h=load('-ascii','solo_dati.fir');
m=length(h);
z=[0:0.001:1];
q=length(z);
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));
end
figure;
plot(z,respa)
title('risposta del filtro in ampiezza');
xlabel('frequenza normalizzata [0...nyquist]');
ylabel('dB');
