p-q 分解法潮流计算程序设计

《p-q 分解法潮流计算程序设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、z=[00.10+0.4*i0.3*i0.12+0.5*i;0000.08+0.4*i;0000;0000];%分别定义所需矩阵的初值，方便编程所需变量%n=length(z);k=zeros(n,n);Y=zeros(n,n);yd=zeros(n,n);y=zeros(n,n);%输入原始数据如下%yd=[00.01528*i00.01920*i;001528*i000.01413*i;0000;0.01920*i0.01413*i00];%其中yd(i,j)表明i,j结点间与i结点接地电阻%%y(i,j)表明i,j结点间

2、正常联接电阻%k(1,3)=1.1;%数据处理%form=1:nforj=1:nifz(m,j)~=0y(m,j)=1/z(m,j);y(j,m)=y(m,j);endendendform=1:nforj=1:nifk(m,j)~=0y(m,j)=k(m,j)/z(m,j);y(j,m)=y(m,j);yd(m,j)=(k(m,j)-1)*k(m,j)/z(m,j);yd(j,m)=(1-k(m,j))/z(m,j);endendendform=1:nforj=1:nifm==jY(m,j)=sum(y(m,:))+sum(y

3、d(m,:));elseY(m,j)=-y(m,j);Y(j,m)=Y(m,j);endendendA=[-0.3-0.550.50;-0.18-0.1300;111.11.05;0000];G=real(Y);B=imag(Y);%修正方程式的系数矩阵就是导纳矩阵的虚部%B1=B([1,2,3],[1,2,3]);B2=B([1,2,],[1,2,]);fork1=0:100form=1:(n-1)sum=0;forj=1:nh=A(m,3)*A(j,3)*(G(m,j)*cos(2*pi/360*(A(m,4)-A(j,4

4、)))+B(m,j)*sin(2*pi/360*(A(m,4)-A(j,4))));sum=sum+h;endDP(m,1)=A(m,1)-sum;endV1=A([1,2,3],[3]);a=DP./V1;a=inv(-B1)*180/pi*a;DS=V1.a;A([1,2,3],[4])=A([1,2,3],[4])+DS;form=1:2sum=0;forj=1:nw=A(m,3)*A(j,3)*(G(m,j)*sin(2*pi/360*(A(m,4)-A(j,4)))-B(m,j)*cos(2*pi/360*(A(m

5、,4)-A(j,4))));sum=sum+w;endDQ(m,1)=A(m,2)-sum;endV2=A([1,2],[3]);b=DQ./V2;b=inv(-B2)*b;V2=V2+b;A([1,2],[3])=A([1,2],[3])+b;ifmax(max(abs(DP)),max(abs(DQ)))<0.00001break;endendsum=0;sum1=0;sum2=0;forj=1:nx=A(4,3)*A(j,3)*(G(4,j)*cos(2*pi/360*(A(4,4)-A(j,4)))+B(4,j)*si

6、n(2*pi/360*(A(4,4)-A(j,4))));sum=sum+x;c=A(4,3)*A(j,3)*(G(4,j)*sin(2*pi/360*(A(4,4)-A(j,4)))-B(4,j)*cos(2*pi/360*(A(4,4)-A(j,4))));sum1=sum1+c;d=A(3,3)*A(j,3)*(G(3,j)*sin(2*pi/360*(A(3,4)-A(j,4)))-B(3,j)*cos(2*pi/360*(A(3,4)-A(j,4))));sum2=sum2+d;endA(4,1)=sum;A(4,2

7、)=sum1;A(3,2)=sum2;fori=1:4P(i)=A(i,1);Q(i)=A(i,2);V(i)=A(i,3);S(i)=A(i,4);endYdisp('节点1节点2节点3节点4');Pdisp('节点1节点2节点3节点4');Qdisp('节点1节点2节点3节点4');Vdisp('节点1节点2节点3节点4');Sdisp('节点1节点2节点3节点4');disp('PQVS');disp(A);