极客工坊

 找回密码
 注册

QQ登录

只需一步,快速开始

查看: 61726|回复: 39

Arduino教程:MPU6050的数据获取、分析与处理

  [复制链接]
发表于 2015-6-27 11:32:52 | 显示全部楼层 |阅读模式
本帖最后由 Devymex 于 2015-6-27 16:25 编辑

摘要

MPU6050是一种非常流行的空间运动传感器芯片,可以获取器件当前的三个加速度分量和三个旋转角速度。由于其体积小巧,功能强大,精度较高,不仅被广泛应用于工业,同时也是航模爱好者的神器,被安装在各类飞行器上驰骋蓝天。

随着Arduino开发板的普及,许多朋友希望能够自己制作基于MPU6050的控制系统,但由于缺乏专业知识而难以上手。此外,MPU6050的数据是有较大噪音的,若不进行滤波会对整个控制系统的精准确带来严重影响。

MPU6050芯片内自带了一个数据处理子模块DMP,已经内置了滤波算法,在许多应用中使用DMP输出的数据已经能够很好的满足要求。关于如何获取DMP的输出数据,我将在以后的文章中介绍。本文将直接面对原始测量数据,从连线、芯片通信开始一步一步教你如何利用Arduino获取MPU6050的数据并进行卡尔曼滤波,最终获得稳定的系统运动状态。

一、Arduino与MPU-6050的通信

为避免纠缠于电路细节,我们直接使用集成的MPU6050模块。MPU6050的数据接口用的是I2C总线协议,因此我们需要Wire程序库的帮助来实现Arduino与MPU6050之间的通信。请先确认你的Arduino编程环境中已安装Wire库。
Wire库的官方文档(http://www.arduino.cc/en/Reference/Wire)中指出:在UNO板子上,SDA接口对应的是A4引脚,SCL对应的是A5引脚。MPU6050需要5V的电源,可由UNO板直接供电。按照下图连线。

Arduino2MPU6050

Arduino2MPU6050

(紫色线是中断线,这里用不到,可以不接)

MPU6050的数据写入和读出均通过其芯片内部的寄存器实现,这些寄存器的地址都是1个字节,也就是8位的寻址空间,其寄存器的详细列表说明书请点击下载:https://www.olimex.com/Products/Modules/Sensors/MOD-MPU6050/resources/RM-MPU-60xxA_rev_4.pdf

1.1 将数据写入MPU-6050

在每次向器件写入数据前要先打开Wire的传输模式,并指定器件的总线地址,MPU6050的总线地址是0x68(AD0引脚为高电平时地址为0x69)。然后写入一个字节的寄存器起始地址,再写入任意长度的数据。这些数据将被连续地写入到指定的起始地址中,超过当前寄存器长度的将写入到后面地址的寄存器中。写入完成后关闭Wire的传输模式。下面的示例代码是向MPU6050的0x6B寄存器写入一个字节0。
  1. Wire.beginTransmission(0x68); //开启MPU6050的传输
  2. Wire.write(0x6B); //指定寄存器地址
  3. Wire.write(0); //写入一个字节的数据
  4. Wire.endTransmission(true); //结束传输,true表示释放总线
复制代码


1.2 从MPU-6050读出数据

读出和写入一样,要先打开Wire的传输模式,然后写一个字节的寄存器起始地址。接下来将指定地址的数据读到Wire库的缓存中,并关闭传输模式。最后从缓存中读取数据。下面的示例代码是从MPU6050的0x3B寄存器开始读取2个字节的数据:
  1. Wire.beginTransmission(0x68); //开启MPU6050的传输
  2. Wire.write(0x3B); //指定寄存器地址
  3. Wire.requestFrom(0x68, 2, true); //将输据读出到缓存
  4. Wire.endTransmission(true); //关闭传输模式
  5. int val = Wire.read() << 8 | Wire.read(); //两个字节组成一个16位整数
复制代码


1.3 具体实现

通常应当在setup函数中对Wire库进行初始化:
  1. Wire.begin();
复制代码

在对MPU6050进行各项操作前,必须启动该器件,向它的0x6B写入一个字节0即可启动。通常也是在setup函数完成,代码见1.1节。

二、 MPU6050的数据格式

我们感兴趣的数据位于0x3B到0x48这14个字节的寄存器中。这些数据会被动态更新,更新频率最高可达1000HZ。下面列出相关寄存器的地址,数据的名称。注意,每个数据都是2个字节。

    0x3B,加速度计的X轴分量ACC_X
    0x3D,加速度计的Y轴分量ACC_Y
    0x3F,加速度计的Z轴分量ACC_Z
    0x41,当前温度TEMP
    0x43,绕X轴旋转的角速度GYR_X
    0x45,绕Y轴旋转的角速度GYR_Y
    0x47,绕Z轴旋转的角速度GYR_Z

MPU6050芯片的座标系是这样定义的:令芯片表面朝向自己,将其表面文字转至正确角度,此时,以芯片内部中心为原点,水平向右的为X轴,竖直向上的为Y轴,指向自己的为Z轴。见下图:
a980247c82e3dba5ff39e3c6bd896ed4_b.jpg

我们只关心加速度计和角速度计数据的含义,下面分别介绍。

2.1 加速度计

加速度计的三轴分量ACC_X、ACC_Y和ACC_Z均为16位有符号整数,分别表示器件在三个轴向上的加速度,取负值时加速度沿座标轴负向,取正值时沿正向。

三个加速度分量均以重力加速度g的倍数为单位,能够表示的加速度范围,即倍率可以统一设定,有4个可选倍率:2g、4g、8g、16g。以ACC_X为例,若倍率设定为2g(默认),则意味着ACC_X取最小值-32768时,当前加速度为沿X轴正方向2倍的重力加速度;若设定为4g,取-32768时表示沿X轴正方向4倍的重力加速度,以此类推。显然,倍率越低精度越好,倍率越高表示的范围越大,这要根据具体的应用来设定。

我们用f表示倍率,f=0为2g,f=3为16g,设定加速度倍率的代码如下:

  1. Wire.beginTransmission(0x68); //开启MPU-6050的传输
  2. Wire.write(0x1C); //加速度倍率寄存器的地址
  3. Wire.requestFrom(0x68, 1, true); //先读出原配置
  4. unsigned char acc_conf = Wire.read();
  5. acc_conf = ((acc_conf & 0xE7) | (f << 3));
  6. Wire.write(acc_conf);
  7. Wire.endTransmission(true); //结束传输,true表示释放总线
复制代码


再以ACC_X为例,若当前设定的加速度倍率为4g,那么将ACC_X读数换算为加速度的公式为: equation.png ,g可取当地重力加速度。

2.2 角速度计

绕X、Y和Z三个座标轴旋转的角速度分量GYR_X、GYR_Y和GYR_Z均为16位有符号整数。从原点向旋转轴方向看去,取正值时为顺时针旋转,取负值时为逆时针旋转。

三个角速度分量均以“度/秒”为单位,能够表示的角速度范围,即倍率可统一设定,有4个可选倍率:250度/秒、500度/秒、1000度/秒、2000度/秒。以GYR_X为例,若倍率设定为250度/秒,则意味着GYR取正最大值32768时,当前角速度为顺时针250度/秒;若设定为500度/秒,取32768时表示当前角速度为顺时针500度/秒。显然,倍率越低精度越好,倍率越高表示的范围越大。

我们用f表示倍率,f=0为250度/秒,f=3为2000度/秒,除角速度倍率寄存器的地址为0x1B之外,设定加速度倍率的代码与2.1节代码一致。

以GYR_X为例,若当前设定的角速度倍率为1000度/秒,那么将GRY_X读数换算为角速度(顺时针)的公式为: equation.png

三、运动数据分析

在读取加速度计和角速度计的数据并换算为物理值后,根据不同的应用,数据有不同的解译方式。本章将以飞行器运动模型为例,根据加速度和角速度来算出当前的飞行姿态。

3.1 加速度计模型

我们可以把加速度计想象为一个正立方体盒子里放着一个球,这个球被弹簧固定在立方体的中心。当盒子运动时,根据假想球的位置即可算出当前加速度的值。想象如果在太空中,盒子没有任何受力时,假想球将处于正中心的位置,三个轴的加速度均为0。见下图:
4947f925172461770a3b067d2a66412b_b.jpg

如果我们给盒子施加一个水平向左的力,那么显然盒子就会有一个向左的加速度,此时盒内的假想球会因为惯性作用贴向盒内的右侧面。如下图所示:

05.png

为了保证数据的物理意义,MPU6050的加速度计是以假想球在三轴上座标值的相反数作为三个轴的加速度值。当假想球的位置偏向一个轴的正向时,该轴的加速度读数为负值,当假想球的位置偏向一个轴的负向时,该轴的加速度读数为正值。

根据以上分析,当我们把MPU6050芯片水平放于地方,芯片表面朝向天空,此时由于受到地球重力的作用, 假想球的位置偏向Z轴的负向,因此Z轴的加速度读数应为正,且在理想情况下应为g。注意,此加速度的物理意义并不是重力加速度,而是自身运动的加速度,可以这样理解:正因为其自身运动的加速度与重力加速度大小相等方向相反,芯片才能保持静止。

3.2 Roll-pitch-yaw模型与姿态计算

表示飞行器当前飞行姿态的一个通用模型就是建立下图所示坐标系,并用Roll表示绕X轴的旋转,Pitch表示绕Y轴的旋转,Yaw表示绕Z轴的旋转。

bb0b0678c9b8ebea10e41331e84d9533_b.jpg

由于MPU6050可以获取三个轴向上的加速度,而地球重力则是长期存在且永远竖直向下,因此我们可以根据重力加速度相对于芯片的指向为参考算得当前姿态。

为方便起见,我们让芯片正面朝下固定在上图飞机上,且座标系与飞机的坐标系完全重合,以三个轴向上的加速度为分量,可构成加速度向量a(x,y,z)。假设当前芯片处于匀速直线运动状态,那么a应垂直于地面上向,即指向Z轴负方向,模长为|a|=g=\sqrt{x^2+y^2+z^2}(与重力加速度大小相等,方向相反,见3.1节)。若芯片(座标系)发生旋转,由于加速度向量a仍然竖直向上,所以Z轴负方向将不再与a重合。见下图。

dd818cd0ef514049fc92614d2adac435_b.jpg

为了方便表示,上图坐标系的Z轴正方向(机腹以及芯片正面)向下,X轴正方向(飞机前进方向)向右。此时芯片的Roll角φ(黄色)为加速度向量与其在XZ平面上投影(x,0,z)的夹角,Pitch角ω(绿色)与其在YZ平面上投影(0,y,z)的夹角。求两个向量的夹角可用点乘公式: equation.png ,简单推导可得:
equation.png 以及 equation.png
注意,因为arccos函数只能返回正值角度,因此还需要根据不同情况来取角度的正负值。当y值为正时,Roll角要取负值,当x轴为负时,Pitch角要取负值。

3.4 Yaw角的问题

因为没有参考量,所以无法求出当前的Yaw角的绝对角度,只能得到Yaw的变化量,也就是角速度GYR_Z。当然,我们可以通过对GYR_Z积分的方法来推算当前Yaw角(以初始值为准),但由于测量精度的问题,推算值会发生漂移,一段时间后就完全失去意义了。然而在大多数应用中,比如无人机,只需要获得GRY_Z就可以了。

如果必须要获得绝对的Yaw角,那么应当选用MPU9250这款九轴运动跟踪芯片,它可以提供额外的三轴罗盘数据,这样我们就可以根据地球磁场方向来计算Yaw角了,具体方法此处不再赘述。

四、数据处理与实现

MPU6050芯片提供的数据夹杂有较严重的噪音,在芯片处理静止状态时数据摆动都可能超过2%。除了噪音,各项数据还会有偏移的现象,也就是说数据并不是围绕静止工作点摆动,因此要先对数据偏移进行校准 ,再通过滤波算法消除噪音。

4.1 校准

校准是比较简单的工作,我们只需要找出摆动的数据围绕的中心点即可。我们以GRY_X为例,在芯片处理静止状态时,这个读数理论上讲应当为0。但它往往会存在偏移量,比如我们以10ms的间隔读取了10个值如下:
  1. -158.4, -172.9, -134.2, -155.1, -131.2, -146.8, -173.1, -188.6, -142.7, -179.5
复制代码
这10个值的均值,也就是这个读数的偏移量为-158.25。在获取偏移量后,每次的读数都减去偏移量就可以得到校准后的读数了。当然这个偏移量只是估计值,比较准确的偏移量要对大量的数据进行统计才能获知,数据量越大越准,但统计的时间也就越慢。一般校准可以在每次启动系统时进行,那么你应当在准确度和启动时间之间做一个权衡。

三个角速度读数GYR_X、GYR_Y和GYR_Z均可通过统计求平均的方法来获得,但三个加速度分量就不能这样简单的完成了,因为芯片静止时的加速度并不为0。

加速度值的偏移来自两个方面,一是由于芯片的测量精度,导至它测得的加速度向量并不垂直于大地;二是芯片在整个系统(如无人机)上安装的精度是有限的,系统与芯片的座标系很难达到完美重合。前者我们称为读数偏移,后者我们称为角度偏移。因为读数和角度之间是非线性关系,所以要想以高精度进行校准必须先单独校准读数偏移,再把芯片固定在系统中后校准角度偏移。然而,由于校准角度偏移需要专业设备,且对于一般应用来说,两步校准带来的精度提升并不大,因此通常只进行读数校准即可。下面介绍读数校准的方法。我们还3.2节的飞机为例,分以下几个步骤:
  • 首先要确定飞机的坐标系,对于多轴飞行器来说这非常重要。如果坐标系原点的位置或坐标轴的方向存在较大偏差,将会给后面的飞控造成不良影响。
  • 在确定了飞机的坐标系后,为了尽量避免读数偏移带来的影响,首先将MPU6050牢牢地固定在飞机上,并使二者座标系尽可能的重合。当然把Z轴反过来装也是可以的,就是需要重新推算一套角度换算公式。
  • 将飞机置于水平、坚固的平面上,并充分预热。对于多轴无人机而言,空中悬停时的XY平面应当平行于校准时的XY平面。此时,我们认为芯片的加速度方向应当与Z轴负方向重合,且加速度向量的模长为g,因此ACC_X和ACC_Y的理论值应为0,ACC_Z的理论值应为-16384(假设我们设定2g的倍率,1g的加速度的读数应为最大值-32768的一半)。
  • 由于ACC_X和ACC_Y的理论值应为0,与角速度量的校准类似,这两个读数偏移量可用统计均值的方式校准。ACC_Z则需要多一步处理,即在统计偏移量的过程中,每次读数都要加上16384,再进行统计均值校准。

4.2 卡尔曼滤波

对于夹杂了大量噪音的数据,卡尔曼滤波器的效果无疑是最好的。如果不想考虑算法细节,可以直接使用Arduino的Klaman Filter库完成。在我们的模型中,一个卡尔曼滤波器接受一个轴上的角度值、角速度值以及时间增量,估计出一个消除噪音的角度值。跟据当前的角度值和上一轮估计的角度值,以及这两轮估计的间隔时间,我们还可以反推出消除噪音的角速度。

实现代码见4.3节。下面介绍卡尔曼滤波算法细节,不感兴趣的可跳过。

(想看的人多了再写)

4.3 实现代码

以下代码在Arduino软件1.65版本中编译、烧写以及测试通过。
  1. // 本代码版权归Devymex所有,以GNU GENERAL PUBLIC LICENSE V3.0发布
  2. // http://www.gnu.org/licenses/gpl-3.0.en.html
  3. // 相关文档参见作者于知乎专栏发表的原创文章:
  4. // http://zhuanlan.zhihu.com/devymex/20082486

  5. //连线方法
  6. //MPU-UNO
  7. //VCC-VCC
  8. //GND-GND
  9. //SCL-A5
  10. //SDA-A4
  11. //INT-2 (Optional)

  12. #include <Kalman.h>
  13. #include <Wire.h>
  14. #include <Math.h>

  15. float fRad2Deg = 57.295779513f; //将弧度转为角度的乘数
  16. const int MPU = 0x68; //MPU-6050的I2C地址
  17. const int nValCnt = 7; //一次读取寄存器的数量

  18. const int nCalibTimes = 1000; //校准时读数的次数
  19. int calibData[nValCnt]; //校准数据

  20. unsigned long nLastTime = 0; //上一次读数的时间
  21. float fLastRoll = 0.0f; //上一次滤波得到的Roll角
  22. float fLastPitch = 0.0f; //上一次滤波得到的Pitch角
  23. Kalman kalmanRoll; //Roll角滤波器
  24. Kalman kalmanPitch; //Pitch角滤波器

  25. void setup() {
  26.   Serial.begin(9600); //初始化串口,指定波特率
  27.   Wire.begin(); //初始化Wire库
  28.   WriteMPUReg(0x6B, 0); //启动MPU6050设备

  29.   Calibration(); //执行校准
  30.   nLastTime = micros(); //记录当前时间
  31. }

  32. void loop() {
  33.   int readouts[nValCnt];
  34.   ReadAccGyr(readouts); //读出测量值
  35.   
  36.   float realVals[7];
  37.   Rectify(readouts, realVals); //根据校准的偏移量进行纠正

  38.   //计算加速度向量的模长,均以g为单位
  39.   float fNorm = sqrt(realVals[0] * realVals[0] + realVals[1] * realVals[1] + realVals[2] * realVals[2]);
  40.   float fRoll = GetRoll(realVals, fNorm); //计算Roll角
  41.   if (realVals[1] > 0) {
  42.     fRoll = -fRoll;
  43.   }
  44.   float fPitch = GetPitch(realVals, fNorm); //计算Pitch角
  45.   if (realVals[0] < 0) {
  46.     fPitch = -fPitch;
  47.   }

  48.   //计算两次测量的时间间隔dt,以秒为单位
  49.   unsigned long nCurTime = micros();
  50.   float dt = (double)(nCurTime - nLastTime) / 1000000.0;
  51.   //对Roll角和Pitch角进行卡尔曼滤波
  52.   float fNewRoll = kalmanRoll.getAngle(fRoll, realVals[4], dt);
  53.   float fNewPitch = kalmanPitch.getAngle(fPitch, realVals[5], dt);
  54.   //跟据滤波值计算角度速
  55.   float fRollRate = (fNewRoll - fLastRoll) / dt;
  56.   float fPitchRate = (fNewPitch - fLastPitch) / dt;

  57. //更新Roll角和Pitch角
  58.   fLastRoll = fNewRoll;
  59.   fLastPitch = fNewPitch;
  60.   //更新本次测的时间
  61.   nLastTime = nCurTime;

  62.   //向串口打印输出Roll角和Pitch角,运行时在Arduino的串口监视器中查看
  63.   Serial.print("Roll:");
  64.   Serial.print(fNewRoll); Serial.print('(');
  65.   Serial.print(fRollRate); Serial.print("),\tPitch:");
  66.   Serial.print(fNewPitch); Serial.print('(');
  67.   Serial.print(fPitchRate); Serial.print(")\n");
  68.   delay(10);
  69. }

  70. //向MPU6050写入一个字节的数据
  71. //指定寄存器地址与一个字节的值
  72. void WriteMPUReg(int nReg, unsigned char nVal) {
  73.   Wire.beginTransmission(MPU);
  74.   Wire.write(nReg);
  75.   Wire.write(nVal);
  76.   Wire.endTransmission(true);
  77. }

  78. //从MPU6050读出一个字节的数据
  79. //指定寄存器地址,返回读出的值
  80. unsigned char ReadMPUReg(int nReg) {
  81.   Wire.beginTransmission(MPU);
  82.   Wire.write(nReg);
  83.   Wire.requestFrom(MPU, 1, true);
  84.   Wire.endTransmission(true);
  85.   return Wire.read();
  86. }

  87. //从MPU6050读出加速度计三个分量、温度和三个角速度计
  88. //保存在指定的数组中
  89. void ReadAccGyr(int *pVals) {
  90.   Wire.beginTransmission(MPU);
  91.   Wire.write(0x3B);
  92.   Wire.requestFrom(MPU, nValCnt * 2, true);
  93.   Wire.endTransmission(true);
  94.   for (long i = 0; i < nValCnt; ++i) {
  95.     pVals[i] = Wire.read() << 8 | Wire.read();
  96.   }
  97. }

  98. //对大量读数进行统计,校准平均偏移量
  99. void Calibration()
  100. {
  101.   float valSums[7] = {0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0};
  102.   //先求和
  103.   for (int i = 0; i < nCalibTimes; ++i) {
  104.     int mpuVals[nValCnt];
  105.     ReadAccGyr(mpuVals);
  106.     for (int j = 0; j < nValCnt; ++j) {
  107.       valSums[j] += mpuVals[j];
  108.     }
  109.   }
  110.   //再求平均
  111.   for (int i = 0; i < nValCnt; ++i) {
  112.     calibData[i] = int(valSums[i] / nCalibTimes);
  113.   }
  114.   calibData[2] += 16384; //设芯片Z轴竖直向下,设定静态工作点。
  115. }

  116. //算得Roll角。算法见文档。
  117. float GetRoll(float *pRealVals, float fNorm) {
  118.   float fNormXZ = sqrt(pRealVals[0] * pRealVals[0] + pRealVals[2] * pRealVals[2]);
  119.   float fCos = fNormXZ / fNorm;
  120.   return acos(fCos) * fRad2Deg;
  121. }

  122. //算得Pitch角。算法见文档。
  123. float GetPitch(float *pRealVals, float fNorm) {
  124.   float fNormYZ = sqrt(pRealVals[1] * pRealVals[1] + pRealVals[2] * pRealVals[2]);
  125.   float fCos = fNormYZ / fNorm;
  126.   return acos(fCos) * fRad2Deg;
  127. }

  128. //对读数进行纠正,消除偏移,并转换为物理量。公式见文档。
  129. void Rectify(int *pReadout, float *pRealVals) {
  130.   for (int i = 0; i < 3; ++i) {
  131.     pRealVals[i] = (float)(pReadout[i] - calibData[i]) / 16384.0f;
  132.   }
  133.   pRealVals[3] = pReadout[3] / 340.0f + 36.53;
  134.   for (int i = 4; i < 7; ++i) {
  135.     pRealVals[i] = (float)(pReadout[i] - calibData[i]) / 131.0f;
  136.   }
  137. }
复制代码

本文为Devymex于知乎原创首发,另发布于极客工坊,未经授权严禁转载,发现必究,欢迎举报,谢谢!
授权知乎日报转载,授权非营利性的公益机构转载,一般营利性机构/平台非特殊情况不予授权。
回复

使用道具 举报

 楼主| 发表于 2015-6-27 16:21:08 | 显示全部楼层
本帖最后由 Devymex 于 2015-6-27 19:30 编辑
wing 发表于 2015-6-27 13:32
这个内容有点深得慢慢看,楼主你的头像是名画么?


不是什么名画,觉得有意思的一幅图。
2000 (3).jpg
回复 支持 0 反对 1

使用道具 举报

发表于 2015-6-27 12:21:11 | 显示全部楼层
谢谢   楼主
回复 支持 反对

使用道具 举报

发表于 2015-6-27 13:32:12 | 显示全部楼层
这个内容有点深得慢慢看,楼主你的头像是名画么?
回复 支持 反对

使用道具 举报

发表于 2015-6-27 16:01:17 | 显示全部楼层
各位工坊达人,我新开了一个贴吧和DIY网站,大家想玩的可以来玩哦,刚开的没什么人
爱自制贴吧:http://tieba.baidu.com/f?kw=%E7%88%B1%E8%87%AA%E5%88%B6&fr=wwwt
52DIY网站:http://www.52diys.com/
回复 支持 反对

使用道具 举报

发表于 2015-7-19 22:42:16 | 显示全部楼层
万能的楼主,您也提到的零漂问题。我的板子也是,MPU6050开始运动之前,零偏是一个值,运动后,零偏马上漂了,逐渐上升。积分误差会很大。 所以不知道,您对这种随时间变化的零漂怎么处理???谢谢
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-7-24 19:31:38 | 显示全部楼层
lfshit 发表于 2015-7-19 22:42
万能的楼主,您也提到的零漂问题。我的板子也是,MPU6050开始运动之前,零偏是一个值,运动后,零偏马上漂了 ...

单靠六轴传感器的数据是不行的,要九轴的才能修偏。
回复 支持 反对

使用道具 举报

发表于 2015-8-5 09:46:18 | 显示全部楼层
我想问一下不用加上MPU6050的头文件吗?读取角度的那些函数是哪个库的?谢谢了
回复 支持 反对

使用道具 举报

发表于 2015-8-8 15:32:24 | 显示全部楼层
库文件可以传一个吗?
回复 支持 反对

使用道具 举报

发表于 2015-8-11 11:00:17 | 显示全部楼层
谢谢楼主,现在我数学水平还达不到,收藏起来以后学习
回复 支持 反对

使用道具 举报

发表于 2015-8-19 15:33:38 | 显示全部楼层
收藏了 ,谢谢
回复 支持 反对

使用道具 举报

发表于 2015-8-19 17:39:26 | 显示全部楼层
Devymex 发表于 2015-6-27 16:21
不是什么名画,觉得有意思的一幅图。

楼主还是换一张图片吧!看着挺瘆人的!!!
回复 支持 反对

使用道具 举报

发表于 2015-8-19 21:09:32 | 显示全部楼层
好东西,介绍的很详细。解决了我的困扰,谢谢
回复 支持 反对

使用道具 举报

发表于 2015-8-20 09:46:09 | 显示全部楼层
为什么看不到啊
回复 支持 反对

使用道具 举报

发表于 2015-8-26 15:40:08 | 显示全部楼层
很详细,顶。求介绍卡尔曼滤波算法细节~
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|联系我们|极客工坊 ( 浙ICP备09023225号 )

GMT+8, 2020-11-27 08:52 , Processed in 0.059871 second(s), 29 queries .

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表