[精品]数值分析试题A09

1、中国石油大学(北京)2008 2009学年第一学期研究生期末考试试题 A(闭卷考试)课程名称：数值分析所有试题答案写在答题纸上，答案写在试卷上无效题号-一一-二二三四五七八总分得分注：计算题取小数点后四位一、填空题(共30分，每空3分)1、 已知Xk(k = 0,1,-, n) 是互异节点，(x)是对应节点的Lagrange插值基函数，P(x)是任意一个首项系数为1的 + 1次多项式，则 P(x) 一 P )()=k=O、.-十 af x 3 + x 2,0 x 12、 设分段多项式S(x) = 322x + bx + ex -1, lx2是以0,1,2为节点的三次样条函数，贝=, c =o3

2、、 如果A是正交矩阵，则Conc2(A)=。4、 用二3.141作为的近似值，贝 Ux有位有效数字，其绝对误差限为o3 r5、 数值积分公式f(x)dxR f(l) + f(2)是否为插值型求积公式：一，其代数jo2精度为O6、下歹matlab程序中s2计算的是，并指明si与s2的区别为其中：aex = axlO x;a,x e R ot=O;s2=lel4;for i=l:le6temp= l/(le3+i); t=t+temp; s2=s2+temp;endsl= t+lel4;(8分)已知函数表0121011试利用重节点 Newt on差商构造满足插值条件P(0) = 1,P(1) =

3、0,P=1,P(2) = 1,的三次多项式P(.r),(要求构造出差商表)三、(8 分)已知向量 x = (2,0,2,l) r，试构造 Householder变换阵，使 Hx = (Q,Q,k,Q)其中k w R。旦=一 1),试利用勒让德正交多项四、(12分)已知勒让德(Legendre)正交多项式 R = l,g = X,式在二次多项式类H= span中求一个多项式 S(x),使其成为f (x) = e x在-1,1 上的最佳平方逼近函数，并计算岀平方误差。五、(10分)写岀求解线性代数方程组X - 2 矛 2 + 2 矛 3 = 5七、(12分)给岀计算x =72 + A/2 + V2

4、 + ?-的迭代格式，讨论迭代格式的收敛性，并证明x =2o八、(8分)求解常微分方程初值问题y = f(x,y)0,故取 k = % 3 o（3 分）y = C3e3 = ke 3 =（ 0,0,-3,0）1T111001-1002H =1 UUTP15-100-10-20-55方法 2:取 k = 3 ,y = (0,0,3,U=x y = (2,0,1,1-02 -13022021-012(3分)(3(22四、(12 分)spanl, x2 = spanP,P2,设 &x) =+c 2P2(30,鸟）（乌的）I丫勺二鸟）（旦旦）5时八。2 j2co =e-e? 2 .350388即：2一

5、 c, =e 7eT ? 0.1431245 -c =?(e eT)Rl.175194解得:1C2 =|(e-7A)?0.357813分）1 一S(x) CqP+c_r = c(e e ) + (s 5e )x(3x 1) = 0.536715., +0.9962893分）2p|= (/, /) - (/, s)=-号疽 + 36 = 0.8135五、（10分）方程组的Gauss-Seidel 迭代格式(5分)过 5 =5 + 2 妲 *)+2 廿 *) 可即)=(+ X (s) /3 / ) = ( 2_2 妒5) /7其迭代矩阵为六、(12 分)方法 1 :4=2,l 2 =224 3 =

6、丫 3 = 2 ()X (1)=,33333100f 100 0、000110、20100341-1一 2 -202-2-130 002_2BG33=420 70047其特征方程为72-2 232-2 32 0 =2U-262 =022 0 72解之得A=?i 2=0,/i3=26(2谱半径p(Bj = 1 ,故迭代发散。(8 由 Ly =5 得 y = 1, !，:,j./曰4 3 2 1 田以=(2七、(12分)由题意可得岀其迭代格式为迎+i = J2 + * 且0叫v 2伊 3)= + x,当 g0, 2吐伊(工)g 0,2(3 分)1-119,二(2 + x).22 A/ 2 + x当时， |9(x)|=00解得：足=一 1, =2.其中x; =-l 0舍去。可得x=2.即解得 x = 2.(3 分)八、(8 分)假定 y, = y(x,), 将改进欧拉公式写成：hy+i =)+a D a” )+ , )+ 仪yn+i = yA + -lyXA + f(x n,y(x 卜 + hfXx n,y(.x nA + hy,(xJfy(x n,y(x nAO(h )=y(x n ) + hyxJ + y x n) + ?(/z )(3 分)另一方面，依 Tay