STM32 PID温度控制实战：参数自整定与PWM实现

1. STM32 PID温度控制实战：从参数整定到PWM输出

做嵌入式开发这些年，PID控制算法绝对是让人又爱又恨的存在。特别是做温度控制这类大惯性系统时，手动调参简直就像在迷宫里转悠。最近用STM32做了个温控项目，把自整定算法和抗积分饱和这些坑都踩了一遍，今天就把完整实现方案和调试心得整理出来。

先说说这个系统的核心需求：用STM32驱动加热片，通过PWM控制功率，将温度稳定在设定值（比如100℃）。难点在于加热系统有较大热惯性，传统手动PID调参很难达到理想效果。我的解决方案是用阶跃响应法实现PID参数自整定，配合抗积分饱和处理，最终实现了±2℃的控制精度。

2. PID控制核心实现

2.1 数据结构设计

先看PID算法的数据结构设计，这是整个系统的核心：

c复制typedef struct {
    float Kp;         // 比例系数
    float Ki;         // 积分系数  
    float Kd;         // 微分系数
    float integral_max; // 积分限幅
    float output_max;   // 输出限幅
    float last_error;   // 上次误差
    float integral;     // 积分累计
} PID_TypeDef;

这个结构体有几个关键设计点：

1. integral_max：这是解决积分饱和的关键。当系统长时间存在误差时，积分项会不断累积导致输出过大，这个参数限制了积分项的最大值
2. output_max：物理系统都有输出限制（如PWM最大100%占空比），这个参数确保输出不会超出硬件限制
3. 分离了last_error和integral：方便在不同中断周期中保持状态

实际调试中发现，对于温度控制系统，integral_max设为output_max的20%-30%效果最佳。太大容易超调，太小则稳态误差难以消除

2.2 中断服务实现

温度控制需要定时采样，我用TIM2定时器触发中断，基准周期设为100ms：

c复制void TIM2_IRQHandler(void) {
    static uint32_t tick = 0;
    if(__HAL_TIM_GET_FLAG(&htim2, TIM_FLAG_UPDATE)) {
        __HAL_TIM_CLEAR_FLAG(&htim2, TIM_FLAG_UPDATE);
        
        float current_temp = DS18B20_Read(); // 读取温度传感器
        float error = target_temp - current_temp;
        
        // 积分项计算与限幅
        pid.integral += error * dt;
        if(pid.integral > pid.integral_max) 
            pid.integral = pid.integral_max;
        else if(pid.integral < -pid.integral_max) 
            pid.integral = -pid.integral_max;
        
        // 微分项计算
        float derivative = (error - pid.last_error) / dt;
        
        // PID输出计算
        float output = pid.Kp * error + pid.Ki * pid.integral + pid.Kd * derivative;
        
        pid.last_error = error;
        __HAL_TIM_SET_COMPARE(&htim3, TIM_CHANNEL_1, (uint32_t)(output));
    }
}

几个关键细节：

1. dt取值：温度系统惯性大，采样周期建议50-200ms。周期太短会导致微分项噪声大，太长则响应迟钝
2. 微分处理：实际使用中可以加一阶低通滤波，避免传感器噪声被放大
3. 输出限幅：虽然结构体中有output_max，但在中断中还要再做一次限制，确保安全

3. 自整定算法实现

3.1 阶跃响应法原理

手动调参费时费力，我采用Ziegler-Nichols阶跃响应法实现自整定。基本原理是：

1. 逐步增大比例增益，使系统产生等幅振荡
2. 记录临界增益Ku和振荡周期Tu
3. 根据公式计算PID参数：
   • Kp = 0.6 * Ku
   • Ki = 1.2 * Ku / Tu
   • Kd = 0.075 * Ku * Tu

3.2 代码实现

c复制void PID_AutoTune(void) {
    float Ku, Tu;
    float output = 0;
    uint8_t oscillation_count = 0;
    
    // 禁用积分项
    float temp_integral = pid.integral;
    pid.integral = 0;
    
    while(1) {
        output += 0.5f; // 每次增加0.5%输出
        __HAL_TIM_SET_COMPARE(&htim3, TIM_CHANNEL_1, (uint32_t)output);
        HAL_Delay(100); // 等待系统响应
        
        float peak1 = 0, peak2 = 0;
        if(DetectOscillation(&peak1, &peak2)) {
            Ku = output / (peak1 - peak2); // 临界增益
            Tu = CalculateOscillationPeriod(); // 振荡周期
            break;
        }
    }
    
    // Ziegler-Nichols参数计算
    pid.Kp = 0.6 * Ku;
    pid.Ki = 1.2 * Ku / Tu;
    pid.Kd = 0.075 * Ku * Tu;
    
    // 恢复积分项
    pid.integral = temp_integral;
}

振荡检测函数需要记录温度极值点：

c复制bool DetectOscillation(float *peak1, float *peak2) {
    static float last_temp = 0;
    static bool rising = false;
    static uint8_t peak_count = 0;
    
    float current_temp = DS18B20_Read();
    
    if(current_temp > last_temp) {
        rising = true;
    } else if(current_temp < last_temp && rising) {
        // 检测到峰值
        if(peak_count == 0) {
            *peak1 = last_temp;
            peak_count++;
        } else {
            *peak2 = last_temp;
            return true;
        }
        rising = false;
    }
    
    last_temp = current_temp;
    return false;
}

实测发现，对于典型温控系统，振荡周期Tu通常在30-60秒之间。自整定时要有足够耐心，整个过程可能需要5-10分钟

4. PWM输出配置

4.1 定时器参数设计

加热控制通过PWM实现，使用TIM3通道1输出：

c复制htim3.Instance = TIM3;
htim3.Init.Prescaler = 84-1;  // 84MHz/84 = 1MHz
htim3.Init.CounterMode = TIM_COUNTERMODE_UP; 
htim3.Init.Period = 1000-1;   // 1MHz/1000 = 1kHz PWM
htim3.Init.ClockDivision = TIM_CLOCKDIVISION_DIV1;
HAL_TIM_PWM_Start(&htim3, TIM_CHANNEL_1);

频率选择考虑：

1. 1kHz频率足够高，能避免可闻噪声
2. 对于电阻丝加热，这个频率也远高于热惯性时间常数
3. 如果驱动继电器，必须降到10Hz以下，否则会快速损坏触点

4.2 动态输出限制

系统在不同阶段需要不同的控制策略：

c复制// 根据误差动态调整输出限制
if(fabs(error) > 20.0f) {
    pid.output_max = 100.0f; // 大温差时全功率加热
} else if(fabs(error) > 5.0f) {
    pid.output_max = 70.0f;  // 中等温差时限制输出
} else {
    pid.output_max = 50.0f;  // 小温差时精细调节
}

这种变限幅策略能显著加快系统响应，实测从25℃升温到100℃的时间缩短了约40%

5. 调试经验与问题排查

5.1 常见问题处理

5.2 传感器处理技巧

1. 软件滤波：DS18B20读数可以取滑动平均

c复制#define FILTER_LEN 5
float temp_history[FILTER_LEN];
float filtered_temp = 0;

for(int i=0; i<FILTER_LEN-1; i++){
    temp_history[i] = temp_history[i+1];
}
temp_history[FILTER_LEN-1] = DS18B20_Read();

for(int i=0; i<FILTER_LEN; i++){
    filtered_temp += temp_history[i];
}
filtered_temp /= FILTER_LEN;

2. 布线注意：DS18B20信号线要远离PWM线，避免干扰

5.3 进阶优化方向

1. 模糊PID：根据误差大小动态调整PID参数
2. 前馈控制：加入环境温度补偿
3. 多段温控：实现升温-保温多段曲线

这套系统在3D打印机热床控制上实测效果很好，从室温加热到100℃约3分钟，超调不超过2℃，稳态误差±0.5℃。最关键的是自整定功能让参数调试变得非常简单，更换不同加热元件后只需重新自整定即可获得理想参数。