每日精学 155｜互谱相位斜率时延估计

2026-09-25

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

整数相关峰只能给出离散时延栅格上的位置。宽带信号的互谱包含连续相位信息，时延可以由相位随频率变化的斜率求得。该方法把亚采样时延测量转化为直线拟合，同时引入相位展开、可靠频带选择及频率单位一致性的要求。

01
模型与符号

选定频带内互谱非零；相位可正确解缠；频点等权、误差近似等方差；未知常相位作为截距处理。

f_m：频率（Hz）

ϕ_m：互谱解缠相位（rad）

τ：相对时延（s）

ϕ₀：常相位（rad）

02
从模型到公式

相位模型

第二路相对第一路延迟时，互谱相位斜率为负。

$G_{yx}(f)=A(f)e^{\mathrm j\phi_0}e^{-\mathrm j2\pi f\tau},\quad\phi_m=\phi_0-2\pi f_m\tau+e_m$

最小二乘目标

常相位与时延共同进入线性拟合。

$J(a,b)=\sum_m(\phi_m-a-bf_m)^2$

消去截距

令截距导数为零，把频率和相位分别中心化。

$a=\bar\phi-b\bar f,\quad J(b)=\sum_m[(\phi_m-\bar\phi)-b(f_m-\bar f)]^2$

斜率解

对斜率求导，得到协方差与频率离差平方和之比。

$\widehat b=\frac{\sum_m(f_m-\bar f)(\phi_m-\bar\phi)}{\sum_m(f_m-\bar f)^2},\quad\widehat\tau=-\frac{\widehat b}{2\pi}$

估计方差

等方差不相关相位误差下，频带离差越大，时延方差越小。

$\operatorname{var}(\widehat\tau)=\frac{\sigma_\phi^2}{4\pi^2\sum_m(f_m-\bar f)^2}$

解缠条件

相邻频点纯时延相位跳变小于半周是简单顺序解缠的必要安全条件。

$2\pi|\tau|\Delta f<\pi\quad\Longrightarrow\quad|\tau|<\frac1{2\Delta f}$

频率平移不改变时延估计

把所有拟合频率改写为 f'_m=f_m−f_c，斜率分母保持不变：

$\sum_m(f'_m-\overline{f'})^2=\sum_m(f_m-\bar f)^2.$

相位与中心化频率的交叉和也不变，因此估计时延不依赖任意载频坐标原点。

改变的是截距，它吸收 −2π f_cτ。对于等间距 N 点网格，可直接求出

$\sum_m(f_m-\bar f)^2=(\Delta f)^2\frac{N(N^2-1)}{12}.$

该式把方差与频率间隔、点数联系起来，但增加高度相关的频点并不等价于增加独立观测。

频带加宽和网格加密是不同动作；若相位噪声协方差不再为对角阵，应以完整协方差计算斜率方差。

03
物理含义

常相位改变相位直线的高度，时延改变其斜率。中心化把这两种效应分开，也避免大频率偏置导致正规方程病态。有效频带越宽，同样的相位噪声对斜率的影响越小；孤立幅度深谷则可能触发整周解缠错误。

04
固定参数算例

用64个频点、已知3.2毫秒时延及固定高斯相位扰动，比较单次拟合、2000次重复方差和闭式方差，并验证无噪声拟合与加入常相位后的时延不变。数值脚本使用固定随机种子；

随机实验的判据采用有限容差，代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差，关键数值写入独立结果文件。

固定参数数值验证结果

量	数值
frequency count	64
true delay seconds	0.0032
single estimate seconds	0.00317254385107
empirical variance	2.03736428925e-10
theory variance	1.84131766427e-10
noiseless error	4.33680868994e-19
phase offset error	4.33680868994e-19

图 1　互谱相位斜率时延估计的固定参数数值验证

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

05
数字实现的边界

SystemVerilog实现频率与相位的乘积、频率平方和两路总和累加；中心化、除法和相位解缠不在子核内。端口为16位有符号码，内部64位，同步复位，valid空拍保持状态；

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

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

06
适用范围

用于双通道群时延标定、滤波链路延迟测量及线性相位一致性分析。频率单位必须在输入输出接口明确，避免将每样本弧度与赫兹混用。

07
工程上需要注意

多径或色散会产生非线性相位，单一时延不能完整解释。

频点噪声相关时，给定方差公式需改成广义最小二乘形式。

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

互谱相位斜率提供连续值时延估计；中心化拟合隔离常相位，但相位解缠与有效带宽决定其可靠性。

内容依据：每日精学第 155 课《互谱相位斜率时延估计》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读8

写评论...

每日精学 155｜互谱相位斜率时延估计

01模型与符号

02从模型到公式