二分法上機(jī)實(shí)驗報告.doc

1中國礦業(yè)大學(xué)(北京)理學(xué)院數(shù)值分析實(shí)驗報告實(shí)驗名稱 二分法求方程的近似解 實(shí)驗時間 2013 年 3 月 13 日組長簽名 1110720011011107200101111072001021110720010311107200119組員簽名班級 信息與計算科學(xué) 一班 學(xué)號11107200120成績一、實(shí)驗?zāi)康?#xff0c;內(nèi)容 二、相關(guān)背景知識介紹 三、代碼四、數(shù)值結(jié)果 五、計算結(jié)果的分析 六、計算中出現(xiàn)的問題，解決方法及體會一、實(shí)驗?zāi)康?#xff0c;內(nèi)容通過用二分法求一個具體的方程 ， 精確到 的近似根，掌握二分法，會編寫相關(guān)代碼。02x1,502、相關(guān)背景知識介紹(1)算法原理或計算公式 二分法的基本思路是通過計算隔根區(qū)間的中點(diǎn)，逐步將隔根區(qū)間縮小，從而可得方程的近似根數(shù)列 。nx1*2kab(2)程序設(shè)計思路二分法原理不妨設(shè)方程 在隔根區(qū)間 上 .0xf ,ba0,bff1將區(qū)間 二分得中點(diǎn) ,計算 在中點(diǎn)的函數(shù)值 ，若 ，則 就是方程的根；,ba2x2baf02baf 2bax否則,若 ,由于 在 內(nèi)不變號，所以隔根區(qū)間變?yōu)?.同理，若 , 隔根區(qū)間變?yōu)?.將新區(qū)間記為2fxf , 0f 2,ba.1ba2將 重復(fù)上述步驟 1.得到一系列隔根區(qū)間1,kb,aa并有 ，且后一區(qū)間的長度都是前一區(qū)間長度的一半，所以 的長度為 ，區(qū)間 的長度趨向0f*kkxb kb,a ,2kabk kb,a于零，即這些區(qū)間最終收縮于一點(diǎn) ， 顯然 就是方程 的根。x0xf三、代碼(Matlab)cleara0;b1;eps10-5;k0;y a - 2-a;while absa-bepskk1xab/2if a-2-a02xbreakelseif y*x-pow2-x0bx;else ax;endendendxfx-pow2-xK四、數(shù)值結(jié)果k 1 2 3 4 5 6 7 8 9x 0.5000 0.7500 0.6250 0.6875 0.6563 0.6406 0.6484 0.6445 0.6426k 10 11 12 13 14 15 16 17x 0.6416 0.6411 0.6414 0.6412 0.6412 0.6412 0.6412 0.6412f -5.4870e-006k 17五、計算結(jié)果的分析這個方程運(yùn)用二分法求得四位有效數(shù)字的解為 x0.6414。圖一6、計算中出現(xiàn)的問題，解決方法及體會1.剛開始設(shè)計程序，出現(xiàn)死循環(huán)，無法得出結(jié)果，要經(jīng)過不斷調(diào)試。2.不熟悉 Matlab 的用法，導(dǎo)致解決問題效率下降，要學(xué)好相關(guān)基礎(chǔ)知識。3.由圖一可以看出二分在區(qū)間很大時，收斂接近根的速度很快，而當(dāng)區(qū)間較小時，靠近要求精度根速度變得相當(dāng)緩慢，計算量也很大，對滿足高精度的要求比較困難。我們想象此法可以和別的后半段接近根較快的方法結(jié)合，先求出適用于別的方法的初值，應(yīng)當(dāng)會優(yōu)化算法。3教 師 評 語指導(dǎo)教師 年 月 日

總結(jié)

以上是生活随笔為你收集整理的计算机科学导论二分法实验报告,二分法上机实验报告.doc的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。