基于Matlab数字信号滤波处理系统设计.docVIP

绪论 课题背景 在数字信号处理中，数字滤波占有极其重要的地位。目前对数字滤波器的设计有多种方法。其Matlab软件已成为设计数字滤波器的强有力工具。该软件是 1984年由美国MathWorks公司推出的一套用于数值计算及图形处理的高性能的可视化软件，它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便友好的界面和用户环境，深受工程技术人员及科技专家的欢迎，并很快成为应用学科计算机辅助分析、设计、仿真、教学等领域不可或缺的基础软件传统的数字滤波器设计过程复杂、计算工作量大、滤波特性调整困难,但利用Matlab信号处理工具箱（signal processing toolbox） XXXXX 数字滤波器简介 数字滤波器是指输入、输出均为数字信号，并对输入信号起到滤波作用的系统 。信号经滤波器处理，在时域就是信号与滤波器的冲击响应相卷积，即输入、输出关系为y(n)=x(n)×h(n)；而在频域则是信号频谱与滤波器频响的乘积，即。数字滤波器实质上是一个由有限精度算法实现的线性移不变离散时间系统，它的基本工作原理是利用离散系统特性对系统输入信号进行加工、处理和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输出。与模拟滤波相比，数字滤波具有很多突出的优点，如:它可以满足滤波器对幅度和相位特性的严格要求，可以避免模拟滤波无法克服的电压漂移、温度漂移和噪声等问题。数字滤波器按其单位样值响应的性质可分为有限冲激响应滤波器FIR finite impulse response和无限冲激响应滤波器IIR infinite impulse response。其中FIR滤波器具有以下优点： ①具有严格的线性相位，同时又具有任意的幅度特性；②单位样值响应是有限长的，因而系统一定是稳定的；③可以用快速傅里叶变换FFT算法来实现FIR滤波，从而可大大提高运算效率。 滤波器设计原理 数字滤波器是数字信号处理的重要组成部分，在许多科学技术领域中广泛应用。滤波器的功能是把输入序列通过一定的运算，变成输出序列。数字滤波器可以用差分方程来表示，即： （1） 对应的系统函数为： （2） 由式2知，若分母多项式的系数,那么该系统就变成FIR系统 ,即: （3） 显然，系数即为该系统的单位抽样响应h(0)，h(1)，，h(M)，且当 n M时，h(n)≡0。 由式3知，FIR系统只有零点，因此系统总是稳定的。另外，FIR系统还有易实现线性相位、允许设计多通带（或多阻带）滤波器、硬件容易实现等优点。 数字滤波器一般可用两种方法实现：一种是根据描述数字滤波器的数字模型或信号流图，用数字硬件构成专用的数字信号处理机；另一种是编写滤波运算程序，在计算机上运行。 滤波器设计与仿真 滤波器设计指标 以低通滤波器为例，设计指标如下： 采样频率: Fs=80KHz 通带：0≤f≤15KHz，幅度失真：±2%，（=0. 02） 阻带：20KHzf，衰减大于40dB，（=0. 01）。 经计算得相应数字滤波器的指标： 通带截止频率：； 阻带最小衰减：，（阻带波纹=0.01）。 滤波器设计方法 调用Matlab信号处理工具箱函数，分别采用窗函数法、频率采样法和切比雪夫最佳一致逼近法（调用 remez函数）进行设计。在用窗函数法和切比雪夫最佳一致逼近法设计时采用相同的阶次，以便于对结果进行比较。 确定滤波器的阶次 调用 remezord函数的程序清单如下： Clear f=[0.375 0.5 ]; A=[1 0]; dev=[0.02 0.01]; [n f0 A0 W]=remezord(f,A,dev) 运行结果如下： n =28, f0 = [0 0. 375 0. 5000 1. 0000],A0 = [1 1 0 0],W=[1 2] 窗函数法 程序清单如下: n=28; Wn=0.375; window=chebwin(n+1,30); b=fir1(n,Wn,window); [h,w]=freqz(b,1,256,1); hr=abs(h); h=20*log10(h); plot(w,h),grid; xlabel(频率（归一化）); ylabel(幅度(dB)); 运行程序后，幅频特性曲线如图 1所示，h(n)数据如下： -0.0067 0.0023 0.0084 0.0043 -0.0105 -0.0179 0.0000 0.0295 0.0291 -0.0205 -0.0713 -0.0383 0.1099 0.2933