《电力系统分析》课程设计-极坐标表示的牛顿拉夫逊法潮流计算程序设计

目录
1任务书 (2)
2.模型简介及等值电路 (3)
3.设计原理 (4)
4.修正方程的建立 (6)
5.程序流程图及MATLAB程序编写 (8)
6.结果分析 (14)
7.设计总结 (17)
8.参考文献 (17)
《电力系统分析》
课程设计任务书
题目极坐标表示的牛顿拉夫逊法潮流计算程序设计
学生姓名学号专业班
级
电气工程及其
自动化
设计内容与要求1. 设计要求
掌握MATLAB语言编程方法；理解和掌握运用计算机进行潮流计算的基本算法原理；针对某一具体电网，进行潮流计算程序设计。

其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解，培养学生运用所学知识分析和解决问题的能力。

2. 内容
1）学习并掌握MATLAB语言。

2）掌握变压器非标准变比概念及非标准变比变压器的等值电路。

掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。

3）掌握电力系统功率方程、变量和节点分类。

4）掌握利用极坐标表示的牛-拉法进行潮流计算的方法和步骤。

5）选择一个某一具体电网，编制程序流程框图。

6）利用MATLAB语言编写该模型的潮流计算程序，并上机调试程序，对计算结果进行分析。

7）整理课程设计论文。

起止时间2013 年7 月 4 日至2013 年7月10日指导教师签名年月日
系（教研室）主任
签名
年月日学生签名年月日
2 模型简介及等值电路 2.1 课程设计模型：模型3
电力网络接线如下图所示，各支路阻抗标幺值参数如下：Z12=0.02+j0.06，Z13=0.08+j0.24， Z23=0.06+j0.18， Z24=0.06+j0.12， Z25=0.04+j0.12， Z34=0.01+j0.03， Z45=0.08+j0.24， k=1.1。

该系统中，节点1为平衡节点，保
持1
1.060V j =+&为定值；节点2、3、4都是PQ 节点，节点5为PV 节点，给定的注入功率分别为:
20.200.20
S j =+，
3-0.45-0.15
S j =，
40.400.05
S j =--，
50.500.00S j =-+，5 1.10V =&。

各节点电压（初值）标幺值参数如下： 节点
1 2 3 4 5
Ui （0）=ei （0）+jfi （0）
1.06+j0.0
1.0+j0.0 1.0+j0.0
1.0+j0.0
1.1+j0.0
计算该系统的潮流分布。

计算精度要求各节点电压修正量不大于10-5。

图2-1
2.2 模型分析
节点1是平衡节点，节点2、3、4是PQ节点，节点5是PV节点。

由题可得等值电路模型中各节点之间的导纳：y12=5.000-j15.000，y13=1.2500-j3.7500，y22=0.2750-j0.8250，y23=1.667-j5.000，y24=3.333-j6.667，y25=2.7500-j8.2500，y34=10.0000-j30.0000，y55=-0.25+j0.75
2.3 等值电路模型
在图2-2中，将图2-1中的编号重新编排，节点2、3、4、5、1替换为1、2、3、4、5。

则各节点之间的导纳变为y12=1.667-j5，y13=3.333-j6.667，y14=2.75-j8.25，y15=5-j15，y52=1.25-j3.75，y23=10-j30，y34=1.25-j3.75，y11=0.275-j0.825，y44=-0.25+j0.75。

3 设计原理
本题采用了题目要求的牛顿－拉夫逊潮流计算的方法。

牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。

3.1潮流计算的定解条件
题中所给图表示一个五节点的简单电力系统，n 个节点电力系统的潮流方程的一般形式是
.
.
1
n i i
j ij j i
P jQ V Y V =-=∑ (1,2,3,...)i n =
或
.
.
1
n j i i i ij j V P jQ V Y =+=∑ (1,2,3,...)i n =
按变量的不同，一般将节点分为三种类型。

（1） PQ 节点
这类节点的有功功率和无功功率是给定的，节点(,)V δ是待求量。

通常变电所都是这一类型节点。

由于没有发电设备，故其发电功率为零。

属于这一类节点的有按给定有功、无功功率发电的发电厂母线和没有其他电源的变电所母线。

（2） PV 节点
这类节点有功功率P 和电压幅值V 是给定的，节点的无功功率Q 和电压的相位δ是待求量。

这类节点必须有足够的可调无功容量，用以维持给定的电压幅值，因此又称为电压控制节点。

一般选择有一定无功储备的发电厂和具有可调无功电源设备的变电所作为PV 节点。

（3） 平衡节点
在潮流分布算出以前，网络中的功率损耗是未知的，因此，网络中至少有一个节点的有功功率P 不能给定，这个节点承担了系统系统的有功功率平衡，故称之为平衡节点。

潮流计算时，一般只设一个平衡节点。

3.2潮流计算的约束条件
(1)所有节点电压必须满足
min max i i V V V ≤≤
(2)所有电源节点的有功功率和无功功率必须满足
min max
min max
}Gi Gi Gi Gi Gi Gi P P P Q Q Q ≤≤≤≤
(3)某些节点之间的电压应满足
max
i j i j
δδδδ-<-
3.3牛顿-拉夫逊的基本原理
设欲求解的非线性代数方程为
()0f x =
设方程的真实解为x ，则必有()0f x =。

用牛顿-拉夫逊法求方程真实解x 的步骤如下：
首先选取余割合适的初始估值(0)
x 作为方程
(0)()0f x =的解，若恰巧有(0)()0f x =，则方程的真实解即为(0)x x =若(0)()0f x ≠，则做下一步。

取
(1)(0)(0)x x x =+∆ 则(1)(0)(0)
()()f x f x x =+∆
其中(0)x ∆为初始估值的增量，即(0)
(1)(0)x
x x ∆=-。

设函数()f x 具有任意阶导数
(1)(0)(0)(0)'(0)(0)''(0)(0)2()()()()[()]/2...f x f x x f x f x x f x x =+∆=+∆+∆+ 若所取的
(0)
x ∆足够小，则含
(0)2()x ∆的项及其余的一切高阶项均可略去，并使其等于零，即：
(1)(0)(0)(0)'(0)(0)
()()()()0f x f x x f x f x x =+∆≈+∆= 故得
(0)(0)
'(0)
()
()f x x
f x ∆=-
可见，只要'(0)()f x ≠0，即可根据上式求出第一次的修正估值(1)x ，若恰巧有
(1)()f x =0，则方程的真实解即为(1)x x =。

若(1)
()0f x ≠，则用上述方法由(1)x 再
确定第二次的修正估值。

如此反复叠代下去，直到求得真实解x 为止。

设第K 次的估值为第（K+1）次的修正估值，则有
()(1)
()
'()
()()k k k k f x x
x
f x +=-
迭代过程的收敛数据为
1
()()k f x ε<或2
()k x ε∆<
其中，1
ε
，
2
ε
为预先给定的小正数。

4 修正方程的建立
极坐标表示的牛拉法修正方程为：
⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢
⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢
⎢⎢⎣
⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆∆∆∆∆∆∆δδδδn P
nn np n n n n pn pp P P P P n p
n p n p
n p n P U U U U H H N H N H H H N H N H J J L J L J H H N H N H J J L J L J H H N H N H P P Q P Q P 222111
22112211
22222221212222222121111212111111121211112211//
式中留出了（n-m ）行空格和（n-m ）列空列。

式中的有功，无功功率不平衡量Q P i
i
∆∆.
分别由式（3-1a ）,式（3-1b ）可得为
∑==+-=∆n
j j ij ij ij ij i i B G U U P P 1j i sin cos ）
（δδ ∑==--=∆n j j ij ij ij ij i i B G U U Q Q 1
j i cos sin ）
（δδ （3-1a ，b ）
而式中雅可比矩阵的各个元素则分别为 j i P H δ∂∂=
ij ；j j i ij U U P N ∂∂=；j i Q J δ∂∂=ij ; j j
i
ij U U Q L ∂∂= （3-2） 修正方程中将i U ∆改为/i i U U ∆只是为使公式（3-2）中个偏导数的表示形式上更相似，为求取这些偏导数，可将i P 、i Q 分别展开如下
()
()
2121cos sin sin cos j n
i i
ii i j ij ij ij ij j j i j n
i i
ii i j ij ij ij ij j j i
P U G U U G B Q U B U U G B δδδδ==≠==≠=++=-+-∑∑ （3-3a,b ）
计及
cos cos()
sin()sin ()sin sin()
cos()cos ()cos cos()
sin()sin ()
sin sin()cos()cos ()ij i j i j ij
j i j ij i j i j ij j i j ij i j i j ij i i j ij i j i j ij i i j δδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδ∂∂-⎫=
=-=∂-∂-∂∂-==--=-∂-∂-⎬
∂∂-==--=-∂∂-∂∂-==-=∂∂-⎪
⎪
⎪⎪
⎪⎪
⎪⎪
⎪
⎪⎭
（3-4）
j i ≠时，由于对特定的j ，只有该特定节点的j δ，从而特定的ij i j δδδ=-是
变量，由式（3-2）到式（3-4）可得
()()sin cos cos sin i ij i j ij ij ij ij j i ij i j ij ij ij ij i P H U U G B Q
J U U G B δδδδδδ∂⎫
=
=-⎪∂⎪
⎬∂⎪
==-+⎪∂⎭
（3-5a ）
相似的，由于对特定的j ，只有该特定节点的j U 是变量，可得
()()cos sin sin cos i ij i j ij ij ij ij i i
ij i j ij ij ij ij j
P N U U G B U Q L U U G B U δδδδ∂⎫
=
=+⎪∂⎪
⎬∂⎪
==-⎪∂⎭ （3-5b ） j=i 时，由于i δ是变量，所有ij i j δδδ=-都是变量，可得
()()1
1
sin cos cos sin j n
i
ii i j ij ij ij ij j i j i j n
i
ii i j ij ij ij ij j i j i P H U U G B Q J U U G B δδδδδδ==≠==≠⎫∂==--⎪
∂⎪⎪⎬∂⎪
==+⎪∂⎪⎭
∑∑ （3-5c ） 相似的，由于i U 是变量，可得
()()21
21
cos sin 2sin cos 2j n
i
ii i j ij ij ij ij i ii j i j i j n
i
ii i j ij ij ij ij i ii j i j i P N U U G B U G U Q L U U G B U B U δδδδ==≠==≠⎫∂==++⎪
∂⎪⎪⎬∂⎪
==--⎪∂⎪⎭∑∑ （3-5d ）
5 设计流程图及程序的编写
输入数据
形成节点导纳矩阵
设定节点起始计算电压δ U
置迭代次数
应用公式计算不平衡量△P△Q
雅克比矩阵J是否已经全部
形成，i>n?
按公式计算雅克比矩阵元素J
增大节点号，i=i+1
解修正方程，由△P,△Q和J计算电压修正量△u和△δ
求出△u和△δ
迭代是否收敛，△u，△δ﹤ε
计算平衡节点功率S和线路的功率S
结束
K=k+1
△u=U+△u △δ=δ+△δ
置节点号i=1
程序中用到的符号所代表的意义：
Y 代表导纳矩阵
JJ 代表雅克比矩阵
∆
pp代表有功功率的不平衡量P i
∆
qq 代表无功功率的不平衡量Q i
∆U
uu 代表各节点电压和相角的不平衡量δ∆
U 代表各节点的电压
S 代表线路的功率
Q4 代表平衡节点的注入无功功率
K 代表迭代次数
N1 代表PQ节点和PV节点的总数
m,n代表系统中的节点总数，把平衡节点标为最大号
S5 平衡节点功率
d 变压器对地导纳
%The following program for load calculation is based on MATLAB6.5 %以下部分为输入原始数据（到标示‘///’标志为止）。

%求取节点导纳矩阵
g(1,1)=0.275;b(1,1)=0.825;
g(1,2)=1.667;b(1,2)=-5;
g(1,3)=3.333;b(1,3)=-6.667;
g(1,4)=2.75;b(1,4)=-8.25;
g(1,5)=5;b(1,5)=-15;
g(2,1)=1.667;b(2,1)=-5;
g(2,2)=-0;b(2,2)=0;
g(2,3)=10;b(2,3)=-30;
g(2,4)=0;b(2,4)=0;
g(2,5)=1.25;b(2,5)=-3.75;
g(3,1)=3.33;b(3,1)=-6.667;
g(3,2)=10;b(3,2)=-30;
g(3,3)=0;b(3,3)=0;
g(3,4)=1.25;b(3,4)=-3.75;
g(3,5)=0;b(3,5)=0;
g(4,1)=2.75;b(4,1)=-8.25;
g(4,2)=0;b(4,2)=0;
g(4,3)=1.24;b(4,3)=-3.75;
g(4,4)=-0.25;b(4,4)=0.75;
g(4,5)=0;b(4,5)=0;
g(5,1)=5;b(5,1)=-15;
g(5,2)=1.25;b(5,2)=-3.75;
g(5,3)=0;b(5,3)=0;
g(5,4)=0;b(5,4)=0;
g(5,5)=0;b(5,5)=0;
d(1,4)=0.275-0.825*j;
d(4,1)=-0.25+0.75*j
d(1,5)=0;d(5,1)=0;
%求取节点导纳矩阵。

for m=1:5
for n=1:5
if m==n
G(m,m)=g(m,1)+g(m,2)+g(m,3)+g(m,4)+g(m,5);
B(m,m)=b(m,1)+b(m,2)+b(m,3)+b(m,4)+b(m,5);
else
G(m,n)=-g(m,n);
B(m,n)=-b(m,n);
end
end
end
Y=G+j*B
%//////////////////////////////////////////////////////////////
%//设定节点起始节点电压，并给出已知功率值。

%//下面将题中节点2、3、4、5、1分别替换为节点1、2、3、4、5即节点4为PV节点，节点5为平衡节点。

delt(1)=0;delt(2)=0;delt(3)=0;delt(4)=0;
u(1)=1.0;u(2)=1.0;u(3)=1.0;
p(1)=0.2;q(1)=0.2;
p(2)=-0.45;q(2)=-0.15;p(3)=-0.4;q(3)=-0.05;p(4)=-0.5;q(4)=0.0;
%置迭代次数k=0，并计算节点功率的不平衡量，。

k=0;precision=1;
k,delt,u
N1=4;
while precision>0.00001
u(4)=1.1; delt(5)=0;u(5)=1.06;
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
pp(m)=p(m)-sum(pt);qq(m)=q(m)-sum(qt);
end
pp,qq
%///////////////////////////////////////////////////////
%求取雅可比矩阵元素（m=n 时）。

for m=1:N1
for n=1:N1+1
h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
end
H(m,m)=sum(h0)-u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m))); N(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))
+B(m,m)*sin(delt(m)-delt(m)));
J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m))); L(m,m)=sum(l0)+2*u(m)^2*B(m,m)+u(m)^2*(G(m,m)*sin(delt(m)-delt(m))
-B(m,m)*cos(delt(m)-delt(m)));
end
for m=1:N1-1
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=N1:N1
JJ(2*m-1,2*m-1)=H(m,m);
end
%///////////////////////////////////////////////////////
%求取雅可比矩阵元素（m不等于n时）
for m=1:N1
for n=1:N1
if m==n
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);
end
end
end
for m=1:N1-1 %求前六行六列元素(非对角)
for n=1:N1-1
if m==n
else
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
for m=N1
for n=1:N1-1 %求取第七行的元素
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);
end
end
for n=N1
for m=1:N1-1 %求的第七列元素
JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m,2*n-1)=J(m,n);
end
end
%///////////////////////////////////////////////////////
%解修正方程式，由，和计算电压修正量和。

for m=1:N1-1
PP(2*m-1)=pp(m);PP(2*m)=qq(m);
end
for m=N1
PP(2*m-1)=pp(m);
end
uu=-inv(JJ)*PP';precision=max(abs(uu));uu
%///////////////////////////////////////////////////////
%若结果不收敛，执行下列语句。

for n=1:N1-1
delt(n)=delt(n)+uu(2*n-1);
u(n)=u(n)+uu(2*n);
end
for n=N1
delt(n)=delt(n)+uu(2*n-1);
end
k=k+1;
k,delt,u
end
%///////////////////////////////////////////////////////
%若结果收敛，计算各节点电压，平衡节点功率，PV节点功率和线路功率。

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);
end
S5=U(5)*sum(conj(I)); %平衡节点功率
%PV节点功率
for n=1:N1+1
q4(n)=u(4)*u(n)*(G(4,n)*sin(delt(4)-delt(n))-B(4,n)*cos(delt(4)-delt(n)));
end
Q4=sum(q4)
%线路功率a
for m=1:N1+1
for n=1:N1+1
S(m,n)=U(m)*(conj(U(m))*conj(d(m,n))+(conj(U(m))-conj(U(n)))*conj(-Y(m,n)));
end
end
%///////////////////////////////////////////////////////
%显示运行结果（至结束）。

Y
JJ
S
B
pp
qq
uu
U
k
Q4
S5
6 程序执行结果及分析
导纳矩阵Y
Y =
13.0250 -34.0920i -1.6670 + 5.0000i -3.3330 + 6.6670i -2.7500 + 8.2500i -5.0000 +15.0000i
-1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i
-3.3300 + 6.6670i -10.0000 +30.0000i 14.5800 -40.4170i -1.2500 + 3.7500i 0
-2.7500 + 8.2500i 0 -1.2400 + 3.7500i 3.7400 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500i
雅克比矩阵JJ
JJ =
-40.9565 -15.9241 5.8963 1.8286 7.9334 3.7268 10.0491
15.5240 -41.3566 -1.8286 5.8963 -3.7268 7.9334 -3.0796
5.8139 2.0756 -44.1859 -14.2290 34.2658 11.2955 0
-2.0756 5.8139 15.1290 -43.8859 -11.2955 34.2658 0
7.7421 4.1061 34.1899 11.5231 -46.3468 -16.3011 4.4148
-4.1061 7.7421 -11.5231 34.1899 17.1011 -46.2468 -1.4718
9.8870 3.5657 0 0 4.4149 1.4597 -14.3020
平衡节点功率
S5=1.2299-1.0678i
PV节点的注入无功功率
Q4=-0.6893
各节点的节点电压
U =
1.0966 - 0.0677i 1.0624 - 0.0881i 1.0663 - 0.0920i 1.0959 - 0.0945i 1.0600
迭代次数
K=6
迭代过程中各节点功率的不平衡量
K △P1 △Q1 △P2 △Q2 △P3 △Q3 △P4 △Q4
0 0.5000 2.7500 -3.750 0.0750 -0.2750 0.3250 -0.6375 -0.4125
1 -0.1013 -0.3711 0.0261 -0.0269 0.0501 -0.0454 0.0208 0.8565
2 -0.1620 -0.0428 0.001
3 0.0018 0.0017 -0.0013 0.0047 0.7059
3 -0.0016 -0.0042 -0.0002 0.000
4 0.0002 0.0001 0.000
5 0.6907
4 -0.152.3
*e-3 -0.0004
*e-3
0.0199
*e-3
0.0001
*e-3
0.0209
*e-3
0.0000
*e-3
0.0433
*e-3
0.6894
*e-3
5 -0.1416
*e-4 -0.0000
*e-4
0.0200
*e-4
0.0000
*e-4
0.0023
*e-4
0.0000
*e-4
0.0378
*e-4
0.6893
*e-4
迭代过程中的各节点电压的修正量
k △δ1△U1 △δ2 △U2 △δ3△U3 △δ4
0 -0.0648 0.1115 -0.0882 0.0757 -0.0918 0.0801 -0.0876
1 0.0031 -0.0114 0.0054 -0.0089 0.0056 -0.0090 0.0016
2 -0.0000 -0.001
3 0.0000 0.0007 -0.0000 -0.0008 -0.0001
3 -0.0004
*e-3 -0.1190
*e-3
-0.0004
*e-3
-0.0475
*e-3
-0.0005
*e-3
-0.5320
*e-3
-0.0010
*e-3
4 -0.0000
*e-4 -0.1073
*e-3
-0.0000
*e-4
-0.2950
*e-4
-0.0000
*e-4
-0.0351
*e-4
-0.0001
*e-4
5 -0.0000
*e-6 -0.9640
*e-6
-0.0000
*e-6
-0.1828
*e-6
-0.0000
*e-6
-0.2303
*e-6
-0.0001
*e-6 迭代过程中的各节点电压的修正量
k △δ1△U1 △δ2 △U2 △δ3△U3 △δ4
0 -0.0648 0.1115 -0.0882 0.0757 -0.0918 0.0801 -0.0876
1 0.0031 -0.0114 0.0054 -0.0089 0.0056 -0.0090 0.0016
2 -0.0000 -0.001
3 0.0000 0.0007 -0.0000 -0.0008 -0.0001
3 -0.0004
*e-3 -0.1190
*e-3
-0.0004
*e-3
-0.0475
*e-3
-0.0005
*e-3
-0.5320
*e-3
-0.0010
*e-3
4 -0.0000
*e-4 -0.1073
*e-3
-0.0000
*e-4
-0.2950
*e-4
-0.0000
*e-4
-0.0351
*e-4
-0.0001
*e-4
5 -0.0000
*e-6 -0.9640
*e-6
-0.0000
*e-6
-0.1828
*e-6
-0.0000
*e-6
-0.2303
*e-6
-0.0001
*e-6
迭代过程中各节点电压
K δ1 U1 δ2U2 δ3U3 δ4U4 δ5U5
0 0 1.0000 0 1.0000 0 1.0000 0 1.1000 0 1.0600
1 -0.0648 1.1115 -0.088
2 1.0757 -0.0918 1.0801 -0.0876 1.1000 0 1.0600
2 -0.0617 1.1002 -0.0828 1.0668 -0.0861 1.0711 -0.0861 1.1000 0 1.0600
3 -0.0617 1.0989 -0.0828 1.0661 -0.0861 1.0703 -0.0860 1.1000 0 1.0600
4 -0.0617 1.0987 -0.0828 1.0660 -0.0861 1.0703 -0.0861 1.1000 0 1.0600
5 -0.0617 1.0987 -0.0828 1.0660 -0.0861 1.0703 -0.0861 1.1000 0 1.0600
6 -0.061
7 1.0987 -0.082
8 1.0660 -0.0861 1.0703 -0.0861 1.1000 0 1.0600
各线路的功率S ij
j
i
1 2 3 4 5
1 0 0.1839 +
0.1398i 0.2968 +
0.1152i
0.5722 +
0.9064i
-0.8529 +
1.0305i
2 -0.1812 -
0.1319i 0 0.0687 -
0.1735i
0 -0.3375 +
0.1554i
3 -0.2917 -
0.1052i -0.0683 +
0.1744i
0 -0.0400 -
0.1193i
4 -0.5407 -
0.8120i 0 0.0407 +
0.1226i
0 0
5 0.8825 -
0.9416i 0.3472 -
0.1262i
0 0 0
7 设计总结
（1）首先要正确理解牛拉法潮流计算具体方法和步骤。

（2）本次课程设计的关键在于建立修正方程和对三类节点的处理以及运用Matlab进行正常的编程。

（3）此模型中变压器用∏形等值电路，潮流计算中迭代了6次得到了满足精度的各节点电压，编程时，将节点做以适当的调整，把平衡节点编号置为最大号，有利于程序的编程和理解。

8 参考文献
（1）《电气工程专业毕业设计指南电力系统分册》陈跃主编中国水利水电出版社。

（2）《MATLAB程序设计与应用》张智星主编清华大学出版社。