oustafod

PURPOSE

OUSTAFOD Oustaloup's approximation of a s^r operator

SYNOPSIS

function G=oustafod(r,N,wb,wh)

DESCRIPTION

 OUSTAFOD Oustaloup's approximation of a s^r operator

  OUSTAFOD(r,N,wb,wh) computes the Oustaloup filter approximation of a
           fractional-order operator s^r of the order N and valid in the
           frequency range of (wb, wh). The function returns a ZPK object
           containing the continuous-time filter. The following equation
           is used to construct the filter:
                       N
                     -----
           s^r = K *  | | (s+w'k)/(s+wk),
                      | |
                     k= -N

           where w'k = wb*(wh/wb)^((k+N+(1/2)*(1-r))/(2N+1))
                 wk  = wb*(wh/wb)^((k+N+(1/2)*(1+r))/(2N+1))
                 K   = wh^r.

CROSS-REFERENCE INFORMATION

This function calls:

tf TF Return Z(s) of the fractance circuit in form of a transfer function
zpk ZPK Return Z(s) of the fractance circuit in form of a zero-pole-gain model

This function is called by:

oustapp OUTSAPP Obtain a integer-order approximation of a fractional-order system.
frac_imp_rc_lowpass_oc1 FRAC_IMP_RC_LOWPASS_OC1 Generate component values for Cauer I RC LP imp.