昆明理工大学 《电力系统计算机辅助分析》上机实验指导书
实验二 电力系统计算中常用的数值算法及电力网络的数
学模型
上机实验目的
1、掌握高斯消去法、因子表法求解线性方程组的基本原理； 4、读懂课本中相应算法的实现例程；
5、 掌握节点导纳矩阵和节点阻抗矩阵的生成算法； 6、 掌握电力系统计算所需原始数据的准备过程； 7、 了解节点优化编号的基本原理和算法。
预习要求
1、 回顾《计算方法》中的基本理论，准备上述算法的计算流程图； 2、 阅读课本中的相应程序，将各算法改编为独立的函数文件。 3、 读懂课本中的算法实现代码，将其改编为独立的函数文件； 4、 整理课本例17.6的原始数据，设计成独立的数据文件；
5、 查阅资料了解某一商用或开源电力系统仿真软件对原始数据文件的要求。
上机实验内容
1、 用因子表法求解电路方程中的节点电压；
u1 I1 10V u210Ω 50Ω 30Ω I2 15Ω u3I3 30Ω
2、将图2.1中的电压源分别改为2V、4V、6V、8V时，用因子表进行求解，体
会该算法带来的益处；
5
昆明理工大学 《电力系统计算机辅助分析》上机实验指导书
3、 上机验证课本例17.6－17.7的计算结果；
4、 设计系统变更时修改导纳矩阵的函数（参考课本P.323），并调试通过。
思考题
1、高斯消去法与因子表法有何异同？ 2、什么是雅可比矩阵？
3、节点导纳阵与节点阻抗阵之间有何关系，二者之间如何转换？
4、 网络矩阵的修改函数对电力系统计算有何作用？
报告要求及其它说明
1、 按实验报告格式进行提交；
2、 提交各个自行封装的算法函数的文件清单、源码和调用规则； 3、 说明原始数据文件的输入要求；
4、 对验证的程序，报告必须包括应用过程和结果； 5、 给出导纳矩阵修正函数的调用规则及其实现代码； 。
6
昆明理工大学 《电力系统计算机辅助分析》上机实验指导书
实验三 短路电流计算
上机实验目的
1、掌握电力系统短路计算的基本原理及算法； 2、掌握计算机求取转移电抗的算法； 3、灵活掌握阻抗矩阵在短路计算中的应用。
预习要求
1、看懂课本的例程，将其划分为相对独立的模块；
2、设计一个新的计算流程，合并课本第18章例程的计算功能； 3、整理故障计算必需的已知条件。
上机实验内容
1、 合并课本第18章例程。即要求实现计算指定节点各种故障条件下，各支路的
故障电流分布情况。其中对于实用算法，可只提供计算电抗。提示：尽量使用先前实验中已封装好的函数。
2、 用新程序验证第18章的各例题。
思考题
1、计算短路电流的目的是什么？
2、在计算短路电流时，是否用到了节点阻抗阵中的所有元素？据此，可以如何简化短路计算程序？
报告要求及其它说明
1、按实验报告格式进行提交；
2、提供合并后计算程序的流程图与源代码； 3撰写合并后计算程序的使用说明；
3、结合思考题完成实验报告。鼓励在报告的讨论部分说明进一步丰富该计算程序的设想。
7
昆明理工大学 《电力系统计算机辅助分析》上机实验指导书
计算电抗算短路电流例18.3，课本P358页程序中用到的短路电流计算曲线数据
Xjs 0s 0.01s 3.897 3.310 0.06s 3.494 3.063 0.1s 3.38 2.902 0.2s 3.022 2.646 0.4s 2.705 2.410 0.5s 2.634 2.360 0.6s 2.589 2.332 1s 2.507 2.299 2s 2.415 2.306 4s 2.378 2.316 0.2731 4.025 0.3227 3.368
实验四 潮流计算
上机实验目的
1、 掌握牛顿拉夫逊法求解非线性方程组的基本原理； 2、 了解欧拉法及改进欧拉法求解微分方程的基本原理； 3、 掌握用N-L法计算潮流的基本原理和程序设计要点； 4、 掌握用PQ分解法计算潮流的基本原理和程序设计要点。
预习要求
1、准备潮流计算前各种已知参数表格；
2、找出课本中两个潮流例程中实现相同功能的代码，设计合理的函数将其封装起来。
上机实验内容
1、用牛顿拉夫逊法求解如下非线性方程组，注意保留迭代过程中的中间值。
22??3??x1?x2?5x1?1?0(0) x???设置精度值，0.01 0.0001，最大迭代次数14 15?223????x1?x2?3x2?3?0次
2、上机验证课本例19.1和例19.2的计算结果；
3、将两种算法的程序计算程序进行合并，要求共用输入、输出部分，网络方程形成函数和线性方程组求解函数，并实现算法的人工选择；
4、根据实际的进度（可选择）进一步丰富潮流计算程序的功能，如网络修改、节点电压越限提示、支路过载Q越线提示等功能。
8
昆明理工大学 《电力系统计算机辅助分析》上机实验指导书
思考题
1、两种算法的收敛过程图中出现差异的原因是什么？出现降阶， 2、PQ分解法是如何利用电力系统特征的？
3、为什么PQ分解法在经过很大的简化后却不会降低计算结果的精度？ 4、潮流计算中的平衡节点起什么作用？
5、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？
答：因为潮流计算节点注入的有功和无功是 如果电压复制越线，系统武功太充足，PQ变为pv节点，把V座位上限制计算；如果是PV节点，无功越线，超过上限，电压幅值，pv变pq，以上线作为武功定值，从新确定PQ结点个数再进行计算。
6、如果要以某一潮流分布情况作为短路计算的初值，如何考将潮流计算程序与短路计算程序联系起来？
报告要求及其它说明
1、按实验报告格式进行提交；
2、提供合并后计算程序的流程图与源代码； 3、撰写合并后计算程序的使用说明；
4、结合思考题完成实验报告。鼓励在报告的讨论部分说明进一步丰富该计算程序的设想。
9