姜启源之代数方程与差分方程模型

,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.02761013 (焦耳)即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. 使 尽量小，如 =0,靠经济实力控制数量不变,结果解释,模型的推广,生产者根据当前时段和前一时段 的价格决定下一时段的产量.,生产者管理水平提高,设供应函数为,需求函数不变,二阶线性常系数差分方程,x0为平衡点,研究平衡点稳定，即k, xkx0的条件,方程通解,(c1, c2由初始条件确定),1, 2特征根，即方程 的根,平衡点稳定，即k, xkx0的条件:,平衡点稳定条件,比原来的条件 放宽了!,模型的推广,6.5 减肥计划节食与运动,背景,多数减肥食品达不到减肥目标，或不能维持.,通过控制饮食和适当的运动，在不伤害身体 的前提下，达到减轻体重并维持下去的目标.,分析,体重变化由体内能量守恒破坏引起.,饮食（吸收热量）引起体重增加.,代谢和运动（消耗热量）引起体重减少.,体重指数BMI=w(kg)/l2(m2). 18.525 超重; BMI30 肥胖.,模型假设,1）体重增加正比于吸收的热量每8000千卡 增加体重1千克；,2）代谢引起的体重减少正比于体重每周每千克 体重消耗200千卡 320千卡(因人而异), 相当于70 千克的人每天消耗2000千卡 3200千卡；,3）运动引起的体重减少正比于体重，且与运动形式 有关；,4）为了安全与健康，每周体重减少不宜超过1.5千克, 每周吸收热量不要小于10000千卡.,某甲体重100千克，目前每周吸收20000千卡热量，体重维持不变。现欲减肥至75千克.,第一阶段：每周减肥1千克，每周吸收热量逐渐减少,直至达到下限（10000千卡）；,第二阶段：每周吸收热量保持下限，减肥达到目标.,2）若要加快进程，第二阶段增加运动，试安排计划.,1）在不运动的情况下安排一个两阶段计划.,减肥计划,3）给出达到目标后维持体重的方案.,确定某甲的代谢消耗系数,即每周每千克体重消耗 20000/100=200千卡,基本模型,w(k) 第k周(末)体重,c(k) 第k周吸收热量, 代谢消耗系数(因人而异),1）不运动情况的两阶段减肥计划,每周吸收20000千卡 w=100千克不变,=1/8000(千克/千卡),第一阶段: w(k)每周减1千克, c(k)减至下限10000千卡,第一阶段10周, 每周减1千克，第10周末体重90千克,吸收热量为,1）不运动情况的两阶段减肥计划,第二阶段：每周c(k)保持Cm, w(k)减至75千克,1）不运动情况的两阶段减肥计划,基本模型,第二阶段：每周c(k)保持Cm, w(k)减至75千克,第二阶段19周, 每周吸收热量保持10000千卡, 体重按 减少至75千克.,运动 t=24 (每周跳舞8小时或自行车10小时), 14周即可.,2）第二阶段增加运动的减肥计划,t每周运动时间(小时),取t=0.003, 即t=24,=1/8000(千克/千卡), =0.025,增加运动相当于提高代谢消耗系数,2）第二阶段增加运动的减肥计划,提高12%,减肥所需时间从19周降至14周,减少25%,这个模型的结果对代谢消耗系数很敏感.,应用该模型时要仔细确定代谢消耗系数 (对不同的人; 对同一人在不同的环境).,3）达到目标体重75千克后维持不变的方案,每周吸收热量c(k)保持某常数C，使体重w不变,不运动,运动(内容同前),6.6 按年龄分组的种群增长,不同年龄组的繁殖率和死亡率不同.,建立差分方程模型，讨论稳定状况下种群的增长规律.,假设与建模,种群按年龄大小等分为n个年龄组，记i=1,2,n,时间离散为时段，长度与年龄组区间相等，记k=1,2,以雌性个体数量为对象.,第i 年龄组1雌性个体在1时段内的繁殖率为bi,第i 年龄组在1时段内的死亡率为di, 存活率为si=1- di,假设 与 建模,xi(k)时段k第i 年龄组的种群数量,按年龄组的分布向量,预测任意时段种群按年龄组的分布,Leslie矩阵(L矩阵),(设至少1个bi0),稳定状态分析的数学知识,L矩阵存在正单特征根1，,若L矩阵存在bi, bi+10, 则,P的第1列是x*,特征向量,解释,L对角化,稳态分析k充分大种群按年龄组的分布, 种群按年龄组的分布趋向稳定，x*称稳定分布, 与初始分布无关。, 各年龄组种群数量按同一倍数增减， 称固有增长率,3）=1时, 各年龄组种群数量不变, 1个个体在整个存活期 内的繁殖数量为1,稳态分析,存活率 si是同一时段的 xi+1与 xi之比,（与si 的定义 比较）,3）=1时,人口模型,连续型人口模型的离散形式,xi(k)k年i 岁的女性人数(模型只考虑女性人口).,bi(k)k年i 岁女性生育率(每人平均生育女儿数).,dii 岁女性死亡率，si=1-di存活率,i1, i2生育区间,k年育龄女性平均生育女儿数,总合生育率(生育胎次),年龄分布向量,hi生育模式,人口模型,存活率矩阵,生育模式矩阵,x(k)状态变量, (k)控制变量,双线性方程(对x(k), (k)线性),原模型,