計(jì)算方法教案

第1章誤差分析與數(shù)值計(jì)算3§1.1引言3§1.2絕對(duì)誤差與相對(duì)誤差、有效數(shù)字9§1.3近似數(shù)的簡(jiǎn)單算術(shù)運(yùn)算12§1.4數(shù)值計(jì)算中誤差分析的一些原則13第2章非線性方程（組）的近似解法15§2.1引言15§2.2根的隔離16§2.3對(duì)分法16§2.3對(duì)分法17§2.4迭代法19§2.6弦截法21§2.6弦截法22§1.7用牛頓法解方程組23本章小結(jié)25第3章線性方程組的解法26§3.1引言26§3.2高斯消去法28§3.3矩陣的LU分解31§3.4對(duì)稱(chēng)矩陣的LDLT分解32§3.5線性方程組解的可靠性33§3.6簡(jiǎn)單迭代法34本章小結(jié)43第4章矩陣特征值與特征向量的計(jì)算44§4.1引言44§4.2冪法和反冪法45§4.3雅可比方法46§4.4QR方法*51本章小結(jié)52第5章插值與擬合53§5.1引言53§5.2插值多項(xiàng)式的存在和唯一性54§5.3拉格朗日插值多項(xiàng)式55§5.4均差插值公式57§5.5差分等距結(jié)點(diǎn)插值公式59§5.6愛(ài)爾米特插值公式61§5.7分段低次插值62§5.8三次樣條函數(shù)63§5.9曲線擬合的最小二乘法67本章小結(jié)70第6章數(shù)值積分和數(shù)值微分71§6.1引言71§6.2牛頓一科特斯型積分公式72§6.3復(fù)合求積公式74§6.4龍貝格求積公式77§6.5高斯求積公式78§6.6二重積分的數(shù)值積分法80§6.7數(shù)值微分81本章小結(jié)83第7章常微分方程的數(shù)值解法84§7.1引言84§7.2歐拉法和改進(jìn)的歐拉法85§7.3龍格-庫(kù)塔方法86§7.4線性多步法89§7.5算法的穩(wěn)定性與收斂性91§7.6微分方程組和高階微分方程解法92本章小結(jié)94第1章誤差分析與數(shù)值計(jì)算§1.1引言1、課程任務(wù)和目的：在第七屆國(guó)際軟件工程學(xué)術(shù)會(huì)議上，“計(jì)算方法”被列入應(yīng)用方法學(xué)的研究領(lǐng)域，強(qiáng)調(diào)了計(jì)算方法的研究應(yīng)用與軟件方法學(xué)的研究密切結(jié)合。這就說(shuō)明了計(jì)算方法與軟件之間的聯(lián)系以及在應(yīng)用軟件研制中的地位與作用，計(jì)算方法是研究各種數(shù)學(xué)問(wèn)題求解的數(shù)值計(jì)算方法。在計(jì)算機(jī)成為數(shù)值計(jì)算的主要工具的今天，則要求研究適合于計(jì)算機(jī)使用的數(shù)值計(jì)算方法。計(jì)算方法就是研究用計(jì)算機(jī)解決數(shù)學(xué)問(wèn)題的數(shù)值方法及其理論，它的內(nèi)容包括函數(shù)的數(shù)值逼近、數(shù)值微分與數(shù)值積分、非線性方程值解、線性方程組數(shù)值解、常微和偏微數(shù)值解等，即都是以數(shù)學(xué)問(wèn)題為研究對(duì)象的。因此，計(jì)算方法是數(shù)學(xué)的一個(gè)分支，只是它不象純數(shù)學(xué)那樣只研究數(shù)學(xué)本身的理論，是把理論與計(jì)算緊密結(jié)合，著重研究數(shù)學(xué)問(wèn)題的數(shù)值方法及其理論，計(jì)算方法是計(jì)算機(jī)應(yīng)用和軟件研制開(kāi)發(fā)的重要組成部分，通過(guò)本課程的學(xué)習(xí)和上機(jī)實(shí)習(xí)，使學(xué)生掌握利用計(jì)算機(jī)進(jìn)行科學(xué)計(jì)算的基本理論和基本方法，并且學(xué)會(huì)將基本理論和基本方法應(yīng)用于軟件開(kāi)發(fā)以及軟件研制。2、本課程基本要求（1）掌握方法的基本原理和思想。（2）掌握方法處理的技巧及與計(jì)算機(jī)的結(jié)合。（3）掌握誤差分析，收斂性及穩(wěn)定性的基本理論。（4）學(xué)會(huì)進(jìn)行可靠的理論分析，對(duì)近似計(jì)算要確保精度要求，要進(jìn)行誤差分析。（5）通過(guò)例子，學(xué)習(xí)使用各種計(jì)算方法解決實(shí)際計(jì)算問(wèn)題。（6）通過(guò)上機(jī)實(shí)踐，能編寫(xiě)算法和實(shí)現(xiàn)算法。（7）掌握數(shù)值計(jì)算中一些最基本、最常用的計(jì)算方法和算法。3、本課程與各課程的關(guān)系：由于本課內(nèi)容包括了微積分、代數(shù)、常微分方程的數(shù)值方法，學(xué)生必須掌握這幾門(mén)課的基本內(nèi)容才能學(xué)好這一課程，同時(shí)，學(xué)習(xí)此課程還必須具備計(jì)算機(jī)系統(tǒng)的初步知識(shí)，掌握一門(mén)常用的高級(jí)語(yǔ)言，如：BASIC、PASCAL、C語(yǔ)言等，并須具備一定的編程能力。4、本課程的特點(diǎn)：（1）面向計(jì)算機(jī)，要根據(jù)計(jì)算機(jī)特點(diǎn)提供實(shí)際可行的有效算法。即算法只能包括加、減、乘、除運(yùn)算和邏輯運(yùn)算，是計(jì)算機(jī)能直接處理的。（2）有可靠的理論分析，能任意逼近并達(dá)到精度要求，對(duì)近似算法要保證收斂性和數(shù)值穩(wěn)定性，還要對(duì)誤差進(jìn)行分析，而且都是建立在相應(yīng)數(shù)學(xué)理論基礎(chǔ)上的。（3）有好的計(jì)算復(fù)雜性。時(shí)間復(fù)雜性好是指節(jié)省時(shí)間；空間復(fù)雜性好是指節(jié)省存儲(chǔ)量。這也是建立算法時(shí)要研究的問(wèn)題，因?yàn)樗P(guān)系到算法能否在計(jì)算機(jī)上完成。（4）要有數(shù)值實(shí)驗(yàn)。即任何一種算法除了從理論上要滿足上述三點(diǎn)外，還要通過(guò)數(shù)值實(shí)驗(yàn)證明是行之有效的。計(jì)算方法最基本的立足點(diǎn)是容許誤差，在誤差容許的范圍內(nèi)對(duì)某一數(shù)學(xué)問(wèn)題進(jìn)行近似計(jì)算，得到能滿足要求的近似結(jié)果?，F(xiàn)實(shí)世界中誤差是普遍存在的，由于世界上沒(méi)有絕對(duì)精確的量具（絕對(duì)精確的量具是沒(méi)有刻度的），因此人類(lèi)通過(guò)量具采集的數(shù)據(jù)都是近似值，另一方面，我們的生產(chǎn)、實(shí)驗(yàn)工具都不是絕對(duì)精確的，這就使得人類(lèi)在生產(chǎn)和科學(xué)實(shí)驗(yàn)中必需容許誤差。計(jì)算機(jī)的應(yīng)用可以分為二個(gè)方面，即數(shù)值計(jì)算和非數(shù)值計(jì)算。利用計(jì)算機(jī)進(jìn)行數(shù)值計(jì)算的過(guò)程如下圖所示：SKIPIF1<0在上圖中，計(jì)算方法的任務(wù)是：由建立的數(shù)學(xué)模型給出可編程并由計(jì)算機(jī)能完成的計(jì)算方法，然后編程和上機(jī)求解。由于計(jì)算方法是編程后可由計(jì)算機(jī)求解的近似計(jì)算方法，如何確保近似解的精度顯得尤為重要，必須深入討論有關(guān)誤差的基本概念和基本理論，為近似計(jì)算的精度分析打下基礎(chǔ)。1、誤差的來(lái)源（種類(lèi)）誤差的來(lái)源主要有以下四種（1）模型誤差：建立數(shù)學(xué)模型時(shí)的誤差。例如：在求重量的數(shù)學(xué)模型G=m*g中，重量G不是僅與質(zhì)量和重力加速度有關(guān)，它還與溫度、測(cè)量地點(diǎn)的海拔、地層結(jié)構(gòu)等眾多因素有關(guān)，為了使模型較為簡(jiǎn)單和實(shí)用，采用抓住主要矛盾的方法，去掉了大量對(duì)重量影響不大的次要因素，建立了上述重量的近似模型，由此產(chǎn)生了模型誤差。（2）觀測(cè)誤差：采集數(shù)據(jù)時(shí)的誤差。采集數(shù)據(jù)時(shí)，通常是依靠?jī)x器和量具，由于沒(méi)有絕對(duì)精確的儀器和量具，因此采集的數(shù)據(jù)有誤差，此誤差稱(chēng)為觀測(cè)誤差。（3）舍入誤差：由于計(jì)算機(jī)字長(zhǎng)有限而產(chǎn)生的誤差。硬件再發(fā)展，計(jì)算機(jī)的字長(zhǎng)總是有限的，在計(jì)算過(guò)程中，當(dāng)數(shù)據(jù)的長(zhǎng)度超過(guò)了計(jì)算機(jī)的字長(zhǎng)時(shí)，計(jì)算機(jī)就會(huì)進(jìn)行四舍五入，由此產(chǎn)生的誤差稱(chēng)為舍入誤差。（4）截?cái)嗾`差：無(wú)限形式的有限化而產(chǎn)生的誤差。在計(jì)算中有時(shí)會(huì)運(yùn)用無(wú)限形式的計(jì)算公式，例如臺(tái)勞公式：SKIPIF1<0顯然此公式無(wú)法進(jìn)行計(jì)算，因此必需根據(jù)實(shí)際需要，從某一項(xiàng)起將后面的各項(xiàng)截?cái)啵碨KIPIF1<0由此產(chǎn)生的誤差稱(chēng)為截?cái)嗾`差?！?.2絕對(duì)誤差與相對(duì)誤差、有效數(shù)字為描述方便，首先約定x*是精確值x的近似值。引入誤差的概念，其目的是為了衡量近似值x*的好壞。（1）絕對(duì)誤差：x*x由于精確值x通常無(wú)法確定，因此絕對(duì)誤差無(wú)法計(jì)算，由此引入絕對(duì)誤差限的概念。絕對(duì)誤差限：絕對(duì)誤差的一個(gè)上界。即：若|x*x|e，則稱(chēng)e為x*的絕對(duì)誤差限。絕對(duì)誤差限的性質(zhì)是：A.不唯一這是因?yàn)閨x*x|的上界是不唯一的。B.可確定只要我們對(duì)x*的實(shí)際背景有一定的了解，就不難確定|x*x|的上界。例如，x*表示身高，則|x*x|的上界可為3米。當(dāng)x*是你求出的，那么為了說(shuō)明你的工作認(rèn)真，你一定會(huì)將|x*x|的上界估計(jì)得盡量小，因此在這種意義上絕對(duì)誤差限可用來(lái)衡量x*的好壞。由于絕對(duì)誤差限沒(méi)有考慮問(wèn)題的規(guī)模，因此有時(shí)它也不能衡量x*的好壞。例如：x是地球與太陽(yáng)的距離，y是分子中二個(gè)原子間的距離，若|x*x|1公里，|y*y|1厘米，則并不能說(shuō)y*比（2）相對(duì)誤差：(x*x)/x*相對(duì)誤差限：相對(duì)誤差絕對(duì)值的一個(gè)上界。3、有效數(shù)字這里我們必須搞清楚什么是有效數(shù)字以及如何確定x*有幾位有效數(shù)字。（1）有效數(shù)字的定義若|x*-x|<x*的某一位的半個(gè)單位，則稱(chēng)x*精確到這一位，并從這一位開(kāi)始，一直到前面第一個(gè)不為零的數(shù)都是x*的有效數(shù)字。此定義實(shí)際上定義了什么叫精確到某一位和什么叫有效數(shù)字。例如：若x*精確到小數(shù)點(diǎn)后第3位，即指|x*x|0.510-3。（2）有效數(shù)字的判定方法方法一：四舍五入此方法首先確定x*是由x的哪一位四舍五入產(chǎn)生的，然后從這一位的前一位開(kāi)始一直到前面第一個(gè)不為零的數(shù)都是x*的有效數(shù)字。例1若x=0.872596,x*=0.87，求x*的有效位數(shù)。解：SKIPIF1<0x*是由x的小數(shù)點(diǎn)后第三位四舍五入產(chǎn)生的，所以x*有二位有效數(shù)字。注意，方法一判定有效數(shù)字很簡(jiǎn)單，但有時(shí)會(huì)失效。例如，若x=0.272987x*=0.273102，此時(shí)無(wú)法用方法一確定x*的有效位數(shù)，原因是x*不是由x四舍五入產(chǎn)生的，在這種情況下，必須用有效數(shù)字的定義來(lái)確定x*的有效位數(shù)。即方法二：用定義此方法首先計(jì)算|x*x|，再判斷它小于等于x*的哪一位的半個(gè)單位，然后從近一位開(kāi)始，一直到第一個(gè)不為零的數(shù)都是有效數(shù)字。例2若x=0.62073，x*=0.6207，確定x*的有效位數(shù)。解：因?yàn)閨x*x|0.00030.5104，x*精確到小數(shù)點(diǎn)后第4位，所以x*有四位有效數(shù)字。例3若x=0.080199，x*=0.802，確定x*的有效位數(shù)。解：因?yàn)閨x*x|=0.000010.5105，所以SKIPIF1<00.5103，推出x*有三位有效數(shù)字。例4若x=6.28936，x*=7.3132，確定x*的有效位數(shù)。解：|x*x|=0.023570.5101，所以x*有二位有效數(shù)字?！?.3近似數(shù)的簡(jiǎn)單算術(shù)運(yùn)算§1.4數(shù)值計(jì)算中誤差分析的一些原則為保證計(jì)算結(jié)果的高精度，在進(jìn)行數(shù)值計(jì)算時(shí)應(yīng)遵循下述幾個(gè)原則。（1）在進(jìn)行除法時(shí)，要避免除數(shù)的絕對(duì)值<<被除數(shù)的絕對(duì)值。①為什么要“避免”？若不“避免”，則除出的結(jié)果很大，由于計(jì)算機(jī)字長(zhǎng)有限，它裝不下，因此會(huì)進(jìn)行四舍五入，一個(gè)很大的數(shù)進(jìn)行四舍五入時(shí)舍去的部分也會(huì)很大，這會(huì)使舍入誤差變大。②怎樣“避免”？因?yàn)橛脩糁魂P(guān)心最后的計(jì)算結(jié)果，當(dāng)中間計(jì)算過(guò)程中出現(xiàn)了除數(shù)的絕對(duì)值<<被除數(shù)的絕對(duì)值時(shí)，就應(yīng)該換一種計(jì)算方法，以避免這種情況的發(fā)生，以后我們將會(huì)針對(duì)具體的計(jì)算問(wèn)題來(lái)討論“避免”的方法。（2）在進(jìn)行減法時(shí)，要避免二個(gè)相近的數(shù)相減。①為什么要“避免”？若不“避免”，就可能失去大量的有效數(shù)字，例如：若a=30001和b=30000都有五位有效數(shù)字，因?yàn)閍-b=1，所以結(jié)果至多有1位有效數(shù)字。②怎么“避免”？“避免”的思路與第1個(gè)原則中“避免”的思路相同，須針對(duì)具體計(jì)算問(wèn)題來(lái)討論。（3）要防止“大數(shù)吃小數(shù)”①什么是“大數(shù)吃小數(shù)”？我們用一個(gè)例子為說(shuō)明。計(jì)算8756294874SKIPIF1<0，其中n=1020，0<ai<106。此題是一個(gè)很大的數(shù)與很多很小的數(shù)相加，若采用將大數(shù)依次與a1,a2,,an相加，由于計(jì)算機(jī)字長(zhǎng)有限，因此在與ai相加時(shí)會(huì)進(jìn)行四舍五入將ai舍去，這樣，最后的結(jié)果仍是大數(shù)，這就是大數(shù)將a1,a2,,an吃掉了。②為什么要“避免”？盡管每個(gè)小數(shù)都很小，但它們很多，可能它們的和比大數(shù)還大，而最后計(jì)算工結(jié)果為大數(shù)，顯然誤差可能很大。③怎樣“避免”？有的同學(xué)提出先將小數(shù)相加，然后再與大數(shù)相加，這個(gè)思路是對(duì)的，但有一個(gè)漏洞，因?yàn)樾?shù)相加到一定程度也會(huì)變成大數(shù)，它也開(kāi)始吃小數(shù)了。可以采取分部相加的方法解決。第2章非線性方程（組）的近似解法§2.1引言方程f(x)=0的解稱(chēng)為方程的根。也叫做函數(shù)f(x)的零點(diǎn)。方程求根大致包括三個(gè)問(wèn)題（1）方程有沒(méi)有根？如果有根，有幾個(gè)根？（2）哪里有根？求有根的區(qū)間，區(qū)間內(nèi)的任意一點(diǎn)作為根的近似值。（3）根的精確化，已知一個(gè)根的近似值后設(shè)法逐步把根精確化，直到足夠精確為止。本課程主要研究問(wèn)題（2）和（3）?！?.2根的隔離求方程f(x)=0的解的近似值時(shí)，首先要確定若干個(gè)區(qū)間，使每個(gè)區(qū)間內(nèi)只有的一個(gè)根，這個(gè)步驟稱(chēng)為根的隔離。對(duì)一般的方程，根的隔離有兩種方法（1）試值法。求出f(x)在若干點(diǎn)上的函數(shù)值，觀察函數(shù)值符號(hào)變化的情況，從而確定隔根區(qū)間。（2）作圖法。畫(huà)出y=f(x)的草圖，觀察曲線y=f(x)與x軸交點(diǎn)的大致位置，從而確定隔根區(qū)間。例1.2.1討論方程f(x)=2x3-4x2+4x+2=0-1012-10-5-1012-10-5051011.52-1-0.500.5§2.3對(duì)分法11.52-1-0.500.5設(shè)有方程f(x)=0在（ab）內(nèi)有且僅有一個(gè)根x*，這時(shí)有f(a)f(b)<0可用對(duì)分法求x*的近似值，方法如下（1）準(zhǔn)備：計(jì)算區(qū)間（ab）兩個(gè)端點(diǎn)的函數(shù)值f(a)，f(b)（2）對(duì)分：取c=(a+b)/2為（ab）的中點(diǎn)，計(jì)算f(c)

（3）判斷：如果f(c)=0，則c為f(x)=0的根，否則檢驗(yàn)： 若f(c)f(a)<0，則方程的根位于[ac]內(nèi)，用c代替b，若f(c)f(b)<0，則方程的根位于[cb]內(nèi)，用c代替a。

（4）檢驗(yàn)：若|b-a|<e(e為精度要求)此時(shí)計(jì)算結(jié)束x*=c，否則轉(zhuǎn)（2）。例1.3.1用對(duì)分法求方程f(x)=x3+2x-5=0在[12]內(nèi)的根，[e=10-5有根區(qū)間f=inline('x^3+2*x-5')f(1),f(2)f=inline('x^3+2*x-5')f(1),f(2)fplot(f,[12]),gridon1.00001.50001.25001.50001.25001.37501.31251.37501.31251.34381.32811.34381.32811.33591.32811.33201.32811.33011.32811.3291方程的解x=1.3286§2.4迭代法設(shè)有方程f(x)=0在[ab]上有且僅有一個(gè)根x*，可用迭代法求x*的近似值，方法如下（1）將方程f(x)=0寫(xiě)成迭代形式x=(x)

（2）在[ab]上任取一個(gè)初始值x0。

（3）計(jì)算x1=(x0)

（4）若|x1x0|<e(e為精度要求)，此時(shí)計(jì)算結(jié)束x*=x1，否則令x0=x1轉(zhuǎn)（3）。例1.4.1用迭代法解方程x=10x-2，x0=1分別采用迭代格式x=10x-2和x=log(x+2)，觀察兩個(gè)計(jì)算過(guò)程的區(qū)別。e=1e-3迭代過(guò)程:1.00000.47710.39390.37910.37640.3759迭代6次x=0.3759f=inline('f=inline('log10(x+2)')x=1x=f(x)例1.4.2用迭代法求方程f(x)=x3+2x-5=0的根，x0=1[SKIPIF1<0]。迭代過(guò)程:1.00001.44221.28371.34491.32201.33061.32741.32861.3281f=inline('(5-2*x)^(1/3)')x=1x=f(x)迭代9次x=1.3281§2.5牛頓迭代法f=inline('(5-2*x)^(1/3)')x=1x=f(x)牛頓法是解方程f(x)=0的重要方法，它也是一種迭代法。設(shè)有方程f(x)=0在[ab]上有且僅有一個(gè)根x*，可用牛頓法求x*的近似值，方法如下（1）求函數(shù)f(x)的導(dǎo)函數(shù)f(x)，牛頓法迭代公式為x=xf(x)/f(x)

（2）在[ab]上任取一個(gè)初始值x0。

（3）計(jì)算x1=x0f(x0)/f(x0)

（4）若|x1x0|<e(e為精度要求)，此時(shí)計(jì)算結(jié)束x*=x1，否則令x0=x例1.5.1用牛頓法解方程f(x)=x3-2x2-4x-7=0在[34]內(nèi)的根[x0=4迭代過(guò)程:4.00003.67863.63293.6320迭代4次x=3.6320f=inline('x^3-2*x^2-4*x-7')fd=inline('3*x^2-4*x-4')x=4x=x-f(x)/fd(x)§2.6弦截法f=inline('x^3-2*x^2-4*x-7')fd=inline('3*x^2-4*x-4')x=4x=x-f(x)/fd(x)弦截法也是一種是解方程f(x)=0的迭代法，它的特點(diǎn)是不需要計(jì)算f(x)的函數(shù)f(x)，且收斂速度也相當(dāng)快，是工程計(jì)算中常用的算法之一。設(shè)有方程f(x)=0在[ab]上有且僅有一個(gè)根x*，可用弦截法求x*的近似值，方法如下（1）求函數(shù)f(x)在區(qū)間[ab]的兩個(gè)端點(diǎn)的函數(shù)值f(x0),f(x1)，其中a=x0,b=x1（2）計(jì)算x2=x1f(x1)[x1x0]/[f(x1)f(x0)]

（3）若|x2x1|<e(e為精度要求)，此時(shí)計(jì)算結(jié)束x*=x2，否則令x0=x1x1=x2轉(zhuǎn)（§1.7用牛頓法解方程組設(shè)有非線性方程組u(x,y)=0，v(x,y)=0，在(xn,yn)按臺(tái)勞級(jí)數(shù)展開(kāi)，取展開(kāi)式的第1，2項(xiàng)得到SKIPIF1<0SKIPIF1<0其中(xn,yn)是根的第n次近似值，如果Jn0方程組的第n+1次近似值(xn+1,yn+1)可用以下公式計(jì)算SKIPIF1<0例1.7.1用牛頓法解方程組u(x,y)=x3+y34=0，v(x,y)=x4+y23=0迭代初值x0=1,y0=1.4。（u/x=3x2u/y=3y2v/x=4x3v/y=2y例1.7.2用牛頓法解方程組u(x,y)=2x3y21=0，v(x,y)=xy3y4=0迭代初值x0=1.2,y0=1.7。（u/x=6x2u/y=2x v/x=6y3v/y=2xy21）例1.7.3用牛頓法解方程組u(x,y)=xcos(y)=0，v(x,y)=ysin(x)=0迭代初值x0=0,y0=0。（u/x=1 u/y=sin(y)v/x=cos(x) v/y=1）本章小結(jié)為了比較各種迭代方法的收斂速度，我們引入收斂階的概念。設(shè)迭代過(guò)程xn+1=(xn)收斂于方程x=(x)的根x*，令en=xnx*，en稱(chēng)為迭代誤差，如果存在實(shí)數(shù)P1和非零常數(shù)K，使得SKIPIF1<0，則稱(chēng)該迭代過(guò)程為P階收斂的。P=1稱(chēng)為線性收斂，P>1稱(chēng)為超線性收斂，P=2稱(chēng)為平方收斂，顯然P越大，迭代過(guò)程收斂的越快?？梢宰C明當(dāng)x*是方程f(x)=0的單根時(shí)，牛頓法是平方收斂的。當(dāng)x*是方程f(x)=0的重根時(shí)，牛頓法僅為線性收斂。弦截法的收斂階P=1.618。對(duì)分法的收斂速度與公比為1/2的等比級(jí)數(shù)相同。牛頓法：收斂速度最快，但要計(jì)算f(x)的導(dǎo)函數(shù)，計(jì)算量大，有發(fā)散問(wèn)題。弦截法：收斂速度次之，不需要計(jì)算f(x)的導(dǎo)函數(shù)計(jì)算量比牛頓法小，有發(fā)散問(wèn)題。對(duì)分法：收斂速度最慢，但簡(jiǎn)單有效，不存在發(fā)散問(wèn)題。它一定收斂到有根區(qū)間[ab]內(nèi)的某個(gè)根。第3章線性方程組的解法§3.1引言在科學(xué)實(shí)驗(yàn)和工程設(shè)計(jì)中，經(jīng)常用到解線性方程組的問(wèn)題。本章討論用計(jì)算機(jī)求解線性方程組的兩類(lèi)主要方法：直接法和迭代法。解線性方程組的一般表達(dá)式SKIPIF1<0根據(jù)矩陣的性質(zhì)可以寫(xiě)成SKIPIF1<0

簡(jiǎn)記為Ax=b其中SKIPIF1<0方程組Ax=b有唯一解的充分必要條件是|A|0。我們只討論這種情況下的解法。解線性方程組的方法可以分為兩類(lèi)：一類(lèi)是直接法，它只包含有限次的四則運(yùn)算，在每次運(yùn)算都無(wú)舍入誤差的情況下，所得到的是方程組的準(zhǔn)確解。由于實(shí)際計(jì)算中總是有舍入誤差，所以實(shí)際得到的也是近似解。令一類(lèi)是迭代法，它首先選擇一組初始值，再運(yùn)用同樣的計(jì)算步驟，重復(fù)計(jì)算，得到近似解。由于這類(lèi)方法中出現(xiàn)了極限過(guò)程，必須研究迭代過(guò)程的收斂性。本章主要介紹：直接法中的高斯消去法和主元高斯消去法。迭代法中的簡(jiǎn)單迭代法和塞德?tīng)枂蔚??！?.2高斯消去法以n=4為例說(shuō)明高斯消去法的計(jì)算過(guò)程，設(shè)有線性方程組SKIPIF1<0SKIPIF1<0SKIPIF1<0SKIPIF1<0經(jīng)過(guò)3次消元步驟，得到以上形式。從最后一個(gè)方程中解出x4,依此回代得到方程組的全部解。6x1+3x2+2x3=6

例2.2.3用高斯消去法解方程組10x1+5x2+6x3=0

8x1+5x2+3x3=0方程組的增廣矩陣[A|b]6326105608530消元6.00003.00002.00006.0000002.6667-10.000001.00000.3333-8.0000方程組系數(shù)矩陣主對(duì)角線元素為零，消元過(guò)程無(wú)法進(jìn)行!例2.2.4用列主元高斯消去法解例2.2.3中的方程組。方程組的增廣矩陣[A|b]6326105608530選主元1056063268530消元10.00005.00006.0000000-1.60006.000001.0000-1.80000選主元10.00005.00006.0000001.0000-1.8000000-1.60006.0000消元10.00005.00006.0000001.0000-1.8000000-1.60006.0000回代得到方程組的解5.6250-6.7500-3.7500§3.3矩陣的LU分解§3.4對(duì)稱(chēng)矩陣的LDLT分解§3.5線性方程組解的可靠性§3.6簡(jiǎn)單迭代法設(shè)有方程組Ax=b，變?yōu)榈问?，x=Mx+f，或x(k+1)=Mx(k)+f，任取初始值x(0)程迭代得到x(0)，x(1)，x(2)，，x(k)，若極限SKIPIF1<0存在，則x*就是原方程組的解。以n=4為例SKIPIF1<0x(k+1) M x(k)f寫(xiě)成分量形式SKIPIF1<0定理1若SKIPIF1<0，則簡(jiǎn)單迭代法對(duì)任意初始值x(0)和f都收斂。定理2若SKIPIF1<0，則簡(jiǎn)單迭代法對(duì)任意初始值x(0)和f都收斂。定理3迭代公式x(k+1)=Mx(k)+f，對(duì)任意初始值x(0)和f都收斂的充分必要條件是矩陣M的各個(gè)特征值的模都小于1?！?.6雅可比迭代法與高斯-塞德?tīng)柕ㄔ诤?jiǎn)單迭代法的基礎(chǔ)上作改進(jìn)x(k+1)=M1x(k+1)+M2x(k)+f，以n=4為例SKIPIF1<0x(k+1)M1x(k+1)M2x(k)f寫(xiě)成分量形式SKIPIF1<0定理1若SKIPIF1<0，則塞德?tīng)柕▽?duì)任意初始值x(0)和f都收斂。定理2若SKIPIF1<0，則塞德?tīng)柕▽?duì)任意初始值x(0)和f都收斂。定理3迭代公式x(k+1)=M1x(k+1)+M2x(k)+f，對(duì)任意初始值x(0)和f都收斂的充分必要條件是矩陣(I-M1)-1M松弛法（SuccessiveOverRelaxationMethod）x(k+1)=x(k)+(b-Ax(k))稱(chēng)為松弛因子，>1超松弛法，>1超松弛法，>1低松弛法。定理3松弛法對(duì)任意初始值x(0)和f都收斂的必要條件是0<<2。例2.6.1分別用雅可比迭代法和塞德?tīng)柕ń夥匠探MSKIPIF1<0誤差e<10-3雅可比迭代法迭代公式x(k+1)=Mx(k)+f,寫(xiě)成分量形式SKIPIF1<0初始值k=0(000),迭代過(guò)程x1(k) x2(k) x3(k) x1(k) x2(k) x3(k)-2.40005.00000.3000-4.00023.00311.9999-4.46124.24952.2802-4.00123.00002.0010-4.55582.74462.4671-4.00022.99922.0002-3.99132.62752.0345-3.99972.99981.9998-3.85792.98491.8865M=[0-0.4-0.2;0.250-0.5;-0.20.30]f=M=[0-0.4-0.2;0.250-0.5;-0.20.30]f=[-2.4;5;0.3]x=[0;0;0]x=M*x+f-4.03033.02372.0219-4.01382.98152.0132-3.99522.99001.9972-3.99543.00261.9960塞德?tīng)柕ǖ絰(k+1)=M1x(k+1)+M2x(k)+f寫(xiě)成分量形式SKIPIF1<0clearM1=[0-0.4-0.2;00-0.5;000]M2=[000;0.2500;-0.20.30]f=[-2.4;5;0.3]x=[0clearM1=[0-0.4-0.2;00-0.5;000]M2=[000;0.2500;-0.20.30]f=[-2.4;5;0.3]x=[0;0;0]B=inv((eye(3)-M1))x=B*M2*x+B*fx1(k) x2(k) x3(k)-2.40005.00000.3000-4.40004.85000.3000-4.00403.07661.8667-3.99962.98712.0238-4.00003.00211.9961-4.00002.99972.0006-4.00003.00002.0000例2.6.2分別用雅可比迭代法和塞德?tīng)柕ń夥匠探MSKIPIF1<0誤差e<10-3（1）雅可比迭代法迭代公式x(k+1)=Mx(k)+fSKIPIF1<0取初始值k=0(000),迭代過(guò)程x1(k) x2(k) x3(k)1.00000-2.00005.0000-0.99600.0080-1.00805.00806.0080-1.00005.00006.0000-1.00005.00006.0000 分析，M的特征方程SKIPIF1<0（2）塞德?tīng)柕ǖ絰(k+1)=M1x(k+1)+M2x(k)+fSKIPIF1<0，不收斂分析，B=(I-M1)-1M2的特征方程SKIPIF1<0例2.6.3分別用雅可比代法和塞德?tīng)柕ń夥匠探MSKIPIF1<0誤差e<10-3（1）雅可比迭代法迭代公式x(k+1)=Mx(k)+fSKIPIF1<0分析，M的特征方程SKIPIF1<0，所以簡(jiǎn)單迭代法不收斂。塞德?tīng)柕ǖ絰(k+1)=M1x(k+1)+M2x(k)+fSKIPIF1<0分析，B=(I-M1)-1M2的特征方程SKIPIF1<0，收斂取初始值k=0(000),迭代過(guò)程x1(k) x2(k) x3(k)2.00001.000001.49800.2500-0.74802.24900.2495-1.43702.59380.4216-1.59392.58610.5039-1.54312.51960.5117-1.49902.49370.5027-1.49172.49450.4986-1.49682.49910.4988-1.50012.50060.4997-1.50062.50040.5001-1.5002本章小結(jié)本章討論了解線性方程組的直接解法和迭代解法。直接解法比較適用與系數(shù)矩陣稠密（既零元素較少）的中、小型線性方程組，但對(duì)系數(shù)矩陣是帶狀或近似帶狀的大型線性方程組也適用。直接解法中的列主元高斯消去法具有精度較高和省時(shí)的優(yōu)點(diǎn)，是計(jì)算機(jī)中常用的算法。迭代解法中主要介紹了雅可比迭代法、高斯-塞德?tīng)柕ê退沙诜?。迭代法具有?jì)算公式簡(jiǎn)單、程序設(shè)計(jì)容易、占用計(jì)算機(jī)內(nèi)存較少的優(yōu)點(diǎn)。適用于解大型稀疏矩陣（既零元素較多）線性方程組。高斯-塞德?tīng)柕ㄊ窃谘趴杀鹊ǖ幕A(chǔ)上改進(jìn)得到，在很多情況下可以加快收斂速度，但它的收斂域與雅可比迭代法不同，因此不能互相取代。松弛法可以加速迭代過(guò)程的收斂速度，但要適當(dāng)選擇松弛因子（0<<2）。在選擇迭代法時(shí)，要特別注意檢驗(yàn)方法的收斂性問(wèn)題。第4章矩陣特征值與特征向量的計(jì)算§4.1引言求矩陣的特征值和特征向量，是代數(shù)計(jì)算中的重要問(wèn)題。在自然科學(xué)和工程中的許多問(wèn)題，例如電磁振蕩、橋梁的振動(dòng)，機(jī)械振動(dòng)等都可以歸結(jié)為求矩陣的特征值和特征向量問(wèn)題。矩陣A的特征值和特征向量是指，如果數(shù)和非零列向量x滿足關(guān)系式Ax=x，則數(shù)稱(chēng)為A的特征值非零列向量x稱(chēng)為A的與特征值對(duì)應(yīng)的特征向量。計(jì)算n階矩陣A的特征值，就是求特征方程|AI|=0的根i(i=1,2,,n)。齊次線性方程組(AiI)x=0的非零解xi，是i對(duì)應(yīng)的特征向量。本章討論一些在計(jì)算機(jī)上計(jì)算矩陣的特征值和特征向量的較為穩(wěn)定的數(shù)值算法?！?.2冪法和反冪法1、冪法：計(jì)算n階矩陣A的模最大的特征值（主特征值）及對(duì)應(yīng)的特征向量。任取n維列向量x(0)，用迭代公式x(k+1)=Ax(k)計(jì)算得到x(0)，x(1)，x(2)，設(shè)x(0)=a1v1+a2v2++anvn，因?yàn)锳vj=ivj所以x(1)=Ax(0)=a11v1+a22v2++annvnx(2)=Ax(1)=a112v1+a222v2++ann2vn一般地有x(k+1)=Ax(k)=a11kv1+a22kv2++annkvn=1k[a1v1+a2(2/1)kv2++an(n/1)kvn]當(dāng)k充分大時(shí)x(k+1)a11k+1v11x(k)向量x(k+1)與x(k)向近似地只差一個(gè)倍數(shù)，這個(gè)倍數(shù)就是模最大的特征值1。1、反冪法：Ax=x，A-1Ax=A-1x，A-1x=-1x，即A的特征值的倒數(shù)-1是A的逆矩陣A-1的特征值。用冪法求A-1的模最大的特征值，它的倒數(shù)就是A的模最小的特征值?！?.3雅可比方法對(duì)于2階方陣SKIPIF1<0SKIPIF1<0令SKIPIF1<0對(duì)于n階方陣（以n=3為例）SKIPIF1<0令SKIPIF1<0作變換矩陣SKIPIF1<0則有SKIPIF1<0令SKIPIF1<0作變換矩陣SKIPIF1<0則有SKIPIF1<0令SKIPIF1<0作變換矩陣SKIPIF1<0則有SKIPIF1<0一般地說(shuō)，令SKIPIF1<0，可以將A中的元素的aij和aji變?yōu)?。在實(shí)際計(jì)算中采用以下公式SKIPIF1<0例4.3.1 用雅可比求對(duì)稱(chēng)矩陣SKIPIF1<0的特征值和特征向量。消去第i行第j列的元素[ij]=[12]A->1.00000.0000-0.70710.00003.0000-0.7071-0.7071-0.70712.0000消去第i行第j列的元素[ij]=[13]A->0.6340-0.32510.0000-0.32513.0000-0.62800.0000-0.62802.3660消去第i行第j列的元素[ij]=[12]A->0.59010.0000-0.08390.00003.0438-0.6223-0.0839-0.62232.3660消去第i行第j列的元素[ij]=[13]A->0.5862-0.02930.0000-0.02933.0438-0.62160.0000-0.62162.3700消去第i行第j列的元素[ij]=[23]A->0.5862-0.0252-0.0150-0.02523.41400.0000-0.01500.00001.9998消去第i行第j列的元素[ij]=[12]A->0.58590.0000-0.01500.00003.41420.0001-0.01500.00011.9998§4.4QR方法*QR方法是求一般矩陣A的全部特征值和特征向量的一種迭代方法。其基本思路是利用矩陣A的QR分解，通過(guò)迭代格式SKIPIF1<0將A化為相似的上三角矩陣（或分塊上三角矩陣），從而求出A的全部特征值和特征向量。例3.4.1 用QR方法求矩陣SKIPIF1<0的特征值和特征向量。特征向量-0.8165-0.6215-0.6760-0.4082-0.6215-0.67600.40820.4770-0.2935特征值 1.00000.69724.3028本章小結(jié)本章介紹了求矩陣的特征值和對(duì)應(yīng)的特征向量的幾種方法。冪法可以求出矩陣的主特征值和對(duì)應(yīng)的特征向量，優(yōu)點(diǎn)是算法簡(jiǎn)單，但當(dāng)|1/2|1時(shí)，收斂速度很慢。反冪法可以求出矩陣的模最小的特征值和對(duì)應(yīng)的特征向量。雅可比方法是利用一系列正交相似變換（即平面旋轉(zhuǎn)變換）把實(shí)對(duì)稱(chēng)矩陣A化為對(duì)角陣（近似），從而求出實(shí)對(duì)稱(chēng)矩陣全部特征值。QR方法是用鏡向反射陣將矩陣A作QR分解，是一種求矩陣的全部特征值的有效方法。第5章插值與擬合§5.1引言已知表格函數(shù)y=f(x)xix0x1x2xn-1xnf(xi)y0y1y2yn-1yn構(gòu)造一個(gè)公式p(x)近似地表示f(x)，解決這個(gè)問(wèn)題的方法有兩類(lèi)：一類(lèi)是插值法，另一類(lèi)是擬合法，又稱(chēng)為逼近法。已知函數(shù)y=f(x)在互異點(diǎn)x0 ,x1,x2,,xn-1,xn上的函數(shù)值y0,y1,y2,,yn-1,yn，構(gòu)造一個(gè)函數(shù)p(x)使得p(xi)=yi這樣的問(wèn)題稱(chēng)為插值問(wèn)題。y=f(x)稱(chēng)為被插值函數(shù)，[x0xn]稱(chēng)為插值區(qū)間，p(x)稱(chēng)為插值函數(shù)，x0 ,x1,x2,,xn-1,xn稱(chēng)為插值點(diǎn)，在插值區(qū)間內(nèi)部用p(x)代替f(x)稱(chēng)為內(nèi)插，在插值區(qū)間外部用p(x)代替f(x)稱(chēng)為外推，R(x)=f(x)-p(x)稱(chēng)為插值函數(shù)p(x)的誤差?！?.2插值多項(xiàng)式的存在和唯一性定理：給出n+1個(gè)插值點(diǎn)及函數(shù)值xix0x1x2xn-1xnf(xi)y0y1y2yn-1yn求一個(gè)n次多項(xiàng)式pn(x)=a0+a1x+a2x2++anxn(x0,y0)(x1,y1)(x2,y2)(x3,y(x0,y0)(x1,y1)(x2,y2)(x3,y3)(x4,y4)§5.3拉格朗日插值多項(xiàng)式1、給出2個(gè)插值點(diǎn)(x0,y0)，(x1,y1)可以得到一次多項(xiàng)式SKIPIF1<02、給出3個(gè)插值點(diǎn)(x0,y0)，(x1,y1)，(x2,y2)可以得到二次多項(xiàng)式SKIPIF1<0不難驗(yàn)證p2(x)滿足插值條件p2(x0)=y0p2(x1)=y1p2(x2)=y23、給出n+1個(gè)插值點(diǎn)(x0,y0)，(x1,y1),,(xn,yn)可以得到一個(gè)n次多項(xiàng)式SKIPIF1<0其中SKIPIF1<0例5.3.1按下列表格求y(-0.5)和y(0.5)的值。x|123 y|7解：插值多項(xiàng)式SKIPIF1<0L2(-0.5)=0.2500L2(0.5)=4例4.3.2按下列表格求y(2.5)、y(4.5)、y(5.5)的值。x|解：插值多項(xiàng)式L4(x)=-0.7917x4+9.25x3-37.21x2+60.75x-32L4(2.5)=0.9180,L4(4.5)=6.1323,L4(5.5)=-8.9637§5.4均差插值公式已知函數(shù)f(x)在互異點(diǎn)x0,x1,,xn上的值為f(x0),f(x1),,f(xn)稱(chēng)SKIPIF1<0為函數(shù)f(x)在點(diǎn)xi,xj處的一階均差。稱(chēng)SKIPIF1<0為函數(shù)f(x)在點(diǎn)xi,xj,xk處的二階均差。一般地稱(chēng)SKIPIF1<0為函數(shù)f(x)在點(diǎn)x0,x1,,xn上的n階均差。牛頓插值公式Nn(x)=f(x0)+(x-x0)f(x0,x1) +(x-x0)(x-x1)f(x0,x1,x2)+ +(x-x0)(x-x1)(x-xn-1)f(x0,x1,,xn)牛頓插值公式具有遞推關(guān)系Nk+1(x)=Nk(x)+(x-x0)(x-x1)(x-xk)f(x0,x1,,xk+1)新增加一個(gè)節(jié)點(diǎn)，只需要增加計(jì)算一項(xiàng)(x-x0)(x-x1)(x-xk)f(x0,x1,,xk+1)§5.5差分等距結(jié)點(diǎn)插值公式1、差分已知函數(shù)f(x)在等距節(jié)點(diǎn)x0,x1,,xn上的值xx0x+hx+2hx+nhf(x)y0y1y2yn稱(chēng)△yk=yk+1yk為函數(shù)f(x)在點(diǎn)xk處的一階差分。稱(chēng)△2yk=△yk+1△yk為函數(shù)f(x)在點(diǎn)xk處的二階差分。一般地稱(chēng)△myk=△m-1yk+1△m-1yk為函數(shù)f(x)在點(diǎn)xk處的m階差分。2、等距結(jié)點(diǎn)插值公式牛頓向前插值公式SKIPIF1<0牛頓向后插值公式SKIPIF1<0§5.6愛(ài)爾米特插值公式定理：給出n+1個(gè)插值點(diǎn)上的函數(shù)值和導(dǎo)數(shù)值xix0x1x2xn-1xnf(xi)y0y1y2yn-1ynf(xi)y0y1y2yn-1yn求一個(gè)2n+1次多項(xiàng)式H(x)滿足Hn(xi)=yi，Hn(xi)=yi(i=0,1,2,,n)。SKIPIF1<0其中SKIPIF1<0SKIPIF1<0§5.7分段低次插值§5.8三次樣條函數(shù)在xoy平面上給定n+1個(gè)點(diǎn)(x0,y0),(x1,y1),,(xn,yn)，構(gòu)造一個(gè)函數(shù)S(x)滿足以下條件（1）S(xi)=yi(i=0,1,2,,n)

（2）在區(qū)間(x0,xn)內(nèi)S(x)具有連續(xù)的二階導(dǎo)數(shù)（3）在每個(gè)子區(qū)間[xi1,xi]上S(x)（表達(dá)式是Si(x)）是一個(gè)三次多項(xiàng)式。滿足以上條件S(x)稱(chēng)為三次樣條插值多項(xiàng)式。三次樣條插值多項(xiàng)式的推導(dǎo)設(shè)在子區(qū)間[xi1,xi]上S(x)=Si(x)(i=1,2,,n)由條件（1）得Si(xi-1)=yi1，Si(xi)=yi設(shè)S(x)在節(jié)點(diǎn)xi1處的二階導(dǎo)數(shù)為Mi1，在節(jié)點(diǎn)xi處的二階導(dǎo)數(shù)為Mi，即Si(xi-1)=Mi1，Si(xi)=Mi，顯然Si(x)是x的線性函數(shù)，根據(jù)拉格朗日插值公式有SKIPIF1<0，記xi-xi1=hi有SKIPIF1<0將上式積分兩次得到SKIPIF1<0利用Si(xi-1)=yi1，Si(xi)=yi定出積分常數(shù)C1和C2SKIPIF1<0解此方程組得到SKIPIF1<0代入上式整理后得到SKIPIF1<0SKIPIF1<0SKIPIF1<0類(lèi)似地有SKIPIF1<0SKIPIF1<0因?yàn)镾i+1(xi)=Si(xi)所以有SKIPIF1<0整理后得到關(guān)于位知數(shù)M0,M1,,Mn的線性方程組aiMi-1+2Mi+biMi+1=di(i=1,2,其中SKIPIF1<0邊界條件：n-1個(gè)方程組，n+1個(gè)位知數(shù)，所有需要補(bǔ)充兩個(gè)條件，常見(jiàn)的是以下兩種（1）給定端點(diǎn)的二階導(dǎo)數(shù)M0=a，Mn=b，特別地，當(dāng)a=b=0時(shí)稱(chēng)為自然樣條。（2）給定端點(diǎn)的一階導(dǎo)數(shù)S(x0)=a，S(xn)=b，這時(shí)有SKIPIF1<0例5.7.1給定插值條件，和兩種邊界條件(1)m0=1,m3=0(2)M0=1,M3=0

x|0123

y|0000

m0=1,m3=0s1(x)=0.73333x3-1.7333x2+x0<x<1s2(x)=-0.2x3+1.0667x2-1.8x+0.933331<x<2s3(x)=0.066667x3-0.53333x(2)M0=1,M3=0s1(x)=-0.21111x3+0.5x2-0.28889x0<x<1s2(x)=0.055556x3-0.3x2+0.51111x-0.26667§5.9曲線擬合的最小二乘法給出表格函數(shù)xix0x1x2xn-1xnf(xi)y0y1y2yn-1yn求一個(gè)m(m<n)次多項(xiàng)式為pm(x)=a0+a1x+a2x2++amxm逼近這組數(shù)據(jù)。pm(x)的系數(shù)的確定：SKIPIF1<0即SKIPIF1<0(i=0,1,2,,n)SKIPIF1<0由微分學(xué)知，使(a0,a1,,am)達(dá)到極小值的a0,a1,,am滿足必要條件SKIPIF1<0(K=0,1,2,,m)寫(xiě)成分量形式SKIPIF1<0正規(guī)方程組解此方程組得到a0,a1,,am即得到所要的多項(xiàng)式。例5.8.1有數(shù)據(jù)表如下，分別用二次和三次多項(xiàng)式逼近這組數(shù)據(jù)。x|-3-2-10123 y|p1(x)=0.1786x+0.5714p2(x)=0.17857x2+0.17857x-0.14286p3(x)=3.2895*10-17x3+0.17857x2+0.17857x-0.14286一次多項(xiàng)式逼近誤差平方和2.8214，二次多項(xiàng)式逼近誤差平方和0.1429，三次多項(xiàng)式逼近誤差平方和0.1429。本章小結(jié)本章介紹了兩種構(gòu)造函數(shù)f(x)的逼近函數(shù)的方法—插值法和曲線擬合。關(guān)于插值法，討論的是多項(xiàng)式插值，它要求所構(gòu)造的多項(xiàng)式函數(shù)嚴(yán)格地通過(guò)給定的所有數(shù)據(jù)點(diǎn)。由Lagrange插值基函數(shù)的討論，導(dǎo)出了Lagrange插值公式及其余項(xiàng)公式。作為對(duì)Lagrange插值的改進(jìn)，建立了具有遞推性的牛頓插值公式。Hermite插值是一種在插值節(jié)點(diǎn)上插值函數(shù)與被插值函數(shù)不僅有相同的函數(shù)值，而且還有相同的一階導(dǎo)數(shù)值的多項(xiàng)式插值。由于高次插值會(huì)產(chǎn)生龍格現(xiàn)象，因此討論了分段插值法，重點(diǎn)介紹了三次樣條插值公式。關(guān)于曲線擬合，討論的是曲線擬合的最小二乘法，它不要求所構(gòu)造的逼近函數(shù)嚴(yán)格地通過(guò)給定的所有數(shù)據(jù)點(diǎn)，只是在多項(xiàng)式中，以使得殘差的平方和最小為標(biāo)準(zhǔn)，選擇多項(xiàng)式，以作為被逼近函數(shù)的近似替代。本章最后介紹了數(shù)值微分，其中5點(diǎn)公式是常用的。第6章數(shù)值積分和數(shù)值微分§6.1引言計(jì)算定積分的牛頓-萊布尼茲公式SKIPIF1<0，其中F(x)是函數(shù)f(x)的原函數(shù)。在實(shí)際應(yīng)用中（1）常遇到某些函數(shù)f(x)的原函數(shù)不能用初等函數(shù)表示，例如sin(x2)，cos(x2),sin(x)/x，1/log(x)等。（2）有些函數(shù)f(x)是用表格形式表示，無(wú)法得到它們的原函數(shù)。（3）有些函數(shù)f(x)的原函數(shù)十分復(fù)雜，不利于工程上的使用。因此，要研究計(jì)算定積分的近似方法：數(shù)值積分?！?.2牛頓一科特斯型積分公式由上一章知，任意函數(shù)f(x)可以用一個(gè)拉格朗日插值多項(xiàng)式n(x)近似表示，因此f(x)的定積分也可以用n(x)的定積分近似表示SKIPIF1<0以此為基礎(chǔ)得到的數(shù)值積分公式稱(chēng)為牛頓一科特斯型積分公式。SKIPIF1<0其中SKIPIF1<0取a=x0<x1<<xn=b為一組等距節(jié)點(diǎn)，令x=a+bt，h=(b-a)/n得到SKIPIF1<0記SKIPIF1<0，于是SKIPIF1<0所以得SKIPIF1<0對(duì)不同的n系數(shù)如下n牛頓一科特斯型積分公式系數(shù)余項(xiàng)（誤差）11/2，1/2 (梯形公式)(1/12)h3f(2)()a<21/6，4/6，1/6 （拋物線公式或辛卜生公式）(1/90)h5f(4)(31/8，3/8，3/8，1/8(3/80)h5f(4)(47/90，32/90，12/90，32/90，7/90 （柯特斯公式）(8/945)h7f(6)(519/288，75/288，50/288，50/288，75/288，19/288(275/12096)h7f(6)(641/840，216/840，27/840，272/840，27/840，216/840，41/840(9/1400)h3f(2)(§6.3復(fù)合求積公式1、復(fù)合梯形公式取a=x0<x1<<xn=b為一組等距節(jié)點(diǎn)將積分區(qū)間[ab]分為n個(gè)子區(qū)間[xkxk-1]xk–xk-1=(b-a)/n=h，在每個(gè)子區(qū)間[xkxk-1]上應(yīng)用梯形公式得：SKIPIF1<0于是得到復(fù)合梯形公式SKIPIF1<02、復(fù)合辛卜生公式在積分區(qū)間[ab]取點(diǎn)a=x0<x1<<x2n-1<x2n=b將[ab]分為2n個(gè)子區(qū)間，令x2k-2x2k=(b-a)/n=2h，在每個(gè)子區(qū)間[x2k-2x2k]上應(yīng)用辛卜生公式得：SKIPIF1<0（k=1,2,,n）于是得到復(fù)合辛卜生公式SKIPIF1<03、復(fù)合柯特斯公式在積分區(qū)間[ab]取點(diǎn)a=x0<x1<<x4n-1<x4n=b將[ab]分為4n個(gè)子區(qū)間，令x4k-4x4k=(b-a)/n=4h，在每個(gè)子區(qū)間[x4k-4x4k]上應(yīng)用柯特斯公式得：SKIPIF1<0（k=1,2,,n）得到復(fù)合柯特斯公式SKIPIF1<0（k=1,2,,n）4、步長(zhǎng)h的自動(dòng)選擇復(fù)合梯形公式|T2nTn|<3e（e表示允許誤差）其中Tn和T2n分別表示取n和2n時(shí)用復(fù)合梯形公式計(jì)算得到的積分近似值。（2）復(fù)合辛卜生公式|S2nSn|<15e（e表示允許誤差）其中Sn和S2n分別表示取n和2n時(shí)用復(fù)合辛卜生公式計(jì)算得到的積分近似值。（3）復(fù)合柯特斯公式|C2nCn|<63e（e表示允許誤差）其中Cn和C2n分別表示取n和2n時(shí)用復(fù)合柯特斯公式計(jì)算得到的積分近似值。§6.4龍貝格求積公式梯形公式與辛卜生公式之間的關(guān)系SKIPIF1<02、辛卜生公式與柯特斯公式之間的關(guān)系SKIPIF1<03、龍貝格積分公式SKIPIF1<0§6.5高斯求積公式一個(gè)求積公式，若對(duì)于任何次數(shù)不超過(guò)m的多項(xiàng)式都準(zhǔn)確成立，則稱(chēng)這個(gè)求積公式的代數(shù)精確度為m。定義：使插值求積公式SKIPIF1<0的代數(shù)精確度為2n-1的節(jié)點(diǎn)x1,x2,,xn稱(chēng)為高斯點(diǎn)，對(duì)應(yīng)的插值求積公式稱(chēng)為高斯求積公式。定理1：若x1,x2,,xn是高斯點(diǎn)，則SKIPIF1<0定理2：高斯求積公式的系數(shù)恒為正，且有SKIPIF1<0n節(jié)點(diǎn)xk(n)系數(shù)Ak(n)余項(xiàng)（誤差）積分區(qū)間[-11]102f(2)()/31<<120.5773503+0.577350311f(4)()/13530.5773503

0+0.57735035/98/9

5/9f(6)()/1575040.8611363

0.3399810

+0.3399810+0.86113630.3478548

0.6521452

0.6521452

0.3478548f(8)()/347287550.90617990.5384693

0

+0.5384693+0.90617990.2369269

0.4786287

0.5688889

0.4786287

0.2369269f(10)()/1237732650§6.6二重積分的數(shù)值積分法§6.7數(shù)值微分給出表格函數(shù)如下表所示，求f(x)在節(jié)點(diǎn)f(xi)處的導(dǎo)數(shù)的近似值。xix0x1x2xn-1xnf(xi)y0y1y2yn-1yn方法：用插值多項(xiàng)式的導(dǎo)數(shù)近似f(x)的導(dǎo)數(shù)。三點(diǎn)公式：SKIPIF1<0SKIPIF1<0三個(gè)相鄰節(jié)點(diǎn)的選法，一般是在所考察的節(jié)點(diǎn)兩側(cè)各選取一個(gè)節(jié)點(diǎn)，如果一側(cè)的無(wú)節(jié)點(diǎn)，則用另一側(cè)的節(jié)點(diǎn)補(bǔ)足。一階導(dǎo)數(shù)的