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

1次阅读

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

Buck 电路作为最基本的 DC-DC 拓扑，其开关管周期性通断会引入强烈的时变非线性特性。新手在尝试分析其动态行为时，常遇到两个典型问题：

开关级仿真效率低下 ：由于需要捕捉 ns 级的开关瞬态，仿真步长必须极小，导致 1ms 的电路动态可能需要数分钟计算
理论分析困难 ：直接求解开关状态的微分方程几乎不可能，难以提取系统的传递函数等关键特性

状态空间平均法通过将开关周期内的变量平均值作为状态量，巧妙地将非线性系统转化为等效线性模型。这种方法不仅能保持关键动态特性，还能将仿真速度提升 100 倍以上。

列出各开关状态方程
开关导通时段（$0<t<DT_s$）：
$$\begin{cases}
L\frac{di_L}{dt} = V_{in} – v_o \
C\frac{dv_o}{dt} = i_L – \frac{v_o}{R}
\end{cases}$$
开关关断时段（$DT_s<t<T_s$）：
$$\begin{cases}
L\frac{di_L}{dt} = -v_o \
C\frac{dv_o}{dt} = i_L – \frac{v_o}{R}
\end{cases}$$
加权平均处理
引入占空比 $D$ 作为权重系数，得到平均模型：
$$\begin{cases}
L\frac{d\langle i_L \rangle}{dt} = D V_{in} – \langle v_o \rangle \
C\frac{d\langle v_o \rangle}{dt} = \langle i_L \rangle – \frac{\langle v_o \rangle}{R}
\end{cases}$$
小信号线性化
在稳态工作点 $(I_L,V_o)$ 附近引入扰动：
$$\begin{cases}
\langle i_L \rangle = I_L + \hat{i}L \
\langle v_o \rangle = V_o + \hat{v}_o \
D = D_0 + \hat{d}
\end{cases}$$
最终得到矩阵形式的状态空间方程：
$$\frac{d}{dt}\begin{bmatrix}
\hat{i}_L \
\hat{v}_o
\end{bmatrix} =
\begin{bmatrix}
0 & -1/L \
1/C & -1/(RC)
\end{bmatrix}
\begin{bmatrix}
\hat{i}_L \
\hat{v}_o
\end{bmatrix} +
\begin{bmatrix}
V/L \
0
\end{bmatrix}\hat{d}$$

% Buck 电路参数
L = 100e-6;     % 电感量 (H)
C = 470e-6;     % 电容量 (F)
R = 5;          % 负载电阻 (Ω)
Vin = 24;       % 输入电压 (V)
D = 0.5;        % 占空比
Ts = 1/100e3;   % 开关周期 (s)

% 状态空间方程定义
A = [0 -1/L; 1/C -1/(R*C)];
B = [Vin/L; 0];
C_mat = eye(2);  % 输出全部状态量
D_mat = zeros(2,1);

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

% 使用 ode45 求解
[t,x] = ode45(@(t,x) A*x + B*D, tspan, x0);

% 结果可视化
figure;
subplot(2,1,1);
plot(t, x(:,1)); 
title('电感电流动态响应');
xlabel('时间 (s)'); ylabel('电流 (A)');

subplot(2,1,2);
plot(t, x(:,2)); 
title('输出电压动态响应');
xlabel('时间 (s)'); ylabel('电压 (V)');