HIPASS

Short fom:

float array[*][data_columns] hipass(float array[*][4] filter_parameters, float amplification, float array[*][data_columns] data)

Returns the data array, with all columns passed through the filter specified by the filter_parameters and multiplied by the amplification factor. Each data column is filtered separately.

Range form:

float array[*][data_columns] hipass(float array[*][4] filter_parameters, float amplification, float array[*][data_columns] data, integer first_column, integer last_column = data_columns)

Returns the data array, with columns first_column through last_column passed through the filter specified by the filter_parameters and multiplied by the amplification factor. Each data column in the specified range is filtered separately.

If last_column is omitted, it defaults to first_column (i.e. only first_column is filtered).

All arguments (including all array elements) are converted to float. If conversion fails, the default value 0.0 is used.

Each row in filter_parameters describes a second order IIR (Infinite Impulse Response) filter implementing the difference equation

where x[n] is the input sequence and y[n] is the filter's output. The filter coefficients are listed in the order {b₀, b₁, a₁, a₂}.

All second order filters are applied to the input sequence in parallell. The function result is obtained by summing their outputs and multiplying the sum with the factor -amplification (note the sign!).

Filter parameters can be specified manually or created with make_hipass (recommended).

Tip: The frequency response of any linear, time-invariant filter is just the Fourier transform of the filter's impulse response (the filter's output when the input sequence is a unit impulse, i.e. a single 1 followed by zeros) so to visualize the filter's frequency response, all you have to do is pass a unit impulse through the filter, apply fftp to the output and plot the result. See the example below.

_data = mop("sin()", makevector(100, 0, 2 * PI / 100)) + mop("sin()", makevector(100, 0, 32 * PI / 100))

contains a sine wave with period 100 (normalized frequency 1) mixed with a sine wave with period 100 / 16 (normalized frequency 0.16). To filter out the longer wave, we can use a high pass filter with normalized cut-off frequency 0.25:

Plugged into graph sources and plotted together, the above expressions look like this:

The input is in black, the high-pass output in blue. Both wave components can be seen to have suffered some amplitude loss, but the longer wave much more so than the shorter ones.

To improve the result, we can reduce the width of the transition band by adding more poles to our filter. With 15 poles,

The longer wave component is now almost completely gone, while the amplitude of the shorter component is almost unaffected.

Tip: In cases like this, with a sharply defined, narrow frequency range to filter out, the notch function is usually a better choice than a generic high pass filter.

We can visualize the difference in frequency response between the two filters by passing a unit impulse through both and computing the Fourier transforms of their output.

(where the "+ 1" is to skip the zeroth, constant "harmonic") yield the following amplitude responses (blue for 5 poles, red for 15 poles):

It should be noted that the better frequency resolution of the higher order filter comes at the price of larger phase delays at the high frequencies which we are trying to isolate. With the graph sources