课程设计---电力系统短路故障的计算机算法程序设计.doc

信息工程学院课程设计（论文）电力系统分析课程设计报告书 题目: 电力系统短路故障的计算机算法程序设计 专 业：电气工程及其自动化班 级： 学 号： 学生姓名： 指导教师： 2012年 3 月 11 日 学生姓名 学 号 成 绩设计题目电力系统对称短路故障的计算计算法程序设计设计内容电力系统故障的计算程序设计及编制和调试。设计要求1在对称短路计算、简单不对称短路计算中任选一种计算作为计算作业。2计算机语言自选。3设计、编制、调试出相关的通用计算程序。4输入输出数据一律以文件格式形成。5要求计算的题目（1）采用所编制的程序进行电力系统分析例6-3题，习题6-10、6-11对称短路计算；（2）采用所编制的程序进行电力系统分析例8-5题、习题8-1简单不对称短路计算；时间安排第10周以前参考资料1、社电力系统分析的计算机算法邱晓燕等 中国电力出版社2、电力系统基础 陈光会 王敏 中国水利电力出版社3、电力系统分析（上册） 何仰赞等 华中理工大学出版社4、电力系统稳态分析 陈珩 中国电力出版社5、matlab在电气工程中应用李维波 中国电力出版社6、matlab电力系统设计与分析吴天明等 国防工业出版目录1目的与原理41.1关于电力系统短路故障的计算机算法程序设计目的41.2原理42总体设计62.1程序设计说明62.2选择计算机语言63程序设计93.1流程图：93.2习题实例104总结18参考文献191目的与原理1.1关于电力系统短路故障的计算机算法程序设计目的电力系统正常运行的破坏多半是由于短路故障引起的，发生短路时，系统从一种状态剧变成另一种状态，并伴随复杂的暂态现象。所谓短路故障，是指一切不正常的相与相之间或相与地发生通路的情况。根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。1.2原理计算机计算原理应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三则是选择合适的计算机语言编制计算程序。建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂。运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用内存少，以利提高计算机的解题规模。选择合适的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错误。所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。最后将所编制出的正确计算程序，用于电力系统的实际计算。电力系统短路计算计算机算法一般在电力系统短路计算中，多数情况下只要计算短路电流、电压的周期分量起始值。因此，电力系统短路电流计算的数学模型，可归结为求解短路故障初始状态下的等值电路稳态解的问题。对于三相对称短路，可建立一相等值电路的计算模型，对于不对称短路，则可应用对称分量法建立系统的正、负、零序网络，从而建立故障计算的序网模型。（1）对称短路计算有关数学模型和原理框图以及已知结果的例题，参见电力系统分析第六章。常用的计算方法为节点导纳矩阵法或节点阻抗矩阵法，其形成方法分别参见电力系统分析第四章。（2）简单不对称短路计算有关数学模型、原理框图以及已知结果的例题参见电力系统分析第八章8-4节。2总体设计2.1程序设计说明根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学们对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。 2.2选择计算机语言matlab软件的基本功能在进行电力系统潮流计算，应用计算机算法进行求解节点导纳矩阵时候，我们要用到matlab的知识，所以应当对matlab的基本指令有所了解，下面是matlab的一些基本功能：进入matlab之后，会看到一个matlab command window，称为命令窗，它是最主要的窗口，既是键入命令也是显示计算结果的地方。另外还有一个编程窗，专门用来编辑应用程序。还有一个主窗口，用来记录已使用过的历史命令和已打开的目录，方便使用者查找。如果绘图还会自动弹出一个绘图窗，专门用来显示绘制的图形。matlab一般有3种进行计算的方法，第1种就如同使用计算器，直接输入数值和运算符，立即从屏幕上获得结果。第2种先对变量赋值，然后再输入由变量构成的表达式，也可立即获得结果。第3种，就是采用编程的方法来解决较复杂的，诸如含有判断、循环、迭代、递归等算法的较复杂的问题。上述方法中，第2和第3包括了数组和矩阵运算，只要定义了数组和矩阵变量，就可以如同普通代数运算一样直接用变量进行数学运算，十分方便。我们使用的是第四代计算机语言的matlab，利用其丰富的函数资源，它的优点如下： 1.语言简洁紧凑，使用方便灵活，库函数极其丰富。matlab程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。可以说，用matlab进行科技开发是站在专家的肩膀上。具有fortran和c等高级语言知识的读者可能已经注意到，如果用fortran或c语言去编写程序，尤其当涉及矩阵运算和画图时，编程会很麻烦。例如，如果用户想求解一个线性代数方程，就得编写一个程序块读入数据，然后再使用一种求解线性方程的算法（例如追赶法）编写一个程序块来求解方程，最后再输出计算结果。在求解过程中，最麻烦的要算第二部分。解线性方程的麻烦在于要对矩阵的元素作循环，选择稳定的算法以及代码的调试动不容易。即使有部分源代码，用户也会感到麻烦，且不能保证运算的稳定性。解线性方程的程序用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会根据矩阵的特性选择方程的求解方法，所以用户根本不用怀疑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的核心文件和工具箱文件都是可读可改的源文件，用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。3程序设计3.1流程图： 输入数据形成节点导纳矩阵选择故障点if输入数据计算节点阻抗矩阵if列元素用公式（6-10）计算短路电流if if=1/（zff+zf）z电流if用公式（6-11）计算各点电压 vi=1-zif/( zff+zf)用公式（6-9）计算指定支路的电路 ipq=(kvp-vq)/zpq输出结果matlab流程图开 机输入节点数n支路数n1数组bk=1k1=b(k,1)k2=b(k,2)y(k1,k1)=y(k1,k1)+1/b(k,3)+b(k,4)/2；y(k1,k2)=y(k1,k2)-1/(b(k,3)*b(k,5)；y(k2,k1)=y(k1,k2);y(k2,k2)=y(k2,k2)+1/(b(k,3)*b(k,5)2)+b(k,4)/2;k=k+1nk=n1？ y输出结果3.2习题实例【例6-3】在如图3-1所示的电力系统中分别在节点1和节点5接入发电机支路，其标幺值参数为：。在节点3发生三相短路，计算短路电流及网络中的电流分布。线路的电阻和电容略去不计，变压器的标幺变比等于1。各元件参数的标幺值如下： 图3-1电力系统等值网络图 图3-2 三相短路时的等值网络图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.9802jz=inv(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.08i运行结果如下：please input jd=3y =columns 1 through 40 -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.4348icolumn 5 0 0 0 0 + 5.4348i 0 - 9.9802ians =columns 1 through 40 + 0.0567i 0 - 0.0086i 0 - 0.0081i 0 - 0.0077i0 - 0.0086i 0 - 0.0146i 0 - 0.0137i 0 - 0.0130i0 - 0.0081i 0 - 0.0137i 0 + 0.0289i 0 + 0.0116i0 - 0.0077i 0 - 0.0130i 0 + 0.0116i 0 + 0.0306i0 - 0.0042i 0 - 0.0071i 0 + 0.0063i 0 + 0.0166icolumn 5 0 - 0.0042i 0 - 0.0071i 0 + 0.0063i 0 + 0.0166i 0 + 0.1093ii = 0 - 5.3763iv1 = 0.5151v2 = 0.1758v3 = 0v4 = 0.1339v5 = 0.5285i54 = 0 - 2.1447ii43 = 0 - 2.6774ii23 = 0 - 2.7047ii12 = 0 - 3.2309ii24 = 0 - 0.5242i【习题6-10】如图3-9所示网络，略去负荷，试用节点阻抗矩阵法求节点5发生三相短路时，短路点的短路电流及线路l-2、l-3的电流。已知各元件参数如下：发电机g-1: sn = 120 mva, x”d = 0.23; g-2: sn = 60 mva, x”d = 0.14变压器t-1: sn = 120 mva, vs = 10.5%; t-2: sn = 60 mva, vs = 10.5%线路参数 x1 = 0.4/km, b1 = 2.810-6 s/km线路长度 l-1: 120 km, l-2: 80 km, l-3: 70 km取sb = 120 mva, vb = vavggg-1g-2t-1t-2l-1l-2l-312345电力系统网络图三相短路时的等值网络图节点导纳矩阵如下所示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=);z=0 + 0.1819i 0 + 0.0616i 0 + 0.1600i 0 + 0.1079i 0 + 0.1326i; 0 + 0.0616i 0 + 0.2078i 0 + 0.0898i 0 + 0.1536i 0 + 0.1242i; 0 + 0.1600i 0 + 0.0898i 0 + 0.2330i 0 + 0.1571i 0 + 0.1932i; 0 + 0.1079i 0 + 0.1536i 0 + 0.1571i 0 + 0.2689i 0 + 0.2174i; 0 + 0.1326i 0 + 0.1242i 0 + 0.1932i 0 + 0.2174i 0 + 0.3427ii=1/0.3427iv1=1-0.1326i*iv2=1-0.1242i*iv3=1-0.1932i*iv4=1-0.2174i*iv5=0i54=(v5-v4)/0.2174ii42=(v4-v2)/0.1536ii21=(v2-v1)/0.0616ii43=(v4-v3)/0.1571ii32=(v3-v2)/0.0898i运行结果如下所示：y = 0 -13.8716i 0 0 + 9.5238i 0 0 0 0 - 8.3333i 0 0 + 4.7916i 0 0 + 9.5238i 0 0 -15.2329i 0 + 2.2960i 0 + 3.4440i0 0 + 4.7916i 0 + 2.2960i -10.9646 0 + 3.9360i0 0 0 + 3.4440i 0 + 3.9360i -7.3569 ans =-0.0085 + 0.1211i -0.0045 - 0.0080i -0.0124 + 0.0713i -0.0079 - 0.0139i -0.0260 - 0.0100i -0.0045 - 0.0080i -0.0226 + 0.1155i -0.0066 - 0.0116i -0.0394 - 0.0078i 0.0096 - 0.0241i -0.0124 + 0.0713i -0.0066 - 0.0116i -0.0180 + 0.1039i -0.0115 - 0.0202i -0.0378 - 0.0146i -0.0079 - 0.0139i -0.0394 - 0.0078i -0.0115 - 0.0202i -0.0685 - 0.0136i 0.0167 - 0.0420i -0.0260 - 0.0100i 0.0096 - 0.0241i -0.0378 - 0.0146i 0.0167 - 0.0420i -0.1066 - 0.0087ii = 3.1263v3 = 0 + 0.4632iv4 = 0 + 0.3890iil2 = 0 + 1.5953iil3 = 0 - 1.5311i【习题6-11】电力系统等值电路如图3-13所示，支路阻抗的标幺值已注明图中。（1）形成节点导纳矩阵（或节点阻抗矩阵），并用以计算节点3的三相短路电流。（2）另选一种方法计算短路电流，并用以验证（1）的计算结果。电力系统等值电路图最终计算等值电路图jd=input(please input jd=);y=-10.0j, 5.0j, 0; 5.0j, -14.5j, 10.0j; 0, 10.0j, -17.6667jinv(y)i1=1.05/0.2ji3=1.05/0.15jv3=0.98434if=v3/0.10714程序运行结果如下所示：please input jd=3y = 0 -10.0000i 0 + 5.0000i 0 0 + 5.0000i 0 -14.5000i 0 +10.0000i 0 0 +10.0000i 0 -1