我目前正在Keras中试验生成对抗网络。 如this论文中所建议,我想使用历史平均损失函数。这意味着我要惩罚网络权重的变化。 我不确定如何巧妙地实现它。
我正在根据this帖子的答案实施自定义损失功能。
def historical_averaging_wrapper(current_weights, prev_weights):
def historical_averaging(y_true, y_pred):
diff = 0
for i in range(len(current_weights)):
diff += abs(np.sum(current_weights[i]) + np.sum(prev_weights[i]))
return K.binary_crossentropy(y_true, y_pred) + diff
return historical_averaging
网络的权重受到惩罚,权重在每批数据之后都会发生变化。
我的第一个想法是在每次批处理之后更新损失函数。 大概是这样的:
prev_weights = model.get_weights()
for i in range(len(data)/batch_len):
current_weights = model.get_weights()
model.compile(loss=historical_averaging_wrapper(current_weights, prev_weights), optimizer='adam')
model.fit(training_data[i*batch_size:(i+1)*batch_size], training_labels[i*batch_size:(i+1)*batch_size], epochs=1, batch_size=batch_size)
prev_weights = current_weights
这合理吗?我认为这种方法似乎有点“混乱”。 还有另一种以“更智能”的方式执行此操作的可能性吗? 就像更新数据生成器中的损失函数并使用fit_generator()一样? 预先感谢。
答案 0 :(得分:0)
损失函数是使用张量在图形上进行的操作。 您可以在损失函数中定义其他张量以保存先前的值。这是一个示例:
import tensorflow as tf
import tensorflow.keras.backend as K
keras = tf.keras
class HistoricalAvgLoss(object):
def __init__(self, model):
# create tensors (initialized to zero) to hold the previous value of the
# weights
self.prev_weights = []
for w in model.get_weights():
self.prev_weights.append(K.variable(np.zeros(w.shape)))
def loss(self, y_true, y_pred):
err = keras.losses.mean_squared_error(y_true, y_pred)
werr = [K.mean(K.abs(c - p)) for c, p in zip(model.get_weights(), self.prev_weights)]
self.prev_weights = K.in_train_phase(
[K.update(p, c) for c, p in zip(model.get_weights(), self.prev_weights)],
self.prev_weights
)
return K.in_train_phase(err + K.sum(werr), err)
变量prev_weights
保存先前的值。请注意,在计算重量误差之后,我们添加了K.update
操作。
用于测试的样本模型:
model = keras.models.Sequential([
keras.layers.Input(shape=(4,)),
keras.layers.Dense(8),
keras.layers.Dense(4),
keras.layers.Dense(1),
])
loss_obj = HistoricalAvgLoss(model)
model.compile('adam', loss_obj.loss)
model.summary()
一些测试数据和目标函数:
import numpy as np
def test_fn(x):
return x[0]*x[1] + 2.0 * x[1]**2 + x[2]/x[3] + 3.0 * x[3]
X = np.random.rand(1000, 4)
y = np.apply_along_axis(test_fn, 1, X)
hist = model.fit(X, y, validation_split=0.25, epochs=10)
在我的测试中,模型损失随时间减少。