• 締切済み

特異値分解法、最小2乗法について

「NUMERICAL RECIPES in C」 (ニューメリカル レシピーズ イン シー) と言う、本に載っている 特異値分解法と最小2乗法のプログラムを完成させたいのですが 誰か、教えていただけないでしょうか? 別に、それでなくても、特異値分解法と最小2乗法のC言語でのプログラム の作り方を教えてください! また、特異値分解法と最小2乗法についてもよくわからないので 何かわかりやすい説明あれば教えていただけると幸いです! よろしく、お願いします!

みんなの回答

回答No.1

管理者より: 同等の質問があるのでそちらをご参照下さい

参考URL:
http://www.okweb.ne.jp/kotaeru.php3?q=169753

関連するQ&A

  • C言語について教えてください。

    「NUMERICAL RECIPES in C」 (ニューメリカル レシピーズ イン シー) と言う、本に載っている 特異値分解法と最小2乗法のプログラムを完成させたいのですが 誰か、教えていただけないでしょうか? 別に、それでなくても、特異値分解法と最小2乗法のC言語でのプログラム の作り方を教えてください! よろしく、お願いします!

  • 最小二乗法について

    最小二乗法では二乗和の誤差 Σ[i=1~n]{Yi-(α+βXi )}^2 (iは添え字です) を最小化するα,βを推定することを考えますが、 これは単純にα,βで偏微分してそれを0とおいて 連立方程式を解くだけでよいのですか? といいますのも、2変数関数の極値を求める場合、 Hessianを計算して判別しますよね? ただ一階偏導関数が0になるからといって、 そこで極値をとるとは限らない気がしたので… それとも最小二乗法の場合は必ずとるようになっているのでしょうか? 手元の本には、 「この二乗和は非負値なので、αとβで偏微分したものを0とするα,βが上式を最小にする値である」 とあるのですが、一般に非負値だとこの ようなことが言えるのでしょうか?

  • 最小二乗法について

     y=ax+b+c/x という式での最小二乗法の求め方を 教えてください。

  • 連立方程式と特異値分解

    連立方程式と特異値分解 (方程式の本数)>(変数の数)という過剰条件の連立方程式における、最小二乗法を用いた最適解の算出について学習しています。 まず簡単に前提知識をまとめます。 (1)行列Aとベクトルb、xを用いた Ax=b という連立方程式について、||Ax-b||^2を最小にするために、 各変数での偏導関数=0を解くと、Aの転置行列A^Tを用いた正規方程式 (A^T)Ax=(A^T)b が得られ、擬似逆行列を用いて最適解xを求められること。 (2)Aの列ベクトルが1次従属である場合に、 (A^T)Aが逆行列を持たず、正規方程式が解けないこと。 (3)(2)の場合に特異値分解を用いると、数値計算精度の観点から 信頼性のある解を得られるらしい?こと。 ここまでは理解しています。 化学的な収支の問題を解く途中過程で特異値分解にぶつかってしまいまして、 線形代数を復習しているのですが、なかなか特異値分解の解説が見当たらず。。 素人的な観点で申し訳ないのですが、 【疑問1】特異値分解によって何故、(2)の状態の最適解を導けるのか? 【疑問2】特異値はどんな指標として扱うことができるのか? といった内容に関して、アドバイスを頂ければ幸いです。 特に【疑問2】については、1次従属の状態を判別するのに役立つとかなんとか、、 といった情報も見てしまい、その真偽を知りたいところです。 また、「連立方程式」「最小二乗法」「特異値分解」のキーワードで学習に向いている 参考書やWebサイト等あれば、是非お願いします。 長々と失礼致しました。

  • 2次の最小二乗法

    1次(ax+b)の最小二乗法は、 http://szksrv.isc.chubu.ac.jp/lms/lms1.html に載っている通りに求めるのですが、 2次(ax^2+bx+c)のa,b,cを求める式はどうなるのですか?

  • logの入った最小二乗法について

    Y=ln(x-b)-ln(x-c) といったように、xが二つに分かれてしまっていて、さらに式にあるようにbやcが入っているときの最小二乗法のやり方を教えてください。お願いします。

  • 円の最小二乗法のプログラム

    今、研究でVBで円の最小二乗法のプログラムを使いたいのですが全然わかりません。誰か詳しい人お願いします。

  • 多変数の最小二乗法

    Y=AX+Bという形で(Y,X)の組から係数A,Bを推定する方法が最小二乗法かと思います。(そのような整理が妥当であるか、意味があるか、についてはデータの分布形状や変動係数などで別途検証する必要がありますがここでは不問として) さて、その次ですが、YY=AX+BY+CZという関数で表示しようと考えたとします。(YY,X,Y,Z)というデータ系列があり、係数(A,B,C)を求めるというような解析です。(前述のようにそれが妥当かどうかは不問として。)係数A,B,Cを最小二乗法のようにできるだけ妥当に求める方法を教えていただきたいのですが。 1)たとえば、YY=AX+CCとして最小二乗法でAを求めたあと、残差CCについてCC=BY+CZとおいて再度最小二乗法でBを求めるというような手法を何度も行うとか? これだと解く順番に依存するようなので全部の順番でやって平均値を取るとか? 下手な考え休むに似たりのようですが。 2)カルマンフィルタ、ニューラルネットワークの手法で解析するということもあるでしょうか。その場合、YY=AX+BY+CZという枠組みではなくなりますが、それでもよいと言えばよいです。 このような問題を解く方法について解説されている本とか分野の名称(数理統計学とか?)などを教えて頂けると探索する手がかりとなります。今は手がかりすらはっきりしていないのです。 最小二乗法についてはその程度だったら専門分野ではなく、各分野で個別に解説されているようで込み入った問題になったときにどこを捜索したらよいかわからないもので質問しました。よろしくお願いします。

  • 最小2乗

    (1) (x,y)=(-3,4),(-2,1),(-1,0),(0,-1),(1,1),(2,2),(3,3)を 方程式 y = ax^2 + bx + c で最小2乗法に当てはめたときの a,b,cを求める (2) 要素(x,y,z)=(1,3,2),(2,1,1),(3,2,4),(4,5,3),(5,4,5)が z = ax + by + p とモデル化されるとき、残差の2乗Σp^2 が 最小になるように定式化して、最小となるa,bを求める という問題です。 これらを行列を用いて解きたいのですが、どのようにすればよいか分かりません。よろしくお願いします。

  • L1ノルムの最小化

    最適化問題で、 特異値分解法は最小2乗解を得ることができるので、L2ノルムの最小値に向かって推定が行われていると思うのですが、一般てきにポピュラーなL1ノルムを最小化するアルゴリズムにはどのような方法が使われているのでしょうか。 よろしく御願いします。