2.3.5牛顿法的特点及牛顿下山法
牛顿法有以下特点：
● 收敛速度快，达到平方收敛；
● 计算量较大，对函数)(x f 解析性质要求较高。
● 具有局部收敛性，即当初值0x 与*x 太远时，收敛速度很慢，甚至不收敛，仅在*x 附近可达到平方收敛；
将牛顿迭代公式修改为以下形式
)
()
(1k k k k x f x f x x '-=+λ
其中参数)1,0(∈λ称为下山因子。用试算的方法选取),2
1
,21,21,1(32 λ，使得
)()(1k k x f x f <+
称满足上述要求的算法为牛顿下山法。
2.4 弦位法
2.4.1 弦位法的思想
弦位法又称弦截法或割线法。是用差商代替牛顿公式中的微商)(k x f '；或者说是用)(x f 在点))(,(11--k k x f x 和))(,(k k x f x 处的割线的零点作为新的迭代点
)()
()()
(111--+---=k k k k k k k x x x f x f x f x x
2.4.2 弦位法的几何意义 略
2.4.3 弦位法的计算步骤 略
2.4.4 弦位法的收敛性及收敛速度
定理2-7 设*x 是方程0)(=x f 的单根，且)(x f 在*x 的某邻域有二阶连续导数，则∈?10,x x ),(*δx N ，弦位法在*x 附近局部收敛，且按阶1.618收敛到*x 。 2.4.4 弦位法的特点
● 收敛速度夹较快，但比牛顿法慢； ● 超线性收敛，收敛阶为1.618；
● 无需计算导数，每步只需计算一次函数值；
● 属于多点迭代，而牛顿法和一般迭代法属于单点迭代。
2.5 迭代法的改善与加速
2.5.1 一般的加速算法
加速算法思想是先求)(1k k x x ?=+，然后选取合适的数m 和n ，令 k k k nx x m x +=++11 作为下一次迭代点。
为选取合适的线性组合系数m,n 。由于)(x ?连续，且)(**x x ?=，从而 ))(()()(**1*k k k x x x x x x -'=-=-+ξ??? 用a 近似代替上式中的)(ξ?'，立即有 )(*1*k k x x a x x -≈-+ 整理得
k k x a
a
x a x +-+≈
+1111* 从而得到加速迭代公式
k k k x a
a
x a x +-+=+1)(111?
2.5.2 埃特金（Aithen ）加速算法
埃特金算法的思想是将一般迭代法和弦位法相结合来实现加速迭代算法的收敛速度的。
为构造埃特金法加速迭代公式，设k x 是)(x x ?=的一个近似解，首先令
)(1k k x x ?=+，)(?1k k x x
?=+ 然后在曲线)(x y ?=上，过点),(1
+k k
x x P 和)?,(?11
++k k x
x P
做连线，该直线的两点式方程为
k
k k k k k x x x y x x x x
--=
--++++1111? 最后考虑到方程)(x x ?=的解是曲线)(x y ?=和直线x y =的交点*P 的横坐标，
可用弦P
P ?与直线x y =的交点1
+k P 代替*P ,即用点1
+k P 的横坐标作为)(x x ?=的一个近似解，其值可通过将上式中的y 用x 代替并解出x 得到
11211?2?+++++--=
k k k k k k x
x x x x
x x 用上式进一步可以构造出两种具有加速收敛的迭代公式
11211112111?2)?(??22?++++++++++---=+--≈k k k k k k k k k k k k k x x x x x
x x x x x x x x 称上式为埃特金算法的迭代公式；
1
12
1112111
?2)(?22?+++++++++---
=+--≈k k k k k k k k k k k k k x x x x x x x x x x x x x 称上式为Steffensen （斯姬芬森）外推迭代公式。实质上，若令
[]x
x x x x x x +---
=)(2))(()()(2
????ψ，则 ,2,1,0)(1==+k x x k k ψ称为Steffensen 迭代公式。
2.5.3 埃特金算法的几何意义 略 p25图2-7 2.5.4 埃特金算法的计算步骤
略 例8p25，原本发散的迭代公式13
1-=+k k x x 用埃特金算法是收敛的。
2.5.4 加速算法的收敛性
设*x 是)(x x ?=的精确解，由迭代法收敛性定理知：若)(*x ?'在*x 的某个邻域上连续，若1)(0*<''x ?，则迭代算法肯定不收敛。但经过改进后的Aitken 算法超线性收敛，Steffensen （斯基芬森）外推算法若
)(x ?'''在*x 处连续，且1)(*≠'x ?，则算法至少有二次收敛性。
定理2-8 Steffensen 迭代法收敛定理 设函数)(x ψ按上式定义的迭代函数：
(1) 若*x 是)(x ?的不动点，)(x ?'在*x 处连续，且1)(≠'x ?，则*x 也是)(x ψ的不动点，反之，若*x 是)(x ψ的不动点，则*x 也是)(x ?的不动点。
(2) 若*x 是)(x ?的不动点，)(x ?'在*x 处连续，且1)(≠'x ?，则Steffensen 迭代公式至少是平方收敛。（证明见6p60）
2.6 代数方程求根
当)(x f 是多项式时，方程)(x f =0称作代数方程。在牛顿迭代法中，要求
)(),(00x f x f '，设 n n n
n n n a x a x a x a x a a x a x a x a x f ++++=++++=---)))((()(12101110
令 ),,2,1(,1000n i b x a b a b i i i =+==-, 则 n b x f =)(0。
为求)(0x f '，考查)(x f 求导，即
)
))((())))((()))((()(()
))((())))(((]))((([)
))(((]))(([)(1210121022101012101210221012101210--------++++++++
+++++=++++++++'+++=++++'+++='n n n n n n n n a x a x a x a x a x a x a x a x a x a x a x a x a x a a x a x a x a x
a x a x a x a x a x a x a x a a x a x a x a x a x a x a x a x f 则
1
020201001020100010201000202010001000)))((()))((())))((()))((()(()(-----++++=++++++++
+++++='n n n n n b x b x b x b x b a x a x a x a x a x a x a x a x a x a x a x a x a x a x f 上述结果用于代数方程)(x f =0的牛顿迭代公式)(/)(1k k k k x f x f x x '-=+ 中，)(),(k k x f x f '的计算公式如下：
??
?
??=≤≤+==-n k i k i i b
x f n i b x a b a b )(1100 ??
?
??='-≤≤+==--110
0)(11n k i k i i c
x f n i c x b c b c 2.7 求方程的m 重根
设*x 是方程0)(=x f 的)2(≥m 重根，则)()()(*x g x x x f m -=，其中
0)(*≠x g 。若)(x f 在*x 的某个邻域内有m 阶连续导数，由于 0)()()(*)1(**==='=-x f x f x f m ，而0)(*)(≠x f m ，由泰勒展式，得
m m x x m f x f )(!
)
()(*1)(-=ξ
1*2)()()!1()
()(---='m m x x m f x f ξ
2*3)()()!
2()
()(---=''m m x x m f x f ξ
其中321,,ξξξ介于x 与*x 之间，由牛顿迭代函数)
()
()(x f x f x x '-
=?可得 *
2)(1)(**
])
()()([lim )(lim )(**x mf f x x x x x m m x x x x =--==→→ξξ?? m f m f f m x f x f x f x x m m m x x x x x x 1
1)]
([)()()1(lim )]([)()(lim )(lim )(22)(3)(1)(2*
***-=-='''='='→→→ξξξ?? 由于**)(x x =?且1)(0*<'为提高收敛阶数，分两种情况进行修改：
（1） 当重数m 已知时，构造迭代函数)
()
()(1x f x mf x x '-
=?，此时 显然**1)(x x =?，利用上面的结果m x f x f x f x
x 1
1)]
([)()(lim 2
*
-='''→，不难推出 0)(*1
='x ?，从而迭代公式)(11k k x x ?=+至少平方收敛。 （2） 当重数m 未知时，令)()()(x f x f x u '==)
()()()
()(*
*x g x x x mg x g x x '-+-， 显然，*x 是)(x u 的单零点。令迭代函数)
()
()(2x u x u x x '-
=?，此时， 0)]([)
()(lim )(lim )(2
*2
**
='''='='→→x u x u x u x x x x x x ??，从而迭代公式
)
()()]([)
()()()()(221k k k k k k k k k k k x f x f x f x f x f x x u x u x x x ''-''-
='-
==+?至少平方收敛。 P27例9 作业：p28 3,4,9,14,15,18
数值分析参考答案(第二章)doc资料
数值分析参考答案(第 二章)
第二章 插值法 1．当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)23537623l x l x x x x x x x =-+=---+-+=+- 2．给出()ln f x x =的数值表 用线性插值及二次插值计算ln0.54的近似值。 解：由表格知， 01234012340.4,0.5,0.6,0.7,0.8;()0.916291,()0.693147()0.510826,()0.356675()0.223144 x x x x x f x f x f x f x f x ======-=-=-=-=- 若采用线性插值法计算ln0.54即(0.54)f ， 则0.50.540.6<<
2 112 1 221 11122()10(0.6)()10(0.5)()()()()() x x l x x x x x x l x x x x L x f x l x f x l x -==----= =---=+ 6.93147(0.6) 5.10826(0.5)x x =--- 1(0.54)0.62021860.620219L ∴=-≈- 若采用二次插值法计算ln0.54时， 1200102021101201220212001122()() ()50(0.5)(0.6) ()() ()() ()100(0.4)(0.6) ()()()() ()50(0.4)(0.5) ()() ()()()()()()() 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 L x f x l x f x l x f x l x --==------==-------= =----=++ 500.916291(0.5)(0.6)69.3147(0.4)(0.6)0.51082650(0.4)(0.5) x x x x x x =-?--+---?--2(0.54)0.615319840.615320L ∴=-≈- 3．给全cos ,090x x ≤≤的函数表，步长1(1/60),h '==若函数表具有5位有效数字，研究用线性插值求cos x 近似值时的总误差界。 解：求解cos x 近似值时，误差可以分为两个部分，一方面，x 是近似值，具有5位有效数字，在此后的计算过程中产生一定的误差传播；另一方面，利用插值法求函数cos x 的近似值时，采用的线性插值法插值余项不为0，也会有一定的误差。因此，总误差界的计算应综合以上两方面的因素。 当090x ≤≤时， 令()cos f x x = 取0110,( )606018010800 x h ππ ===?=
计算方法第二章方程求根上机报告
实验报告名称 班级：学号：姓名：成绩： 1实验目的 1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点。 2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法。 2 实验内容 用牛顿法和割线法求下列方程的根 x^2-e^x=0; x*e^x-1=0; lgx+x-2=0; 3实验步骤 1)根据二分法和牛顿迭代法，割线法的算法编写相应的求根函数； 2)将题中所给参数带入二分法函数，确定大致区间； 3)用牛顿迭代法和割线法分别对方程进行求解； 3 程序设计 牛顿迭代法x0=1.0; N=100; k=0; eps=5e-6; delta=1e-6; while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if k>N disp('Newmethod failed')
break end if(abs(x1-x0)=delta) c=x1; x1=cutnext(x0,x1); x0=c; %x0 x1μYí?μ?μ?x1 x2 è?è?±￡′??úx0 x1 end k=k+1; if k>N disp('Cutline method failed') break; end if(abs(x1-x0)数值分析第二章小结
第2章线性方程组的解法 --------学习小结 一、本章学习体会 通过本章知识的学习我首先了解到求解线性方程组的方法可分为两类：直接法和迭代法。计算机虽然运行速度很快，但面对运算量超级多的问题，计算机还是需要很长的时间进行运算，所以，确定快捷精确的求解线性方程组的方法是非常必要的。 本章分为四个小节，其中前两节Gauss消去法和直接三角分解法因为由之前《线性代数》学习的一定功底，学习起来还较为简单，加之王老师可是的讲解与习题测试，对这一部分有了较好的掌握。第三节矩阵的条件数与病态方程组，我 Ax 的系数矩阵A与左端向量b的元素往往是通首先了解到的是线性方程组b 过观测或计算而得到，因而会带有误差。即使原始数据是精确的，但存放到计算机后由于受字长的限制也会变为近似值。所以当A和b有微小变化时，即使求解过程精确进行，所得的解相对于原方程组也可能会产生很大的相对误差。对于本节的学习掌握的不是很好，虽然在课后习题中对课堂知识有了一定的巩固，但整体感觉没有很好的掌握它。第四节的迭代法，初次接触迭代法，了解到迭代法就是构造一个无线的向量序列，使他的极限是方程组的解向量。迭代法应考虑收敛性与精度控制的问题。三种迭代方法的基本思想我已经掌握了，但是在matlab 的编程中还存在很大的问题。 在本节的学习中我认为我最大的问题还是程序的编写。通过这段时间的练习，虽然掌握了一些编写方法和技巧。相比于第一章是对其的应用熟练了不少，但在程序编写上还存在很多问题。希望在以后的学习中能尽快熟练掌握它，充分发挥它强大的作用。 二、本章知识梳理 2.1、Gauss消去法（次重点） Gauss消去法基本思想：由消元和回代两个过程组成。 a(k=1,2,```,n-1)均不为零的充分必要条件定理顺序Gauss消去法的前n-1个主元素)(k kk 是方程组的系数矩阵A的前n-1个顺序主子式
数值分析课后题答案
数值分析 第二章 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 ''∴= --
数值分析习题集及答案
（适合课程《数值方法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 =按递推公式 11783 100 n n Y Y -=- ( n=1,2,…) 计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加, 而相对误差却减小. 11. 序列{}n y 满足递推关系1101 n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 (21)f =-,取 2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 6 3 11,(322), ,9970 2. (21) (322) --++ 13. 2 ()ln(1)f x x x =- -,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等 价公式 2 2 ln(1)ln(1)x x x x - -=-+ + 计算,求对数时误差有多大? 14. 试用消元法解方程组{ 10 10 12121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin , 2 s ab c = 其中c 为弧度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证 明面积的误差s ?满足 . s a b c s a b c ????≤ ++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令
数值分析第四版习题及答案
第四版 数值分析习题 第一章 绪 论 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 =- 计算,求对数时误差有多大?
数值分析第二章复习与思考题
第二章复习与思考题 1.什么是拉格朗日插值基函数？它们是如何构造的？有何重要性质？ 答：若n 次多项式()),,1,0(n j x l j =在1+n 个节点n x x x <<< 10上满足条件 (),,,1,0,, ,0, ,1n k j j k j k x l k j =?? ?≠== 则称这1+n 个n 次多项式()()()x l x l x l n ,,,10 为节点n x x x ,,,10 上的n 次拉格朗日插值基函数. 以()x l k 为例，由()x l k 所满足的条件以及()x l k 为n 次多项式，可设 ()()()()()n k k k x x x x x x x x A x l ----=+- 110， 其中A 为常数，利用()1=k k x l 得 ()()()()n k k k k k k x x x x x x x x A ----=+- 1101， 故 ()()()() n k k k k k k x x x x x x x x A ----= +- 1101 ， 即 ()()()()()()()()∏ ≠=+-+---=--------=n k j j j k j n k k k k k k n k k k x x x x x x x x x x x x x x x x x x x x x l 0110110)( . 对于()),,1,0(n i x l i =，有 ()n k x x l x n i k i k i ,,1,00 ==∑=，特别当0=k 时，有 ()∑==n i i x l 0 1. 2.什么是牛顿基函数？它与单项式基{ }n x x ,,,1 有何不同？ 答：称()()()(){ }10100,,,,1------n x x x x x x x x x x 为节点n x x x ,,,10 上的牛顿基函数，利用牛顿基函数，节点n x x x ,,,10 上的n 次牛顿插值多项式()x P n 可以表示为 ()()()()10010---++-+=n n n x x x x a x x a a x P 其中[]n k x x x f a k k ,,1,0,,,,10 ==.与拉格朗日插值多项式不同，牛顿插值基函数在增加节点时可以通过递推逐步得到高次的插值多项式，例如 ()()()()k k k k x x x x a x P x P --+=++ 011,
数值分析第四版习题和答案解析
第四版 数值分析习题 第一章绪论 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的多项式中,求在的最佳一致逼近多项式.
数值分析第五版计算实习题复习进程
数值分析计算实习题 第二章 2-1 程序： clear;clc; x1=[0.2 0.4 0.6 0.8 1.0]; y1=[0.98 0.92 0.81 0.64 0.38]; n=length(y1); c=y1(:); for j=2:n %求差商 for i=n:-1:j c(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1)); end end syms x df d; df(1)=1;d(1)=y1(1); for i=2:n %求牛顿差值多项式 df(i)=df(i-1)*(x-x1(i-1)); d(i)=c(i)*df(i); end disp('4次牛顿插值多项式'); P4=vpa(collect((sum(d))),5) %P4即为4次牛顿插值多项式,并保留小数点后5位数pp=csape(x1,y1, 'variational');%调用三次样条函数 q=pp.coefs; disp('三次样条函数'); for i=1:4 S=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1]; S=vpa(collect(S),5) end x2=0.2:0.08:1.08; dot=[1 2 11 12]; figure ezplot(P4,[0.2,1.08]); hold on y2=fnval(pp,x2); x=x2(dot); y3=eval(P4); y4=fnval(pp,x2(dot)); plot(x2,y2,'r',x2(dot),y3,'b*',x2(dot),y4,'co'); title('4次牛顿插值及三次样条'); 结果如下： 4次牛顿插值多项式 P4 = - 0.52083*x^4 + 0.83333*x^3 - 1.1042*x^2 + 0.19167*x + 0.98 三次样条函数
数值分析课后习题答案
第一章 题12 给定节点01x =-，11x =，23x =，34x =，试分别对下列函数导出拉格朗日插值余项： （1） （1） 3 ()432f x x x =-+ （2） （2） 4 3 ()2f x x x =- 解 （1）(4) ()0f x =， 由拉格朗日插值余项得(4)0123() ()()()()()()0 4!f f x p x x x x x x x x x ξ-=----=； （2）(4) ()4!f x = 由拉格朗日插值余项得 01234! ()()()()()() 4! f x p x x x x x x x x x -= ----(1)(1)(3)(4)x x x x =+---. 题15 证明：对于()f x 以0x ，1x 为节点的一次插值多项式()p x ，插值误差 012 10()()()max () 8x x x x x f x p x f x ≤≤-''-≤. 证 由拉格朗日插值余项得 01() ()()()()2!f f x p x x x x x ξ''-= --，其中01x x ξ≤≤， 01 0101max ()()()()()()()() 2!2!x x x f x f f x p x x x x x x x x x ξ≤≤''''-=--≤-- 01210()max () 8x x x x x f x ≤≤-''≤. 题22 采用下列方法构造满足条件(0)(0)0p p '==，(1)(1)1p p '==的插值多项式 ()p x ： （1） （1） 用待定系数法； （2） （2） 利用承袭性，先考察插值条件(0)(0)0p p '==，(1)1p =的插值多项式 ()p x . 解 （1）有四个插值条件，故设230123()p x a a x a x a x =+++，2 123()23p x a a x a x '=++， 代入得方程组001231123010231 a a a a a a a a a =? ?+++=?? =? ?++=? 解之，得01230 021 a a a a =??=?? =??=-?
数值分析第五版答案(全)
第一章 绪论 1．设0x >,x 的相对误差为δ，求ln x 的误差。 解：近似值*x 的相对误差为***** r e x x e x x δ-= == 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈ 进而有(ln *)x εδ≈ 2．设x 的相对误差为2%，求n x 的相对误差。 解：设()n f x x =，则函数的条件数为'()||() p xf x C f x = 又 1'()n f x nx -=, 1||n p x nx C n n -?∴== 又((*))(*)r p r x n C x εε≈? 且(*)r e x 为2 ((*))0.02n r x n ε∴≈ 3．下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指 出它们是几位有效数字：*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =, *57 1.0.x =? 解：*1 1.1021x =是五位有效数字； *20.031x =是二位有效数字； *3385.6x =是四位有效数字； *456.430x =是五位有效数字； *57 1.0.x =?是二位有效数字。 4．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中**** 1234,,,x x x x 均为第3题所给的数。 解：
*4 1*3 2*13*3 4*1 51 ()102 1()102 1()102 1()102 1()102x x x x x εεεεε-----=?=?=?=?=? ***124***1244333 (1)() ()()() 111101010222 1.0510x x x x x x εεεε----++=++=?+?+?=? ***123*********123231132143 (2)() ()()() 1111.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈ **24****24422 *4 33 5 (3)(/)()() 110.0311056.430102256.43056.430 10x x x x x x x εεε---+≈??+??=?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？ 解：球体体积为343V R π= 则何种函数的条件数为 23'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε=%1
数值分析第五版全答案
第四章 数值积分与数值微分 1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度： 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解： 求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。 （1）若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =，则 1012h A A A -=++ 令()f x x =，则 110A h A h -=-+ 令2 ()f x x =，则 3 221123 h h A h A -=+ 从而解得 01 1431313A h A h A h -?=?? ?=?? ?=?? 令3 ()f x x =，则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 故 101()()(0)()h h f x dx A f h A f A f h --=-++? 成立。 令4 ()f x x =，则
第二章_Volterra_方程的求解
第二章 Volterra 方程的求解 §2.1 第二类Volterra 方程求解 积分方程是近代数学的一个重要分支，它与微分方程、泛函分析、计算数学和有机分析等有着紧密的联系.同时，它也是解决力学、数学物理和工程技术等问题的一种重要工具. 本章首先介绍积分方程的基本概念，其次利用压缩映照原理讨论积分方程的可解性及逐次逼近方法，并扼要介绍Fredholm 定理，讨论一些非线性积分方程的解法. 第二类Volterra 方程一般形式为： ()(,)()()x a x K x t t dt f x ?λ?=+? , (2.1.1) A ． 化为常微分方程求解 例2.1.10().x x t e t dt x ?-=? 解 由0 ()x x t e e t dt x ?-=? ，
得0 ()x t x e t dt xe ?--=?, 求导得(),x x x e x e xe ?---=- 即()1x x ?=-. 例2.1.2 0()().x x x t dt e ??=+? 解 求导得()().x x x e ??'=+ 定解条件0 0(0)() 1.t dt e ??=+=? 化为微分方程 , (0) 1. x e ???'?=+? =? 容易得到()(1)x x x e ?=+. 定理2.1.1 如果第二类Volterra 方程(2.1.1)的核(,)K x t 为()x t -的(1)n -次多项式 01(,)()()() K x t a x a x x t =+-2 2()()2!a x x t +-11()()(1)!n n a x x t n --++ -- , 令 1 1()()()(1)!x n a y x x t t dt n ?-=--?，
第二章 数值分析--方程求根
第二章 方程求根 教学内容： 1．二分法 2．基本迭代法 3．牛顿法 4．弦位法 5．埃特金法和斯基芬森法 6．重根的情况 教学重点： 各种算法的思路及迭代公式的构造 教学难点： 各种算法的收敛性、收敛速度及误差估计 计划学时：5-6学时 授课提纲： 方程求根就是求函数)(x f 的零点*x ，即求解方程 0)(=x f 这里，0)(=x f 可以是代数方程，也可以不是，如超越方程。 方程的根既可以是实数，也可以是复数；既可能是单根，也可能是重根；即可能要求求出给定范围内的某个根，也可能要求求出方程全部的根。 本章介绍的方法对两类方程都适用，但大部分都是要求知道根在什么范围内，且在此范围内只有一个单根。若有α使得0)(,0)(≠'=ααf f ，则称α是方程0)(=x f 的单根；若有α使得 0)(,0)()()()()1(≠==='=-ααααm m f f f f ， 则称α是方程0)(=x f 的m 重根。 设)(x f 在区间[a,b]连续，若0)()(
2.1.2 二分法思想 区间对分，去同存异 2.1.3 二分法计算步骤 步1：令2/)(0b a x +=，计算)(0x f ； 步2：若0)(0=x f ，令0*x x =，计算结束； 步3：若)(0x f *)(a f >0，令0x a =；否则令0x b =； 步4：若ε≤-||a b ，令2/)(*b a x +=,计算结束；否则转步1。 2.1.4 二分法误差分析和收敛性 记第k 次区间中点为k x ，则有 2/)(0*a b x x -≤-，21*2/)(a b x x -≤-，1*2/)(,+-≤-k k a b x x 故当∞→k 时，*x x k →。 为使ε≤-k x x *，解不等式ε≤-+12/)(k a b ，得 12ln /]ln )[ln(---≥εa b k 2.1.5 二分法的优缺点 ● 算法简单直观，易编程计算； ● 只需)(x f 连续即可； ● 区间收缩速率相同，收敛速度慢； ● 无法求复根和偶重根。 例2-1 p15例1 2.2 迭代法 2.2.1 迭代法原理 0)(=x f ? )(x x ?= )(x f 的根 )(x ?的不动点 2.2.2 迭代法思路 任取初值],[0b a x ∈，令)(01x x ?=，)(12x x ?=，反复迭代，即得 ),2,1,0(),(1 ==+k x x k k ? 直到满足精度要求的k x 来近似*x 。称)(x x ?=为迭代公式，)(x ?为迭代函数，{k x }为迭代序列。 若{k x }收敛时，称迭代公式是收敛的。此时设=∞ →k k x lim *x ，当)(x ?连续时 )()lim ()(lim lim *1*x x x x k k k k k ???====∞ →∞ →+∞ → 亦即0)(*=x f 。若{k x }不收敛，称迭代公式是发散的。
数值分析习题集及答案
数值分析习题集 （适合课程《数值方法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次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误
数值分析第五版答案
第一章 绪论 3．下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指 出它们是几位有效数字：*1 1.1021x =,*20.031x =, *3385.6x =, * 456.430x =,*57 1.0.x =? 解：*1 1.1021x =是五位有效数字； *20.031x =是二位有效数字； *3385.6x =是四位有效数字； *456.430x =是五位有效数字； *57 1.0.x =?是二位有效数字。 4．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234 ,,,x x x x 均为第3题所给的数。 解： *4 1* 3 2* 13* 3 4* 1 51 ()1021()1021()1021()1021()102 x x x x x εεεεε-----=?=?=?=?=? *** 124***1244333 (1)()()()() 1111010102221.0510x x x x x x εεεε----++=++=?+?+?=? *** 123*********123231132143 (2)() ()()() 111 1.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈
** 24**** 24422 * 4 33 5 (3)(/) ()() 11 0.0311056.430102256.43056.430 10x x x x x x x εεε---+≈ ??+??= ?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？ 解：球体体积为34 3 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε=?≈ 6．设028Y = ，按递推公式1n n Y Y -= （n=1,2,…） 计算到100Y 27.982≈（5位有效数字），试问计算100Y 将有多大误差？ 解：1n n Y Y -=- 10099Y Y ∴=- 9998Y Y = 9897Y Y =-…… 10Y Y =- 依次代入后，有1000100Y Y =- 即1000Y Y = 27.982, 100027.982Y Y ∴=-
数值分析第四版习题及答案
数值分析第四版习题及答案
第四版 数值分析习题 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****1 2 3 4 5 1.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. 设0 28,Y =按递推公式 11 783100 n n Y Y -=( n=1,2,…) 计算到100Y .若取78327.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字78327.982). 8. 当N 充分大时,怎样求2 11N dx x +∞ +?? 9. 正方形的边长大约为100㎝,应怎样测量才能
使其面积误差不超过1㎝2 ? 10. 设212 S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小. 11. 序列{}n y 满足递推关系1 101n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 21)f =,取2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 63 22)70 2. (21)(322)--++ 13. 2 ()ln(1)f x x x =-,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 2 2 ln(1)ln(1)x x x x -=-+ 计算,求对数时误差有多大? 14. 试用消元法解方程组 {101012121010;2. x x x x +=+=假定只用 三位数计算,问结果是否可靠? 15. 已知三角形面积 1 sin ,2 s ab c = 其中c 为弧 度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证明面积的误差s ?满足 .s a b c s a b c ????≤++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令
