哈尔滨工程大学数字信号处理实验六探析.doc

实 验 报 告 课程名称实验项目名称实验类型实验学时班级学号姓名指导教师实验室名称实验时间实验成绩预习部分实验过程 表现实验报告 部分总成绩教师签字日期 实验六 离散时间滤波器设计 预习报告 实验名称 离散时间滤波器设计 实验要求 设计和分析一组滤波器，获得对设计过程的深入认识，并且掌握几种标准设计方法的特点，要求合组讨论、单独实验，4学时。 实验原理 IIR数字滤波器的设计 a、脉冲响应不变法变换原理 将模拟滤波器的s平面变换成数字滤波器的平面，而将模拟滤波器映射成数字滤波器。MATLAB信号处理工具箱中提供了IIR数字滤波器设计的函数，常用的函数： IIR滤波器的阶数选择 Buttord----巴特沃斯滤波器阶数选择 cheb1ord-----切比雪夫I型滤波器阶数选择 cheb2ord ----切比雪夫Ⅱ型滤波器阶数选择 IIR滤波器的设计 Buttrer-----巴特沃斯滤波器设计 cheby1-----切比雪夫I型滤波器设计 cheby2 ----切比雪夫Ⅱ型滤波器设计 maxflat----通用的巴特沃斯低通滤波器设计 b、巴特沃斯滤波器设计 巴特沃斯滤波器是通带，阻带都单调衰减的滤波器。 调用buttord函数确定巴特沃斯滤波器的阶数，格式[N,Wc]=buttord(Wp,Ws,Ap,As) 调用butter函数设计巴特沃斯滤波器，格式[b,a]=butter(N,Wc,options) 利用以上两个函数可以设计出模拟滤波器，格式为 [N,Wc]=buttord(Wp,Ws,Ap,As,’s’) [b,a]=butter(N,Wc,options.’s’) c、切比雪夫I型滤波器的设计 切比雪夫I型滤波器为通带波纹控制器；在通带呈现波纹特性，阻带单调衰减。 [N,Wc]= cheb1ord (Wp,Ws,Ap,As) [b,a]= cheby1 (N,Ap,Wc,options) d、切比雪夫Ⅱ型滤波器的设计 切比雪夫Ⅱ型滤波器为阻带波纹控制器；在阻带呈现波纹特性，通带单调衰减。 [N,Wc]= cheb2ord (Wp,Ws,Ap,As) [b,a]= cheby2 (N,As,Wc,options) 已知模拟滤波器。可以利用脉冲响应不变法转换函数impinvar将其变为数字滤波器，调用格式为[bz,az]=impinvar(b,a,Fs). e、双线性变换原理 采用非线性频率压缩方法，克服了脉冲响应不变法产生频率响应的混叠失真，使是s平面与z平面建立了一一对应的单值关系，消除多值变换性，频谱混叠现象。 已知模拟滤波器，可以利用双线性变换函数bilinear将其变换为数字滤波器，调用格式为 [bz,az]=bilinear(b,a,Fs). 2、窗函数法设计FIR数字滤波器 FIR滤波器设计需使频率响应H(e^jw)逼近所要求的理想频率响应，窗函数法设计FIR数字滤波器是在时域中进行的，用窗函数截取无限长的hd(n),这样H(e^jw)逼近与理想值。一旦选取了窗函数，其指标就是给定的，所以由窗函数设计FIR滤波器就是有阻带衰减指标确定用什么窗，由过渡带宽估计窗函数的长达N。常用的有： hd=boxcar(N) ht=triang(N) hd=hanning(N) hd=hamming(N) hd=blackman(N) hd=kaiser(N,β) MATLAB中提供的fir可以用来设计FIR滤波器，调用格式为h=fir1(M,Wc,’ftype’,window) 实验内容 IIR数字滤波器的设计 要求通带截止频率fp=3KHz,通带最大衰减ap=1dB,阻带截止fs=4.5kHz,阻带最小衰减as=15dB,采样频率fc=30kHz,用脉冲响应不变法设计一个切比雪夫低通滤波器，并图示滤波器的振幅特性，检验wp,ws对应的衰减。 用双线性变换法设计一个切比雪夫I型数字高通滤波器，技术指标为：采样频率fc=2kHz,通带截止频率fp=700Hz,通带最大衰减ap=1dB,阻带边缘频率fs=500Hz,阻带最小衰减as=32dB。 窗函数法设计FIR数字滤波器 窗函数法设计低通滤波器 （1）N=26，分别利用矩形窗、汉宁窗和布莱克曼窗设计该滤波器，且滤波器具有线性相位。绘出脉冲响应h(n)及滤波器的频率响应。 （2）增加N，观察过渡带和最大尖峰值的变化。 2、 利用凯泽窗设计线性相位高通数字滤波器， N=31，且滤波器具有线性相位。 实验六 离散时间滤波器设计 实验报告 实验原理 IIR数字滤波器的设计 a