Курсовая работа: Розрахунок трьохфазного мостового випрямляча

[s,Z]=ode113(@F_first, [0 (Dtx(3)-Dtx(2))+10^-8], Y0);

s=s+t(length(t));

Y=[Y;Z];

t=[t;s];

n=0;

Ox=0;

Px=0;

for x=0:0.000001:(Dtx(3)-Dtx(2))

n=n+1;

Ox(n)=E*sin(w*x+PhaseShift);

Px(n)=x;

end

Px=Px+P(length(P));

P=[P, Px];

O=[O, Ox];

PhaseShift=PhaseShift+(Dtx(3)-Dtx(2))*w;

Dtx=Dtx+Dt(length(Dt));

D=[D, Dx];

Dt=[Dt, Dtx];

end

for x=1:length(t)

Ufb(x)=(Y(x, 2)*5/36-Uref)*Kfb;

if (Ufb(x)<0)

Ufb(x)=0;

end

end

I=Y;

I(:, 2)=[];

subplot(4,1,1);