
function fs=fty(s);
 [n,m]=size(s);
 fs=fft([s(:,ceil(m/2)+1:m),s(:,1:ceil(m/2))].');
 fs=[fs(ceil(m/2)+1:m,:);fs(1:ceil(m/2),:)].';
