MATLAB符号工具箱求解RLC电路微分方程报错排查

报错原因 #

你写的第二个初始条件cond2 = diff(i(t),t) == V0/L;没有指定t=0时刻的取值，dsolve会把它当成一个额外的微分方程，而非初始条件。原微分方程已经是1个未知数i(t)的二阶方程，加上这个错误的"方程"后，方程数（2个）远超未知数数量（1个），直接触发"无法化简方阵系统"的报错。

修正步骤 #

1. 修正初始条件：把导数初始条件明确限定在t=0时刻，MATLAB R2019b及以上版本支持diff(i(t),t)(0)的语法，老版本可以用subs(diff(i(t),t), t, 0)。
2. 拆分通解与特解求解：先不带初始条件求通解，再带入初始条件求特解，逻辑更清晰。

完整可运行代码 #

syms i(t) R L C V0 

% 定义RLC元件参数
R = 10;
L = 1;
C = 1;
V0 = 12;

% 建立二阶微分方程（零输入响应）
eqn = L*diff(i(t), t, 2) + R*diff(i(t), t) + (1/C)*i(t) == 0;

% 定义正确的初始条件
cond1 = i(0) == 0;              % t=0时电流为0
cond2 = diff(i(t), t)(0) == V0/L; % t=0时电流导数为V0/L

% 求解通解（不含初始条件）
general_sol = dsolve(eqn);
disp('i(t)的通解：')
disp(general_sol)

% 求解特解（带入初始条件）
particular_sol(t) = dsolve(eqn, cond1, cond2);
disp('i(t)的特解：')
disp(particular_sol(t))

额外说明 #

你的RLC参数是过阻尼配置（阻尼比ζ=R/(2*sqrt(L/C))=5>1），特解会是两个指数衰减项的和，运行代码后可以看到具体表达式。如果需要可视化，加一行fplot(particular_sol(t), [0, 5])就能画出暂态响应曲线。

内容的提问来源于stack exchange，提问作者BP130