两信号互相关时延计算峰值定位错误问题修正

问题现状 #

• 基于两信号互相关运算计算信号时延，始终无法得到正确计算结果
• 涉及输入数据：i1（1×2097152 复双精度型序列）、i2（1×2097152 复双精度型序列）、ti1（1×2097152 双精度型时间序列）
• 预期目标：互相关结果的标记点准确落在最大幅值点位置
• 当前异常：代码输出的标记点位置与实际幅值极值点存在明显偏差

现有实现代码 #

t = ti1; % 1x2097152 double
t_Step = t(2)-t(1);
t_tmp = [-flip(t(2:end)),t];
x = i1; %1x2097152 complex double
y = i2; % 1x2097152 complex double
corrL=length(x)+length(y)-1;

X=fft(x,corrL); % x的FFT计算
Y=fft(y,corrL); % y的FFT计算

Z=X.*(conj(Y)); % 频域哈达玛积
z=fftshift(ifft(Z)); % 转换为时域互相关结果

Ly=length(z);
tz= t_tmp;

[~,idz] = max(z) % 查找互相关结果最大值索引
TD = tz(idz)*t_Step      % 计算时延值

if(TD<0)
distance = -TD * 204357.1915 % 特定介质中光速为204357.1915，计算传输距离
else(TD<0)
distance = TD * 204357.1915
end
%绘图部分
%绘制原始信号与延迟信号
plot(t,x,'r',t,y,'g'); 
legend('Original','Delayed');
%绘制互相关信号
[z_max, index] = max(imag(z))
figure
%plot(tz,z,tz(index),y(index),'o'); hold on
plot(tz,z,tz(index),z_max,'o')

xlabel('Time');
ylabel('Magnitude');
title('z signal');
grid

异常效果示意 #

图中标记点未对齐互相关曲线的实际峰值位置，存在明显偏移

修正方案 #

代码存在3处核心问题导致定位偏差：

• 峰值检索对象错误：z是复数互相关序列，直接调用max(z)会按实部大小排序，绘图时调用max(imag(z))仅检索虚部最大值，二者都无法定位真实互相关峰值。互相关峰值对应复数序列的模最大值，需统一用abs(z)作为检索对象。
• 条件分支逻辑错误：else分支后多余编写了(TD<0)判断，属于无效语法，距离计算可直接取时延绝对值完成，无需额外分支判断。
• 时间轴匹配校验：互相关结果长度为length(x)+length(y)-1，需保证构造的时间轴tz长度与z完全一致，避免索引错位。

修正后核心代码段：

% 检索互相关模值的峰值索引
[~,idz] = max(abs(z));
TD = tz(idz)*t_Step;
% 直接取时延绝对值计算传输距离
distance = abs(TD) * 204357.1915;
% 绘图部分同步修改峰值检索逻辑
[z_max, index] = max(abs(z));

内容的提问来源于stack exchange，提问作者ASB