csvd.m

来自「这是在网上下的一个东东」· M 代码 · 共 37 行

function [U,s,V] = csvd(A,tst) 
%CSVD Compact singular value decomposition. 
% 
% s = csvd(A) 
% [U,s,V] = csvd(A) 
% [U,s,V] = csvd(A,'full') 
% 
% Computes the compact form of the SVD of A: 
%    A = U*diag(s)*V', 
% where 
%    U  is  m-by-min(m,n) 
%    s  is  min(m,n)-by-1 
%    V  is  n-by-min(m,n). 
% 
% If a second argument is present, the full U and V are returned. 
 
% Per Christian Hansen, IMM, 06/22/93. 
 
if (nargin==1) 
  if (nargout > 1) 
    [m,n] = size(A); 
    if (m >= n) 
      [U,s,V] = svd(full(A),0); s = diag(s); 
    else 
      [V,s,U] = svd(full(A)',0); s = diag(s); 
    end 
  else 
    U = svd(full(A)); 
  end 
else 
  if (nargout > 1) 
    [U,s,V] = svd(full(A)); s = diag(s); 
  else 
    U = svd(full(A)); 
  end 
end