Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,580

ЗАДАЧА ПРЕСЛЕДОВАНИЯ В СИСТЕМЕ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ «MATHCAD»

Дубанов А.А. 1 Эрдынеева Л.И. 1
1 Сибирский Государственный Универститет Телекоммуникаций и Информатики
Описывается модели решения задачи преследования. Преследование производится по пересеченной местности. Ландшафт местности представлен точечным базисом, по которому выполнена полиномиальная регрессия. По результатам математического моделирования созданы анимированные изображения.
задача преследования
траектория кролика
траектория лисы
точечный базис
1. «Задачи по физике», Библиотечка «Квант», выпуск № 5, 1980, И.Ш. Слободецкий, Л.Г. Асламазов.
2. URL: http://blagovest2002.narod.ru/files/Fox_and_Rabbit_Sinusoud.rar.
3. URL: https://www.youtube.com/watch?v=ExIC7n1mlKc , «Кролик» и «Лиса» на плоскости.
4. URL: http://www.youtube.com/watch?v=Xv_q_EgdUX4 Траектория «Кролика».
5. URL: https://www.youtube.com/watch?v=dsfkqSlZGk8 , «Кролик» и «Лиса».
6. URL: http://dubanov.exponenta.ru/2014/Paper/Fox_and_Rabbit_2015.html Сайт «http://dubanov.exponenta.ru» «Геометрическое моделирование в системе MathCAD».
7. URL: http://dubanov.exponenta.ru/2014/Paper/Archive/Fox_and_Rabbit.zip Сайт «http://dubanov.exponenta.ru» «Геометрическое моделирование в системе MathCAD».
8. URL: http://www.youtube.com/watch?v=_iQAFGB6iFM «Лиса» работает на опережение. Проекция на горизонтальную плоскость.
9. URL: http://www.youtube.com/watch?v=cerQXgDUwwM «Лиса» работает на опережение. Трехмерный вид.

В данной статье предлагаются методы решения классической задачи из школьного курса физики, описанной в замечательной книге «Задачи по физике», Библиотечка «Квант», выпуск № 5, 1980, И.Ш. Слободецкий, Л.Г. Асламазов [1]. Позволю себе привести полный текст задачи № 13: «За лисой, бегущей равномерно и прямолинейно со скоростью v1 гонится собака, скорость которой v2 постоянна по абсолютной величине и направлена все время на лису. В момент времени, когда скорости v1 и v2 оказались взаимно перпендикулярными, расстояние между лисой и собакой было равно l. Каково было ускорение собаки в этот момент?»

Мы в системе MathCAD будем решать похожие задачи, имеющие свое происхождение из указанной выше классической задачи. «Кролик» движется по случайной траектории. Его преследует «Лиса». Необходимо найти траекторию движения «Лисы», ее зависимость от времени. Вектор скорости «Лисы» направлен на «Кролика» и по модулю неизменен.

«Кролик» и «Лиса» движутся по плоскости

Пусть траектория движения «кролика» (рис. 1) описывается уравнением

dub03.wmf,

а траектория движения «лисы»:

dub04.wmf, dub05.wmf – абсолютная неизменная величина скорости движения «лисы».

duban1.wmf

Рис. 1. «Кролик» движется по плоскости

Тот факт, что вектор скорости движения «лисы» направлен все время на кролика, и то, что абсолютная величина скорости движения «лисы» равна dub06.wmf, описывается системой уравнений (1.1):

dub07.wmf (1.1)

где dub08.wmf – координаты траектории «кролика» на плоскости,

dub09.wmf – координаты траектории «лисы», dub10.wmf – абсолютная величина скорости «лисы», t – время.

Система уравнений (1.1) относительно переменных dub11.wmf, dub12.wmf имеет решение (1.2):

dub13.wmf (1.2)

Если бы траектория движения «кролика»

dub14.wmf

имела бы явную функциональную зависимость от времени, то данная задача в таком математическом пакете, как «MathCAD 15», реализуется в несколько операторов, подводящих ко встроенным решателям систем дифференциальных уравнений I порядка (rkfixed, Rkadaptи т.д.). По следующей ссылке [2] вы сможете скачать архивированный файл программы, выполненной в системе «MathCAD 15», по следующей ссылке вы сможете посмотреть анимированное изображение результатов работы данной программы [3].

«Кролик» и «Лиса» движутся по пересеченной местности

Произведем постановку задачи преследования «Лисой» «Кролика», двигающихся по пересеченной местности, заданной поверхностью Z = f(X, Y).Траектория движения «Кролика» задается уравнением

dub15.wmf,

где T – время движения «Кролика».

Необходимо найти траекторию движения «Лисы»

dub16.wmf,

при следующих условиях:

1. Горизонтальная проекция скорости движения «Лисы» на плоскость (X, Y) направлена на горизонтальную проекцию точек движения «Кролика» (рис. 1).

2. Скорости движения «Кролика» и «Лисы» неизменны и равны величинам dub17.wmf, dub18.wmf.

Условия для решения поставленной задачи приводят к системе уравнений (2.1)

dub19.wmf. (2.1)

Первое уравнение системы (2.1) означает то, что горизонтальная проекция вектора скорости «Лисы» dub20.wmf на плоскость (X, Y) направлена на горизонтальную проекцию dub21.wmf текущего положения «Кролика». Второе уравнение означает то, что скорость передвижения «Лисы» равна постоянной величине dub22.wmf. Третье уравнение получается из дифференцирования уравнения Z = f(X, Y) по времени T. Фактически, это означает то, что точка траектории «Лисы»

dub23.wmf

принадлежит поверхности Z = f(X, Y). Cистема уравнений (2.1) имеет решение в явном виде относительно переменных

dub24.wmf, dub25.wmf, dub26.wmf:

dub27.wmf (2.2)

При реализации задачи, рассматриваемой в данной статье, в математических пакетах, в частности в системе «MathCAD», мы предприняли следующее:

1. В системе «AutoCAD» была построена модель поверхности Z = f(X, Y) и был произведен импорт по горизонталям в текстовые файлы. Затем точки были обработаны в системе «MathCAD» и была произведена полиномиальная регрессия. На плоскости (X, Y) была построена дополнительная сетка для проведения двумерной сплайн-интерполяции и построения поверхностей dub28.wmf и dub29.wmf, которые используются при решении системы уравнений (2.2).

2. В системе «AutoCAD» формируется массив горизонтальной проекции точек траектории «Кролика» dub30.wmf. Далее, вводится формальный параметр t, чтобы мы могли построить параметрическую функцию

dub31.wmf.

Нами, в качестве формального параметра t, был выбран такой, который равномерно проходит через значения [0, … N – 1], где N – количество точек массива dub32.wmf. После выполнения сплайн-интерполяции функция траектории движения «Кролика» принимает вид

dub33.wmf.

Далее, производится параметризация формального параметра от длины дуги траектории движения «Кролика»:

dub34.wmf.

Поскольку, длина дуги траектории «Кролика» может быть выражена: dub35.wmf, где T – это время движения «Кролика», то мы теперь можем в качестве основного параметра использовать время T.

3. Анимированное изображение движения «Кролика» нами предоставлено здесь [4].

4. При заданных стартовых позициях «Кролика» и «Лисы», при заданных значениях dub36.wmf и dub37.wmf было получено численное решение системы дифференциальных уравнений (2.2). С результатом анимации можно ознакомиться здесь [5]. На анимации видно, что при выбранных значениях «Лиса» догоняет «Кролика» несколько раз.

5. Полный листинг программы, выполненной в системе MathCAD с подробными комментариями можно посмотреть здесь [6].

6. Полный архив программы можно скачать здесь [7].

«Лиса» работает на опережение

Для прогнозирования «Лисой» движения «Кролика», мы предлагаем использовать следующую модель (рис. 2).

dub38.wmf, (3.1)

где dub39.wmf – горизонтальная проекция вектора dub40.wmf (рис. 2). Здесь векторная величина скорости «Кролика» равна

dub41.wmf.

duban2.wmf

Рис. 2. «Лиса» прогнозирует движение «Кролика»

На выборе параметра приращения по времени ∆T остановимся подробнее. Введем следующую весовую функцию:

dub42.wmf,

где

dub43.wmf,

dub44.wmf – выбираемое нами пороговое значение дистанции между «Кроликом» и «Лисой». Измененный листинг программы в системе MathCAD вы можете просмотреть и скачать на сайте [7]. Анимированное изображение ситуации, когда «Лиса» бежит на опережение, вы можете посмотреть здесь (в проекции на плоскость [8] и трехмерный вид [9]).


Библиографическая ссылка

Дубанов А.А., Эрдынеева Л.И. ЗАДАЧА ПРЕСЛЕДОВАНИЯ В СИСТЕМЕ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ «MATHCAD» // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 9-1. – С. 7-11;
URL: http://www.applied-research.ru/ru/article/view?id=10181 (дата обращения: 28.02.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074