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

信息工程学院课程设计（论文）信息工程学院课程设计（论文）;0,15.3846j,-35.3846j,20.000j,0;0,0,12.5000j,0,0,20.000j,-37.9348j,5.4348j;5.4348j,-9.9802j]Z=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运行结果如下：pleaseinputjd=3Columns1through40-16.1905i0+9.5238i0+9.5238i0+37.4084i0+15.3846i0+12.5000i0+15.3846i0-35.3846i0+20.0000i0+12.5000i0+20.0000i0-37.9348iColumn50+5.4348i0-9.9802ians=Columns1through40+5.4348i0+0.0567i0-0.0086i0-0.0081i0-0.0077i0-0.0086i0-0.0146i0-0.0137i0-0.0130i0-0.0081i0-0.0137i0.0289i0.0116i0-0.0077i0-0.0130i0.0116i0;0,15.3846j,-35.3846j,20.000j,0;0,0,12.5000j,0,0,20.000j,-37.9348j,5.4348j;5.4348j,-9.9802j]Z=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运行结果如下：pleaseinputjd=3Columns1through40-16.1905i0+9.5238i0+9.5238i0+37.4084i0+15.3846i0+12.5000i0+15.3846i0-35.3846i0+20.0000i0+12.5000i0+20.0000i0-37.9348iColumn50+5.4348i0-9.9802ians=Columns1through40+5.4348i0+0.0567i0-0.0086i0-0.0081i0-0.0077i0-0.0086i0-0.0146i0-0.0137i0-0.0130i0-0.0081i0-0.0137i0.0289i0.0116i0-0.0077i0-0.0130i0.0116i0.0306i0-0.0042i0-0.0071i0.0063i0.0166iY=[0-16.1905j,9.5238j, 0,0+9.5238j,37.4084j,15.3846j,12.5000j,0;Column50-0.0042i0-0.0071i0+0.0063i0+0.0166i0+0.1093iI=0-5.3763i1=0.51512=0.17583=04=0.13395=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:S=120MVA,x”d=0.23;G-2: S=60MVA,x”d=0.14变压器T-1:S=120MVA,VS=10.5%;T-2: Sn=60MVA,VS=10.5%线路参数Xi=0.4Q/km,b1=2.8x10-6S/km线路长度L-1:120km,L-2:80km,L-3:70km取&=120MVA,VB=Mv

电力系统网络图节点导纳矩阵如下所示Y=-j13.87160 j9.5238 0 00 -j8.33330 j4.76190j9.5238 0 -j15.2329j2.2960j3.44400 j4.7619j2.2960 -j10.9646j3.9360设计程序如下：jd=input('pleaseinputjd=');Z=[0+0.1819i 0+0.0616i 0+0.1600i 0+0.1079i0+0.1326i;0+0.0616i 0+0.2078i 0+0.0898i 0+0.1536i0+0.1242i;

0+0.1600i0.1932i;0+0.1079i0.2174i;0+0.1326i0.3427i]I=1/0.3427iV1=1-0.1326i*I0+0.1600i0.1932i;0+0.1079i0.2174i;0+0.1326i0.3427i]I=1/0.3427iV1=1-0.1326i*IV2=1-0.1242i*IV3=1-0.1932i*IV4=1-0.2174i*IV5=00+0.0898i0+0.1536i0+0.1242i0+0.2330i0+0.1571i0+0.1932i0+0.1571i0+0+0.2689i0+0+0.2174i0+I54=(V5-V4)/0.2174iI42=(V4-V2)/0.1536iI21=(V2-V1)/0.0616iI43=(V4-V3)/0.1571iI32=(V3-V2)/0.0898i运行结果如下所示:Y=0-13.8716i00 00+9.5238i 0+3.4440i0 0Y=0-13.8716i00 00+9.5238i 0+3.4440i0 03.9360i0 0-7.3569ans=0+9.5238i-8.3333i 0 00-15.2329i 0+2.2960i+4.7916i+4.7916i 0 +2.2960i -10.9646 00+3.4440i 0+3.9360i-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.0078i0.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.0136i0.0167-0.0420i-0.0260-0.0100i0.0096-0.0241i-0.0378-0.0146i0.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('pleaseinputjd=');Y=[-10.0j, 5.0j, 0;5.0j, -14.5j, 10.0j;0, 10.0j, -17.6667j]inv(Y)I1=1.05/0.2jI3=1.05/0.15jV3=0.98434If=V3/0.10714程序运行结果如下所示：pleaseinputjd=3最终计算等值电路图jd=input('pleaseinputjd=');Y=[-10.0j, 5.0j, 0;5.0j, -14.5j, 10.0j;0, 10.0j, -17.6667j]inv(Y)I1=1.05/0.2jI3=1.05/0.15jV3=0.98434If=V3/0.10714程序运行结果如下所示：pleaseinputjd=30-10.0000i 0+5.0000i 00+5.0000i0-14.5000i0+10.0000ians=0+0.1394i0+0.0789i0+0.0446i0+10.0000i0-17.6667i0+0.0789i0+0.1577i0+0.0893i0+0.0446i0+0.0893i0+0.1071iI1=0-5.2500iI3=0-7.0000iV3=0.9843If=9.18744总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。同时由于求解过程中用到求节点导纳矩阵，求矩阵的逆等等，又使我对以前所学的知识有了一次很好的温习， 我通过自己看书查资料，自己把以前比较零碎的知识点理论在实践设计中加以整合， 明白了以前学习中的很多不足。通过这次课程设计，既弥补了知识点的疏漏，同时也看到了研究性学习的效果，从研究中学习，理论结合实际，将理论运用到实际，同时在实践中发现问题，然后解决问题。而且在此次课程设计中，我发现了自己的基础知识有很多的不足。这些基础的缺乏给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认真去做，持之以恒，就会有新的发现，有意外的收获。对于这次的课程设计依然有许多不足之处，今后我会把 MATLAB语言学好。这次我在MATLAB电子仿真与应用教程中学到了很多东西