分段三次hermite插值 - 360文档中心

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

C语言课程设计报告

课程名称：程序设计语言

组员：学号：魏文豪135********

汤恒 135********

彭建平135******** 指导教师：李华刚

2015年 12 月 13 日

流程图

代码！

#include

#define m 4

#define n 5

void main()

{

int i,k;

float x[n+1],y[n+1],yy[n+1],h,z[m];

printf("请按行输入一系列的x值:\n");

for(k=0;k

scanf("%f",&x[k]);

printf("请按行输入一系列的y值:\n");

for(k=0;k

scanf("%f",&y[k]);

printf("请输入一系列的y'的值:\n");

for(k=0;k

scanf("%f",&yy[k]);

printf("请按行输入这%d个插值点:\n",m+1);

for(i=0;i

scanf("%f",&z[i]);

printf("%f\n",z[i]);

for(i=0;i

for(k=0;k

if(z[i]>=x[k]&&z[i]<=x[k+1])

{

h=pow((z[i]-x[k+1])/(x[k]-x[k+1]),2.0)*(1+2*(z[i]-x[k])/(x[k+1]-x[k]))*y[k]+pow( (z[i]-x[k])/(x[k+1]-x[k]),2.0)*(1+2*(z[i]-x[k+1])/(x[k]-x[k+1]))*y[k+1]+pow((z[i ]-x[k+1])/(x[k]-x[k+1]),2.0)*(z[i]-x[k])*yy[k]+pow((z[i]-x[k])/(x[k+1]-x[k]),2.0

)*(z[i]-x[k+1])*yy[k+1];

printf("h(%f)=%f\n",z[i],h);

}

运行图