tic%計算運行時間
N=30;%時域采樣信號有N點
ft=600;%時域的采樣頻率
dt=1/ft;%時域采樣間隔dt
%采樣時長為N*dt
n=0:N-1;
t=n*dt;%時域采樣對應的時間
L=512;%DFT點數
ws=2*pi*ft;%頻域下的采樣周期
f1=100;
f2=120;
x=cos（2*pi*f1*t）+cos（2*pi*f2*t）;%時域信號
x（1L）=0;?%補零
a=zeros（1L）;%空數組以保存實部序列
b=zeros（1L）;%空數組以保存虛部序列
c=zeros（1L）;%空數組以保存幅值序列
for?k=0:L-1
????for?i=0:L-1
????????a（k+1）=a（k+1）+x（i+1）*cos（2*pi*k*i/L）;
????????b（k+1）=b（k+1）-x（i+1）*sin（2*pi*k*i/L