Matlab を使用して QR 分解を実装し、行列 A = [1 2; でテストすることになっています。1 1; 2 3; 2 2]。これが今まで書いてみたコードです。しかし、いくつかのエラーが発生しており、どこに問題があるのか わかりません。誰でもそれを見つけることができますか？
A = [1 2; 1 1; 2 3; 2 2];
m = 4
Q=A;
Q(:,1) = A(:,1)/norm(A(:,1));
K = eye(m);
for j=2:n
Q(:,j) = ((K - Q(:,j-1)*Q(:,j-1)')*A(:,j))/norm((K - Q(:,j-1)*Q(:,j-1)')*A(:,j));
K = K - Q(:,j-1)*Q(:,j-1)';
end
R=Q'*A;