2015年数学建模讲座 4月20日课件.ppt

2015年数学建模讲座 微分方程模型讨论：基本思想、方法与数值模拟分析 贺天兰 微分方程=自然定律 常微分方程：联系着自变量、未知函数、及未知函数的某些导数之间的关系式 总是与变化率相联系：如 速度、加速度、曲率、增长率、出生率、死亡率等 本讲座分两个部分 模型举例  模型求解的基本理论 一、模型举例 体重变化模型(一阶线性常微分方程) 减肥计划模型（差分方程） 国家综合实力模型（一阶二维非线性微分方程组） 模型一、体重变化模型 某人食量10467J/Day,假设其中5038J/Day因新陈代谢而自动消耗，其它活动，如健身等所消耗的热量是69J/(kg.Day)与他体重的乘积。假设与脂肪形式存储的热量100%有效，且1kg脂肪含热量41868J.试研究此人体重随时间变化的规律。 1、问题分析：问题涉及的原则、理论或方法做适当分析 要研究此人体重随时间变化的规律，就是要找出体重随时间变化的函数关系式。 注意各变量与它们的微小增量之间的关系。 2、问题假设：提出有利条件，对实际问题做理想化近似 (1)设t时刻某人的体重为W(t)，一天开始时此人的体重为W(0); (2)W(t)是t的光滑函数； (3)体重在某时间段内的变化等于输入与输出之差，其中，其中输入是指扣除基本新陈代谢之后的净食量吸收，输出是指如健身时的消耗； (4)只考虑一天的情况，且不会出现异常现象。 3.模型建立:实际问题---数学语言:函数\方程\方程组等 此人每天的体重变化=输入-输出 输入热量=10467-5038=5429J 输出热量=69×W=69WJ 4、模型求解 对数学翻译进行求解，既避免了对原问题的直接处理，又利用了现有的数学知识。 体重变化曲线 Clear W0=60; t=0:3000; W=5429/69-(5429-69*w0)/69*exp(-23*t/12956); Plot(t,w,’r’),grid,legend(‘体重变化曲线’)； 5、结果结论：将数学结果返回到实际中做相应的解释 给出明确可信的结论！ 6、敏感性分析：讨论参数的变化对结果的影响 讨论：增加或减少热量的摄入，体重将如何变化？ 增加或减少健身活动，体重又将如何变化？ 那个参数的变化对体重比较敏感？ 定性分析 模型三 综合实力的微分方程模型讨论 一个地区物质文明和精神文明的综合称为综合实力 二、相关的理论知识 微分方程模型讨论：基本思想、方法与数值模拟分析 贺天兰 基本概念：1、微分方程的解、通解 例：y=sin(arcsinx+c) 2、初值问题：柯西问题与特解 3、常微分方程的阶、线性、非线性常微分方程 4、问题：哪些常微分方程可以用初等方法求解？ 自治的常微分方程与非自治的常微分方程 可以用初等方法求解的方程（组） 可分离变量的一阶方程及求解 一阶线性微分方程求解：常数变易法 全微分方程(Total/complete/filly)DE M(x,y)dx+N(x,y)dy=0是全微分方程的充要条件是 一阶二维线性方程组的一般理论 二阶常系数线性微分方程 二阶线性常系数非齐次微分方程 二阶线性常系数非齐次微分方程解的结构 5、黎卡提方程可以求解吗？ 6、一阶方程解的存在性、唯一性 7、定性方法与数值方法 1、一阶微分方程的几何意义； 2、斜率场—两种特例； 解析方法与定性方法相结合；举例 3、数值方法：欧拉方法 一阶的自治与非自治方程 相线：自治方程简化的斜率场 利用相线画解的图像简图 例 画出相线，并研究解的趋势 平衡点的分类：汇、源、结点 具有Allee效应的Logistic模型 分歧 自然条件下海洋中鱼的总量： 8、 一阶二维非线性方程组模型的定性、数值讨论 模型举例 Lotka-Volterra 竞争模型 处理非线性自治常微分方程模型的定性方法 在平衡点处线性化，借助线性模型而得到结论 平衡解、线性化定理 原系统与线性化系统在平衡点处 线性系统与非线性系统在平衡点附近 差不多 function dy=ODEfun6(t,y) dy=zeros(2,1); dy(1)=-3*y(1)+y(2);dy(2)=-y(2); return clear [t1,Y1]=ode45(ODEfun6,[-8:0.1:8],[1;0]); [t2,Y2]=ode45(ODEfun6,[-8:0.1:8],[1;3]); [t3,Y3]=ode45(ODEfun6,[-8:0.1:8],[-1;0]); subplot(3,1,1),plot(t1,Y1,-g);subplot(3,1,2),plot(t2,Y2,-b);subplot(3,1,3),plot(t3,Y3,-r);figure; s