Signal Processing Stack Exchange is a question and answer site for practitioners of the art and science of signal, image and video processing. It only takes a minute to sign up.
Sign up to join this community
I want to extract the envelope of music signals with a sample rate of 44.1kHz.
I used the MATLAB command: abs(hilbert(mysignals)), but the resulting signal is similar to the original signal, not the envelope.
Does taking the Hilbert transform extract the envelope of an audio signal? Or is my calculation incorrect?
That really depends on your definition of "enevelope" and what you need it for.
The Hilbert transform calculates the "analytic" signal, i.e. it calculates a matching imaginary part to a real signal by shifting the phase by 90 degrees in the frequency domain. It's reputation of calculating the "envelope" comes mainly from communication technology. It works very well with a narrow band signals, like amplitude or frequency or phase modulated signals. It's based on the simple fact that the analytic signal of a sine wave is a complex exponentional, the magnitude of which is 1. So the hilbert transform (or to be precise calculating the magnitude of the analytic signal) will eliminate the carrier quite well and leaves you with the modulation signal, or "envelope".
For broad band music, this is mostly not the case and the analytic signal isn't particularly helpful. A cymbal crash is basically a noise burst. The hilbert transform of this is just another noise burst and so is the magnitude of the analytic signal. You can't eliminate a carrier, if there isn't one to start with.
IF you are looking for something that roughly gives you "perceived loudness as a function of time", there many different ways to do this from simple lossy peak detectors all the way to sophisticated perceptual algorithms including cochlea models, binaural & temporal masking, nerve excitation models etc. It really depends on your application.
The Hilbert transform does NOT calculate an analytic signal. If we have a real-valued sequence $x(n)$, and the Hilbert transform of $x(n)$ is $x_H(n)$, the analytic signal, $a(n)$, associated with $x(n)$ is:
$$a(n) = x(n) + jx_H(n)$$
Computing the magnitude of $a(n)$ will give you the instantaneous envelope of the original $x(n)$ sequence. This works great for AM demodulation. If you compute the instantaneous phase of $a(n)$ to get $\phi(n)$, and compute the time derivative of $\phi(n)$, you have FM demodulation.
Not sure about the theory, but here's something that seems to work in practice:
If you divide your audio signal into sub-bands of much less than one octave width (which obviously requires a low frequency limit for a finite number of bands), then the weighted sum of the instantaneous amplitudes from all the complex signal plus Hilbert transforms of all the sub-bands seems to approach something that looks like an instantaneous amplitude of the envelope of the whole audio signal (depending on the lowest sub-band used compared to the lowest frequency present in the audio signal, and the sub-band width(s) in fractions of an octave). You will also have to disregard the start-up transient of the longest Hilbert transform filter used (which again requires some high-pass cut-off with a long enough audio snippet).
ADDED:
Could the issue be a psycho-perceptual one, not a complex analytic one?
Take the sum of 2 audio frequency pure sinusoids of the same magnitude but a couple Hz apart. A human will hear a "beat" in the audio, and the Hilbert transform magnitude will also show the beating as the two analytic exponential spirals go in and out of phase. Now make the 2 sinusoids 20 Hz or more apart, and the Hilbert magnitude will still show a beating in the envelope of the sum, but a human might hear a 2 tone interval with a constant envelope instead, depending on the separation in Hz and semitones.
Previous answers (from @Hilmar, @RickLyons, and @hotpaw2) are very informative. Hope I can also provide some.
First of all, there are more than one way to construct an complex signal based on a real one. Think about
$x(t)=\cos(2\pi f_1 t)\cos(2\pi f_2 t)$, where $f_1<f_2$,
one can define either
$z(t)=\cos(2\pi f_1 t)e^{j2\pi f_2 t}$, or $z(t)=\cos(2\pi f_2 t)e^{j2\pi f_1 t}$
However, Hilbert transform only generates one analytic signal (the first one). According to Bedrosian's product theorem, for real signal $x(t)=a(t)\cos\phi(t)$, its Hilbert transform
$H\{a(t)\cos\phi(t)\}=a(t)H\{\cos\phi(t)\}=a(t)\sin\phi(t)$
if $|F\{a(t)\}|=|A(\omega)|>0$ only for $\omega<\omega_{critical}$,
while $|F\{\cos\phi(t)\}|>0$ only for $\omega>\omega_{critical}$.
In this case, $a_z(t)=|x(t)+jH\{x(t)\}|=a(t)$ will be the desirable envelope. Otherwise, $a_z(t)\ne a(t)$ and might be hard to interpret. Please refer to [2] for more details.
[I'm not quite sure about my following reasoning...]
The analytic signal associated with any real bandpass ($\omega_a<\omega<\omega_b$) signal can be expressed as $z(t)=A(t)e^{j(\omega_a+\phi(t))}=A(t)e^{j\phi(t)}e^{j\omega_a}=z_{base}(t)e^{j\omega_a}$,
where the spectrum of $z_{base}(t)$ and thus $A(t)$ is upper bounded by $\Delta\omega=\omega_b-\omega_a$,
while the spectrum of $e^{j\phi(t)}e^{j\omega_a}$ is lower bounded by $\omega_a$
If $\omega_b<2\omega_a$ (or $\Delta\omega<\omega_a$), then $H\{A(t)\cos(\omega_a+\phi(t))\}=A(t)\sin(\omega_a+\phi(t))$, and Hilbert transform will provide reasonable envelope.
I think this is why to bandpass a broadband signal (e.g., music audio) in a series of octaves before applying Hilbert transform results in reasonable outputs.
[2]: Boashash, B. (1992). Estimating and interpreting the instantaneous frequency of a signal. I. Fundamentals. Proceedings of the IEEE, 80(4).
