ФУНКЦИИ VBA EXCEL ДЛЯ ЧИСЛЕННОГО РЕШЕНИЯ ЗАДАЧИ КОШИ СИСТЕМЫ ОДУ 1-ГО ПОРЯДКА | Опубликовать статью ВАК, elibrary (НЭБ)

ФУНКЦИИ VBA EXCEL ДЛЯ ЧИСЛЕННОГО РЕШЕНИЯ ЗАДАЧИ КОШИ
СИСТЕМЫ ОДУ 1-ГО ПОРЯДКА

Научная статья

Сдвижков О.А.*

Российский государственный университет туризма и сервиса, Пушкино МО, Россия

* Корреспондирующий автор (oasdv[at]yandex.ru)

Аннотация

Приведены программные коды пользовательских функций VBA Excel, реализующих численные методы решения задачи Коши для системы обыкновенных дифференциальных уравнений 1-го порядка методами Эйлера и Рунге-Кутта четвертого порядка точности. Параметрами функций являются начальное и конечное значения независимой переменной, массив начальных значений зависимых переменных и число шагов, в программных кодах применяются динамические массивы. Пользовательские функции можно применять к задаче Коши и одного дифференциальных уравнений 1-го порядка. Все пользовательские функции имеют описания, поясняющие их назначения. Применение пользовательских функций показано на конкретных примерах численного решения задачи Коши системы дифференциальных уравнений 1-го порядка и отдельных дифференциальных уравнений 1-го порядка.

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

EXCEL VBA FUNCTIONS FOR NUMERICAL SOLUTION OF A CAUCHY PROBLEM
OF THE FIRST-ORDER ODE SYSTEM

Research article

Sdvizhkov O.A.*

Russian State University of Tourism and Service, Pushkino, Moscow Oblast, Russia

* Corresponding author (oasdv[at]yandex.ru)

Abstract

The current study presents program codes of custom VBA Excel functions that implement numerical methods for solving the Cauchy problem for a system of 1st order ordinary differential equations using the Euler and Runge-Kutta 4th order methods. The parameters of the functions are the initial and final values of the independent variable, an array of initial values of the dependent variables, and the number of steps. Program codes use dynamic arrays. Custom functions can be applied to the Cauchy problem and one of the 1st order differential equations.  All of the custom functions have descriptions that explain their purpose. The article demonstrates the use of custom functions using specific examples of numerical solutions to the Cauchy problem of 1st order differential equations system and individual 1st-order differential equations.

Keywords: system of differential equations, Cauchy problem, numerical methods, code, dialog box.

Введение

Системы обыкновенных дифференциальных уравнений 1-го порядка имеют большое прикладное значение [9]. Однако численные методы решения задачи Коши системы обыкновенных дифференциальных уравнений 1-го порядка отличаются большой трудоемкостью [1], [2], [6]. Поэтому задача информатизации этих методов, доступной самому широкому кругу пользователей ПК, каким является круг пользователей Microsoft Office, несомненно, является актуальной.

Цель данного исследования – разработка на языке программирования VBA Excel [3], [4], [7] программных кодов пользовательских функций [10], которыми численно решается задача Коши системы обыкновенных дифференциальных уравнений 1-го порядка методами Эйлера и Рунге-Кутта.

Программный код функции VBA Excel для решения методом Рунге-Кутта системы двух дифференциальных уравнений приведен в [5]. Пользовательские функции VBA Excel для численного решения задачи Коши обыкновенного дифференциального уравнения 1-го порядка приведены в [9]. 

1. Пользовательская функция ЭЙЛЕР

В методе Эйлера численного решения задачи Коши системы дифференциальных уравнений

    (1.1)

     (1.2)

на отрезке [x₀, x₀ + n·h] применяется формула:

     (1.3)

j = 1, 2, …, m; i= 1, 2, …, n.

Правые части системы (1.1) в VBA Excel можно задать пользовательской функцией СИСТЕМ, листинг которой имеет вид:

Function СИСТЕМ (Массив As Variant)

‘массив состоит из значений:i – номера уравнения,х,y1,y2,…, yn

Select Case Массив(1)

Case 1

СИСТЕМ = f₁(Массив(2), Массив(3), …, Массив(m))

Case 2

СИСТЕМ = f₂(Массив(2), Массив(3), …, Массив(m))

…………………

Case m

СИСТЕМ = f_n(Массив(2), Массив(3), …, Массив(m))

End Select

End Function

Пусть система (1.1) задана пользовательской функцией СИСТЕМ. Тогда пользовательская функция ЭЙЛЕР, код VBA Excel которой приведен в листинге 1, выполняет вычисления по формуле (1.3).

Листинг 1. Код функции ЭЙЛЕР и ее описания

Function ЭЙЛЕР(Нач_x0, Конеч_x, Нач_y0(), Шагов) As Variant

Dim n As Integer, m As Integer, x(), y(), z(), c()

n = Шагов:m = UBound(Нач_y0)

ReDim z(1 To m, 0 To n):ReDim x(0 To n)

ReDim y(1 To m + 2):ReDim c(1 To m)

x(0) = Нач_x0:h = (Конеч_x – x(0)) / n

For i = 1 To m

z(i, 0) = Нач_y0(i)

Next

For j = 1 To n

 For i = 1 To m

y(1) = i: y(2) = x(0) + h * (j – 1)

 For r = 3 To m + 2

y(r) = z(r – 2, j – 1)

 Next

 Массив = y

z(i, j) = z(i, j – 1) + h * СИСТЕМ(Массив)

 Next

Next

For i = 1 To m

c(i) = z(i, n)

Next

ЭЙЛЕР = c

End Function

Sub InstallFunc2()

Application.MacroOptions Macro:=”ЭЙЛЕР”, Description:= _

“Возвращает значения yi, полученные методом Эйлера”

End Sub

В листинг 1 входят массивы:

z(1 To m, 0 To n)– двумерный массив значений функций у₁, у₂, …, y_m в точках x₀, x₀+h, …, x₀+n·h;

x(0 To n) – массив значений x₀, x₀+h, …, x₀+n·h;

y(1 To m + 2) – массив значений (i, x, y₁, y₂, … y_m), где i – номер уравнения системы;

c(1 To m) – массив итоговых значений массива z.

Пример 1.1. Найти приближенные значения решения системы

с начальными условиями . Вычислить значения решения при х=0,1; 0,2; 0,3; 0,4.

1. В стандартный модуль вводится листинг функции СИСТЕМ:

Function СИСТЕМ(Массив As Variant)

Select Case Массив (1)

Case 1

СИСТЕМ = Массив (4)

Case 2

СИСТЕМ = Массив (3)

End Select

End Function

2. В диапазон А1:А5 вводятся значения, принимаемые переменной х: 0; 0,1; 0,2; 0,3; 0,4.
3. В диапазон В1:С1 вводятся начальные значения, принимаемые у₁ и у₂: 0 и 1.
4. Выделяется диапазон В2:С2, вызывается функция ЭЙЛЕР и вводятся данные:

Рис. 1.1 – Применение функции ЭЙЛЕР в примере 1.1

 

5. Команда Ctrl+Shift +Enter возвращает в диапазоне В2:С2 результаты.
6. Копирование данных диапазона В2:С2 в диапазон В3:С5 завершает решение задачи (рис. 1.2).

Рис. 1.2 – Вид решения примера 1.1 функцией ЭЙЛЕР

 

Пример 1.2. Найти приближенные значения решения системы

с начальными условиями . Вычислить значения решения при х=0,1; 0,2; 0,3; 0,4; 0,5.

1. В стандартный модуль вводится листинг функции СИСТЕМ:

Function СИСТЕМ(Массив As Variant)

Select Case Массив(1)

Case 1

СИСТЕМ = Массив(4) + Массив(5)

Case 2

СИСТЕМ = Массив(3) + Массив(5)

Case 3

СИСТЕМ = Массив(3) + Массив(4)

End Select

End Function

2. В диапазон А1:А6 вводятся значения, принимаемые х: 0; 0,1; 0,2; 0,3; 0,4; 0,5.
3. В диапазон В1:D1 вводятся начальные значения, принимаемые у₁, у₂, у₃: 1, 2, 3.
4. Выделяется диапазон В2:D2, вызывается функция ЭЙЛЕР и вводятся данные:

Рис. 1.3 – Применение функции ЭЙЛЕР в примере 1.2

 

5. Команда Ctrl+Shift +Enter возвращает в диапазоне В2:D2 результаты.
6. Копирование данных диапазона В2:D2 в диапазон В3:D6 завершает решение задачи (рис. 1.4).

Рис. 1.4 – Вид решения примера 1.2 функцией ЭЙЛЕР

 

Функция Эйлер применима и к задаче Коши обыкновенного дифференциального уравнения 1-го порядка.

Пример 1.3. Найти приближенные значения решения задачи Коши

    (1.4)

на отрезке [0, 1], разбивая отрезок на 10 равных частей.

В этом случае код функции СИСТЕМ имеет вид:

Function СИСТЕМ(Массив As Variant)

СИСТЕМ = Массив(2) + Массив(3)

End Function

Применение функции ЭЙЛЕР приводит к результатам, показанным на рисунке 1.3.

Рис. 1.5 – Вид решения примера 1.3 функцией ЭЙЛЕР

 

2. Пользовательская функция РУНГЕ

В методе Рунге-Кутта четвертого порядка точности, которым численно решается задача Коши (1.1, 1.2):

   (2.1)

где

    (2.2)

Пусть система (1.1) задана пользовательской функцией СИСТЕМ. Тогда пользовательская функция РУНГЕ, программный код VBA Excel которой приведен в листинге 2, выполняет вычисления по формулам (2.1), (2.2).

Листинг 2. Код функции РУНГЕ и ее описания

Function РУНГЕ(Нач_x, Конеч_x, Нач_y() As Variant, Шагов) _

As Variant

Dim n As Integer, m As Integer, h As Single, x(), y(), z(), c()

n = Шагов:m = UBound(Нач_y)

ReDim z(1 To m, 0 To n):ReDim x(0 To n)

ReDim y(1 To m + 2):ReDim c(1 To m)

Dim K(1 To 4) As Double

x(0) = Нач_x:h = (Конеч_x – x(0)) / n

For i = 1 To m

z(i, 0) = Нач_y(i)

Next

For i = 1 To n

For j = 1 To m

y(1) = j: y(2) = x(0) + h * (i – 1)

For r = 3 To m + 2

y(r) = z(r – 2, i – 1)

Next

Массив = y:K(1) = h * СИСТЕМ(Массив)

y(1) = j: y(2) = x(0) + h * (i – 1) + h / 2

For r = 3 To m + 2

y(r) = z(r – 2, i – 1) + K(1) / 2

Next

Массив = y

K(2) = h * СИСТЕМ(Массив)

y(1) = j: y(2) = x(0) + h * (i – 1) + h / 2

For r = 3 To m + 2

y(r) = z(r – 2, i – 1) + K(2) / 2

Next

Массив = y

K(3) = h * СИСТЕМ(Массив):y(1) = j: y(2) = x(0) + h * i

For r = 3 To m + 2

y(r) = z(r – 2, i – 1) + K(3)

Next

Массив = y

K(4) = h * СИСТЕМ(Массив)

z(j, i) = z(j, i – 1) + (K(1) + 2 * K(2) + 2 * K(3) + K(4)) / 6

Next

Next

For i = 1 To m

c(i) = z(i, n)

Next

РУНГЕ = c

End Function

Sub InstallFunc3()

Application.MacroOptions Macro:=”РУНГЕ”, Description:= _

“Возвращает значения yi, полученные методом Рунге-Кутта”

End Sub

В листинге 2 массивы аналогичные массивам листинга 1, только:

K(1 To 4) – массив значений, вычисляемых по формулам (2.2).

Решение примера 1.1 с помощью функции РУНГЕ показано на рисунке 2.1.

Рис. 2.1 – Вид решения примера 1.1 функцией РУНГЕ

Решение примера 1.2 с помощью функции РУНГЕ показано на рисунке 2.2.

Рис. 2.2 – Вид решения примера 1.2 функцией РУНГЕ

 

Решение примера 1.3 с помощью функции РУНГЕ показано на рисунке 2.3.

Рис. 2.3 – Вид решения примера 1.3 функцией РУНГЕ

 

Пример 2.1. Построить с помощью функции РУНГЕ графики численности популяций, описываемых системой Вольтерра – Лотка

при начальных условиях x₁(0)=3, x₂(0)=1; параметр t принимает значения от 0 до 3.

Расчетная таблица строится по аналогии с таблицей примера 1.1, значения независимой переменной в диапазоне А1:А31, в ячейку В2 вводится формула =РУНГЕ(0;A2;{3;1};A2/0,01). По таблице строятся графики численности популяций (рис. 2.4).

Рис. 2.4 – Графики численности популяций

Заключение

1. Программные коды пользовательских функций VBA Excel, реализующих численные методы решения задачи Коши для системы обыкновенных дифференциальных уравнений 1-го порядка методами Эйлера и Рунге-Кутта четвертого порядка, публикуются впервые.
2. Как показывают примеры, применение рассматриваемых функций не вызывает трудностей, избавляет от больших объемов рутинных вычислений, позволяет получать результаты по такому сложному разделу, как численные методы решения систем дифференциальных уравнений, на современном (цифровом) уровне.
3. Рассматриваемые пользовательские функции будут полезны как в научных исследования, связанных с системами обыкновенных дифференциальных уравнений 1-го порядка, так и в учебном процессе.

Конфликт интересов Не указан.

Conflict of Interest None declared.