每日精学 151|互相关时延的符号约定

2026-09-25

← 技术专题 · 每日精学目录

雷达 · 通信 · 电子战 / 基础精学 151

互相关峰值的正负号取决于两个信号的排列及共轭位置。统一定义可避免数字时延读数在接口间反向。

实际记录只包含有限个采样值。由记录计算得到的数值与总体参数并不相同:前者受记录长度、噪声、边界约定及运算精度影响,后者属于给定模型的属性。建立可复核结论需要先固定输入含义、统计假设和归一化方式,再比较可计算量与理论量。

01
模型与符号

有限实序列在记录外补零;第二路是第一路的正向整数延迟;相关采用非循环定义。

d:第二路相对第一路的延迟(样本)

k:相关扫描滞后(样本)

Cyx[k]:第二路相对第一路的互相关(幅度平方)

02
从模型到公式

信号模型

正时延表示第二路在更晚的样本时刻出现相同内容。

y[n]=a x[n-d],\quad a>0

统一相关定义

将第二路向前取样,与第一路同一基准索引相乘。

C_{yx}[k]=\sum_n y[n+k]x^*[n]

代入延迟模型

以自相关表示互相关,可直接确定峰位。

C_{yx}[k]=a\sum_nx[n+k-d]x^*[n]=a C_{xx}[k-d]

峰值位置

当序列自相关具有唯一零滞后最大值时,峰值对应正延迟。

\widehat d=\arg\max_k|C_{yx}[k]|=d

交换输入

变换求和变量并取共轭,可得交换次序后的镜像关系。

C_{xy}[k]=C_{yx}^*[-k]

频域对应

采用正变换负指数约定时,互谱排列须与时域定义一致。

C_{yx}[k]=\operatorname{IDFT}\{Y[m]X^*[m]\}[k],\quad M\geq N_x+N_y-1

FFT数组索引与物理滞后的换算

上式只列出有效线性滞后,其余零填充索引不参与取峰。直接按 M/2 分割仅适用于另行限定了中心化滞后范围的循环读数,不能覆盖任意不等长记录。采用负指数正变换时,乘积 Y X* 的逆变换把负滞后置于循环数组末端。设变换长为 M,数组索引 r 对应的有符号滞后可写为

k=\begin{cases}r,&0\le r<N_y,\\r-M,&M-N_x<r<M.\end{cases}

但有限线性相关的实际支撑还受到两路记录长度约束,零填充产生的空白区不应解释为新的有效滞后。

一个最小检查是令参考路仅在零处为一,延迟路仅在 d>0 处为一,则

C_{yx}[k]=\delta[k-d],\qquad C_{xy}[k]=\delta[k+d].

该检查同时区分输入次序、共轭约定和循环索引重排。

把相关幅度取绝对值不会消除滞后符号差异;符号来自索引排列,而不是相关值本身为正或为负。

03
物理含义

延迟符号不是信号固有的标签,而是相对参考路和索引方向的约定。时域峰位置、频域相位斜率以及接口说明必须使用同一排列。FFT末尾的负滞后区需要重排,不能直接把原始数组索引解释为物理延迟。

04
固定参数算例

以非周期随机序列和17样本延迟验证两个输入次序的峰位分别为正17与负17,比较直接求和与足长FFT的全部相关样本。数值脚本使用固定随机种子;随机实验的判据采用有限容差,代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差,关键数值写入独立结果文件。

固定参数数值验证结果
量数值
template length80
true delay samples17
first order peak17
swapped peak-17
FFT 长度256
direct fft error2.84217094304e-14
峰值67.9058962166
图 1:互相关时延的符号约定的固定参数数值验证图 1:互相关时延的符号约定的固定参数数值验证

图 1 互相关时延的符号约定的固定参数数值验证

表中数值和图中曲线由 MATLAB 实际运行得到,随机种子为 20261070;运行版本与完整结果另存为独立文件。

05
数字实现的边界

SystemVerilog实现外部对齐两路后的实互相关乘积与能量块累加;相关扫描地址生成不在子核中。端口为16位有符号码,内部64位,同步复位,valid空拍保持状态;

按外部首尾标志进行块累加并在块尾输出。自检覆盖多组数据、空拍、边界和复位,仅验证所述整数子核,不代表完整算法的硬件验证。

独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。

06
适用范围

用于双通道采样一致性检查、接口延迟标定与数字信号链单元测试。跨软件库迁移时,可用单脉冲和已知正延迟作为符号回归用例。

07
工程上需要注意

周期信号可能有多个等高相关峰,不能得到唯一绝对时延。

FFT长度不足会产生循环回绕,改变峰值及符号解释。

有限精度实现必须区分输入量化、内部舍入和最终结果格式,不能仅依据浮点曲线估计硬件误差。

在规定的互相关定义下,延迟路放在第一个输入位置产生正峰位,交换两路则翻转滞后并取共轭。

内容依据:每日精学第 151 课《互相关时延的符号约定》。图表沿用原课固定参数数值结果,属于数值验证,不代表设备实测。

阅读7
分享
写评论...