利用MATLAB实现循环卷积.docxVIP

一、实验目的 利用 MATLAB 实现循环卷积。 比较循环卷积与线性卷积的区别。 二、实验条件 PC 机，MATLAB7.0 三、实验内容 循环卷积的定义：两个序列的 N 点循环卷积定义为： [h(n) ? x(n)] ? ?N ?1 h(m)x(n ? m) (0 ? n ? N ) N N k ?0 利用 MATLAB 实现两个序列的循环卷积可以分三个步骤完成： (1)初始化：确定循环点数 N，测量输入 2 个序列的长度。 循环右移函数：将序列 x(n)循环右移，一共移 N 次（N 为循环卷积的循环次数）， 最后将每次循环成的新序列组成一个矩阵 V。 相乘：将 x(n)移位后组成的矩阵 V 与第二个序列 h(n)对应相乘，即得循环卷积结果。程序如下： 程 序 一 ： clear;close all; N=10; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N); x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,.); title(循环卷积的结果); xlabel(n);ylabel(y1(n)); 运行后所得图形如下： 观察所得的循环卷积结果发现并没有呈现周期性的序列，因此将程序做下列改变。程序二： clear;close all; N=40; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,x2,x2,x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N); x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,.); title(循环卷积的结果); xlabel(n);ylabel(y1(n)); 从图中可以看出循环卷积的结果已经呈循环序列。将程序进一步改进——在 x2 序列之间加一些零矩阵；程序如下： clear;close all; N=50; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,zeros(1,7),x2,zeros(1,7),x2,zeros(1,7),x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N); x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,.); title(循环卷积的结果); xlabel(n);ylabel(y1(n)); 线性卷积与循环卷积的区别： 从循环卷积的定义公式中可以看出，循环卷积和线性卷积的不同之处在于：①两个 N 点序列的 N 点循环卷积的结果仍为 N 点序列，而两个 N 点序列线性卷积的结果的长度则变为为 2N-1；②循环卷积对序列的移位采取循环移位，而线性卷积对序列则是 采取线性位移。而就是这两点不同，导致循环卷积和线性卷积有不同的结果和性质。 然而虽然循环卷积和线性卷积虽然是不同的概念，但是它们之间有一个有意义的 公式联系在一起 y(n) ? [h(n) ? x(n)] N ? ?? ??? r ??? ? ? y(n ? rN )?G ?? N ? (n) ，其中 y(n) ? h(n)* x(n) ; 也就是说，两个序列的N 点循环卷积是他们的线性卷积以 N 为周期的周期延阔。设 序列 h(n)的长度为 N ,序列 x(n)的长度为 N ，此时，线性卷积结果的序列的点数为 1 2 N ? N ? N -1;因此如果循环卷积的点数 N 小于 N ? N ?1 ，那么上述周期性延阔的 1 2 1 2