Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кудряшов Всеволод / 13 лаб / 13 практика

.docx
Скачиваний:
2
Добавлен:
04.12.2024
Размер:
60.58 Кб
Скачать

Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования «Московский технический университет связи и информатики»

Практическая работа

Направление подготовки 10.05.02

«Информационная безопасность телекоммуникационных систем»

______________________________________

Руководитель

Кудряшов Всеволод Владимирович

______________________________________

Работа соответствует предъявляемым

требованиям

Заведующий кафедрой.

___________________________.

«_______» _____________________ 20____г.

Москва – 202 г.

Практическая работа №13

  1. Устанавливаем GNU Octave

  1. Для решения задачи необходимо составить уравнение: функции высоты камня над поверхностью Земли от времени.

Сила тяжести действующая на камень

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)

  1. Теперь занесем полученные данные в 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

График построен верно