数学模型 第06章(第五版)ppt课件.pptx

差分方程若干离散点上未知变量数值的方程.,描述离散时间段上客观对象的动态变化过程.,现实世界中随时间连续变化的动态过程的近似.,差分方程与代数方程都是离散模型的数学表述，二者有着类似的向量-矩阵表达形式，求解过程也存在相互联系.,第六章 差分方程与代数方程模型,6.1 贷款购房6.2 管住嘴迈开腿6.3 市场经济中的物价波动6.4 动物的繁殖与收获6.5 信息传播6.6 原子弹爆炸的能量估计6.7 CT技术的图像重建6.8 等级结构6.9 中国人口增长预测,每月还多少钱,贷款购房需考虑的问题,买多大的房子,一共贷多少钱,网上的房贷计算器,轻击鼠标即得,输入必要信息,贷款购房最简单的差分方程模型,6.1 贷款购房,单利和复利,单利 万元存5年定期, 年利率4.75%, 到期后本息(本金加利息)：10000(1+0.04755)=12375元.,两种计算利息的基本方式,复利 万元存1年定期, 年利率为3%, 到期不取则自动转存, 5年后本息：10000 (1+0.03)5=11593元.,利滚利！,复利本息：(1+r)n,单位本金、同一利率r、同一存期n计算单利和复利：,单利本息：1+nr,1+nr,零存整取 每月固定存额，约定存款期限，到期一次支取本息的定期储蓄.,按单利计算的业务零存整取,单利和复利,方式：5元起存，多存不限，存期1年、3年、5年.,勤俭节约、科学理财,例 每月存入3000元，存期年（年利率3.5%）,零存整取 计算器,累计存入金额180,000元到期本息总额196,012.50元,单利和复利,按单利计算的业务零存整取,a每月存入金额, r 月利率, n 存期(月),xk 存入k个月后的本息,k=n递推至k=,a =3000, r =0.035/12, n =125 (月),x1=a+ar,xk= xk-1+a+akr, k2,3, n,x2= x1+a+a2r,xn= na+ar(1+2+n),= na+ar (+) ,等额本息贷款和等额本金贷款,房贷计算器的选项,贷款类别：商业贷款, 公积金, 组合型,年利率不同,计算方法：根据贷款总额或面积、单价计算.,按揭年数：可选至30年. 选择20年.,银行利率：基准利率、利率上限或下限. 选择商业 贷款的基准利率6.55%.,还款方式：等额本息还款或等额本金还款.,等额本息贷款和等额本金贷款,例1“房贷计算器”选择等额本息还款, 输入: 商业贷款总额100万元, 期限20年, 年利率6.55%.,建立等额本息还款方式的数学模型, 并作数值计算,等额本息还款每月归还本息(本金加利息)数额相同.,等额本金还款每月归还本金数额相同, 加上所欠本金的利息.,点击“开始计算”得: 还款总额1796447.27元, 月均还款7485.2元.,等额本息还款模型,a每月还款金额,x0 贷款总额,r 月利率,n 贷款期限(月),xk 第k月还款后尚欠金额,xk= xk-1(1+r)a, k1,2, n,k=n递推至k=,贷款到期时xn=0,xn= x0(1+r)na1+(1+r)+(1+r)n-1,= x0(1+r)na (1+r)n (1+r)n1,本月欠额=上月欠额的本息还款金额,等额本息还款模型,A1 还款总额,a每月还款金额,x0 贷款总额,r 月利率,n 贷款期限(月),例1 x0 =100(万元), r=0.0655/12, n=1220=240(月),与房贷计算器给出的相同,例“房贷计算器”选择等额本金还款, 输入: 商业贷款总额100万元, 期限20年, 年利率6.55%.,等额本息贷款和等额本金贷款,建立等额本金还款方式的数学模型, 并作数值计算.,点击“开始计算”得到: 还款总额1657729.17元, 每月还款金额由第月的9625元逐月递减, 最后月为4189.41元.,等额本金还款模型,x0 贷款总额,r 月利率,n 贷款期限(月),xk 第k月还款金额,还款金额逐月减少归还本金x0/n所产生的利息x0r/n,每月归还本金x0/n,第1月还款金额,k=n递推至k=2,等额本金还款模型,x0 贷款总额,r 月利率,n 贷款期限(月),xk 第k月还款金额,A2 还款总额,例2 x0 =100(万元), r=0.0655/12, n=1220=240(月),与房贷计算器给出的相同,等额本息与等额本金方式的比较,等额本息方式简单，便于安排收支.,等额本金方式每月还款金额前期高于等额本息方式, 后期低于等额本息方式, 适合当前收入较高人群.,等额本息方式还款总额大于等额本金方式.,等额本息方式前期还款额较少, 所欠本息的利息逐月归还, 所以利息总额较大.,例1 例2： A1=1796447.27(元) , A2=1657729.17(元).,模型适用于任何还款周期(半月、一季度等)将公布的年利率折换为一个还款周期的利率.,小结与评注,贷款购房两种基本还款方式：等额本息、等额本金.,要点: 明确利息计算, 列出差分方程, 利用递推关系.,不同还款周期一次还款金额和还款总额都不一样.,周期越短还款总额越小？,测评体重的标准体重指数(BMI Body Mass Index),BMI=w/l2，w体重(kg), l身高(m).,例.l=1.70m, w=63.5 kg,多数减肥食品达不到减肥效果，或不能维持.,通过控制饮食和适当运动，在不伤害身体的 前提下，达到减轻体重并得以控制的目的.,标准的身材!,6.2 管住嘴迈开腿,模型分析,人体通过食物摄入热量, 通过代谢和运动消耗热量.,二者平衡,体重不变.,分析对热量的吸收和消耗, 建立体重变化规律的模型.,平衡被破坏则体重变化.,减肥计划应以不伤害身体为前提.,增加运动量是加速减肥的有效手段.,以周为时间单位制订减肥计划.,1）体重增加正比于吸收的热量, 平均8000kcal 增加体重1kg.,2）代谢引起的体重减少正比于体重, 每周每千克 体重消耗200 320kcal (因人而异).,3）运动引起的体重减少正比于体重, 且与运动 形式和运动时间有关.,模型假设,4）为了安全与健康, 每周吸收热量10000kcal,且每周 减少量1000kcal; 每周体重减少量 1.5kg.,调查资料,食物每百克所含热量,运动每小时每千克体重消耗热量,模型假设,基本模型,c(k) 第k周吸收热量 (kcal),w(k)第k周(初)体重 (kg) , k=1,2,热量转换系数,平均8000kcal增加体重1kg,代谢系数(因人而异).,由和吸收热量 c(k)决定体重w(k)的变化规律.,减肥计划的提出,某人身高1.70m, 体重100kg, BMI高达34.6. 目前每周吸收20000kcal热量,体重长期未变.,制订减肥计划使体重减至75kg(BMI=26)并维持下去.,1. 在正常代谢情况下安排一个两阶段计划: 第一阶段：吸收热量每周减少1000kcal, 直至达到 安全下限10000 kcal/周；,2. 为加快进程而增加运动，重新安排两阶段计划.,3. 给出达到目标后维持体重不变的方案.,第二阶段：每周吸收热量保持下限, 达到减肥目标 .,减肥计划的制定,1. 确定某人的代谢消耗系数,每周每千克体重消耗 20000/100=200kcal,每周吸收20000kcal, 体重100kg不变.,正常代谢消耗相当弱.,=1/8000,c(k)= c, w(k+1)=w(k)= w,w=100,c=20000, + = +c(k) ,2. 正常代谢情况下的第一阶段计划,吸收热量由20000kcal每周减少1000kcal, 直至达到安全下限10000 kcal/周., + = +c(k) ,c(10)= 10000,第一阶段需10周,3. 正常代谢情况下的第二阶段计划,吸收热量保持下限cmin =10000kcal/周,体重减至75kg., + = +c(k) ,两阶段计划共需32周.,第二阶段需22周.,4. 为加快进程而增加运动,t每周运动时间 (h),运动每小时每千克体重消耗热量,取 t=40,w(11)= 89.3319 kg，w(11+12)= 74.7388kg,第二阶段缩短为12周,两阶段计划共需22周.,(如每周步行7h加乒乓4h),5. 检验“每周体重减少量 1.5kg”, + = . +., =, . +., =, + = . +., =, . +., =,正常代谢,增加运动,编程计算w(k),w(k) w(k+1) 1.5kg,6. 达到目标后维持体重不变的方案,每周吸收热量保持常数c使体重w=75kg不变.,=15000kcal/周,c由20000kcal/周直接减至15000,14000,13000,12000时体重w(k)下降曲线.,c =14000时w(72)=75kg,c =12000时w(40)=75kg,比两阶段计划的时间长，吸收热量突减对身体不利.,75kg,7. 达到目标体重所需时间与每周吸收热量的关系,令目标体重w*=w(n+1), 记初始体重w1= w(1), + = +c(k) ,k=1递推至k=n,w*=75w1=100,小结与评注,减肥科学化、定量化需要研究人体体重变化的规律.,计算中由于增加运动使由0.025提高到0.03时(变化20%), 减肥所需时间从32周减少到22周(变化约30%)体重变化对相当敏感.,体重变化既有普遍规律也与每个人特殊生理条件有关，特别是代谢消耗系数.,消费者在自由竞争的市场经济中常会遇到商品价格的波动现象.,供大于求,商品数量与价格在波动,6.3 市场经济中物价的波动,物价的波动,商品数量和价格主要由供求关系决定.,供求平衡,供求失衡,波动的两种形式,振幅逐渐减小，最终趋向平稳.,振幅越来越大，如不干预将导致经济崩溃.,讨论政府的干预方式,描述波动现象,研究趋向平稳的条件,建立数量价格模型,模型假设,xk第k时段商品数量,yk 第k时段商品价格,时段生产周期 (饲养周期、种植周期),2. yk由消费者需求关系决定,3. xk+1由生产者供应关系决定,4. xk, yk偏离x0, y0不大时, 偏离yk-y0与xk-x0成正比,yky0, 价格过低,xkx0, 供过于求,偏离xk+1-x0与yk-y0成正比.,差分方程模型,xk, yk的差分方程组,消去yk-y0,xk的差分方程模型,x0, y0稳定,xkx0,yky0,xk, yk,x0, y0不稳定,差分方程模型,模型分析,例. 平衡状态：x0=100, y0=10(元).,设x1=110,x0, y0稳定,xkx0,yky0,x0, y0不稳定,例. 平衡状态：x0=100, y0=10(元).,设x1=110,=0.1,=5,模型分析,xk,yk, 商品数量减少1单位, 价格上涨幅度, 价格上涨1单位, (下时段)供应的增量, 消费者对需求的敏感程度, 生产者对价格的敏感程度,小, 有利于经济稳定, 小, 有利于经济稳定,模型分析,消费者需求关系,生产者供应关系,f与g的交点P0(x0,y0) 平衡点,一旦xk=x0，则yk=y0,且 xk+1=xk+2=x0 , yk+1=yk+2= =y0,蛛网模型,xk第k时段商品数量,yk 第k时段商品价格, 0 =( 0 ), +1 0 =( 0 ),差分方程模型的图形表示,x1,设x1偏离x0,P0是稳定平衡点,P0是不稳定平衡点,蛛网模型,斜率取绝对值,差分方程模型与蛛网模型的一致,x0, y0稳定,xkx0,yky0,xk, yk,x0, y0不稳定, 0 =( 0 ), +1 0 =( 0 ),差分方程模型,蛛网模型,P0是稳定平衡点,P0是不稳定平衡点,直线 f 斜率 Kf =,直线 g 斜率 Kg=1/, 1/, 1/,政府的干预办法,1. 使 尽量小，如 =0,以行政手段控制价格不变,2. 使 尽量小，如 =0,靠经济实力控制数量不变,差分方程模型的推广, +1 0 =( 0 ),根据当前和前一时段的价格决定下一时段的产量.,生产者管理水平和素质提高,xk，yk的差分方程组,已知, 及x0, y0, 由初始值x1, x2递推地计算xk, yk .,差分方程模型的推广,x0, y0不稳定,原模型,新模型,=0.24, =5 不变,xk,yk, +1 0 =( 0 ),xkx0,yky0, 0 =( 0 ), =1.21,x0, y0稳定, =1.21,讨论稳定条件,二阶线性常系数差分方程,差分方程模型的推广, 0 =( 0 ),1, 2特征根, 特征方程, 稳定平衡点,1, 21,比原模型的稳定条件 放宽了.,差分方程模型的推广,特征方程,稳定条件(xkx0): 1, 21,生产者管理水平和素质的提高有利于经济稳定！,政府干预措施具有人们熟知的现实背景.,小结与评注,对市场经济中“供不应求价格上涨、供过于求价格下跌”的现象用两种模型描述和解读：,差分方程：便于运算,蛛网模型：直观鲜明,模型参数有明确的经济学含义敏感系数.,差分方程平衡点的稳定性有明显的实际意义，反映了数学与现实的密切关系.,野生动物种群在自然环境下繁殖、成长、死亡,不同年龄动物的数量比例保持平衡.,饲养动物种群在人类控制下,使不同年龄动物的数量比例达到稳定的预期目标.,建立动物种群的自然增长模型.,讨论饲养动物种群的稳定收获.,6.4 动物的繁殖与收获,按年龄分组的动物种群增长模型,不同年龄动物的繁殖率、死亡率差别较大.,建立按年龄分组种群增长的差分方程模型.,时段与年龄组相对应.,种群通过雌性繁殖而增长.,总体数量按性别比计算.,讨论稳定状况下种群的增长规律.,种群按照年龄等间隔地分为n个年龄组.,时间分成与年龄组区间大小相等的时段.,每个年龄组的雌性个体在一个时段内的繁殖率 和死亡率不随时段变化.,模型假设,在稳定环境下和不太长时期内,模型建立,xi(k)第i年龄组第k时段的种群数量, i=1,2, n，k=0,1,2,bi第i年龄组的繁殖率 (每个雌性个体一个时 段繁殖的数量).,di 第i年龄组的死亡率(一个时段内死亡数量 在总量中的比例).,si=1- di 存活率,bi0，至少一个bi 0.,0si1, sn=0,模型建立,第1年龄组(出生婴儿)k+1时段数量=各年龄组k时段繁殖数量之和.,k时段第i年龄组存活的部分到k+1时段演变为第i+1年龄组.,n个变量的差分方程组,已知bi, si 及xi(0),按年龄分组的种群增长模型,k=0,1,2,模型建立,Leslie矩阵(矩阵L), +1 +1 = ,按年龄分组的种群数量,x(k)的归一化向量，按年龄分组的分布向量.,Leslie模型,模型求解,例. 种群分5个年龄组, 繁殖率为 b1=0, b2=0.2, b3=1.8, b4=0.8, b5=0.2, 存活率为s1=0.5, s2=0.8, s3=0.8, s4=0.1, 各年龄组初始数量均为100只.,求任意时段各年龄组数量x(k)及分布向量x*(k).,x(0),x(k)= 1 , 2 , 5 T,模型求解,x(k)= 1 , 2 , 5 T,x*(k)趋向稳定,x(k)仍在增长,k充分大,结果分析,分析k充分大后x(k), x*(k)的变化规律,稳定状态分析的数学知识,对应特征向量 x,矩阵L存在最大特征根 (正单根),c常数,满足,结果分析,k充分大x(k), x*(k)的特性,特征向量,1. 分布向量 稳定分布.,2. 数量,与初始分布无关.,各年龄组数量按同一倍数 (固有增长率)增减.,3. =1时 x(k) cx ,si等于同一时段相邻年龄组的数量比.,结果分析,用算例验证x(k), x*(k)的特性,x*=0.4559, 0.2223, 0.1734, 0.1353, 0.0132T,模型求解中x*(30)近似于x*,1. 由L计算得到=1.0254,2. 模型求解中xi(30)与xi(29)之比约为=1.0254.,3. =1.0254比1略大, xi+1(30)与 xi(30)之比近似于si,饲养动物种群的持续稳定收获模型,同一年龄组种群的收获量在每个时段都相等,实现方法: 每个年龄组每个时段种群的增长量=同一时段的收获量.,控制饲养动物各年龄组的数量, 实现持续稳定收获:,假定自然环境下饲养动物仍服从种群增长模型:,种群数量始终不变.,xi(k)第i年龄组第k时段的种群数量.,hi 第i年龄组种群的收获系数(收获量与总量之比),模型建立,增长量 = 收获量,种群增长模型,实现持续稳定收获种群数量x(k)=x (对k不变),模型建立,持续稳定收获,的最大特征根,给定bi, si，选择收获系数hi,持续稳定收获, 种群数量的稳定分布：,模型建立,的特征向量( ),持续稳定收获,增长量 = 收获量,模型求解,例. 设一个种群分成个年龄组, 各年龄组的繁殖率 为 b1=0, b2=5, b3=2, 存活率为s1=0.8, s2=0.5.,确定各年龄组的收获系数以实现持续稳定收获.,持续稳定收获的条件,1. 取h1=0, h2=0.75, h3=1,求种群及收获量按年龄组的稳定分布.,持续稳定收获,2. 取h1=0.5, h2=0.5, h3=1,模型求解,满足持续稳定收获条件,2. 出售50%的幼畜和成年牲畜及全部老年牲畜.,收获量的稳定分布,种群数量的稳定分布,1. 不出售幼畜, 出售75%成年牲畜及全部老年牲畜.,1. h1=0, h2=0.75, h3=1,2. h1=0.5, h2=0.5, h3=1,1. ，2.,1. ，2.,小结与评注,人口增长与动物种群数量变化规律相同, 类似建立离散型女性人口模型Leslie模型,模型基本假定：种群参数(繁殖率、存活率)只与年龄有关, 与时段无关(稳定环境、时间不长).,Leslie矩阵为常数矩阵L可用特征根方法作稳定性分析.,6.5 信息传播,当今每天都有大量的、正面或负面的信息，甚至谣言，通过各种传统的、近代的、特别是互联网的渠道，在几乎没有限制的人群中传播.,考察总人数一定的封闭环境, 开始极少数人得到了一条信息或制造了一条谣言, 然后通过人与人之间的交流在人群中传播, 使获知的人越来越多.,在合理的简化假设下，建立数学模型来描述信息传播的规律，研究其发展趋势.,模型假设,1. 在封闭环境中人群的总人数不变.,2. 信息通过已获知的人向未获知的人传播.,已获知信息的人数越多 (传播人群), 每天新获知信息的人数越多.,未获知信息的人数越多 (潜在人群), 每天新获知信息的人数越多.,3. 每天新获知信息的人数与已获知信息的人数和未获知信息的人数的乘积成正比.,模型建立,k第k天已获知信息人数 (传播人群)，N pk第k天未获知信息人数(潜在人群).,假设：第k+1天新获知信息人数pk+1 pk=pk与pk和N pk的乘积成正比.,N总人数,对于潜在人群的一位而言,每天新获知信息人数的百分比增量.,c是反映传播速度的参数，c越大传播速度越快,模型求解,pk+1 pk = cpk(N pk),设N=1000，初始获知信息的人数 p0=10.,80天后pk才接近N,k第k天获知信息人数,k接近N只需40天,logistic微分方程的离散形式差分方程,pk+1=(1+N) ( + ),logistic微分方程的离散形式差分方程,无法得到pk的显式表达式,关注k时pk的变化,=+, = + ,一阶非线性差分方程,通过讨论的平衡点及其稳定性研究pk的收敛性质(k).,预备知识6-1 差分方程的类型、求解及稳定性,标准形式k+1=(1)的平衡点及其稳定性,解代数方程 x=f(x)=bx(1 x) 得到非零平衡点 ,为判断x*的稳定性计算f(x)=bx(1 x)在x*的导数,根据x*稳定性条件 ,=+,只需讨论b3时xk的变化规律(k).,不同b值下k+1=(1)的计算结果(初值x0=0.2),xkx* (单调收敛),b=1.7,b=2.9,b=3.15,b=3.525,b=3.565,b=3.7,xkx* (振荡收敛),xk的2个子序列趋向另外2个平衡点,xk的4个子序列趋向另外4个平衡点,xk的8个子序列趋向另外8个平衡点,xk不趋向任何平衡点, 出现混沌现象,拓展阅读6-1 差分方程k+1=(1)的收敛、分岔和混沌,模型讨论,pk+1=(1+N) ( + ),k+1=(1),平衡点 =/,平衡点 p*=N,1b3, 稳定,cN2, p*稳定,=+, = + ,b2, xk单调收敛于x*,cN1, pk单调收敛于p*,出现cN1, pk振荡收敛 、分岔、混沌,23, 出现分岔、混沌,？,标准形式,原方程,模型讨论,pk+1=(1+N) ( + ),设总人数N固定，讨论参数c (传播速度) 的上限.,传播过程中对于任意的k都有pk+1 N.,cN1 ?,c= (pk+1 )/ ,k可以无限接近N.,模型拓展,信息自由传播,a 每天被制止传播谣言的人数比例.,pk+1 pk = cpk(N pk),被制止传播谣言的人又加入到听信谣言的潜在人群中.,=+-a, = +a ,被制止传播谣言的人既不再传播也不会听信，从此退出这个信息传播系统.,增加新的人群被制止传播后退出系统者.,模型拓展,qk+1 qk = apk,qk第k天退出系统者的人数,pk+1 pk = cpk(N pk qk ) apk,a 每天被制止传播谣言的人数比例.,k, qk 联立组成非线性差分方程组.,给定 p0和q0可递推计算pk和qk,小结与评注,讨论方程k+1=(1)的平衡点及稳定性,说明非常简单的非线性差分方程, 也会出现比相应的logistic微分方程有着远为复杂的特性.,信息传播模型的解只能具有单调增的收敛形式，完全符合人们的直观认识.,实际上信息传播速度c不可能保持不变，用随k而变的ck代替c, 仍可递推计算pk, 结果会更好.,受限环境下传染病的蔓延和生物种群的增长都可以建立类似的数学模型.,1945年7月16日美国科学家在新墨西哥州的阿拉莫戈多沙漠试爆了全球第一颗原子弹, 震惊世界!,当时资料是保密的, 无法准确估计爆炸的威力.,英国物理学家泰勒研究了两年后美国公开的录像带, 利用数学模型估计这次爆炸释放的能量为19.2103t.,后来公布爆炸实际释放的能量为21103t,6.6 原子弹爆炸的能量估计,泰勒测量： 时刻t 所对应的“蘑菇云”的半径r,原子弹爆炸的能量估计,爆炸产生的冲击波以爆炸点为中心呈球面向四周传播, 爆炸的能量越大, 在一定时刻冲击波传播得越远.,冲击波由爆炸形成的“蘑菇云”反映出来.,泰勒用量纲分析方法建立数学模型, 辅以小型试验,又利用测量数据对爆炸的能量进行估计.,物理量的量纲,长度 l 的量纲记 L=l,质量 m的量纲记 M=m,时间 t 的量纲记 T=t,动力学中基本量纲 L, M, T,速度 v 的量纲 v=LT-1,导出量纲,加速度 a 的量纲 a=LT-2,力 f 的量纲 f=LMT-2,引力常数 k 的量纲 k,对无量纲量，=1(=L0M0T0),量纲齐次原则,=fl2m-2=L3M-1T-2,在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系.,量纲齐次原则,等式两端的量纲一致,量纲分析利用量纲齐次原则寻求物理量之间的关系.,例：单摆运动,求摆动周期 t 的表达式,设物理量 t, m, l, g 之间有关系式,1, 2, 3 为待定系数，为无量纲量,(1)的量纲表达式,与 对比,对 x,y,z的两组测量值x1,y1,z1 和x2,y2,z2, p1 = f( x1,y1,z1), p2 = f( x2, y2,z2 ),为什么假设这种形式?,设p= f(x,y,z),x,y,z的量纲单位缩小a,b,c倍,量纲齐次原则,单摆运动,单摆运动中 t, m, l, g 的一般表达式,基本解,设 f(q1, q2, , qm) = 0,ys = (ys1, ys2, ,ysm)T , s = 1,2, m-r,F( 1, 2, m-r ) = 0 与 f (q1, q2, , qm) =0 等价, F未定.,定理 (Buckingham),是与量纲单位无关的物理定律，X1,X2, ,Xn 是基本量纲, nm, q1, q2, ,qm 的量纲可表为,量纲矩阵记作,记爆炸能量为E，将“蘑菇云”近似看成一个球形.,时刻 t 球的半径为 r,t, E,空气密度, 大气压强P,基本量纲：L, M, T,原子弹爆炸能量估计的量纲分析方法建模,r与哪些因素有关？,量纲矩阵,y=(1,-2/5,-1/5,1/5,0) y=(0,6/5,-2/5,-3/5,1)T,原子弹爆炸能量估计的量纲分析方法建模,原子弹爆炸能量估计的数值计算,时间 t 非常短能量 E 非常大,泰勒根据一些小型爆炸试验的数据建议,用r, t 的实际数据做平均,空气密度 =1.25 (kg/m3),1103t (TNT能量)= 4.1841012J,实际值21 103t,泰勒的计算,最小二乘法拟合 r=atb,E=8.02761013 (J), 即19.2 103t,取y平均值得c=6.9038,模型检验,b=0.4058,2/5,量纲分析法的评注,物理量的选取,基本量纲的选取,基本解的构造,结果的局限性, () = 0中包括哪些物理量是至关重要的.,基本量纲个数n; 选哪些基本量纲.,有目的地构造 Ay=0 的基本解.,方法的普适性,函数F和无量纲量未定.,不需要特定的专业知识.,CT(计算机断层成像 )技术是20世纪50至70年代由美国科学家科马克和英国科学家豪斯费尔德发明的.,1971年第一代供临床应用的CT设备问世.,螺旋式CT机等新型设备被医疗机构普遍采用.,CT技术在工业无损探测、资源勘探、生态监测等领域也得到了广泛的应用.,背景,什么是CT，它与传统的X射线成像有什么区别？,6.7 CT技术的图像重建,一个半透明物体嵌入5个不同透明度的球,概念图示,单方向观察无法确定球的数目和透明度,让物体旋转从多角度观察能分辨出5个球及各自的透明度,人体内脏,胶片,传统的X射线成像原理,CT技术原理,探测器,X射线,X光管,人体内脏,CT技术: 在不同深度的断面上,从各个角度用探测器接收旋转的X光管发出、穿过人体而使强度衰减的射线;,经过测量和计算将人体器官和组织的影像重新构建.,图像重建,X射线强度衰减与图像重建的数学原理,射线强度的衰减率与强度成正比.,I射线强度,l物质在射线方向的厚度,物质对射线的衰减系数,I0入射强度,射线沿直线L穿行, 穿过由不同衰减系数的物质组成的非均匀物体(人体器官).,X射线强度衰减与图像重建的数学原理,右端数值可从CT 的测量数据得到,多条直线L的线积分,FQ(q)与Q相距q的直线L的线积分Pf(L)对所有q的平均值,Radon变换,Radon逆变换,图像重建,数学原理,实际上只能在有限条直线上得到投影(线积分).,图像重建在数学方法上的进展，为CT技术在各个领域成功的和不断拓广的应用提供了必要条件.,图像重建的代数模型,每个像素对射线的衰减系数是常数,m个像素(j=1, m),n束射线(i=1,n),Li的强度测量数据,j像素j的衰减系数,lj射线在像素j中的穿行长度,J(Li)射线Li穿过的像素j的集合,图像重建的代数模型,常用算法,设像素的边长和射线的宽度均为,中心线法,aij射线Li的中心线在像素j内的长度lij与之比.,面积法,aij射线Li的中心线在像素j内的面积sij与之比.,中心法,aij=1射线Li经过像素j的中心点.,图像重建的代数模型,中心法的简化形式,假定射线的宽度为零, 间距,aij=1 Li经过像素j内任一点,根据A和b, 由 确定像素的衰减系数向量x,m和n很大且m n, 方程有无穷多解,+ 测量误差和噪声,在x和e满足的最优准则下估计x,代数重建技术(ART),社会系统中需要适当且稳定的等级结构.,描述等级结构的演变过程，预测未来的结构.,确定为达到某个理想结构应采取的策略.,引起等级结构变化的因素：,系统内部等级间的转移：提升和降级.,系统内外的交流：调入和退出(退休、调离等).,用马氏链模型描述确定性的转移问题 (将转移比例视为概率).,6.8 等级结构,基本模型,a(t)等级结构,等级 i=1,2,k（如助教、讲师、教授）,数量分布 n(t)=(n1(t), n2(t), ,nk(t)ni(t) t 年属于等级i 的人数， t =0,1,比例分布 a(t)=(a1(t), a2(t), ,ak(t),转移矩阵 Q=pijkk, pij 是每年从i 转至j 的比例,基本模型,ri 每年调入 i 的比例 (在总调入人数中),ij 每年从i 转至j 的比例,基本模型, 基本模型,基本模型,等级结构a(t) 状态概率,P转移概率矩阵,用调入比例进行稳定控制,问题：给定Q, 哪些等级结构可以用合适的调入比例保持不变,用调入比例进行稳定控制,a*,例 大学教师(助教、讲师、教授)等级 i=1,2,3，已知每年转移比例,可行域A,稳定域B,用调入比例进行稳定控制,研究稳定域B的结构,用调入比例进行稳定控制,稳定域B是k维空间中以 si 为顶点的凸多面体,研究稳定域B的结构,例,S1,稳定域B是以si为顶点的三角形,等级结构的演变、预测和控制在社会系统中有 广泛应用.,讨论总人数和内部转移比例不变情况下, 用调入 比例控制级结构的变化.,建立等级结构演变过程的基本方程, 预测未来结构.,各种推广情况：总人数按照一定比例增长； 调入比例有界；调入比例固定而用内部转移比例 控制级结构的变化.,小结与评注,6.9 中国人口增长预测,中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。,2007年公布的国家人口发展战略研究报告提出，如果我国人口总量峰值2033年前后达到峰值15亿人左右，全国总和生育率应保持在1.8左右，过高或过低都不利于人口与经济社会的协调发展.,要求从中国实际情况和人口增长特点出发，参考2005年人口抽样数据，建立中国人口增长的数学模型，并对人口增长中短期和长期趋势做出预测。,以发表在工程数学学报2007年增刊二上学生优秀论文为基本材料, 介绍建模过程.,近年来中国的人口发展出现了一些新的特点, 如老龄化进程加速、出生人口性别比持续升高、乡村人口城镇化等因素, 影响着中国人口的增长.,中国人口增长预测全国大学生数学建模竞赛2007年A 题,2005年人口抽样数据,20岁城市男性占城市总人口的0.69%, 死亡率为20岁城市男性总人口的0.59,20岁城市生育子女的妇女占20岁城市妇女总人口的29.01 ,问题分析,用数学建模预测人口增长的方法：差分方程、微分方程、回归分析、时间序列等.,结合所给数据以差分方程组的Leslie模型为基础.,考虑不同地区、不同性别人口参数的差别及农村人口向城市迁移等因素.,按照地区和性别建立以时间和年龄为基本变量的中国人口增长模型.,利用历史数据估计生育率、死亡率及人口迁移等参数，代入模型求解并作预测.,模型假设,中国人口是封闭系统, 将数据中的市、镇合并为城市, 与农村(乡)作为两个地区; 只考虑农村向城市人口的单向迁移, 不考虑与境外的相互移民.,对中短期人口预测, 生育率、死亡率及人口迁移等参数用历史数据估计; 长期预测考虑总和生育率的控制、城镇化指数的变化趋势等因素.,女性每胎生育一个子女.,模型建立,k年, i岁(不满i+1岁), 地区j(j=1城市, j=2农村), 性别l (l=m男性, l=w女性)的人口数量.,死亡率(死亡人数占总人数的比例).,存活率.,生育率(k年每位i岁女性平均生育婴儿数).,婴儿性别比(男婴比例).,i1,i2育龄区间.,迁移数量(迁入为正, 迁出为负).,模型建立, +1,+1 = , , + , , i=1,2,n-1,l=m,w,j=1,2, 1,+1 = 0, = 1 2 , , + 0, ,=1,2,模型建立,k年每位育龄女性的生育数.,假定所有女性在育龄期间都保持这个生育数.,表述、控制人口增长的重要指标.,生育率 (k年每位i岁女性生育数)分析,生育模式(i岁女性生育数在育龄女性 中的比例).,模型建立,农村向城市迁移的人口数量,假设单向迁移,人口分布向量,存活率矩阵,生育模式矩阵,人口迁移向量,模型建立,按地区和性别划分、以年龄为离散变量、随时段演变的人口发展模型.,因生育过程使女性人口出现在男性人口方程中, 导致方程之间的耦合作用.,4n阶差分方程组,参数估计, +1 = ( , , , , , ), =,=1,2,=0,1,2,人口发展模型表达式,需要先估计参数：存活率 , 婴儿性别比 ,总和生育率 ,生育模式 , 人口迁移 .,再由 递推计算 +1 .,存活率的估计,死亡率与年龄关系大, 与地区、性别和时间的关系小.,中短期预测：将过去若干年不同地区、性别和各年龄人口的死亡率简单地取平均值.,长期预测：用统计方法对历史数据加以处理，并参考发达国家人口死亡率的演变过程给出估计值.,参数估计,中国几十年来死亡率降低较快, 未来趋势仍持续下降.,参数估计,生育率的估计,中国女性生育率已在持续下降后大致保持稳定.,中短期预测：将过去若干年不同地区、性别和各年龄人口的生育率简单地取平均值.,长期预测：设定几个不同水平的总和生育率,生育模式可采用概率论中的分布:,ic生育高峰年龄,婴儿性别比可由数据中0岁男婴和女婴的比率得出.,人口迁移的估计,参数估计,农村迁移的人口总数按照城镇化率的增长估算., k年全国城镇、农村和总人口, 城镇化率,根据历史数据预测,k年农村迁移人口总数,=(k)-(k-1)X(k),迁移人口年龄、性别分布用典型城市的资料代替.,迁移人口比例,模型求解, +1 = ( , , , , , ), =,=1,2,=0,1,2,选定初始年份用人口发展模型递推计算,全面、完整地描述人口的演变过程.,人口指数简明反映一个国家或地区的人口特征.,人口总数, = =1 2 =0 , + , ,平均年龄,模型求解,人口指数,平均寿命,按k年死亡率计算的k年出生人口平均存活时间.,老龄化指数,平均年龄与平均寿命之比,抚养指数(k),每个劳动力平均抚养的无劳动力人数,L(k)劳动力人数,im1, im2, iw1, iw2 劳动力年龄,小结与评注,人口迁移是比较复杂的因素, 历史数据不多、不完整, 模型中人口迁移部分比较粗糙.,生育在一定意义上可以控制, 应该对总和生育率和生育模式设定几种方案,对预测结果进行分析.,差分方程组形式的Leslie模型简明地表述了人口发展过程中生育和死亡两个自然因素的作用.,死亡基本上不可控制, 稳定环境下死亡率变化慢, 用历史数据能较好地对中短期预测做出估计.,