# Cardiopulmonary Signal Detection Based on Magnetic Induction

Authors: Dan Yang; Bin Xu; Hongyu Qiao; Xu Wang
Journal: Journal of Sensors
 (2017)

---

## Abstract

Cardiopulmonary signal contains vital and rich physiological information, which is very useful for clinical diagnosis and home healthcare. Convectional monitoring methods still need uncomfortable sensors or are expensive for nonprofessional common consumers. In this paper, we proposed a cardiopulmonary signal detection method based on electromagnetic induction. Based on thoracic volume variation affecting the biological impedance, which can be detected by magnetic induction, our idea is to develop a not complicated and practical measuring platform and exploit the intrinsic properties of cardiopulmonary signals with respiratory and heart rate. We have some theoretical analysis about the relationship of thoracic volume with biological impedance, sensor-head parameters, and the optimal measuring position. Then a whole measuring system has been designed and evaluated. The respiratory and heart rate obtained from the proposed method are not significantly different from the reference method (noncontact BCG).


---

## Body

## 1. Introduction

According to the World Health Organization (WHO), cardiovascular diseases (CVD) are the leading causes of death and disability in the world. An estimated 17.5 million people died from CVD in 2016, representing 30% of all global deaths [1]. Meanwhile the elderly population has increased to almost 810 million in 2012. In 2050, the number of aged people (60 and above) is about to reach a staggering 2 billion [2]. These cause the needs to monitor patient health status while they are out of the hospital to increase the comfort and wellbeing of the patients which can be in a familiar environment and at the same time decrease the significant cost on healthcare. Effective disease management through continuous monitoring and information fusion of vital physiological signals in patients is viewed as the key mechanism for the mentioned problems.Cardiopulmonary signal includes respiratory and heart rate information that can be monitored to detect acute changes in personal condition (e.g., apnea) and potentially provide an early warning of impending life-threatening deterioration. Current methods require sensors attached to the body or not robust to subject motion. The most widely used method for detecting cardiopulmonary signal is electrocardiography (ECG) which provides useful information about the cardiovascular system. ECG captures the electrical activity of the heart as a time series data with Ag/AgCl electrodes attached to exposed bare skin; this will cause restraint and produce some burdens to subjects [3]. Airflow is another popular measuring vital sign (respiratory) method with common sensors fixed at the nasal cavity or the chest [4]. Video camera methods run the gamut ranging from high-end expensive camera systems that track bright objects placed on the chest and used in CT image gating to low-cost vital signs “apps” for mobile devices that give heart rate by detecting facial flushing with each beat and respiration rate by chest movement [5]. Thermal cameras detect air movement and temperature effects, while radar and Wi-Fi disturbance methods detect chest wall motion [6, 7]. UWB Doppler radar monitors vital signs of a person in a noncontact manner by coherently processing the time history of narrow pulses emitted by the radar and subsequently backscattered by the human chest [8]. These sensor-attached methods are uncomfortable and inconvenient for mobile applications. They cannot be also applied to all patients. In addition, contactless devices monitoring the cardiopulmonary signal are more professional, which are not affordable for home-based use. Therefore, in this paper we prefer to detect cardiopulmonary signals without contact, especially for home, nonburden, ambient-assisted living and long-term monitoring.Magnetic induction is a valuable measurement contactless modality based on eddy current induction, which has been used to monitor human physiological information during the last decade. An alternating magnetic field is generated by excitation coils, which induces eddy currents into the conductive tissue of a body. Usually the diseased tissue of the organism will produce differences in conductivity, and this change could be reflected in the characteristics of the induced magnetic field, including frequency, amplitude, and phase. Then the pathological and physiological condition of biological tissue was analyzed by detecting the induced magnetic field. Previously, several approaches have been proposed to monitor human physiological signals with magnetic induction techniques. In [9], the authors proposed a simple measurement hardware holding with actual excitation and measurement coils. Testing experiments have been on a healthy adult volunteer who was resting on bed. However, the testbed is not equipped in the normal bed, and some manual postprocessing is required. In [10–12], they present a simulation model with a six-channel axial coil gradiometer underneath a neonatal thorax. They analyze the coverage of the whole bed, induction into the body, and distance between measuring points. Some animal trials have been done for evaluating the feasibility of monitoring vital signs by magnetic induction. They just discussed the feasibility by more simulating. Alternatively. In [13, 14], reported the related experiment of using magnetic induction methods to place sensors under the bed for physiological signal detection. The measuring principles are based on phase difference between exciting coil and receiving coil. Teichmann et al. in [15, 16] discussed the capacitive or inductive coupling for recording cardiorespiratory activity. The proposed system contains inductive sensor head and capacitive sensor head. The experiments reveal the inductive sensor is sensitive to tissue conductivity.Inspired by these prior works, we aim to improve and develop a home-used, feasible cardiorespiratory monitoring system based on magnetic induction. We consider the thoracic volume can be varied by physiological activity especially due to heart or lung breathing motion, and this will change the conductivity with the thoracic volume. In our solution, the only sensor-head coil is used to couple biological impedance change. The human body is connected in the oscillation tank circuit by coil inductive sensing thoracic volume variation. The oscillation frequency will change with thoracic equivalent impedance change caused by cardiopulmonary activity. The contribution of our work is as the following: we proposed a noncontact cardiorespiratory signal monitoring technique based magnetic induction, which exploits the intrinsic properties of detecting signal about respiratory or heart rate. And we enable the only sensor head in our designed system to obtain the useful information and simplify the conventional system hardware design. The proposed method is thoroughly evaluated by the theoretical analysis and the laboratory tests and also demonstrates comparable performance against guided breath behavior or the BCG method.


## 2. Magnetic Induction Measuring Model

### 2.1. Materials’ Conductivity VariationΔσ Related to the Impedance Variation ΔZ

Thez-component of magnetic field B from a coil of N turns, carrying AC current Iejωt centered on the x-y-axis, as seen in Figure 1, can be described by the Biot-Savart law:(1)Bzρ,z,t=μ0NIejωt4π∫02πρ-ρcos⁡θρ2+R2+z2-2ρRcos⁡θ3/2dθ,where μ0=4π×10-7 N/A2 is the permeability of free space and ρ is the radius from the z-axis j=-1, ρ=x2+y2, θ=tan-1(y/x)Figure 1
Single coil detection system coordinate.Thez-component of the induced magnetic field by the single turn coil can be denoted as(2)Bsz=C1∫L1L2dz∫0rzdρ∫0ρdρ′∫0πdθ∫0πdαF1ω,where(3)ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε,ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε=Δf+jΔX,F1=ρ′R2-ρ′Rcos⁡θ·ρ-hcos⁡αR2+ρ′2+z2-2ρ′Rcos⁡θ3/2·h2+ρ2+z2-2hρcos⁡α3/2.The magnetic fluxΦs of the circle area with the radius R and symmetry z-axis is(4)Φs=∫0RBsz2πρ′dρ′.The variation of induced electromotive force (EMF) in the coil sensor is(5)ΔE=-NdΦsdt=-jωNΦs.The induced EMF also can be denoted by(6)ΔE=-IejωtΔZ.Then(7)ΔZ=-ΔEIejωt=C2I^,where (8)C2=32π3×10-14N2f2R3σ+jωε,I^=∫L1/RL2/Rdz∫0rzdρ∫0ρdρ′∫01dh∫0πdθ∫0πdα,F=ρ′1-ρ′cos⁡θ·hρ-hcos⁡α1+ρ′2+z2-2ρ′cos⁡θ3/2·h2+ρ2+z2-2hρcos⁡α3/2.The relationship between the coil impedance with the material conductivity and dielectric constant is(9)ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε=ΔR+jΔX,where ω is current frequency and ε and σ are the material relative conductivity and dielectric parameter. For biological tissues, the real part of ΔZ is far greater than the imaginary part of ΔZ; then(10)ΔZ≈ΔR,ΔZ∝Δσ.


## 2.1. Materials’ Conductivity VariationΔσ Related to the Impedance Variation ΔZ

Thez-component of magnetic field B from a coil of N turns, carrying AC current Iejωt centered on the x-y-axis, as seen in Figure 1, can be described by the Biot-Savart law:(1)Bzρ,z,t=μ0NIejωt4π∫02πρ-ρcos⁡θρ2+R2+z2-2ρRcos⁡θ3/2dθ,where μ0=4π×10-7 N/A2 is the permeability of free space and ρ is the radius from the z-axis j=-1, ρ=x2+y2, θ=tan-1(y/x)Figure 1
Single coil detection system coordinate.Thez-component of the induced magnetic field by the single turn coil can be denoted as(2)Bsz=C1∫L1L2dz∫0rzdρ∫0ρdρ′∫0πdθ∫0πdαF1ω,where(3)ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε,ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε=Δf+jΔX,F1=ρ′R2-ρ′Rcos⁡θ·ρ-hcos⁡αR2+ρ′2+z2-2ρ′Rcos⁡θ3/2·h2+ρ2+z2-2hρcos⁡α3/2.The magnetic fluxΦs of the circle area with the radius R and symmetry z-axis is(4)Φs=∫0RBsz2πρ′dρ′.The variation of induced electromotive force (EMF) in the coil sensor is(5)ΔE=-NdΦsdt=-jωNΦs.The induced EMF also can be denoted by(6)ΔE=-IejωtΔZ.Then(7)ΔZ=-ΔEIejωt=C2I^,where (8)C2=32π3×10-14N2f2R3σ+jωε,I^=∫L1/RL2/Rdz∫0rzdρ∫0ρdρ′∫01dh∫0πdθ∫0πdα,F=ρ′1-ρ′cos⁡θ·hρ-hcos⁡α1+ρ′2+z2-2ρ′cos⁡θ3/2·h2+ρ2+z2-2hρcos⁡α3/2.The relationship between the coil impedance with the material conductivity and dielectric constant is(9)ΔZ=32π3×10-14N2f2R3I^Δσ+jωε=32π3×10-14N2f2R3I^Δσ+j64π3×10-14N2f2R3I^ωε=ΔR+jΔX,where ω is current frequency and ε and σ are the material relative conductivity and dielectric parameter. For biological tissues, the real part of ΔZ is far greater than the imaginary part of ΔZ; then(10)ΔZ≈ΔR,ΔZ∝Δσ.


## 3. Hardware Platform

In our design, the hardware platform consists of coil sensor head, signal receiver unit, and signal controlling unit. The coil sensor works for coupling the impedance change. The signal receiver unit includes LC oscillator circuit and postprocessing circuit, which respond for sampling the oscillating frequency related to the target conductivity variation. The signal controlling unit is built with the central control processor and display part used for controlling signal receivers, displaying the results and communicating with PC for all further signal processing steps.


### 3.1. Coil Sensor Head

In our system, we design the circular plane PCB coil for coupling the body conductivity change. The coil turn is 10, the average radius is 2.9 cm, and the width and distance of wire line both are 0.0254 cm, shown in Figure2(a). The theoretical inductance is 12.08 μH. The practical inductance value of our coil is 11.90 μH by high frequency Q meter (shown in Figure 2(b)).Figure 2
The coil sensor head and measuring device.
(a)
The designed circular plane coil sensor
(b)
The inductive value of designed coil byQ meter high frequency


### 3.2. Signal Receiver Unit

In our solution, the detecting signal unit works by detecting the change in coil sensor impedance caused by motion in the heart and lung. The oscillator frequency variation is related to the body conductivity change. The coil sensor is connected in Colpitts oscillating circuit. The schematic diagram is shown in Figure3, where L denotes the coil sensor inductance, ΔZ denotes the human body conductivity change, and R1, R2, R3, C1, and C2 are the circuit parameters. The oscillator frequency ωΔZ can be denoted by the equation(11)ωΔZ=C1+C2R1+ΔZ/R1+C3R1+R2+ΔZ/R1C1LC2+C3.Figure 3
Colpitts oscillating circuit considering inductive coilL and varying body impedance ΔZ.When the circuit is powered, the resonance frequency isω0 just caused by the coil sensor impedance L. Then the whole circuit is placed near the target object. The thoracic variation of target object due to heart and lungs motion will cause the body conductivity change; correspondingly the equivalent impedance of human body will be varied ΔZ. Then the oscillator frequency of the Colpitts oscillator circuit will be ω0+ωΔZ. For the value of varied impedance, ΔZ is very small and the signal of the frequency change ωΔz will be very weak. In order to detect the faint frequency change, the chip OPA690 was selected to be configured as postprocessing circuit for amplifying and filtering. OPA690 is unity-gain stable, voltage-feedback op amps. The amplifier circuit gain is 10, and the output is connected to the other OPA690, designed as Butterworth active filter with cutoff frequency 5.1 Mhz. For counting the frequency, the waveform transfer circuit is designed for shaping the frequency waveform.


### 3.3. Signal Controlling Unit

The signal controlling unit includes primary control processor and LCD. The main control board using chip MSP430F5438A is responsible for data receiving and communication with PC. The control chip clock frequency is set as 24 MHz. Thin Film Transistor (TFT) is used for displaying the detecting frequency signal waveform. The whole board has a dimension of 9.2 cm × 7.4 cm. The power supply is 3.3 V obtained by AMS1117. The whole designed system is shown in Figure4.Figure 4
Hardware platform: (a) main board and (b) TFTLCD.
(a)
(b)


## 3.1. Coil Sensor Head

In our system, we design the circular plane PCB coil for coupling the body conductivity change. The coil turn is 10, the average radius is 2.9 cm, and the width and distance of wire line both are 0.0254 cm, shown in Figure2(a). The theoretical inductance is 12.08 μH. The practical inductance value of our coil is 11.90 μH by high frequency Q meter (shown in Figure 2(b)).Figure 2
The coil sensor head and measuring device.
(a)
The designed circular plane coil sensor
(b)
The inductive value of designed coil byQ meter high frequency


## 3.2. Signal Receiver Unit

In our solution, the detecting signal unit works by detecting the change in coil sensor impedance caused by motion in the heart and lung. The oscillator frequency variation is related to the body conductivity change. The coil sensor is connected in Colpitts oscillating circuit. The schematic diagram is shown in Figure3, where L denotes the coil sensor inductance, ΔZ denotes the human body conductivity change, and R1, R2, R3, C1, and C2 are the circuit parameters. The oscillator frequency ωΔZ can be denoted by the equation(11)ωΔZ=C1+C2R1+ΔZ/R1+C3R1+R2+ΔZ/R1C1LC2+C3.Figure 3
Colpitts oscillating circuit considering inductive coilL and varying body impedance ΔZ.When the circuit is powered, the resonance frequency isω0 just caused by the coil sensor impedance L. Then the whole circuit is placed near the target object. The thoracic variation of target object due to heart and lungs motion will cause the body conductivity change; correspondingly the equivalent impedance of human body will be varied ΔZ. Then the oscillator frequency of the Colpitts oscillator circuit will be ω0+ωΔZ. For the value of varied impedance, ΔZ is very small and the signal of the frequency change ωΔz will be very weak. In order to detect the faint frequency change, the chip OPA690 was selected to be configured as postprocessing circuit for amplifying and filtering. OPA690 is unity-gain stable, voltage-feedback op amps. The amplifier circuit gain is 10, and the output is connected to the other OPA690, designed as Butterworth active filter with cutoff frequency 5.1 Mhz. For counting the frequency, the waveform transfer circuit is designed for shaping the frequency waveform.


## 3.3. Signal Controlling Unit

The signal controlling unit includes primary control processor and LCD. The main control board using chip MSP430F5438A is responsible for data receiving and communication with PC. The control chip clock frequency is set as 24 MHz. Thin Film Transistor (TFT) is used for displaying the detecting frequency signal waveform. The whole board has a dimension of 9.2 cm × 7.4 cm. The power supply is 3.3 V obtained by AMS1117. The whole designed system is shown in Figure4.Figure 4
Hardware platform: (a) main board and (b) TFTLCD.
(a)
(b)


## 4. Experiments and Results

### 4.1. Cardiopulmonary Activity and Impedance Changes

According the physiological data, human thorax is supposed as the sphere. Table1 lists the human thorax volume at various cardiopulmonary activity for male and female. During the simulation study, the distance between the coil and the thorax sphere is 10 cm, coil turn is 1, and the radius of coil is 6 cm. The oscillation frequency is 5 Mhz. The relationship between the changes of thoracic volume and the detection impedance was simulated, shown in Figure 5. During the process of deep inspiration exercise, the change of impedance was increased obviously. The main reason was that cardiopulmonary activity was obvious and the change of electrical conductivity was obvious in deep breathing process. And compared to women, men change a little more, which with their own physiological reasons have a certain relationship.Table 1
Human thorax volume at various cardiopulmonary activity for male and female.
Inspiration
Normal breath
Expiration
Holding breath
Volume change/L
1.25
0.25
0.75
0.75
Equivalent radius/cm
Male
8.85
7.33
5.37
6.94
Female
8.64
7.03
4.74
6.6Figure 5
The effect of the thoracic radius and impedance.


### 4.2. Coil Parameters

In order to measure the cardiopulmonary signal efficiently, the coil sensor parameters are very important such as the radius and measuring position. The saline was used for analyzing the coupling degree between the coil and the object signal. We simulated the testbed for discussing the coil sensor parameters. In the simulated testbed, coil sensor was placed below the cylindrical cylinder, which is used to contain 1 cm/ms saline. There was a small vessel, which is suspended in the middle of the cylindrical cylinder by manual controlling. The cylinder diameter and sample vessel diameter were 5 cm and 4 cm, respectively. The conductivity range in small vessel representing different types of tissue and organsσs= 1~8 cm/ms.Figure6(a) shows the relationship between the signal output Δf of the coil sensor by the object saline σs=8 cm/ms at the distance h. The results are based on different oscillation frequency and coil turns. At the same measuring condition, the higher the oscillation frequency and the greater the number of coil turns, the larger the couple frequency change output. There is a small oscillation in the tail of the curve. This may be caused by the sample vessel moving or manual operation factors (such as moving speed and hand shaking). Figure 6(b) shows the influence of coil sensor radius change for the couple frequency. The greater the coil sensor radius, the higher the simulated system’s sensitivity. But when the distance between the coil and target object is near 2 cm, for the target small vessel radius is 2 cm, the radius r=1.5 cm coil sensor has better sensitivity.Figure 6
The relation of signal output frequency with coil to object distanceh.
(a)
The sample conductivity is fixed; the fundamental frequency and the number of coil turns are varied
(b)
Only change in the coil radiusIn order to improve the system sensitivity, the optimal position of coil sensor has been discussed. The object body tank is represented by the cylinderσtunk=0.3 s/m, the heart is simulated by the sphere σheart=0.57 s/m, and the trapezoid cylinder represented the lungs σlung=0.35 s/m. There are six positions to be considered, shown in Figure 7(a). The sphere with radius 40 cm is the border, the cylinder represents the human body, the two cylindrical cylinders are denoted as lungs, and the small sphere is simulated as a heart. Positions 1 and 2 located in the lungs, position 3 located in the middle of heart, position 4 located in the boundary between the lung and heart, position 5 located in the middle of left lung and right lung, and position 6 is far from the lung, shown as Figure 7(a). Setting the heart as the center, we set 5 points at both sides in horizontal and vertical distribution. Each point is spacing 2 cm. Figure 7(b) shows the imaginary part of magnetic flux density when the coil sensors at different six positions move along the horizontal axis. The coil sensor is placed at various positions, and the variation in the value of the magnetic flux density is very small. The change is much higher when the coil sensor located at various horizontal location. At the distance to the center 6 cm, the magnetic flux density is largest. The value of position 3 is greatest due to heart conductivity, the measuring range of coil sensor moving is greater, and the signal contains much integrity. Then position 3 is selected as best position.Figure 7
Simulated coil sensor position at COMSOL and MATLAB. (a) Simulated coil sensor position at COMSOL. (b) Comparison of different coil positions in MATLAB.
(a)
(b)


### 4.3. System Evaluation

#### 4.3.1. Frequency Change Detecting

The proposed cardiopulmonary signal measuring board can be fixed on target object clothes, or sitting chair, or lying bed. The measured position should be near the human body, especially not too far as the discussed position 3. In our experiment, we evaluate the designed system at sitting state by the testbed. A young normal girl was asked to wear the test board on the backside (sternal left 2 cm), shown as Figure8(a). Regular daily activities were performed: standing and sitting. Figure 7 shows the result of one experiment in which the human subject follows instructions to perform various cardiorespiratory activities. In TFT displaying the x-axis is representing the time unit second, and the y-axis described the frequency unit MHz. At first, the subject was asked to perform normal breath 12.5 s and then deep inspiration 7.5 s when she was sitting on the chair. It clearly shows that the frequency change in normal breathing is about 200 Hz and in deep inspiration is 280 Hz~360 Hz, shown in Figure 8(b). We also observed the cardiorespiratory variation when subject was standing silently. He normally breathed 7 s, held breath for about 7 s, and recovered the normal cardiorespiratory status. The waveforms displayed that the frequency change in normal breath was 180 Hz~220 Hz and the frequency change in holding breath is 48 Hz~60 Hz, shown in Figure 8(c).Figure 8
Cardiopulmonary signal waveform in different states: (a) the subject with the board, (b) prone position, normal breathing, and deep breathing, and (c) sitting, normal breathing, and apnea.
(a)
(b)
(c)In the normal breathing, the changes in thoracic conductivity are greater than the changes in electrical conductivity caused by heart motion. Therefore, the frequency change in normal breath is mainly due to the respiratory activity, and the frequency change in holding breath is caused by heart motion. Compared with the experiment results in Figure6, frequency change at normal breath is about 280~360 Hz, and frequency variation due to holding breath is about 45~75 Hz. According to our system measuring rules, the amplitude change corresponding to the frequency change is about 5.09 mV~6.21 mV in normal breathing, 7.91 mV~10.18 mV at deep inspiration, and 1.27 mV~2.09 mV when the person is holding breath (Figure 8).


#### 4.3.2. Cardiopulmonary Signal Estimation

Another experiment was conducted for evaluating the system performance. We compared frequency variation related to heart and lung activity detected by our system with the other representing cardiopulmonary state physiology signal BCG, which also reflects the respiratory behavior variation [17]. We obtained the BCG data and frequency variation due to electromagnetic induction by our designed system when the target object is at normal breath state and holding breath state.Figures9 and 10 show 10-second recordings of the BCG signal and frequency change variation due to thoracic volume.Figure 9
Normal respiratory signal: (a) BCG signal and (b) system measurement signal.Figure 10
Holding breath signal: (a) BCG signal and (b) system measurement signal.In Figure9(a), the blue line is BCG signal, which can be seen as the ground truth, and red dashed line represented the respiratory waveform envelope. And in Figure 9(b), the blue line is frequency change due to cardiopulmonary activities by the designed system in this paper. Figure 10 shows the frequency variation and BCG signal when the object is holding breath. Obviously, the frequency change represents the heart and lung activity at two states.In order to determine whether to use the frequency change detected by magnetic induction for estimation of respiratory rate and breath rate, we also analyze the signal frequency by FFT. Welch periodogram of the 15-second recording signal was shown in Figures9 and 10. Normally, the respiratory frequency is about 12~20 bmp, and heart rate is about 60–100 bmp for normal adults. In Figure 11, there are two large peaks at 0.297 Hz and 1.12 Hz. The corresponding respiratory and heart rate are 18 bmp and 67 bmp. According to Figure 10, in the 8-second holding breath, the waveform changes 9 times due to heart rate. The heart rate of the object is 67 bmp. The results obtained by frequency analysis are basically the same.Figure 11
Welch period-gram (respiratory rate = 0.297 Hz and heart rate = 1.12 Hz).


## 4.1. Cardiopulmonary Activity and Impedance Changes

According the physiological data, human thorax is supposed as the sphere. Table1 lists the human thorax volume at various cardiopulmonary activity for male and female. During the simulation study, the distance between the coil and the thorax sphere is 10 cm, coil turn is 1, and the radius of coil is 6 cm. The oscillation frequency is 5 Mhz. The relationship between the changes of thoracic volume and the detection impedance was simulated, shown in Figure 5. During the process of deep inspiration exercise, the change of impedance was increased obviously. The main reason was that cardiopulmonary activity was obvious and the change of electrical conductivity was obvious in deep breathing process. And compared to women, men change a little more, which with their own physiological reasons have a certain relationship.Table 1
Human thorax volume at various cardiopulmonary activity for male and female.
Inspiration
Normal breath
Expiration
Holding breath
Volume change/L
1.25
0.25
0.75
0.75
Equivalent radius/cm
Male
8.85
7.33
5.37
6.94
Female
8.64
7.03
4.74
6.6Figure 5
The effect of the thoracic radius and impedance.


## 4.2. Coil Parameters

In order to measure the cardiopulmonary signal efficiently, the coil sensor parameters are very important such as the radius and measuring position. The saline was used for analyzing the coupling degree between the coil and the object signal. We simulated the testbed for discussing the coil sensor parameters. In the simulated testbed, coil sensor was placed below the cylindrical cylinder, which is used to contain 1 cm/ms saline. There was a small vessel, which is suspended in the middle of the cylindrical cylinder by manual controlling. The cylinder diameter and sample vessel diameter were 5 cm and 4 cm, respectively. The conductivity range in small vessel representing different types of tissue and organsσs= 1~8 cm/ms.Figure6(a) shows the relationship between the signal output Δf of the coil sensor by the object saline σs=8 cm/ms at the distance h. The results are based on different oscillation frequency and coil turns. At the same measuring condition, the higher the oscillation frequency and the greater the number of coil turns, the larger the couple frequency change output. There is a small oscillation in the tail of the curve. This may be caused by the sample vessel moving or manual operation factors (such as moving speed and hand shaking). Figure 6(b) shows the influence of coil sensor radius change for the couple frequency. The greater the coil sensor radius, the higher the simulated system’s sensitivity. But when the distance between the coil and target object is near 2 cm, for the target small vessel radius is 2 cm, the radius r=1.5 cm coil sensor has better sensitivity.Figure 6
The relation of signal output frequency with coil to object distanceh.
(a)
The sample conductivity is fixed; the fundamental frequency and the number of coil turns are varied
(b)
Only change in the coil radiusIn order to improve the system sensitivity, the optimal position of coil sensor has been discussed. The object body tank is represented by the cylinderσtunk=0.3 s/m, the heart is simulated by the sphere σheart=0.57 s/m, and the trapezoid cylinder represented the lungs σlung=0.35 s/m. There are six positions to be considered, shown in Figure 7(a). The sphere with radius 40 cm is the border, the cylinder represents the human body, the two cylindrical cylinders are denoted as lungs, and the small sphere is simulated as a heart. Positions 1 and 2 located in the lungs, position 3 located in the middle of heart, position 4 located in the boundary between the lung and heart, position 5 located in the middle of left lung and right lung, and position 6 is far from the lung, shown as Figure 7(a). Setting the heart as the center, we set 5 points at both sides in horizontal and vertical distribution. Each point is spacing 2 cm. Figure 7(b) shows the imaginary part of magnetic flux density when the coil sensors at different six positions move along the horizontal axis. The coil sensor is placed at various positions, and the variation in the value of the magnetic flux density is very small. The change is much higher when the coil sensor located at various horizontal location. At the distance to the center 6 cm, the magnetic flux density is largest. The value of position 3 is greatest due to heart conductivity, the measuring range of coil sensor moving is greater, and the signal contains much integrity. Then position 3 is selected as best position.Figure 7
Simulated coil sensor position at COMSOL and MATLAB. (a) Simulated coil sensor position at COMSOL. (b) Comparison of different coil positions in MATLAB.
(a)
(b)


## 4.3. System Evaluation

### 4.3.1. Frequency Change Detecting

The proposed cardiopulmonary signal measuring board can be fixed on target object clothes, or sitting chair, or lying bed. The measured position should be near the human body, especially not too far as the discussed position 3. In our experiment, we evaluate the designed system at sitting state by the testbed. A young normal girl was asked to wear the test board on the backside (sternal left 2 cm), shown as Figure8(a). Regular daily activities were performed: standing and sitting. Figure 7 shows the result of one experiment in which the human subject follows instructions to perform various cardiorespiratory activities. In TFT displaying the x-axis is representing the time unit second, and the y-axis described the frequency unit MHz. At first, the subject was asked to perform normal breath 12.5 s and then deep inspiration 7.5 s when she was sitting on the chair. It clearly shows that the frequency change in normal breathing is about 200 Hz and in deep inspiration is 280 Hz~360 Hz, shown in Figure 8(b). We also observed the cardiorespiratory variation when subject was standing silently. He normally breathed 7 s, held breath for about 7 s, and recovered the normal cardiorespiratory status. The waveforms displayed that the frequency change in normal breath was 180 Hz~220 Hz and the frequency change in holding breath is 48 Hz~60 Hz, shown in Figure 8(c).Figure 8
Cardiopulmonary signal waveform in different states: (a) the subject with the board, (b) prone position, normal breathing, and deep breathing, and (c) sitting, normal breathing, and apnea.
(a)
(b)
(c)In the normal breathing, the changes in thoracic conductivity are greater than the changes in electrical conductivity caused by heart motion. Therefore, the frequency change in normal breath is mainly due to the respiratory activity, and the frequency change in holding breath is caused by heart motion. Compared with the experiment results in Figure6, frequency change at normal breath is about 280~360 Hz, and frequency variation due to holding breath is about 45~75 Hz. According to our system measuring rules, the amplitude change corresponding to the frequency change is about 5.09 mV~6.21 mV in normal breathing, 7.91 mV~10.18 mV at deep inspiration, and 1.27 mV~2.09 mV when the person is holding breath (Figure 8).


### 4.3.2. Cardiopulmonary Signal Estimation

Another experiment was conducted for evaluating the system performance. We compared frequency variation related to heart and lung activity detected by our system with the other representing cardiopulmonary state physiology signal BCG, which also reflects the respiratory behavior variation [17]. We obtained the BCG data and frequency variation due to electromagnetic induction by our designed system when the target object is at normal breath state and holding breath state.Figures9 and 10 show 10-second recordings of the BCG signal and frequency change variation due to thoracic volume.Figure 9
Normal respiratory signal: (a) BCG signal and (b) system measurement signal.Figure 10
Holding breath signal: (a) BCG signal and (b) system measurement signal.In Figure9(a), the blue line is BCG signal, which can be seen as the ground truth, and red dashed line represented the respiratory waveform envelope. And in Figure 9(b), the blue line is frequency change due to cardiopulmonary activities by the designed system in this paper. Figure 10 shows the frequency variation and BCG signal when the object is holding breath. Obviously, the frequency change represents the heart and lung activity at two states.In order to determine whether to use the frequency change detected by magnetic induction for estimation of respiratory rate and breath rate, we also analyze the signal frequency by FFT. Welch periodogram of the 15-second recording signal was shown in Figures9 and 10. Normally, the respiratory frequency is about 12~20 bmp, and heart rate is about 60–100 bmp for normal adults. In Figure 11, there are two large peaks at 0.297 Hz and 1.12 Hz. The corresponding respiratory and heart rate are 18 bmp and 67 bmp. According to Figure 10, in the 8-second holding breath, the waveform changes 9 times due to heart rate. The heart rate of the object is 67 bmp. The results obtained by frequency analysis are basically the same.Figure 11
Welch period-gram (respiratory rate = 0.297 Hz and heart rate = 1.12 Hz).


## 4.3.1. Frequency Change Detecting

The proposed cardiopulmonary signal measuring board can be fixed on target object clothes, or sitting chair, or lying bed. The measured position should be near the human body, especially not too far as the discussed position 3. In our experiment, we evaluate the designed system at sitting state by the testbed. A young normal girl was asked to wear the test board on the backside (sternal left 2 cm), shown as Figure8(a). Regular daily activities were performed: standing and sitting. Figure 7 shows the result of one experiment in which the human subject follows instructions to perform various cardiorespiratory activities. In TFT displaying the x-axis is representing the time unit second, and the y-axis described the frequency unit MHz. At first, the subject was asked to perform normal breath 12.5 s and then deep inspiration 7.5 s when she was sitting on the chair. It clearly shows that the frequency change in normal breathing is about 200 Hz and in deep inspiration is 280 Hz~360 Hz, shown in Figure 8(b). We also observed the cardiorespiratory variation when subject was standing silently. He normally breathed 7 s, held breath for about 7 s, and recovered the normal cardiorespiratory status. The waveforms displayed that the frequency change in normal breath was 180 Hz~220 Hz and the frequency change in holding breath is 48 Hz~60 Hz, shown in Figure 8(c).Figure 8
Cardiopulmonary signal waveform in different states: (a) the subject with the board, (b) prone position, normal breathing, and deep breathing, and (c) sitting, normal breathing, and apnea.
(a)
(b)
(c)In the normal breathing, the changes in thoracic conductivity are greater than the changes in electrical conductivity caused by heart motion. Therefore, the frequency change in normal breath is mainly due to the respiratory activity, and the frequency change in holding breath is caused by heart motion. Compared with the experiment results in Figure6, frequency change at normal breath is about 280~360 Hz, and frequency variation due to holding breath is about 45~75 Hz. According to our system measuring rules, the amplitude change corresponding to the frequency change is about 5.09 mV~6.21 mV in normal breathing, 7.91 mV~10.18 mV at deep inspiration, and 1.27 mV~2.09 mV when the person is holding breath (Figure 8).


## 4.3.2. Cardiopulmonary Signal Estimation

Another experiment was conducted for evaluating the system performance. We compared frequency variation related to heart and lung activity detected by our system with the other representing cardiopulmonary state physiology signal BCG, which also reflects the respiratory behavior variation [17]. We obtained the BCG data and frequency variation due to electromagnetic induction by our designed system when the target object is at normal breath state and holding breath state.Figures9 and 10 show 10-second recordings of the BCG signal and frequency change variation due to thoracic volume.Figure 9
Normal respiratory signal: (a) BCG signal and (b) system measurement signal.Figure 10
Holding breath signal: (a) BCG signal and (b) system measurement signal.In Figure9(a), the blue line is BCG signal, which can be seen as the ground truth, and red dashed line represented the respiratory waveform envelope. And in Figure 9(b), the blue line is frequency change due to cardiopulmonary activities by the designed system in this paper. Figure 10 shows the frequency variation and BCG signal when the object is holding breath. Obviously, the frequency change represents the heart and lung activity at two states.In order to determine whether to use the frequency change detected by magnetic induction for estimation of respiratory rate and breath rate, we also analyze the signal frequency by FFT. Welch periodogram of the 15-second recording signal was shown in Figures9 and 10. Normally, the respiratory frequency is about 12~20 bmp, and heart rate is about 60–100 bmp for normal adults. In Figure 11, there are two large peaks at 0.297 Hz and 1.12 Hz. The corresponding respiratory and heart rate are 18 bmp and 67 bmp. According to Figure 10, in the 8-second holding breath, the waveform changes 9 times due to heart rate. The heart rate of the object is 67 bmp. The results obtained by frequency analysis are basically the same.Figure 11
Welch period-gram (respiratory rate = 0.297 Hz and heart rate = 1.12 Hz).


## 5. Conclusion

In this paper, we studied the problem of cardiopulmonary signal detection human localization by electromagnetic induction. We proposed a method that uses oscillator circuits to couple impedance change information due to breathing from coil sensors for human cardiopulmonary signal detection. The human impedance change due to thoracic volume variation can be coupled by coil sensors. Then oscillator frequency changed with the equivalent impedance of the couple coil. Our method has the advantage of reducing the obtrusiveness to the target object, while maintaining for detecting signals. We conducted experiments in a hardware platform and the feasibility and validity of the proposed method are evaluated. In the future, we will further improve our testbed and signal analysis method. For example, we need to minimize the board, and also we will also investigate the algorithm suitable for lying in a smart home environment.



---

*Source: 1752560-2017-06-14.xml*