P-Q 分解法潮流计算程序设计.docVIP

下载本文档

11
0
约1.99千字
约 6页
2017-08-12 发布于安徽
举报
版权申诉

P-Q 分解法潮流计算程序设计.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

z=[0 0.10+0.4*i 0.3*i 0.12+0.5*i;0 0 0 0.08+0.4*i;0 0 0 0;0 0 0 0]; %分别定义所需矩阵的初值，方便编程所需变量% n=length(z); k=zeros(n,n); Y=zeros(n,n); yd=zeros(n,n); y=zeros(n,n); %输入原始数据如下% yd=[0 0.01528*i 0 0.01920*i;0 01528*i 0 0 0.01413*i;0 0 0 0;0.01920*i 0.01413*i 0 0]; %其中yd(i,j)表明i,j结点间与i结点接地电阻% %y(i,j)表明i,j结点间正常联接电阻% k(1,3)=1.1; %数据处理% for m=1:n for j=1:n if z(m,j)~=0 y(m,j)=1/z(m,j); y(j,m)=y(m,j); end end end for m=1:n for j=1:n if k(m,j)~=0 y(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); end end end for m=1:n for j=1:n if m==j Y(m,j)=sum(y(m,:))+sum(yd(m,:)); else Y(m,j)=-y(m,j); Y(j,m)=Y(m,j); end end end A=[-0.3 -0.55 0.5 0;-0.18 -0.13 0 0;1 1 1.1 1.05;0 0 0 0]; G=real(Y); B=imag(Y); %修正方程式的系数矩阵就是导纳矩阵的虚部% B1=B([1,2,3],[1,2,3]); B2=B([1,2,],[1,2,]); for k1=0:100 for m=1:(n-1) sum=0; for j=1:n h=A(m,3)*A(j,3)*(G(m,j)*cos(2*pi/360*(A(m,4)-A(j,4)))+B(m,j)*sin(2*pi/360*(A(m,4)-A(j,4)))); sum=sum+h; end DP(m,1)=A(m,1)-sum; end V1=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; for m=1:2 sum=0; for j=1:n w=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,4)-A(j,4)))); sum=sum+w; end DQ(m,1)=A(m,2)-sum; end V2=A([1,2],[3]); b=DQ./V2;b=inv(-B2)*b; V2=V2+b; A([1,2],[3])=A([1,2],[3])+b; if max(max(abs(DP)),max(abs(DQ)))0.00001 break; end end sum=0; sum1=0; sum2=0; for j=1:n x=A(4,3)*A(j,3)*(G(4,j)*cos(2*pi/360*(A(4,4)-A(j,4)))+B(4,j)*sin(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; end A(4,1)=sum;A(4,2)=sum1;A(3,2)=sum2; for i=1:4 P(i)=A(i,1); Q(i)=A(i,2); V(i)=A(i,3); S(i)=A(i,4); end Y disp(节点1 节点2 节点3 节点4); P disp(节点1 节点2 节点3 节点4); Q disp(节点1 节点2 节点3 节点4); V disp(节点1 节点2 节点3 节点4); S disp(节点1 节点2 节点3 节点4); disp(P Q V S); disp(A);