Возведение в квадрат с помощью TensorFlow с классом набора данных

Я хочу написать нейронную сеть, которая ищет распределение x ^ 2 без предопределенной модели. Точнее, ему даются некоторые точки в [-1,1] с их квадратами для обучения, а затем ему нужно будет воспроизвести и предсказать подобное, например. [-10,10]. У меня более-менее получилось - без наборов данных. Но затем я попытался изменить его, чтобы использовать наборы данных и научиться их использовать. Теперь мне удалось запустить программу, но вывод стал хуже, чем раньше, в основном это постоянный 0.

Предыдущая версия была похожа на x^2 в [-1,1] с линейным продолжением, что было лучше. >Предыдущий вывод с синей линией, теперь плоской. И цель была бы совпасть с красным..

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

# square2.py - drugie podejscie do trenowania sieci za pomocą Tensorflow
# cel: nauczyć sieć rozpoznawać rozkład x**2
# analiza skryptu z:
# https://stackoverflow.com/questions/43140591/neural-network-to-predict-nth-square

import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np
from tensorflow.python.framework.ops import reset_default_graph

# def. danych do trenowania sieci
# x_train = (np.random.rand(10**3)*4-2).reshape(-1,1)
# y_train = x_train**2
square2_dane = np.load("square2_dane.npz")
x_train = square2_dane['x_tren'].reshape(-1,1)
y_train = square2_dane['y_tren'].reshape(-1,1) 

# zoptymalizować dzielenie danych
# x_train = square2_dane['x_tren'].reshape(-1,1)
# ds_x = tf.data.Dataset.from_tensor_slices(x_train)
# batch_x = ds_x.batch(rozm_paczki)
# iterator = ds_x.make_one_shot_iterator()

# określenie parametrów sieci
wymiary = [50,50,50,1]
epoki = 500
rozm_paczki = 200

reset_default_graph()
X = tf.placeholder(tf.float32, shape=[None,1])
Y = tf.placeholder(tf.float32, shape=[None,1])

weights = []
biases = []
n_inputs = 1

# inicjalizacja zmiennych
for i,n_outputs in enumerate(wymiary):
    with tf.variable_scope("layer_{}".format(i)):
        w = tf.get_variable(name="W", shape=[n_inputs,n_outputs],initializer = tf.random_normal_initializer(mean=0.0,stddev=0.02,seed=42))
        b=tf.get_variable(name="b",shape=[n_outputs],initializer=tf.zeros_initializer)
        weights.append(w)
        biases.append(b)
        n_inputs=n_outputs

def forward_pass(X,weights,biases):
    h=X
    for i in range(len(weights)):
        h=tf.add(tf.matmul(h,weights[i]),biases[i])
        h=tf.nn.relu(h)
    return h    

output_layer = forward_pass(X,weights,biases)
f_strat = tf.reduce_mean(tf.squared_difference(output_layer,Y),1)
f_strat = tf.reduce_sum(f_strat)
# alternatywna funkcja straty
#f_strat2 = tf.reduce_sum(tf.abs(Y-y_train)/y_train)
optimizer = tf.train.AdamOptimizer(learning_rate=0.003).minimize(f_strat)

with tf.Session() as sess:
    sess.run(tf.global_variables_initializer())
    # trenowanie
    dataset = tf.data.Dataset.from_tensor_slices((x_train,y_train))
    dataset = dataset.batch(rozm_paczki)
    dataset = dataset.repeat(epoki)
    iterator = dataset.make_one_shot_iterator()
    ds_x, ds_y = iterator.get_next()
    sess.run(optimizer, {X: sess.run(ds_x), Y: sess.run(ds_y)})
    saver = tf.train.Saver()
    save = saver.save(sess, "./model.ckpt")
    print("Model zapisano jako: %s" % save)

    # puszczenie sieci na danych
    x_test = np.linspace(-1,1,600)
    network_outputs = sess.run(output_layer,feed_dict = {X :x_test.reshape(-1,1)})

plt.plot(x_test,x_test**2,color='r',label='y=x^2')
plt.plot(x_test,network_outputs,color='b',label='sieć NN')
plt.legend(loc='right')
plt.show()

Я думаю, что проблема с вводом обучающих данных sess.run(optimizer, {X: sess.run(ds_x), Y: sess.run(ds_y)}) или с определением ds_x, ds_y. Это моя первая такая программа. Так что это был вывод для строк (внутри блока «видит»)

with tf.Session() as sess:
    sess.run(tf.global_variables_initializer())
    # trenowanie
    for i in range(epoki):
        idx = np.arange(len(x_train))
        np.random.shuffle(idx)
        for j in range(len(x_train)//rozm_paczki):
            cur_idx = idx[rozm_paczki*j:(rozm_paczki+1)*j]
            sess.run(optimizer,feed_dict = {X:x_train[cur_idx],Y:y_train[cur_idx]})
    saver = tf.train.Saver()
    save = saver.save(sess, "./model.ckpt")
    print("Model zapisano jako: %s" % save)

Спасибо!

PS: меня очень вдохновила нейронная сеть для прогнозирования n-го квадрата< /а>


person whereabouts    schedule 04.12.2017    source источник


Ответы (1)


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

sess.run(optimizer, {X: sess.run(ds_x), Y: sess.run(ds_y)})
  1. Будет выполнен только один шаг обучения, потому что этот код не находится в цикле. Ваш исходный код выполнил len(x_train)//rozm_paczki шагов, что должно привести к большему прогрессу.

  2. Два вызова sess.run(ds_x) и sess.run(ds_y) выполняются на отдельных этапах, что означает, что они будут содержать значения из разных пакетов, которые не связаны между собой. Каждый вызов sess.run(ds_x) или sess.run(ds_y) перемещает Iterator к следующему пакету и отбрасывает любые части элемента ввода, которые вы явно не запросили при вызове sess.run(). По сути, вы получите X из пакета i и Y из пакета i+1 (или наоборот), а модель будет обучаться на недопустимых данных. Если вы хотите получить значения из одного и того же пакета, вам нужно сделать это за один вызов sess.run([ds_x, ds_y]).

Есть еще две проблемы, которые могут повлиять на эффективность:

  1. Dataset не перемешивается. Ваш исходный код вызывает np.random.shuffle() в начале каждой эпохи. Вы должны включить dataset = dataset.shuffle(len(x_train)) перед dataset = dataset.repeat().

  2. Неэффективно извлекать значения из Iterator обратно в Python (например, когда вы делаете sess.run(ds_x)) и возвращать их обратно на шаг обучения. Более эффективно передавать выходные данные операции Iterator.get_next() непосредственно на шаг прямой связи в качестве входных данных.

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

import tensorflow as tf
import matplotlib.pyplot as plt
import numpy as np

# Generate training data.
x_train = np.random.rand(10**3, 1).astype(np.float32) * 4 - 2
y_train = x_train ** 2

# Define hyperparameters.
DIMENSIONS = [50,50,50,1]
NUM_EPOCHS = 500
BATCH_SIZE = 200

dataset = tf.data.Dataset.from_tensor_slices((x_train,y_train))
dataset = dataset.shuffle(len(x_train))  # (Point 3.) Shuffle each epoch.
dataset = dataset.repeat(NUM_EPOCHS)
dataset = dataset.batch(BATCH_SIZE)
iterator = dataset.make_one_shot_iterator()

# (Point 2.) Ensure that `X` and `Y` correspond to the same batch of data.
# (Point 4.) Pass the tensors returned from `iterator.get_next()`
# directly as the input of the network.
X, Y = iterator.get_next()

# Initialize variables.
weights = []
biases = []
n_inputs = 1
for i, n_outputs in enumerate(DIMENSIONS):
  with tf.variable_scope("layer_{}".format(i)):
    w = tf.get_variable(name="W", shape=[n_inputs, n_outputs],
                        initializer=tf.random_normal_initializer(
                            mean=0.0, stddev=0.02, seed=42))
    b = tf.get_variable(name="b", shape=[n_outputs],
                        initializer=tf.zeros_initializer)
    weights.append(w)
    biases.append(b)
    n_inputs = n_outputs

def forward_pass(X,weights,biases):
  h = X
  for i in range(len(weights)):
    h=tf.add(tf.matmul(h, weights[i]), biases[i])
    h=tf.nn.relu(h)
  return h

output_layer = forward_pass(X, weights, biases)
loss = tf.reduce_sum(tf.reduce_mean(
    tf.squared_difference(output_layer, Y), 1))
optimizer = tf.train.AdamOptimizer(learning_rate=0.003).minimize(loss)
saver = tf.train.Saver()

with tf.Session() as sess:
  sess.run(tf.global_variables_initializer())

  # (Point 1.) Run the `optimizer` in a loop. Use try-while-except to iterate
  # until all elements in `dataset` have been consumed.
  try:
    while True:
      sess.run(optimizer)
  except tf.errors.OutOfRangeError:
    pass

  save = saver.save(sess, "./model.ckpt")
  print("Model saved to path: %s" % save)

  # Evaluate network.
  x_test = np.linspace(-1, 1, 600)
  network_outputs = sess.run(output_layer, feed_dict={X: x_test.reshape(-1, 1)})

plt.plot(x_test,x_test**2,color='r',label='y=x^2')
plt.plot(x_test,network_outputs,color='b',label='NN prediction')
plt.legend(loc='right')
plt.show()
person mrry    schedule 05.12.2017