2014级数字信号处理课程设计报告.docxVIP

2014级数字信号处理课程设计报告 题目： FIR Polyphase Realization 1. 题目要求： Task 1: 对3.4kHz的单频正弦信号进行8kHz采样，利用内插零点方式将其采样率升至1MHz； Task 2: 设计内插滤波器 Task 3: 分析用直接型实现结构和多相实现结构分别实现滤波器时对字长的敏感性，设计对滤波器系数进行量化的方案。 Task 4: 得到量化后的系数，画出系数量化后的滤波器的实际幅频响应曲线，和未量化时滤波器幅频响应曲线进行比较，分析字长对幅度响应误差的影响； 2. 课程设计思路： 涉及到的知识点 内插的基本原理及内插之后信号频谱的变化。 高倍数变采样率信号处理中，FIR滤波器的直接实现结构与多相实现结构之间的性能对比。 直接型结构和多相实现结构的有限字长效应。 方案分析 方案一：用matlab自带函数interp实现内插，fdatool实现滤波器，对其进行多相结构实现，对不同字长的滤波器，输入同一个随机序列，观察输出的变化来选择更好的结构和量化方案。 方案二：由内插的原理先进行时域插零，频域滤波来实现内插，fdatool实现滤波器，对其先进行级联型分解，再多相结构展开，由零点图观察字长的选取对幅频响应的影响，最终确定量化方案，画出实际幅频响应。 方案三：由内插的原理先进行时域插零，频域滤波来实现内插，fdatool实现滤波器，对其进行多相结构实现，对不同字长的滤波器，输入同一个随机序列，观察输出的变化来选择更好的结构和量化方案。给出合适的实际幅频响应曲线。 分析上面三种方案，可以看出，方案一没法更直观的观察内插的基本原理和时域变化；方案二中由零点图来分析不同结构对字长的敏感性效果不好；因此采用方案三，既能直观看出内插的原理又能很好的确定量化方案。 设计框架和流程 具体设计过程： 代码 时域内插零点输出 频域观察内插过程 分组实现多相结构 同一输入观察输出的差别 二进制n位截断 生成的图与数据等 内插原理的实现 图（1-1） 图（1-2） 设计内插滤波器： 图（2-1） 图（2-2）  = 3 \* GB3 \* MERGEFORMAT ③直接型实现结构和多相实现结构分别实现滤波器时对字长的敏感性 均方差 N=7N=8N=9N=10N=11N=12N=13N=14N=15N=16N=17N=18多相型0.52550.35140.30390.25950.23270.20690.18030.15120.13890.13250.12920.1276直接型0.52550.35270.30060.23190.17440.11590.06980.03440.01660.00850.00420.0021 表3 图（3-1） 字长n=9 图（3-2） 字长n=12 图（3-3） 字长n=16 图（3-4） 幅频响应曲线 图（4-1） 图（4-2） 对数据的分析  = 1 \* GB3 \* MERGEFORMAT ①内插原理的实现： 从图（1-2）可以看出在时域内插L-1个零点之后，频域上会出现频谱压缩，镜像频率成分也出现，经过fdatool设计的低通滤波器后，频谱只剩下原本想要的成分，同时图（1-1）上看，也得到了升采样率后的波形，跟理想波形很接近，同时分析可得由于设计的是FIR滤波器，冲击响应为有限长的因果序列，因此，进行时域卷积是，输出会相应的有延迟，这跟时域图中很好的对合上（第一张图可以看出内插采样波形要比实际1M输出延迟一定单位）。  = 2 \* GB3 \* MERGEFORMAT ②内插滤波器的设计： 利用fdatool设计等纹波低通滤波器，参数选择分别是 设计的低通滤波器如图（2-1）所示，图（2-2）展现了通带的细节，可以看出，在通带，幅频响应还是很平稳，等纹波振动的幅度也不太大，满足所需要求。  = 3 \* GB3 \* ME