电力系统计算机辅助潮流计算.docx

文档介绍：电力系统计算机辅助潮流计算 实验报告 (
学生姓名:
学 院：电力学院
系 别：电力系
专 业：电气工程及其自动化
二0 年十二月
（2011-2012学年第一学期）
1、 实验目的：
了解计算机潮流分析的基本原理、主要=30;G(4,4) = ;B(4,4) = ;G(4,5)=0;B(4,5)=0;
G(5,1)=0;B(5,1)=0;G(5,2)=-5;B(5,2)=15;G(5,3) = ;B(5,3) = ;
G(5,4)=0;B(5,4)=0;G(5,5) = ;B(5,5) = ;
Y=G+j*B %形成节点导
纳矩阵
delt(1)=0;delt(2)=0;delt(3)=0; delt(4)=0;
u(1) = ;u(2) = ;u(3) = ;u(4) = ; ps(1) = ;qs(1) = ;ps(2) = ;qs(2) = ;
ps( 3) = ;qs ⑶=；ps⑷= ;qs ⑷=;
%设迭代初值
k=1;precision=1
和精度
N1=4;
点数
while precision>
delt(5)=0;u(5) = ;
号
for m=1:N1
for n=1:N1+1
%设迭代次数
%PQ节
%判断是否满足精度要求
%给定平衡节点编
pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)
-delt(n)));
qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)
-delt(n)));
end
pi(m)=sum(pt);qi(m)=sum(qt);
入功率
%计算PQ节点的注
dp(m)=ps(m)-pi(m);
dq(m)=qs(m)-qi(m) ; %计算PQ节点的功率不
平衡量
end
for m=1:N1
for n=1:N1
if m==n
H(m,m)=-qi(m)-u(m)八2*B(m,m); N(m,m)=pi(m)+u(m)”*G(m,m);
J(m,m)=pi(m)-u(m)”*G(m,m); L(m,m)=qi(m)-u(m)”*B(m,m);
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);
else
H(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m )-delt(n)));
J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt( m)-delt(n)));
N(m,n)=-J(m,n);L(m,n)=H(m,n);
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);
JJ(2*m,2*n-1)=J(m,n); JJ(2*m,2*n)=L(m,n);
End
end
end %计算jocbi各项，并放入统一矩阵JJ中，对JJ下标统
—编号
JJ
for m=1:N1
PP(2*m-1)=dp(m);
PP(2*m)=dq(m);
End %按统一矩阵形成功率不
平衡
uu=inv(JJ)*PP';precision=max(abs(uu)); % 判断是
否收敛 for n=1:N1
delt(n)=delt(n)+uu(2*n-1);
u(n)=u(n)+uu(2*n)*u(n); %将结果分解为电压幅值
和角度
end %求解修正方程，得电压幅值变化量(标幺值)和角度
变化量
k=k+1;
end
for n=1:N1+1
U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));
end
for m=1:N1+1
I(m)=Y(5,m)*U(m); %求平衡节点的注入电流
/5=U Y U i=1
end
S5=U(5)*sum(conj(I)) %求平衡节点的注入功率
*
S = V5 I
for m=1:N1+1
for n=1:N1+1
S(m,n)=U(m)*(conj(U(m))-conj(U