第9章 2-非線性方程組的迭代解法

中至少有一個是x的非線性函數(shù)，若其全為線性的則為線性方程組。9.3非線性方程組的迭代解法

含有n個未知數(shù)的n個方程的非線性方程組為

(1)其中為n維列向量，非線性方程組包括：

高次方程組,即代數(shù)方程組超越方程組產(chǎn)生背景:許多科學理論與工程技術都可化為非線性方程組求解的特點：無求解公式，無直接解法，

難求得精確解。求解的方法：間接法即迭代法。迭代法求解的要求：

收斂計算效率(快慢)

數(shù)值穩(wěn)定性(考慮計算機的舍入誤差)初始值好迭代公式合適(好的)一.一般迭代法1.建立

xyy=xxyy=xxyy=xxyy=xx*x*x*x*x0p0x1p1x0p0x1p1x0p0x1p1x0p0x1p12.迭代法的收斂性

(局部收斂性)（收斂定理）證:

1o由微分中值定理

2o~4o注：L越小，收斂越快。定理3.1指出:只要構造的迭代函數(shù)滿足迭代法的收斂階(收斂速度)

定義3.1.:設若有實數(shù)c>0,p≥1,使

則稱是p階收斂,相應的迭代法稱為p階方法.特別,p=1,稱線性收斂;1<p<2,稱超線性收斂

p=2,稱平方收斂。(2)3.非線性方程組的一般迭代法(3)并構造不動點迭代法

(4)把方程組(1)改寫成下面便于迭代的等價形式：對于給定的初始點,若由此生成的序列收斂，(1)的解。是方程組的不動點，是迭。即滿足連續(xù)函數(shù).則的是自變量是連續(xù)的,即.且代函數(shù)例1設有非線性方程組把它寫成等價形式

并由此構造不動點迭代法

取初始點。計算結果列于表1，可見迭代收斂到方程的解表1k012…

18

1900.80.928…0.9999999720.99999998900.80.931…0.9999999720.999999989函數(shù)也稱映射，若函數(shù)的定義域為，則可用映射符號簡便地表示為。為了討論不動點迭代法（4）的收斂性，先定義向量值函數(shù)的映內性和壓縮性。二.牛頓迭代法1.一元方程牛頓法

將在點作Taylor展開:

——Taylor展開線性化（重要思想）近似于解出,

記為,則Newton迭代法的幾何意義

與軸的交點，作為下一個迭代點,即用在點處的切線Newton迭代法需要求每個迭代點處的導數(shù)復雜！(6)這種格式稱為簡化Newton迭代法精度稍低Newton法的改進無論哪種迭代法：Newton迭代法簡化Newton法用Newton迭代法求解:x0=2x1=-3.54x2=13.95x3=-279.34x4=122017是否收斂均與初值的位置有關.例3

x0=1x1=-0.5708x2=0.1169x3=-0.0011x4=7.963110-10x5=0收斂發(fā)散迭代法的局部收斂性2.非線性方程組的Newton法對于非線性方程組，也可以構造類似于一元方程的Newton迭代法。設是方程組（1）的解，是方程組的一個近似解。用點處的一階Taylor展開式近似每一個分量函數(shù)的值

，有其中為的Jacobi矩陣在的值，而寫成向量形式有若矩陣非奇異，則可以用使（7）右端為零的向量作為一個新的的近似值，記為，于是得到Newton迭代法（8)其中是給定的初值向量。如果寫成一般不動點迭代的形式，則Newton迭代函數(shù)為(9)在Newton法實際計算過程中，第k步是先解線性方程組解出后，再令，其中包括了計算向量和矩陣

(10)例4用Newton法解例1的方程組解

對該方程組有取初始向量，解方程組，即求出后，。同理計算結果列于表2。可見，Newton法的收斂速度比例1中的迭代法（4）要快的多。表2k01