python之拉格朗日插值算法完成计算

一、提出问题

拉格朗日插值法，如何计算数值？

二、详解拉格朗日插值算法

插值法就是利用已知的点(x,y)建立合适的插值函数y=f(x),将未知点w带入插值函数 f(x)可以求出函数值 f(w)。对于平面上相异（无两点在一条直线上）的 n个点，我们可以找到一个 n次多项式

使这个多项式函数经过这些点。有n+1个未知数(A)，我只需要知道n+1个点的函数值就可以解出这n+1个未知数,拉格朗日插值法就是求得这个多项式函数。
1、拉格朗日插值法
(1)将n+1个点代入多项式得

(2)方程的系数矩阵为A的行列式为

三、举例说明

1、要求
根据上述基本原理，用python程序语言将已知多点数据画出图像，并利用这组数值计算一组数据[-3,2,3,-1]对应的拉格朗日插值，并将其形成的多个对应点在图像中标出如图1
2、具体代码

import pylab as pl
import numpy as np def lagrange(x,y,w):n = len(x)res=0for i in range(n):temp = 1 for j in range(n):if i!=j: temp = temp * (w-x[j]) / (x[i] - x[j]) res += temp * y[i]return resx = [16,9,25]
y = [4,3,5]
lagrange(x,y,10)   x1=np.linspace(-5,5,20)
y1= x1**2+np.random.uniform(-1,1,20)
x2=[-3,2,3,-1]
y2=[]
for i in x2:y2.append(lagrange(x1,y1,i))
pl.plot(x1,y1, 'r')
pl.scatter(x2,y2)
pl.xlim(-6,6)
pl.ylim(-5,30)
pl.show()

3、运行结果

以上就是本文所有内容，希望能帮到大家！！！

python之拉格朗日插值算法完成计算相关推荐

Python实现拉格朗日插值算法
关于拉格朗日插值算法的定义可以查看维基百科拉格朗日插值法. 这里直接上代码,当前代码不是最优算法 class Algorithm:@staticmethoddef LagrangeInterpolat ...
函数插值的python实现——拉格朗日、牛顿插值
函数插值的python实现--拉格朗日.牛顿插值 1. 拉格朗日(Larange)插值实现代码 2. 牛顿(Larange)插值 2.1 牛顿插值多项式的基本形式 2.2 牛顿均差插值多项式 (1) ...
python 科学计算基础教程电子版-自学Python 编程基础、科学计算及数据分析
自学Python 编程基础.科学计算及数据分析 epub pdf mobi txt 下载自学Python 编程基础.科学计算及数据分析 epub pdf mobi txt 下载 ☆☆☆☆☆ 李金著 ...
全面对比 MATLAB、Julia、Python，谁在科学计算中更胜一筹？
数百种编程语言,各有优劣,各自也都有自己最为适用的场景.那么就科学计算领域而言,主流的 MATLAB.Julia.Python 会有哪些最为独特的优势呢?又存在哪些让开发者无力的缺陷?在本文中,我们将 ...
梯度、梯度法、python实现神经网络的梯度计算
[机器学习]梯度.梯度法.python实现神经网络的梯度计算一.python实现求导的代码: 二.what is 梯度三.使用梯度法寻找神经网络的最优参数四.神经网络的梯度计算一.python ...
Python数据分析学习笔记：计算向量夹角
Python数据分析学习笔记:计算向量夹角通过计算两个向量夹角了解两个向量之间的相关性. # 计算向量夹角import numpy as npdef included_angle(a, b):a_n ...
Python数据分析学习笔记：计算相关系数
Python数据分析学习笔记:计算相关系数 1.相关系数概念相关系数,或称线性相关系数.皮氏积矩相关系数(Pearson product-moment correlation coefficient ...
python怎么算积分_python计算积分
python有多个方法计算积分,下面介绍其中三个,以下式为例: 方法一:直接用numpy计算 start = 1 stop = 2 length = 101 x = np.linspace(start ...
python修改ABAQUS的inp计算文件
python修改ABAQUS的inp计算文件这里的函数可以与上一个博客Matlab-python-ABAQUS数据交互及联合使用的函数接口相互配合和集成,实现matlab与Abaqus进行联合仿真计 ...