1.拉格朗日差值法
python中的scipy
库提供了拉格朗日差值法的函数,直接调用
案例
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy.interpolate import lagrange
cat_sale = pd.read_excel('data/catering_sale.xls')
cat_sale.drop('日期', axis=1, inplace=True)
# 过滤异常值,并置为空值
cat_sale['销量'][(cat_sale['销量'] < 400) | (cat_sale['销量'] > 5000)] = np.NAN
# 自定义列向量插值函数
def ployinterp_columns(s, n, k=4):
if n < k:
y = s[list(range(0, n)) + list(range(n + 1, n + k + 1))]
elif n > len(s) - k - 1:
y = s[list(range(n - k, n)) + list(range(n + 1, len(s)))]
else:
y = s[list