计算方法第一章习题

计算方法-刘师少版第一章课后习题完整答案

1
9000 m=1
9000.00
解（1）∵ 2.0004＝0.20004×10 ,
x − x ∗ = x − 0.20004 ≤ 0.000049 ≤ 0.5 × 10 −4
m-n=-4,m=1 则 n=5,故 x=2.0004 有 5 位有效数字
x1 =2，相对误差限 ε r =
1 1 × 10 −( n −1) = × 101−5 = 0.000025 2 × x1 2× 2
-2
（2）∵ －0.00200= -0.2×10 ,
m=-2
x − x ∗ = x − (−0.00200) ≤ 0.0000049 ≤ 0.5 × 10 −5
m-n=-5, m=-2 则 n=3,故 x=-0.00200 有 3 位有效数字
x1 =2，相对误差限 ε r =
4
1 × 101−3 =0.0025 2× 2
4 3 4 πR − π ( R * ) 3 3 ε r* (V ) = 3 4 3 πR 3 R 3 − (R* )3 ( R − R * )( R 2 + RR * + R * ) = = R3 R3 R − R * R 2 + RR * + R * R − R * R 2 + RR * + RR * = ⋅ ≈ ⋅ R R R2 R2
可以得到计算积分的递推公式：
I n = 1 − nI n −1
1 0
n = 1,2, L
1 0
I 0 = ∫ e x −1 dx = e x −1
则准确的理论递推式实际运算的递推式两式相减有
* *
= 1 − e −1
I n = 1 − nI n −1
* * In = 1 − nI n −1 * * * In − In = −n( I n −1 − I n −1 ) = − ne( I n −1 ) *

计算方法_课后习题答案

(4.5)(0.01172)

0.00879
（2）采用 Newton 插值多项式 y x N2(x) 根据题意作差商表：
i
xi
0
4
1
6.25
f (xi ) 2 2.5
一阶差商 2 9
2
9
3
2 11
二阶差商 4 495
N2 (7) 2 29 (7 4) ( 4 495) (7 4) (7 6.25) 2.6484848

1
e2
则根据二次Lagrange插值公式得：
L2 (x)

(x ( x0

x1)(x x2 ) x1)(x0 x2 )
y0

(x ( x1

x0 )(x x2 ) x0 )(x1 x2 )
y1

(x ( x2

x0 )(x x1) x0 )(x2 x1)
y2
2(x 1)(x 0.5) 2x(x 0.5)e1 4x(x 1)e0.5
8. 求作 f x xn1 关于节点 xi i 0,1, , n 的 Lagrange 插值多项式，并利用
插值余项定理证明
n
n
xin1li 0 1n xi
i0
i0
式中 li x 为关于节点 xi i 0,1, , n 的 Lagrange 插值基函数。
2 02 12 4 23 4 04 14 2 3
1 x2 3x 2 x 4 3x x2 6x 8 23 x x2 5x 4 1 x x2 3x 2
8
4
8

《数值计算方法》习题答案

《数值计算方法》课后题答案详解吉林大学第一章习题答案1. 已知(1)2,(1)1,(2)1f f f −===，求()f x 的Lagrange 插值多项式。

解：由题意知：()01201212001020211012012202121,1,2;2,1,1()()(1)(2)()()6()()(1)(2)()()2()()(1)(1)()()3(1)(2)(1)(2)()2162nj j j x x x y y y x x x x x x l x x x x x x x x x x l x x x x x x x x x x l x x x x x x x x L x y l x ==−=====−−−−==−−−−+−==−−−−−+−==−−−−+−==×+×−∴∑()2(1)(1)131386x x x x +−+×=−+2. 取节点01210,1,,2x x x ===对x y e −=建立Lagrange 型二次插值函数，并估计差。

解11201201210,1,;1,,2x x x y y e y e −−======1)由题意知：则根据二次Lagrange插值公式得：02011201201021012202110.510.520.51()()()()()()()()()()()()()2(1)(0.5)2(0.5)4(1)(224)(43)1x x x x x x x x x x x x L x y y y x x x x x x x x x x x x x x x x e x x e e e x e e x −−−−−−−−−−−−=++−−−−−−=−−+−−−=+−+−−+22)Lagrange 根据余项定理，其误差为(3)2210122()1|()||()||(1)(0.5)|3!61max |(1)(0.5)|,(0,1)6()(1)(0.5),()330.5030.2113()61()0.2113(0.21131)(0.21130.5)0.008026x f R x x e x x x x x x t x x x x t x x x x t x R x ξξωξ−+≤≤==−−≤−−∈′=−−=−+=−==≤××−×−=∴取并令可知当时，有极大值3. 已知函数y =在4, 6.25,9x x x ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。

数值计算方法课后习题答案

第一章绪论（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-5 极限运算法则

x3 2x 2 (4) lim 2 x 2 ( x 2) ( x 2) 2 穷大， 0, 根据无穷小的倒数是无解： lim 3 2 x 2 x 2 x x3 2x2 lim 2 x 2 ( x 2)
5x 4 x x 1 x 1 5 x 4 x 将分子、分母同乘 ( 5 x 4 x )( 5 x 4 x ) 解： lim lim 分子的有理化因子 x 1 x 1 x 1 ( x 1)( 5 x 4 x ) (5) lim
分子、分母 ( x 1)( x 2) x2 lim lim 1 2 2 约去零因子 x 1 (1 x)(1 x x ) x 1 1 x x
(3) lim
(n 1)(n 2)(n 3) 1. 3 n 5n 解：应用第44页定理 3 的(3)时，要求分子、分母的极限各自存在且分母的
解法二： ( 1 2 n 1 2 (n 1) ) lim
n
n(n 1) n(n 1) lim n 2 2 n(n 1) n(n 1) n(n 1) n(n 1) 2 2 2 2 lim n n(n 1) n(n 1) 2 2 n 1 lim lim n n(n 1) n(n 1) n 1 1 1 1 1 1 2 2 2 n 2 n
代数和所求极限是无限个具有 . 极限的函数的代数和相悖处为 .
1 n(n 1) 1 1 2 (n 1) 2 正确的解法： lim lim n n 2 n2 n2 3 1 3 1 (2) lim lim lim 0; 3 3 x 1 1 x 1 x x 1 1 x x 1 1 x 解：教材第44页定理3告诉我们：各自具有极限的函数的和、差、积、

数值计算方法课后习题答案(李庆扬等) (修复的)

，。

第一章绪论（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 。

计算方法_习题第一、二章答案

第一章误差1 问3.142，3.141，722分别作为π的近似值各具有几位有效数字？分析利用有效数字的概念可直接得出。

解 π=3.141 592 65…记x 1=3.142，x 2=3.141，x 3=722.由π- x 1=3.141 59…-3.142=-0.000 40…知3411110||1022x π--⨯<-≤⨯ 因而x 1具有4位有效数字。

由π- x 2=3.141 59…-3.141=-0.000 59…知2231021||1021--⨯≤-<⨯x π因而x 2具有3位有效数字。

由π-722=3.141 59 …-3.142 85…=-0.001 26…知231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。

2 已知近似数x*有两位有效数字，试求其相对误差限。

分析本题显然应利用有效数字与相对误差的关系。

解利用有效数字与相对误差的关系。

这里n=2,a 1是1到9之间的数字。

%5101101|*||)(|1211*=⨯≤⨯≤-=+-+-n rx x x ε3 已知近似数的相对误差限为0.3%，问x*至少有几位有效数字？分析本题利用有效数字与相对误差的关系。

解 a 1是1到9间的数字。

1112*10110113%3.0)(--⨯≤⨯=<=x r ε 设x*具有n 位有效数字，令-n+1=-1，则n=2，从而x*至少具有2位有效数字。

4 计算sin1.2，问要取几位有效数字才能保证相对误差限不大于0.01%。

分析本题应利用有效数字与相对误差的关系。

解设取n 位有效数字，由sin1.2=0.93…，故a 1=9。

411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式411101021-+-≤⨯n a 知取n=4即可满足要求。

5 计算76017591-，视已知数为精确值，用4位浮点数计算。

计算方法第一章绪论(32学时)-2014.2

教材聂玉峰、王振海等《数值方法简明教程》，高等教育出版社，2011作业计算方法作业集（A、B）参考书¾封建湖，车刚明计算方法典型题分析解集（第三版）西北工业大学出版社，2001¾封建湖，聂玉峰，王振海数值分析导教导学导考（第二版）西北工业大学出版社，2006¾车刚明，聂玉峰，封建湖，欧阳洁数值分析典型题解析及自测试题（第二版）西北工业大学出版社，2003西北工业大学理学院欧阳洁2第一章绪论§1 引言§2 误差的度量与传播§3 选用算法时应遵循的原则西北工业大学理学院欧阳洁3§1 引言科学与工程领域中运用计算机求解问题的一般过程：1 实际问题的提出2 建立数学模型3 设计可靠、高效的数值方法4 程序设计5 上机实践计算结果6 数据处理及结果分析西北工业大学理学院欧阳洁4学习算法的意义科学计算（数值模拟）已经被公认为与理论分析、实验分析并列的科学研究三大基本手段之一。

计算方法课程的研究对象具有广泛的适用性，著名流行软件如Maple、Matlab、Mathematica 等已将其绝大多数内容设计成函数，简单调用之后便可以得到运行结果。

但由于实际问题的具体特征、复杂性, 以及算法自身的适用范围决定了应用中必须选择、设计适合于自己特定问题的算法，因而掌握数值方法的思想和内容至关重要。

西北工业大学理学院欧阳洁5鉴于实际问题的复杂性，通常将其具体地分解为一系列子问题进行研究，本课程主要涉及如下几个方面问题的求解算法：¾非线性方程求根¾线性代数方程组求解¾函数插值¾曲线拟合¾数值积分与数值微分¾常微分方程初值问题的数值解法¾矩阵特征值与特征向量计算西北工业大学理学院欧阳洁6§2 误差的度量与传播一误差的来源与分类模型误差：数学模型与实际问题的误差观测误差：观测结果与实际问题的误差截断误差：数学模型的理论解与数值计算问题的精确解之间的误差舍入误差：对超过某有限位数的数据进行舍入所产生的误差西北工业大学理学院欧阳洁75 使用数值稳定性好的公式一个算法，如果初始数据微小的误差仅使最终结果产生微小的误差，或在运算过程中舍入误差在一定条件下能够得到控制，则称该算法（数值）稳定，否则称其为（数值）不稳定．西北工业大学理学院欧阳洁26总结1.数值运算的误差估计2.绝对误差、相对误差与有效数字3.数值运算中应遵循的若干原则西北工业大学理学院欧阳洁30。

计算方法(孙志忠)习题第一章绪论

第一章第一章绪论绪论绪论
一、填空填空（1）为了使计算32)
1(6)1(41310−−−+−+=x x x y 的乘除法运算次数尽量少，应将表达式改写为( )
1220.010s n x t x a a a ∗=×L ()准确值与其有位有效数字的近似值1(0)a x x ∗≠−≤的绝对误差 ( )。

（3）设 2.40315a =是真值 2.40194x =的近似值，则a 有（）位有效数字。

（4）设 2.001,0.8030a b ==−是由真值,x y 经四舍五入得来的近似值，则a b +的误差限为（）。

（5）为了减少舍入误差的影响，应将表达式62002008−改写为（）。

（6）用4位浮点数计算出表达式030.6228100.687610×+×的正确结果是（）
二、选择题
（1）
计算 61) 1.4α−≈，采用下列算式计算，问哪一个
得到的结果最好？（）
3A B C 32
D 、99－70、（－（2）若近似值a 的数字个数越多，则其相对误差越小。

（）
A、正确
B、错误
（4）若误差限为0.5×10－5，那么近似数0.003400有5位有效数字。

（）
A、正确
B、错误。

计算方法教程（第2版）习题答案

《计算方法教程(第二版)》习题答案第一章习题答案1、浮点数系),,,(U L t F β共有 1)1()1(21++---L U t ββ 个数。

3、a ．4097b ．62211101110.0,211101000.0⨯⨯c ．6211111101.0⨯4、设实数R x ∈，则按β进制可表达为：,1,,,3,2,011)11221(+=<≤<≤⨯++++++±=t t j jd d lt t d t t d dd x βββββββ按四舍五入的原则，当它进入浮点数系),,,(U L t F β时，若β211<+t d ，则 l tt d dd x fl ββββ⨯++±=)221()(若 β211≥+t d ，则 l tt d d d x fl ββββ⨯+++±=)1221()(对第一种情况：t l lt l t t d x fl x -++=⨯≤⨯+=-βββββ21)21(1)()(11对第二种情况：t l ltl t t d x fl x -++=⨯≤⨯--=-ββββββ21)21(1)(11 就是说总有： tl x fl x -≤-β21)( 另一方面，浮点数要求 β<≤11d ，故有l x ββ1≥，将此两者相除，便得t x x fl x -≤-121)(β 5、a ． 5960.1 b ． 5962.1 后一种准确6、最后一个计算式：00025509.0原因：避免相近数相减，避免大数相乘，减少运算次数7、a ．]!3)2(!2)2(2[2132 +++=x x x yb ．)21)(1(22x x x y ++=c ．)11(222-++=x x x yd ． +-+-=!2)2(!6)2(!4)2(!2)2(2642x x x x ye ．222qp p q y ++=8、01786.098.5521==x x9、 m )10(m f - 1 233406.0- 3 20757.0- 5 8.07 710计算宜采用：])!42151()!32141()!22131[()(2432+⨯-+⨯-+⨯--=x x x f第二章习题答案1、a ．T x )2,1,3(= b ．T x )1,2,1,2(--= c ．无法解2、a ．与 b ．同上， c ．T T x )2188.1,3125.0,2188.1,5312.0()39,10,39,17(321---≈---=7、a ．⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛---14112111473123247212122123211231321213122 b ． ⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫⎝⎛----333211212110211221213231532223522121⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎪⎭⎫⎝⎛-=111211212130213221219、T x )3415.46,3659.85,1220.95,1220.95,3659.85,3415.46(1= T x )8293.26,3171.7,4390.2,4390.2,3171.7,8293.26(2= 10、T LDL 分解：)015.0,579.3,9.1,10(diag D =⎪⎪⎪⎪⎪⎭⎫⎝⎛=16030.07895.05.018947.07.019.01L Cholesky 分解⎪⎪⎪⎪⎪⎭⎫⎝⎛=1225.01408.10833.15811.18918.12333.12136.23784.18460.21623.3G 解：)1,1,2,2(--=x 12、16,12,1612111===∞A A A611,4083.1,61122212===∞A A A2)(940)()(12111===∞A Cond A Cond A Cond524)(748)()(22221===∞A C o n d A C o n d A C o n d⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=--180.0000180.0000- 30.0000 180.0000- 192.0000 36.0000- 30.0000 36.0000- 9.0000,0.0139 0.1111- 0.0694- 0.1111- 0.0556 0.1111- 0.0694- 0.1111- 0.0139 1211A A1151.372,1666.0212211==--A A15、 1A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代不收敛； 2A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代不收敛； 3A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代收敛；第三章习题答案1、Lagrange 插值多项式：)80.466.5)(20.366.5)(70.266.5)(00.166.5()80.4)(20.3)(70.2)(00.1(7.51)66.580.4)(20.380.4)(70.280.4)(00.180.4()66.5)(20.3)(70.2)(00.1(3.38)66.520.3)(80.420.3)(70.220.3)(00.120.3()66.5)(80.4)(70.2)(00.1(0.22)66.570.2)(80.470.2)(20.370.2)(00.170.2()66.5)(80.4)(20.3)(00.1(8.17)66.500.1)(80.400.1)(20.300.1)(70.200.1()66.5)(80.4)(20.3)(70.2(2.14)(4--------⨯+--------⨯+--------⨯+--------⨯+--------⨯=x x x x x x x x x x x x x x x x x x x x x L Newton 插值多项式：)80.4)(20.3)(70.2)(00.1(21444779.0)20.3)(70.2)(00.1(527480131.0)70.2)(00.1(855614973.2)00.1(117647059.22.14)(4----+------+-+=x x x x x x x x x x x N2、设)(x y y =，其反函数是以y 为自变量的函数)(y x x =，对)(y x 作插值多项式：)1744.0)(1081.0)(4016.0)(7001.0(01253.0)1081.0)(4016.0)(7001.0(01531.0)4016.0)(7001.0(009640.0)7001.0(3350.01000.0)(----+---+--+--=y y y y y y y y y y y N 3376.0)0(=N 是0)(=x y 在]4.0,3.0[中的近似根。

黄云清版数值计算方法习题解答

第一章引论（习题）2．证明：记 x x f =)( ，则)()(***x x x x x xx x f E r +-=-=)(21**x E x x x x x xr ≈-⋅+=.3．证明：令： )()()(b a fl b a fl b a **-*=δ可估计： 1|)(|-≥*c b a fl β （c 为b a *阶码），故： 121||--≤c t c ββδt-=121β 于是： )1()()(δ+*=*b a b a fl .4．解 (1) )21()1(22x x x ++. (2))11(2x x x x x-++.(3) xxx x x x x cos 1sin )cos 1(sin cos 12+≈+=-.6．解 a 的相对误差：由于 31021|)(|-⋅≤-=a x x E . x a x x E r -=)(, 221018110921)(--⋅=⨯≤x E r . （1Th ） )(a f 对于)(x f 的误差和相对误差. |11||)(|a x f E ---==()25.021011321⨯⋅≤-+---ax x a =310-33104110|)(|--⨯=-≤a f E r .9．解递推关系： 1101.100-+-=n n n y y y (1) 取初值 10=y ， 01.01=y 计算可得： 11001.10022-⨯=-y 10001.1-=410-= 6310-=y ， 8410-=y ， 10510-=y ， …(2) 取初值 50101-+=y ， 2110-=y , 记： n n n y y -=ε,序列 {}n ε ，满足递推关系，且 5010--=ε ， 01=ε1101.100-+-=n n n εεε, 于是： 5210-=ε,531001.100-⨯=ε, 55241010)01.100(---⨯=ε,55351002.20010)01.100(--⨯-⨯=ε, 可见随着 n ε 的主项 5210)01.100(--⨯n 的增长，说明该递推关系式是不稳定的.第二章多项式插值 (习题)1. 方法一. 由 Lagrange 插值公式)()()()()(332211003x l f x l f x l f x l f x L ⋅+⋅+⋅+⋅=)1)((31)2)()(1()1)(()(123210---=-----=x x x x x x x l , ))(1(2)1)()(1()(21221211--=--+=x x x x x x l ， x x x x x x l )1()()1()1!()(2382121232--=-⋅⋅-+=, )()1(12)()1()(2121213-+=⋅⋅-+=x x x x x x x l . 可得： )21()(23-=x x x L方法二. 令：)()21()(3B Ax x x x L +-=由 23)1(3-=-L ， 21)1(3=L ，定A ，B （称之为待定系数法）2. 证明(1) 由于 j i j i x l ,)(δ= 故： =)(x L n ∑=ni i k i x l x0)( ，当 j x x = 时有： k j j n x x L =)( ， n j ,,1,0 =)(x L n 也即为 kx 的插值多项式，由唯一性，有：∑==ni k i ki x x l x)( ， n k ,,1,0 =证明(2)：利用Newton 插值多项式)(],[)()(0100x x x x f x f x N n -+=)()(],,[100---++n n x x x x x x f )()()()()()(00101x l x x x x x x x x x f n n =----=差商表：f(x) 一阶二阶 … n 阶差商0x 1 1x 0101x x -)()(11020x x x x --n x 0 0)()(1010n x x x x --代入)(*式有：)()()()()(1)(020*******n n n x x x x x x x x x x x x x x x N -----++--+=- . )(0x l 为n 次代数多项式，由插值多项式的唯一性：有 )()(0x N x l n ≡.4. 解作)(x f 以b a a ,,ε+为节点的Lagrange 插值多项式，有： )()()(22x R x L x f +=，其中：)()()()()()()()()()(2εεεεε+-+--+-----=a fb a b x a x a f b a b x a x x L)()()()()(b f a b a b a x a x εε------+，)()()(!3)()(2b x a x a x f x R ----'''=εζ ， b a <<ζ 令： 0→ε 有 )()(6)()()(22b x a x f x R x R --'''=→ζ，又：)()()()([)()(2a f a b ax a f a b a x x b x L εεεεε----+----= )]()()()()(a f a b a x a f a b a x -------+εεεε )()()()()(b f a b a b a x a x εε------+)()()2()(2a f ab a b x x b --+-→)()()()(a f a b a x x b '---+ )()()()(22x P b f a b a x =--+ 故当 0→ε 时，成立公式： )()()(x R x P x f +=.5. 解：因为34)(3'-=x x f ，2''12)(x x f =)(x f 为凹函数.又从数值表可见：当]5.0,1.0[∈x 时，)(x f 单调下降.有反函数)(1y fx -=)(y f的Newton 插值多项式：)17440.0)(10810.0)(40160.0)(70010.0(01225.0)10810.0)(40160.0)(70010.0(01531.0)40160.0)(70010.0(0096436.0)70010.0(33500.01.0)(4+---+------+--=y y y y y y y y y y y N.337.0)0(4*≈=N x7．解 1)(37++=x x x f .有：=]2,,2,2[71f !7)()7(ξf =1, !8)(]2,,2,2[)8(810ηf f = 0=.9. 证明：(1) =⋅-⋅=⋅∆++i i i i i i g f g f g f 11)(i i i i i i i i g f g f g f g f ⋅-⋅+⋅-⋅++++1111i i i i f g g f ∆+∆=+1.(3) n x n n)1()1(-=∆！)()(nh x h x x h n ++此题可利用数学归纳法：设 k n = 成立，证明 1+=k n 成立.又 1=n 时是成立的.10. 证明：记： 2]2/)1([)(+=n n n f ，33321)(n n g +++=有： 3)1()()1()(+=-+=∆n n f n f n f 故： ∑-=∆=10)()(n k k f n g ∑-=-+=1)]()1([n k k f k f2]2/)1([)0()(+=-=n n f n f .13．解作重节点差商的Newton 插值公式)1(]1,1[)1()(+--+-=x f f x P 22)1(]1,0,1,1[)1(]0,1,1[+--++--+x x f x f )1()1(]1,1,0,1,1[2-+--+x x x f 重节点差商表：i x i f 一阶二阶三阶四阶10-=x 110-=x 1 201=x 1 0 -212=x 1 0 0 112=x 1 2 2 1 0得 22)1()1(2)1(21)(+++-++=x x x x x P 13+-=x x .17．证：取 ,00=x 211=x ， 12=x ， 21=h00=f ， 11=f ， 12=f 记： )(i i x s M ''= ， 2,1,0=i有 hx x M h x x M x S 01101)(-+-=''x M x M 102)21(2+-= )21(2)1(2)(212-+-=''x M x M x S 又三弯矩方程为：( 2],,[210-=x x x f )244210-=++M M M ， )24(41201M M M ++-=.分段积分：⎰⎰+''=''∆1021221)]([)]([dx x s dx x s ⎰''12221)]([dx x s ⎰+-+=21201)]21([4dx x M x M ⎰-+-121221)]21()1([4dx x M x M⎰⎰-+-+-+-=121121221201)]21()1([4)]1()21([4dxx M x M dx x M x M由于 ⎰=-1212241)21(dx x ，⎰=-1212241)1(dx x ，⎰=--121481)1()21(dx x x ，于是：⎰++++=''∆1022212110202]2[61))((M M M M M M M dx x S 又： )24(41201M M M ++-=记 =),(20M M I ⎰∆''12))((dx x S=)()24(41[6120202220M M M M M M +++-+ ])24(81220M M +++由00=∂∂M I， 02=∂∂M I . 得：⎩⎨⎧=+-=-07072020M M M M 即当： 020==M M 时, ),(20M M I 达最小故：⎰=⋅⋅≥''∆102212)24(8161))((dx x S ,由最小模原理： ⎰≥''1212)]([dx x f .20. 解利用三弯矩方法 )(i i x s M ''= ， 2,1,0=i 10=x ， 22=x ， 32=x⎪⎩⎪⎨⎧-=+=++=+542364622121010M M M M M M M解得： 70-=M ， 201=M ， 372-=M]2,1[∈x 72431729)(231-+-=x x x x s ]3,2[∈x 105229367219)(232+-+-=x x x x s .第三章最佳逼近及其实现 (习题)2. 解 (1) ⎰'⋅'=badx x g x f g f )()(),( 不是 ),(b a c '中的内积，事实上容易验证：),(),(f g g f = ， ),(),(g f g f λλ= ),(),(),(w g w f w g f +=+但是 0),(=f f 当且仅当 0)(≡x f . 条件不满足，因为： ⎰='⋅'=badx x f x f f f 0)()(),(推出0)(≡'x f ，0)(≠=const x f . 因而 ),(g f 不是 ),(b a C '中的内积.(2) ),(g f 是 =],[10b a C {}],[)(,0)(:)(b a C x f a f x f '∈'=空间的内积，这是因为： 0),(=f f 推出 0)(='x f , C x f =)(，又],[10b a C f ∈ ，故 0)(=x f .4. 解：由于 0)(],,[2≠''∈x f b a c f ，则)(x f ''于],[b a 上保号，由定理5的推论2可知：)()(1x P x f -的交错点组恰有三个交错点，且 a x =1，b x =3，即： ⎪⎪⎩⎪⎪⎨⎧=-'='-=+-==+-==+-=0)()(,)()()(,)()()(,)()()(122210223103311011αρααρααρααx f x e x x f x e x x f x e x x f x e 故： a b a f b f x f --='=)()()(21α，2)()(2)()(220x a a b a f b f x f a f +⋅---+=α 记 c x =2 ，即证得(1).(2) 若 x x f cos )(= ，]2,0[],[π=b a此时由 ab a f b fc f --=')()()( 得：π2sin =c ， )2sin(πarc c =，πα21-=πππα2)4(2120-+=2)/2sin(2ππarc ⋅+)4(212-+=πππππ)2sin(arc +. 误差估计：)()(10b f b f E -+=-=ααρ)4(212-+=πππ1)2sin(-+ππarc5. 解：选取α ，使得：=)(αI ||max 211x x x α-≤≤ ，达到极小，即要求 x x *)(*αϕ= ，于]1,0[上一致逼近于2x ，如图应选 *α ，使得：x x x *)(2αϕ-=,于 ]1,0[ 上有两个轮流为正负偏差点，其中之一为1，另一个假设为 ζ 于是： )()1(ζϕα-=, 0)(='ζϕ ，（ ζ为)(x ϕ的极值点）得： αζζα+-=-2102=-αζ 解得：ζα2= ，0122=-+ζζ, 212,1±-=ζ取12-=ζ ， 222-=α. 又： α 是唯一的.6. 证明：由最佳一致逼近的特征定理，)(*x P n 为)(x f 的最佳一致逼近多项式，则存在2+n 个点b x x x a n ≤<<<≤+110使得： )()()(*k n k k x P x f x e -==*)1(n kP f --σ.又由于 ],[)(b a C x f ∈ ，于 ),(1+i i x x 中有一个点 i η ，1+<<i i i x x η ，使得： 0)()()(*=-=i n i i P f e ηηη, n i ,,1,0 =即： )(*x P n 为)(x f 满足插值条件： )()(*i i n f P ηη= ， n i ,,1,0 = 的插值多项式.7．解：求C*，使得：C x f C I bx a R C -=≤≤∈)(max min *)(记 C x f x e -=)()(, 依最佳一致逼近的特征定理：应取 )](min )(max [21*],[],[x f x f C b a b a +=*)()(C x f x e -=于 ],[b a 才有两个轮流正负的偏差点，（即 )(x f 于],[b a 上的最大值点和最小值点）1x ，2x )(max )(],[1x f x f b a = ， )(min )(],[2x f x f b a =此时： *)(m a x )1()(],[C x f x e b a ii --=σ即 *C 为)(x f 的零次最佳逼近多项式.8. 解： 436)(23+++=x x x x f 2)(34)3(62031T T T T +++=014T T ++01232112112323T T T T +++= 因为)(413x T 与零偏差最小，故： 012221121123)(T T T x P ++=421132++=x x . 为)(x f 的最佳一致逼近多项式.9. 证明：我们仅证明)(x f 是偶函数时，)(x P n 亦是偶函数.由于)(x P n 为)(x f的最佳一致逼近多项式，有：)()()(max ],[f E x P x f n n a a =--和： [,max ()()()]n n a af x P x E f ----=即： )()()(m a x ],[f E x P x f n n a a =---)(x P n -亦是)(x f 的最佳一致逼近多项式，由最佳一致逼近多项式的惟一性，有： )()(x P x P n n =-即： )(x P n 为偶函数.11．解：设 x a a x P 10*1)(+= ， 2210*2)(x b x b b x P ++= 分别为)(x f 的一次、二次最佳平方逼近多项式。

数值计算方法》习题答案

《数值计算方法》课后题答案详解吉林大学第一章习题答案1. 已知(1)2,(1)1,(2)1f f f −===，求()f x 的Lagrange 插值多项式。

解：由题意知：()01201212001020211012012202121,1,2;2,1,1()()(1)(2)()()6()()(1)(2)()()2()()(1)(1)()()3(1)(2)(1)(2)()2162nj j j x x x y y y x x x x x x l x x x x x x x x x x l x x x x x x x x x x l x x x x x x x x L x y l x ==−=====−−−−==−−−−+−==−−−−−+−==−−−−+−==×+×−∴∑()2(1)(1)131386x x x x +−+×=−+2. 取节点01210,1,,2x x x ===对x y e −=建立Lagrange 型二次插值函数，并估计差。

解11201201210,1,;1,,2x x x y y e y e −−======1)由题意知：则根据二次Lagrange插值公式得：02011201201021012202110.510.520.51()()()()()()()()()()()()()2(1)(0.5)2(0.5)4(1)(224)(43)1x x x x x x x x x x x x L x y y y x x x x x x x x x x x x x x x x e x x e e e x e e x −−−−−−−−−−−−=++−−−−−−=−−+−−−=+−+−−+22)Lagrange 根据余项定理，其误差为(3)2210122()1|()||()||(1)(0.5)|3!61max |(1)(0.5)|,(0,1)6()(1)(0.5),()330.5030.2113()61()0.2113(0.21131)(0.21130.5)0.008026x f R x x e x x x x x x t x x x x t x x x x t x R x ξξωξ−+≤≤==−−≤−−∈′=−−=−+=−==≤××−×−=∴取并令可知当时，有极大值3. 已知函数y =在4, 6.25,9x x x ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。

计算方法习题及答案

第一章绪论一.填空题1.*x 为精确值x 的近似值；()**x f y=为一元函数()x f y =1的近似值；()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-：***rx x e x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅ ()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂2、计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫舍入误差。

3、分别用2.718281，2.718282作数e的近似值，则其有效数字分别有 6 位和 7 位；又取1.73≈（三位有效数字），则-211.73 10 2≤⨯。

4、设121.216, 3.654x x ==均具有3位有效数字，则12x x 的相对误差限为 0.0055 。

5、设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为 0.01 。

6、已知近似值 2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为 0.000021 .7、递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取0 1.41y ≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定不稳定 .8、精确值 14159265.3*=π，则近似值141.3*1=π和1415.3*2=π分别有 3 位和 4 位有效数字。

9、若*2.71828x e x =≈=,则x 有 6 位有效数字，其绝对误差限为1/2*10-5。

10、设x*的相对误差为2％，求(x*)n 的相对误差0.02n 二、计算题1. 有一个长方形水池,由测量知长为(50±0.01)米,宽为(25±0.01)米,深为(20±0.01)米,试按所给数据求出该水池的容积,并分析所得近似值的绝对误差和相对误差公式,并求出绝对误差限和相对误差限. 解：设长方形水池的长为L ，宽为W,深为H ，则该水池的面积为V=LWH当L=50,W=25,H=20时，有 V=50*25*20=25000(米3) 此时，该近似值的绝对误差可估计为()()()()()()()=V V VV L W H L W HWH L HL W LW H ∂∂∂∆≈∆+∆+∆∂∂∂∆+∆+∆ 相对误差可估计为：()()r V V V∆∆=而已知该水池的长、宽和高的数据的绝对误差满足()()()0.01,0.01,0.01L W H ∆≤∆≤∆≤故求得该水池容积的绝对误差限和相对误差限分别为()()()()()()325*20*0.0150*20*0.0150*25*0.0127.5027.501.1*1025000r V WH L HL W LW H V V V -∆≤∆+∆+∆≤++=∆∆=≤=2.已知测量某长方形场地的长a=110米,宽b=80米.若()()**0.1 0.1a a b b -≤-≤米，米试求其面积的绝对误差限和相对误差限. 解：设长方形的面积为s=ab当a=110,b=80时，有 s==110*80=8800(米2) 此时，该近似值的绝对误差可估计为()()()()()=b s ss a b a ba ab ∂∂∆≈∆+∆∂∂∆+∆ 相对误差可估计为：()()r s s s∆∆=而已知长方形长、宽的数据的绝对误差满足()()0.1,0.1a b ∆≤∆≤故求得该长方形的绝对误差限和相对误差限分别为()()()()() 80*0.1110*0.119.019.00.0021598800r s b a a b s s s ∆≤∆+∆≤+=∆∆=≤= 绝对误差限为19.0；相对误差限为0.002159。

初一数学第一章练习题

初一数学第一章练习题在初一的数学学习中，第一章是一个重要的里程碑，涵盖了一些基础的数学知识和技巧。

为了巩固所学的内容，我们需要进行一些练习题。

以下是一些针对初一数学第一章的练习题，希望能够帮助同学们更好地理解和掌握所学的知识。

练习题1：整数的加减1. 计算：（-5）+ 8 - 12 + 32. 计算：7 - 2 + 33. 计算：（-3）- 7 + 4 - 9练习题2：整数的乘除1. 计算：（-2）×（-5）2. 计算：6 ÷ 23. 计算：（-15） ÷ 5练习题3：小数的加减乘除1. 计算：0.6 + 0.352. 计算：2.1 - 0.93. 计算：0.4 × 0.24. 计算：1.5 ÷ 0.51. 计算：（-4）²2. 计算：√163. 计算：√(0.01)练习题5：图形的面积和周长1. 已知正方形的一条边长为4cm，求它的面积和周长。

2. 一个长方形的长为6m，宽为3m，求它的面积和周长。

3. 一个圆的直径为8cm，求它的面积和周长。

练习题6：比例和百分数1. 用比例的方法解决下列问题：某班男生占全班人数的5/8，如果男生有24人，那么全班人数是多少？2. 用百分数的方法解决下列问题：一块售价40元的物品打折50%，打完折后的价格是多少元？练习题7：代数表达式1. 将以下自然语言表示的表达式转换为代数表达式：- 两个数相加的结果- 一个数减去另一个数的两倍- 一个数的平方减去41. 解方程：2x + 3 = 92. 解方程：4x - 7 = 5x + 13. 解不等式：3x + 5 > 2x + 8以上是一些初一数学第一章的练习题，希望同学们能够认真尝试解答，巩固所学的知识。

如果遇到困难，可以寻求同学、老师或家长的帮助。

通过勤奋练习，我们相信大家能够取得进步，为接下来的数学学习打下坚实的基础。

祝大家学习进步！。