Классификация сейсмических сигналов на основе нейросетевых технологий



РУБРИКИ	Классификация сейсмических сигналов на основе нейросетевых технологий	РЕКОМЕНДУЕМ

Главная

Историческая личность

История

Искусство

Литература

Москвоведение краеведение

Авиация и космонавтика

Административное право

ПОИСК

Классификация сейсмических сигналов на основе нейросетевых технологий

Sj – компонента целевого вектора (S={0,1} – взрыв, S={1,0}-землетрясение)

Таким образом после предварительного обучения и формирования кластеров

в слое Кохонена, на фазе вторичного обучения все нейроны каждого

полученного кластера соединяются активными (единичными) синапсами со своим

выходным нейроном, характеризующим данный кластер.

Выход нейронов второго слоя определяется выражением:

(11) [pic]

где: [pic]

Kj - размерность j-ого кластера, т.е. количество нейронов слоя Кохонена

соединённых с нейроном j выходного слоя отличными от нуля

коэффициентами.

R - пороговое значение (0 < r=""><>

Пороговое значение R можно выбрать таким образом, чтобы с одной

стороны не были потеряны значения активированных кластеров, а с другой

стороны - отсекался "шум не активизированных кластеров".

В результате на каждом шаге обработки исходных данных на выходе

получаются значения Sj, которые характеризуют явление, породившее данную

входную ситуацию ([pic] - землетрясение; [pic] - взрыв).

4.5 Выводы по разделу.

Итак, подводя итог данной главе, следует сказать, что это далеко не

полный обзор нейросетевых архитектур, которые успешно справляются с

задачами классификации. В частности ничего не было сказано о вероятностных

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

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

менее известных, но хорошо себя зарекомендовавших. Соответственно проблема

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

сейсмических сигналов вполне актуальна. В идеале, конечно хотелось бы

проверить эффективность хотя бы нескольких из них и выбрать наилучшую. Но

для этого необходимо проводить более масштабные исследования, которые

займут много времени. На данном этапе исследований была сделана попытка

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

чтобы убедиться в эффективности подхода в целом. В главе 6 детально

обсуждаются нейросеть, которая была исследована в рамках настоящей

дипломной работы.

5. Методы предварительной обработки данных.

Если возникает необходимость использовать нейросетевые методы для

решения конкретных задач, то первое с чем приходится сталкиваться – это

подготовка данных. Как правило, при описании различных нейроархитектур, по

умолчанию предполагают что данные для обучения уже имеются и представлены в

виде, доступном для нейросети. На практике же именно этап предобработки

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

обучения нейросети также может решающим образом зависеть от того, в каком

виде представлена информация для ее обучения.

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

понижения размерности исходных данных, позволяющие увеличить

информативность обучающей выборки.

5.1 Максимизация энтропии как цель предобработки.

Рассмотрим основной руководящий принцип, общий для всех этапов

предобработки данных. Допустим, что в исходные данные представлены в

числовой форме и после соответствующей нормировки все входные и выходные

переменные отображаются в единичном кубе. Задача нейросетевого

моделирования – найти статистически достоверные зависимости между входными

и выходными переменными. Единственным источником информации для

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

больше бит информации принесет пример – тем лучше используются имеющиеся в

нашем распоряжении данные.

Рассмотрим произвольную компоненту нормированных (предобработанных)

данных: [pic]. Среднее количество информации, приносимой каждым примером

[pic], равно энтропии распределения значений этой компоненты [pic]. Если

эти значения сосредоточены в относительно небольшой области единичного

интервала, информационное содержание такой компоненты мало. В пределе

нулевой энтропии, когда все значения переменной совпадают, эта переменная

не несет никакой информации. Напротив, если значения переменной [pic]

равномерно распределены в единичном интервале, информация такой переменной

максимальна.

Общий принцип предобработки данных для обучения, таким образом состоит

в максимизации энтропии входов и выходов.

5.2 Нормировка данных.

Как входами, так и выходами могут быть совершенно разнородные величины.

Очевидно, что результаты нейросетевого моделирования не должны зависеть от

единиц измерения этих величин. А именно, чтобы сеть трактовала их значения

единообразно, все входные и выходные величин должны быть приведены к

единому масштабу. Кроме того, для повышения скорости и качества обучения

полезно провести дополнительную предобработку, выравнивающую распределения

значений еще до этапа обучения.

Индивидуальная нормировка данных.

Приведение к единому масштабу обеспечивается нормировкой каждой

переменной на диапазон разброса ее значений. В простейшем варианте это –

линейное преобразование:

[pic]

в единичный отрезок: [pic]. Обобщение для отображения данных в интервал

[pic], рекомендуемого для входных данных тривиально.

Линейная нормировка оптимальна, когда значения переменной [pic] плотно

заполняют определенный интервал. Но подобный «прямолинейный» подход

применим далеко не всегда. Так, если в данных имеются относительно редкие

выбросы, намного превышающие типичный разброс, именно эти выбросы определят

согласно предыдущей формуле масштаб нормировки. Это приведет к тому, что

основная масса значений нормированной переменной [pic] сосредоточится

вблизи нуля [pic] Гораздо надежнее, поэтому, ориентироваться при нормировке

не а экстремальные значения, а на типичные, т.е. статистические

характеристики данных, такие как среднее и дисперсия.

[pic], где

[pic], [pic]

В этом случае основная масса данных будет иметь единичный масштаб, т.е.

типичные значения все переменных будут сравнимы (рис. 6.1)

Однако, теперь нормированные величины не принадлежат гарантированно

единичному интервалу, более того, максимальный разброс значений [pic]

заранее не известен. Для входных данных это может быть и не важно, но

выходные переменные будут использоваться в качестве эталонов для выходных

нейронов. В случае, если выходные нейроны – сигмоидные, они могут принимать

значения лишь в единичном диапазоне. Чтобы установить соответствие между

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

диапазон изменения переменных.

Линейное преобразование, представленное выше, не способно отнормировать

основную массу данных и одновременно ограничить диапазон возможных значений

этих данных. Естественный выход из этой ситуации – использовать для

предобработки данных функцию активации тех же нейронов. Например,

нелинейное преобразование

[pic], [pic]

нормирует основную массу данных одновременно гарантируя что[pic] (рис. 5.2)

Как видно из приведенного выше рисунка, распределение значений после

такого нелинейного преобразования гораздо ближе к равномерному.

Все выше перечисленные методы нормировки направлены на то, чтобы

максимизировать энтропию каждого входа (выхода) по отдельности. Но, вообще

говоря, можно добиться гораздо большего максимизируя их совместную

энтропию. Существуют методы, позволяющие проводить нормировку для всей

совокупности входов, описание некоторых из них приведено в [4].

6.3 Понижение размерности входов.

Поскольку заранее неизвестно насколько полезны те или иные входные

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

число входных параметров, в надежде на то, что сеть сама определит, какие

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

результатам, а к тому же еще и увеличивает сложность обучения. Напротив,

сжатие данных, уменьшение степени их избыточности, использующее

существующие в них закономерности, может существенно облегчить последующую

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

алгоритмов, предназначенных для понижения размерности данных с минимальной

потерей информации:

. Отбор наиболее информативных признаков и использование их в процессе

обучения нейронной сети;

. Кодирование исходных данных меньшим числом переменных, но при этом

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

данных.

Рассмотрим более подробно оба типа алгоритмов.

5.3.1 Отбор наиболее информативных признаков.

Для того, чтобы понять какие из входных переменных несут максимум

информации, а какими можно пренебречь необходимо либо сравнить все признаки

между собой и определить степень информативности каждого из них, либо

пытаться найти определенные комбинации признаков, которые наиболее полно

отражают основные характеристики исходных данных.

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

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

позволяют применять его для более распространенных задач.

Для выбора подходящей комбинации входных переменных используется так

называемые генетические алгоритмы [5], которые хорошо приспособлены для

задач такого типа, поскольку позволяют производить поиск среди большого

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

5.3.2 Сжатие информации. Анализ главных компонент.

Самый распространенный метод понижения размерности - это анализ главных

компонент (АГК).

Традиционная реализация этого метода представлена в теории линейной

алгебры. Основная идея заключается в следующем: к данным применяется

линейное преобразование, при котором направлениям новых координатных осей

соответствуют направления наибольшего разброса исходных данных. Для эти

целей определяются попарно ортогональные направления максимальной вариации

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

размерности, порожденное компонентами с наибольшей вариацией [4]. Один из

недостатков классического метода главных компонент состоит в том, что это

чисто линейный метод, и соответственно он может не учитывать некоторые

важные характеристики структуры данных.

В теории нейронных сетей разработаны более мощные алгоритмы,

осуществляющие “нелинейный анализ главных компонент”[3]. Они представляют

собой самостоятельную нейросетевую структуру, которую обучают выдавать в

качестве выходов свои собственные входные данные, но при этом в ее

промежуточном слое содержится меньше нейронов, чем во входном и выходном

слоях. (рис 5.3). Сети подобного рода носят название – автоассоциативные

сети.

Чтобы восстановить свои входные данные, сеть должна научиться

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

этом случае носит название правило обучения Ойя для однослойной сети.

Учитывая то, что в такой структуре веса с одинаковыми индексами в обоих

слоях одинаковы, дельта-правило обучения верхнего (а тем самым и нижнего)

слоя можно записать в виде:

[pic], где

[pic] [pic],и [pic] [pic],

[pic], j=1,2,…,d – компонента входного вектора;

[pic], выходы сети j=1,…,d;

d - количество нейронов на входном ми выходном слоях (размерность вектора

признаков);

yi - выход с i-го нейрона внутреннего слоя, i=1,…,M

M – количество нейронов на внутреннем слое;

. - коэффициент обучения;

wij=wkj - веса сети , соответственно между входным – скрытым и скрытым –

выходным слоями.

Скрытый слой такой сети осуществляет оптимальное кодирование входных

данных, и содержит максимально возможное при данных ограничениях количество

информации. После обучения внешний интерфейс (wij) (рис.5.4) может быть

сохранен и использован для понижения размерности.

Нелинейный анализ главных компонент.

Главное преимущество нейроалгоритмов в том, что они легко обобщаются на

случай нелинейного сжатия информации, когда никаких явных решений уже не

существует. Можно заменить линейные нейроны в описанных выше сетях –

нелинейными. С минимальными видоизменениями нейроалгоритмы будут работать и

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

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

на сигмоидную в правиле обучения Ойя:

[pic]

приводит к новому качеству.

Таким образом, нейроалгоритмы представляют собой удобный инструмент

нелинейного анализа, позволяющий относительно легко находить способы

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

5.4 Выводы по разделу.

Конечно, описанными выше методиками не исчерпывается все разнообразие

подходов к ключевой для нейро-анализа проблеме формирования пространства

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

главных компонент. Также, большего внимания заслуживают генетические

алгоритмы. Необъятного не объять. Главное, чтобы за деталями не терялся

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

избыточности всеми возможными способами. Это повышает информативность

примеров и, тем самым, качество нейропредсказаний.

6. Реализация нейросетевой модели и исследование ее технических

характеристик.

Ранее было показано, какими средствами нейроинформатики можно пытаться

решить задачу идентификации типа сейсмического источника, какие процедуры

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

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

этой главе представлено решение задачи на базе двухслойного персептрона,

так как именно он был выбран на начальном этапе исследований. Дано также

описание алгоритма обучения и методов его оптимизации.

6.1 Структура нейронной сети.

Итак, для решения задачи идентификации типа сейсмического события

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

многослойный персептрон, а точнее его двухслойную реализацию (рис. 6.1).

Как показали эксперименты, увеличение числа скрытых слоев не приводит к

лучшим результатам, а лишь усложняет процесс обучения, поэтому и была

выбрана именно реализация с одним скрытым слоем нейронов.

На вход сети подается p-мерный вектор признаков {xi, i=1,2,…,p}. Для

определенности будем рассматривать случай, когда p=9, хотя исследования

проводились и для p=5, p=18. Оптимальное количество нейронов на скрытом

слое H подбиралось экспериментально для разных p. Соответственно при p = 9

достаточно брать H равным также 9 или немного больше. Для разбиения

исходных данных на два класса на выходе сети достаточно одного нейрона.

Между входным и скрытым слоями, а также между скрытым и выходным слоями

использовалась полносвязная структура.

С учетом этих дополнений опишем принятые на рисунке 7.1 обозначения:

p – размерность исходных данных (количество признаков используемых для

классификации);

H – число нейронов на скрытом слое;

xi – компонента входного вектора признаков, i = 1,…,p;

x0 ( 1 – постоянное воздействие используемое для работы нейронной сети;

wji – весовые коэффициенты между входным и скрытым слоями, i = 0,1,…,p ,

j = 1,…,H;

vk - весовые коэффициенты между скрытым и выходным слоями, k = 0,1,…,H.

zj – значение выхода j-го нейрона скрытого слоя; z0 ( 1, j = 1,…,H;

[pic]

y – значение выходного нейрона сети (выход сети)

[pic] (12)

f1(x) –функция активации нейронов скрытого слоя;

f2(x) –функция активации нейрона выходного слоя.

В качестве функции активации f1(x) для нейронов скрытого слоя и f2(x)

для единственного нейрона на выходе сети предлагается использовать одну и

ту же функцию, а именно сигмоидную функцию активации, для краткости будем

обозначать ее как f(x):

[pic],

с производной в виде

[pic].

Вид такой функции представлен на рис.6.2

Т.к. значения функции f(x) ограничены в диапазоне [0, 1], результат

сети y(x) может принимать любые действительные значения из этого же

диапазона, в следствии чего логично интерпретировать выходы сети следующим

образом: если y(x) > 0.5, то вектор принадлежит к одному классу (взрывы), в

противном случае к другому (землетрясения).

6.2 Исходные данные.

На вход нейронной сети предлагается подавать вектора признаков

составленные из сейсмограмм. О том, какие признаки были использованы для

этой задачи и как они получены, было рассказано ранее в разделе 3.1. Стоит

отметить, что проблема формирования векторов признаков – это исключительно

проблема сейсмологии. Поэтому для исследования эффективности применения

нейронных сетей в качестве исходных данных были использованы уже готовые

выборки векторов, которые содержали в себе примеры и землетрясений и

взрывов.

Размерность векторов признаков p=9, хотя , как было отмечено в

предыдущем разделе, проводились эксперименты и с другим количеством

признаков.

Для работы с нейросетью рекомендуется использовать исходные данные не в

первоначальном виде, а после предварительной обработки при помощи процедуры

индивидуальной нормировки по отдельному признаку, описанной в разделе 5.2.

Это преобразование состоит в следующем:

[pic]

где

xi – исходное значение вектора признаков, точнее его i-я компонента;

xi,min – минимальное значение по i-му признаку, найденное из всей

совокупности исходных данных, включающей оба класса событий;

xi,max – максимальное значение по i-му признаку …

Выбор именно этой нормировки, а не более универсальных, которые описаны

в разделе 5, в настоящих исследованиях продиктованы тем обстоятельством,

что непосредственно признаки измеренные по сейсмограммам, подвергаются

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

y=Ln(x) и z=(1/7)(y1/7-1),

и уже из этих значений формируются обучающие вектора. Такие преобразования

приводят к большей кластеризации точек в многомерном пространстве, однако

диапазон изменения каждого из признаков не нормирован относительно

интервала [-1, 1], а выбранная нормировка позволяет без потери информации

перенести все входные значения в нужный диапазон.

6.3 Определение критерия качества системы и функционала его оптимизации.

Если через [pic] обозначить желаемый выход сети (указание учителя), то

ошибка системы для заданного входного сигнала (рассогласование реального и

желаемого выходного сигнала) можно записать в следующем виде:

[pic][pic], где

k — номер обучающей пары в обучающей выборке, k=1,2,…,n1+n2

n1 - количество векторов первого класса;

n2 - число векторов второго класса.

В качестве функционала оптимизации будем использовать критерий минимума

среднеквадратической функции ошибки:

[pic]

6.4 Выбор начальных весовых коэффициентов.

Перед тем, как приступить к обучению нейронной сети, необходимо задать

ее начальное состояние. От того насколько удачно будут выбраны начальные

значения весовых коэффициентов зависит, как долго сеть за счет обучения и

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

Как правило, всем весам на этом этапе присваиваются случайные величины

равномерно распределенные в диапазоне [-A,A], например [-1,1], или [-3,3].

Однако, как показали эксперименты, данное решение не является наилучшим и в

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

инициализации, а именно:

1. Присваивать весам случайные величины, заданные не равномерным

распределением, а нормальным распределением с параметрами N[(,(], где

выборочное среднее (=0, а дисперсия ( = 2, или любой другой небольшой

положительной величине. Для формирования нормально распределенной

величины можно использовать следующий алгоритм:

Шаг 1. Задать 12 случайных чисел x1, x2, …,x12 равномерно

распределенных в диапазоне [0,1]. xi ( R[0,1].

Шаг 2. Для искомых параметров ( и ( величина [pic], полученная по

формуле:

[pic]

будет принадлежать нормальному распределению с параметрами N[(,(].

2. Можно производить начальную инициализацию весов в соответствии

с методикой, предложенной Nguyen и Widrow [7]. Для этой методики

используются следующие переменные

[pic] число нейронов текущего слоя

[pic] количество нейронов последующего слоя

[pic] коэффициент масштабирования:

[pic]

Вся процедура состоит из следующих шагов:

Для каждого нейрона последующего слоя[pic]:

Инициализируются весовые коэффициенты (с нейронов текущего слоя):

[pic]случайное число в диапазоне [-1,1] ( или [pic]).

Вычисляется норма [pic]

Далее веса преобразуются в соответствии с правилом:

[pic]

Смещения [pic] выбираются случайным образом из диапазона [pic].

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

результатов, в сравнении со стандартным алгоритмом начальной инициализации

весов.

6.5 Алгоритм обучения и методы его оптимизации.

Приступая к обучению выбранной нейросетевой модели, необходимо было

решить, какой из известных типов алгоритмов, градиентный (обратное

распространения ошибки) или стохастический (Больцмановское обучение)

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

подход, который и представлен в этом разделе.

Обучение нейронных сетей как минимизация функции ошибки.

Когда функционал ошибки нейронной сети задан (раздел 6.3), то главная

задача обучения нейронных сетей сводится к его минимизации. Градиентное

обучение – это итерационная процедура подбора весов, в которой каждый

следующий шаг направлен в сторону антиградиента функции ошибки.

Математически это можно выразить следующим образом:

[pic], или , что то же самое : [pic],

здесь (( - темп обучения на шаге (. В теории оптимизации этот метод

известен как метод наискорейшего спуска.[]

Метод обратного распространения ошибки.

Исторически наибольшую трудность на пути к эффективному правилу

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

функции ошибки [pic]. Дело в том, что ошибка сети определяется по ее

выходам, т.е. непосредственно связана лишь с выходным слоем весов. Вопрос

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

найти производные по соответствующим весам. Нужна была процедура передачи

ошибки с выходного слоя к предшествующим слоям сети, в направлении обратном

обработке входной информации. Поэтому такой метод, когда он был найден,

получил название метода обратного распространения ошибки (error back-

propagation ).

Разберем этот метод на примере двухслойного персептрона с одним

нейроном на выходе.(рис 6.1) Для этого воспользуемся введенными ранее

обозначениями. Итак,

[pic] -Функция ошибки (13)

[pic] -необходимая коррекция весов коррекция весов (14)

для выходного слоя (v записывается следующим образом.

[pic]

Коррекция весов между входным и скрытым слоями производится по формуле:

[pic] (15)

[pic] [pic] [pic]

Подставляя одно выражение в другое получаем

[pic] (16)

Производная функции активации, как было показано ранее (раздел 6.1),

вычисляется через значение самой функции. [pic][pic]

Непосредственно алгоритм обучения состоит из следующих шагов:

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

сети.

2. Вычислить выход сети y(x) по формуле (12).

3. Вычислить разность между выходом сети и требуемым значением для данного

вектора (13).

4. Если была допущена ошибка при классификации выбранного вектора, то

подкорректировать последовательно веса сети сначала между выходным и

скрытым слоями (15), затем между скрытым и входным (16).

5. Повторять шаги с 1 по 4 для каждого вектора обучающего множества до тех

пор, пока ошибка на всем множестве не достигнет приемлемого уровня.

Несмотря на универсальность, этот метод в ряде случаев становится

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

увеличить скорость сходимости функционала ошибки, разработано много

модификаций стандартного алгоритма, в частности две из которых и

предлагается использовать.

Многостраничное обучение.

С математической точки зрения обучение нейронных сетей (НС) – это

многопараметрическая задача нелинейной оптимизации. В классическом методе

обратного распространения ошибки (single-режим) обучение НС рассматривается

как набор однокритериальных задач оптимизации. Критерий для каждой задачи -

качество решения одного примера из обучающей выборки. На каждой итерации

алгоритма обратного распространения параметры НС (синаптические веса и

смещения) модифицируются так, чтобы улучшить решение одного примера. Таким

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

оптимизации.

Из теории оптимизации следует, что при решении многокритериальных задач

модификации параметров следует производить, используя сразу несколько

критериев (примеров), в идеале - все. Тем более нельзя ограничиваться одним

примером при оценке производимых изменений значений параметров.

Для учета нескольких критериев при модификации параметров используют

агрегированные или интегральные критерии, которые могут быть, например,

суммой, взвешенной суммой или квадратным корнем от суммы квадратов оценок

решения отдельных примеров.

В частности, в настоящих исследованиях изменения весов проводилось

после проверки всей обучающей выборки, при этом функция ошибки

рассчитывалась в виде :

[pic]

где,

k - номер обучающей пары в обучающей выборке, k=1,2,…,n1+n2

n1 - количество векторов первого класса;

n2 - число векторов второго класса.

Как показывают тестовые испытания, обучение при использовании пакетного

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

Автоматическая коррекция шага обучения.

В качестве еще одного расширения традиционного алгоритма обучения

предлагается использовать так называемый градиентный алгоритм с

автоматическим определением длины шага (. Для его описания необходимо

определить следующий набор параметров:

. начальное значение шага(0 ;

. количество итераций, через которое происходит запоминание данных сети

(синоптических весов и смещений);

. величина (в процентах) увеличения шага после запоминания данных сети, и

величина уменьшения шага в случае увеличения функции ошибки.

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

Затем происходит заданное число итераций обучения с заданным шагом. Если

после завершения этих итераций значение функции ошибки не возросло, то шаг

обучения увеличивается на заданную величину, а текущие значения весов и

смещений записываются на диск. Если на некоторой итерации произошло

увеличение функции ошибки, то с диска считываются последние запомненные

значения весов и смещений, а шаг обучения уменьшается на заданную величину.

При использовании автономного градиентного алгоритма происходит

автоматический подбор длины шага обучения в соответствии с характеристиками

адаптивного рельефа, и его применение позволило заметно сократить время

обучения сети без потери качества полученного результата.

Эффект переобучения.

Одна из наиболее серьезных трудностей изложенного подхода обучения

заключается в том, что таким образом минимизируется не та ошибка, которую

на самом деле нужно минимизировать, а ошибка, которую можно ожидать от

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

хотелось бы, чтобы нейронная сеть обладала способностью обобщать результат

на новые наблюдения. В действительности сеть обучается минимизировать

ошибку на обучающем множестве, и в отсутствие идеального и бесконечно

большого обучающего множества это совсем не то же самое, что минимизировать

"настоящую" ошибку на поверхности ошибок в заранее неизвестной модели

явления [5]. Иначе говоря, вместо того, чтобы обобщить известные примеры,

сеть запомнила их. Этот эффект и называется переобучением.

Соответственно возникает проблема – каким методом оценить ошибку

обобщения? Поскольку эта ошибка определена для данных, которые не входят в

обучающее множество, очевидным решением проблемы служит разделение всех

имеющихся в нашем распоряжении данных на два множества: обучающее – на

котором подбираются конкретные значения весов, и валидационного – на

котором оцениваются предсказательные способности сети. На самом деле,

должно быть еще и третье множество, которое вообще не влияет на обучение и

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

сети. Ошибки, полученные на обучающем, валидационном и тестовом множестве

соответственно называются ошибка обучения, валидационная ошибка и тестовая

ошибка.

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

подхода:

. Ранняя остановка обучения;

. Прореживание связей (метод от большого к малому);

. Поэтапное наращивание сети ( от малого к большому).

Самым простым является первый метод. Он предусматривает вычисление во время

обучения не только ошибки обучения, но и ошибки валидации, используя ее в

качестве контрольного параметра. В самом начале работы ошибка сети на

обучающем и контрольном множестве будет одинаковой. По мере того, как сеть

обучается, ошибка обучения, естественно, убывает, и, пока обучение

уменьшает действительную функцию ошибок, ошибка на контрольном множестве

также будет убывать. Если же контрольная ошибка перестала убывать или даже

стала расти, это указывает на то, что сеть начала слишком близко

аппроксимировать данные и обучение следует остановить. Рисунок 6.5 дает

качественное представление об этой методике.

Использование этой методики в работе с сейсмическими данными затруднено

тем обстоятельством, что исходная выборка очень мала, а хотелось бы как

можно больше данных использовать для обучения сети. В связи с этим было

принято решение отказаться от формирования валидационного множества, а в

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

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

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

этого условия проводились дополнительные эксперименты, показавшие что при

определенном минимуме ошибки обучения достигался относительный минимум

ошибки на тестовых данных.

Два других подхода для контроля переобучения предусматривают

постепенное изменение структуры сети. Только в одном случае происходит

эффективное вымывание малых весов (weight elimination) ,т.е. прореживание

малозначительных связей, а во втором, напротив, поэтапное наращивание

сложности сети. [3,4,5].

6.6 Формирование обучающей выборки и оценка эффективности обученной

нейросетевой модели.

Из исходных данных необходимо сформировать как минимум две выборки –

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

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

оценки эффективности обученной нейронной сети.

Как правило, используют следующую методику: из всей совокупности данных

случайным образом выбирают около 90% векторов для обучения, а на оставшихся

10% тестируют сеть. Однако, в условиях малого количества примеров эта

процедура становится неэффективной с точки зрения оценивания вероятности

ошибки классификации. В разделе 4.4 был описан другой, наиболее точный

метод расчета ошибки классификации. Это, так называемый, метод скользящего

экзамена (синонимы: cross-validation, “plug-in”-метод).[7,9].

В терминах нейронных сетей основную идею метода можно выразить так:

выбирается один вектор из всей совокупности данных, а остальные

используются для обучения НС. Далее, когда процесс обучения будет завершен,

предъявляется этот выбранный вектор и проверяется правильно сеть распознала

его или нет. После проверки выбранный вектор возвращается в исходную

выборку. Затем выбирается другой вектор, на оставшихся сеть вновь

обучается, и этот новый вектор тестируется. Так повторяется ровно n1+n2

раз, где n1–количество векторов первого класса, а n2 - второго.

По завершению алгоритма общая вероятность ошибки P подсчитывается

следующим образом:

[pic], где

N= n1+n2 - общее число примеров;

E– число ошибочных векторов (сеть неправильно распознала предъявляемый

пример).

Недостатком этого метода являются большие вычислительные затраты,

связанные с необходимость много раз проводить процедуру настройки весовых

коэффициентов, а в следствии этого и большое количество времени, требуемое

для вычисления величины [pic].

Однако в случае с малым количеством данных для определения

эффективности обученной нейронной сети рекомендуется применять именно метод

скользящего экзамена или некоторые его вариации. Стоит отметить, что

эффективность статистических методов классификации сейсмических сигналов

также проверяется методом скользящего экзамена. Таким образом, применяя его

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

экспериментов с возможностями стандартных методов.

7. Программная реализация.

При выборе программного обеспечения для решения определенной задачи с

помощью нейронных сетей возможны два подхода: использование готовых решений

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

программы.

Первый подход позволяет получить быстрое решение, не вдаваясь в

детальное изучение работы алгоритма. Однако, хорошие пакеты, в которых

имеются мощные средства для реализации различных парадигм нейронных сетей,

обработки данных и детального анализа результатов, стоят больших денег. И

это сильно ограничивает их применение. Еще одним недостатком является то,

что несмотря на свою универсальность, коммерческие пакеты не реализовывают

абсолютно все возможности моделирования и настройки нейронных сетей, и не

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

Если же возникает необходимость построить нейросетевое решение,

адаптированное под определенную задачу, и реализованное в виде отдельного

программного модуля, то первый подход чаще всего неприемлем.

Именно такие требования и были выдвинуты на начальном этапе

исследований. Точнее, необходимо было разработать программу,

предназначенную для классификации сейсмических данных при помощи

нейросетевых технологий, а также работающую под операционной системой Unix

(Linux и Sun Solaris SystemV 4.2). В результате была разработана программа,

реализующая основные идеи нейроинформатики, изложенные в разделе 6.

Следует отметить, что базовый алгоритм программы был выполнен под

системой Windows 95, а лишь затем оптимизирован под Unix по той причине,

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

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

программы требуется много времени.

Для большей совместимости версий под различные платформы использовались

возможности языка программирования С.

7.1 Функциональные возможности программы.

В программе “nvclass.с” – (нейро-классификатор векторов данных)

реализована модель двухслойного персептрона, представленная в разделе 6.

Эта программа предназначена для соотнесения тестируемого вектора признаков

сейсмической информации к одному из двух классов. Входные данные

представляют собой предысторию сейсмических явлений конкретного региона, а

также тестируемый вектор признаков, соответствующий сейсмическому событию,

не включенному в предысторию. Эти данные считываются из соответствующих

файлов в виде набора векторов признаков заданной размерности.

Автоматически, в зависимости от размерности входных векторов, определяется

конфигурация нейронной сети т.е. по умолчанию для заданной размерности

входных данных выбирается определенное (рекомендуемое по результатам

предварительных экспериментов) число нейронов в входном и скрытом слоях,

хотя при желании эти параметры легко меняются.

В качестве правила обучения нейронной сети реализован алгоритм

обратного распространения ошибки и некоторые методы его оптимизации. После

завершения процесса обучения тестируемый вектор признаков подается на вход

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

вероятностью соотносит входной вектор к одному из заданных классов.

Для дополнительной настройки нейронной сети в программе реализован ряд

процедур, описанных в разделе 6. Из них можно выделить следующие:

. Различные процедуры начальной инициализации весовых коэффициентов;

. Пакетный режим обучения;

. Алгоритм коррекции шага обучения;

. Процедуры предварительной обработки данных;

. Алгоритм оценки эффективности – cross-validation.

. Процедура многократного обучения.

Последняя процедура (многократное обучение сети) предусмотрена для

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

число раз (Cycle) генерируются матрицы начальных весовых коэффициентов и

выполняется алгоритм обучения и идентификации тестового вектора. По

полученным результатам обучения и тестирования выбирается вариант

наибольшего повторения и итоговое решение принимается исходя из него. Для

исключения неоднозначности это число выбирается положительным , целым,

нечетным.

В программе “nvclass” предусмотрены следующие режимы функционирования:

. «Внешний» режим – идентификации тестового входного вектора признаков;

. «Внутренний» режим – идентификация вектора признаков из набора векторов

предыстории.

«Внешний» режим предназначен для классификации вновь поступивших

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

. Классификация на нейронных сетях, уже обученных на данных из конкретных

регионов регионов,

. Классификация с повторным обучением нейронной сети.

«Внутренний» режим служит для оценки вероятности ошибки идентификации

сети и включает два подрежима – проверки правильности идентификации одного

из векторов набора предыстории и последовательной проверки всех заданных

векторов (“cross-validation”).

Режим работы программы устанавливается в файле настроек.

После завершения работы основные результаты записываются в

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

детального анализа. Пример файла приведен в приложении 4.

7.2 Общие сведения.

1. Программный пакет предназначенный для идентификации типа сейсмического

события включает следующие модули:

. Исходный код программы “nvclass.c” и “nvclass.h”;

. Файл с настройками режима работы программы “nvclass.inp”;

. Файл с обучающей выборкой векторов “vector.txt”;

. Файл с векторами для тестирования сети “vector.tst”;

. Файл, содержащий описание определенной конфигурации сети и весовые

коэффициенты этой уже обученной сети “nor18.net”.

. Файл автоматической компиляции “Makefile” (Только для версии под Unix).

. Файл отчета о результатах работы программы “Report.txt”.

2. В настоящий момент разработано две версии программы. Одна работает под

операционной системой Dos 6.2 и выше, а другая под Unix (Linux, Solaris

V4.2).

3. Необходимое средство компиляции:

. Для Dos (Windows) – любой компилятор Си. Например, Borland C++ 3.1 или

выше.

. Для Unix стандартный компилятор cc, входящий в состав базовой

комплектации любой операционной системы семейства Unix.

7.3 Описание входного файла с исходными данными.

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

файл, в котором хранится информация о размерности векторов, их количестве и

сами вектора данных. Файл должен иметь форму числовой матрицы. Каждая

строка матрицы соответствует одному вектору признаков. Количество признаков

должно совпадать с параметром NDATA. Количество столбцов равно количеству

признаков плюс два. Первый столбец содержит порядковый номер вектора в

общей совокупности данных (соответствует последовательности 1, 2,

3,...,NPATTERN), а в последнем столбце записаны значения указателя

классификатора: 1- для вектора из первого класса, 0 – для вектора из

второго класса. Все числовые параметры разделяются пробелами и записываются

в кодах ASCII. Пример файла приведен в приложении 2.

7.4 Описание файла настроек.

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

“nvclass.inp”. Пример файла приведен в приложении 3. Для настройки

используются следующие переменные:

TYPE - РЕЖИМ РАБОТЫ ПРОГРАММЫ

TYPE=1_1

Это значение соответствует внешнему режиму функционирования

программы без обучения нейронной сети, т.е. тестирование на заранее

обученной нейронной сети. При этом надо задать следующие параметры:

1. NDATA –Размерность входных данных

2. TESTVECTOR – Имя файла с тестируемым вектором

3. NETWORKFILE – Имя файла с матрицами весов предварительно обученной

сети

TYPE=1_2

Это значение соответствует внешнему режиму функционирования

программы с обучением нейронной сети и тестированием на ней заданного

вектора. Необходимо задать следующие параметры:

1. NDATA –Размерность входных данных

2. NPATTERN –Количество векторов признаков

3. PATTERNFILE-Имя файла с набором векторов признаков

4. TESTVECTOR – Имя файла с тестируемым вектором;

5. RESNETFNAME- Имя выходного файла с матрицами весов обученной сети.

TYPE=2_1

Данное значение соответствует внутреннему режиму с проверкой

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

программы необходимо задать следующие параметры:

1. NDATA –Размерность входных данных

2. NPATTERN –Количество векторов признаков

3. PATTERNFILE -Имя файла с набором векторов признаков

4. NUMBERVECTOR -Номер тестового вектора признаков из заданной выборки

TYPE=2_2

При данном значении параметра программа будет функционировать во

внутреннем режиме с последовательной проверкой всех векторов

(“cross_validation”). Необходимо задать следующие параметры :

1. NDATA -Размерность входных данных

2. NPATTERN –Количество векторов признаков

3. PATTERNFILE -Имя файла с набором векторов признаков

NDATA РАЗМЕРНОСТЬ ВЕКТОРОВ ПРИЗНАКОВ

Задается размерность векторов признаков, или количество признаков

в каждом векторе наблюдений. Этой величине должны соответствовать все

входные данные в текущем сеансе работы программы.

NPATTERN КОЛИЧЕСТВО ВЕКТОРОВ ПРИЗНАКОВ

Этот числовой параметр характеризует объем обучающей выборки и

соответствует количеству строк во входном файле PATTERNFILE.

PATTERNFILE ИМЯ ФАЙЛА С НАБОРОМ ВЕКТОРОВ ПРИЗНАКОВ

Имя файла, содержащего наборы векторов признаков предыстории

сейсмических явлений региона с указателями классификатора.

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

Имя файла, содержащего вектор признаков, который необходимо

идентифицировать. Файл должен иметь форму строки (числа разделяются

пробелами). Количество признаков должно соответствовать переменной

NDATA.

NETWORKFILE ИМЯ ФАЙЛА С МАТРИЦАМИ ВЕСОВ ПРЕДВАРИТЕЛЬНО ОБУЧЕННОЙ СЕТИ.

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

предварительно обученной нейронной сети с фиксированной размерностью

входных данных. Файл формируется на предыдущих этапах работы

программы. Необходимо учитывать количество признаков NDATA (явно

указанных в имени файла, под которые проектировалась нейронная сеть

(NDATA соответствует количеству входов сети) и символьную

аббревиатуру региона, из которого получена сейсмическая информация.

RESNETFNAME ИМЯ ВЫХОДНОГО ФАЙЛА С МАТРИЦАМИ ВЕСОВ ОБУЧЕННОЙ СЕТИ

Имя файла, содержащего параметры спроектированной и обученной

нейронной сети в данном сеансе эксплуатации программы. В имени файла

обязательно следует указывать символьную абревиатуру региона, из

которого получена сейсмическая информация и размерность векторов

признаков NDATA обрабатываемой информации, чтобы избежать путаницы в

интерпретации разных моделей. (Например, norv18.net или isrl9.net).

NUMBERVECTOR ПОРЯДКОВЫЙ НОМЕР ВЕКТОРА ПРИЗНАКОВ

Этот параметр соответствует номеру вектора признаков (номеру

строки в первом столбце матрицы) из файла PATTERNFILE. Этот вектор

признаков с указателем классификатора в дальнейшем будет

интерпретироваться как тестовый вектор. Он удаляется из всего набора

, а оставшиеся NPATTERN-1 векторов будут использованы в качестве

обучающей выборки.

REPORTFNAME ИМЯ ФАЙЛА ОТЧЕТА

Имя файла с результатами работы программы.

InitWeigthFunc ФУНКЦИЯ ИНИЦИАЛИЗАЦИИ НАЧАЛЬНЫХ ВЕСОВЫХ КОЭФФИЦИЕНТОВ

СЕТИ.

InitWeigthFunc=Gauss

Начальные матрицы весовых коэффициентов будут выбраны как

нормально распределенные случайные величины с математическим

ожиданием Alfa и среднеквадратическом отклонении Sigma (

N[Alfa,Sigma]).

InitWeigthFunc=Random

Начальные матрицы весовых коэффициентов будут выбраны как

равномерно распределенные случайные величины в диапазоне [-

Constant,Constant].

(Значение по умолчанию – InitWeigthFunc= RandomDistribution[-3,3],

т.е. Constant=3)

Constant ДИАПАЗОН РАВНОМЕРНО РАСПРЕДЕЛЕННЫХ СЛУЧАЙНЫХ ВЕЛИЧИН

Смотри InitWeigthFunc …

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

СЛУЧАЙНЫХ ВЕЛИЧИН

Смотри InitWeigthFunc …

Alfa МАТЕМАТИЧЕСКОЕ ОЖИДАНИЕ НОРМАЛЬНО РАСПРЕДЕЛЕННЫХ СЛУЧАЙНЫХ

ВЕЛИЧИН

Смотри InitWeigthFunc …

WidrowInit NGUYEN-WIDROW ИНИЦИАЛИЗАЦИЯ .

Параметр позволяет сформировать начальные весовые коэффициенты по

методике предложенной Nguyen и Widrow. Возможные варианты: “Yes” –

провести соответствующую инициализацию. “No”- не использовать эту

процедуру.(Значение по умолчанию – “No”)

Shuffle ПЕРЕМЕШИВАНИЕ ВЕКТОРОВ ПРИЗНАКОВ

При значении параметра “Yes” – входные вектора будут

предварительно перемешаны. При “No” – вектора будут подаваться на

вход сети в той последовательности, в которой они расположены во

входном файле (PATTERNFILE). (Значение по умолчанию – “Yes”).

Scaling ПРЕДВАРИТЕЛЬНАЯ ОБРАБОТКА ВЕКТОРОВ ПРИЗНАКОВ.

Этот параметр служит для использования в рамках программы

“nvclass” процедуры масштабирования входных данных. Эта процедура

позволяет значительно ускорить процесс обучения нейронной сети, а

также качественно улучшает результаты тестирования. Возможные

значения параметра: “Yes”,”No”. (Значение по умолчанию – “Yes”).

LearnToleranse ТОЧНОСТЬ ОБУЧЕНИЯ.

Параметр определяющий качество обучения нейронной сети. При

достижении заданной точности ? для каждого вектора признаков из

обучающей выборки настройка весовых коэффициентов сети заканчивается

и сеть считается обученной. (Значение по умолчанию – 0.1).

Eta КОЭФФИЦИЕНТ ОБУЧЕНИЯ НЕЙРОННОЙ СЕТИ.

Значение коэффициента задает скорость и качество обучения

нейронной сети. Используется для алгоритма обратного распространения

ошибки. (Значение по умолчанию–1.0)

MaxLearnCycles МАКСИМАЛЬНОЕ КОЛИЧЕСТВО ИТЕРАЦИЙ ОБУЧЕНИЯ

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

будет автоматически завершен. (Величина по умолчанию- 2000)

Loop КОЛИЧЕСТВО ПОВТОРОВ ОБУЧЕНИЯ.

Параметр задает величину полных циклов функционирования программы

(целое нечетное число). В каждом цикле формируются начальные матрицы

весов производится обучение сети и осуществляется классификация

тестового вектора. Результаты всех циклов обрабатываются, и

формируется итоговое заключение . (Значение по умолчанию=1).

7.5 Алгоритм работы программы.

Алгоритм работы программы зависит от режима, в котором она

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

операций:

1. Инициализация сети;

2. Настройка;

3. Проверка тестовых векторов.

Инициализация

В этом разделе происходит считывание всех данных из соответствующих

файлов (файл с примерами обучающей выборки, файл с конфигурацией обученной

сети, файл с примерами для тестирования). Затем, в зависимости от режима

функционирования, либо происходит инициализация всех весовых коэффициентов

сети заданным образом, либо сразу начинается проверка тестовых векторов на

обученной заранее нейронной сети, конфигурация которой считана из файла.

Настройка.

Если выбранный режим предусматривает выполнение алгоритма обучения

нейронной сети, то программа, после считывания исходных данных, и начальной

инициализации весовых коэффициентов выполняет процедуру их настройки до тех

пор, пока не выполнится одно из условий остановки. Либо значение ошибки

обучения достигнет желаемого уровня и сеть будет считаться обученной, либо

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

число. По мере выполнения алгоритма автоматически формируется полный отчет

о состоянии сети.

Проверка тестовых векторов.

На этом этапе происходит тестирование заданных векторов. Причем

возможны два варианта: тестируемый вектор может быть считан из файла, а

также можно задать номер тестируемого вектора в выборке исходных данных и

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

записываются в файл отчета.

7.6 Эксплуатация программного продукта.

Перед тем, как приступить к эксплуатации программного продукта

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

записаны исходная выборка векторов и с основными переменными файла настроек

программы.

Для корректной работы в дальнейшем желательно придерживаться

определенной последовательности действий:

1. Подготовить исходные данные согласно принятом формату.

2. Изменить в соответствии с требованиями определенные поля в файле

настроек.

3. Запустить программу.

4. Проанализировать результат, записанный в соответствующем файле.

7.7 Результат работы программы.

Для исследований возможностей разработанного программного обеспечения

были проведены различные эксперименты, основная цель которых - подобрать

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

работы содержали наименьшее количество ошибок идентификации. Методика, по

которой оценивалась ошибка классификации, основана на подходе “cross-

validation”.

Эксперименты проводились на данных, полученных из сейсмограмм,

записанных в Норвежской сейсмологической сети. В исходной выборке

насчитывалось 86 событий из разных классов, из них соответственно 50 –

землетрясений и 36 – взрывов. Исследования проводились для разного числа

признаков идентификации, а именно для 18 и 9 размерных векторов признаков.

Первая серия экспериментов была проведена на 18 размерных векторах.

Структура нейронной сети соответствовала , где 18 – количество

нейронов во входном слое, 9- число нейронов на первом скрытом слое , 1-

размерность выхода сети. Увеличение нейронов на скрытом слое не приводило к

улучшению результатов, а при уменьшении возникали дополнительные ошибки, в

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

Далее представлены описание параметров настройки программы во входных

файлах и результаты тестирования.

В качестве начальной конфигурации использовались следующие значения

настраиваемых параметров в файле “nvclass.inp”:

TYPE=2_2

NDATA=18

NPATTERN=86

PatternFile=norv18.pat

NetStructure=[18,9,1]

WidrowInit=No

Shuffle=Yes

Scaling=Yes

Eta=0.7

MaxLearnCycles=1950

Loop=5

Результаты экспериментов отражают количество ошибок идентификации от

различных параметров настройки программы.

Для примера рассмотрим влияние процедуры начальной инициализации

весовых коэффициентов и точности обучения на ошибку классификации. На

рисунках 7.1 и 7.2 едставлены эти результаты.

Отметим, что более стабильные результаты получаются в случае

инициализации весов при помощи нормально распределенных величин. Можно

добиться всего лишь 4-5 ошибок из 86, что соответствует ошибке

идентификации равной 5-6 процентов.

Для 9 размерных векторов признаков была использована следующая

структура нейронной сети , т.е. 5 нейронов на скрытом слое было

достаточно для получения хороших результатов.

В качестве примера приведем исследования аналогичные тем, которые

описаны выше.(Рис. 7.3, 7.4).

Последнюю диаграмму можно представить в виде.

Уже сейчас можно сделать вывод, что при использовании не всего набора

признаков идентификации, а некоторой части признаков результаты заметно

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

начальной инициализации не играет.

Представленные эксперименты не отражают полной картины о возможностях

применения нейронных сетей для идентификации типа сейсмического события, но

они экспериментально подтверждают эффективность нейросетевых технологий для

решения этой задачи.

8. Заключение

Проведенные исследования подтвердили эффективность применения

нейросетевых технологий для идентификации типа источника сейсмических

события. При определенных настройках нейронной сети можно добиться

результатов, когда вероятность правильного распознавания составляет 96.5%.

Ошибки возникают только на 3 векторах из 86. Если сравнивать полученные

результаты с теми, которые можно достичь при использовании стандартных

методов классификации, один из вариантов которых приведен в разделе 4, то

они практически повторяют друг друга. И статистика и нейронные сети

ошибаются одинаковое количество раз, причем на одних и тех же векторах. Из

86 событий статистические методы ошибаются на 3 векторах (1–землетрясение и

2-взрыва), и нейросетевой классификатор также ошибается именно на этих

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

одного метода над другим.

Заметим, что в настоящих исследованиях были использованы довольно общие

и универсальные технологии нейроинформатики (многослойные сети применяются

для решения многих задач, но это не всегда самая оптимальная

нейроструктура), а применение более узких и специализированных нейронных

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

частности, при помощи нейропакетов на тех же данных были поставлены ряд

экспериментов над сетями Кохонена, описанными в разделе 5.4. Результаты

показали, что количество ошибок идентификации в большинстве случаев

составляет 3-4 вектора, т.е. практически совпадают с результатами,

полученными на многослойных сетях и классических методах.

Итак, подводя итог всему выше сказанному, выделим основные результаты

проведенных исследований:

1. Нейронные сети позволяют успешно решать проблему определения типа

источника сейсмического события.

2. Новое решение не уступает по эффективности традиционным методам,

использующимся в настоящее время для решения исследуемой задачи.

3. Возможны улучшения технических характеристик нейросетевого решения.

В качестве дальнейших исследований, направленных на повышение

эффективности нейросетевого решения, можно предложить следующие:

. Для многослойных сетей прямого распространения решить проблему начальной

инициализации весовых коэффициентов. Если предположить, что существует

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

конкретной реализацией выборки данных, предназначенной для обучения

нейронной сети, то можно объяснить те случаи, когда результат несколько

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

учитывающего распределение исходных данных позволит получать более

стабильные результаты.

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

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

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

предназначенной конкретно для решения данной задачи могут привести к

улучшению полученных результатов.

Список литературы.

1. Уоссермен Ф. “Нейрокомпьютерная техника” - М.: Мир,1992.

2. Горбань А.Н., Дубинин-БарковскийВ.Л., Кирдин А.Н. “Нейроинформатика” СП

“Наука” РАН 1998.

3. Горбань А.Н., Россиев Д.А. “Нейронные сети на персональном компьютере”

СП “Наука” РАН 1996.

4. Ежов А.А., Шумский С.А. “Нейрокомпьютинг и его применение в экономике и

бизнесе”.1998.

5. Bishop C.M. “Neural Networks and Pattern Recognition.” Oxford Press.

1995.

6. Goldberg D. “Genetic Algorithms in Machine Learning, Optimization, and

Search.” – Addison-Wesley,1988.

7. Fausett L.V. “Fundamentals of Neural Networks: Architectures, Algorithms

and Applications”, Prentice Hall, 1994.

8. Kohonen T. “Self-organization and Associative Memory”, Berlin: Springer-

Verlag, 1989.

9. Kushnir A.F., Haikin L.M., Troitsky E.V. “Physics of the earth and

planetary interiors” 1998.

10. Копосов А.И., Щербаков И.Б., Кисленко Н.А., Кисленко О.П., Варивода

Ю.В. Отчет по научно-исследовательской работе "Создание аналитического

обзора информационных источников по применению нейронных сетей для задач

газовой технологии"; и др., ВНИИГАЗ, 1995, www.neuralbench.ru

11. Fukunaga K., Kessel D.L., “Estimation of classification error”, IEEE

Trans. Comp. C 20, 136-143. 1971.

12. Деев А.Д., “Применение статистического дискриминационного анализа и его

ассимптотического расширения для сравнения различных размерностей

пространства.”, РАН 195, 759-762. 1970.

Приложение.

1. Пример выборки сейсмограмм.

В левом столбце представлены сейсмограммы, описывающие взрывы, а в

правом – землетрясения.

2. Пример файла с векторами признаков.

Представлена выборка из файла 9_Norv.txt, содержащего 9 размерные вектора

признаков.

NumOfPattern: 86

PatternDimens: 9

1 -14.3104 -13.2561 -13.4705 -13.4306 -14.1015 -13.3503 -13.3805

-13.7369 -0.3494 0

2 -14.6881 -13.6349 -12.9050 -13.4323 -14.2279 -13.4720 -13.2117

-13.5791 -1.2801 0

3 -14.4036 -14.1745 -13.8014 -12.7209 -14.6283 -13.9589 -13.4649

-12.9716 -0.8250 0

…

55 -14.3693 -13.4362 -11.4072 -12.3129 -14.8612 -13.3480 -12.8517

-13.4014 -0.7738 0

56 -14.2856 -12.6858 -13.8215 -13.4282 -14.0982 -13.1587 -13.2792

-13.7852 -1.3442 0

57 -14.4822 -13.1141 -13.7787 -13.4466 -13.6761 -13.2969 -13.6033

-13.9252 -0.6642 1

58 -13.5522 -13.1302 -13.5444 -14.1471 -13.2994 -13.2368 -13.9776

-14.4295 -0.9973 1

59 -14.8524 -11.9846 -13.7231 -14.2496 -13.4809 -13.0515 -13.8950

-14.3923 -1.8284 1

…

85 -14.5994 -13.6920 -12.8539 -13.7629 -14.1699 -13.2075 -13.3422

-13.6788 -11.9537 1

86 -14.3821 -13.6093 -12.8677 -13.7788 -14.1260 -13.3246 -13.2966

-13.6453 -11.4304 1

3. Файл с настройками программы

# Common parameters for programm "NVCLASS"

# # # # # # # # # # # # # # # # # # # #

# 1_1 - OnlyTest mode , 1_2 - TestAfterLearn mode,

# 2_1 - CheckOneVector , 2_2 - CrossValidation mode.

TYPE=2_2

NDATA=9

NPATTERN=86

PatternFile=9_Norv.txt

NTEST=10

TestVector=vector.tst

NetworkFile=9.net

ResNetFname=9.net

NumberVector=57

ReportFile=Report.txt

Debug=Yes

# Next parameters was define in result experiments and if you will

# change it, the any characteristics of Neural Net may be not optimal

# (since may be better then optimal).

# # # # # # # # # # #

# 'NetStructure' must be: [NDATA,NUNIT1,1] (NOUT=1 always)

# value 'AUTO'-'NetStructure' will be define the programm.(See help).

# example : [18,9,1], or [18,18,1], or [9,9,5,1]

NetStructure=[18,12,1]

# may be: [Gauss] or [Random]

InitWeigthFunc=Gauss

Constant=3

Alfa=0

Sigma=1.5

Widrow=No

Shuffle=Yes

Scaling=Yes

LearnTolerance=0.1

Eta=1

MaxLearnCycles=50

Loop=3

#end of list

4. Пример файла отчета.

NVCLASS report - Wed Jun 02 15:58:02 1999

Type = 1_2

Neural Net -

PatternFile - vect.txt

Test Vector(s) - vector.tst

ResNetFname - 12.net

LearnTolerance = 0.10

InitialWeigthFunc = Gauss[ 0.0, 1.5]

< loop="" 1="">

Learning cycle result:

NumIter = 5

NumLE = 3

Error vector(s): 58, 59, 63,

+-----+------+--------+------+

| N | ID | Result |Target|

+-----+------+--------+------+

| 1 | 24 | 0.1064 | 0 |

| 2 | 25 | 0.9158 | 1 |

| 3 | 26 | 0.0452 | 0 |

| 4 | 27 | 0.0602 | 0 |

| 5 | 28 | 0.0348 | 0 |

| 6 | 29 | 0.0844 | 0 |

| 7 | 30 | 0.1091 | 0 |

| 8 | 31 | 0.0821 | 0 |

| 9 | 32 | 0.0298 | 0 |

| 10 | 33 | 0.2210 | 0 |

+-----+------+--------+------+

< loop="" 2="">

Learning cycle result:

NumIter = 5

NumLE = 5

Error vector(s): 33, 34, 55, 58, 63,

+-----+------+--------+------+

| N | ID | Result |Target|

+-----+------+--------+------+

| 1 | 24 | 0.1279 | 0 |

| 2 | 25 | 0.9929 | 1 |

| 3 | 26 | 0.0960 | 0 |

| 4 | 27 | 0.1463 | 0 |

| 5 | 28 | 0.1238 | 0 |

| 6 | 29 | 0.1320 | 0 |

| 7 | 30 | 0.1478 | 0 |

| 8 | 31 | 0.1235 | 0 |

| 9 | 32 | 0.0740 | 0 |

| 10 | 33 | 0.5140 | 1 |

+-----+------+--------+------+

5. Файл описания функций, типов переменных и используемых библиотек

“nvclass.h”.

* --- Neuro classificator---

* Common defines

#include

//#include

#include

#define DefName "nvclass.inp"

#define MAXDEF 100

#define MAXLINE 256

#define NMAXPAT 100

#define NMXINP 20

#define NMXUNIT 20

#define CONT 0

#define EXIT_OK 1

#define EXIT_CNT 2

#define RESTART 911

#define MAXEXP 700 /* Max arg exp(arg) without error 'OVERFLOW'

#define Random 10

#define Gauss 20

#define OK 0

#define Error 1

#define Yes 77

#define No 78

#define Min 0 /* Find_MinMax(...) */

#define Max 1

#define TYPE_ONE 21

#define TYPE_TWO 22

#define TYPE_THREE 23

#define TYPE_FOUR 24

int NDATA = 0;

int NUNIT1 = 0;

int NUNIT2 = 0;

int NUNIT3 = 0;

int NOUT = 1;

int NPATTERN = 0; /* Number of input pattern*/

int NWORK = 0; /* Number of work pattern*/

int NTEST= 0; /* Number of test pattern*/

int result;

int STOP = 0;

int NumOut = 250; /* Number of itteration, after which show result in

debugfile. */

int Num_Iter=10;/* The parameters requred in the procecc of */

float Percent=0.25; /* dinamic lerning with change 'eta' */

float LearnTolerance = 0.10;

float TestTolerance = 0.5;

float MAX_ERR=0.00001; /* min error */

float eta = 1.0; /* learning coefficient*/

float MIN_ETA=0.000001;

float **Array_MinMax;

int *Cur_Number;

float W1[NMXINP][NMXUNIT];

float W2[NMXUNIT];

float PromW1[NMXINP][NMXUNIT];

float PromW2[NMXUNIT];

float PromW1_OLD[NMXINP][NMXUNIT];

float PromW2_OLD[NMXUNIT];

float Err1[NMXUNIT];

float Err2;

float OLD_ERROR;

float GL_Error=0.0;

float Out1[NMXUNIT];

float Out2;

char NetStr[20]="Auto"; /* String with pattern of Net Structure*/

int Type = TYPE_THREE; /* Enter the mode of work of programm */

int InitFunc = Random; /* Random [=10] weigth will RandomDistribution

Gauss [=20] - ... GaussianDistributon */

float Constant = 1; /* RandomDistribution [-Constant,Constant]*/

float Alfa = 0; /* GaussianDistribution [Alfa,Sigma]*/

float Sigma = 1; /* ... */

int Widrow = No; /* Nguyen-Widrow initialization start weigth*/

int Loop = 1; /* Number repeat of Learning cycle */

char *PatternFile; /* File with input patterns*/

char *TestVector;

char *ReportFile="report.txt"; /* name of report file */

char *NetworkFile; /* Name of input NetConfig file */

char *ResNetFname; /* Name of output NetConfig file */

int DEBUG = Yes; /* if 'Yes' then debug info in the DebugFile

char *DebugFile="Logfile.log"; /* Name of the debug file*/

int NumberVector = 0; /* Number of TEST vector */

int Shuffle = Yes; /* Flag - shuffle the input vectors*/

int Scaling = Yes; /* Scaling input vector */

int MaxLearnCycles = 1999; /* Max number of learning iteration */

FILE *Dfp; /* Debug file pointer */

FILE *Rfp; /* Report file pointer*/

typedef struct Pattern {

int ID; /* ID number this vector in all set of pattern */

float *A; /* pattern (vector) A={a[0],a[1],...,a[NDATA]} */

float Target; /* class which this vector is present*/

} PAT;

PAT *Input;

PAT *Work;

PAT *Test;

/* lines in defaults file are in the form "NAME=value" */

typedef struct Default {

char *name; /* name of the default */

char *value; /* value of the default */

} DEF;

/* structure of statistics info about one test vector */

typedef struct Statistic {

int ID; /* Primery number from input file */

float Target;

float TotalRes; /* Total propability */

int Flag; /* Flag = 1, if vector was error and = 0

in over case */

float *result; /* Result of testing vector on current

iteration */

int *TmpFlag; /* analog 'Flag' on current itteration */

int *NumIter; /* Number iteration of learning on which

Learning cycle STOPED */

int **NumLE; /* Error vectors after cycle of learning

was test*/

} STAT;

/* structure of the some result of learning cycle */

typedef struct ResLearning {

int NumIter;

int LearnError[NMAXPAT+1]; /* A[0]-count of error,

A[1]-ID1,

A[2]-ID2,...

A[NMAXRL]-ID?.*/

} RL;

/* function prototypes */

void OnlyTestVector(void);

void TestAfterLearn (void);

void CheckOneVector ( void );

void CrossValidation ( void );

DEF **defbuild(char *filename);

DEF *defread(FILE *fp);

FILE *defopen (char *filename);

char *defvalue(DEF **deflist, const char *name);

int defclose(FILE *fp);

void defdestroy(DEF **, int);

void getvalues(void);

void Debug (char *fmt, ...);

void Report (char *fmt, ...);

void Widrow_Init(void);

int Init_W( void );

float RavnRaspr(float A, float B);

float NormRaspr(float B,float A);

void ShufflePat(int *INP,int Koll_El);

float F_Act(float x);

float Forward (PAT src);

int LearnFunc (void);

int Reset (float ResErr, int Cnt, int N_Err);

void Update_Last (int n, float Total_Out);

void Update_Prom1 (int n);

void Prom_to_W (void);

void Update_All_W (int num, float err_cur );

void Init_PromW(void);

void Prom_to_OLD(void);

int CheckVector(float Res, PAT src);

int *TestLearn(int *src);

RL FurtherLearning(int NumIteration,

float StartLearnTolerans,

float EndLearnTolerans,

RL src);

STAT *definestat (PAT src);

STAT **DefineAllStat (PAT *src,int Num);

void FillStatForm (STAT *st, int iteration, float res, RL lr);

void FillSimpleStatForm (STAT *st, float res);

void destroystat ( STAT *st, int param);

void DestroyAllStat (STAT **st, int Num);

void PrintStatHeader(void);

void printstat(STAT *st);

void PrintStatLearn(RL src);

void PrintTestStat(STAT **st, int len);

void PrintErrorStat (STAT **st,int Len);

int DefineNetStructure (char *ptr);

void getStructure(char buf[20]);

PAT patcpy (PAT dest, PAT src);

PAT* LocPatMemory(int num);

void ReadPattern (PAT *input, char *name,int Len);

void FreePatMemory(PAT* src, int num);

void ShowPattern (char *fname, PAT *src, int len);

void ShowVector(char *fname,PAT src);

float getPatTarget (float res);

PAT* DataOrder (PAT* src,int Len, int Ubit, PAT* dest, PAT* test);

void FindMinMax (PAT *src,int Dimens, int Num_elem, float **Out_Array);

void ConvX_AB_01(PAT src);

int *DefineCN (int len);

int getPosition (int Num, int *src, int Len);

void DestroyCN (int *src);

void ShowCurN (int LEN);

float **LocateMemAMM(void);

void FreeAMM (float **src);

void WriteHeaderNet(char *fname, float **src);

void WriteNet (char *fname,int It);

void ReadHeaderNet(char *fname, float **src);

int ReadNet (char *fname, int It);

FILE *OpenFile(char *name);

int CloseFile(FILE *fp);

/* End of common file */

6. Файл автоматической компиляции программы под Unix -“Makefile”.

CC= cc

LIBS= -lm

OBJ= nvclass.o

nvclass: $(OBJ)

$(CC) -o nvclass $(LIBS) $(OBJ)

nvclass.o: nvclass.c

7. Основной модуль - “nvclass.с”

* Neuron Classificator ver 1.0

#include "common.h"

/* =========================

* MAIN MODULE

* =========================

void main (int argc, char *argv[])

{ int i;

char buf[MAXLINE], PrName[20], *ptr;

time_t tim;

time(&tim);

/* UNIX Module */

Dfp = OpenFile(DebugFile);

strcpy(buf,argv[0]);

ptr = strrchr(buf,'/');

ptr++;

strcpy(PrName,ptr);

Debug ("\n\n'%s' - Started %s",PrName,ctime(&tim));

getvalues();

Rfp = OpenFile(ReportFile);

DefineNetStructure(NetStr); /* NetStr string from input file */

getStructure(buf);

Debug ("\nNeyral net %s",buf);

Input = LocPatMemory(NPATTERN);

Work = LocPatMemory(NPATTERN);

Array_MinMax = LocateMemAMM();

Cur_Number = DefineCN (NPATTERN);

printf("\nMetka - 1");

if (Type == TYPE_ONE)

OnlyTestVector ();

if (Type == TYPE_TWO)

TestAfterLearn ();

if (Type == TYPE_THREE)

CheckOneVector ();

if (Type == TYPE_FOUR)

CrossValidation();

time(&tim);

Debug ("\n\n%s - Normal Stoped %s",PrName,ctime(&tim));

CloseFile(Dfp);

CloseFile(Rfp);

FreeAMM (Array_MinMax);

DestroyCN (Cur_Number);

FreePatMemory(Input,NPATTERN);

FreePatMemory(Work, NPATTERN);

}

* ^OnlyTestVectors - read net from (NetworkFile) and test the

TestVector(s)

void OnlyTestVector(void)

{ char buf[MAXLINE+1];

STAT **st, *stat;

int i,j;

float Res;

Debug ("\nOnlyTestVector proc start");

Debug ("\n NPATTERN = %d",NPATTERN);

Debug ("\n NTEST = %d",NTEST);

Test = LocPatMemory(NTEST);

ReadPattern(Test,TestVector, NTEST);

/* ShowPattern ("1.tst",Test,NTEST);*/

PrintStatHeader();

st = DefineAllStat (Test,NTEST);

ReadHeaderNet(NetworkFile,Array_MinMax);

if (Scaling == Yes)

{ for (i=0;i ",i+1);

ReadNet(NetworkFile,i+1);

for (j=0;j x/32768

* -> x from [0,1]

float RavnRaspr(float A, float B)

{float x;

x = (B-A)*rand()/(RAND_MAX+1.0) + A;

return x;

}

float NormRaspr(float A,float B)

{ float mat_ogidanie=A, Sigma=B;

float Sumx=0.0, x;

int i;

for (i=0;i N[a,sigma]*/

x = Sigma*(Sumx-6) + mat_ogidanie;

return x;

}

int Init_W ( void )

{ int i,j;

float A, B;

time_t t,t1;

t = time(NULL);

t1=t;

/* restart random generator*/

while (t==t1)

srand((unsigned) time(&t));

if (InitFunc == Random)

{ A = -Constant;

B = Constant;

Debug ("\nInit_W () --- Start (%ld))",t);

Debug ("\n InitFunc=Random[%4.2f,%4.2f]",A,B);

for(i=0; i LearnTolerance)

{ NumErr++;

Init_PromW();

Update_Last(K, Res);

Update_Prom1(K);

Prom_to_W();

}

err_cur = err_cur + (ep[n]*ep[n]);

}

err_cur=0.5*(err_cur/NWORK);

result = Reset(err_cur, num, NumErr);

if ((num % NumOut)==0)

Debug("\nStep :%d NumErr :%d Error:%6.4f",num,NumErr,err_cur);

} while (result == CONT || result == RESTART);

Debug("\nStep :%d NumErr :%d Error:%6.4f",num,NumErr,err_cur);

return num;

}

-----------------------

w11

8A. 1.3. 2CEA;>9=K9 ?5@A5?B@>=.

8A. 1.2. 5@540B>G=K5 DC=:F88 =59@>=>2.

(

y=f(s)

8A. 1.1. 0B5=

f(x)

w12

w3p

w40

vРис. 1.3. Двухслойный персептрон.

Рис. 1.2. Передаточные функции нейронов.

(

y=f(s)

Рис. 1.1. Математический нейрон

f(x)

w12

w3p

w40

Рис 1.4. Этапы нейросетевого решения задачи.

Рис 1.5. Процесс обучения нейронной сети.

Рис 1.6. Задача прогнозирования.

Рис 1.7. Задача распознавания.

Рис 1.8. Задача распознавания.

[pic]

Рис 1.9. Применение нейросетей для отыскания зависимостей переменных

Рис 1.10. Задача кластеризации.

Сейсмограмма

Решающее

правило

Землетрясение

Взрыв

Рис. 2.1. Постановка задачи классификации сейсмических сигналов.

Рис.3.1. Пример сейсмограммы.

Рис. 4.1. Одиночный нейрон

y=f(x)

(

[pic]

Рис. 4.2. Линейный дискриминатор.

Рис 4.3. Сети Ворда

Рис.4.4 Сеть Кохонена.

[pic]

Рис.4.5 Сеть Кохонена для классификации сейсмических сигналов.

[pic]

Рис 5.1. Гистограмма значений переменной при наличии редких, но больших по

амплитуде отклонений от среднего.

[pic]

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

активации [pic]

[pic]

Рис 5.3. Автоассоциативная сеть с узким горлом – аналог правила обучения

Ойя.

[pic]

Рис 5.4. Слой линейных нейронов

Рис. 6.1. Двухслойный персептрон.

(

z0(1

x0(1

w10

wH0

wHp

f(x)

Рис. 6.2. Единичная сигмоида, со значениями в диапазоне (0,1).

Ошибка обучения

Ошибка на тестовом наборе

Рис 6.3 Ранняя остановка обучения в момент минимума ошибки валидации.

[pic]

Рис. 7.1 Инициализация весовых коэффициентов случайным образом из интервала

[-a,a].

[pic]

Рис. 7.2. Инициализация весовых коэффициентов с помощью нормального

распределения с параметрами N[-(,(].

[pic]

Рис. 7.3. Инициализация весовых коэффициентов случайным образом из

интервала [-a,a].

[pic]

Рис 7.4. Инициализация весовых коэффициентов с помощью нормального

распределения с параметрами N[-(,(].

[pic]

Страницы: 1, 2