发布时间 : 星期三 文章数字信号处理实验..更新完毕开始阅读0244887a82c4bb4cf7ec4afe04a1b0717ed5b365
程序计算结果如下：
X8e?k??1*?N?k?及x4?n??x8r?n??IDFT[X8e?k?]，正好与图X8?k??X82中x4?n?的16点X4?k?及x4?n?相同。
jX8o?k??1*?N?k?及x5?n??x8i?n??IDFT[X8o?k?]/j，正好与X8?k??X82图中16点的X5?k?及x5?n?相同。
实验三 双线性变换法设计IIR数字滤波器
一、 实验目的
1. 掌握利用双线性变换法设计IIR数字滤波器的原理和具体方法。 2. 深入理解利用双线性变换法设计IIR数字滤波器的优缺点及使用范围。 二、 实验内容
1. 利用巴特沃思模拟滤波器,通过双线性变换法设计巴特沃思数字滤波器,数字
滤波器的技术指标为
采样周期为T=1。 程序代码
T=2; %设置采样周期为2 fs=1/T; %采样频率为周期倒数 Wp=0.25*pi/T;
Ws=0.35*pi/T; 率
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,'s'); %调用butter函数设计巴特沃斯滤波
器
W=linspace(0,pi,400*pi); hf=freqs(B,A,W); subplot(2,1,1);
plot(W/pi,abs(hf)/abs(hf(1))); 频特性曲线
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)/abs(Hz(1))); 的幅频特性曲线
grid on;
title('巴特沃斯数字滤波器'); xlabel('Frequency/Hz');
%指定一段频率值
%计算模拟滤波器的幅频响应 %绘出巴特沃斯模拟滤波器的幅
%调用双线性变换法 %返回频率响应 %绘出巴特沃斯数字低通滤波器
ylabel('Magnitude');
进行试验
巴特沃斯模拟滤波器1Magnitude0.5000.10.20.30.40.50.60.7Frequency/Hz巴特沃斯数字滤波器0.80.911.5Magnitude10.5000.10.20.30.40.50.6Frequency/Hz0.70.80.91输出图像： 实验分析
通过查workspace可以得出脉冲响应为：
滤波器幅度响应：滤波器相位响应：