计算方法与实习第五版-习题答案.
计算机05级计算方法试卷A、B及参考答案

计 算 机 系 05 级计算方法试卷(A)(2008.1)班级_______________姓名______________学号_____________得分______________ 本卷考试时间为90分钟。
1.(10分)已知,0050.10010001≈具有七位有效数字,试从防止误差的角度给出适当的方法计算1000010001-,并说明理由。
2.(10分)简单叙述秦九韶法(或Horner 方法)计算多项式值的方法思想, 并应用该方法计算多项式322)(245+-+-=x x x x x f 在5.0=x 处的值,写出计算过程并说明计算过程中所需的乘法次数。
3. (10分)选用适当的方法求方程032=-x e x 在5.0 附近的一个根,要求所求根的误差不超过210-=ε。
4. (15分)用LU 分解法或高斯消去法解方程组5并利用插值多项式近似计算)(x f 在8.1处的值。
67. (15分)简述龙贝格(Romberg )求积方法的思想, 并选取适当的数值积分方法,求积分dx xe x ⎰102,要求误差不超过21021-⨯。
8. (15分)写出四阶龙格—库塔(Runge-Kutta )方法,并选用适当的方法求解初值问题,取2.0=h , 计算机05级计算方法试卷(B)(2008.1)班级_______________姓名______________学号_____________得分______________ 本卷考试时间为90分钟。
3.(10分)数列{n x }满足递推公式122,n n x x -=- ,若73.130≈=x (有三位有效数字),问 ①从0x 计算到k x 时误差有多大?②上述计算是稳定的?4.(10分)给出计算多项式334)(345+-+-=x x x x x f 在0x 处的值方法,使其所需乘法次数尽可能少.6. (10分)用适当数值方法求方程310x x +-= 在区间[0,1] 上的一个根,精度310-=ε。
最新数值分析(第五版)计算实习题第五章作业

数值分析第五章第一题:LU分解法:建立m文件function h1=zhijieLU(A,b)%h1各阶主子式的行列式值[n n]=size(A);RA=rank(A);if RA~=ndisp('请注意:因为A的n阶行列式h1等于零,所以A不能进行LU分解。
A的秩RA如下:')RA,h1=det(A);returnendif RA==nfor p=1:nh(p)=det(A(1:p,1:p));endh1=h(1:n);for i=1:nif h(1,i)==0disp('请注意:因为A的r阶主子式等于零,所以A不能进行LU分解。
A的秩RA和各阶顺序主子式h1依次如下:')h1;RAreturnendendif h(1,i)~=0disp('请注意:因为A的r阶主子式都不等于零,所以A能进行LU分解。
A的秩RA和各阶顺序主子式h1依次如下:')for j=1:nU(1,j)=A(1,j);endfor k=2:nfor i=2:nfor j=2:nL(1,1)=1;L(i,i)=1;if i>jL(1,1)=1;L(2,1)=A(2,1)/U(1,1);L(i,1)=A(i,1)/U(1,1);L(i,k)=(A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k);elseU(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j);endendendendh1;RA,U,L,X=inv(U)*inv(L)*bendend输入:>> A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2];>> b=[8;5.900001;5;1];>> h1=zhijieLU(A,b)输出:请注意:因为A的r阶主子式都不等于零,所以A能进行LU分解。
A的秩RA和各阶顺序主子式h1依次如下:RA =4U =10.0000 -7.0000 0 1.00000 2.1000 6.0000 2.30000 0 -2.1429 -4.23810 -0.0000 0 12.7333L =1.0000 0 0 0-0.3000 1.0000 0 00.5000 1.1905 1.0000 -0.00000.2000 1.1429 3.2000 1.0000X =-0.2749-1.32981.29691.4398h1 =10.0000 -0.0000 -150.0001 -762.0001列主元高斯消去法:建立m文件function [RA,RB,n,X]=liezhu(A,b)B=[A b];n=length(b);RA=rank(A);RB=rank(B);zhicha=RB-RA;if zhicha>0disp('请注意:因为RA~=RB,所以方程组无解')returnwarning off MATLAB:return_outside_of_loopendif RA==RBif RA==ndisp('请注意:因为RA=RB,所以方程组有唯一解')X=zeros(n,1);C=zeros(1,n+1);for p=1:n-1[Y,j]=max(abs(B(p:n,p)));C=B(p,:);B(p,:)=B(j+p-1,:);B(j+p-1,:)=C;for k=p+1:nm=B(k,p)/B(p,p);B(k,p:n+1)=B(k,p:n+1)-m*B(p,p:n+1);endendb=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n);for q=n-1:-1:1X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);endelsedisp('请注意:因为RA=RB<n,所以方程组有无穷多解') endend输入:>> A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2];>> b=[8;5.900001;5;1];>> [RA,RB,n,X]=liezhu(A,b),H=det(A)输出:请注意:因为RA=RB,所以方程组有唯一解RA =4RB =4n =4X =0.0000-1.00001.00001.0000H =-762.0001第二题:建立列主元高斯消去法m文件(题一中已有)(1)输入:>> format compact>> A=[3.01 6.03 1.99;1.27 4.16 -1.23;0.987 -4.81 9.34];>> b=[1;1;1];>> [RA,RB,n,X]=liezhu(A,b),h=det(A),C=cond(A)输出:请注意:因为RA=RB,所以方程组有唯一解RA =3RB =n =3X =1.0e+03 *1.5926-0.6319-0.4936h =-0.0305C =3.0697e+04(2)输入:>> A=[3.00 6.03 1.99;1.27 4.16 -1.23;0.990 -4.81 9.34]; >> b=[1;1;1];>> [RA,RB,n,X]=liezhu(A,b),h=det(A)输出:请注意:因为RA=RB,所以方程组有唯一解RA =3RB =3n =3X =119.5273-47.1426-36.8403h =-0.4070第三题:输入:>> clear>> A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10];>> b=[32 23 33 31]’;>> dA=det(A),lamda=eig(A),Ac2=cond(A,2)输出:dA =1.0000lamda =0.01020.84313.858130.28872.9841e+03下面分析误差性态:建立m文件:function Acp=pjwc(A,jA,b,jb,p)%Acp矩阵A的p条件数cond%pjwc:p范数解的误差性态分析%jA是A的近似矩阵jA=A+δA,jb=b+δbAcp=cond(A,p);dA=det(A);X=A\b;deltaA=jA-A;pndA=norm(deltaA,p);deltab=jb-b;pndb=norm(deltab,p);if pndb>0jX=A\jb;Pnb=norm(b,p);pnjx=norm(jX,p);deltaX=jX-X;pnjdX=norm(deltaX,p);jxX=pnjdX/pnjX;pnX=norm(X,p);xX=pnjdX/pnX;pndb=norm(deltab,p);xAb=pndb/pnb;pnbj=norm(jb,p);xAbj=pndb/pnbj;Xgxx=Acp*xAb;endif pndA>0jX=jA\b;deltaX=jX-X;pnX=norm(X,p);pnjdX=norm(deltaX,p);pnjX=norm(jX,p);jxX=pnjdX/pnjX;xX=pnjdX/pnX;pnjA=norm(jA,p);pnA=norm(A,p);pndA=norm(deltaA,p);xAbj=pndA/pnjA;xAb=pndA/pnA;Xgxx=Acp*xAb;endif (Acp>50)&(dA<0.1)disp('请注意:AX=b是病态的,A的p条件数Acp,A的行列式值dA,解X,近似解jX,解的相对误差xX,解的相对误差估计Xgxx,b或A的相对误差xAb依次如下:')Acp,dA,X',jX',xX',jxX',Xgxx',xAb',xAbj'elsedisp('请注意:AX=b是良态的,A的p条件数Acp,A的行列式值dA,解X,近似解jX,解的相对误差xX,解的相对误差估计Xgxx,b或A的相对误差xAb依次如下:')Acp,dA,X',jX',xX',jxX',Xgxx',xAb',xAbj'end输入:>> jA=[10 7 8.1 7.2;7.08 5.04 6 5;8 5.98 9.89 9;6.99 5 9 9.98];>> jb=b;p=2;>> Acp=pjwc(A,jA,b,jb,p)输出:请注意:AX=b是良态的,A的p条件数Acp,A的行列式值dA,解X,近似解jX,解的相对误差xX,解的相对误差估计Xgxx,b或A的相对误差xAb依次如下:Acp =2.9841e+03dA =1.0000ans =1.0000 1.0000 1.0000 1.0000 ans =-9.5863 18.3741 -3.2258 3.5240 xX =10.4661jxX =0.9842Xgxx =22.7396xAb =0.0076xAbj =0.0076Acp =2.9841e+03第四题:(1)输入:建立m文件:for n=2:6a=hilb(n);pnH(n-1)=cond(a,inf);endpnHn=2:6;plot(n,pnH);可见条件数随着n的增大而急剧增大(2)输入:>> n=2;H=hilb(n);>> x=(linspace(1,1,n))';>> b=H*x;>> [RA,RB,n,X]=gauss(H,b)输出:请注意:因为RA=RB,所以方程组有唯一解RA =2RB =2n =2X =1.00001.0000输入:>> r=b-H*X,deltax=X-x输出:r =deltax =1.0e-15 *0.4441-0.6661输入:>> n=3;H=hilb(n);>> x=(linspace(1,1,n))';>> b=H*x;>> [RA,RB,n,X]=gauss(H,b)>> r=b-H*X,deltax=X-x输出:X =1.00001.00001.0000r =1.0e-15 *0.2220deltax =1.0e-13 *-0.02000.1221-0.1255>> n=4;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.0000r =1.0e-15 *-0.4441-0.1110deltax =1.0e-12 *-0.02220.2485-0.59800.3886>> n=5;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00001.0000r =1.0e-15 *0.22200.1110deltax =1.0e-11 *-0.00350.0524-0.19370.2591-0.1148>> n=6;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00001.0000r =1.0e-15 *0.22200.1110deltax =1.0e-11 *-0.00350.0524-0.19370.2591-0.1148>> n=7;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00001.00001.0000r =1.0e-15 *0.22200.1110deltax =1.0e-09 *-0.00080.0219-0.14820.3854-0.42540.1677>> n=8;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00001.00001.00001.00001.0000r =1.0e-15 *-0.2220-0.1110-0.1110deltax =1.0e-06 *-0.00000.0018-0.02360.1279-0.34420.4870-0.34660.0978>> n=9;H=hilb(n);>> x=(linspace(1,1,n))'; >> b=H*x;>> [RA,RB,n,X]=gauss(H,b) >> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00001.00001.00001.00001.00001.0000r =1.0e-15 *0.44410.2220-0.22200.22200.2220-0.1110deltax =1.0e-04 *-0.00000.0002-0.00280.0197-0.07220.1471-0.16870.1017-0.0251>> n=10;H=hilb(n);>> x=(linspace(1,1,n))';>> b=H*x;>> [RA,RB,n,X]=gauss(H,b)>> r=b-H*X,deltax=X-xX =1.00001.00001.00001.00000.99991.00030.99961.00040.99981.0000r =1.0e-15 *0.4441-0.22200.2220-0.11100.11100.1110deltax =1.0e-03 *-0.00000.0001-0.00230.0205-0.09740.2669-0.43690.4214-0.22090.0485一、填空题(每空1分,共40分)1、中国共产党是_中国工人阶级___的先锋队,同时是中国人民和中华民族的先锋队,是_______中国特色社会主义事业___ 的领导核心。
计算方法与实习(第五版)期末复习资料

《计算机在材料科学中的应用》习题课第一章 误差等概念1. 误差来源:模型误差、观测误差、截断误差、舍入误差2. 绝对误差(限):e=x*-x ,|e|=|x*-x|≤ε3. 相对误差(限):e r =(x*-x)/x ,|e r |=|x*-x|/|x|≤εr4. 有效数字:|e|≤m-n 11025. 防止误差的危害:避免两相近数相减,多数作乘数或小数作除数,大数“吃”小数第二章 方程求根1. 根的存在及隔离2. 二分法:误差是()k+11b-a 23. 迭代法:'1x (x)|(x)|1 ||k k x x ϕϕε+=<-<, ,4. 加速法:'()L x ϕ≈取, 1111() L 1Lk k k k k k x x x x x x ϕ-+--+++⎧⎪⎨+-⎪⎩-==() 5. 牛顿迭代法:1000''1'111111'f()f()f ()0f ()f() f ()=c f()-f()f()()f ()=f()-f()f() f ()k k k k k k k k k k k k k k k k k k k k k k k x x x x x x x x x x x c x x x x x x x x x x x x x x x x λλ++--+--+->-----g ''=, 选取时使得简化牛顿法:,=拟牛顿法(割线法): ,=牛顿下山法:=, 选取下山因子使得1|f()|<|f()|k k x x +⎧⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎩第三章 方程组求解1. 消去法:高斯消去法,列主元消去法,高斯-约当法,消元因子 ()()k ikik k kka l a =消元公式 (k+1)(k)(k)ij ij ik kj (k+1)(k)(k)i i ik k a =a -l a (i,j=k+1,k+2,...,n)b =b -l b (i=k+1,k+2,...,n)⎧⎪⎨⎪⎩ 回代公式 kjn(k)(k)kjj=k+1k (k)kkb - a x x =(k=n,...,1)a∑2. 矩阵直接分解:紧凑格式3. 追赶法4. 迭代法:收敛条件1||||nii ij j j ia a =≠>∑①雅可比法迭代格式:ji n(k)i ij j=1j i(1)iib -a x x =(i=1,2,...,n) a k ≠+∑②高斯-赛德尔法迭代格式:jji i-1n(k+1)(k)i ij ij j=1j=i+1(1)iib -a x -a x x =(i=1,2,...,n)a k +∑∑第四章 插值法1. 插值多项式2012j j j j (1)n+1 ()()... , (x )= f( x )= y (j=0,1,...,n) x [a,b],() ()=()-()=()(n+1)!n n n n f x P x a a x a x a x P f R x f x P x x ξω+≈=++++=插值条件,插值节点,插值区间插值余项2. 拉格朗日插值: 插值基函数 n 001 () L ()()0 n nji j i i j i j j ix x i j l x x y i jx x ==≠-=⎧==⎨≠-⎩∑∏g ,3. 差商:10011002010122101k-2k 01k-2k-101k k k-1f(x )-f(x )f[x ,x ]=x -x f[x ,x ]-f[x ,x ]f[x ,x ,x ]=x -x f[x ,x ,...,x ,x ]-f[x ,x ,...,x ,x ]f[x ,x ,...,x ]=x -x 一阶差商二阶差商k 阶差商4. 牛顿插值公式f(x)=f(x 0)+f[x 0,x 1](x-x 0)+f[x 0,x 1,x 2](x-x 0)(x-x 1)+… +f[x 0,x 1,…,x n ](x-x 0)(x-x 1)…(x-x n-1) 5. 差分(等间距节点)111122111 = () , () -() -() - - k k k k k k k k k k k k k k m m m k k k x x kh x x f f x f x x h f f f f x x h f f f f x x h f f fm f f f δ+-+---+=+-∆≡∇≡≡∆=∆∆k 0k+1k 等距节点时,(k=0,1,...,n ),h=记则在处以为步长的向前差分:在处以为步长的向后差分:在处以为步长的中心差分:同样也有各自的阶差分111111122- -m m m k k k m m m k k k f f f f f fδδδ-----+-∇=∇∇=6. 牛顿前插公式20000001012nf f f ()=()+(-)()()....()...()()h 2!h n!h n n n f x f x x x x x x x x x x x R x -∆∆∆+--++--+7. 样条插值:三次样条插值,要求光滑、连续第五章 曲线拟合最小二乘原理2012n2i 01m j j j=1n (j=1,2,...,n),[]()...a (i=0,1,..., m),• (a ,...,a )= [P(x ) - y ] (x)(x,y ) m m p x a a x a x a x p n ϕ=++++∑j j 1n 有对数据(x ,y )在x ,x 上求一个m 次多项式适当选取使得,a 为最小值,则称为最小二乘拟合多项式是间的经验公式。
计算机网络第五版课后习题答案---精品管理资料

计算机网络第五版答案第一章概述1-01 计算机网络向用户可以提供那些服务?答: 连通性和共享1-02 简述分组交换的要点。
答:(1)报文分组,加首部(2)经路由器储存转发(3)在目的地合并1-03 试从多个方面比较电路交换、报文交换和分组交换的主要优缺点。
答:(1)电路交换:端对端通信质量因约定了通信资源获得可靠保障,对连续传送大量数据效率高。
(2)报文交换:无须预约传输带宽,动态逐段利用传输带宽对突发式数据通信效率高,通信迅速。
(3)分组交换:具有报文交换之高效、迅速的要点,且各分组小,路由灵活,网络生存性能好。
1—04 为什么说因特网是自印刷术以来人类通信方面最大的变革?答: 融合其他通信网络,在信息化过程中起核心作用,提供最好的连通性和信息共享,第一次提供了各种媒体形式的实时交互能力。
1-05 因特网的发展大致分为哪几个阶段?请指出这几个阶段的主要特点。
答:从单个网络APPANET向互联网发展;TCP/IP协议的初步成型建成三级结构的Internet;分为主干网、地区网和校园网;形成多层次ISP结构的Internet;ISP首次出现。
1-06 简述因特网标准制定的几个阶段?答:(1)因特网草案(Internet Draft)—-在这个阶段还不是RFC 文档。
(2)建议标准(Proposed Standard) ——从这个阶段开始就成为RFC 文档。
(3)草案标准(Draft Standard)(4)因特网标准(Internet Standard)1-07小写和大写开头的英文名字internet 和Internet在意思上有何重要区别?答:(1)internet(互联网或互连网):通用名词,它泛指由多个计算机网络互连而成的网络。
;协议无特指(2)Internet(因特网):专用名词,特指采用TCP/IP 协议的互联网络区别:后者实际上是前者的双向应用1-08 计算机网络都有哪些类别?各种类别的网络都有哪些特点?答:按范围:(1)广域网W AN:远程、高速、是Internet的核心网。
数值分析第五版第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 。
计算机网络第五版答案完整版

计算机网络第五版答案完整版Computers have become an integral part of our lives, and computer networking plays a crucial role in connecting various devices and enabling communication between them. In the fifth edition of the book "Computer Networks," the authors provide comprehensive answers and solutions to various exercises and questions, giving readers a complete understanding of the concepts and principles of computer networking.1. Introduction to Computer NetworksComputer networks are essential for transmitting data and enabling communication between devices. In this chapter, the authors explain the basics of networks, including the OSI model, network topology, and types of networks. They also provide answers to exercises that help readers grasp the fundamental concepts of computer networking.2. Physical LayerThe physical layer is responsible for transmitting raw bits over a communication channel. The authors cover topics such as analog and digital signals, transmission media, and modulation techniques. Through detailed explanations and examples, readers gain a solid understanding of the physical layer's functions and mechanisms.3. Data Link LayerThe data link layer ensures reliable data transfer between two connected nodes. This chapter focuses on topics like error detection and correction, flow control, and media access control. The authors provide accurateanswers to questions related to these concepts, enabling readers to comprehend the data link layer's role in establishing error-free communication.4. Network LayerThe network layer facilitates the delivery of data packets across multiple networks. Addressing, routing algorithms, and internet protocols are among the key topics discussed in this chapter. By examining the provided answers, readers can enhance their knowledge of the network layer's functionalities and protocols.5. Transport LayerThe transport layer provides end-to-end communication between applications running on different hosts. This chapter explores topics like multiplexing, demultiplexing, reliable data transfer protocols, and congestion control. The authors offer complete solutions to exercises, allowing readers to grasp the complexities and mechanisms of the transport layer.6. Application LayerThe application layer enables network applications to communicate with each other. This chapter covers topics such as domain name system (DNS), email protocols, and World Wide Web (WWW) protocols. The authors present accurate and detailed answers, enabling readers to understand the application layer's role in facilitating various network services.7. Network SecurityNetwork security is essential to protect data from unauthorized access and malicious activities. This chapter discusses topics like symmetric and asymmetric encryption, public key infrastructure (PKI), and network security protocols. The authors provide comprehensive answers, helping readers to comprehend the importance of network security and the techniques used to safeguard data.8. Multimedia NetworkingIn the modern era, multimedia applications require robust networking capabilities. This chapter delves into topics such as streaming and real-time applications, multimedia protocols, and quality of service (QoS) mechanisms. By examining the provided answers, readers can gain a deeper understanding of the challenges and solutions in multimedia networking.9. Network ManagementEfficient network management is crucial for ensuring the smooth operation of computer networks. This chapter covers topics like Simple Network Management Protocol (SNMP), network monitoring, and network troubleshooting. The authors provide accurate and comprehensive answers, allowing readers to learn about the tools and techniques used in network management.By providing complete and accurate answers to exercises, the fifth edition of "Computer Networks" equips readers with the knowledge and understanding required to master the field of computer networking. The authors' attention to detail, clear explanations, and concise yet informative solutions make this book an invaluable resource for students, professionals, and anyone interested in computer networks.。
数值分析课程第五版课后习题答案(李庆扬等)

数值分析课程第五版课后习题答案(李庆扬等)数值分析课程第五版课后习题答案(李庆扬等)第一章:数值分析导论1. 解答:数值分析是一门研究如何使用计算机来解决数学问题的学科。
它包括了从数学理论到计算实现的一系列技术。
数值分析的目标是通过近似的方式求解数学问题,其结果可能不是完全精确的,但是能够满足工程或科学应用的要求。
2. 解答:数值分析在实际应用中起着重要的作用。
它可以用于求解复杂的数学方程、计算机模拟及建模、数据的统计分析等等。
数值分析是科学计算和工程计算的基础,对许多领域都有着广泛的应用,如物理学、经济学、生物学等。
3. 解答:数值方法指的是使用数值计算的方式来求解数学问题。
与解析方法相比,数值方法一般更加灵活和高效,可以处理一些复杂的数学问题。
数值方法主要包括了数值逼近、插值、数值积分、数值微分、线性方程组的求解、非线性方程的求根等。
4. 解答:计算误差是指数值计算结果与精确解之间的差异。
在数值计算中,由于计算机的有限精度以及数值计算方法本身的近似性等因素,都会导致计算误差的产生。
计算误差可以分为截断误差和舍入误差两种。
第二章:数值误差分析1. 解答:绝对误差是指实际值与精确值之间的差异。
例如,对于一个计算出的数值近似解x和精确解x_0,其绝对误差为| x - x_0 |。
绝对误差可以衡量数值近似解的精确程度,通常被用作评估数值计算方法的好坏。
2. 解答:相对误差是指绝对误差与精确解之间的比值。
对于一个计算出的数值近似解x和精确解x_0,其相对误差为| (x - x_0) / x_0 |。
相对误差可以衡量数值近似解相对于精确解的精确度,常用于评估数值计算方法的收敛速度。
3. 解答:舍入误差是由于计算机的有限精度而引起的误差。
计算机中使用的浮点数系统只能表示有限的小数位数,因此在进行数值计算过程中,舍入误差不可避免地会产生。
舍入误差会导致计算结果与精确结果之间存在差异。
4. 解答:误差限度是指对于给定的数值计算问题,所能容忍的误差范围。
数值分析第五版课后答案2篇

数值分析第五版课后答案2篇数值分析第五版课后答案(一)第一章1.1 机器精度的数值为2^-52 ≈2.22 × 10^-16。
1.2 Example 1.2设f(x) = (1 - cosx)/sinx,则f(0)的分母为0,无法进行数值计算。
1.3 Example 1.3设f(x) = (1 - cosx)/sinx,则f(0)的分子为0,因此有f(0) = 0。
1.4 Example 1.4(a) 将x的值从1.8改为1.799,则f(x)的值由-0.000000000000159为0.000000000000313,差值为0.000000000000472。
(b) 我们有f'(x) = sinx/(1 - cosx) - 1/sin^2x。
将x的值从1.8改为1.799,利用f(x)和f'(x)的值可以得到下面的近似式:f(x + Δx) ≈ f(x) + f'(x)Δx = -0.000000000000159 + 0.449787416887455×0.001 = -0.000000000000137。
与(a)中的结果相近。
1.5 Example 1.5(a) 当x很接近于0时,函数值的符号取决于cosx的符号,其中cosx接近于1。
因此,函数值为正。
(b) 当x很接近于π时,函数值的大小趋于无穷大,因为分母趋向于0,而分子不为0。
1.6 Example 1.6(a) 因为函数在x = 0处是奇函数,所以它的导数为偶函数。
(b) 首先,我们有f''(0) = -2,因此x = 0是最大值。
其次,我们有f''(x) = -2 - 8sin^2x。
由于-f''(x)在x = 0处是正的,我们有当x越接近0时,f''(x)越小,也就意味着函数在x = 0处是严格的最大值。
1.7 Example 1.7(a) 我们有f(x) = x^3 - 2x^2 - 5x + 6,f'(x) =3x^2 - 4x - 5和f''(x) = 6x - 4。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.《计算方法》课程主要研究以计算 机为工具的 数值 分析方法 ,并评价 该算法的计算误差。 2.近似值作四则运算后的绝对误差限 公式为 ( x1 x2 ) ( x1 ) ( x2 ) ,近似值 1.0341的相对误差限不大于 1 10 2 , 则它至少有三位有效数字。 4
a. 0.41003×106 c. 4.10037×105 b. 0.41004×106 d. 上溢 a.0.00235 c.0.0023 b.0.0023471 d.0.00234711
第一章 绪论 练习
6. 自然数e*=2.718281828459045…,取 e≈2.71828,那么e的有效数字是: ( b )
解:1)求单调区间 f’(x)=-1-cosx,可知在(3.14, 0)区间f’(x)<0,单调递减 2)在(3.14, 0)区间逐步搜索 f(0)=1-0-sin0=1>0,f(1)=1-1-sin1=-sin1<0 ∴方程1-x-sinx=0在[0,1]中有且只有1个根。 ba 1 1 3)求二分次数 k 1 k 1 *10 3
绪论
习题1——10:设 f ( x) 8x5 0.4x4 4x3 9x 1 用秦九韶法求f(3)。 解:
8
0.4
24
4
0
9
1
x3
8
70.8
74.8
224.4
224.4
673.2 664.2
1992.6
1993.6
23.6
∴ f(3)=1993.6
第一章 绪论 练习
根的近似值
0.5000 0.7500 0.8750 0.9375 0.9063
序号
1 2 3 4 5
序号
6 7 8 9 10
根的近似值
0.9219 0.9141 0.9180 0.9200 0.9209
方程求根
习题2——3:用简单迭代法求方程ex-4x=0的 根,并验证收敛性,精确到4位有效数字。 解:1.找出方程的有根区间
(1)单调区间:
令f’(x)=ex-4=0, x=ln4≈1.4,所以有两个单调区间: [- ∞,1.4](递减)和[1.4, ∞](递增)
(2)有根区间:∴ 存在两个有根区间为:[0,1] 和[2,3]
[- ∞,1.4]区间:f(0)=1>0,f(1)=e-4<0,所以有根区间为:[0,1] [1.4,+ ∞]区间:f(2)=e2-8<0,f(3)=e3-12>0,所以有根区间为:[2,3]
方程求根
习题2——3:用简单迭代法求方程ex-4x=0的 根,并验证收敛性,精确到4位有效数字。
第一章 绪论 练习
3.设数据x1,x2的绝对误差限分别为0.05和 0.005,那么两数的乘积x1x2的绝对误差限 (x1x2)= 0.05 x2 0.005x1 。 4. 0.00234711 具有 5 位有效数字的近似值 是: ( b ) 5. 在β=10,t=5,-L=U=5的截断机上, 与数410037对应的规格化浮点数是: ( d )
3.2589
3.2590 4.3820 0.00078925
绪论
习题1——4:已知下列近似值x1=4.8675, x2=4.08675, x3=0.08675,求x1+x2+x3 的误差限。 4 5 5 e ( x ) 0 . 5 * 10 , e ( x ) 0 . 5 * 10 , e ( x ) 0 . 5 * 10 解: 1 2 3
2)判断二分次数 由(b-a)/2k+1=1/2k+1≤1/2*10-3,解得k≥3ln10/ln2≥9.965, 所以需要二分10次,才能满足精度要求。
方程求根
习题2——2:用二分法求方程2e-x-sinx=0在区 间[0,1]内的1个实根,要求3位有效数字。
解:3)迭代计算
∴x ≈0.921
3 ln( 10 ) k
ln( 2)
9.965
2
2
2
∴需二分10次
方程求根——二分法
习题2——2:用二分法求方程2e-x-sinx=0在区 间[0,1]内的1个实根,要求3位有效数字。
解:1)判断是否在该区间有且仅有一个根
f(0)=2>0,f(1)=2/e-sin1≈-0.1<0, f’(x)=-2e-x-cosx,f’=-3,-2/e-cos1<0
计算方法(数值分析)
习题答案——第一、二章
教师:马英杰 成都理工大学 核自学院
绪论
习题1——1:指出下列各数有几位有效数字
4.8675 4.08675 0.08675 96.4730 96*105 5 6 4 6 2
0.00096
2Байду номын сангаас
绪论
习题1——2:对下列各数写出具有5位有效数 字的近似值
3.25894 3.25896 4.382000 0.000789247
e( x1 x2 x3 ) e( x1 ) e( x2 ) e( x3 ) e( x1 x2 x3 ) e( x1 ) e( x2 ) e( x3 )
e( x1 ) e( x2 ) e( x3 )
0.5 *10 4 0.5 *105 0.5 *105 0.6 *10 4
绪论
习题1——6:一台10进制的计算机,4位字长, 阶码p∈[-2,3],可以表示的机器数有多少个? 给出它的最大数、最小数及距原点最近的非 零数,并求fl(x)的相对误差限。
解:β=10,t=4,L=-2,U=3
机器数个数:2*(β-1)*βt-1*(U-L+1)+1=2*9*103*6+1=108001 距原点最近的非零数:±0.1000*10-2 最大的数:0.9999*103 最小的数:-0.9999*103 相对误差限:0.5*10-3(舍入机), 10-3(截断机)
a.5位 b.6位 c.7位 d.8位
7. 数13.013627……的有四位有效数字 的近似值是: ( d )
a.13.00 c.13.014 b.13.02 d.13.013
方程求根
习题2——1:证明方程1-x-sinx=0在[0,1]中有 且只有1个根,用二分法求误差不大于1/2*10-3 的根需要迭代多少次?