dft.m

来自「连续信号离散傅立叶变换」· M 代码 · 共 7 行

function [Xk] = dft(xn,N)
% DFT
n = 0:1:N-1;
k = 0:1:N-1;
xn1 = [xn,zeros(1,N-length(xn))]; 
Xk = xn1*exp(-j*2*pi*n'*k/N);