PX4二次开发——多旋翼姿态控制

本文详细介绍了PX4飞行控制系统中的串级PID加前馈复合控制策略,包括内外环控制逻辑、罗德里格旋转公式及姿态解耦控制原理等关键技术点。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

一、准备知识

1.串级PID

        PX4中多旋翼姿态控制是采用的串级PID加前馈的复合控制策略。分内外两环,外环只有一个比例环节,内环采用PID加前馈的控制方式。外环主要是由control_attitude()函数负责,由期望姿态(使用四元数表示)和当前姿态(也是使用四元数表示)得到角度差,作用于外环的比例环节(对偏航的处理比较特殊),得到期望角速率,传递给内环;内环主要是由control_attitude_rates()函数负责,由期望角速率和当前角速率得到角速率差,作用于内环的PID环节,同时期望角速率通过前馈通道累加到PID输出的控制量上,得到最终的要传递给混控器的控制量(可以理解为转矩)。此时多旋翼姿态控制模块的使命就基本完成了,接下来(src/drivers/px4fmu/fmu.cpp)由混控器根据这些控制量和混控器矩阵(仅对旋翼而言)进行力矩的分配,然后将分配到各个电机上的力矩换算成PWM值,由px4fmu最终落实到硬件上。

        PX4中的串级PID就讲这么多,深入了解请参考其他相关资料。

2.罗德里格旋转公式


3.Roll/Pitch和Yaw解耦控制

首先,简单通俗地介绍一下多旋翼姿态控制的基本原理(以四旋翼为例)。

四旋翼横滚角和俯仰角控制原理:

        电机转动产生轴向的升力,对称位置的电机转速不一样,产生升力差,进而产生改变横滚角或俯仰角的转矩。

四旋翼偏航角控制原理:

        对称位置的电机转速一样,相邻位置的电机转速不一样,而多旋翼对称位置电机转向一致,相邻位置电机转向相反,这样会产生垂直于轴向的扭矩,进而改变偏航角。

现在,我们进入正题。

为什么要解耦控制?

        因为横滚角/俯仰角和偏航角控制方式不一样,横滚角/俯仰角是考靠轴向升力之差产生转矩,而偏航角是靠扭矩,所以它们的响应速度也不一样,偏航角控制的响应速度只有前两者的一半左右,所以分开控制能获得更快的响应速度。

怎么实现解耦控制?

        先【主要】Roll和Pitch联动,使当前B系中的Z轴向期望姿态的Z轴靠拢,后【主要】进行Yaw的旋转,即对齐X和Y轴。

        实际代码中的实现:给Yaw的旋转加一个范围是[0,1]的权重,它是取的(cos<z,z_sp>)^2(两Z轴之间夹角的余弦值的平方),使得Z轴夹角越大,Yaw旋转的权重越小,反之越大。

4.“轴角法”表示误差

        当前姿态可以通过绕一条轴旋转一定的角度来达到期望姿态。用一条旋转轴u(单位向量)和一个旋转角θ来表示旋转的方法就是所谓的“轴角法”,可以表示为方向为u,模为θ的向量  θu = (θ a, θ b, θ c)T。这样乘外环比例系数后直接就可以得到期望机体角速率。(代码中对Roll和Pitch误差采用了这种方法来表示,对Yaw误差则是单独处理,然后附加在误差向量上,差别不大)

二、多旋翼姿态控制模块(mc_att_control)源码剖析(基于1.7.2版本代码)

        前面提到控制分了内外两环,外环由control_attitude()函数负责,内环由control_attitude_rates()函数负责,下面就开始剖析这两个函数。

首先外环:




然后内环:



三、最后

        PX4多旋翼位置控制就先讲到这里,有没有看得头晕眼花啊?偷笑(去年暑假我第一次接触无人机的时候反正是晕了微笑

        这篇文章花了本考研党接近一天时间,所有的文字都是自己动手敲的,所有的注释也都是自己亲手加上去的,如果有问题还请多多包涵,另外尊重原创大笑
        (上传之后图变得不是很清楚了,单击图片可以看大图,需要原图的话发邮件给我,我把图给发过去。e-mail:hello.icq@qq.com)
### PX4 旋翼无人机姿态控制程序解析 #### 姿态控制器架构概述 PX4旋翼姿态控制系统是一个分层设计的闭环反馈系统,主要包括内环角度速率控制和外环角度控制两部分。该系统通过接收来自外部传感器的数据并结合设定点信息来调整电机输出,实现稳定飞行。 #### 内环角速度控制 为了精确地维持期望的姿态,在内部存在一个快速响应的角度速率回路。此环节负责将实际测量得到的身体坐标系下的角速度转换为目标指令所对应的参考值之间的差异,并利用PID算法计算出所需的转矩命令: ```cpp // 计算角速度误差 Vector3f rate_error = _rate_setpoint - gyro; // 应用 PID 控制律获得力矩需求 _vector_thrust_body += _ctrl_rate.getControl(rate_error, dt); ``` 上述代码片段展示了如何基于当前陀螺仪读数`gyro`与预期角速度设定点 `_rate_setpoint` 来求得偏差 `rate_error` ,进而调用 PID 调节器对象 `_ctrl_rate` 获取最终施加于机体上的额外推力矢量[^1]。 #### 外环姿态控制 在外围则构建了一个较慢但更注重精度的态度保持机制。它接受由导航模块产生的航向、俯仰及横滚角作为输入信号,并将其转化为相应的角速度设定点传递给内环处理: ```cpp // 更新姿态估计 _quat.setFromEuler(roll_sp, pitch_sp, yaw_sp); // 将四元数形式的姿态设为新的状态机目标 _attitude_setpoint.q_d = _quat; // 使用四元数差乘法获取旋转矩阵中的前两个轴方向 _matrix_R_to_earth = Matrix3f(_attitude_setpoint.q_f).transpose(); _matrix_R_to_body = Matrix3f(_quat).transpose(); // 得到相对于水平面投影后的期望朝向变化率 Vector2f body_z_axis_proj = rpy_from_rotation(_matrix_R_to_body)[2].xy().normalize(); Vector2f earth_z_axis_proj = rpy_from_rotation(_matrix_R_to_earth)[2].xy().normalize(); // 设置横向偏移补偿项 float lateral_offset_compensation = atan2(-body_z_axis_proj.y(), body_z_axis_proj.x()) - atan2(-earth_z_axis_proj.y(), earth_z_axis_proj.x()); // 构建完整的三维空间内的角速度请求 Vector3f att_rate_setpoint(lateral_offset_compensation * Kp_roll, (pitch_angle_desired - attitude_estimate.pitch) * Kp_pitch, heading_rate_command); ``` 这里实现了从欧拉角表示的姿态目标至角速度域的具体映射过程;其中涉及到四元数运算用于表征刚体转动关系,同时也考虑到了侧风等因素可能引起的附加扰动修正[^3]。 #### 推力分配策略 当完成内外双闭环调节之后,下一步就是依据所得总合力矩大小合理安排各执行机构的工作强度。对于常规布局类型的旋翼平台而言,通常会采用如下方式来进行动力学解耦和平滑过渡: ```cpp for (int i = 0; i < NUM_ACTUATORS; ++i){ float thrust_i = base_thrust + mix_factor[i][ROLL]*moment_x + mix_factor[i][PITCH]*moment_y+ mix_factor[i[YAW]]*moment_z; // 防止出现异常情况导致失控现象的发生 if (!isfinite(thrust_i)){ thrust_i = 0.f; } } ``` 这段伪码描述了按照预定义好的混合系数矩阵 `mix_factor[][]` 对整体作用效果按需拆分成单个螺旋桨应承担的任务份额,并且加入了必要的数值检验措施以保障安全性[^4]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值