实验1拉格朗日插值与牛顿插值

《实验1拉格朗日插值与牛顿插值》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西华大学计算机系学生上机实践报告西华数学与计算机学院上机实践报告课程名称：计算方法A年级：2010级上机实践成绩：指导教师：严常龙姓名：上机实践名称：拉格朗日插值和牛顿插值法学号：上机实践日期：上机实践编号：1上机实践时间：一、目的1．通过本实验加深对拉格朗日插值和牛顿插值法构造过程的理解；2．能对上述两种插值法提出正确的算法描述编程实现。二、内容与设计思想自选插值问题，编制一个程序，分别用拉格朗日插值法和牛顿插值法求解某点的函数近似值。（从课件或教材习题中选题）已知y=f(x)的数据表如下，求t=0.63处的函数值z=f(t)。ixiyiixiyi10.100.90483760.570.56

2、552520.150.86070870.700.49658530.250.77880180.850.42741540.400.67032090.930.39455450.500.606531101.000.367879三、使用环境操作系统：软件环境：四、核心代码及调试过程1.拉格朗日插值代码#includemain(){floatx[8],y[8];floatt,z=0.0,s;intn=0,i=0,tp=0,start,end,j;printf("请输入节点控制数n");printf("");scanf("%d",&n);printf("请输入节点t的值");print

3、f("");第5页共5页西华大学计算机系学生上机实践报告scanf("%f",&t);//节点赋值for(i=0;in-1){end=n-1;start=end-7;}//用所选插值点做拉格朗日插

4、值多项式运算for(i=start;i#definemax20main(){doublex[max],y[max];//插值点xi以及对应函数yi的函数值doublet,fw;//t为待求的参数点fw为对应要求的值doubleft[max

5、][max];//表示各阶差商第二维为阶数intn=0,i=0,j=0;doubleb[max];/*printf("请输入节点控制数n");printf("");scanf("%d",&n);printf("请输入节点t的值");printf("");scanf("%f",&t);//节点赋值for(i=0;i

6、;y[1]=0.860708;x[2]=0.25;y[2]=0.778801;x[3]=0.40;y[3]=0.670320;x[4]=0.50;y[4]=0.606531;x[5]=0.57;y[5]=0.565525;x[6]=0.70;y[6]=0.496585;x[7]=0.85;y[7]=0.427415;第5页共5页西华大学计算机系学生上机实践报告x[8]=0.93;y[8]=0.394554;x[9]=1.00;y[9]=0.367879;for(i=0;i

7、t[i][0]=y[i];//计算0阶差商}for(i=1;i=1;i--){b[i-1]=ft[i-1][i-1]+b[i]*(t-x[i-1]);}fw=b