【文章內(nèi)容簡介】 + j + j + j 23 22 + j + j + j 24 23 + j + j + j 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 26 25 + j + j + j 27 26 + j + j + j 28 26 + j + j + j 29 26 + j + j + j 29 28 + j + j + j 4 0 + j + j + j 5 0 + j + j + j 11 12 + j + j + j 13 12 + j + j + j 30 2 + j + j + j 31 6 + j + j + j 32 10 + j + j + j 34 20 + j + j + j 33 19 + j + j + j 35 22 + j + j + j 36 23 + j + j + j 37 25 + j + j + j 38 29 + j + j + j 20 19 + j + j + j 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 結(jié)果分析：此程序的運行結(jié)果和試驗程序給出的結(jié)果是一致的。說明程序無誤，但在精確度上有微小差異，這主要是和導納矩陣的精確度以及顯示精度有關。 心得：本程序分模塊進行，先是排序，再是求導納陣，然后求雅閣比，再進行迭代運算，程序本身很簡潔明了，運行的時候只需要在 matlab 里輸入 main 就行了，然后打開 BUS 和 line 所在的 .m文件，結(jié)果就會自動存在 result 文件中了，通過編寫牛頓拉夫遜法 matlab 潮流計算程序復習了潮流計算的知識，也實現(xiàn)了計算機算法 附錄： 實驗源程序： Main函數(shù)： clear [dfile,pathname]=uigetfile(39。*.m39。,39。Select Data File39。)。 if pathname == 0 error(39。 you must select a valid data file39。) else lfile =length(dfile)。 % strip off .m eval(dfile(1:lfile2))。 end [nb,mb]=size(bus)。 [nl,ml]=size(line)。 % 計算 bus和 line矩陣的行數(shù)和列數(shù) [bus,line,nPQ,nPV,nodenum] = Num(bus,line)。 % 對節(jié)點重新排序的子程序 Y = y(bus,line) % 計算節(jié)點導納矩陣的子程序 myf = fopen(39。39。,39。w39。)。 fprintf(myf,39。計算結(jié)果 39。)。 fclose(myf)。 % 在當前目錄下生成“ ”文件，寫入節(jié)點導納矩陣 format long EPS = 。 % 設定誤差精度 for t = 1:100 % 開始迭代計算，設定最大迭代次數(shù)為 100，以便不收斂情況下及時跳出 [dP,dQ] = dPQ(Y,bus,nPQ,nPV)。 % 計算功率偏差 dP和 dQ的子程序 J = Jac(bus,Y,nPQ)。 % 計算雅克比矩陣的子程序 UD = zeros(nPQ,nPQ)。 for i = 1:nPQ UD(i,i) = bus(i,2)。 % 生成電 壓對角矩陣 end end dAngU = J\[dP。dQ]。 dAng = dAngU(1:nb1,1)。 % 計算相角修正量 dU = UD*(dAngU(nb:nb+nPQ1,1))。 % 計算電壓修正量 bus(1:nPQ,2) = bus(1:nPQ,2) dU。 % 修正電壓 bus(1:nb1,3) = bus(1:nb1,3) dAng。 % 修正相角 if (max(abs(dU))EPS)amp。(max(abs(dAng))EPS) 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 break end % 判斷是否滿足精度誤差，如滿足則跳出，否則返回繼續(xù)迭代 end bus = PQ(bus,Y,nPQ,nPV)。 % 計算每個節(jié)點的有功和無功注入的子程序 [bus,line] = ReNum(bus,line,nodenum)。 % 對節(jié)點恢復編號的子程序 YtYm = YtYm(line)。 % 計算線路的等效 Yt和 Ym的子程序，以計算線路潮流 bus_res = bus_res(bus)。 % 計算節(jié)點數(shù)據(jù)結(jié)果的子程序 S_res = S_res(bus,line,YtYm)。 % 計算線路潮流的子程序 myf = fopen(39。39。,39。a39。)。 fprintf(myf,39。牛頓－拉夫遜法潮流計算結(jié)果 節(jié)點計算結(jié)果： n節(jié)點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率 \n39。)。 for i = 1:nb fprintf(myf,39。% 39。,bus_res(i,1))。 fprintf(myf,39。% 39。,bus_res(i,2))。 fprintf(myf,39。% 39。,bus_res(i,3))。 fprintf(myf,39。% + j %\n39。,real(bus_res(i,4)),imag(bus_res(i,4)))。 end fprintf(myf,39。n線路計算結(jié)果： n節(jié)點 I 節(jié)點 J 線路功率 S(I,J) 線路功率 S(J,I) 線路損耗 dS(I,J)\n39。)。 for i = 1:nl fprintf(myf,39。% 39。,S_res(i,1))。 fprintf(myf,39。% 39。,S_res(i,2))。 fprintf(myf,39。% + j % 39。,real(S_res(i,3)),imag(S_res(i,3)))。 fprintf(myf,39。% + j % 39。,real(S_res(i,4)),imag(S_res(i,4)))。 fprintf(myf,39。% + j%\n39。,real(S_res(i,5)),imag(S_res(i,5)))。 end fclose(myf)。 % 迭代結(jié)束后繼續(xù)在“ ”寫入節(jié)點計算結(jié)果和線路計算結(jié)果 程序結(jié)束 作用為對節(jié)點重排序，并修改相應的線路數(shù)據(jù) function [bus,line,nPQ,nPV,nodenum] = Num(bus,line) [nb,mb]=size(bus)。 [nl,ml]=size(line)。 nSW = 0。 % number of swing bus counter nPV = 0。 % number of PV bus counter nPQ = 0。 % number of PQ bus counter for i = 1:nb, % nb為總節(jié)點數(shù) type= bus(i,6)。 if type == 3, nSW = nSW + 1。 % increment swing bus counter 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 SW(nSW,:)=bus(i,:)。 elseif type == 2, nPV = nPV +1。 % increment PV bus counter PV(nPV,:)=bus(i,:)。 else nPQ = nPQ + 1。 % increment PQ bus counter PQ(nPQ,:)=bus(i,:)。 end end。 bus=[PQ。PV。SW]。 newbus=[1:nb]39。 nodenum=[newbus bus(:,1)]。 bus(:,1)=newbus。 for i=1:nl for j=1:2 for k=1:nb if line(i,j)==nodenum(k,2) line(i,j)=nodenum(k,1)。 break end end end end 作用為計算節(jié)點導納矩陣 function Y = y(bus,line) [nb,mb]=size(bus)。 [nl,ml]=size(line)。 Y=zeros(nb,nb)。 for k=1:nl I=line(k,1)。 %讀入線路參數(shù) J=line(k,2)。 Zt=line(k,3)+j*line(k,4)。 Yt=1/Zt。 Ym=line(k,5)+j*line(k,6)。 K=line(k,7)。 if (K==0)amp。(J~=0) %