利用matlab软件求解常数e和欧拉常数γ.pdfVIP

利利⽤⽤matlab软软件件求求解解常常数数e和和欧欧拉拉常常数数γ

数学实验报告

利⽤matlab软件求解常数e和欧拉常数γ实验⽬的：

利⽤matlab软件计算常数e和γ，并尝试利⽤不同的算法计算，⽐较计算精和时间，找到较好的算法。

掌握matlab程序求和、求极限的⽅法，学会寻找更优算法。

实验内容：

1、求e

e可以来源于两个数列的极限和，即

en=lim(1+1/x)^x,(x-+∞)（1式）

sn=1/0!+1/1!+1/2!+1/3!+1/4!+1/5!+……（2式），

根据1式，可在matlab上设计如下代码：

forn=1:15

n=10^n;

e=(1+1/n)^n%求常数e的循环语句

end

formatlong%使结果显⽰16位双精数

结果：

e的标准值约为：2.71828182845904523536

由上述结果可知，使⽤1式，有很⼤的缺陷，不仅精连10^-7都没有，⽽且当n=10^9误差开始变⼤。

根据2式，可得如下代码：

sum=0;

t=1;

forn=1:18

t=n*t;

sum=sum+(1/t);

end%求常数e的循环语句

e=1+sum

formatlong%使结果显⽰16位双精数

结果：

e的标准值约为：2.71828182845904523536

如上所⽰，随着n的增⼤，e的计算值越来越接近e的真实值.但是，当n的值⼤于17后，计算的精不再提⾼，原因是双精型

数只能精确到16位，所以结果只有个位以及⼩数点后15位（最后⼀位是近似取的），⽽1/18!=1.56*10^-16，所以n超过18再

往下计算不会更精确。。

在1式代码中，(1+1/n)和n都只能精确到⼩数点后16位，两者相乘，结果精将只能精确到8位。在2式中，,每⼀项都能精确到

⼩数点后16位，⽽e是所有项的和，求和后仍然能够精确到⼩数点后16位。所以，对于某些使⽤数学软件求解的问题，如果对

精有要求，应该尽量使⽤加、减运算，少⽤其他的运算(例如乘、除、乘⽅、对数等)，这样可以提⾼运算精。

2、求γ

如此欧拉常数γ也可以使⽤matlab求出较为精确的值。

可由公式γ=lim(n→∞)[(1+1/2+1/3+…+1/n)-ln(n)]得出。

forn=1:10

s=0;

fori=1:10^n

s=s+(1/i);

end

y=s-log(10^n)%求常数γ的循环语句

end

formatlong%使结果显⽰16位双精数

结果：

γ的权威数值约等于0.57721566490153286060651209。

上述结果精约有10^-8，虽然精还有提⾼空间，然⽽matlab上运⾏时，结果表明，增加n值，精提⾼，运算时间也将⼤

⼤加长（通常n超过10就需要好⼏分钟）

使⽤级数来计算：

forn=1:9

s=0;

fori=1:10^n

s=s+(1/i)-log(1+1/i);

end

y=s%求常数γ的循环语句

end

formatlong%使结果显⽰16位双精数

结果：

γ的权威数值约等于0.57721566490153286060651209。

上述结果精也约有10^-8，n再⼤就计算时间⼤⼤增加了，且可能会出错。

上⾯两个算法都不能算出更为精确的欧拉常数的值，因为调和级数收敛较慢，因此matlab计算很缓慢。这也表明，这个算法还

有待改进，需要更好的算法才能计算出γ更精确的值。

实验总结：

本⽂探索了使⽤数学软件求解常见常数的⽅法，并⽐较了算法的优劣，⽽算法不同，计算精、时间相差很⼤，这表明，数学

不能完全依靠计算机，⼈所编写的算法也⾮常重要，计算机只是⼈的⼯具，⼈的思维能⼒才是最重要的。另外，我们也需要计

算机强⼤计算能⼒的帮助，学会使⽤数学软件，才能更加有效地发展数学。