本文關(guān)鍵詞：高階Runge-Kutta方法的構(gòu)造及研究

  更多相關(guān)文章： 顯式對(duì)稱Runge-Kutta方法 隱式辛Runge-Kutta方法 根樹理論 W-變換 辛方法

【摘要】：常微分方程在自然科學(xué)的很多學(xué)科領(lǐng)域都有著重要的應(yīng)用，如自動(dòng)控制、電子學(xué)裝置的設(shè)計(jì)、彈道計(jì)算、飛機(jī)和導(dǎo)彈飛行的穩(wěn)定性以及化學(xué)反應(yīng)過程穩(wěn)定性的研宄等都可以轉(zhuǎn)化為求解常微分方程或研究常微分方程解的性質(zhì).當(dāng)前，，計(jì)算機(jī)技術(shù)的快速發(fā)展為常微分方程的應(yīng)用及理論研究提供了有力的工具.然而，只有一些特殊類型的常微分方程初值問題能夠得到用解析式表示的精確解，大量的常微分方程初值問題很難得到其精確解的解析式，有的甚至根本無法用解析式來表示，因此我們只能依賴于數(shù)值方法，以獲得常微分方程初值問題的數(shù)值解.Runge-Kutta方法是求解常微分方程初值問題的經(jīng)典方法.本文考慮構(gòu)造求解常微分方程初值問題的高階隱式辛Runge-Kutta方法和8級(jí)6階顯式對(duì)稱Runge-Kutta方法.在第三章中，我們主要研究了高階隱式對(duì)稱辛Runge-Kutta方法的構(gòu)造.對(duì)于高階隱式對(duì)稱辛Runge-Kutta方法的構(gòu)造，我們利用了Hairer和Wanner教授提出的W-變換理論，通過在變換矩陣:中選取不同的參數(shù)α,β和γ的值，我們不僅得到了經(jīng)典的Gauss方法、LobattoⅢA、LobattoⅢB、LobattoⅢC、Lobatto ⅢE、Lobatto ⅢS，而且我們得到了一類新的隱式對(duì)稱辛Runge-Kutta方法(我們稱之為Lobatto ⅢSX方法).我們構(gòu)造的Lobatto ⅢSX方法的辛性質(zhì)可以由Hairer和Wanner教授給出的相關(guān)結(jié)論直接得到.在第四章中，我們主要研究了8級(jí)6階顯式對(duì)稱Runge-Kutta方法的構(gòu)造.對(duì)于高階顯式Runge-Kutta方法構(gòu)造的困難性和復(fù)雜性可以從表格(2-1)輕易地得出.為了克服己有高階顯式方法構(gòu)造的困難性，我們從8級(jí)顯式Runge-Kutta方法的一般形式出發(fā)，借助Runge-Kutta方法的一般伴隨方法(在本文中，我們稱之為對(duì)稱伴隨方法)、辛伴隨方法(新定義的)給出了一種構(gòu)造8級(jí)6階顯式對(duì)稱Runge-Kutta方法的新途徑.首先，我們利用8級(jí)顯式Runge-Kutta方法(4.2.1)及其辛伴隨方法(4.2.4)得到相應(yīng)的8級(jí)對(duì)稱辛方法(4.2.5)(詳見表4-2),并在A*=AS*及積分公式(b,c)對(duì)稱的條件下,我們給出了簡化階條件B(p)(p ≤ 6),C(1)和D(1)的化簡,再利用得到的C(1)和D(1)的等價(jià)性,證明了對(duì)稱辛方法(4.2.5)對(duì)應(yīng)簡化階條件CS*(1),DS*(1),CS*(2)和DS*(2)的等價(jià)關(guān)系.其次,我們利用積分公式(b,c)的對(duì)稱性、A*= C(1)(?)D(1)、Cs*(1)(?)DS*(1)及Cs*(2)(?)DS*(2)的相關(guān)結(jié)論,對(duì)8級(jí)6階顯式對(duì)稱Runge-Kutta方法需要滿足的33個(gè)階條件分三步進(jìn)行了相應(yīng)的簡化.通過這些簡化過程,我們得到8級(jí)6階顯式對(duì)稱Runge-Kutta方法需要滿足階條件(4.2.55)式.然后,我們將c2,c3,c4,b2,b3,a43(或a42)看作自由變量,并通過求解階條件(4.2.55)式,得到8級(jí)6階顯式對(duì)稱Runge-Kutta方法的Butcher表中系數(shù)a32,a42(或a43),a52,a53,a54,a62,a63,a72和b4的表達(dá)式.最后,我們利用(4.2.6),(4.2.11),(4.2.13),(4.2.14),(4.3.7),(4.3.8),(4.3.9),(4.3.10),(4.3.11),(4.3.12),(4.3.13),(4.3.14)和(4.3.20)等式,通過選擇自由變量a43(或a42),c2和c4的值,構(gòu)造了一類具體的8級(jí)6階顯式對(duì)稱Runge-Kutta方法.
【學(xué)位授予單位】：上海師范大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2017
【分類號(hào)】：O241.8

【相似文獻(xiàn)】

中國期刊全文數(shù)據(jù)庫 前10條

1 ;ORDER RESULTS FOR ALGEBRAICALLY STABLEMONO-IMPLICIT RUNGE-KUTTA METHODS[J];Journal of Computational Mathematics;1999年06期

2 李銀山,張年梅,楊桂通;1/3 SUBHARMONIC SOLUTION OF ELLIPTICAL SANDWICH PLATES[J];Applied Mathematics and Mechanics(English Edition);2003年10期

3 蘇凱;王錦紅;張宏偉;王晚生;;顯式和對(duì)角隱式Rung-Kutta方法求解中立型泛函微分方程的非線性穩(wěn)定性[J];數(shù)值計(jì)算與計(jì)算機(jī)應(yīng)用;2011年01期

4 阮保庚;高階雙參量Runge-Kutta方法[J];數(shù)學(xué)物理學(xué)報(bào);1999年03期

5 楊彪,孫樂平;延時(shí)微分方程多步Runge-Kutta方法的P-穩(wěn)定性(英文)[J];上海師范大學(xué)學(xué)報(bào)(自然科學(xué)版);1997年04期

6 王森;關(guān)于高階Runge-Kutta方法的對(duì)稱性及藕對(duì)性[J];數(shù)學(xué)研究與評(píng)論;1999年03期

7 ;NUMERICAL SIMULATION FOR SOLITARY WAVES OF RLW EQUATION[J];Journal of Hydrodynamics(Ser.B);2004年02期

8 甘四清;Runge-Kutta方法的強(qiáng)正則性[J];長沙鐵道學(xué)院學(xué)報(bào);2000年04期

9 鄒巾英,丁效華,劉明珠;二級(jí)二階和三級(jí)三階連續(xù)Runge-Kutta-Nystr(?)m方法[J];黑龍江大學(xué)自然科學(xué)學(xué)報(bào);2003年02期

10 曹陽,李慶揚(yáng);A-stable Explicit Nonlinear Runge-Kutta Methods[J];Tsinghua Science and Technology;1998年04期

中國重要會(huì)議論文全文數(shù)據(jù)庫 前6條

1 ;Spurious Solutions of Runge-Kutta Methods for Delay Differential Equations[A];第九屆全國微分方程數(shù)值方法暨第六屆全國仿真算法學(xué)術(shù)會(huì)議論文集[C];2004年

2 ;H-stability of Runge-Kutta Methods with General Variable Stepsize for Pantograph Equation[A];第九屆全國微分方程數(shù)值方法暨第六屆全國仿真算法學(xué)術(shù)會(huì)議論文集[C];2004年

3 ;Symplectic Runge-Kutta Methods for The Linear Quadratic Regulator Problem[A];保結(jié)構(gòu)算法2004離散變分和上同調(diào)及其在動(dòng)力系統(tǒng)計(jì)算中的應(yīng)用[C];2004年

4 劉明珠;李冬松;;Runge-Kutta方法對(duì)于比例方程的漸近穩(wěn)定性[A];新世紀(jì) 新機(jī)遇 新挑戰(zhàn)——知識(shí)創(chuàng)新和高新技術(shù)產(chǎn)業(yè)發(fā)展（上冊(cè)）[C];2001年

5 孫耿;;辛Runge-Kutta方法和擾動(dòng)的Hamilton系統(tǒng)[A];保結(jié)構(gòu)算法2004離散變分和上同調(diào)及其在動(dòng)力系統(tǒng)計(jì)算中的應(yīng)用[C];2004年

6 ;High-order time discretizations in seismic modeling[A];中國科學(xué)院地質(zhì)與地球物理研究所2007學(xué)術(shù)論文匯編(第三卷)[C];2008年

中國博士學(xué)位論文全文數(shù)據(jù)庫 前4條

1 夏開封;高階Runge-Kutta方法的構(gòu)造及研究[D];上海師范大學(xué);2017年

2 李計(jì)勇;二階振蕩微分方程數(shù)值方法研究[D];南京大學(xué);2012年

3 吳志橋;L-穩(wěn)定格式求解結(jié)構(gòu)動(dòng)力學(xué)方程和多體系統(tǒng)動(dòng)力學(xué)方程[D];國防科學(xué)技術(shù)大學(xué);2009年

4 蔣成香;高階對(duì)角隱式辛Runge-Kutta方法研究及應(yīng)用[D];上海師范大學(xué);2015年

中國碩士學(xué)位論文全文數(shù)據(jù)庫 前10條

1 邱明明;隨機(jī)時(shí)滯微分方程顯式Runge-Kutta方法的漸近穩(wěn)定性[D];上海師范大學(xué);2015年

2 韓明崗;誤差校正和指數(shù)隨機(jī)Runge-Kutta方法[D];哈爾濱工業(yè)大學(xué);2015年

3 孫瑞;非線性中立型泛函微分方程Runge-Kutta方法的穩(wěn)定性分析[D];長沙理工大學(xué);2014年

4 穆薩（Ibrahim Hussein Musa Tahir）;生物微分方程的保結(jié)構(gòu)數(shù)值方法[D];南京農(nóng)業(yè)大學(xué);2014年

5 錢永峰;線性離散型時(shí)滯系統(tǒng)的多導(dǎo)數(shù)Runge-Kutta方法[D];華中科技大學(xué);2015年

6 金永虎;三角指數(shù)擬合Runge-Kutta方法[D];華中科技大學(xué);2011年

7 劉佳蘭;兩類兩步Runge-Kutta方法的代數(shù)穩(wěn)定性與定量誤差分析[D];湘潭大學(xué);2006年

8 陳全發(fā);Runge-Kutta-Nystr(?)m方法的若干新性質(zhì)與階結(jié)果[D];湘潭大學(xué);2007年

9 布春霞;解常微分方程的三步Runge-Kutta方法[D];華中科技大學(xué);2007年

10 康明;求解剛性振蕩問題的單調(diào)隱式Runge-Kutta方法[D];湘潭大學(xué);2011年

本文編號(hào)：1273706