尝试制作这个四旋翼飞控的过程,感触颇多,整理了思绪之后,把重要的点一一记下来;

这个飞控是基于STM32,整合了MPU6050,即陀螺仪和重力加速计,但没有融合电子罗盘;

另外,四旋翼飞行器的运动方式请百度百科,不太复杂,具体不再赘述;

这是飞控程序的控制流程(一个执行周期):

比较重要的地方:

1.i2c通信方式;

  因为我不是学电类专业,最开始对i2c这些是没有一点概念,最后通过Google了解了一些原理,然后发现STM32的开发库是带有i2c通信的相关函数的,但是我最后还是没有用这些函数。

我通过GPIO模拟i2c,这样也能获得mpu6050的数据,虽然代码多了一些,但是比较好的理解i2c的原理。

  STM32库实现的模拟i2c代码(注释好像因为编码问题跪了):

/*******************************************************************************

// file                :    i2c_conf.h

// MCU                : STM32F103VET6

// IDE                : Keil uVision4

// date                £º2014.2.28

*******************************************************************************/

#include "stm32f10x.h"

#define   uchar unsigned char

#define   uint unsigned int

#define      FALSE 0

#define   TRUE  1

void I2C_GPIO_Config(void);

void I2C_delay(void);

void delay5ms(void);

int I2C_Start(void);

void I2C_Stop(void);

void I2C_Ack(void);

void I2C_NoAck(void);

int I2C_WaitAck(void);

void I2C_SendByte(u8 SendByte);

unsigned char I2C_RadeByte(void);

int Single_Write(uchar SlaveAddress,uchar REG_Address,uchar REG_data);

unsigned char Single_Read(unsigned char SlaveAddress,unsigned char REG_Address);

/*******************************************************************************

// file                :  i2c_conf.c

// MCU                : STM32F103VET6

// IDE                : Keil uVision4

// date                £º2014.2.28

*******************************************************************************/

#include "i2c_conf.h"

#define SCL_H         GPIOB->BSRR = GPIO_Pin_6

#define SCL_L         GPIOB->BRR  = GPIO_Pin_6

#define SDA_H         GPIOB->BSRR = GPIO_Pin_7

#define SDA_L         GPIOB->BRR  = GPIO_Pin_7

#define SCL_read      GPIOB->IDR  & GPIO_Pin_6        //IDR:¶Ë¿ÚÊäÈë¼Ä´æÆ÷¡£

#define SDA_read      GPIOB->IDR  & GPIO_Pin_7

void I2C_GPIO_Config(void)

{

GPIO_InitTypeDef  GPIO_InitStructure;

GPIO_InitStructure.GPIO_Pin =  GPIO_Pin_6;

GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_OD;   //¿ªÂ©Êä³öģʽ

GPIO_Init(GPIOB, &GPIO_InitStructure);

GPIO_InitStructure.GPIO_Pin =  GPIO_Pin_7;

GPIO_InitStructure.GPIO_Speed = GPIO_Speed_50MHz;

GPIO_InitStructure.GPIO_Mode = GPIO_Mode_Out_OD;

GPIO_Init(GPIOB, &GPIO_InitStructure);

}

void I2C_delay(void)

{

int i=6; //ÕâÀï¿ÉÒÔÓÅ»¯ËÙ¶È    £¬¾­²âÊÔ×îµÍµ½5»¹ÄÜдÈë

while(i)

{

i--;

}

}

void delay5ms(void)

{

int i=5000;

while(i)

{

i--;

}

}

int I2C_Start(void)

{

SDA_H;                                            //II2ЭÒé¹æ¶¨±ØÐëÔÚʱÖÓÏßΪµÍµçƽµÄÇ°ÌáÏ£¬²Å¿ÉÒÔÈà Êý¾ÝÏßÐźŸıä

SCL_H;

I2C_delay();

if(!SDA_read)

return FALSE;                //SDAÏßΪµÍµçƽÔò×ÜÏßæ,Í˳ö

SDA_L;

I2C_delay();

if(SDA_read)

return FALSE;                //SDAÏßΪ¸ßµçƽÔò×ÜÏß³ö´í,Í˳ö

SDA_L;

I2C_delay();

return TRUE;

}

void I2C_Stop(void)

{

SCL_L;

I2C_delay();

SDA_L;

I2C_delay();

SCL_H;

I2C_delay();

SDA_H;

I2C_delay();

}

void I2C_Ack(void)

{

SCL_L;

I2C_delay();

SDA_L;

I2C_delay();

SCL_H;

I2C_delay();

SCL_L;

I2C_delay();

}

void I2C_NoAck(void)

{

SCL_L;

I2C_delay();

SDA_H;

I2C_delay();

SCL_H;

I2C_delay();

SCL_L;

I2C_delay();

}

int I2C_WaitAck(void)      //·µ»ØΪ:=1ÓÐACK,  =0ÎÞACK

{

SCL_L;

I2C_delay();

SDA_H;

I2C_delay();

SCL_H;

I2C_delay();

if(SDA_read)

{

SCL_L;

I2C_delay();

return FALSE;

}

SCL_L;

I2C_delay();

return TRUE;

}

void I2C_SendByte(u8 SendByte) //Êý¾Ý´Ó¸ßλµ½µÍλ//

{

u8 i=8;

while(i--)

{

SCL_L;

I2C_delay();

if(SendByte&0x80)        // 0x80 = 1000 0000;

SDA_H;

else

SDA_L;

SendByte<<=1;   // SendByte×óÒÆһλ¡£

I2C_delay();

SCL_H;

I2C_delay();

}

SCL_L;

}

unsigned char I2C_RadeByte(void)  //Êý¾Ý´Ó¸ßλµ½µÍλ//

{

u8 i=8;

u8 ReceiveByte=0;

SDA_H;

while(i--)

{

ReceiveByte<<=1;      //×óÒÆһ룬

SCL_L;

I2C_delay();

SCL_H;

I2C_delay();

if(SDA_read)

{

ReceiveByte"=0x01; //дÈë

}

}

SCL_L;

return ReceiveByte;

}

int Single_Write(uchar SlaveAddress,uchar REG_Address,uchar REG_data)

{

if(!I2C_Start())

return FALSE;

I2C_SendByte(SlaveAddress);   //·¢ËÍÉ豸µØÖ·+дÐźŠ   //I2C_SendByte(((REG_Address & 0x0700) >>7) | SlaveAddress & 0xFFFE);    //ÉèÖøßÆðʼµØÖ·+Æ÷¼þµØÖ·

if(!I2C_WaitAck())

{

I2C_Stop();

return FALSE;

}

I2C_SendByte(REG_Address );   //ÉèÖõÍÆðʼµØÖ·

I2C_WaitAck();

I2C_SendByte(REG_data);

I2C_WaitAck();

I2C_Stop();

delay5ms();

return TRUE;

}

unsigned char Single_Read(unsigned char SlaveAddress,unsigned char REG_Address)

{

unsigned char REG_data;

if(!I2C_Start())

return FALSE;

I2C_SendByte(SlaveAddress); //I2C_SendByte(((REG_Address & 0x0700) >>7) | REG_Address & 0xFFFE);//ÉèÖøßÆðʼµØÖ·+Æ÷¼þµØÖ·

if(!I2C_WaitAck())

{

I2C_Stop();

return FALSE;

}

I2C_SendByte((u8) REG_Address);   //ÉèÖõÍÆðʼµØÖ·

I2C_WaitAck();

I2C_Start();

I2C_SendByte(SlaveAddress+1);

I2C_WaitAck();

REG_data= I2C_RadeByte();

I2C_NoAck();

I2C_Stop();

return REG_data;

}

2.mpu6050;

  然后用写好的模拟i2c函数读取mpu6050,根据mpu6050手册的各寄存器地址,读取到了重力加速计和陀螺仪的各分量;

传感器采样率设置为200Hz,这个值是因为我电调频率为200Hz,也就是说,我的程序循环一次0.005s,一般来说,采样率高点没问题,别比执行一次闭环控制的周期长就行了;

陀螺仪量程±2000°/s,加速计量程±2g, 量程越大,取值越不精确;

这里注意,由于我们没有采用磁力计,而陀螺仪存在零偏,所以最终在yaw方向上没有绝对的参考系,不能建立绝对的地理坐标系,这样最好的结果也仅仅是在yaw上存在缓慢漂移。

3.互补滤波;

  融合时,陀螺仪的积分运算很大程度上决定了飞行器的瞬时运动情况,而重力加速计通过长时间的累积不断矫正陀螺仪产生的误差,最终得到准确的机体姿态。

  这里我们采用Madgwick提供的UpdateIMU算法来得到姿态角所对应的四元数,之后只需要经过简单运算便可转换为实时欧拉角。感谢Madgwick大大为开源做出的贡献。

1 #define sampleFreq    512.0f        // sample frequency in Hz

2 #define betaDef        0.1f        // 2 * proportional gain

3

4

5 volatile float beta = betaDef;

6 volatile float q0 = 1.0f, q1 = 0.0f, q2 = 0.0f, q3 = 0.0f;

7

8 float invSqrt(float x);

9

10 void MadgwickAHRSupdateIMU(float gx, float gy, float gz, float ax, float ay, float az) {

11     float recipNorm;

12     float s0, s1, s2, s3;

13     float qDot1, qDot2, qDot3, qDot4;

14     float _2q0, _2q1, _2q2, _2q3, _4q0, _4q1, _4q2 ,_8q1, _8q2, q0q0, q1q1, q2q2, q3q3;

15

16     // Rate of change of quaternion from gyroscope

17     qDot1 = 0.5f * (-q1 * gx - q2 * gy - q3 * gz);

18     qDot2 = 0.5f * (q0 * gx + q2 * gz - q3 * gy);

19     qDot3 = 0.5f * (q0 * gy - q1 * gz + q3 * gx);

20     qDot4 = 0.5f * (q0 * gz + q1 * gy - q2 * gx);

21

22     // Compute feedback only if accelerometer measurement valid (avoids NaN in accelerometer normalisation)

23     if(!((ax == 0.0f) && (ay == 0.0f) && (az == 0.0f))) {

24

25         // Normalise accelerometer measurement

26         recipNorm = invSqrt(ax * ax + ay * ay + az * az);

27         ax *= recipNorm;

28         ay *= recipNorm;

29         az *= recipNorm;

30

31         // Auxiliary variables to avoid repeated arithmetic

32         _2q0 = 2.0f * q0;

33         _2q1 = 2.0f * q1;

34         _2q2 = 2.0f * q2;

35         _2q3 = 2.0f * q3;

36         _4q0 = 4.0f * q0;

37         _4q1 = 4.0f * q1;

38         _4q2 = 4.0f * q2;

39         _8q1 = 8.0f * q1;

40         _8q2 = 8.0f * q2;

41         q0q0 = q0 * q0;

42         q1q1 = q1 * q1;

43         q2q2 = q2 * q2;

44         q3q3 = q3 * q3;

45

46         // Gradient decent algorithm corrective step

47         s0 = _4q0 * q2q2 + _2q2 * ax + _4q0 * q1q1 - _2q1 * ay;

48         s1 = _4q1 * q3q3 - _2q3 * ax + 4.0f * q0q0 * q1 - _2q0 * ay - _4q1 + _8q1 * q1q1 + _8q1 * q2q2 + _4q1 * az;

49         s2 = 4.0f * q0q0 * q2 + _2q0 * ax + _4q2 * q3q3 - _2q3 * ay - _4q2 + _8q2 * q1q1 + _8q2 * q2q2 + _4q2 * az;

50         s3 = 4.0f * q1q1 * q3 - _2q1 * ax + 4.0f * q2q2 * q3 - _2q2 * ay;

51         recipNorm = invSqrt(s0 * s0 + s1 * s1 + s2 * s2 + s3 * s3); // normalise step magnitude

52         s0 *= recipNorm;

53         s1 *= recipNorm;

54         s2 *= recipNorm;

55         s3 *= recipNorm;

56

57         // Apply feedback step

58         qDot1 -= beta * s0;

59         qDot2 -= beta * s1;

60         qDot3 -= beta * s2;

61         qDot4 -= beta * s3;

62     }

63

64     // Integrate rate of change of quaternion to yield quaternion

65     q0 += qDot1 * (1.0f / sampleFreq);

66     q1 += qDot2 * (1.0f / sampleFreq);

67     q2 += qDot3 * (1.0f / sampleFreq);

68     q3 += qDot4 * (1.0f / sampleFreq);

69

70     // Normalise quaternion

71     recipNorm = invSqrt(q0 * q0 + q1 * q1 + q2 * q2 + q3 * q3);

72     q0 *= recipNorm;

73     q1 *= recipNorm;

74     q2 *= recipNorm;

75     q3 *= recipNorm;

76 }

77

78

79 float invSqrt(float x) {

80     float halfx = 0.5f * x;

81     float y = x;

82     long i = *(long*)&y;

83     i = 0x5f3759df - (i>>1);

84     y = *(float*)&i;

85     y = y * (1.5f - (halfx * y * y));

86     return y;

87 }

  四元数转欧拉角的算法可参考 http://www.cnblogs.com/wqj1212/archive/2010/11/21/1883033.html

4.获取期望姿态;

  也就是遥控部分了,让用户介入控制。

  本着拿来主义的原则,用上”圆点博士开源项目”提供的安卓的开源蓝牙控制端。

  圆点博士给出了数据包格式,同过HC-06蓝牙模块接连到STM32串口1,再无线连接到控制端,这样我们就可以获得控制端不断发送的数据包了,并实时更新期望姿态角,这里只需要注意输出的姿态角和实时姿态角方向一致以及数据包的校验就行了。

5.PID控制算法;

  由于简单的线性控制不可能满足四轴飞行器这个灵敏的系统,引入PID控制器来更好的纠正系统。

  简介:PID实指“比例proportional”、“积分integral”、“微分derivative”,这三项构成PID基本要素。每一项完成不同任务,对系统功能产生不同的影响。

以Pitch为例:

  

error为期望角减去实时角度得到的误差;

iState为积分i参数对应累积过去时间里的误差总和;

if语句限定iState范围,繁殖修正过度;

微分d参数为当前姿态减去上次姿态,估算当前速度(瞬间速度);

总调整量为p,i,d三者之和;

这样,P代表控制系统的响应速度,越大,响应越快。

I,用来累积过去时间内的误差,修正P无法达到的期望姿态值(静差);

D,加强对机体变化的快速响应,对P有抑制作用。

PID各参数的整定需要综合考虑控制系统的各个方面,才能达到最佳效果。

输出PWM信号:

PID计算完成之后,便可以通过STM32自带的定时资源很容易的调制出四路pwm信号,采用的电调pwm格式为50Hz,高电平持续时间0.5ms-2.5ms;

我以1.0ms-2.0ms为每个电机的油门行程,这样,1ms的宽度均匀的对应电调的从最低到最高转速。

至此,一个用stm32和mpu6050搭建的飞控系统就算实现了。

剩下的调试PID参数了。

四旋翼飞行器的飞控实现相关推荐

  1. 四旋翼飞行器12——飞控接收机的三种接收模式

    四旋翼飞行器12--飞控接收机的三种接收模式 文章目录 四旋翼飞行器12--飞控接收机的三种接收模式 一.PWM型 二.PPM型 三.SBUS 无线遥控就是利用高频无线电波实现对模型的控制.如天地飞的 ...

  2. STM32单片机:四旋翼飞行器的飞控实现

    ---------------------------------------------------------------------------------------------------- ...

  3. 四旋翼飞行器7——主要开源飞控介绍

    四旋翼飞行器7--主要开源飞控介绍 文章目录 四旋翼飞行器7--主要开源飞控介绍 一.APM 特性包括: 硬件包括: 二 Pixhawk 特性包括: 各种恰到好处的传感器. 各种可扩展I/O接口和专用 ...

  4. 【飞控理论】四旋翼飞行器控制原理

    该篇博客是对锡月科技无人机飞行控制原理教学内容的整理 1.四旋翼飞行器的两种结构: "×"字模式: Pitch和 Roll与1,3.2,4两组电机呈 45°夹角 . "十 ...

  5. 四旋翼飞行器基本知识(四旋翼飞行器结构和原理+四轴飞行diy全套入门教程)

    转载两篇日志: 第一篇<四旋翼飞行器结构和原理> 第二篇<四旋翼飞行diy全套入门教程> =========================================== ...

  6. 四旋翼飞行器5——各类方案设计及参考

    下面几篇文章及资料是这两天看到的比较好的资源,可以多看看,对于新手来说,是个不错的提升自己的一些资料.可以分析每一种方案的特点,主控啊,传感器啊,不过到最后肯定要自己实践,才能得到最好的控制效果. 文 ...

  7. 四旋翼飞行器平衡传感器数据 处理方法探讨

    四旋翼飞行器平衡传感器数据 处理方法探讨 [摘要]四旋翼飞行器的制作与研究日渐火热,而且其应用价值日渐凸显,当今四旋翼飞行器大都使用飞控板控制,其优点是操作简单,稳定性比较好,但好多的复杂功能无法自定 ...

  8. 四旋翼飞行器2——自己设计四旋翼飞行器的硬件准备和基础知识

    四旋翼飞行器2--自己设计四旋翼飞行器的硬件准备和基础知识 文章目录 四旋翼飞行器2--自己设计四旋翼飞行器的硬件准备和基础知识 DIY四轴需要准备什么零件 [基本原理与名词解释] 1.遥控器篇 什么 ...

  9. 四旋翼飞行器9——APM地面站使用详解(校准结束)

    四旋翼飞行器9--APM地面站使用详解 如果你手上有一块apm飞控,那么APM飞控第一步--下载APM地面站,校准传感器.遥控器. 下面的文章,大部分参考这个: http://bbs.loveuav. ...

最新文章

  1. 关于什么事情能做到和不能做到的思考
  2. LeNet试验(四)使用shuffle label观察网络的记忆能力
  3. .NET Core微服务 权限系统+工作流(一)权限系统
  4. c++读取文件–结束条件的判断
  5. php维护session,维护带有cookie的PHP session_start()
  6. 创建下标为1-10的整形数组
  7. 火狐 移动 html 元素,python中的Firefox+Selenium:如何交互式地获取元素html?
  8. 腾讯公司被法院强制执行25元,刚好是QQ超级会员的价格
  9. Linux内存管理:内存描述之内存页面page
  10. android studio编译找不到程序包,AndroidStudio编译提示“程序包R不存在”可能的解决方案,希望你用不到。。。-Go语言中文社区...
  11. jquery multiSelect 多选下拉框
  12. BZOJ 2337 XOR和路径(概率DP)
  13. mkfs.ext3 快速格式化_U盘格式化
  14. Linux安装及管理程序
  15. CAD转图片如何调整输出格式?
  16. Fe3O4纳米颗粒的表面接枝修饰/氨基乙酸|L-半胱氨酸(L-Cys)修饰的Fe3O4包裹TiO2(Fe3O4@TiO2/L-Cys)复合纳米粒子
  17. es6 arry fill
  18. 你想要的二维码美化效果都在这里了(前景图,圆点二维码,液化效果)
  19. C语言--逗号运算符及逗号表达式
  20. 每天干的啥?(2018.09)

热门文章

  1. [市民呼声] 双色球3.5亿大奖爆出惊天骗局!(转载)
  2. 微软CRM,路在何方?
  3. JS_自己写的JQ插件
  4. 兔子繁殖问题与解决方案
  5. HTTP响应头信息泄露
  6. 画数轴的步骤_如何用几何画板画数轴?
  7. 不用辅助空间的归并排序
  8. 辞职的时候要告诉领导已经找好工作了吗
  9. Date时间格式转换
  10. 翻译:Identifying Encrypted Malware Traffic with Contextual Flow Data利用上下文流数据识别加密恶意软件流量