linear controller with anti reset windup implemented by bounding the integral state:
[sim, u] = ld_lin_awup_controller(sim, ev, r, y, Ta, tfR, min__, max__)
e = r-y
u = ld_limited_integrator( e, min__, max__ ) + tfR*e