LLC谐振变换器数字补偿器设计：从状态空间建模到2P2Z实现

1. 项目概述：为什么数字补偿是LLC谐振变换器的“大脑升级”？

在电源设计领域，LLC谐振变换器因其高效率、高功率密度和软开关特性，早已成为中高功率应用的首选拓扑。但很多工程师朋友在实际调试时，常常会遇到一个头疼的问题：环路稳定性。传统的模拟补偿器，比如Type II或Type III，虽然经典，但面对LLC这种高阶、非线性、工作点随负载和输入电压剧烈变化的系统，往往显得力不从心。调好了满载，轻载可能振荡；输入电压一变，环路特性又漂了。这背后的核心原因，是模拟补偿器的零极点位置是固定的，无法自适应地跟随被控对象的变化。

这就是我们这次要深入探讨的“LLC谐振变换器数字补偿器设计”。简单来说，就是给LLC变换器装上一个“数字大脑”——用一颗微控制器（MCU）或数字信号处理器（DSP），实时采样输出电压和电流，通过运行在芯片内部的数字控制算法（补偿器），计算出精确的开关频率或占空比指令，从而实现对输出电压的精准、快速且稳定的控制。相较于模拟补偿，数字补偿的优势是颠覆性的：参数可在线调整、易于实现复杂控制算法（如自适应、前馈）、无温漂和老化问题，并且便于监控和通信。

本次设计的核心路径是“从状态空间建模到2P2Z控制器实现”。这听起来有点学术，但拆解开来就非常清晰：状态空间建模是为了让我们在数字域里，用数学方程精确地描述LLC变换器这个物理系统的动态行为，这是设计控制器的前提。而2P2Z（两极点两零点）控制器，则是一种在数字控制中非常经典且实用的补偿器结构，它足够灵活，能够为LLC这类系统提供所需的相位裕度和增益裕度。整个项目，就是打通从理论模型到实际代码的完整链路，让你不仅能调通环路，更能深刻理解每一个参数背后的物理意义和设计逻辑。

2. 核心思路与方案选型：为何是状态空间模型+2P2Z？

面对LLC的数字控制，市面上有各种方案，比如直接移植模拟补偿器的数字化（如双线性变换），或者用数字PID。我们选择“状态空间平均建模 + 2PZ数字补偿器”这条路径，是经过深思熟虑的，主要基于以下几点考量：

2.1 为什么选择状态空间平均模型？

LLC谐振变换器是一个典型的非线性、时变系统。直接对其非线性方程进行控制器设计极其困难。状态空间平均法，通过在开关周期内对状态变量（通常是电感电流和电容电压）进行平均，巧妙地将其线性化，得到一个能在开关频率以下描述系统低频小信号特性的线性时不变模型。这个模型是设计线性控制器（如2P2Z）的黄金标准。

• 精准性：相比于简单的等效电路法，状态空间平均模型能更准确地反映谐振腔（Lr, Cr, Lm）的动态交互，尤其是增益曲线在谐振点附近的特性，这对环路设计至关重要。
• 通用性：模型以状态矩阵（A, B, C, D）的形式呈现，与负载类型（电阻、电流源等）和变压器模型（考虑磁化电感）可以很好地结合，模型扩展性强。
• 便于数字化：得到的连续域传递函数，可以直接应用成熟的离散化方法（如零阶保持器法）转换为离散域模型，为数字控制器设计铺平道路。

2.2 为什么选择2P2Z结构？

数字补偿器的形式有很多，从最简单的PI到更复杂的PID、PR控制器等。2P2Z结构在这里脱颖而出：

• 结构匹配：LLC电压模式控制（通过调节开关频率来稳压）的开环传递函数，通常呈现双极点特性（由输出LC滤波器主导），并带有由右半平面零点（如果有）或谐振腔引起的相位滞后。一个2P2Z补偿器正好可以提供两个零点来提升中频段增益、补偿相位，同时用两个极点来抑制高频噪声和限制带宽。它在提供足够补偿能力的同时，结构又不过于复杂。
• 实现简单：2P2Z的差分方程形式清晰，在MCU中只需几次乘加运算即可完成一次控制计算，对计算资源要求低，中断服务程序执行时间短。
• 调参直观：零极点的位置直接对应到补偿器的系数，物理意义明确。通过调整零点来补偿对象相位，调整极点来塑造带宽和噪声抑制，工程师的调试直觉可以很好地映射到参数上。

注意：这里说的2P2Z是数字补偿器的实现结构。在设计时，我们通常在连续域（s域）根据环路指标设计好补偿器传递函数，然后通过离散化方法（如双线性变换）将其转换为离散域（z域）的2P2Z形式。切勿混淆设计域和实现域。

2.3 整体设计流程预览我们的设计将遵循一个清晰的流程：1) 建立LLC谐振变换器的状态空间平均模型，推导出控制到输出的传递函数 Gvd(s)。2) 在连续域设计满足环路带宽、相位裕度等指标的补偿器 Gc(s)。3) 将 Gc(s) 离散化为 Gc(z)，并整理成2P2Z的差分方程形式。4) 在MCU/DSP中编程实现该差分方程，并集成到数字脉宽调制（DPWM）闭环中。5) 通过仿真和实验验证环路性能。

3. 从物理电路到数学模型：LLC状态空间平均建模详解

这是整个项目最核心的理论基础。我们以一个最典型的半桥LLC谐振变换器为例，其功率级包括输入电压Vin，谐振电感Lr，谐振电容Cr，变压器励磁电感Lm，理想变压器（匝比n:1），以及输出滤波电容Co和负载电阻Ro。

3.1 状态变量选取与平均模型推导

对于LLC，我们通常关注其低频动态（远低于开关频率），因此输出滤波电容Co两端的电压Vo和流过Co的电流（或负载电流）是关键。但为了得到控制（开关频率fs）到输出（Vo）的传递函数，我们需要建立包含谐振腔在内的模型。一个经典且有效的方法是采用“扩展描述函数法”或“基波近似法”的简化状态空间平均。

这里我们采用一种更直观的方法：将LLC谐振网络在开关频率fs处等效为一个交流阻抗，再通过整流和滤波得到直流输出。但为了状态空间模型，我们定义以下状态变量：

• x1 = Vcr (谐振电容电压)
• x2 = Ilm (励磁电感电流)
• x3 = Vo (输出电压) —— 这是我们的主要控制目标。 实际上，对于电压模式控制，我们常忽略谐振状态的细节，直接对输出滤波器建模。但更精确的模型需要考虑谐振网络产生的等效增益。

一个实用的建模步骤是：

1. 确定工作模式：LLC通常工作在感性区域（fs < fr， fr为谐振频率），以确保主开关管实现零电压开通（ZVS）。
2. 推导直流增益：首先推导出LLC的直流电压增益公式 M = Vo/(Vin/2n) = f(fn, Q, m)。其中fn=fs/fr， Q为品质因数， m=Lm/Lr。这个公式描述了稳态下输入、输出与控制频率的关系。
3. 小信号线性化：对直流增益公式 M 在某个稳态工作点（如额定输入电压、半载）附近进行小信号扰动分析。令 fs = Fs + \hat{fs}， Vo = Vo + \hat{vo}， 其中^表示小信号量。通过泰勒展开并忽略高阶项，可以得到控制量（开关频率扰动 \hat{fs}）到输出量（输出电压扰动 \hat{vo}）的小信号传递函数 Gvd(s) = \hat{vo}(s) / \hat{fs}(s)。
4. 结合输出滤波器：将上一步得到的 Gvd(s) 视为一个“压控电压源”，其后面接上输出LC滤波器（Lo, Co）。这样，整个功率级的小信号模型就包含了LLR谐振增益的特性和输出滤波器的动态。

最终，我们可以将系统表达为标准的状态空间形式：

dx/dt = A * x + B * u y = C * x + D * u

其中，状态向量x通常包含输出电感电流（如果有）和输出电容电压；输入u是开关频率的小信号扰动；输出y是输出电压的小信号扰动。矩阵A, B, C, D的系数由电路参数（Lr, Cr, Lm, Co, Ro等）和工作点决定。

3.2 建模中的关键注意事项与技巧

• 工作点选择：LLC的模型参数（尤其是增益和相位）强烈依赖于工作点（输入电压和负载）。必须选择一个最具代表性的工作点进行建模和初始控制器设计，通常选择额定输入电压和半载或满载。设计好的控制器需要在其他工作点进行鲁棒性校验。
• 模型验证：在深入设计补偿器之前，务必用仿真工具（如PLECS、Simplis或MATLAB/Simulink）验证你推导出的状态空间模型或传递函数。在时域里给一个频率小阶跃，对比详细开关模型和平均模型的输出响应，确保两者在低频段吻合良好。
• 参数提取：模型中需要的Lr, Lm, Cr等参数，必须使用LCR表在合适的测试条件下（如频率接近谐振频率）实际测量得到，而不是仅仅依赖器件数据手册的标称值。变压器的漏感（作为Lr的一部分）尤其需要精确测量。

4. 连续域补偿器设计与离散化

拿到准确的功率级传递函数 Gp(s) = Gvd(s) * Gfilter(s) 后，我们就可以在s域设计补偿器 Gc(s)了。

4.1 环路指标确定与补偿器结构选择

首先明确设计目标：

• 穿越频率 Fc：即环路增益为0dB的频率。对于LLC，一般取开关频率的1/20到1/10。例如，开关频率为100kHz，Fc可取5kHz。太高易受开关噪声影响，太低动态响应慢。
• 相位裕度 PM：通常目标为45°到60°，以保证足够的稳定性和良好的动态响应。
• 增益裕度 GM：通常要求大于10dB。

画出 Gp(s) 的波特图。你会发现，在穿越频率附近，Gp(s) 的相位下降很快（由于谐振腔和输出滤波器的双重作用），相位裕度可能为负（系统不稳定）。因此，我们需要一个相位超前补偿器来提升相位。一个在s域传递函数为以下形式的2P2Z补偿器非常适合：

Gc(s) = Kc * (1 + s/ωz1) * (1 + s/ωz2) / [ (1 + s/ωp1) * (1 + s/ωp2) ]

其中，Kc是直流增益，ωz1和ωz2是两个零点的角频率，ωp1和ωp2是两个极点的角频率。

设计策略：

1. 将两个零点 ωz1, ωz2 放置在功率级主极点（通常是输出滤波器引起的）附近或略低的位置，以抵消其相位滞后，提升中频增益。
2. 将第一个极点 ωp1 放置在ESR零点（如果有）或远高于穿越频率的位置，主要用于积分作用，保证低频高增益以抑制稳态误差。
3. 将第二个极点 ωp2 放置在穿越频率的2-4倍处或接近谐振频率，用于衰减高频开关噪声，防止其被放大。
4. 调整增益 Kc，使得补偿后的开环传递函数 Gc(s)*Gp(s) 在设定的穿越频率 Fc 处增益为0dB。

4.2 离散化：从s域到z域

数字控制器运行在离散时间域，我们需要将设计好的 Gc(s) 转换为 Gc(z)。最常用的方法是双线性变换（Tustin变换）：

s = (2/Ts) * (z-1)/(z+1)

其中，Ts是数字控制系统的采样周期。将上式代入 Gc(s)，经过整理，即可得到z域的传递函数，其一般形式为：

Gc(z) = (b0 + b1*z^-1 + b2*z^-2) / (1 + a1*z^-1 + a2*z^-2)

这正是2P2Z结构！其中，b0, b1, b2, a1, a2 就是我们将要在MCU中使用的系数。

离散化过程中的核心要点：

• 采样频率选择：采样频率 Fs = 1/Ts。为了满足香农采样定理并减少相位延迟，采样频率通常为开关频率的2倍或更高。对于电压模式LLC，常在开关周期中点或结束时采样输出电压，采样频率等于开关频率也是常见的，但需注意由此带来的额外相位滞后。
• 频率预畸变：双线性变换会将s域的频率响应扭曲。对于关心的穿越频率Fc，需要进行预畸变校正：ω_warped = (2/Ts) * tan(ω*Ts/2)，在设计Gc(s)时，针对的零点极点频率应使用ω_warped。很多现代控制设计工具（如MATLAB的c2d函数）会自动处理这一点。
• 系数量化效应：计算出的系数b0, b1, b2, a1, a2通常是浮点数。在定点MCU中，需要将其量化为定点数（如Q格式）。量化会引入误差，可能影响极点位置，甚至导致不稳定。设计时需要评估系数精度要求，并可能在软件中采用更高精度的计算（如32位累加）。

5. 2P2Z数字补偿器的MCU实现与编程

理论模型和系数都已就绪，现在进入实战环节——用代码实现这个2P2Z控制器。

5.1 差分方程与算法实现

由z域传递函数Gc(z) = Y(z)/E(z) = (b0 + b1*z^-1 + b2*z^-2) / (1 + a1*z^-1 + a2*z^-2)，可以写出对应的差分方程：

y[n] = b0 * e[n] + b1 * e[n-1] + b2 * e[n-2] - a1 * y[n-1] - a2 * y[n-2]

其中：

• e[n]是当前采样时刻的误差信号，e[n] = Vref - Vout_ADC[n]。Vref是数字化的电压参考值，Vout_ADC[n]是ADC采样得到的输出电压值。
• y[n]是当前计算出的控制器输出（即控制量，对于LLC频率控制，就是开关周期寄存器值Tpwm的调整量ΔD或直接是频率指令）。
• e[n-1],e[n-2],y[n-1],y[n-2]是过去两个时刻的误差和输出值，需要存储在全局变量或静态变量中。

在MCU的中断服务程序（通常是PWM周期中断或ADC转换完成中断）中，执行以下步骤：

1. 采样：读取ADC结果，得到Vout_ADC[n]。
2. 计算误差：e[n] = Vref - Vout_ADC[n]。
3. 执行2P2Z计算：按照上述差分方程计算y[n]。注意运算顺序和数据类型，防止溢出。
4. 输出限幅：对y[n]进行限幅处理，确保其不超出DPWM模块允许的最大最小值范围，防止积分饱和或产生危险的频率指令。
5. 更新PWM寄存器：将y[n]转换为具体的PWM比较值或周期值，写入相应的硬件寄存器。
6. 数据更新：为下一次计算做准备：e[n-2] = e[n-1];e[n-1] = e[n];y[n-2] = y[n-1];y[n-1] = y[n]。

5.2 编程实战技巧与避坑指南

• 定点数运算与Q格式：如果使用定点MCU（如很多DSP），强烈建议使用Q格式数进行运算。例如，使用Q15格式（1位符号位，15位小数位）来表示小数系数和中间变量。乘法后需要右移15位来对齐。务必注意乘加运算过程中的累加器位宽（通常32位），防止溢出。

  // 示例：Q15格式下的2P2Z计算（伪代码） #define B0_Q15 (int16_t)(0.1234 * 32768) // 将系数转换为Q15整数 #define B1_Q15 (int16_t)(0.5678 * 32768) // ... 其他系数 int32_t acc; // 32位累加器 acc = (int32_t)B0_Q15 * e_n; acc += (int32_t)B1_Q15 * e_n_1; acc += (int32_t)B2_Q15 * e_n_2; acc -= (int32_t)A1_Q15 * y_n_1; acc -= (int32_t)A2_Q15 * y_n_2; y_n = (int16_t)(acc >> 15); // 结果转换回Q15

• 抗积分饱和（Anti-Windup）：当输出因限幅而饱和时，积分项（在2P2Z中体现在递归反馈项上）可能会持续累积（windup），导致系统退出饱和时产生大的超调。简单的抗饱和方法是：在计算y[n]后，如果进行了限幅，则将限幅后的y[n]值回代，反算出对应的“正确的”历史输出值y[n-1]（或冻结积分项），用于下一次迭代。这能有效改善大信号动态响应。
• 采样同步与延迟补偿：从采样到PWM更新存在计算延迟和PWM硬件更新延迟。这个延迟会在环路中引入额外的相位滞后。在建模时，可以将其近似为一个纯延迟环节e^{-sTd}，其中Td约为1.5个采样周期。在设计补偿器时，需要将这个延迟考虑进去，预留更多的相位裕度。或者，采用预测性控制算法来主动补偿这部分延迟。
• 参考值软启动：上电时，不要直接将Vref设置为目标值。应采用一个从0缓慢斜坡上升到目标值的软启动过程，避免启动冲击电流。这可以通过在中断中逐步增加Vref的值来实现。

6. 环路调试、验证与常见问题排查

代码烧录后，真正的挑战才开始。你需要验证环路的实际性能。

6.1 调试工具与方法

• 波特图测试：这是最权威的验证方法。需要一台网络分析仪或具备波特图功能的示波器（如Keysight InfiniiVision系列带MSO的型号）。将一个小幅度的正弦扫频扰动信号注入环路（通常在误差放大器输出或DA注入点），测量注入点前后信号的幅值和相位差，即可绘制出开环波特图。实测的穿越频率、相位裕度、增益裕度应与设计值吻合。
• 动态负载测试：使用电子负载对输出进行阶跃负载变化（如从25%跳到75%负载）。用示波器观察输出电压的瞬态响应。记录下冲/下冲电压、恢复时间、振荡次数。一个好的环路应表现为快速、平滑的恢复，超调小。
• 输入电压阶跃测试：改变输入电压，观察输出电压的调整过程，评估系统的输入抗扰度。

6.2 常见问题与解决方案速查表

6.3 调试心得：模型与现实的桥梁

在实际调试中，你可能会发现，无论模型多么精确，实测结果与理论总有偏差。这非常正常。数字补偿器设计的优势就在于，你可以快速修改几个系数，在线观察效果。

我的个人经验是：将理论计算出的系数作为初始值，而不是最终值。首先确保系统能在所有工作点稳定（不振荡），然后通过波特图测试，微调补偿器零极点位置。例如，如果相位裕度不够，可以尝试同时将两个零点向高频移动一小段（比如提高20%），并相应调整增益。每次修改后，重新进行动态负载测试，观察时域响应是否改善。这个过程是连接理论模型和物理世界的桥梁，也是积累“手感”的关键。

最后，别忘了进行应力测试：在最高温、最低温、最高输入电压、最低输入电压、满载、空载等各种极端组合下测试系统，确保环路在任何情况下都稳定可靠。数字补偿器的参数如果存储在非易失性存储器中，甚至可以针对不同的典型工作条件存储多组参数，实现自适应的最优控制，这正是其超越模拟补偿器的潜力所在。