电力系统稳态分析课程设计---两机五节点网络潮流计算-牛拉法

m)-dt(n))); end

disp('各节点无功率不平衡量为：') qq(m)=q(m)-sum(qt) end for m=2:5

qq1=flipud(rot90(qq(2:5))); qq2=qq1(1:4,:);

ut=-b2*qq2/u4;

disp( '各节点电压的修正量：')

ut

disp( '各节点电压的大小为：') u(m)=u(m)+ut(m-1) end for n=1:5

U(n)=u(n)*(cos(dt(n))+j*sin(dt(n))); end for n=1:5

I(n)=Y(1,n)*U(n); end

disp('各节点功率s为：') S1=U(1)*sum(conj(I))

13

for m=1:5 for n=1:5

S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));

end end

disp('各节点功率sij为：') S

wucha=max(abs(ut)) k=k+1 End

第四节 程序运行结果显示 节点总数为： N = 5 平衡节点为： ans = 1 PQ节点为：

JD = 2 3 4 5

Y =6.2500 -18.7500i -5.0000 +15.0000i -1.2500 + 3.7500i 0 0

-5.0000 +15.0000i 10.8340 -32.5000i -1.6670 + 5.0000i

14

-1.6670 + 5.0000i -2.5000 + 7.5000i

-1.2500 + 3.7500i -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0

0 -1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i -1.2500 + 3.7500i

0 -2.5000 + 7.5000i -1.2500 + 3.7500i 3.7500 -11.2500i 各节点有功率不平衡量为： dp =0 0.5000

各节点有功率不平衡量为： dp = 0 0.5000 -0.3750 各节点有功率不平衡量为：

dp = 0 0.5000 -0.3750 -0.4000 各节点有功率不平衡量为：

dp = 0 0.5000 -0.3750 -0.4000 -0.6000 dtt =-0.0370 -0.0855 -0.0918 -0.1086 各节点电压的相位角dt为：

dt =0 -0.0370 0 0 0 dtt =-0.0370 -0.0855 -0.0918 -0.1086 各节点电压的相位角dt为：

dt = 0 -0.0370 -0.0855 0 0 dtt =-0.0370 -0.0855 -0.0918 -0.1086

15

0 各节点电压的相位角dt为：

dt = 0 -0.0370 -0.0855 -0.0918 0 dtt =-0.0370 -0.0855 -0.0918 -0.1086 各节点电压的相位角dt为：

dt = 0 -0.0370 -0.0855 -0.0918 -0.1086 各节点无功率不平衡量为： qq = 0 1.2119 各节点无功率不平衡量为： qq = 0 1.2119 -0.0773 各节点无功率不平衡量为：

qq = 0 1.2119 -0.0773 -0.1919 各节点无功率不平衡量为：

qq =0 1.2119 -0.0773 -0.1919 -0.3196 各节点电压的修正量： ut =

0.0395 0 0 0 0.0081 0 0 0 0.0064 0 0 0 0.0001 0 0 0 各节点电压的大小为：

u =1.0600 1.0395 1.0000 1.0000 1.0000 各节点电压的修正量：

16