function [m,ff]=risposta(nome)
% leggi il file con i dati del filtro
h=load('-ascii',nome);
m=length(h);
% plotta la risposta del filtro
z=[0:0.001:1];
q=length(z);
% controlla quando scende sotto di 3 db
ff=0;
level=-3;
%
for j=1:q
    d=0;
    for k=1:m
        d=d+h(k)*exp(-i*pi*z(j)*(k-1));
    end
    respa(j)=20*log10(abs(d));
    if( respa(j)> level) 
	ff=z(j);
	end
end
figure;
plot(z,respa)
title('risposta del filtro in ampiezza');
xlabel('frequenza normalizzata [0...nyquist]');
ylabel('dB');
% campionamento a 44100, Nyquist a 22050
ff=ff*22050;
%
X=sprintf('\nfiltro fir di %d termini, -3dB a %d Hertz\n',m,ff);
disp(X);
end
