MATLAB使用fft求取给定音频信号的频率

一段10s立体声音频，采样率位8000Hz，已知频率为1000Hz

 clc;
 clear;
 [data, Fs] = audioread('1khz_stereo_8000.wav');
 fs=Fs;
 N=;
 n=*fs:*fs+N-;     %从第2s结束开始取样
 f=(n-*fs)*fs/N;     %数字频率转换成模拟频率 本来应该是n/N * fs，但由于n是从第三个周期开始的
 temp=data(:,);      %取单个声道样本
 x=temp(n);
 y=fft(x,N);
 mag=abs(y);          %求得Fourier变换后的振幅
 plot(f,mag);
 [max,pos] = max(mag);%找到幅值最高点，记录所在横坐标
 freq = f(pos);       %幅值最高对应基频分量
 disp(freq);
 grid on;

频谱图为