基于極坐標(biāo)牛頓拉夫遜潮流計(jì)算修正版

1、2 模型簡(jiǎn)介及等值電路2.1 課程設(shè)計(jì)模型：模型3電力網(wǎng)絡(luò)接線如下圖所示，各支路阻抗標(biāo)幺值參數(shù)如下：Z12=0.02+j0.06，Z13=0.08+j0.24， Z23=0.06+j0.18， Z24=0.06+j0.12， Z25=0.04+j0.12， Z34=0.01+j0.03， Z45=0.08+j0.24， k=1.1。該系統(tǒng)中，節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，保持為定值；節(jié)點(diǎn)2、3、4都是PQ節(jié)點(diǎn)，節(jié)點(diǎn)5為PV節(jié)點(diǎn)，給定的注入功率分別為:， ，。各節(jié)點(diǎn)電壓（初值）標(biāo)幺值參數(shù)如下： 節(jié)點(diǎn) 12345Ui（0）=ei（0）+jfi（0）1.06+j0.01.0+j0.01.0+j0.01.0+j0

2、.01.1+j0.0計(jì)算該系統(tǒng)的潮流分布。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于10-5。圖2-12.2模型分析節(jié)點(diǎn)類型介紹按變量的不同，一般將節(jié)點(diǎn)分為三種類型。1 節(jié)點(diǎn)這類節(jié)點(diǎn)的有功功率和無(wú)功功率是給定的，節(jié)點(diǎn)是待求量。通常變電所都是這一類型節(jié)點(diǎn)。由于沒(méi)有發(fā)電設(shè)備，故其發(fā)電功率為零。有些情況下，系統(tǒng)中某些發(fā)電廠輸出的功率在一段時(shí)間內(nèi)是固定時(shí)，該發(fā)電廠母線也作為節(jié)點(diǎn)。因此，電力系統(tǒng)中絕大多數(shù)節(jié)點(diǎn)屬于這一類型。2 節(jié)點(diǎn)這類節(jié)點(diǎn)有功功率P和電壓幅值V是給定的，節(jié)點(diǎn)的無(wú)功功率Q和電壓的相位是待求量。這類節(jié)點(diǎn)必須有足夠的可調(diào)無(wú)功容量，用以維持給定的電壓幅值，因此又稱為電壓控制節(jié)點(diǎn)。一般選擇有一定無(wú)功儲(chǔ)備的

3、發(fā)電廠和具有可調(diào)無(wú)功電源設(shè)備的變電所作為節(jié)點(diǎn)。3 平衡節(jié)點(diǎn)在潮流分布算出以前，網(wǎng)絡(luò)中的功率損耗是未知的，因此，網(wǎng)絡(luò)中至少有一個(gè)節(jié)點(diǎn)的有功功率P不能給定，這個(gè)節(jié)點(diǎn)承擔(dān)了系統(tǒng)系統(tǒng)的有功功率平衡，故稱之為平衡節(jié)點(diǎn)。1.2、各節(jié)電參數(shù)：由模型中所給列出下表：各節(jié)點(diǎn)電壓和注入功率（初值）標(biāo)幺值參數(shù)如下： 節(jié)點(diǎn) 12345Ui（0）=ei（0）+jfi（0）1.06+j0.01.0+j0.01.0+j0.01.0+j0.01.1+j0.000.20+j0.20-0.45-j0.151.0+j0.01.1+j0.0 各節(jié)點(diǎn)之間的導(dǎo)納：y12=5.000-j15.000，y13=1.2500-j3.7500，

4、y22=0.2750-j0.8250，y23=1.667-j5.000，y24=3.333-j6.667，y25=2.7500-j8.2500，y34=10.0000-j30.0000，y55=-0.25+j0.751.3 等值電路模型 由于計(jì)算時(shí)一般將平衡節(jié)點(diǎn)放到最大編號(hào)，故在本模型中將節(jié)點(diǎn)2、3、4、5、1分別替換為節(jié)點(diǎn)1、2、3、4、5，也即是4換為PV節(jié)點(diǎn)，5為平衡節(jié)點(diǎn)。將變壓器用等值電路，由此繪制等值電路如下：y12=1.667-j5，y13=3.333-j6.667，y14=2.75-j8.25，y15=5-j15，y52=1.25-j3.75，y23=10-j30，y34=1.2

5、5-j3.75，y11=0.275-j0.825，y44=-0.25+j0.75。2.3 等值電路模型 在圖2-2中，將圖2-1中的編號(hào)重新編排，節(jié)點(diǎn)2、3、4、5、1替換為1、2、3、4、5。則各節(jié)點(diǎn)之間的導(dǎo)納變?yōu)閥12=1.667-j5，y13=3.333-j6.667，y14=2.75-j8.25，y15=5-j15，y52=1.25-j3.75，y23=10-j30，y34=1.25-j3.75，y11=0.275-j0.825，y44=-0.25+j0.75。3 設(shè)計(jì)原理本題采用了題目要求的牛頓拉夫遜潮流計(jì)算的方法。牛頓法是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。解決電力系

6、統(tǒng)潮流計(jì)算問(wèn)題是以導(dǎo)納距陣為基礎(chǔ)的，因此，只要在迭代過(guò)程中盡可能保持方程式系數(shù)距陣的稀疏性，就可以大大提高牛頓法潮流程序的放率。3.1潮流計(jì)算的定解條件題中所給圖表示一個(gè)五節(jié)點(diǎn)的簡(jiǎn)單電力系統(tǒng)，n個(gè)節(jié)點(diǎn)電力系統(tǒng)的潮流方程的一般形式是 或 3.2潮流計(jì)算的約束條件(1)所有節(jié)點(diǎn)電壓必須滿足(2)所有電源節(jié)點(diǎn)的有功功率和無(wú)功功率必須滿足(3)某些節(jié)點(diǎn)之間的電壓應(yīng)滿足3.3牛頓-拉夫遜的基本原理 設(shè)欲求解的非線性代數(shù)方程為設(shè)方程的真實(shí)解為，則必有。用牛頓-拉夫遜法求方程真實(shí)解的步驟如下： 首先選取余割合適的初始估值作為方程的解，若恰巧有，則方程的真實(shí)解即為若，則做下一步。取 則其中為初始估值的增量，即

7、。設(shè)函數(shù)具有任意階導(dǎo)數(shù) 若所取的足夠小，則含的項(xiàng)及其余的一切高階項(xiàng)均可略去，并使其等于零，即： 故得 可見(jiàn)，只要0，即可根據(jù)上式求出第一次的修正估值，若恰巧有 =0，則方程的真實(shí)解即為。若，則用上述方法由再確定第二次的修正估值。如此反復(fù)疊代下去，直到求得真實(shí)解為止。設(shè)第K次的估值為第（K+1）次的修正估值，則有 迭代過(guò)程的收斂數(shù)據(jù)為 或其中，為預(yù)先給定的小正數(shù)。4 修正方程的建立極坐標(biāo)表示的牛拉法修正方程為： 式中留出了（n-m）行空格和（n-m）列空列。式中的有功，無(wú)功功率不平衡量分別由式（3-1a）,式（3-1b）可得為 （3-1a，b）而式中雅可比矩陣的各個(gè)元素則分別為 ；; （3-2）

8、式（4-44）中將改為只是為使公式（4-46）中個(gè)偏導(dǎo)數(shù)的表示形式上更相似，為求取這些偏導(dǎo)數(shù)，可將、分別展開(kāi)如下 （3-3a,b）計(jì)及 （3-4）時(shí)，由于對(duì)特定的j，只有該特定節(jié)點(diǎn)的，從而特定的是變量，由式（4-46）到式（4-48）可得 （3-5a）相似的，由于對(duì)特定的j，只有該特定節(jié)點(diǎn)的是變量，可得 （3-5b）j=i時(shí)，由于是變量，所有都是變量，可得 （3-5c）相似的，由于是變量，可得 （3-5d）5 設(shè)計(jì)流程圖及程序的編寫(xiě)輸入數(shù)據(jù)形成節(jié)點(diǎn)導(dǎo)納矩陣設(shè)定節(jié)點(diǎn)起始計(jì)算電壓 U置迭代次數(shù)應(yīng)用公式計(jì)算不平衡量PQ雅克比矩陣J是否已經(jīng)全部形成，in?按公式計(jì)算雅克比矩陣元素J增大節(jié)點(diǎn)號(hào)，i=i+

9、1解修正方程，由P,Q和J計(jì)算電壓修正量u和求出u和迭代是否收斂，u，計(jì)算平衡節(jié)點(diǎn)功率S和線路的功率S結(jié)束K=k+1u=U+u =+置節(jié)點(diǎn)號(hào)i=1程序中用到的符號(hào)所代表的意義：Y 代表導(dǎo)納矩陣JJ 代表雅克比矩陣pp代表有功功率的不平衡量qq 代表無(wú)功功率的不平衡量uu 代表各節(jié)點(diǎn)電壓和相角的不平衡量U 代表各節(jié)點(diǎn)的電壓S 代表線路的功率Q4 代表PV節(jié)點(diǎn)的注入無(wú)功功率K 代表迭代次數(shù)N1 代表PQ節(jié)點(diǎn)和PV節(jié)點(diǎn)的總數(shù)m,n代表系統(tǒng)中的節(jié)點(diǎn)總數(shù)，把平衡節(jié)點(diǎn)標(biāo)為最大號(hào)S5 平衡節(jié)點(diǎn)功率d 變壓器對(duì)地導(dǎo)納%The following program for load calculation is

10、based on MATLAB6.5%以下部分為輸入原始數(shù)據(jù)（到標(biāo)示/標(biāo)志為止）。%求取節(jié)點(diǎn)導(dǎo)納矩陣g(1,1)=0.275;b(1,1)=0.825;g(1,2)=1.667;b(1,2)=-5;g(1,3)=3.333;b(1,3)=-6.667;g(1,4)=2.75;b(1,4)=-8.25;g(1,5)=5;b(1,5)=-15;g(2,1)=1.667;b(2,1)=-5;g(2,2)=-0;b(2,2)=0;g(2,3)=10;b(2,3)=-30;g(2,4)=0;b(2,4)=0;g(2,5)=1.25;b(2,5)=-3.75;g(3,1)=3.33;b(3,1)=-6.6

11、67;g(3,2)=10;b(3,2)=-30;g(3,3)=0;b(3,3)=0;g(3,4)=1.25;b(3,4)=-3.75;g(3,5)=0;b(3,5)=0;g(4,1)=2.75;b(4,1)=-8.25;g(4,2)=0;b(4,2)=0;g(4,3)=1.24;b(4,3)=-3.75;g(4,4)=-0.25;b(4,4)=0.75;g(4,5)=0;b(4,5)=0;g(5,1)=5;b(5,1)=-15;g(5,2)=1.25;b(5,2)=-3.75;g(5,3)=0;b(5,3)=0;g(5,4)=0;b(5,4)=0;g(5,5)=0;b(5,5)=0;d(1,4

12、)=0.275-0.825*j;d(4,1)=-0.25+0.75*jd(1,5)=0;d(5,1)=0;%求取節(jié)點(diǎn)導(dǎo)納矩陣。for m=1:5 for n=1:5 if m=n G(m,m)=g(m,1)+g(m,2)+g(m,3)+g(m,4)+g(m,5); B(m,m)=b(m,1)+b(m,2)+b(m,3)+b(m,4)+b(m,5); else G(m,n)=-g(m,n); B(m,n)=-b(m,n); end endendY=G+j*B%/%/設(shè)定節(jié)點(diǎn)起始節(jié)點(diǎn)電壓，并給出已知功率值。%/下面將題中節(jié)點(diǎn) 2、3、4、5、1分別替換為節(jié)點(diǎn)1、2、3、4、5即節(jié)點(diǎn)4為PV節(jié)點(diǎn)，節(jié)

13、點(diǎn)5為平衡節(jié)點(diǎn)。delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;u(1)=1.0;u(2)=1.0;u(3)=1.0;p(1)=0.2;q(1)=0.2;p(2)=-0.45;q(2)=-0.15;p(3)=-0.4;q(3)=-0.05;p(4)=-0.5;q(4)=0.0; %置迭代次數(shù)k=0，并計(jì)算節(jié)點(diǎn)功率的不平衡量 ， 。k=0;precision=1;k,delt,uN1=4;while precision0.00001 u(4)=1.1; delt(5)=0;u(5)=1.06;for m=1:N1 for n=1:N1+1pt(n)=u(m)*u(

14、n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);endpp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);endpp,qq%/%求取雅可比矩陣元素 （m=n 時(shí)）。for m=1:N1 for n=1:N1+1h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);n0(n)=-

15、u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);endH(m,m)=sum(h0)-u(m)2*(G(m,m)*sin(delt(m)-delt(m)-B(m,m)*cos(delt(m)-delt(m);N(m,m)=sum(n

16、0)-2*u(m)2*G(m,m)+u(m)2*(G(m,m)*cos(delt(m)-delt(m) +B(m,m)*sin(delt(m)-delt(m);J(m,m)=sum(j0)+u(m)2*(G(m,m)*cos(delt(m)-delt(m)+B(m,m)*sin(delt(m)-delt(m);L(m,m)=sum(l0)+2*u(m)2*B(m,m)+u(m)2*(G(m,m)*sin(delt(m)-delt(m) -B(m,m)*cos(delt(m)-delt(m);end for m=1:N1-1 JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m

17、)=N(m,m); JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m); end for m=N1:N1 JJ(2*m-1,2*m-1)=H(m,m); end %/%求取雅可比矩陣元素 （m不等于n時(shí)） for m=1:N1 for n=1:N1 if m=n elseH(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n);N(m,n

18、)=-J(m,n);L(m,n)=H(m,n); end endendfor m=1:N1-1 %求前六行六列元素(非對(duì)角) for n=1:N1-1 if m=n elseJJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);JJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=L(m,n); end endend for m=N1 for n=1:N1-1 %求取第七行的元素JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n); endend for n=N1 for m=1:N1-1 %求的第七列元素JJ(2

19、*m-1,2*n-1)=H(m,n);JJ(2*m,2*n-1)=J(m,n); endend %/%解修正方程式，由 ， 和 計(jì)算電壓修正量 和 。for m=1:N1-1 PP(2*m-1)=pp(m);PP(2*m)=qq(m);endfor m=N1 PP(2*m-1)=pp(m);enduu=-inv(JJ)*PP;precision=max(abs(uu);uu%/ %若結(jié)果不收斂，執(zhí)行下列語(yǔ)句。 for n=1:N1-1 delt(n)=delt(n)+uu(2*n-1); u(n)=u(n)+uu(2*n);end for n=N1 delt(n)=delt(n)+uu(2*n

20、-1); end k=k+1;k,delt,uend%/ %若結(jié)果收斂，計(jì)算各節(jié)點(diǎn)電壓，平衡節(jié)點(diǎn)功率，PV節(jié)點(diǎn)功率和線路功率 。for n=1:N1+1 U(n)=u(n)*(cos(delt(n)+j*sin(delt(n);endfor m=1:N1+1 I(m)=Y(5,m)*U(m);endS5=U(5)*sum(conj(I); %平衡節(jié)點(diǎn)功率%PV節(jié)點(diǎn)功率for n=1:N1+1q4(n)=u(4)*u(n)*(G(4,n)*sin(delt(4)-delt(n)-B(4,n)*cos(delt(4)-delt(n);end Q4=sum(q4)%線路功率 afor m=1:N1+

21、1 for n=1:N1+1 S(m,n)=U(m)*(conj(U(m)*conj(d(m,n)+(conj(U(m)-conj(U(n)*conj(-Y(m,n); endend%/%顯示運(yùn)行結(jié)果（至結(jié)束）。YJJSBppqquuUk Q4S56 程序執(zhí)行結(jié)果及分析Y = 13.0250 -34.0920i -1.6670 + 5.0000i -3.3330 + 6.6670i -2.7500 + 8.2500i -5.0000 +15.0000i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.75

22、00i -3.3300 + 6.6670i -10.0000 +30.0000i 14.5800 -40.4170i -1.2500 + 3.7500i 0 -2.7500 + 8.2500i 0 -1.2400 + 3.7500i 3.7400 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500ik = 0delt = 0 0 0 0 0u = 1.0000 1.0000 1.0000 1.1000 1.0600pp = 0.5000 -0.3750 -0.2750 -0.6364qq = 2.7500 0.

23、0750 0.3250 -0.4125uu = -0.0648 0.1115 -0.0882 0.0757 -0.0918 0.0801 -0.0876k = 1delt = -0.0648 -0.0882 -0.0918 -0.0876 0 u = 1.1115 1.0757 1.0801 1.1000 1.0600pp = -0.1013 0.0261 0.0501 0.0208qq = -0.3710 -0.0269 -0.0454 0.8566uu =0.0031 -0.0114 0.0054 -0.0089 0.0056 -0.0090 0.0016k = 2delt = -0.06

24、17 -0.0828 -0.0861 -0.0860 0u = 1.1002 1.0668 1.0711 1.1000 1.0600pp = -0.0162 0.0013 0.0017 0.0047qq = -0.0428 0.0018 -0.0013 0.7061uu = -0.0000 -0.0013 0.0000 -0.0007 0.0000 -0.0008 -0.0001k = 3delt =-0.0617 -0.0828 -0.0861 -0.0860 0u = 1.0989 1.0661 1.0703 1.1000 1.0600pp = -0.0016 0.0002 0.0002

25、0.0005qq = -0.0042 0.0004 0.0001 0.6909uu = 1.0e-003 *-0.0004 -0.1190 -0.0004 -0.0475 -0.0005 -0.0532 -0.0010k = 4delt = -0.0617 -0.0828 -0.0861 -0.0861 0u = 1.0987 1.0660 1.0703 1.1000 1.0600pp = 1.0e-003 * -0.1523 0.0199 0.0209 0.0433qq = -0.0004 0.0001 0.0000 0.6896uu = 1.0e-004 *-0.0000 -0.1073

26、-0.0000 -0.0295 -0.0000 -0.0351 -0.0001k = 5delt =-0.0617 -0.0828 -0.0861 -0.0861 0u = 1.0987 1.0660 1.0703 1.1000 1.0600pp =1.0e-004 * -0.1416 0.0200 0.0231 0.0377qq = -0.0000 0.0000 0.0000 0.6895uu =1.0e-006 * -0.0000 -0.9640 -0.0000 -0.1828 -0.0000 -0.2303 -0.0001k = 6delt =-0.0617 -0.0828 -0.086

27、1 -0.0861 0u = 1.0987 1.0660 1.0703 1.1000 1.0600Q4 =-0.6895Y = 13.0250 -34.0920i -1.6670 + 5.0000i -3.3330 + 6.6670i -2.7500 + 8.2500i -5.0000 +15.0000i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i -3.3300 + 6.6670i -10.0000 +30.0000i 14.5800 -40.4170i -1.2500 + 3.7500

28、i 0 -2.7500 + 8.2500i 0 -1.2400 + 3.7500i 3.7400 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500iJJ = -40.9565 -15.9241 5.8963 1.8286 7.9334 3.7268 10.0491 15.5240 -41.3566 -1.8286 5.8963 -3.7268 7.9334 -3.0796 5.8139 2.0756 -44.1859 -14.2290 34.2658 11.2955 0 -2.0756 5.8139 15.1

29、290 -43.8859 -11.2955 34.2658 0 7.7421 4.1061 34.1899 11.5231 -46.3468 -16.3011 4.4148 -4.1061 7.7421 -11.5231 34.1899 17.1011 -46.2468 -1.4718 9.8870 3.5657 0 0 4.4149 1.4597 -14.3020S = 0 0.1839 + 0.1398i 0.2968 + 0.1152i 0.5722 + 0.9064i -0.8529 + 1.0305i -0.1812 - 0.1319i 0 0.0687 - 0.1735i 0 -0

30、.3375 + 0.1554i -0.2917 - 0.1052i -0.0683 + 0.1744i 0 -0.0400 - 0.1193i 0 -0.5407 - 0.8120i 0 0.0407 + 0.1226i 0 0 0.8825 - 0.9416i 0.3472 - 0.1262i 0 0 0 B = -34.0920 5.0000 6.6670 8.2500 15.0000 5.0000 -38.7500 30.0000 0 3.7500 6.6670 30.0000 -40.4170 3.7500 0 8.2500 0 3.7500 -11.2500 0 15.0000 3.

31、7500 0 0 -18.7500pp = 1.0e-004 * -0.1416 0.0200 0.0231 0.0377qq = -0.0000 0.0000 0.0000 0.6895uu =1.0e-006 * -0.0000 -0.9640 -0.0000 -0.1828 -0.0000 -0.2303 -0.0001U = 1.0966 - 0.0677i 1.0624 - 0.0881i 1.0663 - 0.0920i 1.0959 - 0.0945i 1.0600 k = 6Q4 = -0.6895S5 = 1.2297 - 1.0678i d = 0 0 0 0.2750 -

32、 0.8250i 0 0 0 0 0 0 0 0 -0.2500 + 0.7500i 0 0 0 Y =13.0250 -34.0920i -1.6670 + 5.0000i -3.3330 + 6.6670i -2.7500 + 8.2500i -5.0000 +15.0000i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i -3.3300 + 6.6670i -10.0000 +30.0000i 14.5800 -40.4170i -1.2500 + 3.7500i 0 -2.7500

33、+ 8.2500i 0 -1.2400 + 3.7500i 3.7400 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500ik = 0delt = 0 0 0 0u = 1.0000 1.0000 1.0000 1.1000 1.0600pp = 0.5000 -0.3750 -0.2750 -0.6364qq = 2.7500 0.0750 0.3250 -0.4125uu = -0.0648 0.1115 -0.0882 0.0757 -0.0918 0.0801 -0.0876k = 1delt =-0

34、.0648 -0.0882 -0.0918 -0.0876 0u =1.1115 1.0757 1.0801 1.1000 1.0600pp = -0.1013 0.0261 0.0501 0.0208qq = -0.3710 -0.0269 -0.0454 0.8566uu = 0.0031 -0.0114 0.0054 -0.0089 0.0056 -0.0090 0.0016k =2delt =-0.0617 -0.0828 -0.0861 -0.0860 0pu = 1.1002 1.0668 1.0711 1.1000 1.0600p =-0.0162 0.0013 0.0017 0

35、.0047qq = -0.0428 0.0018 -0.0013 0.7061uu =-0.0000 -0.0013 0.0000 -0.0007 0.0000 -0.0008 -0.0001k =3delt = -0.0617 -0.0828 -0.0861 -0.0860 0u = 1.0989 1.0661 1.0703 1.1000 1.0600pp = -0.0016 0.0002 0.0002 0.0005qq = -0.0042 0.0004 0.0001 0.6909uu =1.0e-003 * -0.0004 -0.1190 -0.0004 -0.0475 -0.0005 -

36、0.0532 -0.0010k = 4delt =-0.0617 -0.0828 -0.0861 -0.0861 0u = 1.0987 1.0660 1.0703 1.1000 1.0600pp =1.0e-003 *-0.1523 0.0199 0.0209 0.0433qq =-0.0004 0.0001 0.0000 0.6896uu =1.0e-004 * -0.0000 -0.1073 -0.0000 -0.0295 -0.0000 -0.0351 -0.0001k = 5delt =-0.0617 -0.0828 -0.0861 -0.0861 0u =1.0987 1.0660

37、 1.0703 1.1000 1.0600pp = 1.0e-004 * -0.1416 0.0200 0.0231 0.0377qq = -0.0000 0.0000 0.0000 0.6895uu = 1.0e-006 * -0.0000 -0.9640 -0.0000 -0.1828 -0.0000 -0.2303 -0.0001k = 6delt =0.0617 -0.0828 -0.0861 -0.0861 0u = 1.0987 1.0660 1.0703 1.1000 1.0600Q4 = -0.6895Y =13.0250 -34.0920i -1.6670 + 5.0000i

38、 -3.3330 + 6.6670i -2.7500 + 8.2500i -5.0000 +15.0000i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i -3.3300 + 6.6670i -10.0000 +30.0000i 14.5800 -40.4170i -1.2500 + 3.7500i 0 -2.7500 + 8.2500i 0 -1.2400 + 3.7500i 3.7400 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0

39、0 6.2500 -18.7500iJJ =-40.9565 -15.9241 5.8963 1.8286 7.9334 3.7268 10.0491 15.5240 -41.3566 -1.8286 5.8963 -3.7268 7.9334 -3.0796 5.8139 2.0756 -44.1859 -14.2290 34.2658 11.2955 0 -2.0756 5.8139 15.1290 -43.8859 -11.2955 34.2658 0 7.7421 4.1061 34.1899 11.5231 -46.3468 -16.3011 4.4148 -4.1061 7.742

40、1 -11.5231 34.1899 17.1011 -46.2468 -1.4718 9.8870 3.5657 0 0 4.4149 1.4597 -14.3020S =0 0.1839 + 0.1398i 0.2968 + 0.1152i 0.5722 + 0.9064i -0.8529 + 1.0305i -0.1812 - 0.1319i 0 0.0687 - 0.1735i 0 -0.3375 + 0.1554i -0.2917 - 0.1052i -0.0683 + 0.1744i 0 -0.0400 - 0.1193i 0 -0.5407 - 0.8120i 0 0.0407

41、+ 0.1226i 0 0 0.8825 - 0.9416i 0.3472 - 0.1262i 0 0 0 B =-34.0920 5.0000 6.6670 8.2500 15.0000 5.0000 -38.7500 30.0000 0 3.7500 6.6670 30.0000 -40.4170 3.7500 0 8.2500 0 3.7500 -11.2500 0 15.0000 3.7500 0 0 -18.7500pp =1.0e-004 * -0.1416 0.0200 0.0231 0.0377qq = -0.0000 0.0000 0.0000 0.6895uu = 1.0e

42、-006 *-0.0000 -0.9640 -0.0000 -0.1828 -0.0000 -0.2303 -0.0001U = 1.0966 - 0.0677i 1.0624 - 0.0881i 1.0663 - 0.0920i 1.0959 - 0.0945i 1.0600 k = 6Q4 = -0.6895S5 = 1.2297 - 1.0678i迭代過(guò)程中各節(jié)點(diǎn)功率的不平衡量K P1Q1P2Q2P3Q3P4Q400.50002.7500-3.7500.0750-0.27500.3250-0.6375-0.41251-0.1013-0.37110.0261-0.02690.0501-0.

43、04540.02080.85652-0.1620-0.04280.00130.00180.0017-0.00130.00470.70593-0.0016-0.0042-0.00020.00040.00020.00010.00050.69074-0.152.3*e-3-0.0004*e-30.0199*e-30.0001*e-30.0209*e-30.0000*e-30.0433*e-30.6894*e-35-0.1416*e-4-0.0000*e-40.0200*e-40.0000*e-40.0023*e-40.0000*e-40.0378*e-40.6893*e-4 迭代過(guò)程中的各節(jié)點(diǎn)電壓的修正量k