切比雪夫2型滤波器的设计2
2009-12-24 23:21阅读：
2、切比雪夫2型带通数字滤波器
clc
clear
%IIR-DF的设计
%确定数字滤波器的指标
Fs=10000;
fp=[440,560];fs=[275,900];
wp=2*fp/Fs;
%滤波器边界频率（关于pi归一化）
ws=2*fs/Fs;
Ap=0.1;
As=40;
%确定参数（
N：阶数，wso：阻带截止频率）
[N,wso]=cheb2ord(wp,ws,Ap,As);
%计算切比雪夫2型带通滤波器的系统函数的系数
[b,a]=cheby2(N,As,wso);
%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的符合信号st
st=mstg;
%调用滤波器实现函数filter进行滤波
y2t=filter(b,a,st);
%绘制DF的幅频响应曲线图
[H,w]=freqz(b,a,2000);
plot(w/pi,20*log10(abs(H)/max(abs(H))));
grid on
axis([0,1,-60,0]);
xlabel('w/\pi');
ylabel('幅度/dB');
title('切比雪夫2型带通数字滤波器幅频响应曲线');
%绘图显示切比雪夫2型低通数字滤波器输出波形y2(t)频谱Y2(t)
M=800;
T=1/Fs;t=0:T:(M-1)*T;
Tp=M*T;k=0:M-1;f=k/Tp;
figure;
subplot(2,1,1);
plot(t,y2t);
grid on
axis([0,0.08,-1,1]);
xlabel('t/s');
ylabel('y2(t)');
title('y2(t)的时域波形');
subplot(2,1,2);
Y2f=fft(y2t,M)/M;
stem(f,abs(Y2f),'.');
grid on
axis([0,Fs/5,0,0.25]);
xlabel('f/Hz');
ylabel('幅度');
title('y2(t)的频谱');
3、切比雪夫2型高通数字滤波器
clc
clear
%IIR-DF的设计
%确定数字滤波器的指标
Fs=10000;
fp=890;fs=550;
wp=2*fp/Fs;
%滤波器边界频率（关于pi归一化）
ws=2*fs/Fs;
Ap=0.1;
As=40;
%确定参数（N：阶数，wso：阻带截止频率）
[N,wso]=cheb2ord(wp,ws,Ap,As);
%计算切比雪夫2型高通滤波器的系统函数的系数
[b,a]=cheby2(N,As,wso,'high');
%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的符合信号st
st=mstg;
%调用滤波器实现函数filter进行滤波
y3t=filter(b,a,st);
%绘制DF的幅频响应曲线图
[H,w]=freqz(b,a,2000);
plot(w/pi,20*log10(abs(H)/max(abs(H))));
grid on
axis([0,1,-60,0]);
xlabel('w/\pi');
ylabel('幅度/dB');
title('切比雪夫2型高通数字滤波器幅频响应曲线');
%绘图显示高通滤波输出波形y3(t)
M=800;
T=1/Fs;t=0:T:(M-1)*T;
Tp=M*T;k=0:M-1;f=k/Tp;
figure;
subplot(2,1,1);
plot(t,y3t);
grid on
axis([0,0.08,-1,1]);
xlabel('t/s');
ylabel('y3(t)');
title('y3(t)的时域波形');
subplot(2,1,2);
Y3f=fft(y3t,M)/M;
stem(f,abs(Y3f),'.');
grid on
axis([0,Fs/5,0,0.25]);
xlabel('f/Hz');
ylabel('幅度');
title('y3(t)的频谱');
