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

微分方程式

N/A
N/A
Protected

Academic year: 2021

シェア "微分方程式"

Copied!
18
0
0

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

全文

(1)

微分方程式

2. 同次型

(2)

1 復習・変数分離形

–1–

正規形(

y = · · ·

の形)の

1

階常微分方程式

y = F(x, y)

を考えている

変数分離形の常微分方程式とは

dy

dx = f(x)g(y)

の形のものをいう。このとき次のように変数を分離して

dy

g(y) = f(x)dx

一般解は

dy g(y) =

f(x)dx +C, (C :

定数

)

で与えられる

.

(3)

2 同次形

–2–

正規形の

1

階常微分方程式

y = F(x, y)

で、右辺

F(x, y)

y/x

の関数になっているとき、すなわち

y = f

(y x

)

の形の微分方程式を同次形という

.

同次形の場合には原則として変数変換

u = y x

を行うことにより変数分離形に帰着させることができる

.

(4)

3 同次形の解き方

–3–

y = f (y

x )

において

u = y x

とおくと

y = xu

なので

(

左辺

) : y = (xu) = u + xu (

右辺

) : f

(y x

)

= f(u)

したがって

u +xu = f(u)

よって

xu = f(u) u = u = f(u) u x

という変数分離形の常微分方程式に変形される

.

(5)

3 同次形の解き方・続き

–4–

変数分離形の常微分方程式

du

dx = f(u) u x

は変数分離形の解き方に従うと

du

f(u) u =

dx

x =

du

f(u) u = log|x| +C

左辺を積分したものを

g(u) =

du f(u) u

と表すと一般解は次で与えられる

g (y

x )

= log|x| +C

(6)

4 同次形の解き方・注意

–5–

同次形の方程式

y = f (y

x )

変数分離形の解き方と同様に

f(u) u = 0

を満たす 値

u = m

が存在するとき問題が起こる。

u = y/x

なので

y = mx

(直線)となる。実際にこのとき

y = m, f

(y x

)

= f(m) = m

となって

y = mx

は同次形の方程式の解になっている。

多くの場合は

y = mx

は一般解の中に含まれる

(7)

5 同次形・例

–6–

例1

. y = x2 +y2 2xy

右辺は分子分母を

x2

で割って

x2 + y2

2xy = 1 + (y/x)2

2y/x

となるので 同次形である。

u = y

x

とおくと

y = xu

より

y = xu +u,

右辺は

1 + u2

2u

より

xu +u = 1 + u2

2u = u = 1 −u2

2xu = u2 1 2u ·

(

1 x

)

変数分離形の解き方で

∫ 2u du

u2 1 =

dx x +C

したがって

log|u21| = log|x|+C log|x(u21)| = C x(u21) = ±eC

(8)

5 同次形・例(つづき)

–7–

x(u2 1) = ±eC

において、

D = ±eC

とおくと

D ̸= 0

で 

x(u2 1) = D.

u = y/x

を代入すれば

y2 x2 = Dx.

変数分離形の解き方で分母になる

u2 1 = 0

のとき、すなわち

u = ±1

のときが問題になる

.

u = ±1

のときは

y = ±x

となり

,

もとの微分方程式を満足してお り

,

さきほどの一般解で抜けていた

D = 0

に対応している

.

以上により、求める一般解は

y2 x2 = Dx.

(9)

5 同次形に帰着できる方程式

–8–

同次形の方程式

y = f (y

x )

よくでてくる同次形として

(A, B, a, b

は定数)

y = f

(Ax +By ax+by

)

= f

(A +By/x a +by/x

)

同次形でない方程式

(A, B, C, a, b, c

は定数)

y = f

(Ax +By +C ax+ by + c

)

この方程式は同次形に変換できる

(10)

5 同次形に帰着できる方程式 2

–9–

同次形でない方程式

(A, B, C, a, b, c

は定数)

y = f

(Ax +By +C ax+ by + c

)

もし

C = 0, c = 0

であれば同次形になる

: y = f

(Ax + By ax+ by

)

= f

(A +B y/x a +b y/x

)

したがって元の方程式を変形して同次形になるようにする

(11)

ここで変数変換

x = X +x0, y = Y +y0

を行うと

–10–

dy

dx = dY dX

(定数だけずらしても微分は変わらない)

dy dx = f

(Ax +By +C ax+ by +c

)

の文字を変えると

dY

dX = f

(A(X + x0) + B(Y + y0) + C a(X +x0) + b(Y + y0) + c

)

= f

(AX + BY + (Ax0 +By0 +C) aX + bY + (ax0 +by0 +c)

)

ここで

,

定数項が消えるようにうまく

x0, y0

を選ぶ

(12)

5 同次形に帰着できる方程式 2 ・続き

–11–

0) x0, y0

に関する連立

1

次方程式

Ax0 + By0 +C = 0 ax0 +by0 +c = 0

を満たす

x0, y0

が存在すれば

,

次の同次形を得る

: dY

dX = f

(AX +BY aX +bY

)

連立方程式は解が一意に定まるとき・不能・不定の3通りある

1)

一意

A/a ̸= B/b

のときは一意に解を持つ

(13)

5 同次形に帰着できる方程式 2 ・続き

–12–

2)

不能

A/a = B/b ̸= C/c

のときはは解がない 元の方程式に立ち返って

A/a = B/b = k

とおくと

y = f

(Ax +By +C ax+ by + c

)

= f

(k(ax+ by) + C ax+by +c

)

u = ax+by

と置くと

u = a+by

より

y = (u−a)/b.

したがって

1

b (u −a) = f

(ku + C u +c

)

x

の入らない方程式

u = bf

(ku +C u +c

) + a

になるので変数分離形になっている

(14)

5 同次形に帰着できる方程式 2 ・続き

–13–

3)

不定

A/a = B/b = C/c

のときは解を無数に持つ

A = ka, B = kb, C = kc

を元の方程式に代入すると

y = f

(Ax + By +C ax +by +c

)

= f

(kax +kby +kc ax +by +c

)

= f(k)

と実は右辺は定数になっているので

y = f(k)x +D (D

は定数)が解

(15)

6 まとめ

–14–

(I)

同次形の微分方程式

y = f (y

x )

は、

u = y/x

と置いて変数分離形にできる

(II)

同次形ではないが

y = f

(Ax +By +C ax+ by + c

)

は連立方程式

Ax0 + By0 +C = 0 ax0 +by0 +c = 0

を解いて同次形にできる

(16)

7 同次形・演習問題

–15–

(1) (3x2 +y2)y 2xy = 0. (2) y = 2x +y 3 x + 2y 5. [

追加問題

3.1]

(1) x + yy = 2y (2) (x +y) + (x y)y = 0 (3) y2 +x2y = xyy (4) xy2y = x3 +y3 (5) 2xy (3x2 +y2)y = 0

(6) (

xcos y

x + ysin y x

) y =

(

ysin y

x xcos y x

) xy [

追加問題

3.2]

(1) (5x7y)(x3y+4)y = 0 (2) (3x+y5)(x3y5)y = 0 (3) (4x 6y 1) (2x 3y + 2)y = 0

(4) (6x 2y 7) (3x −y + 4)y = 0 (5) y = 2

( y + 2 x + y + 1

)2

(6) y =

(x y + 3 x y + 1

)2

(17)

演習問題解説 (1) y = 2xy

3x2+y2 より同次形なので y =xu とおくとy =xu+u. よってxu+u= 2u 3 +u2. u について解くとu =u(u2+ 1)

x(u2+ 3). 変数分離形なので

u2+ 3

u(u2+ 1)du=

dx x . よって log

u3 (u2+ 1)

=log|x|+C. となるのでlog xu3

(u2+ 1)

=C より u3

(u2+ 1) = D x (た だしD=±eC とおいた). u=y/xだったので y3

(xy2+x3) = D x

(2) y = 2x+y3

x+ 2y5 において x=X+a, y=Y +b とおくとdY

dX = 2X+Y + 2a+b3 X+ 2Y +a+ 2b5. ここで連立一次方程式 2a+b3 = 0, a+ 2b5 = 0 をといて a= 1/3, b= 7/3.

以下ではdY

dX = 2X+Y

X+ 2Y を考える. Y =XuとおいてXu+u= 2 +u

1 + 2u よりu = 1

X·22u2 1 + 2u と変数分離形になる.

1 + 2u 22u2 du=

∫ (

3

4(u1) 1 4(u+ 1)

)

du=3

4log|u1| −1

4log|u+ 1| より3

4log|u1| − 1

4log|u+ 1|= log|x|+C となるのでx4(u1)3(u+ 1) =D となる。

ただし D=±e4C. u=Y /X = (y7/3)/(x1/3)を代入して y, x の陰関数 27x4(xy+ 2)3(83x3y)

(13x)4 =D を得る。

(18)

[問題] dy

dx = xy1 x+y+ 3

[解説]x=X+a, y =Y +b とおくと dY

dX = XY + (ab1) X+Y + (a+b+ 3).

連立方程式 ab1 = 0, a+b+ 3 = 0 をといてa=1, b=2. ここで u=Y /X とおくと u+Xu = 1u

1 +u. 正規系に直して

u = 1 X ·

(1u 1 +u u

)

= 1

X · 12uu2 1 +u これは変数分離形なので

1 +u

12uu2 du=

1 X dX よって

1

2log|u2+ 2u1|= log|X|+D となってlog|X2(u2+ 2u1)|=2D よりX2(u2+ 2u1) =±eD. C =±eD とおき、u=Y /XよりY2+ 2XY X2 =C.

ここでa=1, b =2なので

(y+ 2)2+ 2(x+ 1)(y+ 2)(x+ 1)2 =C.

参照

関連したドキュメント

[r]

しかし何かを不思議だと思うことは勉強をする最も良い動機だと思うので,興味を 持たれた方は以下の文献リストなどを参考に各自理解を深められたい.少しだけ案

[r]

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

参加方式 対面方式 オンライン方式 使用可能ツール zoom Microsoft Teams. 三重県 鈴鹿市平田中町1-1