数值方法试卷及其答案.docVIP

《数值计算与MATLAB语言》2003-2004 第1学期试卷A (一页开卷考试) ************************** 数值计算试题 *********************************** 取y0=30, 按递推公式 去计算y100, 若取(五位有效数字), 试问计算y100将有多大的误差? [8分] 2. 用对分区间法求解方程x3-2x2-1=0在2.2附近的实根，准确到三位有效数字。[9分] 3. 用LU分解求下列方程组：[8分] 4. 对下述矩阵计算 : [4分×2] 5. ，请计算它的条件数和特征值. [8分] 6. 给定 xi = i +1, i = 0, 1, 2, 3, 4, 5. 下面哪个是 l2(x)的图像，并简要说明理由？[4分] 7. 用拉格郎日插值找经过点(-3,-1),(0,0),(3,1),(6,2)的三次插值公式 [10分] 8、用最小二乘原理求一个形如y=a+bx2的经验公式,使与下列数据相拟合[10分] x 19 25 31 38 44 y 19 32.3 49 73.3 97.8 9、已知数表 x 1.8 2.0 2.2 2.4 2.6 f(x) 3.12014 4.42569 6.04241 8.03014 10.46675 计算积分（12分） 用复化梯形公式； 用复化辛甫生公式； 用牛顿-柯特斯公式。 10、用欧拉法求解下列常微分方程（取h=0.2，在区间[0,0.6]上计算）（8分） ***************************** MATLAB语言试题 ***************************** 11. 传入一个方阵A，编写一个自定义函数 “test1” 完成如下的功能： [8分] 把A的主对角线元素写入文本文件 “abc.txt” 中； 把B中的元素进行随机初始化（B的行列数与A相同）； 从文本 “abc.txt” 中读取元素到矩阵C，把C中每个元素值均加100； 分别计算：D=B*C；E=B.*C； 把矩阵D和E作为传出的参数传出。 12. 编写一个自定义函数求解任一函数的导数, 要求函数形式为 dy=ddf(x,e). [7分] 其中, e为误差精度, dy为x对应的f(x)的导数值. f(x)为被求导的函数. 《数值计算与MATLAB语言》2003-2004 第1学期试卷A参考答案 (一页开卷考试) ************************** 数值计算试题 *********************************** 1. 取y0=30, 按递推公式 去计算y100, 若取(五位有效数字), 试问计算y100将有多大的误差? [8分] 解: 2. 用对分区间法求解方程x3-2x2-1=0在2.2附近的实根，准确到三位有效数字。[9分] 解：对方程左边的函数进行预处理知在2.2附近取x1=2.1，x2=2.3，对应的y1=-0.559，y2=0.587，符合对分区间法的要求 [3分] 。迭代过程如下：[4分] 因此，准确到三位有效数字，得方程的根为x=2.203125。[2分] 3. 用LU分解求下列方程组：[8分] 解：由LU分解的元素计算公式，对系数矩阵进行分解得L和U矩阵如下：[4分] 由LY=b进行回代，得：Y=(1,1.5,4)T 。[2分] 代入方程组UX=Y得：X=(49,-29,-2)T。[2分] 4. 对下述矩阵计算 : [4分×2] 解： （1）6, 7, 5.465 （2）12, 23, 15.6525 5. ，请计算它的条件数和特征值. [8分] 解：考察 A 的特征根： 6. 给定 xi = i +1, i = 0, 1, 2, 3, 4, 5. 下面哪个是 l2(x)的图像，并简要说明理由？[4分] 解：C图像。(略) 7. 用拉格郎日插值找经过点(-3,-1),(0,0),(3,1),(6,2)的三次插值公式 [10分] 解：按照拉格朗日插值公式，得到如下的表达式：（每个li表达式各占2分，结论2分） 8、用最小二乘原理求一个形如y=a+bx2的经验公式,使与下列数据相拟合（10分） x 19 25 31 38 44 y 19 32.3 49 73.3 97.8 解：按照最小二乘原理，得：[原理2分；正规方程4分；系数解2分；最终表达式2分] 所以，表达式为：y=1.01+0.05x2 9、已知数表 x 1.8 2.0 2.2 2.4 2.6 f(x) 3.12014 4.42569 6.04241 8.03014 10.46675 计算积分（12分） 用复化梯形公式； 用复化辛甫生公式； 用牛顿-柯特斯公式。 解：h=0.2