Processing math: 100%
関連するページを見るにはこのグラフ図を利用してください.

階数低減法 (reduction of order)

n 階の微分演算子

L(n)1=ni=0ai(x)Di=an(x)Dn+an1(x)Dn1++a1(x)D+a0(x)     - - - (1)

を用いて表される n 階非同次線形微分方程式 L(n)1[y]=f(x)  を考える.ここで,D=d/dx である(微分演算子).同次方程式 L(n)1[y]=0  の解(独立な解は n 個)が一つでも分かっていれば,この非同次方程式の階数を一つ下げて, n1 階の線形微分方程式の形に直すことができる.その方法を 階数低減法(reduction of order) という.

2階線形微分方程式の場合の例

まず,同次方程式 L(n)1[y]=0  のゼロでない解 y1(x) が分かっているとする.つまり,

L(n)1[y1]=ni=0ai(x)y(i)1=0     - - - (2)

である.この解と,定数ではない x の未知関数 u(x) とを用いて,非同次方程式の一般解を

y(x)=u(x)y1(x)     - - - (3)

とおく.式 (3) に式 (1) を作用させて,

y(i)=diydxi =ij=0Ciju(j)y(ij)1

を用いると,

L(n)1[y]=ni=0ai(x)y(i)

=ni=0ai(x)ij=0iCju(j)y(ij)1

=ni=0ai(x){uy(i)1+ij=1iCju(j)y(ij)1}

=uni=0ai(x)y(i)1+ni=1ai(x)ij=1iCju(j)y(ij)1    - - - (4)

となる.式 (4) の第1項目は式 (2) より消え,第2項目の和の順序を入れ替えて j=0 からの和に変更すると,

L(n)1[y]=nj=1[ni=jai(x)iCjy(ij)1]u(j) =n1j=0[ni=j+1ai(x)iCj+1y(ij1)1]u(j+1)     - - - (5)

が得られる.ここで, w(x)=du/dx を導入すると, w(j)=u(j+1)  であり,式 (5) の [ ]の中を

bj=ni=j+1ai(x)iCj+1y(ij1)1

とおくと,

L(n)1[y]=n1j=0bj(x)w(j)=L(n1)2[w]   ,   L(n)2=ni=0bi(x)Di

と書ける.したがって, n 階の線形微分方程式 L(n)1[y]=f(x)  が n1 階の線形微分方程式 L(n1)2[w]=f(x) の形に直せるのである.この n1 階の線形微分方程式を(解くことができれば)解いて w を求め, u=wdx から u を求めてやれば,式 (3) より一般解 y が求まる.


◎ 2階線形微分方程式の例

階数低減法の例として最もよく知られているのは定数係数の2階同次線形微分方程式の特性方程式の解が実重根の場合(*)であるが,まずはより一般的な場合を考える.

2階非同次線形微分方程式

y+a1(x)y+a0(x)y=f(x)     - - - (6)

について,同次方程式

y+a1(x)y+a0(x)y=0     - - - (7)

の一つの解が y1 と分かっているとする. x の未知関数 u(x) を用いて, y=uy1 とおくと,

y=uy1+uy1   ,   y=uy1+2uy1+uy1

である.これらを式 (6) に代入すると,

uy1+2uy1+uy1+a1(uy1+uy1)+a0uy1=f
⇒   uy1+u(2y1+a1y1)+u(y1+a1y1+a0y1)=f

が得られる. y1 は式 (7) の解なので,上式左辺の第3項目の( )の中の式はゼロである.したがって, w=u  とおいて,両辺を y1 で割って整理すると,

y1w+(2y1+a1y1)w=f   ⇒   w+(2y1y1+a1)w=fy1     - - - (8)

となる.これが階数低減法によって得られた1階微分方程式である.この1階微分方程式の積分因子は

μ=exp[(2y1y1+a1)dx] =exp[2y1y1dx+a1dx] =exp[2log|y1|+a1dx] =y21exp[a1dx]

であり,この積分因子を用いて,一般解として

w=1μ[μfy1dx+C] =y21ea1dx[y1fea1dxdx+C]    ( C :任意定数)

が求まる.あとは, u=wdx から u を求めてやれば,解 y=uy1 が求まる.


◎ 定数係数の2階同次線形微分方程式の特性方程式の解が実重根の場合

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

y+ay+by=0    ( a,b :実定数)     - - - (9)

において, a24b=0  のときの解を求めることを考える.この場合,特性方程式 λ2+aλ+b=0 の解は

λ=a±a24b2=a2   (実重根)

であるので,式 (9) の一つの解は

y1=ea2x     - - - (10)

である. x の未知関数 u(x) を用いて, y=uy1 =uea2x とおいて,式 (9) に代入して整理すると,

uea2xauea2x+a24uea2x+a(uea2xa2uea2x)+buea2x=0

⇒   {u14(a24b)u}ea2x=0   ⇒   uea2x=0

となり, w=u  とおいて(今の場合,特におく必要はないが), ea2x0 であることを考えると,1階微分方程式 w=0 が得られる.この解は w=c1 であり,これより u=wdx=c1x+c2 が求まる( c1,c2 :任意定数).したがって,式 (9) の一般解は

y=uea2x =(c1x+c2)ea2x

となる.このことから,式 (9) のもう一つの独立な解は y2=xea2x であることが分かる.


ホーム>>カテゴリー分類>>微分>>微分方程式>>階数低減法

[ページトップ]

金沢工業大学

利用規約

google translate (English version)