畢業(yè)論文--基于牛頓拉夫遜法潮流計(jì)算的matlab實(shí)驗(yàn)報(bào)告-資料下載頁
2024-11-14 21:41本頁面
【導(dǎo)讀】頓-拉夫遜法潮流計(jì)算程序。程序具有較強(qiáng)通用性。1）系統(tǒng)容量基準(zhǔn)值為100MVA。
【正文】 YtYm k,4 Ym。 YtYm k,5 Ym。 end if K 0 amp。 J 0 % 對(duì)地支路 : K 0,J 0,R X 0 YtYm k,4 Ym。 end if K 0 % 變壓器線路 : Zt 和 Ym 為折算到 i 側(cè)的值 ,K 在 j 側(cè) YtYm k,3 Yt/K。 YtYm k,4 Ym+Yt* K1 /K。 YtYm k,5 Yt* 1K /K/K。 end if K 0 % 變壓器線路 : Zt 和 Ym 為折算到 K 側(cè)的值 ,K 在 i 側(cè) YtYm k,3 Yt*K。 YtYm k,4 Ym+Yt* 1+K 。 YtYm k,5 Yt* K^2+K 。 end end 計(jì)算并返回節(jié)點(diǎn)數(shù)據(jù)結(jié)果 [nb,mb] size bus 。 bus_res zeros nb,4 。 % bus_res 矩陣儲(chǔ)存著節(jié)點(diǎn)計(jì)算結(jié)果 bus_res :,1:2 bus :,1:2 。 bus_res :,3 bus :,3 *180 / pi。 % 相角采用角度制 bus_res :,4 bus :,4 + sqrt 1 *bus :,5 。 % 注入功率 end 計(jì)算并返回線路潮流 [nl,ml] size line 。 S_res zeros nl,5 。 % S_res 矩陣儲(chǔ)存著線路潮流計(jì)算結(jié)果 S_res :,1:2 line :,1:2 。 % 前兩列為節(jié)點(diǎn)編號(hào) for k 1:nl I S_res k,1 。 J S_res k,2 。 if J~ 0 amp。 I~ 0 S_res k,3 bus I,2 ^2* conj YtYm k,3 +conj YtYm k,4 bus I,2 *bus J,2 * cos bus I,3 +j*sin bus I,3 * conj cos bus J,3 +j*sin bus J,3 *conj YtYm k,3 。 S_res k,4 bus J,2 ^2* conj YtYm k,3 +conj YtYm k,5 bus I,2 *bus J,2 * conj cos bus I,3 +j*sin bus I,3 * cos bus J,3 +j*sin bus J,3 *conj YtYm k,3 。 S_res k,5 S_res k,3 + S_res k,4 。 % 利用公式計(jì)算非接地支路的潮流 else if J 0 S_res k,3 bus I,2 ^2*conj YtYm k,4 。 S_res k,5 S_res k,3 +S_res k,4 。 else S_res k,4 bus J,2 ^2*conj YtYm k,5 。 S_res k,5 S_res k,3 +S_res k,4 。 % 利用公式計(jì)算接地支路的潮流 end end end end 作業(yè) y pifun 10000 y clear all tic。 y pifun 100000 。 toc。 Elapsed time is seconds. function y pifun n y 0。 for i 1:n y y+ 1 ^ i1 *1/ 2*i1 。 end y y*4。 end a 0:pi/12:2*pi。 x 5*cos a 。 y 3*sin a 。 plot x,y clear format long 1 100。 eps1 。 eps2 。 x0 [ ]39。 x x0。 for i 1:1, A Fd x 。 b F x 。 dx A\b。 %矩 陣左除，即 b 除以 A x xdx。 sprintf 39。第 %d 次迭代結(jié)果： \n39。,i [x dx F x ] %在屏幕上輸出每次的 x i ,dx i ,F x i if abs dx eps1 amp。 abs F x eps2 break end end ans 第 1 次迭代結(jié)果： ans ans 第 2 次迭代結(jié)果： ans ans 第 3 次迭代結(jié)果： ans ans 第 4 次迭代結(jié)果： ans ans 第 5 次迭代結(jié)果： ans ans 第 6 次迭代結(jié)果： ans ans 第 7 次迭代結(jié)果： ans ans 第 8 次迭代結(jié)果： ans 0 0 x function J Fd x J [3,x 3 *sin x 2 *x 3 ,x 2 *sin x 2 *x 3 。 2*x 1 , 162* x 2 + ,cos x 3 。 0,x 3 *exp x 2 *x 3 ,x 2 *exp x 2 *x 3 +20]。 end function y F x y 1 3*x 1 cos x 2 *x 3 。 y 2 x 1 ^281* x 2 + ^2+sin x 3 +。 y 3 exp x 2 *x 3 +20*x 3 + 10*pi3 /3。 y [y 1 。y 2 。y 3 ]。 end format long 1 100。 eps 10^10。 x0 [ ]39。 x x0。 for i 1:1, A Fd x 。 b F x 。 dx A\b。 %矩陣左除，即 b 除以 A x xdx。 sprintf 39。第 %d 次迭代結(jié)果： \n39。,i [x dx F x ] %在屏幕上輸出每次的 x i ,dx i ,F x i if abs dx eps amp。 abs F x eps break end end ans 第 1 次迭代結(jié)果： ans +002 * ans 第 2 次迭代結(jié)果： ans ans 第 3 次迭代結(jié)果： ans ans 第 4 次迭代結(jié)果： ans ans 第 5 次迭代結(jié)果： ans ans 第 6 次迭代結(jié)果： ans 0 0 0 x 1 5 4 function y F x y 1 x 1 5*x 2 ^2+7*x 3 ^2+12。 y 2 3*x 1 *x 2 +x 1 *x 3 11*x 1 。 y 3 2*x 2 *x 3 +40*x 1 。 y [y 1 。y 2 。y 3 ]。 end function J Fd x J [1,10*x 2 ,14*x 3 。 3*x 2 +x 3 11, 3*x 1 ,x 1 。 40,2*x 3 ,2*x 2 ]。 end 開始 形成節(jié)點(diǎn)導(dǎo)納矩陣 輸入原始數(shù)據(jù) 設(shè)節(jié)點(diǎn)電壓， i 1,2? ,n,is 置迭代次數(shù) 置節(jié)點(diǎn)號(hào) i 1 按式（ 33），（ 34）計(jì)算雅克比矩陣元素 按式（ 32）計(jì)算節(jié)點(diǎn)的，節(jié)點(diǎn)的， 求解修正方程式，得， 雅克比矩陣是否已全部形成？ 計(jì)算平衡節(jié)點(diǎn)及 PV 節(jié)點(diǎn)功率 求， 迭 代次數(shù) k k+1 i i+1 ？ 潮流計(jì)算完成 計(jì)算各節(jié)點(diǎn)電壓的新值：
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1