微分方程模型2

实验07讲评参考答案_微分方程模型(2学时)

实验07讲评、参考答案 讲评未按时交的同学 批改情况：

附参考答案： 实验07 微分方程模型（2学时） （第5章 微分方程模型） 1.（验证）传染病模型2（SI 模型）p136~138 传染病模型2（SI 模型）： 0(1),(0)di k i i i i dt =-= 其中， i (t )是第t 天病人在总人数中所占的比例。 k 是每个病人每天有效接触的平均人数（日接触率）。 i 0是初始时刻（t =0）病人的比例。 1.1 画~di i dt 曲线图p136~138 取k =0.1，画出 i dt di ~的曲线图，求i 为何值时dt di 达到最大值，并在曲线图上标注。 %传染病模型2（SI 模型）的di/dt~i 曲线图 %文件名：p137fig2.m %λ=0.1 clear; clc; 《数学建模实验》 王平

提示：fplot, fminbnd, plot, text, title, xlabel 1）画曲线图 用fplot函数，调用格式如下： fplot(fun,lims) fun必须为一个M文件的函数名或对变量x的可执行字符串。 若lims取[xmin xmax]，则x轴被限制在此区间上。 若lims取[xmin xmax ymin ymax]，则y轴也被限制。 本题可用 fplot('0.1*x*(1-x)',[0 1.1 0 0.03]); 2）求最大值 用求解边界约束条件下的非线性最小化函数fminbnd，调用格式如下： x=fminbnd('fun',x1,x2) fun必须为一个M文件的函数名或对变量x的可执行字符串。 返回自变量x在区间x1

微分方程建模案例

第五章微分方程建模案例 微分方程作为数学科学的中心学科，已经有三百多年的发展历史，其解法和理论已日臻完善，可以为分析和求得方程的解（或数值解）提供足够的方法，使得微分方程模型具有极大的普遍性、有效性和非常丰富的数学涵。微分方程建模包括常微分方程建模、偏微分方程建模、差分方程建模及其各种类型的方程组建模。微分方程建模对于许多实际问题的解决是一种极有效的数学手段，对于现实世界的变化，人们关注的往往是其变化速度、加速度以及所处位置随时间的发展规律，其规律一般可以用微分方程或方程组表示，微分方程建模适用的领域比较广，涉及到生活中的诸多行业，其中的连续模型适用于常微分方程和偏微分方程及其方程组建模，离散模型适用于差分方程及其方程组建模。本章主要介绍几个简单的用微分方程建立的模型，让读者一窥方程的应用。下面简要介绍利用方程知识建立数学模型的几种方法： 1．利用题目本身给出的或隐含的等量关系建立微分方程模型 这就需要我们仔细分析题目，明确题意，找出其中的等量关系，建立数学模 型。 例如在光学里面，旋转抛物面能将放在焦点处的光源经镜面反射后成为平行光线，为了证明具有这一性质的曲线只有抛物线，我们就是利用了题目中隐含的条件——入射角等于反射角来建立微分方程模型的。 2．从一些已知的基本定律或基本公式出发建立微分方程模型

我们要熟悉一些常用的基本定律、基本公式。例如从几何观点看，曲线 y y(x)上某点的切线斜率即函数y y(x)在该点的导数；力学中的牛顿第二运 动定律：F ma ，其中加速度a 就是位移对时间的二阶导数，也是速度对时间 的一阶导数等等。从这些知识出发我们可以建立相应的微分方程模型。 例如在动力学中，如何保证高空跳伞者的安全问题。对于高空下落的物体， 我们可以利用牛顿第二运动定律建立其微分方程模型， 设物体质量为m ，空气阻 力 系数为k ，在速度不太大的情况下，空气阻力近似与速度的平方成正比；设时 刻t 时物体的下落速度为v ，初始条件：v (o ) 0.由牛顿第二运动定律建立其微 分方程模型: 求解模型可得: 体在地面上的投影面积。根据极限速度求解式子，在m,, 一定时，要求落地速 度w 不是很大时，我们可以确定出s 来，从而设计出保证跳伞者安全的降落伞的 直径大小来 3?利用导数的定义建立微分方程模型 dv m 一 dt mg kv 2 ? k(exp[2t 由上式可知，当t 其中，阻力系数k 1) 时，物体具有极限速度: lim v t mg :k ， s ， 为与物体形状有关的常数, 为介质密度，s 为物 、mg(exp[2t 1)

第二章动力学系统的微分方程模型

第二章：动力学系统的微分方程模型 利用计算机进行仿真时，一般情况下要给出系统的数学模型，因此有必要掌握一定的建立数学模型的方法。在动力学系统中，大多数情况下可以使用微分方程来表示系统的动态特性，也可以通过微分方程可以将原来的系统简化为状态方程或者差分方程模型等。在这一章中，重点介绍建系统动态问题的微分方程的基本理论和方法。 在实际工程中，一般把系统分为两种类型，一是连续系统；其数学模型一般是高阶微分方程；另一种是离散系统，它的数学模型是差分方程。 §2.1 动力学系统统基本元件 任何机械系统都是由机械元件组成的，在机械系统中有3种类型的基本机械元件：惯性元件、弹性元件和阻尼元件。 1 惯性元件：惯性元件是指具有质量或转动惯量的元件，惯量可以定义为使加速度（或角加速度）产生单位变化所需要的力（或力矩）。 惯量（质量）= ） 加速度（力（2 /) s m N 惯量（转动惯量）= ） 角加速度（力矩（2/) s rad m N ? 2 弹性元件：它在外力或外力偶作用下可以产生变形的元件，这种元件可以通过外力做功来储存能量。按变形性质可以分为线性元件和非线性元件，通常等效成一弹簧来表示。 对于线性弹簧元件，弹簧中所受到的力与位移成正比，比例常数为弹簧刚度k 。 x k F ?= 这里k 称为弹簧刚度，x ?是弹簧相对于原长的变形量，弹性力的方向总是指向弹 簧的原长位移，出了弹簧和受力之间是线性关系以外，还有所谓硬弹簧和软弹簧，它们的受力和弹簧变形之间的关系是一非线性关系。 3 阻尼元件：这种元件是以吸收能量以其它形式消耗能量，而不储存能量，可以形象的表示为一个活塞在一个充满流体介质的油缸中运动。阻尼力通常表示为： α x c R = 阻尼力的方向总是速度方向相反。当1=α，为线性阻尼模型。否则为非线性阻 尼模型。应注意当α等于偶数情况时，要将阻尼力表示为： ||1--=αx x c R 这里的“-”表示与速度方向相反

第三章-微分方程模型

微分方程模型 1.1微分方程模型简介 对于现实世界的变化，人们关注的往往是变量之间的变化率，或者变化速度、加速度以 及所处的位置随时间的发展规律，之中的规律一般可以写成一个（偏）微分方程或方程组。 所以实际问题中，有大批的问题可以用微分方程来建立数学模型，涉及的领域包括物理学、 化学、天文学、生物学、力学、政治、经济、军事、人口、资源等等。微分方程建模是数学 建模的重要方法，因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色 色的实际问题化成微分方程的定解问题，大体上可以按以下几步： 1?、根据实际要求确定要研究的量（自变量、未知函数、必要的参数等）并确定坐标系； 2?、找出这些量所满足的基本规律（物理的、几何的、化学的或生物学的等等）； 3?、运用这些规律列出方程和定解条件。 2.1微分方程模型运用实例 例1:发射卫星为什么用三级火箭 采用运载火箭把人造卫星发射到高空轨道上运行，为什么不能用一级火箭而必须用多级 火箭系统？ 下面通过建立运载火箭有关的数学模型来回答上述问题。 火箭是一个复杂的系统，为了使问题简单明了，我们只从动力系统和整体结构上分析， 并且假设引擎是足够强大的。 首先解决第一个问题：为什么不能用一级火箭发射人造卫星，下面用三个数学模型回答 这个问题： （1 ）卫星进入600km高空轨道时，火箭必须的最低速度。 首先将问题理想化，假设： （i）卫星轨道是以地球中心为圆心的某个平面上的圆周，卫星在此轨道上以地球引力作为向心力绕地球作平面匀速圆周运动； （ii ）地球是固定于空间中的一个均匀球体，其质量集中于球心； iii）其它星球对卫星的引力忽略不计。 建模与求解：设地球半径为R，质量为M ;卫星轨道半径为r，卫星质量为m。 根据假设（"）和（iii），卫星只受到地球的引力，由牛顿万有引力定律可知其引力大小为 GMm F— (1) r 其中G为引力常数。 为消去常数G，把卫星放在地球表面，则由（1）式得 GMm 亠m2 mg 2 或GM 二R g R 再代入（1）式，得

3.1 微分方程模型的建模步骤

第3章微分方程模型 3.1 微分方程模型的建模步骤 在自然科学以及工程、经济、医学、体育、生物、社会等学科中的许多系统，有时很难找到该系统有关变量之间的直接关系——函数表达式，但却容易找到这些变量和它们的微小增量或变化率之间的关系式，这时往往采用微分关系式来描述该系统——即建立微分方程模型。我们以一个例子来说明建立微分方程模型的基本步骤。 例1 某人的食量是10467（焦/天），其中5038（焦/天）用于基本的新陈代谢（即自动消耗）。在健身训练中，他所消耗的热量大约是69（焦/公斤?天）乘以他的体重（公斤）。假设以脂肪形式贮藏的热量100%地有效，而1公斤脂肪含热量41868（焦）。试研究此人的体重随时间变化的规律。 模型分析 在问题中并未出现“变化率”、“导数”这样的关键词，但要寻找的是体重（记为W ）关于时间t 的 函数。如果我们把体重W 看作是时间t 的连续可微函数，我们就能找到一个含有的dt dW 微分方程。 模型假设 1.以)(t W 表示t 时刻某人的体重，并设一天开始时人的体重为0W 。 2．体重的变化是一个渐变的过程。因此可认为 )(t W 是关于t 连续而且充分光滑的。 3．体重的变化等于输入与输出之差，其中输入是指扣除了基本新陈代谢之后的净食量吸收；输出就是进行健身训练时的消耗。 模型建立 问题中所涉及的时间仅仅是“每天”，由此，对于“每天” 体重的变化=输入-输出。 由于考虑的是体重随时间的变化情况，因此，可得 体重的变化/天=输入/天—输出/天。 代入具体的数值，得 输入/天 = 10467（焦/天）—5038（焦/天）=5429（焦/天）， 输出/天 = 69（焦/公斤?天）×W （公斤）= 69W （焦/天）。 体重的变化/天=t W ??（公斤/天）dt dW t =→?0 考虑单位的匹配，利用 “公斤/天=公斤焦天 焦/41868 /”, 可建立如下微分方程模型

第五章微分方程模型

第五章 微分方程模型 、 某人每天由饮食获取10467焦热量，其中5038焦用于新陈代谢，此外每公斤体重需支付69焦热量作为运动消耗，其余热量则转化为脂肪，已知以脂肪形式贮存的热量利用率为100%，每公斤脂肪含热量41868焦，问此人的体重如何随时间而变化 解： 设此人的体重为w ，则根据题意有，每天获取的热量，减去新陈代谢，减去运动消耗的热量，剩余的按利用率100% 转化为脂肪，即有下列等式成立： 1046750386941868 w dw dt --= 经化简有： 232313956139565429()41868t t w e t e c - =-?+ 假设此人现在的体重为0w ，则此人的体重随时间的变化如下： 2323139561395605429()41868t t w e t e w - =-?+ 、 生活在阿拉斯加海滨的鲑鱼服从Malthus 增长模型)(003.0)(t p dt t dp = 其中t 以分钟计。在0=t 时一群鲨鱼来到此水域定居，开始捕食鲑鱼。鲨鱼捕杀鲑鱼的速率是)(001.02t p ，其中)(t p 是t 时刻鲑鱼总数。此外，由于在它们周围出现意外情况，平均每分钟有条鲑鱼离开此水域。 （1）考虑到两种因素，试修正Malthus 模型。 （2）假设在0=t 是存在100万条鲑鱼，试求鲑鱼总数 )(t p ，并问∞→t 时会发生什么情况 解：

（1），由题可知， 在考虑两种因素后，修正后的Malthus 模型如下： 2()0.003()0.001()0.002dp t p t p t dt =-- （2），假设在0t = 时，存在100万条鲑鱼，即(0)1000000p = ，解下列初值问题 2()0.003()0.001()0.002(0)1000000 dp t p t p t dt p ?=--???=? 解得 0.0010.0012999998()11000001t t ae p t a ae --+==-其中 当t →∞ 时，2p →。 、 根据罗瑟福的放射性衰变定律，放射性物质衰变的速度与现存的放射性物质的原子数成正比，比例系数成为衰变系数，试建立放射性物质衰变的数学模型。若已知某放射性物质经时间21T 放射物质的原子下降至原来的一半（21T 称为该物质的半衰期）试决定其衰变系数。 解： 假设初始时刻该放射性物质的原子数位0N ，在时间t 时，该放射性物质的原子个数为N ，设衰变系数为k ，则有下列微分方程： 0,(0)dN kN N N dt =-= 解得 0()kt N t N e =

微分方程模型

数学建模学习辅导 第三章 微分方程模型 本章重点： 车间空气清洁问题、减肥问题、单种群增长问题与多物种相互作用问题等数学模型的建立过程与所使用的方法 复习要求： 1．进一步理解建模基本方法与基本建模过程，掌握平衡原理与微元法在建模中的用法． 所谓平衡原理是指自然界的任何物质在其变化的过程中一定受到某种平衡关系的支配.注意发掘实际问题中的平衡原理是从物质运动机理的角度组建数学模型的一个关键问题.就象中学的数学应用题中等量关系的发现是建立方程的关键一样. 微元法是指在组建对象随着时间或空间连续变化的动态模型时，经常考虑它在时间或空间的微小单元变化情况，这是因为在这些微元上的平衡关系比较简单，而且容易使用微分学的手段进行处理．这类模型基本上是以微分方程的形式给出的． 例1 设警方对司机饮酒后驾车时血液中酒精含量的规定为不超过80%(mg/ml). 现有一起交通事故,在事故发生3个小时后,测得司机血液中酒精含量是56%(mg/ml), 又过两个小时后, 测得其酒精含量降为40%(mg/ml),试判断: 事故发生时,司机是否违反了酒精含量的规定? 解：模型建立 设)(t x 为时刻t 的血液中酒精的浓度, 则依平衡原理时间间隔],[t t t ?+内, 酒精浓度的改变量 t t x x ??∝?)(, 即 t t kx t x t t x ?-=-?+)()()( 其中k >0为比例常数, 式前负号表示浓度随时间的推移是递减的, 遍除以t ?, 并令0→?t , 则得到 ,d d kx t x -= 且满足40)5(,56)3(==x x 以及0)0(x x =. 模型求解 容易求得通解为kt c t x -=e )(, 代入0)0(x x =,得到 kt x t x -=e )(0 则)0(0x x =为所求. 又由,40)5(,56)3(==x x 代入0)0(x x =可得 17.04056e 40e 56e 25030=?=????==--k x x k k k 将17.0=k 代入得 25.93e 5656e 17.03017 .030≈?=?=??-x x >80

常微分课后答案解析第二章

第一章 绪论 §1、1 微分方程:某些物理过程的数学模型 §1、2 基本概念 习题1、2 1.指出下面微分方程的阶数,并回答方程就是否线性的: (1) y x dx dy -=24; (2)0122 2 2=+??? ??-xy dx dy dx y d ; (3)0322 =-+?? ? ??y dx dy x dx dy ; (4)x xy dx dy dx y d x sin 352 2=+-; (5) 02cos =++x y dx dy ; (6)x e dx y d y =+??? ? ??22sin . 解 (1)一阶线性微分方程; (2)二阶非线性微分方程; (3)一阶非线性微分方程; (4)二阶线性微分方程; (5)一阶非线性微分方程; (6)二阶非线性微分方程. 2.试验证下面函数均为方程022 2=+y dx y d ω的解,这里0>ω就是常数. (1)x y ωcos =; (2)11(cos C x C y ω=就是任意常数); (3)x y ωsin =; (4)22(sin C x C y ω=就是任意常数); (5)2121,(sin cos C C x C x C y ωω+=就是任意常数); (6)B A B x A y ,()sin(+=ω就是任意常数). 解 (1)y x dx y d x dx dy 2 222cos ,sin ωωωωω-=-=-=,所以022 2=+y dx y d ω,故

x y ωcos =为方程的解. (2)y x C y x C y 2 2 11cos , sin ωωωωω-=-=''-=',所以022 2=+y dx y d ω,故x C y ωcos 1=为方程的解. (3)y x dx y d x dx dy 2222sin ,cos ωωωωω-=-==,所以02 2 2=+y dx y d ω,故x y ωsin =为方程的解. (4)y x C y x C y 2 2 22sin , cos ωωωωω-=-=''=',所以022 2=+y dx y d ω,故x C y ωsin 2=为方程的解. (5)y x C x C y x C x C y 2222121sin cos , cos sin ωωωωωωωωω-=--=''+-=',所 以02 2 2=+y dx y d ω,故x C x C y ωωsin cos 21+=为方程的解. (6)y B x A y B x A y 2 2 )sin(, )cos(ωωωωω-=+-=''+=',故0222=+y dx y d ω,因 此)sin(B x A y +=ω为方程的解. 3.验证下列各函数就是相应微分方程的解: (1)x x y sin = ,x y y x cos =+'; (2)212x C y -+=,x xy y x 2)1(2 =+'-(C 就是任意常数); (3)x Ce y =,02=+'-''y y y (C 就是任意常数); (4)x e y =,x x x e ye y e y 2212-=-+'-; (5)x y sin =,0cos sin sin 22 2 =-+-+'x x x y y y ; (6)x y 1- =,12 22++='xy y x y x ; (7)12 +=x y ,x y x y y 2)1(2 2 ++-='; (8))()(x f x g y = ,) () ()()(2x f x g y x g x f y '-'='.

常微分课后答案解析第二章

第一章 绪论 §1.1 微分方程：某些物理过程的数学模型 §1.2 基本概念 习题1.2 1．指出下面微分方程的阶数，并回答方程是否线性的： （1） y x dx dy -=24； （2）0122 2 2=+??? ??-xy dx dy dx y d ； （3）0322 =-+? ? ? ??y dx dy x dx dy ； （4）x xy dx dy dx y d x sin 352 2=+-； （5） 02cos =++x y dx dy ； （6）x e dx y d y =+??? ? ??22sin ． 解 （1）一阶线性微分方程； （2）二阶非线性微分方程； （3）一阶非线性微分方程； （4）二阶线性微分方程； （5）一阶非线性微分方程； （6）二阶非线性微分方程． 2．试验证下面函数均为方程02 2 2=+y dx y d ω的解，这里0>ω是常数． （1）x y ωcos =； （2）11(cos C x C y ω=是任意常数)； （3）x y ωsin =； （4）22(sin C x C y ω=是任意常数)； （5）2121,(sin cos C C x C x C y ωω+=是任意常数)； （6）B A B x A y ,()sin(+=ω是任意常数)．

解 （1）y x dx y d x dx dy 2222cos ,sin ωωωωω-=-=-=，所以02 2 2=+y dx y d ω，故x y ωcos =为方程的解． （2）y x C y x C y 2 2 11cos , sin ωωωωω-=-=''-='，所以0222=+y dx y d ω，故 x C y ωcos 1=为方程的解． （3）y x dx y d x dx dy 2 222sin ,cos ωωωωω-=-==，所以022 2=+y dx y d ω，故x y ωsin =为方程的解． （4）y x C y x C y 2 2 22sin , cos ωωωωω-=-=''='，所以022 2=+y dx y d ω，故x C y ωsin 2=为方程的解． （5）y x C x C y x C x C y 2222121sin cos , cos sin ωωωωωωωωω-=--=''+-='， 所以022 2=+y dx y d ω，故x C x C y ωωsin cos 21+=为方程的解． （6）y B x A y B x A y 2 2 )sin(, )cos(ωωωωω-=+-=''+='，故02 2 2=+y dx y d ω，因此)sin(B x A y +=ω为方程的解． 3．验证下列各函数是相应微分方程的解： （1）x x y sin = ，x y y x cos =+'； （2）212x C y -+=，x xy y x 2)1(2 =+'-（C 是任意常数）； （3）x Ce y =，02=+'-''y y y （C 是任意常数）； （4）x e y =，x x x e ye y e y 2212-=-+'-； （5）x y sin =，0cos sin sin 22 2 =-+-+'x x x y y y ； （6）x y 1- =，12 22++='xy y x y x ； （7）12 +=x y ，x y x y y 2)1(2 2 ++-='；

微分方程模型建模实例

微分方程模型建模实例 1．一个半球状雪堆，其体积融化的速率与半球面面积S成正比，比例系数k > 0。设融化中雪堆始终保持半球状，初始半径为R且3小时中融化了总体积的7/8，问雪堆全部融化还需要多长时间？ 2．从致冰厂购买了一块立方体的冰块，在运输途中发现，第一小时大约融化了1/4 （1）求冰块全部融化要多长时间（设气温不变） （2）如运输时间需要2.5小时，问：运输途中冰块大约会融化掉多少？ 3．一展开角为α的圆锥形漏斗内盛着高度为H的水，设漏斗底部的孔足够大（表面张力不计），试求漏斗中的水流光需要多少时间？ 4．容器甲的温度为60度，将其内的温度计移入容器乙内，设十分钟后温度计读数为70度，又过十分钟后温度计读数为76度，试求容器乙内的温度。 5．一块加过热的金属块初始时比室温高70度，20分钟测得它比室温高60度，问：（1）2小时后金属块比室温高多少？（2）多少时间后，金属块比室温高10度？ 6．设初始时容器里盛放着含净盐10千克的盐水100升，现对其以每分钟3升的速率注入清水，容器内装有搅拌器能将溶液迅时搅拌均匀，并同时以每分钟2升的速率放出盐水，求1小时后容器里的盐水中还含有多少净盐？ 7．某伞降兵跳伞时的总质量为100公斤（含武器装备），降落伞张开前的空气阻力为0.5v,该伞降兵的初始下落速度为0，经8秒钟后降落伞打开，降落 伞打开后的空气阻力约为0.6 试球给伞降兵下落的速度v(t),并求其下落的极限速度。 8． 1988年8月5日英国人Mike McCarthy创建了一项最低开伞的跳伞纪录，它从比萨斜塔上跳下，到离地179英尺时才打开降落伞，试求他落地时的速度。 9．证明对数螺线r=A 上任一处的切线与极径的夹角的正切为一常 数，（）

第二章 微 分 方 程 模 型.

第二章 微 分 方 程 模 型 建立微分方程模型就是把物理、化学、生物科学、工程科学和社会科学中的规律和原理用含有待定函数的导数或微分的数学关系式表示出来。这一章我们由浅入深地介绍一些微分方程模型。 2.1 简单模型 例1 物体在空气中的下落与特技跳伞问题 假设质量为m 的物体在空气中下落，空气阻力与物体的速度平方成正比，阻尼系数为k （>0）,求物体的运动规律。 解 所谓运动规律即下落距离与时间的关系，如图2.1.1, 建立坐标系。设x 为物体下落的距离，于是物体下落的速度为 dx v dt =， 加速度为 22d x a dt =， 根据牛顿第二定律F ma =，可以列出微分方程 2 22d x d x m k m g d t d t ?? =-+ ???, （2.1.1） 负号表示阻力方向与速度方向相反。 例2 单摆的自由振动问题。 如图2.1.2 为一个单摆，上端固定在O 点，M 为一质量为m 的质点，摆杆OM 之长为L （摆杆的质量忽略不计）。单摆的平衡位置为铅垂线'OO 。将质点M 拉开，使OM 与'OO 成一个角度0θ，然后放手任其自由运动，试求摆杆OM 和铅垂线'OO 的夹角θ与时间t 的关系。 解 将重力分解为径向力F 与切向力T ，T 的大小为sin mg θ，M 的切向加速 度为22d a L dt θ =，于是，由牛顿第二定律，列出微分方程 22s i n d m a m L m g dt θ θ== ， 即 22s i n d g dt L θθ=-， （2.1.2）

设初始时刻0t =，摆杆的初始位置为0θ，初始角速度为0，则单摆的运动规律的研究就化为微分方程的初值问题 ()()22 00' 0s i n ,,0.t t d g dt L t t θθθθθ==?=-??? =??=??? （2.1.3） 图2.1.1 图2.1.2 例3 考古和地质学中文物和化石年代的测定问题。 考古、地质学等方面的专家常用14C （碳14）来估计文物或化石的年代。它们的依据是，宇宙射线不断轰击大气层，使之产生中子，中子与氧气作用生成具有放射性的14C 。这种放射性碳可以氧化成二氧化碳。二氧化碳被植物所吸收，而动物又以植物为食物，于是放射性碳就被带到各种动植物体内。由于14C 是放射性的，无论存在于空气中或生物体内它都在不断衰变，活着的生物通过新陈代谢不断地摄取14C ，使得生物体内的14C 与空气中的14C 有相同的百分含量。生物体死后它停止摄取14C ，因而尸体内的14C 由于不断衰变而不断减少。碳定年代法就是根据14C 的衰变减少量的变化情况来判定生物的死亡时间的。 基本假设 （1）现代生物体中14C 的衰变速度与古代生物体中14C 的衰变速度相同（依据是地球周围大气中14C 的百分含量可认为基本不变，即宇宙射线照射大气层的强度自古至今基本不变）； （2）14C 的衰变速度与该时刻14C 的含量成正比（这条假设的根据来自于原子物理学理论）。 下面用微分方程建模。 设在时刻t （年）生物体中14C 的存量为()x t ，由假设（2）知

数学建模实验答案微分方程模型

实验07 微分方程模型（2学时） （第5章 微分方程模型） 1.（验证）传染病模型2（SI 模型）p136~138 传染病模型2（SI 模型）： 0(1),(0)di k i i i i dt =-= 其中， i (t )是第t 天病人在总人数中所占的比例。 k 是每个病人每天有效接触的平均人数（日接触率）。 i 0是初始时刻（t =0）病人的比例。 1.1 画~di i dt 曲线图p136~138 取k =0.1，画出i dt di ~的曲线图，求i 为何值时dt di 达到最大值，并在曲线图上标注。 参考程序：

提示：fplot, fminbnd, plot, text, title, xlabel 1）画曲线图 用fplot函数，调用格式如下： fplot(fun,lims) fun必须为一个M文件的函数名或对变量x的可执行字符串。 若lims取[xmin xmax]，则x轴被限制在此区间上。 若lims取[xmin xmax ymin ymax]，则y轴也被限制。 本题可用 fplot('0.1*x*(1-x)',[0 1.1 0 0.03]); 2）求最大值 用求解边界约束条件下的非线性最小化函数fminbnd，调用格式如下：x=fminbnd('fun',x1,x2) fun必须为一个M文件的函数名或对变量x的可执行字符串。 返回自变量x在区间x1

第五章----微分方程模型

第五章 微分方程模型 5.1、 某人每天由饮食获取10467焦热量，其中5038焦用于新陈代谢，此外每公斤体重需支付69焦热量作为运动消耗，其余热量则转化为脂肪，已知以脂肪形式贮存的热量利用率为100%，每公斤脂肪含热量41868焦，问此人的体重如何随时间而变化？ 解： 设此人的体重为w ，则根据题意有，每天获取的热量，减去新陈代谢， 减去运动消耗的热量，剩余的按利用率100% 转化为脂肪，即有下列等式成立： 1046750386941868 w dw dt --= 经化简有： 232313956139565429()41868t t w e t e c -=-?+ 假设此人现在的体重为0w ，则此人的体重随时间的变化如下： 2323139561395605429()41868t t w e t e w - =-?+ 5.2、 生活在阿拉斯加海滨的鲑鱼服从Malthus 增长模型)(003.0)(t p dt t dp = 其中t 以分钟计。在0=t 时一群鲨鱼来到此水域定居，开始捕食鲑鱼。鲨鱼捕杀鲑鱼的速率是)(001.02t p ，其中)(t p 是t 时刻鲑鱼总数。此外，由于在它们周围出现意外情况，平均每分钟有0.002条鲑鱼离开此水域。 （1）考虑到两种因素，试修正Malthus 模型。 （2）假设在0=t 是存在100万条鲑鱼，试求鲑鱼总数 )(t p ，并问∞→t 时会发生什么情况？ 解： （1），由题可知， 在考虑两种因素后，修正后的Malthus 模型如下： 2()0.003()0.001()0.002dp t p t p t dt =-- （2），假设在0t = 时，存在100万条鲑鱼，即(0)1000000p = ，解下列初值问题

第三章,微分方程模型

第三章 微分方程建模 在许多实际问题的研究中，要直接导出变量之间的函数关系较为困难，但要导出包含未知函数的导数或微分的关系式却较为容易，此时即可用建立微分方程模型的方法来研究实际问题。例如，根据自由落体运动的重力加速度g 为常数及初始条件即可得出自由落体运动的公式、根据单摆的受力分析及牛顿第二定理即可得到单摆运动满足的方程等等就是典型的实例。本章除了介绍一些来自经典力学的物理及一些几何方面的微分方程问题以外，也介绍了一些稍有不同的微分方程应用题。这些模型研究的主要是来自于非物理领域的实际问题，对这些问题，我们将分析其特征，根据具体情况进行类比，提出假设条件并建立微分方程模型加以研究。提出的假设条件不同，将会导出不同的微分方程。最后还要将求解的结果与实际现象进行对比，如果差异较大还应反复修改假设建立新的模型。因此，在这类模型中，微分方程被当成了研究问题的工具。事实上，在连续变量问题的研究中，微分方程或微分方程组还是十分常用的数学工具之一。 §3.1 几个简单实例 例3.1 （理想单摆运动的周期） 本例的目的是建立理想单摆运动满足的微分方程，由该微分方程即可得出理想单摆运动的周期公式。 （图3-1） 从图3-1中不难看出，小球所受的合力为 θsin mg ，根据牛顿第二定律可得： θθ sin mg ml -= 从而得出两阶微分方程： sin 0 (0)0,(0)g l θθθθθ?+=???==? （3.1） 这就是理想单摆运动满足的微分方程。 （3.1）是一个两阶非线性常微分方程，不容易求解。根据微积分知识，当θ很小时，有sin θ≈θ，此时，为简单起见，我们可考察（3.1）的近似线性方程：

微分方程模型2(2课时)

第2章 微分方程模型2（2课时） 教学目的 懂得如何根据实际情况建立微分方程。 教学内容 介绍应用微分方程方法建模。 模型Ⅱ 人口模型 1. 引言 在研究某些实际问题时，经常无法直接得到各变量之间的联系，问题的特性往往会给出关于变化率的一些关系。利用这些关系，我们可以建立相应的微分方程模型。在自然界以及工程技术领域中，微分方程模型是大量存在的。它甚至可以渗透到人口问题以及商业预测等领域中去，其影响是广泛的。 问题：人口问题是当今世界上最令人关注的问题之一，一些发展中国家的人口出生率过高，越来越威胁着人类的正常生活，有些发达国家的自然增长率趋于零，甚至变为负数，造成劳动力紧缺，也是不容忽视的问题。另外，在科学技术和生产力飞速发展的推动下，世界人口以空前的规模增长，统计数据显示： 年 1625 1830 1930 1960 1974 1987 1999 人口（亿） 5 10 20 30 40 50 60 可以看出，人口每增长十亿的时间，由一百年缩短为十二三年。我们赖以生存的地球，已经带着它的60亿子民踏入了21世纪。 长期以来，人类的繁衍一直在自发地进行着。只是由于人口数量的迅速膨胀和环境质量的急剧恶化，人们才猛然醒悟，开始研究人类和自然的关系，人口数量的变化规律，以及如何进行人口控制等问题。 我国是世界第一人口大国，地球上每九个人中就有一个中国人，在20世纪的一段时间内我国人口的增长速度过快，如下表： 年 1908 1933 1953 1964 1982 1990 2000 人口（亿） 3.0 4.7 6.0 7.2 10.3 11.3 12.95 有效地控制人口的增长，不仅是使我国全面进入小康社会、到21世纪中叶建

第5章 微分方程模型

第5章 微分方程模型 一、讨论题 1. 某人每天由饮食获取10467焦热量，其中5038焦用于新陈代谢，此外每公斤体重需支付69焦热量作为运动消耗，其余热量则转化为脂肪，已知以脂肪形式贮存的热量利用率为100%，每公斤脂肪含热量41868焦，问此人的体重如何随时间而变化？ 2. 根据罗瑟福的放射性衰变定律，放射性物质衰变的速度与现存的放射性物质的原子数成正比，比例系数成为衰变系数，试建立放射性物质衰变的数学模型。若已知某放射性物质经时间21T 放射物质的原子下降至原来的一半（21T 称为该物质的半衰期）试决定其衰变系数。 3. 建立耐用消费品市场销售量的模型。如果已知了过去若干时期销售量的情况，如何确定模型的参数。 4. 两棵不同类别的植物种在一起，按比例吸取养料，试建立它们的生长模型。 二、思考题 1、具有什么样特征的数学建模问题需要用微分方程方法建立模型？ 2、用微分方程方法建立数学模型的基本步骤是什么？应注意哪些问题？ 3、某种细菌的增长率不知道，但假设它是常数，试验开始时估计大约有11500个细菌，一时后有2000个，问四时后大约有多少细菌？ 三、习题 1. 生活在阿拉斯加海滨的鲑鱼服从Malthus 增长模型 )(003.0) (t p dt t dp = 其中t 以分钟计。在0=t 时一群鲨鱼来到此水域定居，开始捕食鲑鱼。鲨鱼捕杀鲑鱼的速 率是)(001.02 t p ，其中)(t p 是t 时刻鲑鱼总数。此外，由于在它们周围出现意外情况，平 均每分钟有0.002条鲑鱼离开此水域。 （1）考虑到两种因素，试修正Malthus 模型。 （2）假设在0=t 是存在100万条鲑鱼，试求鲑鱼总数 )(t p ，并问∞→t 时会发生什么情况？ 2. 用具有放射性的14 C 测量古生物年代的原理是：宇宙线轰击大气层产生中子，中子 与氮结合产生14C 。植物吸收二氧化碳时吸收了14C ，动物食用植物从植物中得到14 C 。在 活组织中14C 的吸收速率恰好与14 C 的衰变速率平衡。但一旦动植物死亡，它就停止吸收14C ，于是14C 的浓度随衰变而降低。由于宇宙线轰击大气层的速度可视为常数，既动物刚 死亡时14 C 的衰变速率与现在取的活组织样本（刚死亡）的衰变速率是相同的。若测得古生 物标本现在14C 的衰变速率，由于14 C 的衰变系数已知，即可决定古生物的死亡时间。试建 立用14C 测古生物年代的模型（14 C 的半衰期为5568年）。 3. 试用上题建立的数学模型，确定下述古迹的年代： （1）1950年从法国Lascaux 古洞中取出的碳测得放射性计数率为0.97计数（min ?g ），而活树木样本测得的计数为6.68计数（min ?g ），试确定该洞中绘画的年代； （2）1950年从某古巴比伦城市的屋梁中取得碳标本测得计数率为4.09计数（min ?g ），活数标本为6.68计数（min ?g ），试估计该建筑的年代。 4. 一容器用一薄膜分成容积为A V 和B V 的两部分，分别装入同一物质不同浓度的溶液。

微分方程模型(二)

微分方程模型 1．一个半球状雪堆，其体积融化的速率与半球面面积S成正比，比例系数k > 0。设融化中雪堆始终保持半球状，初始半径为R且3小时中融化了总体积的7/8，问雪堆全部融化还需要多长时间？ 2．从致冰厂购买了一块立方体的冰块，在运输途中发现，第一小时大约融化了1/4 （1）求冰块全部融化要多长时间（设气温不变） （2）如运输时间需要2.5小时，问：运输途中冰块大约会融化掉多少？ 3．一展开角为α的圆锥形漏斗内盛着高度为H的水，设漏斗底部的孔足够大（表面张力不计），试求漏斗中的水流光需要多少时间？ 4．容器甲的温度为60度，将其内的温度计移入容器乙内，设十分钟后温度计读数为70度，又过十分钟后温度计读数为76度，试求容器乙内的温度。 5．一块加过热的金属块初始时比室温高70度，20分钟测得它比室温高60度，问：（1）2小时后金属块比室温高多少？（2）多少时间后，金属块比室温高10度？ 6．设初始时容器里盛放着含净盐10千克的盐水100升，现对其以每分钟3升的速率注入清水，容器内装有搅拌器能将溶液迅时搅拌均匀，并同时以每分钟2升的速率放出盐水，求1小时后容器里的盐水中还含有多少净盐？ 7．某伞降兵跳伞时的总质量为100公斤（含武器装备），降落伞张开前的空气阻力为0.5v,该伞降兵的初始下落速度为0，经8秒钟后降落伞打开，降落伞打 开后的空气阻力约为0.6试球给伞降兵下落的速度v(t),并求其下落的极限速度。 8． 1988年8月5日英国人Mike McCarthy创建了一项最低开伞的跳伞纪录，它从比萨斜塔上跳下，到离地179英尺时才打开降落伞，试求他落地时的速度。 9．证明对数螺线r=A上任一处的切线与极径的夹角的正切为一常数， （） 10．实验证明，当速度远低于音速时，空气阻力正比与速度，阻力系数大约为0.005。现有一包裹从离地150米高的飞机上落下，（1）求其落地时的速度（2）如果飞机高度更大些，结果会如何，包裹的速度会随高度而任意增大吗？11．生态学家估计人的内禀增长率约为0.029，已知1961年世界人口数为30.6 亿（3.06×）而当时的人口增长率则为0.02。试根据Logistic模型计算：

扩散问题的偏微分方程模型_数学建模

实验一SDH网元基本配置 一、实验目的： 通过本实验，了解SDH光传输的原理和系统组成，了解ZXMP S325设备的硬件构成和单板功能，学习ZXONM 300 网管软件的使用方法，掌握SDH 网元配置的基本操作。 二、实验器材： 1、SDH 设备：3 套ZXMP 325； 2、实验用维护终端。 三、实验原理 1、SDH 原理 同步数字体制（SDH）是为高速同步通信网络制定的一个国际标准，其基础在于直接同步复用。按照SDH 组建的网络是一个高度统一的、标准化的、智能化的网络，采用全球统一的接口以实现多环境的兼容，管理操作协调一致，组网与业务调度灵活方便，并且具有网络自愈功能，能够传输所有常见的支路信号，应用于多种领域（如光纤传输，微波和卫星传输等）。 SDH 具有以下特点： （1）接口：接口的规范化是设备互联的关键。SDH对网络节点接口（NNI）作了统一的规范，内容包括数字信号数率等级、帧结构、复接方法、线路接口、监控管理等。 电接口：STM-1是SDH的第一个等级，又叫基本同步传送模块，比特率为155.520Mb/s；STM-N是SDH第N个等级的同步传送模块，比特率是STM-1的N 倍(N=4n=1，4，16，- - -)。 光接口：采用国际统一标准规范。SDH仅对电信号扰码，光口信号码型是加扰的NRZ 码，信号数率与SDH 电口标准信号数率相一致。 （2）复用方式 a）低速SDH----高速SDH，字节间插； b) 低速PDH-----SDH，同步复用和灵活的映射。 （3）运行维护:用于运行维护（OAM）的开销多，OAM功能强——这也是线路编码不用加冗余的原因. （4）兼容性:SDH 具有很强的兼容性，可传送PDH 业务，异步转移模式信号（ATM）及其他体制的信号。 （5）SDH 复用映射示意图

[数学]第二章微分方程模型

第二章微分方程管模型 第一节人口学模型 人口问题是当今世界人类面临的五大问题的首要问题。我国是世界上人口最多的国家，由于20世纪五六十年代人口政策方面的失误，不仅造成人口总数增长过快，而且年龄结构也不合理。人口的过份增长给我国经济发展造成沉重袍袱，严重地影响经济建设。能否有效地控制人口的增长，已成为本世纪中叶我国能否达到中等发达国家以至赶上发达国家的关键。 建立数学模型对人口发展过程进行描述、分析和预测，并进而研究控制人口增长的生育政策，已引起有关人口专家和官员的极大关注和兴趣。以下我们就如何建立人口数学模型作简要的介绍。 一. 马尔萨斯人口增长模型 对于如何预测人口的增长，早在8世纪人们就开始进行了。英国早期经济学家马尔萨斯(1766-1843)他在担任牧师期间，根据教堂100多年人口出生统计资料，他发现人口出生率是一个常数。于是在1798年发表的《人口原理》一书中提出了哄动于世的马尔萨斯人口模型：假设x(t)表示t时刻人口总数，r为人口增长率(常数)，其他影响人口增长的因素均不考虑。则在t到t+△t这段时间内人口总数增长为 ()()() +?-=? x t t x t rx t t

两端同除以t ?，并令0t ?→，得 ()00dx rx dt x t x ?=???=? ……… (1) 其解为 ()() 00r t t x t x e -= (2) 方程(1)称为马尔萨斯模型，它的解是一个以r e 为公比的几何级数。马尔萨斯根据这个模型认为人口的增长是按几何级数增加，而物质的增长只能按算术级数增加。因此，人口必须加以控制。不幸的是马尔萨斯的这一忠告没有引起人们的足够重视。 马尔萨斯模型 (1) 和 (2) 与19世纪以前欧州一些地区的人口统计数据十分吻合。据估计1961年全世界人口总数为93.0610?,而在此之前的十来年间人口按年2％的速率增长。因此有 9001961, 3.0610,0.02t x r ==?= 于是 ()()0.02196193.0610t x t e -=?? (3) 根据统计资料，在1700—1961年间世界人口大约每35年增加一倍。由(3)式可以算出每34.6年人口增加一倍。 事实上，设在0T t t =-内地球上人口增加一倍，即当0t t =时， 90.0209 0 3.0610 3.0610x e ?=??=?,当0 T t t =-时，90.0202 3.0610,T x e =?所以