实验六数字滤波器设计教程.docx

实验六数字滤波器设计 13通信工程 方鹏豪 134173708 1）基于Butterworth 模拟滤波器原型，使用双线性变换法分别设计IIR 数字低通滤波器，参数指标：通带截止频率Wp=0.2 π ，通带波动值Rp=1dB，阻带截止频率Ws=0.3 π ，阻带最小衰减Rs=20dB。滤波器采样频率Fs=2000 Hz。 要求：绘出幅频特性、相频特性图并比较结果，绘出滤波器零极点图，写出滤波器的系统函数。 代码： %数字滤波器指标 Wp=0.2*pi; Ws=0.3*pi; Rp=1; Rs=20; %转换为模拟滤波器指标 Fs=2000; T=1/Fs; Omgp=(2/T)*tan(Wp/2); %双线性变换法中频率预畸变 Omgs=(2/T)*tan(Ws/2); %模拟原型滤波器计算 [N, Omgc]=buttord(Omgp, Omgs, Rp, Rs, s) %计算阶数和截止频率 [z0, p0, k0]=buttap(N); %归一化原型设计 b0=k0*real(poly(z0)) %求原型滤波器分子分母系数 a0=real(poly(p0)); [b, a]=lp2lp(b0, a0, Omgc) %变换为模拟滤波器分子分母系数 %用双线性变换法计算数字滤波器系数 [num, den]=bilinear(b, a, Fs) %数字滤波器设计结果 %方法2：直接求模拟滤波器系数 [b, a]=butter(N, Omgc, s); [num, den]=bilinear(b, a, Fs); %方法3：直接求数字滤波器系数 [N, Wc]=buttord(Wp/pi, Ws/pi, Rp, Rs); [num, den]=butter(N, Wc); %求数字滤波器的频率特性 [H, w]=freqz(num, den); db_H=20*log10((abs(H)+eps)/max(abs(H))); %转换为分贝值 subplot(2, 2, 1); plot(w/pi, abs(H)); ylabel(幅度); xlabel(频率); title(幅度响应);axis([0, 1, 0, 1.1]); set(gca, XTickMode, manual, XTick, [0, 0.2, 0.3, 1]); ripple=10^(-Rp/20); rs=10^(-Rs/20); set(gca, YTickMode, manual, YTick, [0, rs, ripple, 1]); grid subplot(2, 2, 2); plot(w/pi, angle(H)/pi); ylabel(相位); xlabel(频率); title(相位响应);axis([0, 1, -1, 1]); set(gca, XTickMode, manual, XTick, [0, 0.2, 0.3, 1]); set(gca, YTickMode, manual, YTick, [-1, 0, 1]); grid subplot(2, 2, 3); plot(w/pi, db_H); ylabel(幅度(dB)); xlabel(频率); title(幅度响应(dB)); axis([0, 1, -50, 5]); set(gca, XTickMode, manual, XTick, [0, 0.2, 0.3, 1]); set(gca, YTickMode, manual, YTick, [-50, -20, -1, 0]); grid subplot(2, 2, 4); zplane(num, den); axis([-1.1, 1.1, -1.1, 1.1]); title(零极点图); 运算图形结果： 2）基于切比雪夫I 型模拟滤波器原型设计，使用脉冲响应不变法设计一个IIR 数字带通滤波器，参数指标：下通带截止频率Wp1=0.3 π ，上通带截止频率Wp2=0.7 π ，通带最大衰减Rp=1dB；下阻带截止频率Ws1=0.2 π ，上阻带截止频率Ws2=0.8 π ，阻带最小衰减Rs=20dB。滤波器采样频率Fs=2000 Hz。 要求：绘出幅频特性、相频特性图并比较结果，绘出滤波器零极点图，写出滤波器的系统函数。 代码： %数字滤波器指标 Wp1=0.3*pi; Wp2=0.7*pi; Ws1=0.2*pi; Ws2=0.8*pi; Rp=1; Rs=20; %转换为模拟滤波器指标 Fs=2000; T=1/Fs; Omgp1=Wp1*Fs; Omgp2=Wp2*Fs; Omgs1=Ws1*Fs; Omgs2=Ws2*Fs; Omgp=[Omgp1, Omgp2]; Omgs=[Omgs1, Omg