解PDE約束優(yōu)化問題的交替方向迭代法
發(fā)布時(shí)間:2022-01-17 15:08
帶偏微分方程(PDE)約束優(yōu)化問題的數(shù)值求解是應(yīng)用數(shù)學(xué)領(lǐng)域中重要而具有挑戰(zhàn)性的問題之一,其在現(xiàn)代工業(yè)、醫(yī)學(xué)、經(jīng)濟(jì)學(xué)等應(yīng)用領(lǐng)域都具有很重要的應(yīng)用.對傳統(tǒng)的帶L2-控制成本的PDE約束優(yōu)化問題,理論和數(shù)值解法都取得了豐富的研究成果.而對于帶L1-控制成本的PDE約束優(yōu)化問題,研究成果還不多.與有限維l1-正則化一樣,L1-控制成本具有誘導(dǎo)稀疏的特性,因此該類問題在控制裝置的布放問題以及材料和機(jī)械裝置的拓?fù)鋬?yōu)化等領(lǐng)域有重要應(yīng)用.鑒于有限維稀疏優(yōu)化在交替方向類算法上和應(yīng)用上的成功,我們嘗試研究了帶L1-控制成本的PDE約束優(yōu)化問題的交替方向類算法,并取得如下研究成果:1.帶L1-控制成本的稀疏最優(yōu)控制問題的FE-ihADMM算法.若用常規(guī)的分片線性有限元對連續(xù)問題進(jìn)行離散,不同于有限維的l1-范數(shù),離散的L1-范數(shù)不具有可分結(jié)構(gòu),甚至不能通過引進(jìn)人工變量轉(zhuǎn)化成具有可分結(jié)構(gòu)的問題,因此不適合用有限維交替方向乘子法(ADMM)求解.我們提出了一種具有可分結(jié)構(gòu)的有限元離散格式,盡管該離散格式會帶來額外的離散誤差,但我們給出的誤差估計(jì)表明該離散格式仍具有與常規(guī)離散格式一樣的誤差階O(h).進(jìn)一步,我們...
【文章來源】:大連理工大學(xué)遼寧省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:217 頁
【學(xué)位級別】:博士
【部分圖文】:
圖2.3例2.2:最優(yōu)控制叫在/7?=?2_7的三維圖和平面圖.(a)三維圖.(b)平面圖.深紅色和深藍(lán)色區(qū)??域?qū)?yīng)w/z?=?±30,綠色區(qū)域?qū)?yīng)w;z?=?0??
帶控制約束的最優(yōu)控制問題的精確解是未知的.作為精確解的代替,我們采用網(wǎng)格大小??為/f?=?2_1()上的數(shù)值解作為一個(gè)參考解.??圖3.1給出了網(wǎng)格大小為A?=?2_6時(shí)離散最優(yōu)控制叫的三維圖.表3.1給出了在L2模意??義下控制《的離散誤差的數(shù)值結(jié)果.?dāng)?shù)值結(jié)果證實(shí)了定理3.]給出的理論誤差結(jié)果???表3.1中給出了ihADMM、經(jīng)典ADMM算法和LADMM算法求解例3.1的計(jì)算精度、??迭代步數(shù)和CPU時(shí)間的數(shù)值結(jié)果.實(shí)驗(yàn)結(jié)果表明ihADMM算法在計(jì)算時(shí)間上要明顯優(yōu)??于經(jīng)典ADMM算法和LADMM算法,尤其是當(dāng)網(wǎng)格大。练浅Pr(shí).進(jìn)一步,隨著網(wǎng)格??大。璧募蛹(xì),考慮ihADMM算法和LADMM算法的迭代步數(shù),可以看出ihADMM算法??和LADMM算法具有網(wǎng)格獨(dú)立性.相反,經(jīng)典ADMM算法卻不具有網(wǎng)格獨(dú)立性.此外,??為了得到更高精度的解,我們采用了兩階段策略.在表3.2中我們給出了相應(yīng)的數(shù)值結(jié)??果.作為對比,表3.2中我們同樣給出了帶Armijo線搜索的全局化PDAS算法的數(shù)值結(jié)??果.從數(shù)值結(jié)果可以看出,兩階段的算法在CPU計(jì)算時(shí)間以及迭代步數(shù)上要表現(xiàn)地比??帶Armijo線搜索的全局化PDAS算更有優(yōu)勢.這些數(shù)值結(jié)果都說明了我們的ihADMM算??法在求解例3.1到中等精度的解的高效性.??
根據(jù)(4.14),很容易看到這樣的一個(gè)事實(shí):當(dāng)|/?丨</?時(shí)意味著h?=?0,這也就解釋了L1-??范數(shù)為什么可以誘導(dǎo)控制w的稀疏性質(zhì)?此外,由于(4.14)同時(shí)也暗不??了we妒㈦).這樣根據(jù)非光滑方程(4.14),我們在圖4.1中展示了w和p的關(guān)系.??w個(gè)??b——JT??aa-?f3?/?j??? ̄?j? ̄?^?P?ab?+?p?P???i?a??圖4.1控制變量W和伴隨狀態(tài)變量/7的關(guān)系??Fig.?4.1?The?relationship?between?u?and?p??因此,從圖4.1我們可以看出,當(dāng)f充分大時(shí),最優(yōu)控制《將恒等于0,即%?=?0.具體??地,我們有如下引理.??引理4.2若正則參_3滿足條件yS?2灼:=-?A-%,.)丨L,)時(shí),那么問題(P)的唯一??最優(yōu)解為(>¥,¥)?=?(A_1Sy,.,?0).??-99-??
【參考文獻(xiàn)】:
期刊論文
[1]SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY B(?)NARD TYPE[J]. Yanzhen Chang School of Mathematics and System Science,Shandong University,Jinan 250100,China Danping Yang Department of Mathematics,East China Normal University,Shanghai 200062,China. Journal of Computational Mathematics. 2008(05)
本文編號:3594950
【文章來源】:大連理工大學(xué)遼寧省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:217 頁
【學(xué)位級別】:博士
【部分圖文】:
圖2.3例2.2:最優(yōu)控制叫在/7?=?2_7的三維圖和平面圖.(a)三維圖.(b)平面圖.深紅色和深藍(lán)色區(qū)??域?qū)?yīng)w/z?=?±30,綠色區(qū)域?qū)?yīng)w;z?=?0??
帶控制約束的最優(yōu)控制問題的精確解是未知的.作為精確解的代替,我們采用網(wǎng)格大小??為/f?=?2_1()上的數(shù)值解作為一個(gè)參考解.??圖3.1給出了網(wǎng)格大小為A?=?2_6時(shí)離散最優(yōu)控制叫的三維圖.表3.1給出了在L2模意??義下控制《的離散誤差的數(shù)值結(jié)果.?dāng)?shù)值結(jié)果證實(shí)了定理3.]給出的理論誤差結(jié)果???表3.1中給出了ihADMM、經(jīng)典ADMM算法和LADMM算法求解例3.1的計(jì)算精度、??迭代步數(shù)和CPU時(shí)間的數(shù)值結(jié)果.實(shí)驗(yàn)結(jié)果表明ihADMM算法在計(jì)算時(shí)間上要明顯優(yōu)??于經(jīng)典ADMM算法和LADMM算法,尤其是當(dāng)網(wǎng)格大。练浅Pr(shí).進(jìn)一步,隨著網(wǎng)格??大。璧募蛹(xì),考慮ihADMM算法和LADMM算法的迭代步數(shù),可以看出ihADMM算法??和LADMM算法具有網(wǎng)格獨(dú)立性.相反,經(jīng)典ADMM算法卻不具有網(wǎng)格獨(dú)立性.此外,??為了得到更高精度的解,我們采用了兩階段策略.在表3.2中我們給出了相應(yīng)的數(shù)值結(jié)??果.作為對比,表3.2中我們同樣給出了帶Armijo線搜索的全局化PDAS算法的數(shù)值結(jié)??果.從數(shù)值結(jié)果可以看出,兩階段的算法在CPU計(jì)算時(shí)間以及迭代步數(shù)上要表現(xiàn)地比??帶Armijo線搜索的全局化PDAS算更有優(yōu)勢.這些數(shù)值結(jié)果都說明了我們的ihADMM算??法在求解例3.1到中等精度的解的高效性.??
根據(jù)(4.14),很容易看到這樣的一個(gè)事實(shí):當(dāng)|/?丨</?時(shí)意味著h?=?0,這也就解釋了L1-??范數(shù)為什么可以誘導(dǎo)控制w的稀疏性質(zhì)?此外,由于(4.14)同時(shí)也暗不??了we妒㈦).這樣根據(jù)非光滑方程(4.14),我們在圖4.1中展示了w和p的關(guān)系.??w個(gè)??b——JT??aa-?f3?/?j??? ̄?j? ̄?^?P?ab?+?p?P???i?a??圖4.1控制變量W和伴隨狀態(tài)變量/7的關(guān)系??Fig.?4.1?The?relationship?between?u?and?p??因此,從圖4.1我們可以看出,當(dāng)f充分大時(shí),最優(yōu)控制《將恒等于0,即%?=?0.具體??地,我們有如下引理.??引理4.2若正則參_3滿足條件yS?2灼:=-?A-%,.)丨L,)時(shí),那么問題(P)的唯一??最優(yōu)解為(>¥,¥)?=?(A_1Sy,.,?0).??-99-??
【參考文獻(xiàn)】:
期刊論文
[1]SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY B(?)NARD TYPE[J]. Yanzhen Chang School of Mathematics and System Science,Shandong University,Jinan 250100,China Danping Yang Department of Mathematics,East China Normal University,Shanghai 200062,China. Journal of Computational Mathematics. 2008(05)
本文編號:3594950
本文鏈接:http://sikaile.net/kejilunwen/yysx/3594950.html
最近更新
教材專著