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

下载本文档

54
0
约6.36千字
约 9页
2016-12-31 发布于重庆
举报
版权申诉

极坐标法求解潮流计算报告.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文档。上传文档

程序清单： %极坐标法求解潮流计算 clear;close;clc; %清变量、图形窗口、命令窗口。 %format long format short Branch=[1 2 0.1+0.4j 0.01528j 1 0 3 1 0.3j 0 1.1 1 1 4 0.12+0.50j 0.01920j 1 0 2 4 0.08+0.40j 0.01413j 1 0 ]; %Branch矩阵：列1、支路首端号；列2、支路末端号；列3、支路阻抗（R）；列4、支路对地导纳（S）； %列5、支路的变化；列6、支路首端处于K侧为1,1侧为0 Y=zeros(4); %定义节点导纳矩阵 for i=1:4 if Branch(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); else j=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)+Branch(i,4); Y(k,k)=Y(k,k)+Branch(i,5)^2/Branch(i,3)+Branch(i,4); end end %生成节点导纳矩阵 disp(节点导纳矩阵:); Y G=real(Y);B=imag(Y); V=[1;1;1.1;1.05];%给定V的初始计算值 U=abs(V); d=angle(V); disp(节点电压的大小：) U disp(节点电压的相位角：) d disp(节点注入有功功率：) Ps=[-0.3;-0.55;0.5;0] disp(节点注入无功功率：) Qs=[-0.18;-0.13;0;0] %由各节点电压向量（状态变量）可得各节点注入功率： for k=1:4 for i=1:4 det(k,i)=d(k)-d(i); end end for i=1:4 temp1=0; temp2=0; for j=1:4 temp1=temp1+U(j)*(G(i,j)*cos(det(i,j))+B(i,j)*sin(det(i,j))); temp2=temp2+U(j)*(G(i,j)*sin(det(i,j))-B(i,j)*cos(det(i,j))); end P(i)=U(i)*temp1; Q(i)=U(i)*temp2; end del_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_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取全