【正文】 j,3)))。 end end % 分別計(jì)算 N矩陣的對角及非對角元素 if i nPQ+1 if i~=j K(i,j)=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)))。 K L]，并初始化 Qi = zeros(nb1,1)。 % 對 dP和 dQ賦初值 PV節(jié)點(diǎn)不需計(jì)算 dQ 平衡節(jié)點(diǎn)不參與計(jì)算 for i = 1:n1 for j = 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)))。 Y(J,I)=Y(I,J)。 Y(I,J)=Y(I,J)Yt。 Zt=line(k,3)+j*line(k,4)。 nodenum=[newbus bus(:,1)]。 % increment PV bus counter PV(nPV,:)=bus(i,:)。 [nl,ml]=size(line)。 fprintf(myf,39。% 39。,bus_res(i,3))。 for i = 1:nb fprintf(myf,39。 % 計(jì)算線路潮流的子程序 myf = fopen(39。 % 計(jì)算電壓修正量 bus(1:nPQ,2) = bus(1:nPQ,2) dU。 % 設(shè)定誤差精度 for t = 1:100 % 開始迭代計(jì)算，設(shè)定最大迭代次數(shù)為 100，以便不收斂情況下及時(shí)跳出 [dP,dQ] = dPQ(Y,bus,nPQ,nPV)。,39。 you must select a valid data file39。 計(jì)算結(jié)果牛頓－拉夫遜法潮流計(jì)算結(jié)果 節(jié)點(diǎn)計(jì)算結(jié)果： n 節(jié)點(diǎn) 節(jié) 點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率 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 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 36 + j 37 + j 38 + j 39 + j 線路計(jì)算結(jié)果： n 節(jié)點(diǎn) I 節(jié)點(diǎn) 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 13 10 + j + j + j 14 13 + j + j + j 15 14 + j + j + j 16 15 + j + j + j 17 16 + j + j + j 19 16 + j + j + j 21 16 + j + j + j 24 16 + j + j + j 18 17 + j + j + j 27 17 + j + j + j 22 21 + j + j + j 23 22 + j + j + j 24 23 + j + j + j 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 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 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 結(jié)果分析：此程序的運(yùn)行結(jié)果和試驗(yàn)程序給出的結(jié)果是一致的。 32 10 0 0 。