Graphs_kater

clc;clear;
% параметры заданного судна
W=25.5; %водоизмещение в тоннах
N=1000; %мощность ГЭУ в лс
v=30; %скорость в узлах
%per=1852/3600;

% перевод в систему СИ
M=1000*W %масса катера
Nmax=735.5*N
vmax=0.51*v
Fmax=Nmax/vmax %максимальная сила тяги движителя
dFmax=0.1*Fmax %допустимое изменение силы тяги

dt=1; %шаг интегрирования
dP=dt*(dFmax/Fmax)*100 %относительное допустимое изменение силы тяги за время равное 1с
A=Fmax/(vmax^2)

% построение графиков
xi=0; x_pred=0; P=0; t=0;
global time; global x; global vv; global p;
x=[]; time=[]; vv=[]; p=[];

while P<100 % разгон + выход на максимальное значение силы тяги
x_sled=xi+(xi-x_pred)+(P*Fmax*(dt^2)/100-A*(xi-x_pred)*abs(xi-x_pred))/M;
v_sled=(x_sled-xi)/dt;
vv=[vv,v_sled];
p=[p,P];
x_pred=xi; xi=x_sled;
P=P+dP;
time=[time,t];
x=[x,x_sled];
t=t+dt;
end;
while v_sled<0.9999*vmax % разгон на максимальной тяге до максимальной скорости
x_sled=xi+(xi-x_pred)+(P*Fmax*(dt^2)/100-A*(xi-x_pred)*abs(xi-x_pred))/M;
v_sled=(x_sled-xi)/dt;
vv=[vv,v_sled];
p=[p,P];x_pred=xi; xi=x_sled;
time=[time,t];
x=[x,x_sled];
t=t+dt;
end;

while P>-100 % торможение, выход на обратную силу тяги
x_sled=xi+(xi-x_pred)+(P*Fmax*(dt^2)/100-A*(xi-x_pred)*abs(xi-x_pred))/M;
v_sled=(x_sled-xi)/dt;
vv=[vv,v_sled];
x_pred=xi; xi=x_sled;
p=[p,P];
P=P-dP;
time=[time,t];
x=[x,x_sled];
t=t+dt;
end;
while v_sled>=0.001 % торможение, до нуля
x_sled=xi+(xi-x_pred)+(P*Fmax*(dt^2)/100-A*(xi-x_pred)*abs(xi-x_pred))/M;
v_sled=(x_sled-xi)/dt;
vv=[vv,v_sled];
p=[p,P];
x_pred=xi; xi=x_sled;
time=[time,t];
x=[x,x_sled];
t=t+dt;
end;
tinv=time';
pinv=p';
vinv=vv';
xinv=x';
subplot(3,1,1)
plot(time,x);grid on; xlabel 't, c'; ylabel 'x, м'; title 'Расстояние';
subplot(3,1,2)
plot(time,vv);grid on; xlabel 't, c'; ylabel 'v, м/с'; title 'Скорость';
subplot(3,1,3)
plot(time,p);grid on; xlabel 't, c'; ylabel 'P, %'; title 'Относительная сила тяги';
time=time';
x=x';
vv=vv';
p=p';