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

第一章习题解答

3、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。 解：设=()u f x ，

()()()()()()||||||||||()||()||||()||()||||r r r x e u df x e x df x e x e u u dx u dx u x df x x df x x e x x dx u dx u δ=

≈==≤

()||10.2(())||()||ln ln ln r r r r df x x x x f x x x dx u x x x x

δδδδ==??==

4、长方体的长宽高分别为50cm ，20cm 和10cm ，试求测量误差满足什么条件时其表面积的误差不超过1cm 2。

解：设2()S xy yz zx =++

{}[]{}(,,)(,,)(,,)()||()||()||()(,,)(,,)(,,)||||||max (),(),()2()2()2()max (),(),()1

S x y z S x y z S x y z e S e x e y e z x y z

S x y z S x y z S x y z e x e y e z x y z y z z x x y e x e y e z ???≤++????????≤++ ??????

=+++++< {}[]11max (),(),()2()2()2()4()

110.0031254(502010)320

e x e y e z y z z x x y x y z <=+++++++===++

测量误差小于0.00625时其表面积的误差不超过1cm 2。

7

、计算61)

1.414≈。利用以下四种计算格式，试问哪一种算法误差最小。

（1

（2

）3(3- （3

（4

）99- 解：计算各项的条件数'()(())||()

xf x cond f x f x = 111.41461(),(())| 3.5147(1)

x f x c o n d f x x ===+ 3221.414()(32),(())|49.

3256x f x x c o n d f x ==-=

331.41431(),(())| 1.4557

(32)x f x c o n d f x x ===+ 441.414()9970,(())|4949

x f x x c o n d f x ==-= 由计算知，第三种算法误差最小。

n n

118 ∞=∑、考虑无穷级数，它是微积分中的发散级数。在计算机上计算该级数的部分和，会得到怎样的结果？为什么？

解：在计算机上计算该级数的是一个收敛的级数。因为随着n 的增大，会出现大数吃小数的现象。

9、 通过分析浮点数集合F=（10，3，-2，2）在数轴上的分布讨论一般浮点数集的分布情

况。

解：浮点数集合F=（10，3，-2，2）在数轴上离原点越近，分布越稠密；离原点越远，分布越稀疏。一般浮点数集的分布也符合此规律。

11、为减少乘除法运算次数，应将下面算式怎样改写？

3217151318)()(-+-+-+

=x x x y 解：令1

1-=x u ，则18357+++=u u u y ))(( 12、试推导求函数值)(x f y =的绝对误差和相对误差。 解：**(())()()'()()'()()e f x f x f x f x x f x e x ξ=-=-≈

***()()'()()'()()'()(())()()()()()

r r f x f x f x x f x x x f x e f x x x e x f x f x f x x f x ξ---==≈= 13、试推导求函数值),(y x f 的条件数)),((y x f Cond 。 解：(,)(,)(,)(,)()()f x y f x y f x y f x y x x y y x y ??=+-+-+?

?

''(,)(,)(,)(,)()()(,)(,)(,)y x yf x y xf x y f x y f x y x x y y f x y f x y x f x y y

---≈+ ''''(,)(,)(,)(,)()()(,)(,)(,)(,)(,)()()max ,(,)(,)y x y x yf x y xf x y f x y f x y x x y y f x y f x y x f x y y

yf x y xf x y x x y y f x y f x y x y ---≤+????--??≤ +??? ??????? '(,)'(,)((,))(,)(,)

x y xf x y yf x y Cond f x y f x y f x y ∴=+

14、序列{}n x 满足递推公式

??????=-=≈=+)

,,(.1057321310n n x x x n n 求计算到20x 的的误差，并讨论计算过程的稳定性。 解：000x x e -=

000111555e x x x x e =-=-=

021*********e x x x x e =---=-=)()( …

9020201084545?≈=.e e

误差逐渐增大，计算不稳定

15、写出下面Matlab 程序所描述的数学表达式。

（1）for j=1: n

for i=1: m

y ( i )= A (i, j)*x(j)+y(i)

end

end

(2)for j=1:n

y=x(j)*A(:,j)+y

end

解：（1）m n n m R y R x R A y Ax y ∈∈∈+=?,,,

（2）m n n m R y R x R A y Ax y ∈∈∈+=?,,,

数值分析最佳习题(含答案)

第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5 105.0-?,那么近似数有几位有效数字？（有效数字的计算） 解：2*103400.0-?=x ，325* 102 1 1021---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 解：10314159.0?= π，欲使其近似值* π具有4位有效数字，必需 41*1021 -?≤-ππ，3*3102 11021--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算） 解：3* 1021-?≤ -a a ，2*102 1 -?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤-+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2 123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤-+-≤-b b a a a b b a ab 故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 解：已知 δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5* =，已知 cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差限与相对误差 限。（误差限的计算） 解： * 2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 π ππ252.051.02052)5,20(),(2=??+????≤-v r h v

数值分析试题及答案汇总

数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则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 )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )

数值分析习题与答案

第一章绪论 习题一?1．设x＞0，x＊的相对误差为δ,求f(ｘ)=ln x的误差限。 解:求lnx的误差极限就是求f(x)=lnx的误差限，由公式(１.2.４）有 已知x*的相对误差满足，而 ,故 即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字，并给出其误差限与相对误差限。 解:直接根据定义和式（１.2.２)(1.２.3）则得?有５位有效数字，其误差限,相对误差限 有２位有效数字， 有5位有效数字， 3．下列公式如何才比较准确？ (１)?(2） 解:要使计算较准确，主要是避免两相近数相减，故应变换所给公式。

（1)?(2） 4.近似数x*=０.０３１0,是 3 位有数数字。 5.计算取,利用 :式计算误差最小。 四个选项: 第二、三章插值与函数逼近 习题二、三 1. 给定的数值表 用线性插值与二次插值计算ln０．5４的近似值并估计误差限. 解：仍可使用n＝1及n＝２的Laｇrange插值或Neｗｔo ｎ插值，并应用误差估计（5.8)。线性插值时,用0．５及0．6两点，用Newton插值??误差限 ，因,

故? 二次插值时，用0.5，0．6,０.7三点,作二次Ｎewton插值 ?误差限，故? 2. 在－4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过,函数表的步长ｈ应取多少? 解：用误差估计式（５．8), ?令 因?得 3. 若,求和.

解：由均差与导数关系 ?于是 4. 若互异,求 的值,这里p≤n+1. 解:,由均差对称性 可知当有?而当P＝n ＋1时 ?于是得 5. 求证. 解:解:只要按差分定义直接展开得 ? 6. 已知的函数表

华南理工大学数值分析试题-14年下-C

华南理工大学研究生课程考试 《数值分析》试卷C （2015年1月9日） 1. 考前请将密封线内各项信息填写清楚； 所有答案请按要求填写在本试卷上； 课程代码：S0003004； 4. 考试形式：闭卷； 5. 考生类别：硕士研究生； 本试卷共八大题，满分100分，考试时间为150分钟。 一、（12分）解答下列问题： 1）设近似值0x >，x 的相对误差为δ，试证明ln x 的绝对误差近似为δ。 2）利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值（须写出计算形式），并统计乘法次数。 （12分）解答下列问题： 1）设()235f x x =+，求[]0,1,2f 和[]0,1,2,3f 。 2）利用插值方法推导出恒等式： 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。

（1）设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族，其中1)(0=x q ，求1()q x 和2()q x 。 （2）求形如2y a bx =+的经验公式，使它与下列数据拟合： 四、（14分）对积分()10I f x dx = ?，试 （1）构造一个以012113,,424 x x x ===为节点的插值型求积公式； （2）指出所构造公式的代数精度； （3）用所得数值求积公式计算积分1 203x dx ?的精确值； （4）指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。

（1）设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 （2）用列主元Gauss 消去法解方程组： 12312315410030.112x x x ????????????=????????????-?????? 六、（13分）对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? （11220a a ≠ ） （1）证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散； （2）当同时收敛时，试比较它们的收敛速度。

数值分析习题

习题1 1． 填空题 (1) 为便于算法在计算机上实现，必须将一个数学问题分解为 的 运算； (2) 在数值计算中为避免损失有效数字，尽量避免两个 数作减法运算；为避免 误差的扩大，也尽量避免分母的绝对值 分子的绝对值； (3) 误差有四大来源,数值分析主要处理其中的 和 ； (4) 有效数字越多,相对误差越 ； 2. 用例1.4的算法计算10,迭代3次,计算结果保留4位有效数字. 3. 推导开平方运算的误差限公式,并说明什么情况下结果误差不大于自变量误差. 4. 以下各数都是对准确值进行四舍五入得到的近似数,指出它们的有效数位、误差限和相对误差限. 95123450304051104000003346087510., ., , ., .x x x x x -==?===? 5. 证明1.2.3之定理1.1. 6. 若钢珠的的直径d 的相对误差为1.0%,则它的体积V 的相对误差将为多少。（假定钢珠为标准的球形） 7. 若跑道长的测量有0.1%的误差,对400m 成绩为60s 的运动员的成绩将会带来多大的误差和相对误差. 8. 为使20的近似数相对误差小于0.05%,试问该保留几位有效数字. 9. 一个园柱体的工件,直径d 为10.25±0.25mm,高h 为40.00±1.00mm,则它的体积V 的近似值、误差和相对误差为多少． 10 证明对一元函数运算有 r r xf x f x k x k f x εε'≈= () (())(),() 其中 并求出157f x x x ==()tan ,.时的k 值，从而说明f x x =()tan 在2 x π ≈时是病态问题． 11. 定义多元函数运算 1 1 1,,(),n n i i i i i i S c x c x εε====≤∑∑其中 求出S ε()的表达式，并说明i c 全为正数时，计算是稳定的，i c 有正有负时，误差难以控制． 12. 下列各式应如何改进,使计算更准确：

数值分析复习题要答案

第一章 1、ln2＝0.69314718…，精确到 10－3 的近似值是多少？ 解 精确到 10－3＝0.001，即绝对误差限是 e ＝0.05％，故至少要保留小数点后三位才可以。 ln2≈0.693。 2、设115.80,1025.621≈≈x x 均具有5位有效数字，试估计由这些数据计算21x x ， 21x x +的绝对误差限 解：记126.1025, 80.115x x == 则有11232411 10, | 102|||2 x x x x --≤?-≤?- 所以 121212121212211122||||||||||||x x x x x x x x x x x x x x x x x x -=-+-+≤-- 3411 80.11610 6.10102522 0.007057-==??+≤?? 1212112243|()|||11 |10100.0005522 |x x x x x x x x --≤≤?+?=+-+-+- 3、一个园柱体的工件,直径d 为10.250.25mm,高h 为40.00 1.00mm,则它的体 积V 的近似值、误差和相对误差为多少。 解： ()() 22222222 4 314210254000000330064 221025400002510251002436444 3300624362436 0073873833006 , .....; ()()()......, ..().()..% .r d h V d h V mm d h V dh d d h V mm V V V πππππεεεεε= ≈=??===+=???+?==±====第二章： 1、分别利用下面四个点的Lagrange 插值多项式和Newton 插值多项式N 3(x )， 计算L 3(0.5)及N 3(-0.5) x －2 －1 0 1 f (x ) －1 1 2

数值分析习题集及答案[1].(优选)

数值分析习题集 （适合课程《数值方法A 》和《数值方法B 》） 长沙理工大学 第一章 绪 论 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. 利用公式(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 -=( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx 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. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若

数值分析习题集及答案Word版

数值分析习题集 （适合课程《数值方法A 》和《数值方法B 》） 长沙理工大学 第一章 绪 论 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. 利用公式(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. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx 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. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?

(完整版)数值分析第7章答案

第七章非线性方程求根 一、重点内容提要 （一）问题简介 求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为 函数()f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法 设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在 (a,b)内仅有一个根.令00,a a b b ==,计算0001 ()2x a b =+和0()f x .若0()0f x =则*x x =,结束计算;若00()()0f a f x >,则令10,1a x b b ==,得新的有根区间11[,]a b ;若 00()()0 f a f x <,则令 10,10 a a b x ==,得新的有根区间 11[,]a b .0011[,][,]a b a b ?,11001()2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得 出新的有根区间22[,] a b ,如此反复进行,可得一有根区间套 1100...[,][,]...[,] n n n n a b a b a b --????

数值分析1-4习题及答案

1、 0.1％,要取几位有效数字？ （ c ） (a) 2 (b) 3 (c) 4 (d) 5 2、若* 12.30x =是经过四舍五入得到的近似数，则它有几位有效数字？ （ c ） (a) 2 (b) 3 (c) 4 (d) 5 3、已知n ＋1个互异节点(x 0,y 0), (x 1,y 1),…, (x n ,y n )和过这些点的拉格朗日插值基函数l k (x )(k =0,1,2,…,n )，且ω(x )=(x －x 0) (x －x 1)… (x －x n )．则n 阶差商f (x 0,x 1,…, x n )= ( ) (a) ∑=n k k k y x l 0 )( (b) ∑='n k k k k x l y 0)( (c) ∑=n k k k x y 0)(ω (d) ∑='n k k k x y 0)(ω 4、已知由数据（0，0），（0.5，y ），（1，3），（2，2）构造出的三次插值多项式 33()6 P x x y 的 的系数是，则 等于 （ ） (a) -1.5 (b) 1 (c) 5.5 (d) 4.25 5、设(0,1,2,3,4)i x i =为互异结点，()i l x 为拉格朗日插值基函数,则 4 2 () ()i i i x x l x =-∑等于 （ a ） (a) 0 (b) 1 (c) 2 (d) 4 4()[,],()()(),()(),( )(), ' () ' (),22 ()()_________________________f x C a b H x a b a b H a f a H b f b H f H a f a f x H x ∈++====-=设是满足下列插值条件的三次多项式：则插值余项 1、 是以0，1，2为节点的三次样条函数，则b=-2,c=3 2、 已知(1)0,(1)3,(2)4,f f f =-=-=写出()f x 的牛顿插值多项式 2()P x =___2537 623x x +-__,其余项表达式 R(x)=__() (1)(1)(4) [1,4]6 f x x x ξξ'''-+-∈-_______________________ 3、 确定求积公式1 0121 ()(1)(0)'(1)f x dx A f A f A f -≈-++? 中的待定参数，使其代数精度 尽量高,则A 0=_ 29__________, A 1=__169________, A 2=_29 _______,代数精度=__2_________。

数值分析习题

习题一 1.1 求下列各数的具有四位有效数字的近似值, 并指出其绝对误差限和相对误差限 )1.0ln(,121,101 1,1014321== = = x x x x 1.2 下列各数都是对准确值进行四舍五入得到的近似值, 指出它们的绝对误差限、相对误差限和有效数字的位 数。 3 * 5* 4* 3* 2* 1100.5,5000,50.31,3015.0,0315.0?=====x x x x x 1.3 为了使 3 1的近似值的相对误差不超过0.1%, 问应取几位有效数字? 1.4 怎样计算下列各题才能使得结果比较精确? (1) x x sin )sin(-+ε,其中ε充分小 (2) ? ++1 2 1N N x dx ,其中N 是充分大的正数 (3) x x sin cos 1-,其中x 充分小 (4) o 1cos 1- (5) 1001.0-e (6) )11010ln(84-- 1.5 求方程01562=+-x x 的两个根, 使至少具有四位有效数字。 习题二 2.1 证明方程043 =-+x x 在区间[1,2]内有且仅有一个根。如果用二分法求它具有五位有效数字的根，试问需对 分多少次？(不必求根) 2.2 用二分法求方程0134 =+-x x 在[0.3, 0.4]内的一个根, 精度要求2 10 2 1-?= ε。 2.3 找出下列方程的有根区间，选择适当的初始点用二分法求方程的根，精度要求2 10 -=ε。 (1) 02 =--x x ； (2) 06cos 2 =-++-x e x x ； (3) 01tan =--x x ； (4) 0sin 2=--x e x 。 2.4 考虑方程032 =-x e x ，将其改写为3 x e x ± =，取00=x ，用两种迭代公式迭代，分别收敛到1.0和-0.5附 近的两个根(取精度要求3 10-=ε)。

数值分析整理版试题及答案

数值分析整理版试题及答案

例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x 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 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为

[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f 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 ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P 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 x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有

数值分析试题1

数值分析试卷1 一、填空题（每空2分，共30分） 1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字； 2. 设)(x f 可微，求方程)(x f x =根的牛顿迭代格式是_______________________________________________； 3. 对1)(3++=x x x f ，差商=]3,2,1,0[f _________________； =]4,3,2,1,0[f ________； 4. 已知??? ? ??-='-=1223,)3,2(A x ，则=∞||||Ax ________________，=)(1A Cond ______________________ ； 5. 求解线性方程组?????=+=+045 11532121x x x x 的高斯—赛德尔迭代格式为_______________________________________；该迭代格式迭代矩阵的谱半径=)(G ρ_______________； 二、（12分）(1)设LU A =，其中L 为下三角阵，U 为单位上三角阵。已知 ?????? ? ??------=2100121001210012A ，求L ，U 。 (2)设A 为66?矩阵，将A 进行三角分解：LU A =，L 为单位下三角阵，U 为上三角阵，试写出L 中的元素65l 和U 中的元素56u 的计算公式。 三、给定数据表如下 x 0.20.40.60.81 1.2f(x)212523202124 （1） 用三次插值多项式计算f ( 0.7 ) 的近似值； （2） 用二次插值多项式计算f ( 0.95 ) 的近似值： （3） 用分段二次插值计算 f ( x ) )2.12.0(≤≤x 的近似值能保证有几位有

数值分析习题

第一章 绪论 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 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 1 dx 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 。（插值多项式的构造）

东南大学_数值分析_第七章_偏微分方程数值解法

第七章 偏微分方程数值解法 ——Crank-Nicolson 格式 ****（学号） *****（姓名） 上机题目要求见教材P346,10题。 一、算法原理 本文研究下列定解问题(抛物型方程) 22(,) (0,0)(,0)() (0) (0,)(), (1,)() (0)u u a f x t x l t T t x u x x x l u t t u t t t T ?αβ???-=<<≤≤???? =≤≤??==<≤?? (1) 的有限差分法，其中a 为正常数，,,,f ?αβ为已知函数，且满足边界条件和初始条件。关于式(1)的求解，采用离散化方法，剖分网格，构造差分格式。其中，网格剖分是将区域{}0,0D x l t T =≤≤≤≤用两簇平行直线 (0) (0)i k x x ih i M t t k k N τ==≤≤?? ==≤≤? 分割成矩形网格，其中,l T h M N τ==分别为空间步长和时间步长。将式(1)中的偏导数使用不同的差商代替，将得到不同的差分格式，如古典显格式、古典隐格式、Crank-Nicolson 格式等。其中，Crank-Nicolson 格式具有更高的收敛阶数，应用更广泛，故本文采用Crank-Nicolson 格式求解抛物型方程。 Crank-Nicolson 格式推导：在节点(,)2 i k x t τ +处考虑式(1)，有 22(,)(,)(,)222 i k i k i k u u x t a x t f x t t x τττ??+-+=+?? (2) 对偏导数 (,)2 i k u x t t τ ?+?用中心差分展开 []2311+13 1(,)(,)(,)(,) ()224k k i k i k i k i i k i k u u x t u x t u x t x t t t t ττηητ++??+=--<

数值分析复习题及答案

数值分析复习题 一、选择题 1. 3.142和3.141分别作为π的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式()()2 11211()(2)636f x dx f Af f ≈++?，则A ＝（ ） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（ ） A ．() 00l x ＝0，()110l x = B ． ()00l x ＝0，()111l x = C ．()00l x ＝1，()111l x = D ． ()00l x ＝1，()111l x = 4. 设求方程()0f x =的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=??++=??--=? 作第一次消元后得到的第3个方程（ ）. A ．232x x -+= B ．232 1.5 3.5x x -+= C ．2323x x -+= D ．230.5 1.5x x -=- 二、填空 1. 设 2.3149541...x *=，取5位有效数字，则所得的近似值x= . 2.设一阶差商 ()()()21122114,321f x f x f x x x x --= ==---， ()()()322332615,422f x f x f x x x x --===--

则二阶差商 ()123,,______f x x x = 3. 设(2,3,1)T X =--, 则2||||X = ，=∞||||X 。 4．求方程 2 1.250x x --= 的近似根，用迭代公式 1.25x x =+，取初始值 01x =， 那么 1______x =。 5．解初始值问题 00'(,)()y f x y y x y =??=?近似解的梯形公式是 1______k y +≈。 6、 1151A ??= ?-??，则A 的谱半径 ＝ 。 7、设 2()35, , 0,1,2,... , k f x x x kh k =+== ，则[]12,,n n n f x x x ++= 和[]123,,,n n n n f x x x x +++= 。 8、若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都 。 9、解常微分方程初值问题的欧拉（Euler ）方法的局部截断误差为 。 10、为了使计算 23123101(1)(1)y x x x =+ +----的乘除法运算次数尽量的少,应将表达式改写 成 。 11. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 12. 一阶均差()01,f x x = 13. 已知3n =时，科茨系数()()()33301213,88C C C ===，那么 ()33C = 14. 因为方程()420x f x x =-+=在区间[]1,2上满足 ，所以()0f x =在区间内有根。 15. 取步长0.1h =，用欧拉法解初值问题()211y y y x y ?'=+???=?的计算公式 . 16.设 * 2.40315x =是真值 2.40194x =的近似值，则*x 有 位有效数字。

数值分析题库答案

1. 正方形的边长大约为100cm ，应怎样测量才能使面积误差不超过1cm 2? 2. 已测得某场地长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的计算结果能有几位有效数字？

6. 已知33348 7.034.0sin ,314567.032.0sin ==，请用线性插值计算3367.0sin 的值，并估计截断误差. 7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并估计误差. 8. 已知 1 6243sin ,sin π ππ== =请用抛物插值求sin50的值,并估计误差

9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x 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?? ???? ，上的三次埃尔米特

数值分析第七章非线性方程求根习题答案

第七章非线性方程求根 （一）问题简介 求单变量函数方程 ()0f x = (7.1) 的根是指求*x （实数或复数），使得(*)0f x =.称*x 为方程(7.1)的根,也称*x 为函数() f x 的零点.若()f x 可以分解为 ()(*)()m f x x x g x =- 其中m 为正整数,()g x 满足()0g x ≠,则*x 是方程(7.1)的根.当m=1时,称*x 为单根;当m>1时,称*x 为m 重根.若()g x 充分光滑,*x 是方程(7.1)的m 重根,则有 (1)() (*)'(*)...(*)0,(*)0m m f x f x f x f x -====≠ 若()f x 在[a,b]上连续且()()0f a f b <,则方程(7.1)在(a,b)内至少有一个实根,称[a,b]为方程(7.1)的有根区间.有根区间可通过函数作图法或逐次搜索法求得. (二)方程求根的几种常用方法 1.二分法 设()f x 在[a,b]上连续,()()0f a f b <,则()0f x =在(a,b)内有根*x .再设()0f x =在(a,b)内 仅有一个根.令00,a a b b ==,计算0001()2x a b =+和 0()f x .若0()0f x =则*x x =,结束计算;若 00()()0 f a f x >,则令 10,1a x b b ==,得新的有根区间 11[,] a b ;若 00()()0 f a f x <,则令 10,10a a b x ==,得新的有根区间11[,]a b .0011[,][,]a b a b ?,11001 () 2b a b a -=-.再令1111 ()2x a b =+计算1()f x ,同上法得出新的有根区间22[,]a b ,如此反复进行,可得一有根区 间套 1100...[,][,]...[,] n n n n a b a b a b --???? 且110011 *,0,1,2,...,()...() 22n n n n n n a x b n b a b a b a --<<=-=-==-. 故 1 lim()0,lim lim ()* 2n n n n n n n n b a x a b x →∞→∞→∞-==+=

数值计算方法第七章习题 2013

计算方法 第七章 习题 复习与思考题 1．设f ∈C [a , b ]，写出三种常用范数2 1 f f 及∞ f 。 2．f , g ∈C [a , b ]，它们的内积是什么？如何判断函数族{? 0, ? 1, …, ? n }∈C [a , b ]在[a ,b ]上线性无关？ 3．什么是函数f ∈C [a , b ]在区[a , b ]上的n 次最佳一致逼近多项式？ 4．什么是f 在[a , b ] 上的n 次最佳平方逼近多项式？什么是数据{}m i f 0的最小二乘曲 线拟合？ 5．什么是[ a , b ]上带权ρ (x )的正交多项式？什么是[ -1, 1 ]上的勒让德多项式？它有什 么重要性质？ 6．什么是切比雪夫多项式？它有什么重要性质？ 7．用切比雪夫多项式零点做插值得到的插值多项式与拉格朗日插值有何不同？ 8．什么是最小二乘拟合的法方程？用多项式做拟合曲线时，当次数n 较大时为什么不直接求解法方程？ 9．哪种类型函数用三角插值比用多项式插值或分段多项式插值更合适？ 10．判断下列命题是否正确？ （1）任何f (x ) ∈C [a , b ]都能找到n 次多项式P n (x ) ∈ H n ，使| f (x ) - P n (x ) | ≤ ε ( ε 为任给的误差限)。 （2）n n H x P ∈)(* 是f (x )在[ a , b ]上的最佳一致逼近多项式，则)()(lim * x f x P n n =∞ →对 ],[b a x ∈?成立。 （3）f (x ) ∈C [a , b ]在[a , b ]上的最佳平方逼近多项式P n (x ) ∈ H n 则)()(lim x f x P n n =∞ →。 （4）)(P ~ x n 是首项系数为1的勒让德多项式，Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 ? ? --1 1 21 1 2d )(d )](P ~ [x x Q x x n n 。 （5）)(T ~ x n 是[-1 , 1]上首项系数为1的切比雪夫多项式。Q n (x ) ∈ H n 是任一首项系数为1的多项式，则 .)(max )(~ max 1 11 1x Q x T n x n x ≤≤-≤≤-≤ （6）当数据量很大时用最小二乘拟合比用插值好。