Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/grav/russian/lecture/filtr/labs/hank.m Дата изменения: Sat Nov 27 16:26:25 2010 Дата индексирования: Mon Oct 1 23:41:11 2012 Кодировка:

function X=hank(G1,L,N,L_min, K_max)
if L~=L_min X=0; return;
end;

for ii=1:1:L_min
X1(ii)=0;
X1(N-ii+1)=0;
for j=0:1:ii-1
X1(ii)=X1(ii)+G1(ii-j,j+1);
X1(N-ii+1)=X1(N-ii+1)+G1(L_min-(ii-j)+1,K_max-j);
end;
X1(ii)=X1(ii)/ii;
X1(N-ii+1)=X1(N-ii+1)/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)-1)
% 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;