利用STM32实现自平衡机器人功能项目
将机器人整体开源，同时总结一下机器人搭建过程中遇到的坑和未来的改进方向。在分享的文件里包含了结构设计、程序控制、电路设计以及其他模块相关资料供大家参考。
第一：机器人原理分析
首先来看成品图：
如图所示，该机器人根据陀螺仪的位姿数据，通过三个全向轮驱动底部球体调整自己在球上的位置，保持动态平衡的同时实现全向移动。
保持动态平衡过程需要对机器人进行运动学分析，这里参考了平衡小车之家的运动学方程：
自平衡控制问题转化为三步：输入x、y角度—控制器计算—输出a、b、c电机转速的控制模型。
第二：控制器设计
首先考虑参考平衡车控制，球上自平衡机器人本质上依然是一个一阶倒立摆问题。
这里参考了飞思卡尔直立车的控制方法，采用串级pid控制器，外环pd角度环，内环速度pi环。
由于我的驱动方案选择的是42步进电机，在速度闭环的时候有些问题。正常的直流电机+编码器的控制方案可以通过编码器将轮子的真实速度计算出来，从而和控制器的理想转速作差，实现速度控制。
而我这里的速度闭环是通过计算上一个时钟周期时给步进电机的控制量，通过运动学方程分解，得到机器人的虚拟速度，与理想转速作差控制。我认为这种速度闭环方式还是存在一定缺陷的，但是在网上查看论文的时候我发现有很多自平衡机器人都是用42步进电机来实现速度闭环的，不知道是什么方法。
这里还可以好好思考一下为什么角度环要用pd控制，速度环要pi控制，角度环的p部分和d部分对机器人控制有什么影响？在很多csdn调试平衡车的博客中都有解释，这里就留给大家思考了。
第三：硬件及结构设计
自平衡机器人的硬件清单有：
56mm全向轮 45元/个
42步进电机 25/个
42步进闭环模块 59.8元/个
lm2596s降压模块 20元
stm32f103c8t6-4飞控板 59.8元
gy-521六轴陀螺仪 25元
用到的模块大致如上所示，c8t6的价格随着最近芯片涨价直线上升，我白嫖了实验室的两块板子，现在买一块实在太贵，可以等芯片价格稳定一些再买。其余开关排针等常见元件不再赘述。
电路原理图如下所示：
机器人使用solidworks设计整体结构，底板可在某宝定制6050太空铝切割，蓝色件为正常3d打印件。
第四：程序部分
在keil 5中开发stm32。
控制程序采用定时器0.5ms定时中断的方式进行计算，每触发两次中断计算对电机控制一次，这里还是推荐大家采用外部中断读取gy-521上的int引脚的方式，控制计算周期。gy-521上的int引脚每5ms触发一次跳变，采用外部中断的方式可以严格保证读取位姿数据与计算处理同步。
int tim1_up_irqhandler (void) { u8 key_cal; if(tim_getitstatus(tim1,tim_it_update) != reset) { tim_clearitpendingbit(tim1,tim_it_update); flag_target=!flag_target; key_cal=key_scan(0); if(state_flag==1)//矫正结束 { if(flag_target==1)//每读取两次陀螺仪控制一次 { read_dmp(); //===读取倾角 scope(); return 0; } } if(key_cal==1)//矫正按键 { angle_zero_x=angle_balance_x; angle_zero_y=angle_balance_y; key_cal=0; flag_stop=0; } if(key_cal==2||key_cal==3)//矫正按键 { flag_stop=1;//关闭速度环i积分 key_cal=0; } angle_bias_x =angle_balance_x-angle_zero_x; //获取y方向的偏差 angle_bias_y =angle_balance_y-angle_zero_y; //获取y方向的偏差 if(control_mode==0)//pid控制模式 { encoder_analysis(motor_a,motor_b,motor_c); //正运动学分析，得到x y方向的速度 balance_pwm_x= balance_x(angle_bias_x,gyro_balance_x);//x方向的倾角控制 balance_pwm_y=-balance_y(angle_bias_y,gyro_balance_y); //y方向的倾角控制// if(++flag_target_2==4)//速度环频率慢于加速度环 但是还没加速度环 // { velocity_pwm_x=velocity_x(compute_x); //x方向的速度控制 velocity_pwm_y=velocity_y(compute_y); //y方向的速度控制 // flag_target_2=0;// } move_x =balance_pwm_x+velocity_pwm_x; //===x方向控制量累加 move_y =balance_pwm_y+velocity_pwm_y; //===y方向控制量累加 move_z=0; kinematic_analysis(move_x,move_y,move_z);//逆运动学分析得到abc电机控制量 } motor_a=target_a;//直接调节pwm频率 motor_b=target_b;//直接调节pwm motor_c=target_c;//直接调节pwm//以下都是为了速度连续化处理防止突变 if(motor_a==0) motor_a=motor_a_last; if(motor_b==0) motor_b=motor_b_last; if(motor_c==0) motor_c=motor_c_last; xianfu_pwm(2000); set_pwm(motor_a,motor_b,motor_c); gyro_balance_x_last=gyro_balance_x; gyro_balance_y_last=gyro_balance_y; gyro_balance_z_last=gyro_balance_z; angle_balance_x_last=angle_balance_x; angle_balance_y_last=angle_balance_y; angle_balance_z_last=angle_balance_z; motor_a_last=motor_a; motor_b_last=motor_b; motor_c_last=motor_c; } return 0;} 对于电机控制，由于采用的驱动方案是步进电机，调速的方式是改变驱动步进电机的脉冲频率。我这里选择了三个定时器，动态调节定时器的频率，具体方式是在初始化时设定好定时器的预分频系数psc的值，然后在程序里动态更改arr寄存器的值，从而改变定时器的定时频率。//这里以a电机的速度控制为例 输入为 电机方向和电机速度void set_motora_speed(u8 dir,u16 speed){ u32 arr; arr=speed; tim_arrpreloadconfig(tim3,disable); tim3->arr=arr;//计数到10000在归零重新计数 tim3->ccr4=arr/2;//保持占空比为50% tim_arrpreloadconfig(tim3,enable); tim_cmd(tim3,enable); if(dir==0) { gpio_setbits(gpioa,gpio_pin_1); } else { gpio_resetbits(gpioa,gpio_pin_1); }}
小车的运动学分解代码实现如下，参考了平衡小车之家的代码：
/**********************************************************函数功能：小车运动数学模型入口参数：x y z 三轴速度或者位置返回 值：无***********************************************************/void kinematic_analysis(float vx,float vy,float vz){ target_a = vx + l_parameter*vz; target_b = -x_parameter*vx + y_parameter*vy + l_parameter*vz; target_c = -x_parameter*vx - y_parameter*vy + l_parameter*vz;}/*****************************************************************函数功能：小车运动 正运动学分析 入口参数：a b c三个电机的速度返回 值：无******************************************************************/void encoder_analysis(float va,float vb,float vc){ compute_x=(va*2-vb-vc); compute_y=((vb-vc)*sqrt(3)); compute_z=(va+vb+vc); } 其余代码不全放出，可在文末点击“阅读原文”下载查看。
第五：总结与展望
球上自平衡机器人可以作为算法试验平台， 输入输出固定，更换不同控制器，将数据导入matlab进行分析即可比较控制器性能。
个人认为结构有两个改进方向，一方面参考以下论文：余义. 单球驱动自平衡机器人位姿解算与控制系统研究[d].武汉科技大学,2019。论文中采用的四足式驱动结构更有利于机器人自平衡控制。
另一方面可以增加球体和机器人固定装置，利用机械结构将机器人与底部驱动球结合成一个整体防止机器人跳轮等问题。同时驱动球对于机器人平衡的影响较大，最好还是定制空心钢球，然后喷漆增大摩擦力，最有利于机器人自平衡控制。
控制部分的改进，首先是控制原理，本文是针对建立好的运动学方程进行分析，通过串级pid算法来实现自平衡运动。该机器人的控制问题本质上是一阶倒立摆问题，可以采用动力学建模的方式，通过动力学分析算出平衡需要的虚拟力矩，再对电机进行力矩控制。
其次是控制器，pid控制算法应用广泛但也有一定的缺点，可以考虑采用模糊pid，adrc自抗扰控制器，强化学习等智能控制算法对机器人自平衡进行控制。
Energous 和 橙群微电子推出面向物联网工业应用的低维护无线传感器
pcb电路设计的实用小技巧
VR全景虚拟现实线上旅游的优势都有哪些
AI空前的发展热度，推动AI城市名片争夺战进入了白热化阶段
华为P10什么时候出 华为P10和华为P9的区别在哪里
利用STM32实现自平衡机器人功能项目
led点阵8×8滚动显示原理_led点阵屏编程
传音攻守并重，海外新兴市场迎可持续增长新机遇
如何判断光纤收发器是否有问题
学习Raspberry Pi树莓派不可忽略的十个要点
本源量子携手瀚海博兴实现新冠病毒检测领域新突破
电磁兼容暗室是什么
频频缺货,为何MLCC上了天、铝电解电容却涨不起来？
诺基亚3310复刻版最新消息：诺基亚3310复刻版正式开启预约，外观、配置升级游戏回归经典
低压无功补偿装置MSCGD的原理结构和作用概述
进击的苏宁 苏宁欲利用区块链助力智慧零售健康发展
毫米科技的AI布局思路以及智能家居时代变革的经验分享
酷芯携通信SoC AR8030亮相滴水湖中国RISC-V产业论坛
什么是区块链以及加密货币
光载毫米波无线电通信技术的现状与发展
第一：机器人原理分析
首先来看成品图：
如图所示，该机器人根据陀螺仪的位姿数据，通过三个全向轮驱动底部球体调整自己在球上的位置，保持动态平衡的同时实现全向移动。
保持动态平衡过程需要对机器人进行运动学分析，这里参考了平衡小车之家的运动学方程：
自平衡控制问题转化为三步：输入x、y角度—控制器计算—输出a、b、c电机转速的控制模型。
第二：控制器设计
首先考虑参考平衡车控制，球上自平衡机器人本质上依然是一个一阶倒立摆问题。
这里参考了飞思卡尔直立车的控制方法，采用串级pid控制器，外环pd角度环，内环速度pi环。
由于我的驱动方案选择的是42步进电机，在速度闭环的时候有些问题。正常的直流电机+编码器的控制方案可以通过编码器将轮子的真实速度计算出来，从而和控制器的理想转速作差，实现速度控制。
而我这里的速度闭环是通过计算上一个时钟周期时给步进电机的控制量，通过运动学方程分解，得到机器人的虚拟速度，与理想转速作差控制。我认为这种速度闭环方式还是存在一定缺陷的，但是在网上查看论文的时候我发现有很多自平衡机器人都是用42步进电机来实现速度闭环的，不知道是什么方法。
这里还可以好好思考一下为什么角度环要用pd控制，速度环要pi控制，角度环的p部分和d部分对机器人控制有什么影响？在很多csdn调试平衡车的博客中都有解释，这里就留给大家思考了。
第三：硬件及结构设计
自平衡机器人的硬件清单有：
56mm全向轮 45元/个
42步进电机 25/个
42步进闭环模块 59.8元/个
lm2596s降压模块 20元
stm32f103c8t6-4飞控板 59.8元
gy-521六轴陀螺仪 25元
用到的模块大致如上所示，c8t6的价格随着最近芯片涨价直线上升，我白嫖了实验室的两块板子，现在买一块实在太贵，可以等芯片价格稳定一些再买。其余开关排针等常见元件不再赘述。
电路原理图如下所示：
机器人使用solidworks设计整体结构，底板可在某宝定制6050太空铝切割，蓝色件为正常3d打印件。
第四：程序部分
在keil 5中开发stm32。
控制程序采用定时器0.5ms定时中断的方式进行计算，每触发两次中断计算对电机控制一次，这里还是推荐大家采用外部中断读取gy-521上的int引脚的方式，控制计算周期。gy-521上的int引脚每5ms触发一次跳变，采用外部中断的方式可以严格保证读取位姿数据与计算处理同步。
int tim1_up_irqhandler (void) { u8 key_cal; if(tim_getitstatus(tim1,tim_it_update) != reset) { tim_clearitpendingbit(tim1,tim_it_update); flag_target=!flag_target; key_cal=key_scan(0); if(state_flag==1)//矫正结束 { if(flag_target==1)//每读取两次陀螺仪控制一次 { read_dmp(); //===读取倾角 scope(); return 0; } } if(key_cal==1)//矫正按键 { angle_zero_x=angle_balance_x; angle_zero_y=angle_balance_y; key_cal=0; flag_stop=0; } if(key_cal==2||key_cal==3)//矫正按键 { flag_stop=1;//关闭速度环i积分 key_cal=0; } angle_bias_x =angle_balance_x-angle_zero_x; //获取y方向的偏差 angle_bias_y =angle_balance_y-angle_zero_y; //获取y方向的偏差 if(control_mode==0)//pid控制模式 { encoder_analysis(motor_a,motor_b,motor_c); //正运动学分析，得到x y方向的速度 balance_pwm_x= balance_x(angle_bias_x,gyro_balance_x);//x方向的倾角控制 balance_pwm_y=-balance_y(angle_bias_y,gyro_balance_y); //y方向的倾角控制// if(++flag_target_2==4)//速度环频率慢于加速度环 但是还没加速度环 // { velocity_pwm_x=velocity_x(compute_x); //x方向的速度控制 velocity_pwm_y=velocity_y(compute_y); //y方向的速度控制 // flag_target_2=0;// } move_x =balance_pwm_x+velocity_pwm_x; //===x方向控制量累加 move_y =balance_pwm_y+velocity_pwm_y; //===y方向控制量累加 move_z=0; kinematic_analysis(move_x,move_y,move_z);//逆运动学分析得到abc电机控制量 } motor_a=target_a;//直接调节pwm频率 motor_b=target_b;//直接调节pwm motor_c=target_c;//直接调节pwm//以下都是为了速度连续化处理防止突变 if(motor_a==0) motor_a=motor_a_last; if(motor_b==0) motor_b=motor_b_last; if(motor_c==0) motor_c=motor_c_last; xianfu_pwm(2000); set_pwm(motor_a,motor_b,motor_c); gyro_balance_x_last=gyro_balance_x; gyro_balance_y_last=gyro_balance_y; gyro_balance_z_last=gyro_balance_z; angle_balance_x_last=angle_balance_x; angle_balance_y_last=angle_balance_y; angle_balance_z_last=angle_balance_z; motor_a_last=motor_a; motor_b_last=motor_b; motor_c_last=motor_c; } return 0;} 对于电机控制，由于采用的驱动方案是步进电机，调速的方式是改变驱动步进电机的脉冲频率。我这里选择了三个定时器，动态调节定时器的频率，具体方式是在初始化时设定好定时器的预分频系数psc的值，然后在程序里动态更改arr寄存器的值，从而改变定时器的定时频率。//这里以a电机的速度控制为例 输入为 电机方向和电机速度void set_motora_speed(u8 dir,u16 speed){ u32 arr; arr=speed; tim_arrpreloadconfig(tim3,disable); tim3->arr=arr;//计数到10000在归零重新计数 tim3->ccr4=arr/2;//保持占空比为50% tim_arrpreloadconfig(tim3,enable); tim_cmd(tim3,enable); if(dir==0) { gpio_setbits(gpioa,gpio_pin_1); } else { gpio_resetbits(gpioa,gpio_pin_1); }}
小车的运动学分解代码实现如下，参考了平衡小车之家的代码：
/**********************************************************函数功能：小车运动数学模型入口参数：x y z 三轴速度或者位置返回 值：无***********************************************************/void kinematic_analysis(float vx,float vy,float vz){ target_a = vx + l_parameter*vz; target_b = -x_parameter*vx + y_parameter*vy + l_parameter*vz; target_c = -x_parameter*vx - y_parameter*vy + l_parameter*vz;}/*****************************************************************函数功能：小车运动 正运动学分析 入口参数：a b c三个电机的速度返回 值：无******************************************************************/void encoder_analysis(float va,float vb,float vc){ compute_x=(va*2-vb-vc); compute_y=((vb-vc)*sqrt(3)); compute_z=(va+vb+vc); } 其余代码不全放出，可在文末点击“阅读原文”下载查看。
第五：总结与展望
球上自平衡机器人可以作为算法试验平台， 输入输出固定，更换不同控制器，将数据导入matlab进行分析即可比较控制器性能。
个人认为结构有两个改进方向，一方面参考以下论文：余义. 单球驱动自平衡机器人位姿解算与控制系统研究[d].武汉科技大学,2019。论文中采用的四足式驱动结构更有利于机器人自平衡控制。
另一方面可以增加球体和机器人固定装置，利用机械结构将机器人与底部驱动球结合成一个整体防止机器人跳轮等问题。同时驱动球对于机器人平衡的影响较大，最好还是定制空心钢球，然后喷漆增大摩擦力，最有利于机器人自平衡控制。
控制部分的改进，首先是控制原理，本文是针对建立好的运动学方程进行分析，通过串级pid算法来实现自平衡运动。该机器人的控制问题本质上是一阶倒立摆问题，可以采用动力学建模的方式，通过动力学分析算出平衡需要的虚拟力矩，再对电机进行力矩控制。
其次是控制器，pid控制算法应用广泛但也有一定的缺点，可以考虑采用模糊pid，adrc自抗扰控制器，强化学习等智能控制算法对机器人自平衡进行控制。
Energous 和 橙群微电子推出面向物联网工业应用的低维护无线传感器
pcb电路设计的实用小技巧
VR全景虚拟现实线上旅游的优势都有哪些
AI空前的发展热度，推动AI城市名片争夺战进入了白热化阶段
华为P10什么时候出 华为P10和华为P9的区别在哪里
利用STM32实现自平衡机器人功能项目
led点阵8×8滚动显示原理_led点阵屏编程
传音攻守并重，海外新兴市场迎可持续增长新机遇
如何判断光纤收发器是否有问题
学习Raspberry Pi树莓派不可忽略的十个要点
本源量子携手瀚海博兴实现新冠病毒检测领域新突破
电磁兼容暗室是什么
频频缺货,为何MLCC上了天、铝电解电容却涨不起来？
诺基亚3310复刻版最新消息：诺基亚3310复刻版正式开启预约，外观、配置升级游戏回归经典
低压无功补偿装置MSCGD的原理结构和作用概述
进击的苏宁 苏宁欲利用区块链助力智慧零售健康发展
毫米科技的AI布局思路以及智能家居时代变革的经验分享
酷芯携通信SoC AR8030亮相滴水湖中国RISC-V产业论坛
什么是区块链以及加密货币
光载毫米波无线电通信技术的现状与发展