1

我今天刚开始使用 boost，发现这篇文章非常有帮助。我正在尝试使用 boost::bisect 来求解一系列值的参数方程。如果我想解决 0.8 的值，则以下方法有效：

#include <boost/math/tools/roots.hpp>   

//declare constants
const double Y_r = 0.2126;
const double Y_g = 0.7152;
const double Y_b = 0.0722;

struct FunctionToApproximate  {
  double operator() (double x)  {
      return (pow (x, 2.4) * Y_r) + (pow ((x*0.6)+0.4, 2.4) * Y_g) + (1 * Y_b) - 0.8;
  }
};

struct TerminationCondition  {
  bool operator() (double min, double max)  {
    return fabs(min - max) <= t_c;
  }
};

using boost::math::tools::bisect;

std::pair<double, double> result = bisect(FunctionToApproximate(), 0.0, 1.0, TerminationCondition());

double root = (result.first + result.second) / 2;

我想把它包装在一个循环中，这样我就可以解决 0.8 以外的值。我该怎么办？

非常感谢！

4

1 回答 1

1

您可以添加一个状态/数据成员来FunctionToApproximate保存要在每次调用中减去的值：

struct FunctionToApproximate  {
    FunctionToApproximate(double val) :val(val){}
    double operator() (double x)  {
        return (pow (x, 2.4) * Y_r) + (pow ((x*0.6)+0.4, 2.4) * Y_g) + (1 * Y_b) - val;
    }
    double val;
};

然后将计算包装在一个循环中应该是直截了当的。

for(double i = 0.1; i <= 1.0; i+=1.0) {
    std::pair<double, double> result = bisect(FunctionToApproximate(i), 0.0, 1.0, TerminationCondition());

    double root = (result.first + result.second) / 2;
}
于 2018-05-06T01:50:11.457 回答