s2.m

来自「数字信号处理的MATLAB仿真加实验指导」· M 代码 · 共 9 行

%a12.m
N=8;n=0:N-1;
xn=sin(pi/8*n);
Xk=fft(xn,N);
magXk=abs(Xk);
k=(0:length(magXk)'-1)*N/length(magXk);
subplot(4,2,1);stem(n,xn);
subplot(4,2,2);stem(k,magXk);