电气工程自动化毕业论文.doc

《电气工程自动化毕业论文.doc》由会员分享，可在线阅读，更多相关《电气工程自动化毕业论文.doc（36页珍藏版）》请在三一办公上搜索。

1、湖南科技大学毕业论文基于Matlab的电力系统潮流仿真计算作者：系别：成教专业：电气工程自动化学号： 指导老师 ： 时间：2011年2月13日目 录摘要1第一章 电力系统潮流计算简介 2第二章 电力网络的数学模型-节点电压方程4 2.1 电力网络的节点电压方程2 2.2节点导纳矩阵2 2.2.1 自导纳和互导纳的确定方法2 2.2.2 节点导纳矩阵的性质和意义4 2.2.3 非标准变比变压器等值电路5 第三章 潮流计算的数学模型-潮流方程7 3.1 潮流计算的基本方程7 3.2 电力系统节点分类10 3.3 潮流计算的约束条件11第四章 牛顿-拉夫逊法潮流计算12 4.1 牛顿-拉夫逊法基本原

2、理12 4.2 牛顿-拉夫逊法潮流计算方法14 4.2.1 以直角坐标形式表示14 4.2.2 牛顿-拉夫逊法程序框图17 4.2.3 牛顿-拉夫逊法求解过程18 4.2.4 程序需要输入数据19总结20参考文献21附录 牛顿拉夫逊法潮流计算程序及其结果分析22基于Matlab的电力系统潮流仿真计算学员姓名：单位：摘要：Matlab是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序的难度并节省了时间，从而可把主要的精力集

3、中在算法的构思而不是编程上。传统的潮流计算程序缺乏图形用户界面，结果显示不直观，难于与其他分析功能集成，网络原始数据输入工作量大且易于出错。Matlab潮流计算研究近年来得到了长足的发展。针对这一现状结合电力系统的基本特点，以牛顿拉夫逊法潮流计算方法为例，对IEEE-6BUS标准试验系统的潮流计算进行仿真,提出了基于Matlab的电力系统潮流仿真计算。它大大提高了计算速度，占用内存少；计算结果有良好的可靠性和可信性；适应性好，即能处理变压器变比调整，系统元件的不同描述和与其它程序配合的能力强。关键词：电力系统潮流计算;MATLAB仿真;牛顿拉夫逊法潮流计算;面向对象AbstractThe MA

4、TLAB is a kind of handing over with each other type, face to the procedure design language of the object, extensive apply in the industry field and academic circles, used for primarily the matrix carries to calculate, also having the mighty function in the aspects of number analyzing, automaticly co

5、ntrol imitating, arithmetic figure signal handling, development analyzing, painting waiting at the same time.Pass the language of M, can weave to write the calculate way with the way of the similar mathematical formula, lowered consumedly the procedure a difficulty for needing combined to save time,

6、 from but can concentrate the main energy conceive outline in the calculate way but is not to weave the distance top.Traditional current calculation the procedure lacks the sketch customer interface, result manifestation do not keep the view, difficult in gather with other analysis function, network

7、 primitive data input the work has great capacity and apt to come amiss.The MATLAB current calculation research got the substantial development in recent years.This present condition that aim at the basic characteristics joins together the electric power system, with Newton-Raphson method current ca

8、lculation method for a current for, to the standard of BUS IEEE-6 experimenting system calculation proceeds to imitate true, put forward according to the MATLAB electric power system current imitates true calculation.It increased consumedly calculation speed, take up the inside saves little;Calculat

9、ion result contain good and dependable with can letter;The adaptability is good, can handle namely the transformer changes the ratio adjustment, the system an ability for of differently describing with with other procedure matching with is strong.Key Words:Electric power system current calculation;

10、the MATLAB imitates true; Newton-Raphson;Face to the object第一章 电力系统潮流计算的简介所谓电力系统的潮流计算就是采用一定的方法确定系统中各处的电压和功率分布实为功率流（power flow），电力界惯称潮流。电压（包括幅值U和相位）和功率（包括有功功率P和无功功率Q ）是表征电力系统稳定运行的主要物理量。电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需

11、要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行电力系统潮流计算从50年代中期就已经开始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几点：（1）计算方法的可靠性或收敛性；（2）对计算机内存量的要求；（3）计算速度

12、；（4）计算的方便性和灵活性。电力系统的潮流计算和一般交流电路计算的根本差别是：后者已知和待求的是电压和电流，而前者是电压和功率。正是这一差别决定了二者本质上的不同：描述交流电路的方程，如节点电压方程，回路电流方程是线性方程，而描述电力系统稳定运行特性的潮流方程是非线性方程。由此似使得求解方法有了根本不同：线性方程可直接采用消去法求解，而非线性方程只能采用迭代法求解。这就是电力系统潮流问题的特点：已知和待求的是电压和功率，为非线性关系，需迭代求解。牛顿-拉夫逊法是求解非线性方程有效的迭代计算。它能可靠的收敛，并给出正确答案。电力系统中进行电力系统潮流计算的目的是：（1）确定电力系统的运行方式。

13、（2）检查各系统中的各元件是否过压，过载。（3）为电力系统继电保护的整定提供依据。（4）为电力系统的稳定计算提供初值。（5）为电力系统规划和经济运行提供分析基础。利用计算机进行电力系统潮流计算分为以下几步：建立数学模型；确定计算方法；制定计算流程图并编制程序；上机调试并对结果进行分析。第二章 电力网络的数学模型2.1 电力网络的节点电压方程 (2-1)为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点为负。既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者代数之和。为节点电压列向量，由于节点电压是对称于参考节

14、点而言的，因而需先选定参考节点。在电力系统中一般以地为参考节点。如整个网络无接地支路，则需要选定某一节点为参考。设网络中节点数为（不含参考节点），则,均为n*n列向量。为n*n阶节点导纳矩阵。2.2 节点导纳矩阵 2.2.1自导纳和互导纳的确定算法。节电导纳矩阵的节点电压方程：展开为： : (2-2)是一个n*n阶节点导纳矩阵，其阶数就等于网络中除参考节点外的节点数。例如：课程设计所给图中，n=6。节点导纳矩阵的对角元素 (i=1,2,n)成为自导纳。自导纳数值上就等于在i节点施加单位电压，其他节点全部接地时,经节点i注入网络的电流，因此，它可以定义为： (2-3) 图1 图2以图二所示网络为

15、例,取i=1,在节点1接电压源,节点2,3的电压源短接,按定义,可得: ()从而 .由此可得,节点i的自导纳数值上就等于与节点直接连接的所有支路导纳的总和。节点导纳矩阵的非对角元素 (j=1,2,n;i=1,2,.,n;j=i)称互导纳,由此可得互导纳数值上就等于在节点i施加单位电压,其他节点全部接地时,经节点j注入网络的电流,因此可定义为: (2-4)以图2为例,取=2,在节点2接电压源,节点1,3的电压源短接,按定义可得: 从而,。 由此可见,节点j,i之间的互导纳数值上就等于连接节点j,i支路到导纳的负值。显然,恒等于.互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且,由于每个

16、节点所连接的支路数总有一个限度，随着网络中节点数的增加非零元素相对愈来愈少,节点导纳矩阵的稀疏度,即零元素数与总元素的比值就愈来愈高。根据定义直接求取节点导纳矩阵时，注意以下几点:(1)节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。(2)节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。(3)节点导纳矩阵的对角元素就等于各该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。(4)节点导纳矩阵的非对角元素等于连接节点i,j支路导纳的负值。因此，一般情况下，节点导纳矩阵

17、的对角元素往往大于非对角元素的负值。(5)节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求求取这个矩阵的上三角或下三角部分。2.2.2 节点导纳矩阵的性质和意义是n*n阶方阵，其对角元素 (i=1,2,-n)称为自导纳，非对角元素 (i,j=1,2,n, )称为互导纳。将节点电压方程展开为可见, (2-5)表明，自导纳在数值上等于仅在节点i施加单位电压而其余节点电压均为零（即其余节点全部接地）时，经节点i注入网络的电流。其显然等于与节点i直接相连的所有支路的导纳之和。同时可见.表明，互导纳在数值上等于仅在节点j施加单位电压而其余节点电压均为零时，经节点i注入网络的电流，

18、其显然等于()即=。为支路的导纳，负号表示该电流流出网络。如节点ij之间无支路直接相连，则该电流为0，从而=0。注意字母几种不写法的不同意义：粗体黑字表示导纳矩阵，大写字母代表矩阵中的第i行第j列元素，即节点i和节点j之间的互导纳。小写字母i,j支路的导纳等于支路阻抗的倒数数，。节点导纳矩阵的性质：(1)为对称矩阵，=。如网络中含有源元件，如移相变压器，则对称性不再成立。(2)对无接地支路的节点，其所在行列的元素之和均为零，即 。对于有接地支路的节点，其所在行列的元素之和等于该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。(3)具有强对角性：对角元素的值不小于同一行或同一

19、列中任一元素(4)为稀疏矩阵，因节点i ,j 之间无支路直接相连时=0，这种情况在实际电力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元素与全部元素之比，即 , 式中Z 为中的零元素。S 随节点数n 的增加而增加：n=50，S可达92%；n=100，S 可达90%；n=500，S可达99%，充分利用节点导纳矩阵的稀疏性可节省计算机内存，加快计算速度，这种技巧称为稀疏技术。2.2.3 非标准变比变压器等值电路.变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算.双绕组变压器可用阻抗与一个理想变压器串联的电路表示.理想变压器只是一个参数,那就是变比。现在变压器阻抗按实

20、际变比归算到低压侧为例,推导出变压器型等值电路.图3 双绕组变压器原理图图4 变压器阻抗归算到低压侧等值模型流入和流出理想变压器的功率相等 (2-6)式中, 是理想变压器的变比,和 分别为变压器高,低绕组的实际电压.从图四直接可得： （2-7） 从而可得: （2-8）式中,又因节点电流方程应具有如下形式: （2-9）将式（1-8）与（1-9）比较，得： 因此可得各支路导纳为: （2-10） 由此可得用导纳表示的变压器型等值电路:图5 变压器等值电路第三章 潮流计算的数学模型-潮流方程本章从节点电压方程入手导出潮流计算的数学模型-潮流方程。3.1 潮流计算的基本方程在潮流问题中，任何复杂的电力系

21、统都可以归纳为以下元件（参数）组成。（1）发电机（注入电流或功率）（2）负荷（注入负的电流或功率）（3）输电线支路（电阻，电抗）（4）变压器支路（电阻，电抗，变比）（5）母线上的对地支路（阻抗和导纳）（6）线路上的对地支路（一般为线路充电点容导纳）集中了以上各类型的元件的简单网络如图：图6 潮流计算用的电网结构图图7 潮流计算等值网络采用导纳矩阵时，节点注入电流和节点电压构成以下线性方程组 (3-1)其中 可展开如下形式 (3-2)由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中的注入电流用节点注入功率来表示。节点功率与节点电流之间的关系为 (3-3)式中，因此用导纳矩阵时，

22、PQ节点可以表示为把这个关系代入式中 ，得（3-4）式（3-4 ）就是电力系统潮流计算的数学模型-潮流方程。它具有如下特点：（1）它是一组代数方程，因而表征的是电力系统的稳定运行特性。（2）它是一组非线性方程，因而只能用迭代方法求其数值解。（3）由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式-极坐标形式，直角坐标形式和混合坐标形式。a.取 ,得到潮流方程的极坐标形式: (3-5)b. 取 , ,得到潮流方程的直角坐标形式: (3-6)c.取, ,得到潮流方程的混合坐标形式: (3-7)不同坐标形式的潮流方程适用于不同的迭代解法。例如：利用牛顿-拉夫逊迭代法

23、求解，以直角坐标和混合坐标形式的潮流方程为方便；而P-Q解耦法是在混合坐标形式的基础上发展而成，故当然采用混合坐标形式。（4）它是一组n个复数方程，因而实数方程数为2n个但方程中共含4n个变量：P,Q,U和，i=1,2,n，故必须先指定2n个变量才能求解。3.2 电力系统节点分类用一般的电路理论求解网络方程，目的是给出电压源（或电流源）研究网络内的电流（或电压）分布，作为基础的方程式，一般用线性代数方程式表示。 然而在电力系统中，给出发电机或负荷连接母线上电压或电流（都是向量）的情况是很少的，一般是给出发电机母线上发电机的有功功率（P）和母线电压的幅值（U），给出负荷母线上负荷消耗的有功功率（

24、P）和无功功率（Q）。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然的把节点分成三种类型。 1、PQ节点：对这类节点指定P和Q，U和待求。电力系统绝大多数节点均属此类。如变电所母线节点，其无电源功率，负荷功率又已知，故该节点的注入功率已知；又如一些按指定有功和无功功率发电的电厂，其,指定，所带机端负荷也已知，从而节点的注入功率已知。 2、PU节点：对这类节点指定P和U，Q和待求。电力系统中此类节点少数，个别小系统基本没有。设置节点是为了控制该点的电压为一定值，从而保证电力系统的质量，为控制电压必须要有一定的无功功率可供调节。故这类节点是有一定无功

25、储备的发电厂和装有无功电源（电容器，调相机或静止无功补偿器）的变电所，这类节点也称为电压控制节点。3、平衡节点：对这类节点指定U和，其有功功率P和无功功率Q由保证全系统功率平衡的条件决定，因而又称为平衡节点。一般取其=0，电力系统潮流计算中必须有而且只有一U个节点，负责系统频率调整的主调频厂基本上起着平衡节点的作用。3.3 潮流计算的约束条件 电力系统运行必须满足一定技术和经济上的要求.这些要求够成了潮流问题中某些变量的约束条件,常用的约束条件如下:1、 节点电压应满足 (3-8)从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU节点电压幅值必须按上述条件

26、给定。因此，这一约束条件对PQ节点而言。2、 节点的有功功率和无功功率应满足 (3-9) PQ节点的有功功率和无功功率,以及PU节点的有功功率,在给定是就必须满足上述条件,因此,对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验.3、 节点之间电压的相位差应满足 (3-10) 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行

27、方式，重新进行计算。第四章 牛顿-拉夫逊法潮流计算4.1 牛顿拉夫逊法的基本原理设有单变量非线性方程 （4-1）求解此方程式时，先给出解的近似值，它与真解的误差为，则将满足方程（4-1），即 （4-2）将（4-2）式左边的函数在附近展开泰勒级数得：（4-3）式中，分别为函数在处的一阶导数，n 阶导数。如果差值很小，（4-3）式右端的二次及以上阶次的各项均可略去。因此（4-3）可以简化为： （4-4）这是对于变量的修正量的线性方程式，也称修正方程式。由这个方程式可以解出修正量，即 （4-5）用所求得的去修正近似解，便得由于式（3-4）是略去了高次相的简化式，因此所解出的修正量也只是近似值。修正后

28、的近似解同真解仍然有误差。但是，这样的迭代计算可以反复进行下去，迭代计算的通式是 （4-6） 迭代过程的收敛判据为： （4-7）或 （4-8）式中和为预先给定的小正数。牛顿法也可以应用于求解任意多个变量的非线性方程组。假定已给出各变量的初值，令，分别为各变量的的修正量，使其满足方程，即 （4-9）将上式中的n个多元函数在初值附近分别展开成泰勒级数，并略去含有，的二次及以上的阶次的各项，于是便得： （4-10） 方程式（4-10）也可以写成矩阵形式： （4-11）方程式（4-11）是对于修正量，的线性方程组，称为牛顿法的修正方程式。利用高斯消去法或三角分解法 可以解出修正量，。然后就可以对初始近

29、似解进行修正：（4-12）如此反复进行迭代，在进行k+1次迭代时，从求解修正方程式 （4-13）可以得到修正量，。然后按下式 （4-14）对各变量进行修正。式（4-13）和4-14）也可以缩写为： （4-15）和 （4-16）式中，X和分别是由n个变量和修正量组成的n维列向量；是由n个多元函数组成的维列向量；是n*n阶方阵，称为雅可比矩阵，它的第i,j个元素是第i个函数对j个变量的的饿偏导数；上角标（）表示阵的每一个元素都在点处取值。迭代过程一直进行到满足收敛判据 （4-17）或 （4-18）为止。和为预先给定的小正数。4.2牛顿-拉夫逊法潮流计算方法4.2.1 以直角坐标形式表示把牛顿法用于

30、潮流计算，首先应将潮流方程的右端展开，并且分开实部和虚部。采用直角坐标时，节电电压表示为： 导纳矩阵元素则表示为：将上述表示式代入潮流方程式的右端，展开并分出实部和虚部，便得： （4-19）按照节点的分类：（1）PQ节点的有功功率和无功功率是给定的，第个节点的给定功率设为和。假定系统中的第1，2， m 号节点为PQ 为节点，对其中每一个节点可列方程： （4-20）（2）PU节点的有功功率和节点电压幅值是给定的。假定系统中的第m+1,m+2,mn-1为PU为节点，则对其中每一个节点可列方程： （4-21）（3）第n号节点为平衡节点，其电压是给定的，故不需要参加迭代。式（4-20）和（4-21）总

31、共包含了2(n-1)个方程，待求的变量也是2(n-1)个。我们还可以看到，方程（4-20）和（4-21）已经具备了方程组的形式。因此，不难写出如下的修正方程式： (4-22)式中， ; 上述方程中雅可比矩阵的各元素，可以对（4-20）和（4-21）式求偏导数获得。当时， （4-23）当时， (4-24)从表达式（4-23），（4-24）可以看出雅可比矩阵具有以下特点：（1）矩阵中各元素是节点电压的函数，在迭代过程中，这些元素随节点电压的变化而变化。（2）导纳矩阵中的某些非对角元素为零时，雅可比矩阵中对应的元素也为零，若，则必有。（3）雅可比矩阵不是对称矩阵。4.2.2 牛顿-拉夫逊法程序框图

32、程序框图如下图所示： 图8 牛顿-拉夫逊法程序框图4.2.3牛顿-拉夫逊法潮流计算求解过程(1)给出各节点的电压初始值。(2) 将以上的电压初始值代入式（4-20）和（4-21）计算各类节点的不平衡量和。(3) 将电压初始值再代入式 （4-23）和（4-24）,求出修正方程式中雅可比矩阵的各元素。(4) 解修正方程式(4-22),求出修正量 。(5)修正各节点电压:, 。(6)将,再代入式(4-20),(4-21),求。(7)按条件（4-17）校验是否收敛，即 。(8)如果收敛，迭代到此结束，转入计算各线路潮流和平衡节点的功率，并打印输出计算结果,如果不收敛,则转回第(2)进行下一次迭代计算,

33、直到收敛为止。4.2.4程序需要输入的数据（1）n为节点数,nl为支路数，isb为平衡母线节点号（固定为1），pr为误差精度。（2）输入由支路参数形成的矩阵B1矩阵B1的每行是由以下参数构成的：某支路的首端号P；末端号Q，且 PQ; 支路的阻抗（R+jX）;支路的对地容抗；支路的变比K；折算到哪一侧的标志（如果支路的首端P处于高压侧则请输入“1”，否则请输入“0”）。（3）输入各节点参数形成的矩阵B2矩阵B2的每行是由下列参数构成的：节点所接发电机的功率 ；节点负荷的功率；节点电压的初始值；节点电压的给定值；节点所接的无功补偿设备的容量；节点分类标号igl a.PQ节点，b:PV节点,c:平衡

34、节点。 (4)输入由节点号及其对地阻抗形成的矩阵X。总 结电力系统潮流计算是电力系统分析的基本计算，其目标是求取电力系统在给定运行状态的相关参数，即节点电压和功率分布，用以检查系统各元件是否过负荷、各点电压是否满足要求、功率的分布和分配是否合理以及功率损耗等。其计算结果是对复杂电力系统正常和故障条件下稳态运行状态分析的基础，在规划设计中，用于接线方式、电气设备以及导线截面选择；在电网运行中，用于确定运行方式、制订检修计划和确定调整电压的措施；为继电保护、电力系统自动化操作提供设计与整定数据。在电力系统的多个领域都要用到电网潮流计算。传统的手工计算工作量大，易出错，潮流计算费时费力。目前，计算机

35、已成为潮流计算的重要工具。利用牛顿-拉夫逊法、采用Matlab语言进行潮流计算，提高了准确性，使用方便，已得到成功应用。参考文献1.电力系统分析 何仰赞、温增银华中科技大学出版社2.电力系统分析 孟祥萍、高燕高等教育出版社3.电力系统分析 陈怡、蒋平、万秋兰等中国电力出版社4.电力系统 李海燕中国电力出版社5.数值分析 谢伟松天津大学出版社6.Matlab电力系统设计与分析吴天明、彭彬国防工业出版社7.Matlab程序设计及应用 张铮中国铁道出版社8Matlab程序设计与应用 刘卫国、张颖高等教育出版社附录牛顿-拉夫逊法潮流计算程序及其结果分析（六节点）f3=fopen(NodeNum2.tx

36、t);node=fscanf(f3,%d,1 inf);fclose(f3);NodeMax=node(6);NodePVNum=node(1);NodeSlack=node(6);PVNo=node(5);BranchNum=7;%输入线路参数结束f4=fopen(BranData2.txt);Bran=fscanf(f4,%d %d %f %f %f,5 inf);Bran=Bran;v1234=zeros(21,2*(PVNo-1);%tap=0;y=zeros(NodeMax,2*NodeMax);for t=1:size(Bran,1), p=Bran(t,1); q=Bran(t,

37、2); r=Bran(t,3); x=Bran(t,4); k=Bran(t,5); model=x*x+r*r; if k0, y(p,2*p)= y(p,2*p)-k;y(q,2*q)=y(q,2*q)-k;end if k=0, k=1; end y(p,2*p-1)=y(p,2*p-1)+r/(model*k*k); y(p,2*p)=y(p,2*p)+(-1)*x/(model*k*k); y(q,q*2-1)=y(q,q*2-1)+r/(model); y(q,q*2)=y(q,q*2)+(-1)*x/(model); y(p,q*2-1)=(-1)*r/(model*k); y(

38、p,q*2)=x/(model*k); y(q,p*2-1)=y(p,q*2-1); y(q,p*2)=y(p,q*2); end y%for t=2:2:2*NodeMax, % y(:,t)=j*y(:,t); %end%输入PQV参数f5=fopen(PQVData2.txt);pqv=fscanf(f5,%d %f %f %f %f,5 inf);pqv=pqv;fclose(f5);s0=zeros(NodeMax,2);v=zeros(NodeMax,2);v20=zeros(2*(NodeMax-1),1);s0(:,1)=pqv(:,2);s0(:,2)=pqv(:,3);v(:,1)=pqv(:,4);v(:,2)=pqv(:,5);for m=PVNo:(NodeMax-1),