dft.m

连续信号离散傅立叶变换

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