function f=sijie_hanshu() clc clear close all
Fs=48000;%%%采样率 t=0:1/Fs:0.02;%%时长 f=1200; w=2pi40;
x=5sinc(wt).exp(1i2pif*t); %x信号 x=awgn(x,-5); % figure % plot(t,abs(x));
tt=-0.01:1/Fs:0.03; ff=200; x_delay=sinc(wtt).exp(1i2pifftt); %x信号 x_delay=awgn(x_delay,-5); % figure % plot(tt,abs(x_delay));
%%%时差为0.01 频差为1000 xor_len=length(x_delay)-length(x)+1; FFT_len=2^11;
xor_value=zeros(xor_len,FFT_len); %%行列关系 %%%滑动求解二维互模糊 for i=1:xor_len temp=x.conj(x_delay(i:i+length(x)-1)); xor_value(i,:)=fftshift(abs(fft(temp,FFT_len))); %%可以忽略exp（1i2pi(Fd-f)*n）？ end
[a b]=max(xor_value,[],2); %%每行最大值组成a， b [c d]=max(a);
b(d); %% delta_f=Fs*b(d)/FFT_len-Fs/2; %%频差 delta_t=d/Fs; %%时差
f=delta_f
