標(biāo)題: stm32f103rct6卡爾曼濾波源碼 [打印本頁(yè)]
作者: xcj666    時(shí)間: 2018-5-5 23:00
標(biāo)題: stm32f103rct6卡爾曼濾波源碼
卡爾曼濾波源碼
以stm32f103rct6為例
測(cè)試串口通信數(shù)據(jù)
MPU數(shù)據(jù)進(jìn)行處理
檢測(cè)X軸加速度，Y軸角速度
參考網(wǎng)上例程

單片機(jī)源程序如下:
  1. #include "stm32f10x.h"
  2. #include <stdio.h>
  3. #include "delay.h"
  4. #include "sys.h"
  5. #include "Usart.h"
  6. #include "MPU6050.h"
  7. //******角度參數(shù)************
  8. float Gyro_y;        //Y軸陀螺儀數(shù)據(jù)暫存
  9. float Angle_gy;      //由角速度計(jì)算的傾斜角度
  10. float Accel_x;             //X軸加速度值暫存
  11. float Angle_ax;      //由加速度計(jì)算的傾斜角度
  12. float Angle;         //最終測(cè)量角度
  13. void TIM2_Configuration(void)
  14. {
  15.   TIM_TimeBaseInitTypeDef  TIM_TimeBaseStructure;//定義一個(gè)定時(shí)器結(jié)構(gòu)體變量

  17.   RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2, ENABLE);//使能定時(shí)器2
  18.   
  19.   TIM_TimeBaseStructure.TIM_Period = (10000 - 1);        //計(jì)數(shù)100次，因?yàn)閺?開(kāi)始，所以減1
  20.   TIM_TimeBaseStructure.TIM_Prescaler = (72 - 1);        //時(shí)鐘72分頻，因?yàn)?不分頻，所以減1
  21.   TIM_TimeBaseStructure.TIM_ClockDivision = 0;        // 使用的采樣頻率之間的分頻比例
  22.    TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up;        //向上計(jì)數(shù)
  23.   TIM_TimeBaseInit(TIM2, &TIM_TimeBaseStructure);        //初始化定時(shí)器2

  25.   TIM_ClearITPendingBit(TIM2, TIM_IT_Update);                //清除定時(shí)器2中斷標(biāo)志位
  26.   TIM_ITConfig(TIM2, TIM_IT_Update, ENABLE);                //打開(kāi)定時(shí)器2中斷

  28.   TIM_Cmd(TIM2, ENABLE);  //計(jì)數(shù)器使能，開(kāi)始計(jì)數(shù)
  29. }

  31. void NVIC_Configuration(void)
  32. {
  33.   NVIC_InitTypeDef NVIC_InitStructure;
  34.   
  35.   NVIC_SetVectorTable(NVIC_VectTab_FLASH, 0x0000);
  36.   NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);

  38. /********************注冊(cè)定時(shí)器2中斷處理函數(shù)***********************/

  40.   NVIC_InitStructure.NVIC_IRQChannel = TIM2_IRQn;
  41.   NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 0;
  42.   NVIC_InitStructure.NVIC_IRQChannelSubPriority = 0;
  43.   NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE;
  44.   NVIC_Init(&NVIC_InitStructure);

  46. /******************注冊(cè)串口1中斷服務(wù)函數(shù)************************/
  47.         NVIC_InitStructure.NVIC_IRQChannel = USART1_IRQn;//配置串口中斷
  48.         NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 1;//占先式優(yōu)先級(jí)設(shè)置為0
  49.         NVIC_InitStructure.NVIC_IRQChannelSubPriority = 1; //副優(yōu)先級(jí)設(shè)置為0
  50.         NVIC_InitStructure.NVIC_IRQChannelCmd = DISABLE;//中斷禁止
  51.         NVIC_Init(&NVIC_InitStructure);//中斷初始化


  54. }

  56. /*************卡爾曼濾波*********************************/
  57. void Kalman_Filter(float Accel,float Gyro)               
  58. {
  59.                
  60.         static const float Q_angle=0.001;  
  61.         static const float Q_gyro=0.003;
  62.         static const float R_angle=0.5;
  63.         static const float dt=0.01;                          //dt為kalman濾波器采樣時(shí)間;
  64.         static const char  C_0 = 1;
  65.         static float Q_bias, Angle_err;
  66.         static float PCt_0, PCt_1, E;
  67.         static float K_0, K_1, t_0, t_1;
  68.         static float Pdot[4] ={0,0,0,0};
  69.         static float PP[2][2] = { { 1, 0 },{ 0, 1 } };
  70.        
  71.         Angle+=(Gyro - Q_bias) * dt; //先驗(yàn)估計(jì)

  73.         Pdot[0]=Q_angle - PP[0][1] - PP[1][0]; // Pk-先驗(yàn)估計(jì)誤差協(xié)方差的微分

  75.         Pdot[1]= -PP[1][1];
  76.         Pdot[2]= -PP[1][1];
  77.         Pdot[3]=Q_gyro;
  78.        
  79.         PP[0][0] += Pdot[0] * dt;   // Pk-先驗(yàn)估計(jì)誤差協(xié)方差微分的積分
  80.         PP[0][1] += Pdot[1] * dt;   // =先驗(yàn)估計(jì)誤差協(xié)方差
  81.         PP[1][0] += Pdot[2] * dt;
  82.         PP[1][1] += Pdot[3] * dt;
  83.                
  84.         Angle_err = Accel - Angle;        //zk-先驗(yàn)估計(jì)
  85.        
  86.         PCt_0 = C_0 * PP[0][0];
  87.         PCt_1 = C_0 * PP[1][0];
  88.        
  89.         E = R_angle + C_0 * PCt_0;
  90.        
  91.         K_0 = PCt_0 / E;
  92.         K_1 = PCt_1 / E;
  93.        
  94.         t_0 = PCt_0;
  95.         t_1 = C_0 * PP[0][1];

  97.         PP[0][0] -= K_0 * t_0;                 //后驗(yàn)估計(jì)誤差協(xié)方差
  98.         PP[0][1] -= K_0 * t_1;
  99.         PP[1][0] -= K_1 * t_0;
  100.         PP[1][1] -= K_1 * t_1;
  101.                
  102.         Angle        += K_0 * Angle_err;         //后驗(yàn)估計(jì)
  103.         Q_bias        += K_1 * Angle_err;         //后驗(yàn)估計(jì)
  104.         Gyro_y   = Gyro - Q_bias;         //輸出值(后驗(yàn)估計(jì))的微分=角速度
  105. }
  106. void Angle_Calculate(void)
  107. {
  108.         static        uint8_t DataBuffer[2];        //數(shù)據(jù)緩存
  109. /****************************加速度****************************************/
  110.         I2C_ReadBuffer(DataBuffer, ACCEL_XOUT_H, 2);
  111.         Accel_x  = (short)((DataBuffer[0]<<8)+DataBuffer[1]);          //讀取X軸加速度
  112.         Angle_ax = (Accel_x +220) /16384;   //去除零點(diǎn)偏移,計(jì)算得到角度（弧度）
  113.         Angle_ax = Angle_ax*1.2*180/3.14;     //弧度轉(zhuǎn)換為度,

  115. /****************************角速度****************************************/
  116.         I2C_ReadBuffer(DataBuffer, GYRO_YOUT_H, 2);
  117.         Gyro_y = (short)((DataBuffer[0]<<8)+DataBuffer[1]);              //靜止時(shí)角速度Y軸輸出為-18左右
  118.         Gyro_y = (Gyro_y + 18)/16.4;         //去除零點(diǎn)偏移，計(jì)算角速度值
  119. //                Angle_gy = Angle_gy + Gyro_y*0.01;  //角速度積分得到傾斜角度.       

  121. /***************************卡爾曼濾波+角度融合*************************************/
  122.         Kalman_Filter(Angle_ax,Gyro_y);       //卡爾曼濾波計(jì)算傾角

  124. /*******************************互補(bǔ)濾波******************************************/

  126. /*補(bǔ)償原理是取當(dāng)前傾角和加速度獲
  127.         得傾角差值進(jìn)行放大，然后與陀螺
  128.         儀角速度疊加后再積分，從而使傾
  129.         角最跟蹤為加速度獲得的角度0.5
  130.         為放大倍數(shù)，可調(diào)節(jié)補(bǔ)償度;
  131.         0.01為系統(tǒng)周期10ms       
  132. */       
  133. //        Angle = Angle + (((Angle_ax-Angle)*0.5 + Gyro_y)*0.01);*/
  134.                                                                                                                                          
  135. }
  136. int main(void)
  137. {
  138.         Stm32_Clock_Init(6);//系統(tǒng)時(shí)鐘設(shè)置為外部晶振，9倍頻
  139.         delay_init(72);//系統(tǒng)SysTick初始化
  140.         Usart_Configuration();
  141.         TIM2_Configuration();
  142.         NVIC_Configuration();
  143.         I2C_Congiguration();
  144.         MPU6050_Init();
  145.         while (1)
  146.         {
  147.                 printf("%f        %f        ",Angle,Gyro_y);
  148.                 USART1_Send_Enter();       
  149.                 delay_ms(100);
  150.         }
  151. }
  152. /****************定時(shí)器2中斷服務(wù)函數(shù)***************************************/
  153. void TIM2_IRQHandler(void)
  154. {
  155.         if (TIM_GetITStatus(TIM2, TIM_IT_Update) != RESET)
  156.   {
  157.                 TIM_ClearITPendingBit(TIM2, TIM_IT_Update);
  158.                 Angle_Calculate();
  159.         }
  160. }
復(fù)制代碼

所有資料51hei提供下載:
卡爾曼濾波.rar (427.33 KB, 下載次數(shù): 68)






歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1