发布时间 : 星期六 文章数字信号处理实验..更新完毕开始阅读
2. 利用巴特沃思模拟滤波器，通过双线性变换法设计数字带阻滤波器，数字滤
波器的技术指标为
采样周期为T=1。 程序代码
T=1; %设置采样周期为1 fs=1/T; %采样频率为周期倒数 wp=[0.25*pi,0.75*pi]; ws=[0.35*pi,0.65*pi]; Wp=(2/T)*tan(wp/2);
Ws=(2/T)*tan(ws/2); %设置归一化通带和阻带截止频率 Ap=20*log10(1/0.9);
As=20*log10(1/0.18); %设置通带最大和最小衰减 [N,Wc]=buttord(Wp,Ws,Ap,As,'s');%调用butter函数确定巴特沃斯滤波器阶数
[B,A]=butter(N,Wc, 'stop','s'); %调用butter函数设计巴特沃斯滤波器
W=linspace(0,2*pi,400*pi); %指定一段频率值
hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 subplot(2,1,1); plot(W/pi,abs(hf));
%绘出巴特沃斯模拟滤波器的幅频特性
曲线
grid on;
title('巴特沃斯模拟滤波器'); xlabel('Frequency/Hz'); ylabel('Magnitude');
[D,C]=bilinear(B,A,fs); Hz=freqz(D,C,W); subplot(2,1,2);
plot(W/pi,abs(Hz)); 频特性曲线
grid on;
title('巴特沃斯数字滤波器'); xlabel('Frequency/Hz'); ylabel('Magnitude');
进行试验 输出图像：
%调用双线性变换法 %返回频率响应 %绘出巴特沃斯数字带阻滤波器的幅
巴特沃斯模拟滤波器1Magnitude0.5000.20.40.60.811.21.4Frequency/Hz巴特沃斯数字滤波器1.61.821Magnitude0.5000.20.40.60.811.2Frequency/Hz1.41.61.82 实验分析
通过查workspace可以得出脉冲响应为：
滤波器幅度响应：滤波器相位响应：
3. 通过切比雪夫模拟滤波器，利用双线性变换法设计数字高通滤波器，数字滤
波器的技术指标为
采样周期为T=1，并观察T分别为0.01，0.1，0.25，0.5，0.8时数字滤波器的频率响应。 程序代码
T=1; fs=1/T; Wp=(2/T)*tan(0.75*pi/2);
Ws=(2/T)*tan(0.65*pi/2); %设置归一化通带和阻带截止频率
Ap=20*log10(1/0.9);
As=20*log10(1/0.18); %设置通带最大和最小衰减 [N,Wc]=cheb1ord(Wp,Ws,Ap,As,'s');%调用cheb1ord函数确定切比雪夫I型滤波器阶数
[B,A]=cheby1(N,Ap,Wc, 'high','s'); %调用cheby1函数设计切比雪夫I型滤波器
W=linspace(0,7*pi,400*pi); %指定一段频率值
hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 subplot(2,1,1);
%设置采样周期为1
%采样频率为周期倒数