根据零力矩点理论分析机器人行进过程的稳定条件，利用稳定裕度的概念，在支撑多边形中求取最优稳定点来规划零力矩点．可以为避免 walk步态中频繁调整躯干姿态导致的能耗和行进速度损失[1]

此外，“零力矩点”是判定仿人机器人动态稳定运动的重要指标，ZMP落在四足机器人支撑多边形的范围里面，则机器人可以稳定的行走。

前言

以下内容均为论文《基于零力矩点的四足机器人非平坦地形下步态规划与控制》的部分解读，想要跟深入了解请查看原文

一、支撑多边形

支撑多边形指的是四足机器人足底与地面之间所有接触点所构成的多边形区域。由于walk步态在低速行走时具有稳定性好、对越障和不平地形适应性强等优势，本文将主要对walk步态下的zmp进行研究。

首先我们要确定支撑多边形各个顶点的坐标值。因此需要先对四足机器人进行建模．

坐标轴

以基体质心为原点定义机器人坐标系{B}
X 轴指向机器人前进方向
Z 轴向上
Y 轴由右手定则确定

坐标系

将坐标系{B}投影至支撑面，我们得到世界参考坐标系{W}；

定义4条腿顺序编号（左前LF、右前RF、右后RH、左后LH）依次为１～４，坐标系分别为{LF}、{RF}、{RH}、{LH}。

各个坐标系位置如下图所示。根据关节角度值由正运动学求出末端点 $p_{ei}$ ，其中i = 1,2,3,4。在单腿坐标系{i}下的坐标表示 $^iP_{ei}$ 。根据坐标系间的几何关系，先求出{B}中的各个足端坐标 $^BP_{ei}$ ，然后在根据坐标系间的变换求出相对于{W}的足端位置 $^WP_{ei}$ ，最后由各支撑点坐标确定出支撑多边形．

图片来自参考文献[1]

二、零力矩点

在世界参考系{W}下，假设四足机器人的总重量为Ｍ，

所受重为 $g=[0, 0, -g]^T$ ，
中心坐标为 $p_c = [x_c, y_c, z_c]^T$ ，
zmp所在位置 $p=[x_p, y_p, z_p]^T$ ，
支撑点受到地面反作用力合力为 $f$ ，

则地面反作用力绕远点的力 $\tau$ 为：

$\tau = p\times f + \tau_p$

其中 $\tau_p$ 为过零力矩点的力矩。

动量定理：如果一个系统不受外力或所受外力的矢量和为零，那么这个系统的总动量保持不变，公式为： $F_t=mv′－mv=p′－p$

角动量定理：质点系对一点（或一轴）的角动量对时间的导数等于外力系对此点（或此轴）的主矩

根据动量定理和角动量定理，有以下关系：

$\dot{P} = Mg + f$

$\dot{L} = p_c \times Mg + \tau$

综合上式，消去 $\tau$ 和 $f$ 得：

$\tau_p = \dot{L} - p_c \times Mg + (\dot{P}-Mg)\times p$

由ZMP理论，令上式中的水平分量 $\tau_{px}$ 和 $\tau_{py}$ 为０，解得ZMP位置坐标 $x_p$ 和 $y_p$ ，

$\left\{\begin{matrix} x_p = \frac{Mgx_c + z_p\dot{P_x} -\dot{L_y}}{Mg + \dot{P_z}}\\ \\ y_p = \frac{Mg y_c + z_p\dot{P_y} + \dot{L_x}}{Mg + \dot{P_z}} \end{matrix}\right.\tag{1}$

其中 $z_p$ 为支撑面离地高度，当机器人出于平面时， $z_p=0$ 。把四足机器人简化成一个质点，其动量和角动量分别为：

$\left\{\begin{matrix} P = M\dot{p_c} = M[\dot{x_c} \quad \dot{y_c} \quad \dot{z_c} ]^T\\ \\ L = p_c \times M\dot{p_c} = M[\dot{x_c} \quad \dot{y_c} \quad \dot{z_c} ]^T \end{matrix}\right. \tag{2}$

将（2）代入（1）中：

$\left\{\begin{matrix} x_p = x_c - \frac{(z_c - z_p)\ddot{x_c}}{\ddot{z_c}+g}\\ \\ y_p = y_c - \frac{(z_c-z_p)\ddot{y_c}}{\ddot{z_c}+ g} \end{matrix}\right.$