• 沒有找到結果。

甲、 近似函數值 :

在文檔中 顏 國 勇 2012 (頁 96-108)

• 利用 「微分函數」(differential function) , 迅速但不精細.

• 利用 「Taylor 展開式」 :Taylor 定理中, 若函數fNδ(p)n 階可微,

∀ x ∈ Nδ(p), ∃ c ∈ (p, x) 或(x, p) 使得 f (x) = f (p) + f(p)(x − p) + f′′(p)

2! (x − p)2+ · · · +f(n−1)(p)

(n − 1)!(x − p)n−1+f(n)(c)

n! (x − p)n

= Pn−1(x) + Rn(x).

我們將以 Pn−1(x) 做為 f (x) 之一近似值, 而視 Rn(x) 則為誤差. (參見拙著 《 微積分學.電子書 2.1 版》 第四章第一節及第五節. )

[註] : 利用 「微分函數」 方法等於只利用 P1(x) 做為 f (x) 之近似值, 不討論誤差 之控制.

乙、

近似根值

:Newton 迭代法,f (x) = 0之根的近似值, 應有以下諸步驟. 若欲解方程式 f (x) = 0之一近似根至小數第 k.

1 先求 f(x) ;

2g(x) = x − f (x) f(x) ;

3 利用 Bolzano 定理或繪圖法, 估計欲求之根r 的大略位置, 令之為 x1. 其次求 x2 = g(x1) , 至小數k + 1,

x3 = g(x2) , 至小數k + 1, ...

xn+1 = g(xn) , 至小數k + 1, 直至 xn+1xn 之前k + 1 位小數相同時停止 . 4 由上述 xn 以四捨五入法求出 k 位小數之近似值 s .

5 驗算 :ǫ = 10−k/2 ,f (s + ǫ)f (s − ǫ) 異號,s確為吾人所需. 理由 為

s − ǫ < r < s + ǫ ⇔ −ǫ < r − s < ǫ ⇔ |r − s| < ǫ = 10−k/2.

◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ • ◦ ⋆ •

12 單元 近似值問題 97

♠ (12-1) 試利用 「微分函數」(differential function) 以求 8.98 之近似值.

分析 當 |h| 很小時, f(p) · h 可做為 f (p + h) − f(p) 之近似值. 但當 |h| 並非很小時, 這種方法並不精確, 應改以Taylor 定理方法解決.

(不少學者習慣以 dy = ydx 做為 ∆y =

def f (x + dx) − f(x) 之近似值. ) 解 設f (x) =x, p = 9, h = −0.02,則

√8.98 = f (8.98) = f (9 − 0.02)

= f (p + h)

≈ f(p) + f(p) · h

= f (9) + f(9) · (− 0.02)

=√

9 + 1 2√

9 · (− 0.02)

= 2.9967.

♠ (12-2) 試以下列方法求 3之近似值, 使誤差之絕對值小於 5 · 10−5 . (a) 以 Taylor 展開式求之 ;

(b) 以 Newton 迭代法求之 .

分析 以 Taylor 展開式求近似值時, 必須令餘項之絕對值 |Rn(x)| 小於指定誤差 5 · 10−5, 再以逆推法求出 n , 最後再以n − 1Taylor 多項式 Pn−1(x)做為原題之近似值.(a)3 =4 − 1 , 故令 f : (0, +∞) → R : f(x) =√

x , 則 f(x) = 1

2x−1/2, f(4) = 1

2· 1 2, f′′(x) = 1

2· −1

2 x−3/2, f′′(4) = 1 2 · −1

2 · 1 23, f′′′(x) = 1

2 ·−1 2 · −3

2 x−5/2, f′′′(4) = 1 2· −1

2 · −3 2 · 1

25, ...

f(n−1)(x) = 1 2 · −1

2 · −3

2 · · ·−(2n − 5)

2 x−(2n−3)/2, f(n−1)(4) = 1

2· −1 2 · −3

2 · · ·−(2n − 5)

2 · 1

22n−3, f(n)(x) = 1

2 · −1 2 · −3

2 · · ·−(2n − 3)

2 x−(2n−1)/2 = (−1)n−11 · 3 · · · (2n − 3) 2n· x(2n−1)/2 .

是以f 在點 p = 4Taylor 展開式為 f (x) = f (4) + f(4)(x − 4) + f′′(4)

2! (x − 4)2+ · · · +f(n−1)(4)

(n − 1)!(x − 4)n−1+ Rn(x)

= 2 + x − 4

22 −(x − 4)2 2! 25 + · · ·

+ (−1)n−21 · 3 · · · (2n − 5)

(n − 1)! 23n−4 (x − 4)n−1+ Rn(x), 其中

Rn(x) = f(n)(c)

n! (x − 4)n = (−1)n−11 · 3 · · · (2n − 3)

2n· n! c(2n−1)/2 (x − 4)n. c 介於 4x 之間,x = 3 , 餘項為誤差, 由題意應有,

|Rn(3)| < 5 · 10−5

(−1)n−11 · 3 · · · (2n − 3)

2n· n!c(2n−1)/2 (3 − 4)n

< 5 · 10−5

⇐ 1 · 3 · 5 · · · (2n − 3)

2n· n! 3(2n−1)/2 < 5 · 10−5,



∵c ∈ (3, 4) ⇒ 1 c < 1

3



⇐ 1 · 3 · 5 · · · (2n − 3)

2n· n! 3n−1 < 5 · 10−5, 

∵3(2n−1)/2> 3n−1

⇐ n ≥ 7 , (令n = 1, 2, · · · 逐一測試. ) 此時,

√3之近似值 = P6(x) = 2 − 1

1! · 22 − 1

2! · 25 − 1 · 3

3! · 28 −1 · 3 · 5

4! · 211 − 1 · 3 · 5 · 7 5! · 214

−1 · 3 · 5 · 7 · 9 6! · 217

= 1.7321.

(b) 設 f : R → R : f(x) = x2 − 3, 則 f (x) = 0 之一根為 31/2, 誤差之絕對值小於 5 · 10−5 代表近似根至小數第 4. 因為 f(31/2) = 2 · 31/2 6= 0 ,f′′(x) = 2 為連續, 滿足牛頓迭代法之條件.

g(x) = x − f (x)

f(x) = x −x2− 3

2x = x2+ 3 2x , 取 x1 = 2 ,

x2 = g(x1) = 1.75, x3 = g(x2) = 1.73214, x4 = g(x3) = 1.732105, x5 = g(x4) = 1.732105.

12 單元 近似值問題 99 令近似根 s = 1.7321 , 則因 f (1.73205) < 0 , f (1.73215) > 0 , 知 s = 1.7321 為所欲求.

♠ (12-3)

(a) 試求 cosnMaclaurin 展開式, 其中 n 為奇數;

(b) 利用 (a) 以求 cos 5 之近似值至第六位小數(即誤差之絕對值小於 5 · 10−7).

分析 可能有人認為要解決 cos 5 之近似值只要打開計算器或電腦即可, 不需如此大費周章, 其實, 電腦中正是以 Taylor (or Maclaurin) 展開式解決諸特殊函數之近似值.

(a) 令 f (x) = cos x , 則

f (x) = cos x, f (0) = 1,

f(x) = − sin x, f(0) = 0, f′′(x) = − cos x, f′′(0) = −1, f′′′(x) = sin x, f′′′(0) = 0, f(4)(x) = cos x, f(4)(0) = 1,

... ...

f(2k)(x) = (−1)kcos x, f(2k+1)(x) = (−1)k+1sin x, 顯然 cos 在點 0之奇數階導數皆為 0.

n 為奇數, 令 n = 2k + 1, k ∈ N , 則∀x ∈ R, ∃ c ∈ (0, x) or (x, 0) 使得 f (x) = f (0) +f(0)

1! x + · · · + f(2k−1)(0)

(2k − 1)!x2k−1+f(2k)(0)

(2k)! x2k+ R2k+1(x)

= 1 − x2 2! +x4

4! − x6

6! − · · · + (−1)k

(2k)! x2k+(−1)k+1sin c (2k + 1)! x2k+1. 此乃為 sin xn = 2k + 1Maclaurin 展開式.

(b) 由題意 誤差之絕對值 小於 5 · 10−7,x = 5 = 5 · π 180 = π

36, 應有

|R2k+1(x)| < 5 · 10−7

sin c (2k + 1)!(π

36)2k+1

< 5 · 10−7

⇐ 1

(2k + 1)!(π

36)2k+1< 5 · 10−7

⇐ 107 5 (π

36)2k+1< (2k + 1)!

⇐ k ≥ 2 , (註)

(註: 利用計算器,k = 1, 2逐一測試, 此處 k = 2 即已成功. )k = 2,

P4(0.0872665) = 1 − (0.0872665)2

2 + (0.0872665)4

24 = 0.9961947.

cos 5 六位小數之近似值為 0.996195.

♠ (12-4)

(a) 試求 f (x) = sinh x 之 nMaclaurin 展開式;

(b) 試求 sinh 1 之一近似值, 使其誤差之絕對值小於 10−4.

分析 由給定 「誤差之絕對值」 範圍以求項數 n, 常遭遇 |f(n)(c)| 之上界如何求得, 一 旦克服, 困難也可迅速解決.

(a) 由於

f (x) = sinh x, f (0) = 0, f(x) = cosh x, f(0) = 1, f′′(x) = sinh x, f′′(0) = 0, f′′′(x) = cosh x, f′′′(0) = 1, f(4)(x) = sinh x, f(4)(0) = 0,

... ...

f(2k−1)(x) = cosh x, f(2k−1)(0) = 0, f(2k)(x) = sinh x,

f 之 n = 2k, (k ∈ N), 階 Maclaurin 展開式為 ∀ x ∈ R, ∃ c ∈ (0, x) 或 (x, 0) 使得

f (x) = f (0) + f(0)x + f′′(0)

2! x2+ · · · + f2k−1(0)

(2k − 1)!x2k−1+f(2k)(c) (2k)! x2k

= x + x3 3! +x5

5! + · · · + x2k−1

(2k − 1)!+ R2k(x).

R2k(x) = sinh c (2k)!x2k.

(註: n 為奇數時, f Maclaurin 展開式同理, 略之. )

12 單元 近似值問題 101

(b) 由原設誤差之絕對值 小於 10−4

|R2k(1)| < 10−4 ⇐ sinh c

(2k)! < 10−4

⇐ 1.5

(2k)! < 10−4,



∵ 0 < sinh c < sinh 1 = e − e−1 2 < e

2 < 3 2



⇐ k ≥ 4.

亦即 sinh 1 之近似值為 1 + 1 3! + 1

5! + 1

7! = 1.1752.

♠ (12-5)

(a) 試求 f (x) = ln(1 + x) 之nMaclaurin 展開式 ;

(b) 利用 (a)ln 2 之近似值時, 若欲使誤差之絕對值小於 10−5. 試問 n 至少應為 多少 ? 並將近似值表為有限級數 (不必化簡).

分析 題目之用意在於 : 測試讀者使用逆推法以求 「項數」n 之能力.(a) 因為

f(x) = (1 + x)−1, f(0) = 1, f′′(x) = −(1 + x)−2, f′′(0) = −1, f′′′(x) = 2!(1 + x)−3, f′′′(0) = 2!, f(4)(x) = −3!(1 + x)−4, f(4)(0) = −3!,

... ...

f(n−1)(x) = (−1)n(n − 2)!(1 + x)−n+1, f(n−1)(0) = (−1)n(n − 2)!, f(n)(x) = (−1)(n+1)(n − 1)!(1 + x)−n.

對於 x > 0 , (或 −1 < x < 0 ) , 由 Taylor 定理知, 存在 c ∈ (0, x), ( 或 c ∈ (x, 0) ) 使得

ln(1 + x) = f (0) + f(0)

1! x +f′′(0)

2! x2+f′′′(0)

3! x3+ · · · +f(n−1)(0)

(n − 1)! xn−1+ f(n)(c) n! xn

= x −1

2x2+1

3x3− · · · + (−1)n 1

n − 1xn−1+ (−1)n+1 1

n(1 + c)nxn. 此為 f (x) = ln(1 + x) 之nMaclaurin 展開式.

(b) 由題意,x = 1 , 因誤差之絕對值< 10−5, 應有

|Rn(1)| < 10−5 ⇔ 1

n(1 + c)n < 10−5

⇐ 1

n ≤ 10−5, ( ∵ 1/(1 + c)n < 1 )

⇐ n ≥ 105. 取 n = 100, 000 , 此時 ln 2 之近似值為

Pn−1(1) = 1 − 1 2+ 1

3− · · · + 1 99999.

(

: 由上述答案知, 以此一方法求 ln 2 之近似值, 似乎不太高明, 下一題有更

好的方法. )

♠ (12-6) (續上題) 利用ln 1 + x 1 − x

之Maclaurin 展開式, 試求ln 2 之近似值, 使其誤差 之絕對值小於 10−5. (提示 : 令 x = 1/3 ) .

分析 上題為求 ln 2 之近似值, 使用 ln(1 + x) , 需要一個 99999 項之級級, 改用本法只 需不到十項.

解 令f : (−1, 1) → R : f(x) = ln(1 + x) − ln(1 − x) . 則∀x ∈ (−1, 1), f(x) = (1 + x)−1+ (1 − x)−1,

f′′(x) = −(1 + x)−2+ (1 − x)−2, f′′′(x) = 2(1 + x)−3+ 2(1 − x)−3, f(4)(x) = −3!(1 + x)−4+ 3!(1 − x)−4,

...

f(n−1)(x) = (−1)n(n − 2)!(1 + x)−(n−1)+ (n − 2)!(1 − x)−(n−1), f(n)(x) = (−1)n+1(n − 1)!(1 + x)−n+ (n − 1)!(1 − x)−n,

f 在點之Maclaurin 展開式為, ( n 為奇數) f (x) = f (0) + f(0)x + f′′(0)

2! x2+ f′′′(0)

3! x3+ · · · + f(n−1)(0)

(n − 1)!xn−1+ Rn(x)

= 0 + 2x + 0 + 2x3

3 + 0 + · · · + 2xn−2

n − 2 + 0 + Rn(x)

= 2 x + x

3 + · · · + xn−2 n − 2



+ Rn(x), 內Rn(x) = 1

n

 1

(1 + c)n + 1 (1 − c)n

xn, 而c 介於0 與 x.x = 1/3 , 是以

ln 2 = f 1 3

= 2 1 3 + 1

3 · 3 + · · · + 1 (n − 2) · 3n−2

+ Rn

 1 3

.

12 單元 近似值問題 103

Rn

 1 3



< 10−5 ⇔ 1 n

 1

(1 + c)n + 1 (1 − c)n

 1

3n < 10−5

⇐ 1

n3n



1 + 1

(1 −13)n



< 10−5

⇐ n ≥ 13.

n = 13 ,ln 2 之近似值為 P12 = 2 1

3+ 1

3 · 33 + 1

5 · 35 + 1

7 · 37 + 1

9 · 39 + 1 11 · 311



= 0.69315.

♠ (12-7) 設L乃過點(1, 0)且為曲線 y = x2 之法線, 此法線與曲線之交點為(x0, x20) , 試 求x0 之近似值至小數第四位.

分析 本題需先求曲線外之點 (1,0) 的法線, 最後以牛頓迭法解出 x0 之近似值.

解 設 F (x) = x2,F(x) = 2x .(1,0) 和交點為 (x0, x20) 連線之斜率為 x

20

x0− 1, 而過點(x0, x20) 之切線斜率為 2x0, 應有

x20

x0 − 1 = − 1 2x0

, 亦即 2x30+ x0− 1 = 0 . 以牛頓迭代法解x0 之近似值.

1f (x) = 2x3+ x − 1,f(x) = 6x2+ 1 . (注意 : fF 勿相混. )

2

g(x) = x − f (x)

f(x) = x −2x3+ x − 1 6x2+ 1 . 3 由圖中我們可看出2x3+ x − 1 = 0 之根 x0 應在0.5

附近,x1 = 0.5,

x2 = g(x1) = 0.6, x3 = g(x2) = 0.58975, x4 = g(x3) = 0.58975,

4 上述 x3 之小數第五位四捨五入得 s = 0.5898 .

12–1

5 驗算 : 因 f (0.58985) > 0, f (0.58975) < 0 , 故知 x0 之近似值為 0.5898.

♠ (12-8) 設曲線 y = sin x 與y = exp(−x) 於最接近原點之交點為 (x0, y0).

(a) 繪出此二曲線及其交點之簡圖 ; (b) 試求 x0 之近似值至小數第五位 .

分析 繪圖之後得知交點之橫坐標在 0.5附近, 可用以為牛頓迭代法之初始點, 解題時則應 令f (x) = sin x − e−x, 而尋找f (x) = 0 之根.

(a) 此二曲線及其交點之簡圖如下 :

12–2

(b) 設 f (x) = sin x − e−x, 則 f(x) = cos x + e−x, 令 g(x) = x − f (x)

f(x) = x − sin x − e−x cos x + e−x, 利用牛頓迭代法求 f (x) = 0之根 x0,x1 = 0.5 ,

x2 = g(x1) = 0.585644, x3 = g(x2) = 0.588529, x4 = g(x3) = 0.588533, x5 = g(x4) = 0.588533,

令 s = 0.58853 , 則因 f (0.588525) < 0 , f (0.588535) > 0 , 故知 x0 之五位小 數近似值為 0.58853.

♠ (12-9) 設過點 (−1, 0) 且與自然對數函數 ln 相切之切點為 (a, ln a) , 試以 Newton 迭代 法, 求之 a 之近似值至第四位小數.

分析 為了方便, 微積分教師常設計相對容易計算的題目, 然而應用上又常遇到難以處理的 問題, 此時, 近似值就成了最佳工具.

12 單元 近似值問題 105

12–3

解 切點為(a, ln a) , 切線之斜率為 ln(a) = 1/a , 故切線為 y − ln a = 1

a(x − a) , 因其 經過點(−1, 0) , 是以

0 − ln a = 1

a(−1 − a),

亦即a為方程式a ln a − a − 1 = 0 之根, 令f (x) = x ln x − x − 1 , 我們將以Newton 迭代法求 f (x) = 0之近似根.

1g(x) = x − f (x)

f(x) = x −x ln x − x − 1

ln x = x + 1 ln x . 2x1 = 4 ,

x2 = g(x1) = 3.60673, x3 = g(x2) = 3.59115, x4 = g(x3) = 3.59115. 3 上述 x4 四捨五入得 a 之四位小數之近似值為 s = 3.5911 .

4 驗算 : 因 f (3.59105) > 0, f (3.59115) < 0 , 故知 a 之第四位小數近似值為 3.5911.

♠ (12-10) 某君於某年元旦起每月初存款M 元於某銀行, 存入12次後於翌年元旦領回本利

, 假設利息係按月以複利計算 .

(a) 設月利率為 r , 試將本利和A 表為 rM 之函數;

(b) 若 M = 900, A = 12, 000 , 試利用 Newton 迭代法以求 r 之近似值至小數第四 位.

分析 在華人地區許多人參加 「互助會」, 例如某甲為某一互助會之會腳, 自某年元旦起每 個月月初將 「活會款」 900 元交給 「會頭」, 交了十二次, 某甲乃 12 名 「會腳」 最後一 位, 第二年元旦 「會頭」 將收來的會款 (每位 「死會」 付 1000) 共計 12,000 元交給 某甲, 此一情況與 (b) 題相同.

(a) A = M X12 k=1

(1 + r)k = M(1 + r)((1 + r)12− 1) (1 + r) − 1 . (b) 為方便計令a = 1 + r , 則

A = 12000 ⇔ 900a(a12− 1)

a − 1 = 12000

⇔ 3a13− 3a = 40a − 40

⇔ 3a13− 43a + 40 = 0.

f (a) = 3a13− 43a + 40 , 求f (a) = 0 之近似根.f(a) = 39a12− 43 , g(a) = a − f (a)

f(a) = a − 3a13− 43a + 40 39a12− 43 , 並取 a1 = 1.02 ,

a2 = g(a1) = 1.016159, a3 = g(a2) = 1.0161037, a4 = g(a3) = 1.0161037,

令 s = 1.0161 , 則因 f (1.01615) > 0, f (1.01605) < 0, 知 sa 之四位小數近 似值, 即月利率 r 之四位小數近似值為 0.0161 = 1.61% .

[註] : 這麼高的利率只發生在高通膨的年代.

♠ (12-11) 設一平行於X軸之直線與y = cosh x相交而得一 弦, 試求此弦之長度之近似值至第四位小數, 以使弦之 中點 By = cosh x 之頂點 A 之距離等於弦長, 亦 即AB = QP (參見右圖).

[

] : 美國密蘇里州的聖路易大拱門 (St. Louis Gate-way Arch) 即為一倒立之懸索線, 高與寬皆為 630. 本題解出點Px坐標 x0, 先倒立, 再將 x, y 皆 予放大 k = 630

2 · x0

= 127.711 倍, 最後向上平移, 即 可算出拱門函數 :

g : [−315, 315] → R : g(x) = k − k · coshx

k + 630.

12–4

分析 若設點 P 之坐標為 (x0, y0) , 由原設 y0− 1 = 2 x0, 再設法解出x0.

12 單元 近似值問題 107 解 設點 P 之坐標為 (x0, y0) , 由題意知應解方程式 2x0 = cosh x0 − 1 , 此方程式不易 以邏輯方式解出 a , 令 f (x) = cosh x − 2x − 1 , 我們改以 Newton 迭代法解方程式 f (x) = 0 . 因f(x) = sinh x − 2 , 故令

g(x) = x − f (t)

f(t) = x − cosh x − 2x − 1 sinh x − 2 , 令x1 = 3 , 利用計算器或電腦可求得

x2 = g(x1) = 2.61739 x3 = g(x2) = 2.482005 x4 = g(x3) = 2.466680 x5 = g(x4) = 2.466497 x6 = g(x5) = 2.466497.

令s = 2.46650 , 因 f (2.466495) < 0, f (2.466505) > 0 , 故知 x0 之五位小數近似值 為2.46650 , 即弦長之四位小數近似值為 2x0 = 4.9330 .

(註 : x0 近似值多取一位小數, 以免 2x0 之末位造成誤差. )

在文檔中 顏 國 勇 2012 (頁 96-108)

相關文件