第3讲：求导积分与微分方程数值解(第2次课).pptVIP

第三讲 微分方程数值解 内容：本讲首先以单摆微分方程求解过程引入， 介绍基于MATLAB的微分方程求解函数， 然后重点讲解微分方程(组) 图形图像解法， 最后简介 Differential Equation Editor 目的：掌握微分方程数值解的一般思路和方法 要求：能够处理应用类型微分方程数值解问题 掌握Malthus人口模型参数计算 (本讲实验题目) 掌握常用求解函数 dsolve ode23 ode45 ode15s 掌握图形图像求解方法 斜率场 / 相平面 / 等值线 了解基于Simulink的微分方程数值仿真工具 DEE 大多数微分方程无法求解析解? 微分方程是研究函数变化规律的有力工具，在科技、工程、经济管理、及生态、环境、人口、交通等各个领域有着广泛的应用。 建立微分方程可以依据物理的、或其他原理和规律建立的平衡关系，但是！更重要的问题是如何求解这些微分方程(组) 部分微分方程可以求得解析解，但是绝大多数的非线性、变系数微分方程或“难以求解”或“求不出解”，所以对于实际问题，研究微分方程的数值解具有重要意义！ 我们所熟悉的微分方程？ 微分方程初值问题的最简单形式： 单摆微分方程求解：建立方程 引例：单摆微分方程求解，由已知运动规律建立微分方程，在求解过程中采用了两种方法，方法一近似简化，方法二求数值解 (实验室练习) 单摆微分方程求解：求近似解 单摆微分方程求解：求数值解 由简化方程求得近似解为： 单摆微分方程求解：求数值解 然后是主调指令，也可写成主调文件loaddanbai.m warning off ts=0; tf=10; a0=0.1745; cond0=[a0,0]; %初始化变量 [t,x]=ode23(danbai,ts,tf,cond0); %调用ode23函数求解 g=9.8; l=25; w=sqrt(g/l); y=a0*cos(w*t); %近似解 [t,x(:,1),y] %输出t对应的数值解和近似解 subplot(1,2,2); stem(t,x(:,1),ro); title(数值解) subplot(1,2,1); hold on; stem(t,y, bp); plot(t,y, b-); title(近似解) 用dsolve函数求解微分方程 MATLAB求解微分方程解析解的函数dsolve Symbolic solution of ordinary differential equations . Syntax~ r = dsolve(eq1,eq2,...,cond1,cond2,...,v) 题例1：p49-4.4.1/ex1，ex2 dsolve(Dy=1+y^2) dsolve(Dtheta=1+theta^2,theta(0)=1,xi) dsolve(x^2*D2y+x*Dy+(x^2-(1/2)^2)*y=0,y(pi/2)=2,Dy(pi/2)=-2/pi,x) pretty(ans) 提示：一些需要注意的细节… 用dsolve函数求解微分方程 MATLAB求解微分方程组解析解的函数dsolve 题例2：求解 dsolve(Dy+2*x*y=x*exp(-x^2),x) 题例3：求解 dsolve((x^2-1)*Dy+2*x*y-cos(x)=0,y(0)=1,x) 题例4：求解 dsolve(D2y+3*Dy+exp(x)=0,x) 用dsolve函数求解微分方程组 MATLAB求解微分方程组解析解的函数dsolve 题例5： p50-4.4.1/ex3，ex4 [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g) [f,g]=dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1,x) 下面的指令有否区别？ dsolve(Dy=x*sin(x)/cos(y)) dsolve(Dy=x*sin(x)/cos(y),x) 提示： 用dsolve求解存在解析解的微分方程相当方便，在“只要结果，不求过程”的场合，节约了大量时间。 课本引例: Malthus人口模型计算 Malthus认为单位时间内人口净增长率为常数： 微分方程数值的相关方法 求微分方程数值解的方法很多，比如：欧拉法，龙格-库塔法等。Runge-Kutta法本质上是间接使用Taylor级数的一种技术，有兴趣的同学可以参考微分方程数值解相关专著。 首选ode45 ，次选ode15s，低精度可选ode23 基于R