数字信号处理96142new.docVIP

数 字 信 号 处 理 一、实验要求 编程利用FFT进行卷积计算，通过实验比较出快速卷积优越性。写出实验报告。 二、实验原理 利用FFT进行离散卷积的步骤归纳如下： (1)、设x(n)的列长为N1，h(n)的列长为N2，要求 y(n)=x(n)*h(n)=x(k)h((n-k))NRN(n)= [1] (2)、为使两有限长序列的线性卷积可用其圆周卷积来代替而不产生混淆，必须选择N≥N1+ N2-1。为使用基-2FFT来完成卷积计算，故要求N=2v(v是整数)。用补零的办法使x(n)，h(n)具有列长N，即 x(n)= h(n)= (3)为用圆周卷积定理计算线性卷积，先用FFT计算x(n) ，h(n) 的N点离散傅立叶变换 x(n) X(k) [2] h(n) H(k) [3] (4)组成卷积 Y(k)=X(k)H(k) [4] (5)利用IFFT计算Y(k)的离散傅立叶逆变换得到线性卷积y(n)。由于 y(n)=WN-nk=[WNnk]* [5] 可见，y(n)可由求(1/N)Y*(k)的FFT再取共轭得到。 三、实验题目 两个正弦序列的卷积（均为两个周期，256点） 输入序列： 卷积输出： （2）正弦序列与三角序列的卷积（正弦序列为两个周期，256点；三角序列为一个周期，256个点） 输入序列： 卷积输出： （3）两个矩形序列的卷积（均为两个周期，256个点，占空比0.5） 输入序列： 卷积输出： （4）单位冲击与正弦波（单位冲击序列为256个点，正弦序列为1.7个周期，256个点）输入序列： 卷积输出： 任何序列与单位冲击序列的卷积为原序列，所以结果正确。 (5)正弦序列与矩形序列的卷积（两序列均为256个点，正弦序列为两个周期，矩形序列为两个周期，占空比为0.2） 输入序列： 卷积输出： 三、实验代码 （1）两个正弦序列的卷积（均为两个周期，256点） k=1:256; s1=sin(k/64*pi); s2=s1; xk=fft(s1,2*length(k)-1); subplot(3,1,1) stem(k,s1); title(正弦序列1) subplot(3,1,2) stem(k,s2); title(正弦序列2) yk=fft(s2,2*length(k)-1); rm=ifft(xk.*yk); m=(-255):(255); subplot(3,1,3) stem(m,rm) xlabel(m); ylabel(rm); title(卷积输出) （2）正弦序列与三角序列的卷积（正弦序列为两个周期，256点；三角序列为一个周期，256个点） k=1:256; s1=sin(k/64*pi); s2=[1/64*k(1:63),-k(64:192)/64+2,1/64*(k(193:256)-192)-1] xk=fft(s1,2*length(k)-1); subplot(3,1,1) stem(k,s1); title(正弦序列) subplot(3,1,2) stem(k,s2); title(三角序列) yk=fft(s2,2*length(k)-1); rm=ifft(xk.*yk); m=(-255):(255); subplot(3,1,3) stem(m,rm) xlabel(m); ylabel(rm); title(卷积输出) （3）两个矩形序列的卷积（均为两个周期，256个点，占空比0.5） k=1:256; y=ones(size(k)); s1=[y(1:64),-1.*y(65:128),y(129:192),-1.*y(193:256)]; s2=s1; xk=fft(s1,2*length(k)-1); subplot(3,1,1) stem(k,s1); title(矩形序列1) subplot(3,1,2) stem(k,s2); title(矩形序列2) yk=fft(s2,2*length(k)-1); rm=ifft(xk.*yk); m=(-255):(255); subplot(3,1,3) stem(m,rm) xlabel(m); ylabel(rm); title(卷积输出) （4）单位冲击与正弦波（单位冲击序列为256个点，正弦序列为1.7个周期，256个点） k=0:255; y=ones(size