李庆扬数值分析第五版习题答案清华大学出版社
李庆扬 数值分析第五版 习题答案
第2章 复习与思考题01ii i ii kx x x x 的基函数称为主要性质有 0,()1,k i kx i k()1n l x、什么是牛顿基函数？它与单项式基答：牛顿差值基函数为00101x ),(x x )(x x ),...,(x x )(x x )...(x x )}n 牛顿差值基函数中带有常数项01,,...n x x x ，这有单项式基不同。
阶均差？它有何重要性质 01n 2n 01n 2n -11[,,...,,][,,...,,]n n f x x x f x x x x x xk j 0j 0j-1j j+1j -k x x x x x x x （）...()()...()和k 阶均差的性质0101k-10[,,...,][,,...,]k kf x x x f x x x x x （分子前项多xk ）[a,b]上存在阶导数，且节点2n ,[a,b]x ，则1()!f n0()nn n ik k kk k i i ki kx x y l x y x x ，(j 1,2,....,n)个点的牛顿插值多项式01[,,...,]k f x x x ,(k 1,2,....,n)两者的主要差异是未知数不一致。
拉格朗日插值多项式是系数知道，但基函数不知道。
牛顿插值多项式是函数知道，但系数不知道。
与一般多项式基本相同。
y ，其中系数矩阵用下列基底作多项式插值时，120001211112222121...1...1 (1)...n n n n n nnx x x x x x x x x x x x ，无非零元素。
）拉格朗日基底为01{(),(),...,()}n l x l x l x ，已知数为未知数为01{(),(),...,()}n l x l x l x ，则系数矩阵为00101x ),(x x )(x x ),...,(x x )(x x )...(x x )}n ，已，未知数为012{,,,...,}n a a a a ，则系数矩阵为102020211010100...010...01()()...0...............1()()...()n nnnnj j x x x x x x x x x x x x x x x x ，为下三角矩阵，矩阵的上三角元0。
数值分析上机(四)
Ps ：题目均来自数值分析第五版作者：李庆扬，王能超，易大义 编 出 版 社：清华大学出版社误差分析问题：求下列方程的实根 （1）2x 320xx e -+-=（2）322+10x 200xx +-=要求：（1）设计一种不动点迭代法，要使迭代序列收敛，然后再用斯特芬森加速迭代，计算到81|x x |10k k ---<为止。
（2）用牛顿迭代，同样计算到81|x x |10k k ---<，输出迭代初值及各次迭代值和迭代次数k ，比较方法的优劣。
代码部分： /**函数**/function y = fun(x) y=x^3+2*x^2+10*x-20; endfunction y = fun1( x) y=x^2-3*x+2-exp(x); % y=2*log(x)+log(3); endfunction [y,k]= niudun(x0)%NUIDUN Summary of this function goes here % Detailed explanation goes here x(1)=x0; k=1; des=1;while des>1.0e-8x(k+1)=x(k)-fun1(x(k))/dfun1(x(k)); des=abs(x(k+1)-x(k)); k=k+1; end y=x(k); k=k; endfunction [y,k]= sitefensen(x0,f)%SITEFENSEN Summary of this function goes here % Detailed explanation goes here%x0为初值，n 为迭代次数，f 为迭代函数 x(1)=x0;des=1;k=1;while des>1.0e-8y(k)=f(x(k));z(k)=f(y(k));x(k+1)=x(k)-(y(k)-x(k))^2/(z(k)-2*y(k)+x(k));des=abs(x(k+1)-x(k));k=k+1;endy=x(k);k=k;end%fun的导数function y= dfun(x)%DFUN Summary of this function goes here % Detailed explanation goes herey=3*x^2+ 4*x+10;end%fun1的导数function y = dfun1( x )%DFUN1 Summary of this function goes here % Detailed explanation goes herey=2*x-exp(x)-3;endclearclc%不动点迭代法n=100;x0=0.5;%初值k=0;des=1;while des>1.0e-8x=(x0^2+2-exp(x0))/3;% 2*log(x)+log(3)% x=(-x0^3-10*x0+20)/2*(x0+eps);des=abs(x-x0);k=k+1;x0=x;enddisp('不动点迭代解->')fprintf('%f\n',x)disp('迭代次数->')fprintf('%d\n',k)disp('误差->')fprintf('%f\n',abs(0-fun1(x))) %斯蒂芬森加速f='(x^2+2-exp(x))/3';f=inline(f);[yy,kk]=sitefensen(0.5,f);disp('斯蒂芬森加速解->') fprintf('%f\n',yy)disp('迭代次数->')fprintf('%d\n',kk)disp('误差->')fprintf('%f\n',abs(0-fun1(yy))) %牛顿迭代法[yy1,kk1]=niudun(0.5);disp('牛顿迭代法解->')fprintf('%f\n',yy1)disp('迭代次数->')fprintf('%d\n',kk1)disp('误差->')fprintf('%f\n',abs(0-fun1(yy1)))第一个方程的运行结果如下：不动点迭代解->0.257530迭代次数->14误差->0.000000斯蒂芬森加速解->0.257530迭代次数->5误差->0.000000牛顿迭代法解->0.257530迭代次数->5误差->0.000000结论：由上述三个表格可以看出去在迭代初值相同的情况下，斯蒂芬森加速和牛顿加速迭代次数都明显少于不动点迭代。
数值分析课程第五版课后习题答案(李庆扬等)1之欧阳引擎创编
第一章 绪论（12）欧阳引擎（2021.01.01）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 xxx 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 nk k kεεεε；（3）*4*2/x x 。
数值分析课程第五版课后习题答案(李庆扬等)
数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章：数值分析导论1. 解答：数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题，其结果可能不是完全精确的，但是能够满足工程或科学应用的要求。
2. 解答：数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础，对许多领域都有着广泛的应用，如物理学、经济学、生物学等。
3. 解答：数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比，数值方法一般更加灵活和高效，可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答：计算误差是指数值计算结果与精确解之间的差异。
在数值计算中，由于计算机的有限精度以及数值计算方法本身的近似性等因素，都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章：数值误差分析1. 解答：绝对误差是指实际值与精确值之间的差异。
例如，对于一个计算出的数值近似解x和精确解x_0，其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度，通常被用作评估数值计算方法的好坏。
2. 解答：相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0，其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度，常用于评估数值计算方法的收敛速度。
3. 解答：舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数，因此在进行数值计算过程中，舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答：误差限度是指对于给定的数值计算问题，所能容忍的误差范围。
李庆扬-数值分析第五版第5章与第7章习题答案
（3）一个单位下三角矩阵的逆仍为单位下三角矩阵。
答：正确。
（4）如果A非奇异，则Ax = b的解的个数是由右端向量b的决定的。
答：正确。解释：若A|b与A的秩相同，则A有唯一解。若不同，则A无解。
（5）如果三对角矩阵的主对角元素上有零元素，则矩阵必奇异。
（6）范数为零的矩阵一定是零矩阵。
答：正确。
17、大熊座的明显标志就是我们熟悉的由七颗亮星组成的北斗七星，
11.判断下列命题是否正确：
8、铁生锈的原因是什么？人们怎样防止铁生锈？（1）非线性方程（或方程组）的解通常不唯一（正确）
（2）牛顿法是不动点迭代的一个特例（正确）
11、月食：当地球转到月球和太阳的中间，太阳、地球、月球大致排成一条直线时，地球就会挡住太阳射向月球的光，这时在地球上的人就只能看到月球的一部分或全部看不到，于是就发生了月食。（3）不动点迭代法总是线性收敛的（错误）
4从以上可以看出，每次运算后，区间长度减少一半，是线形收敛。
3.什么是函数 的不动点？如何确定 使它的不动点等价于 的零点
P215.
将方程 改写成等价的形式 ，若要求 满足 ，则 ；反之亦然，称 为函数 的一个不动点。
4.什么是不动点迭代法？ 满足什么条件才能保证不动点存在和不动点迭代序列收敛于 的不动点
齐次性
三角不等式
设 为向量，则三种常用的向量范数为：（第3章p53,第5章p165）
7、何谓矩阵范数？何谓矩阵的算子范数？给出矩阵A =(ai j)的三种范数||A||1，||A||2，||A||∞，||A||1与||A||2哪个更容易计算？为什么？
向量范数定义见p162，需要满足四个条件。
正定条件
均是对称正定矩阵的性质。应予以记住。
数值分析课程第五版课后习题答案(李庆扬等)1之欧阳美创编
第一章 绪论（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 n x n 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 n k 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之欧阳与创编
第一章 绪论（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 xxx x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x n x n 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 n k 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 nk k kεεεε；（3）*4*2/x x 。
李庆扬-数值分析第五版第5章与第7章习题答案
3）设 ，则 ，从而 ，所以迭代方法发散。
4）设 ，则 ，从而
，所以迭代方法发散。
第5章
复习与思考题
1、用高斯消去法为什么要选主元？哪些方程组可以不选主元？
答：使用高斯消去法时，在消元过程中可能出现 的情况，这时消去法无法进行；即时主元素 ，但相对很小时，用其做除数，会导致其它元素数量级的严重增长和舍入误差的扩散，最后也使得计算不准确。因此高斯消去法需要选主元，以保证计算的进行和计算的准确性。
4从以上可以看出，每次运算后，区间长度减少一半，是线形收敛。
3.什么是函数 的不动点？如何确定 使它的不动点等价于 的零点
P215.
将方程 改写成等价的形式 ，若要求 满足 ，则 ；反之亦然，称 为函数 的一个不动点。
4.什么是不动点迭代法？ 满足什么条件才能保证不动点存在和不动点迭代序列收敛于 的不动点
因为C的一、二、三阶顺序主子式分别为1，5，1，所以C能够分解为三角阵的乘积，并且分解是唯一的。
12、设
，
计算A的行范数，列范数，2-范数及F-范数。
本题考查的是矩阵范数的定义及求法
行范数0.6+0.5=1.1
列范数0.5+0.3=0.8
2-范数的计算需要用到特征值，特征值的计算可以使用幂法进行计算，也可以直接求。
（1）计算 的递推公式
（2）解Ly=f
（3）解UX=y
10、用改进的平方根法解方程组
。
本题明确要求使用平方根法进行求解。实际考查的LDU分解。见P157
。
11、下列矩阵能否分解为 （其中L为单位下三角阵，U为上三角阵）？若能分解，那么分解是否唯一。
， ， 。
LU分解存在的条件
李庆扬数值分析第五版习题答案清华大学出版社
2
( A*) 2 A * ( x*) .
当 x* 100 时，若 ( A*) 1, 则 ( x*)
1 102 2
2
故测量中边长误差限不超过 0.005cm 时，才能使其面积误差不超过 1cm 10．设 S
1 2 gt ，假定 g 是准确的，而对 t 的测量有 0.1 秒的误差，证明当 t 增加时 S 的 2
1 dx N 1 x2
N 1
arctan(tan( )) tan tan arctan 1 tan tan N 1 N arctan 1 ( N 1) N 1 arctan 2 N N 1
9．正方形的边长大约为了 100cm，应怎样测量才能使其面积误差不超过 1cm ？ 解：正方形的面积函数为 A( x) x
2
13． f ( x) ln( x x 1) ,求 f (30) 的值。若开平方用 6 位函数表，问求对数时误差有多 大？若改用另一等价公式。 ln( x x 1) ln( x
2
x 2 1)
计算，求对数时误差有多大？ 解
f ( x) ln( x x 2 1) , f (30) ln(30 899)
( y10 *) 1010 ( y0 *) 1 1010 102 2 1 108 2
计算到 y10 时误差为
1 108 ，这个计算过程不稳定。 2
6
12．计算 f ( 2 1) ，取 2 ，利用下列等式计算，哪一个得到的结果最好？
1 , ( 2 1)6
x0 1, x1 1, x2 2, f ( x0 ) 0, f ( x1 ) 3, f ( x2 ) 4; l0 ( x) l1 ( x) l2 ( x ) ( x x1 )( x x2 ) 1 ( x 1)( x 2) ( x0 x1 )( x0 x2 ) 2 ( x x0 )( x x2 ) 1 ( x 1)( x 2) ( x1 x0 )( x1 x2 ) 6 ( x x0 )( x x1 ) 1 ( x 1)( x 1) ( x2 x0 )( x2 x1 ) 3
李庆扬数值分析第五版习题答案解析清华大学出版社
又
即计算值比准确值大。
故 在 内至少有三个互异零点，
依此类推， 在 内至少有一个零点。
记为 使
又
其中 依赖于
分段三次埃尔米特插值时，若节点为 ，设步长为 ，即
在小区间 上
16．求一个次数不高于4次的多项式P（x），使它满足
解：利用埃米尔特插值可得到次数不高于4的多项式
设
其中，A为待定常数
从而
17．设 ，在 上取 ，按等距节点求分段线性插值函数 ，计算各节点间中点处的 与 值，并估计误差。
19。观测物体的直线运动，得出以下数据：
时间t(s)
0
0.9
1.9
3.0
3.9
5.0
距离s(m)
0
10
30
50
80
110
求运动方程。
解：
被观测物体的运动距离与运动时间大体为线性函数关系，从而选择线性方程
令
则
则法方程组为
从而解得
故物体运动方程为
20。已知实验数据如下：
19
25
31
38
44
19.0
32.3
将 代入得
由此得矩阵开工的方程组为
求解此方程组，得
又 三次样条表达式为
将 代入得
21．若 是三次样条函数，证明：
若 ，式中 为插值节点，且 ，则
证明：
从而有
第三章 函数逼近与曲线拟合
1． ，给出 上的伯恩斯坦多项式 及 。
解：
伯恩斯坦多项式为
其中
当 时，
当 时，
2．当 时，求证
证明：
若 ，则
3．证明函数 线性无关
解：
采用复化梯形公式时，余项为
数值分析上机(八)
Ps：均来自数值分析第五版作者：李庆扬，王能超，易大义编出版社：清华大学出版社幂法求特征值特征向量问题：给定矩阵A=[5 4 1 14 5 1 11 1 4 21 12 4].(1)用幂法求A的主特征值及对应的特征向量，并用瑞丽商加速法观察加速效果。
(2)利用反幂法迭代，试用不同p值，求A的不同特征值及特征向量，比较结果。
代码部分：/*mifa.m*/clearclc%给定实对称阵A=[5 4 1 14 5 1 11 1 4 21 12 4];n=length(A);%任取初值v0=rand(n,1);k=100;v(:,1)=v0;for i=1:k+1v(:,i+1)=A*v(:,i);v(:,i)=v(:,i+1);end%内部函数求得的最大特征值t1=max(eig(A))%幂法求最大特征值t2=v(:,k+1)./v(:,k)%瑞利商加速法%任意向量xl=rand(n,1);%迭代次数kk=100;for i=1:kku(:,i)=A*xl./max(A*xl);xl=A*xl;endfor k=1:10y(k)=sum(u(:,k)'*(A*u(:,k)))./sum(u(:,k).^2);endy运行结果：t1 =10t2 =10101010y =Columns 1 through 79.7613 9.9487 9.9873 9.9968 9.9992 9.999810.0000Columns 8 through 1010.0000 10.0000 10.0000结果分析：调用内部函数命令eig算得最大特征值和幂法所求一样，瑞丽商加速法在第八步收敛到所求的特征值。
/*fanmifa.m*/%% 反幂法clearclc%给定实对称阵A=[5 4 1 14 5 1 11 1 4 21 12 4];n=length(A);%任取初值v0=rand(n,1);k=100;v(:,1)=v0;for i=1:k+1v(:,i+1)=A*v(:,i);v(:,i)=v(:,i+1);end%内部函数求得的最小特征值t0=min(eig(A))%任取向量u0=rand(n,1);N=100;u(:,1)=u0;for k=2:Nv(:,k)=A\u(:,k-1);u(:,k)=v(:,k)/max(v(:,k)); endtmin=1/max(v(:,N))运行结果：t0 =1.0000tmin =1.0000结果分析：反幂法求得的最小特征值为1.00。
数值分析课程第五版课后习题答案李庆扬等
数值分析课程第五版课后习题答案李庆扬等在学习数值分析这门课程的过程中，课后习题的练习与答案的参考对于我们深入理解和掌握知识点起着至关重要的作用。
李庆扬等编写的《数值分析》第五版教材，其课后习题涵盖了丰富的知识点和多种解题思路。
下面，我将为大家详细解析部分课后习题的答案。
首先，让我们来看一道关于插值法的习题。
题目是：给定函数值$f(0)＝0$，＄f(1)＝1$，＄f(2)＝4$，利用线性插值和抛物插值分别计算$f(15)＄的值。
对于线性插值，我们设直线方程为$L_1(x)＝ax ＋ b$。
将已知的两个点$（0,0)＄和$（1,1)＄代入，可得方程组：＄＼begin{cases}b ＝ 0 ＼＼ a ＋ b ＝ 1\end{cases}＄解得$a ＝ 1$，＄b ＝ 0$，所以$L_1(x) ＝ x$。
则$f(15) ＼approxL_1(15) ＝ 15$。
对于抛物插值，设抛物线方程为$L_2(x)＝ax^2 ＋ bx ＋ c$。
将三个点$（0,0)＄，＄（1,1)＄，＄（2,4)＄代入，得到方程组：＄＼begin{cases}c ＝ 0 ＼＼ a ＋ b ＋ c ＝ 1 ＼＼ 4a ＋ 2b ＋ c ＝4\end{cases}＄解这个方程组，可得$a ＝ 1$，＄b ＝ 0$，＄c ＝ 0$，所以$L_2(x) ＝ x^2$。
则$f(15) ＼approx L_2(15) ＝ 225$。
接下来是一道关于数值积分的题目。
求积分$＼int_{0}＾｛1} x^2 dx$的数值解，分别使用梯形公式和辛普森公式。
梯形公式为：＄T ＝＼frac{b a}｛2} ＼times f(a) ＋ f(b)＄，代入$a ＝ 0$，＄b ＝ 1$，＄f(x) ＝ x^2$，可得：＄T ＝＼frac{1 0}｛2} ＼times 0^2 ＋ 1^2 ＝ 05$辛普森公式为：＄S ＝＼frac{b a}｛6} ＼times f(a) ＋ 4f(＼frac{a ＋ b}｛2}）＋ f(b)＄，代入可得：＄S ＝＼frac{1 0}｛6} ＼times 0^2 ＋ 4 ＼times （＼frac{1}｛2}）＾2 ＋ 1^2 ＝＼frac{1}｛3}＄再看一道关于解线性方程组的习题。
数值分析课程第五版课后习题答案(李庆扬等)(OCR)
根是ｘ，，２…，ｘ－，且Ｖ。ｘ，ｘ…·，ｘ）＝Ｖ，Ｃｘ６，ｘ…·）（ｘ－ｘ）…（ｘ－ｘ）。
Ｖ，（ｘｏ，ｘ，…ｘ－ｘ）＝１１】 －ｘ，）用ａ－ｘ，）
［证明］由
可得求证。
＝Ｖ，（Ｃｘ８，ｘ，…，ｘＸ））１１（ｘ－ｘ）
２、当ｘ＝１－１，２时，ｆ（ｘ）＝０，－３．４，求ｆ（ｘ）的二次插值多项式。
Ｌ，（ｘ）＝ｙ％（（ｘｘ６－－ｘｘ，）（（ｘｘ－２ｘ－ｘ２２））
ｙ＝ｆ（ｘ）＝ｆ０．５）＝－０．６９３１４７，ｙ２＝ｆ（ｘ）＝ｆ（０．６）＝－０．５１０８２６，则
Ｌ２（ｘ）＝ｙ。 （ｘ－ｘ）（ｘ－ｘ２）
（ｘ６－ｘ）ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ２）
（ｘ－ｘｏ）（ｘ－ｘ） （ｘ２－ｘｏ）（ｘ２－ｘ）
＝－０．９１６２９１×．（０（．ｘ４－０－．０５．）５（）ｘ（－００．．６４）－０．６－．
３０—＋２—９．ｘ９５８３ｘ３１￣０２＇＝０．８３３６×１０４
１４、试用消元法解方程ｘ组１＋１０＂ｘ＝１００
ｘ＋ｘ２＝２
，假定只有三位数计算，问结果是否
可靠？
［解］精确解为ｘ１＝０１００－＊１ １０＂－２ ，当使用三位数运 算时，得到
ｘ ＝１，ｘ２＝１，结果可靠。
１５、已知三角形面积ｓ＝ｓ去＝ ａｂｓｉｎｃ，其中ｃ为弧度，０＜ｃ＜ 且测量ａ，ｂ，ｃ
位有效数字；ｘ＝５６．４３０有５位有效数字；ｘ＝７×１０有２位有效数字。 ４、利用公式（３．３）求下列各近似值的误差限，其中ｘ，ｘ；，ｘ，ｘ；均为第３题所给
的数。
（１）ｘ＋ｘ２＋ｘ：
ｅ（ｘ＋ｘ写＋ｘ）＝＞
［解］
Ｅ（ｘ）＝Ｅ（ｘ）＋Ｅ（ｘ）＋Ｅ（ｘ；）
３＋ｔｘ１０＝１．０５×１０３
（２）ｘｘｘ；
数值分析课程第五版课后习题答案李庆扬等
数值分析课程第五版课后习题答案李庆扬等数值分析作为一门重要的数学课程，对于许多理工科学生来说是必须掌握的知识。
李庆扬等编著的《数值分析》第五版教材备受青睐，而课后习题的答案则成为了同学们检验自己学习成果、加深对知识理解的重要参考。
在学习数值分析的过程中，课后习题起到了巩固和拓展知识的关键作用。
通过完成这些习题，我们能够更加深入地理解数值分析中的各种算法和概念，如插值法、数值积分、常微分方程数值解法等。
而准确的答案则能够帮助我们及时发现自己的错误和不足，从而有针对性地进行改进和提高。
以插值法这一章节的习题为例，我们可能会遇到要求用拉格朗日插值多项式、牛顿插值多项式等方法来构造插值函数，并计算给定节点处的函数值。
在解答这类问题时，需要我们熟练掌握插值公式的推导和计算过程，同时要注意误差的分析和控制。
答案中会详细展示每一步的计算过程，让我们能够清晰地看到如何从给定的节点数据得到最终的插值结果。
对于数值积分部分的习题，可能会涉及到梯形公式、辛普森公式等不同的数值积分方法。
在求解过程中，需要准确确定积分区间和节点，计算相应的系数，并最终得到积分的近似值。
答案会给出具体的计算步骤和结果，同时还会对不同方法的精度和误差进行比较和分析，帮助我们更好地理解各种数值积分方法的特点和适用范围。
常微分方程数值解法的习题则通常要求我们运用欧拉方法、改进的欧拉方法、龙格库塔方法等求解给定的初值问题。
这需要我们对这些方法的原理和公式有深入的理解，并能够正确地进行编程实现或手算求解。
答案中会详细讲解每一种方法的应用过程，以及如何根据给定的精度要求选择合适的解法。
在求解课后习题的过程中，我们不能仅仅满足于得到答案的结果，更要注重理解答案背后的思路和方法。
比如，在遇到错误答案时，要认真分析自己的解题过程，找出错误的原因，并通过与正确答案的对比，加深对知识点的理解。
同时，我们还可以尝试对答案进行拓展和延伸，思考如何将所学的知识应用到实际问题中，提高自己解决实际问题的能力。
李庆扬数值分析第五版习题答案清华大学出版社
李庆扬数值分析第五版习题答案清华大学出版社Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】第一章 绪论1．设0x >,x 的相对误差为δ，求ln x 的误差。
解：近似值*x 的相对误差为*****r e x xe 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．利用公式求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解：*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少解：球体体积为343V R π= 则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=故度量半径R 时允许的相对误差限为1(*)10.333r R ε=⨯≈ 6．设028Y =，按递推公式1n n Y Y -=（n=1,2,…） 计算到100Y27.982≈（5位有效数字），试问计算100Y 将有多大误差解：1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后，有1000100Y Y =-即1000Y Y =，若取27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
数值分析课程第五版课后习题答案(李庆扬等)1之欧阳引擎创编
第一章 绪论（12）欧阳引擎（2021.01.01）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 n x n 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 nk k kεεεε；（3）*4*2/x x 。
精品!数值分析 第五版课后习题完整答案(李庆扬等)
第一章 绪论（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章习题答案-李庆扬
第1章复习与思考题习题0，依据定义|X*X|/X*|X*-X|=X***X X：按定义求解E(lnX)=|lnX-LnX*|=ln(X/X*)=ln(1-)X|lnX ln*|/lnX*ln(1)/lnX*ln(1)/lnX*按泰勒展开求解，2(x*)=(x*)(x x*)"()(x x*)f f f1(lnX*)|lnX lnX*|lnX*(X X*)**r(lnX*)|lnX lnX*|/ln */ln *X X X X问题是解法1错了吗？ 很小时，ln(1-)=，求n X 的相对误差*0.02*X X ： 按照定义：{(10.02)X*}*(1.02)n n X (X )(X X*)/X* 1.021(10.02)1n n n n n多项式展开，有100.02n n i i：按照泰勒展开11(X X*)*(X X*)nX*(0.02X*)0.02X*)(XX*)/X*0.02*/*0.02n n n n n nn nn nnX nX X n问题是解法1错了吗？10.02n n i i 收敛于（0.02/(1-0.02)= 0.002004008016032064128256释？应该没有错，按照泰勒展开，相当于将误差限放大了。
3、下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指出它们是几位有效数字* 1.1021x ，0.031x ，385.6，56.430，7 1.0。
解：有4位有效数字；X2有1位有效数字；X3有3位有效数字；X4有1位有效数字**24x x ，*1x x 其中***123x x x ，，，公式（2.3）：1(A*)|()*|(X )n k k kfX 解：******124124()()()()0.510x x x x x x ************1232311321234313()()()()0.031385.60.510 1.1021385.60.510 1.10210.0310.510(0.59768212.48488 1.708255)100.214790815x x x x x x x x x x x x******2*2442244323335(/)1/()/()()1/56.4300.5100.031/(56.430)0.510(10.031/56.430)/56.4300.510=(0.99945064681906787169945064681907)0.5/56.430100.885610x x x x x x x5、计算求体积要使相对误差限为1%，问度量半径R 所允许的相对误差是多少？解：球的体积公式3V R23(V)3R (R)/R 3(R)/R 0.01r有(R)0.01/3R (R)(R)/1/300r R6、设Y0=28，按递推公式11783100n nY Y ，1,2,3.....n 计算到Y100，若取27.982（解：n n n (Y )(Y Y )，所以(Y )0n 。