Error Bounds
Error Bounds
電腦裡面要表達一個數字,像是
所以你也可以知道如果我們要電腦幫我們解 Ax = b,那麼必定會有誤差,因為電腦除了整數之外,大部分的數字是幾乎不可能沒有誤差的
Def
我們要解 Ax = b 這個線性系統的時候,假設我們透過迭代方法找到一個解
如果現在只有
那我們就期待如果
我們看個例子,這邊我們解一個 Ax = b:

A 矩陣的第一個 row 和第二個 row 只有一點點差異,如果我們用手去解你會發現他是可逆的,然後解會是
但如果我們今天用某種迭代方法算出一個解長
那為什麼,原因完全出現在 A 矩陣裡面,我們接 Thm 7.27
Thm 7.27
我們假設
然後如果

我們先看一下上面的在說什麼,他的意思是逼近解跟真實解的絕對誤差會小於等於右邊那個東西,而下面那個則是逼近解與真實解的相對誤差會小於等於右邊那個式子
所以一個是絕對誤差,一個是相對誤差
套用到剛剛的例子,就會長這樣:

也就是說雖然

所以關鍵就是在
證明:

Def 7.28

如果 K(A) 接近 1,我們就說這個系統是良置的,如果遠遠大於 1,就說這個系統是病態的
舉個例子:

Thm 7.29
接下來是另一個主題,就像一開始所說我們會有 rouding error,又稱 perturbation,也就是說我們實際上 Ax = b 丟進去電腦時會跑出個

那我們把它就代回 Ax = b:

這時我們就會發現電腦裡面算出來的
證明:
