matlab-DFT函数

%%DFT函数，输入变量fk为时域离散点。Fs为采样频率
function H=DFT(fk,Fs)
N=length(fk); 
n=0:N-1;
k=0:N-1;
W=exp(-j*2*pi/N);
kn=k'*n;
fn = fk*(W.^kn);
subplot(211)
plot(n,fn);
subplot(212)
L = N;
P2 = abs(fn/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1) 
title('DFT')
xlabel('f (Hz)')
ylabel('|P1(f)|')            %内置DFT变换坐标
end