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

2階の定数係数線形常微分方程式

N/A
N/A
Protected

Academic year: 2021

シェア "2階の定数係数線形常微分方程式"

Copied!
4
0
0

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

全文

(1)

Revised at 18:33, January 31, 2017

数学特論B 第13回

http://my.reset.jp/˜gok/math/advanced/ 1

2階の定数係数線形常微分方程式

1

同次方程式

1.1

特性方程式による一般解の分類

一般に微分方程式:

y 00 (x) + ay 0 (x) + by(x) = 0

a, b

は実数)

の一般解は、その特性2次方程式

w 2 + aw + b = 0

の2つの解の様子で分類すれば次の いずれかに当てはまります(

C 1 , C 2

は任意の定数)

特性方程式の解 微分方程式の一般解 相異なる実数解

w 1 , w 2 y(x) = C 1 e w

1

x + C 2 e w

2

x

実重解

w y(x) = C 1 xe wx + C 2 e wx

共役な2つの複素数解

r ± iq y(x) = C 1 e rx cos qx + C 2 e rx sin qx

1.2

基本的な計算と微分作用素による記述

1階の微分方程式:

f 0 (x) + pf(x) = 0

の両辺に

e px

を掛けると

0 = f 0 (x)e px + pf(x)e px = { f (x)e px } 0

と云う風に積の微分法と見る事が出来、微分が

0

となるのは定数関数でしたから

f (x)e px = C

従って、

f (x) = Ce px

が分かります(

C

は任意の定数)。

これを更に別の見方で見てみましょう。

今やった計算は、要するに

f (x)e px = g(x)

と置けば

g(x)

が定数関数であることが分 かると云うわけなんですが、これはつまり

f (x) = g(x)e px

と置くと云う事であって、

これを微分方程式の左辺に代入してみると

f 0 (x) + pf (x) = ©

g(x)e px0

+ pg(x)e px = g 0 (x)e px

となっています。

一般に関数を食べて関数を吐き出す 関数 の事を作用素と呼ぶ事にすれば、関数を 微分すると云う作用素

dx d

と関数を定数倍すると云う作用素の和を使って

µ d dx + p

h(x) = d

dx h(x) + ph(x) = h 0 (x) + ph(x)

と書く事が出来ますが、更に

D = dx d + p

と書くことにすればさっきの計算は

Df (x) = g 0 (x)e px

と書く事が出来ます。これを繰り返せば次が成り立つことが分かります:

事実

1.1 D = dx d + p

と置けば、以下が成り立ちます:

D n ©

g(x)e px

= g (n) (x)e px .

1.3

相異なる実数解の場合の具体例

1.3.1

因数分解法 〜1階の微分方程式を2回解く事に帰着する〜

微分作用素による記述を採用すれば微分方程式:

f 00 (x) 5f 0 (x) + 6f (x) = 0

はあたかも左辺を因数分解するかの様に

µ d dx 3

∂ Ωµ d dx 2

f (x)

æ

= 0

と変形する事が出来ます。ここで

° d

dx 2 ¢

f (x) = h(x)

と置けば、方程式は

µ d

dx 3

h(x) = 0

と書く事が出来、これは簡単に一般解

h(x) = C 1 e 3x

が求まります(

C 1

は任意の定数)。

そこで元の

f (x)

に戻してやれば

µ d

dx 2

f(x) = C 1 e 3x · · · ( )

(2)

Revised at 18:33, January 31, 2017

数学特論B 第13回

http://my.reset.jp/˜gok/math/advanced/ 2

ですからまたここで

f (x) = g(x)e 2x

と置けば、

g 0 (x)e 2x = C 1 e 3x g 0 (x) = C 1 e x

g(x) = C 1 e x + C 2

f (x) = C 1 e 3x + C 2 e 2x

C 1 , C 2

は任意の定数)

となって一般解が求まります。あるいは

g(x)

は導入せずに直接

( )

式を書き直して

f 0 (x) 2f(x) = C 1 e 3x

f 0 (x)e 2x 2f (x)e 2x = C 1 e x

° f (x)e 2x ¢ 0

= C 1 e x f (x)e 2x = C 1 e x + C 2

f(x) = C 1 e 3x + C 2 e 2x

と計算しても良いでしょう。

1.3.2

平方完成法 〜自明な2階方程式に帰着する〜

同様の記法を使えば同じ微分方程式はあたかも左辺を平方完成するかの様に

0 = f 00 (x) 5f 0 (x) + 6f (x) =

µ d dx 5

2

∂ 2

f (x) 25

4 f (x) + 6f (x) µ d

dx 5 2

∂ 2

f (x) = 1 4 f (x)

と変形する事が出来ます。ここで

f (x) = g(x)e

52

x

と置けば、最初に注意した事実から

µ d

dx 5 2

∂ 2

f (x) = g 00 (x)e

52

x

となっているので、先の方程式は

g 00 (x)e

52

x = 1

4 g(x)e

52

x

すなわち、

g 00 (x) = 1 4 g(x)

と変形する事が出来ます。これは簡単に一般解が求められて

g(x) = C 1 e

12

x + C 2 e

12

x

でした(

C 1 , C 2

は任意定数)から、結局もとの

f (x)

に戻してやれば

g(x)e

52

x = C 1 e

12

x e

52

x + C 2 e

12

x e

52

x

すなわち、

f (x) = C 1 e 3x + C 2 e 2x

という具合に求める一般解が求まりました。

1.4

実重解の場合の具体例

f 00 (x) 6f 0 (x) + 9f (x) = 0

この微分方程式も同様の記法に従えば

µ d dx 3

∂ 2

f (x) = 0

と書く事が出来ますが、ここで

f(x) = g(x)e 3x

と置けば

0 =

µ d dx 3

∂ 2

f (x) = g 00 (x)e 3x ,

従って

g 00 (x) = 0

となってこれは簡単に解が求まって

g(x) = C 1 x + C 2

C 1 , C 2

は任意の定数)

です。最後にこれを

f (x)

に戻してやれば一般解は次の通りです:

f (x) = C 1 xe 3x + C 2 e 3x .

1.5

共役な2つの複素数解の場合の具体例

1.5.1

平方完成法

f 00 (x) 2f 0 (x) + 5f (x) = 0

も、同様の記法によってあたかも左辺を平方完成するかの様に

µ d dx 1

∂ 2

f (x) f (x) + 5f (x) = 0 µ d

dx 1

∂ 2

f (x) = 4f (x)

と変形する事が出来ます。そこで同じ様に

f (x) = g(x)e x

と置けば、

µ d dx 1

∂ 2

(g(x)e x ) = 4g(x)e x g 00 (x)e x = 4g(x)e x

g 00 (x) = 4g(x)

(3)

Revised at 18:33, January 31, 2017

数学特論B 第13回

http://my.reset.jp/˜gok/math/advanced/ 3

となりますが、これは簡単に一般解が分かります(

C 1 , C 2

は任意の定数)

g(x) = C 1 cos 2x + C 2 sin 2x

元に戻してやれば解を得る事が出来ます:

f(x)e x = C 1 cos 2x + C 2 sin 2x f (x) = C 1 e x cos 2x + C 2 e x sin 2x.

2

非同次方程式

2.1

同次の場合と同様な因数分解法で頑張るやり方

対応した同次方程式の特性方程式が実数解(重解でも

OK

)を持つ場合にはこの方法 で解を求める事が出来ます(虚数解の場合も全く同様に可能ではありますが、虚数を含 んだ計算をしなければならない)。

2階の定数係数・非同次微分方程式:

f 00 (x) 5f 0 (x) + 6f (x) = (12x 7)e x

は左辺を因数分解するかの様に

µ d dx 3

∂ Ωµ d dx 2

f (x)

æ

= (12x 7)e x

と変形する事が出来、ここで

° d

dx 2 ¢

f (x) = h(x)

と置けば、方程式は

µ d

dx 3

h(x) = (12x 7)e x

と書く事が出来ます。両辺に

e 3x

を掛けてやる事によって

e 3x h 0 (x) 3e 3x h(x) = (12x 7)e 4x

© e 3x h(x)0

= (12x 7)e 4x e 3x h(x) =

Z

(12x 7)e 4x dx

= (12x 7) 1

4 e 4x Z

12 1

4 e 4x dx

= 12x 7

4 e 4x 3

4 e 4x + C 1

h(x) = ( 3x + 1)e x + C 1 e 3x

と云う具合に一般解

h(x) = ( 3x + 1)e x + C 1 e 3x

が求まります(

C 1

は任意の定数)。

そこでこの結果を元の

f (x)

に戻してやれば

µ d

dx 2

f (x) = ( 3x + 1)e x + C 1 e 3x

f 0 (x) 2f (x) = ( 3x + 1)e x + C 1 e 3x f 0 (x)e 2x 2f (x)e 2x = ( 3x + 1)e 3x + C 1 e x

° f (x)e 2x ¢ 0

= ( 3x + 1)e 3x + C 1 e x f (x)e 2x =

Z

( 3x + 1)e 3x dx + C 1 e x + C 2

= ( 3x + 1) 1

3 e 3x Z

( 3) 1

3 e 3x dx + C 1 e x + C 2

= 3x 1

3 e 3x + 1

3 e 3x + C 1 e x + C 2

f (x) = xe x + C 1 e 3x + C 2 e 2x

となって一般解が求まります(

C 1 , C 2

は任意の定数)。

2.2

解の構造に注目するやり方

一般に非同次方程式の解は

(非同次方程式の一般解)

=

(非同次方程式の1つの解)

+

(同次方程式の一般解)

と云う構造をしていますので、同次方程式の一般解の方は今見た様に求めて、後は1つ の解を何らかの方法で見つけてくれば良いわけです。

未知関数を含まない部分がこの様に『(1次式)×(指数関数)』であるならば、広く 知られているテクニックとして、

f (x) = (ax + b)e x

と仮定して、これがこの非同次方 程式を満たす様に係数

a, b

を決めてしまうと云う方法があります。

実際、

f (x) = (ax + b)e x

を方程式の左辺に代入してみると、

f 00 (x) 5f 0 (x) + 6f (x) = (ax 2a + b)e x 5( ax + a b)e x + 6(ax + b)e x

= (12ax 7a + 12b)e x

ですからこれが左辺

(12x 7)e x

と一致するためには

a = 1, b = 0

です。

従ってこの非同次方程式の1つの解

xe x

が求まりました。これと同次方程式の一般 解(これはさっき求めましたね)を使えば結局非同次方程式の一般解は

f (x) = xe x + C 1 e 3x + C 2 e 2x

C 1 , C 2

は任意の定数)

となる事が分かります(当たり前ですがさっきの解と同じです)。

(4)

Revised at 18:33, January 31, 2017

数学特論B 第13回

http://my.reset.jp/˜gok/math/advanced/ 4

Exercise

基本演習

1 (

阪大基礎工

9)

微分方程式(*)

y 00 2y 0 + 5y = e x + xe x

について、

以下の問いに答えて下さい。

(1)対応した同次方程式の一般解を求めて下さい。

(2)

y 00 2y 0 + 5y = e x

の解を1つ求めて下さい。

(3)

y 00 2y 0 + 5y = xe x

の解を1つ求めて下さい。

(4)(*)の一般解を求めて下さい。

(5)

y(0) = a, y 0 (0) = b

のときの(*)の解を求めて下さい。

基本演習

2 (

阪大基礎工

13) y 00 + ay 0 + by = 0

の解

y(x)

を、

(1)

b = 2a 2

(2)

b = a 4

2 (3)

b = 2a 2

のときにそれぞれ求めて下さい(

a

は実定数とします)。

基本演習

3 (

東大工

11) y 00 2y 0 + 5y = e x , y(0) = p, y 0 (0) = q

の解を求む。

基本演習

4 (

京大工

13) y 00 4y 0 + 3y = 0

の一般解を求めて下さい。

基本演習

5 (

神戸大理数

17) y 00 (a + b)y 0 + aby = 0

の一般解を求めて下さい。

基本演習

6 (

名工大

H9)

(1)

y 00 + 2y 0 3y = xe x

を解くために、

y(x) = z(x)e x

とおくと、微分方程式

z 00 +az 0 + bz = x

が導かれる。定数

a, b

の値を求めて下さい。

(2)上で導かれた微分方程式

z 00 + az 0 + bz = x

の一般解を求めて下さい。

基本演習

7 (

名工大

14)

定数係数2階線形微分方程式

y 00 y 0 2y = 20 cos 2x

初期条件

y(0) = 2, y 0 (0) = 4

を満たす解

y(x)

を求めて下さい。

基本演習

8 (

東工大

15)

微分方程式

y 00 + 2y 0 + y = e x

の解で、

y(0) = 1, y 0 (0) = 4 e

を満たすものを求めて下さい。

発展演習

9 (

京大工

14) y 00 + 9y = 7 sin 3x

において

z(x) = ° d

dx + 3i ¢

y(x)

として 以下の問いに答えて下さい。

(1)

° d

dx 3i ¢

z(x) = 7 sin 3x

z(x)

について解いて下さい。

(2)今求めた

z(x)

に対し

° d

dx + 3i ¢

y(x) = z(x)

y(x)

について解いて下さい。

基本演習

10 (

農工大

26)

微分方程式

y 00 + 4y 0 + 4y = 1+x e

2x2 の解

y = y(x)

のうち で条件

y(0) = 0, y 0 (0) = 1 2

を満たすものを求めて下さい。

基本演習

11 (

農工大

27)

微分方程式

y 00 + 3y = cos

3x

の解

y = y(x)

が、

y(0) = 1, y 0 (0) = 1

を満たすとき

y

を求めてください。

基本演習

12 (

東工大

26)

(1)

y 00 + y = 0

の一般解を求めて下さい。

(2)

y 00 + y = e x

の一般解を求めて下さい。

基本演習

13 (

神戸大

27)

以下の各微分方程式の一般解を求めて下さい。

(1)

2y 00 − 5y 0 + 2y = 0

(2)

2y 00 − 5y 0 + 2y = e x .

基本演習

14 (

三重大

26)

微分方程式

y 00 4y 0 + 4y = 0

の解を求めて下さい。

基本演習

15 (

北大

26)

微分方程式

y 00 + 3y 0 + 2y = 0

の一般解を求めて下さい。

基本演習

16 (

神戸大

26) y 00 2y 0 = xe 2x

について以下の問いに答えて下さい。

(1)この方程式は

y = (Ax 2 + Bx)e 2x (A, B

は定数

)

の形の特殊解を持つこと を示し、

A, B

を決めて下さい。

(2)この方程式の一般解を求めて下さい。

基本演習

17 (

滋賀県立大

26)

(1)

y 00 8y 0 + 16y = 0

の一般解を求めて下さい。

(2)2階定数係数非同次線形常微分方程式

y 00 8y 0 + 16y = 2 cos x

の特殊解 を求めて下さい

(

特殊解を

y(x) = A sin x + B cos x

と仮定してよい)。

(3)上記(2)の非同次線形常微分方程式の一般解を書き下して下さい。

発展演習

18 (

佐賀大

26) f (t) = Ae ∏t sin(ωt + θ)

が解となるような、

t

を独立変 数とする

f

の2階微分方程式を一つ書いて下さい(

A, ∏, ω, θ

は定数)。

発展演習

19 (

千葉大

27)

次の微分方程式の一般解

y(x)

を求め、与えられた初期 条件を満たす解曲線の概形を図示してください。

y 00 + 2y 0 + 17y = 0, x 0

初期条件:

y(0) = 1, y 0 (0) = 1

発展演習

20 (

奈良女

26)

(1)

y 0 = x y+1

2

+2x

の一般解を求めて下さい。

(2)解の形として

x = ae iωt

を仮定し、以下に示す手順で微分方程式

mx 00 (t) = kx ∏x 0 (t)

を解く事にします。ここで

a

は正の実定数で

ω

は複素定数とすします。

(a)

x 0 (t)

を計算し、それを

x

を用いて表して下さい。

(b)

x 00 (t)

を計算し、それを

x

を用いて表して下さい。

(c)

ω

が満たすべき方程式を導いて下さい。

(d)上で求めた方程式を解くことによって

ω

を求めて下さい。

参照

関連したドキュメント

Yamamoto: “Numerical verification of solutions for nonlinear elliptic problems using L^{\infty} residual method Journal of Mathematical Analysis and Applications, vol.

[r]

[r]

この節では mKdV 方程式を興味の中心に据えて,mKdV 方程式によって統制されるような平面曲線の連 続朗変形,半離散 mKdV

Existence of weak solution for volume preserving mean curvature flow via phase field method. 13:55〜14:40 Norbert

ある周波数帯域を時間軸方向で複数に分割し,各時分割された周波数帯域をタイムスロット

In this paper, we consider the discrete deformation of the discrete space curves with constant torsion described by the discrete mKdV or the discrete sine‐Gordon equations, and

(火力発電のCO 2 排出係数) - 調整後CO 2 排出係数 0.573 全電源のCO 2 排出係数