北理工电力系统潮流计算.docxVIP

电力系统作业 班级： 姓名： 学号： 实验名称 用matlab编写潮流计算程序 实验目的 掌握潮流计算中的牛顿－拉夫逊迭代法。 结合第三章所学加深对潮流计算的理解。 3、实验内容 编写潮流计算程序，计算下图等效电力系统网络的潮流。 其中各节点已进行编号，节点1为平衡节点， 运行时其电压为242kV。 4、设计思路及程序设计 4.1 设计思路 (1)形成节点导纳矩阵Y (2)设各节点电压初值e,f (3)将各节点电压的初值代入相关公式求各节点的不平衡量ΔP、ΔQ和ΔU (4)将各节点电压初值代入相关公式求修正方程的系数矩阵—雅可比矩阵的各个元素。 (5)解修正方程式的系数矩阵—求各节点电压的变化量，即修正量Δe、Δf。 (6)计算各节点电压的新值U。 (7)运用各节点电压的新值自第三步开始进入下一次迭代。 (8)计算平衡节点功率和线路功率。 4.2 程序设计 y12=1/(5.9+1i*31.5); y23=1/(0.8+1i*23); y34=1/(65+1i*100); y45=1/(65+1i*100); y51=1/(3+1i*110); Y=[y12+y51,-y12,0,0,-y51; -y12,y12+y23,-y23,0,0; 0,-y23,y23+y34,-y34,0; 0,0,-y34,y34+y45,-y45; -y51,0,0,-y45,y45+y51;]; %定义节点导纳矩阵 U=[0;242;0;220;0;220;0;220;0;220]; %电压以【f1,e1…f5,e5】形式存储 UL=[242+0*1i,220+0*1i,220+0*1i,220+0*1i,220+0*1i]; %电压以【e1+f1*i…e5+f5*i】形式存储便于之后的功率计算 U2=[0;0;0;0;0;0;0;0]; %电压修正量以【Δf1, Δe1…Δf5, Δe5】形式存储 P=[0;10;-180;-100;-50;-30;40;30]; %注入功率以【P2,Q2…P5,Q5】形式存储 P1=[0;0;0;0;0;0;0;0]; %各节点功率以【P2,Q2…P5,Q5】形式存储 P2=[0;0;0;0;0;0;0;0]; %各节点功率不平衡量以【ΔP2,ΔQ2…ΔP5,ΔQ5】形式存储，便于与雅可比矩阵进行运算 S=zeros(5,5); %定义了各线路功率 J=zeros(8,8); %定义了雅可比矩阵 I=[0;0;0;0;0]; %定义了各节点注入电流 a=zeros(5,5); %定义了各节点注入电流实部 b=zeros(5,5); %定义了各节点注入电流虚部 M1=[0,0,0,0]; %用来显示迭代过程中各节点功率不平衡量M2=[0,0,0,0]; %用来显示迭代过程中各节点电压修正量 M3=[0,0,0,0]; %用来显示迭代过程中各节点电压 for m=1:1:10 %求解各节点功率 for k=3:2:9 for j=1:1:5 P1(k-2,1)=P1(k-2,1)+U(k+1,1)*(real(Y((k+1)/2,j))*U(2*j,1)-imag(Y((k+1)/2,j))*U(2*j-1,1))+U(k,1)*(real(Y((k+1)/2,j))*U(2*j-1,1)+imag(Y((k+1)/2,j))*U(2*j,1)); %求解各节点功率实部 P1(k-1,1)=P1(k-1,1)+U(k,1)*(real(Y((k+1)/2,j))*U(2*j,1)-imag(Y((k+1)/2,j))*U(2*j-1,1))-U(k+1,1)*(real(Y((k+1)/2,j))*U(2*j-1,1)+imag(Y((k+1)/2,j))*U(2*j,1)); %求解各节点功率虚部 end end P2=P-P1; %求解各节点功率不平衡量 for j=2:1: