КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ ПРОЦЕССОВ

СРЕДСТВАМИ СИСТЕМЫ КОМПЬЮТЕРНОЙ МАТЕМАТИКИ MAPLE

Голанова Анна Викторовна

канд. пед. наук, доц кафедры информатики и вычислительной математики

Ленинградского государственного университета имени А.С. Пушкина,

РФ, г. Пушкин

E-mail: golanova@yandex.ru

Голикова Екатерина Ивановна

канд. пед. наук, доц. кафедры информатики и вычислительной математики

Ленинградского государственного университета имени А.С. Пушкина,

РФ, г. Пушкин

E-mail: golikova_kat@inbox.ru

COMPUTER SIMULATION OF PHYSICAL PROCESSES

BY MEANS OF COMPUTER ALGEBRA SYSTEM MAPLE

Anna Golanova

candidate of Science, associate professor of department of computer science and calculus mathematics of Pushkin Leningrad State University,

Russia, Pushkin

Ekaterina Golikova

candidate of Science, associate professor of department of computer science and calculus mathematics of Pushkin Leningrad State University,

Russia, Pushkin

АННОТАЦИЯ

В статье описан и реализован в системе компьютерной математики Maple алгоритм решения физических задач, приводящих к дифференциальным уравнениям первого порядка.

ABSTRACT

The article described and implemented in the computer algebra system Maple algorithm for solving physical problems leading to differential equations of the first order.

Ключевые слова: физическая задача, дифференциальное уравнение, компьютерное моделирование, математическая модель.

Keywords: physical problem, differential equation, computer simulation, mathematical model.

Компьютерное моделирование используется для описания и анализа процессов разнообразной природы и является частью научно-исследовательской деятельности. Построение компьютерной модели требует не только знаний в конкретной предметной области, но и знаний в области вычислительной математики и программирования.

Бакалавры, обучающиеся по направлениям 09.03.03 «Прикладная информатика» (профиль подготовки «Прикладная информатика в экономике») и 44.03.05 «Педагогическое образование» (профили подготовки «Математика и информатика», «Информатика»), в рамках дисциплины «Компьютерное моделирование» должны иметь представление о классах программных средств и их назначении, знать инструментарий и владеть технологическими приёмами работы. При изучении этой дисциплины рассматриваются различные физические и геометрические задачи. В процессе решения таких задач возникает необходимость построения их математической и компьютерной моделей. Математическими моделями таких задач очень часто являются обыкновенные дифференциальные уравнения первого порядка с заданными начальными условиями [1; 2].

Процесс решения геометрических и физических задач состоит из последовательности взаимосвязанных этапов. Для геометрической задачи этот процесс был подробно описан в [2]. В данной статье мы остановимся на описании алгоритма решения физической задачи.

Для решения физической задачи, сводящейся к обыкновенному дифференциальному уравнению первого порядка, необходимо:

1)  построить чертёж в декартовой системе координата;

2)  записать второй закон Ньютона в векторной форме и найти проекции на оси координат;

3)  получить дифференциальное уравнение и найти его общее решение;

4)  найти частное решение дифференциального уравнения, используя исходные данные задачи;

5)  используя найденное частное решение, найти искомую величину и ответить на поставленный вопрос задачи.

Для нахождения, численного и аналитического решений полученного дифференциального уравнения удобно использовать онлайн калькуляторы (http://matematikam.ru/calculate-online/differential-equations.php, http://reshit-online.ru/diff-equation.html, http://math.semestr.ru/math/diffur.php и др.) и различные системы компьютерной математики (Mathematica, Maple, MatLab, Mathcad, Maxima и др.). Основная задача онлайн калькуляторов – нахождение аналитического решения. Использование систем компьютерной математики направлено, в первую очередь, на нахождение не только аналитического, но численного решения обыкновенного дифференциального уравнения, а также на построение графиков решений.

Для наглядной иллюстрации вышеописанного алгоритма, рассмотрим решение следующей физической задачи.

Снаряд выпущен вертикально вверх с начальной скоростью v₀ = 100 м/с. Сопротивление воздуха замедляет его движение, сообщая снаряду отрицательное ускорение, равное –kv² (где v – мгновенная скорость снаряда, а k – аэродинамический коэффициент). Определить время достижения снарядом наивысшего положения (по [3, с. 111, № 36]).

Для решения данной задачи:

1.  Построим чертёж к задаче в декартовых координатах (Рис. 1).

Рисунок 1. Иллюстрация задачи

2.  Запишем второй закон Ньютона в векторной форме:

где: m – масса снаряда, F₁ – сила притяжения, F₂ – сила аэродинамического сопротивления, F₃ – сила ветра.

Найдём проекции на оси:

где: g – ускорение свободного падения, V_x – скорость ветра (проекция на ось X).

3.  Проинтегрируем второе уравнение системы и подставим результат в первое уравнение:

В результате получим уравнение первого порядка с разделяющимися переменными. Его общее решение:

4.  Для нахождения частного решения воспользуемся начальным условием v_y(0) = 100:

Выразим постоянную C:

Таким образом, частное решение имеет вид:

5.  В момент достижения снарядом наивысшего положения его скорость равна 0, следовательно:

Зададим значения коэффициентов:

·     аэродинамический коэффициент k = 0,5 (коэффициент лобового сопротивления для конуса 2:1);

·     ускорение свободного падения g = 9,8 м/с²;

·     скорость бокового ветра V_x = 2 м/с.

В результате получим: t_max = 0,6267 c.

Таким образом, время достижения снарядом наивысшего положения составляет 0,6267 c.

Для нахождения, аналитического и численного решений полученного дифференциального уравнения, будем использовать систему компьютерной математики Maple.

1.  Найдём аналитическое решение дифференциального уравнения

2.  Найдём численное решение дифференциального уравнения. Численное решение ищем, используя метод Рунге-Кутта четвертого порядка и метод Эйлера.

3.  Построим графики найденных решений.

Рисунок 2. Графики аналитического и численного решений дифференциального уравнения

Отметим, что использование системы компьютерной математики Maple позволило найти аналитическое и численное решения полученного дифференциального уравнения и построить их графики в одной системе координат, визуально оценить различия аналитического и численного решений.

Оценим погрешность найденного численного решения.

Таблица 1.

Абсолютная и относительная погрешности численных решений

Полученные результаты позволяют сделать вывод о том, что при выбранном количестве разбиений (n = 1000), аналитическое и численное решение отличаются весьма незначительно.

Для более тщательного и подробного исследования полученной компьютерной модели, необходимо изменение исходных данных и сравнение полученных результатов. Например, при постоянном значении аэродинамического коэффициента отследить время достижения снарядом наивысшего положения при изменении начальной скорости или при постоянном значении начальной скорости отследить время достижения снарядом наивысшего положения при изменении аэродинамического коэффициента.

Список литературы:

1.    Голанова А.В., Голикова Е.И. К вопросу об отборе содержания лабораторных работ по дисциплине «Компьютерное моделирование» для бакалавров по направлению «Педагогическое образование». // XVIII Царскосельские чтения: материалы междунар. науч. конф. – СПб: ЛГУ им. А.С. Пушкина, 2014. – Т. III. – С. 130–134.

2.    Голанова А.В., Голикова Е.И. Применение системы компьютерной математики Maple для решения задач дифференциальной геометрии. // Естественные и математические науки в современном мире: Сб. ст. по материалам XXIII междунар. науч. конф. № 23. – Новосибирск: Изд. «Сибак», 2014. – С. 23–29. 

3.    Самойленко А.М., Кривошея С.А., Перестюк Н.А. Дифференциальные уравнения: примеры и задачи. Учеб. пособие. – 2-е изд., перераб. – М.: Высш. шк., 1989. – 383 с.