1. RRT算法快速探索随机树的诞生与局限我第一次接触RRT算法是在2013年做移动机器人项目时。当时需要让机器人在未知环境中自主导航A*这类基于网格的算法在动态环境中显得力不从心直到发现了这个基于随机采样的神奇方法。RRT的核心思想就像盲人摸象你站在起点位置随机朝某个方向伸出手采样找到离你最近的树枝最近邻节点然后朝那个方向生长一小段步长扩展。重复这个过程直到摸到终点。这种生长式探索特别适合高维空间和复杂环境因为它避免了构建完整地图的开销。用Python实现基础RRT只需要不到50行代码。关键步骤包括def rrt(start, goal, map_area, max_iter1000, step_size20): tree {start: None} # 用字典存储树结构 for _ in range(max_iter): rand_point random_sample(map_area) nearest find_nearest(tree.keys(), rand_point) new_point steer(nearest, rand_point, step_size) if not collision_check(nearest, new_point, obstacles): tree[new_point] nearest if distance(new_point, goal) step_size: if not collision_check(new_point, goal, obstacles): return reconstruct_path(tree, goal) return None但实际使用中我发现了三个致命问题路径质量不稳定同样的环境运行10次可能得到10条完全不同长度的路径收敛速度慢在开阔区域会浪费大量采样点在无用区域锯齿现象由于随机扩展路径常常像醉汉走路一样曲折最让我印象深刻的是在一次无人机测试中RRT生成的路径竟然让无人机在空中画了个之字形多飞了将近30%的距离。这促使我开始研究它的改进算法。2. RRT*渐进最优的突破与代价2015年实验室新来的师弟兴奋地告诉我师兄RRT能找最优路径了这个星号()带来的改变确实令人惊艳。RRT*在两方面做了关键改进重布线机制就像城市道路规划不仅考虑新建道路还会检查是否能让周边居民抄近路。具体实现时每个新节点会检查半径r范围内的邻居看看能否通过这个新节点获得更短路径def rewire(tree, new_node, neighbors, obstacles): for neighbor in neighbors: potential_cost tree[new_node][cost] distance(new_node, neighbor) if potential_cost tree[neighbor][cost]: if not collision_check(new_node, neighbor, obstacles): tree[neighbor][parent] new_node update_cost(tree, neighbor, potential_cost)代价函数传播每个节点都记录从起点到该点的路径代价这让我想起Dijkstra算法。在二维平面中我们通常用欧式距离作为代价而在机器人领域可能需要考虑能耗、通过性等因素。但实际部署时发现了新问题在30x30m的仓库环境中要让路径代价收敛到理论最优值的5%以内需要近5000次迭代耗时超过3秒。更糟的是当空间维度增加时如机械臂的6维空间收敛速度呈指数级下降。有次演示时机械臂足足思考了8分钟才开始移动场面一度尴尬。3. Informed-RRT*椭圆采样的智慧飞跃2017年读到的Informed-RRT*论文让我眼前一亮。它巧妙地利用了椭圆的基本性质椭圆上任意点到两个焦点的距离之和相等。将起点和终点作为焦点当前最优路径长度决定椭圆大小就形成了一个天然的优化区域。具体实现时需要在找到初始路径后切换采样策略def informed_sample(c_best, start, goal, map_area): if c_best float(inf): return random_sample(map_area) # 初始阶段仍用全局采样 # 构建旋转矩阵将椭圆对齐起点和终点 c_min distance(start, goal) center (start goal) / 2 angle math.atan2(goal[1]-start[1], goal[0]-start[0]) # 椭圆参数计算 a c_best / 2 b math.sqrt(a**2 - (c_min/2)**2) # 在单位圆内采样后映射到椭圆 while True: r np.random.uniform(-1, 1, 2) if np.linalg.norm(r) 1: x a * r[0] * math.cos(angle) - b * r[1] * math.sin(angle) center[0] y a * r[0] * math.sin(angle) b * r[1] * math.cos(angle) center[1] if 0 x map_area[0] and 0 y map_area[1]: return (x, y)实测效果令人振奋在相同迭代次数下路径优化速度比RRT快3-5倍。特别是在狭窄通道环境中传统RRT容易卡在局部最优而Informed版本能快速调整采样策略。不过要注意椭圆采样时的坐标变换我有次忘记考虑旋转角度结果采样点全部偏移到了地图外。4. 工程实践中的调参经验与陷阱经过数十个项目的实战检验我总结出以下关键参数设置原则步长选择室内环境地图对角线长度的2-5%无人机最大转弯半径的1.5倍机械臂关节运动限度的10-20%邻域半径公式def calculate_radius(dim, cardV, gamma1.5, eta20): return min(gamma * (math.log(cardV)/cardV)**(1/dim), eta)其中dim为空间维度cardV是当前节点数。gamma过大导致计算负担重过小则可能错过优化机会。常见坑点碰撞检测不精确特别是对于非点状机器人需要计算包络体积动态障碍物处理简单的重新规划可能导致抖动浮点误差累积长时间运行后可能出现节点无法连接的情况有次为AGV小车部署时没考虑车体半径结果规划出的路径让小车卡在了货架之间。后来改用膨胀障碍物的方法才解决问题。另一个教训是在ROS中实现时没有做适当的坐标系转换导致所有采样点都偏移了2米。在算法选择上我的经验法则是只需快速可行解用RRT或RRT-Connect静态环境求优解Informed-RRT*高维空间考虑降低采样维度或使用投影方法动态环境结合人工势场或速度障碍法最近在处理一个机械臂拣选项目时我将Informed-RRT*与深度学习结合用神经网络预测优化椭圆的朝向和偏心率使收敛速度又提升了40%。这让我意识到经典算法与现代AI技术的结合还有巨大探索空间。