找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6776|回复: 1
打印 上一主题 下一主题
收起左侧

自己做的双轮自平衡小车

[复制链接]
跳转到指定楼层
楼主
控制器:ATmega16;8MHz;
加速度传感器:MMA2260;陀螺仪:EWTS82;
传感器的融合:卡尔曼滤波;
马达:EN_2342CR(速比64)+双路12脉冲编码器+CD40106对信号整形;
驱动板芯片:CD4001+IR2111+IRF1404(驱动电流可以很大);
制作资料在附件压缩包里面,另外这篇文章http://www.51hei.com/bbs/dpj-47882-1.html 供参考

自己做的自平衡小车,基本达到预期效果,图片(拍摄效果不佳)及文件如下,视频请见:
http://v.youku.com/v_show/id_XMTc3ODg0Mjky.html
http://v.youku.com/v_show/id_XMTgxNTAxOTA4.html
http://v.youku.com/v_show/id_XMTgxNTAzMTI4.html
资料供大家参考;


照片 (原文件名:RII_2.jpg)


照片 (原文件名:RII_1.jpg)


照片 (原文件名:RII_3.jpg)


今天做了取消马达转速传感器的测试,利用PWM的值低通滤波后代替转速信号,也可以稳定控制,但在移动距离偏差较大时不能修正,小范围内和利用传感器闭环控制差不多;源文件如下
  1. //MCU:Mega16;晶振:8MHz;
  2. //PWM:4KHz;滤波器频率:100Hz;系统频率:100Hz;10ms;
  3. //赵国霖:10.05.10;
  4. #include
  5. #include
  6. #include

  7. //#define checkbit(var,bit)     (var&(0x01<<(bit)))     /*定义查询位函数*/
  8. //#define setbit(var,bit)     (var|=(0x01<<(bit)))     /*定义置位函数*/
  9. //#define clrbit(var,bit)     (var&=(~(0x01<<(bit)))) /*定义清零位函数*/


  10. //-------------------------------------------------------
  11. //输出端口初始化
  12. void PORT_initial(void)
  13. {
  14.         DDRA=0B00000000;
  15.         PINA=0X00;
  16.         PORTA=0X00;
  17.         
  18.         DDRB=0B00000000;
  19.         PINB=0X00;
  20.         PORTB=0X00;
  21.         
  22.         DDRC=0B00010000;
  23.         PINC=0X00;
  24.         PORTC=0X00;
  25.         
  26.         DDRD=0B11110010;
  27.         PIND=0X00;
  28.         PORTD=0X00;
  29. }


  30. //-------------------------------------------------------
  31. //定时器1初始化
  32. void T1_initial(void)                                       
  33. {
  34.         TCCR1A|=(1<<com1a1)|(1<<wgm10)|(1<<com1b1); t1:8位快速pwm模式、匹配时清0,top时置位
  35.         TCCR1B|=(1<<wgm12)|(1<<cs11); pwm:8分频:8m="" 8="" 256="4KHz;
  36. }


  37. //-------------------------------------------------------
  38. //定时器2初始化
  39. void T2_initial(void)                        //T2:计数至OCR2时产生中断
  40. {
  41.         OCR2=0X4E;                        //T2:计数20ms(0X9C)10ms(0X4E)时产生中断;
  42.         TIMSK|=(1<<ocie2);
  43.         TCCR2|=(1<<wgm21)|(1<<cs22)|(1<<cs21)|(1<<cs20); ctc模式,1024分频
  44. }


  45. //-------------------------------------------------------
  46. //外部中断初始化
  47. void INT_initial(void)
  48. {
  49.         MCUCR|=(1<<isc01)|(1<<isc00)|(1<<isc11)|(1<<isc10); int0、int1上升沿有效
  50.         GICR|=(1<<int0)|(1<<int1); 外部中断使能
  51. }


  52. //-------------------------------------------------------
  53. //串口初始化;
  54. void USART_initial( void )
  55. {        
  56.         UBRRH = 0X00;
  57.         UBRRL = 51;                //f=8MHz;设置波特率:9600:51;19200:25;
  58.         UCSRB = (1<<rxen)|(1<<txen); 接收器与发送器使能;
  59.         UCSRC = (1<<ursel)|(1<<ucsz0)|(1<<ucsz1); 设置帧格式:="" 8="" 个数据位,="" 1="" 个停止位;
  60.         
  61.         UCSRB|=(1<<rxcie); usart接收中断使能
  62. }


  63. //-------------------------------------------------------
  64. //串口发送数据;
  65. void USART_Transmit( unsigned char data )
  66. {
  67.         while ( !( UCSRA & (1<<udre))); 等待发送缓冲器为空;
  68.         UDR = data;                //将数据放入缓冲器,发送数据;
  69. }


  70. //-------------------------------------------------------
  71. //串口接收数据中断,确定数据输出的状态;
  72. #pragma interrupt_handler USART_Receive_Int:12
  73. static char USART_State;
  74. void USART_Receive_Int(void)
  75. {
  76.         USART_State=UDR;//USART_Receive();
  77. }


  78. //-------------------------------------------------------
  79. //计算LH侧轮速:INT0中断;
  80. //-------------------------------------------------------
  81. static int speed_real_LH;
  82. //-------------------------------------------------------
  83. #pragma interrupt_handler SPEEDLHINT_fun:2
  84. void SPEEDLHINT_fun(void)
  85. {
  86.         if (0==(PINB&BIT(0)))
  87.         {
  88.                 speed_real_LH-=1;
  89.         }
  90.         else
  91.         {
  92.                 speed_real_LH+=1;
  93.         }
  94. }


  95. //-------------------------------------------------------
  96. //计算RH侧轮速,:INT1中断;
  97. //同时将轮速信号统一成前进方向了;
  98. //-------------------------------------------------------
  99. static int speed_real_RH;
  100. //-------------------------------------------------------
  101. #pragma interrupt_handler SPEEDRHINT_fun:3
  102. void SPEEDRHINT_fun(void)
  103. {
  104.         if (0==(PINB&BIT(1)))
  105.         {
  106.                 speed_real_RH+=1;
  107.         }
  108.         else
  109.         {
  110.                 speed_real_RH-=1;
  111.         }
  112. }


  113. //-------------------------------------------------------
  114. //ADport采样:10位,采样基准电压Aref
  115. //-------------------------------------------------------
  116. static int AD_data;
  117. //-------------------------------------------------------
  118. int ADport(unsigned char port)
  119. {
  120.         ADMUX=port;
  121.         ADCSRA|=(1<<aden)|(1<<adsc)|(1<<adps1)|(1<<adps0); 采样频率为8分频;
  122.         while(!(ADCSRA&(BIT(ADIF))));
  123.         AD_data=ADCL;
  124.         AD_data+=ADCH*256;
  125.         AD_data-=512;
  126.         return (AD_data);
  127. }


  128. //*
  129. //-------------------------------------------------------
  130. //Kalman滤波,8MHz的处理时间约1.8ms;
  131. //-------------------------------------------------------
  132. static float angle, angle_dot;                 //外部需要引用的变量
  133. //-------------------------------------------------------
  134. static const float Q_angle=0.001, Q_gyro=0.003, R_angle=0.5, dt=0.01;
  135.                         //注意:dt的取值为kalman滤波器采样时间;
  136. static float P[2][2] = {
  137.                                                         { 1, 0 },
  138.                                                         { 0, 1 }
  139.                                                 };
  140.         
  141. static float Pdot[4] ={0,0,0,0};

  142. static const char C_0 = 1;

  143. static float q_bias, angle_err, PCt_0, PCt_1, E, K_0, K_1, t_0, t_1;
  144. //-------------------------------------------------------
  145. void Kalman_Filter(float angle_m,float gyro_m)                        //gyro_m:gyro_measure
  146. {
  147.         angle+=(gyro_m-q_bias) * dt;
  148.         
  149.         Pdot[0]=Q_angle - P[0][1] - P[1][0];
  150.         Pdot[1]=- P[1][1];
  151.         Pdot[2]=- P[1][1];
  152.         Pdot[3]=Q_gyro;
  153.         
  154.         P[0][0] += Pdot[0] * dt;
  155.         P[0][1] += Pdot[1] * dt;
  156.         P[1][0] += Pdot[2] * dt;
  157.         P[1][1] += Pdot[3] * dt;
  158.         
  159.         
  160.         angle_err = angle_m - angle;
  161.         
  162.         
  163.         PCt_0 = C_0 * P[0][0];
  164.         PCt_1 = C_0 * P[1][0];
  165.         
  166.         E = R_angle + C_0 * PCt_0;
  167.         
  168.         K_0 = PCt_0 / E;
  169.         K_1 = PCt_1 / E;
  170.         
  171.         t_0 = PCt_0;
  172.         t_1 = C_0 * P[0][1];

  173.         P[0][0] -= K_0 * t_0;
  174.         P[0][1] -= K_0 * t_1;
  175.         P[1][0] -= K_1 * t_0;
  176.         P[1][1] -= K_1 * t_1;
  177.         
  178.         
  179.         angle        += K_0 * angle_err;
  180.         q_bias        += K_1 * angle_err;
  181.         angle_dot = gyro_m-q_bias;
  182. }
  183. //*/

  184. /*
  185. //-------------------------------------------------------
  186. //互补滤波
  187. //-------------------------------------------------------
  188. static float angle,angle_dot;                 //外部需要引用的变量
  189. //-------------------------------------------------------        
  190. static float bias_cf;
  191. static const float dt=0.01;
  192. //-------------------------------------------------------
  193. void complement_filter(float angle_m_cf,float gyro_m_cf)
  194. {
  195.         bias_cf*=0.998;                        //陀螺仪零飘低通滤波;500次均值;
  196.         bias_cf+=gyro_m_cf*0.002;
  197.         angle_dot=gyro_m_cf-bias_cf;
  198.         angle=(angle+angle_dot*dt)*0.90+angle_m_cf*0.05;        
  199.         //加速度低通滤波;20次均值;按100次每秒计算,低通5Hz;
  200. }
  201. */        




  202. //-------------------------------------------------------
  203. //AD采样;
  204. //以角度表示;
  205. //加速度计:1.2V=1g=90°;满量程:1.3V~3.7V;
  206. //陀螺仪:0.5V~4.5V=-80°~+80°;满量程5V=200°=256=200°;
  207. //-------------------------------------------------------
  208. static float gyro,acceler;
  209. //-------------------------------------------------------
  210. void AD_calculate(void)
  211. {
  212.         
  213.         acceler=ADport(2)+28;                        //角度校正
  214.         gyro=ADport(3);        
  215.         
  216.         acceler*=0.004069;                //系数换算:2.5/(1.2*512);
  217.         acceler=asin(acceler);
  218.         gyro*=0.00341;                        //角速度系数:(3.14/180)* 100/512=0.01364;        
  219.         
  220.         Kalman_Filter(acceler,gyro);
  221.         //complement_filter(acceler,gyro);
  222. }
  223.         


  224. //-------------------------------------------------------
  225. //PWM输出
  226. //-------------------------------------------------------
  227. void PWM_output (int PWM_LH,int PWM_RH)
  228. {
  229.         if (PWM_LH<0)
  230.         {
  231.                 PORTD|=BIT(6);
  232.                 PWM_LH*=-1;
  233.         }
  234.         else
  235.         {
  236.                 PORTD&=~BIT(6);
  237.         }
  238.         
  239.         if (PWM_LH>252)
  240.         {
  241.                 PWM_LH=252;
  242.         }
  243.         
  244.         
  245.         if (PWM_RH<0)
  246.         {
  247.                 PORTD|=BIT(7);
  248.                 PWM_RH*=-1;
  249.         }
  250.         else
  251.         {
  252.                 PORTD&=~BIT(7);
  253.         }
  254.         
  255.         if (PWM_RH>252)
  256.         {
  257.                 PWM_RH=252;
  258.         }
  259.         
  260.                
  261.         OCR1AH=0;
  262.         OCR1AL=PWM_LH;                        //OC1A输出;
  263.         
  264.         OCR1BH=0;
  265.         OCR1BL=PWM_RH;                        //OC1B输出;
  266.         
  267. }




  268. //-------------------------------------------------------
  269. //计算PWM输出值
  270. //车辆直径:76mm; 12*64pulse/rev; 1m=3216pulses;
  271. //-------------------------------------------------------        
  272. //static int speed_diff,speed_diff_all,speed_diff_adjust;
  273. //static float K_speed_P,K_speed_I;
  274. static float K_voltage,K_angle,K_angle_dot,K_position,K_position_dot;
  275. static float K_angle_AD,K_angle_dot_AD,K_position_AD,K_position_dot_AD;
  276. static float position,position_dot;
  277. static float position_dot_filter;
  278. static float PWM;
  279. static int speed_output_LH,speed_output_RH;
  280. static int Turn_Need,Speed_Need;
  281. //-------------------------------------------------------        
  282. void PWM_calculate(void)        
  283. {
  284.         if ( 0==(~PINA&BIT(1)) )        //左转
  285.         {
  286.                 Turn_Need=-40;
  287.         }
  288.         else if ( 0==(~PINB&BIT(2)) )         //右转
  289.         {
  290.                 Turn_Need=40;
  291.         }
  292.         else        //不转
  293.         {
  294.                 Turn_Need=0;
  295.         }
  296.         
  297.         if ( 0==(~PINC&BIT(0)) )        //前进
  298.         {
  299.                 Speed_Need=-2;
  300.         }
  301.         else if ( 0==(~PINC&BIT(1)) )        //后退
  302.         {
  303.                 Speed_Need=2;
  304.         }
  305.         else        //不动
  306.         {
  307.                 Speed_Need=0;
  308.         }
  309.         
  310.         K_angle_AD=ADport(4)*0.007;
  311.         K_angle_dot_AD=ADport(5)*0.007;
  312.         K_position_AD=ADport(6)*0.007;
  313.         K_position_dot_AD=ADport(7)*0.007;

  314.         position_dot=speed_output_RH*0.04;
  315.         
  316.         position_dot_filter*=0.9;                //车轮速度滤波
  317.         position_dot_filter+=position_dot*0.1;        
  318.         
  319.         position+=position_dot_filter;
  320.         //position+=position_dot;        
  321.         position+=Speed_Need;
  322.         
  323.         if (position<-768)                //防止位置误差过大导致的不稳定
  324.         {
  325.                 position=-768;
  326.         }
  327.         else if  (position>768)
  328.         {
  329.                 position=768;
  330.         }
  331.         
  332.         PWM = K_angle*angle *K_angle_AD + K_angle_dot*angle_dot *K_angle_dot_AD +
  333.                                         K_position*position *K_position_AD + K_position_dot*position_dot_filter *K_position_dot_AD;
  334.         
  335.         
  336.         speed_output_RH = PWM - Turn_Need;
  337.         speed_output_LH = - PWM - Turn_Need ;
  338.         
  339.         /*
  340.         speed_diff=speed_real_RH-speed_real_LH;                        //左右轮速差PI控制;
  341.         speed_diff_all+=speed_diff;
  342.         speed_diff_adjust=(K_speed_P*speed_diff+K_speed_I*speed_diff_all)/2;
  343.         */
  344.                
  345.         PWM_output (speed_output_LH,speed_output_RH);        
  346. }

  347. //-------------------------------------------------------
  348. //定时器2中断处理
  349. //-------------------------------------------------------
  350. static unsigned char temp;
  351. //-------------------------------------------------------
  352. #pragma interrupt_handler T2INT_fun:4
  353. void T2INT_fun(void)                                                
  354. {
  355.         AD_calculate();
  356.         
  357.         PWM_calculate();
  358.         
  359.         if(temp>=4)                        //10ms即中断;每秒计算:100/4=25次;
  360.         {               
  361.                 if (USART_State==0X30)                //ASCII码:0X30代表字符'0'
  362.                 {
  363.                         USART_Transmit(angle*57.3+128);
  364.                         USART_Transmit(angle_dot*57.3+128);
  365.                         USART_Transmit(128);        
  366.                 }
  367.                 else if(USART_State==0X31)                //ASCII码:0X30代表字符'1'
  368.                 {
  369.                         USART_Transmit(speed_output_LH+128);
  370.                         USART_Transmit(speed_output_RH+128);                        
  371.                         USART_Transmit(128);
  372.                 }
  373.                 else if(USART_State==0X32)                //ASCII码:0X30代表字符'2'
  374.                 {
  375.                         USART_Transmit(speed_real_LH+128);
  376.                         USART_Transmit(speed_real_RH+128);
  377.                         USART_Transmit(128);
  378.                 }
  379.                 else if(USART_State==0X33)                //ASCII码:0X30代表字符'3'
  380.                 {
  381.                         USART_Transmit(K_angle+128);
  382.                         USART_Transmit(K_angle_dot+128);
  383.                         USART_Transmit(K_position_dot+128);
  384.                 }                                
  385.                 temp=0;                                
  386.         }
  387.         speed_real_LH=0;
  388.         speed_real_RH=0;        
  389.         temp+=1;        
  390. }



  391. //-------------------------------------------------------
  392. int i,j;
  393. //-------------------------------------------------------
  394. void main(void)
  395. {
  396.         PORT_initial();
  397.         T2_initial();
  398.         INT_initial();
  399.         USART_initial ();
  400.         
  401.         SEI();
  402.         
  403.         K_position=0.8 * 0.209;                //换算系数:(256/10) * (2*pi/(64*12))=0.20944;//256/10:电压换算至PWM,256对应10V;
  404.         K_angle=34 * 25.6;                //换算系数:256/10 =25.6;
  405.         K_position_dot=1.09 * 20.9;                //换算系数:(256/10) * (25*2*pi/(64*12))=20.944;
  406.         K_angle_dot=2 * 25.6;                //换算系数:256/10 =25.6;
  407.         
  408.         for (i=1;i<=500;i++)                //延时启动PWM,等待卡尔曼滤波器稳定
  409.         {
  410.                 for (j=1;j<=300;j++);;
  411.                
  412.         }
  413.         T1_initial();
  414.                
  415.         while(1)
  416.         {
  417.                 ;
  418.         }
  419. }
复制代码

OK:其实控制很简单,利用定时器产生10ms中断,每次中断进行如下工作:
1:AD采样加速度传感器和陀螺仪,然后卡尔曼滤波得出角度与角速度;(滤波模块借鉴老外的)(互补滤波效果感觉不是很好);
2:计算车轮的速度,积分得出位置;
3:利用PD算法得出PWM值=K1*angle + K2*angle_dot + K3*speed + K4 * position;
4:前进后退给定参考速度计算即可;
注:
1:增量式轮速传感器A相中断,读取B相电位判断前进还是后退,在10ms的时间内累加,计算车轮的速度;
2:由于10ms累加的轮速信号不多,直接计算车体会发抖,所以增加了低通滤波,解决问题;
3:我在外面增加4个电位器可以手动调节4个K值,方便调试;
以上;

附件的文件夹里面有matlab建模的资料,可以求出所需的4个K的值,(参考NXT的),大家可自行调试K值看看对车辆的影响;

回复hahahagg
环路控制怎么这么慢啊?
振荡很久才稳定...
-----------------------------------------------------------------------

其实我也在思考,可能原因:
1:对速度做了滤波导致相位延迟;
2:倾斜角速度和车轮速度的K值不够大;
经实践:K值增大后,振荡减小;:)

回复format
敢问楼主的车轮和电机哪里买的?
价格?
-----------------------------------------------------------------------

电机在淘宝上买二手的,大几十元一个,因为考虑到减速器,而且带传感器;
其实后面实践得知可以不要位置传感器,这样选择电机的空间就很大了,买便宜的;
建议:
1:电机扭矩要大;说白了,功率要大;转动惯量大似乎的更稳定,所以不用空心杯的也许更好;
2:减速比可以大一点,用较大的轮子;
3:重心越靠下,抗扰性越强;但是可能会影响稳定性;

只是删除了RS232通讯和多余的注释(包括互补滤波),效果和原来的一样,这个是没有速度传感器的版本;片内8M振荡器: 无速度传感器的精简版程序.rar (2.34 KB, 下载次数: 19)

电路图: protel99格式的电路图_563431QMA9X5.rar (74.88 KB, 下载次数: 16)

全部制作资料: 平衡小车制作资料ProjectRII.rar (1.44 MB, 下载次数: 24)

</aden)|(1<<adsc)|(1<<adps1)|(1<</ursel)|(1<<ucsz0)|(1<</rxen)|(1<</int0)|(1<</isc01)|(1<<isc00)|(1<<isc11)|(1<</wgm21)|(1<<cs22)|(1<<cs21)|(1<</ocie2);
</wgm12)|(1<</com1a1)|(1<<wgm10)|(1<

评分

参与人数 1黑币 +12 收起 理由
YJGG + 12 很给力!

查看全部评分

分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏4 分享淘帖 顶 踩
回复

使用道具 举报

沙发
ID:113207 发表于 2016-4-10 02:58 | 只看该作者
常用计算方法有两种:互补滤波和卡尔曼滤波;都可以输出校正后的角度与角速度;我采用的是卡尔曼滤波,考虑到单片机的运算能力,是经过精简的,当然这些都是老外做的,关于卡尔曼滤波,我找了很多资料,较好的如下:
Kalman Filter Tutorial for Balancing Robot.pdf (76.47 KB, 下载次数: 15) kalman_intro_chinese.pdf (539.84 KB, 下载次数: 10)
互补滤波如下:点击此处下载 filter.rar (1.22 MB, 下载次数: 12)
我觉得卡尔曼的效果很好,但是初期有较大波动,要经过一段时间才能稳定;
http://www.docin.com/p-52169676.html
这个地方的可供参考;

回复

使用道具 举报

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

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

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