ms14_4_1.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 10 行

Ts=0.5;C=[32,64,512];
for r=0:2;
    N=C(r+1);
    xn=[ones(1,11),zeros(1,N-11)];
    D=2*pi/(N*Ts);
    k=floor(-(N-1)/2:(N-1)/2);
    X=fftshift(fft(xn,N));
    subplot(3,2,2*r+1);plot(k*D,abs(X),'k');
    subplot(3,2,2*r+2);stairs(k*D,angle(X),'k'); 
end