极坐标法求解潮流计算报告

1、例题么值如下:在图所示的简单电力系统中，网络各元件参数的标 石g = O.1LO -I- yO.4O -V120 = jO.O528例题么值如下:= Mg = 0-01413系统中节点1、2为PQ节点，节点3为PV 节点，节点4为平衡节点，已给定马十/&一30-= Mg = 0-01413系统中节点1、2为PQ节点，节点3为PV 节点，节点4为平衡节点，已给定马十/&一30-处.18A, I ；e2j= 0.55-/0.13乌=。5,吃=11。=105Z0容许误差 = 10一1试用牛顿法orPQ分解法计算潮流分布。程序清单：%极坐标法求解潮流计算clear;close;clc; %清变量、图形

2、窗口、命令窗口。%format longformat shortBranch=1 2 0.1+0.4j0.01528j 1 03 10.3j01.1 14 0.12+0.50j 0.01920j 1 04 0.08+0.40j 0.01413j 1 0;%Branch矩阵：列1、支路首端号；列2、支路末端号；列3、支路阻抗(R)；列4、支路对 地导纳(S)；%列5、支路的变化；列6、支路首端处于K侧为1,1侧为0Y=zeros(4); %定义节点导纳矩阵for i=1:4if Branch(i,6)=0%不含变压器的支路j=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(

3、j,k)-1/Branch(i,3);Y(kj)=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(kj)=Y(j,k);Y(jj)=Y(jj)+1/Branch(i,3)+Branch(i,4);Y(k,k)=Y(k,k)+Branch(i,5)A2/Branch(i,3)+Branch(i,4);endend%生成节点导纳矩阵

4、disp(节点导纳矩阵:)；YG=real(Y);B=imag(Y);V=1;1;1.1;1.05;%给定V的初始计算值U=abs(V);d=angle(V);disp(*节点电压的大小：)Udisp(*节点电压的相位角：)ddisp(*节点注入有功功率：)Ps=-0.3;-0.55;0.5;0disp(*节点注入无功功率：)Qs=-0.18;-0.13;0;0%由各节点电压向量(状态变量)可得各节点注入功率：for k=1:4for i=1:4det(k,i)=d(k)-d(i);endendfor i=1:4temp1=0;temp2=0;for j=1:4temp1=temp1+U(j)

5、*(G(ij)*cos(det(ij)+B(ij)*sin(det(ij);temp2=temp2+U(j)*(G(iJ)*sin(det(iJ)-B(iJ)*cos(det(ij);endP(i)=U(i)*temp1;Q(i)=U(i)*temp2;enddel_W= -0.30-P(1);-0.55-P(2);0.50-P(3);-0.18-Q(1);-0.13-Q(2)n=0;%辅助循环计数变量while (any(abs(del_W)1e-5)&(n15)n=n+1;disp(迭代次数：)；disp(n)%修正方程式：%Ps(1)-P(1)del_d(1)%Ps(2)-P(2)del

6、_d(2)%Ps(3)-P(3)=-J1 * del_d(3)%Qs(2)-Q (2)del_U(2)/U(2)%Qs(1)-Q(1)del_U(1)/U(1)%雅克比矩阵为：j= H N%M L %其中：H取全矩阵；N取全行，最后PV的节点数列略去；M取全列， 最后PV个节点数行略去；%L少取PV个节点数行和PV个节点数 列！很重要% %求雅可比矩阵的参数%H=zeros(3);N=zeros(3);M=zeros(3);L=zeros(3);for i=1:3for j=1:3 if i=jH(i,j)= U(i*B(i,i)+Q(i);N(ij)=-U(i)人2*G(i,i)-P(i);

7、M(i,j)= U(i)人2*G(i,i)-P(i);L(i,j)= U(i)A2*B(i,i)-Q(i);endif i=jH(iJ)=-U(i)*U(j)*(G(ij)*sin(det(iJ)-B(ij)*cos(det(ij);N(iJ)=-U(i)*U(j)*(G(iJ)*cos(det(iJ)+B(iJ)*sin(det(iJ);M(ij)=-N(i,j);L(iJ)=H(iJ);endendend%-求出雅克比矩阵j%J=zeros(5);J= H N(:,1:2) M(1:2,:) L(1:2,1:2);disp(雅可比矩阵为：)J%解修正方程并得出修正后的电压向量%del_du

8、=-Jdel_W for k=1:3 d(k)=d(k)+del_du(k);endfor k=1:2U(k)=U(k)+del_du(k+3)*U(k);end clear jdisp(节点电压为：) disp(U.*exp(j*d) %计算节点不平衡量%for k=1:4for i=1:4det(k,i)=d(k)-d(i);endendfor i=1:4temp1=0;temp2=0; for j=1:4temp1=temp1+U(j)*(G(ij)*cos(det(ij)+B(ij)*sin(det(ij);temp2=temp2+U)*(G(i,j)*sin(det(i,j)-B(i

9、,j)*cos(det(i,j); endP(i)=U(i)*temp1;Q(i)=U(i)*temp2;enddel_W= -0.30-P(1);-055-P(2);0.50-P(3);-0.18-Q(1);-0.13-Q(2);disp(节点不平衡量为：)disp(del_W)end%-最终结论%clear jdisp(最终各节点电压幅值为：)disp(abs(U.*exp(j*d)disp(最终各节点电压相角(度)为：)disp(d)disp(最终各节点注入功率：)S=P+j*Qdisp(最终各节点注入有功功率为：)Pdisp(最终各节点注入无功功率为：)Q运行结果：节点导纳矩阵：Y =

10、0 + 3.6667i -0.4539 + 1.8911i0-0.4808 + 2.4038i0 - 3.3333i0 + 3.6667i -0.4539 + 1.8911i0-0.4808 + 2.4038i0 - 3.3333i000.9346 - 4.2616i0 + 3.6667i0-0.4539 + 1.8911i -0.4808 + 2.4038i节点电压的大小：U =1.00001.00001.10001.0500节点电压的相位角：d =0000节点注入有功功率:Ps =-0.3000-0.55000.5000 0节点注入无功功率:Qs =-0.1800-0.1300 0 0de

11、l_W =-0.2773-0.52600.5000-0.05100.0196迭代次数：1雅可比矩阵为:-8.37192.35294.0333-1.01940.58822.3529-4.877000.5882-1.04504.03330-4.033301.0648-0.58820-8.11382.3529-0.58821.093002.3529-4.5778del_du =-0.0088-0.10780.1151-0.0065-0.0237节点电压为：0.9935 - 0.0088i 0.9707 - 0.1051i 1.0927 + 0.1264i1.0500节点不平衡量为：-0.0000-0

12、.02050.0045-0.0438 -0.0245迭代次数：2雅可比矩阵为：-8.27252.2148-8.27252.21483.97641.3286-0.79342.3276-4.6118 0-0.34221.54853.9764-0.72870.342200.7934-0.4895-3.9764-0.49550-0.4955-8.00012.327602.2148-4.4009del_du =0.0001-0.00470.0023-0.0088-0.0117节点电压为：0.9847 - 0.0086i 0.9589 - 0.1083i 1.0924 + 0.1289i 1.0500节点

13、不平衡量为:1.0e-003 *0.0995-0.41950.0793-0.4503-0.3181迭代次数：3雅可比矩阵为：-8.17332.28183.9403-0.71050.32442.1660-4.531500.7875-0.44583.94030-3.9403-0.499901.3107-0.3244-0.4999-7.81422.2818-0.78751.545002.1660-4.2721del_du =1.0e-003 *0.0109-0.08920.0440-0.1017 -0.1603节点电压为：0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1290i1.0500节点不平衡量为： 1.0e-006 *0.0417-0.11040.0135-0.0666-0.0666最终各节点电压幅值为：0.98470.96481.1