电力系统课程设计--高压输电网潮流的计算机算法程序设计

pe==3,%求平衡节点数nSW=nSW+1;SW(nSW,:)=bus(I,:);elseiftype==2,%求PV节点数nPV=nPV+1;PV(nPV,:)=bus(I,:);else%求PQ节点数nPQ=nPQ+1;PQ(nPQ,:)=bus(I,:);endendbus=[PQ;PV;SW];%按PQPVSW顺序排列的节点数据矩阵newbus=[1:nb]';nodenum=[newbusbus(:,1)];%形成节点号对应表bus(:,1)=newbus;forI=1:nl%根据节点号对应表重新排列线路矩阵节点编号forJ=1:2fork=1:nbifline(I,J)==nodenum(k,2)line(I,J)=nodenum(k,1);breakendendendend%******************建立节点导纳矩阵*****************functionY=ybus(bus,line)globalnl;globalnb;globalbus;globalline;globalY;globalmyf;Y=zeros(nb,nb);%对导纳矩阵赋初值0fork=1:nlI=line(k,1);%读入线路参数J=line(k,2);Zt=line(k,3)+j*line(k,4);ifJ~=0%接地支路分母为零Yt=1/Zt;%非接地支路为阻抗的倒数endYm=line(k,5)+j*line(k,6);%求YmK=line(k,7);%求Kif(K==0)&(J~=0)%普通线路:K=0;Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J);endif(K==0)&(J==0)%对地支路:K=0,J=0,R=X=0;Y(I,I)=Y(I,I)+Ym;end