【正文】
+ 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é)果是一致的。 38 29 0 0 。 36 23 0 0 。 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 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。 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 % 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??? ? ? ??????? ? ????? ?? ? ? ? ? ??????? ? ? ? ? ????????? ????????? 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 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???????? ? ? ?????? ? ? ? ?????? 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 三、 求解問題及其結(jié)果 開始 形成節(jié)點(diǎn)導(dǎo)納矩陣 輸入原始數(shù)據(jù) 設(shè)節(jié)點(diǎn)電壓 (0) (0)iief, i=1,2? ,n,i? s 置迭代次數(shù) 0k? 置節(jié)點(diǎn)號(hào) i=1 按式( 33),( 34)計(jì)算雅克比矩陣元素 按式( 32)計(jì)算 PQ 節(jié)點(diǎn)的 ()kiP? , ()kiQ? , PV 節(jié)點(diǎn)的 ()kiP? , ()2kiU? 求解修正方程式,得 ()kie? , ()kif? 雅克比矩陣是否已全部形成? 計(jì)算平衡節(jié)點(diǎn)及 PV節(jié)點(diǎn)功率 求 ()max||ke? , ()max||kf? 迭代次數(shù) k=k+1 i=i+1 ( ) ( )m ax m ax| | ,| |kkef ?? ? ?? 潮流計(jì)算完成 計(jì)算各節(jié)點(diǎn)電壓的新值: ( 1) ( ) ( )k k kie e e? ? ?? ( 1) ( ) ( )k k kif f f? ? ?? 基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 IEEE美國新英格蘭 10機(jī) 39節(jié)點(diǎn)測試系統(tǒng) 一、 系統(tǒng)單線圖 G GGGGGG GG G303912253729172693381654182728362435222120342319331011131415831126327基于牛頓拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 二、系統(tǒng)參數(shù) 1)系統(tǒng)容量基準(zhǔn)值為 100MVA。 ( 2)編程語言為 MATLAB?;谂nD拉夫遜法潮流計(jì)算的 matlab 實(shí)驗(yàn)報(bào)告 一、實(shí)驗(yàn)?zāi)康? 應(yīng)用 MATLAB 語言編寫具有一定通用性的牛 頓 拉夫遜法潮流計(jì)算程序。 要求: ( 1)潮流計(jì)算方法為牛頓 拉夫遜法。 ( 3)程序具有較強(qiáng)通用性。 2) 負(fù)荷數(shù)據(jù)見表 D1 表 D1 負(fù)荷 數(shù)據(jù) 節(jié)點(diǎn)號(hào) 有功( MW) 無功( Mvar) 節(jié)點(diǎn)號(hào) 有功( 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ā)電機(jī)數(shù)據(jù)見表 D2 表 D2 發(fā)電機(jī)數(shù)據(jù) 發(fā)電機(jī) 節(jié)點(diǎn)號(hào) 有功( MW) 無功( Mvar) 電壓( .) PV PQ PQ PV PV PQ PV PV PV BUS30 BUS32 BUS33 BUS34 BUS35 BUS36 BUS37 BUS38 BUS39 平衡節(jié)點(diǎn) BUS31 (電壓幅角 ) (幅值) 4)線路參數(shù)見表 D3 表 D3 線路參數(shù) 序號(hào) 節(jié)點(diǎn) I 節(jié)點(diǎn) 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