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