- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
三次样条插值算法
//三次样条插值算法
#includeiostream.h
void main()
{
cout输入节点个数:;
int num;
cinnum;
int n=num-1;
double *x=new double[num];
double *y=new double[num];
cout按序依次输入各处节点及函数值endl;
for(int i=0;i=n;i++)
cinx[i]y[i];
double *h=new double [num];
for(i=0;i=n-1;i++)
h[i]=x[i+1]-x[i];
double *a1=new double[num];
double *b1=new double[num];
cout请输入边界。按1输入第一种边界条件,按2输入第二种边界条件: ;
int pd;
cinpd;
if(pd==1)
{
cout按顺序输入两端点的微商:endl;
int k1,k2;
cink1k2;
a1[0]=0;a1[n]=1;b1[0]=2*k1;b1[n]=2*k2;
}
else
{
cout按顺序输入两端点的二阶微商:endl;
int t1,t2;
cint1t2;
a1[0]=1;a1[n]=0;b1[0]=3*(y[1]-y[0])/h[0];b1[n]=3*(y[n]-y[n-1])/h[n-1];
}
for(i=1;i=n-1;i++)
{
a1[i]=h[i-1]/(h[i-1]+h[i]);
b1[i]=((1-a1[i])*(y[i]-y[i-1])/h[i-1]+a1[i]*(y[i+1]-y[i])/h[i])*3;
}
double*a2=new double[num];
double*b2=new double[num];
a2[0]=-a1[0]/2;b2[0]=b1[0]/2;
for(i=1;i=n;i++)
{
a2[i]=-a1[i]/(2+(1-a1[i])*a2[i-1]);
b2[i]=(b1[i]-(1-a1[i])*b2[i-1])/(2+(1-a1[i])*a2[i-1]);
}
double *m=new double[n+2];
m[n+1]=0;
for(i=n;i=0;i--)
{
m[i]=a2[i]*m[i+1]+b2[i];
}
cout输入所要求的点:endl;
int xx;
cinxx;
for(i=0;i=n-1;i++)
{
if(x[i]=xxxx=x[i+1]) break;
}
if(i=n-1)
{
double kk;
kk=(1+2*(xx-x[i])/(x[i+1]-x[i]))*(xx-x[i+1])*(xx-x[i+1])*y[i]/(x[i]-x[i+1])/(x[i]-x[i+1])+
(1+2*(xx-x[i+1])/(x[i]-x[i+1]))*(xx-x[i])*(xx-x[i])*y[i+1]/(x[i+1]-x[i])/(x[i+1]-x[i])+
(xx-x[i])*(xx-x[i+1])*(xx-x[i+1])*m[i]/(x[i]-x[i+1])/(x[i]-x[i+1])+
(xx-x[i+1])*(xx-x[i])*(xx-x[i])*m[i+1]/(x[i+1]-x[i])/(x[i+1]-x[i]);
coutxx的函数近似值为kkendl;
}
else
{
cout输入的值超出范围.endl;
}
}
用以上算法解问题1:给定y=f (x)的函数表如下,边界条件为s’’(75)=0,s’’(80)=0,用三次样条插值求f(78.3)的近似值。
x 75 76 77 78 79 80
f (x) 2.768 2.833 2.903 2.979 3.062 3.153
程序运行结果如下:
用以上算法解问题2:给定函数y=f (x)的函数表如下,边界条件s’(0.25)=1,s’(0.53)=0.6868,用三次样条插值算法求f (0.35)的近似值。
x 0.25 0.3 0.39 0.45 0.53
y=f (x) 0.5 0.5477 0.6245 0.6708 0.728
程序运行结果如下:
您可能关注的文档
- XXXMSVC方案.doc
- xu学习方法.doc
- xxx煤矿材料回收管理制度.doc
- WORKFLOW流程系统V1.4用户说明书.doc
- XX社区2012年文化工作计划.doc
- YJ材料台帐格式.doc
- zfp第七章(二阶).ppt
- z_C++语言试卷.doc
- WCDMA无线网络容量估算.ppt
- ××县党员志愿者管理办法(试行).doc
- 第十一章 电流和电路专题特训二 实物图与电路图的互画 教学设计 2024-2025学年鲁科版物理九年级上册.docx
- 人教版七年级上册信息技术6.3加工音频素材 教学设计.docx
- 5.1自然地理环境的整体性 说课教案 (1).docx
- 4.1 夯实法治基础 教学设计-2023-2024学年统编版九年级道德与法治上册.docx
- 3.1 光的色彩 颜色 电子教案 2023-2024学年苏科版为了八年级上学期.docx
- 小学体育与健康 四年级下册健康教育 教案.docx
- 2024-2025学年初中数学九年级下册北京课改版(2024)教学设计合集.docx
- 2024-2025学年初中科学七年级下册浙教版(2024)教学设计合集.docx
- 2024-2025学年小学信息技术(信息科技)六年级下册浙摄影版(2013)教学设计合集.docx
- 2024-2025学年小学美术二年级下册人美版(常锐伦、欧京海)教学设计合集.docx
文档评论(0)