ms14_1.m

《数字信号处理实验(MATLAB版)》这

n1=0;n2=7;n0=4;
n=n1:n2; 
N=length(n);
xn=[(n-n0)>=0];
subplot(2,2,1);stem(n,xn,'filled','k');
title('原序列x(n)');
k=0:N-1;
Xk=fft(xn,N);
subplot(2,1,2);stem(k,abs(Xk),'filled','k'); 
title('Xk=DFT(x(n))');
xn1=ifft(Xk,N);
subplot(2,2,2);stem(n,xn1,'filled','k'); 
title('x(n)=IDFT(Xk)');