[理学]ch2_2DFT计算卷积.ppt

fftfilt(h,x,n) h: FIR filter x: input sequence n为DFT点数，一般取2的整数次幂。 利用MATLAB实现分段卷积 * 利用MATLAB实现分段卷积 % Generate the noise sequence N = 64; d = rand(N,1)-0.5; % Generate the uncorrupted sequence and add noise for m = 1:1:N, s(m) = 2*(m-1)*((0.9)^(m-1)); x(m) = s(m) + d(m); end % the length of moving average filter M = 4; % Generate the moving average filter coefficients h = ones(1,M)/M; % Perform the overlap-add filtering operation y = fftfilt(h,x,8); * 4点滑动平均系统去噪结果 * 2. 重叠保留法（overlap save） 方法： (1) 将x[k]长序列分段，每段长度为L。 (2) 各段序列xn[k]与 M点短序列h[k]循环卷积。 (3) 从各段循环卷积中提取线性卷积结果。 前M-1个点不是线性卷积的点， 因yn[k]=xn [k] L h[k] 故分段时，每段与其前一段有M-1个点重叠。 * 2.4.2 长序列和短序列的线性卷积 2. 重叠保留法（overlap save） - - x [ k ( M 1)] M - 1 - - L ( M 1) L - 1 x 0 [ k ] x 1 [ k ] 2 L - M k 第一段前需补M-1个零 * 2.4.2 长序列和短序列的线性卷积 2. 重叠保留法（overlap save） 记：yn[k] =xn [k] L h[k] * 2.4.2 长序列和短序列的线性卷积 2. 重叠保留法（overlap save） y0[k]中的[M-1, L-1]点对应于线性卷积x[k]*h[k]中的[0 , L-M]点。 y1[k]中的[M-1, L-1]点对应于线性卷积x[k]*h[k]中的[ L-(M-1), 2L-M-(M-1)]点。 * 2.4.2 长序列和短序列的线性卷积 y[k]={2, 7, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 41, 14} 解: 重叠相加法 x1[k]={2, 3, 4, 5, 6} x2[k]={7, 8, 9, 10, 11} x3[k]={12,13, 14} y1[k]={2, 7, 12, 16, 20, 17, 6} y2[k]={ 7, 22, 32, 36, 40, 32, 11} y3[k]={12, 37, 52, 41, 14} 例3: 已知序列x[k]=k+2,0?k?12, h[k]={1,2,1}，试分别利用重叠相加法和重叠保留法计算线性卷积, 取L=5 。 * 解: 重叠保留法 y[k]={2, 7, 12, 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 41, 14} x1[k]={0, 0, 2, 3, 4} x2[k]={3, 4, 5, 6 ,7} x3[k]={6 ,7 , 8, 9, 10} y1[k]= x1[k]?h[k]= {11, 4, 2, 7, 12} x4[k]={9, 10 , 11, 12,13} y2[k]= x2[k]?h[k]= {23, 17, 16, 20, 24} y3[k]= x3[k]?h[k]= {35, 29, 28, 32, 36} y4[k]= x4[k]?h[k]= {47, 41, 40, 44, 48} x5[k]={12,13, 14, 0, 0} y5[k]= x5[k]?h[k]= {12, 37, 52, 41, 14} 例3: 已知序列x[k]=k+2,0?k?12, h[k]={1,2,1}，试分别利用重叠相加法和重叠保留法计算线性卷积, 取L=5 。 * 数字信号处理(Digital Signal Processing) 信号与系统系列课程组 国家电工电子教学基地 第2章 离散傅里叶变换(DFT) * 2.1 有限长序列的傅里叶分析 2.2 离散傅里叶变换的性质 2.3 离散傅里叶变换与