本文關鍵詞：高階Runge-Kutta方法的構造及研究

  更多相關文章： 顯式對稱Runge-Kutta方法 隱式辛Runge-Kutta方法 根樹理論 W-變換 辛方法

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

【相似文獻】

中國期刊全文數(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 蘇凱;王錦紅;張宏偉;王晚生;;顯式和對角隱式Rung-Kutta方法求解中立型泛函微分方程的非線性穩(wěn)定性[J];數(shù)值計算與計算機應用;2011年01期

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

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

6 王森;關于高階Runge-Kutta方法的對稱性及藕對性[J];數(shù)學研究與評論;1999年03期

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

8 甘四清;Runge-Kutta方法的強正則性[J];長沙鐵道學院學報;2000年04期

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

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

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

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

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

3 ;Symplectic Runge-Kutta Methods for The Linear Quadratic Regulator Problem[A];保結構算法2004離散變分和上同調及其在動力系統(tǒng)計算中的應用[C];2004年

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

5 孫耿;;辛Runge-Kutta方法和擾動的Hamilton系統(tǒng)[A];保結構算法2004離散變分和上同調及其在動力系統(tǒng)計算中的應用[C];2004年

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

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

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

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

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

4 蔣成香;高階對角隱式辛Runge-Kutta方法研究及應用[D];上海師范大學;2015年

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

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

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

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

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

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

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

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

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

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

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

本文編號：1273706