常微分方程式を解く方法
指数根と を持つ正規特異点を持つ微分方程式のいくつかの解 。
r
=
1
2
{\displaystyle r={\frac {1}{2}}}
−
1
{\displaystyle -1}
数学 において 、 フェルディナント・ゲオルク・フロベニウス にちなんで名付けられた フロベニウス法は 、
および の 形式の
線形2階常 微分方程式の 無限級数 解を求める方法です 。
ず
2
あなた
″
+
p
(
ず
)
ず
あなた
′
+
q
(
ず
)
あなた
=
0
{\displaystyle z^{2}u''+p(z)zu'+q(z)u=0}
あなた
′
≡
d
あなた
d
ず
{\textstyle u'\equiv {\frac {du}{dz}}}
あなた
″
≡
d
2
あなた
d
ず
2
{\textstyle u''\equiv {\frac {d^{2}u}{dz^{2}}}}
正規特異点 の近傍において 。
ず
=
0
{\displaystyle z=0}
で割ると、 形式の微分方程式が得られますが
、 p ( z )/ z または q ( z )/ z 2 のいずれかが z = 0 で解析的でない場合は、通常のべき級数法では解くことができません。フロベニウス法では、 p ( z ) と q ( z ) 自体 が 0 で 解析 的 で
ある か 、 または 他の場所では解析的であり、0 での極限が両方とも存在する(そして有限である)ことを条件に、このような微分方程式の べき 級数 解 を 作成 でき ます 。
ず
2
{\displaystyle z^{2}}
あなた
″
+
p
(
ず
)
ず
あなた
′
+
q
(
ず
)
ず
2
あなた
=
0
{\displaystyle u''+{\frac {p(z)}{z}}u'+{\frac {q(z)}{z^{2}}}u=0}
歴史: フロベニウスの実際の貢献
フロベニウスの貢献 [1]は、関係する級数解の可能なすべての 形式 (下記参照)という点においてはそれほど大きくない。これらの形式はすべて、以前に [2]フックス[ 3] によって確立されていた 。 [4] 指標 多項式 (下記参照)とその役割もフックスによって確立されていた。 [2]
フロベニウスによる理論への最初の貢献は、解析的冪級数に独立変数の任意の冪 rを乗じた形をとる最初の線形独立解に関して、一般化冪級数の係数は 再帰関係 に従い、常に簡単に計算できるということを示したことである (下記参照) 。
フロベニウスの2番目の貢献は、指数方程式の根が整数だけ異なる場合、 2番目の線形独立解の一般 形 (下記参照)は、上記のパラメータ rに関する微分 [5] に基づく手順によって得られることを示したことである。
フロベニウスの1873年の出版物[1] の大部分は、 解に含まれるすべての級数の収束の証明と、これらの級数の収束半径の確立に費やされました。
フロベニウス法の説明:まず、線形独立解
フロベニウスの方法は、次の形式のべき級数解を求めるものである。
あなた
(
ず
)
=
ず
r
∑
け
=
0
∞
あ
け
ず
け
、
(
あ
0
≠
0
)
{\displaystyle u(z)=z^{r}\sum _{k=0}^{\infty }A_{k}z^{k},\qquad (A_{0}\neq 0)}
差別化:
あなた
′
(
ず
)
=
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
1
{\displaystyle u'(z)=\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-1}}
あなた
″
(
ず
)
=
∑
け
=
0
∞
(
け
+
r
−
1
)
(
け
+
r
)
あ
け
ず
け
+
r
−
2
{\displaystyle u''(z)=\sum _{k=0}^{\infty }(k+r-1)(k+r)A_{k}z^{k+r-2}}
上記の微分を元の常微分方程式に代入すると次のようになります。
ず
2
∑
け
=
0
∞
(
け
+
r
−
1
)
(
け
+
r
)
あ
け
ず
け
+
r
−
2
+
ず
p
(
ず
)
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
1
+
q
(
ず
)
∑
け
=
0
∞
あ
け
ず
け
+
r
=
∑
け
=
0
∞
(
け
+
r
−
1
)
(
け
+
r
)
あ
け
ず
け
+
r
+
p
(
ず
)
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
+
q
(
ず
)
∑
け
=
0
∞
あ
け
ず
け
+
r
=
∑
け
=
0
∞
[
(
け
+
r
−
1
)
(
け
+
r
)
あ
け
ず
け
+
r
+
p
(
ず
)
(
け
+
r
)
あ
け
ず
け
+
r
+
q
(
ず
)
あ
け
ず
け
+
r
]
=
∑
け
=
0
∞
[
(
け
+
r
−
1
)
(
け
+
r
)
+
p
(
ず
)
(
け
+
r
)
+
q
(
ず
)
]
あ
け
ず
け
+
r
=
[
r
(
r
−
1
)
+
p
(
ず
)
r
+
q
(
ず
)
]
あ
0
ず
r
+
∑
け
=
1
∞
[
(
け
+
r
−
1
)
(
け
+
r
)
+
p
(
ず
)
(
け
+
r
)
+
q
(
ず
)
]
あ
け
ず
け
+
r
=
0
{\displaystyle {\begin{aligned}&z^{2}\sum _{k=0}^{\infty }(k+r-1)(k+r)A_{k}z^{k+r- 2}+zp(z)\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-1}+q(z)\sum _{k=0}^{\infty }A_{k}z^{k+r}\\={}&\sum _{k=0}^{\infty }(k+r-1)(k +r)A_{k}z^{k+r}+p(z)\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r}+q(z)\sum _{k=0}^{\infty }A_{k}z^{k+r}\\={}&\sum _{k=0}^{\infty }[(k+r-1)(k+r)A_{k}z^ {k+r}+p(z)(k+r)A_{k}z^{k+r}+q(z)A_{k}z^{k+r}]\\={}&\合計 _{k=0}^{\infty }\left[(k+r-1)(k+r)+p(z)(k+r)+q(z)\right]A_{k}z ^{k+r}\\={}&\left[r(r-1)+p(z)r+q(z)\right]A_{0}z^{r}+\sum _{k =1}^{\infty }\left[(k+r-1)(k+r)+p(z)(k+r)+q(z)\right]A_{k}z^{k+r}=0\end{整列した
この式は
、 r の 2 次式である 指数多項式
として知られています。 指数多項式 の一般的な定義は、無限級数における z の最低のべき乗の係数です。この場合、これは r 番目の係数になります が、最小の指数は、与えられた微分方程式に応じて r − 2、 r − 1 などになる可能性があります。この詳細は覚えておくことが重要です。微分方程式のすべての級数を同じインデックス値 (上記の式では k = 1) から開始するように同期させるプロセスでは、複雑な式になる可能性があります。ただし、指数根を解く際には、 z の最低のべき乗の係数のみに注意が向けられます 。
r
(
r
−
1
)
+
p
(
0
)
r
+
q
(
0
)
=
私
(
r
)
{\displaystyle r\left(r-1\right)+p\left(0\right)r+q\left(0\right)=I(r)}
これを用いると、 z k + r の係数の一般的な表現は 次のようになる。
私
(
け
+
r
)
あ
け
+
∑
じゅう
=
0
け
−
1
(
じゅう
+
r
)
p
(
け
−
じゅう
)
(
0
)
+
q
(
け
−
じゅう
)
(
0
)
(
け
−
じゅう
)
!
あ
じゅう
、
{\displaystyle I(k+r)A_{k}+\sum _{j=0}^{k-1}{(j+r)p^{(kj)}(0)+q^{(kj)}(0) \over (kj)!}A_{j},}
これらの係数は微分方程式の解であるためゼロでなければならないので、
私
(
け
+
r
)
あ
け
+
∑
じゅう
=
0
け
−
1
(
じゅう
+
r
)
p
(
け
−
じゅう
)
(
0
)
+
q
(
け
−
じゅう
)
(
0
)
(
け
−
じゅう
)
!
あ
じゅう
=
0
∑
じゅう
=
0
け
−
1
(
じゅう
+
r
)
p
(
け
−
じゅう
)
(
0
)
+
q
(
け
−
じゅう
)
(
0
)
(
け
−
じゅう
)
!
あ
じゅう
=
−
私
(
け
+
r
)
あ
け
1
−
私
(
け
+
r
)
∑
じゅう
=
0
け
−
1
(
じゅう
+
r
)
p
(
け
−
じゅう
)
(
0
)
+
q
(
け
−
じゅう
)
(
0
)
(
け
−
じゅう
)
!
あ
じゅう
=
あ
け
{\displaystyle {\begin{aligned}I(k+r)A_{k}+\sum _{j=0}^{k-1}{(j+r)p^{(kj)}(0)+q^{(kj)}(0) \over (kj)!}A_{j}&=0\\[4pt]\sum _{j=0}^{k-1}{(j+r)p^{(kj)}(0)+q^{(kj)}(0) \over (kj)!}A_{j}&=-I(k+r)A_{k}\\[4pt]{1 \over -I(k+r)}\sum _{j=0}^{k-1}{(j+r)p^{(kj)}(0)+q^{(kj)}(0) \over (kj)!}A_{j}&=A_{k}\end{aligned}}}
上記のA k の級数解
は
、
あなた
r
(
ず
)
=
∑
け
=
0
∞
あ
け
ず
け
+
r
{\displaystyle U_{r}(z)=\sum _{k=0}^{\infty }A_{k}z^{k+r}}
ず
2
あなた
r
(
ず
)
″
+
p
(
ず
)
ず
あなた
r
(
ず
)
′
+
q
(
ず
)
あなた
r
(
ず
)
=
私
(
r
)
ず
r
{\displaystyle z^{2}U_{r}(z)''+p(z)zU_{r}(z)'+q(z)U_{r}(z)=I(r)z^{ r}}
U r ( z ) における r の指示多項式の根の 1 つを選択すると 、微分方程式の解が得られます。根の差が整数でない場合は、もう一方の根で線形独立な別の解が得られます。
例
解決しましょう
ず
2
ふ
″
−
ず
ふ
′
+
(
1
−
ず
)
ふ
=
0
{\displaystyle z^{2}f''-zf'+(1-z)f=0}
全体をz 2 で割ると、 z = 0
で必要な特異点を持つ ようになります
。
ふ
″
−
1
ず
ふ
′
+
1
−
ず
ず
2
ふ
=
ふ
″
−
1
ず
ふ
′
+
(
1
ず
2
−
1
ず
)
ふ
=
0
{\displaystyle f''-{1 \over z}f'+{1-z \over z^{2}}f=f''-{1 \over z}f'+\left({1 \over z^{2}}-{1 \over z}\right)f=0}
シリーズソリューションを使用する
ふ
=
∑
け
=
0
∞
あ
け
ず
け
+
r
ふ
′
=
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
1
ふ
″
=
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
{\displaystyle {\begin{aligned}f&=\sum _{k=0}^{\infty }A_{k}z^{k+r}\\f'&=\sum _{k=0}^ {\infty }(k+r)A_{k}z^{k+r-1}\\f''&=\sum _{k=0}^{\infty }(k+r)(k+r-1)A_{k}z^{k+r-2}\end{整列}}}
さて、代入すると
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
−
1
ず
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
1
+
(
1
ず
2
−
1
ず
)
∑
け
=
0
∞
あ
け
ず
け
+
r
=
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
−
1
ず
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
1
+
1
ず
2
∑
け
=
0
∞
あ
け
ず
け
+
r
−
1
ず
∑
け
=
0
∞
あ
け
ず
け
+
r
=
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
−
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
2
+
∑
け
=
0
∞
あ
け
ず
け
+
r
−
2
−
∑
け
=
0
∞
あ
け
ず
け
+
r
−
1
=
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
−
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
2
+
∑
け
=
0
∞
あ
け
ず
け
+
r
−
2
−
∑
け
−
1
=
0
∞
あ
け
−
1
ず
け
−
1
+
r
−
1
=
∑
け
=
0
∞
(
け
+
r
)
(
け
+
r
−
1
)
あ
け
ず
け
+
r
−
2
−
∑
け
=
0
∞
(
け
+
r
)
あ
け
ず
け
+
r
−
2
+
∑
け
=
0
∞
あ
け
ず
け
+
r
−
2
−
∑
け
=
1
∞
あ
け
−
1
ず
け
+
r
−
2
=
{
∑
け
=
0
∞
(
(
け
+
r
)
(
け
+
r
−
1
)
−
(
け
+
r
)
+
1
)
あ
け
ず
け
+
r
−
2
}
−
∑
け
=
1
∞
あ
け
−
1
ず
け
+
r
−
2
=
{
(
r
(
r
−
1
)
−
r
+
1
)
あ
0
ず
r
−
2
+
∑
け
=
1
∞
(
(
け
+
r
)
(
け
+
r
−
1
)
−
(
け
+
r
)
+
1
)
あ
け
ず
け
+
r
−
2
}
−
∑
け
=
1
∞
あ
け
−
1
ず
け
+
r
−
2
=
(
r
−
1
)
2
あ
0
ず
r
−
2
+
{
∑
け
=
1
∞
(
け
+
r
−
1
)
2
あ
け
ず
け
+
r
−
2
−
∑
け
=
1
∞
あ
け
−
1
ず
け
+
r
−
2
}
=
(
r
−
1
)
2
あ
0
ず
r
−
2
+
∑
け
=
1
∞
(
(
け
+
r
−
1
)
2
あ
け
−
あ
け
−
1
)
ず
け
+
r
−
2
{\displaystyle {\begin{aligned}\sum _{k=0}^{\infty }&(k+r)(k+r-1)A_{k}z^{k+r-2}-{\frac {1}{z}}\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-1}+\left({\frac {1}{z^{2}}}-{\frac {1}{z}}\right)\sum _{k=0}^{\infty }A_{k}z^{k+r}\\&=\sum _{k=0}^{\infty }(k+r)(k+r-1)A_{k}z^{k+r-2}-{\frac {1}{z}}\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-1}+{\frac {1}{z^{2}}}\sum _{k=0}^{\infty }A_{k}z^{k+r}-{\frac {1}{z}}\sum _{k=0}^{\infty }A_{k}z^{k+r}\\&=\sum _{k=0}^{\infty }(k+r)(k+r-1)A_{k}z^{k+r-2}-\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-2}+\sum _{k=0}^{\infty }A_{k}z^{k+r-2}-\sum _{k=0}^{\infty }A_{k}z^{k+r-1}\\&=\sum _{k=0}^{\infty }(k+r)(k+r-1)A_{k}z^{k+r-2}-\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-2}+\sum _{k=0}^{\infty }A_{k}z^{k+r-2}-\sum _{k-1=0}^{\infty }A_{k-1}z^{k-1+r-1}\\&=\sum _{k=0}^{\infty }(k+r)(k+r-1)A_{k}z^{k+r-2}-\sum _{k=0}^{\infty }(k+r)A_{k}z^{k+r-2}+\sum _{k=0}^{\infty }A_{k}z^{k+r-2}-\sum _{k=1}^{\infty }A_{k-1}z^{k+r-2}\\&=\left\{\sum _{k=0}^{\infty }\left((k+r)(k+r-1)-(k+r)+1\right)A_{k}z^{k+r-2}\right\}-\sum _{k=1}^{\infty }A_{k-1}z^{k+r-2}\\&=\left\{\left(r(r-1)-r+1\right)A_{0}z^{r-2}+\sum _{k=1}^{\infty }\left((k+r)(k+r-1)-(k+r)+1\right)A_{k}z^{k+r-2}\right\}-\sum _{k=1}^{\infty }A_{k-1}z^{k+r-2}\\&=(r-1)^{2}A_{0}z^{r-2}+\left\{\sum _{k=1}^{\infty }(k+r-1)^{2}A_{k}z^{k+r-2}-\sum _{k=1}^{\infty }A_{k-1}z^{k+r-2}\right\}\\&=(r-1)^{2}A_{0}z^{r-2}+\sum _{k=1}^{\infty }\left((k+r-1)^{2}A_{k}-A_{k-1}\right)z^{k+r-2}\end{aligned}}}
( r − 1) 2 = 0 から 1 の二重根が得られます。この根を使用して、 z k + r − 2 の係数を 0 に設定します (解となるため)。これにより、次の式が得られます。
したがって、再帰関係が得られます。
(
k
+
1
−
1
)
2
A
k
−
A
k
−
1
=
k
2
A
k
−
A
k
−
1
=
0
{\displaystyle (k+1-1)^{2}A_{k}-A_{k-1}=k^{2}A_{k}-A_{k-1}=0}
A
k
=
A
k
−
1
k
2
{\displaystyle A_{k}={\frac {A_{k-1}}{k^{2}}}}
いくつかの初期条件が与えられれば、再帰を完全に解くことも、べき級数形式で解を得ることもできます。
係数の比は 有理関数 なので 、べき級数は 一般化された超幾何級数 として表すことができます。
A
k
/
A
k
−
1
{\displaystyle A_{k}/A_{k-1}}
「例外的なケース」: 整数で区切られたルート
前の例では、重複根を持つ指数多項式が使用されており、与えられた微分方程式に対して 1 つの解しか得られません。一般に、フロベニウス法では、指数方程式の根が整数 (ゼロを含む) で区切られていない場合、2 つの独立した解が得られます。
根が重複しているか、根が整数だけ異なる場合、2 番目の解は次を使用して見つけることができます。
ここで 、 は最初の解 (不等根の場合は大きい方の根に基づく)、 は小さい方の根、定数 C と係数 を決定する必要があります。 が選択されると (たとえば 1 に設定)、までの C と が 決定されますが、は含まれません 。 は任意に設定できます。これにより、 の残りの部分が決定されます。 場合によっては、定数 C を 0 にする必要があります。
y
2
=
C
y
1
ln
x
+
∑
k
=
0
∞
B
k
x
k
+
r
2
{\displaystyle y_{2}=Cy_{1}\ln x+\sum _{k=0}^{\infty }B_{k}x^{k+r_{2}}}
y
1
(
x
)
{\displaystyle y_{1}(x)}
r
2
{\displaystyle r_{2}}
B
k
{\displaystyle B_{k}}
B
0
{\displaystyle B_{0}}
B
k
{\displaystyle B_{k}}
B
r
1
−
r
2
{\displaystyle B_{r_{1}-r_{2}}}
B
k
.
{\displaystyle B_{k}.}
例 : 次の微分方程式 ( a = 1 および b = 2 の Kummer 方程式 ) を考えます。
指数方程式の根は -1 と 0 です。 2 つの独立した解は であり 、 どの解にも対数は現れないことがわかります。 解は、 べき乗ゼロから始まるべき級数を持ちます。 べき級数で始まる 再帰関係では、項の係数に制限はなく 、任意に設定できます。 係数をゼロに設定すると、この微分方程式では他のすべての係数がゼロになり、解 1/ z が得られます。
z
u
″
+
(
2
−
z
)
u
′
−
u
=
0
{\displaystyle zu''+(2-z)u'-u=0}
1
/
z
{\displaystyle 1/z}
e
z
/
z
,
{\displaystyle e^{z}/z,}
(
e
z
−
1
)
/
z
{\displaystyle (e^{z}-1)/z}
z
−
1
{\displaystyle z^{-1}}
z
0
,
{\displaystyle z^{0},}
例外的な場合の級数係数の直列再帰関係
指数多項式の根が整数(ゼロを含む)だけ異なる場合、2番目の線形独立解に含まれるすべての級数の係数は、 直列再帰関係 から直接計算できます。 [5]これらの直列関係は、パラメータ r について微分するというフロベニウスの独自の発明をさらに発展させ 、このアプローチを使用してすべてのケースで級数の係数を実際に計算することによって構築できます。 [5]
参照
外部リンク
参考文献
^ ab Frobenius、Ferdinand Georg (1968) [原文は Journal für die reine und angewandte Mathematik 76、214-235 (1873)]。 「Uber die Integration der lineen Differentialgleichungen durch Reihen」。 ゲザメルテ・アブハンドルンゲン (ドイツ語)。ベルリン: Springer-Verlag。 84–105ページ。
^ ab グレイ、ジェレミー (1986)。 リーマンからポアンカレまでの線型微分方程式と群論 。ボストン:バークハウザー 。ISBN 0-8176-3318-9 。
^ フックス、ラザロ・インマヌエル (1865)。 "Zur Theorie der lineen Differentialgleichungen mit veranderlichen Coefficienten"。 Gesammelte Mathematische Werke von L. Fuchs (ドイツ語)。ミシガン大学の図書館。
^ フックス、ラザロ・インマヌエル (1866)。 "Zur Theorie der lineen Differentialgleichungen mit veranderlichen Coefficienten"。 数学に関するジャーナル 。 66 : 159–204。
^ abc van der Toorn, Ramses (2022年12月27日). 「正規特異点に関する対数フロベニウス級数解の係数のタンデム再帰関係」. 公理 . 12 (1): 32. doi : 10.3390/axioms12010032 . ISSN 2075-1680.