数值分析期末试题

合集下载

数值分析期末考试题

数值分析期末考试题

数值分析期末考试题一、选择题1. 在数值分析中,用于求解线性方程组的雅可比方法属于以下哪种迭代法?A. 直接迭代法B. 间接迭代法C. 外推法D. 松弛法2. 插值法中,拉格朗日插值多项式的主要特点是?A. 适用于多项式插值B. 适用于函数值已知的情况C. 只适用于单点插值D. 适用于分段插值3. 在数值积分中,辛普森法则是一种?A. 单区间求积公式B. 双区间求积公式C. 三区间求积公式D. 多区间求积公式4. 误差分析中,截断误差通常与以下哪个概念相关?A. 舍入误差B. 舍入误差的补偿C. 条件数D. 病态条件5. 非线性方程求解中,牛顿法的收敛速度通常?A. 较慢B. 较快C. 与初始值有关D. 与方程的性质有关二、填空题1. 在求解三对角线性方程组时,托马斯算法是一种________方法。

2. 多项式插值中,牛顿插值多项式可以通过________法来构建。

3. 数值积分中,高斯求积法是一种________方法。

4. 误差传递的估计通常通过________公式来进行。

5. 非线性方程的求解中,二分法是一种________方法。

三、简答题1. 请简述数值分析中的条件数概念及其在解方程中的应用。

2. 描述线性方程组迭代法中的收敛性判断方法,并给出收敛域的计算公式。

3. 解释插值和拟合的区别,并举例说明各自的应用场景。

4. 阐述数值积分中梯形法则的原理及其误差估计方法。

5. 讨论非线性方程求解中不动点理论和收敛性的关系。

四、计算题1. 给定线性方程组如下,请使用高斯消元法求解未知数x、y、z的值: \[\begin{cases}2x + y + z = 6 \\x + 3y + 2z = 11 \\3x + y + 4z = 17\end{cases}\]2. 假设有一个函数f(x) = sin(x),给定插值节点如下,请使用拉格朗日插值法构造一个三次插值多项式,并计算在x=π/4处的插值误差。

数值分析试卷

数值分析试卷

数值分析期末考试试卷一、填空题:(15分)1、为了使计算球体体积时的相对误差不超过1%,测量半径R 的允许相对误差 为2、设矩阵A=⎪⎪⎭⎫ ⎝⎛3.01.05.06.0,则A 是行范数 ,列范数 ,2范数 , F 范数3、要使积分值dx x ⎰-202)-(sinx πβα最小,则a= , β= 4、若矩阵A 是正交阵,则有cond (A )=5、对于矩阵A=⎪⎪⎪⎭⎫ ⎝⎛111a a a a a a ,若A 为正定矩阵则a的范围是 ,若A 只对J 法收敛,则a的范围是二、选择题:(15分)1、设x>0,x 的相对误差为δ,则lnx 的相对误差为()A.ln δB.δ1 C.δ D.2δ 2、以下哪个公式不需要做变换就可以避免有效数字的损失()A.212-x e B.sinx-siny C.arctanx-arctany D.都要做变换 3、计算过程中取有限位数字进行运算而引起的误差称为()A.模型误差B.截断误差C.舍入误差D.观测误差4、Legendre 多项式前五项根的范围是()A.(-1,1)B.(0,+∞)C.(-∞,0)D.(-∞,+∞)5、三、判断题:(15分)1、复化梯形求积公式和复化Simpson 求积公式都是二阶收敛的 ()2、三点Gauss-Legendre 求积公式具有5次代数精度 ()3、如果矩阵A 的按绝对值最大特征值和最小特征值之比很小,则A 是病态的()4、方程组Ax=b 的SOR 法收敛有20<<ω,则20<<ω时,方程组Ax=b 的SOR 法收敛 ()5、映内性既可保证不动点的存在,也可保证其唯一性 ()6、e=2.71828......,x=e ,则A x =2.7有两位有效数字 ()7、若A n n R ⨯∈矩阵,当B=PA PT 且P 为正交阵,则有F F A B = ()8、若矩阵A 为严格对角占优矩阵,则解方程Ax=b 的J 法和GS 法都收敛 ()9、设P n n R ⨯∈且非奇异,X 为R n 上的一种向量范数,则p X 也是R n 上的一种向量范数。

《数值分析》A卷期末考试试题及参考答案

《数值分析》A卷期末考试试题及参考答案

一、单项选择题(每小题3分,共15分) 1、用Simpson 公式求积分1401x dx +⎰的近似值为 ( ).A.2924 B.2429C.65D. 562、已知(1)0.401f =,且用梯形公式计算积分2()f x dx ⎰的近似值10.864T =,若将区间[0,2]二等分,则用递推公式计算近似值2T 等于( ). A.0.824 B.0.401 C.0.864 D. 0.8333、设3()32=+f x x ,则差商0123[,,,]f x x x x 等于( ).A.0B.9C.3D. 64的近似值的绝对误差小于0.01%,要取多少位有效数字( ). A.3 B.4 C.5 D. 25、用二分法求方程()0=f x 在区间[1,2]上的一个实根,若要求准确到小数 点后第四位,则至少二分区间多少次( ).A.12B.13C.14D. 15二、填空题(每小题4分,共40分)1、对于迭代函数2()=(3)ϕ+-x x a x ,要使迭代公式1=()ϕ+k k x x则a 的取值范围为 .2、假设按四舍五入的近似值为2.312,则该近似值的绝对误差限为 .3、迭代公式212(3)=,03++>+k k k k x x a x a x a收敛于α= (0)α>. 4、解方程4()530f x x x =+-=的牛顿迭代公式为 . 5、设()f x 在[1,1]-上具有2阶连续导数,[1,1]x ∀∈-,有1()2f x ''≤,则()f x 在[1,1]-上的线性插值函数1()L x 在点0处的误差限1(0)R ≤______.6、求解微分方程初值问题2(0)1'=-⎧⎨=⎩y xy yy ,0x 1≤≤的向前Euler 格式为 .7、设310131013A -⎛⎫⎪=-- ⎪ ⎪-⎝⎭,则A ∞= .8、用梯形公式计算积分112-⎰dx x 的近似值为 . 9、设12A 21+⎡⎤=⎢⎥⎣⎦a 可作Cholesky 分解,则a 的取值范围为 . 10、设(0)1,(0.5) 1.5,(1)2,(1.5) 2.5,(2) 3.4f f f f f =====,若1=h ,则用三点公式计算(1)'≈f .三、解答题(共45分) 1、给定数据用复化Simpson 公式计算1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛. (5分)4、已知数据试对数据用最小二乘法求出形如=+y x b的拟合曲线. (8分) 5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (8分) 6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦一、单项选择题(每小题3分,合计15分) 1、A 2、D 3、C 4、C 5、D 二、填空题(每小题3分,合计30分) 1、0<<a ; 2、31102-⨯; 3;4、4135345++-=-+k k k k k x x x x x ; 5、14; 6、1(2)+=+-n n n n n y y h x y y ; 7、5;8、34-; 9、3>a ;10、1.2;三、计算题(合计55分) 1、给定数据用复化Simpson 公式计算 1.381.30()f x dx ⎰的近似值,并估计误差,小数点后保留3位. (8分)解: 401024S [()4()()]6-=++x x f x f x f x ………… 1分 1.38 1.30(3.624 4.20 5.19)6-=+⨯+ 0.341= ………… 2分20422012234S [()4()()][()4()()]66--=+++++x x x xf x f x f x f x f x f x =0.342 ………… 6分2211[]15-≈-I S S S =-⨯40.6710 ………… 8分 2、用直接三角分解法求线性代数方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡432631531321321x x x 的解. (8分) 解:设111213212223313233u u u 123100135l 100u u 136l l 100u ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=*⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦………… 1分 111=u ,212=u ,313=u ,121=l ,131=l 122=u ,223=u ,132=l133=u ,133=l …………6分所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111011001L ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100210321U …………7分 由b Ly =得Ty )1,1,2(=;由y Ux =得Tx )1,1,1(-=. ………… 8分3、求()λx ,使得迭代公式1()()λ+=+k k k k f x x x x 求方程2()31=+-f x x x 的根的相应迭代序列{}k x 具有平方收敛.(6分)解:要使迭代序列具有平方收敛,则()0ϕ'*=x ………… 2分 而()()()ϕλ=+f x x x x ,即 ………… 3分 2()()()()10()λλλ''**-**+=*f x x x f x x …………4分 而()0*=f x 则有()1()λ'*=-*f x x ………… 5分所以()()23λ'=-=--x f x x ………… 6分4、已知数据试对数据用最小二乘法求出形如=+ay x b的拟合曲线. (8分) 解:因为11=+b x y a a ,令0111,,,====b a a y x x a a y……2分 则有法方程01461061410⎛⎫⎛⎫⎛⎫= ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭a a ……5分解出014,1==-a a ,则1,4=-=-a b ……7分 所以1=4-y x……8分5、已知(2)8f -=,(0)4f =,(2)8=f ,试求二次拉格朗日插值多项式. (7分)解:01()(2)8l x x x =- …………2分 211()(4)4l x x =-- …………4分21()(2)8l x x x =+ …………6分 2012()()(2)()(0)()(2)L x l x f l x f l x f =-++24=+x …………7分6、设矩阵A 如下,根据谱半径判断用Jacobi 迭代法求解方程组Ax b =的敛散性.(8分)1102111221012A ⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦解:100010001D ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,00010021002L ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦,10021002000U ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………3分1100211()0221002J B D L U -⎡⎤⎢⎥⎢⎥⎢⎥=+=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦…………5分 2102111()0222102J E B λλλλλλ⎡⎤-⎢⎥⎢⎥⎢⎥-=--=-=⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦…………6分()2J B ρ=…………7分 所以用Jacobi 迭代法求解方程组Ax b =收敛 …………8分。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案试题一:1. 简答题(共10分)a) 什么是数值分析?它的主要应用领域是什么?b) 请简要解释迭代法和直接法在数值计算中的区别。

2. 填空题(共10分)a) 欧拉方法是一种______型的数值解法。

b) 二分法是一种______法则。

c) 梯形法则是一种______型的数值积分方法。

3. 计算题(共80分)将以下函数进行数值求解:a) 通过使用二分法求解方程 f(x) = x^3 - 4x - 9 = 0 的近似解。

b) 利用欧拉方法求解微分方程 dy/dx = x^2 + 2x + 1, y(0) = 1 在 x = 1 处的解。

c) 使用梯形法则计算积分∫[0, π/4] sin(x) dx 的近似值。

试题二:1. 简答题(共10分)a) 请解释什么是舍入误差,并描述它在数值计算中的影响。

b) 请解释牛顿插值多项式的概念及其应用。

2. 填空题(共10分)a) 数值稳定性通过______号检查。

b) 龙格-库塔法是一种______计算方法。

c) 零点的迭代法在本质上是将方程______转化为______方程。

3. 计算题(共80分)使用牛顿插值多项式进行以下计算:a) 已知插值节点 (-2, 1), (-1, 1), (0, 2), (1, 4),求在 x = 0.5 处的插值多项式值。

b) 已知插值节点 (0, 1), (1, 2), (3, 7),求插值多项式,并计算在 x = 2 处的值。

c) 使用 4 阶龙格-库塔法求解微分方程 dy/dx = x^2 + 1, y(0) = 1。

答案:试题一:1. a) 数值分析是研究使用数值方法解决数学问题的一门学科。

它的主要应用领域包括数值微积分、数值代数、插值和逼近、求解非线性方程、数值积分和数值解微分方程等。

b) 迭代法和直接法是数值计算中常用的两种方法。

迭代法通过反复迭代逼近解,直到满足所需精度为止;而直接法则通过一系列代数运算直接得到解。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案一、选择题(每题5分,共20分)1. 在数值分析中,下列哪个算法不是用于求解线性方程组的?A. 高斯消元法B. 牛顿法C. 雅可比法D. 追赶法答案:B2. 插值法中,拉格朗日插值法属于:A. 多项式插值B. 样条插值C. 线性插值D. 非线性插值答案:A3. 以下哪个选项不是数值分析中的误差来源?A. 截断误差B. 舍入误差C. 计算误差D. 测量误差答案:C4. 在数值积分中,梯形法则的误差项是:A. O(h^2)B. O(h^3)C. O(h)D. O(1)答案:A二、填空题(每题5分,共20分)1. 牛顿插值法中,插值多项式的一般形式为:______。

答案:f(x) = a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1) + ...2. 牛顿迭代法求解方程的根时,迭代公式为:x_{n+1} = x_n -f(x_n) / __________。

答案:f'(x_n)3. 在数值分析中,______ 用于衡量函数在区间上的近似积分值与真实积分值之间的差异。

答案:误差4. 线性方程组的解法中,______ 法是利用矩阵的LU分解来求解。

答案:克兰特三、解答题(每题10分,共60分)1. 给定函数f(x) = e^(-x),使用拉格朗日插值法,求x = 0.5时的插值值。

解答:首先选取插值节点x_0 = 0, x_1 = 0.5, x_2 = 1,对应的函数值分别为f(0) = 1, f(0.5) = e^(-0.5), f(1) = e^(-1)。

拉格朗日插值多项式为:L(x) = f(0) * (x-0.5)(x-1) / (0-0.5)(0-1) + f(0.5) * (x-0)(x-1) / (0.5-0)(0.5-1) + f(1) * (x-0)(x-0.5) / (1-0)(1-0.5)将x = 0.5代入得:L(0.5) = 1 * (0.5-0.5)(0.5-1) / (0-0.5)(0-1) + e^(-0.5) * (0.5-0)(0.5-1) / (0.5-0)(0.5-1) + e^(-1) * (0.5-0)(0.5-0.5) / (1-0)(1-0.5)计算得L(0.5) = e^(-0.5)。

数值分析期末复习题答案

数值分析期末复习题答案

数值分析期末复习题答案一、选择题1. 以下哪个算法是用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 共轭梯度法D. 辛普森积分法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的主要区别是什么?A. 插值点的选取不同B. 插值多项式的构造方式不同C. 计算复杂度不同D. 适用的函数类型不同答案:B3. 在数值积分中,梯形法则和辛普森法则的主要区别是什么?A. 精度不同B. 适用的积分区间不同C. 计算方法不同D. 稳定性不同答案:A二、简答题1. 解释什么是数值稳定性,并举例说明。

答案:数值稳定性指的是数值方法在计算过程中对于舍入误差的敏感程度。

例如,在求解线性方程组时,如果系数矩阵的条件数很大,则该方程组的数值解对舍入误差非常敏感,即数值稳定性差。

2. 说明数值微分与数值积分的区别。

答案:数值微分是估计函数在某一点的导数,而数值积分是估计函数在某个区间上的积分。

数值微分通常用于求解函数的局部变化率,而数值积分用于求解函数在一定区间内的累积效果。

三、计算题1. 给定一组数据点:(1, 2), (2, 3), (3, 5), (4, 6),请使用拉格朗日插值法构造一个三次插值多项式。

答案:首先写出拉格朗日插值基函数,然后根据数据点构造插值多项式。

具体计算过程略。

2. 给定函数 f(x) = x^2,使用牛顿-科特斯公式中的辛普森积分法在区间 [0, 1] 上估计积分值。

答案:首先确定区间划分,然后应用辛普森积分公式进行计算。

具体计算过程略。

四、论述题1. 论述数值分析中误差的来源及其控制方法。

答案:误差主要来源于舍入误差和截断误差。

舍入误差是由于计算机在进行浮点数运算时的精度限制造成的,而截断误差是由于数值方法的近似性质导致的。

控制误差的方法包括使用高精度的数据类型、选择合适的数值方法、增加计算步骤等。

五、综合应用题1. 给定一个线性方程组 Ax = b,其中 A 是一个 3x3 的矩阵,b 是一个列向量。

《数值分析》期末复习题(1)

《数值分析》期末复习题(1)

《数值分析》期末复习题一、单项选择题1. 数值x *的近似值x =0.32502×10-1,若x 有5位有效数字,则≤-*x x ( ).(A)21×10-3 (B) 21×10-4 (C) 21×10-5 (D) 21×10-6 2. 设矩阵A =10212104135⎡⎤⎢⎥⎢⎥⎣⎦,那么以A 为系数矩阵的线性方程组A X =b 的雅可比迭代矩阵为( )(A)00.20.10.200.40.20.60--⎡⎤⎢⎥--⎢⎥--⎣⎦(B)10.20.10.210.40.20.61⎡⎤⎢⎥⎢⎥⎣⎦(C) 00.20.10.200.40.20.60⎡⎤⎢⎥⎢⎥⎣⎦ (D)021204130⎡⎤⎢⎥⎢⎥⎣⎦3. 已知(1)1,(2)4,(3)9f f f ===,用拉格朗日2次插值,则(2.5)f =( )(A) 6.15 (B) 6.25 (C) 6.20 (D) 6.10 4. 抛物形求积公式的代数精度是( )A. 1,B. 2 ,C. 3,D. 45. 改进欧拉格式的 局部截断误差是( ). (),A O h 2. (),B O h 3. (),C O h 4. ().D O h二、填空题1、以722作为π的近似值,它有( )位有效数字; 2、经过)1,2( ),2,1( ),1,0(C B A 三个节点的插值多项式为( ); 3、用高斯-赛德尔迭代法解方程组⎩⎨⎧-=+-=+,10,232121x bx bx x 其中b 为实数,则方法收敛的充分条件是b 满足条件( );4、取步长为1.0=h ,用欧拉法计算初值问题22',(0)0,y x y y ⎧=+⎨=⎩的解函数)(x y ,它在3.0=x 的近似值为( );5、已知方程0sin 1=--x x 在)1,0(有一个根,使用二分法求误差不大于41021-⨯的近似解至少需要经过( )次迭代。

数值分析期末考试题及答案

数值分析期末考试题及答案

数值分析期末考试题及答案一、选择题(每题2分,共20分)1. 在数值分析中,下列哪个算法用于求解线性方程组?A. 牛顿法B. 高斯消元法C. 插值法D. 傅里叶变换答案:B2. 以下哪个选项不是数值分析中的误差类型?A. 舍入误差B. 截断误差C. 测量误差D. 累积误差答案:C3. 多项式插值中,拉格朗日插值法的特点是:A. 插值点必须等距分布B. 插值多项式的次数与插值点的个数相同C. 插值多项式是唯一的D. 插值多项式在插值点处的值都为1答案:B4. 在数值分析中,下列哪个方法用于求解非线性方程?A. 辛普森法则B. 牛顿迭代法C. 欧拉法D. 龙格-库塔法答案:B5. 以下哪个是数值稳定性的指标?A. 收敛性B. 收敛速度C. 条件数D. 误差传播答案:C二、简答题(每题10分,共20分)1. 简述高斯消元法求解线性方程组的基本原理。

答案:高斯消元法是一种直接解法,通过行变换将增广矩阵转换为上三角形式,然后通过回代求解线性方程组。

它包括三个基本操作:行交换、行乘以非零常数、行相加。

2. 解释什么是数值稳定性,并举例说明。

答案:数值稳定性是指数值解对输入数据小的扰动不敏感的性质。

例如,某些数值方法在计算过程中可能会放大舍入误差,导致结果不可靠,这样的方法就被认为是数值不稳定的。

三、计算题(每题15分,共30分)1. 给定线性方程组:\[\begin{align*}x + 2y - z &= 4 \\3x - y + 2z &= 1 \\-x + y + z &= 2\end{align*}\]使用高斯消元法求解该方程组,并给出解。

答案:首先将增广矩阵转换为上三角形式,然后回代求解,得到\( x = 1, y = 2, z = 1 \)。

2. 给定函数 \( f(x) = x^2 - 3x + 2 \),使用拉格朗日插值法在\( x = 0, 1, 2 \) 处插值,并求出插值多项式。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数值分析期末试题一、填空题(20102=⨯分)(1)设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=283012251A ,则=∞A ______13_______。

(2)对于方程组⎩⎨⎧=-=-34101522121x x x x ,Jacobi 迭代法的迭代矩阵是=J B ⎥⎦⎤⎢⎣⎡05.25.20。

(3)3*x 的相对误差约是*x 的相对误差的31倍。

(4)求方程)(x f x =根的牛顿迭代公式是)('1)(1n n n n n x f x f x x x +--=+。

(5)设1)(3-+=x x x f ,则差商=]3,2,1,0[f 1 。

(6)设n n ⨯矩阵G 的特征值是n λλλ,,,21 ,则矩阵G 的谱半径=)(G ρi ni λ≤≤1max 。

(7)已知⎥⎦⎤⎢⎣⎡=1021A ,则条件数=∞)(A Cond 9(8)为了提高数值计算精度,当正数x 充分大时,应将)1ln(2--x x 改写为)1ln(2++-x x 。

(9)n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。

(10)拟合三点))(,(11x f x ,))(,(22x f x ,))(,(33x f x 的水平直线是)(3131∑==i i x f y 。

二、(10分)证明:方程组⎪⎩⎪⎨⎧=-+=++=+-12112321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。

证明:Jacobi 迭代法的迭代矩阵为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=05.05.01015.05.00J BJ B 的特征多项式为)25.1(5.05.0115.05.0)det(2+=---=-λλλλλλj B IJ B 的特征值为01=λ,i 25.12=λ,i 25.13-=λ,故25.1)(=J B ρ>1,因而迭代法不收敛性。

三、(10分)定义内积⎰=1)()(),(dx x g x f g f试在{}x Span H ,11=中寻求对于x x f =)(的最佳平方逼近元素)(x p 。

解:1)(0≡x ϕ,x x ≡)(1ϕ,1),(100==⎰dx ϕϕ,21),(101==⎰xdx ϕϕ,31),(1211==⎰dx x ϕϕ,32),(10==⎰dx x f ϕ,52),(11==⎰dx x x f ϕ。

法方程⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡5232312121110c c 解得1540=c ,15121=c 。

所求的最佳平方逼近元素为 x x p 1512154)(+=,10≤≤x 四、(试用三次多项式以最小二乘法拟合所给数据。

解:332210)(x c x c x c c x y +++=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----=84211111000111118421A , ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=130034003401034010001005A A TT T y A )4.14,7,2.4,9.2(=法方程y A Ac A T T =的解为4086.00=c ,39167.01=c ,0857.02=c ,00833.03=c 得到三次多项式3200833.00857.039167.04086.0)(x x x x y +++=误差平方和为000194.03=σ五. (10分) 依据如下函数值表建立不超过三次的Lagrange 插值多项式,用它计算)2.2(f ,并在假设1)()4(≤x f 下,估计计算误差。

解:先计算插值基函数1478781)40)(20)(10()4)(2)(1()(230+-+-=------=x x x x x x x lx x x x x x x l 38231)41)(21)(01()4)(2)(0()(231+-=------=x x x x x x x l -+-=------=2324541)42)(12)(02()4)(1)(0()(x x x x x x x l 12181241)24)(14)(04()2)(1)(0()(233+-=------=所求Lagrange 插值多项式为121445411)(3)(23)(9)()()()(233210303+-+-=+++==∑=x x x x l x l x l x l x l x f x L i i i 从而0683.25)2.2()2.2(3=≈L f 。

据误差公式))()()((!4)()(3210)4(3x x x x x x x x f x R ----=ξ及假设1)()4(≤x f 得误差估计:0396.09504.0!41)42.2)(22.2)(12.2)(02.2(!4)()()4(3=⨯≤----=ξf x R 六. (10分) 用矩阵的直接三角分解法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡71735 30103421101002014321x x x x解 设⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡443433242322434241323121020111113010342110100201u u u u u ul l l l l l由矩阵乘法可求出ij u 和ij l⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡10101211011111434241323121l l l l l l ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡21210102010201443433242322u u u u u u解下三角方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡7173510101211014321y y y y 有51=y ,32=y ,63=y ,44=y 。

再解上三角方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡463521************x x x x 得原方程组的解为11=x ,12=x ,23=x ,24=x 。

七. (10分) 试用Simpson 公式计算积分dx e x⎰211的近似值, 并估计截断误差。

解:0263.2)4(612215.11211=++-≈⎰e e e dx e xx e xx x x f15678)4()2436121(+++=43.198)1()(max )4()4(21==≤≤f x f x截断误差为06890.0)(max 2880)12()4(2152=-≤≤≤x f R x八. (10分) 用Newton 法求方程2ln =-x x 在区间) ,2(∞内的根, 要求8110--<-kk k x x x 。

解:此方程在区间) ,2(∞内只有一个根s ,而且在区间(2,4)内。

设2ln )(--=x x x f则 x x f 11)('-=, 21)(''xx f = Newton 法迭代公式为1)ln 1(112ln 1-+=----=+k k k kk k k k x x x x x x x x , ,2,1,0=k 取30=x ,得146193221.34=≈x s 。

九. (10分) 给定数表求次数不高于5的多项式)(5x H ,使其满足条件⎩⎨⎧====2 ,0),()(3,2 ,1 ,0 ),()('55i x f x H i x f x H i i i i 其中,1i x i +-= 3 ,2 ,1 ,0=i 。

解:先建立满足条件)()(3i x f x p =, 3,2,1,0=i的三次插值多项式)(3x p 。

采用Newton 插值多项式[][]))((,,)(,)()(1021001003x x x x x x x f x x x x f x f x p --+-+=+[]))()((,,,2103210x x x x x x x x x x f ---)1()1(61)1()1(410-+-+-++=x x x x x x326161914x x x --+= 再设 )2)(1()1)(()()(35--+++=x x x x b ax x p x H ,由⎩⎨⎧=-++==-+-+-=-1.0)2)(()1()1(1)6)(()1()1('3'5'3'5b a p H b a p H 得⎪⎩⎪⎨⎧=+=+-6017811b a b a 解得36059-=a ,360161=b 。

故所求的插值多项式)2)(1()59161(36016161914)(2325---+--+=x x x x x x x x H。

相关文档
最新文档