电力系统暂态上机计算课程设计报告--附程序

1、精选优质文档-倾情为你奉上 课程设计报告( 20112012年度第二学期)名 称：电力系统暂态上机 院 系：电气与电子工程学院 班 级：电 气 0914 班 学 号： 学生姓名： 雷 刚 指导教师： 麻 秀 范 老 师 设计周数： 两 周 成 绩： 日期： 2012年 7月3日专心-专注-专业一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、设计正文（详细内容见附录,用A4纸，页数不限）1 对称短路计算过程流程图和计算结果2 不对称短路计算过程

2、流程图和计算结果3 静态稳定性计算过程流程图和计算结果4 暂态稳定性计算过程流程图和计算结果5 思考题三、课程设计总结或结论在这次电力系统暂态上机课程中，我们主要讨论了各种类型短路故障下（三相短路、单相接地短路、两相相间短路、两相接地短路、单相经电阻接地）系统网络中的电压、电流分布的计算以及电力系统的静态稳定、暂态稳定问题的分析。通过本次课程设计，我对之前学习的知识有了一个巩固和提高，同时对电力系统故障分析有了更深刻的理解。电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以

3、简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。这与手算系统短路电流时使用的网络化简方法大大不同。同时也明白了计算机编程方法对于电力系统稳态暂态计算的简便性、快速性和重要性。最后，经过此次课程设计，对于我之前在数学建模中的学习到的Matlab知识起到了升华的作用，可以说这也算是一次比赛把，对较复杂的选择语句、循环语句的使用等有了更深的了解，为今后学习、工作中埋下一个伏笔，相

4、信会受益匪浅。四、参考文献1. 电力系统暂态分析，李光琦，中国电力出版社，2006年，第三版；2. 电力系统分析（上、下），何仰赞，华中科技大学出版社，1998年，第二版。3. 电力系统故障的计算机辅助分析 重庆大学出版社 米麟书等 4. 电力系统潮流计算 天津大学出版社 宋文南等 5. 电力系统故障分析 清华大学出版社 周荣光6. 短路电流实用计算方法 电力工业出版社 西安交通大学等7. 精通Matlab6.5北京航空航天大学出版社，张志涌附录（设计流程图、计算结果、思考题答案）1. 对称短路计算过程流程图和计算结果流程图如下：输入数据根据Ifault值选择故障类型用节点导纳阵求逆得到节点阻

5、抗矩阵选择短路点f计算各序电流的故障分量计算故障点各序电流计算故障点各序电压的故障分量计算故障后各节点电压值形成节点导纳矩阵确定Ifault值开始计算支路电流计算计及发电机端点的各节点电压结束计算结果如下：(1)导纳矩阵:Y1 =Y2=0 -26.6667i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -33.3333i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000iY0= 0 -29.0476i 0 + 5.0000i 0 + 5.0000i 0 + 5.0000i 0 -48.0228i 0 + 5.0000i

6、 0 + 5.0000i 0 + 5.0000i 0 -10.0000iYY1=YY2= 0 -40.0000i 0 +10.0000i 0 +10.0000i 0 +20.0000i 0 0 +10.0000i 0 -60.0000i 0 +10.0000i 0 0 +40.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000i 0 0 0 +20.0000i 0 0 0 -30.0000i 0 0 0 +40.0000i 0 0 0 -60.0000i(2) 节点3发生三相短路故障l 故障点三序电流、三相电流、三相电流有效值：I1= 0 - 9.8592i I2

7、=0 I3=0 ; Ia= 0 - 9.8592i Ib=-8.5383 + 4.9296i Ic=8.5383 + 4.9296iIa_effective = Ib_effective = Ic_effective = 9.8592l 各节点三序电压、三相电压：节点U1U2U0UaUbUc10.4507 000.4507 -0.2254 - 0.3903i -0.2254 + 0.3903i 20.5352000.5352-0.2676 - 0.4635i-0.2676 +0.4635i300000040.6338000.6338-0.3169 - 0.5489i-0.3169 +0.548

8、9i 50.6901000.6901-0.3451 - 0.5977i-0.3451 +0.5977il 各支路三序电流、三相电流：支路I1I2I0IaIbIc1-20 + 0.8451i000 - 0.8451i -0.7319 +0.4225i 0.7319 +0.4225i 1-30 - 4.5070i 000 + 4.5070i 3.9032 - 2.2535i -3.9032 - 2.2535i 2-30 - 5.3521i000 + 5.3521i4.6351 - 2.6761i-4.6351 - 2.6761i2不对称短路计算过程流程图和计算结果不对称短路计算过程流程图如三相短路

9、流程图。 (1) 节点3发生A相短路接地故障l 故障点三序电流、三相电流：I1= 0 - 3.1239i I2= 0 - 3.1239i I3= 0 - 3.1239i ; Ia= 0 - 9.3718i Ib= 0 - 0.0000i Ic= 0 - 0.0000i ；l 各节点三序电压、三相电压：节点U1U2U0UaUbUc10.8286-0.1714-0.07190.5853-0.4004 - 0.8660i-0.4004 + 0.8660i20.8552-0.1448-0.04620.6643-0.4014 - 0.8660i-0.4014 + 0.8660i30.6857-0.314

10、3-0.37140.0000-0.5571 - 0.8660i-0.5571 + 0.8660i40.7704 + 0.4448i-0.0956 + 0.0552i00.6748 + 0.5000i0.0000 - 1.0000i-0.6748 + 0.5000i50.7853 + 0.4534i-0.0807 + 0.0466i00.7046 + 0.5000i0.0000 - 1.0000i-0.7046 + 0.5000il 各支路三序电流、三相电流：支路I1I2I0IaIbIc1-20 + 0.2665i0 + 0.2665i0 + 0.1286i0 + 0.6616i0 - 0.13

11、79i0 - 0.1379i1-30 - 1.4287i0 - 1.4287i0 - 1.4977i0 - 4.3551i0.0000 - 0.0689i-0.0000 - 0.0689i2-30 - 1.6952i0 - 1.6952i0 - 1.6263i0 - 5.0167i0.0000 + 0.0689i-0.0000 + 0.0689i(2) 节点3发生A相经10电阻接地故障l 故障点三序电流、三相电流：I1= 0.0333 - 0.0004i I2= 0.0333 - 0.0004i I3= 0.0333 - 0.0004i ; Ia= 0.1000 - 0.0011i Ib= 0

12、 Ic= 0 ；l 各节点三序电压、三相电压：节点U1U2U0UaUbUc11.0000 - 0.0018i-0.0000 - 0.0018i-0.0000 - 0.0008i1.0000 - 0.0044i-0.5000 - 0.8650i-0.5000 + 0.8671i21.0000 - 0.0015i-0.0000 - 0.0015i-0.0000 - 0.0005i1.0000 - 0.0036i-0.5000 - 0.8650i-0.5000 + 0.8671i31.0000 - 0.0034i-0.0000 - 0.0034i-0.0000 - 0.0040i0.9999 - 0

13、.0107i-0.5000 - 0.8666i-0.5000 + 0.8654i40.8666 + 0.4990i-0.0006 - 0.0010i00.8660 + 0.4980i0.0000 - 1.0000i-0.8660 + 0.5020i50.8665 + 0.4991i-0.0005 - 0.0009i00.8660 + 0.4983i0.0000 - 1.0000i-0.8660 + 0.5017il 各支路三序电流、三相电流：支路I1I2I0IaIbIc1-2-0.0028 + 0.0000i-0.0028 + 0.0000i-0.0014 + 0.0000i-0.0071 +

14、 0.0001i0.0015 - 0.0000i0.0015 - 0.0000i1-30.0152 - 0.0002i0.0152 - 0.0002i0.0160 - 0.0002i0.0465 - 0.0005i0.0007 - 0.0000i0.0007 - 0.0000i2-30.0181 - 0.0002i0.0181 - 0.0002i0.0174 - 0.0002i0.0535 - 0.0006i-0.0007 + 0.0000i-0.0007 + 0.0000i(3) 节点3发生b、c两相短路故障l 故障点三序电流、三相电流：I1= 0 - 4.9697i I2= 0 + 4.9

15、697i I3= 0 ; Ia= 0 Ib= -8.6078 - 0.0000i Ic= 8.6078 + 0.0000i；l 各节点三序电压、三相电压：节点U1U2U0UaUbUc10.72730.272701.0000-0.5000 - 0.3937i-0.5000 + 0.3937i20.76970.230301.0000-0.5000 - 0.4671i-0.5000 + 0.4671i30.50000.500001.0000-0.5000-0.500040.7139 + 0.4122i0.1521 - 0.0878i00.8660 + 0.3244i0.0000 - 0.6487i-

16、0.8660 + 0.3244i50.7376 + 0.4259i0.1284 - 0.0741i00.8660 + 0.3518i0.0000 - 0.7035i-0.8660 + 0.3518il 各支路三序电流、三相电流：支路I1I2I0IaIbIc1-20 + 0.4240i0 - 0.4240i000.7343 -0.7343 1-30 - 2.2729i0 + 2.2729i00-3.9367 3.9367 2-30 - 2.6968i0 + 2.6968i00-4.67114.6711(4) 节点3发生b、c两相短路接地故障l 故障点三序电流、三相电流：I1= 0 - 6.447

17、3i I2= 0 + 3.4921i I3= 0 + 2.9552i ; Ia= 0 Ib= -8.6078 + 4.4327i Ic= 8.6078 + 4.4327i ；l 各节点三序电压、三相电压：节点U1U2U0UaUbUc10.64620.19160.06800.9058-0.3509 - 0.3937i-0.3509 + 0.3937i20.70120.16180.04370.9067-0.3879 - 0.4671i-0.3879 + 0.4671i30.35130.35130.35131.05400040.6687 + 0.3861i0.1069 - 0.0617i00.775

18、6 + 0.3244i0.0000 - 0.6487i-0.7756 + 0.3244i50.6995 + 0.4038i0.0902 - 0.0521i00.7897 + 0.3518i0.0000 - 0.7035i-0.7897 + 0.3518il 各支路三序电流、三相电流：支路I1I2I0IaIbIc1-20 + 0.5500i0 - 0.2979i0 - 0.1217i0 + 0.1304i0.7343 - 0.2477i-0.7343 - 0.2477i1-30 - 2.9486i0 + 1.5971i0 + 1.4167i0 + 0.0652i-3.9367 + 2.0925i

19、3.9367 + 2.0925i2-30 - 3.4987i0 + 1.8950i0 + 1.5384i0 - 0.0652i-4.6711 + 2.3402i4.6711 + 2.3402i3静态稳定性计算过程流程图和计算结果静态稳定性计算过程流程图如下：是否开始结束赋初值（delta,omega,步长h,时段长度Duration等）等输入故障切除时间CutTimei=0i<=round(CutTime/h)求时段初变化率求时段末估计值求时段末变化率求平均变化率求时段末计算值i=i+1显示此时功角同理求故障切除后的功角变化过程画出摇摆曲线计算结果：Ke=0.,delta=91,P=1.

20、,Eq=1.Ke=0.,delta=92,P=1.,Eq=2.Ke=0.,delta=93,P=1.,Eq=2.Ke=0.,delta=95,P=1.,Eq=2.Ke=0.,delta=96,P=1.,Eq=2.Ke=0.,delta=97,P=1.,Eq=2.Ke=0.,delta=98,P=1.,Eq=2.Ke=0.,delta=99,P=1.,Eq=2.Ke=0.,delta=101,P=1.,Eq=2.Ke=1.,delta=102,P=1.,Eq=2.Ke=1.,delta=103,P=1.,Eq=2.Ke=1.,delta=104,P=1.,Eq=2.Ke=1.,delta=105

21、,P=1.,Eq=2.Ke=1.,delta=106,P=1.,Eq=2.Ke=1.,delta=107,P=1.,Eq=2.Ke=1.,delta=108,P=1.,Eq=2.Ke=1.,delta=109,P=1.,Eq=2.Ke=1.,delta=108,P=1.,Eq=2.Ke=1.,delta=105,P=1.,Eq=2.Ke=2.,delta=103,P=1.,Eq=2.Ke=2.,delta=100,P=1.,Eq=2.Ke=2.,delta=98,P=1.,Eq=2.Ke=2.,delta=96,P=1.,Eq=2.Ke=2.,delta=94,P=1.,Eq=2.Ke=2.,

22、delta=92,P=1.,Eq=2.Ke=2.,delta=90,P=1.,Eq=2.Ke=2.,delta=88,P=1.,Eq=2.Ke=2.,delta=86,P=1.,Eq=2.Ke=2.,delta=84,P=1.,Eq=2.Ke=3.,delta=82,P=1.,Eq=2.Ke=3.,delta=81,P=1.,Eq=2.Ke=3.,delta=79,P=1.,Eq=2.Ke=3.,delta=77,P=1.,Eq=2.Ke=3.,delta=76,P=1.,Eq=2.Ke=3.,delta=74,P=1.,Eq=2.Ke=3.,delta=73,P=1.,Eq=2.Ke=3.,

23、delta=71,P=1.,Eq=2.Ke=3.,delta=70,P=1.,Eq=2.Ke=3.,delta=69,P=1.,Eq=2.Ke=4.,delta=67,P=1.,Eq=2.Ke=4.,delta=66,P=1.,Eq=2.Ke=4.,delta=65,P=1.,Eq=2.Ke=4.,delta=63,P=1.,Eq=2.Ke=4.,delta=62,P=1.,Eq=2.Ke=4.,delta=61,P=1.,Eq=2.Ke=4.,delta=60,P=1.,Eq=2.Ke=4.,delta=59,P=1.,Eq=2.Ke=4.,delta=58,P=1.,Eq=2.Ke=4.,

24、delta=56,P=1.,Eq=2.Ke=5.,delta=55,P=1.,Eq=2.Ke=5.,delta=54,P=1.,Eq=2.Ke=5.,delta=53,P=1.,Eq=2.Ke=5.,delta=52,P=1.,Eq=2.Ke=5.,delta=51,P=1.,Eq=2.Ke=5.,delta=50,P=1.,Eq=1.Ke=5.,delta=49,P=1.,Eq=1.Ke=5.,delta=49,P=1.,Eq=1.最终选择放大倍数Ke=2.24暂态稳定性计算过程流程图和计算结果 暂态稳定性计算过程流程图如下：否否否是是是开始清屏设初值：Xd_all,Xdd_all,Xq,X

25、q_all,Eq0,UG0,Xe,Xd,XddKe=0.1Ke_result=Ke,delta_result=delta, K1_result= K1,Kemax_result=KemaxP_result=Eq*U*sin(delta*pi/180)/Xd_all,delta=delta+0.01delta=49.00求Kemax,Kemin求Eq,Eqq求UGd,UGq,UG求K1K6Ke<5.7delta<=110.00(K1>0)&(Ke>=Kemin)&(Ke<=Kemax)结束打印Ke_result，delta_result, K1_re

26、sult,Kemax_result计算结果如下： 0.15s时切除故障的摇摆曲线 0.25s时切除故障的摇摆曲线根据摇摆曲线判断：0.15s时切除故障系统暂态稳定，0.25s时切除故障系统的功角无限增大，系统失去稳定性。另据试探，0.2s时系统临界稳定，为保证系统暂态稳定，切除时间最大值约为0.2s。程序：Matrix导纳阵计算程序clc;clear;%-输入已知条件-%bus_Num1=3; %节点数bus_Num2=5; %包括发电机节点的节点数branch_Num1=3; %线路数branch_Num2=5; %包括发电机支路的支路数branch1_No1=1,1,2;%支路首节点bra

27、nch1_No2=2,3,3;%支路末节点branch2_No1=1,1,2,1,2;branch2_No2=2,3,3,4,5;%输入支路各序阻抗，z1_branch表示支路正序阻抗，z2_branch表示支路负序阻抗，z0_branch表示支路零序阻抗z1_branch(1)=j*0.1;z1_branch(2)=j*0.1;z1_branch(3)=j*0.1;z2_branch=z1_branch;z0_branch(1)=j*0.2;z0_branch(2)=j*0.2;z0_branch(3)=j*0.2;%-第一步：不考虑发电机节点计算节点导纳矩阵-%节点导纳矩阵，Y1表示不计发

28、电机节点的正序网络节点导纳阵，Y2表示不计发电机节点的负序网络节点导纳阵，Y0表示不计发电机节点的零序网络节点导纳阵，Y1=zeros(bus_Num1);Y1(1,1)=1/(j*0.15);Y1(2,2)=1/(j*0.075);%请同学们求正序节点导纳矩阵for i=1:branch_Num1 y=1/z1_branch(i); ii=branch1_No1(i); jj=branch1_No2(i); Y1(ii,ii)=Y1(ii,ii)+y; Y1(ii,jj)=Y1(ii,jj)-y; Y1(jj,ii)=Y1(jj,ii)-y; Y1(jj,jj)=Y1(jj,jj)+y;en

29、dY2=Y1;%负序等于正序Y0=zeros(bus_Num1);Y0(1,1)=1/(j*0.0525);Y0(2,2)=1/(j*0.0263);%请同学们求零序节点导纳矩阵for i=1:branch_Num1 y=1/z0_branch(i); ii=branch1_No1(i); jj=branch1_No2(i); Y0(ii,ii)=Y0(ii,ii)+y; Y0(ii,jj)=Y0(ii,jj)-y; Y0(jj,ii)=Y0(jj,ii)-y; Y0(jj,jj)=Y0(jj,jj)+y;end %-第二步：考虑发电机节点计算节点导纳矩阵-%节点导纳矩阵，YY1表示计及发电机

30、节点的正序网络节点导纳阵，YY2表示计及发电机节点的负序网络节点导纳阵z1_branch(4)=j*0.05;z1_branch(5)=j*0.025;z2_branch=z1_branch;YY1=zeros(bus_Num1);YY1(4,4)=1/(j*0.1);YY1(5,5)=1/(j*0.05);%请同学们求正序节点导纳矩阵for i=1:branch_Num2 y=1/z1_branch(i); ii=branch2_No1(i); jj=branch2_No2(i); YY1(ii,ii)=YY1(ii,ii)+y; YY1(ii,jj)=YY1(ii,jj)-y; YY1(j

31、j,ii)=YY1(jj,ii)-y; YY1(jj,jj)=YY1(jj,jj)+y;endYY2=YY1;%负序等于正序%-第三步：计算节点阻抗矩阵-%Z1=inv(Y1);Z2=inv(Y2);Z0=inv(Y0);ZZ1=inv(YY1);ZZ2=inv(YY2);短路电流计算程序shortcircuitclc;clear;%数据来源于教材电力系统暂态分析P77例（3-4），P143例（5-7）Matrix %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量Fault_Node=input('输入短路点编号 ;n Fault_Node=');Fault_Type=input

32、('输入短路类型 ;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 1 a a2 1;%-第一步：计算短路点的序电流，相电流-%计根据故障类型选择不同的计算公式，计算故障点各序电流if Fault_Type=0 %三相短路 I_Fau

33、lt1=1/Z1(Fault_Node,Fault_Node); I_Fault2=0; I_Fault0=0;elseif Fault_Type=1 %a相接地短路 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); I_Fault2=I_Fault1; I_Fault0=I_Fault1; elseif Fault_Type=2 %a相经10欧姆电阻接地短路 Zf=10 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault

34、_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+3*Zf); I_Fault2=I_Fault1; I_Fault0=I_Fault1; elseif Fault_Type=3 %bc两相相间短路 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node); I_Fault2=-1*I_Fault1; I_Fault0=0; elseif Fault_Type=4 %bc两相短路接地 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,

35、Fault_Node)*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); I_Fault2=-1*I_Fault1*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); I_Fault0=-1*I_Fault1*Z2(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); end;st

36、r='短路电流'strIabc=T*I_Fault1 I_Fault2 I_Fault0.' %相量Iabc_effective=abs(Iabc) %有效值%-第二步：计算各个节点的序电压，相电压-%故障分量 Ifault_node1=zeros(bus_Num1,1);Ifault_node2=zeros(bus_Num1,1);Ifault_node0=zeros(bus_Num1,1);for m=1:bus_Num1 if m=Fault_Node Ifault_node1(m,1)=-1*I_Fault1; %故障电流分量 Ifault_node2(m,1)

37、=-1*I_Fault2; Ifault_node0(m,1)=-1*I_Fault0; else Ifault_node1(m,1)=0; Ifault_node2(m,1)=0; Ifault_node0(m,1)=0; endendUfault_node1=zeros(bus_Num1,1);Ufault_node2=zeros(bus_Num1,1);Ufault_node0=zeros(bus_Num1,1);Ufault_node1=Y1Ifault_node1; %故障电压分量Ufault_node2=Y2Ifault_node2;Ufault_node0=Y0Ifault_no

38、de0;%正常分量 Unormal_node1=ones(bus_Num1,1);Unormal_node2=zeros(bus_Num1,1);Unormal_node0=zeros(bus_Num1,1);%各个节点的序电压Uall_node1=zeros(bus_Num1,1);Uall_node2=zeros(bus_Num1,1);Uall_node0=zeros(bus_Num1,1);Uall_node1=Unormal_node1+Ufault_node1; %叠加原理Uall_node2=Unormal_node2+Ufault_node2;Uall_node0=Unorma

39、l_node0+Ufault_node0;str='各节点序电压'strU120=Uall_node1 Uall_node2 Uall_node0.' % 序分量str='各节点相电压'strUabc=T*Uall_node1 Uall_node2 Uall_node0' %相量('表示求转置矩阵）Uabc_effective=abs(Uabc) %有效值%-第三步：计算各个支路的序电流，相电流-%Ibranch_1=zeros(branch_Num1,1);Ibranch_2=zeros(branch_Num1,1);Ibranch_0

40、=zeros(branch_Num1,1);for i=1:branch_Num1 Ibranch_1(i,1)=(Uall_node1(branch1_No1(i),1)-Uall_node1(branch1_No2(i),1)/z1_branch(i); %支路电压/支路阻抗 Ibranch_2(i,1)=(Uall_node2(branch1_No1(i),1)-Uall_node2(branch1_No2(i),1)/z2_branch(i); Ibranch_0(i,1)=(Uall_node0(branch1_No1(i),1)-Uall_node0(branch1_No2(i),

41、1)/z0_branch(i);endstr='各支路序电流'strIbranch_120=Ibranch_1 Ibranch_2 Ibranch_0.' % 序分量str='各支路相电流'strIbranch_abc=T*Ibranch_1 Ibranch_2 Ibranch_0' %相量('表示求转置矩阵）Ibranch_abc_effective=abs(Ibranch_abc) %有效值%-第四步：计算发电机节点的序电压，相电压-%故障分量 I2fault_node1=zeros(bus_Num2,1);I2fault_node2

42、=zeros(bus_Num2,1);I2fault_node0=zeros(bus_Num2,1);for m=1:bus_Num2 if m=Fault_Node I2fault_node1(m,1)=-1*I_Fault1; I2fault_node2(m,1)=-1*I_Fault2; I2fault_node0(m,1)=-1*I_Fault0; else Ifault_node1(m,1)=0; Ifault_node2(m,1)=0; Ifault_node0(m,1)=0; endendU2fault_node1=zeros(bus_Num2,1);U2fault_node2=

43、zeros(bus_Num2,1);U2fault_node0=zeros(bus_Num1,1);U2fault_node1=YY1I2fault_node1;U2fault_node2=YY2I2fault_node2;%正常分量 U2normal_node1=ones(bus_Num2,1);U2normal_node2=zeros(bus_Num2,1);U2normal_node0=zeros(bus_Num2,1);%各个节点的序电压U2all_node1=zeros(bus_Num2,1);U2all_node2=zeros(bus_Num2,1);U2all_node0=zer

44、os(bus_Num2,1);U2all_node1=U2normal_node1+U2fault_node1;U2all_node2=U2normal_node2+U2fault_node2;U2120=U2all_node1 U2all_node2 U2all_node0' % 序分量U2abc=T*U2all_node1 U2all_node2 U2all_node0' ; %相量('表示求转置矩阵）U2abc_effective=abs(U2abc); %有效值str='发电机节点序电压'strU4120=U2120(:,4)U5120=U212

45、0(:,5)str='发电机节点相电压相量及有效值'strU4abc=U2abc(:,4)U5abc=U2abc(:,5)U4abc=U2abc_effective(:,4)U5abc=U2abc_effective(:,5)%-第五步：计算发电机支路的序电流，相电流-%Ibranch4_1=Ibranch_1(1,1)+Ibranch_1(2,1);Ibranch4_2=Ibranch_2(1,1)+Ibranch_2(2,1);Ibranch4_0=Ibranch_0(1,1)+Ibranch_0(2,1);Ibranch5_1=-Ibranch_1(1,1)+Ibranc

46、h_1(3,1);Ibranch5_2=-Ibranch_2(1,1)+Ibranch_2(3,1);Ibranch5_0=-Ibranch_0(1,1)+Ibranch_0(3,1);str='发电机支路序电流'strIbranch4_120=Ibranch4_1 Ibranch4_2 Ibranch4_0' % 序分量Ibranch5_120=2Ibranch5_1 Ibranch5_2 Ibranch5_0' % 序分量str='发电机支路相电流'strIbranch4_abc=T*Ibranch4_1 Ibranch4_2 Ibranch

47、4_0' %相量('表示求转置矩阵）Ibranch4_abc_effective=abs(Ibranch4_abc) %有效值Ibranch5_abc=2T*Ibranch5_1 Ibranch5_2 Ibranch5_0' %相量('表示求转置矩阵）Ibranch5_abc_effective=abs(Ibranch5_abc) %有效值静态稳定计算stability_smallsignalclear;clc;U=1; %系统电压Eq0=1.972; %空载电动势UG0=1.21; %机端电压Xe=0.504; %线路电抗值Xd=0.982; %同步电抗Xdd

48、=0.344; %暂态电抗Xd_all=1.486; %系统电抗(包含同步电抗)Xdd_all=0.848; %系统电抗(包含暂态电抗)Tj=10; %惯性时间常数Td0=10; %励磁绕组时间常数Ke=0.1;while(Ke<5.7) for delta=49:110 %请同学们求各个参数 Eq1=Xdd_all*Eq0/Xd_all+U*(1-Xdd_all/Xd_all)*cos(delta*pi/180); UGd0=U*sin(delta*pi/180)*Xd/Xd_all; UGq0=Eq0*Xe/Xd_all+Xd*U*cos(delta*pi/180)/Xd_all; UG=sqrt(UGd02+UGq02); a=Ke2*Xe2/(Xd_all2)-1; b=2*(Ke2*Xd*Xe*U*cos(delta*pi/180)/(Xd_all2)+(Eq0+Ke*UG0); c=Ke2*Xd2*U2*cos(delta)*cos(delta)/(Xd_all2)+Ke2*Xd2*U2*sin(delta)*sin(delta)/(Xd_all2)-(Eq0+Ke*UG0)2; Eq=(-b+sqrt(b*b-4*a*c)/(2*a);%求解方程 K1=Eq1*U/Xdd_all*cos