%x=sin(w.*t)+0.3.*sin(3.*w.*t)+0.1.*sin(5.*w.*t)+0.08.*sin(7.*w.*t);
tic%开始计时
%取样
fs=256;%采样频率
N=256;%采样点数
T=N./fs;
dett=1./fs;
t=0:dett:T-dett;
f=10;%频率 hz
w=2.*pi.*f;
X=sin(w.*t)+0.3.*sin(3.*w.*t)+0.1.*sin(5.*w.*t)+0.08.*sin(7.*w.*t);
subplot(2,1,1)
stem(t,X)
title('取样点');
Xn=ones(1,N);
Z=0;
for k=0:N-1
for n=0:N-1
Y=X(n+1).*exp(-j.*2.*pi./N).^(n.*k);
Z=Z+Y;
Xn(k+1)=Z;
end
Z=0;
end
%Xn=fft(X);
Xn=Xn./N.*2;
Xn(1)=Xn(1)./2;
Xn(N)=Xn(N)./2;
subplot(2,1,2)
stem(t,abs(Xn))
title('离散傅里叶变换')
toc%结束计时