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

 找回密碼
 立即注冊

QQ登錄

只需一步，快速開始

搜索
查看: 4111|回復: 0
打印 上一主題 下一主題
收起左側

stm32f103rct6卡爾曼濾波源碼

[復制鏈接]
跳轉到指定樓層
樓主
ID:323149 發表于 2018-5-5 23:00 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
卡爾曼濾波源碼
以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)


分享到:  QQ好友和群QQ好友和群 QQ空間QQ空間 騰訊微博騰訊微博 騰訊朋友騰訊朋友
收藏收藏2 分享淘帖 頂頂 踩踩
回復

使用道具 舉報

您需要登錄后才可以回帖 登錄 | 立即注冊

本版積分規則

手機版|小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術交流QQ群281945664

Powered by 單片機教程網

快速回復 返回頂部 返回列表
主站蜘蛛池模板： 伊人看片| 亚洲品质自拍视频网站 | 欧美日韩久 | 99re视频在线 | 日韩一三区 | 天天搞天天操 | 成人av激情 | 亚洲一区二区三区视频 | 国产精久久久久久久 | 日韩成人在线网站 | 五月激情婷婷网 | 黑人巨大精品 | 91精品国产91久久久久久密臀 | 视频一区二区三区中文字幕 | 日日干夜夜干 | 狠狠狠| 中文字幕av亚洲精品一部二部 | 国产午夜精品一区二区三区嫩草 | 一区二区免费在线 | 亚洲国产精品一区二区三区 | 天天玩天天干天天操 | 成人深夜小视频 | 久久久黄色 | 午夜免费视频 | 美女久久久久 | 国产高清精品一区二区三区 | 欧美做暖暖视频 | 国产精品久久一区二区三区 | 99re在线免费视频 | 亚洲国产区 | 一区二区在线看 | 丝袜美腿一区二区三区 | 国产精品18hdxxxⅹ在线 | 午夜寂寞影院在线观看 | 日日夜夜影院 | 国产免费xxx | 精品久| 影音先锋中文在线 | 国产精品一区二区久久精品爱微奶 | 国产成人免费网站 | 欧美亚洲视频 |