序列线性卷积与圆周卷积相关运算

实验目的

本实验结合理论教材中有关离散傅里叶变换的教学内容，通过编程、上机调试程序，进一步增强使用计算机解决问题的能力。掌握线性卷积与圆周卷积软件实现的方法，并验证两者之间的关系。

实验问题 

实验代码 

1.function y=cirshfit(x,m,N);
2.if length(x)>N
3.error;
4.end
5.x = [x, zeros(1,N - length(x))];
6.n = 0: N-1;
7.n = mod(n-m,N);
8.y = x(n+1);
9.
10.function y=circonvtim(x1,x2,N)
11.n=0:N-1;
12.x1=[x1,zeros(1,N-length(x1))];
13.x2=[x2,zeros(1,N-length(x2))];
14.x3=x2(mod(-n,N)+1);
15.for m=0:N-1
16.x4=cirshfit( x3 ,m,N);
17.x5=x1 .* x4;
18.y(m+1)=sum( x5 );
19.end
20.
21.clc;close all;clear all;
22.x1 = [2,1,1,2];
23.x2 = [1,-1,-1,1];
24.y1=circonvtim(x1,x2,4);
25.y2=circonvtim(x1,x2,7);
26.y3=circonvtim(x1,x2,8);
27.figure;
28.subplot(4,1,1);stem(y1);title('N=4');axis([0,8,-2.5,2.5]);
29.subplot(4,1,2);stem(y2);title('N=7');axis([0,8,-2.5,2.5]);
30.subplot(4,1,3);stem(y3);title('N=8');axis([0,8,-2.5,2.5]);
31.y4=conv(x1,x2);
32.subplot(4,1,4);stem(y4);title('线性卷积');axis([0,8,-2.5,2.5]);

由实验结果可得，当N为7时，在N点区间内有相同的线性卷积和循环卷积，即当圆周卷积长度N>=N1+N2-1时, N点圆周卷积能取代线性卷积。而当N大于7时，需要补零。