【正文】 + 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 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 結(jié)果分析：此程序的運行結(jié)果和試驗程序給出的結(jié)果是一致的。 38 29 0 0 。 36 23 0 0 。 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 33 19 0 0 。 32 10 0 0 。 30 2 0 0 。 11 12 0 0 。 4 0 0 0 0 0。 29 26 0 0。 27 26 0 0。 24 23 0 0。 22 21 0 0。 18 17 0 0。 21 16 0 0。 17 16 0 0。 15 14 0 0。 13 10 0 0。 39 9 0 0。 8 7 0 0。 7 6 0 0。 6 5 0 0。 5 4 0 0。 4 3 0 0。 3 2 0 0。 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 % 1 2 3 4 5 6 7 % line: from bus to bus R, X, G, B/2 K line=[ 2 1 0 0。 二、程序流程圖 所用公式 112 2 2 2[ ( ) ( ) ][ ( ) ( ) ]()jni i i ij j ij j i ij j ij jjjni i i ij j ij j i ij j ij jji i i iP P e G e B f f G f B eQ Q f G e B f e G f B eU U e f?????? ? ? ? ? ?????? ? ? ? ? ???? ? ? ? ?????? ij? 2200iij ij i ij iiiij ij i ij iiiij ij i ij i ijiiij ij i ij i ijiiijiiijiPH B e G ffPN G e B feQJ G e B f NfQL B e G f HeURfUSe??? ? ? ??????? ? ????? ?? ? ? ? ? ??????? ? ? ? ? ????????? ????????? 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 ij?2222iii ij i ij i iiiiii ij i ij i iiiiii ij i ij i iiiiii ij i ij i iiiiii iiiii iiPH B e G f bfPN G e B f aeQJ G e B f afQL B e G f beURffUSee??? ? ? ? ??????? ? ? ????? ?? ? ? ? ??????? ? ? ? ? ????????? ????????? 其中 11()()jnii ii i ii i ij j ij jjjijnii ii i ii i ij j ij jjjia G e B f G e B fb G f B e G f B e???????? ? ? ?????? ? ? ? ?????? 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 三、 求解問題及其結(jié)果 開始 形成節(jié)點導(dǎo)納矩陣 輸入原始數(shù)據(jù) 設(shè)節(jié)點電壓 (0) (0)iief， i=1,2? ,n,i? s 置迭代次數(shù) 0k? 置節(jié)點號 i=1 按式（ 33），（ 34）計算雅克比矩陣元素 按式（ 32）計算 PQ 節(jié)點的 ()kiP? ， ()kiQ? ， PV 節(jié)點的 ()kiP? ， ()2kiU? 求解修正方程式，得 ()kie? ， ()kif? 雅克比矩陣是否已全部形成？ 計算平衡節(jié)點及 PV節(jié)點功率 求 ()max||ke? ， ()max||kf? 迭代次數(shù) k=k+1 i=i+1 ( ) ( )m ax m ax| | ,| |kkef ?? ? ?？ 潮流計算完成 計算各節(jié)點電壓的新值： ( 1) ( ) ( )k k kie e e? ? ?? ( 1) ( ) ( )k k kif f f? ? ?? 基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 IEEE美國新英格蘭 10機 39節(jié)點測試系統(tǒng) 一、 系統(tǒng)單線圖 G GGGGGG GG G303912253729172693381654182728362435222120342319331011131415831126327基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 二、系統(tǒng)參數(shù) 1）系統(tǒng)容量基準(zhǔn)值為 100MVA。 （ 2）編程語言為 MATLAB。基于牛頓拉夫遜法潮流計算的 matlab 實驗報告 一、實驗?zāi)康? 應(yīng)用 MATLAB 語言編寫具有一定通用性的牛 頓 拉夫遜法潮流計算程序。 要求： （ 1）潮流計算方法為牛頓 拉夫遜法。 （ 3）程序具有較強通用性。 2) 負(fù)荷數(shù)據(jù)見表 D1 表 D1 負(fù)荷 數(shù)據(jù) 節(jié)點號 有功（ MW） 無功（ Mvar） 節(jié)點號 有功（ MW） 無功（ Mvar） BUS 3 BUS 4 BUS 7 BUS 8 BUS 12 BUS 15 BUS 16 BUS 18 BUS 20 BUS 21 BUS 23 BUS 24 BUS 25 BUS 26 BUS 27 BUS 28 BUS 29 BUS 39 3）發(fā)電機數(shù)據(jù)見表 D2 表 D2 發(fā)電機數(shù)據(jù) 發(fā)電機 節(jié)點號 有功（ MW） 無功（ Mvar） 電壓（ .） PV PQ PQ PV PV PQ PV PV PV BUS30 BUS32 BUS33 BUS34 BUS35 BUS36 BUS37 BUS38 BUS39 平衡節(jié)點 BUS31 (電壓幅角 ) （幅值） 4）線路參數(shù)見表 D3 表 D3 線路參數(shù) 序號 節(jié)點 I 節(jié)點 J R(.) X(.) B/2(.) LN1 LN2 LN3 LN4 LN5 LN6 LN7 LN8 LN9 BUS2 BUS39 BUS3 BUS25 BUS4 BUS18 BUS5 BUS14 BUS6 BUS1 BUS1 BUS2 BUS2 BUS3 BUS3 BUS4 BUS4 BUS5