多次元行列の小さな値と大きい値の勾配降下

Question

私は勾配降下法を使って推定したい行列Rを持っています。コードは小さな次元行列（例えば100x1）に対して良好に働いていますが、次元を大きくするとNaN値が得られます（例：10x3）。助けてください！

R =[196,242,3;186,302,3;22,377,1;244,51,2;166,346,1;298,474,4;115,265,2;253,465,5;305,451,3;6,86,3]; 
N = length(R); 
M = size(R,2); 
K = 3; 

P = rand(N,K); 
Q = rand(M,K); 
alpha = 0.002; 
beta = 0.02; 


for iter = 1 : 5000 
e = R - P*Q'; 
P_new = P + 2*alpha*(e*Q-beta*P); 
Q_new= Q + 2*alpha*(e'*P-beta*Q); 
mse2(iter) = norm(R - P*Q')/norm(R); 

P=P_new; 
Q=Q_new; 
end 

R_est = P*Q'; 

Answer 1

問題はforループにあります。最大のダブルmatlabには1.7977e+308が格納できます。ループが7回実行された後、PとQの値は10^101のオーダーで、8回の実行後にはinfになります。私はあなたが推定行列によって何を意味するのかは分かりませんが、あなたのnanの値はPのQが急速に成長する方法から来ています。注記の他の1つは、mse2変数は使用されません。その目的は何ですか？