1. 程式人生 > >雙線性二次插值原理解析

雙線性二次插值原理解析

對某個多項式函式,已知有給定的k + 1個取值點:

(x_{0},y_{0}),\ldots ,(x_{k},y_{k})

其中x_{j}對應著自變數的位置,而y_{j}對應著函式在這個位置的取值。

假設任意兩個不同的xj都互不相同,那麼應用拉格朗日插值公式所得到的拉格朗日插值多項式為:

L(x):=\sum _{{j=0}}^{{k}}y_{j}\ell _{j}(x)

其中每個\ell _{j}(x)拉格朗日基本多項式(或稱插值基函式),其表示式為:

\ell _{j}(x):=\prod _{{i=0,\,i\neq j}}^{{k}}{\frac  {x-x_{i}}{x_{j}-x_{i}}}={\frac  {(x-x_{0})}{(x_{j}-x_{0})}}\cdots {\frac  {(x-x_{{j-1}})}{(x_{j}-x_{{j-1}})}}{\frac  {(x-x_{{j+1}})}{(x_{j}-x_{{j+1}})}}\cdots {\frac  {(x-x_{{k}})}{(x_{j}-x_{{k}})}}.[3]

拉格朗日基本多項式\ell _{j}(x)的特點是在x_{j}上取值為1,在其它的點x_{i},\,i\neq j上取值為0

範例

假設有某個二次多項式函式f,已知它在三個點上的取值為:

  • f(4)=10
  • f(5)=5.25
  • f(6)=1

要求f(18)的值。

首先寫出每個拉格朗日基本多項式:

\ell _{0}(x)={\frac  {(x-5)(x-6)}{(4-5)(4-6)}}
\ell _{1}(x)={\frac  {(x-4)(x-6)}{(5-4)(5-6)}}
\ell _{2}(x)={\frac  {(x-4)(x-5)}{(6-4)(6-5)}}

然後應用拉格朗日插值法,就可以得到p的表示式(p為函式f的插值函式):

p(x)=f(4)\ell _{0}(x)+f(5)\ell _{1}(x)+f(6)\ell _{2}(x)
.\,\,\,\,\,\,\,\,\,\,=10\cdot {\frac  {(x-5)(x-6)}{(4-5)(4-6)}}+5.25\cdot {\frac  {(x-4)(x-6)}{(5-4)(5-6)}}+1\cdot {\frac  {(x-4)(x-5)}{(6-4)(6-5)}}
.\,\,\,\,\,\,\,\,\,\,={\frac  {1}{4}}(x^{2}-28x+136)

此時代入數值\ 18就可以求出所需之值:\ f(18)=p(18)=-11

證明

存在性

對於給定的k+1個點:(x_{0},y_{0}),\ldots ,(x_{k},y_{k}),拉格朗日插值法的思路是找到一個在一點x_{j}

取值為1,而在其他點取值都是0的多項式\ell _{j}(x)。這樣,多項式y_{j}\ell _{j}(x)在點x_{j}取值為y_{j},而在其他點取值都是0。而多項式L(x):=\sum _{{j=0}}^{{k}}y_{j}\ell _{j}(x)就可以滿足

L(x_{j})=\sum _{{i=0}}^{{k}}y_{i}\ell _{i}(x_{j})=0+0+\cdots +y_{j}+\cdots +0=y_{j}

在其它點取值為0的多項式容易找到,例如:

(x-x_{0})\cdots (x-x_{{j-1}})(x-x_{{j+1}})\cdots (x-x_{{k}})

它在點x_{j}取值為:(x_{j}-x_{0})\cdots (x_{j}-x_{{j-1}})(x_{j}-x_{{j+1}})\cdots (x_{j}-x_{{k}})。由於已經假定x_{i}兩兩互不相同,因此上面的取值不等於0。於是,將多項式除以這個取值,就得到一個滿足“在x_{j}取值為1,而在其他點取值都是0的多項式”:

\ell _{j}(x):=\prod _{{i=0,\,i\neq j}}^{{k}}{\frac  {x-x_{i}}{x_{j}-x_{i}}}={\frac  {(x-x_{0})}{(x_{j}-x_{0})}}\cdots {\frac  {(x-x_{{j-1}})}{(x_{j}-x_{{j-1}})}}{\frac  {(x-x_{{j+1}})}{(x_{j}-x_{{j+1}})}}\cdots {\frac  {(x-x_{{k}})}{(x_{j}-x_{{k}})}}

這就是拉格朗日基本多項式。


我們的方法是這樣的,根據水平方向上的雙線性二次插值,由f(I,j)和f(i+1,j)求取f(x,j),由

f(I,j+1)和f(i+1,j+1)求取f(x,j+1),然後再根據這兩點的二次插值求取f(x,y)。

  根據前面的例題,我們可以很容易的求取各點插值如下:

                        f(x,j)=(i+1-x)f(I,j)+(x-i)f(i+1,j)               公式1-(4)

                      f(x,j+1)=(i+1-x)f(I,j+1)+(x-i)f(i+1,j+1)           公式1-(5)

                       f(x,y)=(i+1-y)f(x,j)+(y-j)f(x,j+1)               公式1-(6)

  以上三式綜合可以得到:

  f(x,y)=(j+1-y)(i+1-x)f(I,j)+(j+1-y)(x-i)f(i+1,j)+(y-j)(i+1-x)f(I,j+1)+(y-j)(x-i)f(i+1,j+1)     公式1-(7)

  我們令x=i+p,y=j+q得:

  f(i+p,j+q)=(1-q)(1-p)f(I,j)+p(1-q)f(i+1,j)+q(1-p)f(I,j+1)+pqf(i+1,j+1)                公式1-(8)

  上式即為數字影象處理中的雙線性二次插值公式。

參考部落格:https://www.cnblogs.com/ECJTUACM-873284962/p/6833391.html

                   http://blog.csdn.net/trent1985/article/details/45150677