电力系统潮流计算第一题部分课程设计

课程设计(论文)

题目名称潮流计算课程设计

课程名称电力系统稳态分析

学生姓名徐玛丽

学号08４1229186

系、专业电气工程系08电力二班

指导教师黄肇王晓芳

２011年1 月5 日

邵阳学院课程设计(论文）任务书

注：1．此表由指导教师填写,经系、教研室审批,指导教师、学生签字后生效;

2.此表1式3份,学生、指导教师、教研室各1份

指导教师（签字）: 学生（签字)：

邵阳学院课程设计(论文）评阅表

学生姓名徐玛丽学号0841229186

系电气工程及其自动化专业班级08电力二班

题目名称潮流计算课程设计课程名称电力系统分析

二、指导教师评定

注：１、本表是学生课程设计(论文）成绩评定的依据，装订在设计说明书(或论文）的“任务书”页后面;

2、表中的“评分项目”及“权重”根据各系的考核细则和评分标准确定。

目录

第1章潮流计算课题及算法ﻩ错误!未定义书签。

1．1 潮流计算课题 .......................................... 错误!未定义书签。

1.2 极坐标下P-Q法的算法 .............................. 错误!未定义书签。

１.2.１节点导纳矩阵Y ........................... 错误!未定义书签。

１.２.2 简化雅可比矩阵B/和Ｂ/／ ............. 错误!未定义书签。

１.2.3 修正和迭代 .................................... 错误!未定义书签。第2章手工计算ﻩ错误!未定义书签。

第3章程序设计ﻩ错误!未定义书签。

3.1 流程图........................................................... 错误!未定义书签。

3．２潮流计算程序 .......................................... 错误!未定义书签。

3.3 潮流计算程序运行结果 .............................. 错误!未定义书签。结束语......................................................................... 错误!未定义书签。

第１章 潮流计算课题及算法

1.1 潮流计算课题

题目一:在图１所示的简单电力系统中，系统中节点1、2为PQ 节点，节点3为

PV 节点,节点4为平衡节点，已给定3.04.01j s --=,2.03.02j s --=,4.03=P ，

02.13=V ，05.14=V ， 04=θ，网络各元件参数的标幺值如表２所示，给定电压

的初始值如表2所示,收敛系数00001.0=ε。试求：

4

4θ∠V 1

12

jQ

图1 简单电力系统

表1 网络各元件参数的标幺值

表２ 各节点电压(初值)标幺值参数

（３)采用极坐标下的Q P -分解法计算图１网络的潮流分布。

1．２ 极坐标下P-Q 法的算法

1.２.1 节点导纳矩阵Y

根据题目提供的各节点的参数，求得节点导纳矩阵

Y

ii

=ij y j

i y ∑+0

y Y

ik

ik

-=

1.2.2 简化雅可比矩阵Ｂ／和B //

通过上一步的导纳矩阵,形成有功迭代和无功迭代的简化雅可比矩阵B ／和B ／

／

对雅可比矩阵进行三角分解，形成因子表，为后面进行修正方程计算作好准备。 1．2．3 修正和迭代

第一步，给定P Ｑ节点初值和各节点电压相角初值。

第二步，作第一次有功迭代，按公式计算节点有功功率不平衡量。 第三步，做第一次无功迭代,按公式计算无功功率不平衡量，计算时电压相角

最新的修正值。解修正方程式，可得各节点电压幅值的修正量。

第四步，第一轮有功迭代和无功迭代便做完了。

第五步,按公式计算平衡节点功率。直到节点不平衡功率下降到１０-5

以下,

迭代便可以结束。

第2章手工计算

第3章程序设计3.1 流程图

3.2 潮流计算程序

电力系统极坐标下的PＱ分解法潮流计算

ｄisp(＇电力系统极坐标下的PQ分解法潮流计算＇:);

cｌeaｒ

n=iｎｐut(＇请输入结点数：n=＇);

ｎ1=ｉｎｐut('请输入PＶ结点数：n１=');

n２=input('请输入PQ结点数：n２=');

isb=inpuｔ('请输入平衡结点:isb=');

pr=inpuｔ('请输入精确度:ｐr=');

K=inpｕｔ（'请输入变比矩阵看:K=＇);

C=inpuｔ('请输入支路阻抗矩阵：Ｃ=＇)；

ｙ=ｉnｐｕt(＇请输入支路导纳矩阵：y=＇)；

U=inｐut（'请输入结点电压矩阵:U=')；

S=inｐut（'请输入各结点的功率:S＝')；

Ｚ＝ｚerｏｓ(1，n);N=zeｒos(n２,n2+n1);L=zerｏs(n1+n２,n2)；ＱT1=ｚeroｓ（１,ｎ1＋n2);

for m＝1:n

for R=１：n

Ｃ（ｍ,m)=Ｃ(m,m)＋y(ｍ,R);

if K(m,Ｒ)~=０

C(m，m）=C(m,m)＋1/((K(m,R)*C(m,R）)/（K(m，R)－1）);

Ｃ（R,R)=C(Ｒ，Ｒ）+１／（(K(m,R)^2*C(m,R))/(１－K（m,R)));

C(m，R)=C(m，R)*K(m，R)；

C(R，ｍ）=Ｃ（m，R);

eｎd

end

end