《用matlab实现fir滤波器》由会员分享，可在线阅读，更多相关《用matlab实现fir滤波器（8页珍藏版）》请在金锄头文库上搜索。
1、求用matlab实现fir滤波器（低通、高通，带通）浏览次数：1059次悬赏分：10|提问时间：2009-3-24 10:39|提问者：linyuchery520要求输入信号带有噪声，无论用那种设计方法实现，只要滤波效果好就行了。急急急推荐答案窗函数设计低通滤波器：fp=1000;fc=1200;as=100;ap=1;fs=22000;wp=2*fp/fs；wc=2*fc/fs;N=ceil(as-7.95)/(14.36*(wc-wp)/2)+1;beta=0.1102*(as-8.7);win dow=Kaiser(N+1,beta);b=fir1(N,wc,w in dow);freq
2、z(b,1,512,fs);咼通滤波器：fs=22000;Wp=2*5000/fs；Ws=2*4800/fs；Ap=1;As=100;N=ceil(8*pi/(Wp-Ws)+1;N=N+mod(N+1,2)+1;Wc=(Wp+Ws) /2/pi;h=fi门(N,Wc,high);omega=li nspace(0,pi,512);freqz(h,1,omega);带通滤波器：fs=22000;Wp1=2*1200/fs；Wp2=2*3000/fs；Wc1=2*1000/fs；Wc2=2*3200/fs；Ap=1；As=100；W1=(Wp1+Wc1)/2;W2=(Wp2+Wc2)/2;wdth
3、=m in (Wp1-Wc1),(Wc2-Wp2);N=ceil(11*pi/wdth)+1;b = fir1(N,W1 W2);freqz(b,1,512,fs)双线性法设计低通滤波器：fp=1000;fc=1200;as=100;ap=1;fs=22000;wp=2*fp/fs；wc=2*fc/fs;%归一化截止频率n,wn =ellipord(wp,wc,ap,as);%求数字滤波器的最小阶数和归一化截止频率b,a=ellip( n, ap,as,w n);%求传递函数的分子分母系数freqz(b,a,512,fs);基于MATLAB的FIR数字滤波器的设计与实现 期刊门户-中国期刊网2
4、009-2-11来源:黑龙江科技信息2008年9月下供稿文/李洋 洋 江亮亮导读简要介绍了 FIR滤波器的设计基础和利用窗函数设计FIR滤波器的方法，运用MATLAB语言实现了低通FIR滤波器的设计仿真。摘 要：简要介绍了 FIR 滤波器的设计基础和利用窗函数设计 FIR 滤波器的方 法，运用MATLAB语言实现了低通FIR滤波器的设计仿真。关键词：MATLAB； FIR滤波器；窗函数设计1 MATLAB 简介MATLAB是“矩阵实验室”(MATrix LABoratoy)的缩写，它是由美国Mathworks 公司于1984年正式推出的，是一种以矩阵运算为基础的交互式程序 语言，专门针对科学、
5、工程计算及绘图的需求1。随着版本的不断升级，内容 不断扩充，功能更加强大，从而被广泛应用于仿真技术、自动控制和数字信号处 理等领域。全球现有超过五十万的企业用户和上千万的个人用户，广泛分布在航 空航天、金融财务、机械化工、电信、教育等各个行业2。2 数字滤波器简介滤波器的种类很多，从功能上可以分为低通、高通、带通和带阻滤波器， 上述每种滤波器又可以分为模拟滤波器和数字滤波器。如果滤波器的输入输出都 是数字信号，则这样的滤波器称之为数字滤波器，它通常通过一定的运算关系改 变输入信号所含频率成分的相对比例或者滤除某些频率成分来实现滤波。根据数 字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即
6、无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器3。有数字信号处理的一般理论 可知，IIR滤波器的特征是具有无限持续时间的冲激响应，而FIR滤波器使冲激 响应只能持续一定的时间。3 数字滤波器的设计步骤3.1确定指标。在设计一个滤波器之前，必须首先确定一些技术指标。这些 技术指标需要来制定。在很多实际应用中，例如语音或音频信号处理中，数字滤 波器常用来实现选频操作。因此，指标的形式一般确定为频域中幅度和相位的响 应。3.2逼近。确定了滤波器的技术指标后，就可以利用数学和 DSP 的基本原理 提出一个滤波器模型来逼近给定的目标。这一步是滤波器设计的所要研究的主要 问题，是整个设计过
7、程的重中之重。3.3硬件或软件实现。上两步的结果是得到以差分或系统函数或冲激响应描 述的滤波器，根据这个描述就可以用硬件或在计算机上用软件实现。4 窗函数设计法4.1 窗函数的理论根据。窗函数设计滤波器的基本思想，就是根据给定的滤 波器技术指标，选择滤波器的阶数N和合适的窗函数(n)。即用一个有限长 度的窗口函数序列 (n)来截取一个无限长的序列hd (n)获得一个有限长序 列h(n)，即h (n) =w (n) *hd (n)，并且要满足以下两个条件：(1)窗谱 主瓣尽可能地窄，以获得较陡的过渡带；(2)尽量减少窗谱的最大旁瓣的相对幅 度，也就是能量尽量集中于主瓣，使峰肩和纹波减小，就可增多
8、阻带的衰减。这 就给窗函数序列的形状和长度选择提出了严格的要求。4.2用窗函数设计 FIR 滤波器的步骤(1) 根据技术要求确定待求滤波器的单位取样响应hd (n)。如果给出待求滤波 器的频响为Hd (ejw),那么单位取样响应应用下式求出：hd (n)=Hd(ejw)ejwndw(1)如果Hd(ejw)较为复杂，或者不能用封闭公式Hd(ejw)表示时，则不能用上 式求出hd (n)。我们可以对Hd(ejw)从w=0到w=2n采样M点，采样值为 Hd (e)，k=0，1，2，M-1并用代替公式(1)中的dw，公式(1)近似 写成：hM (n)=Hd (e)e(2)根据频率采样定理，hM (n)
9、与hd (n)应满足如下关系：hM (n)=hd (n+rM) 因此，如果M选的较大，可以保证在窗口内hM (n)有效逼近hd (n)。实际计算公式(2)，可以用Hd (ejw)的M点采样值，进行M点的IDFT (IFFT) 得到。如果给出通带阻代衰减和边界频率要求，可选用理想滤波器作为逼近函数，从而 用理想滤波器的特性作傅里叶逆变换，求出hd (n)。例如，理性低通滤波器如 公式( 3)所示，Hd (ejw)=ejwa，w?燮 wcO, wcw?燮 n (3)求出单位取样响应hd (n)如下式：hd (n) =为保证线性相位，取 a=( N-1 ) /2。(2) 根据对过渡带及阻带衰减的要求
10、，选择窗函数的形式，并估计窗口长度N。 设待求滤波器的长度用Aw来表示，它近似等于窗函数主瓣宽度。因过渡带Aw 近似于窗口长度N成反比，N=A/Aw, A决定于窗口形式。例如，矩形窗A=4n 哈明窗A=8n等。按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证 阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。(3) 计算滤波器的单位取样响应h (n)，h (n)=w (n)*hd (n)。式中w (n) 是上面选择好的窗函数。如果要求线性相位，则要求hd (n)和w (n)均对(N-1) /2对称。如要求h (n)对(N-1)/2奇对称，只要保证hd (n)对(N-1)/2奇对 称就可以
11、了。(4) 验算技术指标是否满足要求。设计出的滤波器频率响应应用下式计算： H (ejw)=h (n)e-jwn计算上式时可用FFT算法。如果H (ejw)不满足要求，根据具体情况重复(2)(3) (4)步，直到满足要求。5用各种窗函数设计FIR数字滤波器实例例分别用矩形窗和Hamming窗设计线性相位FIR低通滤波器。要求 通带截止频率wc=n/4，单位脉冲响应h (n)的长度N=21。绘出h (n)及其幅 频响应特性曲线。设计分析(1)用窗函数法设计 FIR 数字滤波器时，先求出相应的理想滤波器(本例应为 理想低通)单位脉冲响应hd(n)，再根据阻带最小衰减选择合适的窗函数w(n), 最后
12、得到 FIR 滤波器单位脉冲响应 h(n) =w(n) *hd(n)。(2) 本题中，wc=n/4, N = 21,所以线性相位理想低通滤波器的单位脉冲响应 为：hd (n) =(3) 为了满足线性相位FIR滤波器条件h (n) =h (N-1-n),要求a= (N-1) /2=10。信号处理工具箱中有窗生成函数boxcar，hamming，hanning和blackman等。Wn=boxcar (m)产生长度为m的矩形窗函数列向量wn,其他窗函数产生工具箱 函数的调用格式相同。(4) MATLAB程序clear； close allN=21 ； wc=pi/4；%理想低通滤波器参数n=0:N
13、-1 ； r= (N-1) /2；hdn=sin (wc* (n-r) /pi./ (n-r)；%计算理想低通单位脉冲响应hd (n)if rem (N，2) =0 hdn (r+1) =wc/pi； end % N 为奇数时，处理 n=r 点的0/0型wn 1=boxcar (N)；% 矩形窗hn 1=hdn.*wn 1；% 加窗% 以上两条语句可代以 fir 函数：hn 1=fir1 (N-1,wc/pi,boxcar (N)；wn2=hamming (N)；% hamming 窗hn2=hdn.*wn2；% 加窗% 以上两条语句可代以 fir 函数：hn2=fir1 (N-1，wc/pi
14、，hamming (N)； k=3； k=1；%绘图函数p751f的位置参数hnwplot (hn 1，k)； title (矩形窗设计的 h (n)%hamming 窗wn2=hamming (N)；hn 2=hd n.*wn 2；k=3；hnwplot (hn2,k)； title (hamming 设计的 h (n)(5) 程序运行结果(6)结果分析对两种窗函数的设计结果分别如上图所示。由图中可以看出，不同的窗 函数生成的过渡带宽和阻带最小衰减是不同的。这就是选择窗函数的根据。FIR 滤波器的应用十分广泛，当今的许多信号处理系统，如图像处理等系统要求具有 线性相位特性。在这方面FIR滤波
15、器有独特的优点，我们能很容易设计出具有 严格线性相位的FIR系统。此外，FIR滤波器的冲激响应应是有限长序列，其系 统函数为一个多项式，它所包含的极点多为原点，所以FIR滤波器永远是稳定 的。结束语FIR数字滤波器在数字处理领域有非常重要的地位，应用MATLAB语 言进行FIR数字滤波器的设计仿真时，可根据设计要求随时改变参数，以使滤 波器达到最优化。参考文献1 陈怀琛.MATLAB及在电子信息课程中的应用（第三版）M电子工业出版社,2006.3.2 中仿信息科技有限公司.MATLAB产品介绍EB/OL.http:/ 丁玉美数字信号处理（第二版）M西安：西安电子科技大学出版社，2000.12. 作者简介：李洋洋，女， 22岁，华中师范大学信息技术系电子信息工程
