子函数Tworank

function yd=Tworank(t,y)

R1=20;

R2=10;

R3=10;

C=100e-6;

L=0.5;

yd=[-y(2)/C;-(R1+R2*R3/(R2+R3))/L*y(2)+y(1)/L];

end

主函数erjie

clc

clear all

tspan=[0,0.2];

y0=[25;5];%y0和y0'的初始值

[t,YY]=ode45('Tworank',tspan,y0);

subplot(2,1,1),plot(t,YY(:,1),'color',[0.8500 0.3250 0.0980],'linewidth',2),

title('电容电压');

legend('uct')

xlabel('t/s')

ylabel('uct')

subplot(2,1,2),plot(t,YY(:,2),'-b','linewidth',2),

title('电感电流');

legend('i2t')

xlabel('t/s')

ylabel('i2t')