matlab各种潮流计算代码,基于MATLAB的潮流计算源程序代码

《基于MATLAB的潮流计算源程序代码》由会员分享，可在线阅读，更多相关《基于MATLAB的潮流计算源程序代码(7页珍藏版)》请在人人文库网上搜索。

1、*电力系统直角坐标系下的牛顿拉夫逊法潮流计算*clearclcload E:dataIEEE014_Node.txtNode=IEEE014_Node;weishu=size(Node);nnum=weishu(1,1); %节点总数 load E:dataIEEE014_Branch.txtbranch=IEEE014_Branch;bwei=size(branch); bnum=bwei(1,1); %支路总数Y=(zeros(nnum);Sj=100;%*节点导纳矩阵*for m=1:bnum;s=branch(m,1); %首节点e=branch(m,2); %末节点R=branch(。

2、m,3); %支路电阻X=branch(m,4); %支路电抗B=branch(m,5); %支路对地电纳k=branch(m,6);if k=0 %无变压器支路情形Y(s,e)=-1/(R+j*X); %互导纳Y(e,s)=Y(s,e); endif k=0 %有变压器支路情形Y(s,e)=-(1/(R+j*X)*k);Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/(R+j*X)*k2); Y(e,e)=-(k-1)/(R+j*X)*k); %对地导纳endY(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形endfor t=1:。

3、nnum; Y(t,t)=-sum(Y(t,:)+Node(t,12)+j*Node(t,13); %求支路自导纳endG=real(Y); %电导B=imag(Y); %电纳%*节点分类*pq=0; pv=0; blancenode=0;pqnode=zeros(1,nnum);pvnode=zeros(1,nnum);for m=1:nnum;if Node(m,2)=3blancenode=m; %平衡节点编号else if Node(m,2)=0pq=pq+1;pqnode(1,pq)=m; %PQ节点编号else if Node(m,2)=2pv=pv+1;pvnode(1,pv)=。

4、m; %PV节点编号endendendend%*设置电压初值*Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化for n=1:nnumUoriginal