Hermite插值报告..docVIP

数值逼近实习报告 题目： 分段三次Hermite插值 班级： 计算111班 学号： 3110811009 姓名： 刘艳平 指导老师： 秦新强 2013-3 题目：分段三次Hermite插值 姓名：刘艳平 学号：3110811009 目的意义 紧接上节的分段线性插值，我们研究分段埃尔米特（Hermite）插值。分段线性插值在特定的每个小区间上只能保证连续性，而不光滑，我们想要得到插值区间上光滑的分段插值多项式，可采用分段埃尔米特（Hermite）插值，再此处我们考虑在整个[a,b]区间上用分段三次Hermite插值多项式来逼近f(x),一般的将带有导数的插值多项式称作Hermite插值多项式。 这样，我们增加了节点处的导数值，这样得到的插值多项式既有很好的准确度，又可以满足一定光滑度的要求。 数学公式 其中，在区间[]上的线性函数为 三、代码 #includeiostream //3次Hermite插值 using namespace std; #define N 2 void main() { double x[N],y[N],yy[N]; int i,j; double X,p1,p2,p3,p4,p,h; cout输入已知点的x坐标:endl; for(i=0;iN;i++) cinx[i]; cout输入已知点的y坐标:endl; for(j=0;jN;j++) ciny[j]; cout输入已知点的导数:endl; for(i=0;iN;i++) cinyy[i]; cout输入X:endl; cinX; for(i=1;i=N;i++) { h=x[i]-x[i-1]; if(Xx[i-1]) { p1=((h+(X-x[i-1])*2)*(X-x[i])*(X-x[i])*y[i-1])/(h*h*h); p2=((h-(X-x[i])*2)*(X-x[i-1])*(X-x[i-1])*y[i])/(h*h*h); p3=((X-x[i-1])*(X-x[i])*(X-x[i])*yy[i-1])/(h*h); p4=((X-x[i])*(X-x[i-1])*(X-x[i-1])*yy[i])/(h*h); p=p1+p2+p3+p4; } } cout结果为:pendl; } 四、数值实例 已知点的横坐标为1、2，纵坐标为2、3，在两点的导数值分别为1和-1，求x=1.5处的y值。 解： 计算结果为2.75， 且运行结果如下： 输入已知点的x坐标: 1 2 输入已知点的y坐标: 2 3 输入已知点的导数: 1 -1 输入X: 1.5 结果为:2.75 Press any key to continue 对计算机结果进行分析评价 分段三次Hermite插值多项式计算的结果比较逼近真实值，而且满足了一定的光滑度。 六、参考文献 [1] 秦新强，数值逼近，西安理工大学出版社，2010  4