Buenas,
He estado mirando y creo que esta bien.
Prueba esto:
% A mi si que se me parecen:
plot(abs(IR))
figure; plot(frec1,abs(Amp));
figure; freqz(
% este es el que menos se parece, probablemente por lo cercano que estan los picos en el espectro
% puedes probar utilizando una frecuencia más alta para
% que esten más separados.
% yo utilizaría la función filter
y= filter(b,1,x);
Espero que te sirva de algo.