常州大学数值分析作业—第三章.docVIP

第一章：9.设，给出计算函数值的一个合适算法，并在字长m给定的，十进制计算机上给出数值计算结果。 解：由 得 字长为5的十进制计算机上计算和，并与的精确值 1.0075376410479比较，说明差异存在理由，其中，。 解：字长为5时的误差很大，这是因为设置的字长有限，就不可避免的使舍入误差不断积累。 把字长改为9时，误差已经大幅度减小。这说明，加大字长可以显著减小误差。 举例介绍数组矩阵常见运算。 解：举例如下 12.对任意给定的实数a、b、c、试编写Matlab程序,求方程的根。 解：利用教材例11的方法： 当b0时，，。 当b0时，，。 利用及，给出一个计算π的方法，根据此方法编写程序，给出π的至少有10位有效数字的近似值。 解：根据题中所给公式，容易得到： 分别利用下式给出计算ln2的近似方法，编写相应程序并比较算法运行情况。 解： 由运行结果可知， 方法二的绝对误差比方法一的误差要小得多。 这是因为方法一给出的计算公式含有相近数相减项，损失了有效数字。 而方法二给出的计算公式避免了相近数相减，具有较好的精度。 分别用Jacobi迭代法、Gauss-seidel迭代法解方程组 解：Jacobi迭代法收敛，Gauss-seidel迭代法不收敛。 编写LU分解法、改进平方根法、追赶法的Matlab程序，并进行相关数值实验。 3.将矩阵进行Doolittle和Crout分解