微分方程模型new.pptVIP

微分方程的基本概念 含有未知函数以及未知函数的导数或微分的方程 常微分方程求解 方程特点：能够解析求解 方程特点：方程比较复杂，不能够得到解析解。 例2：求 在命令窗口中输入命令: [x,y]=ode45(myfun2,[0,12],[0,1,1]); plot(x,y(:,1),*) 小结： 求解ODE的基本过程： （1）根据问题所属学科中的规律、定律、公式，用常 微分方程与初始条件进行描述； （2）运用数学中的变量替换，把高阶的微分方程（组） 写成一阶微分方程组，初始条件也作相应的改变； （3）根据（1）与（2）的结果 编写M-函数文件odefile； （4）将odefile文件与初始条件传递给求解器solver中的 一个，运行后就可得到ODE的在指定区间上的解 列向量（其中包含y及不同阶的倒数）。 令 作业：枯井深度的估计 有一口枯井和一只跑表，用扔下一块石头听回声的方法来估计枯井的深度。假定你能准确的测定时间，你有应当怎样推算古井的深度？请你分析这一问题。 微分方程模型 1 微分方程基本概念 2 微分方程的求解 3 微分方程模型——传染病模型 4 常见微分方程模型 微分方程中出现的未知函数的导数或微分的最高阶数 常微分方程（ODE）：未知函数是一元函数的微分方程； 偏微分方程（PDE）：未知函数是多元函数的微分方程； 微分方程 微分方程的阶 微分方程的分类 常微分方程符号解 （MATLAB）符号解语法:(在命令窗口输入) dsolve(equation,condition) 其中：equation代表常微分方程，方程中 Condition为初始条件。 例1:计算微分方程 在初始条件 下的特解。 dsolve(x*Dy+2*y=exp(x), y(1)=2*exp) 常微分方程数值解 例如，非线性微分方程等，一般要求得到的解在若干点上满足规定精确度的近似值，或得到一个满足精度要求的便于计算的表达式。 （MATLAB）数值解语法： [x,y]=solver(fun,xspan,y0) 求解器solver 有:ode45;ode23;ode113;ode23t等，常 用ode45;ode23; fun代表编写的函数文件名； xspan代表求解区间； y0代表初始条件。 在初始条件 下的数值解。 分析：令 则微分方程变为微分方程组 解：建立M文件myfun2.m function dy=myfun2(x,y) dy=zeros(3,1); dy(1)=y(2); dy(2)=y(3); dy(3)=-(51/100)*y(1)*y(2); 导数的作用 在连续问题中表示瞬时变化率 在离散问题中逼近平均变化率 根据函数及其变化率之间的关系确定函数 微分方程建模 根据建模目的和问题分析作出简化假设 按照内在规律或用类比法建立微分方程 传染病模型 问题 描述传染病的传播过程 分析受感染人数的变化规律 预报传染病高潮到来的时刻 预防传染病蔓延的手段 按照传播过程的一般规律，用机理分析方法建立模型 t 时刻已感染者(病人)所占比例为 y(t). 每个病人每天有效接触(足以使人致病)人数为?，且设此疾病既不导致死亡也不会康复（流传初期） 模型1 假设 若有效接触的是病人，则不能使病人数增加 必须区分已感染者(病人)和未感染者(健康人) 建模 ? 设某地区总人数为N,最初时刻病人所占比例为y0 ， 模型2 区分已感染者(病人)和未感染者(健康人) 假设 1）总人数N不变，病人和健康 人的 比例分别为 2）每个病人每天有效接触人数为?, 且使接触的健康人致病。 建模 ? ~ 日 接触率 SI 模型 每个病人每天可使? 个健康者变为病人。 因为病人数为N 个， 时间内增加的病人数为： 可表示为 模型2 1/2 tm y y0 1 0 t tm~传染病高潮到来时刻 ? (日接触率)? ? tm? Logistic 模型 病人可以治愈！ ? t=tm, dy/dt 最大 模型3 传染病无免疫性——病人治愈成为健康人，健康人可再次被感染 增加假设 SIS 模型 3）病人每天治愈的占病人总数的比例为? ?