数值分析期末复习题

合集下载

山东交通学院数值分析期末复习题

山东交通学院数值分析期末复习题

数值分析A 卷复习题一、填空题(每题3分,共30分)1.已知 3.201,0.57a b ==是经过四舍五入后得到的近似值,则a b +有( 2 )位有效数字。

2.设()i l x 是以,(0,1,,9)k x k k ==为节点的Lagrange 插值基函数,则9()==∑ik kl k ( i )。

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

4.设有矩阵2304A -⎡⎤=⎢⎥⎣⎦,则1A =( 7 )。

5.三次样条函数是在各个子区间上的( 3 )次多项式。

6.用牛顿下山法求解方程303x x -=的下山条件是( 1()()k k f x f x +≤ )。

7.已知3n =时的Newton-Cotes 系数(3)018C =,则(3)3C =( 1/8 )。

8.若线性代数方程组Ax b =的系数矩阵A 为严格对角占优阵,则雅可比迭代和高斯-塞德尔迭代都 ( 收敛 )。

9.用Gauss-Seidel 迭代法解方程组1212423x ax ax x +=⎧⎨+=-⎩,其中a 为实数,方法收敛的充要条件是a 满足a <<。

二、单选题(每题3分,共30分)1.用1+x 近似表示xe 所产生的误差是( C )误差。

A.模型B.观测C.截断D.舍入 2.已知数21234721,0.721,0.700,710x x x x -====⨯是由四舍五入得到的,则它们的有效数字的位数应分别为( A )。

A.3,3,3,1B.3,3,3,3C.3,3,1,1D.3,3,3,2 3.设(1)1,(0)3,(2)4-===f f f ,则抛物线插值多项式中2x 的系数为( A )。

A.-0.5B.0.5C.2D.-24. 为求方程3210x x --=在区间[]1.3,1.6内的一个根,把方程改写成下列形式,并建立相应的迭代公式,迭代公式不收敛的是( A )。

A.2111,:k x x x +==-迭代公式1221111,:k kx x x x +=+=+迭代公式C.32213111/,:()k kx x x x +=+=+迭代公式 D.23212111,:k k kk xx x x x x +-==+++迭代公式5.线性方程组AXB =能用高斯消元法求解的充分必要条件为( D )。

数值分析复习题及答案

数值分析复习题及答案

数值分析复习题一、选择题1. 3.142和3.141分别作为π的近似数具有〔 〕和〔 〕位有效数字.A .4和3B .3和2C .3和4D .4和42. 求积公式()()211211()(2)636f x dx f Af f ≈++⎰,那么A =〔 〕A . 16B .13C .12D .233. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足〔 〕A .()00l x =0,()110l x = B .()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =4. 设求方程()0f x =的根的牛顿法收敛,那么它具有〔 〕敛速。

A .超线性B .平方C .线性D .三次5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=⎧⎪++=⎨⎪--=⎩ 作第一次消元后得到的第3个方程〔 〕.A .232x x -+= B .232 1.5 3.5x x -+= C .2323x x -+= D .230.5 1.5x x -=-二、填空1. 设2.3149541...x *=,取5位有效数字,那么所得的近似值x=.2.设一阶差商()()()21122114,321f x f x f x x x x --===---,()()()322332615,422f x f x f x x x x --===-- 那么二阶差商()123,,______f x x x =3. 设(2,3,1)TX =--, 那么2||||X = ,=∞||||X 。

4.求方程 21.250x x --= 的近似根,用迭代公式x =01x =, 那么 1______x =。

5.解初始值问题 00'(,)()y f x y y x y =⎧⎨=⎩近似解的梯形公式是 1______k y +≈。

数值分析期末考试和答案

数值分析期末考试和答案

数值分析期末考试和答案一、单项选择题(每题2分,共20分)1. 在数值分析中,下列哪个方法用于求解线性方程组?A. 插值法B. 迭代法C. 直接法D. 拟合法答案:C2. 以下哪个数值方法是用于求解非线性方程的?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 拉格朗日插值法答案:B3. 在数值积分中,梯形法则的误差与下列哪个因素无关?A. 被积函数的二阶导数B. 积分区间的长度C. 积分区间的划分数量D. 被积函数的一阶导数答案:D4. 以下哪个数值方法是用于求解常微分方程的?A. 欧拉方法B. 牛顿迭代法C. 拉格朗日插值法D. 高斯消元法答案:A5. 在数值分析中,下列哪个方法用于求解特征值问题?A. 高斯消元法B. 幂迭代法C. 牛顿迭代法D. 梯形法则答案:B6. 以下哪个数值方法是用于求解线性最小二乘问题的?A. 高斯消元法B. 梯形法则C. 正交分解法D. 牛顿迭代法答案:C7. 在数值分析中,下列哪个方法用于求解非线性方程组?A. 高斯消元法B. 牛顿迭代法C. 线性插值法D. 欧拉方法答案:B8. 在数值分析中,下列哪个方法用于求解偏微分方程?A. 有限差分法B. 牛顿迭代法C. 线性插值法D. 梯形法则答案:A9. 在数值分析中,下列哪个方法用于求解优化问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 单纯形法答案:D10. 在数值分析中,下列哪个方法用于求解插值问题?A. 高斯消元法B. 梯形法则C. 牛顿迭代法D. 拉格朗日插值法答案:D二、填空题(每题2分,共20分)1. 在数值分析中,求解线性方程组的直接法包括______消元法和______消元法。

答案:高斯;LU2. 牛顿迭代法的收敛速度是______阶的。

答案:二3. 梯形法则的误差与被积函数的______阶导数有关。

答案:二4. 欧拉方法是一种求解______阶常微分方程的数值方法。

答案:一5. 幂迭代法是求解______特征值问题的数值方法。

数值分析期末考试题

数值分析期末考试题

数值分析期末考试题一、选择题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处的插值误差。

数值分析期末试题及答案

数值分析期末试题及答案

数值分析期末试题及答案试题一: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 \) 处插值,并求出插值多项式。

数值分析期末考卷

数值分析期末考卷

数值分析期末考卷一、选择题(每题4分,共40分)A. 插值法B. 拟合法C. 微分法D. 积分法A. 高斯消元法B. 高斯赛德尔迭代法C. 共轭梯度法D.SOR方法3. 下列哪个算法不是求解非线性方程的方法?A. 二分法B. 牛顿法C. 割线法D. 高斯消元法A. 梯形法B. 辛普森法C. 高斯积分法D. 复化求积法A. 欧拉法B. 龙格库塔法C.亚当斯法D. 高斯消元法A. 幂法B. 反幂法C. 逆迭代法D. QR算法A. 梯度下降法B. 牛顿法C. 共轭梯度法D. 高斯消元法A. 拉格朗日插值法B. 牛顿插值法C. 埃尔米特插值法D. 分段插值法A. 前向差分法B. 后向差分法C. 中心差分法D. 拉格朗日插值法A. 牛顿法B. 割线法C. 雅可比迭代法D. 高斯消元法二、填空题(每题4分,共40分)1. 数值分析的主要任务包括数值逼近、数值微积分、数值线性代数和______。

2. 在求解线性方程组时,迭代法的收敛速度与______密切相关。

3. 牛顿法的迭代公式为:x_{k+1} = x_k f(x_k)/______。

4. 在数值积分中,复化梯形公式的误差为______。

5. 求解常微分方程初值问题,龙格库塔法的阶数取决于______。

6. 矩阵特征值的雅可比方法是一种______方法。

7. 梯度下降法在求解无约束优化问题时,每次迭代的方向为______。

8. 拉格朗日插值多项式的基函数为______。

9. 数值微分中的中心差分公式具有______阶精度。

10. 在求解非线性方程组时,牛顿法的迭代公式为:x_{k+1} =x_k J(x_k)^{1}______。

三、计算题(每题10分,共60分)1. 给定数据点(1,2),(2,3),(3,5),(4,7),求经过这四个数据点的拉格朗日插值多项式。

2. 用牛顿迭代法求解方程x^3 2x 5 = 0,初始近似值为x0 = 2,计算前三次迭代结果。

数值分析期末复习-福大研究生版

数值分析期末复习-福大研究生版

数值分析期末复习题型:一、填空 二、判断 三、解答(计算) 四、证明第一章 误差与有效数字一、 有效数字1、 定义:若近似值x*的误差限是某一位的半个单位,该位到x*的第一位非零数字共有n 位,就说x*有n 位有效数字。

2、 两点理解:(1) 四舍五入的一定是有效数字(2) 绝对误差不会超过末位数字的半个单位eg. 3、 定理1(P6):若x*具有n 位有效数字,则其相对误差限为4、 考点:(1)计算有效数字位数:一个根据定义理解,一个根据定理1(P7例题3)二、 避免误差危害原则 1、 原则:(1) 避免大数吃小数(方法:从小到大相加;利用韦达定理:x1*x2= c / a )(2) 避免相近数相减(方法:有理化)eg. 或(3) 减少运算次数(方法:秦九韶算法)eg.P20习题14三、 数值运算的误差估计 1、 公式:(1) 一元函数:|ε*( f (x *))| ≈ | f ’(x *)|·|ε*(x )|或其变形公式求相对误差(两边同时除以f (x *)) eg.P19习题1、2、5(2) 多元函数(P8)eg. P8例4,P19习题4第二章 插值法一、 插值条件1、 定义:在区间[a,b]上,给定n+1个点,a ≤x 0<x 1<…<x n ≤b 的函数值yi=f(xi),求次数不超过n 的多项式P(x),使 2、 定理:满足插值条件、n+1个点、点互异、多项式次数≤n 的P(x)存在且唯一二、 拉格朗日插值及其余项1、 n 次插值基函数表达式(P26(2.8))2、 插值多项式表达式(P26(2.9))3、 插值余项(P26(2.12)):用于误差估计*(1)11102n r a ε--≤⨯;x εx εx εx ++=-+();1ln ln ln ⎪⎪⎭⎫ ⎝⎛+=-+x εx εx x cos 1-2sin 22x =n i y x P ii n ,,2,1,0)(Λ==4、 插值基函数性质(P27(2.17及2.18))eg.P28例1三、 差商(均差)及牛顿插值多项式 1、 差商性质(P30):(1) 可表示为函数值的线性组合(2) 差商的对称性:差商与节点的排列次序无关 (3) 均差与导数的关系(P31(3.5)) 2、 均差表计算及牛顿插值多项式四、埃尔米特插值(不用背公式) 两种解法:(1) 用定义做:设P 3(x)=ax 3+bx 2+cx+d ,将已知条件代入求解(4个条件:节点函数值、导数值相等各2个)(2) 牛顿法(借助差商):重节点eg.P49习题14 五、三次样条插值定义(1) 分段函数,每段都是三次多项式(2) 在拼接点上连续(一阶、二阶导数均连续) (3)考点:利用节点函数值、导数值相等进行解题第三章 函数逼近与曲线拟合一、 曲线拟合的最小二乘法解题思路:确定ϕi ,解法方程组,列方程组求系数(注意ϕi 应与系数一一对应)eg.P95习题17nj y x S j j ,,1,0,)(Λ==形如y=ae bx 解题步骤: (1) 线性化(2)重新制表(3)列法方程组求解(4)回代第四章 数值积分与数值微分一、 代数精度 1、 概念:如果某个求积公式对于次数不超过m 的多项式准确成立,但对于m+1次多项式不准确成立,则称该求积公式具有m 次代数精度 2、 计算方法:将f(x)=1,x,x 2, …x n 代入式子求解 eg.P100例1二、 插值型的求积公式求积系数定理:求积公式至少具有n 次代数精度的充要条件是:它是插值型的。

数值分析 计算方法 期末试卷3及参考答案

数值分析 计算方法 期末试卷3及参考答案

数值分析计算方法期末试卷3及参考答案数值分析计算方法期末试卷3及参考答案数值分析&lpar;计算方法&rpar;期末试卷3及参考答案一.填空题(每空3分后,共30分后)1.截断误差2.-x(x-2),x(x-1),103.14.=x-x2xk-f(xk)k2x(x5.6,5,26,9k-f'k)1.结构轻节点的差商表中:所以,要求的newton插值为:n3(x)=-5(x-1)+2(x-1)(x-2)+(x-1)(x-2)(x-3)插值余项是:r(x)=f'''(ξ)(x-1)2(x-2)或:r(x)=f[x,1,2,3,4](x-1)(x-2)(x-3)(x-4)=x3-4x2+32.(1)解:f(x)=1时,左=⎰10f(x)dx=1,右=a0+a1,左=右得:a0+a1=1f(x)=x时,左=⎰1f(x)dx=102,右=b10+a1,左=右得:b0+a1=2f(x)=x2时,左=⎰10f(x)dx=13,右=a1,左=右得:a1=13联立上述三个方程,解得:0=3,b0=6,a1=3f(x)=x3时,左=⎰1110f(x)dx=4,右=a1=3,左≠右所以,该求积公式的代数精度是2(2)解:过点0,1构造f(x)的hermite 插值h2(x),因为该求积公式代数精度为2,所以有:⎰h2(x)dx=a0h2(0)+a1h2(0)+b0h2(0)=a0f(0)+a1f(1)+b'0f(0)其求积余项为:r(f)=⎰1f(x)dx-[a0f(0)+a'01f(1)+b0f(0)]f'''(η)0f(x)dx-⎰0h2(x)dx=⎰(x-1)dx=f'''(ζ)123!⎰0x(x-1)dx=-f'''(ζ)72所以,k=-3.解:改进的euler公式是:具体内容至本题中,解的公式就是:⎰n+1=yn+0.2(3xn+2yn)=1.4yn+0.6xn⎰⎰yn+1=yn+0.1[3xn+2yn+3xn+1+2n+1]⎰y(0)=1⎰⎰n+1=yn+hf(xn,yn)⎰yn+1=yn+[f(xn,yn)+f(xn+1,n+1)]⎰⎰2代入求解得:1=1.4,y1=1.542=2.276,y2=2.48324.解:设f(x)=x3+2x-5,则f'(x)=3x2+2,牛顿迭代公式为:xk+1=xk-将x0=1.5代入上式,得f(xk)f'(xk)3xk+2xk-5=xk-23xk+232xk+53xk+2x1=1.34286,x2=1.37012,x3=1.32920,x4=1.32827,x5=1.32826x5-x4=0.00001所以,方程的近似根x5=1.328265.解,jacobi迭代公式是:⎰k+152k1k⎰x1=3-3x2-3x3⎰⎰k+13k⎰x2=-x12⎰k+1⎰x3=4-x1k+1⎰⎰gauss-seidel迭代公式是:⎰k+152k1k⎰x1=3-3x2-3x3⎰⎰k+13k+1⎰x2=-x12⎰k+1⎰x3=4-x1k+1⎰⎰(2)设其系数矩阵是a,将a分解为:a=d-l-u,其中⎰d=300⎰⎰000⎰⎰0-2020⎰⎰,l=-200⎰⎰001⎰⎰,u=00⎰⎰-100⎰⎰⎰00jacobi 运算矩阵就是:⎰⎰0-2-1⎰bj=d-1(l+u)=00⎰⎰-200⎰⎰⎰001⎰⎰⎰-100⎰⎰⎰-100⎰⎰⎰-100⎰⎰⎰gauss-seidel迭代矩阵是:b-l)-1u=300⎰⎰0-2-1⎰220⎰⎰000⎰j=(d⎰101⎰⎰⎰⎰000⎰⎰⎰20=10⎰⎰6-230⎰0-2⎰00⎰-206⎰⎰⎰00二.证明-1⎰0⎰⎰0⎰⎰=1⎰0-2302⎰02-1⎰证明:x1a0>0且xk+1=2x)⇒xk>0所以存有:xak+1=1(xkx)≥12k2即:数列xk有下界;x2kk+1=2(xk+x)≤2(xk+x)=xk由单调递减且有下界的数列极限存在可知序列xk极限存在。

数值分析报告期末考试复习题及其问题详解

数值分析报告期末考试复习题及其问题详解

数值分析期末考试复习题及其答案1. 已知325413.0,325413*2*1==X X 都有6位有效数字,求绝对误差限。

(4分)解:由已知可知,n=65.01021,0,6,10325413.0016*1=⨯==-=⨯=ε绝对误差限n k k X 2分 620*21021,6,0,10325413.0-⨯=-=-=⨯=ε绝对误差限n k k X 2分2. 已知⎢⎢⎢⎣⎡=001A 220- ⎥⎥⎥⎦⎤440求21,,A A A ∞ (6分) 解:{},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 ()A A A T max 2λ= 1分⎢⎢⎢⎣⎡=001A A T 420 ⎥⎥⎥⎦⎤-420⎢⎢⎢⎣⎡001 220- ⎥⎥⎥⎦⎤440=⎢⎢⎢⎣⎡001 080 ⎥⎥⎥⎦⎤3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A3. 设32)()(a x x f -= (6分) ① 写出f(x)=0解的Newton 迭代格式② 当a 为何值时,)(1k k x x ϕ=+ (k=0,1……)产生的序列{}k x 收敛于2解:①Newton 迭代格式为:xa x x x ax a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(22321+=+=---=-=+ϕ 3分②时迭代收敛即当222,11210)2(',665)('2<<-<-=-=a a x a x ϕϕ 3分4. 给定线性方程组Ax=b ,其中:⎢⎣⎡=13A ⎥⎦⎤22,⎥⎦⎤⎢⎣⎡-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α(k=0,1……)求解Ax=b ,问取什么实数α,可使迭代收敛(8分)解:所给迭代公式的迭代矩阵为⎥⎦⎤--⎢⎣⎡--=-=ααααα21231A I B 2分其特征方程为 0)21(2)31(=----=-αλαααλλB I 2分即,解得αλαλ41,121-=-= 2分 要使其满足题意,须使1)(<B ρ,当且仅当5.00<<α 2分5. 设方程Ax=b ,其中⎢⎢⎢⎣⎡=211A 212 ⎥⎥⎥⎦⎤-112,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=765b 试讨论解此方程的Jacobi 迭代法的收敛性,并建立Gauss-Seidel 迭代格式 (9分)解:U D L A ++=⎢⎢⎢⎣⎡--=+-=-210)(1U L D B J 202-- ⎥⎥⎥⎦⎤-012 3分0,03213=====-λλλλλJ B I 2分即10)(<=J B ρ,由此可知Jacobi 迭代收敛 1分 Gauss-Seidel 迭代格式:⎪⎩⎪⎨⎧--=--=+-=++++++)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(12276225k k k k k k k k k x x x x x x x x x (k=0,1,2,3……) 3分6. 用Doolittle 分解计算下列3个线性代数方程组:i i b Ax =(i=1,2,3)其中⎢⎢⎢⎣⎡=222A 331 ⎥⎥⎥⎦⎤421,23121,,974x b x b b ==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= (12分)解:①11b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=9741x A=⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211=LU 3分 由Ly=b1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡974 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 1分 由Ux1=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 得x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 2分 ②22b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 由Ly=b2=x1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 1分 由Ux2=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 得x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 2分 ③33b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0由Ly=b3=x2,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 1分 由Ux3=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 得x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-025.0375.0 2分7. 已知函数y=f(x)有关数据如下:要求一次数不超过3的H 插值多项式,使'11'33)(,)(y x H y x H i i == (6分)解:作重点的差分表,如下:3分21021101011001003))(](,,,[))(](,,[)](,[][)(x x x x x x x x f x x x x x x x f x x x x f x f x H --+--+-+= =-1+(x+1)-x(x+1)+2x.x(x+1)=232x x + 3分8. 有如下函数表:试计算此列表函数的差分表,并利用Newton 前插公式给出它的插值多项式 (7分)解:由已知条件可作差分表,3分i ih x x i =+=0 (i=0,1,2,3)为等距插值节点,则Newton 向前插值公式为: 033210022100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ∆---+∆--+∆-+==4+5x+x(x-1)=442++x x 4分9. 求f(x)=x 在[-1,1]上的二次最佳平方逼近多项式)(2x P ,并求出平方误差 (8分)解:令22102)(x a x a a x P ++= 2分取m=1, n=x, k=2x ,计算得: (m,m)=dx ⎰-111=0 (m,n)=dx x ⎰-11=1 (m,k)= dx x ⎰-112=0(n,k)= dx x ⎰-113=0.5 (k,k)= dx x ⎰-114=0 (m,y)= dx x ⎰-11=1(n,y)=dx x⎰-112=0 (k,y)= dx x ⎰-113=0.5得方程组:⎪⎩⎪⎨⎧==+=5.05.005.011201a a a a 3分解之得c a a c a 2,1,210-=== (c 为任意实数,且不为零)即二次最佳平方逼近多项式222)(cx x c x P -+= 1分 平方误差:32),(22222222=-=-=∑=i i i y a fp f ϕδ 2分10. 已知如下数据:用复合梯形公式,复合Simpson 公式计算⎰+=10214dx x π的近似值(保留小数点后三位) (8分)解:用复合梯形公式:)}1()]87()43()85()21()83()41()81([2)0({1618f f f f f f f f f T ++++++++==3.139 4分用复合Simpson 公式: )}1()]43()21()41([2)]87()85()83()81([4)0({2414f f f f f f f f f S ++++++++==3.142 4分11. 计算积分⎰=20sin πxdx I ,若用复合Simpson 公式要使误差不超过51021-⨯,问区间]2,0[π要分为多少等分?若改用复合梯形公式达到同样精确度,区间]2,0[π应分为多少等分? (10分)解: ①由Simpson 公式余项及x x f x x f sin )(,sin )()4(==得544)4(2041021)1()4(360)(max )4(1802)(-≤≤⨯≤=≤n x f n f R x n πππππ 2分即08.5,6654≥≥n n ,取n=6 2分即区间]2,0[π分为12等分可使误差不超过51021-⨯ 1分②对梯形公式同样1)(''max 20≤≤≤x f x π,由余项公式得51021)2(122)(-⨯≤≤n f R n ππ2分即255,2.254=≥n n 取 2分即区间]2,0[π分为510等分可使误差不超过51021-⨯ 1分12. 用改进Euler 格式求解初值问题:⎩⎨⎧==++1)1(0sin 2'y x y y y 要求取步长h 为0.1,计算y(1.1)的近似值 (保留小数点后三位)[提示:sin1=0.84,sin1.1=0.89] (6分)解:改进Euler 格式为:⎪⎩⎪⎨⎧++=+=+-++-+)],(),([2),(1111n n n n n n n n n n y x f y x f hy y y x hf y y 2分 于是有⎪⎩⎪⎨⎧+++-=+-=+-++-+-+)sin sin (05.0)sin (1.012112121n n n n n n n n n n n n n x y y x y y y y x y y y y (n=0,1,2……) 2分 由y(1)=0y =1,计算得⎪⎩⎪⎨⎧=≈=+-=-838.0)1.1(816.0)1sin 11(1.01121y y y 2分 即y(1.1)的近似值为0.83813. ][],[],,[lim ],[),,(],,[)(0'000000'x f x x f x x f x x f b a x b a C x f x x ==∈∈→证明:定义:设(4分)证明:]['],[],[],[lim ][][lim]['00000000000x f x x f x x f x x f x x x f x f x f x x x x ===--=→→故可证出 4分14. 证明:设nn RA ⨯∈,⋅为任意矩阵范数,则A A ≤)(ρ (6分)证明:设λ为A 的按模最大特征值,x 为相对应的特征向量,则有Ax=λx 1分 且λρ=)(A ,若λ是实数,则x 也是实数,得Ax x =λ 1分而x x ⋅=λλ x A x ,⋅≤⋅⋅≤λ故x A Ax 2分由于A x 0x ≤≠λ得到,两边除以 1分故A A ≤)(ρ 1分 当λ是复数时,一般来说x 也是复数,上述结论依旧成立。

数值分析期末复习题

数值分析期末复习题

一、填空题1.设真值x=983350,则其近似值y=98000的有效数字的位数 ,绝对误差为 , 相对误差为 。

2.x=0.1062,y=0.947,计算x+y 其有效数字的位数为 。

3.对f(x)=x 3+x+1,差商f[0,1,2,3]= ;f[0,1,2,3,4]= 。

4.设f(x)可微,求方程x=f(x)根的牛顿迭代法格式是 。

5.设方程x=ϕ(x)有根x *,且设ϕ(x)在含x *的区间(a,b)内可导,设x 0∈(a,b)则迭代格式x k+1=ϕ(x k )收敛的充要条件为 。

6.求解线性方程组Ax=b 的迭代格式x (k+1)=Jx (k)+f 收敛的充要条件为 。

7.⎪⎪⎭⎫ ⎝⎛=011001001001....A ,||A||∝= ,cond(A)∝= 。

8.n 次Legendre 多项式的最高次项系数为 。

9.中矩形公式:)()2()(a b b a f dx x f b a -+=⎰的代数精度为 。

10.求积公式:)1(21)0()(10f f dx x f '+≈⎰的代数精度为 。

11.在区间[1,2]上满足插值条件⎩⎨⎧==3)2(1)1(P P 的一次多项式P(x)= 。

12.设∑==n k k k n x f A f I 0)()(是函数f(x)在区间[a,b]上的插值型型求积公式,则 ∑=n k k A0= 。

13.梯形公式和改进的Euler 公式都是 阶精度的。

二、计算题1.利用矩阵的高斯消元法,解方程组⎪⎩⎪⎨⎧=++=++=++2053182521432321321321x x x x x xx x x2.设有函数值表试求各阶差商,并写出Newton 插值多项式。

3.求解超定方程组⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛43231211121x x的最小二乘解。

4.给定下列函数值表:求3次自然样条插值函数5.给定x x f =)(在x=100, 121, 144 三点处的值,试以这三点建立f(x)的二次(抛物)插值公式,利用插值公式求115的近似值并估计误差。

数值分析试题及答案汇总

数值分析试题及答案汇总

数值分析试题一、填空题(2 0×2′) 1.⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-=32,1223X A 设x =0.231是精确值x *=0.229的近似值,则x 有 2 位有效数字。

2. 若f (x )=x 7-x 3+1,则f [20,21,22,23,24,25,26,27]= 1 ,f [20,21,22,23,24,25,26,27,28]= 0 。

3. 设,‖A ‖∞=___5 ____,‖X ‖∞=__ 3_____,‖AX ‖∞≤_15_ __。

4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ,则使用该迭代函数的迭代解法一定是局部收敛的。

5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。

6. 当插值节点为等距分布时,若所求节点靠近首节点,应该选用等距节点下牛顿差商公式的 前插公式 ,若所求节点靠近尾节点,应该选用等距节点下牛顿差商公式的 后插公式 ;如果要估计结果的舍入误差,应该选用插值公式中的 拉格朗日插值公式 。

7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是:=∑=ni i x a 0)( 1 ;所以当系数a i (x )满足 a i (x )>1 ,计算时不会放大f (x i )的误差。

8. 要使20的近似值的相对误差小于0.1%,至少要取 4 位有效数字。

9. 对任意初始向量X (0)及任意向量g ,线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。

10. 由下列数据所确定的插值多项式的次数最高是 5 。

11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。

12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的,其中的残差r i = (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ,(i =0,1,…,n )。

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

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

《计篥方法P 实验报告1. 已知X ; =325413, X ; =0.325413都有6位有效数字,求绝对误差限。

(4 分)解:由已知可知,n 二6X : =0.325413x1()6* =6北一n = 0,绝对误差限^ =丄 xl0° =0.522X ; =0・325413xl0°,k=(U—〃 = —6,绝对误差限& =-xl0"62・ -2分心("刃=皿{1,8,32} = 32 1分|H|2 =732=4^23. 设/(x) = (x 2-«)3(6 分)① 写岀f (x)二0解的Newton 迭代格式②当a 为何值时,仏|=卩(忑)(k 二0,1……)产生的序列伉}收敛于、伍【值分析期末考试复习题及其答案1 02.已知4= 02 0 -2解:”州=max{l,4,8} = &分4求IKMJK (6分) 4||^||x =max{l,6,6} = 6,分皿讥如)分_1 0 0 ■ '1A TA = 0 2-20 4 40 0■ 24 二 08 0 -2 4.0 32_w :①Newton 迭代格式为:丿(忑)0(戈)=竺+丄6 6%(屛-°)' _ 5x k a , ・・ , ,6忑(X ;_G )26 6x©⑴三一曲曲以血)卜10—6/~vF<1,即-2<“ <22时迭代收敛 4・给定线性方程组Ax 二b,其中:A =3 -1用迭代公式牙=才「+a(b- Ax (k}) (k=0,1 ........... )求解 Ax 二b, 问取什么实数Q ,可使迭代收敛 (8分)-a 1 一2a.其特征方程为|刀-=八° 一3a}2(X=02分aA-(l -2<z)即,解得=l-a,22 =l-4a2分 要使其满足题意,须使p (B ) < 1,当且仅当0 vav0・52分'12 -2"丁111,b = 6 2 2 1.7.迭代法的收敛性,并建立Gauss-Seidel 迭代格式 (9分)解:A=L+D+UB, =-D~\L + U)= -1-2 -2 0 -22 -1 0|/1/ — By | = A 3= 0,/lj = = Aj = 0即p (B y ) = 0<l,由此可知Jacobi 迭代收敛1一 3a-la所给迭代公式的迭代矩阵为B = I — aA =2分试讨论解此方程的Jacobi5. 设方程Ax 二b,其中A =Gauss-Seidel 迭代格式:X 严=5-2垮)+2宅) <垮+—6-尤严—才 兀严)=7 — 2#申一 2卅Z用Doolittle 分解汁算下列3个线性代数方程组:Ax f =b, (i=l,2,3)其中"2 1 r4' A = 2 3 2 ,s = 7 2 3 4 96. 解: 上2 =兀]上3 =X2 (12分) ① Ax l =b 、x\ =A=由 Ly=bl,由 lxl=y, ②虹=b 2=LU 即y= 4791 10 0 11 12 得xl 二12 0'2 1 rT 2 32 x2= 12 3 4.11 0 o'TT1 1 oy= 1得y= 01 1 11由 Ly=b2=xl,即 (k=0,1,2, 3……)"2 1r丁「OS由Ux2=y,即 02 1 x2= 0得x2二 00 0 2③山3 =仏"2 1r'0.5'2 3 2 x3= 023 4.)0 0"'0.5 '由 Ly=b3=x2,即 1 10 y= 0得y 二 -0.51 110 .0 .'2 1r0.5 ''0.375 -由Ux3二y,即2 1 x3= -0.5得x3二 -0.250 2要求一次数不超过3的H 插值多项式,使 H 3(x /) = y r ,/73(x 1) = y I解:作重点的差分表,如下:H 3(X)= f[x {)] + /[x^x^ ](x-x Q ) + f[x 0,x^x { ](x-x ())(x-x I ) + /[xD ,x p x p x 2](x-x 0)(x-x I )2 二T+(x+l)-x(x+l)+2x. x(x+l)二2x 3 + x 27•已知 函 数 8.有如下函数表:关 数 据(6分)试讣算此列表函数的差分表,并利用Newton前插公式给出它的插值多项式(7分) 解:由已知条件可作差分表,x f = x0 + ih = i (i=0, 1, 2, 3 )为等距插值节点,则Newton向前插值公式为:N3 = f0 + 气評农 + +(—。

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

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

数值分析期末考试复习题及其答案1. 已知325413.0,325413*2*1==X X 都有6位有效数字,求绝对误差限.(4分)解:由已知可知,n=65.01021,0,6,10325413.0016*1=⨯==-=⨯=ε绝对误差限n k k X 2分 620*21021,6,0,10325413.0-⨯=-=-=⨯=ε绝对误差限n k k X 2分2. 已知⎢⎢⎢⎣⎡=001A 220- ⎥⎥⎥⎦⎤440求21,,A A A ∞ (6分)解:{},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 ()A A A T max 2λ= 1分⎢⎢⎢⎣⎡=001A A T 420 ⎥⎥⎥⎦⎤-420⎢⎢⎢⎣⎡001 220- ⎥⎥⎥⎦⎤440=⎢⎢⎢⎣⎡001 080 ⎥⎥⎥⎦⎤3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A3. 设32)()(a x x f -= (6分) ① 写出f(x )=0解的Newton 迭代格式② 当a 为何值时,)(1k k x x ϕ=+ (k=0,1……)产生的序列{}k x 收敛于2解:①Newton 迭代格式为:xa x x x ax a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(22321+=+=---=-=+ϕ 3分②时迭代收敛即当222,11210)2(',665)('2<<-<-=-=a a x a x ϕϕ 3分4. 给定线性方程组Ax=b ,其中:⎢⎣⎡=13A ⎥⎦⎤22,⎥⎦⎤⎢⎣⎡-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α(k=0,1……)求解Ax=b ,问取什么实数α,可使迭代收敛 (8分)解:所给迭代公式的迭代矩阵为⎥⎦⎤--⎢⎣⎡--=-=ααααα21231A I B 2分其特征方程为0)21(2)31(=----=-αλαααλλB I 2分即,解得αλαλ41,121-=-= 2分 要使其满足题意,须使1)(<B ρ,当且仅当5.00<<α 2分5. 设方程Ax=b,其中⎢⎢⎢⎣⎡=211A 212 ⎥⎥⎥⎦⎤-112,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=765b 试讨论解此方程的Jacobi 迭代法的收敛性,并建立Gauss —Seidel 迭代格式 (9分)解:U D L A ++=⎢⎢⎢⎣⎡--=+-=-210)(1U L D B J 202-- ⎥⎥⎥⎦⎤-012 3分0,03213=====-λλλλλJ B I 2分即10)(<=J B ρ,由此可知Jacobi 迭代收敛 1分 Gauss-Seidel 迭代格式:⎪⎩⎪⎨⎧--=--=+-=++++++)1(2)1(1)1(3)(3)1(1)1(2)(3)(2)1(12276225k k k k k k k k k x x x x x x x x x (k=0,1,2,3……) 3分6. 用Doolittle 分解计算下列3个线性代数方程组:i i b Ax =(i=1,2,3)其中⎢⎢⎢⎣⎡=222A 331 ⎥⎥⎥⎦⎤421,23121,,974x b x b b ==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡= (12分)解:①11b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=9741x A=⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211=LU 3分 由Ly=b1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡974 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 1分 由Ux1=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡234 得x1=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 2分 ②22b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 由Ly=b2=x1,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡111 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 1分 由Ux2=y,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001 得x2=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 2分③33b Ax =⎢⎢⎢⎣⎡222 331 ⎥⎥⎥⎦⎤421x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0由Ly=b3=x2,即⎢⎢⎢⎣⎡111 110 ⎥⎥⎥⎦⎤100y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡005.0 得y=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 1分 由Ux3=y ,即⎢⎢⎢⎣⎡002 021 ⎥⎥⎥⎦⎤211x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-05.05.0 得x3=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-025.0375.0 2分7. 已知函数y=f (x)有关数据如下:要求一次数不超过3的H 插值多项式,使'11'33)(,)(y x H y x H i i == (6分)解:作重点的差分表,如下:3分21021101011001003))(](,,,[))(](,,[)](,[][)(x x x x x x x x f x x x x x x x f x x x x f x f x H --+--+-+= =-1+(x+1)-x (x+1)+2x.x(x+1)=232x x + 3分8. 有如下函数表:试计算此列表函数的差分表,并利用Newton 前插公式给出它的插值多项式 (7分)解:由已知条件可作差分表,3分i ih x x i =+=0 (i=0,1,2,3)为等距插值节点,则Newton 向前插值公式为: 033210022100003!3))()((!2))((!1)()(f h x x x x x x f h x x x x f h x x f x N ∆---+∆--+∆-+==4+5x+x (x-1)=442++x x 4分9. 求f (x )=x 在[-1,1]上的二次最佳平方逼近多项式)(2x P ,并求出平方误差 (8分)解:令22102)(x a x a a x P ++= 2分取m=1, n=x , k=2x ,计算得: (m ,m)=dx ⎰-111=0 (m,n )=dx x ⎰-11=1 (m,k)=dx x ⎰-112=0(n,k )=dx x ⎰-113=0。

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