誰かﾏｸﾛｰﾘﾝ展開、教えて!
■ このスレッドは過去ログ倉庫に格納されています
0001うー
NGNGマクローリン展開。
「
世の中凄いことを考える人がいるもので、ありとあらゆる関数が
この形で展開できるというのです。
f(x)=a0 + a1x + a2x^2 + a3x^3 + a4x^4 ････
」
ってあって、上記が正しいこと前提に最後までつっぱしてるんで、納得いかないっす。
ネットでも検索したし、参考書も見たけど、どれも似たように、なぜか証明せず
それ前提で、次の微分へと進んじゃいます。（そもそも無限回数微分可能な関数のみが
展開できるのであって、ありとあらゆる関数ってわけじゃないですよね？）
ちなみに、その後、ｎ回微分して各係数を求める行為そのものは理解できるっす。（教科書通りだし）
0002デフォルトの名無しさん
NGNG￣￣￣￣￣∨￣￣￣ (´´
∧∧ ） (´⌒(´
⊂（ﾟДﾟ⊂⌒｀つ≡≡≡(´⌒;;;≡≡≡
￣￣ (´⌒(´⌒;;
ｽﾞｻﾞｰｰｰｰｰｯ
0003デフォルトの名無しさん
NGNG￣￣￣￣￣∨￣￣￣
,〜（（((((((〜〜､
( ＿（（（((((((＿ )
|/ ~^^＼)／^^~ヽ|
| ＿ 《 ＿ |
(|-(_//_)-(_//_)-|）
| 厶､ |
＼ |||||||||||| ／
⊂ ＼_______⊂ ⌒`つ≡≡≡(´⌒;;;≡≡≡
￣ (´⌒(´⌒;;
0004デフォルトの名無しさん
NGNG| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 数学板でなくてよかったな
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0005デフォルトの名無しさん
NGNG一件だけヒット。↓
http://www.gem.hi-ho.ne.jp/katsu-san/audio/ft_lt.htm
と思ったら
「（これの証明は、マクローリン展開してやれば
出来るのですが、そのマクローリン展開の証明は
かなり長くなるので割愛します。） 」だってさ。（欝）
思うに
マクローリン展開がなぜできる？に関する疑問に応えられる人は
いないんじゃない？盲目的に、次の説明へと進んでたよ、俺の大学の数学の授業。
まぁ、テストで点とれりゃいーんだけどね。
0006デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ まともな条件付でテイラーの定理の帰結だ
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0007デフォルトの名無しさん
NGNG0008デフォルトの名無しさん
NGNGそのものの説明をすっとばしてる。プロセス解説が目的なのに
プロセス抜かして結果だけ載せてる参考書って多いやねぇ・・・
とりあえず、「多項式で表現できるんだ」だから「微分して近似値を
求められるんだ」と頭ごなしに教えられたよ、漏れも。
実際、sin() cos() のソフト演算って、結果さえ知ってれば
大学生でも作れるからねぇ。
0009デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 無限回可微分とテイラー展開可能じゃ全然違うぞ
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0010デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 近似したやつが中心から外れると近似になる保証も無いし
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0011デフォルトの名無しさん
NGNG0012デフォルトの名無しさん
NGNG剰余項による収束の評価でよければ検索すれば
チャカチャカ出てくると思う。
キーワード｢テイラー展開｣｢剰余項｣で検索してみなよ。
001312
NGNG0014デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 無限回可微分とテイラー展開可能とローラン展開可能じゃ全然違うぞ
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0015デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 言いたいことはわからんでもないが
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0017デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 要は、複素数まで広げないと( ﾟДﾟ)ﾊｧ?ってこと？
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
001912
NGNG“テイラー展開では”ではないけどね
0020デフォルトの名無しさん
NGNGなどの多項式に展開できることを図でイメージできないよ、俺。
そもそも、中心からずれれば、誤差が出るんだよね？10のもなーが
言ってるみたいに。まぁ、sin()だとかcos()は幸いにして周期関数
だから、0 から 2πの間に丸め込めるから深刻な問題ではないけど。
（CPUも、内部ではマクローリン展開ではじき出してるので、cos(x) と
cos(x+n*π)とでは、誤差があると明記されてるしね。）
002120
NGNG0022デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ fmodすらしないで露骨に展開してるんですか？
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0023デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ よく読んだら周期の境界のお話でした
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0024げーむ屋
NGNGVC++付属の（というか多くの標準数学ライブラリ）の△関数関連内部
に存在するfmodは無意味であるばかりか、足かせとなります。
いっぽう、某PS2の実行部隊1のチップ内にある△関数命令は
fmodかまさないばかりか、2*π の付近の時点で精度がかなりくたばってます。
Ｘ箱万歳。
0025デフォルトの名無しさん
NGNG曰
| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ >>24 CAD屋ですが、勉強になります
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0026デフォルトの名無しさん
NGNG実際にグラフにして、イメージ掴むしかないか・・・
でも、分母が階乗だから、すぐに小さくなっちゃうんだよね−？
0027デフォルトの名無しさん
NGNG| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ x^nの係数はｆ^(n)(0)/n!だから1/(1-x)とかはそうでもないよ
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||代||Ｌ二⊃ .￣￣＼
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0028デフォルトの名無しさん
NGNG「解析概論」,高木貞治著,岩波出版,ISBN4-00-005171-7,\2600
をお勧めする。
それのP.213によれば、
1)「解析函数は、それが正則なる領域内の任意の点において
Taylor級数に展開される。」
なお、「複素平面上の或る領域Kの各点において微分可能な函数を
Kにおいて正則な解析函数という。あるいは略して単に正則ともいう。」
らしいです。つまり、一回微分可能であればよい。
また、1の証明は、収束円とCauchyの積分公式を用いてされるが
詳細は上記の本みれ。のってるから。上手く説明できましぇん。
おおざっぱなイメージ的には連続で滑らかな一本線になっている関数なら
OKそうって感じかな。
0029デフォルトの名無しさん
NGNG0030
NGNG（複素）関数論、解析などの授業を忘れちまったおいらは
本とかみながらなじゃないとついてけんわい。
0031デフォルトの名無しさん
NGNG0032デフォルトの名無しさん
NGNG0033デフォルトの名無しさん
NGNGsin(x)などの関数のある特定区間(sin(a)の近く,aは定数)で
近似曲線ってかけるよね。例えば 0付近では
一次: f(x) = x
三次: f3(x) = x^3 + 0*x^2 + (-0.917)*x + 0
五次: f5(x) = x^5 + 0*x^4 + (-0.975)*x^3 + 0*x^2 + 0.821*x + 0
みたいに。実際に書いてくらべてみよう。項が多くなればなるほど、
値が近い区間も広くなっていってだんだん sin(x)のグラフに
似てくるのがわかるはずだ。
0034デフォルトの名無しさん
NGNGしまった。上記のサンプルだと、
f(x) = xが一番綺麗にみえるかもしれない。
まぁ、いいか。いわんとすることはわかってもらえると思うし。
0035デフォルトの名無しさん
NGNGというよりも、宇治社中って、数学覚えたて学生の陶酔記事にしか見えんよ。
知識欲旺盛な半端学生に多い傾向なんだけどさ。
論理展開飛躍はないんだけど、肝心なとこの説明を省いてる。
そして、おそらく本人も、省いた意識すらないと思われ。
0036トの名無しさ
NGNG志井マガってそんなんばっかり。
0037デフォルトの名無しさん
NGNG数学板に行って「高校生でも分かるように説明して欲しい」
といえば、懇切丁寧に説明してくれるかもしれない・・・
1なんてそっちのけで、議論が進む可能性も高いが。
そういう私はちゃんと理解してないです。
f(x+h)の近似式 → テーラー展開 → マクローリン展開
という、カンの流れだけで。
数学科じゃなきゃ、こんなもんでしょ。
お役に立てずごめんなさい。
0038デフォルトの名無しさん
NGNG| | ∧＿∧ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
ﾉ__丶（〃´Д｀)＿ ＜ 数学の出来ないプログラマは逝ってよしですか？
||十||／ .| ¢、 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
＿ ||四||| | . ￣丶.)
＼ ||才||Ｌ二⊃ .￣￣＼ ｵﾒｰﾗｶﾞ ﾅﾆ ﾊﾅｼﾃﾙﾉｶ ｻｯﾊﾟﾘﾀﾞ ﾖ
||＼`~~´ プログレ♪ ＼
||＼||￣￣￣￣￣￣￣||￣
. ||￣￣￣￣￣￣￣||
0039デフォルトの名無しさん
NGNGいや、Ｃマガ読むと分かるんだけど、
「みんなが疑問に思っている証明すべき肝心な証明」が抜けてるんだよ。
いや、証明でなくてもいい。曖昧で直感的な説明でもいいから、ちゃっかり抜け
はしてほしくなかった。
誰も厳密な数学的証明を求めるなんて揚げ足とること要求してないの。
まぁ、Ｃマガのあの記事だけじゃないんだけどね。
大抵の数学教師も、あそこはすっ飛ばす。
「順次、微分していくと・・・」なんて感じで。
違うよ、疑問なのは、三角関数が多項式に近似できてしまう根拠を知りたいよー！
って感じで。教科書にも載ってないんだよね。
0040デフォルトの名無しさん
NGNGうーむ。あくまでも、「カン」で言えば…
平面上の2点を決めれば、それらを通る直線が1つ定まる。
3点を決めれば、それらを通る2次関数が1つ定まる。
4点を決めれば、それら全部を通る3次関数が1つ定まる…
というわけで、
なんでもいいから近似したい曲線の上に、n個の点を取れば、
そのn個の点を全部通る、n-1次関数が求まるはずなわけだ。
で、
点の間隔をどんどん細かくして、nを無限まで続ければ、
どんな関数（連続な曲線）でも、
多項式で近似できるということが、直感的に分かる。
そうすると、
あとは、どうやって各項の係数を求めるのか？と言う話になる…
これじゃダメ？
0041デフォルトの名無しさん
NGNG厳密な証明はCマガみたいなプログラミング雑誌に書くべきものでもないと思うんだけど。
大学初等レベルの教科書見ればちゃんと書いてあるわけだし。
図書館行きゃいくらでも調べれるし、2chで聞きたかったら数学板に行くのが一番いい。
0042デフォルトの名無しさん
NGNGおー素晴らしい説明。
いや、実は、僕の場合、就職後に知ってしまったので、
その直感的説明も聞いてしまったんだけどね。
でも本来、その手の説明をＣマガの時点でしてほしかったと。
（あるいは大学の講習の時点で）
>あとは、どうやって各項の係数を求めるのか？と言う話になる…
>これじゃダメ？
あとは、微分して各項の係数を求めていく機械的作業が残ってるのに
学校じゃ、ここばかりスポットライトを当ててると。
それを言いたかっただけ・・・
0043デフォルトの名無しさん
NGNG書いてるの? 厳密な証明じゃないよ。
0044デフォルトの名無しさん
NGNGこれってスプライン曲線とかの、
データを曲線に当てはめる方法と同じなのかな？
0045デフォルトの名無しさん
NGNG>この形で展開できるというのです。
これは大嘘。展開できるものについては展開できます。
0046デフォルトの名無しさん
NGNG当たり前じゃん。いたるところ微分可能なものなら、とか具体的なこと
書かないと何も言ってないのと同じ
0047デフォルトの名無しさん
NGNGワイヤシュトラウスの近似定理は結構すごいこと言ってるのね
0049デフォルトの名無しさん
NGNGその説明だとラグランジュ補間になってるような。
マクローリンだと１点なんじゃないの？
このスレみるまで忘れてたから間違ってるかも。
005040
NGNGそう。あれは補完なんだけど、質問のポイントが、
「連続関数を多項式で表すことができる理由」
だったから、まー、その説明としてはいいんじゃないかと。
1点から行くとなると、イメージ的に難しいんだよなー。
1点を決めて接線を引けば、近傍では近似できる
っていうのはいいだろうけど、その先が…
なぜ、f''/2なのか？
教授と同じで、「順次微分すると…」って説明になる？
005140
NGNGある連続関数f(x)を、多項式で近似したい。
（多項式で近似できるってことは分かってるとして）
x=0から始めることして、まずは、
直線（1次関数）で近似することにする。
すると、直感的にこれは、(0, f(0))でf(x)に接する
接線になるってことは、分かってくれますよね？
というわけで、xが0に近いときには、
f(x) ≒ f'(0)x + f(0) と近似できる。 … (1)
さらに、導関数f'(x)も近似したいということになると、
上でやったf(x)の例を、f'(x)に当てはめて、
f'(x) ≒ f''(0)x + f'(0) となり、 … (2)
(1)と(2)の結果を総合して、
f(x) ≒ (f''(0)/2) x^2 + f'(0) x + f(0) となる。
同様に、f''(x)も近似したいとすれば、
f''(x) ≒ f'''(0) x + f''(0) … (3)
(1)〜(3)より、
f(x) ≒ (f'''(0)/6) x^3 + (f''(0)/2) x^2 + f'(0) x + f(0)
で、これをずっと続けて、
f(x) ≒ f(0) + f'(0) x + … + (f^(n)(0)/n!) x^n + …
これだと、高次の導関数を順次等しくしていくと、
元の関数もそれに応じて、近似がよくなっていく
ということを説明し切れていないような…
0052デフォルトの名無しさん
NGNGってのが正しくないからでしょう
0053デフォルトの名無しさん
NGNG物語などの話の展開の仕方です。小説家がよく使います。
0054デフォルトの名無しさん
NGNGマクローリンならx=0におけるn階までの微分値が等しいn次多項式を
求めてるだけ。微分値が等しいんだから似たような関数になるだろうと。
nが増えていけばより近くなるってのは直感的に分かるし。
x=0から離れるほどずれは大きくなるっていうのも直感的に分かるよね。
0056デフォルトの名無しさん
NGNGこれがなぜってことなんだろうなぁ。
0057ﾏｸﾛｰﾘﾝの公式の形の意味
NGNGとゆう形に展開できたとすると、
f' = b + 2cx + 3dxx + 4exxx + ...
f'' = 2c + 3*2dx + 4*3exx + 5*4fxxx + ,,
f''' = 3*2*d + 4*3*2ex + 5*4*3fxx + ,,
みたいになるから、上の式に x=0 を代入すると
f(0) = a
f'(0) = 1*b
f''(0) = 2*1*c
f'''(0) = 3*2*1*d
みたいになる。だから
a = f(0)
b =f'(0)
c = f''(0) /2!
d = f'''(0) / 3!
となり、以下同様。これらを一番上の式に代入すると
ﾏｸﾛｰﾘﾝの公式になる。
こっちのほうが見やすいかなぁとおもって、
たとえば x^3 を xxx という風に書いてみた。
0058５７
NGNG無限級数の収束性については、ちゃんとした解析学の
教科書を参考にすればよいと思われ。
0059デフォルトの名無しさん
NGNGhttp://cheese.2ch.net/test/read.cgi/math/1010708150/l50
くだらねぇ問題はここへ書け ver.3.14159265
http://cheese.2ch.net/test/read.cgi/math/1010722815/l50
0060デフォルトの名無しさん
NGNG0061元、十四代酔っぱらい
NGNG突っ込みどころがいくつかあったりするYO!
>>40
>平面上の2点を決めれば、それらを通る直線が1つ定まる。
>3点を決めれば、それらを通る2次関数が1つ定まる。
>4点を決めれば、それら全部を通る3次関数が1つ定まる…
>というわけで、
細かいけど、定点のx座標は異なるのでしょう。
数学屋としてはともかく（？）、プログラマとしてはアウト。
>点の間隔をどんどん細かくして、nを無限まで続ければ、
>どんな関数（連続な曲線）でも、
>多項式で近似できるということが、直感的に分かる。
直観的にも論理的にも情緒的にもわかりません。
詳しく教えてください（ｗ
実解析的 = 連続、ですか？？？
>そうすると、
>あとは、どうやって各項の係数を求めるのか？と言う話になる…
そうすると、ならね。
0062デフォルトの名無しさん
NGNGあんたのように数学的に厳密な証明のほうがわかりやすいと
思ってるﾔｼにはこんな説明はそもそも必要ないの
0063元、十四代酔っぱらい
NGNGこんな場で証明をする気もスペースもないし。
っていうか、>>40の主張はウソ。
安易な言葉使いが問題になるのはプログラムも同じだろ？
それだけのこと、別に数学が特別なわけでもない。
数学コンプなプログラマー見苦しい。
0064デフォルトの名無しさん
NGNG…いや、なんとなく
0065デフォルトの名無しさん
NGNG>>24
fdlibm の e_rem_pio2.c みたいなのは、引数がゼロからかなり離れたとこでも
結果の精度を保証する為には、やむを得ないと思う(確か Intel 系の CPU は、
1584bits も持ってない筈)。
もちろん、そんなにゼロから離れた引数では使わない事が事前に分かっている
場合は、そこをパスするのもアリだと思う。
あと、PS2 は… 0〜2π ではなく、-π〜π で使えって事なんでしょう。
0066デフォルトの名無しさん
NGNG0067デフォルトの名無しさん
NGNG>>1が疑問を持ったように「ありとあらゆる関数」に関して
マクローリン展開ができるとは限らんわけだ。
展開可能であるためにはΣ(fn(0)/n!)が収束しなきゃいけない、と。
fn(0)はf(x)のn次導関数に0を代入したものね。
証明は・・・忘れたなぁ。
骨組み的には、収束してほしい式を上から押さえて、
押さえた式が目的の値に収束すればよい、
という感じだったように思うのだが。
記憶だけで書いてるので間違いあるかも。
0068デフォルトの名無しさん
NGNG級数展開して打ち切ったのと少し係数が違うよね？
あれはどういうふうにやってるの？
もしかしてひたすら差が小さくなるように繰返して求める？
006965
NGNG実際には(マクローリンではなく)チェビシェフ展開したやつをベースに、
使う区間等に合わせて係数を補正したのを使ってると思う。
■ このスレッドは過去ログ倉庫に格納されています