【C++】ExprTk：C++数学表达式解析利器

ExprTk：C++数学表达式解析利器

一、ExprTk简介

ExprTk是由Arash Partow开发的轻量级C++数学表达式解析库，支持实时解析和计算数学表达式。使用时只需要包含 .hpp文件即可

二、主要特性

1. 表达式支持：

// 支持复杂表达式
const std::string expr = "clamp(-1, sin(2πx) + cos(y/2π), +1)";

• 算术运算（±*/%^）
• 逻辑运算（&& || !）
• 300+内置函数（sin/cos/clamp等）
• 条件表达式（if-else）
• 字符串处理

1. 变量绑定（如您的示例所示）：

symbol_table.add_variable("x", x); // 绑定变量
symbol_table.add_constants();       // 添加π等常量

2. 自定义函数：

exprtk::function_compositor<double> fc(symbol_table);
fc.add("rect", "x,y", "x*y", "a,b"); // 添加矩形面积函数

三、使用示例

1. 条件表达式（温度控制系统）

std::string expr = 
    "if(temp > 35.0, 0.0, "          // 过热保护
    "if(temp < 20.0, 100.0, "        // 低温加热
    "(20.0 - temp) * 5.0))";         // 比例控制

double temp = 25.0;
symbol_table.add_variable("temp", temp);
expression_t expression;
parser.compile(expr, expression);

temp = 18.0;  // 输出：加热功率=100%
temp = 28.0;  // 输出：加热功率=(20-28)*5= -40（冷却）

2. 自定义函数（3D几何计算）

// 注册球体积计算函数
fc.add("sphere_volume", "r", "(4.0/3.0) * pi * r^3", "radius");

// 使用示例
std::string expr = "sphere_volume(radius) + cylinder_volume(5, height)";
symbol_table.add_variable("radius", 2.0);
symbol_table.add_variable("height", 10.0);

3. 批量计算优化

// 预编译表达式
exprtk::expression<double> expr;
parser.compile("x^2 + y^2", expr);

// 批量计算（10万次）
std::vector<std::pair<double, double>> data(100000);
for(auto& point : data) {
    symbol_table.variable_ref("x") = point.first;
    symbol_table.variable_ref("y") = point.second;
    double result = expr.value();
    //...存储结果
}

4. 错误处理增强

if (!parser.compile(expr, expression)) {
    for(size_t i = 0; i < parser.error_count(); ++i) {
        auto error = parser.get_error(i);
        std::cerr << "错误位置: " << error.token.position << "\n"
                  << "类型: " << error.diagnostic << "\n"
                  << "附近内容: " 
                  << std::string(error.token.value.data(), 20) << "\n";
    }
}
// 示例错误输出：
// 错误位置: 15
// 类型: ERR001 - 未定义变量 'var1'
// 附近内容: ...var1 + sin(x...

5. 物理引擎应用（抛物线运动）

std::string motion_expr = 
    "v0 * t * cos(theta) + 0.5 * ax * t^2;"  // X轴位移
    "v0 * t * sin(theta) - 0.5 * g * t^2";   // Y轴位移

symbol_table.add_constant("g", 9.8);  // 重力加速度
symbol_table.add_variable("v0", 50.0); // 初速度
symbol_table.add_variable("theta", 45.0 * pi / 180.0); // 角度转弧度

四、性能优化技巧

// 1. 表达式缓存
static std::unordered_map<std::string, exprtk::expression<double>> expr_cache;

if (!expr_cache.contains(expr_str)) {
    exprtk::expression<double> expr;
    parser.compile(expr_str, expr);
    expr_cache.emplace(expr_str, std::move(expr));
}
auto& expr = expr_cache[expr_str];

// 2. 内存预分配
symbol_table.reserve(16); // 预分配变量存储空间

// 3. 禁用扩展功能（提升15%速度）
parser_settings.disable_all_control_structures();
parser_settings.disable_all_syntax_checks();