2022秋季學(xué)期 廣東工業(yè)大學(xué) 《數(shù)值計(jì)算》期末習(xí)題集

習(xí)題一1.什么叫數(shù)值方法？數(shù)值方法的基本思想及其優(yōu)劣的評(píng)價(jià)標(biāo)準(zhǔn)如何？答：數(shù)值方法是利用計(jì)算機(jī)求解數(shù)學(xué)問題近似解的方法。2.古代數(shù)學(xué)家祖沖之曾以EQ\*jc3\*hps23\o\al(\s\up7(355),113)作為圓周率π的近似值，問此近似值具有多少位有效數(shù)字？EQ\*jc3\*hps23\o\al(\s\up7(325),133)*xx=π35517該近似值具有7為有效數(shù)字。1.試構(gòu)造迭代收斂的公式求解下列方程：解：k0123xk00.250.250980.25098**局部收斂EQ\*jc3\*hps23\o\al(\s\up30(x),k)00)'局部收斂k0123456789xk**2.方程x3x21=0在x=1.5附近有根，把方程寫成三種不同的等價(jià)形式：13x=3x=2x=2x=;2，對(duì)應(yīng)迭代公式xk+1=;11xk1x1判斷以上三種迭代公式在x0=1.5的收斂性，選一種收斂公式求出x0=1.5附近的根到4位有效數(shù)字。解：1x2（2）Q(x)=（3）Q(x)=(x)=x(1+x2)(x)=(x1)Q(x0)'<1局部收斂局部收斂Q(x0)'<1Q(x0)'>1不是局部收斂迭代公式（10123456781.444441.479291.4569761.471081.462091.467791.44161.4664791.46501.465931.46531.465721.465481.465631.4655341.465595**迭代公式（2k0123456xk**3.用牛頓法求方程f(x)=x3一2x2一4x一7=0在[3,4]中的根的近似值（精確到小數(shù)點(diǎn)后兩位）。f'(x)=3x2一4x4EQ\*jc3\*hps13\o\al(\s\up13(3),k)EQ\*jc3\*hps23\o\al(\s\up8(一),3x)EQ\*jc3\*hps23\o\al(\s\up8(一),4x)EQ\*jc3\*hps23\o\al(\s\up8(4),k)x7k0123xk3.53.643.633.63**4.試證用牛頓法求方程(x一2)2(x+3)=0在[1,3]內(nèi)的根x*=2是線性收斂的。解：y次迭代公式(xk2)(xk+y次迭代公式(xk2)(2xk+1)(xk2)(2xk+1)=*xxk+1(xk2)(xk+3)故k喻2=xk2*xxk從而ekkekek喻2故牛頓迭代公式是線性收斂的的迭代公式，并討論其收斂性。的迭代公式，并討論其收斂性。解：f(x)=x3一af'(x)=3x2 kkx3a kk相應(yīng)的牛頓迭代公式為EQ\*jc3\*hps13\o\al(\s\up2(2),k)相應(yīng)的牛頓迭代公式為迭代函數(shù)(x)=2xEQ\*jc3\*hps23\o\al(\s\up12(3一),3x3)2a，Q''(x)迭代函數(shù)則則故迭代公式收斂，且收斂階為2階。1.設(shè)有方程組EQ\*jc3\*hps13\o\al(\s\up12(1),x)EQ\*jc3\*hps13\o\al(\s\up12(2),x2)(1)考察用Jacobi法，Gauss-Seidal法解此方程組的收斂性；(2)用Jacobi法及Gauss-Seidal法解方程組，要求當(dāng)*(k+1)-*(k)<10-4時(shí)迭代終止。故用雅克比法及高斯-塞德爾法解此方程均收斂。=-x2-x3-3雅克比法：xk+1)=-xk)-xk)-，xk+1)=xk)-xk)+5，取初始向量x0)=x0)=x0)=0，迭代18次有xEQ\*jc3\*hps13\o\al(\s\up4(1),i)8-xEQ\*jc3\*hps13\o\al(\s\up4(1),i)7<10-4（i=1,2,3）高斯-塞德爾法：xk+1)=-xk)-xk)-，xk+1)=xk)-xk)+5取初始向量x0)=x0)=x0)=0，迭代8次有xEQ\*jc3\*hps13\o\al(\s\up4(8),i)-xEQ\*jc3\*hps13\o\al(\s\up4(7),i)<10-4（i=1,2,3）2.用選列主元高斯消去法求解方程組「3-14-12-2L2「3-14-12-2L2-3-2|||7]|(AD)喻|0喻|0|||0解得---|喻|0-1]「5]1-3x2-2x3=0解：|3|3|3|-147|T-153742|-3334-3-43.用三角分解法求解方程組EQ\*jc3\*hps13\o\al(\s\up12(解：系數(shù)矩陣三角分解為：EQ\*jc3\*hps23\o\al(\s\up11(2),4)EQ\*jc3\*hps23\o\al(\s\up11(4),8)EQ\*jc3\*hps23\o\al(\s\up11(8),1)EQ\*jc3\*hps23\o\al(\s\up11(1),2)EQ\*jc3\*hps23\o\al(\s\up11(0),1)EQ\*jc3\*hps23\o\al(\s\up11(0),0)EQ\*jc3\*hps23\o\al(\s\up11(-2),0)EQ\*jc3\*hps23\o\al(\s\up11(4),0)EQ\*jc3\*hps23\o\al(\s\up11(8),3)原方程可表為：EQ\*jc3\*hps23\o\al(\s\up3(2),6)||||01-1|2||||01-1|2(5-2-10)T解|y2||解|||10-32||x2x3」|y1|「|y1|「5]0]Ly3」「-20「-20||EQ\*jc3\*hps27\o\al(\s\up6(2),1)EQ\*jc3\*hps27\o\al(\s\up6(9),9)EQ\*jc3\*hps27\o\al(\s\up6(1),0)EQ\*jc3\*hps27\o\al(\s\up6(2),9)EQ\*jc3\*hps27\o\al(\s\up6(1),5)EQ\*jc3\*hps27\o\al(\s\up6(5),38)-1.5316,0.2211,0.1316TEQ\*jc3\*hps23\o\al(\s\up9(1),1)解：104)1)|計(jì)算cond(A)104)A-1||1104]104-11EQ\*jc3\*hps23\o\al(\s\up7(1),1)EQ\*jc3\*hps23\o\al(\s\up7(0),0)EQ\*jc3\*hps23\o\al(\s\up12(4),4)EQ\*jc3\*hps23\o\al(\s\up9(1),1)1.給出概率積分f(x)=e-x2dx的數(shù)據(jù)表：試用二次插值計(jì)算f(0.472).X0.460.470.480.49f(x)0.48465550.49375420.50274980.5116683解：取插值節(jié)點(diǎn)：x2L2(x)=EQ\*jc3\*hps23\o\al(\s\up11(2),Σ)yili(x) (x-x1)(x-x2)(x-x0)(x-x2)=y0(x0-x1)(x0-x2)+y1(x1-x0)(x1-x2)+y2 x-x0)(x-x1)x2-x0)(x2-x1)f(0.472)=L2(0.472)=0.49556162.已知y=sinx的函數(shù)表Xsinx0.997490.999570.99166試構(gòu)造出差商表，利用二次Newton插值公式計(jì)算sin(1.609)(保留5位小數(shù))，并估計(jì)其誤差.解：由題意得如下差商表kxkf(xk)fx0,xkfx0,x1,xk11.60.999570.0221.70.99166-0.02915-0.49950又f(x)=sinx,f(3)(x)=-cosxf(3)(ξ)<maxcosx<0.128843.用最小二乘法求一個(gè)形如y=a+bx2的經(jīng)驗(yàn)公式，使與下列數(shù)據(jù)相擬合X25313844Y32.349.073.397.8解：依題意x2解：依題意xiiEQ\*jc3\*hps13\o\al(\s\up10(2),i)EQ\*jc3\*hps13\o\al(\s\up8(4),i)C0=yiQ0ii正則方程為〈01正則方程為〈01故擬合曲線為y=0.973+0.05x21．試確定下面求積公式1f(x)dx~C[f(x0)+f(x1)+f(x2)]使其具三次代數(shù)精度.解：要公式有3次代數(shù)精度，需有EQ\*jc3\*hps13\o\al(\s\up5(2),0)EQ\*jc3\*hps13\o\al(\s\up5(2),1)EQ\*jc3\*hps13\o\al(\s\up5(2),2)EQ\*jc3\*hps13\o\al(\s\up6(3),0)EQ\*jc3\*hps13\o\al(\s\up6(3),1)EQ\*jc3\*hps13\o\al(\s\up6(3),2)22故求積公式為1f(x)dx=[f(0)+f()+f(-)]2．分別用復(fù)合梯形公式及復(fù)合Simpson公式計(jì)算dx,(取步長(zhǎng)h=1/6).解1）用復(fù)合梯形公式a=1,b=2,h=故N=6f(x)=f(a)==1.4427f(a+)= f(a+)=f(a+)=f(a+)= ln(713) ln(512) ln(813) f( f(b)==1.8205f(xn)=16.3582dx~[f(a)+f(b)+2f(xn)]=1.6351167（2）用復(fù)合Simpson公式：f(x12)=f(a+)=1.4760f(x32)=f(a+)==1.5414f(x52)=f(a+)=1.6055f(x72)=f(a+)==1.6683f(x9)=f(a+9)=714=1.72992f(x11)=f(a+11)=23112=1.79052f(xn+1)=39.2464n=02∫f(x)dx心[f(a)+f(b)+4f(xn+12)+2f(xn)]=1.63521673．試構(gòu)造兩點(diǎn)Gauss公式w1f(x)dx心A0f(x0)+A1f(x1),并由此計(jì)算積分(精確到10-4)解：二次Lagendre多項(xiàng)式：2-33dxEQ\*jc3\*hps23\o\al(\s\up22(A),A)EQ\*jc3\*hps13\o\al(\s\up19(0),1)EQ\*jc3\*hps13\o\al(\s\up15(1),1)