发布时间 : 星期二 文章(强烈推荐)两机五节点网络潮流计算方法牛拉法和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