地大北京数值分析.docVIP

对于积分（n=0,1,2…）为何用递推公式，计算不稳定，应如何改进？（计算过程保留六位有效数字） 解：利用分部积分 容易算出，取六位有效数字0.367879推出： …… 故绝对错误，现分析误差传播： 依此类推，误差急剧放大。 改进算法：，若任取为出发值，误差可变为原来的倍，计算结果如下： … 对很大的正数N，如何构造算法，使的计算过程稳定 解： 令，，则 对大的N值计算 解： 按展开式 得 构造两种递推公式计算（n=0,1,2…）并说明两种方法的优劣 解：，从而得出两种递推公式： 算法1：，取，有误差，到误差放大了倍，算法不好 算法2：，任取，即使有误差，到的时候，缩小到了原来的，算法很好 用平方根法解方程组 解： ， ，故A对称正定 ，，， ， 由Cholesky分解求得 由解得，再由解得 用分解法解方程组 解：设 由矩阵乘法得 ， ， ， ， ， 由得：，， 由得： 为什么在用Gauss消去法解线性方程组时，选主元法稳定，试举例说明。 解：主元是Gauss消去法消元与回代过程中的分母，分母过小会放大误差，故应选绝对值最大者为主元。现举例说明： 计算过程保留两位有效数字 解 增广矩阵为 回代， 与准确解，相差很远 现改为选主元的方法： 回代， 可得到比较接近准确的解。 8、设，若用Gauss-Seidel方法解方程时，是否收敛？ 解：对Gauss-Seidel方法迭代矩阵 ， ，故Gauss-Seidel法收敛。 9、设方程组，试考察用Jacobi法求解的收敛性 解： 其中 ， Jacobi法迭代矩阵 故， 故用Jacobi法解此方程组不收敛。 10、用追赶法求解线性方程组（过程可用分数计算） 解： ， ， ， 计算出后回代 x 0 1 2 3 y 5 -9 -5 35 11、已知函数的观测数据为右表： 利用Newton均差求解插值多项式 解： 其中 ， ， x 0.5 0.6 0.7 y -0.693 -0.511 -0.357 12、由的函数表用线性差值 用点0.5和0.6及二次插值计算ln0.54 （计算过程保留三位有效数字） 解： x 0 1 2 3 y 5 -9 -5 35 13、已知函数的观测数据为右表： 求解Lagrange插值多项式 解： 14、设已知函数的数值如下： ， ， ， 求不超过3次的多项式，使得满足插值条件。 解： ， ， ， ，，现构造 求出 ， 从而 由得： 故此有 15、设x与y之间的经验公式为，试根据下列的数据确定其中的常数a和b（计算过程保留3位有效数字） x 2 3 4 7 10 15 y 6 8 10 11 10 11 解：，令，，重新转换表格数据 ？ ？ ？ ？ ？ ？ ？ ？ ？ ？ ？ ？ 计算下表： 1 2 3 4 5 ∑ 由 得 16、给定数据 （i=0,1,2,3,4） x 1.00 1.25 1.50 1.75 2.00 y 5.10 5.79 6.53 7.45 8.46 求的最小二乘拟合曲线 解：不是多项式，但两端取对数得 若令，，则，它是线性最小二乘拟合问题 为求得A、b，先将化为，转化的数据为： 1.00 1.25 1.50 1.75 2.00 1.629 1.756 1.876 2.008 2.135 根据最小二乘法列表求方程组系数 故有法方程 解得：A=1.122 ，b=0.5056 ，a==3.071 最小二乘拟合曲线为 17、用最小二乘法求一个形如的经验公式，使他拟合下列数据并计算均方误差 19 25 31 38 44 19.0 32.3 49.0 73.3 97.8 解：本题给出拟合曲线