functionFS=dft1(A,a,B,b,C)fs=10000;%采样点频率N=1000;%采样点个数n=0:(N-1);x=A*cos(a*pi*n)+B*sin(b*pi*n)+C;%定义一般性输入信号形式y=x;%定义一个数组s=0;%记录最大峰值FS=[0,0,0]

;%将返回值定义为数组用于返回多个数fork=1:N%实现离散傅里叶变换y(k)=0;n=1;while(n