张扬 2010212420数值计算方法实验.pdfVIP

  1. 1、本文档共23页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

实验一插值法

一、算法设计思想:

在插值计算中,取样点的多少往往会影响所得插值函数优化程度。一般情况

下,取样点越多所得插值函数越优化,对应的函数值与标准函数值越接近。

在MATLAB中实现分段线性插值,分段二次插值,拉格朗日插值的命令为

interp1,其调用格式为:Y1=interp1(X,Y,X1,’method’)

函数根据X,Y的值,计算函数在X1处的值。X,Y是两个等长的已知向量,分

别描述采样点和样本值,X1是一个向量或标量,描述欲插值点,Y1是一个

与X1等长的插值结果。

二、程序清单:

1.分段线性插值

#includestdio.h

#includemath.h

doubleLagrange2(double*x,double*y,doubleinput)

{

doubleoutput;

inti;

for(i=0;i5;i++)

{

if(x[i]=inputx[i+1]=input)

{

output=y[i]+(y[i+1]-y[i])*(input-x[i])/(x[i+1]-x[i]);

break;

}

}

returnoutput;

}

2.分段二次插值

doubleLagrange3(double*x,double*y,doubleu)

{

inti,k=0;

doublev;

for(i=0;i6;i++)

{

if(ux[1])

{

k=0;

v=y[k]*(u-x[k+1])*(u-x[k+2])/((x[k]-x[k+1])*(x[k]-x[k+2]))+y[k+1]*(u-

x[k])*(u-x[k+2])/((x[k+1]-x[k])*(x[k+1]-x[k+2]))+y[k+2]*(u-x[k])*(u-x

[k+1])/((x[k+2]-x[k])*(x[k+2]-x[k+1]));

}

if((x[i]uu=x[i+1])(fabs(u-x[i])=fabs(u-x[i+1])))

{

k=i-1;

v=y[k]*(u-x[k+1])*(u-x[k+2])/((x[k]-x[k+1])*(x[k]-x[k+2]))+y[k+1]*(u-

x[k])*(u-x[k+2])/((x[k+1]-x[k])*(x[k+1]-x[k+2]))+y[k+2]*(u-x[k])*(u-x

[k+1])/((x[k+2]-x[k])*(x[k+2]-x[k+1]));

}

if((x[i]uu=x[i+1])fabs(u-x[i])fabs(u-x[i+1]))

{

k=i;

v=y[k]*(u-x[k+1])*(u-x[k+2])/((x[k]-x[k+1])*(x[k]-x[k+2]))+y[k+1]*(u-

x[k])*(u-x[k+2])/((x[k+1]-x[k])*(x[k+1]-x[k+2]))+y[k+2]*(u-x[k])*(u-x

[k+1])/((x[k+2]-x[k])*(x[k+2]-x[k+1]));

}

if(ux[4])

{

k=3;

v=y[k]*(u-x[k+1])*(u-x[k+2])/((x[k]-x[k+1])*(x[k]-x[k+2]))+y[k+1]*(u-

x[k])*(u-x[k+2])/((x[k+1]-x[k])*(x[k+1]-x[k+2]))+y[k+2]*(u-x[k])*(u-x

[k+1])/((x[k+2]-x[k])*(x[k+2]-x[k+1]));

}

}

returnv;

}

voidmain()

{

doublex[6]={0.0,0.1,0.195,0.3,0.401,0.5},y[6]=

{0.39894,0.39695,0.39142,0.38138,0.36812,0.35206};

doubleu;

scanf(%lf,u);

拉格朗日插值:

doubleLagrange1(double*x,double*y,doublexx)

{

inti,j;

double*a,yy=0.000;

a=newdouble[6];

for(i=0;i

文档评论(0)

130****6553 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档