不同方案收斂速度的比較

1、數(shù)值分析實驗報告實驗名稱不同方案收斂速度的比較實驗時間 13年10 月 9日姓名班級學(xué)號成績一、 實驗?zāi)康?通過實驗體會數(shù)值計算中算法選擇的重要地位。2了解數(shù)值計算方法的必要性。3體會數(shù)值計算的收斂性與收斂速度。二、 實驗內(nèi)容 三種求2的算法比較三、算法描述 方案一：利用級數(shù)， 設(shè) ，則方案二：對上述，按生成新數(shù)列，則方案三：利用級數(shù) 設(shè)，則四、程序流程圖由于實驗方案明顯、簡單，實現(xiàn)步驟及流程圖省略。五、實驗結(jié)果計算結(jié)果如下： 運算可知 ln2= 由方案一利用matlab編輯，可以得到運算結(jié)果為： S = k = 100001由方案二利用matlab編輯，可以得到運算結(jié)果為： S = n =3

2、1由方案三可以得到運算結(jié)果為： k =15六、實驗結(jié)果分析 由方案一可以知道當(dāng)k=100001時，S和ln2的誤差計算可知：|S-ln2|= 4.999974963726395e-06 ,是在誤差允許的范圍內(nèi)得到的一個相對準(zhǔn)確的結(jié)果，S = ，通過循環(huán)得到的結(jié)果可知，k=100001，相對循環(huán)的次數(shù)很多，說明它收斂的速度很慢，相對較耗費時間，工程量比較大！ 由方案二可知n=31時，S和ln2的誤差計算可知：|S-ln2|= 4.862077586875735e-06，也是在誤差允許的范圍內(nèi)通過運用方案二得到的計算結(jié)果，S =,由這個循環(huán)得到的結(jié)果可知，n=31，由于是從3開始的，也即是循環(huán)的次

3、數(shù)為（31-3）=28，較方案一來說，運算次數(shù)大量減少，說明它的收斂速度較快，耗時很少，便于平常的運算和實驗操作！ 由方案三可知n=15時，S和ln2的誤差計算可知：|S-ln2|= 3.840474466665356e-06，也是在誤差允許的范圍內(nèi)得到的結(jié)果，S=，這個循環(huán)的次數(shù)為方案一的一半左右，而且誤差也是相對方案一和二來說比較小，次數(shù)為15反映出它的收斂速度很快，有可操作性！ 綜合比較上面三種方案，針對方案一，對于漸進式的選擇和舍入誤差使得收斂速度變慢，運算量增加；而方案二，每次運算的Sn，Sn-1，Sn-2都會在計算機內(nèi)自己舍入誤差，使得運算的結(jié)果都不是在精確值的條件下運行的，由于這

4、些誤差會影響收斂速度；對于方案三，由于漸進式的選擇，使得運算的結(jié)果較為接近ln2，而且這個運算式的舍入過程很少，所以由于舍入誤差引起的誤差也較少。綜上所述可以發(fā)現(xiàn)，一個漸進式即算法的選擇對于一個計算過程很重要，誤差估計和遞推計算的穩(wěn)定性都會直接影響運算結(jié)果。選擇好的算法不僅可以在時間上，也可以在復(fù)雜程度上縮小工作量，進而得到一個相對準(zhǔn)確的數(shù)值。教 師 評 語指導(dǎo)教師： 年 月 日數(shù)值分析上機實驗原始記錄實驗名稱：不同方案收斂速度的比較 實驗時間： 2013 年 10 月 9 日姓名： 學(xué)號： 班級： %方案一S=0;s=0;k=1; while abs(S-log(2)>=1/2*10e

5、-6 s=(-1)(k-1)/k; S=S+s; k=k+1; endS kS =k = 100001>> %方案二S=0;s=0;n=3;while abs(S-log(2)>=0.5*10(-5)A=0;B=0;C=0; for i=1:n s=(-1)(i-1)/i; A=A+s; i=i+1; end s; A; for j=1:(n-1) s=(-1)(j-1)/j; B=B+s; j=j+1; end s;B; for k=1:(n-2) s=(-1)(k-1)/k; C=C+s; k=k+1; end s; C;S=A-(A-B)2/(A-2*B+C);n=n+1;endS;n;>> >> SS =>> nn = 31%方案三S=0;s=0;k=1; whil