文档介绍：.
1 / 8
学生学号
实验课成绩
理工大学
学生实验报告书
实验课程名称 数值分析
开 课 学 院 计算机科学与技术学院
指导老师
学 生 姓 名
学生专业班级
2010—2
cout<<"请输入坐标："<<endl;
while<count1<count>
{
cin>>point[count1][0]>>point[count1][1];count1++;
.
5 / 8
}
cout<<"计算f<x>请输入x："<<endl;
double x;
cin>>x;
double d[100];
for<int i=0;i<count;i++>
d[i]=point[i][1];
for<int j=1;j<count;j++>
for<i=count-1;i>=j;i-->
d[i]=<d[i]-d[i-1]>/<point[i][0]-point[i-j][0]>;
double f=d[0],Lu=1,L;
for<i=1;i<count;i++>
{
Lu=Lu*<x-point[i-1][0]>;
L=Lu*d[i];
f=f+L;
}
cout<<"f<x>的值为："<<f<<endl;
return 0;
}
[埃米尔特插值法源程序]
#include<iostream>
using namespace std;
struct point
{
double x;
double y;
double d;
}point[100];
int main<>
{
cout<<"请输入坐标点个数："<<endl;
int count;
cin>>count;
int count1=0;
cout<<"请依次输入坐标的横纵坐标和对应的导数："<<endl;
while<count1<count>
{
cin>>point[count1].x>>point[count1].y>>point[count1].d;
count1++;
}
cout<<"计算f<x>请输入x："<<endl;
double x;
cin>>x;
.
6 / 8
double f=0,Lu=1,Laa=0,La,Lb;
for<int i=0;i<count;i++>
{
for<int j=0;j<count;j++>
{
if<j==i>continue;
Lu=Lu*<x-point[j].x>/<point[i].x-point[j].x>;
Laa=Laa+/<point[i].x-point[j].x>;
}
La=point[i].y*<1-2*<x-point[i].x>*Laa>*Lu*Lu;