clear
close all

Nim = 50;
Nv = 480;
Nh = 640;

fs = Nv*Nh;
t = (1:Nim*Nv*Nh)'/fs;
%%

f0 = fs/(Nh*Nv);   % f0 = Tv;
f0 = fs/Nh* (1 - 1e-4) ;        % f0 = Th;
 
filt0 = exp(-(1:Nh*10)/fs*2e2);
filt0 = filt0/sum(abs(filt0));
f = sin(2*pi*f0*t);
f(f<0) = 0;

fout = zeros(Nv,Nh,Nim);
for n = 1 : Nim
    fout(:,:,n) = reshape(f( (n-1)*Nv*Nh+1:n*Nv*Nh)',Nh,Nv)';
end
%%
for j = 1 : Nim
imshow(uint8(255*squeeze(fout(:,:,j))));
end