數值分析试卷10计科专升本(B)卷参考答案.docVIP

数值分析2012－2013第1学期10计算专升期末试卷B 参考答案及评分标准 一、判断题（每小题2分，共20分） 1、T 2、F 3、T 4、T 5、T 6、T 7、T 8、T 9、T 10、T 二、计算题（每题8分，共40分） 设有微分方程。设步长为0.1，用Euler方法，计算的近似值：解：记步长为。Euler方法是以为起点，以为切线，构造直线，并以所构造直线在点处的值作为的近似，写成表达式有 （4分）依次计算的结果 （8分） 设，已知试用线性插值及二次插值计算的近似值，并估计误差。解：（1）因为0.54介于0.5与0.6之间，为了进行内插值，所以选取为插值节点，构造插值基函数： 插值函数为：并有余项：所以 （4分）（2）选取节点，构造二次插值基函数 有插值函数：有余项：所以 （8分） 设有实验数据 试求与的函数关系。解：由图上可以看出与大致呈线性关系。设  记，现在的目标是确定使达到最小。为此，令 写成矩阵的形式有 （5分）解之，得，即与的函数关系大致为 （8分） 用Newton迭代法求方程在上的一个根，并用做４次迭代计算。解：对求导，有，构造Newton迭代格式 （4分）作4次迭代的结果为 （8分） 用LU分解的方法求下列两个方程组的解与解：记，则原方程组可表示为 对作分解有 （4分）记，解得，求解，得。记，解得，求解，解得。 （8分）三、综合计算题（每题16分，共16分） 设有线性方程组试求（1）给出解线性方程组的Jacobi迭代、Gauss-Seidel迭代矩阵（2）判断解线性方程组的Jacobi迭代、Gauss-Seidel迭代的收敛性；（3）选取收敛速度较快的一种迭代方法，取进行四次迭代计算解：（1）分别记Jacobi迭代与Gauss-Seidel迭代矩阵为，则  （6分）（2）系数矩阵是严格对角占优的，所以Jacobi迭代、G-S迭代都是收敛的。 （10分）（3）的三个特征值为： 所以的谱半径为的三个特征值为 所以的谱半径为由于所以Gauss-Seidel迭代收敛速度快于Jacob迭代的收敛速度，取Gauss-Seidel迭代，以为初始向量，计算四次的结果为 （16分） 四、应用题（每题12分，共24分） 尝试只用加、减、乘、除四种运算，计算的近似值，并以实例说明。解：（1）因为是方程的根，所以可以用求根的方法求的近似值，记，则，构造迭代格式 这个算法中只含有乘法、除法与加法运算。另外，由此构造的序列，满足 其中，显然当时，总有，所以 （8分）（2）取，用上述迭代法，计算5次的结果为 （12分） 试用数值积分的方法计算的近似值。解：（1）因为，于是可以利用数值积分算法计算的近似值。如果用龙贝格算法，其算法过程为：将区间分为等分，计算： 计算 并以作为的近似值。 （8分）（2）以为例，采用复化梯型求积公式，其过程如下： ，记，，做4等分，节点为，对应的函数值为 所以 其误差为 （12分）蚅肂 试卷第 6 页 共 6 页 试卷第 1 页 共 6 页