计算方法的课后答案
计算方法-刘师少版第一章课后习题完整答案
分， 试给出此递推公式误差的传播规律， 计算 I 10 时误差被放大了多少倍？这个算法是数值稳定的 吗？ 解： I =
∫x
0 1 0
1
n
e x −1 dx ， n = 0,1,2,L,10 ，由分部积分法有
1 0
n −1 x −1 I n = ∫ x n e x −1 dx = x n e x −1 1 e dx 0 − n∫ x
er ( x n ) =
e( x n ) nx n −1 ( x − x * ) x − x* = = n = n ⋅ er ( x) = αn% x xn xn
x n 的相对误差为 an%
1.10 设 x＞0，x 的相对误差为 δ ，求 ln x 的误差。 解： e(ln x) ≈
1 ( x − x * ) = er ( x) = δ x
N +1
N
1 dx = arctan( N + 1) − arctan N 1+ x2 1 = arctan 1 + N ( N + 1) 1 2 gt ，假定 g 是准确的，而对 t 的测量有±0.1s 的误差，证明当 t 增加时，s 的绝对误差 2
1.12 设 s =
增加，而相对误差减少。 解：由题意知， e( s ) = s − s = gt (t − t ) = gt ⋅ e(t ) = 0.1gt
5
计算方法
于是
* * * * e( I 10 ) = −10e( I 9 ) = 10 ⋅ 9e( I 8 ) = L = 10!e( I 0 )
计算 I 10 时的误差被扩大了 10 倍，显然算法是数值不稳定的 1.14 设 f ( x) = 8 x − 0.4 x + 4 x − 9 x + 1 ，用秦九韶算法求 f (3)
计算方法习题第一、二章答案
第一章 误差1 问，，722分别作为π的近似值各具有几位有效数字？分析 利用有效数字的概念可直接得出。
解 π= 592 65… 记x 1=，x 2=，x 3=722.由π- x 1= 59…= 40…知3411110||1022x π--⨯<-≤⨯ 因而x 1具有4位有效数字。
由π- x 2= 59…= 59…知2231021||1021--⨯≤-<⨯x π因而x 2具有3位有效数字。
由π-722= 59 … 85…= 26…知231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。
2 已知近似数x*有两位有效数字，试求其相对误差限。
分析 本题显然应利用有效数字与相对误差的关系。
解 利用有效数字与相对误差的关系。
这里n=2,a 1是1到9之间的数字。
%5101211021|*||*||)(|1211*=⨯⨯≤⨯≤-=+-+-n ra x x x x ε3 已知近似数的相对误差限为%，问x*至少有几位有效数字？分析 本题利用有效数字与相对误差的关系。
解 a 1是1到9间的数字。
1112*10)1(2110)19(21102110003%3.0)(--⨯+≤⨯+⨯=⨯<=a x r ε 设x*具有n 位有效数字，令-n+1=-1，则n=2，从而x*至少具有2位有效数字。
4 计算，问要取几位有效数字才能保证相对误差限不大于%。
分析 本题应利用有效数字与相对误差的关系。
解 设取n 位有效数字，由=…，故a 1=9。
411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式411101021-+-≤⨯n a 知取n=4即可满足要求。
5 计算76017591-，视已知数为精确值，用4位浮点数计算。
解 =-76017591 8×10-2-0.131 6×10-2=×10-5结果只有一位有效数字，有效数字大量损失，造成相对误差的扩大，若通分后再计算：56101734.0105768.01760759176017591-⨯=⨯=⨯=- 就得到4位有效数字的结果。
数值计算方法课后习题答案
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
计算方法 课后习题答案
得到方程组
3。举例说明一个非奇异矩阵不一定存在LU分解。
例如：设
与题设相矛盾，所以一个非奇异矩阵不一定存在LU分解。
4。下列矩阵能否分解为LU（其中L为单位下三角矩阵，U为上三角矩阵）？若能分解，那么分解是否唯一？
解：
设 B可以进行LU分解，则B=
计算得
其中。 。
解：(1)由题意，可设 ，由Lagrange插值余项公式得
(2)由(1)式可知，
15．给定数据表：
1
0
2
3
构造出函数 的差商表，并写出它的三次 插值多项式.
解：利用Newton插值公式：
先作出差商表
一阶差商
二阶差商
三阶差商
0
1
3
1
3/2
13/4
1/2
2
0
3
1/6
1/3
3
2
5/3
-2/3
-5/3
证明：据题4可知，
令 ，则有 。注意到
（证明见王能超数值简明教程145页题6）
令 即有 。
9.已知 ，求差商 和 。
解：根据差商与微商的关系，有
10.已知 互异，求 。其中 。（此题有误。）（见王能超《教程》P149－题2）
解：因为 ，则
由差商性质 可知，
11.设首项系数为1的n次式 有n个互异的零点 ，证明
解：1）用梯形公式有：
事实上，
2）Simpson公式
事实上，
3）由Cotes公式有：
事实上，
2．证明Simpson公式 具有三次代数精度。
证明：
而当 时
左侧：
右侧：
计算方法习题答案
计算方法习题答案计算方法习题答案计算方法是一门重要的学科，它在我们的日常生活中扮演着重要的角色。
无论是在工作中还是在学习中，我们都需要运用计算方法来解决各种问题。
然而，有时我们可能会遇到一些困难，无法找到正确的答案。
在这篇文章中，我将为你提供一些计算方法习题的答案，并解释一些基本的计算方法概念。
第一题：计算两个数的和答案：要计算两个数的和，只需将这两个数相加即可。
例如，如果给定的两个数是3和5，那么它们的和就是3+5=8。
第二题：计算一个数的平方答案：要计算一个数的平方，只需将这个数乘以自己。
例如，如果给定的数是4，那么它的平方就是4*4=16。
第三题：计算一个数的百分比答案：要计算一个数的百分比，需要将这个数乘以百分比的表示形式，并将结果除以100。
例如，如果要计算50的10%，则计算方法是50*(10/100)=5。
第四题：计算一个数的平均值答案：要计算一组数的平均值，需要将这些数相加，然后将结果除以数的个数。
例如，如果给定的一组数是3、4和5，那么它们的平均值就是(3+4+5)/3=4。
第五题：计算一个数的阶乘答案：要计算一个数的阶乘，需要将这个数与比它小1的数相乘，并继续乘以比前一次乘积小1的数，直到乘到1为止。
例如，如果要计算5的阶乘，则计算方法是5*4*3*2*1=120。
以上是一些常见的计算方法习题的答案。
通过这些例子，我们可以看到计算方法在解决实际问题中的应用。
无论是在日常生活还是在工作中，掌握计算方法是非常重要的。
除了以上习题的答案，还有一些更复杂的计算方法可以应用于更高级的问题。
例如，线性规划是一种常用的优化方法，可以用于解决最大化或最小化目标函数的问题。
数值积分是一种用于计算曲线下面积的方法，可以应用于物理学、经济学等领域。
这些方法需要更深入的学习和理解，但它们在解决实际问题中起到了重要的作用。
总结起来，计算方法是一门重要的学科，它在我们的日常生活中扮演着重要的角色。
通过掌握基本的计算方法概念和解题技巧，我们可以更好地解决各种问题。
计算方法课后习题集规范标准答案
习 题 一3.已知函数y =4, 6.25,9x x x ===处的函数值，试通过一个二次插值函解：0120124, 6.25,9;2, 2.5,3y x x x y y y =======由题意 （1） 采用Lagrange插值多项式220()()j j j y L x l x y ==≈=∑27020112012010*********()|()()()()()()()()()()()()(7 6.25)(79)(74)(79)(74)(7 6.25)2 2.532.255 2.25 2.75 2.7552.6484848x y L x x x x x x x x x x x x x y y y x x x x x x x x x x x x ==≈------=++------------=⨯+⨯+⨯⨯-⨯⨯= 其误差为(3)25(3)25(3)2[4,9]2()(7)(74)(7 6.25)(79)3!3()83max |()|40.0117281|(7)|(4.5)(0.01172)0.008796f R f x x f x R ξ--=---==<∴<=又则（2）采用Newton插值多项式2()y N x =≈ 根据题意作差商表：224(7)2(74)()(74)(7 6.25) 2.64848489495N =+⨯-+-⨯-⨯-≈4. 设()()0,1,...,k f x x k n ==，试列出()f x 关于互异节点()0,1,...,i x i n =的Lagrange 插值多项式。
注意到：若1n +个节点()0,1,...,i x i n =互异，则对任意次数n ≤的多项式()f x ，它关于节点()0,1,...,i x i n =满足条件(),0,1,...,i i P x y i n ==的插值多项式()P x 就是它本身。
可见，当k n ≤时幂函数()(0,1,...,)kf x x k n ==关于1n +个节点()0,1,...,i x i n =的插值多项式就是它本身，故依Lagrange 公式有()00(),0,1,...,nn n k kk i j j j j j i j ii jx x x l x x x k n x x ===≠-=≡=-∑∑∏特别地，当0k =时，有()0001nn n ij j j i j ii jx x l x x x ===≠-=≡-∑∑∏而当1k =时有()000nnn ij j j j j i j ii jx x x l x x x x x ===≠⎛⎫- ⎪=≡ ⎪- ⎪⎝⎭∑∑∏ 5.依据下列函数表分别建立次数不超过3的Lagrange 插值多项式和Newton 插值多项式，并验证插值多项式的唯一性。
计算方法引论课后答案
计算方法引论课后答案第一章误差1.什么是模型误差，什么是方法误差？例如，将地球近似看为一个标准球体，利用公式 $A=4\pi r$ 计算其表面积，这个近似看为球体的过程产生的误差即为模型误差。
在计算过程中，要用到 $\pi$，我们利用无穷乘积公式计算 $\pi$ 的值：pi=2\cdot\frac{2}{1}\cdot\frac{2}{3}\cdot\frac{4}{3}\cdot\f rac{4}{5}\cdot\frac{6}{5}\cdot\frac{6}{7}\cdot\frac{8}{7}\cdot\ frac{8}{9}\cdot\cdots我们取前9项的乘积作为 $\pi$ 的近似值，得$\pi\approx3.xxxxxxxx5$。
这个去掉 $\pi$ 的无穷乘积公式中第9项后的部分产生的误差就是方法误差，也称为截断误差。
2.按照四舍五入的原则，将下列各数舍成五位有效数字：816.956，76.000，.322，501.235，.182，130.015，236.23.解：816.96，76.000，.501.24，.130.02，236.23.3.下列各数是按照四舍五入原则得到的近似数，它们各有几位有效数字？81.897，0.008，136.320，050.180.解：五位，三位，六位，四位。
4.若 $1/4$ 用 0.25 表示，问有多少位有效数字？解：两位。
5.若 $a=1.1062$，$b=0.947$，是经过舍入后得到的近似值，问：$a+b$，$a\times b$ 各有几位有效数字？已知 $da<\frac{1}{2}\cdot10^{-4}$，$db<\frac{1}{2}\cdot10^{-3}$，又 $a+b=0.\times10$。
begin{aligned}d(a+b)&=da+db\leq da+db=\frac{1}{2}\cdot10^{-4}+\frac{1}{2}\cdot10^{-3}=0.55\times10^{-3}<\frac{1}{2}\cdot10^{-2}end{aligned}所以 $a+b$ 有三位有效数字；因为 $a\timesb=0.xxxxxxxx\times10$。
计算方法答案王能超
计算方法答案王能超【篇一：计算方法习题集及实验指导书】s=txt>计算机科学与技术系檀明2008-02-10课程性质及目的要求（一）课程性质自计算机问世以来，科学计算一直是计算机应用的一个重要领域，数值计算方法是解决各种复杂的科学计算问题的理论与技术的基础。
《计算方法》课程讨论用于科学计算中的一些最基本、最常用的算法，不但具有数学的抽象性与严密的科学性的特点，而且具有应用的高度技术性的特点。
它对于培养从事计算机应用的科技人才有着重要的作用，是计算机应用专业（本科段）的一门重要的技术基础课程。
（二）目的要求通过本课程的学习和上机实验，了解用计算机解决科学计算问题的方法特点，掌握计算方法中的一些基本概念、基本公式和相应的算法流程，提高根据算法描述设计高级语言程序并进行验证的技能。
在学习过程中，应注重理解和应用，在搞清基本原理和基本概念的基础上，通过习题、编程和上机等环节，巩固和加深已学的内容，掌握重要的算法及其应用。
注重理论与算法的学习和应用相结合，强调编程及上机计算的技能培养，是本课程不同于一般数学课程的重要特点。
（三）学习方法指导1．循序渐进逐章学习本课程从第二章开始，每章都讨论一个大类的算法。
虽然各算法是相对独立的，但是也存在相互联系与前后继承的关系。
前面的概念和算法学好了，后面的内容也就容易学，越学越感到容易。
前面的内容没有学好，后面就会感到难学，甚至会出现越来越感到困难、失去学习信心的情况。
2．稳扎稳打融会贯通学习要扎实、要讲求实效。
每一个重要的概念和公式，都会搞清楚，做到融会贯通。
只有这样，才能取得学习的学习效果。
3．多学练勤做习题教材及本习题集中的每一章都附有适量的习题，可以帮助考生巩固和加深理解所学的知识，提高解题能力。
因此，在学习过程中，应当适合习题进行思考，应当尽可能多做习题，遇到某些不会做的题，应三思之后再请老师给予提示。
4．抓住特点前后联系本课程只讲了五大类算法。
每类算法都是针对一类特定的计算问题，都有其自身的特点。
数值计算方法课后习题答案(李庆扬等) (修复的)
，。
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x xx x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++；[解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
计算方法_课后习题答案
L3 x 的最高次项系数是 6，试确定 y1 。
解： l0 (x)
x x1 x0 x1
x x2 x0 x2
x x3 x0 x3
x 0.5 0 0.5
x 1 0 1
x2 02
= x3
7 2
x2
7 2
x 1
l1 ( x)
x x0 x1 x0
(2 2e1 4e0.5 )x2 (4e0.5 e1 3)x 1
2)根据Lagrange余项定理，其误差为
| R2 (x) ||
f
(3) ( 3!
)
21
(
x)
||
1 6
e
x(
x
1)(
x
0.5)
|
1 max | x(x 1)(x 0.5) |, (0,1) 6 0x1
x2 02
x4= 04
x3
7x2 14x 8 8
l1 ( x)
x x0 x1 x0
x x2 x1 x2
x x3 x1 x3
x0 1 0
x2 1 2
x4 1 4
=
x3
6x2 3
8x
l2 (x)
x x0 x2 x0
i j
而当 k 1时有
n
x jl j
j0
x
n
n
j0 i0 i j
x xi x j xi
x
j
计算方法引论课后答案
第一章 误差1. 试举例,说明什么是模型误差,什么是方法误差.解: 例如,把地球近似看为一个标准球体,利用公式24A r π=计算其表面积,这个近似看为球体的过程产生的误差即为模型误差.在计算过程中,要用到π,我们利用无穷乘积公式计算π的值: 其中我们取前9项的乘积作为π的近似值,得这个去掉π的无穷乘积公式中第9项后的部分产生的误差就是方法误差,也成为截断误差.2. 按照四舍五入的原则,将下列各数舍成五位有效数字:816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 2363. 下列各数是按照四舍五入原则得到的近似数,它们各有几位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位 三位 六位 四位4. 若1/4用0.25表示,问有多少位有效数字? 解: 两位5. 若 1.1062,0.947a b ==,是经过舍入后得到的近似值,问:,a b a b +⨯各有几位有效数字?解: 已知4311d 10,d 1022a b --<⨯<⨯, 又0.2053210a b +=⨯,()433211110100.551010222d a b da db da db ----+=+≤+=⨯+⨯=⨯<⨯,所以a b +有三位有效数字;因为0.1047571410a b ⨯=⨯,所以a b ⨯有三位有效数字.6. 设120.9863,0.0062y y ==,是经过舍入后作为12,x x 的近似值.求1211,y y 的计算值与真值的相对误差限及12y y ⋅与真值的相对误差限.解: 已知-4-41112221211d ,d ,d =10,d 1022x y x x y x x x =+=+⨯=⨯, ()44111111110d d 12dr dr 0.50100.9863x xx x x y --⨯⎛⎫==≈=≈⨯ ⎪⎝⎭;()42222222110d d 12dr dr 0.81100.0062x xx x x y --⨯⎛⎫==≈=≈⨯ ⎪⎝⎭;()()()4221212dr dr dr 0.50100.81100.8210x x x x ---⋅=+≈⨯+⨯≈⨯.7. 正方形的边长约为100cm,应该怎样测量,才能使其面积的误差不超过1cm 2. 解: 设正方形面积为S,边长为a,则S=a 2.所以要使:2d d 2d 1s a a a ==≤,则要求211d 0.5102200a a -≤==⨯.所以边长的误差不能超过20.510-⨯cm.8. 用观测恒星的方法求得某地维度为4502'''(读到秒),试问:计算sin ϕ将有多大误差?解: ()()1d sin cos d cos 45022ϕϕϕ*''⎛⎫'''== ⎪⎝⎭.9 . 真空中自由落体运动距离s 与时间的关系由公式212s gt =确定,g 是重力加速度.现在假设g 是准确的,而对t 的测量有0.1s ±的误差,证明t 增加时,距离的绝对误差增加而相对误差却减小. 证明: 因为:221d d d d d d d ;2.122s gt t gt t t s gt gt t s s t gt ⎛⎫=====⎪⎝⎭ d s 与t成正比,d s s与t 成反比,所以当d t 固定的时候, t增加时,距离的绝对误差增加而相对误差却减小.10. 设0x >,x 的相对误差为δ,求ln x 的绝对误差. 解: 已知d x x δ=,所以ln x 的绝对误差()d d ln x x xδ==. 11. 设x 的相对误差为%α,求nx 的相对误差.解: 1d d d %n n n nx nx x n xn x x xα-===. 12. 计算球的体积,为了使相对误差限为1%,问度量半径R 时允许的相对误差限如何? 解: 已知343V R π=,设()d dr R R a R ==,则要使得 ()()3d dr dln d ln 3d ln 3d ln 3dr 31%V V V R R R R a V ========,则11%3a =⋅. 第二章 插值法与数值微分1.设y =在100,121,144x =三处的值是很容易求得的,试以这三个点建立y =的二次插值多项式,并用,且给出误差估计.用其中的任意两点,构造线性插值函数,用得到的三个线性插值函数,,并分析其结果不同的原因.解: 已知012012100,121,144;10,11,12x x x y y y ======,建立二次Lagrange 插值函数可得:()211510.7228L ≈=.误差()()()()()()2012012,,,,3!f R x x x x x x x x x x ξξξ'''=---∈,所以利用前两个节点建立线性插值函数可得:()111510.7143L ≈=.利用后两个节点建立线性插值可得:()111510.7391L ≈=.利用前后两个节点建立线性插值可得:()111510.6818L ≈=.与,二次插值比线性插值效果好,利用前两个节点的线性插值比其他两个线性插值效果好.此说明,二次插值比线性插值效果好,内插比外插效果好.2. 利用(2.9)式证明 证明: 由(2.9)式当01x x x <<时,()()01max x x x f f x ξ≤≤''''≤,()()()01201101max 4x x x x x x x x x ≤≤--≤- 所以3. 若()0,1,...,j x n 为互异节点,且有 证明 证明: 由于且()0nk j j j x l x =∑和kx都为k 次多项式,而且在k+1个不同的节点处的函数值都相同0,1,...,k n =, 所以马上有()0,0,1,...,nk kj j j x l x xk n =≡=∑.4. 设给出sin x 在[],ππ-上的数值表,用二次插值进行计算,若希望截断误差小于510-,问函数表的步长最大能取多少?解: 记插值函数为p(x),则所以()()()()11cos max sin 3!i i i x x p x x x x x x ππξ-+-≤≤--=---()cos 1ξ-≤;令()()()()11i i i g x x x x x x x -+=---,设1i x x th -=+,得又()()()[]12,0,2t t t t t ϕ=--∈的最大值为10.3849ϕ⎛= ⎝⎭,所以有 所以0.0538h ≤.5. 用拉格朗日插值和牛顿插值找经过点()()()()3,1,0,2,3,2,6,10---的三次插值公式. 解: Lagrange 插值函数:牛顿插值: 首先计算差商也可以利用等距节点构造,首先计算差分 可得前插公式 和后插公式6. 确定一次数不高于4的多项式()x ϕ,使()()()()()00,00,111,21ϕϕϕϕϕ''=====. 解: 利用重节点计算差商则可构造Hermite 插值函数满足题设条件:7. 寻找过1n +个点01,,...,n x x x 的21n +次多项式()21n H x +,满足条件: 解: 和Lagrange 插值函数的构造类似,可将插值函数写成其中,基函数满足条件 (1)()()(),,,21n i n i h x h x P n ∈+;(2)()()()(),,,,,0;,0n i n i n ij ij n i j j ijj h x h x h x h x δδ''====则可由已知条件,可得()()()()2,,,12n i n i i i n i h x l x x x l x '⎡⎤=--⎣⎦;()()()2,,n i i n i h x x x l x '=-.所以可得8. 过0,1两点构造一个三次Hermite 插值多项式,满足条件: 解: 计算重节点的差商马上可得9. 过给定数组(1) 作一分段线性插值函数.(2) 取第二类边界条件,作三次样条插值多项式.(3) 用两种插值函数分别计算75.5,78.3x x ==的函数值. 解: (1)做分段线性插值函数可得:其中, ()[][]076 75,76;0 75,76.x x l x x ⎧-∈⎪=⎨∉⎪⎩()[][][]175 75,7677 76,77;0 75,77.x x l x x x x ⎧-∈⎪=-∈⎨⎪∉⎩ (2)把已知节点值带入M 关系式可得: 由边界条件可得050M M ==,所以上面方程组变为可求解方程组解得12340.0058,0.0067,0.0036,0.0071M M M M ====.所以可得在每个区间上的三次样条函数的表达式: (3)当75.5x =时,()()()50175.5 2.76875.5 2.83375.5 2.8005I l l =+=;()()()()()30.00580.005875.575.576 2.7687675.5 2.83375.575 2.79966s ⎛⎫=-+-+--= ⎪⎝⎭当78.3x =时,()()()53475.5 2.97978.3 3.06278.3 3.0039I ll =+=;10. 若给出sin ,cos ,tan x x x 的函数表:用表上的数据和任一插值公式求: (1) 用tan x 表格直接计算tan1.5695.(2) 用sin1.5695和cos1.5695来计算tan1.5695.并讨论这两个结果中误差变化的原因. 解: 利用Lagrange 插值直接用tan 表计算得tan1.5695819.0342874999274≈;利用Lagrange 插值计算sin 得sin1.56950.99999917500000≈;利用Lagrange 插值计算cos 得cos1.56950.00129630000000≈;最后利用sin/cos 计算tan 得tan1.5695771.4257309264500≈.出现小除数,误差被放大.11. 求三次样条函数()s x ,已知和边界条件解: 把表中数据带入M 关系式可得由边界条件还可得到两个方程: 联立两个方程组可解得:带入M 表达式便可得所求三次样条函数.12. 称n 阶方阵()ij A a =具有严格对角优势,若 (1) 试证明:具有严格对角优势的方阵必可逆. (2) 证明:方程组(2.62)解存在唯一.证明: (1)设矩阵A 按行严格对角占优,如果A 奇异,则存在非零向量x 使得Ax=0,写成分量形式为令指标0i 使得00i x x∞=≠,则因此0000010n i i i i j j j i x a a =≠⎛⎫⎪-≤ ⎪ ⎪⎝⎭∑ 即000010ni i i j j j i a a =≠-≤∑上式与矩阵按行严格对角占优矛盾,因此矩阵非奇异. (2)方程组(2.62)由于该方程组系数矩阵为严格对角占优的方阵,所以由克拉默法则可知方程组存在唯一解.。
数值计算方法课后习题答案(李庆扬等)
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
计算方法 课后习题答案
计算方法课后习题答案计算方法课后习题答案计算方法是一门重要的学科，它涉及到数值计算、算法设计和数据处理等方面的内容。
在学习计算方法的过程中，课后习题是不可或缺的一部分。
通过解答习题，我们可以巩固所学的知识，提高自己的计算能力。
下面是一些计算方法课后习题的答案，希望对大家的学习有所帮助。
1. 矩阵的转置矩阵的转置是将矩阵的行和列互换得到的新矩阵。
对于一个m×n的矩阵A，它的转置记作A^T。
转置后的矩阵A^T的行数和列数分别为原矩阵A的列数和行数。
例如，对于一个3×2的矩阵A，它的转置A^T是一个2×3的矩阵。
2. 矩阵的加法和减法矩阵的加法和减法是对应位置上的元素进行相加或相减得到的新矩阵。
对于两个相同大小的矩阵A和B，它们的和记作A+B，差记作A-B。
加法和减法的运算规则是相同位置上的元素进行相应的运算。
3. 矩阵的乘法矩阵的乘法是指将两个矩阵相乘得到一个新矩阵的运算。
对于两个矩阵A和B，它们的乘积记作AB。
矩阵乘法的运算规则是矩阵A的行与矩阵B的列进行相乘，并将结果相加得到新矩阵的对应位置上的元素。
4. 矩阵的逆矩阵的逆是指对于一个可逆矩阵A，存在一个矩阵B，使得AB=BA=I，其中I是单位矩阵。
如果一个矩阵A存在逆矩阵，则称其为可逆矩阵或非奇异矩阵。
求解矩阵的逆可以使用伴随矩阵和行列式的方法。
5. 线性方程组的求解线性方程组是指由一组线性方程组成的方程组。
求解线性方程组的方法有很多，包括高斯消元法、LU分解法、迭代法等。
其中，高斯消元法是一种常用的求解线性方程组的方法，它通过消元和回代的过程，将线性方程组转化为上三角形矩阵或对角矩阵，从而求解出方程组的解。
6. 数值积分的方法数值积分是指通过数值计算的方法来求解定积分的近似值。
常用的数值积分方法包括梯形法则、辛普森法则和龙贝格法则等。
这些方法都是基于将定积分转化为离散求和的形式，通过计算离散点上的函数值来估计定积分的近似值。
计算方法课后习题答案
计算方法课后习题答案在计算方法课程中，学生通常会接触到各种数学问题的求解方法，包括但不限于数值分析、线性代数、微分方程等。
以下是一些课后习题的解答示例：习题一：求解线性方程组设线性方程组为：\[ \begin{align*}a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n &= b_1, \\a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n &= b_2, \\\vdots \quad \quad & \ \vdots \\a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n &= b_m,\end{align*} \]解答：使用高斯消元法或矩阵分解法求解上述方程组。
首先将系数矩阵转换为行简化阶梯形式，然后回代求解未知数 \( x_1, x_2,\ldots, x_n \)。
习题二：数值积分给定函数 \( f(x) \)，需要在区间 \( [a, b] \) 上进行数值积分。
解答：可以使用梯形法、辛普森法等数值积分方法。
例如，使用梯形法的公式为：\[ \int_a^b f(x)dx \approx \frac{h}{2} \left( f(a) + 2f(a+h) + 2f(a+2h) + \cdots + 2f(b-h) + f(b) \right), \]其中 \( h = \frac{b-a}{n} \) 是区间的等分宽度，\( n \) 是等分数。
习题三：常微分方程的数值解给定一个常微分方程 \( y' = f(x, y) \)，初始条件为 \( y(x_0) = y_0 \)。
解答：使用欧拉法或龙格-库塔法求解。
以欧拉法为例，其迭代公式为：\[ y_{n+1} = y_n + h f(x_n, y_n), \]其中 \( h \) 是步长，\( x_{n+1} = x_n + h \)。
最优化计算方法课后习题集答案解析
解：取 ， 时，DFP法的第一步与最速下降法相同
， ，
，
以下作第二次迭代
，
其中，
，
所以
令 ， 利用 ，求得
所以 ，
以下作第三次迭代
，
，
所以
令 ， 利用 ，求得
所以 ， 因为 ，于是停止
即为最优解。
习题四
包括题目： P95页 3；4；8；9(1)；12选做；13选做
3题解如下
3.考虑问题 ，其中
X1,x2,x3≥0 （3）
求出点（1,1,0）处的一个下降可行方向.
解：首先检查在点（1,1,0）处哪些约束为有效约束。检查易知（1），X3≥0为有效约束。设所求可行方向d=(d1,d2,d3)T。根据可行方向d的定义，应存在a>0，使对∀t∈（0，a）能有
X+td=(1+td1,1+td2,0+td3)T
（1）
s.t.
（2）
s.t.
（1）解：非线性规划的K-T条件如下：
（1）
（2）
（3）
再加上约束条件 （4）
为求出满足（1）~（4）式的解，分情况考虑：
①若（4）式等号不成立，即 ，那么由（2）式得 ，将 代入（1）式解得 ， ，所得值不满足 的条件，故舍去。
②若（4）式等号成立，由（1）式可以解得 ， ，代入（4）式有：
JBi
1
2
3
4
5
6
7
8
9
di0
1
1
0
-5/6
-1/6
1
10/6
4
0
0
38/6
2
0
1
-9/6
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《计算方法》习题答案第一章 数值计算中的误差1．什么是计算方法？（狭义解释）答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。
2．一个实际问题利用计算机解决所采取的五个步骤是什么？答：一个实际问题当利用计算机来解决时，应采取以下五个步骤： 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值，并编程获得解。
解：400)(2345-+⋅+-⋅+=x x x x x x P ，从而 1 0 -1 0 1 -4 -3 -3 9 -24 72 -2191-38-2473-223所以，多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。
5．叙述误差的种类及来源。
答：误差的种类及来源有如下四个方面：（1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。
（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。
（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。
（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。
这样引起的误差称为舍入误差。
6．掌握绝对误差（限）和相对误差（限）的定义公式。
答：设*x 是某个量的精确值，x 是其近似值，则称差x x e -=*为近似值x 的绝对误差（简称误差）。
若存在一个正数ε使ε≤-=x x e *，称这个数ε为近似值x 的绝对误差限（简称误差限或精度）。
把绝对误差e 与精确值*x 之比***x x x x e e r -==称为近似值x 的相对误差，称*x εη=为近似值x 的相对误差限η≤r e ，由于真值*x 是未知的，所以常常用xe x x x e r =-=*来表示相对误差，于是相对误差可以从绝对误差求出。
7．近似值的规格化表示形式如何？答：一般地，对于一个精确值*x ，其近似值x 的规格化形式为m p x x x x 10.021⨯±= ，其中{}),2,1(9,2,1,0,01p i x x i =∈≠，p 为正整数，m 为整数。
8．有效数字的概念是什么？掌握有效数字与误差的关系。
答：若近似值x 的（绝对）误差限是它的某一位的半个单位，也就是说该近似值准确到这一位，且从该位起直到前面第一个非零数字为止的所有数字都称为有效数字。
若近似值x 的（绝对）误差限为n m x x e -⨯≤-=1021*，则称x 为具有n 位有效数字的有效数，或称它精确到nm -10位，其中的每一位数字n x x x ,,21都是x 的有效数字。
设精确值*x 的近似值x 的规格化形式为mp x x x x 10.021⨯±= ，若x 具有n 位有效数字，则其相对误差限为n r x e -⨯≤111021；反之，若x 的相对误差限为n r x e -⨯+≤1110)1(21，则x 至少有n 位有效数字。
9．下列各数都是对真值进行四舍五入后获得的近似值，试分别写出它们的绝对误差限，相对误差限和有效数字的位数。
（1）024.01=x （2）4135.02=x （3）50.573=x （4）600004=x （5）55108⨯=x ；解：（1）0005.01*1≤-=x x e ；0021.0*≤=-=xe x x x e r ；有三位有效数字。
（2）00005.02*2≤-=x x e ；000121.0*≤=-=xex x x e r ；有四位有效数字。
（3）005.03*3≤-=x x e ；000087.0*≤=-=xex x x e r ；有四位有效数字。
（4）5.04*4≤-=x x e ；0000084.0*≤=-=xex x x e r ；有五位有效数字。
（5）5.05*5≤-=x x e ；000000625.0*≤=-=xex x x e r ；有六位有效数字。
10．为了使19的相对误差≤0.1%，问至少应取几位有效数字？解：由19的首位数是4.设近似数*x 有n 位有效数字，由定理4.1可知，相对误差001.010421)(1*≤⨯⨯≤-n r x e ，解得097.3≥n ，即取4位有效数字，近似数的相对误差不超过0.1%。
11．已知33,3100,1150)(*2==-+==x x x x x P y ，计算)3100(*p y =及)33(P y =，并求x 和y 的相对误差。
解： 55555.51150)3100()3100()3100(2*-≈-+==p y 281150)33()33()33(2-=-+==P y333.0)(*≈-=x x x e 0101.0)()(≈=xx e x e r44444.22)(*≈-=y y y e 801587.0)()(≈=yy e y e r 12．写出误差估计的一般公式（以二元函数),(y x f z =为例）。
解：二元函数),(y x f z =的绝对误差： )(|)(|)(),(),(y e yfx e x f z e y x y x ⋅∂∂+⋅∂∂≈二元函数的相对误差： z y e y f z x e x f z z e z e y x y x r )(|)(|)()(),(),(⋅∂∂+⋅∂∂≈=)(|)(|),(),(y e yfz y x e x f z x r y x r y x ⋅∂∂⋅+⋅∂∂⋅=13．用电表测得一个电阻两端的电压和流过的电流范围分别为V V 2220±=，A I 1.010±=，求这个电阻的阻值R ，并估算其绝对误差和相对误差。
解：2)(≤V e ，1.0)(≤I e ，又2,1,IV I R I V R I V R -=∂∂=∂∂=。
所以： 42.01.01002202101)(|)(|)(|)(|)(),(),(),(),(=⨯+⨯=⋅∂∂+⋅∂∂≤⋅∂∂+⋅∂∂≈I e I R V e V R I e IRV e V R R e I V I V I V I V21099.1)()(-⨯≈=RR e R e r 。
14．若01.045.0,01.003.1*2*1±=±=x x ，计算22121x e x y +=的近似值，并估计)(y e 及其上界。
解：45.0221)03.1(e y +≈ )(21))(()21()21()(2*22*21*11*11*1*x x x x e e x x x x e x e x y y y e -++-=+-+=-=),(,01.0211006.2)(21))((*2221*11*12*2x x e e e x x x x x x ∈⨯⨯+⨯=-++-≤-ξξ15．已测得某场地长为m l 110=，宽d 的值为m d 80=，已知m l l l e 2.0)(*≤-=，m d d d e 1.0)(*≤-=，试求面积ld s =的绝对误差限和相对误差限。
解：由ld s =,l ds d l s =∂∂=∂∂,，m l l l e 2.0)(*≤-=，m d d d e 1.0)(*≤-=。
可得：301.0802.0110)(|)(|)(|)(|)(),(),(),(),(=⨯+⨯=⋅∂∂+⋅∂∂≤⋅∂∂+⋅∂∂≈d e ds l e l s d e d s l e l s s e d l d l d l d l 3104.3)()(-⨯≈=ss e s e r 。
16．掌握二元函数的加、减、乘、除和开方运算的绝对误差和相对误差估计公式。
解：（1）加、减运算： 由于()1/=∂+∂x y x ()()(),1/,1/,1/-=∂-∂=∂-∂=∂+∂y y x x y x y y x ，所以()()()()()()()()()()()()()()()()()()()()()|||/|||/|||,//,,//,y e y x y x e y x x y x e y e y x y x e y x x y x e y e x e y x e y e y x y x e y x x y x e y e x e y x e r r r r r r r r r ⋅-+⋅-≤-⨯--⨯-≈--≈-⨯++⨯+≈++≈+从而有（2）乘法运算： 由于()()，x yxy y x xy =∂∂=∂∂,所以()()()()()()y e x e xy e y xe x ye r r r +≈+≈,xy e ，从而()()()|||||||||y e x x e y xy e ⋅+⋅≤（3）除法运算： 由于2)(,1)(yx y y x y x yx-=∂∂=∂∂，所以)()(1)(2y e yxx e y y xe -≈，)()()(y e x e yxe r r r -≈（4）乘方及开方运算：由于()1-=∂∂n nnx xx ，所以()()()()x ne x e x e nx x e r n r n n ≈≈-,1 17．求方程01562=+-x x 的两个根，使它至少具有4位有效数字（982.27783≈）。
解：782.55982.272812114)56(5621=+≈⨯⨯⨯--+=x017863.0782.55112≈==x c x 19．求方程01162=+-x x 的较小正根，要求有3位有效数字。
解：937.15937.7812114)16(1621=+≈⨯⨯⨯--+=x062747.0937.15112≈==x c x 所以较小正根为062747.02≈x 。
20．设4110,,2,1,0, ==⎰n dx e x I x n n 。
（1）证明：4110,,2,1,0, =-=-n nI e I n n ；（2）给出一个数值稳定的算法，并证明算法的稳定性。
（1）证明：11111---=-===⎰⎰⎰n x n xn xn n nI e x d e nx e e d x dx e x I（2）)(11n n I e nI -=-设n n n I I e -=*，则n nnn n n n n n n e n I I e e nI I e e n I I e 1110*0022*221*11=-==-==-=------当n 无限大时，n e 越小，所以该算法稳定。