Method and device for multichannel multifrequency analysis of an object

• US 8,773,151 B2
• Filed: 04/24/2007
• Issued: 07/08/2014
• Est. Priority Date: 04/24/2006
• Status: Expired due to Fees

• Current Assignee: Ou Eliko Tehnoloogia Arenduskeskus
• Original Assignee: Ou Eliko Tehnoloogia Arenduskeskus
• Inventors: Min, Mart, Annus, Paul, Kuusik, Alar, Land, Raul, Parve, Toomas, Ronk, Ants, Haapalainen, Antti
• Primary Examiner(s): NATALINI, JEFF WILLIAM
• Application Number:
  US12/298,389
• Publication Number:
  US 20100225303A1
• Time in Patent Office: 2,632 Days
• Field of Search: 324/707, 324/709, 324/76.24
• US Class Current: 324/707
• CPC Class Codes: A61B 5/053 Measuring electrical impedance or conductance of a portion of the body
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: G01R 23/00 Arrangements for measuring frequencies; Arrangements for analysing frequency spectra

Abstract

A method for multichannel multifrequency analysis of an object, applying a set of excitation signals to the object and sampling the response signal from the object, using uniform and non-uniform undersampling. Non-uniform sampling of the response signal is performed, i.e., the sampling is performed for two or more different frequencies in one observation time slot. Also, uniform sampling of the response signal is performed, i.e., the sampling of a signal, corresponding to one frequency, is performed for two or more channels within one observation time slot and then sampling the same signal for another frequency for two or more channels within the next observation time slot.

Claims

1. Method for multichannel multifrequency analysis of an object, comprising:

generating a set of excitation signals wherein all excitation signals are periodic signals with preselected frequencies;

applying said set of excitation signals simultaneously to the object;

receiving from the object, in response to said set of excitation signals, a response signal corresponding to said set of excitation signals, modulated by the object'"'"'s properties, said response signal comprising a set of components with different frequencies;

calculating a set of phase-related sampling instances based on said set of excitation signals, for each frequency of interest and for at least one predetermined observation time slot;

demodulating said response signal to separate the response to each excitation signal by quantizing and processing of said response signal at each sampling instances of said set of sampling instances, wherein said quantizing is performed for each frequency of interest of said response signal in predetermined mode; and

obtaining a set of main frequency and time tagged estimates, each of said estimates comprising a plurality of sub-estimates, said sub-estimates representing a signal value, for said observation time slot, wherein the length of said observation time slot is selected so that it contains integer numbers of periods of each excitation signal of said set of excitation signals, wherein each of said estimates are sorted separately so the estimates at each frequency of interest are obtained while impacts of other estimates are cancelled out.

2. A method as in claim 1, wherein said response signal is sampled with predetermined sampling pattern within said observation time slot.

3. A method as in claim 2, wherein said response signal is sampled uniformly within said observation time slot.

4. A method as in claim 2, wherein undersampling is applied to demodulate at least one of said components of said response signal.

5. A method as in claim 2, wherein separating measurement channels in time domain is performed by uniformly sampling said response signal for one frequency of interest for every measurement channel in one observation time slot and for each additional frequency of interest for every measurement channel in corresponding additional observation time slot.

6. A method as in claim 2, wherein separating measurement channels is performed by non-uniformly sampling said response signal for each frequency of interest for one measurement channel in first observation time slot and for each frequency of interest for each additional measurement channel in corresponding additional observation time slot.

7. A method as in claim 6, wherein said response signal is quantized more than once at predetermined sampling rate for at least one sampling instance of said set of sampling instances.

8. A method as in claim 7, wherein said sampling instances are calculated, using said preselected frequencies of each excitation signal of said set of excitation signals and for each frequency of interest.

9. A method as in claim 8, wherein said sampling instances are estimated for said set of excitation signals and for each frequency of interest.

10. A method as in claim 1 wherein at least one frequency of interest of said response signal is compensated by a predetermined signal.

11. A method as in claim 1 wherein at least one frequency of interest of said response signal is compensated by an estimated signal.

12. A method as in claim 1 wherein at least one frequency of interest of said response signal is compensated using analogue compensation.

13. A method as in claim 1 wherein at least one frequency of interest of said response signal is compensated using digital compensation.

14. A method as in claim 1, wherein an inphase sampling is applied for said frequency of interest, wherein obtained inphase samples are used to separate inphase estimates.

15. A method as in claim 1, wherein a quadrature sampling is applied for said frequency of interest, wherein obtained quadrature samples are used to separate quadrature estimates.

16. A method as in claim 1, wherein off the phase sampling is applied for each frequency of interest, whereas obtained off the phase samples are used for extracting an additive signal component removal of an offset, or a common-mode noise by summing cancellation technique.

17. A method as in claim 1, wherein said preselected frequencies f_p of said excitation signals of said set of excitation signals are chosen according to formula

f_p=J_p*2^p−

1*f_R,

where J_p is an integer, p is an integer from 1 to P, where P is the number of excitation signals in said set of excitation signal, and f_R is a repetition rate of an excitation and measurement cycle.