《微分差分方程模型幻灯片.ppt》由会员分享，可在线阅读，更多相关《微分差分方程模型幻灯片.ppt（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、微分差分方程模型第1页，共45页，编辑于2022年，星期六 在研究实际问题时在研究实际问题时,我们常常不能直接得出变量之间的我们常常不能直接得出变量之间的关系关系,但却能容易得出包含变量导数在内的关系式但却能容易得出包含变量导数在内的关系式,这就是微这就是微分方程分方程.在现实社会中在现实社会中,又有许多变量是离散变化的又有许多变量是离散变化的,如人口数、生如人口数、生产周期与商品价格等产周期与商品价格等,而且离散的运算具有可操作性而且离散的运算具有可操作性,差分差分正是联系连续与离散变量的一座桥梁正是联系连续与离散变量的一座桥梁.不管是微分方程还是差分方程模型，有时无法得到其解析不管是微分方
2、程还是差分方程模型，有时无法得到其解析解解(必要时必要时,可以利用计算机求其数值解可以利用计算机求其数值解),),既使得到其解析解既使得到其解析解,尚有未知参数需要估计尚有未知参数需要估计(这是可利用第二章参数估计方法这是可利用第二章参数估计方法).).而在实际问题中而在实际问题中,讨论问题的解的变化趋势很重要，因讨论问题的解的变化趋势很重要，因此，以下只对其平衡点的稳定性加以讨论此，以下只对其平衡点的稳定性加以讨论.第2页，共45页，编辑于2022年，星期六3.1 微分方程模型微分方程模型 如果如果则称平衡点则称平衡点x0是是稳定稳定的的.称代数方程称代数方程 f(x)=0 的实根的实根x=
3、x0为方程为方程(3-1)的的平衡点平衡点(或奇点或奇点).它也是方程它也是方程(3-1)的解的解.设设第3页，共45页，编辑于2022年，星期六稳定性判别方法稳定性判别方法由于由于在讨论方程在讨论方程(3-1)的的来代替来代替.稳定性时，可用稳定性时，可用 易知易知 x0也是方程也是方程(3-2)的平衡点的平衡点.(3-2)的通解为的通解为关于关于x0是否稳定有以下结论：是否稳定有以下结论：若若则则x0是稳定的；是稳定的；若若则则x0是不稳定的是不稳定的.这个结论对这个结论对于于(4-1)也是也是成立的成立的.第4页，共45页，编辑于2022年，星期六 关于常微分方程组的平衡点及其稳定性关于
4、常微分方程组的平衡点及其稳定性,设设代数方程组代数方程组的实根的实根x=x0,y=y0称为方程称为方程(3-3)的的平衡点平衡点,记作记作P0(x0,y0).它也是方程它也是方程(3-3)的解的解.第5页，共45页，编辑于2022年，星期六如果如果则称平衡点则称平衡点P0是是稳定稳定的的.下面给出判别平衡点下面给出判别平衡点P0是否稳定的是否稳定的判别准判别准则则.设设 则当则当p0且且q0时时,平衡点平衡点P0是稳定的；当是稳定的；当p0或或q0时时,平衡点平衡点P0是不稳定的是不稳定的.第6页，共45页，编辑于2022年，星期六3.2 3.2 差分方程模型差分方程模型 对于对于k阶差分方程
5、阶差分方程F(n;xn,xn+1,xn+k)=0 (3-6)若有若有xn=x(n),满足满足F(n;x(n),x(n+1),x(n+k)=0,则称则称xn=x(n)是差分方程是差分方程(3-6)的的解解,包含个任意常数包含个任意常数的解称为的解称为(3-6)的的通解通解,x0,x1,xk-1为已知时称为为已知时称为(3-6)的的初始条件初始条件,通解中的任意常数都由初始条件确定后的通解中的任意常数都由初始条件确定后的解称为解称为(3-6)的的特解特解.若若x0,x1,xk-1已知已知,则形如则形如xn+k=g(n;xn,xn+1,xn+k-1)的差分方程的解可以在计算机上实现的差分方程的解可以
6、在计算机上实现.第7页，共45页，编辑于2022年，星期六 若有常数若有常数a是差分方程是差分方程(3-6)的解的解,即即F(n;a,a,a)=0,则称则称 a是差分方程是差分方程(3-6)的的平衡点平衡点.又对差分方程又对差分方程(3-6)的任意由初始条件确定的解的任意由初始条件确定的解 xn=x(n)都有都有xna(n),则称这个平衡点则称这个平衡点a是是稳定稳定的的.一阶常系数线性差分方程一阶常系数线性差分方程 xn+1+axn=b,(其中其中a,b为常数为常数,且且a-1,0)的通解为的通解为xn=C(-a)n+b/(a+1)易知易知b/(a+1)是其平衡点是其平衡点,由上式知由上式知
7、,当且仅当当且仅当|a|1时时,b/(a+1)是稳定的平衡点是稳定的平衡点.第8页，共45页，编辑于2022年，星期六 二阶常系数线性差分方程二阶常系数线性差分方程xn+2+axn+1+bxn=r,其中其中a,b,r为常数为常数.当当r=0时时,它有一特解它有一特解x*=0；当当r 0,且且a+b+1 0时时,它有一特解它有一特解x*=r/(a+b+1).不管是哪种情形不管是哪种情形,x*是其平衡点是其平衡点.设其特征方程设其特征方程 2+a +b=0的两个根分别为的两个根分别为 =1,=2.第9页，共45页，编辑于2022年，星期六 当当 1,2是两个不同实根时是两个不同实根时,二阶常系数线
8、性差二阶常系数线性差分分方程的通解为方程的通解为xn=x*+C1(1)n+C2(2)n;当当 1,2=是两个相同实根时是两个相同实根时,二阶常系数线性二阶常系数线性差分差分方程的通解为方程的通解为xn=x*+(C1+C2 n)n;当当 1,2=(cos +i sin )是一对共轭复根时是一对共轭复根时,二阶常系数线性差分二阶常系数线性差分方程的通解为方程的通解为xn=x*+n(C1cosn +C2sinn ).易知易知,当且仅当特征方程的任一特征根当且仅当特征方程的任一特征根|i|1时时,平衡点平衡点x*是稳定的是稳定的.则则第10页，共45页，编辑于2022年，星期六对于一阶非线性差分方程对
9、于一阶非线性差分方程xn+1=f(xn)其平衡点其平衡点x*由代数方程由代数方程x=f(x)解出解出.为分析平衡点为分析平衡点x*的稳定性的稳定性,将上述差分方程近似将上述差分方程近似为一阶常系数线性差分方程为一阶常系数线性差分方程时时,上述近似线性差分方程与上述近似线性差分方程与原原非线性差分方程的非线性差分方程的稳定性相同稳定性相同.因此因此当当时时,x*是稳定的；是稳定的；当当时时,x*是不稳定的是不稳定的.当当第11页，共45页，编辑于2022年，星期六3.3 3.3 观众厅地面设计观众厅地面设计1 问题的提出在影视厅或报告厅，经常会为前边观众遮挡住自己的视线而苦恼。显然，场内的观众都
10、在朝台上看，如果场内地面不做成前低后高的坡度模式，那么前边观众必然会遮挡后面观众的视线。试建立数学模型设计良好的报告厅地面坡度曲线。第12页，共45页，编辑于2022年，星期六建立坐标系oo处在台上的设计视点bb第一排观众的眼睛到x轴的垂 直距离xyadda第一排观众与设计视点的水平距离d相邻两排的排距视线升高标准x表示任一排与设计视点的水平距离求任一排x与设计视点o的竖直距离函数使此曲线满足视线的无遮挡要求。问题第13页，共45页，编辑于2022年，星期六2 问题的假设(建模初等变量假设)1)观众厅地面的纵剖面图一致，只需求中轴线上地面的起伏曲线即可。2)同一排的座位在同一等高线上。3)每个
11、坐在座位上的观众的眼睛与地面的距离相等。4)每个坐在座位上的观众的头与地面的距离也相等（等变量假设）。5)所求曲线只要使观众的视线从紧邻的前一个座位的人的头顶擦过即可。第14页，共45页，编辑于2022年，星期六3 建模设眼睛升起曲线应满足微分方程初始条件obxyadd1）从第一排起，观众眼睛与o点的连线的斜率随排数的增加而增加，而眼睛升起曲线显然与这些直线皆相交，故此升起曲线是凹的。第15页，共45页，编辑于2022年，星期六2）选择某排和相邻排oyx-dC(x,0)C2(x+d,0)MM2M1xN1ABN相似于D第16页，共45页，编辑于2022年，星期六再计算相似于第17页，共45页，编
12、辑于2022年，星期六4 模型求解 微分不等式（比较定理）设函数定义在某个区域上，且满足1）在D上满足存在唯一性定理的条件；2）在D上有不等式则初值问题与的解在它们共同存在区间上满足第18页，共45页，编辑于2022年，星期六第19页，共45页，编辑于2022年，星期六所求曲线的近似曲线方程（折衷法）折衷法第20页，共45页，编辑于2022年，星期六5 总结与讨论有时只需求近似解。方法利用微分不等式建模；模型讨论obxyadd1）视点移动时升起曲线如何求得？2）怎样减少地面的坡度？调整参数、相邻排错位。3）衡量经济的指标？座位尽量多、升起曲线占据的空间尽量少等。第21页，共45页，编辑于202
13、2年，星期六 3.4 碳定年代法碳定年代法考古、地质学等方面的专家常用14C测定法（通常称碳定年代法）来估计文物或化石的年代。第22页，共45页，编辑于2022年，星期六v 14C的蜕变规律的蜕变规律v14C是一种由宇宙射线不断轰击大气层，使大气层产生中子，中子与氮气作用生成的具有放射性的物质。这种放射性碳可氧化成二氧化碳，二氧化碳被植物所吸收，而植物又作为动物的食物，于是放射性碳被带到各种动植物体内。v14C是放射性的，无论在空气中还是在生物体内他都在不断蜕变，这种蜕变规律我们可以求出来。通常假定其蜕变速度与该时刻的存余量成正比。第23页，共45页，编辑于2022年，星期六v设在时刻t（年）
14、，生物体中14C的存量为x(t),生物体的死亡时间记为t0=0,此时14C含量为x0,由假设，初值问题 v （1.1）的解为v （1.2）v其中，为常数，k前面的符号表示14C的存量是递减的。（1.2）式表明14C是按指数递减的，而常数k可由半衰期确定，第24页，共45页，编辑于2022年，星期六v若14C的半衰期为T,则有v (1.3)v将（1.3）代入（1.2）得v v即有 （1.4）第25页，共45页，编辑于2022年，星期六v碳定年代法的根据碳定年代法的根据v 活着的生物通过新陈代谢不断摄取14C,因而他们体内的14C与空气中的14C含量相同，而生物死亡之后，停止摄取14C，因而尸体内
15、的14C由于不断蜕变而不断减少。碳定年代法就是根据生物体死亡之后体内14C蜕变减少量的变化情况来判断生物的死亡时间的。第26页，共45页，编辑于2022年，星期六v碳定年代法的计算碳定年代法的计算v由（1.4）解得 v (1.5)v由于x(0),x(t)不便于测量，我们可把(1.5)作如下修改.v对(1.2)式两边求导数，得v (1.6)v而 (1.7)第27页，共45页，编辑于2022年，星期六v(1.6)和(1.7)两式相除，得 将上式代入(1.5)，得 v (1.8)这样由(1.8)可知，只要知道生物体在死亡时体内14C的蜕变速度 和现在时刻t的蜕变速度 ,就可以求得生物体的死亡时间了，
16、在实际计算上，都假定现代生物体中14C的蜕变速度与生物体死亡时代生物体中14C的蜕变速度相同。第28页，共45页，编辑于2022年，星期六v马王堆一号墓年代的确定马王堆一号墓年代的确定v马王堆一号墓于1972年8月出土，其时测得出土的木炭标本的14C平均原子蜕变数为29.78/s,而新砍伐木头烧成的木炭中14C 平均原子蜕变数为38.37/s,又知14C的半衰期为5568年，这样，我们可以把 ,，T=5568 年代入(1.8)，得 v这样就估算出马王堆一号墓大约是在2000多年前。第29页，共45页，编辑于2022年，星期六v两个注记两个注记v（1）马王堆中的古代科技之谜）马王堆中的古代科技之
17、谜v素纱蝉衣素纱蝉衣：两件轻薄的衣服，丝绸，极轻且两千年不腐，南京云锦研究所接受国家科技攻关，用了二十年时间，于1990年成功研制出类似素纱蝉衣的复制品，但该复制品比汉代的还重50克，已不可能再轻了。v女尸千年不腐：病理知识女尸千年不腐：病理知识：女尸解剖显示患有非常严重的冠心病；肺部有肺结核的钙化，肺部钙化是肺结核痊愈后的表现。2000年后的今天，要想控制肺结核，除自身的第30页，共45页，编辑于2022年，星期六v抵抗力要强外，还要有好的营养，要想痊愈是很困难的。两处胆结石，其一在胆总管，有蚕豆大，胆道被堵得水泄不通。三种寄生虫，其中竟有血吸虫，其症状应为腹胀如鼓，骨瘦如柴，但该女子皮下脂
18、肪异常丰满，显然血吸虫被有效的控制住了。该西汉贵妇生前病魔缠身，但从其遗体上未发现长期卧床养病的迹象。一个同时患有这么多疾病的人，能够长期稳定控制病情，在今天也是一个奇迹，说明汉代医术已达到了相当高的水平。第31页，共45页，编辑于2022年，星期六v（2）碳定年代法的不足）碳定年代法的不足 现在，14C年代测定法已受到怀疑，在2500-10000年前这段时间中与其他断代法的结果有差异。1966年，耶鲁实验室的Minze Stuiver 和加利福尼亚大学圣地亚哥分校的Hans E.Suess在一份报告中指出了这一时期使14C年代测定产生误差的根本原因。在那个年代，宇宙射线的放射强度减弱了，偏差
19、的峰值发生在大约6000年以前。第32页，共45页，编辑于2022年，星期六v这两位研究人员的结论出自对Brist/econe松树所作的14C年代测定的结果，因为这种松树同时还提供了精确的年轮断代。他们提出了一个很成功的误差公式，用来校正根据14C断代定出的2300-6000年前这期间的年代：v真正的年代=14C年1.4900。第33页，共45页，编辑于2022年，星期六3.4 范范.梅格伦伪造名画案梅格伦伪造名画案 第二次世界大战比利时解放后，荷兰保安机关开始搜捕纳粹分子的合作者，发现一名三流画家H.A.Vanmeegren曾将17世纪荷兰著名画家Jan.Vermeer的一批名贵油画盗卖给德
20、寇，于1945年5月29日通敌罪逮捕了此人。Vanmeegren被捕后宣称他从未出卖过荷兰的利益，所有的油画都是自己伪造的，为了证实这一切，在狱中开始伪造Vermeer的画耶稣在学者中间。当他的工作快完成时，又获悉他可能以伪造罪被判刑，于是拒绝将画老化，以免留下罪证。第34页，共45页，编辑于2022年，星期六 为了审理这一案件，法庭组织了一个由化学家、物理学家、艺术史学家等参加的国际专门小组，采用了当时最先进的科学方法，动用了X-光线透视等，对颜料成份进行分析，终于在几幅画中发现了现代物质诸如现代颜料钴蓝的痕迹。这样，伪造罪成立，Vanmeegren被判一年徒刑。1947年11月30日他在狱
21、中心脏病发作而死去。但是，许多人还是不相信其余的名画是伪造的，因为，Vanmeegren在狱中作的画实在是质量太差，所找理由都不能使怀疑者满意。直到20年后，1967年，卡内基梅隆大学的科学家们用微分方程模型解决了这一问题。第35页，共45页，编辑于2022年，星期六原理原理著名物理学家卢瑟夫（Rutherford）指出：物质的放射性正比于现存物质的原子数。设 时刻的原子数为 ，则有为物质的衰变常数。初始条件第36页，共45页，编辑于2022年，星期六半衰期碳-14铀-238镭-226铅-210能测出或算出，只要知道 就可算出这正是问题的难处，下面是间接确定 的方法。断代。第37页，共45页，
22、编辑于2022年，星期六油画中的放射性物质油画中的放射性物质 白铅（铅的氧化物）是油画中的颜料之一，应用已有2000余年，白铅中含有少量的铅(Pb210)和更少量的镭(Ra226)。白铅是由铅金属产生的，而铅金属是经过熔炼从铅矿中提取来出的。当白铅从处于放射性平衡状态的矿中提取出来时，Pb210的绝大多数来源被切断，因而要迅速蜕变，直到Pb210与少量的镭再度处于放射平衡，这时Pb210的蜕变正好等于镭蜕变所补足的为止。第38页，共45页，编辑于2022年，星期六铀238镭226铅210钋210铅206（放射性）（无放射性）第39页，共45页，编辑于2022年，星期六假设假设（1）镭的半衰期为
23、1600年，我们只对17 世纪的油画感兴趣，时经300多年，白铅中镭至少还有原量的90%以上，所以每克白铅中每分钟镭的衰变数可视为常数，用 表示。（2）钋的半衰期为138天容易测定，铅210的半衰期为22年，对要鉴别的300多年的颜料来说，每克白铅中每分钟钋的衰变数与铅210的衰变数可视为相等。第40页，共45页，编辑于2022年，星期六建模建模设 时刻每克白铅中含铅210的数量为 ，为制造时刻 每克白铅中含铅210的数量。为铅210的衰变常数。则油画中铅210含量第41页，共45页，编辑于2022年，星期六求解求解均可测出。可算出白铅中铅的衰变率 ，再于当时的矿物比较，以鉴别真伪。矿石中铀的
24、最大含量可能 23%，若白铅中铅210每分钟衰变超过3 万个原子，则矿石中含铀量超过 4%。第42页，共45页，编辑于2022年，星期六测定结果与分析测定结果与分析画名画名钋钋210衰变原子数衰变原子数镭镭226衰变原子数衰变原子数Emmaus的信徒们8.50.82洗足12.60.26读乐谱的妇人10.30.3弹曼陀林的妇人8.20.17做花边的人1.51.4欢笑的女孩5.26.0第43页，共45页，编辑于2022年，星期六若第一幅画是真品，铅210每分钟每克衰变不合理，为赝品。第44页，共45页，编辑于2022年，星期六同理可检验第2，3，4幅画亦为赝品，而后两幅画为真品。第45页，共45页，编辑于2022年，星期六