|
|
1行目: |
1行目: |
− | {{Differential equations}} | + | {{テンプレート:20180815sk}} |
− | '''線型微分方程式'''<ref group="注">'''線形'''等の用字・表記の揺れについては[[線型性]]を参照。</ref>(せんけいびぶんほうていしき、{{lang-en-short|''linear differential equation''}})は、[[微分]]を用いた[[線型写像|線型作用素]](線型[[微分作用素]]){{mvar|L}} と[[未知関数]] {{mvar|y}} と既知関数 {{mvar|b}} を用いて
| |
− | : {{math|1=''Ly'' = ''b''}}
| |
− | の形に書かれる[[微分方程式]]のこと。
| |
− | | |
− | == 概要 ==
| |
− | 線型微分方程式
| |
− | :<math>Ly = b</math>
| |
− | は、{{math|''b'' ≠ 0}} の場合、2 つの解 {{math|''s''<sub>1</sub>, ''s''<sub>2</sub>}} を任意に取り、その差 {{math|''d'' {{=}} ''s''<sub>1</sub> − ''s''<sub>2</sub>}} を考えると、{{mvar|L}} が線型作用素であることから
| |
− | : <math>\begin{align}
| |
− | Ld &= L\left(s_1 - s_2\right)\\
| |
− | &=Ls_1 - Ls_2\\
| |
− | &=b - b\\
| |
− | &=0
| |
− | \end{align}</math>
| |
− | となり、{{math|''b'' {{=}} 0}} の場合に帰着する。この {{math|''b'' {{=}} 0}} の場合の線型微分方程式は(もとの方程式に属する)'''斉次'''あるいは'''同次'''な {{en|(homogeneous)}}{{efn2|name="homogeneous"|ここでいう {{en|homogeneous}} は[[斉次函数]]のような次数に関する語ではなく、解函数あるいは解空間のある種の「等質性」を表すために用いられており、むしろ[[等質空間]]などでの語法が近い。しかし、斉次(形、方程式)・同次(形、方程式)と訳すのが定訳であり、等質方程式や非等質形のように呼ぶことはないかあってもかなり稀。}}方程式と呼ばれる。{{math|''s''<sub>1</sub> {{=}} ''d'' + ''s''<sub>2</sub>}} であることを考えれば線型微分方程式 {{math|''Ly'' {{=}} ''b''}} のすべての解は {{math|''Ly'' {{=}} ''b''}} の[[特殊解]]と、元の方程式に対応する斉次方程式
| |
− | :<math>Ly = 0</math>
| |
− | の解の和となる。したがって、線型微分方程式を解くことは特殊解を見つける問題と、斉次方程式を解く問題に分けることができる。また、{{mvar|L}} が線型作用素であることから、斉次方程式の解は[[線型性]]を持ち、解同士の和や、解の定数倍も解になる。
| |
− | | |
− | 関数の代わりに[[数列]]を(同時に、[[微分]]の代わりに[[有限差分|差分]]を)考えると、類似の概念として[[漸化式]](差分方程式)を捉えることができる(離散化)。線型差分方程式と線型微分方程式の間で、特性方程式を用いる解法など、いくつかの手法を共通に用いることができる。
| |
− | | |
− | == 定義 ==
| |
− | === 高階単独型 ===
| |
− | {{mvar|x}} の関数 {{mvar|y}} の高階微分 {{math|{{sfrac|''d <sup>j</sup>y''|''dx <sup>j</sup>''}}}} および、可微分関数 {{math|''a<sub>j</sub>''(''x'') (1 ≤ ''j'' ≤ ''n''), ''b''(''x'')}} により
| |
− | :<math>\frac{d^n y}{dx^n} + a_{n-1}(x)\frac{d^{n-1}y}{dx^{n-1}} + \cdots + a_1(x)y = b(x)</math>
| |
− | で表される微分方程式を'''単独高階型の線型微分方程式'''という。{{math|''b'' {{=}} 0}} であるとき'''斉次'''{{efn2|name="homogeneous"}}あるといい、
| |
− | :<math>\frac{d^n y}{dx^n} + a_{n-1}(x)\frac{d^{n-1}y}{dx^{n-1}} + \cdots + a_1(x)y = 0</math>
| |
− | を元の方程式に'''属する斉次方程式'''という。
| |
− | | |
− | 微分作用素 {{math|''L''({{sfrac|''d''|''dx''}})}} を
| |
− | :<math>L\left(\frac{d}{dx}\right) = \frac{d^n}{dx^n} + a_{n-1}(x)\frac{d^{n-1}}{dx^{n-1}} + \cdots + a_1(x)</math>
| |
− | で定めると、未知関数 {{mvar|y}} への作用 {{math|''L''({{sfrac|''d''|''dx''}})''y''}} は {{mvar|y}} に関して線型性を持つ。
| |
− | :<math>\begin{align}
| |
− | &L\left(\frac{d}{dx}\right)\left(y_1(x) + y_2(x)\right) = L\left(\frac{d}{dx}\right)y_1(x) + L\left(\frac{d}{dx}\right)y_2(x),\\
| |
− | &L\left(\frac{d}{dx}\right)\lambda y(x) = \lambda L\left(\frac{d}{dx}\right)y(x).
| |
− | \end{align}</math>
| |
− | | |
− | === 1 階連立型 ===
| |
− | 各成分が変数 {{mvar|x}} の(適当な階数の)可微分関数である {{mvar|n}} 次元縦ベクトル {{math|'''y'''(''x'')}}, {{mvar|m}} 次元縦ベクトル {{math|'''b'''(''x'')}} および {{math|''m'' × ''n''}} 行列 {{math|''A''(''x'')}} に対し、
| |
− | :<math>\frac{d\mathbf{y}}{dx} = A(x)\mathbf{y} + \mathbf{b}</math>
| |
− | で定義される微分方程式系を、{{math|''A''(''x'')}} を'''[[係数行列]]''' {{en|(coefficient matrix)}} とする''' 1 階連立型線型微分方程式'''などと呼ぶ。
| |
− | {{math|'''b'''(''x'') {{=}} '''0''' (for all ''x'')}} である場合、方程式は'''斉次'''<ref group="注" name="homogeneous" />であるといい、
| |
− | :<math>\frac{d\mathbf{y}}{dx} = A(x)\mathbf{y}</math>
| |
− | を元の方程式に'''属する斉次方程式'''という。右辺の {{math|''A''(''x'')'''y'''}} は {{math|'''y'''}} に関して線型性を持つ。
| |
− | :<math>\begin{align}
| |
− | &A(x)\left(\mathbf{y}_1(x) + \mathbf{y}_2(x)\right) = A(x)\mathbf{y}_1(x) + A(x)\mathbf{y}_2(x),\\
| |
− | &A(x)\lambda \mathbf{y}(x) = \lambda A(x)\mathbf{y}(x).
| |
− | \end{align}</math>
| |
− | | |
− | 高階単独型線型微分方程式は、変換
| |
− | :<math>y_i := \frac{d^{i-1}y}{dx^{i-1}}, \qquad i = 1, 2,\dots, n.</math>
| |
− | により 1 階連立型の線型微分方程式に変形できる。従って、1 階連立型の線型微分方程式について成り立つ性質は、そのまま高階単独型の線型微分方程式にも適用できる。
| |
− | | |
− | == 解と解空間 ==
| |
− | === 基本解 ===
| |
− | 斉次な線型微分方程式に対し、関数の集合 {{math|''B'' {{=}} {'''y'''<sub>1</sub>(''x''), '''y'''<sub>2</sub>(''x''), ..., '''y'''<sub>''n''</sub>(''x'')}}} がその微分方程式の解空間の[[基底]]となるならば、{{mvar|B}} に属する関数 {{math|'''y'''<sub>''j''</sub>(''x'') (''j'' {{=}} 1, 2, ..., ''n'')}} のことを、その微分方程式の'''基本解'''という。つまり、斉次な線型微分方程式の'''一般解'''はすべて基本解の[[線型結合]]として得られる。また、一般の線型微分方程式では、その方程式の 1 つの特殊解と、その方程式に属する斉次方程式の一般解<ref group="注">つまり基本解の[[線型結合]]</ref>の線型結合が一般解を与える。
| |
− | | |
− | === ロンスキー行列式 ===
| |
− | {{main|ロンスキー行列式}}
| |
− | 斉次方程式の解としていくつかの関数が得られたとき、特に係数行列の形が {{math|''n'' × ''n''}} 成分の[[正方行列]]で、{{mvar|n}} 個の解 {{math|'''y'''<sub>1</sub>(''x''), '''y'''<sub>2</sub>(''x''), ..., '''y'''<sub>''n''</sub>(''x'')}} が得られたとき、それが基本解であるかどうかは次の行列式
| |
− | :<math> W(x) = \begin{vmatrix} | |
− | y_{11}(x) & y_{12}(x) & \cdots & y_{1n}(x) \\
| |
− | y_{21}(x) & y_{22}(x) & \cdots & y_{2n}(x) \\
| |
− | \vdots & \vdots & \ddots & \vdots \\
| |
− | y_{n1}(x) & y_{n2}(x) & \cdots & y_{nn}(x)
| |
− | \end{vmatrix}\quad
| |
− | (\mathbf{y}_j(x) = \begin{pmatrix}
| |
− | y_{1j}(x) \\ y_{2j}(x) \\ \vdots \\ y_{nj}(x)
| |
− | \end{pmatrix})
| |
− | </math>
| |
− | が常に 0 でないことを確認することによって判定できる(実際には任意の 1 点で 0 でないといえば十分である)。
| |
− | | |
− | また、単独高階型の場合には、既に述べた方法でこれを 1 階連立型に帰着すると、解は {{math|'''y'''<sub>''j''</sub> {{=}} (''y''<sub>''j''</sub>, {{sfrac|''dy<sub>j</sub>''|''dx''}}, ..., {{sfrac|''d'' <sup>''n''−1</sup>''y<sub>j</sub>''|''dx'' <sup>''n''−1</sup>}})}} の形で出てくるから、上の行列式は次のように書き換えられる:
| |
− | :<math> W(x) = \begin{vmatrix}
| |
− | y_1 & y_2 & \cdots & y_n \\[5pt]
| |
− | \cfrac{dy_1}{dx} & \cfrac{dy_2}{dx} & \cdots & \cfrac{dy_n}{dx} \\[5pt]
| |
− | \vdots & \vdots & \ddots & \vdots \\[3pt]
| |
− | \cfrac{d^{n-1} y_1}{dx^{n-1}} & \cfrac{d^{n-1} y_2}{dx^{n-1}} & \cdots & \cfrac{d^{n-1} y_n}{dx^{n-1}}
| |
− | \end{vmatrix}.</math>
| |
− | これを'''ロンスキー行列式''' {{en|(Wronski determinant)}} または'''ロンスキアン''' {{en|(Wronskian)}} という。
| |
− | | |
− | == 定数係数の斉次常微分方程式の解法 ==
| |
− | {{mvar|a<sub>k</sub>}} を既知の定数とする斉次線型常微分方程式
| |
− | :<math>\frac{d^n y}{dx^n} + a_{n-1}\frac{d^{n-1}y}{dx^{n-1}} + \cdots + a_0y = 0</math>
| |
− | の左辺に対し、各 {{math|{{sfrac|''d <sup>k</sup>y''|''dx <sup>k</sup>''}}}} を {{mvar|t<sup>k</sup>}} に置き換えて得られる[[多項式]]
| |
− | | |
− | :<math> F(t) = \sum_{k=0}^n a_k t^k = t^n + a_{n-1}t^{n-1} + \dots + a_0</math>
| |
− | | |
− | をこの常微分方程式の'''特性多項式''' {{en|(characteristic polynomial)}}、更に {{mvar|t}} の[[代数方程式]] {{math|''F''(''t'') {{=}} 0}} をこの常微分方程式の'''[[特性方程式]]''' {{en|(characteristic equation)}} という。
| |
− | | |
− | {{mvar|ω}} を代数方程式 {{math|''F''(''t'') {{=}} 0}} の根とすれば、[[指数関数#一般化|指数関数]] {{math|exp(''ωx'')}} は {{math|{{sfrac|''d <sup>k</sup>''exp(''ωx'')|''dx <sup>k</sup>''}} {{=}} ''ω<sup>k</sup>''exp(''ωx'')}} を満たすから、
| |
− | :<math>
| |
− | F\left(\frac{d}{dx}\right){\rm exp}(\omega x)
| |
− | = F(\omega){\rm exp}(\omega x) = 0
| |
− | </math>
| |
− | となり、{{math|''y'' {{=}} exp(''ωx'')}} は元の常微分方程式の解である。ただし、{{math|''f'' ({{sfrac|''d''|''dx''}})}} は、多項式 {{math|''f'' (''t'')}} の {{mvar|t<sup>k</sup>}} を {{math|{{sfrac|''d <sup>k</sup>''|''dx <sup>k</sup>''}}}} に置き換えた[[微分作用素]]である。
| |
− | | |
− | 特性多項式 {{math|''F''(''t'')}} が[[重根]]を持たなければ、[[線型代数学]]でよく知られた事実により集合 {{math|{exp(''ωx'') {{!}} ''ω'' は ''F''(''t'') の根}}} は元の常微分方程式の解を生成する<ref group="注">つまり、基本解になる。</ref>。重根を持つならば {{math|''x''exp(''ωx'')}} などがさらに必要となる。
| |
− | | |
− | == 関数係数の斉次常微分方程式の解法 ==
| |
− | | |
− | 1960年以降の研究で,定数係数ではない関数係数<ref name="sugit-2007ngze" />の斉次常微分方程式の解法が報告されている<ref name="diffequ-2005ngzd" />。
| |
− | | |
− | 主に,求積法による解法が多く、2 階線型常微分方程式をはじめ、多くの非線型常微分方程式がある<ref name="diffequ-2005ngzd" />。
| |
− | これらの中に、一般の陰関数型の常微分方程式があるので、この陰関数型の関数に線型の関数型を与えれば、線型の常微分方程式が得られる。
| |
− | | |
− | 以下に,求積法で解ける主な関数係数の 2 階線型常微分方程式の例を記述する<ref name="diffequ-2005ngzd" /><ref name="diffeq-198205">長島 隆廣[常微分方程式134例とその解]丸善出版サービスセンター,1982年5月発行,国立国会図書館・請求記号 MA117-111,全国書誌番号 82049441</ref>。
| |
− | | |
− | {| class="wikitable" style="text-align:center"
| |
− | |+ 求積法で解ける 2 階線型常微分方程式の例<ref group="注">{{mvar|a}} と {{mvar|b}} および {{mvar|α}} と {{mvar|β}} は定数、{{math|''C''<sub>1</sub>, ''C''<sub>2</sub>}} は[[積分定数]]である。ただし、{{math|''a'' ≠ 0, ''b'' ≠ 0, ''β'' ≠ 0}} を満たす。</ref>
| |
− | ! 常微分方程式 !! 一般解<ref name="diffequ-2005ngzd" />
| |
− | |-
| |
− | | <math>\frac{{d}^2y}{{d}x^2}-xP(x)\frac{{d}y}{{d}x}+P(x)y=0</math>||<math>y=x \Bigl\{C_1 +C_2 \int \frac{1}{\,x^2 \,}\exp \Bigl( \int\! x P(x) \,dx \Bigr)\, dx \Bigr\}</math>
| |
− | |-
| |
− | | <math>\frac{{d}^2y}{{d}x^2}+P(x)\frac{{d}y}{{d}x}-a(a+P(x))y=0</math>||<math>y=e^{ax}\Bigl\{C_1 +C_2\! \int \exp \Bigl(\! -2ax -\!\! \int\! P(x)\,dx \Bigr)\, dx \Bigr\}</math>
| |
− | |-
| |
− | | <math>P(x)\frac{{d}^2 y}{{d}x^2}+(a+bx)\frac{{d}y}{{d}x}-by=0</math>||<math>y= C_1\!\! \int \! \! \int \!\! \frac{1}{\,P(x)\,}\exp \Bigl(\! -\!\! \int \! \frac{\,a+bx\,}{P(x)}\,dx \Bigr)\, dx\, dx +C_2\Bigl(x+\frac{a}{\,b\,}\Bigr)</math>
| |
− | |-
| |
− | | <math>\frac{{d}^2y\,}{{d}x^2}-\left(\frac{1}{2P(x)}\cdot\frac{{d}P(x)}{{d}x}\right)\frac{{d}y}{{d}x}+P(x)y=0</math>||<math>y=C_1\sin\left(\int\sqrt{P(x)}\,{d}x\right)+C_2\cos\left(\int\sqrt{P(x)}\,{d}x\right)</math>
| |
− | |-
| |
− | | <math>\frac{{d}^2y\,}{{d}x^2}-\left(\frac{1}{P(x)}\cdot\frac{{d}P(x)}{{d}x}\right)\frac{{d}y}{{d}x}-\left(P(x)\right)^2 y=0</math>||<math>y=C_1\exp\left(\int P(x)\,{d}x\right)+C_2\exp\left(-\int P(x)\,{d}x\right)</math>
| |
− | |-
| |
− | | <math>x\frac{{d}^2 y\,}{{d}x^2}+(\alpha + \beta x)\frac{{d}y}{{d}x}+\beta y = 0</math>||<math>y=x^{1-\alpha}e^{-\beta x} \left( C_1 \int{}x^{\alpha-2} e^{\beta x}\,{d}x + C_2 \right)</math>
| |
− | |}
| |
− | | |
− | == 注釈 ==
| |
− | {{notelist2}}
| |
− | | |
− | == 参考文献 ==
| |
− | {{reflist|refs=
| |
− | <ref name="sugit-2007ngze">{{cite book|和書
| |
− | |editor=日本数学会
| |
− | |title=岩波・数学辞典
| |
− | |edition=第 4 版
| |
− | |year=2007
| |
− | |publisher=岩波書店
| |
− | |isbn=978-4-00-080309-0
| |
− | }}</ref>
| |
− | <ref name="diffequ-2005ngzd">{{cite book|和書
| |
− | |last=長島
| |
− | |first=隆廣
| |
− | |title=常微分方程式 80 余例とその厳密解
| |
− | |publisher=近代文芸社
| |
− | |year=2005
| |
− | |isbn=4-7733-7282-6. 国立国会図書館蔵書, 請求記号:MA117-H55(東京 本館書庫)
| |
− | }}</ref>
| |
− | }}
| |
− | | |
− | == 関連項目 ==
| |
− | * [[関数解析学]]
| |
− | * [[微分方程式]]
| |
− | ** [[常微分方程式]]
| |
− | ** [[偏微分方程式]]
| |
− | * [[ロンスキー行列式]]
| |
− | * [[求積法]]
| |
− | * [[変数分離]]
| |
− | * [[線型代数学]]
| |
− | ** [[線型性]]
| |
− | ** [[線型方程式]]
| |
− | | |
− | {{DEFAULTSORT:せんけいひふんほうていしき}}
| |
− | | |
− | {{sci-stub}}
| |
− | | |
− | [[Category:微分方程式]]
| |
− | [[Category:線型代数学]]
| |
− | [[Category:解析学]]
| |
− | [[Category:数学に関する記事]]
| |