数值分析复习题答案.docx

《数值分析复习题答案.docx》由会员分享，可在线阅读，更多相关《数值分析复习题答案.docx（16页珍藏版）》请在课桌文档上搜索。

1、数值分析复习题一、填空Chapterl绪论近似数x*=0.4231关于真值x=0.4229有位有效数字.用IoOO.1近似真值IO(X)时，其有效数字有4位，准确值X*与其有t位有效数字的近似值=21O(1O)的绝对误差为x*-x10vz设三=2.40315是真值X=2.40194的近似值，那么(有3位有效数字。-J-104=-104设一近似数X*=25231具有5位有效数字，那么其相对误差限是2x24,其绝对误差-104限是2xl-7x=.1-=当X很大时，为防止损失有效数字，应该使x+l+VxChapter2插值方法s(x)=3x6+6-5x2+1j那么九一3,-2,-1,0,1,2,3=

2、3假设f(x)=24+-3,那么fl,2,3,4,5,6=。对f(x)=x3+3x2-x+5，差商f0,l,2,3,4二0。设F(X)=X6-3/+/-5,那么差商/04,2,3,4,5,6=1丫=3)的均差几。/，%=5,/x4,x0,x2=9fl43,x2=i4,f,x3,x2=8,.那么均差fx4,x2,x=9O(交换不变性)X112-32(x+1)(%-2)+(x+l)(x-1)设有数据那么其2次Lanmge插值多项式为23,2次拟合多项式为(最正确平方逼近可求)。？以n+1个整数点k(k=O2,n)为节点的Lagrange插值基函数为0)(k力Idk(X)=T=0,1,2,.,n),

3、那么k=oXo?(注：九一，那么有拉格朗日插值公式：y4()=EyJk(X)k=0,x = 0,1,2,，2； y=0,l,2x3-lOxlS(X)=4132-(x-l)3+a(x-l)2+b(x-l)+clx2假设2是三次样条函数，那么：a=_3_,b=_3_,c=O。三次样条函数S(X)满足：S(X)在区间a,b内二阶连续可导，S(xk)=yk(),k=0,l,2,n,且满足S(X)在每个子区间xk,xk+l上是不超过三次的多项式。(1.5x+lX=0,2过(0,1),(2,4),(3,1)点的分段线性插值函数P(x)J-3x+lx=2,3i. Xxxlx2x3x4ii. yyyly2y3

4、y4设有函数表如：W丫m0mlrn2m3m4,那么可利用分段三次Hermite插值，其插值多项式的次方为三次.？Chapter3函数的最正确平方逼近无Chapter4数值积分与数值微分牛顿柯特斯求积公式的系数和=。积分区间的长度(b-a)o(验证梯形、辛普森、科特斯公式满足)？r311If(x)dx=-f(一)+-f(l)12数值求积公式品)434的代数精度为：2次代数精度。(依次将函数IKK-代入验证是否满足，可得代数精度)rill3的代数精度为：3次代数精度。求积公式求积分 f(x)dx的近似值，其辛卜生公式为b - ar q/、+ /0) + 4/(-)Jo2Jo/(X心-2/(一)-/

5、(一)+2/(一)求积分I的近似值,其复化梯形公式叱9+25/(/)=,l+dr+/S)=设J。，那么用梯形公式得近似值为22Chapters线性方程组的直接解法能用高斯消元法求解Ar=S的充要条件是A的各阶顺序主子式不为零(P113),当满足条件7且3时(各阶顺序主子式不为零),A可作LU分解，当满足条件。3时(A为n阶对称正定矩阵)，必有分解式A=LL，其中L是对角元素为正的下三角阵。Chapter6线性方程组的迭代解法-548，那么皿=17 ,设A=112,那么阿=20。7设45 o510 JAII 产3-233-,那么:ML=-8_,|隗=_3IAe=_9IAcMLM=24。P(A)方

6、阵A的谱半径是指maxA. in I矩阵A的条件数是指。MA)=MMAll非奇异矩阵A的条件数Cond(A)=? ? , A是病态是指条件数数值很大。？ ？则条件数Cond8（八）=ChaPter8非线性方程的数值解法解方程f(x)=0的简单迭代法的迭代函数(x)满足在有根区间内l(x)LVL那么在有根区间内任意取一点作为初始值，迭代解都收敛。R=。一0U)-2利用二分法求f(X)=0在外切上根的近似值，误差限为2及设f(X)可微，那么求方程2=f(x)根的牛顿迭代格式为求后的近似值，其牛顿迭代格式为M5厂W3的近似值，其牛顿迭代格式是5%求解方程/*)二的 Newton迭代公式为f()割线公

7、式为Xr.1=xk(RJt-Xrl)k=1,2,3./UJ-Z(V1)序列KJn=O满足递推关系：丫=10丫-1，（11=1，2，一）,假设丫0有误差，这个计算过程不稳定。Chapter9常微分方程初值问题的数值解法微分方程数值解的几何意义是指用直线代替曲线。？求解常微分方程处值问题y = (x,y), ax+)”0-4242纥=！T41426max26012-124=4=，4=；77P(BJ=须11.u12，求得：21,计算：13I21所以，该迭代公式不收敛(即：发散)。5x1-1Ix2+x3=18-x1+X2+4巧=6二、对下述方程组I12xi+x2-直接应用高斯一塞德尔迭代法求解是否收敛

8、？如果不收敛试设法给出收敛的迭代公式，并简述理由。解：1、迭代公式的矩阵形式：“二&0+九其中：=-(DL)-,=-(DL)-B,P(BJ =2.22.228.6max-0.2-4.2-6.6.l3 ，* 俎 =O,L= -2.2-10.0379/,4=-2.2 +10.0379z4 笆，水w :，”舁：10.2762 1,所以该迭代公式不收敛(即：发散)。2、构造收敛的迭代公式？ ？5x 11/ + / = 18-xl+x2+ 4x3 = 612司+乙-七=9化为:12x1 + x2 - X3 = 9* 5x1 -1 Ix2 + X3 = 18-x, +z +M =6那么可得到新的:A为严格

9、对角占优矩阵，所以采用雅可比和高斯塞德尔迭代都收敛!,b 二三、0、31用迭代公式”=产+次从一,其中：伏=0,1,2,)。求解Ar二问取什么实数。可使迭代收敛，什么。可使迭代收敛最快。解:k将心叫)化为标准形式册”=(E+A*令B=(E+A),f=-h,可得：=+/由已经条件可得:l + 20O2aOl + 5OaO1 + 30解得.4 =l+5,% =l + 4,4 = l + P(B) = 根据迭代法收敛的充要条件：max. 1 lz3 ”可得关于的不等式:l + 5all + 4all6r 2一一Ca05- 02-2 a a(-,0)2八a (-,0)/1,所以在 5时，夕(切1,即迭

10、代收敛。2、求解可使迭代收敛最快:题三示意图2分别将AB)Tl+5,p(8)=l+4,p(B)=l+作出曲线图,如上图所示。在丁)的区间P(B) 内，max九1Z3 I的曲线为黑色粗线，那么夕mini)为折线的最低点(红点)，即为曲线夕(3) = l + 5同和z3) = l + 的交点，求得:1Ct =3,使得(3)最小。1 a 判断夕(8),越小收敛精度越高。当 3时,0(B) = Pmin(B),所以迭代收敛最快。-2x1 + Ix2 + 33 = 12 +L)Ta = ( +L)Tb其过程同下题六(2)!五、给定线性方程组(1)写出GaUSS-Seidel迭代格式；(2)分析该迭代格式

11、是否收敛。检验高斯-赛德尔迭代，”川=8产+工，其中：民=(。+L)W.=(O+L)-%其过程同下题六(2)!10.5。50.510.5六、给定线性方程组Ax=b,其中A=L551J,证明雅可比迭代法发散，而高斯一赛德尔迭代法收敛。证明：迭代公式的矩阵形式：二瓦+,分别检验2(3),进行敛散性判断。I、检验雅可比迭代，幻十(其中：B=-DU=-DbiO-0.5-0.5B=-0.5O-0.5-05-0.5O,求解得：4=05,4=5,P(B)=I所以雅可比迭代发散！2、检验高斯-赛德尔迭代，”=用/)+工，其中：s=(D+L),Uffs=(D+L)lb0-0.5-0.5Bx=00.25-0.25

12、00.1250.375J求解得.A=0=0.3125+0.1654z,=0.3125-0.1654z0(B)=O.35361,所以雅可比迭代收敛！20x(M)=(E-aA)x(k)-ab=x(k+i)=Bx(k+f,那么：B=(E-aA)J=abf求解B的特征值4,可求得P(B)只需证明夕(8)1即可证明收敛。过程同上！3213A=,b=八、给定线性方程组Ar=A,其中L12JL-1J,用迭代公式(=+S-Ah)(k=0,l,2)求解，问取什么实数可使迭代收敛，什么。可使迭代收敛最快。解：同题三！ChaPter8非线性方程的数值解法一、在求非线性f(x)=O根的近似值时，论证简单迭代法一般为线

13、性收敛，而牛顿迭代法为平方收敛。P200证明：？1、一般迭代法：Xk+产g(Xk)，k=0,12，由于/_/+=g(4)(9一五),g*,*,_畀=IgC)Ig()k8所以,U4)4,因而同假设Iga)I且Iga)IQ(X)=SinX1、证明：由条件，迭代函数为332I2(x)=-sin.d-lXeR可得：3I3,所以，对于任意的初始值不，该方法收敛。2、令：/U)=2cosx-3x+12,那么其导数f(x)=-2SinX-3f(x,)2cosxl-3xl.+12=K=4+二白一由牛顿迭代公式J(Z)可得：ZSInZ+J三、给定方程丁+4炉-10=分析该方程存在几个根，并构造求近似根的迭代公式

14、，证明所用的迭代公式是收敛的。解：=_81、令：/(x)=x3+4x2-10t令/()=3+8x=0,解得不(-OO,-,0/(X)在3为增函数，3为减函数，0，+0。)为减函数，具体函数形状如图a所示，又由于QQQ/(-2)=-0.515,/(0)=-10x1=-f(-)-0.51503,建立坐标系，从图中可以看出3,3为局部-8,0)最大值，/(0)=-100j图中可知该方程有一个根。图a图b2、由于/=一50,/(2)=140,/(x)=6+80构造牛顿迭代公式，f(xk)+4x-10xk+l=xk卜=Xk_7Ef(xky34+8%证明：验证迭代公式是否满足以下条件：？在xwl,2上，/

15、(幻、八X)存在且F(X)O(x)O符号不变，满足条件；/(2)0%,XC,由于/(X)O,应使f(o),比方XO=1.5,/(1.5)=2.375O,即可假设满足条件；综上，可知该迭代公式收敛！五、给定方程X-L优一2二。(1)分析该方程存在几个根，找出每个根所在的区间；(2)构造求近似根的迭代公式，并证明所用的迭代公式是收敛的。解：方法同上下题！四、给定方程/(幻=(%一1)-1=。分析该方程存在几个根；用迭代法求出这些根，精确至四位有效数；证明所试用的格式是收敛的。解：1、分析方程存在几个根：/(x)=(x-l)2x2Z0j所以/3在(，向上为增函数，同时y)o,所以f(X)存在一个根。

16、2、用迭代法求解：迭代格式为：加=e+1,由于/0,八1取Xo=I5,计算过程如下：Ixp-X11I-103=-1014*22,具有四位有效数字，所以X=12238。3、证明：迭代函数g(%)=Xz=+Ig(x)=-2x检验是否满足一下条件，g(x)=e+1在xl,2上存在且连续，满足条件；(2)g()在XWu2上为减函数,因而，可得：g(x)2,xU,2,满足条件：M(X)I=-2xe21=L02-ah2 + tz二、假设用梯形公式求)=N的近似解，其中,试证明:X=(1)(2)对固定的玉二，当7-时，”收敛于准确解。解：(1)根据梯形公式(改良的欧拉法公式)：由条件得：/，y)=y=y),)()二%h.rz、2-ah利用:那么:证毕!XfJ/(XD=X=W-2-ah.2-ah.1-.2-ah.i.2-ah.iX=二前I=R)NT=(否)X-o3五、初值问题1武)二的解为)一3。假设=九y是用改良欧拉公式得到的y()在巧y(%)-K=2,Z=1,2,3,处的近似值，证明：6。XO)=0h+=-lf(k，然)+,m),%=0,L2,一1证明：根据改良的欧拉法公式：2f(xty)=x2=f(x)fxi=ih9yif代入得到迭代过程，如下：必T)*)+4尸将，式相加可得：62那么:命题得证。