二阶电路零状态响应的另一种实现方式

RLC1.m文件

%RLC状态方程
function dydt = rlc1(t,y)
dydt=zeros(size(y));
us=1500;
r=0.1;
l=0.01;
c=0.003;
%y(1)=uc,y(2)=ic=il,ic=c*duc/dt,ul=l*dil/dt=us-r*ic-uc
%dydt =[y(2)/c;(us-r*y(2)-y(1))/l];
dydt(1)=y(2)/c;
dydt(2)=(us-r*y(2)-y(1))/l;

RLC2.m文件

clear all
clc
close
t0=0.01;
t1=1;
y0=[0;0];
[t,y]=ode45('rlc1',t0,t1,y0);
figure(1);
subplot(211);
plot(t,y(:,1));
title('RLC二阶电路零状态响应 电压波形');
xlabel('时间t（s）');
ylabel('电压U（V）');
grid on;
subplot(212);
plot(t,y(:,2));
title('RLC二阶电路零状态响应 电流波形');
xlabel('时间t（s）');
ylabel('电流I（A）');
grid on;