1 of 80

微分方程初步

李柏翰

國立師大附中/leepohan@gmail.com

22 April, 2022

*

1

2 of 80

基礎知識

*

2

BuzzNote ( https://tg.pe/buzzBuzzNote ( https://tg.pe/buzz ; 李柏翰老師自建網站)

李柏翰老師 教學影片

1.Ch1 三角函數(Math0) (https://youtu.be/A050tIZZDvg)

微積分入門(Math1(Math1) (https://youtu.be/Z2MqFSzl8GQ)

微積分教學與應用(Math1a) (https://youtu.be/lHItY_X_m38)

2.Ch2 向量,斜拋運動(Math2(Math2) (https://youtu.be/A1WHBb3e4PM)

3.Ch3 靜力平衡 (Phy1)(https://youtu.be/hUTu2Iugehs )

4. Ch4 圓周運動與簡諧運動 (Phy2(Phy2)(https://youtu.be/rVqLp5ovMj0)

5. Ch7 圓周運動,向心力,星球重力位能能量轉換(Phy2a) (https://youtu.be/3F3d7UyUf_w )

向心加速度推導(Phy2b) (https://youtu.be/fXQnEgDPeDo ),行星重力位能(Phy2c)(https://youtu.be/LCfYtdjHkyE)

6.Ch8 能量轉換(Phy3(Phy3) (https://youtu.be/9xYXY3BG-YM)

7.Ch9 碰撞(Phy4(Phy4) ( https://youtu.be/1BFmF7ajhsk)

3 of 80

基礎知識

*

3

BuzzNote ( https://tg.pe/buzzBuzzNote ( https://tg.pe/buzz ; 李柏翰老師自建網站)

李柏翰老師 教學影片

8.Ch 10 PV=nRT (Phy5(Phy5) (https://youtu.be/k7I5EgglbNc)

Ch 10 氣體動力論 (Phy5a) (https://youtu.be/6keK5xB5H_Q)

9.波動方程初步 (Phy6(Phy6) ( https://youtu.be/e_i86p-OMd0)

10.干涉與繞射 (Phy7 (Phy7) (https://youtu.be/bKiztPwPYG4 (Phy7) (https://youtu.be/bKiztPwPYG4),雙狹縫干涉(Phy7a) (https://youtu.be/WPmmNne1CZU)。

11. Ch13-1 光波的干涉,Maxwell (Phy8(Phy8) (https://youtu.be/ny0TBjWSAWo )

Ch13-2 光波的干涉,Maxwell (Phy9(Phy9) (https://youtu.be/sUxOD8bJiRQ )

Ch13 光波的干涉波動,視深 (Phy9a) (https://youtu.be/zb4rDF22ASw )

12.微分方程 (Math3 (Math3) (https://youtu.be/qJbA_sVJ2w4 (Math3) (https://youtu.be/qJbA_sVJ2w4) (網站 (Math3) (https://youtu.be/qJbA_sVJ2w4) (網站) (講義)。

13微分方程初步(Math4(Math4)(https://youtu.be/Jxx70ws4WGQ )(ppt)

14 波動光學 (Phy9)(ppt)

16 Ch18 Maxwell equations, 重力位能,磁能轉換 (Phy10(Phy10) (https://youtu.be/qmF41SEkX7)

Ch18 Maxwell equations 推導 (Phy11(Phy11) (https://youtu.be/TXL2oLZczio)

Ch21-4 Bohr氫原子能階 (Phy15) (https://youtu.be/l2Z8p5XnQJs)

4 of 80

*

4

y=f(x)

f(x) 為一函數,

Q(x,f(x)) 是圖形上異於 P 的任意一點﹐

割線

Q 點逐漸向 P 點靠近,

P(a,f(a))是 y=f(x) 圖形上的一個定點,

P(a,f(a))

Q(x,f(x))

切線

1. 斜率與微分

微積分-part 1

5 of 80

x

y

O

x=1

y=f(x)=x2

x

y

O

x=1

y=f(x)=x2

x

y

O

x=1

y=f(x)=x2

函數 y=f(x)=x2 在區間 [0,1]所圍的綠色區域面積為何?

2. 面積與積分

n 條紅色長條面積和

所求綠色面積

n 條藍色長條面積和

6 of 80

一個無窮數列〈an〉,當 n 趨近於無限大時( 記做 n → ∞ )

an 的值能夠趨近一個定值 α

此時稱數列〈an〉為收斂數列,其極限值為 α

n → ∞ 時,

注意:n 越大,

1. 收斂數列:

極限

7 of 80

若數列〈an〉與〈bn〉為收斂數列,

2. 數列極限的基本性質:

則:數列〈an+bn〉、〈anbn〉、〈anbn〉都收斂,

8 of 80

3. 範例 求下列各極限值:

9 of 80

4. 範例

10 of 80

5. 範例

11 of 80

6. 範例

12 of 80

7. 範例

畫出 y=f(x) 與 y=g(x) 的圖形,試問兩者是否相等?

x

O

y

4

2

x

O

y

y=x+2

4

2

解:(1) y=f(x)的圖形如左下圖,且f(x)的定義域為 R

(2) y=g(x)的圖形如右下圖,且g(x)的定義域為{xx≠2} 。

f(x) 與 g(x) 不相等

13 of 80

y=f (x)

x

1.999

3.999

1.99

3.99

1.8

3.8

1.9

3.9

4

2

2.001

4.001

2.01

4.01

2.1

4.1

2.2

4.2

x左側逐漸接近 2

x右側逐漸接近 2

f(x)隨之接近 4

f(x)隨之接近 4

(1) 考慮函數 y=f (x)=x+2,

觀察函數 f (x) 的值的變化。

8. 觀察函數的極限

x

O

y

y=x+2

4

2

稱「當 x 趨近 2 時,

x 愈接近 2 時,

函數 f(x)隨之愈接近 4

函數 f(x)的極限為 4

(但是 x 2)

x 的值愈來愈接近 2

14 of 80

y=f (x)

x

1.999

3.999

1.99

3.99

1.8

3.8

1.9

3.9

2.001

4.001

2.01

4.01

2.1

4.1

2.2

4.2

x左側逐漸接近 2

x右側逐漸接近2

f(x)隨之接近 4

f(x)隨之接近 4

x

O

y

y=x+2

4

2

x 趨近 2 時,

函數 f(x)的極限為 4

,當x 趨近2 時,

2

×

觀察 f (x) 值的變化。

15 of 80

x

O

y

L

a

x

O

y

L

a

y=f(x)

y=f(x)

當函數 f(x) 滿足

,稱函數 f(x) 在 x=a 連續

x

O

y

L

a

y=f(x)

f(x) 在 x=a 連續

f(x) 在 x=a 不連續

9. 函數的連續

說明:

f(x) 在 x=a 不連續

16 of 80

若函數 f(x) 在 x=a 連續

10. 連續函數與極限

(1) 在 x=a 連續:

(2) 在 x=a 不連續:

例:f(x)=x+3x=2 連續

f(x)=x2x=3 連續

f(x)=3x2+4在 x=1 連續

x=a 代入 f(x) 即可

須另外化簡求極限

17 of 80

11. 範例 已知點 P(2,4)在二次函數 f(x)=x2 的圖形上,

y

x

O

y=x2

P(2,4)

切線 L

Q(x,f(x))

求以 P 點為切點的切線方程式。

解:切線 L 的斜率

割線

18 of 80

12. 範例 已知點 P(1,2) 在三次函數 f(x)=x3+1 的圖形上,

y=x3+1

P(1,2)

Q(x,f(x))

切線 L

求以 P 點為切點的切線方程式。

解:切線 L 的斜率

割線

19 of 80

13. 函數的夾擠定理:

若三個函數 fgh 滿足下列兩條件:

(1) 在 a 附近且 x≠a x 值,恆有 f(x) g(x) h(x)

h(x)

f(x)

g(x)

α

說明:

20 of 80

14. 範例:

C

B

A

O

x

1

tanx

1

21 of 80

22 of 80

x 很小時,sinxx

注意:

23 of 80

15. 範例 試利用夾擠定理證明:半徑為 r 之圓的面積=πr2

A

O

B

r

r

證明:(1) 圓內接正 n 邊形面積

= n ∙ ΔAOB面積

24 of 80

P

O

Q

∴ 由夾擠定理得:圓面積=πr2

圓內接正 n 邊形面積 圓面積圓外切正 n 邊形面積

r

H

(2) 圓外切正 n 邊形面積

= n ∙ ΔPOQ面積

25 of 80

若函數 f(x) 在定義域中 x=a

例如: 求函數 f(x)=x2 x=5 處的導數。

1. 導 數 的 定 義

微分

26 of 80

2. 範例

27 of 80

(1) 當函數 f(x) 在 x=a 處的導數 f′(a) 存在時,

稱函數 f(x) 在 x=a 可微分

3. 導 數 與 微 分

(2) 若函數 f(x) 在 x=a 處可微分,

則函數圖形 y=f(x) 在點 P(a, f(a)) 的切線斜率為 f′(a),

否則稱 f(x) 在 x=a 處不可微分。

P(a,f(a))

Q(x,f(x))

切線

y=f(x)

28 of 80

設函數 f(x) 表示一個運動質點在時刻 x 的位移函數,

4. 平均速度與瞬時速度

29 of 80

5. 範例 已知一物體從離地面 50 公尺高處自由落下,

(1)求時刻t=2到t=3之間的平均速度 (2)求時刻t=2的瞬時速度。

30 of 80

(2) 若函數 f(x) 在定義域中的每一個 x=a 處都有導數

如此對應構成一個函數關係,此函數以 f ′(x) 表示之。

我們稱 f ′(x)f(x) 的導函數,並稱 f(x) 為可微分函數。

(1) 若函數 f(x) 在定義域中 x=a

6. 導 函 數

31 of 80

例:函數 f(x)=x3x=a 處的導數

32 of 80

7. 範例

33 of 80

34 of 80

8. 範例

35 of 80

36 of 80

r 是任意實數,且 f(x)=xr

9. 微分公式

則導函數 f ′(x)= rxr−1 。

37 of 80

10. 範例 已知 P(2,−2) 為二次函數 f(x)=x2−2x+2 的圖形外一點,

求通過 P 點的切線方程式。

x

O

y

P(2,− 2)

y=x2−2x+2

解:設切點A(a , a2−2a+2)

過 P 點的切線 L:y−(−2)=(2a−2)(x−2)

A(a,a2−2a+2)

f ′(x)=2x−2

切線 L

a2−2a+2−(−2)=(2a−2)(a−2)

整理得 a2−4a=0

(1)切點(4,10), f ′(4)=6

(2)切點(0,2), f ′(0)= −2

∴ 切線斜率=f ′(a)=2a−2

A在 L上:

a=4 或 0 。

⇒ 切線 L:6xy=14

⇒ 切線 L:2x+y=2

38 of 80

1. 範例 f(x)=x2 的圖形與直線 y=0,x=0 及 x=1 所圍成的

區域 R ,將閉區間 [0,1] 等分成 n 個線段後,

(1) 求上和Un與下和Ln。 (2) 求區域 R 的面積。

x

y

O

x=1

解:(1)上和 Un(藍色長條面積和 )

y=f(x)=x2

積分-part 2

39 of 80

x

y

O

x=1

y=f(x)=x2

下和 Ln (紅色長條面積和 )

40 of 80

x

y

O

x=1

y=f(x)=x2

x

y

O

x=1

y=f(x)=x2

x

y

O

x=1

y=f(x)=x2

n 條紅色長條面積和

所求綠色面積

n 條藍色長條面積和

41 of 80

x

y

O

x=1

y=f(x)=x2

x

y

O

x=1

y=f(x)=x2

2. 積分的符號表示法

總和

說明:

42 of 80

x

y

O

y=f(x)

2

x

y

O

y=f(x)

4

x

y

O

y=f(x)

3

x

y

O

y=f(x)

1

設函數 f(x) 在區間 [0,b] 上連續,且 b>0,f(x)≥0,則:

3、積分與面積的關係

43 of 80

x

y

O

y=f(x)

f(t)

f(c)

t

c

4. 微積分基本定理

若函數 f(x) 在區間 [0,b] 上連續,且 0 < t < b

說明:

44 of 80

積分

微分

5. 積分即為反導數

此時 f(x) 為 g(x) 的導函數。

我們稱 g(x) 為 f(x) 的反導函數

例如:f(x)=x2 的反導函數可為

45 of 80

f(x)=x3 的反導函數為

f(x)=x4 的反導函數為

f(x)=12x2+10x 的反導函數為

說明:

f(x)=x3+2x2−4 的反導函數為

46 of 80

6. 定積分的求法

x

y

O

y=f(x)

f(b)

f(a)

b

a

= 紅色區域面積

說明:

若函數 f(x) 在區間 [a,b] 上連續,

47 of 80

x

y

O

5

f(x)=x3

f(x)

dx

例如:

48 of 80

7. 定積分的性質

若函數 f(x) 與 g(x) 在區間 [a,b] 上連續,則:

49 of 80

說明:

(4)f(x)≥0,

x

y

O

y=f(x)

f(b)

f(a)

b

a

(5)f(x)≤0,

x

y

O

a

b

f(a)

f(b)

y=f(x)

說明:

f(x) 的圖形與直線 y=0,x=ax=b 所圍面積為 R,則:

50 of 80

8. 範例 右下圖為 f(x)=x2−3x+2 的圖形,求灰色區域的面積。

x

O

y

y=f(x)

f(x)

dx

f(x)

dx

1

2

解:在區間 [0,1] 上,f(x)≥0

在區間 [1,2] 上,f(x)≤ 0

51 of 80

9. 範例 求兩拋物線 y=x2 y= −x2+2x+4 所圍成之區域面積。

x

O

y

(−1,1)

(2,4)

y= −x2+2x+4

y=x2

dx

解:所求面積

(−x2+2x+4)x2

52 of 80

設函數 f(x) 在區間 [a,b] 上連續

y=f(x) 與直線 y=0x=a,x=b

所圍成之區域繞 x 軸旋轉,

x

y=f(x)

a

b

x

旋轉半徑R(x)=f(x)

10. 旋轉體的體積

所得旋轉體的體積為

53 of 80

y=f(x)

a

b

x

說明:

薄片

旋轉體

很細的長條形

旋轉一圈

堆 疊

f(x)

π (f(x))2

x

54 of 80

11. 範例 證明:直圓錐體的體積等於其底面積與高

之乘積的三分之一。

解:所求體積

x

y

O

x

(h,r)

r

h

55 of 80

(r,0)

x

y

(0,r)

O

(0,−r)

(−r,0)

x

解:所求體積

12. 範例

56 of 80

Sec2: 微分方程初步

物理學與動力學?

🡺 連結常常是微分方程(Differential Equations)

什麼是微分方程?

🡺一些變數微分後與若干變數導數的函數組合,所形成的方程式,所有的變數導數綁定在此方程式中,不再彼此獨立無關。

*

56

(牛頓運動定律F = ma ,時間為變數,F = -k x )

57 of 80

微分方程定義

微分方程

(1)常微分方程(ODE,Ordinary Differential Equations),只含一個獨立變數以及它的常導數,則稱為常微分方程。

(2)偏微分方程(PDE,Partial Differential Equations)

含有兩個或兩個以上獨立變數以及它的偏導數,則可定義此為偏微分方程。

Schrodinger equations 就是偏微分方程,變數 t,x

*

57

(free-particle)

58 of 80

微分方程 order & degree

微分方程

🡺階(order,微分次數)

🡺次(degree,微分項幾次方)

ODEPDE 都有線性非線性之分。

二階一次

線性常微分方程

*

58

59 of 80

(1)一階一次

線性常微分方程

(2)二階一次

線性常微分方程

(3)一階一次

非線性常微分方程

(4)一階二次

非線性常微分方程

(5)二階一次

線性偏微分方程

*

59

幾階

幾次

線性?

60 of 80

微分方程的解

*

60

微分方程求解

  • 若求得解為函數y,而cy 仍是其解,c 是常數,可判斷為齊次微分方程(homogeneous DE)。

  • n 階微分方程,可以得到n 個任意常數解的通解。

  • 如果再加入邊界條件,或初始值的判斷,可以得到特解,如此特定的解不是通解,而是異解(singular solution)。

61 of 80

微分方程 Q1

Answer:

  1. (分離變數)

如果ODE 能夠將變數分隔開來,用積分處理,問題求解將會是非常容易。比如說此題,等號分隔開變數與常數。

*

61

Q1 :

62 of 80

微分方程Q1求解

(2). 兩邊同時積分

*

62

得到通解 c1,c2 為常數。

63 of 80

微分方程求解

(3). 若有起始條件

*

63

特解完成

64 of 80

微分方程 Q2

Answer:

*

64

Q2 :彈簧拉動物體運動,F = ma,若恢復力為彈簧彈力,

F = -kx, 如何列出ODE,k > 0?

ODE完成

65 of 80

微分方程 Q3

*

65

Q3 :彈簧拉動物體運動,F = ma,若恢復力為彈簧彈力,

F = -kx, 受到阻力 f= -bv,如何解,k > 0?

66 of 80

微分方程Q3求解

Answer:

*

66

 

67 of 80

微分方程Q3求解

*

67

 

Answer:

68 of 80

微分方程Q3求解

*

68

討論: r和w值對震盪圖形的關係

 

(Fig.2: 以r=2,w=1為例)

69 of 80

微分方程Q3求解

*

69

討論: r和w值對震盪圖形的關係

 

(Fig.3: 以r=2,w=2為例)

70 of 80

微分方程Q3求解

*

70

討論: r和w值對震盪圖形的關係

 

Fig.4(a):X-t,以r=1,w=2。

Fig.4(b): V-t 圖

71 of 80

微分方程Q3求解

*

71

討論: r和w值對震盪圖形的關係

Fig.5(b): V-t 圖

Fig.5(a): X-t,以r=0,w=2。

 

72 of 80

微分方程Q3求解

*

72

討論: r和w值對震盪圖形的關係

Fig.5(c): x-t,V-t,V-x 相圖

 

73 of 80

微分方程Q3求解

*

73

討論: r和w值對震盪圖形的關係

 

Fig.6: x-t,V-t,V-x 相圖

74 of 80

微分方程Q3求解

*

74

Fig.7: x-t 振盪相圖

 

 

(6)整理上述情況,(a) r>w,overdamping(b) r =w critically damped(c) r<w,underdamping

x 是位移,v 是速度,b= 2r = 阻尼

w = 振盪頻率

改寫成兩條一階一次微分方程

75 of 80

一階常微分方程

一階常微分方程在日常生活中有很多運用,常見的表示方式有:

*

75

又或者是表示為:

76 of 80

一階常微分方程

一階常微分方程如果可以寫成:

*

76

兩邊同時積分,則可以得到:

77 of 80

微分方程Q4求解

*

77

Q4:

Answer:

  1. (分離變數)

目前變數x,y 可以分隔開來,如果能夠將變數分隔開來,用積分處理,問題求解將會是非常容易。

78 of 80

微分方程Q4求解

*

78

Answer:

(2) 利用連鎖律

兩邊同時積分

79 of 80

微分方程Q4求解

*

79

(3) 對照連鎖律

80 of 80

微分方程Q4求解

*

80

(4) 整體解

k 為常數,DE 詳解完成。

這便是用連鎖律變數分隔處理,這個連鎖律的使用是

一階微分方程解常見的技巧。Sec1完!