matlab和dspFIR数字滤波器设计方法讨论.docVIP

基于MATLAB和DSP的FIR数字滤波器设计方法的讨论 摘要：介绍了基于Matlab和DSP设计FIR数字滤波器的四种方法，并经过实验，将一个多频率成分的信号通过该滤波器并进行比对,实验结果表明,这四种方法都能有效地实现FIR滤波器的设计任务。 关键词：FIR； DSP；Matlab； CCSLink Abstract: This paper introduced four methods to design FIR digital filter based on MATLB and DSP. Through the experiment, a signal with multiple frequencies was fed into this filter, and a contrast was drawn. The result of this experiment demonstrates the efficiency of the four methods to fulfill the designing task of the FIR digital filter. Keywords: FIR; DSP; Matlab；CCSLink 1 前言 随着信息和数字技术的发展，数字信号处理已成为当今极其重要的的学科和技术领域之一，而数字滤波作为之中的一个重要环节，也一直是人们研究的重点。数字滤波器是一种用来过滤离散时间信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。与IIR滤波器相比，FIR的实现是非递归的，总是稳定的；更重要的是，FIR滤波器在满足幅频响应要求的同时，可以获得严格的线性相位特性。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。FIR滤波器的设计方法有许多种，如窗函数设计法、频率采样设计法和最优化设计法等。窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列FIR_TASK: STM #INPUT,INBUF_P STM #OUTPUT,OUTBUF_P STM #KS-1,BRC RPTBD LOOP-1 STM #N,BK LD *INBUF_P+,A FIR_FILTER: STL A,*DATA_P+% RPTZ A,#N-1 MAC *DATA_P+0%,*COEF_P+0%,A STH A,*OUTBUF_P+ LOOP: 程序中运用一个乘加指令MAC实现两列数据的卷积过程，该指令能在1个机器周期内完成对2个操作数的间接寻址(Graph(Time/Frequency…，按图1所示填入相应参数，看频域图要把Display type选择为FFT Magnitude。输入信号和滤波后输出信号的时域图和频域图如图2所示。 图1 观测信号波形的属性设置 图2 CCS环境下的输入输出信号时域频域图 2) Matlab环境下编写M脚本文件实现 MATLAB信号处理工具箱提供了各种窗函数、滤波器设计函数和滤波器实现函数。本文的滤波器设计及滤波程序如下：shuru(1:256)=sin((0:255)/fs*2*pi*500)+sin((0:255)/fs*2*pi*1000)+ sin((0:255)/fs*2*pi*2500); %产生长度为256的输入信号，包含500Hz、1000Hz和2500Hz的正弦波 xishu=fir1(16,1500/8000*2); %用Hamming窗得到16阶的滤波器系数 shuchu=conv(shuru,xishu); %用卷积指令conv对输入信号进行滤波 t=(0:255)/fs; w=8000*(0:127)/256; shuru1=fft(shuru,256); %把输入信号进行256点fft变换，用于后续画信号频谱图 shuchu1=fft(shuchu,256); %把输出信号也进行256点fft变换 subplot(2,2,1) %画出输入信号的时域图 plot(t,shuru);grid on;title(输入信号时域图); 以上是滤波器的设计以及滤波过程主要程序，画图程序简略。 程序执行结果得到的输入输出信号的时域频域图如图3所示： 图3 Matlab编程实现滤波前