数值分析实验五(分段线性插值).docVIP

《数值分析》实验报告 实验编号：实验五 课题名称：分段线性插值 一、算法介绍 给定的函数为f(x)=1/(1+25x*x)，-1=x=1，将给定区间分成10分，得到11个节点：x[0],x[1],…,x[10]，对于[x[0],x[1]]区间上，插值基函数L[0]=(x-x[1])/(x[0]-x[1])，其他区间上为零。对于[x[i-1],x[i]]区间上，插值基函数L[i] = (x-x[i-1])/(x[i]-x[j-1])。对于(x[i],x[i+1]]上，插值基函数L[i] = (x-x[i+1])/(x[i]-x[i+1])，其中i=1，2，…，10。 对于(x[n-1],x[n]]区间上，L[i] = (x-x[10])/(x[9]-x[10]),其他区间上为零。构造得到的多项式g(x)=∑f(x[i])*l[i],其中i=0,1,2,…,10。 二、程序代码 // testView.cpp : implementation of the CTestView class // #include stdafx.h #include test.h #include testDoc.h #include testView.h #ifdef _DEBUG #define new DEBUG_NEW #undef THIS_FILE static char THIS_FILE[] = __FILE__; #endif ///////////////////////////////////////////////////////////////////////////// // CTestView IMPLEMENT_DYNCREATE(CTestView, CView) BEGIN_MESSAGE_MAP(CTestView, CView) //{{AFX_MSG_MAP(CTestView) // NOTE - the ClassWizard will add and remove mapping macros here. // DO NOT EDIT what you see in these blocks of generated code! //}}AFX_MSG_MAP // Standard printing commands ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint) ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint) ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview) END_MESSAGE_MAP() ///////////////////////////////////////////////////////////////////////////// // CTestView construction/destruction CTestView::CTestView() { // TODO: add construction code here } CTestView::~CTestView() { } BOOL CTestView::PreCreateWindow(CREATESTRUCT cs) { // TODO: Modify the Window class or styles here by modifying // the CREATESTRUCT cs return CView::PreCreateWindow(cs); } ///////////////////////////////////////////////////////////////////////////// // CTestView drawing void CTestView::OnDraw(CDC* pDC) { CTestDoc* pDoc = GetDocument(); ASSERT_VALID(pDoc); // TODO: add draw code for native data here int i,j,k; double x,y,p_x,p_y,l,xx[100],f[100],F[100],sum,p_sum; CPen MyPen,*OldPen; pDC-SetViewportOrg(400,400); //定义坐标原点 for(i=-500;i500;i++) { pDC-SetPixel(0,i,RGB(0,0,0)); pDC-S