MATLAB极坐标求解潮流计算

1、电力系统分析大作业潮流计算班级：XXX姓名：XXX学号：XXX例题在图所示的简单电力系统中，网络各元件参数的标% %= =o.ioo.io十八.4.4。= =J7uoJ7uo= =yO.O1yO.O1528528q q学= =_/O-3_/O-39 9k k= =1.11.1= =O-O-1 12 2H-H-JO.SQJO.SQ，i4o3410i4o3410= =_ _/ /.0192.0192z z2424= =KO8KO8+ +yO4yO4= =JO.014.3JO.014.3系统中节点 1、 2 为 P。 节点， 节点 3 为叫节点，节点 4 为平衡节点，已给定P Pu uI I鸡，=V

2、.3O=V.3O一用.18.18P P”必=。55-55-刊1313心=0.5=0.5,%,%=1.10=1.10乙=1.05X0*=1.05X0*容许误差二1010一；试用牛顿法。rPQrPQ分解法计算潮流分布程序清单：帆角坐标法求解潮流计算clear;close;clc;Branch=120.1+0.4j0.01528j10310.3jinf1.111 40.12+0.50j0.01920j10240.08+0.40j0.01413j10%Branch矩阵：1、支路首端号；2、支路末端号；3、支路阻抗；4、支路对地导纳；5、支路的变化；6、支路首端处于K侧为1,1侧为0Y=zeros(4)

3、;%节点导纳矩阵fori=1:4ifBranch(i,6)=0%不含变压器的支路j=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-1/Branch(i,3);Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3)+Branch(i,4);Y(k,k)=Y(k,k)+1/Branch(i,3)+Branch(i,4);elsej=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-Branch(i,5)/Branch(i,3);Y(k,j)=Y(j,k);Y(j，j尸Y(j，j)+1/Branch(i,3)

4、;Y(k,k)=Y(k,k)+Branch(i,5)A2/Branch(i,3);endenddisp(节点导纳矩阵,);YG=real(Y);B=imag(Y);V=1;1;1.1;1.05;%给定V的初始计算值e=real(V);f=imag(V);disp(节点电压的实部：,)edisp(节点电压的虚部：,)fdisp(节点注入有功功率：,)Ps=-0.3;-0.55;0.5;0disp(节点注入无功功率：,)Qs=-0.18;-0.13;0;0岫各节点电压向量(状态变量)可得各节点注入功率:P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G

5、*f+B*e);del_W=-0.30-P(1);-0.18-Q(1);-0.55-P(2);-0.13-Q(2);0.5用(3);1.1八2-(3)人2片(3)人2;n=0;%辅助循环计数变量while(any(del_W1e-5)&(n5)n=n+1;disp(迭代次数：);disp(n)%修正方程式:%-0.30-P(1)%-0.18-Q(1)%-0.55-P(2)%-0.13-Q(2)%0.5-P(3)del_e(1)del_f(1)del_e(2)+J*del_f(2)=0del_e(3)%1.1A2-e(3)A2-f(3)A2del_f(3)%雅克比矩阵为：J=P1_e1P

6、1_f1P1_e2P1_f2P1_e3P1_e3;%Q1_e1Q1_f1Q1_e2Q1_f2Q1_e3Q1_f3%P2_e1P2_f1P2_e2P2_f2P2_e3P2_f3%Q2_e1Q2_f1Q2_e2Q2_f2Q2_e3Q2_f3%P3_e1P3_f1P3_e2P3_f2P3_e3P3_f3%V3_e1V3_f1V3_e2V3_f2V3_e3V3_f3%求雅可比矩阵的参数GeBf=G*e-B*f;%辅助计算函数GfBe=G*f+B*e;%辅助计算函数fori=1:2forj=1:3ifi=jP_e(i,j)=-GeBf(i)-G(i,j)*e(i)-B(i,j)*f(i);P_f(i,j

7、)=-GfBe(i)-G(i,j)*f(i)+B(i,j)*e(i);Q_e(i,j)=GfBe(i)+B(i,j)*e(i)-G(i,j)*f(i);Q_f(i,j)=-GeBf(i)+G(i,j)*e(i)+B(i,j)*f(i);elseP_e(i,j)=-G(i,j)*e(i)-B(i,j)*f(i);Q_f(i,j尸-P_e(i,j);P_f(i,j)=B(i,j)*e(i)-G(i,j)*f(i);Q_e(i,j)=P_f(i,j);endendendforj=1:2P_e(3,j)=-G(3,j)*e(3)-B(3,j)*f(3);P_f(3,j)=B(3,j)*e(3)-G(3

8、,j)*f(3);V_e(3,j)=0;V_f(3,j)=0;endP_e(3,3)=-GeBf(3)-G(3,3)*e(3)-B(3,3)*f(3);P_f(3,3)=-GfBe(3)-G(3,3)*f(3)+B(3,3)*e(3);V_e(3,3)=-2*e(3);V_f(3,3)=-2*f(3);%求出雅克比矩阵%J=zeros(6,6);fori=1:3forj=1:3J(2*i-1,2*j-1)=P_e(i,j);J(2*i-1,2*j)=P_f(i,j);endendfori=1:2forj=1:3J(2*i,2*j-1)=Q_e(i,j);J(2*i,2*j)=Q_f(i,j);

9、endendforj=1:3J(6,2*j-1)=V_e(3,j);J(6,2*j)=V_f(3,j);enddisp(,雅可比矩阵为：,)J%解修正方程并得出修正后的电压向量%del_V=-Jdel_W;cleari;fork=1:3V(k)=V(k)+del_V(2*k-1)+i*del_V(2*k);enddisp(节点电压为：,)disp(V)%计算节点不平衡量%e=real(V);f=imag(V);P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G*f+B*e);del_W=-0.30-P(1);-0.18-Q(1);-0.55-P(2

10、);-0.13-Q(2);0.5-P(3);1.1A2-e(3)A2-f(3)A2;disp(节点不平衡量为：)disp(del_W)end%最终Z论%disp(最终各节点电压幅值为：,)disp(abs(V)disp(最终各节点电压相角(度)为：)disp(180*angle(V)/pi)disp(最终各节点注入功率：,)S=P+i*Qdisp(最终各节点注入有功功率为：,)Pdisp(最终各节点注入无功功率为：,)Q运行结果:Branch=1.00002.00000.1000+0.4000i0+0.0153i1.000003.00001.00000+0.3000iInf1.10001.00

11、001.00004.00000.1200+0.5000i0+0.0192i1.000002.00004.00000.0800+0.4000i0+0.0141i1.00000节点导纳矩阵Y=1.0421-8.2429i-0.5882+2.3529i0+3.6667i-0.4539+1.8911i-0.5882+2.3529i1.0690-4.7274i0-0.4808+2.4038i0+3.6667i00-3.3333i0-0.4539+1.8911i-0.4808+2.4038i0节点电压的实部：e=1.00001.00001.10001.0500节点电压的虚部：f=0000节点注入有功功率:

12、Ps=-0.3000-0.55000.500000.9346-4.2616i节点注入无功功率:Qs=-0.1800-0.130000迭代次数：1雅可比矩阵为：-1.0194-8.37190.58822.352903.6667-8.11381.06482.3529-0.58823.666700.58822.3529-1.0450-4.8770002.3529-0.5882-4.57781.09300004.0333000-3.66670000-2.20000节点电压为：0.9935-0.0088i0.9763-0.1078i1.1000+0.1267i1.0500节点不平衡量为：-0.0013-

13、0.0028-0.0135-0.05470.0030-0.0160迭代次数：2雅可比矩阵为：-0.8091-8.36130.60522.33250.03243.6429-7.99921.40712.3325-0.60523.6429-0.03240.82802.2338-1.0190-4.6364002.2338-0.8280-4.36412.087800-0.46444.033300-0.0324-3.64290000-2.2000-0.2533节点电压为：0.9847-0.0086i0.9590-0.1084i1.0924+0.1289i1.0500节点不平衡量为：-0.0000-0.0000-0.0003-0.00110.