如何更精准计算12.8V LiFePO4电池剩余电量百分比？

现有实现说明 #

Grafana数据 #

现有电量计算代码 #

static const uint8_t lifepo4[NUM_BATT_READINGS][2] = {
  {144, 100}, {136, 100}, {134,  99}, {133,  90},
  {132,  70}, {131,  40}, {130,  30}, {129,  20},
  {128,  17}, {125,  14}, {120,   9}, {100,   0}
};

int16_t measurements_get_batt_perc(int16_t mv)
{
  int16_t pct = 0;
  int32_t x1, x2, vD, pD;
  int32_t theta;
  if (mv >= lifepo4[0][0]*100)
      pct = lifepo4[0][1]; // cap value at 100%
  else
  {
    for (int i=0; i<NUM_BATT_READINGS-1; i++)
    {
      // put voltages into mV instead of 100s of mV
      x1 = (int32_t)(lifepo4[i][0])*100L;   // Higher value
      x2 = (int32_t)(lifepo4[i+1][0])*100L; // Lower value
      if (mv <= x1 && mv > x2)
      {
        vD = x1-x2;
        pD = (int32_t)(lifepo4[i][1] - lifepo4[i+1][1]);
        if (pD != 0)
          theta = pD*10000/vD;
        else
          theta = 0;
        pct = lifepo4[i][1] - (x1-mv)*theta/10000;
        break;
      }
    }
  }

  return pct;
}

参考电压曲线 #

上述代码基于LiFePO4电池电压-电量查找表计算剩余电量，对应的电压曲线如下：

问题 #

现有方法未考虑负载变化对电压的影响，导致电量估算不准确。请问是否有更优的计算方案？是否存在适用于变负载场景的电量计算公式？或如何将电流积分与电压数据进行融合以提升精度？

优化方案 #

1. 引入开路电压（OCV）校准 #

负载变化会导致电池端电压出现压降（放电）或抬升（充电），直接用端电压估算误差极大，核心解决思路是获取开路电压：

• 静置法：当电池电流接近0且持续30分钟以上时，读取端电压作为OCV，此时用原有LUT查表即可得到准确电量；
• 动态估算：无法静置时，用Thevenin等效电路模型计算OCV，公式示例：
  OCV = V + I*(R0 + R1*(1 - e^(-t/τ1)))
  其中V是端电压，I是负载电流（放电为正、充电为负），R0为欧姆内阻，R1为极化内阻，τ1为极化时间常数，这些参数需通过电池实测获取。

2. 电流积分（库仑计数）+ 电压校准 #

库仑计数通过积分电流计算电量消耗/补充，精度高但存在累计误差，需结合电压定期校准：

• 基础计算：剩余容量 = 初始容量 - ∫(I*dt)（放电时I为正，充电时为负）
• 校准机制：
  • 静置状态下，用OCV对应的电量修正库仑计数的累计值；
  • 充电至满电时，重置剩余容量为额定容量，清零累计误差；
  • 每日一次结合OCV与LUT结果，微调库仑计数的剩余容量。

3. 融合算法简化实现示例 #

// 全局变量
float remaining_capacity = BATTERY_RATED_CAPACITY; // 初始额定容量（Ah）
float last_current = 0;
uint32_t last_time_ms = 0;

int16_t get_batt_perc(int16_t mv, int16_t ma, uint32_t current_time_ms) {
    // 1. 更新库仑计数
    if (last_time_ms != 0) {
        float delta_time_h = (current_time_ms - last_time_ms) / 3600000.0f;
        float delta_cap = ma * delta_time_h / 1000.0f; // 毫安转安时
        remaining_capacity -= delta_cap;
        // 限制容量范围，防止溢出
        remaining_capacity = max(0.0f, min(remaining_capacity, BATTERY_RATED_CAPACITY));
    }
    last_current = ma;
    last_time_ms = current_time_ms;

    // 2. 估算OCV并校准
    int16_t ocv_batt_perc;
    if (abs(ma) < 10) { // 电流极小，视为静置，直接用端电压当OCV查表
        ocv_batt_perc = measurements_get_batt_perc(mv);
        // 校准库仑计数
        remaining_capacity = (ocv_batt_perc / 100.0f) * BATTERY_RATED_CAPACITY;
    } else {
        // 动态负载下用库仑计数结果，每日校准一次
        static uint32_t last_calib_time = 0;
        if (current_time_ms - last_calib_time > 86400000) {
            float estimated_ocv = estimate_ocv(mv, ma); // 用Thevenin模型估算OCV
            ocv_batt_perc = lookup_ocv_perc((int16_t)(estimated_ocv * 1000)); // 转mV查表
            float calib_factor = (ocv_batt_perc / 100.0f) / (remaining_capacity / BATTERY_RATED_CAPACITY);
            remaining_capacity *= calib_factor;
            last_calib_time = current_time_ms;
        }
    }

    // 计算最终百分比
    int16_t final_perc = (remaining_capacity / BATTERY_RATED_CAPACITY) * 100;
    return clamp(final_perc, 0, 100);
}

4. 额外优化点 #

• 温度补偿：LiFePO4的OCV随温度变化，需在LUT中加入温度维度，或对OCV做温度修正；
• 老化补偿：定期通过满充放循环更新电池额定容量，适配电池衰减；
• 内阻监测：长期监测内阻变化，辅助判断电池健康状态，修正估算模型。

内容的提问来源于stack exchange，提问作者Jdelmore