节点导矩阵的建立.docVIP

下载本文档

39
0
约 4页
2017-06-14 发布于河南
举报
版权申诉

节点导矩阵的建立.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文档。上传文档

节点导矩阵的建立

如上图所示的简单电力系统中，网络各元件参数的标幺值如下： z12=0.10+j0.40 y120=y210=j0.01528 z13=j0.3,k=1.1 z14=0.12+j0.50 y140=y410=j0.01920 z24=0.08+j0.40 y240=y420=j0.01413 系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点。节点导纳矩阵的运行程序如下： clc Clear disp(网络各元件参数用标幺值表示！！！); N0=input(请输入节点数:N0=); n1=input(请输入支路数:n1=); l=input(请输入PQ节点的个数=); for m=1:l c(m)=input([请输入第,num2str(m),个PQ节点的节点号为：]); end t=input(请输入PV节点的个数=); for m=1:t ??? c(m)=input([请输入第,num2str(m),个PV节点的节点号为：]); end b=input(请输入平衡节点号:b=); %%由支路参数形成矩阵B1 disp(各支路连接情况：) ?i=1; for m=1:n1 ????? syms Y N ??? p=input([第,num2str(m),条支路的起始节点]); ??? q=input([第,num2str(m),条支路的终止节点]); ??? mn=input([第,num2str(m),条支路是否有变压器（请输入‘Y’或‘N’）]); ??? y=0;k=1; ??? if mn==Y; ??????? k=input(请输入变压器变比(标幺值):); ??????? z=input([请输入第,num2str(m),条支路的线路阻抗]); ??? else ???????? z=input([请输入第,num2str(m),条支路的线路阻抗：]); ???????? y=input([请输入第,num2str(m),条支路线路的对地阻抗：]); ??? end ?? ??? B1(i,1)=p;B1(i,2)=q;B1(i,3)=z;B1(i,4)=y;B1(i,5)=1/k; ??? i=i+1; end disp(由支路参数形成的矩阵B1) B1 %求节点导纳矩阵 Y=zeros(N0); e=zeros(1,N0); f=zeros(1,N0); for i=1:n1??????? ??? ?????? p=B1(i,1);q=B1(i,2); ??? Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); ??? Y(q,p)=Y(p,q);????????????????????? ??? Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;? ??? Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;??? End disp(导纳矩阵 Y=); disp(Y) 程序运行结果：