三次埃尔米特插值.docVIP

《计算方法》课程设计报告 学生姓名: 学 号： 陈 洋 1009300109 刘 睿 1009300122 学 院: 班 级: 题 目: 指导教师： 职称 讲师 朱秀丽 讲师 尚宝欣 讲师 2012年月日目 录目 录 I 一、 摘要 1 二、 算法设计 2 2.1分段线性插值 2 2.2 分段三次埃尔米特插值 2 2.3 功能框图 3 三、 例题计算 4 四、误差及结果分析 9 4.1 例题误差分析 9 4.2 结点个数对插值结果的影响 9 五、总结及心得体会 12 参考文献 13 源程序 14 摘要 分段线性插值与分段定义的线性插值,在相邻插值节点的区间上对应的是同一个线性函数由于它们的表现形式不一样从而产生为两种不同的计算方法,相应的误差表现形式也不一样.拉格朗日插值余项利用f(x)的二阶导数,要f(x)的二阶导数存在,对于二阶导数不存在的情况不能估算出它的误差,所以适用范围比较小.现在我们可以利用一阶导数就估算出误差,给计算带来许多的方便为了避免高次插值可能出现的大幅度波动现象，在实际应用中通常采用分段低次插值来提高近似程度，比如可用分段线性插值或分段三次埃尔米特插值来逼近已知函数，但它们的总体光滑性较差。为了克服这一缺点，一种全局化的分段插值方法——成为比较理想的工具。 二、算法设计 2.1 分段线性插值 分段线性插值就是通过插值点用折线段连接起来的逼近函数。设已知节点上的函数值记求一折线函数满足： （1） （2） （3）在每个小区间上是线性函数. 则称为分段线性插值函数。 由定义可知在每个小区间上可表示为 分段线性插值的误差估计可利用插值余项公式 得到 或 其中由此还可得到 在上一致成立，故在上一致收敛到。 2.2 分段三次埃尔米特插值 分段线性插值函数的导数是间断的，若在节点上除已知函数值外还给出导数值这样就可构造一个导数连续的分段插值函数它满足条件： （1） （2）; （3）在每个小区间上是三次多项式。 根据两点三次插值多项式可知，在区间上的表达式为 上式对于成立。 利用三次埃尔米特插值多项式的余项可得误差估计 其中 2.3功能框图 程序可以通过输入插值函数，插值结点个数，来实现插值，还可以选择插值方式，具体框图如下 三、例题计算 插值函数为 输入6个节点的分段线性插值 输入6个节点的分段三次埃尔米特插值 请输入插值函数：x^2+1 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 1 请输入插值节点数N=6 是否继续（1为继续）：1 请输入插值函数：x^2 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 2 请输入插值节点数N=6 是否继续（1为继续）：2 2.插值函数为 输入节点为8的分段线性插值 输入节点为8的分段三次埃尔米特插值 请输入插值函数：cos(x) 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 1 请输入插值节点数N=8 是否继续（1为继续）：1 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 2 请输入插值节点数N=8 3.插值函数为 输入节点为10的分段线性插值 输入节点为10的分段线三次埃尔米特插值 请输入插值函数：exp(x) 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 1 请输入插值节点数N=10 是否继续（1为继续）：1 请选择插值法（1为分段线性插值，2为分段三次埃尔米特插值）： 2 请输入插值节点数N=10 四、结果分析 误差 n值取8,由分段线性插值余项公式 可得误差限为 由插值图像和误差可知误差结果较大 由分段三次埃尔米特插值余项公式 可得误差限为 误差为零因此分段三次埃尔米特插值较为精确 2. 误差 n值取10，分段线性插值误差限为 分段三次埃尔米特插值误差限为 由误差结果可知分段三次埃尔米特插值得到的结果更为精确 3. 误差 n值取12，分段线性插值误差限为 分段三次埃尔米特插值误差限为 由误差结果可知分段三次埃尔米特插值得到的结果更为精确 综上可知，同等条件下分段三次埃尔米特插值得到的结果更为精确 4.2结点个数对插值结果的影响 为了进一步研究误差与插值节点的关系，我们以例题3的分段三次埃尔米特插值为例 n值取4时，误差为15.0974，插值图像如下 n值取5时，误差为6.1839，插值图像如下 n值取6时，误差为2.9822，插值图像如下 n值取8时，误差为0.9436，插值图像如下 通过误差限和得到的