设计的fir带阻滤波器.docVIP

PAGE 《DSP应用系统设计》综合设计报告 题目： 带阻数字滤波器的设计 姓名 XXX 学号 ＿XXXXXXX 手机 XXXXXXX Email HYPERLINK mailto:@ XXXXXXXXX 专 业 ＿ 自动化 班级 ＿XXXXXX 开课学期 2011 至 2012 学年 2＿学期 提交时间 2012 年 7 月 5 日 引言 数字滤波是现代数字信号处理的重要研究内容，在信号分析、图像处理、模式识别、自动控制等领域得到了广泛应用。按照数字滤波器的特性，它可以被分为线性与非线性、因果与非因果、无限长冲击响应（IIR）与有限长冲击响应（FIR）等等。与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。 基于DSP设计的带阻滤波器实现了对信号的滤波，使其在某一频率带内的信号被滤去而使其他频率的信号通过。Matlab语言是一种用于科学计算的高级语言，随着Matlab信号处理工具箱的不断完善和发展，使数字滤波器的计算机辅助设计得以实现。 二、实验目的 利用DSP的原理和方法设计FIR数字带阻滤波器（采用Hamming窗），具体设计参数如下： 采样频率Fs=8000Hz。 （2）滤波器阶数 L=128。 （3）输入信号为三个正弦波的叠加，其频率分别为f1=500Hz，f2=1500Hz，f3=2800Hz，幅值均为0.333。 （4）带阻滤波器的边界频率fp=[1000 3000]，则对应的归一化数字角频率为Wn=fp/(Fs/2)=[0.25，0.75]。 三、实验场地及仪器、设备和材料： （1） PC兼容机一台，操作系统为Windows2000。 （2）安装Code Composer Studio 4.2软件与Matlab软件。 （3）DSP系统实验箱（硬件连接使用，本实验未使用）。 四、实验内容 1 、系统方案 一个典型的DSP系统除了数字信号处理部分外，还包括A/D和D/A两部分。这是因为自然界的信号，如声音、图像等大多是模拟信号，因此需要将其数字化后进行数字信号处理，模拟信号的数字化即称为A/D转换。数字信号处理后的数据可能需还原为模拟信号，这就需要进行D/A转换。一个仅包括A/D和D/A两部分的简化数字信号处理系统功能如图1所示。 图1 简化数字信号处理系统功能框图 DSP系统的输入信号多种多样，但一般都要转换为模拟信号，经过抗混叠滤波后由模数转换芯片将模拟信号变换为数字信号。然后，DSP根据实际需要对其进行相应的处理，如FFT、卷积等;处理后的数字信号再由数模转换芯片将其转换为模拟值，在对其进行平滑滤波后最终输出模拟信号。需要说明的是，对输入模拟信号的采样频率要大于其本身频率的两倍，因为根据奈奎斯特采样定律，只有采样频率至少为输入信号频率的两倍时输入信号的信息才保证不会丢失。 设计带阻滤波器，其中N+1为滤波器阶数，b为返回的长度为N+1的系数向量。Wn为归一化边界频率，0 Wn 1.0，Wn 处增益为-6dB。Wn =1.0对应的非归一化频率为Fs/2，其中Fs为采样频率。 本设计中，设Fs=8000Hz，设计带通阻波器，格式为：b = fir1(N,Wn)或b = fir1(N,Wn, stop) 其中Wn = [W1 W2], W1、W2为边界频率，当带阻滤波器的边界频率fp=[1000， 3000]，则对应的归一化数字角频率为Wn=fp/(Fs/2)=[0.25，0.75]。执行b = fir1(N,Wn, stop)与b=round(b*32768) 指令得到的系数为： b=[0,0,-28,0,0,0,36,0,0,0,-52,0,0,0,77,0,0,0,-113,0,0,0,161,0,0,0,-223,0,0,0,304,0,0,0,-407,0,0,0,541,0,0,0,-718,0,0,0,965,0,0,0,-1336,0,0,0,1974,0,0,0,-3411,0,0,0,10415,0,16397,0,10415,0,0,0,-3411,0,0,0,1974,0,0,0,-1336,0,0,0,965,0,0,0,-718,0,0,0,541,0,0,0,-407,0,0,0,304,0,0,0,-223,0,0,0,161,0,0,0,-113,0,0,0,77,0,0,0,-52,0,0,0,36,0,0,0,-28,0,0] 注意: 上述滤波器系数为浮点格式，当DSP程序