こんばんは、satoです。
目次 非表示
はじめに…
この記事では専門的な数学の話をします。
どちらかといえば「数学に興味がある」とか「数学を勉強している」人向けになります。
その辺りをご了承くださいませ(;’∀’)
今日の話は「微分方程式をどうやって解くのか?」ということについて、よく書かれている公式がどうして成り立つのかを考えてみたことをまとめてみました。
その中で「接続」という概念も出てくるので、その解説をします。
微分方程式の解法
たとえば、次の微分方程式を考えます。
このように方程式の中で高々2階微分で表されるものを2階微分方程式と言います。
2階微分すると
![Rendered by QuickLaTeX.com -1](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-7b34c01098c83fa602de54e9d74d63a9_l3.png)
これ、答えは
![Rendered by QuickLaTeX.com \cos x, \sin x](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-089c0a1fad09c0d30288c2e01b90adcf_l3.png)
教科書によくある解法
大学一年生とかで習う解法は次のようになります。
まず方程式をとすると、
を
に置き換えてできる
を解くと
となります。
これを使うと微分方程式の解がとなります。ただし、
は任意の定数です。
ここで、オイラーの公式を用いると
となり、
![Rendered by QuickLaTeX.com \cos x, \sin x](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-089c0a1fad09c0d30288c2e01b90adcf_l3.png)
この解法の中で一つ方程式を解いたところがありますが、これを「特性方程式」と言います。
特性方程式を解いたら、その解から微分方程式の解が解けてしまうのです。
もう少し詳しく
ところで、どうしてこのように解答できるのでしょうか?それは次のような計算をしているからです。
最後の方程式は「微分するとが出てくる関数」ということで簡単に求められます。
特性方程式は「微分方程式を因数分解」するために使ったわけです。
また、最後の話は「重ね合わせの原理」というものを使っています。先ほどの簡単な微分方程式を満たす二つの解はどちらも、元々の微分方程式の解を満たします。
ここから、二つの解を足し合わせたものも解となるというのが重ね合わせの原理なのです。
特性方程式が重解を持つ場合
こうして、2階微分方程式を解けるのですが、この解法には問題が一つありまして…特性方程式が重解を持つ時にうまく行かないということです。
の特性方程式は
![Rendered by QuickLaTeX.com \lambda^2-2\lambda+1=0](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-d0a7900c91e8910373e0d507706f2c00_l3.png)
![Rendered by QuickLaTeX.com \lambda=1](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-3fad11e39fc84afed3ac93a27d13fbaa_l3.png)
と一つの解しか出てきません。こういうとき教科書には二つ目の解として
![Rendered by QuickLaTeX.com y=C_2xe^x](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-1892744c0b5b1928549ef7f40739b1ba_l3.png)
が…どうしてこの解が出てきたのでしょうか?
疑問
そもそも、特性方程式がどうして出てきたのか?特性方程式と微分方程式にはどういう関係があるのか?
また、重ね合わせの原理がどうして成り立つのか?
など、いろいろ疑問が出てきます。今日はこれについて自分なりに理解したことをまとめてみます。
行列の話が分かればある程度は分かるはず。大学一年生の線形代数の知識を前提としています。
微分方程式を行列で表す
というのが、これから話す解き方のコンセプトです。
記号の解説
二つの関数があったときに、
というように2つの関数によるベクトルを考えます。そして、このベクトルに関して
というようなベクトルに対する作用素を定義します。ようするに
![Rendered by QuickLaTeX.com \nabla](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-d8b9de357424b30554f6eb9a1c65efb7_l3.png)
を行列で表す
さて、今
を満たすような関数
![Rendered by QuickLaTeX.com y](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-0af556714940c351c933bba8cf840796_l3.png)
![Rendered by QuickLaTeX.com \displaystyle\begin{pmatrix}y\\y'\end{pmatrix}](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-7a638759b133a3f2be81f9b977b6ab86_l3.png)
![Rendered by QuickLaTeX.com \nabla](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-d8b9de357424b30554f6eb9a1c65efb7_l3.png)
微分方程式から
が成り立ちます。最初は当たり前、二番目の式は微分方程式の言いかえです。
これを使うと
となります。
ここまで何をしたかと言いますと微分方程式を満たす関数とその導関数の組ががどのように変換されるかを行列で表現したというわけです。
行列と言えば…
大学一年生の時に「固有値」を嫌というほど計算させられたと思います。というわけで、さっき出てきた行列(とします)の固有値を求めましょう。
なので、
となります。(
は単位行列)
ここで「特性方程式」が出てきましたね!
さて、ついでに固有ベクトルを求めますとのときは
、
![Rendered by QuickLaTeX.com \lambda=-i](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-c5456e00d0a1458dbdbf11d26c44ee9d_l3.png)
となります。ここで
と置くと、
が成り立ちます。対角化ですね。ここから
となります。
微分方程式を変形する
さて、上を用いると
となります。ここで、ベクトル
![Rendered by QuickLaTeX.com v](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-ef71511c70f0e4b25cc6bd69f3bc20c2_l3.png)
とすると、
となります。つまり、
となるのです!この方程式から
が得られるので、あとは
![Rendered by QuickLaTeX.com P](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-650eb7688af6737ac325425b5c9a5982_l3.png)
![Rendered by QuickLaTeX.com y=C_1e^{ix}+C_2e^{-ix}](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-9fc574e7077a6fbea382a43245378406_l3.png)
いやー長かったですね…。
で、結局何をしたの?
先ほどの解法に比べると分かりにくいし面倒くさい気がするこの「微分方程式を行列に置き換える」という解法ですが、結局何をしたのでしょうか?まとめると
与えられた微分方程式の解とその導関数の組によるベクトルを、という作用がどのように変換するのかを考えます。
そうして、という変換によって動かない固有ベクトルを求めます。すると、この組に対しては微分方程式が簡単になります。
その簡単な微分方程式を解くことで元の微分方程式の解を得ました。
という流れです。で、このというのが接続と言われているものです。微分方程式とその導関数の組は一つの「幾何学」的なものがあるのですが、その間の変換を表しています。
で、その行列での表現を「接続行列」と言ったりします。
長いようですが、一つ一つの変形がどのようになっているのかをより深く理解できます。
特性方程式は接続行列の固有方程式
微分方程式の因数分解の正体は固有ベクトルによるでの変換
重ね合わせの原理は固有ベクトルから元の微分方程式の解を出す計算
に対応しているわけです。
さらに、特性方程式の解が重解の時にどうして二つ目の解がああなるのかについても説明できます。
特性方程式が重解を持つ場合
上と同様に
とすると、微分方程式の解とその導関数の組からなるベクトルに対して
となります。ここで出てきた行列を
![Rendered by QuickLaTeX.com B](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-770fd1447ccf2fc229801b486b0d8f8a_l3.png)
![Rendered by QuickLaTeX.com B](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-770fd1447ccf2fc229801b486b0d8f8a_l3.png)
![Rendered by QuickLaTeX.com (\lambda-1)^2=0](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-1fa4a3047facbf81f786dc0a1717b335_l3.png)
![Rendered by QuickLaTeX.com 1](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-4868771cbc422b5818f85500909ce433_l3.png)
この固有ベクトルを求めると
の一種類だけになります。あれ、これだと対角化できないぞ…(;’∀’)
広義固有ベクトル
そこで、広義固有ベクトルというものを使います。固有ベクトルはとなるベクトルですが、広義な固有ベクトルというのは適当な自然数
に対して
となるベクトルを表します。
求めるときには固有ベクトルを使って
となる
を求めて、あとは順次…という流れです。
この場合は計算すると、たとえば
が当てはまります。さらに、
とすると、次が成り立ちます。
このように、対角化できない行列も上三角行列に変えることができます。
微分方程式を変形
これを用いると
が成り立ちます。よって、ベクトル
![Rendered by QuickLaTeX.com u](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-43fe27dc3e528266a619764d90fce60b_l3.png)
とすると、
となります。
これで、
という微分方程式が得られました。このうち二番目はすぐに計算できて
![Rendered by QuickLaTeX.com u_2=C_2e^x](https://setsuri-nihon.net/wp-content/ql-cache/quicklatex.com-448068b6454ede4d97581c76c4bd36bc_l3.png)
非線形微分方程式
さて、これを使うと一つ目の微分方程式が
となりますが、これは実を言いますと非線形微分方程式というものになります。今まで扱っていた「線形微分方程式」と違って、これを解くことは簡単ではありません。
この場合はが解になります。これを「特殊解」と言います。実際にこれを求めるときにはラプラス変換とかを使えばいいのですが…私はまだ理解できてません(´・ω・`)
というわけで、一つ目の微分方程式の解はとなりました。
であることからこれが元の微分方程式の解となるわけです。
つまり、重解の時に出てくるもう一つの解の正体はこの非線形微分方程式の特殊解だったのでした。
まとめ
というわけで、2階線形微分方程式の解法を詳しく考えてみた結果を書きました。
まとめると、微分方程式を「接続」という概念を通して行列に置き換えて、固有方程式を解くことで固有値を求め簡単な形に帰着させるということでした。
ちなみに、微分方程式とその導関数の組は一つの「幾何学」的なものと書いた話は正確に表現しようとすると「-加群」の話に繋がります。
とは微分作用素環のことです。
ただ数式を当てはめるだけで解ける微分方程式。
だけど、その解法の奥底には線型代数の理論を駆使した世界、そして幾何学の世界があるのです。
このように、御言葉も奥深い世界があります。単に言葉だけ見るとなにも分からないですが、その奥にある霊の世界を見ると…とても深い次元の言葉だとわかるのです。
…初めて数式を交えた専門的な記事を書きましたが、とても長くなりましたね(;’∀’)ここまで読んでいただきありがとうございました。
今後もこういう研究メモを書いてみようと思うので、数学好きの方はチェックしてみてください(*‘∀‘)
この記事を書いたブロガー
-
「素直に、深く、面白く」がモットーの摂理男子。霊肉ともに生粋の道産子。30代になりました。目指せ数学者。数学というフィールドを中心に教育界隈で色々しています。
軽度の発達障害(ADHD・PD)&HSP傾向あり。
最近の投稿
日常生活2024年7月26日自信って大事だなと実感した一週間
日常生活2024年7月25日天でも地でも求められる「誠実さ」
日常生活2024年7月24日どうして、私に聖書が、御言葉が必要なのか?
日常生活2024年7月23日摂理人の『明け方』を巡る戦争。