用Mathematica研究自然对数的底数e.docVIP

用Mathematica研究自然对数的底数e 作 者：陈 龙 摘要：是一个奇妙有趣的无理数，它取自瑞士数学家欧拉的英文字头。与被认为是数学中最重要的两个超越数，、及（为虚数单位）三者间存在的关系。本文利用Mathematica软件研究了自然对数的底数，介绍了的一些相关知识、与自然对数的关系以及的值的计算方法等。 关键词：Mathematica，，自然对数 一、引言 远在公元前，圆周率就被定义为“周长与直径之比”。自古以来，的近似值一直取为3.14或。通过许多数学家的努力，的近似值位数不断增加。目前用电脑计算圆周率。由于电脑速度等功能不断改进，今后的近似值位数会越来越多。 另外一个奇妙有趣的无理数是，它取自瑞士数学家欧拉（Euler，1707-1783）的英文字头。欧拉首先发现此数并称之为自然数。但是，这种所谓的自然数与常见正整数1，2，3，……截然不同。确切地讲，应称为“自然对数的底数”。 与被认为是数学中最重要的两个超越数（transcendental number，若一数为之根，其中为某一至少一次的整系数多项式，则此数称为代数数（algebraic number），否则称为超越数）。、及（为虚数单位）三者间存在的关系。本文主要介绍的一些知识以及用Mathematica软件来计算。 二、欧拉数 考虑数列，==，，其中=，，，应用下述关于级数收敛的基本定理之一可证明出其极限存在。 定理1.设数列为单调且有界，则当时，（为一有限数）。 首先，对=，显然为单调递增数列。其次，=2，=，而时， =1+1+ 1+1+ = 1+3， 即数列以3为一上界。故有定理1知，数列收敛至一实数，由于此极限值与圆周率一样在许多数学的公式中出现，所以不可避免的需要给它一个特别的符号。欧拉似乎是第一个体会到此数之重要性的数学家，他并以来表示此数。后来符号就被广为采用，后人并称为欧拉数（Euler’s number）以纪念他。由于为时之极限，故可表示为 （1） = 。 以下说明如何以来求之近似值，事实上收敛至的速度极快。这里借助一几何级数，对任意， = + + + = + 故对， （2） + 。 若令，则上式为 （3） + 。 即对，与之差最多为。由于随着增长速度极快，故为的一个很好的估计值。例如，若=10，则与之差小于，因此经由计算，得到=2.718281…。 N[a[10],50] 2.7182818011463844797178130511463844797178130511464 N[E,50] 2.7182818284590452353602874713526624977572470937000 N[a[10]-E,50] True 当然若取大一些便可再更精确些，如=2.71828182845904523536028…。这是欧拉用笔算得到的之小数前23位。欧拉22岁时，在一篇论文中写着“这个数的对数是1，以命名之，它的值为2.71828…，它的常用对数为0.4342944…”。 是无理数的证明（这是欧拉在1737年所证出），可利用前述（3）对的估计式。设=为一有理数，其中，为二互质正整数。易见，此因介于2与3之间，故不可能为整数。现由（3）式知 + 。 将上式每项各乘以得 ++。 而由之定义知，为一整数，如此则得整数介于两相邻整数及+之间的矛盾结果。故不是有理数。 下面我们来看另一种常见的引进的方法。考虑数列 = ，。 则由二项式定理（Binomial Theorem）可得 = = = 1+1+ = 3 。 又由上面第三个等号的右侧可看出，的每一项对递增，且比多一正的项，故为一单调递增且有界数列必有极限。故得证存在。 接着证明。对，仍由前述第三个等号之右侧可得 。 若先固定，而令，则上式左侧趋近于，而右侧趋近于。即此时有，而又有，因此 ， 。 令，由夹逼定理，便得。也就是我们得到下述重要的极限结果： （4）