problems in QR function
显示 更早的评论
function [Q,R] = cgs(A)
[m,n] = size(A);
Q = A; R = zeros(n);
for k = 1:n
R(1:k-1,k) = Q(:,1:k-1)'*A(:,k);
Q(:,k) = A(:,k)- Q(:,1:k-1)*R(1:k-1,k);
R(k,k) = norm(Q(:,k));
Q(:,k) = Q(:,k)/R(k,k);
end
end
please how to call this function 1000times?
采纳的回答
Walter Roberson
2019-4-3
A = .... whatever is appropriate
for iteration = 1 : 1000
[Q,R] = cgs(A);
end
This would not appear to make much sense to do unless you are changing A between iterations, or unless you are doing this in order to time the execution. If you are doing this for timing purposes then we would recommend that you instead use
timing = timeit(@() cgs(A), 2)
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
