每日精学 171｜Burg反射系数与双向误差

2026-09-25

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

Burg一阶反射系数通过同时最小化前向与后向预测误差得到，分母对两端数据能量作对称处理。

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

01
模型与符号

实值有限记录；一阶前向残差为 x_n+kx_n−1，后向残差为 x_n−1+kx_n；共同索引上的总输入能量为正。

f_n,b_n：前向与后向输入误差（幅度）

k：反射系数（无量纲）

J(k)：双向误差目标（幅度平方）

02
从模型到公式

两方向残差

同一个实反射系数同时作用于前后两个方向。

$f_n^{+}=f_n+kb_n,\quad b_n^{+}=b_n+kf_n$

联合能量

在共同有效索引上相加，避免单边归一化。

$J(k)=\sum_n[(f_n+kb_n)^2+(b_n+kf_n)^2]$

目标展开

定义交叉乘积和与总能量后，目标为简单二次函数。

$C=\sum_nf_nb_n,\quad D=\sum_n(f_n^2+b_n^2),\quad J(k)=D(1+k^2)+4kC$

最小点

正分母保证严格凸性。

$\frac{\mathrm dJ}{\mathrm dk}=2Dk+4C=0,\quad\widehat k=-\frac{2C}{D}$

模长边界

逐样本平方非负可给出交叉积绝对值的界。

$2|C|\leq\sum_n(f_n^2+b_n^2)=D\quad\Longrightarrow\quad|\widehat k|\leq1$

一阶输入

初始前后误差就是相邻样本，预测反馈系数等于反射系数的负值。

$f_n=x_n,\quad b_n=x_{n-1},\quad\widehat a=-\widehat k$

双向能量给出的反射系数界

实数前后向误差序列 f_n,b_n 的 Burg 系数由交叉和与双向能量构成。逐点不等式 2|f_nb_n|≤ f_n²+b_n² 给出

$\left|\frac{-2\sum_n f_nb_n}{\sum_n(f_n^2+b_n^2)}\right|\le1.$

等号发生于两列误差以相同符号关系成比例的退化情况；

严格稳定需要严格小于一。若分母为零，两路残差均为零，系数不能通过直接除法确定。

$\sum_n(f_n^2+b_n^2)=0\quad\Longrightarrow\quad f_n=b_n=0\ \text{对全部样本成立}.$

分子、分母使用的有效样本区间必须相同，否则上述能量界不再保证成立。

边界裁剪、空拍计数和符号约定都应先固定，再比较不同实现的反射系数。

03
物理含义

双向目标不偏向某一端的能量，尤其适合短记录。反射系数的模长界来自数据上的确定性不等式，而不是渐近统计保证；若前后序列完全成比例，边界可能取等。

04
固定参数算例

对0.7系数AR序列求Burg一阶反射系数，扫描双向目标确认最小点，并检查解析梯度为零、目标下降和模长不超过一。数值脚本使用固定随机种子；随机实验的判据采用有限容差，代数关系的判据采用浮点残差。图中同时展示参数变化趋势与对应误差，关键数值写入独立结果文件。

固定参数数值验证结果

量	数值
样本数	2000
true AR coefficient	0.7
estimated reflection	-0.708800514913
estimated AR coefficient	0.708800514913
objective gradient	0
unpredicted total energy	8053.2261683
minimum total energy	4007.30007922

图 1　Burg反射系数与双向误差的固定参数数值验证

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

05
数字实现的边界

SystemVerilog实现双向误差交叉积、两路平方和及样本计数累加；反射系数除法由软件完成。端口为16位有符号码，内部64位，同步复位，valid空拍保持状态；

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

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

06
适用范围

用于短记录AR谱建模与格型预测器初始化。数据均值和趋势应先按模型要求处理，避免低频确定性成分被误解释为极点。

07
工程上需要注意

高阶Burg需要逐阶收缩的有效索引，不能直接重复一阶公式。

模长接近一时有限精度可能导致稳定性边界问题。

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

一阶Burg估计是双向平方误差的唯一二次最小点，其幅度界由交叉积与能量关系给出。

内容依据：每日精学第 171 课《Burg反射系数与双向误差》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读8

写评论...

每日精学 171｜Burg反射系数与双向误差

01模型与符号

02从模型到公式