用matlab实现fir滤波器.doc

求用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); window=Kaiser(N+1,beta); b=fir1(N,wc,window); freqz(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=fir1(N,Wc,high); omega=linspace(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=min((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,wn);%求传递函数的分子分母系数 freqz(b,a,512,fs); 基于MATLAB的FIR数字滤波器的设计与实现 期刊门户-中国期刊网2009-2-11来源:《黑龙江科技信息》2008年9月下供稿文/李洋洋 江亮亮 [导读]简要介绍了FIR滤波器的设计基础和利用窗函数设计FIR滤波器的方法，运用MATLAB语言实现了低通FIR滤波器的设计仿真。 摘 要：简要介绍了FIR滤波器的设计基础和利用窗函数设计FIR滤波器的方法，运用MATLAB语言实现了低通FIR滤波器的设计仿真。关键词：MATLAB；FIR滤波器；窗函数设计1 MATLAB简介 MATLAB是“矩阵实验室”（MATrix LABoratoy）的缩写，它是由美国Mathworks公司于1984年正式推出的，是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求[1]。随着版本的不断升级，内容不断扩充，功能更加强大，从而被广泛应用于仿真技术、自动控制和数字信号处理等领域。全球现有超过五十万的企业用户和上千万的个人用户，广泛分布在航空航天、金融财务、机械化工、电信、教育等各个行业[2]。2 数字滤波器简介 滤波器的种类很多，从功能上可以分为低通、高通、带通和带阻滤波器，上述每种滤波器又可以分为模拟滤波器和数字滤波器。如果滤波器的输入输出都是数字信号，则这样的滤波器称之为数字滤波器，它通常通过一定的运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分来实现滤波。根据数字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器[3]。有数字信号处理的一般理论可知，IIR滤波器的特征是具有无限持续时间的冲激响应，而FIR滤波器使冲激响应只能持续一定的时间。3 数字滤波器的设计步骤 3.1确定指标。在设计一个滤波器之前，必须首先确定一些技术指标。这些技术指标需要来制定。在很多实际应用中，例如语音或音频信号处理中，数字滤波器常用来实现选频操作。因此，指标的形式一般确定为频域中幅度和相位的响应。 3.2逼近。确定了滤波器的技术指标后，就可以利用数学和DSP的基本原理提出一个滤波器模型来逼近给定的目标。这一步是滤波器设计的所要研究的主要问题，是整个设计过程的重中之重。 3.3硬件或软件实现。上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器，根据这个描述就可以用硬件或在计算机上用软件实现。4 窗函数设计法 4.1 窗函数的理论根据。窗函数设计滤波器的基