【正文】 作用為計(jì)算雅克比矩陣 [nb,mb] size bus 。 % 對(duì) 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 。 % 總節(jié)點(diǎn)個(gè)數(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。 end if 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: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 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 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ò)編寫(xiě)牛頓拉夫遜法 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 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