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