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

实 验 报 告 课程名称 实验项目名称 实验类型 实验学时 班级 学号 姓名 指导教师 实验室名称 实验时间 实验成绩 预习部分 实验过程 表现 实验报告 部分 总成绩 教师签字 日期 实验六 离散时间滤波器设计 预习报告 实验名称 离散时间滤波器设计 实验要求 设计和分析一组滤波器，获得对设计过程的深入认识，并且掌握几种标准设计方法的特点，要求合组讨论、单独实验，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，且滤波器具有线性相位。 实验六 离散时间滤波器设计