拉格朗日插值——基于scipy实现
来源:互联网 发布:权限管理的软件 编辑:程序博客网 时间:2024/05/30 19:35
拉格朗日插值数学原理
在节点上给出节点基函数,然后做基函数的线性组合,组合系数为节点函数值,这种插值多项式称为拉格朗日插值公式
已最简单的两点插值为例,已知函数y = f (x)在给定互异点x0, x1上的值为y0= f (x0),y1=f (x1)线性插值就是构造一个一次多项式
P1(x) = ax + b
使它满足条件
P1 (x0) = y0 ,P1 (x1) = y1
其几何解释就是一条直线,通过已知点A (x0, y0),B(x1, y1)。
线性插值计算方便、应用很广,但由于它是用直线去代替曲线,因而一般要求[x0, x1]比较小,且f(x)在[x0, x1]上变化比较平稳,否则线性插值的误差可能很大。为了克服这一缺点,有时用简单的曲线去近似地代替复杂的曲线,最简单的曲线是二次曲线,用二次曲线去逼近复杂曲线的情形。
具体可见:
http://blog.csdn.net/jiangnanyouzi/article/details/5698122
python 实现
from scipy.interpolate import lagrangex = [3, 6, 9]y = [10, 8, 4]lagrange(x,y)#poly1d([ -0.11111111, 0.33333333, 10. ])
以上 lagrange(x,y) 的输出值 poly1d([−0.11111111,0.33333333,10.]) 值的是多项式的三个系数
即: a0=−0.11111,a1=0.3333333,a2=10.
如果要进行插值操作,可以:
lagrange(x, y)(10)# 2.222222
0 0
- 拉格朗日插值——基于scipy实现
- Scipy 基础 —— 稀疏矩阵
- SciPY——数值计算库
- SciPy——数值计算库
- 拉格朗日插值 python scipy
- SciPy
- SciPy
- SciPy
- SciPy
- scipy.special —— 排列、组合与阶乘
- Python——解决scipy导入出错问题
- scipy —— 丰富的子包(io、cluster)
- 【Python】三、Scipy——《用Python做科学计算》
- Python笔记——Python基础、numpy基础、SciPy、Matplotlib
- python科学计算之scipy——optimize
- Windows、Linux安装pandas,scipy,sklearn(scikit—learn)
- scipy.cluster模块实现kmeans聚类
- 最小二乘拟合-python scipy实现
- HTML 内容隐藏显示以及display和visibility的对比
- JS数组的常用方法总结
- 学习Siki学堂C#初级教程、中级教程过程中的笔记和最后完成的自定义列表MyList
- 草根程序员转型做项目管理走过的点点滴滴之"经验总结"
- 不删该任何文件解决 “/usr/lib64/libstdc++.so.6: version `GLIBCXX_3.4.14' not found “
- 拉格朗日插值——基于scipy实现
- C++异常与继承
- Piggy-Bank [dp][完全背包]
- hdu3853 LOOPS 【概率DP+期望问题】
- JS字符串的一些常用方法
- 简单的单向链表的创建(总结)
- 第一次写博客,记录一下今天遇到的一个错误
- 递归之算24
- mahout in action推荐系统代码解释