這篇文章主要介紹matlab如何實現(xiàn)二分法方程求根,文中介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們一定要看完!
創(chuàng)新互聯(lián)是一家以網(wǎng)站設計建設,小程序開發(fā)、網(wǎng)站開發(fā)設計,網(wǎng)絡軟件產品開發(fā),企業(yè)互聯(lián)網(wǎng)推廣服務為主的民營科技公司。主要業(yè)務涵蓋:為客戶提供網(wǎng)站策劃、網(wǎng)站設計、網(wǎng)站開發(fā)、主機域名、網(wǎng)站優(yōu)化排名、賣友情鏈接等服務領域。憑借建站老客戶口碑做市場,建設網(wǎng)站時,根據(jù)市場搜索規(guī)律和搜索引擎的排名收錄規(guī)律編程,全力為建站客戶設計制作排名好的網(wǎng)站,深受老客戶認可和贊譽。
假設函數(shù) 在區(qū)間 上連續(xù),并且 ,此時就可以用二分法求解。
求解偽代碼:
a1 = a;
b1 = b;
計算中點
如果 ,那么方程的解為 ,終止
如果
如果 ,
如果 ,
重復上述步驟2到4,直到滿足誤差,停止迭代。
, 是第n次迭代的結果,p為真實解。
方程 在區(qū)間【0,2】上的近似解。
;
所以區(qū)間縮放到【1,2】,
利用matlab計算如下:
% -------------- inputs -------------------f = @(x) 3*x^2-x-3;a = 0;b = 2;% tolerance / max iterTOL = 1e-4; NI = 50;% -------------------------------------------------------% STEP 1: initializationi = 1;fa = f(a);converge = false; % convergence flag% STEP 2: iterationwhile i<=NI% STEP 3: compute p at the i's stepp = a+(b-a)/2;fp = f(p);% STEP 4: check if meets the stopping criteriaif (abs(fp)0a = p; fa = fp;elseb = p;endendend
以上是“matlab如何實現(xiàn)二分法方程求根”這篇文章的所有內容,感謝各位的閱讀!希望分享的內容對大家有幫助,更多相關知識,歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道!