1、 數(shù)學(xué)物理方程反問題數(shù)學(xué)物理方程反問題 麥宏晏麥宏晏2021-12-2822021-12-282n什么是反問題n應(yīng)用背景n具體實(shí)例2021-12-283什么是反問題n耳朵能耳朵能“聽出聽出”鼓的形狀嗎？鼓的形狀嗎？ 僅僅通過鼓的聲音能否判斷出鼓的形狀？問題最早由丹僅僅通過鼓的聲音能否判斷出鼓的形狀？問題最早由丹麥著名物理學(xué)家麥著名物理學(xué)家Lorentz在在1910年的一次講演中提出，它的年的一次講演中提出，它的背景來自于射線理論。背景來自于射線理論。 一個(gè)物體的音色可以由一串譜一個(gè)物體的音色可以由一串譜 來確定，它們?cè)趤泶_定，它們?cè)谖锢砩蠈?duì)應(yīng)著物體的固有頻率。物理上對(duì)應(yīng)著物體的固有頻率。“盲人

2、聽鼓盲人聽鼓”即是要求通過即是要求通過已知的譜來確定一個(gè)鼓面的形狀。經(jīng)過數(shù)學(xué)家們的演算給已知的譜來確定一個(gè)鼓面的形狀。經(jīng)過數(shù)學(xué)家們的演算給出了否定的答案。但是，從鼓聲中我們確實(shí)能得到相當(dāng)多出了否定的答案。但是，從鼓聲中我們確實(shí)能得到相當(dāng)多的形狀信息：我們能夠的形狀信息：我們能夠“聽聽”出鼓的面積有多大、周邊有出鼓的面積有多大、周邊有多長(zhǎng)甚至鼓的內(nèi)部是否有洞、有幾個(gè)洞。多長(zhǎng)甚至鼓的內(nèi)部是否有洞、有幾個(gè)洞。 2021-12-2842021-12-284n什么是反問題什么是反問題？ 反問題是相對(duì)于正問題而言的。以前面所舉的反問題是相對(duì)于正問題而言的。以前面所舉的“盲人聽鼓盲人聽鼓”反問題為例，它的正

3、問題就是要在已知鼓的形狀的條件下，研反問題為例，它的正問題就是要在已知鼓的形狀的條件下，研究其發(fā)聲規(guī)律，這在數(shù)學(xué)物理歷史上已經(jīng)研究在先，而且比較究其發(fā)聲規(guī)律，這在數(shù)學(xué)物理歷史上已經(jīng)研究在先，而且比較成熟。此時(shí)鼓的所有譜都能通過一套算法利用計(jì)算機(jī)算出來。成熟。此時(shí)鼓的所有譜都能通過一套算法利用計(jì)算機(jī)算出來。我們可以這樣理解：世間的事物或現(xiàn)象之間往往存在著一定的我們可以這樣理解：世間的事物或現(xiàn)象之間往往存在著一定的自然順序，如時(shí)間順序、空間順序、因果順序，等等。所謂正自然順序，如時(shí)間順序、空間順序、因果順序，等等。所謂正問題，一般是按著這種自然順序來研究事物的演化過程或分布問題，一般是按著這種自然

4、順序來研究事物的演化過程或分布形態(tài)，起著由因推果的作用。反問題則是根據(jù)事物的演化結(jié)形態(tài)，起著由因推果的作用。反問題則是根據(jù)事物的演化結(jié)果，由可觀測(cè)的現(xiàn)象來探求事物的內(nèi)部規(guī)律或所受的外部影果，由可觀測(cè)的現(xiàn)象來探求事物的內(nèi)部規(guī)律或所受的外部影響，由表及里，索隱探秘，起著倒果求因的作用。可以看出，響，由表及里，索隱探秘，起著倒果求因的作用。可以看出，正、反兩方面都是科學(xué)研究的重要內(nèi)容。正、反兩方面都是科學(xué)研究的重要內(nèi)容。 2021-12-2852021-12-285n什么是反問題？什么是反問題？ 盡管一些經(jīng)典反問題的研究可以追溯很早，反問題這一學(xué)科盡管一些經(jīng)典反問題的研究可以追溯很早，反問題這一學(xué)科

5、的興起卻是近幾十年來的事情。在科學(xué)研究中經(jīng)常要通過間接觀的興起卻是近幾十年來的事情。在科學(xué)研究中經(jīng)常要通過間接觀測(cè)來探求位于不可達(dá)、不可觸之處的物質(zhì)的變化規(guī)律；生產(chǎn)中經(jīng)測(cè)來探求位于不可達(dá)、不可觸之處的物質(zhì)的變化規(guī)律；生產(chǎn)中經(jīng)常要根據(jù)特定的功能對(duì)產(chǎn)品進(jìn)行設(shè)計(jì)，或按照某種目的對(duì)流程進(jìn)常要根據(jù)特定的功能對(duì)產(chǎn)品進(jìn)行設(shè)計(jì)，或按照某種目的對(duì)流程進(jìn)行控制。這些都可以提出為某種形式的反問題。可見，反問題的行控制。這些都可以提出為某種形式的反問題。可見，反問題的產(chǎn)生是科學(xué)研究不斷深化和工程技術(shù)迅猛發(fā)展的結(jié)果，而計(jì)算技產(chǎn)生是科學(xué)研究不斷深化和工程技術(shù)迅猛發(fā)展的結(jié)果，而計(jì)算技術(shù)的革命又為它提供了重要的物質(zhì)基礎(chǔ)。術(shù)的

6、革命又為它提供了重要的物質(zhì)基礎(chǔ)。 現(xiàn)在，反問題的研究已經(jīng)遍及現(xiàn)代化生產(chǎn)、生活、研究的各現(xiàn)在，反問題的研究已經(jīng)遍及現(xiàn)代化生產(chǎn)、生活、研究的各個(gè)領(lǐng)域。我們下面具體介紹一些常見的反問題應(yīng)用背景，希望大個(gè)領(lǐng)域。我們下面具體介紹一些常見的反問題應(yīng)用背景，希望大家能夠?qū)λ幸粋€(gè)概括的了解。家能夠?qū)λ幸粋€(gè)概括的了解。2021-12-286定向設(shè)計(jì)物性探測(cè)掃描成像逆時(shí)反演及其他應(yīng)用背景2021-12-2872021-12-287 工業(yè)生產(chǎn)離不開產(chǎn)品設(shè)計(jì)，如何設(shè)計(jì)出優(yōu)質(zhì)產(chǎn)品使之更好工業(yè)生產(chǎn)離不開產(chǎn)品設(shè)計(jì)，如何設(shè)計(jì)出優(yōu)質(zhì)產(chǎn)品使之更好地實(shí)現(xiàn)其功能，是關(guān)系到廠家信譽(yù)和企業(yè)生存的大問題。在這地實(shí)現(xiàn)其功能，是關(guān)系到廠家

7、信譽(yù)和企業(yè)生存的大問題。在這方面，反問題研究可以為企業(yè)家出謀劃策。事實(shí)上，最早的反方面，反問題研究可以為企業(yè)家出謀劃策。事實(shí)上，最早的反問題研究就是起源于定向設(shè)計(jì)問題。我們知道，單擺的等時(shí)性問題研究就是起源于定向設(shè)計(jì)問題。我們知道，單擺的等時(shí)性只是在小角度的假設(shè)下才近似成立。能不能找到一種特殊軌線只是在小角度的假設(shè)下才近似成立。能不能找到一種特殊軌線的擺，使它嚴(yán)格滿足等時(shí)性？的擺，使它嚴(yán)格滿足等時(shí)性？Huygens于于1673年提出并解決了年提出并解決了這一問題，這種特殊的軌線就是旋輪線，它的方程為這一問題，這種特殊的軌線就是旋輪線，它的方程為 當(dāng)代工業(yè)產(chǎn)品的極大豐富為反問題的研究提供了廣闊的

8、用當(dāng)代工業(yè)產(chǎn)品的極大豐富為反問題的研究提供了廣闊的用武之地，許多工業(yè)設(shè)計(jì)問題是相當(dāng)困難的，需要用到高深的數(shù)武之地，許多工業(yè)設(shè)計(jì)問題是相當(dāng)困難的，需要用到高深的數(shù)學(xué)手段。例如，國(guó)外的光學(xué)儀器廠家提出：能否設(shè)計(jì)一種光柵學(xué)手段。例如，國(guó)外的光學(xué)儀器廠家提出：能否設(shè)計(jì)一種光柵，利用其非線性衍射效應(yīng)產(chǎn)生出高能量的單色射線？這就是一，利用其非線性衍射效應(yīng)產(chǎn)生出高能量的單色射線？這就是一 2021-12-288個(gè)定向設(shè)計(jì)問題，它要求數(shù)學(xué)家利用推導(dǎo)和計(jì)算手段構(gòu)造出所個(gè)定向設(shè)計(jì)問題，它要求數(shù)學(xué)家利用推導(dǎo)和計(jì)算手段構(gòu)造出所需要的曲面（光柵）形狀。需要的曲面（光柵）形狀。 定向設(shè)計(jì)的應(yīng)用相當(dāng)廣泛。比如說：一個(gè)城市的

9、某條街道定向設(shè)計(jì)的應(yīng)用相當(dāng)廣泛。比如說：一個(gè)城市的某條街道車流量很大，不堪負(fù)荷，怎樣通過鋪設(shè)新的路段來進(jìn)行分流？車流量很大，不堪負(fù)荷，怎樣通過鋪設(shè)新的路段來進(jìn)行分流？在軍事行動(dòng)中如何對(duì)不同種類的炮火進(jìn)行分布以達(dá)到特定的轟在軍事行動(dòng)中如何對(duì)不同種類的炮火進(jìn)行分布以達(dá)到特定的轟炸效果？這類問題往往涉及各種事物的組合、分配布局，要求炸效果？這類問題往往涉及各種事物的組合、分配布局，要求在各種相互制約、相互影響的因素中尋找出最佳方案，為領(lǐng)導(dǎo)在各種相互制約、相互影響的因素中尋找出最佳方案，為領(lǐng)導(dǎo)的決策提供依據(jù)。的決策提供依據(jù)。2021-12-2892021-12-289 給你一只管子，不允許直接進(jìn)入內(nèi)部

10、測(cè)量，你能算出里給你一只管子，不允許直接進(jìn)入內(nèi)部測(cè)量，你能算出里面的形狀嗎？如果管子是軸對(duì)稱的，這時(shí)只需要知道內(nèi)部的面的形狀嗎？如果管子是軸對(duì)稱的，這時(shí)只需要知道內(nèi)部的截面半徑就可以了。美國(guó)貝爾電話實(shí)驗(yàn)室的截面半徑就可以了。美國(guó)貝爾電話實(shí)驗(yàn)室的Sondhi和和Gophinath提供了一個(gè)方法：在管子的一邊發(fā)出聲音，用儀提供了一個(gè)方法：在管子的一邊發(fā)出聲音，用儀器測(cè)量管口的位移速度和壓力。通過測(cè)量結(jié)果就可以推知管器測(cè)量管口的位移速度和壓力。通過測(cè)量結(jié)果就可以推知管內(nèi)的截面半徑。內(nèi)的截面半徑。 這個(gè)例子，它實(shí)際上暗示了許多不能直接測(cè)量的物性探這個(gè)例子，它實(shí)際上暗示了許多不能直接測(cè)量的物性探測(cè)問題可

11、以通過類似的間接方法來解決。我們通常說測(cè)問題可以通過類似的間接方法來解決。我們通常說“上天上天入地入地”都是很困難的事情，可是在一些情況下似乎必須都是很困難的事情，可是在一些情況下似乎必須“入地入地”才能解決問題，比如說石油勘探。石油通常埋在幾千米的才能解決問題，比如說石油勘探。石油通常埋在幾千米的地下，無法直接觀察油田的位置和儲(chǔ)量，靠試打井的辦法來地下，無法直接觀察油田的位置和儲(chǔ)量，靠試打井的辦法來探測(cè)不但費(fèi)用昂貴（一口井的代價(jià)要上千萬元），而且效率探測(cè)不但費(fèi)用昂貴（一口井的代價(jià)要上千萬元），而且效率2021-12-2810極低（只能探測(cè)到井附近的局部信息）。一個(gè)可行的辦法是通極低（只能探測(cè)

12、到井附近的局部信息）。一個(gè)可行的辦法是通過地面爆炸向地下發(fā)射地震波，同時(shí)接收地層的反射波信號(hào)。過地面爆炸向地下發(fā)射地震波，同時(shí)接收地層的反射波信號(hào)。可以想象，地面接收到的反射信號(hào)中含有地下的物性結(jié)構(gòu)信息可以想象，地面接收到的反射信號(hào)中含有地下的物性結(jié)構(gòu)信息（地層的密度、聲速等等），利用數(shù)學(xué)手段將這些信息提取出（地層的密度、聲速等等），利用數(shù)學(xué)手段將這些信息提取出來，就可以對(duì)地下的油儲(chǔ)及其分布作出科學(xué)的判斷。來，就可以對(duì)地下的油儲(chǔ)及其分布作出科學(xué)的判斷。 類似的探測(cè)方法可以應(yīng)用于許多方面，如：農(nóng)用土壤分析類似的探測(cè)方法可以應(yīng)用于許多方面，如：農(nóng)用土壤分析、地下水勘查，甚至于在考古發(fā)現(xiàn)上也有應(yīng)用。

13、位于三峽庫(kù)區(qū)、地下水勘查，甚至于在考古發(fā)現(xiàn)上也有應(yīng)用。位于三峽庫(kù)區(qū)的四川省云陽縣故陵鎮(zhèn)有一個(gè)大土包，相傳為楚國(guó)古墓，但是的四川省云陽縣故陵鎮(zhèn)有一個(gè)大土包，相傳為楚國(guó)古墓，但是歷經(jīng)三千余年的變遷，已經(jīng)難以確認(rèn)了。科技工作者在地表利歷經(jīng)三千余年的變遷，已經(jīng)難以確認(rèn)了。科技工作者在地表利用地震波法、高精度磁法、電場(chǎng)巖性探測(cè)和地化方法四種手段用地震波法、高精度磁法、電場(chǎng)巖性探測(cè)和地化方法四種手段進(jìn)行探測(cè)，不但確認(rèn)了古墓的存在，而且得到了關(guān)于古墓的埋進(jìn)行探測(cè)，不但確認(rèn)了古墓的存在，而且得到了關(guān)于古墓的埋藏深度、形狀、大小甚至墓道的準(zhǔn)確信息，為搶救和保護(hù)文物藏深度、形狀、大小甚至墓道的準(zhǔn)確信息，為搶救和保

14、護(hù)文物作出了貢獻(xiàn)。作出了貢獻(xiàn)。2021-12-28112021-12-2811 如果把下落的物體用掃描射線替代，從另一個(gè)角度來看它如果把下落的物體用掃描射線替代，從另一個(gè)角度來看它為我們提供了從射線的走時(shí)響應(yīng)反推其傳播軌跡的方法，將不為我們提供了從射線的走時(shí)響應(yīng)反推其傳播軌跡的方法，將不同軌跡射線的反演結(jié)果組合起來就能得到傳播介質(zhì)的內(nèi)部形態(tài)同軌跡射線的反演結(jié)果組合起來就能得到傳播介質(zhì)的內(nèi)部形態(tài)信息。本世紀(jì)初，信息。本世紀(jì)初，Hebglotz和和Wiechebt應(yīng)用應(yīng)用Abel型反演方法型反演方法解決了在一定對(duì)稱條件下通過地震波的走時(shí)曲線來反推地層內(nèi)解決了在一定對(duì)稱條件下通過地震波的走時(shí)曲線來反

15、推地層內(nèi)部形貌的方法。據(jù)此部形貌的方法。據(jù)此Mohobovic（1909年）發(fā)現(xiàn)了地殼與地年）發(fā)現(xiàn)了地殼與地幔之間的斷層。現(xiàn)在，利用地震波的接收信號(hào)通過成像來考察幔之間的斷層。現(xiàn)在，利用地震波的接收信號(hào)通過成像來考察地層地貌形態(tài)已經(jīng)成為地球物理勘探最為重要的手段。例如，地層地貌形態(tài)已經(jīng)成為地球物理勘探最為重要的手段。例如，通過走時(shí)成像，可以得到地震波在不同深度的傳播速度；而在通過走時(shí)成像，可以得到地震波在不同深度的傳播速度；而在已知速度的前提下，利用聲波方程或其單程波方程偏移成像方已知速度的前提下，利用聲波方程或其單程波方程偏移成像方法，又可以得到反射界面的位置和形狀。法，又可以得到反射界面的

16、位置和形狀。 2021-12-2812 成像的另一個(gè)重要應(yīng)用是醫(yī)學(xué)上的計(jì)算機(jī)層析成像（成像的另一個(gè)重要應(yīng)用是醫(yī)學(xué)上的計(jì)算機(jī)層析成像（CT），這是），這是X光射線自光射線自Roentgen發(fā)明（獲發(fā)明（獲1900年諾貝爾獎(jiǎng)）以來在醫(yī)療診斷上年諾貝爾獎(jiǎng)）以來在醫(yī)療診斷上的重大進(jìn)展，其發(fā)明人的重大進(jìn)展，其發(fā)明人Hounsfield和和Cormack因此獲得了因此獲得了1979年的年的諾貝爾醫(yī)學(xué)獎(jiǎng)。諾貝爾醫(yī)學(xué)獎(jiǎng)。CT技術(shù)是醫(yī)學(xué)、電子技術(shù)、計(jì)算機(jī)技術(shù)和反演數(shù)學(xué)技術(shù)是醫(yī)學(xué)、電子技術(shù)、計(jì)算機(jī)技術(shù)和反演數(shù)學(xué)相結(jié)合的產(chǎn)物，它利用計(jì)算機(jī)來對(duì)穿越人體的相結(jié)合的產(chǎn)物，它利用計(jì)算機(jī)來對(duì)穿越人體的X射線信號(hào)進(jìn)行處理，射線

17、信號(hào)進(jìn)行處理，來重建體內(nèi)的結(jié)構(gòu)信息，生成透視圖象供醫(yī)療診斷參考，其核心算來重建體內(nèi)的結(jié)構(gòu)信息，生成透視圖象供醫(yī)療診斷參考，其核心算法的數(shù)學(xué)基礎(chǔ)是二維法的數(shù)學(xué)基礎(chǔ)是二維Radon變換。繼之而起的是基于三維變換。繼之而起的是基于三維Radon變變換的核磁共振成像，在診斷效果和無傷害性方面更為優(yōu)越。事實(shí)上換的核磁共振成像，在診斷效果和無傷害性方面更為優(yōu)越。事實(shí)上，類似的方法也可以借助于聲波、光波、電磁波在無損探傷、雷達(dá)，類似的方法也可以借助于聲波、光波、電磁波在無損探傷、雷達(dá)偵察、射電望遠(yuǎn)鏡探測(cè)、環(huán)境監(jiān)測(cè)等多方面有廣泛應(yīng)用。偵察、射電望遠(yuǎn)鏡探測(cè)、環(huán)境監(jiān)測(cè)等多方面有廣泛應(yīng)用。2021-12-28132

18、021-12-2813 我們經(jīng)常遇到這樣的問題：知道了某個(gè)事物的現(xiàn)在狀我們經(jīng)常遇到這樣的問題：知道了某個(gè)事物的現(xiàn)在狀態(tài)，希望了解它的過去，即通常所說的態(tài)，希望了解它的過去，即通常所說的“恢復(fù)歷史的本來面恢復(fù)歷史的本來面目。這往往可以提為逆時(shí)反問題。它所研究的對(duì)象一般要滿目。這往往可以提為逆時(shí)反問題。它所研究的對(duì)象一般要滿足某種類型的演化方程或數(shù)學(xué)模式。例如，通過遠(yuǎn)程測(cè)得的足某種類型的演化方程或數(shù)學(xué)模式。例如，通過遠(yuǎn)程測(cè)得的某次爆炸產(chǎn)生的輻射波，如何確定爆炸的位置和初始能量？某次爆炸產(chǎn)生的輻射波，如何確定爆炸的位置和初始能量？這是波動(dòng)方程的逆時(shí)反問題；又如，根據(jù)近來的溫度變化能這是波動(dòng)方程的逆時(shí)

19、反問題；又如，根據(jù)近來的溫度變化能否確定過去某個(gè)時(shí)間的溫度狀態(tài)？這就成為熱傳導(dǎo)方程的逆否確定過去某個(gè)時(shí)間的溫度狀態(tài)？這就成為熱傳導(dǎo)方程的逆時(shí)反問題。時(shí)反問題。 2021-12-2814 反問題的研究起源于數(shù)理方程，反問題的研究也促進(jìn)了人們反問題的研究起源于數(shù)理方程，反問題的研究也促進(jìn)了人們對(duì)世界的識(shí)。一個(gè)著名的例子是反散射方法在孤立子發(fā)現(xiàn)中的作對(duì)世界的識(shí)。一個(gè)著名的例子是反散射方法在孤立子發(fā)現(xiàn)中的作用：反散射問題是量子物理學(xué)研究中的一個(gè)問題，通過譜和譜函用：反散射問題是量子物理學(xué)研究中的一個(gè)問題，通過譜和譜函數(shù)在無窮遠(yuǎn)處的散射性態(tài)反推一維數(shù)在無窮遠(yuǎn)處的散射性態(tài)反推一維Schordinger方程

20、的位勢(shì)函數(shù)。方程的位勢(shì)函數(shù)。它由前蘇聯(lián)數(shù)學(xué)家它由前蘇聯(lián)數(shù)學(xué)家Gelfand和和Levitan（1955年）一舉解決。在年）一舉解決。在此基礎(chǔ)上引發(fā)了一系列突破性進(jìn)展，最為著名的是利用這個(gè)結(jié)果此基礎(chǔ)上引發(fā)了一系列突破性進(jìn)展，最為著名的是利用這個(gè)結(jié)果Lax（1968年）得到了關(guān)于年）得到了關(guān)于KDV方程的巧妙解法，從而發(fā)現(xiàn)了方程的巧妙解法，從而發(fā)現(xiàn)了非線性方程中的孤立子現(xiàn)象。這是近代非線性科學(xué)研究的重要事非線性方程中的孤立子現(xiàn)象。這是近代非線性科學(xué)研究的重要事件。件。2021-12-2815具體實(shí)例例例1 金融學(xué)中的一個(gè)簡(jiǎn)單計(jì)息反問題金融學(xué)中的一個(gè)簡(jiǎn)單計(jì)息反問題 在通常利率的連續(xù)復(fù)利率模型中，投資

21、變化的百分率是一個(gè)給定的常利率，用微分方程表示為其中 是 時(shí)刻投資價(jià)值，由此導(dǎo)出價(jià)值的指數(shù)增長(zhǎng)模型 在變化利率模型中，利率與時(shí)間有關(guān)，即 ，于是有 /dudlnuurdtdt( )u tt( ) =(0)rtu tue( )r = r t2021-12-2816 ，給定與時(shí)間有關(guān)的利率和初值 ，求價(jià)值記錄 的問題稱為正問題；另一方面，由價(jià)值記錄求變化率的問題稱為反問題。注意這里正問題的求解是積分運(yùn)算，這是一個(gè)穩(wěn)定的過程，而反問題的求解需要微分運(yùn)算，是不穩(wěn)定的。我們通常感興趣的是反問題，如果已知價(jià)值記錄的解析表達(dá)式，求 是一個(gè)簡(jiǎn)單的微分運(yùn)算，對(duì)時(shí)間進(jìn)行離散 ，用差商代替導(dǎo)數(shù)，則用來近似 ，其中

22、。 直接從原始模型出發(fā)可以建立求反問題近似解的另一種d0( )= (0)tr(t) tu tue(0)u()ut()r=r t012,ktkh,k= , ,1k+kklnulnurh( )kr t= ( )kkuu t2021-12-28172021-12-2817方法將方程兩邊在 上積分得到 ，利用梯形公式近似右端積分即得由此解出這里的問題中都假定利率和價(jià)值記錄是嚴(yán)格正的函數(shù)。durudt1kk+tt，11()( )=k+ktk+ktu tu trudt111= ()/2k+kkkk+k+uuh ru +r u1112()/=k+kkkk+k+uuhruru2021-12-2818例例2 熱

23、傳導(dǎo)問題熱傳導(dǎo)問題 與外界存在熱交換的物體表面溫度不能及時(shí)達(dá)到熱平衡，如果物體比外界熱，熱量從物體往外流出從而使物體冷卻，這一現(xiàn)象最簡(jiǎn)單的模型是牛頓冷卻定律，即表面溫度的變化與物體表面和外界溫度之差成比例，如果時(shí)刻 物體的表面溫度為 ，而外界溫度為常數(shù) ，由牛頓冷卻定律得其中熱交換系數(shù)為正常數(shù) ，這是一個(gè)經(jīng)典的指數(shù)衰減模型，確定表面溫度的正問題有唯一解t( )u tA()duAudt( ) =+ ( (0)tu tAuA e2021-12-28192021-12-2819這個(gè)解依賴外界溫度 、初始溫度 和熱交換系數(shù) 三個(gè)參數(shù)。當(dāng)然在適當(dāng)時(shí)候表面溫度觀測(cè)值確定了識(shí)別參數(shù) 、 和 的反問題的解。

24、熱是一種形式的能量，物體的熱含量是物體分子動(dòng)能的一種測(cè)量標(biāo)準(zhǔn)，溫度計(jì)測(cè)量得的溫度與物體分子的平均動(dòng)能有關(guān)。物體的熱含量不僅依賴于它的溫度，而且依賴于它的質(zhì)量，如3千克的鐵球在給定的溫度下，它的熱能是相同溫度的1千克的鐵球的3倍。熱還與材料的種類有關(guān)，1千克的棉花球在給定溫度下，其熱能比相同溫度的1千克鐵球的少。這一思想可描述為：其中 是物體的溫度， 是它的質(zhì)量， 是物質(zhì)的比熱。A(0)uA(0)uQ = cmuumc2021-12-28202021-12-2820 我們討論具有最簡(jiǎn)單幾何形狀的物體的內(nèi)部溫度：一個(gè)橫截面積為單位面積的單位長(zhǎng)度的棒子，并假設(shè)它正好處于 軸的單位坐標(biāo)上。假設(shè)棒子的質(zhì)

25、量密度和比熱分別為 和 ，并假設(shè)棒子的側(cè)面是隔熱的，所以溫度空間上僅為單變量 的函數(shù)。這樣的棒子在 上的點(diǎn)在 時(shí)刻的溫度函數(shù)為 。考慮棒子上的一小段區(qū)域 ，這一薄層上的熱量變化速度近似為 ，由能量守恒原理，這個(gè)量等于流入這一薄層的熱量速度加上這一薄層單位時(shí)間內(nèi)可能產(chǎn)生的熱量。如果定義（ 時(shí)刻 點(diǎn)）單位體積產(chǎn)生的熱量為 ，那么單位時(shí)間內(nèi)，橫截面積為單位面積的薄層產(chǎn)生的熱量近似為 。 熱只能從 的左側(cè)或 的右側(cè)流入（或流出）這一薄層。最后，由傅立葉定律，通過表面的熱流速度與表面c01x ，0t( , )u x tx,x+xc uxttx( , )f x tfxxx+ xxx2021-12-2821

26、上溫度梯度的負(fù)數(shù) 成比例。（這里的負(fù)號(hào)是因?yàn)闊崃渴菑臒崃飨蚶洌Ｒ虼藛挝粫r(shí)間內(nèi)通過表面流入薄層 的凈熱流加上單位時(shí)間內(nèi)部生成的熱量，我們得到這一薄層中熱量變化速度近似為由能量守恒原理，這與前面計(jì)算的比率相等，當(dāng)區(qū)間 收縮到點(diǎn) 時(shí)這一模型變成這就是著名的傅立葉熱方程。熱方程的正問題包括給定邊界條件，即端點(diǎn)的溫度 和初始溫度 分布及參數(shù)值 ，求解 時(shí)刻所有點(diǎn)的溫度，一般來說這些參數(shù)x+ xxuukkfxxxu/ x x0()x+ xxxuukkc uuxxlimf =kftxxx(0, ), (1, )ut u t( ,0)u xc,k0tx,x+ xx,x+ x2021-12-28222021-

27、12-2822是時(shí)間、空間和溫度的函數(shù)。 我們研究識(shí)別和估計(jì)熱模型中的分布參數(shù)，特別的，我們考慮由觀測(cè)得到的棒子中點(diǎn)的溫度變化過程，確定與時(shí)間有關(guān)的參數(shù) 的問題練習(xí)：練習(xí)：下面給出測(cè)量得到的不同時(shí)刻按牛頓定律冷卻的物體表面溫度，求解周圍溫度、初始溫度和熱交換系數(shù)。 ( ) = (0.5, )h tut( )a t22( ),01,0(0, )(1, )0( ,0) = sinuua txttxututu xx 5 10 15 72 62 54tu2021-12-28232021-12-2823例例3 線性代數(shù)中的反問題線性代數(shù)中的反問題 線性代數(shù)主要關(guān)注的是線性方程組解的存在性、唯一性和穩(wěn)定性

28、，線性代數(shù)的正問題指給定 矩陣和一個(gè) 維向量 ，求解 維向量 ，而尋求滿足 的解 的原因反問題受到更多的關(guān)注，還有一類識(shí)別反問題即確定矩陣使得對(duì)給定的輸入-輸出滿足 。 首先考慮原因反問題，求解向量 滿足 ，這個(gè)問題的解存在當(dāng)且僅當(dāng) 屬于子空間 ，這一子空間就是所有 的列向量的線性方程組合形成的 的子空間。確定 是否屬于 也就是解是否存在，尋找所有的解是通過高斯消去法來實(shí)現(xiàn)的。 m nnxmb=AxAx=bxAx=bnxRAx=bb( )=:nRAAx x RAmRb( )R A2021-12-28242021-12-2824唯一性問題是通過 的零空間 來解決的，高斯消去法是一個(gè)刻畫零空間的有

29、效手段。 的解 關(guān)于右端項(xiàng) 的擾動(dòng)的穩(wěn)定性，可以用矩陣的條件數(shù)來量化。我們假設(shè)矩陣 為可逆矩陣，我們想知道 在怎樣一個(gè)相對(duì)誤差內(nèi)會(huì)導(dǎo)致解的相對(duì)大的改變，假設(shè)是右端項(xiàng) 的一個(gè)擾動(dòng)，擾動(dòng)的大小可以用范數(shù)來衡量，令 是右端項(xiàng)為 的解，那么所以可以得到如下的相對(duì)誤差估計(jì)A( )=:0nN AxRAx =Ax=bxbAbbbbb / |b |x b111| | |xx = A bA bAb b11| | | | | |b bb bxxAbAAxbb2021-12-28252021-12-2825于是其中 稱為矩陣 的條件數(shù)，因此條件數(shù)給出了由給定的右端項(xiàng)相對(duì)誤差造成的解的相對(duì)誤差的上界。對(duì)于大條件數(shù)的矩

30、陣，即病態(tài)矩陣，右端相對(duì)小的擾動(dòng)會(huì)引起解的相對(duì)大的變化，在這種意義下，病態(tài)線性方程組是不穩(wěn)定的。 考慮無解或有無窮多解的線性方程組。當(dāng) 不屬于矩陣的值域時(shí)， 無解。零空間、值域和轉(zhuǎn)置之間的關(guān)系使得我們可以求它的一類廣義解：最小二乘解，即解向量 使得所有的 中的范數(shù)達(dá)到最小。如果 是最小二乘解，那么對(duì)于任意向量 ，函數(shù)|( )|xxbbcond Axb1( ) =|cond AAAbAnuRnxRnvRAx=bnxRA2021-12-28262021-12-2826在 時(shí)達(dá)到最小值，由極值的必要條件 得 ，所以對(duì)于所有的 成立 ，也就是說如果 是最小二乘解，那么 其中 是 的轉(zhuǎn)置矩陣。相反的，如

31、果 ，那么對(duì)于任何 ，2222( ) | ( + )|2(,)|g tA u tvbAu bAv Au b tAvt0t=(0) 0g(,)0Av AubnvR( ,)0v A AuA bA AuA bAAA AuA bnx R2222| | ()| ()|2(,) |Ax bA x uAu bA x ux u A AuA bAu bAu bnuR2021-12-28272021-12-2827即 是 最小二乘解。 所有的 最小二乘解就是對(duì)稱問題 的普通解，由于 ,這一對(duì)稱問題總是有一個(gè)解。事實(shí)上,如果 是 最小二乘解，那么對(duì)于任意 ， 也是最小二乘解，即最小二乘解的集合構(gòu)成平行于零空間的超平

32、面，因此如果 有一個(gè)非平凡的零空間，那么 有無限多個(gè)最小二乘解，其中有一個(gè)最小二乘解能夠區(qū)別其它的解，即這一解與零空間正交，這樣的最小二乘解至多只有一個(gè)。我們接受這一廣義解的概念，那么每個(gè)線性系統(tǒng)都有唯一的廣義解。 下面考慮識(shí)別反問題，即給定關(guān)于 的向量對(duì)確定矩陣 的反問題。通過詢問適當(dāng)?shù)?觀察輸出的 來識(shí)別黑箱，我們可以控制輸入，安排他們?yōu)榫€性無關(guān)的，用一個(gè) 矩陣表示輸入，其中矩陣的每一列分別輸入向量nu RAx=bA AuA b( )= ()R A R A Au( )vN Au vAx=bAx=bAx=bAx=bxbnpA1x ,AA2021-12-28282021-12-2828 ，記為

33、 ，相應(yīng)的輸出向量可以用 矩陣 來表示。如果存在唯一的矩陣 滿足 ，那么稱可以通過矩陣對(duì) 來識(shí)別 。考慮三種情況，首先 是不可能的，因?yàn)橄蛄拷M 是線性無關(guān)的；第二種情況 ，那么 是可逆的， 是可識(shí)別的。第三種情況 ，在這種情況下存在一個(gè)向量 與 正交，令 是 矩陣，它的第一行是 ，其他行都是零向量，那么 是 的零矩陣，因此如果 就有 ，此時(shí)從信息 無法識(shí)別 2,pxxXAAx=b( , )X BApn12,px ,xxp=nXApnnqR12,px ,xxCm nqCXm pAx=b()A+C X=B( , )X BAmpB2021-12-28292021-12-2829例例4 混合溶液的流動(dòng)

34、混合溶液的流動(dòng) 典型的正問題為給定一個(gè)初始濃度的溶液，以及溶液流入和流出的速度，計(jì)算將來時(shí)刻溶液的濃度。這里給出最簡(jiǎn)單的模型：已知容器的容積為 ，濃度為 的溶液以一個(gè)給定速度 流入，而這被充分?jǐn)嚭偷娜芤河忠韵嗤乃俣?從容器中排出。 這一模型的建立僅僅依賴于速度的平衡，如果 代表容器中時(shí)刻 溶質(zhì)的質(zhì)量，那么 隨時(shí)間變化的速度就是溶液流入容器的速度和它流出的速度之差，因此或者可以給出容器中溶液濃度的微分方程Var( )q ttqdqqarrdtVr2021-12-28302021-12-2830這一微分方程與通過牛頓得到的冷卻過程模型有相同的形式，有唯一的解其中 是容器中溶液的初始濃度，給定參數(shù)

35、 ，這一正問題的解給出了該溶液任意時(shí)刻的濃度。這一簡(jiǎn)單模型蘊(yùn)藏了許多有趣的反問題，例如，假設(shè)容器是個(gè)有污染的地下水滲入的地下蓄水池，通過抽樣這個(gè)蓄水池的溶液（通過預(yù)測(cè)的探測(cè)器），可以測(cè)量得到這一蓄水池任意時(shí)刻的液體濃度，這些測(cè)量結(jié)果可以用來反演這一模型中的參數(shù)，這一模型也可以推廣，例如流出和流入的速度不同時(shí)，這個(gè)模型中()d cracd tV0( ) =+()rtVc taca e0c0,a cV,r2021-12-28312021-12-2831體積 就是一個(gè)依賴于時(shí)間的參數(shù)。反問題也可以這樣提出，流速或流入液體的濃度是依賴于時(shí)間的。練習(xí)：練習(xí)：含有未知（常數(shù)）濃度的污染物質(zhì)的地下水以一個(gè)未

36、知（常數(shù)）速度滲入到一個(gè)容積為1000加侖的蓄水池中，這一被充分?jǐn)嚭偷娜芤阂酝瑯拥乃俣认蛲鉂B漏，測(cè)量結(jié)果表明蓄水池中污染物質(zhì)的初始濃度為 ，一天以后污染物濃度為 ，兩天后為 ，求地下水中污染物質(zhì)的濃度是多少？地下水滲入蓄水池的速度是多少？ V1%1.1%1.2%2021-12-28322021-12-2832例例5 利用萬有引力尋寶利用萬有引力尋寶 考慮質(zhì)量為 的質(zhì)點(diǎn)，位于單位深度的平靜的河底，假定任何其他的重力源都是完全均勻的，從而由質(zhì)點(diǎn)源產(chǎn)生的重力異常是唯一的重力影響源。確定河底的金塊對(duì)河面單位質(zhì)量質(zhì)點(diǎn)的引力是一個(gè)正問題，牛頓萬有引力定律認(rèn)為該作用力等于兩質(zhì)點(diǎn)質(zhì)量乘積的萬有引力常數(shù)倍，再除以兩質(zhì)點(diǎn)距離的平方。我們考慮基于河底測(cè)得數(shù)據(jù)來確定質(zhì)點(diǎn)質(zhì)量和位置的反問題，測(cè)得的數(shù)據(jù)包括河面基準(zhǔn)點(diǎn)到測(cè)量點(diǎn)的距離以及河底的金塊產(chǎn)生的對(duì)河面 處單位質(zhì)量質(zhì)點(diǎn)的引力的垂直分量的估計(jì)值 （例如利用靈敏彈簧秤測(cè)得）， mx2021-