matlabで最小二乗法を解く

Question

前提・実現したいこと

matlabで最小二乗法にL1ノルムの正則化項を加えた式を解いています．
前回質問したことで，lasso関数を用いた部分のエラーは解決されました．
Xa:27×256　single
Xb:60×256　single
XaとXbは上記のようなサイズで，lasso関数を用いた解TLは
TL:1620×100　single
になりました．
最小二乗法の後に，
D＝sqrt(sum(sum((Xa-TL*Xb).^2))/256);
を解かなければならいのですが，TLのサイズが合わないため解くことができません

発生している問題・エラーメッセージ

配列の次元が一致していません

該当のソースコード

D＝sqrt(sum(sum((Xa-TL*Xb).^2))/256);

Answer 1

Y=rand(27,256);
X=rand(60,256);

n1=size(Y,1);
n2=size(X,1);
b=Y(:);
A=kron(X',eye(n1));
W=lasso(A,b);
[~,n3]=size(W);
w=reshape(W,n1,n2,n3);

D=sqrt(norm(A*W-b)^2/256)

Comments

[kakikukekonbu]

実行したのですが，エラーが返されました．
引数は 2 次元であるか、少なくとも 1 つの引数がスカラーでなければなりません。

サイズは以下のようでした．
b:6912×1
A:6912×1620
ｗ:1620×100

[WathMorks]

wではなくWです。

[kakikukekonbu]

訂正して実行したのですが同様のエラーが返されました

[WathMorks]

sqrt(sum(A*W-b).^2/256)
でWの列それぞれに対する誤差が計算できます

[kakikukekonbu]

解きたい式が
D＝sqrt(sum(sum((Xa-TL*Xb).^2))/256);
上記の式になるのですが，
回答を頂いた式とは同意でしょうか．

[kakikukekonbu]

W:1620×100
w：27×60×100

[WathMorks]

同じです。

[kakikukekonbu]

行列の次元が違うのでエラーが返されるのですが，解決法はありますでしょうか

[WathMorks]

sqrt(sum(A*W-b).^2/256)
でWの列それぞれに対する誤差が計算できます.

これでできませんか？

[kakikukekonbu]

やってみますね．
ちなみになんですけど，
w=reshape(W,n1,n2,n3)；
ここでreshapeしたwは後に使わないのですが，
reshapeする理由は何なのでしょうか．素人質問で申し訳ありません

[WathMorks]

前回のコードをコピペした結果です。
各正則化パラメータそれぞれに対するwが欲しいなら、それを使えばよい、という意味合いです。

[kakikukekonbu]

なるほどです．

実行したのですが，
同様のエラーが出ます

[WathMorks]

W=lasso(A,b);
が実行できるなら
A*W-b
も問題なく計算できると思うのですが。

[kakikukekonbu]

すみません．完全にこちらのミスでした，．
無事解決できました．ありがとうございます．