MPU6050数据轻松分析.docx
《MPU6050数据轻松分析.docx》由会员分享，可在线阅读，更多相关《MPU6050数据轻松分析.docx（14页珍藏版）》请在冰点文库上搜索。
MPU6050数据轻松分析
MPU6050数据轻松分析
这个文章是根据自己学习，查资料的汇总，同时把一些自己的心得加进去。
如果有什么不对的，欢迎请大家指正、交流。
邮箱：
z
最近看到加速度计和陀螺仪比较火，而且也有很多人都在研究。
于是也在网上淘了一个mpu6050模块，想用来做自平衡小车。
可是使用起来就发愁了。
网上关于mpu6050的资料的确不少，但是大家都是互相抄袭，然后贴出一段程序，看完之后还是不知道所以然。
经过翻阅各个方面的资料，以及自己的研究在处理mpu6050数据方面有一些心得，在这里和大家分享一下。
在处理加速度计和陀螺仪用到的方法都是比较简单的，这里的简单并不是不需要任何基础知识，只是这些基本知识都是最基本的，比如简单的三角函数，数学计算，物理知识，c语言以及基本的arduino知识（如果不会arduino会其它单片机也是一样的，本文实践是使用arduino），如果还不具备这些知识那就快去补课吧。
1、加速度和陀螺仪原理
当然，在开始之前至少要弄懂什么是加速度计，什么是陀螺仪吧，否则那后边讲的都是没有意义的。
简单的说，加速度计主要是测量物体运动的加速度，陀螺仪主要测量物体转动的角速度。
这些理论的知识我就不多说了，都可以在网上查到。
这里推荐一篇讲的比较详细的文章《AGuideTousingIMU（AccelerometerandGyroscopeDevices）inEmbeddedApplications》，在网上可以直接搜索到。
2、加速度测量
在开始之前，不知大家是否还记得加速度具有合成定理？
如果不记得可以先大概了解一下，其实简单的举个例子来说就是重力加速度可以理解成是由x,y,z三个方向的加速度共同作用的结果。
反过来说就是重力加速度可以分解成x,y,z三个方向的加速度。
加速度计可以测量某一时刻x,y,z三个方向的加速度值。
而自平衡小车利用加速度计测出重力加速度在x,y,z轴的分量，然后利用各个方向的分量与重力加速度的比值来计算出小车大致的倾角。
其实在自平衡小车上非静止的时候，加速度计测出的结果并不是非常精确。
因为大家在高中物理的时候都学过，物体时刻都会受到地球的万有引力作用产生一个向下的重力加速度，而小车在动态时，受电机的作用肯定有一个前进或者后退方向的作用力，而加速度计测出的结果是，重力加速度与小车运动加速度合成得到一个总的加速度在三个方向上的分量。
不过我们暂时不考虑电机作用产生的运动加速度对测量结果的影响。
因为我们要先把复杂的事情分解成一个个简单的事情来分析，这样才能看到成果，才会有信心继续。
下边我们就开始分析从加速度得到角度的方法。
如下图，把加速度计平放，分别画出xyz轴的方向。
这三个轴就是我们后边分析所要用到的坐标系。
把mpu6050安装在自平衡车上时也是这样的水平安装在小车底盘上的，假设两个车轮安装时车轴和y轴在一条直线上。
那么小车摆动时，参考水平面就是桌面，并且车轴（y轴）与桌面始终是平行的，小车摆动和移动过程中y轴与桌面的夹角是不会发生变化的，一直是0度。
发生变化的是x轴与桌面的夹角以及z轴与桌面的夹角，而且桌面与x轴z轴夹角变化度数是一样的。
所以我们只需要计算出x轴和z轴中任意一个轴的夹角就可以反映出小车的倾斜的情况了。
为了方便分析，由于y轴与桌面夹角始终不变，我们从y轴的方向俯看下去，那么这个问题就会简化成只有x轴和z轴的二维关系。
假设某一时刻小车上加速度计（mpu6050）处于如下状态，下图是我们看到简化后的模型。
在这个图中，y轴已经简化和坐标系的原点o重合在了一起。
我们来看看如何计算出小车的倾斜角，也就是与桌面的夹角a。
上图g是重力加速度，gx、gz分别是g在x轴和z轴的分量。
由于重力加速度是垂直于水平面的，得到：
角a+角b=90度
X轴与y轴是垂直关系，得到：
角c+角b=90度
于是轻松的就可以得出：
角a=角c
根据力的分解，g、gx、gz三者构成一个长方形，根据平行四边形的原理可以得出：
角c=角d
所以计算出角度d就等效于计算出了x轴与桌面的夹角a。
前边已经说过gx是g在x轴的分量，那么根据正弦定理就可以得出：
Sind=gx/g
得到这个公式可是还是得不到想要的角度，因为需要计算反正弦，而反正弦在单片机里不是很好计算。
为了得到角度，于是又查了相关资料，原来在角度较小的情况下，角度的正弦与角度对应的弧度成线性关系。
先看看下边的图：
这个图x轴是角度，取值范围是0~90度，有三个函数曲线，分别是：
Y=sinx正弦曲线
Y=x*3.14/180弧度
Y=0.92*x*3.14/180乘以一个0.92系数的弧度
从图上可以看出，当角度范围是0~29度时：
sinx=x*3.14/180
对于自平衡车来说，小车的摆动范围在-29~29度之内，如果超过这个范围，小车姿态也无法调整，所以对于自平衡小车sinx=x*3.14/180基本上是成立的。
当然有时候也会担心-29~29度的摇摆范围还是无法满足需求。
那可以给上边的公式乘一个系数。
得到如下公式：
Sinx=k*x*3.14/180
从上边的函数对比图可以看出，当系数取0.92时，角度范围可以扩大到-45~45度。
经过这一系列的分析，终于得到角度换算方法：
由
Sind=gx/g
Sind=k*d*3.14/180
得到：
gx/g=k*d*3.14/180
那么角度就可以通过如下公式计算出：
d=180*gx/（k*g*3.14）
而gx可以从加速度计里读出来，所以这下角度就可以轻松得到了。
而之前也说过这个角度不是很精确，但是至少可以反映出角度变化的趋势。
不过可以通过卡尔曼滤波等算法把加速度计读出的角度和陀螺仪读出的角度结合起来，使小车的角度更加准确。
3、陀螺仪
通过陀螺仪来测量角度就很简单了，因为陀螺仪读出的是角速度，大家都知道，角速度乘以时间，就是转过的角度。
把每次计算出的角度做累加就会等到当前所在位置的角度。
先看下图：
假设最初陀螺仪是与桌面平行，单片机每tms读一次陀螺仪的角速度，当读了三次角速度以后z轴转到上图的位置，则在这段时间中转过的角度为x：
角x=角1+角2+角3
假设从陀螺仪读出的角速度为w，那总角度为：
X=（w1*t1+w2*t2+w3*t3）/1000
假设经过n次，那么总的角度如下：
X=（w1*t1+w2*t2+w3*t3+…+wn*tn）/1000
实际上这就是一个积分过程。
其实这种计算出来的角度也存在一定的误差，而且总的角度是经过多次相加得到的，这样误差就会越积累越大，最终导致计算出的角度与实际角度相差很大。
于是也可以使用卡尔曼滤波把加速度计读出的角度结合在一起，使计算出的角度更准确。
4、mpu6050与arduino连接
为了方便大家学习和复制，同时也应大家的要求，详细列出从哪里可以找到，对比了很多觉得这个是最划算和实用的，供大家参考
名称
来源
Mpu6050模块
Arduino模块
其它
Mpu6050可以直接买现成的模块，模块使用i2c接口和arduino通信，连接方式如下：
Mpu6050
Arduino
VCC
VCC
GND
GND
SCL
AD5
SDA
AD4
5、Arduino代码
在开始之前先介绍一下用到的组件和算法。
Mpu6050是通过i2c接口和arduino连接，所以需要i2c库，另外mpu6050也有现成的库。
由于陀螺仪本身就有偏差，多次积累之后偏差会越来越大越来越不准确，所以需要用一些滤波算法来校正，这里用的是卡尔曼滤波算法。
卡尔曼滤波在网上讲的很多，这里就不多说了。
网上有一个测量温度的通俗版，另外维基百科里的介绍比较简单。
有兴趣的可以去研究一下。
另外也可以用其它算法来滤波，这里也只是做演示所以不多介绍。
//ArduinoWirelibraryisrequiredifI2CdevI2CDEV_ARDUINO_WIREimplementation
//isusedinI2Cdev.h
#include"Wire.h"
//I2CdevandMPU6050mustbeinstalledaslibraries,orelsethe.cpp/.hfiles
//forbothclassesmustbeintheincludepathofyourproject
#include"I2Cdev.h"
#include"MPU6050.h"
//classdefaultI2Caddressis0x68
//specificI2Caddressesmaybepassedasaparameterhere
//AD0low=0x68（defaultforInvenSenseevaluationboard）
//AD0high=0x69
MPU6050accelgyro;
int16_tax,ay,az;
int16_tgx,gy,gz;
doubletotal_angle=0;
#defineLED_PIN13
/*把mpu6050放在水平桌面上，分别读取读取2000次，然后求平均值*/
#defineAX_ZERO（-1476）/*加速度计的0偏修正值*/
#defineGX_ZERO（-30.5）/*陀螺仪的0偏修正值*/
voidsetup（）{
//joinI2Cbus（I2Cdevlibrarydoesn'tdothisautomatically）
Wire.begin（）;
//initializeserialcommunication
//（38400chosenbecauseitworksaswellat8MHzasitdoesat16MHz,but
//it'sreallyuptoyoudependingonyourproject）
Serial.begin（38400）;
//initializedevice
Serial.println（"InitializingI2Cdevices..."）;
accelgyro.initialize（）;
//verifyconnection
Serial.println（"Testingdeviceconnections..."）;
Serial.println（accelgyro.testConnection（）?
"MPU6050connectionsuccessful":
"MPU6050connectionfailed"）;
}
/*通过卡尔曼滤波得到的最终角度*/
floatAngle=0.0;
/*由角速度计算的倾斜角度*/
floatAngle_gy=0.0;
floatQ_angle=0.9;
floatQ_gyro=0.001;
floatR_angle=0.5;
floatdt=0.01;/*dt为kalman滤波器采样时间;*/
charC_0=1;
floatQ_bias,Angle_err;
floatPCt_0=0.0,PCt_1=0.0,E=0.0;
floatK_0=0.0,K_1=0.0,t_0=0.0,t_1=0.0;
floatPdot[4]={0,0,0,0};
floatPP[2][2]={{1,0},{0,1}};
/*卡尔曼滤波函数,具体实现可以参考网上资料,也可以使用其它滤波算法*/
voidKalman_Filter（floatAccel,floatGyro）
{
Angle+=（Gyro-Q_bias）*dt;
Pdot[0]=Q_angle-PP[0][1]-PP[1][0];
Pdot[1]=-PP[1][1];
Pdot[2]=-PP[1][1];
Pdot[3]=Q_gyro;
PP[0][0]+=Pdot[0]*dt;
PP[0][1]+=Pdot[1]*dt;
PP[1][0]+=Pdot[2]*dt;
PP[1][1]+=Pdot[3]*dt;
Angle_err=Accel-Angle;
PCt_0=C_0*PP[0][0];
PCt_1=C_0*PP[1][0];
E=R_angle+C_0*PCt_0;
if（E!
=0）
{
K_0=PCt_0/E;
K_1=PCt_1/E;
}
t_0=PCt_0;
t_1=C_0*PP[0][1];
PP[0][0]-=K_0*t_0;
PP[0][1]-=K_0*t_1;
PP[1][0]-=K_1*t_0;
PP[1][1]-=K_1*t_1;
Angle+=K_0*Angle_err;
Q_bias+=K_1*Angle_err;
}
voidloop（）{
//readrawaccel/gyromeasurementsfromdevice
doubleax_angle=0.0;
doublegx_angle=0.0;
unsignedlongtime=0;
unsignedlongmictime=0;
staticunsignedlongpretime=0;
floatgyro=0.0;
if（pretime==0）
{
pretime=millis（）;
return;
}
mictime=millis（）;
accelgyro.getMotion6（&ax,&ay,&az,&gx,&gy,&gz）;
/*加速度量程范围设置2g16384LSB/g
*计算公式：
*前边已经推导过这里再列出来一次
*x是小车倾斜的角度,y是加速度计读出的值
*sinx=0.92*3.14*x/180=y/16384
*x=180*y/（0.92*3.14*16384）=
*/
ax-=AX_ZERO;
ax_angle=ax/262;
/*陀螺仪量程范围设置250131LSB//s
*陀螺仪角度计算公式:
*小车倾斜角度是gx_angle,陀螺仪读数是y,时间是dt
*gx_angle+=（y/（131*1000））*dt
*/
gy-=GX_ZERO;
time=mictime-pretime;
gyro=gy/131.0;
gx_angle=gyro*time;
gx_angle=gx_angle/1000.0;
total_angle-=gx_angle;
dt=time/1000.0;
Kalman_Filter（ax_angle,gyro）;
Serial.print（ax_angle）;Serial.print（","）;
Serial.print（total_angle）;Serial.print（","）;
Serial.println（Angle）;
pretime=mictime;
}
Mpu6050取值范围调整的方法,修改MPU6050库中的mpu6050.cpp,修改initialize函数中标成红色的行。
voidMPU6050:
:
initialize（）{
setClockSource（MPU6050_CLOCK_PLL_XGYRO）;
setFullScaleGyroRange（MPU6050_GYRO_FS_250）;/*陀螺仪取值范围*/
setFullScaleAccelRange（MPU6050_ACCEL_FS_2）;/*加速度计取值范围*/
setSleepEnabled（false）;//thankstoJackElstonforpointingthisoneout!
}
下别贴出分析出的图供大家参考,图是使用SericalChart，通过收到的串口数据绘制的。
红色的是陀螺仪读出的角度，绿色的是加速度计读出的角度。
红色是陀螺仪读出的角度，绿色是加速度计读出的角度，黄色是卡尔曼滤波融合后的角度，由于参数调整的不是很好，可以看出黄色曲线有滞后的现象。
升级会员 