计算方法试卷(I).doc

一、 (共20分) 已知函数的四组数据如下：试作一个三次Newton插值多项式 ；利用 求 的近似值，并作误差分析。已知函数的四组数据如下：试作一个三次Newton插值多项式 ；利用 求 的近似值，并作误差分析。解 1 ）由 Newton 公式 　　于是可得 　　2 ）审视所给的函数表，可以看出它是由函数 ，对 造出的。因此，求 ，即求 之值。现利用 求 的近似值，有 看来误差比较大。我们作误差分析。利用 Newton 余项公式 　　先补充求均差 ，这里是 ，计算结果如下（对应接上表的最末一行）： 　　于是 与 　　比较，可见理论分析相当吻合。此外，这里由于我们看出 ，因此，计算均差时使用 的精确值 。如果 的具体形式未给出，也可以 的近似值（如这里的 2 ）作误差估计用。 证明高斯－塞德尔法收敛； 写出高斯－塞德尔法迭代公式； 取初始值，求出。 因为严格对角占优矩阵，所以高斯－塞德尔迭代法收敛。 高斯－塞德尔法迭代公式为： 初值，计算得 三、用三角分解法解方程组（20分） （1）完成分解 ，，， ，，，，， 所以矩阵的三角分解 （2）解方程组，，， （3）解方程组，，， 所以 六、解微分方程初值问题。（共20分） 拟用四阶经典 方法及取 求解下列初值问题 　　试写出求解过程的算法描述。 拟用四阶经典 方法及取 求解下列初值问题 　　试写出求解过程的算法描述。 解 四阶 公式用于所给问题的计算公式及计算过程是　　其中， 　　因此，求解过程的算法描述如下（公式略作加工）： 　　打印表头：“ ”， “ ”， “ ” 　　打印输出： ， ， 　　 　　for do 　　{ 　　 　　 　　 　　 　　 　　打印输出： ， ， } y(0.1)=y1 K1=3y0/(1+x0)=3 K2=3(y0+0.05K1)/(1.05+x0)=3*1.15/1.05=3.2857 K3=3(y0+0.05K2)/(1.05+x0)=3*(1+0.05*3.2857)/1.05=3.3265 K4=3(y0+0.1K3)/(1.1+x0)=3(1+0.1*3.3265)/1.1=3.6345 y1=y0+(K1+2K2+2K3+K4)/6=3.3098 班 级_________________姓 名_________________学 号_________________ ……………………………………………装 订 线…………………………………………………… 常 州 工 学 院 试 卷 I 卷 共 4 页 第 4 页 09 / 10 学年第 二 学期 计算方法 试卷 考试类型 开卷 课程代码 0302004 班 级_________________姓 名_________________学 号_________________ ……………………………………………装 订 线…………………………………………………… 常 州 工 学 院 试 卷 I 卷 共 4 页 第 1 页 一 二 三 四 五 六 七 八 九 十 十一 十二 总分 KS002-1