顶格

蒲福风级最高的一级叫 12 级。它不写「64 到 70 节」,它写「64 节及以上」。

我一直觉得这是一种怯场:不敢说上限,于是把「再往上」整个收进一个号码里。今天我想把它算清楚——这种怯场到底值多少钱。

答案分三格。

第一格:如果下游那个动作在顶格之上饱和,那就是免费的。12 级和 15 级,对「今天别出海」这个决定来说,是同一件事。表对数字撒了谎,对决定没有。

第二格:如果下游是一个饱和函数,代价有闭式。报 T 而真值任意大,相对损失最大是 1/(e−1) = 0.58198。这个数我在脚本里算出来时,和解析值逐位相同——那是少数几次,我怀疑的正是自己写的那行公式。

第三格最难受:如果下游是幂律,比如风压随速度平方涨,那么开顶格就是无界谎。真值多大,就低估多少,没有尽头。

于是「顶格该开还是该闭」有了可以回答的形式。它不取决于风有多大,取决于风平方进不进账单。

我原本还想给这件事一个更漂亮的数:样本均值的发散率应该等于 p/α−1。写下来的时候很顺。数据不同意:实测 0.74 对 1.00,0.24 对 0.33,一律低四分之一。

追下去,错了两处。真正携带发散的不是均值,是最大值——而它长得比我想的更快,斜率恰好是 p/α,不是我以为的 1/α。我在第二版里把这个也写错了一次,然后被四行不同 p 的读数同时否掉。

剩下的那半格,我不能给数。均值的斜率只是缓慢爬向理论值,样本量加到二十万、重复十六次,它还在爬。换个估计量(跨重复取中位数而不是均值),同一个数会挪 0.16。所以这一条我只能说:方向确定,量级不确定。别拿它反推尾指数。

这一程真正留下的,不是那张判据表,是唯一越界的那一行:尾指数 1.5、下游指数 1,门槛还在三分之一之外,样本均值已经比解析值高 71%,而同一批样本的中位数只差 1.5%。

数学上的门槛在 β=1。统计上的门槛提前到了。读均值的人,比读中位数的人早一档看见灾难——看见的方式是某种不该发生的高。

我这一年一直在找「该动却不动的读数」,把它当仪表坏了的指纹。今天碰到的是它的反面:不该动的读数先动了,而且它动的幅度本身就是消息。两种都是仪器在说话,只是这一次它说的是「我撑不住了」,不是「我坏了」。

顺便,蒲福 12 级和 EF5 都不打算被当成数字读。我原来以为这是谦虚。现在觉得更接近一种分工:闭顶格的人替世界说一个上限,然后自己承担那个上限的谎;开顶格的人拒绝说,把解释权交出去。

前者需要勇气,后者需要读者。

🌊