【正文】 end function J Fd x J [1,10*x 2 ,14*x 3 。第 %d 次迭代結(jié)果： \n39。 eps 10^10。 0,x 3 *exp x 2 *x 3 ,x 2 *exp x 2 *x 3 +20]。 b F x 。 x 5*cos a 。 S_res k,5 S_res k,3 +S_res k,4 。 J S_res k,2 。 bus_res zeros nb,4 。 J 0 % 對(duì)地支路 : K 0,J 0,R X 0 YtYm k,4 Ym。 if Zt~ 0 Yt 1/Zt。 % 重新賦值回 bus，完成 bus 矩陣的重新編號(hào) for i 1:nl for j 1:2 for k 1:nb if line i,j nodenum k,1 line i,j nodenum k,2 。 bus n,5 bus n,5 +bus n,2 *bus j,2 * real Y n,j *sin bus n,3 bus j,3 imag Y n,j *cos bus n,3 bus j,3 。 else K i,j Pi i,1 +real Y i,j * bus i,2 ^2 。 Pi zeros nb1,1 。 if i nPQ+1 dQ i,1 dQ i,1 bus i,2 *bus j,2 * real Y i,j *sin bus i,3 bus j,3 imag Y i,j *cos bus i,3 bus j,3 。 end if K 0 % 變壓器線路 : Zt和 Ym為折算到 K側(cè)的值 ,K在 i 側(cè) Y I,I Y I,I +Yt+Ym。 Y J,I Y I,J 。 Yt 1/Zt。 bus :,1 newbus。 else nPQ nPQ + 1。 nSW 0。% + j % 39。,S_res i,1 。 fprintf myf,39。% 39。39。 % 修正電壓 bus 1:nb1,3 bus 1:nb1,3 dAng。 % 計(jì)算功率偏差dP 和 dQ 的子程序 J Jac bus,Y,nPQ 。w39。 else lfile length dfile 。說(shuō)明程序無(wú)誤，但在精確度上有微小差異，這主要是和導(dǎo)納矩陣的精確度以及顯示精度有關(guān)。 34 20 0 0 。 29 28 0 0。 27 17 0 0。 14 13 0 0。 8 5 0 0。 39 1 0 0。畢業(yè)論文 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 一、實(shí)驗(yàn)?zāi)康? 應(yīng)用 MATLAB 語(yǔ)言編寫(xiě)具有一定通用性的牛 頓 拉夫遜法潮流計(jì)算程序。 3 2 0 0。 7 6 0 0。 15 14 0 0。 22 21 0 0。 4 0 0 0 0 0。 33 19 0 0 。 心得：本程序分模塊進(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。 % strip off .m eval dfile 1:lfile2 。 。 % 計(jì)算雅克比矩陣的子程序 UD zeros nPQ,nPQ 。 % 修正相角 if abs dU EPS amp。,39。,bus_res i,1 。% + j %\n39。 fprintf myf,39。,real S_res i,4 ,imag S_res i,4 。 % number of swing bus counter nPV 0。 % increment PQ bus counter PQ nPQ,: bus i,: 。 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 。 Ym line k,5 +j*line k,6 。 end if K 0 amp。 Y J,J Y J,J +K*K*Yt。 end end end % 利用循環(huán)計(jì)算求取 dP 和 dQ 作用為計(jì)算雅克比矩陣 [nb,mb] size bus 。 for i 1:nb1 for j 1:nb Pi i,1 Pi 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 。 end % 分別計(jì)算 K 矩陣的對(duì)角及非對(duì)角元素 if j nPQ+1 if i~ j L i,j bus i,2 *bus j,2 * real Y i,j *sin bus i,3 bus j,3 imag Y i,j *cos bus i,3 bus j,3 。 end % 計(jì)算平衡節(jié)點(diǎn)的無(wú)功及有功注入 作用為對(duì)節(jié)點(diǎn)和線路數(shù)據(jù)恢復(fù)編號(hào) [nb,mb] size bus 。 break end end end end % 恢復(fù) line 的編號(hào) 作用為計(jì)算線路的等效 Yt 和 Ym，以計(jì)算線路潮流 [nl,ml] size line 。 end Ym line k,5 +j*line k,6 。 end if K 0 % 變壓器線路 : Zt 和 Ym 為折算到 i 側(cè)的值 ,K 在 j 側(cè) YtYm k,3 Yt/K。 % bus_res 矩陣儲(chǔ)存著節(jié)點(diǎn)計(jì)算結(jié)果 bus_res :,1:2 bus :,1:2 。 if J~ 0 amp。 % 利用公式計(jì)算接地支路的潮流 end end end end 作業(yè) y pifun 10000 y clear all tic。 y 3*sin a 。 dx A\b。 end function y F x y 1 3*x 1 cos x 2 *x 3 。 x0 [ ]39。,i [x dx F x ] %在屏幕上輸出每次的 x i ,dx i ,F x i if abs dx eps amp。 3*x 2 +x 3 11, 3*x 1 ,x 1 。y 3 ]。 sprintf 39。 end format long 1 100。 2*x 1 , 162* x 2 + ,cos x 3 。 for i 1:1, A Fd x 。 end a 0:pi/12:2*pi。 else S_res k,4 bus J,2 ^2*conj YtYm k,5 。 % 前兩列為節(jié)點(diǎn)編號(hào) for k 1:nl I S_res k,1 。 end end 計(jì)算并返回節(jié)點(diǎn)數(shù)據(jù)結(jié)果 [nb,mb] size bus 。 end if K 0 amp。