kron_svd.m

curves ti si s a nice code.

function [U, s, V] = kron_svd(A)
%
%   [U, s, V] = kron_svd();
%
%  Given a kronMatrix object, this function computes
%  an "SVD" factorization that is to be used for a
%  psfMatrix object.
%
%  On Entry:
%     A  -  kronMatrix obect
%
%  On Exit:
%     U, V  -  kronMatrix objects, orthogonal
%        s  -  column vector containing the singular values
%              of A.  Note that these are not sorted from largest
%              smallest.
%
[U, S, V] = svd(A);
s = kron(diag(S.a), diag(S.b));