第五章常微方程初值問題

2023-01-06 16:22:59 字數 3086 閱讀 7033

我們討論的導數問題。實際計算及討論函式在具體某點的導數值。

我們自然想到微積分中的導數定義:

用作為的近似。

理論上這種近似,h越小越好,但實際計算相近的數相減會損失有效數字,實際計算中要注意這些問題。

二階導數的三點公式:

精度高乙個量級

往往利用插值方法先求出插值函式,在對函式求導,用之代替原函式的導數。

數值積分我們也曾這樣做。

本章我們主要解決初值問題的數值解法。即:

數值解就是要求出在

一系列點上的近似值

。假設在所討論問題中連續。由leibniz公式(或對方程積分)有:

依據上式可得到不同的數值解法。

§1 尤拉法和改進尤拉法、預估-校正法

由積分可用左矩形公式,得:

顯式用右矩形公式,得:

隱式梯形公式,得:

改進尤拉公式

精度比上兩個高乙個量級。

隱式和改進尤拉法都是乙個方程。如果方程解時比較複雜,則用迭代發近似求解。

先用顯式尤拉公式計算初值,再用改進尤拉公式迭代求解。

(數學證明當在所討論的範圍有界在h取得很小時,這種迭代收斂)。

公式為:

程式結構:

f(x,y)=y*y

write(*,*) 'input x0,y0'

read(*,*) x0,y0

write(*,*) 'input [x0,b] b'

read(*,*) b

write(*,*) 'input h'

read(*,*)h

y=y0

f0=f(x0,y0)

do 100 x=x0, b, h

k1=k2=k3=

k4=y=y+

y=y0+h*f0

do 200 k=1,10

y=y0+h/2.0*(f0+f(x,y))

200 continue

write (*,*) x+h,y, 1.0/(1-x-h)

y0=y

f0=f(x,y)

100 continue

stop

end注意:所計算的範圍必須是從x0到某值b。因為是逐點遞推的。

當計算範圍為[b,x0]時,步長h取負值。

作業:數值解初值問題:

用改進公式解方程,並和精確解作圖比較。

注意:本題中x的取值範圍不能接近1,為什麼?

將方程換成,精確解為計算[-2,2]間曲線

將方程換成,精確解為計算[-2,2]間曲線

改進尤拉法在判斷迭代收斂時很不好求(為什麼?未知!)。

往往計算中用預估-校正法,公式如下:

左矩形預估

右矩形校正§2 龍格-庫塔法(runge-kutta)

尤拉法計算簡單,精度不高。

預估校正法可寫為:

將上式加以推廣

**為待定係數,當時就是預估-校正法。

我們現在要選取上面的三個係數,使得上式的誤差為。

我們將在展開:

****

其中:另一方面,將 **式中的在已知點展開

二元函式展開

帶入**式得:[, , , , , , , , ]

與****式比較

****

我們取時** 式的誤差為

三個未知數兩個方程,所以具體的取法很多。

當時,就是尤拉預估-校正法。

另乙個常取,此時公式為:

中點公式

以上是二階龍格-庫塔公式,按照相似的方法做出三階公式

%%%為待定係數,選取這五個係數,可使得上式的誤差為。

方法和前面一樣,將在展開到三次項;

同時將%%%式也展開到三次項;

在比較對應項的係數。

最後可得

六個未知數,四個方程,常取如下兩個公式:

四階公式

上式為古典龍格庫塔公式,最常用

程式很簡單。

寫時如果要用k作變數,將其定義成實型。

作業1、將尤拉改進法程式改為四階古典龍格庫塔法,並比較結果。

2、將以前的作業改成用龍格庫塔法來做。

3、考慮如下問題計算方法

4、自己選一系統(氫原子,線性諧振子,無限深勢阱。。。。。),設定某種圍繞,計算相關量子態機率隨時間變化,畫出曲線。注意:表示式並沒有考慮自發躍遷。

類似於:設有五位同學a1、a2、a3、a4、a5,開始a1有十萬元錢,其他人沒錢。沒人每天花費其所有的20%從其他四位同學那裡買東東(相當於圍繞矩陣元)。

計算五位同學擁有的錢數隨時間變化曲線。

§5.4含時間微擾理論

定態微擾理論不含時間,求下的定態近似波函式和能級

本節討論哈密算符含有與時間相關的微擾

即其中與時間無關,僅與時間有關,因含時,

求解體系的波函式要由含時的薛方程準確解出,通常時很困難的

含時間微擾理論,由的定態波函式近似地計算可以無微擾體系()在微擾作用下由乙個量子態躍遷到另乙個量子態的躍遷機率。

設體系哈密頓算符

是微擾已精確求解,本徵函式

將按的定態波函式展開

[, ]

代入薛方程

t': 'span', 'c': '因為', 'r': 'r_32'}, ]

所以上式化為

以左乘上式兩邊,積分

其中上式給出了隨時間的變化規律

即給出了波函式

薛方程的另一種表示方式—矩陣表示

求解方程形式解 *

[, , ]

設微擾在時開始,且這時體系處於的第k個本徵態,則

即當時,

但很小時,或作用時間很短時,

t': 'span', 'c': '零級近似', 'r': 'r_13'}]

代入*得t': 'span', 'c': '一級近似', 'r': 'r_13'}]

在時刻體系處於態的機率,

即從初態躍遷到的機率

①一級近似,波函式用零級近似,積分中

②躍遷的機率

等於的機率

因 ③多次迭代精度更高,數值解最適於此類問題

[, ]如按能級

度簡併初態度簡併1 末能級簡併度則躍遷機率要對末態求和

②初能級簡併即度簡併

初未都簡併

八年級下第五章分式與分式方程

第四節分式方程 二 班級 數學班姓名 李增盛學號 學習目標 1 體會分式方程到整式方程的轉化思想,掌握分式方程的解法 2 了解分式方程產生增根的原因,會檢驗根的合理性 3 培養學生的數學轉化思想和觀察 模擬 探索的能力 學習方法 自主 總結與小組合作交流相結合 學習重難點 重點 掌握分式方程的解法解...

第五章一元一次方程形成性評價

時間 40分鐘滿分 100分 班級姓名得分 一 選擇題 每小題5分,滿分30分 1 下列方程是一元一次方程的是 a x 2y 9 b x2 3x 1 c d 2 方程的解是,則等於 a 8 b 0c 2d 8 3 方程的解是 a b c d 4 方程2 去分母得 a 2 2 2x 4 x 7 b 1...

第五章《二元一次方程組》水平測試

八年級數學 上 第二次月考試題 時間 100分鐘滿分 120分 班級 學號姓名 一 選一選,看完四個選項再做決定!每小題3分,共30分 1 下列不是二元一次方程組的是 a b c d 2 由,可以得到用表示的式子是 a b c d 3 方程組的解是 a b cd 4 方程是二元一次方程,則的取值為 ...