MATLAB与工程应用函数2.pptVIP

5.4.1.3 常微分方程初值问题的数值解法 一. 龙格－库塔法简介 一般高阶微分方程总可以写成一阶常微分方程组, 可以应用龙格－库塔法求解. 如: =fj(t,y1,y2 ,…,yn)j=1 , 2 , … , n 四阶龙格－库塔法的计算公式为: 其中K与方程及步长有关. 二. 龙格－库塔法的实现 基于龙格－库塔法，MATLAB提供了求常微分方程数值解的函数，一般调用格式为： [t,y]=ode23(fname,tspan,y0) [t,y]=ode45(fname,tspan,y0) 其中fname是定义f(t,y)的函数文件名，该函数文件必须返回一个列向量。tspan形式为[t0,tf],表示求解区间。y0是初始状态列向量。t和y分别给出时间向量和相应的状态向量。 例5-17 见MATLAB –HELP---ODE45---EXAMPLE1 例5-18 求解著名的Van der Pol方程。 y(1)’’ - 1000*(1 - y(1)^2)*y(1)’ +y(1)=0; 见MATLAB –HELP---ODE45---EXAMPLE2 例5-19 画Van der Pol方程的极限环。 画位移速度曲线，将1000改为0.01，初始条件改为[0，2] 5.4.1.4 函数极值 MATLAB提供了基于单纯形算法