实验五：分段线性插值.docVIP

《数值分析》实验报告 实验序号：实验五 题目名称: 分段线性插值? 学号: 姓名: 任课教师: 马季骕 专业班级:计算机科学与技术（非师范） 1、 实验目的： 设在区间[a,b]上，给定n+1个插值节点a=X0X1X2。。。Xn=b和相应的函数值y0，y1，y2，。。。，yn ，求一个插值函数φ（x），具有下面性质： （1）φ（x）=yj（j=0,1,2,…,n）? （2）φ（x）在每一个小区间[xj，yj+1]上是线性函数。 2、? 算法分析： 分段线性插值的算法思想： 分段线性插值需要在每个插值节点上构造分段线性插值基函数 ，然后再作它们的线性组合。分段线性插值基函数的特点是在对应的插值节点上函数值取 1，其它节点上函数值取0。插值基函数如下： 有了基函数就可以直接写出分段线性插值函数的表达式 。 具体程序设计： for(i=0;i=20;i++)?????? //选取节点 { ?????? ax[i]=-1+((2/20.0)*i);? //选取 上的21个对称的节点 ????? ?ay[i]=1.0/(1+25*ax[i]*ax[i]); } ?x1=-1; while(x1=1) { ????? m=0; ????? x1=x1+0.00001; ????? for(i=0;i=20;i++) ????? { ????????? if(i==0x1=-1x1=-0.9) ????????? { ??????????? n=ay[0]*((x1-ax[1])/(ax[0]-ax[1])); //计算 并和相应的函数值组合 ??????????? m=n+m; ?????????? } ????????? else ????????? { ?? ??????????if(x1=ax[i-1]x1=ax[i]) //计算 并和相应的函数值组合 ???????????? { ??????????????? n=ay[i]*((x1-ax[i-1])/(ax[i]-ax[i-1])); ??????????????? m=n+m; ???????????? } ???????????? else ???????????? { ??????????????? if(x1ax[i]x1=ax[i+1]) ??????????????? { ?????????????????? n=ay[i]*((x1-ax[i+1])/(ax[i]-ax[i+1])); ?????????????????? m=m+n; ???????????????? } ??????????????? else ??????????????? { ?????????????????? if(i==19x1ax[19]x1=ax[20]) ????????????????? { ??????????????????? n=ay[20]*((x1-ax[19])-(ax[20]-ax[19]));//计算 并和相应的函数值组合 ??????????????????? m=n+m; ?????????????????? } ???????????????? } ????????????? } ?????????? } ?????? } ?? } 3、?实验结果截图： 在[-1,1]区间上选取了21个等分节点的分段线性插值函数的图像如下： ?? 4、? 程序代码 // SHIYAN456View.cpp : implementation of the CSHIYAN456View class // ? #include stdafx.h #include SHIYAN456.h ? #include SHIYAN456Doc.h #include SHIYAN456View.h ? #ifdef _DEBUG #define new DEBUG_NEW #undef THIS_FILE static char THIS_FILE[] = __FILE__; #endif ? ///////////////////////////////////////////////////////////////////////////// // CSHIYAN456View ? IMPLEMENT_DYNCREATE(CSHIYAN456View, CView) ? BEGIN_MESSAGE_MAP(CSHIYAN456View, CView) ????? //{{AFX_MSG_MAP(CSHIYAN456View) ????? ON_COMMAND(ID_FFunction, OnFFunction) ????? ON_COM