自动控制中的基础数学 微分方程与差分方程(王翼编著)思维导图

常微分方程在数学建模中的应用(免费版)

常微分方程在数学建模中的应用 这里介绍几个典型的用微分方程建立数学模型的例子. 一、人口预测模型 由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去,则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 例1( 马尔萨斯 (Malthus ) 模型) 英国人口统计学家马尔萨斯(1766—1834)在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是:在人口自然增长过程中,净相对增长(出生率与死亡率之差)是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型. 解 设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理(因人口总数很大,可近似地这样处理,此乃离散变量连续化处理),据马尔萨斯的假设,在t 到t t ?+时间段内,人口的增长量为 t t rN t N t t N ?=-?+)()()(, 并设0t t =时刻的人口为0N ,于是 ?????==. , 00)(d d N t N rN t N 这就是马尔萨斯人口模型,用分离变量法易求出其解为 )(00e )(t t r N t N -=, 此式表明人口以指数规律随时间无限增长. 模型检验:据估计1961年地球上的人口总数为9 1006.3?,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3?=N ,02.0=r ,于是 ) 1961(02.09 e 1006.3)(-?=t t N . 这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人 口大约每35年翻一番,而上式断定34.6年增加一倍(请读者证明这一点). 但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地(事实上,地球表面还有80%被水覆盖),我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改. 例2(逻辑Logistic 模型) 马尔萨斯模型为什么不能预测未来的人口呢?这主要是地

数学建模作业、微分方程实验、北京工业大学

2微分方程实验 1、微分方程稳定性分析 绘出下列自治系统相应的轨线,并标出随 t 增加的运动方向,确定平■衡点, 并按稳定的、渐近稳定的、或不稳定的进行分类: 解:(1)由 f (x ) =x=0, f (y ) =y=0;可得平衡点为(0,0), ___ 1 0 系数矩阵A ,求得特征值入1=1,入2=1; 0 1 p=-(入1+入2)=-2<0 , q=入1入2=1>0;对照稳定性的情况表,可知平■衡点(0, 0) 是 不稳定的。 图形如下: (2)如上题可求得平衡点为(0,0 ),特征值入1=-1,入2=2; p=-(入1+入2)=-1<0 , q-入1入2=-2<0;对照稳定性的情况表,可知平■衡点(0, 0) 是 不稳定的。 其图形如下: dx ⑴dt dt x, y; dx dt dy dt dx x, ⑶尸 2y ;晋 dx y , (4) ? 2x;也 dt x+1, 2y.

(3) 如上题可求得平■衡点为(0,0 ),特征值入1=0 + 1.4142i,入2=0 -1.4142i; p=-(入1+入2)= 0, q-入1入2=1.4142>0;对照稳定性的情况表,可知平■衡点(0, 0)是不稳定的。 其图形如下: (4) 如上题可求得平衡点为(1,0 ),特征值入1=-1,入2=-2; p=-(入1+入2)= 3>0, q=入1入2=2>0;对照稳定性的情况表,可知平■衡点(1, 0) 是稳定的。 其图形如下:

2、种群增长模型 一个片子上的一群病菌趋向丁繁殖成一个圆菌落.设病菌的数目为N,单位 成员的增长率为r1,则由Malthus生长律有竺r1 N,但是,处丁周界表面的dt 那些病菌由丁寒冷而受到损伤,它们死亡的数量与N2成比例,其比例系数为r2, 求N满足的微分方程.不用求解,图示其解族.方程是否有平衡解,如果有,是否为稳定的? 解:由题意很容易列出N满足的微分方程:坐r1N r2N; f(N) dt 令f(N)=O,可求得方程的两个平■衡点N1=0,N2=「22/r i2 1 1 d2N 1 5 5 2 (r1 r2N 2) (r1N r2N 2) dt 2 进而求得 A d2N 令r dt 2 2 0可求得N=r2 /4r〔 则N=N1 N=N2 N=r22/4r i2可以把第一象限划为三部分,且从下到上三部分中分 0,冬dt2 .2 2 c dN cdN c dN cdN 0, ;—0, —r 0; —0, ―r dt dt dt dt 则可以画出N (t) 的图形,即微分方程的解族,如下图所示:

常微分方程在数学建模中的应用.

微分方程应用 1 引言 常微分方程的形成与发展和很多学科有着密切的联系,例如力学、天文学、物理学等.数学的其他分支的快速发展,产生出很多新兴学科,这些新兴学科的产生都对常微分方程的发展有着深刻的影响,而且当前计算机的快速发展更是为常微分方程的应用及理论研究提供了非常有力的工具. 数学解决实际问题就必须建立模型,而数学建模就是把数学语言描述实际现象的过程.利用数学去解决各类实际问题时,建立数学模型是十分重要的一步,但是也是最困难的一步.建立数学模型的过程,是把错综复杂的实际问题简化、抽象为合理的数学结构的过程.要通过大量调查、收集相关数据资料,观察和研究实际对象的固有特征和内在规律,抓住问题的主要矛盾,建立起反映实际问题的数量关系,然后利用数学的理论和方法去分析和解决问题. 因此本文先简要介绍了如何建立微分方程模型,并通过具体的实例来简单地介绍了微分方程在数学建模中的应用. 2 数学模型简介 通常我们把现实问题的一个模拟称为模型.如交通图、地质图、航空模型和建筑模型等.利用字母、数学及其它数学符号建立起来的等式或不等式以及图表、图象、框图等来模拟现实的模型称为数学模型.数学模型在实际生活中经常碰到,如求不规则图形的面积,可建立定积分的数学模型,求变化率的问题可建立导数模型,统计学中抽样调查,买彩票中奖的概率问题等等.学会建立数学模型对解决实际生活问题会有很大的帮助. 建立数学模型是沟通摆在面前的实际问题与数学工具之间联系的一座必不可少的桥梁.随着科学技术的进步,特别是电子计算机技术的迅速发展,数学已经渗透到从自然科学技术到工农业生产建设,从经济生活到社会生活的各个领域.一般地说,当实际问题需要我们对所研究的现实对象提供分析、预报、决策、控制等方面的定量结果时,往往都离不开数学的应用,而建立数学模型则是这个过程的关键环节. 3 常微分方程模型 3.1 常微分方程的简介

高等数学第七章微分方程习题

第七章 微分方程与差分方程 习题7-1(A ) 1. 说出下列微分方程的阶数: ;02)()1(2=+'-'x y y y x ;0)2(2=+'+'''y y x y x .0)32()67()3(=++-dy y x dx y x 2. 下列函数是否为该微分方程的解: x e x y y y y 2; 02)1(==+'-'' )(2; 0)()2(2为任意常数C x x C y xdy dx y x -==++ ),(cos sin ; 0) 3(212122 2为任意常数C C ax C ax C y y a dx y d +==+ )(ln ; 02)()4(2xy y y y y y x y x xy =='-'+'+''+ 3. 在下列各题中,确定函数关系式中所含的参数,写出符合初始条件的函数: ;5, )1(0 22==-=x y C y x ;1,0,)()2(0 221=' =+===x x x y y e x C C y . 0,1, )(sin )3(21='=-===ππx x y y C x C y 4. 写出下列条件确定的曲线所满足的微分方程: 点横坐标的平方。 处的切线的斜率等于该曲线在点),()1(y x 轴平分。被,且线段轴的交点为处的法线与曲线上点y PQ Q x y x P ),()2( 习题7-1(B ) 1.在下列各题中,对各已知曲线族(其中 C 1, C 2, C 3 都是任意常数)求出相应的微分方程: ; 1)()1(22=+-y C x . )2(21x x e C e C xy -+= 2.用微分方程表示下列物理问题: 平方成反比。温度的成正比,与的变化率与气压对于温度某种气体的气压P T P )1( 。 速度成反比(比例系数同时阻力与, 成正比(比例系数与时间用在它上面的一个力的质点作直线运动,作一质量为)))2(11k k t m 习题7-2(A ) 1.求下列微分方程的通解: ;0ln )1(=-'y y y x ;0553)2(2='-+y x x ; )()3(2y y a y x y '+='-'

高数公式大全(全)

高数公式大全 1.基本积分表: 三角函数的有理式积分: 2 22212211cos 12sin u du dx x tg u u u x u u x += =+-=+=, , ,  一些初等函数: 两个重要极限: ? ?????????+±+=±+=+=+=+-=?+=?+-==+==C a x x a x dx C shx chxdx C chx shxdx C a a dx a C x ctgxdx x C x dx tgx x C ctgx xdx x dx C tgx xdx x dx x x )ln(ln csc csc sec sec csc sin sec cos 222 22 22 2C a x x a dx C x a x a a x a dx C a x a x a a x dx C a x arctg a x a dx C ctgx x xdx C tgx x xdx C x ctgxdx C x tgxdx +=-+-+=-++-=-+=++-=++=+=+-=????????arcsin ln 21ln 211csc ln csc sec ln sec sin ln cos ln 2 2222222? ????++-=-+-+--=-+++++=+-= ==-C a x a x a x dx x a C a x x a a x x dx a x C a x x a a x x dx a x I n n xdx xdx I n n n n arcsin 22ln 22)ln(221 cos sin 22 2222222 2222222 22 2 22 2 ππx x arthx x x archx x x arshx e e e e chx shx thx e e chx e e shx x x x x x x x x -+=-+±=++=+-==+= -=----11ln 21)1ln(1ln(:2 :2:22)双曲正切双曲余弦双曲正弦...590457182818284.2)11(lim 1sin lim 0==+=∞→→e x x x x x x

数学建模——微分方程的应用

第八节 数学建模——微分方程的应用举例 微分方程在物理学、力学、经济学和管理科学等实际问题中具有广泛的应用,本节我们将集中讨论微分方程的实际应用,尤其是微分方程经济学中的应用. 读者可从中感受到应用数学建模的理论和方法解决实际问题的魅力. 分布图示 ★衰变问题 ★逻辑斯谛方程 ★价格调整问题 ★人才分配问题 内容要点: 一、衰变问题 镭、铀等放射性元素因不断放射出各种射线而逐渐减少其质量, 这种现象称为放射性物质的衰变. 根据实验得知, 衰变速度与现存物质的质量成正比, 求放射性元素在时刻t 的质量. 用x 表示该放射性物质在时刻t 的质量, 则 dt dx 表示x 在时刻t 的衰变速度, 于是“衰变速度与现存的质量成正比”可表示为 .kx dt dx -= (8.1) 这是一个以x 为未知函数的一阶方程, 它就是放射性元素衰变的数学模型, 其中0>k 是比例常数, 称为衰变常数, 因元素的不同而异. 方程右端的负号表示当时间t 增加时, 质量x 减少. 解方程(8.1)得通解.kt Ce x -=若已知当0t t =时, ,0x x =代入通解kt Ce x -=中可得,00kt e x C -= 则可得到方程(8.1)特解 ,)(00t t k e x x --= 它反映了某种放射性元素衰变的规律. 注: 物理学中, 我们称放射性物质从最初的质量到衰变为该质量自身的一半所花费的时间为半衰期, 不同物质的半衰期差别极大. 如铀的普通同位素( U 238)的半衰期约为50亿年;通常的镭( Ra 226)的半衰期是上述放射性物质的特征, 然而半衰期却不依赖于该物质的初始量, 一克Ra 226 衰变成半克所需要的时间与一吨Ra 226衰变成半吨所需要的时间同样都是1600年, 正是这种事实才构成了确定考古发现日期时使用的著名的碳-14测验的基础.

(完整版)高等数学微分方程试题

第十二章 微分方程 §12-1 微分方程的基本概念 一、判断题 1.y=ce x 2(c 的任意常数)是y '=2x 的特解。 ( ) 2.y=(y '')3是二阶微分方程。 ( ) 3.微分方程的通解包含了所有特解。 ( ) 4.若微分方程的解中含有任意常数,则这个解称为通解。 ( ) 5.微分方程的通解中任意常数的个数等于微分方程的阶数。 ( ) 二、填空题 1. 微分方程.(7x-6y)dx+dy=0的阶数是 。 2. 函数y=3sinx-4cosx 微分方程的解。 3. 积分曲线y=(c 1+c 2x)e x 2中满足y x=0=0, y ' x=0=1的曲线是 。 三、选择题 1.下列方程中 是常微分方程 (A )、x 2+y 2=a 2 (B)、 y+0)(arctan =x e dx d (C)、22x a ??+22y a ??=0 (D ) 、y ''=x 2+y 2 2.下列方程中 是二阶微分方程 (A )(y '')+x 2y '+x 2=0 (B) (y ') 2+3x 2y=x 3 (C) y '''+3y ''+y=0 (D)y '-y 2=sinx 3.微分方程2 2dx y d +w 2 y=0的通解是 其中c.c 1.c 2均为任意常数 (A )y=ccoswx (B)y=c sinwx (C)y=c 1coswx+c 2sinwx (D)y=c coswx+c sinwx 4. C 是任意常数,则微分方程y '=3 23y 的一个特解是 (A )y-=(x+2)3 (B)y=x 3+1 (C) y=(x+c)3 (D)y=c(x+1)3 四、试求以下述函数为通解的微分方程。 1.2 2 C Cx y +=(其中C 为任意常数) 2.x x e C e C y 3221+=(其中21,C C 为任意常数) 五、质量为m 的物体自液面上方高为h 处由静止开始自由落下,已知物体在液体中受的阻力与运动的速度成正比。用微分方程表示物体,在液体中运动速度与时间的关系并写出初始条件。

数学建模之微分方程建模与平衡点理论

微分方程 列微分方程常用的方法: (1)根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律来建 立微分方程模型。 (2)微元分析法 利用已知的定理与规律寻找微元之间的关系式,与第一种方法不同的是对 微元而不是直接对函数及其导数应用规律。 (3)模拟近似法 在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有 所了解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建立能 近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性 质,再去同实际情况对比,检验此模型能否刻画、模拟某些实际现象。 一、模型的建立与求解 1.1传染病模型 (1)基础模型 假设:t时刻病人人数() x t连续可微。每天每个病人有效接触(使病人治病的接触)的人数为λ,0 t=时有0x个病人。 +?病人人数增加 建模:t到t t

()()()x t t x t x t t λ+?-=? (1) 0,(0)dx x x x dt λ== (2) 解得: 0()t x t x e λ= (3) 所以,病人人数会随着t 的增加而无限增长,结论不符合实际。 (2)SI 模型 假设:1.疾病传播时期,总人数N 保持不变。人群分为两类,健康者占总人数的比例为s(t),病人占总人数的比例为i(t)。 2.每位病人每天平均有效接触λ人,λ为日接触率。有效接触后健康者变为病人。 依据:患病人数的变化率=Ni(t)(原患病人数)* λs(t)(每个病人每天使健康人变为病人的人数) 建模: di N Nsi dt λ= (4) 由于 ()()1s t i t += (5) 设t=0时刻病人所占的比例为0i ,则可建立Logistic 模型 0(1),(0)di i i i i dt λ=-= (6)

高等数学 微分方程

第十二章 微分方程 § 1 微分方程的基本概念 1、由方程x 2-xy+y 2=C 所确定的函数是方程( )的解。 A. (x-2y)y '=2-xy '=2x-y C.(x-2)dx=(2-xy)dy D.(x-2y)dx=(2x-y)dy 2、曲线族y=Cx+C 2 (C 为任意常数) 所满足的微分方程 ( ) 4.微分方程y '=y x 21-写成以 y 为自变量,x 为函数的形式为( ) A.y x 21dx dy -= B.y x 21dy dx -= '=2x-y D. y '=2x-y §2 可分离变量的微分方程 1.方程P(x,y)dx+Q(x,y)dy=0是( ) A.可分离变量的微分方程 一阶微分方程的对称形式, C.不是微分方程 D.不能变成 ) y ,x (P ) y ,x (Q dy dx -= 2、方程xy '-ylny=0的通解为( ) A y=e x B. y=Ce x cx D.y=e x +C 3、方程满足初始条件:y '=e 2x-y , y|x=0=0的特解为( ) A. e y =e 2x +1 2 1 e ln x 2+= C. y=lne 2x +1-ln2 D. e y =21e 2x +C 4、已知y=y(x)在任一点x 处的增量α+?+=?x x 1y y 2 ,且当?x →0时,α是?x 高阶无穷小,y(0)=π,则y(1)=( ) A. 2π B. π C. 4 e π 4e ππ 5、求特解 cosx sinydy=cosy sinxdx , y|x=0=4 π 解:分离变量为tanydy=tanxdx ,即-ln(cosy)=-ln(cosx)-lnC ,cosy=ccosx 代入初始条件:y|x=0= 4π 得:2 2C =特解为:2cosy=cosx 6、求微分方程()2 y x cos y x 2 1cos dx dy +=-+满足y(0)=π的特解。

最新31微分方程与微分方程建模法汇总

31微分方程与微分方 程建模法

第三章微分方程模型 3.1微分方程与微分方程建模法 一、微分方程知识简介 我们要掌握常微分方程的一些基础知识,对一些可以求解的微分方程及其方程组,要求掌握其解法,并了解一些方程的近似解法。 微分方程的体系:(1)初等积分法(一阶方程及几类可降阶为一阶的方程) ?Skip Record If...?(2)一阶线性微分方程组(常系数线性微分方程组的解法) ?Skip Record If...?(3)高阶线性微分方程(高阶线性常系数微分方程解法)。其中还包括了常微分方程的基本定理。 0.常数变易法:常数变易法在上面的(1)(2)(3)三部分中都出现过,它是由线性齐次方程(一阶或高阶)或方程组的解经常数变易后求相应的非齐次方程或方程组的解的一种方法。 1.初等积分法:掌握变量可分离方程、齐次方程的解法,掌握线性方程的解法,掌握全微分方程(含积分因子)的解法,会一些一阶隐式微分方程的解法(参数法),会几类可以降阶的高阶方程的解法(恰当导数方程)。 分离变量法:(1)可分离变量方程: ?Skip Record If...? (2) 齐次方程:?Skip Record If...? 常数变易法:(1) 线性方程,?Skip Record If...??Skip Record If...?

(2) 伯努里方程,?Skip Record If...??Skip Record If...? 积分因子法:化为全微分方程,按全微分方程求解。 对于一阶隐式微分方程?Skip Record If...?有 参数法:(1) 不含x或y的方程:?Skip Record If...? (2) 可解出x或y的方程:?Skip Record If...? 对于高阶方程,有 降阶法:?Skip Record If...? 恰当导数方程 一阶方程的应用问题(即建模问题)。 2.一阶线性微分方程组:本部分主要内容有:一是一阶线性微分方程组的基本理论(线性齐次、非齐次微分方程组的通解结构,刘维尔公式等),二是常系数线性微分方程组的解法(求特征根,单根与重根[待定系数法]),三是常数变易法。本部分内容与线性代数关系密切,如线性空间,向量的线性相关与线性无关,基与维数,特征方程、特征根与特征向量,矩阵的若当标准型等。3.高阶线性微分方程:了解高阶线性微分方程的基本理论(线性齐次、非齐次微分方程的通解结构,刘维尔公式等); n阶线性常系数微分方程解法:(1)求常系数齐次线性微分方程基本解组的待定指数函数法;(2)求一般非齐次线性方程解的常数变易法;(3)求特

高等数学第七章微分方程试题及答案

第七章 常微分方程 一.变量可分离方程及其推广 1.变量可分离的方程 (1)方程形式: ()()()()0≠=y Q y Q x P dx dy 通解() ()? ?+=C dx x P y Q dy (注:在微分方程求解中,习惯地把不定积分只求出它的一个原函数,而任意常数另外再加) (2)方程形式:()()()()02211=+dy y N x M dx y N x M 通解()()()() C dy y N y N dx x M x M =+??1221 ()()()0,012≠≠y N x M 2.变量可分离方程的推广形式 (1)齐次方程 ?? ? ??=x y f dx dy 令 u x y =, 则()u f dx du x u dx dy =+= ()c x c x dx u u f du +=+=-?? ||ln 二.一阶线性方程及其推广 1.一阶线性齐次方程 ()0=+y x P dx dy 它也是变量可分离方程,通解()?-=dx x P Ce y ,(c 为任意常数) 2.一阶线性非齐次方程 ()()x Q y x P dx dy =+ 用常数变易法可求出通解公式 令()()?-=dx x P e x C y 代入方程求出()x C 则得 ()()()[] ?+=??-C dx e x Q e y dx x P dx x P 3.伯努利方程 ()()()1,0≠=+ααy x Q y x P dx dy 令α-=1y z 把原方程化为()()()()x Q z x P dx dz αα-=-+11 再按照一阶线性非齐次方程求解。 4.方程: ()()x y P y Q dx dy -=1可化为()()y Q x y P dy dx =+ 以y 为自变量,x 为未知函数 再按照一阶线性非齐次方程求解。

高等数学微分方程练习题

(一)微分方程的基本概念 微分方程:含未知函数的导数或微分的方程,称为微分方程、 微分方程的阶:微分方程所含未知函数的最高阶导数或微分的阶数称为微分方程的阶数、 1、不就是一阶微分方程. A、正确 B、不正确 2、不就是一阶微分方程. A、正确 B、不正确 一阶线性微分方程:未知函数及其导数都就是一次的微分方程d ()() d y P x y Q x x +=称为一阶 线性微分方程、 微分方程的解:如果一个函数代入微分方程后,方程两边恒等,则称此函数为微分方程的解、通解:如果微分方程的解中所含独立任意常数C的个数等于微分方程的阶数,则此解称为微分方程的通解、 特解:在通解中根据附加条件确定任意常数C的值而得到的解,称为特解、 1、就是微分方程的解. A、正确 B、不正确 2、就是微分方程的解. A、正确 B、不正确 3、就是微分方程的通解. A、正确 B、不正确 4、微分方程的通解就是( ). A、 B、 C、 D、

(二)变量可分离的微分方程:()()dy f x g y dx = 一阶变量可分离的微分方程的解法就是: (1)分离变量:1221()()()()g y f x dy dx g y f x =;(2)两边积分:1221()()()()g y f x dy dx g y f x =?? 左边对y 积分,右边对x 积分,即可得微分方程通解、 1、微分方程 的通解就是( ). A 、 B 、 C 、 D 、 2、微分方程的通解就是( ). A 、 B 、 C 、 D 、 3、微分方程的通解就是( ). A 、 B 、 C 、 D 、 4、微分方程 的通解就是( ). A 、 B 、 C 、 D 、 5、微分方程 的通解就是( ). A 、 B 、 C 、 D 、 6、微分方程的通解( ). A 、 B 、 C 、 D 、 7、微分方程 的通解就是( ). A 、 B 、 C 、 D 、 8、 x y dy e dx -=就是可分离变量的微分方程. A 、正确 B 、不正确

数学建模之微分方程建模与平衡点理论

微分方程 列微分方程常用的方法: (1)根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律来建立微分方程模型。 (2)微元分析法 利用已知的定理与规律寻找微元之间的关系式,与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 (3)模拟近似法 在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有所了解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建立能近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性质,再去同实际情况对比,检验此模型能否刻画、模拟某些实际现象。 一、模型的建立与求解 1.1传染病模型 (1)基础模型 假设:t 时刻病人人数()x t 连续可微。每天每个病人有效接触(使病人治病的接触)的人数为λ,0t =时有0x 个病人。 建模:t 到t t +?病人人数增加 ()()()x t t x t x t t λ+?-=?(1) 0,(0)dx x x x dt λ==(2) 解得: 0()t x t x e λ=(3) 所以,病人人数会随着t 的增加而无限增长,结论不符合实际。 (2)SI 模型

假设:1.疾病传播时期,总人数N 保持不变。人群分为两类,健康者占总人数的比例为s(t),病人占总人数的比例为i(t)。 2.每位病人每天平均有效接触λ人,λ为日接触率。有效接触后健康者变为病人。 依据:患病人数的变化率=Ni(t)(原患病人数)*λs(t)(每个病人每天使健康人变为病人的人数) 建模: di N Nsi dt λ=(4) 由于 ()()1s t i t +=(5) 设t=0时刻病人所占的比例为0i ,则可建立Logistic 模型 0(1),(0)di i i i i dt λ=-=(6) 解得: 01()111kt i t e i -= ??+- ??? (7) 用Matlab 绘制图1()~i t t ,图2 ~di i dt 图形如下, 结论:在不考虑治愈情况下

第三讲微分方程的理论与数学建模

第三讲 微分方程的理论与数学建模 一、微分方程模型的建立 函数是事物的内部联系在数量方面的反映,如何寻找变量之间的函数关系,在实际应用中具有重要意义。在许多实际问题中,往往不能直接找出变量之间的函数关系,但是根据问题所提供的情况,有时可以列出含有要找的函数及其导数的关系式。这就是所谓的微分方程,从而得出微分方程模型。 例1 物体冷却过程的数学模型 将物体放置于空气中,在时刻0=t 时,测量得它的温度为1500=u C ,10分钟后测量得温度为 C u 1001=。我们要求此物体的温度u 和时间t 的关系,并计算20分钟后物体的温度。这里我们假定 空气温度保持为C u a 24=。 解 为了解决上述问题,需要了解有关热力学的一些基本规律。例如,热量总是从温度高的物体向温度低的物体传导的;在一定的温度范围内,一个物体的温度变化速度与这一物体的温度和其所在介质温度的差值成正比。这是已为实验证实了的牛顿(Newton )冷却定律。 设物体在时刻t 的温度为)(t u u =,则温度的变化速度以 dt du 来表示。注意到热量总是从温度高的物体向温度低的物体传导的,因而a u u >0。所以温度差a u u -恒正;又因物体将随时间而逐渐冷却,故温度变化速度dt du 恒负。故有: dt du )(a u u k --= (1.1) 这里0>k 是比例常数。方程(1.1)就是物体冷却过程的数学模型,它含有未知函数u 及它的一阶导数dt du ,这样的方程称为一阶微分方程。 为了解出物体的温度u 和时间t 的关系,我们要从方程(1.1)中解出u 。注意到a u 是常数,且0>-a u u ,可将(1.1)改写成 kdt u u u u d a a -=--)( (1.2) 这样u 和t 就被分离开了。两边积分,得到 c kt u u a ~)ln(+-=- (1.3) 这里c ~是任意常数。上式可写成 c kt a e u u ~+-=- 令c e c ~=,则有 kt a ce u u -+= (1.4) 再根据初始条件: 当0=t 时,0u u = (1.5) 可得a u u c -=0,于是 kt a a e u u u u --+=)(0 (1.6) 如果k 的数值确定了,(1.6)就完全决定了温度u 和时间t 的关系。 根据条件10=t 时,1u u =,得到 k a a e u u u u 1001)(--+= 由此得到a a u u u u k --=10ln 101051.066.1ln 10 1≈=。从而 t e u 051.012624-+= (1.7)

数学建模作业求解常微分方程和人口模型问题

实验报告 课程名称:数学建模 课题名称:求解常微分方程与人口模型 专业:信息与计算科学 姓名:胡家炜 班级: 123132 完成日期: 2016 年 6 月 10 日

一.求解微分方程的通解 (1). dsolve('2*x^2*y*Dy=y^2+1','x') ans = (exp(C3 - 1/x) - 1)^(1/2) -(exp(C3 - 1/x) - 1)^(1/2) i -i (2). dsolve('Dy=(y+x)/(y-x)','x') ans = x + 2^(1/2)*(x^2 + C12)^(1/2) x - 2^(1/2)*(x^2 + C12)^(1/2) (3). dsolve('Dy=cos(y/x)+y/x','x') ans = (pi*x)/2-x*log(-(exp(C25 + log(x)) - i) /(exp(C25 + log(x))*i - 1))*i (4). dsolve('(x*cos(y)+sin(2*y))*Dy=1','x') ans = -asin(x/2 + lambertw(0, -(C30*exp(- x/2 - 1))/2) + 1) (5). dsolve('D2y+3*Dy-y=exp(x)*cos(2*x)','x') ans = C32*exp(x*(13^(1/2)/2 - 3/2)) + C33*exp(-x*(13^(1/2)/2 + 3/2)) + (13^(1/2)*exp(x*(13^(1/2)/2-3/2))*exp((5*x)/2(13^(1/2)*x)/2)* (2*sin(2*x) - cos(2*x)*(13^(1/2)/2 - 5/2)))/(13*((13^(1/2)/2 - 5/2)^2 +4))-(13^(1/2)*exp(x*(13^(1/2)/2+3/2))*exp((5*x)/2 +(13^(1/2)*x)/2)*(2*sin(2*x)+cos(2*x)*(13^(1/2)/2+5/2))) /(13*((13^(1/2)/2 + 5/2)^2 + 4)) (6)dsolve('D2y+4*y=x+1+sin(x)','x') ans = cos(2*x)*(cos(2*x)/4 - sin(2*x)/8 + sin(3*x)/12 - sin(x)/4 + (x*cos(2*x))/4 - 1/4) + sin(2*x)*(cos(2*x)/8 - cos(3*x)/12 + sin(2*x)/4 + cos(x)/4 + (x*sin(2*x))/4 + 1/8) + C35*cos(2*x) + C36*sin(2*x)

高等数学第九章微分方程试题及答案

第九章 常微分方程 一.变量可分离方程及其推广 1.变量可分离的方程 (1)方程形式: ()()()()0≠=y Q y Q x P dx dy 通解() ()? ?+=C dx x P y Q dy (注:在微分方程求解中,习惯地把不定积分只求出它的一个原函数,而任意 常数另外再加) (2)方程形式:()()()()02211=+dy y N x M dx y N x M 通解()()()() C dy y N y N dx x M x M =+??1221 ()()()0,012≠≠y N x M 2.变量可分离方程的推广形式 (1)齐次方程 ?? ? ??=x y f dx dy 令 u x y =, 则()u f dx du x u dx dy =+= ()c x c x dx u u f du +=+=-?? ||ln 二.一阶线性方程及其推广 1.一阶线性齐次方程 ()0=+y x P dx dy 它也是变量可分离方程,通解()?-=dx x P Ce y ,(c 为任意常数) 2.一阶线性非齐次方程 ()()x Q y x P dx dy =+ 用常数变易法可求出通解公式 令()()?-=dx x P e x C y 代入方程求出()x C 则得 ()()()[] ?+=??-C dx e x Q e y dx x P dx x P 3.伯努利方程 ()()()1,0≠=+ααy x Q y x P dx dy 令α-=1y z 把原方程化为()()()()x Q z x P dx dz αα-=-+11 再按照一阶线性非齐次方程求解。 4.方程: ()()x y P y Q dx dy -=1可化为()()y Q x y P dy dx =+ 以y 为自变量,x 为未知函数 再按照一阶线性非齐次方程求解。 三、可降阶的高阶微分方程

常微分方程在数学建模中的应用论文正稿

毕业论文 论文题目:常微分方程在数学建模中的应用姓名: 学科专业: 指导教师: 完成时间:

常微分方程是数学理论(特别是微积分)联系实际的重要工具,它不仅与几何学、力学、电子技术、自动控制、星际航行、甚至和化学、生物学、农业以及经济学都有着密切的联系。本文结合实践背景,建立数学模型,并利用所得结果去解释某些实际问题。 关键字常微分方程、人口预测模型、市场价格模型、混合溶液的数学模型、震动模型

第一章人口预测模型 第二章市场价格模型 第三章混合溶液的数学模型第四章震动模型

绪论 当我们描述实际对象的某些特性随时间(或空间)而演变的过程、分析它的变化规律、预测它的未来性态,研究它的控制手段时,通常要建立对象的动态模型。建模时首先要根据建模目的和对问题的具体分析作出简化假设,然后按照对象在的或可以类比的其他对象的规律列出微分方程,求出方程的解并将结果翻译回实际对象,就可以进行描述、分析、预测或控制了。 事实上在微分方程课程中,解所谓应用题时我们遇到简单的建立动态模型问题,例如“一质量为m的物体自高h处自由下落,初速度是零,设阻力与下落速度的平方成正比,比例系数为k,求下落速度随时间的变化规律。”又如“容器有盐水100L,含盐10kg,令以3L/min的速度从一管放进净水,以2L/min的速度从另一管抽出盐水,设容器盐水浓度始终是均匀的,求容器含盐量随时间变化规律。”本文讨论的是常微分方程在数学建模中的应用。

第一章 人口预测模型 由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去,则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 例1(马尔萨斯(Malthus )模型) 英国人口统计学家马尔萨斯(1766—1834)在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是:在人口自然增长过程中,净相对增长(出生率与死亡率之差)是常数,即单位时间人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型. 解 设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理(因人口总数很大,可近似地这样处理,此乃离散变量连续化处理),据马尔萨斯的假设,在t 到t t ?+时间段,人口的增长量为 t t rN t N t t N ?=-?+)()()(, 并设0t t =时刻的人口为0N ,于是 ?????==. ,00)(d d N t N rN t N 这就是马尔萨斯人口模型,用分离变量法易求出其解为 )(00e )(t t r N t N -=, 此式表明人口以指数规律随时间无限增长. 模型检验:据估计1961年地球上的人口总数为91006.3?,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3?=N ,02.0=r ,于是 )1961(02.09e 1006.3)(-?=t t N . 这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间

常微分方程在数学建模中的应用

北方民族大学学士学位论文 论文题目:常微分方程在数学建模中的应用 院(部)名称:信息与计算科学学院 学生姓名:马木沙 专业:信计学号:20093490 指导教师姓名:魏波 论文提交时间: 论文答辩时间: 学位授予时间: 北方民族大学教务处制

摘要 本文利用常微分方程和数学建模二者之间的联系,了解微分方程的一般理论、微分方程解的存在惟一性、微分方程的稳定性问题、通过几个典型的数学模型如:人口模型、减肥的数学模型、化工车间通风模型、传染病的传播模型及定性分析等例子来体现微分方程在数学建模中的应用. 用数学理论解决实际生活中的问题.微分方程的出现以及运用微分方程在数学建模中的应用,就是为了更好地使更多的人理解并运用数学理论,更好的解决实际生活中的问题.努力在各个领域利用并渗透数学知识的广泛运用. 关键词:常微分方程,数学建模,数学模型

Abstract In this paper, ordinary differential equations and mathematical modeling contact between the two, understand the general theory of differential equations, stability problems of the existence and uniqueness of differential equations, differential equations, several typical mathematical models such as: demographic model,example of the mathematical model of weight loss, chemical plant ventilation model, spread of infectious diseases, model and qualitative analysis to reflect the application of differential equations in mathematical modeling. found that the application of mathematical theory to study and solve problems in the actual process of the emergence of ordinary differential equations andOrdinary Differential Equations in Mathematical Modeling widely used, in order to better enable ordinary people to understand and use mathematical theory, solving real-world problems. sublimation theory by the knowledge-based transformation to the ability to type, highlight the differential equationsand differential equations in mathematical modeling efforts made outstanding and significant contribution in various fields. Keywords: ordinary differential equations, mathematical modeling, mathematical model.

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

实验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

相关文档
最新文档