第6讲微分方程模型导论.ppt

把形形色色的实际问题化成微分方程的定解问题，大体上可以按以下几步 1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。 2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。 3. 运用这些规律列出方程和定解条件。 数学建模与数学实验 山东工商学院数学与信息科学学院 微分方程模型 涉及“改变”、“变化”、“增加”、“减少”、“衰变”、“边际”、“速度”、 “运动”等词语的确定性连续问题。 微分方程建模的对象 微分方程建模是数学建模的重要方法，因为在很多 实际问题中，经常要涉及各变量的变化率问题，这 些问题的数学描述将导致微分方程模型。所谓微分方 程，就是表示未知函数、未知函数的导数与自变量之 间的关系的方程。未知函数是一元函数的微分方程， 称之为常微分方程；未知函数是多元函数的微分方程， 称之为偏微分方程. 微分方程建模的手段 （1） 按照内在规律或用类比法建立微分方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。 在生物、经济等学科的实际问题中，许多现象的规律性不很清楚，即使有所了解也是极其复杂的，建模时在不同的假设下去模拟实际的现象，建立能近似反映问题的微分方程，然后从数学上求解或分析所建方程及其解的性质，再去同实际情况对比，检验此模型能否刻画、模拟某些实际现象。 （3）模拟近似法 （2） 微元分析法 利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 在实际的微分方程建模过程中，也往往是上述方法的综合应用。不论应用哪种方法，通常要根据实际情况，做出一定的假设与简化，并要把模型的理论或计算结果与实际情况进行对照验证，以修改模型使之更准确地描述实际问题并进而达到预测预报的目的。 称代数方程 f (x)=0 的实根x = x0为方程(9.1)的平衡点(或奇点). 它也是方程的解. 设 常微分方程的稳定性理论 定义1. 如果存在某个邻域，使得微分方程(9.1)的解x(t)从这个邻域的某个x(0)出发，满足 则称 平衡点x0是 渐进稳定的，否则称平衡点x0是不渐进稳定的。 由于 在讨论方程(9.1)的 来代替. 稳定性时，可用 易知 x0也是方程(9.2)的平衡点. (9.2)的通解为 关于x0是否稳定有以下结论： ① 若 则x0是稳定的； ② 若 则x0是不稳定的. 微分方程（组）的Matlab求解 求微分方程（组）的解析解命令: dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’) 结果：u = tan(t+c) 1. 微分方程的解析解 解 输入命令: 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 2. 用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、在解n个未知函数的方程组时，x0和x均为n维向量，m-文件中的待解方程组应以x的分量形式写成. M文件应接受两个参数（尽管m文件中未必用到这两个参数），返回值为一个列向量。