close

Вход

Забыли?

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

?

Решение задачи Дирихле для уравнения колебаний струны методом сопряженных градиентов.

код для вставкиСкачать
В. И. Васильев, А. М. Кардашевский, В. В. Попов. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ МЕТОДОМ
СОПРЯЖЕННЫХ ГРАДИЕНТОВ
УДК 519.63
В. И. Васильев, А. М. Кардашевский, В. В. Попов
РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ
ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ
МЕТОДОМ СОПРЯЖЕННЫХ ГРАДИЕНТОВ
Рассматривается неклассическая задача Дирихле для уравнения колебания струны, описываемое начально-краевой
задачей для гиперболичеcкого типа второго порядка. Для ее численного решения используется итерационное уточнение
начального условия для производной по времени, предложенное в работе [1] для решения ретроспективной обратной задачи
теплопроводности. В данной работе предлагается наиболее быстро сходящийся итерационный метод сопряженных градиентов. Приведены примеры расчетов для модельных задач, в том числе и со случайными погрешностями во входных данных.
Ключевые слова: обратные задачи, гиперболическое уравнение второго порядка, задача Дирихле, конечно-разностный
метод, явная разностная схема, метод сопряженных градиентов, случайные погрешности.
V. I. Vasilyev, A. M. Kardashevsiy, V. V. Popov
The Dirichlet Problem
for Oscillation Equation of String
by Conjugate Gradient Method
The article considers a non-classical Dirichlet problem for the oscillation equation of string described by an initial-boundary value
problem for the second-order hyperbolic equation. For numerical solution there is used the iterative improvement of the initial condition
for the time derivative, which proposed in [6] for the solution of a retrospective inverse heat conduction problem. The author proposes the
most rapidly converging iterative conjugate gradient method. Examples of calculations for model problems including one with random
errors in input data are presented.
Keywords: inverse problems, second-order hyperbolic equation, Dirichlet problem, finite-difference method, explicit difference
scheme, conjugate gradient method, random errors.
ВАСИЛЬЕВ Василий Иванович – д. ф.-м. н., профессор,
зав. каф. вычислительных технологий ИМИ СВФУ им.
М. К. Аммосова.
E-mail: vasvasil@mail.ru
VASILYEV Vasiliy Ivanovich – D. Sc., Professor, Head of
the department «Computational Technologies» of Institute of
Mathematics and Informatics, North-Eastern Federal University
n. a. M. K. Ammosov.
E-mail: vasvasil@mail.ru
КАРДАШЕВСКИЙ Анатолий Михайлович – зав.
кабинетом каф. «Вычислительные технологии» ИМИ СВФУ
им. М. К. Аммосова.
E-mail: kardam123@gmail.com
KARDASHEVSKIY Anatoliy Michaylovich – Head of Cabinet
of the department «Computational Technologies» of Institute of
Mathematics and Informatics, North-Eastern Federal University
n. a. M. K. Ammosov.
E-mail: kardam123@gmail.com
ПОПОВ Василий Васильевич – доцент-исследователь каф.
вычислительных технологий ИМИ СВФУ им. М. К. Аммосова.
E-mail: pvvchoxoon@gmail.com
POPOV Vasily Vasilievich – Associate professor and researcher
of the department «Computational Technologies» of Institute of
Mathematics and Informatics, North-Eastern Federal University
n. a. M. K. Ammosov.
E-mail: pvvchoxoon@gmail.com
Введение
Задача Дирихле для уравнения колебаний струны
относится к классу условно корректных задач математической физики [2-3]. Доказательство корректности
постановки задачи Дирихле для уравнений гиперболического типа второго порядка можно найти в книгах [24]. Общий подход к решению таких некорректных задач
базируется на использовании градиентных итерационных методов. Основу этих методов составляет нахождение приближенного численного решения некорректных
задач из итерационной минимизации соответствующего функционала. В цикле работ С. И. Кабанихина и его
учеников [5-6] численно исследовано решение задачи
Дирихле для волнового уравнения методами итераций
Ландвебера и наискорейшего спуска, проведено теоретическое исследование их устойчивости и сходимости
приближенного решения к решению исходной дифференциальной задачи для одномерного и двумерного
уравнений гиперболического типа второго порядка.
Для задачи восстановления начального условия
уравнения теплопроводности в работе [1] предложен
быстро сходящийся итерационный метод, продемонстрирована его высокая эффективность при решении
модельных задач. В работе [7] нами предложен итера43
ВЕСТНИК СВФУ, 2015, том 12, № 2
ционный метод решения неклассической задачи для
уравнения колебаний струны, в которой по заданному
распределению смещения струны в конечный момент
времени определяется скорость струны в начальный
момент времени.
В данной работе для численного решения задачи
Дирихле для гиперболического уравнения второго порядка предлагается использовать аналогичный подход
с итерационным уточнением начального условия. Для
поставленной неклассической задачи по заданному
решению на конечный момент времени методом сопряженных градиентов уточняется начальное условие для
производной по времени искомого решения. Приведены
примеры расчетов для модельных задач, в том числе
со случайными погрешностями во входных данных.
Результаты вычислительного эксперимента демонстрируют высокую эффективность предлагаемого итерационного метода.
Решение прямой задачи для уравнения колебаний
струны
Колебание тонкой струны описывается дифференциальным уравнением с частными производными
гиперболического типа второго порядка
2
2
∂ u ∂ u
(1)
−
= 0, x ∈ (0, l ) , 0 < t ≤ T .
∂t 2 ∂x 2
В предположении, что концы колеблющейся струны
x=0, x=l закреплены, выполняются однородные граничные условия первого рода
u (0, t ) = 0, u (l , t ) = 0, 0 < t ≤ T .
(2)
В начальный момент времени заданы смещение и
скорость точек струны
∂u
(3)
( x,0) = ν ( x), 0 ≤ x ≤ l.
∂t
При достаточной гладкости входных данных φ(x),
ν(x) и выполнении условий согласования начальных и
граничных условий задача (1)-(3) поставлена корректно
[8].
Для ее численного решения используем метод
конечных разностей. В области определения задачи (1)(3), представляющей собой прямоугольник [0,l]×[0,T],
введем равномерную прямоугольную сетку
u ( x,0) = ϕ ( x ) ,
ωhτ = ωh ×ωτ ,
где
ωh = { xi = ih, i =0,1,…, n; h = l / n} ,
ωτ = {t j = jτ , j =0,1,…, m;τ = T / m} .
Задаче (1)-(3) на пространственно-временной сетке
ωhτ поставим в соответствие симметричную разностную схему с весовым множителем σ:
44
yij +1 − 2 yij + yij−1
2
=
yi +1 − 2 yi + yi−1
τ
i = 1, n −1, j = 1, m −1,
h2
,
(4)
y0j +1 = 0, ynj +1 = 0, j = 1,m −1,
(5)
yi0 = ϕ ( xi ) , yi1 = yi0 + τν ( xi ) , i = 0,n,
(6)
здесь использовано обозначение
yi = σ yij +1 + (1− 2σ ) yij + σ yij−1 .
На множестве сеточных функций y Î H таких, что
y(x)=0, x Ï ωh , определим сеточный оператор A соотношением
Ay = − y xx , x ∈ ωh .
Хорошо известно [9], что для решения разностной
схемы (4)-(6) при выполнении условия устойчивости
1
1
σ≤ − 2
4 τ || A ||
(7)
справедлива следующая априорная оценка:
|| Y ||∗ ≤|| Y ||∗ ≤…≤|| Y (τ ) ||∗ , t ∈ ωτ ,
норма || Y ||* определяется следующим образом:
(8)
2
1
1
Y * = || y + y ||2A +τ 2 || yt ||2 1 , R = E + σ A.
R− A
4
4
4
Тем самым норма решения задачи Коши со временем не возрастает. Отметим, что для явной разностной
схемы (σ=0) условие устойчивости (7) принимает простейший вид τ≤h. Таким образом, на решении исходной
дифференциальной задачи (1)-(3) разностная схема (4)(6) обладает вторым порядком аппроксимации, и при
выполнении условия
1 h2
σ≤ − 2
4 4τ
ее решение при достаточной гладкости входных данных
сходится к решению дифференциальной задачи (1)-(3)
со скоростью порядка 0(τ2+h2). Приведем результаты
вычислительного эксперимента. Сначала проведем расчеты, показывающие влияние шагов сетки на точность
определения приближенного решения прямой задачи
для уравнения колебаний струны. Для этого рассмотрим задачу с точным решением
u ( x, t ) = sin ( x ) cos (t ) , x ∈  0, l ], t ∈[ 0, T  ,
где l=π, T=π/2.
На рис. 1 приведен график зависимости от шагов
пространственно-временной сетки точности решения
явной разностной схемы на решении исходной задачи в
финальный момент времени
В. И. Васильев, А. М. Кардашевский, В. В. Попов. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ МЕТОДОМ
СОПРЯЖЕННЫХ ГРАДИЕНТОВ
Рис. 1. Решение прямой задачи для уравнения колебаний струны
m
|| z ||=
n−1
∑
i =1
( yim
2
− u ( xi , tm )) h .
Расчеты проводились на последовательности сеток,
на графике количество узлов пространственной сетки
через k вычисляется по формуле n=20∙2k , т. е. при n=20,
40, 80, 160, 320, 640. Сплошная линия (1) соответствует
расчету с максимально допустимым временным шагом,
обеспечивающим устойчивость разностной схемы τ=h,
следовательно, m принимает те же значения, что и n.
Прерывистой линией (2) представлена та же зависимость при τ=h/2, а штрих-пунктирной линией (3) – при
τ=h/4. То, что графики представляют собой прямые
линии, подтверждает теоретическое утверждение о втором порядке сходимости по h и τ решения разностной
схемы к решению исходной дифференциальной задачи.
Представленные графики также показывают, что при
решении задачи Дирихле критерий выхода из итерационного цикла на данной конкретной пространственно-временной сетке следует согласовывать с точностью
решения прямой задачи.
Численную реализацию разностной схемы (4)-(6)
проведем для задачи (1)-(3) c однородными граничными
условиями, начальным условием для скорости смещения струны µ1(t)≡0, µ2(t)≡0, 0<t≤T, φ(x)≡0, 0≤x≤l и при
неоднородном начальном условии для смещения
2
ϕ ( x) = 12e−50( x−l / 2) , 0 ≤ x ≤ 1.
На рис. 2 приведены результаты решения прямой
задачи по явной разностной схеме (σ=0) при l=1, T=0,99
на последовательности сгущающихся сеток c (n, m)=(25,
25); (50, 50); (100, 100); (200, 200). Здесь слева даны графики функции φ(x), а в середине – график решения в
финальный момент времени yim , i = 0,n . Отметим, что
пространственно-временная сетка построена таким
обра­зом, что условие устойчивости τh выполняется с
минимальным запасом. Справа на рис. 2 представлено
решение yij , i = 0,n, j = 0,m разностной схемы (n=200,
m=200).
Результаты численного эксперимента показывают,
что достаточно хорошую точность получаем на достаточно грубых сетках с n=50, m=50. В дальнейшем для
задания значения решения в конечный момент времени
для неклассической задачи будем использовать решение
прямой задачи в финальный момент времени, полученное на достаточно подробной сетке.
Постановка задачи Дирихле для уравнения колебаний струны
В прямоугольнике ищем функцию u(x,t) – решение
гиперболического уравнения второго порядка
∂ 2 u ∂ 
∂u 
− k ( x )  = 0, x ∈ (0, l ) , 0 < t < T .
2

∂x 
∂x 
∂t
(9)
Предположим, что концы колеблющейся струны
x=0, x=l закреплены, т. е. выполняются граничные условия первого рода
u (0, t ) = 0, u (l , t ) = 0, 0 < t < T .
(10)
В начальный и финальный моменты времени задаются смещения струны
u ( x,0) = ϕ ( x ) , u ( x, T ) = φ ( x ) , 0 ≤ x ≤ l.
(11)
Пусть коэффициент k(x) дифференциального уравнения (7) удовлетворяет условиям 0<k1≤k(x)≤k2. Задача
Дирихле относится к классу неклассических задач математической физики, является условно корректной
[3]. Она может иметь не единственное решение. Как
показывает нижеследующий простейший пример. Пусть
45
ВЕСТНИК СВФУ, 2015, том 12, № 2
Рис. 2. Решение прямой задачи для уравнения колебаний струны
φ(x)≡0, ϕ(x)≡0, тогда неклассическая задача (9)-(11) при
коэффициенте k(x)≡0, помимо тривиального решения
u(x,t)≡0, имеет семейство решений, задаваемых формулой
u ( x, t ) = c ⋅ sin (k p x )⋅ sin (npt ) , c ∈ R, k , l , T ∈ N .
Разностная задача Дирихле
Обозначим через ωhτ множество внутренних узлов
пространственно-временной сетки
ωhτ = ωh ×ωτ ,
где
ωh = { xi = ih, i =1,2,…, n −1} ,
ωτ = {t j = jτ , j =1,2,…, m −1} .
На множестве сеточных функций y Î H таких, что
y(x)=0, x Ï ωh , определим сеточный оператор A соотношением
Ay = −(a ( x) y x ) x , x ∈ ωh ,
положив, например, a(x)=k(x–0,5h).
46
В сеточном гильбертовом пространстве H скалярное
произведение и норму введем соотношениями
( y, v) = ∑ yvh, ||y ||= ( y, y )1/ 2 .
x∈ω
В H оператор A является самосопряженным, положительно определенным A=A*>0 при k(x)≥k1, k1>0 и
ограниченным ||A||≤4k2/h2.
Сначала от задачи Дирихле перейдем к задаче
определения решения дифференциально-операторного
уравнения
d2y
+ Ay = 0, x ∈ ω , 0 < t < T
dt 2
(12)
при заданных дополнительных условиях, обеспечивающих единственность ее решения
y ( x , 0) = ϕ ( x ) , y ( x , T ) = φ ( x ) , x ∈ ω h .
(13)
При использовании симметричной трехслойной
разностной схемы с весовым множителем σ дискретный
В. И. Васильев, А. М. Кардашевский, В. В. Попов. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ МЕТОДОМ
СОПРЯЖЕННЫХ ГРАДИЕНТОВ
аналог задачи (10)-(11) имеет вид:
∨
y − 2 y + y
∨

+ Aσ y + (1− 2σ ) y + σ y  = 0, ( x, t ) ∈ ωhτ ,
2



τ
(14)
y (0) = ϕ ( x ) , y (tn ) = φ ( x ) , x ∈ ωh ,
(15)
где 0≤σ≤1 и использована безиндексная система обозначений, введенная А. А. Самарским [9]:
∨
y = y j +1 , y = y j , y = y j−1 , x ∈ ω .
h
∂u
( x , 0) = ν ( x ) , x ∈ ω .
∂t
Обозначим через yj разностное решение на момент
времени tj =jτ, где τ>0 – шаг по времени, причем mτ=T.
Дискретный аналог прямой задачи имеет вид:
∨
∨
y  = 0,

y (0) = ϕ ( xi ) , y (τ ) = y (0) + τν ( x ) , x ∈ ωh .
(16)
(17)
Итерационный метод
Для нахождения решения дискретной задачи Дирихле
(12)-(13) используем наиболее быстро сходящийся итерационный метод вариационного типа – метод сопряженных градиентов [10], основанный на последовательном уточнении искомого начального условия ν(x) с дальнейшим решением на каждой итерации дискретного
аналога прямой задачи (16)-(17). Придадим этой задаче
соответствующую операторную формулировку. Последовательно исключая промежуточные значения сеточной
функции y (t ) , t ∈ ωt при заданных φ и ν на конечный
момент времени, получим
y N = Aν + Bϕ ,
где ,  – соответственно операторные полиномы от
положительно определенного и самосопряженного оператора A. С учетом этого факта приближенному решению обратной задачи естественно сопоставить решение
следующего операторного уравнения
Aν = φ ( x ) − Bϕ = φ .
ν k +1 = α k +1ν k + (1− α k +1 )ν k −1 − α k +1 β k +1rk ,
k = 0,1,...,
Для численной реализации системы линейных
алгебраических уравнений (12)-(13) воспользуемся
итерационным методом сопряженных градиентов,
связанным с уточнением производной по времени в
начальный момент времени. На каждой итерации решается прямая задача для уравнения колебаний струны с
использованием стандартной симметричной трехслойной разностной схемы, обладающей вторым порядком
аппроксимации по h и τ [9].
Итак, пусть вместо обратной задачи рассматривается прямая задача для этого же уравнения, когда вместо
второго из условий используется начальное условие
y − 2 y + y

+ Aσ y + (1− 2σ ) y + σ
2

τ
( x, t ) ∈ ωhτ ,
В силу самосопряженности оператора A самосопряженным является и оператор , . Однако мы ничего не
можем сказать об однозначной разрешимости данного
сеточного уравнения.
Для решения этого уравнения используем трехслойный итерационный метод сопряженных градиентов,
записанный в каноническом виде [9]
(18)
где rk = ν k − φ – невязка.
Для вычисления итерационных параметров βk+1, αk+1
используются следующие рекуррентные формулы
β k +1 =
(rk , rk )
, k = 0,1,...,
( Ark , rk )
−1
 β
(rk , rk ) 1 
k +1

α k +1 = 1−
 ,

β k (rk −1 , rk −1 ) α k 
(19)
(20)
k = 1, 2,...., α1 = 1.
В этом итерационном методе вычисление ykn осуществляется численной реализацией прямой задачи для
очередного приближения искомой сеточной функции νk:
∨
y − 2 y + y
k
k
k
τ2
( x, t ) ∈ ωhτ
∨ 

+ Aσ y k + (1− 2σ ) yk + σ y k  = 0,


(21)
с начальными условиями
yk0 = ϕ , y1k = ϕ + τν k , x ∈ ωh .
(22)
Следовательно, невязка вычисляется по формуле
rk = ykm − φ , x ∈ ωh . Аналогично организуется вычисление вектора zk =Ark:
∨
z k − 2 zk + z k
2
τ
( x, t ) ∈ ωhτ
∨ 

+ Aσ z k + (1− 2σ ) zk + σ z k  = 0,


(23)
с начальными условиями
zk0 = 0, z1k = τ rk , x ∈ ωh .
(24)
Таким образом, итерационный процесс реализуется
в следующем порядке:
1. Полагаем k=0 и задаем начальное приближение
искомой функции ν k , x Î ωh .
2. Запускаем счетчик итераций k=k+1, последовательно решая прямые задачи (21)-(22) и (23)-(24), определяем, соответственно, невязку rk = ykm − φ ,� � x ∈ ωh и
вектор zk = rk .
3. По рекуррентным формулам (19)-(20) определяем
значения итерационных параметров βk+1, αk+1.
4. По формуле (18) находим очередное приближение
искомого начального условия ν k +1 ( x), x ∈ ωh .
47
ВЕСТНИК СВФУ, 2015, том 12, № 2
5. Этот процесс продолжается до тех пор, пока не
выполнится критерий остановки итераций rk<ε, иначе
возвращаемся к пункту 2.
Примеры расчетов
Вычислительную эффективность предложенного
метода удобно проиллюстрировать на примере численного решения простейшей задачи Дирихле для уравнения колебания струны. Будем искать приближенное
решение задачи Дирихле для уравнения
∂ 2u ∂ 2u
−
= 0, 0 < x < 1, 0 < t < T ,
∂t 2 ∂x 2
с однородными граничными условиями
u (0, t ) = 0, u (1, t ) = 0, 0 < t < T ,
и условиями в начальный и финальный моменты времени
u ( x, 0) = 0, u ( x, T ) = φ ( x ) , 0 ≤ x ≤ 1.
Пусть h – шаг равномерной сетки по пространственной переменной
ωh = { x | x = xi = ih, i = 0,1,..., n; nh = 1}.
В рамках квазиреального вычислительного эксперимента ограничимся примером численного решения
обратной задачи Дирихле, которое соответствует
решению прямой задачи для уравнения колебаний с
теми же граничными условиями и двумя начальными
условиями
u ( x, 0) = 0,
∂u ( x , 0 )
= ν ( x ) , 0 ≤ x ≤ 1.
∂t
Из решения прямой задачи определяется функция
φ ( x) , которая и присутствует в постановке задачи Дирихле:
u ( x, T ) = φ ( x ) , 0 ≤ x ≤ 1.
В качестве начального приближения искомого
начального условия возьмем функцию, отдаленно напоминающую искомую функцию
ν 0 = 5 sin (π x ) , 0 ≤ x ≤ 1.
Приведем результаты по решению обратной задачи
Дирихле в условиях, когда искомая функция ϕ(x) имеет
разный вид, на разных сетках, при разных T. В расчетах
задавали ε=0,0001.
На рис. 3 приведены результаты в разные моменты
времени m=100, l=1, T=0,24, 0,48, 0,72, 0,96, 1,2; n=25, 50,
100, 200. Решение в разные моменты времени. Соответственно, итераций 9, 8, 9, 5, 13, а время счета 0,332, 0,563,
0,899, 0,445, 1,436 c.
На рис. 4 приведены аналогичные результаты, когда
искомая функция задана в виде
48
2
2
2
ν ( x) = 7(−e−100( x−l / 4) + e−20( x−l / 2) − e−100( x−3l / 4) ,
0 ≤ x ≤ 1.
В этом случае итерации сходятся значительно медленнее, соответственно, количество итераций 55, 48, 43,
39, 63, а время счета составляет, соответственно, 1,736,
2,012, 2,692, 3,261, 6,592 c.
Наконец, в численных экспериментах сеточная
функция φ ( x ) ,� � x ∈ ωh возмущалась, как в работе [1],
следующим образом:
φδ ( x) = φ + δσ ( x ) , x ∈ ωh ,
где σ(x) – равномерно распределенные на отрезке [–1, 1]
случайные величины. Итерационный процесс сопряженных градиентов обрывался по достижению невязки
значения δ, т. е. при
rn(δ ) < δ .
Проведены расчеты для уровня погрешности во
входных данных, определяемых величиной δ= 0,001
(25). Здесь приведены данные расчетов при использовании предлагаемого итерационного метода без сглаживающего оператора. Для выделения более гладкого решения, принадлежащего вместо L2(ω) сеточному пространству W21 (ω ) , зададим сглаживающий оператор в
виде
Dy = −5 y xx + y, x ∈ ω
сеточных функций, обращающихся в нуль в граничных
узлах.
Графики функций ν(x) и ϕ(x) показаны слева на рис.
5. Здесь приведены результаты счета на пространственной сетке M=200, h=0,005 при помощи явной схемы для
волнового уравнения σ=0 с временным шагом, удовлетворяющим условию устойчивости: τ=0,00495h, T= 0,99,
N=200.
В данном случае итерационный процесс сопряженных градиентов сходится достаточно быстро, всего за
5-15 итераций, поскольку каждый раз генератор случайных чисел выдает новый набор случайных чисел.
В расчетах брали ε=0,007. Очевидно, рассчитывать на
точное восстановление начального ν(x) нельзя в силу
пониженной гладкости входных данных. С хорошей
точностью приближенное решение на конечный момент времени мы имеем от более гладкого начального
условия, которое находится при итерационном решении
обратной задачи при сглаживании входных данных. В
двумерном случае количество итераций значительно
растет. В этом случае для существенного ускорения
скорости сходимости итерационного процесса возникает необходимость использования в методе сопряженных
градиентов предобуславливателя.
Заключение
Приведенные результаты вычислительного эксперимента подтверждают: работоспособность предлагаемого
В. И. Васильев, А. М. Кардашевский, В. В. Попов. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ МЕТОДОМ
СОПРЯЖЕННЫХ ГРАДИЕНТОВ
Рис. 3. Решение в разные моменты времени
Рис. 4. Решение в разные моменты времени
Рис. 5. Восстановленное начальное условие
49
ВЕСТНИК СВФУ, 2015, том 12, № 2
метода решения задачи Дирихле для гиперболического
уравнения, ее достаточно высокую эффективность; несущественное развитие возмущений в исходных данных с
естественным увеличением амплитуды при возрастании
уровня шумов.
Авторы благодарны профессору П. Н. Вабищевичу
за конструктивные замечания и советы.
Работа выполнена при финансовой поддержке
РФФИ (Грант № 13-01-00719).
Литература
1. Самарский А. А., Вабищевич П. Н., Васильев В. И.
Итерационное решение ретроспективной обратной задачи теплопроводности // Математическое моделирование. – 1997. – Т.
9, № 5. – C. 119-127.
2. Samarskii A. A., Vabishchevich P. N. Numerical Methods for
Solving Inverse Problems of Mathematical Physics. De Gruyter,
2007.
3. Kabanikhin S. I. Inverse and ill – posed problems: theory and
applications. De Gruyter, 2012.
4. Денисов А. М. Введение в теорию обратных задач. – М.:
Изд-во Московского ун-та, 1994.
5. Kabanikhin S. I., Bektemesov M. A., Nurseitov D. B.,
Krivorotko O. I., Alimova A. N. An optimization method in the
Dirichlet problems for the wave equation // Journal of inverse and
ill – posed problems. 2012. – № 2 (20). – C. 193-211.
6. Кабанихин С. И., Криворотько О. И. Численный метод
решения задачи Дирихле для волнового уравнения // Сиб. журн.
индустр. матем. – 2012, – Т. 15, № 4. – С. 90-101
7. Васильев В. И., Попов В. В., Еремеева М. С., Кардашевский А. М. Итерационное решение одной неклассичесой задачи
для уравнения колебаний струны // Вестник МГТУ им. Н. Э.
Баумана. Сер. Естественные науки. – 2015. – № 3. – С. 77-87.
8. Тихонов А. Н., Самарский А. А. Уравнения математичес­
50
кой физики. – М: Наука, 1972.
9. Самарский А. А. Теория разностных схем. – М.: Наука,
1989.
10. Самарский А. А., Николаев Е. С. Методы решения сеточных уравнений. – М.: Наука, 1978.
References
1. Samarskii A. A., Vabishchevich P. N., Vasil’ev V. I. Iteratsionnoe
reshenie retrospektivnoi obratnoi zadachi teploprovodnosti //
Matematicheskoe modelirovanie. – 1997. – T. 9, № 5. – C. 119-127.
2. Samarskii A. A., Vabishchevich P. N. Numerical Methods for
Solving Inverse Problems of Mathematical Physics. De Gruyter,
2007.
3. Kabanikhin S. I. Inverse and ill – posed problems: theory and
applications. De Gruyter, 2012.
4. Denisov A. M. Vvedenie v teoriiu obratnykh zadach. – M.:
Izd-vo Moskovskogo un-ta, 1994.
5. Kabanikhin S. I., Bektemesov M. A., Nurseitov D. B.,
Krivorotko O. I., Alimova A. N. An optimization method in the
Dirichlet problems for the wave equation // Journal of inverse and
ill – posed problems. 2012. – № 2 (20). – C. 193-211.
6. Kabanikhin S. I., Krivorot’ko O. I. Chislennyi metod resheniia
zadachi Dirikhle dlia volnovogo uravneniia // Sib. zhurn. industr.
matem. – 2012, – T. 15, № 4. – S. 90-101
7. Vasil’ev V. I., Popov V. V., Eremeeva M. S., Kardashevskii
A. M. Iteratsionnoe reshenie odnoi neklassichesoi zadachi dlia
uravneniia kolebanii struny // Vestnik MGTU im. N. E. Baumana.
Ser. Estestvennye nauki. – 2015. – № 3. – S. 77-87.
8. Tikhonov A. N., Samarskii A. A. Uravneniia matematicheskoi
fiziki. – M: Nauka, 1972.
9. Samarskii A. A. Teoriia raznostnykh skhem. – M.: Nauka,
1989.
10. Samarskii A. A., Nikolaev E. S. Metody resheniia setochnykh
uravnenii. – M.: Nauka, 1978.
Документ
Категория
Без категории
Просмотров
11
Размер файла
1 805 Кб
Теги
решение, методов, уравнения, струна, градиент, сопряженное, дирихле, колебания, задачи
1/--страниц
Пожаловаться на содержимое документа