姓名: 学号: 班级:
一、目的
根据所给的电力系统,编制短路电流计算程序,通过计算机进行调试,最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解,同时加强计算机实际应用能力的训练。 二、上机内容
电力系统故障的计算程序设计及编制和调试。采用所编制的程序进行《电力系统分析》例6-3题的对称短路计算。有关数学模型和原理框图以及已知结果的例题,参见《电力系统分析》第六章。常用的计算方法为节点导纳矩阵法或节点阻抗矩阵法,其形成方法分别参见《电力系统分析》第四章。
3、选择所用计算机语言的理由
我们使用的是第四代计算机语言的MATLAB,利用其丰富的函数资源,它的优点如下:
1. 语言简洁紧凑,使用方便灵活,库函数极其丰富。MATLAB程序书写形式自由,利用起丰富的库函数避开繁杂的子程序编程任务,压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写,用户不必担心函数的可靠性。可以说,用MATLAB进行科技开发是站在专家的肩膀上。
2.运算符丰富。由于MATLAB是用C语言编写的,MATLAB提供了和C语言几乎一样多的运算符,灵活使用MATLAB的运算符将使程序变得极为简短。
3.MATLAB既具有结构化的控制语句(如for循环,while循环,break语句和if语句),又有面向对象编程的特性。
4.程序限制不严格,程序设计自由度大。例如,在MATLAB里,用户无需对矩阵预定义就可使用。
5.程序的可移植性很好,基本上不做修改就可以在各种型号的计算机和操作系统上运行。
6.MATLAB的图形功能强大。在FORTRAN和C语言里,绘图都很不容易,但在MATLAB里,数据的可视化非常简单。MATLAB还具有较强的编辑图形界面的能力。
7.MATLAB的缺点是,它和其他高级程序相比,程序的执行速度较慢。由于MATLAB的程序不用编译等预处理,也不生成可执行文件,程序为解释执行,所以速度较慢。
8.功能强大的工具箱是MATLAB的另一特。MATLAB包含两个部分:核心部分和各种可选的工具箱。核心部分中有数百个核心内部函数。其工具箱又分为两类:功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能,图示建模仿真功能,文字处理功能以及与硬件实时交互功能。功能性工具箱用于多种学科。而学科性工具箱是专业性比较强的,如control,toolbox,signl proceessing toolbox,commumnication toolbox等。这些工具箱都是由该领域内学术水平很高的专家编写的,所以用户无需编写自己学科范围内的基础程序,而直接进行高,精,尖的研究。 9.源程序的开放性。开放性也许是MATLAB最受人们欢迎的特点。除内部函数以外,所有MATLAB的核心文件和工具箱文件都是可读可改的源文件,用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。
四、程序主框图、子框图及主要数据变量说明
流程图:
用公式(6-10)计算短路电流If If=1/(Zff+zf)z电流If
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
用公式(6-10)计算短路电流If If=1/(Zff+zf)z电流If
用公式(6-10)计算短路电流If If=1/(Zff+zf)z电流If
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
用公式(6-9)计算指定支路的电路 Ipq=(Kvp-Vq)/zpq
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)
输出结果
用公式(6-9)计算指定支路的电路 Ipq=(Kvp-Vq)/zpq
用公式(6-9)计算指定支路的电路 Ipq=(Kvp-Vq)/zpq
MATLAB流程图
五、关键程序及程序注释
Y=
-j13.8716 0 j9.5238 0 0
0 -j8.3333 0 j4.7619 0
j9.5238 0 -j15.2329 j2.2960 j3.4440
0 j4.7619 j2.2960 -j10.9646 j3.9360
程序如下:
jd=input('please input jd=’);
Y=[0-16.1905j, 9.5238j, 0, 0 , 0;
0+9.5238j, 37.4084j, 15.3846j, 12.5000j, 0;
0, 15.3846j, -35.3846j, 20.000j, 0;
0, 12.5000j, 20.000j, -37.9348j, 5.4348j;
0, 0, 0, 5.4348j, -9.9802j]
Z=inv(Y)
I=1/0.1860i
V1=1-0.0902i*I
V2=1-0.1533i*I
V3=0
V4=1-0.1611i*I
V5=1-0.0877i*I
I54=(V5-V4)/0.184i
I43=(V4-V3)/0.05i
I23=(V2-V3)/0.065i
I12=(V1-V2)/0.105i
I24=(V2-V4)/0.08i
运行结果如下:
please input jd=3
Y =
Columns 1 through 4
0 -16.1905i 0 + 9.5238i 0 0
0 + 9.5238i 0 +37.4084i 0 +15.3846i 0 +12.5000i
0 0 +15.3846i 0 -35.3846i 0 +20.0000i
0 0 +12.5000i 0 +20.0000i 0 -37.9348i
0 0 0 0 + 5.4348i
Column 5
0
0
0
0 + 5.4348i