潮流计算实验

1、电力系统分析实验报告 实验一：潮流计算的计算机算法 clear; n=10; nl=10; isb=1; pr=0.00001; B 1 =1 2 0. 0 35 1 2 + 0. 0 83 0 6i 0. 1 3455 i 1 0; 2 3 0. 0068 + 0. 1 8375i 0 1. 0 2 38 1 1 ; 1 4 0. 0 56 2 0+ 0. 1 3 2 89 i 0 . 05 382 i 1 0; 4 5 0. 00 8 1 1 + 0. 2 4 54 9i 0 1. 02 38 1 1; 1 6 0. 0 56 2 0+ 0. 1 3 2 89 i 0 . 05 382 i

2、 1 0; 4 6 0. 042 1 5+ 0. 09 9 67i 0 . 04 037 i 1 0; 6 7 0. 0 0 68 + 0. 1 8 3 75i 0 1. 0 2 38 1 1; 6 8 0. 0 28 1 0+ 0. 06 645i 0.1 0 764 i 1 0; 8 1 0 0. 00 8 1 1 + 0. 2454 9i 0 1 1; 8 9 0. 0 35 1 2 + 0. 0 8 30 6i 0. 1 345 5i 1 0 B2 =0 0 1. 1 1. 1 0 1; 0 0 1 0 0 2; 0 0. 343 + 0. 21256i1 0 0 2; 0 0 1 0

3、 0 2; 0 0. 204 + 0. 1 2638i1 0 0 2; 0 0 1 0 0 2; 0 0. 306 + 0. 1 8962 i1 0 0 2; 0 0 1 0 0 2; 0. 5 0 1. 1 1. 1 0 3; 0 0. 343 + 0. 2 1 256i1 0 0 2 Y=zeros (n );e=zeros(1, n);f=zeros(1, n);V=zeros(1, n);sida=zeros(1, n);S仁zeros( nl); for i=1: nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1

4、); end Y(p ,q)=Y( p,q)-1./(B1(i,3)*B1(i,5); Y(q, p)=Y( p,q)； Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)A2)+B1(i,4)./2; Y( p,p )=Y( p,p )+1./B1(i,3)+B1(i,4)./2; end dis PC导纳矩阵丫=); dis p(Y) % G=real(Y);B=imag(Y); for i=1: n e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4); end for i=1: n S(i)=B2(i,1)-B2(i,2)

5、; B(i,i)=B(i,i)+B2(i,5); end P=real(S);Q=imag(S); ICT1=0;IT2=1;N0=2* n;N=N0+1;a=0; while IT2=0 IT2=0;a=a+1; for i=1: n if i=isb C(i)=0;D(i)=0; for j1=1: n 工(Gij*ej-Bij*fj) 2 (Gij*fj+Bij*ej) C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);% end P 仁C(i)*e(i)+f(i)*D(i); Q1

6、=C(i)*f(i)-e(i)*D(i); V2=e(i)A2+f(i)A2; if B2(i,6)=3 DP=P (i)-P1; DQ=Q(i)-Q1; for j1=1: n if j1=isb X2=B(i,j1)*e(i)-G(i,j1)*f(i); X3=X2; X4=-X1; p=2*i-1;q=2*j1-1; J( p,q)=X3;J( p, N)=DQ;m=p+1; J(m,q)=X1;J(m,N)=D P;q=q+1; J(p ,q)=X4;J(m,q)=X2; elseif j1=i X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i); X3=D(i)+B(i

7、,i)*e(i)-G(i,i)*f(i); X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i); p=2*i-1;q=2*j1-1;J( p,q)=X3;J( p,N)=DQ; m=p+1; J(m,q)=X1;q=q+1;J( p,q)=X4;J(m,N)=D P; J(m,q)=X2; end end else DP=P (i)-P1; DV=V(i)A2-V2; for j1=1: n if j1=isb X2=B(i,j1)*e(i)-G(i,j1)*f(i); X5=0;X6=0; p=2*i-1;q=2*j1-1;J( p,q)=X5;J( p,N)=DV; m=p+

8、1; J(m,q)=X1;J(m,N)=D P;q=q+1;J( p,q)=X6; J(m,q)=X2; elseif j1=i X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i); X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J( p,q)=X5;J( p,N)=DV; m=p+1; J(m,q)=X1;J(m,N)=D P;q=q+1;J( p,q)=X6; J(m,q)=X2; end end end end end for k=3:N0 k仁k+1;N1=N; for k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k);

9、 end J(k,k)=1; if k=3 k4=k-1; for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end if k=N0 break; end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2); end J(k3,k)=0; end end end

10、 for k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); k1=k+1; f(L)=f(L)-J(k1,N); end for k=3:N0 DET=abs(J(k,N); if DET=pr IT2=IT2+1; end end ICT2(a)=IT2; ICT1=ICT1+1; end disp(迭代次数：); disp( ICT1); disp(没有达到精度要求的个数：); dis p( ICT2); for k=1: n V(k)=sqrt(e(k)A2+f(k)A2); sida(k)=ata n(f(k)./e(k)*180./pi; E(k)=e

11、(k)+f(k)*j; end disp(各节点的实际电压标幺值E为(节点号从小到大排列):); disp( E); dispC); disp(各节点的电压大小V为(节点号从小到大排列)：)； disp (V); dispC); disp(各节点的电压相角sida为(节点号从小到大排列)：)； dis p( sida); for p=1: n C(p )=0； for q=1: n C(p )=C( p)+conj(Y( p,q)*conj(E(q); end S(p )=E( p)*C( p); end disp(各节点的功率S为(节点号从小到大排列)：); disp (S); dispC)

12、; disp(各条支路的首端功率Si为(顺序同输入B1时一致)：)； for i=1: nl p=B1(i,1);q=B1(i,2); if B1(i,6)=0 Si( p, q)=E( p)*(conj(E( p)*conj(B1(i,4)./2)+(conj(E( p)*B1(i,5). -conj(E(q)*conj(1./(B1(i,3)*B1(i,5); Siz(i)=Si (p ,q); else Si( p, q)=E( p)*(conj(E( p)广conj(B1(i,4)./2)+(conj(E( p) ./B1(i,5). -conj(E(q)*conj(1./(B1(i,

13、3)*B1(i,5); Siz(i)=Si (p ,q); end dis p(Si( p,q); SSi( p, q)=Si (p ,q); ZF=S(, nu m2str( p), nu m2str(q),)= ,n um2str(SSi( p,q); dis p( ZF); disp(); end disp(各条支路的末端功率Sj为(顺序同输入B1时一致)：)； for i=1: nl p=B1(i,1);q=B1(i,2); if B1(i,6)=0 Sj(q, p)=E(q)*(conj(E(q)*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5). -conj

14、(E( p)*conj(1./(B1(i,3)*B1(i,5); Sjy(i)=Sj(q, p); else Sj(q, p)=E(q)*(conj(E(q)*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5). -conj(E( p)*conj(1./(B1(i,3)*B1(i,5); Sjy(i)=Sj(q, p); end dis p(Sj(q, p); SSj(q, p)=Sj(q, p); ZF=S(, nu m2str(q), ,n um2str( p),)=, nu m2str(SSj(q, p); dis p( ZF); disp(); end disp(各

15、条支路的功率损耗DS为(顺序同输入B1时一致)：); for i=1: nl p=B1(i,1);q=B1(i,2); DS(i)=Si( p, q)+Sj(q, p); dis p(DS(i); DDS(i)=DS(i); ZF=DS(, nu m2str( p), nu m2str(q),)= ,n um2str(DDS(i); dis p( ZF); disp(); end figure(1); sub plot(1,2,1); plot(V); xIabelC节点号);ylabel(电压标幺值); grid on; sub plot(1,2,2); pl ot(sida); xIabe

16、lC节点号);ylabel(电压角度); grid on; figure(2); sub plot(2,2,1); P=real(S);Q=imag(S); bar( P); xlabel(节点号);ylabel(节点注入有功); grid on; sub plot(2,2,2); bar(Q); xlabel(节点号);ylabel(节点注入无功); grid on; sub plot(2,2,3); P1=real(Siz);Q1=imag(Siz); bar(P 1); ); ); xlabel(支路号);ylabel(支路首端注入有功 grid on; sub plot(2,2,4);

17、 bar(Q1); xlabel(支路号);ylabel(支路首端注入无功 grid on; B1 = 1.0000 2.0000 0.0351 + 0.0831i 0 + 0.1346i 1.0000 ( 2.0000 3.0000 0.0068 + 0.1838i 0 1.0238 1.0000 1.0000 4.0000 0.0562 + 0.1329i 0 + 0.05 38i 1.0000 ( 4.0000 5.0000 0.0081 + 0.2455i 0 1.0238 1.0000 1.0000 6.0000 0.0562 + 0.1329i 0 + 0.0538i 1.0000

18、 ( 4.0000 6.0000 0.0422 + 0.0997i 0 + 0.0404i 1.0000 1 6.0000 7.0000 0.0068 + 0.1838i 0 1.0238 1.0000 6.0000 8.0000 0.0281 + 0.0664i 0 + 0.1076i 1.0000 ( 8.0000 10.0000 0.0081 + 0.2455i 0 1.0000 1.0000 8.0000 9.0000 0.0351 + 0.0831i 0 + 0.1346i 1.0000 ( 0 0 0 0 0 0 B2 = 0 0 1.1000 1.1000 0 1.0000 0

19、0 1.0000 0 0 2.0000 0 0.3430 + 0.2126i 1.0000 0 0 2.0000 0 0 1.0000 0 0 2.0000 0 0.2040 + 0.1264i 1.0000 0 0 2.0000 0 0 1.0000 0 0 2.0000 0.5000 导纳矩阵 0.3060+ 0.1896i 0.3430+ 0.2126i Y= Colu mns 1 through 6 9.7177 -22.8591i -2.6996 + 6.3834i -4.3185 +10.2135i -2.6996 + 6.3834i 1.0000 1.0000 1.1000 1.

20、0000 -4.3185 +10.2135i 4.5104 -15.3311i -0.1964 -0.1964 5.3083i 0.2011 2.0000 2.0000 1.1000 3.0000 2.0000 5.3083i -5.4347i -2.6996 + 6.3834i 6.4271 -18.7292i -0.1313 + 3.9744i -3.5993 + 8.5110i -0.1313 + 3.9744i0.1344 - 4.0690i -2.6996 + 6.3834i -3.5993 + 8.5110i 11.8891 -32.7444i -0.1964 + 5.3083i

21、-5.3984 +12.7660i Colu mns 7 through 10 -0.1964 + 5.3083i -5.3984 +12.7660i 4 0.2011 - 5.4347i 9.8514 -26.9275i -4.3185 +10.2135i -0.1344 + 4.0690i -4.3185 +10.2135i 4.3185 -10.1462i -0.1344 + 4.0690i 0.1344 - 4.0690i 迭代次数: 没有达到精度要求的个数: 18 各节点的实际电压标幺值 E为（节点号从小到大排列）: Colu mns 1 through 6 1.1000 1.075

22、7 - 0.0207i1.0050 - 0.0780i1.0772 - 0.0175i 1.0171 - 0.0631i1.0762 - 0.0152i Colu mns 7 through 10 1.0112 - 0.0666i1.0778 - 0.0051i1.0996 + 0.0304i 1.0177 - 0.0814i 各节点的电压大小 V为（节点号从小到大排列）： 1.1000 1.0759 1.0080 1.0773 1.0191 1.0763 1.0134 1.0778 1.1000 1.0209 各节点的电压相角 sida为（节点号从小到大排列 ）： 0-1.1046-4.43

23、73-0.9283 -3.5503 -0.8106 -3.7665 -0.2718 S(2,3)=0.34409+0.24201i 1.5822-4.5707 各节点的功率S为（节点号从小到大排列）: Colu mns 1 through 6 0.7165 + 0.2587i0.0000 + 0.0000i-0.3430 - 0.2126i -0.0000 0.0000i -0.2040 0.1264i-0.0000 + 0.0000i Colu mns 7 through 10 -0.3060 - 0.1896i-0.0000 + 0.0000i0.5000 + 0.0089i-0.3430

24、 - 0.2126i 各条支路的首端功率Si为（顺序同输入B1时一致）： 0.3485 + 0.0932i S(1,2)=0.3485+0.093157i 0.3441 + 0.2420i 0.1904 + 0.0760i S(1,4)=0.19038+0.07599i 0.2044 + 0.1400i S(4,5)=0.20445+0.13999i 0.1777 + 0.0895i S(1,6)=0.17767+0.089525i -0.0163 - 0.0055i S(4,6)=-0.016305-0.0054856i 0.3069 + 0.2128i S(6,7)=0.30686+0.2

25、1281i -0.1477 - 0.0234i S(6,8)=-0.14767-0.02338i 0.3443 + 0.2509i S(8,10)=0.34427+0.25091i -0.4925 - 0.1508i S(8,9)=-0.49251-0.15077i 各条支路的末端功率 Sj为（顺序同输入B1时一致）: -0.3441 - 0.2420i S(2,1)=-0.34409-0.24201i -0.3430 - 0.2126i S(3,2)=-0.343-0.21256i -0.1881 - 0.1345i S(4,1)=-0.18815-0.13451i -0.2040 - 0.1264i S(5,4)=-0.204-0.12638i -0.1755 - 0.1482i S(6,1)=-0.17551-0.14815i 0.0163 - 0.0413i S(6,4)=0.0163