電力系統(tǒng)分析課程設(shè)計(jì)---復(fù)雜網(wǎng)絡(luò)牛頓—拉夫遜法潮流分析與計(jì)算的設(shè)計(jì)-預(yù)覽頁(yè)
【正文】 學(xué)出版社， 2021 年版； 《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩 編，水利電力出版社； 復(fù)雜網(wǎng)絡(luò) NR 法潮流分析與計(jì)算的設(shè)計(jì) 電力系統(tǒng)的潮流計(jì)算是電力系統(tǒng)分析課程基本計(jì)算的核心部分之一。 其次，利用計(jì)算機(jī)編程對(duì)電力系統(tǒng)穩(wěn)態(tài)運(yùn)行的各參數(shù)進(jìn)行解析和計(jì)算；編程完成復(fù)雜網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣的形成；電力系統(tǒng)支路改變、節(jié)點(diǎn)增減的程序變化；編程完成各元件的功率損耗、各段網(wǎng)絡(luò)的電壓損耗、各點(diǎn)電壓、功率大小和方向的計(jì)算。 首先，根據(jù)給定的電 力系統(tǒng)簡(jiǎn)圖，通過(guò)手算完成計(jì)算機(jī)算法的兩次迭代過(guò)程，從而加深對(duì)牛頓 拉夫遜法的理解，有助于計(jì)算機(jī)編程的應(yīng)用。 四、進(jìn)程安排 第一周： 第一天上午：選題，查資料，制定設(shè)計(jì)方案； 第一天下午 —— 第五天：復(fù)習(xí)潮流計(jì)算的計(jì)算機(jī)算法，完成給定網(wǎng)絡(luò)的潮流計(jì)算； 第二周： 第一 —— 四天：利用 matlab 編程完成潮流計(jì)算，并對(duì)照手算結(jié)果，分析誤差 第五天下午：答辯，交設(shè)計(jì)報(bào)告。運(yùn)用以直角坐標(biāo)表示的牛頓 拉夫遜計(jì)算如圖 1所示系統(tǒng)中的潮流分布。 表 1 各節(jié)點(diǎn)電壓標(biāo)幺值參數(shù) U1 U2 U3 U4 U5 U6 表 2 線路、變壓器阻抗標(biāo)幺值 線路 L2 L3 L4 L5 T1 T2 Y/2 阻抗 + + + + 表 3 節(jié)點(diǎn)輸出功率 節(jié)點(diǎn) ② ③ ④ ⑤ ⑥ 功率 2+j1 + + + 5 注：各 PQ 節(jié)點(diǎn)的電 壓取 1是為了方便計(jì)算和最后驗(yàn)證程序的正確性。這個(gè)收斂條件比較直觀，它可以直接顯示出最終結(jié)果的功率誤差 如果收斂，就進(jìn)一步計(jì)算各段電力線路潮流和平衡節(jié)點(diǎn)功率，并打印出計(jì)算結(jié)果；如果不收斂，轉(zhuǎn)回 3步進(jìn)行下一次迭代計(jì)算，直到收斂為止。YKK1Y11111yT2T16645255020T22T155453450,4044342340,30,3325233020211221111?????????????????????????????????????KyjyyKyjyyyyyjyyyyyyyyYKKYKyjYKKYKTTTT 非對(duì)角線上的元素為： 00。0。同時(shí) 整個(gè) 求解過(guò)程中用到求節(jié)點(diǎn)導(dǎo)納矩陣，求矩陣的逆陣等等，又讓我們遇到了很多問(wèn)題，我們經(jīng)過(guò)查資料問(wèn)同學(xué) 、問(wèn)老師，最后才把問(wèn)題解決。這些基礎(chǔ)的缺乏給我的設(shè)計(jì)造成了不小的障礙 ，今后我會(huì)努力學(xué)習(xí)來(lái)彌補(bǔ)這些不足 ?？蓪?duì)于大系統(tǒng)來(lái)說(shuō)手工計(jì)算是不可完成的任務(wù)，因此用計(jì)算機(jī)計(jì)算潮流是必不可少的，只要編好程序計(jì)算機(jī)就會(huì)給你一個(gè)自己想要的滿意答案。 y12=1/(*)。 y16=0。 y24=0。 y32=y23。 y36=0。 y44=+1/(+)+1/(+)。 y52=y25。 y56=1/(*)。y64=0。 y21 y22 y23 y24 y25 y26。y61 y62 y63 y64 y65 y66]%形成導(dǎo)納矩陣 %*********牛頓 — 拉夫遜法潮流計(jì)算 ************ E(2)=。E(6)=。F(5)=0。 %提取導(dǎo)納實(shí)部 ,虛部分別為 G,B S(2)=21j。 S(6)=5。precision=。%輸入平衡節(jié)點(diǎn)電壓實(shí)部 ,虛部 e,f for m=2:N1+1 for n=1:N1+1 Pt(n)=(E(m)*(G(m,n)*E(n)B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)))。 %計(jì)算功率誤差 end for m=2:N1+1 for n=1:N1+1 Ai(n)=G(m,n)*E(n)B(m,n)*F(n)。%計(jì)算 dP/df L(m,m)= sum(Bi)B(m,m)*E(m)+G(m,m)*F(m)。 JJ(2*m2,2*m3)=L(m,m)。 %計(jì)算 dP/df L(m,n)=B(m,n)*E(m)+G(m,n)*F(m)。 JJ(2*m2,2*n3)=L(m,n)。 %修正雅可比矩陣第十行 9 元素為 2e JJ(10,10)=2*F(6)。%修正 PQ 矩陣第十個(gè)元素為 du^2 dU=inv(JJ)*PQ39。 F(n+1)=F(n+1)+dU(2*n)。 break end end %************功率及功率損耗的計(jì)算 *********** for m=1:N1+1 I(m)=YB(1,m)*U(m)。 y30=。 Y=[y10 y20 y30 y40 y50 y60]。 %輸出電力線路功率總損耗 dataS及平衡節(jié)點(diǎn)功率 s1及節(jié)點(diǎn)間線路上的流動(dòng)功率矩陣 s
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1