М-файл функции parashut.m:
%Функция моделирования движения парашютиста
function dhdt=parashut(t,h)
global k1 k2 g m
% система ДУ первого порядка
dhdt(1,1)= -h(2);
dhdt(2,1)=(m*g-k1*h(2)-k2*h(2)*h(2))/m
М-файл вывода результатов parashutist.m:
% Моделирование движения парашютиста
% Васильцов С. В.
clc
global h0 g m k1 k2 a
% k1-линейный коэффициент пропорциональности, определяющийся свойствами среды и формой тела. Формула Стокса.
k1=6*0.0182*0.4;
%k2-квадратичный коэффициент пропорциональности, пропорционален площади сечения тела, поперечного по
%отношения к потоку, плотности среды и зависит от формы тела.
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')
Экранная форма вывода графиков.