微分方程-精选（公开课件）.pptVIP

课后作业 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 疾病传播的微分方程模型 目前甲型H1N1流感肆虐，试建立传染病流行的微分方程数学模型，对疾病流行情况进行模拟和分析 我们采用SIR模型，即将人群分成三种人：S健康者，I感染者，R治愈者。 假设健康者与感染者接触就会感染流感； 感染者经过一段时间治疗后就变成治愈者 治愈者不再会被感染 我们设健康者、感染者和治愈者在第t天的比例分别是s(t),i(t)和r(t)。 试写出描述s(t),i(t)和r(t)三者变化的微分方程组，并用Matlab求解。 讨论不同参数对方程解的影响 * * * Matlab基础及应用——微分方程 张江 北京师范大学管理学院 系统科学系 今日内容 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 微分方程 降落伞问题 Dsolve进行数学求解 微分方程组 lotka-volterra方程组 利用ode23命令进行数值求解 将高阶微分方程化为微分方程组 降落伞模型 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 设人+降落伞的总重量为m=200kg， 下落过程中受到的阻力f=bv，其中v为降落伞的速度，b为一个常数b=196N.s/m 设降落伞从1000m的高空落下 求下落位置和时间t的关系，经过多少时间到达地面？下落过程中的最大速度多少？ 降落伞模型 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 mg f=bv 完整的数学模型 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 mg f=bv 微分方程 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 含有未知函数，及其未知函数导数的方程称为微分方程 根据未知函数中所含变元的个数可以把微分方程分为： 常微分方程 偏微分方程 微分方程与微分方程组 含有多个未知函数的方程组为微分方程组 Matlab中求解微分方程的命令 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 [y1,y2,…,y12]=dsolve(eq1,eq2,...,cond1,cond2,...,v) 输入有3个部分： 符号化的微分方程，符号化的初始条件，要求解的自变量（默认为t） 如果没有初始条件，则返回的函数中将包含常数C1,C2,… Matlab求解 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 v=dsolve(200*Dv+196*v-200*9.8=0,v(0)=0,t) ezplot(v,[0,100]); 求解最大速度 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 syms t; vmax=limit(v,t,inf) 观察v(t)的图形发现，v存在一个渐进的最大值，利用极限的方法来求该最大值 求解距离函数x(t) 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 x=int(v,0,t) figure; ezplot(x,[0,100]); 对v(t)的积分，得到距离函数x(t) 求解距离函数x(t) 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 eqn=strcat(char(x), - 1000); fzero(inline(eqn),100) 通过求解方程x(t)=1000，得到下落的时间t* 求解下列方程 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 y=dsolve(‘D2y+y=1-t^2’) 用dsolve求解微分方程组 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 求解下列方程组 [u v]=dsolve(Du=3*u-2*v,Dv=2*u-v) 练习1 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 弹簧振子 m 物块质量m=1,弹簧的倔强系数k=1,在光滑的水平面上，初始位移为1cm， 初速度为0，求物块的运动方程 x 0 捕食与被捕食者模型 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 设狼种群在任意时刻的个数为Y(t)，兔子种群在任意时刻的个数为X(t) 关于两个未知函数X(t)和Y(t)的方程组 X Y aX kXY cY pXY 捕食与被捕食者模型 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 X Y aX kXY cY pXY 尝试用dsolve求解 北京师范大学管理学院 系统科学系 张江 MATLAB 基础及应用 X Y aX kXY cY pXY [x y]=dsolve(Dx=x-0.01*x*y,Dy=0.02*x*y-y,x(0)=20,y(0)=20,t) 采用数