ニュートンラフソン法
のバックアップソース(No.8)
Unity学習帳2冊目
ニュートンラフソン法
のバックアップソース(No.8)
[
トップ
] [
差分
|
バックアップ
|
リロード
] [
新規
|
一覧
|
検索
|
最新
|
ヘルプ
]
[ ]
差分
を表示
現在との差分
を表示
ニュートンラフソン法
へ行く。
« Prev
Next »
#jsmath **ニュートンラフソン法 [#x75a0b2c] 追加資料:[[はじめMath! Javaでコンピュータ数学_第67回 微分・積分の数学 ニュートン・ラフソン法 [前編]:http://gihyo.jp/dev/serial/01/java-calculation/0067]] ①\(x+{ e }^{ x }=0\) ②\({ x }^{ 2 }=5\) この二つの方程式の未知数\(x\)をニュートンラフソンを用いて解いてみようと思う(この「仕組み」自体は何気にゲームのいろんな場面で利用されている) ニュートンラフソンの式は非常にシンプルで以下のようになっている \({ x }_{ n+1 }={ x }_{ n }-\frac { f\left( { x }_{ n } \right) }{ f'\left( { x }_{ n } \right) } \) この式がどのようにして導かれるのか、はじまりは直線の方程式から来ている 点 \(\left( \alpha ,\beta \right)\) を通る傾き \(m\) の直線の方程式は \(y=m\left( x-\alpha \right) +\beta \) これを変形すると \(y-\beta =m\left( x-\alpha \right) \) 例えば点 \(\left( \alpha =3,\beta =5 \right)\) を通る、傾き \(m=\frac { 1 }{ 2 }\) の直線は \( y=\frac { 1 }{ 2 } \left( x-3 \right) +5\) となる &ref(grp1.png); \(y=\frac { 1 }{ 2 } \left( x-3 \right) +5\quad\) と \(y=0\quad\) との連立方程式を\(x\)に対して解くことにより… \(0=\frac { 1 }{ 2 } \left( x-3 \right) +5\quad \Leftrightarrow \quad \frac { 1 }{ 2 } x-\frac { 3 }{ 2 } +5=0\quad \Leftrightarrow \quad \frac { 1 }{ 2 } x=\frac { 3 }{ 2 } -5\quad \Leftrightarrow \quad \frac { 1 }{ 2 } x=-\frac { 7 }{ 2 } \quad \Leftrightarrow \quad x=-7\) グラフ上の両グラフ線の交点が求まる。これを踏まえたうえで &ref(oneday.png); &font(Red){''この「直線の方程式」を「微積分を利用した関数」で考えると&br;''}; \(y=f\left( x \right) のx=\alpha\) における接線の傾きは \(m=f'\left( \alpha \right) \) 関数 \(y=f\left( x \right)\) 上の点 \(\left( \alpha ,f\left( \alpha \right) \right)\) における接線の方程式は \(y=f'\left( \alpha \right) \left( x-\alpha \right) +f\left( \alpha \right) \) これを変形すると \(y-f\left( \alpha \right) =f'\left( \alpha \right) \left( x-\alpha \right) \quad \quad \cdots ③ \) &font(Red){となる。この③式は「接線の方程式」と呼ばれるもので非常に有名な公式なので微積分と一緒に覚えておいた方が良い&br;}; この式を約分して式変形すると \(\frac { y-f\left( \alpha \right) }{ x-\alpha } =f'\left( \alpha \right) \) となる ここで \(y=f\left( x \right)\) とするならば、この式は \(\frac { f\left( x \right) -f\left( \alpha \right) }{ x-\alpha } =f'\left( \alpha \right) \) となり、この\(x\)を\(b\)に変えると \(\frac { f\left( b \right) -f\left( \alpha \right) }{ b-\alpha } =f'\left( \alpha \right) \) になる さらに、この式に「極限」\(\displaystyle \lim _{ b\rightarrow \alpha }{ } \)を加えると \(\displaystyle f'\left( \alpha \right) =\lim _{ b\rightarrow \alpha }{ \frac { f\left( b \right) -f\left( \alpha \right) }{ b-\alpha } } \) という微分で見慣れた式となる。つまり③式は微分の式を極限なしに変形しているものだという事を強く意識した方が良いという事になっている ここで上記の微分の式を利用して①の問題を解くことを考えてみる。これからやりたい事は以下の動画イメージになる。青い線が求めたい関数が描くグラフ。赤い線が接線の方程式で描いたグラフ 赤い点が青い点に近づいていき、重なった後、またゆっくり離れていく様子に注目してほしい 目的は求めたい関数のy=0とのxに対する連立方程式の解。つまり青い点の\(x\)位置が知りたい &ref(Animation2.gif); これはニュートンラフソンの式を漸化式として数回演算することで求められる \({ x }_{ n+1 }={ x }_{ n }-\frac { f\left( { x }_{ n } \right) }{ f'\left( { x }_{ n } \right) } \) この式は、実は微分の式の変形とy=0との連立方程式から導かれている。上記の微分の式の\(\alpha\)を\(a\)に変えて作ってみる \( その前に以下の違いを確認しておく\\ \\ \lim _{ b\rightarrow a }{ \frac { f\left( b \right) -f\left( a \right) }{ b-a } =f'\left( a \right) } \\ \\ \frac { f\left( b \right) -f\left( a \right) }{ b-a } \simeq f'\left( b \right) \\ \\ わかるだろうか?これは「到達している傾きの位置」が違うという事だ。極限を使うとbはaに限りなく近づくのでf'\left( a \right) となる\\ しかし、極限を利用しないとf'\left( b \right) となる。これはよくよく考えると当たり前のことを言っている\\ でも、これはよくやりやすい勘違いなので注意してほしい。極限なしで平均値の算出で考えを進めると以下になる\\ \\ f'\left( b \right) \simeq \frac { f\left( b \right) -f\left( a \right) }{ b-a } \\ \\ \quad \Leftrightarrow \quad f'\left( b \right) \left( b-a \right) \simeq f\left( b \right) -f\left( a \right) \quad \Leftrightarrow \quad \quad f'\left( b \right) b-f'\left( b \right) a\simeq f\left( b \right) -f\left( a \right) \\ \Leftrightarrow \quad \quad f'\left( b \right) b\simeq f\left( b \right) -f\left( a \right) +f'\left( b \right) a\quad \Leftrightarrow \quad \quad b\simeq \frac { f\left( b \right) -f\left( a \right) +f'\left( b \right) a }{ f'\left( b \right) } \quad \Leftrightarrow \quad \quad b\simeq \frac { f\left( b \right) }{ f'\left( b \right) } -\frac { f\left( a \right) }{ f'\left( b \right) } +a\\ \Leftrightarrow \quad \quad b-\frac { f\left( b \right) }{ f'\left( b \right) } \simeq a-\frac { f\left( a \right) }{ f'\left( b \right) } \\ \\ ここで、f(a)=0との連立方程式をaに対して解くと考える(bがaに近づいていく)\\ \\ b-\frac { f\left( b \right) }{ f'\left( b \right) } \simeq a-\frac { 0 }{ f'\left( a \right) } \quad \Leftrightarrow \quad \quad a\simeq b-\frac { f\left( b \right) }{ f'\left( b \right) } \) &ref(grp.png);
« Prev
Next »
ニュートンラフソン法 のバックアップ一覧
ニュートンラフソン法 のバックアップソース(No. All)
1: 2016-10-06 (木) 22:27:47
osinko
2: 2016-10-07 (金) 00:43:52
osinko
3: 2016-10-07 (金) 01:53:51
osinko
4: 2016-10-07 (金) 09:27:14
osinko
5: 2016-10-07 (金) 13:47:41
osinko
6: 2016-10-24 (月) 23:06:28
osinko
7: 2016-11-09 (水) 02:15:02
osinko
8: 2016-11-09 (水) 09:47:21
osinko
現: 2016-11-09 (水) 18:38:42
osinko