DADiSP Online Help
Click here to see this page in full context

DADiSP Worksheet Functions > Function Categories > Fourier Transforms and Signal Processing > BLACKMANHARRIS

 

BLACKMANHARRIS

Purpose:

Multiplies a series with a 4 term Blackman-Harris window.

Syntax:

BLACKMANHARRIS(series, ampflag, "sym")

series

-

A series or array.

ampflag

-

Optional. An integer, the amplitude correction flag:

0

:

do not correct amplitude (default)

1

:

correct amplitude

2

:

correct RMS amplitude

3

:

correct mean squared amplitude

"sym"

-

Optional. A string, the symmetry flag:

"symmetric"

:

Starting and ending points are equal (default).

"periodic"

:

Periodically extended window where the symmetric end point is removed. Conforms to the ISO standard.

"iso"

:

Same as "periodic".

Alternate Syntax:

BLACKMANHARRIS(N, ampflag, "sym")

N

-

An integer, the length of the window.

ampflag

-

Optional. An integer, the amplitude correction flag:

0

:

do not correct amplitude (default)

1

:

correct amplitude

2

:

correct RMS amplitude

3

:

correct mean squared amplitude

"sym"

-

Optional. A string, the symmetry flag:

"symmetric"

:

Starting and ending points are equal (default).

"periodic"

:

Periodically extended window where the symmetric end point is removed. Conforms to the ISO standard.

"iso"

:

Same as "periodic".

Returns:

A series or array.

Example:

W1: gsin(1000, .001, 45)

W2: spectrum(blackmanharris(W1))

W3: spectrum(blackmanharris(W1, 1))

 

The MAX of W2 == 0.359 and the MAX of W3 == 1.0. The amplitude of the spectrum in W3 has been corrected to take into account amplitude effects of the symmetric 4 term Blackman-Harris window.

 

 

where n is the nth point (1 <= n <= N) and N is the number of points to generate.

Example:

W4: blackmanharris(1000, "periodic")

 

Creates a 1000 point periodic Blackman-Harris window that conforms to the ISO 18431-1 standard.

 

 

where n is the nth point (1 <= n <= N) and N is the number of points to generate.

Remarks:

If ampflag == 1, the correction factor is the mean of the spectral window. This assures that the spectrum of a sinusoid of amplitude A has a peak of A.

 

If ampflag == 2, the correction is applied as follows:

 

 w = blackmanharris(s) * rms(s) / rms(blackmanharris(s))

 

This assures that:

 

sqrt(area(psd(w))) == rms(s) approximately

 

If ampflag == 3, the correction is applied as follows:

 

w = winfun(s) / sqrt(mean(win * win) 

 

where win is the windowing function.

 

The "sym" flag controls the window symmetry as follows:

 

"Symmetric" sets the last point to be the same value as the first point. For an N point symmetric window, a N-1 point periodic window is effectively created and the Nth point is set to the same value as the first point.

 

"Periodic" or "iso" creates a periodic window function useful in spectrum analysis applications. "Periodic" or "iso" conforms to the ISO 18431-1 standard for windowing functions.

 

The Hamming, Hanning, Flattop and Blackman windows are part of the family of cosine window functions. The ISO 18431-1 standard periodic form of these windowing functions are defined by:

 

 

where K is the number of window coefficients and N is the length of the window. The symmetric form of the window can be constructed by setting N to N-1.

 

For the Blackman-Harris window:

 

 

See GBLACKMANHARRIS to generate a Blackman-Harris window.

 

See BLACKMAN to multiply a series with a 3 term Blackman window.

See Also:

BLACKMAN

FLATTOP

GBLACKMAN

GBLACKMANHARRIS

HAMMING

HANNING

KAISER

PSD

SPECTRUM

WINFUNC