for i=1:1:n-window
%对从当前点开始的window长度的数据进行傅里叶变换
y =fft(data_bit(i:i+window-1));
y =abs(y);%得到目标频率傅里叶变换结果中对应的index
index_impulse =round(f/fs*window);%考虑到声音通信过程中的频率偏移,我们取以目标频率为中心的5个频率采样点中最大的一个来代表目标频率的强度
impulse_fft_5k(i)=max(y(index_impulse-2:index_impulse+2));
end