clear;      % Azzera il workspace di MATLAB
Fs=8000;
NBITS=16;
Ts=1/Fs;

Fc=400;     %Frequenza di taglio del filtro IIR da implementare come FIR
a=1-exp(-(2*pi*Fc/Fs));
b=1-a;
N=8000;
F=1600;


for n=0:N-1
    x(n+1)=sin(2*pi*F*n*Ts);
end

y(1)=0;
for n=2:N-1
    y(n)=a*x(n)+b*y(n-1);
end
plot(y(1:200));

