制約付き最小二乗問題に関しての質問です。
minimize |Vector(x)-Σ_[k=1,n]w_k*Vector(x_k)|^2 (変数はw_kです)
subject to Σ_[k=1,n]w_k=1
(Vector(x)の次元=Vector(x_k)の次元>n)
(w_k(1<=k<=n)はスカラーです)
上記の制約付き最小二乗問題をラグランジュで解きたいんですが、うまく解けません。
答えをclosed-formでかけるらしいのですが、導出過程がわかる方いらっしゃいましたらぜひご教授ください。
よろしくお願い致します。
No.1ベストアンサー
- 回答日時:
書きやすいように、次のように定めます。
y = Vector(x)
X = (Vector(x_1), ... ,Vector(x_n))
w = (w_1, ... , w_n)'
i = (1, ... , 1)'
j_k = (0, ... , 0, 1, 0, ... , 0)' 第k要素だけが1で後はゼロ
J = |1, -1, 0, ...| (n-1)×n の行列
|1, 0, -1, ...|
|.............|
|1....., 0, -1|
とするとこの問題は
min_w {y - Xw}'{y - Xw}
s.t. w'i = 1
と書けるはずです。したがって
e = y - Xw
とおけば、
L = e'e + 2R (1-w'i)
とラグランジュ関数が書けるので、最小化から
-X'e = R i
i'w = w'i = 1
という連立方程式が書けます。したがってn+1本の方程式とn+1個の未知数ですので、この式は単一の解を持つ、つまりclosed-formで書けることが分かります。
j_g' X'e = j_h' X'e
(j_g - j_h)' X'e = 0
これをたてにn-1本並べると、
J X'e = J X'X w - J X'y = 0
と書ける。したがって
J'J X'X w = J'J X'y
w = (J'J X'X)^{-1} J'J X'y
この J のあたり、三角行列やなんかを使ってもう少し簡単に出来るような気がするが、これでどうでしょう。
この回答への補足
丁寧な回答でとても参考になりました。
ありがとございました。
とある論文を見てみると、下記のように答えをかけるみたいです。
(at9_am様の表記をお借りしてます)
G=(y*i'-X)'(y*i'-X)
w=G^{-1}*i/(i'*G^{-1}*i)
がんばって解いてみましたが、この通りでません。
こちらの導出過程もわかる方がいらっしゃいましたら、ぜひご教授ください。どうぞよろしくお願い致します。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 【高1 数学Ⅰ 二次関数】 二次関数 f(x)=x^2-4ax+8a がある。ただし、aは正の定数と 3 2022/07/23 15:46
- 数学 数学 2時間数に関わる問題について教えてください。 x≧1 y≧-1 2x+y=5 であるとき、xy 7 2022/10/29 10:57
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- 数学 最大エントロピー原理をpythonで実装したい 2 2022/06/21 13:10
- 小学校 通分と約分について 通分と約分について教えてください。 子どもの宿題で、2/8と3/12を通分をし大 3 2022/11/28 21:07
- 数学 【 数I 2次関数 最大・最小 】 問題:関数y=x²+2x+c (-2≦x≦2)の最大値 が5であ 3 2022/06/19 08:41
- 統計学 t統計量とF統計量について 9 2023/01/05 14:23
- C言語・C++・C# このプログラミング誰か教えてください 9 2022/04/22 18:50
- 大学受験 ある大学の過去問なのですが、回答に解説がなく困っています。誰かこの問題の解説をつけて欲しいです(тт 1 2022/11/03 22:44
- 数学 数学の問題です。この問題が分かりません!解答によるとx=-1で最小値3となっています。私は最小値が- 3 2022/08/12 12:02
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
2進数のバイアス表現について
-
全員と同じグループを経験でき...
-
3次元での点群に対する最小二...
-
東大文系数学の問題です!
-
n!が10の40乗で割り切れるとき...
-
1/x+1/y≦1/2 , 2<x,2<yのとき、...
-
最大元と最小元をもつことの証...
-
中学受験用の小5算数の問題です
-
Gnuplot 最小二乗フィッティン...
-
2次関数の問題の場合分けで理解...
-
最小ハミング距離とは?
-
Xの二次関数 y=x ²ーmx+m(mは...
-
0は公約数?
-
最小分解体と拡大次数について
-
2次関数の問題です
-
mを実数とする。xの関係式f(x)=...
-
ある数字を割り切れる最小の数...
-
3変数の場合の最小値の求め方‐...
-
最大元と最大値 最小元と最小値...
-
大学入試の過去問です。 解き方...
マンスリーランキングこのカテゴリの人気マンスリー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つの放物線間の最短距離
-
間違いの理由を教えてください...
-
アルキメデスの公理について 任...
おすすめ情報