файлы матлаб / lab2

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);