Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Dual-optimized adaptive Kalman filtering algorithm based on BP neural network and variance compensation for laser absorption spectroscopy

Open Access Open Access

Abstract

A dual-optimized adaptive Kalman filtering (DO-AKF) algorithm based on back propagation (BP) neural network and variance compensation was developed for high-sensitivity trace gas detection in laser spectroscopy. The BP neural network was used to optimize the Kalman filter (KF) parameters. Variance compensation was introduced to track the state of the system and to eliminate the variations in the parameters of dynamic systems. The proposed DO-AKF algorithm showed the best performance compared with the traditional multi-signal average, extended KF, unscented KF, KF optimized by BP neural network (BP-KF) and KF optimized by variance compensation (VC-KF). The optimized DO-AKF algorithm was applied to a QCL-based gas sensor system for an exhaled CO analysis. The experimental results revealed a sensitivity enhancement factor of 23. The proposed algorithm can be widely used in the fields of environmental pollutant monitoring, industrial process control, and breath gas diagnosis.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

This work focuses on filtering techniques for laser spectroscopy, which has been extensively used in atmospheric environmental monitoring, industrial process control, and medical diagnosis [13]. Exhaled breath analysis is an emerging field that offers a non-invasive, painless, and online means for medical diagnosis with a fast response time. However, the sensitivity and resolution of measurements are affected by intrinsic noises and interferences induced by electrical components and background changes [45]. Therefore, a highly precise processing of the measured spectra data is crucial.

Several filtering techniques based on hardware have been recently developed to process the measured spectra [67]. However, their practical application has been limited by their complex system structure, high cost, and poor compatibility. Moreover, the intrinsic noises, including flicker, shot, and detector noises, are difficult to remove [89]. Therefore, a filtering technique based on software has become an ideal choice due to its simple system structure and low cost. Although numerous filtering techniques based on software have been applied to process gas absorption spectra, only few techniques have been considered effective for the online filtering of these spectra. Although widely used, the multi-signal averaging filtering (MAF) is a time-consuming method that can only suppress white noise [1011]. Meanwhile, the wavelet transform denoising technique allows for an efficient online filtering [1213], but depends on a large number of parameters and involves a complex denoising calculation. The discontinuity of wavelet coefficients near the threshold also limits the filtering performance of this approach [1415]. An adaptive Savitzky–Golay filter algorithm has been recently proposed for real-time infrared gas detection and has demonstrated an excellent performance in optimizing the filtering process. However, this approach can only filter out high-frequency noise [1617].

Kalman filtering (KF), as an efficient tool for reducing the noise of polluted signals, has recently attracted the interest of researchers and has been successfully applied in the fields of GPS positioning and navigation, satellite measurement and control, biomedicine, and signal processing [1820]. KF does not need to store a large number of measured data and demonstrates a fast processing speed. Therefore, KF is a suitable filtering technique for real-time exhaled breath detection [2122]. However, in a nonlinear system, KF causes a serious deviation in the spectral absorption part. To address this problem, an extended KF (EKF) algorithm has been proposed. This approach can achieve an optimal prediction and control of the nonlinear model, thereby ensuring its applicability [2324]. Although this approach obtains a good effect in nonlinear systems, EKF has an uncertain system model and filtering parameters for a fuzzy system. Recently, an adaptive KF (AKF) algorithm that combines system identification and filtering estimation has been proposed. This approach continuously evaluates the dynamics of the system and then estimates and modifies the model parameters and noise statistical characteristics to reduce the actual error of filtering [2526]. AKF shows an excellent performance in dynamic systems and has been applied to measure non-stationary signals, such as instantaneous power and time-varying spectra [2728]. The filtering deviation of the nonlinear system in AKF is usually corrected by selecting the sub-optimal filtering results that reduce the efficiency of the algorithm.

In this paper, a dual-optimized adaptive KF (DO-AKF) algorithm is developed based on back propagation (BP) neural network and variance compensation to enhance the measurement system performance of KF. The BP neural network is used to optimize the KF parameters, and variance compensation is applied to obtain an accurate system state equation and filtering result. Compared with KF, DO-AKF can realize “self-adjustment” and “tracking feedback” to achieve optimal signal-to-noise ratios (SNRs) of the spectra. The construction of DO-AKF is discussed in section 2. The simulation of the absorption spectrum filtering process and the experimental setup of the exhaled CO gas detection are described in section 3. The denoising results of several spectral filtering algorithms are compared in section 4, and the proposed DO-AKF demonstrates the best performance.

2. Mathematical model and algorithm theory

KF is an optimal state estimation algorithm that is widely used in linear systems. However, this approach may cause a serious deviation in nonlinear signal processing. The system adaptability can be improved by introducing the BP neural network.

The BP neural network has a hierarchical feed forward network architecture. The output of each layer is sent directly to each neuron in the above layer. The entire process requires a minimum of three layers, namely, the input layer, hidden layer, and output layer [29]. The structure of the BP neural network is shown in Fig. 1.

 figure: Fig. 1.

Fig. 1. Schematic diagram of the BP neural network structure and information flow.

Download Full Size | PDF

The error and correlation of the dynamic system may influence the stability of KF. To improve filtering stability and accuracy, the filtering variance matrix needs to be optimized by using variance compensation.

Variance compensation is coupled with the spectral vector and prediction residual [30]. The spectral noise variance matrix is estimated in real time. To compensate for the limitations in spectral noise variance and covariance estimation, the spectral filtering vector needs to be predicted and corrected continuously. The proposed DO-AKF is constructed by combining KF theory with the BP neural network and variance compensation algorithm, for which flow chart is shown in Fig. 2.

 figure: Fig. 2.

Fig. 2. Flow chart of the developed DO-AKF algorithm.

Download Full Size | PDF

2.1. Spectral system process model

In laser spectroscopy, the measurement equation and the state equation in the state space model can be characterized by the spectral state and measurement vector in the spectral filtering system. The spectral filtering system and measurement vector are defined as

$$X(k) = AX(k - 1) + BU(k) + W(k)\,\mbox{and}$$
$$Z(k) = HX(k) + V(k),$$
where k is the discrete time index, $X(k)$ is the spectral state vector, $U(k)$ is the process control vector of the system, $Z(k)$ is the measurement vector. A, $B$, and $H$ are the spectral state transition matrixes. $W(k)$, $V(k)$ are the process and measurement noises, respectively. All of these parameters A, $B,\, H,\, W(k)$ and $V(k)$ are the hyperparameters of the state space model in the meanwhile.

The spectral state vector $X(k)$ and measurement vector $Z(k)$ can be substituted by the estimated and detected values of spectral data, respectively. As for the hyperparameters in the state space model, it is generally assumed that the parameters of the state space model are known and constants for the convenience of calculation. However in the proposed spectral filtering system, BP neural network is structured to traverse the optimal values of these hyperparameters or the state transition matrixes A, $B$, and H instead of defining them as constants as usual. The process and measurement noises are represented by $W(k)$ and $V(k)$, where the expected value, $E[{W(k)} ]= E[{V(k)} ]= 0$. And their covariance matrices are represented by $Q(k)$ and $R(k)$, respectively. The covariance $Q(k)$ represents the process noise introduced by the actual optical path, while $R(k)$ represents the measurement noise introduced by the sensor.

The spectral vectors are not constant in the actual system. Therefore, the KF algorithm cannot accurately reflect the time-varying characteristics of the spectral model. The BP neural network is used to optimize the real-time system states and vectors, and each matrix is updated in each step of KF.

1) Time recursive state variables

According to the system model, the predicted state vector can be obtained recursively as

$$X(k|k - 1) = AX(k - 1|k - 1) + BU(k),$$
where $X(k|k - 1)$ is the predicted spectral state vector, and $X(k - 1|k - 1)$ is the previous state vector.

2) Forward calculation of error covariance

The corresponding predicted covariance matrix P is obtained as

$$P(k|k - 1) = AP(k - 1|k - 1){A^T} + Q(k - 1),$$
where ${A^T}$ means the transpose operation of A. $P(k|k - 1)$ is the corresponding predicted covariance matrix, and $P(k - 1|k - 1)$ is the previous covariance matrix.

2.2. Spectral parameters update equation of KF

  • 1) Update the spectral filtering state variables
    $$X(k|k) = X(k|k - 1) + Kg(k)({Z(k) - HX(k|k - 1)} )$$
  • 2) Calculate the Kalman gain
    $$Kg(k) = P(k|k - 1){H^T}{(HP(k|k - 1){H^T} + R)^{ - 1}}$$
  • 3) Update the error covariance
    $$P(k|k) = ({I - Kg(k)H} )P(k|k - 1),$$
    where $X(k|k)$ and $P(k|k)$ are the optimized state vector and error covariance matrix, respectively, and $Kg$ is the Kalman gain.

The original absorption spectra includes two parts, i.e. the non-absorption baseline part (corresponding to laser intensity background) and the molecular absorption peak part. In this case, the non-linear laser intensity background makes large errors and even diverging problems when a common KF algorithm model is used. Therefore, BP neural network is innovatively proposed to traverse their optimal values, instead of defining them as constants in a common KF algorithm model. The variance compensation is applied to compensate for the inadequacy of dynamic noise equation and covariance estimation. Combining with the precise state transition parameters and the variance compensation, the improved KF algorithm model can adapt to the instability of the nonlinear system (i.e. the original spectroscopy signal).

2.3. BP neural network and variance compensation optimization

Based on the variance compensation principle, the prediction residual $F(k)$ is defined as

$$F(k) = L(k) - \hat{L}(k|k - 1),$$
where $L(k)$ and $\hat{L}(k|k - 1)$ denote the measured and best predicted vectors, respectively. The prediction residuals are used to estimate and correct the covariance vector of dynamic noise and to compensate for the deficiencies of the dynamic noise equation or covariance estimation in the filtering process.

The BP neural network is constructed through the spectral state vector, and the following parameters are selected for the input layer of the BP neural network: (a). the difference between the spectral state and predicted state vectors ($X(k) - X(k|k - 1)$), (b). the gap between the measured vector and forecast measured data ($Z(k) - HX(k|k - 1)$), and (c). KF gain filtering $Kg$. After optimizing the BP neural network, the filtering compensation vector $\Delta X(k)$ is defined as

$$\Delta X(k) = D - X(k),$$
where D is the theoretical vector after filtering. By optimizing the BP neural network and variance compensation, DO-AKF can be used to improve the spectral filtering accuracy.

3. Experimental details

The spectral data is usually polluted by various noises, which are categorized into two sources: interference noises due to spectral and optical interferences, and electronic noises from detectors, lasers, and additional electronics. Therefore, it is significant to perform data preprocessing to obtain high SNR data. In this study, a dual-optimized adaptive DO-AKF algorithm based on BP neural network and variance compensation was developed for high-sensitivity trace gas detection. To demonstrate the applicability for this technique, the simulated spectra by referring to the HITRAN database [31] and experimental CO data were used to evaluate the proposed algorithm, and detailed comparison with other common algorithms are also presented.

3.1 Optimization of the BP neural network

Similar to other digital signal processing techniques, the effectiveness of KF largely depends on the system state equation. The selection of an appropriate system state equation is essential to achieve a trade-off between reducing noise and avoiding bias. In this work, a reasonable BP neural network is constructed to optimize the system state equation and correct the filtering result.

The number of neurons in the input and output layers of the BP neural network is determined by the dimensional data of the input and output vectors. Given that the input and output data are composed of wave number and transmittance, the number of neurons in the input and output layers is 2. The BP neural network with a single hidden layer can theoretically approximate any continuous function in a closed interval. The number of hidden layer units can be determined by $P\mbox{ = }\sqrt {N + M} \mbox{ + }a$, where $N$ and $M$ denote the number of input and output layer units, respectively [32]. However, this equation only provides a rough estimation, and both convergence and learning efficiency need to be improved. In practical applications, the optimization of the hidden layer number and number of hidden layer units is crucial to achieve the best SNR. The simulated spectra shown in section 4.1 are used for constructing the optimal network structure to find the optimal SNR and optimizing the filter. The number of hidden layers and hidden layer units ranges between 1 and 10, and the filtering results of different hidden layer and unit numbers are compared. The SNRs under different layer and unit numbers are shown in Fig. 3, and the results for the optimal number of hidden layer units under each hidden layer are shown in Table 1. Figure 3 shows that the effectiveness of final denoising largely depends on the hidden layers and hidden layer units. The optimal effect is achieved under the combination of four hidden layers and seven hidden layer units.

Tables Icon

Table 1. Optimal hidden layer unit number and filtered SNR statistics for different hidden layer numbers (ranging from 1 to 10).

 figure: Fig. 3.

Fig. 3. SNR distribution of filtered signals with different combinations of hidden layer numbers and hidden layer units (ranging from 1 to 10).

Download Full Size | PDF

3.2 Experimental setup

To demonstrate the performance of DO-AKF, this algorithm is applied to process the absorption spectroscopy of the exhaled CO. The experimental setup is illustrated in Fig. 4. A thermoelectrically cooled CW room temperature QCL (ALPES, Switzerland) emitting at ∼4.6 µm is used as the optical source with a maximum peak output power of 30 mW. The laser output power and wave number are controlled by temperature (Starter Kit, Alpes) and current controllers (LDX 3232, ILX Light wave), respectively. The linewidth of laser is estimated to be less than 0.001 cm−1, and the beam of the QCL laser passes through a CaF2 mirror to combine the tracer laser with the trace laser. The astigmatic Herriott multi-pass cell with an effective optical path length of 62.5 m is used as gas cell. A ZnSe beam splitter is used to co-align a visible red light (632.8 nm) for beam tracing purposes. The exiting beam is then focused into a thermoelectric cooled mercury cadmium telluride detector (Vigo, VI-4TE-5) with a 75 mm focal length plano-convex lens. This detector does not require liquid nitrogen cooling and can be used for long-term operations. The laser frequency is scanned from 2188.5 cm−1 to 2191.5 cm−1 by using a triangular wave at a typical frequency of 100 Hz. The detector signal is subsequently acquired by using a DAQ card (National Instruments, USB-6259), and the data are processed and saved by a LabVIEW program. A polytetrafluoroethylene filter (0.22 µm) for dust retention and an integrating mass flow controller (Alicat, USA) are connected to the inlet of the multi-pass cell. The system pressure is maintained by using a diaphragm pump and pressure controller to avoid potential upstream pressure transients.

 figure: Fig. 4.

Fig. 4. Schematic diagram of the experimental device for exhaled breath detection.

Download Full Size | PDF

4. Analysis and results

4.1 Results and analysis of simulated data

In the range of 2188.5 cm−1 to 2191.5 cm−1 and pressure of 1 atm, the spectrum of CO satisfying the Viogt line-shape function at a temperature of 296 K is obtained by using the HITRAN database. An optical path length of 62.5 m is used and the CO mixing ratio is 1.4 ppm. The spectral data are superimposed with the Gaussian noise (white noise) and non-Gaussian noise, respectively. Under the background of Gaussian noise (white noise, the mean value is 0, the variance is 1 and the standard deviation is 0.009), the SNR of the spectral data is 7.9162 dB. Noise pollution greatly reduces the accuracy of the spectrum. The MAF, EKF, unscented KF (UKF), variance compensation optimized KF (VC-KF), BP neural network optimized KF (BP-KF), and DO-AKF are used to process the simulated spectrum. The filtered results are shown in Fig. 5, and the SNRs of different filtering methods are shown in Table 2.

 figure: Fig. 5.

Fig. 5. Simulation of the CO gas absorption spectra of the six filtering algorithms (under the background of Gaussian noise).

Download Full Size | PDF

Tables Icon

Table 2. Comparison of the SNR and RMSE for the filtered signal obtained by different filtering methods.

As shown in Fig. 5 and Table 2, the SNRs have increased from 7.9162 dB to 20.2685 dB, and the root mean square error (RMSE) has decreased from 5.255×10−5 to 2.401×10−6. The filtering effect of MAF is not obvious, whereas the spectra after the application of EKF, UKF, VC-KF, BP-KF, and DO-AKF are relatively smooth. The VC-KF shows a more effective filtering effect than EKF, and UKF has higher filtering SNR than BP-KF. The DO-AKF shows the best performance and obtains the highest SNR among the six filtering algorithms.

The frequency response of the filter is an important property of filtering method. The residual calculated from the difference between the filtered data and the origin data without noise is used to evaluate the performance of the proposed filter algorithm. The frequency domain distribution of the residual noise is analyzed by Fast Fourier Transformation (FFT) before and after applying various filters, as shown in Fig. 6. As can been seen, MAF and EKF have poor noise filtering performance, which still retain the frequency distribution characteristics of the original noise. UKF, VC-KF, BP-KF, DO-AKF can remove more high frequency noise. DO-AKF shows the best filtering performance, and the amplitude of low-frequency noise can be reduced to one third compared to UKF, VC-KF and BP-KF algorithm.

 figure: Fig. 6.

Fig. 6. Frequency domain distribution of raw noise and filtered residual noise.

Download Full Size | PDF

Indeed, the simulated noised in Fig. 5 is Gaussian noise (white noise). In order to demonstrate the influence of the non-Gaussian noise on the proposed DO-AKF algorithm. The non-Gaussian background noise (exponential distribution noise, mean parameter = 0.005) is simulated for test, and the filtered results are shown in Fig. 7, and the SNRs of different filtering methods are summarized in Table 3.

 figure: Fig. 7.

Fig. 7. Simulation of the CO gas absorption spectra of the six filtering algorithms. (under the background of non-Gaussian noise).

Download Full Size | PDF

Tables Icon

Table 3. Comparison of the SNR and RMSE for the filtered signal obtained by different filtering methods.

As shown in Fig. 7 and Table 3, in the case of non-Gaussian noise, a spectral data with a SNR of 8.2634 dB is simulated. After the application of DO-AKF filtering algorithm, the SNR has been improved from 8.2634 dB to 21.1207 dB, and the RMSE has decreased from 4.913×10−5 to 2.286×10−6. Although UKF and EKF are both non-linear filtering methods, UKF has better filtering effect (the second highest filtering SNR) than EKF. And DO-AKF also shows the best performance and obtains the highest SNR under these non-Gaussian noise among the six filtering algorithms.

The transmission spectrum and residual of EKF and VC-KF are compared in Fig. 8. VC-KF obtains a lower residual than EKF, whereas the SNR significantly increases from 14.2457 dB (EKF) to 16.1512 dB (VC-KF). Variance compensation is used to optimize the filtering process in real time. The variance is more flexible in noise reduction and can compensate for the filtering variance based on the current filtering state.

 figure: Fig. 8.

Fig. 8. Comparison of transmission spectra and residual distributions for the EKF and VC-KF algorithms. (a) EKF and residual of the CO gas absorption spectra; and (b) VC-KF and residual of the CO gas absorption spectra

Download Full Size | PDF

Given the aforementioned factors, DO-AKF is identified as the most suitable filtering algorithm for extracting a gas absorption signal from a contaminated signal. Depending on the optimal traversal of BP neural network to the filtering parameters and the optimal compensation of variance to the filtering covariance, DO-AKF can be adapted to both Gaussian and non-Gaussian noise and maintain the optimal filtering performance.

4.2 Results and analysis of experimental data

To demonstrate the effectiveness of DO-AKF in processing the spectral signal, the measured and filtered spectra of exhaled CO at a pressure of 1atm and temperature of 296 K are shown in Fig. 9(a). In our case, measurement time resolution of 1 Hz is commonly selected for practical application, which including data acquisition, signal average, filtering process and spectral fitting, data storage. We restricted the spectral data records to 1024 sample points, sampling rate of data acquisition board is set to 200 kHz. Higher responsivity can also be realized by using high-speed data acquisition card. The absorption characteristic is partially submerged in noise. The filtered transmission spectrum has been significantly smoothed, thereby facilitating the identification and fitting of the absorption spectrum. The measured transmission spectrum of exhaled CO is converted into a gas absorption spectrum and fitted by using the Voigt function as shown in Fig. 9(b). The Doppler broadened HWHM ωG is a function of temperature T and molecular weight M, which can be written as:

$${\omega _G} = \frac{{{v_0}}}{c}\sqrt {\frac{{2k\mbox{T}\ln 2}}{M}} = 3.5812 \times {10^{ - 7}}{v_0}\sqrt {\frac{T}{M}} .$$
The Doppler width is fixed during the fitting procedure. The peak position is fixed in the wavelength obtained from the database. The peak position is fixed in the wavelength obtained from the database. The fitted result shows that the R-square value of the fitting is as high as 0.99992, thereby suggesting that the filtered absorption spectrum shows good agreement with the theoretical shape. The corresponding enhancement factors (the ratio of the peak absorption coefficient and standard deviation of the fitted spectrum) is 23 times higher than the fitted result without a filter, thereby highlighting the great potential of DO-AKF in trace gas detection in various fields, including environmental pollution, biomedical, and industrial testing.

 figure: Fig. 9.

Fig. 9. (a) Typical example of experimentally determined and filtered CO absorption spectra at a pressure of 1 atm; and (b) absorption spectra of CO and fitted spectra with a Voigt function centered at ∼2190.017 cm−1.

Download Full Size | PDF

The bulk matrix of breath is a mixture of N2, O2, CO2, H2O and other compounds, such as CO and N2O. The main constitutes of human breath, N2 and O2, usually do not interfere with the spectra of CO. Therefore, the predominant interference specie for breath CO detection is H2O and CO2. The line intensity of CO is several orders of magnitude higher than H2O and CO2. Generally, by reducing sample gas pressure, the spectral interference of H2O and other species can be effectively reduced. To verify the universality of the DO-AKF algorithm for spectral data at different pressures, a series of spectra are measured at different pressures (30, 50, 80, and 100 mbar) near 2190.01 cm−1, and the filtered results of the measured spectral data after DO-AKF filtering process are shown in Fig. 10. At the same time, the top-left panels and top-right panels show the filtering details near absorption peaks. The algorithm can effectively eliminate the interference of background noise.

 figure: Fig. 10.

Fig. 10. The filtered results of the real-time measured exhaled CO gas spectral data (at four different pressures) after the filtering process of DO-AKF.

Download Full Size | PDF

Allan variance is used to investigate the performance of the sensor system. The custom CO gas is injected into the sample cell for a time scale of several minutes, and a pressure of 150 mbar is maintained. The CO concentrations are calculated from the fitted peak area of the acquired spectra. The real-time measured CO concentration levels are shown in Fig. 11(a), and the Allan deviation is shown in Fig. 11(b). The average and standard deviations of the measured data are 457.45 ppb and 18.20 ppb, respectively, while the average and standard deviations of the filtered data are 457.27 ppb and 0.86 ppb, respectively. As shown in Fig. 11(b), the optimum sensitivity of 0.81 ppbv is reached at 181 s.

 figure: Fig. 11.

Fig. 11. (a) Real-time measured CO concentrations; and (b) the Allan deviation plot of the measured and filtered data.

Download Full Size | PDF

5. Conclusion

In this paper, the DO-AKF algorithm was developed based on BP neural network and variance compensation to achieve a high-sensitivity trace gas detection in laser spectroscopy. The BP neural network was used to optimize the KF parameters, and variance compensation was introduced to eliminate the parameter variation and divergence of dynamic systems. The developed algorithm was then applied to a QCL based gas sensor system to measure exhaled CO, and a sensitivity enhancement factor of 23 was obtained. DO-AKF showed two main advantages over the other algorithms, including MAF, EKF, UKF, VC-KF, and BP-KF. First, DO-AKF showed a significant filtering effect. Second, in dynamic systems, DO-AKF demonstrated universality by using variance compensation and BP neural network. Therefore, DO-AKF can be widely used in the fields of environmental pollutant monitoring, industrial process control, and breath gas diagnosis.

Funding

University Natural Science Research Project of Anhui Province (KJ2018A0034); National Basic Research Program of China (973 Program) (2016YFC03022022016YFC0301900); National Natural Science Foundation of China (41875158, 61705002, 61905001); Natural Science Foundation of Anhui Province (1808085QF198, 1908085QF276).

References

1. Y. F. Ma, Y. Tong, Y. He, X. G. Jin, and F. K. Tittel, “Compact and sensitive mid-infrared all-fiber quartz-enhanced photoacoustic spectroscopy sensor for carbon monoxide detection,” Opt. Express 27(6), 9302–9312 (2019). [CrossRef]  

2. H. P. Wu, L. Dong, and H. D. Zheng, “Beat frequency quartz-enhanced photoacoustic spectroscopy for fast and calibration-free continuous trace-gas monitoring,” Nat. Commun. 8(1), 15331 (2017). [CrossRef]  

3. G. P. Asner, R. E. Martin, D. E. Knapp, R. Tupayachi, C. B. Anderson, F. Sinca, N. R. Vaughn, and W. Llactayo, “Airborne laser-guided imaging spectroscopy to map forest trait diversity and guide conservation,” Science 355(6323), 385–389 (2017). [CrossRef]  

4. P. Liu, W. Huang, W. Zhang, and L. Fang, “An emd-sg algorithm for spectral noise reduction of fbg-fp static strain sensor,” IEEE Photonics Technol. Lett. 29(10), 814–817 (2017). [CrossRef]  

5. A. Janani and M. Sasikala, “Investigation of different approaches for noise reduction in functional near-infrared spectroscopy signals for brain–computer interface applications,” Neural Comput. Appl. 28(10), 2889–2903 (2017). [CrossRef]  

6. A. Chabuda, P. Durka, and J. Żygierewicz, “High frequency SSVEP-BCI with hardware stimuli control and phase-synchronized comb filter,” IEEE T. Neur. Sys. Reh 26(2), 344–352 (2018). [CrossRef]  

7. J. Li, B. Yu, W. Zhao, and W. Chen, “A review of signal enhancement and noise reduction techniques for tunable diode laser absorption spectroscopy,” Appl. Spectrosc. Rev. 49(8), 666–691 (2014). [CrossRef]  

8. L. Zhang, G. Tian, J. Li, and B. Yu, “Applications of absorption spectroscopy using quantum cascade lasers,” Appl. Spectrosc. 68(10), 1095–1107 (2014). [CrossRef]  

9. M. Chen, C. Wang, J. Wang, H. Luo, and Z. Meng, “53-dB phase noise suppression and hz-range linewidth emission in compact brillouin/erbium fiber laser,” Opt. Express 25(16), 19216 (2017). [CrossRef]  

10. C. Li, L. Shao, L. Jiang, X. Qiu, J. Wei, and W. Ma, “Simultaneous measurements of CO and CO2 employing wavelength modulation spectroscopy using a signal averaging technique at 1.578 µm,” Appl. Spectrosc. 72(9), 1380–1387 (2018). [CrossRef]  

11. C. Eichner, S. F. Cauley, J. Cohen-Adad, H. E. Möller, R. Turner, K. Setsompop, and L. L. Wald, “Real diffusion-weighted MRI enabling true signal averaging and increased diffusion contrast,” NeuroImage, S1053811915006941 (2015). [CrossRef]  

12. Z. Bo, L. Sun, H. Yu, X. Yong, and Z. Cong, “A method for improving wavelet threshold denoising in laser-induced breakdown spectroscopy,” Spectrochim. Acta, Part B 107, 32–44 (2015). [CrossRef]  

13. R. X. Yi, L. B. Guo, X. H. Zou, J. M. Li, Z. Q. Hao, X. Y. Yang, X. Y. Li, X. Y. Zeng, and Y. F. Lu, “Background removal in soil analysis using laser- induced breakdown spectroscopy combined with standard addition method,” Opt. Express 24(3), 2607 (2016). [CrossRef]  

14. G. Wang, S. J. Kim, G. S. Kil, and S. W. Kim, “Optimization of wavelet and thresholding for partial discharge detection under hvdc,” IEEE Trans. Dielectr. Electr. Insul. 24(1), 200–208 (2017). [CrossRef]  

15. Z. Zainuddin and P. Ong, “Optimization of wavelet neural networks with the firefly algorithm for approximation problems,” Neural Comput. Appl. 28(7), 1715–1728 (2017). [CrossRef]  

16. B. Zimmermann and A. Kohler, “Optimizing savitzky-golay parameters for improving spectral resolution and quantification in infrared spectroscopy,” Appl. Spectrosc. 67(8), 892–902 (2013). [CrossRef]  

17. N. R. Koluguri, G. N. Meenakshi, and P. K. Ghosh, “Spectrogram enhancement using multiple window savitzky-golay (MWSG) filter for robust bird sound detection,” IEEE/ACM T. Audio Spe. 25(6), 1183–1192 (2017). [CrossRef]  

18. Y. V. Kistenev, A. V. Borisov, D. A. Kuzmin, O. V. Penkova, N. Y. Kostyukova, and A. A. Karapuzikov, “Exhaled air analysis using wideband wave number tuning range infrared laser photoacoustic spectroscopy,” J. Biomed. Opt. 22(1), 017002 (2017). [CrossRef]  

19. F. Caron, E. Duflos, D. Pomorski, and P. Vanheeghe, “GPS/IMU data fusion using multisensor kalman filtering: introduction of contextual aspects,” Inform. Fusion 7(2), 221–230 (2006). [CrossRef]  

20. H. Zhang, R. Ayoub, and S. Sundaram, “Sensor selection for kalman filtering of linear dynamical systems: complexity, limitations and greedy algorithms,” Automatica 78, 202–210 (2017). [CrossRef]  

21. Y. Huang, Y. Zhang, and X. Wang, “Kalman-filtering-based in-motion coarse alignment for odometer-aided sins,” IEEE Trans. Instrum. Meas. 66(12), 3364–3377 (2017). [CrossRef]  

22. D. G. Abrecht, J. M. Schwantes, R. K. Kukkadapu, B. S. Mcdonald, G. C. Eiden, and L. E. Sweet, “Real-time noise reduction for mössbauer spectroscopy through online implementation of a modified kalman filter,” Nucl. Instrum. Methods Phys. Res., Sect. A 773, 66–71 (2015). [CrossRef]  

23. Z. Yin, G. Li, Y. Zhang, and L. Jing, “Symmetric-strong-tracking-extended-kalman-filter-based sensorless control of induction motor drives for modeling error reduction,” IEEE T. Inform 15(2), 650–662 (2019). [CrossRef]  

24. K. Kumar, D. Yadav, and B. V. Srinivas, “Adaptive noise models for extended kalman filter,” J. Guid. Control & Dynam. 14(2), 475–477 (1991). [CrossRef]  

25. Y. Huang, Y. Zhang, Z. Wu, L. Ning, and J. Chambers, “A novel adaptive kalman filter with inaccurate process and measurement noise covariance matrices,” IEEE Trans. Autom. Control 63(2), 594–601 (2018). [CrossRef]  

26. L. Zanni, J. Y. L. Boudec, R. Cherkaoui, and M. Paolone, “A prediction-error covariance estimator for adaptive kalman filtering in step-varying processes: application to power-system state estimation,” IEEE Trans. Contr. Syst. Technol. 25(5), 1683–1697 (2017). [CrossRef]  

27. X. Fan, J. Wang, D. Zhu, and T. Qi, “Speckle noise reduction technique for lidar echo signal based on self-adaptive pulse-matching independent component analysis,” Opt. & Laser. Eng. 103, 92–99 (2018). [CrossRef]  

28. A. O. Salam, R. E. Sheriff, S. R. Al-Araji, K. Mezher, and Q. Nasir, “Adaptive interacting multiple model-kalman filter for multitaper spectrum sensing in cognitive radio,” Electron. Lett. 54(5), 321–322 (2018). [CrossRef]  

29. Z. Li and X. Zhao, “Bp artificial neural network based wave front correction for sensor-less free space optics communication,” Opt. Commun. 385, 219–228 (2017). [CrossRef]  

30. T. Nose, “Efficient implementation of global variance compensation for parametric speech synthesis,” IEEE/ACM T. Audio Spe. 24(10), 1694–1704 (2016). [CrossRef]  

31. L. S. Rothman, I. E. Gordon, Y. Babikov, A. Barbe, D. Chris Benner, F. Bernath, M. Birk, L. Bizzocchi, V. Boudon, L. R. Brown, A. Campargue, K. Chance, E. A. Cohen, L. H. Coudert, V. M. Devi, B. J. Drouin, A. Fayt, J.-M. Flaud, R. R. Gamache, J. J. Harrison, J.-M. Hartmann, C. Hill, J. T. Hodges, D. Jacquemart, A. Jolly, J. Lamouroux, R. J. Le Roy, G. Li, D. A. Long, O. M. Lyulin, C. J. Mackie, S. T. Massie, S. Mikhailenko, H. S. P. Müller, O. V. Naumenko, A. V. Nikitin, J. Orphal, V. Perevalov, A. Perrin, E. R. Polovtseva, C. Richard, M. A. H. Smith, E. Starikova, K. Sung, S. Tashkun, J. Tennyson, G. C. Toon, V. l. G. Tyuterev, and G. Wagner, “The HITRAN, Molecular spectroscopic database,” J. Quant. Spectrosc. Radiat. Transfer 130, 4–50 (2013). [CrossRef]  

32. N. J. Guliyev and V. E. Ismailov, “A single hidden layer feedforward network with only one neuron in the hidden layer can approximate any univariate function,” Neural Comput. 28(7), 1289–1304 (2016). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (11)

Fig. 1.
Fig. 1. Schematic diagram of the BP neural network structure and information flow.
Fig. 2.
Fig. 2. Flow chart of the developed DO-AKF algorithm.
Fig. 3.
Fig. 3. SNR distribution of filtered signals with different combinations of hidden layer numbers and hidden layer units (ranging from 1 to 10).
Fig. 4.
Fig. 4. Schematic diagram of the experimental device for exhaled breath detection.
Fig. 5.
Fig. 5. Simulation of the CO gas absorption spectra of the six filtering algorithms (under the background of Gaussian noise).
Fig. 6.
Fig. 6. Frequency domain distribution of raw noise and filtered residual noise.
Fig. 7.
Fig. 7. Simulation of the CO gas absorption spectra of the six filtering algorithms. (under the background of non-Gaussian noise).
Fig. 8.
Fig. 8. Comparison of transmission spectra and residual distributions for the EKF and VC-KF algorithms. (a) EKF and residual of the CO gas absorption spectra; and (b) VC-KF and residual of the CO gas absorption spectra
Fig. 9.
Fig. 9. (a) Typical example of experimentally determined and filtered CO absorption spectra at a pressure of 1 atm; and (b) absorption spectra of CO and fitted spectra with a Voigt function centered at ∼2190.017 cm−1.
Fig. 10.
Fig. 10. The filtered results of the real-time measured exhaled CO gas spectral data (at four different pressures) after the filtering process of DO-AKF.
Fig. 11.
Fig. 11. (a) Real-time measured CO concentrations; and (b) the Allan deviation plot of the measured and filtered data.

Tables (3)

Tables Icon

Table 1. Optimal hidden layer unit number and filtered SNR statistics for different hidden layer numbers (ranging from 1 to 10).

Tables Icon

Table 2. Comparison of the SNR and RMSE for the filtered signal obtained by different filtering methods.

Tables Icon

Table 3. Comparison of the SNR and RMSE for the filtered signal obtained by different filtering methods.

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

X ( k ) = A X ( k 1 ) + B U ( k ) + W ( k ) and
Z ( k ) = H X ( k ) + V ( k ) ,
X ( k | k 1 ) = A X ( k 1 | k 1 ) + B U ( k ) ,
P ( k | k 1 ) = A P ( k 1 | k 1 ) A T + Q ( k 1 ) ,
X ( k | k ) = X ( k | k 1 ) + K g ( k ) ( Z ( k ) H X ( k | k 1 ) )
K g ( k ) = P ( k | k 1 ) H T ( H P ( k | k 1 ) H T + R ) 1
P ( k | k ) = ( I K g ( k ) H ) P ( k | k 1 ) ,
F ( k ) = L ( k ) L ^ ( k | k 1 ) ,
Δ X ( k ) = D X ( k ) ,
ω G = v 0 c 2 k T ln 2 M = 3.5812 × 10 7 v 0 T M .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.