信号的频谱分析实验matlab,数字信号处理 Matlab实验四 FFT的频谱分析

XXXX大学XXXX学院

实验名称 FFT的频谱分析

实验目的：

加深理解FIR数字滤波器的时域特性和频域特性，掌握FIR数字滤波器的设计原理与设计方法，以及FIR十足滤波器的应用。

实验内容：

FIR 数字滤波器总是稳定的系统，且可以设计成具有线性相位的。其在数据通信、图像处理、语音信号处理的实际应用领域中得到广泛的应用。N阶有限冲击响应(FIR)数字滤波器的转移函数为：H(z)= h(n)z-n，系统的单位

n=0N-1冲击响应h(n)是长度为N的有限长因果序列。当满足好h(n)=h(N-n-1)的对称条件时，改FIR数字滤波器具有线性相位。FIR数字滤波器的设计方法主要有窗函数法和频率采样法。

实验步骤

1. 窗函数法

FIR滤波器的冲击响应就是系统函数各次项的系数，所以设计FIR滤波器的方法之一就是：从时域出发，截取有限长的一段冲击响应作为H(z)系，冲击响应长度N就是函数H(z)的阶数。主要N足够长，并且截取的方法合理，总能够满足频域的要求。这就是FIR滤波器的窗函数设计法。

2. 频率采样法

频域采样法是先对理想频响Hd(ej )进行采样，得到样值H(k),再利用插值公式直接求出系统转换函数H(z)，以便实现;或者求出频响Hd(ej )，以便与理想频响进行比较。在[0,2 ]区间上对Hd(ej )进行N点采样，等效于时域以N为周期延拓。 实验结果及分析：

例4-1

实验代码

wp=0.4*pi;ws=0.6*pi;deltaw=ws-wp;%过渡带宽计算

N0=ceil(6.6*pi/deltaw)+1;%按海明窗计算所需的滤波器阶数N0