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

Я пытаюсь разработать 1D сверточную нейронную сеть с остаточными связями и нормализацией партии на основе определения уровня аритмии на уровне кардиолога с сверточными нейронными сетями с использованием keras. Это код до сих пор:

# define model
x = Input(shape=(time_steps, n_features))

# First Conv / BN / ReLU layer
y = Conv1D(filters=n_filters, kernel_size=n_kernel, strides=n_strides, padding='same')(x) 
y = BatchNormalization()(y)
y = ReLU()(y)

shortcut = MaxPooling1D(pool_size = n_pool)(y)

# First Residual block
y = Conv1D(filters=n_filters, kernel_size=n_kernel, strides=n_strides, padding='same')(y) 
y = BatchNormalization()(y)
y = ReLU()(y)
y = Dropout(rate=drop_rate)(y)
y = Conv1D(filters=n_filters, kernel_size=n_kernel, strides=n_strides, padding='same')(y) 
# Add Residual (shortcut)
y = add([shortcut, y])

# Repeated Residual blocks   
for k in range (2,3): # smaller network for testing

    shortcut = MaxPooling1D(pool_size = n_pool)(y)
    y = BatchNormalization()(y)
    y = ReLU()(y)
    y = Dropout(rate=drop_rate)(y)
    y = Conv1D(filters=n_filters * k, kernel_size=n_kernel, strides=n_strides, padding='same')(y)    
    y = BatchNormalization()(y)
    y = ReLU()(y)
    y = Dropout(rate=drop_rate)(y)
    y = Conv1D(filters=n_filters * k, kernel_size=n_kernel, strides=n_strides, padding='same')(y) 
    y = add([shortcut, y])

z = BatchNormalization()(y)
z = ReLU()(z)    
z = Flatten()(z)
z = Dense(64, activation='relu')(z)
predictions = Dense(classes, activation='softmax')(z)

model = Model(inputs=x, outputs=predictions)

# Compiling 
model.compile(loss='categorical_crossentropy', optimizer='adam', metrics=['categorical_accuracy'])

# Fitting 
model.fit(train_x, train_y, epochs=n_epochs, batch_size=n_batch)

И это график упрощенной модели того, что я пытаюсь построить.

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

Сеть состоит из 16 остаточных блоков с 2 сверточными слоями на блок. Все сверточные слои имеют длину фильтра 16 и имеют 64k-фильтры, где k начинается с 1 и увеличивается каждый 4-й остаточный блок. Каждый альтернативный остаточный блок подбирает свои входы в 2 раза, поэтому исходный вход в конечном итоге подвыбор в 2 раза. Когда остаточный блок подвыражает ввод, соответствующие соединения быстрого доступа также подвыражают их ввод, используя операцию Max Pooling с тем же коэффициентом подвыборки.

Но я могу заставить его работать, если я использую одинаковое количество фильтров в каждом слое Conv1D, с k = 1, strides = 1 и padding = same, без применения MaxPooling1D. Любые изменения этих параметров приводят к несоответствию размера тензора и невозможности компиляции со следующей ошибкой:

ValueError: Operands could not be broadcast together with shapes (70, 64) (70, 128)

Кто-нибудь знает, как устранить несоответствие этого размера и заставить его работать?

Кроме того, если вход имеет более одного канала (или функций), несоответствие даже хуже! Есть ли способ иметь дело с несколькими каналами?