Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 286590, 9 pages
http://dx.doi.org/10.1155/2015/286590

Research Article

Compressive Sensing Approach in the Hermite Transform Domain

Srdjan Stanković, Ljubiša Stanković, and Irena Orović

University of Montenegro, Faculty of Electrical Engineering, Dzordza Vasingtona bb, 81000 Podgorica, Montenegro

Received 20 August 2015; Accepted 16 November 2015

Academic Editor: Yuqiang  Wu

Copyright © 2015 Srdjan Stanković et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Compressive sensing has attracted significant interest of researchers providing an alternative way to sample and reconstruct the signals. This approach allows us to recover the entire signal from just a small set of random samples, whenever the signal is sparse in certain transform domain. Therefore, exploring the possibilities of using different transform basis is an important task, needed to extend the field of compressive sensing applications. In this paper, a compressive sensing approach based on the Hermite transform is proposed. The Hermite transform by itself provides compressed signal representation based on a smaller number of Hermite coefficients compared to the signal length. Here, it is shown that, for a wide class of signals characterized by sparsity in the Hermite domain, accurate signal reconstruction can be achieved even if incomplete set of measurements is used. Advantages of the proposed method are demonstrated on numerical examples. The presented concept is generalized for the short-time Hermite transform and combined transform.

1. Introduction

The Hermite polynomials and Hermite functions have attracted the attention of researchers in various fields of engineering and signal processing [1–7], such as in quantum mechanics (harmonic oscillators), ultra-high band telecommunication channels, and ECG data compression using Hermite functions representation of the QRS complexes. A set of Hermite functions forming an orthonormal basis is suitable for approximation, classification, and data compression tasks [3]. Since the Hermite functions are eigenfunctions of the Fourier transform, time and frequency spectra are simultaneously approximated. Here, we are especially interested in a class of signals that are sparse in the Hermite transform domain. Note that, generally, such signals are not sparse in the Fourier transform domain. In the light of compressive sensing (CS) theory [8–12], we propose the method for efficient reconstruction of signals from its incomplete set of samples using the Hermite transform. The number of Hermite functions used in this approach is much smaller compared to the original length of the signal. The proposed approach is useful in the applications where the significant information is missing and the total signal reconstruction is required. Note that the large amount of missing signal samples may occur as a consequence of the compressed sampling strategy but also as a consequence of discarding damaged signal parts [13–15]. The theory is illustrated through examples showing that the Hermite transform based CS for certain types of signals can outperform the Fourier transform related reconstructions. Furthermore, in analogy with the time-frequency analysis based on the Fourier transform, the short-time Hermite transform is defined as a linear representation that reveals the local behavior of windowed signal parts. If the signal components are of short duration, then the short-time Hermite transform is more suitable for compressive sensing than the standard Hermite transform. Finally, the possibility of combining different transforms depending on the signal characteristics is explored.

The paper is organized as follows. The theory behind the Hermite transform and the fast method for Hermite coefficients calculation is given in Section 2. The formulation of compressive sensing approach in the Hermite transform domain is given in Section 3. The possibilities to exploit other sparsity domains based on the short-time Hermite transform and combined transform are presented in Section 4. The experimental evaluation of the proposed approach is given in Section 4, while the concluding remarks are given in Section 5.

2. Hermite Transform

The Hermite functions provide good localization and the compact support in both time and frequency domain [4]. The th order Hermite function is defined as follows:

The Hermite functions provide an orthonormal basis set for an optimal representation of different signals using the fewest number of basis functions. Signal expansion into Hermite functions, known as the Hermite transform, has been used for both 1D and 2D signals in various applications. The Hermite expansion for a signal can be defined as follows:

where are the Hermite functions and is the number of functions used for the approximation. The number of Hermite functions could be usually much smaller than the number of signal samples (). The Hermite coefficients can be defined by using the Hermite polynomials as follows:

where

represents the Hermite polynomial. An efficient procedure for calculation of Hermite coefficients can be done by applying the Gauss-Hermite quadrature [5]:

where are zeros of Hermite polynomials. By using the Hermite functions instead of polynomials, a simplified expression is obtained:

The constants are calculated as follows:

3. Compressive Sensing Formulation in the Hermite Transform Domain

Generally, the compressive sensing scenarios are focused on the new sampling strategy, which results in a large number of randomly missing samples compared to the standard sampling methods [8]. Hence, based on a small set of acquired measurements, the entire signal needs to be reconstructed. The missing samples in compressive sensing generally cannot be recovered using standard interpolation methods due to the complexity of nonstationary signals in real applications. Namely, the interpolation methods such as polynomial fit, cubic spline interpolation, or similar usually assume certain model function, which is mostly inappropriate for time domain signal modeling. Therefore, the compressive sensing reconstruction is formulated in the literature as an optimization problem (rather than interpolation) which reconstructs the signal by finding the sparsest transform domain solution corresponding to the available small set of samples.

In the CS context, we are dealing with a small set of randomly chosen samples of . Let us assume that we have only out of available samples ( is the total number of signal samples and ). The vector of available measurements is denoted as . Now we may write

where is original full signal (of length ) written in the vector form, while () is the measurement matrix. The original signal can be expressed using the Hermite transform as follows:

where is the vector of Hermite transform coefficients, while is the inverse Hermite transform matrix of size ():

The direct Hermite transform matrix is given by . In the extended form, (9) can be written as follows:

The Hermite basis functions are calculated using the fast recursive realization defined as follows:

According to (8) and (9), we have

The reconstructed signal can be obtained as a solution of linear equations with unknowns. The system is undetermined and can have infinitely many solutions. Now we assume that the signal is sparse in the Hermite transform domain. It means that the observed signal can be efficiently represented by a very small number of Hermite expansion coefficients, such that . Therefore, the optimization based mathematical algorithms should be employed to search for the sparsest solution. A near optimal solution is achieved by using the norm based minimization as follows:

where is the Hermite transform vector of signal .

In order to solve the previous minimization problem, first we need to calculate the initial Hermite transform using the available set of samples with the time support . Therefore, we observe the signal in the following form:

The initial vector of Hermite transform coefficients can be then calculated using Hermite functions () as follows:

where contains only the columns of that correspond to instants . Alternatively, we can write

In order to determine the signal support in the Hermite transform domain, the initial vector of the Hermite transform coefficients is compared by the threshold :

The exact values of Hermite coefficients at positions selected in vector are obtained as a solution of the CS problem:

The CS matrix is obtained from the inverse Hermite transform matrix , using columns that correspond to frequencies and rows corresponding to measurements with support . The system is solved in the least square sense as follows:

where () denotes the conjugate transpose operation.

Analysis of Components Reconstruction. Let us consider the isometry property of Hermite transform matrix :

The previous equation can be written as follows:

Now, observe the first term on the right side given by , particularly the sum of squared values of different Hermite functions. Figure 1 illustrates the sum of squared values of different Hermite functions calculated in the zeros of Hermite polynomials for = (128, 100, 70 and 50). Unlike in the case of the Fourier transform basis, where the sum of absolute values of complex exponential basis would be constant for any frequency , from Figure 1 we can note an approximate low-pass characteristic of the curves, meaning that the lower-order coefficients are favored compared to the higher order coefficients. Consequently, when applying the threshold for components detection, it would be easier to detect a set of low-order coefficients than the high-order ones.

Figure 1: Sum of squared values of Hermite functions for different .

4. Introducing the Short-Time Hermite Transform and Short-Time Combined Transform

4.1. Short-Time Hermite Transform

Let us assume that in (10) and define an Hermite matrix:

The short-time Hermite transform (STHT) can be defined as a composition of Hermite transform matrices whose size is defined by the window width. Without loss of generality, we may assume the nonoverlapping unit rectangular windows (with the width of samples). The total transform matrix for the STHT can be created as follows:

where is identity matrix of size (), is zero matrix, and is the total length of the signal, while denotes the Kronecker product. The STHT can be defined using transform matrix as follows:

where is a time domain signal vector of size and is a Hermite transform vector of size , while STHT is a column vector containing all STHT vectors , ( is the number of windows). Furthermore we may write

where is matrix of size that maps the global transform domain information in into local information in . In the case of windows with variable length over time (time-varying windows), the smaller matrices within will be of different sizes. Particularly, for a set of rectangular windows of sizes instead of we will have .

In the context of compressive sensing, the STHT allows us to define two types of CS problem. A common CS problem is defined by assuming that the missing samples appear in the time domain, while the sparsity is exploited in the transform domain. However, the missing samples may also appear in the STHT domain, after applying certain filter forms such as the -estimate filtering in the presence of strong noisy pulses. The two CS optimization problems are discussed in the sequel.

(1) Assume that CS procedure is done in the time domain, such that the measurements are taken from each windowed signal part: . Based on the definitions given in Section 3, the minimization problem is given by

for , where is the total number of windows, while is the inverse Hermite transform of size , where only the rows corresponding to the set of positions are kept.

(2) Let us observe the missing values in the STHT domain, where the measurements vector is denoted by . The CS problem formulation is defined using the linear relationship (26). The CS matrix, denoted by , is formed by omitting the rows from corresponding to the positions of missing values in . Then, a simple linear relationship can be established between the reduced observations in and the sparse Hermite transform vector C (corresponding to the entire signal). Hence, the CS problem is given in the following form:

It means that the missing samples in the can be reconstructed such as to provide the best concentrated . The above CS optimization problem can be solved using various existing CS reconstruction algorithms.

4.2. Short-Time Combined Transform

On the basis of the STHT, we can also define the short-time combined transform as follows:

where the transform matrix is made as a combination of Hermite transform matrices and other transforms:

where can be either the Hermite transform or some other suitable transform matrix of the same size. For instance, we may observe one example of the combination of Hermite transform and Fourier transform as a special case:

The combined matrix is thus composed of Hermite transform matrix of size , two sequential Fourier transform matrices of size , and again one Hermite transform matrix of the same size. In the case of windows with variable length over time, the short-time combined transform can be written as:

where represents the vector containing first samples in and vector contains the next samples from , while contains the last samples from such that .

We can again observe the set of measurements corresponding to different signal parts: , where . Therefore, the CS minimization problem is given by

for , where in superscript denotes that we use only the rows corresponding to available samples defined by the set . The inverse transform denoted by can be changed in each window depending on the signal characteristics.

5. Experimental Evaluation

Example 1. In order to illustrate efficiency of the proposed method, let us observe a signal that is sparse in the Hermite transform domain. The time domain signal is shown in Figure 2(a). The Hermite transform of the observed signal consists of ten components with unit amplitudes, as shown in Figure 2(b). The available number of samples is 60% of the total signal length, while the maximal order of Hermite functions used in the expansion is limited to .
After calculating the initial Hermite transform using available signal samples, the threshold is applied to select signal support in the Hermite domain. The threshold is set empirically to , , after performing a large number of experiments. Moreover, the results are not very sensitive on threshold settings, and even lower thresholds (e.g., , ) can be used, because all additional (false) components selected by such a threshold would obtain their true zero values afterwards. The suitable threshold can be even theoretically derived (under certain assumptions). Namely, the effects caused by missing samples in time domain can be modeled by the noise in the Hermite transform domain. The noise can be described using two random variables: the first one corresponds to the noise appearing at nonsignal components (noise alone), while the second one corresponds to the noise appearing at the signal components positions. The two random variables can be described by the corresponding mean values and variances that could be further employed to define the threshold level. Namely, with certain probability, we may define the threshold value that is just above the level of noise components. Consequently, such a threshold would select only the signal components determining the signal support used in the reconstruction procedure. Since the threshold derivation would require significant space and analyses, it could be a topic of some further work.
The selected components in the Hermite transform domain are shown in Figure 3. Using the detected positions of components, the exact values of Hermite coefficients are calculated according to (20) (Figure 3). The reconstructed signal and the absolute error between original and reconstructed signal are shown in Figures 4(a) and 4(b), respectively. For the comparison, the reconstruction based on the Fourier transform is considered. However, it is shown that the Fourier transform based reconstruction cannot be used in this case, because it produces a significant error (Figure 4(c), 10 largest Fourier transform components are used for reconstruction). The reason is in the fact that the signals that are sparse in Hermite transform domain do not exhibit sparsity in the Fourier transform domain (or some other transform domains). This means that the exact results, in this case, can be obtained only using the Hermite transform providing the sparse signal representation. Similar conclusion holds for other transform domains where the signal is not sparse (Fractional Fourier, Wavelet, etc.).

(a)

(b)

Figure 2: The original signal: (a) time domain representation, (b) Hermite transform of signal calculated using Hermite functions.

Figure 3: Hermite coefficients of signal with missing samples -○, and the reconstructed Hermite coefficients -×.

(a)

(b)

(c)

Figure 4: The reconstructed signal: (a) time domain representation, (b) the absolute error between original and reconstructed signal, (c) for comparison: reconstructed signal based on the Fourier transform reconstruction approach (-°) and the corresponding error (solid line).

Example 2. In this example we have observed different sparse signals (having 10 components) in the Hermite domain through the 1000 repetitions of the previously described procedure. The aim is to test how the accuracy of the proposed method changes for different number of missing samples. Namely, the number of missing samples was increased in steps of 10% starting from zero up to 80% of missing samples. For each instance on the -axes, we performed 1000 iterations with different 10-component signals, and the mean square errors between original and reconstructed signals are calculated and shown in Figure 5(a). We can observe that, up to 70% of missing samples, the MSE is still low, but it starts to increase significantly afterwards (for values higher than 70%). This means that the reliable reconstruction cannot be guaranteed for such a large number of missing samples (such as 80% of missing samples).
Additionally, the proposed approach is tested in the presence of external additive Gaussian noise. Namely, the MSE is calculated for different values of SNR (for each SNR, we assume 1000 realizations of random noise). The results are presented in Figure 5(b), showing that the proposed approach can be efficiently used for SNR > 8 dB. For SNR < 8 dB the MSE starts to increase rapidly.

(a)

(b)

Figure 5: (a) MSE calculated for 1000 realizations and for different number of missing samples, (b) MSE calculated for different values of SNR in 1000 realizations of external noise.

Example 3. This example illustrates the concept of combined Hermite-Fourier transform and CS reconstruction. Observe the signal with the total length of samples composed of two parts: first part exhibits sparsity in the Hermite transform domain with components (HT components) at positions 7 and 12 with amplitudes 0.8 and 0.65; second part exhibits sparsity in the Fourier transform domain with components (FT components) at positions −3 and 4 with amplitudes 1 and 0.8, respectively. Now consider the compressive sensing scenario with only 50% of samples available (original full length signal and available measurements are shown in Figure 6). The original short-time combined transform obtained using windows width equal to 16 samples is shown in Figure 7(a), where the white fields represent 0 value. Thus, the signal represented by a full set of samples is sparse meaning that it is fully concentrated on a few nonzero components. In the case when we deal with 50% of missing samples, the corresponding short-time combined transform is shown in Figure 7(b). Note that in this case the sparsity is disturbed due to the missing samples. Instead of zero values at nonsignal position, the noise appears as a consequence of missing samples.

Figure 6: Missing samples and available measurements in time domain.

(a)

(b)

Figure 7: (a) Original short-time combined transform of full length signal (FT: Fourier transform, HT: Hermite transform), (b) short-time combined transform of available measurements.

The threshold based components selection and signal reconstruction (in analogy to (18) and (20)) is applied to both parts of the signal, FT components, and HT components (Figures 8(a) and 8(b), resp.). Based on the identified positions of components, the exact signal reconstruction is done using the procedure presented in Section 3. The reconstructed signal components are shown in Figure 9.

(a)

(b)

Figure 8: Selecting the components of interest by applying threshold to the following: (a) FT components and (b) HT components.

Figure 9: Reconstructed signal components.

6. Conclusion

The possibility of using the Hermite transform in compressive sensing applications was explored in this work. The compressive sensing setup and signal reconstruction approach in the Hermite transform domain were defined. A simple and fast procedure for the total reconstruction of signals using selected Hermite transform coefficients provides the results very close to the original signal even when we deal with significant missing information. It is important to emphasize that the proposed concept in the Hermite transform domain can be also combined with other known compressive sensing solvers. The entire concept is generalized and extended by defining the short-time Hermite transform as well as the short-time combined transform. These two transforms open more possibilities to apply the compressive sensing approach in different scenarios.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the Montenegrin Ministry of Science, Project Grant “New ICT Compressive Sensing Based Trends Applied to: Multimedia, Biomedicine and Communications” (ACRONYM: CS-ICT).

References

A. Sandryhaila, S. Saba, M. Püschel, and J. Kovacevic, “Efficient compression of QRS complexes using hermite expansion,” IEEE Transactions on Signal Processing, vol. 60, no. 2, pp. 947–955, 2012.

M. Lagerholm, G. Peterson, G. Braccini, L. Edenbrandt, and L. L. Sörnmo, “Clustering ECG complexes using hermite functions and self-organizing maps,” IEEE Transactions on Biomedical Engineering, vol. 47, no. 7, pp. 838–848, 2000.

B. Beliczynski, “Approximation of functions by multivariable hermite basis: a hybrid method,” in Adaptive and Natural Computing Algorithms, vol. 6593 of Lecture Notes in Computer Science, pp. 130–139, Springer, 2011.

Y. K. Alp and O. Arıkan, “Time-frequency analysis of signals using support adaptive Hermite-Gaussian expansions,” Digital Signal Processing, vol. 22, no. 6, pp. 1010–1023, 2012.

A. Krylov and D. Kortchagine, “Fast Hermite projection method,” in Proceedings of the International Conference on Image Analysis and Recognition (ICIAR '06), pp. 329–338, Póvoa de Varzim, Portugal, September 2006.

D. Kortchagine and A. Krylov, “Projection filtering in image processing,” in Proceedings of the International Conference on Computer Graphics and Applications, pp. 42–45, Moscow, Russia, August-September 2000.

I. Orovi, S. Stankovi, T. Chau, C. M. Steele, and E. Sejdić, “Time-frequency analysis and hermite projection method applied to swallowing accelerometry signals,” EURASIP Journal on Advances in Signal Processing, vol. 2010, Article ID 323125, 2010.

R. G. Baraniuk, “Compressive sensing,” IEEE Signal Processing Magazine, vol. 24, no. 4, pp. 118–124, 2007.

D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006.

E. J. Candes, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489–509, 2006.

D. Vukobratovic and A. Pizurica, “Compressed sensing using sparse adaptive measurements,” in Proceedings of the Symposium on Information Theory in the Benelux (SITB '14), Eindhoven, The Netherlands, May 2014.

D. Angelosante, G. B. Giannakis, and E. Grossi, “Compressed sensing of time-varying signals,” in Proceedings of the 16th International Conference on Digital Signal Processing (DSP '09), pp. 1–8, Santorini-Hellas, Greece, July 2009.

L. Stankovic, S. Stankovic, and M. Amin, “Missing samples analysis in signals for applications to L-estimation and compressive sensing,” Signal Processing, vol. 94, no. 1, pp. 401–408, 2014.

S. Stanković, I. Orović, and L. Stanković, “An automated signal reconstruction method based on analysis of compressive sensed signals in noisy environment,” Signal Processing, vol. 104, pp. 43–50, 2014.

S. Stanković, L. Stanković, and I. Orović, “Relationship between the robust statistics theory and sparse compressive sensed signals reconstruction,” IET Signal Processing, vol. 8, no. 3, pp. 223–229, 2014.