本福特定律和齐夫定律是一回事吗

Posted dog250

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了本福特定律和齐夫定律是一回事吗相关的知识,希望对你有一定的参考价值。

关于本福特定律的简单解释和推导,参见:
https://zhuanlan.zhihu.com/p/440462854

思考本福特定律,与齐夫定律对照,它们之间似乎可以相互推导,是真的吗?

本福特定律说首数为 n n n的概率: P ( n ) = log ⁡ 10 n + 1 n = log ⁡ 10 ( n + 1 ) − log ⁡ 10 n P(n)=\\log_10\\dfracn+1n=\\log_10(n+1)-\\log_10n P(n)=log10nn+1=log10(n+1)log10n

写成连续的形式: P ( x ) = log ⁡ 10 ( x + 1 ) − log ⁡ 10 x P(x)=\\log_10(x+1)-\\log_10x P(x)=log10(x+1)log10x

从这个形式上看,它是一个定积分 ∫ x x + 1 log ⁡ 10 e n d n \\displaystyle\\int_x^x+1\\dfrac\\log_10endn xx+1nlog10edn 。设不定积分式为 F ( x ) F(x) F(x),则:

F ( x ) = ∫ log ⁡ 10 e x d x F(x)=\\displaystyle\\int\\dfrac\\log_10exdx F(x)=xlog10edx

积分 F ( x ) F(x) F(x)实际上就是所有首数字概率的积累分布函数,其概率密度函数为一个反比例函数:

f ( x ) = log ⁡ 10 e x f(x)=\\dfrac\\log_10ex f(x)=xlog10e

从本福特定律的概念上讲,首数字为 n n n的概率可以写成两种形式:

  • 定积分的形式: P i n t ( n ) = F ( x ) ∣ n n + 1 P_int(n)=F(x)|_n^n+1 Pint(n)=F(x)nn+1
  • 概率密度的形式: P p r o b ( n ) = f ( n ) P_prob(n)=f(n) Pprob(n)=f(n)

连续化是为了拟合微积分计算,回到离散的形式:

P p r o b ( n ) = f ( n ) = log ⁡ 10 e n P_prob(n)=f(n)=\\dfrac\\log_10en Pprob(n)=f(n)=nlog10e

换一种写法:

P p r o b ( n ) × n = log ⁡ 10 e P_prob(n)\\times n=\\log_10e Pprob(n)×n=log10e

这看起来符合齐夫定律。来看下是不是。

经过了连续~离散变换,连续情况下的反比例形式不能用于离散情况的计算,只能直观理解 P ( n ) × n = 常 数 C P(n)\\times n=常数C P(n)×n=C。现在直接从本福特定律的结论入手,实际计算一下:

P ( n ) × n = n × log ⁡ 10 n + 1 n = log ⁡ 10 ( n + 1 n ) n P(n)\\times n=n\\times \\log_10\\dfracn+1n=\\log_10(\\dfracn+1n)^n P(n)×n=n×log10nn+1=log10(nn+1)n

设:

g ( n ) = ( n + 1 n ) n g(n)=(\\dfracn+1n)^n g(n)=(nn+1)n

g ( n ) g(n) g(n)快速逼近 e e e,但仅在 n n n取1~9时, g ( x ) g(x) g(x)有意义,分别为:

g ( 1 ) = 2 , g ( 2 ) = 2.25 , g ( 3 ) = 2.37 , g ( 4 ) = 2.44 , g ( 5 ) = 2.48 , g ( 6 ) = 2.52 , g ( 7 ) = 2.54 , g ( 8 ) = 2.56 , g ( 9 ) = 2.58 g(1)=2, g(2)=2.25,g(3)=2.37,g(4)=2.44,g(5)=2.48,g(6)=2.52,g(7)=2.54,g(8)=2.56,g(9)=2.58 g(1)=2,g(2)=2.25,g(3)=2.37,g(4)=2.44,g(5)=2.48,g(6)=2.52,g(7)=2.54,g(8)=2.56,g(9)=2.58

log ⁡ 10 x \\log_10x log10x单调递增,计算 log ⁡ 10 g ( 1 ) \\log_10g(1) log10g(1) log ⁡ 10 g ( 9 ) \\log_10g(9) log10g(9)的值,分别为:

log ⁡ 10 g ( 1 ) = 0.301 \\log_10g(1)=0.301 log10g(1)=0.301
log ⁡ 10 g ( 9 ) = 0.411 \\log_10g(9)=0.411 log10g(9)=0.411

它们相差非常小,可近似为符合齐夫定律。

这是为什么?

通过上述推导, P i n t P_int Pint P p r o b P_prob Pprob是可以相互转换的,只要可以将事情抽象成 P i n

以上是关于本福特定律和齐夫定律是一回事吗的主要内容,如果未能解决你的问题,请参考以下文章

本福特定律

齐普夫定律的起源

本福特定律(Benford‘s law)的直观解释

概率论杂记

墨菲定律的典故及详细内容

定律-墨菲定律:百科