[en] The description of complex systems requires to select and to compact the relevant information. The wavelet theory constitutes an appropriate framework for defining adapted representation bases obtained from a controlled hierarchy of approximations. The optimization of the wavelet analysis depend mainly on the chosen analysis method and wavelet family. Here the analysis of the harmonic oscillator wave function was carried out by considering a Spline bi-orthogonal wavelet base which satisfy the symmetry requirements and can be approximated by simple analytical functions. The goal of this study was to determine a selection criterion allowing to minimize the number of elements considered for an optimal description of the analysed functions. An essential point consists in utilization of the wavelet complementarity and of the scale functions in order to reproduce the oscillating and peripheral parts of the wave functions. The wavelet base representation allows defining a sequence of approximations of the density matrix. Thus, this wavelet representation of the density matrix offers an optimal base for describing both the static nuclear configurations and their time evolution. This information compacting procedure is performed in a controlled manner and preserves the structure of the system wave functions and consequently some of its quantum properties