clear;
N=1000;
fs=5000;
n=0:N-1;
xn=1+0.5*sin(30*pi*n/fs)+2*sin(80*pi*n/fs);
subplot(2,1,1);%2行1列显示
plot(n,xn);
title('f(t)=1+0.5*sin(30*pi*t)+2*sin(80*pi*t)')
xlabel('时间(s)');
N=40;
fs=200;
n=0:N-1;
xn=1+0.5*sin(30*pi*n/fs)+2*sin(80*pi*n/fs);
k=0:N-1;
nk=n'*k;
WN=exp(-j*2*pi/N);
WNnk=WN.^nk;
Fn=xn*WNnk;
subplot(2,1,2);
amp=abs(Fn(1,1:N/2))*2/N;
f=n*fs/N;
plot(f(1,1:N/2),amp);
title('DFS')
xlabel('频率(Hz)');
ylabel('幅度');
grid on;
运行结果: