#### Abstract

Sub- and ultraharmonics generation by ultrasound contrast agents makes possible sub- and ultraharmonics imaging to enhance the contrast of ultrasound images and overcome the limitations of harmonic imaging. In order to separate different frequency components of ultrasound contrast agents signals, nonlinear models like single-input single-output (SISO) Volterra model are used. One important limitation of this model is its incapacity to model sub- and ultraharmonic components. Many attempts are made to model sub- and ultraharmonics using Volterra model. It led to the design of mutiple-input singe-output (MISO) Volterra model instead of SISO Volterra model. The key idea of MISO modeling was to decompose the input signal of the nonlinear system into periodic subsignals at the subharmonic frequency. In this paper, sub- and ultraharmonics modeling with MISO Volterra model is presented in a general framework that details and explains the required conditions to optimally model sub- and ultraharmonics. A new decomposition of the input signal in periodic orthogonal basis functions is presented. Results of application of different MISO Volterra methods to model simulated ultrasound contrast agents signals show its efficiency in sub- and ultraharmonics imaging.

#### 1. Introduction

Medical diagnostic using ultrasound imaging was greatly improved with the introduction of ultrasound contrast agents. In ultrasound imaging, contrast agents are microbubbles [1]. Historically, the important difference between the acoustic impedance of the tissue and the gas encapsulated within the microbubbles was the first step to improve the contrast of echographic images. However, the contrast was still improved by taking into account the nonlinear behavior of microbubbles. In fact, when microbubbles were insonified by a sinusoidal excitation, they respond by generating harmonic components [2]. For example, second harmonic imaging (SHI) [3] consists in transmitting a signal at frequency and receiving echoes at twice the transmitted frequency . However, harmonic generation during the propagation of ultrasound in the nonperfused tissue limits the contrast [4].

Many years ago, experimental studies have shown the existence of subharmonics at [5] and ultraharmonics at [6] in the microbubble response under specific conditions of frequency and pressure. The absence of these components in the backscattered signal by the tissue has enabled the introduction of sub- and ultraharmonics as an alternative of the harmonic imaging in order to enhance the contrast. Sub- and ultraharmonic imaging consists of transmitting a signal of frequency and extracting components at .

Many models have been developed to understand the dynamics of the microbubble [2]. Microbubble oscillation can be accurately described using models such as Rayleigh-Plesset modified equation [7–9]. However, to enable optimal separation of harmonic components, other nonlinear models like single-input single-output (SISO) Volterra model have been preferred [10]. A well known limitation of SISO Volterra model is its capacity to model exclusively harmonic components sub- and ultraharmonics are not modeled [11].

To overcome this difficulty, Boaghe and Billings [12] have proposed a multiple-input single-output (MISO) Volterra-based method. Input signals are specified by having sub-harmonic component at frequency . This approach has been applied in ultrasound medical imaging [13].

However, neither Boaghe and Billings [12] nor Samakee and Phukpattaranont [13] have clearly justified the required conditions to design a MISO Volterra decomposition able to model sub- and ultraharmonics.

To answer this untreated point, we propose a more general framework which firstly gives a clear justification regarding the choice of the model and secondly can offer interesting alternatives.

This paper is organized as follows: after recalling Volterra model and presenting the general framework of MISO Volterra methods, simulations of contrast ultrasound medical imaging followed by results are presented. Finally, a discussion completed by a conclusion closes the paper.

#### 2. SISO Volterra Model

Volterra series were introduced like Taylor series with memory [10]. Let and be, respectively, the input and the output signals in the discrete time domain of the nonlinear system (see Figure 1). The output of Volterra model of order and memory is given in [14]. Note that, in our study focused on ultrasound imaging, a third-order Volterra model is sufficient for the available transducers bandwidths. The output of SISO Volterra model of order and memory is given by where is the kernel of order of the filter, with .

Equation (1) could be rewritten as follows where the output signal is: where is the length of the signal , and the vector of kernels is where the input matrix is with vector with .

The vector of kernels is calculated to minimize the mean square error (MSE) between and according to the equation where is the symbol of the mathematical expectation.

Vector is calculated using the least squares method if is invertible. Otherwise, regularization techniques can be used.

Nevertheless, as reported in [12], it is not possible to model sub- and ultraharmonics with SISO Volterra model under this formulation. This is due to the fact mentioned in [12] that SISO Volterra model can only model frequencies at integer multiples of the input frequency.

To overcome this limitation, Boaghe and Billings [12] proposed a MISO Volterra-based solution and not any more a SISO Volterra. This point is discussed in Section 3.

#### 3. General Framework of MISO Volterra Model

According to Boaghe and Billings’ claims [12], it is possible to model sub- and ultraharmonic components of the signal if the excitation signal to Volterra model has the sub-harmonic component at . The solution proposed by Boaghe and Billings [12] to show up the sub-harmonic component at frequency is to decompose the input signal into multiple input signals , each signal having frequency components at and . From our point of view, Boaghe and Billings’ approach [12] claimed two conditions that are intrinsically coupled by the choice of the decomposition method as follows:(i)the input signal to Volterra model has sub-harmonic frequency component at ;(ii)Volterra system is a MISO system described by

The block diagram of MISO Volterra model is presented in Figure 2.

A third condition that is not really explained in [12], however, it is a crucial condition to carry out this modeling procedure. It is the orthogonality condition between each multiple input of MISO Volterra model. Taking into account this third condition makes it possible to generalize Boaghe and Billings’ approach presented in [12] as follows: where are coefficients to be adjust and is the periodic orthogonal basis functions having a spectral component at . Different bases could be proposed. In this study, two bases are presented as follows.(1)In [12] a first periodic basis of orthogonal functions is proposed as follows: where is the sampling period, represents the convolution product, and is the rectangular function equal to when and equal to zero otherwise. Note that this approach is MISO1. (2)In the present work, a new periodic basis of orthogonal functions is presented as follows: where is the Hilbert transform of and . Note that this second is MISO2 approach.

For our application in contrast medical imaging, the sub-harmonic frequency is [5–7], so .

As an illustration, when and , the decomposition is written: (1)for the first basis, as follows: with , and (2)and for the second basis, as follows: with , and where is the Dirac function. Finally, and can be simply rewritten as follows:

The two signals and for the two previous bases are represented in Figure 3.

**(a)**

**(b)**

It is obvious to show that for the two bases, the signals and are orthogonal because (From a statistical point of view, the two signals and are orthogonal if and only if . If and are stationary and ergodic, then .). The algebraic area of the signal , shown in Figure 3, is equal to zero.

Finally, if the components are orthogonal to each other, then this also means that the output of Volterra model can be decomposed as follows: where the components are also orthogonal to each other. A proof of this propriety is given in Appendix A.

The consequence of this statement is that MISO Volterra model can be considered as parallel SISO Volterra models as depicted in Figure 4.

#### 4. Simulations

To validate the different proposed bases and to quantify its performances for application in contrast ultrasound medical imaging, realistic simulations are proposed. To carry out the simulations, the free simulation program BUBBLESIM developed by Hoff [7] was used to calculate the oscillations and scattered echoes for a specified contrast agent and excitation pulse. A modified version of Rayleigh-Plesset equation was chosen. The model presented by Church [15] and then modified by Hoff [7] is based on the theoretical description of microbubbles as air-filled particles with surface layers of elastic solids. In order to simulate the mean behavior of a microbubble cloud, we hypothesized that the response of a cloud of microbubbles was times the response of a single microbubble with the mean properties.

The incident burst to the microbubble is a sinusoidal wave of frequency MHz (The resonance frequency of a microbubble of *μ*m is about MHz. Therefore, the emission frequency at MHz is nearly the double of the resonance frequency.) To ensure the presence of sub- and ultraharmonics with moderate destruction of microbubbles, Forsberg et al. have proposed in [16] a pressure range from MPa to MPa. To limit the destruction of microbubbles, we set the pressure level to the lowest value at MPa. The burst consists of 18 cycles. The sampling frequency is MHz. The parameters of the microbubble are given in the Table 1 [13].

#### 5. Results

In this research, the performances of different modeling methods are evaluated qualitatively and quantitatively.

##### 5.1. Qualitative Evaluation

To evaluate qualitatively the two MISO methods, MISO1 (with the basis proposed in [12]) and MISO2 (with the new basis proposed in the present work) with respect to SISO Volterra method, temporal representations of and , and spectral representations and of the nonlinear system backscattered by the contrast agent in nonlinear mode are presented in Figure 5.

**(a)**

**(b)**

Results presented in Figure 5 are obtained for a signal to tissue ratio and using Volterra model of order and memory .

To better distinguish the different harmonic components of ultrasound signal, six cycles of *μ*s are presented in Figure 5(a), and a bandwidth of MHz covering the 3 harmonics potentially accessible in ultrasound imaging is presented in Figure 5(b). For both types of representations, the fundamental frequency, harmonics, sub- and ultraharmonics are well apparent. In Figure 5(a) (top), only harmonic components at , , and are modeled by SISO Volterra. This result confirms that SISO Volterra system is unable to correctly model sub- and ultraharmonics at frequencies and . In Figure 5 (middle, bottom), all the spectral components are correctly modeled validating the two MISO approaches.

##### 5.2. Quantitative Evaluation

To determine accurately the performances of the two methods and to know which Volterra approach provides the best performances a quantitative study is necessary. The relative mean square error (RMSE) defined as follows: is evaluated for different noise levels at the system output. The noise level, adjusted as a function of SNR, is Gaussian and white. Ten realizations are made to evaluate the fluctuations of RMSE. RMSE for and dB is reported in Figure 6. A zoom in Figure 6(d) shows the fluctuations of the EQMR around a mean value.

**(a)**

**(b)**

**(c)**

**(d)**

The main result of these simulations shows that regardless the SNR values, MISO Volterra methods provide a much better RMSE than SISO Volterra method. In fact, a gap between SISO Volterra method and the two methods MISO1 and MISO2 going from 5 to dB can be obtained depending on the SNR conditions. These results confirm that SISO Volterra method is not suitable for sub- and ultraharmonic modeling. A zoom on Figure 6(d) emphasizes the small fluctuations of the RMSE. This result shows the robustness of the two MISO Volterra approaches towards noise.

Note that the RMSE obtained with the two MISO Volterra approaches are similar and follows the same trend. However, a small advantage in favor of MISO2 method with respect to MISO1 method for memory values smaller than 6 is noted.

Finally, the more the memory increases, the more the RMSE decreases, indicating that the different methods tend asymptotically toward the optimal solution.

#### 6. Discussions and Conclusions

In the present research, we proposed a general framework that describes harmonic, sub-, and ultraharmonics modeling using Volterra decomposition. This framework allowed us to highlight three essential criteria instead of two, to accurately model sub- and ultraharmonics:(i)as suggested in [12], the basis should be periodic of period ; (ii)as suggested in [12], Volterra system should be a MISO system; (iii)as suggested in this work, the decomposition of the input signal to Volterra model must be done with an orthogonal basis.

This general framework has also justified the different steps of the decomposition thus allowing to propose new periodic orthogonal bases more efficient. It is the same for the choice of the order of Volterra model, which was limited to three. In fact, for more or less severe constraints on the ultrasound transducers bandwidth, the order can be reduced or increased.

This more general formulation provides a methodological basis for optimal sub- and ultraharmonics contrast imaging and opens a new research axis for more efficient periodic orthogonal bases of MISO Volterra systems and also for new MISO systems based on Hammerstein models or Wiener models.

#### Appendix

#### A. Decomposition of MISO Volterra Model of 2 Inputs to 2 SISO Volterra Models

A MISO Volterra model with inputs is equivalent to SISO Volterra models if and only if the mean square error to be minimized between and is the same in both cases. We will determine the conditions that must be satisfied by the inputs and of MISO Volterra when , to have this equivalence.

Volterra kernels are calculated using the least squares method by minimizing the mean square error (MSE) between and the modeled signal :

For MISO Volterra model, the decomposition of into and such that requires that . It follows that . The error to be minimized is

For the 2 SISO Volterra models of inputs and and outputs and , respectively, the error to be minimized is

A MISO Volterra model could be seen as 2 SISO Volterra models if (A.2) and (A.3) are equal. This equality gives

One possible solution is that each term of the equation is equal to zero:

Therefore where means orthogonal. Elsewhere, and are calculated according to (1). That implies that

One possible solution is that each term of the equation is equal to zero. For the first term, we get

The last equation implies that . For the other terms in (A.7), we obtain the same conclusion . Therefore, and are orthogonal if and are also orthogonal.

Elsewhere, if and are the estimations of and , then and . This means that the orthogonality of and implies the orthogonality of each couple formed by the four signals presented in (A.5). This is true if and only if .

Therefore, a MISO Volterra model with two inputs could be treated as two SISO Volterra models if the two inputs are orthogonal. This result could be generalized for MISO Volterra model with inputs.

#### Acknowledgment

The authors would like to thank the Lebanese council of scientific research (CNRSL) for financing this work.