ms16_3.m

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

wc=3;
Tmax=0.1*pi/wc
Ts=input('(Ts<Tmax)Ts=');
Nmin=10*2*pi/wc/Ts
N=input('(N>Nmin)N=');
D=2*pi/(Ts*N);
M=floor(wc/D);
Xa=[ones(1,M+1),zeros(1,N-2*M-1),ones(1,M)];
n=-(N-1)/2:(N-1)/2;
xa=abs(fftshift(ifft(Xa/Ts,N)));
plot(n*Ts,xa,'k');