hermite插值c语言程序,张艳-埃尔米特Hermite 插值逼近的C语言程序.doc
张艳-埃尔米特Hermite 插值逼近的C语言程序.doc
论文题目:埃尔米特(Hermite) 插值逼近的C语言程序
院 系: 数学科学学院
专 业: 数学与应用数学
姓 名: 张 艳
学 号:
指导教师: 侯 存 贵
完成时间: 2007-5-15
埃尔米特(Hermite) 插值逼近的C语言程序
张艳
包头师范学院数学科学学院
摘要:本文主要探讨埃尔米特(Hermite) 插值逼近的C语言程序算法,着重分析其推导过程,并给出了其C语言程序以及埃尔米特(Hermite)插值逼近的简单应用.
关键词:Hermite插值多项式;插值条件; Hermite插值基函数.
一、Hermite插值多项式定义
定义: 设基点互异.给定,,.要求插值多项式满足 (1)
则称为二重密切Hermite插值多项式,简称为Hermite插值多项式.称为二重插值基点.
(1)式共有个条件,因此Hermite插值多项式通常次数不超过,故可将记为.
二、埃尔米特(Hermite)插值多项式的存在唯一性
定理:关于互异基点满足条件(1)的二重密切次Hermite插值多项式存在且唯一.
证明:设有2n+1次多项式
(2)
满足条件(1)即 (3)
由(2)式知(3)式是一个关于的阶线性方程组.的存在唯一性决定于(3)式为齐次线性方程系组,即当=0,=0()时,(3)式仅有平凡解 .现用反证法证明:若齐次方程组有非平凡解,则表示存在一个次数不高于的多项式满足
为的互异的二重零点,即次多项式有个零点(包括重数),这和代数基本定理相矛盾.
三、埃尔米特(Hermite)插值多项式的构造
由定理知存在且唯一,我们用类似于拉格朗日插值多项式的构造方法来构造Hermite插值多项式
设,,,分别满足插值条件
(4)
(5)
(其中表示克罗内克(Kronecker)符号.当, 时=1;当,时=0.)的次多项式,于是次数不超过次多项式 能够满足插值条件 .因而就是所要求的次多项式.因此只要构造出满足条件(4)(5)的和即可.我们把满足插值条件(4)(5)的次多项式和()称为Hermite插值基函数.下面构造和
由于关于基点的拉格朗日基函数满足 ( ,)且是次多项式,结合插值条件(4)设为,则依据条件(4)要求应有
当时,由于 故
当时应有
从中解出 可得
从而得到
设为 则
依据条件(5)要求应有
当时,由于 故
当时,故也应有 而
故 即
因此,我们得到埃尔米特插值函数的基函数为
根据插值条件,利用二重密切的Hermite插值基函数的性质,Hermite插值多项式可简单地表示为
四、埃尔米特(Hermite)插值多项式误差
在求解某些数学问题时,用有限的过程代替无限过程所产生的误差称为截断误差(或方法误差).
定理:a:设的导数于[a,b]连续,于(a,b)内存在, ()互异;
b:为Hermite插值多项式;
则 .
其中与有关。
五、埃尔米特(hermite)插值程序
功能:给定个基点,()上的函数值及一阶导数值,,用埃尔米特(Hermite)插值公式计算出给定插值点处的函数近似值.
程序如下:
#include
#include
#include
#include
#define EPSILON 1.0e-12
#define N 3
double hermite(double x,double xi[N],double yi[N],double dyi[N])
{ int i,j;
static double li,sum,y,gix[N],hix[N];
for(i=0;i
{ li=1.0;sum=0.0;
for(j=0;j
if(j!=i)
{ if(fabs(xi[i]-xi[j])
{ printf("The interpolation base points overlapping!\n");
printf("Strike any key to exit!\n");
getch();
exit(1);
hermite插值c语言程序,张艳-埃尔米特Hermite 插值逼近的C语言程序.doc相关推荐
- 埃尔米特插值c语言编程,131张艳-埃尔米特(Hermite)插值逼近的C语言程序.doc
PAGE PAGE 1 论文题目:埃尔米特(Hermite) 插值逼近的C语言程序 院 系: 数学科学学院 专 业: 数学与应用数学 姓 名: 张 艳 学 号: 指导教师: 侯 存 贵 完成时间: 2 ...
- 插值与拟合 (一) : 拉格朗日多项式插值 、Newton插值 、分段线性插值、Hermite插值 、样条插值、 B 样条函数插值、二维插值
插值:求过已知有限个数据点的近似函数. 拟合:已知有限个数据点,求近似函数,不要求过已知数据点,只要求在某种意义下它在这些点上的总偏差最小. 插值和拟合都是要根据一组数据构造一个函数作为近似,由于近似 ...
- c语言埃尔米特插值思路,【数学建模算法】(26)插值和拟合:埃尔米特(Hermite)插值和样条插值...
1.埃尔米特(Hermite)插值 1.1.Hermite插值多项式 如果对插值函数,不仅要求它在节点处与函数同值,而且要求它与函数有相同的一阶.二阶甚至更高阶的导数值,这就是 Hermite 插值问 ...
- c语言上机指导实验7,[C语言程序设计上机指导与同步训练(刘迎春张艳霞)]实验7.pdf...
您所在位置:网站首页 > 海量文档  > 计算机 > C/C++资料 [C语言程序设计上机指导与同步训练(刘迎春张艳霞 ...
- 埃尔米特(Hermite)插值
Hermite插值满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值.对于高阶导数的情况,Hermite插值多项式比较复杂,在实际情况中,常常遇到的是函数值与一阶导数给定的情况.在此情 ...
- 艾尔米特插值的MATLAB实现,埃尔米特(Hermite)插值
Hermite插值满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值.对于高阶导数的情况,Hermite插值多项式比较复杂,在实际情况中,常常遇到的是函数值与一阶导数给定的情况.在此情 ...
- c语言设计基础张小毅答案,C语言程序设计基础
<c语言程序设计上机指导与习题集>是<c语言程序设计>(全丽莉.明平象主编,重庆大学出版社出版)的配套教材.本书内容分为两部分.第1部分是c语言程序设计环境介绍和c语言程序设计 ...
- 数学建模准备 插值(拉格朗日多项式插值,牛顿多项式插值,分段线性插值,分段三次样条插值,分段三次Hermite插值)
文章目录 摘要(必看) 0 基础概念 什么是插值 插值用途 什么是拟合 插值和拟合的相同点 插值和拟合的不同点 1 常用的基本插值方法 1.1 多项式插值法 1.1.1 拉格朗日多项式插值法 多项式插 ...
- c语言程序与设计教学设计,浅谈C语言程序设计课程的教学设计
C语言程序设计课程是一门计算机专业程序设计语言类的基础课程,如何将枯燥复杂的语法.单调抽象的函数让学生们有兴趣学习.研究.文章通过在C语言程序设计这门课教学的设计,浅谈了自己的见解和方法. 大多数计算 ...
- 【数学与算法】【分段三次Hermite插值】和【分段三次样条插值】
光滑曲线在数学上的定义是什么?? 原文链接:光滑曲线在数学上的定义是什么? 回答1: 定义:切线随切点的移动而连续转动. 若函数f(x)f(x)f(x)在区间(a,b)(a,b)(a,b)内具有一阶连 ...
最新文章
- 为什么catch了异常,但事务还是回滚了?
- 矩阵的三维变换(转)
- cachecloud:Redis云管理平台
- 纵观计算机网络发展历程,人工智能在计算机网络技术中的应用分析
- 著名的用户界面设计准则
- linux指定运行so文件,linux – 强制二进制文件使用特定(较新)版本的共享库(.so)...
- 双物块弹簧阻尼系统模型,李雅普诺夫稳定性分析,matlab仿真验证
- EZEMC测试软件_AR EMC测试软件EMCWARE
- 使用cmd命令行查看Windows系统激活信息
- cmake-自动识别新增子模块
- bind dns mysql,linux下bind9.8+dlz+mysql 的dns服务器局域网配置
- 西恩科技更新招股书:IPO前大手笔分红“套现”, 赵志安为实控人
- Algorithm negotiation fail解决
- 深度学习入门(三十七)计算性能——硬件(TBC)
- 用html语言编写勾股定理,HTML5 勾股定理动画演示
- [4G5G专题-131]:流程 - LTE的功率控制
- 计算机科学速成课笔记(详解)
- oracle19c克隆PDB
- NOI Linux备忘(竞赛)
- parseObject和fromObject的区别
热门文章
- php168 shell,php168 v6/v7鸡肋Getshell安全问题
- xp怎么看计算机是多少位的,WinXP系统怎么看电脑是32位还是64位?
- C# TextBox 自动换行问题解决
- Matlab无法使用GBK编码保存文件,改用UTF-8编码
- Git问题解决方案:不建议在没有为偏离分支指定合并策略时执行pull操作(Pulling without specifying how to reconcile divergent branches)
- 服务器30hz显示器240hz,显示器刷新率上不去,这锅到底让谁背
- 小米格式化fastboot_小米fastboot刷机教程
- CPU的DataSheet和UserManual的区别
- 小游戏学习--获取已发布微信小游戏源码 creator cocos 编译出微信平台代码 喵喵切水果 跑酷
- python生成excel文件二维码_Python实现读取Excel表内容批量生成二维码