我如何在Keras中创建一个自定义损失函数，包括神经网络输出相对于输入的导数？

Question

我想训练一个神经网络来满足常微分方程。我用keras序列模型构建了一个这样的神经网络

import keras
from keras.models import Sequential
from keras.layers import Dense
from keras.layers.advanced_activations import LeakyReLU, PReLU
model = Sequential()
model.add(Dense(output_dim = 20, init = 'uniform', activation = 'relu', input_dim = 1,use_bias=True))
model.add(Dense(output_dim = 20, init = 'uniform', activation = 'relu',use_bias=True))
model.add(Dense(output_dim = 20, init = 'uniform', activation = 'relu',use_bias=True))
model.add(Dense(output_dim = 1,init='uniform',activation='sigmoid'))`

现在，我希望神经网络的输出满足以下微分方程：

dy/dx+y=x

其中x位于0和1之间。 假设N代表构造的神经网络的输出，我想要下面的损失函数

∑(dN/dx_i+N-x_i)^2

我如何使用Keras序列模型来实现这一点