add3_1b.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 14 行

n=0:10;x=10*(0.8).^n;
[xep,xop]=circevod(x);
X=dft(x,11);
Xep=dft(xep,11);
Xop=dft(xop,11);
subplot(2,2,1);stem(n,real(X));axis([-0.5,10.5,-5,50]);
title('real(dft(x))');xlabel('k');
subplot(2,2,2);stem(n,imag(X));axis([-0.5,10.5,-20,20]);
title('imag(dft(x))');xlabel('k');
subplot(2,2,3);stem(n,real(Xep));axis([-0.5,10.5,-5,50]);
title('dft[xep(n)]');xlabel('k');
subplot(2,2,4);stem(n,imag(Xop));axis([-0.5,10.5,-20,20]);
title('dft(xop)');xlabel('k');