Способ оперативной идентификации морских целей по их информационным полям на базе нейро-нечетких моделей
METHOD FOR OPERATIONAL IDENTIFICATION OF MARINE TARGETS BY THEIR INFORMATION FIELDS BASED ON NEURO-FUZZY MODELS

Изобретение относится к гидроакустике и может быть использовано для реализации операций нейросетевого распознавания классов целей (надводный или подводный объект), обнаруженных по признакам амплитудно-фазовой модуляции низкочастотных сигналов накачки морской среды излучениями и полями объектов.
Принцип работы параметрических антенн как с высокочастотной (десятки-сотни кГц), так и с низкочастотной (десятки-сотни Гц) накачкой морской среды, основан на использовании естественных нелинейных свойств морской среды.
Известны способы и реализующие их параметрические антенны, использующие высокочастотную накачку морской среды (см. Новиков Б.К., Тимошенко В.И. Параметрические антенны в системах гидролокации. - Л.: Судостроение. - 1990. - С. 17-40, 203-225; Мироненко М.В., Малашенко А.Е., Карачун Л.Э., Василенко A.M. Низкочастотный просветный метод дальней гидролокации гидрофизических полей морской среды: монография. - Владивосток: СКБ САМИ ДВО РАН, 2006. - 173 с.; Пятакович В.А., Василенко A.M., Филиппова А.В. Технология создания автоматизированной системы дальнего приема и нейросетевой классификации гидрофизических полей морских акватории. - Известия Тульского государственного университета. Технические науки. Вып. 7. Тула: Изд-во ТулГУ, 2017. - С. 247-258; Пятакович В.А., Василенко A.M., Рычкова В.Ф. Интеллектуальная система нейросетевой классификации морских целей. - Морские интеллектуальные технологии. - СПб.: 2018. №2 (40). том 2. - С. 115-120; Пятакович В.А., Рычкова В.Ф., Филиппов Е.Г. Нейронные сети как вариант вычислительной структуры системы классификации морских целей. - Морские интеллектуальные технологии. - СПб.: 2020. №1 (47) Том 2. - С. 163-174; Пятакович В.А., Суров А.Б., Рычкова В.Ф. Расчет эффективности классификации целей интеллектуальной системой ВМФ, использующей комплекс вычислительных операций нейронных сетей. - Морские интеллектуальные технологии. - СПб.: 2020. №1 (47) Том 2. - С. 175-185; Пятакович В.А., Николаев А.В., Костиков Е.А. Автоматизация обработки информации в интеллектуальной системе морского мониторинга. - Проблемы машиностроения и автоматизации. - М.: 2020. №4. - С. 72-79).
Известны способы приема упругой волны в морской среде, повышающие величину параметра нелинейности среды в зоне приема сигналов. Для чего в зону параметрического приема сигналов, кроме упругой волны, вводят дополнительный электромагнитный сигнал, подвергнутый частотно-временной модуляции, или модулированный гидродинамический поток жидкости. Модуляцию гидродинамического потока жидкости осуществляют путем изменения плотности за счет управляемого насыщения газом, или изменения температуры, или изменения ее химического состава, (пат. №2158029 РФ. Способ приема упругой волны в морской среде (варианты); Опубл. 20.10.2000. Бюл. №29.; пат. №2452041 РФ. Способ параметрического приема волн различной физической природы в морской среде. Опубл. 27.05.2012. Бюл. №15.). Недостатками способов являются низкая чувствительность и, как следствие, ограниченная дальность (единицы километров) параметрического приема информационных волн различной физической природы в инфразвуковом и дробном (единицы-доли герц) диапазонах частот.
Наиболее близким по технической сущности к заявляемому изобретению является «Способ обнаружения и нейросетевого распознавания признаков полей различной физической природы, генерируемых морскими целями» (пат. №2682088 С1 (РФ, МПК G01S 15/02 (2006.01); Опубл. 14.03.2019. Бюл. №8), заключающийся в том, что сначала формируют в морской среде зону нелинейного взаимодействия и параметрического преобразования волн накачки с сигналами объектов, для чего излучатель и приемную антенну размещают на противоположных границах контролируемого участка морской среды, затем волны накачки, промодулированные с сигналами объекта, принимают и усиливают в полосе параметрического преобразования, переносят их частотно-временной масштаб в высокочастотную область, проводят узкополосный спектральный анализ, выделяют параметрические составляющие суммарной или разностной частоты, по которым с учетом временного и параметрического преобразования волн восстанавливают характеристики сигналов объекта, которые подают в тракт нейросетевого распознавания и классификации, состоящий из блока распознавания класса цели по амплитудно-частотным характеристикам, реализующего вычислительные операции искусственной нейронной сети и охваченного обратной связью с блоком обучения, в память которого записаны данные математически обработанных образов спектрограмм морских целей, причем на первый вход блока распознавания класса цели по амплитудно-частотным характеристикам поступают данные с выхода спектроанализатора тракта приема, обработки и регистрации сигналов, а на второй его вход поступают данные с блока обучения тракта нейросетевого распознавания и классификации, после чего искусственную нейронную сеть настраивают по классификационным признакам полей, генерируемых морскими целями, запускают вычислительные операции, по результатам которых корректируют ее весовые коэффициенты и формируют вывод о степени принадлежности исследуемой области спектра объекту классификации (надводный или подводный объект).
Известно, что результатом параметрического преобразования взаимодействующих волн является их взаимная амплитудно-фазовая модуляция. Малое отличие частот (в пределах одного порядка) просветных волн и волн, генерируемых объектом, обеспечивает наиболее интенсивное их взаимодействие. Амплитуда взаимодействующих волн и индекс фазовой модуляции могут быть представлены в следующем виде
n
c
n
c
0
0
где γ - коэффициент нелинейности морской среды; ω, ω - частота волны накачки и полезного сигнала, соответственно; P, P - затухание волны накачки и полезного сигнала, соответственно; V - объем среды нелинейного взаимодействия и параметрического преобразования волн; R - расстояние от точки излучения до точки расположения объекта; ρ - плотность, c - скорость звука в морской среде.
Сформированные в результате преобразования просветных волн параметрические составляющие суммарной и разностной частоты при обработке широкополосных сигналов выделяются, как признаки амплитудно-фазовой модуляции, что обосновано математическими зависимостями и подтверждено результатами морских экспериментов (см. Мироненко М.В., Малашенко А.Е., Карачун Л.Э., Василенко A.M. Низкочастотный просветный метод дальней гидролокации гидрофизических полей морской среды: монография. - Владивосток: СКБ САМИ ДВО РАН, 2006. - 173 с.; Пятакович В.А., Василенко A.M., Мироненко М.В.; Технологии нелинейной просветной гидроакустики и нейро-нечетких операций в задачах распознавания морских объектов: - монография. - Владивосток: Дальневост. федерал, ун-т, 2016 190 с. ISBN 978-5-7444-3790-9.; Мироненко М.В., Пятакович В.А., Василенко A.M. Результаты экспериментальных исследований способа определения профиля морского объекта и реализующей его системы. - Мониторинг. Наука и технологии. 2017. - №2 (31) - С. 64-69.).
Спектр взаимодействующих волн состоит из бесконечного числа боковых составляющих, частоту и амплитуду которых можно найти из известного выражения
c
n
n
p
где  P - результирующее и мгновенное значения давления модулированной волны, соответственно; 2ω - удвоенная частота модулированной волны; Ω - волна, генерируемая объектом; t - время; J - функции Бесселя n-го порядка; A - амплитуда модулированной волны; m - коэффициент модуляции. Как видно из выражения, значения частот боковых составляющих отличаются от удвоенной центральной частоты 2ωо (равной сумме частот взаимодействующих волн) на величину ±n⋅Ω, где n - любое целое число. Амплитуды боковых составляющих для соответствующих частот (2ω±n Ω) определяются величиной множителя
p
При малых значениях коэффициента модуляции m спектр взаимодействующих волн приближенно состоит из удвоенной центральной частоты 2ω и ее боковых частот 2ω+Ω и 2ω-Ω.
Недостатком способа-прототипа является отсутствие операций, которые обеспечивают регулирование параметров формирования и редукции выборок блоком обучения, комплексное сокращение размерности данных при автоподстройке базы правил, исходя из выборки значений параметров объекта классификации, формирование сигнала с номером продукционного правила и типа функции принадлежности типу цели для оптимизации вычислительных процессов, выполняемых в тракте нейросетевого распознавания и классификации, обеспечивающем конечное классификационное решение по обнаруженным морским целям (надводный или подводный объект), что ограничивает функциональные возможности системы-прототипа.
Задача, на решение которой направлено заявляемое изобретение, заключается в расширении функциональных возможностей способа-прототипа, т.е. его реализации как способа оперативной идентификации морских целей по их информационным полям на базе нейро-нечетких моделей. Конечное классификационное решение по обнаруженным морским целям (надводный или подводный объект) принимается обученной нейронной сетью по степени принадлежности исследуемой области спектра объекту классификации.
Техническим результатом предлагаемого изобретения является автоматизация процесса распознавания классов морских целей (надводный или подводный объект), обнаруженных по признакам амплитудно-фазовой модуляции низкочастотных сигналов накачки морской среды излучениями и полями объектов, комплексное сокращение размерности данных при автоподстройке базы правил за счет формирования и редукции выборки эталонных образцов математических моделей морских целей, осуществляемой при помощи тракта регулирования параметров формирования и редукции выборок, необходимых для реализации конечного процесса классификации в тракте нейросетевого распознавания и классификации, что обеспечивает повышение вероятности достоверной оперативной классификации морской цели на 7-9% большую, чем при использовании прототипа.
Указанный технический результат достигается путем формирования и редукции выборки эталонных образцов математических моделей морских целей трактом регулирования параметров формирования и редукции выборок, самостоятельно производящим автоподстройку своей базы правил и ее нейро-нечеткую коррекцию, с применением вычислительных операций модифицированной комбинированной распознающей сети, состоящей из сетей Кохонена - Гросберга и адаптивной нейро-нечеткой сети ANFIS, использующей гибридный алгоритм обучения с генетической настройкой параметров, для оперативно обновляемой библиотеки спектрограмм морских целей блока обучения тракта нейросетевого распознавания и классификации, обеспечивающего конечное классификационное решение по обнаруженным морским целям.
Для решения поставленной задачи разработан способ оперативной идентификации морских целей по их информационным полям на базе нейро-нечетких моделей, заключающийся в том, что сначала формируют в морской среде зону нелинейного взаимодействия и параметрического преобразования волн накачки с сигналами объекта, для чего излучатель и приемную антенну размещают на противоположных границах контролируемого участка морской среды, затем волны накачки, промоделированные сигналами объекта, принимают и усиливают в полосе параметрического преобразования, переносят их частотно-временной масштаб в высокочастотную область, проводят узкополосный спектральный анализ, выделяют параметрические составляющие суммарной или разностной частоты, по которым с учетом временного и параметрического преобразования волн восстанавливают характеристики сигналов объекта, которые подают в тракт нейросетевого распознавания и классификации, состоящий из блока распознавания класса цели по амплитудно-частотным характеристикам, реализующего вычислительные операции искусственной нейронной сети и охваченного обратной связью с блоком обучения, в память которого записаны данные математически обработанных образов спектрограмм морских целей, причем на первый вход блока распознавания класса цели по амплитудно-частотным характеристикам поступают данные с выхода спектроанализатора тракта приема, обработки и регистрации сигналов, а на второй его вход поступают данные с блока обучения тракта нейросетевого распознавания и классификации.
Принципиальным отличим от прототипа является то, что данные подают в дополнительно введенный тракт регулирования параметров формирования и редукции выборок на вход, охваченного обратной связью с блоком фаззификации, блока нечетких правил, где формируют сигнал номера нового продукционного правила и нового типа функции принадлежности типу цели и далее подают на вход логического устройства формирования и редукции выборок, функцию которого выполняет модифицированная комбинированная распознающая сеть, состоящей из сетей Кохонена - Гросберга и адаптивной нейро-нечеткой сети ANFIS, использующей гибридный алгоритм обучения с генетической настройкой параметров, в котором осуществляют компенсацию функций принадлежности нового вида правил, определяют номер правила, необходимого для замены в основной базе правил, а также новый тип функции принадлежности для данного правила, и далее на блок фаззификации, выход которого соединен с входом блока дефаззификации и автоматического регулирования базы правил, в котором производят автоподстройку своей базы правил, исходя из выборки образцов математических моделей морских целей, формирование и редукцию выборки эталонных образцов математических моделей морских целей и коррекцию данных оперативно обновляемой библиотеки математически обработанных образов спектрограмм морских целей для блока обучения тракта нейросетевого распознавания и классификации, после чего искусственную нейронную сеть настраивают и формируют вывод о степени принадлежности исследуемой области спектра объекту классификации (надводный или подводный объект).
Как известно, извлечение полезной информации из гидроакустических сигналов определяет основы алгоритмизации обработки данных в способе обнаружения и нейросетевого распознавания признаков полей различной физической природы, генерируемых морскими целями. Для формирования вектора признаков, являющегося входным информационным массивом распознающей сети, используется метод масок. Процесс формирования информационных массивов необходим для решения двух задач, первая из которых представляет собой процесс формирования эталонных образцов, необходимых для реализации процесса обучения распознающей сети, и вторая для распознавания целей (см. Пятакович В.А., Богданов В.И., Назаренко П.К. Принцип автоматического распознавания образа цели: материалы Международной конференции «Математическое моделирование физических, экономических, технических, социальных систем и процессов». - Ульяновск: УГУ, 2003. - С. 31, 32; Пятакович В.А., Василенко A.M., Хотинский О.В. Распознавание и классификация источников формирования полей различной физической природы в морской среде: монография. - Владивосток: Мор. гос. ун-т, 2017. - 255 с.; Пятакович В.А., Василенко A.M., Хотинский О.В. Нейросетевые технологии в интеллектуальных системах обнаружения и оперативной идентификации морских целей: монография. - Владивосток: Мор. гос. ун-т, 2018. - 263 с.; Пятакович В.А., Василенко A.M., Пашкеев С.В. Автоматизированная система мониторинга морской среды для решения задач обнаружения технического объекта. - Двойные технологии. - М: 2018. №4 (85). - С. 85-88; Пятакович В.А., Рычкова В.Ф. Параметрическая оптимизация нейросетевой системы классификации морских целей по критерию надежности. - Морские интеллектуальные технологии. - СПб.: 2018. №4 (42) Том 5. - С. 153-162; Пятакович В.А. Опытная модель разработанного образца распознающего модуля системы классификации морских целей на базе нейросетевых технологий. Проблемы и методы разработки и эксплуатации вооружения и военной техники ВМФ: Сб. статей. - Владивосток: ТОВВМУ, 2018. - Вып. 98. - С. 176-181; Пятакович В.А. Решение задачи классификации морской цели с использованием нейронных нечетких сетей на основе модели вывода Мамдани-Заде. Проблемы и методы разработки и эксплуатации вооружения и военной техники ВМФ: Сб. статей. - Владивосток: ТОВВМУ, 2019. - Вып. 101. - С. 275-284; Пятакович В.А., Филиппов Е.Г. Метод формирования и редукции выборок для решения задач классификации морских целей нейроклассификатором. Проблемы и методы разработки и эксплуатации вооружения и военной техники ВМФ: Сб. статей. - Владивосток: ТОВВМУ, 2019. - Вып. 101. - С. 312-320; Пятакович В.А. Обучение глубоких нейронных сетей системы оперативной идентификации морских целей ВМФ РФ. Проблемы и методы разработки и эксплуатации вооружения и военной техники ВМФ: Сб. статей. - Владивосток: ТОВВМУ, 2020. - Вып. 104. - С. 177-186; Пятакович В.А., Рычкова В.Ф., Филиппов Е.Г. Нейронные сети как вариант вычислительной структуры системы классификации морских целей. - Морские интеллектуальные технологии. - СПб.: 2020. №1 (47) Том 2. - С. 163-174; Пятакович В.А., Суров А.Б., Рычкова В.Ф. Расчет эффективности классификации целей интеллектуальной системой ВМФ, использующей комплекс вычислительных операций нейронных сетей. - Морские интеллектуальные технологии. - СПб.: 2020. №1 (47) Том 2. - С. 175-185; Пятакович В.А., Николаев А.В., Костиков Е.А. Автоматизация обработки информации в интеллектуальной системе морского мониторинга. - Проблемы машиностроения и автоматизации. - М: 2020. №4. - С. 72-79).
Идея состоит в том, что для каждой маски ищется максимальное амплитудное значение, которое и является ортом вектора классификационных признаков. Для автоматизации процесса поиска экстремума в зоне одной маски использовалась сеть поиска максимума MAXNET. Итерации сети завершаются после того, как выходные нейроны сети перестают меняться. Тип элементов входных сигналов - целые или действительные числа, тип элементов выходных сигналов - действительные числа. Размерности входных и выходных сигналов совпадают. Тип активационной функции - линейная с насыщением (используется линейный участок). Число синапсов в сети равно N (N-1). Формирование синаптических весов происходит согласно формуле
ij
j
j
где W - i-й синаптический вес j-го нейрона; N - число элементов входного сигнала (количество нейронов в сети). Функционирование сети задается выражением y(0)=x,  i=1, …, N; j=1, …, N.
j
i
где x - элемент (орт) входного сигнала сети; y - выход j-го нейрона.
Нормализация входного вектора признаков, полученного после анализа масок сетью MAXNET, производится согласно выражения
Границы диапазона значений  известны и определяются моделью входного гидроакустического сигнала.
*
*
Обучение распознающей сети производится на основе алгоритма обратного распространения ошибки, реализующего градиентный метод оптимизации функционала вида: , где Т - вектор синаптических весов сети; (X, Y) - обучающие пары;  - норма вектора (см. Пятакович В.А., Василенко A.M., Хотинский О.В. Распознавание и классификация источников формирования полей различной физической природы в морской среде: монография. - Владивосток: Морской гос. ун-т им. Г.И. Невельского, 2017. - 255 с.; Пятакович В.А., Василенко A.M. Перспективы и ограничения использования геометрических методов распознавания акустических образов морских объектов применительно к задаче управления нейросетевой экспертной системой. - Фундаментальные исследования. - М: 2017. - №7. - С. 65-70; Пятакович В.А., Василенко A.M. Мироненко М.В. Обучение нейронной сети как этап разработки экспертной системы для классификации источников физических полей при мониторинге акваторий. - Вестник Инженерной школы Дальневосточного федерального университета. - Владивосток: Дальневост. федерал, ун-т, 2017. №3 (32). - С. 138-149. DOI.org/10.5281/zenodo.897021; Пятакович В.А. Система классификации морских целей на базе нейросетевых технологий. - Морские интеллектуальные технологии. - СПб.: 2018. №4 (42) Том 5. - С. 169-176).
Для автокоррекции и регулирования алгоритма обратного распространения ошибки при обучении распознающей сети используется тракт регулирования параметров формирования и редукции выборок, сочетающий в себе адаптивные подходы самообучения и опыт эксперта, что при наличии неопределенных параметрических возмущений, позволяет оперативно осуществлять коррекцию значений, соответствующих новым условиям классификации морской цели.
Изобретение поясняется чертежами, где на фиг. 1 показана функциональная схема системы оперативной идентификации морских целей по их информационным полям на базе нейро-нечетких моделей, содержащая следующие элементы:
1. Излучающий преобразователь (подводный звуковой маяк марки ПЗМ-400 излучающий сигналы на частоте около 400 Гц).
2. Приемный преобразователь.
3. Морская среда.
4. Рабочая зона нелинейного взаимодействия и параметрического преобразования волн накачки и информационных волн.
5. Объекты (морские цели, генерирующие акустические, электромагнитные и гидродинамические излучения).
6. Тракт излучения сигналов накачки.
6.1. Генератор сигналов накачки стабилизированной частоты.
6.2. Усилитель мощности.
6.3. Блок согласования.
7. Тракт приема, обработки и регистрации информационных сигналов.
7.1. Широкополосный усилитель.
7.2. Преобразователь частотно-временного масштаба.
7.3. Спектроанализатор.
7.4. Регистратор.
8. Тракт нейросетевого распознавания и классификации.
8.1. Блок распознавания класса цели по амплитудно-частотным характеристикам.
8.2. Блок обучения.
9. Тракт регулирования параметров формирования и редукции выборок.
9.1. Блок нечетких правил и функций.
9.2. Блок фаззификации.
9.3. Логическое устройство формирования и редукции выборок.
9.4. Блок дефаззификации и автоматического регулирования базы правил.
Общая структура модифицированной комбинированной распознающей сети, состоящей из сетей Кохонена и Гросберга представлена на фиг. 2.
Модификация сети состоит в добавлении к сети Кохонена и Гросберга сети MAXNET, что весьма важно для решаемой задачи.
ij
(1)
На каждый нейрон первого слоя через синапсы с весами {T}, i=1, 2, 3; j=1, 2, 3, подаются все компоненты входного вектора  Число нейронов во втором (скрытом) слое определяется взаимным расположением и формой разделяемых множеств.
ij
(2)
На каждый нейрон второго слоя через синапсы с весами {T}, i=1, 2, 3; j=1, 2, 3, подаются выходные сигналы первого слоя. Число нейронов третьего (выходного) слоя определяется числом рассматриваемых классов, подлежащих распознаванию.
ij
2n
n
n
i
(3)
(i)
(i)
(i)
На каждый нейрон третьего слоя через синапсы с весами {T}, i=1, 2, 3; j=1, 2, 3, подаются выходные сигналы второго слоя. Значения выходных сигналов третьего слоя образуют вектор  решений. Нейроны, составляющие сеть, одинаковы и имеют функцию активации известного типа  где x, y и I - значения r-го входного сигнала, выходного сигнала и внешнего смещения n-го нейрона i-го слоя; N - число нейронов в i-м слое; i=1, 2, 3.
Предварительная обработка входных векторов производится нормализацией входного вектора признаков, полученного после анализа масок сетью MAXNET или после получения статистических оценок согласно выражения  Границы диапазона значений  известны и определяются моделью входного гидроакустического сигнала.
Предварительная обработка данных, один из самых трудозатратых шагов, необходима, чтобы последующие алгоритмы обучения смогли извлечь из выборки максимум знаний.
Процедура архитектурной самонастройки и настройки весовых коэффициентов (синапсов) сети Кохонена включает в себя следующие этапы.
1. Вводится произвольное число нейронов L=L0 со случайным образом нормализованными векторами синапсов, равномерно распределенными на интервале [-1, 1].
2. На слой подается один из обучающих векторов входных сигналов (очевидно, что предварительно обработанный), определяются потенциалы на выходах всех нейронов и номер L*-нейрона «победителя».
3. Определяется угол β* между обучающим вектором признаков и вектором синапсов (весовых коэффициентов) нейрона "победителя".
4. Если выполнено условие β*<β, то осуществляется настройка синапсов нейрона "победителя" путем усреднения по всем шагам обучения, на которых он оказался нейроном "победителем" и последующей нормализации.
Если β*>β, то в слой директивным порядком вводится L+1 нейрон, синапсом которого (весам связей) присваиваются значения компонентов соответствующего обучающего вектора.
5. Выбирается очередной входной вектор обучающей выборки и повторяются процедуры 2, 3, 4.
Обучение слоя Гросберга является традиционным обучением с учителем и может производиться как одновременно с архитектурной самоорганизацией и настройкой слоя Кохонена по каждому входному вектору выборки, так и автоматно (после обучения слоя Кохонена). Во всех случаях правило обучения можно представить следующим алгоритмом:  l=1, 2. где  - l-й компонент желаемого выходного вектора классификатора;  - выходной сигнал j-го нейрона слоя Кохонена, при обучении S-му входному вектору признаков.
Функционирование рассматриваемого распознающего устройства. Предварительно обработанный вектор входных признаков наблюдаемого объекта {xi} подается на вход сети и распределяется через весовые коэффициенты связей (синапсы)  на нейроны слоя Кохонена. Входные потенциалы нейронов  i=1, …, L.
После этого слой Кохонена начинает функционировать как конкурирующая сеть с латеральными связями. В результате на выходе слоя формируется вектор с одним единичным компонентом, соответствующим нейрону «победителю» и с нулевыми остальными коэффициентами. Через синапсы  вектор входа поступает на нейроны выходного слоя Гросберга, которые функционируют согласно следующего алгоритма:  l=1, 2, где 10 ( ) - функция единичного скачка.
На фиг. 3 представлена нейросетевая интерпретация алгоритма многомерной классификации.
На фиг. 4 приведена таблица интерпретации трехмерного выходного вектора распознавания гидроакустических сигналов по амплитудно-частотной характеристике.
На фиг. 5 представлен метод масок, используемый для распознавания по амплитудно-частотной характеристике.
1
2
j
k
В каждой маске по реальной характеристике определяется максимальное амплитудное значение сигнала A, А, …, A …, A. Выбор величины Δ, а, следовательно, и числа масок определяется возможностями распознающей сети (реально 10÷100). Увеличение числа масок приводит к повышению достоверности входной информации и к повышению сложности (увеличению числа нейронов входного слоя) распознающего устройства, то есть имеет место классический конфликт между качеством и сложностью. Возможен вариант исследования шумового портрета по частям, то есть низкочастотной, среднечастотной и высокочастотной составляющих отдельно.
На фиг. 6 представлена сравнительная характеристика методов классификации морской цели по времени обучения и классификации.
Наиболее скоростным методом по качеству классификации является модифицированный метод потенциальных функций на основе НС (фиг. 6. Табл 2). НС обладают более сильными аппроксимационными способностями, чем другие методы классификации целей, а алгоритм многомерной классификации, в отличие от методов потенциальных функций, учитывает значимость признаков посредством учета частных значимостей двупризнаковых классификаций.
На фиг. 7 и фиг. 8 представлены результаты вычислительного эксперимента по определению коэффициента распознавания (классификации), определяемого как отношение числа распознанных объектов к общему числу испытаний в процентах, для надводных и подводных объектов в условиях зашумления сигнала в диапазоне от - 10 до 20 дБ. Как видно из рисунков, распознавание и классификация морских целей с помощью вычислительных операций модифицированной комбинированной распознающей сети, состоящей из сетей Кохонена (конкурирующая сеть) и Гросберга и адаптивной нейро-нечеткой сети ANFTS, использующей гибридный алгоритм обучения позволяет повысить вероятность достоверной оперативной классификации как надводных, так и подводных целей на 7-9%.
Способ оперативной идентификации морских целей по их информационным полям на базе нейро-нечетких моделей
Излучающий преобразователь 1 и приемный преобразователь 2 размещают в морской среде 3 с учетом закономерностей многолучевого распространения волн в протяженном гидроакустическом канале, что обеспечивает формирование и эффективное использование пространственно-развитой рабочей зоны 4 нелинейного взаимодействия и параметрического преобразования просветных волн и волн различной физической природы, генерируемых объектами 5 (см. Свидетельство о государственной регистрации программы для ЭВМ «Расчет лучевой картины» №2016616822. 21.06.2016; Свидетельство о государственной регистрации программы для ЭВМ «Программа имитационного моделирования процесса распространения гидроакустических сигналов» №2017664296. 20.12.2017; Свидетельство о государственной регистрации программы для ЭВМ «Программно-вычислительный комплекс имитационного моделирования морской информационной ситуации при идентификации целей» №2018612944 01.03.2018. Бюл. №3; Свидетельство о государственной регистрации программы для ЭВМ «Специализированный нейросетевой комплекс классификации зашумленных сигналов морских целей» №2018619739 РФ. 10.08.2018. Бюл. №8; Свидетельство о государственной регистрации программы для ЭВМ «Программа проектирования и обучения искусственных нейронных сетей персептронного типа» №2019611559 РФ. 29.01.2019. Бюл. №2).
Сформированный генератором 6.1 сигнал накачки стабилизированной частоты поступает на вход усилителя мощности 6.2, тракта излучения сигналов накачки 6, затем на вход блока согласования 6.3, выход которого соединен с подводным кабелем, соединяющим выход тракта излучения сигналов накачки 6 и вход излучающего преобразователя 1. Излучатель 1 озвучивает среду сигналами накачки стабилизированной частоты в диапазоне десятки-сотни герц. Накачка контролируемой морской среды сложными частотно- или фазово-модулированными сигналами обеспечивает повышение эффективности параметрического приема волн акустически слабозаметных объектов.
На различных режимах движения объекты 5 генерируют излучения, приводящие к изменению величины характеристик проводящей жидкости (плотности и (или) температуры и (или) теплоемкости и т.д.), которые в зависимости от их физической сущности модулируют низкочастотные сигналы накачки морской среды. В спектре нелинейно преобразованных волн (информационных волн) появляются низкочастотные и высокочастотные составляющие, как результат модуляции амплитуды и фазы низкочастотной волны накачки излучениями и полями объектов 5. Являясь неразрывно связанной компонентой просветной волны модуляционные составляющие переносятся на большие расстояния и обнаруживаются в блоках тракта приема, обработки и регистрации информационных сигналов 7.
Информационные волны принимают антенной 2, усиливают в полосе параметрического преобразования (блок 7.1), переносят частотно-временной масштаб сигналов в высокочастотную область (блок 7.2), проводят узкополосный спектральный анализ (блок 7.3) и выделяют из них параметрические составляющие суммарной и разностной частоты, по которым восстанавливают признаки проявления полей, генерируемых объектами.
Операция преобразования частотно-временного масштаба сигнала в блоке 7.2 тракта приема, обработки и регистрации сигналов, обеспечивает увеличение концентрации энергии нелинейно преобразованных сигналов и эффективность выделения из них признаков полей, генерируемых морскими целями.
Операции спектрального анализа в блоке 7.3 позволяют выделять дискретные составляющие суммарной или разностной частоты в узкополосных спектрах нелинейно преобразованных сигналов, по которым восстанавливают характеристики волн объектов 5.
Амплитудно-частотные характеристики сигналов объекта (надводный или подводный объект), полученные с помощью узкополосного спектрального анализа в тракте приема, обработки и регистрации сигналов 7 подают на первый вход блока 8.1 распознавания класса цели по амплитудно-частотным характеристикам тракта нейросетевого распознавания и классификации 8.
На второй вход блока 8.1 распознавания класса цели по амплитудно-частотным характеристикам поступают данные с блока обучения 8.2 тракта нейросетевого распознавания и классификации 8.
Далее по амплитудно-частотным характеристикам распознают класс (надводный или подводный объект) обнаруженной морской цели в тракте нейросетевого распознавания и классификации 8. Анализ низкочастотной, среднечастотной и высокочастотной составляющих амплитудно-частотной характеристики производят раздельно, так как генеральные признаки для различных типов объектов могут находиться в различных частотных диапазонах.
Процесс обучения многослойной нейронной сети (МНС) является итерационным и, в общем случае, достаточно длительным, поскольку заранее нельзя определить количество итераций, необходимых для обучения НС.
ij
Для нейросетевой реализации сравнения расстояний и определения значения K используется следующее выражение:
где ψ(x) - логистическая функция.
Если функция ψ(х) будет дискретной, например, пороговой:
ij
то K будет принимать значение 0 или 1.
Если функция ψ(x) будет вещественной, например, сигмоидной:
ij
то K будет принимать значения на интервале [0, 1], чем ближе значение этой функции будет к 0, тем ближе экземпляр будет к классу, которому сопоставлено значение 0, и, соответственно, наоборот, чем ближе значение этой функции будет к 1, тем ближе экземпляр будет к классу, которому сопоставлено значение 1.
Использование сигмоидной функции может быть более предпочтительным на практике, поскольку она позволяет не только определить к какому классу ближе экземпляр, но и на сколько ближе.
Для вычисления разности расстояний  подставляются соответствующие выражения:
После математических преобразований получаем:
i
j
Выражения для  и  вычисляются на основе формального нейрона, имеющего один вход, на который подается значение x или x, вес которого равен  или  соответственно.
Порог нейрона (нулевой вес) в этом случае будет равен  или  соответственно.
(μ, k)
Правила вычисления параметров алгоритма многомерной классификации в этом случае останутся неизменными, а параметры и функции активации нейронной сети (НС) необходимо определить на их основе по следующим правилам. Функция активации ψ k-го нейрона μ-го слоя:
Весовой коэффициент  p-го входа k-го нейрона μ-го слоя:
Значения выходных сигналов третьего слоя образуют вектор решений  элементы которого представлены в табл. 1. на фиг. 4.
Распознающую сеть блока 8.1 распознавания класса цели по амплитудно-частотным характеристикам тракта нейросетевого распознавания и классификации 8, настраивают по классификационным признакам полей, генерируемых морскими целями, с помощью блока обучения 8.2 тракта нейросетевого распознавания и классификации 8, в который входят априори полученные шумовые портреты вероятных объектов распознавания и устройство подготовки данных.
Корректируют весовые коэффициенты распознающей сети блока 8.1, распознавания класса цели по амплитудно-частотным характеристикам с помощью блока обучения 8.2, тракта нейросетевого распознавания и классификации 8. Далее данные подают в дополнительно введенный тракт регулирования параметров формирования и редукции выборок 9 на вход блока нечетких правил и функций 9.1, где формируют сигнал номера нового продукционного правила и нового типа функции принадлежности типу цели и далее подают на вход логического устройства формирования и редукции выборок 9.3, тракта регулирования параметров формирования и редукции выборок, функцию которого выполняет модифицированная комбинированная распознающая сеть, состоящая из сетей Кохонена - Гросберга и адаптивная нейро-нечеткая сеть ANFIS, использующая гибридный алгоритм обучения, в котором осуществляют компенсацию функций принадлежности нового вида правил, определяют номер правила, необходимого для замены в основной базе правил, а также новый тип функции принадлежности для данного правила, и далее на блок фаззификации 9.2, выход которого соединен с входом блока дефаззификации и автоматического регулирования базы правил 9.4, в котором производят автоподстройку своей базы правил, исходя из выборки образцов математических моделей морских целей, формирование и редукцию выборки эталонных образцов математических моделей морских целей и коррекцию данных оперативно обновляемой библиотеки математически обработанных образов спектрограмм морских целей для блока обучения 8.2 тракта нейросетевого распознавания и классификации 8, после чего искусственную нейронную сеть настраивают и формируют вывод о степени принадлежности исследуемой области спектра объекту классификации (надводный или подводный объект).
Функцию логического устройства формирования и редукции выборок 9.3 в системе нечеткого управления тракта регулирования параметров формирования и редукции выборок 9 выполняет модифицированная комбинированная распознающая сеть, состоящая из сетей Кохонена - Гросберга и адаптивная нейро-нечеткая сеть ANFIS, использующая гибридный алгоритм обучения. Исходная выборка, для работы логического устройства формирования и редукции выборок, будучи отображенной в одномерное пространство, позволит выделить на одномерной оси интервалы ее значений, соответствующие кластерам разных классов в исходном N-мерном пространстве. Определив границы интервалов на одномерной оси, можно найти ближайшие к ним экземпляры, которые и составят формируемую подвыборку.
Алгоритм работы логического устройства формирования и редукции выборок.
Этап инициализации. Задать исходную выборку данных X=<x, y>.
Этап анализа характеристик выборки. Определить  и  - соответственно, минимальное и максимальное значения j-го признака, j=1, 2, …, N. Определить число интервалов для каждого признака: k=K 1n S, а также длины интервалов:
j
Этап расчета обобщенных признаков. Для каждого s-го экземпляра, s=1, 2,…, S: определить k - номер интервала значений по каждому j-му признаку, j=1, 2, …, N, в который попадает s-й экземпляр
рассчитать координату s-го экземпляра по обобщенной оси
Это позволяет отобразить исходную выборку на одномерную обобщенную ось I (заметим, что при этом произойдет потеря части информации вследствие неявного квантования пространства признаков при преобразовании).
s
s
s
s
q
q
q
q
q
I
Этап анализа обобщенной оси. Сформировать набор кортежей I={<I, y, s>}. Упорядочить набор I в порядке неубывания значений I. Просматривая обобщенную ось в порядке увеличения ее значений, определить граничные значения ее интервалов <l,r>, в которых номер класса y остается неизменным, где l, r - соответственно левое и правое граничные значения q-го интервала обобщенной оси. Обозначим: K - номер класса, соответствующий q-му интервалу обобщенной оси; k - число интервалов обобщенной оси.
q
I
q
Этап анализа характеристик интервалов. Для каждого q-го интервала обобщенной оси, k=1, 2, …, k, определить S - число попавших в него экземпляров, а также номера этих экземпляров.
*
Этап формирования обучающей выборки. Среди экземпляров q-го интервала включить в обучающую выборку X все экземпляры: его класса, находящиеся на одной из границ интервала
его класса, ближайшие к одной из границ интервала:
q
q
где α - пороговый коэффициент, регулирующий близость экземпляров к границам интервала (например, можно задать: α=0,l(r-l));
интервалов с малым числом экземпляров:
i
s=1, 2, …, S, q=1, 2, …, k,
где β - некоторый пороговый коэффициент, 0<β<1 (например, можно задать: β=0,1);  - среднее число экземпляров в интервале обобщенной оси.
Этап устранения избыточности обучающей выборки. Определить расстояния между всеми экземплярами, вошедшими в сформированную обучающую выборку, сформировав матрицу расстояний R (для упрощения и ускорения вычислений будем оперировать квадратами расстояний):
Заметим, что R(s, р)=R(p, s), a R(s, s)=0.
До тех пор, пока , g≠р, выполнять в цикле действия: найти в матрице расстояний два экземпляра с наименьшим расстоянием между собой:
если два ближайших экземпляра принадлежат к одному и тому же классу, то оставить в обучающей выборке только тот из них, который находится ближе к экземплярам других классов, а другой исключить из нее:
*
Скорректировать соответствующим образом элементы матрицы R, установив: R(s, q)=R(q, s)=-1; если два ближайших экземпляра принадлежат к разным классам, то перейти к выполнению этапа дополнения (уточнения) обучающей выборки. Этап дополнения (уточнения) обучающей выборки. Определить разность исходной и сформированной выборок X'=Х/Х.
s'
s'
*
*
Последовательно для каждого s'-го экземпляра <x, y> выборки X', s'=1, 2, …, S' относительно экземпляров сформированной выборки X найти расстояние (квадрат расстояния) от него до каждого экземпляра выборки X:
*
если ближайший к s'-му экземпляру экземпляр сформированной выборки принадлежит к другому классу, то включить его в выборку X:
В результате выполнения данного метода для исходной выборки X получаем сформированную обучающую выборку X', которая будет обладать основными топологическими свойствами исходной выборки. При этом также из исходной выборки можно получить также тестовую выборку как разность исходной и сформированной обучающей выборок. Разбиение признакового пространства для выборки эмпирических наблюдений необходимо для определения нечетких термов признаков как проекций соответствующих блоков на координатные оси. Формирование разбиения предлагается осуществлять путем их выполнения в приведенной ниже последовательности.
Шаг 1. Инициализация. Задать обучающую выборку <x, y>.
Шаг 2. По оси каждого признака j=1, 2, …, N определить одномерные расстояния между экземплярами: .
Среди полученных расстояний найти минимальное расстояние, большее нуля:
Шаг 3. Для каждого признака определить количество интервалов разбиения диапазона его значений:
j
j
j
а также определить длину интервала наблюдаемых значений каждого признака: r=max(x)-min(x), j=1, 2, …, N.
j
Шаг 4. Разбить ось j-го признака на n интервалов. Определить координаты левых и правых границ для каждого p-го интервала j-го признака по формулам:
Шаг 5. Сформировать блоки-кластеры и задать номера их классов путем выполнения шагов 5.1-5.8.
q, j
Шаг 5.1. Сформировать прямоугольные блоки  в N-мерном пространстве признаков на пересечении соответствующих интервалов значений признаков. Занести в B номер интервала j-го признака, который соответствует q-му блоку.
Шаг 5.2. Определить номера классов для прямоугольных блоков в N-мерном пространстве признаков:
Установить коэффициент уверенности классификации для блоков:
q
Шаг 5.3. Для тех блоков, у которых K=-1, q=G+1, G+2, …,
G+Q, установить:
где  - количество экземпляров k-го класса, попавших в q-й блок-кластер.
q
Шаг 5.4. Для тех блоков, у которых номер класса K=0, q=G+1, G+2, …, G+Q, определить расчетный номер класса, для чего предлагается использовать модифицированный нерекуррентный метод потенциальных функций.
Шаг 5.5. Вычислить расстояние между q-м и р-м блоками, q=G+1, G+2, …, G+Q, p=q+1, …, G+Q, как:
q
p
p
q
Или  где  При этом R(B, B)=R(B, B).
Шаг 5.6. Определить потенциал, наводимый совокупностью блоков, принадлежащих к k-му классу, на р-й блок с неизвестной классификацией:
q=G+1, G+2,..., G+Q, р=q+1.....G+Q.
k
q
где L - количество блоков, принадлежащих к к-му классу, S - количество экземпляров обучающей выборки, попавших в q-й блок.
р
Шаг 5.7. Установить номер класса для р-го блока с неизвестной классификацией (К=0) по формуле:  р=G+1, G+2, …, G+Q.
Шаг 5.8. Модифицировать значения коэффициентов уверенности для блоков: q =G+1, G+2, …, G+Q.
Шаг 6. Выполнить объединение смежных блоков-кластеров.
q
q
p
q, i
p, i
Выполнить объединение смежных блоков, принадлежащих к одному и тому же классу: для ∀q, р=G+1, G+2, …, G+Q; q≠р: если K>0, K=K и  ∀i≠j:B=B, i=1, 2, …, N, j=1, 2, …, N; тогда объединить блоки q и p по j-му признаку:
- установить:
p
p
p, i
- удалить р-й блок: K=0, α=0, B=0, i=1, 2, …, N.
*
*
*
*
*
Шаг 7. Из объединенного множества (ОМ) выделить подмножество экземпляров, относящихся к блокам-кластерам, номера классов которых не совпадают с номерами классов экземпляров. Применить для полученного разбиения и выделенного подмножества процедуру уточнения разбиения и дообучения модели. Шаг 8. Остановить. Уточнение разбиения и дообучение модели. Если имеется разбиение признакового пространства, которое нужно уточнить (дообучить) на основе новых наблюдений <х, y>,  j=1, 2, …, N; s=1, 2, …, S; то необходимо из набора новых наблюдений исключить те наблюдения, которые попадают в блоки имеющегося разбиения и соответствуют им по номеру класса, скорректировав соответствующим образом S. Для тех наблюдений, которые не совпадают с классами блоков, целесообразно сформировать отдельные точечные кластеры.
Для каждого нового наблюдения сформировать интервалы по признакам и занести в  номера интервалов для каждого j-го признака, соответствующие новому кластеру, а также определить:
где γ - некоторая константа, γ∈.
Для определения целесообразности применения предложенного метода для конкретной задачи на практике используем нотацию Ландау в так называемом «мягком виде» и оценим сложность этапов предложенного метода. Для этапа инициализации вычислительной сложностью можно пренебречь, а пространственная сложность может быть оценена как O(NS).
2
4
2
2
Для этапа анализа характеристик выборки вычислительная сложность составит O(2NS), а пространственная - O(4N). Для этапа расчета обобщенных признаков вычислительная сложность может быть оценена как O(9NS), а пространственная - O(N+S). Для этапа анализа обобщенной оси вычислительная сложность может быть оценена как O(2S+2S), а пространственная - O(3S+3K1nS). Для этапа анализа характеристик интервалов вычислительная сложность может быть оценена как O(3KS1nS), а пространственная - O(K1nS). Для этапа формирования обучающей выборки вычислительная сложность может быть оценена как O(20KS1nS), а пространственная в виде O(0,2K1nS+0,2S). Для этапа устранения избыточности ОВ вычислительная и пространственная сложность могут быть оценены соответственно как O(0,0016S+(0,44+0,08N)S) и O(0,04S).
2
4
2
2
2
Для этапа дополнения (уточнения) обучающей выборки вычислительная и пространственная сложность оцениваются соответственно как O((0,48N+0,32)S) и O(0,8S). Таким образом, общая сложность метода может быть оценена как: вычислительная - O(0,0016S+0,56NS+2.76S+11NS+2S++23KS1nS); пространственная - O(0,04S+NS+5N+5S+4,2K1nS).
i
j
j
k
Алгоритм коррекции значений настроек регулятора тракта регулирования параметров формирования и редукции выборок, реализован для нахождения номера нечеткого правила вида: (ЕСЛИ Tpi Условие И G Условие ТО N «Заключение»)
k=1, 2, …, N,
1
1
N
N
1
1
N
N
при наличии обучающего множества ((D, N), …, (D, N)), ((D, ФП), …, (D, ФП)).
i:
1
i1
р1
i2
1
in
i,
ij
i
Для моделирования неизвестного отображения f используется алгоритм нечеткого вывода, применяются предикатные правила: П ЕСЛИ е есть A И Т есть A И G есть A, ТО N=S i=1, 2, …, m, где A - нечеткие множества описывающие высказывания: «отрицательная», «нулевая», «положительная», «малое», «среднее», «большое» и т.д. S - вещественные числа (номер правила). Степень истинности μ правила i - определяется с помощью операции умножения Ларсена (Larsen):
k
k
k
k
i
для моделирования логического оператора «И», выход нечеткой системы Uопределяется методом центра тяжести:  Функция ошибки для i-го предъявленного значения вида: E=0.5(U-N) позволяет использовать градиентный метод для подстройки параметров заданных предикатных правил, а величина S корректируется по соотношению:
i=1, 2, …, m,
где η - константа, характеризующая скорость обучения сети. Аналогичным образом определяются параметры функции принадлежности.
В ситуации неопределенности, т.е. влияния случайных внешних и параметрических возмущений, блок дефаззификации и автоматического регулирования базы правил 9.4 осуществляет их компенсацию в базе правил, а также компенсацию функций принадлежности нового вида (с другими универсумом Е и E'.
Таким образом, тракт регулирования параметров формирования и редукции выборок определяет номер правила (N), необходимого для замены в основной базе правил, а также новый тип функции принадлежности для данного правила. Блок дефаззификации и автоматического регулирования базы правил 9.4 тракта регулирования параметров формирования и редукции выборок самостоятельно производит автоподстройку своей базы правил, исходя из выборки образцов математических моделей морских целей, осуществляя формирование и редукцию выборки эталонных образцов математических моделей морских целей и коррекцию данных оперативно обновляемой библиотеки математически обработанных образов спектрограмм морских целей тракта нейросетевого распознавания и классификации 8, обеспечивающего конечное классификационное решение по обнаруженным морским целям (надводный или подводный объект).
Поскольку базы нечетких правил часто характеризуются большим объемом, то наиболее актуальной является задача объединения и преобразования нечетких правил, суть которой заключается в том, чтобы на основе начального набора нечетких правил сформировать новую базу правил меньшего объема, которая в достаточной мере представляла бы начальную базу и была бы менее избыточной.
Для решения задачи формирования и редукции нечетких правил необходимо привести этапы метода, в соответствие с особенностями решаемой задачи, который представлен в виде последовательности шагов 1-17.
Шаг 1. Инициализация. Задаются статические параметры работы метода: коэффициенты α, β, p. Для каждого из возможных классов выходных значений создается свое пространство поиска, и, соответственно, свое отдельное множество агентов, отдельный граф поиска, представляющий собой лингвистические термы, которые могут быть включены в правила, а также для каждого пространства поиска рассчитываются эвристические значения значимости отдельного терма.
p
q
init
где  - значение эвристической значимости лингвистического терма p для описания класса q; о - экземпляр входной выборки, содержащей N экземпляров; μ(o), μ(o) - значение функции принадлежности объекта о терму р и классу q, соответственно; Т - количество лингвистических термов; K - количество классов. В каждом пространстве поиска каждому узлу графа поиска ставится в соответствие начальное значение количества дискретных составляющих (ДС) τгде  - значение количества ДС для р-го терма в пространстве поиска для q-го класса на первой итерации поиска.
Шаг 2. Установить: t=1. Шаг 3. Установить: i=1. Шаг 4. Установить: j=1. Шаг 5. Установить: k=1.
Шаг 6. Выбор терма для добавления в правило j-го агента в пространстве поиска i-го класса.
j
Шаг 6.1. Для j-го агента на основе правила случайного выбора рассчитывается вероятность включения k-го лингвистического терма в правило, описывающего i-й класс выходного значения  где  - вероятность добавления k-то терма в правило j-го агента в пространстве поиска для i-го класса; R - множество термов, которые могут быть добавлены в правило j-го агента. Формирование данного множества определяет вид правил, которые могут составляться в процессе поиска, то есть предполагается, что правило может включать выражения типа ИЛИ. После добавления терма из данного множества исключается только данный терм, если же предполагается, что правило не может включать выражения типа ИЛИ, то кроме выбранного терма, исключаются и все термы, описывающие данный атрибут.
Шаг 6.2. Проверить условие  где rand (1) - случайное число из интервала [0; 1]. Если условие выполняется, тогда лингвистический терм k добавляется в правило j-го агента, удаляется из множества возможных термов для данного агента и выполняется переход к шагу 7. В противном случае - переход к шагу 6.3.
Шаг 6.3. Установить k=k+1.
Шаг 6.4. Если были рассмотрены все термы, то установить: k=1. Переход к шагу 6.1.
Шаг 7. Проверка завершения перемещения j-го агента.
Шаг 7.1. Если множество термов, которые j-й агент может добавить в формируемое правило, пусто, то выполняется переход к шагу 8.
Шаг 7.2. Определяется, сколько экземпляров i-го класса покрывает правило j-го агента.
j
j
j
j
Шаг 7.2.1. Для экземпляра о, относящегося к классу i, рассчитывается степень соответствия сформированного правила R и экземпляра о  где match(R, o) - степень соответствия между правилом j-го агента R и экземпляромо;  - мера соответствия между р-м атрибутом в правиле R и соответствующим атрибутом экземпляра о. Данная мера рассчитывается следующим образом:
p
где q - отдельный терм, относящийся к области описания атрибута р; Q - количество термов, относящихся к области описания атрибута р.
j
j
j
Шаг 7.2.2. Проверить условие: match(R, о)≥inMatchMin, где inMatchMin - заданный параметр, который определяет, какое минимальное значения соответствия является достаточным, чтобы считать, что правило R в достаточной степени описывает объект о. Если условие выполняется, то считается, что данный объект о покрывается правилом R.
j
Шаги 7.2.1 и 7.2.2 выполняются для всех экземпляров, относящихся к классу i, и на основании получаемых данных увеличивается счетчик cntMatch, в котором хранится количество экземпляров, покрываемых правилом R.
Шаг 7.3. Проверить условие: cntMatch≥nCntMatchMini, где inCntMatchMini - предельное минимальное количество экземпляров i-го класса, которое должно покрываться правилом. Если указанное условие выполняется, то считается, что правило покрывает необходимое количество экземпляров, и j-й агент завершил свое перемещение, выполняется переход к шагу 8. В противном случае - выполняется переход к шагу 5.
Шаг 8. Если j<cntAgents, то установить: j=j+1 и выполнить переход к шагу 5. В противном случае - переход к шагу 9.
Шаг 9. Если i<k, то установить: i=i+1 и выполнить переход к шагу 4. В противном случае - переход к шагу 10.
Шаг 10. Случайным образом формируются базы правил. Создастся nBases баз правил, при этом для описания каждого класса выходного значения выбирается одно правило из соответствующего пространства поиска.
Шаг 11. Оценка качества сформированных баз правил. Для оценки качества баз правил используется входная обучающая выборка, для каждого экземпляра которой выбирается правило с самой высокой степенью совпадения, на основании которого определяется расчетный класс, к которому относится данный экземпляр, основываясь на соответствующей базе правил. Основываясь на данных, полученных при помощи базы знаний, исходя из заданной обучающей выборки, рассчитываем оценку качества базы правил:  где cntNotMatch - количество экземпляров, для которых класс был определен неверно с помощью заданной базы правил; Q - качество прогнозирования класса экземпляров на основе соответствующей базы правил.
high
threshold
high
threshold
Шаг 12. Проверить условие: Q≥Q, где Q - качество прогнозирования базы знаний, которая характеризуется наилучшей точностью прогнозирования; Q - приемлемое качество прогнозирования. Если указанное условие выполняется, то производится переход к шагу 17, в противном случае - переход к шагу 13.
Ω
high
Ω
high
_
Шаг 13. Добавление ДС. Добавление ДС выполняется с целью повышения приоритетности тех термов, включение которых в правила способствует повышению качеству прогнозирования результирующих баз правил. В связи с этим количество добавляемого коэффициента приоритетности прямо пропорционально качеству прогнозирования базы правила, в которое входит заданное нечеткое правило. При этом добавление ДС предлагается выполнять только для тех термов, входящих в правила баз нечетких правил, для которых выполняется условие: Q≥δ⋅Q, где δ коэффициент, определяющий, насколько близко качество прогнозирования базы правил Q должно приближаться к лучшему качеству прогнозирования Q, чтобы можно было применять процедуру добавления ДС для правил, входящих в данную базу правил Ω. Таким образом, добавление ДС выполняется для каждого терма, входящего в правило, которое, в свою очередь, входит в базу правил Ω.  ∀p∈R, ∀R⊂Ω, где  - количество ДС для терма р в пространстве поиска для класса q, который определяется с помощью соответствующего правила.
Шаг 14. Исключение ДС. Для исключения худших термов, то есть таких, которые при включении их в правила, понижают качество прогнозирования с помощью соответствующего правила, применяют процедуру исключения ДС, которая выполняется в конце каждой итерации и применяется для всех узлов во всех графах поиска. Исключение ДС морской цели (МЦ) выполняется в соответствии с формулой:  ,,  где ρ - коэффициент исключения, который задается при инициализации.
max
Шаг 15. Если t<t, тогда установить: t=t+1 и выполнить переход к шагу 16, в противном случае - считается, что выполнено максимально допустимое количество итераций, и выполняется переход к шагу 17.
Шаг 16. Перезапуск агентов. Все данные обновляются, агенты размещаются в случайные точки пространств поиска. Переход к шагу 3. Шаг 17. Останов.
На основе проведенных экспериментов были получены базы нечетких правил, которые характеризовались точностью классификации МЦ данным методом в 85,6%. Следовательно, предложенный метод обеспечивает формирование редуцированной базы нечетких правил, которая характеризуется достаточно точными результатами, сокращая при этом временные затраты на расчет. Существенным достоинством данного метода, является то, что поиск производится одновременно на нескольких графах, что позволяет контролировать влияние термов на качество прогнозирования распознавания разных классов МЦ.
В процессе работы способа постоянно происходит коррекция данных оперативно обновляемой библиотеки математически обработанных образов спектрограмм морских целей, статистики классификационной характеристики, уточняется описание классов морских целей для текущих гидролого-акустических и помехосигнальных условий, на основании чего формируется актуальный порог для классификации.
Таким образом, обнаружив цель по признакам амплитудно-фазовой модуляции низкочастотных сигналов накачки морской среды излучениями и полями объекта, и используя для регулирования параметров формирования и редукции выборок адаптивную нейро-нечеткую коррекцию с комплексным сокращением размерности данных при автоподстройке базы правил за счет формирования и редукции выборки эталонных образцов оперативно обновляемой библиотеки математически обработанных образов спектрограмм морских целей, а также архитектуру распознающей нейронной сети в виде модифицированной комбинированной распознающей сети, состоящей из сетей Кохонена (конкурирующая сеть) и Гросберга и адаптивной нейро-нечеткой сети ANFIS, использующей гибридный алгоритм обучения, можно в автоматизированном режиме распознавать класс цели по амплитудно-частотным характеристикам и делать вывод о степени принадлежности исследуемой области спектра объекту классификации (надводный или подводный объект). По результатам экспериментальных данных, эффективность классификации морских целей в заявленном способе увеличивается на 7-9% относительно прототипа и других способов классификации, использующих модуляции низкочастотных сигналов накачки морской среды излучениями и полями объектов.
Способ промышленно применим, так как для его реализации используются распространенные компоненты и изделия радиотехнической промышленности и вычислительной техники.