numObservations = 128;
X = rand(numInputChannels,numObservations);
X = dlarray(X,"CB");
T = rand(numOutputChannels,numObservations);
T = dlarray(T,"CB");
function [loss,gradients,lap] = modelLoss(net,X,T)
Y = forward(net,X);
loss = l1loss(Y,T);
X = stripdims(X);
Y = stripdims(Y);
lap = dllaplacian(Y,X,1);
gradients = dlgradient(loss,net.Learnables);
end