16t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 46 行

clc;
N=31;
n=[0:N];
x=n.*(stepseq(0,0,N)-stepseq(16,0,N));
y=stepseq(0,0,N)-stepseq(8,0,N);
X=fft(x);
Y=fft(y);
Z=X.*Y;
z=ifft(Z);
subplot(3,1,1);
stem(n,z);
title('线性卷积');
axis([0,25,0,100]);
N1=15;
n1=[0:31];
x1=n1.*(stepseq(0,0,31)-stepseq(16,0,31));
y1=stepseq(0,0,N1)-stepseq(8,0,N1);
X1=fft(x1,16);
Y1=fft(y1);
Z1=X1.*Y1;
z1=ifft(Z1);
subplot(3,1,2);
n1=[0:15];
stem(n1,z1);
title('16点圆周卷积');
axis([0,20,0,100]);
N=31;
n=[0:N];
x=n.*(stepseq(0,0,N)-stepseq(16,0,N));
y=stepseq(0,0,31)-stepseq(8,0,31);
X=fft(x);
Y=fft(y);
Z=X.*Y;
z=ifft(Z);
subplot(3,1,3);
stem(n,z);
title('32点圆周卷积');
axis([0,25,0,100]);