第六章--代数方程与差分方程模型

6.1投入产出模型6.2CT技术的图像重建6.3原子弹爆炸的能量估计6.4市场经济中的蛛网模型6.5减肥计划——节食与运动6.6按年龄分组的种群增长第六章代数方程与差分方程模型•国民经济各个部门之间存在着相互依存和制约关系，每个部门将其他部门的产品或半成品经过加工（投入）变为自己的产品（产出）.•根据各部门间投入和产出的平衡关系，确定各部门的产出水平以满足社会的需求.•20世纪30年代由美国经济学家列昂节夫提出和研究.•从静态扩展到动态，与数量经济分析方法日益融合，应用领域不断扩大.6.1投入产出模型背景建立静态投入产出数学模型，讨论具体应用.投入产出表国民经济各部门间生产和消耗、投入和产出的数量关系产出投入农业工业建筑业运输邮电批零餐饮其他服务外部需求总产出农业464788229131271312842918工业499860514444035571223408316814建筑业593202312426912875运输邮电62527128163671464771570批零餐饮79749140431302739272341其他服务146128527222521954227255414初始投入1663485165970312183093总投入2918168142875157023415414中国2002年投入产出表（产值单位：亿元）直接消耗系数表产出投入农业工业建筑业运输邮电批零餐饮其他服务农业0.1590.0470.0800.0080.0540.002工业0.1710.5120.5020.2570.2380.226建筑业0.0020.0010.0010.0130.0100.023运输邮电0.0210.0310.0450.1040.0290.027批零餐饮0.0270.0450.0490.0270.0560.050其他服务0.0500.0760.0950.1430.0940.100一个部门的单位产出对各个部门的直接消耗中国2002年直接消耗系数表由投入产出表直接得到农业每1亿元产出直接消耗0.159亿元农业产品直接消耗0.171亿元工业产品反映国民经济各个部门之间的投入产出关系投入产出的数学模型xi~第i部门的总产出di~对第i部门的外部需求xij~第i部门对第j部门的投入aij~直接消耗系数——第j部门单位产出对第i部门的直接消耗xij~第j部门总产出对第i部门的直接消耗每个部门的总产出等于总投入xj~第j部门的总投入injijidxx1jijijxxa/ijnjijidxax1设共有n个部门nnijaA)(Tnxxx),(1Tnddd),(1dAxx技术水平没有明显提高模型应用假设直接消耗系数不变问题1如果某年对农业、工业、建筑业、运输邮电、批零餐饮和其他服务的外部需求分别为1500,4200,3000,500,950,3000亿元,问这6个部门的总产出分别应为多少？d=(1500,4200,3000,500,950,3000)TA由直接消耗系数表给出6个部门的总产出x=(3277,17872,3210,1672,2478,5888)（亿元）.dAxx(),IAxddAIx1)(求解模型应用dAxx(),IAxddAIx1)(dAIx1)(总产出对外部需求线性Δd~d增加1个单位x的增量若农业的外部需求增加1单位Td)0,0,0,0,0,1(1)(AIΔx为的第1列6个部门的总产出分别增加1.2266，0.5624，0.0075，0.0549，0.0709，0.1325单位.问题2如果6个部门的外部需求分别增加1个单位,问它们的总产出应分别增加多少？求解其余外部需求增加1单位1)(AIΔx为的其余各列6.2CT技术的图像重建•CT(计算机断层成像)技术是20世纪50至70年代由美国科学家科马克和英国科学家豪斯费尔德发明的.•1971年第一代供临床应用的CT设备问世.•螺旋式CT机等新型设备被医疗机构普遍采用.•CT技术在工业无损探测、资源勘探、生态监测等领域也得到了广泛的应用.背景什么是CT，它与传统的X射线成像有什么区别？光源人眼光源人眼一个半透明物体嵌入5个不同透明度的球概念图示单方向观察无法确定球的数目和透明度让物体旋转从多角度观察能分辨出5个球及各自的透明度人体内脏胶片传统的X射线成像原理CT技术原理探测器X射线X光管人体内脏CT技术:在不同深度的断面上,从各个角度用探测器接收旋转的X光管发出、穿过人体而使强度衰减的射线;经过测量和计算将人体器官和组织的影像重新构建.图像重建X射线强度衰减与图像重建的数学原理•射线强度的衰减率与强度成正比.I~射线强度l~物质在射线方向的厚度μ~物质对射线的衰减系数IdldII0~入射强度•射线沿直线L穿行,穿过由不同衰减系数的物质组成的非均匀物体(人体器官).leII0LIIdlyx0ln),(LdlyxII)),(exp(0I0L0yxμ(x,y)Ldlyxl)),(X射线强度衰减与图像重建的数学原理LIIdlyx0ln),(右端数值可从CT的测量数据得到多条直线L的线积分Ldlyx),(被积函数μ(x,y)LfdlyxfLP),()(FQ(q)~与Q相距q的直线L的线积分Pf(L)对所有q的平均值拉东变换0)(1)(qqdFQfQ拉东逆变换图像重建反映人体器官大小、形状、密度的图像数学原理实际上只能在有限条直线上得到投影(线积分).图像重建在数学方法上的进展，为CT技术在各个领域成功的和不断拓广的应用提供了必要条件.图像重建的代数模型Δlj每个像素对射线的衰减系数是常数m个像素(j=1,…,m),n束射线(i=1,…,n)iII)/ln(0~Li的强度测量数据μj~像素j的衰减系数Δlj~射线在像素j中的穿行长度J(Li)~射线Li穿过的像素j的集合LIIdlyx0ln),(像素j射线LiniIIliLJjjji,,2,1,)/ln(0)(σjLiσlij图像重建的代数模型常用算法设像素的边长和射线的宽度均为σ,)/ln(0iibIIjjx中心线法aij~射线Li的中心线在像素j内的长度lij与σ之比.iLJjjjIIli)/ln(0)(bAx面积法aij~射线Li的中心线在像素j内的面积sij与σ之比.sijnibxamjijij,,2,1,1中心法aij=1~射线Li经过像素j的中心点.图像重建的代数模型中心法的简化形式假定射线的宽度为零,间距σaij=1~Li经过像素j内任一点987654321L4L8L7L6L5L3L2L1σ100100000010110100001011010000000001011000000110011001100110011000000100AbAx根据A和b,由确定像素的衰减系数向量xm和n很大且mn,方程有无穷多解+测量误差和噪声在x和e满足的最优准则下估计x代数重建技术(ART)beAx6.3原子弹爆炸的能量估计1945年7月16日美国科学家在新墨西哥州阿拉莫戈多沙漠试爆了全球第一颗原子弹,震惊世界!当时资料是保密的,无法准确估计爆炸的威力.英国物理学家泰勒研究了两年后美国公开的录像带,利用数学模型估计这次爆炸释放的能量为19.2千吨.后来公布爆炸实际释放的能量21千吨t(ms)r(m)t(ms)r(m)t(ms)r(m)t(ms)r(m)t(ms)r(m)0.1011.10.8034.21.5044.43.5361.115.0106.50.2419.90.9436.31.6546.03.8062.925.0130.00.3825.41.0838.91.7946.94.0764.334.0145.00.5228.81.2241.01.9348.74.3465.653.0175.00.6631.91.3642.83.2659.04.6167.362.0185.0泰勒测量：时刻t所对应的“蘑菇云”的半径r原子弹爆炸的能量估计爆炸产生的冲击波以爆炸点为中心呈球面向四周传播,爆炸的能量越大，在一定时刻冲击波传播得越远.冲击波由爆炸形成的“蘑菇云”反映出来.泰勒用量纲分析方法建立数学模型,辅以小型试验,又利用测量数据对爆炸的能量进行估计.物理量的量纲长度l的量纲记L=[l]质量m的量纲记M=[m]时间t的量纲记T=[t]动力学中基本量纲L,M,T速度v的量纲[v]=LT-1导出量纲221rmmkf加速度a的量纲[a]=LT-2力f的量纲[f]=LMT-2引力常数k的量纲[k]对无量纲量，[]=1(=L0M0T0)量纲齐次原则=[f][l]2[m]-2=L3M-1T-2在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系.量纲齐次原则等式两端的量纲一致量纲分析~利用量纲齐次原则寻求物理量之间的关系.例：单摆运动)1(321glmt321][][][][glmtlmgm求摆动周期t的表达式设物理量t,m,l,g之间有关系式1,2,3为待定系数，为无量纲量2/12/10321glt(1)的量纲表达式glt2与对比33212TLMT12003321对x,y,z的两组量测值x1,y1,z1和x2,y2,z2,p1=f(x1,y1,z1),p2=f(x2,y2,z2)2121pppp为什么假设这种形式?设p=f(x,y,z)),,(),,(),,(),,(222111222111czbyaxfczbyaxfzyxfzyxfx,y,z的量纲单位缩小a,b,c倍zyxzyxf),,(p=f(x,y,z)的形式为),,(),,,(22221111czbyaxfpczbyaxfp量纲齐次原则321glmt单摆运动0002010010101004321)()()()(TMLTMLTMLTMLTMLyyyy000241243TMLTMLyyyyy201001010100][][][][TMLgTMLlTMLmTMLt单摆运动中t,m,l,g的一般表达式0),,,(glmtf020041243yyyyyglt12)/(glt4321yyyyglmty1~y4为待定常数,为无量纲量0)(FT)1,1,0,2(Tyyyyy),,,(4321基本解设f(q1,q2,,qm)=0mjXqniaijij,,2,1,][1ys=(ys1,ys2,…,ysm)T,s=1,2,…,m-rF(1,2,…,m-r)=0与f(q1,q2,,qm)=0等价,F未定.Pi定理(Buckingham)是与量纲单位无关的物理定律，X1,X2,Xn是基本量纲,nm,q1,q2,qm的量纲可表为,}{mnijaA量纲矩阵记作rrankA若线性齐次方程组0Ay有m-r个基本解，记作mjyjssjq1为m-r个相互独立的无量纲量,且则记爆炸能量为E，将“蘑菇云”近似看成一个球形.时刻t球的半径为rt,E空气密度ρ,大气压强P基本量纲：L,M,T21322][;][;][;][;][MTLPMLMTLETtLr),,,(PEtr原子弹爆炸能量估计的量纲分析方法建模r与哪些因素有关？rtEρP20210111001320153ALMT量纲矩阵0),,,,(PEtrfy=(1,-2/5,-1/5,1/5,0)y=(0,6/5,-2/5,-3/5,1)T原子弹爆炸能量估计的量纲分析方法建模rARank3ARank5/125/15/15/21EtrErt5/132565/35/25/6