牛拉潮流程序.docxVIP

clc; %%%+++++++++++++++++++++++++++++++++++输++入带有变压器的支路矩阵中各节点对应各变比 % function[] [NODE,Branch]=OpDF_; %打开矩阵 文件 Node=NODE; N=Node(:,1); %节点号 Type=Node(:,2); BR=Branch K=Branch(:,6); n=length(N); nbr=length(K); %将支路信息保存在B冲%% %支路变压器变比， 0代表没有 %节点数 %支路数 Total_of_Bus1=size(NODE); %%取%节点矩阵的行和列 Total_of_Bus=Total_of_Bus1(1,1) %%bu 矩阵的行数即 节点数 Total_of_Branch1=size(Branch); %%%取%支路矩阵的行和列 Total_of_Branch=Total_of_Branch1(1,1); %%支路 branch 矩阵行数即 支路数 Z=zeros(Total_of_Bus1); format short %排序标志位 %卩节点标志位 %P\节点标志位 b=1; pq=0; pv=0; ph=0; % % %平衡节点标志位 for a=1:Total_of_Bus if NODE(a,2) == 0 Z(b,:)=NODE(a,:); b=b+1; pq=pq+1; end end 统计PV节点数2代表pv节点 for a=1:Total_of_Bus if NODE(a,2) == 2 Z(b,:)=NODE(a,:); b=b+1; pv=pv+1; end end 统计平衡节点数 3 代表平衡节点 for a=1:Total_of_Bus if NODE(a,2) == 3 Z(b,:)=NODE(a,:); b=b+1; ph=ph+1; end end Z Z2=Z; %%%%%%%%%%%%%%%%%%%%%将%节%点%%进%行%重%新%排%序%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% mm=zeros(n,1); for i=1:n mm(i,1)=i; end Z1(:,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 end end Branch(:,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(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 end end %%%%%%%%%%%%%%%%%求%自%导%纳%%%%%%%%%%%%%%%%%%%%%%%%%%