课程设计任务书
学生姓名：专业班级：
指导教师：工作单位：
题目: 简单闭式网络潮流计算程序设计
初始条件：
节点数：4 支路数：4 计算精度：0.00010
支路1： 3 (0.0200+j0.0800)
支路2： 4 （0.0400+j0.1200)
支路3： 4 （0.0500+j0.1400）
支路4： 4 （0.0400+j0.1200)
节点1：PQ节点，S(1)=-0.6000-j0.2500
节点2：PQ节点，S(2)=-0.8000-j0.3500
节点3：PV节点，P(3)=0.4000 V(3)=0.9500
节点4：平衡节点，U(4)=1.0000 0.0000
要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）
（1）编写计算原理。
（2）编写计算程序。
（3）编写程序并计算结果。
（4）结果分析。
时间安排：
第16周，6月6—10号，指导老师下达任务书，审题、查阅相关资料；
第16周，6月11—16号，分析、计算，编写程序，撰写报告；
第17周，6月17号，论文答辩。
指导教师签名：年月日
系主任（或责任教师）签名：年月日
本文运用matlab软件对给定课设题目进行潮流计算。潮流计算是电力系统课程最基本最常用的计算。潮流计算是指对电力系统正常运行状况的分析和计算。通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、PQ分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。
通过潮流计算，可以确定各母线的电压幅值和相角，各元件流过的功率，整个系统的功率损耗。潮流计算是实现安全经济发供电的必要手段和重要工作环节。因此潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有广泛的运用。
本课程设计采用PQ分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。
关键词：PQ分解法matlab 潮流计算
1.题目原始数据及化简 (1)
2.PQ分解法 (2)
2.1PQ分解法基本思想 (2)
2.2 PQ分解法潮流计算基本步骤 (5)
3编程及运行 (6)
3.1 PQ分解法潮流计算程序框图 (6)
3.2源程序代码 (7)
3.3运行程序及结果分析： (19)
4.小结 (22)
5.参考文献 (23)
1
1.题目原始数据及化简
原始数据：
编程计算潮流
节点数：4 支路数：4 计算精度：0.00010
支路1：
3 (0.0200+j0..0800)
支路2
： 4 （0.0400+j0.1200) 支路3：
4 （0.0500+j0.1400） 支路4： 4 （0.0400+j0.1200)
节点1：PQ 节点，S(1)=-0.6000-j0.2500
节点2：PQ 节点，S(2)=-0.8000-j0.3500
节点3：PV 节点，P(3)=0.4000 V(3)=0.9500
节点4：平衡节点，U(4)=1.0000 0.0000
根据原始数据所画电路简化图如图1 1 3
4
2
图1电路简化图
2 2.PQ 分解法
2.1PQ 分解法基本思想
PQ 分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标形式，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功和无功分开进行迭代其主要特点是以一个（n-1）阶和一个m 阶不变的、对称的系数矩阵B ,B '''代替原来的（n+m-1）阶变化的、不对称的系数矩阵M ，以此提高计算速度，降低对计算机贮存容量的要求。
牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，修正方程式为：
⎥⎦⎤⎢⎣⎡∆∆⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡∆∆V V L N J H Q P /δ
或展开为：
V V L J Q V V N H P //∆⋅+∆⋅=∆∆⋅+∆⋅=∆δδ （1）
以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。
电力系统中有功功率主要与各节点电压向量的角度有关，无功功率则主要受各节点电压幅值的影响。大量运算经验也告诉我们，矩阵N 及J 中各元素的数值相对是很小的，因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代，即将式(4)化简为：
V V L Q H P /∆⋅=∆∆⋅=∆δ （2）
这样，由于我们把2n 阶的线性方程组变成了二个n 阶的线性方程组，因而计算量和内存方面都有改善。但是，H ，L 在迭代过程中仍然不断变化，而且又都是不对称矩阵。对牛顿法的第二个化简，也是比较关键的一个化简，即把式
(5)中的系数矩阵简化为在迭代过程中不变的对称矩阵。
众所周知，一般线路两端电压的相角差是不大的(通常不超过10~20度)，因此可以认为：
cos 1sin ij ij ij ij G B θθ≈ （3）
