每日精学 164｜复正弦频率的CRLB与时间中心化

2026-09-25

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

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

复正弦频率与未知初相位相互耦合。消去初相位后，频率信息由采样时刻的中心离差平方和决定。

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

01
模型与符号

单个非零复正弦，幅度已知、相位与角频率未知；独立圆对称复高斯噪声的每复样本功率为sigma平方；使用局部无偏界。

ω：角频率（rad/sample）

ϕ：未知初相（rad）

A：正弦幅度（幅度）

t_n：中心化时间索引（样本）

02
从模型到公式

均值模型

噪声方差按复样本总功率定义。

$\mu_n=Ae^{\mathrm j(\phi+\omega n)},\quad\mathrm E|v_n|^2=\sigma^2$

均值导数

相位导数与频率导数只相差时间索引。

$\frac{\partial\mu_n}{\partial\phi}=\mathrm j\mu_n,\quad\frac{\partial\mu_n}{\partial\omega}=\mathrm jn\mu_n$

信息矩阵

复高斯模型产生二倍实部系数。

$\boldsymbol I=\frac{2A^2}{\sigma^2}\begin{bmatrix}N&\sum n\\\sum n&\sum n^2\end{bmatrix}$

消去相位

Schur补扣除由未知相位吸收的信息。

$I_{\omega\mid\phi}=\frac{2A^2}{\sigma^2}\left[\sum n^2-\frac{(\sum n)^2}{N}\right]$

中心离差求和

均匀索引的二阶中心矩有闭式。

$t_n=n-\frac{N-1}{2},\quad\sum_nt_n^2=\frac{N(N^2-1)}{12}$

最终下界

若频率用每样本周期表示，还要除以四倍圆周率平方。

$\operatorname{var}(\widehat\omega)\geq\frac{6\sigma^2}{A^2N(N^2-1)},\quad\operatorname{var}(\widehat f)\geq\frac{\operatorname{CRLB}_\omega}{4\pi^2}$

未知相位导致时间原点信息被消去

复正弦频率导数与时间索引成正比，未知常相位导数则为常数方向。消去相位参数后，频率信息只依赖中心化时间，而不依赖任意选择的时间原点：

$J_{\omega\mid\phi}=\frac{2|A|^2}{\sigma_c^2} \sum_{n=0}^{N-1}(n-\bar n)^2.$

这里 σ_c² 是复噪声模平方的均值。

等间隔索引的平方和为

$\sum_{n=0}^{N-1}(n-\bar n)^2=\frac{N(N^2-1)}{12}.$

该结果解释长记录的信息增长与约三次幂的关系，但它是局部无偏估计界，不包含频率混叠、错误峰选择和低信噪比阈值效应。

把复噪声总方差误写为每个实维方差会产生二倍尺度错误。

03
物理含义

时间原点改变初相但不改变可获取的频率信息。把时间置于记录中心使信息矩阵的相位频率交叉项为零，清晰显示频率精度依赖整个观测跨度，而不仅是样本数。

04
固定参数算例

对长度8至128比较信息矩阵求逆与闭式下界，检查任意时间原点平移后的Schur补不变，并展示已知相位与未知相位的差异。数值脚本使用固定随机种子；随机实验的判据采用有限容差，代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差，关键数值写入独立结果文件。

固定参数数值验证结果

量	数值
信号幅度	1.2
复噪声功率	0.3
N64 angular crlb	4.76953601954e-06
N64 cycle crlb	1.20813758731e-07
matrix closed relative error	5.55111512313e-16
time shift information error	0
N64 known phase crlb	1.22055055618e-06

图 1　复正弦频率的CRLB与时间中心化的固定参数数值验证

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

05
数字实现的边界

SystemVerilog实现外部中心化时间索引的平方和、索引和及计数累加；信息矩阵求逆在软件中完成。端口为16位有符号码，内部64位，同步复位，valid空拍保持状态；

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

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

06
适用范围

用于单频参数精度预算、记录长度选择及频率估计器蒙特卡洛结果的局部基准。噪声功率是复样本总功率还是单实分量方差必须明确。

07
工程上需要注意

低信噪比的频率跳峰产生阈值效应，实际误差可能远离局部界。

多个频率分量时信息矩阵不能直接沿用单频结果。

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

未知初相位要求扣除时间均值项，最终频率信息与中心化时间的平方和成正比。

内容依据：每日精学第 164 课《复正弦频率的CRLB与时间中心化》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读7

写评论...

每日精学 164｜复正弦频率的CRLB与时间中心化

01模型与符号

02从模型到公式