
Кудряшов Всеволод / 13 лаб / 13 практика
.docxОрдена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования «Московский технический университет связи и информатики»
Практическая работа
Направление подготовки 10.05.02
«Информационная безопасность телекоммуникационных систем»
______________________________________
Руководитель
Кудряшов Всеволод Владимирович
______________________________________
Работа соответствует предъявляемым
требованиям
Заведующий кафедрой.
___________________________.
«_______» _____________________ 20____г.
Москва – 202 г.
Практическая работа №13
Устанавливаем GNU Octave
Для решения задачи необходимо составить уравнение: функции высоты камня над поверхностью Земли от времени.
Сила тяжести действующая на камень
Fт=mg
Где m - масса камня
g - ускорение свободного падения
Т.к на камень не действуют второстепенные силы, и он движется вертикально вниз, то движение происходит по 1 координате (пусть у - прямая относительно которой движется тело)
Ускорение направлено противоположенно т.е.
ma= -F
F= Fт /// => ma = - mg => g= -a (масса сокращается , значит не учитывается в данном уравнении движения)
Используем формулу равноускоренного движения vn=Vn0 + gt (n0=0 ), тогда получим приращение скорости
А через производную dv/dt= -g >>
V(t)= -g(t), а относительно прошлого приращения
V(t)=V(t-1) -g(t)
Теперь занесем полученные данные в GNU Octave
3.1 Указываем функцию
Через производную - скорость
И через нее же - ускорение
function dydt = F(y, t)
g = 10
dydt(1) = y(2)
dydt(2) = -g
Endfunction
3.2 указываем диапазон времени
>> t0 = 0
t0 = 0
>> tend = 5
tend = 5
>> deltaT = 0.1
deltaT = 0.1000
3.3 указываем границы функции, высоту
>> y0 = [100; 0]
3.4 запускаем уравнение с введенными данными
y = lsode("F", y0, t)
3.5 Строим график относительно т
plot(y(:,1), t )
3.6 Получим параболу, которая соответствует уравнению h=at^2/2
График построен верно