【正文】 4)。 for k=1:nl I=line(k,1)。 [nl,ml]=size(line)。 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)。 nodenum=[newbus bus(:,1)]。SW]。 bus=[PQ。 % increment PQ bus counter PQ(nPQ,:)=bus(i,:)。 % increment PV bus counter PV(nPV,:)=bus(i,:)。 % increment swing bus counter 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 SW(nSW,:)=bus(i,:)。 % number of PQ bus counter for i = 1:nb, % nb為總節(jié)點(diǎn)數(shù) type= bus(i,6)。 % number of swing bus counter nPV = 0。 [nl,ml]=size(line)。 end fclose(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。)。 end fprintf(myf,39。% + j %\n39。,bus_res(i,3))。 fprintf(myf,39。% 39。,bus_res(i,1))。 for i = 1:nb fprintf(myf,39。牛頓－拉夫遜法潮流計(jì)算結(jié)果 節(jié)點(diǎn)計(jì)算結(jié)果： n節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率 \n39。)。,39。 % 計(jì)算線路潮流的子程序 myf = fopen(39。 % 計(jì)算線路的等效 Yt和 Ym的子程序，以計(jì)算線路潮流 bus_res = bus_res(bus)。 % 計(jì)算每個(gè)節(jié)點(diǎn)的有功和無(wú)功注入的子程序 [bus,line] = ReNum(bus,line,nodenum)。 % 修正相角 if (max(abs(dU))EPS)amp。 % 計(jì)算電壓修正量 bus(1:nPQ,2) = bus(1:nPQ,2) dU。 dAng = dAngU(1:nb1,1)。 % 生成電 壓對(duì)角矩陣 end end dAngU = J\[dP。 % 計(jì)算雅克比矩陣的子程序 UD = zeros(nPQ,nPQ)。 % 設(shè)定誤差精度 for t = 1:100 % 開(kāi)始迭代計(jì)算，設(shè)定最大迭代次數(shù)為 100，以便不收斂情況下及時(shí)跳出 [dP,dQ] = dPQ(Y,bus,nPQ,nPV)。 fclose(myf)。計(jì)算結(jié)果 39。)。,39。 % 對(duì)節(jié)點(diǎn)重新排序的子程序 Y = y(bus,line) % 計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣的子程序 myf = fopen(39。 [nl,ml]=size(line)。 % strip off .m eval(dfile(1:lfile2))。 you must select a valid data file39。)。,39。 心得：本程序分模塊進(jìn)行，先是排序，再是求導(dǎo)納陣，然后求雅閣比，再進(jìn)行迭代運(yùn)算，程序本身很簡(jiǎn)潔明了，運(yùn)行的時(shí)候只需要在 matlab 里輸入 main 就行了，然后打開(kāi) BUS 和 line 所在的 .m文件，結(jié)果就會(huì)自動(dòng)存在 result 文件中了，通過(guò)編寫牛頓拉夫遜法 matlab 潮流計(jì)算程序復(fù)習(xí)了潮流計(jì)算的知識(shí)，也實(shí)現(xiàn)了計(jì)算機(jī)算法 附錄： 實(shí)驗(yàn)源程序： Main函數(shù)： clear [dfile,pathname]=uigetfile(39。 計(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