4. 实时仿真建模基础

在电源控制器 HIL 验证中，被控对象（如 Buck、Boost、Buck-Boost 变换器）不能直接以真实开关器件模型接入实时仿真器。原因很直接：若按开关频率 100 kHz 建模，实时步长需远小于 10 μs，FPGA 或 CPU 难以在有限资源下完成多路变换器的高保真解算。因此，工程上普遍采用开关周期平均模型，把开关纹波“平均掉”，只保留低频包络动态，再用状态空间平均法建立连续时间模型，最后通过离散化变成实时求解器可执行的差分方程。

一、状态空间平均法的基本思路

以理想 Buck 变换器为例，设电感电流 i_L、电容电压 v_C 为状态变量，输入电压 V_g，占空比 d。在一个开关周期 T_s 内，开关管导通时段状态方程为：

L·di_L/dt = V_g − v_C
C·dv_C/dt = i_L − v_C/R

关断时段为：

L·di_L/dt = −v_C
C·dv_C/dt = i_L − v_C/R

将两段按占空比 d 加权平均，得到平均状态方程：

L·di_L/dt = d·V_g − v_C
C·dv_C/dt = i_L − v_C/R

写成标准形式 ẋ = A·x + B·u，其中 x = [i_L, v_C]^T，u = [V_g, d]^T。注意 d·V_g 是非线性乘积项，在稳态工作点附近做小信号线性化，令 d = D + d̂、V_g = V_g + v̂_g，忽略二阶小量后得到小信号模型，这才是控制器设计常用的对象。

二、开关周期平均模型与实时实现

平均模型在仿真中通常以“受控源 + 平均开关”形式实现：用占空比 d 控制一个电压源或电流源，替代真实 MOSFET 和二极管。这样模型在连续域是平滑的，求解器可以用较大步长。但要注意，平均模型不包含开关纹波、死区、体二极管反向恢复等细节，因此它适合验证环路稳定性、负载阶跃响应、输入扰动抑制，不适合验证开关尖峰、EMI、驱动时序。

在 HIL 中，平均模型常与 PWM 生成器配合：控制器输出占空比指令，HIL 侧用该指令驱动平均模型，同时把电感电流、输出电压反馈回控制器。这样闭环就能在实时步长下运行。

三、离散化方法的选择

连续模型必须离散成差分方程才能被实时求解器执行。常用方法有三种：

• 前向欧拉：x[k+1] = x[k] + T_s·(A·x[k] + B·u[k])。实现最简单，但稳定性差，步长稍大就发散。
• 后向欧拉：x[k+1] = x[k] + T_s·(A·x[k+1] + B·u[k+1])。需要解隐式方程，但数值稳定性好，适合刚性系统。
• Tustin（双线性）：s ≈ 2/T_s·(z−1)/(z+1)。频域保真度高，常用于控制器离散化，但可能引入频率预畸变。

以 Buck 参数 L = 100 μH、C = 470 μF、R = 5 Ω、T_s = 10 μs 为例，LC 谐振频率约 f₀ = 1/(2π√(LC)) ≈ 734 Hz。前向欧拉要求 T_s < 2/|λ_max|，而系统极点实部约 −1000 rad/s 量级，10 μs 步长下前向欧拉仍可稳定，但若步长增至 100 μs 就会发散。后向欧拉在同样步长下仍稳定，代价是相位滞后略大。

四、求解器选择与数值稳定性

实时仿真器通常提供定步长求解器：显式欧拉、RK2/RK4、后向欧拉、Tustin。选择依据是：

• 若模型刚性不强、步长小，RK4 精度最高，但计算量大；
• 若模型含快速动态（如寄生参数、电流环带宽高），后向欧拉或 Tustin 更稳；
• 若要求严格实时，定步长显式方法必须满足稳定性条件，否则宁可降步长或换隐式方法。

数值稳定性还要关注量化误差和代数环。HIL 中 ADC 量化、PWM 分辨率、定点运算都会引入极限环。若控制器输出直接反馈到平均模型输入，可能形成代数环，需插入单位延迟或使用隐式求解。

五、工程实践中的关键参数与判断

在 HIL 项目中，建议按以下顺序确定模型参数：先根据实物测量或数据手册确定 L、C、R、V_g；再计算 LC 谐振频率 f₀；然后选择实时步长 T_s，一般要求 T_s ≤ 1/(20·f₀)，同时满足求解器稳定性条件；最后用阶跃响应对比平均模型与详细开关模型的包络，若两者在低频段吻合，说明平均模型有效。

若发现 HIL 输出出现持续振荡或发散，优先检查：步长是否超过稳定性边界、离散化方法是否与模型刚性匹配、是否存在代数环、定点量化是否引入极限环。必要时改用后向欧拉或 Tustin，并在反馈路径插入一个采样周期延迟。

总结来说，状态空间平均法解决“建模”问题，开关周期平均模型解决“实时性”问题，离散化方法解决“可执行”问题，求解器选择与数值稳定性解决“可信”问题。四者环环相扣，缺一不可。