基于满足动力学约束下的路径规划(State Lattice Planning,Kinodynamic-RRT*)

总结课程《深蓝学院移动机器人路径规划》

Kinodynamic : Kinematic + Dynamic

1.State Lattice Planning

汽车质点模型不满足动力学约束,需要满足两点之间是可行的运动连接。
1.离散机器人控制空间
在这里插入图片描述

2.离散机器人状态空间
在这里插入图片描述
考虑动力学约束后,建立一个机器人模型,由状态和输入描述。 S ˙ = f ( s , u ) \dot{S}=f(s,u) S˙=f(s,u)
当给系统不同的输入,持续一定的时间T。(离散机器人控制空间)
在这里插入图片描述(没有任务导向性,效率相对低。)
当选择两个状态,找到两个状态之间的输入控制(u,T)。(离散机器人状态空间)
在这里插入图片描述(导向性,难以实现)


Sample in control space
在这里插入图片描述在这里插入图片描述
当确定末状态后,解算中间状态。状态转移方程:
在这里插入图片描述


2.1 Boundary Value Problem

在这里插入图片描述

2.2 Optimal Boundary Value Problem (OBVP)

从一个状态到另一个状态,使得其 J k J_k Jk最小值。
J ∑ = ∑ k = 1 3 J k , J k = 1 T ∫ 0 T j k ( t ) 2 d t J_{\sum}=\sum_{k=1}^3{J_k},J_k={1 \over T}\begin{array}{cc} \int_0^T{ j_k(t)^2 dt }\end{array} J=k=13Jk,Jk=T10Tjk(t)2dt状态: S k = ( p k , v k , a k ) S_k=(p_k,v_k,a_k) Sk=(pk,vk,ak) ,输入 u k = j k u_k=j_k uk=jk
系统模型: S k ˙ = f s ( s k , u k ) = ( v k , a k , j k ) \dot{S_k}=f_s(s_k,u_k)=(v_k,a_k,j_k) Sk˙=fs(sk,uk)=(vk,ak,jk)

利用极小值原理:
在这里插入图片描述
引入协态costate: λ = ( λ 1 , λ 2 , λ 3 ) \lambda=(\lambda_1,\lambda_2,\lambda_3) λ=(λ1,λ2,λ3),由系统变量决定。
定义Hamiltonian function: H ( s , u , λ ) = 1 T j 2 + λ T f s ( s , u ) H(s,u,\lambda)={1\over T}j^2+\lambda^Tf_s(s,u) H(s,u,λ)=T1j2+λTfs(s,u)得到最优的状态量 s s s,最优输入 u u u

3.Planning in Frenet-serret Frame

在这里插入图片描述

4.混合A* (Hybrid A*)

在这里插入图片描述
在每个栅格中保证只存在一个节点。
启发函数的设计:
在这里插入图片描述

5.基于运动学约束的RRT*

1.(Sample)在整个状态空间采样,包括位置和速度。
2.(Near)利用optimal control将新找到的节点和邻近节点进行OBVP连接。
3.(ChooseParent)以最优控制cost指标为半径的球,范围内的节点。
4.(find near nodes efficiently)后向可达集
在这里插入图片描述
5.(Rewire)前向可达集

这是一道比较经典的计数问题。题目描述如下: 给定一个 $n \times n$ 的网格图,其中一些格子被标记为障碍。一个连通块是指一些被标记为障碍的格子的集合,满足这些格子在网格图中连通。一个格子是连通的当且仅当它与另一个被标记为障碍的格子在网格图中有公共边。 现在,你需要计算在这个网格图中,有多少个不同的连通块,满足这个连通块的大小(即包含的格子数)恰好为 $k$。 这是一道比较经典的计数问题,一般可以通过计算生成函数的方法来解决。具体来说,我们可以定义一个生成函数 $F(x)$,其中 $[x^k]F(x)$ 表示大小为 $k$ 的连通块的个数。那么,我们可以考虑如何计算这个生成函数。 对于一个大小为 $k$ 的连通块,我们可以考虑它的形状。具体来说,我们可以考虑以该连通块的最左边、最上边的格子为起点,从上到下、从左到右遍历该连通块,把每个格子在该连通块中的相对位置记录下来。由于该连通块的大小为 $k$,因此这些相对位置一定是 $(x,y) \in [0,n-1]^2$ 中的 $k$ 个不同点。 现在,我们需要考虑如何计算这些点对应的连通块是否合法。具体来说,我们可以考虑从左到右、从上到下依次处理这些点,对于每个点 $(x,y)$,我们需要考虑它是否能够与左边的点和上边的点连通。具体来说,如果 $(x-1,y)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们就是连通的;同样,如果 $(x,y-1)$ 和 $(x,y)$ 都在该连通块中且它们在网格图中有公共边,那么它们也是连通的。如果 $(x,y)$ 与左边和上边的点都不连通,那么说明这个点不属于该连通块。 考虑到每个点最多只有两个方向需要检查,因此时间复杂度为 $O(n^2 k)$。不过,我们可以使用类似于矩阵乘法的思想,将这个过程优化到 $O(k^3)$ 的时间复杂度。 具体来说,我们可以设 $f_{i,j,k}$ 表示状态 $(i,j)$ 所代表的点在连通块中,且连通块的大小为 $k$ 的方案数。显然,对于一个合法的 $(i,j,k)$,我们可以考虑 $(i-1,j,k-1)$ 和 $(i,j-1,k-1)$ 这两个状态,然后把点 $(i,j)$ 加入到它们所代表的连通块中。因此,我们可以设计一个 $O(k^3)$ 的 DP 状态转移,计算 $f_{i,j,k}$。 具体来说,我们可以考虑枚举连通块所包含的最右边和最下边的格子的坐标 $(x,y)$,然后计算 $f_{x,y,k}$。对于一个合法的 $(x,y,k)$,我们可以考虑将 $(x,y)$ 所代表的点加入到 $(x-1,y,k-1)$ 和 $(x,y-1,k-1)$ 所代表的连通块中。不过,这里需要注意一个细节:如果 $(x-1,y)$ 和 $(x,y)$ 在网格图中没有相邻边,那么它们不能算作连通的。因此,我们需要特判这个情况。 最终,$f_{n,n,k}$ 就是大小为 $k$ 的连通块的个数,时间复杂度为 $O(n^2 k + k^3)$。 参考代码:
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

秃头队长

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值