time_delay.txt

来自「用多个正弦产生仿真信号,两个阵元接受到信号后,估计两信号的延迟」· 文本 代码 · 共 50 行

广义互相关时延估计程序


SN1=20; 
nn=4096; 
tt=1:nn;
fs=1000; 
S=3*sin(2*pi*17*tt/fs)+2.5*sin(2*pi*35*tt/fs)+1.5*sin(2*pi*51*tt/fs)+sin(2*pi*102*tt/fs); %采样频率为1000

Ps=S*S'/nn;
ps=diag(Ps);
refp=2*10.^(SN1/10);
tmp=sqrt(refp./ps);
SS=diag(tmp)*S;
%SS=SS+randn(size(SS)); 

for j=1:nn
S2(j+549)=1/dist(2)*SS(j);
end
%S2=S2+randn(size(S2));

for j=1:nn
S3(j+580)=1/dist(3)*SS(j); 
end%
%S3=S3+randn(size(S3));

n=200;
m=580-549;
SS2=[S2(tao(2)+1:tao(2)+1+n),zeros(1,m)];
SS3=[zeros(1,m),S3(tao(3)+1:tao(3)+1+n)]; 

SS2=[SS2,zeros(1,length(SS2))];
SS3=[zeros(1,length(SS3)),SS3];

fx2=fft(SS2);fx3=fft(SS3);
Pxy=fx2.*conj(fx3)./(abs(fx2).*abs(fx3)+eps); 

R=ifft(Pxy);

R=real(R);
c=max(abs(R));
a=find(abs(R)==c) 
m=(a-length(R)/2)*1/fs 
figure(3)
plot(R)



另外 4 处求的是两信号的互功率谱，请问互功率谱和功率谱密度的区别，如果去掉 4 处的分母，可以吗？