别再只用电量百分比了！用STM32的ADC和电流传感器，手把手教你实现更准的电池SOC估算

张

张建站

2026/5/29 4:55:40

10分钟阅读

别再只用电量百分比了！用STM32的ADC和电流传感器，手把手教你实现更准的电池SOC估算

嵌入式实战基于STM32的电池SOC精准估算系统设计在物联网设备和便携式电子产品中电池状态监测SOC的准确性直接影响用户体验和设备可靠性。许多开发者习惯使用简单的电压百分比法但这种方法在动态负载下误差可能高达20%-30%。本文将分享如何利用STM32内置ADC和外接电流传感器构建一个误差控制在5%以内的实用SOC估算系统。1. 硬件架构设计与关键元件选型1.1 电流检测方案对比选择适合的电流传感器是SOC估算的基础常见方案有传感器类型精度成本适用电流范围接口方式分流电阻运放±1%低10A差分ADCACS712霍尔效应±1.5%中30A单端ADCINA219数字芯片±0.5%较高3.2AI2C对于大多数嵌入式设备推荐使用50mΩ/1%精度的分流电阻配合AD8217差分放大器成本约$0.5典型电路如下// 电流检测电路校准参数 #define SHUNT_RESISTOR 0.05 // 50mΩ #define GAIN 20 // 放大器增益 #define ADC_REF 3.3 // ADC参考电压 float read_current(void) { uint16_t adc_raw HAL_ADC_GetValue(hadc1); float voltage (adc_raw / 4095.0) * ADC_REF; return (voltage / GAIN) / SHUNT_RESISTOR; // 单位安培 }1.2 电压采样优化技巧电池电压采样需注意使用1%精度的分压电阻在分压电路后添加0.1uF电容滤除高频噪声定期校准ADC基准电压void ADC_Calibration(void) { HAL_ADCEx_Calibration_Start(hadc1, ADC_SINGLE_ENDED); vrefint *(__IO uint16_t*)(0x1FFF75AA); // 读取内部参考电压 }2. 核心算法实现与优化2.1 库仑计实现关键细节电流积分法需要解决三个核心问题时间基准精度使用硬件定时器触发采样// 配置TIM2每100ms触发ADC采样 htim2.Instance TIM2; htim2.Init.Prescaler 1600-1; // 16MHz/1600 10kHz htim2.Init.CounterMode TIM_COUNTERMODE_UP; htim2.Init.Period 1000-1; // 10kHz/1000 10Hz HAL_TIM_Base_Start(htim2);累积误差处理采用滑动窗口平均滤波#define WINDOW_SIZE 10 float current_window[WINDOW_SIZE]; float coulomb_count 0.0; void update_soc(float current) { static uint8_t index 0; current_window[index] current; index (index 1) % WINDOW_SIZE; float avg_current 0; for(int i0; iWINDOW_SIZE; i) { avg_current current_window[i]; } avg_current / WINDOW_SIZE; coulomb_count avg_current * 0.1; // 100ms时间间隔 }容量温度补偿根据NTC读数调整额定容量float temp_compensation(float nominal_capacity, float temp_c) { // 锂离子电池典型温度补偿曲线 if(temp_c 0) return nominal_capacity * 0.7; if(temp_c 10) return nominal_capacity * 0.9; if(temp_c 45) return nominal_capacity * 0.95; return nominal_capacity; }2.2 电压-SOC曲线拟合实践不同电池化学体系的电压特性差异显著建议通过实验获取实际放电曲线。以18650锂离子电池为例// 基于实测数据的分段线性插值 float voltage_to_soc(float voltage) { const float curve[] { // 电压V, SOC% 3.00, 0, 3.50, 10, 3.70, 30, 3.85, 50, 4.00, 80, 4.20, 100 }; for(int i0; i5; i) { if(voltage curve[i*2] voltage curve[(i1)*2]) { float slope (curve[(i1)*21] - curve[i*21]) / (curve[(i1)*2] - curve[i*2]); return curve[i*21] slope * (voltage - curve[i*2]); } } return (voltage 3.0) ? 0 : 100; }3. 混合算法设计与实现3.1 动态权重调整策略结合库仑计和电压法的优势设计自适应混合算法typedef struct { float soc_coulomb; float soc_voltage; float weight; // 0-1, 偏向库仑计的程度 } soc_estimator; void update_estimator(soc_estimator *est, float current, float voltage) { // 更新库仑计部分 est-soc_coulomb (current / BAT_CAPACITY) * 0.1 * 3600; // 100ms间隔 // 更新电压法部分 est-soc_voltage voltage_to_soc(voltage); // 动态调整权重 if(fabs(current) 0.1) { // 小电流时信任电压法 est-weight 0.2; } else { // 大电流时信任库仑计 est-weight 0.8; } } float get_final_soc(soc_estimator *est) { return est-weight * est-soc_coulomb (1-est-weight) * est-soc_voltage; }3.2 充电状态特殊处理充电阶段需要特别处理检测充电器连接状态满电判断电流降至C/10以下充电效率补偿通常取95%-98%#define CHARGER_DETECT_PIN GPIO_PIN_4 bool is_charging(void) { return HAL_GPIO_ReadPin(GPIOA, CHARGER_DETECT_PIN) GPIO_PIN_SET; } void handle_charging(soc_estimator *est) { if(is_charging() current 0) { est-soc_coulomb * 0.97; // 充电效率补偿 } if(is_charging() current (BAT_CAPACITY/10)) { est-soc_coulomb 100.0; // 满电复位 } }4. 系统校准与维护机制4.1 工厂校准流程建议在生产环节执行三步校准零点校准void calibrate_zero_offset(void) { float sum 0; for(int i0; i100; i) { sum HAL_ADC_GetValue(hadc1); HAL_Delay(10); } zero_offset sum / 100; EEPROM_Write(ZERO_OFFSET_ADDR, zero_offset); }满量程校准# 使用标准电流源输入1A电流 $ apply_current 1.000 $ calibrate_full_scale容量标定void capacity_calibration(void) { // 完全放电后记录库仑计数 BAT_CAPACITY fabs(coulomb_count / 1000); // 转换为mAh EEPROM_Write(CAPACITY_ADDR, BAT_CAPACITY); }4.2 运行期自维护实现三种自修正机制OCV自学习当检测到静止状态电流0.05C持续5分钟时用当前电压修正SOCif(fabs(current) BAT_CAPACITY*0.05/1000) { idle_counter; if(idle_counter 3000) { // 5分钟 soc_coulomb voltage_to_soc(voltage); idle_counter 0; } } else { idle_counter 0; }容量衰减跟踪记录每次完整循环的放电量void update_capacity_degradation(void) { static float last_full_charge 0; if(soc_coulomb 95) last_full_charge coulomb_count; if(soc_coulomb 5) { float actual_capacity last_full_charge - coulomb_count; BAT_CAPACITY 0.9*BAT_CAPACITY 0.1*actual_capacity; } }温度补偿系数自适应根据历史数据优化温度曲线参数float learn_temp_coeff(float measured_soc, float temp_c) { static float coeff_table[6] { /* 初始值 */ }; int index (int)(temp_c / 10); coeff_table[index] 0.9*coeff_table[index] 0.1*measured_soc; return coeff_table[index]; }