微分方程数值解修改.ppt

作案时间的确定 一场笔墨官司（放射性废物的处理问题） 结果 ：T(10)=18+42 =25.870， 该物体温度降至300c 需要8.17分钟. T(t)=18+42 , t ≥0. 美国原子能委员会（现为核管理委员会）处理浓缩放射性废物，是将废物放入密封性能很好的圆桶中，然后扔到水深300英尺的海里.他们这种做法安全吗？ 分析：可从各个角度去分析造成危险的因素，这里仅考虑圆桶泄露的可能. 联想：安全 、危险 微分方程模型2 问题的关键 *圆桶至多能承受多大的冲撞速度？(40英尺/秒); *圆桶和海底碰撞时的速度有多大？ 新问题：求这一种桶沉入300英尺的海底时的末速度.（原问题是什么?） 可利用的数据条件： 圆桶的总重量 W=527.327（磅） 圆桶受到的浮力 B =470.327（磅） 圆桶下沉时受到的海水阻力 D=Cv，C=0.08 可利用牛顿第二定律，建立圆桶下沉位移满足的微分方程： 方程的解为 计算碰撞速度，需确定圆桶和海底的碰撞时间t0 ? 分析：考虑圆桶的极限速度 ≈713.86（英尺/秒）＞＞40（英尺/秒） 原问题得到解决了吗? * 目的 内容 2．学会用MATLAB求微分方程的数值解． 1．学会用MATLAB求简单微分方程的解析解． 1．求简单微分方程的解析解． 2．求微分方程的数值解． 3． 数学建模实例． 求微分方程的数值解 （一）常微分方程数值解的定义 （二）建立数值解法的一些途径 （三）用MATLAB软件求常微分方程的数值解 返 回 1．目标跟踪问题一：导弹追踪问题 2．目标跟踪问题二：慢跑者与狗 3．地中海鲨鱼问题 返 回 数学建模实例 微分方程的解析解 求微分方程（组）解析解的命令: dsolve(‘方程1’,‘方程2’,…,‘方程n’,‘初始条件’,‘自变量’) To MATLAB（ff1） 结 果：u = tg(t-c) 解 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x) 结 果 为 : y =3e-2xsin（5x） To MATLAB（ff2） 解 输入命令 ： [x,y,z]=dsolve(Dx=2*x-3*y+3*z, 　　Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t)； x=simple(x) % 将x化简 y=simple(y) z=simple(z) 结 果 为：x = (c1-c2+c3+c2e -3t-c3e-3t)e2t y = -c1e-4t+c2e-4t+c2e-3t-c3e-3t+c1-c2+c3)e2t z = (-c1e-4t+c2e-4t+c1-c2+c3)e2t To MATLAB（ff3） 返 回 微分方程的数值解 （一）常微分方程数值解的定义 在生产和科研中所处理的微分方程往往很复杂,且大多得不出一般解．而实际中的对初值问题，一般是要求得到解在若干个点上满足规定精确度的近似值，或者得到一个满足精确度要求的便于计算的表达式． 因此，研究常微分方程的数值解法是十分必要的． 返 回 为什么要学习微分方程数值解 微分方程是研究函数变化规律的重要工具，有着广泛的应用。如: 物体的运动, 电路的电压, 人口增长的预测 许多微分方程没有解析解，数值解法是求解的重要手段，如 基本内容 3. 实际问题用微分方程建模，并求解 2. 龙格-库塔方法的MATLAB实现 *4. 数值算法的收敛性、稳定性与刚性方程 两个最常用的数值算法: 欧拉（Euler）方法 龙格-库塔（Runge-Kutta）方法 （三）用MATLAB软件求常微分方程的数值解 [t，x]=solver(’f’,ts,x0,options) ode45 ode23 ode113ode15sode23s 由待解方程写成的M文件名 ts=[t0,tf]，t0、tf为自变量的初值和终值 函数的初值 ode23：组合的2/3阶龙格–库塔–费尔贝格算法 ode45：运用组合的4/5阶龙格–库塔–费尔贝格算法 自变量值 函数值 用于设定误差限(缺省时设定相对误差10-3, 绝对误差10-6), 命令为：options=odeset（’reltol’,rt,’abstol’,at）, rt，at：分别为设定的相对误差和绝对误差． 1．在解含