三次樣條插值課件

引例2.7.1三次樣條插值函數(shù)的概念一背景二、樣條函數(shù)的定義

例2.13

定理2.8（3次樣條插值函數(shù)存在唯一）2.7.2三彎矩法邊界條件1（固支邊界）邊界條件2（簡支邊界）邊界條件3（周期邊界）例2.14,2.152.7.3m關(guān)系式2.7.4三次樣條插值函數(shù)的性質(zhì)§2.7三次樣條插值引例2.7.1三次樣條插值函數(shù)的概念一背景二、樣條函1引例:y=sinx

在區(qū)間[0，]上的插值逼近

1.二次插值

2.兩點埃爾米特插值3.分段埃爾米特插值x 0 /2

Sinx 0 1 0Cosx 1 0 –1x 0

Sinx 0 0Cosx 1 –1引例:y=sinx在區(qū)間[0，]上的插值逼近2高次插值出現(xiàn)龍格現(xiàn)象L-插值（牛頓插值）Hermite插值分段插值但分段線性插值在節(jié)點處不一定光滑分段Hermite插值但導(dǎo)數(shù)值不容易提?。ㄕ业剑駷榈玫焦饣雀?、應(yīng)用方便的插值函數(shù),我們引入樣條插值函數(shù)?！皹訔l”名詞來源于工程中船體、汽車、飛機等的外形設(shè)計：給出外形曲線上的一組離散點(樣點)，如(xi,yi)，i=0,1,2,…,n,將有彈性的細(xì)長木條或鋼條(樣條)在樣點上固定，使其在其它地方自由彎曲，這樣樣條所表示的曲線，稱為樣條曲線(函數(shù))。一背景2.7.1三次樣條插值函數(shù)的概念高次插值出現(xiàn)龍格現(xiàn)象L-插值（牛頓插值）Hermite插值分3x=-5:5;y=1./(1+x.^2);plot(x,y,x,y,'o')x=-5:5;y=1./(1+x.^2);xi=-5:.05:5;yi=spline(x,y,xi);plot(xi,yi,'b',x,y,'ro')被插值函數(shù):-5≤x≤53/18x=-5:5;x=-5:5;被插值函數(shù):-5≤x≤54x=[0,0.0155,0.1485,0.3493,0.6480,1.0547,2.0];y=[0,0.1242,0.3654,0.4975,0.5472,0.4781,0];n=length(x);t=0:n-1;tt=0:.25:n-1;xx=spline(t,x,tt);yy=spline(t,y,tt);plot(xx,yy,x,y,'o')x=[0,0.0155,0.1485,0.3493,5相同數(shù)據(jù)3次樣條插值與Lagrange插值效果比較CubicSplineInterpolationLagrangeInterpolation三次樣條插值課件6下面介紹應(yīng)用最廣且只有二階連續(xù)導(dǎo)數(shù)的三次樣條函數(shù)．在數(shù)學(xué)上，三次樣條曲線表現(xiàn)為近似于一條分段的三次多項式，它要求在節(jié)點處具有一階和二階連續(xù)導(dǎo)數(shù)。二、樣條函數(shù)的定義定義2.8（三次樣條函數(shù)）在每一個小區(qū)間上是次數(shù)多項式。，即具有連續(xù)的一階，二階導(dǎo)數(shù)。滿足下述條件：如果函數(shù)設(shè)有對[a,b]的剖分的一個3次樣條函數(shù)。為關(guān)于剖分則稱下面介紹應(yīng)用最廣且只有二階連續(xù)導(dǎo)數(shù)的三次樣條函7定義2.8*給定區(qū)間[a,b]上的一個分劃:a=x0<x1<…<xn=b已知

f(xj)=yj(j=0,1,···,n),如果滿足:(1)S(x)在[xj，xj+1]上為三次多項式;(2)S”(x)在區(qū)間[a，b]上連續(xù);(3)S(xj)=yj

(j=0,1,···,n).則稱

S(x)為三次樣條插值函數(shù).定義2.8*給定區(qū)間[a,b]上的一個分劃:滿足:(8注：三次樣條與分段Hermite插值的根本區(qū)別在于S(x)自身光滑，不需要知道f的導(dǎo)數(shù)值（除了在2個端點可能需要）；而Hermite插值依賴于f在所有插值點的導(dǎo)數(shù)值。f(x)H(x)S(x)注：三次樣條與分段Hermite插值的根本區(qū)別在于S(x9三次樣條插值課件10①插值條件:S(xj)=yj

(j=0,1,···,n)n+1個②連續(xù)性條件:S(xj+0)=S(xj-0)

(j=1,···,n-1)S'(xj+0)=S'

(xj-0)

(j=1,···,n-1)S'

'

(xj+0)=S'

'

(xj-0)

(j=1,···,n-1)3(n-1)個共可建立方程(4n-2)個?。》匠虜?shù)少于未知數(shù)個數(shù)？？①插值條件:S(xj)=yj(j=0,1,11

共有個條件,要唯一確定,還必須附加2個條件這兩個條件常在插值區(qū)間[a,b]的邊界點a,b處給出，稱為邊界條件。邊界條件的類型很多，常見的有：③附加2個條件，有多種給法.最常見的給法是:(a)固支邊界(b)簡支邊界特別地,(自然邊界,三次自然樣條);(1)(2)注：一般不取一端是一階導(dǎo)數(shù)而另一端是二階導(dǎo)數(shù)。共有個條件,要唯一確定,還必須附加212第3種邊界條件（周期邊界條件）：注意：上述②給出的個條件是問題本身隱含的，①和③共個獨立條件須提供，故節(jié)點三次樣插值問題只有個自由度.(請與分段三次Hermite插值比較!)第3種邊界條件（周期邊界條件）：注意：上述②給出的13且(1)如果是定義在上函數(shù)且已知函數(shù)表定理2.8（3次樣條插值函數(shù)存在唯一）唯一3次樣條插值函數(shù),且滿足(2)給定邊界條件，則于存在且(1)如果是定義在上函數(shù)且已知函數(shù)表定理2.14例2.13已知f(–1)=1,f(0)=0,f(1)=1.求[–1，1]上的三次自然樣條(滿足自然邊界條件).解設(shè)

則有:S(-1)=–a1+b1–c1+d1=f(-1)=1,S(0)=d1=f(0)=0,S(1)=a2+b2+c2+d2=f(1)=1,

S(0-0)=d1=S(0+0)=d2,S'-(0)=c1=S'+(0)=c2,S''-(0)=b1=S''+(0)=b2

由自然邊界條件:S''(0)=–6a1+2b1=0,S'(1)=6a2+2b2=0解方程組,得a1=-a2=1/2,b1=b2=3/2,c1=c2=d1=d2=0例2.13已知f(–1)=1,f(0)=15問題的解

x=[-1,0,1];y=[1,0,1];f1=inline('0.5*x.^3+1.5*x.^2');f2=inline('-0.5*x.^3+1.5*x.^2');t1=-1:.1:0;t2=0:.1:1;p1=f1(t1);p2=f2(t2);plot(x,y,'o',[t1,t2],[p1,p2],’r’)Holdon,plot([t1,t2],[t1,t2].^2)y=x2問題的解x=[-1,0,1];y=[1,0,1];y=x16三次樣條插值函數(shù)可以有多種表達式，有時用二階導(dǎo)數(shù)值表示時，使用更方便。在力學(xué)上解釋為細(xì)梁在處的彎矩，并且得到的彎矩與相鄰兩個彎矩有關(guān)，故稱用表示的算法為三彎矩算法。2.7.2

構(gòu)造三次樣條插值函數(shù)的三彎矩法

------三次樣條插值函數(shù)的二階導(dǎo)數(shù)表示三次樣條插值函數(shù)可以有多種表達式，有17由兩點拉格朗日插值可表示為參數(shù)對上式積分,得再積分,得由兩點拉格朗日插值可表示為參數(shù)對上式積分,得再積分,得18由條件，確定積分常數(shù)由條件，確定積分常數(shù)19將上式代入(2.48)得到三次樣條插值函數(shù)的表達式由上討論可知,只要確定Mj(j=0,1,…n)這n+1個值,就可定出三次樣條插值函數(shù)S(x)。為了確定Mj(j=0,1,…n),對S(x)求導(dǎo)得將上式代入(2.48)得到三次樣條插值函數(shù)的表達式由上20三次樣條插值課件21三次樣條插值課件22（2.55）

上式兩邊同乘以,即得方程

若記

（2.56）（2.55）上式兩邊同乘以,即得方程若記23所得方程可簡寫成（2.58）

即

（2.57）——三彎矩方程所得方程可簡寫成（2.58）即（2.57）——三彎矩方24這是一個含有n+1個未知數(shù)、n-1個方程的線性方程組.要完全確定Mi(i=0,1,…,n)的值還需要補充兩個條件,這兩個條件通常根據(jù)實際問題的需要，根據(jù)插值區(qū)間a,b的兩個端點處的邊界條件來補充。這是一個含有n+1個未知數(shù)、n-1個方程的線性方程組25由(2.53),得由(2.54),得(1)若已知，則令j=0，令j=n，邊界條件1（固支邊界）-----由(2.53),得由(2.54),得(1)若已知，則令j26對角占優(yōu)的三對角帶狀矩陣對角占優(yōu)的三對角帶狀矩陣27(2)若已知，代入方程(2.58)，只需解n-1個方程邊界條件2（簡支邊界）-----對角占優(yōu)的三對角帶狀矩陣(2)若已知，代入方程(2.58)，只需解n-1個方程邊28(3)對第三類邊界條件：兩邊同除以(j=n)(j=n)(j=0)邊界條件3（周期邊界）-----(3)對第三類邊界條件：兩邊同除以(j=29令得又由,三彎矩方程可寫為令得又由,三彎矩方程可寫為30小結(jié)：在三個邊界條件下的三彎矩方程小結(jié)：在三個邊界條件下的三彎矩方程31說明：（1）方程組(2.59)~(2.61)系數(shù)矩陣都是嚴(yán)格對角占優(yōu)矩陣，因此方程組(2.59)~(2.61)有唯一解（2）Mj在力學(xué)上為細(xì)梁在xj處截面處的彎矩,且彎矩與相鄰的兩個彎矩有關(guān),故方程組(2.59)~(2.61)稱為三彎矩方程。Mj在數(shù)學(xué)上稱為曲率。

實際上,方程組(2.59)~(2.61)的系數(shù)矩陣是一類特殊的矩陣，在后面線性方程組的解法中，將專門介紹這類方程組的解法和性質(zhì)。

說明：（1）方程組(2.59)~(2.61)系32.三次樣條插值課件33三次樣條插值課件34三次樣條插值課件35在本例中，將代入整理后可得：故所求三次樣條插值函數(shù)為：在本例中，將故所求三次樣條插值函數(shù)為：36例2.15

已知的函數(shù)值如下：

x1245

f(x)1342在區(qū)間1,5上求三次樣條插值函數(shù)S(x),使它滿足邊界條件

解:這是在第二種邊界條件下的插值問題，故確定

的方程組形如（2.60）所示，由已知邊界條件,有

則得求解的方程組為

例2.15已知的函數(shù)值如下：在區(qū)間1,5上求三次樣條37根據(jù)給定數(shù)據(jù)和邊界條件算出與根據(jù)給定數(shù)據(jù)和邊界條件算出與38則得方程組解得

又

即得S(x)在各子區(qū)間上的表達式,由式（2.51）知,S(x)在上的表達式為代入式(2.50)將代入上式化簡后得

則得方程組解得又即得S(x)在各子區(qū)間上的表達式代入式39同理S(x)在上的表達式為

S(x)在上的表達式為同理S(x)在上的表達式為S(x)在40故所求的三次樣條插值函數(shù)S(x)在區(qū)間上的表達式為故所求的三次樣條插值函數(shù)S(x)在區(qū)間上的表達式41練習(xí)設(shè)在節(jié)點上，函數(shù)的值為，。試求三次樣條插值函數(shù)，滿足條件

解（1）是固支邊界，先求，再求解，可知練習(xí)設(shè)在節(jié)點42對第一類邊界條件代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有對第一類邊界條件代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)43（２）是簡支條件，不過要注意的不同。由于和已知，故可以化簡得代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有（２）是簡支條件，不過要注意44由此解得。

將代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有由此解得45練習(xí)已知離散點:

(1.1,0.4000),(1.2,0.8000),(1.4,1.6500),(1.5,1.8000),取自然邊界條件

M0=Mn=0,構(gòu)造三次樣條插值函數(shù)，并計算

f(1.25).解

n=3.

∵

h0=x1-x0=0.1,h1=0.2,h2=0.1,因此，分段的三次樣條插值函數(shù)為,由（2.50）由(2.56)計算得練習(xí)已知離散點:(1.1,0.4000),(1.246

上述三次樣條插值的基本思想和特點是：先利用一階導(dǎo)數(shù)在內(nèi)節(jié)點上的連續(xù)性以及邊界條件,列出確定二階導(dǎo)數(shù)的線性方程組（力學(xué)上稱為三彎矩方組），由此解出，再用來表達S(x)。實際上，還可以通過別的途徑來求取三次樣條插值函數(shù)。如：可以先利用二階導(dǎo)數(shù)在內(nèi)節(jié)點上的連續(xù)性及邊界條件，列出確定一階導(dǎo)數(shù)的線性方程組（力學(xué)上稱為三轉(zhuǎn)角方程組），由此解出，再用表達S(x),在某些情況下，這種方法比前者更簡單適用。上述三次樣條插值的基本思想和特點是：實際上，還可以通過別的途472.7.3m關(guān)系式——用一階導(dǎo)數(shù)表示的樣條插值函數(shù)

給定插值點

(xi,yi),設(shè)S'(xi)=mi,

i=0,1,2,…,n,則[xi,xi+1]上的三次Hermite插值為

令

hi=xi+1-xi,∵S(x)∈C2[a,b],對(2.62)求二階導(dǎo)數(shù)2.7.3m關(guān)系式——用一階導(dǎo)數(shù)表示的樣條插值函數(shù)48

令

xi+=xi+0,在[xi,xi+1]上得到xi點的右導(dǎo)數(shù)，

同理，在[xi-1,xi]上構(gòu)造三次樣條插值

S(x),在[xi-1,xi]上得點xi的左導(dǎo)數(shù)，令xi+=xi+0,在[xi,xi+1]49三種邊界條件：

三種邊界條件：50由此可解得m1,m2,…,mn-1，從而得S(x)的表達式.（2.66）①對于邊界條件(1),兩個方程則m1,m2,…,mn-1滿足方程組

由此可解得m1,m2,…,mn-1，從而得S(x)的表51①對于邊界條件(2),可導(dǎo)出兩個方程:(2.67)①對于邊界條件(2),可導(dǎo)出兩個方程:(2.67)52若令則（2.65）和（2.67）可合并成矩陣形式（2.68）可解出從而得S(x)的表達式.若令則（2.65）和（2.67）可合并成矩陣形式（2.68）53由(2.65)和(2.6)可解出，方程組的矩陣形式為③對于邊界條件（3），可得（2.69）其中（2.70）由(2.65)和(2.6)可解出54在實際應(yīng)用中，如果不需要規(guī)定內(nèi)節(jié)點處的一階導(dǎo)數(shù)值，那么使用三次樣條插值函數(shù)會得到很好的效果。三次樣條插值函數(shù)不僅在內(nèi)節(jié)點處的二階導(dǎo)數(shù)是連續(xù)的，而且逼近具有很好的收斂性，也是數(shù)值穩(wěn)定的。下面給出三次樣條插值函數(shù)的一些重要性質(zhì)。2.7.4三次樣條插值函數(shù)的性質(zhì)在實際應(yīng)用中，如果不需要規(guī)定內(nèi)節(jié)點處的一階導(dǎo)數(shù)值，那么使55值函數(shù),則有估計式定理2.9設(shè)函數(shù)記則對任意滿足邊界條件（2.44）或（2.45）的三次樣條插（2.69）其中由于誤差估計與收斂性定理的證明比較復(fù)雜，下面只給出誤差估計的結(jié)論。值函數(shù),則有估計式定理2.9設(shè)函數(shù)記56

誤差估計式（2.69）除可以用于誤差估計外,它進一步表明，當(dāng)時，在插值區(qū)間上，對于滿足邊界條件（2.44）或（2..45）的插值函數(shù)，不僅一致收斂于，而且一致收斂于，一致收斂于。用三次樣條繪制的曲線不僅有很好的光滑度，而且當(dāng)節(jié)點逐漸加密時，其函數(shù)值在整體上能很好地逼近被插函數(shù)，相應(yīng)的導(dǎo)數(shù)值也收斂于被插函數(shù)的導(dǎo)數(shù)，不會發(fā)生龍格現(xiàn)象。因此三次樣條在計算機輔助設(shè)計中有廣泛的應(yīng)用。誤差估計式（2.69）除可以用于誤差估計外,它進一步57

插值名稱

插值條件

構(gòu)造方法優(yōu)點簡便計算方法

缺點L-插值先求基函數(shù)再求插值函數(shù)對稱秦九韶算法(1)計算量大(2)數(shù)值不穩(wěn)定、不收斂(Runge現(xiàn)象)本章內(nèi)容：H-插值一、同上二、利用牛頓插值

收斂同上高次有Runge現(xiàn)象分段插值基函數(shù)法收斂

同上

光滑性差三次樣條插值基函數(shù)法收斂同上光滑性有改進改進方法：列維爾算法、埃特金算法、牛頓法注：分段插值中乍看上去沒有構(gòu)造插值基函數(shù)，但實際上線性插值用了線性L-插值基函數(shù)，三次樣條插值的基函數(shù)可為：插值插值構(gòu)造優(yōu)點簡便計算方法58引例2.7.1三次樣條插值函數(shù)的概念一背景二、樣條函數(shù)的定義

例2.13

定理2.8（3次樣條插值函數(shù)存在唯一）2.7.2三彎矩法邊界條件1（固支邊界）邊界條件2（簡支邊界）邊界條件3（周期邊界）例2.14,2.152.7.3m關(guān)系式2.7.4三次樣條插值函數(shù)的性質(zhì)§2.7三次樣條插值引例2.7.1三次樣條插值函數(shù)的概念一背景二、樣條函59引例:y=sinx

在區(qū)間[0，]上的插值逼近

1.二次插值

2.兩點埃爾米特插值3.分段埃爾米特插值x 0 /2

Sinx 0 1 0Cosx 1 0 –1x 0

Sinx 0 0Cosx 1 –1引例:y=sinx在區(qū)間[0，]上的插值逼近60高次插值出現(xiàn)龍格現(xiàn)象L-插值（牛頓插值）Hermite插值分段插值但分段線性插值在節(jié)點處不一定光滑分段Hermite插值但導(dǎo)數(shù)值不容易提?。ㄕ业剑駷榈玫焦饣雀摺?yīng)用方便的插值函數(shù),我們引入樣條插值函數(shù)?！皹訔l”名詞來源于工程中船體、汽車、飛機等的外形設(shè)計：給出外形曲線上的一組離散點(樣點)，如(xi,yi)，i=0,1,2,…,n,將有彈性的細(xì)長木條或鋼條(樣條)在樣點上固定，使其在其它地方自由彎曲，這樣樣條所表示的曲線，稱為樣條曲線(函數(shù))。一背景2.7.1三次樣條插值函數(shù)的概念高次插值出現(xiàn)龍格現(xiàn)象L-插值（牛頓插值）Hermite插值分61x=-5:5;y=1./(1+x.^2);plot(x,y,x,y,'o')x=-5:5;y=1./(1+x.^2);xi=-5:.05:5;yi=spline(x,y,xi);plot(xi,yi,'b',x,y,'ro')被插值函數(shù):-5≤x≤53/18x=-5:5;x=-5:5;被插值函數(shù):-5≤x≤562x=[0,0.0155,0.1485,0.3493,0.6480,1.0547,2.0];y=[0,0.1242,0.3654,0.4975,0.5472,0.4781,0];n=length(x);t=0:n-1;tt=0:.25:n-1;xx=spline(t,x,tt);yy=spline(t,y,tt);plot(xx,yy,x,y,'o')x=[0,0.0155,0.1485,0.3493,63相同數(shù)據(jù)3次樣條插值與Lagrange插值效果比較CubicSplineInterpolationLagrangeInterpolation三次樣條插值課件64下面介紹應(yīng)用最廣且只有二階連續(xù)導(dǎo)數(shù)的三次樣條函數(shù)．在數(shù)學(xué)上，三次樣條曲線表現(xiàn)為近似于一條分段的三次多項式，它要求在節(jié)點處具有一階和二階連續(xù)導(dǎo)數(shù)。二、樣條函數(shù)的定義定義2.8（三次樣條函數(shù)）在每一個小區(qū)間上是次數(shù)多項式。，即具有連續(xù)的一階，二階導(dǎo)數(shù)。滿足下述條件：如果函數(shù)設(shè)有對[a,b]的剖分的一個3次樣條函數(shù)。為關(guān)于剖分則稱下面介紹應(yīng)用最廣且只有二階連續(xù)導(dǎo)數(shù)的三次樣條函65定義2.8*給定區(qū)間[a,b]上的一個分劃:a=x0<x1<…<xn=b已知

f(xj)=yj(j=0,1,···,n),如果滿足:(1)S(x)在[xj，xj+1]上為三次多項式;(2)S”(x)在區(qū)間[a，b]上連續(xù);(3)S(xj)=yj

(j=0,1,···,n).則稱

S(x)為三次樣條插值函數(shù).定義2.8*給定區(qū)間[a,b]上的一個分劃:滿足:(66注：三次樣條與分段Hermite插值的根本區(qū)別在于S(x)自身光滑，不需要知道f的導(dǎo)數(shù)值（除了在2個端點可能需要）；而Hermite插值依賴于f在所有插值點的導(dǎo)數(shù)值。f(x)H(x)S(x)注：三次樣條與分段Hermite插值的根本區(qū)別在于S(x67三次樣條插值課件68①插值條件:S(xj)=yj

(j=0,1,···,n)n+1個②連續(xù)性條件:S(xj+0)=S(xj-0)

(j=1,···,n-1)S'(xj+0)=S'

(xj-0)

(j=1,···,n-1)S'

'

(xj+0)=S'

'

(xj-0)

(j=1,···,n-1)3(n-1)個共可建立方程(4n-2)個！！方程數(shù)少于未知數(shù)個數(shù)？？①插值條件:S(xj)=yj(j=0,1,69

共有個條件,要唯一確定,還必須附加2個條件這兩個條件常在插值區(qū)間[a,b]的邊界點a,b處給出，稱為邊界條件。邊界條件的類型很多，常見的有：③附加2個條件，有多種給法.最常見的給法是:(a)固支邊界(b)簡支邊界特別地,(自然邊界,三次自然樣條);(1)(2)注：一般不取一端是一階導(dǎo)數(shù)而另一端是二階導(dǎo)數(shù)。共有個條件,要唯一確定,還必須附加270第3種邊界條件（周期邊界條件）：注意：上述②給出的個條件是問題本身隱含的，①和③共個獨立條件須提供，故節(jié)點三次樣插值問題只有個自由度.(請與分段三次Hermite插值比較!)第3種邊界條件（周期邊界條件）：注意：上述②給出的71且(1)如果是定義在上函數(shù)且已知函數(shù)表定理2.8（3次樣條插值函數(shù)存在唯一）唯一3次樣條插值函數(shù),且滿足(2)給定邊界條件，則于存在且(1)如果是定義在上函數(shù)且已知函數(shù)表定理2.72例2.13已知f(–1)=1,f(0)=0,f(1)=1.求[–1，1]上的三次自然樣條(滿足自然邊界條件).解設(shè)

則有:S(-1)=–a1+b1–c1+d1=f(-1)=1,S(0)=d1=f(0)=0,S(1)=a2+b2+c2+d2=f(1)=1,

S(0-0)=d1=S(0+0)=d2,S'-(0)=c1=S'+(0)=c2,S''-(0)=b1=S''+(0)=b2

由自然邊界條件:S''(0)=–6a1+2b1=0,S'(1)=6a2+2b2=0解方程組,得a1=-a2=1/2,b1=b2=3/2,c1=c2=d1=d2=0例2.13已知f(–1)=1,f(0)=73問題的解

x=[-1,0,1];y=[1,0,1];f1=inline('0.5*x.^3+1.5*x.^2');f2=inline('-0.5*x.^3+1.5*x.^2');t1=-1:.1:0;t2=0:.1:1;p1=f1(t1);p2=f2(t2);plot(x,y,'o',[t1,t2],[p1,p2],’r’)Holdon,plot([t1,t2],[t1,t2].^2)y=x2問題的解x=[-1,0,1];y=[1,0,1];y=x74三次樣條插值函數(shù)可以有多種表達式，有時用二階導(dǎo)數(shù)值表示時，使用更方便。在力學(xué)上解釋為細(xì)梁在處的彎矩，并且得到的彎矩與相鄰兩個彎矩有關(guān)，故稱用表示的算法為三彎矩算法。2.7.2

構(gòu)造三次樣條插值函數(shù)的三彎矩法

------三次樣條插值函數(shù)的二階導(dǎo)數(shù)表示三次樣條插值函數(shù)可以有多種表達式，有75由兩點拉格朗日插值可表示為參數(shù)對上式積分,得再積分,得由兩點拉格朗日插值可表示為參數(shù)對上式積分,得再積分,得76由條件，確定積分常數(shù)由條件，確定積分常數(shù)77將上式代入(2.48)得到三次樣條插值函數(shù)的表達式由上討論可知,只要確定Mj(j=0,1,…n)這n+1個值,就可定出三次樣條插值函數(shù)S(x)。為了確定Mj(j=0,1,…n),對S(x)求導(dǎo)得將上式代入(2.48)得到三次樣條插值函數(shù)的表達式由上78三次樣條插值課件79三次樣條插值課件80（2.55）

上式兩邊同乘以,即得方程

若記

（2.56）（2.55）上式兩邊同乘以,即得方程若記81所得方程可簡寫成（2.58）

即

（2.57）——三彎矩方程所得方程可簡寫成（2.58）即（2.57）——三彎矩方82這是一個含有n+1個未知數(shù)、n-1個方程的線性方程組.要完全確定Mi(i=0,1,…,n)的值還需要補充兩個條件,這兩個條件通常根據(jù)實際問題的需要，根據(jù)插值區(qū)間a,b的兩個端點處的邊界條件來補充。這是一個含有n+1個未知數(shù)、n-1個方程的線性方程組83由(2.53),得由(2.54),得(1)若已知，則令j=0，令j=n，邊界條件1（固支邊界）-----由(2.53),得由(2.54),得(1)若已知，則令j84對角占優(yōu)的三對角帶狀矩陣對角占優(yōu)的三對角帶狀矩陣85(2)若已知，代入方程(2.58)，只需解n-1個方程邊界條件2（簡支邊界）-----對角占優(yōu)的三對角帶狀矩陣(2)若已知，代入方程(2.58)，只需解n-1個方程邊86(3)對第三類邊界條件：兩邊同除以(j=n)(j=n)(j=0)邊界條件3（周期邊界）-----(3)對第三類邊界條件：兩邊同除以(j=87令得又由,三彎矩方程可寫為令得又由,三彎矩方程可寫為88小結(jié)：在三個邊界條件下的三彎矩方程小結(jié)：在三個邊界條件下的三彎矩方程89說明：（1）方程組(2.59)~(2.61)系數(shù)矩陣都是嚴(yán)格對角占優(yōu)矩陣，因此方程組(2.59)~(2.61)有唯一解（2）Mj在力學(xué)上為細(xì)梁在xj處截面處的彎矩,且彎矩與相鄰的兩個彎矩有關(guān),故方程組(2.59)~(2.61)稱為三彎矩方程。Mj在數(shù)學(xué)上稱為曲率。

實際上,方程組(2.59)~(2.61)的系數(shù)矩陣是一類特殊的矩陣，在后面線性方程組的解法中，將專門介紹這類方程組的解法和性質(zhì)。

說明：（1）方程組(2.59)~(2.61)系90.三次樣條插值課件91三次樣條插值課件92三次樣條插值課件93在本例中，將代入整理后可得：故所求三次樣條插值函數(shù)為：在本例中，將故所求三次樣條插值函數(shù)為：94例2.15

已知的函數(shù)值如下：

x1245

f(x)1342在區(qū)間1,5上求三次樣條插值函數(shù)S(x),使它滿足邊界條件

解:這是在第二種邊界條件下的插值問題，故確定

的方程組形如（2.60）所示，由已知邊界條件,有

則得求解的方程組為

例2.15已知的函數(shù)值如下：在區(qū)間1,5上求三次樣條95根據(jù)給定數(shù)據(jù)和邊界條件算出與根據(jù)給定數(shù)據(jù)和邊界條件算出與96則得方程組解得

又

即得S(x)在各子區(qū)間上的表達式,由式（2.51）知,S(x)在上的表達式為代入式(2.50)將代入上式化簡后得

則得方程組解得又即得S(x)在各子區(qū)間上的表達式代入式97同理S(x)在上的表達式為

S(x)在上的表達式為同理S(x)在上的表達式為S(x)在98故所求的三次樣條插值函數(shù)S(x)在區(qū)間上的表達式為故所求的三次樣條插值函數(shù)S(x)在區(qū)間上的表達式99練習(xí)設(shè)在節(jié)點上，函數(shù)的值為，。試求三次樣條插值函數(shù)，滿足條件

解（1）是固支邊界，先求，再求解，可知練習(xí)設(shè)在節(jié)點100對第一類邊界條件代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有對第一類邊界條件代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)101（２）是簡支條件，不過要注意的不同。由于和已知，故可以化簡得代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有（２）是簡支條件，不過要注意102由此解得。

將代入三次樣條插值函數(shù)的表達式（2.50），經(jīng)化簡有由此解得103練習(xí)已知離散點:

(1.1,0.4000),(1.2,0.8000),(1.4,1.6500),(1.5,1.8000),取自然邊界條件

M0=Mn=0,構(gòu)造三次樣條插值函數(shù)，并計算

f(1.25).解

n=3.

∵

h0=x1-x0=0.1,h1=0.2,h2=0.1,因此，分段的三次樣條插值函數(shù)為,由（2.50）由(2.56)計算得練習(xí)已知離散點:(1.1,0.4000),(1.2104

上述三次樣條插值的基本思想和特點是：先利用一階導(dǎo)數(shù)在內(nèi)節(jié)點上的連續(xù)性以及邊界條件,列出確定二階導(dǎo)數(shù)的線性方程組（力學(xué)上稱為三彎矩方組），由此解出，再用來表達S(x)。實際上，還可以通過別的途徑來求取三次樣條插值函數(shù)。如：可以先利用二階導(dǎo)數(shù)在內(nèi)節(jié)點上的連續(xù)性及邊界條件，列出確定一階導(dǎo)數(shù)的線性方程組（力學(xué)上稱為三轉(zhuǎn)角方程組），由此解出，再用表達S(x),在某些情況下，這種方法比前者更簡單適用。上述三次樣條插值的基本思想和特點是：實際上，還可以通過別的途1052.7.3m關(guān)系式——用一階導(dǎo)數(shù)表示的樣條插值函數(shù)

給定插值點

(xi,yi),設(shè)S'(xi)=mi,

i=0,1,2,…,n,則[xi,xi+1]上的三次Hermite插值為

令

hi=x