当前位置：文档之家› 数值分析作业答案(第5章)part2

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

.证明：

(1).如果A 是对称正定矩阵，则1-A 也是对称正定矩阵

(2).如果A 是对称正定矩阵，则A 可以唯一地写成L L A T =，其中L 是具有正对角元的下三角矩阵。

证明：

(1).因A 是对称正定矩阵，故其特征值i λ皆大于0，因此1-A 的特征值1

-i λ也皆大于0。因此1

-i λ也皆大于0，故A 是可逆的。又

111)()(---==A A A T T

则1-A 也是对称正定矩阵。

(2).由A 是对称正定，故它的所有顺序主子阵均不为零，从而有唯一的杜利特尔分解

U L A ~

=。又

022211111

11111DU u u u u u u u u u U n n nn =??

??????

??????????=?????????

其中D 为对角矩阵，0U 为上三角矩阵，于是

~DU L U L A ==

由A 的对称性，得

T T

L D U A A ==

由分解的唯一性得

L U T =

从而

L D L A =

由A 的对称正定性，如果设),,2,1(n i D i =表示A 的各阶顺序主子式，则有

011>=D d ，01

-i i

i D D d ，n i ,,3,2 =

故

121D

D d d d d d d d d d D n n n =??????

??????

??????????

?????=????????????=

因此

T T

LL D L D L L D D L A ===)(21~

121~

，

其中2

D L L =为对角元素为正的下三角矩阵。

.用列主元消去法解线性方程组

???

??=++-=-+-=+-6

1531815331232

1321321x x x x x x x x x 并求出系数矩阵A 的行列式(即A det )的值。

解

??????----?→?-=????

??????----?→???

???

?????----??→?-

=-=?113/110053/7101513

186

76/3118/176/7053/7101513

186111153312151318)(323

21312

1m b A m m r r

所以解为33=x ，22=x ，11=x ，66det -=A 。

.用追赶法解三对角方程组b Ax =，其中

????????????????--------=2100012100012100012100012A ，???????

?????????=00001b 。

解设A 有分解

???????

?????????????????

???----=????????????????--------11

111111211211

211211243

154321

ββββααααα，由公式

???

??===+===-4

,3,2,,5,4,3,2,,

,111111i c i b c b i i i

i i i i βααβαβαα 其中)5,,2,1( =i b i ，)4,,2,1( =i c i 分别是系数矩阵的主对角元素及其下边和上边的次对角线元素。

具体计算，可得

21=α，232=α，343=α，454=α，5

5=α，

211-=β，322-=β，433-=β，5

44-=β。

由

??????

????????=?????????????????????????????????????

?----00001561

451341231254321y y y y y ，

得211=

y ，312=y ，413=y ，514=y ，6

5=y ；再由

???

???????

?????????????=??

????????????????????????

?????????

?--

-61514131211541

2121154321x x x x x ，得615=

x ，314=x ，213=x ，322=x ，6

1=x 。

.下述矩阵能否分解为LU (其中L 为单位下三角矩阵，U 为上三角矩阵)若能分解，那么分解是否唯一

??????????=764142321A ，??????????=133122111B ，??

?????=461561552621C 。

解 A 中02=?，故不能分解。但由于010det ≠-=A ，所以若交换A 的第1行与第3行，则可以分解且分解是唯一的。

在B 中，032=?=?，故不能分解。但B 可以分解为

?????-??????????=33320010011113121u l B ，

其中32l ，33u 为任意常数，且U 奇异，故分解不唯一。

对于C ，)3,2,1(0=≠?i i ，故C 可以分解且分解唯一。

???????????????=131621136121C 。

.求证：(1).∞∞

≤≤x n x x

1；(2).F F A A A n

≤≤21

。证明 (1).由定义知

∞=∞

=≤≤=≤≤∞

==≤=≤=∑∑∑x n x

x x x x x

i n i i n

i n i i i n

i 1

11max max

故∞∞

≤≤x n x x

1。

(2).由范数定义，有

∑∑∑∑∑=======+++==+++≤=n

i F

j ij n

i in

i i n

i i T

T n T T T A a a a a A A tr A A A A A A A A A 1212

21max 2

2)()

()()()( λλλλ

又

221max 2

21)]()()([1)(F T n T T T A n

A A A A A A n A A A =+++≥=λλλλ

所以F F A A A n

≤≤21

。

.设n

n R

P ?∈且非奇异，又x 设为n

R 上一向量范数，定义

Px x P =

试证明P x 是n

R 上向量的一种范数。

证明只需证明P x 满足向量范数的三个条件。

(1).因P 非奇异，故对任意0≠x ，有0≠Px ，故0≥=Px x

，当且仅当0=x 时，

有0==Px x

。

(2).对任意R ∈α，有

P P x Px x P x αααα===。

(3).对任意n

R y x ∈,，有

P P P y x Py Px Py Px y x P y x +=+≤+=+=+)(，故P x 是n R 上的向量

范数。

.设A 为对称正定矩阵，定义()2

1,x Ax x A

=，试证明P x 是n R 上向量的一种范数。

证明只需证明A x 满足向量范数的三个条件。

(1).因A 正定对称，故当0=x ， ()0,2

1==x Ax x

；而当0≠x 时，

()0,2

1>=x Ax x

。

(2).对任意R ∈α，有

()A T T A x Ax x x A x x x A x ααααααα====()()(,2

1。

(3).因A 正定，故有分解T

LL A =，因而

))()(()()(x L x L x L x LL x Ax x x

T T T T T

T T A

====

对任意n

R y x ∈,，由2?的三角不等式有

)(T T

T T T T T A

L L y L x L y L x L y x L y

x +=+≤+=+=+，

故A x 是n

R 上的向量范数。

数值分析试题及答案汇总

数值分析试题一、填空题（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 )

数值分析课后题答案

数值分析第二章 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、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。解因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分）这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。（1分）而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ；解设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分）由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性；解 1）Jacoib 迭代格式为

数值分析作业答案

数值分析作业答案插值法 1、当x=1，-1，2时，f(x)=0，-3，4，求f(x)的二次插值多项式。（1）用单项式基底。（2）用Lagrange插值基底。（3）用Newton基底。证明三种方法得到的多项式是相同的。解：（1）用单项式基底设多项式为: ，所以：所以f(x)的二次插值多项式为：（2）用Lagrange插值基底 Lagrange插值多项式为：所以f(x)的二次插值多项式为： (3) 用Newton基底: 均差表如下： xk f(xk) 一阶均差二阶均差 1 0 -1 -3 3/2 2 4 7/ 3 5/6 Newton插值多项式为：所以f(x)的二次插值多项式为：由以上计算可知，三种方法得到的多项式是相同的。 6、在上给出的等距节点函数表，若用二次插值求ex的近似值，要使截断误差不超过10-6，问使用函数表的步长h应取多少？解：以xi-1,xi,xi+1为插值节点多项式的截断误差，则有式中令得插值点个数

是奇数，故实际可采用的函数值表步长 8、，求及。解：由均差的性质可知，均差与导数有如下关系：所以有： 15、证明两点三次Hermite插值余项是并由此求出分段三次Hermite插值的误差限。证明：利用[xk，xk+1]上两点三次Hermite插值条件知有二重零点xk和k+1。设确定函数k(x): 当或xk+1时k(x)取任何有限值均可；当时，，构造关于变量t的函数显然有在[xk,x][x,xk+1]上对g(x)使用Rolle定理，存在及使得在，，上对使用Rolle定理，存在，和使得再依次对和使用Rolle定理，知至少存在使得而，将代入，得到推导过程表明依赖于及x 综合以上过程有：确定误差限：记为f(x)在[a,b]上基于等距节点的分段三次Hermite插值函数。在区间[xk，xk+1]上有而最值进而得误差估计： 16、求一个次数不高于4次的多项式，使它满足，，。

数值分析试题答案

数值分析试题答案 1、构造拉格朗日插值多项式(X)p 逼近3 (x)f x =，要求（1）取节点011,1x x =-=作线性插值（2）取节点0121,0,1x x x ===作抛物插值答案：（1）代入方程得 0110 10010 1,1(x)y (x x )x y y y y p x x =-=-=+ -=- （2）代入方程得 1202011220120102101220210.1(x x )(x x )(x x )(x x )(x x )(x x ) (x)y x (x x )(x x )(x x )(x x )(x x )(x x )y y p y y ==------= ++=------ 2、给出数据点：01234 39 61215 i i x y =?? =? 用1234,,,x x x x 构造三次牛顿插值多项式3 () N x ，并计算 1.5x =的近似值3(1.5)N 。 33333133.15()93(1) 4.5(1)(2)2(1)(2)(3)(1.5) 5.6250, ()36 4.5(1)3(1)(2)(1.5)7.5000, 1.54 (1.5)(1.5)((1.5)(1.5)) 1.17194 N x x x x x x x N N x x x x x x x N R f N N N =+-+------==+--+--=-=-≈ -=四（分） 3、已知分别用拉格朗日插值法和牛顿插值法求)(x f 的三次插值多项式)(3x P ，并求)2(f 的近似值（保留四位小数）。答案： )53)(43)(13() 5)(4)(1(6 )51)(41)(31()5)(4)(3(2 )(3------+------=x x x x x x x L

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

第一章绪论 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 +=+=假定只用三位数计算,问结果是否可靠?

数值分析第一次作业及参考答案

数值计算方法第一次作业及参考答案 1. 已测得函数()y f x =的三对数据：（0，1），（－1，5），（2，－1），（1）用Lagrange 插值求二次插值多项式。（2）构造差商表。（3）用Newton 插值求二次插值多项式。解：(1)Lagrange 插值基函数为 0(1)(2)1 ()(1)(2)(01)(02)2 x x l x x x +-= =-+-+- 同理 1211 ()(2),()(1)36 l x x x l x x x = -=+ 故 2 20 2151 ()()(1)(2)(2)(1) 23631 i i i p x y l x x x x x x x x x =-==-+-+-++=-+∑ （2）令0120,1,2x x x ==-=，则一阶差商、二阶差商为 011215 5(1) [,]4, [,]20(1) 12 f x x f x x ---= =-= =----- 0124(2) [,,]102 f x x x ---= =- 实际演算中可列一张差商表：（3）用对角线上的数据写出插值多项式 2 2()1(4)(0)1*(0)(1)31P x x x x x x =+--+-+=-+ 2. 在44x -≤≤上给出()x f x e =的等距节点函数表，若用二次插值求x e 的近似值，要使截断误差不超过6 10-，问使用函数表的步长h 应取多少解： ()40000(), (),[4,4],,,, 1.x k x f x e f x e e x x h x x h x x th t ==≤∈--+=+≤考察点及

(3) 2000 4 43 4 3 () ()[(()]()[()] 3! (1)(1) (1)(1) 3!3! .(4,4). 6 f R x x x h x x x x h t t t e t h th t h e h e ξ ξ =----+ -+ ≤+??-= ≤∈- 则 4 36 ((1)(1) 100.006. t t t h - -+± << Q在点得 3.求2 () f x x =在［a,b］上的分段线性插值函数() h I x，并估计误差。解： 22 22 11 1 111 22 11 11 1 () () k k k k h k k k k k k k k k k k k k k k k k k x x x x x x I x x x x x x x x x x x x x x x x x x x x x ++ + +++ ++ ++ + --- =+= --- ?-? -=+- - [] 2 11 22 11 ()()()[()] 11 ()() 44 h h k k k k k k k k R x f x I x x x x x x x x x x x x x h ++ ++ =-=-+- =--≤-= 4.已知单调连续函数() y f x =的如下数据用插值法计算x约为多少时() 1. f x=（小数点后至少保留4位）解：作辅助函数()()1, g x f x =-则问题转化为x为多少时，()0. g x=此时可作新的关于() i g x的函数表。由() f x单调连续知() g x也单调连续，因此可对() g x的数值进行反插。的牛顿型插值多项式为 1()0.110.097345( 2.23)0.451565( 2.23)( 1.10) 0.255894( 2.23)( 1.10)(0.17) x g y y y y y y y - ==-+++++ -++-

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

例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 ρ=，这样，有

数值分析第五版全答案chap1

第一章绪论 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 =,*5 7 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 π= 则何种函数的条件数为 2 3 '4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε=

数值分析作业答案part

6.4.设??? ? ? ??=5010010a b b a A ，0det ≠A ，用a ，b 表示解线性方程组f Ax =的雅可比迭代与高斯—塞德尔迭代收敛的充分必要条件。解雅可比迭代法的迭代矩阵 ? ??? ??? ? ??----=???? ? ??----????? ??=-050100100100000001010101 a b b a a b b a B J ， ?? ? ?? -=-1003||2ab B I J λλλ，10||3)(ab B J = ρ。雅可比迭代法收敛的充分必要条件是3 100 ||

数值分析试卷及其答案1

1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分）解：由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 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分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ

1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： x a x x x a x 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)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使迭代收敛（8分）解：所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分

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

第五章习题解答 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 =进行

数值分析试题及答案.

一、单项选择题（每小题3分，共15分） 1. 和分别作为π的近似数具有（）和（）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ()()2 1 121 1()(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， ()110 l x = B ． () 00l x ＝0， ()111 l x = C ． () 00l x ＝1， ()111 l x = D ． () 00l x ＝1， ()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛，则它具有（）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组12312312 20223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程（）. A ．232 x x -+= B ．232 1.5 3.5 x x -+= C ． 2323 x x -+= D ． 230.5 1.5 x x -=- 单项选择题答案

二、填空题（每小题3分，共15分） 1. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时，科茨系数 ()()() 33301213,88C C C === ，那么() 3 3C = 4. 因为方程 ()420 x f x x =-+=在区间 []1,2上满足，所以 ()0 f x =在区间内有根。 5. 取步长0.1h =，用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公式 . 填空题答案

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

第四版数值分析习题第一章绪论 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 =- 计算,求对数时误差有多大?

数值分析试题及答案

数值分析试题及答案一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（）和（）位有效数字. A．4和3 B．3和2 C．3和4 D．4和4 2. 已知求积公式，则＝（） A． B．C．D． 3. 通过点的拉格朗日插值基函数满足（） A．＝0，B．＝0， C．＝1，D．＝1， 4. 设求方程的根的牛顿法收敛，则它具有（）敛速。 A．超线性B．平方C．线性D．三次 5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程（）. A．B． C．D．单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷人二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案

1. 9和 2. 3. 4. 5. 得分评卷人三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解，，所以分段线性插值函数为 2. 已知线性方程组（1）写出雅可比迭代公式、高斯－塞德尔迭代公式；（2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为雅可比迭代公式为高斯－塞德尔迭代法公式用雅可比迭代公式得用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根（1）请指出为什么初值应取2？（2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案

《数值计算方法》试题集及答案

《计算方法》期中复习试题一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为，拉格朗日插值多项式为。答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )；答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度为( 5 )；

12、为了使计算 32)1(6)1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式1999 2001-改写为 199920012 + 。 13、用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。 14、计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 0.4268 ，用辛卜生公式计算求得的近似值为 0.4309 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。 15、设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿插值多项式为 )1(716)(2-+=x x x x N 。 16、求积公式?∑=≈b a k n k k x f A x x f )(d )(0的代数精度以( 高斯型 )求积公式为最高，具有 ( 12+n )次代数精度。 17、已知f (1)=1,f (3)=5,f (5)=-3,用辛普生求积公式求?5 1 d )(x x f ≈( 12 )。 18、设f (1)=1， f (2)=2，f (3)=0，用三点式求≈')1(f ( 2.5 )。 19、如果用二分法求方程043 =-+x x 在区间]2,1[内的根精确到三位小数，需对分（ 10 ）次。 20、已知?????≤≤+-+-+-≤≤=31)1()1()1(2110)(2 33x c x b x a x x x x S 是三次样条函数，则 a =( 3 )，b =（ 3 ），c =（ 1 ）。 21、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数，则 ∑== n k k x l 0)(( 1 )，∑== n k k j k x l x 0 )(( j x )，当2≥n 时 = ++∑=)()3(20 4x l x x k k n k k ( 32 4++x x )。 22、区间[]b a ,上的三次样条插值函数)(x S 在[]b a ,上具有直到_____2_____阶的连续导数。

数值分析作业答案(第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

数值分析试卷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 的近似值能保证有几位有

文档之家

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

数值分析试题及答案汇总

数值分析课后题答案

数值分析试卷及其答案

数值分析作业答案

数值分析试题答案

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

数值分析第一次作业及参考答案

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

数值分析第五版全答案chap1

数值分析作业答案part

数值分析试卷及其答案1

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

数值分析试题及答案.

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

最新数值分析课程第五版课后习题答案(李庆扬等)1

数值分析试题及答案

《数值计算方法》试题集及答案

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

数值分析试题1