ms16r3.m

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

wc=3.14;wm=2*wc;
Ts=0.1*pi/wm;
N=floor(10*2*pi/wm/Ts);
D=2*pi/(Ts*N);
M=floor(wc/D);
w=[-(N-1)/2:(N-1)/2]*D;
Xa=[zeros(1,M+1),1,zeros(1,N-2*M-1),1,zeros(1,M)];
n=0:N-1;
xa=abs(fftshift(ifft(Xa/Ts,N)));
plot(n*Ts,xa,'k');