看zzq大佬的博客,看到了这个看似很深奥的东西,实际很简单(反正比FFT简单,我是一个要被FFT整疯了的孩子)
拉格朗日插值法
是什么
可以找到一个多项式,其恰好在各个观测点取到观测到的值。这样的多项式称为拉格朗日(插值)多项式
数学上来说,拉格朗日插值法可以给出一个恰好穿过二维平面上若干个点的多项式函数。
定义
概念
一般地,若已知y=f(x)在互不相同n+1个点x0,x1,...xn处的函数值y0,y1,...yn(即该函数(x0,y0),(x1,y1),...(xn,yn)这n+1个点)
则可考虑构造一个过这n+1个点的、次数不超过n的多项式y=Pn(x),使其满足:
Pn(xk)=yk,k=0,1,...,n(*)
//就是说,对于很多很多个点,我们可以找到一个满足所有点的函数解析式。
要估计任一点v,v!=xi,i=0,1,...,n,则可以用Pn(v)的值作为准确值f(v)的近似值,此方法叫“插值法”。
称式(*)为插值条件(准则),含xi(i=0,1,...,n)的最小区间[a,b],其中a=min{x0,x1,...,xn},b=max{x0,x1,...,xn}.
定义
满足插值条件的、次数不超过n的多项式是存在而且唯一的。
一般形式运用方法
在平面上有
共n个点,现作一条函数
使其图像经过这n个点。
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekF1WW1SemRHRjBhV011WTI5dEx5MDBiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNekl2YzJsbmJqMWhOell6WVRnMlpEWTFaRGxtTW1Rek1qUXhNVEl5WldSaE9HVmpPRFF3Tnk4eVpUSmxZamt6T0RsaU5UQTBabU15TURsbU5USTNPRGxsTW1Sa1pUY3hNVGt4WldZMlpEa3dMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
使其图像经过这n个点。
作法:设集合
是关于点
,都有
使得
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjek11WW1SemRHRjBhV011WTI5dEx6ZFFiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNVGt2YzJsbmJqMDBPVFF6T0RobVptSTVNV000TnpBeFpESmlObUkyWldZeU56ZG1ZamN4TVM4d1pHWTBNekZoWkdOaVpXWTNOakE1TW1Nd05tRm1PRFF5TkdSa1lUTmpZemRqWkRrNVpUbG1MbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekl1WW1SemRHRjBhV011WTI5dEx6bG1iek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNek12YzJsbmJqMHdNbU0yTWpabFptRmtaV1pqWlRGaVpXVXlZbU5sWXpsaFpqVXhNMlF3TUM5a09ETXpZemc1TldReE5ETmhaRFJpTkRNellUSXpaRGs0T0RBeU5XRmhabUUwTUdZd05qazJMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
的角标的集合,
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjek11WW1SemRHRjBhV011WTI5dEx5MVFiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNVE0zTDNOcFoyNDlaakE0WkdJME5tVTBPV0U1T0RJeU5tSmpZekV5WmpJMFltUTRNbUk1TjJFdk5qTmtPV1l5WkRNMU56SmpNVEZrWm1JeVlUVTRPVE00TmpreU56WXlaREJtTnpBell6SmlNQzVxY0djPS5qcGc%3D.jpg?w=700&webp=1)
,作n个多项式
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekl1WW1SemRHRjBhV011WTI5dEx5MW1iek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVPRE12YzJsbmJqMDBOelprTXpZd01XTmpOV014TURNNE1qQTNaV016WXpGaU1qRXhZemt4WVM4NU5HTmhaREZqT0dFM09EWmpPVEUzWkRFd016WTBNakpqTXpOa056Qmpaak5pWXpjMU4yRmpMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
。对于任意
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjek11WW1SemRHRjBhV011WTI5dEx6ZFFiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVORFF2YzJsbmJqMDROalkyWW1KaU1tRXpNREUwWXpBNE1XUXpZakk1WVRFd1lUZGlPR014WXk4M1pHUTVPR1F4TURBeFpUa3pPVEF4TkdVM1pHSmlNalUzTVdWak5UUmxOek0yWkRFNU5tRTVMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjek11WW1SemRHRjBhV011WTI5dEx6ZFFiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNVGd3TDNOcFoyNDlObVF4T1dZM016UmlOREE1Tm1JMk16ZzFNVGsxWVRVNE0yTXpNamczTXpNdlkyRmxaamMyTURrMFlqTTJZV05oWmpRNU9ESmhObVUzTnpaa09UaGtNVEF3TVdVNU9XTTBaUzVxY0djPS5qcGc%3D.jpg?w=700&webp=1)
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekV1WW1SemRHRjBhV011WTI5dEx5MTJiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNell2YzJsbmJqMW1NMlJoTlRGbFlqYzJaRGs0WkRFd056SmtOREJoTXpjeU1UTm1aREEyTUM5aFl6WmxaR1JqTkRVeFpHRTRNV05pWlRSa1l6UTBaR0UxT0RZMlpEQXhOakE1TWpRek1XVmlMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
是n-1次多项式,且满足
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjek11WW1SemRHRjBhV011WTI5dEx5MVFiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVNVE0zTDNOcFoyNDlNRGRqT0RjNVpqZ3lOVEpsWWprek9HVTRObVEzWldZeFpUSTJNamcxWm1VdlkyUmlaalpqT0RFNE1EQmhNVGxrT0dWbU1URm1aVEZqTXpsbVlUZ3lPR0poTmpGbE5EWTVZeTVxY0djPS5qcGc%3D.jpg?w=700&webp=1)
并且
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekF1WW1SemRHRjBhV011WTI5dEx6azBiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVOamt2YzJsbmJqMDFZbVUyTkRJMU9XTXdPREEyTlRNNE4yWmxZV0UzTVdFNU5tUmtZekkxTnk4NVpUTmtaamhrWTJReE1EQmlZV0V4TURWbE5EZ3dZV00wWkRFd1lqa3hNbU01Wm1NeVpXUTVMbXB3Wnc9PS5qcGc%3D.jpg?w=700&webp=1)
。最后可得 。形如上式的插值多项式
称为拉格朗日(Lagrange)插值多项式。
例如:当n=4时,上面的公式可简化为:
![集训DAYn——拉格朗日插值法 集训DAYn——拉格朗日插值法](https://image.shishitao.com:8440/aHR0cHM6Ly9iYnNtYXguaWthZmFuLmNvbS9zdGF0aWMvTDNCeWIzaDVMMmgwZEhCekwyZHpjekV1WW1SemRHRjBhV011WTI5dEx5MTJiek5rVTJGblgzaEpOR3RvUjJ0d2IxZExNVWhHTm1ob2VTOWlZV2xyWlM5ekpUTkVPRGcwTDNOcFoyNDlabVF3WVRJeU16UTVabU5oWkRGak9HUTBZbUptTVRKbU5HSXpaVFkzWXpRdllqZG1aRFV5Tmpaa01ERTJNRGt5TkRZMk5XTTVNREF3WkdVd056TTFabUZsTm1Oa016UmhZaTVxY0djPS5qcGc%3D.jpg?w=700&webp=1)
这是一个过4个点的唯一的三次多项式。
其实很简单?!