常微分方程数值解与matlab精选.ppt

谢谢同学们！ * * 实验4-常微分方程数值解 1. 求解常微分方程数值方法介绍 （1）一阶微分方程 求方程(1)的数值解, 就是计算(精确)解在一系列离散点 的近似值. 通常取相等的步长h,于是xn=x0+nh(n=1,2,…). (a) 欧拉方法 基本思想是在小区间[xn,xn+1]上用差商 代替方程(1)左端的导数 而方程右端函数f(x,y(x))中的x取[xn,xn+1]上得某一点, 公式为 （2） 实验4-常微分方程数值解 (b) Runge-Kutta方法 基本思想是用小区间[xn,xn+1]上的若干个点的导数的线性组合代替方程(2)右端的 , 一般形式为 (3) 满足 并使(3)的局部截断误差 -------L级p阶Runge-Kutta公式 实验4-常微分方程数值解 (2) 常微分方程组和高阶方程的数值方法 欧拉方法和Runge-Kutta方法可直接推广到求常微分方程组, 如对 欧拉公式为 Runge-Kutta公式有类似的形式. 对高阶方程 (5) 需先降阶化为一阶常微分方程组, 降阶方法不唯一. 简单、常用的方法是令y1=y,将(5)化为 实验4-常微分方程数值解 2. Runge-Kutta方法的MatLab实现 对微分方程(组)的初值问题 Runge-Kutta方法用MatLab命令实现: [t,x]=ode23(@f,ts,x0,options) %用3级2阶Runge-Kutta公式 [t,x]=ode45(@f,ts,x0,options) %用5级4阶Runge-Kutta公式 命令的输入f是待解方程写成的函数M文件: function dx=f(t,x) Dx=[f1;f2;…;fn]; 实验4-常微分方程数值解 2. Runge-Kutta方法的MatLab实现 举例：仿真模拟著名的Lorenz系统混沌图 其中， 先建立一个函数M文件 function xdot=lorenz(t,x) sigma=10;r=28;row=8/3; xdot=[-sigma*x(1)+sigma*x(2); (r-x(3))*x(1)-x(2); x(1)*x(2)-row*x(3)]; 实验4-常微分方程数值解 2. Runge-Kutta方法的MatLab实现 画出Lorenz系统图 clear all;clf; options=odeset(RelTol,1e-5,AbsTol,1e-5); tspan=[0,100];x0=[1,2,3]; [t,x]=ode45(lorenz,tspan,x0,options); l=length(x(:,1)); a=1;b=l; figure(1) plot3(x(a:b,3),x(a:b,1),x(a:b,2),‘b’);grid on;%画出三维相图 xlabel(z);ylabel(x);zlabel(y); figure(2) subplot(311);plot(t,x(a:b,1)) ;%画三分量演化图 subplot(312);plot(t,x(a:b,2)) subplot(313);plot(t,x(a:b,3)) 实验4-常微分方程数值解 2. Runge-Kutta方法的MatLab实现 作业报告：著名的Duffing系统（描述弹簧系统性质） 其中 类似的，分别画出F=1，2，3，4，6等时的相图 翻阅一些参考书，你能得到一些什么结论？ 实验4-常微分方程数值解 3. 实例 问题 缉私艇追击走私船