《电力系统计算机辅助分析》上机实验指导书4Word下载.docx

1、1、利用图书馆和互联网资源，了解MATLAB软件的基本信息。列举出MATLAB软件的主要特点和应用领域。2、列出实现三种结构的和程序流程控制的关键字和语法规则；3、学习MATLAB函数的编写方法。上机实验内容1、 启动MATLAB软件，识别出常用的四个窗口命令窗口（Command Window）、历史命令窗口（Command History）、工作空间浏览器（Workspace Browser）。2、 在用户盘上以自己的学号建立一个文件夹，并将该文件夹设置为当前工作路径。在命令窗口中可以执行与Dos命令一样的文件操作。3、 在命令窗口中执行“helpdesk”命令，仔细阅读MATLAB的联机帮

2、助，了解该帮助系统的使用方法。4、上机求解下面的电路问题，脚本文件已经给出。已知某个电路的端电压的表达式为：区间时，试绘制电压的波形。提示：在MATLAB的M文件编辑器中键入以下命令语句，并保存为ex.m：Clear,clc,closet=0:0.5:10;u=6*exp（-2*t）;plot（t,u）;title（RC电路的电压响应曲线）;xlabel（时间/s）;ylabel（电压/v）grid%5、 对于如下电路图，已知V，自己编写程序，画出和的波形（四个周期）。写出的向量表达式用于的计算（写出串联段和并联段的阻抗，利用分压公式），写出的幅值和相角就能写出波形表达式，设定波形的长度（四个

3、周波）和采样频率（两点间隔）后作图。注意画图时使用标题，横纵坐标及网格线等命令。%mUsm=200; %的幅值aUsm=（pi/180）*（+90）;的相角Usm=mUsm*exp（j*aUsm）;用于计算mUcm=abs（Ucm）; %求幅值aUcm=angle（Ucm）; %求相角T=4*2*pi/w; %用w表示，T表示四个周波N=100; %画图的点数dt=T/N; %两点间隔，即采样频率dt:T;for k=1:101 %画波形us（k）=mUsm*cos（w*t（k）+aUsm）;uc（k）=mUcm*cos（w*t（k）+aUcm）;endplot（） %注意区分两条曲线的线型和

4、颜色grid ontitle（）xlabel（,fontsize,16） % fontsize,16表示横坐标的字体大小为16ylabel（）legend（us,uc）%思考题及说明（选作3题）1、 MATLAB软件是什么？有什么特点？2、 怎样在命令窗口中进行简单的计算？3、 如何获得命令或函数的帮助文档？4、 MATLAB的数据类型与其它高级语言相比有何特点？5、 什么时候用for循环、什么时候用while循环？6、 MATLAB函数必须遵循的规则有哪些？报告要求及其它说明1、按实验报告格式进行提交；2、对验证的程序，报告必须包括应用过程和结果；3、对编制的程序，报告必须包括文件清单、流程

5、图、源码、计算实例及结果。4、结合思考题完成实验报告。鼓励在报告的讨论部分提出对例程的修改方案。实验二电力系统计算中常用的数值算法及电力网络的数学模型1、掌握高斯消去法、因子表法求解线性方程组的基本原理；4、读懂课本中相应算法的实现例程；5、 掌握节点导纳矩阵和节点阻抗矩阵的生成算法；6、 掌握电力系统计算所需原始数据的准备过程；7、 了解节点优化编号的基本原理和算法。1、 回顾计算方法中的基本理论，准备上述算法的计算流程图；2、 阅读课本中的相应程序，将各算法改编为独立的函数文件。3、 读懂课本中的算法实现代码，将其改编为独立的函数文件；4、 整理课本例17.6的原始数据，设计成独立的数据文

6、件；5、 查阅资料了解某一商用或开源电力系统仿真软件对原始数据文件的要求。1、 用因子表法求解电路方程中的节点电压；2、将图2.1中的电压源分别改为2V、4V、6V、8V时，用因子表进行求解，体会该算法带来的益处；3、 上机验证课本例17.617.7的计算结果；4、 设计系统变更时修改导纳矩阵的函数（参考课本P.323），并调试通过。思考题1、高斯消去法与因子表法有何异同？2、什么是雅可比矩阵？3、节点导纳阵与节点阻抗阵之间有何关系，二者之间如何转换？4、 网络矩阵的修改函数对电力系统计算有何作用？1、 按实验报告格式进行提交；2、 提交各个自行封装的算法函数的文件清单、源码和调用规则；3、

7、说明原始数据文件的输入要求；4、 对验证的程序，报告必须包括应用过程和结果；5、 给出导纳矩阵修正函数的调用规则及其实现代码；。实验三短路电流计算1、 掌握电力系统短路计算的基本原理及算法；2、 掌握计算机求取转移电抗的算法；3、 灵活掌握阻抗矩阵在短路计算中的应用。1、 看懂课本的例程，将其划分为相对独立的模块；2、 设计一个新的计算流程，合并课本第18章例程的计算功能；3、 整理故障计算必需的已知条件。1、 合并课本第18章例程。即要求实现计算指定节点各种故障条件下，各支路的故障电流分布情况。其中对于实用算法，可只提供计算电抗。尽量使用先前实验中已封装好的函数。2、 用新程序验证第18章的

8、各例题。1、 计算短路电流的目的是什么？2、 在计算短路电流时，是否用到了节点阻抗阵中的所有元素？据此，可以如何简化短路计算程序？2、提供合并后计算程序的流程图与源代码；3撰写合并后计算程序的使用说明；3、 结合思考题完成实验报告。鼓励在报告的讨论部分说明进一步丰富该计算程序的设想。计算电抗算短路电流例18.3，课本P358页程序中用到的短路电流计算曲线数据Xjs0s0.01s0.06s0.1s0.2s0.4s0.5s0.6s1s2s4s0.27314.0253.8973.4943.383.0222.7052.6342.5892.5072.4152.3780.32273.3683.3103.0

9、632.9022.6462.4102.3602.3322.2992.3062.316实验四潮流计算1、 掌握牛顿拉夫逊法求解非线性方程组的基本原理；2、 了解欧拉法及改进欧拉法求解微分方程的基本原理；3、 掌握用N-L法计算潮流的基本原理和程序设计要点；4、 掌握用PQ分解法计算潮流的基本原理和程序设计要点。1、 准备潮流计算前各种已知参数表格；2、 找出课本中两个潮流例程中实现相同功能的代码，设计合理的函数将其封装起来。1、 用牛顿拉夫逊法求解如下非线性方程组，注意保留迭代过程中的中间值。设置精度值，0.01 0.0001，最大迭代次数14 15次2、 上机验证课本例19.1和例19.2的计

10、算结果；3、 将两种算法的程序计算程序进行合并，要求共用输入、输出部分，网络方程形成函数和线性方程组求解函数，并实现算法的人工选择；4、 根据实际的进度（可选择）进一步丰富潮流计算程序的功能，如网络修改、节点电压越限提示、支路过载Q越线提示等功能。1、 两种算法的收敛过程图中出现差异的原因是什么？出现降阶，2、 PQ分解法是如何利用电力系统特征的？3、 为什么PQ分解法在经过很大的简化后却不会降低计算结果的精度？4、 潮流计算中的平衡节点起什么作用？5、 为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的

11、？哪些可调？哪些不可调？答：因为潮流计算节点注入的有功和无功是 如果电压复制越线，系统武功太充足，PQ变为pv节点，把V座位上限制计算；如果是PV节点，无功越线，超过上限，电压幅值，pv变pq，以上线作为武功定值，从新确定PQ结点个数再进行计算。6、 如果要以某一潮流分布情况作为短路计算的初值，如何考将潮流计算程序与短路计算程序联系起来？3、撰写合并后计算程序的使用说明；实验五静态稳定和暂态稳定计算1、 掌握简单电力系统稳定计算的基本原理；2、 掌握分段计算法的实现过程；3、 部分了解MATLAB矩阵分析、数值积分算法应用。1、 看懂课本的例程；2、 用手工推导例20.1与例20.2的计算过程

12、。1、 上机验证例20.1与例20.2的计算结果；2、 用MATLAB数值工具箱中的函数（解微分方程、求特征值）替代例程中的相应部分，重复验证过程。1、 简单电力系统的稳定性计算方法可能应用于什么样的实际系统中？2、 MATLAB提供的算法有何特点？2、提交计算原理的详细说明；3、提交修改后的源程序；鼓励在报告的讨论部分进一步探讨更如何在单机无穷大系统中使用更详细的发电机及其控制系统的模型，或经典模型下的多机电力系统稳定性分析方法。 学生实验报告表头格式昆明理工大学 电力工程 学院学生实验报告实验课程名称： 电力系统计算机辅助分析 开课实验室： 年 月 日年级、专业、班学号姓名成绩实验项目 名 称电力系统计算中常用的数值算法指 导 教 师教 师 评 语 教师签名：年 月 日一、上机实验目的二、算法功能与原理三、程序流程四、过程记录五、源码及仿真结果六、讨论