No.2ベストアンサー
- 回答日時:
Excelを用いるのであれば,a,b,c,dをパラメータとして
ソルバーを用いるのが良いのではないでしょうか.
まずは,ツール=>アドインでソルバーが入っているか確認して
入っていなければ,ソルバー アドインのところをチェック下さい.
シート上にパラメータa,b,c,d,をきまったセルに入れます.
xの値とyの値とax^d+bx+cの値を横に並べてy-(ax^d+bx+c)も並べます.
取得したx,yの分,同じ用に縦に並べます.
y-ax^d+bx+cのセルをsumsq()で平方和します.
この状態でツール=>ソルバーを起動させると,
目的セルは平方和を入れたセルで,目標値は最小値とします.
変化させるセルはパラメータa,b,c,dになります.
(各パラメータについて正のみなどの制約条件も入れられます.)
うまくいかない場合は,オプション設定で制限時間や反復回数を変化させると良いと思います.
この回答への補足
回答有難う御座います。
早速ソルバーで計算させてみました。(詳しい説明ありがとう御座いました。)
残差平方和がより小さくなるように、オプションで
制限時間や探索法など色々変えて見たのですが、
測定値となかなか一致しません。(全体的に右に少しシフトした感じに
なっています。)
現在、残差平方和が3.47719E-05ほどなのですが、
通常はどのくらいの残差平方和を一致の目安としているのでしょうか。
度々申し訳御座いませんが、回答のほど宜しくお願い致します。
No.4
- 回答日時:
もともとの値がいくらか分かりませんが,
1程度の数値だとすると,5桁落ちはたいしたものと思います.
最小二乗法などの回帰曲線は近似曲線なので,測定点を通ることは
ほとんどありえません.
全体的に右に少しシフトした感じ=>これはx軸上に平行移動
したようなイメージなのでしょうか.
そうであれば,y=ax^d+bx+cを
y=a(x-e)^d+b(x-e)+c=a(x-e)^d+bx+c'
と平行移動させれば良いかと思います.
パラメータを1個増やすことになりますが...
ちなみに,最小二乗法ではパラメータ数は
データ数の半分程度が限界と何かで読んだことがあります.
回答有難う御座います。
>最小二乗法などの回帰曲線は近似曲線なので,測定点を通ることは
ほとんどありえません.
そうだったんですね。勉強不足でした。有難う御座います。
>全体的に右に少しシフトした感じ=>これはx軸上に平行移動
したようなイメージなのでしょうか.
そうです。傾きも少し大きくなっていますが・・。
>そうであれば,y=ax^d+bx+cを
y=a(x-e)^d+b(x-e)+c=a(x-e)^d+bx+c'
と平行移動させれば良いかと思います.
パラメータを1個増やすことになりますが...
有難う御座います。そうやればいいんですね。
しかし、そうすると求めたい式の形が変わってしまう
気がするのですが・・・。
>ちなみに,最小二乗法ではパラメータ数は
データ数の半分程度が限界と何かで読んだことがあります.
それは初耳でした。大変勉強になりました。
フィッティングするように色々と試して見たいと思います。
有難う御座いました。
No.3
- 回答日時:
最も大きな次数の項の寄与が大きい(i.e. y≒ax^d)のなら、Logを取って
log y = d・log x + a
として最小二乗フィットさせれば、最大次数dを推定できると思います。
普通、実験結果からザっと次数を決定するときは、結果をLog-Logプロットして
その勾配を見る、ってしませんか?
最大次数が推定できたら、後は次数に応じて普通に最小二乗法を適用すれば
良いのでは。
この回答への補足
回答有難う御座います。
教えて頂いた通り
log y = d・log x + a
として最小二乗フィットさせ、
dを求めましたところ、d=0.00679となったので、この次数を用いて
y = ax^d + bx + c
に最小二乗フィットさせてみました。
その結果、測定値と大きくずれてしまったので、探索方法などを変えて
フィットさせて行くうちに、dはd=2に近づいて行きました。
また、未知パラメータの初期値は測定値のグラフにおいて、excelの
近似曲線の追加で出た値を入力しています。
間違えている点などありましたら、ご指摘頂けると有難いです。
No.1
- 回答日時:
dの値で大きくグラフが変化するでしょう。
このような場合には、行列ではちょっと無理です。ただし、
Y=aiX^i + a(i-1)X^(i-1) +... + a(1)x + a(0)
の形に直せれば、可能です。
なおせないとして
作図法で、エイヤーと、a,b,c,d を決定してください。
a,b,c,d のいずれかひとつを少し(普通は1/100くらい)変化させて、計算してみて、残差が減ったらば、変化後の値が正しい
という作業を飽きれるまで(山さが十分少なくなったと判断できるまで)つづけてくだ最。
今のパソコンですと早いですから、2-10秒でできるでしょう。
8086無wait , 無コプロで10分ぐらいの計算ですから。
早速の回答、有難うございます。
excelのソルバーで図を描かせて、そのような作業をしてみたいと思います。
有難う御座いました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
首吊りどこ締めるの
-
彼女のことが好きすぎて彼女の...
-
白血球が多いとどんな心配があ...
-
EXCELで条件付き書式で空白セル...
-
勃起する時って痛いんですか? ...
-
小数点以下を繰り上げたものを...
-
2つの数値のうち、数値が小さい...
-
イタリアから帰国する際、肉製...
-
エクセル指定した範囲からラン...
-
検便についてです。 便は取れた...
-
エクセルのラベルの値(文字列...
-
Excel 0目標に対して数字があ...
-
エクセルで数式の答えを数値と...
-
甲状腺が腫れているが血液検査...
-
風俗店へ行く前のご飯
-
excelのIF関数 A,Bの大きいほ...
-
テスターで断線を調べる方法教...
-
EXCELの条件付き書式で数式を空...
-
VLOOKUP関数を使用時、検索する...
-
ワードのページ番号をもっと下...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
首吊りどこ締めるの
-
彼女のことが好きすぎて彼女の...
-
白血球が多いとどんな心配があ...
-
勃起する時って痛いんですか? ...
-
精子が黄色?
-
検便についてです。 便は取れた...
-
これって喉仏ですか? 私は女性...
-
小数点以下を繰り上げたものを...
-
甲状腺が腫れているが血液検査...
-
精液の落とし方を教えてください
-
エクセル指定した範囲からラン...
-
EXCELで条件付き書式で空白セル...
-
2つの数値のうち、数値が小さい...
-
精子に血が・・・
-
エクセルのラベルの値(文字列...
-
イタリアから帰国する際、肉製...
-
ある範囲のセルから任意の値を...
-
リンク先のファイルを開かなく...
-
中出しをするとお腹が痛い・・・。
-
MIN関数で空白セルを無視したい...
おすすめ情報