CN115153580B - A skin impedance measurement method and system for fetal electrocardiogram monitoring - Google Patents

A skin impedance measurement method and system for fetal electrocardiogram monitoring Download PDF

Info

Publication number
CN115153580B
CN115153580B CN202210909613.1A CN202210909613A CN115153580B CN 115153580 B CN115153580 B CN 115153580B CN 202210909613 A CN202210909613 A CN 202210909613A CN 115153580 B CN115153580 B CN 115153580B
Authority
CN
China
Prior art keywords
electrode
fetal
impedance measurement
human body
skin impedance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210909613.1A
Other languages
Chinese (zh)
Other versions
CN115153580A (en
Inventor
陆泽宏
李玮琛
赵妍
张博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GBA National Institute for Nanotechnology Innovation
Original Assignee
GBA National Institute for Nanotechnology Innovation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GBA National Institute for Nanotechnology Innovation filed Critical GBA National Institute for Nanotechnology Innovation
Priority to CN202210909613.1A priority Critical patent/CN115153580B/en
Publication of CN115153580A publication Critical patent/CN115153580A/en
Application granted granted Critical
Publication of CN115153580B publication Critical patent/CN115153580B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/344Foetal cardiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0531Measuring skin impedance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A61B5/1118Determining activity level
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analogue processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00Evaluating a particular growth phase or type of persons or animals
    • A61B2503/02Foetus

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Dermatology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Mathematical Physics (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种用于胎儿心电监测的皮肤阻抗测量方法和系统,利用三轴陀螺仪采集孕妇的人体姿态信息,在孕妇处于平静状态下才启动阻抗测量环节,启动交流激励电流源,通过P电极对孕妇皮肤表面施加预置交流电流,通过差分放大器获取P电极和N电极之间的差分信号作为激励响应信号,对激励响应信号进行预处理,得到滤除母胎生物电信号后的激励响应信号,再对滤除母胎生物电信号后的激励响应信号进行离散傅里叶变换,求解出基波频率实部和虚部的复数结果,得到人体皮肤阻抗,解决了现有的皮肤与电极界面处的阻抗测量方法存在的忽略激励响应信号的虚部信息以及采用模拟开关器件带来噪声干扰,导致阻抗测量准确度下降的技术问题。

The present invention discloses a skin impedance measurement method and system for fetal electrocardiogram monitoring. A three-axis gyroscope is used to collect human posture information of a pregnant woman. The impedance measurement link is started only when the pregnant woman is in a calm state. An AC excitation current source is started, and a preset AC current is applied to the skin surface of the pregnant woman through a P electrode. A differential amplifier is used to obtain a differential signal between the P electrode and the N electrode as an excitation response signal. The excitation response signal is preprocessed to obtain an excitation response signal after filtering out the maternal and fetal bioelectric signals. The excitation response signal after filtering out the maternal and fetal bioelectric signals is then subjected to a discrete Fourier transform to solve the complex results of the real part and the imaginary part of the fundamental frequency to obtain the human skin impedance. The technical problems of ignoring the imaginary part information of the excitation response signal and using analog switch devices to cause noise interference, resulting in reduced impedance measurement accuracy, are solved in the existing impedance measurement method at the interface between the skin and the electrode.

Description

Skin impedance measurement method and system for fetal electrocardiograph monitoring
Technical Field
The invention relates to the technical field of electrocardiograph monitoring, in particular to a skin impedance measurement method and system for fetal electrocardiograph monitoring.
Background
The non-invasive fetal electrocardiosignal (NI-FECG) equipment can collect maternal and fetal electrocardiosignals through an electrode arranged on the abdomen of a pregnant woman, separate the maternal and fetal electrocardiosignals to obtain the Maternal Electrocardiosignals (MECG) and the Fetal Electrocardiosignals (FECG), and calculate the fetal heart rate and uterine contraction according to the maternal and fetal electrocardiosignals. However, the FECG signal is very weak, usually 10-50 microvolts, and is easy to submerge in interference signals of human bodies, environments and circuits. The low and stable impedance at the interface of the skin and the electrode is the key of the acquired high-fidelity biological signals, and can provide a quality evaluation basis for the maternal-fetal electrocardiosignal and correct the result of separating the maternal-fetal signals. However, on the one hand, the complexity of human skin causes great signal quality change, and on the other hand, when the electrode is attached to the skin, the problems of wrinkling, falling off and the like are easily caused, so that a method for accurately measuring the interface impedance of the skin and the electrode is needed, when a good fetal signal cannot be acquired due to overlarge interface impedance, impedance information can be fed back to a user, and the user is prompted to perform inspection so as to ensure the accuracy of acquiring fetal electrocardiosignals.
The existing impedance measurement method at the interface of the skin and the electrode generally adopts a method of carrying out ohm law conversion on the amplitude of the excitation response signal, solves the impedance information, ignores the imaginary part information of the excitation response signal, cannot effectively combine the capacitive impedance information to calculate the impedance value, easily causes the measurement error to be amplified, is not matched with the impedance model of the human body, and leads to the reduction of the measurement accuracy of the human body impedance. Meanwhile, in the prior art, an analog switching device is adopted in the impedance measurement method at the interface of the skin and the electrode, and signals of the electrode pair are routed to an impedance signal measurement module or an electrocardiosignal measurement module, so that the impedance signal and the electrocardiosignal are distinguished on a physical link, but the analog switching device has thermal noise and is easy to introduce external noise, and is easy to interfere with weak bioelectricity signals, so that the impedance measurement precision is reduced. Therefore, the invention provides a skin impedance measurement method and a skin impedance measurement system for fetal electrocardio monitoring, which are used for solving the technical problems that the accuracy of impedance measurement is reduced due to the fact that the existing impedance measurement method at the interface of skin and an electrode ignores the imaginary part information of an excitation response signal and adopts an analog switching device to bring noise interference.
Disclosure of Invention
The invention provides a skin impedance measurement method and a skin impedance measurement system for fetal electrocardio monitoring, which are used for solving the technical problems that the accuracy of impedance measurement is reduced due to the fact that the existing impedance measurement method at the interface of skin and an electrode ignores the imaginary part information of an excitation response signal and adopts an analog switching device to bring noise interference.
In view of this, a first aspect of the present invention provides a skin impedance measurement method for fetal electrocardiography monitoring, comprising:
Acquiring human body posture information through a triaxial gyroscope arranged on a pregnant woman;
judging whether the pregnant woman is in a calm state or not according to the human body posture information, if so, starting an impedance measurement link, and measuring the skin impedance of the human body;
The impedance measurement link includes:
Starting an alternating current excitation current source, applying preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, and obtaining a differential signal between the P electrode and an N electrode through a differential amplifier to serve as an excitation response signal, wherein the alternating current excitation current source, the P electrode and the N electrode are connected in a ring shape to form a loop, and the P electrode and the N electrode are connected with two input ends of the differential amplifier;
preprocessing the excitation response signal to obtain an excitation response signal after the maternal biological electrical signal is filtered, wherein the preprocessing comprises analog-to-digital conversion processing and filtering processing;
performing discrete Fourier transform on the excitation response signal after the maternal biological electrical signal is filtered, and solving a complex result of a real part and an imaginary part of fundamental frequency to obtain the skin impedance of the human body.
Optionally, the method further comprises:
after the measurement of the skin impedance of the human body is completed, an alternating current excitation current source is turned off, a fetal electrocardiosignal measuring link is started, the maternal-fetal electrocardiosignal is collected through a P electrode, an N electrode and other collecting electrodes, and the maternal-fetal electrocardiosignal is subjected to maternal-fetal electrocardiosignal separation by combining the measured skin impedance of the human body, so that the fetal electrocardiosignal is obtained.
Optionally, acquiring the human posture information by a triaxial gyroscope mounted on the pregnant woman includes:
Acquiring an instantaneous angular velocity value of a human body in a X, Y, Z dimension direction acquired by a triaxial gyroscope arranged on a pregnant woman in real time, and calculating a vector module length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment;
Caching the vector modular length at each moment to obtain a vector modular length time sequence;
carrying out local spectrum transformation of unit window data on the vector modulus long-time sequence by adopting short-time Fourier transformation, and converting the local spectrum data into power spectrum data;
And summing the power spectrum data, calculating a power total value, judging that the current human body posture is in an active state if the power total value exceeds a threshold value, and judging that the current human body posture is in a static state if the power total value exceeds the threshold value.
Optionally, calculating the vector modular length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment includes:
The instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment is represented as a vector space, the Euclidean norm of the vector space at each moment is calculated, and the vector modular length at each moment is obtained.
Alternatively, a square wave excitation current of 24nA at 250Hz is preset as the alternating current.
Optionally, performing discrete fourier transform on the excitation response signal after filtering the maternal-fetal bioelectric signal, and solving a complex result of a real part and an imaginary part of fundamental frequency to obtain skin impedance of the human body, including:
constructing an excitation response signal with the filtered maternal-fetal bioelectric signals into a finite-length complex sequence with the imaginary number of 0, and performing discrete Fourier transform on the finite-length complex sequence with the imaginary number of 0;
The complex information with the fundamental wave frequency component of 250Hz is taken out from the discrete Fourier transform result, and the human skin impedance is calculated by combining the differential amplification factor of the differential amplifier, wherein the calculation formula is as follows:
Wherein Z is the human skin impedance, b is the differential amplification factor of the differential amplifier, R 250 is the real value of the complex information with the fundamental frequency component of 250Hz extracted from the discrete Fourier transform result, and I 250 is the imaginary value of the complex information with the fundamental frequency component of 250Hz extracted from the discrete Fourier transform result.
The invention provides a skin impedance measurement system for fetal electrocardio monitoring, which comprises a gyroscope data analysis module, a P electrode, an N electrode, a three-axis gyroscope, an alternating current excitation current source, a differential amplifier, an analog-to-digital converter, a right leg driving circuit module, an information interaction interface and an impedance measurement module, wherein the gyroscope data analysis module is used for analyzing the data of the fetus;
the three-axis gyroscope is connected with the gyroscope data analysis module;
the gyroscope data analysis module is connected with an alternating current excitation current source;
The positive end of the alternating current excitation current source is connected with the P electrode, the negative end of the alternating current excitation current source is connected with the N electrode, the P electrode and the N electrode are connected with two input ends of the differential amplifier, the output end of the differential amplifier is connected with the input end of the analog-to-digital converter and the right leg driving circuit module, the right leg driving circuit module comprises a right leg driving circuit and an E electrode, and the triaxial gyroscope, the P electrode and the N electrode are attached to the abdomen of a pregnant woman;
the output end of the analog-to-digital converter is connected with the impedance measurement module, and the output end of the impedance measurement module is connected with the information interaction interface;
When the skin impedance measuring system for fetal electrocardio monitoring is used for measuring the skin impedance of a human body, the three-axis gyroscope collects the human body posture information and sends the human body posture information to the gyroscope data analysis module, the gyroscope data analysis module judges whether the pregnant woman is in a calm state according to the human body posture information, if so, an impedance measuring link is started, and the skin impedance of the human body is measured;
The impedance measurement link includes:
Starting an alternating current excitation current source, applying preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, and obtaining a differential signal between the P electrode and an N electrode through a differential amplifier to serve as an excitation response signal, wherein the alternating current excitation current source, the P electrode and the N electrode are connected in a ring shape to form a loop, and the P electrode and the N electrode are connected with two input ends of the differential amplifier;
Preprocessing the excitation response signal through an impedance measurement module to obtain an excitation response signal after the maternal biological electrical signal is filtered, wherein the preprocessing comprises analog-to-digital conversion processing and filtering processing;
And performing discrete Fourier transform on the excitation response signal after the maternal biological electrical signal is filtered by an impedance measurement module, and solving complex results of a real part and an imaginary part of fundamental frequency to obtain the skin impedance of the human body.
Optionally, the device also comprises a fetal electrocardiograph measurement module and other acquisition electrodes;
the other acquisition electrodes are connected with the input end of the differential amplifier, the input end of the fetal electrocardiograph measurement module is connected with the output end of the analog-to-digital converter, and is also connected with the impedance measurement module, and the output end of the fetal electrocardiograph measurement module is connected with the information interaction interface;
The fetal electrocardiosignal module is used for measuring fetal electrocardiosignals after finishing the impedance measurement of the skin of a human body;
After the measurement of the skin impedance of the human body is completed, the impedance measurement module sends a control signal to the alternating current excitation current source, so that the alternating current excitation current source is turned off, the maternal-fetal electrocardiosignals are collected through the P electrode, the N electrode and other collecting electrodes, and the maternal-fetal electrocardiosignals are subjected to maternal-fetal electrocardiosignal separation by combining the skin impedance of the human body measured by the impedance measurement module, so that the fetal electrocardiosignals are obtained.
Optionally, the gyroscope data analysis module is specifically configured to:
Acquiring the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction acquired by the triaxial gyroscope in real time, and calculating the vector module length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment;
Caching the vector modular length at each moment to obtain a vector modular length time sequence;
carrying out local spectrum transformation of unit window data on the vector modulus long-time sequence by adopting short-time Fourier transformation, and converting the local spectrum data into power spectrum data;
And summing the power spectrum data, calculating a power total value, judging that the current human body posture is in an active state if the power total value exceeds a threshold value, and judging that the current human body posture is in a static state if the power total value exceeds the threshold value.
Optionally, calculating the vector modular length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment includes:
The instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment is represented as a vector space, the Euclidean norm of the vector space at each moment is calculated, and the vector modular length at each moment is obtained.
From the above technical scheme, the skin impedance measurement method and system for fetal electrocardiograph monitoring provided by the invention have the following advantages:
The skin impedance measurement method for fetal electrocardiosignal monitoring provided by the invention comprises the steps of collecting human body posture information of a pregnant woman by using a triaxial gyroscope, judging whether the pregnant woman is in a calm state currently, starting an impedance measurement link when the pregnant woman is in the calm state, starting an alternating current excitation current source, applying preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, obtaining a differential signal between the P electrode and an N electrode through a differential amplifier as an excitation response signal, preprocessing the excitation response signal to obtain the excitation response signal after filtering the maternal-fetal bioelectric signal, and performing discrete Fourier transform on the excitation response signal after filtering the maternal-fetal bioelectric signal to solve complex results of a fundamental frequency real part and an imaginary part, thereby obtaining the skin impedance of the human body.
The principle and effect of the skin impedance measurement system for fetal electrocardiograph monitoring provided by the invention are the same as those of the skin impedance measurement method for fetal electrocardiograph monitoring provided by the invention, and the skin impedance measurement system for fetal electrocardiograph monitoring is not repeated here.
Drawings
For a clearer description of embodiments of the invention or of solutions according to the prior art, the figures which are used in the description of the embodiments or of the prior art will be briefly described, it being obvious that the figures in the description below are only some embodiments of the invention, from which, without the aid of inventive efforts, other relevant figures can be obtained for a person skilled in the art.
FIG. 1 is a flow chart of a skin impedance measurement method for fetal electrocardiograph monitoring provided in the present invention;
FIG. 2 is a schematic diagram of a skin impedance measurement system for fetal electrocardiograph monitoring according to the present invention;
Fig. 3 is a schematic diagram of a fetal electrocardiograph flow path of a skin impedance measurement system for fetal electrocardiograph monitoring provided in the present invention.
Detailed Description
In order to make the present invention better understood by those skilled in the art, the following description will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings, and it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
For ease of understanding, referring to fig. 1, an embodiment of a skin impedance measurement method for fetal electrocardiographic monitoring is provided in the present invention, comprising:
step 101, acquiring human body posture information through a triaxial gyroscope arranged on a pregnant woman.
It should be noted that, in the embodiment of the present invention, the three-axis gyroscope on the pregnant woman may be mounted on the abdomen of the pregnant woman in a pasting manner, and the three-axis gyroscope adopts the three-axis gyroscope of the BMI160 inertial navigation component. The sampling rate of the three-axis gyroscope is set to 1600Hz, and the human body posture information of the pregnant woman is acquired through the three-axis gyroscope. Specifically, the triaxial gyroscope can collect the instantaneous angular velocity values of the human body in X, Y, Z three dimensions in real time, which are respectively recorded as g 1、g2 and g 3, and the vector space of the instantaneous angular velocity values of the human body in X, Y, Z three dimensions corresponding to the time t is represented by g t, so that g t corresponding to each time is represented as:
gt={g1,g2,g3}
The vector modulo length, l, is calculated by solving the Euclidean norm for three elements in g t:
The vector modulus length l is proportional to the human activity intensity and is used for measuring the instantaneous intensity of human activity at the current time t. In the embodiment of the invention, the vector module length l of the current time t is calculated in real time, and the newly generated vector module length l is continuously cached and combined into an unlimited time sequence, namely a vector module length time sequence d (x):
d(x)={l0,l1,......,∞}
Where l 0 is the vector modulo length of the first cache and l 1 is the vector modulo length of the second cache.
According to the sensor characteristics of the three-axis gyroscope, when the human body posture is in a calm state, the angular velocity values of the sensor on three axes X, Y, Z tend to be 0, otherwise, the greater the intensity of irregular movement of the human body posture is, the greater the absolute value of the corresponding angular velocity output by the sensor is.
Because the single measurement result of the gyroscope is an instantaneous value, but the instantaneous value cannot effectively represent the current activity state of the human body, finite-length time sequence data are required to be constructed by sliding a time window or rolling the time window, and the activity state of the human body is estimated by analyzing the data characteristics in each unit time window
At 1600Hz sampling rate, using the buffer 2048 data points as unit window data, adopting Short Time Fourier Transform (STFT) to perform local spectrum transformation of the unit window data according to the time sequence of the vector modulus long time sequence, combining the unit window data with a window function to perform spectrum conversion, and solving the local spectrum of the vector modulus long time sequence d (x). The Short Time Fourier Transform (STFT) formula is defined as:
m∈[0,1,......,∞]
n∈[0,1,......,1024]
Where H is the frame shift length, preferably 512, m is the number of frame shifts, d (mXH+k) is the vector modulo long time series, w (k) is the window function, j is the imaginary unit, and n is the argument of the windowed Fourier transform.
In order to mitigate the influence of spectrum leakage generated in the short-time Fourier transform process, a Hanning window is selected for calculation, the length of the Hanning window is 2048 points, and the window function w (k) is defined as:
When the independent variable m increases, the hanning window function is represented as sliding along the d (x) sequence, the sliding step is h=512 data points, after each sliding of the window function, the unit window data of the subscript corresponding to d (x) is subjected to one-dimensional matrix multiplication, the multiplied sequence result is solved, the multiplied sequence result is converted into a complex sequence with the imaginary number of 0, discrete Fourier transformation is carried out, and the single local spectrum transformation of d (x) is completed.
Compared with the traditional Fourier Transform (FT), the short-time Fourier transform (STFT) can reflect not only the spectrum characteristics of the original signal, but also the time-varying law of the spectrum of the original signal from a higher dimension. Therefore, by means of the frequency domain characteristics and the time domain characteristics of the STFT, the change intensity of the gyroscope data in the window time can be accurately counted, and the time boundary of the human body entering the calm state can be found. And then the activity state of the human body at the current moment is judged.
After the local spectrum transformation is completed, the local spectrum data are converted into power spectrum data, energy corresponding to all frequencies in the power spectrum data is summed, a power total value is calculated, if the power total value exceeds a threshold value, the current human body posture is judged to be in an active state, otherwise, the current human body posture is judged to be in a static state, and the pregnant woman is regarded as being in a calm state. When the current human body posture is judged to be in a static state, the impedance measurement of the human body skin is carried out, otherwise, the new data is cached by the sliding unit window continuously, the local frequency spectrum transformation of the next round of unit window data is carried out, and the corresponding human body posture judgment result is output.
Step 102, judging whether the pregnant woman is in a calm state according to the human body posture information, if so, starting an impedance measurement link, and measuring the skin impedance of the human body.
When the current human body posture is judged to be in a static state according to the total power value in the step 101, namely, when the pregnant woman is in a calm state, an impedance measurement link is started, an alternating current excitation current source is started, preset alternating current is applied to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, a differential signal between the P electrode and an N electrode is obtained through a differential amplifier to serve as an excitation response signal, wherein the alternating current excitation current source, the P electrode and the N electrode are connected in a ring mode to form a loop, the P electrode and the N electrode are connected with two input ends of the differential amplifier, the excitation response signal is preprocessed, the excitation response signal after the maternal biological electrical signal is filtered is obtained, the preprocessing comprises analog-digital conversion processing and filtering processing, discrete Fourier transformation is carried out on the excitation response signal after the maternal biological electrical signal is filtered, and a complex result of a fundamental frequency real part and an imaginary part is solved, so that the human skin impedance is obtained.
In typical human skin impedance measurement applications, the impedance measurement range is usually from the left upper limb to the right lower limb of the human body, most of the excitation signals applied to the human body are alternating current signals with the frequency being more than 50KHz and the current amplitude being 300 microamps, and the alternating current excitation frequency applied to human tissues is high, the relative intensity of the current is high, and the risk of human electrolyte flocculation is caused. When fetal electrocardiographic detection is carried out, the electrode pair is usually attached to the abdomen of a human body, the attaching positions of the electrode pair are similar, and the impedance generated by the skin between the P electrode and the N electrode and the interface between the electrode and the skin is measured. In the embodiment of the invention, a 250Hz low-frequency square wave excitation signal which is not easy to penetrate human skin tissue is adopted, and a 24nA weak current signal which is harmless to the skin tissue is adopted. By applying a constant frequency, constant current square wave excitation current signal to human skin, a signal response is induced by the complex impedance of human skin and the electrode-skin interface impedance characteristics, providing a response signal with complex impedance information for subsequent impedance measurements.
The skin impedance measurement method for fetal electrocardio monitoring provided by the embodiment of the invention utilizes a triaxial gyroscope to acquire human body posture information of a pregnant woman, judges whether the pregnant woman is in a calm state at present, starts an impedance measurement link when the pregnant woman is in the calm state, starts an alternating current excitation current source, applies preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, acquires a differential signal between the P electrode and an N electrode as an excitation response signal through a differential amplifier, preprocesses the excitation response signal to obtain an excitation response signal after filtering a maternal biological electrical signal, and then carries out discrete Fourier transform on the excitation response signal after filtering the maternal biological electrical signal to solve complex results of a real part and an imaginary part of fundamental frequency to obtain the skin impedance of the human body. The current of the alternating current excitation current source flows out from the N electrode after being led into the human skin from the P electrode, and under the condition of constant current, the response signal is caused between the P electrode and the N electrode due to the complex impedance characteristic of the electric signal of the human skin and the impedance characteristic of the gel bonding degree between the electrode pair and the skin, and the response intensity of the signal is inversely proportional to the complex impedance of the human skin and inversely proportional to the impedance of the interface of the electrode and the skin. The skin impedance measurement method for fetal electrocardio monitoring provided by the invention considers the imaginary part information of the excitation response signal, reduces the error of human skin impedance measurement, avoids the risk of thermal noise interference of the analog switch device and external noise interference, improves the accuracy of skin impedance measurement for fetal electrocardio monitoring, and solves the technical problems of ignoring the imaginary part information of the excitation response signal and noise interference caused by the analog switch device in the existing impedance measurement method at the interface of skin and electrode, thereby reducing the accuracy of impedance measurement.
In one embodiment, in the impedance measurement link, an ac excitation current source applies an excitation signal to human skin through a P-electrode pair, thereby causing a response signal. The response signal belongs to an aliasing signal and comprises an excitation response signal and a maternal bioelectric signal, wherein the excitation response signal and the maternal bioelectric signal are caused by complex impedance of human skin to an excitation signal. The aliasing signals are subjected to analog-to-digital conversion, excitation response signals are screened out through a digital band-pass filter, the digital band-pass filter filters 1024 sampling values cached in a rolling window, the conducting range is 240 Hz-260 Hz, and the filtering output result is a one-dimensional real sequence a (r) with the length equal to 1024. Constructing a finite-length complex sequence s (r) with an imaginary number of 0 by using a (r), and performing discrete Fourier transform on the finite-length complex sequence s (r), wherein a formula of the discrete Fourier transform is defined as follows:
where F (q) is the result of the discrete Fourier transform and s (r) is the finite length complex sequence of the input.
The result F (q) of the discrete fourier transform is in complex form, so that each frequency component contains real and imaginary information. Since the frequency of the signal applied to the human body by the ac excitation current source is 250Hz, complex information whose fundamental frequency component is 250Hz is extracted from the spectrum result F (q):
Z250=R250+jI250
Wherein j is an imaginary unit, R 250 is a real value of complex information with a fundamental wave frequency component of 250Hz extracted from a discrete Fourier transform result, the real value is characterized as an electrical impedance component in the measured human body complex impedance, I 250 is a virtual value of complex information with the fundamental wave frequency component of 250Hz extracted from the discrete Fourier transform result, and the real value is characterized as a capacitance component in the measured human body complex impedance.
Setting the differential amplification factor of the differential amplifier as b, and calculating the final human skin impedance by combining the differential amplification factor as b as follows:
Wherein Z is the impedance of human skin.
The spectrum conversion is carried out on the excitation response signal with limited length, and the impedance calculation is carried out through real number and imaginary number information, so that the accuracy and stability of the system for measuring the impedance of the human skin are improved.
In one embodiment, after the measurement of the skin impedance of the human body is completed, an alternating current excitation current source is turned off, a fetal electrocardiograph link is started, a maternal electrocardiograph signal is collected through a P electrode, an N electrode and other collecting electrodes, and maternal electrocardiograph signals are separated by combining the measured skin impedance of the human body, so that fetal electrocardiograph signals are obtained.
After the impedance measurement of the skin of the human body is completed, the alternating current excitation current source is turned off, the P electrode and the N electrode respectively collect human body potential signals, and the differential signals of the P electrode and the N electrode are mother fetal electrocardiosignals. After the alternating current excitation current source is turned off, the impedance measurement function is turned off, and excitation response signals caused by an excitation source on the electrode pairs do not exist, so that signals output by the analog-to-digital converter only comprise mother tyre mixed bioelectric signals, the mother tyre mixed bioelectric signals are filtered through an IIR digital filter, direct current components and power frequency interference in the signals are removed, and QRST signals of mother tyre electrocardiosignals are extracted. And according to the characteristics of the maternal-fetal electrocardiosignals, completing the maternal-fetal electrocardiosignal separation through a principal component analysis algorithm.
The circuit structure used for the skin impedance measurement method for fetal electrocardio monitoring provided by the embodiment of the invention is shown in figure 2, and comprises a gyroscope data analysis module, a P electrode, an N electrode, a triaxial gyroscope, an alternating current excitation current source, a differential amplifier, an analog-to-digital converter, a right leg driving circuit module, an information interaction interface and an impedance measurement module, wherein the triaxial gyroscope is connected with the gyroscope data analysis module, the gyroscope data analysis module is connected with the alternating current excitation current source, the positive end of the alternating current excitation current source is connected with the P electrode, the negative end is connected with the N electrode, the P electrode and the N electrode are connected with two input ends of the differential amplifier, the output end of the differential amplifier is connected with the input end of the analog-to-digital converter and the right leg driving circuit module, the right leg driving circuit module comprises a right leg driving circuit and an E electrode, the three-axis gyroscope, the P electrode and the N electrode are attached to the abdomen of a pregnant woman, the output end of the analog-to-digital converter is connected with the impedance measurement module, the output end of the impedance measurement module is connected with the information interaction interface, the output end of the impedance measurement module is also connected with the input end of the alternating current excitation current source, and the impedance measurement module sends a control signal to the alternating current excitation current source after completing impedance measurement, so that the alternating current excitation current source is closed. The circuit structure of fig. 2 operates on the principle that:
The gyroscope is used for detecting the posture change of the human body, when the human body posture is detected to be in an active state, impedance measurement is not carried out, and otherwise, the impedance measurement is carried out.
The alternating current excitation current source is used for generating an alternating current excitation signal harmless to a human body.
The system input electrode comprises at least one set of electrode pairs, each pair comprising a P electrode and an N electrode. During impedance measurement, current is applied to human skin from the P electrode and from the N electrode loop. When fetal electrocardio detection is carried out, an alternating current excitation power supply is turned off, and the P electrode and the N electrode respectively collect human body potential signals.
The differential amplifier is used for amplifying differential signals formed between each group of P electrodes and N electrodes. During impedance measurement, the differential signals of the P electrode and the N electrode are excitation response signals due to the complex impedance characteristic of human skin, and during fetal electrocardio detection, the differential signals of the P electrode and the N electrode are mother fetal mixed electrocardio signals.
The analog-to-digital converter is used for carrying out quantization sampling on the signals after differential amplification.
The impedance measurement module comprises a digital band-pass filter and an impedance analysis unit, and the fetal electrocardiograph measurement module comprises an IIR digital filter and a mother fetal electrocardiograph separation module.
The digital band-pass filter is used for filtering when detecting impedance, so that interference of electrocardiosignals of a human body is reduced.
The IIR digital filter is used for filtering when detecting fetal electrocardio, and reducing interference of direct current and high-frequency signals.
The impedance analysis unit is used for analyzing a response signal formed by a human body to the excitation source and calculating an impedance value from the response signal;
The number of the P-electrode and the N-electrode and other collecting electrodes (a plurality of electrodes are used in fetal electrocardiographic detection, including the P-electrode and the N-electrode in fig. 2, and further including other collecting electrodes, and the number of the other collecting electrodes is selected according to an actual monitoring scene, which is not limited herein) are used for collecting maternal electrocardiographic signals in fetal electrocardiographic detection.
The maternal-fetal electrocardiosignal separation module is used for analyzing electrocardiosignals acquired from a human body, separating maternal-fetal electrocardiosignals and extracting fetal electrocardiosignals;
the system comprises a right leg driving circuit, wherein the right leg driving circuit feeds common mode signals to a human body through an E electrode and is used for inhibiting common mode interference caused by human body functions during electrocardiograph measurement.
The information interaction interface is used for issuing the calculated real-time impedance data or fetal electrocardio data outwards. The information interaction interface adopts a wireless mode of BLE low-power consumption Bluetooth to form a man-machine interaction channel.
When measuring the impedance formed by the electrode pair (P electrode and N electrode) and the skin, the impedance measurement is divided into a front-stage system and a rear-stage system, wherein the front-stage system is an alternating current excitation source plus electrode (front input), and the rear-stage system is amplification plus analog-digital conversion plus signal analysis (rear-stage processing). The front stage system introduces an alternating current excitation current source, the excitation current source generates square wave excitation signals with constant on-off frequency and constant current, the excitation signals are applied to the P electrode, and the current flows out of the N electrode after passing through the skin of a human body and forms a current signal loop with the excitation current source. The differential amplifier in the post-stage system is used for amplifying signals of the P electrode and the N electrode, sampling and quantizing are carried out by the analog-to-digital converter, and after the quantized result is filtered by the digital band-pass filter, the impedance analysis unit carries out data analysis on the filtered result. The current flows out from the N electrode after being led into the human skin from the P electrode, and under the condition of constant current, a response signal is caused between the P electrode and the N electrode due to the complex impedance characteristic of the electric signal of the human skin and the impedance characteristic of the gel bonding degree between the electrode pair and the skin, and the response intensity of the signal is inversely proportional to the complex impedance of the human skin and inversely proportional to the impedance of the interface of the electrode and the skin.
The latter stage system carries out differential signal amplification and analog-to-digital conversion sampling on the electrode pair attached to the skin of the human body. The common mode interference of bioelectric signals can be effectively reduced by using a differential input mode, and the differential signals are amplified and then are subjected to signal quantization sampling by using an analog-to-digital converter. According to the rolling window buffer mechanism, a first-in first-out rule is followed, the length of 1024 sampling values are buffered in a rolling way, and the buffered time sequence data not only comprises excitation response signal data caused by the application of excitation current to the skin, but also comprises electrocardiosignal data of a human body.
Typically, the spectral energy range of the human body electrocardiosignal is less than 30Hz, and the spectral energy of the excitation signal applied by the front stage is concentrated at 250Hz and the frequency multiplication thereof. In order to demodulate and separate response data with complex impedance information from time sequence data of the aliasing electrocardiac, in this embodiment, a digital band-pass filter is constructed by adopting a digital signal processing technology, and the response data with complex impedance information is screened out from the sampled aliasing digital signal by the filter. The passband of the digital band-pass filter is 240 Hz-260 Hz.
For ease of understanding, referring to fig. 2 and 3, an embodiment of a skin impedance measurement system for fetal electrocardiograph monitoring is provided in the present invention, which includes a gyroscope data analysis module, a P-electrode, an N-electrode, a tri-axis gyroscope, an ac excitation current source, a differential amplifier, an analog-to-digital converter, a right leg driving circuit module, an information interaction interface, and an impedance measurement module;
the three-axis gyroscope is connected with the gyroscope data analysis module;
the gyroscope data analysis module is connected with an alternating current excitation current source;
The positive end of the alternating current excitation current source is connected with the P electrode, the negative end of the alternating current excitation current source is connected with the N electrode, the P electrode and the N electrode are connected with two input ends of the differential amplifier, the output end of the differential amplifier is connected with the input end of the analog-to-digital converter and the right leg driving circuit module, the right leg driving circuit module comprises a right leg driving circuit and an E electrode, the E electrode is attached to the right leg of a pregnant woman, and the three-axis gyroscope, the P electrode and the N electrode are attached to the abdomen of the pregnant woman;
the output end of the analog-to-digital converter is connected with the impedance measurement module, and the output end of the impedance measurement module is connected with the information interaction interface;
When the skin impedance measuring system for fetal electrocardio monitoring is used for measuring the skin impedance of a human body, the three-axis gyroscope collects the human body posture information and sends the human body posture information to the gyroscope data analysis module, the gyroscope data analysis module judges whether the pregnant woman is in a calm state according to the human body posture information, if so, an impedance measuring link is started, and the skin impedance of the human body is measured;
The impedance measurement link includes:
Starting an alternating current excitation current source, applying preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, and obtaining a differential signal between the P electrode and an N electrode through a differential amplifier to serve as an excitation response signal, wherein the alternating current excitation current source, the P electrode and the N electrode are connected in a ring shape to form a loop, and the P electrode and the N electrode are connected with two input ends of the differential amplifier;
Preprocessing the excitation response signal through an impedance measurement module to obtain an excitation response signal after the maternal biological electrical signal is filtered, wherein the preprocessing comprises analog-to-digital conversion processing and filtering processing;
And performing discrete Fourier transform on the excitation response signal after the maternal biological electrical signal is filtered by an impedance measurement module, and solving complex results of a real part and an imaginary part of fundamental frequency to obtain the skin impedance of the human body.
Also included are fetal electrocardiograph modules and other acquisition electrodes (not labeled in fig. 2 and 3);
the other acquisition electrodes are connected with the input end of the differential amplifier, the input end of the fetal electrocardiograph measurement module is connected with the output end of the analog-to-digital converter, and is also connected with the impedance measurement module, and the output end of the fetal electrocardiograph measurement module is connected with the information interaction interface;
The fetal electrocardiosignal module is used for measuring fetal electrocardiosignals after finishing the impedance measurement of the skin of a human body;
After the measurement of the skin impedance of the human body is completed, the impedance measurement module sends a control signal to the alternating current excitation current source, so that the alternating current excitation current source is turned off, the maternal-fetal electrocardiosignals are collected through the P electrode, the N electrode and other collecting electrodes, and the maternal-fetal electrocardiosignals are subjected to maternal-fetal electrocardiosignal separation by combining the skin impedance of the human body measured by the impedance measurement module, so that the fetal electrocardiosignals are obtained.
The gyroscope data analysis module is specifically used for:
Acquiring the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction acquired by the triaxial gyroscope in real time, and calculating the vector module length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment;
Caching the vector modular length at each moment to obtain a vector modular length time sequence;
carrying out local spectrum transformation of unit window data on the vector modulus long-time sequence by adopting short-time Fourier transformation, and converting the local spectrum data into power spectrum data;
And summing the power spectrum data, calculating a power total value, judging that the current human body posture is in an active state if the power total value exceeds a threshold value, and judging that the current human body posture is in a static state if the power total value exceeds the threshold value.
Calculating the vector modular length of each moment according to the instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment, including:
The instantaneous angular velocity value of the human body in the X, Y, Z dimension direction at each moment is represented as a vector space, the Euclidean norm of the vector space at each moment is calculated, and the vector modular length at each moment is obtained.
The short-time Fourier transform formula for carrying out local frequency spectrum transformation of unit window data on the vector modulus long-time sequence by adopting short-time Fourier transform is as follows:
m∈[0,1,......,∞]
n∈[0,1,......,1024]
Where H is the frame shift length, m is the number of frame shifts, d (mxh+k) is the vector modulo long time series, w (k) is the window function, j is the imaginary unit, and n is the argument of the windowed fourier transform.
The preset alternating current is a square wave excitation current of 250Hz24 nA.
Performing discrete Fourier transform on the excitation response signal after the maternal biological electrical signal is filtered, solving complex results of a real part and an imaginary part of fundamental frequency, and obtaining skin impedance of a human body, wherein the method comprises the following steps:
constructing an excitation response signal with filtered maternal-fetal bioelectric signals into a finite-length complex sequence with an imaginary number of 0, and performing discrete Fourier transform on the finite-length complex sequence with the imaginary number of 0, wherein the formula of the discrete Fourier transform is as follows:
wherein F (q) is the result of discrete Fourier transform, s (r) is the input finite-length complex sequence;
The complex information with the fundamental wave frequency component of 250Hz is taken out from the discrete Fourier transform result, and the human skin impedance is calculated by combining the differential amplification factor of the differential amplifier, wherein the calculation formula is as follows:
Wherein Z is the human skin impedance, b is the differential amplification factor of the differential amplifier, R 250 is the real value of the complex information with the fundamental frequency component of 250Hz extracted from the discrete Fourier transform result, and I 250 is the imaginary value of the complex information with the fundamental frequency component of 250Hz extracted from the discrete Fourier transform result.
The skin impedance measurement system for fetal electrocardio monitoring provided by the invention utilizes the three-axis gyroscope to acquire human body posture information of a pregnant woman, judges whether the pregnant woman is in a calm state at present, starts an impedance measurement link when the pregnant woman is in the calm state, starts an alternating current excitation current source, applies preset alternating current to the skin surface of the pregnant woman through a P electrode connected with the alternating current excitation current source, acquires a differential signal between the P electrode and an N electrode as an excitation response signal through a differential amplifier, preprocesses the excitation response signal to obtain the excitation response signal after filtering the maternal biological electrical signal, and then carries out discrete Fourier transform on the excitation response signal after filtering the maternal biological electrical signal, so as to solve complex results of a real part and an imaginary part of fundamental frequency, and obtain human skin impedance.
The principle of the skin impedance measurement system for fetal electrocardiograph monitoring provided in the embodiment of the present invention is the same as that of the skin impedance measurement method for fetal electrocardiograph monitoring in the embodiment of the skin impedance measurement method for fetal electrocardiograph monitoring, and the skin impedance measurement system for fetal electrocardiograph monitoring is not described in detail herein.
While the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those skilled in the art that the foregoing embodiments may be modified or equivalents may be substituted for some of the features thereof, and that the modifications or substitutions do not depart from the spirit and scope of the embodiments of the invention.

Claims (9)

1.一种用于胎儿心电监测的皮肤阻抗测量方法,其特征在于,包括:1. A skin impedance measurement method for fetal electrocardiogram monitoring, comprising: 通过安装在孕妇身上的三轴陀螺仪采集人体姿态信息;The body posture information is collected by a three-axis gyroscope installed on the pregnant woman; 根据人体姿态信息判断孕妇是否处于平静状态,若是,则启动阻抗测量环节,测量人体皮肤阻抗;Judging whether the pregnant woman is in a calm state according to the human posture information, if so, starting the impedance measurement link to measure the human skin impedance; 阻抗测量环节包括:The impedance measurement process includes: 启动交流激励电流源,通过与交流激励电流源连接的P电极对孕妇皮肤表面施加预置交流电流,通过差分放大器获取P电极和N电极之间的差分信号作为激励响应信号,其中,交流激励电流源、P电极和N电极环形连接形成回路,P电极和N电极与差分放大器的两个输入端连接;The AC excitation current source is started, and a preset AC current is applied to the skin surface of the pregnant woman through a P electrode connected to the AC excitation current source, and a differential signal between the P electrode and the N electrode is obtained as an excitation response signal through a differential amplifier, wherein the AC excitation current source, the P electrode, and the N electrode are connected in a ring to form a loop, and the P electrode and the N electrode are connected to two input terminals of the differential amplifier; 对激励响应信号进行预处理,得到滤除母胎生物电信号后的激励响应信号,预处理包括模数转换处理和滤波处理;Preprocessing the stimulus response signal to obtain the stimulus response signal after filtering out the maternal and fetal bioelectric signals, the preprocessing including analog-to-digital conversion processing and filtering processing; 对滤除母胎生物电信号后的激励响应信号进行离散傅里叶变换,求解出基波频率实部和虚部的复数结果,得到人体皮肤阻抗;Perform discrete Fourier transform on the excitation response signal after filtering out the maternal and fetal bioelectrical signals, solve the complex results of the real and imaginary parts of the fundamental frequency, and obtain the human skin impedance; 对滤除母胎生物电信号后的激励响应信号进行离散傅里叶变换,求解出基波频率实部和虚部的复数结果,得到人体皮肤阻抗,包括:After filtering out the maternal and fetal bioelectric signals, the excitation response signal is subjected to discrete Fourier transform to solve the complex results of the real and imaginary parts of the fundamental frequency and obtain the human skin impedance, including: 对滤除母胎生物电信号后的激励响应信号构造成虚数为0的有限长复数序列,对虚数为0的有限长复数序列进行离散傅里叶变换;The excitation response signal after filtering out the maternal and fetal bioelectrical signals is constructed into a finite length complex sequence with an imaginary number of 0, and the finite length complex sequence with an imaginary number of 0 is subjected to discrete Fourier transform; 从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息,结合差分放大器的差分放大倍数,计算人体皮肤阻抗,计算公式为:The complex information with the fundamental frequency component of 250Hz is taken out from the discrete Fourier transform result, and the human skin impedance is calculated by combining the differential gain of the differential amplifier. The calculation formula is: 其中,Z为人体皮肤阻抗,b为差分放大器的差分放大倍数,R250为从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息的实数值,I250为从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息的虚数值。Wherein, Z is the human skin impedance, b is the differential gain of the differential amplifier, R 250 is the real value of the complex information with the fundamental frequency component of 250 Hz taken out from the discrete Fourier transform result, and I 250 is the imaginary value of the complex information with the fundamental frequency component of 250 Hz taken out from the discrete Fourier transform result. 2.根据权利要求1所述的用于胎儿心电监测的皮肤阻抗测量方法,其特征在于,还包括:2. The skin impedance measurement method for fetal electrocardiogram monitoring according to claim 1, further comprising: 当完成人体皮肤阻抗测量后,关闭交流激励电流源,启动胎儿心电测量环节,通过P电极、N电极以及其他采集电极采集母胎心电信号,结合测量的人体皮肤阻抗对母胎心电信号进行母胎心电分离,得到胎儿心电信号。After completing the human skin impedance measurement, turn off the AC excitation current source, start the fetal ECG measurement link, collect maternal and fetal ECG signals through the P electrode, N electrode and other collection electrodes, and perform maternal and fetal ECG separation on the maternal and fetal ECG signals in combination with the measured human skin impedance to obtain the fetal ECG signal. 3.根据权利要求1所述的用于胎儿心电监测的皮肤阻抗测量方法,其特征在于,通过安装在孕妇身上的三轴陀螺仪采集人体姿态信息,包括:3. The skin impedance measurement method for fetal electrocardiogram monitoring according to claim 1 is characterized in that the human body posture information is collected by a three-axis gyroscope installed on the pregnant woman, comprising: 实时获取安装在孕妇身上的三轴陀螺仪采集的人体在X、Y、Z维度方向上的瞬时角速度值,根据每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值计算每个时刻的向量模长;The instantaneous angular velocity values of the human body in the X, Y, and Z directions collected by the three-axis gyroscope installed on the pregnant woman are obtained in real time, and the vector modulus at each moment is calculated according to the instantaneous angular velocity values of the human body in the X, Y, and Z directions at each moment; 对每个时刻的向量模长进行缓存,得到向量模长时间序列;Cache the vector modulus length at each moment to obtain a long vector modulus time series; 采用短时傅里叶变换对向量模长时间序列进行单位窗口数据的局部频谱变换,将局部频谱数据转换为功率谱数据;Short-time Fourier transform is used to perform local spectrum transformation of unit window data on the vector modulus long time series, and the local spectrum data is converted into power spectrum data; 对功率谱数据进行求和,计算功率总值,若功率总值超过阈值,则判定当前人体姿态为活动状态,否则,判定当前人体姿态为静止状态。The power spectrum data is summed to calculate the total power value. If the total power value exceeds the threshold, the current human body posture is determined to be in an active state, otherwise, the current human body posture is determined to be in a static state. 4.根据权利要求3所述的用于胎儿心电监测的皮肤阻抗测量方法,其特征在于,根据每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值计算每个时刻的向量模长,包括:4. The skin impedance measurement method for fetal electrocardiogram monitoring according to claim 3 is characterized in that the vector modulus at each moment is calculated according to the instantaneous angular velocity value of the human body in the X, Y, and Z dimensional directions at each moment, comprising: 将每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值表征成向量空间,计算出每个时刻向量空间的Euclidean范数,得到每个时刻的向量模长。The instantaneous angular velocity values of the human body in the X, Y, and Z dimensions at each moment are represented as a vector space, the Euclidean norm of the vector space at each moment is calculated, and the vector modulus at each moment is obtained. 5.根据权利要求1所述的用于胎儿心电监测的皮肤阻抗测量方法,其特征在于,预置交流电流为250Hz24nA的方波激励电流。5. The skin impedance measurement method for fetal electrocardiogram monitoring according to claim 1, characterized in that the preset alternating current is a square wave excitation current of 250Hz24nA. 6.一种用于胎儿心电监测的皮肤阻抗测量系统,其特征在于,包括陀螺仪数据分析模块、P电极、N电极、三轴陀螺仪、交流激励电流源、差分放大器、模数转换器、右腿驱动电路模块、信息交互接口和阻抗测量模块;6. A skin impedance measurement system for fetal electrocardiogram monitoring, characterized in that it includes a gyroscope data analysis module, a P electrode, an N electrode, a three-axis gyroscope, an AC excitation current source, a differential amplifier, an analog-to-digital converter, a right leg drive circuit module, an information interaction interface and an impedance measurement module; 三轴陀螺仪与陀螺仪数据分析模块连接;The three-axis gyroscope is connected to the gyroscope data analysis module; 陀螺仪数据分析模块与交流激励电流源连接;The gyroscope data analysis module is connected to the AC excitation current source; 交流激励电流源的正端与P电极连接,负端与N电极连接,P电极和N电极与差分放大器的两个输入端连接,差分放大器的输出端连接模数转换器的输入端与右腿驱动电路模块,右腿驱动电路模块包括右腿驱动电路和E电极,三轴陀螺仪、E电极、P电极和N电极贴覆在孕妇腹部;The positive end of the AC excitation current source is connected to the P electrode, and the negative end is connected to the N electrode. The P electrode and the N electrode are connected to two input ends of a differential amplifier. The output end of the differential amplifier is connected to the input end of an analog-to-digital converter and a right leg driving circuit module. The right leg driving circuit module includes a right leg driving circuit and an E electrode. The three-axis gyroscope, the E electrode, the P electrode, and the N electrode are attached to the abdomen of the pregnant woman. 模数转换器的输出端连接阻抗测量模块,阻抗测量模块的输出端连接信息交互接口;The output end of the analog-to-digital converter is connected to the impedance measurement module, and the output end of the impedance measurement module is connected to the information interaction interface; 当用于胎儿心电监测的皮肤阻抗测量系统用于测量人体皮肤阻抗时,三轴陀螺仪采集人体姿态信息发送给陀螺仪数据分析模块,陀螺仪数据分析模块根据人体姿态信息判断孕妇是否处于平静状态,若是,则启动阻抗测量环节,测量人体皮肤阻抗;When the skin impedance measurement system for fetal electrocardiogram monitoring is used to measure human skin impedance, the three-axis gyroscope collects human posture information and sends it to the gyroscope data analysis module. The gyroscope data analysis module determines whether the pregnant woman is in a calm state based on the human posture information. If so, the impedance measurement link is started to measure the human skin impedance. 阻抗测量环节包括:The impedance measurement process includes: 启动交流激励电流源,通过与交流激励电流源连接的P电极对孕妇皮肤表面施加预置交流电流,通过差分放大器获取P电极和N电极之间的差分信号作为激励响应信号,其中,交流激励电流源、P电极和N电极环形连接形成回路,P电极和N电极与差分放大器的两个输入端连接;The AC excitation current source is started, and a preset AC current is applied to the skin surface of the pregnant woman through a P electrode connected to the AC excitation current source, and a differential signal between the P electrode and the N electrode is obtained as an excitation response signal through a differential amplifier, wherein the AC excitation current source, the P electrode, and the N electrode are connected in a ring to form a loop, and the P electrode and the N electrode are connected to two input terminals of the differential amplifier; 通过阻抗测量模块对激励响应信号进行预处理,得到滤除母胎生物电信号后的激励响应信号,预处理包括模数转换处理和滤波处理;Preprocessing the excitation response signal through the impedance measurement module to obtain the excitation response signal after filtering out the maternal and fetal bioelectric signals, the preprocessing including analog-to-digital conversion processing and filtering processing; 通过阻抗测量模块对滤除母胎生物电信号后的激励响应信号进行离散傅里叶变换,求解出基波频率实部和虚部的复数结果,得到人体皮肤阻抗;The impedance measurement module performs discrete Fourier transformation on the excitation response signal after filtering out the maternal and fetal bioelectric signals, and solves the complex results of the real and imaginary parts of the fundamental frequency to obtain the human skin impedance; 对滤除母胎生物电信号后的激励响应信号进行离散傅里叶变换,求解出基波频率实部和虚部的复数结果,得到人体皮肤阻抗,包括:After filtering out the maternal and fetal bioelectric signals, the excitation response signal is subjected to discrete Fourier transform to solve the complex results of the real and imaginary parts of the fundamental frequency and obtain the human skin impedance, including: 对滤除母胎生物电信号后的激励响应信号构造成虚数为0的有限长复数序列,对虚数为0的有限长复数序列进行离散傅里叶变换;The excitation response signal after filtering out the maternal and fetal bioelectrical signals is constructed into a finite length complex sequence with an imaginary number of 0, and the finite length complex sequence with an imaginary number of 0 is subjected to discrete Fourier transform; 从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息,结合差分放大器的差分放大倍数,计算人体皮肤阻抗,计算公式为:The complex information with the fundamental frequency component of 250Hz is taken out from the discrete Fourier transform result, and the human skin impedance is calculated by combining the differential gain of the differential amplifier. The calculation formula is: 其中,Z为人体皮肤阻抗,b为差分放大器的差分放大倍数,R250为从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息的实数值,I250为从离散傅里叶变换结果中取出基波频率分量为250Hz的复数信息的虚数值。Wherein, Z is the human skin impedance, b is the differential gain of the differential amplifier, R 250 is the real value of the complex information with the fundamental frequency component of 250 Hz taken out from the discrete Fourier transform result, and I 250 is the imaginary value of the complex information with the fundamental frequency component of 250 Hz taken out from the discrete Fourier transform result. 7.根据权利要求6所述的用于胎儿心电监测的皮肤阻抗测量系统,其特征在于,还包括胎儿心电测量模块和其他采集电极;7. The skin impedance measurement system for fetal ECG monitoring according to claim 6, characterized in that it also includes a fetal ECG measurement module and other acquisition electrodes; 其他采集电极与差分放大器的输入端连接,胎儿心电测量模块的输入端与模数转换器的输出端连接,还与阻抗测量模块连接,胎儿心电测量模块的输出端与信息交互接口连接;The other acquisition electrodes are connected to the input end of the differential amplifier, the input end of the fetal electrocardiogram measurement module is connected to the output end of the analog-to-digital converter and the impedance measurement module, and the output end of the fetal electrocardiogram measurement module is connected to the information interaction interface; 胎儿心电测量模块用于在完成人体皮肤阻抗测量后测量胎儿心电信号;The fetal ECG measurement module is used to measure the fetal ECG signal after completing the human skin impedance measurement; 在完成人体皮肤阻抗测量后,阻抗测量模块向交流激励电流源发送控制信号,使得交流激励电流源关闭,通过P电极、N电极以及其他采集电极采集母胎心电信号,结合阻抗测量模块测量的人体皮肤阻抗对母胎心电信号进行母胎心电分离,得到胎儿心电信号。After completing the human skin impedance measurement, the impedance measurement module sends a control signal to the AC excitation current source to turn off the AC excitation current source, and collects maternal and fetal ECG signals through the P electrode, N electrode and other collection electrodes. The maternal and fetal ECG signals are separated by combining the human skin impedance measured by the impedance measurement module to obtain the fetal ECG signal. 8.根据权利要求6所述的用于胎儿心电监测的皮肤阻抗测量系统,其特征在于,陀螺仪数据分析模块具体用于:8. The skin impedance measurement system for fetal electrocardiogram monitoring according to claim 6, wherein the gyroscope data analysis module is specifically used for: 实时获取三轴陀螺仪采集的人体在X、Y、Z维度方向上的瞬时角速度值,根据每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值计算每个时刻的向量模长;Acquire the instantaneous angular velocity values of the human body in the X, Y, and Z directions collected by the three-axis gyroscope in real time, and calculate the vector modulus at each moment according to the instantaneous angular velocity values of the human body in the X, Y, and Z directions at each moment; 对每个时刻的向量模长进行缓存,得到向量模长时间序列;Cache the vector modulus length at each moment to obtain a long vector modulus time series; 采用短时傅里叶变换对向量模长时间序列进行单位窗口数据的局部频谱变换,将局部频谱数据转换为功率谱数据;Short-time Fourier transform is used to perform local spectrum transformation of unit window data on the vector modulus long time series, and the local spectrum data is converted into power spectrum data; 对功率谱数据进行求和,计算功率总值,若功率总值超过阈值,则判定当前人体姿态为活动状态,否则,判定当前人体姿态为静止状态。The power spectrum data is summed to calculate the total power value. If the total power value exceeds the threshold, the current human body posture is determined to be in an active state, otherwise, the current human body posture is determined to be in a static state. 9.根据权利要求8所述的用于胎儿心电监测的皮肤阻抗测量系统,其特征在于,根据每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值计算每个时刻的向量模长,包括:9. The skin impedance measurement system for fetal electrocardiogram monitoring according to claim 8, characterized in that the vector modulus at each moment is calculated according to the instantaneous angular velocity value of the human body in the X, Y, and Z dimensional directions at each moment, comprising: 将每个时刻的人体在X、Y、Z维度方向上的瞬时角速度值表征成向量空间,计算出每个时刻向量空间的Euclidean范数,得到每个时刻的向量模长。The instantaneous angular velocity values of the human body in the X, Y, and Z dimensions at each moment are represented as a vector space, the Euclidean norm of the vector space at each moment is calculated, and the vector modulus at each moment is obtained.
CN202210909613.1A 2022-07-29 2022-07-29 A skin impedance measurement method and system for fetal electrocardiogram monitoring Active CN115153580B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210909613.1A CN115153580B (en) 2022-07-29 2022-07-29 A skin impedance measurement method and system for fetal electrocardiogram monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210909613.1A CN115153580B (en) 2022-07-29 2022-07-29 A skin impedance measurement method and system for fetal electrocardiogram monitoring

Publications (2)

Publication Number Publication Date
CN115153580A CN115153580A (en) 2022-10-11
CN115153580B true CN115153580B (en) 2025-05-27

Family

ID=83476521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210909613.1A Active CN115153580B (en) 2022-07-29 2022-07-29 A skin impedance measurement method and system for fetal electrocardiogram monitoring

Country Status (1)

Country Link
CN (1) CN115153580B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024103214A1 (en) * 2022-11-14 2024-05-23 深圳市韶音科技有限公司 Signal collection system
CN117243585B (en) * 2023-11-10 2024-04-16 深圳市奋达智能技术有限公司 Noninvasive heart displacement measurement method, noninvasive heart displacement measurement system, computer equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105232034A (en) * 2014-06-04 2016-01-13 中国科学院半导体研究所 Electrode and circuit multiplex structure for detecting electrocardiograph signals and human impedance
KR101849857B1 (en) * 2016-12-12 2018-04-18 김준호 Wearable living body diagnosis device
CN109793516A (en) * 2019-01-21 2019-05-24 浙江大学 A kind of skin electrical impedance flexible detection device and skin electrical impedance detection method
CN111012312A (en) * 2019-12-25 2020-04-17 中国科学院合肥物质科学研究院 Portable Parkinson's disease bradykinesia monitoring intervention device and method
CN111166308A (en) * 2020-01-23 2020-05-19 北京津发科技股份有限公司 A multi-parameter physiological signal acquisition method of finger parts

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101194834A (en) * 2006-12-05 2008-06-11 重庆博恩富克医疗设备有限公司 Bio-electrical impedance measuring method and apparatus
WO2013075270A1 (en) * 2011-11-25 2013-05-30 Yang Chang-Ming Object, method, and system for detecting heartbeat or whether or not electrodes are in proper contact
CN103196591B (en) * 2013-03-07 2015-10-21 同济大学 A kind of structural loads recognition methods based on regularization and svd
CN109222934B (en) * 2018-10-30 2019-06-28 北京康加科技有限公司 System and method for bio-electrical impedance coherence measurement
KR102810731B1 (en) * 2019-12-11 2025-05-21 삼성전자주식회사 Method for processing biometric signal, the electronic device and storage medium therefor
CN113827247B (en) * 2021-09-27 2022-07-15 中国科学院深圳先进技术研究院 Spread spectrum modulation electrode contact impedance online measurement device and method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105232034A (en) * 2014-06-04 2016-01-13 中国科学院半导体研究所 Electrode and circuit multiplex structure for detecting electrocardiograph signals and human impedance
KR101849857B1 (en) * 2016-12-12 2018-04-18 김준호 Wearable living body diagnosis device
CN109793516A (en) * 2019-01-21 2019-05-24 浙江大学 A kind of skin electrical impedance flexible detection device and skin electrical impedance detection method
CN111012312A (en) * 2019-12-25 2020-04-17 中国科学院合肥物质科学研究院 Portable Parkinson's disease bradykinesia monitoring intervention device and method
CN111166308A (en) * 2020-01-23 2020-05-19 北京津发科技股份有限公司 A multi-parameter physiological signal acquisition method of finger parts

Also Published As

Publication number Publication date
CN115153580A (en) 2022-10-11

Similar Documents

Publication Publication Date Title
Berwal et al. Motion artifact removal in ambulatory ECG signal for heart rate variability analysis
Nayak et al. Filtering techniques for ECG signal processing
CN115153580B (en) A skin impedance measurement method and system for fetal electrocardiogram monitoring
JP7277970B2 (en) Electrocardiogram measurement method and system using wearable device
KR20120094870A (en) System for measure of electromyogram by suit type electromyogram sensors and analisys method of rehabilitation using thereof
CN110897631B (en) Real-time pregnancy monitoring device and method
CN106264504A (en) Noninvasive Blood Pressure Measurement System based on finger arteriogram and method
CN202179534U (en) Subminiature atrial fibrillation detector
KR20120094857A (en) Apparatus for measure of bio signal and method for analysis of rehabilitation training by suit type bio sensors
CN113974576B (en) Sleep quality monitoring system and monitoring method based on magnetocardiogram
CN112450899A (en) Self-calibration continuous blood pressure measuring device and method
CN113100776B (en) A fatigue monitoring system and method for fusing myoelectric and electrocardiographic signals
CN201847678U (en) Waist belt type wireless heart rate detection device
CN115192005A (en) A portable electrocardiogram and seismogram combined acquisition system and data processing method
CN109602395B (en) Noninvasive multichannel arterial system detection method and device
CN111714123A (en) A system and method for detecting electromyographic signals on the surface of human lower back
Yama et al. Development of a wireless capacitive sensor for ambulatory ECG monitoring over clothes
Jung et al. Development of wearable ECG measurement system using EMD for motion artifact removal
CN103876731B (en) A kind of fetal electrocardiogram signal extraction device and method
CN111134641A (en) A sleep monitoring chip system and sleep monitoring chip
DongseoUniversity The deployment of novel techniques for mobile ECG monitoring
CN211022624U (en) Continuous blood pressure measuring device based on fingertip pulse wave conduction time
JP2002078695A (en) Electrocardiogram measuring instrument
CN114795235B (en) Single-lead electrocardiogram monitoring method and system based on morphological contour algorithm
CN112535483B (en) Heart sympathetic activity detection device based on high-frequency electrocardio

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant