牛拉潮流程序.docx
《牛拉潮流程序.docx》由会员分享,可在线阅读,更多相关《牛拉潮流程序.docx(9页珍藏版)》请在咨信网上搜索。
1、编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第9页 共9页clear;clc;%+输入带有变压器的支路矩阵中各节点对应各变比% function%=NODE,Branch=OpDF_; %打开矩阵 (test.m)文件Node=NODE;N=Node(:,1); %节点号Type=Node(:,2); %节点类型BR=Branch %将支路信息保存在BR中%K=Branch(:,6); %支路变压器变比,0代表没有变压器n=length(N); %节点数nbr=length(K); %支路数Total_of_Bus1=size(NODE);%取节点矩阵的行和列Total_
2、of_Bus=Total_of_Bus1(1,1)%bus矩阵的行数即 节点数Total_of_Branch1=size(Branch);%取支路矩阵的行和列Total_of_Branch=Total_of_Branch1(1,1);% 支路branch矩阵行数即 支路数Z=zeros(Total_of_Bus1); %将节点排序重新存储节点信息%定义为 节点数的方阵format shortb=1; %排序标志位pq=0; %PQ节点标志位pv=0; %PV节点标志位ph=0; %平衡节点标志位%-按照PQ,PV,平衡节点的次序排序各种节点%-统计PQ节点数 0代表是pq节点for a=1:T
3、otal_of_Bus if NODE(a,2) = 0 Z(b,:)=NODE(a,:); b=b+1; pq=pq+1; endend%-统计PV节点数2代表pv节点for a=1:Total_of_Bus if NODE(a,2) = 2 Z(b,:)=NODE(a,:); b=b+1; pv=pv+1; endend%-统计平衡节点数 3代表平衡节点for a=1:Total_of_Bus if NODE(a,2) = 3 Z(b,:)=NODE(a,:); b=b+1; ph=ph+1; endendZZ2=Z;%将节点进行重新排序%mm=zeros(n,1);for i=1:n m
4、m(i,1)=i;endZ1(:,1)=mm(:,1);Branch1=zeros(nbr,2);for i=1:n if Z(i,1)=Z1(i) for j=1:nbr if Branch(j,1)=Z(i,1) Branch1(j,1)=Z1(i); end if Branch(j,2)=Z(i,1) Branch1(j,2)=Z1(i); end end else for j=1:nbr if Branch(j,1)=Z(i,1) Branch1(j,1)=Z(i,1); end if Branch(j,2)=Z(i,1) Branch1(j,2)=Z(i,1); end end en
5、dendBranch(:,1)=Branch1(:,1);Branch(:,2)=Branch1(:,2);Z(:,1)=Z1(:,1); j=sqrt(-1); %-矩阵已经完成按照PQ,PV,平衡节点的顺序排列起来-YSNODE=Z; %保存排序后的原始节点数据%= Y=zeros(n,n);%求互导纳%for i=1:n for t=1:nbr if (Branch(t,1)=i|Branch(t,2)=i)& Branch(t,6)=0 %非变压器支路% Y(Branch(t,1),Branch(t,2)=-1/(Branch(t,3)+j*Branch(t,4); Y(Branch(
6、t,2),Branch(t,1)=Y(Branch(t,1),Branch(t,2); else if (Branch(t,1)=i|Branch(t,2)=i)&Branch(t,6)=0 %变压器支路% Y(Branch(t,1),Branch(t,2)=(-1/(j*Branch(t,4)/Branch(t,6); Y(Branch(t,2),Branch(t,1)=Y(Branch(t,1),Branch(t,2); end end endend%求自导纳% for i=1:n for t=1:nbr if (Branch(t,1)=i|Branch(t,2)=i)& Branch(t
7、,6)=0 %非变压器支路% Y(i,i)=Y(i,i)+1/(Branch(t,3)+j*Branch(t,4)+(1/2)*j*Branch(t,5); else if Branch(t,1)=i&Branch(t,6)=0 %变压器支路且i为首节点% Y(i,i)=Y(i,i)+1/(j*Branch(t,4); else if Branch(t,2)=i&Branch(t,6)=0%_变压器支路且i为末节点% Y(i,i)=Y(i,i)+(1/(j*Branch(t,4)/(Branch(t,6)*Branch(t,6); end end end endend %若有并联电容器组,则自
8、导纳要加上并联电容器的导纳%for i=1:n if NODE(i,13)=0 Y(i,i)=Y(i,i)+j*NODE(i,13) endendY n=length(N);G=real(Y); %实部,即电导B=imag(Y); %虚部,即电纳%给定初始的电压值与相位值%U_first=Z(:,3); %初始电压幅值phase_first=Z(:,4); %初始相位值e=U_first.*cos(phase_first);f=U_first.*sin(phase_first);%计算Delta_P初始功率量%P=Z(:,5); %节点负荷有功分量Q=Z(:,6); %节点负荷无功分量PG=Z
9、(:,7); %发电机发出的有功QG=Z(:,8); %发电机发出的无功U0=Z(:,9); %节点电压都的初始值Delta_P=zeros(1,n-1);for i=1:n-1 for j=1:n Delta_P(i)=Delta_P(i)-e(i)*(G(i,j)*e(j)-B(i,j)*f(j)+f(j)*(G(i,j)*f(j)+B(i,j)*e(j); endendfor i=1:n-1 Delta_P(i)=Delta_P(i)-(P(i)-PG(i);endDelta_P%计算Delta_Q初始功率量%m=0;for i=1:n; if Type(i)=2; %计算PV节点的个数
10、 m=m+1; %m代表pv节点个数 endendDelta_Q=zeros(1,n-m-1);for i=1:n-m-1 for j=1:n Delta_Q(i)=Delta_Q(i)-f(i)*(G(i,j)*e(j)-B(i,j)*f(j)+e(i)*(G(i,j)*f(j)+B(i,j)*e(j); endendfor i=1:n-m-1 Delta_Q(i)=Delta_Q(i)-(Q(i)-QG(i);endDelta_QDelta_V=zeros(1,m);for i=1:m for j=1:n if Type(j)=2 Delta_V(i)=U0(j)2-(e(j)2+f(i)
11、2); end endendDelta_Vnum=0;disp(第,num2str(num),次时的Delta总的失配量为:) %-进入循环体判断是否满足条件-%-先算出最大值,作为判断是否收敛的依据-%DEL=Delta_P Delta_Q; %Delta_P Delta_Q_%MAX =max(abs(DEL); MAXTheta_first=zeros(1,n);U_f=U_first;Delta_F_E1=Theta_first(1:n-1) U_f(1:n-m-1);Delta_F=Delta_F_E1;Delta_Cor=Delta_F_E1; %_Delta_the Delta_
12、u_%disp(第一次最大失配量误差:,num2str(MAX)%-循环判断-% if MAX1e-004 % 判断依据 disp(-下面开始下一次迭代过程!-)endwhile MAX1e-004 num=num+1;%形成雅克比矩阵%-先求非对角元素-(H)-%Hik=zeros(n-1,n-1);for i=1:n-1 for k=1:n-1 if i=k theik=Theta_first(i)-Theta_first(k); Hik(i,k)=-U_first(i)*U_first(k)*(G(i,k)*sin(theik)-B(i,k)*cos(theik); end endend
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 程序
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【Fis****915】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【Fis****915】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。