第一章:插值方法

寫在前面
插值問題講的是啥呢脱柱?實際上就是求原函數(shù)佣蓉,但是我們這里只有幾個離散點的值兆解,本章主要研究的就是不同的方法求出不同的插值多項式馆铁。

1.1 拉格朗日插值多項式

線性插值:線性插值是最簡單的插值函數(shù),就是兩點決定一條直線锅睛,用兩點式表示了而已埠巨。


那我們換一種表示方法,把兩點式中的兩個式子給他換個名字:

這個l?l?就是傳說中的線性插值基函數(shù)

接下來講一個貫穿整個插值函數(shù)的例題:


現(xiàn)在看這道題可能有點簡單现拒,到了后面再做的時候就明白這道題的霸道啦~


拋物插值:拋物插值就是線性插值的進階版辣垒,給出三個點,求出來的插值多項式就是所謂的拋物插值具练。
從線性插值中推理過來:


這里面的l?乍构,l?l? 展開后扛点,就是介個樣子:

那我們再把這個例題掏出來:


很明顯哥遮,這個結(jié)果比線性插值精度高一些~


一般的:
對于我們的l?:


那一般的插值基函數(shù)就是這樣的:
可以看到,分母上減數(shù)沒有l?陵究,被減數(shù)是l?眠饮,這是個用來記憶的規(guī)律~

所以根據(jù)基函數(shù),可以得到n階插值多項式的表達(dá)方式如下:

為了表示方便铜邮,引入w來做一個標(biāo)記仪召,w實際上就是分母:

那么插值多項式就可以表示為:

實際上吧,我覺得這里沒有太大必要這樣表示松蒜,但是為了后續(xù)的理解扔茅,姑且先記一下叭~繞了個小彎。

1.2 拉格朗日插值多項式余項

啥是余項?
通俗來講秸苗,余項就是誤差召娜,所以插值多項式的余項可以表示成:

其中f(x)是精確解,L?(x)是所求的多項式惊楼。

再進一步寫:

快看這里的w眼熟不玖瘸!就是上面那個展開奇奇怪怪的東西~

具體的證明不需要記秸讹,但是要記住,余項表達(dá)式只有在f(x)高階導(dǎo)數(shù)存在時才能用雅倒。
通常璃诀,我們求函數(shù)的n+1階導(dǎo)數(shù)max|f??1(x)| = M???,從而將誤差放縮:

所以蔑匣,我們直接拿來結(jié)果劣欢!

  • n=1時,即線性插值余項:
  • n=2時裁良,即拋物插值余項:

那我們再把上面那道例題掏出來氧秘,這里的運算比較復(fù)雜, 考試估計不會這樣趴久,但求誤差估計的方法很典型,看明白就好

這里方法有一個局限性搔确,就是必須要知道導(dǎo)函數(shù)的上界彼棍,屬于事前誤差估計,那如果上界不知道呢膳算?


事后誤差估計方法

事后誤差估計是個怎么回事兒呢座硕?通俗來講,就是多算一位涕蜂,分別把L?L???的式子算出來华匾,近似相等,可以得到結(jié)果和誤差:

1.2 差商

定義: 一階差商就是机隙,函數(shù)值之差比上自變量之差:


二階差商就是一階差商的差商:

那么一般地蜘拉,K階差商:

性質(zhì)

  • K階差商可以表示為f(x0)···f(xk)的線性組合
  • 差商與節(jié)點的排列次序無關(guān),具有對稱性

計算:使用差商表最方便

另外還有一個題型無法使用差商表:


這時候要想起來一個和計算有關(guān)的性質(zhì):

則有:

一般這樣的題型求導(dǎo)之后有鹿,導(dǎo)數(shù)值都是一個常數(shù)旭旭,可以直接計算出結(jié)果。

1.3 牛頓插值多項式

實際上牛頓和拉格朗日插值是等價的葱跋,拉格朗日插值有高度的對稱性持寄;牛頓插值多項式來自于差商,其意義在于具有承襲性娱俺,即增加一項可以從上一項推出來稍味。

定義式:

遞推式:

Newton插值余項

這個東西順便證明了差商的第三條性質(zhì)。

又雙叒叕掏出來上面那道例題:


在使用Newton插值多項式的時候荠卷,先根據(jù)變量值和函數(shù)值計算出差商表模庐,再結(jié)合公式帶入即可:

1.4 分段插值

龍格現(xiàn)象:所謂龍格現(xiàn)象,就是當(dāng)插值多項式的次數(shù)隨著節(jié)點個數(shù)增加時僵朗,有可能產(chǎn)生激烈的震蕩從而不符合原函數(shù)赖欣。
分段插值:分段插值就是將被插值函數(shù)分成一小段一小段屑彻,在每個小段里面逼近,從而達(dá)到比較好的效果顶吮。

1.4.1 分端Lagrange插值

分段線性插值:將一個區(qū)間化為n個小區(qū)間社牲,記h是所有區(qū)間長度的最大值,則Ih在[a,b]上連續(xù)悴了、存在且在每一段上都是線性多項式搏恤,即為分段線性差值函數(shù)

實際上就是在每個小區(qū)間上的折線湃交,可以用拉格朗日插值多項式表示:

重點要記一下余項的算法:

分段線性插值:

分段二次插值:

1.4.2 分段Hermit插值

為了克服拉格朗日插值中熟空,分段點處不可導(dǎo)的問題

這邊的證明太太太長了,我們記幾個關(guān)鍵的公式即可:

三次Hermit插值公式:
其中的四個式子就是三次Hermit插值基函數(shù)

還是用一道例題熟悉吧

方法一:基函數(shù)法

方法二:待定系數(shù)法
我覺得待定系數(shù)好理解一點emm

1.5 樣條函數(shù)

樣條函數(shù)的特點是搞莺。充分光滑息罗,即導(dǎo)數(shù)連續(xù);又有一定的間斷性才沧,即分段的特性迈喉。

三次樣條插值

計算三次樣條函數(shù)時,需要的邊界條件:

  • 端點處的一階導(dǎo)數(shù):
  • 端點處的二階導(dǎo)數(shù):
    當(dāng)二階導(dǎo)數(shù)的兩端點值都為0時温圆,稱為自然邊界條件挨摸,樣條函數(shù)稱為自然樣條函數(shù)
  • f(x)是周期函數(shù),且xn-x0是一個周期岁歉,則要求S(x)也是周期函數(shù):

整個例題理解一下:

這種題的解法圍繞定義下手得运,S(x)需滿足在作用域內(nèi)二階連續(xù)可導(dǎo),且一階和二階導(dǎo)函數(shù)連續(xù)锅移,帶入即可熔掺。


接下來講講三次樣條插值函數(shù)的計算方法

三轉(zhuǎn)角方法(題目中給的是端點一階導(dǎo)數(shù)):
具體的公式推導(dǎo)實在太麻煩,直接上干貨帆啃,先記幾個公式:

這是由分段Hermit推來的式子瞬女,放在這幫助理解

其中h是劃分每一段的長度,為簡便計算努潘,引入三個已知量:

這里的mi是樣條函數(shù)的一階導(dǎo)數(shù)诽偷,結(jié)合上面三個式子,就能得到mi

還是整個例題看看吧:


這里把區(qū)間分成三份疯坤,首先計算λ报慕、μg压怠,由公式代入可得:
這里不要被那個矩陣唬住了眠冈,想不明白為什么就帶入上面公式的最后一個求mi,通過λμ蜗顽、g解二元一次方程組布卡。


三彎矩方法(題目中給了端點的二階導(dǎo)數(shù))
Mi為S的二階導(dǎo)數(shù),則有:


再次引入λi雇盖、μi忿等、di

從而求解Mi:

再來個類似的例題:


同樣分為三份區(qū)間,帶入公式:

1.6 數(shù)據(jù)擬合的最小二乘法

這里就不過多證明崔挖,直接上例題尋找考點吧

解法過程如下:

  • 先描點畫一個大概的草圖贸街,判斷函數(shù)的次數(shù),本題可以看出是一次
  • 用冪函數(shù)擬合曲線狸相,即用g(x) = 1+x+x2+x3+····作為擬合多項式
  • 結(jié)合給出的離散數(shù)據(jù)寫出向量薛匪,然后兩兩做內(nèi)積

    其中系數(shù)分別為φ0φ0內(nèi)積=8,φ0φ1內(nèi)積=4脓鹃、φ1φ1內(nèi)積44逸尖,φ0f內(nèi)積=3.9,φ1f內(nèi)積=46

再來一個題型:


這個題型的主要方法就是瘸右,把解方程組冷溶,轉(zhuǎn)化為求G(x,y)的最小值,即每一項都為零的時候成立尊浓,進一步轉(zhuǎn)化成求偏導(dǎo):

最終求得的是近似解而非精確解。

習(xí)題

1.插值多項式的次數(shù)與插值節(jié)點的個數(shù)有關(guān)系纯衍。正確
2.若n+1個插值節(jié)點互不相同栋齿,則滿足插值條件且不大于n次的插值多項式唯一存在
3.拉格朗日插值基函數(shù)滿足插值條件襟诸,lk(xi)= a,(i=k); lk(xi)= b,(i≠k)瓦堵。則a=1b=0歌亲。
4.用拉格朗日插值法求插值多項式就是對應(yīng)節(jié)點xk的基函數(shù)lk(x)與相應(yīng)節(jié)點函數(shù)值yk乘積之和菇用。正確
5.插值多項式余項中的ζ與x無關(guān)。錯誤

1.Newton插值法的插值基函數(shù)既有與節(jié)點相關(guān)又有升冪的特點陷揪,從而改進了拉格朗日插值不具有繼承性的不足惋鸥。正確
2.差商f[x0,x1,x2] = f[x1,x0,x2] = f[x1,x2,x0]。正確
3.插值節(jié)點從x0到xn悍缠,對節(jié)點重新排列之后卦绣,對應(yīng)的Newton插值是否相等?相等
4.n+1個節(jié)點的拉格朗日插值多項式與牛頓插值多項式只是表現(xiàn)形式不同飞蚓,實質(zhì)上是等價的滤港。正確
5.Newton插值多項式中,每增加一個插值節(jié)點趴拧,所有的差商值都需要重新計算溅漾。錯誤

1.Runge現(xiàn)象產(chǎn)生的原因是插值節(jié)點不多山叮。錯誤
2.插值多項式的余項隨著節(jié)點的增多而在某些點可能產(chǎn)生激烈的震蕩
3.分段線性插值克服了高次插值多項式誤差可能產(chǎn)生震蕩的不足,但分段線性插值函數(shù)在整個插值區(qū)間上只能保證連續(xù)添履,但不連續(xù)可導(dǎo)屁倔。
4.三次樣條插值函數(shù)S(x)是分段函數(shù)。
5.三次樣條插值要求插值函數(shù)在整個插值區(qū)間上都是三階連續(xù)可導(dǎo)缝龄。錯誤
6.已知相同的離散數(shù)據(jù)汰现,插值和逼近(擬合)可獲得相同的函數(shù)表達(dá)式。錯誤
7.用最小二乘法進行數(shù)據(jù)擬合時叔壤,獲得的正則線性方程組的系數(shù)矩陣是對稱矩陣瞎饲。正確

書后題:


后記
這一章太難了,太難了炼绘,加油兄弟們

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末嗅战,一起剝皮案震驚了整個濱河市咽笼,隨后出現(xiàn)的幾起案子阀蒂,更是在濱河造成了極大的恐慌,老刑警劉巖玖翅,帶你破解...
    沈念sama閱讀 206,378評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件脚曾,死亡現(xiàn)場離奇詭異东且,居然都是意外死亡,警方通過查閱死者的電腦和手機本讥,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,356評論 2 382
  • 文/潘曉璐 我一進店門珊泳,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人拷沸,你說我怎么就攤上這事色查。” “怎么了撞芍?”我有些...
    開封第一講書人閱讀 152,702評論 0 342
  • 文/不壞的土叔 我叫張陵秧了,是天一觀的道長。 經(jīng)常有香客問我序无,道長验毡,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 55,259評論 1 279
  • 正文 為了忘掉前任帝嗡,我火速辦了婚禮米罚,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘丈探。我一直安慰自己录择,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 64,263評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著隘竭,像睡著了一般塘秦。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上动看,一...
    開封第一講書人閱讀 49,036評論 1 285
  • 那天尊剔,我揣著相機與錄音,去河邊找鬼菱皆。 笑死须误,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的仇轻。 我是一名探鬼主播京痢,決...
    沈念sama閱讀 38,349評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼篷店!你這毒婦竟也來了祭椰?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 36,979評論 0 259
  • 序言:老撾萬榮一對情侶失蹤疲陕,失蹤者是張志新(化名)和其女友劉穎方淤,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體蹄殃,經(jīng)...
    沈念sama閱讀 43,469評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡携茂,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,938評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了诅岩。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片邑蒋。...
    茶點故事閱讀 38,059評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖按厘,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情钱慢,我是刑警寧澤逮京,帶...
    沈念sama閱讀 33,703評論 4 323
  • 正文 年R本政府宣布,位于F島的核電站束莫,受9級特大地震影響懒棉,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜览绿,卻給世界環(huán)境...
    茶點故事閱讀 39,257評論 3 307
  • 文/蒙蒙 一策严、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧饿敲,春花似錦妻导、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,262評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽术浪。三九已至,卻和暖如春寿酌,著一層夾襖步出監(jiān)牢的瞬間胰苏,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,485評論 1 262
  • 我被黑心中介騙來泰國打工醇疼, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留硕并,地道東北人。 一個月前我還...
    沈念sama閱讀 45,501評論 2 354
  • 正文 我出身青樓秧荆,卻偏偏與公主長得像倔毙,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子辰如,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,792評論 2 345