7,539
社区成员
发帖
与我相关
我的任务
分享请发表友善的回复…
发表回复
avisnet 2006-10-30
- 打赏
- 举报
x,y是两个double类型的数组，size是数组中的元素个数，val是要插值的点
avisnet 2006-10-30
- 打赏
- 举报
拉格朗日插值
double lagrange(const double* x, const double* y, unsigned int size, double val)
{
double sum = 0.0;
for(unsigned int k = 0; k < size; k++)
{
double mul = 1.0;
for(unsigned int j = 0; j < size; j++)
{
if(j != k)
{
mul *= (val - x[j]) / (x[k] - x[j]);
}
}
sum += y[k] * mul;
}
return sum;
}
double lagrange(const double* x, const double* y, unsigned int size, double val)
{
double sum = 0.0;
for(unsigned int k = 0; k < size; k++)
{
double mul = 1.0;
for(unsigned int j = 0; j < size; j++)
{
if(j != k)
{
mul *= (val - x[j]) / (x[k] - x[j]);
}
}
sum += y[k] * mul;
}
return sum;
}
