ARTICLE DETAIL

建站实战干货

来自一线的建站与推广经验沉淀,每一条都经过真实交付验证。

不手写推导模型和 Jacobian,如何把 OCP/MPC 模型生成可部署 C++ SDK

2026/8/14 17:08:58 拓冰建站 浏览量
不手写推导模型和 Jacobian,如何把 OCP/MPC 模型生成可部署 C++ SDK

很多做机器人、自动驾驶、工业运动控制、能源调度的团队,都会遇到类似的问题:

  • 算法原型在 Python、MATLAB 或 CasADi 里能跑。

  • 真正上车、上机器人、上工业控制器时,又要改成 C++。

  • 模型、约束、代价函数一变,矩阵、导数、接口代码就要跟着改。

  • 最后工程里堆了很多手写 Jacobian、手写稀疏矩阵、手写参数映射和手写调试工具。

这篇文章分享一个实时 OCP/MPC 工程化思路:用 Python/CasADi 定义模型,自动生成 C++ 模型库和 runtime SDK,让目标端只运行 C++,不依赖 Python 和 CasADi。

这不是在讨论某个控制理论公式,而是讨论一个更工程化的问题:

如何把一个优化控制模型,从可运行原型,快速变成可部署、可验证、可诊断的 C++ SDK。

为什么手写 MPC 工程化很痛苦

MPC/OCP 问题本身通常不难描述。一个典型问题会包含:

  • 状态变量

  • 控制变量

  • 时变参数

  • 动力学模型

  • 阶段代价和终端代价

  • 等式约束和不等式约束

  • 软约束和 slack

  • warm start

  • 求解状态诊断

但工程落地时,真正消耗时间的往往不是“写出数学公式”,而是下面这些事情:

  • 把模型公式翻译成 C++。

  • 推导和维护 Jacobian / Hessian / 残差。

  • 保证稀疏结构和变量索引不出错。

  • 适配 x86、ARM64、QNX、嵌入式 Linux。

  • 输出 benchmark,确认平均耗时、最大耗时和迭代次数。

  • 求解慢或失败时,可以保存现场并离线 replay。

  • 给业务工程提供稳定的 C++ API 和 CMake 接入方式。

如果模型还在频繁变化,这些维护成本会不断放大。

一个更直接的工程化方式

一种更直接的方式是:

Python/CasADi 定义模型 -> 自动代码生成 -> C++ 模型库 -> runtime solver SDK -> CMake package -> 用户 C++ 工程接入

目标端只运行 C++,不依赖 Python/CasADi。

这条链路的核心价值是:

  • 不需要手动推导模型、矩阵和导数。

  • 自动生成模型专用 C++ 代码。

  • 支持从零建模,也支持已有 MPC/OCP 模型迁移。

  • 生成模型 SDK 和 runtime SDK。

  • 用户侧只接入 C++ API 和 CMake。

  • 支持 static/shared library。

  • 支持 toolchain 交叉编译。

  • 支持 benchmark、replay dump 和慢求解诊断。

一个最小例子应该是什么样

用户侧不应该关心底层求导和矩阵展开。比较理想的使用方式应该接近下面这样:

from rtocp_codegen import OcpModel ocp = OcpModel("MinimalLineTrackingModel") x = ocp.state("x") y = ocp.state("y") theta = ocp.state("theta") v = ocp.control("v") omega = ocp.control("omega") ref_x = ocp.parameter("ref_x") ref_y = ocp.parameter("ref_y") ocp.dynamics([ v * cos(theta), v * sin(theta), omega, ]) ocp.stage_cost((x - ref_x) ** 2 + (y - ref_y) ** 2 + 0.1 * omega ** 2)

然后生成 C++ SDK:

rtocp-codegen export-static-sdk model.py --install-dir ./model_sdk

在用户 C++ 工程中,通过 CMake 接入:

find_package(RTOCP REQUIRED) find_package(MinimalLineTrackingModel REQUIRED) target_link_libraries(app PRIVATE MinimalLineTrackingModel::minimal_line_tracking_model RTOCP::rtocp_solver )

最终业务代码只需要设置参数、调用求解、读取结果。

为什么不是只用通用非线性优化库

通用非线性优化库非常成熟,比如 Ceres、g2o、GTSAM、IPOPT、acados 等,在各自领域都很强。

但实时控制和产品交付里,很多时候问题不只是“求一个优化解”,还包括:

  • 模型如何快速变化。

  • 求导和稀疏结构如何维护。

  • 如何生成可部署 C++ 代码。

  • 如何接入客户 C++ 工程。

  • 如何在 ARM/QNX/嵌入式 Linux 上构建。

  • 如何记录每帧耗时、迭代次数和失败原因。

  • 出现慢求解或失败时如何 replay。

所以这个方向更像是“实时优化控制工程化 SDK”,而不是单纯的通用 solver。

典型适用场景

这套思路不局限于自动驾驶。

比较适合的场景包括:

  • AMR / AGV / 无人叉车的路径跟踪和动态障碍物避让。

  • 机器人局部运动控制。

  • 工业多轴运动控制。

  • 激光 / CNC 加工轨迹速度规划。

  • 半导体设备和精密平台运动控制。

  • 光储充 / 微电网 / 充电站功率分配。

  • 无人机安全着陆和能量约束控制。

这些问题的共同点是:

有动态系统 + 有约束 + 需要在线或准实时决策 + 手写规则和手写矩阵难维护

Demo 可以验证什么

为了验证这条链路,可以设计几类 demo:

动态障碍物绕行

展示车辆或机器人在滚动优化中避让障碍物,同时输出每帧迭代次数和求解耗时。

这个 demo 主要验证:

  • 动态约束建模。

  • 障碍物安全距离约束。

  • 连续帧 warm start。

  • 每帧耗时统计。

  • 慢求解诊断。

工业运动控制

展示复杂加工轮廓下的轨迹平滑、速度变化、加速度约束和求解耗时。

这个 demo 主要验证:

  • 非车、非机器人场景也能使用。

  • 工业轨迹的速度和加速度约束。

  • 曲线、拐角、复杂轮廓下的实时优化能力。

无人机安全着陆

展示无人机在低电量或受限能量条件下,根据高度、速度、电量和推力约束,生成安全着陆轨迹。

这个 demo 主要验证:

  • 高度、速度、推力和能量状态建模。

  • 近地速度限制和安全边界约束。

  • 能量 reserve 约束。

  • 推力变化率约束。

  • 安全控制类问题的实时优化能力。

充电站功率分配

展示多个充电车辆在总功率约束、电价变化和 SOC 目标下的功率分配。

这个 demo 主要验证:

  • 资源调度类问题也可以用 OCP/MPC 形式表达。

  • 不同阶段参数可以变化。

  • 可解释地输出功率、SOC、成本和约束状态。

更关注的不是“炫技”

对很多工程团队来说,最重要的不是一个 solver 名字,而是下面这些问题:

  • 半天能不能完成从 0 到 1 的模型 C++ SDK 闭环。

  • 能不能不用手写推导和矩阵展开。

  • 能不能接进已有 C++ 工程。

  • 能不能在目标平台构建。

  • 能不能看到平均耗时、最大耗时和迭代次数。

  • 求解异常时能不能复现。

所以一个更实际的验收方式是:

建模 -> 代码生成 -> C++ 接入 -> benchmark 验证

如果这个闭环能快速跑通,后续再讨论模型精度、约束设计、性能优化和目标平台部署。

适合交流的问题

如果你的团队正在做以下事情,可能适合交流:

  • 正在把 Python/MATLAB/CasADi 控制原型迁移到 C++。

  • 手写 MPC/OCP 工程维护成本高。

  • 需要在 ARM64、QNX、嵌入式 Linux 上部署实时优化。

  • 需要把控制模型交付成 SDK。

  • 需要 benchmark 和 replay 诊断能力。

  • 模型、约束、代价函数频繁变化。

如果只是 PID、LQR 或简单规则控制已经足够,这类工具可能并不是必须的。

结语

实时优化控制的难点,不只在求解器本身,也在从模型到工程交付的整条链路。

如果能把模型定义、自动求导、代码生成、C++ SDK、交叉编译、benchmark 和 replay 诊断串起来,MPC/OCP 的工程落地成本会明显下降。

如果你也在做 MPC/OCP 工程化、嵌入式部署或实时优化控制,欢迎在评论区交流。

相关产品:RTOCP,实时优化自动代码生成 SDK
官网:https://qiwentech.cn
支持:OCP/MPC 建模、C++ SDK 生成、嵌入式部署、性能验证
试用:可在官网提交 试用申请