三章微分方程和差分方程模型研究报告.pptVIP

第三章 微分方程和差分方程模型; 在研究实际问题时, 我们常常不能直接得出变量之间的关系,但却能容易得出包含变量导数在内的关系式,这就是微分方程. 在现实社会中,又有许多变量是离散变化的,如人口数、生产周期与商品价格等, 而且离散的运算具有可操作性, 差分正是联系连续与离散变量的一座桥梁. 不管是微分方程还是差分方程模型，有时无法得到其解析解(必要时,可以利用计算机求其数值解),既使得到其解析解,尚有未知参数需要估计(这是可利用第二章参数估计方法). 而在实际问题中,讨论问题的解的变化趋势很重要，因此，以下只对其平衡点的稳定性加以讨论.;3.1 微分方程模型 ;稳定性判别方法;如果;3.2 差分方程模型 ; 若有常数a是差分方程(3-6)的解, 即; 二阶常系数线性差分方程 xn+2 + axn+1 + bxn = r, 其中a, b, r为常数.; ① 当?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*+ ? n (C1cosn? + C2sinn? ). 易知,当且仅当特征方程的任一特征根 |?i |＜1时, 平衡点x*是稳定的. ;对于一阶非线性差分方程 xn+1 = f (xn );3.3 观众厅地面设计;建立坐标系;2 问题的假设;3 建模;2）选择某排;再计算;4 模型求解 ;;所求曲线的近似曲线方程（折衷法）;5 总结与讨论; 3.4 碳定年代法; 14C的蜕变规律 14C是一种由宇宙射线不断轰击大气层，使大气层产生中子，中子与氮气作用生成的具有放射性的物质。这种放射性碳可氧化成二氧化碳，二氧化碳被植物所吸收，而植物又作为动物的食物，于是放射性碳被带到各种动植物体内。 14C是放射性的，无论在空气中还是在生物体内他都在不断蜕变，这种蜕变规律我们可以求出来。通常假定其蜕变速度与该时刻的存余量成正比。;设在时刻t（年），生物体中14C的存量为x(t),生物体的死亡时间记为t0=0,此时14C含量为x0,由假设，初值问题 （1.1）的解为 （1.2） 其中，为常数，k前面的符号表示14C的存量是递减的。（1.2）式表明14C是按指数递减的，而常数k可由半衰期确定，;若14C的半衰期为T,则有 (1.3) 将（1.3）代入（1.2）得 即有 （1.4） ; 碳定年代法的根据 活着的生物通过新陈代谢不断摄取14C,因而他们体内的14C与空气中的14C含量相同，而生物死亡之后，停止摄取14C，因而尸体内的14C由于不断蜕变而不断减少。碳定年代法就是根据生物体死亡之后体内14C蜕变减少量的变化情况来判断生物的死亡时间的。;碳定年代法的计算 由（1.4）解得 (1.5) 由于x(0),x(t)不便于测量，我们可把(1.5)作如下修改. 对(1.2)式两边求导数，得 (1.6) 而 (1.7); (1.6)和(1.7)两式相除，得 将上式代入(1.5)，得 (1.8) 这样由(1.8)可知，只要知道生物体在死亡时体内14C的蜕变速度 和现在时刻t的蜕变速度 ,就可以求得生物体的死亡时间了，在实际计算上，都假定现代生物体中14C的蜕变速度与生物体死亡时代生物体中14C的蜕变速度相同。 ;马王堆一号墓年代的确定 马王堆一号墓于1972年8月出土，其时测得出土的木炭标本的14C平均原子蜕变数为29.78/s,而新砍伐木头烧成的木炭中14C 平均原子蜕变数为38.37/s,又知14C的半衰期为5568年，这样，我们可以把 , ， T=5568 年代入(1.8)，得 这样就估算出马王堆一号墓