Численно-аналитическое моделирование восстановления оптических сигналов и изображений тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Исаев, Юсуп Ниязбекович

• Исаев, Юсуп Ниязбекович
• доктор физико-математических наук
• 2004

Оглавление диссертации доктор физико-математических наук Исаев, Юсуп Ниязбекович

Глава 1. Методы моделирования полиномиального восстановления случайной фазы на основе измерений локальных наклонов волнового фронта.

§1.1 Преобразование Радона в задаче фазового оптического контроля.

§1.2 Представление аберраций в базисе Карунена-Лоэва-Обухова.

§1.3 Погрешность алгоритмов и численный эксперимент.

§ 1.4 Синтез базиса для оптимального разложения фазы оптической волны с учётом внешнего масштаба турбулентности.

§ 1.5 Аналитический метод определения векторных ортогональных полиномов для градиента фазы при произвольной геометрии приёмной апертуры.

Результаты главы -1.

Глава-2. Простые и быстрые базисы для восстановления случайных волновых полей.

§ 2.1 Представление функций Карунена-Лоэва-Обухова в базисах Уолша и Хаара.

§2.2 Численное моделирование восстановления случайных волновых полей с использованием быстрых преобразований.

§2.3 Критерии эффективности адаптивных оптических систем при различных базисах разложения фазы случайной волны.

Результаты главы - 2.

Глава 3 Представление 1-D оптических сигналов в ортогональных вэйвлет базисах.

§3.1 Синтез ортогональных вэйвлет базисов.

§3.2 Разложение сигналов в ортогональный вэйвлет базис.

§3.3 Графическое представление скейлинговых функций и вэйвлетов.

§3.4 Примеры разложения сигналов в D1- вэйвлет базис.

Результаты главы - 3.

Глава 4 Биортогональный и комплексный вэйвлет базисы.

§4.1 Скейлинговые функции биортогоналыюго вэйвлет базиса.

§4.2 Критерии синтеза вэйвлет базиса.

§4.3 Комплексные вэйвлет базисы.

§4.4 20-Вэйвлет базисы.

Результаты главы - 4.

Глава 5 Представление дифференциальных и обратных им операторов в вэйвлет базисах. Вэйвлет - пакеты.

§5.1 Выбор базиса для представления функций с особенностями.

§5.2 Представление дифференциальных и обратных им операторов в базисе скейлинговых функций.

§5.3 Вэйвлет-пакеты.

Результаты главы -5.

Глава 6. Восстановление распределения интенсивности лазерного излучения по температуре поверхности мишени. Неоднородный поток.

§6.1 Постановка задачи. Выводы основных соотношений на основе метода инвариантного погружения.

§6.2 Восстановление параметров лазерного пучка по температуре нагретой поверхности.

§6.3 Регуляризация решения многомерной обратной задачи теплопроводности. Численные результаты.

§6.4 Восстановление интенсивности лазерного пучка по зашумлённому температурному полю мишени.

§6.5 Обработка модельного и лабораторного экспериментов.

Результаты главы -6.

Введение диссертации (часть автореферата) на тему «Численно-аналитическое моделирование восстановления оптических сигналов и изображений»

АКТУАЛЬНОСТЬ РАБОТЫ

Уникальные свойства лазерного излучения обусловили широкое применение лазеров в различных областях науки и техники. В их числе нелинейная оптика, технологические процессы, вычислительная техника, медицина, и т.п. В особую проблему выливаются задачи концентрации лазерной энергии на значительные расстояния в атмосфере, задачи использования оптического излучения для зондирования атмосферных параметров, дальнометрии, локации, навигации и связи. Применение оптического излучения для этих целей сталкивается с необходимостью учета влияния атмосферы на параметры волны. Для подавления атмосферных искажений необходимо использовать адаптивное управление волновым фронтом излучения, а следовательно, необходимо измерять интенсивность и фазу фронта волны в цепи обратной связи адаптивных оптических систем (АОС) в реальном масштабе времени[1-10]. Интенсивность и фаза лазерного излучения, прошедшего атмосферу со случайными неоднородностями используются для извлечения информации об атмосфере или передачи информации через неё. В результате взаимодействия с турбулентной средой распространения фаза и интенсивность становятся сложными топологическими объектами. Происходит преобразование гладкого рельефа волнового фронта в прерывистую изломанную структуру со степенными особенностями и сингулярностями. Для удобства математического моделирования (АОС) и упрощения анализа эти сложные математические объекты, условно будем называть их сигналами, обычно представляются в виде сумм ортогональных составляющих бесчисленным количеством способов. И поскольку каждый раз система ортогональных функций, применяемая для разложения бывает известна, то интенсивность и фаза оптического излучения полностью определяются наборами весовых чисел для этих функций.

Такие наборы чисел - это спектры оптических сигналов.

Спектр - это единственно возможная форма аналитического выражения сигнала в рамках линейной теории. И весь вопрос сводится лишь к выбору наиболее подходящей базисной системы функций удобной для решения той или иной практической задачи.

Насколько адекватным реальности окажется модель анализируемого сигнала зависит от удачи выбора базиса. При решении (моделировании) задач восстановления оптического поля возникает необходимость хранить большие массивы данных изображения, так как оптические сигналы прошедшие атмосферу являются случайными, и следовательно необходимо хранить ансамбль реализаций. Возникает необходимость устранять шумы оптического изображения, сокращать время обработки данных и тем самым приблизиться к реальным временным масштабам физического процесса. Отслеживать эволюцию частоты сигнала, вызванную различными временными и пространственными масштабами неоднородности.

СОСТОЯНИЕ ПРОБЛЕМЫ

Существующие методы восстановления мощного лазерного излучения заключаются в восстановлении распределения интенсивности по болометрическим измерениям, по результатам измерений с помощью матриц пироэлектрических, пиромагнитных детекторов излучения, матриц фотоприёмников [11-17]. Однако указанные измерители не всегда удовлетворяют требованиям по пространственному разрешению и сложны в конструктивном исполнении. Измерения распределения интенсивности мощного лазерного излучения становятся особенно сложными в натурных условиях на наклонных трассах распространения. Возможный способ определения энергетических параметров падающего на мишень излучения в этом случае заключается в восстановлении распределения интенсивности в пучке из результатов дистанционного определения температурного поля нагреваемой им поверхности мишени или измерения рассеянного аэрозолем или мишенью лазерного излучения [18]. Перспективным для таких измерений является использование тепловизоров [19-22]. Способных обеспечить значительное пространственное разрешение и временное разрешение при измерении. Однако вопросы математического моделирования и алгоритмы функционирования подобных методов реконструкции параметров лазерных пучков до появления работ автора не рассматривались.

Измерения распределения фазы оптической волны в силу квадратичности оптических детекторов, реагирующих на интенсивность излучения, всегда требуют организации специальных приёмов для извлечения информации о фазе из измерений её интенсивности, то есть решение «фазовой проблемы» в том или ином варианте её постановки [23-27]. Датчики волнового фронта измеряют, как правило, локальные наклоны волнового фронта в дискретном множестве точек [23-25, 28-30], после чего эти данные преобразуются в значения фазы. Чаще всего при моделировании фазовые возмущения раскладываются в ряд по элементарным компонентам, в качестве которых принято использовать круговые полиномы Цернике, описывающие классические аберрации оптических систем. При этом коэффициенты разложения отвечающие за относительный вклад отдельных мод (наклонов, кривизны, дисторсии .) в фазовые искажения. Существующие алгоритмы определения полиномиального разложения фазы основаны на представлении производных - локальных наклонов волнового фронта -конечными разностями с последующим решением системы линейных алгебраических уравнений[23-25]. Такой подход в известной мере затрудняет оценку качества восстановления из-за присутствия шумов в измерениях наклонов фазы, накладывает ограничения на размер, форму и положения датчиков.

Указанные недостатки могут быть устранены на базе аналитических соотношений связывающих фазу волны с её производными по координатам [30, 31]. Однако соответствующие алгоритмы не были разработаны. Не была решена задача построения статистически оптимальных модовых разложений фазы для устройств когерентной оптики, позволяющих адекватно моделировать как восстановление волнового фронта, так и функционирование устройств когерентной оптической адаптивной техники, предназначенных для функционирования в турбулентной атмосфере. Одним из возможных подходов к решению подобных задач моделирования может рассматриваться использование разложения Карунена-Лоева-Обухова (КЛО), если речь идёт о низкочастотных аберрациях. Если в атмосфере присутствуют высокочастотные аберрации, вызывающие появления сингулярностей и скачков в волновом фронте, в этом случае необходимо использовать вэйвлет преобразования. Теория вэйвлет преобразований за последние годы получила сильное развитие, главным образом усилиями зарубежных математиков [32-37]. Тем не менее, во первых, для большой армии отечественных исследователей в этой области до сих пор присутствует определённый вакуум. Во вторых, большая часть зарубежных работ слишком поверхностна и при большой ширине обхвата даёт лишь наброски узловых моментов. В третьих, в большей части работ по вэйвлетам гарантируется возможность разложения исследуемой функции в один из рассматриваемых базисов без объяснений критерия его выбора, то есть, нет ответа на вопрос какой базис нужно выбирать для моделирования конкретной задачи. Оптическое излучение - это, как правило, различного рода векторные и скалярные поля, которые при моделировании представляются в виде изображений, поэтому в зависимости от того хотите вы выделить скрытые периодичности, контур изображения, повысить его контрастность, выделить сингулярности или осуществить его сжатие, нужно выбирать различные вэйвлет базисы.

Таким образом, целью диссертационной работы является разработка методов численно-аналитического моделирования восстановления оптических сигналов и изображений.

ЦЕЛЬ РАБОТЫ

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

Разработка математических методов моделирования и алгоритмов восстановления распределения фазы и коэффициентов его разложения в оптимальный базис по измерениям локальных наклонов волнового фронта.

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

Первая глава посвящена разработке методов моделирования и алгоритмов восстановления фазы на базе аналитического интегрального представления, связывающего фазу с ее производной [30, 31].

Для описания фазовых искажений световых полей на атмосферных трассах обычно используют их разложение на моды в системе базисных функций. Среди наиболее известных методов решения этой задачи можно выделить метод разложения фазы по полиномам Цернике [23-25, 28,29,38]. Использование аналитического соотношения [30,31] позволяет, как показано в диссертации, восстанавливать коэффициенты разложения фазы по полиномам Цернике без предварительного этапа ее восстановления. При этом положение и форма датчиков могут быть произвольными. Такой подход позволяет быстро оценивать аберрации контролируемой оптической системы.

Самодостаточность представления Цернике нарушается, если речь идет об аберрациях, обусловленных турбулентностью атмосферы или вносимых самой адаптивной системой. В таких случаях более эффективным оказывается универсальное разложение, отвечающее ряду условий оптимальности, разложение Карунена-Лоэва-Обухова. Оно характеризуется следующими свойствами:

1. Минимальность среднеквадратической ошибки при удержании заданного числа членов в бесконечном ряде разложения фазы.

2. Наибольшее (по сравнению с любым другим разложением) количество информации о представляемой усеченным рядом функции фазы, какое бы число членов рядов не удерживалось.

3. Некоррелированность коэффициентов разложения, что упрощает дальнейшее использование результатов и их анализ.

В то же время до работ автора существовало устоявшееся мнение о невозможности реализации данного разложения потому, что "собственные функции разложения характеристик искаженного поля имеют весьма сложную структуру, а практическая реализация их в виде корректирующих устройств с переменным базисом функций оказывается затруднительным"[38]. Отсутствие достаточной теоретической проработки указанной проблемы приводило к тому, что на практике ограничивались близким к оптимальному разложением Цернике. В связи с этим, в третьей главе диссертационной работы значительное внимание было уделено созданию метода получения аналитического представления базиса Карунена-Лоэва-Обухова для разложения фазы и установлению связи этого представления с разложением в базисе Цернике. В главе описываются алгоритмы и результаты построения оптимального модового разложения случайной фазы оптической волны в турбулентной атмосфере с учётом внешнего масштаба турбулентности. Приводится аналитический вид основных мод Карунена-Лоэва-Обухова в представлении базиса Цернике.

Во второй главе рассматривается связь статически оптимального базиса Карунена-Лоева-Обухова (КЛО) с преобразованиями Хаара и Уолша, которые характеризуются простотой технической реализации в виде корректирующих устройств адаптивной оптической системы. Эти функции принимают только два значения -1 и 1, в силу чего удаётся избежать операции умножения при реализации алгоритмов на ЭВМ. Время выполнения преобразования Уолша в 10 раз быстрее быстрого преобразования Фурье того же массива данных[39,40], а преобразование Хаара во столько же раз быстрее преобразования Уолша[40]. Таким образом, связь между быстрыми базисами и оптимальными функциями KJIO позволяют оптимизировать работу АОС, и сокращает время отклика адаптивного контура. В этой же главе рассчитываются критерии качества функционирования АОС в различных базисах и приводится явный вид матриц преобразования от одного базиса к другому.

В третьей и четвёртой главах описывается синтез оптимальных вэйвлет базисов для моделирования сигналов, прошедших атмосферу в условиях сильной турбулентности. Из большого класса вэйвлет-базисов рассматривается синтез вэйвлетов, подчиняющихся условию кратно-масштабного анализа [41-45]. Так как именно такие вэйвлеты сочетают в себе свойства иерархической структуры турбулентности, её фрактальности, они подчиняются свойствам быстрых преобразований на порядок превышающие скорость вычислений БПФ, они близки к статистически оптимальному базису KJIO [42,46,47], и, наконец, они позволяют представлять степенные особенности и сингулярности оптических пучков претерпевших сильное искажение в атмосфере.

В третьей главе рассматривается синтез ортогональных, симметричных вэйвлет-базисов, вопреки установившемуся мнению о невозможности их существования [43,48]. Их симметричность и ортогональность сокращает быстродействие численных алгоритмов. На основе энтропийного критерия конструируется базис, наилучшим образом отображающий форму анализируемого сигнала.

Одной из важных проблем видения изображений, наблюдаемых в рассеивающих средах и описываемых в терминах уравнения сверки, является задача восстановления замутненных изображений. Класс задач, связанных с этой проблемой, например, в спутниковом мониторинге подстилающей поверхности Земли, необозримо велик. Сложность решения этой задачи определяется тем, что функция размытия точки (ФРТ) на момент регистрации изображений, как правило, неизвестна, кроме того, необходимо учитывать некорректность операции обращения свертки. Среди многообразия реальных атмосферно - оптических ситуаций нередки случаи, когда мы наблюдаем достаточно четкое, хотя и слабоконтрастное изображение на фоне общего «размытия». Это соответствует такому виду ФРТ, которая содержит четко выраженную дельта - составляющую и «размытое основание». Если построить преобразование, «усиливающее» роль высокочастотных составляющих наблюдаемого размытого изображения и подавить фоновую составляющую, то получим восстановленное высококонтрастное изображение. Технология использования вейвлет - преобразований для решения этой задачи заключается в следующем. Синтезируется специальное вейвлет -преобразование, позволяющее выделять градиенты изображения. Поэтому в четвёртой главе рассматривается синтез биортогональных вэйвлетов. Преимущество биортогональных вэйвлетов перед ортогональными в их повышенной гладкости. В возможности синтезировать вэйвлет наперёд заданными числом нулевых моментов, что существенно повышает оптимальность разложения [42,47,49]. Разложения по биортогональным вэйвлетам позволяет визуализовать скрытые периодичности сигналов и гасить низкочастотные тренды в изображениях для повышения контрастности. В этой же главе приводятся комплексные вэйвлеты, позволяющие выделять контуры изображения. При разложении спутниковых изображений в комплексный вэйвлет-базис происходит «перекачка» сингулярностей изображения в мнимую часть разложения. Сингулярности - это резкие изменения градиента интенсивности изображения [47,49]. После фильтрации мнимой части изображения от шума получается ярко выраженная картина контуров.

При моделировании восстановления фазы и интенсивности оптического поля приходится использовать дифференциальные и обратные им операторы. Необходимо решать задачи оптимального гашения шумов изображения. Как было упомянуто выше, вэйвлеты не имеют аналитического представления. Поэтому в пятой главе приводятся представления дифференциальных и обратных им операторов в вэйвлет - базисах, а для оптимального гашения шумов описывается алгоритм синтеза адаптивных вэйвлет пакетов, которые наилучшим образом «подстраиваются» под шумы [47,50].

В шестой главе описываются методы моделирования и алгоритмы восстановления интенсивности волны по температурному полю поверхности мишени.

Полагая тепловые потери незначительными, а коэффициент отражения R известным, можно считать, что интенсивность падающего пучка с точностью до множителя (1 - R) равна тепловому потоку на поверхности мишени, за счет которого формируется температурное поле в образце. Для восстановления потока необходимо решить обратную задачу теплопроводности (задачу пересчета граничных условий).

Значительное число работ [15-17, 31, 32] посвящено решению прямых задач, где определяется температура внутри облучаемой мишени в зависимости от ее конфигурации и интенсивности падающего на нее теплового потока. При этом достаточно решить дифференциальное уравнение теплопроводности с граничными и краевыми условиями.

Задача восстановления интенсивности падающего излучения требует решения задачи теплопроводности в ее обратной постановке. В этом случае решения имеют вид интегральных соотношений, связывающих температуру на поверхности мишени с падающим на нее тепловым потоком [33-35]. Эти интегральные соотношения имеют особенности и требуют специальных методов аналитической регуляризации, при численной реализации которых возникают большие трудности. Решение одномерной задачи (падающий поток должен быть однородным) для бесконечной в поперечном направлении мишени приводится в работах [5, 6, 51, 52]. Однако для реальных ситуаций эти задачи не представляют большого интереса. Поэтому автором разработан метод решения обратной тепловой задачи в многомерной постановке [51-55], без ограничений на продольные размеры мишени и граничные режимы. Такой подход позволяет учесть неоднородность распределения интенсивности падающего излучения и предоставляет свободу выбора размеров мишени и граничных условий на ее поверхности при реализации такого способа измерений.

Процесс переноса тепла в твердых телах характеризуется значительным сглаживанием особенностей граничных функций по мере удаления рассматриваемой точки внутрь тела от теплообменной поверхности, темп изменения температуры в удаленных внутренних точках может оказаться существенно ниже темпа изменения температуры внешней поверхности. "Такая физика распространения тепла приводит к известной патологической особенности обратных задач, заключающейся в отсутствии непрерывной зависимости результатов от входных температурных данных (нарушения условия корректности по Адамару)"[56].

Классический метод решения задач восстановления распределения интенсивности лазерного пучка по температурному полю на поверхности произвольной мишени, заключается в следующем: решается прямая задача теплопроводности с краевыми условиями (определяется температура в образце), после решения прямой задачи решается задача Коши для приграничной к мишени области (определяется поток). Метод инвариантного погружения [57-61], заключающийся в введении дополнительной пространственной переменной (параметром погружения является продольный размер мишени), позволяет связать краевую задачу с задачей Коши, что освобождает от необходимости решать прямую краевую задачу. В заключительном уравнении погружения имеется прямая связь потока с температурой на поверхности мишени и содержится параметр погружения, варьируя который могут быть получены решения для произвольных размеров мишени (бесконечно-тонкая, ограниченная и полуограниченная). Результатом решения данной задачи является обобщенное интегральное соотношение, связывающее поток с температурой на поверхности мишени. Для численной реализации такого решения необходимо избавиться от сингулярностей в интегральных соотношениях. Для устранения сингулярностей определяется порядок их полюсов, в соответствии с которым применяются соответствующие правила регуляризации. Таким образом, на основе метода погружения удается решить обратную задачу теплопроводности.

Полученные соотношения позволили получить аналитические соотношения и разработать численные алгоритмы:

1.Для интенсивности пучка через температуру поверхности нагретой мишени;

2. Для вектора координат центра тяжести интенсивности;

3. Для эффективного радиуса пучка;

4. Для функционала фокусировки, равного мощности пучка в пределах заданной апертуры.

С помощью данного метода проведена обработка результатов лабораторного эксперимента по восстановлению интенсивности пучка с помощью тепловизора.

В заключении диссертации приводятся основные результаты проведённых исследований.

НАУЧНАЯ НОВИЗНА

Впервые на основе вариационного алгоритма сконтруированы симметричные ортогональные вэйвлеты, подчиняющиеся быстрым алгоритмам преобразования и позволяющие сократить время работы адаптивного контура.

Впервые получен алгоритм позволяющий синтезировать симметричные биортогональные и комплексные вэйвлеты высокого порядка позволяющие визуализовать скрытые периодичности оптических сигналов и выделять контуры спутниковых оптических изображений.

Впервые на основе метода инвариантного погружения разработан метод расчета интенсивности лазерного пучка по температурному полю на поверхности мишени для двумерного и многомерного случаев при произвольных размерах мишени и произвольных граничных условиях на ее обратной поверхности. С помощью этого метода определены интегральные моменты и энергетические характеристики пучка. Впервые применен способ канонической регуляризации в процедуре построения многомерных обобщенных функций, возникающих при решении обратной задачи теплопроводности.

На основе представления дифференциальных операторов в вэйвлет базисах предложен метод восстановления коэффициентов разложения фазы, позволяющий эффективно оценивать относительный вклад различных типов аберраций контролируемых оптических систем без этапа восстановления фазы.

ПРАКТИЧЕСКАЯ ЗНАЧИМОСТЬ

Методы и алгоритмы, приведенные в работе, применимы для определения пространственно-временного распределения интенсивности лазерного излучения в турбулентной атмосфере, его интегральных моментов и энергетических функционалов на основе данных дистанционного измерения температуры поверхности, нагретой лазерным излучением. Созданные на основе этих методов пакеты прикладных программ позволили провести обработку результатов реального лабораторного эксперимента по тепловизионному измерению температуры нагреваемой поверхности и провести восстановление падающего теплового потока.

Разработанные алгоритмы и пакеты программ позволяют синтезировать ортогональные, биортогональные и комплексные вэйвлеты различных порядков и различных фрактальных размерностей, которые позволяют визуализовать масштабы неоднородностей и определять фрактальную размерность среды распространения оптического сигнала для последующего его моделирования.

Разработаны алгоритмы, позволяющие синтезировать адаптивные вэйвлет - пакеты, подстраивающиеся под особенности анализируемого сигнала, что весьма эффективно при оптимальной фильтрации сигнала.

Разработанные алгоритмы и пакеты программ позволяют из измерений наклонов волнового фронта непосредственно рассчитывать вклад отдельных модовых составляющих фазы в суммарные аберрации волновых полей, что отвечает многим практическим потребностям.

Найденные соотношения перехода от разложений в базисе Карунена-Лоэва-Обухова к разложению по полиномам Цернике, функциям Хаара, функциям Адамара (Уолша) и разработанные алгоритмы и программы позволяют осуществить статистически оптимальную модовую компенсацию турбулентных искажений фазы, оптимизировать и сократить время работы адаптивных систем, функционирующих в перечисленных базисах.

ДОСТОВЕРНОСТЬ результатов диссертации подтверждается их совпадением в частных случаях с результатами расчетов, выполненными другими авторами с помощью других подходов, удовлетворительным согласием результатов расчетов по разработанным алгоритмам и программам с данными лабораторного эксперимента.

Автором выносятся на защиту следующие

ОСНОВНЫЕ ЗАЩИЩАЕМЫЕ ПОЛОЖЕНИЯ

1. Статистически оптимальный базис КЛО минимизирующий дисперсию коэффициентов разложения оптических сигналов и учитывающий внешний масштаб турбулентности среды распространения, позволяющий оптимизировать работу корректирующих устройств и сократить время работы адаптивного оптического контура.

2.Вариационный алгоритм синтеза симметричных квазиортогональных вэйвлетов с желаемой гладкостью и обладающих свойством быстрых преобразований, позволяющих существенно сократить время синтеза и анализа нестационарных оптических сигналов.

3. Алгоритм синтеза биортогональных, симметричных вэйвлетов, обладающих высокой локальностью, как в пространственной, так и в частотной областях, с высоким числом нулевых моментов, и обеспечивающих визуализацию скрытых периодичностей и сингулярностей оптических сигналов, прошедших турбулентную атмосферу.

4.Алгоритм синтеза комплексных, ортогональных, симметричных вэйвлетов, обладающих высоким числом нулевых моментов и позволяющих осуществлять анализ и синтез комплексных оптических сигналов и выделять контуры спутниковых оптических изображений.

5.Дифференциальные и обратные им операторы, в том числе и обобщенные, в представлении вэйвлет-базисов, позволяющие восстанавливать волновой фронт лазерного пучка по его локальным наклонам и интенсивность пучка по температуре нагретой им мишени с учетом сингулярностей волнового фронта, при этом обеспечивается повышенная точность и устойчивость представления и сокращается время численной реализации.

6.Алгоритм синтеза адаптивных вэйвлет-пакетов обладающих высокой локальностью и оптимальных в смысле минимума энтропии.

7.Алгоритм решения методом погружения некорректной пространственно-временной обратной задачи теплопроводности.

Основные материалы диссертации опубликованы в следующих работах :

1. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Восстановление распределения интенсивности лазерного излучения по температуре поверхности секционированной мишени // Оптика атмосферы, 1991, т. 4, №2, с. 166172.

2. Аксенов В.П., Исаев Ю.Н. Преобразование Радона в задаче фазового оптического контроля // Оптика атмосферы. 1991, т. 4, № 12, с. 166-172.

3. Аксенов В.П., Исаев Ю.Н. Восстановление параметров лазерного пучка по температурному полю нагретой поверхности // Оптика атмосферы и океана, 1992, т. 5, № 5, с. 509-516.

4. Аксенов В.П., Исаев Ю.Н. Оптимальное модовое разложение фазы, восстановленной по измерениям наклонов волнового фронта в турбулентной атмосфере. I. Представление аберраций в базисе Карунена-Лоэва-Обухова // Оптика атмосферы и океана, 1994, т. 7, № 7, с. 947-954.

5. Аксенов В.П., Банах В.А., Захарова Е.В., Исаев Ю.Н. Оптимальное модовое разложение фазы, восстановленной по измерениям наклонов волнового фронта в турбулентной атмосфере. И. Погрешность алгоритмов и численный эксперимент // Оптика атмосферы и океана, 1994, т. 7, № 7, с. 955-959.

6. Исаев Ю.Н. Инженерный метод расчета распределения интенсивности излучения по температуре поверхности мишени // Оптика атмосферы и океана, 1994, т. 7, № 10, с. 1433-1436.

7. Исаев Ю.Н., Захарова Е.В. Восстановление интенсивности лазерного пучка по температурному полю тонкой мишени // Оптика атмосферы и океана, 1995, т. 8, № 6, с. 841-846.

8. Аксенов В.П., Банах В.А., Исаев Ю.Н., Захарова. Е.В., Тихомирова О.В. Дифракционный томографический датчик волнового фронта // Оптика атмосферы и океана, 1995, т. 8, № 12, с.1884-1888.

9. Исаев Ю.Н. Аналитический метод определения векторных ортогональных полиномов для градиента фазы при произвольной геометрии приемной апертуры // Оптика атмосферы и океана, 1995, т. 8, № 10, с. 1539-1541.

10. Исаев Ю.Н., Захарова Е.В. Восстановление интенсивности лазерного пучка по зашумленному температурному полю мишени. 4.1. Построение регуляризирующего алгоритма на основе БПФ // Оптика атмосферы и океана 1996, т. 9, № 10, с. 1353-1358.

11. Исаев Ю.Н., Захарова Е.В., Аксенов В.П., Исаков А.В., Рейно В.В, Цвык Р.Ш. Восстановление интенсивности лазерного пучка по зашумленному температурному полю мишени. 4.II. Обработка модельного и лабораторного экспериментов // Оптика атмосферы и океана 1996, т. 9, № 10, с. 1359-1366.

12. Исаев Ю.Н., Захарова Е.В. Представление функций Карунена-Лозва-Обухова в базисах Уолша и Хаара // Оптика атмосферы и океана, 1997, т. 10, №8, с. 959-966.

13. Исаев Ю.Н., Захарова Е.В. Синтез оптимального базиса для восстановления случайных волновых полей // Оптика атмосферы и океана, 1998, т. 11, № 5, с. 451-454.

14. Исаев Ю.Н., Захарова Е.В. Критерии эффективности адаптивных оптических систем при различных базисах разложения фазы случайной волны // Оптика атмосферы и океана, 1999, т. 12, № 8, с. 708-711.

15. Исаев Ю.Н., Синтез вэйвлет-базиса для анализа оптических сигналов Ч. I. Ортогональный вэйвлет-базис // Оптика атмосферы и океана, 2002, т. 15, № 11, с. 974-981.

16. Исаев Ю.Н. Синтез вэйвлет-базиса для анализа оптических сигналов Ч. II Биортогональный и комплексный вэйвлет-базисы // Оптика атмосферы и океана, 2003, т. 16, № 4, с. 329-336.

17. Исаев Ю.Н. Синтез вэйвлет-базиса для анализа оптических сигналов. Ч. III. Представление дифференциальных и обратных им операторов в вэйвлет-базисах. Вэйвлет-пакеты // Оптика атмосферы и океана 2003, т. 16, №8, с. 688-694.

18. Aksenov V.P., Isaev Yu.N. Analitical representation of the phase and its mode components reconstructed according to the wave front slopes // Optics Letters, 1992, v. 17, No 17, p. 1180-1182.

20. Аксенов В.П., Захарова E.B., Исаев Ю.Н. Измерение теплового потока по температурному полю нагретой поверхности. Ч. I. Однородный поток // Инженерно-физический журнал, 1994, т. 64, № 3-4, с. 275-280.

21. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Измерение теплового потока по температурному полю нагретой поверхности. Ч. II. Неоднородный поток // Инженерно-физический журнал, 1995, т. 68, № 4, с. 622-628.

22. Исаев Ю.Н., Колчанова В.А., Хохлова Т.Е. Определение параметров двухполюсника как эквивалентной схемы замещения электрического разряда при воздействии импульсного напряжения // Электричество, 2003, № 11, с. 63-69.

23. Тартаковский В.А, Исаев Ю.Н., Несвитайло В.Д., Волков Ю.В., Попов В.Н. Монотонность, причинность и демодуляция годичных колец деревьев // Автометрия, 2003, № 5, с. 118-127.

24. Яворовский Н.А., Исаев Ю.Н., Корнев Я.И., Чен-Бен Нам, Хасельберг М.Б. Определение параметров двухполюсника как эквивалентной схемы замещения электрического разряда при воздействии импульсного напряжения // Известия вузов «Физика», 2003, № 10, с. 3-7.

25. Исаев Ю.Н. Конструирование биортогональных вэйвлет базисов для оптимального представления сигналов // Известия Томского политехнического университета, 2004, т. 307, № 1, с. 37-42.

26. Исаев Ю.Н. Конструирование вэйвлет базисов для оптимального представления оптических изображений // Известия Томского политехнического университета, 2004, т. 307, № 2, с. 31-40.

27. Aksenov V.P., Isaev Yu.N. Reconstruction of laser radiation intensity distribution from temperature along target surface // SPIE's Intern. Conf. on Industrial Photonics. Thermosence XV. Orlando, USA, 1993, p. 298-307.

28. Aksenov V.P., Isaev Yu.N. Retreiving the laser beam intensity distribution from temperature field along of the heated target surface // CLEO Europe/EQEC'94 Abstract, Amsterdam, Holland, 1994, p.109-110.

29. Aksenov V.P., Isaev Yu.N., Zakharova E.V. Spatial-temporal reconstruction of laser beam intensity distribution from the temperature along surface of the heated target // SPIE's Intern. Conf. on Thermal Sensing and Imaging Diagnostic Applications XVIII, Orlando, USA, 1996, v. 2766, p. 336-358.

30. Isaev Yu.N., Zakharova E.V. Expansion of random wave phase in terms of eigenfunctions of phase correlation function // VIII International Symposium "Remote Sensing", Toulouse, France, 2001, p. 1011-1012.

31. Isaev Yu.N., Banakh V.A., Zakharova E.V. Orthogonal expansions of phase of random wave fields // SPIE Proceedings, Barcelona, Spain, 1998, V. 3494, p. 1234-1237.

32. Isaev Yu.N., Banakh V.A., Zakharova E.V. Representation of Karhunen-Loeve-Obukhov functions in bases of Walsh and Haar // SPIE Proceedings San Diego, California, USA, 1997, v. 3126, p. 507-514.

33. Isaev Yu.N., Zakharova E.V, Representation of phase of distorted optical wave through the orthonormal bases including the outer scale of turbulence //

Numerical experiment. VII International Symposium "Remote Sensing", Barcelona, Spain, 2000, p.107-108.

34. Аксенов В.П., Банах B.A, А.А.Землянов, Исаев Ю.Н. Определение моментов интенсивности и функционала фокусировки лазерного пучка по рассеянному излучению // IV Всесоюзный симпозиум по вычислительной томографии. Тез. докл. - Ташкент, 1989, с. 31.

35. Аксенов В.П., Исаев Ю.Н. Определение энергетического центра эффективного размера лазерного пучка по трем проекциям, зарегистрированным в пределах острого угла // IV Всесоюзный симпозиум по вычислительной томографии. Тез. докл. - Ташкент, 1989, с. 41.

36. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Восстановление тепловых потоков по температуре нагретой поверхности // Оптические методы измерений и способы обработки данных теплофизических и нейтронно-физических процессов в элементах энерготехники. Тез. докл. -Севастополь, 1990, с. 62.

37. Аксенов В.П., Исаев Ю.Н. Аналитический метод восстановления фазового распределения в лазерном пучке по измерениям локальных наклонов // Фотометрия и ее метрологическое обеспечение. 8 Всес. научно-техн. конференция. Тез. докл. - Москва, 1990, с. 62.

38. Аксенов В.П., Исаев Ю.Н. Преобразование Радона в задачах фазового контроля оптических полей // XI Всес. симпоз. по распространению лазерного излучения в атмосфере и водных средах. Тез. докл. - Томск, 1991, с.123.

39. Аксенов В.П., Исаев Ю.Н. Фазовый и модовый контроль в когерентной оптике на основе преобразования Радона // XI Всес. симпоз. по распространению лазерного излучения в атмосфере и водных средах. Тез. докл. - Томск, 1991, с.132-133.

40. Аксенов В.П., Исаев Ю.Н. Аналитическое решение многомерной задачи восстановления теплового потока по температуре нагретой поверхности Условно-корректные задачи математической физики и анализа. Всес. конф. Тез. докл. - Новосибирск, 1992, с. 44-45.

41. Аксенов В.П., Исаев Ю.Н. Томографические методы измерения фазы световых пучков в системах лазерной, астрономической оптики, оптического контроля // Симпоз. «Прикладная оптика-94». Тез. докл. -СПб, 1994, с. 231.

42. Исаев Ю.Н. Синтез вэйвлет базиса для анализа сигналов прошедших турбулентную атмосферу // Тез.док. Международная конференция ТГУ, "Сопряжённые задачи механики, информатики и экологии". Томск, 2003, с. 23.

43. Исаев Ю.Н. Представления дифференциальных и обратных им операторов в вэйвлет базисах для обработки оптических сигналов // X Юбилейный международный симпозиум «Оптика Атмосферы и Океана. Физика Атмосферы». - Томск, 2003, с. 70.

44. Исаев Ю.Н. Фирсов К.М. Вэйвлет анализ спектров поглощения атмосферных газов // X Юбилейный международный симпозиум «Оптика Атмосферы и Океана. Физика Атмосферы». - Томск, 2003, с. 49.

45. Исаев Ю.Н. Артамонов Е.С., Протасов К.Т. Адаптивное восстановление замутненных аэрокосмических изображений с использованием вейвлет преобразований // X Юбилейный международный симпозиум «Оптика Атмосферы и Океана. Физика Атмосферы». - Томск, 2003, с. 92.

46. Исаев Ю.Н. Захарова Е.В. Вэйвлет-анализ атмосферных метео-пара-метров // Х-Юбилейный международный симпозиум «Оптика Атмосферы и Океана. Физика Атмосферы». - Томск, 2003, с. 73.

47. Исаев Ю.Н., Колчанова В.А. Вэйвлет базис для анализа и синтеза сигналов // X Международная научно-практическая конференция студентов, аспирантов и молодых ученных, Т.2. - Томск, 2004, с. 149150.

АПРОБАЦИЯ РЕЗУЛЬТАТОВ

Основные результаты, полученные в диссертации, представлены в опубликованных работах [1-^$]. Докладывались на X, XI Всесоюзном симпозиуме по распространению лазерного излучения (Томск, 1989, 1991), IV Всесоюзном симпозиуме по вычислительной томографии (Ташкент, 1989), VIII Всесоюзной научно-технической конференции "Фотометрия и ее метрологическое обеспечение" (Москва, 1990), Всесоюзной конференции "Оптические методы измерений и способы обработки данных теплофизических и нейтронно-физических процессов в элементах энерготехники" (Севастополь, 1990), Всесоюзной конференции "Условно-корректные задачи математической физики и анализа" (Новосибирск, 1992), SPIE's International Conference of Industrial Photonics Thermosence XV (Orlando, USA, 1993), I Межреспубликанском симпозиуме "Оптика атмосферы и океана" (Томск, 1994), CLEO Europe/EQEC'94 (Amsterdam, Holland, 1994), Симпозиуме "Прикладная оптика-94" (С.-Петербург, 1994). II Межреспубликанском симпозиуме "Оптика атмосферы и океана" (Томск, 1995), III Межреспубликанском симпозиуме "Оптика атмосферы и океана" (Томск, 1996), SPIE's Intern. Conf. on Thermal Sensing and Imaging Diagnostic Applications XVIII, Orlando,(USA, 1996), CLEO Europe/EQEC'96 (Hamburg, Gergany,

1996), IV Симпозиуме "Оптика атмосферы и океана" (Томск,

1997), SPIE Proceedings San Diego, California (USA, 1997),VII Международном симпозиуме "Оптика атмосферы и океана" (Томск, 2000),VII International Symposium "Remote Sensing", Barcelona (Spain, 2000), VIII International Symposium "Remote Sensing", Toulouse, (France, 2001), Международной конференции ТГУ "Сопряжённые задачи механики, информатики и экологии"

Томск, 2002), IX-Международном научно-практической конференции студентов аспирантов и молодых ученных, ТПУ, (Томск, 2003), X Международном симпозиуме " Оптика атмосферы и океана. Физика Атмосферы ", (Томск, 2003). V-Сибирском совещании по климато - экологическому мониторингу, (Томск, 2003).

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Исаев, Юсуп Ниязбекович

ЗАКЛЮЧЕНИЕ

Приведём основные результаты проведённых исследований.

1. На основе аналитического соотношения, связывающего фазу с ее пространственным градиентом, получены формулы, позволяющие определять коэффициенты полиномиального разложения фазы непосредственно из измерения ее локальных наклонов (минуя этап восстановления фазы). Предложен полуаналитический подход к проблеме получения оптимальных в среднем по ансамблю модовых разложений случайной фазы в системах адаптивной оптики, на основе которого получены представления функций Карунена-Лоэва-Обухова.

2. Предложен способ оптимизации классического разложения Цернике, заключающийся в коррекции коэффициентов этого разложения на основе информации о пространственной корреляции в случайно неоднородной среде. Приведен явный вид матрицы преобразования полиномов Цернике в оптимальный базис. Показано, что такие оптимальные коэффициенты могут быть найдены непосредственно из измерений наклонов волнового фронта.

3. Исследована эффективность предлагаемого полуаналитического подхода к проблеме построения оптимального для турбулентной атмосферы базиса в представлении случайной фазы в системах когерентной адаптивной оптики. Описаны алгоритмы численной реализации, оценена погрешность, допускаемая при замене бесконечных рядов в представлении функции Карунена-Лоэва-Обухова рядом с конечным числом членов. Приведены результаты численного эксперимента и аналитический вид основных мод Карунена-Лоэва-Обухова в представлении через базис Цернике. Разработанные подходы и алгоритм могут служить основой для построения программных комплексов, обеспечивающих управление корректорами волнового фронта, функционирующими по принципам модовой или зональной коррекции.

4. Предложен алгоритм позволяющий вычислять матрицы перехода из оптимального базиса КЛО в быстрые базисы Уолша и Хаара, что позволяет существенно сократить время отклика адаптивной оптической системы.

5. На основе свойств локальности симметричных, ортогональных и биортогональных вэйвлетов позволяющих локально сужать функцию рассеяния точки изображения разработан новый алгоритм повышения контрастности изображения.

6. На основе свойств комплексных вэйвлетов позволяющих выделять сингулярности сигналов предложен новый алгоритм выделения контуров изображения.

7. На основе представления дифференциальных операторов в вэйвлет-базисах получен алгоритм восстановления волнового фронта с сингулярностями по его локальным наклонам.

8. Разработанные алгоритмы по моделированию и восстановлению волнового фронта в условиях сильной турбулентности могут служить основой для создания программного продукта, обеспечивающего функционирование адаптивной оптической системы в условиях сильной турбулентности, вызывающей появления вихревых сингулярностей в волновом фронте.

9. На основе метода инвариантного погружения получены формулы для восстановления интенсивности лазерного пучка по температуре нагретой поверхности бесконечно протяженной в поперечном направлении пластины для различных продольных длин и краевых условий. Аналитическое решение пространственно-временной обратной задачи теплопроводности в многомерной постановке получено впервые.

10. Решен вопрос регуляризации обобщенного решения многомерной обратной задачи пересчета граничных условий. Получены расчетные формулы для восстановления распределения интенсивности, интегральных моментов и энергетических функционалов лазерного пучка по температурному распределению на поверхности нагреваемой мишени.

11. Построены алгоритмы, учитывающие особенности решения обратной некорректной задачи, проведены численные и обработан лабораторный эксперименты по восстановлению тепловых потоков по зашумленным значениям температуры. Оценены точностные характеристики алгоритмов для характерных моделей временного хода интенсивности и показано, что использование сглаживающих сплайнов и алгоритмов, не содержащих дифференцирования значений температуры, повышает точность восстановления теплового потока. Разработанные алгоритмы могут служить основой для создания программного продукта, обеспечивающего функционирование соответствующих измерительных систем.

Список литературы диссертационного исследования доктор физико-математических наук Исаев, Юсуп Ниязбекович, 2004 год

1. Применение адаптивной оптики для коррекции изображения наблюдаемого объекта // Радиоэлектроника за рубежом, 1978. Вып.5. С.20-31.

2. Пузырев В.А., Данилевич А.Б. Системы автоматической коррекции волнового фронта // Зарубежная радиоэлектроника, 1980. № 6. С. 45-63.

3. Рознер Р.Д. Гетеродинный прием оптического сигнала, распространяющегося в атмосфере // ТИИЭР, 1968. Т.5. № 1. С. 144-146.

4. Татарский В.И. Адаптивные системы и когерентность I, II // Изв. ВУЗов. Радиоэлектроника, 1981. Т.24. № 7. С.861-883.

5. Фрид Д. Гетеродинный прием оптического сигнала при атмосферных искажениях волнового фронта // ТИИЭР, 1978. Т.66. № 6. С.31-85.

6. MullerR.A., BufflngtonA. Real time correcting of atmospherical telescope image through image sharpening. // J. Opt. Soc. Am. 1974. V.64. № 9. P. 1200-1210.

7. T. Stewart McKechnie. Light propagation through the atmosphere and the properties of images formed by of large ground-based telescopes // J. Opt. Soc. Am. 1991. V.8. № 2. P. 346-364.

8. T. Stewart McKechnie. Atmospheric turbulence and the resolution limits of large ground-based telescopes // J. Opt. Soc. Am. 1992. V.9. № 11. P. 19371954.

9. Tatarskii V.I. and Zavorotny V.U. Atmospheric turbulence and the resolution limits of large ground-based telescopes: comment // J. Opt. Soc. Am. 1993. V.10. № 11. P. 2410-2417.

10. J.H. Tarng and Chih-Ming Chen. Recovery of laser beam propagation through a turbulent amplifier using phase cojugation // Opt. Engineering. 1993. V.32. № l.P. 121-130.

11. Латынин Ю.М., КолтокЮ.В., Кузьмичев B.M. Пиромагнитные детекторы излучения // Квантовая электроника, 1980. Вып. 19. 88-97 с.

12. Колток Ю.В., Кузьмичев В.М., Латынин Ю.М., Манита О.Ф. Пиромагнитный детектор излучения // Квантовая электроника, 1978. Вып. 14. 82-84 с.

13. Аполлонов В.В., Бункин Ф.В., Хомич В.Ю., Четкин С.А. Термодеформационный метод измерения распределения интенсивности мощного лазерного пучка // Письма в ЖТФ, 1978. Т.4. Вып.17. 10171019 с.

14. Измерение энергетических параметров и характеристик лазерного излучения //Под ред. А.Ф.Котюка. М.: Радио и связь, 1981. 288с.

15. Григорьев А. Б. Импульсный нагрев излучениями. Ч. I. М.: Наука, 1974. 320 с.

16. Григорьев А. Б. Импульсный нагрев излучениями. Ч. II. М.: Наука, 1979. 727 с.

17. Абильшитов Г.Л., Велихов Е.П., Голубев B.C. и др. Мощные газоразрядные СС>2-лазеры и применение в технологии. М.: Наука, 1984.312с.

18. Аксенов В.П., Пикалов В.В. Томографическое восстановление пространственно-энергетических параметров лазерных пучков // Квантовая электроника, 1990. Т.14. № 2. С. 167-172.

19. Ллойд Дж. Системы тепловидения. М.: Мир, 1978. 414 с.

20. Кузиковский А. В., Погодаев В. А., Рождественский А. Е. // II Совещание по атмосферной оптике. Томск: ИОА СО АН СССР, 1980. Ч. И. С. 177.

21. Мурин Г.А. Теплотехнические измерения. М.: Энергия, 1979. 423с.

22. Госсорг Ж. Инфракрасная термография. М.: Мир, 1988. 400с.

23. Воронцов М.А., КорябинА.В., Шмальгаузен В.И. Управляемые оптические системы. М.: Наука, 1988. 272 с.

24. Тараненко В.Г., Шанин О.И. Адаптивная оптика. М.: Радио и связь, 1990. 112с.

25. Воронцов М.А., Шмальгаузен В.И. Принципы адаптивной оптики. М.: Наука, 1985. 336 с.

26. ФервердаХ.А. Обратные задачи в оптике. М.: Машиностроение, 1984. 213с.

27. Barakat R., Sandler В.Н. Determination of the wave-front aberration function from measured values of the point-spread function: a two-dimensional phase retreival problem// J. Opt. Soc. Am. 1992. V.9. № 10. P. 1715-1723.

28. Адаптивная оптика // Под ред. Д.Фрида. М.: Мир, 1980. 456с.

29. Харди Дж.У. Активная оптика: новая техника управления световым пучком // ТИИЭР. 1978. Т.66. № 6. С. 31-85.

30. Аксенов В.П. Аналитическое представление фазы, восстановленной по производным волнового фронта // В сб. "XVI Всесюзн. конф. по распространению радиоволн". Харьков. 1990. 4.2.

31. Исаев Ю.Н. Аналитический метод определения векторных ортогональных полиномов для градиента фазы при произвольной геометрии приемной апертуры //Оптика атмосферы и океана. 1995, т.8, N 10. с.1539-1541

32. Beylkin G., Coifman R., Daubechies I., Mallat S., Meyer Y., Raphael L., and Ruskai B. Wavelets and Their Applications, Jones and Bartlett, Cambridge, MA, 1992.

33. Daubechies I. Wavelets, CBMSNSF Series in Appl. Math. SIAM Publ., Philadelphia, 1992.

34. Meyer Y. Omdelettes et Operateurs, in two volumes, Hermann, Paris, 1990.

35. Chui С. K. (ed.), Approximation Theory and Functional Analysis, Academic Press, Boston, 1991.

36. M.Antonini, M.Barlaud, P.Mathieu, and I.Daubechies.Image coding using wavelet transforms, IEEE Trans. Image Process., 1. pp. 205-220. 1992.

37. A.Cohen and I.Daubechies. A stability criterion for biortogonal wavelet bases and their related subband coding schemes, Duke Math.J., 68. pp. 313-335.

38. Лукьянов Д. П., Корниенко А. А., Рудницкий Б.Е. Оптические адаптивные системы. М.: Радио и связь, 1989. 238с.

39. Голубов Б.И., Ефимов А.В., Скворцов В.А. Ряды и преобразования Уолша: Теория и применения. М.: Наука, 1987. 344 с.

40. Залманзон Л.А. Преобразование Фуье, Уолша, Хаара и их применения в управлениии, связи и других областях М.: Наука 1989. 496 с.

41. Новиков И.Я., Стечкин С.Б. Основы теории всплесков// Успехи Математических Наук. 1998. Т.53. № 6. С.53-128.

42. Петухов А.П. Введение в теорию базисов всплесков. СПб: Издательство СПбГТУ, 1999. 132с.

43. Добеши И. Десять лекций по вейвлетам. Пер. с англ. Е.В. Мищенко. Под ред. А.П. Петухова. М.:РХД, 2001

44. Чуй К. Введение в вэйвлеты Пер. с англ. под. ред. Я.М. Жилейкина. М.:Мир, 2001

45. R.A.DeVore, B.Jawerth and B.J.Lucier. "Image compression through wavelet transform coding", IEEE Transactions on Information Theory, 1992, Vol.38. N2 pp.719-746.

46. Исаев Ю.Н. Синтез вэйвлет-базиса для анализа оптических сигналов Часть I. Ортогональный вэйвлет базис.// Оптика атмосферы и океана 2002, Т.15, N11. с.974-981

47. Исаев Ю.Н. Синтез вэйвлет-базиса для анализа оптических сигналов Часть II. Биотогональный и комплексный вэйвлет базисы.// Оптика атмосферы и океана 2003, Т. 16, N4. с.329-336

48. Э. Столниц, Т.ДеРоуз, Д.Салезин. Вейвлеты в компьютерной графике. РХД Москва-Ижевск. 2003. 271с.

49. Исаев Ю.Н. Артамонов Е.С., Протасов К.Т. Адаптивное восстановление замутненных аэрокосмических изображений с использованием вейвлет преобразований , Х-Юбилейный международный симпозиум «Оптика Атмосферы и Океана. Физика Атмосферы», 2003, 24-28 Июня

50. Исаев Ю.Н. Синтез вэйвлет-базиса для анализа оптических сигналов. Часть III. Представление дифференциальных и обратных им операторов в вэйвлет базисах. Вэйвлет пакеты// Оптика атмосферы и океана, 2003, N8, 688-694с.

51. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Измерение теплового потока по температурному полю нагретой поверхности. I. Однородный поток // ИФЖ, 1994. Т. 64. № 3-4. С. 275-280.

52. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Измерение теплового потока по температурному полю нагретой поверхности. II. Неоднородный поток // ИФЖ, 1995. Т. 68. № .4 С.622-628

53. Аксенов В.П., Исаев Ю.Н. Аналитическое решение многомерной задачи восстановления теплового потока по температуре нагретой поверхности // В кн.: Условно-корректные задачи математической физики и анализа. Всес. конф. Тез. докл. Новосибирск. 1992.

54. Aksenov V.P., Isaev Yu.N. Retreiving the laser beam intensity distribution from temperature field along of the heated target surface // CLEO Europe/EQEC'94 Abstract. Amsterdam, Holland. 1994.

55. Алифанов O.M. Прямые методы решения линейных обратных задач теплопроводности // ИФЖ. 1973. Т. 24. N6. С. 11241128.

56. Aksenov V.P., Isaev Yu.N. Reconstruction of laser radiation intensity distribution from temperature along target surface // SPIE's Intern. Conf. on Industrial Photonics. Thermosence XV. Orlando, USA. 1993.

57. Аксенов В.П., Исаев Ю.Н. Восстановление параметров лазерного пучка по температурному полю нагретой поверхности // Оптика атмосферы и океана. 1992. Т. 5. № 5. С. 509-516.

58. Касти Дж., Калаба Р. Методы погружения в прикладной математике.М.: Мир, 1976. 233с.

59. H.Kagiwada, R.Kalaba and B.Vereeke. "The Invariant Imbedding Numericaal Method for Fredholm Integral equations with Generate Kernels" // Journal of Approximation Theory. V.l. №1. 1970. P.355-364.

60. Кляцкин В.И., Метод погружения в теории распространения волн. М.: Наука, 1986. 256с.

61. Витреченко Э.А., Лукин В.П., Пушной Л.А., Тартаковский В.А. Проблемы оптического контроля. Новосибирск: Наука, 1990. 351 с.

62. Аксенов В.П., Исаев Ю.Н. Аналитический метод восстановления фазового распределения в лазерном пучке по измерениям локальных наклонов // В кн.: Фотометрия и ее метрологическое обеспечение. 8 Всес. научно-техн. конференция. Тез. докл. Москва. 1990.

63. Аксенов В.П., Исаев Ю.Н. Преобразование Радона в задачах фазового контроля оптических полей // XI Всес. симпоз. по распространению лазерного излучения в атмосфере и водных средах. Тез. докл. Томск. 1991.

64. Аксенов В.П., Исаев Ю.Н. Фазовый и модовый контроль в когерентной оптике на основе преобразования Радона // XI Всес. симпоз. по распространению лазерного излучения в атмосфере и водных средах. Тез. докл. Томск. 1991.

65. Аксенов В.П., Исаев Ю.Н. Томографические методы измерения фазы световых пучков в системах лазерной, астрономической оптики, птического контроля // Симпоз. Прикладная оптика-94. Тез. докл. С.-Петербург. 1994.

66. Аксенов В.П., Банах В.А., Захарова Е.В., Исаев Ю.Н. Аналитическое представление аберраций волнового фронта в турбулентной атмосфере в базисе Карунена-Лоэва-Обухова. // В кн.: Оптика атмосферы и океана. I Межресп. симпоз. Тез. докл. Томск. 1994.

67. Aksenov V.P., Isaev Yu.N. Analitical representation of the phase and its mode components reconstructed according to the wave front slopes // Optics Letters. 1992. V. 17. № 17. P.

68. Аксенов В.П., Исаев Ю.Н. Преобразование Радона в задаче фазового оптического контроля // Оптика атмосферы. 1991. Т. 4. № 12. С. 166-172.

69. Пикалов В.В., Преображенский Н.Г. Реконструктивная томография в газовой динамике и физике плазмы. Новосибирск: Наука, 1987. 226 с.

70. Тихонов А.Н., Арсенин В.Я., Тимонов А.А. Математические задачи компьютерной томографии. Москва: Наука, 1987. 160 с.

71. Воронцов М.А., Корябин А.В., Шмальгаузен В.И. Управляемые оптические системы. М. : Наука, 1988. 272 с.

72. Борн М., Вольф Э. Основы оптики. М.: Наука, 1973. 719 с.

73. Wang J. Y., Markey J. K.//J. Opt. Soc. Amer. 1978. V. 68. N. 1. P. 78

74. Ван Трис Г. Теория обнаружения оценок и модуляции. М.: Советское радио, 1972. 744с.

75. Обухов А. М. Турбулентность и динамика атмосферы. Ленинград: Гидрометеоиздат, 1988. 414с.

76. Воронцов М. А., Шмальгаузен В. И. Принципы адаптивной оптики. М.: Наука, 1985. 334с.

77. Татарский В. И. Распространение волн в турбулентной атмосфере. М.: Наука, 1967. 548с.

78. ЦлафЛ.Я. Вариационные исчисления и интегральные уравнения. Москва: Наука, 1966. 176с.

79. Верлань А.Ф., Сизиков B.C. Интегральные уравнения. Методы, алгоритмы, программы. Киев: Наукова думка. 1986. 543 с.

80. Виленкин Н. Я. Специальные функции и теория представлений групп. М.: Наука, 1965. 588с.

81. Никифоров А.Ф., Уваров В.Б. Специальные функции математической физики. Москва: Наука, 1984. 343 с.

82. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. М.: Наука, 1987. 598с.

83. Исаев Ю.Н., Захарова Е.В. Представление функций Карунена-Лоэва-Обухова в базисах Уолша и Хаара // Оптика атмосферы и океана. 1997. Т. 10. № 8. С. 959-966.

84. Исаев Ю.Н., Захарова Е.В. Синтез оптимального базиса для восстановления случайных волновых полей// Оптика атмосферы и океана 1998 т 11, N 5. С.451-454.

85. Voitsekhovich V.V. Outer scale of turbulence: comparison of different models// JOSA A. 1995. Vol. 12. N 6. P.1346-1353.

86. Winker D.M. Effect of a finite outer scale on the Zernike decomposition of atmospheric optical turbulence// JOSA A. 1991. Vol.8. N 10. P. 1586-1517.

87. Брычков Ю.А., Маричев О.И., Пруднико А.П. Интегралы и . ряды. Специальные функции. Москва: Наука, 1983. 437с.

88. Безуглов Д.А., Мищенко Е.Н., Серпенинов О.В. // Оптикаатмосферы. 1991. Т.4. N 2. С.161.

89. Безуглов Д.А. // Оптика атмосферы. 1991. Т.4. N 12 С.1303.

90. Edward P.Wallner //JOSA. V. 73. N 12. P.l771-1776.

91. D.P.Peterson, K.H.Cho // JOSA, 1986. V. 3. N 6. P.818-825.

92. Бицадзе А.В. "Основы теории аналитических функций комплексного переменного". М.: Наука, 1984. 320с.

93. Аксенов В.П., Исаев Ю.Н. //Оптика атмосферы. 1991. Т.4. N 12 С.1321-1325.

94. Исаев Ю.Н. "Разработка методов и численно-аналитическое моделирование восстановления параметров лазерных пучков". Канд. дисс. Томск: ИОА. 1995. С.135.

95. Gavrielides // Optics Letters, 1982, V. 7. N 11. P.526-528.

96. Тараненко В.Г., Шанин О.И. Адаптивная оптика. М.: Радио и связь, 1990. 112с.

97. Астафьева Н.М. Вейвлет-анализ: спектральный анализ локальных возмущений (основы теории и примеры применения). // Известия вузов "Проблемы нелинейной динамики". 1996. Т. 4. №2. С.3-39.ч

98. Mallat. S.G. Theory for multiresolution signal decomposition: The wavelet representation. //IEEE Transactions on Pattern Analysis and Machine Intelligence. 1989 V.31(7). P. 674-693.

99. Collineau S., Brunet Y. Detection on turbulent coherent motions in a forest canopy. Part I. Wavelet analysis. //Boundary-Layer Meteorlogy.1993. V.65. P. 357-379.

100. Noll R.J. Zernike and atmospheric turbulence. // J. Opt. Soc. Amer. 1976. V. 66.N.3. P. 207-211.

101. Ярославский Л.П. Цифровая обработка сигналов в оптике и голографии: Введение в цифровую оптику. М.: Радио и связь, 1987. 296с.

102. А.М.Трахтман, В.А.Трахтман. Основы теории дискретных сигналов на конечных интервалах. М.: Советское радио, 1975. 208с.

103. Калиткин Н.Н. Численные методы. М.: Наука, 1978. 512с.

104. Шишаков К.В., Шмальгаузен В.И. Полиномиальное разложение атмосферных аберраций // Оптика атмосферы и океана. 1990. Т.З. №12. С. 1244-1248.

105. Dai Guang-ming. Modal compensation of atmospheric turbulence with the use of Zernike polynomials and Karhunen-Loeve functions // JOSA A. 1995. V. 12. N. 10. P. 2182-2193.

106. Wang J.Y. Phase-compensated optical beam propagation throught atmospheric turbulence // Appl. Opt. 1978. 17. N 16. P. 2580

107. Астафьева Н.М. Вейвлет-анализ: спектральный анализ локальных возмущений (основы теории и примеры применения). // Известия вузов "Проблемы нелинейной динамики". 1996. Т. 4. №2. С.3-39.

108. Ш.Новиков И .Я., Стечкин С.Б. Основы теории всплесков// Успехи Математических Наук. 1998. Т.53. № 6. С.53-128.

109. Дрёмин И.М., Иванов О.В., Нечитайло В.А. Вейвлеты и их использование// Успехи Физических Наук. 2001. Т.171. № 5. С.465-501.

110. Wickerhauser M.V. Two fast approximate wavelet algorithms for image processing, classification, and recognition.// Optical engineering. 1994. V. 33, N 7. P.2225-2235.

111. Френке JI. Теория сигналов. М: Советское радио, 1974. 343 с.

112. Ахмед Н., Рао К.Р. Ортогональные преобразования при обработке цифровых сигналов. М: Связь, 1980. 247 с.

113. Mingui Sun., Robert J. Sclabassi, Ching-Chung Li Yingping Zhang, Harold H.Szu. Symmetrical wavelet transforms for edge localization. .// Optical engineering. 1994. 1994. V. 33, N 7, P.2272-2281.

114. Gilbert Strang and Vasily Strela. Ortogonal multiwavelets with vanishing moments.// Optical engineering. 1994. 1994. V. 33, N 7, P.2104-2107.

115. DeVore R.A., Jawerth B. and LucierB.J. "Image compression through wavelet transform coding", IEEE Transactions on Information Theory, 1992, Vol.38. N2 pp.719-746.

116. Хургин Я.И Яковлев В.П. Финитные функции в физике и технике. М.:Наука, 1971,с.406.

117. Минкович Б.М., Хургин Я.И. Теория синтеза антенн, М: Советское радио, 1969, с.236.

118. Кеч В., Теодореску П. Введение в теорию обобщённых функций с приложениями в технике. Москва: Мир, 1978. 518 с.

119. Гельфанд И.М., Шилов Г.Е. Обобщенные функции и действия над ними. Москва: Физматгиз, 1959, 470 с.

120. Владимиров B.C. Обобщённые функции в математической физике. М: Наука, 1976.370 с.

121. Jian Lu, Dennis М. Healy, John В. Weaver. Contrast enhancement of medical imadiges using multiscale edge representation. // Optical Engineering. 1994, V.33, N7, p.2151-2161.

122. Канев Ф.Ю., Лукин В.П., Макенова Н.А. Регистрация фазового профиля когерентного излучения и реализация адаптивного управления лазерным пучком при наличии особых точек в волновом фронте./Юптика атмосферы и океана 2002, T.15,N11. с.1018-1026.

123. Тартаковский В.А., Сенников В.А. Условие появления и плотность действительных нулей флуктуирующей волны.//2002. Т.15, N 11. С. 969993.

124. Gerald Kaiser. A friendly guide to wavelets. Cambridge, MA 1994.

125. Gerald Kaiser. An algebraic theory of wavelets, Part I: Complex structure and operational calculus, SIAM J.Math. Anal. 23 (1992), 222-243, 1992a

126. Тихонов A.H., Гончарский A.C., Степанов B.B., ЯголаА.Г. Регуляризирующий алгоритмы и априорная информация. М.: Наука, 1983. 200 с.

127. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука. 1986. 286с.

128. Попов Д.А. О применении гладких регуляризаторов для вычисления свертки. ДАН СССР, 1984. Т.276, N 1, с.38-42.

129. Coifman R.R., Meyer Y., Wickerhauser V.Wavelet Analysis and Signal Processing. Wavelets and their Applications, Jones and Barlett, Boston. 1992. P. 153-178.

130. Chui Ch.K., Li Chun. Nonortogonal wavelet wavelet packets, SIMAJ. Math. Anal V. 24 (1993). № 3. P. 712-738.

131. Теоретические основы и конструирование численных алгоритмов задач математической физики. Под редакцией К.И. Бабенко. М.: Наука. 1997.295.C.

132. Диткин В.А., Прудников А.П. Интегральные преобразования и операционное исчисление. М.: Наука, 1974. 544 с.

133. Сороко Л.М. Основы голографии и когерентной оптики. М.: Наука. 1971. 616 с.

134. Папулис А. Теория систем и преобразований в оптике. М.: Мир. 1971, 496 с.

135. Мамфорд А.В. Лекции о тета-функциях. М.: Мир, 1988, 448 с.

136. Привалов И.И. Введение в теорию функцию комплексного переменного. М.: Наука. 1884,432 с.

137. Прудников А.Д., БрычковЮ.А., Маричев О.И. Интегралы и ряды. Дополнительные главы. М.: Наука, 1986. 800 с.

138. Лаврентьев М.А., Шабат Б.В. Методы теории функции комплекного переменного М.: Наука, 1987, 688 с.

139. Алифанов О.М. Идентификация процессов теплообмена летательных аппаратов. М.: Машиностроение, 1979. 216 с.

140. Шилов Г.Е. Математический анализ. Второй специальный курс. М.: Наука, 1965. 327 с.

141. Нуссенцвейг Х.М. Причинность и дисперсионные соотношения. М.: Мир, 1976 с.

142. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука. 1989. 623 с.

143. Бохнер С. Лекции об интегралах Фурье. М.: Мир, 1962, 360 с.

144. Ю.Н.Исаев. "Разработка методов и численно-аналитическое моделирование восстановления параметров лазерных пучков". Канд. дисс. Томск: ИОА. 1995. С. 135.

145. Г.И.Василенко, А.М.Тараторин. Восстановление изображений. М.: Радио и связь, 1986. 302с.

146. Ж.Макс. Методы и техника обработки сигналов при физических измерениях. М.: Мир, 1983. Т. 1. 311с.

147. В.А.Морозов. Регулярные методы решения некорректно поставленных задач. М.: Наука, 1987. 240с.

148. В.А.Морозов. // ЖВФ и МФ. 1966. Т.6. №1.

149. В.А.Морозов. // ЖВФ и МФ. 1974. Т.14. №2.

150. В.К.Иванов, В.В.Васин, В.П.Танана. Теория линейных некорректных задач и ее приложения. М.: Наука, 1978.

151. В.П.Танана. Методы решения операторных уравнений. М.: Наука, 1981.

152. Доботкин А.А., Исаков А.В., Ильин А.Б., Петренко А.П., Рейно В.В., Цвык Р.Ш., Шерстобитов М.В. // Оптика атмосферы и океана, 1994. Т.7. №5. С.673-679.

153. Исаков А.В., Ильин А.Б., Петренко А.П., Рейно В.В., Цвык Р.Ш., Шерстобитов М.В. // В кн.: XII Межресп. симпоз. по распространению лазерного излучения в атмосфере и водных средах. Тез. докл. Томск. 1993. С.188.

154. Аксенов В.П., Захарова Е.В., Исаев Ю.Н. Восстановление распределения интенсивности лазерного излучения по температуре поверхности секционированной мишени // Оптика атмосферы. 1991. Т. 4. № 2. С. 166172.