A 回答 (3件)
- 最新から表示
- 回答順に表示
No.3
- 回答日時:
Excelソルバではどうでしょうか(最近こればっかり)。
非線形の最小二乗法は数値的に解くのがよろしいかと思います(その解くもののことをソルバと言います)。このようなソルバを使えば、最適化の常道として目的関数f=Σ[yi-(axi^2*log(b)+cxi+d)]^2を最小にするようなa,b,c,dを求めることができます。参考URL:http://office.microsoft.com/ja-jp/excel/HA011245 …
No.2
- 回答日時:
あなたが本気で
f(a,b,c,d) = y - (ax^2 * log(b) + cx +d)
に関して,Σf^2 を最小化しようとしているなら,
log(b)が存在するために微分して零云々の方法はそのまま使えません.
評価関数を微分して零にして云々は誤差項がパラメターの二次関数であるときだけ有効だからです.
二次関数以外の関数では微分して零の点は極大値か極小値であり,最小値である保証が何処にもないからです.
この場合は非線型最小二乗法と呼ばれ,一発で解を求める方法はありません.
この場合,Gauss-Newton法などが良く使われます.
これは適当な初期解から始めて反復的に解を改良してゆく方法です.
a*log(b)を一つの係数としてまとめて良いなら,
線型最小二乗法なので二次関数への当てはめで良いでしょう.
No.1
- 回答日時:
こんばんは。
1つの点(xk、yk)の、求める曲線からのずれをεkと置けば、
εk = axk^2*logb + cxk + d - yk
二乗誤差は、
εk^2 = (axk^2*logb + cxk + d - yk)^2
二乗誤差の合計をSとすると、
S =Σ[k=1→n]εk^2
= Σ[k=1→n](axk^2*logb + cxk + d - yk)^2
ここで、x1~xn、および、y1~yn はデータ、つまり、既知の数値なので定数です。
逆に、a、b、c、dの4つは、変数です!!!
二乗誤差の合計を最小にするということは、二乗誤差の合計が極小になる条件を求めるのと同じです。
極小とは言っても、変数が4つありますから、それぞれの変数について極小の条件を立てることになります。
そして、極小の条件とは、a~dについて、偏微分したものが、それぞれゼロになることです!!!
ですから、
∂S/∂a = 0
∂S/∂b = 0
∂S/∂c = 0
∂S/∂d = 0
を満たすa、b、c、d を求めれば、最小二乗法ができたことになります。
∂S/∂a = ∂/∂a・Σ[k=1→n](axk^2*logb + cxk + d - yk)^2
= Σ[k=1→n]xk^2*logb・2(axk^2*logb + cxk + d - yk)
= 0
∂S/∂b = Σ[k=1→n]axk^2*1/b・2(axk^2*logb + cxk + d - yk)
= 0
∂S/∂c = Σ[k=1→n]xk・2(axk^2*logb + cxk + d - yk)
= 0
∂S/∂d = Σ[k=1→n]2(axk^2*logb + cxk + d - yk)
= 0
という連立方程式になります。
(どっか、計算を間違えていたら、ごめんなさい。)
以上で考え方は終わりです。
あとは、表計算ソフトなどを利用して計算するとよいでしょう。
ところで、
alogb って、一体で1つの定数として扱えるので、
単なる二次関数への回帰でよいと思うんですが。
(いずれにしろ、最小二乗法の考え方は、上記のとおりです。)
この回答へのお礼
お礼日時:2008/10/19 16:29
わかりやすい解説ありがとうございます。
数学の苦手な私でもなんとか理解できました。
これを用いて最小二乗法の適用をしてみたいと思います。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 log底10真数1/75 ただし、 log底10真数2=0.3 log底10真数3=0.5とする 式 2 2022/05/30 22:51
- 数学 乗法公式の問題についてです。 (x-y)(2x+y)??? 2 2022/10/18 19:50
- Visual Basic(VBA) vbaで条件付き書式を設定したときの適用範囲について 1 2023/07/17 23:14
- 統計学 統計検定2級を取ろうと勉強中なのですが分からないことがあったので質問させていただきます。 スタージェ 6 2023/01/01 23:02
- 数学 logの含まれてる式からすなわち〜となっている式にどのように変形したらこうなるのか教えてください。 2 2023/04/18 22:16
- 法学 会社法 (監査役の退任) 第480条について 3 2022/05/20 21:09
- 国産車 車高調の選択について 3 2023/02/03 20:07
- その他(スマートフォン・携帯電話・VR) スマホの任意整理について。 auかんたん決済、スマホ料金が重なり、(機種の分割料金も含めて)支払い出 1 2023/01/10 18:54
- 数学 特定の座標点を通る回帰を行う方法について。 2 2022/10/10 10:27
- 数学 中3乗法公式の利用について。 ①103の二乗 ②99の二乗 ③101×99 ④52×48 この4問の 2 2022/04/20 22:14
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報