第四章 非線性方程求根--高等工程數(shù)學(xué)_第1頁
第四章 非線性方程求根--高等工程數(shù)學(xué)_第2頁
第四章 非線性方程求根--高等工程數(shù)學(xué)_第3頁
第四章 非線性方程求根--高等工程數(shù)學(xué)_第4頁
第四章 非線性方程求根--高等工程數(shù)學(xué)_第5頁
已閱讀5頁,還剩73頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

1、高等工程數(shù)學(xué)高等工程數(shù)學(xué)教學(xué)要求教學(xué)要求 1 1)教材:高等工程數(shù)學(xué),鄭洲順編)教材:高等工程數(shù)學(xué),鄭洲順編2 2)考核:平時(shí))考核:平時(shí)30%+30%+期末考試期末考試70%70%平時(shí)平時(shí)30%30%:平時(shí)登記平時(shí)登記+ +每章論文報(bào)告。每章論文報(bào)告。論文報(bào)告要求:論文報(bào)告要求:以本章介紹的數(shù)學(xué)方法為基礎(chǔ),自主選擇有一定以本章介紹的數(shù)學(xué)方法為基礎(chǔ),自主選擇有一定實(shí)際意義、背景知識大眾化的實(shí)際問題,進(jìn)行計(jì)實(shí)際意義、背景知識大眾化的實(shí)際問題,進(jìn)行計(jì)算分析,寫出學(xué)習(xí)論文。包括:(算分析,寫出學(xué)習(xí)論文。包括:(1 1)問題描述,)問題描述,(2 2)建立模型(包括必要的假設(shè)),()建立模型(包括必要

2、的假設(shè)),(3 3)算法)算法選擇(必須是本章的算法)和設(shè)計(jì),(選擇(必須是本章的算法)和設(shè)計(jì),(4 4)計(jì)算)計(jì)算流程與流程與MatlabMatlab程序,(程序,(5 5)計(jì)算結(jié)果與分析。)計(jì)算結(jié)果與分析。高等工程數(shù)學(xué)高等工程數(shù)學(xué)中南大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院中南大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院非線性方程求根的數(shù)值解法非線性方程求根的數(shù)值解法第四章第四章 養(yǎng)老保險(xiǎn)問題養(yǎng)老保險(xiǎn)問題 非線性方程求根的數(shù)值解法非線性方程求根的數(shù)值解法養(yǎng)老保險(xiǎn)問題養(yǎng)老保險(xiǎn)問題4.1非線性方程求根的數(shù)值方法非線性方程求根的數(shù)值方法4.2養(yǎng)老保險(xiǎn)模型的求解養(yǎng)老保險(xiǎn)模型的求解4.34.1.1 4.1.1 問題的引入問題的引入 養(yǎng)老保險(xiǎn)是保險(xiǎn)

3、中的一種重要險(xiǎn)種,保養(yǎng)老保險(xiǎn)是保險(xiǎn)中的一種重要險(xiǎn)種,保險(xiǎn)公司將提供不同的保險(xiǎn)方案供以選擇,險(xiǎn)公司將提供不同的保險(xiǎn)方案供以選擇,分析保險(xiǎn)品種的實(shí)際投資價(jià)值。也就是說,分析保險(xiǎn)品種的實(shí)際投資價(jià)值。也就是說,如果已知所交保費(fèi)和保險(xiǎn)收入,則按年或如果已知所交保費(fèi)和保險(xiǎn)收入,則按年或按月計(jì)算實(shí)際的利率是多少?或者說,保按月計(jì)算實(shí)際的利率是多少?或者說,保險(xiǎn)公司需要用你的保費(fèi)實(shí)際至少獲得多少險(xiǎn)公司需要用你的保費(fèi)實(shí)際至少獲得多少利潤才能保證兌現(xiàn)你的保險(xiǎn)收益?利潤才能保證兌現(xiàn)你的保險(xiǎn)收益?4.1 4.1 養(yǎng)老保險(xiǎn)問題養(yǎng)老保險(xiǎn)問題4.1.2 4.1.2 模型分析模型分析 假設(shè)每月交費(fèi)假設(shè)每月交費(fèi)200200元至

4、元至6060歲開始領(lǐng)取歲開始領(lǐng)取養(yǎng)老金,養(yǎng)老金,投保人投保人若若2525歲起投保,屆時(shí)養(yǎng)歲起投保,屆時(shí)養(yǎng)老金每月老金每月22822282元;如元;如3535歲起保,屆時(shí)月歲起保,屆時(shí)月養(yǎng)老金養(yǎng)老金10561056元;試求出保險(xiǎn)公司為了兌元;試求出保險(xiǎn)公司為了兌現(xiàn)保險(xiǎn)責(zé)任,每月至少應(yīng)有多少投資收現(xiàn)保險(xiǎn)責(zé)任,每月至少應(yīng)有多少投資收益率?這也就是投保人的實(shí)際收益率。益率?這也就是投保人的實(shí)際收益率。4.1.3 4.1.3 模型假設(shè)模型假設(shè) 這是一個(gè)過程分析模型問題。過程的結(jié)果在條件一這是一個(gè)過程分析模型問題。過程的結(jié)果在條件一定時(shí)是確定的。定時(shí)是確定的。 整個(gè)過程可以按月進(jìn)行劃分,因?yàn)榻毁M(fèi)是按月進(jìn)行

5、整個(gè)過程可以按月進(jìn)行劃分,因?yàn)榻毁M(fèi)是按月進(jìn)行的。假設(shè)投保人到第的。假設(shè)投保人到第 k 月止所交保費(fèi)及收益的累計(jì)總月止所交保費(fèi)及收益的累計(jì)總額為額為 Fk ,每月收益率為,每月收益率為 r ,用,用 p, , q 分別表示分別表示6060歲之歲之前和之后每月交費(fèi)數(shù)和領(lǐng)取數(shù),前和之后每月交費(fèi)數(shù)和領(lǐng)取數(shù),N N表示停交保險(xiǎn)費(fèi)的月表示停交保險(xiǎn)費(fèi)的月份,份,M M表示停領(lǐng)養(yǎng)老金的月份。表示停領(lǐng)養(yǎng)老金的月份。4.1.4 4.1.4 模型建立模型建立v 在整個(gè)過程中,離散變量在整個(gè)過程中,離散變量 Fk 的變化規(guī)律滿足:的變化規(guī)律滿足:v 在這里在這里 Fk 實(shí)際上表示從保險(xiǎn)人開始交納保險(xiǎn)費(fèi)以后,實(shí)際上表示

6、從保險(xiǎn)人開始交納保險(xiǎn)費(fèi)以后,保險(xiǎn)人賬戶上的資金數(shù)值。保險(xiǎn)人賬戶上的資金數(shù)值。 11(1),0,1,.,1(1), N+1 .kkkkFFrp kNFFrq kNM,4.1.4 4.1.4 模型建立模型建立v 我們關(guān)心的是,在第我們關(guān)心的是,在第 M M 月時(shí),月時(shí), FM 能否為非負(fù)能否為非負(fù) 數(shù)?如果為正,則表明保險(xiǎn)公司獲得收益;如為負(fù),則數(shù)?如果為正,則表明保險(xiǎn)公司獲得收益;如為負(fù),則表明保險(xiǎn)公司出現(xiàn)虧損。當(dāng)為零時(shí),表明保險(xiǎn)公司最后表明保險(xiǎn)公司出現(xiàn)虧損。當(dāng)為零時(shí),表明保險(xiǎn)公司最后一無所一無所獲獲,所有的收益全歸保險(xiǎn)人,把它作為保險(xiǎn)人的,所有的收益全歸保險(xiǎn)人,把它作為保險(xiǎn)人的實(shí)際收益。實(shí)際收

7、益。 v 從這個(gè)分析來看,引入變量從這個(gè)分析來看,引入變量 FK ,很好地刻畫了整,很好地刻畫了整個(gè)過程中資金的變化關(guān)系;特別是引入收益率個(gè)過程中資金的變化關(guān)系;特別是引入收益率 r r,雖然,雖然它不是我們所求的保險(xiǎn)人的收益率,但從問題系統(tǒng)環(huán)境它不是我們所求的保險(xiǎn)人的收益率,但從問題系統(tǒng)環(huán)境中來看,必然要考慮引入另一對象中來看,必然要考慮引入另一對象保險(xiǎn)公司的經(jīng)營保險(xiǎn)公司的經(jīng)營效益,以此作為整個(gè)過程中各量變化的表現(xiàn)基礎(chǔ)。效益,以此作為整個(gè)過程中各量變化的表現(xiàn)基礎(chǔ)。4.1.5 4.1.5 模型求解模型求解 在在(4.1.4)(4.1.4)中兩式,中兩式,(可(可取初始值取初始值F0=0 ),我

8、們可以得到:,我們可以得到:MNkrrqrFFNkrrprFFNkNkNkkkk,.,1,1)1()1(,.,2 , 1 , 0,1)1()1 (0 再分別取,再分別取,k=Nk=N和和k=Mk=M,并利用,并利用F FM M=0=0可以可以導(dǎo)導(dǎo)出:出:0)1)(1 ()1 (pqrpqrNMM它是一個(gè)非線性方程。它是一個(gè)非線性方程。如何求解此方程?如何求解此方程? 代數(shù)方程求根問題是一個(gè)古老的數(shù)學(xué)問題。早在代數(shù)方程求根問題是一個(gè)古老的數(shù)學(xué)問題。早在1616世紀(jì)就找到了三次、四次方程的求根公式。世紀(jì)就找到了三次、四次方程的求根公式。 但直到但直到1919世紀(jì)才證明了世紀(jì)才證明了 次的一般代數(shù)方

9、程式是次的一般代數(shù)方程式是不能用代數(shù)公式求解的,因此需要研究用數(shù)值方法求得不能用代數(shù)公式求解的,因此需要研究用數(shù)值方法求得滿足一定精度的代數(shù)方程式的近似解。滿足一定精度的代數(shù)方程式的近似解。 在工程和科學(xué)技術(shù)中許多問題常歸結(jié)為求解非線性在工程和科學(xué)技術(shù)中許多問題常歸結(jié)為求解非線性方程式問題。正因?yàn)榉蔷€性方程求根問題是如此重要和方程式問題。正因?yàn)榉蔷€性方程求根問題是如此重要和基礎(chǔ),因此它的求根問題很早就引起了人們的興趣,并基礎(chǔ),因此它的求根問題很早就引起了人們的興趣,并得到了許多成熟的求解方法。得到了許多成熟的求解方法。 下面下面 我們首先了解一下非線性方程的基本概念。我們首先了解一下非線性方程

10、的基本概念。5n 定義定義4.2.14.2.1 設(shè)有一個(gè)非線性方程設(shè)有一個(gè)非線性方程 ,其中,其中 為實(shí)變?yōu)閷?shí)變量量 的非線性函數(shù)。的非線性函數(shù)。 (1 1)如果有)如果有 使使 ,則稱則稱 為為方程的根方程的根,或?yàn)?,或?yàn)?的零點(diǎn)。的零點(diǎn)。 (2 2)當(dāng))當(dāng) 為多項(xiàng)式,即為多項(xiàng)式,即 則稱則稱 為為 次次代數(shù)方程代數(shù)方程, 包含指數(shù)函數(shù)或者三角函包含指數(shù)函數(shù)或者三角函數(shù)等特殊函數(shù)時(shí),則稱數(shù)等特殊函數(shù)時(shí),則稱 為為特殊方程特殊方程。 (3 3)如果)如果 ,其中,其中 。 為正整數(shù),則為正整數(shù),則稱稱 為為 的的 重根重根。當(dāng)。當(dāng) 時(shí)稱時(shí)稱 為為 的的單根單根。4.2.1 根的搜索相關(guān)定義根的

11、搜索相關(guān)定義 0f x ( )f xxx()0f xx f xn f x110( )00nnnnnf xa xaxaxaa( ) 0f x ( )f x( )0f x ( ) ()( )mf xx xg x( )0g xmx( ) 0f x m1mx( ) 0f x 4.2 4.2 非線性方程求根的數(shù)值方法非線性方程求根的數(shù)值方法 定理定理4.2.14.2.1 設(shè)設(shè) 為具有復(fù)系數(shù)的為具有復(fù)系數(shù)的 次代數(shù)方程,則次代數(shù)方程,則 在復(fù)數(shù)域上恰有在復(fù)數(shù)域上恰有 個(gè)根個(gè)根( 重根計(jì)算重根計(jì)算 個(gè))。如果個(gè))。如果 為實(shí)系數(shù)方程,則復(fù)數(shù)根成對出現(xiàn),即當(dāng)為實(shí)系數(shù)方程,則復(fù)數(shù)根成對出現(xiàn),即當(dāng): : 為為 的

12、復(fù)根,則的復(fù)根,則 亦是亦是 的根。的根。 定理定理4.2.24.2.2 設(shè)設(shè) 在在 連續(xù)連續(xù), ,且且 ,則存,則存在在 ,使得,使得 ,即,即 在在 內(nèi)存在實(shí)零點(diǎn)。內(nèi)存在實(shí)零點(diǎn)。( )f x,a b( )( ) 0f a f b,xa b( )0f x f x, )a b(( ) 0f x n( ) 0f x nrr( )0f x 0i ( ) 0f x ( ) 0f x i4.2.2 4.2.2 逐步搜索法逐步搜索法v 對于方程對于方程 , ,為明確起見,為明確起見,設(shè)設(shè) , , ,從區(qū)間左端點(diǎn),從區(qū)間左端點(diǎn) 出發(fā)按某個(gè)預(yù)出發(fā)按某個(gè)預(yù)定步長定步長 (如取(如取 , 為正整數(shù)),一步一步為

13、正整數(shù)),一步一步地向右跨,每跨一步進(jìn)行一次根的收索。即檢查節(jié)點(diǎn)地向右跨,每跨一步進(jìn)行一次根的收索。即檢查節(jié)點(diǎn) 上的函數(shù)值上的函數(shù)值 的符號,若的符號,若 ,則,則 即即為方程解。若為方程解。若 ,則方程根在區(qū)間,則方程根在區(qū)間 中,其中,其寬度為寬度為 。 bahN 0f x ,xa b 0f a ( )0f b 0 xahNkxa kh kf x 0kf x kx 0kf x 1,kkxxh4.2.2 4.2.2 逐步搜索法逐步搜索法 例例4.2.14.2.1 考察方程考察方程 由于由于 則則 在在 內(nèi)至少有一個(gè)根,內(nèi)至少有一個(gè)根,設(shè)從設(shè)從 出發(fā),以出發(fā),以 為步長向右進(jìn)行根的搜索。列為步

14、長向右進(jìn)行根的搜索。列表記錄各節(jié)點(diǎn)函數(shù)值的符號??梢娫诒碛涗浉鞴?jié)點(diǎn)函數(shù)值的符號。可見在 內(nèi)必有一根。內(nèi)必有一根。 表表4.2.14.2.1 的符號的符號x00.51.01.5 的符號-+ 31 0f xx x f x f x 01 0,25 0ff 0,20 x 0.5h1.0,1.5( )f x 易見此方法應(yīng)用關(guān)鍵在步長易見此方法應(yīng)用關(guān)鍵在步長 的選擇上。的選擇上。很明顯,只要步長很明顯,只要步長 取得足夠小,利用此法就取得足夠小,利用此法就可以得到任意精度的根,但可以得到任意精度的根,但 縮小,搜索步數(shù)縮小,搜索步數(shù)增多,從而使計(jì)算量增大,用此方法對高精度增多,從而使計(jì)算量增大,用此方法對

15、高精度要求不合適。要求不合適。hhh4.2.3 4.2.3 二分法二分法 對非線性方程:對非線性方程: 其中其中 在在 連續(xù)且連續(xù)且 無妨設(shè)無妨設(shè) 在在 內(nèi)僅有一個(gè)零點(diǎn)。內(nèi)僅有一個(gè)零點(diǎn)。 求方程(求方程( )的實(shí)根)的實(shí)根 的二分法過程,的二分法過程,就是將就是將 逐步分半,檢查函數(shù)值符號的變化,逐步分半,檢查函數(shù)值符號的變化,以便確定包含根的充分小區(qū)間。以便確定包含根的充分小區(qū)間。 0 f x 4.2.1 f x,a b 0f a f b f x,a b4.2.1x,a b 二分法的步驟如下二分法的步驟如下:記記 , 第第1 1步:分半計(jì)算步:分半計(jì)算 ,將,將 分半。計(jì)算中點(diǎn)分半。計(jì)算中點(diǎn)

16、 及及 。若。若 ,則根必在,則根必在 內(nèi),否則必在內(nèi),否則必在 內(nèi),內(nèi),(若(若 ,則,則 ),于是得到長度一半的區(qū)間),于是得到長度一半的區(qū)間 含根含根, ,即即 , ,且且 。 第第 步(步( k 次次分半計(jì)算)重復(fù)上述過程。分半計(jì)算)重復(fù)上述過程。1aa1bb1112abx1122,a xab:1122, ,xbab:22111()2baba1k11a,b 1f x11( )( ) 0f af x1( ) 0f x 1xx22,a b22() ( )0f af bk 設(shè)已完成第設(shè)已完成第1 1步步第第 步,分半計(jì)算得到含根步,分半計(jì)算得到含根區(qū)間區(qū)間 ,且滿,且滿足足 , ,即即 ,

17、即即 , 則第則第k k步的分半計(jì)算:步的分半計(jì)算: ,且有:,且有:1122,kkabababL , kkxa b11()2kkkbaba2kkkabx1 22kkkkbaxxba 4 .2 .21k( ) ( ) 0kkf a f b 確定新的含根區(qū)間確定新的含根區(qū)間 ,即如果,即如果 ,則根必在則根必在 內(nèi),否則必內(nèi),否則必在在 內(nèi),且有:內(nèi),且有: 。總之,由上總之,由上述二分法得到序列述二分法得到序列 , 由由(4.2.2)(4.2.2)有:有: 。 可用二分法求方程可用二分法求方程 的實(shí)根的實(shí)根 的近似值的近似值到任意指定的精度,這是因?yàn)椋涸O(shè)到任意指定的精度,這是因?yàn)椋涸O(shè) 為給定為

18、給定精度要求,則由精度要求,則由 ,可得分,可得分半計(jì)算次數(shù)半計(jì)算次數(shù) k 應(yīng)滿足:應(yīng)滿足:11,kkab( ) ( ) 0kkf a f x 11,kkkkabax:11,kkkkabxb:111()2kkkbaba kxlimkkxx0( ) 0f x x2kkbaxx 二分法的優(yōu)點(diǎn)是方法簡單,且只要求二分法的優(yōu)點(diǎn)是方法簡單,且只要求 連續(xù)即可,可用二分法求出連續(xù)即可,可用二分法求出 在在 內(nèi)全部實(shí)根,但二分法不能求復(fù)內(nèi)全部實(shí)根,但二分法不能求復(fù)根及偶數(shù)重根,且收斂較慢,函數(shù)值計(jì)算次數(shù)根及偶數(shù)重根,且收斂較慢,函數(shù)值計(jì)算次數(shù)較多。較多。 lnlnln 2bak 4.2.3( )f x()0

19、fx,ab 例例4.2.24.2.2 用二分法求用二分法求 在在1,21,2內(nèi)一內(nèi)一個(gè)實(shí)根,且要求精確到小數(shù)點(diǎn)后第三位。個(gè)實(shí)根,且要求精確到小數(shù)點(diǎn)后第三位。(即(即 ) 解解 由由 代入公式代入公式(4.2.3) (4.2.3) ,可確定所需分半次數(shù)為可確定所需分半次數(shù)為 ,計(jì)算結(jié)果部分如下表:,計(jì)算結(jié)果部分如下表:(顯然(顯然 )。)。6( )1f xxx*3121 0kxx30 .51 01,2)ab(11k(1)10,(2)0ff K81.132813 1.140625 1.136719 0.020619 91.132813 1.136719 1.134766 0.4268415 101

20、.132813 1.134766 1.133789 111.133789 1.134766 1.134277 表表4.2.2 4.2.2 部分計(jì)算結(jié)果部分計(jì)算結(jié)果kakbkx()kf x00959799. 00045915. 04.2.4 4.2.4 迭代法迭代法 迭代法迭代法是一種逐次逼近法。它是求解代數(shù)方程,超是一種逐次逼近法。它是求解代數(shù)方程,超越方程及方程組的一種基本方法,但越方程及方程組的一種基本方法,但存在收斂性及收斂存在收斂性及收斂快慢的問題快慢的問題。 用迭代法求解用迭代法求解 的近似根,首先需將此方的近似根,首先需將此方程化為等價(jià)的方程:程化為等價(jià)的方程: 然而將然而將 化為

21、等價(jià)方程化為等價(jià)方程 的方法是的方法是很多的。很多的。 ( ) 0f x (4.2.4)( )0f x ( ) x g x (4.2.4) 例例4.2.34.2.3 對方程對方程 可用不同的方法將其化為等價(jià)方程:可用不同的方法將其化為等價(jià)方程: (1 1) (2 2)()sin0.50fxxx1sin0.5( )xxgx12sin0.5( )xxgx定義定義4.2.2 (迭代法)(迭代法)設(shè)方程為設(shè)方程為 , ,取方程根的一個(gè)初始近取方程根的一個(gè)初始近似似 ,且按下述逐次代入法,構(gòu)造一個(gè)近似解序列:,且按下述逐次代入法,構(gòu)造一個(gè)近似解序列: 這種方法稱為這種方法稱為迭代法迭代法(或稱為單點(diǎn)迭代

22、法),(或稱為單點(diǎn)迭代法), 稱為稱為迭代函數(shù)迭代函數(shù)。若由迭代法產(chǎn)生序列若由迭代法產(chǎn)生序列 有極限存在,即有極限存在,即 ,稱,稱 為為收斂或迭代過程收斂或迭代過程 收斂,否則稱迭代法不收斂。若收斂,否則稱迭代法不收斂。若 連續(xù),連續(xù),且且 ,則,則 , 即即 為方程為方程 的解(稱的解(稱 為函數(shù)為函數(shù) 的不動點(diǎn)),顯然的不動點(diǎn)),顯然在由方程在由方程 轉(zhuǎn)化為等價(jià)方程轉(zhuǎn)化為等價(jià)方程 時(shí),選擇不同的迭代時(shí),選擇不同的迭代函數(shù)函數(shù) ,就會產(chǎn)生不同的序列,就會產(chǎn)生不同的序列 (即使初值(即使初值 選擇一樣)且選擇一樣)且這些序列的收斂情況也不會相同。這些序列的收斂情況也不會相同。( )xg x0

23、 x 10211, kkxg xxg xxg x (4.2.5)( )g x kx*limkkxxkx (4.2.5)*limkkxx1limlim ( )lim( )kkkkkkxxg xgxg x g x (4.2.4)( )g xxx( )0f x ( )xg x( )g x kx0 x例例4.2.44.2.4 對例對例4.2.14.2.1中方程考查用迭代法求根中方程考查用迭代法求根 由計(jì)算可以看出,我們選取的兩個(gè)函數(shù)由計(jì)算可以看出,我們選取的兩個(gè)函數(shù) ,分別構(gòu)造序列分別構(gòu)造序列 收斂情形不一樣(初值都取為收斂情形不一樣(初值都取為1 1),),在在 中中 收斂且收斂且 ,在,在 中計(jì)算

24、出中計(jì)算出 無定義。無定義。 111sin0.5,0,1,2,sin0.5 ,0,1,2,kkkkaxxkbxxk 12,g x g x kx( )a kx1.497300 x( )b114sin0.5sin1.987761x0 1.0 1.0 11.341471 0.523599 21.473820 0.023601 31.049530 -0.496555 41.497152 -1.487761 51.497289 61.497300 71.497300 k( )kax( ) kbx ( )kaf x73.6 10表表4.2.3 4.2.3 部分計(jì)算結(jié)果部分計(jì)算結(jié)果 因此對用迭代法求方程因此

25、對用迭代法求方程 的近似根,的近似根,需要研究下述問題:需要研究下述問題: (1 1)如何選取迭代函數(shù))如何選取迭代函數(shù) 使迭代過程使迭代過程 收斂。收斂。 (2 2)若)若 收斂較慢時(shí),怎樣加速收斂較慢時(shí),怎樣加速 收斂。收斂。 ( )0f x ( )g x 1kkxg x kx kx 迭代法的幾何意義:迭代法的幾何意義: 從幾何意義看,求方程從幾何意義看,求方程 根的問題,是求曲線根的問題,是求曲線 與直線與直線 交點(diǎn)的橫坐標(biāo)交點(diǎn)的橫坐標(biāo) ,當(dāng)?shù)瘮?shù),當(dāng)?shù)瘮?shù) 的導(dǎo)數(shù)函數(shù)的導(dǎo)數(shù)函數(shù) 在根在根 處滿足下述幾種條件時(shí),從幾何上來看迭代過程處滿足下述幾種條件時(shí),從幾何上來看迭代過程 的收斂情

26、況如圖的收斂情況如圖4.2.14.2.1。 從曲線從曲線 上一點(diǎn)上一點(diǎn) 出發(fā),沿著平行于出發(fā),沿著平行于x x軸方向前進(jìn)交軸方向前進(jìn)交 于一點(diǎn)于一點(diǎn) 再從點(diǎn)再從點(diǎn) 沿平行于沿平行于y y軸方向前軸方向前進(jìn)交進(jìn)交 于于 點(diǎn),顯然點(diǎn),顯然 的橫坐標(biāo)就是的橫坐標(biāo)就是 ,繼續(xù),繼續(xù)這過程就得到序列這過程就得到序列 ,且從幾何上觀察知道在(,且從幾何上觀察知道在(1 1),(),(2 2)情)情況下況下 收斂于收斂于 ,在(,在(3 3),(),(4 4)情況)情況 不收斂于不收斂于 。 ( )xg x( )yg xyxx( )g x xgx1kkxg x( )yg x000,px g x0Q0Q( )

27、y gx1p1p10 xg xkx kx*x kx*xyx圖圖4.2.1 4.2.1 迭代法的幾何意義圖迭代法的幾何意義圖 由迭代法的幾何由迭代法的幾何意義意義知,為了保證迭代過程收斂,應(yīng)該要求迭代函知,為了保證迭代過程收斂,應(yīng)該要求迭代函數(shù)的導(dǎo)數(shù)滿足條件數(shù)的導(dǎo)數(shù)滿足條件 。當(dāng)。當(dāng) 時(shí),原方程在時(shí),原方程在 中可能中可能有幾個(gè)根或迭代法不收斂,為此有關(guān)于迭代收斂性的定理有幾個(gè)根或迭代法不收斂,為此有關(guān)于迭代收斂性的定理4.2.34.2.3。 定理定理4.2.34.2.3 設(shè)有方程設(shè)有方程 , (1) (1) 設(shè)設(shè) 于于 一階導(dǎo)數(shù)存在,一階導(dǎo)數(shù)存在, (2) (2) 當(dāng)當(dāng) 時(shí),有時(shí),有 , (3

28、) (3) 滿足條件:滿足條件: 則有:則有: 在在 上有唯一解上有唯一解 , 對任意選取初始值對任意選取初始值 ,迭代過程,迭代過程 收斂即收斂即 , 誤差估計(jì)誤差估計(jì)()1gx , xab , a bg( )xxg( ) xa,ba,bxg( ) a,bx ( )1, , ,g xLxa b g( ) x1g( )xx , a b*x3*11 1kkkxxxxL4*10 ,(1,2,.)1kkLxxxxkL20 , xab*klim xx1(),k0,1,.kkxg x證明證明 只證只證 , , 由定理?xiàng)l件由定理?xiàng)l件 ,當(dāng)取,當(dāng)取 時(shí),則有時(shí),則有 記誤差記誤差 ,由中值定理有:,由中值定

29、理有: ,其中其中 在在 與與 之間,即之間,即 ,又由條件,又由條件有:有: ,由此遞推可,由此遞推可得:得: ,由,由 故故 。 由迭代公式由迭代公式 有:有: ,其中其中c c在在 與與 之間,于是:之間,于是: 即即 。 234*111*() -L x(1-L) kkkkkkkkkxxxxxxxxxxxxxxx*11111kkkkkLxxxxxxLL*klimxx20 , xa b , ,(1,2,.)kxa bk*kkexx*1()()( )()kkkxxg xg xg c xxc*xkx , ca b*2*120.kkkkxxL xxL xxL xx01L*1( )kkkxxg c

30、 xxL xx321()kkxg x1111( )()( )()kkkkkkkkxxg xg xg c xxLxx1kxkx3 由上面由上面 反復(fù)利用代入上式中有反復(fù)利用代入上式中有: 由定理由定理 結(jié)果可知,當(dāng)計(jì)算得到的相鄰兩次迭代滿足條件結(jié)果可知,當(dāng)計(jì)算得到的相鄰兩次迭代滿足條件 時(shí),則誤差時(shí),則誤差 。 因此在計(jì)算機(jī)上可利用因此在計(jì)算機(jī)上可利用 來控制算法終止,但要來控制算法終止,但要注意注意 時(shí),即使時(shí),即使 很小,誤差很小,誤差 仍然可能很大。仍然可能很大。 另外,當(dāng)已知另外,當(dāng)已知 及及 及給定精度要求及給定精度要求 時(shí),利用定時(shí),利用定理理 結(jié)果可確定使誤差達(dá)到給定精度要求時(shí)所需

31、要迭代次數(shù)結(jié)果可確定使誤差達(dá)到給定精度要求時(shí)所需要迭代次數(shù)k k,事,事實(shí)上,由實(shí)上,由 解得:解得: 411kkkkxxLxx*1121210111L .11kkkkkkkkLxxxxxxLLLxxxxLL31kkxx*11kxxL*kxx1kkxx1L 1kkxx01,x x(1)L L4*101kkLxxxxL10(lnln)/ln 1xxkLL 定理?xiàng)l件定理?xiàng)l件 ,在一般情況下,可能對大,在一般情況下,可能對大范圍的含根區(qū)間不滿足,而在根的鄰近是成立的,為此有范圍的含根區(qū)間不滿足,而在根的鄰近是成立的,為此有如下迭代過程的局部收斂性結(jié)果。如下迭代過程的局部收斂性結(jié)果。 定理定理4.2.

32、44.2.4 (迭代法的局部收斂性)設(shè)給定方程(迭代法的局部收斂性)設(shè)給定方程 (1 1)設(shè))設(shè) 為方程的解,為方程的解, (2 2)設(shè))設(shè) 在在 的鄰域內(nèi)連續(xù)可微,且有的鄰域內(nèi)連續(xù)可微,且有 ,則對任意初值則對任意初值 (在(在 的鄰域內(nèi)),迭代過程的鄰域內(nèi)),迭代過程 , 收斂于收斂于 。( )1, , g xLxa bg( )xx*xg( ) x*x*()1gx0 x*x1( )kxg x0,1,.k *x 例例4.2.54.2.5 用用迭代法解方程迭代法解方程( )ln(2)0f xxx 解解 (1)(1)顯然有顯然有 即知方程于即知方程于0,20,2及及-1.9,-1-1.9,-1內(nèi)

33、有根記為內(nèi)有根記為 。 (2) (2)考察取初值考察取初值 迭代過程迭代過程 的收斂性,的收斂性,其中迭代函數(shù)為其中迭代函數(shù)為 ,顯然,顯然 , ,及,及 為增函為增函數(shù),則當(dāng)數(shù),則當(dāng) 時(shí),時(shí), ,又由,又由 則有則有 。 于是由定理于是由定理4.2.44.2.4可知,當(dāng)初值可知,當(dāng)初值 時(shí),迭代過時(shí),迭代過程程 收斂,如果要求收斂,如果要求 的近似根準(zhǔn)確到小數(shù)的近似根準(zhǔn)確到小數(shù)點(diǎn)后第點(diǎn)后第6 6位(即要求位(即要求 )由計(jì)算結(jié)果可)由計(jì)算結(jié)果可知知 。且。且 則則 , 。 (0) (2)0, ( 1.9) ( 1)0ffff 12,xx00,2x 1ln(2)kkxx1( ) ln(2)g

34、xx1(0) ln(2) 0.693 0g1(2) ln(4) 1.368 2g1( )g x02x10( )2g x1()2gxx111( )(0)1(0,2)22g xgxx 00,2x 611102kx x 1ln(2)kkxx1x7151410 xx12L *1141.1461931xx714()0.8 10f x 表表4.2.4 4.2.4 部分計(jì)算結(jié)果表部分計(jì)算結(jié)果表k1ln(2)kkxx0 00.00.01 10.693147180.693147182 20.990710460.9907104614141.14619311.146193115151.14619321.146193

35、2 (3) (3)為了求為了求-1.9,-1-1.9,-1內(nèi)方程的根。由迭代方程內(nèi)方程的根。由迭代方程 ,顯,顯然然 ,所以迭代過程,所以迭代過程 (初值(初值 )不能保證收斂于)不能保證收斂于 。 (4) (4)若將方程轉(zhuǎn)化為等價(jià)方程若將方程轉(zhuǎn)化為等價(jià)方程 或或 則則 ,且,且 ( 時(shí))時(shí)) 所以當(dāng)選取所以當(dāng)選取 時(shí)迭代過程時(shí)迭代過程 收斂。如取收斂。如取 ,則迭代則迭代1212次有次有 ,且,且 。 由此例可見,對于方程由此例可見,對于方程 ,迭代函數(shù),迭代函數(shù) 取不同形式,取不同形式,相應(yīng)的迭代法產(chǎn)生的相應(yīng)的迭代法產(chǎn)生的 收斂情況也不一樣,因此,我們應(yīng)該選擇收斂情況也不一樣,因此,我們應(yīng)

36、該選擇迭代函數(shù)時(shí)構(gòu)造的迭代過程迭代函數(shù)時(shí)構(gòu)造的迭代過程 收斂,且收斂較快。收斂,且收斂較快。1ln(2)kkxx*121( )()12gxgxx1ln(2)kkxx*002 1.9, 1,xxx *2xe2xx 22( )xxegx2g ( )xxe22g ( )g ( 1)0.386 1x 1.9, 1x2( ) 1.9, 1g x 0 1.9, 1x k 12kxxe01x *2121.841405660 xx 812()0.2 10f x ( ) 0f x ( )g x kx1()kkxg x關(guān)于迭代公式的加工:關(guān)于迭代公式的加工: 對于收斂的迭代過程,只要迭代足夠多次,總可以使結(jié)果達(dá)到

37、對于收斂的迭代過程,只要迭代足夠多次,總可以使結(jié)果達(dá)到任意的精度。但有時(shí)迭代收斂緩慢,從而使計(jì)算量變得很大,因此任意的精度。但有時(shí)迭代收斂緩慢,從而使計(jì)算量變得很大,因此迭代過程的加速是一個(gè)很重要的課題。迭代過程的加速是一個(gè)很重要的課題。 設(shè)設(shè) 為根為根 的某個(gè)預(yù)測值,用迭代公式校正一次得:的某個(gè)預(yù)測值,用迭代公式校正一次得: 由中值定理:由中值定理: , 介于介于 之間,若之間,若 改改變不大。近似地取某常數(shù)變不大。近似地取某常數(shù) ,則由,則由 可以期望按上式右端求得的可以期望按上式右端求得的 是比更好的近似值。是比更好的近似值。 若將每得到一次改進(jìn)值算作一步,并用若將每得到一次改進(jìn)值算作一

38、步,并用 和和 分別表示第分別表示第 步的校正值和改進(jìn)值,則步的校正值和改進(jìn)值,則加速迭代計(jì)算方案加速迭代計(jì)算方案如下:如下:預(yù)測預(yù)測:校正校正: 2101101()111LLxxxxxxLLLxkxkk1 ()kxg xk 111k ( )1kkLxxxxL10()xg x*10( )()x xgxx *0 xx,( )g x*10101()11Lx xLxxxxxLL 0 x*xL0111150 0.5()0.60.6 ()1.60.5310 kxxkkkkkxexexxxxx由于在附近時(shí),則預(yù)測校正迭代公式為:預(yù)測:校正:只須由迭代 次即可得精確到的結(jié)果。18若直接用迭代法求,達(dá)到這一結(jié)

39、果須迭代次。-* e0.56714.xxx例: 求解方程。其解為: 設(shè)設(shè) 的某近似值的某近似值 ,將校正值,將校正值 再校正一次再校正一次得:得: ,由,由 與與 得:得: 由此得:由此得: 。這樣將上式右端作為。這樣將上式右端作為改進(jìn)公式就不再含有導(dǎo)數(shù)信息了。但需要用到兩次迭代的結(jié)果進(jìn)行改進(jìn)公式就不再含有導(dǎo)數(shù)信息了。但需要用到兩次迭代的結(jié)果進(jìn)行加工。如果仍將得到一次改進(jìn)值作為一步,則計(jì)算過程如下:加工。如果仍將得到一次改進(jìn)值作為一步,則計(jì)算過程如下: 上述處理過程稱為(上述處理過程稱為(埃特金埃特金AitkenAitken)方法。)方法。*x0 x10()xg x21( )xg x*21()

40、xxLx x *10()()x xLx x*01*21xxxxxxxx2*02121201201222x xxxxxxxxxxxx2()k11k1k1k1k12k1k1k1 () () 2kkkxg xxg xxxxxxxx校 正 :再 校 正 :()改 進(jìn) :由于使用參數(shù)由于使用參數(shù) ,這在實(shí)際應(yīng)用中不方便,下面進(jìn)行改進(jìn)計(jì)算。,這在實(shí)際應(yīng)用中不方便,下面進(jìn)行改進(jìn)計(jì)算。3: Aitken 10 xx 例用方法求解方程。30 1( )(1.5),:xxg xxAitken 解:若將方程改為等價(jià)方程:則直接用迭代法 取是不收斂的 現(xiàn)在以這種迭代法公式為基礎(chǔ)形成算法如下.32472. 1105*5x

41、的近似解次可得精度為則迭代k 1k 1k 13k 13k 1102k 11k 1 1 1 (1.5) =2kkkkxxxxxxxxxxxx()4.2.5 Newton4.2.5 Newton公式公式 對于方程對于方程 ,應(yīng)用迭代法時(shí)先要改寫成,應(yīng)用迭代法時(shí)先要改寫成 ,即需要針對即需要針對 構(gòu)造不同的合適的迭代函數(shù)構(gòu)造不同的合適的迭代函數(shù) 。 顯然可以取迭代函數(shù)為顯然可以取迭代函數(shù)為 ,相應(yīng)迭代公式,相應(yīng)迭代公式為為 。 一般地,這種迭代公式不一定收斂,或者速度很慢。對此公一般地,這種迭代公式不一定收斂,或者速度很慢。對此公式應(yīng)用前面的加速技術(shù)具體格式為:式應(yīng)用前面的加速技術(shù)具體格式為: (

42、)0f x ( )xg x( )f x( )g x( )( )g xxf x 1()kkkxxf x1111()()1kkkkkkkxxf xLxxxxL 記記 ,則上二式可合并寫為:,則上二式可合并寫為: 。此公式稱為簡單的此公式稱為簡單的NewtonNewton公式,迭代函公式,迭代函數(shù)數(shù)為:為: 。 又由于又由于 為為 的近似值,而的近似值,而 ,因此,因此 實(shí)際上是實(shí)際上是 的近似值,故用的近似值,故用 代替上式中的代替上式中的 即得到下面的迭代函數(shù):即得到下面的迭代函數(shù): 。 相應(yīng)的迭代公式為:相應(yīng)的迭代公式為: , 即為即為NewtonNewton公式。公式。1ML1()kkkf

43、xxxM( )( )fxg xxML( )g x( )( )g xxfx1ML( )f x( )fxM( )( )( )f xg xxfx1()()kkkkf xxxfx4.2.6 Newton4.2.6 Newton法的幾何意義法的幾何意義牛頓法牛頓法的基本思想的基本思想就是將非線性方程就是將非線性方程 逐步線性化求解,設(shè)逐步線性化求解,設(shè) 有近似的根有近似的根 ,將,將 在在 處處 展開得:展開得: 從而從而 近似地表為:近似地表為: 。方程。方程 的根的根 即為曲線即為曲線 與與 軸交軸交點(diǎn)的橫坐標(biāo)。設(shè)點(diǎn)的橫坐標(biāo)。設(shè) 為為 近似近似值值,過曲線,過曲線 上橫坐標(biāo)上橫坐標(biāo) 為為 的點(diǎn)的點(diǎn)

44、作曲線的切線,該切線作曲線的切線,該切線 與與 軸軸交交點(diǎn)的橫坐標(biāo)即為點(diǎn)的橫坐標(biāo)即為 的新近似的新近似 值值 ,它與,它與 軸交點(diǎn)的橫坐標(biāo)為:軸交點(diǎn)的橫坐標(biāo)為: ,因此,因此 NewtonNewton法亦稱切線法法亦稱切線法。( )f x( )0f x ( )0f x kxTaylor( )( )( )()kkkf xf xf x x xkxkx( )0f x ()()()0kkkf xfxxx( )0f x ( )yf xx*x( )yf xkxkpx*x1kxx1()/()kkkkxxf xf x4.2.7 Newton4.2.7 Newton法的局部收斂性法的局部收斂性 定義定義4.2.

45、34.2.3 設(shè)迭代過程設(shè)迭代過程 收斂于方程收斂于方程 的的根根 ,如果迭代誤差,如果迭代誤差 ,當(dāng),當(dāng) 時(shí)有:時(shí)有: 則稱該迭代過程則稱該迭代過程為為 階收斂階收斂的。的。 定理定理4.2.54.2.5 對迭代過程對迭代過程 如果如果 在在 附近連續(xù),附近連續(xù),且:且: 且且 ,則該迭代過程在,則該迭代過程在 附近是附近是 階收斂的。階收斂的。1()kkxg x( )x gx*x*kkexxk 1 (0,)kpkecce為常數(shù)pk 1( ),kxg x( )( )pgx*x*(p-1)*g( ) g( ) . g( ) 0 xxx (p)*g ( ) 0 x *xp 證明證明 由于由于 ,

46、則有前面關(guān)于迭代法的局部收斂性定理知:此迭代,則有前面關(guān)于迭代法的局部收斂性定理知:此迭代過程過程 具有局部收斂性。即具有局部收斂性。即 。將。將 在在 處處 展開,展開,并注意到并注意到 有:有: 而,而, 從而上式化為:從而上式化為: *g( )0 1x 1( )kkxg x0kx ( )kg x*xTaylor*(p-1)*g( ). g( )0 xx( )*()()()() , , !ppkkkkkgg xg xxxx xp*k 1( ), ( )kxg xxg x()*k1()()!ppkkgxxxxp 即:即: 故知迭代過程具有故知迭代過程具有 階收斂性。階收斂性。 *()()*1

47、k1p*k()e()e()!ppkkpkxxggxxxppp 定理定理4.2.54.2.5表明迭代過程的收斂速度依賴于迭代函數(shù)表明迭代過程的收斂速度依賴于迭代函數(shù) 的選取,的選取,如果如果 時(shí)時(shí) 。則迭代過程只可能是線性收斂的。則迭代過程只可能是線性收斂的。( )g x , xa b( ) 0g x 對于對于NewtonNewton法,由迭代函數(shù)為:法,由迭代函數(shù)為: 則則 , 若若 為為 的一個(gè)單根。即的一個(gè)單根。即 ,則由上式知,則由上式知 。 由上面定理可知由上面定理可知NewtonNewton法在根法在根 的鄰域內(nèi)是平方收斂的(二階收斂的鄰域內(nèi)是平方收斂的(二階收斂 的)。的)。()

48、()-()fxgxxfx222( )( )( )( )( )g ( )1( )( )fxf x fxf x fxxfxfx*x( )f x *()00f xfx,*( ) 0g x *x 特別地考察特別地考察NewtonNewton公式:設(shè)公式:設(shè) 二次連續(xù)可微,二次連續(xù)可微,則則 , 在在 之間,特別之間,特別地取地取 ,注意,注意 ,則,則 設(shè)設(shè) 。兩邊同除以。兩邊同除以 ,得:,得: (注:(注: ),利用),利用NewtonNewton公式,即有:公式,即有: 當(dāng)當(dāng) ,則,則 , 或或( )f x2( )( )()()()()2!kkkkff xf xfxxxxx,kx x*xx*(

49、) 0f x *2( )0( )( )( )()()2!kkkkff xf xf xxxxx( )0kf x( )kf x*2()()()()2()kkkkkfxfxxxxfxfx1()()kkkkf xxxfx*1*2()()2()kkkkxxfmxxfx k *()()2()2()kffxfxfx *2211( )()2 ( )kkkk kkfxxexxmef x 可見可見 (誤差)與(誤差)與 的誤差的誤差 的平方成比例。當(dāng)初始誤差的平方成比例。當(dāng)初始誤差 充分小時(shí),以后迭代的誤差將減少得非??臁7粗浞中r(shí),以后迭代的誤差將減少得非常快。反之 ,則放大。則放大。NewtonNewton

50、法每計(jì)算一步,需要計(jì)算一次函數(shù)值法每計(jì)算一步,需要計(jì)算一次函數(shù)值 和一次導(dǎo)數(shù)和一次導(dǎo)數(shù)值值 。 例例4.2.64.2.6 用用NewtonNewton法求解法求解 。 解解 顯然顯然 。則在。則在0,20,2內(nèi)方程有一個(gè)根,求導(dǎo)內(nèi)方程有一個(gè)根,求導(dǎo) 則則NewtonNewton公式為:公式為: 取取 ,迭代,迭代6 6次得近似根為次得近似根為 , , 。這表明,當(dāng)初值這表明,當(dāng)初值 取值靠近取值靠近 時(shí),時(shí),NewtonNewton法收斂且收斂較快,否則法收斂且收斂較快,否則NewtonNewton法可能不收斂。法可能不收斂。(0) (2)0ff4( )(6) / 4xfxex414()(2)

51、1()(6) / 4kkxkkkkkxkkf xexxxxfxex1kekxke*00exx01e ()kf x( )kf x4( )(2)10 xf xex01.0 x *0.783596x *6( )3.8 10f x0 x*x 下面考慮下面考慮NewtonNewton法的誤差估計(jì),由中值定理有:法的誤差估計(jì),由中值定理有: ,當(dāng)當(dāng) 充分接近充分接近 時(shí),有時(shí),有 因此,用因此,用NewtonNewton法求方程單根法求方程單根 的近似根的近似根 的誤差的誤差 可用可用 來估計(jì)。來估計(jì)。 *( )( )( )( )()kkkkf xf xf xfxxkx*1( )( )( )( )Newt

52、onkkkkkkkf xf xxxxxff x*x*xkx*kkexx 1kkxx4.2.8 Newton4.2.8 Newton法應(yīng)用舉例法應(yīng)用舉例1. 1. 對給定的正數(shù)對給定的正數(shù) ,應(yīng)用,應(yīng)用NewtonNewton法解二次方程法解二次方程 可導(dǎo)出求可導(dǎo)出求開方值開方值 的計(jì)算格式:的計(jì)算格式: 可證明公式可證明公式 對任意函數(shù)初值對任意函數(shù)初值 都是收斂的。這是因?yàn)椋憾际鞘諗康?。這是因?yàn)椋?C20 xCC11() 2kkkCxxx (4.2.7)(4.2.7)00 x 21211()21()2kkkkkkxCxCxxCxCx 兩式相除得:兩式相除得: 利用此式遞推可得:利用此式遞推可

53、得: ( 由由 可知:可知: ,則:,則: )而)而 , 故由公式知故由公式知 即迭代法恒收斂。即迭代法恒收斂。211()kkkkxCxCxCxC2220202()1kkkkkkkxCxCqqxCCxCxCq200()kkkxCxCxCxC2211kkkqxCq2222()1kkkkkqxCxC qCq00 x 001xCqxC()kxC k 例例4.2.74.2.7 求求 的近似值,要求的近似值,要求 終止迭代。終止迭代。 解解 取取 經(jīng)經(jīng)6 6次迭代后:次迭代后: , , ,故故 。 例例4.2.4.2.8 8 對給定正數(shù)對給定正數(shù) ,應(yīng)用,應(yīng)用NewtonNewton法求解法求解 ,由此

54、式可導(dǎo)出求,由此式可導(dǎo)出求 而不用除法的計(jì)算程序:而不用除法的計(jì)算程序: 。 這個(gè)算法對于沒有設(shè)置除法操作的電子計(jì)算機(jī)是有用的??梢宰C明,這個(gè)算法對于沒有設(shè)置除法操作的電子計(jì)算機(jī)是有用的??梢宰C明,此算法初值滿足此算法初值滿足 時(shí)是收斂的,這是因?yàn)椋簳r(shí)是收斂的,這是因?yàn)椋?即:即: ,令,令 ,有遞推公式:有遞推公式: ,反復(fù)遞推得:反復(fù)遞推得: 。 當(dāng)當(dāng) ,即,即 時(shí),有時(shí),有 即即 ,從,從而迭代法收斂。而迭代法收斂。 106110kkxx1110()2kkkxxx01.0 x 53.16227767x 63.16227766x 7650.1 10 xx103.16227766C10Cx1

55、C1(2)kkkxxCx020 xC21111(2)()kkkkxxCxC xCCC211(1)kkCxCx1kkrC x21kkrr20kkrr0011rCx020 xC0kr 1kxC4.2.9 Newton4.2.9 Newton下山法下山法 NewtonNewton法收斂性依賴于法收斂性依賴于 初值的選取,如果初值的選取,如果 偏離偏離 較遠(yuǎn),較遠(yuǎn),則則NewtonNewton法可能發(fā)散。法可能發(fā)散。 例如,對方程例如,對方程 。求在。求在 附近的一個(gè)根附近的一個(gè)根 。若取。若取初值初值 ,則由,則由NewtonNewton法:法: 計(jì)算得計(jì)算得 ,僅迭代,僅迭代3 3次即得有次即得有

56、6 6位有效數(shù)字的近似值位有效數(shù)字的近似值 。但若取初值。但若取初值 則由同一則由同一NewtonNewton公公式計(jì)算得式計(jì)算得 ,這反而比,這反而比 更遠(yuǎn)離所求根更遠(yuǎn)離所求根 ,迭代迭代發(fā)散。為防止發(fā)散,對迭代過程加一下降要求:發(fā)散。為防止發(fā)散,對迭代過程加一下降要求: 滿足這項(xiàng)要求的算法稱為滿足這項(xiàng)要求的算法稱為下山法下山法。0 x0 x*x31 0 xx 1.5x *x01.5x 312131kkkkkxxxxx1231.34783,1.32520,1.32472xxx3x00.6x 117.9x 00.6x*1.32472x 1()()kkfxfx 將將NewtonNewton法與下

57、山法結(jié)合,即在下山法保證函數(shù)下降法與下山法結(jié)合,即在下山法保證函數(shù)下降條件下,用條件下,用NewtonNewton法加速收斂。為此,可將法加速收斂。為此,可將NewtonNewton計(jì)計(jì) 算結(jié)果算結(jié)果 與每一步近似值與每一步近似值 作加權(quán)均:作加權(quán)均: ,其中,其中 ( )稱稱為為下山下山 因子因子。選擇下山因子。選擇下山因子 以以保證下降性保證下降性。 的的選擇方法選擇方法是:由是:由 反復(fù)減半的試探法,若能反復(fù)減半的試探法,若能找到找到 使下降性成立,則下山成功,否則下山失敗使下降性成立,則下山成功,否則下山失敗,改改變初值變初值 重新開始。重新開始。 1()()kkkkf xxxfxkx

58、11(1)kkkxxx0110 x4.2.10 4.2.10 弦截法與拋物法弦截法與拋物法 NewtonNewton法法 每迭代一次計(jì)算函數(shù)值每迭代一次計(jì)算函數(shù)值 ,導(dǎo),導(dǎo)數(shù)值數(shù)值 各一次,當(dāng)各一次,當(dāng) 函數(shù)本身比較復(fù)雜時(shí),求導(dǎo)數(shù)值更加困難。函數(shù)本身比較復(fù)雜時(shí),求導(dǎo)數(shù)值更加困難。 下面方法多利用以前各次計(jì)算的函數(shù)值下面方法多利用以前各次計(jì)算的函數(shù)值 來回避導(dǎo)來回避導(dǎo)數(shù)值數(shù)值 的計(jì)算,導(dǎo)出這種求根方法的基本原理是插值法。的計(jì)算,導(dǎo)出這種求根方法的基本原理是插值法。 設(shè)設(shè) 是是 的一組近似值,利用對應(yīng)的函數(shù)的一組近似值,利用對應(yīng)的函數(shù)值值 ,構(gòu)造插值多項(xiàng)式,構(gòu)造插值多項(xiàng)式 ,適當(dāng)選取,適當(dāng)選取 的

59、一個(gè)根作為的一個(gè)根作為 的新的近似根的新的近似根 。這樣就確定了一個(gè)迭代。這樣就確定了一個(gè)迭代過程,記迭代函數(shù)為過程,記迭代函數(shù)為 ,則,則 ,下面具體考,下面具體考察察 (弦截法),(弦截法), (拋物法)兩種情形。(拋物法)兩種情形。 1()()kkkkf xxxfx( )kf x( )kf xf1( ), (),kkf xf x()kfx1,kkk rx xx( ) 0f x 1( ), (), ()kkk rf xf xf x( )rp x( ) 0rp x ( )0f x 1kxg11( ,)kkkk rxg x xx1r 2r 4.2.11 4.2.11 弦截法弦截法 設(shè)設(shè) 為為

60、的近似根,過點(diǎn)的近似根,過點(diǎn) , 構(gòu)造構(gòu)造一次插值多項(xiàng)式一次插值多項(xiàng)式 ,并用,并用 的根作為的根作為 的新的近似的新的近似根根 。由于。由于 則由則由 可得:可得: 另外,公式另外,公式(4.2.9)(4.2.9)也可以用導(dǎo)數(shù)也可以用導(dǎo)數(shù) 的差商的差商 近近似取代似取代NewtonNewton公式中的公式中的 ,同樣得公式,同樣得公式 。1,kkxx( ) 0f x ( , ( )kkx f x11(, ()kkxf x1( )p x1( ) 0p x ( )0f x 1kx111()()( )()() kkkkkkf xf xp xf xxxxx(4.2.8)1( )0p x 111()(

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論