数字信号处理实验报告黎美琪 201300800610 13通信2实验一名称：周期卷积、循环卷积和线性卷积比较 一、实验目的1.理解周期卷积、循环卷积、线性卷积的定义2.用图像显示上述几种卷积并对其进行直观的比较 二、实验步骤 自行设定：）它们的线性卷积（）求它们的循环卷积（求它们的周期卷积（两个有限长序列3)8(2)8)1(2012,81,1129,1)(,2012,81,0129,8)(21==⎩⎨⎧≤≤≤≤-≤≤=⎩⎨⎧≤≤≤≤≤≤-=N N n n n n x n n n n n x实验代码：（大部分语句为图像显示处理）%循环卷积&线性卷积&周期卷积 %%线性卷积 figure(1);set(gcf, 'color', 'w')%将图的背景设置为白色x1=[zeros(1,8),[1:4],zeros(1,4),zeros(1,8)];%原有限长序列x1（n ） x2=[zeros(1,8),ones(1,4),zeros(1,4),zeros(1,8)] ; %原有限长序列x2（n ） L=length(x1)%长度L M=length(x2)%长度My1=conv(x1,x2) %线性卷积 subplot(311) stem(x1);title('有限长序列x1（n ）') axis([1 L 0 5]) subplot(312) stem(x2);title('有限长序列x2（n ）') axis([1 M 0 1]) subplot(313) stem(y1);grid on ; title('线性卷积')axis([1 L+M-1 0 11]) %%循环卷积(圆周卷积) figure(2);set(gcf, 'color', 'w')%将图的背景设置为白色%x11=[[1:4],zeros(1,4),[1:4],zeros(1,4),[1:4],zeros(1,4)];x11=[[1:4],zeros(1,2),[1:4],zeros(1,2),[1:4],zeros(1,2),[1:4],zeros(1,2)];y2=conv(x2,x11)P=length(x22)%长度Psubplot(311);stem(x11);title('有限长序列x1的周期延拓x11（n）')axis([1 L 0 5])subplot(312)stem(x2);title('有限长序列x2（n）')axis([1 M 0 1])subplot(313)stem(y2);grid on;title('循环卷积')axis([1 P+M-1 0 11])%%周期卷积figure(3);set(gcf, 'color', 'w')%将图的背景设置为白色x22=[ones(1,4),zeros(1,4),ones(1,4),zeros(1,4),ones(1,4),zeros(1,4)]; y2=conv(x1,x22)Q=length(x22)%长度Qsubplot(311)%stem(x11);stem(x11);%title('有限长序列x1（n）')title('有限长序列x1的周期延拓x11（n）')axis([1 L 0 5])subplot(312);stem(x22);title('有限长序列x2的周期延拓x2（n）')axis([1 Q 0 1])subplot(313)stem(y2);grid on;title('周期卷积')%axis([1 L+Q-1 0 15])axis([1 P+Q-1 0 11])(一)线性卷积1.线性卷积步骤1）将序列x2（n）翻褶2）平行向右移位3）被卷积两序列对应序号值相乘，再相加X2(-m)00001111X2(1-m)0000111 1 Y(8)=1X2(2-m)000011 11 Y(9)=3X2(3-m)00001 111 Y(10)=6X2(4-m)0000 1111 Y(11)=10X2(5-m)000 01111 Y(12)=9X26-m)00 001111 Y(13)=7X2(7-m) 0 0001111 Y(14)=4X2(8-m) 00001111 Y(15)=0X2(9-m) 0000111 1 Y(6)=0X2(10-m) 000011 11 Y(17)=0X2(11-m) 00001 111 Y(18)=0X2(12-m) 0000 1111 Y(19)=0X2(13-m) 000 01111 Y(20)=0X2(14-m) 00 001111 Y(21)=0X2(15-m) 0 0001111 Y(22)=0注意：为方便比较几种不同卷积的结果，设定的序列的初始位置在n=9。

因为前面的平移相乘结果都为0，所以前面省略了一部分，这里列出的是主要部分，且x2（n-m）中的n是在8的基础上向右平移的位数。

3.线性卷积图像：(二)周期卷积基本原理：将h(n) 进行周期延拓，周期为N：∑∞-∞=+=rrNnhnh)()(~计算)(~n x 与)(~n h 的周期卷积)(~n y N ：∑∑∑∑∑∑∑∞-∞=∞-∞=-=∞-∞=-=-=-=+=-+=+-=-=-=r r N m r N m N m N m N rN n y m rN n h m x rN m n h m x m n h m x m n h m x n y )()]()([)()()(~)()(~)(~)(~110111.周期卷积步骤1）将两个主值序列都进行周期延拓得到x11（n ）和x22（n ） 2）对应序号相乘并相加求和 3）周期性重复 2.周期卷积列表X1(m) 12340000 y （n ） X2(m)11110000X11((m))8 12340000 12340000 12340000X22((m))8 111100001111000011110000X11(-(m))8 00004321 00004321 00004321 X11(1-(m))8 10000432 10000432 10000432 Y(9)=1 X11(2-(m))8 21000043 21000043 21000043 Y(10)=3 X11(3-(m))8 32100004 32100004 32100004 Y(11)=6 X11(4-(m))8 43210000 43210000 43210000 Y(12)=10 X11(5-(m)8 04321000 04321000 04321000 Y(13)=9 X11(6-(m))8 00432100 00432100 00432100 Y(14)=7 X11(7-(m))8 00043210 00043210 00043210Y(15)=4X11(8-(m))800004321 00004321 00004321（周期性重复） Y(16)=03.周期卷积的图像：基本原理：对于有限长序列x(n)和y(n)( 0<=n<=N-1 ) DFT[()]()DFT[()]()x n X k y n Y k ==若()()()F k X k Y k =10()IDFT[()]()(())()N N N m f n F k x m y n m R n -===-∑x(n)和y(n)的N 点循环卷积，记作()()n x n y ⊗，这个卷积可以看作是周期序列x （n ）和y （n ）做周期卷积后再取主值序列。

1.循环卷积步骤1）补零（如果两虚列长度不同，需要补零使两序列长度相同） 2）其中一个序列x1（n ）周期延拓为x2（n ） 3）x11（n ）翻褶，截取计算区域 4）循环移位5）被卷积两序列对应序号值相乘，再相加 6）取主值序列循环卷积长度N(8)>=N1(4)+N2(4)-1 循环卷积长度N(6)<=N1(4)+N2(4)-1三、分析总结1.对比N=8和N=6两种情况下的循环卷积结果：2.对比周期卷积、循环卷积、线性卷积的结果：周期卷积)(~n y N 是x(n)与h(n)的线性卷积y(n) 的周期延拓。

由于)(~n x 与)(~n h 的周期都为N ，因此它们的周期卷积)(~n y N 的周期也为N ，正好等于y(n)的长度，即上式中以N 为周期的周期延拓没有发生混叠，线性卷积y(n)正好是周期卷积)(~n y N 的一个周期。

循环卷积又是周期卷积的主值序列，因此，此时循环卷积yN(n)与线性卷积y(n)完全相同，即：∑-=-≤≤-===⊗=110)()()()()(~)()()(N m N N N N n m n h m x n y n R n y n h n x n y 四、学习体会通过此次实验深入了解了周期卷积、循环卷积、线性卷积三者之间的关系，且对其原理也有了更加深刻的理解。

通过这次实验为学会了一种新的思想：从比较中找出相同点和不同点，这样对概念的理解会更加深刻。

此次实验还遇到了一个问题：stem 图形都是从n=1开始画图的？尝试了多种方法也没能达到目的效果，虽然这个对实验结果没有很大的影响，但是用了多种方法没能成功，且花费了较多时间，没能抓住重点。

⎩⎨⎧-+≤≤-+≥*=*-+≥-+≥-+201)()()()(1)(11)()()(21212121211212111N N n N N N n x n x n x n x N N N N N n y N N N N N n y N n f n f N 能代表线性卷积点循环卷积时，即当循环卷积的长度。

周期延拓才无混叠现象为周期进行以时，所以只有当的长度为序列。

的周期延拓序列的主值为周期以是线性卷积点循环卷积可见，。