• 沒有找到結果。

常系数非齐次线性微分方程 第九节

N/A
N/A
Protected

Academic year: 2021

Share "常系数非齐次线性微分方程 第九节"

Copied!
16
0
0

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

全文

(1)

常系数非齐次线性微分方程 第九节

型 )

( )

(

x e P x f

x m

x x

P e

x

f

( )  x[ l ( ) cos

型 ]

sin )

~ (

x x

P

n

 一、

二、

(2)

) (x f y

q y

p

y      ( p , q 为常数 )

二阶常系数线性非齐次微分方程 :

根据解的结构定理 , 其通解为

Y

y   y *

非齐次方程特解 齐次方程通解

求特解的方法

根据

f (x)

的特殊形式

,

*

给出特解 y

的待定形式 , 代入原方程比较两端表达式以确定待定系数 .

待定系数法

(3)

) ( [ Q x

e

x

  ( 2   p ) Q  ( x )  ( 

2

p   q ) Q ( x ) ] )

(x P

e

x m

一、

f ( x ) e

x

P

m

( x )

 为实数 ,

) (x P

m

设特解为

y *  e

x

Q ( x ) ,

其中 为待定多项式 ,

) (x Q

] ) ( )

( [

* e Q x Q x

y  

x

  

] ) ( )

( 2

) ( [

* e

2

Q x Q x Q x

y  

x

     

代入原方程 , 得

)

(x Q 

(1) 若  不是特征方程的根 ,

,

2

p   q  0

即 

则取

), (x

Q

m 从而得到特解 形式为

y *  e

x

Q

m

( x ) .

) ( )

2

(  p Qx

  (

2

p q ) Q ( x )P

m

(x )

为 m 次多项式

.

Q (x)

m

次待定系数多项

(4)

(2) 若 是特征方程的单根

,

2

p   q  0 , 2 p 0 , )

(x

Q

为 m 次多项式 ,故特解形式为

y *  x Q

m

( x ) e

x (3) 若  是特征方程的重根

,

2

p   q  0 , 2   p  0 , )

(x

Q 

是 m 次多项式 ,

故特解形式为

y *  x

2

Q

m

( x ) e

x 小结 对方程① ,

) 2 , 1 , 0 (

) (

*  x Q x e k

y

k m x

此结论可推广到高阶常系数线性微分方程 .

)

(x

Q   ( 2   p ) Q  ( x )  ( 

2

p   q ) Q ( x ) P

m

(x )

当 是特征方程的

k

重根

时 , 可设

特解

(5)

例 1.

求方程 y   2 y   3 y  3 x  1

的一个特解 . 解 : 本题 而特征方程为

r

2

 r 2  3  0 ,

不是特征方程的根

设所求特解为.

y *  b

0

xb

1

,

代入方程 :

1

3 2

3

3

0

1

0

 

b x b b x

比较系数 , 得

 

 b 3

0

 3

1 3

2

0

1

b b 3

, 1 1

1

0

  bb

于是所求特解为

.

3

*  x   1 y

 0

,

 0

(6)

例 2.

求方程 y   5 y   6 yx e

2x的通解 .

解 : 本题 特征方程为

r

2

 r 5  6  0 ,

其根为

对应齐次方程的通解为

YC

1

e

2x

C

2

e

3x 设非齐次方程特解为

y *  x ( b

0

xb

1

) e

2 x

比较系数 , 得

 

 b 2

0

 1 0

2 b

0

 b

1

 , 1

2 1

1

0

  b   b

因此特解为

y *  x ( 

12

x  1 ) e

2x

. 3

,

2

2

1

r

r

代入方程得

 2 b

0

xb

1

 2 b

0

x

所求通解为

yC

1

e

2x

C

2

e

3x

 (

12

x

2

x ) e

2x

. ,

 2

(7)

二、

f ( x )  e

x

P

l

( x ) cos  xP ~

n

( x ) sin  x  型

P

m

x e

i x

x

f ( ) ( )

( )

P

m

( x ) e

(i)x 第二步 求出如下两个方程的特解

x m

x e

i

P y

q y

p

y      ( )

( )

 

  p y q y y

分析思路 :

第一步 将 f (x) 转化为

第三步 利用叠加原理求出原方程的特解 第四步 分析原方程特解的特点

x m

x e

i

P ( )

( )

(8)

第一步 利用欧拉公式将 f (x) 变形

 

e

x

x

f ( )

 

  

i

x P

x

P

l n

2 )

~ ( 2

)

(

i x

e

( )

 

  

i

x P

x

P

l n

2 )

~ ( 2

)

(

i x

e

( )

P

m

x e

i x

x

f ( ) ( )

( )

P

m

( x ) e

(i) x

P

m

( x ) e

(i) x

P

m

( x ) e

(i)x

 

m  max n , l , ) (x P

l

2

x i x

i

e

e

)

~ x ( P

n

 

 

i e

e

i x i x

2

(9)

第二步 求如下两方程的特解

  i

是特征方程的

k

重根 ( k = 0, 1)

, i x

k

Q

m

x e x

y

1

 ( )

( )

( Q

m

( x ) 为m 次多项式 )

( y

1

)   p ( y

1

)   q y

1

P

m

( x ) e

(i) x

等式两边取共轭 :

x m

x e

i

P y

q y

p

y

1 1

 

1

 ( )

( )

 

1

这说明 y

为方程 ③ 的特解 .

x m

x e

i

P y

q y

p

y      ( )

( )

x m

x e

i

P y

q y

p

y      ( )

( )

则 ②

特解 :

(10)

第三步 求原方程的特解

利用第二步的结果 , 根据叠加原理 , 原方程有特 解 :

y *  y

1

y

1

x k

e x

Q

m

e

i x

Q

m

e

i x

原方程

y   p y   q ye

x

P

l

( x ) cos x P ~

n

( x ) sin x

x k

e x

Q

m

(cos  xi sin  x )

) sin

(cos x i x

Q

m

  

x k

e x

R

m

cos  xR ~

m

sin  x

m m

R R , ~

其中

均为

m

次多项式 .

(11)

第四步 分析

y

的特点

R x R x

e x

y y

y

m x m

k

cos  ~ sin 

1 1

y

1

y

1

*

y y

所以 因此 R

m

, R ~

m 均为

m

多项式 .

y

1

y

1

y

本质上为实函数 ,

y

1

y

1

(12)

小 结 :

P x x P x x

e

x l

( ) cos   ~

n

( ) sin 

对非齐次方程

y q y

p

y    

) ,

( p q 为常数

R x R x

e x

y * 

k x m

cos   ~

m

sin 

则可设特解 :

其中

为特征方程的

k

重根 ( k = 0, 1) ,

  i

n l

m  max ,

上述结论也可推广到高阶方程的情形 .

(13)

例 4.

求方程 y   yx cos 2 x

的一个特解 . 解 : 本题

特征方程

, 2 ,

0 

 

故设特解为

x d

x c x

b x

a

y *  (  ) cos 2  (  ) sin 2

不是特征方程的根 ,

i i   2

 

代入方程得

x x

x a

d x

c x

c b

x

a 3 4 ) cos 2 ( 3 3 4 ) sin 2 cos 2 3

(       

0

2

 1  r

, )

( x x

P

l

P ~

n

( x )  0 ,

比较系数 , 得 9

4 3

1

, 

a

d

. 2 sin 2

cos

*

31

x x

49

x

y

于是求得一个特解

1 3 

 a

0 4

3  

b c

0 3 

 c 3  4  0

d a b  c  0

(14)

例 5.

求方程 y   9 y  18 cos 3 x  30 sin 3 x

的通解 . 解 : 特征方程为

r

2

 9  0 ,

其根为

对应齐次方程的通解为

YC

1

cos 3 xC

2

sin 3 x

) 3 sin 3

cos (

* x a x b x

y  

比较系数 , 得

a  5 , b  3 ,

因此特解为

y *  x ( 5 cos 3 x  3 sin 3 x )

i r

1,2

  3

代入方程 :

6 b cos 3 x  6 a sin 3 x

所求通解为

x C

x C

y

1

cos 3 

2

sin 3

为特征方程的单根 ,

i

 3

) 3 sin 3

3 cos 5

( x x

x

x x 30 sin 3 3

cos

18 

因此设非齐次方程特解为

(15)

内容小结

m

x e

x

P y

q y

p

y ( )

.

1     

为特征方程的 k ( = 0, 1, 2) 重根 ,

m x

k

Q x e x

y *  ( )

则设特解为

] sin

)

~ ( cos

) ( [

.

2 y   p y   q ye

x

P

l

xxP

n

xx

为特征方程的 k ( = 0, 1 ) 重

根 ,

  i

x k

e x

y * 

则设特解为

]

sin )

~ ( cos

) (

[ R

m

xxR

m

xx

l n

m  max ,

3. 上述结论也可推广到高阶方程的情形 .

(16)

思考与练习

时可设特解为

x

x x

f ( ) cos )

1 当 

e

x

x x

x

f ( ) cos 2

2

)

2 当  

x

y *  ( a xb ) cos x

*

y ( a xb ) cos 2 x  ( c xd ) sin 2 xk e

2x

)

(x f y

y   

时可设特解为

P x x P x x

e x

f ( ) 

x l

( ) cos   ~

n

( ) sin 

x k

e x

y * 

n l

m  max ,

提示 :

x d

cx ) sin ( 

1 . (

填空 )

] sin

)

~ ( cos

) (

[ R

m

xxR

m

xx

參考文獻

相關文件

线性拟合与二次拟合 数据拟合的线性模型 一次多项式拟合公式..

超定方程组QR分解算法 数据拟合确定常微分方程..

酸鹼滴定、一次微分、二次微分曲線圖..

对于二元函数的自变量,我

笛卡儿企图通过坐标系给几何引进新方 法, 他的成就远远超出他的期望. 坐标系是数 学中的双刃剑, 使得几何的目的可以通过代 数达到, 反过来,

为此, 我们需要建立函 数的差商与函数的导数间的基本关系式, 这些关系式称为“微分学中值定理”...

Sample a point in the plane, iterates Newton's Method on f(z) using that point as the initial value, and colors the point based on which of the four roots the iterates converge

• 孤波解的性质:1- x 趋向无穷时 u 趋于0;2- 波是单向传播的;3- 波幅与波速成正比;4- KdV方程的色散与非线性刚好抵消,导致