怎样用python实现牛顿插值?
Admin 2021-09-08 群英技术资讯 783 次浏览
今天给大家分享的是关于用python实现牛顿插值的内容,用牛顿插值算法能帮助更快速的实现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]]
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
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和Pygame实现24点小游戏,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
这篇文章主要介绍了Python中schedule模块关于定时任务使用方法,文章围绕主题展开详细的内容介绍,具有一定的参考价值,需要的小伙伴可以参考一下
这篇文章主要介绍了如何利用 K-Means 聚类进行色彩量化,以减少图像中颜色数量。文中的代码具有一定的学习价值,感兴趣的小伙伴可以关注一下
变量是拥有匹配对象的名字(标识符),命名空间是一个包含了变量名称们(键)和它们各自相应的对象们(值)的字典,下面这篇文章主要给大家介绍了关于Python基础教程之名称空间以及作用域的相关资料,需要的朋友可以参考下
Python中字符串反转常用的五种方法:使用字符串切片、使用递归、使用列表reverse()方法、使用栈和使用for循环。
成为群英会员,开启智能安全云计算之旅
立即注册Copyright © QY Network Company Ltd. All Rights Reserved. 2003-2020 群英 版权所有
增值电信经营许可证 : B1.B2-20140078 粤ICP备09006778号 域名注册商资质 粤 D3.1-20240008