电力系统潮流计算课程设计

1、课 程 设 计 报 告学生姓名:学 号：学 院:电气工程学院班 级:题 目:电力系统潮流计算初 壮指导教师： 职称: 副教授 指导教师： 李翠萍 职称: 副教授 2014年 01月10日摘 要潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压幅值和相角，各元件流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种，牛顿拉夫逊法是数学上解非线性方程

2、组的有效方法，有较好的收敛性。将牛顿拉夫逊法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性，稀疏性及节点编号顺序优划等技巧，使牛顿拉夫逊法在收敛性，占用内存，计算速度等方面的优点都超过了阻抗法。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。首先，画出系统的等效电路图，在计算出各元件参数的基础上，应用牛顿拉夫逊法以及matlab软件进行计算对给定系统图进行了四种不同负荷下的潮流计算，经过调节均得到符合电压限制及功率限制的潮流分布。其次，牛顿拉夫逊法具有较好的收敛性，上述计算过程经过四到五次迭代后均能收敛。根据运算结果，分析各支路损

3、耗和系统总损耗。关键词：牛顿-拉夫逊法 matlab 潮流计算 matpow一、 题目原始资料1系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。变电所1变电所2变电所3变电所435kv母线10kv母线35kv母线10kv母线一次侧电压220kv一次侧电压220kv母线1母线3母线2线路长为90km线路长为80km线路长为90km线路长为100km2*qfq-50-22*qfs-50-2tqn-100-22*tqn-100-2。线路长为80km。电厂一电厂二线路长为70km2、发电厂资料：母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300mw ），母线3为机压母线，机压母线上

4、装机容量为100mw，最大负荷和最小负荷分别为20mw和50mw,发电厂二装机容量为200mw。3、变电所资料：(一) 变电所1、2、3、4低压母线的电压等级分别为：35kv 10kv 35kv 10kv(二) 变电所的负荷分别为：60mw 40mw 70mw 50mw(三) 每个变电所的功率因数均为cos=0.85；(四) 变电所1和变电所3分别配有两台容量为75mva的变压器，短路损耗414kw，短路电压（%）=16.7；变电所2和变电所4分别配有两台容量为63mva的变压器，短路损耗为245kw，短路电压（%）=10.5。4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于

5、图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。二、 课程设计内容及要求课设内容：1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。1) 4个变电所的负荷同时以2%的比例增大；2) 4个变电所的负荷同时以2%的比例下降3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的调整。（变电所

6、低压母线电压10kv要求调整范围在9.5-10.5之间；电压35kv要求调整范围在35-36之间）5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次）6. 最终形成课程设计成品说明书。课设要求：1. 在读懂程序的基础上画出潮流计算基本流程图2. 通过输入数据，进行潮流计算输出结果3. 对不同的负荷变化，分析潮流分布，写出分析说明。4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。5. 打印利用ddrts进行潮流分析绘制的系统图，以及潮流分布图。三、 题目分析3.1牛顿拉夫逊法概要首先对一般的牛顿拉夫逊法作一简单的说明。已知一个变量x

7、函数为：到此方程时，由适当的近似值出发，根据：反复进行计算，当满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿拉夫逊法。这一方法还可以做下面的解释，设第次迭代得到的解语真值之差，即的误差为时，则：把在附近对用泰勒级数展开上式省略去以后部分的误差可以近似由上式计算出来。比较两式，可以看出牛顿拉夫逊法的休整量和的误差的一次项相等。用同样的方法考虑，给出个变量的个方程：对其近似解得修正量可以通过解下边的方程来确定：式中等号右边的矩阵都是对于的值。这一矩阵称为雅可比（jacobi）矩阵。按上述得到的修正向量后，得到如下关系这比更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算

8、满足为预先规定的小正数，是第n次迭代的近似值。3.2牛顿法的框图及求解过程3.2.1用牛顿法计算潮流时，有以下的步骤：a) 给这各节点电压初始值；b) 将以上电压初始值代入公式，求修正方程的常数项向量 ；c) 将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。d) 解修正方程式；e) 修正各节点电压，；f) 将，在带入方程式，求出；g) 检验是否收敛，即h) 如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回（2）进行下次迭代计算，直到收敛为止。3.2.2程序框图如下3.3 matpower的功能与使用方法matpower所用的所有数据文件均

9、为matlab的m文件或者mat文件，用来定义和返回变量、basemva、bus、branch、gen等。 basemva变量是一个标量，用来设置基准容量。对于计算中采用有名值，可以根据需要设置，如ioomva；对于计算中采用标么值，一般设置为1。 bus变量是一个矩阵，用来设置电网中各母线参数，其格式为：bus_i、type、pd、qd、gs、bs、83ea、vm、va、basekv、zone、vmax、vmin。bus-i用来设置母线编号，范围为129997。type用来设置母线类型，1为pq节点母线，2为pv节点母线，3为平衡（参考)节点母线。pd和qd用来设置母线注入负荷的有功和无功功

10、率。gs、bs用来设置与母线并联电导和电纳。basekv用来设置该母线基准电压。vm和va用来设置母线电压的幅值、相位初值。vmax和vmin用来设置工作时母线最高、最低电压幅值。area和zone用来设置电网断面号和省耗分区号，一般都设置为1，前者可设置范围为l100，后者可设置范围为1999。 branch变量是一个矩阵，用来设置电网中各支路参数，其格式为：fbus、tbus、r、x、b、ratea、mteb、mtec、ratio、ajlgle、status。fbus和tbus用来设置该支路由起始节点(母线)编号和终止节点(母线)编号。r、x和b用来设置该支路的电阻、电抗和充电电纳。rat

11、ea、 rateb和rate(；分别用来设置该支路长期、短期和紧急允许功率。ratio用来设置该支路的变比，如果支路元件仅仅是导线为o；如果支路元件为变压器，则该变比为fl，us侧母线的基准电压与thus侧母线的基准电压之比。angle用来设置支路的相位角度，如果支路元件为变压器，就是变压器的转角；如果支路元件不是变压器，相位角度为0度。status用来设置支路工作状态，1表示投入运行，0表示退出运行。 gen变量也是一个矩阵，用来设置接入电网中的发电机(电源)参数，其格式为：bus、魄、qg、qmax、qmin、vg、mbase、status、pmax、pmin。bus用来设置接人发电机(电

12、源)的母线编号。pg和qg用来设置接人发电机(电源)的有功和无功功率。pmax和pmin用来设置接入发电机(电源)的有功功率最大、最小允许值。9max和qmin用来设置接入发电机(电源)的无功功率最大、最小允许值。vg用来设置接人发电机(电源)的工作电压。mbase用来设置接入发电机(电源)的功率基准，如为缺省值，就是baselva变量的值。status用来设置发电机(电源)工作状态，l表示投入运行，0表示退出运行。3.4 问题求解电压是衡量电力系统电能质量的标准之一。电压过高或过低，都将对人身及其用电设备产生重大的影响。保证用户的电压接近额定值是电力系统调度的基本任务之一。当系统的电压偏离允

13、许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值范围内。本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整范围。常规的电力系统潮流计算中一般具有三种类型的节点：pq、pv及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：pq节点（注入有功功率及无功功率），pv节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。3.4.1变量的分类：负荷消耗的有功、无功功率、电源发出的有功、无功功率、母线或节点的电压大小和相位、 、在这十二个变量中，负荷消耗的有功和无功功率无法

14、控制，因它们取决于用户，它们就称为不可控变量或是扰动变量。电源发出的有功无功功率是可以控制的自变量，因此它们就称为控制变量。母线或节点电压的大小和相位角是受控制变量控制的因变量。其中, 、主要受、的控制, 、主要受、的控制。这四个变量就是简单系统的状态变量。为了保证系统的正常运行必须满足以下的约束条件：对控制变量 对没有电源的节点则为对状态变量的约束条件则是 对某些状态变量还有如下的约束条件 3.4.2节点设置及分类 根据系统图可知此系统为两端供电网路，将母线1，2设为节点1，6，将变电所1、2、3、4的高低压侧分别设为节点2、3、4、5、7、8、9、10。并且，将节点1设为平衡节点，将节点6

15、设为pv节点，其余节点设为pq节点。3.4.3参数求取设定基准值，所以根据题目原始资料，计算发电厂、变压器及线路的参数。（1）运用下列公式计算变压器参数：采用变压器参数为折算至高压侧的数值，其变比k1，其中， 75mva的变压器： 两台变压器并联，得阻抗假设我们所采用的变压器有5个抽头，电压调节范围为2*2.5%， 对应的分接头开始时设变压器高压侧主分接头,设四个变电所变压器的非标准变比，降压变压器5个分接头的非标准变比如下，以备调压时选用变压器非标准变比按照上面数据输入后到得到各个接点电压标幺值，具体检验过程如下：可得出变压器二次侧的实际电压：形成节点参数矩阵时，还需要我们计算出各节点所接发

16、电机功率，节点负荷功率（可根据负荷有功功率，发电机有功功率及已知的功率因数计算）变电所1 变电所2 变电所3 变电所4 平衡节点为节点，所设节点电压的初始值为，给定值为1.05pv节点为节点，给定电压值为1.05，。（2）计算线路参数支路阻抗 支路对地电容 支路70km输电线路：支路80km输电线路：支路90km输电线路：支路90km同杆架设双回线：支路80km同杆架设双回线：支路100km同杆架设双回线：3.5 等值电路的计算3.6代入程序并运行从结果中可以观察到各个节点电压标幺值分别为:节点12345678910电压1.0180.9620.9911.0021.0401.0181.0251.

17、0321.0501.050折算到有名值以后，可以观察到节点1、4的电压都与题目给定的范围相比偏高，因此调节变压器分接头和发电厂电压，通过对系统进行多次调整，最终得到合理结果，调整结果如下: 调节方法：调节方法分接头1分接头2分接头3分接头4未调整1111调整后1.0251.01.01.05调节结果： 节点号12345678910未调整1.0180.9620.9911.0021.0401.0181.0251.0321.0501.050调整后0.9920.9620.9910.9511.0401.0181.0251.0321.0501.050 具体的潮流分布见下表：各条支路的首端功率p各条支路的首端

18、功率q各条支路的末端功率p各条支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.03-1.40-0.87s(7,3)1.200.83-1.20-0.74s(8,4)0.800.55-0.80-0.50s(5,9)-0.80-0.530.81-0.10s(6,9)-0.22-1.090.240.87s(6,7)-1.180.061.19-0.19s(7,9)0.62-1.32-0.601.14s(7,8)-3.020.683.08-1.03s(8,10)-3.880.484.00-0.84变压器分接头的调高，使最终调节使得电压在规定范围内时，有功损耗为11.

19、95mw。3.6.1 4个变电所负荷同时以2%的比例增大由于变电所负荷以2%的比例增大，所以各个变电所负荷变为变电所1 =0.816+j0.5057变电所2 =1.428+j0.8850变电所3 =1.224+j0.7585变电所4 =0.816+j0.5057变电所4 =0.816+j0.5057将程序中相应矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：调节方法分接头1分接头2分接头3分接头4调整前1111调整后1.0251.01.01.05 调节结果： 节点电压12345678910调整前1.0170.9600.9891.0001.0391.0171.0

20、241.0311.0501.050调整后0.9910.9600.9890.9501.0391.0171.0241.0311.0501.050 由上表观察到，进行最后一次调节时节点1,4的电压均在题目允许的范围内，变压器分接头的调高，使最终调节使得电压在规定范围内时，有功损耗为12mw。由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下： 各条支路的首端功率p各条支路的首端功率q各条支路的末端功率p各条支路的末端功率qs(5,1)0.820.54-0.82-0.51s(6,2)1.431.06-1.43-0.88s(7,3)1.230.85-1

21、.22-0.76s(8,4)0.820.56-0.82-0.51s(5,9)-0.82-0.540.82-0.08s(6,9)-0.25-1.110.270.89s(6,7)-1.180.051.19-0.18s(7,9)0.58-1.33-0.571.15s(7,8)-3.000.663.06-1.02s(8,10)-3.880.454.00-0.813.7.2 4个变电所负荷同时以2%的比例减小由于变电所负荷以2%的比例减小，所以各个变电所负荷变为： 变电所1 =0.784+j0.4859变电所2 =1.372+j0.8502变电所3 =1.1776+j0.7287变电所4 =0.784+

22、j0.4859将程序中相应矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：调节方法分接头1分接头2分接头3分接头4调整前1111调整后1.0251.01.01.05 调节结果： 节点电压12345678910调整前1.0190.9640.9921.0031.0401.0191.0261.0321.0501.050调整后0.9930.9640.9920.9521.0401.0191.0261.0321.0501.050 由上表观察到，进行最后一次调节时节点1,4的电压均在题目允许的范围内，变压器分接头的调高，使最终调节使得电压在规定范围内时，有功损耗为11.95

23、mw。由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：各条支路的首端功率p各条支路的首端功率q各条支路的末端功率p各条支路的末端功率qs(5,1)0.790.52-0.78-0.49s(6,2)1.381.01-1.37-0.85s(7,3)1.180.81-1.18-0.73s(8,4)0.790.54-0.78-0.49s(5,9)-0.79-0.520.79-0.11s(6,9)-0.19-1.080.200.86s(6,7)-1.190.071.20-0.20s(7,9)0.66-1.30-0.641.13s(7,8)-3.030

24、.703.10-1.05s(8,10)-3.880.514.00-0.873.8 断开环网回线的潮流的分布（1）分析断开6、9支路时的潮流分析当断开6、9支路时，导致支路数变为9程序矩阵改变，其他参数不变。用matlab软件进行潮流计算进行调整。节点号12345678910调整前1.0180.9160.9710.9911.0400.9751.0061.0211.0501.050调整后0.9920.9160.9710.9651.0400.9751.0061.0211.0501.050 由上表观察到，进行最后一次调节后节点1,2,3,4的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后

25、各个支路的有功损耗，记录于下表： 功率损耗pqs(5,1)0.0010.04s(6,2)0.0060.18s(7,3)0.0030.09s(8,4)0.002 0.05s(5,9)0.0050.01s(6,9)0.0000.00s(6,7)0.0250.06s(7,9)0.0380.09s(7,8)0.0590.14总损耗0.2560.93具体的潮流分布如下表：各条支路的首端功率p各条支路的首端功率q各条支路的末端功率p各条支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.05-1.40-0.87s(7,3)1.200.83 -1.20-0.74s(8,4

26、)0.800.55-0.80-0.50s(5,9)-0.80-0.530.81-0.10s(6,9)0.00 0.000.000.00s(6,7)-1.41-1.051.430.96s(7,9)0.39-1.92-0.351.79s(7,8)-3.020.123.08-0.47s(8,10)-3.88-0.084.00-0.27 分析断开6、7支路时的潮流分析各条支路的首端功率p各条支路的首端功率q各条支路的末端功率p各条支路的末端功率qs(5,1) 0.800.53-0.80-0.50s(6,2)1.411.04-1.40-0.87s(7,3)1.200.83-1.20-0.74s(8,4)

27、0.80 0.55-0.80-0.50s(5,9)-0.80-0.530.81 -0.10s(6,9)-1.41-1.041.440.87s(6,7)0.000.00 0.000.0s(7,9)1.81 -1.62-1.751.54 s(7,8)-3.020.793.08-1.14s(8,10)-3.880.604.00-0.95分析断开7、9支路时的潮流分析各条支路的首端功率p各条支路的首端功率q各条支路的末端功率p各条支路的末端功率qs(5,1)0.800.53-0.80-0.50s(6,2)1.411.04 -1.40-0.87s(7,3)1.200.83 -1.20-0.74s(8,4

28、)0.800.55-0.80 -0.50s(5,9)-0.80-0.53 0.81 -0.10s(6,9)0.38-1.81-0.341.65s(6,7)-1.790.771.82-0.85s(7,9)0.000.000.000.00s(7,8)-3.020.013.08-0.36s(8,10)-3.88-0.194.00-0.16 四、总结从本次电力系统潮流计算课程设计中，通过计算在不同的负荷下各支路节点潮流分布的变化，分析负荷变动时，对功率的影响，我们掌握了常用的调整电压的措施。在本次课程设计过程当中，锻炼了自己实际操作分析能力，理论联系实际，对运行中的电力系统，通过潮流计算可以预知各种负

29、荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。实际的情况远比我们计算的情况复杂，这让我深刻了解了潮流计算的重要性。精准的潮流计算不仅可以使电网处于稳定且平衡的状态，这对电网的安全运行起到关键性的作用，还可以实现运行的经济性。所以我认为学好电力系统的关键在于学好潮流计算，电力系统的安全运行不是儿戏，一个很小的失误就可能造成整个电网的崩溃，这给国家和社会带来的危害和损失之大是可想而知的。学好电力系统这门课、学好潮流计算是我们电专业学生的当务之急。没有坚实的基础就不

30、可能建造出高楼大厦，所以，要想在以后的工作岗位上有所成就，就必须打好基础，努力提高自己的专业素养。 参考文献1、电力系统分析课程设计综合实验，祝淑萍 等，中国电力出版社;2、电力系统稳态分析，陈珩，中国电力出版社;3、matlab实用指南，苏金明，王永利，电子工业出版社。4、陈珩.电力系统稳态分析.北京：中国电力出版社，20055、/matpower/附 录function mpc=caseks5% matpower case format:version 2mpc.version=2;% - power flow data - % sys

31、tem mva basempc.basemva=1;1% bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.4 0.8676 0 0 1 1 0 1 1 1.1 0.94 ; 3 1 1.2 0.7436 0 0 1 1 0 1 1 1.1 0.94; 4 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94 ; 5 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 6 1 0

32、0 0 0 1 1 0 1 1 1.1 0.95; 7 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 8 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 9 3 0 0 0 0 1 1 0 1 1 1.1 0.95; 10 2 0 0 0 0 1 1 0 1 1 1.1 0.95; ; % generator data% bus pg qg qmax qmin vg mbase status pmax pminmpc.gen= 9 0 0 10 -10 1.05 1 1 7 0; 10 4 0 10 -10 1.05 1 1 4 0; ; % branch data

33、% fbus tbus r x b ratea rateb ratec ratio angle status angmin angmax mpc.branch= 5 1 1.543*10-3 4.165*10-2 0 0 0 0 1.025 0 1 -360 360; 6 2 1.84*10-3 5.565*10-2 0 0 0 0 1 0 1 -360 360; 7 3 1.543*10-3 4.165*10-2 0 0 0 0 1 0 1 -360 360; 8 4 1.84*10-3 5.565*10-2 0 0 0 0 1.025 0 1 -360 360; 5 9 0.8034*10

34、-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; 6 9 1.2854*10-2 3.0397*10-2 0.2353 0 0 0 0 0 1 -360 360; 6 7 0.8034*10-2 1.8998*10-2 0.1471 0 0 0 0 0 1 -360 360; 7 9 1.1248*10-2 2.6597*10-2 0.2059 0 0 0 0 0 0 -360 360; 7 8 0.6427*10-2 1.5199*10-2 0.4706 0 0 0 0 0 1 -360 360; 8 10 0.8034*10-2 1.8998*10-2

35、 0.5882 0 0 0 0 0 1 -360 360; ; return; function mpc=caseks2% matpower case format:version 2mpc.version=2;% - power flow data - % system mva basempc.basemva=1;1% bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.4 0.8676 0 0 1 1 0

36、 1 1 1.1 0.94 ; 3 1 1.2 0.7436 0 0 1 1 0 1 1 1.1 0.94; 4 1 0.8 0.4958 0 0 1 1 0 1 1 1.1 0.94 ; 5 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 6 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 7 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 8 1 0 0 0 0 1 1 0 1 1 1.1 0.95; 9 3 0 0 0 0 1 1 0 1 1 1.1 0.95; 10 2 0 0 0 0 1 1 0 1 1 1.1 0.95; ; % generato

37、r data% bus pg qg qmax qmin vg mbase status pmax pminmpc.gen= 9 0 0 10 -10 1.05 1 1 7 0; 10 4 0 10 -10 1.05 1 1 4 0; ; % branch data % fbus tbus r x b ratea rateb ratec ratio angle status angmin angmax mpc.branch= 5 1 1.543*10-3 4.165*10-2 0 0 0 0 1.025 0 1 -360 360; 6 2 1.84*10-3 5.565*10-2 0 0 0 0

38、 1 0 1 -360 360; 7 3 1.543*10-3 4.165*10-2 0 0 0 0 1 0 1 -360 360; 8 4 1.84*10-3 5.565*10-2 0 0 0 0 - 0 1 -360 360; 5 9 0.8034*10-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; 6 9 1.2854*10-2 3.0397*10-2 0.2353 0 0 0 0 0 1 -360 360; 6 7 0.8034*10-2 1.8998*10-2 0.1471 0 0 0 0 0 1 -360 360; 7 9 1.1248*10

39、-2 2.6597*10-2 0.2059 0 0 0 0 0 1 -360 360; 7 8 0.6427*10-2 1.5199*10-2 0.4706 0 0 0 0 0 1 -360 360; 8 10 0.8034*10-2 1.8998*10-2 0.5882 0 0 0 0 0 1 -360 360; ; return; function mpc=caseks31% matpower case format:version 2mpc.version=2;% - power flow data - % system mva basempc.basemva=1; % bus data % bus_i type pd qd gs bs area vm va basekv zone vmax vminmpc.bus= 1 1 0.816 0.5057 0 0 1 1 0 1 1 1.1 0.94; 2 1 1.428 0.8850 0 0 1 1 0 1 1 1.1