6.1 投入产出模型投入产出模型6.2 CT技术的图像重建技术的图像重建6.3 原子弹爆炸的能量估计原子弹爆炸的能量估计6.4 市场经济中的蛛网模型市场经济中的蛛网模型6.5 减肥计划减肥计划节食与运动节食与运动6.6 按年龄分组的种群增长按年龄分组的种群增长第六章第六章 代数方程与差分方程模型代数方程与差分方程模型梁俗桨超灶克粟阳腰腹挖庚滦喷铁俏诧踩蛀坯校碉吾扁缚冬蛇憋测达瓤傈投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 国民经济各个部门之间存在着国民经济各个部门之间存在着相互依存和制约相互依存和制约关系，关系，每个部门将其他部门的产品或半成品经过加工（每个部门将其他部门的产品或半成品经过加工（投入投入）变为自己的产品（变为自己的产品（产出产出）. 根据各部门间根据各部门间投入和产出的平衡投入和产出的平衡关系，确定各部关系，确定各部门的产出水平以满足社会的需求门的产出水平以满足社会的需求 . 20世纪世纪30年代由美国经济学家列昂节夫提出和研究年代由美国经济学家列昂节夫提出和研究. 从静态扩展到动态，与数量经济分析方法日益融合，从静态扩展到动态，与数量经济分析方法日益融合，应用领域不断扩大应用领域不断扩大 .6.1 投入产出模型投入产出模型背景背景建立静态投入产出数学模型，讨论具体应用建立静态投入产出数学模型，讨论具体应用. 丙妇舍香臆缔嗅肮言哇雅得粕魔汐瘩棋奉扳咏黄均蛮宪扩粮衔踪旭未躲澡投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出表投入产出表 国民经济各部门间生产和消耗、投入和产出的数量关系国民经济各部门间生产和消耗、投入和产出的数量关系 产出出投入投入农业工工业建筑建筑业运运输邮电批零批零餐餐饮其他其他服服务外部外部需求需求总产出出农业464788229131271312842918工工业499860514444035571223408316814建筑建筑业593202312426912875运运输邮电62527128163671464771570批零餐批零餐饮79749140431302739272341其他服其他服务146128527222521954227255414初始投入初始投入1663485165970312183093总投入投入2918168142875157023415414中国中国2002年投入产出表（产值单位：亿元）年投入产出表（产值单位：亿元） 募畅谋巩噬绞贝盂甚殊黍钵驮鞍这地泌辣砚炊鸭锹对融汁坪驳侈涪钧订跟投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量直接消耗系数表直接消耗系数表 产出出投入投入农业工工业建筑建筑业运运输邮电批零餐批零餐饮其他服其他服务农业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亿元工业产品亿元工业产品 反映国民经济各个部门之间的投入产出关系反映国民经济各个部门之间的投入产出关系绘秩香矽稳夷轰员张讽凋殴拈诸蒋茬齿踩浅迅沮蜜句啃幢标优腥搞印矽熔投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出的数学模型投入产出的数学模型 xi第第i部门的总产出部门的总产出di对第对第i部门的部门的外部需求外部需求xij第第i部门对第部门对第j部门的投入部门的投入aij直接消耗系数直接消耗系数第第j部部门单位产出单位产出对对第第i部门的直接消耗部门的直接消耗 xij第第j部部门总产出出对第第i部门的直接消耗部门的直接消耗每个部门的总产出等于总投入每个部门的总产出等于总投入 xj第第j部门的总投入部门的总投入 设共有共有n个个部门部门体自莉挺或墓牛蹄粪败坤粉趾憋莱票隔连松喳怕察凸田溉葫莫苦怕犬仰狮投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量技术水平没有明显提高技术水平没有明显提高模型应用模型应用 假设直接消耗系数不变假设直接消耗系数不变 问题问题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)（亿元）（亿元）.求解求解精卉饰擎迄乱况禹刻贪椿意休阁活那绑汗伶衍羌根占奔些燃股猫棺孝轰铃投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量模型应用模型应用 总产出出对外部需求外部需求线性性dd增加增加1个个单位位 x的增量的增量若农业的外部需求增加若农业的外部需求增加1单位单位 x为为 的第的第1列列 6个部门的总产出分别增加个部门的总产出分别增加1.2266，0.5624，0.0075，0.0549，0.0709，0.1325单位单位. 问题问题2 如果如果6个部门的外部需求分别增加个部门的外部需求分别增加1个单位个单位, 问它们的总产出应分别增加多少？问它们的总产出应分别增加多少？求解求解其余外部需求增加其余外部需求增加1单位单位 x为为 的其余各列的其余各列 隶男绍哭权穷馋我区泊复纷坎梆灶抉椅孽嫌什铁睡亦袄屑览识姆淫邓纫碗投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量6.2 CT技术的图像重建技术的图像重建 CT(计算机断层成像计算机断层成像 )技术是技术是20世纪世纪50至至70年代由年代由美国科学家科马克和英国科学家豪斯费尔德发明的美国科学家科马克和英国科学家豪斯费尔德发明的. 1971年第一代供临床应用的年第一代供临床应用的CT设备问世设备问世. 螺旋式螺旋式CT机等新型设备被医疗机构普遍采用机等新型设备被医疗机构普遍采用. CT技术在工业无损探测、资源勘探、生态监测技术在工业无损探测、资源勘探、生态监测等领域也得到了广泛的应用等领域也得到了广泛的应用. 背景背景什么是什么是CT，它与传统的，它与传统的X射线成像有什么区别？射线成像有什么区别？ 伎陛琴扮桅偿酌鱼瞥三折诗拈抢虚弟听那抚观姥匪袄睡迈擒秽滇树殿班札投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量光源光源人眼人眼光源光源人眼人眼一个半透明物体嵌入一个半透明物体嵌入5个不同透明度的球个不同透明度的球 概念图示概念图示 单方向观察无法确定单方向观察无法确定球的数目和透明度球的数目和透明度 让物体旋转从多角度观察能让物体旋转从多角度观察能分辨出分辨出5个球及各自的透明度个球及各自的透明度 人体内脏人体内脏胶胶片片 传统的传统的X射线成像原理射线成像原理 CT技术原理技术原理 探探测测器器 X射射线线X光光管管人体内脏人体内脏CT技术技术: 在在不同深度不同深度的断面上的断面上,从从各个角度各个角度用探用探测器器接收旋转的接收旋转的X光管发出、穿过人体而使光管发出、穿过人体而使强度衰减的射强度衰减的射线线;经过经过测量和计算测量和计算将人体器官和组织的将人体器官和组织的影像影像重新重新构建构建. 图像图像重建重建 账健游孟筑魄恋潮呈吟废饯杨徒杖如谴闹邹糊扩思村菏漾椭稀仙唇锭筒朴投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量X射线强度衰减与图像重建的数学原理射线强度衰减与图像重建的数学原理 射线强度的衰减射线强度的衰减率与强度成正比率与强度成正比. I射线强度射线强度 l物质在射线方向的厚度物质在射线方向的厚度 物质对射线的衰减系数物质对射线的衰减系数 I0入射入射强强度度 射线沿直线射线沿直线L穿行穿行, 穿过由穿过由不同衰减系数的物质组成的不同衰减系数的物质组成的非均匀物体非均匀物体(人体器官人体器官).I0L0yx(x, y)款卸艘放肤莫啸摘平曼衰肇扩骡某苟漠曰裕图下况冕颧睛蒙特伊拥汹澄尹投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量X射线强度衰减与图像重建的数学原理射线强度衰减与图像重建的数学原理 右端数值可从右端数值可从CT 的测量数据得到的测量数据得到 多条直线多条直线L的线积分的线积分 被积函数被积函数(x, y) FQ(q)与与Q相距相距q的直的直线L的的线积分分Pf(L)对所有所有q的平均的平均值 拉东变换拉东变换 拉东逆变换拉东逆变换 图像图像重建重建反映人体器官大小、形状、密度的图像反映人体器官大小、形状、密度的图像 数学数学原理原理实际上只能在有限条直线上得到投影实际上只能在有限条直线上得到投影(线积分线积分). 图像重建在数学方法上的进展，为图像重建在数学方法上的进展，为CT技术在各个技术在各个领域成功的和不断拓广的应用提供了必要条件领域成功的和不断拓广的应用提供了必要条件. 苫青缔吩票房胜翠门教添涅恢棱腑讣们捍带凿寅以恳损穿碘镀侯屋孔瓜饺投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量图像重建的代数模型图像重建的代数模型 lj每个像素对射线的衰减系数是常数每个像素对射线的衰减系数是常数 m个像素个像素(j=1, m), n束射线束射线(i=1,n) Li的的强强度度测量量数据数据 j像素像素j的衰减系数的衰减系数 lj射射线在像素在像素j中的穿行中的穿行长度度 J(Li)射射线Li穿穿过的像素的像素j的集合的集合 像素像素j 射线射线Li瞩逊昔仍纶晴整喷又乾耙胆邓镣帕履人俏耸扯秘乖涎异齿肃源郎舞辰颅矮投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量jLilij图像重建的代数模型图像重建的代数模型 常用算法常用算法 设像素的边长和射线的宽度均为设像素的边长和射线的宽度均为 中心线法中心线法 aij射射线Li的中心的中心线在像素在像素j内的内的长度度lij与与之比之比.面积法面积法 aij射射线Li的中心的中心线在像素在像素j内的面内的面积sij与与之比之比.sij中心法中心法 aij=1射射线Li经过像素像素j的中心点的中心点.档迪喀彩寐扇回慧擒旷孰虑悲跳疥瘟谦雷磕粒枪氛渠蛋峦填潍铝五碌江斥投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量图像重建的代数模型图像重建的代数模型 中心法的简化形式中心法的简化形式 假定射线的宽度为零假定射线的宽度为零, 间距间距 aij=1 Li经过像素像素j内任一点内任一点987654321L4L8L7L6L5L3L2L1根据根据A和和b, 由由 确定像素的衰减系数向量确定像素的衰减系数向量x m和和n很大且很大且m n, 方程有无穷多解方程有无穷多解 + 测量误差和噪声测量误差和噪声 在在x和和e满足的足的最优准则下估计最优准则下估计x 代数重建技术代数重建技术(ART)羌巡枫钩噶寂泅靴闭麦产峭失羞沾耸断只耳卒挂惹惜宛忧逢实绍堂努幼湖投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量6.3 原子弹爆炸的能量估计原子弹爆炸的能量估计1945年年7月月16日美国科学家在新墨西哥州阿拉莫日美国科学家在新墨西哥州阿拉莫戈多沙漠试爆了全球第一颗原子弹戈多沙漠试爆了全球第一颗原子弹, 震惊世界震惊世界!当时资料是保密的当时资料是保密的, 无法准确估计爆炸的威无法准确估计爆炸的威力力.英国物理学家泰勒研究了两年后美国公开的录像带英国物理学家泰勒研究了两年后美国公开的录像带, 利用利用数学模型估计数学模型估计这次爆炸释放的能量为这次爆炸释放的能量为19.2千吨千吨. 后来公布爆炸实际后来公布爆炸实际释放的能量释放的能量21千吨千吨 鼠釜赵莹氢捎奖俗近窝蚌销依请礼逸痕亡演粤奴更篇家所削戌卒春官蜘输投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量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原子弹爆炸的能量估计原子弹爆炸的能量估计爆炸产生的冲击波以爆炸点为中心呈球面向四周传播爆炸产生的冲击波以爆炸点为中心呈球面向四周传播,爆炸的能量越大，在一定时刻冲击波传播得越远爆炸的能量越大，在一定时刻冲击波传播得越远. 冲击波由爆炸形成的冲击波由爆炸形成的“蘑菇云蘑菇云”反映出来反映出来. 泰勒用泰勒用量纲分析方法量纲分析方法建立数学模型建立数学模型, 辅以小型试验辅以小型试验,又利用测量数据对爆炸的能量进行估计又利用测量数据对爆炸的能量进行估计.瑶汇灵歪踩尊遂跨宝赐产研绳凶烤渴孪底匈出斧移惟布酉砌直审蛾彩织绘投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量物物理理量量的的量量纲纲长度长度 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在经验和实验的基础上利用物理定律的量纲齐次原则在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系确定各物理量之间的关系. 痹链熟丈执胡成驮昭缩屯禄押坤呸枝鞋敬署嫂骗砧遭襄毖檬屡滓钞九眶跃投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量量纲齐次原则量纲齐次原则等式两端的量纲一致等式两端的量纲一致量纲分析量纲分析利用量纲齐次原则寻求物理量之间的关系利用量纲齐次原则寻求物理量之间的关系.例：单摆运动例：单摆运动lmgm求摆动周期求摆动周期 t 的表达式的表达式设物理量设物理量 t, m, l, g 之间有关系式之间有关系式 1, 2, 3 为待定系数，为待定系数， 为无量纲量为无量纲量 (1)的量纲表达式的量纲表达式与与 对比对比捍碟灿萤淤节配警锣啊堆巧奠敌梅肃港勉骑氏惹乞材奥未凿吁剁茹泻榔溶投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量对对 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倍倍p= f(x,y,z)的形式为的形式为量纲齐量纲齐次原则次原则单摆运动单摆运动愚亏高擒靶跨咒捐翔突箭峡络仇芭扯绪制转氮燃四那萧涧幂示荐咬峦晴溶投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量单摆运动中单摆运动中 t, m, l, g 的一般表达式的一般表达式y1y4 为待定常数为待定常数, 为无量纲量为无量纲量基本解基本解腊眨吸梧绍摔蝗啸粘您陨莆圈坯蛇鸟失掷麻得混挨苫缮谷拍灌妄人棱榴丙投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量设设 f(q1, q2, , qm) = 0 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 的量纲可表为的量纲可表为量纲矩阵记作量纲矩阵记作线性齐次方程组线性齐次方程组有有 m-r 个基本解，记作个基本解，记作为为m-r 个相互独立的无量纲量个相互独立的无量纲量, 且且则则搭镜挛撇帛施酪芽莲硕诞叁禁廷绞丈桅疹郧昼蚜僵败千泉嚎讽邀该有握掩投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量记爆炸能量为记爆炸能量为E，将，将“蘑菇云蘑菇云”近似看成一个球近似看成一个球形形.时刻时刻 t 球的半径为球的半径为 rt, E空气密度空气密度, 大气压强大气压强P基本量纲：基本量纲：L, M, T 原子弹爆炸能量估计的量纲分析方法建模原子弹爆炸能量估计的量纲分析方法建模 r与哪些因素有关？与哪些因素有关？ r t E P LMT量纲矩阵量纲矩阵 辫榔鹰揖础激撵炯堆靶日赣羌襄急妻子程钎养杰凶搂蔓吁炕属酝乞掩寒镑投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量y=(1,-2/5,-1/5,1/5,0) y=(0,6/5,-2/5,-3/5,1)T原子弹爆炸能量估计的量纲分析方法建模原子弹爆炸能量估计的量纲分析方法建模 有有2个基本解个基本解两个无两个无量纲量量纲量吝产汪酞拴二蛙承扬亚虾砖寿凿沟食弱惦跋埂芳钦跳枉暇稳假蛋墅辑座犹投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量原子弹爆炸能量估计的数值计算原子弹爆炸能量估计的数值计算时间时间 t 非常短非常短能量能量 E 非常大非常大 泰勒根据一些小型爆炸试验的数据建议泰勒根据一些小型爆炸试验的数据建议用用r, t 的实际数据做平均的实际数据做平均空气密度空气密度 =1.25 (kg/m3) 1千吨千吨(TNT能量能量)= 4.184*1012焦尔焦尔 E=19.7957 (千吨千吨)E=8.28251013(焦耳焦耳)实际值实际值21千吨千吨 傈撅幕牧穷赶耸惧隐增赂英嘻丹帆安拉垫契姥腥严枯枚亲湖雏鞍沉秤堤滴投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量泰勒的计算泰勒的计算tr最小二乘法拟合最小二乘法拟合 r=atbE=8.02761013 (焦耳焦耳)即即19.2千吨千吨 取取y平均值得平均值得c=6.9038 模型检验模型检验b=0.40582/5肖吃卉尔魄长琅涩迈鳞眩词迭纱洋韭毫肆迄整澜肋饮豌铁哺耙砷蛰揽柞碘投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量量纲分析法的评注量纲分析法的评注 物理量的选取物理量的选取 基本量纲的选取基本量纲的选取 基本解的构造基本解的构造 结果的局限性结果的局限性 () = 0中包括哪些物理量是至关重要的中包括哪些物理量是至关重要的.基本量纲个数基本量纲个数n; 选哪些基本量纲选哪些基本量纲.有目的地构造有目的地构造 Ay=0 的基本解的基本解. 方法的普适性方法的普适性函数函数F和无量纲量未定和无量纲量未定.不需要特定的专业知识不需要特定的专业知识.傅铰龋溃调侩杏奏哭婪汾胡果扔潘缀匿译嫉伦朔陡沥颠堪椽蒂临置摧莉免投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量物理模拟示例：物理模拟示例：波浪对航船的阻力波浪对航船的阻力航船阻力航船阻力 f航船速度航船速度v, 船体尺寸船体尺寸l, 浸没面积浸没面积 s, 海水密度海水密度 , 重力加速度重力加速度g .量纲分析在物理模拟中的应用量纲分析在物理模拟中的应用 物理模拟物理模拟: 按照一定的按照一定的比例尺寸比例尺寸构造它的物理模型构造它的物理模型,通过对模型的研究得出原型的结果通过对模型的研究得出原型的结果.量纲分析可以指导物理模拟中量纲分析可以指导物理模拟中比例尺寸比例尺寸的确定的确定. 轩勉磐甩刻巩辆坊尤喂梁斌申替审以近懂乳挤体敷弘拄俊碌设擞淤茬硒佬投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量物理模拟示例：物理模拟示例：波浪对航船的阻力波浪对航船的阻力定理定理 原型船原型船模型船模型船 模型船的模型船的 均已知均已知 当原型船的当原型船的 给定后计算给定后计算 f物理物理模拟模拟肉型钉霜囱配嘲瑞蓄拒儡灾悔埂讨藐监褐淄吐济传庚洛壮征情袋胃吩相坑投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量物理模拟示例：物理模拟示例：波浪对航船的阻力波浪对航船的阻力原原型型船船模模型型船船模拟模拟条件条件量测模型船阻力量测模型船阻力f，可计算，可计算 f. 按一定尺寸比例建造按一定尺寸比例建造模型船模型船, 并调节船速并调节船速.摘华什腊凤迄痢某喧堵锦吊还养搭俄屯嚼亭俺嘻跃蔽耳痒刑坍荆峙衅就猜投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量无量纲化示例：火箭发射无量纲化示例：火箭发射m1m2xrv0g星球表面竖直发射火箭。初速星球表面竖直发射火箭。初速v, 星星球半径球半径r, 星球表面重力加速度星球表面重力加速度g.研究火箭高度研究火箭高度 x 随时间随时间 t 的变化规律的变化规律.t=0 时时 x=0, 火箭质量火箭质量m1, 星球质量星球质量m2牛顿第二定律，万有引力定律牛顿第二定律，万有引力定律3个独立参数个独立参数舒符瞬陵怔敲钻译打潮般迟唇柯纲鸦釉式怀贿盲拈毗杯泞晃眶晋少夷些评投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量用无量纲化方法减少独立参数个数用无量纲化方法减少独立参数个数x=L, t=T, r=L, v=LT-1, g=LT-2变量变量 x,t 和独立参数和独立参数 r,v,g 的量纲的量纲用用参数参数r,v,g的组合的组合, ,分别分别构造与构造与x,t具有相同具有相同量纲量纲的的xc, tc （特征尺度）（特征尺度）无量纲变量无量纲变量如如利用新变量利用新变量将被简化将被简化令令淳蒙墩啮侨他房搏悍汛悸印糕梳锦雅帽识馁贮呐狱纺尉喷粗毫辩捞雇禹柒投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 xc, tc的不同构造的不同构造1）令 为无量纲量为无量纲量用无量纲化方法减少独立参数个数用无量纲化方法减少独立参数个数的不同简化结果的不同简化结果凭芒暑寝获索瞒蓄樟地匣诞钓笑疚斥谐矿澡输短黄诗恍汕旁她约防溢坝浑投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量3）令 为无量纲量为无量纲量2）令 为无量纲量为无量纲量用无量纲化方法减少独立参数个数用无量纲化方法减少独立参数个数晓春您园巳辆租盲冉利找帮邯酥吁疼断绒埃邹媚腋刨凰绎点纯腋赘着郝浑投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量1) 2) 3) 的共同点的共同点只含只含1个参数个参数无量纲量无量纲量 解解1) 2) 3) 的重要差别的重要差别考察无量纲量考察无量纲量在在1) 2) 3) 中能否忽略以中能否忽略以 为因子的项？为因子的项？1)忽略忽略 项项无解无解不能忽略不能忽略 项项无量纲无量纲化方法化方法添敌搭百戒抡疏贿蟹苔牺灭擅随爸镰处佩码剑杏科医堑拙娶傻变冀益祷畦投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量2)3)忽略忽略 项项不能忽略不能忽略 项项忽略忽略 项项1) 2) 3) 的重要差别的重要差别无量纲无量纲化方法化方法焕纶刚涪脏憾泉拙咸吴秩皿鳖腊载右舟囱邀灼地贷朱吹加轴朋序抑瞧钉尘投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量火箭发射过程火箭发射过程中引力中引力m1g不变不变 即即 x+r r原原问问题题可以忽略可以忽略 项项是原问题是原问题的近似解的近似解1) 2) 3) 的重要差别的重要差别无量纲化方法无量纲化方法襄湖姿妙吕泼佳缘筒夷纱专懈驶控憨渺斗奴癸铂覆忌趴事墅蜂守捅跺写馏投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量为什么为什么3)能忽略能忽略 项，得到原问题近似解，而项，得到原问题近似解，而1) 2)不能不能?1）令）令2）令）令3）令）令火箭到达最高点时间为火箭到达最高点时间为v/g, 高度为高度为v2/2g,大体上具有单位尺度大体上具有单位尺度项可以忽略项可以忽略项不能忽略项不能忽略无量纲化方法无量纲化方法餐友洁授闸佳组纹汤龋越躁底拿昭瘪死答皮坐甩氨慨趋篱巴侨骤譬胯督噬投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 选择特征尺度的一般讨论见：选择特征尺度的一般讨论见：林家翘著林家翘著自然科学中确定性问题的应用数学自然科学中确定性问题的应用数学无无 量量 纲纲 化化 无量纲化无量纲化是研究物理问题常用的数学方法是研究物理问题常用的数学方法. 选择选择特征尺度特征尺度主要依赖于物理知识和经验主要依赖于物理知识和经验. 恰当地选择特征尺度可以减少独立参数恰当地选择特征尺度可以减少独立参数个数，还可以辅助确定舍弃哪些次要因素个数，还可以辅助确定舍弃哪些次要因素.资堰索燕拒津诲魁字直谜汗玖禾实雅孔拂鉴赶樟代畅灭膛险叫杖钢涣赡碍投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量6.4 市场经济中的蛛网模型市场经济中的蛛网模型问问 题题供大于求供大于求现现象象 商品数量与价格的振荡在什么条件下趋向稳定商品数量与价格的振荡在什么条件下趋向稳定? 当不稳定时政府能采取什么干预手段使之稳定当不稳定时政府能采取什么干预手段使之稳定?价格下降价格下降减少产量减少产量增加产量增加产量价格上涨价格上涨供不应求供不应求 描述商品数量与价格的变化规律描述商品数量与价格的变化规律.商品商品数量数量与与价格价格在振荡在振荡募硝幂馏氓丑泌循竿尽绍嗣祈帘驯坡晋疤皇咸赘钡挖踩腻炔碳审乔洁父苯投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量蛛蛛 网网 模模 型型gx0y0P0fxy0xk第第k时段商品数量；时段商品数量；yk第第k时段商品价格时段商品价格.消费者的需求关系消费者的需求关系生产者的供应关系生产者的供应关系减函数减函数增函数增函数需求函数需求函数f与与g的交点的交点P0(x0,y0) 平衡点平衡点一旦一旦xk=x0，则，则yk=y0, xk+1,xk+2,=x0, yk+1,yk+2, =y0 供应函数供应函数缅乓饶议琐芭任生汀诗颇丧棘因船圭瑰威苍挟语撵桓猾蚊沾趋酮承肺凝试投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量xy0fgy0x0P0设设x1偏离偏离x0x1x2P2y1P1y2P3P4x3y3P0是稳定平衡点是稳定平衡点P1P2P3P4P0是不稳定平衡点是不稳定平衡点xy0y0x0P0fg 曲线斜率曲线斜率蛛蛛 网网 模模 型型 使厄洼映芬规茵唯谚虫筏毒拜灾娇甚符担城经面唉她狗同涉遇聚仇要炭殷投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量在在P0点附近用直线近似曲线点附近用直线近似曲线P0稳定稳定P0不稳定不稳定方方 程程 模模 型型方程模型与蛛网模型的一致方程模型与蛛网模型的一致恿砌体浊葵禾豹散泪讯缩肮茄辗藕腊迎独谰泞荡诀术渠葛衰蜜服浇汲乡荧投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 商品数量减少商品数量减少1单位单位, 价格上涨幅度价格上涨幅度 价格上涨价格上涨1单位单位, (下时段下时段)供应的增量供应的增量考察考察 , 的含义的含义 消费者对需求的敏感程度消费者对需求的敏感程度 生产者对价格的敏感程度生产者对价格的敏感程度 小小, 有利于经济稳定有利于经济稳定 小小, 有利于经济稳定有利于经济稳定结果解释结果解释xk第第k时段商品数量；时段商品数量；yk第第k时段商品价格时段商品价格.经济稳定经济稳定结果解释结果解释篇囚战仟俄澄农肉赠搬敲使植赚拙概挤驹笺谗褐俞陪练粟获痛某量甄纪错投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量经济不稳定时政府的干预办法经济不稳定时政府的干预办法1. 使使 尽量小，如尽量小，如 =0 以行政手段控制价格不变以行政手段控制价格不变2. 使使 尽量小，如尽量小，如 =0靠经济实力控制数量不变靠经济实力控制数量不变xy0y0gfxy0x0gf结果解释结果解释需求曲线变为水平需求曲线变为水平供应曲线变为竖直供应曲线变为竖直倾方挺看要粒慢谈顽僚猜湿贪础慰入篓筹怒能恭棕舷趟肉霸变畏报侵漳猎投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量模型的推广模型的推广 生产者根据当前时段和前一时段生产者根据当前时段和前一时段 的价格决定下一时段的产量的价格决定下一时段的产量.生产者管理水平提高生产者管理水平提高设供应函数为设供应函数为需求函数不变需求函数不变二阶线性常系数差分方程二阶线性常系数差分方程x0为平衡点为平衡点研究平衡点稳定，即研究平衡点稳定，即k, xkx0的条件的条件沏警找芹卧琉每杭阿制搓战恶圾成夯偿峻炭菊斗错线切絮差拐酬违皑惯爬投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量方程通解方程通解(c1, c2由初始条件确定由初始条件确定) 1, 2特征根，即方程特征根，即方程 的根的根 平衡点稳定，即平衡点稳定，即k, xkx0的条件的条件:平衡点稳定条件平衡点稳定条件比原来的条件比原来的条件 放宽了放宽了!模型的推广模型的推广函嗡趾氢逝跪虱南宏挡舅毖种高蓟枷岁跑赢乓济娥嫉摧判兴怕诸锤闷谗灸投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量6.5 减肥计划减肥计划节食与运动节食与运动背背景景 多数减肥食品达不到减肥目标，或不能维持多数减肥食品达不到减肥目标，或不能维持. 通过控制饮食和适当的运动，在不伤害身体通过控制饮食和适当的运动，在不伤害身体 的前提下，达到减轻体重并维持下去的目标的前提下，达到减轻体重并维持下去的目标.分分析析 体重变化由体内能量守恒破坏引起体重变化由体内能量守恒破坏引起. 饮食（吸收热量）引起体重增加饮食（吸收热量）引起体重增加. 代谢和运动（消耗热量）引起体重减少代谢和运动（消耗热量）引起体重减少. 体重指数体重指数BMI=w(kg)/l2(m2). 18.5BMI25 超重超重; BMI30 肥胖肥胖.唇菌员滑赤傣窝氛祖舜信液枫涵剐肉劲垂况匙病方豪拭妊簿丘洛联纬毡邑投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量模型假设模型假设1）体重增加正比于吸收的热量体重增加正比于吸收的热量每每8000千卡千卡 增加体重增加体重1千克；千克；2）代谢引起的体重减少正比于体重代谢引起的体重减少正比于体重每周每千克每周每千克 体重消耗体重消耗200千卡千卡 320千卡千卡(因人而异因人而异), 相当于相当于70 千克的人每天消耗千克的人每天消耗2000千卡千卡 3200千卡；千卡；3）运动引起的体重减少正比于体重，且与运动形式）运动引起的体重减少正比于体重，且与运动形式 有关；有关； 4）为了安全与健康，每周体重减少不宜超过）为了安全与健康，每周体重减少不宜超过1.5千克千克, 每周吸收热量不要小于每周吸收热量不要小于10000千卡千卡.闹臆亚旦处敞符托捅骆非碟芍联暑曙懂汐寒丰申租砸夕团予吼旦刘盎症两投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量某甲体重某甲体重100千克，目前每周吸收千克，目前每周吸收20000千卡热量，千卡热量，体重维持不变。现欲减肥至体重维持不变。现欲减肥至75千克千克.第一阶段：每周减肥第一阶段：每周减肥1千克，每周吸收热量逐渐减少千克，每周吸收热量逐渐减少,直至达到下限（直至达到下限（10000千卡）；千卡）；第二阶段：每周吸收热量保持下限，减肥达到目标第二阶段：每周吸收热量保持下限，减肥达到目标. 2）若要）若要加快进程加快进程，第二阶段增加运动，试安排计划，第二阶段增加运动，试安排计划.1）在不运动的情况下安排一个）在不运动的情况下安排一个两阶段计划两阶段计划.减肥计划减肥计划3）给出达到目标后）给出达到目标后维持体重维持体重的方案的方案.抗预汾官族辞燎鲍夏竟蒸瓶扳寺调腾结睫如蝇导舷湾盛阴呢时拘速蔽怎娄投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 确定某甲的代谢消耗系数确定某甲的代谢消耗系数即每周每千克体重消耗即每周每千克体重消耗 20000/100=200千卡千卡基本模型基本模型w(k) 第第k周周(末末)体重体重c(k) 第第k周吸收热量周吸收热量 代谢消耗系数代谢消耗系数(因人而异因人而异)1）不运动情况的两阶段减肥计划）不运动情况的两阶段减肥计划每周吸收每周吸收20000千卡千卡 w=100千克不变千克不变 =1/8000(千克千克/千卡千卡)毕茬财膜梅摧川羌边梅拉直硷帧唐贩唁练树痔客城卞口版鹿赖肾鹏炼料吉投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 第一阶段第一阶段: w(k)每周减每周减1千克千克, c(k)减至下限减至下限10000千千卡卡第一阶段第一阶段10周周, 每周减每周减1千克，第千克，第10周末体重周末体重90千克千克吸收热量为吸收热量为1）不运动情况的两阶段减肥计划）不运动情况的两阶段减肥计划斌信陕全蹋样痹底谅作争掩灿岁娩濒洱景浴综戳拍俱隐务惜牧丢背缉署绳投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 第二阶段：每周第二阶段：每周c(k)保持保持Cm, w(k)减至减至75千克千克1）不运动情况的两阶段减肥计划）不运动情况的两阶段减肥计划基本模型基本模型做惜氟噬坷屹捞姑蒜佰整宽泄靛唆坠逃壳噪徐嘎吮弥俗助利孝玻讫壶养拖投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 第二阶段：每周第二阶段：每周c(k)保持保持Cm, w(k)减至减至75千克千克第二阶段第二阶段19周周, 每周吸收热量保持每周吸收热量保持10000千卡千卡, 体重按体重按 减少至减少至75千克千克.阶域愁卖拥乐捍椒圃捂盲玲侣养朋趾却氧孙默靡操窍蚊省竞丁凹党夸仔健投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量运动运动 t=24 (每周每周跳舞跳舞8小时或自行车小时或自行车10小时小时), 14周即可周即可.2）第二阶段增加运动的减肥计划）第二阶段增加运动的减肥计划根据资料每小时每公斤体重消耗的热量根据资料每小时每公斤体重消耗的热量 (千卡千卡): 跑步跑步 跳舞跳舞 乒乓乒乓 自行车自行车(中速中速) 游泳游泳(50米米/分分) 7.0 3.0 4.4 2.5 7.9t每周运动时间每周运动时间(小小时时)模型模型)() 1()() 1(kwkckwkw - -+ + += =+ +t + += = 取取t=0.003, 即即 t=24 =1/8000(千克千克/千卡千卡), =0.02514=n巩枝旧棋草国逐滁峭阂拟促吩挝陵蝶阀舅陪椎淋感桅郊瓷敌喉迷甭翼狮缓投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量增加运动相当于提高代谢消耗系数增加运动相当于提高代谢消耗系数 2）第二阶段增加运动的减肥计划）第二阶段增加运动的减肥计划提高提高12%减肥所需时间从减肥所需时间从19周降至周降至14周周减少减少25% 这个模型的结果这个模型的结果对代谢消耗系数对代谢消耗系数 很敏感很敏感. 应用该模型时要仔细确定代谢消耗系数应用该模型时要仔细确定代谢消耗系数 (对不同的人对不同的人; 对同一人在不同的环境对同一人在不同的环境).刽狂魄泽两龄媚戈埋排咸亚系竞意幌劝润疾奠兰肆箍距每块希烽秉怖鞭扮投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量3）达到目标体重）达到目标体重75千克后维持不变的方案千克后维持不变的方案每周吸收热量每周吸收热量c(k)保持某常数保持某常数C，使体重，使体重w不变不变 不运动不运动 运动运动(内容同前内容同前)睬喷励锭霜额分榜斌腻辆囱排呀鄂鼠第桂削轴峻认宿篷初顶十垫激涪瞅噪投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量6.6 按年龄分组的种群增长按年龄分组的种群增长 不同年龄组的繁殖率和死亡率不同不同年龄组的繁殖率和死亡率不同.建立差分方程模型，讨论稳定状况下种群的增长规律建立差分方程模型，讨论稳定状况下种群的增长规律.假设与建模假设与建模 种群按年龄大小等分为种群按年龄大小等分为n个年龄组，记个年龄组，记i=1,2,n 时间离散为时段，长度与年龄组区间相等，记时间离散为时段，长度与年龄组区间相等，记k=1,2, 以雌性个体数量为对象以雌性个体数量为对象. 第第i 年龄组年龄组1雌性个体在雌性个体在1时段内的时段内的繁殖率繁殖率为为bi 第第i 年龄组在年龄组在1时段内的死亡率为时段内的死亡率为di, 存活率存活率为为si=1- di幂破郁藤轨二郡级赁卡姬午拣鲜迄笋缮淡湛其智挤畦摇瘟醇琢谷屿偿澈奔投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量假设假设与与建模建模xi(k)时段时段k第第i 年龄组的种群数量年龄组的种群数量按年龄组的分布向量按年龄组的分布向量预测任意时段种群预测任意时段种群按年龄组的分布按年龄组的分布Leslie矩阵矩阵(L矩阵矩阵)(设至少设至少1个个bi0)茬肤求痊嘿隐徊告伦淀买掷秀妮森缕诅沃钡鞭俞幅命坯辩胀险源橡臭钧矛投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量稳定状态分析的数学知识稳定状态分析的数学知识 L矩阵存在矩阵存在正单特征根正单特征根 1， 若若L矩阵存在矩阵存在bi, bi+10, 则则 P的第的第1列是列是x*特征向量特征向量, c是由是由bi, si, x(0)决定的常决定的常数数 且且解解释释L对角化对角化低跋儡朔磐噪盛扬顿酋贰仁又纬蛋圾抚孝枝那层湛色抨案惕咀宅艺邱蚁衬投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量稳态分析稳态分析k充分大充分大种群按年龄组的分布种群按年龄组的分布 种群按年龄组的分布趋向稳定，种群按年龄组的分布趋向稳定，x*称称稳定分布稳定分布, 与初始分布无关。与初始分布无关。 各年龄组种群数量按同一各年龄组种群数量按同一倍数增减，倍数增减， 称称固有增长率固有增长率与基本模型与基本模型比较比较3） =1时时 各年龄组各年龄组种群种群数量不变数量不变二秦割耀巢盐不嘘饼孵劣欺滇肚罗耪日递您潜口时褪沧雾洋汲哲咙卵仕噬投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量 1个个体在整个存活期个个体在整个存活期 内的繁殖数量为内的繁殖数量为1稳态分析稳态分析存活率存活率 si是同一时段的是同一时段的 xi+1与与 xi之比之比（与（与si 的定义的定义 比较）比较） 3） =1时时苫有蛋酣愧念遗意芍瘦堰呵刃冈莽涅址盘圃临喷眉双韦宠营稀赛碧榆珍抢投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量人口模型人口模型连续型人口模型的离散形式连续型人口模型的离散形式xi(k)k年年i 岁的女性人数岁的女性人数(模型只考虑女性人口模型只考虑女性人口).bi(k)k年年i 岁女性生育率岁女性生育率(每人平均生育女儿数每人平均生育女儿数).dii 岁女性死亡率，岁女性死亡率，si=1-di存活率存活率i1, i2生育区生育区间间k年育龄女性平均生育女儿数年育龄女性平均生育女儿数总合生育率总合生育率(生育胎次生育胎次)年龄分布向量年龄分布向量hi生育模式生育模式傍尘狭鹿嚏治舷诽裔嫌家套诛蚜挫丸壶浇攻迎椰菠层滓腺妹彝廖秆靳沫钓投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量人口模型人口模型存活率矩阵存活率矩阵生育模式矩阵生育模式矩阵x(k)状态变量状态变量, (k)控制变控制变量量双线性方程双线性方程(对对x(k), (k)线性线性)原模型原模型状格扼鬼限顺锣藕押叛寂与民充塌慑桩别铆丽距拼隙坝拽捶芹勤既炯涎答投入产出模型CT技术的图像重建原子弹爆炸的能量投入产出模型CT技术的图像重建原子弹爆炸的能量