由于一直在调试本项目,好久没有发文章,最近本项目的pid调试初见成效!开始正文前首先感谢各位粉丝的支持,以及对本项目技术上支持的老师以及师兄,谢谢你们!
基于stm32的多旋翼无人机
一、多旋翼无人机飞行原理
四轴飞行器基本原理是通过飞控控制四个电机旋转带动浆叶产生升力,分别控制每一个电机和浆叶产生不同的升力从而控制飞行器的姿态和位置。四轴在空中可以实现八种运动,分别为垂直上升、垂直下降、向前运动、向后运动、向左运动、向后运动、顺时针改变航向、逆时针改变航向。
1.1 多旋翼无人机的组成
四轴飞行器主要是由电机、电调、电池、浆叶、机架、遥控器、飞控组成。下面以我们四轴及市场上常见的 diy 大四轴来介绍这些部件。
1.1.1 电机
电机根据目前市场的供给以及本项目的需求,采用的是空心杯无刷电机8520,如下图:
其转速可达到12000~15000转/min,并且其价格较低,性能较稳定,更适合本项目的使用。
1.1.2 电调
电调即为电子调速器,控制电机转动、停止及转速。有刷电机电调通常只需要一个 mos 管,飞控输出 pwm 即可控制电机,电调所采用的mos管如下图所示:无刷电机电调模块内部通常由一个 mcu 和三相桥电路组成,mcu 通过控制三相桥来实现无刷电机换相。同样,无刷电机电调模块也只需飞控输出 pwm 即可控制电机
1.1.3 航模电池
航模所用电池为3.7v锂电池(可充电),由于受到多旋翼无人机自身动力的问题,故采用800mh锂电池来供电,遥控器也采用800mh供电,飞行时间大概20min左右。电池如下图:
电池插座采用空对空接头。
电池容量 mah:锂电池的容量,如 2000mah 的电池,以 2000ma 放电,可持续放电 1
小时;以 1000ma 放电,可持续放电 2 小时。
电池节数:电池 2s、3s、4s 代表锂电池节数。锂电池 1 节的标准电压为 3.7v,3s 代
表有 3 节 3.7v 的电池在里面,电压为 11.1v。
1.1.4 正反浆
浆液采用的是75mm直径的正反浆,如下图所示:
浆叶旋转时会产生自旋力导致四轴自旋,为了抵消自旋力相隔电机的浆叶旋转方向要不一样,但是浆的风都是要往下吹,这就出现了正反浆的说法。通常顺时针转的叫正浆,逆时针转的是反浆。
1.2 垂直上升及下降
当四轴飞行在空中自稳后,m1、m2、m3、m4 四个电机同时转速增大或同时转速减小,即可发生垂直上升运动或垂直下降运动,如下图所示:
1.3 向前飞行及向后飞行
当四轴飞行在空中自稳后,m2、m3 转速增大 m1、m4 转速不变或减小即可实现向前运动。相反,m2、m3 转速减小或不变 m1、m4 转速增加,即可实现向后运动,如下图所示:
向左和向右飞行同理可以实现,这里就不做介绍。
1.4 顺时针改变航向和逆时针改变航向
当四轴飞行在空中自稳后,m1、m3 转速增大 m2、m4 转速不变或减小即可实现顺时针改变航向。m1、m3 转速减小或不变 m2、m4 转速增加即可实现逆时针改变航向。如下图所示:
二、多旋翼无人机飞控电路设计
2.1 系统框架设计
根据minfly的系统框架图,来进行设计,主控芯片采用stm32103c8t6,三轴陀螺仪采用mpu6050芯片,收发无线模块采用nrfl2401模块(亿百特)等。
2.2 主控mcu模块
主控采用的是stm32103c8t6芯片,电路图如下:
主控 mcu 为四轴飞行器的大脑,对飞行器稳定飞行起着至关重要的作用。它同时承担
着多种责任,包括:传感器数据读取、数据融合、pid 控制、电机控制、无线通信和 usb
通信等。
主控 mcu 连接了一个 usb 接口,此接口可以用作与上位机通信,也可以用作固件升级,是一个非常方便适用的接口。
2.3 三轴加速陀螺仪模块
三轴加速陀螺仪模块采用的是mpu6050芯片作为主控,mpu6050 imu 在单芯片上集成了一个 3 轴加速度计和一个 3 轴陀螺仪。以及一个可扩展的数字运动处理器 dmp(digital motion processor)。它也被称为六轴运动跟踪设备或 6 dof 设备,因为它有 6 个输出,即 3 个加速度计输出和 3 个陀螺仪输出。以当前地面为水平面检测x,y,z三轴方向上的加速度,并转换为电信号来进行输出。同时进行三轴陀螺仪传感器进行使用,得到三轴方向上的倾角,即姿态角。如下图所示:
mpu6050主控芯片的控制原理如图:
设计电路图如下:
2.4 电机驱动模块
本项目采用微型高速 8520空心杯电机,电机转速高达 15000r/min,能够为飞行器提供充沛的动力。电机采用 nmos 管 si2302,3v 门级驱动电压下,导通电阻只有几十毫欧,驱动电流高达 3a,轻松驱动 8520空心杯电机,从而带动飞行器飞行。电路设计如下图:
2.5 无线通信模块
无线通信模块所采用的是nrf24l01模块来进行通信,nrf24l01是一款新型单片射频收发器件,工作与2.4ghz~2.5ghz ism频段。内置频率合成器、功率放大器、晶体振荡器、调制器等功能模块,并融合了增强型shockburst技术,其中输出功率和通信频道可通过程序进行配置。
nrf24l01有工作模式有四种:
- 收发模式
- 配置模式
- 空闲模式
- 关机模式
如图所示(e01-ml01s):
芯片方案:nrf24l01p
工作频率:2.4~2.525ghz
发射功率:0dbm
通信距离:0.1km
接口类型:spi
产品重量:0.5±0.1g
产品简介:采用挪威nordic公司原装进口nrf24l01p射频芯片,收发 一体;全进口工业级元器件,邮票孔贴片型,内置pcb板载天线;体积小易嵌入,性能优异,数据传输快,适合近距离传输。
要实现通信还需要e01-ml01d模块来实现交互通信,模块的图片如下:
芯片方案:nrf24l01p
工作频率:2.4~2.525ghz
发射功率:0dbm
通信距离:0.1km
通信接口:spi
产品重量:0.9±0.1g
产品简介:挪威进口nrf24l01+芯片,日本进口阻容感,美国进口晶振,无铅焊接工艺; 目前已经大量用于电力线,品质相当可靠,使用寿命长。
无线通信模块电路图如下:(分别为接受端和发送端)
2.6 升压、稳压电路模块
考虑本项目所需8520电机供电为3.7v,而本项目还需要5v供电,电池的供电为3.7v,因此需要设计升压电路来进行从3.7v到5v的升压,根据模拟电子技术以及电路原理自主设计了3.7v----5v的升压电路。主控为dc–dc电源芯片,采用肖特基二极管来进行外部稳压,具体电路如下图所示:(篇幅影响不进行详细的讲解,感兴趣的粉们可以私信我讨论相关知识)
最终得到5v的升压电源,方便项目其他模块的供电。
相反可以设计出稳压电路从5v—3.3v的降压电路,具体原理不做说明,具体见下图所示:
2.7 led电路模块
此部分为led指示灯电路设计,主要由3个led直插式灯以及两枚rgb贴片灯组成,具体i单路如下图所示:
后续逻辑功能将在程序中进行设计。
三、多旋翼无人机飞控算法设计
3.1 飞行姿态表示法
飞行器姿态有多种表示方式,常见的是四元数,欧拉角,矩阵和轴角。他们各自有其自身的优点,在不同的领域使用不同的表示方式。在四轴飞行器中使用到了四元数和欧拉角。
3.1.1 欧拉角
用来确定定点转动刚体位置的 3 个一组独立角参量,由章动角 θ、旋进角(即进动角)ψ 和自转角 j 组成,为莱昂哈德·欧拉首先提出而得名。对于在三维空间里的一个参考系,任何坐标系的取向,都可以用三个欧拉角来表现。参考系又称为实验室参考系,是静止不动的。而坐标系则固定于刚体,随着刚体的旋转而旋转。
如下图所示:
设定 xyz-轴为参考系的参考轴。称 xy-平面与 xy-平面的相交为交点
线,用英文字母(n)代表。
zxz 顺规的欧拉角可以静态地这样定义:
α 是 x-轴与交点线的夹角,β 是 z-轴与 z-轴的夹角,γ 是交点线与 x-轴的夹角。
3.1.2 四元数
四元数是复数的不可交换延伸。如把四元数的集合考虑成多维实数空间的话,四元数就代表着一个四维空间,相对于复数为二维空间。
四元数的计算,四元数可以理解为一个实数和一个向量的组合,也可以理解为四维的向量。
其中的q为一个四元数,其模的长度为:
对四元数进行单位化,与线性代数中的单位话相似,可得到:
再由创造出来一个变量q关于旋转角得到的一个变量,即可表示为:
由于“四元数表示”转“欧拉角表示”。(这个地方跳过了复杂的换算步骤,我也不太理解)
3.1.3 pid控制
当今的闭环自动控制技术都是基于反馈的概念以减少不确定性。反馈理论的要素包括三个部分:测量、比较和执行。测量关键的是被控变量的实际值,与期望值相比较,用这个偏差来纠正系统的响应,执行调节控制。在工程实际中,应用最为广泛的调节器控制规律为比例、积分、微分控制,简称 pid 控制,又称 pid 调节。
pid 控制器(比例-积分-微分控制器)是一个在工业控制应用中常见的反馈回路部件,由比例单元 p、积分单元 i 和微分单元 d 组成。pid 控制的基础是比例控制;积分控制可消除稳态误差,但可能增加超调;微分控制可加快大惯性系统响应速度以及减弱超调趋势。如下图所示:
3.2 飞控软件框架设计
主要程序设计框图参考了minfly的设计框图:
主要任务关系如下:
说明:此处不包含app设计。
radiolinktask:无线通信任务。该任务主要负责接收从 nrf51822 发送(串口方式)过来的数据,然后对数据进行打包和校验,打包成 atkp 格式并校验无误后发送到atkprxanltask 的接收队列里,同时回传一帧数据给 nrf51822。
usblinkrxtask:usb 通信接收任务。该任务主要负责接收上位机发下来(usb 虚拟串口方式)的数据,然后对数据进行打包和校验,打包成 atkp 格式并校验无误后发送到atkprxanltask 的接收队列里。
atkprxanltask:atkp 数据包接收处理任务。该任务主要是处理遥控器和上位机发下来的数据包,解析到的控制指令则发送到 stabilizertask 中去。
stabilizertask:四轴平衡控制任务。该任务运行的内容比较多,也是比较关键的内容。包括传感器数据读取,数据融合,获取控制数据,空翻检测,异常检测,pid 控制,pwm输出控制等。
wifilinktask:手机控制任务。该任务主要是接收 wifi 摄像头模块的串口数据,然后按照 wifi 摄像头模块通讯协议解析成对应的控制指令,并将控制指令发送到 stabilizertask。
atkptxtask:atkp 数据包发送任务。该任务主要是获取 stabilizertask 中的传感器数据、姿态数据、电机 pwm 输出数据等数据以定周期发送给 radiolinktask 和 usblinktxtask,由这两个任务分别发送飞遥控器和上位机。
usblinkrxtask:usb 通信发送任务。该任务主要负责发送atkptxtask 发送过来的数据包,这些数据包主要是传感器数据、姿态数据等。
3.3 飞控软件开发
3.3.1 姿态解算与pid算法
算法流程图如下:
关于姿态解算,采用互补滤波算法进行姿态解算,更新周期 500hz。mcu 通过iic(模拟 iic)读取加速计和陀螺仪数据寄存器,然后对加速计数据 iir 低通滤波,对陀螺仪数据加偏置调整,然后对加计数据和陀螺数据进行融合,输出姿态数据(roll/pitch/yaw)。
角度环 pid 控制器,更新周期 500hz,z 轴高度 pid 控制器,更新周期 250hz。得到实际油门值和姿态控制量数据,我们就可以把油门值和姿态控制量数据整合,整合周期 1000hz,然后通过控制 pwm 控制电机,从而控制四轴。
目前常见的飞控系统中只使用一个姿态传感器芯片,这个芯片集成了加速度计、陀螺仪以及磁传感器。mpu6050算法主要代码如下:
#include "mpu6050.h"
#include "iic.h"
#include "systick.h"
#include "acc_cal.h"
s16_xyz accraw = {0}; //加速度计原始数据
s16_xyz gyroraw = {0}; //陀螺仪原始数据
si_f_xyz accbutterworthdata = {0}; //加速度计巴特沃斯低通滤波后的数据
si_f_xyz gyrobutterworthdata = {0}; //陀螺仪巴特沃斯低通滤波后的数据
si_f_xyz acc_att_lpf = {0};
si_f_xyz acc_fix_lpf = {0};
si_f_xyz acc_1_lpf = {0};
si_f_xyz acc_butter_lpf = {0};
si_f_xyz gyro_lpf = {0};
si_f_xyz gyro_offset = {0,0,0}; //陀螺仪零偏数据
_mpu6050_data mpu = {0};
//mpu初始化
void mpu6050_init(void)
{
iic_write_one_byte(0xd0,pwr_mgmt_1, 0x80);
delay_ms(100);
iic_write_one_byte(0xd0,pwr_mgmt_1, 0x00); //唤醒mpu
/* when dlpf is disabled( dlpf_cfg=0 or 7),陀螺仪输出频率 = 8khz;
when dlpfis enabled,陀螺仪输出频率 = 1khz
fs(采样频率) = 陀螺仪输出频率 / (1 + smplrt_div)*/
iic_write_one_byte(0xd0,smplrt_div, 0x00); //sample rate. fsample= 1khz/(<this value>+1) = 1000hz
iic_write_one_byte(0xd0,mpu_config, 0x03); //内部低通 acc:44hz gyro:42hz
iic_write_one_byte(0xd0,gyro_config, 0x18); // gyro scale :+-2000deg/s
iic_write_one_byte(0xd0,accel_config, 0x10); // accel scale :+-8g (65536/16=4096 lsb/g)
}
//两字节数据合成
static int getdata(unsigned char reg_address)
{
unsigned char h,l;
h = iic_read_one_byte(0xd0,reg_address);
l = iic_read_one_byte(0xd0,reg_address+1);
return ((h<<8)+l);
}
//get id
uint8_t get_mpu_id(void)
{
u8 mpu_id;
mpu_id = iic_read_one_byte(0xd0,who_am_i);
return mpu_id
//读取陀螺仪三轴数据量
void getgyroraw(void)
{
gyroraw.x = getdata(gyro_xout_h) - gyro_offset.x; //原始数据
gyroraw.y = getdata(gyro_yout_h) - gyro_offset.y;
gyroraw.z = getdata(gyro_zout_h) - gyro_offset.z;
gyrobutterworthdata.x = (float)butterworth_lpf(((float)gyroraw.x),&gyrobutterdata[0],&gyro_30hz_parameter); //巴特沃斯低通滤波后的数据
gyrobutterworthdata.y = (float)butterworth_lpf(((float)gyroraw.y),&gyrobutterdata[1],&gyro_30hz_parameter);
gyrobutterworthdata.z = (float)butterworth_lpf(((float)gyroraw.z),&gyrobutterdata[2],&gyro_30hz_parameter);
}
//求取iir滤波因子
void get_iir_factor(float *out_factor,float time, float cut_off)
{
*out_factor = time /( time + 1/(2.0f * pi * cut_off) );
}
//加速度iir低通滤波
void acc_iir_lpf(si_f_xyz *acc_in,si_f_xyz *acc_out,float lpf_factor)
{
acc_out->x = acc_out->x + lpf_factor*(acc_in->x - acc_out->x);
acc_out->y = acc_out->y + lpf_factor*(acc_in->y - acc_out->y);
acc_out->z = acc_out->z + lpf_factor*(acc_in->z - acc_out->z);
}
//加速度计滤波参数
_butterworth_parameter acc_5hz_parameter =
{
1, -1.778631777825, 0.8008026466657,
0.005542717210281, 0.01108543442056, 0.005542717210281
};
_butterworth_data acc_butter_data[3];
//加速度计巴特沃斯低通
void acc_butterworth_lpf(si_f_xyz *accin,si_f_xyz *accout)
{
accout->x = butterworth_lpf(accin->x,&acc_butter_data[0],&acc_5hz_parameter);
accout->y = butterworth_lpf(accin->y,&acc_butter_data[1],&acc_5hz_parameter);
accout->z = butterworth_lpf(accin->z,&acc_butter_data[2],&acc_5hz_parameter);
}
//原始加速度量转为 g
void accdatatranstog(si_f_xyz *accin,si_f_xyz *accout)
{
accout->x = (float)(accin->x * acc_raw_to_g);
accout->y = (float)(accin->y * acc_raw_to_g);
accout->z = (float)(accin->z * acc_raw_to_g);
}
//滤波后的数据转成(弧度/秒)单位
void radtransform(si_f_xyz *gyroin,si_f_xyz *gyroradout)
{
gyroradout->x = (float)(gyroin->x * gyro_raw_to_radian_s);
gyroradout->y = (float)(gyroin->y * gyro_raw_to_radian_s);
gyroradout->z = (float)(gyroin->z * gyro_raw_to_radian_s);
}
//滤波后的数据转成(度/秒)单位
void degtransform(si_f_xyz *gyroin,si_f_xyz *gyrodegout)
{
gyrodegout->x = (float)(gyroin->x * gyro_raw_to_deg_s);
gyrodegout->y = (float)(gyroin->y * gyro_raw_to_deg_s);
gyrodegout->z = (float)(gyroin->z * gyro_raw_to_deg_s);
}
3.3.2 无线通信软件开发
这里根据上面介绍的初始化程序然后根据需要发送的数据将数据传送到发送和接受缓冲区进行发送与接受。
#include "nrf24l01.h"
#include "spi.h"
#include "systick.h"
#include "led.h"
#include "imath.h"
#include "pair_freq.h"
const u8 tx_address[tx_adr_width]={0x1f,0x2e,0x3d,0x4c,0x5b};
const u8 rx_address[rx_adr_width]={0x1f,0x2e,0x3d,0x4c,0x5b};
void nrf24l01init(void)
{
gpio_inittypedef gpio_initstructure;
rcc_apb2periphclockcmd(rcc_apb2periph_gpiob|rcc_apb2periph_gpioa, enable);
//ce
gpio_initstructure.gpio_pin = gpio_pin_1;
gpio_initstructure.gpio_mode = gpio_mode_out_pp;
gpio_initstructure.gpio_speed = gpio_speed_50mhz;
gpio_init(gpiob, &gpio_initstructure);
//irq
gpio_initstructure.gpio_pin = gpio_pin_4;
gpio_initstructure.gpio_mode = gpio_mode_ipd;
gpio_initstructure.gpio_speed = gpio_speed_50mhz;
gpio_init(gpioa, &gpio_initstructure);
nrf_ce_l;
spi_csn_h;
}
//无线是否在位检测
u8 nrf24l01_check(void)
{
u8 buf[5]={0x18,0x18,0x18,0x18,0x18};
u8 i;
spi_write_buf(nrf_write_reg+tx_addr,buf,5);
spi_read_buf(tx_addr,buf,5);
for(i=0;i<5;i++){
if(buf[i]!=0x18)
break;
}
if(i!=5)
return 1;
return 0;
}
//向寄存器写入值
u8 spi_write_reg(u8 reg,u8 value)
{
u8 status;
spi_csn_l;
status = spi_rw_byte(reg);
spi_rw_byte(value);
spi_csn_h;
return(status);
}
//读取寄存器值
u8 spi_read_reg(u8 reg)
{
u8 reg_val;
spi_csn_l;
spi_rw_byte(reg);
reg_val = spi_rw_byte(0xff);
spi_csn_h;
return(reg_val);
}
//读出寄存器中连续len个字节长度的值
u8 spi_read_buf(u8 reg,u8 *pbuf,u8 len)
{
u8 status,u8_ctr;
spi_csn_l;
status = spi_rw_byte(reg);
for(u8_ctr=0;u8_ctr<len;u8_ctr++)
pbuf[u8_ctr]=spi_rw_byte(0xff);
spi_csn_h;
return status;
}
//向寄存器写入连续len个字节的值
u8 spi_write_buf(u8 reg, u8 *pbuf, u8 len)
{
u8 status,u8_ctr;
spi_csn_l;
status = spi_rw_byte(reg);
for(u8_ctr=0; u8_ctr<len; u8_ctr++)
spi_rw_byte(*pbuf++);
spi_csn_h;
return status;
}
//接收模式
void nrf24l01receivemode(void)
{
nrf_ce_l;
spi_write_reg(setup_aw, 0x03); // 设置地址宽度为 5bytes
spi_write_buf(nrf_write_reg+rx_addr_p0,(u8*)pair.addr,rx_adr_width);//设置接收地址(rx)
spi_write_reg( nrf_write_reg+feature, 0x06 );//使能动态负载长度及ack应答
spi_write_reg(nrf_write_reg+dynpd, 0x01); //使能接收管道0动态负载长度
spi_write_reg(nrf_write_reg+en_aa,0x01); //使能通道0的自动应答
spi_write_reg(nrf_write_reg+en_rxaddr,0x01); //使能通道0的接收地址
spi_write_reg(nrf_write_reg+rf_ch,pair.freq_channel); //设置频点(rf通道)
spi_write_reg(nrf_write_reg+rx_pw_p0,rx_pload_width); //设置接收数据通道0有效数据宽度为11
spi_write_reg(nrf_write_reg+rf_setup,0x07); //设置射频数据率为1mhz,发射功率为7dbm
spi_write_reg(nrf_write_reg+config, 0x0f); //配置基本工作模式的参数;开启crc,配置为接收模式,开启所有中断
nrf_ce_h;
}
//接收数据包
u8 nrf24l01_rxpacket(u8 *rxbuf)
{
u8 sta;
sta = spi_read_reg(nrf_read_reg+status); //状态标志位
spi_write_reg(nrf_write_reg+status,sta);
if(sta&rx_ok) //接收成功
{
spi_read_buf(rd_rx_pload,rxbuf,rx_pload_width);
spi_write_reg(flush_rx,0xff);
return 0;
}
return 1;
}
//该函数初始化nrf24l01到tx模式
//设置tx地址,写tx数据宽度,设置rx自动应答的地址,填充tx发送数据,选择rf频道,波特率和lna hcurr
//pwr_up,crc使能
//当ce变高后,即进入rx模式,并可以接收数据了
//ce为高大于10us,则启动发送.
void nrf24l01_tx_mode(void)
{
nrf_ce_l;
spi_write_reg(setup_aw, 0x03); // 设置地址宽度为 5bytes
spi_write_buf(nrf_write_reg+tx_addr,(uint8_t*)pair.addr,tx_adr_width); //写tx节点地址
spi_write_buf(nrf_write_reg+rx_addr_p0,(uint8_t*)pair.addr,rx_adr_width); //设置tx节点地址,主要为了接收ack
//nrf24l01_write_reg(nrf_write_reg+feature, 0x02 );//使能动态负载长度及带负载的ack应答
//nrf24l01_write_reg(nrf_write_reg+dynpd, 0x01); //使能接收管道0动态负载长度
spi_write_reg(nrf_write_reg+en_aa,0x01); //使能通道0的自动应答
spi_write_reg(nrf_write_reg+en_rxaddr,0x01); //使能通道0的接收地址
spi_write_reg(nrf_write_reg+rf_ch,pair.freq_channel); //设置rf通道
spi_write_reg(nrf_write_reg+setup_retr,0x1a); //设置自动重发间隔时间:500us;最大自动重发次数:10次
spi_write_reg(nrf_write_reg+rf_setup,0x07); //设置射频数据率为1mhz,发射功率为7dbm
spi_write_reg(nrf_write_reg+config,0x0e); //配置基本工作模式的参数;开启crc,配置为发射模式,开启所有中断
nrf_ce_h; //ce为高,10us后启动发送
}
//启动nrf24l01发送一次数据
//sendbuff:待发送数据首地址
//返回值:发送完成状况
uint8_t nrf24l01_txpacket(uint8_t *sendbuff)
{
uint8_t state;
nrf_ce_l;
//nrf24l01_write_buf(spi_write_reg+rx_addr_p0,(uint8_t*)pair.addr,rx_adr_width);
spi_write_buf(wr_tx_pload,sendbuff,tx_pload_width);
nrf_ce_h; //启动发送
while(nrf_irq!=0); //等待发送完成
state=spi_read_reg(nrf_write_reg+status); //读取状态寄存器的值
spi_write_reg(nrf_write_reg+status,state); //清除tx_ds或max_rt中断标志
if(state&max_tx){ //达到最大重发次数
spi_write_reg(flush_tx,0xff); //清除tx fifo寄存器
return max_tx;
}
if(state&tx_ok){ //发送完成
return tx_ok;
}
return 0xff; //其他原因发送失败
}
3.3.3 角度环 pid 和角速度环 pid
pid 更新函数,pid 采用的标准 pid,其数学公式如下:
将数学公式转换为 c 代码,pid 更新函数是这样的:
float pidupdate(pidobject* pid, const float error)
{
float output;
pid->error = error;
pid->integ += pid->error * pid->dt;
if (pid->integ > pid->ilimit)
{
pid->integ = pid->ilimit;
}
else if (pid->integ < pid->ilimitlow)
{
alientek
minifly
28 / 48
atk-minifly 开发指南
pid->integ = pid->ilimitlow;
}
pid->deriv = (pid->error - pid->preverror) / pid->dt;
pid->outp = pid->kp * pid->error;
pid->outi = pid->ki * pid->integ;
pid->outd = pid->kd * pid->deriv;
output = pid->outp + pid->outi + pid->outd;
pid->preverror = pid->error;
return output;
}
pidobject 为 pid 对象结构体数据类型,第一个参数为将被更新的 pid 结构体对象,第二个参数则是偏差(期望值-测量值),积分项为偏差对时间的积分,微分项则是偏差对时间的微分,然后函数里面有三个参数 pid->kp,pid->ki,pid->kd 分别指的是该 pid 对象的比例项,积分项和微分项系数,每个 pid 对象都有属于自己的 pid 系数,pid 初始化 pid 对象的时候会设定一组默认的系数,同时这组系数是可以调整的,我们常说的 pid 参数整定,其实就是调整这组系数,让它满足你的系统。
其函数原型如下:
void attitudeanglepid(attitude_t *actualangle, attitude_t *desiredangle, attitude_t *outdesiredrate)
{/* 角度环 pid */
outdesiredrate->roll = pidupdate(&pidangleroll, desiredangle->roll - actualangle->roll);
outdesiredrate->pitch = pidupdate(&pidanglepitch, desiredangle->pitch - actualangle->pitch);
float yawerror = desiredangle->yaw - actualangle->yaw ;
if (yawerror > 180.0f)
yawerror -= 360.0f;
else if (yawerror < -180.0)
yawerror += 360.0f;
outdesiredrate->yaw = pidupdate(&pidangleyaw, yawerror);
}
attitude_t 是一个姿态数据结构类型,函数参数 actualangle 是一个结构体指针,指向实际角度结构体变量(数据融合输出值)state->attitude, desiredangle 指向期望角度结构体变量(设置的角度)attitudedesired,outdesiredrate 则是角度环的输出,指向期望角速度结构体变量 ratedesired。
然后是角速度环 pid,其函数原型如下:
void attituderatepid(axis3f *actualrate, attitude_t *desiredrate, control_t *output)
{/* 角速度环 pid */
output->roll = pidoutlimit(pidupdate(&pidrateroll, desiredrate->roll - actualrate->x));
output->pitch = pidoutlimit(pidupdate(&pidratepitch, desiredrate->pitch - actualrate->y));
alientek
minifly
29 / 48
atk-minifly 开发指南
output->yaw = pidoutlimit(pidupdate(&pidrateyaw, desiredrate->yaw - actualrate->z));
}
3.3.4 姿态控制量和油门值整合
设置x 模式飞行,电机转向和姿态解算正方向(箭头指示正方向):
control->thrust 为油门控制量,这个值增大四轴升高,减小则下降。control->roll,control->pitch,control->yaw 为 pid 输出的姿态控制量。油门控制量和姿态控制量整合后控制电机,整合代码在 power_control.c 文件的函数 powercontrol ()中实现,代码如下:
void powercontrol(control_t *control) /*功率输出控制*/
{
s16 r = control->roll / 2.0f;
s16 p = control->pitch / 2.0f;
motorpwm.m1 = limitthrust(control->thrust - r - p + control->yaw);
motorpwm.m2 = limitthrust(control->thrust - r + p - control->yaw);
alientek
minifly
31 / 48
atk-minifly 开发指南
motorpwm.m3 = limitthrust(control->thrust + r + p + control->yaw);
motorpwm.m4 = limitthrust(control->thrust + r - p - control->yaw);
if (motorsetenable)
{
motorpwm=motorpwmset;
}
motorssetratio(motor_m1, motorpwm.m1); /*控制电机输出百分比*/
motorssetratio(motor_m2, motorpwm.m2);
motorssetratio(motor_m3, motorpwm.m3);
motorssetratio(motor_m4, motorpwm.m4);
}
roll 方向受外力向左旋转(向右为正),为了恢复平衡,则 m3 和 m4 同侧出力,m1和 m2 反向出力(m1 和 m2 的 roll 为-,m3 和 m4 的 roll 为+);
pitch 方向受外力向后旋转(向前为正),为了恢复平衡,则 m2 和 m3 同侧出力,m1和 m4 反向出力(m1 和 m4 的 pitch 为-,m2 和 m3 的 pitch 为+);
yaw 方向受外力顺时针旋转(逆时针为正),为了恢复平衡,则 m1 和 m3 同侧出力,m2 和 m4 反向出力,(m2 和 m4 的 yaw 为-,m1 和 m3 的 yaw 为+);
bool 型变量 motorsetenable为 true,使能手动设置电机占空比,这样可以方便单独调试某几个电机,默认不使能。
motorssetratio()当然就是设定对应电机定时器通道占空比的函数了,设定的占空比作用到 mos 管,然后控制电机,从而控制四轴。
3.3.5 4d 空翻算法
4d 空翻实现原理是只使用内环 pid–角速度环 pid 控制器,姿态角度期望值直接作为角速度环的期望值,测量值使用 3 轴陀螺仪数据,这样我们控制的不是四轴的角度,而是四轴的转动角速度。知道了如何控制四轴转动,当然就能控制翻滚了。
4d 空翻源码比较多,我就不贴出来了,空翻具体实现过程自行去 flip.c 文件查看空翻实现函数 flyerflipcheck(),空翻过程有好几个状态,flipstate 指示空翻的当前状态,其定义如下:
static enum
{
flip_idle = 0,
flip_set,
flip_speed_up,
flip_slow_down,
alientek
minifly
32 / 48
atk-minifly 开发指南
flip_period,
flip_finished,
rever_speed_up,
flip_error,
}flipstate = flip_idle;
flip_idle 为空翻空闲状态,在此状态下,四轴实时检测是否要执行空翻命令。如果检测到空翻指令,则状态切换到 flip_set,在此状态下,我们设置一些四轴翻滚用到的参数,参数设置完成后切换到加速上升状态 flip_speed_up,因为空翻特技会有掉高问题,所以我们在真正 4d 翻滚之前先加速一段时间,当 z 轴速度达到一定值后,进入减速状态,为什么翻滚之前要这个减速状态呢,答案是为了更好的空翻。减速到设定值后才进入真正的翻滚状态 flip_period,前面状态都是为空翻做准备的。
3.4 限制于篇幅,其他模块的算法不做说明(太多了……)
四、飞控的制作
制作之前需要将原理图导入pcb并设计出板子形状。
4.1 pcb板图
为方便后续封装的配置,将原理图导入到ad软甲中进行设计,板图如下:分别为飞机和遥控器板图:
将设计图纸发给厂家,最终得到上面板子的板子:
4.2 贴片器件的焊接
因为以前没有焊接贴片式器件的经验,所以去工作室练习了芯片灯复杂器件的焊接。完成这些之前,你需要有自己的加热台,有条件的可以买个热风枪。
将器件按照原理图中的位置进行摆放,我采购的器件的封装为0603封装(仅限电阻电感),如下图:
4.3 实物图
最终可以得到实物如下:
五、飞控程序烧录及调试
焊接好器件就可以进行程序的烧录以及+
器件的调试工作了。
在这里我想说一下我遇到的一些问题,以及解决的方法。在调试的过程,我重复了上面3次的焊接步骤,具体原因如下:
-
芯片焊接出现问题,有了连锡的现象
-
mpu6050封装加热时出现,由于加热时间过长,导致焊盘脱落
-
遥控器焊接完毕,oled显示屏不能显示,最终发现是程序出现了问题,最终导致芯片烧坏,被迫换板
如下图,分别是上述三种原因导致的焊接失败的板子:
经过三次的失败,我得到了自己的焊接芯片的技术方法,最终焊接成功。如图:
程序的烧录采用的是keil软件进行烧录:
此处显示编译无误,即可调试进行烧录程序,如果芯片焊接无误这里显示的是这样的:
或者采用stm32 unity进行验证,单击这里:
如果出现如图所示,证明芯片焊接没有问题:
之后就在keil中进行烧录程序:单机load进行烧录,等待完成。遥控器也是同样的操作,完成步骤后如图:
然后成品就完成了,后面我根据需求,自己加装了天线。
六、 总结
装箱后的飞机就是这样,目前还没有试飞,有空时候去外边飞一飞看,整机花销100以内,还是很亲民的,哈哈哈哈哈
创作不易,这个飞机花费了太多的时间,希望粉丝朋友们喜欢,谢谢大家!!!
发表评论