非負条件の最小2乗法のプログラムを作成
したいのですが,参考文献やプログラムが
あれば,教えて下さい。考え方のヒントでも
よいです。
具体的には,関数f[{x}]=|[A]{x}-{b}|^2・・・(1)
がxi>=0(i=1,2,..,n)・・・(2)の条件で,
最小となる{x}を求めるという問題を解くプログラムを
作りたいと考えています。
ここで,{x}=(x1,x2,...,xn)
{b}=(b1,b2,...,bm)
[a]は サイズm x nのマトリクス(m>n)
であり,[A]および{b}は既知です。
条件(2)が無い場合の最小二乗のプログラムは
作成できますが,条件(2)を満足させるという条件がある
場合には,その条件を具体的にどのようにプログラム化
するのか不明であり質問した次第です。
宜しくお願いします。
No.1ベストアンサー
- 回答日時:
何も考えずに回答してみる。
普通に最小値を求めて領域 { x_i >= 0 | i = 1,2,...,n } に含まれなければ、この領域のエッジで最小値をとるしかないだろうから。
特に制限を設けずに {x} を求める
x が非負ならそれが答え
それ以外の場合
x_1 = 0 として {x} を求める、x_1 = 0 として {x} を求める... x_n = 0 として {x} を求める。
非負の解となったものの内で、f({x}) が最小であるものが答え
上記 n 個の解すべてが NG だった場合
x_1 = x_2 = 0 として {x} を求める、(以下 n(n-1)/2 通りの組合せ)
どんどん、= 0 の数を増やしていって、すべてに失敗した場合には (0, 0, ... 0) が最小値を与える。
No.2
- 回答日時:
これは線形の最小二乗法じゃあ使えないですよね。
非線形の場合は初期値によって停留点が複数ある場合があるから、
初期値をいろいろ変えて、xiが全部正になる解をを探すんですかね・・・・
なんとなくですが、xiが正にならなければならない問題で、
負の解が出てくるっていうのはモデルそのものを見直したほうがいいような気がするです。
ご解答ありがとうございます。
あれから,いろいろさがしまして,非負条件の最小二乗法の
コード自体は存在するようです。
しかし,ご指摘のとおり,余り素直な解き方ではないようなので,
私の問題の場合,最終的には2次計画法のコードを利用しました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# c言語 1 2022/11/09 13:05
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- C言語・C++・C# numpyスライス機能を使った数値計算 2 2023/05/08 16:01
- 統計学 主成分分析について 1 2022/06/14 14:53
- C言語・C++・C# C言語 3 2022/11/09 13:27
- その他(プログラミング・Web制作) pythonのプログラムについての質問です。 1 2023/05/26 10:31
- その他(プログラミング・Web制作) プログラムの勉強のおすすめは 7 2022/12/09 20:09
- 数学 写真の問題について質問なのですが、(別解)のやり方で、②(正の解1つと負の異なる解2つを持つ)の条件 2 2023/08/11 16:38
- C言語・C++・C# C言語 3 2022/10/04 15:07
- JavaScript 最小二乗法 2 2023/01/01 20:57
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
2進数のバイアス表現について
-
おしどり遊び(テイトの飛び石...
-
全員と同じグループを経験でき...
-
問題文は解答欄に載せます。 四...
-
3次元での点群に対する最小二...
-
数学?算数の問題です どのよう...
-
n!が10の40乗で割り切れるとき...
-
5406を13で割ったときの絶対値...
-
最小ハミング距離とは?
-
楕円のxの最大値と最小値
-
y=x^xの最小値
-
なぜ自然数を平方した数の約数...
-
次の問題を解いてください。 実...
-
高一の数学です。以下の2問の回...
-
領域の最大・最小の問題です!
-
なぜ、y"で上に凸、下に凸がわ...
-
高校数学、線分長の最小値
-
至急!!二次関数について aは...
-
分解能について
-
最小分解体と拡大次数について
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
全員と同じグループを経験でき...
-
2進数のバイアス表現について
-
3次元での点群に対する最小二...
-
数学2です x>0のとき、x + 16/(...
-
至急!!二次関数について aは...
-
高校数学1の問題集に、2次関数...
-
中学受験用の小5算数の問題です
-
5406を13で割ったときの絶対値...
-
x.>0ときγ(x)が最小値となるxの...
-
最大元と最大値 最小元と最小値...
-
おしどり遊び(テイトの飛び石...
-
1/x+1/y+2/z=1を満たす自然数解
-
3で割ると2余り、7で割ると4余...
-
1/x+1/y≦1/2 , 2<x,2<yのとき、...
-
次の問題を解いてください。 実...
-
整式について。
-
Xの二次関数 y=x ²ーmx+m(mは...
-
2つの放物線間の最短距離
-
間違いの理由を教えてください...
-
アルキメデスの公理について 任...
おすすめ情報