Fastest way to compute spectral norm of a matrix?
35 次查看(过去 30 天)
显示 更早的评论
Hello
I can't find any mention of the spectral norm in the documentation. I want to calculate
max(|Ax|)/x for any vector x, given a matrix A.
There's a method via the eigenvalues of A'A but that sounds very memory intensive, I was wondering if anyone knows a faster in-built option.
Thanks
Mike
0 个评论
采纳的回答
Matt J
2013-6-11
See the NORM command. It accepts matrix inputs.
3 个评论
Matt J
2013-6-11
编辑:Matt J
2013-6-11
Hmmm. Looks like NORM is not as fast as the direct calculation,
A=rand(1000);
tic;
n1=sqrt( max( eig(A.'*A) ) );
toc; %Elapsed time is 0.150860 seconds.
tic;
n2=norm(A);
toc; %Elapsed time is 0.319899 seconds.
However, I'm pretty sure it's because NORM uses SVD instead of EIG, which is probably more numerically stable.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!