Discretization schemes for fractional-order differentiators and integrators

For fractional-order differentiator s^r where r is a real number, its discretization is a key step in digital implementation. Two discretization methods are presented. The first scheme is a direct recursive discretization of the Tustin operator. The second one is a direct discretization method using the Al-Alaoui operator via continued fraction expansion (CFE). The approximate discretization is minimum phase and stable. Detailed discretization procedures and short MATLAB scripts are given. Examples are included for illustration