卡尔曼滤波用的是非矩阵还是矩阵的?

subingkiun

在网上找到的非矩阵的卡尔曼滤波的程序,非矩阵的,有人用他直立过么?

volatile float QingJiao = 0;
volatile float Gyro_Data = 0;
void Kalman(void)
{
   float Q =1,R = 300;
   static float RealData = 0,RealData_P = 0;
   float NowData = 0,NowData_P = 0;
   float Kg = 0,gyroscope_rate = 0,accelerometer_angle=0;
   float Acc_x = 0,Acc_z = 0, Gyro = 0;
   static float his_acc = 0.0,his_accx = 0.0,his_accz = 0.0;
   
   while(!ATD0STAT0_SCF);   //等待转换完成
        
    Acc_x = (float)ATD0DR1;
    Acc_z = (float)ATD0DR2;
    Gyro  = (float)ATD0DR0;
   
   if(Gyro > 4090) Gyro += 1000;
   else
   if(Gyro > 4084) Gyro += 500;
   
   if(Gyro < 55)   Gyro -= 1000;
   else
   if(Gyro < 60)   Gyro -= 500;
      
    Acc_x = Acc_x - 2342.0;
    Acc_z = Acc_z - 2076.0;
    Gyro  = Gyro  - 2048.0;
    Gyro_Data = Gyro;
   
    OutData[0] = Gyro_Data;
    accelerometer_angle = atan2f(-Acc_x,Acc_z);
    OutData[1] = accelerometer_angle*1000;

    gyroscope_rate = Gyro*0.0023;         //参考电压3.3v 12位ADC 放大9.1倍 enc-03 0.67mv/deg./sec.
   
                                          //(3300/4096)/(0.67*9.1)*(3.14/180) =  0.0023
   
    NowData = RealData + gyroscope_rate*0.01;                 
    //1.预估计 X(k|k-1) = A(k,k-1)*X(k-1|k-1) + B(k)*u(k)
    NowData_P = sqrt(Q*Q+RealData_P*RealData_P);            
    //2.计算预估计协方差矩阵   P(k|k-1) = A(k,k-1)*P(k-1|k-1)*A(k,k-1)'+Q(k)
    Kg = sqrt(NowData_P*NowData_P/(NowData_P*NowData_P+R*R));
    //3.计算卡尔曼增益矩阵 K(k) = P(k|k-1)*H(k)' / (H(k)*P(k|k-1)*H(k)' + R(k))
    RealData = NowData + Kg*(accelerometer_angle - NowData);
    //4.更新估计 X(k|k) = X(k|k-1)+K(k)*(Z(k)-H(k)*X(k|k-1))
    RealData_P = sqrt((1-Kg)*NowData_P*NowData_P);           
    //5.计算更新后估计协防差矩阵 P(k|k) =（I-K(k)*H(k)）*P(k|k-1)
   
    QingJiao =  RealData;  
   
   
    OutData[2] = QingJiao*1000;
    OutPut_Data();
}
他用的AD转换是12位精度的。但是输出的波形是这样的
那我要调什么参数？
是不是陀螺仪的比例大了？
真心求指教！！！
红色的是陀螺仪
黄色是加速度的角度
蓝色的是融合后的。

帅气、灰太郎

   Acc_x = (float)ATD0DR1;
    Acc_z = (float)ATD0DR2;
    Gyro  = (float)ATD0DR0;   

请问下 Acc_x,Acc_z是x轴和z轴加速度吗？
还有平衡车只要一个加速度计 就行了。。。。。
不理解。。。。。

huangjianwu

Idenktsie

这个float Acc_x = 0,Acc_z = 0, Gyro = 0;是定义什么？？

Idenktsie

Acc_x = Acc_x - 2342.0;
     Acc_z = Acc_z - 2076.0;
     Gyro  = Gyro  - 2048.0;
     Gyro_Data = Gyro;
这些数据怎么来的啊？？

huayuc19

还有看LZ 贴的程序 用到了X轴。想问下 呢个求atan2f求的角度有什么用？

huayuc19

想请教下楼主 卡尔曼里面的dt怎么改？

subingkiun

jeyran 发表于 2013-3-16 18:01
村田的哪一款

是哪家公司做的模块吧?
我们用的是自己做的

jeyran

村田的哪一款

subingkiun

jeyran 发表于 2013-3-16 17:55
楼主用哪一款的陀螺仪?

官方不是规定了用村田公司的么?