每日精学 260｜序贯概率比检验的对数累积

2026-10-01

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

01
摘要与关键词

固定样本检验先收集规定数量的数据；序贯检验则随证据累积决定继续或停止。对数似然比的加法递推便于数字实现，但 Wald 门限公式是基于越界近似的错误率配置，需记录实际越界行为。

技术分析在明确的模型和符号约定下给出基本关系、推导与适用边界，并用固定参数数值算例和独立整数接口自检分别核对模型运算与数字子核。两种验证层次的输入尺度和覆盖范围分别记录，有限算例不外推为完整系统性能。

关键词：序贯概率比检验的对数累积；条件模型；数值验证；数字实现

02
问题背景

该问题位于由模型关系到有限精度数据通路的连接处。数学等式首先规定各量的语义；随后需要确认计算采用的输入、范围与尺度符合等式前提。公开基础信号处理中的观测证据、约束集合和数值表示具有不同作用，混淆这些作用会使输出虽落在合法范围，却不再对应所声明的模型。

03
假设与符号表

独立 Bernoulli 观测，在假设零下参数 p0=0.3，在假设一下 p1=0.6。错误预算 α=0.05、β=0.1；最大 500 个样本，超限记为未决。

所有对数未特别注明时采用自然对数。向量和矩阵的行列方向由公式明确指定；概率质量、条件密度和未归一化权重分别按定义使用。整数码只表示某个中间量，不默认具有物理单位。实际数值算例保留所列固定参数，不能据其有限取值范围推断一般输入的性能。

S_n：累计似然比（自然对数单位）

A,B：上与下比值门限（正数）

p₀,p₁：两假设参数（概率）

04
基础公式与完整推导

单样本增量

两 Bernoulli 密度相除再取自然对数，得到仅有两个可能值的增量。

$\ell(x)=x\ln(p_1/p_0)+(1-x)\ln[(1-p_1)/(1-p_0)].$

累计递推

独立观测的联合似然比为乘积，因此自然对数按样本相加。

$S_n=S_{n-1}+\ell(x_n),\quad S_0=0.$

停止区域

介于两门限之间时继续收集观测；跨界后输出并停止当前序列。

$S_n\ge\ln A\Rightarrow H_1,\quad S_n\le\ln B\Rightarrow H_0.$

Wald近似门限

忽略跨越门限时的超调得到近似配置。离散增量通常产生超调，所以实际错误率需单独评估。

$A\approx\frac{1-\beta}{\alpha},\qquad B\approx\frac\beta{1-\alpha}.$

漂移方向

期望增量与相应 Bernoulli KL 散度关联，使证据平均向正确边界移动。

$\mathbb E_0\ell=-D(p_0\Vert p_1)<0,\quad\mathbb E_1\ell=D(p_1\Vert p_0)>0.$

05
一致性与敏感性分析

超调的定义

上界跨越时的实际累计证据高于门限，近似门限推导省略了该随机差值。

$O_+=S_\tau-\ln A\ge0$

停止时间是首次离开继续区间的样本编号；无穷或超限状态必须在实际接口中单独记录。

$\tau=\inf\{n:S_n\notin(\ln B,\ln A)\}$

模型关系与数值误差的区分

推导中的恒等关系在其假设下精确成立；计算时观测模型误差、输入量化和中间舍入属于不同来源。设理想运算为 y=f(x)，输入量化扰动为 e_x，具体数值实现误差为 e_a，则输出可写为

$\widehat y=f(x+e_x)+e_a.$

若 f 在所讨论的连续参数点可微且扰动足够小，一阶项为 Df(x)e_x。

地址、候选掩码、分支标签和状态版本是离散量，不满足这种局部微分假设；改变一个控制位可能换掉整个合法集合或数据来源。因此连续误差界不能用于替代控制语义检查。

只有将同一输入送入独立的数学参考与被测实现，差值才隔离数值实现误差。若两个程序各自生成随机输入，输出不同可能来自不同样本；若参考重复使用相同错误的符号约定，逐点一致也不能证明定义正确。固定参数算例、结构不变量和边界事务提供互补证据。

06
物理含义

停止时间由观测证据决定，平均样本量不是固定采样数。误判、未决和停止样本数是不同统计量。

数学输出需要按定义解释。无量纲比值、自然对数度量、概率权重和整数地址分别回答不同问题；即使都存放在定点寄存器内，也不能互相替代。输出幅度变化可能来自公共尺度、观测变化或可行集合变化，只有回到模型定义才能确定其含义。

07
数值验证

每个假设各 10000 条序列，最大 500 步，记录错误率、未决数和平均停止样本数，绘制停止样本数的经验分布。

配套 MATLAB 语法脚本实际在 GNU Octave 10.3.0 中执行，完成断言、数据导出与图表生成。MATLAB R2026b 启动返回许可证到期错误，因此本次没有取得 MATLAB 原生运行结果。

该边界只限定运行环境，不将历史 MATLAB 日志作为本次数值证据。完整实际输出、运行版本及绘图数据与脚本一起保存。

固定参数算例的实际关键数值

结果字段	实际值
false alarm observed	0.0407
miss observed	0.0711
mean samples h0	12.5875
mean samples h1	14.4157
undecided sequences	0

图 1　序贯概率比检验的对数累积的固定参数数值关系

图中横轴为脚本中规定的扫描参数或样本索引，纵轴为实际计算值；多条曲线按脚本列顺序对应各计算序列。曲线连接用于显示离散扫描结果，不意味着中间所有参数已经仿真。确定性穷举仅覆盖声明的有限域；随机向量用于补充不同组合，不构成连续参数上的证明。

08
数字实现与验证边界

八位有符号累计码加增量，宽位结果与上下门限比较，输出两边界标志。

整数接口及有效周期

独立模块统一使用同步高有效复位和输入输出有效标志。输入总线分成八个十六位槽，各槽低八位承载本子核实际使用的码字；控制总线提供对应掩码、标签或地址。未使用槽不参与数学运算。

输出分成四个三十二位槽，负结果采用二补码；保护位用于保存中间精确整数结果，不代表算法必须使用三十二位物理精度。

被测时钟上升沿先判断复位；复位优先于有效输入并清除输出数据与输出有效位。非复位的有效输入在该边沿形成输出，输出有效位同步置一。无效输入周期仅清除输出有效位，数据保持上一次有效结果。调用者只能在输出有效时消费数据，不能把保持的旧值当作新的计算结果。

量化尺度与中间位宽

若整数码 q 对应实际量 x=q2^−F，小数位数 F 应在接口中固定。两个同尺度量相加保持尺度，乘积的尺度为 2^−2F。

模块中的宽位整数运算先保存完整结果；若外部需要重新缩放，必须明确舍入方向和饱和策略。

$x_1+x_2=(q_1+q_2)2^{-F},\qquad x_1x_2=(q_1q_2)2^{-2F}.$

八位有符号输入包含最小码字 −128，它的正幅值 128 不能用八位有符号数表示。

绝对值、符号取反、乘积及多项加和必须先扩展位宽。精确整数参考直接使用宽整数计算，输出编码再转换为对应二补码，这样可区分算法差异与位宽回绕。

实际自检覆盖

Icarus Verilog 实际编译并运行独立自检，确认 4448 个有效向量、262 个无效周期和 5 次同步复位全部通过。输入包括零、一、127、128、254 和 255 的码字边界，多级组合与固定种子的随机向量。

自检同时核对有效标志、无效周期数据保持、复位与有效同时出现以及中途复位后继续接收。

实际仿真事件摘录

事件	输入低三十二位	输出低三十二位
0	00000000	00000000
1	00010001	00000002
2	007f007f	000000fe
3	00800080	ffffff00

事件摘录以十六进制保存接口位型，完整输入和全部输出槽见独立追踪文件。低位摘录不承担全部数值结果的解释；完整测试通过状态来自对每个有效事务的全总线比较。参考向量由独立 Python 整数计算生成，不依赖被测 RTL 内的比较函数。

上述验证针对所声明的整数子核。它没有运行器件综合、静态时序、功耗分析或板级测量，也没有验证包含缓存、存储仲裁和完整算法控制器的系统。模块采用固定界限循环、组合算术与时钟寄存器构成，可作为可综合源代码；具体资源数和最高频率仍需目标工具与器件约束。

09
适用条件与工程应用

适用于两简单假设与独立同分布观测。未知参数、时变噪声和样本相关需要修改模型。

工程接入应先确认输入是否正是公式中的量，再确定尺度、位序与事务边界。若某个量已由上游近似或饱和处理，应把该处理列入误差来源；后续子核的整数精确性不能恢复上游丢失的信息。并行和流水线可以改变吞吐与延迟，但不得改变同一事务中的控制与数据配对关系。

10
限制与非理想因素

有限样本上限、门限超调和定点累积会改变实际错误率。没有声明最大样本数的数字实现还需防止计数和累加器溢出。

模型不匹配与有限精度并不具有相同修正方法。前者需要重新给出观测关系、约束或统计假设；后者需要增加保护位、选择合适缩放或改进数值运算。有限样本的误差为零仅说明当前输入与参考一致，不能自动推出未测试输入、完整链路或实际设备的误差为零。

11
配套脚本文件

独立文件及用途

文件	用途
matlab/lesson260_verify.m	固定参数数值计算与图表生成
matlab/numerical_results.txt	本次实际运行环境及关键数值
matlab/plot_data.csv	图表逐点数据
systemverilog/lesson260_core.sv	所声明的可综合整数子核
systemverilog/tb_lesson260_core.sv	有效周期、边界与复位自检
systemverilog/reference_vectors.txt	独立整数参考向量
systemverilog/trace_results.txt	实际事务追踪

12
结论

期望增量与相应 Bernoulli KL 散度关联，使证据平均向正确边界移动。停止时间由观测证据决定，平均样本量不是固定采样数。误判、未决和停止样本数是不同统计量。

固定参数数值复验和整数子核自检分别支持模型层关系与数字接口行为。结论的适用范围受所列假设、候选或状态集合、量化尺度以及实际运行环境限制；完整系统性能需在扩展模型与验证条件下另行评估。

13
参考资料

MathWorks. Constant False Alarm Rate Detection. 在线原始资料。

GNU Octave. Arithmetic and Linear Algebra. 公开运行环境文档。

内容依据：每日精学第 260 课《序贯概率比检验的对数累积》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读4

写评论...

每日精学 260｜序贯概率比检验的对数累积

01摘要与关键词

02问题背景

03假设与符号表

04基础公式与完整推导