[2018年最新整理]实验一插值计算-1240111118-武元甲.docVIP

实验一插值方法 一 实验目的 通过本次上机实习，能够进一步加深对各种插值算法的理解；学会使用用三种类型的插值函数的数学模型、基本算法，结合相应软件（如VC/VB/Delphi/Matlab/JAVA/Turbo C）编程实现数值方法的求解。并用该软件的绘图功能来显示插值函数，使其计算结果更加直观和形象化。 二 实验内容 通过程序求出插值函数的表达式是比较麻烦的，常用的方法是描出插值曲线上尽量密集的有限个采样点，并用这有限个采样点的连线，即折线，近似插值曲线。取点越密集，所得折线就越逼近理论上的插值曲线。本实验中将所取的点的横坐标存放于动态数组中，通过插值方法计算得到的对应纵坐标存放于动态数组中。 本实验将Lagrange插值、Newton插值和三次样条插值实现为一个C++类CInterpolation，并在Button单击事件中调用该类相应函数，得出插值结果并画出图像。CInterpolation类为 class CInterpolation { public: CInterpolation();//构造函数 CInterpolation(float *x1, float *y1, int n1);//结点横坐标、纵坐标、下标上限 ~ CInterpolation();//析构函数 ………… ………… int n, N;//结点下标上限，采样点下标上限 float *x, *y, *X;//分别存放结点横坐标、结点纵坐标、采样点横坐标 float *p_H,*p_Alpha,*p_Beta,*p_a,*p_b,*p_c,*p_d,*p_m;//样条插值用到的公有指针，分别存放，，，，，，和 }; 其中，有参数的构造函数为 CInterpolation(float *x1, float *y1, int n1) { //动态数组x1，y1中存放结点的横、纵坐标，n1是结点下标上限（即n1+1个结点） n=n1; N=x1[n]-x1[0]; X=new float[N+1]; x=new float[n+1]; y=new float[n+1]; for(int i=0;i=n;i++) { x[i]=x1[i]; y[i]=y1[i]; } for(int i=0;i=N;i++) X[i]=x[0]+i; } 2.1 Lagrange插值 ，其中 对于一个自变量，要求插值函数值，首先需要计算对应的Lagrange插值基函数值 float l(float xv,int i) //求插值基函数的值 { float t=1; for(int j=0;j=n;j++) if(j!=i) t=t*(xv-x[j])/(x[i]-x[j]); return t; } 调用函数l(float x,int i)，可求出 float p_l(float x) //求在一个点的插值结果 { float t=0; for(int i=0;i=n;i++) t+=y[i]*l(x,i); return t; } 调用p_l(float x)可实现整个区间的插值 float *Lagrange() //求整个插值区间上所有采样点的插值结果 { float *Y=new float[N+1]; for(int k=0;k=N;k++) Y[k]=p_l(x[0]+k*h); return Y; } Newton插值 ， 其中， 对于一个自变量，要求插值函数值，首先需要计算出和 float *f() { //该函数的返回值是一个长度为n＋1的动态数组，存放各阶差商 } float w(float x, int i) { //该函数计算 } 在求的函数中调用*f()得到各阶差商，然后在循环中调用w(float x)可得出插值结果 float p_n(float x) { //该函数计算在一点的值 } 调用p_n(float x)可实现整个区间的插值 float *Newton() { //该函数计算出插值区间内所有点的值 } 2.3 三次样条插值 三次样条插值程序可分为以下四步编写： 计算结点间的步长、、； 利用、、产生三对角方程组的系数矩阵和常数向量； 通过求解三对角方程组，得出中间结点的导数； 对自变量x，在对应区间上，使用Hermite插值； 调用上述函数，实现样条插值。 将每步写成函数： (1) void GetH(void) { //该函数