【數(shù)學(xué)建模算法】(25)插值與擬合:分段線性插值

1.插值多項(xiàng)式的震蕩

用 Lagrange 插值多項(xiàng)式L_{n}(x)近似f(x)(a \leq x \leq b),雖然隨著節(jié)點(diǎn)個(gè)數(shù)的增加,L_{n}(x)的次數(shù)n變大,多數(shù)情況下誤差\left|R_{n}(x)\right|會(huì)變小。但是n增大時(shí),L_{n}(x)的光滑性變壞,有時(shí)會(huì)出現(xiàn)很大的振蕩。理論上,當(dāng)n \rightarrow \infty,在[a, b]內(nèi)并不能保證L_{n}(x)處處收斂于f(x)。Runge給出了一個(gè)有名的例子:

f(x)=\frac{1}{1+x^{2}}, \quad x \in[-5,5]

對(duì)于較大的|x|,隨著n的增大,L_{n}(x)震蕩越來(lái)越大,事實(shí)上可以證明,僅當(dāng)x \leq 3.63時(shí),才有\lim _{n \rightarrow \infty} L_{n}(x)=f(x),而在此區(qū)間外,L_{n}(x)是發(fā)散的。

高次插值多項(xiàng)式的這些缺陷,促使人們轉(zhuǎn)而尋求簡(jiǎn)單的低次多項(xiàng)式插值。

2.分段線性插值

簡(jiǎn)單地說(shuō),將每?jī)蓚€(gè)相鄰的節(jié)點(diǎn)用直線連起來(lái),如此形成的一條折線就是分段線性插值函數(shù),記做I_{n}(x),它滿足I_{n}\left(x_{i}\right)=y_{i},且I_{n}(x)在每個(gè)小區(qū)間\left[x_{i}, x_{i+1}\right]是線性函數(shù)(i=0,1, \cdots, n)。

I_{n}(x)可表示為
I_{n}(x)=\sum_{i=0}^{n} y_{i} l_{i}(x)
l_{i}(x)=\left\{\begin{array}{l}{\frac{x-x_{i-1}}{x_{i}-x_{i-1}}},x \in\left[x_{i-1}, x_{i}\right](i=0時(shí)舍去) \\ {\frac{x-x_{i+1}}{x_{i}-x_{i+1}}} ,x \in\left[x_{i}, x_{i+1}\right](i=n時(shí)舍去)\\ {0},其他\end{array}\right.

I_{n}(x)有良好的收斂性,即對(duì)于x \in[a, b]有,
\lim _{n \rightarrow \infty} I_{n}(x)=f(x)
I_{n}(x)計(jì)算x點(diǎn)的插值時(shí),只用到x左右的兩個(gè)節(jié)點(diǎn),計(jì)算量與節(jié)點(diǎn)個(gè)數(shù)n 無(wú)關(guān)。但n越大,分段越多,插值誤差越小。實(shí)際上用函數(shù)表作插值計(jì)算時(shí),分段線性插值就足夠了,如數(shù)學(xué)、物理中用的特殊函數(shù)表,數(shù)理統(tǒng)計(jì)中用的概率分布表等。

3.用Matlab實(shí)現(xiàn)分段線性插值

用 Matlab 實(shí)現(xiàn)分段線性插值不需要編制函數(shù)程序,Matlab 中有現(xiàn)成的一維插值函數(shù) interp1。

y=interp1(x0,y0,x,'method')
method 指定插值的方法,默認(rèn)為線性插值。其值可為:

method值 功能
'nearest' 最近項(xiàng)插值
'linear' 線性插值
'spline' 逐段3次樣條插值
'cubic' 保凹凸性3次插值

所有的插值方法要求 x0 是單調(diào)的。
當(dāng) x0 為等距時(shí)可以用快速插值法,使用快速插值法的格為'*nearest'、'*linear'、
'*spline'、'*cubic'。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容