数字信号处理大作业报告

数字信号处理大作业序列谱分析及FFT快速卷积

学院（系）：软件学院

专业：软件工程

学生姓名：许利志

学号：201492300

班级：软网1404

完成日期：2016.4.28

大连理工大学

Dalian University of Technology

报告内容

一、设计题目

1.序列的谱分析及特性

2.用FFT实现快速卷积

二、设计目的及意义

序列的谱分析及特性

对于三种典型序列------单位采样序列、实指数序列、矩形序列，要求：

(1). 画出以上序列的时域波形图；

(2). 求出以上序列的傅里叶变换；

(3). 画出以上序列的幅度谱及相位谱，并对相关结果予以理论分析；

(4). 对以上序列分别进行时移，画出时移后序列的频谱图，验证傅里叶变换的时移性质；

(5). 对以上序列的频谱分别进行频移，求出频移后频谱所对应的序列，并画出序列的时域波形图，验证傅里叶变换的频移性质。

用FFT实现快速卷积

(1).加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好地利用FFT 进行数字信号处理。

(2).理解快速卷积在实际生活中的重要作用，加深对知识的理解。

三．设计题目的描述以及要求

序列的谱分析及特性

（1）画出序列的时域波形图；

（2）求出序列的傅里叶变换；

（3）画出序列的幅度谱及相位谱，并对相关结果予以理论分析；

（4）对序列进行时移，画出时移后序列的频谱图，验证傅里叶变换的时移性质；

用FFT实现快速卷积

（1）FFT的出现，使DFT得到广泛应用。在实际应用中，为了分析时域离散LTI系统或者序列滤波时，需要计算两个序列的线性卷积。为了提高运算速度，可以利用FFT来实现。

（2）设计并编写程序来实现线性卷积运算，可以在重叠相加法和重叠保留法中任选一种方法实现

（3）给出输入信号和输出信号的图形描述，以及简要的说明以及计算中间过程的图形描述及简要说明，每计算完一段后输出分段卷积段号、该段卷积结果、该段图像并暂停，按任意键继续下一段。

四、设计原理

序列的谱分析及特性

1.时移原理在这个序列运算中，x[n]的每一个样本都移动（即延迟）k个采样周期，移位后的序列为y(n)。当k >0时每一个样本向右移动，称为x(n)的延时序列；当k<0时，每一个样本向左移动，称为x(n)的超前序列：y(n)=x(n-k) 。

2.频移原理

3.傅里叶变换（DFT)原理

离散傅里叶变换的结果为有限长和离散的，它实质上是对序列傅里叶变换在频域均匀离散的结果，因而使数字信号处理可以在频域采用数字运算的方法进行，大大增加傅里叶变换的灵活性和使用性。离散傅里叶变换的定义如下

其中为旋转因子NW，N为变换区间长度。

用FFT实现快速卷积

MATLAB 中计算序列的离散傅里叶变换和逆变换是采用快速算法，利用fft 和ifft 函数实现。

（1）.

输入参数：为待计算DFT 的序列，N 为序列的长度。

输出参数：为序列的IDFT。

（2）.

输入参数：xk为待计算IDFT 的序列，N 为序列xk的长度。

输出参数：Xn 为序列xk的IDFT。

五、设计内容以及实现方法

序列的谱分析及特性：

1.单位采样序列的谱分析及特性实现

2.实指数序列的谱分析及特性实现

3.矩阵序列的的谱分析及特性实现

以下列出第一部分代码（其他见电子文档）：

clear

n=1:50;

x=zeros(1,50);

x(1)=1;

subplot(3,1,1);

stem(x);

axis([0 50 -1 1]);

title('单位采样序列');

N=25;

k=-N:N;

X=x*(exp(-j*pi/25)).^(n'*k);

magX=abs(X);

subplot(3,1,2);

stem(magX);

title('单位采样序列的幅度谱');

angX=angle(X);

subplot(3,1,3);stem(angX);

title ('单位采样序列的相位谱');

n=1:50; x=zeros(1,50); t=10; x(t)=1;

figure

subplot(3,1,1);

stem(x);title('单位采样序列的时移');

k=-25:25; X=x*(exp(-j*pi/25)).^(n'*k);

magX=abs(X);

subplot(3,1,2);

stem(magX);title('单位采样序列时移的幅度谱');

angX=angle(X);

subplot(3,1,3);

stem(angX);

title ('单位采样序列时移的相位谱');

n=1:50; x=zeros(1,50);

x(1)=1; l = 5; y = exp(-j*pi/25*l).^n;

z=x.*y;

figure

subplot(3,1,1);

stem(z);

title('频移后单位采样序列');

k=-25:25; X=z*(exp(-j*pi/25)).^(n'*k);

magX=abs(X);

subplot(3,1,2);stem(magX);

title('频移后单位采样序幅度谱');

angX=angle(X); subplot(3,1,3);

stem(angX);

title ('频移后单位采样序列的相位谱');

用FFT实现快速卷积（代码及结果见第六部分）

六、设计结果及改进建议

序列的谱分析及特性：

单位采样序列的生成及傅立叶变换的幅度谱单位采样序列的时移特性如下图4.2所示相位谱图，如上图（1）单位采样序列，