《MATLAB潮流程序(IEEE14直角坐标牛拉法)》由会员分享，可在线阅读，更多相关《MATLAB潮流程序(IEEE14直角坐标牛拉法)（10页珍藏版）》请在金锄头文库上搜索。
1、%功率基值%节点号%节点类型%负载有功%负载无功%发电机发出有功%发电机发出无功%电压设定点%并联电容电纳标幺值%支路端点号%线路类型%两点间电阻%两点间电抗MATLAB潮流程序（IEEE14直角坐标牛拉 法） clearbaseMVA=100;%读Datal中数 据load Datal .txtBus=Data1(:,1);Vtype=Data1(:,5);Pload=Data1(:,8);Qload=Data1(:,9);Pgen=Data1(:,10);Qgen=Data1(:,11);Vset=Data1(:,13);Qsh=Data1(:,17);%读Data2中数 据load Da
2、ta2.txtll=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(ll);% 支路数for x=1 Jin switch Ltype(x) case 1 y1(ll(x),JJ(x)=1/(R(x)+i*X(x);y1(JJ(x),ll(x)=y1(ll(x),JJ(x);y3(ll(x),JJ(x)=i
3、*B(x);y3(JJ(x),ll(x)=i*B(x);case 2 y1(ll(x),JJ(x)=1/(R(x)+i*X(x)*K(x);y1(JJ(x),ll(x)=y1(ll(x),JJ(x);y2(ll(x)JJ(x)=(1-K(x)/(R(x)+i*X(x)*K(xF2);y2(JJ(x)JI(x)=(K(x)-1)/(R(x)+i*X(x)*K(x);end endclear xY=zeros(14);for x=1:14Y(x,x)=sum(y1 (x,:)+sum(y2(x5:)+sum(y3(x,:)+i*Qsh(x);end clear x;Y=Y-y1;G=real(Y)
4、;B=imag(Y);%设电压初值U=Vset;e=real(U);f=imag(U);%Ps=zeros(1,14);Qs=zeros(1,14);D=ones(26,1);for x=1:14Ps(x)=(Pgen(x)-Pload(x)/baseM VA;Qs(x)=(Qgen(x)-Qload(x)/baseMVA;endclear x;N=0;Jacbi=zeros(26);while max(abs(D)0.000001for x=2:14电压不平衡量节点功率及switch Vtype(x)case 1%PQ 节点D(2*x-3)=Ps(x)-e(x)*(G(x;)*e-B(x3:
5、)*f)-f(x)*(G(x5:)*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(x5:)*e-B(x,:)*f)-f(x)*(G(xJ:)*f+B(x5:)*e);D(2*x-2)=Vset(x)广Vset(x)-(e(x)八 2+f(x)八 2);e ndendclear mfor 1=2:14克比求雅矩阵 for J=2:14if l=J%非对角元素Jacbi(2*l-3),(2*J-3)=B(lJJ)*e(l)-
6、G(l,J)*f(l);Jacbi(2*l-3),(2*J-2)=-(G(IJ)*e(l)+B(l,J)*f(l);switch Vtype(l)case 1%PQ 节点Jacbi(2*l-2),(2*J-3)=G(l3J)*e W+B(l,J)pi);Jacbi(2%2)；(2J2)=B(l,J)P(l)-G(l,J)*f(l);case2%PV 节点Jacbi(2*l-2),(2*J-3)=0;Jacbi(2*l-2),(2*J-2)=0;endelse%对角元素Jacbi(2*l-3J2*J-3)=-(G(l5:)*f+B(l,:)*e)+B(l,l)*e(l)-G(IJ)*f(l);J
7、acbi(2*l-352*J-2)=-(G(l5:)*e-B(l,:)*f)-G(IJ)*e(l)-B(l,l)*f(l);switch Vtype(l)case 1%PQ 节点 Jacbi(2%232*J-3)=-(G(l/)*e-B(l/)*f)+G(IJ)*e(l)+B(IJ)*f(l);Jacbi(2*l-2,2*J- 2)=(G(l/)*f+B(l,:)*e)+B(IJ)*e(l)-G(IJ)*f(l);case22*f(l)；Jacbi(2*l-2,2*J-2)=-2*e(l);endendendendclear I J;Deta=-inv(Jacbi)*D;正方程forx=2:1
8、4初值 f(x)=f(x)+Deta(2*x-3),1); e(x)=e(x)+Deta(2*x-2),1);endclear x;%PV 节点 Jacbi(2%2,2J-3)=%修%新电压U=e+i*f:N=N+1;endN=N-1;S0=U(1)*(conj(Y(1,:)*conj(U); 衡节点功率S1=zeros(20,1);%始端功率%末S2=zeros(20,1);端功率 for x=1:20S1(x)=U(ll(x)*(conj(U(ll(x)*(conj(y2(ll(x)+y3(ll(x)+i*Qsh(ll(x)+(conj(U(ll(x)- conj(U(JJ(x)*conj
9、(y1(ll(x)5JJ(x);S2(x)=U(JJ(x)*(conj(U(JJ(x)*(conj(y2(JJ(x)+y3(JJ(x)+i*Qsh(JJ(x)+(conj(U(JJ(x) )-conj(U(ll(x)*conj(y1(ll(x),JJ(x);end detaS=S1+S2;%线路损耗功率%电%相Vabs=abs(U);压幅值Angle=atan(f./e)*180/Pi ;角%显示数据dispC迭代次数N=f);disp(N);dis PC各节点电压dis PC节点幅值相角)；disP(Bus Vabs Angle);dis PC平衡节点功率disp(SO);dis PC线路功率J;disPf HJJ始端功率末端功率线路损耗功率)；disP(ll JJS1 S2 detaS);Datal及Data2数据文件也已上传文件名就是Datal和Data2
