def func(x,y,X,infor=True):

list2=[y[0]] # 差商表的对角线的第一个元素始终是y0

count=1

while(True):

if len(y)>1:

list=[] # 空列表用来保存,每次计算后差商表的行

for i in range(len(y)-1):

n=x[i+count]-x[i]

m=y[i+1]-y[i]

l=m/n

list.append(l)

list2.append(list[0]) # list2用来记录差商表的对角线元素,每计算一次,取行的第一个元素

count += 1

y = list

else:

break

if infor: # 判断是否要继续计算,结果

W=0

for i in range(len(list2)):

if i==0:

w=list2[i]

else:

w = list2[i]

for j in range(i):

w*=(X-x[j])

W+=w

print(‘牛顿插值:‘, W)

return ‘牛顿差商表对角线列:%s‘ %list2

ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],‘‘,infor=False)

print(ret)

ret=func([0.32, 0.34, 0.36],[0.314567, 0.333487, 0.352274],0.3367)

print(ret)

运行结果:

原文:http://blog.51cto.com/13747953/2308551

牛顿插值法python代码_Python实现牛顿插值法(差商表)相关推荐

  1. 牛顿迭代法Python代码,全网最详细,教学向

    代码功能包括函数图像展示,初始值选取收敛区间判断,迭代结果输出,迭代过程图像输出. 因讲解过于冗长,先将完整代码直接放在这里,只是想抄个模板方便修改的可以直接拿去用啦,有不了解的地方可以再翻下去看. ...

  2. postman生成python代码_python代码直接生成可执行exe文件

    前言: python中有一个很厉害的库可以直接将写好的python代码打包成可执行的.exe文件,生成好的文件不需要python运行环境和pycharm等IDE就可以直接双击运行,是不是很厉害呢!这样 ...

  3. 初学者怎样看懂python代码_Python 初学者,如何快速的入门

    随着疫情的到来,不少人开启了退休之前最长的一次度假,唯一不同的是,这次的度假是在家里完成的.同时还面临着低薪,高额的房租房贷卡贷等. 所以就有了油然而生的想法,何不在这唯一一次长时间的假期来学习一项技 ...

  4. word转换pdf的python代码_python word转pdf代码实例

    原理 使用python win32 库 调用word底层vba,将word转成pdf 安装pywin32 pip install pywin32 python代码 from win32com.clie ...

  5. 一段有趣的python代码_Python 五个有趣的彩蛋,你都知道吗?

    阅读文本大概需要 6.6 分钟. 当一门编程语言是开源的时候,往往会有产生一些搞笑和有趣的东西.通常,这意味着社区的贡献者会为该语言添加一些有趣和特别的彩蛋以及隐藏的特性(当然前提是不会增加在生产环境 ...

  6. 海绵宝宝python代码_Python那些事——python还能这么玩?帅呆了

    当我们还是孩子的时候,总是迷恋漫画书,当翻到我们的英雄们开始行动时会激动不已. 大家总是争论谁是最厉害的超级英雄,认真地讨论他们的家族历史,或者梦想自己拯救高谭市.我们很多人用自己的童年创造了这样的难 ...

  7. 拉格朗日插值python代码_Python实现的拉格朗日插值法示例

    本文实例讲述了Python实现的拉格朗日插值法.分享给大家供大家参考,具体如下: 拉格朗日插值简单介绍 拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法. 许多实际问题中都 ...

  8. 缺失值处理 - 拉格朗日插值法 - Python代码

    目录 缺失值处理 拉格朗日差值法的理论基础 拉格朗日插值法代码实现 其他数据预处理方法 缺失值处理 处理缺失值常用的办法可分为三类:删除记录.数据插补.不处理. 其中常见的数据插补法有: 如果通过删除 ...

  9. 拉格朗日插值法——python代码实现

    一.拉格朗日的基本思想: 二.线性插值 三.多个点 四.代码实现 def lagrange(xx,y):l=len(y)l_n = 0for k in range(l):xxx=xx.copy()x_ ...

  10. 高斯牛顿迭代法matlab代码,优化算法--牛顿迭代法

    简书同步更新 牛顿法给出了任意方程求根的数值解法,而最优化问题一般会转换为求函数之间在"赋范线性空间"的距离最小点,所以,利用牛顿法去求解任意目标函数的极值点是个不错的思路. 方程 ...

最新文章

  1. rsync 远程同步——安全高效的异地备份策略
  2. hdu5823 (附带数的二进制子集)
  3. js符号输入不可用_JS 控制非法字符的输入代码
  4. java演练 数组的逆序文字玩法 你是猪才怪
  5. android6.0原生brower_Android Browser学习一 application的初始化
  6. 二叉树:一入递归深似海,从此offer是路人
  7. Redis的高级应用-安全性和主从复制
  8. 2016西安教师职称计算机考试,2016年教师职称计算机考试题库.doc
  9. 微信小程序分享/转发功能方法
  10. 米家电磁炉显示e10_MIJIA 米家 电磁炉
  11. 如何建立风险分析模型
  12. oracle数据库学习笔记(二)(创建表、DDL、DML、单行插入数据、多行插入数据)
  13. Latex改变页边距
  14. INS/GNSS组合导航(十二)如何读懂MEMS惯性器件精度?
  15. 用cd 命令进入和退出D盘文件夹的操作步骤。
  16. 微电子以后什么方向前途好?FPGA吗?
  17. 多读书多看报...... http://www.it-ebooks.info/
  18. c语言把字符串复制到剪贴板,如何把位图拷贝到剪贴板中
  19. 小法师系列-8-番外-失恋的并发症
  20. openjudge 1.8.5 计算鞍点

热门文章

  1. think-cell 数据表无法打开怎么办
  2. 电机控制park变换公式推导
  3. 单系统 台电x80pro_台电X80HD安装Win8单系统教程
  4. android通电自动开机,【图】改造导航仪,实现通电自动开机!
  5. 2000年建模b题matlab,2013数学建模B题matlab代码
  6. 短视频APP管理系统源码 直播系统源码
  7. selenium 谷歌驱动 下载网址
  8. android gps free,GPS HUD Speedometer Free
  9. CopyU!v2 已经收录到腾讯软件管家!
  10. 苹果CoreFoundation源代码