数学模型(第四版)幻灯片-第六章

,6.1 投入产出模型 6.2 CT技术的图像重建 6.3 原子弹爆炸的能量估计与量纲分析 6.4 市场经济中的蛛网模型 6.5 减肥计划节食与运动 6.6 按年龄分组的种群增长,第六章 代数方程与差分方程模型,国民经济各个部门之间存在着相互依存和制约关系，每个部门将其他部门的产品或半成品经过加工（投入）变为自己的产品（产出）.,根据各部门间投入和产出的平衡关系，确定各部门的产出水平以满足社会的需求 .,20世纪30年代由美国经济学家列昂节夫提出和研究.,从静态扩展到动态，与数量经济分析方法日益融合，应用领域不断扩大 .,6.1 投入产出模型,背景,建立静态投入产出数学模型，讨论具体应用.,投入产出表,国民经济各部门间生产和消耗、投入和产出的数量关系,中国2002年投入产出表（产值单位：亿元）,直接消耗系数表,一个部门的单位产出对各个部门的直接消耗,中国2002年直接消耗系数表,由投入产出表直接得到,农业每1亿元产出直接消耗0.159亿元农业产品,直接消耗0.171亿元工业产品,反映国民经济各个部门之间的投入产出关系,投入产出的数学模型,xi第i部门的总产出,di对第i部门的外部需求,xij第i部门对第j部门的投入,aij直接消耗系数第j部门单位产出对第i部门的直接消耗,xij第j部门总产出对第i部门的直接消耗,每个部门的总产出等于总投入,xj第j部门的总投入,设共有n个部门,技术水平没有明显提高,模型应用,问题1 如果某年对农业、工业、建筑业、运输邮电、批零餐饮和其他服务的外部需求分别为1500, 4200, 3000, 500, 950, 3000亿元, 问这6个部门的总产出分别应为多少？,d=(1500, 4200, 3000, 500, 950, 3000)T,A由直接消耗系数表给出,6个部门的总产出 x=(3277, 17872, 3210, 1672, 2478, 5888)（亿元）.,求解,模型应用,总产出对外部需求线性,dd增加1个单位,x的增量,若农业的外部需求增加1单位,x为 的第1列,6个部门的总产出分别增加1.2266，0.5624，0.0075，0.0549，0.0709，0.1325单位.,问题2 如果6个部门的外部需求分别增加1个单位, 问它们的总产出应分别增加多少？,求解,其余外部需求增加1单位,x为 的其余各列,6.2 CT技术的图像重建,CT(计算机断层成像 )技术是20世纪50至70年代由美国科学家科马克和英国科学家豪斯费尔德发明的.,1971年第一代供临床应用的CT设备问世.,螺旋式CT机等新型设备被医疗机构普遍采用.,CT技术在工业无损探测、资源勘探、生态监测等领域也得到了广泛的应用.,背景,什么是CT，它与传统的X射线成像有什么区别？,一个半透明物体嵌入5个不同透明度的球,概念图示,单方向观察无法确定球的数目和透明度,让物体旋转从多角度观察能分辨出5个球及各自的透明度,人体内脏,胶片,传统的X射线成像原理,CT技术原理,探测器,X射线,X光管,人体内脏,CT技术: 在不同深度的断面上,从各个角度用探测器接收旋转的X光管发出、穿过人体而使强度衰减的射线;,经过测量和计算将人体器官和组织的影像重新构建.,图像重建,X射线强度衰减与图像重建的数学原理,射线强度的衰减率与强度成正比.,I射线强度,l物质在射线方向的厚度,物质对射线的衰减系数,I0入射强度,射线沿直线L穿行, 穿过由不同衰减系数的物质组成的非均匀物体(人体器官).,X射线强度衰减与图像重建的数学原理,右端数值可从CT 的测量数据得到,多条直线L的线积分,FQ(q)与Q相距q的直线L的线积分Pf(L)对所有q的平均值,拉东变换,拉东逆变换,图像重建,数学原理,实际上只能在有限条直线上得到投影(线积分).,图像重建在数学方法上的进展，为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.3 原子弹爆炸的能量估计,1945年7月16日美国科学家在新墨西哥州阿拉莫戈多沙漠试爆了全球第一颗原子弹, 震惊世界!,当时资料是保密的, 无法准确估计爆炸的威力.,英国物理学家泰勒研究了两年后美国公开的录像带, 利用数学模型估计这次爆炸释放的能量为19.2千吨.,后来公布爆炸实际释放的能量21千吨,泰勒测量： 时刻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未定.,Pi定理 (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),1千吨(TNT能量) = 4.184*1012焦尔,实际值21千吨,泰勒的计算,最小二乘法拟合 r=atb,E=8.0276×1013 (焦耳)即19.2千吨,取y平均值得c=6.9038,模型检验,b=0.4058,2/5,量纲分析法的评注,物理量的选取,基本量纲的选取,基本解的构造,结果的局限性, () = 0中包括哪些物理量是至关重要的.,基本量纲个数n; 选哪些基本量纲.,有目的地构造 Ay=0 的基本解.,方法的普适性,函数F和无量纲量未定.,不需要特定的专业知识.,物理模拟示例：波浪对航船的阻力,航船阻力 f,航船速度v, 船体尺寸l, 浸没面积 s, 海水密度, 重力加速度g .,量纲分析在物理模拟中的应用,物理模拟: 按照一定的比例尺寸构造它的物理模型,通过对模型的研究得出原型的结果.,量纲分析可以指导物理模拟中比例尺寸的确定.,物理模拟示例：波浪对航船的阻力,定理,原型船,模型船,模型船的 均已知,当原型船的 给定后计算 f,物理模拟,物理模拟示例：波浪对航船的阻力,原型船,模型船,模拟条件,量测模型船阻力f'，可计算 f.,无量纲化示例：火箭发射,星球表面竖直发射火箭。初速v, 星球半径r, 星球表面重力加速度g.,研究火箭高度 x 随时间 t 的变化规律.,t=0 时 x=0, 火箭质量m1, 星球质量m2,牛顿第二定律，万有引力定律,3个独立参数,用无量纲化方法减少独立参数个数,用参数r,v,g的组合,分别构造与x,t具有相同量纲的xc, tc （特征尺度）,无量纲变量,如,令,xc, tc的不同构造,1）令,为无量纲量,用无量纲化方法减少独立参数个数,3）令,2）令,用无量纲化方法减少独立参数个数,1) 2) 3) 的共同点,1) 2) 3) 的重要差别,考察无量纲量,在1) 2) 3) 中能否忽略以为因子的项？,1),无解,无量纲化方法,2),3),1) 2) 3) 的重要差别,无量纲化方法,原问题,是原问题的近似解,1) 2) 3) 的重要差别,无量纲化方法,为什么3)能忽略项，得到原问题近似解，而1) 2)不能?,3）令,火箭到达最高点时间为v/g, 高度为v2/2g,大体上具有单位尺度,无量纲化方法,选择特征尺度的一般讨论见：林家翘著自然科学中确定性问题的应用数学,无 量 纲 化,无量纲化是研究物理问题常用的数学方法.,选择特征尺度主要依赖于物理知识和经验.,恰当地选择特征尺度可以减少独立参数个数，还可以辅助确定舍弃哪些次要因素.,6.4 市场经济中的蛛网模型,问 题,供大于求,现 象,商品数量与价格的振荡在什么条件下趋向稳定?,当不稳定时政府能采取什么干预手段使之稳定?,描述商品数量与价格的变化规律.,商品数量与价格在振荡,蛛 网 模 型,xk第k时段商品数量；yk第k时段商品价格.,消费者的需求关系,生产者的供应关系,减函数,增函数,f与g的交点P0(x0,y0) 平衡点,一旦xk=x0，则yk=y0,xk+1,xk+2,=x0, yk+1,yk+2, =y0,设x1偏离x0,x1,P0是稳定平衡点,P0是不稳定平衡点,曲线斜率,蛛 网 模 型,在P0点附近用直线近似曲线,P0稳定,P0不稳定,方 程 模 型,方程模型与蛛网模型的一致, 商品数量减少1单位, 价格上涨幅度, 价格上涨1单位, (下时段)供应的增量,考察 , 的含义, 消费者对需求的敏感程度, 生产者对价格的敏感程度,小, 有利于经济稳定, 小, 有利于经济稳定,结果解释,xk第k时段商品数量；yk第k时段商品价格.,结果解释,经济不稳定时政府的干预办法,1. 使 尽量小，如 =0,以行政手段控制价格不变,2