[Python] Hermite 插值

　　　　　　　　

 

 

# -*- coding: utf-8 -*-
#Program 0.5 Hermite Interpolation

import matplotlib.pyplot as plt
import numpy as np

#计算基函数的导数值
def dl(i, xi):
	result = 0.0
	for j in range(0,len(xi)):
		if j!=i:
			result += 1/(xi[i]-xi[j])
	result *= 2
	return result

#计算基函数值
def l(i, xi, x):
	deno = 1.0
	nu = 1.0

	for j in range(0, len(xi)):
		if j!= i:
			deno *= (xi[i]-xi[j])
			nu *