Evaluating svd() on slice of matrix array.

4 次查看（过去 30 天）
Morten Nissov
Morten Nissov 2020-1-20
编辑： 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?

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

回答（1 个）

David Goodmanson
David Goodmanson 2020-1-23
编辑：David Goodmanson 2020-1-23
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.

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by