RLC.m
clear all;
clc;
u=dsolve('0.00003*D2y+0.0003*Dy+y=1500','y(0)=0,Dy(0)=0','t'); %求RLC二阶电路零状态响应,LC(d2u/dt)+RC(du/dt)+u=us
syms t;
i=diff(u);
i=inline(i);
t=linspace(0.01,1,1000);
for k=1:1000
y(k)=0.003*i(t(k));
end
plot(t,y);
RLC.m
clear all;
clc;
u=dsolve('0.00003*D2y+0.0003*Dy+y=1500','y(0)=0,Dy(0)=0','t'); %求RLC二阶电路零状态响应,LC(d2u/dt)+RC(du/dt)+u=us
syms t;
i=diff(u);
i=inline(i);
t=linspace(0.01,1,1000);
for k=1:1000
y(k)=0.003*i(t(k));
end
plot(t,y);