何必手算?Matlab求解电力电子电路微分方程|Matlab符号计算|电力电子


clc %清屏
close all %关闭所有窗口
clear %清空工作空间
%定义符号变量
syms L R Em E I1(t) I10 T I2(t) I20
%定义微分方程本体1
ode = L*diff(I1, t) + R*I1(t) + Em == E;
conds = I1(0)==I10; %定义初始条件
res1 = dsolve(ode, conds); %求解微分方程
res1 = subs(res1, L/R, T); %变量代换
res = collect(res1, I10); %合并同类项
symdisp(res) %显示微分方程结果
%定义微分方程本体2
ode = L*diff(I2, t) + R*I2(t) + Em == 0;
conds = I2(0) == I20; %定义初始条件
res2 = dsolve(ode, conds); %求解微分方程
res2 = subs(res2, L/R, T); %变量代换
res = collect(res2, [I20 Em]); %合并同类项
symdisp(res) %显示微分方程结果