N = 40;
t = linspace(0,2*pi,N);t = t(2:end);
x = sin(t+1)+cos(2*t+2);				% neporušený signál
x1 = sin(t+1)+cos(2*t+2)+rand(1,N-1)-0.5;	% signál s šumem
y = fft(x1);
y(5:end-3) = 0;  	% odstranìní vysokých frekvencí - šumu
x2 = ifft(y);
plot(t,x,'k',t,x1,'k.',t,x2,'k--');












