每日精学 166｜单复正弦频率的集中似然

2026-09-25

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

单复正弦的未知复幅度可在给定频率下解析消去，将联合最大似然问题化为一维周期图峰值搜索。

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

01
模型与符号

独立圆对称复白高斯噪声；单个正弦频率和复幅度未知；只在一个不混叠频率区间内搜索。

c：未知复系数（复幅度）

ω：角频率（rad/sample）

J_c：消去复幅度后的残差（幅度平方）

02
从模型到公式

联合拟合

高斯负对数似然除去常数后等价于残差能量。

$J(c,\omega)=\sum_{n=0}^{N-1}|x_n-ce^{\mathrm j\omega n}|^2$

固定频率投影

复指数向量能量为N，最小二乘复系数由正交投影给出。

$\widehat c(\omega)=\frac1N\sum_nx_ne^{-\mathrm j\omega n}$

平方展开

残差分为数据总能量及与候选正弦的投影。

$J=\sum_n|x_n|^2-2\operatorname{Re}\!\left[c^*\sum_nx_ne^{-\mathrm j\omega n}\right]+N|c|^2$

集中目标

代入最优复系数，两个交叉项与平方项合并。

$J_c(\omega)=\sum_n|x_n|^2-\frac1N\left|\sum_nx_ne^{-\mathrm j\omega n}\right|^2$

频率解

最小残差等价于最大化连续频率周期图。

$\widehat\omega=\arg\max_{\omega\in[-\pi,\pi)}\left|\sum_nx_ne^{-\mathrm j\omega n}\right|^2$

离散搜索误差

有限网格只近似连续最大点。仅在局部光滑单峰且网格足够细时可用半网格量级估计；它不是一般严格界，更不是统计下界。

$\omega_q=\omega_{\min}+q\Delta\omega,\quad |\omega_{q^*}-\omega_*|\lesssim\Delta\omega/2$

集中似然作为一维子空间投影

固定试探频率 ω 后，复幅度的最小二乘解为

$\widehat c=\boldsymbol v^{\mathrm H}\boldsymbol y/N$

其中

$v_n=e^{\mathrm j\omega n}$

消去幅度后残差为

$J(\omega)=\|\boldsymbol y\|^2- \frac{|\boldsymbol v^{\mathrm H}\boldsymbol y|^2}{N}.$

因此最大化谱投影能量与最小化残差完全等价。

根据投影不等式，得分满足

$0\le\frac{|\boldsymbol v^{\mathrm H}\boldsymbol y|^2}{N}\le\|\boldsymbol y\|^2.$

该界可检测模板范数或归一化遗漏。

FFT 只提供离散网格上的候选，连续频率最大值通常不在栅格上；局部细化也不保证从错误主瓣跳到全局最优。多分量数据不满足单正弦模型时，最优投影仍可定义，但参数解释需谨慎。

03
物理含义

集中似然消除了不需要数值搜索的参数，只留下频率这一非线性参数。FFT可以快速给出一组候选值，但网格更密不等于记录更长：它减小数值搜索误差，不改变观测信息。

04
固定参数算例

用80点复正弦及噪声计算高密度候选频率，验证集中残差与直接重构残差一致，检查无噪声真频率残差为零及频率峰值位置。数值脚本使用固定随机种子；随机实验的判据采用有限容差，代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差，关键数值写入独立结果文件。

固定参数数值验证结果

量	数值
记录长度	80
true angular frequency	0.283
estimated angular frequency	0.2824
search step	0.0002
concentrated identity error	7.81597009336e-14
noiseless true residual	2.71479084961e-30
minimum noisy residual	8.80083340412

图 1　单复正弦频率的集中似然的固定参数数值验证

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

05
数字实现的边界

SystemVerilog实现候选复振荡器与输入的共轭乘积实虚块累加；候选频率生成和峰搜索由软件完成。端口为16位有符号码，内部64位，同步复位，valid空拍保持状态；

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

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

06
适用范围

用于单频仪器校准和参数估计基准。工程实现可粗网格定位后局部优化，但必须检查全局峰与局部峰，避免把优化器收敛等同于正确频率。

07
工程上需要注意

多音、调频及时间变化幅度违反单分量模型。

有限网格、窗口加权和频率区间限制都会改变数值最大点。

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

未知复幅度可解析消去；单复正弦频率最大似然因此等价于连续频率周期图峰值。

内容依据：每日精学第 166 课《单复正弦频率的集中似然》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读8

写评论...

每日精学 166｜单复正弦频率的集中似然

01模型与符号

02从模型到公式