ВУЗ:
Составители:
93
G_p(i-1)=0;
else G_p(i-1)=mu1*sqrt(p(i-1)^2-p_pot^2);
end
G_pi=G_p(i-1)+(M_p(i)-M_p(i-1))/dt;
if G_pi<=0
G_pi=0;
end
M_v(i)=M_v(i-1)+dt*(G_v-G_pi);
G_sm=mu2*sqrt(p_t^2+p_k^2);
if p_t < p_k
G_sm=0;
end
t_k(i)=(K*F*t_v(i-1)+G_sm*(c_k*t_t_vx+r_t))/(K*F+G_sm*c_k);
t_v(i)=t_v(i-1)+dt*(K*F*(t_k(i)-t_v(i-1))-G_v*c_v*(t_v(i-1)-t_v_vx)-
G_pi*Lamda)/(c_v*M_v(i));
end
% построение графиков зависимостей
figure(1)
plot(t,t_v);
figure(2)
plot(t,t_k);
figure(3)
p(end)=p(end-1);
plot(t,p);
figure(4)
G_p(end)=G_p(end-1);
plot(t,G_p);
figure(5)
plot(t,M_v);