C语言进行潮流计算.doc

1、电力系统课程设计C 语言潮流计算学院：电气工程班级：电 092 班学号： 0912002020学生姓名：闵 凯2013.3.7电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤，是指在一定的系统结构和运行条件下，确定系统运行状态的计算，也即是对各母线（节点）电压，各元件（支路） 传输电线或功率的计算。 通过计算出的节点电压和功率分布用以检查系统各元件是否过负荷，各点电压是否合理，以及功率损耗等。即使对于一个简单的电力系统，潮流计算也不是一件简单就可以完成的事，其运算量很大，因此如果对于一个大的、复杂的电网来说的话，由于其节点多，分支杂，其计算量可想而知， 人工对其计算也更是难上加难

2、了。特别是在现实生活中，遇到一个电力系统不会像我们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率，而是只知道它的首端电压和末端功率， 更是使计算变的头疼万分。 为了使计算变的简单， 我们就可以利用计算机，用 C 语言编程来实现牛顿 -拉夫逊（ Newton-Raphson ）迭代法，最终实现对电力系统潮流的计算。一用牛顿 -拉夫逊迭代法进行电力系统潮流计算的相关概念1.节点导纳矩阵.?如图所示的电力网络,将节点 i 和 j 的电压用 U i 和 U j 表示 ,它们之间的支路导纳表示为.?yij ,那么有基尔霍夫电流定律可知注入接点I的电流 Ii (设流入节点的电流为正)等于

3、离开节点I 的电流之和 ,因此有j?I ijI?I in ?n.?I iI ij(U iU j )(1-1)j0yijj 0ii?nn?I iUyijyij U(1-2)j 0j0iinyijYiiyijYij 则可将 (1-2)改写为 :如令j0i?n?I iYijU ijI=1,2, ,n.(1-3)j1i上式也可以写为: I=YU(1-4)其中 Y 为节点导纳矩阵，也称为稀疏的对称矩阵，它是nn 阶方阵。对角元 Y ii 称为自导纳，它等于与该节点 I直接相连的所有支路导纳总和；非对角元Y ij（ i j）称为互导纳或转移导纳，它等于连结节点I ，j 支路导纳的负数，且有Yij =Y j

4、i ，当节点 I， j 之间没有支路直接相连时， Y ij =Y ji =0 。电力系统的分析计算中，往往要作不同的运行方式下的潮流计算，如果系统发生变化，如投切一条线路或一台变压器， 由于改变了一条支路的状态或参数只影响该支路两端节点的自导纳和他们之间的互导纳， 因而对每一种运行方式不必重新形成导纳矩阵， 只需对原有导纳矩阵作相应的修改即可。2潮流计算的功率方程在实际的电力系统中，已知的条件往往不是节点的注入电流而是负荷和发电机功率，这些功率一般不随节点的电压变化而变化，而节点的电流则是随电压的变化而变化的，在已知节点导纳矩阵的情况下，必须用已知的节点功率来替代未知的节点注入电流，且因此才能

5、求出节点电压，每一个节点的注入功率方程式为：?IU i U i iPGijQGiPLijQ LiSGiSLiGPiPGiPLi（ 1-5）QiQGiQLi（ 1-6）?n?I iYij U j(1-7)j 1节点注入电流用功率和电压表示为：?(PGiPLi )j (QGi QLi )PijQiI iSi（ 1-8）?U iU iU i功率方程可以表示为：PijQ in?Yij U j（ 1-9）?U ij13节点分类对于有 n 个节点的电力网络，可以列出n 个功率方程，由图可知一个节点有四个变量：注入有功功率Pi，注入无功功率Qi，节点电压幅值 U i 和相角i 。 n 个节点有4n 个变量

6、,但只有 2n 个关系式 ,所以为了使潮流有确定解,必须给定其中 2n个变量。根据给定节点变量的不同，可以有以下三种类型的节点：（ 1）PQ 节点：给定注入功率Pi， Qi，即已知 PGi ，PLi ， QGi， QLi ，待求 U i， i。例如：降压变电所母线（负荷节点） ，固定出力的发电厂母线。（ 2）PV 节点：给定了注入有功功率Pi （PGi ，PLi ）， U i 和 QLi ，待求 QGi（ Qi），i。例如：有一定无功电源的降压变电所母线，有一定储备的发电厂母线。（ 3）平衡节点：给定了 Ui， i和 PLi ， QLi ，待求 PGi ， QGi，即 Pi， Qi，用来平衡全

7、电网的功率，通常在一个独立的电力系统中只设一个平衡节点。4牛顿 -拉夫逊迭代法牛顿 -拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的求解过程。对于一个 n 维非线性方程组：f i ( x1 , x2 ,., xn )yin=1,2,3,n假 定 其 初 值 为 x1(0) ,x2(0), ,xn(0) , 也 即 其 近 似 解 , 它 与 真 值 之 间 的 误 差 为 x1( 0) , x2 ( 0) ,., xn ( 0) 也即各变量与真解之间的修正量。将这 n 个方程式都在初值的附近展开成Taylor 级数且忽略二次项及高次项,则可得修正方程fi (x 1(0)

8、, x 2(0) , , xn(0) )f 1x1( 0).f nxn(0)yi , I=1,2, n. (1-10)x10x n0将修正方程写成矩阵形式:(0)(0)(0)f1.f 1( 0)x n 0y1f1 (x 1, x 2, x n)x10x1:.(1-11)(0)(0)(0)f nf n( 0)ynfn (x 1, x 2, x n).xnx1x n 00f 1.f 1x1x n00其中令 J=. .,称之为雅可比 (Jacobi) 方阵。f n.f nx1x n00它的第 I 行，第 j列交点的元素为第I个函数 f i ( x1 , x2 ,. xn ) 对第 j个变量 xj 的

9、偏导数在点(0)(0)(0)）的值，所以方程组是线性方程，可用于求出x1(0), x2( 0),.,(0 )（ x1,x2, ,xnxn ，从而得到新的近似解，xi(1)xi ( 0)xi (0 )（ 1-12）于是得到一般迭代式：(k)(k)(k)f 1f1.y1 f1 (x 1, x 2, , x n)x1 kx n kx1(k ):.（1-13 ）yn fn (x 1(k) , x 2(k) , x n (k) )f n.f nx1 kx n k于是得到近似解：xn(k )xi (k 1)xi ( k )xi ( k)迭代一直进行到Max|y i-f i(x 1(0),x2(0),xn(

10、0)| 或 Max|（ 1-14）x i (k)|km?增大迭代次数， k-k+1计算各节点电压新值：ei (k 1)ei (k )ei ( k )( k 1)(k )( k)fif if i否停止四用编程方法求解实际问题如图所示的一个电力网络，0.45+j0.150.4+j0.0513460.08+j0.240.01+j0.030.40.180.18j2+6+j0.06+j.20.0000j.+080.020.04+j0.125- (0.2+j0 .2)0.6+j0.1G?已知： U 1 1.06j 0 为定值， 其余四个节点都是PQ 节点，且给定的注入功率分别为：0.400.60j 0.1

11、0S20.20 j 0.20, S30.45 j 0.15, S4j0.05, S5由上图可得相应的节点导纳矩阵Y=6.25-18.75i-5+15i-1.25+3.75i00;-5+15i10.834-32.5i-1.667+5i-1.667+5i-2.5+7.5i;-1.25+3.75i-1.667+5i12.917-38.75i-10+30i0;0-1.667+5i-10+30i12.917-38.75i-1.25+3.75i;0-2.5+7.5i0-1.25+3.75i3.75-11.25i 五程序清单#include #include float divRe(b1,b2,b3,b4)

12、float b1,b2,b3,b4;float a1r;a1r=(b1*b3+b2*b4)/(b3*b3+b4*b4);return(a1r);float divIm(b1,b2,b3,b4)float b1,b2,b3,b4;float a1i;a1i=(b2*b3-b1*b4)/(b3*b3+b4*b4);return(a1i);float mulRe(b1,b2,b3,b4)float b1,b2,b3,b4;float a2r;a2r=b1*b3-b2*b4;return(a2r);float mulIm(b1,b2,b3,b4)float b1,b2,b3,b4;float a2i;

13、a2i=b2*b3+b1*b4;return(a2i);float Max(float a,int n)int i;float max;for(i=0;iai+1)max=ai;ai=ai+1;ai+1=max;return(max);main()int i,j,k,n,km;float eps,sumpi1,sumpi2,sumqi1,sumqi2,max,sumir,sumii,I1r,I1i;float pi05,qi05,detpi5,detqi5,Iir05,Iii05,J088,detsi8,detui8, u88,l88,y8,ui18,H44,N44,J44,L44,ei15,

14、fi15;static float ybr55=6.250,-5.000,-1.250,0,0,-5.000,10.834,-1.667,-1.667,-2.500, -1.250,-1.667,12.917,-10.000,0,0,-1.667,-10.000,12.917,-1.250, 0,-2.500,0,-1.250,3.750;static float ybi55=-18.750,15.000,3.750,0,0,15.000,-32.500,5.000,5.000,7.500, 3.750,5.000,-38.750,30.000,0,0,5.000,30.000,-38.750

15、,3.750, 0,7.500,0,3.750,-11.250;float ei05=1.06,1.0,1.0,1.0,1.0;float fi05=0,0,0,0,0;float pi5=0,0.2,-0.45,-0.4,-0.6;float qi5=0,0.2,-0.15,-0.05,-0.1;k=0;km=6;eps=0.00001;dok+=1;printf(Now start.n);printf(The %d timesn,k);for(i=1;i5;i+)printf(pi%d=%-14.6f,i,pii);sumpi2=0;sumqi2=0;for(i=1;i5;i+)for(j

16、=0;j5;j+)sumpi1=(ei0i*(ybrij*ei0j-ybiij*fi0j)+fi0i*(ybrij*fi0j+ybiij*ei0j); sumpi2+=sumpi1;pi0i=sumpi2;printf(pi0%d=%-13.6f,i,pi0i);sumpi2=0;for(i=1;i5;i+)for(j=0;j5;j+)sumqi1=(fi0i*(ybrij*ei0j-ybiij*fi0j)-ei0i*(ybrij*fi0j+ybiij*ei0j); sumqi2+=sumqi1;qi0i=sumqi2;printf(qi0%d=%-13.6f,i,qi0i);sumqi2=0

17、;for(i=1;i5;i+)detpii=pii-pi0i;detqii=qii-qi0i;printf(detpi%d=%-21.6f,i,detpii);printf(detqi%d=%-21.6fn,i,detqii);for(i=1;i5;i+)Iir0i=divRe(pi0i,-qi0i,ei0i,-fi0i);Iii0i=divIm(pi0i,-qi0i,ei0i,-fi0i);printf(Iir0%d=%-22.6f,i,Iir0i);printf(Iii0%d=%-22.6fn,i,Iii0i);for(i=0;i4;i+)for(j=0;j4;j+)if(i=j) Hij

18、=-ybii+1j+1*ei0i+1+ybri+1j+1*fi0i+1+Iii0i+1;Nij=ybri+1j+1*ei0i+1+ybii+1j+1*fi0i+1+Iir0i+1;Jij=-ybri+1j+1*ei0i+1-ybii+1i+1*fi0i+1+Iir0i+1;Lij=-ybii+1j+1*ei0i+1+ybri+1j+1*fi0i+1-Iii0i+1;else Hij=ybri+1j+1*fi0i+1-ybii+1j+1*ei0i+1;Nij=ybri+1j+1*ei0i+1+ybii+1j+1*fi0i+1;Jij=-ybii+1j+1*fi0i+1-ybri+1j+1*ei0

19、i+1;Lij=ybri+1j+1*fi0i+1-ybii+1j+1*ei0i+1;for(i=0;i8;i+)for(j=0;j8;j+)if(i%2=0&j%2=0) J0ij=Hi/2j/2;else if(i%2=0&j%2!=0) J0ij=Ni/2(j-1)/2;else if(i%2!=0&j%2=0) J0ij=J(i-1)/2j/2;else J0ij=Li/2(j-1)/2;printf( 输出雅可比矩阵：n);for(i=0;i8;i+)for(j=0;j8;j+)printf(%-10.4f,J0ij);for(i=0;i8;i+)if(i%2=0) detsii=de

20、tpi(i+2)/2;else detsii=detqi(i+1)/2;printf(detsi%d=%-11.6f,i,detsii);for(i=0;i8;i+) uii=1.000;for(n=0;n8;n+)for(i=n;i8;i+)lin=J0in;for(j=0;j=n-1;j+)lin-=(lij*ujn);for(j=n+1;j8;j+)unj=J0nj;for(i=0;i=n-1;i+)unj-=(lni*uij);unj/=lnn;for(i=0;i8;i+)yi=detsii;for(j=0;j=0;i-)detuii=yi;for(j=i+1;jn;j+)detuii

21、-=(uij*detuij);for(i=0;i8;i+)printf(detui%d=%-11.6f,i,detuii);for(i=0;i8;i+)if(i%2=0) ui1i=detuii+fi0i/2+1;else ui1i=detuii+ei0(i+1)/2;printf(ui1%d=%-13.6f,i,ui1i);for(i=1;i5;i+)ei1i=ui12*i-1;fi1i=ui12*i-2;for(i=1;i5;i+)printf(ei1%d=%-13.6f,i,ei1i);printf(fi1%d=%-13.6f,i,fi1i);max=Max(detui,8);print

22、f(max=%fn,max);for(i=1;i5;i+)ei0i=ei1i;fi0i=fi1i;for(i=1;ieps&kkm);printf(All do %d timesn,k);sumir=0;sumii=0;for(i=0;i5;i+)I1r=mulRe(ybr0i,-ybi0i,ei0i,-fi0i);I1i=mulIm(ybr0i,-ybi0i,ei0i,-fi0i);sumir+=I1r;sumii+=I1i;pi0=mulRe(ei00,fi00,sumir,sumii);qi0=mulIm(ei00,fi00,sumir,sumii);printf(S1=%f+j%fn,

23、pi0,qi0);ei10=ei00;fi10=fi00;for(i=0;i5;i+)printf(u%d=%f%fn,i+1,sqrt(ei1i*ei1i+fi1i*fi1i),atan(fi1i/ei1i)*180/3.14159);六运行结果：Now start.The 1 timespi1=0.200000pi2=-0.450000pi3=-0.400000pi4=-0.600000pi01=-0.300000pi02=-0.075001pi03=0.000000pi04=0.000000qi01=-0.900000qi02=-0.225000qi03=0.000000qi04=0.0

24、00000detpi1=0.500000detqi1=1.100000detpi2=-0.374999detqi2=0.075000detpi3=-0.400000detqi3=-0.050000detpi4=-0.600000detqi4=-0.100000Iir01=-0.300000Iii01=0.900000Iir02=-0.075001Iii02=0.225000Iir03=0.000000Iii03=0.000000Iir04=0.000000Iii04=0.000000输出雅可比矩阵：33.400010.5340-5.0000-1.6670-5.0000-1.6670-7.500

25、0-2.5000-11.134031.60001.6670-5.00001.6670-5.00002.5000-7.5000-5.0000-1.667038.975012.8420-30.0000-10.00000.00000.00001.6670-5.0000-12.992038.525010.0000-30.00000.00000.0000-5.0000-1.6670-30.0000-10.000038.750012.9170-3.7500-1.25001.6670-5.000010.0000-30.0000-12.917038.75001.2500-3.7500-7.5000-2.500

26、00.00000.0000-3.7500-1.250011.25003.75002.5000-7.50000.00000.00001.2500-3.7500-3.750011.2500detsi0=0.500000detsi1=1.100000detsi2=-0.374999detsi3=0.075000detsi4=-0.400000 detsi5=-0.050000 detsi6=-0.600000 detsi7=-0.100000detui0=-0.047295detui1=0.042961detui2=-0.086292detui3=0.015391detui4=-0.092225de

27、tui5=0.014105detui6=-0.107605detui7=0.009342ui10=-0.047295ui11=1.042961ui12=-0.086292ui13=1.015391ui14=-0.092225ui15=1.014105ui16=-0.107605ui17=1.009342ei11=1.042961fi11=-0.047295ei12=1.015391fi12=-0.086292ei13=1.014105fi13=-0.092225ei14=1.009342fi14=-0.107605max=0.042961Now start.The 2 timespi1=0.2

28、00000pi2=-0.450000pi3=-0.400000pi4=-0.600000pi01=0.277045pi02=-0.449244pi03=-0.410254pi04=-0.616366qi01=0.222040qi02=-0.118305qi03=-0.013816qi04=-0.036371detpi1=-0.077045detqi1=-0.022040detpi2=-0.000756detqi2=-0.031695detpi3=0.010254detqi3=-0.036184detpi4=0.016366detqi4=-0.063629Iir01=0.255454Iii01=

29、-0.224478Iir02=-0.429431Iii02=0.153007Iir03=-0.400000Iii03=0.050001Iir04=-0.600000Iii04=0.100000输出雅可比矩阵：33.159413.0920-5.1360-1.9751-5.1360-1.9751-7.7040-2.9621-12.581133.60831.9751-5.13601.9751-5.13602.9621-7.7040-4.9331-2.124138.384816.0302-29.5988-12.74270.00000.00002.1241-4.9331-16.889038.078812

30、.7427-29.59880.00000.0000-4.9168-2.1516-29.5009-12.907838.155316.2729-3.6876-1.61352.1516-4.916812.9078-29.5009 -17.072938.05531.6135-3.6876-7.3011-3.33040.00000.0000-3.6505-1.665211.05164.39563.3304-7.30110.00000.00001.6652-3.6505-5.595610.8516detsi0=-0.077045 detsi1=-0.022040 detsi2=-0.000756 dets

31、i3=-0.031695detsi4=0.010254detsi5=-0.036184detsi6=0.016366detsi7=-0.063629detui0=-0.000435detui1=-0.007504detui2=0.001742detui3=-0.010067detui4=0.002076detui5=-0.010761detui6=0.003195detui7=-0.013070ui10=-0.047730ui11=1.035457ui12=-0.084550ui13=1.005324ui14=-0.090149ui15=1.003344ui16=-0.104410ui17=0

32、.996272ei11=1.035457fi11=-0.047730ei12=1.005324fi12=-0.084550ei13=1.003344fi13=-0.090149ei14=0.996272fi14=-0.104410max=0.003195Now start.The 3 timespi1=0.200000pi2=-0.450000pi3=-0.400000pi4=-0.600000pi01=0.200525pi02=-0.449917pi03=-0.400025pi04=-0.599997qi01=0.200204qi02=-0.149683qi03=-0.049610qi04=

33、-0.099159detpi1=-0.000525detqi1=-0.000204detpi2=-0.000083detqi2=-0.000317detpi3=0.000025detqi3=-0.000390detpi4=-0.000003detqi4=-0.000841Iir01=0.184354Iii01=-0.201847Iir02=-0.431957Iii02=0.185219Iir03=-0.391092Iii03=0.084584Iir04=-0.585382Iii04=0.160879输出雅可比矩阵：32.933412.9537-5.0977-1.9648-5.0977-1.9648-7.6466-2.9466-12.585033.33711.9648-5.09771.