迭代法的設(shè)計和運行

1、數(shù)值實驗報告：迭代法的設(shè)計和運行姓 名：孫 闖學(xué) 號：P200909071實驗報告數(shù)值實驗1.1 迭代法的設(shè)計和運行 姓名：孫闖 學(xué)號：P200909071 一、實驗?zāi)康牡ㄊ墙飧鞣N方程(含方程組)的基本方法, 它常通過構(gòu)造一個定常的迭代格式，重復(fù)計算而產(chǎn)生一個收斂的解的序列（標量序列，向量序列等），逐步逼近問題的真解。對同一個問題，?？稍O(shè)計不同的迭代格式，這些格式的計算效果可能相差很大，對初始值的選擇的要求不同，收斂速度也不同。二、實驗題求方程在1，1.5內(nèi)的根, 設(shè)計下列3個迭代格式：(1) ； (2) ；(3) 三、實驗步驟和要求 (1) 選取若干個可行的初始值, 分別用3種迭代格式進

2、行計算, 記錄每步的計算結(jié)果, 并在坐標紙上畫出, 分析其變化趨勢,判別其收斂性;(2) 觀察和分析迭代法結(jié)果關(guān)于初始值的依賴關(guān)系; (3) 你還能設(shè)計出新的迭代格式嗎?四、實驗結(jié)果（1）取初始值x0=1.10，迭代結(jié)果如下迭代步數(shù)n迭代格式1 4.9289999999999991.4721582795338281.4002800840280102-2.020014210889998e+0021.3047475779514131.36079234494813138.079171661661471e+0061.3945287527281701.3657949658048234-5.2735215

3、20393093e+0201.3498190038521361.36515814053526051.466567876063306e+0621.3730089642613391.3652391578493736-3.154325477210480e+1861.3612188643181131.36522884997455271.3672760056485991.36523016143783981.36418057074659291.365766768990587101.364955079617933111.365370733269326121.365157960973024131.365266

4、899027172141.365211128768512151.365239681452983161.365225063708751171.365232547459919發(fā)散收斂收斂（2）取初始值x0=1.25，迭代結(jié)果如下迭代步數(shù)n迭代格式123456789101112131415163.046875000000000-52.3724174499511721.326368855327903e+005-2.333490535726542e+0151.270627151578478e+046-2.051419097135232e+138Inf1.4183507147387771.33666350

5、22588991.3794767424336421.3578398600116361.368987870356384 1.3632994708712001.3662166229369891.3647244488432561.3654887225372861.3650975328383201.3652978300701201.3651952917475411.3652477889942691.3652209128569421.3652346725077521.3652276281284381.3801311186847091.3633380946752001.3654707847825841.3

6、651993812478311.3652339107429841.365229517559543發(fā)散收斂收斂（3）取初始值x0=1.30，迭代結(jié)果如下：迭代步數(shù)n迭代格式123456789101112131415 2.343000000000000-22.4778436070000019.323516446902724e+003-8.108219573345890e+011 5.330605009160471e+035-1.514710057873552e+107 Inf1.3966925216381741.348647661864871 1.373591233274327 1.3609163

7、08043531 1.367429720635338 1.364101566687551 1.365807135043137 1.364934392575762 1.365381318675585 1.365152540199034 1.365269673871691 1.3652097080555861.365240408776732 1.3652246913402981.3652327380963381.373605639486890 1.3641656335521121.365365453979880 1.365212781718142 1.3652322057977511.365229

8、734478521發(fā)散收斂收斂 （4）取初始值x0=1.35，迭代結(jié)果如下：迭代步數(shù)n迭代格式123456789101112131.599624999999998 -2.728696237447233 -2.194548379150916 -0.889680276517144 6.648405259779069 -4.540248702677044e+0029.276704496659131e+007 -7.983276783143062e+023 5.087958496701220e+071 -1.317136182179706e+2151.372918879613796 1.3612656

9、44820378 1.367252230974365 1.364192788092528 1.365760526179413 1.364958278834917 1.365369096217715 1.365158799295786 1.365266469894964 1.365211348482401 1.365239568971871 1.3652251212957411.3652325179778781.367171854049327 1.364983021234603 1.365261439166450 1.3652260151519791.365230522110880發(fā)散收斂收斂（

10、5）取初始值x0=1.40，迭代結(jié)果如下：迭代步數(shù)n迭代格式12345678910111213140.8160000000000027.609237503999985 -6.545713651042345e+002 2.787455514286298e+008 -2.165827389138137e+025 1.015948104902566e+076 -1.048611397131342e+2281.346848172586651 1.374483304275448 1.360452145192081 1.367665374147659 1.363980405747051 1.3658690

11、29182214 1.364902669703743 1.365397550251570 1.365144227823156 1.365273928846953 1.365207529505614 1.365241524066006 1.3652241203438951.3652330304213571.360827634879543 1.365790470333855 1.365158712401375 1.365239085089733 1.3652288592316941.365230160260058發(fā)散收斂收斂 五、實驗分析選取迭代初值為1.35，迭代結(jié)果如圖1所示。由圖可明顯知道，

12、第一種迭代方式發(fā)散，第二種和第三種收斂。但是第二種迭代方式迭代的步數(shù)多，因為其迭代所得的值是忽大忽小的，所以迭代所需的步數(shù)多。第三種迭代方式步數(shù)最少。所以選擇合適的迭代方式非常重要。圖1 初始值為1.35的迭代結(jié)果圖取不同的迭代初值，迭代的步數(shù)是不一樣的。如果取得初值越接近方程的根，則迭代所需的步數(shù)越少。比如本題方程的根是1.36523，取初始值1.35迭代時，第二種和第三種收斂的迭代格式步數(shù)均是最少的。但是對于第一種發(fā)散的迭代格式，不管迭代初值如何，仍然是發(fā)散的。該方程的根為1.36523。迭代誤差小于迭代初值迭代步數(shù)1.101771.251661.301561.351351.40146六、

13、實驗拓展1設(shè)計新的迭代格式和，取迭代初值為x0=1.40迭代步數(shù)n迭代格式11.2926608140191311.35816618911174821.4912246256274611.36669133333052431.0338407095517601.36492882562862141.7888918319131611.36529213786356050.704774834213240+ 1.220705820753264i1.36521720135122762.468416038638107- 0.379515051007131i1.36523265576040771.3652294684

14、63052發(fā)散收斂2用Newton法（切線法）設(shè)計迭代格式：由于f(x)=x3+4x2-10，所以f(x)=3x2+8x,根據(jù)Newton法迭代 得到新的迭代格式為 取迭代初值x0=1.10，結(jié)果如下迭代步數(shù)n迭代格式11.40804505229284021.36610134082686531.36523038537962241.36523001341416556收斂，且且收斂速度較快在1，1.5內(nèi)取不同迭代初值，Newton迭代法的迭代步數(shù)相差不大，而且收斂速度都很快。迭代初值迭代步數(shù)1.1041.2541.3031.3531.403附錄：MATLAB程序如下主程序：function A=d

15、iedai(x,max)%輸入量-x0是初始值，k是迭代步數(shù)；輸出量-矩陣A，包含迭代步數(shù)、偏差和迭代值xkformat longA=zeros(max,3)；%預(yù)置MAX*3零矩陣k=1;x(k)=x; piancha=abs(x(k)-x);A(k,:)=(k-1) piancha x;while (k<=max) x(k+1)=fun4(x(k)；%程序中調(diào)用迭代函數(shù)y=(x) piancha= abs(x(k+1)-x(k); k=k+1; xk=x(k); A(k,:)=(k-1) piancha xk; if (piancha <10(-5) disp(祝賀你！此迭代格式收斂，且收斂速度較快) break endendif(piancha >=10(-5) disp(請用戶注意：此迭代格式發(fā)散，請