Derive the discrete transfer function for the parallel form of a PID controller.
G_{c}(s)= K_{c}(1+ \frac {1}{τ_{1}s}+ τ_{D}s) (1)
using the BD substitution for s. Compare the result with the velocity form of the PID algorithm given in Equation
\Delta p_{k} = p_{k}-p_{k-1}=K_{c}[(e_{k}-e_{k-1})+\frac{\Delta t}{\tau _{I}}e_{k} +\frac{\tau _{D}}{\Delta t}(e_{k}-2e_{k-1}+e_{k-2}) ] (*)