clc %清屏

close all %關(guān)閉所有窗口

clear %清空工作空間

%定義符號(hào)變量

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); %合并同類項(xiàng)

symdisp(res) %顯示微分方程結(jié)果

%定義微分方程本體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]); %合并同類項(xiàng)

symdisp(res) %顯示微分方程結(jié)果