- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验4插值法与数值积分
实验4 插值法与数值积分思考问题:插值多项式是否阶次越高越好?数值积分与插值的关系是什么?逐次分半梯形公式求积如何判断误差是否满足要求?答:插值多项式不一定阶次越高越好;数值积分与插值的关系是:二者都是用已知给出的点通过构造新的函数来估计预测未知点,求取我们想知道的内容。数值积分与插值有密切联系,可以先对求积节点做插值,再对插值函数求积分,如果插值函数与被积函数足够接近,就可以把插值函数的积分近似为被积函数的积分;逐次分半梯形公式求积可以用相邻两次求得的结果的差的绝对值来间接判断误差是否满足要求。4.1 用拉格朗日插值法求2的平方根。提示:可以用抛物线插值,,,。运行结果为:4.2 用牛顿插值法求2的平方根。提示:可以用抛物线插值,,,。运行结果为:4.3用逐次分半梯形公式求积计算。运行结果为:程序代码4.1#include stdio.h#include math.h#include conio.h#define MAXSIZE 50void input(double x[MAXSIZE],double y[MAXSIZE],long n);void main(void){double x[MAXSIZE],y[MAXSIZE],_x,_y,t;long n,i,j;printf(\nPlease input the number of interpolation node:);scanf(%ld,n);input(x,y,n);printf(\nPlease input the interpolation point:);scanf(%lf,_x);_y=0;for(i=0;i=n-1;i++){ t=1;for(j=0;j=n-1;j++)if(j!=i) t*=(_x-x[j])/(x[i]-x[j]);_y+=t*y[i]; }printf(\nThe interpolation point(x,y)=(%lf,%lf).,_x,_y);getch();}void input(double x[MAXSIZE],double y[MAXSIZE],long n){long i;for(i=0;i=n-1;i++){printf(\nPlease input the interpolation node x[%ld],y[%ld]:,i,i);scanf(%lf,%lf,x[i],y[i]);}}4.2#include stdio.h#include math.h#include conio.h#define MAXSIZE 50void input(double x[MAXSIZE],double f[MAXSIZE][MAXSIZE],long n);void main(void){double x[MAXSIZE],f[MAXSIZE][MAXSIZE],_x,_y;long n,i,j;printf(\nPlease input the number of interpolation node:);scanf(%ld,n);input(x,f,n);printf(\nPlease input the interpolation point:);scanf(%lf,_x);for(j=1;j=n-1;j++)for(i=j;i=n-1;i++)f[i][j]=(f[i][j-1]-f[i-1][j-1])/(x[i]-x[i-j]);_y=f[n-1][n-1];for(i=n-2;i=0;i--)_y=f[i][i]+(_x-x[i])*_y;printf(\nThe interpolation point(x,y)=(%lf,%lf).,_x,_y);getch();}void input(double x[MAXSIZE],double f[MAXSIZE][MAXSIZE],long n){long i;for(i=0;i=n-1;i++){printf(\nPlease input the interpolation node x[%ld],y[%ld]:,i,i);scanf(%lf,%lf,x[i],f[i][0]);}}4.3#include stdio.h#include math.h#include conio.hdouble f(double x);void main(void){double x,a,b,h,tn,t2n,epsilon;long n,i;printf(\nPlease input the integral interval boundary a,b:);scanf(%lf,%lf,a,b);printf(\nPlease input
文档评论(0)