1、 课程设计报告( 20142015 年度第二学期)名 称：电力系统暂态上机计算院 系： 电气与电子工程学院 班 级： 电气 1211 学 号： 1121210205 学生姓名： 郝阳 指导教师： 陈艳波 设计周数： 两周 成 绩： 日期： 2015 年 7 月 4 日课程 电力系统暂态上机课程设计报告1一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、设计正文（详细内容见附录,用 A4 纸，页数不限）1 对称短路计算过程流程图和计算结果2 不对
2、称短路计算过程流程图和计算结果3 静态稳定性计算过程流程图和计算结果4 暂态稳定性计算过程流程图和计算结果5 思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的
3、分析方法是采用对称分量法将不对称问题转化为对称问题来求解。同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。这与手算系统短路电流时使用的网络化简方法大大不同。在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于 matlab 软件的计算。利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在 matlab 软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用
4、起来简便易懂。利用改进欧拉法计算最大切除角或切除时间虽然比较方便，但是由于计算机有效位数的限制而引起的舍入误差随步长的减小以及运算次数的增多而增大。在上机编程中体会比较深的是自己犯的几个错误，比如说在算短路电流时，发电机之路的电流和电压是否要经过相位的变换，在静态稳定计算中的角度是弧度制，在计算中使用的都是标幺值等等，这些看似细微的地方，往往考验的就是自己对电力系统分析中的知识的基础，虽然小但是也非常的重要，往往就是这些错误会令程序得不到正确的结果，值得我注意和反思。总的来说，结合我们上学期潮流上机编程的经验，虽然这次用的是 matlab，但是明显感觉到对计算机处理电力系统问题的简便，程序的编
5、写也显得更加的轻松和得心应手，体会到 matlab 软件的强大。经过两个星期的回顾与学习，我巩固了电力系统暂态分析理论知识，在编程的过程中很深刻的感受到要想得出正确的结果，就必须认真地理解课本上对应的例题，知道这些原理，不能仅仅照着书抄公式；同时我也初步掌握了采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础，更好地建设我国的电力系统。四、参考文献1. 电力系统暂态过程 ，常鲜戎、赵书强编，机械工业出版社，2010 年 1 月，第一版；2. 电力系统分析基础 ，李庚银，机械工业出版社，2011 年，第一版；课程
6、电力系统暂态上机课程设计报告2附录（设计流程图、计算结果、思考题答案）1. 对称短路计算过程流程图和计算结果流程图：输入数据选择故障类型用节点导纳阵求逆得到节点阻抗矩阵计算各节点的序电压、相电压根据故障类型计算故障点序电流、相电流计算各支路的序电流、相电流计算发电机节点的序电压、相电压形成节点导纳矩阵选择短路点开始计算发电机支路的序电流、相电流结束课程 电力系统暂态上机课程设计报告3计算结果：（1） 导纳矩阵：正序阻抗： 10. -4.i0. +1.0i. +1.0i. +20.i.0 +.i+-6 4. .i. .i.Y -2.i. .i. .i2 -30. .i0. 4.0i. .0i.
7、0.i.0 -6.i 负序阻抗： 2. -4.i. +1.i. +1.i. +2.i. +.i+1-6 40. .0i. .0i.Y -20.i.0 .i.0 .i-3. .i. 4.i. .i. .i. -6.i 零序阻抗： 0. -3.0i. +5.0i. +5.0i+5-. .i. .i. -1.iY （2） 节点 3 发生三相短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 9.8592i Ia = 9.8592I2 = 0 Ib = 9.8592I0 = 0 Ic = 9.8592b. 各节点三序电压、三相电压有效值：U1 U2 U0 Ua Ub Uc1 0.
8、4507 0 0 0.4507 0.4507 0.45072 0.5352 0 0 0.5352 0.5352 0.53523 0.0000 0 0 0.0000 0.0000 0.00004 0.6338 0 0.6338 0.6338 0.63385 0.6901 0 0.6901 0.6901 0.6901c. 各支路三序电流、三相电流有效值：电压节点课程 电力系统暂态上机课程设计报告4I1 I2 I0 Ia Ib Ic1-20.0000 + 0.8451i0 0 0.8451 0.8451 0.84511-30.0000 - 4.5070i0 0 4.5070 4.5070 4.507
9、02-30.0000 - 5.3521i0 0 5.3521 5.3521 5.35211-40.0000 + 3.6620i0 0 3.6620 3.6620 3.66202-50.0000 + 6.1972i0 0 6.1972 6.1972 6.1972G10.0000 - 3.6620i0 3.6620 3.6620 3.6620G20.0000 - 6.1972i0 6.1972 6.1972 6.19722. 不对称短路计算过程流程图和计算结果（1） 节点 3 发生 A 相短路接地故障：a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 3.1167i Ia = 9.
10、3500I2 = 0.0000 - 3.1167i Ib = 0.0000I0 = 0.0000 - 3.1167i Ic = 0.0000b. 各节点三序电压、三相电压有效值：U1 U2 U0 Ua Ub Uc1 0.8264 -0.1736 -0.0686 0.5841 0.9518 0.95182 0.8531 -0.1469 -0.0436 0.6625 0.9526 0.95263 0.6839 -0.3161 -0.3678 0.0000 1.0268 1.02684 0.8842 -0.1158 0.8324 1.0000 0.83245 0.9020 -0.0980 0.857
11、3 1.0000 0.8573c. 各支路三序电流、三相电流有效值：电流支路节点电压课程 电力系统暂态上机课程设计报告5I1 I2 I0 Ia Ib Ic1-20.0000 + 0.2671i0.0000 + 0.2671i0.0000 + 0.1247i0.6589 0.1425 0.14251-30.0000 - 1.4248i0.0000 - 1.4248i0.0000 - 1.4960i4.3455 0.0712 0.07122-30.0000 - 1.6919i0.0000 - 1.6919i0.0000 - 1.6207i5.0045 0.0712 0.07121-40.0000
12、+ 1.1576i0.0000 + 1.1576i0.0000 + 1.3713i3.6866 0.2137 0.21372-50.0000 + 1.9590i0.0000 + 1.9590i0.0000 + 1.7453i5.6634 0.2137 0.2137G10.0000 - 1.1576i0.0000 - 1.1576i 2.0050 0.0000 2.0050G20.0000 - 1.9590i0.0000 - 1.9590i 3.3932 0.0000 3.3932（2） 节点 3 发生 A 相经 10 电阻接地故障a. 故障点三序电流、三相电流有效值：I1 = 0.9793 -
13、 2.7705i Ia = 8.8154I2 = 0.9793 - 2.7705i Ib = 0.0000I0 = 0.9793 - 2.7705i Ic = 0.0000b. 各节点三序电压、三相电压有效值：U1 U2 U0 Ua Ub Uc10.8456 - 0.0546i-0.1544 - 0.0546i-0.0610 - 0.0215i0.6437 0.9269 0.986720.8694 - 0.0462i-0.1306 - 0.0462i-0.0388 - 0.0137i0.7080 0.9281 0.986930.7190 - 0.0993i-0.2810 - 0.0993i-0
14、.3269 - 0.1156i0.3333 1.0375 1.010040.8971 - 0.0364i-0.1029 - 0.0364i 0.8149 1.0000 0.8888节点电压支路 电流课程 电力系统暂态上机课程设计报告650.9129 - 0.0308i-0.0871 - 0.0308i 0.8432 1.0000 0.9043c. 各支路三序电流、三相电流有效值：I1 I2 I0 Ia Ib Ic1-2-0.0839 + 0.2375i-0.0839 + 0.2375i-0.0392 + 0.1108i0.6213 0.1343 0.13431-30.4477 - 1.2665
15、i0.4477 - 1.2665i0.4701 - 1.3298i4.0971 0.0672 0.06722-30.5316 - 1.5040i0.5316 - 1.5040i0.5093 - 1.4406i4.7183 0.0672 0.06721-4-0.3638 + 1.0290i-0.3638 + 1.0290i-0.4309 + 1.2190i3.4758 0.2015 0.20152-5-0.6156 + 1.7414i-0.6156 + 1.7414i-0.5484 + 1.5515i5.3396 0.2015 0.2015G10.3638 - 1.0290i0.3638 - 1
16、.0290i 1.8904 0.0000 1.8904G20.6156 - 1.7414i0.6156 - 1.7414i 3.1992 0.0000 3.1992（3） 节点 3 发生 b、c 两相相间短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 4.9296i Ia = 0I2 = 0.0000 + 4.9296i Ib = 8.5383I0 = 0 Ic = 8.5383b. 各节点三序电压、三相电压有效值：U1 U2 U0 Ua Ub Uc1 0.7254 0.2746 0 1.0000 0.6343 0.63432 0.7676 0.2324 0 1.00
17、00 0.6818 0.68183 0.5000 0.5000 0 1.0000 0.5000 0.50004 0.8169 0.1831 0.9222 0.6338 0.92225 0.8451 0.1549 0.9322 0.6901 0.9322支路 电流节点电压课程 电力系统暂态上机课程设计报告7c. 各支路三序电流、三相电流有效值：I1 I2 I0 Ia Ib Ic1-20.0000 + 0.4225i0.0000 + 0.4225i0 0.0000 0.7319 0.73191-30.0000 - 2.2535i0.0000 - 2.2535i0 0.0000 3.9032 3.9
18、0322-30.0000 - 2.6761i0.0000 - 2.6761i0 0.0000 4.6351 4.63511-40.0000 + 1.8310i0.0000 + 1.8310i0 0.0000 3.1714 3.17142-50.0000 + 3.0986i0.0000 + 3.0986i0 0.0000 5.3669 5.3669G10.0000 - 1.8310i0.0000 - 1.8310i 1.8310 3.6620 1.8310G20.0000 - 3.0986i0.0000 - 3.0986i 3.0986 6.1972 3.0986（4） 节点 3 发生 b、c
19、两相接地短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 6.4114i Ia = 0.0000I2 = 0.0000 + 3.4478i Ib = 9.6262I0 = 0.0000 + 2.9636i Ic = 9.6262b. 各节点三序电压、三相电压有效值：U1 U2 U0 Ua Ub Uc1 0.6428 0.1921 0.0652 0.9001 0.5258 0.52582 0.6977 0.1625 0.0415 0.9018 0.6049 0.60493 0.3497 0.3497 0.3497 1.0491 0.0000 0.00004 0.7619
20、0.1281 0.8333 0.6338 0.83335 0.7985 0.1084 0.8578 0.6901 0.8578支路 电流节点电压课程 电力系统暂态上机课程设计报告8c. 各支路三序电流、三相电流有效值：I1 I2 I0 Ia Ib Ic1-20.0000 + 0.5495i0.0000 - 0.2955i0.0000 - 0.1185i0.1355 0.7719 0.77191-30.0000 - 2.9309i0.0000 + 1.5761i0.0000 + 1.4225i0.0677 4.4322 4.43222-30.0000 - 3.4805i0.0000 + 1.87
21、17i0.0000 + 1.5411i0.0677 5.1947 5.19471-40.0000 + 2.3814i0.0000 - 1.2806i0.0000 - 1.3040i0.2032 3.6737 3.67372-50.0000 + 4.0300i0.0000 - 2.1672i0.0000 - 1.6596i0.2032 5.9596 5.9596G10.0000 - 2.3814i0.0000 + 1.2806i 2.0643 3.6620 2.0643G20.0000 - 4.0300i0.0000 + 2.1672i 3.4934 6.1972 3.4934支路 电流课程 电
22、力系统暂态上机课程设计报告93. 静态稳定性计算过程流程图和计算结果流程图：否否否是是是开始清屏设初值：Xd_all,Xdd_all,Xq,Xq_all,Eq0,UG0,Xe,Xd,XddKe=0.1Ke_result=Ke,delta_result=delta,P_result=Eq*U*sin(delta*pi/180)/Xd_all, delta=49求 A 阵特征值求 Eq,Eqq求 UGd,UGq,UG求 K1K6Ke5.7delta=110(V10)end%-计算节点导纳矩阵-%Y1=zeros(nbus_1);for m=1:nbranch_1from_bus=branch_1(
23、m).from;to_bus=branch_1(m).to; %write your programme hereNo1=from_bus;No2=to_bus;if No2=0Y1(No1,No1)=Y1(No1,No1)+1/branch_1(m).x;elseY1(No1,No1)=Y1(No1,No1)+1/branch_1(m).x;Y1(No2,No2)=Y1(No2,No2)+1/branch_1(m).x;Y1(No1,No2)=-1/branch_1(m).x;Y1(No2,No1)=-1/branch_1(m).x;endenddisp(正序阻抗);Y1Y2=Y1; %负序
24、等于正序disp(负序阻抗);Y2Y0=zeros(nbus_0);for m=1:nbranch_0from_bus=branch_0(m).from;to_bus=branch_0(m).to; %write your programme hereNo1=from_bus;No2=to_bus;if No2=0Y0(No1,No1)=Y0(No1,No1)+1/branch_0(m).x0;elseY0(No1,No1)=Y0(No1,No1)+1/branch_0(m).x0;Y0(No2,No2)=Y0(No2,No2)+1/branch_0(m).x0;课程 电力系统暂态上机课程设计
25、报告17Y0(No1,No2)=-1/branch_0(m).x0;Y0(No2,No1)=-1/branch_0(m).x0;endenddisp(零序阻抗);Y0%-计算节点阻抗矩阵-%Z1=inv(Y1);Z2=inv(Y2);Z0=inv(Y0);return;（2） shortcircuit.mfunction shortcircuit()%数据来源于教材电力系统分析基础P77 例（3-4），P182 例（7-3）%-输入已知条件-%bus_Num_1=5; %正负序网络的节点数bus_Num_0=3; %零序网络的节点数bus(1) = struct(idx,1,shift,0);
26、bus(2) = struct(idx,2,shift,0);bus(3) = struct(idx,3,shift,0);bus(4) = struct(idx,4,shift,pi/6);bus(5) = struct(idx,5,shift,pi/6);branch_Num_1=7; %正负序网络的支路数branch_Num_0=5; %零序网络的支路数branch_1(1) = struct(from,1,to,2,x,j*0.1,shift,0);branch_1(2) = struct(from,1,to,3,x,j*0.1,shift,0);branch_1(3) = struc
27、t(from,2,to,3,x,j*0.1,shift,0);branch_1(4) = struct(from,1,to,4,x,j*0.05,shift,0);branch_1(5) = struct(from,2,to,5,x,j*0.025,shift,0);branch_1(6) = struct(from,4,to,0,x,j*0.1,shift,pi/6);branch_1(7) = struct(from,5,to,0,x,j*0.05,shift,pi/6);branch_0(1) = struct(from,1,to,2,x0,j*0.2);branch_0(2) = st
28、ruct(from,1,to,3,x0,j*0.2);branch_0(3) = struct(from,2,to,3,x0,j*0.2);branch_0(4) = struct(from,1,to,0,x0,j*0.05);branch_0(5) = struct(from,2,to,0,x0,j*0.025);%-计算正序和负序导纳矩阵-%Z1,Z0=formAdmitanceMatrix(bus_Num_1,bus_Num_0,branch_Num_1,branch_Num_0,branch_1,branch_0); %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量Z2=Z1;课程 电力
29、系统暂态上机课程设计报告18Fault_Node=input(输入短路点编号：;n Fault_Node=);Fault_Type=input(输入短路类型：;n(1)Fault_Type=0 为三相短路;n(2)Fault_Type=1 为 a相接地短路;n(3)Fault_Type=2 为 a 相经 10 欧姆电阻接地短路n(4)Fault_Type=3 为 bc 两相相间短路n(5)Fault_Type=4 为 bc 两相短路接地nFault_Type=);a=-0.5+j*sqrt(3)/2;T=1 1 1 % T 为对称分量法的变换矩阵，见 P87 公式（4-4）a2 a 1a a2
30、 1;%-第一步：计算短路点的序电流，相电流-%计根据故障类型选择不同的计算公式，计算故障点各序电流if Fault_Type=0 %三相短路I_Fault1=1/Z1(Fault_Node,Fault_Node); %三相短路只有正序电流I_Fault2=0;I_Fault0=0;elseif Fault_Type=1 %a 相接地短路%write your programme hereI_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); %3 个序网阻抗串联I_F
31、ault2=I_Fault1;I_Fault0=I_Fault1;elseif Fault_Type=2 %a 相经 10 欧姆电阻接地短路%write your programme hereZf=10*50/(1152);I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+3*Zf); %3 个序阻抗和电阻串联I_Fault2=I_Fault1;I_Fault0=I_Fault1;elseif Fault_Type=3 %bc 两相相间短路%write your p
32、rogramme hereI_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node); %正序阻抗和负序阻抗并联I_Fault2=-I_Fault1;I_Fault0=0;elseif Fault_Type=4 %bc 两相短路接地%write your programme here%三个序阻抗并联课程 电力系统暂态上机课程设计报告19I_Fault1=1/(Z1(Fault_Node,Fault_Node)+(Z2(Fault_Node,Fault_Node)*Z0(Fault_Node,Fault_Node)/(Z2(Fa
33、ult_Node,Fault_Node)+Z0(Fault_Node,Fault_Node);I_Fault2=-I_Fault1*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node);I_Fault0=-I_Fault1*Z2(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node);end;disp(短路点序电流);I1=I_Fault1I2=I_Fault2I0=I_Fault0disp(短路
34、电流有效值);Iabc=T*I_Fault1;I_Fault2;I_Fault0; %相量Iabc_effective=abs(Iabc) %有效值%-第二步：计算各个节点的序电压，相电压-%write your programme here%故障分量 Ifault_node1=zeros(bus_Num_1,1);Ifault_node2=zeros(bus_Num_1,1);Ifault_node0=zeros(bus_Num_0,1);for m=1:bus_Num_1 if m=Fault_NodeIfault_node1(m,1)=-1*I_Fault1;Ifault_node2(m
35、,1)=-1*I_Fault2;elseIfault_node1(m,1)=0;Ifault_node2(m,1)=0;endendfor m=1:bus_Num_0 if m=Fault_NodeIfault_node0(m,1)=-1*I_Fault0;elseIfault_node0(m,1)=0;endendUfault_node1=zeros(bus_Num_1,1);Ufault_node2=zeros(bus_Num_1,1);Ufault_node0=zeros(bus_Num_0,1);课程 电力系统暂态上机课程设计报告20Ufault_node1=Z1*Ifault_nod
36、e1;Ufault_node2=Z2*Ifault_node2;Ufault_node0=Z0*Ifault_node0;%正常分量 Unormal_node1=ones(bus_Num_1,1);Unormal_node2=zeros(bus_Num_1,1);Unormal_node0=zeros(bus_Num_0,1);U_node1=zeros(bus_Num_1,1);U_node2=zeros(bus_Num_1,1);U_node0=zeros(bus_Num_0,1);U_node1=Ufault_node1+Unormal_node1; %正常分量加故障分量U_node2=
37、Ufault_node2+Unormal_node2;U_node0=Ufault_node0+Unormal_node0;disp(各节点序电压);U1=U_node1U2=U_node2U0=U_node0Uabc=;for m=1:bus_Num_1if m=bus_Num_0tmp=T*U1(m)*exp(j*bus(m).shift);U2(m)*exp(-j*bus(m).shift);U0(m);Uabc=Uabc,tmp;elsetmp=T*U1(m)*exp(j*bus(m).shift);U2(m)*exp(-j*bus(m).shift);0;Uabc=Uabc,tmp;
38、endenddisp(各节点相电压有效值);Uabc_effective=abs(Uabc)% -第三步：计算各个支路的序电流，相电流-%write your programme hereI_branch1=zeros(bus_Num_1,1);I_branch2=zeros(bus_Num_1,1);I_branch0=zeros(bus_Num_0,1);Iabc1=;for m=1:branch_Num_0from_bus = branch_1(m).from;to_bus=branch_1(m).to;if(from_bus=0)I_branch1(m,1)=(0-U1(to_bus)
39、/branch_1(m).x;%支路正序电流I_branch2(m,1)=(0-U2(to_bus)/branch_1(m).x;%支路负序电流 课程 电力系统暂态上机课程设计报告21elseif(to_bus=0)I_branch1(m,1)=(U1(from_bus)-0)/branch_1(m).x;%支路正序电流I_branch2(m,1)=(U2(from_bus)-0)/branch_1(m).x;%支路负序电流 elseI_branch1(m,1)=(U1(from_bus)-U1(to_bus)/branch_1(m).x;%支路正序电流I_branch2(m,1)=(U2(f
40、rom_bus)-U2(to_bus)/branch_1(m).x;%支路负序电流 endfrom_bus = branch_0(m).from;to_bus=branch_0(m).to;if(from_bus=0)I_branch0(m,1)=(0-U0(to_bus)/branch_0(m).x0;%支路零序电流elseif(to_bus=0)I_branch0(m,1)=(U0(from_bus)-0)/branch_0(m).x0;%支路零序电流elseI_branch0(m,1)=(U0(from_bus)-U0(to_bus)/branch_0(m).x0;%支路零序电流end%
41、-考虑变压器支路造成的相移-%tmp=T*I_branch1(m,1)*exp(j*branch_1(m).shift);I_branch2(m,1)*exp(-j*branch_1(m).shift);I_branch0(m,1);%由序分量合成相分量Iabc1=Iabc1,tmp;end;for m=branch_Num_0+1:branch_Num_1from_bus = branch_1(m).from;I_branch1(m,1)=(1-U1(from_bus)/branch_1(m).x;%支路正序电流I_branch2(m,1)=(0-U2(from_bus)/branch_1(
42、m).x;%支路负序电流 %-考虑变压器支路造成的相移-%tmp=T*I_branch1(m,1)*exp(j*branch_1(m).shift);I_branch2(m,1)*exp(-j*branch_1(m).shift);0;%由序分量合成相分量Iabc1=Iabc1,tmp;end;disp(各支路序电流);I1= I_branch1I2= I_branch2I0= I_branch0disp(各支路相电流);Iabc1_effective=abs(Iabc1)return;（3） stability_smallsignal.mfunction stability_smallsig
43、nal()U=1; %系统电压Eq0=1.972; %空载电动势UG0=1.21; %机端电压课程 电力系统暂态上机课程设计报告22Xe=0.504; %线路电抗值Xd=0.982; %同步电抗Xdd=0.344; %暂态电抗Xd_all=1.486; %系统电抗(包含同步电抗)Xdd_all=0.848; %系统电抗(包含暂态电抗)Tj=10;Td0=10;Ke=0.1;Ke_result=;delta_result=;P_result=;while(Ke5.7)for delta=49:110%write your programme herea=(Ke*Ke*Xe*Xe)/(Xd_all
44、*Xd_all)-1;b=2*(Ke*Ke*Xd*Xe*U*cos(delta*pi/180)/(Xd_all*Xd_all)+Eq0+Ke*UG0);c=Ke*Ke*Xd*Xd*U*U/(Xd_all*Xd_all)-(Eq0+Ke*UG0)*(Eq0+Ke*UG0);Eq=(-b+sqrt(b*b-4*a*c)/(2*a);Eqq=Eq*Xdd_all/Xd_all+(Xd_all-Xdd_all)*U*cos(delta*pi/180)/Xd_all;UGd=U*Xd*sin(delta*pi/180)/Xd_all;UGq=Eq*Xe/Xd_all+Xd*U*cos(delta*pi/
45、180)/Xd_all;UG=sqrt(UGd*UGd+UGq*UGq); K1=Eqq*U*cos(delta*pi/180)/Xdd_all+U*U*(Xdd_all-Xd_all)*cos(2*delta*pi/180)/(Xd_all*Xdd_all);K2=U*sin(delta*pi/180)/Xdd_all;K3=Xdd_all/Xd_all;K4=(Xd_all-Xdd_all)*U*sin(delta*pi/180)/Xdd_all;K5=UGd*Xd*cos(delta*pi/180)/(UG*Xd_all)-UGq*Xdd*sin(delta*pi/180)/(UG*Xd
46、d_all);K6=UGq*(Xd_all-Xd)/(UG*Xdd_all); A=zeros(3);A(1,1)=0;A(1,2)=100*pi;A(1,3)=0;A(2,1)=-K1/Tj;A(2,2)=0;A(2,3)=-K2/Tj;A(3,1)=-(K4+Ke*K5)/Td0;A(3,2)=0;A(3,3)=-(1/K3+Ke*K6)/Td0;V=eig(A);V1=real(V(1);课程 电力系统暂态上机课程设计报告23V2=real(V(2);V3=real(V(3);%write your programme hereif (V10)elseKe_result=Ke_resul
47、t;Ke;delta_result=delta_result;delta;P=Eq*U*sin(delta*pi/180)/Xd_all;P_result=P_result;P;s=sprintf(Ke=%f, delta=%d, P=%fn,Ke,delta,P);disp(s);break;end;end; Ke=Ke+0.1; end;subplot(2,1,1),plot(Ke_result,P_result);xlabel(Ke);ylabel(Pmax);subplot(2,1,2),plot(Ke_result,delta_result);xlabel(Ke);ylabel(deltamax);return;（4） stability_transient.mfunction stability_transient(CutTime)if nargin 1CutTime=input(输入切除时间：;n CutTime=);endf=50; %系统额定频率Tj=8.47; %归算后的发电机惯性时间常数PT=1; %正常运行时发电机向无穷大系统传输的有功功率P2M=0.48; %故障存在时发电机的最大功率P3M=1.38; %故障切除后发
