User:Esmehwiki

matlab:

%%%%%% fft of a rectangular pulse * LPF

% time domain

% 1 x 201

f = -10.01:0.1:9.99;               % freq vector

Wf = sin(pi*f)./(pi*f);

Hf=[zeros(1,95), ones(1,11), zeros(1,95)];

figure(1);

plot(f,Wf)

title('Wf')

xlabel('freq')

grid on;

figure(2);

plot(f,Hf)

title('Hf')

xlabel('freq')

grid on;

Yf=Wf.*Hf;

figure(3);

plot(f,Yf)

title('Yf')

xlabel('freq')

grid on;

wt=fft(Wf);

figure(4)

plot(abs(wt))

ht=fft(Hf);

figure(5)

plot(abs(ht))

yt=conv(wt,ht);

figure(6)

plot(abs(yt))

Yf2=fft(yt)

figure(7)

plot(abs(Yf2))