%% Signaalin suodatus
close all
% alustetaan aikaväli
t = linspace(-pi,pi,600);

%% luodaan signaali, ja häiriöitetty signaali
ft = .3*sin(3*t)+sin(t);
fh = ft+rand(size(ft))-.5;

%% piiretään signaalit
plot(t,ft,t,fh);

%% Siirrytään taajuustasoon

ffh  = fft(fh);
figure
plot(fftshift(abs(ffh)))

%% Tehdään suodatus
ffh(abs(ffh)<30) = 0;

%% Siirrytään takaisiin aikatasoon
fk = ifft(ffh);
figure
plot(t,ft,t,fk)
