QR Factorization using Householder Transformation
156 次查看(过去 30 天)
显示 更早的评论
Hello!
I am trying to compute a function householder that taking a matrix A, returns the matrix Q R such that
This is my code so far but it seems that I have a mistake since ~= A:
Any help please?
function [Q,R] = householder(A)
% Write your code here.
[m,n] = size(A);
Q = A;
for j = 1:m
z = Q(j:end,j);
v = [ -sign(z(1))*norm(z)-z(1); -z(2:end)];
P = eye(size(z,1), size(z,1)) - (2/(v'*v))*(v*v');
Q(j:end, j:end) = P* Q(j:end, j:end);
end
R = Q' * A;
end
0 个评论
回答(1 个)
Monisha Nalluru
2021-3-10
Refer to the following blog post which gives example and clear understanding on Household Refelection and Household QR factorization with sample code.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Operators and Elementary Operations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!