电力系统分析课程设计—

1、2009-2010学年度下学期电力系统分析课程设计 电力系统短路故障的计算机 算法程序设计2010年 6 月 1 日目录第1章 课程设计目的.2第2章 课程设计内容.2第3章 课程设计原理.2第4章 选择语言理由.3第5章 程序主框图及主要数据变量说明.5第6章 程序及说明.9第7章 程序结果及分析.12第8章 参考文献.17第1章 课程设计目的根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学们对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。第2章 课程设计内容电力系统故障的计算程序设计及

2、编制和调试。第3章 课程设计原理1计算机计算原理应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三则是选择合适的计算机语言编制计算程序。建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂。运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用内存少，以利提高计算机的解题规模。选择合适的语言编写程序，就是首先确定用什么计

3、算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错误。所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。最后将所编制出的正确计算程序，用于电力系统的实际计算。2电力系统短路计算计算机算法一般在电力系统短路计算中，多数情况下只要计算短路电流、电压的周期分量起始值。因此，电力系统短路电流计算的数学模型，可归结为求解短路故障初始状态下的等值电路稳态解的问题。对于三相对

4、称短路，可建立一相等值电路的计算模型，对于不对称短路，则可应用对称分量法建立系统的正、负、零序网络，从而建立故障计算的序网模型。（1）对称短路计算有关数学模型和原理框图以及已知结果的例题，参见电力系统分析第六章。常用的计算方法为节点导纳矩阵法或节点阻抗矩阵法，其形成方法分别参见电力系统分析第四章。（2）简单不对称短路计算有关数学模型、原理框图以及已知结果的例题参见电力系统分析第八章8-4节。第4章 选择语言理由我们使用的是第四代计算机语言的MATLAB，利用其丰富的函数资源，它的优点如下： 1）语言简洁紧凑，使用方便灵活，库函数极其丰富。MATLAB程序书写形式自由，利用起丰富的库函数避开繁杂

5、的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。可以说，用MATLAB进行科技开发是站在专家的肩膀上。具有FORTRAN和C等高级语言知识的读者可能已经注意到，如果用FORTRAN或C语言去编写程序，尤其当涉及矩阵运算和画图时，编程会很麻烦。例如，如果用户想求解一个线性代数方程，就得编写一个程序块读入数据，然后再使用一种求解线性方程的算法（例如追赶法）编写一个程序块来求解方程，最后再输出计算结果。在求解过程中，最麻烦的要算第二部分。解线性方程的麻烦在于要对矩阵的元素作循环，选择稳定的算法以及代码的调试动不容易。即使有部分源代码，用户也会感

6、到麻烦，且不能保证运算的稳定性。解线性方程的程序用FORTRAN和C这样的高级语言编写，至少需要四百多行，调试这种几百行的计算程序可以说很困难。以下用MATLAB编写以上两个小程序的具体过程。MATLAB求解下列方程，并求解矩阵A的特征值。Ax=b,其中：A= 32 13 45 67 23 79 85 12 43 23 54 65 98 34 71 35b= 1 2 3 4解为：x=Ab;设A的特征值组成的向量e，e=eig（A）。可见，MATLAB的程序极其简短。更为难能可贵的是，MATLAB甚至具有一定的智能水平，比如上面的解方程，MATLAB会根据矩阵的特性选择方程的求解方法，所以用户根

7、本不用怀疑MATLAB的准确性。2）运算符丰富。由于MATLAB是用C语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短。3）MATLAB既具有结构化的控制语句（如for循环，while循环，break语句和if语句），又有面向对象编程的特性。4）程序限制不严格，程序设计自由度大。例如，在MATLAB里，用户无需对矩阵预定义就可使用。5）程序的可移植性很好，基本上不做修改就可以在各种型号的计算机和操作系统上运行。6）MATLAB的图形功能强大。在FORTRAN和C语言里，绘图都很不容易，但在MATLAB里，数据的可视化非常简单。MATLA

8、B还具有较强的编辑图形界面的能力。7）MATLAB的缺点是，它和其他高级程序相比，程序的执行速度较慢。由于MATLAB的程序不用编译等预处理，也不生成可执行文件，程序为解释执行，所以速度较慢。8）功能强大的工具箱是MATLAB的另一特色。MATLAB包含两个部分：核心部分和各种可选的工具箱。核心部分中有数百个核心内部函数。其工具箱又分为两类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能，图示建模仿真功能，文字处理功能以及与硬件实时交互功能。功能性工具箱用于多种学科。而学科性工具箱是专业性比较强的，如control,toolbox,signl proceessing too

9、lbox,commumnication toolbox等。这些工具箱都是由该领域内学术水平很高的专家编写的，所以用户无需编写自己学科范围内的基础程序，而直接进行高，精，尖的研究。9）源程序的开放性。开放性也许是MATLAB最受人们欢迎的特点。除内部函数以外，所有MATLAB的核心文件和工具箱文件都是可读可改的源文件，用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。第5章 程序主框图及主要数据变量说明1.电力系统节点方程的建立等值电路的制定发 电 机：电势源支路 电流源支路（含同步调相机）输电线路：PI型等值电路，R、X、B变 压 器： 一般负荷：恒定阻抗，电动机负荷：电势源支路 电流

10、源支路（同步电动机、感应电动机、以电动机为主的综合负荷，起始次暂态电流计算）2.利用节点阻抗矩阵计算短路电流基本原理3利用节点阻抗矩阵计算短路电流近似计算方法4利用节点阻抗矩阵计算短路电流计算流程图用例6-3的数据，线路的电阻和电容略去不计，电压器的标么变比等于参数：z12=J0.105, z45=J0.184,z24=J0.08,z23=J0.065,z34=J0.05Y11=-J16.1905,Y12=J905238,Y21=J905238,Y22-J37.4084,Y23=J15.3846,Y24=J12.500,Y32=J15.3846,Y33=-J35.3846,Y34=J20.00

11、0,Y42=J15.5000,Y43=-J37.9348,Y45=J5.4348,Y54=J5.4348,Y55=-J9.9802 电路图：流程图： 输入数据形成节点导纳矩阵选择故障点If输入数据计算节点阻抗矩阵If列元素用公式（6-10）计算短路电流If If=1/（Zff+zf）z电流If用公式（6-11）计算各点电压 Vi=1-Zif/( Zff+zf)用公式（6-9）计算指定支路的电路 Ipq=(Kvp-Vq)/zpq输出结果Y=0 -16.1905i 0 + 9.5238i 0 0 0 0 + 9.5238i 0 -37.4084i 0 +15.3846i 0 +12.5000i 0

12、 0 0 +15.3846i 0 -35.3846i 0 +20.0000i 0 0 0 +12.5000i 0 +20.0000i 0 -37.9348i 0 + 5.4348i 0 0 0 0 + 5.4348i 0 - 9.9802i 采用4-3节的方法Z=0 + 0.1181i 0 + 0.0958i 0 + 0.0902i 0 + 0.0858i 0 + 0.0467i 0 + 0.0958i 0 + 0.1629i 0 + 0.1533i 0 + 0.1459i 0 + 0.0794i 0 + 0.0902i 0 + 0.1533i 0 + 0.1860i 0 + 0.1611i

13、0 + 0.0877i 0 + 0.0858i 0 + 0.1459i 0 + 0.1611i 0 + 0.1729i 0 + 0.0941i 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0.0941i 0 + 0.1515i 计算电流及网络中的电流分布If=Vi/Zii=1/Zii V1=V2=V3=V4=V5=1 (注意：Vi=0)V1=V1-Z1i*IfV2=V2-Z2i*IfV3=V3-Z3i*IfV4=V4-Z4i*IfV5=V5-Z5i*IfI54=(V5-V4)/z45I43=(V4-V3)/z43I23=(V2-V3)/z23I12=(V1-

14、V2)/z12I24=(V2-V4)/z42第6章 程序及说明 程序如下：jd=input('please input jd=');if jd=1Y=-16.1905j, 9.5238j, 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.9802jinv(Y)I=1/0.1181iV1=0V2=1-0.0958i*IV3=1-0.090

15、2i*IV4=1-0.0858i*IV5=1-0.0467i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=2Y=-16.1905j, 9.5238j, 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.

16、9802jinv(Y)I=1/0.1629iV1=1-0.0958i*IV2=0V3=1-0.1533i*IV4=1-0.1459i*IV5=1-0.0794i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=3Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.500

17、0j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1860iV1=1-0.0902i*IV2=1-0.1533i*IV3=0V4=1-0.1611i*IV5=1-0.0877i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=4Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j

18、, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1729iV1=1-0.0858i*IV2=1-0.1459i*IV3=1-0.1611i*IV4=0V5=1-0.0941i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=5Y=-16.1905j, 9.5238

19、j, 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.9802jinv(Y)I=1/0.1515iV1=1-0.0467i*IV2=1-0.0794i*IV3=1-0.0877i*IV4=1-0.0941i*IV5=0I54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.1

20、05iI24=(V2-V4)/0.08i else c=errorEnd第7章 程序结果及分析节点1：please input jd=1 y11 =0 -16.1905i y12 =0 + 9.5238i y13 = 0y14 = 0y15 =0y21 =0 + 9.5238iy22 =0 -37.4084iy23 =0 +15.3846iy24 =0 +12.5000iy25 =0y31 =0y32 =0 +15.3846iy33 =0 -35.3846iy34 =0 +20.0000iy35 =0y41 =0y42 =0 +12.5000iy43 =0 +20.0000iy44 =0 -3

21、7.9348iy45 =0 + 5.4348iy51 =0y52 =0y53 =0y54 =0 + 5.4348iy55 =0 - 9.9802iy0 = Columns 1 through5 0 -16.1905i 0 + 9.5238i 0 0 0 0 + 9.5238i 0 -37.4084i 0 +15.3846i 0 +12.5000i 0 0 0 +15.3846i 0 -35.3846i 0 +20.0000i 0 0 0 +12.5000i 0 +20.0000i 0 -37.9348i 0 + 5.4348i 0 0 0 0 + 5.4348i 0 - 9.9802ians

22、=Columns 1 through 4 0 + 0.1181i 0 + 0.0958i 0 + 0.0902i 0 + 0.0858i 0 + 0.0958i 0 + 0.1629i 0 + 0.1533i 0 + 0.1459i 0 + 0.0902i 0 + 0.1533i 0 + 0.1860i 0 + 0.1611i 0 + 0.0858i 0 + 0.1459i 0 + 0.1611i 0 + 0.1729i 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0.0941iColumn 5 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0.0941i 0 + 0.1515ii =0 - 8.4674iv1 =0v2 =0.1888v3 =0.2362v4 =0.2735v5 =0.6046i54 =0 - 1.7993ii43 =0 - 0.7451ii23 =0 + 0.7295ii