2010数值分析练习题选讲.pptVIP

例21. 试确定下面积分公式中的参数使其代数精确度尽量高. 解: 因此 所以该积分公式具有3次代数精确度 . . 若 为 以 为节点的线性插值，则 由引理5.1得 的余项 注:Simpson公式 导出梯形公式 的余项 例22 解 可得各节点的值如下表 0 1 0.125 00.25 00.375 00.5 00.625 00.75 00.875 0 1 0例23. 为简单起见,依次使用8阶复合梯形公式、4阶 复合Simpson公式和2阶复合Cotes公式 解: * * 例1. 解: 故 注： 不是唯一的. 例2. 解: 可见,经四舍五入取近似值,其绝对误差限将 不超过其末位数字的半个单位 5位有效数字 1位有效数字 5位有效数字 例3. 解： 例4. 解线性方程组 解: 这个方程组若用Gauss消去法计算会有小数作除数的现象,若采用换行的技巧,则可避免. 绝对值最大 不需换行 经过回代后可得 事实上,方程组的准确解为 用Doolittle法解方程组 解 例5. 用平方根法解线性方程组 方程组的系数矩阵A对称正定. 设 得 由 得 例6. 解 分析线性方程组的性态 系数矩阵 这是一个病态方程组. 例7. 例8. 用迭代法求方程的近似解,精确到小数点后6位. 解: 本题迭代函数有两种构造形式 因此采用迭代函数 d1 = 0.1000000 d2 = -0.0105171 d3 = 0.1156e-002 d4 = -0.1265e-003 d5 = 0.1390e-004 d6 = -0.1500e-005 d7 = 0.1000e-006 由于|d7| =0.1000e-0061e-6 因此原方程的解为 x7 = 0.090525. x1 = 0.1000000 x2 = 0.0894829 x3 = 0.0906391 x4 = 0.0905126 x5 = 0.0905265 x6 = 0.0905250 x7 = 0.0905251 例9. 为线性收敛. 证: 所以 证: 令 则 所以 故该迭代法至少是平方收敛的. 至少是平方收敛的. 例10. 例11. 用Jacobi迭代法求解方程组,误差不超过1e-4 解: 依此类推,得方程组满足精度的解为x12 迭代次数 为12次 x4 = 3.0241 1.9478 0.9205 d = 0.1573 x5 = 3.0003 1.9840 1.0010 d = 0.0914 x6 = 2.9938 2.0000 1.0038 d = 0.0175 x7 = 2.9990 2.0026 1.0031 d = 0.0059 x8 = 3.0002 2.0006 0.9998 d = 0.0040 x9 = 3.0003 1.9999 0.9997 d = 7.3612e-004 x10 = 3.0000 1.9999 0.9999 d = 2.8918e-004 x11 = 3.0000 2.0000 1.0000 d = 1.7669e-004 x12 = 3.0000 2.0000 1.0000 d = 3.0647e-005 例12. 用Gauss-Seidel迭代法求解例11. 解: 通过迭代,至第7步得到满足精度的解x7 x1 =2.5000 2.0909 1.2273 d =3.4825 x2=2.9773 2.0289 1.0041 d =0.5305 x3 =3.0098 1.9968 0.9959 d =0.0465 x4 =2.9998 1.9997 1.0002 d =0.0112 x5 =2.9998 2.0001 1.0001 d =3.9735e-004 x6 =3.0000 2.0000 1.0000 d =1.9555e-004 x7 =3.0000 2.0000 1.0000 d =1