dtft.asv

数字信号处理（机械工业出版社）的源码

function [X,w] = dtft(x,n,k)
% Computes Discrete Time Fourier Transform
% -----------------------------------
% [X,w] = dtft(x,n,k)
% x = x(n) n point finite-duration sequence
% n = sample point row vector
% k = -M:M, (2M+1) point row vecor
% X = DTFT sequence for w frequency points
% w = frequency point row vector for [0,pi] duration
%

w=(pi/abs(max(k)))*k;
X=x*(exp(-j*pi/abs(max(k)))).^(n'*k);