第二章 非線性方程

1、第二章 非線性方程求根 代數(shù)方程求根問題是一個古老的數(shù)學問題，早在16世紀時找到了三次、四次方程的求根公式。但直到19世紀才證明 次的一般代數(shù)方程式不能用代數(shù)公式求解。因此需要研究用數(shù)值方法求得滿足一定精度的代數(shù)方程式的近似解。 在工程和科學技術中許多問題常常歸結為求解非線性方程式問題，例如在控制系統(tǒng)的設計領域，研究人口增長率等。 本章將介紹幾種求解非線性方程近似解的數(shù)值方法。最后一節(jié)簡單介紹求解非線性方程組的方法。 §1 二 分 法 設有非線性方程 (1.1) 其中，為上連續(xù)函數(shù)且設（不妨設方程（1.1）于內(nèi)僅有一個實根）。 求方程（1.1）實根的二分法過程，就是將含根區(qū)間逐步分半

2、，檢查函數(shù)值符號的變化，以便確定含根的充分小區(qū)間。 二分法敘述如下：記。第1步分半計算：將分半，計算中點及，如果 則根一定在區(qū)間內(nèi)，否則根一定在區(qū)間內(nèi)（若 則）。于是得到長度縮小一半的含根區(qū)間，即 ，且第步分半計算：重復上述過程，設已完成第1步，第步分半計算得到含根區(qū)間 且滿足： （1）， 即； （2）；現(xiàn)進行第k步分半計算： （3） 計算 且有 （1.2）（4） 確定新的含根區(qū)間，即如果，則根一定在內(nèi)，否則根一定在區(qū)間內(nèi)，且有 總之，由上述二分法得到一序列，由（1.2），則有 可用二分法求方程實根的近似值到任意指定的精度。事實上，設為給定精度要求，試確定分半次數(shù)k使 由，兩邊取對數(shù)，即得 （

3、1.3）例1 用二分法求于內(nèi)一個實根，且要求精確到小數(shù)后第3位（即要求）。 解 由和公式（1.3）可確定所需分半次數(shù)。計算結果如下表。 表 6-1k11.02.01.58.89062521.01.51.251.56469731.01.251.125-0.09771341.1251.251.18750.61665351.1251.18751.156250.23326961.1251.156251.1406250.061577871.1251.1406251.13281319575681.1328131.1406251.1367190.020619091.1328131.1367191.13476

4、64.307101.1328131.1347661.133789111.1337891.1347661.134277二分法優(yōu)點是方法簡單，且對只要求連續(xù)即可。可用二分法求出于內(nèi)全部實根。但二分法不能求復根及偶數(shù)重根。二分法：設有方程，其中于連續(xù)，且滿足條件（且設于內(nèi)只有一個實根）。 （1） 計算，； （2） 如果 或 ,則輸出 ； （3） 如果，則；否則其中表示給定的最大分半次數(shù)，當 或時分半終止。 §2 迭 代 法迭代法是一種逐次逼近法。它是求解代數(shù)方程，超越方程及方程組的一種基本方法，但存在收斂性及收斂快慢問題。 為了用迭代法求非線性方程的近似根，首先需要將此方程轉(zhuǎn)化為等價的方程

5、 (2.1)顯然，將轉(zhuǎn)化為等價方程(2.1)的方法是很多的。 例2 方程可用不同方法轉(zhuǎn)化為等價方程 ( a) (b) 定義：（迭代法）設方程為。 （1）選取方程根的一個初始近似，且按下述逐次代入法，構造一近似解序列： （2.2）這種方法稱為迭代法（或稱為單點迭代法）。稱為迭代函數(shù)。 （2）如果由迭代法產(chǎn)生的序列有極限存在，即，則稱為收斂或稱迭代過程（2.2）收斂。否則稱不收斂。 設為連續(xù)函數(shù)，且 ，則有，即為方程（2.1）的解（稱為函數(shù)的不動點）。 事實上，由迭代過程（2.2）兩邊取極限，則有 顯然，在由方程轉(zhuǎn)化為等價的方程時，選擇不同的迭代函數(shù)，就會產(chǎn)生不同的序列（即使初始值選擇一樣），且這

6、些序列的收斂情況也不會相同。 例3 對例2中方程，考察用迭代法求根： （a）， (b) ， 表 6 2 k01.01.011.3414710.52359921.4738200.02360131.49530141.49715251.49728961.49730071.497300 由計算看出，我們選取的兩個迭代函數(shù)和，分別構造序列 收斂情況不一樣（初始值都取為1.0），在（a）種情況收斂且。在（b）種情況出現(xiàn)計算無定義。 因此，對于用迭代法求方程近似根需要研究下述問題： （1） 如何選取迭代函數(shù)使迭代過程 收斂。 （2） 若 收斂較慢時，怎樣加速收斂。 迭代法的幾何意義： 從幾何上解釋，求方程根

7、的問題，是求曲線與直線交點的橫坐標。當?shù)瘮?shù)的導數(shù)在根處滿足下述幾種條件時，從幾何上來考查迭代過程的收斂情況如圖(1)-(4)。從曲線上一點出發(fā)，沿著平行于軸方向前進交于一點，再從點沿平行于軸方向前進交于點，顯然，的橫坐標就是。繼續(xù)這過程就得到序列，且從幾何上觀察知道在（1），（2）情況下收斂于，在（3）、（4）情況不收斂于。 由迭代法的幾何意義可知，為了保證迭代過程收斂，應該要求當時，迭代函數(shù)的導數(shù)滿足條件；否則方程于可能有幾個根或迭代法不收斂，為此有下述關于迭代法收斂性定理。 定理1 設有方程 （1） 設于一階導數(shù)存在； （2） 當時，有； （3） 當時，滿足條件：。 則 （1） 在上有

8、唯一解； （2） 對任意選取初始值，迭代過程收斂，即；（3）； （2.3）（4）誤差估計 （2.4）只證明（2），（3），（4）。 證（2）由定理假設條件（2），當取時，則有，。記誤差，由中值公式有 其中c在與之間，即。又利用假設條件（3）得到誤差的遞推關系 （2.5）反復利用（3.5），得到 （當）即 。證（3） 由迭代公式，顯然有 （2.6）其中 在與之間，于是 即 證（4） 反復利用（2.6），可得 。 由定理3結果（2.3）可知，當計算得到的相鄰兩次迭代滿足條件 （2.7）時，則誤差 。所以在電算時可利用來控制迭代終止，但是要注意，當時，即使很小，但誤差還可能較大。當已知，及給定精度要

9、求時，利用（2.4）可確定使誤差達到給定精度要求所需要迭代次數(shù)。事實上，由可得 （2.8）定理1中的假設條件：當時，。在一般情況下，可能對于大范圍的含根區(qū)間不滿足，而在根的鄰域是成立的，為此有下述迭代過程局部收斂性結果。定理 2（迭代法的局部收斂性）設給定方程 （1）設為方程的解； （2）設在的鄰近連續(xù)可微且有則對任意取初值，迭代過程（）收斂于（稱迭代過程具有局部收斂性）。 證： 取，于是只要驗證定理3中條件（2）成立，定理4即得證。 事實上，設，則滿足 ，其中，。這說明。 例4 試用迭代法解方程：。 解 （1）由于當時，連續(xù)，且顯然有 及 即知，方程于及內(nèi)有根，分別記為及。 （2）考查取初值

10、時，迭代過程的收斂性，其中迭代函數(shù)為。 顯然，及為增函數(shù)，則當時，。又由，則當時，于是由定理3可知，當初值時迭代過程收斂。結果見表63表63 00.010.6931471820.99071046141.1461931151.1461932如果要求近似根準確到小數(shù)后第6位（即要求），則由上表可知，且。所以。于是 ，。（3）為了求內(nèi)方程的根，考查迭代過程 （2.9） 顯然，當時，。所以，迭代過程（2.9）不 收斂于。 （4）可將方程轉(zhuǎn)化等價方程 ，從而，且有。因為當時，所以，當選取初值時，迭代過程 收斂。如取，則迭代12次有 ，且有。由上例可見，對于方程，迭代函數(shù)選取不同，相應由迭代法產(chǎn)生的收斂情

11、況也不一樣。因此，我們應該選擇迭代函數(shù)，使構造的迭代過程 收斂且較快。迭代法：求解方程 （1） 選取解的初始估計；給定精度要求（2）對于，計算，其中為給定的最大迭代次數(shù)。（3）終止條件：當時（或 ）迭代終止。 §3 牛 頓 法 解非線性方程牛頓方法是一種將非線性函數(shù)線性化的方法。牛頓方法的最大優(yōu)點是在方程單根附近具有較高的收斂速度，牛頓方法可用來計算的實根，還可計算代數(shù)方程的復根。它是求解方程應用最廣泛的迭代法之一。 3.1 牛頓法公式及誤差分析 設有非線性方程 （3.1） 其中，假設在上一階連續(xù)可微，且 ；又設是的一個零點的近似值（設）。現(xiàn)考慮用過曲線上點 的切線近似代替函數(shù)，即用

12、線性函數(shù)代替。且用切線的零點，作為方程（3.1）根的近似值，即 （3.2）一般，若已求得，將（3.2）中換為，重復上述過程，即得求方程根的牛頓方法的計算公式 （3.3）下面利用的泰勒公式進行誤差分析。設已知的根的第次近似，于是在點泰勒公式為（設二次連續(xù)可微）： （3.4）其中C在與之間。 如果用線性函數(shù) 近似代替，其誤差為。且用根作為的根的近似值，又得到牛頓公式.現(xiàn)在（3.4）中取，則有 于是 .利用牛頓公式（3.3）即得誤差關系式 （3.5）誤差公式（3.5）說明的誤差是與的誤差的平方成比例的。當初始誤差充分小時，以后迭代的誤差將非?？斓臏p少。 由計算公式（3.3）可知，用牛頓方法求方程根，

13、每計算一步需要計算一次函數(shù)值及一次導數(shù)。 例5 用牛頓法求 根。 解 顯然，方程在，內(nèi)各有一根。牛頓法計算公式為 以為初值求解如下：，由此可見，經(jīng)過3次迭代，即得誤差小于的平方根。同理可得另一平方根-5。注：當初值選取靠近根時，牛頓法收斂且收斂較快；當初值不是選取接近方程根時，牛頓法可能會給出發(fā)散的結果。 3.2 牛頓法的局部收斂性 設有方程，由于牛頓法是一種迭代法，其中迭代函數(shù)為 可用迭代法理論來考查牛頓方法的收斂性。 定理3 設方程有根，而在根鄰近具有連續(xù)二階導數(shù)，且，則存在的一個鄰域 ，使得對于任意初值， 由牛頓法產(chǎn)生的序列收 斂 于，且有 （3.6）證 由于牛頓法是一個迭代法， 其迭代函數(shù)為 故 從而 于是由定理2知，牛頓法迭代法為局部收斂。且由（3.5）取極限即得（3.6）。 注： 在用牛頓法求的單根時，一般可用來估計的誤差，即當充分接近時，若 ，意味著。 §4 弦 割 法 如果函數(shù)比較復雜，求導可能有困難，這時可將牛頓公式中近似用差商來代