数值计算答案.doc

写在前面的话 各位同学：自《计算方法》开课以来，大家的学习热情很高，但也有一部分同学反映做习题有些困难。我们《计算方法》教学组的老师们商量以后，觉得可以考虑借助于网络，提供一部分辅导，通过演示和讲解部分练习题，使大家加深对教材内容的理解，学会如何用理论来解决有关题目和问题。 这只是一个尝试，效果如何还有待实践来检验。有一点是肯定的，那就学数学不能照葫芦画瓢，只有理解了的才能有解题的思路，而成功地解了题，反过来也会加深对理论的理解和认识，这是一个反复提高的过程。 对吗？ 在这里，我们还要感谢那些打印文稿的同学，正是由于他们的劳动，才使大家能看到这几章的内容。 数学系 《计算方法》教学组 2008.11.10 习题一 数值分析引论 提要： 了解数学问题与数值计算问题的联系与区别。 为什么要解数值计算问题？ 知道有哪些基本的数值计算问题。 求值。 方程求解 数值逼近 了解数值计算的基本数学思想和方法 基本思想： 1. 等价变换思想 2. 逐次迭代思想 3. 逐步迭代思想 4. 化整为零：数值积分，求微分方程数值解 5. 化曲为直：牛顿法 6. 转变问题类型 7. 外推思想 基本方法： 1. 直接方法 2. 间接方法：迭代，递归 误差分析 误差来源 绝对误差，相对误差，有效数字。 数值分析 6. 算法性态分析 误差分析中的重要概念与公式 1. 有效数字的位数n： 其中m为整数，为0-9中的一个数字， 如的绝对误差不超过末位的半个单位，即，则称具有n位有效数字。 因此，如的绝对误差界是某一位的半个单位，该位到的第一位非零数字共计有n位。则有n位有效数字。有效数字不仅给出了近似值的大小，还给出了它的绝对误差界。 2. 定理3.1 1）有效数字相对误差界：有n位有效数字相对误差界 2）相对误差界有效数字：的相对误差界满足，则它有n位有效数字。 3、 对 绝对误差： 相对误差： 对 绝对误差： 相对误差： 思考题1.1 设，（1）用左矩形公式，右矩形公式，梯形公式，simpson公式分别计算其近似值；（2）把原积分区间对分为，分别在两子区间用上述公式近似计算，再求和，与（1）的计算结果进行比较，你有何看法？ 解： 不能积分求出精确值。记， 左矩形公式 右矩形公式 梯形公式 Simpson公式 其中， 所以， 把 左矩形公式 右矩形公式 梯形公式 Simpson公式 2. 按定义计算行列式。一个N维行列式有项，每一项为数的乘积，总共要作次乘法和次加法，由公式stirling公式，，计算一个25阶行列式的值，计算量多大？（只计乘除法的次数），并假定计算机每秒可作次乘除法的计算，试求总共要计算多少时间？由此，计算解25阶线性方程组的Gremer法则要作多少时间？ 解：。 计算一个行列式的乘除法次数 。 需要时间： 即总共需要2.740亿年。 解一个25阶线性方程组： Gremer法则，要作26个行列式计算，故需要 用差分法离散化下列边值问题 对照教科书上的例1.6，写出离散后的线性方程组。 解： ， 。把区间作9等分，得： 记 故在内点处有 代入，得方程组： 此为三对角方程组。 思考题 1.2 涉及一个计算的保证收敛的不动点迭代算法。 解： 又 记 故有 迭代格式为 ， 设，有： ，，，，，， 2. 在外推法的介绍的算例中，我们采取了步长取半的方式，实际上只需要步长序列单调递减就可实现外推。试推出这样一般情况下的外推公式，并根据诸值，来计算圆周率的近似值 。 解：设圆直径为1，则其内接正n边形的周长为圆周长。 记，（即把圆周n等分），把n理解为步长。 又设 把原式Taylor展开： 故 其中， 选择当取时， （*） （**） 这里 用圆内接正边形的周长近似的误差，即用离散近似值 近似极限值的误差称为离散化误差。记 用