* @param Omega double specifying the value discrete-time frequency [0 pi] at which the group delay is evaluated.
* @return double containing the resulting group delay in samples.
*/
public double discreteTimeGroupDelay( double Omega ) {
Complex c = Complex.exp( new Complex(0.0, -Omega ) );
Complex N = new Complex( a[order]*order );
for ( int i = order-1; i >= 0; i-- ) {
N = N.times(c).plus( a[i]*i );
}
Complex D = evaluate( c );
return ( N.over(D) ).real();
}