У меня есть кусок недокументированного кода, который я должен понять, чтобы исправить ошибку. Следующий метод называется optimization
и он должен найти максимум очень сложной функции f
. К сожалению, в некоторых случаях он терпит неудачу (т.е. Достигает линии "Максимальная итерация").
Я уже пытался написать некоторые модульные тесты, но это мало помогло.
Поэтому я хочу понять, как этот метод действительно работает, и если он реализует конкретный и хорошо известный алгоритм оптимизации. Возможно, тогда я смогу понять, если он подходит для решения требуемых уравнений.
public static double optimization(double x1, double x2, double x3, Function<Double, Double> f, double epsilon) {
double y1 = f.apply(x1);
double y2 = f.apply(x2);
double y3 = f.apply(x3);
double a = ( x1*(y2-y3)+ x2*(y3-y1)+ x3*(y1-y2)) / ((x1-x2)*(x1-x3)*(x3-x2));
double b = (x1*x1*(y2-y3)+x2*x2*(y3-y1)+x3*x3*(y1-y2)) / ((x1-x2)*(x1-x3)*(x2-x3));
int i=0;
do {
i=i+1;
x3=x2;
x2=x1;
x1=-1.*b/(2*a);
y1=f.apply(x1);
y2=f.apply(x2);
y3=f.apply(x3);
a = ( x1*(y2-y3)+ x2*(y3-y1)+ x3*(y1-y2))/((x1-x2)*(x1-x3)*(x3-x2));
b = (x1*x1*(y2-y3)+x2*x2*(y3-y1)+x3*x3*(y1-y2))/((x1-x2)*(x1-x3)*(x2-x3));
} while((Math.abs(x1 - x2) > epsilon) && (i<1000));
if (i==1000){
Log.debug("Max iteration reached");
}
return x1;
}