a7.m

数字滤波仿真程序

Fs=10;  fs=Fs/2;    %Frec. de muestreo
Fpb1=2;  fpb1=Fpb1/fs;  %Frec. de borde pasa banda

Fpb2=4; fpb2=Fpb2/fs; % Frec. de borde rechaza banda

Rp=1;
n=3  %n=2*N

wn=[fpb1,fpb2];
[b,a]=cheby1(n,Rp,wn);
[H,w]=freqz(b,a,512,1);


%------------------------------------------------------------

%Trazado de la respuesta en Magnitud

subplot(2,1,1);
plot(w,20*log10(abs(H)));
grid on;
title (['Filtro pasa-bajos, Respuesta en magnitud, orden=', num2str(n)]);
xlabel('frecuencia');
ylabel('H(f) db')
axis([0 0.5 -100 0])


subplot(2,1,2);
plot(w,angle(H));
grid on;
title (['Filtro pasa-bajos, Respuesta en magnitud, orden=', num2str(n)]);
xlabel('frecuencia')
ylabel('醤gulo de H rad')

%Resp al impulso

[y,t]= impz(b,a,60);
figure;
stem(t,y);
title (['Filtro pasa-bajos, Respuesta al impulso, orden=', num2str(n)]);

%Ploteo de los polos y ceros

z= roots(b);  %Zeros
p = roots(a);  % Polos

figure 
zplane(z,p)
title(['Trazade polos y ceros para Butter'])