В данной статье мы рассмотрим, что такое рекуррентные нейронные сети и как создать нейронную сеть с нуля в Python.
Содержание
Зачем нужны рекуррентные нейронные сети
Создание рекуррентной нейронной сети на примере
Поставление задачи для рекуррентной нейронной сети
Составление плана для нейронной сети
Предварительная обработка рекуррентной нейронной сети RNN
Фаза прямого распространения нейронной сети
Фаза обратного распространения нейронной сети
Параметры рассматриваемой нейронной сети
Тестирование рекуррентной нейронной сети
Рекуррентные нейронные сети (RNN) — это тип нейронных сетей, которые специализируются на обработке последовательностей. Зачастую их используют в таких задачах, как обработка естественного языка (Natural Language Processing) из-за их эффективности в анализе текста. В данной статье мы наглядно рассмотрим рекуррентные нейронные сети, поймем принцип их работы, а также создадим одну сеть в Python, используя numpy.
Данная статья подразумевает наличие у читателя базовых знаний о нейронных сетях. Будет не лишним прочитать от том как создать нейронную сеть в Python, в которой показаны простые примеры использования нейронов в Python.
Приступим!
Зачем нужны рекуррентные нейронные сети
Один из нюансов работы с нейронными сетями (а также CNN) заключается в том, что они работают с предварительно заданными параметрами. Они принимают входные данные с фиксированными размерами и выводят результат, который также является фиксированным. Плюс рекуррентных нейронных сетей, или RNN, в том, что они обеспечивают последовательности с вариативными длинами как для входа, так и для вывода. Вот несколько примеров того, как может выглядеть рекуррентная нейронная сеть:
Входные данные отмечены красным, нейронная сеть RNN — зеленым, а вывод — синим.
Способность обрабатывать последовательности делает рекуррентные нейронные сети RNN весьма полезными. Области использования:
Машинный перевод (пример Google Translate) выполняется при помощи нейронных сетей с принципом «многие ко многим». Оригинальная последовательность текста подается в рекуррентную нейронную сеть, которая затем создает переведенный текст в качестве результата вывода;
Анализ настроений часто выполняется при помощи рекуррентных нейронных сетей с принципом «многие к одному». Этот отзыв положительный или отрицательный? Такая постановка является одним из примеров анализа настроений. Анализируемый текст подается нейронную сеть, которая затем создает единственную классификацию вывода. Например — Этот отзыв положительный.
Мы собрали ТОП Книг для Python программиста которые помогут быстро изучить язык программирования Python. Список книг: Книги по Python
Далее в статье будет показан пример создания рекуррентной нейронной сети по схеме «многие к одному» для анализа настроений.
Создание рекуррентной нейронной сети на примере
Представим, что у нас есть нейронная сеть, которая работает по принципу «многое ко многим«. Входные данные — x0, х1, … xn, а результаты вывода — y0, y1, … yn. Данные xi и yi являются векторами и могут быть произвольных размеров.
Рекуррентные нейронные сети RNN работают путем итерированного обновления скрытого состояния h, которое является вектором, что также может иметь произвольный размер. Стоит учитывать, что на любом заданном этапе t:
Следующее скрытое состояние ht подсчитывается при помощи предыдущего ht - 1 и следующим вводом xt;
Следующий вывод yt подсчитывается при помощи ht.
Рекуррентная нейронная сеть RNN многие ко многим
Вот что делает нейронную сеть рекуррентной: на каждом шаге она использует один и тот же вес. Говоря точнее, типичная классическая рекуррентная нейронная сеть использует только три набора параметров веса для выполнения требуемых подсчетов:
Wxh используется для всех связок xt ? ht
Whh используется для всех связок ht-1 ? ht
Why используется для всех связок ht ? yt
Для рекуррентной нейронной сети мы также используем два смещения:
bh добавляется при подсчете ht
by добавляется при подсчете yt
Вес будет представлен как матрица, а смещение как вектор. В данном случае рекуррентная нейронная сеть состоит их трех параметров веса и двух смещений.
Следующие уравнения являются компактным представлением всего вышесказанного:
Разбор уравнений лучше не пропускать. Остановитесь на минутку и изучите их внимательно. Помните, что вес — это матрица, а другие переменные являются векторами.
Говоря о весе, мы используем матричное умножение, после чего векторы вносятся в конечный результат. Затем применяется гиперболическая функция в качестве функции активации первого уравнения. Стоит иметь в виду, что другие методы активации, например, сигмоиду, также можно использовать.
Не знаете, что такое функция активации? Вы можете ознакомиться с ними в вводной статье о нейронных сетях. Для оптимальной работы это важно.
Поставление задачи для рекуррентной нейронной сети
К текущему моменту мы смогли реализовать рекуррентную нейронную сеть RNN с нуля. Она должна выполнить простой анализ настроения. В дальнейшем примере мы попросим сеть определить, будет заданная строка нести позитивный или негативный характер.
Вот несколько примеров из небольшого набора данных, который был собран для данной статьи:
Текст
Позитивный?
Я хороший
Да
Я плохой
Нет
Это очень хорошо
Да
Это неплохо
Да
Я плохой, а не хороший
Нет
Я несчастен
Нет
Это было хорошо
Да
Я чувствую себя неплохо, мне не грустно
Да
Составление плана для нейронной сети
В следующем примере будет использована классификация рекуррентной сети «многие к одному». Принцип ее использования напоминает работу схемы «многие ко многим», что была описана ранее. Однако на этот раз будет задействовано только скрытое состояние для одного пункта вывода y:
Рекуррентная нейронная сеть RNN многие к одному
Каждый xi будет вектором, представляющим определенное слово из текста. Вывод y будет вектором, содержащим два числа. Одно представляет позитивное настроение, а второе — негативное. Мы используем функцию Softmax, чтобы превратить эти значения в вероятности, и в конечном счете выберем между позитивным и негативным.
Приступим к созданию нашей рекуррентной нейронной сети.
Предварительная обработка рекуррентной нейронной сети RNN
Упомянутый ранее набор данных состоит из двух словарей Python:
Python
1
2
3
4
5
6
7
8
9
10
11
train_data={
'good':True,
'bad':False,
# ... больше данных
}
test_data={
'this is happy':True,
'i am good':True,
# ... больше данных
}
True = Позитивное, False = Негативное
Для получения данных в удобном формате потребуется сделать определенную предварительную обработку. Для начала необходимо создать словарь в Python из всех слов, которые употребляются в наборе данных:
print('%d unique words found'%vocab_size)# найдено 18 уникальных слов
vocab теперь содержит список всех слов, которые употребляются как минимум в одном учебном тексте. Далее присвоим каждому слову из vocab индекс типа integer (целое число).
Python
1
2
3
4
5
6
# Назначить индекс каждому слову
word_to_idx={w:ifori,winenumerate(vocab)}
idx_to_word={i:wfori,winenumerate(vocab)}
print(word_to_idx['good'])# 16 (это может измениться)
print(idx_to_word[0])# грустно (это может измениться)
Теперь можно отобразить любое заданное слово при помощи индекса целого числа. Это очень важный пункт, так как:
Рекуррентная нейронная сеть не различает слов — только числа.
Напоследок напомним, что каждый ввод xi для рассматриваемой рекуррентной нейронной сети является вектором. Мы будем использовать веторы, которые представлены в виде унитарного кода. Единица в каждом векторе будет находиться в соответствующем целочисленном индексе слова.
Так как в словаре 18 уникальных слов, каждый xi будет 18-мерным унитарным вектором.
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
importnumpy asnp
defcreateInputs(text):
'''
Возвращает массив унитарных векторов
которые представляют слова в введенной строке текста
- текст является строкой string
- унитарный вектор имеет форму (vocab_size, 1)
'''
inputs=[]
forwintext.split(' '):
v=np.zeros((vocab_size,1))
v[word_to_idx[w]]=1
inputs.append(v)
returninputs
Мы используем createInputs() позже для создания входных данных в виде векторов и последующей их передачи в рекуррентную нейронную сеть RNN.
Фаза прямого распространения нейронной сети
Пришло время для создания рекуррентной нейронной сети. Начнем инициализацию с тремя параметрами веса и двумя смещениями.
Обратите внимание: для того, чтобы убрать внутреннюю вариативность весов, мы делим на 1000. Это не самый лучший способ инициализации весов, но он довольно простой, подойдет для новичков и неплохо работает для данного примера.
Для инициализации веса из стандартного нормального распределения мы используем np.random.randn().
Затем мы реализуем прямую передачу рассматриваемой нейронной сети. Помните первые два уравнения, рассматриваемые ранее?
Эти же уравнения, реализованные в коде:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
classRNN:
# ...
defforward(self,inputs):
'''
Выполнение передачи нейронной сети при помощи входных данных
Возвращение результатов вывода и скрытого состояния
Вывод - это массив одного унитарного вектора с формой (input_size, 1)
'''
h=np.zeros((self.Whh.shape[0],1))
# Выполнение каждого шага в нейронной сети RNN
fori,xinenumerate(inputs):
h=np.tanh(self.Wxh@x+self.Whh@h+self.bh)
# Compute the output
y=self.Why@h+self.by
returny,h
Довольно просто, не так ли? Обратите внимание на то, что мы инициализировали h для нулевого вектора в первом шаге, так как у нас нет предыдущего h, который теперь можно использовать.
Давайте попробуем следующее:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
# ...
defsoftmax(xs):
# Применение функции Softmax для входного массива
returnnp.exp(xs)/sum(np.exp(xs))
# Инициализация нашей рекуррентной нейронной сети RNN
rnn=RNN(vocab_size,2)
inputs=createInputs('i am very good')
out,h=rnn.forward(inputs)
probs=softmax(out)
print(probs)# [[0.50000095], [0.49999905]]
Наша рекуррентная нейронная сеть работает, однако ее с трудом можно назвать полезной. Давайте исправим этот недочет.
Фаза обратного распространения нейронной сети
Для тренировки рекуррентной нейронной сети будет использована функция потери. Здесь будет использована потеря перекрестной энтропии, которая в большинстве случаев совместима с функцией Softmax. Формула для подсчета:
Здесь pc является предсказуемой вероятностью рекуррентной нейронной сети для класса correct (позитивный или негативный). Например, если позитивный текст предсказывается рекуррентной нейронной сетью как позитивный текст на 90%, то потеря составит:
При наличии параметров потери можно натренировать нейронную сеть таким образом, чтобы она использовала градиентный спуск для минимизации потерь. Следовательно, здесь понадобятся градиенты.
Обратите внимание: следующий раздел подразумевает наличие у читателя базовых знаний об многовариантном исчислении. Вы можете пропустить несколько абзацев, однако мы рекомендуем все пробежаться по ним глазами. По мере получения новых данных код будет дополняться, и объяснения станут понятнее.
Оригиналы всех кодов, которые использованы в данной инструкции, доступны на GitHub.
Готовы? Продолжим!
Параметры рассматриваемой нейронной сети
Параметры данных, которые будут использованы в дальнейшем:
y — необработанные входные данные нейронной сети;
р — конечная вероятность: р = softmax(y);
с — истинная метка определенного образца текста, так называемый «правильный» класс;
L — потеря перекрестной энтропии: L = -ln(pc);
Wxh, Whh и Why — три матрицы веса в рассматриваемой нейронной сети;
bh и by — два вектора смещения в рассматриваемой рекуррентной нейронной сети RNN.
Установка
Следующим шагом будет настройка фазы прямого распространения. Это необходимо для кеширования отдельных данных, которые будут использоваться в фазе обратного распространения нейронной сети. Параллельно с этим можно будет установить основной скелет для фазы обратного распространения. Это будет выглядеть следующим образом:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
classRNN:
# ...
defforward(self,inputs):
'''
Выполнение фазы прямого распространения нейронной сети с
использованием введенных данных.
Возврат итоговой выдачи и скрытого состояния.
- Входные данные в массиве однозначного вектора с формой (input_size, 1).
'''
h=np.zeros((self.Whh.shape[0],1))
self.last_inputs=inputs
self.last_hs={0:h}
# Выполнение каждого шага нейронной сети RNN
fori,xinenumerate(inputs):
h=np.tanh(self.Wxh@x+self.Whh@h+self.bh)
self.last_hs[i+1]=h
# Подсчет вывода
y=self.Why@h+self.by
returny,h
defbackprop(self,d_y,learn_rate=2e-2):
'''
Выполнение фазы обратного распространения нейронной сети RNN.
- d_y (dL/dy) имеет форму (output_size, 1).
- learn_rate является вещественным числом float.
'''
pass
Градиенты
Настало время математики! Начнем с вычисления . Что нам известно:
Здесь используется фактическое значение , а также применяется дифференцирование сложной функции. Результат следующий:
К примеру, если p = [0.2, 0.2, 0.6], а корректным классом является с = 0, то конечным результатом будет значение = [-0.8, 0.2, 0.6]. Данное выражение несложно перевести в код:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
# Цикл для каждого примера тренировки
forx,yintrain_data.items():
inputs=createInputs(x)
target=int(y)
# Прямое распространение
out,_=rnn.forward(inputs)
probs=softmax(out)
# Создание dL/dy
d_L_d_y=probs
d_L_d_y[target]-=1
# Обратное распространение
rnn.backprop(d_L_d_y)
Отлично. Теперь разберемся с градиентами для Why и by, которые используются только для перехода конечного скрытого состояния в результат вывода рассматриваемой нейронной сети RNN. Используем следующие данные:
Здесь hn является конечным скрытым состоянием. Таким образом:
Аналогичным способом вычисляем:
Теперь можно приступить к реализации backprop().?
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
classRNN:
# ...
defbackprop(self,d_y,learn_rate=2e-2):
'''
Выполнение фазы обратного распространения нейронной сети RNN.
- d_y (dL/dy) имеет форму (output_size, 1).
- learn_rate является вещественным числом float.
'''
n=len(self.last_inputs)
# Подсчет dL/dWhy и dL/dby.
d_Why=d_y@self.last_hs[n].T
d_by=d_y
Напоминание: мы создали self.last_hs в forward() в предыдущих примерах.
Наконец, нам понадобятся градиенты для Whh, Wxh, и bh, которые использовались в каждом шаге нейронной сети. У нас есть:
Изменение Wxh влияет не только на каждый ht, но и на все у , что, в свою очередь, приводит к изменениям в L. Для того, чтобы полностью подсчитать градиент Wxh, необходимо провести обратное распространение через все временные шаги. Его также называют Обратным распространением во времени, или Backpropagation Through Time (BPTT):
Обратное распространение во времени
Wxh используется для всех прямых ссылок xt ? ht, поэтому нам нужно провести обратное распространение назад к каждой из этих ссылок.
Приблизившись к заданному шагу t, потребуется подсчитать :
Производная гиперболической функции tanh нам уже известна:
Используем дифференцирование сложной функции, или цепное правило:
Аналогичным способом вычисляем:
Последнее нужное значение — . Его можно подсчитать рекурсивно:
Реализуем обратное распространение во времени, или BPTT, отталкиваясь от скрытого состояния в качестве начальной точки. Далее будем работать в обратном порядке. Поэтому на момент подсчета значение будет известно. Исключением станет только последнее скрытое состояние hn:
Теперь у нас есть все необходимое, чтобы наконец реализовать обратное распространение во времени ВРТТ и закончить backprop():
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
classRNN:
# ...
defbackprop(self,d_y,learn_rate=2e-2):
'''
Выполнение фазы обратного распространения RNN.
- d_y (dL/dy) имеет форму (output_size, 1).
- learn_rate является вещественным числом float.
'''
n=len(self.last_inputs)
# Вычисление dL/dWhy и dL/dby.
d_Why=d_y@self.last_hs[n].T
d_by=d_y
# Инициализация dL/dWhh, dL/dWxh, и dL/dbh к нулю.
d_Whh=np.zeros(self.Whh.shape)
d_Wxh=np.zeros(self.Wxh.shape)
d_bh=np.zeros(self.bh.shape)
# Вычисление dL/dh для последнего h.
d_h=self.Why.T@d_y
# Обратное распространение во времени.
fortinreversed(range(n)):
# Среднее значение: dL/dh * (1 - h^2)
temp=((1-self.last_hs[t+1]**2)*d_h)
# dL/db = dL/dh * (1 - h^2)
d_bh+=temp
# dL/dWhh = dL/dh * (1 - h^2) * h_{t-1}
d_Whh+=temp@self.last_hs[t].T
# dL/dWxh = dL/dh * (1 - h^2) * x
d_Wxh+=temp@self.last_inputs[t].T
# Далее dL/dh = dL/dh * (1 - h^2) * Whh
d_h=self.Whh@temp
# Отсекаем, чтобы предотвратить разрыв градиентов.
fordin[d_Wxh,d_Whh,d_Why,d_bh,d_by]:
np.clip(d,-1,1,out=d)
# Обновляем вес и смещение с использованием градиентного спуска.
self.Whh-=learn_rate*d_Whh
self.Wxh-=learn_rate*d_Wxh
self.Why-=learn_rate*d_Why
self.bh-=learn_rate*d_bh
self.by-=learn_rate*d_by
Моменты, на которые стоит обратить внимание:
Мы объединили в для удобства;
Мы постоянно обновляем переменную d_h, которая держит самую последнюю версию , что требуется для подсчета ;
Закончив с обратным распространением во времени ВРТТ, мы используем np.clip() на значениях градиента ниже -1 или выше 1. Это поможет избавиться от проблемы со взрывными градиентами. Такое случается, когда градиенты становятся слишком большими из-за огромного количества умноженных параметров. Взрыв, а также исчезновение градиентов не считается редкостью для классических рекуррентных нейронных сетей. Более сложные рекуррентные нейронные сети, например LSTM, лучше подойдут для их обработки.
Когда все градиенты подсчитаны, мы обновляем параметры веса и смещения, используя градиентный спуск.
Мы сделали это! Наша рекуррентная нейронная сеть готова.
Тестирование рекуррентной нейронной сети
Наконец настал тот момент, которого мы так долго ждали — протестируем готовую рекуррентную нейронную сеть.
Для начала, напишем вспомогательную функцию для обработки данных рассматриваемой рекуррентной нейронной сети:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
importrandom
defprocessData(data,backprop=True):
'''
Возврат потери рекуррентной нейронной сети и точности для данных
- данные представлены как словарь, что отображает текст как True или False.
- backprop определяет, нужно ли использовать обратное распределение
'''
items=list(data.items())
random.shuffle(items)
loss=0
num_correct=0
forx,yinitems:
inputs=createInputs(x)
target=int(y)
# Прямое распределение
out,_=rnn.forward(inputs)
probs=softmax(out)
# Вычисление потери / точности
loss-=np.log(probs[target])
num_correct+=int(np.argmax(probs)==target)
ifbackprop:
# Создание dL/dy
d_L_d_y=probs
d_L_d_y[target]-=1
# Обратное распределение
rnn.backprop(d_L_d_y)
returnloss/len(data),num_correct/len(data)
Теперь можно написать цикл для тренировки сети:
Python
1
2
3
4
5
6
7
8
9
10
# Цикл тренировки
forepoch inrange(1000):
train_loss,train_acc=processData(train_data)
ifepoch%100==99:
print('--- Epoch %d'%(epoch+1))
print('Train: Loss %.3f | Accuracy: %.3f'%(train_loss,train_acc))
Вот и все, пошаговое руководство по рекуррентным нейронным сетям на этом закончено. Мы узнали, что такое RNN, как они работают, почему они полезны, как их создавать и тренировать. Это очень малый аспект мира нейронных сетей. При желании вы можете продолжить изучение темы самостоятельно, используя следующие ресурсы:
Подробнее ознакомьтесь с LTSM. Это долгая краткосрочная память, которая характерна более мощной архитектурой рекуррентных нейронных сетей. Будет не лишним ознакомиться с управляемым рекуррентными блоками, или GRU. Это наиболее популярная вариация LTSM;
Поэкспериментируйте с более крупными и сложными RNN. Для этого используйте подходящие ML библиотеки, например, Tensorflow, Keras или PyTorch;
Прочтите о двунаправленных нейронных сетях, которые обрабатывают последовательности как в прямом, так и в обратном направлении. Это позволяет получить больше информации на уровне вывода;
Познакомьтесь поближе с Natural Language Toolkit (NLTK), популярной библиотекой Python для работы с данными на языках, которые используют люди, а не машины.