数字信号处理课程设计基于重叠相加法的圆周卷积

数字信号处理课程设计基于重叠相加法的圆周卷积内容摘要：

nji()函数计算圆周卷积 （ 6）各段 重叠相加 （ 7） 取出实际的输出序列 武汉理工大学《数字信号处理》 课程设计 说明书 6 程序设计流程图 函数 juanji()流程图 开始 对 x1(n)进行 FFT运算 对 x(2)进行 FFT运算 频域相乘 Yk=X1k*X2k 对 Yk 进行反变换得到时域结果y(n) 结束 武汉理工大学《数字信号处理》 课程设计 说明书 7 chongdie（ x,h,N） 函数流程图 方案 1: 开始 x(n) 按 长 度N 进行分段 确定圆周卷积的周期 L 将 序列 按要求填充 确定分段数 K 调用函数 juanji() 将各段处理的结果重叠相加输出最终结果 结束 武汉理工大学《数字信号处理》 课程设计 说明书 8 方案 2: 方案比较 方案 1 和方案 2 的区别在于：在方案 1 中由已知的长序列 x(n)分段大小 N、短序列长度 M 依据公式计算 L N+M1 圆周卷积周期 L。 而在方案 2 中由已知的圆周卷积周期L、短序列长度 M 依据公式计算长序列 x(n)的分段大小。 比较这两种方案可以知道，方案 1 中如果 L 直接取 L=N+M1，则 L 的值不一定是2 的 n 次方，而接下来的要进行的 L 点快速傅里叶变换中， L 的值需要为 2 的 n 次方才能正确计算。 所以方案 1 的必须要求取恰当的 N 值以使 L 点满足要求。 方案 2 中可以直接取 L 为恰当的值从而满足快速傅里叶变换的要求，所以方案 2 比方案 1 的适用范围更大。 开始 取圆周卷积的周期 L 计算每一分段的大小 N 将 序列 按要求填充 计算 分段数 K 调用函数 juanji() 将各段处理的结果重叠相加输出最终结果 结束 武汉理工大学《数字信号处理》 课程设计 说明书 9 3程序代码 主程序： x1=[1,2,3,4,5,5,4,3,2,1]。 x2=[8,5,4]。 L=8。 y=chongdie(x1,x2,L)。 s=1:10。 subplot(2,2,1)。 stem(s,x1) r=1:3。 subplot(2,2,2)。 stem(r,x2) u=1:13。 subplot(2,2,3)。 stem(u,y) 函数 juanji() function y=juanji(x1,x2,L) if length(x1)L %如果 x1 长度大于 L 则产生错误 error(39。 L must not be less than length。