![]() |
![]() |
![]() |
||
![]() |
![]() |
|
Phase Demodulation using the Hilbert transform in the frequency domain
The general ideaA phase modulated signal is a type of signal which contains information in the variation of its phase, an example of a phase modulated signal, in its simplest form, is a single sine wave modulated by another sine wave, such as:
Evidently phase demodulation of a signal involves reconstructing a signal such that one can characterise how the modulated signal's phase changes with time. Phase demodulation is therefore based on this simple idea of setting out to measure how the phase of the signal varies with time. For the above simple phase modulated signal, a pragmatic approach might lead you to consider that the measurement of the phase as in fact being trivial by taking the inverse cosine of the time series Without going into detail, which can be proven to oneself by an interested reader, the inversion of the trigonometric function in the time domain results in an erroneous solution essentially due to the ambiguity of the trigonometric function. For instance, the cosine trigonometric function is ambiguous in that the phase angles cannot be distinguished between being in the 1st and 4th quadrant on the unit circle or similarly between the 2nd and 3rd quadrants. However if we express the above example as a complex exponential,
then characterising the phase at any instant in time could be simply obtained by observing the angle between the real and imaginary value of the complex signal at that same instant in time. Thus expressing a real signal in a complex form, of which the real part is the original signal, is the aim of frequency domain Hilbert transform phase demodulation. This complex signal representation is often referred to as the analytic signal. Therefore it needs to be set about seeking how to change our real signal into its complex form. How to implement (the maths)By utilising Euler's formula,
the above simple example signal, equation (1), can be representing in its analytic form by adding the original signal with the sine of the instantaneous phase (the instantaneous phase being
where Nothing more needs to be discussed about the Hilbert transform itself, suffice to say that it is the technical name of the process to be used here. We will now show a convenient way of constructing the analytic signal for our example signal with a judicious use of the Fourier transform. If we first represent a cosine in its complex form:
Also given that the Fourier transform of a signal is:
Therefore the Fourier transform of (2) is:
where If we now set the negative frequencies of equation (3) above to zero, multiply by 2, then inverse Fourier transform we get a new function
It is seen that How to implement (matlab example)(All steps can be cut and pasted into matlab's desktop window) Create a modulated signal in the same form as in equation (1)
Figure 1. Time series of equation (1). Note that the phase modulation is not readily visible Transform the signal into the frequency domain
Figure 2. Spectrum of phase modulated signal as in equation (1) Set negative frequencies to zero, and double all positive frequencies (remember not to double the zero frequency) and inverse transform back to the time domain to create the analytic signal.
(Note the analytic signal can be created in one step with the use of the matlab command 'hilbert')
Figure 3. One sided spectrum, or spectrum of the so called analytic signal Now calculate the instantaneous angle of the analytic signal and the unwrapped instantaneous phase of the original signal (as the matlab command angle only gives a value between 0:2pi, use the unwrap command to give the non-bound limited phase)
As the instantaneous phase is given by If the carrier frequency is known this can be done by multiplying the carrier frequency by the inverse of the sampling frequency and subtracting from the unwrapped instantaneous phase, or if the carrier frequency is unknown then the linear fit of the unwrapped phase will give the estimate of
Now observe the spectrum of the modulating signal
Figure 4. Spectrum of the phase modulating signal with linear fit of carrier frequency The mass line can be seen to be due to the linear fit not being able to exactly match the linear phase increase. If the actual linear phase is used then the mass line can be seen to be removed:
Figure 5. Spectrum of the phase modulating signal using known carrier frequency If we now observe the phase and amplitude of the modulation signal, they are found to be very close to the actual values
Table 1. Actual and estimated signal parameters after demodulation
Bandwidth considerationsSome considerations on the bandwidth of both the modulating signal and the carrier frequency relationship will now be discussed. Three general rules of 'thumb' will be given for bandwidth's which will commonly result in a phase modulated signal which can be demodulated with conventional techniques.Despite the discussion of the various bandwidth considerations that will be developed below, first and foremost the major requirement for accurate demodulation is the separation of the negative and positive frequency sidebands in the signals spectrum. This criteria is the application of the initial condition of Bedrosian's Theorem [1]. This theorem states in its most succinct form, when applied to a phase modulated signal, that the respective frequency domains of the carrier and modulating functions are non-intersecting and that the frequency of the carrier is higher than the modulating frequency for the general solution of the Hilbert transform to hold [2]. Whether this separation is present is often evident from simply viewing the spectrum, as can be seen in Figure 2 where the positive and negative frequency sideband regions are clearly separated. However when this is not evident from simply viewing the spectrum some rules of 'thumb' which can be applied to help achieve accurate demodulation will be discussed. Fundamentally, for conventional demodulation, the maximum modulating frequency must be at least less than the carrier frequency, however this constraint alone does not provide a signal which is able to be accurately demodulated as only one set of sidebands may only be able to be used in the demodulation. We can in investigate this limitation by firstly expressing a phase modulated signal as its expanded Bessel series.
As it can be seen the spectrum of a modulated signal will have an infinite set of frequencies located at the carrier frequency plus and minus the modulating frequency. Looking at the spectrum it will have components at
Figure 6.(left) This limitation can be discussed with the introduction of a term for the ratio of carrier to modulation frequency being:
It can be seen that first set of sidebands is corrupted by this frequency wrapping and that the sidebands equal or greater than the ratio These limitations and indeed the whole phase demodulation theory can be applied to a modulating signal which is broadband in nature. In that case these bandwidth limitations apply to the highest frequency in the modulating signal.
Figure 7.
The value of Two good rules of thumb that can be derived from these results, that should be adhered to for accurate demodulation are; Firstly
Secondly the modulation amplitude should be limited by:
This now brings up the question of how many sidebands should be included to obtain an accurate demodulation without undue computation cost? The number required for any given degree of accuracy can be obtained from observation of Figure 8 however a general criteria which is often used is only including sidebands up to when any higher sidebands are less than 100 times smaller than the highest amplitude sideband. With the minimum lower limit on the number of sidebands included being 3 pairs. So the last rule of thumb for demodulation is that the number of sidebands that should be included be greater than 3 pairs of sidebands for a modulation amplitude of less than 2.
Figure 8. Estimates of modulation amplitude with increasing number of sidebands used in the demodulation for different modulation amplitudes as shown
Figure 9. Wrapping due to a large modulation amplitude. Lastly filtering should generally be done whenever the full bandwidth is not being used for demodulation, which is normally the case. The signal should be band-pass filtered around the carrier and sidebands, that are to be used in the demodulation, to avoid distortion from out of band frequencies. For example it is shown in Figure 10 the band-pass filtered spectrum for
Figure 10. (left) Spectrum and overlaid band-pass filter. (right) Band-pass filtered spectrum [1] Bedrosian, E.A. [2] Cerejeiras, P., Q. Chen and U. Kahler, [3] Forbes, G.L. and R.B. Randall, ![]() Phase Demodulation by Gareth Forbes is licensed under a Creative Commons Attribution-Non-Commercial 2.5 Australia License. |