Содержание к диссертации
Введение
1. Задача интегральной геометрии с возмущением в трехмерном слое 25
1.1. Задача интегральной геометрии на семействе параболоидов. Формула обращения 26
1.2. Постановка задачи интегральной геометрии с возмущением. Вспомогательные сведения 31
1.3. Вспомогательные утверждения 35
1.4. Теорема единственности решения задачи с возмущением 1.2 43
2. Задачи интегральной геометрии вольтерровского типа 50
2.1. Постановка задачи интегральной геометрии вольтерровского типа на плоскости. Канонический вид 51
2.2. Единственность решения задачи на плоскости 59
2.3. Задача интегральной геометрии вольтерровского типа в трехмерном слое 71
3. Задача Радона с ограниченным диапазоном углов и возмущением 79
3.1. Постановка задачи Радона с ограниченным диапазоном углов и возмущением на плоскости 80
3.2. Задача интегральной геометрии вольтерровского типа на плоскости и задача Радона с ограниченным диапа зоном углов и возмущением 82
3.3. Единственность решения задачи III.1 94
3.4. Полисингулярные интегральные уравнения с возмуще нием в трехмерном пространстве 97
3.5. Преобразование Радона с ограниченным диапазоном углов и возмущением в трехмерном пространстве 112
4. Восстановление поверхностейпо контурам теней 119
4.1. Постановка задачи и вспомогательные утверждения. Восстановление поверхности типа шапки 120
4.2. Совместное восстановление двух поверхностей 129
Литература 138
- Постановка задачи интегральной геометрии с возмущением. Вспомогательные сведения
- Единственность решения задачи на плоскости
- Задача интегральной геометрии вольтерровского типа на плоскости и задача Радона с ограниченным диапа зоном углов и возмущением
- Совместное восстановление двух поверхностей
Введение к работе
Интегральная геометрия занимается изучением преобразований, ставящих в соответствие функциям на многообразии X их интегралы по подмногообразиям некоторого семейства М [28]. Эта актуальная и бурно развивающаяся область современной математики тесно связана с теорией дифференциальных уравнений и математической физикой, геометрическим анализом и имеет многочисленные приложения при математическом исследовании проблем сейсморазведки, интерпретации данных геофизических и аэрокосмических наблюдений, при решении обратных задач астрофизики и гидроакустики [30, 32, 53, 70, 72, 114]. Разработанные здесь методы являются основой для решения проблем из области медицинской и промышленной томографии [60, 61, 79, 109, 120, 122, 123].
Интегральная геометрия является одним из крупных направлений в теории некорректных задач математической физики и анализа. Основы теории некорректных задач были заложены в работах А. Н. Тихонова, М. М. Лаврентьева, В. К. Иванова (см. [35, 53, 77, 78]). С ее дальнейшим развитием и многочисленными приложениями
можно ознакомиться по монографиям и статьям Ю. Е. Аниконова, В. Я. Арсенина, А. Л. Бухгейма, В. В. Васина, В. И. Дмитриева, В. К. Иванова, С. И. Кабанихина, М. М. Лаврентьева, И. В. Мельниковой, В. А. Морозова, А. И. Прилепко, В. Г. Романова, В. П. Тананы,
A. Н. Тихонова, А. М. Федотова, В. Г. Чередниченко, А. Г. Яголы,
B. Г. Яхно и др. [4, 23, 25, 35, 36, 53, 56, 64, 67, 70, 80, 81, 89, 105].
К задачам восстановления функции по известным интегралам от
нее сводится ряд обратных задач фотометрии, т.е. определения некоторых характеристик источника по результатам оптических измерений [38, 54]. Укажем также на связь задач интегральной геометрии с задачами геометрической томографии — восстановления поверхности по заданному набору ее проекций [112].
Обратные задачи для кинетических уравнений как по постановке, так и по методам исследования близки к задачам интегральной геометрии, в частности к задаче обращения преобразования Радона (см., например, [2, 5, 71]).
Задача восстановления функции по известным интегралам от нее вдоль всевозможных гиперплоскостей в евклидовом пространстве рассматривалась в статье Радона [124]. В этой ставшей классической работе были получены явные формулы обращения для четно- и нечетномерных пространств, разработаны методы решения задачи, а также было указано на возможность рассмотрения различных аналогов этого интегрального преобразования. В вышедшей за год до этого работе Функа [111] изучалась проблема восстановления четной функции на сфере, если известны ее интегралы по большим кругам.
Связь задач интегральной геометрии с дифференциальными уравнениями исследовалась в работах Ф. Иона [115, 116] (см. также [37]). В статье [116] рассматривалось преобразование (впоследствии
названное лучевым), которое функции в трехмерном пространстве ставило в соответствие ее интегралы по всевозможным прямым.
Отправляясь от теории представлений групп Ли, И. М. Гельфанд с соавторами в цикле работ рассмотрел задачи интегральной геометрии на линейных многообразиях. Это стало толчком для изучения интегральной геометрии на широких классах однородных пространств, здесь были получены интересные явные формулы обращения (см. [26, 28-30, 41, 74, 75, 86, 87]).
Особый интерес к задачам интегральной геометрии вызвало обнаружение их важности для решения задач компьютерной томографии. Работы Кормака [106, 107] положили начало применению методов компьютерной томографии в медицине. На основе полученных результатов по задачам интегральной геометрии, в том числе для семейств прямых и fc-мерных плоскостей в n-мерном пространстве были созданы эффективные с вычислительной точки зрения алгоритмы решения соответствующих задач томографии (см. [60, 122]).
В связи с необходимостью изучения разнообразных схем сканирования в компьютерной томографии (веерной, конусной и др.) были получены различные формулы обращения для преобразования Радона, а также лучевого преобразования (см. [60, 87, 109, 122]).
Связь задач интегральной геометрии с многомерными обратными задачами для дифференциальных уравнений с частными производными была выявлена в работе М.М.Лаврентьева и В.Г.Романова [51]. Многообразия, которые возникают при сведении обратных задач к задачам интегральной геометрии, являются либо сечениями характеристических коноидов исходного дифференциального уравнения, либо проекциями бихарактеристик на пространство, ортогональное временной координате [53]. В случае обратных задач для уравнения с
переменными коэффициентами соответствующие задачи интегральной геометрии рассматриваются на геометрически достаточно сложных многообразиях.
Задачами интегральной геометрии волыперровского типа называются задачи, которые могут быть сведены к исследованию операторных уравнений Вольтерра [53].
В.Г.Романов исследовал вопросы единственности и устойчивости решения задач интегральной геометрии в случае, когда многообразия инвариантны относительно группы всех движений, параллельных (п — 1)- мерной гиперплоскости [67,70]. Все весовые функции также предполагались инвариантными относительно данной группы.
Достаточно общие результаты по единственности и устойчивости решения получены в работах А.Ю. Аниконова и А.Л. Бухгей-ма для следующих классов задач интегральной геометрии: решение ищется в классах кусочно-аналитических по части переменных функций; многообразия, по которым ведется интегрирование, и весовые функции аналитическим образом зависят от части переменных (параметров) [4, 23].
Важные результаты по теории задач интегральной геометрии были получены также в работах [1, 3, 21, 24, 88, 90].
Следует отметить работы [1, 57-59]), в которых были получены значимые результаты по единственности и устойчивости задач не вольтерровского типа на плоских кривых в ограниченной области и на геодезических римановой метрики в n-мерном пространстве, а также работу [49] , в которой изучались вопросы единственности и устойчивости "в малом" задач интегральной геометрии не вольтерровского типа в ограниченной области в n-мерном пространстве.
Отметим также группу результатов по интересной и имеющей важные приложения в томографии постановке задачи интегральной геометрии на различных семействах прямых, пересекающих некоторое заданное множество точек в пространстве (кривую) [19, 27, 39, 110, 126].
Значительный прогресс достигнут в изучении классической задачи Радона по семействам гиперплоскостей с постоянной весовой функцией, а также для задачи Радона с полными данными. Однако даже в случае восстановления функции на плоскости по ее интегралам на прямых с весовой функцией, монотонной вдоль прямых (плоская задача Радона с поглощением) еще многие вопросы остаются невыясненными. Это показывает и построенный в работе [104] важный пример неединственности. Некоторые существенные результаты для задачи Радона с поглощением были получены в самое последнее время [91, 117, 118].
Таким образом, все большее значение приобретает исследование задач интегральной геометрии для более широких классов многообразий, по которым ведется интегрирование, и весовых функций (см. также [108], сборник трудов [125]). Не меньшее теоретическое и практическое значение имеют задачи интегральной геометрии с неполными данными [60]. В такой постановке наиболее принципиальным, как указано в [53], является вопрос об однозначности определения искомой функции по интегралам от нее.
Приведем краткое изложение основных результатов работы.
В первой главе изучается задача восстановления функции в трехмерном слое, если известны суммы интегралов от нее по поверхности параболоидов, опирающихся на плоскость х$ = 0, и по объемам,
ограниченным параболоидами и плоскостью х$ = 0, с заданной весовой функцией.
Введем обозначения: х Є R3, Є R3; х = (хих2), f= (6,6); Л Є R2;
S = {x Є R+ : 0 < x3 < h, h < oo}; 5 = {x Є R+ : 0 < 3 < h}.
«
В первом параграфе рассматривается задача интегральной геометрии в слое S на семействе параболоидов {V(x)} с вершинами в точ-
xv.cU? Зи *
V(x) = { : *з - Ь = \х - |2; 0 < 6 < М; р(х) — проекция параболоида V(x) на плоскость #з = 0.
Задача 1.1. Восстановить функцию и, если известны интегралы от нее по поверхностям семейства {V(x)} :
Пусть < , > — скалярное произведение в R2; й(Л, ^з) = х- / ei
M*s) = ^- / ei
Получено аналитическое представление для образа Фурье искомой функции по переменным (х\,Х2). Из вида этой формулы вытекает сильная неустойчивость решения сформулированной задачи.
Теорема 1.1. Пусть функция f(x) известна для всех х Є S. Тогда, решение задачи 1.1 в классе Cq(S) единственно, и имеет место представление
з о Здесь /о — функция Бесселя мнимого аргумента.
В 1.2 приводится постановка задачи интегральной геометрии на семействе параболоидов с возмущением, а также некоторые вспомогательные сведения.
В слое S рассмотрим семейство параболоидов {V(x)} с вершинами в точках х и обозначим: Р(х) — часть слоя S , ограниченная поверхностью параболоида V(x) и плоскостью #з = 0.
Задача 1.2. Решить операторное уравнение относительно функции и:
J u(0d+ / W(x,Qu(t)dt = f(x),
р(х) Р(г)
где W — заданная весовая функция.
В параграфе 1.3 содержатся формулировки и доказательства вспомогательных утверждений — лемм 1.2-1.6, которые также использо-вуются при доказательстве основного результата главы — теоремы единственности решения задачи интегральной геометрии с возмущением в трехмерном слое.
В параграфе 1.4 приводятся формулировка и доказательство теоремы единственности решения задачи интегральной геометрии с возмущением в трехмерном слое при достаточно общих предположениях о весовой функции возмущения.
Теорема 1.2. Пусть f(x) известна для всех х Є S; весовая функция W имеет все непрерывные производные до второго порядка включительно и обращается в 0 вместе со своими производными на поверхности параболоида V{x).
Тогда решение задачи 1.2 в классе Cq(S) единственно.
Во второй главе исследуются задачи интегральной геометрии вольтерровского типа общего вида по кривым на плоскости и поверхностям в трехмерном пространстве.
В параграфе 2.1 приводится постановка плоской задачи. Пусть f = (&,&), х = (xi,x2), R+ = {х = (хъх2) : х2 > 0},
L = {х Є R+ : 0 < х2 < 1,1 < оо}. Рассмотрим семейство кривых {Г(ат)} в полосе
L = {х Є R+ : 0 < х2 < /}.
Кривая семейства Г (х) определяется уравнением
6 = ^(я,&)> Ф(х,хі ~hi) = ^(x,x1 + h2) = 0,
где hi = hi(x), hi = h2(x).
Задача II. 1. Восстановить функцию двух переменных u(x), если известны интегралы от нее с заданной весовой функцией р(х,і) по кривым семейства {Т(х)}:
/р(*,6М0<% = /(*). (і)
Т(х)
Приведены условия, при которых сформулированная задача приводится к каноническому виду, и описана процедура такой редукции.
Во втором параграфе главы приводятся формулировка и доказательство теоремы единственности решения задачи II. 1.
Условие А. Для любого х из L и любого угла а Є [<*i,af2], где — 7г/2 < с*і < 0 < ot< 7г/2, существует единственная кривая семейства Т(х) , проходящая через точку х под углом а с осью ОХ\.
Предположим также, что произвольная кривая семейства симметрична относительно прямой, проходящей через вершину кривой перпендикулярно оси OXi.
При этих предположениях исследуем задачу II. 1.
Теорема 2.2. Пусть функция f(x) — правая часть уравнения (1) — известна для всех х из полосы L; функции р Є C8(L), ф Є C9(L) и удовлетворяют условиям:
p(*,ti)>*i>0; (2)
Ф(х,і) < х2 пряі фхи ф(х,хі) = х2; (3)
-оо < S3 < Ш < 62 < -2, (4)
где 5і, 52, ^з ($з < <^) — постоянные.
Тогда решение задачи ІІ.1 в классе Cq(L) единственно.
Таким образом, условия (2)-(4) вместе с естественным условием отделенности весовой функции от 0, условиями на гладкость весовой функции и функции, описывающей структуру кривой семейства, приводят к единственности решения задачи интегральной геометрии вольтерровского типа общего вида на плоскости.
В параграфе 2.3 рассматривается задача интегральной геометрии вольтерровского типа общего вида на семействе поверхностей в трехмерном слое
L = {х Є R3 : 0 < хг < I, / < оо}.
Пусть {Г(#)} — семейство поверхностей, гладко заполняющих слой L и однозначно параметризующихся с помощью своих вершин х Є L.
Задача II.2. Определить функцию и[х), если известны ее интегралы по поверхностям Т(х) с заданной весовой функцией д(х,)
j g(x,i)u{t,)ds = f(x)
Г(х)
для всех X из слоя L.
Здесь также получена аналогичная теорема единственности решения задачи.
В третьей главе рассматривается задача обращения преобразования Радона с ограниченным диапазоном углов и возмущением на плоскости и в трехмерном пространстве.
В параграфе 3.1 приводятся постановка задачи на плоскости и некоторые вспомогательные сведения. Введем обозначения: х Є R2, у Є R2;
L = {х Є R2 : 0 < хч < I, I < сю};
1(х, о) — прямая, проходящая через точку х Є L под углом а к оси OXh а Є Л, АС [0;тг].
Задача III. 1.
Восстановить функцию v, если при всех х Є L и а Є Л известны суммы интегралов от нее по прямой /(я, а) и по R2 с заданной весовой функцией W
J v(y)ds + J W(x,y,a)v(y)dy = V(x,a). (5)
l(x,a) К2
Первое слагаемое в левой части уравнения (5)
J v(y)ds = V0(x:a)
l(x,a)
представляет собой оператор Радона на плоскости с ограниченным диапазоном углов.
В параграфе 3.2 устанавливается связь между задачей интегральной геометрии вольтерровского типа общего вида на плоскости и плоской задачей Радона с ограниченным диапазоном углов и возмущением.
Пусть
= Кь&),У = Ы,У2), х = (хих2),
R+ = ix — (^1^2) : %2 > 0}. Рассмотрим семейство кривых {Г(#)} в R+, кривая семейства Т(х) определяется уравнением 2 = ^(#>i)- Как и во второй главе, рассмотрим задачу нахождения функции двух переменных и(х), если известны интегралы от нее с заданной весовой функцией р(х,\) на семействе кривых {Г(#)}, т. е. задачу решения операторного уравнения относительно функции и(х):
/p(*,i)««Ki = /(*) (б)
Т(Х)
Вместо того, чтобы параметризовать кривую семейства Т(х) при помощи координат ее вершины (х\,Х2), данную кривую можно определить однозначно с помощью трех параметров: у\,У2 и а- Здесь
(УііУ2) — координаты некоторой точки кривой Г(#), а = а(х,у) — направление касательной к кривой в точке у.
Постановка задачи интегральной геометрии с возмущением. Вспомогательные сведения
В слое S рассмотрим семейство параболоидов { (#)} с вершинами в точках х и обозначим: Р{х) — часть слоя S , ограниченная поверхностью параболоида V(x) и плоскостью х% = 0. . Задача 1.2. Решить операторное уравнение относительно функ- ции и: где W — заданная весовая функция. Предположим, что функция и трижды непрерывно дифференциру ема и финитна с носителем в S . Функция W имеет все непрерывные производные до второго порядка включительно и вместе со своими производными обращается в 0 на поверхности параболоида V{x). Пусть f(x) — правая часть уравнения (1.2.1) — тождественно равна 0. Тогда Здесь через /о и /і обозначены соответственно первое и второе слагаемое в левой части (1.2.1). Рассмотрим функцию где p определено выше, у Є [0; 1]. Из определения функции ф следует, что справедливы следующие соотношения: Здесь и всюду далее Введем в рассмотрение функцию Из (1.2.3) и (1.2.6) с учетом соотношения (1.2.4) следует, что функция q удовлетворяет уравнению Лемма 1.1. Функция qo(x) может быть представлена в следующем виде: где весовая функция W\ непрерывна и выражается через функцию W и ее производные. Доказательство. Проинтегрируем функцию q(x, у) по у от 0 до 1. Из формул (1.2.4)-(1.2.6) следует, что Уравнение (1.2.9) является операторным уравнением Вольтерра 2-го рода (см. [53]) относительно функции и. Его регпение имеет вид: где W$ — непрерывная функция, выражающаяся через функцию W и ее производные. Используя выражение для возмущения /і, т.е. второго слагаемого в левой части уравнения (1.2.1), получаем По нашему предположению, весовая функция W вместе со своими производными обращается в 0 на параболоидах V{x). Из определения функции qo(x) (см. (1.2.8)) и (1.2.10) вытекает, что где функция W\ зависит только от функции W и ее производных. Лемма 1.1 доказана. 1.3. Вспомогательные утверждения В этом параграфе приведем формулировки и доказательства лемм, которые будут использоваться при доказательстве основного результата главы — теоремы единственности решения задачи интегральной геометрии с возмущением в трехмерном слое. Введем следующие обозначения: В пространстве переменных (#3? 2/) рассмотрим полуполосу и условию Здесь до и 9i определяются формулами (1.3.2) и (1.3.3) соответственно.
При фиксированных (АьАг) рассмотрим для функции Q задачу Коши (1.3.4), (1.3.5). Лемма 1.2. Функция q может быть продолжена из R2 х 1 ва всю область R2 х ІЇ как решение задачи Коши (1.3.4), (1.3.5). Доказательство. Решение задачи Коши (1.3.4), (1.3.5) имеет вид где JQ и J\ — функции Бесселя нулевого и первого порядков соответственно. Действительно, непосредственной подстановкой в уравнение (1.3.4) нетрудно убедиться, что функция Q, определенная формулой (1.3.6), удовлетворяет этому уравнению. Далее, легко показать, что задача Коши (1.3.5) для эволюционного уравнения (1.3.4) равномерно корректна, и в силу этого ее решение единственно (см. [42]). Применим к обеим частям уравнения (1.2.7) преобразование Фурье по (х\,Х2). Исходя из известных свойств преобразования Фурье, заключаем, что функция q удовлетворяет уравнению дх3 "" 4 0 Таким образом, Из (1.3.7) следует справедливость утверждения леммы. Пусть параметр Ъ 0. Рассмотрим функции Доказательство. Справедливость утверждения леммы вытекает из (1.3.4) и (1.3.8). Введем следующие обозначения: Рассмотрим первое слагаемое в разложении [L#(A, х3,у; Ь)]2. Поскольку и Є CQ(S), С ПОМОЩЬЮ интегрирования по частям легко по- казать, что после несложных, хотя и довольно громоздких преобразований убеждаемся в справедливости утверждения леммы. Аналогично можно показать, что справедлива Лемма 1.5. Имеет место равенство Как следует из условий теоремы и леммы 1.2, при фиксированных А, #з и Ь функцию можно рассматривать как элемент гильбертова пространства % функций, определенных на полуоси 0 у со, со скалярным произведением: легко видеть, что сопряженный оператор определяется следующим образом: Лемма 1.6. Пусть операторы I и I определяются формулами (1.3.16) и (1.3.17) соответственно. Тогда справедливо соотношение Доказательство: Лемма 1.6 доказывается аналогично доказательству леммы 3 из [48]. Действуя по изложенной там схеме, при фиксированных Л, х$ и Ь подействуем на д сначала интегральным оператором J, а затем сопряженным к нему. Из (1.3.15)-(1.3.17) имеем
Единственность решения задачи на плоскости
Изучается задача восстановления функции в полосе по известным интегралам от нее с заданной весовой функцией вдоль кривых. Кри вые, по которым ведется интегрирование, однозначно параметризу ются координатами своих вершин. Рассматриваемая задача относит ся к вольтерровскому типу и имеет достаточно общий вид. В отли чие от указанных во введении работ [4, 23, 67], на кривые семейства и весовую функцию не накладываются ограничения типа инвариант ности относительно группы всех движений плоскости вдоль фикси рованной прямой или аналитичности (кусочной аналитичности) по каким-либо из переменных или параметров задачи. Кривые семей ства предполагаются строго выпуклыми вверх и имеют конечную «і гладкость, весовая функция положительна, отделена от нуля и также имеет конечную гладкость по всем переменным. Описана процедура приведения задачи интегральной геометрии на плоскости к каноническому виду. Доказана теорема единственности решения задачи общего вида в классе трижды непрерывно дифференцируемых функций с носителем в полосе. Далее рассматривается аналогичная постановка задачи интегральной геометрии в трехмерном слое. Здесь также получена теорема единственности ее решения. 2.1. Постановка задачи интегральной геометрии вольтерровского типа на плоскости. Канонический вид Введем обозначения: = ( ъСг)? х — {x\ix2)\ Рассмотрим семейство кривых {Г(х)} в полосе L, произвольная кривая семейства Т(х) определяется уравнением где hi = hi(x), hi = h2(x). Точка x является вершиной кривой Г(#), т.е. V 6 Г(#) \ фх выполняется условие & = Ф{ХЛ\) х2, Задача II. 1. Восстановить функцию двух переменных и(х), если известны интегралы от нее с заданной весовой функцией р(х,і) на кривых семейства {Г(я)}: Следуя идее М.М. Лаврентьева, опишем процедуру приведения задачи интегральной геометрии II. 1 к каноническому виду. Теорема 2.1. Пусть функция f{x) — правая часть уравнения (2.1.1) — известна для всех х из полосы L; функции р Є С6, ф Є С7 и удовлетворяют условиям: Тогда исходную задачу интегральной геометрии можно привести к следующему виду: и функция WQ(,X2) — непрерывна. Доказательство. Без ограничения общности можно считать, что произвольная кривая семейства {Г(ж)} симметрична относительно прямой, проходящей через вершину кривой х перпендикулярно оси ОХ\, весовая функция р(х,\) также удовлетворяет этому условию. Нетрудно проверить, что все приведенные ниже выкладки остаются верными и без этих ограничений. Уравнение (2.1.1) можно записать в следующем виде: Из условий, которым удовлетворяют функции , ф, вытекает, что верно следующее представление: где tpo, до — положительные функции, отделенные от 0, и имеющие такие же порядки гладкости, как и функции р, ф соответственно.
Введем функцию и применим к обеим частям уравнения (2.1.3) оператор интегрирования по х\ с бесконечными пределами и весом к{х)\ Сходимость интегралов в (2.1.7) вытекает из финитности и, гладкости функций и и д и определения функции к. Далее в подобных случаях очевидные рассуждения приводить не будем. В первом интеграле правой части (2.1.7) сделаем замену Равенство (2.1.8) определяет взаимно-однозначное соответствие между переменными жі и 1. Разрешим его относительно Х\ . Аналогично, сделав во втором интеграле замену 6 = х\ + (р(х,2)- получим: Как следует из (2.1.7), (2.1.9) и (2.1.10), функция может быть представлена в виде Как вытекает из (2.1.4)-(2.1.6) и (2.1.12), функция Gi(,X2) имеет вид: где функция GioK,x2) — гладкая, т.е. Gw(,x2) Є C6(L); Таким образом, мы получили операторное уравнение типа Воль-терра 1-го рода (см. [53]) с весовой функцией, имеющей интегрируемую особенность. Применяя оператор дробного дифференцирования порядка 1/2, рассмотрим функцию: Проделав несложные преобразования, получим: Здесь функция G2, как и функция G\, имеет вид у/х2 - 6 где функция G2Q также достаточно гладкая. Рассмотрим функции а также Как вытекает из (2.1.13)-(2.1.15), функция F4(x2) имеет вид где функция 7з(,#2) — гладкая. Проделав эту операцию еще один раз, получим: и функция Сг4о(, #2) также гладкая. Из условий, наложенных на весовую функцию р и функцию гр, описывающую структуру кривой, вытекает, что функция Gio в окрестности плоскости х2 — 2 непрерывно дифференцируема по переменной Р = Vx2 — 2 до шестого порядка включительно и, следовательно, может быть представлена в виде где Уб(,х2) — непрерывна.
Таким образом, отправляясь от исходной задачи интегральной геометрии и продолжая описанную выше процедуру, можно получить функцию и функция Wo(,x2) — непрерывна. Замечание 2.1. Как следует из приведенного доказательства, при р Є С8, ф Є С9 можно привести исходную задачу интегральной геометрии II. 1 к виду (2.1.16) с весовой функцией, выражающейся по формуле (2.1.17), где функция WQ(,X2) будет уже из класса С3. 2.2. Единственность решения задачи на плоскости В полосе L рассмотрим определенное в первом параграфе главы семейство кривых {Г(х)}. По-прежнему будем считать, что произвольная кривая семейства Г(х) определяется уравнением где h\ = hi(х), hi = h2(x). Точка х является вершиной кривой Т(х). Будем предполагать, что выполняется следующее существенное Условие А. Для любого х из L и любого угла а Є [с і,а2] (— 7г/2 «і 0 с 2 п/2) существует единственная кривая семейства Т(х), проходящая через точку х под углом а с осью OXi. Предположим также, что произвольная кривая семейства симме- трична относительно прямой, проходящей через вершину кривой пер пендикулярно оси OXi. При этих предположениях исследуем задачу II. 1, то есть задачу восстановления функции двух переменных и(х) по известным интегралы от нее с заданной весовой функцией р(#,і) на кривых семейства {Г(ж)}:
Задача интегральной геометрии вольтерровского типа на плоскости и задача Радона с ограниченным диапа зоном углов и возмущением
Применив преобразования Фурье-Лапласа к (2.3.4), аналогично предыдущему параграфу заключаем, что задача интегральной геометрии вольтерровского типа общего вида II.2 и задача интегральной геометрии 1.2 в трехмерном слое на семействе параболидов с возмущением переходят друг в друга путем интегральных преобразований с помощью аналитических функций. Отсюда следует, что весовая функция во втором слагаемом из левой части последнего уравнения достаточно гладкая и вместе со своими производными до второго порядка включительно обращается в 0 на параболоидах V(x0). Таким образом, мы попадаем в условия теоремы 1.2, откуда следует единственность решения исходной задачи интегральной геометрии Замечание 2.3. Как указано в параграфе 1.1, сильная неустойчивость задачи 1.1 вытекает из явного вида формулы обращения. Соответствующая задача интегральной геометрии с возмущением будет, вообще говоря, также сильно неустойчивой. То же самое следует из полученных результатов для задачи интегральной геометрии вольтерровского типа II.2. Тем не менее не только в полупространстве R+, но и в трехмерном слое можно, по-видимому, подобрать разрывные весовые функции специального вида таким образом, чтобы решение задачи было устойчивым в пространствах, в определении нормы которых участвует конечное число производных. Изучается задача обращения преобразования Радона на плоскости с ограниченным диапазоном углов и возмущением. Доказана теорема единственности ее репіения. Подобные операторные уравнения при других ограничениях на ядро возмущения рассматривались в работах [45, 49]. Рассматривается также плоская задача интегральной геометрии вольтерровского типа общего вида. Получены условия, при которых вопрос о единственности решения задачи интегральной геометрии "в малом" сводится к исследованию единственности решения операторного уравнения типа задачи Радона с возмущением. Иследована аналогичная постановки задачи Радона с ограниченным диапазоном углов и возмущением в трехмерном пространстве. Получена теорема единственности ее решения. При этом возникают особые многомерные сингулярные уравнения с возмущением, названные в работе по аналогии с бисингулярными интегральными уравнениями (см. [46]) полисингулярными.
Отметим, что, в отличие от плоского случая, задача решения этих уравнений может быть как сильно, так и слабо неустойчивыми. Доказаны теоремы единственности, а также существования решения подобных уравнений. 3.1. Постановка задачи Радона с ограниченным диапазоном углов и возмущением на плоскости Введем обозначения: 1{х,а) — прямая, проходящая через точку х Є L под углом а к оси ОХъ а Є Л, Л С [0; тг]. Задача ІІІ.1. В полосе L рассмотрим операторное уравнение относительно функции v Отметим, что подобные операторные уравнения рассматривались в работах М. М. Лаврентьева и А. Л. Бухгейма [45, 49]. Нами получена теорема единственности решения уравнения (3.1.1) при других ограничениях, накладываемых на весовую функцию возмущения. Легко заметить, что первое слагаемое в левой части уравнения (3.1.1) удовлетворяет диференциальному уравнению 1-го порядка 3.2. Задача интегральной геометрии вольтерровского типа на плоскости и задача Радона с ограниченным диапазоном углов и возмущением Пусть Рассмотрим семейство кривых {Г(:г)} в R/j_, кривая семейства Г(х) определяется уравнением Рассмотрим задачу восстановления функции двух переменных и(х), если известны интегралы от нее с заданной весовой функцией р{х,і) по кривым семейства {Г(а:)}, т. е. задачу решения операторного уравнения относительно функции и(х): га — единичный вектор нормали к прямой 1(х,а) в точке х, Н(х, а) — полуплоскость, ограниченная прямой /(#, а) и содержащая вектор т. Теорема 3.1. Пусть функция f(x) — правая часть уравнения (3.2.1) — известна для всех х яз полосы Lh , функции риф дважды непрерывно дифференцируемы и удовлетворяют условиям: Тогда: 1) Уравнение (3.2.1) может быть приведено к виду где функция В непрерывно дифференцируема по всем переменным; 2) Единственность решения уравнения (3.2.5) в классе Cg(Lft) влечет за собой единственность решения уравнения (3.2.1) в том же классе. Доказательство: Как вытекает из постановки задачи и (3.2.2)-(3.2.4), кривые семейства удовлетворяют следующим условиям: a) они однозначно параметризуются с помощью своих вершин х; b) гладко заполняют полосу L/,; c) любая кривая семейства {Г(х)} может быть продолжена в Rl = {х — (х\,х2) : Х2 0} с сохранением гладкости; d) для любой кривой семейства {Г(а:)} и любой ее точки угол, который образует касательная к кривой Г(х) в точке с положительным направлением оси 0Х\, принадлежит множеству Л — {а : -7г/2 а\ а а2 7г/2 }.
Очевидно также, что ф{х,і) х2, & фхх и найдутся гладкие функции q\ = q\(x),q2 = q2(x) такие, что ф(х, xi - qi) = ф(х, Xi + q2) = 0. Тогда уравнение (3.2.1) можно представить в виде В уравнениии (3.2.6) х\—фі(х, 2) есть обратная функция к функции Ф{хіі) по переменной і при фиксированных (хі,х2)1 если і Є [#і — qi,xi). Аналогично 1 + 2( , 2) есть обратная функция к функции ift(x,{) при фиксированных (#1,2 ), если і Є {хі,х\ + 92]- Ясно, что фі(х,х2) = ф2(х,х2) 0. Как следует из отмеченных выше условий, где А; = 1,2; (ж, ) и 0 (#,2) — гладкие функции, отделенные от нуля. Рассмотрим угол с , который касательная к кривой Т(х) в некоторой ее точке у образует с положительным направлением оси ОХ\. Как отмечено выше, Легко видеть, что При фиксированных (xi,#2) каждому значению у2 из промежутка [0; Х2) отвечают два значения аг, соответствующие левой и правой частям кривой Г(х) . Далее подробно рассмотрим только одно из этих значений, отвечающее, например, левой части кривой. Это равносильно дополнительному условию а Є (0;о2]. Для or, отвечающего правой части кривой Т(х), будут кратко отмечаться необходимые изменения в ходе доказательства.
Совместное восстановление двух поверхностей
В настоящей главе изучается задача восстановления поверхностей, опирающихся на заданную плоскость, по контурам их теней. Эта задача, как отмечено во введении, является задачей геометрической томографии [31, 92, 112, 113]. Последняя же, как отмечено в [112], имеет многочисленные связи с интегральной геометрией, в том числе с задачей Радона. Вводятся понятия оптической проекции точки поверхности, а также контуров теней, отвечающих перпендикулярным и оптическим проекциям. Приведены необходимые и достаточные условия однозначного восстановления плоской кривой по полному набору оптических и перпендикулярных проекций ее точек. Получен результат по единственности решения задачи в случае гладкой поверхности типа шапки. Доказана теорема единственности решения задачи одновременного восстановления двух непересекающихся кусочно-гладких поверхностей типа шапки по полному набору их теней. 4.1. Постановка задачи и вспомогательные утверждения. Восстановление поверхности типа шапки Пусть векторы а(а) лежат в плоскости x i = 0 и образуют угол а с положительным направлением оси ОХ і, а Є [—тг; 0]. Зададим систему ориентированных прямых {/"} следующим образом: 1) при любом г = (г\,Г2) направление прямой /" совпадает с направлением вектора а (а); 2) при любом а Є [—7г; 0] прямые {/"} заполняют все пространство 3) при а Є (—7г;0) прямая /" проходит через точку (гі,Г2,0), а при а = — 7г или а = 0 — через точку (0, п, гг). Множество {/г(/?)} С {/"} состоит из всех прямых вида / , где /З Є [—7г; 0] фиксировано, г Є Поверхностью типа шапки назовем строго выпуклую поверхность, которая пересекается с любой прямой, перпендикулярной плоскости #з = 0, не более чем в одной точке. Рассмотрим гладкую поверхность Н типа шапки, лежащую в R-+ = {х Є R3 : хз 0} и опирающуюся на плоскость х$ = 0. Подчеркнем, что всюду в работе под гладкостью будем понимать принадлежность классу С1. Обозначим: — некоторая точка поверхности Н, Т() — касательная плоскость к поверхности Н в этой точке, /() — плоскость, проходящая через точку параллельно плоскости Введем направления "+" и "-" на прямой /() = Т() П /() и обозначим через Т (Н) множество {/_;/+}.
Пусть Тогда прямую 1@ назовем (і-касательной к поверхности Я в точке . Угол (3 и индекс q 6 Q С R2 однозначно определяются заданием точки Я;(3 = f3(H,0;q = q(H,0-ОбраТНО, Существуют ОГі,ОГ2, — 7Г «і —7г/2 Qf2 R2, Q = ф(Я, /3) такие, что произвольная точка Є Я (за исключением точек, являющихся вершинами кривых вида ЯП/()) однозначно определяется заданием (3 и q: Легко видеть, что тогда /Зі = — тг, /32 = 0, точка является вершиной кривой ЯП/() и в этой точке существуют — тг-касательная и 0-касательная к поверхности Я. Наконец, пусть (3 6 (а а ). Тогда Итак, пусть прямая 1& является /3-касательной к поверхности Я в точке . Тогда точку назовем точкой касания прямой № с поверхностью Н; точку пересечения прямой № с плоскостью #з = О — оптической проекцией точки . Предположим, что при некотором фиксированном индекс q принимает все возможные значения из множества Кривую Sp{H) на поверхности Н, состоящую из точек касания прямых множества {lq(@)} с этой поверхностью, назовем {3-линией тени на поверхности Н. Очевидно, что при / — --71-,/ = 0 (3\- и (3%-линии тени на поверхности Н совпадают; при любых / ф / из (—7г; OL{\ U [с 2; 0) /Зі- и /#2- линии тени на поверхности Н имеют только две общие точки hi и 2, лежащие на плоскости х$ = 0. Теперь можно ввести определения контуров теней гладкой поверхности типа шапки. Рр(Н)-контур тени — кривая, образованная перпендикулярными проекциями точек кривой Sp{H) на плоскость х$ = 0. Lp(H)-контур тени — кривая, образованная оптическими проекциями точек кривой Sp(H). Задача IY.1. Восстановить поверхность Н по контурам теней. В плоском случае естественно говорить о восстановлении кривой по границам тени, то есть по оптическим и перпендикулярным проекциям точек этой кривой. Приведем вначале следующий пример. Пусть заданы функции Найти кривую Р, границы тени которой определяются функциями а(а) и 6(a). Покажем, что кривая Р восстанавливается единственным образом и представляет собой часть параболы х — с — х\, лежащую в полуплоскости Х2 0. Для этого рассмотрим семейство прямых где параметры а и Ь связаны соотношением Кривая Р является огибающей семейства прямых (4.1.3), (4.1.4) и совпадает с его дискриминантой (см. [34]).
Она может быть найдена как решение системы (4.1.5) которое существует и единственно. Исключая из (4.1.5) параметры а и Ь и учитывая условия на а, получаем: х2 = с х1; х\ [—у/с; \/с\. Легко проверить, что для любого а(а) и Ь(а), определенные в (4.1.1) и (4.1.2), являются соответственно перпендикулярной и оптической проекциями точки этой кривой, отвечающей углу а, то есть точки, ориентированная касательная к которой образует угол а; с положительным направлением оси ОХ\. Теперь рассмотрим более общую постановку. Кривой типа шапки на плоскости OX\X i назовем строго выпуклую кривую, которая пересекается с любой прямой, перпендикулярной оси ОХ\ , не более чем в одной точке. Обозначим через p(ct) перпендикулярную, а через ф(а) оптическую проекцию точки кривой, в которой ориентированная касательная к кривой образует угол а с положительным направлением оси ОХ і; Предложение 4.1. Пусть Г Є С1 — кривая типа шапки на плоскости, опирающаяся на ось ОХ\. Тогда для любого a Є Л выполняется следующее условие согласования Доказательство: Запишем уравнение касательной 1(a) к кривой Г, проходящей через точку (ф(а);0): Дифференцируя (4.1.7) поаи учитывая, что получившееся выражение имеет место для всех (хі;хг) Є /(»), приходим к требуемому соотношению между р(а) и ф(а). Предложение доказано.