数学实验三new.docVIP

尚国荣 数学实验三 最佳分数逼近 实验环境:mathematica 实验目的: 1、用“连分数展开”的方法计算圆周率的近似值； 2、通过实验来体会“连分数展开”的方法与其他方法的区别，比较各种方法的优劣； 3、尝试用“连分数展开”的方法对其他的数进行展开。 实验方法： 1、Mathematica中常用的展开数与多项式的函数的使用； 2、计算圆周率“连分数展开”方法，并且利用特定的函数来展开其他数 实验步骤： 一．多项式的展开与化简 多项式是表达式的一种特殊的形式，所以多项式的运算与表达式的运算基本一样，表达式中的各种输出形式也可用于多项式的输出。Mathematica提供一组按不同形式表示代数式的函数。如： 1. 对进行分解，使用的函数为Factor： 代码： 结果： 2. 展开多项式与，使用的函数为Expand: 代码： 结果： 3. 连个多项式相除，总能写成一个多项式和一个有理式相加， Mathematic中提供两个函数PolynomialQuotient和PolynomialRemainder分别返回商式和余式： 代码： 结果： 4. 化简与，使用的函数为Pimplify: 代码： 结果： 二. 的连分数展开 的求解方法之前我们已经有许多种，但都比较繁琐而且误差较大，如何找到误差较小的的近似值求解方法，我们在所得整数3的基础上进行分析，有了整数3，则=3+，其中是3的误差，。只要能找到的最佳分数逼近值，再加3就得到的最佳分数近似值。从而我们使用一种方法“连分数展开“，其原理是 为了寻找与接近的分数，先找接近的整数，显然是7.于是 这是祖冲之的效率。 在此基础上，我们可以再用上述方法，要找到比误差更小的分数近似值，只需要找到比整数7更接近的分数来作为的近似值。由于，其中。先找的最佳整数近似值，显然是16.于是，从而，这就得到祖冲之的密度。 如果还要进一步提高精确度，就应当在考虑的整数近似值16的误差，取的整数近似值294，则可得到的更好的近似值。将这样的过程无限的进行下去，得到的的越来越精确的近似值。 对于 的连续分式的项有理数二次无理数 (循环连续分式）的连续分式的 n 项是很有规律的连续分式前 1000 项的几何平均数