dtft.m

来自「数字信号处理中的几个有用的变换的实现」· M 代码 · 共 24 行

function y=DTFT(x)
%
%x是序列幅度向量
%画出任意序列的傅立叶变换,模及相位
%频率取值于[0，2PI]
%可以取实部，可以取绝对值，也可以取相位
  % x=[1 1 1 1];
    m=size(x,2);
   %n=n1:n2
    n=0:m-1;
    k=-500:500;
    S=size(k);
    M=S(2);
    w=((4*pi)/M)*k;
    X=x*(exp(-j*4*pi/M)).^(n'*k);
%     X1=real(X);
     X2=abs(X);
    subplot(211)
    plot(w,X2);
    fa=angle(X);
    subplot(212)
    plot(w,fa);