解释代码%循环卷积 x1=[1 2 3 2]; x2=[2 1 2 1]; xn1=length(x1); xxn1=0:xn1-1; subplot(2,1,1); stem(xxn1,x1,'.','b'); title('序列x1'); axis([-1,4,0,4]);grid; xn2=length(x2); xxn2=0:xn2-1; subplot(2,1,2); stem(xxn2,x2,'.','b'); title('序列x2'); axis([-1,4,0,4]);grid; figure(2) N=4; x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; m=0:N-1; x=zeros(N,N); for n=0:N-1 x(:,n+1)=x2(mod((n-m),N)+1); end; yn=x1*x; subplot(3,1,1); stem(m,yn,'b','.'); title('序列x1和序列x2的4点循环卷积结果'); axis([-1,4,0,15]);grid; N=7; x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; m=0:N-1; x=zeros(N,N); for n=0:N-1 x(:,n+1)=x2(mod((n-m),N)+1); end; yn=x1*x; subplot(3,1,2); stem(m,yn,'b','.'); title('序列x1和序列x2的7点循环卷积结果'); axis([-1,7,0,15]);grid; N=8; x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; m=0:N-1; x=zeros(N,N); for n=0:N-1 x(:,n+1)=x2(mod((n-m),N)+1); end; yn=x1*x; subplot(3,1,3); stem(m,yn,'b','.'); title('序列x1和序列x2的8点循环卷积结果'); axis([-1,8,0,15]);grid;
时间: 2024-04-05 22:30:16 浏览: 100
03_代码实现卷积操作1
这段代码实现了两个序列x1和x2的循环卷积,并将结果绘制在三个子图中。首先,定义了两个序列x1和x2,然后计算它们的长度,并使用stem函数绘制它们的离散数据点。接下来,定义了一个长度为N的卷积长度,并使用零填充将x1和x2的长度扩展为N。然后,定义了一个长度为N的序列m,以及一个大小为N×N的矩阵x。在for循环中,使用循环矩阵将x2序列进行循环移位,并将移位后的序列赋值给矩阵x的第n+1列。最后,使用x1和x的矩阵乘法得到卷积结果yn,并使用stem函数将结果绘制成离散的数据点。在三个子图中分别展示了x1和x2的4点、7点和8点循环卷积结果,同时使用title函数给子图添加了标题,使用axis函数设置了坐标轴范围,使用grid函数添加了网格线。
阅读全文