基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告

基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告

ID:41050631

大小:531.00 KB

页数:13页

时间:2019-08-15

基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告_第1页
基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告_第2页
基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告_第3页
基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告_第4页
基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告_第5页
资源描述:

《基于.matlab和FPGA联合仿真的FIR低通滤波器设计报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、FIR低通滤波器设计报告1.设计内容本设计是基于FPGA的一个FIR低通滤波器设计,给定一段有高频干扰的信号Hnoise.wav,要求使用matlab设计出一个低通滤波器滤除其干扰频率,并取出10000到10160点进行时频分析。然后使用verilog语言编写出滤波器,联合modelsim进行编译仿真,并将结果与matlab结果进行对比。2.设计原理FIR滤波器响应(简称FIR)系统的单位脉冲响应为有限长序列,系统函数在有限z平面上不存在极点,其运算结构中不存在反馈支路,即没有环路。如果的长度为N,则它的

2、系统函数和差分方程一般具有如下形式:根据差分方程直接画出FIR滤波器的结构,称为直接型结构。如图所示:图2.1FIR滤波器直接结构FIR滤波器的特点:单位脉冲响应序列为有限个;可快速实现;可得到线性相位;滤波器阶数较高。对线性时不变系统保持线性相位的条件是:单位脉冲响应为偶对称或奇对称。即:为设计线性滤波器,应保证h(n)为对称的。1)若N为偶数,其线性相位FIR滤波器的对称结构流图:图2.2若N为偶数线性相位FIR滤波器的对称结构流图图中:“+1”对应偶对称情况,“-1”对应奇对称情况。当n为奇数时,最

3、后一个支路断开。2)若N为奇数,其线性相位FIR滤波器的对称结构流图:图2.3N为奇数线性相位FIR滤波器的对称结构流图在本设计中,我们采用线性FIR低通滤波器,所采用的阶数N=8,所以是偶对称的,估采取图2.2的结构,其中“±1“取“+1”。1.设计思路首先要用matlab对给定的信号进行时频分析来确定干扰信号的中心频率,然后再设计出相应的数字滤波器进行滤波,最后要在FPGA上实现FIR滤波器。首先要确定滤波器的抽头系数。其系数的确定,我们可以通过两种办法来实现:第一种就是通过matlab编写FIR滤波

4、器程序,然后直接导出抽头系数“h(n)”,另外一种办法就是使用matlab自带的FDATOOL简便地设计一个FIR滤波器,然后导出系数。考虑到要更直观地描述FIR滤波器的设计,我采用了第二种方法,用fdatool设计滤波器并导出参数。设计好滤波器后,接着要用matlab对给定的信号进行采样,保存,量化,导出以供FPGA的FIR滤波器模块使用。做好准备工作后,就可以在QUARTUSII里面编写Verilog代码了。由于是硬件描述语言,所以设计的思路很简单,就是通过把输入序列移位,然后首位对称相加再乘以抽头系

5、数,然后把相乘结果再相加最后给输出。其中涉及的难点是FPGA对有符号小数的乘法处理部分。在QUARTUSII编写好模块之后,就要用到modelsim来对设计进行仿真。对于仿真信号的输入,我们直接使用系统函数$readmemh()来读取matlab产生的信号采样值,然后通过时钟的驱动一个一个输入到FIR滤波器模块,最后在modelsim中显示滤波器的输出结果。验证的最后一步,就是要把modelsim输出的仿真结果用matlab进行绘图,和之前matlab程序的输出结果进行比较,从而验证滤波器的设计是否成功。

6、1.设计过程1)对信号进行时频分析语音信号是非平稳的随机过程,频率会随着时间而变化。在语音信号处理研究领域,通常用短时间的傅里叶变换(STFT)来分析其时域和频率特性,也就是在20ms~50ms以内,可以将语音信号看成是平稳的随机过程。平稳的随机过程是功率信号,可以用功率谱来描述。所加的噪声干扰是单音干扰,单音干扰是周期信号。这里提取了20ms的语音信号来分析其时域和频域特性,由于采样频率是8000Hz,因此我们利用了160个采样点,起始点在第10001点,结束点在第10160点。图4.1滤波前的时域图图

7、4.2滤波前的频谱图通过上图可以看出信号在3800HZ左右的地方受到严重的干扰。1)使用matlab设计FIR滤波器用matlab设计线性FIR滤波器,首先要确定其指标,在本设计里,我们规定滤波器的指标如下:阶数N=34,抽样频率fs=8000Hz,截止频率为3300用窗函数设计FIR滤波器,窗函数选择汉明窗。其幅频响应如下:图4.3滤波器幅频响应可以导出滤波器抽头系数h(n)为:h=0.000780227960895479-0.001771093775299170.00203379258009801-0

8、.000494245608574501-0.003344650598584810.00754651554336203-0.00762305727075547-0.0002838373332175980.0145004055957633-0.02556763039525060.02040310781520230.00783243571565760-0.05010586343181080.0793218675548351-0.05

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。