哈工大研究生數(shù)值分析報告精彩試題及問題詳解

1、word1 . x3,2分別是方程x3 x3xn 422x0 3x0 63xo 4 2xi 3xi 63x1 42x； 3x2 63x2 41 / 13 8x 12 0的根；討論用Newton迭代法求它們近似值的收斂階。取初值%2計算根x3的近似值，要求迭代3次。結(jié)果保存4位小數(shù)解： 設(shè) f (x) x3 x2 8x 122f (x)3x22x8f (x)6x2f( 3)0,f (3)0, f (2) 0, f (2) 0, f (2) 10 0如此：3是f(x) 0的單根，故Newton迭代在3附近是平方收斂； 2是f (x) 0的二重根，故Newton迭代在2附近是線性收斂;取 x02 ,

2、 Newton 迭代:xn 1xngf (xn)x3 x2 84 123x2 2x 822xn 3xn 6word2 .設(shè)常數(shù)a 0 ,求出a的取值圍使得解方程組a21x1D2a3x2b213 ax3b3的Jacobi迭代法收斂。解：Jacobi迭代：(k i)(k)xBjxgi0211-203a021Bja 2 03a 1 30iabig a b2a b3迭代矩陣Bj的特征方程:0EBj1 2 a1即：(a)3 14( a) 0特征根：0,14.ia.14譜半徑：(BJ)1時Jacobi迭代收斂1a故：a 44 / 13word23 .設(shè)1用Crout三角分解法求解方程組103Xix2X32

3、用乘幕法求方程組系數(shù)陣的按摸最大的特征值和對應(yīng)的特征向量。解:Crout三角分解:21032101232114LU取v° (0,0,1)T ,計算迭代三次的值2101232,U114Ax bLyUx求解Ly b得I10求解UxTTy得 x 1,1,02v0(0,0,1)T, U0V0max(v0)0,0,1 TV2 Auiv3 Au2V1AU0T2,4,1 ,max(v1)0.5,1, 0.25 T 4U2max(v2)0.5,1,0,8611T9U3max(v3)T0.5,1,0,7306 ,11.443 / 13word4 .試?yán)貌逯刀囗検阶C明：對k 0,1,|,n 2包有等式

4、i 1 (i1)(i證明:- ki0i 1)(i i 1)|(in)i, i1,2,|,nf(x)0,1,|n 2由插值多項式的唯一性，比擬 Lagrange與Newton插值最高項系數(shù)得:ma)rn f,|l|,Xn(Xi X1) (Xi X 1)(X Xi1)“|(x Xn)111由差商與導(dǎo)數(shù)關(guān)系，有fXi,|,Xnf(n1)()(n 1)!1,n將Xii, (i 1,2, |,n), f (x)xk, (k0,1,|n 2）代入上面兩等式，有(i 1) (i i 1)(i i 1) (i n)(i 1) (i i 1)。 i 1) (i n)fX1>|>Xnf(n1)( )

5、0 (n 1)!4 / 13word5 .求4次Hermit插值多項式H (x),滿足:H(0) H (0) 0, H(1) H (1) 1, H (2) 1并寫出誤差表達(dá)式。方 H(1) H (1) 1,H (2) 1a b c 12a 3b 4c 1a 2b 4c 1得 a -, b4122H(x) -x (x 3)4誤差:E(x)解：方法一：因 H(0) H (0) 0,故設(shè)：H(x) x2(a bx cx2),行31一，c24f / 13( ) 22f (x) H(x)x (x 1) (x 2),(0,2)5!方法一：滿足H(0) 0,H (1) H(2) 1的插值多項式為：31 2P

6、2(x)-x -x2 2設(shè)：H(x) p2(x) (A Bx)(x 0)(x 1)(x 2)3H (0) 2B 0,41H (1) (A B) 12得：由 A 1,B344311H (x) -x -x -(x 3)(x 0)(x 1)(x 2)224122x (x 3)4、f ()誤差：E(x) f (x) H(x) (-x2(x 1)2(x 2),(0,2)5!word6.試求求積公式2 f(x)dx Aof (拽)Af (述)的求積系數(shù)Ao, A ,使得其有盡可能高233的代數(shù)精度，是否是Gauss型的？并用此公式計算積分 02 sin xdx結(jié)果保存5位小數(shù)。解： 令f(x) 1,x求積

7、公式準(zhǔn)確成立，有：Ao A 4M穹)A(穹)。得：A。 a 2求積公式:令f(x) x2,x3求積公式準(zhǔn)確成立的，f(x) x4求積公式不是準(zhǔn)確成立的,求積公式代數(shù)精度為3,是Gauss型的;作變換 x (t 2), t 2,28o2sinxdx2)dt2222sin (t28822)dt sin-(t 8 282)0.998486 / 13word7.用最小二乘法求一個形如y ax2 b的經(jīng)驗公式，使它與如下數(shù)據(jù)擬合xi1925313844V解： 取 o(x) 1, 1(x) x2 ,擬合函數(shù)為 y b 0(x) a 1(x) b ax2法方程為：5b 5327a 271.45327a 72

8、77699b 369321.5得： a 0.050351, b 0.9726045擬合函數(shù)為 y 0.050035僅2 0.97260457 / 138.用共腕梯度方法解方程組:2 1 x11 3 x2word取初值x(0)(0,0)T。Po共腕梯度方法：x(k 1) kr(0) bAx(0),(k)kPk,r(k1)(r(k)，產(chǎn))(Pk,Apk)(k)/rkApk(r(k 1),r(k1)(r(k),r(k)Pk ir (k 1)kPk一2 11, 一解： A是對稱正定陣;1 3Po r(0) b Ax(0) (5,5)t(r(0),r(0)(P0, AP0)(1) x(0) x(1) r

9、r(0)10 10(一,一 7755 T0 Ap0(7,7)0P0Pi(2) x(r,r)1(r(0),r(0) 49/4030、r0P0 (,)4949(r,r)(P1, APi)710(1)x 1P1(2,1)Tr(1)r(0)0AP0 (0,0)T解為：x(2,1)t8 / 13word9.應(yīng)用Heun方法:h _ _yn 1 yn(K1 3K2)4Kif(Xn,yn),22K2 f(Xn -h,yn -hKi) 33并在解初值問題5y 8y 0時，問步長h應(yīng)如何選取方能保證方法的絕對穩(wěn)定性?y(0) 2解:yn 11,2中選取數(shù)值穩(wěn)定的步長計算將Heun方法應(yīng)用到方程(1 h)yn,其

10、中h5y5hy(2)的近似值.8y 0上，有:1.6h故取h8,方法是絕對穩(wěn)定的5yn 1y1y21725 yn,1725 yo17一y125342517251.36,345780.9248,25 625(2,0)時，方法是絕對穩(wěn)定的,5(0 , -) (0 ,1.25)時萬法是絕對穩(wěn)定的;451 (0,5) (0,1.25),即 h49 / 13word10.求解常微分方程初值問題yyb的兩步方法:8ynhyn 1yn(5yn 1121求出局部截斷誤差；2討論方法的收斂性；3討論方法的絕對穩(wěn)定性。1125 .8 ,斛.a0 1,a1 0, b 1, b0, b11212（1）把局部截斷誤差T

11、n在4處Taylor展開:TnCoy(Xn) Cihy(Xn)| Crhry(r)(xn)IHcoC4Tn3) IH324n), n(xn , Xn 1 )2C0G 0,方法是相容的;第一特征多項式：（r）(r)0兩根為:01,10,i 1, r11是單根，方法滿足根條件;2-h )r 3h12由收斂的充分必要條件知方法是收斂的。一 、，一一6 一 c（2）穩(wěn)止多項式： （r;h） （1 一h ）r （112由絕對穩(wěn)定性要求知h0,故 1 h12由參考定理知:(r;h)0的兩根ro,i(h)13h5 h12h12_5一 h122 .(1 -h)35(1 h)12h1212/ 5 ,1 h12(

12、6,0),h125 .1 h12即當(dāng)h6,0）時方法是絕對穩(wěn)定的。10 / 13word應(yīng)用1.試確定0是方程f(x) e2x 1 2x 2x2 0的幾重根；取初值比0.25用改良的具有二階收斂速度的Newton迭代法求f (x) 0的根0的近似值。要求迭代2次結(jié)果保存4位小數(shù)。解： f(x) e2x 1 2x 2x2,f (x) 2e2x 2 4xf (x) 4e2x 4f (x) 8e2xf(0) f (0) f (0) 0, f (0) 8 00是方程f (x) 0的3重根;改良的具有二階收斂速度的 Newton迭代法:Xn 1xn3Xn3f(Xn)f (Xn)e24 1 2xn 2x2

13、2e2" 2 4xnXie2x01 2x0 2x02x0 3 2e2% 2 4x0X2e1 1 2x1 2x2Xi 32；2e2x1 2 4x111 / 13word應(yīng)用4.假如用復(fù)化梯形公式計算積分3.e'sinxdx,要求截斷誤差不超過104舍入誤差不計，問需要計算多少個節(jié)點上的函數(shù)值?解:f (x) exsinx, f (x)ex(sinxcosx), f (x) 2ex cosx, f (x) 3ex(cosx sin x)復(fù)化求積公式余項為:En(f)(b a)其中:12h J2h f (),a,b假如517.5cosx 1,有 |f ()2e3En(f)h 3.86518,104,得:10 3h23 10 4故至少需519個節(jié)點才能保證截斷誤差不超過10 4。12 / 13word應(yīng)用9.寫出經(jīng)典4階Runge-Kutta方法求解初值問題y 8 3yy(0) 2的計算公式，并取步長h 0.2 ,計算y(0.4)的近似值.小數(shù)點后至少保存4位解：f(x