polynomial interpolate

最新推荐文章于 2024-04-17 11:30:00 发布

一轮子

最新推荐文章于 2024-04-17 11:30:00 发布

阅读量803

点赞数

分类专栏：数据科学文章标签： analyse python 差值定理

本文链接：https://blog.csdn.net/ddyyll988/article/details/54913286

版权

数据科学专栏收录该内容

8 篇文章 0 订阅

订阅专栏

polynomial interpolate

Wiki Polynomial interpolation\

scipy

solve this system of linear equations:

p (x) = a n x n + a n - 1 x n - 1 + \dots + a 2 x 2 + a 1 x + a 0

$p(x)=a_nx^n+a_{n-1}x^{n-1}+\dots+a_2x^2+a_1x+a_0$

We have to solve this system for $a_k$ to construct the interpolant $p(x)$ . The matrix on the left is commonly referred to as a Vandermonde matrix.

Lagrange polynomial

WIki

$Lagrange polynomials$

example

x	f(x)
1	1
2	8
3	27

∴ L (1.5) = 3

$\therefore L(1.5)=3$

code

from scipy.interpolate import lagrange
L = lagrange([1,2,3],[1,8,27])
print(L)
#poly1d([  6., -11.,   6.])
L(1.5)
#3.0

#check
import numpy as np
vandermonate = np.vander([1,2,3])
A = L.coeffs
y = np.dot(vandermonate,A)
print(y)

Vandermonde:

numpy.vander

polyvander

deal with null numbers using interpolate

pandas

filled = df.interpolate(method='spline', order=4)

#view
df.plot()
filled.plot()

method: {‘linear’, ‘time’, ‘index’, ‘values’, ‘nearest’, ‘zero’,‘slinear’, ‘quadratic’, ‘cubic’, ‘barycentric’, ‘krogh’, ‘polynomial’, ‘spline’ ‘piecewise_polynomial’, ‘pchip’}

see: scipy