数值分析考试题解析
数值分析试题及答案
数值分析试题及答案一、选择题1. 下列哪个方法不适合用于求解非线性方程的根？A. 二分法B. 牛顿法C. 弦截法D. 正割法2. 当使用二分法求解非线性方程的根时，需要满足的条件是：A. 函数f(x)在区间[a, b]上连续B. 函数f(x)在区间[a, b]上单调递增C. 函数f(x)在区间[a, b]上存在根D. 函数f(x)在区间[a, b]上可导3. 数值积分是通过将定积分转化为求和的方法来近似计算积分值的过程。
下列哪个方法是常用的数值积分方法？A. 矩形法则B. 辛普森规则C. 梯形规则D. 高斯-勒让德法则4. 龙格-库塔法是常用于求解常微分方程的数值解法。
以下哪个选项是描述龙格-库塔法的特点？A. 该方法是一种多步法B. 该方法是一种多项式插值法C. 该方法是一种单步法D. 该方法是一种数值积分法5. 用有限差分法求解偏微分方程时，通常需要进行网格剖分。
以下哪个选项是常用的网格剖分方法？A. 多边形剖分法B. 三角剖分法C. 矩形剖分法D. 圆形剖分法二、解答题1. 将函数f(x) = e^x 在区间[0, 1]上用复化梯形规则进行数值积分，分为6个子区间，求得的近似积分值为多少？解：将区间[0, 1]等分为6个子区间，每个子区间的长度为h = (1-0)/6 = 1/6。
根据复化梯形规则的公式，近似积分值为：I ≈ (1/2) * h * [f(0) + 2f(1/6) + 2f(2/6) + 2f(3/6) + 2f(4/6) + 2f(5/6) +f(1)]≈ (1/2) * (1/6) * [e^0 + 2e^(1/6) + 2e^(2/6) + 2e^(3/6) + 2e^(4/6) +2e^(5/6) + e^1]2. 使用二分法求解方程 x^3 - 3x + 1 = 0 在区间[1, 2]上的根。
要求精确到小数点后三位。
解：首先需要判断方程在区间[1, 2]上是否存在根。
数值分析习题和答案解析(最新整理)
(1)
要使
应满足（）.
(2) 已知方程组
，则解此方程组的
Jacobi 迭代法是否收敛（）.它的渐近收敛速度 R(B)=
公式（6.13）直接计算即可。
对
，取 n=8,在分点处计算 f(x)的值构造函数表。
按式（6.11）求出
，按式（6.13）求得
，
积分
2. 用 Simpson 公式求积分 ，并估计误差 解：直接用 Simpson 公式（6.7）得
由（6.8）式估计误差，因
，故
3. 确定下列求积公式中的待定参数，使其代数精确度尽量 高，并指明求积公式所具有的代数精确度.
专业知识分享
WORD 格式
编辑整理
11. 填空题
(1) 满 足 条 件
的插值多项式
p(x)=( ).
(2)
,则 f［1,2,3,4］=( )，f［1,2,3,4,5］
=( ).
(3) 设
为互异节点， 为对应的四次插值基函
数，则
＝( )，
＝( ).
(4) 设
是区间［0,1］上权函数为 ρ(x)=x 的最
WORD 格式
编辑整理
误差估计由公式（5.19）得
这里 仍为 0.565 8． 求 一 个 次 数 不 高 于 四 次 的 多 项 式 p(x),使 它 满 足
解：这种题目可以有很多方法去做，但应以简单为宜。此处
可先造 使它满足
，显然 p(x)=x2(2-x)+Ax2(x-1)2
，再令
由 p(2)=1 求出 A＝ ，于是
5.计算
取 ，利用 ：
式计算误差最小。
四个选项： 第二、三章 插值与函数逼近
习题二、三
数值分析试题与答案
一. 填空题（本大题共4小题，每小题4分，共16分）1.设有节点012,,x x x ，其对应的函数()y f x =的值分别为012,,y y y ，则二次拉格朗日插值基函数0()l x 为 。
2.设()2f x x =，则()f x 关于节点0120,1,3x x x ===的二阶向前差分为 。
3.设110111011A -⎡⎤⎢⎥=--⎢⎥⎢⎥-⎣⎦，233x ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则1A ＝ ，1x = 。
4. 1n +个节点的高斯求积公式的代数精确度为 。
二．简答题（本大题共3小题，每小题8分，共24分）1. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？2. 什么是不动点迭代法？()x ϕ满足什么条件才能保证不动点存在和不动点迭代序列收敛于()x ϕ的不动点？3. 设n 阶矩阵A 具有n 个特征值且满足123n λλλλ>≥≥≥，请简单说明求解矩阵A 的主特征值和特征向量的算法及流程。
三．求一个次数不高于3的多项式()3P x ，满足下列插值条件：i x 1 2 3 i y 2 4 12 i y '3并估计误差。
（10分）四．试用1,2,4n =的牛顿-科特斯求积公式计算定积分1011I dx x=+⎰。
（10分） 五．用Newton 法求()cos 0f x x x =-=的近似解。
（10分） 六．试用Doolittle 分解法求解方程组：12325610413191963630x x x -⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦⎣⎦ （10分）七．请写出雅可比迭代法求解线性方程组123123123202324812231530x x x x x x x x x ++=⎧⎪++=⎨⎪-+=⎩ 的迭代格式，并判断其是否收敛？（10分）八．就初值问题0(0)y yy y λ'=⎧⎨=⎩考察欧拉显式格式的收敛性。
（10分）《数值分析》（A ）卷标准答案（2009－2010－1）一． 填空题（每小题3分，共12分） 1. ()1200102()()()()x x x x l x x x x x --=--； 2.7；3. 3，8；4. 2n+1。
东北大学数值分析考试题解析
数值分析提供了许多实用的算法， 这些算法可以解决各种实际问题， 如线性方程组、微分方程、积分 方程等。这些算法在科学计算、 工程仿真、数据分析等领域都有 广泛的应用。
数值分析在解决实际问题时具有 高效、精确和可靠的特点。通过 数值分析，我们可以快速地得到 问题的近似解，并且可以通过误 差分析来控制解的精度。这使得 数值分析成为解决实际问题的重 要工具。
详细描述
数值分析是一门应用广泛的学科，它通过数学方法将实际问题转 化为可计算的数学模型，并寻求高效的数值计算方法来求解这些 问题。数值分析在科学计算、工程、经济、金融等领域中发挥着 重要的作用，为实际问题的解决提供了有效的工具。
数值分析的应用领域
总结词
数值分析的应用领域非常广泛，包括科学计算、工程、经济、金融等。
非线性方程组的求解精度和速 度取决于所选择的方法和初值 条件。
非线性方程组的求解在科学计 算、工程技术和计算机图形学 等领域有广泛应用。
最优化方法
最优化方法是寻找使某个 函数达到最小或最大的参 数值的方法。
最优化方法的效率和精度 取决于所选择的算法和初 始参数值。
常用的最优化方法包括梯 度下降法、牛顿法和拟牛 顿法等。
数值分析在人工智能领域的应用
总结词
数值分析在人工智能领域的应用关键，涉及深度学习、神经 网络等领域。
详细描述
数值分析为人工智能提供了理论基础和算法支持，特别是在 深度学习和神经网络方面。通过数值分析的方法，可以优化 神经网络的参数和结构，提高人工智能的性能和准确性。
数值分析在金融领域的应用
总结词
常见的迭代法有雅可比迭代法 、高斯-赛德尔迭代法等。
牛顿法
牛顿法是一种基于泰勒级数 的迭代方法，用于求解非线 性方程的根。
数值分析考试卷及详细答案解答汇总
姓名 __________ 班级 ___________ 学号 _____________一、选择题i.F （2,5,-3,4）表示多少个机器数（C ）.A 64B 129C 257D 256 2. 以下误差公式不正确的是（D ）A ・ £（迎 *一七 *）« 5（Xj*）+£（£ *） c ，£（“*•£ *）«|^2 *k （-'l*） + |时住2 *）3. 设° =（、任_1）6,从算法设计原则上定性判断如下在数学上等价的表达式，哪一个在数值计算上将给出°较好的近似值？ （D ）A ———B 99-70V2C （3-2V2）3D —— （V2 +1）6 （3 + 204. 一个30阶线性方程组，若用Crammer 法则来求解，则有多少次乘法？（A ） A31X29X30! B 30X30X30! C31X30X31! D 31X29X29!5. 用一把有亳米的刻度的米尺来测量桌子的长度，读出的长度1235mm,桌子的精确长度 记为（D ） A 1235mm B 1235-0.5mm C 1235+0.5nun D 1235±0.5mm二、填空1. 构造数值算法的基本思想是 近似替代、离散化、递推化 。
2. 十进制123.3转换成二进制为1111011.0而1。
3. 二进制110010.1001转换成十进制为 50.5625 。
4. 二进制o.ioi 转换成十进制为-o75.已知近似数X *有两位有效数字，则其相对误差限 5%。
6.1112=0.69314718...,精确到 10一’的近似值是 0.693。
* *7. x = ;r = 3.1415926・・・，则“ =3.1416 , =3.141的有效数位分别为5 和 3 __________ o8. 设卅=2.001,严=-0.8030是由精确值x 和y 经四舍五入得到的近似值，则兀* +y *的误差限____________________ o9.设x = 2.3149541•…，取5位有效数字，则所得的近似值卅二2.3150 。
数值分析题库答案(含详细解题步骤)
第 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 迭代格式，并判断是否收敛。
《数值分析》所有参考答案
等价三角方程组
， ，
11.设计算机具有4位字长。分别用Gauss消去法和列主元Gauss消去法解下列方程组，并比较所得的结果。
解：Gauss消去法
回代
列主元Gauss消去
15.用列主元三角分解法求解方程组。其中
A= ,
解：
等价三角方程组
回代得
， ， ，
16.已知 ，求 ， ， 。
解：
, ,
17.设 。证明
，(II)
，
当 时
当 时
迭代格式(II)对任意 均收敛
3) ，
构造迭代格式 (III)
，
当 时
当 时
迭代格式(III)对任意 均收敛
4)
取格式(III)
， ， ，
4.用简单迭代格式求方程 的所有实根，精确至有3位有效数。
解：
当 时， ，
1 2
当 时
，
，
， ，
1)
迭代格式 ，
,
当 时， ，
任取 迭代格式收敛于
是中的一种向量范数。
解：
当 时存在 使得
，
,
所给 为 上的一个范数
18.设 。证明
(1) ;
(2) ;
(3) 。
解：(1)
(2)
(3)
19.设
A=
求 ， ， 及 ， 。
解： ,
Newton迭代格式
,
20.设 为 上任意两种矩阵(算子)范数，证明存在常数
， 使得
对一切 均成立。
解：由向量范数的等价性知道存在正常数 使得
,
=0.187622
[23.015625 , 23.015625+0.187622]
数值分析习题与解析
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x -==≈⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯ 所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x --==≈-⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.14＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯ 所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈- 相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯ 有效数字： 因为π＝3.14159265...=0.314159265 (10)22 3.1428571430.3142857143107==⨯，m=1。
《数值分析》复习题(14)解析
《数值分析》一、填空题1. 已知近似数 1.28y *=-，则其绝对误差限为 0.005 ，相对误差限是 0.0039 。
2. 设 1.28057x *=是准确值 1.28367x =的近似值，则x 有 3 位有效数字，相对误差限是 0.0024 。
3. 求方程sin x x =的Newton 迭代公式是 。
5. 在数值计算中，计算1cos 2.6-应变成 来计算。
6. 由秦九韶算法，543262975y x x x x x =-++--应改写为 。
8. 已知3132A -⎡⎤=⎢⎥⎣⎦，则谱半径()A ρ= 32. 测量一支铅笔长是16cm ， 那么测量的绝对误差限是 0.5cm ，测量的相对误差限是3.1% 。
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 迭代格式： 。
《数值分析》练习题及答案解析
《数值分析》练习题及答案解析一、单选题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 。
(完整版)数值分析第四版习题和答案解析
第四版数值分析习题第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====⨯4. 利用公式(3.3)求下列各近似值的误差限:********12412324(),(),()/,i x x x ii x x x iii x x ++其中****1234,,,x x x x 均为第3题所给的数.5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少?6. 设028,Y =按递推公式1n n Y Y -=…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2? 10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小.11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好?3--13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式ln(ln(x x =-计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xxx ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且 101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj j j x l x xk n =≡=∑ii)()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()x f x e =的等距节点函数表,若用二次插值求x e 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n y =,求4n y ∆及4n y δ. 10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)k f x k m ∆≤≤是m k -次多项式,并且()0(m l f x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i) 若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦及0182,2,,2f ⎡⎤⎣⎦.17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限. 19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()nx ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差.试求三次样条插值并满足条件 i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii) (0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n ==,式中i x 为插值节点,且01n a x x x b =<<<=,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式. 8. 如何选取r,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一?9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式. 13. 设()xf x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差.15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005. 16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数. 17. 求a 、b 使[]22sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义()(,)()();()(,)()()()();bbaaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x=在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.sin (1)arccos ()nn x u x +=是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.2y a bx =+.用最小二乘拟合求.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度:(1)101()()(0)()h h f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)1,4n =⎰; (4),6n =.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分10x e dx-⎰并计算误差. 5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2ba f f x dxb a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰. 6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8.1xedx-,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是S a =θ,这里a 是椭圆的半长轴,c 是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长. 10. 证明等式3524sin3!5!n n nnππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计()f x第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
数值分析整理版试题及答案
例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[。
数值分析试题库与答案解析
模 拟 试 卷（一）一、填空题（每小题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.()2534F,,,-表示多少个机器数(C ).A 64B 129C 257D 2562. 以下误差公式不正确的是( D)A ．()()()1212x *x *x *x *εεε-≈+B ．()()()1212x *x *x *x *εεε+≈+C ．()()()122112x *x *x *x *x x *εεε?≈+ D ．()()()1212x */x *x *x *εεε≈-3. 设)61a =, 从算法设计原则上定性判断如下在数学上等价的表达式，哪一个在数值计算上将给出a 较好的近似值？(D )A6)12(1+ B 27099- C 3)223(- D3)223(1+4. 一个30阶线性方程组, 若用Crammer 法则来求解, 则有多少次乘法? ( A )A 31×29×30!B 30×30×30!C 31×30×31!D 31×29×29!5. 用一把有毫米的刻度的米尺来测量桌子的长度, 读出的长度1235mm, 桌子的精确长度记为（ D ）A 1235mmB 1235-0.5mmC 1235+0.5mmD 1235±0.5mm二、填空1.构造数值算法的基本思想是近似替代、离散化、递推化。
2.十进制123.3转换成二进制为1111011.01001。
3.二进制110010.1001转换成十进制为 50.5625 。
4. 二进制0101.转换成十进制为57。
5.已知近似数x*有两位有效数字，则其相对误差限5% 。
6. ln2=0.69314718…，精确到310-的近似值是 0.693 。
7.31415926x .π==，则131416*x .=，23141*x .=的有效数位分别为5 和3 。
8.设200108030x*.,y*.==-是由精确值x y 和经四舍五入得到的近似值，则x*y*+的误差限0.55×10-3 。
数值分析试题(卷)与答案解析
数值分析试题一、 填空题（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 。
x 0 0.5 1 1.5 2 2.5 y =f (x )-2-1.75-10.2524.2511. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。
(汇总)东北大学-数值分析--考试题解析.ppt
构造函数(t)=(t)-H3(t)-C(x)t(t-1)2(t-2) 于是,存在x,使(4)(x)=0,即(4)(x)-4!C(x)=0
R(x) f (4) ( x ) x(x 1)2 (x 2)
4!
五、(12分)试确定参数A,B,C及,使数值积分公式
2
2
f
(x)dx
Af
( )
Bf
(0)
Cf
( )
有尽可能高的代数精度,并问代数精度是多少?它是否是
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=1精0品/文9档,B=16/9,=(12/5)1/2 7
令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),
于是
H3(x)==-x(3x--21.)5x2(2x+-22.)5-x3+x2(精x品-2文)档+2.5x(x-1)2
–0.5x(x-1)(x-2) 6
由于,R(0)=R(1)=R(2)=R(1)=0, 故可设 R(x)=C(x)x(x-1)2(x-2)
(3)因为0<</2,所以() cos / 2 1 sin 0
故,此迭代法线性收敛(收敛阶为1).
三、(14分)设线性方程组
4x1 x2 2x3 1 x1 5x2 x3 2 2x1 x2 6x3 3
(1)写出Jacobi法和SOR法的迭代格式(分量形式);
(2)讨论这两种迭代法的收敛性.
数值分析试题库与答案解析
A1 f (1)的求积公式，并求出
3．用 Newton 法求方程 x ln x 2 在区间 ( 2, ) 内的根 , 要求 xk xk 1 xk
4．用最小二乘法求形如 y a bx 2 的经验公式拟合以下数据：
10 8 .
xi
19
25
30
38
yi
19.0
32.3
49.0
73.3
5．用矩阵的直接三角分解法解方程组
， l 32 2
3
3
解方程组
1
y1
11
y2
5
y3
3 21
10
4
16 得 y1 10, y2 6, y3
，
3
30
再解方程组
5 1 1 3 x1
1 2 x2 1 x3
d1 1
10
d2 1
6 得 x1 1, x2
d3 1 4
3
1, x3 2 .
1 4 解 令 Y ，则 Y a bx 容易得出正规方程组
y
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
由于 ( x) [ x f ( x)] 1 f ( x) ，所以 | ( x) | |1 f ( x) | 1
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
三、(14分)设线性方程组 4x1 x2 2x3 1 x1 5x2 x3 2 2x1 x2 6x3 3
(1)写出Jacobi法和SOR法的迭代格式(分量形式);
七、(6分)设n阶矩阵A=(aij)nn,试证实数
A
n
max
1i, jn
aij
为矩阵A的一种范数.
证明 对任意n阶方阵A,B和常数,有
A
n
max
1i, jn
aij
0,且仅当A 0时 A
0。
A
n
max
1i, jn
aij
n max 1i, jn
aij
A
[
2 fn x2
4 9
h2
2 fn xy
8 9
h2
fn
2 fn y 2
4 9
h2
f
2 n
O(h3)
于是有
yn1
yn
hf n
1 2
h2 ( fn x
fn y
fn)
1 6
h3[
2 fn x2
2 2 fn xy
fn
2 fn y 2
f
2 n
]
AB
n
max
1i, jn
aij
bij
n max[ aij | | bij ]
1i, jn
A
B
n
n
AB
n
max
1i, jn
k
1
aik
bkj
n max 1i, jn
aij
(max 1 i n
aik
k 1
)
n
max
1i, jn
aij
n
max
所以此单步方法为二阶方法.
(2)此单步方法用于方程y=-10y,则有
yn1
yn
h 4
[10
yn
30( yn
20 3
hyn )]
[110h
50h2 ]yn
当h=0.25时,有
110h 50h2 1 2.5 3.125 1.625 1
所以,所得数值解是不稳定的.
因此方程(x)=0有唯一正根,且在区间(1,2)内.
(2)构造迭代格式: xk1 1 sin xk
k 0,1,2,...
由于|(x)|=| cos x / 2 1 sin x|<1,故此迭代法收敛.
取初值x0=1.5, 计算得x1=1.41333, x2=1.40983,由于 |x2-x1|=0.0035<10-2 , 故可取根的近似值x2=1.40983.
考试题解析
一、填空题（每空3分,共30分)
1.设矩阵A= 1
2
2 3
,则(A)=____7___,Cond(A)1=___2_5/_7__.
解 由于
1 A E
2 2 4 7 0
2 3
得特征值: 1 2 3i, 2 2 3i
1 5
x(k 1
)
1 5
x(k 3
)
2 5
x(k 3
1)Fra bibliotek1 3x(k 1
)
1 6
x(k) 2
1 2
x1(
k
1)
x (k ) 1
(x1(k )
1 4
x (k ) 2
1 2
x(k 3
)
1) 4
x
(k 2
1)
x (k ) 2
(1 5
x(k 1
容易验证公式对(x)=x5仍精确成立，故其代数精度为5， 是Gauss公式。 六、(12分)设初值问题
y f (x, y)
y(a)
(1)试证单步法
a xb
K1 f (xn , yn ),
K2
f
(xn
2 3
h,
yn
2 3
hK1
)
y
n
1
yn
h 4
(
K
1
解 只要取(x)=x3-a ,或(x)=1-x3/a. 5.设(x)=x3+x2-3,则差商[3,32,33,34]=__1_____.
6.设l0(x),l1(x),l2(x),l3(x)是以x0,x1,x2,x3为互异节点
的三次插值基函数,则
3
l j (x)(x j
2)3=___(_x_-_2_)_3____.
O(h
4
)
而
y(xn1 )
y(xn
)
hy ( x n
)
1 2
h2
y ( x n
)
1 6
h3
y ( x n
)
O(h 4
)
yn hfn
1 h 2 [f n 2 x
f n y
f
n
]
1 6
h3
y ( x n
)
O(h
4
)
所以有
y(xn1 ) yn1 O(h3 )
j0
7.设S(x)=
x3 x2
2x
3
bx 2
cx 1
0 x 1 是以0,1,2为节 1 x 2
点的三次样条函数,则b=__-__2____c=___3______.
解 由2=b+c+1,5=6+2b+c,8=12+2b,可得
二、(13分)设函数(x)=x2-sinx-1
又A-1= 1 3
72
2 1
,所以‖A‖1=5,‖A-1‖1=5/7.
2.设矩阵A=
1 a
a 1
a 0
,当a取______值时,A可以唯一分解
a 0 1
为GGT,其中G为下三角矩阵.
1aa
解
令
1 a
a 1 a 2 0, a 1
a
1 0
0 1 2a 2 0, 得： 1 a 1
于是 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
由于,R(0)=R(1)=R(2)=R(1)=0, 故可设 R(x)=C(x)x(x-1)2(x-2)
构造函数(t)=(t)-H3(t)-C(x)t(t-1)2(t-2) 于是,存在x,使(4)(x)=0,即(4)(x)-4!C(x)=0
(2)讨论这两种迭代法的收敛性.
(3)取初值x(0)=(0,0,0)T,若用Jacobi迭代法计算时,
预估误差x*-x(10) (取三位有效数字).
解 (1)Jacobi法和SOR法的迭代格式分别为
x1(
k
1)
1 4
x(k) 2
1 2
x(k 3
)
1 4
x
(k 2
1)
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,得 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),
3K 2 )
n 0,1,2,...
y0
是二阶方法.
(2)以此法求解y=-10y, y(0)=1时,取步长h=0.25,所得
数值解yn是否稳定?为什么? 解 (1)由于
K2
f
(xn
2 3
h,
yn
2 3
hK1
)
fn
fn x
2h 3
fn y
2 3
hf
n
1 2
1
2
2
3.向量x=(x1,x2,x3)T,试问|x1|+|2x2|+|x3|是不是一种向 量范数__是____,而|x1|+|2x2+x3|是不是一种向量范数_不_是___.
4.求 3 a 的Newton迭代格式为_x_k_1 __xk__x_3k3_xk_2 a_或_x_k_1__32_x_k __a3_x_k2_.
得x(1)=(1/4,-2/5,1/2)T,于是x(1)-x(0)=1/2,所以有
x* x (10)
Bk
x (1) x (0)
0.7510 0.5 0.113