三次埃尔米特插值.docVIP

《计算方法》课程设计报告 学生姓名: 张学阳 学 号： 1009300132 陈 洋 1009300109 刘 睿 1009300122 学 院: 理学院 班 级: 数学101 题 目: 分段线性及三次埃尔米特插值通用程序 指导教师： 宋云飞 职称： 讲师 朱秀丽 讲师 尚宝欣 讲师 2012年 目 录 TOC \o 1-1 \h \z \u HYPERLINK \l _Toc344060456 目 录 PAGEREF _Toc344060456 \h I HYPERLINK \l _Toc344060457 一、 摘要 PAGEREF _Toc344060457 \h 1 HYPERLINK \l _Toc344060457 二、 算法设计 PAGEREF _Toc344060457 \h 2 HYPERLINK \l _Toc344060457 2.1分段线性插值 PAGEREF _Toc344060457 \h 2 HYPERLINK \l _Toc344060457 2.2 分段三次埃尔米特插值 PAGEREF _Toc344060457 \h 2 HYPERLINK \l _Toc344060457 2.3 功能框图 PAGEREF _Toc344060457 \h 3 HYPERLINK \l _Toc344060457 三、 例题计算 PAGEREF _Toc344060457 \h 4 HYPERLINK \l _Toc344060458 四、误差及结果分析 PAGEREF _Toc344060458 \h 9 HYPERLINK \l _Toc344060457 4.1 例题误差分析 PAGEREF _Toc344060457 \h 9 HYPERLINK \l _Toc344060457 4.2 结点个数对插值结果的影响 PAGEREF _Toc344060457 \h 9 HYPERLINK \l _Toc344060459 五、总结及心得体会 PAGEREF _Toc344060459 \h 12 HYPERLINK \l _Toc344060460 参考文献 PAGEREF _Toc344060460 \h 13 HYPERLINK \l _Toc344060461 源程序 PAGEREF _Toc344060461 \h 14 摘要 分段线性插值与分段定义的线性插值,在相邻插值节点的区间上对应的是同一个线性函数。由于它们的表现形式不一样从而产生为两种不同的计算方法,相应的误差表现形式也不一样.拉格朗日插值余项利用f(x)的二阶导数,要f(x)的二阶导数存在,对于二阶导数不存在的情况不能估算出它的误差,所以适用范围比较小.现在我们可以利用一阶导数就估算出误差,给计算带来许多的方便。 为了避免高次插值可能出现的大幅度波动现象，在实际应用中通常采用分段低次插值来提高近似程度，比如可用分段线性插值或分段三次埃尔米特插值来逼近已知函数，但它们的总体光滑性较差。为了克服这一缺点，一种全局化的分段插值方法——三次样条插值成为比较理想的工具。 在代数插值过程中，人们为了获得较好的近似效果，通常情况下是增加插值节点数.由于二次插值比线性插值近似效果好，因此容易错误地认为插值多项式次数越高越好.事实上,随着插值节点的增多，插值多项式不一定收敛到被插值函数.。通过分段低次插值或样条插值可以得到较好的近似逼近函数，分段低次插值具有公式简单、运算量小、稳定性好、收敛性有保证等优点.随着子区间长度h取得足够小，分段低次插值总能满足所要求的精度.因此分段低次插值应用十分广泛.。分段线性插值是分段低次插值中常见的方法之一，在本文中对函数在(-5，5)上进行分段线性插值，取不同节点个数n，得到不同分段线性插值函数.并用MATLAB编写分段线性插值函数，最后比较用不同节点数所得插值函数与真实函数的误差，从而得出节点数与插值效果的关系。 二、算法设计 2.1 分段线性插值 分段线性插值就是通过插值点用折线段连接起来的逼近函数。设已知节点上的函数值记求一折线函数满足： （1） （2） （3）在每个小区间上是线性函数. 则称为分段线性插值函数。 由定义可知在每个小区间上可表示为 分段线性插值的误差估计可利用插值余项公式 得到 或 其中由此还可得到 在上一致成立，故在上一致收敛到。 2.2 分段三次埃尔米特插值 分段线性插值函数的导数是间断的，若在节点上除已知函数值外还给出