Wide Band Integrators – What Are They?

For some time now it has been conventional ‘wisdom’ that using time based digital integration may cause amplitude errors in the result and that these get worse as the frequency increases. As a result of this, integration using Omega arithmetic has been prevalent by using Fourier Transforms of the signal. This, of course, remains a valid approach and is particularly useful if the data is already in the frequency domain, which was its prime purpose.

It has also been conventional wisdom that a higher order integrator is ‘better’. Now, a rectangular integrator is a zero order integrator and a trapezium integrator is a first order integrator. So according to the ‘second wisdom’ the trapezoidal integrator should produce better results than the rectangular integrator. This is not the case. As will become evident, further studies into wide band integrators have shown this amplitude degredation is limited to what are refered to as ‘odd order’ integrators, in particular integrators such as the Trapeziodal integrator. Even order integrators, such as the rectangle integrator, do not have this frequency related degredation. What remains valid is that using a higher even order integrator is more accurate than that of a lower even order integrator.

First let us illustrate the problem using a trapezoidal integrator. Figure 1 below is the result of using a trapezium integrator on a 100Hz sinewave of amplitude 2\pi and sampled at 1024 samples/second which is well below half the Nyquist frequency. The second graph (Figure 2) is the same except at 300Hz. In both case the blue curve is the mathematically exact result and the black curve is the trapezoidal integration. In both cases the Trapezoidal integration has a loss of amplitude with this loss being worse as the frequency increases towards the Nyquist frequency (Sample Rate/2).

Trapezium integrator on a 100Hz sinewave
Figure 1: Trapezium integrator on a 100Hz sinewave
Trapezium integrator on a 300Hz sinewave
Figure 2: Trapezium integrator on a 300Hz sinewave

Recently, a new set of wide band even order integrators up to 6th order has been introduced into Prosig’s DATS software. These are accurate over the entire frequency range and are comparable to the omega arithmetic algorithm in terms of accuracy. The examples in Figure 3 below show the same two cases as above except that a 4th order wide band integrator has been used in place of the trapezium integrator. There is excellent agreement. The slight apparent discrepancies in the higher frequency case are in fact due to an effective delay as discussed later.

Wide Band integrator on 100Hz & 300Hz sinewaves
Figure 3: Wide Band integrator on 100Hz & 300Hz sinewaves

Ideal Integrator

As a start consider the frequency response characteristics of an ideal integrator. It is readily shown that the desired frequency response of an ideal integrator is

H(i{\omega}) = \frac{1}{i{\omega}} = -\frac{i}{\omega} = \frac{1}{\omega}e^{-\frac{i\pi}{2}} \ldots [\text{Recall }e^{-\frac{i\pi}{2}}=cos\frac{\pi}{2}-isin\frac{\pi}{2} = -i]

That is, the modulus is given by \left|H(i{\omega})\right| = \frac{1}{\omega} and the phase by\phi(i{\omega}) = -\frac{\pi}{2}

For example

Ideal Integrator Transfer Function Modulus
Figure 4: Ideal Integrator Transfer Function Modulus


{\int}cos{\omega}t = \frac{1}{\omega}sin{\omega}t = \frac{1}{\omega}cos\left({\omega}t-\frac{\pi}{2}\right)

In the following it is helpful to know the relationship between the phase and any time delay. Now the quantity known as the phase delay,  {\tau}_{\varphi}(\omega), is simply the time delay equivalent to the particular phase at each frequency, for example

sin({\omega}t - {\varphi}) = sin\{{\omega}(t+{\tau}_{\varphi})\} = sin\left\{{\omega}\left(t-\frac{\varphi}{\omega}\right)\right\}[latex]</p> <p>That is, the phase delay, [latex]{\tau}_{\varphi}(\omega), in seconds is defined as {\tau}_{\varphi}(\omega) = -\frac{\varphi(\omega)}{\omega} seconds, where \varphi(\omega) is in radians and \omega is in radians/sec. Alternatively in degrees and Hz we have \tau_{\varphi}(f) = -\frac{\varphi(f)}{360f}  seconds. The phase delay gives the time delay at each frequency. If the phase delay is linear, that is, if it is a straight line versus frequency, then it means that all frequencies are delayed by the same amount. Specifically, if the phase is in the form \varphi(\omega) = a{\omega}+b  then the phase consists of the constant angle b plus a time delay to the whole signal of a seconds.  Note that the special case of constant phase may be expressed as \varphi(\omega) = 0{\omega} + b.

For completeness, there is a similar phase related delay known as the group delay, \tau_g(\omega), which is defined by \tau_g(\omega) = -\frac{d{\varphi}(\omega)}{d\omega} seconds when \varphi is in radians and \omega is in radians/sec or \tau_g(f) = -\frac{1}{360}\frac{d\varphi(f)}{df} when \varphi is in degrees and $f$ is in Hz. This is most useful when dealing with modulated signals, but in the special case where the phase is linear, that is \varphi(\omega) = a\omega + b, then the group delay also gives the delay of the signal.

The transfer function for the Trapezium integrator, which is an order 1 integrator, is shown in Figure 5. The transfer function of all odd order integrators is similar to that of the trapezium integrator. For these the modulus function is like the ideal integrator up to about the quarter sample rate. After that it increasingly diverges.

Transfer function for the Trapezium integrator
Figure 5: Transfer function for the Trapezium integrator

Ignoring the zero frequency phase, which is artificially set to zero by the software, then the phase is a straight line. The phase delay is half of a sample period of the signal together with a constant minus 900. That is the result is the integrated signal delayed by half a sample time. Obviously the trapezium integrator is limited in its frequency range.

The Rectangular integrator (Figure 6) is an order 0 integrator and does not suffer the amplitude error. It has the same phase delay of half of a sample period together with a constant minus 900.

Rectangular integrator
Figure 6: Rectangular integrator

The effects are shown in the integrated sine waves in Figure 7 & 8. With a sample rate of 1000 samples/second at 60Hz (12% of Nyquist frequency) the results of either a rectangular or a trapezoidal integrator are identical for practical purposes, but at 350Hz (70% of Nyquist frequency), the results are very different. Note that the apparent beating in the 350Hz case is just that, apparent. It is the effect of the lower number of samples per period and that the software draws straight lines between points.

600Hz example
Figure 7: 600Hz example
350Hz example
Figure 8: 350Hz example

To illustrate the half order delay a sine wave was digitally integrated using rectangular integration and is shown in Figure 9 (red curve) compared to direct generation of the integrated sine wave (blue curve). That is, the blue curve is the ‘correct’ one in amplitude and time, the integrated signal is ‘delayed’.

Direct generation v rectangular integration
Figure 9: Direct generation v rectangular integration

To compensate for the delay the integrated curve was ‘advanced’ by adjusting its start time ('base' in DATS terms) by one half of a sample period. The result is shown in Figure 10 without the computer drawn straight lines joining successive points. It is clear that the digitally integrated signal produces the integrated values at the midpoints.

Direct generation v rectangular integration (time adjusted)
Figure 10: Direct generation v rectangular integration (time adjusted)

In summary, if we seek to integrate the time history or similar then we should use even order integrators. These will give a result that lags the ‘true’ integration by half a sample time. If we have a signal in the frequency domain then we should use the Omega arithmetic option, which does not have the minor delay.

The following two tabs change content below.

Dr Colin Mercer

Chief Signal Processing Analyst (Retired) at Prosig
Dr Colin Mercer was formerly at the Institute of Sound and Vibration Research (ISVR), University of Southampton where he founded the Data Analysis Centre. He then went on to found Prosig in 1977. Colin retired as Chief Signal Processing Analyst at Prosig in December 2016. He is a Chartered Engineer and a Fellow of the British Computer Society.

Latest posts by Dr Colin Mercer (see all)

0 0 vote
Article Rating
Notify of

This site uses Akismet to reduce spam. Learn how your comment data is processed.

Newest Most Voted
Inline Feedbacks
View all comments
Luiz Eduardo T. Ferreira
Luiz Eduardo T. Ferreira
8 years ago

Dear Dr. Colin,
Please, forgive my poor English.
I’m just trying to introduce myself into FFT calculations and I have an elementary question to you.
As an example, suppose I have a time history of velocities signals (after integration of accelerations) and I take a FFT of these signals, say, up to the Nyquist frequency.
The real parts of my imaginary vector will return to me the amplitudes and the imaginary ones will return me the phases. So now, I will have the amplitudes related with the several frequencies. However, I would like to know how to obtain the transformation of these amplitudes in terms of my field (primary) variable, that is, in terms of velocities, so that I can get the spectrum in terms of velocities. In this case It will be possible to perform inverse FFTs for single values, that is for each frequency of the spectrum?
Luiz Eduardo

Luiz Eduardo
Luiz Eduardo
8 years ago

Dear Dr. Colin,
First of all, I would like to thank you very much for your attention.
Of course you are right. In fact, the FORTRAN routine I´m using gives me the result back by means of two vectors. One of then brings the amplitudes. The other one, the phases. It makes use of the preliminary vectors where real and imaginary parts are stored during data manipulation. This is only a particularity of the routine I´m using.
I´m sorry to neglect clarity by the time I sent you the question.
In fact, the matter is related with the back transformation of each value of the power spectrum obtained while running FFT routine, into primary variable value (velocity, in this case).
I´m just back transforming each value of the power spectrum, as soon as they are obtained, by means of the IFFT equation (from that point).
I would like to know if the association I´m doing is correct.
Best regards,
Luiz Eduardo

Would love your thoughts, please comment.x