DSP课程设计CLL.docVIP

目录 1.原理分析 3 1.1圆周卷积原理 3 1.2重叠相加法 6 1.3线性卷积、圆周卷积、重叠相加法、DFT、FFT之间的联系 7 2．程序设计 8 2.1程序设计思路 8 2.2程序设计流程图 9 2.3方案比较 10 3．程序代码 12 4.程序运行结果与分析 15 5.心得体会 16 6.参考文献 17 1.原理分析 1.1圆周卷积原理 对两个N点序列和，除了可以做线性卷积外，还有一种很重要的卷积运算，还有一种很重要的卷积运算，就是圆周卷积。 令 则圆周卷积结果长度不变, 为N. 由上式可以得出圆周卷积与周期卷积的关系，就是有限长序列圆周卷积结果的周期延拓，等于它们周期延拓后的周期卷积。也就是说，周期卷积的主值序列，是各周期序列主值序列的圆周卷积。 圆周卷积的实现步骤 若、分别是长度为N、M的序列则与线性卷积至多M+N-1个非零值，如果LM+N-1则周期延拓时必然会有一部分非零值发生混叠；只有当LM+N-1时，周期延拓才不会发生混叠。 之所以讨论用圆周卷积来计算线性卷积的条件，是因为圆周卷积可在频域下利用DFT求得，从而可采用DFT的快速算法FFT来计算，这样就可以利用FFT来计算线性卷积，大大提高运算效率。 1.2重叠相加法 在实际应用中利用FFT来计算两个序列的圆周卷积从而实现计算其线性卷积，但是常遇到的问题是参加卷积的两个序列的长度相差较大，这样长度小的序列就需要补很多的零点，这样就需要打的存储量，运算时间也会变长。所以常用到的解决方法有两种，其中一种就是重叠相加法。 h(n)长度为N，x(n)长度为无限长， x(n)取M点，且与N尽量接近 重叠相加法的卷积示意图 重叠相加法的步骤如下 （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)的重叠部分相加，最后得到结果为 1.3线性卷积、圆周卷积、重叠相加法、DFT、FFT之间的联系 由时域与频域的关系可知，两序列和在时域下进行线性卷积的结果等于这两个序列在频域下相乘后进行反变换回时域的结果。圆周卷积在一定条件下（LM+N-1）与线性卷积得到的结果相同，而圆周卷积可在频域下利用DFT求得，从而可采用DFT的快速算法FFT来计算，这样就可以利用FFT来计算线性卷积，大大提高运算效率。而在利用FFT来计算圆周卷积的过程中当两序列的长度相差较大时采用重叠相加法来进行计算可有效提高计算的效率，减小存储空间的消耗。 2．程序设计 由以上的原理分析可以设计程序利用FFT算法来实现基于重叠相加法的圆周卷积。 2.1程序设计思路 函数juanji(x1,x2,L)设计 x1(n)进行N点快速傅里叶变换得X1k x2(n)进行N点快速傅里叶变换得X2k 进行频域相乘Yk=X1k*X2k 对Yk进行反变换得到时域卷积结果y(n) 函数chongdie（x,h,N）设计 方案1： （1）首先取长序列x(n)进行分段的长度N，以使其分段后的长度与较短的相近 （2）确定圆周卷积的周期L （3）填充序列使得循环中对序列的索引不会超出范围 （4）确定分段数K （5）对序列进行分段调用juanji()函数计算圆周卷积 （6）各段重叠相加 （7）取出实际的输出序列 方案2： （1）首先取圆周卷积的周期L（即进行L点的快速傅里叶变换） （2）计算每一分段的大小N （3）填充序列使得循环中对序列的索引不会超出范围 （4）计算分段数K （5）对序列进行分段调用juanji()函数计算圆周卷积 （6）各段重叠相加 （7）取出实际的输出序列 2.2程序设计流程图 函数juanji()流程图 chongdie（x,h,N）函数流程图 方案1： 方案2： 2.3方案比较 方案1和方案2的区别在于：在方案1中由已知的长序列x(n)分段大小N、短序列长度M依据公式计算LN+M-1圆周卷积周期L.而在方案2中由已知的圆周卷积周期L、短序列长度M依据公式计算长序列x(n)的分段大小。 比较这两种方案可以知道，方案1中如果L直接取L=N+M-1，则L的值不一定是2的n次方，而接下来的要进行的L点快速傅里叶变换中，L的值需要为2的n次方才能正确计算。所以方案1的必须要求取恰当的N值以使L点满足要求。方案2中可以直接取L为恰当的值从而满足快速傅里叶变换的要求，所以方案2比方案1的适用范围更大。 3．程序代码 主程序： x1=[1,2,3,4,5,5,4,3,2,1]; x2=[1,0,1]; L=8; y=chongdie(x1,x2,L); s=1:10;subplot(2