幾類反問題的有效數(shù)值解法
發(fā)布時間:2021-12-22 22:52
反問題在許多領域有著非常廣泛的重要應用,例如醫(yī)學成像、現(xiàn)代工業(yè)、氣象預報等等.也因此,數(shù)值求解反問題引起了廣泛關注.在本論文中,我們考慮以下幾類反問題:帶箱型狀態(tài)約束的橢圓型最優(yōu)控制問題;帶積分型狀態(tài)梯度約束和箱型控制約束的橢圓型最優(yōu)控制問題;圖像恢復中的病態(tài)的線性反問題.受到一階算法成功用于求解大規(guī)模有限維優(yōu)化問題的啟發(fā),我們試圖研究如何利用一些有效的一階算法并結合這些問題的結構來更有效地求解它們,取得的主要成果如下:1.帶箱型狀態(tài)約束的橢圓型最優(yōu)控制問題.對于逐點狀態(tài)約束,相應的拉格朗日乘子是Borel測度函數(shù),這使得最優(yōu)性條件中的互補性條件無法表示成逐點的形式.Lavrentiev正則化被用來解決這一困難,然后我們用分片線性有限元全離散方法來離散正則化的問題.我們關于正則化和離散化帶來的誤差給出了誤差分析,由此可以看到全離散的誤差階并不比變分離散的誤差階差.為了數(shù)值求解離散問題,我們提出了一種異構交替方向乘子法(hADMM)以及一種兩階段策略,它們分別用來得到中等精度和高精度的解.2.帶積分型狀態(tài)梯度約束的橢圓型最優(yōu)控制問題.我們首先證明了該問題的最優(yōu)性條件并提出了一種有限元對偶...
【文章來源】:大連理工大學遼寧省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:126 頁
【學位級別】:博士
【部分圖文】:
圖2.1例2.1中A?=¥時的理想狀態(tài)%.??2" ̄
?幾類反問題的有效數(shù)值解法???取一個相對大的值即可,所以S只需要取一個很小的值就可以保證C(^<£.從這個意義??上說,sHSTA的收斂性受第一項支配,并且其收斂行為與FISTA相似.??作為一個例子,我們在圖4.1中畫出測試圖像‘hst’(有關這一圖像的更多細節(jié)可參見??第4.4部分中的例4.1)中矩陣P和P的奇異值?很容易看到這兩個矩陣的奇異值都很快地??下降到0?例如,Z3的第六大奇異值與最大奇異值之比僅為6.47e-2,P的第十大奇異值??與最大奇異值之比降低到4.27e?-?2.這些結構也可以在其他測試示例中觀察到.??0?!?'?'?' ̄,.?—I?5,?,?,?,?, ̄-?―,????l〇g(^i)[?。乂?1?°??-10?■?-5?-?-??<D?15?■?〇>?10?■????1?20-?5?-15?-?.??^?25?-?-?3?.20?-?_??E?W??W?30-???-.S?-?.??-35?-?-30?-?■??‘???\im?—^?-35?-??+??o??45??'?1???—■-?1?1??-40??1?1?1?1?1???〇?50?1?00?150?200?250?300?0?50?100?1?50?200?250?300??i?(index?of?singular?value)?i?(index?of?singular?value)??(a)尸的奇異值?(b)戶的奇異值??圖4.1例4.1:測試圖像‘hst’中矩陣Z5和戶的奇
?幾類反問題的有效數(shù)值解法???真實圖像?模糊噪聲圖像??jj^J|i!?Q|;??FISTA得到的圖像?sFISTA得到的圖像??圖4.4例4.1中從函數(shù)PRblurdefocus提取的‘satellite’的圖像.??Fig.?4.4?Figures?for?‘satellite’?extracted?from?PRblurdefocus?for?Example?4.1??0.65??s?1?1?!?1?1??■■■!?-i?■■?_i?-i?0.3??1?1?1?1?!?1??i?i?i????—-sFISTA?(s=5)?一?■?sFISTA?(s=5)??〇6?-\?——sFISTA?s=1?■?——sFISTA?(s=1)??\?—FISTA?^〇25?l?—FISTA??E?0.45?■?-g?0-15?■?\?-??I..?V???I?\??i-?------------?rA^??03?■?——????0.25??:?1?1?1?1?1?1?1?1??0??1?1?1?1?1?1?1?1?1???0?5?10?15?20?25?30?35?40?45?50?0?5?10?15?20?25?30?35?40?45?50??iteration?number?iteration?number??相對誤差77?相對殘差y??圖4.5例4.1中從函數(shù)PRblurdefocus提取的‘satellite’的77和7關于迭代次數(shù)的變化趨勢.??Fig.?4.5?t]
【參考文獻】:
期刊論文
[1]An FE-Inexact Heterogeneous ADMM for Elliptic Optimal Control Problems with L1-Control Cost[J]. SONG Xiaoliang,YU Bo,WANG Yiyang,ZHANG Xuping. Journal of Systems Science & Complexity. 2018(06)
本文編號:3547257
【文章來源】:大連理工大學遼寧省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:126 頁
【學位級別】:博士
【部分圖文】:
圖2.1例2.1中A?=¥時的理想狀態(tài)%.??2" ̄
?幾類反問題的有效數(shù)值解法???取一個相對大的值即可,所以S只需要取一個很小的值就可以保證C(^<£.從這個意義??上說,sHSTA的收斂性受第一項支配,并且其收斂行為與FISTA相似.??作為一個例子,我們在圖4.1中畫出測試圖像‘hst’(有關這一圖像的更多細節(jié)可參見??第4.4部分中的例4.1)中矩陣P和P的奇異值?很容易看到這兩個矩陣的奇異值都很快地??下降到0?例如,Z3的第六大奇異值與最大奇異值之比僅為6.47e-2,P的第十大奇異值??與最大奇異值之比降低到4.27e?-?2.這些結構也可以在其他測試示例中觀察到.??0?!?'?'?' ̄,.?—I?5,?,?,?,?, ̄-?―,????l〇g(^i)[?。乂?1?°??-10?■?-5?-?-??<D?15?■?〇>?10?■????1?20-?5?-15?-?.??^?25?-?-?3?.20?-?_??E?W??W?30-???-.S?-?.??-35?-?-30?-?■??‘???\im?—^?-35?-??+??o??45??'?1???—■-?1?1??-40??1?1?1?1?1???〇?50?1?00?150?200?250?300?0?50?100?1?50?200?250?300??i?(index?of?singular?value)?i?(index?of?singular?value)??(a)尸的奇異值?(b)戶的奇異值??圖4.1例4.1:測試圖像‘hst’中矩陣Z5和戶的奇
?幾類反問題的有效數(shù)值解法???真實圖像?模糊噪聲圖像??jj^J|i!?Q|;??FISTA得到的圖像?sFISTA得到的圖像??圖4.4例4.1中從函數(shù)PRblurdefocus提取的‘satellite’的圖像.??Fig.?4.4?Figures?for?‘satellite’?extracted?from?PRblurdefocus?for?Example?4.1??0.65??s?1?1?!?1?1??■■■!?-i?■■?_i?-i?0.3??1?1?1?1?!?1??i?i?i????—-sFISTA?(s=5)?一?■?sFISTA?(s=5)??〇6?-\?——sFISTA?s=1?■?——sFISTA?(s=1)??\?—FISTA?^〇25?l?—FISTA??E?0.45?■?-g?0-15?■?\?-??I..?V???I?\??i-?------------?rA^??03?■?——????0.25??:?1?1?1?1?1?1?1?1??0??1?1?1?1?1?1?1?1?1???0?5?10?15?20?25?30?35?40?45?50?0?5?10?15?20?25?30?35?40?45?50??iteration?number?iteration?number??相對誤差77?相對殘差y??圖4.5例4.1中從函數(shù)PRblurdefocus提取的‘satellite’的77和7關于迭代次數(shù)的變化趨勢.??Fig.?4.5?t]
【參考文獻】:
期刊論文
[1]An FE-Inexact Heterogeneous ADMM for Elliptic Optimal Control Problems with L1-Control Cost[J]. SONG Xiaoliang,YU Bo,WANG Yiyang,ZHANG Xuping. Journal of Systems Science & Complexity. 2018(06)
本文編號:3547257
本文鏈接:http://sikaile.net/shoufeilunwen/jckxbs/3547257.html
最近更新
教材專著