资源描述：
《牛拉潮流程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、clear;clc;%%%+++++++++++++++++++++++++++++++++++++输入带有变压器的支路矩阵中各节点对应各变比%function[]%==========================================================================[NODE,Branch]=OpDF_;%打开矩阵(test.m)文件Node=NODE;N=Node(:,1);%节点号Type=Node(:,2);%节点类型BR=Branch%将支路信息保存在BR中%%K=Branch(:,6);%支路变压器变比，0代表没有变压器n=len
2、gth(N);%节点数nbr=length(K);%支路数37/37Total_of_Bus1=size(NODE);%%%取节点矩阵的行和列Total_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);%将节点排序重新存储节点信息%%%%定义为节点数的方阵formatshortb=1;%排序标志位pq=0;%PQ节点标
3、志位pv=0;%PV节点标志位ph=0;%平衡节点标志位37/37%---------------------按照PQ,PV,平衡节点的次序排序各种节点%-------------------------------------统计PQ节点数0代表是pq节点fora=1:Total_of_BusifNODE(a,2)==0Z(b,:)=NODE(a,:);b=b+1;pq=pq+1;endend%-------------------------------------统计PV节点数2代表pv节点fora=1:Total_of_BusifNODE(a,2)==237/37Z(b,:)=N
4、ODE(a,:);b=b+1;pv=pv+1;endend%-------------------------------------统计平衡节点数3代表平衡节点fora=1:Total_of_BusifNODE(a,2)==3Z(b,:)=NODE(a,:);b=b+1;ph=ph+1;endend37/37ZZ2=Z;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%将节点进行重新排序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%mm=zeros(n,1);fori=1:nmm(i,1)=i;endZ1(:,1)=mm(:,1);Branch1=ze
5、ros(nbr,2);fori=1:nifZ(i,1)~=Z1(i)forj=1:nbrifBranch(j,1)==Z(i,1)37/37Branch1(j,1)=Z1(i);endifBranch(j,2)==Z(i,1)Branch1(j,2)=Z1(i);endendelseforj=1:nbrifBranch(j,1)==Z(i,1)Branch1(j,1)=Z(i,1);endifBranch(j,2)==Z(i,1)Branch1(j,2)=Z(i,1);37/37endendendendBranch(:,1)=Branch1(:,1);Branch(:,2)=Branch
6、1(:,2);Z(:,1)=Z1(:,1);j=sqrt(-1);%--------矩阵已经完成按照PQ，PV，平衡节点的顺序排列起来------YSNODE=Z;%保存排序后的原始节点数据%%%%=======================================================================37/37Y=zeros(n,n);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%求互导纳%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fori=1:nfort=1:nbrif(Branch(t,1
7、)==i
8、
9、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(t,2),Branch(t,1))=Y(Branch(t,1),Branch(t,2));elseif37/37(Branch(t,1)==i
10、
11、Branch(t,2)==i)&&Branch(t,
