资源描述：
《《matlab潮流计算》word版》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、%线路数据Line=[1-节点编号2-线路首端节点号3-线路末端节点号4-支路电阻5-支路电抗6-支路电纳(注意：此处取的是B/2)]Line=[1120.040.250.252130.10.3503230.080.300.25];%变压器数据transform=[1-支路编号2-支路首节点编号3-支路末节点编号4-支路电阻(p.u.)5-支路电抗(p.u.)6-变压器变比(p.u.)]transform=[12400.0151.0523500.031.05];%数据预处理Nbus=5;%节点数nline=size(Line,1);%线路个数ntrans=size(transform,1);%
2、变压器个数slack=5;%平衡节点号Npq=4;%PQ节点的个数Ss=0;%计算节点导纳矩阵Y=zeros(Nbus);ifnline>=1%判断是否存在线路fork=1:nline%以下处理线路t1=Line(k,2);t2=Line(k,3);b2=Line(k,6);%分别取出线路的首端节点编号t1、末端节点编号t2和对地电纳b2Yl=1/(Line(k,4)+j*Line(k,5));%计算线路的支路电导YlY(t1,t1)=Y(t1,t1)+Yl+j*b2;%修正第k条线路首端节点的自导纳Y(t1,t2)=Y(t1,t2)-Yl;%修正第k条线路首端节点与末端节点之间的互导纳Y(t
3、2,t1)=Y(t2,t1)-Yl;%修正第k条线路末端节点与首端节点之间的互导纳Y(t2,t2)=Y(t2,t2)+Yl+j*b2;%修正第k条线路末端节点的自导纳endendifntrans>=1%判断是否存在变压器fork=1:ntrans%以下处理变压器t1=transform(k,2);t2=transform(k,3);t3=transform(k,6);%分别取出变压器的首端节点编号t1、末端节点编号t2和变比t3Yt=1/(transform(k,4)+j*transform(k,5));Yt1=Yt/t3;Yt2=Yt*(1-t3)/(t3*t3);Yt3=Yt*(t3-1)
4、/t3;Y(t1,t1)=Y(t1,t1)+Yt1+Yt2;Y(t1,t2)=Y(t1,t2)-Yt1;Y(t2,t1)=Y(t2,t1)-Yt1;Y(t2,t2)=Y(t2,t2)+Yt1+Yt3;endendG=real(Y);B=imag(Y);%区分节点导纳矩阵的实部和虚部GB%赋初值delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;u(1)=1;u(2)=1;u(3)=1;u(4)=1;p(1)=-0.30;q(1)=-0.18;p(2)=-0.55;q(2)=-0.13;p(3)=0;q(3)=0;p(4)=0.5;q(4)=1.10;p(5)=0.8
5、;q(5)=0.50;k=0;precision=1;Npq=4;%Npq分别是网络中的PQ节点数%[Unbalance]=-[Jacobi][Correction]whileprecision>0.00001%设定误差上限，判断是否继续迭代u(5)=1.06;delt(5)=0;%设定平衡节点电压相角与幅值k;u;delt;form=1:Npqforn=1:Nbuspt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));%由节点电压求得的PQ节点注入有功功率qt(n)=u(m)*u(n)*(G(m,n
6、)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));%由节点电压求得的PQ节点注入无功功率endUnbalance(2*m-1)=p(m)-sum(pt);%计算PQ节点有功功率不平衡量Unbalance(2*m)=q(m)-sum(qt);%计算PQ节点无功功率不平衡量end%[Unbalance]是节点不平衡量矩阵form=1:Npqforn=1:Nbush0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));n0(n)=-u(m)*u(n)*(G
7、(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt
