====== Differences ====== This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
gibson:teaching:fall-2014:iam961:svddemo [2014/09/25 12:28] gibson created |
gibson:teaching:fall-2014:iam961:svddemo [2014/09/25 12:30] (current) gibson |
||
---|---|---|---|
Line 1: | Line 1: | ||
+ | ====== IAM 961 SVD demo ====== | ||
+ | |||
+ | A numerical demo of how uniqueness works for the SVD of real-valued matrices, conducted in Matlab. | ||
<code> | <code> | ||
% Form "random" matrix A with known SVD: A = U1 S1 V1' | % Form "random" matrix A with known SVD: A = U1 S1 V1' | ||
Line 6: | Line 9: | ||
% unitary Q and upper-triangular R. The QR decomp of a matrix with | % unitary Q and upper-triangular R. The QR decomp of a matrix with | ||
% normally distributed random elements with produce a unitary matrix | % normally distributed random elements with produce a unitary matrix | ||
- | % with its orthogonal columns in a random orientation. | + | % with its orthogonal columns in an isotropically random orientation. |
[U1,tmp] = qr(randn(5,5)); | [U1,tmp] = qr(randn(5,5)); |