潮流计算(matlab)实例计算

《潮流计算(matlab)实例计算》由会员分享，可在线阅读，更多相关《潮流计算(matlab)实例计算（11页珍藏版）》请在装配图网上搜索。

1、潮流计算实例潮流例题：根据给定的参数或工程具体要求（如图）,收集和查阅资料；学习相关软件（软件自选:1本设计选择Matlab进行设计）。0 4+j(),O50.040 ifO.O8+j0.24-(0200.2)0.6+jOJ不统接旗盟中节点1为平衡节点.节点，3、4+ 5为PQ节点;2 .在给定的电力网络上画出等值电路图。3 .运用计算机进行潮流计算。4 .编写设计说明书。一、设计原理1.牛顿-拉夫逊原理牛顿迭代法是取x0之后，在这个基础上，找到比 x0更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0的单根附近

2、具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然 后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流 计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写 成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额 定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不潮流计算实例平衡量(未知的)构成了误差方程

3、，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡 量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的 电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。牛顿一拉夫逊迭代法的一般步骤：(1)形成各节点导纳矩阵 丫。(2)设个节点电压的初始值 U和相角初始值e还有迭代次数初值为 0。(3)计算各个节点的功率不平衡量。(4)根据收敛条件判断是否满足，若不满足则向下进行。(5)计算雅可比矩阵中的各元素。(6)修正方程式个节点电压(7)利用新值自第(3)步开始进入下一次迭代，直至达到精度退出循环。(8)计算平衡节点输出功率和各线路功率

4、2 .网络节点的优化1)静态地按最少出线支路数编号这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按 出线支路数有少到多的节点顺序编号。当由n个节点的出线支路相同时，则可以按任意次序对这n个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也2)动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出 来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连 的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出

5、线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动 态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目 的变动情况。3 . MATLA叫程应用Matlab是“Matrix Laboratory ”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、 系统控制、优化和图形显示等应用程序。由于使用Matlab编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输 出结果和精美的图形拷贝，所以它的确为一高效的科研助手。二、设计内容1.设计流程图启

6、动-；输入原始数据潮流计算实例112.程序clear;clc,其中1-4号为PQ节点，5号为平%重新编号，把原题中的节点1,2,3,4,5重新依次编号为5,1,2,3,4衡节点y=0;%输入原始数据，求节点导纳矩阵y (1,2)=1/(0.06+0.18i); y (1,3)=1/(0.06+0.18i); y (1,4)=1/(0.04+0.12i);y(1,5)=1/(0.02+0.06i);y(2,3)=1/(0.01+0.03i);y(2,5)=1/(0.08+0.24i);y(3,4)=1/(0.08+0.24i);y(4,5)=0;for i=1:5for j=i:5 y(j,i)

7、=y(i,j);end end Y=0;%求互导纳 for i=1:5 for j=1:5 if i=j Y(i,j)=-y(i,j); end end end %求自导纳 for i=1:5 Y(i,i尸sum(y(i,:); endY %Y为导纳矩阵 G=real(Y);B=imag(Y);%原始节点功率 S(1)=0.2+0.2i;S(2)=-0.45-0.15i;S(3)=-0.4-0.05i;S(4)=-0.6-0.1i;S(5)=0;P=real(S);Q=imag(S); %就初值 U=ones(1,5);U(5)=1.06;e=zeros(1,5);ox=ones(8,1);f

8、x=ones(8,1);count=0%计算迭代次数while max(fx)1e-5 for i=1:4for j=1:4H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0; endendfor i=1:4for j=1:5oP(i)=oP(i)-U(i)*UO)*(G(i,j)*cos(e(i)-eQ)+B(i,j)*sin(e(i)-e(j); oQ(i)=oQ(i)-U(i)*UO)*(G(i,j)*sin(e(i)-eO)-B(i,j)*cos(e(i)-e(j); endoP(i)=oP(i)+P(i); oQ(i)=oQ(i)+

9、Q(i);endfx=oP,oQ;%求雅克比矩阵i=j 时候求H,N,M,L 如下：for i=1:4for j=1:4if i=j H(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j); N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);L(i,j)=H(i,j);M(i,j)=-N(i,j);end end end H,N,M,L %当i=j时H,N,M,L 如下：for i=1:4 for j=1:5 if i=jH(i,i)=H(i,i)+U(i)

10、*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i, j)*cos (e(i)-e(j); N(i,i)=N(i,i)-U(i)*U(j)*(G(i, j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j); M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j); L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j); end endN(i,i)=N(i,i)-2*(U(i)A2*G(i,i);L(i,i

11、尸L(i,i)+2*(U(i)2*B(i,i);endJ=H,N;M,L %J为雅克比矩阵ox=-(inv(J)*fx);for i=1:4oe(i)=ox(i); oU(i)=ox(i+4)*U(i);endfor i=1:4e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);endcount=count+1;endox,U,e,count%求节点注入的净功率 i=5;for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j)+P(i);Q(i)=U(i)*UQ)*(G(i,j)*sin(e(i)-eQ)

12、-B(i,j)*cos(e(i)-eO)+Q(i); endS(5)=P(5)+Q(5)*sqrt(-1);S%求节点注入电流I=Y*U3.运行结果Y值：Command Window4 i) Hev; to MATLA.fi? .Vatcli this 丫3.see Demes, or read Starts。.10.S333-32.5000i-1.6667+ 5., OOOOi-1.666： + 5.00001-2.5000 + 7.5000i-5.0000 +15. OOQOi-1.6657+ 5.OOOQi12.9167-38-7500i-10.0000 +3Q. 000010-1.25

13、00 + 3.7500i-L6667+ 5.00001-10. 0000+30.OOQOi12.9167 38.75001-1.2500 + 3. 750010-2.5000+ 7.5000iG-1.2500 + 3. 7500i3.7500 -11. 2500i0-5.0000+15* OOQOi-1. 2500+ 3. 7500i0Q6.2500 -18.75001迭代过程：count =QH =05, OQOQ5. 00007, 50005.DOOO03d QQQQO5. 000030.000003. 75007.500003. 7500QQ.66571,66672. 5000I- 66

14、67010+ 000001. 6067IO.000001. 25002.500001.2600QM =0-1.6667-1,6667-2.5000-1,S6670一I口.00000-10,oaoo07. 2500-2.6000-1. 2500O03. 00005.00007. 50005. 0000030.000005.000030. aooo03. 75007.500003. 75000J =-33.4000E.00005.00007. 5000-10r53331.6667L. 66672. 50005. 0000-384975030. 000001,6667-12.841710,00000

15、5. 000030.aooo-33. 75003. 75001.66S7ia.oooo-12. 91671. 25007. 500003. 7500-11. 25002.500001. 2500一工 750011.1333-1.5567-1.6667-2. 5000-3U60005.00005.0000二 5000-1,566712+9917-10.000005+0000-38-525030.00000-L 6667-10. 000012,9167-L 25003.000030.0000-38.7500X 7500-2. 50000-1. 25003. 75007.500003.7500-11

16、. 250005. 35995. 30228. O3S6S, 三芝少2o30.9520O5,203830.829S03. 85767, 722303. 3133Q一aU 55T21B 5235Z, 15111.970110.11370-1,10.48031, 22033. 1Q2301. 33S4L00-1-5572-1.S235-2. 1S 1 1一1.970101011370 I. 9985-ia_ 4EO3a-1. 2203一 3. 1029-1.33840L =1Q5. 35995,35228. 03965. 2222030,9520Q5. 2Q3B30.329603,断67. 722

17、303. 81S3Q-35. 06485. 35995.35223.0396-12.0310L. 55721. 52352. 151L5. 2222-40. 079430,95200I,9701-12,862710.113705, 203830.829879 89133,85781.998510, 4B03-13.8682L 22037. 722303.81S3-11.54063.102901.3384-3. 199IL, 5373-1.5572-1, 5235-2.1511-35. 64005. 35995.36228.0396-1, 97013,7720-10,113705.2222-39

18、, S24730.95200-1.9985-10.480313.6991-1.22035.203830.629379. 31053.8576-3. 10290-1.33944.4412:.722303. 81E3-11.3818O5. 290 45. 291274 9246lessOO5.139 7SO_ 4.26T口3, S01ST.=17O3. 76-14ON 二a 1.S333 1. 9550 3P 口胃百口1.5435 1O.3355 口1. 51129. 9S7S O1.31712.13&8 CL 2050 QM =Q一 1 7 543(5-1,Cl1Z2, 1 35S-1-93&S

19、O-S. 9STAa-1.9658-1O.3355a-1.2050-3, O4B3O-1.3171c05.29045.29127.92465.15S5030.542C05.139730,4267Q3.8018工 621703, 75440一34.71635.2S045. 2B127, 9246-11.33941.5436L 51 122. L36B5,5S5-39.584930, 5420.9393-12, BS509, S87B05, 39730.4267-39.36813, BOLS,965SQ.3355-12,7059).20507. 621703. 7644-1 L, 38613. 0

20、45301. 3171-3. 911-43S4-1. 5436-U 5112-2.136S-35. 1 L735.29045. 29127. 9246-1-9393】工5950-9, 537B05. 1585 39, 23&330.5426Q-1.965g-1。. 335513.5062-1.20505. 139730.425二-39.26833. 8018-3.0453Q-1. 31714. 35237, 21703.7644-11, 186605.28995- 29077, 52375-1580030.539505. 139230.423503. 80137.6210Q3.7640QQ1,

21、54361.51112.136S1.9391Q9.98S701. 965510-33131B 20193. 044S1. 3169a0一L 5436-1. 5111-2- 1 368一1.S3S10一9. 9857a-1,9655-10.33430-1, 2049一工 04480-1.31S90L =05.28995, 29077, 92375. 1560030.5395Q5-130230.423S03. 801 37.621003. 7540a-34.71365.2SS95. 29077. 9237-1 1.83791.54361.51112.13685.15SQ-39.581230, 53

22、9501. 9391-12. 69379.9867Q5.13923QH23G-39.3642180131. 965510.3343-12. 7047L 214g7.021003. 764。-11.38503.044301.3169*3, 161711,1379-L 5436-I. 5111-3.1366-35.11365.28995.29077.9237-1.939113. 5S37-9.986705.1530-39.2S1230. 53950-1,9655-10.334313, 5047-1. 20495. 139230,4236-39,26423.8013-3. 0448Q-1. 3169

23、4436177.621003. 7640-11.1850电压值:1 0e- OO 小一CL 00040. 08350- 11B3CL 21S T0. 2645一。.3080一 Cl . 3 215-0. 39281, 03651. 008S1.00731, 00161 . 0600O, 0461-0, 0839-O, 0896-0, 1044c：Qurrt =平衡节点注入功率及电流:0.200。+ 0+ 2000i -0.4500 - Q, 1500i -0. 400C - 0. 05001 -0.6000 - 0. lOOOi 1. 2982 + 0, 24451I 二0-0645 Oi1935i-0L0953 + 0.2E581-0. 0565 + 0. 16961-0. 0944 + 0. 283210. 1817 - 0.54521