以前も似たような内容の質問をしたのですが、
またわらからない問題が出てきたので投稿しました。
微少量δxを含む次の式を近似したいと思ってます。
f(δx)=δx*exp(-δx)(1 + 2*δx)
ここで解析のために、できれば右辺のカッコ内のδxについては残したままで、
その手前のexp(-δx)のみを1で近似してやりたいです
(展開してでてくるδx^2はのこしてexp(-δx)を消去したい)。
このような近似は行っても大丈夫でしょうか。δxはだいたい10^-2のオーダーです。
自分で考えろと突っ込まれそうですが、
可能かどうか、理由とあわせて教えてもらえると助かります。
No.6ベストアンサー
- 回答日時:
#4です。
A#5の補足について
既に適当な(不適当な)近似をした式に対して、
既に大きな誤差が入っている式を更に近似しても
その近似は、誤差の入った近似式を更に近似するだけで
大元の関数の近似式になるわけではありません。
そういった大元の関数の式を伏せて、質問される事自体
質問として投稿されるのは適当ではありません。
質問者に迷惑をかけるだけです。回答の自体の意味がなくなってしまうからです。
>次の式を一度テイラー展開して出てきたものです。
>g(x)=(1-exp(-x))exp(-x) + 2*(1-exp(-x))^2*exp(-x)
>+ 3*(1-exp(-x))^3*exp(-x)…(◆)
最初から、この式を示して、この式の近似式はどうなるか、といった質問にしていただきたかった。今後はそうして欲しいということです。
近似の微小な変数xをδxと表現することもしません。
0<x<10^(-3)などと適用するxの範囲を示せば十分です。
(◆)の式のg(x)を直接マクローリン展開して求めるのが
正しい近似式を求める方法です。
(◆)のg(x)のマクローリン展開は
g(x)=x+(1/2)x^2+(1/6)x^3-(95/24)*x^4)+ …
となるのでこの式の
第1項から2項目まで、または、3項目までの和の式
g(x)≒x+(1/2)x^2
または
g(x)≒x+(1/2)x^2+(1/6)x^3
が、本来の適切な近似式といえます。
この回答への補足
改めて見返してみると、少し勝手な物言いでした。
上記のコメントは訂正させてください。
解答していただく以上、きちんとした配慮は必要かもしれません。
ともかくありがとうございました。
わかりやすい解説ありがとうございます。
投稿者の数学知識が未熟であるために回答者の
皆さんを混乱させてしまって申し訳ない限りです。
私はテイラー展開による近似手法についてまだ
よく理解できてないようです。
ただ質問自体は
<f(δx)=δx*exp(-δx)(1 + 2*δx)
<ここで解析のために、できれば右辺のカッコ内のδxについては残したままで、
<その手前のexp(-δx)のみを1で近似してやりたいです
これがすべてですので、この式自体がどのように導かれ、
それが正しいかについては回答を望んでいませんし、
式に含まれる近似が不適切であったからといって、それが回答者さんにとって迷惑とまでいわれますと少々気分が悪いです。そのような不利益(未熟な質問にいらいらする)が生じることも理解したうえで個人の趣味で回答者さんは解答されているからと考えるからです。
ただ私自身に関してはもっと数学を勉強する必要があると再確認されました。
No.5
- 回答日時:
数値については既に回答が出ています。
ちょっと別の面からです。
1番前についているδxが混乱のもとですね。
このδx自体が何かの近似で出てきたものだろうと思います。
そうであればf(δx)をδxの1次で近似するということしか意味を持たない事になります。
f(δx)~δx
です。他の関数と比較して関数fの寄与がどの程度かというときにはこれで判断します。
f(δx)の第二項目を問題にするのであれば別の関数g(δx)と一項目だけでは違いが出てこない場合です。その場合、もし最初のδxが何かの近似で出てきているのであればその2項目も考慮しないと正しい評価は出来ません。
exp(-δx)についても2項目を考慮する必要があります。
そのまま使うか1-δxとして計算するかです。
(1+2δx)という式も普通は何かの近似式です。3項目(δxの2乗の項)が小さいとして既に省略されている式です。したがってexp(-δx)を展開して3項目まで計算してもあまり意味はありません。
どこらあたりの桁に影響が現れて来るかという確認のためには意味があるでしょうが数値自体には意味がないのです。
この回答への補足
お礼とあわせてお伺いしたいのですが、
質問の式は、No.5さんがおっしゃられるように、
次の式を一度テイラー展開して出てきたものです。
g(x)=(1-exp(-x))exp(-x) + 2*(1-exp(-x))^2*exp(-x)
+ 3*(1-exp(-x))^3*exp(-x)
これから得られた結果に再度近似を適応するなんて意味がないんじゃないかと思いながら質問してしました。
この式を2次まで評価したうえでexp(-x)を簡単にしてやりたかったのですが
実際にマクローリン展開して2次までとると近時が得られるのだと思うのですが、伺いたいのは近似でexp(-x)を1とおいておきながら
x^2を評価することが意味を持つのかということです。
皆さんご解答ありがとうございます。
返事が送れて申し訳ありません。
代表してNo.5のかたにレスさせてもらいました。
No.1、2、3さんもありがとうございます。
No.4
- 回答日時:
#2です。
A#2のより精度の高い3項目まで採用した式
>f(δx)≒δx+(δx)^2-1.5*(δx^3)
で最後の項の「( )」の「 )」の位置がずれていましたので
以下と差し替えて下さい。
f(δx)≒δx+(δx)^2-1.5*(δx)^3
おまけにこの3項目まで採用したときの計算精度は
x=10^(-2)と置くと
f(10^(-2))≒0.010098500
真の値f(10^(-2))=0.010098508304242
誤差は0.00000001未満ですね。
なお、#1さんのA#1の近似式は、#3さんも指摘されているように間違いですね。
多分、exp(-x)とすべき所をexp(x)と早とちりされたのでしょう。
数値を入れて、近似式と真の値を比較すれば誤差が大きく出ますので分かるかと思います。
No.3
- 回答日時:
exp のテーラー展開は
exp(-ε) = 1 - ε + (1/2)ε^2 - (1/6)ε^3 + …
ですから、
f(ε) = ε・(1 + 2ε)・exp(-ε) = ε・(1 + 2ε)・(1 - ε + εの二次以上の式)
という変形から、
一次近似なら f(ε) = ε + (εの二次以上の式)
二次近似なら f(ε) = ε + ε^2 + (εの三次以上の式)
となります。
f(ε) ≒ ε または
f(ε) ≒ ε + ε^2 であって、
f(ε) ≒ ε + 2 ε^2 はナイでしょう。
No.2
- 回答日時:
f(x)=x*exp(-x)*(1+2x)
をマクローリン展開してやると
f(x)=x+x^2-(3*x^3)/2+(5*x^4)/6 - …
だから,2次の項までで近似すれば
f(x)≒x+x^2
したがって
f(δx)≒δx+(δx)^2
となります。
δx=10^-2とおくと
2次の項までの近似値=0.0101000
真の値=0.010098508304242
となり、真の値の小数点以下6桁目を四捨五入してやれば
0.010100
と小数点以下6桁までは一致しますね。
この程度の誤差でよければ「(δx)^2」の項まででいいですね。
それ以上の小さい誤差が要求されるなら「(δx)^3」の項までとって
f(δx)≒δx+(δx)^2-1.5*(δx^3)
とすればいいでしょう。
No.1
- 回答日時:
テイラー展開してみてはどうでしょう?
exp(x)=1+x+(1/2)x^2+…+(1/n!)x^n+… でしたが今回は
exp(x)=1+(xの1次以上の式) とでもしておけば十分
すると
f(δx)=δx(1 + 2*δx)(1+(δxの1次以上の式))
=δx+2δx^2+(δxの3次以上の式)
δxが十分小さければ(3次以上の式)の部分は2次以下の式と比べて無視できるはずなので
f(δx)≒δx+2δx^2
と思っていいでしょう。
もうちょい精度を上げたければ
exp(x)=1+x+(1/2)x^2+(xの3次以上の式)
などとして次数を上げてみてもいいと思います。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・「I love you」 をかっこよく翻訳してみてください
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・昔のあなたへのアドバイス
- ・かっこよく答えてください!!
- ・あなたが好きな本屋さんを教えてください
- ・スタッフと宿泊客が全員斜め上を行くホテルのレビュー
- ・【大喜利】【投稿~8/27】 こんなガソリンスタンド二度と来るか!なぜそう思った?
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・【お題】動物のキャッチフレーズ
- ・【お題】甲子園での思い出の残し方
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・「それ、メッセージ花火でわざわざ伝えること?」
- ・自分用のお土産
- ・人生で一番お金がなかったとき
- ・一番好きなみそ汁の具材は?
- ・泣きながら食べたご飯の思い出
- ・ちょっと先の未来クイズ第1問
- ・ゴリラ向け動画サイト「ウホウホ動画」にありがちなこと
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
3:4:5の直角三角形の角度
-
テイラー展開とマクローリン展...
-
テイラー展開 1/(1-x^3), 1/(1-...
-
x=√2+√3+√5+√7の整数部分aは?
-
どの参考書にも一次近似式で、...
-
3x3行列の固有値の求め方
-
固有値の値について
-
行列の2分の1乗の計算の仕方
-
ベクトルの内積の表記の仕方、...
-
数学について質問です。 nを正...
-
パソコンで行列はどう書けばいい?
-
行列の積の固有値
-
三角関数の「1/3倍角の公式...
-
行列と級数の和
-
日本数学オリンピック2000年予...
-
【数学】「sinθ=√1-cosθ^2」が...
-
6分の11πはこの範囲(−3分のπか...
-
円環の体積 断面積が半円の内側...
-
arctan(y/x)でxとyを極座標で...
-
行列の積、どうしてそのように...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
テイラー展開とべき級数展開の...
-
テイラー展開 1/(1-x^3), 1/(1-...
-
1/x^2のような近似曲線の求め方
-
テイラー展開とマクローリン展...
-
どの参考書にも一次近似式で、...
-
テイラー展開とローラン展開は...
-
お茶碗の展開図ってどう書くの?
-
n次近似式とテイラー展開について
-
x=√2+√3+√5+√7の整数部分aは?
-
1/sinh xの近似
-
よい例が思い付きませんが、 極...
-
三角関数 角度を求める方法
-
微少量を含む式の近似について
-
テイラー展開の幾何学的意味に...
-
この関数の最大最小、またその...
-
逆三角関数
-
ピカールの逐次近似について。
-
1の近傍での増減の近似
-
微少量の近似について
-
e^0.1 の評価
おすすめ情報