• 沒有找到結果。

提要 33:認識 Euler-Cauchy 方程式的解法(三)--複數根

N/A
N/A
Protected

Academic year: 2021

Share "提要 33:認識 Euler-Cauchy 方程式的解法(三)--複數根"

Copied!
6
0
0

加載中.... (立即查看全文)

全文

(1)

提要 33:認識 Euler-Cauchy 方程式的解法(三)--複數根

Euler-Cauchy 方程式係定義為:

2 0

2 2 + +by=

dx axdy dx

y

x d (1)

上式係變係數之齊性常微分方程式,然而只要作適當之變數變換,即可發覺其在本質上 亦為常係數之齊性常微分方程式,說明如下。首先,令:

ez

x= (2) 故

x

z ln= (3) 而







 

 −

=



 

=



 

 +



 

=

 

= 

=

=

=

2 2

2 2 2

2 1 1 1

1 1

1 1

x dz dy dz

y d x dz dy x dx dz dz

dz d dy

dx d x dz dy x dx

dz d dy

x dz dy dx

d dx dy dx

d dx

y d

x dz dy dx dz dz dy dx dy

故式(1)可改寫為:

(

1

)

0

2

2 + − +by=

dz a dy dz

y

d (4)

上式即為常係數之齊性常微分方程式,其解析方法已於之前介紹過,今再簡述如下。

式(4)係自然界中之問題的化身,大自然的問題之解應與大自然的函數有關,此一 問題所牽涉之大自然的函數為自然指數函數(Exponential Function),因此可考慮問題之 解為:

( )

z

( )

z e z

y = expλ = λ (5) 其中λ 為待解之未知數,稱為特徵根(Characteristic Root)。再將式(5)代回式(4):

( )

2

(

1

) ( ) ( )

0

2 z + − z +be z =

dz e a d

dz e

d λ λ λ

(2)

( )

[

λ2 + a1λ+b

]

eλz =0

因為eλz表所假設的解y,而我們不希望所研討出之解為零,亦即eλz ≠0,故:

(

1

)

0

2 + a− λ+b=

λ (6) 上式係由特徵根所構成之方程式,稱為特徵方程式(Characteristic Equation)。特徵方程 式係一元二次方程式,其解可利用因式分解法或以下所示之公式解法研討出:

( ) ( )

2

4 1

1 a 2 b

a− ± − −

= −

λ (7)

因所考慮之情況為複數根,故上式中根號內之運算結果應為

(

a−1

)

2 −4b<0。可令所得 出之兩個特徵根分別為λ 與1 λ : 2

( ) ( )

2

1 4

1 2

1

− +

= − a i b a

λ 、

( ) ( )

2

1 4

1 2

2

= − a i b a

λ (8)

故所研討出之問題的解分別為:

( )

z e z

y = λ1y

( )

z =eλ2z (9) 或

( )

e λ1 xλ1

y= z = 、y =

( )

ez λ2 = xλ2 (9’) 由重疊原理(Superposition Principle)知,線齊性微分方程式之解可以疊加方式加以表示:

( )

x C1xλ1 C2xλ2

y = + (10) 式(10)稱為問題之通解(General Solution)。

若考慮

2

−1

= a

α 、

( )

2 1 4 − − 2

= b a

β ,並引用尤拉公式(Euler Formula):

x i x

e±ix =cos ± sin (11) 則式(10)可進一步改寫為:

( )

x C xα iβ C xα iβ

y = 1 + + 2 (10’) 亦即:

(3)

( ) ( )

( ) ( )

[ ]

( ) ( )

[ ]

( ) ( )

[ ] [ ( ) ( ) ]

{ }

( ) ( ) ( ) ( )

[

C C x i C C x

]

x

x i

x C

x i

x C

x

x i C

x i C x

x C

x C

x

x C x C x x y

i i

i i

ln sin ln

cos

ln sin ln

cos ln

sin ln

cos

ln exp

ln exp

ln exp ln

exp

2 1 2

1

2 1

2 1

2 1

2 1

β β

β β

β β

β β

α α α

β β

α

β β

α

− + +

=

− +

+

=

− +

=

+

=

+

=

( ) ( )

[

A x B x

]

x

y= α cos βln + sin βln (12) 其中A=C1 +C2B=i

(

C1C2

)

另解

式(1)之解析,亦可直接令:

xm

y= (13) 再代回式(1),則

( )

2

( ) ( )

0

2 2 m + m +b xm =

dx x axd dx

x x d

( )

[

1 2

] (

1

) ( )

0

2 m mxm +axmxm +b xm = x

上式可進一步改寫為:

( )

[

m m1 +am+b

]

xm =0 (14) 因為y =xm ≠0,所以

(

1

)

0

2 + am+b=

m (15) 故

( ) ( )

2

4 1

1 2

1

b a

ma− + − −

= 、

( ) ( )

2

4 1

1 2

2

b a

ma− − − −

= (16)

y =xm1y =xm2均為問題之解,由重疊原理知,問題之通解可表為:

(4)

( )

x C1xm1 C2xm2

y = + (17)

因考慮複數根,即

(

a−1

)

2 −4b<0,若假設

2

−1

= a

α 、

( )

2 1 4 − − 2

= b a

β ,則:

β α i

m1 = + 、m2 =α −iβ (18) 故式(17)亦可表為:

( ) ( )

( ) ( )

[ ]

( ) ( )

[ ]

( ) ( )

[ ] [ ( ) ( ) ]

{ }

( ) ( ) ( ) ( )

[

C C x i C C x

]

x

x i

x C

x i

x C

x

x i C

x i C x

x C

x C

x

x C x C x

x C x

C x y

i i

i i

i i

ln sin ln

cos

ln sin ln

cos ln

sin ln

cos

ln exp

ln exp

ln exp ln

exp

2 1 2

1

2 1

2 1

2 1

2 1

2 1

β β

β β

β β

β β

α α α

β β

α

β β

α

β α β

α

− + +

=

− +

+

=

− +

=

+

=

+

=

+

=

+

( ) ( )

[

A x B x

]

x

y= α cos βln + sin βln (19) 其中A=C1 +C2B=i

(

C1C2

)

附註:

( )

xm =mxm1

( )

m ″ =

(

−1

)

m2

x m m

x

(5)

範例一

試解出Euler-Cauchy 方程式之解:x2y′′+xy′+y =0。

【解答】

y

( )

x =xm,代回原式,則

( ) ( ) ( )

0

2 m ″ + m ′+ m =

x x

x x

x

上式可化簡為:

( )

[

mm−1 +m+1

]

xm =0

因為不希望所研討出之解為零,亦即xm ≠0,故:

0

2 +1= m 上式可因式分解為:

(

m+i

)(

mi

)

=0 故問題之解的特徵根(Characteristic Root)m為:

i

λ = 、λ =i

由重疊原理(Superposition Principle)知,問題之通解(General Solution)為:

( ) ( ) ( )

( ) ( )

( ) ( )

[ ] [ ( ) ( ) ]

( ) ( ) ( ) ( )

( )

x B

( )

x A

x C

C i x C

C

x i

x C

x i

x C

x i C x i C

x C

x C

x C x C x y

i i

i i

ln sin ln

cos

ln sin ln

cos

ln sin ln

cos ln

sin ln

cos

ln exp ln

exp

ln exp ln

exp

2 1 2

1

2 1

2 1

2 1

2 1

+

=

− + +

=

+ +

=

+

=

+

=

+

=

(6)

習題 1. Please solve the following initial value problem.

0

2y′′+xy′+4y=

x with y

( )

1 =2, y′

( )

1 = −1 and x>0.【95 暨大研究所】

2. Solve

(

3x4

)

2 y′′+3 3

(

x4

)

y+36y= .【92 雲科電機所 15%】 0 3. Solve

(

x2 2x 1

)

d y22

(

5x 5

)

dy 5y 0

dx dx

+ + + + + = .【91 北科通訊所 15%】

參考文獻

相關文件

在這一節裡會提到,即使沒辦法解得實際的解函數,我們也 可以利用方程式藉由圖形(方向場)或者數值上的計算(歐拉法) 來得到逼近的解。..

如果函數是由基本函數所組成,至少需要注意:分式函 數分母會等於 0

 MATLAB 程式使用 pass-by-value 的方 式,進行程式與函式間的溝通聯絡,當 程式呼叫函式時, MATLAB

This kind of algorithm has also been a powerful tool for solving many other optimization problems, including symmetric cone complementarity problems [15, 16, 20–22], symmetric

[r]

學習範疇 主要學習成果 級別 級別描述 學習成果. 根據學生的認知發展、學習模式及科本

For R-K methods, the relationship between the number of (function) evaluations per step and the order of LTE is shown in the following

To improve the convergence of difference methods, one way is selected difference-equations in such that their local truncation errors are O(h p ) for as large a value of p as