2015年电工学院数值分A_(有答案版).docVIP

电子科技大学2014-2015学年第 2 学期期 末 考试 A 卷 课程名称：MATLAB与数值分析 考试形式：开卷 考试日期：2015年 6月 日 考试时长：120分钟 课程成绩构成：平时 35 %， 期中 %， 实验 15 %， 期末 50 % 本试卷试题由 4 部分构成，共 7 页。 题号 一 二 三 四 合计 得分 一、填空题（共30分，共每题3分） 1.和Matlab执行命令 结果为 ，执行命令 结果为： ，执行命令结果 。 2. 写出Matlab程序中的命令： ，估计值的有效位数为估计值的有效位数为估计值的有效位数为 4．已知矩阵的holesky分解为，则方程的解为： 。 5. 已知矩阵，则矩阵行范数为为矩阵的谱半径为 6.设有方程组对于迭代求解方法收敛的充分必要条件是 插值函数必须经过插值点，而拟合函数不必经过拟合点 。 ．已知多项式插值节点为则拉格朗日插值为 。 9.采用(Simpson)求积公式求解的结果为 。 10. 取步长为，的前向欧拉迭代公式： 。（共0分，共每题分） × ） 3．采用Lagrange多项式插值时，插值多项式的阶数越高，插值效果越好。 （ × ） 4. 无论问题好坏，好的稳定算法都可得到其近似解。 （ × ） 5. 数值积分中Simpson求积中节点固定，而高斯求积中节点是不固定的。 （ √ ） 三、简答题（共12分） 阅读下列matlab程序，回答问题： 对于矩阵和向量阐述atlab代码所实现的功能，并分析算法的性 寻找上述代码中存在的错误，并改正。(6分) clear all close all clc A = [4, 2, 0, 2, 4, -1; 0, -1, 4]; b = [24; 30;24]; e = 0.001; D = [A(1,1), 0, 0; 0, A(2,2), 0; 0 0 A(3,3)]; M = eye(3) - inv(D)* A; g = inv(D) * b; x = [1;1;1]; for i = 1 : 1000 xtmp = x; x = M * x + g; if norm(x - xtmp, 2) e break; end 答：（1）阐述上述matlab代码的功能，并分析其性能 是利用Jacobi迭代法求解线性方程组方法（3分） 输入矩阵A严格对角占优，该迭代法能够收敛（3分） （2）寻找上述代码中存在的错误，并改正。 A赋值有误，括号第“，”改为分号“；”（3分） 迭代终止条件有误，“”改为“”（3分） （注：指出一处错误，分，共3分，改正正确1个，给分，共计3分。 四、计算题（共48分） 试确定求积公式 中参数、、，使其具有尽可能高的代数精度，并给出其代数精度。() 解：若令，则有 （1）（2分） 若令，，则有 （2）（2分） 若令，，则有 （3）（2分） 若令，，则有 （4）（2分） 于是，联立（1）、（2）、（3）、（4）式，即可解得方程组 解方程得到 （2分） 由上可知，确定的求积公式的代数精度，令，代入该求积公式验证可知，不成立。因此求积公式的代数精度。求积公式的代数精度为3。 （2分） 已知节点: 0 1 2 3 y 1 1.3499 2.226 4.4817 利用最小二乘法拟合模型中的参数和（分）取底数，得到 X 0 1 2 3 Y 1 1.3499 2.226 4.4817 ln(Y) 0 0.3 0.8 1.5 （3分） 构造方程组，写出拟合矩阵 ， （3分） 计算如下公式： ， （3分） 计算系数a、b分别为： ， 得到 （3分） 已知微分方程组， 1) 已知步长，写出该微分方程的前向欧拉迭代公式和梯形迭代公式（分） 写出该微分方程的前向欧拉方法前4步的迭代结果（分）解微分方程的前向欧拉迭代公式和梯形迭代公式前向欧拉公式：（分）梯形迭代公式：（分） （6分对方程组 该Jaco