（
二 〇 一 一 年 十 二 月
（2011-2012学年第一学期）
电力系统计算机辅助潮流计算
实验报告
学 号： 2007
学生姓名：
学 院：电力学院 系 别：电力系
专 业：电气工程及其自动化 班 级： 授课教师：
1、实验目的：
了解计算机潮流分析的基本原理、主要步骤；掌握节点导纳矩阵形成和修改的方法，掌握数据处理的基本方法；熟悉Matlab运行环境，了解Matlab基本编程语句和语法；运用潮流分析程序对给定网络的运行方式做潮流分析，并初步分析计算结果
2、实验要求：
通过预习，对计算机潮流分析基本理论有深入了解；为程序准备必要的、准确的原始数据；熟悉Matlab运行环境，输入潮流程序，上机独立完成程序的调试，给出潮流分析的结果并按要求绘制潮流分布图
3、实验内容：
输入网络参数，包括节点号、节点导纳矩阵、节点功率等；输入潮流程序、调试并输出结果，绘制潮流分布图
4、实验步骤：
1、熟悉原始资料：根据计算要求，整理数据，包括：计算网络中线路、变压器的参数、形成节点导纳矩阵；表示各节点的注入功率。（以上数据均采用有名值计算）
2、读通潮流程序：完成程序的解释和说明，必要时附加对应的公式和程序语言的说明
3、上机调试：熟悉Matlab的运行环境，准确输入原始数据、节点编号、节点注入功率等信息
4、整理计算结果：根据计算结果作电网潮流分布图原始网络：
5、实验数据及处理：
一、实验程序：
clear
G(1,1)=;B(1,1)=;G(1,2)=;B(1,2)=;G(1,3)=0;B(1,3)=0;
G(1,4)=;B(1,4)=;G(1,5)=0;B(1,5)=0;
G(2,1)=;B(2,1)=;G(2,2)=;B(2,2)=;G(2,3)=;
B(2,3)=5;G(2,4)=;B(2,4)=5;G(2,5)=-5;B(2,5)=15;
G(3,1)=0;B(3,1)=0;G(3,2)=;B(3,2)=5;G(3,3)=;
B(3,3)=;G(3,4)=-10;B(3,4)=30;G(3,5)=;B(3,5)=;
G(4,1)=;B(4,1)=;G(4,2)=;B(4,2)=5;G(4,3)=-10;
B(4,3)=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(3)=; ps(4)=;qs(4)=; %设迭代初值
k=1;precision=1 %设迭代次数和精度N1=4; %PQ节点数while precision> %判断是否满足精度要求
delt(5)=0;u(5)=; %给定平衡节点编号
for m=1:N1
for n=1:N1+1
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)^2*G(m,m);
J(m,m)=pi(m)-u(m)^2*G(m,m); L(m,m)=qi(m)-u(m)^2*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)-del t(n)));
J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-de lt(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
551j j i I Y U ==∑
end
S5=U(5)*sum(conj(I)) %求平衡节点的注入功率*
555S V I = for m=1:N1+1 for n=1:N1+1
S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));% end end
end %求节点i,j 节点之间的功率，方向为由i 指向j, *
ij ij i S V I = S %显示支路功率