Курсовая работа: Моделирование движения парашютиста

k2=0.5*1.2*0.4*1.225

g=9.81; % ускорение свободного падения

m=50; % масса манекена

h0=5000; % высота

[t h]= ode45(@parashut,[0 200],[h0 0] )

r=find(h(:,1)>=0);

s=length(r);

b=length(t);

h(s+1:b,:)=[];

t(s+1:b,:)=[];

a=g-(k1*-h(:,2)+k2*h(:,2).*h(:,2))/m % вычисляемускорение

% Построение графика зависимости высоты от времени

subplot(3,1,1), plot(t,h(:,1),'LineWidth',1,'Color','r'),grid on;

xlabel('t, c'); ylabel('h(t), m');

title('Графикзависимостивысотыотвремени', 'FontName', 'Arial','Color','r','FontWeight','bold');

legend('m=50 kg')

% Построение графика зависимости скорости от времени

subplot(3,1,2), plot(t,h(:,2),'LineWidth',1,'Color','b'),grid on;

xlabel('t, c');

ylabel('V(t), m/c');

Title('Графикзависимостискоростиотвремени', 'FontName', 'Arial','Color','b','FontWeight','bold');

legend('m=50 kg')

% Построение графика зависимости ускорения от времени

subplot(3,1,3), plot(t,a,'-','LineWidth',1,'Color','g'),grid on;

text (145, 0,'t, c');

ylabel('a(t), m/c^2');

Title('Графикзависимостиускоренияотвремени', 'FontName', 'Arial','Color','g','FontWeight','bold');

legend('m=50 kg')

Экранная форма вывода графиков.

Список использованных источников