Next Article in Journal
Spatial-Temporal Dynamics of Cropping Frequency in Hubei Province over 2001–2015
Previous Article in Journal
Computational Experiments on the Step and Frequency Responses of a Three-Axis Thermal Accelerometer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circular Regression in a Dual-Phase Lock-In Amplifier for Coherent Detection of Weak Signal

1
Laboratoire de Physico-Chimie de l’Atmosphère, Université du Littoral Côte d’Opale, 59140 Dnkerque, France
2
Laboratoire d’Informatique, Signal et Image de la Côte d’Opale (LISIC), Université du Littoral Côte d’Opale, 62228 Calais, France
*
Authors to whom correspondence should be addressed.
Sensors 2017, 17(11), 2615; https://doi.org/10.3390/s17112615
Submission received: 23 September 2017 / Revised: 9 November 2017 / Accepted: 11 November 2017 / Published: 14 November 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
Lock-in amplification (LIA) is an effective approach for recovery of weak signal buried in noise. Determination of the input signal amplitude in a classical dual-phase LIA is based on incoherent detection which leads to a biased estimation at low signal-to-noise ratio. This article presents, for the first time to our knowledge, a new architecture of LIA involving phase estimation with a linear-circular regression for coherent detection. The proposed phase delay estimate, between the input signal and a reference, is defined as the maximum-likelihood of a set of observations distributed according to a von Mises distribution. In our implementation this maximum is obtained with a Newton Raphson algorithm. We show that the proposed LIA architecture provides an unbiased estimate of the input signal amplitude. Theoretical simulations with synthetic data demonstrate that the classical LIA estimates are biased for SNR of the input signal lower than −20 dB, while the proposed LIA is able to accurately recover the weak signal amplitude. The novel approach is applied to an optical sensor for accurate measurement of NO 2 concentrations at the sub-ppbv level in the atmosphere. Side-by-side intercomparison measurements with a commercial LIA (SR830, Stanford Research Inc., Sunnyvale, CA, USA ) demonstrate that the proposed LIA has an identical performance in terms of measurement accuracy and precision but with simplified hardware architecture.

1. Introduction

Lock-in amplifier (LIA) is an effective device capable of recovering weak signal buried in high noise level. This technique relies on locking an input signal to measure at a specific reference frequency. The signal is extracted from noise using phase sensitive detection (PSD) at this reference frequency [1,2]. In a LIA, the weak input signal is multiplied by a sine-wave signal that acts as reference signal [3]. This multiplication is processed with an analog multiplier, a digital switch or a digital multiplier [4]. A low-pass filter, with an appropriate cut-off frequency, is used to extract the mean value of the product between these two signals. The output filter value is proportional to the amplitude of the input signal and depends on a phase delay. This phase delay between the input signal and the sine-wave reference signal results from frequency fluctuations [5]. In this context it is important to set the reference phase to the value of the input phase signal in order to correctly recover the amplitude [6].
Phase locked loop (PLL) can be implemented in LIA to actively track the phase of the input signal. In this approach an internal reference generator is adjusted to match the phase of the input signal [7]. Andrea De Marcellis et al. [8] added a “phase alignment channel” and a “90 phase shifter channel” to the LIA. These channels generate a reference in-phase signal and a quadrature signal with a 90 shifted phase. The reference signal of the shifted channel and the input signal are provided to the PSD of the LIA. This novel LIA realizes an automatic phase alignment but needs, however, complex analog components and has a long response time.
Dual-phase LIA is an architecture widely used to remove the effects of the phase delay ϕ k . In this architecture, two multiplications are performed with respectively a reference signal and the same signal in quadrature. A low pass filter is applied to the signals to produce the in-phase component I = A cos ( ϕ k ) and the quadrature component Q = A sin ( ϕ k ) . A is the amplitude of the input noisy signal, and ϕ k = arctan ( Q / I ) the phase delay. The in-phase and the quadrature components are squared and summed in order to remove the phase delay and to provide the input signal amplitude: R = I 2 + Q 2 = A [9]. However, R is a nonlinear function of I and Q. It can be shown that for the classical hypothesis of an additive Gaussian noise on the two components in quadrature, the random observations of R are distributed according to a Rician distribution [10]. When the signal-to-noise ratio (SNR) is high, the Rician distribution is similar to a Gaussian distribution and A can be estimated with observations of R. However, the Rician distribution deviates significantly from the Gaussian distribution at low SNR and the estimation of A with observations of R is biased. The bias is a function of the noise noise level [10].
In the present work, we propose to determine the signal amplitude with observations R c = I cos ( ϕ ^ k ) + Q sin ( ϕ ^ k ) , where ϕ ^ k is an estimate of ϕ k . We show in Figure 1 that the observations R c are more accurate than the observations R provided by a classical coherent detection method. Furthermore we show in the next section that these observations are unbiased even at low SNR. This new approach requires, however, a precise estimate of the phase delay ϕ k . This phase delay between the reference and the input signal is mainly due to frequency instability in the input noisy signal and is assumed to vary slowly [11]. It results in a linear phase delay variation. In this article, we introduce a linear-circular regression model to accurately estimate phase delays of the input noisy signal [12].
Signal processing of angular data (here the phase delay) relies on the modeling of angles in a circular Bayesian framework [13,14,15]. The noise is assumed to follow a circular normal von Mises distribution. This distribution acts in the circular domain as the normal distribution in the linear domain for observations defined on the real line [16,17]. Several nonlinear techniques have been proposed in the published works for angular data. These techniques include circular state estimation [18,19] and recursive filtering [20,21,22]. In the present work the circular regression of angular data is used to estimate the phase shift of the input signal [23].
The paper is organized as follows: the Section 2 describes the architecture of the proposed circular regression-based dual-phase digital LIA. The linear-circular model and method used to accurately estimate the phase delay are presented in the Section 3. In the Section 4, the proposed LIA is assessed using synthetic and real experimental data. Finally, high-sensitivity measurements of NO 2 trace gas in the atmosphere are presented and the performance of the proposed LIA is evaluated in comparison with a commercial LIA (SR830, Stanford Research Inc.).

2. Digital Signal Detection

A typical dual-phase LIA (Figure 2a) and the proposed circular regression-based dual-phase LIA (Figure 2b) are schematically depicted in Figure 2. We assume that an input signal s ( t ) contains a carrier of amplitude A at a frequency f. This signal is digitized at a sampling frequency f s by an analog–to-digital convertor (ADC). This discrete sampled signal s n can be expressed as a function of the sampling period T s with the following expression:
s n = A s i n 2 π f n T s + η n
where η n is a Gaussian noise. The general aim of lock-in detection is to extract the amplitude value A, buried in noise. We define the SNR of the input signal as follows:
S N R = A σ g
where σ g is the standard deviation of the Gaussian noise η n . The SNR in dB is:
S N R d B = 10 log 10 A 2 σ g 2
The in-phase i n component and the quadrature q n component are respectively defined as the multiplication of the input signal with the reference signal (at frequency f r ) and the multiplication of the input signal with the reference signal shifted by π / 2 . i n and q n are expressed as follows:
i n = A 2 cos 2 π ( f r f ) n T s + cos 2 π ( f r + f ) n T s + η n i
q n = A 2 sin 2 π ( f r f ) n T s + sin 2 π ( f r + f ) n T s + η n q
where η n i and η n q are independent Gaussian noises with a variance of σ g 2 / 2 . In order to filter out the frequency component f r + f , T w f s samples of the signals i n and q n are summed on a working window T w . This sum acts as a low-pass filter. The integration time T w corresponds to the time constant in classic dual-phase LIA (Figure 2a). The choice of T w will be investigated in a next section (Section 4.3). Two components values, I k and Q k , are processed at each time instant k, multiple of the period T w . Expressions of the two components after integration by the low-pass filter are given below:
I k = A f s T w 2 cos ϕ k + η k I
Q k = A f s T w 2 sin ϕ k + η k Q
where η k I and η k Q are Gaussian noises with a variance of f s T w σ g 2 / 2 . ϕ k is the phase delay between the input signal and the reference. This phase delay evolves linearly with time when the frequency difference f r f is constant. In a classical dual-phase LIA (Figure 2a), the noisy observations of I k and Q k are processed with the following expressions to observe noisy values of the signal amplitude and the phase delay:
R = I k 2 + Q k 2
ϕ ˜ k = arctan Q k I k = ϕ k + Δ ϕ k
where Δ ϕ k is the phase noise related to Gaussian noises η k I and η k Q . Figure 2b shows the proposed LIA architecture. In this architecture, we notice ϕ ^ k the phase delay obtained with the linear-circular regression estimator. With the proposed approach, noisy observations of I k and Q k are processed with the following expression to observe the signal amplitude:
R k c = I k cos ( ϕ ^ k ) + Q k sin ( ϕ ^ k )
With a linear-circular regression approach the estimated phase delay ϕ ^ k is between ϕ ˜ k and ϕ k . The estimate ϕ ^ k tends to ϕ k when the number of observations used to process the circular estimate is high. In the opposite, the estimate ϕ ^ k tends to ϕ ˜ k when the number of observations is low. The corresponding expressions for R k c are given below.
  • When ϕ ^ k = ϕ ˜ k
    (11) R k c = I k cos ( ϕ ˜ k ) + Q k sin ( ϕ ˜ k ) (12) = I k I k I k 2 + Q k 2 + Q k Q k I k 2 + Q k 2 (13) = I k 2 + Q k 2 = A f s T w 2 + η k C
    where η k C is a Rician noise. In this case, R k c = R k and the proposed approach is equivalent to the classical dual-phase approach. The noise on the observed signal is distributed according to a Rician distribution. This approach is called incoherent detection (ID) in the rest of the article.
  • When ϕ ^ k = ϕ k
    (14) R k c = I k cos ( ϕ k ) + Q k sin ( ϕ k ) (15) = A f s T w 2 cos 2 ϕ k + A f s T w 2 sin 2 ϕ k + η k I cos 2 ϕ k + η k Q sin 2 ϕ k (16) = A f s T w 2 + η k C
    In this case, η k C = η k I cos 2 ϕ k + η k Q sin 2 ϕ k is a Gaussian noise. This is the proposed approach. In this case the noise on the observed signal amplitude is distributed according to a Gaussian distribution. This approach is called coherent detection (CD).
We show in Figure 3 the mean value of A obtained with observations of R k c . The x-axis is the SNR of the components I k and Q k . The black line shows the values of the input signal amplitude A. The black line is also the mean value of A obtained with the CD method for a Gaussian noise on the observations. The red line represents the mean value of A obtained with the ID method for a Rician noise on the observations. These curves show that the ID method (incoherent detection) is biased in case of weak signal detection, especially when the SNR of the I k and Q k components is less than 5 dB. These results are obtained for a Rician probability distribution given by the following equation:
f ( x | A , σ ) = x σ 2 exp x 2 A 2 2 σ 2 I 0 x A σ 2
where σ 2 is the variance of the Gauss noises η k I and η k Q . I 0 is the modified Bessel function of first kind and order zero. The mean value μ of the Rician distribution is obtained with the following expression:
μ = σ π 2 L 1 / 2 A 2 2 σ 2
where L q ( x ) is the Laguerre polynomial.

3. Circular Regression

Circular regression has been used in several applications such as wind and wave direction estimation [23,24], GNSS altimetry [25]... In this article we use a linear-circular regression to estimate the phase delay between a reference signal and the signal to be measured. In this model the abscissae is linear and the ordinate is an angle. We can find in the published works three different circular regressions; the circular-linear regression, the circular-circular regression and the linear-circular regression [16].
The main difficulty in a linear-circular regression is the circular nature of the variable response. One way to address the circularity is to assume that the noise is distributed according to a von Mises distribution. We can find in the published works different maximum likelihood approaches for the regression parameters estimation [26,27,28,29]. In this article we derived a similarity function from the likelihood and we propose to maximize this similarity function with a Newton Raphson algorithm.
In this section, we derive an accurate estimate of the phase delay ϕ k of the input signal. We consider a working window T p and observations obtained with a period T w . The number of observations used to process the linear-circular regression thus is n = T p / T w . The choice of the working window T p is discussed in the next section, dedicated to simulation assessment of the proposed LIA approach. The observed phase delay ϕ ˜ k is a circular variable which takes its values on the circumference of a circle. This variable is assumed to follow a linear-circular model with an additive angular noise following a circular von Mises distribution. The noisy linear-circular model is described as follows:
y k = α + β x k + ϵ k ( m o d 2 π )
where y k models the observations of the phase delay ϕ ˜ k at x k . ϵ k is a zero-mean noise distributed according to a von Mises distribution with concentration parameter κ . α and β are the y-intercept and slope of the linear-circular regression. In our application the slope is equal to β = 2 π ( f r f ) . m o d 2 π is the modulo of 2 π . The von Mises distribution is defined by the following equation:
f ( y k / α , β , κ ) = 1 2 π I 0 ( κ ) e x p ( κ cos ( y k ( α + β x k ) ) )
The likelihood of n samples y k is defined by:
L ( y 1 , , y n ) = f ( y 1 / α , β , κ ) f ( y n / α , β , κ )
The likelihood of n samples for the von Mises distribution is given by:
L ( y 1 , , y n ) = 1 2 π I 0 ( κ ) n e x p κ k = 1 n cos ( y k ( α + β x k ) )
The log-likelihood function is given by:
l ( y 1 , , y n ) = log L ( y 1 , , y n ) = n log 2 π I 0 ( κ ) + κ k = 1 n cos ( y k ( α + β x k ) )
In order to derivate the maximum likelihood estimate α ^ of α , the log-likelihood function is derived with respect to α , and the result is given in Equation (24):
l ( y 1 , , y n ) α = κ k = 1 n sin ( y k ( α + β x k ) )
Equation (24) is set to zero to obtain the estimate of α :
l ( y 1 , , y n ) α = 0 cos ( α ) k = 1 n sin ( y k β x k ) = sin ( α ) k = 1 n cos ( y k β x k )
Finally the estimate of α is obtained:
α ^ = arctan * k = 1 n sin ( y k β x k ) k = 1 n cos ( y k β x k )
where A r c t a n * ( ) is the “quadrant-specific” inverse of the tangent. It is shown in [19] that A r c t a n * ( ) is the maximum likelihood estimate. In order to compute the maximum likelihood estimate β ^ of β , the 1st and 2nd order derivatives of the log-likelihood with respect to β are calculated:
l ( y 1 , , y n ) β = κ k = 1 n x k sin ( y k ( α + β x k ) )
2 l ( y 1 , , y n ) β 2 = κ k = 1 n x k 2 cos ( y k ( α + β x k ) )
In this context the second order Taylor approximation of l ( ) at y 1 , , y n can be obtained with the Newton Raphson algorithm and the values of α and β that maximize the log-likelihood are recursively obtained by:
α i = arctan * k = 1 n sin ( y k β i 1 x k ) k = 1 n cos ( y k β i 1 x k )
β i = β i 1 + 1 2 k = 1 n x k sin ( y k ( α i + β i 1 x k ) ) k = 1 n x k 2 cos ( y k ( α i + β i 1 x k ) )
α i and β i converge to the maximum likelihood estimates α ^ and β ^ . The circular regression is equivalent to the linear regression when y k is unwrapped. In this context the variance of the estimate β ^ is defined by:
v a r ( β ^ ) = σ l 2 k = 1 n k T w t ¯ 2
The link between the standard deviation σ l in the linear domain and the parameter of concentration κ in the circular domain is defined by [17]:
σ l 2 = 2 ln ( A ( κ ) )
with A ( κ ) = I 1 ( κ ) I 0 ( κ )
where I n ( ) is the modified Bessel function of the first kind of order n. Finally the relation between κ and the SNR of s ( t ) can be approximated using the following linear expression:
κ = f s T w 2 S N R

4. Experimental Assessments

In order to characterize the ID and CD approaches and to optimize the typical parameters such as T w , T p , synthetic data are simulated with MATLAB in a first step. In a second step, the proposed CD approach is applied to an optical NO 2 sensor for experimental validation.

4.1. Simulation Assessment Using Synthetic Data

The ID and CD methods are assessed through MATLAB simulations. We show in Figure 4a the ID and CD estimates of the input signal amplitude A. These estimates are plotted in Figure 4a as a function of the SNR. The amplitude A is fixed to 1.0 V. The difference of frequency between the signal and the local reference is f r f = 2 Hz (a typical frequency shift of a mechanical chopper usually used for amplitude modulation in optical sensors [11]). The linear-circular regression is processed with phase delay observations obtained in a working window T p = 1 s. The sampling frequency is f s = 44 kHz. The observations of A provided by the ID and CD methods are averaged within an observation window T o = 1 s to provide an estimate A ^ :
A ^ = T w T o k = 1 T o T w A k
As can be seen in Figure 4a, for weak signals with a SNR lower than 5 dB, the ID estimate deviates significantly from the real signal amplitude when the integration period is T w = 0.001 s (black line) or T w = 0.01 s (red line). This effect is due to the biased observations provided by the ID method. For a period T w = 0.01 s (grey line), the CD estimate does not deviate significantly from the real signal amplitude even for a SNR lower than 5 dB. In this case, the observations provided by the CD method are indeed unbiased.
In order to improve the ID estimate and to remove the bias, the integration period T w is increased. In this case we can observe in Figure 4a that the amplitude is underestimated for T w = 0.1 s (blue line) or T w = 0.25 s (green line). This underestimation is explained in Figure 4b with the study of the lock-in detection process in the frequency domain.
In lock-in detection, the incoming signal is down converted at frequency f r f and then filtered out by a low-pass filter (sum function). The transfer function of the low pass filter is given by the following equation:
H ( f ) = sin π f N T s sin π f T s exp j 2 π f ( N 1 ) T s
where N = f s T w , is the number of observations in an integration period T w . Figure 4b shows the transfer function of low-pass filters with different integration periods . As f r f = 2 Hz, the amplitude output of the low-pass filter with integration periods 0.25 s, 0.1 s, 0.01 s and 0.001 s are respectively 0.60, 0.93, 0.99 and 1. These amplitudes are close to the amplitude estimated with the ID method and it demonstrates that a period T w larger than 0.25 s results in an underestimated signal amplitude with a bias ratio of 0.60. In addition, larger T w results in a longer response time of the LIA detection. As a conclusion, using the ID method, the amplitude is either biased at low signal to noise ratio (for large T o ) or underestimated at any signal to noise ratio (for large T w ). In this context, the only way to improve the estimation process is to estimate the signal amplitude with CD observations processed with long observation window T o and short integration time T w .
As a conclusion, using the ID method, the amplitude is either biased at low signal to noise ratio (for large T o ) or underestimated at any signal to noise ratio (for large T w ).
The amplitude observation window T o influences the performance of the CD approach for recovery of weak signals amplitude. We show in Table 1 the minimum SNR required for the estimation of the input signal amplitude with an uncertainty lower than 10 % . In this simulation, the integration time T w is set to 10 ms, T p = 1 s and the frequency difference is fixed to 2 Hz. The minimum SNR of the input signal varies from −14 dB to −22 dB when T o varies from 0.25 s to 1 s. It tends to be constant (−24 dB) with a period T o larger than 2 s. Increasing the size of T o from 0.25 s to 1 s with a T w = 0.01 s can enhance the estimation accuracy, but for an observation window superior to 2 s the performance remains the same.
The estimation accuracy of the CD method can be enhanced with the phase observation window T p . The dependence of the standard deviation (SD) of the estimated slope (that reflects the frequency difference) over the phase observation window T p is investigated Figure 5a. We present in Figure 5a the standard deviation of the estimates α ^ and β ^ as a function of the phase observation window T p for a SNR of −22 dB. We show in Figure 5a that the SD of α ^ and β ^ decrease when the period of observation T p increases. In practice a period of observation T p = 1 s and an integration period T w = 10 ms are used. This choice is based on the simulations results. We indeed observe in Figure 5a that the SD of the estimated slope with T p = 1 s is 0.37 rad/s. This SD value corresponds to a frequency difference of 0.06 Hz which is sufficiently accurate for a typical frequency shift of 2 Hz.
Figure 5b shows a measurement example obtained with simulated data. The phase delay of the input signal is estimated with the circular regression of the CD method discussed in Section 3. For this experiment the amplitude A changes from 0.05 V to 0.10 V and to 0.20 V and the corresponding SNR of the input signals are, respectively, −26 dB, −20 dB and −14 dB. The difference of frequency between the signal and the local replica is fixed to f l f = 2 Hz. For the CD and ID approaches, an integration time T w =10 ms and a phase observation window T p = 1 s are applied. The signal amplitudes estimated by ID approach for a SNR of −26 dB, −20 dB, −14 dB are respectively 0.10 ± 0.01 V, 0.13 ± 0.01 V and 0.21 ± 0.02 V. The signal amplitudes estimated by CD approach are 0.05 ± 0.02 V, 0.10 ± 0.02 V and 0.20 ± 0.02 V and are well consistent with the true values. As expected, the estimated signal amplitudes with the ID approach are biased for a SNR less than −20 dB. The same conclusions are observed with LIA detection of a real signal, as presented in the next paragraph.

4.2. Application to High-Sensitivity and Accurate Measurements of NO 2 at ppbv Concentration Level in the Environment

We implement the proposed “circular phase processing” LIA approach in an optical sensor for accurate measurements of low trace gas concentration of NO 2 in the atmosphere. The optical sensor is based on photoacoustic spectroscopy (PAS). In this sensor a pressure wave (acoustic signal) is generated by thermal expansion resulting from absorption of a power-modulated light by target gas [30]. The resulting acoustic signal, proportional to the target gas concentration, is detected with microphones.
Figure 6a describes a classical LIA architecture of a PAS sensor. A laser beam, modulated with a mechanic chopper at a “reference frequency”, is used to probe absorption of NO 2 . The reference frequency is then used to demodulate the signal with a traditional lock in amplifier. The photoacoustic (PA) cell is a single-pass cylinder which has a length of 26 mm and a diameter of 6 mm [31]. The proposed LIA architecture is described in Figure 6b. The proposed approach uses a fix frequency and no tracking of the reference frequency is necessary.
The detailed experimental set-up is shown in Figure 7. A high power laser diode (NDB7875, Nichia, Tokushima, JAPAN) operating at a wavelength of 444 nm is used as light source. The laser light, driven by a laser diode controller (6340, Arryo instrument, San Luis Obispo, CA, USA), is amplitude-modulated using a mechanical chopper (3501, New Focus, Inc., San Jose, CA, USA). The modulation frequency is matched to the fundamental longitudinal resonance frequency of the cylinder (in order to be an acoustic resonator) at 6260 Hz. Electric condenser microphones (EK-23329-P07, Knowels, Itasca, IL, USA) are used to record the acoustic signal in the PA cell.
Photoacoustic signals are sampled and digitized with a data acquisition system (MX411-P, HBM, Darmstadt, Germany). The following configurations are adopted: a 24-bit ADC at a sampling rate of 96 kHz. As a reference instrument used for inter-comparison, a classical dual-phase LIA SR830 (Stanford Research Inc.) is used with the following configuration: 16-bit DAC and a sampling rate of 256 kHz. The proposed CD approach is implemented in real time with MATLAB.
For NO 2 concentration calibration, we use a standard reference concentration of 10 ppmv (parts per million by volume) NO 2 , diluted with pure nitrogen (N 2 ). We obtain different NO 2 concentrations varying from 6 ppbv to 100 ppbv (parts per billion by volume) using a dilutor (Modele PPA 2000M). The photoacoustic signal is calibrated with measurement of NO 2 concentration provided by an analyzer (AC-31M). The ID and the proposed CD approaches are used to measure different NO 2 concentrations.
Figure 8a shows the estimated NO 2 concentrations obtained with the ID and CD methods as a function of the NO 2 concentration measured by the NOx analyzer AC-31M. For this experimentation the period of observation of the phase delay is T p = 1 s. For both ID and CD measurements the integration time is T w = 0.01 s. The estimated concentration is the average of the measurements obtained over a period T o = 1 s. In Figure 8a the error bar is the root mean square error between the NO 2 concentration measured by the AC-31M analyzer and the concentration estimates obtained with the CD and ID methods. At low concentration, the error bars of the estimates obtained with the ID method are larger than the error bar of the estimates obtained with the CD approach. The inaccuracy of the ID approach in this case is due to the estimation bias of the method. It is obvious that the ID approach can not accurately recover NO 2 concentration less than 20 ppbv. The CD method overcomes this limitation. These results are in accordance with the previously presented theoretical study and assessment using synthetic data.
Figure 8b,c show the observed and estimated phase delay for a NO 2 concentration of respectively 76 ppbv and 6 ppbv. We can observe in Figure 8b,c that the observations of phase delay are more noisy when the concentration of NO 2 decreases. We can observe in Figure 8c that the linear-circular regression is not affected by the 2 π transitions due to the periodic nature of the angular observations. Figure 8b,c, the cycle of the observed phase is 0.56 s, which corresponds to a difference frequency ( f r f ) of 1.79 Hz and an estimated slope β ^ of −11.24 rad/s.

4.3. Comparison Measurements with Commercial LIA

An inter-comparison of NO 2 measurement at low concentrations is carried out using a commercial LIA (SR830) and the proposed LIA architecture. In order to assess the proposed LIA architecture, the SR830 LIA is implemented in two different operation modes. In the first implementation (called external mode SR830E), a link between the chopper and the LIA (as shown in Figure 7, the green arrow) provides an accurate measurement of the PA signal frequency. In the second implementation using an internal mode SR830I, the reference frequency is fixed to 6260 Hz and the above link is removed. These two implementations are compared with the proposed CD method when the integration time is increased.
Figure 9 shows the experimental results for the measurement of a NO 2 concentration of 24.5 ppbv. We show in Figure 9 the assessments of the three different methods and several time constants. The mean value and the standard deviation (error bars) of each estimate are processed in an experimentation of 60 s. In this experiment the time constant Figure 9 is the period T w for SR830E and SR830I LIA. For CD implementation T w is fixed to 0.01 s, T p is fixed to 1 s and the time constant (the abscissa Figure 9) is the period T o .
The results presented in Figure 9 show that the LIA SR830I in internal mode underestimates the concentration when the time constant is >100 ms. This underestimation is due to the frequency difference between the reference frequency and the modulation frequency of the mechanic chopper.
Figure 10a,b show the distribution histogram of the estimated NO 2 concentrations by CD and SR830E LIA with a time constant of 1000 ms. The distribution histograms of the measured concentration are fitted to a Gaussian profile. The mean concentration of 23.9 ppbv obtained by CD method results in an accuracy of 0.6 ppbv (Figure 10a) with a precision of 1.8 ppbv. The mean concentration obtained by the SR830E LIA shown in Figure 10b is 25.2 pppv. Its measurement accuracy and precision are found to be 0.7 ppbv and 1.4 ppbv, respectively.
As can be seen in Figure 9 and Figure 10, the proposed CD method and the LIA SR830E working in external mode have close performance. The advantage of the proposed CD method is that the linear-circular regression accurately estimates the phase delay and removes the noisy phase fluctuations. The implementation of the proposed CD approach is easier and more robust because it only requires a fix reference frequency and simplified hardware components, while SR830 adopts complicated hardware such as low noise amplifier and notch filter to reduce the noise in the input signal and a phase locked loop to track the phase delay.

5. Conclusions

We introduce in the present work a new dual-phase Lock-In Amplification (LIA) processing for coherent detection. The proposed approach relies on precise estimation of the phase delay between the input signal and a reference signal. In our implementation we propose a linear-circular regression estimate of the phase delay. We show with theoretical results and synthetic experimentation that the proposed coherent detection estimate is more accurate than the classical incoherent detection estimate. We show that the proposed approach accurately recovers weak signal intensity at low SNR and outperforms the classical incoherent detection approach for a SNR of the input signal lower than −20 dB.
Coherent detection and incoherent detection approaches are then applied to a PAS-based optical sensor for NO 2 concentration measurements in the atmosphere. The estimated concentrations are assessed with a reference instrument of NO 2 analyzer. The measurements carried out at different NO 2 concentrations indicate that the estimated concentration from incoherent detection approach is biased when NO 2 concentration is less than 20 ppbv. We show that the proposed coherent detection approach overcomes this limitation. This result is in accordance with the theoretical study and assessment using synthetic data.
The measurements comparison, for a constant NO 2 concentration of 24.5 ppbv and different time constants of integration, indicates that the proposed coherent detection LIA and the commercial LIA SR830 have close performances in terms of measurement accuracy and precision. However, the proposed LIA architecture is simpler and the processing is more robust. In future work, the proposed LIA approach using linear-circular regression will be implemented in a field programmable gate array for atmospheric measurements of trace gases pollutants.

Acknowledgments

The authors acknowledge financial support from the French national research agency (ANR) under the CaPPA (ANR-10-LABX-005) and the MABCaM (ANR-16-CE04-0009) contracts. The authors thank the financial support from the CPER CLIMIBIO program. G. W. thanks the “Pôle Métropolitain de la Côte d’Opale” (PMCO) and “the Région Hauts-de-France” for the PhD fellowship support.

Author Contributions

Serge Reboul and Weidong Chen supervised and managed the experimental work and manuscript. Gaoxuan Wang developed the instrument, managed the measurement campaign, processed the data and prepared the manuscript. Jean-Bernard Choquel and Eric Fertein participated in the instrumental development or measurement campaign.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Masciotti, J.; Lasker, J.; Hielscher, A. Digital lock in detection for discrimi-nating multiple modulation frequency with high accuracy and computational efficiency. IEEE Trans. Instrum. Meas. 2008, 57, 182–189. [Google Scholar] [CrossRef]
  2. SRS. About Lock-In Amplifiers. Available online: http://www.thinksrs.com/downloads/PDFs/ApplicationNotes/AboutLIAs.pdf (accessed on 20 September 2017).
  3. Barragan, L.A.; Artigas, J.I.; Alonso, R.; Villuendas, F. A modular low cost digital signal processor-based lock-in card for measuring optical attenuation. Rev. Sci. Instrum. 2001, 72, 247–251. [Google Scholar] [CrossRef]
  4. Burdett, R. Handbook of Measuring System Design; John Wiley and Sons Inc.: Chichester, UK, 2005. [Google Scholar]
  5. Son, H.H.; Jung, I.I.; Hong, N.P.; Kim, D.G.; Choi, Y.W. Signal detection technique utilising ‘lock-in’ architecture using 2 ωc harmonic frequency for portable sensors. Electron. Lett. 2010, 46, 891–892. [Google Scholar] [CrossRef]
  6. Barone, F.; Calloni, E.; DiFiore, L.; Grado, A.; Milano, L.; Russo, G. High performance modular digital lock in amplifier. Rev. Sci. Instrum. 1995, 66. [Google Scholar] [CrossRef]
  7. Hu, A.; Chodavarapu, V. CMOS Optoelectronic Lock-In Amplifier With Integrated Phototransistor Array. IEEE Trans. Biomed. Circuits Syst. 2010, 5, 274–280. [Google Scholar] [CrossRef] [PubMed]
  8. De Marcellis, A.; Ferri, G.; D’Amico, A.; Di Natale, C.; Martinelli, E. A Fully-Analog Lock-In Amplifier with Automatic Phase Alignment for Accurate Measurements of ppb Gas Concentrations. IEEE Sens. J. 2012, 12, 1377–1383. [Google Scholar] [CrossRef]
  9. Macias-Bobadilla, G.; Rodríguez-Reséndiz, J.; Mota-Valtierra, G.; Soto-Zarazúa, G.; Méndez-Loyola, M.; Garduño-Aparicio, M. Dual-Phase Lock-In Amplifier Based on FPGA for Low-Frequencies Experiments. Sensors 2016, 16, 379. [Google Scholar] [CrossRef] [PubMed]
  10. Sijbers, J.; den Dekker, A.J.; Scheunders, P.; Dyck, D.V. Maximum-likelihood estimation of Rician distribution parameters. IEEE Trans. Med. Imaging 1998, 17, 357–361. [Google Scholar] [CrossRef] [PubMed]
  11. Alonso, R.; Villuendas, F.; Borja, J.; Barragán, L.; Salinas, I. Low-cost, digital lock-in module with external reference for coating glass transmission/reflection spectrophotometer. Meas. Sci. Technol. 2003, 14, 551–557. [Google Scholar] [CrossRef]
  12. Kempter, R.; Leibold, C.; Buzsáki, G.; Diba, K.; Schmidt, R. Quantifying circular associations: Hippocampal phase precession. J. Neurosci. Methods 2012, 207, 113–124. [Google Scholar] [CrossRef] [PubMed]
  13. Ghosh, K.; Jammalamadaka, S.R.; Vasudaven, M. Change-point problems for the von Mises distribution. J. Appl. Stat. 1999, 26, 423–434. [Google Scholar] [CrossRef]
  14. Sengupta, A.; Laha, A.K. A Bayesian analysis of the change-point problem for directional data. J. Appl. Stat. 2008, 35, 693–700. [Google Scholar] [CrossRef]
  15. Routtenberg, T.; Tabrikian, J. Bayesian Parameter Estimation Using Periodic Cost Functions. IEEE Trans. Signal Process. 2012, 60, 1229–1240. [Google Scholar] [CrossRef]
  16. Jammalamadaka, S.R.; SenGupta, A. Topics in Circular Statistics; Series on Multivariate Analysis; World Scientific: Singapore, 2001; Volume 5, 336p. [Google Scholar]
  17. Mardia, K.; Jupp, P. Directional Statistics; Wiley Series in Probability and Statistics; John Wiley and Sons Inc.: Chichester, UK, 1999. [Google Scholar]
  18. Stienne, G.; Reboul, S.; Azmani, M.; Choquel, J.B.; Benjelloun, M. GNSS dataless signal tracking with a delay semi-open loop and a phase open loop. Signal Process. 2013, 93, 1192–1209. [Google Scholar] [CrossRef]
  19. Stienne, G.; Reboul, S.; Azmani, M.; Choquel, J.; Benjelloun, M. A multi-temporal multi sensor circular fusion filter. Inform. Fusion 2014, 18, 86–100. [Google Scholar] [CrossRef]
  20. Marković, I.; Ćesić, J.; Petrović, I. Von Mises Mixture PHD Filter. IEEE Signal Process. Lett. 2015, 22, 2229–2233. [Google Scholar] [CrossRef]
  21. Kurz, G.; Gilitschenski, I.; Julier, S.; Hanebeck, U. Recursive Bingham filter for directional estimation involving 180 degree symmetry. J. Adv. Inform. Fusion 2014, 9, 90–105. [Google Scholar]
  22. Kurz, G.; Gilitschenski, I.; Hanebeck, U.D. Recursive Bayesian filtering in circular state spaces. IEEE Aerosp. Electron. Syst. Mag. 2016, 31, 70–87. [Google Scholar] [CrossRef]
  23. Hussin, A.G. Hypothesis Testing of Parameters for Ordinary Linear Circular Regression. Pak. J. Stat. Oper. Res. 2006, 2, 79–86. [Google Scholar] [CrossRef]
  24. Jona-Lasinio, G.; Gelfand, A.; Jona-Lasinio, M. Spatial analysis of wave direction using wrapped Gaussian processes. Ann. Appl. Stat. 2012, 6, 1478–1498. [Google Scholar] [CrossRef]
  25. Kucwaj, J.C.; Reboul, S.; Stienne, G.; Choquel, J.B.; Benjelloun, M. Circular Regression Applied to GNSS-R Phase Altimetry. Remote Sens. 2017, 9, 651. [Google Scholar] [CrossRef]
  26. Gould, A.L. A Regression Technique for Angular Variates. Biometrics 1969, 25, 683–700. [Google Scholar] [CrossRef] [PubMed]
  27. Johnson, R.A.; Wehrly, T.E. Some Angular-Linear Distributions and Related Regression Models. J. Am. Stat. Assoc. 1978, 73, 602–606. [Google Scholar] [CrossRef]
  28. Fisher, N.I.; Lee, A.J. Regression models for an angular response. Biometrics 1992, 48, 665–677. [Google Scholar] [CrossRef]
  29. Esmaieeli Sikaroudi, A.; Park, C. A Mixture of Linear-Linear Regression Models for Linear-Circular Regression. arXiv, 2016; arXiv:stat.ME/1611.00030. [Google Scholar]
  30. Miklós, A.; Hess, P.; Bozóki, Z. Application of acoustic resonators in photoacoustic trace gas analysis and metrology. Rev. Sci. Instrum. 2001, 72, 1937–1955. [Google Scholar] [CrossRef]
  31. Wang, G.; Yi, H.; Hubert, P.; Deguine, A.; Petitprez, D.; Maamary, R.; Fertein, E.; Rey, J.M.; Sigrist, M.W.; Chen, W. Filter-Free Measurements of Black Carbon Absorption Using Photoacoustic Spectroscopy. In Proceedings of the SPIE Quantum Sensing and Nano Electronics and Photonics XIV, San Francisco, CA, USA, 28 January–2 February 2017. [Google Scholar]
Figure 1. Principle of the determination of R. Abscissae represents the in phase component I = A cos ( ϕ k ) + η I , and ordinate represents the quadrature component Q = A sin ( ϕ k ) + η Q , where η I and η Q are Gaussian noise terms. R = O C = I 2 + Q 2 is an observation of A in conventional dual-phase LIA method. Let Δ I and Δ Q be the uncertainty in I and Q, caused by η I and η Q . ϕ 1 and ϕ 2 are the observed phase delays of the input noisy signal. The red area describes the uncertainties using traditional dual-phase LIA. In this article we consider the following observation R c = O B + B C = Q sin ( ϕ k ) + I cos ( ϕ k ) of A. In this case the area of uncertainties is a line between C 1 and C 2 . This area of uncertainties is smaller than in the classical case and it shows that the proposed observation R c is more accurate than R.
Figure 1. Principle of the determination of R. Abscissae represents the in phase component I = A cos ( ϕ k ) + η I , and ordinate represents the quadrature component Q = A sin ( ϕ k ) + η Q , where η I and η Q are Gaussian noise terms. R = O C = I 2 + Q 2 is an observation of A in conventional dual-phase LIA method. Let Δ I and Δ Q be the uncertainty in I and Q, caused by η I and η Q . ϕ 1 and ϕ 2 are the observed phase delays of the input noisy signal. The red area describes the uncertainties using traditional dual-phase LIA. In this article we consider the following observation R c = O B + B C = Q sin ( ϕ k ) + I cos ( ϕ k ) of A. In this case the area of uncertainties is a line between C 1 and C 2 . This area of uncertainties is smaller than in the classical case and it shows that the proposed observation R c is more accurate than R.
Sensors 17 02615 g001
Figure 2. (a) Basic dual-phase lock-in amplification (LIA); (b) Linear-Circular regression-based dual-phase LIA. ADC: analog-digital convertor. ∑: sum function working as low-pass filter. Arctan: the inverse tangent function.
Figure 2. (a) Basic dual-phase lock-in amplification (LIA); (b) Linear-Circular regression-based dual-phase LIA. ADC: analog-digital convertor. ∑: sum function working as low-pass filter. Arctan: the inverse tangent function.
Sensors 17 02615 g002
Figure 3. Mean estimate of A in the case of incoherent detection (ID).
Figure 3. Mean estimate of A in the case of incoherent detection (ID).
Sensors 17 02615 g003
Figure 4. (a) Estimates of the input signal amplitude A as a function of the signal-to-noise ratio (SNR). Assessment of coherent detection (CD) and incoherent detection (ID) methods for different T w : 0.001 s, 0.01 s, 0.1 s and 0.25 s; (b) Frequency response of the Low-Pass Filters (LPF) with T w of 0.001 s, 0.01 s, 0.1 s and 0.25 s.
Figure 4. (a) Estimates of the input signal amplitude A as a function of the signal-to-noise ratio (SNR). Assessment of coherent detection (CD) and incoherent detection (ID) methods for different T w : 0.001 s, 0.01 s, 0.1 s and 0.25 s; (b) Frequency response of the Low-Pass Filters (LPF) with T w of 0.001 s, 0.01 s, 0.1 s and 0.25 s.
Sensors 17 02615 g004
Figure 5. (a) Standard deviation of β ^ and α ^ for a SNR of 22 dB; (b) Example of LIA detection using coherent and incoherent detections.
Figure 5. (a) Standard deviation of β ^ and α ^ for a SNR of 22 dB; (b) Example of LIA detection using coherent and incoherent detections.
Sensors 17 02615 g005
Figure 6. Architectures of LIA involved in a PAS-based NO 2 sensor.
Figure 6. Architectures of LIA involved in a PAS-based NO 2 sensor.
Sensors 17 02615 g006
Figure 7. Experimental set up of photoacoustic spectroscopy (PAS)-based NO 2 sensor.
Figure 7. Experimental set up of photoacoustic spectroscopy (PAS)-based NO 2 sensor.
Sensors 17 02615 g007
Figure 8. (a) Estimated NO 2 concentrations using ID and CD approaches with regards to the NO 2 concentration measured by AC-31M analyzer. The error bars are the root mean square error between the AC-31M measurements and the estimates obtained with the ID or CD methods; (b) Observed and estimated phase delay for a NO 2 concentration of 76 ppbv; (c) Observed and estimated phase delay for a NO 2 concentration of 6 ppbv.
Figure 8. (a) Estimated NO 2 concentrations using ID and CD approaches with regards to the NO 2 concentration measured by AC-31M analyzer. The error bars are the root mean square error between the AC-31M measurements and the estimates obtained with the ID or CD methods; (b) Observed and estimated phase delay for a NO 2 concentration of 76 ppbv; (c) Observed and estimated phase delay for a NO 2 concentration of 6 ppbv.
Sensors 17 02615 g008
Figure 9. Performance comparisons of the estimated NO 2 concentration obtained with a SR830 lock-in amplifier and the proposed CD approach. The error bars are the standard deviation of the concentration estimates.
Figure 9. Performance comparisons of the estimated NO 2 concentration obtained with a SR830 lock-in amplifier and the proposed CD approach. The error bars are the standard deviation of the concentration estimates.
Sensors 17 02615 g009
Figure 10. Histogram showing the distribution of the estimated NO 2 concentration (of 24.5 ppbv measured with a reference NOx analyser) by CD (a) and SR830E (b).
Figure 10. Histogram showing the distribution of the estimated NO 2 concentration (of 24.5 ppbv measured with a reference NOx analyser) by CD (a) and SR830E (b).
Sensors 17 02615 g010
Table 1. Minimum SNR obtained for the estimation of the input signal amplitude with an uncertainty less than 10 % . The SNR is provided as a function of the observation window T o .
Table 1. Minimum SNR obtained for the estimation of the input signal amplitude with an uncertainty less than 10 % . The SNR is provided as a function of the observation window T o .
T o (s)0.250.501.001.502.002.503.003.504.00
SNR (dB)−14−19−22−23−24−24−24−24−24

Share and Cite

MDPI and ACS Style

Wang, G.; Reboul, S.; Choquel, J.-B.; Fertein, E.; Chen, W. Circular Regression in a Dual-Phase Lock-In Amplifier for Coherent Detection of Weak Signal. Sensors 2017, 17, 2615. https://doi.org/10.3390/s17112615

AMA Style

Wang G, Reboul S, Choquel J-B, Fertein E, Chen W. Circular Regression in a Dual-Phase Lock-In Amplifier for Coherent Detection of Weak Signal. Sensors. 2017; 17(11):2615. https://doi.org/10.3390/s17112615

Chicago/Turabian Style

Wang, Gaoxuan, Serge Reboul, Jean-Bernard Choquel, Eric Fertein, and Weidong Chen. 2017. "Circular Regression in a Dual-Phase Lock-In Amplifier for Coherent Detection of Weak Signal" Sensors 17, no. 11: 2615. https://doi.org/10.3390/s17112615

APA Style

Wang, G., Reboul, S., Choquel, J. -B., Fertein, E., & Chen, W. (2017). Circular Regression in a Dual-Phase Lock-In Amplifier for Coherent Detection of Weak Signal. Sensors, 17(11), 2615. https://doi.org/10.3390/s17112615

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop