何必手算？MATLAB解析法求解微分方程|Matlab符号计算

clc

close all

clear

syms L C G iL(t) Is

eqn = L*C*diff(iL, t, 2)+G*L*diff(iL, t)+iL == Is;

DiL = diff(iL, t);

cond = [iL(0)==0, DiL(0)==0];

iL = dsolve(eqn, cond);

iL = subs(iL, [C G L Is], [1e-6 2e-3 1.001 1]);

fplot(iL, [0 0.008], 'LineWidth', 2);

grid on;

diL = diff(iL, t);

fplot(diL, [0 0.008], 'LineWidth', 2);

grid on;

ddiL = diff(diL, t);

fplot(ddiL, [0 0.008], 'LineWidth', 2);

grid on;