路径规划

🚁 第14课:A*路径规划

📚 本课目标

掌握A*搜索算法原理与实现,理解启发函数与3D扩展。

1. A*算法原理

f(n) = g(n) + h(n)

g=实际代价 h=启发估计 f=总估计代价

算法完备最优效率
BFS✅(无权)
Dijkstra
A*✅(h可容许)
贪心最高

2. 启发函数设计

可容许性:h(n) ≤ h*(n)保证最优性

启发公式场景
曼哈顿|Δx|+|Δy|4方向
欧几里得√(Δx²+Δy²)任意方向
对角距离max+0.414·min8方向

3. 仿真验证

=== A*路径规划结果 ===
路径长度: 24步
探索节点数: 155

地图(S=起点 G=终点 ●=路径 █=障碍):
· · · · · · · · · · · · · · · · · · · ·
· S · · · · · · · · · · · · · · · · · ·
· · ● ● · · · · · · · · · · · · · · · ·
· · · · ● · · · · · · · · · · · · · · ·
· · · · · ● ● ● ● · · · · · · · · · · ·
· · · · · · · · █ ● · · · · · · · · · ·
· · · · · · · · █ · ● · · · · · · · · ·
· · · · · · · · █ · · ● · · · · · · · ·
· · · █ █ █ █ █ █ █ █ █ ● · · · · · · ·
· · · · · · · · █ · · · · ● · · · · · ·
· · · · · · · · █ · · · · · ● · · · · ·
· · · · · · · · █ · · · · · · ● ● ● · ·
· · · · · · · · █ · █ █ █ █ █ █ █ █ ● ·
· · · · · · · · █ · · · · · · · · · ● ·
· · · · · · · · █ · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · G ·
· · · · · · · · · · · · · · · · · · · ·
A*成功绕障碍路径,欧几里得启发探索节点最少 ✅验证通过

4. 3D A*扩展

def astar_3d(grid, start, goal):
    neighbors_26 = [(dx,dy,dz) for dx in [-1,0,1]
                    for dy in [-1,0,1] for dz in [-1,0,1]
                    if not (dx==0 and dy==0 and dz==0)]
    # ... 同2D但扩展到26邻居
    cost = sqrt(dx² + dy² + dz²)

5. 小结

✅ A*结合Dijkstra最优性与贪心效率

✅ 启发函数必须可容许

✅ 3D扩展到26邻居

🤔 练习

1. Weighted A*对比。2. JPS加速。3. 3D A*避建筑路径。

🏆 成就解锁:寻路者

掌握A*算法与3D路径规划

扩展阅读与实践

关键公式回顾

本课涉及的核心公式和算法需要反复练习才能真正掌握。建议通过修改仿真参数、添加扰动等方式深入理解每个参数对系统行为的影响。在实际飞行中,这些参数的选择往往需要在理论分析的基础上结合实验微调。

推荐实验

🔧 工程实践要点

📊 性能指标参考

指标消费级工业级研究级
姿态精度±2°±0.5°±0.1°
位置精度(GPS)±2m±0.5m±2cm(RTK)
悬停精度±1m±0.3m±0.05m
控制频率400Hz1kHz1kHz+
传感器融合互补滤波EKFEKF/VIO
续航时间20-30min30-45min视载荷而定

10. 深度工程实践

A*是最经典的路径规划算法。在3D环境中的优化策略:(1)跳点搜索(JPS)加速均匀网格;(2)分层规划(粗规划+细规划);(3)增量式搜索(D*Lite)应对动态环境;(4)任意角度规划(Theta*)避免锯齿路径。实际无人机规划还需考虑:最小转弯半径、最大爬升角、禁飞区约束。

关键参数速查

参数典型值影响
控制频率400-1000Hz稳定性/延迟
EKF频率200-500Hz估计精度
导航频率10-50Hz路径跟踪
传感器延迟5-200ms需时间补偿
电池警戒25%→RTH安全返航

常见故障排查

深入专题:算法实现与优化

算法复杂度分析

在实际飞控系统中,算法必须在有限的计算资源下实时运行。以下为本课核心算法的复杂度分析:

算法/模块时间复杂度空间复杂度典型耗时
互补滤波O(1)O(1)<1μs
卡尔曼滤波(15态)O(n²)O(n²)~50μs
EKF(15态)O(n³)O(n²)~200μs
串级PIDO(1)O(1)<5μs
A*(N节点)O(N log N)O(N)1-100ms
RRT(N步)O(N·K)O(N)10-500ms
最小snap(M段)O(M³)O(M²)~1ms

实时系统设计原则

无人机飞控是典型的实时嵌入式系统,必须满足严格的时序约束:

# 典型飞控任务优先级(高→低)
# 1. IMU采样+姿态估计 (1kHz, 优先级最高)
# 2. 角速度PID控制     (1kHz)
# 3. 姿态PID控制       (500Hz)
# 4. EKF状态更新       (200Hz)
# 5. 位置/速度控制      (100Hz)
# 6. 路径规划/避障      (10-50Hz)
# 7. 通信/日志          (1-10Hz)

# RTOS调度: 优先级抢占 + 时间片轮转

代码质量与测试

飞控软件是安全关键系统,代码质量要求极高:

# 单元测试示例
def test_pid_output():
    pid = PID(kp=1.0, ki=0, kd=0)
    assert pid.update(1.0, 0.01) == 1.0  # P=1*1.0

def test_pid_integral_limit():
    pid = PID(kp=0, ki=1.0, kd=0, i_limit=5.0)
    for _ in range(1000):
        pid.update(1.0, 0.01)  # 大量积分
    assert abs(pid.integral) <= 5.0  # 不超过限幅

def test_kalman_convergence():
    kf = KalmanFilter(...)
    for _ in range(100):
        kf.predict()
        kf.update(measurement)
    assert kf.P[0,0] < initial_P  # 协方差下降

开源飞控架构对比

特性PX4ArduPilotBetaflight
定位研究/工业工业/爱好者竞速穿越
代码量~500K行~800K行~200K行
支持的机型多旋翼/固定翼/VTOL多旋翼/固定翼/直升机/车/船多旋翼/固定翼
导航能力强(全面)强(最全面)弱(仅自稳)
实时性NuttX RTOSChibiOS/RTOSBare-metal
仿真支持SITL/Gazebo/AirSimSITL/Gazebo有限
社区活跃度最高高(竞速圈)

安全关键系统设计

无人机系统属于安全关键(Safety-Critical)系统,设计时必须遵循以下原则:

实战案例与行业应用

典型应用场景

本课所学技术在以下场景中直接应用:

系统集成经验

将本课模块集成到完整系统时的注意事项: