數值計算_插值_單項式基底_牛頓基底_拉格朗日基底算法實例

《數值計算_插值_單項式基底_牛頓基底_拉格朗日基底算法實例》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。

1、數值計算實驗報告學號：201000301220姓名：魏源班級：泰山學堂計算機日期:20120512實驗題目(1)用單項式基函數對上述數據進行插值(2)對上述數據進行Lagrange插值(3)對上述數據進行Newton插值并畫出對應的插值曲線。實驗目的多項式插值編程環(huán)境采用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縮?形?式o?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];}//初?始o?化?￥三¨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的ì?乘?數oyfor(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