数值计算答案北京大学

习题一解答22 3551.取3.14,3.15,—,空作为n的近似值，求各自的绝对误差，相对7 113误差和有效数字的位数。分析：求绝对误差的方法是按定义直接计算。求相对误差的一般方法是先求出绝对误差再按定义式计算。注意，不应先求相对误差再求绝对误差。有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。有了定理2后，可以根据定理2更规范地解答。根据定理2,首先要将数值转化为科学记数形式，然后解答。解：(1)绝对误差：e(x)=n-3.14=3.14159265 3.14=0.00159-^0.0016。相对误差：八x3.14有效数字：因为n=3.14159265…=0.314159265…X10,3.14=0.314X10,m=lo而n-3.14=3.14159265 3.14=0.00159-所以|jt-3.14|=0.00159-^0.005=0.5X10-2=-xlO-2=-xl0'_32 2所以，3.14作为n的近似值有3个有效数字。(2)绝对误差：e(x)=n-3.15=3.14159265 3.14=-0.008407—^-0.0085o相对误差：,、e(x)-0.0085“r,a-2er(x)= = =-0.27xlOx3.15有效数字：因为n=3.14159265…=0.314159265…XlO,3.15=0.315X10,m=lo而n-3.15=3.14159265 3.15=-0.008407-所以In-3.15|=0.008407 ^0.05=0.5X10-'=-xl0-1=-xl0'"22 2所以，3.15作为n的近似值有2个有效数字。(3)绝对误差：22e(x)=兀——=3.14159265••--3.142857143=-0.001264493•••=-0.0013相对误差:/、e(x)-0.0013 in.3er(x)='=———=-0.41x103x22T有效数字：因为n=3.14159265…=0.314159265-X10,22—=3.142857143=0.3142857143x10,m=l。722而万——=3.14159265…一3.142857143=-0.001264493…7所以2?n——=|3.14159265…-3.142857143|=0.001264493••<0.005=0.5xl0-2=-xl0-2=-xl01-32 2所以，三作为“的近似值有3个有效数字。7(4)绝对误差：355e(x)=% =3.14159265••--3,14159292=-0.0000002705•••=-0.000000271113相对误差：^(%)=£00=-0.000000271_0863)<10,7x355U3有效数字：因为n=3.14159265…=0.314159265-X10,355—=3.14159292=0.314159292x10，m=l0113355而% =3.14159265…一3.14159292=-0.0000002705…113所以355ZT--=|3.14159265----3.14159292|=0.0000002705•<0.0000005=0.5x10^=-xl0-6=-xl01-72 2所以，垩作为n的近似值有7个有效数字。113指出：①实际上，本题所求得只能是绝对误差限和相对误差限，而不是绝对误差和相对误差。2、用四舍五入原则写出下列各数的具有五位有效数字的近似数。346.7854,7.000009,0.0001324580,0.600300

解：346.7854-346.79,7.000009弋7.0000,0.0001324580^0.00013246,0.600300-0.60030o指出：注意0o只要求写出不要求变形。3、下列各数都是对准确数进行四舍五入后得到的近似数，试分别指出他们的绝对误差限和相对误差限和有效数字的位数。演=0.0315,x2=0.3015,x3=31.50,x4=5000。分析：首先，本题的准确数未知，因此绝对误差限根据四舍五入规则确定。其次，应当先求绝对误差限，再求相对误差限，最后确定有效数字个数。有效数字由定义可以直接得出。解：由四舍五入的概念，上述各数的绝对误差限分别是f（x,）=0.00005,£（x2）=0.00005,£（x3）=0.005,6*（x4）=0.5由绝对误差和相对误差的关系，相对误差限分别是£区）0.00005d（x）= -= 。0.16%,% 0.0315次々）==0.02%,£（x2）_0.00005x2次々）==0.02%,0.55000=0.01%.有效数字分别有3位、4位、4位、4位。指出：本题显然是直接指出有效数位、直接写出绝对误差，用定义求出相对误差。4.计算而的近似值，使其相对误差不超过0.1%。解：设取n个有效数字可使相对误差小于0.1%,则—xlO'"n<0.1%,2a1而3wji544,显然％=3,此时，—xl01-n=—xl01-n<0.1%,2a} 2x3即LxlOi<10-3,6

也即6x10">1(/所以，n=4o此时，V10=3.162o5、在计算机数系F(10,4,-77,77)中，对x,=0.14281x103与々=-0.314159x10、试求它们的机器浮点数1(xj(i=l,2)及其相对误差。解：y7(x1)=0.1428xl03,e(/Z(x1))=x1-/Z(xl)=0.14281xl03-0.1428xl03=0.00001xl03,y7(x2)=-0.3142xlOl,e(77(x2))=x2-/7(x2)=-0.314159x10'-(-0.3142x10')=0.00041x10'其相对误差分别是_0.00001X103-_0.00001X103-0.1428X103=0.007%,e2_0.000041x1()1--0.3142x10'=-0.013%。6、在机器数系F(10,8,L,U)中，取三个数x=0.23371258x10^,y=0.33678429xl02,z=-0.33677811X102,试按(x+y)+z,x+(y+z)两种算法计算x+y+z的值，并将结果与精确结果比较。解：fl((x+y)+z)=(0.23371258x10"1+0.33678429xlO2)-0.336778HxlO2=(0.00000023x102+0.33678429X102)一0.336778HxlO2=0.33678452xl02-0.3367781IxlO2=0.0000064lx102y7(x+(y+z))=0.23371258xl0-4+(0.33678429xl02-0.3367781IxlO2)=0.23371258x10-4+0.00000618x102=0.00000023xlO2+0.00000618xlO2=0.0000064lx102精确计算得：x+y+z=0.23371258X10-4+0.33678429X102-0.3367781IxlO2=(0.00000023371258xlO2+0.33678429xlO2)-0.336778HxlO2=0.33678452371258X102-0.3367781IxlO2=0.0000641371258xl02第一种算法按从小到大计算，但出现了两个数量级相差较大的数相加,容易出现大数吃小数.而第二种算法则出现了两个相近的数相减，容易导致有效数位的减少。计算结果证明，两者精度水平是相同的。***在机器数系F(10,8,L,5中，取三个数x=0.23371258xl0"',y=0,33678429xl0-2,z=-0.33677811xl02,试按(x+y)+z,x+(y+z)两种算法计算x+y+z的值，并将结果与精确结果比较。解：力((x+y)+z)=(0.23371258x10^+0.33678429xlO-2)-0.336778HxlO2=(0.00233713X10-2+0.33678429x10-2)-0.3367781IxlO2=0.33912142X10-2-0.3367781IxlO2=0.00003391X102-0.3367781IxlO2=-0.3367442x102fl(x+(y+z))=0.23371258x10^4-(0.33678429xlO-2-0.336778HxlO2)=0.23371258x10-4+(0.00003368xl02-0.336778HxlO2)=0.23371258x10^-0.33674742xl02=0.00000023xlO2-0.33674742x102=-0.33674719xl02第一种算法是按从小到大的顺序计算的，防止了大数吃小数，计算更精确。精确计算得：x+y+z=0.23371258x10-4+o33678429x10-2-0.3367781IxlO2=0.000023371258+0.0033678429-33.677811=0.003391214158-33.677811=-33.674419785842=-0.33674419785842xl02显然，也是第一种算法求出的结果和精确结果更接近。7、某计算机的机器数系为F(10,2,L,U),用浮点运算分别从左到右计算及从右到左计算1+0.4+0.3+0.2+0.04+0.03+0.02+0.01试比较所得结果。解：从左到右计算得1+0.4+0.3+0.2+0.04+0.03+0.02+0.01=0.1x10+0.04x10+0.03x10+0.02x10+0.00x10+0.00x10+0.00x10+0.00x10=0.19x10=1.9从右到左计算得

1+0.4+0.3+0.2+0.04+0.03+0.02+0.01=0.01+0.02+0.03+0.04+0.2+0.3+0.4+1=0.1x101+0.2X10T+0.3x10-+0.4xl0-1+0.2+0.3+0.4+l=0.1+0.2+0.3+0.4+1=0.1x10+1=0.1x10+0.1x10=0.2x10=2从右到左计算避免了大数吃小数，比从左到右计算精确。8、对于有效数％=-3.105,/=0.001,9=0.100,估计下列算式的相对误差限X=$+x2+x3,y2=X1X2x3,y3=口X3分析：求和差的相对误差限采取先求出和差的绝对误差限再求相对误差限的方法。求积商的相对误差限采取先求每一个数的相对误差限再求和的方法。解：因为百=一3.105,9=0.001,覆=。1°°都是有效数，所以£（X）=0.0005,£（x2）=0.0005,£（x3）=0.0005g）=0.00053.1050.00050.001=50%,3（项）=0.00050.100=0.5%则E（xl+x2+x3）=£（百）+£（々）+£（&）=0.0005+0.0005+0.0005=0.0015C/ 、 +X2C/ 、 +X2+x3）+X2+X3）=-：! = y-|X1+x2+x310.0015 0.0015|-3.1O5+O.OO1+O.1OO|-3.004=4.99xl（）T=0.05%演尤也演）=6（X）+演马）+仪天）=0.16%+50%+0.5%=50.66%r/这）=3（X2）+S（x3）=50%+0.5%=50.5%当指出：如果简单地用有效数字与误差的关系计算，则不够精确。注意是相对误差限的讨论。符号要正确，商的误差限是误差限的和而不是差。9、试改变下列表达式，使其计算结果比较精确（其中国口1表示x充分接近0,|x|Q1表示x充分大）。(1)InXj-InXj,^=x2；/\1一COSX nIIr-ii,xw0且x口1;X—cotx9x^0且国口1o分析：根据算法设计的原则进行变形即可。当没有简单有效的方法时就采用泰勒展开的方法。解：(1)In玉一InZ=In土；x21—X1+x—(1—x)21-x1+x(l-x)(l+x)_l+x-(l-2x+x2)_3x-x2一~(l-x)(l+x)~-(l-x)(l+x)或(2〃)!+…)x2! +(-l)n+14!(5)1——cotx(2〃)!+…)x2! +(-l)n+14!(5)1——cotxX—< x X3 453B“45,2"R―…)(2〃)！指出：①采用等价无穷小代换的方法一般不可行。近似计算中的误差并不是无穷小量，利用无穷小量等价代换，两个量的差别可能恰恰是影响精度的因素。采用等价无穷小代换，可能只会得到精度水平比较低的结论。例如l-cosx2sin2l(4) 41—(1 (4) 41—(1 1 *,+(-1)/，1-cosx2!4!x xr2 / 2n +(一I)"1—_十.・・2! 4! (2〃)!-2

1cosxsinx—xcosx-cotx=-e .x xsinxxsinxx-xcosx,l . 、工. (xU1,sinx«x)xsinx_1-cosxsinx1-1= (xUl,cosx«l)sinx=0试与上例比较。有时候这种方法可以使用，例如因为cos(x+b)=cosxcosb-sinxsinb,当冏口1时,cosb=l,sinb=0cos(x+^)=cosxcosJ-sinxsin8=cosx-sinxQJ在这个计算中，由于X是常数，X的函数值实际上放大了每一项的计算结果，使得相近的数相减的问题不很突出。而利用一阶的泰勒展开/(x+<5)=/(x)+S[C)(x<J<x+S)，当同口1时，就有f(x+»=f(x)+3f\x),因此cos(x+<5)=cosx-Ssinx和上面的结果一样。但显然，用泰勒展开的方法具有一般性并能得到精度更高的结果，而且不会有方法上出错的可能。②采用洛必达法则也是不可以的。实际上，无论是等价无穷小还是洛必达法则都是极限方法，而因为近似计算中的误差虽然可以近似地看作是微分，但本质上却是一个确定的可能极小的小数而不是无穷小（趋于零的变量），因此近似计算是不能采用极限方法的。③转化的结果要化简，比如化繁分式为简分式，但不能取极限。取极限就违背的了数值计算的本意。所以，1 1-x 1 1-0,,1—x1+x1-1—x1+x1-01+0是错误的。④极小的数做除数，实际上是9型的不定型，要转化为非不定型。10、用4位三角函数表，怎样算才能保证1-COS2。有较高的精度？解：根据l-cos2。=2sin?1。，先查表求出sinl。再计算出要求的结果精度较高。指出：用度数就可以。不必化为弧度。11、利用7^=27.982求方程/-56%+1=0的两个根，使它们至少具有4位有效数字。解：由方程的求根公式，本方程的根为56±V^.56±2R=28土历— 2 2因为=27.982,则=28+^783=28+27.982=55.982如果直接根据求根公式计算第二个根，则因为两个相近的数相减会造成有效数字的减少，误差增大。因此根据韦达定理工/2=1，在求出％=55.982后这样计算％：x2=—=―5—=0.01786=0.1786x10'-%55.982这样就保证了求出的根有四位有效数字。12、试给出一种计算积分0近似值的稳定算法。1解：当n=0时，/0=e~}^x)exdx— —=l-e~ i(j i(jexdx=^v|=e-l)o0 0b b对L运用分部积分法(Judv=犹-J"u)得a a1 1 1In=e~}Jxnexdx=e"(x〃e1—nJxn~lexcbc)= xn~]exdx)0 0 01=1-ne-1jxn]exdx=1-nln_yo由此得到带初值的递推关系式[0=1—eT=1-〃/“T(〃=1,2,3,...)由递推公式i„=i-mn_1解得小=%_/“)，这是逆向的递推公式，对nL的值作估计，有e~ljxnexdx<e~le]^xncbc=o o另有=/Jxnexdx>e~'^xndx=e~]0 0(取e的指数为最小值0,将e'取作e°=1作为常数即可简化公式)则"七”"不那么，我们可以取其上下限的平均值作为其近似值。即取2/1+1可以看出，n越大，这个近似值越精确地接近于准确值。(n越大，L的上限和下限就越接近，近似值区间的长度就越短，近似值和精确值就越接近)此时，en-l=L-i*—In-1=——(In*—L)=—a,|Go|=--|©n|»计算是稳n n 〃！定的。实际上，如果我们要求19,可以先求出1卯，这样求出的19的误差是比120的误差小得多的，而L的误差本身也并不大。实际上，这样求出的L比直接计算出来的精确得多。补充题(一)1、给出数系F(10,4,-5,5)中的最大数、最小数和最小整数。解：最大数：0.9999X105；最小数:-0.9999X105；最小正数：0.0001X10-5o2、已知e=2.71828182845904523536028747…，求它在F(10,5,-5,5)和F(10,8,-5,5)中的浮点数。解:在F(10,5,-5,5)中，/7(e)=O.27183xlO'在F(10,8,-5,5)中，/7(e)=0.27182818x10'3、已知数e的以下几个近似数，它们分别有几位有效数字？相对误差是多少？%=2.7182,x=2.7183,%=2.7182818分析：题目没有说明近似数是通过哪种途径取得的，也就没有明确每个近似数和准确数之间的误差关系。所以，本题的解答应当从求近似数的误差开始。解：因为|e-xo|=O.OOOO8181<yXlO_3=^xl0^4,=0.00002<1xl0-4=^xl01-5, ,|e-x,|=0.00000003<1x10-7=1x101-8.所以，-=2.7182,3=2.7183,飞=2.7182818分别有4、5、8个有效数字。其相对误差分别是L_x||xl0--41er(x0)=। <2 <-xlO-3,’0闾2.71824|e-x1|<-xlO'4,|e-x2|<；x]0-74、数生留与下述各式在实数的意义上是相等的，(3+两3(1)(17-6^)3,(2)[(17-6V8)3]-),(3)(3-V8)6,(4)[(3+78)6]-1,(5)19601-6920的，(6)(I9601+6920指尸。试说明在浮点数系F(10,4,-8,8)中，用哪个公式计算出的结果误差最小。分析：本题实际上是一个算法分析与设计问题，也就是说要应用算法设计的基本原则进行分析讨论。解：在本例中，显然3和唬在浮点数系中是相近的数。进一步地，17和6小、19601和6920次也是相近的数。因此：①为避免相近的数相减，不应采用(1)、(3)、(5)三种计算方法。②在余下的三种计算方法中，(2)需要进行4次乘除法，(4)需要进行7次乘除法，(6)需要进行1次除法。从减少运算次数来说，应采用(6)„所以，采用(6)计算，计算结果误差最小。

X5、/(幻=m/+31-幻］/丁，当国口1时，如何计算才能获得准确的结果？解：当国口1(即很小时)，f(x)的分子是两个相近的小数相减，而分母也是一个小数，因此应避免简单地按原计算顺序直接计算，而应进行变形。由泰勒展开得2 2! 3! 4!2 2! 3! 4!TOC\o"1-5"\h\zr2 r3 xnln(l-x)=-x 2 3 n1

16x24--)x1

16x24--)x5+---]/x3f(X)=[( )/+( )*4+(.83 484397 x397 x1920 x 2448-此处最后略去部分的第一项为-2=-丝/120x326 3840当W口1时，这一部分是相当小的值，可以略去。指出：如果要提高计算精度，就可以考虑保留更多的项。补充题(二)

(―)1、计算e的近似值，使其误差不超过106O2、利用1 rM+,f(x)=-^—=l+x+x2+•••+/'+--——-(0<e<l,x<l)1-x (\-0x)n+2计算H0.1)的近似值，其误差不超过10-2,求no3、3.142和3.141分别作为71的近似数，各有几位有效数字？4、已知近似数x的相对误差限为0.3%,问x至少有几个有效数字?5、已知x的下列3个近似数的绝对误差限都是0.005,问它们的有效数字各有几位？a=138.00,b=-0.0132,c=-0.86x1046、设近似值x=1.234,且绝对误差界为0.0005,则它至少有几位有效数字？7、某校有学生6281人，通常说有6000人。下面哪个式子表示6000这个近似数合适？0.6X1040.60x1040.600X1040.6000x104分析与解答1、解：令f(x尸e;而i)(x尸ex,f(k)(o尸e0=L由麦克劳林公式，可知„2 ti「，xxeC=1+Xd F…-I 1 Xn\(n+1)!(0<8<l)2!n\(n+1)!n\(n+1)!(0<8<l)当x=l时，e=1+1H 1-•,2!J3故此⑴=——<一一。n(n+1)!(n+1)!当n=9时，Rn(l)＜10Y，符合要求。此时，e-2.718285解决这类问题其实很简单。只要知道了泰勒展开式，余下的就只是简单的计算了。泰勒(Taylor)中值定理：若函数f(x)在[a,b]上存在直至n阶的连续导函数，在(a,b)上存在n+1阶导函数，则对任意给定的x,x0G[a,,b],至少存在一点审(a,,b)，使得f(x)=f(X0)+f(xo)(X_/)+于3「，(X-Xo)2+•••+—~~(x-xo)"+- (x-%)，，+l2! 〃！ (n+1)!其中，&a)=▲—向叫做拉格朗日型余项。(n+1)!当xo=O时彳导到麦克劳林公式。/(x)=/(0)+fXoyJx+ +…+—瓜"+.["旬(%)及"+1 (0＜。＜I)2! n\ (n+1)!2、解:xn+]。一6x)"+2xn+]。一6x)"+20.1n+l

(l-0.1)n+20『二(1产0.9"2—9xl0<10-2,9-S+2)<|q-39"+2>[03所以,n=2o3、n=3.14159265...=0.314159265...x10,3.142=0.3142x10,m=1o因为n-3.142=3.14159265...-3.142=-0.00040...

所以，|n-3.142|=0.00040...40.0005=0.5x10-3=-xlO-3=-xlO1-42 2所以，3.142作为ti的近似值有4个有效数字。兀=3.1415926…，3.1415926••--3.1411=0.00059<0.005=0.5x10-2=ixlO-2=-xlO1-31 2 2小数点后几个0,10的指数的绝对值就是几。4、解：设x有n位有效数字，其第一位有效数字按最不利情况3=0.3%3

looo1

2(9+1)xlOl-n1

2x10xlO1-"=-xl0-n=

212xlOM由上可得6x10"=1000,n=2.2,所以取n=2o5、解I|x—iz|=|x—/?|=|x—c|0.005=-x10',所以mn=-2oa=138.00=0.13800x103,则m=3,所以n=3-(-2)=5,即a有5位有效数字;b=-0.0132=-0.132x1O1,贝I」nr-1,所以n=-1-(-2)=1,所以b有1位有效数字。c—0.86x104,则nr-4,所以n=-4*(-2)=-2<0,所以c没有有效数字。6、解：因为近似数x=1.234的绝对误差界为0.0005,所以k*一本0.0005=3x10-3,则mn-3O而x=1.234=0.1234x1Q1,则m=1,所以n=1-(-3)=4,所以，x=1.234有4位有效数字。7、解：哪个式子表示6000这个近似数合适实际上要看近似数6000有多少个有效数字。6281近似到十位、百位，千位分别是

6281=62806281=63006281=6000写成科学记数的形式分别是6281=6280=0.628x1046281=6300=0.63x1046281=6000=0.6X104可见，上述写法中，第一种是合适的。实际上，6281=0.6281x1O”,6000=0.6000x104所以nr4,而|6281-6000|=281=O.281xlO3而|6281-6000|=281=2所以mn=3,则n=m3=4-3=1,即近似数6000只有一个有效数字，所以，只有0.6x104这种写法是合适的。（二）1、已知测量某长方形场地的长为a=110米，宽为b=80米。若|a*-a|WO.1（米），|b*-b|WO.1（米），试求其面积的绝对误差限和相对误差限。2、已知三角形的两个内角的测量误差都不超过0.1。，则计算第三个角时，绝对误差不超过多少。3、若xi=1.03±0.01,X2=0.45±0.01,计算y= *的近似值并估计误差。24、已知测量某长方形场地的长为a=110米，宽为b=80米。若|a*-a|<0.2（米），|b*-b|WO.1（米），试利用多元函数的误差分析方法求其面积S=ab的绝对误差限和相对误差限，并与四则运算的误差分析比较。5、如果用电表测得一个电阻两端的电压和通过的电流分别是V=110±2(V),I=20±0.5(A)试运用欧姆定律/?=求这个电阻值试运用欧姆定律/?=求这个电阻值R的近似值，并估计所求出的近似值的绝对误差和相对误差。6、已知近似值a[=2.21,a2=4.63,a3=7.98是由四舍五入得到的，它们的绝对误差界都是0.005试估计始+%和q-4的相对误差界。。2

分析与解答

1、S=ab,£（S）=£（ab）<a£（b）+be（a）=\9A（m2）S（S）=5（ab）=丝2=19P(2)=-9=0.00217=0.217%P(2)=-92、提示：内角和为180°,而且180是准确数，没有误差。3、由已知，xi=1.03,Axj=0.01,X2=0.45,Ax2=0.01©所以,£a，X2）=2X1=2.06,北=0.7842,€（X1）=Ax]=0.01,8（X2）=Ax2=O.Olo所以，丫的绝对误差限为£（y）K<；（百/2）£（%）+力;（百,9）£（X2）=2.06x0.01+0.7842x0.01=0.028将有关数据代入函数表达式，可以求出函数值的近似值为y=X；=1.845,则y的相对误差限为、£（y）0.028..0/J（y）=—^.= =1.5%y1.845进一步地，本题的绝对误差限可以看作是0.05，那么计算结果中只需要保留到百分位就可以了，即最终结果取1.8，那么计算过程中各数只需要取到千分位。）4、（o6、略解。/（尤1，%2，七）=丑a+七 %f（xl,x2,x3）-f（a1,a2,a3）=（^-+x3）-（^-+a3）x3 a3所以，a2%a2£(/(a1,a2,a3))=|e(/(a1,a2,a3))|=lf(x,,x2,x3)-f(a,,a2,a3a2%a2e(a,)+—e(a2)+-^^-+1e(a3)e(%)ae(%)e(q)+—e(a>)+—--+1a3a；《(空+旦+等+I)xki(r2a3%a3 2=e(a2)=e(a3)=则相对误差限为

下略。解二：根据函数y=/（X，X2,…,七）的函数值的绝对误差e(y)=6(/(百/2，毛,…,天))=相对误差，(y)=簿(7(芭,*2，演,…，怎))=Xv-— 7^(玉)1=1ox,./(x1,x2,x3,---,x„)公式计算。1、用秦九韶算法的多项式格式乘法计算多项式P（x）=x7-2x6-3x4+4x3-x*2x2(12x2(1+2x)(1+x)在x=2处的值p（2）o2、利用等价变换使下面式子的计算结果比较精确。3、指出下列各题的合理计算途径（对给出具体数据的，请算出结果）[1]I-cosl°（三隼函数值取四位有效数字）[2]ln（3O->/3O2-l）（对数函数值取六位有效数字）[3]匕土（其中x的绝对值很小）sinx[4]x127100 1[5]y—!—占〃（〃+1）4、设近似值T0=So=35.7O具有四位有效数字，计算中无舍入误差，试分析分别用递推式Tm=57；-142.8?0S,+1=15,-142.8计算T20和S20所得结果是否可靠。5、计算p6（x）= —3/+x?-4x+l的值凡（3）o1、2、1、2、3、[1]1-cosx=2sin2 6 15 45 136 404 1213所以〃6(3)=12132 6 15 45 136 404 1213所以〃6(3)=1213。习题二解答1.用二分法求方程x3-2x2-4x-7=0在区间［3,4］内的根,精确到10月即误差不超过Lx10一3。分析：精确到与误差不超过IO,不同。30-V302-l= J =0.01667,30+V302-lln(3O-73O2-l)=-4.09414—cosx sinxx—； = =tan—sinx 1+cosx 2[4]x,27=x-x2-x解：因为f(3)=-10V0,f(4)=9>0,所以，方程在区间［3,4］上有根。由-x8-x,6解：因为f(3)=-10V0,f(4)=9>0,所以，方程在区间［3,4］上有根。由[5]由小到大依次相加。鸟1鸟4 1、， 1 100仁心+1)占〃72+1 1011014、设计算Ti的绝对误差为e(Ti)=Ti*-Ti,其中计算T。的误差为£,那么计算T20的误差为e(T2O)=T2o*-T20=(5T19*-142.8)-(5T19-142.8>5(T19*•T19)=5e(T19)=52e(Tl8)=52Oe(To)显然误差被放大，结果不可靠。同理，e(S2o)=[；)e(S0),误差缩小，结果可靠。5、解：将所给多项式的系数按降哥排列，缺项系数为0o2 0 -3 0 1 -4 1|x*fKbn-an|x*fKbn-an_b-a_4-3=—<-xl0-3

T2有2*1>1000,又为21。1024>1000,所以n=ll,即只需要二分11次即可。列表讨论如下:nanbnXnAW的符号1343.500—23.50043.750+33.5003.7503.625—43.6253.7503.688+53.6253.6883.657+63.6253.6573.641+73.6253.6413.633+83.6253.6333.629—93.6293.6333.631一103.6313.6333.632+113.6313.6323.632—x*-xn=3.632o指出：(1)注意精确度的不同表述。精确到103和误差不超过1J3是不同的。(2)在计算过程中按规定精度保留小数，最后两次计算结果相同。如果计算过］程中取4位小数，结果取3位，则如下表：n%Xn跖)的符号1343.5000——23.500043.7500+33.50003.75003.6250—43.62503.75003.6875+53.62503.68753.6563+63.62503.65633.6407+73.62503.64073.6329+83.62503.63293.6290—93.62903.63293.6310—103.63103.63293.6320+113.63103.63203.6315—(3)用秦九韶算法计算f(xn)比较简单。

1*.求方程x22x2-4x-7=0的隔根区间。解：令y=V—2/—4尤一7,贝ij/=3x2-4x-4=3r+2K-2、TOC\o"1-5"\h\z, 2当y，=3x2-4x-4=a+2K-2—0时，有％=——x2=2函数单调区间列表分析如下：X2(-8,一§)_22"32(2,+8)y+0—0+y——弛27-15TOC\o"1-5"\h\z2 149 一， 、 2.因为y——〜 <0 —15<0>所以方程在区间――2、匕无根；3 27 32 149 ,,, 2因为y—― <0,而函数在^8—、上单调增，函数值不可能变号，所以3 27 3方程在该区间上无根；因为丁冬三一15<0,函数在(2,+8)上单调增，所以方程在该区间上最多有一个根，而(3)=-10<0,y(4)=9>0,所以方程在区间(3,4)有一个根。所以，该方程有一个根，隔根区间是(3.4)。2.证明l-x-sinx=0在［0,1］内有一^1^根，使用二分法求误差不大于gxlO"的根，需要迭代多少次？分析：证明方程在指定区间内有一个根，就是证明相应的函数在指定区间有至少一个零点。解：令/(x)=1-x-sinx,g]^/(O)=l-O-sinO=l>O,/(l)=l-l-sinl=-sinl<O>由零点定理，函数f(x)在［0,1］区间有一个根。由I*i,b“一a“b—a1—0 1 1 _4In\2 2" 2" 2" 2W2n''>10000,又为2i°=1024,213=8192<10000,214=16384>10000所以n=15,即需要二分15次。指出：要证明的是有一个解而不是唯一解，因此不必讨论单调性。

20 ， ..试用迭代公式 ,x0=l,求方程/+2/+10%-20=0的X/.+2/+10根，要求精确到IO"。分析：精确到10-5即误差不超过LxlO-52解：令/3)=丁+2/+10工-20列表进行迭代如下：Xk01-711.538463.7596421.29502-1.5238031.401820.7031141.35421-0.3066751.375300.1372161.36593-0.0606771.370090.0270581.36824-0.0119891.369060.00531101.36870-0.00228111.368860.00110121.36879-0.00038131.368820.00025141.368813992x10-5151.3688139'92xl0-5指出：精确到10”可以从两个方面判定。第一，计算过程中取小数到IO一位，最后两个计算结果相同，终止计算。第二，计算过程中取小数到I。-',当k+Lx/vgxlO-终止计算。本题采用第一种方法。

.将■元非线性方程2cosx-e*=0写成收敛的迭代公式，并求其在受二。9附近的根，要求精确到IO-。TOC\o"1-5"\h\z解：2cosx-e*=0改写为2cosx=e*=幺吧=1=幺凶-1=0,则e ex2cosx,、几X=XH 19攻, ,-2sinxe*—2cosxe*gi—A—2期x+cosl2&sinr+-

, ,-2sinxe*—2cosxe*gi—A—在/=05,处，因为2>/2sin-&5'+—gN5y1 ——=09-615<1e所以迭代法8多+广4+幺答-1在/=05.的邻域内收敛。ek列表迭代如下：王00.510.7120.6930.69此时2COS06'9-e°69=o0'0614。.为求方程丁-炉_1=。在/=15.附近的一个根，设将方程改为下列等价形式，并建立相应的迭代公式：d5=1+4■迭代公式X*+]=1+]<X XkTh3=1+x2迭代公式冬«=一!一迭代公式x*+[= -气一共试分析每种迭代公式的收敛性，并取一种公式求出具有4位有效数字的近似值。解：（1）因为x=l+-l，所以迭代函数为g字〜1+4,贝IJX Xg‘个〜,夕65牛m-2x1寸卜二=—\—<1满足局部15 3375收敛性条件，所以迭代公式=1+3具有局部收敛性。1\_（2）因为工=4/3,所以迭代函数为gj-小/。,则gA-—A^-X212x=—X-++X2J=-—彳,3 3 -3改％2d7v1<|g45T 2=0456<l满足局部收敛性条件,所以迭代公式3小15232七+1=*片3具有收敛性。（3）因为》=—L丁，所以迭代函数为g产L丁，则停一1d 字fTOC\o"1-5"\h\z, 1 4-1 1 4P 二—**%—1 ———1，2 2.2 1仁人5TL1=—!~丁=1414>1不满足收敛性条件，所以迭代公式2x05、xk+]=——--；■不具有收敛性。哮一产用迭代公式X』=1+4•列表计算如下:4Xk01.511.44421.48031.45741.47151.46261.46871.46481.46791.465101.466111.465所以，方程的近似根为x*=1465。6.设'x+O*2-3、，应如何取C才能使迭代公式x*+]=e/、具有局部收敛性？解：设C为常数，因为_3~,所以―-术、l+2Cx,要使迭代公式具有局部收敛性，需忸'咻卜|1+2飙|<1,此时即有-1<1+25<1，也即-l<Cxn<Oo即只要C取满足如上条件的常数，就可以使得迭代公式具有局部收敛性。指出：本题的一般形式为：设晒x）=x+cy（x）,应如何取c才能使迭代公式x*+i=of、具有局部收敛性？显然，奴工）=*+守*）是迭代格式4+1=0^、相应的迭代函数，因此该迭代格式要求解的方程是工二以好二工=彳+守*）=/（x）=0o也就是说，这是如何选择C,构造一个求解方程f（x）=O的收敛的迭代格式的问题。因为以x）=x+cy（x）,所以夕'（x）=i+cy'（x）,要使迭代格式收敛，需I"⑸=|1+仃'（何<1解之得-2<OV）<0,即C与/'（X）异号，且|仃'"）|<2。下面的讨论利用了本题的特殊条件，求出了具体的结果：因为。仔-x+Cf2_3、，所以当x=x+C(x2-3)时，WC(x2-3)=0,则X=±G,即函数-x+C<*2—3、的不动点为x*=±G。而+2cx,根据局部收敛性定理，当忖网=1+2。叫<1=>W<c<0时，迭代格式收敛到百；当忖(-百)卜卜+2。(-6)|<1=>0<。<\时，迭代格式收敛到-G。7.用牛顿法求方程Y-3x-1=0在初始值%=2邻近的•个正根，要求k+[-x/<ioo解：因为X3一31一1二0所以有f(x)=d—3x—i,相应的迭代公式为X：— —1 2x；+1**m—%一-3--3―3片-3取x0=2为迭代的初始近似值。迭代的结果列表如下：k0123Xk21.88891.87951.8794因为|七-司=0.0001<|xio-3,符合计算的精度要求,所以x=x3=1.8794o.用牛顿法解方程」-c=0,导出计算数c的倒数而不用除法的一种简单X的迭代公式。用此公式求0.324的倒数，设初始值%=3,要求计算有5位有效数字。解：对于方程^•-c=0,有/'cyL-c,相应的迭代公式为X X应用该迭代公式求0.324的倒数，列表计算如下%0313.08423.086433.0864所以一^"30・864。0324指出：如果将方程，-c=0改写为等价的5-1=0,则有f^cx-\,相应的迭X代公式为无法展开迭代。.设a为已知数，试用牛顿法导出求标的迭代公式，并求极限..标一/+]hm—I00解：设a为正实数，n为自然数，由牛顿法，方程/小。=0的解为x[-a_nx^一%,+4Xkn—\ n-Inxknxk_(〃-1)x(+a;而L=匕(〃-1出+普]

« 4此即求标的迭代公式。由此，贝IJ

lim&T8标一%+1

(4a-xk)2加一[5_1)/+y/a--[(n-l)xklim&T8标一%+1

(4a-xk)2=lim n1— ———=lim 产 』 (五-々)2 - (^/a-xk)2—[(〃―1)+。(1_ 〃]=lim 7= --[(n-l)+a(l--=lim --[(n-l)+a(l--=lim f= 一 —2(喝一七)-1。(1-〃)(-〃)X”]=lim- =limk—^oo-2x(-1)a(l-n)_a(\—n)_a(l-n)l-n2球”-2lim球”一2(标)"+i―研♦—>8指出：本题中，表面上是Z->8的问题，但实际上却是X*->标的问题，引,乙+1才是极限过程中实际的变量。本质上。本题实际上是求极限布—X 标一”〃-1)々+第] 布-1[("-%.+ox：-"]lim—7=———=lim 产 ——=lim 产 5(1"/)2-(Va-Xj2 i>[a——[(〃-1)1+奴『”]=蚂—由于讨论的是:型不定式，且不定式的分母上有2次的"因子，因此两次应用罗必塔法则。解二：首先证明一个定理：定理：设/(x*)=OJ'(x*)HO,又设f(x)在X*的某个邻域内具有连续的二阶导数，则牛顿迭代法具有局部收敛性，且有。吧(—3证明：因为g。)“明所以g'a所以g'a)=HS/(x)/"(x)

(八X))2因为f(x)在邻域内具有连续的二阶导数，所以g'(x)在邻域内连续，且g'(x*)=/(x*)/"(x*)

g'(x*)=/(x*)/"(x*)

(/V))2=0由局部收敛性定理，牛顿迭代法具有局部收敛性。对g，(x)=(x_累)=；；；/；；?求导,根据条件有g"(x*)=f*(x)f\x)由收敛阶定理，若/则8〃。*)=孕2。0,牛顿迭代法二阶收f(x)敛，若/(x*)=0,则g"(x*)=，2=0,牛顿迭代法有更高的收敛阶。/(X)因为牛顿迭代法有二阶收敛性，所以■/V)H =g"(x*)=/'「*)=/,a*)J8(x*")2—2!— 2!—2/'(x*)显然如果x*是方程f(x)=0的单根，贝1」/(犬)=(犬-％”)双工)，且到x,)h0。此时f\x)=(p{x}+(x-x)<pr(x),则f\x)=奴x")H0,可见定理中的条件"/(■/)=0,/(/)。0"可以等价替换为“丁是方程£&)=0的单根”对本题来说，f(x)=xn-a,工*=标是方程的单根,所以f\x)=nxn-',f\^a)=〃(标厂tH0f*(x)=n[n-V)xn~2,/*(Va)=“(〃-1)(布)"'2yfa—xk+i_..xk+i—\[a_n{n—\)(y[a.y'~_n—\_1—

丽-x*)2-J8(“布)2- 2〃(标尸-2加一2标指出：应用分组分解法进行因式分解，分子、分母约去"0"因子，就可以按连续函数的极限性质求解了。10.用快速弦截求方程丁-3*-1=0在初始值凝=2邻近的实根(取％=1.9,要求精确到IO-)。解：因为彳3-3犬一1=0

所以有F(x)=x3-3x-1,相应的迭代公式为.—xk " (xk— )7(演)-/(Xj)取xo=2为迭代的初始近似值。迭代的结果列表如下:kXkXk-Xk-1f(Xk)f(Xk)-f(xk.i)021I1.9-0.10.159-0.84121.8811-0.01890.0130-0.14631.8794-0.00170.0001-0.012941.8794因为卜4-天|=。0°°°<3*10-3,符合计算的精度要求,所以x~x4=1.8794o指出：本教程所说快速弦截法是通常所说的弦截法(割线法)，而它所说弦截法是通常的单点弦截法。11、分别用下列方法求方程4cosx=e，在邻近的根，要求有三位有效数字。7T(1)用牛顿法，取/=彳；(2)用弦截法，取/=?,X]='；(3)用快速弦截法，取解：方程4cosx=e*变形为e*-4cosx=0,贝ij/(x)=e'-4cosx,f\x)=e*+4sinx。牛顿法、弦截法、快速弦截法公式分别为(1)牛顿法_ve"-4cos%

f\xk)\eXi+4sinxt.(2)弦截法/(xj+l.81(占/(xj+l.81(占-0.785)；(3)快速弦截法f(xk)取3位有效数字，分别计算得kxk牛顿法弦截法快速弦截法00.7850.7850.78511.591.571.5721.411.331.3331.391.401.3841.391.381.4051.391.3961.381.39补充题(一)1、确定方程x5+x-10=0的根的个数，找出隔根区间。2、用二分法求方程f(x)=X3+2x-5=0在［2,3］的根的近似值，要求误差不超过0.005o3、用二分法求方程f(x)=x3-2x-5=0在［2,3］的根的近似值，要求误差不超过0.05。24、用二分法求方程/(x)=sinx-工=0的非零实近似根，使误差不超过104-2O5、分析方程/(x)=sinx-］=0的根的分布情况，并用二分法求正根的近似值，使误差不超过102O6、估计用二分法求方程f(x)=x3+4x2-10=0在［1,2］内的根的近似值，为使误差不超过104时所需要的二分次数。分析与解答1、■$,y=x5+x-10,y=5x4+1,显然y'>0,而且函数没有不可导点，所以,函数在区间(-叫8)上是单调增的，故方程最多有一个根。因为因0)=-10<0,y(2)=24>0,所以方程在(0,2)区间有一个根，(0,2)即为方程的隔根区间。

2、因为f(2)=7>0,f(3)=28>0,实际上本方程在指定范围内无根。但如果不加判定，也可以计算出一个值来。所以，用二分法求方程的根必须先行判定。要特别注意的是，完整的二分法的过程是，第一步代入初值，第二步判断是否有解，第三步在有解的前提下求出解来。不进行判断就形式地套用二分法的过程是不可以的，同样地，如果因为无解就放弃讨论也是不正确的。3、因为f(2)=-1<0,f(3)=16>0,所以方程在区间上有解。x*-x|<^ZS=^Z^=2z2=_L<o.O5，所以，2。>20,n=5o“I2 2"T2" 4x*-2.104、画出y=sinx和y=土的曲线，可以看出，4两条曲线除了原点外，在第一象限有且只有一个交点。交点的横坐标介于1.5与2之间(显然，/万、2(幸71/2-1.5,sin(兀/2)=1,——<1,所以在兀/2点，4f(x)>0,而当x=2时，一=1,sinx<1,所以在2点，4f(x)<0o5、画出产sinx和y=3的曲线，可以看出，两条曲线除了原点外，在第一象限有且只有一个交点。交点的横坐标介于1.8与1.9之间(根据图象，用计算器计算估计，当sinx的值从大于二的值变为小于时，隔根区间就找到了)o2要求|X*一Xn|W0.01,可以求出用二分法计算的次数。在区间［1.8,1.9］上，因为। I —G b—a 1.9—1.8 0.1X*区_il= = =—<0011 "I 2 2" 2" 2"所以，n=4o具体计算过程如下nan%XnRXn)的符号11.81.91.85+21.851.91.875+31.8751.91.8875+41.88751.91.89+所以，x*-4=1.89指出：确定求根区间和根的初始近似值，应用M\TLAB工具，用交轨法是重要的途径，可以先确定大致范围，再缩小区间重新画图精细化。在用普通的手工画草图的方法画交轨图的时候，可以借助于计算器使得隔根区间更短，但这种方法只对简单问题有效。6、|x*-xn|<10-5,BP^=-!-<10-s，所以2905。1 2 2因为-5=32768,2^=65537,217=131072,所以n=17。(二)1、对于方程3x2-ex=0,为求最大正根与最小正根的近似值，试分别确定迭代函数g(x)及区间［a,b］,使得当x()G［a,b］时，相应的迭代过程Xk+尸g(x0收敛到要求的根。2、证明：当x0=1.5时，迭代法都收敛于方程f(x)=x3+4x2-10=0在区间［1,2］内的唯一实根X”,分别用上述迭代法求满足精度|Xk+i-Xk|WO-的近似根。3、为求方程f(x)=x3-x2-1=0在xo=L5附近的一个根，可将方程改写成下列等价形式，并建立相应的迭代公式口］改写成x=l+4,迭代公式为%=1+斗；X xk［2］改写成x3=l+x2,迭代公式为x,+1=而W；［3］改写成产=」二，迭代公式为%产二=。・'-1 Jx*-1试分析每一种迭代公式的收敛性。分析与解答1、根据3x2和e'的图像可知，方程3x2-eX=0在实数域上有三个根，分别在区间(-1,0),(0,1),(3,4)内。其最大正根在［3,4］区间，最小正根在［0,1］区间。取迭代函数g(x)=ln3x2,可以得到最大正根，而取迭代函数g(x)=7?不，可以得到最小正根。2、两种迭代法的迭代函数分别在区间［1,2］和［1,1.5］上满足定理2(不动点原理)的条件，故当x0=1.5时两种迭代法都收敛，且分别迭代9次和25,都可得到近似根1.36523。我们讨论第一种迭代法，用定理2证明。它的迭代函数为g(x)=、户。V4+x首先，g(x)是一个减函数，当x=l时，g(1)=>/2，当x=2时，g(2)=Jg所以当x£[l,2]时，l<g(2£g(x)Wg(l)<2,即g(x)w[l,2]。其次，g\x)=-^.1 ,显然这是一个增函数，当x=2时，其函数2J(4+x)3值为全，所以，g()</⑵<1。指出：只给出了含根区间，就只能用定理2证明。3、[1]给出了初始近似值，也即知道了精确根的大致位置，可以用定理4(局部收敛性定理)证明。由题意，方程有实根。下面证明g/(x)连续和g/(x*)<l(x*是方程的精确根)o方程g(x)=l+jg'(x)=5，可见g/(x)在1.5及其附近是连续减函数，因为X Xg(1.5)= -0.59,1.5又在x*的邻域内，由函数g'(x)的连续性，g'(x*)<l,所以此迭代法具有局部的收敛性。指出：一般地说，用定理2(不动点原理)证明只要利用函数的单调性与区间上的最值就可以讨论，而用定理4(局部收敛性定理)则需要用到函数的连续性。习题三解答1、用高斯消元法解下列方程组。TOC\o"1-5"\h\z2%]-x2+3x3=1 ①<4X]+2x2+5x3=4 ②X]+2x2 =7 (3)解：②+(-1)X①，③+(-；)X①消去第二、三个方程的X,,得：2xt-x2+3x3= 1④4x2-x3=2 ⑤-x2--x3=—⑥I2-2 25再由⑥+(-1)x⑤消去此方程组的第三个方程的三，得到三角方程组:2%-x2+3x3=14x2-x3=27 21-x3=8 4回代，得：x3=-6,x2=-1,—=9所以方程组的解为x=(9,-l-6)r注意：①算法要求，不能化简。化简则不是严格意义上的消元法，在算法设计上就多出了步骤。实际上，由于数值计算时用小数进行的，化简既是不必要的也是不能实现的。无论是顺序消元法还是选主元素消元法都是这样。②消元法要求采用一般形式,或者说是分量形式，不能用矩阵，以展示消元过程。要通过练习熟悉消元的过程而不是矩阵变换的技术。矩阵形式错一点就是全错，也不利于检查。一般形式或分量形式：TOC\o"1-5"\h\z2X1-x2+3x3=1 ①4X1+2x2+5x3=4 (2)%+2x2 =7 (3)矩阵形式'2-1'2-142J2向量形式③必须是方程组到方程组的变形。三元方程组的消元过程要有三个方程组，不能亦形出单一的方程④消元顺序占7%二…，不能颠倒。按为支援在方程组中的排列顺序消元也是存储算法的要求。实际上，不按顺序消元是不规范的选主元素。⑤不能化简方程，否则系数矩阵会变化，也不利于算法设计。TOC\o"1-5"\h\z11%)—3x2—2x3-3 ①--23^+1lx2+x3=0 ②X[+2x2+2Xj=—1 ③解：②+(停X①，③+(W)x①消去第二、三个方程的占，得:TOC\o"1-5"\h\z1lXj-3x2— 2x3= 3 (4)52 35 69 0, 1产，产=u⑤25x2+2\3=-14 ⑥I 11-113 1125再由⑥+(-聂)x⑤消去此方程组的第三个方程的三，得到三角方程组:IT1 —3x2— 2x3= 352 35 69< —羽 占=—11 2 11 3 11193 223 X-,- I 52 52回代，得：_223 _1_0£ __4_1_/-~193，%2-T93，X'-797，所以方程组的解为_41106 223rA—197,797，-?972、将矩阵’1020、0111A=20-11(001 1)作LU分解。解：设<1020、00“I(u\\U\2〃13W14A=01I1—*21100°U22M23M24=LU20一I1，32IO°0U3\〃32<001bJa1M2“431000“44J根据矩阵乘法，先求u的第一行，由ulj9得Uy।—1,Uy2=0,〃]3=2,〃]4=0。再求L的第一列，由矩阵乘法，因为％=/“%，所以/“=%,而%=1,所W1I0o以/升=%，所以==0,4=2,/41再求U的第二行，得0o/2(XW12+1X〃22=1，贝U

u22=1—/21xw12=1—0x0=1,/21xw13+lxw23+0xw33=1,贝Ij优23=1—Z21XW13=]_0x2=],Z21xu14+lxw24+0xwj4d-Oxw^=1,则u24=1-/2]x«14=1—0x0=1,再求L的第二列，得xw12+Z32xw22+1x0+0x0=0,贝ij/32=0—/31xw12=0—2x0=0/4,xWj2+Z42xw22+Z43x0+0x0=0,则/42=0—/41xul2=0—0x0=0再求u的第三行，得，31X"13+’32X“23+1XW33=-1，则W33=-1一131X43132X“23二-1-2X2+0X1——5/31xw]4+/32xu24+lxw34+Oxw44=1,贝Iju34=1-/3Ixw14-/32xw24=1-2x0-0x1=1再求L的第三列，得/4IXW134-/42XW23-I-/43XW334-1X0=1,则/43=-1x(l-0x2+0xl)=-1再求U的第四行，得加xw14+Z42xw24+Z43xm34+lx%=1,则W44=1-，4]X/4+I42X〃24+，43XM34=1-0x0-0x1— —Xl)=—所以，矩阵A的LU分解为：’1000、’1020、01000111L—20 10,u=00-5100--1000-1 5)1 5)指出:用分数而表示元素，不能化成近似小数也不化成小数表示。3、用LU分解紧凑格式分解法解方程组。91089108668710A7解一，用-一般格式求解：将系数矩阵作LU分解得:1000、’561000571,u=I0052103I051Ly=b方程组为’1000、792_45~50-50 0-31710-ioo(y^5 „3 1乂J10-1I5)解之得(1、Iio同样地，解方程组取=丫得20、20、-12-53；解二，用LU紧凑格式分解法求解:对增广矩阵三角分解：7910810910877655796_2_45一5一571-552310567\51、11L56575I1010105657517251~2010 1、'5 7 9 10 1、、1—3—6_2_4_355-5-5一-5_17_1T7」 _17-2-25~2~一万~2,c 3 1 35 11 0 - ————71 5 1010J7651610-3原方程组化成同解的上三角方程组为:5x1+7x2+9x3+10x4=11 3~X4=~~1010回代得x=(20,-12,-5,3)ro指出・紧凑格式是直接应用公式进行计算，计算结果保存在A的相应元素位置。从算法的角度，紧凑格式实际体现在数据的存储方法上。由于紧凑格式计算时不再需要A的前面的元素，因此可以进行。4、用列主元的三角分解法解线性方程组。{一+2x?—2七二—13x,-x2+4x3=72x}-3x2—2x3=0解一，列选主元素消元法：先选第一列主元为=3,将第一个方程与第二个方程交换，消去不得：3西-x2+4x3=75 2 4〈 —X) =—3 2 33 37 14 14I3 33再选第二列主元为七3交换第二、三两个方程，消去得三角形方程组:TOC\o"1-5"\h\z3%j—x24-4x3= 77 14 14〈 X、 X,= 3 2 3 3 312 6I 3k3回代求得方程组的解玉=5，%=1，%=2所以方程组的解为17x=(2,l,-)ro解二，列主元素三角分解法：(A,b)=3.22-2-24-3-2-2(A,b)=3.22-2-24-3-2-242-2-3-2-24 72—2—1-3-20-2_73_574143-214-1-2_73_574 714_L4T-H-4 -2同解的三角形方程组为3X1-x2+4x3=7 14 143 3 3-4x3="2回代求得方程组的解演=g，x2=l,石=2所以方程组的解为ITx=(2,l,-)\说明：用矩阵讨论中，矩阵元素进行了化简。5.用追赶法解方程组-1-12-1-12-1,b=T

000

⑼分析：三对角矩阵r,.可以分解如下形式的两个矩阵：a,YiAa23I»2141Am3utAw,4u3由矩阵乘法规则，有w,=%,=~^-(i=2,3,…〃) ,kw,=%-仍t(i=2,3,…,〃)这样可以求出矩阵L和U的所有元素。设有系数矩阵为A的方程组：Ax=b,b=(bt,b2,---bny,这样的方程组称为三对角方程组。三对角方程组经LU分解分解为Ly=b,Ux=y,求解之”=瓦yt.=bi-liyi_l,i=2,3,---,n'x“=%.X,=(^-笈苞+])ui,i=n~l,n-2,---,\这就是所谓追赶法。解：由公式Wj=fXj—2,/?— —二-——77,TOC\o"1-5"\h\zu}2 21 3m2=4一/#\=2-(--)x(-1)=—/产=14\o"CurrentDocument"m22 322 4M3~a3~kfll=2-(_§)X(-1)=]3)7一彳二3TOC\o"1-5"\h\z3 5M4=(X4_/血-2-(--)X(-1)=—4 4,一21一3」«4A544 6“5~a5~4A=2-(-y)X(-l)=y由此得下三角方程组(1~2(1~2必%>4晃1‘I、000O和上三角方程组2-132解上三角方程组145145代入并解上三角方程组2-132-1-1-12-132-1-1-1‘X、“2*4\x57J_2J_3]_45>34566、用改进的Cholesky分解法解方程组[4xt-2x2-4x3=10{-2%+17%+10x3=3[-4X]+10x2+9x3=-7解：设此方程组的系数矩阵为A,右端向量为b,则'10、,b=3'10、,b=37-21710-4109,矩阵A是对称正定矩阵，可以进行乔累斯基分解。4一2-4、“11U\2/3,-21710二U12U22U22“23-410%kW13