久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標題: stm32f103rct6卡爾曼濾波源碼 [打印本頁]
作者: xcj666    時間: 2018-5-5 23:00
標題: stm32f103rct6卡爾曼濾波源碼
卡爾曼濾波源碼
以stm32f103rct6為例
測試串口通信數據
MPU數據進行處理
檢測X軸加速度，Y軸角速度
參考網上例程

單片機源程序如下:
  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. //******角度參數************
  8. float Gyro_y;        //Y軸陀螺儀數據暫存
  9. float Angle_gy;      //由角速度計算的傾斜角度
  10. float Accel_x;             //X軸加速度值暫存
  11. float Angle_ax;      //由加速度計算的傾斜角度
  12. float Angle;         //最終測量角度
  13. void TIM2_Configuration(void)
  14. {
  15.   TIM_TimeBaseInitTypeDef  TIM_TimeBaseStructure;//定義一個定時器結構體變量

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

  25.   TIM_ClearITPendingBit(TIM2, TIM_IT_Update);                //清除定時器2中斷標志位
  26.   TIM_ITConfig(TIM2, TIM_IT_Update, ENABLE);                //打開定時器2中斷

  28.   TIM_Cmd(TIM2, ENABLE);  //計數器使能，開始計數
  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. /********************注冊定時器2中斷處理函數***********************/

  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. /******************注冊串口1中斷服務函數************************/
  47.         NVIC_InitStructure.NVIC_IRQChannel = USART1_IRQn;//配置串口中斷
  48.         NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 1;//占先式優先級設置為0
  49.         NVIC_InitStructure.NVIC_IRQChannelSubPriority = 1; //副優先級設置為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濾波器采樣時間;
  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; //先驗估計

  73.         Pdot[0]=Q_angle - PP[0][1] - PP[1][0]; // Pk-先驗估計誤差協方差的微分

  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-先驗估計誤差協方差微分的積分
  80.         PP[0][1] += Pdot[1] * dt;   // =先驗估計誤差協方差
  81.         PP[1][0] += Pdot[2] * dt;
  82.         PP[1][1] += Pdot[3] * dt;
  83.                
  84.         Angle_err = Accel - Angle;        //zk-先驗估計
  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;                 //后驗估計誤差協方差
  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;         //后驗估計
  103.         Q_bias        += K_1 * Angle_err;         //后驗估計
  104.         Gyro_y   = Gyro - Q_bias;         //輸出值(后驗估計)的微分=角速度
  105. }
  106. void Angle_Calculate(void)
  107. {
  108.         static        uint8_t DataBuffer[2];        //數據緩存
  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;   //去除零點偏移,計算得到角度（弧度）
  113.         Angle_ax = Angle_ax*1.2*180/3.14;     //弧度轉換為度,

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

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

  124. /*******************************互補濾波******************************************/

  126. /*補償原理是取當前傾角和加速度獲
  127.         得傾角差值進行放大，然后與陀螺
  128.         儀角速度疊加后再積分，從而使傾
  129.         角最跟蹤為加速度獲得的角度0.5
  130.         為放大倍數，可調節補償度;
  131.         0.01為系統周期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);//系統時鐘設置為外部晶振，9倍頻
  139.         delay_init(72);//系統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. /****************定時器2中斷服務函數***************************************/
  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. }
復制代碼

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






歡迎光臨 (http://m.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 中文字幕在线免费播放 | 久久视频免费观看 | 91免费视频 | 国产伦精品一区二区三区视频黑人 | av观看网站| 午夜视频免费观看 | 亚洲www.| 玖玖色资源 | 狠狠艹狠狠干 | 精品欧美一区二区精品久久 | 精品久久久久久久久久久久久 | 欧美精品在线观看 | 国产欧美日韩一区二区三区 | 精品国产乱码久久久久久蜜臀网站 | 四虎成人影视 | 午夜精品影院 | 色综合色综合 | 亚洲一区欧美 | 国产又黄又猛 | 精品久久网站 | 可以免费看的av | 日韩精品视频免费 | 国产欧美精品一区 | 97视频在线 | 日日操视频 | 国产精品免费一区 | 欧美一区二区在线播放 | 中文字幕在线观看一区二区 | 国产一区二区三区久久 | 国产福利网| 欧美日本在线观看 | 欧美大片黄 | 亚洲高清在线观看 | 免费一级黄色 | 日本精品国产 | 日韩城人网站 | av久草| 国产精品成人一区二区网站软件 | 好了av在线 | 真实的国产乱xxxx在线 | 日韩精品在线一区 |