dsp带通滤波器设计程序.docVIP

设计要求 试用DSP设计FIR滤波器，分别实现带通的功能，具体要求如下： ① 滤波器的阶数≥5，截止频率自行选定，滤波系数用MATLAB确定。 ② 编制C54XDSP实现FIR滤波器的汇编源程序。 ③ 用软件仿真器完成上述程序的模拟调试。 ④ 以数据文件形式自行设定滤波器输入数据，以数据文件形式输出滤波结果，并与输入数据进行比较分析。用软件仿真器有关工具显示FIR滤波器的输入输出波形，以证明滤波器滤波性能。 FIR滤波器的基本原理 数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所期望的输出序列。一个线性位移不变系统的输出序列和输入序列之间的关系，应满足常系数线性差分方程。FIR滤波器的差分方程为： FIR滤波器的传递函数为： 直接由差分方程得出的实现结构如图１所示： 图１　横截型 直接型﹑卷积型 FIR滤波器的结构图 由上面的公式和结构图可知，FIR滤波算法实际上时一种乘法累加运算。它不断地从输入端读入样本值，经过（）后做乘法累加，输出滤波结果 FIR滤波器的设计 FIR滤波器的设计方法主要有窗函数法和频率采样法。其中，窗函数法是最基本的方法。本次设计采用窗函数法。 一般是先给定所要求的理想滤波器频率响应，由导出，我们知道理想滤波器的冲击响应是无限长的非因果序列，而我们要设计的是是有限长的FIR滤波器，所以要用有限长序列来逼近无限长序列，设： ３-1 常用的方法是用有限长度的窗函数w n 来截取 即： ３-２ 这里窗函数就是矩形序列RN n ,加窗以后对理想低通滤波器的频率响应将产生什么样的影响呢?根据在时域是相乘关系,在频域则是卷积关系： ３-３ 其中,为矩形窗谱, 是FIR滤波器频率响应通过频域卷积过程看的幅度函数H ω 的起伏现象,可知,加窗处理后，对理想矩形的频率响应产生以下几点影响： 使理想频率特性不连续点处边沿加宽，形成一个过渡带，其宽度等于窗的频率响应的主瓣宽度。 在截止频率的两边的地方即过渡带的两边，出现最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，则取决于旁瓣的多少。改变N，只能改变窗谱的主瓣宽度，改变ω的坐标比例以及改变的绝对值大小，但不能改变主瓣与旁瓣的相对比例 此比例由窗函数的形状决定 。对窗函数的要求a、窗谱主瓣尽可能窄，以获取较陡的过渡带；b、尽量减小窗谱的最大旁瓣的相对幅度；即能量集中于主瓣，使肩峰和波纹减小，增大阻带的衰减。 clear all; f1 1500; f2 4000; f3 7000; fs 30000; N 256; T 1/fs; n [0:N-1]; x1 0.08*sin 2*pi*f1*n*T ; x2 0.08*sin 2*pi*f2*n*T ; x3 0.08*sin 2*pi*f3*n*T ; xn ceil 32768* x1+x2+x3 ; figure 1 ; plot xn figure 2 yff abs fft xn ; df n* fs/N ; plot df,yff ; fid fopen indata.inc,w ; fprintf fid, .word%5.0f\n,xn ; fclose fid ; 部分输入信号数据为 共256个采样点 ： .word 0 .word 6097 .word -431 .word 1541 .word -1019 .word 1 .word -2493 .word -5507 ... 2.基于Matlab工具箱Signal中的fir1函数设计FIR带通滤波器，选择Hamming窗函数法，进而获取FIR带通滤波器的系数b，并保存为firin.inc数据文件。这里FIR滤波器的阶数N 32。 实现代码如下： w1 3000/15000*2; w2 5000/15000*2; b fir1 32,[w1,w2],bandpass ; fp fopen firin.inc,wt ; fprintf fp, .word%20.0f\n,b*32768 ; fclose fp ; 滤波器系数为： B0 -4,B1 -0,B2 7, B3 183,B4 -109,B5 -488,B6 328,B7 495,B8 -195,B9 260, B10 -1180,B11 -1433,B12 4086,B13 1879,B14 -7292,B15 -876,B16 8709, B17 -876,B18 -7292,B19 1879,B20 4086,B21 -1433,B22 -1180,B23 260, B24 -195,B25 495,B26 328,B27 -488,B28 -109,B29 183,B30 7,B31