信号处理实验六报告.docxVIP

实验六 离散时间滤波器设计实验原理IIR数字滤波器设计（一）、脉冲响应不变法变换原理脉冲响应不变法将模拟滤波器的平面变换成数字滤波器的平面，从而将模拟滤波器映射成数字滤波器。IIR滤波器的系统函数为(或)的有理分式，即一般满足。1、转换思路：若模拟滤波器的系统函数只有单阶极点，且假定分母的阶次大于分子的阶次，表达式：2、s平面与z平面之间的映射关系。IIR数字滤波器设计的重要环节是模拟低通滤波器的设计，典型的模拟低通滤波器有巴特沃思和切比雪夫(I型和II型)等滤波器。由模拟低通滤波器经过相应的复频率转换为，由经过脉冲响应不变法就得到所需要的IIR数字滤波器。（二）、巴特沃思滤波器设计巴特沃思滤波器是通带、阻带都单调衰减的滤波器。(1) 调用buttord函数确定巴特沃思滤波器的阶数，格式[N,Wc]=buttord(Wp,Ws,Ap,As)其中：Wp，Ws为归一化通带和阻带截止频率；Ap，As为通带最大和最小衰减，单位为dB； N为滤波器阶数，Wc为3dB截止频率，对于带通和带阻滤波器，Wc=[W1，W2]为矩阵，W1和W2分别为通带的上下截止频率。(2) 调用butter函数设计巴特沃思滤波器，格式[b,a]=butter(N,Wc,options)其中：options=’low’, ‘high’, ‘bandpass’, ‘stop’，默认情况下，为低通和带通。 b和a为设计出的IIR数字滤波器的分子多项式和分母多项式的系数。注意，利用以上两个函数也可以设计出模拟滤波器，格式为[N,Wc]=buttord(Wp,Ws,Ap,As,’s’)[b,a]=butter(N,Wc,options,’s’)其中：Wp、Ws和Wc均为模拟频率。（三）、切比雪夫I型滤波器的设计切比雪夫I型滤波器为通带波纹控制器：在通带呈现纹波特性，在阻带单调衰减。[N,Wc]=cheb1ord(Wp,Ws,Ap,As)[b,a]=cheby1(N,Ap,Wc,options)其中的参数含义和巴特沃思的相同。（四）、切比雪夫II型滤波器的设计切比雪夫II型滤波器为阻带波纹控制器：在阻带呈现纹波特性，在通带单调衰减。[N,Wc]=cheb2ord(Wp,Ws,Ap,As)[b,a]=cheby2(N,As,Wc,options)其中的参数含义和巴特沃思的相同。已知模拟滤波器，可以利用脉冲响应不变法转换函数impinvar将其变换为数字滤波器，调用格式为[bz,az]=impinvar(b,a,Fs)其中b,a分别为模拟滤波器系统函数分子、分母多项式系数；Fs为采样频率；bz、az为数字滤波器系统函数的分子、分母多项式系数。（五）、双线性变换法变换原理为克服脉冲响应不变法产生频率响应的混叠失真，可以采用非线性频率压缩方法，使平面与平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，这就是双线性变换法。1、转换思路：写出微分方程差分方程写出由于双线性变换法中，s到z之间的变换是简单的代数关系，得到数字滤波器的系统函数和频率响应，即设模拟系统函数的表达式为应用双线性变换得到的表达式2、s平面与z平面之间的映射关系。用不同的方法选择c可使模拟滤波器频率特性与数字滤波器频率特性在不同频率处有对应的关系。(1) 采用使模拟滤波器与数字滤波器在低频处有较确切的对应关系，即在低频处有。当较小时，。(2) 采用数字滤波器的某一特定频率(例如截止频率)与模拟原型滤波器的一个特定频率严格相对应，则有。已知模拟滤波器，可以利用双线性变换函数bilinear将其变换为数字滤波器，调用格式为[bz,az]=bilinear(b,a,Fs)其中b,a分别为模拟滤波器系统函数分子、分母多项式系数；Fs为采样频率；bz、az为数字滤波器系统函数的分子、分母多项式系数。设计时要注意模拟原型低通频率预畸，否则衰减指标不能满足设计要求。FIR数字滤波器设计1、设计原理FIR滤波器的设计问题，就是要使所设计的FIR滤波器的频率响应逼近所要求的理想滤波器的频率响应。逼近可在时域进行，也可在频域进行。窗函数法设计FIR数字滤波器是在时域进行的，用窗函数截取无限长的，这样得到的频率响应逼近于理的频率响应。2、设计流程(1)给定希望逼近的频率响应函数；(2) 求单位脉冲响应；()(3)由过渡带宽及阻带最小衰减的要求，可选定窗形状，并估计窗口长度N。设待求滤波器的过渡带用表示，它近似等于窗函数主瓣宽度。因过渡带近似与窗口长度成反比，，决定于窗口形式；(4) 计算所设计的FIR滤波器的单位脉冲响应；(5) 由求FIR滤波器的频率响应，检验是否满足设计要求。一旦选取了窗函数，其指标（过渡带宽、阻带衰减）就是给定的。所以由窗函数设计FIR滤波器就是由阻带衰减指标确定用什么窗，由过渡带宽估计