数值分析复习题

合集下载

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

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

数值分析期末考试复习题及其答案1.已知都有6位有效数字,求绝对误差限.(4分)解:由已知可知,n=62分2分2.已知求(6分)解:1分1分1分= 2分1分3.设(6分)①写出f(x)=0解的Newton迭代格式②当a为何值时,(k=0,1……)产生的序列收敛于解:①Newton迭代格式为: 3分② 3分4.给定线性方程组Ax=b,其中:,用迭代公式(k=0,1……)求解Ax=b,问取什么实数,可使迭代收敛(8分)解:所给迭代公式的迭代矩阵为2分其特征方程为2分即,解得2分要使其满足题意,须使,当且仅当2分5.设方程Ax=b,其中,试讨论解此方程的Jacobi迭代法的收敛性,并建立Gauss—Seidel迭代格式(9分)解:3分2分即,由此可知Jacobi迭代收敛1分Gauss-Seidel迭代格式:(k=0,1,2,3 (3)6.用Doolittle分解计算下列3个线性代数方程组:(i=1,2,3)其中,(12分)解:①A= =LU 3分由Ly=b1,即y= 得y= 1分由Ux1=y,即x1= 得x1= 2分②x2=由Ly=b2=x1,即y= 得y= 1分由Ux2=y,即x2= 得x2= 2分③x3=由Ly=b3=x2,即y= 得y= 1分由Ux3=y,即x3= 得x3= 2分7.已知函数y=f(x)有关数据如下:要求一次数不超过3的H插值多项式,使(6分)解:作重点的差分表,如下:3分=-1+(x+1)-x(x+1)+2x。

x(x+1)= 3分8.有如下函数表:试计算此列表函数的差分表,并利用Newton前插公式给出它的插值多项式(7分)解:由已知条件可作差分表,3分(i=0,1,2,3)为等距插值节点,则Newton向前插值公式为:=4+5x+x(x—1)= 4分9.求f(x)=x在[—1,1]上的二次最佳平方逼近多项式,并求出平方误差(8分)解:令2分取m=1,n=x,k=,计算得:(m,m)==0 (m,n)= =1 (m,k)= =0(n,k)= =0。

数值分析期末考试题

数值分析期末考试题

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

数值分析考试复习题

数值分析考试复习题
12/14
第五章 思考题 1. 代数插值问题的存在唯一性定理是如何叙述的 拉格朗日插值和牛顿插值方法各有何特点? 2. 拉格朗日插值和牛顿插值方法各有何特点? 3. Runge反例主要说明一个什么样的问题? 反例主要说明一个什么样的问题? 反例主要说明一个什么样的问题 第六章 思考题 多项式拟合与代数插值问题有何差异? 1. 多项式拟合与代数插值问题有何差异?拟合函数 有何特点? 有何特点? 曲线拟合的最小二乘法有何特点? 2. 曲线拟合的最小二乘法有何特点? 求一个超定方程组的最小二乘解有哪些主要方法? 3. 求一个超定方程组的最小二乘解有哪些主要方法?
11/14
第三章 思考题 高斯消元法消元过程的目标是什么? 1. 高斯消元法消元过程的目标是什么?消元过程需 用多少次乘除法? 用多少次乘除法?有何数学理论支持 解三对角方程组的消元过程有何特点? 2. 解三对角方程组的消元过程有何特点? 矩阵的范数和向量的范数有何联系, 3. 矩阵的范数和向量的范数有何联系,条件数是如 何定义的 第四章 思考题 1. 解线性方程组的迭代法有何特点?它与解方程组 解线性方程组的迭代法有何特点? 的直接法有何不同? 的直接法有何不同? 2. 解线性方程组的迭代法收敛定理对迭代产生的向 量序列的误差是如何估计的? 量序列的误差是如何估计的? 迭代法求解线性方程组的本质是什么? 3. 迭代法求解线性方程组的本质是什么?
6/14
0 yk
1 yk+1
2 yk+2
Ex16.推导左矩形求积公式 推导左矩形求积公式


b
a
f ′(η ) f ( x )dx = (b − a ) f (a ) + (b − a )2 2
a+b f ′′(η ) f ( x )dx = (b − a ) f ( )+ (b − a )3 2 24

数值分析复习试题及参考答案

数值分析复习试题及参考答案

1、设115.80,1025.621≈≈x x 均具有5位有效数字,试估计由这些数据计算21x x 的绝对误差限。

解:由有效数字的定义得,()104121-⨯=x ε,()103221-⨯=x ε()07057.00005.0115.80005.01025.621=⨯+⨯≈x x ε2、设430.56,1021.12≈≈x x均具有5位有效数字,试估计由这些数据计算21x x +的绝对误差限。

解:由有效数字的定义得,()104121-⨯=x ε,()103221-⨯=x ε0055.0)()()(2121=+=+x x x x εεε3、简答题 (1)已知12622)(256+-+-=x xxxx f ,求]1,0[f 及]6,5,4,3,2,1,0[f 。

解:由f(0)=1,f(1)=5得 []()()41011,0=-=f f f因为最高阶差商只出现在最高次,所以[]26,5,4,3,2,1,0=f(2)求积公式[])1()0(121)]1()0([21)(1f f f f dx x f '-'++≈⎰的代数精度为多少? 解:令()xx f =,则()21211021==⎰xdx x f ,右边=21,左边=右边同理令()2xx f =,()3xx f =均准确成立,()4xx f =时,左边≠右边所以,上式具有3阶精度4、求满足下表条件的Hermit 插值多项式。

x0 1)(x f -1 0 )(x f '-210解:使用重节点差商表法x y 一阶二阶 三阶 0 -1 0 -1 -2 1 0 1 3 1 010 9 6()()1236163212322---=-++--=x x x x xx x x H5、已知函数)(x f y =的数据如下:x1 2 4 -5 )(x f3 4 1 0(1)求3次Lagrange 插值多项式; (2)求3次Newton 插值多项式; (3)写出插值余项。

(完整)数值分析试题库与答案解析,推荐文档

(完整)数值分析试题库与答案解析,推荐文档

模 拟 试 卷(一)一、填空题(每小题3分,共30分)1.有3个不同节点的高斯求积公式的代数精度是 次的.2.设,,则=.,= ______.152210142-⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦A 342⎛⎫⎪=- ⎪ ⎪⎝⎭x ∞A1x3.已知y =f (x )的均差(差商),,,01214[,,]3f x x x =12315[,,] 3f x x x =23491[,,]15f x x x =, 那么均差=.0238[,,] 3f x x x =423[,,]f x x x 4.已知n =4时Newton -Cotes 求积公式的系数分别是:则,152,4516,907)4(2)4(1)4(0===C C C = .)4(3C 5.解初始值问题的改进的Euler 方法是阶方法;0(,)()y f x y y x y '=⎧⎨=⎩6.求解线性代数方程组的高斯—塞德尔迭代公式为,123123123530.13260.722 3.51x x x x x x x x x --=⎧⎪-++=⎨⎪++=⎩若取, 则.(0)(1,1,1)=- x(1)=x 7.求方程根的牛顿迭代格式是 .()x f x =8.是以整数点为节点的Lagrange 插值基函数,则01(), (),, ()n x x x 01, ,, ,n x x x =.()nk jk k x x =∑9.解方程组的简单迭代格式收敛的充要条件是.=Ax b (1)()k k +=+x Bx g 10.设,则的三次牛顿插值多项式为(-1)1,(0)0,(1)1,(2)5f f f f ====()f x ,其误差估计式为 .二、综合题(每题10分,共60分)1.求一次数不超过4次的多项式满足:,,()p x (1)15p =(1)20p '=(1)30p ''=,.(2)57p =(2)72p '=2.构造代数精度最高的形式为的求积公式,并求出10101()()(1)2xf x dx A f A f ≈+⎰其代数精度.3.用Newton 法求方程在区间内的根, 要求.2ln =-x x ) ,2(∞8110--<-kk k x x x 4.用最小二乘法求形如的经验公式拟合以下数据:2y a bx=+i x 19253038iy 19.032.349.073.35.用矩阵的直接三角分解法解方程组.⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡71735 30103421101002014321x x x x 6 试用数值积分法建立求解初值问题的如下数值求解公式0(,)(0)y f x y y y '=⎧⎨=⎩,1111(4)3n n n n n hy y f f f +-+-=+++其中.(,),1,,1i i i f f x y i n n n ==-+三、证明题(10分)设对任意的,函数的导数都存在且,对于满足x ()f x ()f x '0()m f x M '<≤≤的任意,迭代格式均收敛于的根.20Mλ<<λ1()k k k x x f x λ+=-()0f x =*x 参考答案一、填空题1.5; 2. 8, 9 ; 3.; 4. ; 5. 二; 911516456. , (0.02,0.22,0.1543)(1)()()123(1)(1)()213(1)(1)(1)312(330.1)/5(220.7)/6(12)*2/7k k k k k k k k k x x x x x x x x x ++++++⎧=++⎪=+-⎨⎪=--⎩7. ; 8. ; 9. ;1()1()k k k k k x f x x x f x +-=-'-j x ()1B ρ<10.32(4)11,()(1)(1)(2)/24(1,2)66x x x f x x x x ξξ+-+--∈-二、综合题1.差商表:11122151515575720204272152230781233234()1520(1)15(1)7(1)(1)(2)5432p x x x x x x x x x x =+-+-+-+--=++++其他方法:设233()1520(1)15(1)7(1)(1)()p x x x x x ax b =+-+-+-+-+令,,求出a 和b.(2)57p =(2)72p '=2.取,令公式准确成立,得:()1,f x x =,, , .0112A A +=011123A A +=013A =116A =时,公式左右;时,公式左, 公式右2()f x x =14=3()f x x =15=524=∴ 公式的代数精度.2=3.此方程在区间内只有一个根,而且在区间(2,4)内。

数值分析期末复习题答案

数值分析期末复习题答案

数值分析期末复习题答案一、选择题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 是一个列向量。

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

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

数值分析期末考试题及答案一、选择题(每题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. 以下误差公式不正确的是( D )A .()1212x x x x ∆-≈∆-∆B .()1212x x x x ∆+≈∆+∆C .()122112x x x x x x ∆≈∆+∆D .1122()x x x x ∆≈∆-∆ 2. 已知等距节点的插值型求积公式()()352kkk f x dx A f x =≈∑⎰,那么3kk A==∑( C )A .1 B. 2 C.3 D. 4 3.辛卜生公式的余项为( c )A .()()32880b a f η-''-B .()()312b a f η-''-C .()()()542880b a f η--D .()()()452880b a f η--4. 用紧凑格式对矩阵4222222312A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦进行的三角分解,则22r =( A ) A .1 B .12C .–1D .–25. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( D ) A .()00l x =0,()110l x = B . ()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =6. 用二分法求方程()0f x =在区间[],a b 上的根,若给定误差限ε,则计算二分次数的公式是n ≥( D )A .ln()ln 1ln 2b a ε-++ B. ln()ln 1ln 2b a ε-+-C.ln()ln 1ln 2b a ε--+ D. ln()ln 1ln 2b a ε--- 7.若用列主元消去法求解下列线性方程组,其主元必定在系数矩阵主对角线上的方程组是( B )A .123123123104025261x x x x x x x x x -+=⎧⎪-+=⎨⎪-+=-⎩ B 。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、选择题(每题2分,共20分)1. 以下哪个算法是数值分析中用于求解线性方程组的直接方法?A. 牛顿法B. 高斯消元法C. 梯度下降法D. 蒙特卡洛方法答案:B2. 插值法中,拉格朗日插值法和牛顿插值法的共同点是:A. 都是多项式插值B. 都使用差商C. 都只适用于等距节点D. 都需要预先知道所有数据点答案:A3. 在数值积分中,辛普森(Simpson)公式比梯形公式的误差:A. 更大B. 更小C. 相同D. 无法比较答案:B4. 以下哪个是数值稳定性分析中常用的方法?A. 条件数B. 收敛性C. 收敛速度D. 误差分析答案:A5. 在求解常微分方程的数值解时,欧拉方法属于:A. 单步法B. 多步法C. 隐式方法D. 显式方法答案:A6. 以下哪个是数值分析中求解非线性方程的迭代方法?A. 高斯-约当消元法B. 牛顿-拉弗森方法C. 雅可比迭代法D. 高斯-赛德尔迭代法答案:B7. 线性插值公式中,如果给定两个点\( (x_0, y_0) \)和\( (x_1, y_1) \),插值多项式是:A. \( y = y_0 + \frac{y_1 - y_0}{x_1 - x_0}(x - x_0) \)B. \( y = y_0 + \frac{y_1 - y_0}{x_0 - x_1}(x - x_0) \)C. \( y = y_0 + \frac{x - x_0}{x_1 - x_0}(y_1 - y_0) \)D. \( y = y_1 + \frac{x_1 - x}{x_1 - x_0}(y_0 - y_1) \)答案:C8. 以下哪个是数值分析中用于求解特征值问题的算法?A. 幂法B. 共轭梯度法C. 牛顿法D. 欧拉法答案:A9. 在数值微分中,使用有限差分法来近似导数时,中心差分法的误差:A. 与步长成正比B. 与步长的平方成正比C. 与步长的立方成正比D. 与步长的四次方成正比答案:B10. 以下哪个是数值分析中用于求解线性最小二乘问题的算法?A. 梯度下降法B. 牛顿法C. 奇异值分解法D. 共轭梯度法答案:C二、简答题(每题10分,共30分)1. 简述数值分析中病态问题的特点及其对算法的影响。

数值分析复习题

数值分析复习题

数值分析复习题⼀、填空1.近似数x*=0.4231关于真值x=0.4229有位有效数字.2.设f(x)可微,则求⽅程x 2=f(x)根的⽜顿迭代格式为 .3.对f(x)=x 3+3x 2-x+5,差商f[0,1,2,3,4]= .4.⽅阵A 的谱半径是指 .5.求积分?ba dx x f )(的近似值,其⾟⼘⽣公式为 .⼆、已知观测数据(1,-5),(2,0),(4,5),(5,6),试⽤最⼩⼆乘法求形如xb ax x +=)(?的经验公式。

(10分)三、求⼀个次数不⾼于4的多项式p 4(x),满⾜下列插值条件 x 0 1 2f(x) 0 1 1)(x f '0 1四、写出计算线性⽅程组=+-=+-=+-272135223121321x x x x x x x 的⾼斯⼀赛德尔迭代格式,并分析此格式的敛散性.五、⽤预估⼀校正法求初值问题=≤≤-='1)0(102y x y x y y在x=0.2处的数值解,步长取h=0.1。

(要求保留⼩数点后4位)六、把区间分成两等份,⽤复化⾟⼘⽣公式计算dxx+1七、在求⾮线性f(x)=0根的近似值时,论证简单迭代法⼀般为线性收敛,⽽⽜顿迭代法为平⽅收敛.⼀填空1.近似数x*=0.4231关于真值x=0.4229有位有效数字.2.设643()35f x x x x =-+-,则差商[0,1,2,3,4,5,6]f = 3.求积分()ba f x dx ?的近似值,其复化梯形公式为4.5点⾼斯求积公式,其代数精度为5.设f(x)可微,则求⽅程x 2=f(x)根的近似值的⽜顿迭代格式为 6.利⽤⼆分法求()0f x =在[,]a b 上根的近似值,误差限为 7.⽅阵A 的谱半径是指 8.矩阵A 的条件数是指 9.能⽤⾼斯消元法求解A x b =的充要条件是 10.设215314278A -??=,则1||||A = ⼆给定线性⽅程组1231232231242122316x x x x x x x x x -++=??-++=??++=? 1. ⽤列主元消元法求解所给线性⽅程组。

数值分析考试题

数值分析考试题

数值分析考试题一、选择题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. 给定函数f(x) = sin(x),在区间[0, π]上使用梯形法则计算积分的近似值,取4个等分点。

2. 设线性方程组如下:\[\begin{cases}2x + y + z = 6 \\x + 2y + 4z = 14 \\3x + y + 2z = 10\end{cases}\]使用高斯消元法求解该方程组的解。

数值分析复习题

数值分析复习题

一、判断题1. 区间[a,b]上,若f(a)f(b)<0,则方程f(x)=0在[a,b]内一定有实根。

2. 22/7作为π=3.1415926……近似值,它有3位有效数字。

3. 设P(x)和Q(x)都是n 次多项式,如果在n +1 个不同的节点x i 上都有P(x i )=Q(x i ),则P(x)≡Q(x) 。

4. 取节点01231, 0, 2 ,4x x x x =-===作2()f x x =的插值多项式()p x ,则()p x 次数为2,插值基函数的次数为3。

5. 插值多项式严格通过所有的节点(x i ,y i )。

6. 若k<=n ,P(x)和Q(x)分别是 x k的通过n +1 个不同的节点的牛顿插值多项式和拉格朗日插值多项式则P(x)≡Q(x)≡x k。

7. 插值多项式次数越高,逼近效果越好。

8. 任何一组互异数据,逼近它们的多项式插值函数仅有一个。

9. 插值多项式次数与拟合曲线都严格通过所给定的数据点。

10. 求积公式:⎰30)(dx x f ≈。

f f f f 是插值型的)]3()2(3)1(3)0([83+++11. 牛顿-科特斯求积公式中的求积节点是等分的。

12. 牛顿法求方程ƒ(x)=0的单根, 在ƒ(x)可导的情况下, 至少二阶收敛。

13. 高斯型求积公式是插值型的。

14. 一阶亚当姆斯格式是单步法。

15. 显式的亚当姆斯公式:+-=+-()n n n n h y y f f 1132是单步法。

16. 求初值问题数值解的四阶亚当姆斯公式是多步法。

17. 如果有一常微分方程数值解法的局部截断误差3111()()n n n T y x y O h +++=-=,则该方法是3阶的。

18. 用一般迭代法求方程()0f x =的根,如其迭代过程()1k k x x ϕ+=发散,则方程()0f x = 的无解。

19. 牛顿法求方程ƒ(x)=0的根, 在ƒ(x)可导的情况下, 至少二阶收敛。

数值分析试题及答案

数值分析试题及答案

数值分析试题及答案一、单项选择题(每题3分,共30分)1. 线性代数中,矩阵A的逆矩阵记作()。

A. A^TB. A^-1C. A^+D. A*答案:B2. 插值法中,拉格朗日插值多项式的基函数是()。

A. 多项式B. 指数函数C. 正弦函数D. 余弦函数答案:A3. 在数值积分中,梯形规则的误差是()阶的。

A. O(h^2)B. O(h^3)C. O(h)D. O(1/h)答案:A4. 求解线性方程组时,高斯消元法的基本操作不包括()。

A. 行交换B. 行乘以非零常数C. 行加行D. 行除以非零常数答案:D5. 非线性方程f(x)=0的根的迭代法中,收敛的必要条件是()。

A. f'(x)≠0B. f'(x)=0C. |f'(x)|<1D. |f'(x)|>1答案:C6. 利用牛顿法求解非线性方程的根时,需要计算()。

A. 函数值B. 函数值和导数值C. 函数值和二阶导数值D. 函数值、一阶导数值和二阶导数值答案:B7. 矩阵的特征值和特征向量是()问题中的重要概念。

A. 线性方程组B. 特征值问题C. 线性规划D. 非线性方程组答案:B8. 在数值分析中,条件数是衡量矩阵()的量。

A. 稳定性B. 可逆性C. 正交性D. 稀疏性答案:A9. 利用龙格现象说明,高阶插值多项式在区间端点附近可能产生()。

A. 振荡B. 收敛C. 稳定D. 单调答案:A10. 雅可比迭代法和高斯-塞德尔迭代法都是求解线性方程组的()方法。

A. 直接B. 迭代C. 精确D. 近似答案:B二、填空题(每题4分,共20分)11. 线性代数中,矩阵A的行列式记作________。

答案:det(A) 或 |A|12. 插值法中,牛顿插值多项式的基函数是________。

答案:差商13. 在数值积分中,辛普森规则的误差是________阶的。

答案:O(h^4)14. 求解线性方程组时,迭代法的基本思想是从一个初始近似解出发,通过不断________来逼近精确解。

《数值分析》复习题(14)

《数值分析》复习题(14)

《数值分析》复习题一、填空题1. 已知近似数 1.28y *=-,则其绝对误差限为 -0.005 ,相对误差限是 0.39% 。

2. 测量一支铅笔长是16cm , 那么测量的绝对误差限是 0.5cm ,测量的相对误差限是3.125% 。

3. 度量一根杆子长250厘米,则其绝对误差限为 0.5cm ,相对误差限是 0.20% 。

4. 在数值计算中,当a1/(√(a+1) +√a ) 5. 在数值计算中,计算356-应变成3561+来计算。

6. 在数值计算中,计算1cos3-应变为2)5.1(sin 2⨯来计算。

7. 若543()2792100f x x x x x =-+-+,则12345[1,4,4,4,4,4]f =____2__________,123456[1,3,3,3,3,3,3]f = 0 。

8. 函数()f x 关于三个节点012,,x x x 的拉格朗日二次插值多项式为 f(x)=f(x0)[(x-x1)(x-x2)/(x0-x1)(x0-x2)] ,9. 当()f x x =时,(,)n B f x =∑f (k/n )Pk(x)=x 。

10. 代数式222236()66x xR x x x +=++ ______________,323222122()23x x R x x x ++=++ __________________. 11. 已知方程组123123123103127322115x x x x x x x x x --=-⎧⎪-++=⎨⎪+-=-⎩,那么收敛的Jacobi 迭代格式为:,收敛的G S -迭代格式为:收敛理由是方程组的系数矩阵为严格对角占优阵12. 已知线性方程组1233111193234184x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦,那么收敛的Jacobi 迭代格式:收敛的G-S 迭代格式: 。

收敛理由是 严格对角占优矩阵 ,13. 求积公式0()nn kk k I Af x ==∑至少有n 次代数精度的充要条件是__严格对角占优矩阵 __________________;当n 是偶数时,牛顿-柯特斯公式()0()()nn n kk k I b a Cf x ==-∑至少有___n+1________次代数精度;高斯求积公式()()()nbk k ak f x x dx A f x ρ=≈∑⎰至少有___2n+1_______次代数精度。

(完整word版)数值分析复习题及答案

(完整word版)数值分析复习题及答案

数值分析复习题一、选择题1.3.142和3.141分别作为二的近似数具有()和()位有效数字A . 4 和 3B . 3 和 2C . 3 和 4D . 4 和 421 r 「2 11 f xdx『1 A f (-)-f (2) 2.已知求积公式6 3 6,则 A =() 11 1 2A .6B .32D . 33•通过点X 0,y 。

,x 1,y 1的拉格朗日插值基函数l ox ,h x 满足()4.设求方程f X =o的根的牛顿法收敛,则它具有( )敛速。

A .超线性B .平方C .线性D .三次X j 2x 2 x 3 = o “ 2% +2x 2 +3x 3 =35.用列主元消元法解线性方程组 厂X | — 3x 2 = 2 作第一次消元后得到的第3个方程 -X 2 X 3 =2 -2X 2 1.5X 3 二 3.5 C .—2x 2 x 3 = 3x 2 - o.5x 3 二、填空1.设-2.3149541...,取5位有效数字,则所得的近似值 x=f X 1,X 2 =2.设一阶差商 f X 2 - f X x 2 一為 1 -4 2-1 f X 1,X 2,X 3= ___________f X 2,X 3=f X 3 - f X 2X 3 -X 2().=-1.56-1 _ 5 4-2 2I 。

* )= 0, h(X 1 )=0Io (X 。

)= 0,)=1C . I o (X o )= 1 h (X 1 ) = 1Io(Xo ) = 1 h(X 1)=18、若线性代数方程组 AX=b 的系数矩阵A 为严格对角占优阵,则雅可比迭代和高斯 9、解常微分方程初值问题的欧拉(Euler )方法的局部截断误差为y=10 +丄 + 2 310、为了使计算 x-1 (x-1) (x-1)的乘除法运算次数尽量的少,应将表达式改写12•—阶均差f x ),x1-13.已知n =3时,科茨系数=8,C/心=8,那么C 33 =18.设 X=(2,一3,7)T,则 ||X|1 厂3•设 X =(2, -3,-1),则 ||X ||2 二,l|X Id4 •求方程x 2-X-仁5" 的近似根,用迭代公式 x 「x425,取初始值x=1 ,那么x1二5 •解初始值问题 y 、f(x,y)y (x0)= y °近似解的梯形公式是yk 16、 一5 1丿,则A 的谱半径Q(A)= 7、设f(x) =3x +5, X k=kh, k =0,1,2,...,则f 1人,焉 1,X n.2】 =-塞德尔迭代都11•设 X=(2,3T T ,则 ||X|"l|X||2 =14.因为方程f x =^4 2 =0在区间1,2I 上满足,所以f x =0在区间内有根。

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

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

数值分析期末考试复习题及其答案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。

数值分析复习题及答案

数值分析复习题及答案

数值分析复习题P56 16.1716.试就f (x )=2x^3+5求商f (1,2,3,4),f (1,2,3,4,5)的值 解:根据n 阶差商定义递推展开式01102110),...,,(),...,,(),...,,(x x x x x f x x x f x x x f n n n n --=-原函数f (x )=2x^3+5,;令0x =1,1x =2,2x =3,...3)3,2,1()4,3,2()4,3,2,1(x x f f f --==2(n=3)同理对f (12345)进行递推可计算得出等于017.给定函数f (x )x^3—4x ,试建立关于节点Xi=i+1(i=0,1……,5)的差商表,并列出关于节点x0,x1,x2,x3的插值多项式p (x ) 解:因为x x x f i i x i 4)()5...,1,0(13-==+= 列出差商表如下: i x i y 一阶差商 二阶差商 三阶差商 四阶差商 五阶差商 1 -3 2 0 3 3 15 15 6 4 48 33 9 1 5 105 57 12 1 0 619287151插值多项式:p(x)=-3+3(x-1)+6(x-1)(x-2)+(x-1)(x-2)(x-3) P94 2.3.42.试判断下列求积公式的代数进度:)(41)31(43)(1x f f dx x f +≈⎰解:当f (x )=1时,左边=1,右边3/4+1/4=1,左边=右边当f (x )=x 时,左边=1/2,右边3/4X1/3+1/4X1=1/2,左边=右边 当f (x )=x^2时,左边=1/3,右边3/4x1/9+1/4x1=1/3,右边=右边当f (x )=x^3时,左边=1/6,右边3/4x1/27+1/4x1=5/18,左边于右边不相等根据代数精度定义求积公式对于mx 次多项式成立,对1+m x不成立则该求积公式具有m次代数精度,所以该求积公式的代数精度为2 次。

数值分析复习题

数值分析复习题

γ1 β2
α3
1
γ2 β3
α4
-5-
= LU γ3 β4
由矩阵乘法得
β1 = 2, γ 1 = 1, γ 2 = 1, γ 3 = 1, α 2 = , β 2 = , α 3 = , β 3 = , α4 =
从而
2 1 5 1 1 2 2 3 ,U = 1 1 5 5 10 7 1 − 3 3 3 7 14 (2)解方程 LY = b 得 Y = (1, , , − ) Τ 2 5 3 1 1 2 L=
解: (1)设有分解
1 0 1 0 0 2 0 1 l 1 0 1 = 21 2 4 3 l 31 1 0 3 l 41 1 l32 l 42 1 l 43
由矩阵乘法可逐行、逐列分别求出 u ij 和 lij
1 l 21 l31 l 41 1 1 0 0 1 u 22 = , 1 2 1 1 0 1 0 1 2 u 23 u 33 0 1 0 2 0 u 24 = 1 0 1 u 34 2 1 u 44 2
3 3 5 1 3 5 9 = l 21 5 9 17 l 31 1 l32 d1 1 d2 d3 1 l 21 1 l31 l32 1
由矩阵乘法可逐行求得 2 5 d1 = 3, d 2 = 2, d 3 = , l 21 = 1, l31 = , l 32 = 2 3 3 (2)解方程组
k →∞
三、各章计算证明题 第二章 代数方程组的直接解法
1 0 1.用杜利特分解法解方程组 1 0 0 2 0 x1 5 1 0 1 x 2 3 = 2 4 3 x3 17 1 0 3 x 4 7 1 1 0 u 22 2 u 23 u 33 0 u 24 u 34 u 44
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数值方法复习题 第一章 绪论
例 1. 已知数 x=2.718281828...,取近似值 x*=2.7182,x*相对于 x 具有几位有效 数字? 点评;考查的有效数字的概念。 e* x x* 2.718281828 2.7182 0.00008182 解; 1 1 0.0005 103 101 4 2 2 故有四位有效数字。 例 2.近似数 x* 0.01999 关于真值 x* 0.02000 有几位有效数字 e* x x* 0.01999 0.02000 0.00001 解: 1 1 0.00005 104 1013 2 2 故有三位有效数字。 例 3. n x* 的相对误差为 x* 的相对误差的 倍。
xk 0.88 0.88469 0.88468 0.88468
点后第四位。 解:因为 f (0) 1 0 , f (1) 3 0 f ( x) 3 x 2 4 x 1 0, x [0,1] 所以方程在[0,1]内有唯一根。 又 f ''( x) 6 x 4 0, x [0,1] 取 x0 0.4 时,虽然 f ( x0 ) f ''( x0 ) 0 ,但是迭代一步后有 f ( x1 ) f ''( x1 ) 0 ,从此以
解:由于系数矩阵按行严格对角占优,故雅克比迭代法和高斯――赛得尔迭代法 都收敛。 1)雅克比迭代法的迭代格式为:
1 (k ) 1 (k ) 7 ( k 1) x2 x3 x1 9 9 9 7 ( k 1) 1 ( k ) x2 x1 8 8 8 ( k 1) 1 ( k ) x3 9 x1 9
1 .显然 ( x ) 在 [0 , 1] 上满足 2 x
封闭性:当 0 x 1 时, 0 ( x) 1 ;
压缩性:当 0 x 1 时, ( x)
1 1 1. 2 (2 x ) 2
故对任意初值 x0 [0 , 1] ,迭代格式 xk 1 不动点 x . (3)经计算得: x 0.414
1 , k 0, 1, 2, 必收敛到 [0 , 1] 上的惟一 2 xk
(请自己补充计算过程)
例7
当常数 c 取合适的值时,两条抛物线 y x 2 x c 与 y 2 x 就在某点相
切,试取 x0 0.3 ,用牛顿迭代法求切点横坐标。误差小于 10 4 时迭代结束。 解 两条曲线求导
0.07 0.005 0.00002
例 6:设 x
1 2 2 1 1 2
.试写出求 x 的迭代格式,讨论该格式的收敛性,并由该迭代
格式求 x 的近似值,使迭代误差不超过 0.001 . 解: (1)可建立迭代格式
xk 1
(2)收敛性讨论:设 ( x )
1 , k 0, 1, 2, 2 xk
y' 2 x 1 和 y' x

1 2
1 2
切点横坐标一定满足 2 x 1 = x 将等式变形为

f ( x) 4 x 3 4 x 2 x 1
然后仿照例 5 构造牛顿迭代公式并求解。 牛顿迭代法 结果为 0.34781 例 8.用迭代法求解方程 f ( x) e x x 2 0 的根的近似值,要求 xn xn 1 102 。 解: (1)验证区间 [0, 0.7] 隔根区间:
取 x (0) (0, 0, 0)T ,迭代结果如下:
k 0 1 2 3 4 5 x1( k ) 0 0.7779 0.9738 0.9942 0.9993 0.9998 x2 ( k ) 0 0.8750 0.9722 0.9967 0.9993 0.9999 x3( k ) 0 0.8889 0.9753 0.9971 0.9994 0.9999
故满足收敛条件,取 x0 = 0.25,计算得 x1 = (0.25) = 0.3385416 x2 = (x1) = 0.3462668 x3 = (x2) =0.3471725 x4 = (x3) =0.3472814 x5 = (x4) =0.3472945
x6 = (x5) =0.3472961 x7 = (x6) =0.3472963 取近似根为 x* = 0.347296 例 4: 用牛顿迭代法建立求平方根 c (c >0)的迭代公式,并用以上公式求
x 1 3 x 1 2 x 3 3 x 4 7
1 0 U 0 0 1 2 3 2 1 2 0 1 1 0 0 2
0 0 0 0 , 1 0 1 1
0.78265 解:设 f ( x ) x 2 c , (x >0)则 c 就是 f (x) =0 的正根。由为 f ’ (x) = 2x,所以得 迭代公式 x k 1 x k xk2 c 2 xk (2.6)

x k 1
1 c xk 2 xk
由于 x >0 时,f ’ (x) >0,且 f (x) > 0,根据定理 3 知:取任意初值 x0 c ,所确 定的迭代序列{xk}必收敛于 c 。 取初值 x = 0.88,计算结果见表 k 0 1 2 3 故可取 0.78265 0.88468
3 2 例 5:用牛顿迭代法求解 f ( x) x 2 x x 1 在 0.4 附近的根。要求精确到小数
由 Ly b 由 T 得: x (1, 0, 1, 0)T
例 3:用雅克比迭代法和高斯――赛得尔迭代法解线性方程组,要求先说明收敛 性,若收敛,取 x (0) (0, 0, 0)T ,迭代至 x ( k 1) x ( k )

103
9 1 1 x1 7 0 x2 7 1 8 1 0 9 x3 8
由 Ly b 由 Ux= y
得:y = (14, -10, -72)T 得: x (1, 2, 3)T
例 2:用矩阵三角分解法解方程组
1 1 0 2 1 1 2 2
解:由 A LU 得:
1 0 0 1 L 1 1 2 0
2 3 1 2 2 2 5 9
x f ( x) 0
根据牛顿迭代格式
x k 1 xk
f ( xk ) f ' ( xk ) xk f ( xk ) 1 f '( xk )
(k 0, 1, 2,)
则相应的得到
xk 1 xk
(k 0, 1, 2,)
例 2:求方程 f ( x) x 10 x 2 0 的一个根 解:因为 f (0) = 1>0 f (1) = -7 <0,知方程在[0, 1]中必有一实根,现将原方 程改为同解方程
n f 解:根据误差传播公式 er ( y * ) ( )* er ( xi* ) xi* i 1 xi
y*
则有 er ( n x *) ( n x *)' er ( x* ) x* / n x * 1 n
第二章 例 1.设 f ( x ) 可微,求 x f ( x ) 根的牛顿迭代公式----。 解;化简得到
此时已有: x ( k 1) x ( k )

103 ,故取 x* (0.9998, 0.9999, 0.9999)T
2)高斯――赛得尔迭代法:
1 (k ) 1 (k ) 7 ( k 1) x2 x3 x1 9 9 9 7 ( k 1) 1 ( k 1) x2 x1 8 8 8 ( k 1) 1 ( k 1) x3 9 x1 9
取 x (0) (0, 0, 0)T ,迭代结果如下: k 0 1 2 3 4 x1( k ) 0 0.7778 0.9942 0.9998 1.0000 x2 (k ) 0 0.9722 0.9993 1.0000 1.0000 x3( k ) 0 0.9753 0.9994 1.0000 1.0000
后迭代法收敛。 x n 1 牛顿迭代公式为: 列表计算 n 0 1 2 3 得根的近似值为 0.4656。 xn 0.4 0.47013 0.46559 0.46557
x n x n 1
3 2 f ( xn ) xn 2 xn xn 1 xn xn 2 f ( x n ) 3 xn 4 xn 1
取初始值 x0 = 1,可逐次算得 x1 = 0.4771 x2 = 0.3939 … x6 = 0.3758 x7 =0.3758 例 3:用迭代法求方程 x 3 3x 1 0 在[0, 0.5]内的根,精确到 10-5。 解:将方程变形 1 x ( x 3 1) 3 1 即 ( x ) ( x 3 1) , 3
第三、四章 例 1:用直接三角分解法解
1 2 3 x1 14 2 5 2 x 2 18 3 1 5 x 20 3
解:由 A LU 得:
1 L 2 1 , 3 5 1 1 2 3 U 1 4 24
10 x x 2
由此得迭代格式 x k 1 lg( x k 2)
x lg( x 2)
收敛性判断;当 x (0,1) 时, ( x) lg( x 2) (0,1) ,且由于
'( x)
1 1 0.2171 1 , 故迭代格式收敛 ( x 2) ln10 2 ln10
易知 ( x) ln(2 x), ( x)
1 x2 ,
则迭代函数在区间 [0, 0.7] 一阶导数存在;迭代函数在区间上单调递减,因此
相关文档
最新文档