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

доктора технических наук
Закасовская, Елена Владимировна
город
Владивосток
год
2011
специальность ВАК РФ
05.13.18
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Нейросетевое и алгебраическое моделирование в распределенных волоконно-оптических измерительных системах в условиях неполноты данных»

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

на правах рукописи

4851343

ЗАКАСОВСКАЯ Елена Владимировна

НЕЙРОСЕТЕВОЕ И АЛГЕБРАИЧЕСКОЕ МОДЕЛИРОВАНИЕ В РАСПРЕДЕЛЕННЫХ ВОЛОКОННО-ОПТИЧЕСКИХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМАХ В УСЛОВИЯХ НЕПОЛНОТЫ ДАННЫХ

Специальность: 05.13.18 - Математическое моделирование, численные методы и комплексы программ

Автореферат

диссертации на соискание ученой степени доктора технических наук

3 О [ДОН 2011

ВЛАДИВОСТОК-2011

4851343

Работа выполнена в лаборатории прецизионных оптических методов измерений Института автоматики и процессов управления ДВО РАН

Научный консультант: член-корреспондент РАН,

доктор физико-математических наук, профессор Кульчин Юрий Николаевич

Официальные оппоненты:

Ведущая организация: Научно-исследовательский институт системных исследований Российской академии наук (НИИСИ РАН), г. Москва

Защита состоится " 4 " июля 20 И г. в 12 час. на заседании диссертационного совета Д 005.007.01 в Институте автоматики и процессов управления ДВО РАН по адресу: 690041, г.Владивосток, ул. Радио, 5.

С диссертацией можно ознакомиться в библиотеке Института автоматики и процессов управления ДВО РАН.

Автореферат разослан »2011г.

Ученый секретарь диссертационного совета

член-корреспондент РАН,

доктор физико-математических наук, профессор

Федотов Анатолий Михайлович

доктор технических наук, профессор Еремин Евгений Леонидович

доктор технических наук, профессор Девятисильный Александр Сергеевич

Д 005.007.01, к.т.н.

А.В. Лебедев

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Широкие возможности в подобных задачах открывает использование распределенных волоконно-оптических измерительных систем (РВОИС) на основе волоконно-оптических измерительных сетей (ВОИС) [1]. Такие системы представляют собой эффективный инструмент для исследования и контроля внешних физических воздействий на протяженные объекты в режиме реального времени.

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

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

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

в выборе схемы укладки ИЛ - обычно проекции известны только в нескольких направлениях при малом числе ИЛ в каждом из них.

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

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

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

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

Следует также заметить, что НС не подходят для решения большинства классических томографических задач по причине значительного числа направ-

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

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

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

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

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

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

Достижение указанной цели требует решения следующих задач:

1. Разработка алгебраических моделей и методов восстановления параметров ФП для нерегулярных распределенных ВОИС томографического типа, основанных на обобщенных теоремах дискретизации и аппроксимации функции проекции на нерегулярной сетке.

2. Разработка метода алгебраического синтеза дополнительных проекционных данных с помощью аппроксимации функции проекции на объединении смежных классов с целью устранения глобальных артефактов для малоракурсных схем укладки ИЛ;

3. Разработка моделей и методов обработки сигналов нерегулярных распределенных ВОИС томографического типа с помощью полнообразных НС радиально-базисного типа, обеспечивающих реконструкцию пространственных распределений параметров ФП в реальном времени;

4. Разработка моделей и методов нейросетевой предобработки сигналов распределенных ВОИС томографического типа в пространстве Радона, обеспечивающих реконструкцию функции проекции в условиях неполноты данных и дальнейшего применения методов аппроксимации для нерегулярных схем, а также создание на их основе нового комбинированного нейро-алгебраического алгоритма восстановления параметров ФП;

5. Создание нового алгоритма восстановления параметров ФП, сочетающего в себе алгебраическую оптимизацию параметров ВОИС с целью локализации мест внешнего воздействия на ВОИС и аппроксимацию функции проекции на нерегулярной сетке.

6. Разработка и исследование новых специализированных нейросете-вых конструкций - комплексов нейронных сетей (КНС):

- для восстановления параметров ФП,

- для синтеза синограмм,

- для классификации измерительной информации;

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

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

Основные положения, выдвигаемые на защиту:

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

2. Разработан и исследован новый метод устранения глобальных артефактов для малоракурсных схем укладки ИЛ в ВОИС, основанный на алгебраическом синтезе синограмм с использованием аппроксимации функции проекции на нерегулярной сетке.

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

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

5. Впервые введена и исследована модель комплекса нейронных сетей для обработки данных РВОИС. Предложены различные по своим функциям модели комплексов нейронных сетей: для стандартизации измерительных данных, синтеза синограмм и восстановления параметров физических полей.

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

Научная новизна.

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

2. Разработан метод устранения глобальных артефактов с помощью алгебраического синтеза синограмм для ультрамалоракурсных схем укладки ИЛ в ВОИС.

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

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

ВОИС томографического типа, обеспечивающих реконструкцию функции проекции на примере НС радиально-базисного типа.

5. Предложен новый комбинированный алгоритм обработки проекционных данных для реконструкции функций распределения физических полей. Этот алгоритм заключается в последовательном выполнении двух процессов: предобработки измерительной информации с помощью НС ра-диально-базисного типа и дальнейшей аппроксимации полученных данных на множестве дискретизации специального вида - на объединении классов смежности.

6. Впервые введено и исследовано понятие комплекса нейронных сетей. Исследованы различные по своим функциям комплексы нейронных сетей для стандартизации измерительных данных, синтеза синограмм, восстановления параметров физических полей.

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

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

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

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

Изучение принципов применения нейросетевых вычислительных методов для обработки сигналов распределенных ВОИС показало возможность реализации параллельных алгоритмов обработки информации при реконструкции пространственных распределений ФП.

Полученные автором результаты были применены для создания моделей нейроподобных вычислительных систем обработки сигналов распределенных ВОИС и могут быть в дальнейшем использованы для создания систем мониторинга протяженных ФП в режиме реального времени.

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

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

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

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

Апробация работы. Изложенные в диссертационной работе материалы докладывались на международных конференциях: "International workshop on Optical Beam Transformation (IWBT)" Vladivostok, Russia, 2001; "International workshop on Bio-Signal Analysis and its Applications", Tokyo, Japan, 2002; "Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics", Vladivostok, Russia, 2002; World Exhibition of Invention, Research and Industrial Innovation "Eureka-2002" in Brussels, 12-17 November, 2002; "IEEE, Physics and Control", 2003; "Optoelectronic information-power technologies", Vinnitsa, Ukraine, 2003; "Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics", Russia, Khabarovsk, 2004;Second International Conference "Physics and control", 2005; "7-th International Conference on Digital Signal Processing and its Applications ", Moscow , Russia, 2005; "6-th Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics" 2007; "7-th Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics" 2008; "8-th Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics", 2009; "9-th Asia-Pacific Conference on Fundamental Problems of Opto- and Microelectronics", Kangnam University, Yongin, Republic of Korea, 2010; "The First Russia and Pacific Conference on Computer Technology and Applications", Vladivostok, Russia, 2010 и

на всероссийских конференциях:

ХЬУ1 Всероссийская межвузовская научно-техническая конференция, Владивосток, ТОВМИ, 2003 ;"Х Всероссийская научно-техническая конференция Нейроинформатика-2009"; "XI Всероссийская научно-техническая конференция Нейроинформатика-2010"; "XII Всероссийская научно-техническая конференция Нейроинформатика-2011".

Публикации. По теме диссертации опубликовано 33 статьи, из которых 13 входят в Перечень ВАК РФ.

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

Структура и объём диссертации. Диссертация состоит из пяти глав, заключения и приложения. Объем диссертации составляет 245 страниц, включает 53 рисунка, 11 таблиц и список литературы из 154 названий.

СОДЕРЖАНИЕ РАБОТЫ

Глава 1 носит вводный характер. В ней рассматривается использование современных измерительных технологий на основе волоконной оптики для создания информационно-измерительных систем, способных осуществлять оперативный сбор и обработку информации о многомерных пространственно-временных функциях распределения параметров различных ФП. Обсуждаются возможности использования РВОИС на основе ВОИС для исследования и контроля внешних физических воздействий на протяженные объекты в режиме реального времени.

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

витальную поверхность сети и отслеживать его перемещения по контролируемой поверхности. В §1.1.3 приведена краткая классификация основных методов классической томографии, применяемых для восстановления параметров ФП по информации, поступающей с РВОИС. В §1.1.4 дается обоснование применения нейросетевых технологий в РВОИС.

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

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

Глава 2 посвящена современным алгебраическим методам восстановления информации, поступающей по ВОИЛ распределенной измерительной сети. При решении большинства таких задач схемы укладки ИЛ в ВОИС далеки от идеальных равномерных схем. Часто встречаются ограничения, в виде уменьшенного шага сканирования, более редкого диапазона углов (рис. 1 а-г), неравномерных отсчетов в разных направлениях укладки ИЛ (рис. \в-г), которые и определяют вид массивов измерительной ин-

Рис. 1. Параллельные схемы укладки измерительных линий в ВОИС

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

нове [4, 5]. В главе 2 рассматривается задача восстановления функций ФП с помощью ВОИС [1] Для неравномерных параллельных схем укладки ИЛ с малым числом направлений сканирования.

Одним из способов сбора информации о состоянии ВОИС является получение синограммы путем сканирования области набором параллельных ИЛ (рис.1с>). Рассмотрим воздействие ФП на распределенную измерительную сеть, периферийная часть которой составлена из интегрирующих ИЛ. Предположим, что ФП с пространственным распределением параметров, описываемым функцией распределения Дхьхг), реализовано на плоской поверхности П. Если через область О распределения ФП проложить прямую ИЛ, чувствительную к воздействию данного поля по всей ее длине, полярные координаты которой задаются как (<р.у), то интегральный сигнал на ее выходе можно описать зависимостью

+00 +00

где 3(г) - дельта-функция Дирака, а з) - функция проекции.

Если индекс / определяет направление сканирования <р i (рис. 13), индекс7 определяет отсчет йу в выбранном / - м направлении, то паре индексов (/, /) соответствует прямая Ьу, вдоль которой сканируется область Г2, а значение проекции вдоль прямой ¿у можно записать в виде

г„ = Я/ {<Рп »«, )= |/(х„х2)<Н. (2)

ч

где К - преобразование Радона функции / а <// - приращение вдоль прямой Ьу.

При разбиении исследуемой области йсй2 на малые элементарные ячейки (считаем функцию / постоянной на и равной /¡¡) матрица/ соответствующая этому разбиению, и вектор-столбец неизвестных 7% примут вид:

/ =

/ А - Л' /- /.« - А.

= /г - Л /«, - Л. - и- (3)

1)+1 /.С-1Н2 ••• /«.,

Для выбранной схемы укладки ИЛ и соответствующей ей матрице А интегральные уравнения (1) преобразуются в систему вида:

АР= в. (4)

В правой части матричного уравнения (4) стоит столбец проекционных данных размера М=у(1)+...+ ){;')+...+ у(р), имеющий вид:

^ = (ЕП.....■■■'&!> ' ' ' > —. ёр^) ^, (5)

где КО _ число ИЛ в /-м направлении укладки.

Пары чисел (де, ^у), определяющие параллельную схему сканирования на плоскости (рис. 1д), могут быть рассмотрены как элементы группы С = Т х 1,

где интервал Т = [0, 2л) выступает в роли циклической группы относительно

сложения по модулю 2тг, а К - в роли аддитивной группы поля К..

В плоскости Оцк точки (<р„ образуют множество дискретизации. Для стандартной параллельной схемы укладки ИЛ в р направлениях с шагом /г матрица решетки дискретизации и двойственная к ней имеют диагональный вид:

'2Р

о %Г

ш - \?> 0 о л

Н1=11¥12п - решетки, порожденные матрицами V/ и V/1 над 2".

где Н =

Математически, задача состоит в том, чтобы восстановить функцию распределения параметра ФП /, по измерениям функции проекции g= 31/ на некоторой решетке.

Если функция / имеет ограниченную существенную ширину спектра, равную Ъ, то для ее восстановления с разрешением порядка 2п/Ъ, требуются стандартные регулярные решетки, для которых должны выполняться стандартные условия дискретизации, т.е. А < л/Ъ, р > Ъ [6].

Преобразование Фурье для функции проекции g = У?/ [6, 7], определим преобразование Фурье и обратное к нему равенствами

I (к, сг) = — ¡¡е(<Р, *)е^9*")Ж{1<р, где (к,а)

О к

Аппроксимация функции проекции g = 91/ на произвольной решетке Н определяется равенством: ~~~^д{у)Хк{2~у)> где г = (<р, 5)6 ГхМ, К

Руч1.

- подходящее компактное множество в ЖхЕ, обратное преобразование Фурье характеристической функции у Соответствующая теорема отсчетов дает оценку ошибки аппроксимации в виде [¿-(г) — ¿'¿■(г)!<тт"' j

(г*1<)\к

где г € ТхК, £ = (к, а) 6 ZxR, при условии, что трансляции множества К элементами решетки Н1 не пересекаются. Известно, что значение функции g = ЧЯ/ мало при 9: 0<5<1 вне множества

= ] (к,а) е Ъ х К: |сг| < Ь,\к\ < шах

9 ' 9

[6]. Таким образом, можно ут-

верждать, что множество К является «существенным» носителем функции g.

На практике часто встречаются схемы укладки измерительных линий, которым не соответствуют решетки, но их можно представить в виде объедине-

/

ния нескольких решеток или их сдвинутых копий: М= и(И,+ ■*/)• Так, например, для схемы укладки ИЛ в четырех направлениях (рис. 1в), имеет место равномерное распределение углов, но неравномерность отсчетов вдоль прямых. Эту схему сканирования можно рассматривать как объединение двух стандартных подсхем с двумя направлениями сканирования, т.е. Аф = я/2. При этом в плоскости 0<рз каждой из них соответствует своё множество дискретизации. Первой схеме соответствует множество дискретизации Л/;, а второй - . Очевидно, что множество М) является решеткой, т.е. М/ = Нь где Н] - подгруппа,

я/2

, т.е. Множество М2

порожденная над I? матрицей

0

решеткой не является, но его можно представить как смежный класс дискрет, (к! 2 0)

ной группы С по подгруппе И^Шг^ , = . Таким образом, объе-

I 0 М

динение М = М|иМ2 = ИlU{]Hl2+(я,/4 0)Г} представляет собой множество

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

Пусть Н: и Нг - произвольные решетки, Л/и Л;- фундаментальные области решеток Н,1 и Н2 соответственно. Пусть для имеет место включение Я\ С /?2, т.е. решетка Н2 является более плотной, чем Щ, а соответственно, и

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

Пусть К = /?2и(7?1+у) спектральное множество, такое, что Л21~| (Я, + у) = 0 и преобразование Фурье функции g = 91/ обращается в ноль вне множества К и пусть известны значения функции проекционных данных g на множествах дискретизации М\ = Ш\ и Мг = Н2 + XI. Предположим, что указанные множества дискретизации обладают следующим свойством:

функции gl и g2 , преобразования Фурье которых, обращаются в ноль вне множеств Л1 + V и Яг, соответственно, могут быть реконструированы от их выборок на множествах дискретизации М\ и М2 соответственно. Тогда основная идея, лежащая в основе представленного здесь подхода, следующая. Функцию я) ищем как разложение вида

8(<Р> = 8г(.Ч>> а)+Р{>Р, s)g^{<p, э), (б)

где Р(<р, 5) - неизвестная функция, обращающаяся в ноль всюду на М2, но не обращающаяся в ноль нигде на М\. Функции gl и g2 обращаются в ноль вне множеств и Я2 , соответственно. Разложение (6) можно использовать для реконструкции функции g от ее выборок на объединении множеств М\ и М2. Точно восстановление функции g(<p, я) описывается следующей теоремой.

Теорема. Пусть известны значения функции проекционных данных g = Щ на множествах дискретизации М\ и Мг, £ = 0 вне спектрального множества АГ = Л2и(й, + г). Пусть любая непрерывная функция И(х) е ¿2((?), преобразование Фурье которой к = 0 вне множества , может быть восстановлена по отсчетам на решетке М\. Тогда функция проекции g-<Rf может быть восстановлена по отсчетам на всем множестве М = Л/, и М2 по формуле:

*М = + (1 -^-'^с, Е ъ{*-у)> (7)

уШ, 1 е

где Хц- обратное преобразование Фурье характеристической функции Х/^

множества Su^Xlg{z) = CгY^g(x1+y)x^(z-x1~ у), для любогогвС,

уеш2

С, = ( 27Г)-'|^.|,/ = 1,2.

Из формулы (7) следует способ восстановления дискретных данных по их проекциям вдоль различных направлений сканирования.

Во-первых, для любого Иг е Н2 + и Уб Неочевидно, что (й2 - хг, V) = 0.

Следовательно, выражение (7) принимает вид g{h2) = ■ Это зна-

чит, что функция аппроксимации не меняет значения проекций вдоль указанных направлений сканирования. Таким образом, автоматически сформирована первая часть данных, которые обозначим как у)2.

Во-вторых, для любого А, € М1 = Н,, из формулы (7) следует, что

еМ = ^¿КМ + ^А (8)

где Д - некоторая корректирующая функция, которая вычисляется во втором слагаемом (7). Это значит, что проекционные данные, соответствующие множеству дискретизации М\ должны быть последовательно подвергнуты действию функции аппроксимации, а затем скорректированы. Обозначим через (5Я])\

результат аппроксимации проекционных данных, соответствующих множеству дискретизации М\\

(«/),(O=ÍHJ+x,g(0 = W|W21 Е ^(y)^ (í - y), Vi € Af, = Н,. (9)

veHj +x2

Из формулы (9), следует, что проекционные данные, соответствующие множеству дискретизации М\, приближаются значениями более плотной решетки Mi = Н2+Х2. Дальнейшая коррекция функции (УЯ j)\ заключается в добавлении к нему слагаемого вида

Е ~XK¡{z_y)> (10)

уе н, 1-е которое влияет только на g(y), у е М = Hi.

Заметим, что основной корректирующей частью функции 5Н А является разность g(y) -S„¡¡+x¡ g(y), y e Hi, стоящая в числителе под знаком суммирования и характеризующая погрешность аппроксимации данных на М\. Обозначим через (Щ/)з результат обработки проекционных данных, соответствующих множеству дискретизации М\, после аппроксимации и коррекции. На следующем этапе обработки происходит комбинирование измерительных (SR/)2 и ап-проксимационных данных Таким образом, описан следующий алгоритм UQC (union of quotient classes) для восстановления функции проекции g = 5R/ исследуемого параметра ФП.

Алгоритм UQC 1.

1. Создание массивов проекционных данных, соответствующих каждому из множества дискретизации;

2. Создание массива измерительных данных. Массив (9?/)2 формируется из проекционных данных по наиболее плотному множеству дискретизации Мг ~

Mi + х2:

3. Аппроксимация проекционных данных, соответствующих менее плотному множеству дискретизации М\ в соответствии с (9). Результат обозначается через (91/) i;

4. Коррекция проекционных данных(5К/)1 в соответствии с (10). Результат обозначается через (9t/)3;

5. Комбинирование измерительных (4Rj)i и аппроксимационных данных (9?/)з. Результат обозначается через 9?/

Аппроксимирующие функции (7) - (10), содержат значения функций вида Xr, где Xr - обратное преобразование Фурье характеристической функции

Xr множества R. Для того чтобы избежать сложного непосредственного вычисления интерполирующих функций вида Хт следует выполнить действия:

1. к функции SMg(x) применить преобразование Фурье, т.е. вычислить зна-

чениефункции (££)(*) = *,(£) Е ¿(4 +у)ег"М,

ve И1

2. затем применить формулу Пуассона [б]:

yzH VüH1

3. вычислить обратное преобразование Фурье.

Проблема неполноты данных возникает, когда значения функции s) недоступны по определенному диапазону аргументов (<р, s). В случае 2D томографии наиболее часто встречается проблема ограниченного набора углов. Если

известно значение аппроксимации SKg{% s), то можно вычислить значение Sa+Xg((p, s). Более того, если известно приближенное значение функции проекции g(<p,s) = 5?y{x,_v) на множестве дискретизации Н, то можно вычислить приближенное значение этой функции на множестве дискретизации И+х. Для этих целей использован алгоритм перехода от значений функции проекции g(<p,s) = Щху) на множестве М\ = Н) к значениям этой функции на множестве дискретизации Mj = ШЬ + хг, использующий формулу (4) и приведенный выше способ вычисления SMg[x) с помощью преобразований Фурье. При этом учитывается, что фундаментальные области R¡uR2 используемых решеток И^и , соответственно, имеют специфическую (непрямоугольную) форму и поэтому к ним нельзя применять стандартные Фурье-алгоритмы непосредственно без угрозы получения серьезных артефактов.

Опишем кратко технику, позволяющую восстановить функцию/ используя процедуру удвоения числа проекционных данных.

Алгоритм UOC 2. Удвоение числа проекционных данных

1. Создание массивов проекционных данных;

2. Создание массива измерительных данных по наиболее плотному множеству дискретизации;

3. Аппроксимация и коррекция проекционных данных по наименее плотному множеству дискретизации;

4. Комбинация измерительных и проекционных данных;

5. Создание массива дополнительных данных, используя сдвиги данных из п. 2 и 3 на половинные углы;

6. Комбинация данных из п. 4 и п.5.

На рис. 2 приведен сравнительный анализ восстановления эталонного двойного воздействия (а) на ВОИС с помощью стандартного метода БВР (б), а также с помощью алгоритмов ЦХ^С 1 (в) и ЦС>С 2 (г). Для каждого из представленных методов над реконструированной функцией указана нормированная среднеквадратичная ошибка (шБе) и коэффициент корреляции(согг).

Образцовое РВР, т<;е - 0.0068127 ирС1, тве = 0.0057664 иС)52, тве = 0.0048831

воздействие согг = 0.15274 согг = 0.13904 согг = 0.20656

Рис. 2. Сравнительный анализ восстановления функции распределения параметра ФП с помощью и(2С1,2 и РВР

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

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

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

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

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

В $ 3.2 проводится исследование особенностей применения НС типа однослойный персептрон для решения томографической задачи восстановления характеристик протяженных ФП по интегральным данным. Отмечается, что использование персептрона позволяет производить параллельную обработку информации за один проход от входа к выходу нейронной сети, что открывает перспективу обработки информации с высокой скоростью, вплоть до скорости света для оптической информации. Решена задача выбора оптимального параметра скорости обучения НС. Получена явная формула для вычисления этого параметра для первого этапа обучения. Для проведения численного анализа нейросетевого метода решения была рассмотрена томографическая задача восстановления информации с ВОИС повышенной размерности (п = 8) с укладкой информационных линий в трех направлениях, а также для ВОИС 3x3 с укладкой ИЛ в двух направлениях. Было показано, что результаты, полученные нейросетевым методом приблизительно те же, что в случае алгебраического моделирования методом Гревиля и для более сложных видов воздействия на ВОИС оба эти метода не дают высокого качества восстановления. Это говорит о том, что однослойные персептроны характеризуются ограниченными возможностями.

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

В § 3.4 представлены основные принципы обработки сигналов нерегулярных распределенных ВОИС томографического типа с помощью полнообразных нейронных ЛВР-сетей, обеспечивающих реконструкцию пространственных распределений параметров ФП.

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

решение задачи восстановления функции fix,у) ищется как выход нейронной

информационный массив, представленный вектором О, является набором томографических данных (4), (5), по которым нейронная сеть должна реконструировать вектор Р. Таким образом, нейросеть должна выполнить преобразование вида Р = А (С), предварительно пройдя обучение на множестве обучающих пар {(С, Р)}. Для создания обучающей страницы использовался способ выбора обучающих пар "с подкреплением", при котором рассматривались пары вида (Gi.Fi), связанные соотношением (?,.

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

_/ ч -Ж-^ЧУ-л,)2) Функции гауссовского типа ¿(х,у) = е 1 ' могут быть использованы в

качестве таковых, так как они имеют по своей природе локальный характер и принимают ненулевые значения только в зоне вокруг определённого центра. При создании обучающих пар ЯВИМЫ использовались функции гауссовского типа, причем параметры Ь„ с„ подбирались как узлы решетки соответствующей схемы укладки ИЛ.

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

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

Тип II. Это пары Гауссиан с непересекающимися носителями. Аналитически эти функции могут быть представлены в виде

при условии, что носители находятся на расстоянии не менее 2idb. Тип III. Эталонные распределения этого типа были получены с помощью процесса рандомизации с нормированием. Каждому целочисленному случайному набору аи аг, b\, ci, Ъг, сг ставилась в соответствие функция (11).

На рис. 3 приводятся результаты тестирования нейронной сети на функциях, как принадлежащих, так и не принадлежащих обучающей странице. В первых столбцах рис. 3 представлены оригиналы функций, а во вторых - результаты восстановления. Экспериментально показано, что данная сеть позволяет осуществлять восстановление функций пространственного распределения исследуемой физической величины с погрешностью в отдельной точке не более 1%. Ни один из приведенных для сравнения традиционных методов (FBP, ART) такого результата не дает.

сети в виде

Формируемый измерительной сетью

(П)

шее» 1.3329е-010 согг = 1

О О

тэе = 0.009461 согг = 0.98341

О 0

шее = 0.0036243 согг = 0.97997

0 0

газе= 0.020325 согг = 0.94654

тве = 7.3567е-011 согг = 1

0 0 5 тве= 0.0012959 согг = 0.98659

О О 5 ШБе = 0.0014293 согг = 0.97071

О 0 5

шее = 0.0060375 согг = 0.96452

Рис. 3. Тестирование 11ВР1Ш при п=5 (слева) и при п=10 (справа).

Построенная сеть ЯВРМЫ обладает хорошими прогнозирующими возможностями и может с успехом использоваться для распределенных вычислений в задачах с применением ВОЛС больших размерностей. Для получения точной информации о функционировании сложных систем, как например, в задачах о нахождении места и величины воздействия на ФП с помощью ВОИС, метод нейросетевой обработки информации ЯВЕМЫ вполне подходит. Однако необходимо отметить, что при этом способе восстановления в задачах с применением ВОИС повышенных размерностей возникают серьезные трудности при обучении сети ввиду слишком большого объема обучающих страниц. Поэтому возникает необходимость поиска оптимальных путей при применении нейронных сетей.

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

ментов и полученными путем нейросетевого моделирования. Представлен подход, позволяющий использовать RBFNN для предобработки проекционных данных с целью дальнейшего применения методов аппроксимации для нерегулярных схем укладки ИЛ более простого вида. В jj 4.1 представлен способ восстановления значений одного неизвестного направления сканирования. Обозначим через Ji, набор томографических данных по г'-му направлению сканирования: 1Я, =

(gn, ga, ■••• gu/iо)Т> 1 Предположим, что практически возможно получить значения данных только для р-1 направления сканирования, а для корректного восстановления функции классическими или специальными методами (например, ART, FBP , UQC 1, 2) необходимо иметь значения проекционных данных в р направлениях.

а) Образцовое воздействие q 5

.......Т... О

А А

д) mse = corr =

- 0.01339 0.56605

е) mse = 0.011061 corr - 0.67981

1

"" Ю О

5 и

0 0 I

Рис. 4а. Комбинированный нейро-алгебраический метод восстановления функции, принадлежащей ОС: а) Эталон; б) - г) восстановление одной проекции и график ошибки по известным значениям трех проекций: б) оригинал проекции, в) Г-Ш РЫМ-восстановление проекции, г) ошибка восстановления проекции; д) восстановление методом РВР для р=4; е) восстановление с помощью нейросетевой генерации одной проекции и дальнейшего применения алгебраического метода и<ЗС 2.

Рис. 46. Комбинированный нейро-алгебраический метод восстановления функции, не принадлежащей ОС (обозначения аналогичны используемым для рис. 4а).

Для решения этой задачи предложено использовать НС радиально-базисного типа, которая должна выполнить преобразование вида: 9?,, (Л2, ,.., Для проведения анализа нейросетевого метода решения

Я) mse = 0.019946 corr - 0.56839

е) mse = 0.014172 corr = 0.74112

задачи с помощью ЯВЕМ^ были рассмотрены томографические задачи восстановления информации с ВОИС повышенной размерности. При ЯВР-восстановлении проекции для функции, принадлежащей обучающей странице (или близкой к функции из ОС), ошибка равна нулю (рис. 4а). Для функции, не принадлежащей обучающей странице (рис. 4 б), значения rn.se не нулевые. Однако следует заметить, что проекция, построенная с помощью НС похожа на исходную (рис.4 б), но имеет более сглаженный вид, что вполне соответствует общей концепции обучаемости НС. В результате применения НС происходит восстановление проекционных данных с недостающих направлений сканирования, что предоставляет возможность применения методов аппроксимации данных на нерегулярной сетке специального вида. Предложен и исследован новый комбинированный алгоритм обработки проекционных данных, сочетающий в себе нейросетевую генерацию проекций с последующим применением методов алгебраической аппроксимации на нерегулярной сетке, представляющей собой объединение классов смежности. Применение комбинированного алгоритма позволяет восстанавливать информацию с полей вдвое более высоких размерностей по сравнению с методами, использованными в главе 3.

В § 4.2 представлен способ удвоения числа направлений сканирования. Предположим, что для корректного восстановления функции / необходимо иметь значения проекционных данных в р направлениях 9?,, У12, 9? , а

известны значения данных только для половины от необходимого числа направления сканирования Л,, 9?3, ...,

Рис. 5. а) Эталон, Ь) иС)С 2 - окончательное восстановление эталонной функции после применения нейро-алгебраического метода, с) оригинал проекции № 2 (под углом 45°), с1) восстановление проекции № 2 с помощью нейронной сети КВК№>1, е) ошибка восстановления функции проекции второго направления, /) оригинал проекции № 4 (под углом 135°), восстановление проекции № 4 с помощью нейронной сети КВРК\Г, И) ошибка восстановления функции проекции четвертого направления.

Для решения этой задачи НС должна выполнить преобразование вида 5Н,, ..., 5Я4, ..., На рис.5 приводится исследование при р=4

для эталонной функции из ОС (слева) и для функции, которая незначительно отличается от функций из ОС (справа).

Очевидно, что два направления укладки требуют значительного (как минимум в два раза) увеличения объёма обучающей страницы по сравнению с тремя направлениями, 2) существенно снижается прогнозирующая способность нейронной сети МЛ-ТСН по сравнению со случаем укладки в трех направлениях. Все эти недостатки являются следствием некорректности поставленной задачи. Одним из способов для решения возникших проблем является повышение вычислительной мощность НС в сочетании с оптимизацией вычислений.

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

Глава 5 посвящена оптимизации вычислительного процесса при обработке информации РВОИС, на основе применения специальных алгоритмов структурирования массивов данных, поступающих с ВОИС, созданию новых нейросетевых конструкций - комплексов нейронных сетей (КНС) и новых комбинированных алгоритмов обработки данных, сочетающих в себе специализированные КНС, аналитические и алгебраические методы.

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

следующих трех процессов: предобработки данных с целью оптимизации геометрии ВОИС, обработки нейросетевыми или алгебраическими методами, восстановления первоначальной геометрии ВОИС. В данном контексте оптимизация геометрии ВОИС заключается в ее "обрезке", т.е. в удалении строк и столбцов по краям матрицы / (3), сумма элементов которых имеет значения, равные (или близкие) нулю. Далее модифицируются сама матрица/ и столбец проекционных данных (5).

эталон

а) б) в) г)

Рис. 6. а) Исследуемая функция, 5), в) проекции вертикального и горизонтального направлений, г) подобласть после оптимизации геометрии ("обрезки") измерительной сети.

Таким образом, при выполнении описанного выше алгоритма производится выделение областей-кандидатов, в которых находятся искомые "объекты" (рис.бг). Эти области в подавляющем большинстве случаев имеют меньшую размерность, чем исходная область. На следующем этапе информация, поступающая с ранее оптимизированной измерительной сети, обрабатывается с помощью как обычных традиционных процедур восстановления, таких как FBP, ART, так и специальных разработанных автором алгебраических алгоритмов UQC (рис. 7), а также специализированных нейросетевых методов.

UQC без предобработки UQC с предобработкой

Образцовое воздействие mse =0.0051873, corr = 0.33157 mse =0.00077268, corr = 0.89099

Рис. 7. Восстановление функции распределения исследуемого параметра ФП с помощью алгебраического алгоритма и()С без предобработки и с предобработкой.

проекция №1

проекция №2

подобласть

30

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

Особое внимание в главе 5 уделено созданию и исследованию специализированных нейросетевых конструкций. Так в §5.4 первые введено и исследовано понятие комплекса нейронных сетей, который представляет собой множество, состоящее из нескольких заранее обученных НС:

SN= {NN(nh mi),NN(n„ mi),..., NN(nK, mK)}, (12)

ni< ...<ni<...<nic,mj< ...<ти,<... < т&,

(nh mi) Ф (nj, mj) при / tj, 1< i,j < K. (13)

Для обработки проекционных данных поступающих от ВОИС размера п'хт', п'<п, т'<т в SN вида (12) выбирается НС вида NN(nh mi), для которой

я/./ < п < п„ т,.!< т<т„ 1 <i,j<K. (14)

Из условий (13), (14) следует, что нейронная сеть NNfa, mi) является сетью наименьшей размерности в SN, с помощью которой можно обработать проекционные данные измерительной сети размера п'хт'.

В данной работе используются три различных комплекса НС 5.5-5.7):

1) SNf = { NNf(nh mi), 1< i < К} - комплекс НС для восстановления параметров ФП. Здесь через NN/(n,, mi) обозначается НС для восстановления параметра ФП по проекционным данным, поступающим с ВОИС размера и,х тя,-.

Если Kj - число ИЛ в соответствующей ВОИС, то сеть NNf{nh mi) выполняет преобразование пространств: Ж*' —> М"'™1.

2) SNg - комплекс НС для восстановления функции проекции. Через NNg(ni, mi) обозначается НС для восстановления функции проекции g(<p, s), а SNg = {NNg(n„ mi), 1< / < К}.Нейронная сеть NNg(n„ mi) выполняет преобразование пространств: Rw -> Мг, где М=у{ 1)+...+ y(j)+-.+ у(р) - общее число ИЛ по всем р известным направлениям сканирования, у - число ИЛ в восстанавливаемом направлении укладки ИЛ.

Комплекс нейронных сетей SNg может использоваться, например, для того, чтобы сделать неполные проекционные данные равномерными по углу или для восстановления синограммы в случае критически малого (р= 2, 3) числа направлений укладки ИЛ.

3) SNj - комплекс НС для стандартизации проекционных данных. Через NNdjij, mi) обозначается НС для преобразования к необходимому виду входных данных в ВОИС, а через SNj обозначается множество из К обученных НС: SNj = {NN^fii, mi), 1 < i < К }. НС NNJirti, т,) выполняет преобразование пространств: —» R^', где y(J)- первоначальное число ИЛ в j-м направлении, а fiij) - требуемое число ИЛ в_/-м направлении.

Комплекс 5'N¡1 может использоваться, например, для того, чтобы сделать неполные проекционные данные равномерными по отсчетам, в частности, для восстановления информации с поврежденной ИЛ.

В §§ 5.8-5.9 представлены новые комбинированные алгоритмы обработки информации, поступающей с распределенных измерительных систем. Представлены различные схемы восстановления функций ФП, сочетающие в себе нейросетевые, аналитические и алгебраические методы. Все указанные комбинированные алгоритмы содержат общую часть - нейросетевую генерацию проекций в тех направлениях, где они отсутствуют. Для применения классического метода РВР этого недостаточно, т.к. данные остаются неравномерными по отсчетам. Для регуляризации по отсчетам применяется комплекс в соответствии со схемой, представленной на рис. 8б.

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

Важной задачей вычислительной томографии является уменьшение времени реконструкции объекта контроля. Поэтому далее вместо затратных по времени классических методов обратного проецирования можно использовать нейросетевое восстановление функции по проекциям (рис. 8г). В соответствии со схемой, представленной на рис. 8г, сначала проекционные данные обрабатываются комплексом нейронных сетей 57У8 , а затем - комплексом &\Гу. При этом проекционные данные сначала приобретают регулярный по углу вид, а затем по восполненным проекционным данным производится восстановление распределения исследуемого параметра ФП.

о) Образцовое 6) РВР, Ш5е = 0.0086244 «) 1X58, пас = 0.0079531 г) 8>|, шее» 0.0071258

воздействие сшт = 0.27348 согг = 0.31233 согг= 0.50508

Рис. 8. а) Образцовое воздействие на измерительную сеть; б) - г) обработка проекционных данных в соответствии с изображенными внизу схемами.

На рис. 9 в виде графиков приведен сравнительный анализ методов, предложенных в главе 5. На рис. 9а-в приведены графики зависимостей нормированных среднеквадратичных ошибок (rn.se) от величины п, определяющей

размер ВОИС в случае одного, двух и трех воздействий на ВОИС для методов, предложенных в главе 5 и схематически изображенных на рис. 8. Первый из предложенных методов - классический алгоритм РВР при р=4. Следующие алгоритмы были применены после обработки комплексом нейронных сетей для восстановления параметра ФП на нерегулярной сетке без удвоения числа проекционных данных (ЪТ<ЗС 1) и с удвоением (ЬТ()С 2). И наконец, нейросете-вой метод представляет собой последовательное применение к неполным нерегулярным данным двух нейросетевых комплексов и Рис. 9 г отражает зависимость времени восстановления параметра ФП от размера ВОИС для конкурирующих между собой нейросетевого и алгебраического методов.

а) Одно воздействие

б) Два воздействия

ф «

Е

10 15 20 25

5 < п < 30 в) Три воздействия

10 15 20 25 30 5 < п < 30 г) Время восстановления функции ФП 0.8

—Ч— ШС2 БЛ/ ,

10 15 20 25 5 < п < 30

10 15 20 25 30 5 < п < 30

Рис.9. Графики зависимостей ошибок (тве) при восстановлении одного (а), двух (б) и трех (в) воздействий на измерительную сеть и времени восстановления параметров ФП (г) от размерности ВОИС.

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Разработан и исследован новый метод реконструктивной томографии для распределенных ВОИС, представляющих собой набор волоконно-оптических измерительных линий, уложенных в соответствии с определенной схемой на исследуемом объекте. Метод основан на максимально обобщенных теоремах дискретизации на теоретико-групповой основе. Он может быть применен для распределенных ВОИС томографического типа в тех случаях, когда не выполняются классические условия корректности выборки, т.е. когда нарушается равномерность схемы сканирования по отсчетам вдоль выбранных направлений. Представлен способ для аппроксимации функции проекции на нерегулярной сетке, для случая, когда множество дискретизации представляет собой объединение нескольких различных решеток.

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

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

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

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

6. Разработаны основные модели и методы применения нейросетевых технологий обработки данных распределенных ВОИС томографического типа, обеспечивающих реконструкцию пространственных распределений параметров ФП в реальном времени для полнообразных НС радиально-базисного типа. При реализации 11ВР1ЧН предложен новый эффективный способ создания функций для обучения нейронной сети, которые строятся на основе тех функций, которые моделируют конкретный физический процесс.

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

8. Предложен новый комбинированный нейро-алгебраический алгоритм обработки проекционных данных для реконструкции функций распределения физических полей. Этот алгоритм заключается в последовательном выполнении двух процессов: нейросетевой предобработки измерительной информации с помощью ЯВРЬШ и дальнейшей аппроксимации полученных данных на множестве дискретизации специального вида - на объединении нескольких решеток или их сдвинутых копий.

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

10. Впервые введено и исследовано понятие специализированной нейросетевой конструкции - комплекса нейронных сетей. Комплекс НС представляет собой множество, состоящее из нескольких заранее обученных нейронных сетей. Каждая НС, принадлежащая комплексу, предназначена для специальной обработки данных, поступающих с ВОИС соответствующего размера. Для обработки проекционных данных поступающих от ВОИС среди таких НС выбирается сеть наименьшей размерности,

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

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

ЦИТИРУЕМАЯ ЛИТЕРАТУРА

1. Кульчин Ю.Н. Распределенные волоконно-оптические измерительные системы. М. : Физматлит. - 2001. - 272 С.

2. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С. Многомерная обработка сигналов с использованием волоконно-оптической измерительной сети // Квантовая электроника. - 1993. - Т. 20, №5.- стр. 711-714.

3. Кульчин Ю.Н., Витрик О.Б., Ромашко Р.В., Петров Ю.С., Кириченко О.В., Каменев О.Т. Томографические методы для исследования векторных полей при помощи волоконно-оптических измерительных систем// Квантовая электроника. - 1997.- Т. 24, №5.- стр.46.

4. Petersen D.P., Middleton D. Sampling and reconstruction of wave-number-limited functions in N-dimensional Euclidean space // Inf. Control. -1962.-Vol. 5. -pp. 279-323.

5. Kluv'anek I. Sampling theorem in abstract harmonic analysis // Mat. Casopis Sloven. Akad. Vied. - 1965. - Vol. 15. - pp. 43-48.

6. Natterer, F. The Mathematics of Computerized Tomography. Stuttgart : B.G. Teubner and John Wiley and Sons. - 1986.

7. Faridani A. Introduction to the mathematics of computed tomography. Inverse Problems and Applications. - 2003. - pp. 1-46

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Kulchin Yu.N., Denisov I.V., Denisova E.V., Milovanov V.I. Application of algebraic methods for restitution of cumulative distribution functions of physical fields // Proc. of International Workshop on Optical Beam Transformation. - IWBT, Vladivostok. - 2001.

2. Kulchin Yu.N., Denisov I.V., Denisova E.V., Sedov V.A., Drozdov R.S. Fiber-Optic Distributed Measuring System with Neural Network Technology of Signal Array Processing // International workshop on Bio-Signal Analysis and its Applications. - Tokyo, Japan. - 2002. - Vol. 2. - pp. 3-6.

3. Денисов И. В., Денисова Е. В. Применение приближенных алгебраических и нейросетевых методов решения томографической задачи // Электронный журнал «Исследовано в России». - 2002. - Т. 201.

- С. 2222-2228.

4. Kulchin Yu. N., Denisova E.V., Denisov I. V., Efimov M.B. Selection of optimal parameter of speed of training of neural network perceptron type // Fundamental Problems of Optoelectronics and Microelectronics, Proceedings of the SPIE. - 2003. - Vol. 5129. - pp. 162-167.

5. Кульчин Ю.Н., Денисов И.В., Денисова E.B. Принципы организации матриц связей оптической нейронной сети на объемных оптических элементах // Нейрокомпьютеры: разработка, применение.1 - 2003. - № 7.-С. 21-27.

6. Kulchin Yu. N., Denisova E.V., Denisov I. V. Synthesis of approximate algebraic and neural-like methods for the solution of the tomography problem //Proc. of SPIE.-2003.-Vol. 5134.-pp. 83-89.

7. Denisov I., Kulchin Y., Kirichenko O., Sedov V., Drozdov R., Denisova E. Model of optoelectronic measuring intelligent system // Proc. of International Conference Physics and Control. - 2003. - Vol. 1. - pp. 172-175.

8. Кульчин Ю. H, Денисова E. В. и др. Макет оптоэлектронной нейроподобной измерительной системы // Нано- и микросистемная техника2. - 2003. - № Ю. - стр. 40 - 42.- ISSN 1813-8586.

9. Кульчин Ю. Н., Денисов И. В., Седов В. А., Кириченко О. В., Денисова Е. В., Дроздов Р. С. Принцип организации волоконно-оптической измерительной сети повышенной размерности // Оптико-электронные информационно-энергетические технологии. - Украина, Винница - 2003.

- № 2 (4). - С. 205-211.

Ю.Денисов И.В., Денисова Е.В., Рыбальченко Н.А. Восстановление двоичных воздействий на волоконно-оптической сети произвольной размерности // Материалы XLVI Всероссийской межвузовской научно-технической конференции. - 2003. - Т. 1. - С. 63-65.

' Жирным шрифтом выделены издания, входящие в перечень ВАК РФ

2 Журнал «Микросистемная техника» с 1999 г. по 2004 г.

11. Kulchin Yu. N., Denisov I. V., Denisova E. V., Piskunov E. N. Prismatic neural chip for distributed measuring networks // Optical Memory & Neural Networks (Information Optics). - 2003. - Vol. 12, № 3. - pp. 237-242.

12. Kulchin Yu. N., Denisova E. V., Denisov I. V. Application of algebraic and neural-like methods for reconstruction of distribution functions of physical fields II Optical Memory & Neural Networks (Information Optics). - 2003. - Vol. 12, № 4. - pp. 283-297.

13. Kulchin Y. N., Denisova E.V., Denisov I. V., Rybalchenko N.A. Principal of reconstruction of the single influences on fiber-optical measuring network // Pacific Science Review. - 2003. - Vol. 5 (1). - pp. 32-37.

14. Denisov I. V., Denisova E.V., Rybalchenko N.A.. Reconstruction of twice influences on the fiber-optical measuring network // Труды II Международной конференции «Параллельные вычисления и задачи управления». - 2004. - С. 531-534.

15. Denisov I.V., Denisova E.V., Rybalchenko N.A. Principal of reconstruction of the ternary influences on the fiber-optical measuring network // Proc. of Fourth Asia-Pacific Conference "Fundamental Problems of Optoelectronics and Microelectronics". - 2004. - Vol. 1. - pp. 238-240.

16. Гридин A.A., Денисова E.B., Денисов И.В. Обработка информации с распределенного температурного поля // Proc. of 7-th International Conference Fundamental Problems of Optoelectronics and Microelectronics , Moscow , Russia. - 2005. - pp. 365-368.

17. Denisov I.V., Denisova E.V., Rybalchenko N.A., Sedov V.A. Definition of the fiber-optical tomography problem // Proc. of International Conference Physics and Control. - 2005. - Vol. 1. - pp. 826 - 829.

18. Zakasovskaya E.V., Fadeev V.V. Restoration of point influences by the fiber-optical network in view of a priori information // Proc. of APCOM. -2007. - CD.

19. Kulchin Yu.N., Zakasovskaya E.V. Application of Radial Basis Function Neural Network for Information Processing in Fiber Optical Distributed Measuring Systems // Optical Memory & Neural Networks (Information Optics). - 2008. - Vol. 17, № 4. - pp. 317-327.

20. Кульчин Ю.Н., Закасовская E.B. Восстановление информации волоконно-оптическими измерительными системами с использованием радиально-базисных нейронных сетей // Нейроинформатика-2009. - Сб. науч. тр. М.: МИФИ.-2009.- С. 289-298.

21. Kulchin Yu. N., Zakasovskaya E. V. Artifacts suppression in limited data problem for parallel fiber optical measuring systems // Optical Memory & Neural Networks (Information Optics). - 2009. - Vol.18, № 3. -pp.171-180.

22. Кульчин Ю.Н., Закасовская E.B. Неравномерные схемы укладки измерительных линий в распределенных волокнно-оптических системах // Информатика и системы управления. - 2009. - № 3(21) . - С. 61-71.

23.Кульчин Ю.Н., Закасовская Е.В. Нейросетевое и алгебраическое моделирование параллельного 2d проецирования в волоконно-оптической томографии при ограниченном числе направлениий сканирования // Компьютерная оптика. - 2009. - Т. 33, № 3. - С. 318- 324.

24. Kulchin Yu.N., Zakasovskaya E.V. Artifacts suppression in limited data problem for parallel fiber optical measuring systems // Proc. of APCOM -

2009. - CD.

25.Кульчин Ю.Н., Закасовская Е.В. Нейросетевое подавление артефактов при малом числе направлений укладки измерительных линий в волоконно-оптических системах // Нейроинформатика-2010 , Сб. науч. тр. МИФИ.-2010.-часть 1. - С. 110-120.

26. Кульчин Ю.Н., Закасовская Е.В. Моделирование параллельного 2D-проецирования в волоконно-оптической томографии для малого числа направлений сканирования II Информатика и системы управления. -

2010.-Т. 1(23).-С. 104-114.

27. Kulchin Yu. N., Zakasovskaya E. V. Optimizing algebraic and neural methods for information processing in distributed fiber-optical measuring systems // Optical Memory & Neural Networks (Information Optics). -2010. -Vol. 19, № 3. - pp. 237-247.

28. Кульчин Ю.Н., Закасовская Е.В. Обработка информации комплексом нейронных сетей в распределенных волоконно-оптических измерительных системах // Компьютерная оптика. - 2010. - Т. 34, № 3. -С. 363 -369.

29. Kulchin Yu.N., Zakasovskaya E.V. Neural-like specialized construction for fiber optical distributed measuring systems in limited-data conditions // Pacific Science Review. - 2010. - Vol. 12 (1). - pp. 80-86.

30. Kulchin Yu.N., Zakasovskaya E.V. Complexes of neural networks for information processing in distributed fiber-optical measuring systems // The First Russia and Pacific Conference on Computer Technology and Applications" (RPC 2010). -Vladivostok, Russia. - 2010. - C. 125-129.

31. Кульчин Ю.Н., Закасовская Е.В. Оптимизация алгоритма обработки информации в распределенных волоконно-оптических измерительных системах // Информатика и системы управления. - 2010. - № 4 (26). -С. 50 - 60 . - ISSN 1814-2400.

32. Кульчин Ю.Н., Закасовская Е.В. Нейросетевые конструкции для распределенных волоконно-оптических измерительных систем // Нейроинформатика-2011, Сб. науч. тр. МИФИ. -2011. - ч. 2. - С. 171-180.

33. Кульчин Ю.Н., Закасовская Е.В. Комплекс нейронных сетей для синтеза синограмм в распределенных волоконно-оптических системах // Нейрокомпьютеры: разработка, применение. - 2011. - № 7. - С. 55-62.

Закасовская Елена Владимировна

НЕЙРОСЕТЕВОЕ И АЛГЕБРАИЧЕСКОЕ МОДЕЛИРОВАНИЕ В РАСПРЕДЕЛЕННЫХ ВОЛОКОННО-ОПТИЧЕСКИХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМАХ В УСЛОВИЯХ НЕПОЛНОТЫ ДАННЫХ

Автореферат

Подписано к печати 21.03.2011 Усл.п.л. 2,0. Уч.-изд.л. 1,7

Формат 60 х 84/16 Тираж 100 экз. Заказ 17

Издано ИАПУ ДВО РАН, г. Владивосток, ул. Радио, 5 Отпечатано участком оперативной печати ИАПУ ДВО РАН г. Владивосток, ул. Радио, 5

Оглавление автор диссертации — доктора технических наук Закасовская, Елена Владимировна

1. СОВРЕМЕННЫЕ ТЕНДЕНЦИИ ОРГАНИЗАЦИИ

РАСПРЕДЕЛЕННЫХ ВОЛОКОННО-ОПТИЧЕСКИХ ИЗМЕРИТЕЛЬНЫХ

СИСТЕМ.

1.1. Распределенные волоконно-оптические измерительные сети томографического типа.

1.1.1. Распределенные волоконно-оптические измерительные системы.

1.1.2. Оптоэлектронная распределенная сигнальная система.

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

1.1.4. Применение нейронных сетей в РВОИС.

1.2. ЦЕЛЬ И ЗАДАЧИ РАБОТЫ. ОСНОВНЫЕ ЗАЩИЩАЕМЫЕ ПОЛОЖЕНИЯ.

1.3. ОБЗОР ПО

ГЛАВАМ.

2. СОВРЕМЕННЫЕ АЛГЕБРАИЧЕСКИЕ МЕТОДЫ ОБРАБОТКИ ДАННЫХ В ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫХ СИСТЕМАХ.

2.1. ОБОЗНАЧЕНИЯ И ОПРЕДЕЛЕНИЯ.

2.1.1. Стандартная параллельная схема сканирования.

2.1.2. Множества дискретизации для параллельной схемы сканирования.

2.2. ОБОБЩЕННЫЕ ТЕОРЕМЫ ДИСКРЕТИЗАЦИИ.

2.2.1. Классическая теорема отсчетов.

2.2.2. Теоретико-групповая терминология.

2.2.3. Свойства групп С и С.

2.2.4. Фундаментальные ^бласти в п-мерном Евклидовом пространстве.

2.2.5. Обобщение теоремы Шеннона на Мп.

2.2.6. Обобщение теоремы Шеннона на случай локально компактной абелевой группы.

2.3. НОСИТЕЛЬ СПЕКТРА ФУНКЦИИ.

2.3.1. Преобразование Фурье функции проекции.

2.3.2. Носитель спектра функции проекции.

2.3.3. Неравномерная схема укладки ИЛ в четырех направлениях.

2.3.4. Восстановление функции проекции на объединении двух множеств дискретизации.

2.3.5. Основная теорема.

2.4. АНАЛИЗ АППРОКСИМАЦИОННОЙ ФОРМУЛЫ И АЛГОРИТМ ВОССТАНОВЛЕНИЯ ФУНКЦИИ ПО ПРОЕКЦИЯМ.

2.4.1. Алгоритм I (и<ЗС 1). Восстановление функции проекций в соответствии с формулой аппроксимации (2.11).

2.4.2. Пример комбинирования измерительных и аппроксимационных данных для схемы укладки измерительных линий в четырех направлениях.

2.4.3. Вычисление аппроксимирующей функции.

2.4.4. Пополнение множества проекций.

2.4.5. Алгоритм II (ХЛ^С 2). Удвоение числа проекционных данных.

2.4.6. Пример схемы получения дополнительных проекций для четырех направлений сканирования.

ВЫВОДЫ.

3. НЕЙРОСЕТЕВЫЕ ПРИНЦИПЫ ОБРАБОТКИ ИНФОРМАЦИИ РАСПРЕДЕЛЕННЫМИ ВОИС. НЕЙРОННЫЕ СЕТИ ПОЛНОГО ОБР АЗА.

3.1. ОБЩИЕ ПРИНЦИПЫ ИСПОЛЬЗОВАНИЯ НС В РВОИС.

3.1.1. Применение нейросистем в РВОИС.

3.1.2. Дискретизация области.

3.2. НЕЙРОННЫЕ СЕТИ ТИПА ОДНОСЛОЙНЫЙ ПЕРСЕПТРОН.

3.2.1. Принципы организации НС типа однослойный персептрон.

3.2.2. Строение однослойного персептрона.

3.2.3. Обучение однослойного персептрона.

3.2.4. Обучение персептрона на множестве изр обучающих пар.

3.2.5. Выбор оптимальньного параметра скорости обучения.

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

3.2.7. Численное моделирование для однослойного персептрона* в случае укладки измерительных линий в двух направлениях.

3.3. ПРИНЦИПЫ ОРГАНИЗАЦИИ НЕЙРОННОЙ СЕТИ ТИПА МНОГОСЛОЙНЫЙ ПЕРСЕПТРОН.

3.4. ИСПОЛЬЗОВАНИЕ НЕЙРОННЫХ СЕТЕЙ РАДИАЛЬНО

БАЗИСНОГО ТИПА.

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

3.4.2. Радиальные функции.

3.4.3. Радиально-базисные функции с компактными носителями.

3.4.4. Радиальные функции гауссовского типа.

3.4.5. Алгоритм построения обучающей страницы.

3.4.6. Эталоны для создания обучающей страницы.

3.4.7. Тестирование нейронной сети радиально-базисного типа.

3.4.8. Сравнительный анализ.

ВЫВОДЫ.

4. ИСПОЛЬЗОВАНИЕ НЕЙРОННЫХ СЕТЕЙ ДЛЯ ПРЕДОБРАБОТКИ ПРОЕКЦИОННЫХ ДАННЫХ.

4.1. ДООПРЕДЕЛЕНИЕ ЗНАЧЕНИЙ ОДНОЙ ПРОЕКЦИИ.

4.1.1. Постановка задачи.

4.1.2. Построение обучающей страницы.

4.1.3. Восстановление функции проекции на объединении двух множеств дискретизации.

4.1.4. Пополнение множества проекционных данных для случая укладки измерительных линий в трех направлениях.

4.1.5. Эталонные распределения для ВОИС при п=10.

4.1.6. Выбор оптимального числа функций в обучающей странице.

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

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

4.1.9. Комбинированный метод восстановления функций ФП.

4.2. УДВОЕНИЕ ЧИСЛА ПРОЕКЦИЙ С ПОМОЩЬЮ НЕЙРОСЕТЕВОГО МЕТОДА.

4.2.1. Случай укладки ИЛ в четырех направлениях.

4.2.2. Эталоны для создания обучающей страницы.

4.2.3. Тестирование НС радиально-базисного типа при п=Т0.

ВЫВОДЫ.

5. ОПТИМИЗАЦИЯ НЕЙРОСЕТЕВЫХ, АЛГЕБРАИЧЕСКИХ И КОМБИНИРОВАННЫХ АЛГОРТИМОВ ОБРАБОТКИ ИНФОРМАЦИИ, ПОСТУПАЮЩЕЙ С РАСПРЕДЕЛЕННЫХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМ

5.1. МАТРИЦА СХЕМЫ СКАНИРОВАНИЯ.

5.2. ОПТИМИЗАЦИЯ ГЕОМЕТРИИ ВОИС.

5.3. АЛГЕБРАИЧЕСКИЕ АЛГОРИТМЫ и<)С И БМирС.

5.4. ОСНОВНЫЕ НЕЙРОСЕТЕВЫЕ КОНСТРУКЦИИ.

5.5. КОМПЛЕКС НЕЙРОННЫХ СЕТЕЙ ДЛЯ ВОССТАНОВЛЕНИЯ ПАРАМЕТРА ФИЗИЧЕСКОГО ПОЛЯ.

5.6. КОМПЛЕКС НС ДЛЯ ВОССТАНОВЛЕНИЯ ФУНКЦИИ ПРОЕКЦИИ

5.7. КОМПЛЕКС НЕЙРОНННЫХ СЕТЕЙ ДЛЯ СТАНДАРТИЗАЦИИ ПРОЕКЦИОННЫХ ДАННЫХ.

5.8. КОМБИНИРОВАННЫЕ МЕТОДЫ С ПРИМЕНЕНИЕМ КОМПЛЕКСОВ НЕЙРОННЫХ СЕТЕЙ И АЛГЕБРАИЧЕСКОГО СИНТЕЗА СИНОГРАММ.

5.9. РЕЗУЛЬТАТЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ.

ВЫВОДЫ.

Заключение диссертация на тему "Нейросетевое и алгебраическое моделирование в распределенных волоконно-оптических измерительных системах в условиях неполноты данных"

Ее основные результаты, сводятся в основном к следующему:

1. Разработан и исследован новый метод реконструктивной томографии для распределенных ВОИС, представляющих собой набор волоконно-оптических измерительных линий, уложенных в соответствии с определенной схемой на исследуемом объекте. Метод основан на максимально обобщенных теоремах дискретизации на теоретико-групповой основе. Он может быть применен для распределенных ВОИС томографического типа в тех случаях, когда не выполняются классические условия корректности выборки, т.е. когда нарушается равномерность схемы сканирования по отсчетам вдоль выбранных направлений. Представлен способ для аппроксимации функции проекции на нерегулярной сетке, для случая, когда множество дискретизации представляет собой объединение нескольких различных решеток.

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

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

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

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

6. Разработаны основные принципы применения нейросетевых методов обработки сигналов распределенных ВОИС томографического типа, обеспечивающих реконструкцию пространственных распределений параметров ФП в реальном времени для полнообразных НС радиально-базисного типа. При реализации КВРЫЫ предложен новый эффективный способ создания функций для обучения нейронной сети, которые строятся на основе тех функций, которые моделируют конкретный физический процесс.

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

8. Предложен новый комбинированный нейро-алгебраический алгоритм обработки проекционных данных для реконструкции функций распределения физических полей. Этот алгоритм заключается в последовательном выполнении двух процессов: нейросетевой предобработки измерительной информации с помощью КВИчПЧ и дальнейшей аппроксимации полученных данных на множестве дискретизации специального вида - на объединении нескольких решеток или их сдвинутых копий.

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

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

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

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

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

ЗАКЛЮЧЕНИЕ

Решение задачи оценки параметров распределенных ФП, подверженных априорно неизвестным внешним воздействиям, имеет исключительно большое значение. Особое значение эта задача имеет для контроля состояния структурных элементов технических и инженерных конструкций, в случаях когда необходим постоянный мониторинг крупномасштабных и пространственно неоднородных ФП. В связи с этим возникает необходимость создания специализированных измерительных средств, обладающих распределенной в пространстве чувствительностью, способных осуществлять измерение параметров различных ФП на значительных площадях пространства, а также обработку формируемых при этом многомерных информационных массивов в реальном масштабе, времени. Широкие возможности в подобных задачах открывает использование измерительных систем на основе распределенных ВОИС. В таких измерительных системах оптические сигналы, поступающие по волоконно-оптическим измерительным линиям, изменяются пропорционально величине внешнего физического воздействия. Сигнал на выходе измерительной линии представляет собой линейный интеграл от функции распределения исследуемого параметра физического поля. Таким образом, возникает задача восстановления функций распределения физических полей по их интегральным характеристикам. Математически — это томографическая задача.

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

Интегральный образ обладает сильной дискретизацией отсчетов не только по углу, так и по числу проекций.

В таких измерительных системах зачастую используются нестандартные схемы укладки измерительных линий, в том числе нерегулярные, а также имеет место сложность геометрии исследуемых областей. Следствием «протяженности» ВОИС является большой объем обрабатываемых данных, приемлемые алгоритмы реконструкции ФП являются итерационными, что не позволяет создавать быстродействующие информационно-измерительные системы даже при применении параллельных принципов организации вычислительных сетей.

Эти и другие особенности накладывают серьезные ограничения на возможные применения РВОИС.

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

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

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

Представлена явная формула для получения дополнительных значений проекционных данных (отсутствующих или искаженных), использующая аппроксимацию функции проекции на нерегулярной сетке. В дальнейшем к этим проекционным данным можно применять стандартные алгоритмы восстановления. Предложены два метода и(^С 1, иС>С 2 восстановления проекционных данных. Приведен сравнительный анализ эффективности работы и<ЗС 1, и(^С 2 со стандартным аналитическим методом волоконно-оптической томографической реконструкции БВР.

Представленные алгоритмы иС^С 1,2 позволяют производить восстановление информации, поступающей от ВОИС в случае неравномерных параллельных схем укладки измерительных линий с малым числом направлений сканирования, определять пространственные распределения исследуемых физических полей или контролируемых поверхностей по информации, поступающей от измерительных сетей при малом числе ИЛ, безошибочно восстанавливать место и величину воздействия, осуществлять восстановление функций с погрешностью не более 3%, оптимизировать принятие компромиссных решений между разрешением и устранением глобальных артефактов (иС^С 2), развивать исследования, направленные на создание нового класса адаптивных контрольно-измерительных устройств, существенно расширять круг решаемых некорректных задач, решение которых на практике является сложной задачей, требующей привлечения специальных методов обработки.

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

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

Ко многим задачам, обработки информации^ может быть-применена НС типа однослойный персептрон, основным преимуществом которого является простота устройства и параллельная архитектура. Качество функционирования нейронной сети* зависит от качества ее обучения, которое в свою очередь зависит от выбора параметров обучения. Поэтому рассмотрение принципов оптимизации^ обучения представляется довольно важным- для разработки высокоэффективные нейроподобных компьютерных сетей. Поэтому был предложен» [84] метод выбора параметра скорости обучения, основанного на комбинации оптимальных параметров: оптимального параметра, выбранного перед началом; обучения и оптимальным параметром после первого цикла обучения.

Нелинейное преобразование может быть реализовано с использованием многослойного персептрона, содержащего внутренние, слои нелинейных нейронов. Как, показали.результаты исследований'[169], для реконструкции? распределения физической величины по томографическим- данным, поступающими от РВОИС с нелинейной функцией преобразования, необходимо* использовать- нейронную сеть типа двухслойный персептрон.

Выбор/ в качестве НС сети радиально-базисного^ типа связан, прежде всего, с возможностью нахождения' решения- в аналитической, форме. Так использование функций гауссовского типа в качестве радиально-базисных позволяет получить решение в виде бесконечно дифференцируемой* функции, быстро убывающей на бесконечности. Кроме того, КВР1чПЧ способна выполнять, локальную аппроксимацию нелинейного отображения, в отличие от ранее исследуемого многослойного персептрона, который осуществляет глобальную аппроксимацию. Таким образом, КВРКЫ способна выполнять нелинейные преобразования входного пространства. НВРТчПЧ обладает высокой скоростью обучения. Это объясняется отсутствием этапа обучения в общепринятом смысле: на первом этапе выделяются компактные группы кластеров, что делает начальные условия обучения оптимальными, а сам алгоритм обучения упрощенным и выигрышным во времени.

Предложен способ создания функций для обучения Е.ВР№чГ, моделирующих конкретный физический процесс. Экспериментально показано, что данная сеть позволяет осуществлять восстановление функций пространственного распределения исследуемой физической величины с погрешностью в отдельной точке не более 1%. Ни один из приведенных для сравнения традиционных методов такого результата не дает. Построенные сети обладают хорошими прогнозирующими возможностями и могут с успехом использоваться для распределенных вычислений в задачах с применением ВОИС.

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

Для этого ЯВРКК были использованы в процессе предобработки проекционных данных, который в свою очередь заключается в последовательном выполнении нейросетевой предобработки с помощью 11ВР№<Г с последующей аппроксимацией на множестве дискретизации специального вида, а именно, на объединении классов смежности.

В результате применения Е-ВРЫИ происходит восстановление проекционных данных с недостающих направлений сканирования, что предоставляет возможность применения методов аппроксимации данных на нерегулярной сетке специального вида.

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

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

Комбинированный алгоритм был рассмотрен и для случаев предельно малого (р=2,3) числа направлений укладки ВОИЛ. При этом нужно заметить, что

1) два направления укладки требуют значительного (как минимум в два раза) увеличения объёма обучающей страницы по сравнению с тремя направлениями,

2) значительно снижается прогнозирующая способность нейронной сети 11ВРЫК по сравнению со случаем укладки в трех направлениях.

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

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

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

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

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

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

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

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

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

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

Изучение принципов применения нейросетевых вычислительных методов для обработки сигналов распределенных волоконно-оптических измерительных сетей показало возможность реализации параллельных алгоритмов обработки информации при' реконструкции пространственных распределений физических полей. Полученные результаты применены* для. создания моделей нейроподобных вычислительных систем обработки сигналов распределенных измерительных сетей и могут быть использованы для создания^ информационно-измерительных систем мониторинга протяженных физических полей в реальном времени.

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

Разработке указанных новых подходов и методов и посвящена настоящая работа.

Библиография Закасовская, Елена Владимировна, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Кульчин Ю.Н. Распределенные волоконно-оптические измерительные системы. М. : Физматлит, 2001. -272 с.

2. Аш Ж. Датчики измерительных систем. Кн.1. М. : Мир, 1992. - 480 с.

3. Ristic L. Sensor technology and devices. Chapter one. Boston: Artech House. 1994.-520 p.

4. Kersey A.D. A review of recent developments in fiber optic sensor technology // Opt. Fiber Technol.- 1996.-Vol. 2, 3.- pp. 291-317 .

5. Марков Н.Г. Автоматические системы сбора и регистрации сейсмической информации. М.: Недра, 1992. - 219 с.

6. Кульчин Ю.Н., Обух В.Ф. Пространственная фильтрация излучения многомодового световода при измерении гидроакустического давления //Квантовая электроника. 1986. - Т. 13.- С. 650.

7. Кульчин Ю.Н., Беловолов М.И., Витрик О.Б., Дианов Е.М., Обух В.Ф. Исследование модуляции фазы и состояния поляризации в маломодовом волоконном световоде при аксиальных деформациях // Квантовая электроника. 1989.- Т. 16, № 11.- С. 2301-2304.

8. Кульчин Ю.Н., Витрик О.Б., Витрик Я.И., Борисенко Л.К., Петров Ю.С., Шестопалов Е.Г. Волоконно-оптический интерферометрический метод для исследования деформаций строительных конструкций // Изв. ВУЗов Сер. Строительство. 2001.- Т. 6.- С. 113-117.

9. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С. Квазираспределенный волоконно-оптический датчик // Измерительная техника,- 1994.- Т. 1.-С. 16-17.

10. Кульчин Ю.Н., Быковский Ю.А., Витрик О.Б. Статистические характеристики когерентного излучения в многомодовых волоконных световодах // Известия вузов. Серия Радиофизика. 1990. - Т. 33, № 11. -С. 1310-1311.

11. Кульчин Ю.Н., Быковский Ю.А., Витрик О.Б., Ларкин А.И. Голографическая согласованная фильтрация сигналов в интерференционных датчиках на многомодовых волоконных световодах// Квантовая электроника.- 1990.- Т. 17, №1.- С. 95-98.

12. Кульчин Ю.Н., Быковский Ю.А., Витрик О.Б. Запись голограмм Френеля излучением, прошедшим через многомодовые волоконные световоды // Оптика и спектроскопия. 1990. - Т. 68, № 5. - С. 11601169.

13. Быковский Ю.А., Кульчин Ю.Н., Обух В.Ф., Смирнов В.Л. Коррелированная перестройка спеклов в интерферометре на многомодовом световоде // Квантовая электроника. 1990. - Т. 17, №8. - С. 1080-1083.

14. Кульчин Ю.Н., Быковский Ю.А., Витрик О.Б. Амплитудная пространственная фильтрация в обработке сигналов одноволоконного многомодового интерферометра // Квантовая электроника. 1990. - Т. 17, № 10.-С. 1377-1378.

15. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С., Каменев О.Т. Метод обработки сигналов двухмодового волоконного интерферометра//Автометрия. 1995.-Т. 5.-С. 32-35.

16. Кульчин Ю.Н., Витрик О.Б., Максаев О.Г., Кириченко О.В., Каменев О.Т. Метод электронной корреляционной обработки спекл-картин для выделения полезного сигнала одноволоконных многомодовых интерферометров // ЖТФ. 1996. - Т. 66, № 12. - С. 137140 .

17. Kulchin Yu.N., Vitrik О.В., Kirichenko O.V., Kamenev O.T., Petrov Yu.S., Maksaev O.V. Method of single-fiber multimode interferometer speckle-signal processing // Optical Engineering. 1997. - Vol. 36, № 5. - pp. 1494-1499.

18. Бусурин Б.И., Носов Ю.Р. Волоконно-оптические датчики. М. : Энергоатомиздат, 1990.

19. Евтихиев H.H., Засовин Э.А., Мировицкий Д.И. Волоконная и интегральная оптика в информационных системах. М. : Изд.МИРЭА, 1987.

20. Кульчин Ю.Н., Витрик О.Б., Обух В.Ф., Петров Ю.С. Исследование интерферометра сдвига в схеме волоконно-оптического датчика давления // Измерительная техника. 1992. - Т. 10. - С. 24-26.

21. Снайдер А., Лав Дж. Теория оптических волноводов. М. : Радио и связь, 1987.-656 с.

22. Бутусов M. М., Галкин С.Л., Оробинский С.П., Пал Б.П. Волоконная оптика и приборостроение. Л. : Машиностроение, 1987.

23. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С. Многомерная обработка сигналов с использованием волоконно-оптической измерительной сети // Квантовая электроника. 1993. - Т. 20, №5.-С. 711-714.

24. Кульчин Ю.Н., Витрик О.Б., Ромашко Р.В., Петров Ю.С., Кириченко О.В., Каменев О.Т. Томографические методы для исследования векторных полей при помощи волоконно-оптических измерительных систем // Квантовая электроника. 1997 . -Т. 24, № 5. - стр. 46.

25. Мировицкий Д.И. Мультиплексированные системы волоконно -оптических датчиков // Изм. техника. 1992. - Т. 1. - С. 40-42.

26. Мировицкий Д.И. Распределенные и квазираспределенные волоконно-оптические датчики // Изм.техника. -1991.- N11.- С. 43-44.

27. Merzbacher С., A. Kersey and Е. Friebele. Fiber optic sensors in concrete structures: a review // Smart Mater. Struct. 1996. - Vol. 5, No.2. - pp. 196208.

28. Senior J.M., Moss S.E., Cusworth S.D. Multiplexing techniques for noninterferometric optical point-sensor networks: A review // Fiber and Integrated Optics. 1998.-Vol. 17, № l.-pp. 3-20.

29. Кульчин Ю.Н., Витрик О.Б., Кириченко O.B. Восстановление физических полей с использованием двухмерной волоконно-оптической измерительной сети // Измерительная техника. 1999. - Т. 3. -С. 24-30.

30. Кульчин Ю.Н., Витрик О.Б., Горбачев К.П., Аносов А.П., Кириченко О.В., Петров Ю.С., Каменев О.Т. Волоконно-оптическая измерительная сеть для регистрации параметров колебательных процессов // Измерительная техника. 1995. - Т. 3. - С. 32-35.

31. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С., Каменев О.Т. Восстановление векторных физических полей методом оптической; томографии// Квантовая электроника. 1995. - Т. 22, №' 10. - С. 10091012.

32. Кульчин Ю.Н., Витрик О.Б., Каменев OiT., Ромашко P.B. Восстановление векторных физических полей с использованием двухмерной волоконно-оптической измерительной сети // Измерительная техника. 1999. - Т. 6. - С. 21-28 .

33. Натеррер, Ф. Математические аспекты компьютерной томографии. -М. : Мир, 1990. -288 с.

34. Луис А., Наттерер Ф. Математические проблемы реконструктивной вычислительной томографии//ТИИЭР.- 1983.-Т. 71, № З.-С. 111-125.

35. Хелгасон С. Преобразование Радона. М. : Мир, 1983 . - 152 с.

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

37. Куницын В.Е., Терещенко Е.Д. Томография ионосферы. М. : Наука, 1991.

38. Левин Г. Г., Вишняков Г. Н. Оптическая томография. М. : Радио и связь, 1989.

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

40. Deans, S.R. The Radon transform and some of its applications. N.Y. : John Wiley and Sons, 1983.

41. Малеханов А.И. О волоконно-оптической томографии акустических полей//Изв. ВУЗ "Радиофизика". 1988,- Т. 31, №1. - С. 1388-1393.

42. Кульчин Ю.Н., Каменев О.Т. Томск: ТПУ. Обучающаяся нейросеть для обработки томографических данных // Кибернетика и вуз. Томск: ТПУ. 1994, № 28. - С. 3-7.

43. Витрик О.Б. Проблема чувствительной кожи и волоконно-оптические измерительные системы // Соросовский образовательный журнал. 2001.- Т. 7, № 1,- С. 108-115.

44. Пикалов В.В., Мельникова Т.С. Томография плазмы. Н. : Наука, 1995.

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

46. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С., Воробьев Ю.Д. Квазираспределенный волоконнооптический датчик // Измерительная техника. 1993, Т. 1, С. 16-17.

47. Кульчин Ю.Н., Витрик О.Б., Петров Ю.С., Кириченко О.В., Каменев О.Т. Восстановление векторных физических полей оптическимтомографическим методом // Квантовая электроника. 1995. - Т. 22, № 10.- С. 1009-1012.

48. Denisov I.V., Kulchin Y.N., Kirichenko O.V., Sedov V.A., Drozdov R.S., Denisova E.V. Model of optoelectronic measuring intelligent system // Proc. of International Conference Physics and Control. -2003. Vol.1, - pp. 172-175.

49. Кульчин Ю. H, Денисова E. В. и др. Макет оптоэлектронной нейроподобной измерительной системы // Нано- и микросистемная техника (Микросистемная техника с 1999 г. по 2004 г.). 2003. - № 10. -стр. 40 - 42.- ISSN 1813-8586.

50. Kulchin Y. N., Denisova E.V., Denisov I. V., Rybalchenko N.A. Principal of reconstruction of the single influences on fiber-optical measuring network // Pacific Science Review. 2003. - Vol. 5 (1). - pp. 3237.

51. Патлах A. JI., Семенов А. С. Светопропускание изогнутых многомодовых оптических волокон // Квантовая электроника. 1983. -Т. 7.-С. 868-870.

52. Губарени Н.М. Вычислительные методы и алгоритмы малоракурсной компьютерной томографии. 1997. - 328 с.

53. Пикалов В.В., Преображенский Н.Г. Вычислительная томография и физический эксперимент// Успехи физических наук. 1983. - Т. 141, № 3.- С. 469-498.

54. Гелъфанд И.М., Граев М.И., Виленкин Н.Я. Интегральная геометрия и связанные с ней вопросы теории представлений. 1962. - 656 с.

55. Хермен Г. Восстановление изображений по проекциям: Основы реконструктивной томографии. М. : Мир, 1983.

56. Norton S.J. . Iterative reconstruction algorithms: Convergence as a function of spatial frequency I I J. Optical Sci. Amer. 1985, Vol. 2. - pp. 613.

57. Natterer F. Fourier Reconstruction in Tomography // Numer. Math. -1985.-Vol. 47.-pp. 343-353.

58. Shepp L.A., Logan B.F. The Fourier reconstruction of a head section // IEEE Trans. Nucl. Sci. 1974. - Vol. NS-21. - pp. 21-43.

59. Ауслендер A. JI., Вешняков Г. H., Левин Г. Г. Решение интегрального уравнения Радона в оптическом процессоре // Оптика и спектроскопия. 1980. - Т. 49, № 5. - стр. 946—951.

60. Рубашов И.Б., Тимонов А.А., Пестряков А.В. О вычислительной томографии//ДАН СССР. 1981,- Т. 258, №4.-С. 846-850.

61. Shepp L.A., Logan B.F. The Fourier reconstruction of a head section //. IEEE Trans. Nucl. Sci. 1974.-Vol. 21, № 3.- pp. 21-43.

62. Штейн И.Н. О применении преобразования Радона в голографической интерферометрии // Радиотехника и радиоэлектроника. 1972. - Т. 17, № 11. - С. 2436-2437.

63. Хорн Б.К.П. Восстановление внутренней структуры объектов с помощью различных схем многолучевого просвечивания // ТИИЭР. -1976. Т. 66, № 5. - С. 27-40.

64. Ценсор Я. Методы реконструкции изображений, основанные на разложении в конечные ряды // ТИИЭР. 1983. - Т. 71, № 3. - С. 148160.

65. Gordon R., Bender R., Herman G.T. Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography II //J. Theor. Bio.-1970.-Vol. 29.-pp. 471 -481.

66. Zakasovskaya E.V., Fadeev V.V. Restoration of point influences by the fiber-optical network in view of a priori information // SPIE Proc. APCOM. -2007.-Vol. 6675.

67. Gilbert P. Iterative methods for the three-dimentional reconstruction of an object from projection // J. Theor. Biol. 1972. - Vol. 36. - C. 105-117.

68. Herman G., Lent A. Iterative reconstruction algoritms // Comput. Biol, and Med. 1976. - Vol. 6. - pp. 273-294.

69. Марчук Г.И. Методы вычислительной математики. М. : Наука, 1980.

70. Матушкин Н.Н., Южаков А.А. Измерительные системы на основе нейронных технологий // Кибернетика и вуз. 1994. - Т. 28. - С. 92-97.

71. Каллан Р. Основные концепции нейронных сетей. Издательский дом «Вильяме», 2001. -287 с.

72. Горбань А.Н., Россиев Д.А. . Нейронные сети на персональном компьютере. Н. : Сибирская издательская фирма РАН, 1996. -276 с.

73. Muller В., Reinhardt J. Neural Networks. Berlin: Springer-Veterad, Heindelberg. -1990.

74. Осовский С. Нейронные сети для обработки информации. -Финансы и статистика, 2002. 344 с.

75. Haykin, S. Neural Networks: a Comprehensive Foundation. New Jersey, Prentice Hall : s.n., 1999. - 842 p.

76. Kulchin Yu.N., Denisov I.V., Denisova E.V., Milovanov V.I. Application of algebraic methods for restitution of cumulative distribution functions of physical fields // International workshop on Optical Beam Transformation. 2001.

77. Денисов И. В., Денисова Е. В. Применение приближенных алгебраических и нейросетевых методов решения томографическойзадачи // Электронный журнал «Исследовано в России». 2002. - Т. 201.- С. 2222-2228.

78. Кульчин lOtH., Денисов И.В., Денисова Е.В. Принципы организации, матриц связей оптической нейронной сети на объемных оптических элементах // Нейрокомпьютеры: разработка, применение. 2003. -№ 7. -С. 21-27.

79. Денисов. И.В., Денисова E.B., Рыбальченко H.A. Восстановление двоичных воздействий на волоконно-оптической, сети произвольной размерности // Материалы XLVI Всероссийской* межвузовской научно-технической конференции.-2003.-Т. 1.-G. 63-65.

80. Denisov I.V., Denisova E.V., Rybalchenko N.A. Reconstruction of twice influences on the fiber-optical measuring network // Труды II Международной конференции^ «Параллельные вычисления и задачи управления». 2004. - стр. 531-534.

81. Kulchin Yu. N., Denisova E.V., Denisov I. V. Synthesis of approximate algebraic and neural-like methods for. the solution of the tomography problem // Proc. of SPIE. 2003. - Vol. 5134. - pp. 83-89.

82. Kulchin Yu. N., Denisov I. V., Denisova E. V., Piskunov E. N. Prismatic neural chip for distributed measuring networks // Optical Memory & Neural Networks. -2003. Vol. 12, № 3. - pp. 237-242.

83. Kulchin Yu. N., Denisova E. V., Denisov I. V. Application of algebraic and neural-like methods for reconstruction of distribution functions ofphysical fields // Optical Memory & Neural Networks. 2003. - Vol. 12, № 4.-pp. 283-297.

84. Denisov I.V., Denisova E.V., Rybalchenko N.A., Sedov V.A. Definition of the fiber-optical tomography problem // Proc. of International Conference Physics and Control. 2005. - Т. 1. - C. 826 - 829.

85. Гридин A.A., Денисова E.B., Денисов И.В. Обработка информации с распределенного температурного поля // Proc. of 7-th International Conference Digital Signal Processing and its Applications , Moscow, Russia. 2005.-pp. 365-368.

86. Kulchin Yu.N., Zakasovskaya E.V. Application of Radial Basis Function Neural Network for Information Processing in Fiber Optical Distributed Measuring Systems // Optical Memory & Neural Networks (Information Optics).-2008.-Vol. 17, №4.-pp. 317-327.

87. Kulchin Yu. N., Zakasovskaya E. V. Artifacts suppression in limited data problem for parallel fiber optical measuring systems // Optical Memory & Neural Networks (Information Optics). 2009. - Vol. 18, №> 3. - pp. 171180.

88. Кульчин Ю.Н., Закасовская E.B. Неравномерные схемы укладки измерительных линий в распределенных волокнно-оптических системах // Информатика и системы управления. 2009, № 3(21) . - С. 61-71.

89. Кульчин Ю.Н., Закасовская Е.В. Нейросетевое и алгебраическое моделирование параллельного 2d проецирования в волоконно-оптической томографии при ограниченном числе направлениий сканирования // Компьютерная оптика. 2009. - Т. 33, № 3. - С. 318324.

90. Кульчин Ю.Н., Закасовская Е.В. Восстановление информации волоконно-оптическими измерительными системами с использованием радиально-базисных нейронных сетей // Нейроинформатика-2009. — Сб. науч. тр. М.: МИФИ.-2009. С.289-298.

91. Кульчин Ю.Н., Закасовская Е.В. Оптимизация алгоритма обработки информации в распределенных волоконно-оптических измерительных системах // Информатика и системы управления. 2010. - № 4 (26). - С. 50 - 60. - ISSN 1814-2400.

92. Кульчин Ю.Н., Закасовская Е.В. Моделирование параллельного 2D-проецирования в волоконно-оптической томографии для малого числа направлений сканирования // Информатика и системы управления. -2010, Т. 1(23).-С. 104-114.

93. Kulchin Yu. N., Zakasovskaya E. V. Optimizing algebraic and neural methods for information processing in distributed fiber-optical measuring systems // Optical Memory & Neural Networks (Information Optics). -2010. Vol. 19, № 3. - C. 237-247.

94. Кульчин Ю.Н., Закасовская Е.В. Обработка информации комплексом нейронных сетей в распределенных волоконно-оптических измерительных системах // Компьютерная оптика. 2010. - Т. 34, № 3. -С. 363 - 369.

95. Кульчин Ю.Н., Закасовская Е.В. Нейросетевое подавление артефактов при малом числе направлений укладки измерительных линий в волоконно-оптических системах // Нейроинформатика-2010, Сб. науч. трудов МИФИ.-2010, часть 1.-С. 110-120.

96. Kulchin Yu.N., Zakasovskaya E.V. Neural-like specialized construction for fiber optical distributed measuring systems in limited-data conditions. Pacific Science Review. 2010, Vol. 12, № 1, pp. 80-86.

97. Kulchin Yu. N., Zakasovskaya E. V. Complexes of neural networks for information processing in distributed fiber-optical measuring systems. The First Russia and Pacific Conference on Computer Technology and Applications (RPC 2010). 2010, pp. 125-129.

98. Кульчин Ю.Н., Закасовская Е.В. Нейросетевые конструкции для распределенных волоконно-оптических измерительных систем //

99. МИФИ, Нейроинформатика-2011. Сб. науч. тр. М.: МИФИ. - 2011. -С. 171-180.

100. Кульчин Ю.Н., Закасовская Е.В. Комплекс нейронных сетей для синтеза синограмм в распределенных волоконно-оптических системах // Нейрокомпьютеры: разработка, применение. 2011. - № 7. - С. 55-62.

101. Kulchin Yu.N., Vitrik О.В. , Kirichenko O.V., Petrov Yu.S. Distribute fiber-optic sensor for seismoacoustic investigation // Proc. 3-rd International Russian Fiber Optic Conference. 1993.-Vol. 2.-pp. 291-294.

102. Kulchin Yu.N., Vitrik O.B., Kirichenko O.V., Petrov Yu.S., Kamenev O.T. Distributed fiber-optic acoustic sensor // Proc. of Distributed and multiplexed fiber optic sensor IV Conf. 1994. - Vol. 2294. - pp. 129-132.

103. Shannon С. Е. Communication in the presence of noise // Proc. Institute of Radio Engineers.- 1949.-T. 37, № 1,-pp. 10-21.

104. Наймарк M.A. Теория представлений групп.-M. : Наука, 1976.

105. Petersen D.P., Middleton D. Sampling and reconstruction of wave-number-limited functions in N-dimensional Euclidean space // Inf. Control. -1962. Vol. 5. -pp. 279-323.

106. Faridani A. An application of a multidimensional sampling theorem to computed tomography // Integral Geometry and Tomography, Contemporary Mathematics. AMS, Providence, RI. 1990. - Vol. 113.

107. Kluv'anek, I. Sampling theorem in abstract harmonic analysis // Mat. Casopis Sloven. Akad. Vied. 1965. - Vol. 15. - pp. 43-48.

108. Faridani A. Introduction to the mathematics of computed tomography // Inverse Problems and Applications. -2003. -pp. 1-46.

109. Behmard H. and Faridani A. Sampling of bandlimited functions on unions of shifted lattices // Fourier Anal. Appl. 2002. - Vol. 8, № 1. - pp. 43-58.

110. McCalloch W.A., Pitts W. A logical calculus of ideas immanents in nervous activity//Bull. Math. Biophys. 1943. -T. 5.-стр. 115-133.

111. Розенблатт Ф. Принципы нейродинамики. Перцептрон и теория механизмов мозга. М.: Мир, 1965. - 480 с.

112. Крюков В.И. Что такое нейрокомпьютеры? Нейрокомпьютер как основа мыслящих ЭВМ. - стр. 84-86.

113. Natterer F. The Mathematics of Computerized Tomography. B.G. Teubner and John Wiley and Sons, Stuttgart. - 1986.

114. Кульчин Ю. Н., Денисов И. В., Каменев О. Т. Оптоэлектронная нейроподобная система обработки выходных данных волоконно-оптической измерительной сети // ПЖТФ. 1999 . - Т. 25, № 6. - С. 6570.

115. Кульчин Ю.Н., Витрик О.Б., Кириченко О.В., Петров Ю.С., Каменев О.Т. Самообучающаяся нейронная сеть для обработки томографических данных // Сборник трудов "Кибернетика и высшая школа".- 1994 .-Т. 28.-С. 3-7.

116. Kulchin Yu., Kamenev О. Self-training neural network model for real time tomography data processing // Laser Biology. 1995. - Vol. 4, № 2. -pp. 625-629.

117. Kamenev O.T. Training two-layer neural network model for tomography data processing // Proc. of Int. Conf. OCEAN'95. San-Diego. - 1995. - pp. 2086-2087.

118. Кульчин Ю.Н., Денисов И.В., Каменев О.Т. Оптоэлектронная нейроподобная система обработки выходных данных волоконно-оптической измерительной сети // Письма в ЖТФ. 1999. - Т. 25, № 6. -С. 65-70.

119. Галушкин А. И. Нейрокомпьютеры. Кн.З. Учеб. пособие для вузов. - 2000. - 528 с.

120. Калан Р. Основные концепции нейронных сетей. М.: Издательский дом "Вильяме". 2001.

121. Muller В., Reinhardt J. Neural Networks. Berlin : Heidelberg, 1990.

122. Wasserman, P.D. Advanced Methods in Neural Computing. New York : Van Nostrand Reinhold, 1993.

123. Yu Hen Hu, Jeng-Neng Hwang. Handbook of Neural Network signal processing. London, New York, Washington, D.C. : CRC Press, Boca Raton, 2002.

124. Rosenblatt F. The Perceptron: A Probabilistic Model for Information Storage and Organization in the Brain // Cornell Aeronautical Laboratory, Psychological Review. 1958. - Vol. 65, № 6. - pp. 386 - 408.

125. Минский M., Пайперт С. Персептроны. — Москва : Мир, 1971.

126. Kulchin Yu., Denisov I., Obuh V., Kamenev O., Vitrik O., Petrov Yu., Romashko R. Computer neural networks for processing of optical tomography information // Pacific Science Review. 1999. - Vol. 1. - pp. 1-4.

127. Denisov I., Kulchin Yu., Obuh V., Kamenev O. Principales of organization of neural-like system on the basis of a matrix of photoelectric cells Optoelectronic Information // Systems and Processing. 2001. - Vol. 4513.-pp. 52-57.

128. Kulchin Yu., Kamenev O., Denisov I. Neural processing system for optical information measuring systems Distributed fiber optical sensors and measuring networks//Proc. of SPIE.-2001.- Vol. 4357.-pp. 109-117.

129. Беклемишев Д.В. Дополнительные главы линейной алгебры.— М. : Наука, 1983.

130. Уоссерман Ф. Нейрокомпьютерная техника. Теория и практика. — М.: Мир, 1992.-240 с.

131. Blair L. Т., S. A. Cassidy. Wavelength division multiplexed sensor network using Bragg fibre reflection gratings // Electron. Lett. 1992. - Vol. 28, № 18.-pp. 1734-1735.

132. Haykin S. Neural networks, a comprehensive foundation. N.Y. : Macmillian College Publishing Company, 1994.

133. Kulchin Yu. N., Panov A. V. Neural network for reconstruction of signal from distributed measuring system of optical amplitude sensors // Pacific Science Review.-2001.-Vol. 3.-pp. 1-4.

134. Denisov I. V., Kamenev О. Т., Kim A. Yu., Kulchin Yu. N., Panov A. V. Neural data processing method for fiber-optic distributed measuring systems // Optical Memory & Neural Networks. 2003. - Vol. 12, № 3. - pp. 165172.

135. Hutchins D. A., Mottram J. T. ,Hines E. L. , Corcoran P., Anthony D. M. A neural network approach to ultrasonic tomography // Proc. IEEE. 1992. -Vol. l.-pp. 365-368.

136. Powell M.J.D. Radial Basis Function approximation for multivariable interpolations: A review // RMCS, Shrivenham, England. 1985. - pp. 143167.

137. Micchelli C.A. Interpolator of scattered data: Distance matrices and cjnditionally positive definite functions // Constructive Aproximation. 1986. -Vol. 2.-pp. 11-22.

138. Hardy R. Multiquadric equations of topography andother irregular surfaces //J. Geophys. Res. 1971.-Vol. 76.-pp. 1905-1915.

139. Kansa E.J. Multiquadrics A scattered data approximation scheme with applications to fluid dynamics I. Surface approximations and partial derivative estimates // Computers and Mathematics with Applications. -1990.-T. 19.-C. 127-145.

140. Wendland H. Scattered Data Approximation. Cambridge : Cambridge University Press, 2005.

141. Schoenberg I. Metric spaces and completely monotone functions // Ann. of Math. 1938. -Vol. 39.-pp. 811-841.

142. Buhmann M. D. Radial basis functions: theory and implementations // Cambridge Monographs on Applied and Computational Mathematics. -2003.-Vol. 12.

143. Wendland H. Piecewise polynomial, positive definite and compactly supported radial basis functions of minimal degree // Advances in Computational Mathematics. 1995. - Vol. 4. - pp. 389-396.

144. Ильин В.А., Поздняк Э.Г. Линейная алгебра. -М. : Наука, 1974.