Evaluating svd() on slice of matrix array.

Question

Morten Nissov 2020-1-20

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/501077-evaluating-svd-on-slice-of-matrix-array

编辑： David Goodmanson 2020-1-23

I have an array of matrices such that

size(A) == [3, 3, 1e3]

for example. These 1000 [3x3] matrices must be orthonormal so I am attempting to project each one to the nearest orthonormal basis, using svd() and the approximation

This function is made to work on a single [3x3] matrix at a time however. A workaround could be using for loops like

[~,~,np] = size(A);
for i=1:np
    [U,~,V] = svd(A(:,:,i));
    A(:,:,i) = U*V';
end

but this function will be called very often with high numbers of matrices so I am attempting to make efficient. Is there a better way to do this?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Goodmanson 2020-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/501077-evaluating-svd-on-slice-of-matrix-array#answer_411569

编辑：David Goodmanson 2020-1-23

在 MATLAB Online 中打开

Hi Morten

I take it that your matrices are close to being orthonomal already. Try

[Q, ~] = qr(A(:,:,i));
A(:,:,i) = Q;

which is about three times faster, not counting overhead to read and write to A(:,:,i). The result is slightly different, but of course still orthogonal.