用雅克比譜配置方法求解分?jǐn)?shù)階最優(yōu)控制問題及收斂性分析
發(fā)布時間:2020-12-09 17:18
本文針對分?jǐn)?shù)階最優(yōu)控制問題,提出Jacobi譜配置方法。首先根據(jù)Hamiltoni-an 量將最優(yōu)控制問題對應(yīng)的極小泛函進(jìn)行改進(jìn),得到分?jǐn)?shù)階最優(yōu)控制問題的必要條件,結(jié)合Caputo和Riemann-Liouville導(dǎo)數(shù)的定義,推導(dǎo)出等價于原問題的方程組;再利用Jacobi譜配置法,求出非線性系統(tǒng)方程組的解,并通過理論分析,得到相應(yīng)的收斂性結(jié)論。最后,通過數(shù)值實驗驗證理論的正確性和方法的有效性。
【文章來源】:湘潭大學(xué)湖南省
【文章頁數(shù)】:42 頁
【學(xué)位級別】:碩士
【部分圖文】:
圖3.1:?2⑴的真解和不同a的數(shù)值解(左);吣)的真解和不同《的數(shù)值解(右)
?t(2<?N<?20)??圖3.2:?a;⑷(左)和A⑷(右)的|丨.||L=〇和||?.?||紀(jì)誤差比較。??的。圖3.2給出了解在||?.?|U〇〇和||?.?意義下的誤差,證實了我們的理論預(yù)測是??可靠的。??例3.4.2我們考慮一個單輸入標(biāo)量系統(tǒng),如下??1?Z*1??Min?J?=?^?(x2(t)+u2(t))dt,??z?Jo?(〇.〇7)??〇D^x(t.)?=u(t),??自由終端條件和初始條件為??x(0)?=?10.??上述對■于a?=?1問題的解析解是??x(t)?=?0.10061268426^?+?9.899387312e-^,??u⑷=0.1006126842(\/^+?-9.899387312(力一l)e_v^_??根據(jù)(3.8)我們有??tD^\{t)?=x(t),??〇D^x(t)?=?-A(t),??x(0)?=?10,?A(l)?=?0.??其中藍(lán)二〇
?t??圖3.3:?:r⑷的真解和不同a的數(shù)值解(左);噸)的真解和不同a的數(shù)值解(右)。??圖3.3中分別給出了真解(a?=?1),及數(shù)值解(a?=?0.9999,?〇;=?0.9,?a?=?0.8,??a?=?0.7)的狀態(tài)變量;r⑷和控制變量w⑴的圖像,可以看出Jacobi譜配置法??在a?=?0.9999的數(shù)值解和真解a?=?1幾乎重合,說明了我們的求解方法是有效的。??例3.4.3我們考慮一個時間變化的線性系統(tǒng),如下??1?r1?〇??Min?J?=?-?(x2⑷?+?ix2⑷)成??2?J〇?(3.68)??〇D^x(t)?=?tx(t)?+?u(t),??自由終端條件和初始條件為??x(0)?=?1.??根據(jù)(3.8)我們有??tD^\{t)?=?x(t)?+?t\{t),??〇Dfx(t)?=?tx(t)?-X(t),??x(0)?=?1,?A(l)?=?0???其中f?=?〇??m⑴=—A⑷.??圖3.4中給出了?〇不同時的狀態(tài)變量rc⑷和控制變量u(t)的變化趨勢
本文編號:2907205
【文章來源】:湘潭大學(xué)湖南省
【文章頁數(shù)】:42 頁
【學(xué)位級別】:碩士
【部分圖文】:
圖3.1:?2⑴的真解和不同a的數(shù)值解(左);吣)的真解和不同《的數(shù)值解(右)
?t(2<?N<?20)??圖3.2:?a;⑷(左)和A⑷(右)的|丨.||L=〇和||?.?||紀(jì)誤差比較。??的。圖3.2給出了解在||?.?|U〇〇和||?.?意義下的誤差,證實了我們的理論預(yù)測是??可靠的。??例3.4.2我們考慮一個單輸入標(biāo)量系統(tǒng),如下??1?Z*1??Min?J?=?^?(x2(t)+u2(t))dt,??z?Jo?(〇.〇7)??〇D^x(t.)?=u(t),??自由終端條件和初始條件為??x(0)?=?10.??上述對■于a?=?1問題的解析解是??x(t)?=?0.10061268426^?+?9.899387312e-^,??u⑷=0.1006126842(\/^+?-9.899387312(力一l)e_v^_??根據(jù)(3.8)我們有??tD^\{t)?=x(t),??〇D^x(t)?=?-A(t),??x(0)?=?10,?A(l)?=?0.??其中藍(lán)二〇
?t??圖3.3:?:r⑷的真解和不同a的數(shù)值解(左);噸)的真解和不同a的數(shù)值解(右)。??圖3.3中分別給出了真解(a?=?1),及數(shù)值解(a?=?0.9999,?〇;=?0.9,?a?=?0.8,??a?=?0.7)的狀態(tài)變量;r⑷和控制變量w⑴的圖像,可以看出Jacobi譜配置法??在a?=?0.9999的數(shù)值解和真解a?=?1幾乎重合,說明了我們的求解方法是有效的。??例3.4.3我們考慮一個時間變化的線性系統(tǒng),如下??1?r1?〇??Min?J?=?-?(x2⑷?+?ix2⑷)成??2?J〇?(3.68)??〇D^x(t)?=?tx(t)?+?u(t),??自由終端條件和初始條件為??x(0)?=?1.??根據(jù)(3.8)我們有??tD^\{t)?=?x(t)?+?t\{t),??〇Dfx(t)?=?tx(t)?-X(t),??x(0)?=?1,?A(l)?=?0???其中f?=?〇??m⑴=—A⑷.??圖3.4中給出了?〇不同時的狀態(tài)變量rc⑷和控制變量u(t)的變化趨勢
本文編號:2907205
本文鏈接:http://sikaile.net/kejilunwen/yysx/2907205.html
最近更新
教材專著