www.gusucode.com > Matlab动力系统和时间序列分析工具箱 > Matlab动力系统和时间序列分析工具箱/lab432/toolbox/SSA/ssa_get_elementary_vector.m

    function Y=ssa_get_elementary_vector(singular_value,U,V,delay);

% last modified 17.03.04

X=singular_value*U*V';
s=size(X);
N=s(2)+(s(1)-1)*delay;

Y=zeros(N,1);
D=Y;
for i=1:s(1)
    Y(1+delay*(i-1):delay*(i-1)+s(2))=Y(1+delay*(i-1):delay*(i-1)+s(2))+X(i,:)';
    D(1+delay*(i-1):delay*(i-1)+s(2))=D(1+delay*(i-1):delay*(i-1)+s(2))+1;
end
Y=Y./D;