实验二Newton插值多项式(4页).doc

《实验二Newton插值多项式(4页).doc》由会员分享，可在线阅读，更多相关《实验二Newton插值多项式(4页).doc（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-实验二 Newton差值多项式1. 实习目的（1） 巩固Newton差值多项式的算法，进一步理解Newton差值问题。（2） 熟练掌握C语言程序设计相关知识。2. 班级：计算092，姓名：薛藏朋，学号30908110723.目的意义4.数学建模（1）差商fx,x,x=；（2） 5.算法Step1:输入n值及（Step2:for i=0,1,n 计算差商fx,x,x=,其中Step3:对于给定的x由 Step4：计算出的值.Step5:输出。6.（1）流程图开始 输入n,xi,yi,xtmp=1,newton=0pointj=(pointj-pointj-1)/(pointsj.x-point

2、sj-1-i.x)i=n?N Ytmp=tmp*(x-pointsi.x)； newton=newton+tmp*pointi+1 i=n? N Y输出Nowton结束 (2)程序 #include #define N 100struct POINT_type double x; double y;void main() int i,j,n; struct POINT_type pointsN; double pointN; double x,tmp=1,newton=0; printf(please input n,0=n=100: n); scanf(%d,&n); printf(Now

3、input (xi,yi),%d: n,n); for(i=0;in;i+) scanf(%lf,%lf,&pointsi.x,&pointsi.y); printf(Now input x: ); scanf(%lf,&x); for(i=0;in;i+) pointi=pointsi.y; for(i=0;ii;j-) pointj=(pointj-pointj-1)/(pointsj.x-pointsj-1-i.x); newton=point0; for(i=0;in;i+) for(j=0,tmp=1;ji;j+) tmp=tmp*(x-pointsi.x); newton=newton+tmp*pointi+1; printf(newton(%lf)=%lfn,x,newton); 7.数值算例8.对计算结果进行分析和评价 结果不是很准确，有待改进9.参考文献【1】张毅坤，曹锰，张亚玲，C语言程序设计教程，西安交通大学出版社，2003.【2】姚全珠，李薇，王晓帆，C+面向对象程序设计，北京：电子工业出版社，2010.【3】秦新强，数值逼近，西安理工大学【4】王萼芳，石声明，高等代数，北京：高等教育出版社。-第 4 页-