論文導讀：MIT Cheetah 3 凸優化 MPC——1 毫秒內解出地面反作用力

Jared Di Carlo, Patrick M. Wensing, Benjamin Katz, Gerardo Bledt, Sangbae Kim · 2018原文連結

論文資訊卡

標題：Dynamic Locomotion in the MIT Cheetah 3 Through Convex Model-Predictive Control
作者：Jared Di Carlo、Patrick M. Wensing、Benjamin Katz、Gerardo Bledt、Sangbae Kim（MIT）
發表：IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), 2018

問題背景

見模型預測控制入門一文討論過的核心困難：四足機器人的步態涉及腳何時觸地、觸地瞬間的衝擊力、摩擦力錐約束，這些接觸動力學本質上是高度非線性甚至非光滑的，完整的非線性 MPC（NMPC）雖然理論上能處理得更精確，但計算成本高，很難在四足機器人需要的高控制頻率下即時求解。這篇論文要解決的問題是：能不能在保留足夠物理真實性的前提下，把這個問題簡化成計算成本低廉的凸優化（convex optimization）問題，讓 MPC 真正能在實體機器人上以夠高的頻率即時運行。

方法核心：簡化模型換取可即時求解的凸問題

論文的核心做法是把機器人動力學簡化成一個「單一剛體」（single rigid body）模型——忽略四條腿本身的慣性與動力學細節，只把機身視為一個受到各腳底地面反作用力作用的剛體，用這個簡化模型去規劃每隻腳在觸地期間該輸出多大的地面反作用力，才能讓機身沿著期望的軌跡運動。這個簡化雖然犧牲了腿部動力學的精確度，但關鍵的好處是：在這個簡化模型下，力規劃問題可以被表述成一個凸優化問題（優化變數是各腳的地面反作用力，約束包含摩擦力錐這類凸集合限制）——凸問題保證只要有解就是全域最佳解，而且有極為高效、成熟的求解器可用，這是論文能達到極端求解速度的根本原因。

論文報告在最長 0.5 秒的預測時域下，這個地面反作用力規劃問題能在 1 毫秒以內求解到最佳解，並以 20 到 30 Hz 的頻率持續重新求解——這正是遞推時域（receding horizon）策略在實體硬體上真正落地所需要的求解速度等級，遠快於完整非線性動力學模型的求解時間。

實驗結果

儘管用的是簡化模型，MIT Cheetah 3 憑藉這套凸優化 MPC 框架展現出涵蓋多種步態的穩健運動能力：站立、小跑（trot）、飛奔小跑（flying-trot）、彈跳（pronk）、綁跳（bound）、側步（pace）、三腳步態，以及完整 3D 的奔馳（gallop）——同一套控制框架能透過調整步態時序參數，涵蓋這麼廣泛的運動模式範圍，展現了地面反作用力規劃這個抽象層級，確實抓住了四足運動控制裡足夠通用的核心結構。

影響與定位

這篇論文確立的「單一剛體模型 + 凸優化 MPC」框架，是理解現代四足機器人商用產品（包括 Boston Dynamics Spot、Unitree 系列等）控制系統設計思路的重要參照——它證明了不需要完整求解高度非線性的全身動力學，用一個經過刻意簡化、但仍保留關鍵物理結構的模型，就能在真實硬體的計算與時間預算內做到可靠的動態運動控制。這個「用簡化換取即時可解性」的設計哲學，直接呼應模型預測控制入門一文中討論線性 MPC 與非線性 MPC 取捨的核心論點，也是後續大量凸 MPC 足式機器人控制研究工作的共同起點。