python计算对角线有理函数插值的方法


Posted in Python onMay 07, 2015

本文实例讲述了python计算对角线有理函数插值的方法。分享给大家供大家参考。具体实现方法如下:

''' p = rational(xData,yData,x)
  Evaluates the diagonal rational function interpolant p(x)
  that passes through he data points
'''  
from numpy import zeros
def rational(xData,yData,x):
  m = len(xData)
  r = yData.copy()
  rOld = zeros(m)
  for k in range(m-1):
    for i in range(m-k-1):
      if abs(x - xData[i+k+1]) < 1.0e-9:
        return yData[i+k+1]
      else:
        c1 = r[i+1] - r[i]
        c2 = r[i+1] - rOld[i+1]
        c3 = (x - xData[i])/(x - xData[i+k+1])
        r[i] = r[i+1] + c1/(c3*(1.0 - c1/c2) - 1.0)
        rOld[i+1] = r[i+1]
  return r[0]

希望本文所述对大家的Python程序设计有所帮助。

Python 相关文章推荐
Python中使用PyHook监听鼠标和键盘事件实例
Jul 18 Python
快速入手Python字符编码
Aug 03 Python
python 中的int()函数怎么用
Oct 17 Python
Windows环境下python环境安装使用图文教程
Mar 13 Python
详解Python中的分组函数groupby和itertools)
Jul 11 Python
pyinstaller打包多个py文件和去除cmd黑框的方法
Jun 21 Python
Python 点击指定位置验证码破解的实现代码
Sep 11 Python
python爬虫开发之使用python爬虫库requests,urllib与今日头条搜索功能爬取搜索内容实例
Mar 10 Python
Python CSS选择器爬取京东网商品信息过程解析
Jun 01 Python
PySide2出现“ImportError: DLL load failed: 找不到指定的模块”的问题及解决方法
Jun 10 Python
python实现二分查找算法
Sep 18 Python
如何创建一个Flask项目并进行简单配置
Nov 18 Python
python计算方程式根的方法
May 07 #Python
python计算牛顿迭代多项式实例分析
May 07 #Python
用C++封装MySQL的API的教程
May 06 #Python
使用wxPython获取系统剪贴板中的数据的教程
May 06 #Python
用Python遍历C盘dll文件的方法
May 06 #Python
使用Python压缩和解压缩zip文件的教程
May 06 #Python
Python发送以整个文件夹的内容为附件的邮件的教程
May 06 #Python
You might like
全国FM电台频率大全 - 25 云南省
2020/03/11 无线电
php+ajax实现无刷新数据分页的办法
2015/11/02 PHP
php支付宝APP支付功能
2020/07/29 PHP
一个JavaScript获取元素当前高度的实例
2014/10/29 Javascript
JavaScript模拟实现继承的方法
2015/03/30 Javascript
微信jssdk用法汇总
2016/07/16 Javascript
移动端使用localStorage缓存Js和css文的方法(web开发)
2016/09/20 Javascript
jQuery右下角悬浮广告实例
2016/10/17 Javascript
AngularJS中如何使用echart插件示例详解
2016/10/26 Javascript
jQuery Ajax前后端使用JSON进行交互示例
2017/03/17 Javascript
jQuery实现打开网页自动弹出遮罩层或点击弹出遮罩层功能示例
2017/10/19 jQuery
Node.js创建HTTP文件服务器的使用示例
2018/05/11 Javascript
Easyui 关闭jquery-easui tab标签页前触发事件的解决方法
2019/04/28 jQuery
解决python2.7用pip安装包时出现错误的问题
2017/01/23 Python
Python自定义函数实现求两个数最大公约数、最小公倍数示例
2018/05/21 Python
pandas实现将日期转换成timestamp
2019/12/07 Python
Python中zip()函数的解释和可视化(实例详解)
2020/02/16 Python
Django使用Celery加redis执行异步任务的实例内容
2020/02/20 Python
Python坐标轴操作及设置代码实例
2020/06/04 Python
如何用python 操作zookeeper
2020/12/28 Python
jupyter notebook远程访问不了的问题解决方法
2021/01/11 Python
科颜氏加拿大官方网站: Kiehl’s加拿大
2016/08/16 全球购物
PUMA官方商城:世界领先的运动品牌之一
2016/11/16 全球购物
ALDI奥乐齐官方海外旗舰店:德国百年超市
2017/12/27 全球购物
告诉你怎样写创业计划书
2014/01/27 职场文书
校园环保标语
2014/06/13 职场文书
军训拉歌口号
2014/06/13 职场文书
工会趣味活动方案
2014/08/18 职场文书
夫妻分居协议书范文
2014/11/26 职场文书
优秀护士事迹材料
2014/12/25 职场文书
大学生考试作弊被抓检讨书
2014/12/27 职场文书
2015年房产经纪人工作总结
2015/05/15 职场文书
党校团干班培训心得体会
2016/01/06 职场文书
心得体会格式及范文
2016/01/25 职场文书
预备党员的思想汇报,你真的会写吗?
2019/06/28 职场文书
《敬重卑微》读后感3篇
2019/11/26 职场文书