7-第七章-常微分方程1.pdfVIP

第七章 常微分方程试验  常微分方程数值解方法  向量场雨箭图绘制方法 1/28 一、数值方法求常微分方程   y  f t y  ( , )  y t y  ( 0 ) 0  命令： [T,Y] = ode23(F,Tspan,y0) 这里，Tspan = [t t ]是常微分方程求解区域，y0是初始值， 0 ，N ‘F’ 是包括函数文件名字的符串。 [T,Y]是求解区域内离散数据以及对应的数值解。 步骤： (1)用函数文件定义一阶微分方程(或方程组)右端函数； (2)用MATLAB命令ode23()求数值解或绘积分曲线。 2/28 例7.1 马尔萨斯模型 以1994 年我国人口为12亿为初值，求解 常微分方程。 分析：N(t)表示人口数量,取人口变化率r=0.015,微分方程 dN 18 0.015N dt 16 N (1994) 12 14 12 1990 1995 2000 2005 2010 2015 2020 编辑窗口 function z=fun1(t,N) z=0.015*N; 命令窗口 ode23(fun1,[1994,2020],12) [T,N]=ode23(fun1,[1994,2020],12) 3/28 例7.2 捕食者与被捕食者问题 海岛上有狐狸和野兔,当野兔数量增多时，狐狸捕食野兔导致 狐群数量增长;大量兔子被捕食使狐群进入饥饿状态其数量下 降;狐群数量下降导致兔子被捕食机会减少,兔群数量回升。 微分方程模型如下  dx  x 0.015xy dt   x (0) 100  , dy y 0.01xy y(0) 20    dt 计算 x(t) ，y (t) 当t ∈[0，20]时的数据。绘图并分析捕食者 和被捕食者的数量变化规律。 4/28 % 创建MATLAB的函数文件 function z=fox