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

实验原理 IIR数字滤波器 脉冲响应不变法变换原理 将模拟滤波器的s平面变换成数字滤波器的平面，而将模拟滤波器映射成数字滤波器。MATLAB信号处理工具箱中提供了IIR数字滤波器设计的函数，常用的函数： IIR滤波器的阶数选择 Buttord----巴特沃斯滤波器阶数选择 cheb1ord-----切比雪夫I型滤波器阶数选择 cheb2ord ----切比雪夫Ⅱ型滤波器阶数选择 IIR滤波器的设计 Buttrer-----巴特沃斯滤波器设计 cheby1-----切比雪夫I型滤波器设计 cheby2 ----切比雪夫Ⅱ型滤波器设计 maxflat----通用的巴特沃斯低通滤波器设计 巴特沃斯滤波器设计 巴特沃斯滤波器是通带，阻带都单调衰减的滤波器。 调用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’) 切比雪夫I型滤波器的设计 切比雪夫I型滤波器为通带波纹控制器；在通带呈现波纹特性，阻带单调衰减。 [N,Wc]= cheb1ord (Wp,Ws,Ap,As) [b,a]= cheby1 (N,Ap,Wc,options) 切比雪夫Ⅱ型滤波器的设计 切比雪夫Ⅱ型滤波器为阻带波纹控制器；在阻带呈现波纹特性，通带单调衰减。 [N,Wc]= cheb2ord (Wp,Ws,Ap,As) [b,a]= cheby2 (N,As,Wc,options) 已知模拟滤波器。可以利用脉冲响应不变法转换函数impinvar将其变为数字滤波器，调用格式为[bz,az]=impinvar(b,a,Fs). 双线性变换原理 采用非线性频率压缩方法，克服了脉冲响应不变法产生频率响应的混叠失真，使是s平面与z平面建立了一一对应的单值关系，消除多值变换性，频谱混叠现象。 已知模拟滤波器，可以利用双线性变换函数bilinear将其变换为数字滤波器，调用格式为 [bz,az]=bilinear(b,a,Fs). FIR窗函数法设计FIR数字滤波器 FIR滤波器设计需使频率响应H(ejw)逼近所要求的理想频率响应，窗函数法设计FIR数字滤波器是在时域中进行的，用窗函数截取无限长的hd(n),这样H(ejw)逼近与理想值。一旦选取了窗函数，其指标就是给定的，所以由窗函数设计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) 实验内容 一 要求通带截止频率fp=3KHz,通带最大衰减ap=1dB,阻带截止fs=4.5kHz,阻带最小衰减as=15dB,采样频率fc=30kHz,用脉冲响应不变法设计一个切比雪夫低通滤波器，并图示滤波器的振幅特性，检验wp,ws对应的衰减。 （1） wp=6*pi*10^3;ws=9*pi*10^3;ap=1;as=15; Fs=30*10^3; wp1=wp/Fs;ws1=ws/Fs; %参数设置 [N,WC]=cheb1ord(wp,ws,ap,as,s); [b,a]=cheby1(N,ap,WC,s); %采用切比雪夫I型滤波器设计 [bz,az]=impinvar(b,a,Fs); %采用脉冲响应不变法 w0=[wp1,ws1]; Hx=freqz(bz,az,w0); [H,W]=freqz(bz,az); dbHx=-20*log10(abs(Hx)/max(abs(H))); %显示幅频特性 plot(W,abs(H)); xlabel(相对频率);ylabel(幅频); grid 得到的结果如下 bz bz = -0.0000 0.0054 0.0181 0.0040 0 az az = 1.0000 -3.0591 3.8323 -2.2919 0.5495 dbHx dbHx = 1.0005