Aspects of Using Chirp Excitation for Estimation of Bioimpedance Spectrum

Short frequency swept signals, known as chirps, are widely used as excitation or stimulus signals in various areas of engineering as radar and sonar techniques, acoustics and ultrasonics, optical and seismological studies, but also in biomedical investigations, including bioimpedance measurement and impedance spectroscopy (Müller & Massarani, 2001; Misaridis & Jensen, 2005; Barsoukov & Macdonald, 2005; Nahvi & Hoyle, 2009).


Introduction
Short frequency swept signals, known as chirps, are widely used as excitation or stimulus signals in various areas of engineering as radar and sonar techniques, acoustics and ultrasonics, optical and seismological studies, but also in biomedical investigations, including bioimpedance measurement and impedance spectroscopy (Müller & Massarani, 2001;Misaridis & Jensen, 2005;Barsoukov & Macdonald, 2005;Nahvi & Hoyle, 2009).
Signal processing in the chirp-based applications is often combined with pulse compression via cross-correlation procedure and Fourier analysis.In this chapter, a similar approach is proposed for estimation of the frequency response (the impedance spectrum) of electrical bioimpedance.An advantage of the chirp-based method is that the characteristics of a biological object can be obtained in a wide frequency range during a very short measurement cycle, which nearly eliminates the influence of low-frequency biological processes (heart beating, breathing, pulsation of blood) to the result of measurement.
The changes of a spectrum monitored at sequent time intervals by means of the Fourier Transform (known as spectrogram) is an informative base for interpreting the processes in biological objects.Furthermore, the signal treatment by the cross-correlation yields a better noise immunity for the measurement system (Barsoukov & Macdonald, 2005), and adds some alternatives for estimation of the bio-impedance properties and behavior.
On the other hand, the bioimpedance is not a constant value, but a function, which is changing in time due to numerous biological processes in the living tissue, and as a matter of fact, is a function of frequency and time altogether as Ż=Z (jω,t).
Injection of the excitation current (stimulus) I exc (t) with known parameters into biological object courses the response voltage V z (t), analysis of which enables to estimate the impedance spectrum Ż(ω) of the object being under the investigation.
Fig. 1 illustrates the path of excitation current through the cells of a tissue, where r ext corresponds to the extracellular resistance, and the resistive components r int together with intercellular capacitances C c constitutes the intracellular impedance (Grimnes & Martinsen, 2008;Min & Parve, 2007).

Fig. 1. Formation of the electrical bioimpedance of tissue
Very often, the response signal is analyzed using the Fourier Transform F(V z (t)) to get information about the frequency-dependent state and changes of the biological matter.Here, the impedance spectrum of EBI manifests as Ż(jω)=F(V z (t))/F(I exc (t)), i.e., the Fourier Transform of the response voltage determines the impedance spectrum of the object one-toone thanks to the predetermined parameters of the excitation signal.
Usually, in theoretical considerations and simulations, the bioimpedance is substituted by a certain RC-circuit.Naturally, the accuracy of such approximation depends on the number and configuration of components.In Fig. 2a, a 5-element circuit is shown, where R 0 corresponds to the extracellular resistance r ext of a tissue, while R 1 , C 1 , R 2 and C 2 stand for the intracellular parameters.The respective Bode diagram has two real poles f p1 , f p2 , and two zeros f z1 , f z2 (Fig. 2b), spread typically over the frequency range from some kHz up to several MHz (Nebuya et al 1999;Pliquett et al, 2000;Grimnes & Martinsen, 2008).
The Laplace transform of the 5-element EBI (Fig. 2a) can be expressed as (Paavle et al 2008 ) For adequate estimation of the bioimpedance, the spectrum of excitation signal should cover the frequency range of the object Ż as much as possible.For that reason, several types of broadband excitation signals (e.g., Maximal Length Sequence (Gawad et al, 2007), multisine (Sanchez et al 2011)) are preferable.However, in this chapter, we will focus on the use of chirp excitation due to several advantages of it: wide and flat amplitude spectrum (Nahvi & Hoyle, 2009) together with independent scalability in the time and frequency domain.

Variety of chirps
A sine-wave based chirp with current phase θ(t) can be described mathematically as where A is the amplitude, and f(t) = (dθ(t)/dt)/2π is instantaneous frequency of the chirp signal.
A definite class of chirps has the instantaneous frequency f(t), which changes accordingly to some power function of the n th order (power chirps).Their specific quantity is chirping rate β=(f fin -f st )/T ch n , where f st and f fin are the initial and final frequencies, respectively, and T ch is the duration of the pulse.Chirps of this type can be expressed as The instantaneous frequency of a chirp can be increasing (up-chirps, β>0) or decreasing (down-chirps, β<0) quantity.Besides, sometimes it is practical to generate chirps with a symmetrical bidirectional frequency change (bidirectional or double chirps).In this case, the actual duration of the pulse is 2T ch , and the sign of chirping rate alters at t =T ch , causing mirrored waveform of the pulse against that moment.
Fig. 3 sketches waveforms of different chirps with the equal pulse duration and almost equal frequency range.A very basic linear (n=1) chirp with f st =0 can be described in accordance with the expression (4) as V ch (t)=Asin(2πf fin t 2 /2T ch ).It is depicted in Fig. 3a.
Fig. 3b shows the waveform of quadratic (n=2) down-chirp.However, the rule of frequency change can be arbitrary.For example, in some specific measurements, excitation with the exponential chirping rate β =( f fin /f st ) t/Tch is appropriate (Darowicki & Slepski, 2004).The waveform of sinusoidal exponential chirps (see Fig. 3c) is described as Generation of a perfect sine-wave chirp requires quite complicated hardware, which can cause problems, especially in on-chip solutions.That is why so-called signum-chirps (known also as pseudo-, binary-or Non-Return-to-Zero (NRZ) chirps) are often implemented (Figs.3d and 3e; the latter one depicts a binary chirp with bidirectional change of frequency).This kind of chirps can be defined by the signum-function of respective sinewave chirps as V ch (t)=V sgn (t) = sign(V sin (t)), which have binary values +A and -A only.Rectangular waveforms simplify signal processing: both generation of excitation and processing of the response will be substantially simpler.Especially simple is calculating of the correlation function or deconvolution in the time domain -shifting and multiplication with a reference signal having only {+1, -1}, or {+1, 0, -1} values (Rufer et al, 2005).
Additional advantage of signum-chirps is their unity crest factor (i.e., the peak amplitude ratio to the root-mean-square (RMS) value of the signal), and major energy compared with the sine-wave chirps of the same length.Unfortunately, using of rectangular signals causes the worse purity of the spectrum due to the accompanying higher harmonic components.

Short-time chirps (titlets)
Commonly, referring to chirps, signals of many cycles (rotations by 2π of the chirp generating vector) are considered (multi-cycle chirps).It does not need to be so, and the chirps with a single cycle or even less (θ(T ch ) ≤ 2π) can be generated and used, too.For such kind of ultra short chirps, the neologism "titlets" or Minimal Length Chirps (MLC) have been used (Min et al, 2011a).
The titlets can be very effective excitation signals in applications, where broadband excitation is necessary, but the minimal power consumption or extremely fast measurement are required at the same time.Because of the prospective use of titlets, a special attention to their properties will be paid below.The diagram in Fig. 4 explains the forming of a singlecycle linear sine-wave chirp.
non-uniformly rotating phasor Similarly, considering (5), one can show that for exponential chirps the duration of a pulse is ( ) Simulated waveforms of some typical titlets with L ≤ 1 are shown in Fig. 5.

Signal analysis and bioimpedance measurement by using Fourier Transform
Spectral analysis is the irreplaceable method as for the study of the behavior of excitation signals in the frequency domain as well as for processing the response signal in the bioimpedance measurements.
It is obligatory to remember that in the modeling and processing of chirp signals the selected sampling rate f s must follow the Nyqvist criterion over the whole frequency range, i.e., the condition f s ≥ 2f fin must be fulfilled.The sampling rate determines the total number of samples N ch during a chirped pulse as N ch = T ch f s .Yet, the simulation or real processing time T tot can be longer than T ch .In this case, the zero padding is used at t>T ch .The length of the input array for the Fast Fourier Transform (FFT) is N=T tot f s ≥ N ch , and it determines the frequency resolution Δf (difference between two successive frequency bins) of the FFT processing, whereas Δf = f s /N according to the uncertainty principle (Vaseghi, 2006).The acquired positive frequency range of the FFT processing in Hz is 0, …, NΔf/2 by steps of Δf.

Principles of bioimpedance measurements
There are several methods to calculate the frequency spectrum of the bioimpedance response.A typical way is the implementation of two FFT-channels for the response and excitation signals separately (Min et al, 2011a).Nevertheless, in the following discussion we will focus mainly on the structure and modeling of a specific bioimpedance measurement system, which incorporates the cross-correlation procedure together with following FFTprocessing (Paavle et al, 2008, Min et al, 2009, see Fig. 6).).In the latter case, the system works as the matched filter, enabling detection of mismatches between the impedance vectors Ż and Ż ref more precisely together with somewhat better noise reduction.The source noise n s (t) and object noise n z (t) can be taken into account in simulations as shown in Fig. 6.
In general, calculation of the CCF proceeds as ( ) where τ is a variable delay (lag) and the overline denotes averaging.
Using broadband chirp excitation, the waveform r xy (τ) of CCF is similar to the sin(x)/x type sinc-function affected by the nature of V z (t) and is strongly compressed close to τ = 0.
The CCF includes both the information about the amplitude level and phase shift of the bioimpedance vector.Thanks to this fact, only a single FFT-block is necessary to obtain the complete information about the object.Otherwise, when using the direct Fourier Transform of the response, another FFT channel is required to establish the basis for phase evaluations.
In accordance with the cross-correlation theorem (generalized Wiener-Khinchin theorem), the correlation function is equivalent to the inverse Fourier Transform of the complex crosspower spectrum density P xy (f) with the magnitude and phase spectral components

Modeling of the measurement system
Accordingly to the main architecture of measurement system (Fig. 6) and model of EBI (Fig. 2a), a special PC-model (Fig. 7) was developed for verifying the theoretical conceptions (Paavle et al, 2008).
where k = 0 … N ch -1, and N ch is the number of samples per length of the chirp pulse (McGhee et al, 2001).
Chirp generator If the uniform delay step is Δτ, then for the overall delay interval (observation time) t obs =τ max -τ min (τ min ≤τ m ≤ τ max ) the number of necessary computing cycles is M + 1 with M = t obs /Δτ.As a result, we acquire an array of correlation values r xy [m] with m = 0 … M. Supposing that the selected Δτ and τ min are integer multiples of the sampling interval 1/f s , the m th delay expresses as τ m = f s (τ min +mΔτ).
In the following Fourier Transform, the array r xy of length M+1 represents the input stream for the FFT-block, while in accordance with the uncertainty principle (see above), the frequency resolution of spectra becomes Δf =1/t obs over the frequency range f = 0 … 1/(2Δτ).For some specific applications, these properties must coincide with minor power consumption and with shortness of signal to ensure quick measurement.It is obvious that long chirps assure the more fitting shape of spectra, except a certain rippling caused by the Gibbs effect.The latter can be suppressed by using a kind of windowing in the time domain (often a boxcar-type window function is used for this purpose), but in this case, the total energy of the signal decreases.According to the Parseval's theorem, the total energy in the frequency domain is the same as the energy in the time domain (Vaseghi, 2006).For chirps with the finite length, a certain part of signal energy falls outside the useful (generated) bandwidth caused by higher frequency components.For that reason, it is necessary to distinguish the total energy of a generated chirp pulse, which varies proportionally with T ch , and the useful energy E exc , falling inside the chirp bandwidth B exc .Typically, E exc < E tot , while absolute values of both quantities depend on the chirp length, waveform and spectral nature.To characterize the percentage of the useful energy, we employ the term of energy-efficiency as δ E =E exc /E tot .For sine-wave chirps with L→∞ the ratio δ E →1.   4) or ( 5), the analytical description of the chirp energy becomes very complicated.So it is reasonable to calculate the respective quantities numerically.Nevertheless, a good approximation of total energy and energy-efficiency of arbitrary chirp signal can be obtained using the results from the Fourier Transform as follows (Paavle et al, 2010): where |V(f i )| is the value of the amplitude spectrum at i th frequency bin, N st and N fin are the numbers of frequency bins, corresponding to the f st and f fin , respectively.N max is the total number of frequency bins, and the divisor in ( 12) corresponds to the total energy E tot .Of course, this method enables to calculate the partial energy for any frequency interval inside the range of 0 to N max Δf.

Spectral features and energy of chirps
The curves in Fig. 11 present the dependence of the energy and energy-efficiency of linear sine-wave chirps on the number of chirp cycles graphically.A selection of power and energy parameters for different chirps and titlets, obtained from the FFT-processing by using the expression ( 12), are converged into Table 1, where the average power ch 2 avg ch tot ch ch 0

() / T PV t d t E T T == ∫
presumes the 1 kΩ load.

Binary signum-chirps
A specific feature of signum-chirps is the gradually decreasing amplitude spectrum by step width of 2B exc as shown in Fig. 12 (Min et al, 2009).Let us analyze this phenomenon.
The fundamental harmonic of a regular rectangular signal with amplitude A 1 = (4/π) A has a root-mean-square (RMS) value 1 /2 4 / ( 2 ) AA =π , energy E 1 = (A 1 2 /2) T ch = (8/π 2 ) A 2 T ch , and power W 1 = E 1 /T ch = (8/π 2 ) A 2 , which creates a constant value power spectral density (PSD) w 1 =W 1 /B exc , V 2 /Hz, within the bandwidth B exc of fundamental harmonic: ( ) The amplitudes of k th higher harmonics are of A k = A 1 /k.Therefore, the power of every k th higher odd harmonic (k =3, 5, 7, …) is equal to W 1 /k 2 , being spread over the frequency range B k =kf fin − kf st =kB exc .Due to the fact that higher harmonics have k times wider bandwidth than the fundamental (first) harmonic has, the PSD for higher harmonics can be expressed through (12) as w k = (w 1 /k 2 )/k = w 1 /k 3 (Min et al, 2009).
The total power of generated signal V ch (t) is gradually distributed over its whole frequency range, theoretically from f st to ∞, see Fig. 12.The PSD p h of every gradual h th level (h= 1, 2, 3, 4, …) of the spectrum, beginning from the first one p 1 , is the sum of power spectral densities of the fundamental and higher harmonics, w 1 and w k , which contribute into the given level h.Into the PSD p 1 of the first level contribute all the signal components p 1 = w 1 + w 3 + w 5 + w 7 + …, but into the second level mere the higher harmonics p 2 = w 3 + w 5 + w in which i is an integer beginning from h, that is: i = h, (h+1), ( h+2), (h+3), …, (h+m), where m is the number of higher odd harmonics taken into account.
The partial power P 1 = p 1 B exc = P exc of the first level (h=1) is the useful excitation power ( ) The power P out falling outside B exc is lost.The lost power can be found as a summed up power of higher harmonics within the next levels of spectrum (h= 2, 3, 4, etc., see Fig. 10): ( ) () where W 1 = w 1 B exc in ( 15) and ( 16) is the power of fundamental harmonic.As the minimal frequency for every level is kf st , then the equations (14-16) are approximate ones, which become absolutely exact only for f st =0.In practice, when the excitation bandwidth B exc exceeds one decade significantly, e.g., f fin / f st > 30, these equations are exact enough for engineering calculations.
The total chirp energy is E tot = E exc + E out = ( P exc +P out ) T ch and the useful excitation energy in it is E exc = P exc T ch .From ( 15) and ( 16) follows the role of useful energy δ E = E exc /E tot .
Considering f st =0, the ratio for the energy-efficiency expresses as Taking into account all the higher harmonics (m→∞ and (max h)→∞), the sums in (17) will obtain limit values, which can be found via Riemann zeta function (Dwight, 1961): www.intechopen.com Riemann's mathematics gives us the following limit values for the sums in ( 17): where ζ(3)≈1.202 is termed the Apéry's constant (Dwight, 1961).Now, using equations ( 19) and ( 20 When f fin / f st >> 30 and m >> 1, the value of δ E = 0.852 and the spectrum has practically uniform energy distribution over the B exc .The relative decrement of every h th step of the averaged amplitude can be calculated from ( 14) as δ h = (E exc /E h ) 1/2 =A 1 /A h .The values for some lower order levels are δ 2 ≈ 4.51, δ 3 ≈ 8.44, and δ 4 ≈ 12.47 − see the approximations given in Fig. 12.For the next levels, a rough approximation δ h ≈ 4(h−1) +1 is appropriate for engineering calculations (Min et al, 2009).The shape of amplitude spectrum becomes more complicated, if the initial frequency f st ≠0.
Though, the spectrum retains the gradual character, the levels do not attenuate monotonically.An example of this kind of spectrum is shown in Fig. 13a, where f st = 50 kHz and f fin =100 kHz.In general, the height and location of spectral levels can be various, depending on the ratio between f st and f fin .
The forming of spectra as the sum of its harmonic components is explained by the draft diagram in Fig. 13b.Fulfillment of the condition 3f st > f fin produces the frequency area B zero = 3f st -f fin next to the excitation bandwidth, where the spectrum has almost zero amplitude.In this particular case, the spectral density within the excitation bandwidth is determined only by the fundamental harmonics from f st to f fin .It can be shown, considering the Eqs.( 13) -( 16), that here the energy-efficiency δ E =8/π 2 .

Ternary chirps
Caused by the absence of particular harmonics (see Sect. 2.2.1), the position and span of amplitude levels of ternary chirps (RZ chirps) differ from the levels of the respective binary chirp.Fig. 14a enables to compare the amplitude spectra of several rectangular chirps with L=1000 at various durations of the zero-level state α d (shortening).
Let us pay attention to the changed stretch of the amplitude levels.For example, in the case of α d =30º, the average amplitude of 2 nd level (h=2) stretches uniformly from f fin up to the 5f fin .As the 3 rd harmonic are absent, this 2 nd level is formed by the amplitudes of the 5 th and higher harmonics up to the frequency 5f fin .Average power of ternary chirps is less than it of respective binary chirps, but the percentage of their useful energy is surprisingly high -over 90% mostly (see also Table 1).Actually, it is possible to generate ternary chirps with any value of the zero-state, but only a few values produce the removing of particular harmonics (see Eq. 6 in Sect.2.2.1).Nevertheless, we can analyze the energy properties of arbitrary ternary chirp.Employing the equation ( 12), the maximal energy-efficiency δ E ≈93.4% was observed at α d ≈22.5º, which does not cause disappearance of any odd harmonic component.Dependence of the energyefficiency on the shortening α d is plotted in Fig. 14b.

Windowing of chirps
In the use of some short-time chirps, the problem of flatness of the amplitude spectrum arises (see Fig. 9b).For the single-cycle sine-wave chirp, the maximum overshoot of normalized voltage spectral density is about +7.4 dB inside the chirp bandwidth.To improve the flatness of spectrum, a kind of additional windowing of chirp pulses should be used (in fact, every finite chirp pulse can be dealt as one inside the rectangular window, however we consider this case as the unwindowed one).As a rule, the windowing accompanies with some loss of total energy and power of signals, but still the δ E can be rather high.On the other hand, frequently the spectral density attains steeper drop-off outside the chirp bandwidth due to windowing.Hence, the optimal choice of windows presumes a certain trade-off always.
In this work, several typical (Hanning, Hamming, Nuttall, etc.) and some specific window functions F win (t) were under study (Barsoukov & Macdonald, 2005).For example, a convenient shaping of amplitude spectra was achieved by implementing window functions F win (t)=sin 2 (πt /T ch ), which can be dealt as a particular case of the Tukey window with the squared-sine lobes and with the tapering time T ch /2.Almost perfect shaping was attained using a non-symmetrical windowing in the form of F win (t)=(t/T ch ) a with the selectable exponent a (usually a =2…8) (Paavle et al, 2010).Some windowing results both in the time and frequency domain are shown in Fig. 15.Usually, the deviation about ± 3 dB of the amplitude spectrum inside the B exc is considered as satisfactory for spectral flatness.There are several ways to achieve such the requirement.
Using the non-symmetric windowing makes the deviation even less than ± 1 dB accessible (see Fig. 15b), but the quantity of useful energy reduces noticeably due to the plain slope of the spectrum at f >f fin : for the case, as shown in Fig. 15, the energy-efficiency δ E =73.1% was observed in simulations.Substantially higher δ E was achieved by using of Nuttall (δ E =92.3%) and Hanning (δ E =95.3%) windows, but the deviation of amplitude spectra inside the B exc is from 0 to -6 dB and from +0.2 to -6 dB, respectively (Paavle et al, 2010).The impact of the squared-sine windowing was almost the same as it of the Hanning window.
As a rule, windowing of the generated signal demands some additional power consumption and reduces the total energy-efficiency.However, in the case of very short signals this drawback can be overcome using the look-up table, which stores the externally calculated values of windowed titlets and which is loaded into the FPGA (Min et al, 2011a).

Simulation examples
The following examples refer to the linear and quadratic sine-wave chirps and are the results of simulation using the modeling structure in Fig. 7.In all the cases, the chirp excitation from f = 0 to 100 kHz and amplitude A=1 V was implemented.The basic model of EBI was a 5-element impedance Ż 1 , consisting of R 0 =1 kΩ, R 1 =200 Ω, R 2 =100 Ω, C 1 =30 nF, C 2 =20 nF and having the following corner frequencies in the Bode diagram (see Fig. 2b): f p1 =2.9, f z1 =26.5, f p2 =45.2, and f z2 =79.6 kHz.Intensive fluctuation of spectra, caused by the concurrent higher harmonics, can be more disturbing for signum-chirp excitation (see Figs. 10b,12,and 13a).This problem can be diminished by proper selection of correlation parameters, e.g., the shorter t obs with rougher frequency resolution smoothes spectrograms essentially (Paavle et al, 2008).

V(t)
Time, ms   For interpreting the variations of frequency responses, the phase spectra should be preferred.First, it is less affected by the spectral fluctuations.Secondly, the phase spectra enable to distinguish changing of object parameters somewhat more clearly than the amplitude ones -compare the respective curves in Figs.16c and 16d.Nevertheless, substantially better resolution of the EBI for the diagnosing purposes can be achieved by using matched filtering, where instead of Ż ref =1 in the reference channel, a predetermined Ż ref ≈Ż x is used, approximately equal to the impedance Ż x under study (see Fig. 6).Joining matched filtering together with analysis of phase spectra allows applying the full scale of some degrees only and permits very good sensitivity for detection of small changes of the object parameters.
Simulation results in implementing of the bidirectional quarter-cycle quadratic chirp excitation for detection of tiny deviations of the EBI vector are presented in Fig. 17  Another advantage of the matched filtering is reducing the impact of additive noise, especially the impact of source noise n s (t).Fig. 18a depicts the noisy input signal, where the generated chirp is affected by the Gaussian noise with the zero mean and variance σ N 2 -see Fig. 6.The source noise causes some error in the phase spectra at the higher frequencies (Fig. 18b).This error depends on the deviation of component parameters of the equivalent circuit (here C 1 was changed), and on the signal-to-noise ratio (SNR=P avg /σ N 2 ), but as a rule, the error remains small even at considerably high level of the source noise.As the noise from object (n z (t) in Fig. 6) affects the cross-correlation procedure nonsymmetrically, then suppression of noise effect is of importance and serves special attention even at the higher SNR values (Min et al, 2011b).

Conclusion
The advantages of using chirps as excitation signals in bioimpedance measurement are their wide and almost flat amplitude spectrum together with the independent scalability both in the time and frequency domain -we can choose the frequency range and duration of the excitation pulse almost independently from each other.These features enable to accommodate the generation of excitation signals with the expected properties of the object to be estimated comparatively simply.
It was shown that the shortening of chirp pulses retains the general benefits of chirps -their flat amplitude spectrum within the predetermined frequency range, which permits their implementation in energy-efficient measuring instruments.More than 90% of the generated excitation energy falls into the desired bandwidth even in the case of very short excitation pulses required for providing ultra quick measurement and analysis of dynamic objects.These requirements are obligatory for investigation of objects with rapidly changing parameters (e.g., for identification of fast moving bioparticles such as cells and droplets in high-throughput micro-fluidic systems), and in the devices, in which the low power consumption is important (e.g., wearable units and medical implants).However, shortening of chirp excitation and measurement time should not be excessive.It must be as short as possible to avoid significant impedance changes during the Fourier analysis, but as long as possible to enlarge the excitation energy and to obtain a better signal-to-noise ratio.www.intechopen.com The described measurement method, based on the sequentially performed cross-correlation and Fourier Transform, enables the joint time-frequency spectral analysis of bioimpedance.
As the time domain cross-correlation function includes the full information about the complex vector Ż(ω), only a single FFT-block is necessary for the Fourier analysis and the complete frequency domain determination of the object.Moreover, the cross-correlation procedure assures substantial suppression of noise introduced by the object.
Applying the modification of system architecture with the cross-correlation based matched filtering enables to avoid the impact of noise in much higher degree then a simple crosscorrelation.In addition, such the matched filtering permits to increase the sensitivity of measurement, especially through the analysis of phase spectra.Identification and interpretation of relative deviations of the object parameters through the detected tiny phase shifts is of great importance, which in turn can give valuable information about the state and processes in the biological objects.

Fig. 6 .
Fig. 6.Basic structure of the measurement system In such a system, the cross-correlation function (CCF) is calculated: (a) between the response V z (t) and excitation V ch (t) at the unity gain in the reference channel (V ref (t)= V ch (t), Ż ref =1), or (b) between the response and predefined reference signal V ref (t) (reference channel includes the known impedance Ż ref ≠1).In the latter case, the system works as the matched filter, enabling detection of mismatches between the impedance vectors Ż and Ż ref more precisely together with somewhat better noise reduction.The source noise n s (t) and object noise n z (t) can be taken into account in simulations as shown in Fig.6.

Fig. 7 .
Fig. 7. Model of the measurement system The cross-correlation function is calculated by the cumulative adder for every m th lag τ m as Next, let us take a look at the results of direct Fourier Transform of distinct chirp pulses.The direct FFT of different excitation signals enables to compare their amplitude and power spectra together with estimation of the energy properties keeping in view the requirements for the bioimpedance measurement.Important and desired spectral properties of excitation signals to improve the quality of wideband measurement are: • flat amplitude spectrum with minimal fluctuation (ripple) together with the absence of overshoots inside the generated (excitation) bandwidth B exc =f finf st ; • steep drop-down of the amplitude spectrum outside the bandwidth B exc ; • maximal energy-efficiency, i.e., the ratio between the energy lying within the generated (excitation) bandwidth B exc and total energy of the signal.

Fig. 8 .
Fig. 8. Voltage spectral density of linear chirps 0…100 kHz (Δf=50 Hz): (a) linear scaling at the unity amplitude of chirp pulses; (b) normalized logarithmic scaling Fig. 8 depicts the amplitude spectra of linear sine-wave chirps of different length in the linear and in the normalized logarithmic scales.Thereby, the base of normalizing is the amplitude |V(f 0 )|, i.e., the magnitude at the lowest frequency bin f 0 of the Fourier Transform.It is obvious that long chirps assure the more fitting shape of spectra, except a certain rippling caused by the Gibbs effect.The latter can be suppressed by using a kind of windowing in the time domain (often a boxcar-type window function is used for this purpose), but in this case, the total energy of the signal decreases.

Fig. 9
Fig.9shows the frequency run and respective amplitude spectra of different titlets.It appears that the desired shape of the spectrum, including flatness, satisfying drop-down outside the bandwidth (slopes from -20 dB/dec to -80 dB/dec were observed) and admissible overshoots can be achieved by a proper selection of the type and length of titlet.Moreover, additional correction and shaping of spectra can be attained by windowing of titlets in the time domain (see below in Sect.3.2).In Fig.10, waveforms and amplitude spectra of some rectangular chirps (see Sect. 2.2.1) with L=10 are shown.Naturally, for this kind of chirps, the rippling of spectra is noticeable.

Fig. 16 .
Fig. 16.Examples of the time-frequency analysis of the EBI (sine-wave chirp excitation with L=1000): (a) waveforms of excitation and response; (b) normalized correlation functions; (c, d) cross-power spectral density and phase spectra at variations of the EBI components Fig. 16 demonstrates simulation results in the case of multi-cycle chirp excitation with pulse duration of T ch =20 ms (1000 cycles).Fig.16ashows the response voltage V z (t) on the background of the non-windowed excitation signal.Fig.16bshows cross-correlation functions in a stretched time scale at the Ż= Ż 1 , and at the unity value impedance Ż= Ż 0 =1 (in fact, the latter case corresponds to the autocorrelation function of the excitation signal).Results of the Fourier Transform of the CCF are shown in Figs.16c and 16d, where one can watch the swing of spectral curves, if the capacitive components of EBI change by +50%.In addition, Fig.16cdemonstrates, how the boxcar-type windowing (here squared-sine Tukey window with the tapering time of 0.05T ch ) reduces the influence of the Gibbs effect.

Fig. 17
Fig. 17.The use of double quarter-cycle quadratic chirp excitation (L=¼, n=2, T ch =2×7.5 μs, f st =0, f fin =100 kHz): (a) signal waveform; (b) variation of instantaneous frequency; (c) amplitude spectrum of the chirp pulse; (d) phase spectra at matched filtering (Ż ref = Ż 1 ) and at changing of the EBI components +5% Fig. 17.The use of double quarter-cycle quadratic chirp excitation (L=¼, n=2, T ch =2×7.5 μs, f st =0, f fin =100 kHz): (a) signal waveform; (b) variation of instantaneous frequency; (c) amplitude spectrum of the chirp pulse; (d) phase spectra at matched filtering (Ż ref = Ż 1 ) and at changing of the EBI components +5% . The waveform of excitation pulse and the corresponding frequency run are shown in Figs.17a and 17b, respectively.Fig. 17c depicts the normalized amplitude spectrum of the excitation signal and Fig. 17d manifests deviation of phase spectra, when the value of single components increases +5% over their initial quality.

Fig. 18 .
Fig. 18.Linear single-cycle chirp (T ch =20 μs) at the presence of additive source noise: (a) noisy input signal; (b) spectra of phase differences at 1.05C 1 and at different input signalto-noise ratio (SNR) levels

Table 1 .
Energy and average power of different chirp pulses with f st ≈0 and f Type of chirp (titlet)L f st , Hz T ch , μs P avg , mW E tot , nJ δ E , fin =100 kHz Considering expressions (11) and (