8.1t.m

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

N=12;w1=2*pi/N;k=0:N-1;
x=[1,2,3,4,5,6,6,5,4,3,2,1];
X=dft(x,N);
magX=abs(X),phaX=angle(X)*180/pi
subplot(2,1,1);
axis([-0.1,12.1,0,50]);
hold on
stem(k,magX);
ylabel('|x(k)|');title('DFT的幅度:N=12');
hold off
subplot(2,1,2);
axis([-0.1,12.1,-400,400]);
hold on
stem(k,phaX);
ylabel('度');title('DFT的相角:N=12');