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中assert用法实例分析
Apr 30 Python
Python正则抓取新闻标题和链接的方法示例
Apr 24 Python
利用Python破解斗地主残局详解
Jun 30 Python
python 实现上传图片并预览的3种方法(推荐)
Jul 14 Python
CentOS7.3编译安装Python3.6.2的方法
Jan 22 Python
python实现淘宝秒杀聚划算抢购自动提醒源码
Jun 23 Python
Django外键(ForeignKey)操作以及related_name的作用详解
Jul 29 Python
在python中使用pymysql往mysql数据库中插入(insert)数据实例
Mar 02 Python
pycharm远程连接vagrant虚拟机中mariadb数据库
Jun 05 Python
python文件及目录操作代码汇总
Jul 08 Python
Python实现淘宝秒杀功能的示例代码
Jan 19 Python
python批量提取图片信息并保存的实现
Feb 05 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
PHP设计模式之代理模式的深入解析
2013/06/13 PHP
php输出金字塔的2种实现方法
2014/12/16 PHP
php用正则判断是否为数字的方法
2016/03/25 PHP
thinkPHP的表达式查询用法详解
2016/09/14 PHP
php实现36进制与10进制转换功能示例
2017/01/10 PHP
php+redis在实际项目中HTTP 500: Internal Server Error故障排除
2017/02/05 PHP
php 使用ActiveMQ发送消息,与处理消息操作示例
2020/02/23 PHP
IE JS无提示关闭窗口不提示的方法
2010/04/29 Javascript
javascript instanceof 内部机制探析
2010/10/15 Javascript
location.href用法总结(最主要的)
2013/12/27 Javascript
判断某个字符在一个字符串中是否存在的js代码
2014/02/28 Javascript
Nodejs进程管理模块forever详解
2014/06/01 NodeJs
基于jquery的文字向上跑动类似跑马灯的效果
2014/09/22 Javascript
基于jQuery.Hz2Py.js插件实现的汉字转拼音特效
2015/05/07 Javascript
关于在Servelet中如何获取当前时间的操作方法
2016/06/28 Javascript
JS实现点击网页判断是否安装app并打开否则跳转app store
2016/11/18 Javascript
jQuery图片轮播(二)利用构造函数和原型创建对象以实现继承
2016/12/06 Javascript
JS实现密码框的显示密码和隐藏密码功能示例
2016/12/26 Javascript
Node.js服务器开启Gzip压缩教程
2017/08/11 Javascript
在Angular中使用JWT认证方法示例
2018/09/10 Javascript
微信小程序实现横向滚动导航栏效果
2019/12/12 Javascript
vue-cli4项目开启eslint保存时自动格式问题
2020/07/13 Javascript
vue下载二进制流图片操作
2020/10/26 Javascript
深入了解Vue3模板编译原理
2020/11/19 Vue.js
[05:48]DOTA2英雄梦之声vol21 屠夫
2014/06/20 DOTA
python套接字流重定向实例汇总
2016/03/03 Python
python3 实现验证码图片切割的方法
2018/12/07 Python
Python如何把字典写入到CSV文件的方法示例
2020/08/23 Python
日本著名的平价时尚女性购物网站:Fifth
2016/08/24 全球购物
新闻专业应届生求职信
2013/10/31 职场文书
文明礼仪事迹材料
2014/01/09 职场文书
创意广告词
2014/03/17 职场文书
党风廉正建设个人工作总结
2015/03/06 职场文书
2015年机关党委工作总结
2015/05/23 职场文书
2015年幼儿园教育教学工作总结
2015/05/25 职场文书
2016学习依法治国心得体会
2016/01/15 职场文书