版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、为铅210的衰变常数一、不同领域中的微分方程模型1 社会及市场经济中的微分方程模型1综合国力的微分方程模型2诱发投资与加速开展的微分方程模型3经济调整的微分方程模型4广告的微分方程模型5价格的微分方程模型 2 战争中的微分方程模型1军备竞赛的微分方程模型2战争的微分方程模型3战斗中生存可能性的微分方程模型4战争的预测与评估模型3 人口与动物世界的微分方程模型1单种群模型及进展开发的单种群模型2弱肉强食模型3两个物种在同一生态龛中的竞争排斥模型4无管理的鱼类捕捞模型5人口预测与控制模型4 疾病的传染与诊断的微分方程模型 1艾滋病流行的微分方程模型2糖尿病诊断的微分方程模型3人体内碘的微分方程模型
2、4药物在体内的分布与排除模型5 自然科学中的微分方程模型1人造卫星运动的微分方程模型2航空航天器翻滚控制的微分方程模型3非线性振动的微分方程模型4PLC电路自激振荡的微分方程模型 5盯梢与追击问题的微分方程模型 1 根据规律列方程二、建立微分方程模型的一般方法 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。2 微元分析法 利用的定理与规律寻找微元之间的关系式,与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。3 模拟近似法 在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有所理解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建
3、立能近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性质,再去同实际情况比照,检验此模型能否刻画、模拟某些实际现象。三、微分方程的解法求解微分方程有三种方法:1求准确解;2求数值解近似解;3定性理论方法。四、典型案例分析问题:案例一 范. 梅格伦Van Meegren伪造名画案 第二次世界大战比利时解放后,荷兰保安机关开场搜捕纳粹分子的合作者,发现一名三流画家曾将17世纪荷兰著名画家的一批名贵油画盗卖给德寇,于1945年5月29日以通敌罪逮捕了此人。Vanmeegren被捕后声称他从未出卖过荷兰的利益,所有的油画都是自己伪造的,为了证实这一切,在狱中开场伪造Vermeer的画?耶
4、稣在学者中间?。当他的工作快完成时,又得悉他可能以伪造罪被判刑,于是回绝将画老化,以免留下罪证。 为了审理这一案件,法庭组织了一个由化学家、物理学家、艺术史学家等参加的国际专门小组,采用了当时最先进的科学方法,动用了X-光线透视等,对颜料成份进展分析,终于在几幅画中发现了现代物质诸如现代颜料钴蓝的痕迹。这样,伪造罪成立, Vanmeegren被判一年徒刑。1947年11月30日他在狱中心脏病发作而死去。 但是,许多人还是不相信其余的名画是伪造的,因为, Vanmeegren在狱中作的画实在是质量太差,所找理由都不能使疑心者满意。直到20年后,1967年,卡内基梅隆大学的科学家们用微分方程模型解
5、决了这一问题。 著名物理学家卢瑟夫Rutherford指出:物质的放射性正比于现存物质的原子数。问题的分析:设 t 时刻的原子数为N(t),那么有其中 为物质的衰变常数,初始条件为 这是一阶常微分方程初值问题,其解为 从中解得 半衰期公式为 , 而碳-14、镭-226、铀-238、铅-210等的半衰期T分别为5568年、1600年、45亿年、22年,N(t)可以测定出,只要知道 就可算出断代。这正是问题的难处,下面是间接确定 的方法。 油画中的放射性物质 白铅铅的氧化物是油画中的颜料之一,应用已有2000余年,白铅中含有少量的铅(Pb210)和更少量的镭(Ra226)。白铅是由铅金属产生的,而
6、铅金属是经过熔炼从铅矿中提取来出的。当白铅从处于放射性平衡状态的矿中提取出来时, Pb210的绝大多数来源被切断,因此要迅速蜕变,直到Pb210与少量的镭再度处于放射平衡,这时Pb210的蜕变正好等于镭蜕变所补足的为止。以以下图是铀裂变示意图 (放射性)(无放射性)铀238镭226铅210钋210铅2061镭的半衰期为1600年,我们只对17 世纪的油画感兴趣,时经300多年,白铅中镭至少还有原量的90%以上,所以每克白铅中每分钟镭的衰变数可视为常数,用r表示。建模: 2钋的半衰期为138天容易测定,铅210的半衰期为22年,对要鉴别的300多年的颜料来说,每克白铅中每分钟钋的衰变数与铅210
7、的衰变数可视为相等。模型假设: 设t时刻每克白铅中含铅210的数量为y(t), 为制造时刻 每克白铅中含铅210的数量。 为铅210的衰变常数。则油画中铅210含量为模型求解: 这是一个一阶线性微分方程的初值问题,由常数变易法得其解为 变形为均可测出,可算出白铅中铅的衰变率 ,再与当时的矿物比较,以鉴别真伪。测定分析与结论:测定结果如以下图 画名铅210衰变原子数镭226衰变原子数Emmaus的信徒们8.50.82洗足12.60.26读乐谱的妇人10.30.3弹曼陀林的妇人8.20.17做花边的人1.51.4欢笑的女孩5.26.0对第一幅画,各量为 可见铅210每分钟每克衰变不合理,故为赝品。
8、同理可检验第2,3,4幅画亦为赝品,而后两幅画为真品。注:类似于上例,在考古、地质学等方面专家常用碳定年代法也是依靠放射性物质的性质来估计文物或化石的年代。代入 计算得而矿石中铀的最大含量为 23%,假设白铅中铅210每分钟衰变超过3 万个原子,那么矿石中含铀量超过 4%。问题:案例二放射性废料的处理问题美国原子能委员会以往处理浓缩的放射性废料时,一直采用把它们装入密封的圆桶里扔到水深约为91米海底的方法。对此,科学家们表示担忧,怕圆桶下沉到海底时与海底碰撞发生破裂而造成核污染。原子能委员会分辩说不会发生这种情况。为此,工程师们进展了碰撞试验,发现当圆桶下沉速度超过米/秒与海底碰撞时,圆桶就可
9、能发生破裂。这样,为防止圆桶碰裂,需要计算一下圆桶下沉到海底时速度是多少。圆桶重量为千克,体积为立方米,海水浮力为千克/立方米。于是,假设圆桶下沉速度小于米/秒,说明原处理放射性废料的方法是平安可靠的,否那么,应该禁用原方法处理放射性废料。大量试验说明圆桶下沉时的阻力与圆桶的方位大致无关,而与下沉的速度成正比,比例系数为。你能判断美国原子能委员会以往处理浓缩的放射性废料方法是否合理吗?符号说明:w: 圆桶重量,这里为千克V: 圆桶体积,这里为立方米B:海水浮力,这里为千克k: 圆桶下沉时的阻力系数,这里为v: 圆桶下沉时的速度D: 圆桶下沉时的阻力,这里为kvt: 圆桶分开海平面下沉的时间,单
10、位为秒y (t):圆桶在t时刻下沉的深度,单位为米问题分析与求解: 因为圆桶下沉满足牛顿第二定律:运动物体受到的力等于该物体的质量与其运动加速度的乘积,即F=ma。在本问题中有,F=W-B-D=W-B-kv,y(0)=0,v(0)=0,于是可以得到如下微分方程:(1)重力浮力阻力求解问题,Matlab命令为用Matlab软件求解。求解出函数y(t)后,求出圆桶下落到水深为91米海底的时间,即求满足y(t)=91的时间t1,然后求出在时间t1的速度即可以得到问题的解答。这里的速度关系由直接求函数y(t)对t的导数得到。syms t;syms y;y=dsolve(m*D2y+k*Dy-w+b=0
11、,y(0)=0,Dy(0)=0)运行结果为y =-m/k2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k2即运行结果如图 b=1025.94*0.208;k=0.12;w=239.456;m=w/9.8; % 这里重力加速度取为米/秒2 t=0:0.1:20;y=-m/k2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k2-91;plot(t,y)求解方程y(t)=91,先画出图形,命令为可见方程的根在t=13附近, 为求解方程,建立M文件:function y=my0(t) b=1025.94*0.208; k=0.12;
12、 w=239.456;m=w/9.8; y=-m/k2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k2-91;结果为t=fsolve(my0, 13) 用如下的命令求解:求y(x)的导数,Matlab代码为b=1025.94*0.208;k=0.12;w=239.456;m=w/9.8;syms t;y1=diff(-m/k2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k2-91)运行结果为 y1 = 21/4276*t)+81439/375即计算结果说明圆桶下沉到水深约为91米海底时的速度约为米/秒,这个速度大于米/秒的
13、速度,因此,计算结果说明美国原子能委员会以往处理浓缩的放射性废料方法是不平安的。故t时的速度为v问题:案例三 弱肉强食微分方程模型生活在同一环境中的各类生物之间,进展着残酷的生存竞争。设想一海岛,居住着狐狸与野兔,狐吃兔,兔吃草,青草如此之丰富,兔子们无无食之忧,于是大量繁殖;兔子一多,狐易得食,狐量亦增,而由于狐狸数量增加吃掉大量兔子,狐群又进入饥饿状态而使其总数下降,这时兔子相对平安,于是兔子总数上升。就这样,狐兔数目交替地增减,无休止的循环,遂形成生态的动态平衡。如何用建立数学模型描绘并预测下一阶段情况?变量说明: 模型各项意义其中a, b, c, d 均为正常数。 x(t)t时刻兔子数
14、目y(t)t时刻狐狸数目数学模型: ax表示兔子的繁殖速度与现存兔子数成正比;-bxy表示狐兔相遇,兔子被吃掉的速度;-cy表示狐狸因同类争食造成的死亡速度与狐狸总数成正 比;dxy表示狐兔相遇,对狐狸有好处而使狐狸繁殖增加的速 度。 这是意大利著名物理学著名生物数学家Volterra建立的微分方程模型,这是达尔文主义的数学表达,称为Volterra-Lotka模型。 设该问题的初始条件为 令 模型求解与分析:得系统的两个稳定点(奇点)为 和取x(t)=0,求解问题()与(),得其解为 取y(t)=0,求解问题()与(),得其解为 说明x轴和y轴都是方程组的轨线,它们的实际意义是:假设没有兔子
15、,狐狸因为没有食物将趋于死亡;而假设没有狐狸,兔子将无限制增长。 当x , y 0 时,方程组等价于这是一个可别离变量的方程,积分得通解两边取指数得 令 证明:当初值 时,方程()确定了xoy面上位于第一象限的一族封闭曲线。,由于 易证f(y)在 处取得极大值 ; 同理g(x)在 处取得极大值 ; 所以当 时,方程无 x 0 , y 0 的解; 当 时,方程有唯一解; 当 时,令F(x , y)=g(x) f (y),则F(x , y)在第一象限内连续,从而对第一象限内任意一条连接原点O与点S的连续曲线L,由于根据连续函数的介值定理,在曲线L上必有一点 ,使得 ,即曲线F(x , y)= C与L相交。由L的任意性知,当 时,曲线()是封闭的。由于这些封闭曲线都不经过奇点S,可见任何从x(0)0, y(0) 0出发的轨线都具有周期性,即函数x(t) , y(t)满足而 这说明假设初始时既有兔子又有狐狸,那么随着时间的增加,狐兔数量呈周期性变化,无休止的形成动态的生态平衡。 此外,当 时,由()式知 ,即y随着时间递增,从而(x(t),y(t)沿曲线逆向运动。 对于模型(),可用定积分计算的平均值。记平均值为 ,则由于 ,故同理可得 结果说明,当开场时,既有兔子又有狐狸,那么随着时间的推移,狐兔数目
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 全国各地区教育信息化发展现状及对策研究考试及答案
- 直播电商社交行为研究论文
- 基于建构主义的单元教学设计研究-以“空间向量与立体几何”单元为例
- 蒙医辩证治疗琪苏希拉偏盛型过敏性紫癜性肾炎的疗效研究及白介素-6的相关性研究
- 数字化赋能下敏捷组织文化建设机制研究
- 快消品行业数字化渠道变革研究
- 新质生产力背景下就业结构变化研究及对策
- 高中生成长型语言思维、二语坚毅对自我调节学习策略使用的影响研究
- 高中生英语词汇学习策略使用与词汇深度知识水平的相关性研究
- 幼儿园中班智育教育知识清单
- 驾校安全职责考核制度
- GB/T 47005-2026增材制造激光定向能量沉积钛合金制件技术规范
- 无人机复杂气象飞行作业规范手册
- 青年红色筑梦之旅课件
- GB/T 3033-2025船舶与海上技术管路系统内含物的识别颜色
- 皮带输送机巡检培训
- 邮政企业会计核算操作流程
- 物业服务酬金制合同范本及风险提示
- 防汛抗旱基本知识培训课件
- 液氨制冷车间拆除方案(3篇)
- 学校交通安全知识培训课件
评论
0/150
提交评论