数值分析第五版第5章与第7章习题答案
数值分析第五章答案
数值分析第五章答案【篇一：数值分析第五版计算实习题】第二章2-1程序：clear;clc;x1=[0.2 0.4 0.6 0.8 1.0];y1=[0.98 0.92 0.81 0.64 0.38];n=length(y1);c=y1(:);or j=2:n %求差商for i=n:-1:jc(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1));endendsyms x df d;df(1)=1;d(1)=y1(1);for i=2:n %求牛顿差值多项式df(i)=df(i-1)*(x-x1(i-1));d(i)=c(i)*df(i);enddisp(4次牛顿插值多项式);p4=vpa(collect((sum(d))),5) %p4即为4次牛顿插值多项式,并保留小数点后5位数 pp=csape(x1,y1, variational);%调用三次样条函数 q=pp.coefs;disp(三次样条函数);for i=1:4s=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1];s=vpa(collect(s),5)endx2=0.2:0.08:1.08;dot=[1 2 11 12];figureezplot(p4,[0.2,1.08]);hold ony2=fnval(pp,x2);x=x2(dot);y3=eval(p4);y4=fnval(pp,x2(dot));plot(x2,y2,r,x2(dot),y3,b*,x2(dot),y4,co);title(4次牛顿插值及三次样条);结果如下：4次牛顿插值多项式p4 = - 0.52083*x^4 + 0.83333*x^3 - 1.1042*x^2 + 0.19167*x + 0.98 三次样条函数x∈[0.2,0.4]时， s = - 1.3393*x^3 + 0.80357*x^2 - 0.40714*x + 1.04 x∈[0.4,0.6]时，s = 0.44643*x^3 - 1.3393*x^2 + 0.45*x +0.92571 x∈[0.6,0.8]时，s = - 1.6964*x^3 + 2.5179*x^2 - 1.8643*x + 1.3886 x∈[0.8,1.0]时，s =2.5893*x^3 - 7.7679*x^2 + 6.3643*x - 0.80571 输出图如下2-3（1）程序：clear;clc;x1=[0 1 4 9 16 25 36 49 64];y1=[0 1 2 3 4 5 6 7 8];%插值点n=length(y1);a=ones(n,2);a(:,2)=-x1;c=1;for i=1:nc=conv(c,a(i,:));endq=zeros(n,n);r=zeros(n,n+1);for i=1:n[q(i,:),r(i,:)]=deconv(c,a(i,:));%wn+1/(x-xk)enddw=zeros(1,n);for i=1:ndw(i)=y1(i)/polyval(q(i,:),x1(i));%系数endp=dw*q;syms x l8;for i=1:nl8(i)=p(n-i+1)*x^(i-1);enddisp(8次拉格朗日插值);l8=vpa(collect((sum(l8))),5)xi=0:64;yi=polyval(p,xi);figureplot(xi,yi,x1,y1,r*);hold ontitle(8次拉格朗日插值);结果如下：8次拉格朗日插值l8 =- 3.2806e-10*x^8 + 6.7127e-8*x^7 - 5.4292e-6*x^6 +0.00022297*x^5 - 0.0049807*x^4 + 0.060429*x^3 - 0.38141*x^2 +1.3257*x输出图如下：第五章4-1（3）程序：clc;clear;y= @(x) sqrt(x).*log(x);a=0;b=1;tol=1e-4;p=quad(y,a,b,tol);fprintf(采用自适应辛普森积分结果为： %d \n, p);结果如下：采用自适应辛普森积分结果为： -4.439756e-01第九章9-1(a)程序：clc;clear;a=1;b=2;%定义域h=0.05;%步长n=(b-a)/h;y0=1;%初值f= @(x,y) 1/x^2-y/x;%微分函数xn=linspace(a,b,n+1);%将定义域分为n等份 yn=zeros(1,n);%结果矩阵yn(1)=y0;%赋初值%以下根据改进欧拉公式求解for i=1:nxn=xn(i);xnn=xn(i+1);yn=yn(i);yp=yn+h*f(xn,yn);yc=yn+h*f(xnn,yp);yn=(yp+yc)/2;yn(i+1)=yn;endxn=yn;%以下根据经典四阶r-k法公式求解for i=1:nxn=xn(i);yn=yn(i);k1=f(xn,yn);k2=f(xn+h/2,yn+h/2*k1);k3=f(xn+h/2,yn+h/2*k2);k4=f(xn+h,yn+h*k3);yn=yn+h/6*(k1+2*k2+2*k3+k4);yn(i+1)=yn;enddisp(改进欧拉法四阶经典r-k法); disp([xn yn])结果如下：改进欧拉法四阶经典r-k法 110.998870.998850.99577 0.99780.991140.996940.985320.996340.978570.996030.971110.996060.963110.996450.95470.997230.945980.998410.9370510.92798 1.0020.91883 1.00440.90964 1.00730.90045 1.01060.89129 1.01430.88218 1.01840.87315 1.02290.86421 1.02780.85538 1.03310.84665 1.0388（b）程序：clc;clear;a=0;b=1;%定义域h=[0.1 0.025 0.01];%步长y0=1/3;%初值f= @(x,y) -50*y+50*x^2+2*x;%微分函数 xi=linspace(a,b,11);y=1/3*exp(-50*xi)+xi.^2;%准确解 ym=zeros(1,11);for j=1:3【篇二：数值分析(第五版)计算实习题第五章作业】题：lu分解法：建立m文件function h1=zhijielu(a,b)%h1各阶主子式的行列式值[n n]=size(a);ra=rank(a);if ra~=ndisp(请注意：因为a的n阶行列式h1等于零，所以a不能进行lu 分解。
李庆扬_数值分析第五版第5章与第7章习题答案
所以2-范数为0.6074
F-范数0.8426
13、求证：
（a） ；
（b） 。
根据定义求证。
。
14、设 且非奇异，又设 为 上一向量范数，定义 。试证明 是 上向量的一种范数。
根据向量范数的定义来证明：
要求就有正定性，齐次性，三角不等式等性质。
显然 ， 、
，从而 是 上向量的一种范数。
（9）不动点迭代法 ，其中 ，若 则对任意处置x0迭代都收敛。（对）
（10）弦截法也是不动点迭代法的特例（正确）
习题
1、用二分法求方程 的正根，要求误差 。
[解]令 ，则 ， ，所以有根区间为 ；
又因为 ，所以有根区间为 ；
，所以有根区间为 ；
，所以有根区间为 ；
，所以有根区间为 ；
，所以有根区间为 ；
（1） 为对称正定矩阵；
（2）
，所以 为对称正定矩阵。
由于 为对称正定矩阵，所以
则
第7章
复习与思考题
1.什么是方程的有根区间？它与求根有何关系？
P213，若 且 ，根据连续函数性质可知 在 内至少有一个实根，这时称 为 的有根区间。
2.什么是二分法？用二分法求 的根， 要满足什么条件？
P213
一般地，对于函数 如果存在实数c,当x=c时，若 ,那么把x=c叫做函数 的零点。解方程即要求 的所有零点。
从而 ，
又当 时，
。
当 时，
。
综上所述， 时最小，这时 ，即 。
18、设 ，计算A的条件数
由 可知， ，从而ቤተ መጻሕፍቲ ባይዱ
，
由 ，
，
由 ，
可得 ，从而
。
数值分析课程第五版课后习题答案
=
1 = 1.7863 × 10 − 2 。 55.982
8、当 N 充分大时，怎样求 ∫ [解]因为 ∫
N +1 N
1 dx ？ 1+ x2
1 dx = arctan( N + 1) − arctan N ，当 N 充分大时为两个相近数相 1+ x2
减，设 α = arctan( N + 1) ， β = arctan N ，则 N + 1 = tan α ， N = tan β ，从而 tan(α − β ) = 因此 ∫
5、计算球体积要使相对误差限为 1%，问度量半径 R 允许的相对误差是多少？ 4 ε * ( π (R* )3 ) 4 3 [解]由 1% = ε r* ( π ( R * ) 3 ) = 可知， 4 3 * 3 π (R ) 3 ′ 4 4 4 ε * ( π ( R * ) 3 ) = 1% × π ( R * ) 3 = π ( R * ) 3 ε * ( R * ) = 4π ( R * ) 2 × ε * ( R * ) ， 3 3 3
ε * ( y n ) = 10ε * ( y n −1 ) = 10 n ε * ( y 0 ) ，
1 1 从而 ε * ( y10 ) = 1010 ε * ( y 0 ) = 1010 × × 10 − 2 = × 10 8 ，因此计算过程不稳定。 2 2 12、计算 f = ( 2 − 1) 6 ，取 2 ≈ 1.4 ，利用下列公式计算，哪一个得到的结果最 好？ 1 ( 2 + 1)
* r
x= x
*
ε ( x * ) = n( x * ) n −1 2% x * = 2n% ⋅ x * ，
(完整版)数值分析第五版答案(全)(最新整理)
第一章 绪论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．设的相对误差为2%，求的相对误差。
x 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 εε≈⋅ 且为2(*)r e x ((*))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．利用公式(2.3)求下列各近似值的误差限：(1) ,(2) ,(3) .***124x x x ++***123x x x **24/x x 其中均为第3题所给的数。
****1234,,,x x x x 解：*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 ππ===A A (*)(*)3(*)r p r r V C R R εεε∴≈=A 又%1(*)1r V ε=故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006．设，按递推公式 （n=1,2,…）028Y =1n n Y Y -=-计算到（5位有效数字），试问计算将有多大误差？100Y 27.982≈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 εεε-∴=+=⨯的误差限为。
数值分析课程第五版课后习题答案(李庆扬等)
数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章：数值分析导论1. 解答：数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题，其结果可能不是完全精确的，但是能够满足工程或科学应用的要求。
2. 解答：数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础，对许多领域都有着广泛的应用，如物理学、经济学、生物学等。
3. 解答：数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比，数值方法一般更加灵活和高效，可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答：计算误差是指数值计算结果与精确解之间的差异。
在数值计算中，由于计算机的有限精度以及数值计算方法本身的近似性等因素，都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章：数值误差分析1. 解答：绝对误差是指实际值与精确值之间的差异。
例如，对于一个计算出的数值近似解x和精确解x_0，其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度，通常被用作评估数值计算方法的好坏。
2. 解答：相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0，其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度，常用于评估数值计算方法的收敛速度。
3. 解答：舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数，因此在进行数值计算过程中，舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答：误差限度是指对于给定的数值计算问题，所能容忍的误差范围。
数值分析课后习题解答
课后习题解答第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式(1.2.2)(1.2.3)则得有5位有效数字，其误差限，相对误差限有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？（1）（2）解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。
（1）（2）4.近似数x*=0.0310,是 3 位有数数字。
5.计算取，利用：式计算误差最小。
四个选项：第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。
线性插值时，用0.5及0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值误差限，故2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h应取多少?解：用误差估计式（5.8），令因得3. 若，求和.解：由均差与导数关系于是4. 若互异，求的值，这里p≤n+1.解：，由均差对称性可知当有而当P＝n＋1时于是得5. 求证.解：解：只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式（5.17）得其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析第五版习题答案清华大学出版社
若，则
则
则法方程组为
从而解得
故
均方误差为
21。在某佛堂反应中，由实验得分解物浓度与时间关系如下：
时间
0 5 10 15 20 25 30 35 40 45 50 55
浓度
0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.62 4.64
用最小二乘法求。
计算到。若取（5位有效数字），试问计算将有多大误差？
解：
……
依次代入后，有
即，
若取,
的误差限为。
7．求方程的两个根，使它至少具有4位有效数字（）。
解：，
故方程的根应为
故
具有5位有效数字
具有5位有效数字
8．当N充分大时，怎样求？
解
设。
则
9．正方形的边长大约为了，应怎样测量才能使其面积误差不超过？
解：正方形的面积函数为
证明：
若
令，可得
当时，
当时，
又，故
得证。
10。证明切xx多项式满足微分方程
证明：
切xx多项式为
从而有
得证。
11。假设在上连续，求的零次最佳一致逼近多项式？
解：
在闭区间上连续
存在，使
取
则和是上的2个轮流为“正”、“负”的偏差点。
由切xx定理知
P为的零次最佳一致逼近多项式。
12。选取常数，使达到极小，又问这个解是否唯一？
解：
令
则在上为奇函数
又的最高次项系数为1，且为3次多项式。
与0的偏差最小。
从而有
13。求在上的最佳一次逼近多项式，并估计误差。
解：
于是得的最佳一次逼近多项式为
数值分析第五版第5章习题答案
第5章
）矩阵行列式的值很小。
）矩阵的范数小。
）矩阵的范数大。
（7）奇异矩阵的范数一定是零。
答：错误，
∞
•可以不为0。
（8）如果矩阵对称，则|| A||1 = || A||∞。
答：根据范数的定义，正确。
（9）如果线性方程组是良态的，则高斯消去法可以不选主元。
答：错误，不选主元时，可能除数为0。
（10）在求解非奇异性线性方程组时，即使系数矩阵病态，用列主元消去法产生的误差也很小。
答：错误。
对于病态方程组，选主元对误差的降低没有影响。
（11）|| A ||1 = || A T||∞。
答：根据范数的定义，正确。
（12）若A是n n的非奇异矩阵，则
)
(
cond
)
(
cond1-
=A
A。
答：正确。
A是n n的非奇异矩阵，则A存在逆矩阵。
根据条件数的定义有：
1
111111 cond()
cond()()
A A A
A A A A A A A
-
------
=•
=•=•=•
习题
如有侵权请联系告知删除，感谢你们的配合！。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第5章
）矩阵行列式的值很小。
）矩阵的范数小。
）矩阵的范数大。
（7）奇异矩阵的范数一定是零。
答：错误，
∞
•可以不为0。
（8）如果矩阵对称，则|| A||1 = || A||∞。
答：根据范数的定义，正确。
（9）如果线性方程组是良态的，则高斯消去法可以不选主元。
答：错误，不选主元时，可能除数为0。
（10）在求解非奇异性线性方程组时，即使系数矩阵病态，用列主元消去法产生的误差也很小。
答：错误。
对于病态方程组，选主元对误差的降低没有影响。
（11）|| A ||1 = || A T||∞。
答：根据范数的定义，正确。
（12）若A是n n的非奇异矩阵，则
)
(
cond
)
(
cond1-
=A
A。
答：正确。
A是n n的非奇异矩阵，则A存在逆矩阵。
根据条件数的定义有：
1
111111 cond()
cond()()
A A A
A A A A A A A
-
------
=•
=•=•=•
习题
第7章复习与思考题
习题
如有侵权请联系告知删除，感谢你们的配合！。