电力系统分析课设 C语言.doc

《电力系统分析课设 C语言.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、电力系统分析课程设计目录1课程设计说明12编程语言的选择及理由13程序设计23.1主体流程图23.2详细流程图33.2.1创建系统（create）流程图33.2.2加载系统函数（load）流程图33.2.3计算子函数（calculate）流程图43.2.4改变短路节点（change）流程图43.3数据及变量说明53.4程序代码及注释63.5测试算例64设计体会13参考文献13附录13i1课程设计说明本课程设计根据电力系统分析课程中所讲的基于节点方程的三相短路计算的原理和方法。采用c语言编程。实现了三相短路计算的计算机实现。本程序输入含有网络拓扑信息和电器元件的参数，输出三相短路的短

2、路电流和短路后网络各节点的电压参数和各支路的电流参数。并以文件的形式保存，方便用户的读取，打印。实现了设计的要求。三相短路在电力系统中虽然发生的概率很小，但一旦发生产生的影响非常严重，电力系统中设备的动稳定和热稳定校验都要以三相短路电流为依据，因此我选择做是电力系统三相短路的计算程序，觉得更有意义。应用计算机进行电力系统计算，首先要建立电力系统相应的数学模型，如建立电力系统等效网络和节点方程；其次是运用合理的计算方法计算结果，如三角分解法，LDU分解法等；第三则是选择合适的计算机语言编制计算程序，实现程序的准确性，提高精确度。本设计利用节点方程来实现电力系统三相短路计算。首先，根据

3、给定的电力系统运行方式制定系统的一相等值网络，计算出各元器件的参数，其次根据网络方程利用变压器和线路参数形成不含发电机和负荷的节点导纳矩阵Y；然后，利用系统形成的节点导纳矩阵求解系统阻抗矩阵Z；第三则是应用节点阻抗矩阵计算短路电流。电力系统短路电流的实际工程计算中，许多实际问题的解决并不需要十分精确的结果。因此采用近似计算法，在建立系统节点方程时，采用网络模型和参数做简化处理，等值电路的制定通常将发电机作为含源支路，表示为节点注入电流源和节点电抗的形式；系统中存在同步发电机时，按发电机支路处理；同时，忽略发电机、变压器和输电线路的电阻，不计线路的电容，略去变压器的励磁电流，负荷忽略

4、不计；在标幺参数计算时，选取各级平均额定电压作为基准电压，忽略元件的额定电压和平均电压的差别，认为变压器变比的标幺值为1；此外，还假定所有的发电机的电势具有相同的相位，所有的元件仅用电抗表示，避免复数运算。三相对称短路为金属性短路。本课程设计的关键部分是利用系统节点导纳矩阵形成节点阻抗矩阵，本程序采用LDU分解法，然后利用节点导纳矩阵和阻抗矩阵的乘积为单位矩阵列写方程，节点导纳矩阵经LDU分解后，回代求解系统阻抗矩阵。利用节点阻抗矩阵进行短路计算，首先利用故障前的节点注入电流和系统节点阻抗矩阵求取故障前系统的各节点电压；然后，利用金属短路时短路点残压为0,对地电阻为0，计算短路点电

5、流；最后，利用故障前节点电压、节点阻抗和短路点电流求得短路后各节点电压和各支路电流。2编程语言的选择及理由能够实现程序编写的语言有多种，本课程设计中采用C语言进行电力系统对称短路计算程序的编写。个人对c语言和metlab相对较熟悉。之所以选择C语言编程，从主观上说，是因为我学过C语言课程，相对于其他的编程语言而言，C语言的应用更为熟悉和了解。另外，从客观上说，C语言是目前世界上最为流行、使用最为广泛的高级程序设计语言。，C的效率MATLAB比高多了，MATLAB虽然语法简单，但执行效率不高。MATLAB是用于特定的方面的，比如说矩阵运算方面、DSP就很有优势，而C语言不同。很多软件

6、的底层都是C编写的。MATLAB提供的是一种基于解释的语言，虽然也是一种高级语言，但是还是很简单的，它的目标在于用户方便，开发简单。而实际上这些方便也是有代价的，编译器承受了这些代价，程序执行效率方面也有损失。另一方面MATLAB这么做也是基于其专用于科学计算的目的来考虑的，无法应用于通用的程序设计。c语言这么做保持了其编译的速度，便于大型程序设计，执行效率也很高，因此我选择用c语言编写程序。1程序设计1.1主体流程图开始欢迎使用用户选择操作，按键123进入创建系统子函数（create）进入加载系统子函数（load）进入退出子函数进入计算子函数（calculate）改短路点子函数（

7、change）是否改路点？结束是图3-1主体流程图1.1详细流程图1.1.1创建系统（create）流程图输入文件名输入电抗支路参数结束输入节点数、支路数、故障节点输入发电机支路参数输入电容支路参数图3-2create函数流程图1.1.1加载系统函数（load）流程图输入系统文件名输出查看系统参数选择进入下一个子函数图3-3节点导纳矩阵形成流程图1.1.2计算子函数（calculate）流程图文件读入电路参数求取LDU分解的因子表矩阵下三角求取DUZj向量求解完毕？否