PQ分解法计算大电网潮流程序.docxVIP

下载本文档

16
0
约7.32千字
约 14页
2017-05-20 发布于浙江
举报
版权申诉

PQ分解法计算大电网潮流程序.docx

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

PQ分解法计算大电网潮流程序

function PQ %用PQ分解法计算大电网潮流 % %bus数组 1.节点编号 2.节点电压 3.节点电压角度 4.注入有功 5.注入无功 6.节点类型(1PQ 2PV 3平衡) %line数组 1.始端节点编号 2.末端节点编号 3.电阻 4电抗 5电导G 6电纳B 7.变比 %打开数据文件 clear clc bus=load(ieee14bus.txt); line=load(ieee14line.txt); linenum(:,[1,2])=line(:,[1,2]); [nb,~]=size(bus); [nl,~]=size(line); % nodenum=[(1:nb) bus(:,1)]; %带入子函数数据处理 [bus,line,nPQ,nPV,nSW,nodenum] =change1_busline( bus,line );%对节点重新编号 Y = admittance(bus,line,1 );%生成节点导纳矩阵 Y1= admittance(bus,line,2 );%生成化简条件3的矩阵B1 Y2=admittance(bus,line,3 );%生成化简条件3的矩阵B2 %----------------------------------------------------- % %临时添加的测试数据 % nPQ=4; nPV=0;nSW=1;nb=5; % Y=[10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i -5+15i % -1.667+5i 12.917-38.75i -10+30i 0 -1.25+3.75i % -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i 0 % -2.5+7.5i 0 -1.25+3.75i 3.75-11.25i 0 % -5+15i -1.25+3.75i 0 0 6.25-18.75i]; % % Y1=[10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i -5+15i % -1.667+5i 12.917-38.75i -10+30i 0 -1.25+3.75i % -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i 0 % -2.5+7.5i 0 -1.25+3.75i 3.75-11.25i 0 % -5+15i -1.25+3.75i 0 0 6.25-18.75i]; % % bus=[1 1 0 0.2 0.2 1 % 2 1 0 -0.45 -0.15 1 % 3 1 0 -0.4 -0.05 1 % 4 1 0 -0.6 -0.1 1 % 5 1.06 0 0 0 3]; % % line=[5 2 1.25 -3.75 0 0 0 % 2 3 10 -30 0 0 0 % 3 4 1.25 -3.75 0 0 0 % 4 1 2.5 -7.5 0 0 0 % 1 2 1.667 -5 0 0 0 % 1 3 1.667 -5 0 0 0 % 1 5 5 -15 0 0 0]; %------------------------------------------------------- bus_PV0=bus((nPQ+1):end,2);%1.05*ones(1,nPV+nSW); bus_U=[ones(1,nPQ) bus_PV0];%电压幅值 bus_e=zeros(nb,1); %电压角度 delta_P=zeros(nPQ+nPV,1); delta_Q=zeros(nPQ,1); % delta_e=zeros(nb-1,1); %delta_U=zeros(nPQ,1); c=0;KP=1;KQ=1;%KP KQ用来判断有功、无