根據有限長序列離散傅立葉變換公式計算:

f1=50; %10Hz

f2=100; %100Hz

%抽樣頻率

Fs=1000; %100Hz

%抽樣點數N

L=10;

N=2^L;

%抽樣脈衝序列

n = 0:N-1;

t = n./Fs;

% f2 一個週期的取樣數

M = floor(Fs/f2);

%被取樣訊號

x = cos(2*pi*f1.*t)+sin(2*pi*f2.*t);

%取樣序列

subplot(311);

stem(t(1:2*M),x(1:2*M));

hold off;

%傅立葉變換

%根據有限長序列的離散傅立葉變換公式計算DFT

n = 0:N-1;

k = 0:N-1;

F = x * exp(-j*2*pi/N).^(n"*k);

subplot(312);

plot(n,abs(F));

subplot(313);

plot(k,angle(F));