15t.m

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

n=0:10;
x=10*(0.8.^n);
x1=fft(x);
k=0:10;
y1=x1.*(exp(8*j*pi*k/11));
y=ifft(y1);
subplot(2,2,1);
stem(n,x);
title('原序列x(n)');
xlabel('n');
axis([0,10,0,12]);
subplot(2,2,2);
stem(n,y);
title('移位序列y(n)');
axis([0,10,0,12]);
n=0:10;
y=10*(0.8.^n);
x=[y(1:1:11) zeros(1,4)];
n1=0:14;
subplot(2,2,3);
stem(n1,x);
title('15点序列x(n)');
xlabel('n');
axis([0,14,0,12]);
x1=fft(x);
k=0:14;
y1=x1.*(exp(8*j*pi*k/15));
y=ifft(y1);
subplot(2,2,4);
stem(n1,y);
title('15点移位序列y(n)');
axis([0,14,0,12]);