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

《数值计算_插值_单项式基底_牛顿基底_拉格朗日基底算法实例》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数值计算实验报告学号：201000301220姓名：魏源班级：泰山学堂计算机日期:20120512实验题目(1)用单项式基函数对上述数据进行插值(2)对上述数据进行Lagrange插值(3)对上述数据进行Newton插值并画出对应的插值曲线。实验目的多项式插值编程环境采用C++或JAVA实现8实验内容（算法、程序、步骤和方法）#include"stdafx.h"#include"iostream"#include"iomanip"usingnamespacestd;classMatrix{pub

2、lic:Matrix();Matrix(intm,intn);Matrix(Matrix&another);voiddisplay();voidclear();intlu();double**mat;introw;intcol;};Matrix::Matrix(){row=0;col=0;}Matrix::Matrix(intm,intn){row=m;col=n;mat=newdouble*[row];for(inti=0;i

3、r(intj=0;j

4、]={11,29,65,125};intcol=0,row=0,p[4];//p为a排?列¢D阵¨®压1缩?形?式º?doublex[4]={},u[4][4]={},l[4][4]={},temp,y[4]={};for(inti=0;i<4;i++){l[i][i]=1;p[i]=i;for(intj=0;j<4;j++)u[i][j]=this->mat[i][j];}//初?始º?化¡¥三¨y角?阵¨®for(col=0;col<4;col++){p[col]=col;for(row=c

5、ol+1;row<4;row++){if(fabs(u[row][col])>fabs(u[(p[col])][col]))p[col]=row;}if(p[col]!=col)for(intk=0;k<4;k++){temp=u[col][k];u[col][k]=u[p[col]][k];u[p[col]][k]=temp;}//选?主¡Â元aif(u[col][col]==0)break;for(row=col+1;row<4;row++){l[row][col]=u[row][col]/

6、u[col][col];}//某3一°?列¢D每?一°?行D的Ì?乘?数ºyfor(intj=col;j<4;j++)8for(inti=col+1;i<4;i++){u[i][j]=u[i][j]-l[i][col]*u[col][j];}}//Ufor(intk=1;k<4;k++){for(col=0;col

7、emp=b[i];b[i]=b[p[i]];b[p[i]]=temp;}//Pbfor(col=0;col<4;col++){if(l[col][col]==0){break;}y[col]=b[col]/l[col][col];for(row=col+1;row<4;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]

8、/u[col][col];for(row=0;row