RU2669484C1 - Method of determining components of impedance of bio-object - Google Patents
Method of determining components of impedance of bio-object Download PDFInfo
- Publication number
- RU2669484C1 RU2669484C1 RU2017122573A RU2017122573A RU2669484C1 RU 2669484 C1 RU2669484 C1 RU 2669484C1 RU 2017122573 A RU2017122573 A RU 2017122573A RU 2017122573 A RU2017122573 A RU 2017122573A RU 2669484 C1 RU2669484 C1 RU 2669484C1
- Authority
- RU
- Russia
- Prior art keywords
- characteristic
- impedance
- values
- calibration
- measured
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 239000003814 drug Substances 0.000 abstract description 3
- 238000010606 normalization Methods 0.000 abstract description 3
- 239000000126 substance Substances 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 2
- 244000309464 bull Species 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
Description
Предлагаемое изобретение относится к медицине и может быть использовано для оценки функционального состояния организма.The present invention relates to medicine and can be used to assess the functional state of the body.
Известен способ измерения электрических величин активного сопротивления, емкости и индуктивности [А.с. 1797079 СССР, МКИЗ G01R 27/18], согласно которому на последовательную активно-емкостную или активно индуктивную цепь подают напряжение постоянного тока. При этом один элемент цепи известен. После подачи напряжения через определенные промежутки времени Δt измеряют первое и второе мгновенные значения напряжения на средней точке измерительной цепи. Неизвестные элементы определяют соответственно по формулам для активно-емкостной и индуктивно-емкостной цепей.A known method of measuring electrical quantities of active resistance, capacitance and inductance [A.S. 1797079 USSR, MKIZ G01R 27/18], according to which a DC voltage is applied to a series active-capacitive or active inductive circuit. In this case, one element of the chain is known. After applying voltage at certain time intervals Δt, the first and second instantaneous voltage values are measured at the midpoint of the measuring circuit. Unknown elements are determined respectively by the formulas for active-capacitive and inductive-capacitive circuits.
Недостаток такого способа измерений сопротивлений в том, что он не позволяет измерить с достаточной точностью значения активной и реактивной составляющих комплексного сопротивления.The disadvantage of this method of measuring resistance is that it does not allow measuring with sufficient accuracy the values of the active and reactive components of the complex resistance.
По способу определения составляющих импеданса биообъекта [см. А.с. СССР №1397024, МПК А61В 5/05, публ. 1988 г., Бил. №19] на биообъект накладывают электроды, через которые подается импульс тока определенной полярности и с амплитудой I0. Так как составляющая импеданса имеет емкостной характер, происходит переходной процесс нарастания напряжения. В моменты времени t1 и t2 измеряют напряжения U1 и U2. Измерение в момент t2 происходит тогда, когда емкость тканей заряжена полностью, т.е. переходной процесс закончился. Величина I0 выбирается такой, чтобы за время действия импульса тока произошел полный заряд емкости тканей. Тогда напряжение на биообъекте пропорционально величине активной составляющей импеданса биообъекта.By the method of determining the components of the impedance of a biological object [see A.S. USSR No. 1397024, IPC А61В 5/05, publ. 1988, Beale. No. 19] electrodes are applied to the biological object through which a current pulse of a certain polarity and with an amplitude of I 0 is supplied. Since the component of the impedance is capacitive in nature, a transient process of increasing voltage occurs. At times t 1 and t 2 measure the voltage U 1 and U 2 . Measurement at time t 2 occurs when the tissue capacity is fully charged, i.e. The transition process is over. The value of I 0 is chosen so that during the duration of the current pulse there is a full charge of the tissue capacity. Then the voltage on the biological object is proportional to the value of the active component of the impedance of the biological object.
Недостатками являются: наличие динамической и методической погрешности и низкая оперативность, вызванная необходимостью ожидания установившегося режима.The disadvantages are: the presence of dynamic and methodological errors and low efficiency, due to the need to wait for the steady state.
Наиболее близким к заявленному техническому решению является способ определения составляющих импеданса биообъекта [Пат. РФ №2509531, МПК8 А61В 5/05, публ. 2014, Бюл. №5.], заключающийся в том, что на биообъект подается импульс тока I0 и измеряют напряжения в моменты времени t2=2t1. По измеренным значениям напряжения и моментам времени регистрируют информативные параметры: потенциал Е и постоянную времени Т, по которым определяют значение активного сопротивления и эквивалентную емкость тканей биообъекта.Closest to the claimed technical solution is a method for determining the components of the impedance of a biological object [Pat. RF №2509531, IPC8 А61В 5/05, publ. 2014, bull. No. 5.], Which consists in the fact that a current pulse I 0 is supplied to the biological object and voltages are measured at time t 2 = 2t 1 . Informative parameters are recorded from the measured voltage values and time instants: potential E and time constant T, which determine the value of active resistance and the equivalent tissue capacity of the biological object.
Недостатком прототипа является то, что он рассчитан на случай, когда оба информативных параметра известны, но, как правило, на практике один из информативных параметров неизвестен.The disadvantage of the prototype is that it is designed for the case when both informative parameters are known, but, as a rule, in practice one of the informative parameters is unknown.
Технической задачей является определение составляющих импеданса биообъекта при неизвестном информативном параметре - максимальном значении потенциала Е.The technical task is to determine the components of the impedance of a biological object with an unknown informative parameter - the maximum value of potential E.
Данная техническая задача достигается тем, что: в способе определения составляющих импеданса биологического объекта, заключающемся в подаче на биообъект импульса стабилизированного тока I0 и измерении напряжения u в момент времени t после начала импульса тока, в качестве составляющих импеданса биообъекта определяют активное сопротивление R и эквивалентную емкость C тканей биообъекта, в отличие от прототипа максимальное значение потенциала Е определяют по калибровочной характеристике, калибровку проводят априори для измеренного Ui и известного Uэi значений напряжения (i=1, 2) в два момента времени t2=2t1, калибровочной характеристикой служит функция максимального значения потенциала Е0i, компенсирующая неопределенность постоянной времени T выбранной произвольно T*, и связывающая эталонную Uэi и измеренную Ui характеристики определения импеданса за счет нормирования измеренных значений известными, по калибровочной характеристике Е0i находят действительные значения постоянной времени T и максимальную величину потенциала E, по которым последовательно строят калибровочную характеристику Е0i и действительную U∂ характеристику определения импеданса.This technical problem is achieved by the fact that: in the method for determining the components of the impedance of a biological object, which consists in applying a stabilized current pulse I 0 to the biological object and measuring the voltage u at time t after the start of the current pulse, the active resistance R and equivalent the capacity C of the tissues of the biological object, in contrast to the prototype, the maximum value of the potential E is determined by the calibration characteristic, calibration is carried out a priori for the measured U i and from of the known U ei voltage values (i = 1, 2) at two time instants t 2 = 2t 1 , the calibration characteristic is the function of the maximum potential value E 0i , which compensates for the uncertainty of the time constant T chosen arbitrarily T *, and relating the reference U ei and the measured U i determining the impedance characteristics due to normalization of the measured values known, the calibration characteristic E 0i are the actual values of the time constant T and the maximum value of the potential E, at which a calibration series E 0i characteristic and actual characteristic U ∂ determining impedance.
Сущность предлагаемого способа поясняется на фиг. 1÷2. Предлагаемый способ включает следующие этапы:The essence of the proposed method is illustrated in FIG. 1 ÷ 2. The proposed method includes the following steps:
1. Определяют максимальное значение потенциала E по калибровочной функции E0i=E0i(t).1. Determine the maximum value of the potential E by the calibration function E 0i = E 0i (t).
2. Калибровку проводят априори для известных эталонной Uэi ((фиг. 1, 1) и Ui (фиг. 1, 2) измеренной значений напряжения.2. Calibration is carried out a priori for the known reference U ei ((Fig. 1, 1) and U i (Fig. 1, 2) measured voltage values.
3. Калибровочной характеристикой служит функция E0i (фиг. 2, 3) максимальной величины потенциала E (фиг. 1, 3), компенсирующая неопределенность постоянной времени T, (фиг. 1, 4) выбранной произвольно T* (фиг. 1, 4а), и связывающая эталонную Uэ и измеренную U зависимости за счет нормирования измеренных значений известными:3. The calibration characteristic is the function E 0i (Fig. 2, 3) of the maximum value of the potential E (Fig. 1, 3), compensating for the uncertainty of the time constant T, (Fig. 1, 4) of a randomly selected T * (Fig. 1, 4a ), and linking the reference U e and the measured U dependences due to the normalization of the measured values known:
По калибровочной характеристике Е0i (фиг. 2, 3) восстанавливают действительную характеристику U∂ (фиг. 2, 4):From the calibration characteristic E 0i (Fig. 2, 3) restore the actual characteristic U ∂ (Fig. 2, 4):
которая максимально приближена к эталонной Uэi (фиг. 2, 1):which is as close as possible to the reference U ei (Fig. 2, 1):
Эталонная характеристика Uэi (фиг. 2, 1) и характеристика, ей тождественная, U∂ (фиг. 2, 4) получены из экспоненциальной динамической характеристики с искомыми информативными параметрами T, Е:The reference characteristic U ei (Fig. 2, 1) and the characteristic identical to it, U ∂ (Fig. 2, 4) are obtained from the exponential dynamic characteristic with the desired informative parameters T, E:
где Т - постоянная времени (фиг. 1, 4) процесса и Е - максимальная величина потенциала (фиг. 1, 3). Физический смысл информативных параметров следует из предельных соотношений:where T is the time constant (Fig. 1, 4) of the process and E is the maximum value of the potential (Fig. 1, 3). The physical meaning of informative parameters follows from the limit relations:
т.е. Е - максимальная величина потенциала для t=∞.those. E is the maximum value of the potential for t = ∞.
т.е. Т - постоянная времени, т.к.those. T is the time constant, because
На практике один из информативных параметров исследуемой характеристики, как правило, неизвестен. В этом случае один параметр выбираем произвольно T* (фиг. 1, 4а), а второй принимает вид функции E0i (фиг. 2, 3), которая компенсирует незнание первого информативного параметра T (фиг. 1, 4). По калибровочной функции E0i нормируется измеренная кривая Ui до тождественного эквивалента Uэi=U∂ (фиг. 2, 1 и 4).In practice, one of the informative parameters of the studied characteristics is usually unknown. In this case, one parameter is chosen arbitrarily T * (Fig. 1, 4a), and the second takes the form of a function E 0i (Fig. 2, 3), which compensates for the ignorance of the first informative parameter T (Fig. 1, 4). According to the calibration function E 0i, the measured curve Ui is normalized to the identical equivalent U ei = U ∂ (Fig. 2, 1 and 4).
Задаем произвольно параметр T*=const (фиг. 1, 4а) вместо неизвестного действительного значения постоянной времени T (фиг. 1, 4). Для компенсации произвольности константы T* (фиг. 1, 4а) максимальное значение потенциала E (фиг. 1, 3) превратится в характеристику E0i (фиг. 2, 3), компенсирующую незнание постоянной времени T (фиг. 1, 4), где i - число измерений (i=1, 2).We arbitrarily set the parameter T * = const (Figs. 1, 4a) instead of the unknown real value of the time constant T (Figs. 1, 4). To compensate for the arbitrariness of the constant T * (Fig. 1, 4a), the maximum value of the potential E (Fig. 1, 3) will turn into a characteristic E 0i (Fig. 2, 3), compensating for the ignorance of the time constant T (Fig. 1, 4), where i is the number of measurements (i = 1, 2).
Калибровочной функцией для неизвестных параметров T, Е служит динамическая характеристика E0i (фиг. 2, 3).The calibration function for the unknown parameters T, E is the dynamic characteristic E 0i (Fig. 2, 3).
Калибровочную характеристику E0i выразим из системы уравнений с известными параметрами T, Е характеристики Uэi, являющейся эталонной (получено путем аппроксимации экспериментальных данных), и характеристики Ui, (фиг. 1, 2), измеренной с произвольной константой T* (фиг. 1, 4а) и характеристикой E0i:We express the calibration characteristic E 0i from the system of equations with known parameters T, E of the characteristic U ei , which is the reference (obtained by approximating the experimental data), and the characteristic U i (Fig. 1, 2), measured with an arbitrary constant T * (Fig. 1, 4a) and characteristic E 0i :
В соответствии с закономерностями калибровки Uэi=Ui, следует калибровочная характеристика E0i (фиг. 2, 3), связывающая между собой эталонную Uэi (фиг. 1, 1) и измеренную Ui, (фиг. 1, 2) характеристики определения импеданса:In accordance with the laws of calibration U ei = U i , follows the calibration characteristic E 0i (Fig. 2, 3), connecting the reference U ei (Fig. 1, 1) and the measured U i (Fig. 1, 2) characteristics impedance definitions:
Следовательно, калибровочной характеристикой служит функция максимального значения потенциала Е0i (фиг. 2, 3), компенсирующая неопределенность значения постоянной времени Т (фиг. 1, 4), выбранной произвольно Т* (фиг. 1, 4а).Therefore, the calibration characteristic is a function of the maximum potential value E 0i (Fig. 2, 3), compensating for the uncertainty of the value of the time constant T (Fig. 1, 4), chosen arbitrarily T * (Fig. 1, 4a).
4. По калибровочной характеристике E0i (фиг. 2, 3) находят действительные значения потенциала E (фиг. 1, 3) и постоянной времени T (фиг. 1, 4), которые являются информативными параметрами, доставляющими оптимум калибровочной характеристике. Из уравнения (5) составим систему уравнений для i=1, 2:4. From the calibration characteristic E 0i (Fig. 2, 3), the actual values of the potential E (Fig. 1, 3) and the time constant T (Fig. 1, 4) are found, which are informative parameters that deliver the optimum to the calibration characteristic. From equation (5) we compose a system of equations for i = 1, 2:
Поделив одно уравнение системы (6) на другое и проэкспоненцировав, при условии t2=2t1, определяют алгоритм оптимизации постоянной времени Т:Dividing one equation of system (6) into another and exponentially, provided t 2 = 2t 1 , determine the optimization algorithm for the time constant T:
Выразим Е из первого уравнения системы (6), подставив найденное Т:Express E from the first equation of system (6), substituting the found T:
5. По полученным информативным параметрам (7) и (8) строят характеристику U∂ (5) (фиг. 2, 4), по которой находят действительную характеристику определения составляющих импеданса, тождественную эквиваленту (4) (фиг. 1, 1).5. Using the obtained informative parameters (7) and (8), we construct the characteristic U ∂ (5) (Fig. 2, 4), by which the actual characteristic for determining the components of the impedance is found, which is identical to the equivalent (4) (Fig. 1, 1).
Адекватность предлагаемого способа физике эксперимента доказывает математическое моделирование действительной характеристики U∂ (фиг. 2, 4), относительно эквивалента экспериментальной характеристики Uэ (фиг. 2, 1), по полученным значениям.The adequacy of the proposed method to the physics of the experiment is proved by mathematical modeling of the actual characteristic U ∂ (Fig. 2, 4), relative to the equivalent experimental characteristic U e (Fig. 2, 1), according to the obtained values.
Проводят оценку адекватности полученных зависимостей по формуле определения относительной погрешности:The adequacy of the obtained dependencies is assessed by the formula for determining the relative error:
ее оценка представлена на фиг. 3.its assessment is presented in FIG. 3.
Относительная погрешность моделирования не превышает 1.2⋅10-15.The relative modeling error does not exceed 1.2⋅10 -15 .
Динамическая погрешность δизм (фиг. 4, 1) измеренной характеристики Ui увеличивается с течением времени с 0 до 75%:The dynamic error δ ISM (Fig. 4, 1) of the measured characteristic U i increases over time from 0 to 75%:
Динамическая погрешность δд (фиг. 4, 2) действительной характеристики U∂ не превышает 0,01.The dynamic error δ d (Fig. 4, 2) of the actual characteristic U ∂ does not exceed 0.01.
Отклонения измеренных характеристик прототипа приводят к большой динамической погрешности относительно нулевой погрешности теоретических значений, что доказывает эффективность предлагаемого способа.Deviations of the measured characteristics of the prototype lead to a large dynamic error relative to the zero error of theoretical values, which proves the effectiveness of the proposed method.
Таким образом, определение составляющих импеданса биообъекта по калибровочной характеристике, компенсирующей неопределенность постоянной времени, выбранной произвольно, по которой определяют действительные значения информативных параметров, действительную характеристику определения составляющих импеданса, в отличие от известных решений, повышает эффективность в несколько раз.Thus, the determination of the components of the biological object impedance by a calibration characteristic that compensates for the uncertainty of the time constant, chosen arbitrarily, by which the actual values of the informative parameters are determined, the actual characteristic of the determination of the impedance components, in contrast to the known solutions, increases the efficiency several times.
Claims (7)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2017122573A RU2669484C1 (en) | 2017-06-27 | 2017-06-27 | Method of determining components of impedance of bio-object |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2017122573A RU2669484C1 (en) | 2017-06-27 | 2017-06-27 | Method of determining components of impedance of bio-object |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2669484C1 true RU2669484C1 (en) | 2018-10-11 |
Family
ID=63862308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2017122573A RU2669484C1 (en) | 2017-06-27 | 2017-06-27 | Method of determining components of impedance of bio-object |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2669484C1 (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6339722B1 (en) * | 1995-09-26 | 2002-01-15 | A. J. Van Liebergen Holding B.V. | Apparatus for the in-vivo non-invasive measurement of a biological parameter concerning a bodily fluid of a person or animal |
UA65068C2 (en) * | 2003-05-26 | 2007-03-15 | Viktor Oleksiiovych Yaruta | Method for measuring electrical parameters of living tissues |
US20090248118A1 (en) * | 2001-12-04 | 2009-10-01 | Boston Scientific Neuromodulation Corporation | Apparatus and method for determining the relative position and orientation of neurostimulation leads |
RU2509531C1 (en) * | 2012-07-04 | 2014-03-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО ТГТУ | Method for determining bioobject impedance components |
RU2586457C1 (en) * | 2015-02-25 | 2016-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО ТГТУ | Method of determining components of impedance of biological effect |
-
2017
- 2017-06-27 RU RU2017122573A patent/RU2669484C1/en not_active IP Right Cessation
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6339722B1 (en) * | 1995-09-26 | 2002-01-15 | A. J. Van Liebergen Holding B.V. | Apparatus for the in-vivo non-invasive measurement of a biological parameter concerning a bodily fluid of a person or animal |
US20090248118A1 (en) * | 2001-12-04 | 2009-10-01 | Boston Scientific Neuromodulation Corporation | Apparatus and method for determining the relative position and orientation of neurostimulation leads |
UA65068C2 (en) * | 2003-05-26 | 2007-03-15 | Viktor Oleksiiovych Yaruta | Method for measuring electrical parameters of living tissues |
RU2509531C1 (en) * | 2012-07-04 | 2014-03-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО ТГТУ | Method for determining bioobject impedance components |
RU2586457C1 (en) * | 2015-02-25 | 2016-06-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ФГБОУ ВПО ТГТУ | Method of determining components of impedance of biological effect |
Non-Patent Citations (1)
Title |
---|
НИКОЛАЕВ Д.В. и др. Биоимпедансный анализ состава тела человека, М., Наука, 2009, с.21-25, 44-54. * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101144684B1 (en) | Apparatus for evaluating characteristics of cell | |
KR102037954B1 (en) | System, controller, and method for determining conductance of an object | |
US10677828B2 (en) | Low frequency impedance measurement with source measure units | |
JP6410986B1 (en) | Battery impedance evaluation apparatus and battery impedance evaluation method | |
Czaja | A microcontroller system for measurement of three independent components in impedance sensors using a single square pulse | |
Mazloomi et al. | Analysis of the frequency response of a water electrolysis cell | |
CN109425833A (en) | Temperature estimating device | |
RU2669484C1 (en) | Method of determining components of impedance of bio-object | |
US11705894B2 (en) | Pulsed high current technique for characterization of device under test | |
JP5717671B2 (en) | Model constant acquisition method and model constant acquisition apparatus | |
RU2509531C1 (en) | Method for determining bioobject impedance components | |
RU152732U1 (en) | ELECTROMETRIC ION METER | |
RU2548780C1 (en) | Method for assessing functional state of haemostasis system | |
RU160681U1 (en) | PULSE RELAXATION DEVICE FOR THE EVALUATION OF NICKEL-CADMIUM BATTERIES | |
CN111989581B (en) | Battery monitoring system | |
RU2586457C1 (en) | Method of determining components of impedance of biological effect | |
KR102037159B1 (en) | A contactless impedance readout system using simulated inductor and switched capacitor | |
US20200400751A1 (en) | Estimation of unknown electronic load | |
RU2594376C1 (en) | Method of capacitors self-discharge time constant measuring | |
RU2672533C1 (en) | Device for measurement of frequency error of thermoelectric converters | |
Kowalewski et al. | Fast high-impedance spectroscopy method using sinc signal excitation | |
JP2017198547A (en) | Insulation degradation testing device and method for testing insulation deterioration | |
RU2616871C1 (en) | Method of determining current localization voltage in powerful hf and uhf bipolar transistors | |
RU2167429C1 (en) | Method measuring thermal resistance of two-terminal networks with well-known temperature coefficient of resistance | |
US6842014B2 (en) | Methods for determining inductance and resistance of an inductor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20200628 |