Например, у меня есть точки
Y X
100 50
90 43
80 32
необходимо решить для y = 50
или
Y X
1/1/2009 100
1/3/2009 97
1/4/2009 94
1/5/2009 92
1/6/2009 91
1/7/2009 89
нужно решить для y = 1/23/2009
Например, у меня есть точки
Y X
100 50
90 43
80 32
необходимо решить для y = 50
или
Y X
1/1/2009 100
1/3/2009 97
1/4/2009 94
1/5/2009 92
1/6/2009 91
1/7/2009 89
нужно решить для y = 1/23/2009
Я использую численный компонент Math.NET http://numerics.mathdotnet.com/
Он содержит "различные методы интерполяции, включая барицентрические подходы и сплайны".
Но, как говорится, есть ложь, проклятая ложь и бикубические сплайновые интерполяции.
Посмотрите, можете ли вы найти то, что хотите, ALGLIB. Конечно, вам все равно придется принимать решения о соответствующем типе интерполяции/экстраполяции для вашей проблемы.
Я не знаю о библиотеках, но здесь простой решатель Secant:
class SecantSolver
{
private int _maxSteps= 10;
private double _precision= 0.1;
public SecantSolver(int maxSteps, double precision)
{
_maxSteps= maxSteps;
_precision= precision;
if (maxSteps <= 0)
throw new ArgumentException("maxSteps is out of range; must be greater than 0!");
if (precision <= 0)
throw new ArgumentException("precision is out of range; must be greater than 0!");
}
private double ComputeNextPoint(double p0, double p1, Func<Double,Double> f)
{
double r0 = f(p0);
double r1 = f(p1);
double p2 = p1 - r1 * (p1-p0) / (r1-r0); // the basic secant formula
return p2;
}
public double Solve( double lowerBound, double upperBound, Func<Double,Double> f, out String message)
{
double p2,p1,p0;
int i;
p0=lowerBound;
p1=upperBound;
p2= ComputeNextPoint(p0,p1,f);
// iterate till precision goal is met or the maximum
// number of steps is reached
for(i=0; System.Math.Abs(f(p2))>_precision &&i<_maxSteps;i++) {
p0=p1;
p1=p2;
p2=ComputeNextPoint(p0,p1,f);
}
if (i < _maxSteps)
message = String.Format("Method converges in " + i + " steps.");
else
message = String.Format("{0}. The method did not converge.", p2);
return p2;
}
}
Использование:
SecantSolver solver= new SecantSolver(200, // maxsteps
0.00000001f/100 // tolerance
);
string message;
double root= solver.Solve(0.10, // initial guess (lower)
1.0, // initial guess (upper)
f, // the function to solve
out message
);
ILNumerics Interpolation Toolbox содержит все стандартные функции интерполяции. Вы найдете различные функции интерполяции/экстраполяции с общим простым интерфейсом. Особым преимуществом перед другими решениями является скорость: эти функции тщательно оптимизируются настолько быстро, насколько это возможно, для многоядерных аппаратных средств и для больших данных.