こんばんは、satoです。
目次 非表示
はじめに…
この記事では専門的な数学の話をします。
どちらかといえば「数学に興味がある」とか「数学を勉強している」人向けになります。
その辺りをご了承くださいませ(;’∀’)
今日の話は「微分方程式をどうやって解くのか？」ということについて、よく書かれている公式がどうして成り立つのかを考えてみたことをまとめてみました。
その中で「接続」という概念も出てくるので、その解説をします。
微分方程式の解法
たとえば、次の微分方程式を考えます。
![\[\frac{d^2y}{dx^2}=-y\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-b83f7a85c10ea787e9559c5b593e0b1f_l3.png "Rendered by QuickLaTeX.com")
このように方程式の中で高々2階微分で表されるものを2階微分方程式と言います。
2階微分すると
これ、答えは
教科書によくある解法
大学一年生とかで習う解法は次のようになります。
まず方程式を
これを使うと微分方程式の解が
ここで、オイラーの公式
![\[\frac{e^{ix}+e^{-ix}}{2}=\cos x, \frac{e^{ix}-e^{-ix}}{2i}=\sin x\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-f9cba604958b70bbd7d5f2fd373c0683_l3.png "Rendered by QuickLaTeX.com")
となり、
この解法の中で一つ方程式を解いたところがありますが、これを「特性方程式」と言います。
特性方程式を解いたら、その解から微分方程式の解が解けてしまうのです。
もう少し詳しく
ところで、どうしてこのように解答できるのでしょうか？それは次のような計算をしているからです。
最後の方程式は「微分すると
特性方程式は「微分方程式を因数分解」するために使ったわけです。
また、最後の話は「重ね合わせの原理」というものを使っています。先ほどの簡単な微分方程式を満たす二つの解はどちらも、元々の微分方程式の解を満たします。
ここから、二つの解を足し合わせたものも解となるというのが重ね合わせの原理なのです。
特性方程式が重解を持つ場合
こうして、2階微分方程式を解けるのですが、この解法には問題が一つありまして…特性方程式が重解を持つ時にうまく行かないということです。
![\[\frac{d^2y}{dx^2}-2\frac{dy}{dx}+y=0\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-5781984667349e0c3f6421bc16104909_l3.png "Rendered by QuickLaTeX.com")
の特性方程式は
![\[(\frac{d}{dx}-1)y=0\Leftrightarrow y=C_1e^{x}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-79a120bd3f84d7613f79df3fd4ccfd66_l3.png "Rendered by QuickLaTeX.com")
と一つの解しか出てきません。こういうとき教科書には二つ目の解として
が…どうしてこの解が出てきたのでしょうか？
疑問
そもそも、特性方程式がどうして出てきたのか？特性方程式と微分方程式にはどういう関係があるのか？
また、重ね合わせの原理がどうして成り立つのか？
など、いろいろ疑問が出てきます。今日はこれについて自分なりに理解したことをまとめてみます。
行列の話が分かればある程度は分かるはず。大学一年生の線形代数の知識を前提としています。
微分方程式を行列で表す
というのが、これから話す解き方のコンセプトです。
記号の解説
二つの関数
![\[\begin{pmatrix}f \\g\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-0f8823a11adc58e0bfe2cf65f31cd515_l3.png "Rendered by QuickLaTeX.com")
というように2つの関数によるベクトルを考えます。そして、このベクトルに関して
![\[\nabla\begin{pmatrix}f \\g\end{pmatrix}=\begin{pmatrix}\frac{df}{dx} \\\frac{dg}{dx}\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-b04b947a4a26b7218f3aa867838a5cec_l3.png "Rendered by QuickLaTeX.com")
というようなベクトルに対する作用素を定義します。ようするに
さて、今
![\[\frac{d^2y}{dx^2}+y=0\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-6f2010a75c5b40ae7834d71544ee699f_l3.png "Rendered by QuickLaTeX.com")
を満たすような関数
微分方程式から
![\[\frac{dy}{dx}=y', \frac{dy'}{dx}=-y\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-942af2cb8c7008cb42a194c084419406_l3.png "Rendered by QuickLaTeX.com")
が成り立ちます。最初は当たり前、二番目の式は微分方程式の言いかえです。
これを使うと
![\[\nabla\begin{pmatrix}y \\y'\end{pmatrix}=\begin{pmatrix}\frac{dy}{dx} \\\frac{dy'}{dx}\end{pmatrix}=\begin{pmatrix}0&1 \\-1&0\end{pmatrix}\begin{pmatrix}y\\y'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-ba1f01cd1d4c2317b34c3a8105321f78_l3.png "Rendered by QuickLaTeX.com")
となります。
ここまで何をしたかと言いますと微分方程式を満たす関数とその導関数の組が
行列と言えば…
大学一年生の時に「固有値」を嫌というほど計算させられたと思います。というわけで、さっき出てきた行列(
ここで「特性方程式」が出てきましたね！
さて、ついでに固有ベクトルを求めますと
![\[c\begin{pmatrix}1\\i\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-7feafe40b2c9c842201968073b63b854_l3.png "Rendered by QuickLaTeX.com")
、
![\[c\begin{pmatrix}1\\-i\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-94d6050d5e0ae66f150828384dc5a305_l3.png "Rendered by QuickLaTeX.com")
となります。ここで
![\[P=\begin{pmatrix}1&1\\i&-i\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-20af6b0150d5d3b26871334650320cc1_l3.png "Rendered by QuickLaTeX.com")
と置くと、
![\[P^{-1}AP=\begin{pmatrix}i&0\\0&-i\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-88907c4c9e13f3610b68adb615fdd49b_l3.png "Rendered by QuickLaTeX.com")
が成り立ちます。対角化ですね。ここから
![\[AP=PP^{-1}AP=P\begin{pmatrix}i&0\\0&-i\end{pmatrix},\ A=APP^{-1}=P\begin{pmatrix}i&0\\0&-i\end{pmatrix}P^{-1}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-b45b74ac9bac77c8e77ff09e5b7184e0_l3.png "Rendered by QuickLaTeX.com")
となります。
微分方程式を変形する
さて、上を用いると
![\[\nabla\begin{pmatrix}y \\y'\end{pmatrix}=A\begin{pmatrix}y\\y'\end{pmatrix}=P\begin{pmatrix}i&0\\0&-i\end{pmatrix}P^{-1}\begin{pmatrix}y\\y'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-2aba2e14b27588411b1a129d8c452bca_l3.png "Rendered by QuickLaTeX.com")
となります。ここで、ベクトル
![\[v=\begin{pmatrix}v_1\\v_2\end{pmatrix}=P^{-1}\begin{pmatrix}y\\y'\end{pmatrix}=\begin{pmatrix}y+y'\\iy-iy'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-7a2dd25d4bb0189fef795e5abbf6e962_l3.png "Rendered by QuickLaTeX.com")
とすると、
![\[\nabla v=\begin{pmatrix}i&0\\0&-i\end{pmatrix}v\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-2fc580288fb9dde0a82e638edb007239_l3.png "Rendered by QuickLaTeX.com")
となります。つまり、
![\[\frac{d}{dx}v_1=iv_1, \frac{d}{dx}v_2=-iv_2\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-c2ddd6490173101efb9090c5d9448b0d_l3.png "Rendered by QuickLaTeX.com")
となるのです！この方程式から
![\[v_1=C_1e^{ix}, v_2=C_2e^{-ix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-c889b245355b0743144478e71dd342e6_l3.png "Rendered by QuickLaTeX.com")
が得られるので、あとは
いやー長かったですね…。
で、結局何をしたの？
先ほどの解法に比べると分かりにくいし面倒くさい気がするこの「微分方程式を行列に置き換える」という解法ですが、結局何をしたのでしょうか？まとめると
与えられた微分方程式の解とその導関数の組によるベクトルを、
そうして、
その簡単な微分方程式を解くことで元の微分方程式の解を得ました。
という流れです。で、この
で、その行列での表現を「接続行列」と言ったりします。
長いようですが、一つ一つの変形がどのようになっているのかをより深く理解できます。
特性方程式は接続行列の固有方程式
微分方程式の因数分解の正体は固有ベクトルによる
重ね合わせの原理は固有ベクトルから元の微分方程式の解を出す計算
に対応しているわけです。
さらに、特性方程式の解が重解の時にどうして二つ目の解がああなるのかについても説明できます。
特性方程式が重解を持つ場合
上と同様に
![\[\nabla\begin{pmatrix}f \\g\end{pmatrix}=\begin{pmatrix}\frac{df}{dx} \\\frac{dg}{dx}\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-a569dcfb031f9a1867df3a4dcc4b44b9_l3.png "Rendered by QuickLaTeX.com")
とすると、微分方程式の解とその導関数の組からなるベクトルに対して
![\[\nabla\begin{pmatrix}y \\y'\end{pmatrix}=\begin{pmatrix}0&1 \\-1&2\end{pmatrix}\begin{pmatrix}y\\y'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-c7d9eea9a199a33eb758b878313d7604_l3.png "Rendered by QuickLaTeX.com")
となります。ここで出てきた行列を
この固有ベクトルを求めると
![\[c\begin{pmatrix}1\\1\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-d1fbe4799d4587686ac55a6ee5964f09_l3.png "Rendered by QuickLaTeX.com")
の一種類だけになります。あれ、これだと対角化できないぞ…(;’∀’)
広義固有ベクトル
そこで、広義固有ベクトルというものを使います。固有ベクトルは
求めるときには固有ベクトル
この場合は計算すると、たとえば
![\[\begin{pmatrix}0 \\1\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-db89bb311c2cc9ac7b0caf3ad289fd38_l3.png "Rendered by QuickLaTeX.com")
が当てはまります。さらに、
![\[Q=\begin{pmatrix}1&0 \\1&1\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-24ffc16b668077d808591941667770f3_l3.png "Rendered by QuickLaTeX.com")
とすると、次が成り立ちます。
![\[Q^{-1}BQ=\begin{pmatrix}1&1 \\0&1\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-8f009a3ad537ab68734aba7f5c94df0f_l3.png "Rendered by QuickLaTeX.com")
このように、対角化できない行列も上三角行列に変えることができます。
微分方程式を変形
これを用いると
![\[\nabla\begin{pmatrix}y \\y'\end{pmatrix}=Q\begin{pmatrix}1&1 \\0&1\end{pmatrix}Q^{-1}\begin{pmatrix}y\\y'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-d0e450d0b6250df19e18ebb2d89044b2_l3.png "Rendered by QuickLaTeX.com")
が成り立ちます。よって、ベクトル
![\[u=\begin{pmatrix}u_1\\u_2\end{pmatrix}=Q^{-1}\begin{pmatrix}y\\y'\end{pmatrix}=\begin{pmatrix}y\\-y+y'\end{pmatrix}\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-9d13b61f98709448867000f457c4a9b1_l3.png "Rendered by QuickLaTeX.com")
とすると、
![\[\nabla u=\begin{pmatrix}1&1 \\0&1\end{pmatrix}u\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-1a4466ad0b4deba52f64314c11fe2275_l3.png "Rendered by QuickLaTeX.com")
となります。
これで、
![\[\frac{d}{dx}u_1=u_1+u_2, \frac{d}{dx}u_2=u_2\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-158e4b823e6c79c0d53f77b883c97fd4_l3.png "Rendered by QuickLaTeX.com")
という微分方程式が得られました。このうち二番目はすぐに計算できて
非線形微分方程式
さて、これを使うと一つ目の微分方程式が
![\[\frac{d}{dx}u_1=u_1+C_2e^x\]](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-9533496af10258cb6e505bb404702635_l3.png "Rendered by QuickLaTeX.com")
となりますが、これは実を言いますと非線形微分方程式というものになります。今まで扱っていた「線形微分方程式」と違って、これを解くことは簡単ではありません。
この場合は
というわけで、一つ目の微分方程式の解は
つまり、重解の時に出てくるもう一つの解の正体はこの非線形微分方程式の特殊解だったのでした。
まとめ
というわけで、2階線形微分方程式の解法を詳しく考えてみた結果を書きました。
まとめると、微分方程式を「接続」という概念を通して行列に置き換えて、固有方程式を解くことで固有値を求め簡単な形に帰着させるということでした。
ちなみに、微分方程式とその導関数の組は一つの「幾何学」的なものと書いた話は正確に表現しようとすると「
ただ数式を当てはめるだけで解ける微分方程式。
だけど、その解法の奥底には線型代数の理論を駆使した世界、そして幾何学の世界があるのです。
このように、御言葉も奥深い世界があります。単に言葉だけ見るとなにも分からないですが、その奥にある霊の世界を見ると…とても深い次元の言葉だとわかるのです。
…初めて数式を交えた専門的な記事を書きましたが、とても長くなりましたね(;’∀’)ここまで読んでいただきありがとうございました。
今後もこういう研究メモを書いてみようと思うので、数学好きの方はチェックしてみてください(*‘∀‘)
この記事を書いたブロガー
-
「素直に、深く、面白く」がモットーの摂理男子。霊肉ともに生粋の道産子。30代になりました。目指せ数学者。数学というフィールドを中心に教育界隈で色々しています。
軽度の発達障害(ADHD・PD)&HSP傾向あり。
最近の投稿
