Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/grav/english/lecture/filtering/Hank.m Дата изменения: Wed Feb 18 22:26:17 2009 Дата индексирования: Mon Oct 1 23:38:45 2012 Кодировка:

%hankelization for singular spectral analysis
% program is written 17.02.2009 by L.V. Zotov

function X=hank(G1,L,N,L_min, K_max)
for ii=1:1:L_min
X1(ii)=0;
for j=0:1:ii-1
X1(ii)=X1(ii)+G1(ii-j,j+1);
end;
X1(ii)=X1(ii)/ii;
end;

for ii=(L_min+1):1:K_max
X1(ii)=0;
for j=0:1:L_min-1
if(L_min==L)
X1(ii)=X1(ii)+G1(L_min-j,(ii-L_min)+j+1);
else
X1(ii)=X1(ii)+G1(ii-j,j+1);
end;
end;
X1(ii)=X1(ii)/L_min;
end;

for ii=(K_max+1):1:N
X1(ii)=0;
for j=0:1:(L_min-(ii-K_max))
if(L_min==L)
X1(ii)=X1(ii)+G1(L_min-j,(ii-L_min)+j);
else
X1(ii)=X1(ii)+G1(K_max-j,j+1);
end;
end;
X1(ii)=X1(ii)/(L_min-(ii-K_max)+1);
end;

X=X1;