何必手算？MATLAB解析法求解微分方程|Matlab符號計(jì)算
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;
標(biāo)簽：
