Frenet 坐标下的车辆运动学

以后轴中心的最近点投影定义车辆 Frenet 状态,在无侧滑运动学假设下分别推导时间域三状态模型与参考路径弧长域二状态模型,并说明弧长换算、坐标有效域及倒车时的符号约定。

首次发布
最后更新

本文目录

基本术语与 Frenet 状态定义

本文只讨论车辆的运动学:描述位置和航向如何随运动变化,不建立轮胎力、载荷转移或车辆质量等动力学模型。车辆在平面内运动,并以后轴中心作为位置参考点。

本文所说的状态,是描述车辆相对于参考路径的几何位置与朝向、并随自变量演化的一组变量;输入是由模型外部给定、用于决定状态变化率的量。以时间 tt 为自变量的模型称为时间域模型,以参考路径弧长 ss 为自变量的模型称为参考路径弧长域模型,也常简称为空间域模型。

参考路径是一条以弧长 ss 参数化的 C2C^2 平面曲线

r:I⟶R2,s⟼r(s).\mathbf r:I\longrightarrow\mathbb R^2, \qquad s\longmapsto\mathbf r(s).

这里 II 是弧长参数 ss 的取值区间。采用直角坐标平面的标准定向,令 JJ 表示逆时针旋转 90∘90^\circ 的线性算子。参考路径的单位切向量、左单位法向量、有符号曲率与连续展开的参考航向角分别记为 T\mathbf T、N\mathbf N、κr\kappa_r 与 θr\theta_r。本文使用下列平面有符号 Frenet 关系:

T(s)=drds,N(s)=JT(s),dTds=κrN,dNds=−κrT,κr=dθrds.\begin{aligned} \mathbf T(s) &= \frac{\mathrm{d}\mathbf r}{\mathrm{d}s},& \mathbf N(s) &= J\mathbf T(s),\\[0.6em] \frac{\mathrm{d}\mathbf T}{\mathrm{d}s} &= \kappa_r\mathbf N,& \frac{\mathrm{d}\mathbf N}{\mathrm{d}s} &= -\kappa_r\mathbf T,\\[0.6em] \kappa_r &= \frac{\mathrm{d}\theta_r}{\mathrm{d}s}. \end{aligned}

曲率符号与左法向约定保持一致:沿 ss 增大的方向经过参考路径时,左转对应 κr>0\kappa_r>0,右转对应 κr<0\kappa_r<0。弧长参数、标架、标准定向和上述公式的完整构造见 《平面曲线的有符号 Frenet 标架》;本文只将其作为车辆运动学的几何基础。

z(t)\mathbf z(t) 表示车辆后轴中心在固定平面坐标系中的位置向量,车辆航向角记为 θ(t)\theta(t),相应的航向单位向量定义为

eθ(t):=[cos⁡θ(t)sin⁡θ(t)].\mathbf e_\theta(t):= \begin{bmatrix} \cos\theta(t)\\ \sin\theta(t) \end{bmatrix}.

沿 eθ\mathbf e_\theta 方向取正的有符号纵向速度记为 v(t)v(t)。无侧滑是指后轴中心的瞬时速度没有垂直于车身纵轴的分量,即

z˙=veθ.\dot{\mathbf z}=v\mathbf e_\theta.

前进时 v>0v>0,倒车时 v<0v<0,停车时 v=0v=0。因此 eθ\mathbf e_\theta 始终描述车身朝向;车辆运动时的实际速度方向还取决于 vv 的符号。

以 ℓ\ell 表示后轴中心轨迹的有向弧长坐标。给定初值 ℓ(t0)\ell(t_0) 后,它定义为

ℓ(t)=ℓ(t0)+∫t0tv(τ) dτ,\ell(t)=\ell(t_0)+\int_{t_0}^{t}v(\tau)\,\mathrm{d}\tau,

因而

dℓdt=v.\frac{\mathrm{d}\ell}{\mathrm{d}t}=v.

前进时 ℓ\ell 增大,倒车时 ℓ\ell 减小。它不是只增不减的累计路程;后者的增量应由 ∣v∣ dt|v|\,\mathrm{d}t 计算。在 v≠0v\ne0 的轨迹段上,ℓ\ell 可作为局部参数,链式法则给出 dz/dℓ=eθ\mathrm{d}\mathbf z/\mathrm{d}\ell=\mathbf e_\theta,因此车辆轨迹相对于这一取向的有符号曲率定义为

κv:=dθdℓ.\kappa_v := \frac{\mathrm{d}\theta}{\mathrm{d}\ell}.

对给定的车辆位置,最近点投影是在参考路径上寻找与后轴中心距离最小的点。本文假设该投影点唯一并位于参考路径参数区间的内部,因此后轴中心与投影点的连线垂直于参考路径。投影点的弧长坐标记为 s(t)s(t);从投影点沿左单位法向量 N\mathbf N 测得的有符号距离记为横向偏差 d(t)d(t),其中 d>0d>0 表示车辆位于参考路径 ss 增大方向的左侧,d<0d<0 表示位于右侧。车辆航向角与投影点处参考航向角之差定义为

φ(t):=θ(t)−θr(s(t)).\varphi(t) := \theta(t)-\theta_r\bigl(s(t)\bigr).

这里的 θ\theta 与 θr\theta_r 都取沿所讨论轨迹段连续展开的角度值,使 φ\varphi 在跨越 ±π\pm\pi 的表示边界时不发生人为跳变。时间域模型以时间 tt 为自变量、以 (s,d,φ)(s,d,\varphi) 为状态;参考路径弧长域模型也称空间域模型,它改用 ss 作为自变量,因此只需以 (d,φ)(d,\varphi) 为状态。图 1 表示这些量之间的几何关系:T(s)\mathbf T(s) 与 N(s)\mathbf N(s) 是投影点处的 Frenet 基向量,紫色双向尺寸线表示从参考路径到后轴中心的横向偏差 dd,后轴中心处的箭头 eθ\mathbf e_\theta 表示车辆航向单位向量。

本文以直立的 d\mathrm{d} 表示微分符号,以区别于斜体的横向偏差变量 dd;d˙\dot d 表示横向偏差对时间的导数。

T(s)N(s)eθd

‖T‖ = ‖N‖ = 1

图 1. 车辆 Frenet 状态的几何定义。蓝点是车辆在参考路径上的投影点,T(s)\mathbf T(s) 与 N(s)\mathbf N(s) 是该点处的切向和左法向基向量;后轴中心满足 z=r(s)+dN(s)\mathbf z=\mathbf r(s)+d\mathbf N(s),紫色双向尺寸线表示横向偏差 dd,eθ\mathbf e_\theta 表示车辆航向单位向量。

时间域 Frenet 运动学

位置坐标运动学

根据图 1 中的法向坐标表示,车辆位置为

z(t)=r(s(t))+d(t)N(s(t)).(1)\mathbf z(t) = \mathbf r\bigl(s(t)\bigr) + d(t)\mathbf N\bigl(s(t)\bigr). \tag{1}

对式 (1) 关于时间 tt 求导,并使用 dr/ds=T\mathrm{d}\mathbf r/\mathrm{d}s=\mathbf T 与 dN/ds=−κrT\mathrm{d}\mathbf N/\mathrm{d}s=-\kappa_r\mathbf T:

z˙=drdss˙+d˙ N+d dNdss˙=Ts˙+d˙ N−dκrTs˙=(1−dκr)s˙ T+d˙ N.\begin{aligned} \dot{\mathbf z} &= \frac{\mathrm{d}\mathbf r}{\mathrm{d}s}\dot s + \dot d\,\mathbf N + d\,\frac{\mathrm{d}\mathbf N}{\mathrm{d}s}\dot s\\[0.6em] &= \mathbf T\dot s + \dot d\,\mathbf N - d\kappa_r\mathbf T\dot s\\[0.6em] &= (1-d\kappa_r)\dot s\,\mathbf T + \dot d\,\mathbf N. \end{aligned}

这里及下文未显式写出的 T\mathbf T、N\mathbf N 与 κr\kappa_r 均在投影点 s(t)s(t) 处取值。无侧滑假设给出 z˙=veθ\dot{\mathbf z}=v\mathbf e_\theta;前进时 v>0v>0,倒车时 v<0v<0,这一关系保持不变。由于 (T,N)(\mathbf T,\mathbf N) 是一组正交归一基,eθ\mathbf e_\theta 可按它在两个基向量方向上的正交投影展开:

eθ=⟨eθ,T⟩T+⟨eθ,N⟩N.\mathbf e_\theta = \langle\mathbf e_\theta,\mathbf T\rangle\mathbf T + \langle\mathbf e_\theta,\mathbf N\rangle\mathbf N.

将

T=[cos⁡θrsin⁡θr],N=[−sin⁡θrcos⁡θr]\mathbf T = \begin{bmatrix} \cos\theta_r\\ \sin\theta_r \end{bmatrix}, \qquad \mathbf N = \begin{bmatrix} -\sin\theta_r\\ \cos\theta_r \end{bmatrix}

代入两个投影系数,得到

⟨eθ,T⟩=cos⁡θcos⁡θr+sin⁡θsin⁡θr=cos⁡(θ−θr)=cos⁡φ,⟨eθ,N⟩=−cos⁡θsin⁡θr+sin⁡θcos⁡θr=sin⁡(θ−θr)=sin⁡φ.\begin{aligned} \langle\mathbf e_\theta,\mathbf T\rangle &= \cos\theta\cos\theta_r + \sin\theta\sin\theta_r\\ &= \cos(\theta-\theta_r) = \cos\varphi,\\[0.7em] \langle\mathbf e_\theta,\mathbf N\rangle &= -\cos\theta\sin\theta_r + \sin\theta\cos\theta_r\\ &= \sin(\theta-\theta_r) = \sin\varphi. \end{aligned}

因此,车辆航向单位向量在 Frenet 标架中的表示为

eθ=cos⁡φ T+sin⁡φ N.\mathbf e_\theta = \cos\varphi\,\mathbf T + \sin\varphi\,\mathbf N.

车辆速度由此写成

z˙=veθ=vcos⁡φ T+vsin⁡φ N.\dot{\mathbf z} = v\mathbf e_\theta = v\cos\varphi\,\mathbf T + v\sin\varphi\,\mathbf N.

图 2 在相同的 Frenet 标架和相同的合速度箭头下,分别画出位置坐标求导与无侧滑速度表示所给出的两个分量。这样的并列使后续的系数比较成为直接的图形对应。

(a)
z˙\dot{\mathbf z}(1−dκr)s˙ T(1-d\kappa_r)\dot{s}\,\mathbf Td˙ N\dot{d}\,\mathbf N
(b)
v eθv cosφ Tv sinφ N

图 2. 同一车辆速度向量在同一 Frenet 标架中的两种分解。(a)对位置坐标求导得到切向分量 (1−dκr)s˙ T(1-d\kappa_r)\dot s\,\mathbf T 与法向分量 d˙ N\dot d\,\mathbf N;(b)无侧滑约束得到切向分量 vcos⁡φ Tv\cos\varphi\,\mathbf T 与法向分量 vsin⁡φ Nv\sin\varphi\,\mathbf N。由于标准正交基下的分解唯一,对应分量的标量系数必须分别相等。

比较同一速度向量在基 (T,N)(\mathbf T,\mathbf N) 下的系数,得到

(1−dκr)s˙=vcos⁡φ,d˙=vsin⁡φ.(1-d\kappa_r)\dot s = v\cos\varphi, \qquad \dot d = v\sin\varphi.

当 1−dκr≠01-d\kappa_r\ne0 时,时间域位置运动学为

s˙=vcos⁡φ1−dκr,d˙=vsin⁡φ.\dot s = \frac{v\cos\varphi}{1-d\kappa_r}, \qquad \dot d = v\sin\varphi.

式 (1) 的法向坐标与速度分解见 [1, §2.1–2.2, Eqs. (14)–(20)],等价的位置表示见 [2, §3.1, Eq. (1)];相同坐标约定下的纵向与横向状态方程见 [3, §4.3, Eqs. (4.33a)–(4.33b)]。

航向误差运动学

航向误差定义为

φ(t)=θ(t)−θr(s(t)).\varphi(t) = \theta(t) - \theta_r\bigl(s(t)\bigr).

关于时间使用链式法则,并代入 κr=dθr/ds\kappa_r=\mathrm{d}\theta_r/\mathrm{d}s:

φ˙=θ˙−dθrdss˙=θ˙−κrs˙.\begin{aligned} \dot\varphi &= \dot\theta - \frac{\mathrm{d}\theta_r}{\mathrm{d}s}\dot s\\[0.6em] &= \dot\theta - \kappa_r\dot s. \end{aligned}

在 v≠0v\ne0 的运动区间,由 κv=dθ/dℓ\kappa_v=\mathrm{d}\theta/\mathrm{d}\ell 以及 dℓ/dt=v\mathrm{d}\ell/\mathrm{d}t=v,

θ˙=dθdℓdℓdt=vκv.\begin{aligned} \dot\theta &= \frac{\mathrm{d}\theta}{\mathrm{d}\ell}\frac{\mathrm{d}\ell}{\mathrm{d}t}\\[0.6em] &= v\kappa_v. \end{aligned}

当 v=0v=0 时,ℓ\ell 不能在该时刻作为 tt 的局部替代参数;若曲率输入 κv\kappa_v 有界,运动学关系 θ˙=vκv\dot\theta=v\kappa_v 可连续延拓到停车状态,此时 θ˙=0\dot\theta=0。

因此航向误差运动学为

φ˙=vκv−κrs˙.\dot\varphi = v\kappa_v - \kappa_r\dot s.

三状态时间域模型

将航向误差运动学与位置坐标运动学合并,得到状态 (s,d,φ)(s,d,\varphi) 的三状态时间域模型:

s˙=vcos⁡φ1−dκr,d˙=vsin⁡φ,φ˙=vκv−s˙ κr=v(κv−κrcos⁡φ1−dκr).(2)\begin{aligned} \dot s &= \frac{v\cos\varphi}{1-d\kappa_r},\\[0.65em] \dot d &= v\sin\varphi,\\[0.65em] \dot\varphi &= v\kappa_v - \dot s\,\kappa_r\\[0.65em] &= v\left( \kappa_v-\frac{\kappa_r\cos\varphi}{1-d\kappa_r} \right). \end{aligned} \tag{2}

这里 κr=κr(s)\kappa_r=\kappa_r(s) 由参考路径给定。本文把有符号纵向速度 vv 与车辆轨迹曲率 κv\kappa_v 作为运动学模型的输入;在具体的转向模型中,κv\kappa_v 通常再由转向几何确定。模型有效要求最近点投影唯一且 1−dκr≠01-d\kappa_r\ne0:前者使车辆位置对应唯一的 (s,d)(s,d),后者是法向坐标映射在切向方向上不退化的局部正则条件。本文将同时满足这些条件的区域称为 Frenet 坐标有效域。与式 (2) 对应的标准 Frenet 车辆模型见 [3, §4.3, Eq. (4.41)]。

时间域模型中的倒车

式 (2) 同样适用于倒车,无须针对倒车重新推导。前述时间域推导已经把 vv 定义为沿航向单位向量 eθ\mathbf e_\theta 取正的有符号纵向速度,并使用统一的无侧滑运动学关系

z˙=veθ,θ˙=vκv.\dot{\mathbf z}=v\mathbf e_\theta, \qquad \dot\theta=v\kappa_v.

这两个关系没有要求 v>0v>0。前进时 v>0v>0,速度向量 z˙\dot{\mathbf z} 与车辆航向单位向量 eθ\mathbf e_\theta 同向;倒车时 v<0v<0,z˙\dot{\mathbf z} 与 eθ\mathbf e_\theta 反向。位置坐标关系 (1) 对时间求导、在 Frenet 标架中比较分量以及对航向误差使用链式法则的过程都不依赖 vv 的正负,因此倒车时仍然得到同一个三状态模型。

记 sgn⁡(x)\operatorname{sgn}(x) 为符号函数:x>0x>0、x=0x=0 与 x<0x<0 时,其值分别为 11、00 与 −1-1。由式 (2),参考路径弧长坐标的变化方向一般满足

sgn⁡(s˙)=sgn⁡(v)sgn⁡ ⁣(cos⁡φ1−dκr).\operatorname{sgn}(\dot s)= \operatorname{sgn}(v) \operatorname{sgn}\!\left(\frac{\cos\varphi}{1-d\kappa_r}\right).

在通常的 Frenet 跟踪区域 1−dκr>01-d\kappa_r>0 且 ∣φ∣<π/2|\varphi|<\pi/2 内,前进对应 s˙>0\dot s>0,倒车对应 s˙<0\dot s<0。更一般的情形应直接由上式判断,不能脱离 φ\varphi 和 1−dκr1-d\kappa_r 的符号仅凭“倒车”断定 s˙\dot s 的符号。

为进行严格对比,固定同一状态 (s,d,φ)(s,d,\varphi)、同一车辆曲率 κv\kappa_v 和同一速度大小 u>0u>0,只把纵向速度从前进的 v=uv=u 改为倒车的 v=−uv=-u。各量的变化如下:

量前进:v=uv=u倒车:v=−uv=-u定义或模型是否改变
车辆速度 z˙\dot{\mathbf z}ueθu\mathbf e_\theta−ueθ-u\mathbf e_\theta无侧滑关系不变
有向弧长变化率 dℓ/dt\mathrm{d}\ell/\mathrm{d}tuu−u-u始终有 dℓ/dt=v\mathrm{d}\ell/\mathrm{d}t=v
参考路径坐标变化率 s˙\dot sucos⁡φ1−dκr\dfrac{u\cos\varphi}{1-d\kappa_r}−ucos⁡φ1−dκr-\dfrac{u\cos\varphi}{1-d\kappa_r}由同一状态方程决定
横向偏差变化率 d˙\dot dusin⁡φu\sin\varphi−usin⁡φ-u\sin\varphidd 的定义不变
航向角速度 θ˙\dot\thetauκvu\kappa_v−uκv-u\kappa_vκv\kappa_v 不人为变号
航向误差变化率 φ˙\dot\varphiu(κv−κrcos⁡φ1−dκr)u\left(\kappa_v-\dfrac{\kappa_r\cos\varphi}{1-d\kappa_r}\right)−u(κv−κrcos⁡φ1−dκr)-u\left(\kappa_v-\dfrac{\kappa_r\cos\varphi}{1-d\kappa_r}\right)φ\varphi 的定义不变
参考曲率 κr\kappa_r 与车辆曲率 κv\kappa_v使用既定符号约定使用同一符号约定均不修改
式 (2)使用原方程使用原方程完全不变

因此,在上述比较条件下,仅改变 vv 的符号会使式 (2) 的三个状态变化率同时反号。时间 tt 仍按正方向推进;与参考路径弧长域模型不同,时间域倒车通过 v<0v<0 表示,而不是通过反向积分时间实现。停车时 v=0v=0,式 (2) 给出 s˙=d˙=φ˙=0\dot s=\dot d=\dot\varphi=0。

参考路径弧长域 Frenet 运动学

本节直接从弧长参数化的位置坐标关系与航向误差关系建立模型,不调用时间域状态方程,也不使用关于时间的导数。

车辆轨迹的有向弧长坐标

变量 ss 是参考路径上投影点的弧长坐标,ℓ\ell 是沿车辆航向单位向量 eθ\mathbf e_\theta 取正的后轴中心轨迹有向弧长坐标。在 Frenet 坐标有效且车辆非静止的轨迹段上,ss 和 dd 都可视为 ℓ\ell 的函数,式 (1) 相应写成

z(ℓ)=r(s(ℓ))+d(ℓ)N(s(ℓ)).\mathbf z(\ell) = \mathbf r\bigl(s(\ell)\bigr) + d(\ell)\mathbf N\bigl(s(\ell)\bigr).

直接关于 ℓ\ell 求导:

dzdℓ=drdsdsdℓ+d ddℓN+ddNdsdsdℓ=Tdsdℓ+d ddℓN+ddNdsdsdℓ.\begin{aligned} \frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell} &= \frac{\mathrm{d}\mathbf r}{\mathrm{d}s}\frac{\mathrm{d}s}{\mathrm{d}\ell} + \frac{\mathrm{d}\,d}{\mathrm{d}\ell}\mathbf N + d\frac{\mathrm{d}\mathbf N}{\mathrm{d}s}\frac{\mathrm{d}s}{\mathrm{d}\ell}\\[0.8em] &= \mathbf T\frac{\mathrm{d}s}{\mathrm{d}\ell} + \frac{\mathrm{d}\,d}{\mathrm{d}\ell}\mathbf N + d\frac{\mathrm{d}\mathbf N}{\mathrm{d}s}\frac{\mathrm{d}s}{\mathrm{d}\ell}. \end{aligned}

由 Frenet 关系 dN/ds=−κrT\mathrm{d}\mathbf N/\mathrm{d}s=-\kappa_r\mathbf T 与 ⟨T,N⟩=0\langle\mathbf T,\mathbf N\rangle=0,

⟨dNds,N⟩=−κr⟨T,N⟩=0.\left\langle \frac{\mathrm{d}\mathbf N}{\mathrm{d}s}, \mathbf N \right\rangle =-\kappa_r\langle\mathbf T,\mathbf N\rangle=0.

将位置导数投影到 N\mathbf N 方向,并使用 ⟨T,N⟩=0\langle\mathbf T,\mathbf N\rangle=0 与 ⟨N,N⟩=1\langle\mathbf N,\mathbf N\rangle=1,得到

⟨dzdℓ,N⟩=d ddℓ.\left\langle \frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell}, \mathbf N \right\rangle = \frac{\mathrm{d}\,d}{\mathrm{d}\ell}.

由 ℓ\ell 的取向和无侧滑约束,车辆轨迹相对于 ℓ\ell 的单位切向量就是航向单位向量:

dzdℓ=eθ=cos⁡φ T+sin⁡φ N.\frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell} = \mathbf e_\theta = \cos\varphi\,\mathbf T + \sin\varphi\,\mathbf N.

其法向分量为

⟨dzdℓ,N⟩=sin⁡φ.\left\langle \frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell}, \mathbf N \right\rangle = \sin\varphi.

比较同一个法向分量,得到

d ddℓ=sin⁡φ=sin⁡(θ−θr).\frac{\mathrm{d}\,d}{\mathrm{d}\ell} = \sin\varphi = \sin(\theta-\theta_r).

这一结果直接来自以 ℓ\ell 参数化的位置关系,没有借助时间导数反推。它最终不显含参考路径曲率,因为 dN/ds=−κrT\mathrm{d}\mathbf N/\mathrm{d}s=-\kappa_r\mathbf T 在 N\mathbf N 方向上的投影为零。最近点映射及其管状邻域中的一般正则性背景见 [4, §2.1, Eq. (4), Proposition 2 and Theorem 3]。

参考路径弧长与车辆轨迹有向弧长的局部换算

将 dN/ds=−κrT\mathrm{d}\mathbf N/\mathrm{d}s=-\kappa_r\mathbf T 代入上一小节的位置导数,可得

dzdℓ=(1−dκr)dsdℓT+d ddℓN.\frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell} = (1-d\kappa_r)\frac{\mathrm{d}s}{\mathrm{d}\ell}\mathbf T + \frac{\mathrm{d}\,d}{\mathrm{d}\ell}\mathbf N.

同时,车辆轨迹相对于有向弧长 ℓ\ell 的单位切向量为

dzdℓ=cos⁡φ T+sin⁡φ N.\frac{\mathrm{d}\mathbf z}{\mathrm{d}\ell} = \cos\varphi\,\mathbf T + \sin\varphi\,\mathbf N.

比较切向分量:

(1−dκr)dsdℓ=cos⁡φ.(1-d\kappa_r)\frac{\mathrm{d}s}{\mathrm{d}\ell} = \cos\varphi.

当 1−dκr≠01-d\kappa_r\ne0 且 cos⁡φ≠0\cos\varphi\ne0 时,上式给出 ds/dℓ≠0\mathrm{d}s/\mathrm{d}\ell\ne0,因而 ss 与 ℓ\ell 在局部互为参数。这里“局部互为参数”是指在当前点附近,同一段车辆轨迹既可以用 ss 描述,也可以用 ℓ\ell 描述,并且两种参数之间存在可微的局部逆变换。两种弧长的局部换算关系为

dℓds=1−dκrcos⁡φ.(3)\frac{\mathrm{d}\ell}{\mathrm{d}s} = \frac{1-d\kappa_r}{\cos\varphi}. \tag{3}

相应的逆关系为

dsdℓ=cos⁡φ1−dκr.\frac{\mathrm{d}s}{\mathrm{d}\ell} = \frac{\cos\varphi}{1-d\kappa_r}.

关于式 (3) 中 dd、κr\kappa_r 与 φ\varphi 如何共同影响弧长局部换算关系,以及坐标映射退化和参数化失效所对应的 Frenet 几何构型,参见交互文章《参考路径弧长与车辆轨迹有向弧长的局部换算:交互分析》。

横向偏差方程

在 ds/dℓ≠0\mathrm{d}s/\mathrm{d}\ell\ne0 的区间内,链式法则给出

d ddℓ=d ddsdsdℓ.\frac{\mathrm{d}\,d}{\mathrm{d}\ell} = \frac{\mathrm{d}\,d}{\mathrm{d}s}\frac{\mathrm{d}s}{\mathrm{d}\ell}.

代入 d d/dℓ=sin⁡φ\mathrm{d}\,d/\mathrm{d}\ell=\sin\varphi 与式 (3) 的逆关系:

d dds=sin⁡φcos⁡φ/(1−dκr)=(1−dκr)tan⁡φ.\begin{aligned} \frac{\mathrm{d}\,d}{\mathrm{d}s} &= \frac{\sin\varphi} {\cos\varphi/(1-d\kappa_r)}\\[0.6em] &= (1-d\kappa_r)\tan\varphi. \end{aligned}

关系 d d/dℓ=sin⁡φ\mathrm{d}\,d/\mathrm{d}\ell=\sin\varphi 不显含参考路径曲率,因为它描述每单位车辆轨迹有向弧长产生的法向变化;而 d d/ds\mathrm{d}\,d/\mathrm{d}s 显含 κr\kappa_r,因为改用参考路径弧长坐标作为自变量时必须使用式 (3) 的逆关系完成尺度换算。

航向误差方程

在式 (3) 的成立条件下,ss 与 ℓ\ell 在局部互为参数。沿同一车辆轨迹,将航向误差直接写成

φ(s)=θ(ℓ(s))−θr(s).\varphi(s)=\theta\bigl(\ell(s)\bigr)-\theta_r(s).

关于参考路径弧长 ss 求导,并使用 κv=dθ/dℓ\kappa_v=\mathrm{d}\theta/\mathrm{d}\ell、κr=dθr/ds\kappa_r=\mathrm{d}\theta_r/\mathrm{d}s 以及式 (3):

dφds=dθdℓdℓds−dθrds=κv1−dκrcos⁡φ−κr.\begin{aligned} \frac{\mathrm{d}\varphi}{\mathrm{d}s} &= \frac{\mathrm{d}\theta}{\mathrm{d}\ell}\frac{\mathrm{d}\ell}{\mathrm{d}s}-\frac{\mathrm{d}\theta_r}{\mathrm{d}s}\\[0.65em] &= \kappa_v\frac{1-d\kappa_r}{\cos\varphi}-\kappa_r. \end{aligned}

第一项是车辆航向角每单位参考路径弧长的变化率,第二项 κr=dθr/ds\kappa_r=\mathrm{d}\theta_r/\mathrm{d}s 是参考航向角在同一自变量下的变化率。

二状态参考路径弧长域模型

将横向偏差方程与航向误差方程合并,得到以 ss 为自变量、以 (d,φ)(d,\varphi) 为状态的空间域模型:

d dds=(1−dκr)tan⁡φ,dφds=(1−dκr)κvcos⁡φ−κr.(4)\begin{aligned} \frac{\mathrm{d}\,d}{\mathrm{d}s} &= (1-d\kappa_r)\tan\varphi,\\[0.7em] \frac{\mathrm{d}\varphi}{\mathrm{d}s} &= \frac{(1-d\kappa_r)\kappa_v}{\cos\varphi}-\kappa_r. \end{aligned} \tag{4}

这里 κr=κr(s)\kappa_r=\kappa_r(s) 由参考路径给定,κv\kappa_v 是车辆轨迹曲率输入。式 (4) 完全由弧长参数化的几何关系得到,不依赖时间域模型。除最近点投影唯一且 1−dκr≠01-d\kappa_r\ne0 外,它还要求 cos⁡φ≠0\cos\varphi\ne0;否则 ds/dℓ=0\mathrm{d}s/\mathrm{d}\ell=0,参考路径弧长坐标 ss 在该点无法作为车辆轨迹的局部自变量。

弧长域模型中的倒车

式 (4) 同样适用于倒车,无须针对倒车重新推导。原因是前述弧长域推导使用的是沿车辆航向单位向量取正的有向弧长坐标 ℓ\ell,并且只依赖无侧滑条件下的几何微分关系

dz=eθ dℓ,dθ=κv dℓ.\mathrm{d}\mathbf z=\mathbf e_\theta\,\mathrm{d}\ell, \qquad \mathrm{d}\theta=\kappa_v\,\mathrm{d}\ell.

这两个关系没有要求 dℓ>0\mathrm{d}\ell>0。前进时 dℓ>0\mathrm{d}\ell>0,位移 dz\mathrm{d}\mathbf z 与车辆航向单位向量 eθ\mathbf e_\theta 同向;倒车时 dℓ<0\mathrm{d}\ell<0,位移 dz\mathrm{d}\mathbf z 与 eθ\mathbf e_\theta 反向,但 dz/dℓ=eθ\mathrm{d}\mathbf z/\mathrm{d}\ell=\mathbf e_\theta 仍然成立。因此,前述位置微分的分量比较、弧长换算和链式法则在倒车时均保持不变,最终得到的仍是同一个二状态模型。

式 (3) 的逆关系给出

dsdℓ=cos⁡φ1−dκr.\frac{\mathrm{d}s}{\mathrm{d}\ell}=\frac{\cos\varphi}{1-d\kappa_r}.

因此一般情况下,ds\mathrm{d}s 与 dℓ\mathrm{d}\ell 的符号关系为

sgn⁡(ds)=sgn⁡(dℓ)sgn⁡ ⁣(cos⁡φ1−dκr).\operatorname{sgn}(\mathrm{d}s)= \operatorname{sgn}(\mathrm{d}\ell) \operatorname{sgn}\!\left(\frac{\cos\varphi}{1-d\kappa_r}\right).

在通常的 Frenet 跟踪区域 1−dκr>01-d\kappa_r>0 且 ∣φ∣<π/2|\varphi|<\pi/2 内,ds/dℓ>0\mathrm{d}s/\mathrm{d}\ell>0,所以前进对应 ds>0\mathrm{d}s>0,倒车对应 ds<0\mathrm{d}s<0。在给定同一车辆曲率函数 κv(s)\kappa_v(s) 并从同一状态出发时,前进与倒车由同一空间域方程确定同一条 (d,φ)(d,\varphi) 状态曲线,只是沿 ss 的相反方向经过它,式 (4) 的右端不需要修改。

在上述常用跟踪区域内,并保持相同的状态、初始条件与曲率输入,各量的对应变化如下:

量前进倒车是否修改定义或模型
车辆轨迹有向弧长坐标增量 dℓ\mathrm{d}\elldℓ>0\mathrm{d}\ell>0dℓ<0\mathrm{d}\ell<0不修改定义
车辆位移 dz=eθdℓ\mathrm{d}\mathbf z=\mathbf e_\theta \mathrm{d}\ell与 eθ\mathbf e_\theta 同向与 eθ\mathbf e_\theta 反向几何关系不变
参考路径弧长增量 ds\mathrm{d}sds>0\mathrm{d}s>0ds<0\mathrm{d}s<0由式 (3) 决定
航向增量 dθ=κvdℓ\mathrm{d}\theta=\kappa_v \mathrm{d}\ell按正向弧长累积对相同 κv\kappa_v,增量方向反转κv\kappa_v 不因倒车而人为变号
状态 dd、φ\varphi由模型沿 ss 正向演化由模型沿 ss 反向演化状态定义不变,也不预先取反
参考曲率 κr\kappa_r 与车辆曲率 κv\kappa_v使用既定符号约定使用同一符号约定均不修改
式 (4) 的两个右端原形式原形式完全不变

相应地,在上述常用跟踪区域内使用 ss 直接求解时,前进从较小的 ss 积分到较大的 ss,倒车则从较大的 ss 积分到较小的 ss。在从 sstarts_{\mathrm{start}} 开始且 ss 单调减小的倒车轨迹段上,如果数值实现另外定义递增进度 σ=sstart−s\sigma=s_{\mathrm{start}}-s,则由 d/dσ=−d/ds\mathrm{d}/\mathrm{d}\sigma=-\mathrm{d}/\mathrm{d}s 可知两个状态方程的右端会整体乘以 −1-1;这是更换自变量造成的结果,不是倒车改变了式 (4)。

参考文献

  1. Jan Eilbrecht and Olaf Stursberg,

    Challenges of Trajectory Planning with Integrator Models on Curved Roads

    , IFAC-PapersOnLine, 53(2), 2020, doi: 10.1016/j.ifacol.2020.12.2492, §2.1–2.2, Eqs. (14)–(20). 本文使用该文献核对 Frenet 位置坐标、速度分解及时间域状态方程。

  2. Keisuke Yoneda et al.,

    Trajectory Optimization and State Selection for Urban Automated Driving

    , Artificial Life and Robotics, 23, 474–480, 2018, §3.1, Eq. (1). 本文使用该文献核对车辆位置的法向坐标表示 z=r(s)+dN(s)\mathbf z=\mathbf r(s)+d\mathbf N(s)。

  3. Oskar Ljungqvist,

    Motion Planning and Feedback Control Techniques with Applications to Long Tractor-Trailer Vehicles

    , Linköping Studies in Science and Technology, Dissertation No. 2070, Linköping University, 2020, doi: 10.3384/diss.diva-165246, §4.3. 本文使用该文献核对 Frenet 坐标正则域及三状态车辆运动学模型。

  4. Catherine Kublik and Richard Tsai,

    Integration over Curves and Surfaces Defined by the Closest Point Mapping

    , Research in the Mathematical Sciences, 3, Article 3, 2016, doi: 10.1186/s40687-016-0053-1, §2.1, Eq. (4), Proposition 2 and Theorem 3. 本文使用该文献说明最近点映射及其管状邻域中的一般正则性背景。