掌握A*搜索算法原理与实现,理解启发函数与3D扩展。
f(n) = g(n) + h(n)
g=实际代价 h=启发估计 f=总估计代价
| 算法 | 完备 | 最优 | 效率 |
|---|---|---|---|
| BFS | ✅ | ✅(无权) | 低 |
| Dijkstra | ✅ | ✅ | 中 |
| A* | ✅ | ✅(h可容许) | 高 |
| 贪心 | ✅ | ❌ | 最高 |
可容许性:h(n) ≤ h*(n)保证最优性
| 启发 | 公式 | 场景 |
|---|---|---|
| 曼哈顿 | |Δx|+|Δy| | 4方向 |
| 欧几里得 | √(Δx²+Δy²) | 任意方向 |
| 对角距离 | max+0.414·min | 8方向 |
=== A*路径规划结果 ===
路径长度: 24步
探索节点数: 155
地图(S=起点 G=终点 ●=路径 █=障碍):
· · · · · · · · · · · · · · · · · · · ·
· S · · · · · · · · · · · · · · · · · ·
· · ● ● · · · · · · · · · · · · · · · ·
· · · · ● · · · · · · · · · · · · · · ·
· · · · · ● ● ● ● · · · · · · · · · · ·
· · · · · · · · █ ● · · · · · · · · · ·
· · · · · · · · █ · ● · · · · · · · · ·
· · · · · · · · █ · · ● · · · · · · · ·
· · · █ █ █ █ █ █ █ █ █ ● · · · · · · ·
· · · · · · · · █ · · · · ● · · · · · ·
· · · · · · · · █ · · · · · ● · · · · ·
· · · · · · · · █ · · · · · · ● ● ● · ·
· · · · · · · · █ · █ █ █ █ █ █ █ █ ● ·
· · · · · · · · █ · · · · · · · · · ● ·
· · · · · · · · █ · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · ● ·
· · · · · · · · · · · · · · · · · · G ·
· · · · · · · · · · · · · · · · · · · ·
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²)✅ 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 |
| 控制频率 | 400Hz | 1kHz | 1kHz+ |
| 传感器融合 | 互补滤波 | EKF | EKF/VIO |
| 续航时间 | 20-30min | 30-45min | 视载荷而定 |
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 |
| 串级PID | O(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 # 协方差下降
| 特性 | PX4 | ArduPilot | Betaflight |
|---|---|---|---|
| 定位 | 研究/工业 | 工业/爱好者 | 竞速穿越 |
| 代码量 | ~500K行 | ~800K行 | ~200K行 |
| 支持的机型 | 多旋翼/固定翼/VTOL | 多旋翼/固定翼/直升机/车/船 | 多旋翼/固定翼 |
| 导航能力 | 强(全面) | 强(最全面) | 弱(仅自稳) |
| 实时性 | NuttX RTOS | ChibiOS/RTOS | Bare-metal |
| 仿真支持 | SITL/Gazebo/AirSim | SITL/Gazebo | 有限 |
| 社区活跃度 | 高 | 最高 | 高(竞速圈) |
无人机系统属于安全关键(Safety-Critical)系统,设计时必须遵循以下原则:
本课所学技术在以下场景中直接应用:
将本课模块集成到完整系统时的注意事项: