数字信号处理课程设计实验内容.docVIP

PAGE 华 北 电 力 大 学 实 验 报 告 201003010308黄日辉 PAGE 第 PAGE 18 页 华 北 电 力 大 学 实 验 报 告 实验环境 MATLAB 6.5 实验名称 实验一： FFT的应用 实 验 目 的 1、熟悉MATLAB在数字信号处理中的应用。 2、掌握利用FFT计算序列线性卷积的基本原理及编程实现。 3、掌握对连续信号进行采样的基本原理和方法，并利用FFT对信号进行频谱分析。 实 验 原 理 1.线型卷积和圆周卷积 设x(n)为L点序列，h(n)为M点序列，x(n)和h(n)的线性卷积为 y yl x(n)和h(n)的N点圆周卷积为 y 圆周卷积与线性卷积相等而不产生交叠的必要条件为 N≥L+M-1 圆周卷积定理：根据DFT的性质，x(n)和h(n)的N点圆周卷积的DFT等于它们DFT的乘积 DFT 2.快速卷积 快速卷积算法用圆周卷积实现线性卷积，根据圆周卷积定理利用FFT算法实现圆周卷积。可以将快速卷积的步骤归纳如下： (1)为了使线性卷积可以用圆周卷积来计算，必须选择 N≥L+M-1; 同时为了能使用基-2FFT完成卷积运算，要求N=2v 实 验 原 理 和h(n)的长度均为N。 (2)计算x(n)和h(n)的N点FFT x h (3)组成乘积 Y (4)利用IFFT计算Y(K)的IDFT，得到线性卷积y(n) Y 3.分段卷积 我们考察单位取样响应为h(n)的线性系统，输入为x(n),输出为y(n)，则 y 当输入序列x(n)极长时，如果要等x(n)全部集齐时再开始进行卷积，会使输出相对输入有较大的延时，再者如果序列太长，需要大量存贮单元。为此我们把x(n)分段，分别求出每段的卷积，合在一起得到最后总的输出。这种方法称为分段卷积分段卷积可细分为重叠保留法和重叠相加法。 重叠保留法：设x(n)的长度为Nx，h(n)的长度为M。我们把序列x(n)分成多段N点序列x 利用FFT实现重叠保留发的步骤如下： 在x(n)前面填充M-1个零，扩大以后的序列为 xn={0,0,…,0,M-1 将x(n)分为若干N点子段，设L=N-M+1为每一段的有效数据长度，则第 实 验 原 理 i段xi(n)(0nN-1 x 计算每一段与h(n)的N点圆周卷积，利用FFT计算圆周卷积： x h Y Y 舍去每一段卷积结果的前M-1个样本，连接剩下样本，得到卷积结果y(n)。 重叠相加法：设h(n)长度为M，将信号x(n)分解成长为L的子段，建议L选择与的M数量级相同，以xi x x x 每一段卷积yi 利用FFT实现重叠保留法的步骤如下： (1)将x(n)分为若干L点子段xi 实 验 原 理 (2)计算每一段与h(n)的卷积yi (3)将各段yi y 4．利用FFT进行频谱分析 若信号本身是有限长的序列，计算序列的频谱就是直接对序列进行FFT运算求得，就代表了序列在之间的频谱值。 幅度谱 相位谱 在MATLAB中，提供了计算模值（或绝对值）的函数abs和计算相角的函数angle，其用法如下： abs(x)；%当x为实数时，计算x的绝对值；当x为复数时，计算x的模值； angle(x); %计算复矢量或复矩阵x的相角，结果为介于之间的值。 若信号是模拟信号，用FFT进行谱分析时，首先必须对信号进行采样，使之变成离散信号，然后就可按照前面的方法用FFT来对连续信号进行谱分析。按采样定理，采样频率应大于2倍信号的最高频率，为了满足采样定理，一般在采样之前要设置一个抗混叠低通滤波器。用FFT对模拟信号进行谱分析的方框图如下所示。 抗混叠低通滤波器 抗混叠低通滤波器 采样T=1/fs N点FFT 实 验 内 容 1.对于两个序列：x(n)=nR16(n),h(n)=R8(n) （1）在同一图形窗口中绘出两序列的时域图形。 （2）利用FFT编程计算两序列的线性卷积，绘出的时域图形。 2．利用FFT对信号进行谱分析 对于连续信号xa(t)=cos(2πf1t)+5cos(2πf2t)+cos(2πf3t) ，其中f1=6.5kHz, f2=7kHz, f3=9kHz, 以采样频率fs=32 kHz对其进行采样， （1）对xa(t)信号采集16点样本，分别作16点和补零到256点的FFT，并分别绘出对应的幅频特性曲线。 （2）对xa(t)信号采集256点样本，分别作256点和512点的FFT，并分别绘出对应的幅频特性曲线。 （3）比较（1）和（2）中的结果，分析采样点数和傅里叶变换点数对FFT的影响，说明