数学实验“等距节点插值Hermite插值分段插值”实验报告.docx

西京学院数学软件实验任务书 课程名称学号 数学软件实验 班级 0912020107 姓名 数 0901 李亚强 实验课题 实验目的 实验要求 实验内容 等距节点插值，Hermite 插值，分段插值（线性，二次，三次） 熟悉等距节点插值，Hermite 插值，分段插值（线性， 二次，三次） 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成 等距节点插值，Hermite 插值，分段插值（线性，二次，三次） 成绩 教师 - - PAGE 1 - 实验十六实验报告 一、实验名称：等距节点插值，Hermite 插值，分段插值（线性， 二次，三次）。 二、实验目的：进一步熟悉等距节点插值，Hermite 插值，分段插值（线性，二次，三次）。 三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 四、实验原理： 等距节点插值： 差分分为前向差分、后向差分和中心差分三种，它们的记法及 定义如下所示: n 阶前向差分公式?n f (x ) ? ?n?1 f (x ) ? ?n?1 f (x ) i i ?1 i n 阶后向差分公式?n f (x ) ? ?n?1 f (x ) ? ?n?1 f (x ) i i i?1 11n 阶中心差分公式? n f (x ) ? ? n?1 f (x ) ?? n?1 f (x ) 1 1 i i? i? 2 2 其中： ? -前向差分； ? -后向差分；? -中心差分。 假设?0 f (x ) ? ?0 f (x ) ? ? 0 f (x ) ? f (x ) ，为了方便计算，构造 差分表( f i i i i i ? f (x ) )。 i 这里只说明前向牛顿插值，其多项式可表示为如下形式： N (x) ? N (x0 ? th) ? t ? ? t ? ? t ? ? f (x ) ? ? ? ?f (x ) ? ? ? ?2 f (x ) ? ? ? ? ?n f (x ) 0 ?1? 0 ? 2 ? 0 ? n ??0 其中h 为步长， h ? x 1 埃尔米特插值： x ，且的取值范围为0 ? t ? n 。 0 埃尔米特插值法满足在节点上等于给定函数值，而且在节点上的导数值也等于给定的导数值，对于有高阶导数的情况，埃尔米特插值多项式比较复杂，在实际应用中，常常遇到的是函数值与一阶导数值给定的情况，在这种情况下，n 个节 点 x , x , x 的埃尔米特插值多项式H (x) 的表达形式如下所示： 1 2 n H (x) ? ?n h [(x ? x)(2a y ? y ) ? y ] i i i?1 i i i i 其中 y ? y(x ), y ? y (x ) i i i h ? ?n i x ? x ( j  )2 , a ? ?n 1 i 分段插值： j ?1 j ?i x ? x i j i j ?1 j ?i x ? x i j 给定插值节点 x i 、 节点函数值 y i 及对应的导数值 y (i ? 0,1,2, , N ) ，则满足下面条件 i p(x ) ? y , p (x ) ? y i i i i 的分段埃尔米特插值函数 p(x) 的表达式如下所示： p (x) ? y (1? 2 x ? x i )( x ? x i?1 )2 ? y (1? 2 x ? x i?1 )( x ? x i )2 ? i i h h i i x ? x i?1 ?h h i i x ? x y (x ? x )( i?1 )2 ? y (x ? x )( i )2 i i h i i?1 i?1 h i h ? x ? x ,( i ? 0,1,2, N ?1), x ?[x , x ] i i?1 i i i ?1 五、实验内容： %等距节点插值 function [f,f0]= dengjujiedian(x,y,x0) syms t; if(length(x) == length(y)) n = length(x); c(1:n) = 0.0; else disp(x和y的维数不相等！); return; end f = y(1); y1 = 0; xx =linspace(x(1),x(n),(x(2)-x(1))); if(xx ~= x) disp(节点之间不是等距的！); return; end for(i=1:n-1) for(j=1:n-i) y1(j) = y(j+1)-y(j); end c(i) = y1(1); l = t; for(k=1:i-1) l = l*