Програмирование физически процесов

Задача о теле брошенном под углом к горизонту дан угол альфа начальная скорость и сопротивление среды масса обьекта как не решая дифференциальные уравнения тоесть не имея общего решения построить траекторию?

0

t -время , его в цикл с мелким шагом от 0 до того сколько времени надо
alpha - нач угол в радианах,
v - начальная скорость

1

так это уже решение диф уравнения я спросил на счет вычисления производной не решая уравнения

0

@aLarman, а где

Сообщение от qsv1991 сопротивление среды

?

0

это мои наработки

Миниатюры

 

0

И вот еще

Миниатюры

 

0

Ребята есть идеи?

0

@qsv1991, как люди должны что-то понять в этом? Перепечатайте, в расширенном режиме ввода сообщения и редактор формул есть, т.ч. все возможно написать вполне аккуратно.

0

так в начале составим диф. уравнение c-сопротивление m- масса g -ускорение свободного падения
m*x"=-c*x` (1)
m*y"=-m*g-c*y`
решение этого диф уравнения есть(v0-начальная скорость)
x(t)=(m/c)*v0*Cos(alfa)(1-exp(-(c*t)/m)) (2)
y(t)=(m/c)*(v0*Sin(alfa)-m*g)*(1-exp(-(c*t)/m)-(m*t*t)/c)
но задача которую поставил передо мной преподователь другая.Он поставил задачу чтобы найти общий метод решения диф ур численно не находя обшего решения Вот.
Зная что
Vx=x`
Vy=y` (3)
подставим (3) в (1) получим
Vx`=-(c*Vx)/m (1)
Vy`=-g-(c*Vy)/m
теперь по правилу производно (T - Delta t)
(Vx(t+T)-Vx(t))/T=-(c*Vx(t))/m
(Vy(t+T)-Vy(t))/T=-g-(c*Vy(t))/m в рез-тате преобразований получим

Vx(t+T)-Vx(t)=T*(-c*Vx(t))/m
Vy(t+T)-Vy(t)=T*(-g-(c*Vy(t))/m)

Vx(t+T)=T*(-c*Vx(t))/m +Vx(t)
Vy(t+T)=T*(-g-(c*Vy(t))/m) +Vy(t)

Vx(t+T)=Vx(t)*(1-T*(c/m) )
Vy(t+T)=-T*g +Vy(t)*(1-T*(c/m))

а если загнать в цикл последнюю систему при малом T и при 0<t< касания с поверхностью мы получим меняющуюся со временем скорости по Х и У а как получить координаты?

Добавлено через 19 часов 6 минут
Ребята ну хоть прокоментируйте а

0