close

Вход

Забыли?

вход по аккаунту

?

Численные решения системы линейных уравнений с дробной степенью матрицы дифференциального оператора.

код для вставкиСкачать
Вестник КРАУНЦ. Физ.-мат. науки. 2013. № 2 (7). C. 7-11. ISSN 2079-6641
МАТЕМАТИКА
УДК 519.622+519.633
ЧИСЛЕННЫЕ РЕШЕНИЯ СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ
С ДРОБНОЙ СТЕПЕНЬЮ МАТРИЦЫ ДИФФЕРЕНЦИАЛЬНОГО
ОПЕРАТОРА ∗
И.А. Ильин1, 2 , Д.С. Нощенко2 , А.С. Пережогин1, 2
1
Институт космофизических исследований и распространения радиоволн ДВО РАН,
684034, Камчатский край, c. Паратунка, ул. Мирная, 7
2 Камчатский государственный университет имени Витуса Беринга, 683031,
г. Петропавловск-Камчатский, ул.Пограничная д.4
E-mail: d72156@gmail.ru
В работе исследуются численные решения линейных дифференциальных уравнений и уравнений в частных производных с дробной степенью разностного оператора. Показаны особенности решения задачи Коши для обыкновенных дифференциальных уравнений с дробной степенью разностного оператора. Введено
условие регуляризации для разностого аналога дробной степени оператора.
Ключевые слова: численный метод, оператор дробной степени, дифференциальные уравнения.
© Ильин И.А., Нощенко Д.С., Пережогин А.С., 2013
MATHEMATICA
MSC 65L12
NUMERICAL SOLUTION OF A LINEAR SYSTEM WITH A FRACTIONAL
POWER
I.A. Il’in1, 2 , D.S. Noshchenko2 , A.S. Perezhogin1, 2
1
Institute of Cosmophysical Researches and Radio Wave Propagation Far-Eastern Branch,
Russian Academy of Sciences, 684034, Kamchatskiy Kray, Paratunka, Mirnaya st., 7,
Russia
2 Kamchatka state university by Vitus Bering, 683031, Petropavlovsk-Kamchatsky,
Pogranichnaya st., 4, Russia
E-mail: d72156@gmail.ru
We investigate numerical solutions of linear ordinary and partial differential equations.
Cauchy’s problem for ordinary equations of first and second order are generalized
with fractional power of finite operator. Condition of regularization for the difference
operator analogue of a fractional power is represented.
Key words: fractional calculus, numerical method
© Il’in I.A., Noshchenko D.S., Perezhogin A.S., 2013
∗ Работа
выполнена при поддержке Министерства образования и науки РФ в рамках программы
стратегического развития ФГБОУ ВПО «Камчатский государственный университет имени Витуса
Беринга» на 2012-2016 гг.
7
ISSN 2079-6641
Ильин И.А., Нощенко Д.С., Пережогин А.С.
Введение
Математическое моделирование физических процессов выполняется с помощью
дробного исчисления в тех случаях, когда физическая система обладает «памятью»
[1]. Это означает, что значения искомой функции зависят от значений функции в
предыдущие или следующие моменты времени. В связи с этим классические модели
физических систем можно обобщить с помощью введения в рассмотрение дробной
степени дифференциальных операторов, которые приставлют собой свертку со степенний функцией. Предполагается, что подобное математическое описание наиболее
точно описывает физические процессы, которые происходят во фрактальных средах.
В настоящей работе рассматриваются простейшие примеры обыкновенных дифференциальных уравнений, для которых вычислительная схема имеет в качестве
оператора дробную степень матрицы. Коэффициенты матрицы согласуются с формальным обобщением разностных оператор путем замены количества сочетаний на
Гамма-функцию.
В одном из примеров мы покажем вычислтельную схемму для дробного диффузионного уравнения. Параметр дробности вводится в систему через временнную
переменную.
Задача I
Рассмотрим уравнение
d α x(t)
+ x(t) = 0,
dt α
α ∈ (0, 1]
x(0) = 1
(1)
(2)
В этой задаче мы можем использовать численную схему, заменив дифференциальный оператор на матрицу. В случае α = 1, мы получаем обыкновенное дифференциальное уравнение первого порядка:
dx xi+1 − xi
≈
, i = 0..n, x0 = 1
dt
τ
Система линейных уравнений для уравнения (1) с учетом начальных условий
будет иметь вид:

 
 1
 
1
0 .. 0
x0
x0
− τ 0 .. 0
−1 1 .. 0 x1 
 1 0 .. 0 x1 
 

 
1
 0 −1 .. 0 x2  = −  0 1 .. 0 x2 
(3)
 

 
τ
 ..
 .. .. .. ..  .. 
.. .. ..  .. 
0
0 .. 1
xn
xn
0 0 .. 0
Вычислим дробную степень α матрицы А в левой части уравнения, которая аппроксимирует разность первого порядка. В нашем случае параметр α выберем равным 12 . Используем рекурсивную схему (4) для нахождения дробной степени матрицы.
Bn+1 =
Bn + A · B−1
n
2
8
(4)
Численные решения системы линейных уравнений . . .
√
Если B0 = A и n → ∞, мы получаем Bn ≈ A. Таким
системе уравнений (5) для численного моделирования.

 
 1
− √τ
1
0
0 0
x0





1
−.5
1
0 0 x1 
1
√ 
= −





−.125
−.5
1
0
x
τ
0
2
x3
−.0625 −.125 −.5 1
0
ISSN 2079-6641
образом, мы приходим к
0
0
1
0
0
0
0
1
 
0
x0


0 x1 

0 x2 
x3
0
(5)
Когда α = 21 , решение имеет особенность в начальной точке t = 0 (рис. 1 a). Модифицируем систему линейных уравнений (5). Для этого выполним замену вектора
(x0 , x1 , . . . , xn ) на вектор (x0 , x1 − x0 , . . . , xn − x0 ) в левой части системы. Таким образом
решение системы уравнений не имеет особенность в начальной точке (рис. 1 б).
Рис. 1. a) решения (3) и (5); 1 – α = 1, 2 – α = 12 ; b) регуляризованное решение (5)
Рассмотрим коэффициенты матрицы
разности порядка ν:
1
H f (x) = lim ν
τ→0 τ
ν
√
A. Мы рассмотрим приближение дробной
∞
Γ(ν + 1)
∑ k!Γ(ν − k + 1) f (x − kτ),
(6)
k=0
где Γ(x) – гамма-функция.
Вычисление
коэффициентов при f (x−kτ) дает совпадение со значением элементов
√
матрицы A. Таким образом, мы получаем совпадение вычислительного метода и
разностным оператором Грюнвальда-Летникова.
Задача II
Рассмотрим формальную запись уравнения с дробной степенью:
dα x
+ x(t) = 0,
dt α
α ∈ (1, 2]
x0 (0) = 0
x(0) = 1,
(7)
(8)
Представим аналог дифференциальног оператора в левой части как √
произведение
матрицы А самой на себе, а в другом случае – как произведение А на A. Матрица
A является первой разностью в системе (1).
9
ISSN 2079-6641
Ильин И.А., Нощенко Д.С., Пережогин А.С.
Для дробной степени α = 32 система имеет решение в виде, аналогичном для
осциллятора с трением. Однако в этом случае отсутствует характерный масштаб
затухания. Это основная особенность для динамических систем с дробными операторами. Приведенные численные решения совпадают с решениями для дробного
осциллятором [2].
Задача III
Пусть u(t, x, y) задана в области (x, y) ∈ [0, 1] × [0, 1] и t ∈ [0, ∞). Функция удовлетворяет уравнению (9) с начальными и граничными условиями.
∂ u ∂ 2u ∂ 2u
= 2+ 2
∂t
∂x
∂y
u(0, x, y) = sin(πx) sin(πy),
u(t, x, y)|Γ = 0
(9)
(10)
Решение этой системы представляется в виде:
2
u(t, x, y) = e−2π t sin(πx) sin(πy)
(11)
Численное решение можно представить явной схемой:
k
uk+1
i j − ui j
(uki+1 j + uki−1 j + ukij+1 + ukij−1 − 4ukij )
(12)
τ
h2
Левая часть (12) может быть представлена в виде матрицы системы линейных
уравнений. Дробная степень этой матрицы дает новое приближение для разностного
аналога дифференциального оператора.
Рассмотрим классическое решение (11) и решение с дробной степенью матрицы в
левой части (12). Отличия решений в случаях α = 1 и α = 21 представлено на рис. 2.
=
Рис. 2. Решение системы (5) при значениях: 1 – α = 2, 2 – α = 32
В численном методе, когда в качестве приближения по переменной t выбирается
первая разность, то есть обычная матрица A, мы получаем приближение классического решения (11). В противном случае решение замедляет свое спадание со времением
для дробной степени матричного оператора (рис. 3). Таким образом динамика такой
задачи совпадает с решением для примера системы в задаче I.
10
Численные решения системы линейных уравнений . . .
ISSN 2079-6641
Рис. 3. Численное решение. Классическая схема – сплошная линия; дробная степень
разностного оператора – сплошная линия с крестами
Заключение
Численные решения обыкновенных дифференциальных уравнений и уравнений
в частных производных с дробной степенью оператора показывает наличие «эффекта памяти». Для физических систем это означает, что процесс имеет длительную
релаксацию процесса. Такую методику моделирования можно применять в задачах
расчета напряженно-деформированного состояния сред.
Библиографический список
1. Тарасов В.Е. Модели теоретической физики с интегро-дифференцированием дробного порядка.
М.: Ижевск: РХД, 2011. 568 с.
2. Нахушев А.М. Дробное исчисление и его применение. М.: Физматлит, 2003. 272 c.
Поступила в редакцию / Original article submitted: 05.11.2013
11
Документ
Категория
Без категории
Просмотров
7
Размер файла
296 Кб
Теги
численные, решение, уравнения, дифференциальной, степенью, система, матрица, оператора, линейный, дробной
1/--страниц
Пожаловаться на содержимое документа