插值

20 0 8年第 2期

建电

三次 H r i插值的扩展及其应用 em t e

洪柏娣

(江传媒学院浙江杭州 3 0 1 浙 10 8)

【摘要1:在总结分析三次 H r i em t e插值多项式的基础上,对三次 H r i em t e插值公式进行了推广和扩展。通过改变三次 H r t插值的初始条件。到了扩展的插值多项式计算公式。给出了扩展的三次 H r i emi e得 em t值格式的有理函数的近似表示 e插

方法以及有理函数的数值积分算法。 【关键词】:三次 Hem t插值近似表示数值积分 ri e

1引言 .

(), (= () () O=1。1 0= ) 0=0

插值法作为数值分析的理论基础。着十分丰富的内容。有它包括 l rne插值、 e tn值、 em t值、分段低次插值、 a ag g N wo插 H ri e插 样条插值等。它们均在数值计算领域发挥着不同程度的作用。

(), ( ) ( ) ( ) l0= 0= 0=0 l 1=1

,

(),。0= o1= () 0=1 () ( 0=0妒 )

7 (),o0:7(=, () 7 0=l7() 7 1 7 0=0: 7。) o经典的插值方法以 T y r值和 L gag al插 o arne插值为代表。 满足这些条件的三次多项式容易构造出来。比如,妒由 T o插值利用函数在定义域内某点处的 0阶至 n阶导数信息 alr 0知妒,缸 B d利用 卜 J卢D a -,=。所以, 得= Id l给出复杂函数或未知函数的近似多项式表达式,arne插值 L gag利用多个离散点的函数值信息给出函数的近似多项式表达式, % () x+1 x=一 () 8进一步根据插值结果对与复杂函数或未知函数相关的理论和应同理可得妒f x 3 j () 9用问题作出讨论。因此 T yo插值和 Lga g插值有着紧密的 alr arn e丸() - =-+ X (0 1)

联系 . a l插值可以看做 L ga g插值的极限形式; a r g Ty r o ar e n L ga e n

1

1 .

插值则是 T y r值的离散化形式。 al插 o 在某些实际应用问题中 .为了保证插

值函数能更好的密合原来的函数 .但要求插值函数”点”即插值函数和被插函数不过,在节点上具有相同的函数值 .且要求”切”即两者在节点处而相 .还具有相同的导数值,类插值称作切触插值,者称为 He.这或 r

(一 )

+ r

( 1 1)

访照三次 H r i插值多项式的推导过程可得: em t e

( x )

)

)嘞

)

( 2 1)

m t插值。由于 T yo插值利用的是”点” i e alr一的各阶导数信息, L gag ar e插值利用的” n多点”数值信息 . H r t值即利用函而 emi e插函数值信息又利用导数信息。所以 H r i插值是 T y r值 em t e al插 o和 L g n e值的综合和推广。这里提出的一种新型的三次 ar g插 a

H r i插值多项式是对三次 H r t值公式的推广和扩展。 em t e emi e插

将 ( ) ( 1式代入 (2式得: 8一 1) 1)

由此得到扩展的三次 Hem t插值多项式。 ri e

结合实际数值计算验证了这种新型的三次 Hem t值多项式 3扩展的三次 h r i ri e插 . em t e插值多项式的应用 31理函数的近似表示 .有的正确性。

2三次 H r i . em t值的扩展 e插

例:定函数给

21次 He i插值多项式及其初始条件 .三 r t m e令三次多项式与() %() ()唬+萌 ( )=++ ) () 1式中相应的妒咖为待定的三次式,别满足条件: D 分 应用扩展的三次 H r t emi e插值多项式求对应的插值多项式 P。 j

厂)南 ( x

妒o 0= l妒 ( )妒 ( ) ( ), o 1= 0=妒 ( ) 0 1=

妒(=l I0=妒() ) I )。 ():=妒(=0 1妒 0 1

庐 ( ) . o 0=庐 () 1=0 0=1庐 ( ) o1=庐 () (), ( ) () 。 ) 1=l 0= 1= (=0 0

() 2

…隅㈤ 离 南

取— J,,到相应的值如下表:,1得 0

1 0 1

-

根据文献 l的叙述 .易推导出:容

) Yg+ ̄ t

)+

)+

() 3

. ,∽

l

0

l

其中 D J任意两点,= J勘。是 h x—

由( ) ( )知: 2及 3式

2

2

/ ) O

0

1

0

P( ) Y, ), x y, ̄ o= o x (= g( 0=,

)=

f、 4

厂(

一 l

2

0

表 1

一一 l

2

我们把()称为三次 H ri 4式 em t值多项式 B的初始条 e插件。 2 .展的三次 He t插值多项式 2扩 r e mi改变三次 H ri em t值的初始条件 .重新构造三次 H r i e插 em t。插值多项式 )其满足,使

( ) Y ( ) Y ( ) Y, ( ) Y O, O, O l ( )

分别代入 (3式得: 1)

'

令尸(铴()锕( hoo ) h 3= x++ y ̄(+  ̄x

式中%妒

(

() 6

为待定的三次多项式,别满足条件分

三次hermite插值matlab,三次hermite插值相关推荐

  1. matlab循环语句从小到大排列,三个例子讲解MATLAB三种循环

    三个例子讲解MATLAB三种循环 FOR循环 在for和end语句之间的{commands}按数组中的每一列执行一次.在每一次迭代中,x被指定为数组的下一列,即在第n次循环中,x=array(:, n ...

  2. matlab三次拟合多项式,matlab三次多项式拟合

    暂无评价|0人阅读|0次下载|举报文档 DSC曲线的三次样条拟合与Matlab实现_能源/化工_工程科技_专业资料.DSC曲线拟合与Matlab实现第... 人阅读|次下载 曲线拟合的最小二乘法mat ...

  3. matlab中三点定圆,MATLAB三点确定一个圆

    1 给定三个不在同一直线上的三个点的坐标,确定圆的轨迹方程 2 3 4 下面是MATLAB中的函数文件 5 function[a,b,c]=ThreePointFitCircle(x,y,flag) ...

  4. opengl绘制三次hermite曲线,三次cardinal曲线

    引言 上篇 上一篇绘制了最基础插值曲线,讲到除了插值经过给定点之外,还有别的约束条件.有时为了所得曲线在端点处具有给定的N阶导数,亦即所称的参数连续性条件和几何连续性条件,给约束矩阵添加了导数约束. ...

  5. 艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值

    Hermite插值满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值.对于高阶导数的情况,Hermite插值多项式比较复杂,在实际情况中,常常遇到的是函数值与一阶导数给定的情况.在此情 ...

  6. hermit插值 matlab,埃尔米特(Hermite)插值及其MATLAB程序

    %hermite.m %求埃尔米特多项式和误差估计的MATLAB主程序 %输入的量:X是n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标向量,Y是纵坐标向量, %以f'(x ...

  7. matlab三次多样式对函数拟合,三次多样式样条函数债券

    1.如何用编程实现三次样条函数 用MatLab 2.如何用Excel进行三次样条函数拟合 安照要求先制作数据表,再插入图表----曲线(也就是折线)系统会自动生成图表的,再设置一下效果! 3.三次样条 ...

  8. 图形图像处理-之-高质量的快速的图像缩放 中篇 二次线性插值和三次卷积插值

    from:http://blog.csdn.net/housisong/article/details/1452249 图形图像处理-之-高质量的快速的图像缩放 中篇 二次线性插值和三次卷积插值    ...

  9. matlab测量液体液位,基于MATLAB三容水箱液位控制系统.doc

    基于MATLAB三容水箱液位控制系统 基于MATLAB三容水箱液位控制系统 计算机技术+沈瑞清+12011112262 本文通过对三容水箱的液位定值控制来说明SIMULINK在自动控制系统中的应用. ...

  10. 卫星轨道的估计问题(Matlab)(三):标准重采样粒子滤波(SIR)对新问题的尝试

    SIR滤波器 关于粒子滤波的基本知识可以参加下面的博客: 粒子滤波,讲的很通俗易懂 基本粒子滤波算法过程: SIR算法伪代码: Matlab代码实现 求解如下所示的滤波问题: %PF的应用 clc,c ...

最新文章

  1. CUDA support is required (OpenCV CMake parameter ‘WITH_CUDA‘ must be true).
  2. JavaScript的格式--从格式做起,做最严谨的工程师
  3. LwIP应用开发笔记之九:LwIP无操作系统TELNET服务器
  4. 一个div 上下两行_web前端工程师如何理解 CSS 布局和块级格式化上下文
  5. Linux内核RCU(Read Copy Update)锁简析-前传
  6. DataSet/DataFrame性能比RDD高?
  7. deepin linux 怎么安装软件,Linux Deepin 从 Backports 安装软件包
  8. Oracle 购买价格 和 服务费 计算方式
  9. 教你制作专属个性二维码!将照片制作成二维码图片!给二维码添加背景图片!
  10. TeamViewer 如何注册账户?
  11. 考勤记录查询java_kaoqin 一个员工考勤软件,用java和sql写的,实现基本的签到和查询功能 Develop 238万源代码下载- www.pudn.com...
  12. TCP/IP 主要报文头格式
  13. Vue 学习——表格列表数据显示百分号%、显示数字改为自定义内容、数据表头下拉、点击筛选数据
  14. html转pdf 图片跨域问题解决
  15. 一周新论文 | 2020年第12周 | 自然语言处理相关
  16. 动手学数据分析-数据可视化
  17. python运算结果是fasls的代码_《利用Python进行数据分析》--整理学习
  18. linux日志切割命令,Linux 服务器log日志切割三种方法【附命令行】
  19. 如何获取新浪股票行情查询接口?
  20. 2021最新Android知识体系总结,实战篇

热门文章

  1. 职业学校计算机教学,职业学校计算机专业教学初探
  2. 修改PC端微信提示音
  3. 卡方拟合优度检验怎么做?
  4. javascript毫秒计时器_JavaScript原生秒表、计时器
  5. 【虚拟光驱怎么安装Win7系统】
  6. 0基础学习3dmax游戏建模有这4个技巧
  7. 网络编程-UDP编程
  8. java hsqldb_HSQLDB: java程序使用hsqldb 入门教程 java启动hsqldb
  9. sql if语句实例
  10. 语音信号处理(赵力)作业答案第8章——语音合成