FIR和IIR数字滤波器的设计方法.docVIP

FIR和IIR数字滤波器的设计方法 数字滤波器和模拟滤波器差别： 1、数字滤波器 　　数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。 　　数字滤波器对信号滤波的方法是:用数字计算机对数字信号进行处理，处理就是按照预先编制的程序进行计算。它的核心是数字信号处理器。 2、模拟滤波器 模拟滤波器分为无源滤波器和有源滤波器，其中无源滤波器由R、L、C组成，有源滤波器由集成运放和R、C组成，不需要使用电感。集成运放的开环电压增益和输入阻抗均很高，输出电阻小，构成有源滤波电路后还具有一定的电压放大和缓冲作用。但集成运放带宽有限，所以目前的有源滤波电路的工作频率难以做得很高。 模拟滤波器的设计，就是用模拟系统的系统函数Ha(s)去逼近所需要的理想特性。标准的模拟低通滤波器的设计公式有巴特沃斯和切比雪夫等，它们都是根据幅度平方函数来确定的。 为逼近理想低通滤波器，其模拟理想低通滤波器的幅度特性可用幅度平方函数表示，即 式中，Ha(s)为所设计的模拟滤波器的系统函数，它是s的有理函数；Ha(jΩ)是其稳态响应，即滤波器频率特性| Ha(jΩ)|为滤波器的稳态振幅特性。 FIR滤波器的设计方法 FIR滤波器的单位脉冲激励响应h(n)是有限长(0≤n≤N-1)，其Z变换为： H(z)=sum(h(n)z-m,0≤m≤N-1) FIR滤波器一般采用间接法设计，常用的方法有窗函数法、频率采样法和切比雪夫等波纹逼近法。用窗函数法设计FIR滤波器的具体步骤为： 根据对阻带衰减及过渡带的指标要求，选择窗函数的类型（矩形窗、三角窗、汉宁窗、汉明窗、凯塞窗等），并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择相同阶数下主瓣窄的窗函数； 构造希望逼近的频率响应函数； 计算h(n)； 加窗得到设计结果。 IIR滤波器的设计方法 设计IIR数字滤波器一般采用直接法和间接法，间接法有脉冲响应不变法和双线性变换法，应用最广泛的是双线性变换法。其设计方法是，先设计过渡模拟滤波器得到系统函数H(s)，然后将H(s)按某种方法转换成数字滤波器的系统函数H(z)，具体步骤如下： 确定所需类型数字滤波器的技术指标，通带截至频率、通带最大衰减、阻带截至频率、阻带最小衰减； 将所需类型的数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Q=2/T tan(0.5w)； 将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标； 设计模拟低通滤波器； 通过频率变换将模拟低通转换成相应类型（高通，带通，带阻）的过渡模拟滤波器； 采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Matlab实现方法： FIR滤波器 以低通滤波器为例，其常用的设计指标有： 通带边缘频率fp(数字频率为Ωp) 阻带边缘频率fs(数字频率为Ωs) 通带内最大纹波衰减δp=-20log10(1-αp)，单位为dB， 阻带最小衰减δs=-20log10(αs)，单位为dB 通带峰值起伏αp=1-10-0.05δp 阻带起伏αs=10-0.05δs 在matlab中，可以利用函数fir1设计FIR滤波器，利用函数freqz画出各滤波器的频率响应。采用[n,Wn,beta,ftype] = kaiserord(f,a,dev)函数来估计滤波器阶数等，得到凯塞窗滤波器。 　　这里的函数kaiserord(f,a,dev)或者kaiserord(f,a,dev,fs)，f为对应的频率，fs为采样频率；当f用数字频率表示时，fs则不需要写。（比如低通滤波器，f=[fp fs]，a=[1 0]为由f指定的各个频带上的幅值向量，dev=[αp αs]；对于带通滤波器，f=[fs1 fp1 fp2 fs2]，a=[0 1 0]，dev=[αs αp αs]）。ωp=2*fp/fs，ωs=2*fs/fs，Wn=(ωp+ωs)/2，由此设计FIR滤波器得： b=fir1(n,Wn,’ftype’,window,’normalization’)设定滤波器和窗函数的类型，其中参数’ftype’是 high for a highpass filter with cutoff frequency Wn. stop for a bandstop filter, if Wn = [w1 w2]. The stopband frequency range is specified by this interval. DC-1 to make the first band o