Abstract

Bioimpedance spectroscopy (BIS) has become an important clinical indicator for monitoring the pathological status of biological tissues, and multifrequency simultaneous measurement of BIS may provide more accurate diagnostic information compared with the traditional frequency-sweep measurement technology. This paper proposes a BIS multifrequency simultaneous measurement method based on multifrequency mixed (MFM) signal excitation and a Nuttall windowed interpolation FFT algorithm. Firstly, the excitation source adopts the nine-frequency MFM signal , which has excellent spectral characteristic and is very suitable for BIS measurement. On this basis, a Nuttall window is adopted to truncate sample data, and an interpolation FFT algorithm based on Nuttall window is built to perform spectral analysis, in which the parameter correction formula is provided based on polynomial approximation. A BIS measurement simulation experiment is performed on an RC three-element equivalent circuit, and results on the 9 primary harmonic frequencies ranging from 3.9 kHz to 1 MHz show a high accuracy with the impedance amplitude relative error , and the phase absolute error . This paper validates the feasibility of BIS multifrequency simultaneous measurement method and establishes an algorithm foundation for the development of practical broadband BIS measurement system.

1. Introduction

Bioimpedance spectroscopy (BIS), which performs measurement of complex electrical impedance of biomaterials over a certain frequency range, has been extensively used in many electrochemical and biomedical applications since it has been shown that BIS can provide physiological and pathological information of biological tissues [1]. In recent years, BIS technology has been increasingly and widely studied and adopted in clinical diagnoses on tissue ischemia [2], mammary cancer [3], and lung cancer [4].

Traditionally, the most commonly used BIS measurement technique is the traditional frequency-sweep (FS) approach which performs measurement of impedance at a single frequency point in a time and sweeps the specific frequencies in the range of interest [5, 6]. To date, the FS method remains popular for its simplicity and ability to describe the stationary properties of linear time-invariant bioimpedance [7]. The main advantage of the FS approach lies in the fact that it can usually obtain high signal-to-noise ratio (SNR) at the cost of the long measuring time, which usually lasts one second to tens of seconds for a complete sweep [8, 9]. However, the physiological status of living body, such as the cardiovascular system, is time-varying and dynamic because of its blood flow, heartbeat, and other factors. Therefore, in the situation of time-varying bioimpedance measurement, the main drawback of the FS method is its sweeping cycle, which can be much longer than the bioimpedance time variations [8]. So the BIS data based on the FS measurement technology cannot accurately represent the instantaneous impedance spectra and may lost important diagnostic information [10, 11].

In recent years, the multifrequency simultaneous (MFS) measurement technique, which applies a broadband excitation and gains its resulting frequency response by means of spectral analysis using the fast Fourier transform (FFT) algorithm, has been becoming popular for time-varying bioimpedance [8, 1215]. Thanks to multiple frequencies being excited at the same time, the time to acquire a complete BIS is drastically reduced to tens of milliseconds [4, 16]. Compared to the FS approach, this MFS approach can obtain the impedance frequency response at different frequencies simultaneously and is highly desirable for time-varying biological systems, such as the respiratory or cardiovascular systems [17, 18].

Proper spread-spectrum excitation is crucial for MFS measurement of BIS. Sanchez et al. [15] had a review on several broadband excitations such as maximum length binary sequences (MLBS), chirp. As a binary signal, the MLBS has the advantage of implementation simplicity and has a dense (continuous) amplitude spectrum with plateau envelope where most of the energy is concentrated. But in most BIS applications, only a discrete set of frequencies is required to fit the impedance data to a model like Cole equation [19], and BIS devices usually need several to tens of isolated frequencies (usually with large frequency intervals) over a wide range of frequency (often from 5 kHz to 1 MHz) [20]. Consequently, broadband excitation with sparse spectrum where energy is distributed on finite frequency points with large frequency spacing is ideal for BIS measurements.

The authors previously proposed a multifrequency mixed (MFM) signal whose majority (65% or more) of energy is concentrated on several expected 2nth primary harmonics [21], which provides a proper broadband excitation for MFS measurement BIS. With the broadband excitation, the complex impedance spectrum then can be just calculated as the ratio of the complex Fourier coefficients of the response voltage to the complex Fourier coefficients of the excitation current signal after FFT operations [15, 22].

However, FFT operations usually bring spectral leakage since rigorous integer-period sampling is usually hard to achieve. The digital Fourier transform (DFT) theory always assumes that the input sequence is periodic (assumed period extension), so if the input sequence finishes on an integer number of periods everything is fine. Otherwise, discontinuity occurs in the assumed extended sequences and brings spectral leakage and picket fence effect in FFT operation [23, 24]. Unfortunately, integer-period sampling is usually hard to achieve, since sampling device is hard to begin and ends exactly at the signal’s head and tail, respectively. Moreover, as a rectangular wave, the MFM signal contains a great number of undesired harmonic components according to the Fourier series theory [25]. These undesired frequency components, especially those whose frequencies are higher than half of the sampling rate, will lead to spectrum aliasing when performing FFT operations and ultimately degrade the measurement precision [26].

Fortunately, the undesired effects of the spectral leakage can be minimized by weighing the samples by a suitable time window, and the picket fence effect can be reduced by adopting interpolation algorithms [27]. It is proved that windows with big side-lobe attenuation and high side-lobe roll-off rate can sufficiently reduce the spectral leakage [28].

This paper synthesizes a nine-frequency MFM signal according to the previously proposed method [21] as the broadband excitation, and realizes multifrequency simultaneous (MFS) measurement of BIS based on a Nuttall windowed interpolation FFT Algorithm. The Nuttall window [29], which has excellent characteristics on side-lobe and main-lobe, is adopted to suppress spectral leakage before FFT, and a Nuttall windowed double-spectral-line interpolation FFT algorithm is developed to reduce the error of picket fence effects. A BIS measurement experiment on an RC three-element equivalent circuit is simulated to evaluate the performance of the proposed algorithm.

2. Synthesis of the MFM Signal

According to the MFM signal synthesis principle introduced by the literature [21], let ; then the nine-frequency MFM signal is synthesized and the time domain waveform of within one period is shown in Figure 1. The harmonic amplitude spectra and power percentage spectra of are shown in Figures 2(a) and 2(b), respectively, where the nine expected 2nth harmonics (the red solid points in Figure 1), namely, the 1st, 2nd, 4th, 8th, 16th, 32nd, 64th, 128th, and 256th harmonics (we called them primary harmonics thereafter) are obviously prominent. The amplitude spectra , power spectra , and initial phase of are shown in Table 1, respectively.

As shown in Figure 2 and Table 1, the synthesized MFM signal has nine 2nth primary harmonic components which occupy up to 65.52% of the total power. The period of one code element width of is 0.5 μs, and its total 512 code elements last a fundamental period μs, which means that the fundamental frequency of is . The nine primary harmonic frequencies of are also listed in Table 1, which expand from 3.90625 kHz to 1 MHz, and cover the main frequency range in most BIS measurements. As a rectangular wave, the MFM signal has a crest factor (CF) of 1.235 ( , lower than the CF of the sinusoid 1.414) according to the computation method mentioned in [30], which is desirable for BIS of living biological systems because a low crest factor may limit the applied current peaks and set biological systems within the linear regime [12, 31]. Furthermore, as a binary-valued function, the MFM signal can be generated by FPGA easily, and the bandwidth of its harmonics can be variable by changing the FPGA’s operating clock, and the number of the expected primary harmonics is also adjustable. Based on the features mentioned above, the newly designed wideband excitation source establishes a good foundation for fast measurement of BIS.

3. Nuttall Windowed Interpolation FFT Algorithm

3.1. Nuttall Window

The Nuttall window [29] is a combination of cosine windows and can be expressed in the time domain as follows: where is the number of terms of the cosine function; denotes the length of the window and ; the coefficient should meet the constraint conditions: , .

The coefficient groups of typical Nuttall windows shown in Tables 2 and 3 list the side-lobe characteristics of typical Nuttall windows and other cosine-combined windows [29].

The side-lobe attenuation level and asymptotic decay rate directly affect the FFT-based spectral analysis results [28]. Among the windows shown in Table 2, the 4-term 3-order Nuttall window has the fastest side-lobe asymptotic decay rate 30 dB/oct and an acceptable side-lobe attenuation level −82.6 dB, which is adopted to truncate sampled data and suppress spectrum leakage in FFT.

3.2. Nuttall Windowed Interpolation FFT Algorithm

If denotes an analog signal with multiple harmonic components, its discrete-time form can be expressed as where denotes the number of harmonic components and is an integer representing the harmonic order. When , , , are the frequency, amplitude, and initial phase of the fundamental harmonic, respectively. When , , , denote the frequency, amplitude, and initial phase of the th harmonic respectively. is the digital sampling rate.

If is truncated by Nuttall window , then the windowed signal after discrete Fourier transform (DFT) can be expressed in frequency domain as where denotes the continuous spectrum of , the frequency resolution, and . is the truncated data length. is the peak spectral of the th harmonic and .

When signal is sampled asynchronously, its peak frequency of the th harmonic can hardly be exactly on the discrete spectral line. As shown in Figure 3, suppose that the largest and the second largest spectral amplitude lines near the peak point are and , respectively; they have the relationship And the amplitudes of these two spectral lines, and , can be expressed as

Set And according to (4), set Then is in the range [−0.5, 0.5] and can be expressed as

Let be the inverse function of (8); namely, ; then the frequency can be corrected by : And set the amplitude as the weighted average of the amplitudes of the and spectral lines:

According to (3) and (7), the initial phase of the th harmonic can be corrected as follows:

According to (1), the spectrum amplitude function of Nuttall window is where denotes the spectrum function of rectangular window and .

Let ; then (12) can be expressed as

Because , , and is generally very large, (13) can be approximately expressed as

Substituting (14) into (8), we can obtain the approximation of

Similarly, substituting (14) into (10), we can obtain the approximation of (10): where denotes a function of and can be expressed as follows:

From (11) and (14), the phase can be corrected as follows:

As shown in (9), (16), and (18), , , are all related to , but has a complicated relationship with in (8) and with in (17). However, the value of and could hardly be figured out by direct analytical method.

This paper proposes a simple and efficient method based on polynomial approximation to estimate the values of and . Let be a series of values in the range ; then the corresponding values of and can be computed according to (14), (8), and (17). Afterwards, based on the known mapping relationship between and , the inverse function of , namely, , can be approximated by a quintic polynomial of :

Also, based on the known mapping relationship between and , can be approximated by a quintic polynomial of :

For every harmonic, we always can find the largest and the second largest spectral lines and and their amplitudes and . Then must have a concrete value according to (6), and the values of and consequent can be calculated by (19) and (20), respectively. With the determined and , the frequency, amplitude, and initial phase of the harmonic, , , , are corrected according to (9), (16), and (18).

As to the MFM signal , which contains multiple harmonic components, set , 2, 4, 8, 16, 32, 64, 128, and 256, and the corresponding primary harmonic parameters , , , can be computed, respectively.

4. Measurement Experiment

On the basis of Nuttall windowed interpolation FFT algorithm mentioned above, we can establish a BIS multifrequency synchronous measurement scheme as shown in Figure 4.

In Figure 4, the excitation current source is the MFM signal , flows through the load , and the voltage drop occurs. In biological tissue, the extracellular and intracellular fluids can be simply equivalent to a resistance, and the cell membrane can be equivalent to a capacitance, which supports the hypothesis that an RC three-element model can represent equivalently the impedance of a biological tissue [1]. This paper adopts an RC three-element equivalent circuit model (dashed box in Figure 4) as the load . In this experiment, we use the typical parameters , , and , since the human body impedance range is 300~1500 Ω, the capacitance value range is 2 pF~22 nF [32].

In order to prove the correctness of the scheme proposed above, we design a BIS measurement simulation experiment. When the RC three-element equivalent circuit is driven by the periodical current source signal , the response voltage signal at any time can be simulated according to the method mentioned in the literature [21]. The simultaneous samplings on and are simulated and the corresponding discrete sequence and are obtained, respectively. Figure 5(a) shows the simulated in a half-period, and Figure 5(b) gives an example of the simultaneous sampling of the excitation and the response .

Like most broadband BIS measurements using periodic excitations, the complex impedance spectrum of the tested RC model can be just calculated as the division of the voltage and current complex Fourier coefficients obtained from their respective FFT operations [15, 22]. and are firstly truncated using the Nuttall window; then the Nuttall windowed interpolation FFT algorithm is performed on the two sequences, respectively. Since rigorous integer-period sampling is hard to achieve in practice, ten or more sampling cycles are usually needed for the windowed FFT algorithm. In this paper, we complete the sampling in about ten cycles, which lasts about 2.56 milliseconds (the period of ). The complex Fourier coefficients of and , namely, and   ( ) are obtained after FFT. Finally, the tested impedance amplitudes and phases are calculated according to the following formula:

The BIS measurement results at the nine primary harmonics are shown in Table 4, in which the impedance amplitude relative error , and the phase absolute error .

5. Conclusion

Multifrequency simultaneous (MFS) measurement of BIS can greatly reduce measurement time and grasp the transient physiological status of living body compared with the traditional frequency-sweep measurement technology. This paper proposes a BIS multifrequency simultaneous measurement approach based on an MFM excitation and a Nuttall Windowed interpolation FFT algorithm. On the basis of the nine-frequency MFM signal excitation, a Nuttall window is adopted to truncate the sample data, and an interpolation FFT algorithm based on Nuttall window is developed to perform spectral analysis, in which the parameter correction formula is provided based on polynomial approximation. A BIS measurement simulation experiment is performed on a RC three-element equivalent circuit, and results show a high accuracy with the impedance amplitude relative error , and the phase absolute error . This paper establishes an algorithm foundation for development of practical broadband BIS measurement systems.

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

This study has been partly supported by two grants from the National Natural Science Foundation of China (30900317, 61273271), a grant from the Scientific Research Plan of Education Bureau of Shaanxi Province, China (12JK0527), and a grant from China Postdoctoral Science Foundation (20110491674).