今天给大家分享的是关于用python实现牛顿插值的内容，用牛顿插值算法能帮助更快速的实现n阶差商的算术，本下文有实例和详细注释供大家参考，需要的朋友参考，接下来跟随小编来学习一下吧。

    1. n阶差商实现

def diff(xi,yi,n):
    """
    param xi:插值节点xi
    param yi:插值节点yi
    param n: 求几阶差商
    return: n阶差商
    """
    if len(xi) != len(yi):  #xi和yi必须保证长度一致
        return
    else:
        diff_quot = [[] for i in range(n)]
        for j in range(1,n+1):
            if j == 1:
                for i in range(n+1-j):
                    diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1]))
            else:
                for i in range(n+1-j):
                    diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j]))
    return diff_quot

    测试一下：

xi = [1.615,1.634,1.702,1.828]
yi = [2.41450,2.46259,2.65271,3.03035]
n = 3
print(diff(xi,yi,n))

    返回的差商结果为：

[[2.53105263157897, 2.7958823529411716, 2.997142857142854], [3.0440197857724347, 1.0374252793901158], [-9.420631485362996]]

    2. 牛顿插值实现

def Newton(x):
    f = yi[0]
    v = []
    r = 1
    for i in range(n):
        r *= (x - xi[i])
        v.append(r)
        f += diff_quot[i][0] * v[i]
    return f

    测试一下：

x = 1.682
print(Newton(x))

    结果为：

2.5944760289639732

    3.完整Python代码

def Newton(xi,yi,n,x):
    """
    param xi:插值节点xi
    param yi:插值节点yi
    param n: 求几阶差商
    param x: 代求近似值
    return: n阶差商
    """
    if len(xi) != len(yi):  #xi和yi必须保证长度一致
        return
    else:
        diff_quot = [[] for i in range(n)]
        for j in range(1,n+1):
            if j == 1:
                for i in range(n+1-j):
                    diff_quot[j-1].append((yi[i]-yi[i+1]) / (xi[i] - xi[i + 1]))
            else:
                for i in range(n+1-j):
                    diff_quot[j-1].append((diff_quot[j-2][i]-diff_quot[j-2][i+1]) / (xi[i] - xi[i + j]))
    print(diff_quot)
    
    f = yi[0]
    v = []
    r = 1
    for i in range(n):
        r *= (x - xi[i])
        v.append(r)
        f += diff_quot[i][0] * v[i]
    return f

    关于用python实现牛顿插值法的实例就分享到这，上述代码仅供参考，感兴趣的朋友可以看看，希望能对大家有帮助，想要了解更多python插值法的内容，大家可以关注群英网络其它相关文章。

文本转载自脚本之家

免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：mmqy2019@163.com进行举报，并提供相关证据，查实之后，将立刻删除涉嫌侵权内容。

猜你喜欢

• 用Python怎样绘制动态的地铁图，实现过程是什么

  这篇文章主要介绍了用python制作北上广深——地铁线路动态图，文中的示例代码讲解详细，对我们的工作或学习都有一定的价值，感兴趣的同学可以学习一下

• 用Python如何爬取疫情数据，并进行数据可视化操作

  这篇文章主要介绍了Python利用爬虫爬取疫情数据并进行可视化的展示，文中的示例代码讲解清晰，对工作或学习有一定的价值，需要的朋友可以参考一下

• torchtext库如何实现导入NLP数据集，流程是什么

  torchtext在文本数据预处理方面特别强大，但我们要知道ta能做什么、不能做什么，并如何将我们的需求用torchtext实现。虽然torchtext是为pytorch而设计的，但是也可以与keras、tensorflow等结合使用。

• 如何用matplotlib画XKCD风格的图表

  内容介绍XKCDmatplotlib对XKCD风格的支持创建XKCD风格的图表使用中文字体创建XKCD风格的图表1.安装字体2.更新中文字体XKCDxkcd是兰道尔门罗（Randall

• pytorch visdom怎样安装，用法是什么？

  pytorch Visdom可视化,是一个灵活的工具,用于创建,组织和共享实时丰富数据的可视化。这篇文章我们就来了解pytorch visdom怎样安装，以及visdom的用法，感兴趣的接下来就跟随小编来了解看看吧！