潮流计算的运行结果

系统接线图：S6c0j2o程序：%开始clccleardisp（节点总数为:）;N=5disp （平衡节点为:）;1disp（PQ节点为：）;JD=2,3,4,5Y=6.25-18.75i-5+15i-5+15i10.834-32.5i-1.25+3.75i-1.667+5i0-1.667+5i-1.25+3.75i00;-1.667+5i-1.667+5i-2.5+7.5i;12.917-38.75i-10+30i0;-10+30i12.917-38.75i-1.25+3.75i;-2.5+7.5i-1.25+3.75i3.75-11.25idisp(导纳矩阵为:Y=);%求出导纳矩阵ei=1.06;1;1;1;1;fi=0;0;0;0;0;P=0;,-0.2;0.45;0.40;0.60;Q=0;-0.20;0.1 5;0.05;0.10;disp(Y);wucha=0.0001;k=0;G=real(Y);B=imag(Y);N1=4;while wucha0.00001for m=2:5pp=0;qq=0;for n=1:5pp=pp+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n);%Pi(0)qq=qq+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)-ei(m)*(G(m,n)*fi(n)+B(m,n)*ei(n);I(m)= (pp-qq*i)/conj(ei(m)+fi(m)*i);%节点注入电流 Ii (0)endDP(m)=P(m)-pp; %有功功率的不平衡量DQ(m)=Q(m)-qq; %无功功率的不平衡量H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+imag(I(m);N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+real(I(m);J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+real(I(m);L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)-imag(l(m);%公式 P131JJ(2*m-3,2*m-3)=H(m,m);JJ(2*m-3,2*m-2)=N(m,m);JJ(2*m-2,2*m-3)=J(m,m);JJ(2*m-2,2 *m-2)=L(m,m);endfor m=2:Nfor n=2:Nif n=melseH(m,n)=-B(m,n)*ei(m)+G(m,n)*fi(m);N(m,n)=G(m,n)*ei(m)+B(m,n)*fi(m);J(m,n)=-B(m,n)*fi(m)-G(m,n)*ei(m);L(m,n)=G(m,n)*fi(m)-B(m,n)*ei(m);JJ(2*m-3,2*n-3)=H(m,n);JJ(2*m-3,2*n-2)=N(m,n);JJ(2*m-2,2*n-3)=J(m,n);JJ(2*m -2,2*n-2)=L(m,n);%公式 P130endendenddisp(雅克比矩阵JJ:);disp(JJ);for m=2:NDW(2*m-3)=DP(m);DW(2*m-2)=DQ(m);%形成 DP 和 DQ 不平衡量的矩阵endDY=JJDW;%解修正方程式wucha=max(abs(DY);disp (第M次修正方程的解DY:);disp(DY);for n=2:Nfi(n)=fi(n)+DY(2*n-3);ei(n)=ei(n)+DY(2*n-2);%计算各节点电压新值endu=ei+fi*i;disp(节点电压的第C(k)次近似值:);disp(u);disp(各点的电压实部ei(单位:V)为(节点号从小到大排列)：)；disp(ei);disp (各点的电压虚部fi(单位:V)为(节点号从小到大排列)：)；disp(fi)；k=k+1;disp(迭代次数:)；disp(k);endfor n=1:N1+1I(n)二Y(l,n)*u(n);%(公式 P139)enddisp(平衡节点的功率)；Sl=u(l)*sum(conj(I)%计算平衡节点的功率for m=1:N1+1for n=1:N1+1S(m,n)二u(m)*(conj(u(m)-conj(u(n)*conj(-Y(m,n);%计算 m-n 各支路功率S2(n,m)二u(n)*(conj(u(n)-conj(u(m)*conj(-Y(n,m);%计算 n-m 各支路功率DS(m,n)=S(m,n)+S2(n,m);%计算各支路功率的损耗(公式P135 4-51a)endenddisp(各支路 m-n 功率)；disp(S)disp(各支路 n-m 功率)；disp(S2)disp (各支路功率损耗)；disp(DS) %结束运行结果：节点总数为：N 二平衡节点为：ans -PQ节点为：JD 二2345二-5. ODDO +1&. DDDDi-1. E5DD + 3. 7&DDi导纳袒阵为汽二6. E5DD -IE. 75DDi0DDD +15. DDDDi 飞-5. DODD +1&. DDDDi10. B34D -3Z. &DDDi-1. 6670 + &. DDDDi-1. 6670-1. E50D + 3. 7&DDi-1. 6670 + 5. DDDDiIE. 5170 -3B. 7&DDi+ &. DDDDi-E. 50DD + 7. &DDDiDODD +15. DDDDi10. B34D -3E. 5DD0i-1. 6670 + 5. DDDDi-1. 6670 + 5. DDDDi-E. &ODD + 7. 5DD0i-10. DODD +30. DDDDiD-1. E&OD + 3. 7500i-1. 6670 + 5. DDDDiIE. 5170 -3B. 7500i-10. DODD +30. DDDDiDD-1. 6670 + 5. DDDDi-10. DODD +30. DDDDiIE. 5170 -3B. 7&DDi-1. E50D + 3. 7&DDiD-1. 6670 + 5. DDDDi-10. DODD +30. DDDDiIE. 5170 -3B. 7500i-1. E&OD + 3. 7500i0-E. 5000 + 7. 5DDDi0-1. 500 + 3. 7&DDi3. 7500 -11. E&DDi0-E. 5000 + 7. 5D00i0-1. E&00 + 3. 7500i3. 7&00 -11. E5D0i雅克比走阵JJ:33, 400010. 5340-5. 0000-1. 6670-5. 0000-1. 6670-7. 5000-E. 5000-11. 134031. 60001. 6670-5. 00001.6670-5. 0000E. 5000-7. 5000-5. 0000-1. 66703B. 9750IE. B4E0-30.0000-10. 0000001. 6670-5. 0000-IE. 99E03B. 5E5010.0000-30. 000000-5. 0000-1. 6670-30. 0000-10. 00003B.75001E. 9170-3. 7500-1. E5001. 6670-5. 000010. 0000-30. 0000-1E.91703B. 75001. E500-3. 7500-7. 5000-E. 500000_3,7500-1. E50011. E5003. 7500E. 5000-7. 5000001.E500-3, 7500-3, 750011. E500第間次修正方程斷解D. 04730. 0E47D. DEB3D. 11E30. 0BEED. 113BD. 10760. UB3节点电压的第f S蛮近忖值:1. D6DD1, DB47 4 D. 047J11. 113 4 D. D&閃i1, H36 4 D.阳亞丄1. 11B3 4 D. 1D761吝直的电压实鄱畀单位：心为节点号从小到大排列）：1.06001,DB471. 1131. 11361. 11E3善点的电压虚IfHil单位汕）为I节直号从小到大排列:o0, 04730. 0B63D, DDES0. 1076迭代E1雅克比走阵JJ：35, B7479. 99T4-5, 5DE3-1. 571750E1-1. 5717-E. 2535-E. 3570-10.431715. 65&51.5717-5. 50211. 5717-5. 50232. 3570-B. 25357D&2-1.4EE744.033311.4624-34. 2310-B. &330001.4E7-5. 70&E-10. 5B4444. 355EE. 533&-34.007E1&-1,3D&E-34,3EB&-B, 36BE44, EB1&11, 101-4, E&11-1, 04611.3S&Z-5. 7E15B. 36EB-34. 3E凹-10. 410144. 39151. 0461-4. Z911-B. 65&4-1.BBEBD0-4. 3E吃-D. 0B44IE. BB463. 5B311.0BEG-G. 65&4D00. 9944-4. 32B2-2. 3G3113. DG46第間次修疋方程的解DY:0. 000E-0. 0116-0. 0010-0. 0171-0. 0015-0. 0177-0. 00E7-0. 019B节点电压的第CG庶近似值：1. 06DD1. 0731 + 0. 0475i1. 0B51 + 0. 0B53i1. 0B5B + 0. O0O7i1. 0SB5 + 0. 1049i諄点的电压奕韶卵单位：心为节点号从小到丈排列）：1. D6D01. 07311. D0&11. D95E1. 0SB5善点的电压虚揶单位汕対节点号尿力到女排列：00. D47E0. DB53D. 0S070. 1D4B迭代诙数:2雅克比走阵JJ：3&. &6G1&. GG64-&. 444G-1. &11-&. 444G-1. &11-G. 1672-2, 3262-10. E7533-5. E1411. 5511-5. 444B1. 5511-5. 444BE. 3E6E-B. 167E-5. 617B-1. 399143. 43E111. E600-33. 7065-E. 39開001.潮1-&. 617G-10.421243. 6431G.旧開-33. 706&00-5. 6303-1. 3731-33. 7B17-B. E36643. 61B911. 0061-4. EEE7-1. 0ED61. 3731-5. 6303B. E366-33. 7B17-10. E7E543. 65061. 0E96-4. EEE7-E. 50U-1.呢规00-4. E507-C.軒軒IE. 71E43. 4E571. 9593-B. 5013000. 9797-4. E507-E. 3RR31E. 7917第屈次慘正另程断解DY二1. 0&-0D3 申0. 0017-0. 1747-0. DOB0D, Z6Z1-0. D14D-0. 2722-0. D29D-0. 314B节点电圧的第C(k)伙近似值:1. 06001. 0720 + 0. D475i1. 0049 + 0. DE53i1. 0955 + 0. DSD7i1. DDES + 0. lD49i各点的电压实鄱訂（单位:：门为（节点号从小到大排列）：1. 06001. D7E&1. DS4&1. D05&1, 09此各点的电压虚鄱门（单位:：门为（节点号从小到大排列）：D0. 047&D. DB53D. D0D70, 1049迭代次数：3雅克比电阵JJ：35, 5633艮 EE46-5. 443&-1. 55013-5. 443&-1. 55013-E. 165&-E. 3E57-10. 273235, E0771. 550E-5. 44391. 550E-5. 4439E. 3257-E. 1659-5. 6165-1. 39R743. 4E3011. 2566-33, 69E6-E. 3900001. 39R7-5. 6165-10. 41E343. 6317E. 3900-33, 69E600-5. 6EE9-1. 3727-33, 7734-E. 234343. 60EE11. 00E7-4. 2217-1. D2931. 3727-5. 6EE9E. 2343-33, 7734-10. E70043. 63941. D293-4 EE17-E. 49E9-1. 95E700-4. 2495-0. 979412. 70993. 4EE11. 95E7-E. 49E9000. 9794-4. 2495-E. 3EE1IE. 7E69第M枚侈正方程的解DY:1.08-007 辛0. 0009-0. 4141-0. 0037-0. 6240-0. 0187-0. 6&12-0. 0485-0. 7B96爷点电圧的第C(k)枚近似值:1. 06001.0729 0. 0475114949 0. 0B53i1.0&55 0. 0907i1.0&B2 0.1049i善点的电压夹郤（单位:町为（搐点呂从小到大排列）:1. 06001. 072&1. 034&1. 09551. 0982善点的电圧戒郤单位:心次（怖血号从小列大排列）:00. 047&0. 0853-0. 09070. 1049迭代次数：451 二-1. 2097 O.OZOBi各支路oi-n功率0-0. 824510. 0464i-0. 3852-0. 0256i000. 8367 - 0. 0099i0-0. 2394-0. 060Bi-0. 2700 - 0. 05621-0. 5272 -0. 0731i0. 3958 0. 0574i0. 242610. 0703i0-0. 1884 0. 02231000. 274010. 06B0i0. 1BB7-0. 021410-0. 0627 +0. 0033i00. 537010. 1025100. 0630 - 0. 002510茗支路n-oi功率0-0. B245+ 0. 04641-0. 3B52 -0. Q256i000. B367 - 0. DD991D-0. 2394 -D. 060Bi-0. 2700 -0. D5621-0. 5272 - D. 073110. 395B + 0. D57410. 2426+ 0. 070310-0. 1BB4 +0. D223i000. 2740+ 0. OBECu0. 1BB7 -D. 021410-0. 0627 + 0. 0033100. 5370+ 0. 1025100. 0630 -0. 00Z510CD. 0121+ 0. 036410. 0106 +D. D31Bi000. 01E1 + 0. 03-64100. 0032 +0. CD95i0. 0040 + D. 0119iD. DD9B + 0. 029510. 0106 + 0. 03-1 BiD. DD32+ 0. 0095100. 0003- + 0. OOCSi0CD. DD40+ 0. 011510. 0003 +0. CDD9i0D. DDD3 + 0. 00081CD. DD9B+ 0. 0295100. 0003- + 0. 000B10