ls=32;lT=16;lk=3;
[t,k]=meshgrid(1:ls,1:lk);s=-(k-1).*(k-3).*sin(2*pi*t/lT);
s=s';st=(s(:))';size([t,k]),size(s),size(st),
for i=1:((lk-1)*ls),S(i,:)=st(i:(i+ls-1));end;ST=S';size(ST),
F=zeros((lk-1)*ls,ls);
for i=1:((lk-1)*ls), F(i,:)=abs(fft(ST,ls));end;