东北大学数值分析课后习题的答案.pptVIP

习题1(第10页) 1-1.下列各数都是经过四舍五入得到的近似值,试分 别指出它们的绝对误差限,相对误差限和有效数字的位数 x1-5.420,x2=0.5420,x2=0.00542,x4=6000,x=0.6×105 解绝对误差限分别为:1-0.5×10-3,82=0.5×10- E30.5×105,E1=0.5,85-0.5×104 相对误差限分别为:ε1=0.5×103/5.420=0.00923%, C12-0.00923%,6:3-0.0923%,84-0.0083%,85=8.3% 有效数位分别为:4位,4位,3位,4位,1位 1-2.下列近似值的绝对误差限都是0.005,试问它们有 几位有效数字.a-1.00031,b=0.042,c=0.00032 解有效数位分别为:3位,1位,0位 1-3.为了使101/2的相对误差小于0.01%,试问应取几位 有效数字? 解因为101/2=3.162.=0.3162.×10,若具有n位有效 数字,则其绝对误差限为0.5×101mn,于是有 E10.5×101/3.162.0.5×101m/30.01% 因此只需n=5.即取1012=3.1623 1-4.求方程x2-56x+1=0的两个根,使它们至少具有四 位有效数字(783≈27982 解x1=28+27.982-55.982,x2=1/x1=0.017863 二习题2(第50页) 2-2(1).用列主元 Gauss消元法解方程组 4 10-70x2|=7 解 10-707 707 0-707 3264 0-0.166.1 5-156 156 02.552.5 10-707 0-707 02.552.5 +02.5525 0-0.166.1 006.26.2 回代得解 2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中 A=132 b=6 解 A=132→33,所以A=}1 122 yI 465 DI y 4 y y3 21 再解3|x 2-4.对矩阵A进行LDM分解和 Crout分解,其中 A=456 61515 解 A=456→43 61515 612 故得 Crout分解:A=43 612 IDM分解为:A=21 341 2-5.对矩阵A进行LDLT分解和G分解,并求解方程组 Ax=b,其中 1648 A=45 b=2 解 1648 A=45 3 8-4 412 故得GG7分解:A= LDLT分解为:A= 9 4 y 0.25 解12 y y 0.875 2-33人y3(3 V3 1.7083 0.25 0.5451 再1 2-3x,|=0.875得x2 1.2916 1.7083 0.5694 2-6(1).给定方程组 0-x+y x y a.用Cra皿er法则求其精确解.b.用 Gauss消元法和列 主元 Gauss消元法求解,并比较结果.(用两位浮点计算) 解a.x=-1/-0.99=1.010101,y=0.98/-0.99=0.989899 b.用 Gauss消元法 102x+y 102x+y=1 x+ 100y=-100 回代得解:y=1,x=0. 再用列主元 Gauss消元法 02x+y=1 x+y=2 回代得解:y=1,x=1 2-8.用追赶法求解方程组 4 14-1 xxxxx 000 14 解 4 14-1 3-56 y 6.6667 解 0,得y3|=1.7857 0 0.47847 200 53.718 27.051 6.6667 8.2052 再解 1.7857,得 5.7693 0.47847 14.872 53.718 53.718 2-10.证明下列不等式 (1)x-yx-z+1z-yll :( 2) x-Ilylsx-yll 证明(1)xyl=‖(xz)+(z-y)‖s‖xl‖+z-y (2)因为‖x|=‖(xy)+y|s‖xyl+lyl 所以‖划|ysxⅶ,同理可证ly| x- yll 于是有|y‖s|xyl