电力系统自动化课程设计内蒙古科技大学.docVIP

电力系统自动化课设 题 目：网络潮流计算(极坐标) 指导教师：杨培宏 专 业：电气工程及其自动化 班 级： 学 号： 学生姓名： 等值电路图表示） MATLAB程序 G(1,1)=2.9050;B(1,1)=-11.5031; G(1,2)=0;B(1,2)=5.3182; G(1,3)=-1.6600;B(1,3)=3.1619; G(1,4)=-1.2450;B(1,4)=2.3714; G(2,1)=0;B(2,1)=5.3182; G(2,2)=0;B(2,2)=-4.7138; G(2,3)=0;B(2,3)=0; G(2,4)=0;B(2,4)=0; G(3,1)=-1.6600;B(3,1)=3.1619; G(3,2)=0;B(3,2)=0; G(3,3)=2.4900;B(3,3)=-4.7040; G(3,4)=-0.8300;B(3,4)=1.5809; G(4,1)=-1.2450;B(4,1)=2.3714; G(4,2)=0;B(4,2)=0; G(4,3)=-0.8300;B(4,3)=1.5809; G(4,4)=2.0750;B(4,4)=-3.9092; Y=G+j*B delt(1)=0;delt(2)=0;delt(3)=0; u(1)=1.0;u(2)=1.0;u(3)=1.05; p(1)=0;q(1)=0;p(2)=-0.5;q(2)=-0.3;p(3)=0.2;u(3)=1.05;p(4)=-0.15;q(4)=-0.1; %置迭代次数k=0，并计算节点功率的不平衡量 ， 。 k=0;precision=1; k,delt,u N1=3; while precision0.00001 delt(4)=0;u(4)=1.05; for m=1:N1 for n=1:N1+1 pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); end pp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt); end pp,qq %/////////////////////////////////////////////////////// %求取雅可比矩阵元素 （m=n 时）。 for m=1:N1 for n=1:N1+1 h0(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(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(n))); end H(m,m)=sum(h0)-u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m))); N(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m)) +B(m,m)*sin(delt(m)-delt(m))); J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m))); L(m,m)=sum(l0)+2*u(m)^2*B(m,m)+u(m)^2*(G(m,m)*sin(delt(m)-delt(m)) -B(m,m)*cos(delt(m)-delt(m))); end for m=1:N1-1 JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m); JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m); end for m=N1:N1 JJ(2*m-1,2*m-1)=H(m,m); end