国赛前第二次讲座资料-1.pptVIP

2014年数学建模讲座 微分方程模型讨论：基本思想、方法与数值模拟分析 贺天兰 微分方程=自然定律 常微分方程：联系着自变量、未知函数、及未知函数的某些导数之间的关系式 总是与变化率相联系：如 速度、加速度、曲率、增长率、出生率、死亡率等 模型举例 体重变化模型 国家综合实力模型 模型一、体重变化模型 某人食量10467J/Day,假设其中5038J/Day因新陈代谢而自动消耗，其它活动，如健身等所消耗的热量是69J/(kg.Day)与他体重的乘积。假设与脂肪形式存储的热量100%有效，且1kg脂肪含热量41868J.试研究此人体重随时间变化的规律。 一、问题分析：问题涉及的原则、理论或方法做适当分析 要研究此人体重随时间变化的规律，就是要找出体重随时间变化的函数关系式。 注意各变量与它们的微小增量之间的关系。 二、问题假设：提出有利条件，对实际问题做理想化近似 (1)设t时刻某人的体重为W(t)，一天开始时此人的体重为W(0); (2)W(t)是t的光滑函数； (3)体重在某时间段内的变化等于输入与输出之差，其中，其中输入是指扣除基本新陈代谢之后的净食量吸收，输出是指如健身时的消耗； (4)只考虑一天的情况，且不会出现异常现象。 三.模型建立:实际问题---数学语言:函数\方程\方程组等 此人每天的体重变化=输入-输出 输入热量=10467-5038=5429J 输出热量=69×W=69WJ 四、模型求解 对数学翻译进行求解，既避免了对原问题的直接处理，又利用了现有的数学知识。 体重变化曲线 Clear W0=60; t=0:3000; W=5429/69-(5429-69*w0)/69*exp(-23*t/12956); Plot(t,w,’r’),grid,legend(‘体重变化曲线’)； 五、结果结论：将数学结果返回到实际中做相应的解释 给出明确可信的结论！ 六、敏感性分析：讨论参数的变化对结果的影响 讨论：增加或减少热量的摄入，体重将如何变化？ 增加或减少健身活动，体重又将如何变化？ 那个参数的变化对体重比较敏感？ 定性分析 模型二 基于周期扰动的综合实力的 微分方程模型讨论 相关的理论知识 微分方程模型讨论：基本思想、方法与数值模拟分析 贺天兰 基本概念：1、微分方程的解、通解 例：y=sin(arcsinx+c) 2、初值问题：柯西问题与特解 3、常微分方程的阶、线性、非线性常微分方程 混合问题：溶液箱（溶质与溶液的混合，如水中溶解盐，且有流入流出） 已知在时刻t=0时，x(0)=x0,计算在时刻t的盐 总量x(t). 假设盐浓度为c1g/L溶液以不变的速率r1L/s流 入箱内，并且箱内的溶液始终通过搅拌来混 合，以不变的速率r0L/s流出。 混合问题的模型 人口增长 放射性物质的衰减：放射性碳定年代 关键：任何生物体中碳原子的比例是不变的 注意：校正因子表 药效的消失 指数增长方程 药效消失举例：麻醉狗 在狗的每公斤体重中，狗的血液中至少含有 45mg的药量，狗才会昏迷。假设这种药的药性在狗的血液中5h消失一半，问为了麻醉一条 50kg的狗1h,需要用药多少? 5、问题：哪些常微分方程可以用初等方法求解？ 自治的常微分方程与非自治的常微分方程 可以用初等方法求解的方程（组） 可分离变量的一阶方程及求解 一阶线性微分方程求解：常数变易法 全微分方程(Total/complete/filly)DE M(x,y)dx+N(x,y)dy=0是全微分方程的充要条件是 6、黎卡提方程可以求解吗？ 7、一阶方程解的存在性、唯一性 8、定性方法与数值方法 1、一阶微分方程的几何意义； 2、斜率场—两种特例； 解析方法与定性方法相结合；举例 3、数值方法：欧拉方法 一阶的自治与非自治方程 相线：自治方程简化的斜率场 利用相线画解的图像简图 例 画出相线，并研究解的趋势 平衡点的分类：汇、源、结点 具有Allee效应的Logistic模型 分歧 自然条件下海洋中鱼的总量： 9、一阶二维线性方程组的一般理论 10、二阶线性常系数微分方程 方程形式 二阶线性常系数非齐次微分方程 二阶线性常系数非齐次微分方程解的结构 11、 一阶二维非线性方程组模型的定性、数值讨论 模型举例 Lotka-Volterra 竞争模型 注意： 数学模型中，等式的两端必须有相同的量纲---物理量的属性 无量纲化---定性讨论---数值讨论 7个基本物理量：长度L 、质量M 、时间T 、电流I 、热力学温度Q 、物质的量N 、发光强度J 个体的数目：S q的量纲为dim(q) A=B ,A+B,A-B: dim(A)=dim(B) 1/A: dim(A)-1 A*B: dim(A