数值分析参考答案第二章
数值分析第二章答案
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)()()3x 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 ==-===-=--==-+-----==------==-+-- 则二次拉格朗日插值多项式为220()()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 =-+=---+-+=+- 5设[]2(),f x Ca b ∈且()()0,f a f b ==求证： 21m ax ()()m ax ().8a x b a x bf x b a f x ≤≤≤≤''≤- 解：令01,x a x b ==，以此为插值节点，则线性插值多项式为10101010()()()x x x x L x f x f x x x x x --=+-- =()()x bx af a f b a b x a --=+--1()()0()0f a f b L x ==∴= 又 插值余项为1011()()()()()()2R x f x L x f x x x x x ''=-=--011()()()()2f x f x x x x x ''∴=--[]012012102()()1()()21()41()4x x x x x x x x x x b a --⎧⎫≤-+-⎨⎬⎩⎭=-=- 又 ∴21m ax ()()m ax ().8a x b a x bf x b a f x ≤≤≤≤''≤- 16．求一个次数不高于4次的多项式P （x ），使它满足(0)(0)0,(1)(1)0,(2)0P P P P P ''=====解：利用埃米尔特插值可得到次数不高于4的多项式0101010,10,10,1x x y y m m ======11300201001012()()()()(12)()(12)(1)j j j j j j H x y x m x x x x xx x x x x x x αβα===+--=---=+-∑∑210110102()(12)()(32)x x x x x x x x x x x α--=---=-2021()(1)()(1)x x x x x xββ=-=-22323()(32)(1)2H x x x x x x x ∴=-+-=-+设22301()()()()P x H x A x x x x =+--其中，A 为待定常数3222(2)1()2(1)P P x x x Ax x =∴=-++-14A ∴= 从而221()(3)4P x x x =-19．求4()f x x =在[,]a b 上分段埃尔米特插值，并估计误差。
数值分析课程第五版课后习题答案(李庆扬等)
数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章：数值分析导论1. 解答：数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题，其结果可能不是完全精确的，但是能够满足工程或科学应用的要求。
2. 解答：数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础，对许多领域都有着广泛的应用，如物理学、经济学、生物学等。
3. 解答：数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比，数值方法一般更加灵活和高效，可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答：计算误差是指数值计算结果与精确解之间的差异。
在数值计算中，由于计算机的有限精度以及数值计算方法本身的近似性等因素，都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章：数值误差分析1. 解答：绝对误差是指实际值与精确值之间的差异。
例如，对于一个计算出的数值近似解x和精确解x_0，其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度，通常被用作评估数值计算方法的好坏。
2. 解答：相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0，其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度，常用于评估数值计算方法的收敛速度。
3. 解答：舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数，因此在进行数值计算过程中，舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答：误差限度是指对于给定的数值计算问题，所能容忍的误差范围。
数值分析(清华大学出版社)第二章课后答案
1.用Gauss 消去法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤---⎢⎢⎢⎢⎣⎡-551631011411014211264321x x x x 解：第一步：交换第三行和第一行，得到如下矩阵⎥⎥⎥⎥⎦⎤----⎢⎢⎢⎢⎣⎡-56153101111402411621做运算()22121E E E →⎪⎭⎫ ⎝⎛+-，()33161E E E →⎪⎭⎫⎝⎛+-，()()441E E E →+，得到增广矩阵 ⎥⎥⎥⎥⎦⎤------⎢⎢⎢⎢⎣⎡0249525213237414210001 第二步：再做运算()3322E E E →+，()44221E E E →⎪⎭⎫⎝⎛+-，得到如下矩阵 ⎥⎥⎥⎥⎦⎤-----⎢⎢⎢⎢⎣⎡94295292113377400210001第三步：做运算()4433713E E E →⎪⎭⎫⎝⎛+，得到 ⎥⎥⎥⎥⎦⎤------⎢⎢⎢⎢⎣⎡21342951919210377400210001利用回代公式求得.790576.0,361257.0,863874.0,115183.11234=-==-=x x x x2、解 2.51 1.48 4.531.480.93 1.302.68 3.041.48⎡⎤⎢⎥-⎢⎥⎢⎥-⎣⎦123x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦=0.051.030.53⎡⎤⎢⎥⎢⎥⎢⎥-⎣⎦ 做两次换行()()()()↔↔3132;E E E E 得2.683.04 1.42.511.48 4.531.480.931.30⎡⎤-⎢⎥⎢⎥⎢⎥-⎣⎦123x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦=0.530.051.03⎡⎤-⎢⎥⎢⎥⎢⎥⎣⎦ 计算()()()()-+→-+→1221330.93657;0.55224;E E E E E E2.683.04 1.481.3672 5.916100.748810.48269⎡⎤-⎢⎥-⎢⎥⎢⎥--⎣⎦123x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦=0.530.546381.3227⎡⎤-⎢⎥⎢⎥⎢⎥⎣⎦计算()()-+→2330.54770;E E E2.683.04 1.4801.36725.9161003.7229⎡⎤-⎢⎥-⎢⎥⎢⎥-⎣⎦123x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦=0.530.546381.0235⎡⎤-⎢⎥⎢⎥⎢⎥⎣⎦ 换行和消去到此结束，经回代计算得到x =()1.440360, 1.577963,0.27494T--3.用Doolittle 三角分解方法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----551631011411014211264321x x x x解：首先对系数矩阵A 做分解LUA =解出：解b y L=，计算出Ty ⎪⎭⎫ ⎝⎛--=74213,521,1,6解y x U=，计算出()T x 115183.1,863874.0,361257.0,790576.0--=4.设][,ij n n a A R A =∈⨯,011≠a ，b Ax =经过高斯消去法一步后变为)2()2(b x A =，其中=)2(A⎥⎦⎤⎢⎣⎡21110A a a T ，(2)A =()(2),2n ij i j a =为(n-1)⨯(n-1)矩阵.其元素为(2)ija =(1)ij a -(1)(1)11i j a a /(1)11a , ,i j =2,3, n. 证明：（1）若A 对称正定，则2A 是对称矩阵。
数值分析第二章答案
∑
n
i=1
ln x i = 0
θ
∧
= −
n
∑ ∑
n
n
i=1
ln x i n
θ
= =
解之得：
i=1
ln x i
（2）母体 X 的期望
E (x) =
∫
+∞ −∞
xf ( x ) d x =
∫
1 0
θ xθ dx =
θ θ +1
而样本均值为：
1 n X = ∑ xi n i =1 令E ( x) = X 得 θ =
x e 2σ 1 n
d x = 2 x ) =
∫
+ ∞ 0
x 2σ
e
−
x σ
d x = − x e ) = 1 ⋅ nσ n
−
x σ
+ ∞
+
0
∫
+ ∞ 0
e
−
x σ
d x =
E (σ ) = E (
∑
n
i=1
i
1 n
∑
n
E ( x
i=1
i
= σ
所以
σ=
∧
1 n ∑ xi σ n i=1 为 的无偏估计量。
∧
X 1− X
5.。解：其似然函数为：
L (σ ) = ∏
i =1
n
1 ⋅e 2σ
−
xi σ
=
1 ⋅e (2σ ) n 1 σ
n i =1
−
1 σ
∑ xi
i =1
n
ln L (σ ) = − n ln(2σ ) − 得： σ =
∧
数值分析第四版课后习题答案
第一章习题解答1、 在下列各对数中，x 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈−=−=aee x a e r π （2）0011.00143.0143.07/1≈=≈−=−=a ee x a e r （3）0127.000004.00031.01000/≈=≈−=−=aee x a e r π （4）001.00143.03.147/100≈=≈−=−=aee x a e r2、已知四个数：001.0,25.134,0250.0,3.264321====x x x x 。
试估计各近似数的有效位数和误差限，并估计运算3211x x x =μ和1431/x x x =μ的相对误差限。
解：21111121101901.0,1021,3,10263.06.23−−⨯≈=⨯==⨯==x x x x n x r δδδ22214212102.0,1021,3,10250.00250.0−−−⨯≈=⨯==⨯==x x x x n x r δδδ 43332333103724.0,1021,5,1013425.025.134−−⨯≈=⨯==⨯==x x x x n x r δδδ 5.0,1021,1,101.0001.04443424==⨯==⨯==−−x x x x n x r δδδ 由相对误差限公式：i r ini n in ni i ir x x fx x f x x x f x x f u δδδ∂∂=∂∂=∑∑==1111),,(),,()(所以有：232123113211103938.0)(1)(−⨯≈++=x x x x x x x x x r δδδμμδ4971.0)(1)(4133141214311≈++−=x x x x x x x x x x r δδδμμδ 3、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
数值分析计算方法第二章作业
1.当x=1，-1，2时，f(x)=0，-3,4，求f(x)的二次差值多项式 （1）用单项式基底 （2）用拉格朗日插值基底
（1）解：设 f(x)abxcx2 则a+b+c=0 a-b+c=-3 a+2b+4c=4
解得
a7,b3,c5 326
所以 f(x)73x5x2
解：由p(0)=0,p(1)=1,p(2)=1,我们可以得出
P 2 ( x ) ( x ( 1 1 ) ) ( ( x 2 ) 2 ) 0 ( 1 ( x ) 0 ( ) x ( 1 2 2 ) ) 1 ( ( 2 x ) ) ( ( 2 x 1 1 ) ) 1 1 2 x 2 3 2 x
将 p'(0)0,p'(1)1 代入到上式中，得出
a 3 ,b 1
4
4
从而有 P4(x)1 4x43 2x39 4x2
p ( x 0 ) f ( x 0 ) , P '( x 0 ) f '( x 0 ) , P ''( x 0 ) f ''( x 0 ) ,p ( x 1 ) f ( x 1 )
解：设 P ( x ) f( x 0 ) f'( x 0 ) ( x x 0 ) f''2 ( x ! 0 )( x x 0 ) 2 a ( x x 0 ) 3
解：设P（x）= ax3bx2cxd
则 P'(x)3ax22bxc
d 0 代入已知条件，得到： c 1
abcd 1 3a 2b c 2
解得a=1,b=-1,c=1,d=0
所以P（x）= x3 x2 x
数值分析课后习题及答案
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
数值分析答案第二章参数估计习题
f(x)= （） { > − ex λ ) λ 0λ ( x解： λe , x ≥ 0
第二章 参数估计 1.设母体X具有负指数分布，它的分布密度 −λ x 为 λe , x ≥ 0 f(x)= 0, x < 0 其中 λ > 0 。试用矩法求的估计量。 解：x e(λ ) f(x)=
0
1
θ −1
dx =
θ θ +1
X 估计EX
X ∴θ = 1− X
1 e 5.设母体X的密度为 f ( x) = 2σ
−
x
σ
, −∞ < x < ∞
试求 σ 的最大似然估计;并问所得估计量是 否的无偏估计. ∑x x n 解： n 1 −σ 1 n − σ
i
L = ∏ f ( xi ) = ∏
i =1 i =1
ln L = n ln θ + (θ − 1)∑ ln xi
i
0, 其他 n
i =1
( θ >0 )
n i =1
d ln L n ＾= − n = + ∑ ln xi = 0,∴θ θ i dθ ∑ ln xi
i
2矩法估计
EX =
用
X 用估计EX
+∞
−∞
∫ x ⋅ f ( x)dx = ∫ x ⋅θ ⋅ x
2
给定置信概率1−α 即
P ( x − uα
2
σ/ n
，有 uα ，使
2
P{ u ≤ uα } = 1 − α
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
若
则步长
在小区间 上，分段线性插值函数为
各节点间中点处的 与 的值为
当 时，
当 时，
当 时，
当 时，
当 时，
误差
又
令
得 的驻点为 和
18．求 在 上分段线性插值函数 ，并估计误差。
解：
在区间 上，
函数 在小区间 上分段线性插值函数为
误差为
19．求 在 上分段埃尔米特插值，并估计误差。
解：
在 区间上，
（2）若 ，则
证明：
（1）
得证。
+
得证。
14． 求 及 。
解：
若
则
15．证明两点三次埃尔米特插值余项是
解：
若 ，且插值多项式满足条件
插值余项为
由插值条件可知
且
可写成
其中 是关于 的待定函数，
现把 看成 上的一个固定点，作函数
根据余项性质，有
由罗尔定理可知，存在 和 ，使
即 在 上有四个互异零点。
根据罗尔定理， 在 的两个零点间至少有一个零点，
则
若采用二次插值法计算 时，
3．给全 的函数表，步长 若函数表具有5位有效数字，研究用线性插值求 近似值时的总误差界。
解：求解 近似值时，误差可以分为两个部分，一方面，x是近似值，具有5位有效数字，在此后的计算过程中产生一定的误差传播；另一方面，利用插值法求函数 的近似值时，采用的线性插值法插值余项不为0，也会有一定的误差。因此，总误差界的计算应综合以上两方面的因素。
第二章插值法
1．当 时， ,求 的二次插值多项式。
解：
则二次拉格朗日插值多项式为
2．给出 的数值表
X
0.4
0.5
0.6
0.7
0.8
lnx
-0.916291
-0.693147
-0.510826
-0.356675
-0.223144
用线性插值及二次插值计算 的近似值。
解：由表格知，
若采用线性插值法计算 即 ，
令
函数 在区间 上的分段埃尔米特插值函数为
误差为
又
20．给定数据表如下：
Xj
0.25
0.30
0.39
0.45
0.53
Yj
0.5000
0.5477
0.6245
0.6708
0.7280
试求三次样条插值，并满足条件：
解：
由此得矩阵形式的方程组为
21 M0
2 M1
2 M2
2 M3
12 M4
求解此方程组得
三次样条表达式为
当 时，
令
取
令
则
当 时，线性插值多项式为
插值余项为
又 在建立函数表时，表中数据具有5位有效数字，且 ，故计算中有误差传播过程。
总误差界为
4．设为互异节点，求证：
（1）
（2）
证明
（1）令
若插值节点为 ，则函数 的 次插值多项式为 。
插值余项为
又
由上题结论可知
得证。
5设 且 求证：
解：令 ，以此为插值节点，则线性插值多项式为
将 代入得
由此得矩阵开工的方程组为
求解此方程组，得
又 三次样条表达式为
将 代入得
21．若 是三次样条函数，证明：
若 ，式中 为插值节点，且 ，则
证明：
从而有
=
插值余项为
6．在 上给出 的等距节点函数表，若用二次插值求 的近似值，要使截断误差不超过 ，问使用函数表的步长h应取多少？
解：若插值节点为 和 ，则分段二次插值多项式的插值余项为
设步长为h，即
若截断误差不超过 ，则
7．若 ，
解：根据向前差分算子和中心差分算子的定义进行求解。
8．如果 是m次多项式，记 ，证明 的k阶差分 是 次多项式，并且 （ 为正整数）。
解：函数 的 展式为
其中
又 是次数为 的多项式
为 阶多项式
为 阶多项式
依此过程递推，得 是 次多项式
是常数
当 为正整数时，
9．证明
证明
得证
10．证明
证明：由上题结论可知
得证。
11．证明
证明
得证。
12．若 有 个不同实根 ，
证明：
证明： 有个不同实根
且
令
则
而
令
则
又
得证。
13．பைடு நூலகம்明 阶均差有下列性质：
（1）若 ，则
故 在 内至少有三个互异零点，
依此类推， 在 内至少有一个零点。
记为 使
又
其中 依赖于
分段三次埃尔米特插值时，若节点为 ，设步长为 ，即
在小区间 上
16．求一个次数不高于4次的多项式P（x），使它满足
解：利用埃米尔特插值可得到次数不高于4的多项式
设
其中，A为待定常数
从而
17．设 ，在 上取 ，按等距节点求分段线性插值函数 ，计算各节点间中点处的 与 值，并估计误差。