Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач
В настоящей статье рассматривается проблема приближенного решения жестко некорректных задач с возмущенными правыми частями. Проводится анализ аппроксимационных свойств комбинации конечномерного варианта тихоновского метода регуляризации и апостериорного выбора параметра регуляризации по принципу рав...
Gespeichert in:
Datum: | 2011 |
---|---|
Hauptverfasser: | , |
Format: | Artikel |
Sprache: | Russian |
Veröffentlicht: |
Інститут прикладної математики і механіки НАН України
2011
|
Schriftenreihe: | Український математичний вісник |
Online Zugang: | http://dspace.nbuv.gov.ua/handle/123456789/124430 |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Zitieren: | Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач / С.Г. Солодкий, А.В. Грушевая // Український математичний вісник. — 2011. — Т. 8, № 3. — С. 447-457. — Бібліогр.: 14 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-124430 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1244302017-09-27T03:03:01Z Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач Солодкий, С.Г. Грушевая, А.В. В настоящей статье рассматривается проблема приближенного решения жестко некорректных задач с возмущенными правыми частями. Проводится анализ аппроксимационных свойств комбинации конечномерного варианта тихоновского метода регуляризации и апостериорного выбора параметра регуляризации по принципу равновесия. Установлено, что в рамках этого подхода удается достичь оптимальный порядок точности на исследуемом классе уравнений. Сравнительный анализ описанного метода с известными ранее подтверждает эффективность предлагаемого подхода. 2011 Article Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач / С.Г. Солодкий, А.В. Грушевая // Український математичний вісник. — 2011. — Т. 8, № 3. — С. 447-457. — Бібліогр.: 14 назв. — рос. 1810-3200 2010 MSC. 47A52, 65R30. http://dspace.nbuv.gov.ua/handle/123456789/124430 ru Український математичний вісник Інститут прикладної математики і механіки НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
Russian |
description |
В настоящей статье рассматривается проблема приближенного решения жестко некорректных задач с возмущенными правыми частями. Проводится анализ аппроксимационных свойств комбинации конечномерного варианта тихоновского метода регуляризации и апостериорного выбора параметра регуляризации по принципу равновесия. Установлено, что в рамках этого подхода удается достичь оптимальный порядок точности на исследуемом классе уравнений. Сравнительный анализ описанного метода с известными ранее подтверждает эффективность предлагаемого подхода. |
format |
Article |
author |
Солодкий, С.Г. Грушевая, А.В. |
spellingShingle |
Солодкий, С.Г. Грушевая, А.В. Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач Український математичний вісник |
author_facet |
Солодкий, С.Г. Грушевая, А.В. |
author_sort |
Солодкий, С.Г. |
title |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
title_short |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
title_full |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
title_fullStr |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
title_full_unstemmed |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
title_sort |
об апостериорном выборе параметра регуляризации при решении жестко некорректных задач |
publisher |
Інститут прикладної математики і механіки НАН України |
publishDate |
2011 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/124430 |
citation_txt |
Об апостериорном выборе параметра регуляризации при решении жестко некорректных задач / С.Г. Солодкий, А.В. Грушевая // Український математичний вісник. — 2011. — Т. 8, № 3. — С. 447-457. — Бібліогр.: 14 назв. — рос. |
series |
Український математичний вісник |
work_keys_str_mv |
AT solodkijsg obaposteriornomvyboreparametraregulârizaciiprirešeniižestkonekorrektnyhzadač AT gruševaâav obaposteriornomvyboreparametraregulârizaciiprirešeniižestkonekorrektnyhzadač |
first_indexed |
2025-07-09T01:25:28Z |
last_indexed |
2025-07-09T01:25:28Z |
_version_ |
1837130663565721600 |
fulltext |
Український математичний вiсник
Том 8 (2011), № 3, 447 – 457
Об апостериорном выборе параметра
регуляризации при решении жестко
некорректных задач
Сергей Г. Солодкий, Анна В. Грушевая
(Представлена А. М. Самойленко)
Аннотация. В настоящей статье рассматривается проблема при-
ближенного решения жестко некорректных задач с возмущенными
правыми частями. Проводится анализ аппроксимационных свойств
комбинации конечномерного варианта тихоновского метода регуля-
ризации и апостериорного выбора параметра регуляризации по прин-
ципу равновесия. Установлено, что в рамках этого подхода удае-
тся достичь оптимальный порядок точности на исследуемом классе
уравнений. Сравнительный анализ описанного метода с известными
ранее подтверждает эффективность предлагаемого подхода.
2010 MSC. 47A52, 65R30.
Ключевые слова и фразы. Жестко некорректная задача, апо-
стериорный выбор параметра регуляризации, принцип равновесия,
оптимальный порядок точности.
1. Постановка задачи
Рассмотрим проблему приближенного решения некорректной за-
дачи, представимой в виде операторного уравнения I рода
Ax = y, (1.1)
где A : X → Y — линейный компактный инъективный оператор,
действующий между гильбертовыми пространствами X и Y , причем
множество Range(A) незамкнуто в Y . Для сокращения выкладок обо-
значим скалярные произведения в обоих пространствах через (·, ·) и
соответствующие им нормы через ‖ · ‖. Более того, тем же символом
‖ · ‖ будем обозначать стандартную операторную норму. Из контекс-
та будет ясно, какое именно пространство или норма имеются ввиду.
Статья поступила в редакцию 30.06.2011
ISSN 1810 – 3200. c© Iнститут математики НАН України
448 Об апостериорном выборе параметра...
Кроме того, будем предполагать, что вместо y доступно лишь неко-
торое приближение yδ ∈ Y , ‖y − yδ‖ ≤ δ, δ > 0.
Обычно уравнение (1.1) называют жестко некорректной задачей,
если его решение x0 = A−1y имеет конечную в некотором смысле
гладкость, а A — оператор бесконечной гладкости. Характерной чер-
той таких задач является то, что решение x0 принадлежит некото-
рому подпространству V , непрерывно вложенному в X, причем син-
гулярные значения оператора канонического вложения JV из V в X
стремятся к нулю с полиномиальной скоростью, в то время как син-
гулярные значения {σl}∞l=1 оператора A стремятся к нулю экспонен-
циально. Следуя [6, 14], в данном случае будем предполагать, что x0
принадлежит множеству
Mp,ρ(A) :=
{
x : x = ln−p(A∗A)−1v, ‖v‖ ≤ ρ
}
(1.2)
при некотором неизвестном значении p > 0 и известном ρ > 0, где
операторная функция ln−p(A∗A)−1 определяется спектральным ра-
зложением оператора
A∗A =
∞
∑
k=1
σ2
k(Ψk, ·)Ψk,
то есть
ln−p(A∗A)−1v =
∞
∑
l=1
ln−p σ−2
l (Ψl, v)Ψl.
В дальнейшем без потери общности будем считать, что ‖A‖ ≤ θ ≤
e−1/2, σl ≤ θ ≤ e−1/2, l = 1, 2, . . ..
Пример 1.1. В качестве примера жестко некорректной задачи рас-
смотрим задачу из теории гравитационной градиентометрии спутни-
ка. Если принять поверхность земли и орбиту спутника за сферы с
радиусами r1 < r2, соответственно, Ωri
= {u ∈ R3, |u| = ri}, i = 1, 2,
то одна из задач в этой теории может быть описана уравнением (1.1)
с оператором
Ax(u) :=
1
4πr1
∫
Ωr1
d2
dr2
2
(
r2
2 − r2
1
|u − v|3
)
x(v) dΩr1
(v), u ∈ Ωr2
. (1.3)
В градиентометрии спутника обычно считается, что точное реше-
ние уравнения (1.1) с оператором (1.3) является элементом сфериче-
ского пространства Соболева
Hs :=
{
f ∈ L2(Ωr1
) : ‖f‖2
s =
∞
∑
l=0
2l+1
∑
k=1
(
l +
1
2
)2s
|〈Y (1)
l,k , f〉|2 < ∞
}
С. Г. Солодкий, А. В. Грушевая 449
при некотором положительном индексе s. С другой стороны, так как
для сингулярных значений {σl} оператора (1.3) справедливо следую-
щее соотношение:
lnσ−2
l ≍
(
l +
1
2
)
,
то найдутся положительные константы c2 > c1 > 0 такие, что для
любого элемента f ∈ Hs соболевского пространства выполняются
двусторонние оценки
c1‖f‖s ≤ ‖ lns(A∗A)−1f‖ ≤ c2‖f‖s.
Это, в частности, означает, что каждый элемент из Hs принадлежит
множеству (1.2) при p = s. Более подробную информацию об этом
можно найти в [4, 9].
Поскольку Range(A) предполагается незамкнутым в Y , то не-
прерывной зависимости решения x0 от входных данных нет. Следо-
вательно, численное решение задачи (1.1)–(1.2) требует применения
специальных методов регуляризации.
Отметим, что проблема устойчивого решения жестко некорре-
ктных задач стала интенсивно изучаться, начиная со статьи [6], в
которой для регуляризации был использован метод Тихонова. Далее
жестко некорректные задачи рассматривались в работах [1, 2, 5, 10].
Так, в частности, в [10] регуляризация жестко некорректных задач
(1.1) с решениями из (1.2) проводилась методом Тихонова, а выбор
параметра регуляризации α осуществлялся по принципу невязки Мо-
розова. Этот подход позволил достичь оптимальную по порядку (в
логарифмической шкале) точность O(ln−p 1
δ ) восстановления реше-
ний из указанного множества для любого p : p > p0 > 0. Идея исполь-
зования тихоновской регуляризации в сочетании с принципом невяз-
ки нашла свое дальнейшее применение в работе [12], где было уста-
новлено, что указанный метод обеспечивает оптимальный порядок
точности также в случае возмущенного оператора в (1.1). Позднее
исследования, инициированные в [10], были распространены в [13] на
более широкий класс экспоненциально некорректных задач, для ко-
торого также удалось добиться оптимального порядка точности вос-
становления решений O((ln . . . ln 1
δ )−p).
В настоящей статье, в отличие от всех упомянутых выше работ,
с целью выбора параметра регуляризации α (в комбинации с мето-
дом Тихонова) будет использовано другое правило останова, а имен-
но, принцип равновесия, который для решения некорректных задач
впервые был предложен в [8]. Ниже будет установлена оптимальность
450 Об апостериорном выборе параметра...
предлагаемого подхода на классе жестко некорректных задач (1.1) с
решениями из множества (1.2).
Схема настоящей статьи следующая. Параграф 2 содержит вспо-
могательные утверждения, необходимые для анализа аппроксимаци-
онных свойств предлагаемого нами подхода. В параграфе 3 описыва-
ется принцип равновесия, а в параграфе 4 доказывается оптималь-
ность по порядку комбинации метода Тихонова с принципом равно-
весия в качестве правила выбора параметра регуляризации.
Напомним, что в рамках метода Тихонова регуляризованное ре-
шение xδ
α определяется как решение вариационной задачи
Iα(X) := ‖Ax − yδ‖2 + α‖x‖2 → min .
Поскольку для численной реализации тихоновского метода необходи-
мо выполнять все вычисления с конечномерным приближением An
вместо A, то вариационная задача Iα(X) → min заменяется конечно-
мерным аналогом
Iα,n(x) := ‖Anx − yδ‖2 + α‖x‖2 → min,
где An — некоторое конечномерное приближение к A, rank(An) = n.
Нахождение приближенного решения требует в этом случае решения
линейного операторного уравнения
αx + A∗
nAnx = A∗
nyδ, (1.4)
то есть приближение ищется в виде xδ
α,n = (αI + A∗
nAn)−1A∗
nyδ.
Отметим, что в данной статье мы не ограничиваемся рамками
какой-либо конкретной схемы дискретизации, а накладываем на дис-
кретизированный оператор An только одно, предельно общее, усло-
вие. А именно, конечномерное приближение An должно выбираться
таким образом, чтобы уровень дискретизации и уровень возмущения
правой части были согласованными по порядку, то есть чтобы выпол-
нялось неравенство
‖A − An‖ ≤ δρ−1. (1.5)
Таким образом, на базе предлагаемого нами общего подхода (1.4),
(1.5) к решению жестко некорректных задач, может быть разработа-
на целая совокупность численных алгоритмов в зависимости от выбо-
ра конкретной схемы дискретизации. Примеры проекционных схем
дискретизации некорректных задач, которые позволяют сохранить
оптимальную точность решения при минимально возможном объеме
вычислений, можно найти, например, в работах [3, 11].
С. Г. Солодкий, А. В. Грушевая 451
2. Вспомогательные утверждения
Приведем необходимые определения и факты.
Общее условие истокопредставимости задается следующим обра-
зом:
Aϕ(ρ) := {x ∈ X, x = ϕ(A∗A)v, ‖v‖ ≤ ρ}, (2.1)
где ϕ(t) — произвольная возрастающая функция такая, что ϕ(0) = 0.
При этом ϕ называется индексной функцией.
В дальнейшем нам понадобится следующее вспомогательное ут-
верждение (см. работу [7, предложение 1]).
Лемма 2.1. Предположим, что ϕ — индексная функция и ϕ(t)/t
невозрастающая. Тогда для всех 0 < s, t ≤ e−1 имеем
sup
0<t≤e−1
∣
∣
∣
s
s − t
ϕ(t)
∣
∣
∣
≤ ϕ(s).
Напомним, что в нашем случае индексная функция имеет вид
ϕ(t) = ln−p 1
t . Очевидно, что она удовлетворяет условию леммы 2.1.
Хорошо известно (см., например, [7]), что для любого ограничен-
ного оператора B выполняются соотношения
B(αI + B∗B)−1 = (αI + BB∗)−1B,
‖(αI + B∗B)−1‖ ≤ α−1, ‖(αI + B∗B)−1B∗‖ ≤ 1
2
√
α
,
‖B(αI + B∗B)−1B∗‖ ≤ 1.
(2.2)
Теорема 2.1. Пусть ‖A‖ ≤ e−1/2 и x0 = A−1y ∈ Mp,ρ. Тогда имеет
место следующее неравенство:
‖x0 − xδ
α,n‖ ≤ ρ ln−p 1
α
+
3δ
2
√
α
. (2.3)
Доказательство. Прежде всего воспользуемся правилом треуголь-
ника
‖x0 − xδ
α,n‖ ≤ ‖x0 − xα‖ + ‖xα − xα,n‖ + ‖xα,n − xδ
α,n‖,
где xα = (αI + A∗A)−1A∗y, xα,n = (αI + A∗
nAn)−1A∗
ny.
Оценим первое слагаемое. Используя лемму 2.1, получим
‖x0 − xα‖ = ‖ ln−p(A∗A)−1v − (αI + A∗A)−1A∗A ln−p(A∗A)−1v‖
452 Об апостериорном выборе параметра...
≤ ρ sup
0<λ≤e−1
∣
∣
∣
ln−p 1
λ
− λ
α − λ
ln−p 1
λ
∣
∣
∣
= ρ sup
0<λ≤e−1
∣
∣
∣
α
α − λ
ln−p 1
λ
∣
∣
∣
≤ ρ ln−p 1
α
.
Для оценки второго слагаемого нам потребуются предварительные
выкладки. Принимая во внимание (2.2), имеем
(αI + A∗A)−1A∗y − (αI + A∗
nAn)−1A∗
ny
= A∗(αI + AA∗)−1y − (αI + A∗
nAn)−1A∗
ny
= [A∗ − (αI + A∗
nAn)−1A∗
n(αI + AA∗)](αI + AA∗)−1y
= (αI + A∗
nAn)−1[(αI + A∗
nAn)A∗ − A∗
n(αI + AA∗)](αI + AA∗)−1y
= (αI + A∗
nAn)−1[αA∗ + A∗
nAnA∗ − αA∗
n − A∗
nAA∗](αI + AA∗)−1y
= (αI + A∗
nAn)−1[α(A∗ − A∗
n) + A∗
n(An − A)A∗](αI + AA∗)−1Ax0.
Отсюда, с учетом неравенств (2.2), получаем оценку второго слагае-
мого
‖xα − xα,n‖ = ‖(αI + A∗
nAn)−1[α(A∗ − A∗
n)
+ A∗
n(An − A)A∗](αI + AA∗)−1Ax0‖ ≤ I1 + I2,
где
I1 := α‖(αI + A∗
nAn)−1(A∗ − A∗
n)(αI + AA∗)−1Ax0‖
≤ ‖A − An‖
2
√
α
ρ ≤ δ
2
√
α
,
I2 := ‖(αI + A∗
nAn)−1A∗
n(An − A)A∗(αI + AA∗)Ax0‖ ≤ δ
2
√
α
.
И, наконец,
‖xα,n − xδ
α,n‖ = ‖(αI + A∗
nAn)−1A∗
ny − (αI + A∗
nAn)−1A∗
nyδ‖
≤ ‖(αI + A∗
nAn)−1A∗
n(y − yδ)‖ ≤ δ
2
√
α
.
Суммируя найденные выше оценки, получим искомую оценку.
3. Принцип равновесия
Будем минимизировать правую часть (2.3), выбирая α согласно
принципу равновесия. Суть этого принципа состоит в выборе параме-
тра регуляризации α таким образом, чтобы уравновесить две функ-
ции, задающие оценку погрешности. В рассматриваемом случае эти
С. Г. Солодкий, А. В. Грушевая 453
функции (см. (2.3)) имеют вид
Φ(α) := ρ ln−p 1
α
,
Ψ(α) :=
3δ
2
√
α
.
Поскольку ϕ(t) = ln−p 1
t — монотонно возрастающая функция, то
с ростом α Φ(α) возрастает, а Ψ(α) убывает. Теперь (2.3) можно
переписать в виде
‖x0 − xδ
α,n‖ ≤ Φ(α) + Ψ(α). (3.1)
В силу поведения функций Φ и Ψ (а именно, их монотонности и во-
гнутости) выбор значения параметра регуляризации α = α̂, миними-
зирующего правую часть (3.1), будет уравновешивать величины Φ(α)
и Ψ(α), т.е. Φ(α̂) = Ψ(α̂) и, следовательно,
‖x0 − xδ
α̂,n‖ ≤ 2Φ(α̂). (3.2)
Однако, если функция ϕ неизвестна, то такой априорный выбор те-
оретически наилучшего значения α̂ невозможен. Поэтому в нашей
ситуации необходимо использовать какое-нибудь апостериорное пра-
вило выбора α. С этой целью воспользуемся принципом равновесия.
Чтобы его применить, рассмотрим дискретное множество возможных
значений параметра регуляризации
△N = {αi = (q2)iα0, i = 1, 2, . . . , N}, q > 1. (3.3)
Здесь α0 = nδ2, N : αN ≍ 1.
Используя принцип равновесия, рассмотрим множество
M+(△N ) =
{
αi ∈ △N : ‖xδ
αi,n − xδ
αj ,n‖ ≤ 4Ψ(αj), j = 1, . . . , i
}
и выберем значение параметра регуляризации по правилу
α = α+ := max{α ∈ M+(△N )}. (3.4)
Ниже (см. теорему 4.1) будет показано, что выбор α = α+ дает оцен-
ку погрешности, которая отличается от наилучшей возможной (3.2)
только на множитель 3q.
Чтобы доказать этот факт, рассмотрим вспомогательное множе-
ство
M(△N ) := {αi ∈ △N : Φ(αj) ≤ Ψ(αi), j = 1, . . . , i}
454 Об апостериорном выборе параметра...
и вспомогательную величину
α∗ := max{α ∈ M(△N )}.
Без потери общности будем считать, что
M(△N ) 6= ∅, △N�M(△N ) 6= ∅.
Теперь можно оценить близость точного и приближенного решений
для параметра регуляризации α+.
4. Основные результаты
Теорема 4.1. Пусть параметр регуляризации выбирается согласно
правилу (3.4). Тогда имеет место следующая оценка:
‖x0 − xδ
α+,n‖ ≤ 6qΦ(α̂).
Доказательство. Начнем доказательство с установления неравенст-
ва α∗ ≤ α+. С учетом (3.1), поведения функций Φ(α), Ψ(α) и опре-
деления множества M(△N ), при любом αj < α∗ имеем
‖xδ
α∗,n − xδ
αj ,n‖ ≤ ‖x0 − xδ
α∗,n‖ + ‖x0 − xδ
αj ,n‖
≤ Φ(α∗) + Ψ(α∗) + Φ(αj) + Ψ(αj) ≤ 2Φ(α∗) + Ψ(α∗) + Ψ(αj)
≤ 3Ψ(α∗) + Ψ(αj) ≤ 4Ψ(αj).
Таким образом, включение α∗ ∈ M+(∆N ) доказано. А значит, выпол-
няется α∗ ≤ α+.
Вновь используя (3.1) при α = α∗, а также определения множеств
M+(△N ) и M(△N ), заключаем
‖x0 − xδ
α+,n‖ ≤ ‖x0 − xδ
α∗,n‖ + ‖xδ
α∗,n − xδ
α+,n‖ ≤ 6Ψ(α∗). (4.1)
В силу монотонности Ψ легко видеть, что
Ψ(q2α∗) =
3δ
2
√
q2α∗
=
1
q
· 3δ
2
√
α∗
=
1
q
Ψ(α∗). (4.2)
С другой стороны, очевидно, что α∗ ≤ α̂ ≤ q2α∗. Вместе с (4.1) и (4.2)
это дает
‖x0 − xδ
α+,n‖ ≤ 6qΨ(q2α∗) ≤ 6qΨ(α̂) = 6qΦ(α̂) = 6qρ ln−p 1
α̂
.
Теорема 4.1 доказана.
С. Г. Солодкий, А. В. Грушевая 455
Теорема 4.2. Пусть выполняются условия теоремы 4.1. Тогда спра-
ведлива оценка
‖x0 − xδ
α+,n‖ ≤ c̃ ln−p 1
δ
,
где константа c̃ в общем случае зависит только от q, ρ и p, а при
ρ ≥ 3/2 c̃ = 6qρ. Более того, метод (1.4), (1.5), (3.4) гарантирует
оптимальную по порядку точность приближения на классе задач
(1.1)–(1.2).
Доказательство. Так как при α = α̂ справедливо Φ(α̂) = Ψ(α̂), то
ρ ln−p 1
α̂
=
3δ
2
√
α̂
,
√
α̂ =
3δ
2ρ
lnp 1
α̂
,
α̂ =
9δ2
4ρ2
ln2p 1
α̂
.
Поскольку для любых x > 0 и ν > 0 выполняется lnx < xν , то
α̂ =
9δ2
4ρ2
ln2p 1
α̂
≤ 9δ2
4ρ2
1
α̂
.
Следовательно, для значения параметра регуляризации, теоретиче-
ски минимизирующего погрешность, справедлива оценка
α̂ ≤ 3δ
2ρ
.
В итоге получаем
‖x0 − xδ
α+,n‖ ≤ 6qρ ln−p 2ρ
3δ
= c̃ ln−p 1
δ
.
С другой стороны, в работах [6, 14] было установлено, что погре-
шность восстановления произвольным приближенным методом то-
чного решения x0 из множества (1.2) по yδ не может быть меньше
величины O(ln−p 1
δ ). Отсюда немедленно следует, что наш метод яв-
ляется оптимальным по порядку на классе жестко некорректных за-
дач (1.1)–(1.2).
Замечание 4.1. Проведем сравнительный анализ наших результа-
тов с уже известными. Так, в статье [10] для решения уравнений (1.1)
с x0 ∈ Mp,ρ(A) также использовался стандартный метод Тихонова,
456 Об апостериорном выборе параметра...
при этом в качестве правила останова рассматривался принцип не-
вязки. Полученные в обеих работах ([10] и настоящей) порядки то-
чности совпадают, в то же время наш подход имеет существенное
преимущество. А именно, при построении метода из [10] задействова-
на нижняя граница для возможных значений p: p > p0, где величина
p0 > 0 предполагается известной. То же самое можно сказать и о ме-
тодах из работ [12, 13]. В нашей же работе это ограничение снято и
все результаты справедливы при любых p > 0.
Литература
[1] А. Б. Бакушинский, М. Ю. Кокурин, Итеративные методы для решения
некорректных операторных уравнений с гладкими операторами, М.: Едито-
риал УРСС, 2002.
[2] М. Ю. Кокурин, Н. А. Юсупова, О необходимых и достаточных условиях
медленной сходимости методов решения линейных некорректных задач //
Изв. вузов. Матем., (2003), No. 2, 81–84.
[3] С. В. Переверзев, С. Г. Солодкий, Оптимальная дискретизация некорре-
ктных задач // Укр. мат. журн., 52 (2000), No. 1, 106–121.
[4] W. Freeden, F. Schneider, Regularization wavelets and multiresolution // Inverse
Problems, (1998), No. 14, 225–243.
[5] T. Hohage, Regularization of exponentially ill-posed problems // Numer. Funct.
Anal. Optim., (2000), No. 21, 439–464.
[6] B. A. Mair, Tikhonov regularization for finitely and infinitely smoothing
operators // SIAM J. Math. Anal., 1 (1994), No. 25, 135–147.
[7] P. Mathe, S. Pereverzev, Regularization of some linear ill-posed problems with di-
scretized random noisy data // Mathematics of Computation, 75 (2006), No. 256,
1913–1929.
[8] S. Pereverzev, E. Schock, On the adaptive selection of the parameter in regulari-
zation of ill-posed problems // SIAM J. Numer. Anal., (2005), No. 43, 2060–2076.
[9] R. Rummel, O. L. Colombo, Gravity field determination from satellite gradi-
ometry // Bull. Geod., (1985), No. 59, 233–246.
[10] E. Schock, S. V. Pereverzev, Morozov’s discrepancy principle for Tikhonov
regularization of severely ill-posed problems in finite-dimensional subspaces //
Numer. Funct. Anal. Optim., (2000), No. 21, 901–916.
[11] S. G. Solodky, On a quasi-optimal regularized projection method for solving
operator equations of the first kind // Inverse Problems, 21 (2005), No. 4, 1473—
1485.
[12] S. G. Solodky, A. V. Mosentsova, Morozov’s discrepancy principle for the Ti-
khonov regularization of exponentially ill-posed problems // Comp. Meth. Appl.
Math., 8 (2008), No. 1, 86–98.
[13] S. G. Solodky, A. V. Mosentsova, Unsaturable methods for solving severely ill-
posed problems // Int. J. Comput. Sci. Math., 2 (2009), No. 3, 229–242.
[14] U. Tautenhahn, Optimality for ill-posed problems under general source conditi-
ons // Numer. Funct. Anal. and Optimiz., 1 (1998), No. 19, 377–398.
С. Г. Солодкий, А. В. Грушевая 457
Сведения об авторах
Cергей
Григорьевич
Солодкий,
Анна Викторовна
Грушевая
Институт математики НАН Украины
ул. Терещенковская 3
01601, Киев
Украина
E-Mail: solodky@imath.kiev.ua,
anna_mos@imath.kiev.ua
|