反问题概述.ppt

1偏微分方程反问题数值解法偏微分方程反问题数值解法 韩 波主要参考书主要参考书：1、刘继军，不适定问题的正则化方法及应用，科学出版社，20052、肖廷延，于慎根，王彦飞，反问题的数值解法，科学出版社，2003，3、韩波，李莉，非线性不适定问题的求解方法及其应用，科学出版社，2011大家好2第一章第一章 反问题概述反问题概述大家好3从墨西哥湾的漏油事件谈起从墨西哥湾的漏油事件谈起 2010年年4月月20日晚，在距美国路易斯安娜州日晚，在距美国路易斯安娜州61公里的墨公里的墨西哥湾海面上，由西哥湾海面上，由BP（英国石油公司）租用的（英国石油公司）租用的“深海地平深海地平线线”(DeepwaterHorizon)钻井平台在从事钻井平台在从事Macondo油田开油田开发作业时发生爆炸，并引发火灾，发作业时发生爆炸，并引发火灾，11名工作人员遇难。经名工作人员遇难。经过约过约36个小时的剧烈燃烧，个小时的剧烈燃烧，“深海地平线深海地平线”于于22日沉入海日沉入海底。两天后，在事故地点出现了严重的石油泄漏，破裂的底。两天后，在事故地点出现了严重的石油泄漏，破裂的油井管道每天向墨西哥湾海域倾倒油井管道每天向墨西哥湾海域倾倒5-10万桶原油（约万桶原油（约8000-16000立方米），历时近三个月，造成了极为严重的环境立方米），历时近三个月，造成了极为严重的环境污染。直至污染。直至7月月15日，油井才被成功封堵。这是美国历史日，油井才被成功封堵。这是美国历史上迄今最严重的海上石油泄露事件。上迄今最严重的海上石油泄露事件。大家好4当关注的焦点聚集到英国石油公司时，我们发当关注的焦点聚集到英国石油公司时，我们发现这个公司在墨西哥湾地区寻找油气藏方面取得了现这个公司在墨西哥湾地区寻找油气藏方面取得了巨大的成功。巨大的成功。2009年年9月月2日，英国石油公司在墨西哥湾距休日，英国石油公司在墨西哥湾距休斯敦市东南斯敦市东南400公里处发现了巨型油田公里处发现了巨型油田Tiber。外界。外界估计其总储量可达估计其总储量可达40-60亿桶，相当于约亿桶，相当于约100个个Macondo。英国石油公司在官方网站上例举了他们近几年英国石油公司在官方网站上例举了他们近几年在墨西哥湾深海地区开发的重要油田，预测产量都在墨西哥湾深海地区开发的重要油田，预测产量都在在1.5亿桶至亿桶至6亿桶之间。亿桶之间。石油巨头们是怎么在墨西哥湾找到这些大石油巨头们是怎么在墨西哥湾找到这些大油田的？油田的？大家好5深海石油开采的深海石油开采的费用高的惊人，即使从事深海石用高的惊人，即使从事深海石油开油开发的公司的公司拥有雄厚的有雄厚的资本和先本和先进的技的技术，也，也不可能在茫茫无不可能在茫茫无际的墨西哥湾到的墨西哥湾到处打井打井窥探。探。ThunderHorse油田：水上油田：水上钻井平台有井平台有3个足球个足球场大，建筑在海水深度大，建筑在海水深度约为1800米的海面上。据估米的海面上。据估计英国石油公司投入英国石油公司投入该项目的建造目的建造费用高达用高达50亿美元！美元！对资源的定位和估算有一套源的定位和估算有一套较准确的科学方法，准确的科学方法，尽量减少决策失尽量减少决策失误。大家好6秘密就在回声中数学物理史上有趣的数学物理史上有趣的问题：不用眼睛看，不用眼睛看，仅仅通通过聆听鼓的声音能否判断聆听鼓的声音能否判断出鼓的形状？（盲人听鼓）出鼓的形状？（盲人听鼓）1910年丹麦著名物理学家年丹麦著名物理学家劳伦兹(Lorentz)在哥廷根的系列演在哥廷根的系列演讲“物理学中物理学中的新、旧的新、旧问题”中提出中提出。大家好7当物体的材料确定后，它的音色和其形状密切相关。在数学上，一个物体的音色可以当物体的材料确定后，它的音色和其形状密切相关。在数学上，一个物体的音色可以由一串由一串谱来确定，它来确定，它们对应着物体的固有着物体的固有频率。率。“盲人听鼓盲人听鼓”即是要求通即是要求通过已知的已知的谱来确定一个来确定一个鼓面的形状。鼓面的形状。Lorentz在他的在他的讲演中猜演中猜测鼓的面鼓的面积可以由极限确定可以由极限确定据据说，当初，当初Hilbert（希（希尔伯特）伯特）认为在他的有生之年不可能看到在他的有生之年不可能看到这个公式的个公式的严格格证明。但是一代宗明。但是一代宗师Hilbert这次作出了次作出了错误的的预测。不到两年。不到两年时间，鼓面，鼓面积的公式就被他的的公式就被他的得意得意门生生Weyl证明了。而且明了。而且证明方法采用的正是明方法采用的正是Hilbert此前不久修此前不久修炼出的独出的独门绝技技-积分理分理论。1954年，年，Pleijel证明了从鼓声中可以明了从鼓声中可以“听听”出鼓的周出鼓的周长。1967年，年，McKean和和Singer证明了从鼓声中可以明了从鼓声中可以“听听”出鼓的内部是否有洞、有几个洞。直到出鼓的内部是否有洞、有几个洞。直到1992年，年，Gordon等人构造等人构造出了两面奇怪的出了两面奇怪的“同声鼓同声鼓”：它：它们的形状不同，却有着相同的音色，的形状不同，却有着相同的音色，单凭耳凭耳朵无法无法鉴别！大家好8大家好9 严格说来，严格说来，“盲人听鼓盲人听鼓”问题问题的答案是否定的。但是，对这个问题的研的答案是否定的。但是，对这个问题的研究启发了我们。当不能用眼睛直接观测时，究启发了我们。当不能用眼睛直接观测时，以耳代目也能够获得关于物体形状的很多以耳代目也能够获得关于物体形状的很多有用信息。举一个生活中的例子，夏天人有用信息。举一个生活中的例子，夏天人们挑西瓜，总是把瓜放在耳边，用手拍一们挑西瓜，总是把瓜放在耳边，用手拍一拍，有经验的人就知道瓜瓤熟不熟。拍，有经验的人就知道瓜瓤熟不熟。大家好10 深海区的石油探测就是应用了类似的原理。勘探地球物深海区的石油探测就是应用了类似的原理。勘探地球物理学家希望能够叩问地球，用耳朵理学家希望能够叩问地球，用耳朵“听听”出地下的地质构出地下的地质构造，从而判断出油藏的造，从而判断出油藏的“准确准确”位置和储量。位置和储量。大家好11 数据采集船上带有气枪。当压缩空气被突数据采集船上带有气枪。当压缩空气被突然释放时，气枪会产生剧烈的爆炸声波。然释放时，气枪会产生剧烈的爆炸声波。声波向地下传播，遇到构造变化会产生反声波向地下传播，遇到构造变化会产生反射、散射和折射。这些回声中携带了地下射、散射和折射。这些回声中携带了地下的地质信息，被海面采集船拖带的检波器的地质信息，被海面采集船拖带的检波器接收，记录为地震数据。海底宝藏的秘密接收，记录为地震数据。海底宝藏的秘密就隐藏在这些数据里。就隐藏在这些数据里。大家好12大家好13大家好14现行的勘探技术主要分为三步：现行的勘探技术主要分为三步：1.数据处理数据处理：海上采集到的数据有很多噪音干扰，海上采集到的数据有很多噪音干扰，比如海浪、背景噪声、相邻采集船传来的气枪回比如海浪、背景噪声、相邻采集船传来的气枪回声，等等。另外，潮汐现象周期性地改变海水的声，等等。另外，潮汐现象周期性地改变海水的深度，造成反射时间的偏差错落；气枪的震源效深度，造成反射时间的偏差错落；气枪的震源效应和海面多次反射的混声也会引起反射信号的畸应和海面多次反射的混声也会引起反射信号的畸变。所以我们要对数据进行变。所以我们要对数据进行“整容整容”，使之更清，使之更清晰、更干净、更可信。晰、更干净、更可信。大家好152.速度建模速度建模：利用多次覆盖的试验数据来构建一利用多次覆盖的试验数据来构建一个统一的地下模型。由于反射波的传播时间对速个统一的地下模型。由于反射波的传播时间对速度的短周期变化不敏感，一步到位求解度的短周期变化不敏感，一步到位求解v(x,y,z)太太困难。退而求其次，我们可以利用数据中的冗余困难。退而求其次，我们可以利用数据中的冗余信息构造一个粗眉大眼的近似速度。这如同绘画，信息构造一个粗眉大眼的近似速度。这如同绘画，先勾勒轮廓，敷设底色，留待下一步皴擦点染，先勾勒轮廓，敷设底色，留待下一步皴擦点染，刻画雕琢。刻画雕琢。大家好163.全波形反演成像全波形反演成像 在较好的初始速度模型基础上，直接利用经过预在较好的初始速度模型基础上，直接利用经过预处理的勘探资料，求取波动方程中的速度模型，处理的勘探资料，求取波动方程中的速度模型，给出地质构造解释。给出地质构造解释。目前，国际上已经能够进行三维速度全波形反演。目前，国际上已经能够进行三维速度全波形反演。弹性波全波形反演研究已经成为国际潮流。弹性波全波形反演研究已经成为国际潮流。大家好17正问题正问题，一般是按着自然顺序来研究事物的一般是按着自然顺序来研究事物的演化过程或分布形态，起着由因推果的作用。演化过程或分布形态，起着由因推果的作用。反问题反问题，是根据事物的演化结果，由可观测的现是根据事物的演化结果，由可观测的现象来探求事物的内部规律或所受的外部影响，起象来探求事物的内部规律或所受的外部影响，起着倒果求因的作用。着倒果求因的作用。大家好18线性方程线性方程 Ax=yAx=yy y=y+n=y+n反问题的提法反问题的提法反问题的提法反问题的提法：1.1.已知已知y y，分离信号和噪音，求出，分离信号和噪音，求出y y；2.2.已知已知y y和和A A，求出输入数据，求出输入数据x x；3.3.已知已知y y和和x x，求出模型机制，求出模型机制A A；4.4.已知已知y y，求出输入，求出输入x x和模型和模型A A，这时问题变成非线性，这时问题变成非线性。大家好19反问题通常体现了一种逆向思维。冯康先反问题通常体现了一种逆向思维。冯康先生在上世纪八十年代初曾经著文生在上世纪八十年代初曾经著文数学物理数学物理中的反问题中的反问题，较早地介绍了这个新的研究，较早地介绍了这个新的研究方向。他将反问题的功能概括为方向。他将反问题的功能概括为“由表及里由表及里”、“索隐探秘索隐探秘”、“倒果求因倒果求因”。在中国。在中国的传统文化中，只有智者高人才能透过现象的传统文化中，只有智者高人才能透过现象看清本质，甚至参透因果，一语破的。科学看清本质，甚至参透因果，一语破的。科学化的反问题研究为我们在解决问题，增长智化的反问题研究为我们在解决问题，增长智慧方面提供了很好的案例和方法论。慧方面提供了很好的案例和方法论。大家好20反问题发展简史反问题发展简史数学物理反问题最早出现在地球物理领域数学物理反问题最早出现在地球物理领域18461846年，法国人年，法国人Le VerrierLe Verrier发现了海王星发现了海王星18801880年，美国学者年，美国学者J.A.EwingJ.A.Ewing等人发明了近代地震仪，提等人发明了近代地震仪，提出了地震记录的分析问题出了地震记录的分析问题19071907年，年，HerglogHerglog提出了地震走时数据的反演提出了地震走时数据的反演19091909年，年，A.MohorovicicA.Mohorovicic发现了莫霍面发现了莫霍面19121912年，年，Beno GutenbegBeno Gutenbeg发现了古登堡面发现了古登堡面19351935年，年，LehmannLehmann发现了地球外核和内核的分界面发现了地球外核和内核的分界面这些工作都对地球物理反演学术思想的形成和发展起到了巨这些工作都对地球物理反演学术思想的形成和发展起到了巨大的推动作用。大的推动作用。大家好21但在第一台数字计算机诞生之前，反问题的发展非常缓慢，但在第一台数字计算机诞生之前，反问题的发展非常缓慢，反演方法只有选择法和量版法反演方法只有选择法和量版法1967-19701967-1970年，美国地球物理学家年，美国地球物理学家BackusBackus和应用数学家和应用数学家GilbertGilbert连续发表了三篇关于平均核法的文章，奠定了反连续发表了三篇关于平均核法的文章，奠定了反演理论的基础演理论的基础TikhonovTikhonov（吉洪诺夫）（吉洪诺夫）2020世纪世纪4040年代，提出了正则化方法年代，提出了正则化方法（19771977，Solutions of ill-posed problems,Solutions of ill-posed problems,美国。中译美国。中译本，本，19791979）7070年代初，英国学者年代初，英国学者G.HonsfieldG.Honsfield研制出了第一台医用研制出了第一台医用机以及他和美国学者机以及他和美国学者A.M.CormackA.M.Cormack共同获得了共同获得了19791979年度年度生理性和医学诺贝尔奖，大大推动了有关不可见物体层析生理性和医学诺贝尔奖，大大推动了有关不可见物体层析成像的研究热潮，也极大地推动了反问题数学理论、数值成像的研究热潮，也极大地推动了反问题数学理论、数值方法以及应用的发展方法以及应用的发展大家好22现代反问题开始于现代反问题开始于7070年代末、年代末、8080年代初，蓬勃发展至今，年代初，蓬勃发展至今，4040多年的时间。多年的时间。广泛应用于石油勘探、工程物探、无损探伤、航空航天、广泛应用于石油勘探、工程物探、无损探伤、航空航天、地下找水、光学、电子、控制等等领域，可以说无处不在。地下找水、光学、电子、控制等等领域，可以说无处不在。国际上四种反问题杂志：国际上四种反问题杂志：Inverse Problems,Inverse Inverse Problems,Inverse Problems and Imaging,Inverse Problems in Science Problems and Imaging,Inverse Problems in Science and Engineering,Journal of inverse and ill-posed and Engineering,Journal of inverse and ill-posed problemsproblems我国反问题的研究最早由计算数学家冯康倡导（我国反问题的研究最早由计算数学家冯康倡导（19821982）。）。他把反问题列为计算数学四大问题之一（正问题、反问题、他把反问题列为计算数学四大问题之一（正问题、反问题、逼近问题和代数问题）逼近问题和代数问题）大家好23科学史上的著名的案例科学史上的著名的案例1781年，天王星被确认为太阳系的第年，天王星被确认为太阳系的第7颗大行星。颗大行星。40年后，年后，法国天文学家法国天文学家Bouvard搜集了一个多世纪来的全部观测资料，搜集了一个多世纪来的全部观测资料，包括了包括了1781年之前的旧数据和之后的新数据，试图用牛顿的天年之前的旧数据和之后的新数据，试图用牛顿的天体力学原理来计算天王星的运动轨道。他发现了一个奇怪的现体力学原理来计算天王星的运动轨道。他发现了一个奇怪的现象：用全部数据计算出的轨道与旧数据吻合得很好，但是与新象：用全部数据计算出的轨道与旧数据吻合得很好，但是与新数据相比误差远超出精度允许的范围；如果仅以新数据为依据数据相比误差远超出精度允许的范围；如果仅以新数据为依据重新计算轨道，得到的结果又无法和旧数据相匹配。重新计算轨道，得到的结果又无法和旧数据相匹配。Bouvard的治学态度非常严谨，他在论文中指出：的治学态度非常严谨，他在论文中指出：“两套数据的不符究两套数据的不符究竟是因为旧的观测记录不可靠，还是来自某个外部未知因素对竟是因为旧的观测记录不可靠，还是来自某个外部未知因素对这颗行星的干扰？我将这个谜留待将来去揭示。这颗行星的干扰？我将这个谜留待将来去揭示。”大家好24首先，首先，Bouvard等天文学家核查了等天文学家核查了1750年以后英国格林尼年以后英国格林尼治天文台对各个行星所作的全部观测记录。结果发现，除天王治天文台对各个行星所作的全部观测记录。结果发现，除天王星以外，对于其它行星的观测记录与理论计算结果都符合得相星以外，对于其它行星的观测记录与理论计算结果都符合得相当好。似乎没有理由怀疑旧的天文观测唯独对天王星失准。既当好。似乎没有理由怀疑旧的天文观测唯独对天王星失准。既然如此，天文学家就需要对天王星的不规律运动作出科学的解然如此，天文学家就需要对天王星的不规律运动作出科学的解释。释。摆在天文学家面前的有两条路：摆在天文学家面前的有两条路：摆在天文学家面前的有两条路：摆在天文学家面前的有两条路：第一条路是质疑牛顿力学的普适性，或许万有引力定律不第一条路是质疑牛顿力学的普适性，或许万有引力定律不适用于距离太阳遥远的天王星，需要对之进行修正；适用于距离太阳遥远的天王星，需要对之进行修正；第二条路是寻找第二条路是寻找Bouvard所猜测的所猜测的“未知因素未知因素”。于是人们。于是人们提出了提出了“彗星撞击彗星撞击”、“未知卫星未知卫星”和和“未知行星未知行星”等多种可等多种可能。能。大家好251841年的暑期，还是英国剑桥大学二年级学生的年的暑期，还是英国剑桥大学二年级学生的Adams就就定下计划，不仅要确认天王星的轨道异常是否来自未知行星的定下计划，不仅要确认天王星的轨道异常是否来自未知行星的引力作用，还要尽可能地确认这颗新行星的轨道，以便通过观引力作用，还要尽可能地确认这颗新行星的轨道，以便通过观测来发现之。这不仅是一个新问题，而且是一个反问题。因为测来发现之。这不仅是一个新问题，而且是一个反问题。因为过去总是已知一颗行星的质量和轨道，根据万有引力定律计算过去总是已知一颗行星的质量和轨道，根据万有引力定律计算出它对另一颗行星产生的轨道摄动。而现在则相反，出它对另一颗行星产生的轨道摄动。而现在则相反，Adams要要假定已知天王星轨道的摄动，来计算出产生这一摄动的未知行假定已知天王星轨道的摄动，来计算出产生这一摄动的未知行星的质量和轨道。由于未知因素很多，实际计算起来是相当复星的质量和轨道。由于未知因素很多，实际计算起来是相当复杂和困难的。杂和困难的。Adams于于1845年彻底解决了这个反问题。他所运用的方法年彻底解决了这个反问题。他所运用的方法在当时是空前新颖的。令人遗憾的是，英国天文学家在当时是空前新颖的。令人遗憾的是，英国天文学家Airy先入先入为主地认为天王星的轨道问题是引力定律不再适用的结果，没为主地认为天王星的轨道问题是引力定律不再适用的结果，没有重视有重视Adams向他提交的新行星的轨道计算结果。向他提交的新行星的轨道计算结果。大家好26几乎与此同时，法国人几乎与此同时，法国人LeVerrier独立地解决了同样独立地解决了同样的反问题。的反问题。1846年年9月月23日，柏林天文台的日，柏林天文台的Galle按照按照LeVerrier提交的计算轨道着手观测，当晚就在偏离预言位提交的计算轨道着手观测，当晚就在偏离预言位置不到置不到1度的地方发现了一颗新的八等星。连续观测的数度的地方发现了一颗新的八等星。连续观测的数据都与据都与LeVerrier的预测结果吻合得很好，证实这是一颗的预测结果吻合得很好，证实这是一颗新行星。这时英国天文台才想起了新行星。这时英国天文台才想起了Adams的工作，悔之晚的工作，悔之晚矣。矣。案子破了。干扰天王星正常运行的那颗神秘天体正是案子破了。干扰天王星正常运行的那颗神秘天体正是太阳系的第太阳系的第8颗大行星颗大行星海王星！不仅长期困扰天文界海王星！不仅长期困扰天文界的天王星轨道异常问题在牛顿力学框架内得到了完满解释，的天王星轨道异常问题在牛顿力学框架内得到了完满解释，而且海王星的发现进一步验证了牛顿力学的正确性。而且海王星的发现进一步验证了牛顿力学的正确性。大家好27大家好28 正问题正问题正问题正问题-“知易行难知易行难”反问题反问题反问题反问题的研究重在探测和发现的研究重在探测和发现-“知之非艰、行知之非艰、行之惟艰之惟艰”。反问题的求解往往违背了事物发展过程的自然顺反问题的求解往往违背了事物发展过程的自然顺序，从而使正问题中的许多良好性质不再满足。序，从而使正问题中的许多良好性质不再满足。更何况我们搜集到的资料经常是真伪交杂、缺失更何况我们搜集到的资料经常是真伪交杂、缺失含混的，这就更增加了求解的困难。尽管如此，含混的，这就更增加了求解的困难。尽管如此，解决反问题仍然像破案猜谜一样引人入胜。解决反问题仍然像破案猜谜一样引人入胜。大家好29什么是反问题？什么是反问题？任何事物都充满着辩证关系，都有正、反两个方面。任何事物都充满着辩证关系，都有正、反两个方面。例如：例如：正问题（正问题（directproblem）：求一个多项式的零点）：求一个多项式的零点反问题（反问题（inverseproblem）：给定零点）：给定零点求求次多项式次多项式正问题正问题（directproblem）：计算已知多项式）：计算已知多项式在在的值的值反问题反问题（inverseproblem）：已知多项式）：已知多项式在点在点的值的值（或近似值），（或近似值），求求（Lagrange插值问题）。插值问题）。大家好30。正问题正问题（directproblem）：已知）：已知，求，求反问题反问题（inverseproblem）：已知）：已知0,1上的连续可微函数上的连续可微函数到底哪一个是反问题？到底哪一个是反问题？先前提出的、人们熟悉的是正问题，而另一个是反问题；先前提出的、人们熟悉的是正问题，而另一个是反问题；稳定的问题是正问题，不稳定的是反问题。稳定的问题是正问题，不稳定的是反问题。大家好31按系统科学观点按系统科学观点l正问题：由因求果正问题：由因求果l反问题：由果求因，或由因果求关系反问题：由果求因，或由因果求关系数学物理反问题：从工程物理中提出的各类反问题，数学物理反问题：从工程物理中提出的各类反问题，都有明确的数学形式和定量要求，都有明确的数学形式和定量要求，故称之为数学物理反问题。故称之为数学物理反问题。微分方程反问题：由微分方程描述的反问题微分方程反问题：由微分方程描述的反问题因讯号因讯号（输入）（输入）关系（模型）关系（模型）果信号果信号（输出）（输出）大家好32大家好33例如例如k，需要提供一些可测量的信息，例如：，需要提供一些可测量的信息，例如：（物体的外界接触表面上的各点流量已知），称（物体的外界接触表面上的各点流量已知），称为附加条件，就构成了热传导方程反问题模型。为附加条件，就构成了热传导方程反问题模型。大家好34大家好35按问题分类：按问题分类：(1)Computerizedtomography(CT)(2)Inversescattering(3)Inverseheatconductionproblems(4)Geophysicalinverseproblems(5)Inverseproblemsinimaging(6)Identificationofparametersindifferentialequations大家好363.非线性与不适定性非线性与不适定性(1)反问题通常是非线性的，而且常常是反问题通常是非线性的，而且常常是highlynonlinear，即使相应的正，即使相应的正问题是线性的。问题是线性的。(2)反问题一般来说是不适定的（反问题一般来说是不适定的（ill-posed）。）。适定与不适定的概念，最初是从研究数学物理方程定解问题提出的。适定与不适定的概念，最初是从研究数学物理方程定解问题提出的。20世纪世纪20年代，法国数学家年代，法国数学家Hadamard（阿达玛）对数学物理方程定解问题（阿达玛）对数学物理方程定解问题提出了如下适定性定义，即定解问题是适定的，是指它具备如下三个条提出了如下适定性定义，即定解问题是适定的，是指它具备如下三个条件：件：问题的解存在；问题的解存在；问题的解唯一；问题的解唯一；解连续依赖于定解条件，即定解条件（数据）改变很小时，相应的解解连续依赖于定解条件，即定解条件（数据）改变很小时，相应的解也改变很小，也就是解具有稳定性。也改变很小，也就是解具有稳定性。不满足上述三个条件之一的问题称为不满足上述三个条件之一的问题称为不适定问题不适定问题。大家好37粗略地说，粗略地说，凡是解不连续地依赖数据的一切数凡是解不连续地依赖数据的一切数学问题都称为不适定的学问题都称为不适定的，不仅限于定解问题。，不仅限于定解问题。符合适定性的数学物理问题很多如：符合适定性的数学物理问题很多如：Laplace方程的方程的Dirichlet问题；问题；波动方程的空间波动方程的空间Cauchy问题；问题；热传导方程的初值问题。热传导方程的初值问题。大家好38大家好39大家好40大家好41大家好42而方程组而方程组的解是的解是大家好43大家好44大家好45取其解取其解 称称为原方程原方程组的正的正则化解，化解，称称为正正则化参数，它化参数，它在在稳定性定性与精度之间进行调节。与精度之间进行调节。大家好大家好46大家好47例 积分方程是唯一可解的，当，解用梯形公式近似积分得其中，从而得到线性系统的近似。大家好48下表表明在时，真解和近似解的误差。表 真解与近似解的误差从表中可以看出，近似解与真解没有关系。而且离散的程度越好，得到的结果越差。大家好49大家好50大家好51大家好52大家好53大家好54大家好55大家好56大家好57大家好58594.4.反反问题应用用举例例 尽管一些尽管一些经典反典反问题的研究可以追溯很早，反的研究可以追溯很早，反问题这一学科的一学科的兴起却是近几十年来的事情。起却是近几十年来的事情。在科学研究中在科学研究中经常要通常要通过间接接观测来探求位于不可达、不可触之来探求位于不可达、不可触之处的物的物质的的变化化规律；生律；生产中中经常要根据特定的功能常要根据特定的功能对产品品进行行设计，或按照某种目的，或按照某种目的对流程流程进行控制。行控制。这些都可以提出些都可以提出为某种形式的反某种形式的反问题。可。可见，反，反问题的的产生是科学研究不断深化和工程技生是科学研究不断深化和工程技术迅猛迅猛发展的展的结果，果，而而计算技算技术的革命又的革命又为它提供了重要的物它提供了重要的物质基基础。大家好60 现在，反问题的研究已经遍及现代化生产、生现在，反问题的研究已经遍及现代化生产、生活、科学研究的各个领域。简单的概括不足以说明活、科学研究的各个领域。简单的概括不足以说明问题，我们下面具体介绍一些常见的反问题类型，问题，我们下面具体介绍一些常见的反问题类型，希望大家能够对它有一个大概的了解。希望大家能够对它有一个大概的了解。大家好61物性探测物性探测给你一只管子，不允许直接进入内部测量，你能算出里面给你一只管子，不允许直接进入内部测量，你能算出里面的形状吗？如果管子是轴对称的，这时只需要知道内部的形状吗？如果管子是轴对称的，这时只需要知道内部的截面半径就可以了。美国贝尔电话实验室的的截面半径就可以了。美国贝尔电话实验室的Sondhi和和Gophinath提供了一个方法：在管子的一边发出声音，提供了一个方法：在管子的一边发出声音，用仪器测量管口的位移速度和压力。通过测量结果就可用仪器测量管口的位移速度和压力。通过测量结果就可以推知管内的截面半径。理论计算与实验结果吻合得很以推知管内的截面半径。理论计算与实验结果吻合得很好。好。不要小看了这个例子，它实际上暗示了许多不能直接测量不要小看了这个例子，它实际上暗示了许多不能直接测量的物性探测问题可以通过类似的间接方法来解决。我们的物性探测问题可以通过类似的间接方法来解决。我们通常说通常说“上天入地上天入地”都是很困难的事情，可是在一些情都是很困难的事情，可是在一些情况下似乎必须况下似乎必须“入地入地”才能解决问题。才能解决问题。大家好62比如比如说石油勘探。石油通常埋在几千米的地下，无法石油勘探。石油通常埋在几千米的地下，无法直接直接观察油田的位置和察油田的位置和储量，靠量，靠试打井的打井的办法来探法来探测不但不但费用昂用昂贵（一口井的代价要上千万元），而且效率极低（一口井的代价要上千万元），而且效率极低（只能探（只能探测到井附近的局部信息）。一个可行的到井附近的局部信息）。一个可行的办法是通法是通过地面爆炸向地下地面爆炸向地下发射地震波，同射地震波，同时接收地接收地层的反射波信的反射波信号。可以想象，地面接收到的反射信号中含有地下的物性号。可以想象，地面接收到的反射信号中含有地下的物性结构信息（地构信息（地层的密度、声速等等），利用数学手段将的密度、声速等等），利用数学手段将这些信息提取出来，就可以些信息提取出来，就可以对地下的油地下的油储及其分布作出科学及其分布作出科学的判断。的判断。这很象在夏天人很象在夏天人们挑西瓜，把瓜放在耳挑西瓜，把瓜放在耳边拍一拍，拍一拍，有有经验的人就知道瓜瓤熟不熟，不需要切开来看，不会破的人就知道瓜瓤熟不熟，不需要切开来看，不会破坏西瓜的完整。坏西瓜的完整。大家好63 类似的探测方法可以应用于许多方面，如：农类似的探测方法可以应用于许多方面，如：农用土壤分析、地下水勘查，甚至于在考古发现上也用土壤分析、地下水勘查，甚至于在考古发现上也有应用。位于三峡库区的四川省云阳县故陵镇有一有应用。位于三峡库区的四川省云阳县故陵镇有一个大土包，相传为楚国古墓，但是历经三千余年的个大土包，相传为楚国古墓，但是历经三千余年的变迁，已经难以确认了。科技工作者在地表利用地变迁，已经难以确认了。科技工作者在地表利用地震波法、高精度磁法、电场岩性探测和地化方法四震波法、高精度磁法、电场岩性探测和地化方法四种手段进行探测，不但确认了古墓的存在，而且得种手段进行探测，不但确认了古墓的存在，而且得到了关于古墓的埋藏深度、形状、大小甚至墓道的到了关于古墓的埋藏深度、形状、大小甚至墓道的准确信息，为抢救和保护文物做出了贡献。准确信息，为抢救和保护文物做出了贡献。大家好64扫描成像扫描成像本世纪初，本世纪初，Hebglotz和和Wiechebt应用应用Abel型反演方法解决了型反演方法解决了在一定对称条件下通过地震波的走时曲线来反推地层内部形在一定对称条件下通过地震波的走时曲线来反推地层内部形貌的方法。据此貌的方法。据此Mohobovic（1909年）发现了地壳与地幔之年）发现了地壳与地幔之间的断层。现在，利用地震波的接收信号通过成像来考察地间的断层。现在，利用地震波的接收信号通过成像来考察地层地貌形态已经成为地球物理勘探最为重要的手段。例如，层地貌形态已经成为地球物理勘探最为重要的手段。例如，通过走时成像，可以得到地震波在不同深度的传播速度；而通过走时成像，可以得到地震波在不同深度的传播速度；而在已知速度的前提下，利用声波方程或其单程波方程偏移成在已知速度的前提下，利用声波方程或其单程波方程偏移成像方法，又可以得到反射界面的位置和形状。像方法，又可以得到反射界面的位置和形状。大家好65 成像的另一个重要应用是医学上的计算机层析成像成像的另一个重要应用是医学上的计算机层析成像（CT），这是），这是X光射线自光射线自Roentgen发明（获发明（获1900年诺贝尔奖）年诺贝尔奖）以来在医疗诊断上的重大进展，其发明人以来在医疗诊断上的重大进展，其发明人Hounsfield和和Cormack因此获得了因此获得了1979年的诺贝尔医学奖。年的诺贝尔医学奖。CT技术是医学、技术是医学、电子技术、计算机技术和反演数学相结合的产物，它利用计算电子技术、计算机技术和反演数学相结合的产物，它利用计算机来对穿越人体的机来对穿越人体的X射线信号进行处理，来重建体内的结构信射线信号进行处理，来重建体内的结构信息，生成透视图像供医疗诊断参考，其核心算法的数学基础是息，生成透视图像供医疗诊断参考，其核心算法的数学基础是二维二维Radon变换。继之而起的是基于三维变换。继之而起的是基于三维Radon变换的核磁共变换的核磁共振成像，在诊断效果和无伤害性方面更为优越。事实上，类似振成像，在诊断效果和无伤害性方面更为优越。事实上，类似的方法也可以借助于声波、光波、电磁波在无损探伤、雷达侦的方法也可以借助于声波、光波、电磁波在无损探伤、雷达侦察、射电望远镜探测、环境监测等多方面有广泛应用。察、射电望远镜探测、环境监测等多方面有广泛应用。大家好66逆时反演逆时反演 在科学研究中，我们经常遇到这样的问题：知道了某个事物在科学研究中，我们经常遇到这样的问题：知道了某个事物的现在状态，希望了解它的过去，即通常所说的的现在状态，希望了解它的过去，即通常所说的“恢复历史恢复历史的本来面目。的本来面目。”这往往可以提为逆时反问题。当然，反问题这往往可以提为逆时反问题。当然，反问题研究不是历史学，它所研究的对象一般要满足某种类型的演研究不是历史学，它所研究的对象一般要满足某种类型的演化方程或数学模式。例如，通过远程测得的某次爆炸产生的化方程或数学模式。例如，通过远程测得的某次爆炸产生的辐射波，如何确定爆炸的位置和初始能量？这是波动方程的辐射波，如何确定爆炸的位置和初始能量？这是波动方程的逆时反问题；又如，根据近来的温度变化能否确定过去某个逆时反问题；又如，根据近来的温度变化能否确定过去某个时间的温度状态？这就成为热传导方程的逆时反问题。时间的温度状态？这就成为热传导方程的逆时反问题。大家好67孤立子的发现孤立子的发现一个著名的例子是反散射方法在孤立子一个著名的例子是反散射方法在孤立子发现中的作用：反散中的作用：反散射射问题是量子物理学研究中的一个是量子物理学研究中的一个问题，通，通过谱和和谱函数在函数在无无穷远处的散射性的散射性态反推一反推一维Schordinger（薛定鄂）方程的（薛定鄂）方程的位位势函数。它由前函数。它由前苏联数学家数学家Gelfand和和Levitan（1955年）年）一一举解决。在此基解决。在此基础上引上引发了一系列突破性了一系列突破性进展，最展，最为著名著名的是利用的是利用这个个结果果Lax（1968年）得到了关于年）得到了关于KDV方程的巧方程的巧妙解法，从而妙解法，从而发现了非了非线性方程中的孤立子性方程中的孤立子现象。象。这是近代是近代非非线性科学研究的重要事件。性科学研究的重要事件。大家好68其它问题其它问题飞行器外形设计飞行器外形设计公路、桥梁、水坝等工程质量检测问题公路、桥梁、水坝等工程质量检测问题沙漠找水沙漠找水钢炉检测钢炉检测大家好69 前面介前面介绍了反了反问题的几种的几种类型，它型，它们在研究和在研究和应用上用上经常是相互常是相互联系的，分系的，分门别类只是只是为了叙述了叙述方便。另外，反方便。另外，反问题与其它数学学科之与其它数学学科之间并没有一并没有一个个严格的界限，而是互格的界限，而是互为补充，互相促充，互相促进。反。反问题的研究起源于数理方程，其反演算法中包含了微分的研究起源于数理方程，其反演算法中包含了微分方程数方程数值解法、最解法、最优化方法和概率化方法和概率统计等方面的等方面的许多思想和技巧。另一方面，反多思想和技巧。另一方面，反问题的研究也促的研究也促进了了人人们对世界的世界的认识，使得研究更全面、深化。，使得研究更全面、深化。大家好70 5.5.反问题研究的难点及对策反问题研究的难点及对策 与正问题相比，反问题的研究起步较晚，发展还远不成熟。与正问题相比，反问题的研究起步较晚，发展还远不成熟。从本质上来说，反问题的研究的难度一般比相应的正问题要从本质上来说，反问题的研究的难度一般比相应的正问题要大。这是因为反问题的求解往往违背了物理过程的自然顺大。这是因为反问题的求解往往违背了物理过程的自然顺序，从而使正问题中的许多良好性质不再满足。这种现象在序，从而使正问题中的许多良好性质不再满足。这种现象在许多学科的研究中都是普遍存在的。比如说：曹雪芹创作了许多学科的研究中都是普遍存在的。比如说：曹雪芹创作了古典名著古典名著红楼梦红楼梦，这是人所共知的，但是要从现存的史，这是人所共知的，但是要从现存的史料和文物料和文物“碎片碎片”来恢复这位伟大作家的人生经历和创作历来恢复这位伟大作家的人生经历和创作历程则是一件万分艰辛的事情，更何况这些程则是一件万分艰辛的事情，更何况这些“碎片碎片”信息真伪信息真伪交