# Discrete-level broadband excitation signals: binary/ternary chirps/Diskretaus lygio placiajuosciai dvigubi arba trigubi cirskimo suzadinimo signalai.

IntroductionChirps (fast frequency sweeps) and chirp-like signals have become popular excitation signals in measurement techniques, particularly as the stimulus in bioimpedance measurement [1, 2]. The main advantage of chirps is the independent scalability in time and frequency domains: the duration of a pulse (usually short) and its frequency range (usually wide) can be adjusted separately. This feature is important in biomedical investigations for quick monitoring of the changing state of living tissue and cells.

The classical description of chirps presumes a sine-wave signal, instantaneous frequency f(t) of which changes linearly in time. Consequently, the current phase [theta](t)=[integral]f(t)dt of this chirp is changing in proportion to the square of time. The waveform of chirp expresses as [V.sub.ch](t)=sin([theta](t))=sin(2[pi][integral]f(t)dt).

Generation of perfect sinusoidal chirps is an intricate task. Therefore, for simplifying generation and signal processing, applying chirps of the rectangular waveform with a small number of discrete amplitude levels should be preferred, e.g., for miniaturized measurement devices. Below we deal with such kinds of discrete-level chirp-like signals, known also as pseudo-chirps [2].

Bioimpedance measurement with chirp stimulus

The impedance of arbitrary biological matter or bioimpedance (BI), shortly, can be characterized by its electrical equivalent (EBI). The EBI is described mathematically by the frequency-dependent complex vector [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Adjusting the parameters of excitation current [I.sub.exc] through the bio-object, and measuring the response voltage [V.sub.z], we can find out the impedance spectrum of it as [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where F denotes the Fourier Transform.

A possible architecture of the system for the EBI measurement is presented in Fig. 1 [3]. It executes the cross-correlation procedure (Corr) between the response and reference signals, and thereupon calculates Fourier transform (FFT) of the cross-correlation function, i.e., the measurement of the unknown bioimpedance Z (or any complex impedance) is accomplished in succession Corr([V.sub.z], [V.sub.ref])[??]FFT[??]{Re(Z(j[omega])); Im(Zj[omega]))}. The result is a broadband impedance spectrum, which characterizes the difference between the EBI under study and the predetermined reference impedance [Z.sub.ref].

As a rule, the informative frequency bandwidth of the BI lies within the range from some kHz up to several MHz [1, 4]. The benefit of chirps is the simplicity of adjusting the sufficiently wide excitation range [B.sub.exc], which covers the frequency bandwidth of the object under measurement.

[FIGURE 1 OMITTED]

In general, the manner of frequency change of chirps can be various: by a power function, exponentially, piecewise linearly, etc. We will pay the main attention to the first mentioned ones or so-called power chirps, the instantaneous frequency of which is changing according to the n-th order power function with arbitrary (incl. fractional) power n. This kind of sine-wave chirps with the amplitude A can be described mathematically as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where [beta]=[B.sub.exc]/[T.sub.ch] characterizes the rate of frequency change (chirping rate), [T.sub.ch] is the pulse duration, and [B.sub.exc]=[f.sub.fin]-[f.sub.o] with the final[f.sub.fin] and initial [f.sub.0] frequencies.

Usually, chirps of many cycles (rotations of chirp-forming vector on the complex plane) are considered. However, it was shown in [5] that chirps of a single cycle or even less ([theta]([T.sub.ch]) [less than or equal to] 2[pi]) can have the practical importance, if the very low power consumption or the extremely short measurement time is required (implantable devices, lab-on-chips). The end frequencies, time duration and number of cycles of a chirp are strictly related. For the n-th order power chirps this relationship can be expressed as [5]

[T.sub.ch] = L(n + 1))(n[f.sub.0] + [f.sub.fin]), (2)

where L is the number of cycles.

The basic discrete-level pseudo-chirp is the binary one (called also signum-chirp or Non-Return-to-Zero (NRZ) chirp), waveform of which represents the signum-function of the respective sine-wave chirp as A-sgn([V.sub.ch](t)), thus having values of { [+ or -] A} only - see Fig. 2a.

[FIGURE 2 OMITTED]

Another important class of pseudo-chirps is ternary (trinary) chirps or Return-to-Zero (RZ) chirps with the possible values of {+A; 0; -A}--see Fig. 2b.

Unlike sine-wave chirps, the spectrum of the discretelevel ones is rippling intensively. Besides, the percentage of their in-band energy is somewhat less at the equal amplitude and duration of the pulse. In practice, these drawbacks should be compensated by the simplicity of the signal processing. An additional benefit of binary chirps is their unity crest factor, and consequently, major energy compared with the sine-wave chirps of the same length.

Next, let us examine some specific spectral features of both the mentioned modifications of linear pseudochirps (n=1).

Binary chirps with [f.sub.0]= 0

A specific feature of binary chirps is the gradually decreasing amplitude and power spectra by plateaus h (h=1, 2, ...; Fig. 3). This effect and the energy distribution of binary chirps along the frequency spectrum, considering [f.sub.0]=0, was discussed by the authors in [6] as follows.

The fundamental harmonic (k=1) of a rectangular signal with levels [+ or -]A has the root-mean-square (RMS) value 4 A /([pi][square root of 2]) and (presuming the unity load) the power W1= (8/[[pi].sup.2])[A.sup.2].

In the case of rectangular-wave chirps the fundamental harmonics from [f.sub.0] to [f.sub.fin] are spread over the whole bandwidth [B.sub.exc], creating a constant value power spectral density (PSD) [w.sub.1]=[W.sub.1]/[B.sub.exc], [V.sub.2]/Hz. The power of every [k.sup.th] higher odd (k=3, 5, .) harmonic Wk=W1/k2 is spread over the frequency range [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] with PSD of [w.sub.k]= [W.sub.1]/([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) - see Fig. 3a.

All the odd harmonic components contribute a certain part of their power for the first plateau (h=1) of the power spectrum with the width of [B.sub.exc] (excitation power [P.sub.exc]). It can be expressed as the sum

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

[FIGURE 3 OMITTED]

The power of every next plateau (h>1) with the width of 2[B.sub.exc] is formed by the harmonics from (2h-1)th and higher, which expresses as (Fig. 3b)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

The total power of all spectral plateaus outside the generated bandwidth is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

For the first spectral plateau, the PSD [p.sub.1]=[P.sub.1]/[B.sub.exc], while for all the next plateaus [p.sub.h]=[P.sub.h]/(2[B.sub.exc]). The values of [p.sub.1]/[p.sub.h]=2[P.sub.1]/[P.sub.h] can be calculated by (3)-(5)--some of the results are shown in Table 1 and in Fig. 3c.

The total energy of the generated chirp signal is [E.sub.tot]=([P.sub.exc]+[P.sub.out])[T.sub.ch]. However, it is essential to know not only the full energy, but also the amount of the energy [E.sub.exc], within the generated bandwidth (useful energy). Energy-efficiency of the generated excitation signal can be expressed as the ratio of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Using (3)-(5) yields [6]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

Where [zeta] denotes Riemann's zeta-function

The theoretical percentage of the useful energy by (6) is [[delta].sub.E][approximately equal to] 0.853. In practice, the energy-efficiency of a signal can be evaluated from the results of FFT-processing as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

where [absolute value of [V.sub.ch]([[phi].sub.i])] is the value of the amplitude spectrum at [i.sup.th] frequency bin [[phi].sub.i], [N.sub.0] and [N.sub.fin] are the numbers of frequency bins, corresponding to the [f.sub.0] and [f.sub.fin], respectively. [N.sub.max] is the total number of frequency bins of the FFT processing.

Some specific values for the chirps with length of 1000 cycles are presented in Table 2.

Binary chirps with [f.sub.0] >0

When the initial frequency of a chirp [f.sub.0][not equal to]0, then the shape of spectra becomes more complicated. Though the spectrum retains the gradual character, the shrinking of plateaus is not always monotonic. It depends on the ratio between [f.sub.0] and [f.sub.fin]. In a certain case, one or more frequency sub-bands [B.sub.z] with zero magnitude appear into a spectrum.

The forming of such a spectrum as the sum of its harmonics is explained by the draft diagram in Fig. 4a.

[FIGURE 4 OMITTED]

It is evident that fulfillment of the condition 3[f.sub.0] >[f.sub.fin] produces the frequency area [B.sub.z1] = 3[f.sub.0] -fin next to the excitation bandwidth, where the spectrum has almost zero magnitude. In this specific case, the spectral density within the excitation bandwidth is determined only by the fundamental harmonic components from [f.sub.0] to [f.sub.fin]. In general, the frequency area above [f.sub.fin] can include more than a single zero sub-band, which separate spectral blocks of particular harmonics. Fulfilling of the condition

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

causes the zero-valued area from the frequency (2i-1)[f.sub.fin] up to (2i+1)[f.sub.0], i.e., [B.sub.z,I] =(ffn+[f.sub.0]) -2i([f.sub.fin]-[f.sub.0]) (Figs. 4b, 4c).

The condition (8) with i=1 has the more practical importance, assuring moderately wide, but distinctly separated excitation bandwidth--see simulated and smoothed spectrum in Fig. 4b, where [f.sub.0]=50 kHz and [f.sub.fin]=100 kHz. This condition with i [greater than or equal to] 1 simplifies evaluation of the energy-efficiency because the energy within the excitation bandwidth is always equal to the energy of first harmonic components from [f.sub.0] to [f.sub.fin] only (see the discussion in the previous chapter)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

Thus, in the case of binary chirps with 3[f.sub.0] >[f.sub.fin] the theoretical energy-efficiency is constantly [[delta].sub.E][approximately equal to]0.81, i.e., it is equal to that of the common rectangular signal. Let us note that at i[right arrow][infinity] is equivalent to the [f.sub.0]=[f.sub.fin].

Ternary chirps

A method to improve the spectral properties of rectangular signals (incl. rectangular chirps) is shortening their duty cycle by a certain amount, thus modifying the binary signal to a ternary one (see Figs. 1b and 4) [7]. The shortening can be characterized by the relative time t/[T.sub.ch] or by the equivalent phase angle [alpha] (for binary signals [alpha] =0) per a quarter of a period (or cycle), within of which the value of signal remains on the zero state (Fig. 5).

[FIGURE 5 OMITTED]

Let us explain the effect of shortening by Fourier series, which expresses for the rectangular-wave signal with the magnitude [+ or -]A as the sum of odd harmonics

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

As it follows from (10), the kth harmonic (k=2i-1 with i=1, 2, 3, ...) is absent from the series F([omega]t), when k[alpha] = [+ or -](2i+1)[pi]/2, and i=0, 1, 2, ..., because of cos(k[alpha])=0. For example, [alpha] =[pi]/6, causes removing of the [(3+6i).sup.th] harmonics, [alpha]=[pi]/10 removes the [(5+10i).sup.th] harmonics, etc. This is valid for the rectangular chirps, too, but in this case the time duration of every quarter-cycle is changing in accordance with the instantaneous frequency (Fig. 5b). Hence, in the generation of ternary chirps [V.sub.ch](t)=sgn(sin([theta](t))) one must follow the current phase: if [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], then [V.sub.ch](t)=0 (Fig. 6). Assuming numerical processing with sampling frequency [f.sub.samp], the current phase [[theta].sub.i] (i=0, 1, 2, ..., [T.sub.ch][f.sub.samp]-1) in degrees according to (1) expresses as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

[FIGURE 6 OMITTED]

Fractional-power (sub-linear) chirps

The previous examples indicated almost constant (on the average) spectral density within the excitation bandwidth (Figs. 3, 4, 6), but in biomedical studies it is sometimes practical to apply the stimulus with the spectral density rising with frequency. This effect can be achieved by using of sub-linear power chirps with the exponent n, slightly less than unity, for example, n=0.8 to 0.9.

Conclusions

The chirp-wave excitation signal is highly convenient in bioimpedance measurements and impedance spectroscopy due to the possibility of fast and broadband identification of objects. Though the best measurement quality can be achieved by implementing the sine-wave chirp stimulus, the advantage of rectangular-wave signals is the simplicity of generation and processing of signals.

Aside the simpler hardware and software, another essential factor is energy consumption, which is extremely important in the development of miniaturized devices. That is why the energy-efficiency of rectangular-wave chirps were analyzed in this work. It was shown, that the useful energy of discrete-level chirps is sufficiently high, exceeding 90% of the generated energy for ternary ones.

Depending on the aim and practicable needs, the shape of spectra can be simply modified by changing the initial frequency and duty cycle of generated chirps or selecting an appropriate mathematical basis of the signal.

The drawback of rectangular-wave chirps is their intensely rippling spectra, which can disturb the interpreting of measurement results. A way to overcome this issue is to focus on the phase spectra, in which the effect of rippling is diminished largely. This opportunity is more closely treated in the earlier works of authors [3].

Acknowledgements

This work was supported by the EU through the European Regional Development Fund, Estonian target-financed project SF0140061s12 and by Enterprise Estonia through the ELIKO Competence Center.

Received 2012 03 19

Accepted after revision 2012 04 17

References

[1.] Barsoukov E., MacDonald J. R. Impedance Spectroscopy: Theory, Experiment, and Applications--New York: Wiley, 2005-595 p.

[2.] Pollakowski M., Ermert H. Chirp signal matching and signal power optimization in pulse echo mode ultrasonic nondestructive testing // IEEE Trans. Ultrason., Ferroelect. Freq. Contr., 1994--No. 41(5).--P. 655-659.

[3.] Min M., Paavle T., Ojarand J. Time-frequency analysis of biological matter using short-time chirp excitation // Proc. European Conf. on Circuit Theory and Design (ECCTD2011), 2011.--P. 585-588.

[4.] Grimnes S., Martinsen 0. G. Bioimpedance and Bioelectricity Basics--Elsevier-Academic Press, 2008.--471 p.

[5.] Paavle T., Min M., Ojarand J., Parve T. Short-time chirp excitations for using in wideband characterization of objects: an overview // Proc. 12th Biennial Baltic Electronics Conf. (BEC2010), 2010.--P. 253-256.

[6.] Min M., Paavle T., Annus P., Land R. Rectangular wave excitation in wideband bioimpedance spectroscopy // Proc. IEEE 4th Int. Workshop on Medical Measurements and Applications (MeMeA2009), 2009.--P. 268-271.

[7.] Min M., Kink A., Land R., Parve T, Ratsep I Modification of pulse wave signals in electrical bioimpedance analyzers for implantable medical devices // Proc. 26th Annual Intern. Conf. EMBS'2004.--IEEE, 2004.--P. 2263-2266.

T. Paavle, M. Min

Th. J. Seebeck Dept. of Electronics, Tallinn University of Technology, Ehitajate tee 5, 19086 Tallinn, Estonia, phone: +372 6201156, e-mail toivo@elin.ttu.ee

http://dx.doi.org/10.5755/j01.eee.122.6.1815

Table 1. The ratio of the PSD plateaus [p.sub.1]/[p.sub.h] [p.sub.1]/[p.sub.2] [p.sub.1]/[p.sub.3] [p.sub.1]/[p.sub.4] 20.3 71.2 155.2 [p.sub.1]/[p.sub.5] [p.sub.1/[p.sub.6] 273.4 424.9 Table 2. Energy and power of chirps (A=1 V, [R.sub.load]=1 k[OMEGA]) [E.sub.tot], [[delta]. Type of chirp P, mW [micro]J sub.e], % sine-wave 0.5 10 99.8 binary 1.0 20 85.1 ternary 18 [degrees] 0.8 16 93.1 ternary 30 [degrees] 0.66 13.2 92.1

Printer friendly Cite/link Email Feedback | |

Title Annotation: | SIGNAL TECHNOLOGY/SIGNALU TECHNOLOGIJA |
---|---|

Author: | Paavle, T.; Min, M. |

Publication: | Elektronika ir Elektrotechnika |

Article Type: | Report |

Geographic Code: | 4EXES |

Date: | Jun 1, 2012 |

Words: | 2716 |

Previous Article: | Generating unit tests for floating point embedded software using compositional dynamic symbolic execution/Vienetu testu generavimas slankiojo... |

Next Article: | Feature extraction and reduction of wavelet transform coefficients for EMG pattern classification/Pozymiu isskyrimas ir vilneliu transformacijos... |

Topics: |