Implementación DSP del filtro digital FIR

[n,Wn,beta,ftype]=kaiserord([7 13 17 23],[0 1 0],[0.01 0.01 0.01],25);

%Encontrar el orden del filtro n= 37, beta=3.4

w1=2*fc1/fs; w2=2*fc2/fs; %Convertir especificaciones de filtro analógico a especificaciones de filtro digital

window=kaiser (n+ 1,beta);%Usar la función de ventana Kaiser

b=fir1(n,[w1 w2],window);Usar ventana aditiva para diseñar la respuesta de frecuencia estándar de la función fir1

freqz(b,1,512);% respuesta de frecuencia del filtro digital

t = (0:100)/Fs;

s = sin(2*pi*t* 5)+sin (2*pi*t*15)+sin(2*pi*t*30);% señal sinusoidal mixta

sf = filter(b,1,s);% filter Señales