2014数值分析(A)

第 1 页 共 6 页 （A 卷）

考试方式：闭卷 太原理工大学 数值分析 试卷（A ） 适用专业：2012级数学专业 考试日期：2015.1.5 时间：120 分钟 共 6页

一．判断题：每小题2分，共20分.

1．两个相近数相减必然会使有效数字损失. （ ） 2．() (0,1,,)i l x i n = 是关于节点i x 的拉格朗日插值基函数,则对于任何次数不大于n 的多项式 ()P x 都有0()()()n i i i l x P x P x ==∑. （ ） 3．任何()[,]f x C a b ∈都能找到n 次多项式()n n P x H ∈，使()()n f x P x ε-≤（ε为给定的误差限. （ ） 4．1n +个点的插值型求积公式的代数精度至少n ，最多可达21n +次. （ ） 5．数值求积公式总是稳定的. （ ） 6. 如果矩阵对称，则1A A ∞=. （ ） 7. A 正定对称，则雅克比迭代与高斯—赛德尔迭代法都收敛. （ ） 8. 实矩阵的特征值一定是实数. （ ） 9. 改进的欧拉法是二级二阶的龙格—库塔方法. （ ） 10．幂法是求矩阵按模最大特征值及相应特征向量的方法. （ ）

第 2 页 共 6 页 （A 卷）

二、填空题：每小题2分，共分10分.

11. 设123456734A ?? ?= ? ???

，则1A =______________,A ∞=_______________. 12. 已知

则012[,,]f x x x =_____________.

13. 1112

A ??=

?-??, 则()A ρ=_________,1()cond A =_____________. 14. * 5.00000x =是经过四舍五入得到的近似数，*x 有_______________位有效数字,误差限为 ____________________.

15. 设i x 为互异的节点(0,1,,)i n = ，则

0()()n k i i i x l x ==∑_______________,(0,1,2,,)k n = 0

()()n k i i

i x x l x =-=∑__________________________.(1,2,,)k n =

三．本题共1小题， 共10分.

16. 计算球体积要使相对误差限为0.01，问度量半径R 所允许的相对误差限是多少？

四．本题共3小题，共30分.

17.（10分）求一个次数小于等于3的多项式

3()

P x，使其满足下列插值条件：

并估计误差.

第 3 页共6 页（A卷）

第 4 页 共 6 页 （A 卷） 18. （10分）设2()21f x x x =++，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{1,}span x Φ=

的最佳平方逼近多项式.

19. （10分）推导求积公式：

3()()()()()224

b a a b f f x dx b a f b a η''+=-+-?

第 5 页 共 6 页 （A 卷）

五．本题共3小题， 共30分.

20. （10分） 对线性方程组12211121x x ??????= ? ? ???????，若用迭代法(1)()()()k k k x x Ax b α+=+-求解，问α在什么范围内取值可使迭代法收敛. 21. （10分）利用格什戈林圆盘定理证明n 矩阵211211141112n n n A n n n n n n n ?? ?

? ? ?= ? ? ? ??? 能与对角矩阵相似，且A 的特征值都是实数.

22. （10分）用改进的欧拉法解初值问题

2

(0)0

y x y

y

'

?=-

?

=

?

，取步长0.1,

h=计算到（保留到小数后4

位）。

第 6 页共6 页（A卷）

数值分析试题(08研)

数值分析试题 一. 填空题： 1. 设A=?? ????4311，则 ||A||1 = ，||A||∞ = _______，()A ρ=_________； 2. 已知函数()y f x =的观测数据为(0,1),(1,2),(2,3，则二次Lagrange 插值多项式22()L x a bx cx =++中a = ， b =_____ ， c =_____； 3. 为使求积公式012()()(0)()h h f x dx A f h A f A f h -≈-++?的代数精度尽量高，则0A =_____，1A =______，2A =______，其具有代数精度为_____次； 4. 设给出(1)2,(0)1,(1)0,(0)2f f f f '-====-，可求得其三次插值多项式 233()H x a bx cx d x =+++中a =____，b =_____ ，c =______ ，d =_____； 5．对3()31f x x x =++，差商[0,1,2,3]f = ；[0,1,2,3,4]f = 。 二．已知函数()y f x =的观测数据为： 1．构造差商表，并写出Newton 插值多项式（按降幂排列）； 2．用最小二乘法求形如 2y a bx cx =++的经验公式使与题目数据拟合； 3．用复化梯形公式计算4 1()f x dx ?的近似值。 三．分别用下列方法求方程3310x x +-=在[0,1]内的根使误差小于110-： 1. Newton 法（取00.4x =）； 2. 试证明用简单迭代格式3/)1(31k k x x -=+求其在[0.2,0.4]内的根是收敛的。 四. 用下列各种方法求解方程组Ax b =，即 ??????????-122111221????????321x x x =???? ??????-001 1．Gauss 消元法； 2．Doolittle 分解法； 3. 写出求Ax b =的解的Jacobi 迭代格式，并取(0)(0,0,0)T x =求(3)x ； 4. 判定矩阵A 对Jacobi 迭代的收敛性，并证明你的结论。 五．1．用2段Simpson 公式（5节点）计算?511dx x 的近似值（计算中取五位有效数字）； 2．若使误差不超过610-，用复化梯形公式计算上述积分至少应取多少个节点？

数值分析课后题答案

数值分析 第二章 2．当1,1,2x =-时，()0,3,4f x =-,求()f x 的二次插值多项式。 解： 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23 537623 l x l x x x x x x x =-+=---+ -+= +- 6．设,0,1,,j x j n =L 为互异节点，求证： （1） 0()n k k j j j x l x x =≡∑ (0,1,,);k n =L （2）0 ()()0n k j j j x x l x =-≡∑ (0,1,,);k n =L 证明 （1） 令()k f x x = 若插值节点为,0,1,,j x j n =L ，则函数()f x 的n 次插值多项式为0 ()()n k n j j j L x x l x == ∑。 插值余项为(1)1() ()()()()(1)! n n n n f R x f x L x x n ξω++=-= + 又,k n ≤Q

(1)()0 ()0 n n f R x ξ+∴=∴= 0()n k k j j j x l x x =∴=∑ (0,1,,);k n =L 0 000 (2)()() (())()()(()) n k j j j n n j i k i k j j j i n n i k i i k j j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑ 0i n ≤≤Q 又 由上题结论可知 ()n k i j j j x l x x ==∑ ()()0 n i k i i k i k C x x x x -=∴=-=-=∑原式 ∴得证。 7设[]2 (),f x C a b ∈且()()0,f a f b ==求证： 21 max ()()max ().8 a x b a x b f x b a f x ≤≤≤≤''≤- 解：令01,x a x b ==，以此为插值节点，则线性插值多项式为 10 101010 ()() ()x x x x L x f x f x x x x x --=+-- =() () x b x a f a f b a b x a --=+-- 1()()0()0 f a f b L x ==∴=Q 又 插值余项为1011 ()()()()()()2 R x f x L x f x x x x x ''=-= -- 011 ()()()()2 f x f x x x x x ''∴= --

数值分析第三版课本习题及答案

第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位 有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1 234,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y . (五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字 . 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 2 12S gt = 假定g 是准确的,而对t 的测量有±秒的误差,证明当t 增加时S 的绝对误差增加,而相 对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…), 若0 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一 等价公式 ln(ln(x x =- 计算,求对数时误差有多大? 14. 试用消元法解方程组 { 101012121010;2. x x x x +=+=假定只用三位数计算,问结果是否可靠?

研究生数值分析试卷

2005~2006学年第一学期硕士研究生期末考试试题（A 卷） 科目名称：数值分析 学生所在院： 学号： 姓名： 注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。 一、（15分）设求方程 0cos 2312=+-x x 根的迭代法 k k x x cos 3 2 41+=+ (1) 证明对R x ∈?0,均有*lim x x k k =∞ →,其中*x 为方程的根. (2) 此迭代法收敛阶是多少? 证明你的结论. 二、（12分）讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的收敛性。 ??? ??=++-=++=-+. 022,1, 122321 321321x x x x x x x x x 三、（8分）若矩阵??? ? ? ??=a a a a A 000002，说明对任意实数0≠a ，方程组b AX =都是非病态的。（范数用∞?） 四、（ 求)(x f 的Hermite 插值多项式)(3x H ，并给出截断误差)()()(3x H x f x R -=。 五、（10分）在某个低温过程中，函数 y 依赖于温度x （℃）的试验数据

为 已知经验公式的形式为 2bx ax y += ，试用最小二乘法求出 a ，b 。 六、（12分）确定常数 a ，b 的值，使积分 [ ] dx x b ax b a I 2 1 1 2 ),(?--+= 取得最小值。 七、（14分）已知Legendre(勒让德)正交多项式)(x L n 有递推关系式： ?? ? ? ???=+-++===-+),2,1()(1)(112)()(, 1)(1110 n x L n n x xL n n x L x x L x L n n n 试确定两点的高斯—勒让德（G —L ）求积公式 ? -+≈1 1 2211)()()(x f A x f A dx x f 的求积系数和节点，并用此公式近似计算积分 ?=2 11 dx e I x 八、（14分）对于下面求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy 的单步法： ??? ? ??? ++==++=+) ,() ,()2 121(1 21211 hk y h x f k y x f k k k h y y n n n n n n

2015年研数值分析A卷

武 汉 大 学 2015-2016第一学期硕士研究生期末考试试题（A 卷） 科目： 数值分析 学生所在院： 学号： 姓名： 一、（12分）设方程230x x e -=，为求其最大正根与最小正根的近似值，试分别确定两个含根区间[,]a b 和两个迭代函数()g x ，使当0[,]x a b ?时，迭代格式1()n n x g x +=分别收敛于最大正根与最小正根。 二、（12分）用杜利特尔（Doolittle ）分解算法求解方程 b Ax =，其中 211625608A ????=?????? 226768b ????=?????? 三、（14分）设方程组 123121113a a x a a x a a x 轾轾轾犏犏犏犏犏犏=-犏犏犏犏犏犏臌臌臌 其中a 为常数。 （1）分别写出Jacobi 迭代格式及 Gauss-Seidel 迭代格式； （2）导出Gauss-Seidel 迭代格式收敛的充分必要条件。 四、（12分）已知 )(x f y = 的数据如下： 求)(x f 的Hermite 插值多项式)(3x H 及其余项。 五、（12分）确定常数 a ，b 的值，使积分 2 1 320(,)I a b x ax bx dx 轾=--犏臌ò 取得最小值。

六、（12 求形如 y bx x =+ 的拟合曲线。 七、（14分）（1）对初值问题 00(,)[,]()dy f t y t a b dt y t y ì??= ??í??=?? 验证改进欧拉方法（也称预估-校正法）与微分方程是相容的； （2） 用改进欧拉方法求下面方程的数值解（取步长5.0=h ）： (0)1 dy dt y ?=???=? [0,1]t ∈ （取5位有效数字计算） 八、（12分）设求积公式 ∑?=≈n k k k b a x f A dx x f 1)()(为高斯型求积公式， 并记 )())(()(21n n x x x x x x x ---= ω （1）问给定的求积公式的代数精度是多少次？ （2）证明: 对任意次数小于等于1-n 的多项式)(x q ，必有?=b a n dx x x q 0)()(ω； （3）证明：n k A k ,,2,1,0 =>

第五章习题解答_数值分析

第五章习题解答 1、给出数据点:0134 19156 i i x y =?? =? (1)用012,,x x x 构造二次Lagrange 插值多项式2()L x ,并计算15.x =的近似值215(.)L 。 (2)用123,,x x x 构造二次Newton 插值多项式2()N x ,并计算15.x =的近似值215(.)N 。 (3)用事后误差估计方法估计215(.)L 、215(.)N 的误差。 解: (1)利用012013,,x x x ===,0121915,,y y y ===作Lagrange 插值函数 2 20 2 1303011915 01031013303152933 ()()()()()() ()()()()()()()() i i i x x x x x x L x l x y x x =------== ?+?+?-------++= ∑ 代入可得2151175(.).L =。 (2)利用 134,,x x x ===,9156,,y y y ===构造如下差商表： 229314134196()()()()()N x x x x x x =+-+---=-+- 代入可得215135(.).N =。 (3)用事后误差估计的方法可得误差为 ()()()02222 03-x 150 x x x -=117513506563-04.()()()(..).x f L R L x N x x x --≈= -≈- ()()()3222203-154 x x -=1175135-1.0938-04 .()()()(..)x x f N R x L x N x x x --≈=-≈- 2、设Lagrange 插值基函数是 0012()(,,,,)n j i j i j j i x x l x i n x x =≠-==-∏ 试证明:①对x ?,有 1()n i i l x ==∑ ②00110001211()()(,,,)()()n k i i i n n k l x k n x x x k n =?=?==??-=+? ∑ 其中01,,,n x x x 为互异的插值节点。 证明： ①由Lagrange 插值多项式的误差表达式10 1()()()()()!n n i i f R x x x n ξ+==-+∏知，对于函数1()f x =进行

07(研)数值分析

数值分析试题 2007.12 一、简答下列各题：（每题4分，共20分） 1．为了提高计算精度，求方程x 2－72x+1=0的根，应采用何种公式，为什么？ 2．设??? ? ??=2112A ，求)(A ρ和2)(A Cond 。 3．设??? ? ? ??=131122321A ，求A 的LU 分解式。 4．问23221)2(x x x x ++=是不是3R 上的向量范数,为什么？ 5．求数值积分公式?-≈b a a b a f dx x f ))(()(的截断误差R[?]。 二、解答下列各题：（每题8分，共56分） 1．已知线性方程组??? ??=-+=-+=-+3 53231 4321 321321x x x x x x x x x ，问能用哪些方法求解？为什么？ 2.解线性方程组b Ax =的Gauss-Seidel 迭代法是否收敛？为什么？其中： ???? ? ??--=211111112A 3.设]2,0[)(4C x f y ∈=，且0)0(,0)2(,2)1(,1)0(='===f f f f ，试求)(x f 的三次插值多项式)(3x H ，并写出余项)()()(33x H x f x R -=。 4．给定离散数据 试求形如3bx a y +=的拟合曲线。 5．求区间[0，1]上权函数为x x =)(ρ的正交多项式)(0x p ，)(1x p 和)(2x p 。 6．证明求积公式： ? +++-≈3 1 ) 5 3 2(5)2(8)532(5[91)(f f f dx x f

是Gauss 型求积公式。 7. 利用2=n 的复化Simpson 公式计算计算定积分 ，并估计误差][f R 。 三、（12分）已知方程0cos 2=-x x ， 1.证明此方程有唯一正根α； 2.建立一个收敛的迭代格式，使对任意初值]1,0[0∈x 都收敛，说明收敛理由和收敛阶。 3.若取初值00=x ，用此迭代法求精度为510-=ε的近似根，需要迭代多少步？ 四、（12分）已知求解常微分方程初值问题： ?? ?∈=='] ,[,)(),(b a x a y y x f y α 的差分公式： ?? ??????? =++==++=+α 0121211) 32 ,32() ,()3(4y hk y h x f k y x f k k k h y y n n n n n n 1．证明：此差分公式是二阶方法； 2．用此差分公式求解初值问题1)0(,10=-='y y y 时，取步长h=0.25,所得数值解是否稳定,为什么？ ?1 0sin xdx

数值分析第四版习题及答案

第四版 数值分析习题 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式求下列各近似值的误差限: ********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 .(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字. 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±秒的误差,证明当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 =- 计算,求对数时误差有多大?

研究生《数值分析》教学大纲

研究生《数值分析》教学大纲 课程名称：数值分析 课程编号：S061005 课程学时：64 学时 课程学分： 4 适用专业：工科硕士生 课程性质：学位课 先修课程：高等数学，线性代数，计算方法，Matlab语言及程序设计 一、课程目的与要求 “数值分析”课是理工科各专业硕士研究生的学位课程。主要介绍用计算机解决数学问题的数值计算方法及其理论。内容新颖，起点较高，并加强了数值试验和程序设计环节。通过本课程的学习，使学生熟练掌握各种常用的数值算法的构造原理和过程分析，提高算法设计和理论分析能力，并且能够根据数学模型，提出相应的数值计算方法编制程序在计算机上算出结果。力求使学生掌握应用数值计算方法解决实际问题的常用技巧。 二、教学内容、重点和难点及学时安排： 第一章? 数值计算与误差分析( 4学时) 介绍数值分析的研究对象与特点，算法分析与误差分析的主要内容。 第一节数值问题与数值方法 第二节数值计算的误差分析 第三节数学软件工具----MATLAB 语言简介 重点：误差分析 第二章? 矩阵分析基础( 10学时) 建立线性空间、赋范线性空间、内积空间的概念，为学习以后各章打好基础。矩阵分解是解决数值代数问题的常用方法，掌握矩阵的三角分解、正交分解、奇异值分解，并能够编写算法程序。 第一节? 矩阵代数基础

第二节? 线性空间 第三节? 赋范线性空间 第四节? 内积空间和内积空间中的正交系 第五节矩阵的三角分解 第六节矩阵的正交分解 第七节矩阵的奇异值分解 难点：内积空间中的正交系。矩阵的正交分解。 重点：范数，施密特(Schmidt) 正交化过程，正交多项式，矩阵的三角分解, 矩阵的正交分解。 第三章? 线性代数方程组的数值方法( 12学时) 了解研究求解线性代数方程组的数值方法分类及直接法的应用范围。高斯消元法是解线性代数方程组的最常用的直接法，也是其它类型直接法的基础。在此方法基础上加以改进，可得选主元的高斯消元法、按比例增减的高斯消元法，其数值稳定性更高。掌握用列主元高斯消元法解线性方程组及计算矩阵的行列式及逆，并且能编写算法程序。掌握矩阵的直接三角分解法：列主元LU 分解，Cholesky分解。了解三对角方程组的追赶法的分解形式及数值稳定性的充分条件。掌握矩阵条件数的定义，并能利用条件数判别方程组是否病态以及对方程组的直接方法的误差进行估计。 迭代解法是求解大型稀疏方程组的常用解法。熟练掌握雅可比迭代法、高斯- 塞德尔迭代法及SOR 方法的计算分量形式、矩阵形式，并能在计算机上编出三种方法的程序用于解决实际问题。了解极小化方法：最速下降法、共轭斜量法。迭代法的收敛性分析是研究解线性代数方程组的迭代法时必须考虑的问题。对于上述常用的迭代法，须掌握其收敛的条件。而对一般的迭代法，掌握其收敛性分析的基本方法和主要结果有助于进一步探究新的迭代法。 第一节求解线性代数方程组的基本定理 第二节高斯消元法及其计算机实现 第三节矩阵分解法求解线性代数方程组 第三节? 误差分析和解的精度改进 第四节? 大型稀疏方程组的迭代法 第五节? 极小化方法 难点：列主元高斯消元法，直接矩阵三角分解。迭代法的收敛性，雅可比迭代法，高斯-塞德尔迭代法，SOR 迭代法。

数值分析(研)试题答案

沈阳航空航天大学研究生试卷（A） 2011-2012 学年第一学期课程名称：数值分析出题人: 王吉波审核人: 一、填空题（本题40 分每空 4 分） 1．设l j (x) ( j 0 ,1, ,n) 为节点x0 , 1 , , x 的n 次基函数，则 l j ( x i ) x n 1, 0, i i j j 。 2．已知函数(x) x 1 f 2 x ，则三阶差商 f [1, 2, 3, 4] = 0 。 3．当n=3 时，牛顿- 柯特斯系数 1 (3) 3 (3) (3) C0 , C C ，则 1 2 8 8 (3) C 3 1 8 。 ( ) Bx( k) f k k 1 收敛的 4．用迭代法解线性方程组Ax=b时，迭代格式, 0,1,2 , x 充分必要条件是(B) 1或B 的谱半径小于 1 。 5．设矩阵 1 2 A ，则A 的条件数 Cond (A)2 = 3 。 2 1 6.正方形的边长约为100cm，则正方形的边长误差限不超过0.005 cm 才能使 其面积误差不超过1 2 cm 。 1 1 7.要使求积公式(0) ( ) 8. f (x)dx f A1 f x1 具有 2 次代数精确度，则 4 x 2/3 ，A1 3/4 。 1 9 18 9 - 27

18 45 0 - 45 其 中， A 8. 用杜利特尔（Doolittle ）分解法分解 A LU ， 9 0 126 9 27 -45 9 135 则 1 1 1 2 3 1 2 L ， 1 - 2 1 0 3 U 9 18 9 9 -18 81 - 27 9 54 9

研究生数值分析习题

1. 五个节点的Newton-Cotes 求积公式的代数精度为______，五个节点的求积公式最高代数精度为___________。（即Gauss 型求积公式） 2. 已知数值求积公式为3 11 ()[(1)4(2)(3)]3 f x dx f f f ≈++? ， 则其代数精度为______。 3. 数值积分公式1 '12 ()[(1)8(0)(1)]9 f x dx f f f -≈-++?的代数 精度为_________。 4. 要使求积公式1 110 1 ()(0)()4 f x dx f A f x ≈ +?具有2次代数精度，则1x =___,1A =___。 5. 在Newton-Cotes 求积公式：() ()()()n b n i i a i f x dx b a C f x =≈-∑? 中，当系数()n i C 是负值时，公式的稳定性不能保证，所以实际应用中，当___________时的Newton-Cotes 求积公式不能使用。 ()8()7()10()6A n B n C n D n ≥≥≥≥ 6. 若用复化梯形公式计算1 0x e dx ?，要求误差不超过6 10-，利 用余项公式估计，至少用______个求积节点。 7. 对于Gauss 型求积公式3 1 ()()()b k k a k f x x dx A f x ρ=≈∑?，其中 ()x ρ为权函数，下列说法错误的是_________。

（A ）该求积公式一定是稳定的； （B ）3 1()k k k A f x b a ==-∑； （C ）该求积公式的代数精度为5； （D ）2 (35)()()0b a x x x x dx ωρ-=? ，其中3 1 ()()k k x x x ω==∏-。 8. 0{()}k k x ?∞ =是区间[0,1]上权函数 ()x x ρ=的最高系数为1的正交多项式族，其中0()1x ?=，则1 40()_______x x dx ?=?。 9. 构造代数精度最高的如下形式的求积公式，并求出其代数精度： 1 010 1 ()()(1)2 xf x dx A f A f ≈+? 10. 数值积分公式形如 1 ()()(0)(1)(0)(1)xf x dx S x Af Bf Cf Df ''≈=+++? （1）试确定参数A 、B 、C 、D ，使公式的代数精度尽量高； （2）设4 ()[0,1]f x C ∈，推导余项公式1 0()()()R x xf x dx S x =-?， 并估计误差。 11. 用8n =的复化梯形公式和复化Simpson 公式计算 1 x e d x -? 时， （1）试用余项估计其误差； （2）计算积分的近似值。

数值分析第四版习题和答案解析

第四版 数值分析习题 第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2％,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝ 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而 相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大这个计算过程 稳定吗 12.计算,取,利用下列等式计算,哪一个得到的结果最好 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大若改用另一等价公式 计算,求对数时误差有多大 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足 第二章插值法 1.根据定义的范德蒙行列式,令 证明是n次多项式,它的根是,且 . 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.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函 数表的步长应取多少 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是 并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误 差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误 差做比较. 2.求证: (a)当时,. (b)当时,. 3.在次数不超过6的多项式中,求在的最佳一致逼近多项式.

《数值分析》第五章答案

习题5 1．导出如下3个求积公式，并给出截断误差的表达式。 (1) 左矩形公式：?-≈b a a b a f dx x f ))(()( (2) 右矩形公式：))(()(a b b f dx x f b a -≈? (3) 中矩形公式：?-+≈b a a b b a f dx x f ))(2 ( )( 解：(1) )()(a f x f ≈， )()()()(a b a f dx a f dx x f b a b a -=≈?? (2) )()(b f x f ≈，??-=≈b a b a a b a f dx b f dx x f ))(()()( )()(2 1)()()()(2 ηηξf a b dx b x f dx b x f b a b a '--=-'=-'=??，),(,b a ∈ηξ (3) 法1 )2 ( )(b a f x f +≈ ， 法2 可以验证所给公式具有1次代数精度。作一次多项式 )(x H 满足 )2()2( b a f b a H +=+，)2 ()2(b a f b a H +'=+'，则有 2 )2 )((!21)()(b a x f x H x f +-''= -ξ， ),(b a ∈ξ 于是 2．考察下列求积公式具有几次代数精度： （1） ?'+ ≈1 )1(2 1 )0()(f f dx x f ； （2） )3 1()31()(1 1f f dx x f +- ≈?-。 解： （1）当1)(=x f 时，左=1，右=1+0=1，左=右； 当x x f =)(时，左21= ，右=2 1 210=+，左=右； 当2 )(x x f =时，左=3 1 ，右=1，左≠右，代数精度为1。

数值分析作业答案(第5章)

5.1.设A 是对称矩阵且011≠a ，经过一步高斯消去法后，A 约化为 ?? ????21 110 A a a T 证明2A 是对称矩阵。 证明 由消元公式及A 的对称性，有 ,,,3,2,,)2(111 11111 )2(n j i a a a a a a a a a a ji i j ji j i ij ij ==-=- = 故2A 对称。 5.2.设n ij a A )(=是对称正定矩阵，经过高斯消去法一步后，A 约化为 ?? ????21 110 A a a T 其中1)2(2)(-=n ij a A 。证明： (1).A 的对角元素;,,2,1,0n i a ii => (2).2A 是对称正定矩阵。 证明 (1).因为A 对称正定，所以 n i e Ae a i i ii ,,2,1,0),( =>=， 其中T i e )0,,0,1,0,,0( =为第i 个单位向量。 (2).由A 的对称性及消元公式，有 ,,,3,2,,)2(111 11111 )2(n j i a a a a a a a a a a ji i j ji j i ij ij ==-=- = 故2A 也对称。 又由A L A a a T 121110=????? ?，其中

??? ?????- =? ????? ? ?????????--=-111 1 11111 21101 1011n n I a a a a a a L ， 可见1L 非奇异，因而对任意0≠x ，由A 的正定性，有 ,0),(),(,011111>=≠x AL x L x AL L x x L T T T T 故T AL L 11正定。 由,000110211 111121111 1?? ? ?? ?=????????-??????=-A a I a a A a a AL L n T T T 而011>a ，故知2A 正定

研究生《数值分析》练习题

硕士研究生 《数值分析》练习题 一、判断题 1、用Newton 切线法求解非线性线性方程可以任选初值。 （ ） 2、求解非线性线性方程，Newton 切线法比弦截法迭代次数多。 （ ） 3、若n n A R ?∈非奇异，用Jacobi 迭代法求解线性方程组Ax b =必收敛。（ ） 4、Lagrange 插值法与Newton 插值法得到同一个插值多项式。 （ ） 二、填空题 1、近似数 3.14108937a =关 于π具 位有效数字。 2、双点弦截法具有 阶收敛速度。 3、求方程x x e =根的单点弦截法迭代公式是 。 4、设2112A ?? = ? ?? ? ，则()A ρ= 。 5、若(),0,1,2,3i l x i =是以01231,3,,x x x x ==为插值节点的Lagrange 插值基函数，则()()3 3012i i i x l =-=∑ 。 6、由下数据表确定的代数插值多项式的不超过 次。 7、若()8754321f x x x x =+-+，则差商[]0,1,2,,8f = 。 8、拟合三点()()()0,1,1,3,2,2A B C 的 直线是y = 。 三、分析与计算题 1、设()14,2,3515T A x -??==-?? -?? ，求∞=,2,1,,p x A p p 和()1A cond 。

2、1001012,20253A x -???? ? ? == ? ? ? ?-???? ，试计算p p x A ,，p=1,2,∞，和1)(A c o n d 。 3、线性方程组,0Ax b b =≠，用Jacobi 迭代法是否收敛，为什么？其中 122111221A -?? ?=-- ? ?--?? 。 4、线性方程组,0Ax b b =≠，用Jacobi 迭代法是否收敛，为什么？其中 2-11=11111-2A ?? ???? ???? 。 5、已知函数表如下: ⑴ ()111.75ln11.75L ≈、估计截断误差并说明结果有几位有效数字； ⑵ ()211.75ln11.75N ≈、估计截断误差并说明结果有几位有效数字。 6、已知函数表 如下: ⑴用Lagrange 插值法求ln 0.55的近似值()10.55N 、估计截断误差并说明结果的有效数字； ⑵用Newton 插值法求ln 0.55的近似值()20.55N 、估计截断误差并说明结果的有效数字。 7、已知数据如下，求满足条件的Hermite 插值多项式。

数值分析习题集及答案

数值分析习题集 （适合课程《数值方法A》和《数值方法B》） 长沙理工大学 第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2％,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少? 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差? 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求? 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝? 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而 相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大?这个计算过程 稳定吗? 12.计算,取,利用下列等式计算,哪一个得到的结果最好? 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 计算,求对数时误差有多大? 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠? 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足 第二章插值法 1.根据定义的范德蒙行列式,令 证明是n次多项式,它的根是,且 .

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.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函 数表的步长应取多少? 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是 并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误 差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误

数值分析习题

第一章 绪论 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5 105.0-?,那么近似数0.003400有几位有效数字？（有效数字的计算） 2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算） 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5* =，已知 cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2 π=的绝对误差限与相对误差 限。（误差限的计算） 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大？（函数误差的计算） 8 设? -=1 1dx e x e I x n n ，求证： （1）)2,1,0(11 =-=-n nI I n n （2）利用（1）中的公式正向递推计算时误差逐步增大；反向递推计算时误差逐步减小。（计算方法的比较选择）

第二章 插值法 习题主要考察点：拉格朗日插值法的构造，均差的计算，牛顿插值和埃尔米特插值构造，插值余项的计算和应用。 1 已知1)2(,1)1(,2)1(===-f f f ，求)(x f 的拉氏插值多项式。（拉格朗日插值） 2 已知9,4,10=== x x x y ，用线性插值求7的近似值。（拉格朗日线性插值） 3 若),...1,0(n j x j =为互异节点，且有 ) ())(())(()())(())(()(11101110n j j j j j j j n j j j x x x x x x x x x x x x x x x x x x x x x l ----------= +-+- 试证明 ),...1,0()(0 n k x x l x n j k j k j =≡∑=。 （拉格朗日插值基函数的性质） 4 已知352274.036.0sin ,333487.034.0sin ,314567.032.0sin ===，用抛物线插值计 算3367.0sin 的值并估计截断误差。（拉格朗日二次插值） 5 用余弦函数x cos 在00=x ，4 1π =x ，2 2π = x 三个节点处的值，写出二次拉格朗日插值 多项式, 并近似计算6 cos π 及其绝对误差与相对误差，且与误差余项估计值比较。（拉格朗 日二次插值） 6 已知函数值212)6(,82)4(,46)3(,10)1(,6)0(=====f f f f f ，求函数的四阶均差 ]6,4,3,1,0[f 和二阶均差]3,1,4[f 。（均差的计算） 7 设)())(()(10n x x x x x x x f ---= 求][1,0p x x x f 之值，其中1+≤n p ，而节点 )1,1,0(+=n i x i 互异。（均差的计算） 8 如下函数值表 建立不超过三次的牛顿插值多项式。（牛顿插值多项式的构造） 9求一个次数小于等于三次多项式)(x p ，满足如下插值条件：2)1(=p ，4)2(=p ， 3)2(='p ，12)3(=p 。（插值多项式的构造）