数值分析 考试题解析
数值分析习题(含答案)
第一章 绪论姓名 学号 班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字？（有效数字的计算） 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。
3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字？（有效数字的计算）解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 解：已知δ=-**xx x ，则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。
数值分析期末试题及答案
数值分析期末试题及答案试题一：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) 迭代法和直接法是数值计算中常用的两种方法。
迭代法通过反复迭代逼近解，直到满足所需精度为止；而直接法则通过一系列代数运算直接得到解。
东北大学数值分析考试题解析
数值分析提供了许多实用的算法， 这些算法可以解决各种实际问题， 如线性方程组、微分方程、积分 方程等。这些算法在科学计算、 工程仿真、数据分析等领域都有 广泛的应用。
数值分析在解决实际问题时具有 高效、精确和可靠的特点。通过 数值分析，我们可以快速地得到 问题的近似解，并且可以通过误 差分析来控制解的精度。这使得 数值分析成为解决实际问题的重 要工具。
详细描述
数值分析是一门应用广泛的学科，它通过数学方法将实际问题转 化为可计算的数学模型，并寻求高效的数值计算方法来求解这些 问题。数值分析在科学计算、工程、经济、金融等领域中发挥着 重要的作用，为实际问题的解决提供了有效的工具。
数值分析的应用领域
总结词
数值分析的应用领域非常广泛，包括科学计算、工程、经济、金融等。
非线性方程组的求解精度和速 度取决于所选择的方法和初值 条件。
非线性方程组的求解在科学计 算、工程技术和计算机图形学 等领域有广泛应用。
最优化方法
最优化方法是寻找使某个 函数达到最小或最大的参 数值的方法。
最优化方法的效率和精度 取决于所选择的算法和初 始参数值。
常用的最优化方法包括梯 度下降法、牛顿法和拟牛 顿法等。
数值分析在人工智能领域的应用
总结词
数值分析在人工智能领域的应用关键，涉及深度学习、神经 网络等领域。
详细描述
数值分析为人工智能提供了理论基础和算法支持，特别是在 深度学习和神经网络方面。通过数值分析的方法，可以优化 神经网络的参数和结构，提高人工智能的性能和准确性。
数值分析在金融领域的应用
总结词
常见的迭代法有雅可比迭代法 、高斯-赛德尔迭代法等。
牛顿法
牛顿法是一种基于泰勒级数 的迭代方法，用于求解非线 性方程的根。
数值分析题库答案(含详细解题步骤)
第 1 页/共 22 页1. 正方形的边长大约为100cm ，应怎样测量才干使面积误差不超过1cm 22. 已测得某场地长l 的值为110=*l m ，宽d 的值为80=*d m ，已知 2.0≤-*l l m,1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.3.为使π的相对误差小于0.001%,至少应取几位有效数字？4.设x的相对误差界为δ，求n x的相对误差界.5.设有3个近似数a=2.31，b=1.93，c=2.24，它们都有3位有效数字，试计算p=a+bc的误差界和相对误差界，并问p的计算结果能有几位有效数字？第 3 页/共 22 页6. 已知333487.034.0sin ,314567.032.0sin ==，请用线性插值计算3367.0sin 的值，并预计截断误差.7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并预计误差.8. 已知16243sin ,sin πππ===请用抛物插值求sin50的值,并预计误差9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x第 5 页/共 22 页10. 已知4)2(,3)1(,0)1(=-=-=f f f ， 求函数)(x f 过这3点的2次牛顿插 值多项式.11. 设x x f =)(，并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ，试用二次牛顿插值多项式计算(2.15)f 的近似值，并研究其误差12. 设],[)(b a x f 在上有四阶延续导数，试求满意条件)2,1,0()()(==i x f x P i i 及)()(11x f x P '='的插值多项式及其余项表达式.13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944⎡⎤⎢⎥⎣⎦，上的三次埃尔米特插值多项式()P x ,使它满意11()()(0,1,2),()(),i i P x f x i P x f x ''===并写出余项第 7 页/共 22 页表达式.14. 设],1,0[,23)(2∈++=x x x x f 试求)(x f 在]1,0[上关于,,1{,1)(x span x =Φ=ρ}2x 的最佳平方逼近多项式15.已知实验数据如下：用最小二乘法求形如y=a+bx2的拟合曲线，并计算均方误差.16.已知数据表如下第 9 页/共 22 页x i 1 2 3 4 5 y iωi4 4.56 8 8.5 2 1 3 1 1试用最小二乘法求多项式曲线与此数据组拟合17. .1)(},1{span ,1]41[)(的最佳平方逼近多项式中的关于上的在在求==Φ=x x x x f ρ18. 决定求积公式⎰++≈10110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A , 使其代数精度尽量高，并指出所决定的求积公式的代数精度.19. 用复化辛普森公式计算积分⎰=10dx e I x , 问区间[0,1]应分多少等分才干使截断误差不超过?10215-⨯第 11 页/共 22 页20. 利用下表中给出的数据，分离用复化梯形公式和复化辛甫生公式计算定积分dx x I ln 21⎰=的近似值（要求结果保留到小数点后六位）21. 用复化梯形公式和复化辛甫生公式计算积分⎰=6.28.1)(dx x f I ,函数)(x f 在某些节点上的值如下图：（本题共14分）22. 决定公式⎰+≈101100)()()(x f A x f A dx x f x 的系数1010,,,x x A A ,使其具有最高代数精度23. 决定求积公式⎰++≈1110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A ,使其代数精度尽量高，并指出所决定的求积公式的代数精度第 13 页/共 22 页24.用LU 分解法求解以下方程组 (10分)123123142521831520x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭25.用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛8892121514131615141321x x x26. 用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛542631531321321x x x27. 设方程组b Ax =，其中⎪⎪⎪⎭⎫⎝⎛-=220122101A ，Tb ⎪⎭⎫ ⎝⎛-=32,31,21, 已知它有解Tx ⎪⎭⎫⎝⎛-=0,31,21，若右端有小扰动61021-∞⨯=bδ，试预计由此引起的解的相对误差.第 15 页/共 22 页28. 设方程组b Ax =，其中212 1.0001A -⎛⎫= ⎪-⎝⎭，11.0001b -⎛⎫= ⎪⎝⎭,当右端向量b 有误差00.0001δ⎛⎫= ⎪⎝⎭b 时，试预计由此引起的解的相对误差(用∞范数计算)29. 给定b Ax =，其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a a a a a a A 证实：(1) 当121<<-a 时，A 对称正定，从而GS 法收敛. (2) 惟独当2121<<-a 时，J 法收敛.30. 对于线性方程组⎪⎩⎪⎨⎧-=+-=-+=+1242043 16343232121x x x x x x x ，列出求解此方程组的Jacobi 迭代格式，并判断是否收敛。
(完整)数值分析试题库与答案解析,推荐文档
模 拟 试 卷（一）一、填空题（每小题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）内。
数值分析试卷及答案
数值分析试卷及答案数值分析试卷一、选择题（共10题，每题2分，共计20分）1. 数值分析的研究内容主要包括以下哪几个方面？A. 数值计算方法B. 数值误差C. 数值软件D. 数学分析答：A、B、C2. 下列哪种方法不属于数值积分的基本方法？A. 插值法B. 微积分基本公式C. 数值微积分D. 数值积分公式答：A3. 数值积分的目的是求解什么？A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答：D4. 数值微分的目的是求解什么？A. 函数的导数B. 函数的原函数C. 函数的极值D. 函数的积分答：A5. 数值微分的基本方法有哪几种？A. 前向差分B. 后向差分C. 中心差分D. 插值法答：A、B、C6. 用数值方法求解方程的基本方法有哪几种？A. 迭代法B. 曲线拟合法C. 插值法D. 数值积分法答：A、B、C7. 用迭代法求方程的根时，当迭代结果满足何条件时可停止迭代？A. 当迭代结果开始发散B. 当迭代结果接近真实解C. 当迭代次数超过一定阈值D. 当迭代结果在一定范围内波动答：B8. 下列哪种插值方法能够确保经过所有给定数据点？A. 拉格朗日插值B. 牛顿插值C. 三次样条插值D. 二次插值答：A、B、C9. 数值解线性方程组的基本方法有哪几种？A. 直接法B. 迭代法C. 插值法D. 拟合法答：A、B10. 下列哪种方程求解方法适用于非线性方程？A. 直接法B. 迭代法C. 插值法D. 曲线拟合法答：B二、填空题（共5题，每题4分，共计20分）1. 数值积分的基本公式是_________。
答：牛顿-科特斯公式2. 数值微分的基本公式是_________。
答：中心差分公式3. 数值积分的误差分为_________误差和_________误差。
答：截断、舍入4. 用插值法求解函数值时，通常采用_________插值。
答：拉格朗日5. 数值解线性方程组的常用迭代法有_________方法和_________方法。
数值分析期末考试题及答案
数值分析期末考试题及答案一、选择题（每题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. 以下误差公式不正确的是（ 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、 已知函数表求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。
解：（1）插值基函数分别为()()()()()()()()()()1200102121()1211126x x x x x x l x x x x x x x ----===--------()()()()()()()()()()021*******()1211122x x x x x x l x x x x x x x --+-===-+---+-()()()()()()()()()()0122021111()1121213x x x x x x l x x x x x x x --+-===-+--+-故所求二次拉格朗日插值多项式为()()()()()()()()()()()2202()11131201241162314121123537623k k k L x y l x x x x x x x x x x x x x ==⎡⎤=-⨯--+⨯-+-+⨯+-⎢⎥⎣⎦=---++-=+-∑（2）一阶均差、二阶均差分别为[]()()[]()()[][][]010*********011201202303,11204,41234,,52,,126f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----===----===---故所求Newton 二次插值多项式为()()[]()[]()()()()()20010012012,,,35311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-++++-=+-例2、 设2()32f x xx =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{}span 1,x Φ=的最佳平方逼近多项式。
《数值分析》练习题及答案解析
《数值分析》练习题及答案解析第一章 绪论主要考查点：有效数字，相对误差、绝对误差定义及关系；误差分类；误差控制的基本原则；。
1. 3.142和3.141分别作为π的近似数具有（ ）和（ ）位有效数字.A ．4和3B ．3和2C ．3和4D ．4和4 答案：A2. 设 2.3149541...x *=，取5位有效数字，则所得的近似值x=___________ .答案：2.31503.若近似数2*103400.0-⨯=x 的绝对误差限为5105.0-⨯,那么近似数有几位有效数字 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
4 . 14159.3=π具有4位有效数字的近似值是多少？解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需！41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（ , ）之间的任意数，都具有4位有效数字。
第二章 非线性方程求根 主要考查点：二分法N 步后根所在的区间，及给定精度下二分的次数计算；非线性方程一般迭代格式的构造，（局部）收敛性的判断，迭代次数计算； 牛顿迭代格式构造；求收敛阶；1.用二分法求方程012=--x x 的正根，要求误差小于0.05。
（二分法）解：1)(2--=x x x f ，01)0(<-=f ，01)2(>=f ，)(x f 在[0，2]连续，故[0，2]为函数的有根区间。
"（1）计算01)1(<-=f ，故有根区间为[1，2]。
（2）计算041123)23()23(2<-=--=f ，故有根区间为]2,23[。
（3）计算0165147)47()47(2>=--=f ，故有根区间为]47,23[。
（4）计算06411813)813()813(2>=--=f ，故有根区间为]813,23[。
数值分析习题及答案解析
第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式(1.2.2)(1.2.3)则得有5位有效数字，其误差限，相对误差限有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？（1）（2）解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。
（1）（2）4.近似数x*=0.0310,是 3 位有数数字。
5.计算取，利用：式计算误差最小。
四个选项：第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。
线性插值时，用0.5及0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值误差限，故2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h 应取多少?解：用误差估计式（5.8），令因得3. 若，求和.解：由均差与导数关系于是4. 若互异，求的值，这里p≤n+1.解：，由均差对称性可知当有而当P＝n＋1时于是得5. 求证.解：解：只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式（5.17）得其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析试题库与答案解析
模 拟 试 卷（一）一、填空题（每小题3分，共30分）1．有3个不同节点的高斯求积公式的代数精度是 次的.2．设152210142-⎡⎤⎢⎥=-⎢⎥⎢⎥-⎣⎦A ，342⎛⎫ ⎪=- ⎪ ⎪⎝⎭x ，则 ∞A = .， 1x = ______.3．已知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．解初始值问题00(,)()y f x y y x y '=⎧⎨=⎩的改进的Euler 方法是 阶方法；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．1(), (),, ()n x x x 是以整数点01, ,, ,n x x x 为节点的Lagrange 插值基函数，则()n kjk k xx =∑= .9．解方程组=Ax b 的简单迭代格式(1)()k k +=+xBx 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 =+的经验公式拟合以下数据：5．用矩阵的直接三角分解法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡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.9115； 4. 1645； 5. 二； 6. (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 ++++++⎧=++⎪=+-⎨⎪=--⎩, ，，0.1543)7. 1()1()k k k k k x f x x x f x +-=-'-; 8. j x ; 9. ()1B ρ<; 10.32(4)11,()(1)(1)(2)/24(1,2)66x x x f x x x x ξξ+-+--∈-二、综合题 1．差商表：233234()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 =+-+-+-+-+ 令(2)57p =，(2)72p '=，求出a 和b. 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(∞内只有一个根s ，而且在区间（2，4）内。
数值分析——邵新慧 考试题解析
是二阶方法. (2)以此法求解y=-10y, y(0)=1时,取步长h=0.25,所得
数值解yn是否稳定?为什么? 解 (1)由于
2 K 2 f ( xn 2 h , y n 3 3 hK 1 ) f n 2 f n 2 fn h hf n x 3 y 3 2 fn 4 2 2 1 2 fn 4 2 2 fn 8 2 [ 2 h h f n 2 h f n O(h 3 ) 2 x 9 xy 9 y 9
3
3 xk a 2 a 2 x k 1 x k 或x k 1 x k x k 2 3 3 3x k 的Newton迭代格式为_______________________.
6.设l0(x),l1(x),l2(x),l3(x)是以x0,x1,x2,x3为互异节点
(x-2)3 的三次插值基函数,则 l j ( x)( x j 2) 3=____________.
H3(x)=20(x)+31(x)+52(x)+0.51(x)
令0(x)=c(x-1)2(x-2),可得0(x)=-0.5(x-1)2(x-2), 令1(x)=x(x-2)(ax+b),可得1(x)=-x(x-2), 令2(x)=cx(x-1)2,可得2(x)=0.5x(x-1)2; 令1(x)=cx(x-1)(x-2),可得1(x)=-x(x-1)(x-2),
(2)讨论这两种迭代法的收敛性.
(3)取初值x(0)=(0,0,0)T,若用Jacobi迭代法计算时,
预估误差x*-x(10) (取三位有效数字).
解
(1)Jacobi法和SOR法的迭代格式分别为
1 (k ) 1 (k ) 1 ( k 1) (k ) (k ) x1 ( x1 x 2 x 3 ) x1 4 2 4 1 ( k 1) 1 (k ) 2 ( k 1) (k ) (k ) x x ( x x x3 ) 2 2 1 2 5 5 5 1 ( k 1) 1 ( k 1) 1 ( k 1) (k ) (k ) x x ( x x x ) 3 3 1 2 3 3 6 2
数值分析报告期末考试复习题及其问题详解
数值分析期末考试复习题及其答案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 也是复数，上述结论依旧成立。
数值分析试题及答案
数值分析试题及答案一、单项选择题（每题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. 求解线性方程组时，迭代法的基本思想是从一个初始近似解出发，通过不断________来逼近精确解。
数值分析试题库与答案解析
y1
y0
h( 2
k1
k2)
2 0. 1 ( 0. 5 0. 5 7 1 4 2 9 )
2. 1 0 7 1 4 2 9
33 5 3 解 设3 5 9
5 9 17
1
d1
1 l21 l31
l21 1
d2
1 l32
l31 l32 1
d3
1
利用矩阵乘法可求得
2
5
d1 3 ， d 2 2 ， d3
， l 21 1 ， l31
7. xk 1 xk xk f (xk ) ; 8. x j ; 9. 1 f (xk)
(B) 1;
10. 1 x3
x2
1 x,
f (4) ( )( x 1)x( x 1)(x 2) / 24
6
6
( 1,2)
二、综合题
1．差商表：
1 15
20
1 15
15
20
7
1 15
22
1
42
8
2 57
30
72
2 57
＝
.
5．解初始值问题
y f ( x, y)
的改进的 Euler 方法是
y(x0) y0
阶方法；
5 x1 3 x2 0.1x3 3
6．求解线性代数方程组
2 x1 6 x2 0.7 x3 2 的高斯—塞德尔迭代公式为
，
x1 2 x2 3.5x3 1
若取 x (0) (1, 1,1) , 则 x(1)
.
7．求方程 x f ( x) 根的牛顿迭代格式是
， l 32 2
3
3
解方程组
1
y1
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
( 4)
五、(12分)试确定参数A,B,C及,使数值积分公式
( x ) x( x 1) 2 ( x 2) 4!
2 f ( x)dx Af ( ) Bf (0) Cf ( ) 有尽可能高的代数精度,并问代数精度是多少?它是否是
2
Gauss公式?
解 令公式对(x)=1,x,x2,x3,x4都精确成立,则有 4=A+B+C, 0=A-C, 16/3=A2+C2, 0=A3-C3 64/5=A4+C4 ,解得:A=C=10/9,B=16/9,=(12/5)1/2
(1)试证方程(x)=0有唯一正根; (2)构造一种收敛的迭代格式xk+1=(xk),k=0,1,2,… 计算精度为=10-2的近似根; (3)此迭代法的收敛阶是多少?说明之. 解 (1)因为0<x1时,(x)<0,x2时,(x)>0,所以(x)
仅在(1,2)内有零点,而当1<x<2时,(x)>0,故(x)单调.
3
3 xk a 2 a 2 x k 1 x k 或x k 1 x k x k 2 3 3 3x k 的Newton迭代格式为_______________________.
6.设l0(x),l1(x),l2(x),l3(x)是以x0,x1,x2,x3为互异节点
(x-2)3 的三次插值基函数,则 l j ( x)( x j 2) 3=____________.
B
k
1 B
x (1) x ( 0 )
0.7510 0.5 0.113 1 0.75
四、(13分)已知(0)=2,(1)=3,(2)=5,(1)=0.5, (1)试建立一个三次插值多项式H3(x),使满足插值条件: H3(0)=2,H3(1)=3,H3(2)=5,H3(1)=0.5; (2)设y=(x)在[0,2]上四次连续可微,试确定插值余项 R(x)=(x)-H3(x). 解 (1)由y0=2,y1=3,y2=5,y1=0.5,得
|x2-x1|=0.0035<10-2 , 故可取根的近似值x2=1.40983. (3)因为0<</2,所以() cos / 2 1 sin 0 故,此迭代法线性收敛(收敛阶为1).
三、(14分)设线性方程组 4 x1 x 2 2 x3 1 x1 5 x 2 x3 2 2 x x 6 x 3 1 2 3 (1)写出Jacobi法和SOR法的迭代格式(分量形式);
1 10h 50h 2 1 2.5 3.125 1.625 1
所以,所得数值解是不稳定的. 七、(6分)设n阶矩阵A=(aij)nn,试证实数
A n max a ij
为矩阵A的一种范数.
1 i , j n
证明
对任意n阶方阵A,B和常数,有
A n max a ij 0, 且仅当A 0时 A 0。
所以有
y ( x n 1 ) y n 1 O(h 3 )
所以此单步方法为二阶方法.
(2)此单步方法用于方程y=-10y,则有
h 20 yn 1 yn [10 yn 30( yn hyn )] [1 10h 50h 2 ] yn 4 3
当h=0.25时,有
3 2 x x 7.设S(x)= 3 2 2 x bx cx 1
j 0
0 x 1 1 x 2
是以0,1,2为节
-2 3 点的三次样条函数,则b=________c=_________.
解
由2=b+c+1,5=6+2b+c,8=12+2b,可得
二、(13分)设函数(x)=x2-sinx-1
( k 1) 1 ( k ) 1 ( k ) 1 x 2 x3 x1 4 2 4 1 (k ) 1 (k ) 2 ( k 1) x x1 x 3 2 5 5 5 1 (k ) 1 (k ) 1 ( k 1) x x1 x 2 3 3 6 2
3.向量x=(x1,x2,x3)T,试问|x1|+|2x2|+|x3|是不是一种向
是 量范数______, 而|x1|+|2x2+x3|是不是一种向量范数不是 _____.
4.求 3 a 解 只要取(x)=x3-a ,或(x)=1-x3/a. 1 5.设(x)=x3+x2-3,则差商[3,32,33,34]=_______.
H3(x)=20(x)+31(x)+52(x)+0.51(x)
令0(x)=c(x-1)2(x-2),可得0(x)=-0.5(x-1)2(x-2), 令1(x)=x(x-2)(ax+b),可得1(x)=-x(x-2), 令2(x)=cx(x-1)2,可得2(x)=0.5x(x-1)2; 令1(x)=cx(x-1)(x-2),可得1(x)=-x(x-1)(x-2),
于是有
1 2 f n f n yn 1 yn hf n h ( fn ) 2 x y 2 fn 2 fn 2 1 3 2 fn h [ 2 2 f n 2 f n ] O(h 4 ) 6 x xy y
而
1 2 1 3 y ( x n 1 ) y ( x n ) hy ( x n ) h y ( x n ) h y ( x n ) O(h 4 ) 2 6 1 2 f n f n 1 3 y n hf n h [ f n ] h y ( x n ) O(h 4 ) 2 x y 6
AB n max
1i , j n
1i , j n
a
k 1
ik kj
b n max aij (max
1i , j n 1i n
b
k 1
n
ik
)
n max aij n max bik A B
1i , k n
所以,实数‖A‖是矩阵A的范数.
考试题解析
一、填空题（每空3分,共30分)
1 2 7 25 / 7 1.设矩阵A= ,则(A)=_______,Cond(A) 1=_______. 2 3
解
1 2 2 4 7 0 A E 2 3
由于
得特征值: 1 2 3i, 2 2 3i
因此方程(x)=0有唯一正根,且在区间(1,2)内. (2)构造迭代格式: x k 1 1 sin x k k 0,1,2,... 由于|(x)|=| cos x / 2 1 sin x |<1,故此迭代法收敛.
取初值x0=1.5, 计算得x1=1.41333, x2=1.40983,由于
是二阶方法. (2)以此法求解y=-10y, y(0)=1时,取步长h=0.25,所得
数值解yn是否稳定?为什么? 解 (1)由于
2 K 2 f ( xn 2 h , y n 3 3 hK 1 ) f n 2 f n 2 fn h hf n x 3 y 3
2 fn 4 2 2 1 2 fn 4 2 2 fn 8 2 [ 2 h h f n 2 h f n ] O (h3 ) 2 x 9 xy 9 y 9
1 i , j n
A n max a ij n max a ij A
1 i , j n 1 i , j n
A B n max aij bij n max[ aij | | bij ] A B
1 i , j n
n
1 i , j n
2 ,所以‖A‖1=5,‖A-1‖1=5/7. 7 2 1 1 a a 2.设矩阵A= a 1 0 ,当a取______值时,A可以唯一分解 a 0 1 1 3 -1 又A =
为GGT,其中G为下三角矩阵.
Hale Waihona Puke 解令1 a a 1 a 1 1 1 a 2 0, a 1 0 1 2a 2 0, 得： a a 1 2 2 a 0 1
(2)因为A是严格对角占优矩阵,但不是正定矩阵,故
Jacobi法收敛,SOR法当0<1时收敛. (3)由(1)可见B=3/4,且取x(0)=(0,0,0)T,经计算可 得x(1)=(1/4,-2/5,1/2)T,于是x(1)-x(0)=1/2,所以有
x * x (10 )
(2)讨论这两种迭代法的收敛性.
(3)取初值x(0)=(0,0,0)T,若用Jacobi迭代法计算时,
预估误差x*-x(10) (取三位有效数字).
解
(1)Jacobi法和SOR法的迭代格式分别为
1 (k ) 1 (k ) 1 ( k 1) (k ) (k ) x1 ( x1 x 2 x 3 ) x1 4 2 4 1 ( k 1) 1 (k ) 2 ( k 1) (k ) (k ) x x ( x x x3 ) 2 2 1 2 5 5 5 1 ( k 1) 1 ( k 1) 1 ( k 1) (k ) (k ) x x ( x x x ) 3 3 1 2 3 3 6 2
容易验证公式对(x)=x5仍精确成立，故其代数精度为5，
是Gauss公式。 六、(12分)设初值问题
y f ( x, y ) y (a) a xb
(1)试证单步法
2 K 2 f ( xn 2 h , y K 1 f ( xn , y n ) , n 3 3 hK 1 ) h y y n 0,1,2,... n 1 n 4 ( K 1 3K 2 ) y 0
于是 H3(x)=-(x-1)2(x-2)-3x(x-2)+2.5x(x-1)2 –0.5x(x-1)(x-2) =x3-2.5x2 +2.5x+2