数字信号处理实验五.docxVIP

附录1： 贵州大学实验报告 学院： 电气工程学院 专业： 测控技术与仪器 班级：测仪131 姓名 杨凯 学号 1308040019 实验组 实验时间 2016.5.20 指导教师 成绩 实验项目名称 FIR数字滤波器设计与软件实现 实验目的 （1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。 掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器 实验要求 （1）对两种设计FIR滤波器的方法（窗函数法和等波纹最佳逼近法）进行分析比较，简述其优缺点。 （2）附程序清单、打印实验内容要求绘图显示的曲线图。 （3）分析总结实验结果。 （4）简要回答思考题。 实验原理 对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率（）为周期进行周期延拓。公式为： 采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的 频谱不产生频谱混叠。 利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。 理想采样信号和模拟信号之间的关系为： 对上式进行傅立叶变换，得到： 在上式的积分号内只有当时，才有非零值，因此： 上式中，在数值上＝，再将代入，得到： 上式的右边就是序列的傅立叶变换，即 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用代替即可。 频域采样定理的要点是： 对信号x(n)的频谱函数X(ejω)在[0，2π]上等间隔采样N点，得到 则N点IDFT[]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为： 由上式可知，频域采样点数N必须大于等于时域离散信号的长度M(即N≥M)，才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x(n),即=x(n)。如果NM，比原序列尾部多N-M个零点；如果NM，z则=IDFT[]发生了时域混叠失真，而且的长度N也比x(n)的长度M短，因此。与x(n)不相同。 在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。 对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。 实验仪器 PC机，MATLAB软件 实验步骤 （1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理； （2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示； 图10.5.1 具有加性噪声的信号x(t)及其频谱如图 （3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。 （4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。 （4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。 eq \o\ac(○,2)采样频率Fs=1000Hz，采样周期T=1/Fs； eq \o\ac(○,3)根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz，阻带截至频率fs=150Hz，换算成数字频率，通带截止频率，通带最大衰为0.1dB，阻带截至频率，阻带最小衰为60dB。] eq \o\ac(○,4)实验程序框图如图10.5.2所示，供读者参考。 Fs=1000 Fs=1000，T=1/Fs xt=xtg 产生信号xt, 并显示xt及其频谱 用窗函数法或等波纹最佳逼近法 设计FIR滤波器hn 对信号xt滤波：yt=fftfilt(hn,xt) 1、计算并绘图显示滤波器损耗函数 2、绘图显示滤波器输出信号yt End 图10.5.2 实验程序框图 实验内容 %《数字信号处理(第三版）学习指导》第