Buck电路状态空间平均模型仿真实战：从理论到MATLAB实现

1次阅读

没有评论

共计 1947 个字符，预计需要花费 5 分钟才能阅读完成。

刚开始学习电力电子时，我最头疼的就是 Buck 电路的动态分析。直接用开关器件搭建仿真模型，虽然能反映真实情况，但每次仿真都要等上好几分钟。特别是做参数优化时，改个参数就要重新跑一遍，效率太低。

后来才知道，这是因为开关管高频切换导致系统具有时变非线性特性。传统开关模型需要非常小的时间步长（通常要小于开关周期的 1 /100），才能准确捕捉开关瞬态，这就带来了巨大的计算量。

状态空间平均法的妙处在于，它通过数学方法将开关电路的非线性特性 ” 线性化 ” 处理。具体来说，就是把一个开关周期内电路的工作状态取平均值，用连续模型来近似描述离散开关行为。

这个方法的数学基础其实很简单：

分别写出开关导通（$0<t<DT_s$）和关断（$DT_s<t<T_s$）时的状态方程
将两个方程按占空比 $D$ 加权平均
引入小信号扰动（$\hat{d}$ 表示占空比扰动）

最终得到的平均模型状态方程：

$$\frac{d}{dt}\begin{bmatrix} i_L \ v_C \end{bmatrix} = \begin{bmatrix} -\frac{r_L}{L} & -\frac{1}{L} \ \frac{1}{C} & -\frac{1}{RC} \end{bmatrix}\begin{bmatrix} i_L \ v_C \end{bmatrix} + \begin{bmatrix} \frac{V_{in}}{L} \ 0 \end{bmatrix}D + \begin{bmatrix} \frac{V_{in}}{L} \ 0 \end{bmatrix}\hat{d}$$

下面我用 MATLAB R2021b 演示具体实现过程。建议先创建两个文件：buck_ss_model.m 定义模型方程，run_simulation.m 运行仿真和绘图。

function dxdt = buck_ss_model(t, x, Vin, L, C, R, rL, D)
    % 状态变量: x(1)=iL, x(2)=vC
    % 电路参数
    A = [-rL/L, -1/L; 1/C, -1/(R*C)];
    B = [Vin/L; 0];

    % 状态方程
    dxdt = A * x + B * D;
end

% 电路参数
Vin = 24;     % 输入电压 (V)
L = 100e-6;   % 电感 (H)
C = 470e-6;   % 电容 (F)
R = 5;        % 负载电阻 (Ω)
rL = 0.1;     % 电感寄生电阻 (Ω)
D = 0.5;      % 占空比

% 仿真设置
tspan = [0 0.01];  % 仿真时间 10ms
x0 = [0; 0];       % 初始状态 [iL; vC]

% 调用 ode45 求解
[t, x] = ode45(@(t,x) buck_ss_model(t,x,Vin,L,C,R,rL,D), tspan, x0);

% 绘制结果
figure;
subplot(2,1,1);
plot(t, x(:,1)); title('电感电流'); xlabel('时间 (s)'); ylabel('i_L(A)');
subplot(2,1,2);
plot(t, x(:,2)); title('输出电压'); xlabel('时间 (s)'); ylabel('v_C(V)');

运行后会看到电感电流和输出电压从零开始的启动波形。如果想观察动态响应，可以修改代码在仿真中途改变占空比。

为了验证平均模型的准确性，我在 Simulink 中搭建了实际的开关电路模型作为基准。对比结果发现：