3. 开关电源建模基础
开关电源的环路补偿设计,本质上是在与一个动态系统打交道。如果只凭经验试凑,不仅效率低下,而且难以保证在全输入电压、全负载范围内的稳定性。本章的核心任务,是为Buck、Boost、Buck-Boost等基本拓扑建立统一的数学描述框架。我们将从状态空间平均法入手,推导出小信号模型,并最终得到控制到输出与输入到输出的传递函数。这套方法构成了数字补偿器设计的“被控对象”基础,也是后续所有频域分析和参数整定的前提。
3.1 状态空间平均法的核心原理
状态空间平均法由R.D. Middlebrook与Slobodan Ćuk于1976年提出,其基本思想是:在一个开关周期内,用占空比作为权重,对开关管导通与关断两个子拓扑的状态方程进行加权平均,从而得到一个连续时间、非线性但时不变的大信号模型。该模型忽略了开关频率附近的纹波分量,只保留直流与低频扰动分量,使得我们可以用经典控制理论进行分析。
以CCM模式下的Buck变换器为例,设电感电流iL和电容电压vC为状态变量,输入电压vg为输入变量,占空比d为控制变量。开关管导通期间(0 < t < dTs),状态方程为:
L·di_L/dt = v_g - v_C
C·dv_C/dt = i_L - v_C/R
开关管关断期间(dTs < t < Ts),续流二极管导通,状态方程为:
L·di_L/dt = -v_C
C·dv_C/dt = i_L - v_C/R
将两个子区间按占空比加权平均,得到大信号平均模型:
L·d⟨i_L⟩/dt = d·⟨v_g⟩ - ⟨v_C⟩
C·d⟨v_C⟩/dt = ⟨i_L⟩ - ⟨v_C⟩/R
这里⟨·⟩表示开关周期平均算子。值得注意的是,平均化操作将不连续的开关网络等效为一个连续系统,但保留了占空比对能量传递的控制作用。这一步骤的物理意义在于:电感电流和电容电压在开关频率处虽有纹波,但其低频包络才是我们关心的动态行为。
3.2 小信号模型推导:线性化的关键步骤
大信号平均模型是非线性的,因为存在d·⟨vg⟩这样的乘积项。为了应用线性控制理论,我们需要在静态工作点附近进行小信号线性化。设每个变量由直流分量与交流小信号分量组成:
⟨v_g⟩ = V_g + v̂_g
d = D + d̂
⟨i_L⟩ = I_L + î_L
⟨v_C⟩ = V_C + v̂_C
其中大写字母表示直流工作点,带“^”的小写字母表示小信号扰动。将上述表达式代入大信号平均模型,并忽略二阶小项(即扰动量的乘积项),同时利用直流工作点关系式进行化简,即可得到小信号线性状态方程。对于Buck变换器,直流工作点满足:
V_C = D·V_g
I_L = V_C/R
小信号方程整理为矩阵形式:
d/dt [î_L; v̂_C] = A·[î_L; v̂_C] + B·[v̂_g; d̂]
v̂_o = C·[î_L; v̂_C] + D·[v̂_g; d̂]
其中系统矩阵A、输入矩阵B、输出矩阵C和前馈矩阵D的具体表达式为:
A = [0, -1/L; 1/C, -1/(R·C)]
B = [D/L, V_g/L; 0, -I_L/C]
C = [0, 1]
D = [0, 0]
这一线性化过程的核心假设是扰动幅度远小于直流分量,即小信号条件。在实际电源系统中,输入电压扰动通常控制在额定值的10%以内,负载阶跃引起的输出电压变化也应在几个百分点量级,这保证了线性化模型的准确性。
3.3 传递函数的建立与物理意义
对小信号状态方程进行拉普拉斯变换,并求解输入到输出的传递函数,可以得到两个关键传递函数:控制到输出传递函数Gvd(s) = v̂_o(s)/d̂(s)|v̂_g=0,以及输入到输出传递函数Gvg(s) = v̂_o(s)/v̂_g(s)|d̂=0。对于理想Buck变换器(忽略电感寄生电阻和电容ESR),推导结果为:
G_vd(s) = V_g / (1 + s·L/R + s²·L·C)
G_vg(s) = D / (1 + s·L/R + s²·L·C)
从表达式可以看出,两个传递函数具有相同的二阶分母,其自然谐振频率为ω₀ = 1/√(LC),品质因数为Q = R·√(C/L)。这一对共轭极点是由LC滤波器决定的,是功率级固有的动态特性。控制到输出传递函数的直流增益为Vg,物理意义是占空比每变化1%,输出电压变化Vg·1%;输入到输出传递函数的直流增益为D,反映了输入电压扰动通过占空比直接传递到输出的比例。
在实际工程中,电容的ESR(等效串联电阻)会在高频引入一个零点,电感的寄生电阻会降低低频增益并增加阻尼。考虑这些寄生参数后,传递函数修正为:
G_vd(s) = V_g · (1 + s·R_ESR·C) / [1 + s·(L/R + R_L·C + R_ESR·C) + s²·L·C·(1 + R_L/R)]
其中RL为电感直流电阻,RESR为电容等效串联电阻。ESR零点通常位于fESR = 1/(2π·RESR·C),如果该零点频率低于穿越频率,会显著影响相位裕度,需要在补偿器设计中重点考虑。
| 传递函数 | 表达式(理想Buck) | 直流增益 | 主要极点/零点 | 设计关注点 |
|---|---|---|---|---|
| Gvd(s) | Vg / (1 + sL/R + s²LC) | Vg | 双极点 @ ω₀ = 1/√(LC) | 决定补偿器零点位置 |
| Gvg(s) | D / (1 + sL/R + s²LC) | D | 双极点 @ ω₀ | 输入电压前馈设计 |
| Zout(s) | sL / (1 + sL/R + s²LC) | 0 | 零点@0,双极点@ω₀ | 负载瞬态响应评估 |
3.4 实战案例:12V转3.3V Buck变换器建模
我们以一个典型的数字电源设计案例来完整走一遍建模流程。设计规格:输入电压Vin = 12V(允许波动±10%),输出电压Vout = 3.3V,最大负载电流Imax = 5A,开关频率fsw = 500kHz,电感L = 4.7μH(DCR = 15mΩ),输出电容C = 2×22μF(MLCC,ESR ≈ 3mΩ)。
第一步:计算直流工作点。稳态占空比D = Vout/Vin = 3.3/12 = 0.275。满载时电感电流平均值IL = Iout = 5A。电感电流纹波ΔIL = (Vin-Vout)·D/(L·fsw) = (12-3.3)×0.275/(4.7μ×500k) ≈ 1.02A,峰值为5.51A,确认CCM模式。
第二步:代入参数计算传递函数。考虑DCR和ESR后,控制到输出传递函数为:
G_vd(s) = 12·(1 + s·3m·44μ) / [1 + s·(4.7μ/1.32 + 15m·44μ + 3m·44μ) + s²·4.7μ·44μ·(1+15m/1.32)]
= 12·(1 + s·1.32e-7) / [1 + s·3.56e-6 + s²·2.07e-10]
由此得到:LC谐振频率f₀ = 1/(2π√(4.7μ×44μ)) ≈ 11.1kHz,ESR零点频率fESR = 1/(2π×3m×44μ) ≈ 1.2MHz。由于ESR零点远高于穿越频率(我们计划设计穿越频率为20kHz),该零点对环路影响可忽略,但需在仿真中确认。
第三步:验证模型准确性。使用SIMPLIS或LTspice搭建开关级仿真电路,在相同工作点下施加占空比阶跃扰动(如从0.275阶跃到0.280),记录输出电压响应。将仿真结果与传递函数的时域响应(通过MATLAB的step函数)对比,重点关注超调量和振荡频率。在本案例中,开关级仿真测得谐振频率约11.3kHz,与模型预测的11.1kHz偏差小于2%,验证了模型的准确性。
第四步:模型应用。基于Gvd(s)的幅频和相频特性,我们设计数字补偿器。在20kHz穿越频率处,功率级增益约为-12dB,相位滞后约-100°。因此补偿器需要在穿越频率处提供+12dB增益和约+50°相位提升,这决定了PID补偿器的零点与极点位置。后续章节将详细展开这一设计过程。
3.5 模型验证与工程注意事项
建立传递函数模型后,必须通过仿真或实测进行验证。推荐的验证流程为:首先在标称输入电压和半载条件下,测量控制到输出的频率响应(使用网络分析仪或扫频仪);然后对比实测Bode图与模型预测的Bode图,重点关注穿越频率附近的增益和相位;最后在不同输入电压(10.8V、12V、13.2V)和不同负载(10%、50%、100%)下重复测量,确认模型在全工况范围内的有效性。
对于数字控制实现,还需注意以下工程细节:
- 采样与保持效应:数字控制器在每个开关周期采样一次电感电流或输出电压,这引入了零阶保持器(ZOH)效应,等效于在控制环路中增加了一个延迟环节。在连续域建模时,通常用e-sTsw近似,或用Pade近似转换为有理函数。
- PWM调制器增益:数字PWM的增益为1/NPWM,其中NPWM为PWM计数器的峰值。例如,若PWM分辨率为12位(NPWM=4096),则调制器增益为1/4096 ≈ 0.000244,需在环路增益计算中计入。
- ADC量化效应:ADC分辨率限制了可检测的最小电压变化,量化噪声在低频段可视为白噪声,但在高频段可能影响环路稳定性。建议ADC分辨率不低于12位,且采样窗口应避开开关噪声尖峰。
最后需要强调的是,状态空间平均法得到的传递函数是设计补偿器的起点,而非终点。实际系统还受到元件公差、温度漂移、寄生参数分布等影响,因此设计完成后必须进行蒙特卡洛分析或最坏情况分析,确保在所有边界条件下环路均稳定且具有足够的相位裕度(通常要求≥45°)。