1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109

%indent
////////////////////////////////////////////////////////////////
* 凌宮読取術：$$ \ddd{^2y}{x^2} $ + $ a\ddd{y}{x} $ + $ by $ = $ f $$ ⇒ $$ \begin{cases} \Big(\ddd{}{x} + \alpha\Big) g = f \\ \Big(\ddd{}{x} + \beta\Big) y = g \end{cases} $$ [#of088eec]
 
;,定数係数の２階線形常微分方程式は運動力学や電磁気学で支配方程式として現れるため、
;,「特性方程式((EMANの物理学/物理数学/定数係数2階線形同次微分方程式： http://homepage2.nifty.com/eman/math/differential09.html))
  ＆定数変化法((EMANの物理学/物理数学/定数係数線形非同次微分方程式： http://homepage2.nifty.com/eman/math/differential11.html))」
  または、
  「特性方程式((Matsuda's Web Page/（高専生のための）微分方程式解法ノート/2. 簡単な線形微分方程式>http://www.tsuyama-ct.ac.jp/matsuda/d-eq/bi2.pdf))
  ＆未定係数法((Matsuda's Web Page/（高専生のための）微分方程式解法ノート/3. 線形微分方程式の特殊解の求め方>http://www.tsuyama-ct.ac.jp/matsuda/d-eq/bi3.pdf))」
  など、
;,いわゆる定番な解法が大学入学早々叩き込まれる。
 
;,問題は、この方法は「解のパターンが分かる前提で答えを組み立てる」手法であること
  ((反面、これらの手法は積分を飛ばすために非常に速い。テスト対策のみならず、実用の面においても覚えた方が良い。))、
;,および、場合分けが複雑で一から覚える知識が多いことで、消化不良となりやすい。
;,一応それぞれの背後には大きな理論があるものの、短時間に習得できるものではない。
 
;,これに対し、凌宮数学では、高校で習う因数分解を利用して微分方程式を分解し、
;,直前に習うはずの１階線形常微分方程式に還元して解く方法を使う。
#ceq(e)
    $$ \ddd{^2y}{x^2} $ + $ a $ \ddd{y}{x} $ + $ b $ y $ = $ f $$
#ceq(e)
    ⇒$$ \bigg( $ \ddd{^2}{x^2} $ + $ a\ddd{}{x} $ + $ b $ \bigg) $ y $ = $ f $$
#ceq(a)
    ２階線形常微分演算子化
#ceq(e)
    ⇒ $$ \bigg( $ \ddd{}{x} $ + $ \alpha $ \bigg) $ \bigg( $ \ddd{}{x} $ + $ \beta $ \bigg) $ y $ = $ f $$  
#ceq(a)
    ２つの１階線形常微分演算子に分解
#ceq(e)
    ⇒ $$ \begin{cases} \Big(\ddd{}{x} + \alpha\Big) g = f \\ \Big(\ddd{}{x} + \beta\Big) y = g \end{cases} $$
#ceq(a)
    ２つの１階線形常微分方程式に分離
#ceq(end)
 
%bodynote
////////////////////////////////////////////////////////////////
* 線形常微分演算の因数分解 [#t3af9f02]
最初の読み替えから演算子部を取り出すと、因数分解した形の式が現れる：
#ceq(e)
    $$ \ddd{^2}{x^2} $ + $ a\ddd{}{x} $ + $ b $$ ⇒ $$ \Big(\ddd{}{x} + \alpha\Big) $ \Big(\ddd{}{x} + \beta\Big) $$
#ceq(end)
 
読み替えが可能な理由は常微分の線形性によるが、右側の展開で簡単に確認できる。
#ceq(e)
    $$ \Big(\ddd{}{x} + \alpha\Big) $ \Big(\ddd{}{x} + \beta\Big) $$
#ceq(e)
    ＝ $$ \Big(\ddd{}{x}\Big)^2 $ + $ \alpha $ \ddd{}{x} $ + $ \beta $ \ddd{}{x} $ + $ \alpha $ \beta $$
#ceq(a)
    $$ \beta $$が定数係数であるため、$$ \beta $$と$$ \ddd{}{x} $$は交換可能
    ((逆に言うと、定数と常微分の可換性を利用しているため、$$ \alpha $$も$$ \beta $$も非定数係数の場合は適応できない。))
#ceq(e)
    ＝ $$ \ddd{^2}{x^2} $ + $ (\alpha + \beta) $ \ddd{}{x} $ + $ \alpha $ \beta $$ 
#ceq(e)
    ⇒ $$ \ddd{^2}{x^2} $ + $ a $ \ddd{}{x} $ + $ b $$ 
#ceq(a)
    $$ a $ \equiv $ \alpha $ + $ \beta $$、$$ b $ \equiv $ \alpha $ \beta $$
#ceq(end)
 
//;,ここで、重要なのは$$ a $ = $ \alpha + \beta $$かつ$$ b $ = $ \alpha $ \beta $$という関係であり、
//;,山勘で因数分解が出来ない場合でも、$$ \lambda $$についての二次方程式$$ \alpha^2 $ - $ a $ \alpha $ + b $ = $ 0 $$を解けば良い$$ \alpha $$と$$ \beta $$が求まる
//  ((この$$ \alpha^2 $ - $ a $ \alpha $ + b $ = $ 0 $$は、定番解法の特性方程式$$ \lambda^2 $ + $ a $ \lambda $ + b $ = $ 0 $$とは係数の符号違いである。&br;$$ (\lambda_1 - \alpha) $ (\lambda_2 - \alpha_2) $ = $ 0 $$という関係に対し$$ \lambda $$を求めるか、$$ \alpha $$を求めるかで符号が反転しているだけで、やっていることは同じである。))。
//;,$$ b $ \equiv $ \alpha $ \beta $$より$$ \beta $ = $ \ffd{b}{\alpha} $$。
//  $$ a $ \equiv $ \alpha + \beta $$に代入して$$ a $ = $ \alpha $ + $ \ffd{b}{\alpha} $$。
//  分母を払って$$ a\alpha $ = $ \alpha^2 $ + $ b $$。
//  移行して$$ \alpha^2 $ - $ a $ \alpha $ + $ b $ = 0 $$。
 
%bodynote
////////////////////////////////////////////////////////////////
* 連立１階線形常微分方程式 [#t23de4a5]
;,因数分解形から連立一階形への変形は、$$ \Big(\ddd{}{x} + \beta\Big) $ y $$を$$ g $$と置いただけである。
;,そうすれば、置き換えられた式も、置く$$ g $$の定義式自体も、自ずと１階線形常微分方程式となる。
#ceq(e)
    $$ \Big(\ddd{}{x} + \alpha\Big) $ \iro[ao]{\Big(\ddd{}{x} + \beta\Big)} $ \iro[ao]{y} $ = $ f $$
    　　⇒　　$$ \begin{cases} \Big(\ddd{}{x} + \alpha\Big) \iro[ao]{g} = f \\ \iro[ao]{\Big(\ddd{}{x} + \beta\Big) y} = \iro[ao]{g} \end{cases} $$
#ceq(end)
 
;,実際、連立も名ばかりのもので、上の式に$$ y $$が無いため、先に$$ g $$から解けば良い。
 
* 解の公式 [#a939cf77]
;,１階線形常微分の解より、連立した２式の解はそれぞれ次のようになる：
#ceq(e)
    $$ \Big( $ \ddd{}{x} $ + $ \alpha $ \Big) $ g $ = $ f $$　　⇒　　$$ f $ = $ e^{-{\alpha}x} $ \int $ e^{{\alpha}x} $ g $ dx $$
#ceq(e)
    $$ \Big( $ \ddd{}{x} $ + $ \beta  $ \Big) $ y $ = $ g $$　　⇒　　$$ g $ = $ e^{-{\beta}x} $ \int $ e^{{\beta}x} $ y $ dx $$
#ceq(end)
 
;,下式を上式に代入して$$ g $$を消せば、解の統一公式が得られる：
#ceq(e)
    ''統一公式''： $$ f $ = $ e^{-{\alpha}x} \!\!\int\!\! e^{{\alpha}x} $ \cdot $ e^{-{\beta}x} \!\!\int\!\! e^{{\beta}x} $ \cdot $ y $ \cdot $ dx^2 $$
#ceq(end)
 
;,この公式を解けば、２つの不定積分で計２つの積分定数が現われて２つの基本解を作る。残りの本体が特殊解となる。
;,定番の解法に対する最大の特徴は、場合分けが一切なく、一本の式に統一されている点である。
;,このため、定番解法に比べると、覚える内容が激減する。
 
%bodynote
////////////////////////////////////////////////////////////////
* 適応例【執筆中】 [#b7f8b8ca]
以下に、積和形の公式と定番手法の場合分けの対応例を示す。
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 3 $ y $ = $ 0 $$};（固有値が実数２つ、同次形）>./Dx43y=0]]
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 4 $ y $ = $ 0 $$};（固有値が実数１つ、同次形）>./Dx44y=0]]
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 5 $ y $ = $ 0 $$};（固有値が虚数２つ、同次形）>./Dx45y=0]]
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 3 $ y $ = $ e^{-4x} $$};（固有値が実数２つ、非同次形で非共鳴）>./Dx43y=-4]]
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 3 $ y $ = $ e^{-3x} $$};（固有値が実数２つ、非同次形で共鳴）>./Dx43y=-3]]
- [[&font(black){$$ \ddd{^2y}{x^2} $ + $ 4 $ \ddd{y}{x} $ + $ 4 $ y $ = $ e^{-2x} $$};（固有値が実数１つ、非同次形で共鳴）>./Dx44y=-2]]
 
- [[./@memo]]
    数学 一覧 検索 最新 バックアップ リンク元   ヘルプ   最終更新のRSS