Obtaining the envelope of AM signal using hilbert transform
56 views (last 30 days)
I have an amplitude modulated signal and want to extract the envelope
I tried using the hilbert transform and use the absolute value to obtain the amplitude.
The figure shows the original signal (wich is also the real part of the hilbert function) in red the imaginary part in blue and the envelope. I obtain the envelope as
As you can see the envelope is not smooth and shows part of the carrier signal.
I'm wondering how to obtain a better estimate of the amplitude (and phase) of the signal.
David Goodmanson on 26 Mar 2021
Edited: David Goodmanson on 26 Mar 2021
I don't think there is much doubt that the effect is due to a DC offset in the signal. In your figure, both the real and imaginary parts are shifted upward slightly as shown by the position of the positive and negative peaks. This code
x = (0:9999)/10000;
offset = .04;
y = cos(2*pi*10*x) + offset;
plot(x,y); grid on
h = hilbert(y) + i*offset; % imaginary offset post hilbert
plot(x,real(h),x,imag(h),x,abs(h)); grid on
shows the effect, with real in blue and imaginary in red. Setting offset = 0 gives a flat line for the envelope as expected.
The problem is, it is easy enough to presume an offset in the original real signal. But for hilbert, the imaginary part of the analytic signal has no DC offset. In order to reproduce your plot I had to apply a DC offset to the imaginary part after using hilbert.
Paul Hoffrichter on 26 Mar 2021
In general, if you increase the carrier frequency, you will get a better envelope.
Specific to your figure, to remove the ripple, you can take the FFT of your resultant envelope signal in the middle region avoiding the AM amplitude shift from low to high and high to low (which will show up in the fft as the highest frequencies). Between the DC and the highest frequencies should be a fairly distince ripple frequency. You can now apply a notch filter to remove the ripple.