资源描述:
《应用数值分析试卷-总.pdf》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、2004年数值分析试题2004.注:计算题取小数点后四位。1.(10分)利用Gauss-Legendre求积公式:1f(x)dx0.5556f(0.7746)0.8889f(0)0.5556f(0.7746)10导出求积分fxdx()的三点高斯型求积公式。32.(15分)写出求解线性代数方程组x2x2x5123xx12312xx7213的Gauss-Seidel迭代格式,并分析此格式的敛散性。210110303.(15分)设矩阵A,02013010(1)试计算
2、
3、A
4、
5、。(
6、2)用Householder变换阵H将A相似约化为上Hessenberg阵,即HAH为上Hessenberg阵。4.(10分)求关于点集1,2,3,4的正交多项式0(),(),x1x2()x。5.(10分)用最小二乘法确定一条经过原点的二次曲线,使之拟合下列数据xi1.02.03.04.0y0.81.51.82.0ix0134i6.(20分)给出数据点:y19156i(1)用xxx,,构造二次Lagrange插值多项式Lx(),并计算x1.5的近似值L(1.5)。01222(2)用xxx,,构造二次Newton插值多
7、项式Nx(),并计算x1.5的近似值N(1.5)。12322(3)用事后误差估计方法估计L(1.5)、N(1.5)的误差。2217.(10分)设矩阵A可逆,A为A的误差矩阵,证明:当A时,1AAA也可逆。8.(10分)设fx()四阶连续可导,xix0ihi,0,1,2.试建立如下数值微分''fx()2()0fx1fx()2公式fx()12h并推导该公式的截断误差数值分析答案33t1、(10分)解:作变换x,则xt[3,0][1,1],于是原积分可化为20113333fxdx()f(t)dtgtdt
8、()22223113[0.5556(0.7746)0.8889(0)0.5556(0.7746)]ggg20.8334(2.6619)1.3334(1.5)0.8334(0.3381)fff即为所求Gauss求积公式。2、(15分)解:方程组的Gauss-Seidel迭代格式为(k1)()k()kx52x2x123(kk1)(1)xx(1)/321xx(kk1)(22(1))/731其迭代矩阵为11022022B1300022G33
9、207004477其特征方程为22323021260207解之得260,1232126谱半径(B)1,故迭代发散。G2143、(15)解:(1)Aamaxij414ij1TT(2)设xy2,1,0,3,2,,0,0由xy且sgn()sgn(x)1可得2222TTu31设uxy0,3,0,3,w0,,0,u22100001300T13T22HI2wwI0,3,0,3
10、6000103310022131320222030从而HAH即为所求上Hessenberg阵。23123002200104、(10分)解:定义fxgx(),()在点集1,2,3,4上的内积为44(,)fgfxgx()()iifigi()()ii11设()1x,由schmitt正交化方法可得0(,x)0()xxx2.510(,)002222(,x0)(,x1)()xxx5x5201(,)(,
11、)00115、(10分)解:2设()xx,()xx,所求曲线为sx()a()xb()x1212110.8241.5,,Y,12391.84162.0则法方程为(,)(,)a(,Y)11121(,)(,)b(,Y)21222即30100a17.2,100354b55解之得a0.9497b0.1129于是所求曲线为2sx()a(
12、)xb()0.9497xx0.1129x。126、(20分)解:(1)由Lagrange插值得22Lx2()yl