MATLAB实验电力系统暂态稳定分析

《MATLAB实验电力系统暂态稳定分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验三电力系统暂态稳定分析电力系统暂态稳定计算实际上就是求解发电机转子运动方程的初值问题，从而得出S-t和3-1的关系曲线。每台发电机的转子运动方程是两个一阶非线性的常微分方程。因此，首先介绍常微分方程的初值问题的数值解法。一、常微分方程的初值问题（一）问题及求解公式的构造方法我们讨论形如式（3・1）的一阶微分方程的初值问题（3-1）y'M=/（X,y）,a

2、xi（20,l,rn-l）为步长。在等步长的情况下，步长为fb-ah=n用比表示在节点兀处解的准确值）‘，（"）的近似值。设法构造序列｛”•｝所满足的一个方程（称为差分方程）畑=为+力・0（毎，）'昇）（3-2）作为求解公式，这是一个递推公式，从（心，儿）出发，采用步进方式，自左相右逐步算出y（x）在所有节点门上的近似值兀（7=1,2,…,n）o在公式（3-2）>

3、>,为求xw只用到前面一步的值x，这种方法称为单步法。在公式（3-2）中的儿+1由片明显表示出,称为显式公式。而形如（3-3）畑=X+h•屮g,y（,治1Ji）（3-3）的公式

4、称为隐式公式，因为其右端妙中还包括兀+

5、0如果由公式求y田时，不止用到前一个节点的值，则称为多步法。由式（3・1）可得dy=f（x,y）dx（3-4）两边在［石，石+J上积分，得)=yg)+f(x.y(x))dx(3-5)Jxi由此可以看出，如果想构造求解公式，就要对右端的积分项作某种数值处理。这种求解公式的构造方法叫做数值积分法。(二)一般的初值问题的解法1.欧拉法和改进欧拉法对于初值问题(3-1),采用数值积分法，从而得到(3・5)。对于(3-5)右端的积分用矩形公式(取左端点)，则得到fgy(x))d^«h・f(Xi,y(Xi))J

6、勺进而得到(3・1)的求解公式(3-2)X+i+〃•/("，)】)(7二()，1，2,n-1)(3-6)此公式称为欧拉(Euler)格式。如果对式(3-5)右端的积分用梯形公式3+]h/(圮y(x))dx«-(/(xPy(xf))^f(xM,y(xM)))Jxj2则可以得到初值问题(3・1)的梯形求解公式如式(3-7)X+i=儿+£・[/(兀，刀)+/■(可+i，X=i)](20,1,2,n-1)(3-7)式(3-7)是个隐式公式。可以采取先用欧拉格式求一个的初步近似值，记作可+「称之为预报值，然后用预报值歹屮替代式(3-7)右端的)冶，

7、再计算得到)冶，称之为校正值，这样建立起来的预报一校正方法称为改进欧拉格式九=兀+〃•/(»，X)"加=儿+1•[/a,兀)+/(和，九)]<3_8)2.龙格一库塔方法在单步法中，应用最广泛的是龙格一库塔(Runge-kutta)法，简称R—K法。下面直接给出一种四阶的龙格一库塔法的计算公式(3-9)畑+2©+2©+KJoK、=h-f(xi,yi)(3-9)K2=hf(xi+£'%+*KJh1K3=h'f(xi+~O+3K2)K4=〃•/(»•+/?,%+KJ它也称为标准(古典)龙格一库塔法。例3・1研究下列微分方程的初值问题/y=]1

8、+x2y(o)=0解:这是一个特殊的微分方程，其解的解析式可以给出，为x应用龙格一库塔法，取力=0.25,根据式(3・9)编写一段程序，由零开始自左相右逐步算出)

9、题及解答表现在图形上。因此，我们可以不用根据式(3-9)编写较复杂的程序，而只需应用MATLAB提供的常微分方程解题器來解决问题。下面给出用MATLAB编写的解题程序。首先编写描述常微分方程的ODE文件，文件名为myfun',便于解题器调用它。functiondy=myfun(x,y)dy=zeros(1,1);dy=l/(l+x^2)-2*y^2;再编写利用解题器指令求解y的程序。clearx0=0;fori=l:4xm=2*i;yO=O;[x,y]=ode45(fmyfun1z[xOxm],[yO]);formatlongy(leng

10、th(y))endplot(x,y,*-')运行上述程序，在得到几个点的函数值的同时，也得到函数y的曲线，如图3-1所示。图3-1根据运算结果画出y的曲线二、简单电力系统的暂态稳定性(一)物理