matlab 应求pq分解法进行潮流计算（Matlab should be solved by PQ decomposition method for power flow calculation）.docVIP

下载本文档

11
0
约5.06千字
约 19页
2017-09-05 发布于河南
举报
版权申诉

matlab 应求pq分解法进行潮流计算（Matlab should be solved by PQ decomposition method for power flow calculation）.doc

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

matlab 应求pq分解法进行潮流计算（Matlab should be solved by PQ decomposition method for power flow calculation）

matlab 应求pq分解法进行潮流计算（Matlab should be solved by PQ decomposition method for power flow calculation） % of this procedure is the function decomposition method for power flow calculation using PQ Clear CLC N=input (please input node number: n=); Nl=input (please input branch number: nl=); Isb=input (please input balance bus number: isb=); Pr=input (please input error: pr=); Disp (please input by the branch parameter matrix formed: B1=); Input B1% [fname pname]= uigetfile (*.xls); File=[pname fname]; X=xlsread (file); B1=[x (:, 1), X (:, 2), (x (: 3) (+x: *i, 4), X (5): X (*i), X (: 6),:, 7)]; Disp (please input by the node parameter matrix formed: B2=); Input B2% [fname pname]= uigetfile (*.xls); File=[pname fname]; Y=xlsread (file); B2=[(Y (: 1) (+y: 2) (*i), y (: 3) (+y: 4), y (*i): 5), y (:, 6), y (:, 7), y (:, 8)]; Disp (please input by the node number and the impedance matrix formed: X=); Input X% [fname pname]= uigetfile (*.xls); File=[pname fname]; Z=xlsread (file); X=z; Na=input (please enter the PQ node number na=); Y=zeros (n); YI=zeros (n); e=zeros (1, n); f=zeros (1, n); V=zeros (1, n); O=zeros (1, n); For i=1:n If X (I, 2) ~=0; P=X (I, 1); Y (P, P) =1./X (I, 2); End End For i=1:nl If B1 (I, 6) ==0 P=B1 (I, 1); q=B1 (I, 2); Else p=B1 (I, 2); q=B1 (I, 1); End Y (P, q) =Y (P, q) -1./ (B1 (I, 3) *B1 (I, 5)); YI (P, q) =YI (P, q) -1./B1 (I, 3); Y (Q, P) =Y (P, Q); YI (Q, P) =YI (P, Q); Y (Q, q) =Y (Q, q) +1./ (B1 (I, 3) *B1 (I, 5) ^2) +B1 (I, 4)./2; YI (Q, q) =YI (Q, q) +1./B1 (I, 3); Y (P, P) =Y (P, P) +1./B1 (I, 3) +B1 (I, 4)./2; YI (P, P) =YI (P, P) +1./B1 (I, 3); End% admittance matrix G=real (Y); B=imag (YI); BI=imag (Y); For i=1:n S (I) =B2 (I, 1) -B2 (I, 2); BI (I, I) =BI (I, I) +B2 (I, 5); End P=real (S); Q=imag (S); For i=1:n E (I) =real (B2 (I, 3)); F (I) =imag (B2 (I, 3)); V (I) =B2 (I, 4); End For i=1:n If B2 (I, 6) ==2 V (I) =sqrt (E (I) ^2+f (I) ^2); O (I) =atan (f (I)./e (I)); End End For i=2:n If i==n B (I, I) =1./B (I, I); Else IC1=i+1; For j1=