第7章 微分方程实验与计算机模拟.pptVIP

*/16 * */16 * */16 微分方程与计算机模拟 常微分方程数值解方法 捕食者与被捕食者问题 有阻力抛射曲线问题 卫星轨道模拟问题 ? ? ? ? 数值方法求常微分方程初值问题 求解步骤: (1)用函数文件定义一阶微分方程(或方程组)右端函数； (2)用MATLAB命令ode23()求数值解或绘积分曲线。 使用格式：[T,Y] = ode23(F,Tspan,y0) 其中,Tspan = [t0，tN]是常微分方程求解区域，y0是初始值,‘F’ 是包括函数文件名字的符串。 返回值(T,Y) 是求解区域内离散数据及对应数值解。 例7.1 马尔萨斯模型,以1994 年我国人口为12亿为初值，求解常微分方程 N(t)表示人口数量,取人口变化率r =0.015,微分方程 function z=fun1(t,N) z=0.015*N; ode23(fun1,[1994,2020],12) [T,N]=ode23(fun1,[1994,2020],12) 命令窗口 ? 编辑窗口 ? 常微分方程组初值问题 一阶常微分方程组初值问题数值求解方法 [T,y] = ode23( F ,Tspan,y0) 其中, F是函数文件, 表示 微分方程右端函数 Tspan = [t0 Tfinal] —— 求解区域; y0 ——