clear,clc,
u=0;
u0=0;
i=0;
t=0;
ph0=0.57;
R=5;
f=50;
L=0.05;
Um=110;
dt=1e-4;
N=3000;
t0=0.1 ;

regim=2;
for k=1:N
t=[t t(k)+dt];
if t(k) >=t0
if regim==1
u= [u Um];
else
u=[u Um*sin(2*pi*f*t(k)+ph0)];
end
else
u=[u u0];
end
di=(u(k)-i(k)*R)/L*dt;
i=[i i(k)+di];
end
plot(t,i);
Соседние файлы в папке файлы матлаб