CC BY
21
ЧЕБЫШЕВСКИЙ СБОРНИК
Том 23. Выпуск 5.
УДК 517 DOI 10.22405/2226-8383-2022-23-5-117-129
Метод Ритца решения дифференциальных уравнений в частных производных с применением теоретико-числовых сеток1
А. В. Родионов
Родионов Александр Валерьевич — старший преподаватель, Тульский государственный педагогический университет им. Л. Н. Толстого (г. Тула). e-mail: rodionovalexnandr@mail. ru
Аннотация
Рассмотрим задачу
Lu(x) = f( х), u(x )|aGs = д(х ),
где f(x), д(х) € E", L — линейный дифференциальный оператор с постоянными коэффициентами, Gs — единичный куб [0; 1]s.
Её решение сводится к отысканию минимума функционала
v(u(х)) = j ... JF (x,u, uXl,. .., uXs) dx\ . .. dxs
as
при заданных граничных условиях.
Значения функционала v(u(x)) в методе Ритца рассматриваются не на множестве всех допустимых функций и(х), а на линейных комбинациях
п
и(х) = Wo(х) wkWk(x),
к = 1
где Wfc (х) — некоторые базисные функции, которые будем находить с помощью теоретико-числовой интерполяции, причём Wo (ж) — функция, удовлетворяющая заданным граничным условиям, а остальные Wk (х) удовлетворяют однородным граничным условиям.
На этих полиномах данный функционал превращается в функцию ip(w) от коэффициентов w1,..., wn. Эти коэффициенты выбираются так, чтобы функция p(w) достигала экстремума. При некоторых ограничениях на функционал v(u(х)) и базисные функции Wfc (х) получим приближённое решение краевой задачи.
Ключевые слова: теоретико-числовой метод, дифференциальные уравнения в частных производных, вариационные методы.
Библиография: 10 названий.
1 Исследование выполнено при финансовой поддержке гранта Министерства образования и науки РФ на развитие молодежных лабораторий, в рамках реализации ТГПУ им. Л. Н. Толстого программы «Приоритет 2030» по Соглашению №073-03-2022-117/7 по теме «Теоретпко-числовые методы в приближенном анализе и их приложения в механике и физике»
Для цитирования:
А. В. Родионов. Метод Ритца решения дифференциальных уравнений в частных производных с применением теоретико-числовых сеток // Чебышевский сборник, 2022, т. 23, вып. 5, с. 117 129.
CHEBYSHEVSKII SBORNIK Vol. 23. No. 5.
UDC 517 DOI 10.22405/2226-8383-2022-23-5-117-129
The Ritz method for solving partial differential equations using
number-theoretic grids
A. V. Rodionov
Rodionov Alexander Valer'evich — senior lecturer, Tula State Lev Tolstoy Pedagogical
University (Tula).
e-mail: rodionovalexnandr@mail. ru
Consider the problem
Abstract
Lu(x) = f( x), u(x )|aGs = g(x ),
where f(x), g(x) G L is a linear differential operator with constant coefficients, Gs is the unit cube [0; 1]s.
Its solution is reduced to finding the minimum of the functional
v(u(x)) = j ... JF (x,u, uXl, .. ., uXs) dx\ ldotsdxs
Gs
under given boundary conditions.
The values of the functional v(u(x)) in the Ritz method are considered not on the set of all admissible functions u(x), but on linear combinations
n
u(x) = Wo (x) +^2wkWk(x),
k = l
where Wk(x) are some basic functions that we will find using number-theoretic interpolation, and W0 (x) is a function that satisfies the given boundary conditions, and the rest Wk (x) satisfy homogeneous boundary conditions.
On these polynomials, this functional turns into a function <f(W) of the coefficients wi,..., wn. These coefficients ^e ^osen so that the function <f(W) reaches an extremum. Under some restrictions on the functional v(u(x)) md the b^is functions Wk(x), we obtain an approximate solution of the boundary value problem.
Keywords: number-theoretic method, partial differential equations, variational methods.
Bibliography: 10 titles.
For citation:
A. V. Rodionov, 2022, "The Ritz method for solving partial differential equations using number-theoretic grids", Chebyshevskii sbornik, vol. 23, no. 5, pp. 117-129.
1. Введение
В середине прошлого века Н. М. Коробов ввёл в рассмотрение класс Е£, а > 1. Это класс функций ¡(х\,...,х3), которые имеют период равный единице по каждой из переменных и для коэффициентов Фурье этой функции С(т\,..., т3) выполняется оценка
С
|С(т\,... ,т3)\ ^ -—
(т\ ■ ... ■ т3)а'
где т = тах(|х|, 1).
Класс периодических функций Е" относительно нормы
|| ¡(х)\\Еа = \С (гп)(т1...т3)а | (1)
является несепарабельным банаховым пространством, изоморфным пространству 1о — пространству всех ограниченных последовательностей комплексных чисел.
В 1959 году в работе [2] Н. М. Коробов предложил квадратурные формулы с параллеле-пипедальными сетками
щ = ({}{%})
^ , к = 1,...,М,
где целые а\,... ,а3 — специально выбранные числа — оптимальные коэффициенты.
Для этих формул на классе Е" выполняется оценка погрешности [/]| = О ('д^); гДе 7 зависит только от а и
Точность найденных формул численного интегрирования значительно превосходит точность как классических, так и вероятностных формул при определенном соотношении между величинами % а и Более того, полученная оценка тем точнее, чем больше гладкость рассматриваемой функции.
Использование параллелепипедальных сеток не ограничивается только численным интегрированием.
В 1961 году В. С. Рябенький в работе [8] предложил численный метод решения задачи Коши для следующего класса дифференциальных уравнений с частными производными:
% = 1 ™
0 ^ Ь ^ Т, -ж <х„ < ж (и = 1,..., в),
и(0,х) = <р(х), х = (х\,... ,х3), (3)
где
п (А А ^ _ ^ ^ ^ т
^дх1,...,дх^ = - ¿0 охп ... дхи (4)
— дифференциальный оператор порядка п(О) = п\ + ... + п3, с максимальным порядком по отдельным переменным, не превосходящим т(О) = тах(щ,..., п3), а <^>(х) = <р(х\,..., х3) — периодическая с периодом единица по каждому из своих аргументов функция из класса Е" (а > т(0>) + 1).
В своей работе В. С. Рябенький предложил некоторый общий подход численного решения задачи Коши с использованием произвольных сеток, для которых выполнены специальные
условия, и показал, что его конструкция применима для многомерных кубических сеток, которые ещё называют равномерными, и для параллелепипедальных сеток Н. М. Коробова.
В работе [6] Н. М. Коробов рассмотрел решение частной задачи Дирихле с нулевым граничным условием для уравнения Пуассона с правой частью из класса Е":
д2и д2и .. „
Щ + ■■■ + ощ = (5)
и(х), ¡(х) еЕ?, (6)
Цхх,..., —х^, ...,х3) = —/(х\, ...,ху ,...,х3) (з = 1,..., в),
и(х) = 0 при х\(1 — х\) ■ ... ■ х3(1 — х3) = 0. (7)
Таким образом,
и(х) \дСв =0, С3 = [0; 1)в.
Метод Н. М. Коробова состоял в получении с помощью параллелепипедальных сеток М^ = ,..., { (к = 0,..., N — 1) приближенного решения указанной задачи Дири-
хле на основании точного решения в виде ряда Фурье, которое легко выписать по ряду Фурье периодической функции /(х). Необходимость приближенного решения обусловлена тем, что, как правило, явного вида ряда Фурье неизвестно, а потому, точное решение имеет только теоретическое значение.
Другой подход к решению дифференциальных уравнений в частных производных использовали Я. М. Жилейкин [4, 5] и В. Т. Стоянцев [9]. В своих работах они сводили задачу решения дифференциального уравнения к задаче вычисления интегралов, которые, в свою очередь, находили с помощью параллелепипедальных сеток.
В 1980 году Л. А. Книжнерман [1] применил вариационный метод решения задач для уравнений в частных производных. Суть его заключается в следующем.
Рассматривается задача Ьи = /, где и — элемент линейного нормированного пространства решений и, / — элемент гильбертова пространства правых частей Е, Ь : и ^ Е — линейный
Ь
Ь и =
и' С и. В качестве приближённого решения берётся функция и, минимизирующая квадрат нормы невязки ||Ьи — f\\'p.
В данной работе будем использовать прямой вариационный метод — метод Ритца. Рассмотрим задачу
Ьи(х) = /(х), (8)
и(х)\эсв =9(х), (9)
где /(х), д(х) е Е", Ь — линейный дифференциальный оператор с постоянными коэффициентами, С3 — единичный куб [0; 1]5.
Её решение сводится к отысканию минимума функционала
ь(и(х)) = У Е (х,и,иХ1 ,...,иХв )йх\...йх3 (10)
при заданных граничных условиях (см. напр. [10]).
Значения функционала ь(и(х)) в методе Ритца рассматриваются не на множестве всех допустимых функций и(х), а на линейных комбинациях
п
и(х) = Шо(х) + (х),
к=1
где (х) — некоторые базисные функции, которые будем находить с помощью теоретико-числовой интерполяции, причём Шо(х) — функция, удовлетворяющая заданным граничным условиям, а остальные (х) удовлетворяют однородным граничным условиям.
На этих полиномах функционал (10) превращается в функцию ^({о) от коэффициентов ..., {п. Эти коэффициенты выбираются так, чтобы функция ^(ги) достигала экстремума. При некоторых ограничениях на функционал ь(и(х)) и базисные функции (х) получим приближённое решение задачи (8).
2. Интерполяционные полиномы по параллелепипедальным сеткам
В 2005 году в работе [3] Н. М. Добровольский, А. Р. Есаян, О. В. Андреева и Н. В. Зайцева предложили способ построения интерполяционного полинома периодической функции многих переменных, заданной в узлах параллелепипедальной сетки.
Их подход заключается в выборе полной решётки Л С . Данной решётке ставятся в соответствие параллелепипедальная сетка М (Л) и конечное множество целочисленных векторов М*(Л) — полная система вычетов фундаментальной решётки Zíi относительно Л. Тогда, если известны значения периодической функции /(х) в узлах сетки М(Л), то тригонометрический полином
1а №= Е с* = N Е ¡(У)е ~2т{гЛ,1})
т ем *(Л) нем (Л)
является интерполяционным для функции /(х).
Данный раздел посвящён описанию фактов, необходимых для построения интерполяционных многочленов по параллелепипедальной сетке.
Пусть Л1,..., Л 8 — линейно независимая система векторов вещественного арифметического пространства М5. Совокупность Л всех векторов вида
+ ... + а3Л3,
где аз независимо друг от друга пробегают все целые рациональные числа, называется решеткой в М8, а сами векторы Л1,..., Л5 — базисом этой решетки.
Л Л*
Л* = {х е М-ГО е Л (х,у) е %3}. (11)
Л
кой М(Л) называется множество М(Л) = Л* ПС8.
Пусть базис Ли = (аи 1,..., аиз) (и = 1,..., в) решётки Л задан матрицей
А =
( аи ■ ■■ аь \
\ аз1 ■ ■ ■ а.3.з )
А = 0.
Величина det А называется детерминантом решётки det Л и не зависит от выбора базиса решётки.
Тогда базис взаимной решётки Л* задаётся матрицей
(
А-1 =
Алл
detA
_А_
detA
\
I- л 1
А31
\ detA
detA
det Л = 0,
/
где — алгебраическое дополнение к элементу Тогда сетка М (Л) имеет вид
к1А1 V + ... + к3А3и
М (Л) = { X
0 =
< 1;и = 1,..., s; к G Z;
}.
det Л
Рассмотрим теперь линейное сравнение
х\ + а2х2 + ... + asxs = 0 ( (mod N)), (12)
где (aj, N) = 1 для всех j = 1,..., s — 1.
Л
/ N 0 0 ••• 0 \
А =
—fl>2 1 0 —as 0 1
\ — as 0 0
1
Детерминант этой решётки равен % а соответствующая ей сетка М(Л) имеет вид
М (Л)
= j ^ к j ka2 ^ J kas
N J ' { N
к = 0,..., N — 1
Усеченной нормой вектора х называется величина q(x) = Х\ ■ ... ■ xs, где для вещественного х обозначаем х = max(1, |х|). Гиперболический параметр д(Л) решетки Л определяется равенством
q(A) = min ^q(x).
х^Л, х=0
Он имеет простой геометрический смысл: гиперболический крест К(Т) не содержит ненулевых точек решетки Л при Т < д(Л).
Гиперболическим крестом называется область
К(Т) = {х | q(x) < Т},
Т
г-и компонентой гиперболического креста К (Т) называется подмножество
Кг(Т) = {х | q(x) ^ Т, ровно г координат х отличны от 0}. Ясно, что справедливо следующее разбиение гиперболического креста:
К (Т ) =
Если М с Л — подрешетка решетки Л, т0 величипа D = det М/ det Л называется индексом подрешетки М решетки Л. Два вектopa х и у решетки Л сравнимы то подрешетке М (находятся в одном классе относительно подрешетки М), если х — у £ М. В этом случае пишем х = у (mod М). Индекс D подрешетки М решетк и Л равен числу классов решетки Л относительно М
Лемма 1. Если решётка Л является решёткой линейного сравнения (12), то два вектора XX,уу € Zs сравнимы по подрешётке Л тогда и только тогда, когда
х\ + а2Х2 + ... + asXs = yi + а2у2 + ... + asys ( (mod N)).
Л
решетки М называется полной системой вычетов решетки Л относительно М.
Обычно полную систему вычетов фундаментальной решетки Zs относительно целочисленной решетки Л будем обозначать через М*(Л), хотя она определена неоднозначно. Ниже будут сформулированы дополнительные условия для выбора М*(Л).
Определение 2. Полная система вычетов фундаментальной решетки Zs относитель-
Л
вычетов, если минимальный гиперболический крест,, содержащий эту полную систему вычетов, имеет минимальное значение своего параметра для всех полных систем вычетов фундаментальной решетки Zs относительно целочисленной решетки Л.
Определение 3. Полная система вычетов фундаментальной решетки Zs относитель-
Л
усеченной нормой среди всех элементов класса вычетов, называется абсолют,но минимальной гиперболической полной системой вычетов.
Такая полная система вычетов обозначается через М* (Л). Вообще говоря, полная система вычетов М* (Л) определена неоднозначно. Это видно на примере решетки NZs при четном N. Действительно, N2 = —N2 (mod N). Уже в одномерном случае две полные системы вычетов {—Ni,..., 0,..., N2} ш {—N2,..., 0,..., Ni} удовлетворяют определению 3. В s-мерном случае таких систем будет 2s. Для однозначности выбора М* (Л) можно еще ввести лексикографический линейный порядок на Zs. Тогда из нескольких возможных элементов с одинаковым значением усеченной нормы выберем наименьший в смысле лексографического упорядочивания. Тем самым М* (Л) будет определено однозначно.
Теорема 1. Для любой функции f(X) на М(Л) полином
Ы f(X))= Е Crne2m^, (13)
т ем *(А)
где
ст = N Е ^(У)е ~2т{т'У1 (14)
уем (А)
является интерполяционным для функции f(X).
Доказательство. См. [3] □
Лемма 2. Пусть значения функции f(X) € R при X € М(Л) и множество М*(Л) симметрично относительно начала, координат,. Тогда, интерполяционный многочлен 1л(f(X)) — действительнозначная функция.
Доказательство. Пользуясь тем, что f(X) € R и (е-2жг(т,'у)) 1 = е~2^г(т,у) получим
с~т = N Е f(y)e-2т(-т'уУ) = N Е f(y) (е~2т(т,уУ)У1 = N Е f(у)е~2жг(т'уУ) = ~ст.
Уем (А) уем(А) уем (Л)
Из равенства = ст следует утверждение леммы. □
Тригонометрический полином (13) зависит от полной системы вычетов М* (Л) решетки Z3 по подрешетке Л. Выбор М* (Л) зависит от класса функций, для которого рассматривается данная задача.
Л
называется наибольшее натуральное число дз(Л), такое, что все целые точки гиперболического креста К(дз(Л)) содержатся в полной системе вычетов фундаментальной решетки Ъ3
Л
Л
Для класса Е" верна следующая теорема.
Л
болической полной системы вычетов М* (К) фундаментальной решетки Ъ3 относительно подрешетки Л для любой функции ¡'(х) € Е" справедлива оценка, сверху погрешности интерполирования:
И^ г II (231п3-\дз(Л)) , ^ 3-2, ЛЛЛ
||¡(х) - 1л(/(хШс < ЫЛ))а-1{ (8 - ща- 1) + О (1п (. (15) □
3. Решение задачи в общем случае
Опишем построение функций Шк(х) (к = 0,... ,п).
Пусть 5 С Ъ3 — произвольное конечное множество целочисленных векторов гп, симметричное относительно начала координат.
Обозначим через Т( 5) пространство всех тригонометрических полиномов с комплексными коэффициентами
Ст. ^ с
Т(5) = \ Р(х) = £ сте2™^ ст € с! .
I тея )
Также обозначим через То(5) С Т( 5) пространство всех функций, принимающих значение пуль на границе «-мерного единичного куба С3 = [0; 1]3.
Найдём сначала на пространстве Т( 5) множество функций удовлетворяющих граничному условию
и(х)\дс: = 9(й), 9(х) = Е Лгпе2тт). (16)
т ея
Обозначим через 73*1 множество всех целочисленных векторов = (]\,..., ]3), каждый из которых имеет координаты, образующие перестановку чисел 1, 2, ..., в такую, что 1 ^ Л < ... < Зг , 1 ^ Л+1 < ... < ]3 ^ 5. Таким образом |73*4| = С3 и 7*0 = 7*3 = {(1, 2,..., 8)}. Также обозначим
5 (¿3,?) = {(т,.. .,тк )|(т1, ...,т3) € Б}.
Для дальнейшего потребуется явный вид ,7* 3_ 1
7*3-1 = {(2,...,8, 1), (1, 3,...,8, 2),..., (1,...,8 - 2,8,8 - 1), (1,..., 5)} ,
|73*,3-1| = в.
В работе [7] доказана следующая лемма.
Лемма 3. Для и(х),(р(х) е Т(5):
¡,(х) = Е сше2™^, <р(Х) = Е (1ше2™^
и(
т ея rn.es
граничное условие
и(х) = р(х) (X е дЩ (17)
выполняется тогда и только тогда, когда для любого ¿3,¿¡^ е J* 3_ 1 и
(,...,т^а_1) е Б (¿3,3-1)
выполняются соотношения
Е = Е (18)
Заметим, что коэффициенты функции <^>(х) являются решением этой системы, а решение соответствующей однородной системы
Е = 0 ¿3,3-1 е Кз-ъ (тл.., тзэ-1) е 5 (¿8,8-1) (19)
яев,
зависит только от множества 5.
Теорема 3. Множество функций на пространстве Т( 5) удовлетворяющих граничному условию (16) представляет собой линейное многообразие и имеет вид
г
и(р,х) = Е ст(Р)е2™т*>, ст(р) = Е Р*-® + (т, (20)
тея г=1
где г — размерность пространства решений системы (19), — некоторый его базис
(Ь = 1,... ,г, т е Б), рг е С — произвольные коэффициенты.
Доказательство. Обозначим
^(х) = Е т^е2™^, г = 1,...,г
т ея
— базис пространства То(5).
Тогда множество функций, удовлетворяющих граничному условию (16) будет иметь вид
и(р, х) = Е р№„(х) + ч>(х) = Е р* Е т^е2™^ + Е (те2™^ = г=1 г=1 тея тея
= Е (¿Р*-® + ат) е2™^, (21)
тея Vг=1 )
где р— произвольный вектор комплексыых чисел. □
Обозначим 5 + С 5 — подмножество векторов в 5, первая ненулевая координата которых
положительна.
Пусть теперь через ТМ(5) С ¥(<5) — пространство вещественнозначных функций пространства Т( 5). Поскольку для таких функций выполнено условие с-г^ = с^, т € 5, то
IР(*) = £ [■
[ т ея+
ТМ( 5) = < Р(х) = V сов (2к(т,х)) + ^ вт (2к(т,х))
Яш ,Рт € ^
где а^ = 2И,е с^, = -21тс
Рассмотрим задачу минимизации функционала
ь(и(х)) = У (х,и,иХ1 ,...,иХв) йх\ ...йх3 (22)
св
на пространстве функций ТМ( 5).
Подставим в (22) тригонометрические полиномы и(р, х), найденные в предыдущем разделе. В этом случае и(и(х)) превращается в функцию <р(р). Далее задача сводится к отысканию коэффициентов при которых достигается минимум этой функции. Метод нахождения этих коэффициентов зависит от вида функционала и(и(х)).
4. Задача Дирихле для уравнения Пуассона
В качестве примера использования описанного метода рассмотрим задачу Дирихле для уравнения Пуассона:
Аи(х) = ¡(х), (23)
и(х)\эсв =9(х). (24)
Решение этой задачи сводится к отысканию минимума функционала
/ 0;/Й (^ )'"
с граничным условием и(х)\дСз = д(х).
Выберем решётку Л и построим интерполяционные полиномы /д ¡(х) и (х) функций ¡(х) и д(х) соответственно. Тогда множество 5 = М*(Л) и эти интерполяционные полиномы будут принадлежать пространству ТМ( М*(Л)). Заменим задачу поиска минимума функционала (25) на пространстве всех допустимых функций задачей поиска минимума данного функционала на пространстве ТМ( М*(Л)).
Теорема 4. Для функций и(х), ¡(ох) € ТМ(5); заданных тригонометрическими многочленами
¡,(х) = ^ с^е /х = ^ ЬгЛ
и(
rn.es rn.es
верно равенство
(и(х)) = 4тт2 £ |с^|2т2 + 2 ^ с^Ь*. (26)
гп ев rn.es
Доказательство. Так как
ди , . Г)тГ1 ( гт ГУ I { ди\ .2
С т С
2т Е с^тье2™(д^) = —4тг2 £ Сгйснткпке2т(т+пп), к = 1,...,з, т ея ^ к' т,пея
(х) / (х) = Е Ьтсле2т,
и(
тлея
то из равенства
/ ... / е 2-(т +П'Х)йх1 ..^ = -т
) ) о, при п = — т,
^
(27)
т е (28)
найдем интеграл:
(и(Х)) =! I —4тг2 Е ст сптк пье2™ (т+Н'£) + 2 £ Ьт<ле2т (гЛ+Л'^) I йхх... йх3 =
к=1гП'ПеЯ ГП'Н ея
= Е стс-ттк(—тк) + 2^2 °ть-т.
к=1т ея т ея
(29)
Поскольку функции и(х), /(х) е ТМ(5), то т е Б с-г^ = ст и Ь-т = Ьт, из чего следует утверждение теоремы. □
Так как, из условия (20) коэффициенты ст = ст(рР), то, обозначив рк = <к + получим многочлен
<р(а,(3) = |ст(р)\2т2 + 2^2 ст(р)Ьт, ст(р) = ЕРг^^ + йт.
т ея т ея г=1
Коэффициенты а, Р, при которых достигается данного многочлена находятся из системы линейных уравнений
* д^ = о, д^ = 0 к = 11,..,г.
д<к ' дРк ' ' '
5. Заключение
В работе рассмотрен подход к решению первой краевой задачи на классах переодических функций Е1^ в мерном кубе. Дело в том, что параллелепипедальные сетки наиболее эффективны именно для таких функций. Однако, от условия периодичности можно избавиться, если воспользоваться одним из вариантов периодизации, предложенных Н. М. Коробовым (см. напр. [6]).
Теоретико-числовые сетки также используются для вычисления интегралов на областях отличных от куба. Отдельный интерес представляет использование данного подхода к решению дифференциальных уравнений в таких областях.
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1. Книжнерман, Л. А. Приложение метода оптимальных коэффициентов к численному решению уравнений в частных производных: Дис. ... канд. физико-математические науки: 01.01.06, 01.01.07. - М.: РГБ, 2006.
2. Коробов, Н. М. О приближённом вычислении кратных интегралов // ДАН СССР. 1959. Т. 124. № 6. С. 1207-1210.
3. Добровольский Н. \!.. Есаян А. Р., Андреева О. В., Зайцева Н. В. Многомерная теоретико-числовая Фурье интерполяция // Чебышевский сборник, 2004 Т. 5. Вып. 1(9). Тула, Из-во ТГПУ им. Л. Н. Толстого. С. 122 - 143.
4. Жилейкин Я. М. О приближённом решениии задачи Дирихле для уравнения Лапласа.-докл. АН СССР, 1964, т. 155, № 5, с. 999-1002.
5. Жилейкин Я. М. О методе приближённого решения задачи Дирихле для уравнения Лапласа в прямоугольном параллелепипеде.- Журн. вычислительной математики и матем. физики, 1965, т. 5, № 2, с. 345-347.
6. Коробов Н. М. Теоретико-числовые методы в приближенном анализе. М.: Физматгиз, 1963.
7. А. В. Родионов. О методе Н. М. Коробова приближенного решения задачи Дирихле // Чебышевский сборник, 15:3 (2014), 48-85.
8. Рябенький В. С. Об одном способе получения разностных схем и об использовании теоре-тикочисловых сеток для решения задачи Коши методом конечных разностей // Тр. матем. ин-та им. В. А. Стеклова. 1961. Т. 60. С. 232-237.
9. Стоянцев В. Т. Решение задачи Коши для параболического уравнения методом квази Монте-Карло.- Журн. вычислительной математики и матем. физики, 1973, т. 13, № 5, с. 1153-1160.
10. Эльсгольц Л. Э. Дифференциальные уравнения и вариационное исчисление. М.: Наука, 1969. 425 с.
11. W. Ritz. Uber eine neue Methode zur Losung gewisser Variationsprobleme der mathematischen Phvsik // Journal fur die Reine und Angewandte Mathematik, 1909, vol. 135, pages 1—61.
REFERENCES
1. Knizhnerman, L. A. "Application of the method of optimal coefficients to the numerical solution of partial differential equations": Dis. ... cand. physical and mathematical sciences: 01.01.06, 01.01.07. - M.: RGB, 2006.
2. N. M. Korobov, 1959, "On the approximate calculation of multiple integrals", DAN SSSR, vol. 124, no. 6, pp. 1207-1210.
3. N. M. Dobrovolskv, A. R. Yesavan, О. V. Andreeva, N. V. Zaitseva, 2004, "Multidimensional number-theoretic Fourier interpolation", Chebvshevskii sbornik, vol. 5, iss. 1(9), pp. 122-143.
4. Zhileikin Ya. M. "On the approximate solution of the Dirichlet problem for the Laplace equation". - Dokl. AN SSSR, 1964, vol. 155, no. 5, p. 999-1002.
5. Zhileikin Ya. M. "On the method of approximate solution of the Dirichlet problem for the Laplace equation in a rectangular parallelepiped".- Zh. Computational Mathematics and Mathematics. Physics, 1965, vol. 5, no. 2, p. 345-347.
6. N. M. Korobov, 1963, "Number-theoretic methods in approximate analysis", Moscow: Fizmatgiz.
7. А. V. Rodionov, 2014, "About N. М. Korobov's method of approximate solution of the Dirichlet problem", vol. Chebvshevskii sbornik, vol. 15:3, pp. 48-85.
8. V. S. Rvabenkv, 1961, "On one method for obtaining difference schemes and on the use of number-theoretic grids for solving the Cauchv problem by the finite difference method", Tr. matem. V. A. Steklov Institute, vol. 60, p. 232-237.
9. V. T. Stovantsev, "Solution of the Cauchv problem for a parabolic equation by the quasi Monte Carlo method", Zh. Computational Mathematics and Mathematics. Physics, 1973, vol. 13, no. 5, p. 1153-1160.
10. Elsgol'ts L. E. Differential Equations and the Calculus of Variations. M.: Nauka, 1969. 425 p.
11. W. Ritz. "Uber eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik", Journal für die Reine und Angewandte Mathematik, 1909, vol. 135, pages 1—61.
Получено: 24.07.2022 Принято в печать
