数字信号课程设计报告用FFT实现快速卷积.docVIP

课 程 设 计 报 告 课程名称 __________数字信号处理_____________ 设计题目 _______用FFT实现快速卷积__________ 专 业 _____________通信工程______________  课 程 设 计 任 务 书 设计题目：_________用FFT实现快速卷积__________________ _________________________________________________________ 设计内容与要求： FFT的出现，使DFT在数字通信、语音信号处理、图像处理、功率谱估计、系统分析与仿真、雷达信号处理、光学、地震及数值分析等各个领域都得到广泛应用。然而，各种应用一般都以卷积和相关运算为依据。在实际应用中，为了分析时域离散LTI系统或者序列滤波时，需要计算两个序列的线性卷积。为了提高运算速度，可以利用FFT来实现。 要求：参考课本上第90页的内容（3.4.1 用DFT计算线性卷积），设计并编写程序来实现重叠相加法计算线性卷积。 课 程 设 计 评 语 成绩： 指导教师：_______________ 年 月 日 目 录? 第1章 设计任务 3 1.1 原理 4 1.1.2 算法基本思想 5 第2章 程序设计 7 2.1 程序设计思路 8 2.1.1 函数juanji(x1,x2,L)设计 8 2..1.2函数chongdie（x,h,N）设计 8 2.2 程序流程图 8 第3章 分析与测试 9 3.1循环卷积设计 10 3.1.1程序 10 3.1.2 测试 10 3.1.3 分析 12 3.2线性卷积设计 12 3.2.1程序 12 3.2.2测试 12 第4章 实验结果 14 4.1实验结果 14 第5章 总结 15 第6章 参考文献 16 第1章 设计任务 计算1个给定序列与输入序列的卷积。 功能： 对给定的数据进行卷积运算，要求分段卷积由循环卷积实现。要求设计有数据导入界面，各种参数从软件界面可以输入，其中给定序列可以由界面输入，对运算前后的数据绘制曲线。 要求： 初步完成总体设计，搭好框架，确定人机对话的界面，确定函数功能，控制参数的输入方法； 设计线性卷积的实现方案； 编写两序列作循环卷积的程序； 通过直接做线性卷积来检验最后结果。 设计步骤： 用结构化设计方法。一个程序划分成若干模块，每一个模块的函数功能要划分好，总体设计应画出流程图； 输入输出界面要友好； 源程序书写要规范，加必要的注释； 要提供通过直接卷积进行检验的结果； 程序一定要要能运行起来。 1.1 原理 1.1.1 算法产生背景 DFT 是连续傅里叶变换在时域和频域上都离散的形式，将时域信号的采样变换为在离散时DFT 具 对于线性非移变离散系统，可由线性卷积表示时域输入输出关系，即 y(n)x(n)*h(n) ? 通常采用循环卷积降低运算量，但实际中往往无法满足对信号处理的实时性要求。因此，产生了重叠相加法和重叠保留法两种典型的算法，用以快速计算线性卷积，成为了DFT 的一个 1.1.2 算法基本思想 重叠相加法是将待过滤的信号分割成长为 N 的若干段，，每一段都可以和有限时宽单位取样回应作卷积，再将过滤后的各段重叠相加。 在实际应用中利用FFT来计算两个序列的圆周卷积从而实现计算其线性卷积，但是常遇到的问题是参加卷积的两个序列的长度相差较大，这样长度小的序列就需要补很多的零点，这样就需要大的存储量，运算时间也会变长。所以常用重叠相加法来解决。 如以下情况： h(n)长度为N，x(n)长度为无限长 x(n)取M点，且与N尽量接近 可采用如下方法来解决 图1-1 重叠相加法的卷积示意图 重叠相加法的步骤如下 （1）将h(n)补零延长到L =M+ N -1，并计算长为L的FFT，得到 H(k)。 （2）分别将xk(n)补零延长到L =M+ N -1，并计算长为L的FFT，得到 Xk(k) （3）计算，并求长为L的反变换，即 （4）将yk(n)的重叠部分相加，最后得到结果为 第2章 程序设计 2.1 程序设计思路 2.1.1 函数juanji(x1,x2,L)设计 x1(n)进行N点快速傅里叶变换得X1k x2(n)进行N点快速傅里叶变换得X2k 进行频域相乘Yk=X1k*X2k 对Yk进行反变换得到时域卷积结果y(n) 2..1.2函数chongdie（x,h,N）设计 （1）首先取圆周卷积的周期L（即进行L点的快速傅里叶变换） （2