数值计算-插值-单项式基底-牛顿基底-拉格朗日基底算法实例.docVIP

数值计算实验报告 (1)用单项式基函数对上述数据进行插值 (2)对上述数据进行Lagrange插值 (3) 对上述数据进行Newton插值 并画出对应的插值曲线。 实验目的 多项式插值 编程环境 采用C++或JAVA实现 实验内容（算法、程序、步骤和方法） #include stdafx.h #include iostream #include iomanip using namespace std; class Matrix { public: Matrix(); Matrix(int m,int n); Matrix(Matrix another); void display(); void clear(); int lu(); double **mat; int row; int col; }; Matrix::Matrix(){ row=0; col=0; } Matrix::Matrix(int m,int n){ row=m; col=n; mat=new double*[row]; for(int i=0;irow;i++){ mat[i]=new double[row]; for(int j=0;jcol;j++) mat[i][j]=0; } } void Matrix::display(){ for(int i=0;irow;i++) { for(int j=0;jcol;j++) coutsetprecision(17)mat[i][j] ; coutendl; } } void Matrix::clear(){ for(int i=0;irow;i++) for(int j=0;jcol;j++) mat[i][j]=0; } int Matrix::lu(){ double b[4]={11,29,65,125}; int col=0,row=0,p[4];//p为a排?列￠D阵¨?压1缩?形?式o? double x[4]={},u[4][4]={},l[4][4]={},temp,y[4]={}; for(int i=0;i4;i++){ l[i][i]=1; p[i]=i; for(int j=0;j4;j++) u[i][j]=this-mat[i][j]; }//初?始o?化?￥三¨y角?阵¨? for(col=0;col4;col++){ p[col]=col; for(row=col+1;row4;row++){ if(fabs(u[row][col])fabs(u[(p[col])][col])) p[col]=row; } if(p[col]!=col) for(int k=0;k4;k++){ temp=u[col][k]; u[col][k]=u[p[col]][k]; u[p[col]][k]=temp; }//选?主??元a if(u[col][col]==0) break; for(row=col+1;row4;row++){ l[row][col]=u[row][col]/u[col][col]; }//某3一°?列￠D每?一°?行D的ì?乘?数oy for(int j=col;j4;j++) for(int i=col+1;i4;i++){ u[i][j]=u[i][j]-l[i][col]*u[col][j]; } }//U for(int k=1;k4;k++){ for(col=0;colk;col++){ temp=l[p[k]][col]; l[p[k]][col]=l[k][col]; l[k][col]=temp; } }//L for(int i=0;i4;i++){ temp=b[i]; b[i]=b[p[i]]; b[p[i]]=temp; }//Pb for(col=0;col4;col++){ if (l[col][col]==0){ break; } y[col]=b[col]/l[col][col]; for(row=col+1;row4;row++){ b[row]=b[row]-l[row][col]*y[col]; } }//Ly=Pb; for(col=4-1;col=0;col--){ if(u[col][col]==0){ break; } x[col]=b[col]/