北京科技大学数学模型M05微分模型.pptVIP

在19世纪中叶,通过刘维尔的工作,人们已经知道绝大多数的微分方程不能用初等积分方法求解.这个结果对于微分方程理论的发展产生了极大影响,使微分方程的研究发生了一个转折.既然初等积分法有着不可克服的局限性,那么是否可以不求微分方程的解,而是从微分方程本身来推断其解的性质呢?定性理论和稳定性理论正是在这种背景下发展起来的.前者由法国数学家庞加莱(Poincaré,1854-1912)在19世纪80年代所创立，后者由俄国数学家李雅普罗夫(Liapunov,1857-1918)在同年代所创立.它们共同的特点就是在不求出方程的解的情况下,直接根据微分方程本身的结构和特点,来研究其解的性质.由于这种方法的有效性,近一百多年以来它们已经成为常微分方程发展的主流.李雅普诺夫创立了处理稳定性问题的两种方法:第一方法要利用微分方程的级数解,在他之后没有得到大的发展;第二方法是在不求方程解的情况下,借助一个所谓的李雅普诺夫函数 和通过微分方程所计算出来的导数的符号性质,就能直接推断出解的稳定性,因此又称为直接法. * 第五章 微分方程模型 第五章 微分方程模型 5.1 传染病模型 5.2 烟雾的扩散与消失 在研究实际问题时，常常会联系到某些变量的变化率或导数， 这样所得到变量之间的关系式就是微分方模型。微分方程模型反映的是变量之间的间接关系，因此，要得到直接关系，就得求微分方程。 微分方程建模 根据建模目的和问题分析作出简化假设 （1）根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。 （2）微元分析法 利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 例1 求 的通解. 求解微分方程有三种方法： 1）求精确解；2）求数值解（近似解）；3）定性 理论方法。 微分方程的解析解 求微分方程（组）的解析解命令: dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’) 结 果：u = tg(t-c) 解 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x) 结 果 为 : y =3e-2xsin（5x） 解 输入命令 ： [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 [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：分别为设定的相对误差和绝对误差. 求解常微分方程 [t,y]=ode23(‘fun’,tspan,y0) 需要把微分方程化为一阶微分方程组 tspan=[t0 tf]表示积分的起始值和终止值； fun是定义函数的文件名； [t,y]=ode45(‘fun’,tspan,y0) y0是初始状态列向量。 常 微 分 方 程 的 数 值 解 法 解: 令 y1=x，y2=y1’ 1、建立m-文件vdp1000.m如下： function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1); 2、取t0=0，tf=3000，输入命令： [T,Y]=ode15s(vdp1000,[0 3000],[2 0]); plot(T,Y(:,1),-) 3、结果如图 解 1、建立m-文件rigid.m如下： function d