电力系统分析大作业-..doc

电力系统分析 程序设计 专业班级 姓 名 学 号 日 期 要求计算所给系统的潮流，设发电机G1的端电压为1 p.u.，发出的有功、无功可调；发电机G2的端电压为1 p.u.，按指定的有功P=0.5 p.u.发电，取ε=10-4。SB=100MVA，UB=Uav。 一牛顿—拉夫逊法 程序代码 clear; clc; %节点编号，5号为平衡节点，4号为PV节点，1、2、3号为PQ节点 %输入原始数据，求节点导纳矩阵 y=0; y(1,1) = 1/(0.03+0.1i)+1/(0.025+0.08i)+0.09i+0.07i; y(1,2) = -1/(0.025+0.08i); y(1,3) = -1/(0.03+0.1i); y(2,1) = -1/(0.025+0.08i); y(2,2) = 1/(0.02+0.06i)+1/(0.025+0.08i)+0.05i+0.07i+1/1.0522^2*1/(0.1905i); y(2,3) = -1/(0.02+0.06i); y(2,4) = 4.9889i; y(3,1) = -1/(0.03+0.1i); y(3,2) = -1/(0.02+0.06i); y(3,3) = 7.7523-28.7757i; y(3,5) = 4.9889i; y(4,2) = 4.9889i; y(4,4) = -5.2493i; y(5,3) = 4.9889i; y(5,5) = -5.2493i; Y=y; display(导纳矩阵为: ); Y %Y 为导纳矩阵 G=real(Y); B=imag(Y); %原始节点功率 Pg=zeros(1,5); Qg=zeros(1,5); Pd=zeros(1,5); Qd=zeros(1,5); S(1)=-0.8055-0.5320i; S(2)=-0.18-0.12i; S(3)=0; S(4)=0.5; S(5)=0; P=real(S); Q=imag(S); %赋初值 U = ones(1,5); a = zeros(1,5); dx=ones(7,1); fx=ones(7,1); count=0; %计算迭代次数 while (max(fx)1e-4) %收敛判据 for i=1:4 for j=1:4 H(i,j)=0; dP(i)=0; end end for i=1:3 for j=1:3 L(i,j)=0; dQ(i)=0; end end for i=1:4 for j=1:3 N(i,j)=0; end end for i=1:3 for j=1:4 K(i,j)=0; end end %P %失配功率表达式 for i=1:4 x=0; for j=1:5 %dP(i)=dP(i)-U(i)*U(j)*(G(i,j)*cos(a(i)-a(j)) + B(i,j)*sin(a(i)-a(j)) ); x=x+U(i)*U(j)*(G(i,j)*cos(a(i)-a(j))+B(i,j)*sin(a(i)-a(j))); end dP(i)=P(i)-x; end for i=1:3 y=0;