dtft.asv

来自「数字信号处理（机械工业出版社）的源码」· ASV 代码 · 共 13 行

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);