电力系统课程设计报告

1、东南大学电气工程学院2015年电力系统综合课程设计报告班级160124姓名齐济学号16012429指导老师陈中小组成员王武森石应隆目录线路数据 2一、网络拓扑 2二、发电机稳态数据 2三、变电站数据 2四、线路数据 3五、发电机暂态数据 3典型方式电力系统潮流计算 4一、发电厂变压器选择及参数计算 4二、形成导纳矩阵 4三、求出系数矩阵B'和B' ' 5四、设置各节点电压初值 5五、求功率不平衡量 5六、解修正方程 5七、求新的初值 6八、判断是否已收敛 6九、计算结果 6变电站变压器选择和主接线设计 7一、变压器选择 71. 变电站主变压器 72. 发电厂变压器 11

2、3. 变压器等效电路 11二、主接线设计 121. 发电厂主接线 122. 变电站主接线 13三、总结与提升 16线路断路器选择 18仿真曲线计算 20感想与建议 23附录一： PQ 解耦法潮流计算 Matlab程序 24附录二：改进欧拉法计算功角特性 Matlab 程序 25詩负荷2T1T228793G2G1网络拓扑线路数据负荷1T3G3负荷341发电机稳态数据名称额定电压机端电压典型方式输出有功G110KV1.05平衡机G210KV1.051.60G310KV1.050.90三、变电站数据负荷额定电压为10KV变电站编号负荷1负荷2负荷3最大负荷（100MVA1.451.11.2典型方式负

3、荷（100MVA1.250.91.0最小负荷（100MVA1.050.60.7Tmax（小时）500040004500功率因数cos ©0.850.80.85低压侧电压（kV）101010对备用要求30%40%30%对调压要求逆调压顺调压常调压负荷重要性重要重要一般四、线路数据名称电压电阻电抗对地导纳（B）7-8110KV0.010.0850.0448-9110KV0.020.0920.0407-5110KV0.030.1020.0765-4110KV0.010.1620.0899-6110KV0.020.1750.0376-4110KV0.030.1830.042负序和正序相等、零

4、序为正序的 3倍五、发电机暂态数据名称D轴电抗D轴暂态电抗阻尼系数Tj （惯性时间常数）G10.140.072.047.2G20.900.221.912.8G21.320.132.16.02典型方式电力系统潮流计算发电厂变压器选择及参数计算三个发电厂的变压器选择如下：T1:单台 SFP7-120000/110T2:双台 SFP7-90000/110T3:单台 SFP7-90000/110经计算得到三台变压器的阻抗、导纳参数为:形成导纳矩阵名称电阻R电抗X对地电导G对地电纳BT10.00290.08750.001060.006T20.00210.05830.00170.0108T30.00420

5、.1170.000850.0054根据线路数据和发电机变压器数据，可以算出电力系统网络中各节点的自导纳和各节点之间的互导纳。由此可得电力系统网络的节点导纳矩阵为厂-8.545i008.54i00000-11.429i000011.423i000-17.16i000008.54i00-19.945i6.15i5.32i000006.15i-15.09i09.02i00005.32i0-10.92i00011.423i009.02i0-31.983i11.6i00000011.6i-21.983i0017.15i005.64i010.3 8iJ12 / 28三、求出系数矩阵B和B节点2作为平衡节点

6、，节点1、3为PV节点，其余节点为PQ节点，由系数B =厂f-8.54508.54000000-17.160000017.158.540-19.9456.155.32000006.15-15.0909.0200005.320-10.92005.640009.020-31.98311.600000011.6-21.98310.83017.15005.64010.38-33.13矩阵的定义可得B''=19.9456.155.320006.15-15.0909.02005.320-10.92005.6409.020-31.98311.6000011.6-21.98310.83005.

7、64010.38-33.13四、设置各节点电压初值五、求功率不平衡量由公式:可求得有功不平衡量和无功不平衡量，从而得到和六、解修正方程求解矩阵方程得出各节点电压相角的修正量和各节点电压模值的修正量七、求新的初值求出修正后的值八、判断是否已收敛由收敛判断依据和（其中&取0.00001）判断是否收敛，若不满足收敛要求，则回到第五步开始下一次迭代，若满足收敛 要求，则退出循环。九、计算结果用Matlab编写PQ解耦法的程序代码见附录，潮流计算的结果如下平衡节点的功率PV节点的注入无功功率程序见附录。变电站变压器选择和主接线设计变压器选择1. 变电站主变压器1.1 变压器选择负荷1:因负荷1为

8、重要负载，故采用两台主变压器为其供电，实现双路供电，提高供 电可靠性。考虑到变压器容量应能满足最大负荷需求，故主变压器所需总容量最少为同时考虑到备用需要，则主变压器所需总容量最少为根据相关设计原则变电所的两台主变压器总安装容量可取为根据以上数据，在 110KV变压器型号表2中选取可用型号为沈阳变压器厂的 SFP7-120000/110型变压器，两台并联使用。根据相关设计原则，变电所其中一台变压器停运后，剩余变压器能承担全部负荷 的60%以上。在当前方案下，一台停运后，剩余变压器能承担负荷百分比为 ，满足设计要求。实际上，考虑变压器的过负荷能力为，变压器的容量远远满足设计要求。负荷2：因负荷2为

9、重要负载，故采用两台主变压器为其供电，实现双路供电，提高供 电可靠性。考虑到变压器容量应能满足最大负荷需求，故主变压器所需总容量最少为同时考虑到备用需要，则主变压器所需总容量最少为根据相关设计原则1，变电所的两台主变压器总安装容量可取为根据以上数据，在 110KV变压器型号表2中选取可用型号为沈阳变压器厂的SFP7-90000/110型变压器，两台并联使用。在当前方案下，一台停运后，剩余变压器能承担负荷百分比为一，满足设计要求。负荷3:因负荷1为一般负载，故采用单台主变压器为其供电。考虑到变压器容量应能满足最大负荷需求，故主变压器所需总容量最少为同时考虑到备用需要，则主变压器所需总容量最少为根

10、据相关设计原则，变电所的两台主变压器总安装容量可取为根据以上数据，在 110KV变压器型号表2中选取可用型号为沈阳变压器厂的SFP7-120000/110 型变压器。在当前方案下，考虑变压器的过负荷能力为1,变压器可承载的最大负荷为，容量满足要求。另考虑到可以进行需求侧管理，变压器满足设计要求。1.2 变压器参数计算SFP7-120000/110（ YN，d11）根据变压器型号表2可知：归算到一次侧，电阻电抗 电导电纳 SFP7-90000/110（ YN，d11）根据变压器型号表2可知：归算到一次侧，电阻电抗 电导电纳1.3 变压器分接头选择变电站变压器分接头选择需考虑负载对调压的要求以及母

11、线电压随负荷变化的电压变化值。经过潮流计算，分别得到每一负载最大最小的情况下相应母线的电压值：最小负荷最大负荷典型负荷负荷1104.71 kv100.39 kv102.64 kv负荷2108.89 kv105.53 kv106.92 kv负荷3106.45 kv100.66 kv103.13 kv负荷1：归算至高压侧的二次侧电压最小负荷时变压器阻抗上的电压降为归算至高压侧的二次侧电压由于二次侧电压要求满足逆调压要求，即最大负荷时，低压侧电压不得低于10.5KV分接头选择为最小负荷时，低压侧电压不得高于10KV。分接头选择为取平均值:选择最接近的分接头负荷2：归算至高压侧的二次侧电压最小负荷时变

12、压器阻抗上的电压降为归算至高压侧的二次侧电压由于二次侧电压要求满足逆调压要求，即最大负荷时，低压侧电压不得低于10.25KV分接头选择为最小负荷时，低压侧电压不得高于10.75KV。分接头选择为取平均值：选择最接近的分接头负荷3：归算至高压侧的二次侧电压最小负荷时变压器阻抗上的电压降为10.2KV归算至高压侧的二次侧电压由于二次侧电压要求满足逆调压要求，即最大负荷时，低压侧电压约为分接头选择为最小负荷时，低压侧电压约为10.5KV。分接头选择为取平均值:选择最接近的分接头2. 发电厂变压器发电厂与电力系统经变压器相连的目的一般是使发电机发出的功率能送入电力系 统1。因此根据发电机稳态数据表可选

13、择合适容量的升压变压器。名称额定电压机端电压典型方式输出有功G110KV1.05平衡机G210KV1.051.60G310KV1.050.90参考110KV变压器型号表 ，综合考虑各发电机输出功率及变压器的过负荷能力 后，选择如下：G1:单台 SFP7-120000/110G2 :双台 SFP7-90000/110G3 :单台 SFP7-90000/110发电机升压变压器分接头选择均为121/10.5。变压器参数计算见1.2节。3. 变压器等效电路图2线路阻抗没有在图中绘出。主接线设计1. 发电厂主接线1.1发电厂主接线设计110KV电压级，出线回路 3回。为保证出现断路器检修时不停电，应采用

14、单母线 分段带旁路母线或双母线带旁路接线，以保证其供电的可靠性和灵活性。但考虑到发 电机单机容量很大，远大于有关设计规程4对选用单母线分段接线每段上不宜超过12MW的规定，且双母线的可靠性和灵活性大幅优于另一者，对重要用电负荷有更好 的保证，故确定为双母线带旁路母线接线方式。负荷3110KV负荷2 补充：主接线设计还应考虑经济费用、可靠性、灵活性等方面的比较。但因无法估 算其施工总投资和年运行费用等数据，在此略去其比较过程。1.2 电网的电能损耗 5根据三个负荷的最大负荷利用小时数 、 、 ，以及每个负荷的功率因数 、 、 ，可计算得到三段线路平均的 和 :根据最大负荷损耗小时数 与的关系表

15、5可知：因此可以计算得到线路的总电能损耗为2. 变电站主接线2.1 主接线设计 变电站的主接线应满足供电可靠，运行灵活，操作检修方便、节约投资和便 于扩建等要求。负荷 1、负荷 2：变电站10KV负荷125MVA (90MVA),变压器为2台120MVA(90MV A)。考 虑到桥型接线可靠性不高，只适用于小容量变电站；多角形接线可以节省费用，并 具有双断路器双母线接线的可靠性6，但不利于后期扩建。根据相关资料 7，当变电站装有两台变压器时， 610KV 侧宜采用分段单母线。 线路为 12 回以上时亦可采用双母线，当不允许停电检修断路器时，可设置旁路母 线。考虑到该负荷为重要用电负荷，因此确定

16、负荷1、负荷 2 的主接线为分段单母线带旁路母线的设计。1-lb4WPJ QSp1QSQSTDT.1QSd(QFp)QSQS-10KV丁11 .11WB图430 / 28上图所示的主接线增设了一组旁路母线WP及各出线回路中相应的旁路隔离开关QSp,分段断路器 QSd兼作旁路断路器 QFp，并设有分段隔离开关 QSd。平时旁路母线不带电，QS1、QS2及QFp合闸，QS3、QS4及QSd断开，主接线系统按单母线分段方式运行。当需要检修某一出线断路器时，可通过闸操作，由分段断路器代替旁路断路器，使旁路断路器经QS4、QFP、QS1接至1段母线，或经QS2、QFP、QS3接至2段母线而带电运行，并经

17、过被检修断路器所在回路的 旁路隔离开关及其两侧的隔离开关进行检修，而不中断其所在线路的供电。此时， 两段工作母线既可通过分段隔离开关QSd并列运行也可分列运行。所以，这种接线方式具有相当高的可靠性和灵活性。但是考虑到Qsd需要检修，且负载的重要性，故采用独立的母线分段断路器。为限制发电厂内部故障和出线故障时的短路电流，以便选用轻型的断路器，故在母线分段断路器上和出线上分别串接母线电抗器和线路电抗器。因此将主接线图修改如下：ULJQSpQFpWP10KVWB负荷3：该负载为一般负载，故可选用成本较低的的简单接线方案。根据以上分析,我们选择无旁路母线的分段单母线接线方式，提高了供电的可靠性和灵活性

18、。并为以后的扩建打下基础。*QLi,QL3.QL4：!丨）III10KV图62.2主变压器的电能损耗根据相关资料，令为电能损耗，则双线圈变压器的电能损耗可以用下式计算:式中一一相同的变压器（并列运行）的台数；每台变压器的额定容量（ KV A ）;n台变压器的总负荷（KV A）;负荷S的使用时间小时数（h）;、一一每台变压器的空载有功损耗（ KW）,无功损耗（Kvar）、一一每台变压器的短路有功损耗（ KW）,短路损耗（Kvar）无功当量，通常取 0.10.15电能损耗（KW h）带入相关参数，可以计算得到：负荷1:负荷2 :负荷3 :三、总结与提升报告确定了每个负载的变压器和其主接线，同时计算

19、出变压器的相应参数以及主接 线的年电能损耗，基本完成了设计要求。在设计过程中参考各种文献等，有效的锻炼了 学生资料查阅能力及自学能力，进一步融会贯通课堂上所学的内容，加深了学生对相关 知识的理解。但该设计仍存在一定的提升空间。经过变电站变压器的分接头的检验可以看出，当前选择的分接头无法很好的满足负 载的调压要求。因此需要更多的电压调整方式对其进行进一步补偿。进一步电压调整方 法有：1发电机调压2改变电力网无功功率分布 3改变输电线路参数进行调压。对于重要负载1和重要负载2，可以增加变压器台数，采取备用变压器的接线图, 但或许会进一步增加经费。如下图表示具有三台工作变压器与一台本备用变压器的变电

20、站。正常时高压母线分 段断路器QF1断开，当任一进线停电时分段断路器均可以自动闭合。三台工作变压器正 常时均经过其高低压侧断路器在各自的分段上工作。备用变压器低压方向经过不同的低 压断路器接到三个低压母线上，这些断路器平时都断开。这种电路可以保证可靠地向重 要用户供电。F2110KV图7参考文献：电能系统基础P207-20822电力系统课程设计及毕业设计参考资料P108电力系统课程设计及毕业设计参考资料P33发电厂电气部分 P106电能系统基础 P108 发电厂电气部分 P115-116电气工程专业毕业设计指南电力系统分册P67线路断路器选择1. 断 路器种类和型式的选择考虑安装调试和运行维护

21、的方便，可考虑少油式断路器，其中的油只做灭弧和触头间弧 隙的绝缘介质，断路器的带电导体与接地部件之间的绝缘主要采用瓷件，油量少、占地少、 价廉，已有长期的运行经验，多用在 110220KV电压等级配电装置中。2. 额 定电压和额定电流的选择高压断路器的额定电压和额定电流的选择应满足：其中Un和Usn分别为断路器和电网的的额定电压；In和Imax分别为断路器的额定电流和电网的最大负荷电流。O 1有原始资料知电网 5-7线路的额定电压为110kv，即UsN=110kv,Un> 110kv。Imax为最大持续工作电流，将负荷1、2、3都设为最大值，由PQ解耦法潮流计算的程序计算出各节点电压（黑

22、色字体为相量表示）Ui=1.050 / -7.789U2=1.050 / 0U3=1.050 / -6.031U4=0.968 / -13.740U5=0.898 / -15.841U6=0.904 / -17.849U7=0.975 / -9.716U8=0.949 / -12.492U9=0.995 / -11.154由导纳矩阵得 Y 57=2.654-9.02il max =(U5 - U7 )?Y57计算得 lmax= 1.1866Z -213.924 ,则 Imax=1.1866。电流基准值 Ib=Sb/ Ub =100MVA/ X110KA =524.86A。所以有 In >

23、 622.8A。3. 开 断电流的选择高压断路器的额定开断电流INbr 是指在额定电压下能保证正常开断的最大短路电流，是表征高压断路器的开断能力的重要参数。本设计中的变电站高压侧为110kv，属于中型变电站，一般采用的是中速断路器，开断时间较长（计非周期分量的影响，采用起始次暂态电流0.1s）短路电流中的非周期分量衰减较多，可不 |"校验，即I也是实际开断瞬间短路电流周期分量有效值。|"的计算：由附录程序，设置短路点为5，计算出短路时5-7之间的电流I” = -22.71 -4.76i,有效值|" =23.20。由此计算出次暂态电流的有名值l" =23.

24、20 X 524.86A=12.177KA，所以iNbr > 12.177KA。4. 短路关合电流的选择若在断路器合闸时已经存在短路故障，则在合闸过程中动静触头在未接触时就有很大的 短路电流流过（预击穿），更容易损坏断路器。且断路器在关合短路电流时不可避免的在接通 后又自动跳闸，此时还需要能切断短路电流。为了保证断路器在关合短路时的安全，断路器 的额定关合电流Inci不应小于短路电流最大冲击值Ish，即Inci > Ish。在高压线路进出线中：Ish =2.55 l" （ l"是起始次暂态电流）在发电机端及母线侧：Ish =2.69 l"所以：Inci

25、 > 2.55X 12.177KA=31.05KA由上述讨论可选择 SW3-110G/1200少油户外断路器。额定电压（KV）额定电流（A）额定开断 电流（KA ）短路关合电流（KA）热稳定电 流（KA ）固有分闸时间（s）合闸时间（s）计算数据110622.812.仃731.05SW3-110G/1200110120015.84115.8(4s)0.070.45. 短路热稳定校验短路电流通过电气设备时，电气设备发热效应不应超过允许值，满足热稳定的条件为：2it t > Qk式中Qk为短路电流产生的热效应；it , t分别为电气设备允许通过的热稳定电流的和时间。由于数据不全，无法校

26、验。6. 短路动稳定校验电气设备承受短路电流机械效应的能力称为动稳定，满足条件为les > ish 其中les为允许通过的动稳定电流的有效值。由于数据不全，无法校验。仿真曲线计算1.由潮流计算结果求初值发电机的电动势可表示为求得各发电机电动势为各发电机的机械功率为2.计算负荷的等值导纳负荷的等值导纳可以表示为求得各负载的等值导纳为3. 修改网络节点的导纳矩阵在原潮流计算用网络导纳矩阵的基础上形成一个包含负荷等值导纳以及 增加发电机电动势节点的导纳矩阵。新增加发电机电动势节点的自导纳为，它们只和相应的发电机端电压节点之间有互导纳，而发电机端电压节点的自导纳也要相应地增加。同时，由于线路5-

27、7发生三相短路故障，相当于在线路5-7上并联一个对地的附加电抗 ，且 值为 0,因此得到修改后的网络节点导纳矩阵为其中A为3*3的方阵，E为3*9的矩阵，C为9*3的矩阵，D为9*9的方阵。经过网络化简可得三相短路故障发生时三台发电机之间的导纳矩阵为得到三台发电机之间的导纳矩阵为2.287741 + 0.401924i0.195099 + 0.2591310.580867 + 0.854405i0.195099 + 0.2591310.297002 - 4.931720i0.133391 + 0.397078i0.580867 + 0.854405i0.133391 + 0.397078i0.

28、363343 - 2.312614i4. 计算发电机功角和电磁功率 三台发电机的初始功角为由发电机功率表达式可计算出各台发电机的电磁功率。5. 求发电机的运动方程及功角特性求解发电机的运动方程，使用改进欧拉法对方程组进行迭代，求出三台发电机功角的变化、电磁功率的变化以及功角特性。功角特性曲线:感想与建议本次课程设计研究的是一个三机九节点的环网系统， 本系统相比课堂上 学习的内容而言更为复杂，不确定因素也更多，更多的是需要我们自行寻找 设计方案，并在我们设计的方案上进行进一步的计算和分析。在这次课程设计中， 我负责系统的潮流计算以及变压器的选择和主接线 的设计。曾经在电力系统稳态的学习过程中接触

29、过 matpower 计算系统的潮流， 但matpower是圭寸装好的程序，这次我们自己编写matlab中PQ解耦法的代码， 让我们更加深刻的理解了相关计算过程，加深了课程和实践的结合。在变压器和主接线的选择时， 主要是查找多个文献的相关设计要求和设 计方案，这有效的提高了我文献查阅和综合资料的能力，对以后的学习和研 究有极大的帮助。对于变压器、主接线相关参数的计算，让我对课本上以前 学习过的知识加深了理解，对曾经不那么重视的部分有了更清晰的了解。其 计算与潮流计算的结合，更是让我们所学的东西融会贯通。通过这次课程设计， 我学到了很多电力系统分析计算的相关知识， 对电 力系统潮流计算和电力系统

30、故障后的暂态分析有了更深一步的了解。* (-0.38+ 6.15+ 5.32+ 11.6* (-0.645dp8/u(8)附录一： PQ 解耦法潮流计算 Matlab 程序clearclcload('tide_B.mat')theta = 0 0 0 0 0 0 0 0 0' u = 1.05 1.05 1.05 1 1 1 1 1 1'dtheta = 0 0 0 0 0 0 0 0 0'du = 0 0 0 1 1 1 1 1 1'dd = 0.00001;dp_u = 1 1 1 1 1 1 1 1' udtheta = 1 1 1

31、 1 1 1 1 1'dq_u = 1 1 1 1 1 1' dutemp = 1 1 1 1 1 1'i = 0;while(abs(dtheta(1)>dd)|(abs(dtheta(9)>dd)|(abs(dtheta(3)>dd)|(abs(dtheta(4)>dd)|(abs(dtheta(5)>dd)|(abs(dtheta(6)>dd)|(abs(dth eta(7)>dd)|(abs(dtheta(8)>dd)|(abs(du(1)>dd)|(abs(du(9)>dd)|(a bs(du(3)&

32、gt;dd)|(abs(du(4)>dd)|(abs(du(5)>dd)|(abs(du(6)>dd)|(abs (du(7)>dd)|(abs(du(8)>dd)dp1 = 0.9 - 1.05 * u(4) * 8.54 * sind(theta(1) - theta(4);dp3 = 1.6 - 1.05 * 17.15 * u(9) * sind(theta(3) - theta(9);dp4 = 0 - u(4) * (1.05 * 8.54 * sind(theta(4) - theta(1) + u(5) * cosd(theta(4) - thet

33、a(5) + 6.15 * sind(theta(4) - theta(5) + u(6) * (-0.872 * cosd(theta(4) - theta(6) + 5.32 * sind(theta(4) - theta(6) + u(4) * 1.56);dp5 = -1.0625 - u(5) * (u(4) * (-0.38 * cosd(theta(5) - theta(4) * sind(theta(5) - theta(4) + u(7) * (-2.654 * cosd(theta(5) - theta(7)+ 9.02 * sind(theta(5) - theta(7)

34、 + u(5) * 3.034);dp6 = -0.85 - u(6) * (u(4) * (-0.872 * cosd(theta(6) - theta(4) * sind(theta(6) - theta(4) + u(9) * (-0.645 * cosd(theta(6) - theta(9) + 5.64 * sind(theta(6) - theta(9) + u(6) * 1.52);dp7 = 0 - u(7) * (11.423 * 1.05 * sind(theta(7) - theta(2) + u(5) * (-2.654 * cosd(theta(7) - theta

35、(5) + 9.02 * sind(theta(7) - theta(5) + u(8) * (-1.365 * cosd(theta(7) - theta(8) + 11.6 * sind(theta(7) - theta(8) + u(7) * 4.4);dp8 = -0.72 - u(8) * (u(7) * (-1.365 * cosd(theta(8) - theta(7) * sind(theta(8) - theta(7) + u(9) * (-2.256 * cosd(theta(8) - theta(9) + 10.38 * sind(theta(8) - theta(9)

36、+ u(8) * 3.621);dp9 = 0 - u(9) * (17.15 * 1.05 * sind(theta(9) - theta(3) + u(6) * cosd(theta(9) - theta(6) + 5.64 * sind(theta(9) - theta(6) + u(8) * (-2.256 * cosd(theta(9) - theta(8) + 10.38 * sind(theta(9) - theta(8) + u(9) * 3.52);dp_u = dp1/u(1) dp3/u(3) dp4/u(4) dp5/u(5) dp6/u(6) dp7/u(7) dp9

37、/u(9)'udtheta = -B1_ * dp_u;dtheta(1) = udtheta(1) / u(1);dtheta(3) = udtheta(2) / u(3);dtheta(4) = udtheta(3) / u(4);dtheta(5) = udtheta(4) / u(5);dtheta(6) = udtheta(5) / u(6);dtheta(7) = udtheta(6) / u(7);dtheta(8) = udtheta(7) / u(8);dtheta(9) = udtheta(8) / u(9);theta = theta + dtheta;dq4 =

38、 0 - u(4) * (-1.05 * 8.54 * cosd(theta(4) - theta(1) + u(5) * sind(theta(4) - theta(5) - 6.15 * cosd(theta(4) - theta(5) + u(6) * (-0.872 * sind(theta(4) - theta(6) - 5.32 * cosd(theta(4) - theta(6) + 19.945 * u(4);dq5 = -0.1875 - u(5) * (u(4) * (-0.38 * sind(theta(5) - theta(4) * cosd(theta(5) - th

39、eta(4) + u(7) * (-2.654 * sind(theta(5) - theta(7)- 9.02 * cosd(theta(5) - theta(7) + 15.09 * u(5);dq6 = -0.15 - u(6) * (u(4) * (-0.872 * sind(theta(6) - theta(4) * cosd(theta(6) - theta(4) + u(9) * (-0.645 * sind(theta(6) - theta(9)- 5.64 * cosd(theta(6) - theta(9) + 10.92 * u(6);dq7 = 0 - u(7) * (

40、-1.05 * 11.423 * cosd(theta(7) - theta(2) + u(5) * (-2.654 * sind(theta(7) - theta(5) - 9.02 * cosd(theta(7) - theta(5) + u(8) * (-1.365 * sind(theta(7) - theta(8) - 11.6 * cosd(theta(7) - theta(8) + 31.983 * u(7);dq8 = -0.18 - u(8) * (u(7) * (-1.365 * sind(theta(8) - theta(7) * cosd(theta(8) - thet

41、a(7) + u(9) * (-2.256 * sind(theta(8) - theta(9)- 10.38 * cosd(theta(8) - theta(9) + 21.938 * u(8);dq9 = 0 - u(9) * (-1.05 * 17.15 * cosd(theta(9) - theta(3) + u(6) * (-0.645 * sind(theta(9) - theta(6) - 5.64 * cosd(theta(9) - theta(6) + u(8) * (-2.256 * sind(theta(9) - theta(8) - 10.38 * cosd(theta

42、(9) - theta(8) + 33.13 * u(9);dq_u = dq4/u(4) dq5/u(5) dq6/u(6) dq7/u(7) dq8/u(8) dq9/u(9)'dutemp = -B2_ * dq_u;du(4) = dutemp(1); du(5) = dutemp(2);du(6) = dutemp(3); du(7) = dutemp(4); du(8) = dutemp(5); du(9) = dutemp(6);u = u + du;i = i + 1;* (-0.38- 6.15- 5.32- 11.6附录二：改进欧拉法计算功角特性 Matlab程序e

43、nd u theta iclearclcload('Yn.mat' )G11= real(Yn(1,1);G12= real(Yn(1,2);G13= real(Yn(1,3);G21= real(Yn(2,1);G22= real(Yn(2,2);G23= real(Yn(2,3);G31= real(Yn(3,1);G32= real(Yn(3,2);G33= real(Yn(3,3);B12 =imag(Yn(1,2);B13 =imag(Yn(1,3);B21 =imag(Yn(2,1);B23 =imag(Yn(2,3);B31 =imag(Yn(3,1);B32 =

44、imag(Yn(3,2);Tj1 =6.02;Tj2 =47.2;Tj3 = 12.8;E1 = abs(1.082 + 0.04i);E2 = abs(1.069 + 0.105i);E3 = abs(1.098 + 0.298i);PT1 = 0.9;PT2 = 1.57;PT3 = 1.6; delta1 = 5.9; delta2 = 5.61; delta3 = 17.34; omega1 = 1; omega2 = 1; omega3 = 1; t = 0;%1?U ?u ?1|? %2?u ? Cy ?u ?1|?%3?u ? Cy ?Uy ?1|? %1?Uy ? Cy ? U

45、 ? u ? e %2?Uy ? Cy ? U ? u ? e %3?Uy ? Cy ? U ? u ? eh = 0.01;omg1 = omega1; omg2 = omega2; omg3 = omega3; dlt1 = delta1;dlt2 = delta2; dlt3 = delta3;?X a?u ?X a?u ?X a?u ?1|?1|?1|?%1?uy ? C%2?Uy ? C%3?uy ? C%1?uy ? C%2?Uy ? C%3?uy ? Cy ?U 1 e ? X y ?U 1 e ?X y ?u 1 e ?X y ?u 1 e ?X y ?u 1 e ?X y ?

46、u 1 e ?XPE1 = E1 * E1 * G11 + E1 * (E2 * (abs(Yn(1,2) * sind(delta1 - delta2 + atand(G12/B12) + E3 * (abs(Yn(1,3) * sind(delta1 - delta3 + atand(G13/B13);PE2 = E2 * E2 * G22 + E2 * (E1 * (abs(Yn(2,1) * sind(delta2 - delta1 + atand(G21/B21) + E3 * (abs(Yn(2,3) * sind(delta2 - delta3 + atand(G23/B23);

47、PE3 = E3 * E3 * G33 + E3 * (E1 * (abs(Yn(3,1) * sind(delta3 - delta1 + atand(G31/B31) + E2 * (abs(Yn(3,2) * sind(delta3 - delta2 + atand(G32/B32);P1 = PE1;P2 = PE2;P3 = PE3; while t <= 0.99d_delta1_0 = (omega1 - 1) * 360 * 50;d_delta2_0 = (omega2 - 1) * 360 * 50;d_delta3_0 = (omega3 - 1) * 360 *

48、50;d_omega1_0 = (1 / Tj1) * (PT1 - PE1);d_omega2_0 = (1 / Tj2) * (PT2 - PE2);d_omega3_0 = (1 / Tj3) * (PT3 - PE3);delta1_t = delta1 + d_delta1_0 * h;delta2_t = delta2 + d_delta2_0 * h;delta3_t = delta3 + d_delta3_0 * h;omega1_t = omega1 + d_omega1_0 * h;omega2_t = omega2 + d_omega2_0 * h;omega3_t =

49、omega3 + d_omega3_0 * h;* sind(delta1_t* sind(delta2_t* sind(delta3_t- delta2_t- delta1_t- delta1_tPE1_t = E1 * E1 * G11 + E1 * (E2 * (abs(Yn(1,2) + atand(G12/B12) + E3 * (abs(Yn(1,3) * sind(delta1_t - delta3_t + atand(G13/B13);PE2_t = E2 * E2 * G22 + E2 * (E1 * (abs(Yn(2,1) + atand(G21/B21) + E3 * (abs(Yn(2,3) * sind(delta2_t - delta3_t + atand(G23/B23);PE3_t = E3 * E3 * G33 + E3 * (E1 * (abs(Yn(3,1) + atand(