数值计算实验5.docVIP

本科实验报告 课程名称： 计算机数值方法B 实验项目： 代数插值 实验地点： 逸夫楼302 专业班级： 学号： 学生姓名： 指导教师： 2013年 4月 20日 一、实验目的和要求 学习使用拉格朗日插值法或牛顿插值法求解 二、实验内容和原理 使用拉格朗日插值法或牛顿插值法求解：已知f(x)在6个点的函数值如下表所示，运用插值方法，求f(0.596)的近似值。 X 0.40 0.55 0.65 0.80 0.90 1.05 f(x) 0.41075 0.57815 0.69675 0.88811 1.02652 1.25386 原理： 设函数在区间[a,b]上n+1互异节点x0,x1,…,xn上的函数值分别为y0,y1,…,yn,求n次插值多项式Pn(x),满足条件 Pn(xj)=yj, j=0,1,…,n 令 Ln(x)=y0l0(x)+y1l1(x)+…+ynln(x)= ∑yili(x) 其中l0(x),l1(x),…, ln(x) 为以x0,x1,…,xn为节点的n次插值基函数，则Ln(x)是一次数不超过n的多项式，且满足 Ln(xj)=yj, L=0,1,…,n 再由插值多项式的唯一性，得 Pn(x)≡Ln(x) 三、主要仪器设备 惠普ProBook6460b VC++运行环境 四、操作方法与实验步骤 #include stdio.h #include math.h #include stdafx.h #includeiostream using namespace std; typedef struct data { float x; float y; }Data; //变量x和函数值y的结构 Data d[20]; //最多二十组数据 float f(int s,int t) //牛顿插值法，用以返回插商 { if(t==s+1) return (d[t].y-d[s].y)/(d[t].x-d[s].x); else return (f(s+1,t)-f(s,t-1))/(d[t].x-d[s].x); } float Newton(float x,int count) { int n; while(1) { cout请输入n值(即n次插值):;//获得插值次数 cinn; if(n=count-1)// 插值次数不得大于count－次 break; else system(cls); } float t=1.0; float y=d[0].y; float yt=0.0; for(int j=1;j=n;j++) { t=(x-d[j-1].x)*t; yt=f(0,j)*t; y=y+yt; } return y; } float lagrange(float x,int count) { float y=0.0; for(int k=0;kcount;k++)//这儿默认为count－次插值 { float p=1.0;//初始化p for(int j=0;jcount;j++) {//计算p的值 if(k==j)continue;//判定是否为同一个数 p=p*(x-d[j].x)/(d[k].x-d[j].x); } y=y+p*d[k].y;//求和 } return y;//返回y的值 } void main() { float x,y; int count; while(1) { cout请输入x[i],y[i]的组数，不得超过组:;//要求用户输入数据组数 cincount; if(count=20) break;//检查输入的是否合法 system(cls); } //获得各组数据 for(int i=0;icount;i++) { cout请输入第i+1组x的值:; cind[i].x; cout请输入第i+1组y的值:; cind[i].y; system(cls); } cout请输入x的值:;//获得变