1、非線性方程 二分法迭代,科學(xué)計(jì)算方法2, ,GPS (Global Positioning System),Reference: Timothy Sauer, Numerical Analysis,非線性科學(xué)是當(dāng)今科學(xué)發(fā)展的一個(gè)重要研究方向,其中非線性方程的求根是一個(gè)不可或缺的研究?jī)?nèi)容。然后非線性方程的求根非常困難。,高于次的代數(shù)方程,不存在通用的求根公式,超越方程一般很難直接求出其準(zhǔn)確解,例:,求 f (x)= 0 根的問題包括如下方面:,根的存在性,即f(x)=0有沒有根? 若有,有 幾個(gè)根?,根的范圍,確定有根區(qū)間,根的精度,已知一個(gè)根的近似值后,能否 將它精確到足夠精度?,Where

2、is Flight MH370?,定理2(零點(diǎn)定理) 若函數(shù)f在閉區(qū)間a, b上連續(xù),且f(a)與f(b)異號(hào), 則至少存在一點(diǎn) x0屬于(a, b),使得f(x0)=0,即方程f(x)=0在(a, b)內(nèi)至少有一個(gè)根。,零點(diǎn)定理回答了根的存在性和根的范圍的 問題,如何去尋找足夠精度的根？,游戲:模仿幸運(yùn)52請(qǐng)同學(xué)們猜測(cè)這部電腦的價(jià)格？,化大為小,a0,b0,a1,b1=a0,x0 a1,b1=x0,b0,已知f(x)=0在a0,b0內(nèi)有一根,且f(a0)f(b0)0,(2)判斷若f(xn-1) f(an-1)0,則anan-1,bn xn-1 否則 anxn-1, bnbn-1,x0 x1

3、,二分法迭代,(1)計(jì)算 xn-1=0.5(an-1+bn-1)和f(xn-1) 判斷若f(xn-1)=0,則xn-1是根,否則轉(zhuǎn)下一步;,例 1 二分法求方程 在區(qū)間 0, 1內(nèi)的根。,函數(shù)在0,1內(nèi)有唯一零點(diǎn), 故0,1是隔根區(qū)間。,解: 令,x4=0.5(0.4375+0.5)=0.4688作為近似值, 誤差不超過1/(25) 。 精確值0.443573534104293。,二分法迭代將得到一系列區(qū)間序列,定理2.2 設(shè)x*是 f(x)=0在a0, b0內(nèi)的唯一根,且 f(a 0)f(b 0)0,則二分過程中各區(qū)間的中點(diǎn)數(shù)列,性質(zhì)1. f(an)f(bn)0; 2. bn an = (b

4、0 a0)/ 2n,滿足 | xn x*| (b 0 a 0)/ 2n+1,注記: 若要,思想簡(jiǎn)單, 僅需函數(shù)值符號(hào), 容易實(shí)現(xiàn)。,二分法的反思:,收斂速度慢和不能求復(fù)根等不足。,function root = bisect(func,x1,x2,n) for i = 1:n x3 = 0.5*(x1 + x2); f3 = feval(func,x3); if f3 = 0 root = x3; return end if f2*f3 0 x1 = x3; f1 = f3; else x2 = x3; f2 = f3; end end,練習(xí) 1,試試實(shí)現(xiàn)二分法？,標(biāo)準(zhǔn): 1. 代碼美 2. 文檔美,作業(yè) 1,試試設(shè)計(jì)新的求根方法？,標(biāo)準(zhǔn): 文檔(word, ppt 或 pdf)和代碼,提示: 黃金分割方法(華羅庚),設(shè)f(x) = 0的根為 x*,通過迭代計(jì)算,產(chǎn)生序列: x0 x1 x2 xn,只須,迭代的思想,迭代法是一種逐次逼近法, 這種方法使用某個(gè)固定公式反復(fù)校正根的近似值, 使之逐步精確化, 最后得到滿足精度要求的結(jié)果。,蘋果成功的原因,Comes out with industry-sha