《电力系统潮流计算课程设计—牛顿-拉夫逊法的应用课程设计》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计—牛顿-拉夫逊法的应用课程设计（26页珍藏版）》请在金锄头文库上搜索。
1、 电力系统潮流计算课程设计电力系统潮流计算课程设计牛顿-拉夫逊法的应用学院学院: :电气工程学院电气工程学院班级班级: :电电 125125学号学号:1212002047:1212002047姓名姓名: :王衡佳王衡佳教师教师: :张新松张新松日期日期:2015.11.15:2015.11.15南通大学 目目 录录一、课程设计内容简介一、课程设计内容简介.1.11.1.课程设计课题课程设计课题.1.12.2.电力系统潮流计算电力系统潮流计算.1.1二、电力系统潮流计算方法二、电力系统潮流计算方法.1.11.1.高斯高斯- -赛德尔迭代法赛德尔迭代法 . .1.12.2.牛顿牛顿- -拉夫逊法拉
2、夫逊法.1.13.P-Q3.P-Q 分解法分解法.2.2三、程序设计三、程序设计.3.31.1.程序设计环境程序设计环境.3.32.2.主程序清单主程序清单.3.33.3.程序运行结果程序运行结果.22.22四、个人小结四、个人小结.23.231.1.课程设计问答课程设计问答.23.232.2.个人心得体会个人心得体会.23.23五、参考文献五、参考文献.24.24南通大学 0一、课程设计内容简介一、课程设计内容简介1.1.课程设计课题课程设计课题设计一个基于高斯赛德尔法或牛顿拉夫逊法的电力系统潮流计算程序，要 求计算 IEEE14 节点，可采用 VC+，MATLAB 或其他工具实现，了解节点
3、导纳矩 阵内各元素的意义以及 PV 节点与 PQ 节点的区别。2.2.电力系统潮流计算电力系统潮流计算电力系统潮流计算是研究电力系统稳态运行运行情况的一种基本电气计算。 它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母 线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮 流计算的结果是电力系统稳定计算和故障分析的基础。 二、电力系统潮流计算方法二、电力系统潮流计算方法1.1.高斯高斯- -赛德尔迭代法赛德尔迭代法以导纳矩阵为基础，并应用高斯-塞德尔迭代的算法是在电力系统中最早 得到应用的潮流计算方法，目前高斯一塞德尔法已很少使用。 将所求方程 f(x) =
4、0 改写为 x=(x)，不能直接得出方程的根，给一个猜测值 x0 得：x1=(x0),又 可取 x1 为猜测值，进一步得：优点： (1)原理简单，程序设计十分容易。 (2)导纳矩阵是一个对称且高度稀疏的矩阵，因此占用内存非常节省。 (3)就每次迭代所需的计算量而言，是各种潮流算法中最小的，并且和网络所包 含的节点数成正比关系。 缺点： (1)收敛速度很慢。 (2)对病态条件系统，计算往往会发生收敛困难：如节点间相位角差很大的重负 荷系统、包含有负电抗支路(如某些三绕组变压器或线路串联电容等)的系统、 具有较长的辐射形线路的系统、长线路与短线路接在同一节点上，而且长短线 路的长度比值又很大的系统
5、。 (3)平衡节点所在位置的不同选择，也会影响到收敛性能。2.2.牛顿牛顿- -拉夫逊法拉夫逊法牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。自从 20 世纪 60 年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计 算速度方面都超过了其他方法，成为直到目前仍被广泛采用的方法。南通大学 1优点： (1)收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭 代 45 次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规 模基本无关。 (2)具有良好的收敛可靠性，对于前面提到的对以节点导纳矩阵为基础的高斯一 塞德尔法呈病态的系统，牛顿法均能可靠地收敛
6、。 (3)牛顿法所需的内存量及每次迭代所需时间均较前述的高斯一塞德尔法为多， 并与程序设计技巧有密切关系。 缺点： 牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有 可能根本不收敛或收敛到一个无法运行的解点上。 解决方法： 对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且 各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为 “平直电压”)，“平直电压”法假定： 这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质 量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问 题。可以先用高斯一塞德尔法迭代 1-2 次
7、；以此迭代结果作为牛顿法的初值， 也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法 迭代。3.3.P-QP-Q 分解法分解法电力系统中常用的 P-Q 分解法派生于以极坐标表示的牛顿-拉夫逊法,其基 本思想是把节点功率表示为电压向量的极坐标形式，以有功功率误差作为修正 电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功和无 功分开进行迭代其主要特点是以一个（n-1）阶和一个 m 阶不变的、对称的系数 矩阵 B,B 代替原来的（n+m-1）阶变化的、不对称的系数矩阵 M，以此提高计算南通大学 2速度，降低对计算机贮存容量的要求。P-Q 分解法在计算速度方面有显著的
8、提 高，迅速得到了推广。本课程设计中不采用此方法，因此不再过多赘述。 3 3、程序设计程序设计1.1.程序设计环境程序设计环境(1)采用方法:牛顿-拉夫逊法 (2)设计工具：VC+6.0 程序 (3)运行环境：win8 64 位操作系统2.2.主程序清单主程序清单#include #include #include #include #include #include #include #include const double CalculateError=0.00001;/定义计算收敛条件 static int AllNodeNumber;/节点总数 static int BranchNu
9、mber;/支路总数 static int GroundNumber;/接地支路总数class BUS public: int OldNumber;/旧的节点编号 int NodeType;/节点类型(0pq 节点,2pv 节点,3平衡节点) double LoadP;/负荷有功 double LoadQ;/负荷无功 double GenP;/发电机有功 double GenQ;/发电机无功 double VoltageVal;/电压幅值 double VoltageAngle;/电压相角 double NodeP;/节点净有功 double NodeQ;/节点净无功 double NodeE
10、;/电压 e 分量 double NodeF;/电压 f 分量BUS():OldNumber(0),NodeType(0),LoadP(0),LoadQ(0),GenP(0),GenQ(0),VoltageVal (1.0),VoltageAngle(0.0),NodeP(0.0),NodeQ(0.0),NodeE(1.0),NodeF(0.0) ; class BRANCH 南通大学 3public: int NodeI;/节点 i 侧 int NodeJ;/节点 j 侧 int Circuit;/判断是否是双回路(0 单回路，1 双回路) double R;/支路电阻 double X;/
11、支路电抗 double Y;/线路变压器变比或接地导纳 BBRANCH():NodeI(0),NodeJ(0),Circuit(1),R(0.0),X(0.0),Y(0.0) ; class Yii/存导纳矩阵对角元素 public: double G,B; Yii():G(0.0),B(0.0) ; class Yij/存导纳矩阵非对角元素和列号 public: double G,B; int j; Yij():G(0),B(0),j(0) ; class GROUND/存接地支路数据 public:int Node; double G; double B; ; class CURRENT/存支路电流 public: double E,F; ; class POWER/存支路功率 public: double P1,Q1; double P2,Q2; ; class POWERFLOW 南通大学 4public: double maxerror; char *FileName, *ResultName; void InputData();/ 读入系统数据 void Tinny2();/ 半动态节点优化编号 void FormY();/
