MATLAB潮流程序IEEE14直角坐标牛拉法.doc

MATLAB潮流程序(IEEE14直角坐标牛拉法)clearbaseMVA=100;%功率基值%%=Data1(:,1);%节点号Vtype=Data1(:,5);%节点类型Pload=Data1(:,8);%负载有功Qload=Data1(:,9);%负载无功Pgen=Data1(:,10);%发电机发出有功Qgen=Data1(:,11);%发电机发出无功Vset=Data1(:,13);%电压设定点Qsh=Data1(:,17);%并联电容电纳标幺值%%=Data2(:,1);JJ=Data2(:,2);%支路端点号Ltype=Data2(:,5);%线路类型R=Data2(:,6);%两点间电阻X=Data2(:,7);%两点间电抗B=Data2(:,8)/2;%线路对地电纳K=Data2(:,14);%变压器非标准变压比%%求导纳矩阵Yy1=zeros(14);y2=zeros(14);y3=zeros(14);lin=length(II);%支路数forx=1:linswitchLtype(x)case1y1(II(x),JJ(x))=1/(R(x)+i*X(x));y1(JJ(x),II(x))=y1(II(x),JJ(x));y3(II(x),JJ(x))=i*B(x);y3(JJ(x),II(x))=i*B(x);case2y1(II(x),JJ(x))=1/((R(x)+i*X(x))*K(x));y1(JJ(x),II(x))=y1(II(x),JJ(x));y2(II(x),JJ(x))=(1-K(x))/((R(x)+i*X(x))*K(x)^2);y2(JJ(x),II(x))=(K(x)-1)/((R(x)+i*X(x))*K(x));endendclearxY=zeros(14);forx=1:14Y(x,x)=sum(y1(x,:))+sum(y2(x,:))+sum(y3(x,:))+i*Qsh(x);endclearx;Y=Y-y1;G=real(Y);B=imag(Y);%%设电压初值U=Vset;e=real(U);f=imag(U);%%Ps=zeros(1,14);Qs=zeros(1,14);D=ones(26,1);forx=1:14Ps(x)=(Pgen(x)-Pload(x))/baseMVA;Qs(x)=(Qgen(x)-Qload(x))/baseMVA;endclearx;N=0;Jacbi=zeros(26);whilemax(abs(D))>=2:14%节点功率及电压不平衡量switchVtype(x)case1%PQ节点D(2*x-3)=Ps(x)-e(x)*(G(x,:)*e-B(x,:)*f)-f(x)*(G(x,:)*f+B(x,:)*e);D(2*x-2)=Qs(x)-f(x)*(G(x,:)*e-B(x,:)*f)+e(x)*(G(x,:)*f+B(x,:)*e);case2%PV节点D(2*x-3)=Ps(x)-e(x)*(G(x,:)*e-B(x,:)*f)-f(x)*(G(x,:)*f+B(x,:)*e);D(2*x-2)=Vset(x).*Vset