(强烈推荐)两机五节点网络潮流计算方法牛拉法和pq法_电力系统稳态分析毕业论文

电力系统稳态分析课程设计

2.2.1 程序的编写

%说明：为了使节点按照先PQ，再PV节点，最后平衡节点的次序编号，以便与公式对照，节点1与节点5对调。 %节点阻抗矩阵

Z=[0,0.04+0.12i,0,0.08+0.24i,0;0.04+0.12i,0,0.06+0.18i,0.06+0

.18i,0.02+0.06i;0,0.06+0.18i,0,0.01+0.03i,0.08+0.24i;0.08+0.24i,0.06+0.18i,0.01+0.03i,0,0;0,0.02+0.06i,0.08+0.24i,0,0]; %求互导纳 for m=1:5 for n=1:5

if Z(m,n)==0 y(m,n)=0; else

y(m,n)=1Z(m,n); end end end %求导纳 for m=1:5 for n=1:5 if m~=n

Y(m,n)=-y(m,n); else

Y(m,n)=sum(y(m,:)); end end end %导纳矩阵

G=real(Y);B=imag(Y);

%计算各节点功率的不平衡量设U=E+jF ；Y=G+Bj；

Y

E=real(U); F=imag(U); %设S=P+Qj；

S=[-0.60-0.10i;0.20+0.20i;-0.45-0.15i;-0.40-0.05i;0] P=real(S);Q=imag(S);

k=0;C=1;

while C > 0.00001 E(5)=1.06;F(5)=0; for m=1:4 for n=1:5

%计算Pi,Qi,设Pi=Pt;Qi=Qt,按照书上的公式：

Pt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)));

Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n))-E(m)*(G(m,n)*F(n)+B(m,n)*E(n))); end

%设P,Q的改变量为dP,dQ dP(m)=P(m)-sum(Pt); dQ(m)=Q(m)-sum(Qt); end

%计算Hii,Nii,Jii,Lii,由公式4-41b 左侧公式实现,sum(Ai)，sum(Bi)用于实现公式中的sigerma从j到n的求和;

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

Bi(n)=G(m,n)*F(n)+B(m,n)*E(n);

Ai(n)=G(m,n)*E(n)-B(m,n)*F(n); end

H(m,m)=sum(Bi)-(B(m,m)*E(m)+G(m,m)*F(m))+2*G(m,m)*F(m); N(m,m)=sum(Ai)-(G(m,m)*E(m)-B(m,m)*F(m))+2*G(m,m)*E(m); J(m,m)=-2*B(m,m)*F(m)+sum(Ai)-(G(m,m)*E(m)-B(m,m)*F(m));

L(m,m)=-2*B(m,m)*E(m)-(sum(Bi)-(B(m,m)*E(m)+G(m,m)*F(m))); end

%设雅可比矩阵为JJ，以下语句用来实现雅可比矩阵中对角线上元素H N J L 的排列

for m=1:4

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=1:4 for n=1:4 if m==n