Another way to add smoothness constraint is to add -norm of the derivative to the objective:
Note that the norm is sensitive to all the derivatives, not just the largest.
We can formulate an LP problem by adding a vector of optimization parameters which bound derivatives:
In matrix form,
The objective function becomes
See Fig.3.41 and Fig.3.42 for example results.
L-Infinity Norm of Derivative Objective