Уравнение движения тела с парашютом. Траектория и график - Matlab - Страница 2

@Fabien159 · 20.04.2020, 12:26 **[ТС]**

Студворк — интернет-сервис помощи студентам

SSC, понял. Спасибо.
Подскажите, а через dsolve можно решить эту систему, что вы написали?

Matlab M
1
2
3
4
dhdt(1,1)= x(2);
dhdt(2,1)=(-m*g-k*x(2)*Vnrm)/m;
dhdt(3,1)= x(4);
dhdt(4,1)=(-k*x(4)*Vnrm)/m;

@SSC · 20.04.2020, 13:15

Сообщение от Fabien159

а через dsolve можно решить эту систему,

Скорее всего никак, уравнения слишком сложны, чтобы получить аналитическое решение.
Однако если требуется получить функцию зависимости расстояния от точки сброса до точки падения от высоты и скорости, то поступают след. образом …..(дальнейшая информация не приводится, что бы меня не обвинили в теоретической подготовке разработчиков бомбоприцелов )

@Fabien159 · 20.04.2020, 17:24 **[ТС]**

SSC, я понял. Я как-то создавал тему, я там пытался решить эту задачу через цикл.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
function raschett
clear; clc;
delt=0.1; %дельта времени
v=200; %начальная скорость
m=13; % вес с парашютом
Cx=0.2; %кф
g=9.81; %ускорение
%h=2000; %высота сброса
s=0; %начальное растояние пролета
ma=m*g-Cx*v^2
t=0;
a=g-((Cx*v^2)/m)
 
n = 200;
 
for i = 1:n-1
    v(i+1)=v(i)+a(i)*delt; %скорость 
    a(i+1)=g-((Cx*v(i)^2)/m); %ускорение
    s(i+1)=s(i)+v(i)*delt+((a(i)*delt^2)/2); %горизонтальный снос
    t(i+1)=t(i)+delt; %время 
  
        
end
plot(t, s)

как думаете, таким методом можно построить траекторию? Если доработать.

@SSC · 21.04.2020, 07:47

Сообщение от Fabien159

Если доработать.

Если доработать, то да.
Вы пытаетесь реализовать метод Эйлера первого порядка.
У него точность при том же размере шага расчета значительно меньше чем у Рунге-Кутта 4-го порядка.
Он будет всегда давать завышенные результаты по дальности.
Но конечно начинать (учиться) надо действительно с Эйлера

@Fabien159 · 21.04.2020, 16:23 **[ТС]**

SSC, спасибо!
Единственное не понимаю еще, как получить высоту таким же путем...

в другой теме вы мне отвечали уже.

Расчет надо вести для 4-ех переменных - высота (h), вертикальная скорость (Vh), горизонтальное перемещение (x), горизонтальная скорость (Vx).
В каждый момент времени вам известны эти величины.
На каждом шаге Вы вычисляете
1. Скорость тела V=sqtr(Vh^2+Vx^2)
2. Силу трения F=k*V^2
3. Раскладываете силу F на составляющие по горизонтали и вертикали.
4. Рассчитываете изменение Vh и Vx
5. Рассчитываете изменение h и x

Вы написали, что в каждый момент времени мне известны эти величины. Но они же не известны?

@SSC · 22.04.2020, 07:09

Сообщение от Fabien159

Но они же не известны?

Известны величины параметров в текущий момент времени:
1 высота
2 скорость падения(подъема)
3 горизонтальное смещение
4 горизонтальная скорость
Для текущего момента времени по представленным формулам ( например как в сообщении 21) Вы рассчитываете производные этих величин. Зная производные находите величины параметров в следующий момент времени.
В методе Эйлера это просто в явную видно. В более сложных методах в сущности это тоже так, но выполняется еще ряд действий для повышения точности метода.

Новые блоги и статьи Все статьи Все блоги /
Запустил конкурс "тем и промптов для текстовых квестов созданных почти чисто ИИ" Adler 06.10.2026 Всем привет! За последние три-четыре дня я создал более 16 текстовых квестовых игр используя преимущественно по одному запросу к ИИ на игру. Мне так понравилось смотреть все ветки/ сцены во всех. . .	ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .	Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .
Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .	Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .	Запрет дублирования строк в табличной части Maks 13.09.2026 Реализация из решения ниже выполнена на нетиповом справочнике "Нормы ТО" с табличной часть "Виды ТО", разработанного в КА2, со следующими реквизитами: - ВидТО (СправочникСсылка. ВидыТО); - ВидГСМ. . .	Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Jin X 06.09.2026 Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Работая с форумом и нейросетями в браузере часто хочется что-то подкорректировать или добавить какого-то функционала. Ниже прикреплён. . .