文档介绍：课程设计任务书
题目
直角坐标表示的牛顿拉夫逊法潮流计算程序设计
学生姓名
张蔚
学号
200805080123
专业班级
电气0803
设
计
内
容
与
要
求
1. 设计要求
掌握MATLAB语言编程方法;理解和掌握运用计算机进行潮流计算的基本算法原理;针对某一具体电网,进行潮流计算程序设计。
其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解,培养学生运用所学知识分析和解决问题的能力。
2. 内容
1)学****并掌握MATLAB语言。
2)掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。
3)掌握电力系统功率方程、变量和节点分类。
4)掌握利用直角坐标表示的牛-拉法进行潮流计算的方法和步骤。
5)选择一个某一具体电网,编制程序流程框图。
6)利用MATLAB语言编写该模型的潮流计算程序,并上机调试程序,对计算结果进行分析。
7)整理课程设计论文。
起止时间
2011 年 6 月 21日至 2011 年6月 25日
指导教师签名
2011 年 6 月 21 日
系(教研室)主任签名
2011 年 6 月日
学生签名
2011 年 6 月日
已知一网络接线如下图所示, 试运用以直角坐标表示的牛顿—拉夫逊法计算该系统的潮流分布。实现下图所示的12节点电力系统潮流计算的程序设计。计算精度要求各节点电压修正量不大于。
牛拉法计算潮流的核心问题是修正方程的建立和求解。未说明这一修正方程的建立过程,首先对网络中各类节点的编号作如下约定:
⑴网络中共有n各节点,编号为1,2,…,n,其中包含一个平衡节点,编号为s;
⑵网络中有n-m个PV节点,编号为m+1,m+2,…,n。据此,我们可以建立如下的修正方程:
式中的,,分别为注入功率和节点电压的不平衡量。则有
式中的雅可比矩阵的各元素分别为
为求取这些偏导数,可将、、分别展开如下
时,由于对待定的j,只有特定节点的fj和ei是变量,则由上面两式可得:
j=i时,为使这些偏导的表示式更简洁,先引入节点注入电流的表示式如下
=[()+]+j[()+]
=a
由上式可得
四、程序流程图
N-L极坐标法:
启动
输入原始数据
形成节点导纳矩阵Y
设节点电压e,f;i=1,2,…,n,i≠s
置迭代次数 k=0
对PQ节点:按公式计算ΔP、ΔQ
对PV节点:按公式计算ΔP、ΔU
置节点号i=1
雅可比矩阵J是否形成,i>n?
计算雅可比矩阵元素J
增大节点号i=i+1
解修正方程式,由ΔP, ΔQ, ΔU, J求各节点电压变量
Δe,Δf;i=1,2,…,n,i≠s
求出、
迭代是否收敛,、≤ε?
计算平衡节点功率S和线路功率S
停止止
计算各节点电压新值:
e=e+,i=1,2,…,n,
f=f+;i≠s
增大迭代次数,k→k+1
五、程序
六、计算结果及分