跳到主要内容

信号分析

突变点检测

在未知时刻搜索一条定量轨迹的均值或方差换挡位置,并用各段占用长度判断新体制是否站稳。

更新日期 2026-08-29 信号分析

这篇能回答什么

突变点检测回答的是:「这条按时间排好的定量轨迹,均值或方差在尚未指定的某个时刻有没有换挡」。算法在整条记录上搜索代价下降最大的切分,而不是先听你报一个日期再做 F 检验。教学结果把采样点折进分段占用:前段(类别 A)、中段(类别 B)、后段(类别 C)各有多少个点。占用是段长,不是「检出了几个变点」。

适合:工艺指标突然抬升、质量均值换班、宏观序列在未知年份改平台。教学示例用时间序列表里已排序的收益或产值。已经知道新政从某月实施、并且有处理组 / 对照组,主识别应是双重差分或事件研究,不要让搜索算法替你「发现」那个已知日期。

什么时候不该用

断点日期事先写在研究设计里时,用回归里的 Chow 分割或事件窗,而不是本页的盲搜。水平缓慢爬坡、没有台阶,先ADF或差分,搜索会在斜坡上切出假台阶。只要「谁领先谁」,去格兰杰因果。单列外推仍用 Auto ARIMA。

如何选择方法

先问断点是已知日历还是未知位置。已知政策日、要点估计处理效应:Chow(同一套回归系数在切点两侧是否相同)或 DID / 事件研究。未知何时换均值、只想把一条波切成几段:本页。Chow 的原假设写在预先指定的 τ;本页的 τ 是搜索出来的,p 值不能按「只检了一次」来读。

  • 突变点检测:1 条有序定量列;输出切分位置,以及前 / 中 / 后段占用。

  • Chow 分割:回归里预先指定一个(或少数)日期,检系数是否跳跃。

  • 双重差分 / 事件研究:已知干预时点,并且有对照轨迹。

  • I-MR / EWMA 控制图:在线报警「这一刻是否出界」,不负责回溯切历史。

输入、输出与限制

输入

分析对象是一条按时间升序的定量列。行序就是搜索轴。不要把截面户或乱序编号当成时间。

  • 分析序列:最少 1 个,最多 1 个。

输出

占用表把采样点打进三段互斥体制:类别 A 前段、类别 B 中段、类别 C 后段,给出频数、占比与合计。图在轨迹上标出切分。读「何时换挡」看图上的竖线;占用只说明各段有多长。

产品限制

  • 占用合计是分段点数,可以接近行数,但缺失行不进入任何一段。

  • 最少段长有下界,贴边的「半个点」不会单独成段。

  • 类别 A/B/C 是按时间顺序的段,不是「均值低 / 中 / 高」三档。

  • 一次运行默认按均值(或方差)损失搜索,不自动拟合带协变量的回归。

  • 检出多个切点时,占用仍折成三档展示,不要把三档理解成「恰好两个变点」。

案例

研究问题。 一条产值或收益轨迹看起来后半段抬高了,但没有事先写进备忘录的换挡日。要先估计切在何处,以及新体制是否已经比旧体制更长,再决定要不要分段建模。

教学示例的分段占用合计 234 个点:前段(类别 A)54(23.08%),中段(类别 B)70(29.91%),后段(类别 C)110(47.01%)。后段最长,提示搜索把较新的一段判成单独体制。下载时间序列表后,可用「annual」产值或「returns」收益练习。234 与工作表行数不必相等。

下载教学示例表(时间序列)

在 StatsX 中操作

  1. 步骤 1 · 上传表格

    打开工作台,将 Excel / CSV / JSON 拖入上传区。练习台阶用「annual」产值;练习噪声里的弱切分用「returns」。单文件不超过 20 MB。

  2. 步骤 2 · 查看画像与质量

    目标列须为定量、已按时间升序。先看轨迹:只有平滑斜坡、没有台阶,搜索容易切在任意处。事件 0/1 列是已知日期的标记,不要当本页的 Y。

  3. 步骤 3 · 选择「突变点检测」

    在算法目录搜索「突变」「变点」或「changepoint」,点选突变点检测。它是回溯切分,不是控制图的在线规则。

  4. 步骤 4 · 指定变量并运行

    分析序列只勾一列。不要同时塞进处理组虚拟变量——那是 DID 的角色,不是搜索轴。

  5. 步骤 5 · 阅读结果

    先看三段长度是否合理(有没有一段短到像噪声切分),再在图上读竖线对应的期次。后段 47.01% 不是「变点个数 = 3」。

StatsX 工作台上传时间序列以便做突变点检测
将教学示例表拖入工作台上传区。
列画像确认有序定量轨迹
列画像给出类型与缺失,便于确认只切一条按时间排列的列。
在目录中选择突变点检测
在算法目录中选定「突变点检测」。
为突变点检测指定序列
分析序列勾选一列定量指标后运行。

如何阅读输出

占用表把 234 个采样点折进三段。前段 54(23.08%)最短,中段 70(29.91%),后段 110(47.01%)最长。这是「新体制已经持续更久」的长度证据,还不是「发生了三次跳跃」。没有逐点 τ̂ 列表时,不要把 47.01% 写成变点日历。

突变点检测结果表
项目类别统计量取值
突变点检测类别 A频数54
突变点检测类别 B频数70
突变点检测类别 C频数110
突变点检测合计有效样本 N234
突变点检测类别 A占比(%)23.08
突变点检测类别 B占比(%)29.91
突变点检测类别 C占比(%)47.01

左右滑动可查看完整表格

教学示例表在 StatsX 工作台中的计算结果口径。

突变点检测分段占用表
工作台中的前/中/后段占用,数字应与上表一致。
带切分标记的轨迹或分段占用柱
换挡时刻在轨迹竖线或附加位置表上读;占用表只提供段长。

前段 23.08%、中段 29.91%、后段 47.01%。后段领先只说明搜索把较新的点划进同一体制。若轨迹仍是斜坡,这三段可能是算法在单位根路径上的过切。先 ADF,再决定要不要相信竖线。Chow 不会出现在这张表里:它要求你先指定日期,再比较两侧回归。

解读顺序:序列是否已排序 → 目视有无台阶 → 最短段是否短到像假阳性 → 图上 τ̂ 的业务含义 → 才写「在第 t 期换挡」。分段之后的模型(各段各自 ARIMA,或加虚拟变量的回归)在本页之外完成。

理论与公式

均值变点把记录切成若干段,每段内部用自己的均值拟合,搜索使残差平方和(或似然)下降最多的位置。单变点可写成:

多变点常用 PELT 或二元分割,用惩罚项阻止切得过碎。方差变点则比较各段样本方差。这些 τ 都是数据依赖的。相比之下,Chow 检验预先固定 τ₀,在回归 y=Xβ+ε 上比较「全样本一套 β」和「τ₀ 两侧两套 β」:

把搜索出来的 τ̂ 再塞进 Chow,F 会被「已经挑过最大切分」抬高,名义 p 值偏小。本页占用表在切分之后清点各段点数。它是段长直方图,不是 F 统计量,也不是变点个数。后段占比高,只说明新体制覆盖了更多点。

注意事项

  • 已知政策日时不要用盲搜代替识别策略;搜索会「恰好」切在你已经知道的日子上,看起来像发现。

  • 单位根或强漂移会制造假台阶,先检验平稳性。

  • 三档加总是 100%,不能把占比解释成「发生了三次突变」。

  • 最短段若只剩十几个点,优先怀疑惩罚太小,而不是写「体制极不稳定」。

  • 控制图报警的是当前点,变点检测回溯整段历史,两者日期不必相同。

  • 多列联动换挡应在 VAR 或方程组里讨论,不要对每一列各切一次再强行对齐。

常见问题

参考文献

  • Chow, G. C. Tests of Equality Between Sets of Coefficients in Two Linear Regressions[J]. Econometrica, 1960, 28(3): 591–605.

  • Killick, R., Fearnhead, P., & Eckley, I. A. Optimal Detection of Changepoints With a Linear Computational Cost[J]. JASA, 2012, 107(500): 1590–1598.

  • Bai, J., & Perron, P. Estimating and Testing Linear Models with Multiple Structural Changes[J]. Econometrica, 2003, 71(1): 47–78.

  • 陈强. 高级计量经济学及 Stata 应用[M]. 高等教育出版社.