常微分方程初值問題的數(shù)值解法第1頁，課件共62頁，創(chuàng)作于2023年2月實際問題中遇到的微分方程通常很復雜，多數(shù)情況下無法求出解的解析表達式，即使求出解，也常常由于計算量太大而不實用。然而實際問題本身又往往只要求給出其解在一系列點上的近似值，這就要依靠數(shù)值解法。其中稱為李氏常數(shù)。從而保證上面的初值問題的解存在并且唯一。所謂數(shù)值解法，就是對于解存在的區(qū)間上一系列的點，不妨假定第2頁，課件共62頁，創(chuàng)作于2023年2月上面給定的初值問題的數(shù)值解法有個基本特點，稱作“步進式”，即求解的過程是按照節(jié)點的排列次序一步步地向前推進。描述這類算法，只須在已知的前提下給出計算的遞推公式。逐個求出的近似值。稱為給定的微分方程初值問題的數(shù)值解。相鄰兩個節(jié)點的間距稱為步長。一般我們總假定，即節(jié)點間是等距的。第3頁，課件共62頁，創(chuàng)作于2023年2月其中為的已知函數(shù)，是給定的常數(shù)，求(1.1)、(1.2)的數(shù)值解。一、方法（一）、方法給定初值問題(1.1)(1.2)方法是解初值問題（1.1）、(1.2)最簡單的數(shù)值解法。由于它的精確度不高，實際計算中已不被采用，然而它在某種程度上卻反映了數(shù)值解法的基本思想。第4頁，課件共62頁，創(chuàng)作于2023年2月這種方法是借助于幾何直觀得到的。由于表示解的曲線通過點，并且在該點處以為切線斜率，于是設想在附近，曲線可以用該點處的切線近似代替，切線方程為第5頁，課件共62頁，創(chuàng)作于2023年2月圖6.1第6頁，課件共62頁，創(chuàng)作于2023年2月也就是說，時，可用近似代替，記這個值為，即于是給出了一種當時，獲得函數(shù)值的近似值的方法。重復上面的作法，在處，就可以得到的近似值第7頁，課件共62頁，創(chuàng)作于2023年2月依此下去，當已經(jīng)得到，則這就是著名的方法的計算格式。

由于方法是用一條折線近似地代替曲線，所以方法也叫折線法。一種計算格式，當在計算時，僅僅用到它前一步的信息，稱它為單步法?？梢姺椒ň褪菃尾椒?。第8頁，課件共62頁，創(chuàng)作于2023年2月將方程(1.1)在區(qū)間上求積分，便得到(1.4)式中右端的積分，可以用數(shù)值積分法計算它的近似值。例如，使用左矩形公式則有（二）改進的方法(1.4)上式右端就是用方法得到的，即第9頁，課件共62頁，創(chuàng)作于2023年2月一般地有這就是公式(1.3)。由此可見，方法也可以看成用矩形公式近似計算某個相應的定積分而得到的。因此可以說，方法之所以精確度不高，正是由于它在計算定積分時，采用矩形公式的緣故。倘若使用較為精確的梯形公式來計算(1.4)式中右端的積分，即第10頁，課件共62頁，創(chuàng)作于2023年2月將它代入(1.4)式的右端，便得到的近似值，用同樣的方法可以得到。一般地有，(1.5)第11頁，課件共62頁，創(chuàng)作于2023年2月這就是改進的Euler方法的計算格式。值得注意的是，Euler方法與改進的Euler方法在計算上有一個明顯的區(qū)別，Euler方法中是由已知的或已經(jīng)算出的量來表達的，得到它不需要解方程，這類方法通常稱為顯示方法；而在改進的Euler方法中，未知數(shù)也隱含在方程右端之中，對于每一個的值都需要通過解方程才能得到，這類方法通常稱為隱式格式。在多數(shù)情況下，要從隱式格式(1.5)中解出是很困難的。因此，通常采用如下的迭代方法來求解。即先用Euler方法算出一個結果，作為(1.5)式的初值，進行迭代，其計算格式為第12頁，課件共62頁，創(chuàng)作于2023年2月(1.6)由可知，當時，迭代格式收斂。也就是說，只要取得充分小，就可能保證迭代序列第13頁，課件共62頁，創(chuàng)作于2023年2月收斂，而且越小，收斂得越快。容易看出，改進的方法雖然提高了精度，然而每一步的計算量卻增加很大，每迭代一次，都要重新計算函數(shù)值，而且迭代需要反復進行若干次。為了簡化算法，通常只迭代一次。具體地講，先用方法求得一個初步的近似值，稱為預估值，再將它代入(1.5)式中作一次校正，這樣處理后，計算格式為(1.7)第14頁，課件共62頁，創(chuàng)作于2023年2月稱它為預估校正格式。可用其中的第一式算出一個預估值，再代入第二式做校正。例1用方法和預估校正法求解初值問題取步長。解分別使用格式與預估校正格式計算，格式的具體形式為第15頁，課件共62頁，創(chuàng)作于2023年2月計算結果見下表。預估校正格式格式第16頁，課件共62頁，創(chuàng)作于2023年2月1.73211.73791.78481.01.41421.41641.43510.51.67331.67821.71780.91.34161.34341.35820.41.61251.61531.64980.81.26491.26621.27740.31.54921.55251.58030.71.18321.18411.19180.21.48321.48601.50900.61.09541.09591.10000.1準確解預校方法方法準確解預校方法方法上面給出的初值問題有解析解，按該式算出的準確值與近似值一起列在上表中，通過比較可以看出方法的精度是較低的。第17頁，課件共62頁，創(chuàng)作于2023年2月二、展開法與截斷誤差利用展開法可以得到初值問題的任意高精度的計算格式。設初值問題有解，且，足夠光滑，則在點處的展開式為展開法（一）第18頁，課件共62頁，創(chuàng)作于2023年2月其中值問題中的函由于足夠光滑，則當時，，式中的各階導數(shù)可由初數(shù)來表達，即第19頁，課件共62頁，創(chuàng)作于2023年2月我們在式右端截取項，即舍去余項，則算得的近似值，即此式稱為階的公式。第20頁，課件共62頁，創(chuàng)作于2023年2月（二）局部截斷誤差及其“階”這個截斷誤差被稱為是階的，即當時，是關于的階無窮小量。在考察計算公式的精度時，我們常常假定第步的結果是精確的，即，在這一前提下，來估計第步計算結果的誤差，即，這一誤差稱為局部截斷誤差。例如，階的公式的第步的局部截斷誤差為第21頁，課件共62頁，創(chuàng)作于2023年2月定義1如果一種方法的局部截斷誤差是階的，則稱該方法是階的。由定義1，階公式是階方法，當時，式變?yōu)檫@正是格式，故知格式是一階方法，其局部截斷誤差為，即為二階的。第22頁，課件共62頁，創(chuàng)作于2023年2月例2證明改進的格式是2階方法。對于方法的“階”和局部截斷誤差的“階”,我們可以這樣來理解：如果式的局部截斷誤差是階的，這說明公式的前步的計算結果都是精確的，即式右端關于次的多項式與左端的在處的級數(shù)的次數(shù)不超過的項，完全重合，而兩端超過次的項不重合。因此，我們稱此方法為階的。第23頁，課件共62頁，創(chuàng)作于2023年2月將左端的與右端的在處作展開，有證明設是初值問題、的精確解，即有，由改進的格式有第24頁，課件共62頁，創(chuàng)作于2023年2月將它們代入式，并將右端稍加整理，有可見，該式兩端的前三項，即的次數(shù)不超過的項完全重合，而從的次方的項開始就不重合了。于是，由定義可知，改進的格式是階方法，而其局部截斷誤差是階的。第25頁，課件共62頁，創(chuàng)作于2023年2月解直接求導數(shù)，有例用展開法求解例中的初值問題。第26頁，課件共62頁，創(chuàng)作于2023年2月1.26941.26490.31.18321.18320.21.09541.09540.1用階公式，取步長，部分計算結果列于表中。表中表示準確值，與比較，可見用階公式得到的數(shù)值解是令人非常滿意的。第27頁，課件共62頁，創(chuàng)作于2023年2月三、方法方法（簡稱方法）是一種構造高精度計算公式的方法。前面我們看到，用展開法確實可以得到高精度的計算公式。然而，方法每提高一階，都要增加很大的計算導數(shù)的工作量，而方法，避開了導數(shù)的計算，采用了另外一種構造格式的途徑。第28頁，課件共62頁，創(chuàng)作于2023年2月首先，從微分中值定理及方程得。這里稱為方程的積分曲線在區(qū)間上的平均斜率。由此可見，只要對此平均斜率提供一種算法，就可以得到一個相應的計算公式。下面，我們來觀察格式和改進的格式，將它們分別寫成方法的基本思想(一)第29頁，課件共62頁，創(chuàng)作于2023年2月前一式是用點處的斜率的平均值來代替平均斜率的，后一式是用兩點上的斜率的平均值來代替平均斜率的。我們已經(jīng)知道格式是1階方法，而改進的格式是2階方法。由此看來，如果在區(qū)間內(nèi)多預報幾個點的斜率值，然后將它們加權平均，以代替上述的平均斜率，就可以構造出更高階的計算公式來。因此，方法的關鍵就在于選擇哪些點上的斜率值，以及如何構造它們的線性組合。第30頁，課件共62頁，創(chuàng)作于2023年2月與格式與改進的格式可以改寫成下面的形式級公式（二）第31頁，課件共62頁，創(chuàng)作于2023年2月舍去誤差項，便得到顯然，若在區(qū)間內(nèi)取個不同的點，記積分曲線在這個點上的斜率分別為，于是我們可以設第32頁，課件共62頁，創(chuàng)作于2023年2月這就是所謂的級階的公式。其中都是待定系數(shù)，并且有待定系數(shù)可用比較系數(shù)的方法求得。即將中的和各都在處展成級數(shù)，然后令兩端關于別的不超過次的同次項的系數(shù)相等，便可求得這些待定系數(shù)。下面以為例，說明待定系數(shù)的求法。當時，由式有第33頁，課件共62頁，創(chuàng)作于2023年2月將式中的與、分別在處作展開，有第34頁，課件共62頁，創(chuàng)作于2023年2月稱為修正的梯形公式。注意，這里用到了二元展開式。將上面的三個展開式代入中，并令兩端的次數(shù)不超過的項的系數(shù)相等，于是得到若取，則可算得，這時，由式得第35頁，課件共62頁，創(chuàng)作于2023年2月稱為修正的矩形公式。以上兩個公式，都是在及的前提之下構造出來的。因此，它們都是級階的公式。若取，則可算得，由式得第36頁，課件共62頁，創(chuàng)作于2023年2月注意，在上面求待定系數(shù)的方程組中，有一個自由參數(shù)，故級階的公式有無窮多個。但是，在這些級公式中，不可能存在高于階的方法。下面，我們給出級公式可以達到的最高階數(shù)：第37頁，課件共62頁，創(chuàng)作于2023年2月標準級階公式依照級階公式的構造過程，我們可以得到更高級高階的公式，其中最常用的就是標準的級階公式，其形式為：例用標準級階公式求解例中給出的初值問題，取。（三）第38頁，課件共62頁，創(chuàng)作于2023年2月解計算公式如下：第39頁，課件共62頁，創(chuàng)作于2023年2月1.7320511.7321401.7378691.051.6124521.6125131.6164740.841.4832401.4832811.4859650.631.341641.3416671.3433600.421.1832161.1832921.1840960.2111100(精確值)(4階R-K方法)(2階R-K方法)將計算結果列于表。第40頁，課件共62頁，創(chuàng)作于2023年2月將表與表的結果相比較，盡管這里步長放大了，但計算的精度卻很高，從而出可以看出選擇方法的重要意義。第41頁，課件共62頁，創(chuàng)作于2023年2月四、線性多步法

前面介紹的幾種方法都是單步法，即在計算時，僅用到它前面一步得到信息。設想，當通過單步法已經(jīng)算出，如何充分地利用這些信息，在計算時獲得較高精度，這就是多步法的基本思想。

假定仍討論本章開始給出的一階微分方程的初值問題第42頁，課件共62頁，創(chuàng)作于2023年2月與其等價的積分方程是前面我們曾使用梯公式，計算式右端的積分，而得到了改進的方法。其實，這也可以理解為是用插值點和的線性插值函數(shù)代替函數(shù)而得到的。由于通常插多項式的次數(shù)越高越精確，所以使我們試圖用高次插值多項式代替，來得到高精度的計算方法。

第43頁，課件共62頁，創(chuàng)作于2023年2月今取和為插值節(jié)點，這時的插值多項式為

第44頁，課件共62頁，創(chuàng)作于2023年2月用代替，便得到的近似值，即

令，并注意到則得

第45頁，課件共62頁，創(chuàng)作于2023年2月第46頁，課件共62頁，創(chuàng)作于2023年2月上式中右端的用代替，就有顯然，這是一個隱式方法，稱式為內(nèi)插公式。第47頁，課件共62頁，創(chuàng)作于2023年2月上面之所以得到的是隱式方法，其原因在于選用了作為插值節(jié)點。例如我們?nèi)『妥鳛椴逯倒?jié)點。這時的插值多項式成為第48頁，課件共62頁，創(chuàng)作于2023年2月用上式代替式右端積分中的，也將得到的近似值，與推導隱式方法的過程類似，可得到如下的顯式公式：我們稱式為外插公式。在討論它們的截斷誤差時，不僅要假定，還要假定‘和，容易證明它們的截斷誤差均為。第49頁，課件共62頁，創(chuàng)作于2023年2月可以單獨使用外插公式，用它每計算一個的值，只需要計算一次的值，計算量小于方法，而它們的截斷誤差為同階。但該方法的明顯不足是開始的幾個值和不能用它算，必須采用其它方法。通常把上面給出的兩個公式聯(lián)合