【正文】 1:n dP(i,1) = dP(i,1)bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)bus(j,3))+imag(Y(i,j))*sin(bus(i,3)bus(j,3)))。 % 總節(jié)點個數(shù)dP = bus(1:n1,4)。Y(J,I)=Y(I,J)。Y(J,J)=Y(J,J)+K*K*Yt。Y(J,I)=Y(I,J)。Y(J,J)=Y(J,J)+Yt/K/K。Y(I,I)=Y(I,I)+Ym。endif (K==0)amp。Y(I,J)=Y(I,J)Yt。Y(I,I)=Y(I,I)+Yt+Ym。if (K==0)amp。Ym=line(k,5)+j*line(k,6)。Zt=line(k,3)+j*line(k,4)。for k=1:nlI=line(k,1)。[nl,ml]=size(line)。for i=1:nlfor j=1:2for k=1:nbif line(i,j)==nodenum(k,2)line(i,j)=nodenum(k,1)。nodenum=[newbus bus(:,1)]。SW]。bus=[PQ。 % increment PQ bus counterPQ(nPQ,:)=bus(i,:)。 % increment PV bus counterPV(nPV,:)=bus(i,:)。 % increment swing bus counterSW(nSW,:)=bus(i,:)。 % number of PQ bus counterfor i = 1:nb, % nb為總節(jié)點數(shù)type= bus(i,6)。 % number of swing bus counternPV = 0。[nl,ml]=size(line)。endfclose(myf)。% + j%\n39。,real(S_res(i,4)),imag(S_res(i,4)))。fprintf(myf,39。% + j % 39。,S_res(i,2))。fprintf(myf,39。% 39。)。endfprintf(myf,39。% + j %\n39。,bus_res(i,3))。 fprintf(myf,39。% 39。,bus_res(i,1))。for i = 1:nbfprintf(myf,39。牛頓－拉夫遜法潮流計算結(jié)果 節(jié)點計算結(jié)果：n節(jié)點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率\n39。)。,39。 % 計算線路潮流的子程序 myf = fopen(39。 % 計算線路的等效Yt和Ym的子程序，以計算線路潮流 bus_res = bus_res(bus)。 % 計算每個節(jié)點的有功和無功注入的子程序[bus,line] = ReNum(bus,line,nodenum)。 % 修正相角 if (max(abs(dU))EPS)amp。 % 計算電壓修正量bus(1:nPQ,2) = bus(1:nPQ,2) dU。dAng = dAngU(1:nb1,1)。 % 生成電壓對角矩陣 endenddAngU = J\[dP。 % 計算雅克比矩陣的子程序UD = zeros(nPQ,nPQ)。 % 設(shè)定誤差精度for t = 1:100 % 開始迭代計算，設(shè)定最大迭代次數(shù)為100，以便不收斂情況下及時跳出[dP,dQ] = dPQ(Y,bus,nPQ,nPV)。fclose(myf)。計算結(jié)果39。)。,39。 % 對節(jié)點重新排序的子程序Y = y(bus,line) % 計算節(jié)點導納矩陣的子程序myf = fopen(39。[nl,ml]=size(line)。 % strip off .m eval(dfile(1:lfile2))。 you must select a valid data file39。)。,39。心得：本程序分模塊進行，先是排序，再是求導納陣，然后求雅閣比，再進行迭代運算，程序本身很簡潔明了，運行的時候只需要在matlab里輸入main就行了，結(jié)果就會自動存在result文件中了，通過編寫牛頓拉夫遜法matlab潮流計算程序復習了潮流計算的知識，也實現(xiàn)了計算機算法附錄：實驗源程序：Main函數(shù)：clear[dfile,pathname]=uigetfile(39。計算結(jié)果牛頓－拉夫遜法潮流計算結(jié)果 節(jié)點計算結(jié)果： n節(jié)點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率 1 + j 2 + j 3 + j 4 + j 5 + j 6 + j 7 + j 8 + j 9 + j 10 + j 11 + j 12 + j 13 + j 14 + j 15 + j 16 + j 17 + j 18 + j 19 + j 20 + j 21 + j 22 + j 23 + j 24 + j 25 + j 26 + j 27 + j 28 + j 29 + j 30 + j 31 + j 32 + j 33 + j 34 + j 35 + j 36 + j 37 + j 38 + j 39 + j 線路計算結(jié)果：n節(jié)點I 節(jié)點J 線路功率S(I,J) 線路功率S(J,I) 線路損耗dS(I,J) 2 1 + j + j + j 39 1 + j + j + j 3 2 + j + j + j 25 2 + j + j + j 4 3 + j + j + j 18 3 + j + j + j 5 4 + j + j + j 14 4 + j + j + j 6 5 + j + j + j 8 5 + j + j + j 7 6 + j + j + j 11 6 + j + j + j 8 7 + j + j + j 9 8 + j + j + j 39 9 + j + j + j 11 10 + j + j + j