spectrum_level

Spectrum level of a signal x.


[SL,f]=spectrum_level(x,nfft,fs,w,nov) % Matlab & Octave
list <- spectrum_level(x, nfft, fs, w, nov) # R

This replaces Matlab's psd function and returns units in dB re root-Hz.

Input var Description Default value
x is a vector containing the signal to be processed. For signals with multiple channels, each channel should be in a column of x.N/AN/A
nfft is the length of the fft to use. Choose a power of two for fastest operation. Default value is 512.N/A512
fs is the sampling rate of x in Hz. Default value is 1.Hz1
w is the window length. The default value is nfft. If w<nfft, each segment of w samples is zero-padded to nfft.N/Anfft
nov is the number of samples to overlap each segment. The default value is half of the window length.N/A.5 * w
Output var Description Units
SL is the spectrum level at each frequency in dB RMS re root-Hz. The spectrum is single-sided and extends to fs/2. The reference level is 1.0 (i.e., white noise with unit variance will have a spectrum level of 3-10*log10(fs). The 3dB is because both the negative and positive spectra are added together so that the total power in the signal is the same as the total power in the spectrum.dB
fs is the vector of frequencies at which SL is calculated.Hz
  1. Matlab: Use [ ] in any argument to access the default value or just don't specify the trailing arguments if all of the defaults are to be used

Matlab & Octave

The example below uses data from the file testset1.nc. You can download it from the animaltags website's example data sets. If the file is saved in your current working directory, load it via:

testset1 = load_nc(‘testset1.nc’)
[SL,f] = spec_lev(testset1.P.data, 4, testset1.P.sampling_rate)

R

The example below uses data from the file testset1.nc. You can download it from the animaltags website's example data sets. If the file is saved in your current working directory, load it via:

testset1 <- load_nc(‘testset1.nc’)
list <- spectrum_level(x = testset1$P$data, nfft = 4, fs = testset1$P$sampling_rate)

bugs@animaltags.org Last modified: 27 July 2017

  • tagwiki/tools/processing/spec_lev.txt
  • Last modified: 2021/06/23 19:51
  • by ort2