Меню сайта
Элементы теории обобщенных функций в применении к задачам обращения лучевого преобразования
Обобщенная функция это непрерывный линейный функционал на пространстве К всех функций a (x), имеющих производные всех порядков и финитный носитель (свой для каждой из функций α (x)). Любая регулярная интегрируемая функция f(x) задает линейный функционал (f, a ):
. (2.2.1)
Однако на пространстве функций K существуют непрерывные линейные функционалы, которые не могут быть заданы с помощью регулярных интегрируемых функций, наиболее известными примерами таких функционалов являются δ-функция и ее производные. Другим широко известным примером является функционал, основанный на функции (1/x)dx. Функция 1/x x является регулярной, однако она не является интегрируемой. При задании соответствующего функционала интеграл
(2.2.2)
понимается в смысле главного значения:
.
Такое понимание интеграла используется при определении преобразования Гильберта от функции α (x) как свертки с функцией 1/xx.
.
Преобразование Гильберта используется, в частности, в одной из формул обращения преобразования Радона в двумерном пространстве. Эта формула обычно приводится в руководствах по компьютерной рентгеновской томографии. Однако метод свертки и обратного проецирования, часто используемый при построении численных алгоритмов томографической реконструкции, основан на несколько другом виде формулы обращения преобразования Радона. В этом методе по существу используется свертка проекционных данных последовательностью функций сходящихся к 1/xx2 в смысле обобщенных функций.
Линейный функционал, соответствующий функции 1/xx2, или, что то же самое, обобщенная функция 1/xx2 определяется формулой [19]
(2.2.3)
Интеграл в (2.2.3) сходится в обычном смысле для любой функции a (x) из пространства основных, и даже из более широкого класса, функций.
В формулах обращения преобразования Радона используется свертка данных с функцией 1/xx2. Свертка обобщенных функций определяется следующим образом.
Пусть заданы два функционала f и g . Действие функционала f *g являющегося их сверткой, на функцию a из пространства основных задается формулой
(f *g, a )= (fx, gy, a (x + y))). (2.2.4)
Здесь gy означает, что функционал действует на функцию a , как функцию переменной y, а функционал f действует на полученную функцию переменной x. Если функционалы f и g можно задать регулярными функциям, то функционал свертки определенный формулой (2.2.4) можно задать функцией, являющейся сверткой соответствующих функций в обычном смысле.
Здесь следует сделать одно замечание. Даже если функция одной переменной a (t ) имеет финитный носитель, функция двух переменных a (x + y) не является функцией с финитным носителем. Это означает, что существование функционала f *g для конкретных функционалов f и g или необходимо доказывать. Известно, что для существования функционала свертки, достаточно, чтобы один из функционалов имел финитный носитель.
Если рассматривать задачи томографии, то там с функцией 1/xx2 сворачиваются исходные данные, которые регулярны и имеют финитный носитель. Можно показать также, что необходимая свертка выражается формулой:
S(r, j ) = I(r, j ) * (-1/p r2 ) =
(2.2.5)
В реальных ситуациях функция I(r, j ) известна в некотором дискретном множестве точек. Для того, чтобы использовать формулу (2.2.4) нужно построить аппроксимацию функции I(r, j ), такую что интеграл в правой части имеет смысл. Интеграл (2.2.4) заведомо сходится, если функция I(r, j ) принадлежит множеству K, то есть имеет финитный носитель и является бесконечно дифференцируемой.
Однако аппроксимация данных бесконечно дифференцируемой функцией может оказаться громоздкой при построении численных алгоритмов. Кроме того, использование бесконечно дифференцируемых функций может приводить к заглаживанию границ областей с резко отличающимися плотностями. Для сходимости интеграла в (2.2.5) достаточно, чтобы функция I(r, j ) имела в каждой точке конечные односторонние производные первого порядка по переменной r. Это позволяет, в частности, использовать кубические сплайны для построения аппроксимации функции I(r, j ).
Читайте больше >>>
Современные лекарственные препараты растительного происхождения
Использование растений в качестве лекарственных средств
пришло в наш век из глубокой древности и до сих пор играет значительную роль в
арсенале лекарственных средств современной медицины.
Это
обусловле ...
Производные бензодиазепина
I.
дикалий клоразепата (транксен)
b. бензодиазепиновые агонисты БД-Р
B. серотониновые агонисты (транквилизаторы)
a. бус ...
Средства визуализации изображений в компьютерной томографии и цифровых рентгенографических системах
Древняя латинская поговорка гласит:
«Diagnosis cetra - ullae therapiae fundamentum» («Достоверный диагноз -
основа любого лечения»). На протяжении многих веков усилия врачей были
направлены на решен ...
Сестринский процесс в работе участковых медсестер при язвенной болезни
«В кругу жертв язвенной болезни все чаще оказываются люди
молодые, и даже подростки. Результаты профилактики и лечения этого недуга не
удовлетворяют ни врачей, ни больных. Социальная цена заболевания все еще
...