• 検索結果がありません。

微 分 方 程 式

N/A
N/A
Protected

Academic year: 2021

シェア "微 分 方 程 式"

Copied!
14
0
0

読み込み中.... (全文を見る)

全文

(1)

微 分 方 程 式

8. 非同次微分方程式

(2)

1. 非同次微分方程式 –1–

非同次線形微分方程式

y

′′

+ P (x)y

+ Q(x)y = R(x) ( )

の一般解を求める.

定理.

(i) ( )

1

つの特殊解

u(x)

および同次微分方程式

y

′′

+ P (x)y

+ Q(x)y = 0 ( ∗∗ )

2

つの

1

次独立な解

y

1

(x), y

2

(x)

が求められたとき

, ( )

の一般 解は

y = u(x) + c

1

y

1

(x) + c

2

y

2

(x)

である.

(ii)

同次微分方程式

( ∗∗ )

1

次独立な

2

つの解

y

1

(x), y

2

(x)

が求 まったとき,

( )

の一般解ははつぎのように与えられる:

y = c

1

y

1

(x)+ c

2

y

2

(x)+ y

1

(x)

R(x)y

2

W (y

1

, y

2

) dx + y

2

(x)

R(x)y

1

W (y

1

, y

2

) dx

(3)

[証明]

(i) y = z + u(x)

とおくと

–2–

L(y) = L(z + u) = L(z) + L(u) = R(x)

L(u) = R

だから

L(z) = 0. L(z) = 0

の一般解は

c

1

y

1

(x) + c

2

y

2

(x)

だから

( )

の一般解は

y = u(x) + c

1

y

1

(x) + c

2

y

2

(x).

ii

[

定数変化法

] y = c

1

(x)y

1

(x) + c

2

(x)y

2

(x)

とおくと

y

= c

1

(x)y

1

+ c

2

(x)y

2

+ c

1

(x)y

1

+ c

2

(x)y

2

.

ここで

c

1

(x)y

1

+ c

2

(x)y

2

= 0

とおく.このとき

y

= c

1

(x)y

1

+ c

2

(x)y

2 となるから微分して

y

′′

= c

1

(x)y

′′1

+ c

2

(x)y

2′′

+ c

1

(x)y

1

+ c

2

(x)y

2 よって

L(y) = c

1

(x)L(y

1

) + c

2

(x)L(y

2

) + c

1

(x)y

1

+ c

2

(x)y

2

= R(x).

(4)

L(y

1

) = 0, L(y

2

) = 0

だから

–3–

c

1

(x)y

1

+ c

2

(x)y

2

= R(x)

ここで

c

1

(x)y

1

+ c

2

(x)y

2

= 0 c

1

(x)y

1

+ c

2

(x)y

2

= R(x)

c

1

(x), c

2

(x)

の連立方程式と思うと

W (y

1

, y

2

) ̸ = 0

より

c

1

(x) = R(x)y

2

W (y

1

, y

2

) c

2

(x) = R(x)y

1

W (y

1

, y

2

) x

について積分すると

c

1

(x) =

R(x)y

2

W (y

1

, y

2

) dx + c

1

, c

2

(x) =

R(x)y

1

W (y

1

, y

2

) dx + c

2

. (c

1

, c

2 は定数).したがって,一般解は

y = c

1

y

1

(x)+ c

2

y

2

(x)+ y

1

(x)

R(x)y

2

W (y

1

, y

2

) dx + y

2

(x)

R(x)y

1

W (y

1

, y

2

) dx

(5)

例題

1

–4–

同次微分方程式

( ∗∗ )

1

つの解

y

1

(x) ̸ = 0

が求まったとき,

(i) y

1

∫ ( e

P dx

y

12

Ry

1

e

P dx

dx )

dx

は非同次微分方程式

( )

1

つの特殊解である

.

(ii) y

1

e

P dx

y

12

dx

y

1

(x)

1

次独立な

( ∗∗ )

の解である

(iii)

非同次微分方程式

( )

の一般解は

c

1

y

1

+ c

2

y

1

e

P dx

y

12

dx + y

1

∫ ( e

P dx

y

12

Ry

1

e

P dx

dx )

dx.

同次方程式の解がひとつわかれば非同次方程式は解ける!

(6)

[証明]

(i) y = uy

1

( )

へ代入すると

–5–

(y

1′′

+ P y

1

+ Qy

1

)u + y

1

u

′′

+ (2y

1

+ P y

1

)u

= R.

y

1′′

+ P y

1

+ Qy

1

= 0

だから

y

1

u

′′

+ (2y

1

+ P y

1

)u

= R.

u

= v

とおくと

y

1

v

+ (2y

1

+ P y

1

)v = R.

これは

1

階線形方程式であり,一般解は

v = e

∫ 2y

1

+ P y

1

y

1

dx

 

R y

1

e

2y

1

+ P y

1

y

1 dx

dx + c

1

  (1)

= e

P dx y

12

(∫

Ry

1

e

P dx

dx + c

1

)

. (2)

(7)

さらに積分すると

–6–

u =

e

P dx y

12

(∫

Ry

1

e

P dx

dx + c

1

)

dx + c

2

. (3)

ここで

c

1

= c

2

= 0

として

y = y

1

u

とすれば

(i)

を得る

(8)

(ii) (3)

R 0, c

1

= 1, c

2

= 0

とすれば解であることはわかる.

1

次独立であることを証明する. 任意の定数

c

1

, c

2に対して

c

1

y

1

+ c

2

y

1

e

P dx

y

12

dx = 0

が成立したとする.

y

1 で割ると

c

1

+ c

2

e

P dx

y

12

dx = 0

さらに微分すると

c

2

e

P dx y

12

= 0

したがって,

c

2

= 0

となりさらに

c

1

= 0

にもなる.

(9)

(iii)

定理

1

より明らか.

–7–

1

かっこ内に与えられた関数が

1

つの特殊解であることを知ってつ ぎの微分方程式を解け.

(1) (1 + x)y

′′

+ xy

y = 0 (y = x).

(2) xy

′′

y

+ (1 x)y = 0 (y = e

x

).

2

かっこ内の関数が対応する同次方程式の解であることを知ってつ ぎの微分方程式を解け.

(1) y

′′

x

x 1 y

+ 1

x 1 y = x 1 (y = e

x

).

(2) (x + 1)y

′′

(2x + 3)y

+ 2y = xe

x

(y = e

2x

).

(10)

定数係数線形常微分方程式 –9–

2

階線形微分方程式

y

′′

+ ay

+ by = 0 (1)

において

a, b

は定数とする.

y = e

λx を代入すると

eλx(λ

2

+ + b) = 0

となる

.

2

次方程式

λ

2

+ + b = 0

を特性方程式といい、その解

λ

を特性根という

λ

が特性根のとき

, y = e

λxは解になっている.

(11)

定理

1

特性根を

λ

1

, λ

2 とする,

–10–

(1) λ

1

, λ

2 が相異なる実根のとき.

e

λ1x

, e

λ2x

1

次独立な(

1

)の 解であり,一般解は

c

1

e

λ1x

+ c

2

e

λ2x である

.

(2) λ

1

= λ

2 のとき.

e

λ1x

, xe

λ2x

1

次独立な

(1

)の解であり,一 般解は

(c

1

+ +c

2

x)e

λ1xである.

(3) λ

1

, λ

2 が共役な複素根であるとき,

λ

1

= α +

α, β)

は実数

β ̸ = 0)

ならば,

e

αx

cos βx, e

αx

sin βx

1

次独立な(

1

) の解で あり,一般解は

e

αx

(c

1

cos βx + c

2

sin βx)

である.

[証明]

(1

) 特性根の性質から

e

λ1x

, e

λ2x が解であることは明 らか.

ロンスキアンをつくると

W (e

λ1x

, e

λ2x

) =

e

λ1x

e

λ2x

λ

1

e

λ1x

λ

2

e

λ2x

= (λ

2

λ

1

)e

12)x

̸ = 0

であるから

1

次独立 である.

(12)

ii) e

λ1xが解であることは明らか。

等根であることから

1

= a

を用いると

L(xe

λ1x

) = xe

λ1x

21

+

1

+ b) + e

λ1x

(2λ

1

+ a) = 0.

ロンスキアンをつくると,

W (e

λ1x

, xe

λ1x

) =

e

λ1x

xe

λ1x

λ

1

e

λ1x

(1 + λ

1

x)e

λ1x

= e

1x

̸ = 0

だから

1

次独立である.

(13)

(iii)

解と係数の関係により 

–12–

λ

1

+ λ

2

= 2α = a, λ

1

· λ

2

= α

2

+ β

2

= b.

より

2α + a = 0, α

2

β

2

+ b = 0 e

αx

cos βx, e

αx

sin βx

(1)

に代入して

L(eαxcosβx) = eαxcosβx(α2 −β2 −aα+b)−eαxsinβx(2α+ a) = 0 L(eαxsinβx) = eαxsinβx(α2 −β2 −aα +b) +eαxcosβx(2α+a) = 0

ロンスキアンをつくると

W(eαxcosβx, eαxsinβx) =

eαxcosβx eαxsinβx

eαx(αcosβx−βsinβx) eαx(αsinβx+βcosβx)

= βeαx ̸= 0

だから

1

次独立である.

(14)

定数係数非同次線形常微分方程式 –13–

非同次方程式

y

′′

+ ay

+ by = R(x)

の一般解については

,

同次微分方程式の解が容易にわかるので定 数変化法を用いればよい.

参照

関連したドキュメント

力学では運動方程式が代表例で、ma=fのaは速度の

変数分離形微分方程式(separation of

自ら問題を解くこ とで初めて微分方程式を解くことができるようになるので , できるだけ多くの問題を 解いてもらいたい...

[r]

[r]

[r]

[r]

[r]