example3_4.m

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

n=[0:1:3];
x=[1,1,1,1];
N=4;
X=dft(x,N);
magX=abs(X),phaX=angle(X)*180/pi
subplot(2,1,1);stem(n,magX);
xlabel('k');ylabel('|X(k)|');title('amplitude of DFT');
subplot(2,1,2);stem(n,phaX);
xlabel('k');ylabel('degree');title('phase of DFT');