无源定位
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