Newton 插值
来源:互联网 发布:淘宝怎样延迟久点收货 编辑:程序博客网 时间:2024/05/16 18:18
define MAXTIMES 10 //最高构造10次插值函数
typedef struct {
double x;
double y;
}Point; //插值节点类型
int Times; //实际要构造的次数
double DCList[(MAXTIMES + 1)*(MAXTIMES + 2) / 2]; //差商列表
Point knotList[MAXTIMES + 1];
void initKnotList(); //插值节点数据输入
void initDCList(); //构造差商表
void displayDCList();
double Newton(double u); //Newton 插值函数, 计算点x处的函数值
int main()
{
double u; //计算该点的插值函数值
initKnotList();initDCList();printf("%f \n", Newton(u));while (getchar() != '\n');getchar();return 0;
}
void initKnotList()
{
int i;
printf("input Times:");scanf_s("%d", &Times, sizeof(int));for (i = 0; i <= Times; i++){ scanf_s("%lf", &knotList[i].x, sizeof(double)); scanf_s("%lf", &knotList[i].y, sizeof(double));}
}
void initDCList()
{
int i, j;
int prePosi, currPosi;
currPosi = 0;for (i = 0; i <= Times; i++) //第i阶差商{ for (j = 0; j <= Times - i; ++j) //第i阶的第j个差商 { if (i == 0) //0阶差商 { DCList[currPosi + j] = knotList[j].y; } else { DCList[currPosi + j] = (DCList[prePosi + j] - DCList[prePosi + j + 1]) / (knotList[j].x - knotList[j + i].x); } } prePosi = currPosi; currPosi = (i + 1)*Times + 1;}
}
void displayDCList()
{
int i, j;
int currPosi;
currPosi = 0;for (i = 0; i <= Times; i++){ for (j = 0; j <= Times - i; ++j) { printf("%f ", DCList[currPosi+j]); } currPosi = ((i+1)*Times) + 1; printf("\n");}
}
double Newton(double u)
{
int i, j;
int DCNumber=0;
double sum = 0;
double w;
for (i = 0; i <= Times; i++){ w = 1; for (j = 0; j < i; ++j) { w *= (u - knotList[j].x); } sum += DCList[DCNumber] * w; DCNumber = ((i + 1)*Times + 1);}return sum;
}
- Newton插值函数
- Newton插值
- Newton 插值
- newton插值
- matlab Newton向前插值
- Matlab Newton向后插值
- newton插值 lagrange 插值
- 牛顿(Newton)插值及其MATLAB程序
- Newton均差插值性质证明
- Newton
- Newton
- 插值
- 插值
- 插值
- 插值
- 插值
- 插值
- 插值
- L-LSP和E-LSP介绍
- js有哪些内置对象????
- 盒图介绍
- javascript表单(正则表达式)验证大全
- IOS推送证书设置
- Newton 插值
- ubuntu下安装命令的区别(dpkg、apt-get、aptitude)
- Learning Swift: Optional Types(详解Optional)
- 设计模式之建造型-单例模式(1)
- 随手记
- Java基础第4讲--逻辑控制语句
- MySql--练习
- 永恒之蓝 端口关闭与开启 bat脚本
- 深度学习中的数学--lecture 1(2)