lagrange.cpp

来自「如何使用MATLAB编程」· C++ 代码 · 共 38 行

#include "StdAfx.h"
#include "lagrange.h"
CLagrange::CLagrange(void)
{
  y = 0;
}
CLagrange::~CLagrange(void)
{
}
double CLagrange::Compute(void)
{
 int i,j;
 double tempY;
 y = 0.0f;
 for (i = 0; i < nCount; i ++) 
 {
  tempY = 1.0f;
  for (j = 0; j < nCount; j ++) 
	{
		if (j != i) 
		{
			tempY *= (x-(pArr->GetAt(j)).x)/((pArr->GetAt(i)).x-(pArr->GetAt(j)).x);
		}
	}
	y += tempY * (pArr->GetAt(i)).y;
 }
	return y;
}
void CLagrange::SetX(double mx)
{
	x = mx;
}
void CLagrange::SetArray(fArray *arr)
{
	pArr   = arr;
	nCount = pArr->GetSize();
}