Скорость объекта, который падает свободно вследствие земного притяжения, может быть смоделирована уравнением: mdvdt=-mg+kv2 где m –

Скорость объекта, который падает свободно вследствие земного притяжения, может быть смоделирована уравнением: mdvdt=-mg+kv2 где m – масса объекта, g=9,81 м/с2 и k – постоянная. Решите это уравнение относительно v для случая m=5 кг, k=0,05 кг/м, 0≤t≤15 сек и v0=0 м/с. Создайте график v как функции времени.

Перепишем уравнение в следующем виде:
dvdt=-g+kv2m
Используемые команды:
@ - создание анонимной функции
Ode45 - Решение нежестких дифференциальных уравнений
Plot – создание графика
Код матлаб:
clear;clc %ochistka sredi
m=5; g=9.81; k=0.05; %zadanie constant
[email protected](t,v)-g+k*v^2/m; %zapicb uravneniya s prisvaivaniem anonimnoi fucnkcii
[t, v]=ode45(dvdt,[0:0.1:15],0); %integrirovanie diff uravneniya v diapazone t s nachalnim usloviem v0
plot(t,v) % sozdanie grafika
xlabel('Time (s)') %podpisb x
ylabel('Velocity (m/s)') %podpisb y
Результат выполнения кода:
Рис