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

PAGE11.课程设计目的及意义：根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。所谓短路，是指电力系统正常运行情况以外的相与相之间或相与地之间的“短接”。在电力系统正常运行时，除中性点外，相与相或相与地之间是绝缘的。如果由于某种原因使其绝缘破坏而构成了通路，我们就称电力系统是发生了短路故障。在三相系统中，短路故障可分为两大类：即对称短路（三相短路）和不对称短路（两相短路、两相接地短路、单相接地短路）。其中三相短路虽然发生的机会较少，但情况严重，又是研究其它短路的基础。所以我们先研究最简单的三相短路电流的暂态变化规律。2.课程设计内容电力系统故障的计算程序设计及编制和调试。3.计算机计算的原理我使用的是第四代计算机语言的MATLAB，MATLAB求解下列方程，并求解矩阵A的特征值。Ax=b,其中：A=32134567237985124323546598347135b=1234解为：x=A\b;设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,signlproceessingtoolbox,commumnicationtoolbox等。这些工具箱都是由该领域内学术水平很高的专家编写的，所以用户无需编写自己学科范围内的基础程序，而直接进行高，精，尖的研究。9）源程序的开放性。开放性也许是MATLAB最受人们欢迎的特点。除内部函数以外，所有MATLAB的核心文件和工具箱文件都是可读可改的源文件，用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。短路电流的构成及求解A.（1）.周期分流电路求解：（2）.非周期分量的求解：B.电机的过渡过程理论及短路电流的求解5．例题的求解：例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.000,Y42=J15.5000,Y43=-J37.9348,Y45=J5.4348,Y54=J5.4348,Y55=-J9.9802(1)电路图：(2)流程图：输入数据输入数据形成节点导纳矩阵形成节点导纳矩阵选择故障点If选择故障点If输入数据计算节点阻抗矩阵If输入数据计算节点阻抗矩阵If列元素用公式（6-10）计算短路电流IfIf=1/（Zff+zf）z用公式（6-10）计算短路电流IfIf=1/（Zff+zf）z电流If用公式（6-11）计算各点电压Vi=1-Zif/(Zff+zf)用公式（6-11）计算各点电压Vi=1-Zif/(Zff+zf)用公式（6-9）计算指定支路的电路Ipq=(用公式（6-9）计算指定支路的电路Ipq=(Kvp-Vq)/zpq输出结果输出结果Y=0-16.1905i0+9.5238i0000+9.5238i0-37.4084i0+15.3846i0+12.5000i000+15.3846i0-35.3846i0+20.0000i000+12.5000i0+20.0000i0-37.9348i0+5.4348i0000+5.4348i0-9.9802i采用4-3节的方法Z=0+0.1181i0+0.0958i0+0.0902i0+0.0858i0+0.0467i0+0.0958i0+0.1629i0+0.1533i0+0.1459i0+0.0794i0+0.0902i0+0.1533i0+0.1860i0+0.1611i0+0.0877i0+0.0858i0+0.1459i0+0.1611i0+0.1729i0+0.0941i0+0.0467i0+0.0794i0+0.0877i0+0.0941i0+0.1515i计算电流及网络中的电流分布If=Viˊ/Zii=1/ZiiV1ˊ=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-V2)/z12I24=(V2-V4)/z426．程序及说明程序如下：jd=input('pleaseinputjd=');ifjd==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.9802j]inv(Y)I=1/0.1181iV1=0V2=1-0.0958i*IV3=1-0.0902i*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.08ielseifjd==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.9802j]inv(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.08ielseifjd==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.5000j,20.000j,-37.9348j,5.4348j;0,0,0,5.4348j,-9.9802j]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.08ielseifjd==4Y=[-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.9802j]inv(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.08ielseifjd==5Y=[-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.9802j]inv(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.105iI24=(V2-V4)/0.08ielsec=errorEnd7.程序结果及分析节点1：pleaseinputjd=1y11=0-16.1905iy12=0+9.5238iy13=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-37.9348iy45=0+5.4348iy51=0y52=0y53=0y54=0+5.4348iy55=0-9.9802iy0=Columns1through50-16.1905i0+9.5238i0000+9.5238i0-37.4084i0+15.3846i0+12.5000i000+15.3846i0-35.3846i0+20.0000i000+12.5000i0+20.0000i0-37.9348i0+5.4348i0000+5.4348i0-9.9802ians=Columns1through40+0.1181i0+0.0958i0+0.0902i0+0.0858i0+0.0958i0+0.1629i0+0.1533i0+0.1459i0+0.0902i0+0.1533i0+0.1860i0+0.1611i0+0.0858i0+0.1459i0+0.1611i0+0.1729i0+0.0467i0+0.0794i0+0.0877i0+0.0941iColumn50+0.0467i0+0.0794i0+0.0877i0+0.0941i0+0.1515ii=0-8.4674iv1=0v2=0.1888v3=0.2362v4=0.2735v5=0.6046i54=0-1.7993ii43=0-0.7451ii23=0+0.7295ii12=0+1.7983ii24=0+1.0584i节点2：pleaseinputjd=2i=0-6.1387iv1=0.4119v2=0v3=0.0589v4=0.1044v5=0.5126i54=0-2.2186ii43=0-0.9085ii23=0+0.9066ii12=0