<注意:このページの書いてある事はかなり好き勝手にやってます> memo:常微分方程式 4乗根の実数値を求める漸化式を作る\({ x }^{ 4 }-C\)のような4次方程式の関数のグラフは以下のようになる。このグラフの\(C\)には\(10\)を代入している。この式のグラフは2次方程式の放物線のグラフと似ているが若干頂点付近の曲線が潰れたような形になっている グラフのY軸\(0\)のヨコ線と、この4次方程式の描くグラフの交点が\(\sqrt [ 4 ]{ C }\)、つまり\({ C }^{ \frac { 1 }{ 4 } }\)となる \(\begin{cases} { y=f(x)=x }^{ 4 }-C \\ f'(a)=\lim _{ b\rightarrow a }{ \frac { f(b)-f(a) }{ b-a } } \\ f(a)=0 \end{cases}\) この式を\(b={x}_{0}\)、\(a={x}_{1}\)として定数から変数にして連立させてまとめると \(f'({ x }_{ 0 })=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ \frac { f({ x }_{ 0 })-0 }{ { x }_{ 0 }-{ x }_{ 1 } } } \) となり、これに関数を適用すると \(4{ x }_{ 0 }^{ 3 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ \frac { { x }_{ 0 }^{ 4 }-C-0 }{ { x }_{ 0 }-{ x }_{ 1 } } } \) になる。この式を変形させて漸化式にしていく \( \rightarrow \quad \quad 4{ x }_{ 0 }^{ 3 }({ x }_{ 0 }-{ x }_{ 1 })=\lim _{ { x }_{ 0 } \rightarrow { x }_{ 1 } }{ { x }_{ 0 }^{ 4 }-C } \\ \rightarrow \quad \quad 4{ x }_{ 0 }^{ 4 }-4{ x }_{ 0 }^{ 3 }\cdot { x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ { x }_{ 0 }^{ 4 }-C } \\ \rightarrow \quad \quad -4{ x }_{ 0 }^{ 3 }\cdot { x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ -3{ x }_{ 0 }^{ 4 }-C } \\ \rightarrow \quad \quad 4{ x }_{ 0 }^{ 3 }\cdot { x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ 3{ x }_{ 0 }^{ 4 }+C } \\ \rightarrow \quad \quad 4{ x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ \frac { 3{ x }_{ 0 }^{ 4 }+C }{ { x }_{ 0 }^{ 3 } } } \\ \rightarrow \quad \quad 4{ x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ 3{ x }_{ 0 }+\frac { C }{ { x }_{ 0 }^{ 3 } } } \\ \rightarrow \quad \quad { x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ \frac { 3 }{ 4 } { x }_{ 0 }+\frac { C }{ 4{ x }_{ 0 }^{ 3 } } } \\ \rightarrow \quad \quad { x }_{ 1 }=\lim _{ { x }_{ 0 }\rightarrow { x }_{ 1 } }{ \frac { 1 }{ 4 } \left( { 3x }_{ 0 }+\frac { C }{ { x }_{ 0 }^{ 3 } } \right) } \\ \rightarrow \quad \quad { x }_{ n+1 }=\frac { 1 }{ 4 } \left( { 3x }_{ n }+\frac { C }{ { x }_{ n }^{ 3 } } \right) \) これをC#のコードにすると以下になる(このコードではCに10を代入して4乗根を計算している)
結果出力された値は 1.778279 となり、これは\(\sqrt [ 4 ]{ 10 } \simeq 1.778279\)を表しているので漸化式が正常に機能している事がわかる memo: 本来、計算機で計算するなら元の関数(元の関数は導関数を積分したものだから)から逆算した方が早い 微積分のメリットは例えるなら「速度の情報しかない時に位置を想定して時間を計算できる」という部分にある 累乗根や対数など計算機に任せられる部分があるなら積分して位置の関数から値を計算する方が良い 物理的な値を求める漸化式を作ってみる思考実験として物理的な値を求める漸化式を作る <問題> この問題を解くため微分積分にあてはめると \(v=9.8t\) を形式的積分の式である \(a{ x }^{ n }\quad \mapsto \quad a\frac { 1 }{ n+1 } { x }^{ n+1 }\) に通すと \(\begin{cases} y=4.9{ t }^{ 2 } \\ y=700 \end{cases}\quad \rightarrow \quad 4.9{ t }^{ 2 }=700\quad \rightarrow \quad { t }=\sqrt { \frac { 700 }{ 4.9 } } =11.952...\)・・・① 11.952秒後に通り過ぎる。この計算をする時、電卓に\(\sqrt { }\)キーや\(\log { } \)キーがない。対数表も持ってない \(f'(a)=\lim _{ b\rightarrow a }{ \frac { f(b)-f(a) }{ b-a } } \quad \mapsto \quad 9.8{ t }_{ 0 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 4.9{ t }_{ 0 }^{ 2 }-700 }{ { t }_{ 0 }-{ t }_{ 1 } } } \) \(9.8{ t }_{ 0 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 4.9{ t }_{ 0 }^{ 2 }-700 }{ { t }_{ 0 }-{ t }_{ 1 } } } \quad \rightarrow \quad 9.8{ t }_{ 0 }\left( { t }_{ 0 }-{ t }_{ 1 } \right) =\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ 4.9{ t }_{ 0 }^{ 2 }-700 } \quad \rightarrow \quad 9.8{ t }_{ 0 }^{ 2 }-9.8{ t }_{ 0 }\cdot { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ 4.9{ t }_{ 0 }^{ 2 }-700 } \\ \) \(\rightarrow \quad -9.8{ t }_{ 0 }\cdot { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ -4.9{ t }_{ 0 }^{ 2 }-700 } \quad \rightarrow \quad 9.8{ t }_{ 0 }\cdot { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ 4.9{ t }_{ 0 }^{ 2 }+700 } \quad \rightarrow \quad { t }_{ 0 }\cdot { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 4.9 }{ 9.8 } { t }_{ 0 }^{ 2 }+\frac { 700 }{ 9.8 } } \\ \) \(\rightarrow \quad { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 4.9 }{ 9.8 } { t }_{ 0 }+\frac { 700 }{ 9.8{ t }_{ 0 } } } \quad \rightarrow \quad { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 1 }{ 2 } { t }_{ 0 }+\frac { 700 }{ 9.8{ t }_{ 0 } } } \quad \rightarrow \quad { t }_{ 1 }=\lim _{ { t }_{ 0 } \rightarrow { t }_{ 1 } }{ \frac { 1 }{ 2 } \left( { t }_{ 0 }+\frac { 700 }{ 4.9{ t }_{ 0 } } \right) } \\ \) \(\rightarrow \quad { t }_{ n+1 }=\frac { 1 }{ 2 } \left( { t }_{ n }+\frac { 700 }{ 4.9{ t }_{ n } } \right) \\ \)・・・② この最期の漸化式をコードにすると以下になる
t0の値はどんな値を入れても良い。このコードでは5000回計算を繰り返しているが、微分も積分も無限回数この計算を繰り返している事になる。その答えは必ず、\(11.952...\)へと収束する。微分積分はこの仕組みを利用してコンパクトに値を求めているらしい |