圆周率的表达式的一种数学推导
题设:求圆周率的一种有理数表达式(近似值)
根据三角函数知识,有:
tan(π4)=1
那么有:
π=arctan(1)×4
至此,我们已经把圆周率表示为了一个连续可导的函数的值的形式。接下来尝试用有理数来表示出它的近似值。
对于连续可导函数,可以使用泰勒级数,将其展开为幂级数的形式,以求取近似值。
泰勒级数的表达式是这样的:
f(x)=∑n=0∞1n!f(n)(x0)(x−x0)n
关于泰勒级数的收敛性此处不再赘述,它在高等数学的教材中将是重要的一章。泰勒级数的思想是用(连续可导)函数在相邻点处的函数值以及各阶导数值来逼近所要求取的位置的值。
令:
f(x)=arctan(x)
则有:
f(1)(x)=11+x2
使用泰勒级数时,我们希望所建立的函数的各阶导函数呈现出比较简约的形式,而上方的函数的多次求导明显将是一件较为复杂的事情,因此进行换元。构造如下函数:
g(t)=11+t=(1+t)(−1)
则有:
g(n)(t)=(−1)n⋅n!⋅(t+1)−n−1
所构造的函数的高阶导函数有很简约的形式,而且由 0 处展开明显可以更好地简化其导数值。
即有:
g(n)(0)=(−1)n⋅n!⋅(0+1)−n−1=(−1)n⋅n!
从而得到了函数g(t)的泰勒级数展开:
g(t)=∑n=0∞(−1)n⋅tn
带回得:
f(1)(x)=g(x2)=∑n=0∞(−1)n⋅x2n
对该式求定积分得原函数:
f(x)−f(0)=∫x0∑n=0∞(−1)n⋅z2ndz
定积分和求和可交换顺序:
f(x)=∑n=0∞(−1)n⋅∫x0z2ndz
计算得:
f(x)=∑n=0∞(−1)n⋅x2n+12n+1
最后,将 1 代入,得:
f(1)=∑n=0∞(−1)n2n+1
即:
f(1)=11−13+15−17+19−111⋯
综上:
π=4×(11−13+15−17+19−111⋯)
这便是人们所熟知的圆周率的无穷级数表达式。当然,这个表达式的产生要远远早于微积分的诞生,膜拜伟大的往圣先贤们…
笔者使用Matlab计算到了式中 n=1000000 时的情形,得到 3.141592653,即小数点后第9位,验证了结果。
小结:
本文简单运用了微积分和泰勒级数的知识,得出了圆周率的表达式的一种推导方法,这是笔者大学时学习高等数学(非数学专业)时产生的一些联想,因本人水平有限,难免存在错漏,欢迎广大朋友批评指正~
本文仅参考了同济大学的《高等数学》第七版下册,灵感来源于笔者脑洞,这种入门级的数学知识应该不会牵扯到侵权吧…缺少相关前置知识的朋友可参阅该书。