Patent Application
20230389814
The present invention relates to a technique regarding the inference of a biological condition, and in particular, relates to a biological signal analysis device using an index relating to cardiac apex beat, a computer program, and a recording medium.
Cardiac apex beat reflects the movement of the apex of the left ventricle and is used in the evaluation of the function of the heart, mainly, the left ventricle. Normally, cardiac apex beat is objectively represented by an ACG (Apex cardiogram) recorded by a MCG (Mechanocardiogram) method, but the apex cardiogram is likely to be affected by recording skill because it captures the movement of the chest wall through various tissues between the heart and the chest wall. Further, mechanocardiogram recording devices are often relatively large and have restrictions that their measurement rooms have to be soundproof.
Considering this, Patent Document 1 discloses a bioinformation measurement device for apex cardiogram, including: an ultrasonic transmitting part causing an ultrasonic signal to be incident toward the vicinity of the cardiac apex from the body surface; and a receiving sensor that receives a reflected ultrasonic wave. The device is capable of easily creating an apex cardiogram even when it is not in a special environment such as a soundproof room.
The present inventors have proposed, in Patent Documents 2 to 5 and so on, arts for capturing vibration generated on the dorsal body surface of a person in a non-constraining manner and inferring the state of the person by analyzing the vibration. The vibration generated on the dorsal body surface of a person is vibration propagated from a human body inner part such as the heart and the aorta and contains information on atrial and ventricular systoles and diastoles, information on vascular wall elasticity which serves as an auxiliary pump for circulation, and information on reflected waves.
In Patent Document 2, slide calculation is performed in which a predetermined time width
is set in a time-series waveform of a dorsal body surface pulse wave of around 1 Hz extracted from vibration (biological signal) propagated through the body surface, to find a frequency slope time-series waveform, and from the tendency of its variation, for example, based on whether its amplitude is on the increase or on the decrease, a biological condition is inferred. It is also disclosed that power spectra of frequencies respectively corresponding to a function regulation signal, a fatigue reception signal, and an activity regulation signal that belong to a predetermined range from ULF band (ultra-low-frequency band) to VLF band (very-low-frequency band) are found through frequency analysis of the biological signals, and the state of a person is determined from time-series variations of the respective power spectra.
Patent Documents 3, 4 disclose a means for determining a homeostasis function level. Further, Patent Document 5 discloses a sound/vibration information collection mechanism including a resonance layer including a natural oscillator having a natural frequency corresponding to sound/vibration information of a biological signal.
Further, Non-patent Document 1 teaches that heart sound is composed of random vibrations. Non-patent Document 2 discloses that waveforms of the aforesaid apex cardiogram are classified and the result is used for diagnosis. Non-patent Document 3 gives an explanation regarding the fluctuation of cardiac cycle.
Patent Document 1: WO2012/165660
Non-patent Document 1: Nobuhiro, Y. et al. Development of Physical Condition Fluctuation Prediction Model Using Trunk Biosignals, Proceeding of Asia-Pacific Vibration Conference 2019 (2019)
Non-patent Document 2: Yoshikawa, J. et al. Physical Examination in Cardiology, Bunkodo Shoten Co., Ltd., 79-95 (2005)
Non-patent Document 3: Kobayashi M., Musha T.: 1/f fluctuation of heart beat period. IEEE Transactions on biomedical engineering. BME, 29, 45-457 (1982)
What Patent Document 1 discloses is a device using 192 ultrasonic oscillators, out of which twelve ultrasonic oscillators compose each set, and it scans while shifting the sets in sequence, to create an apex cardiogram. Further, for creating a more accurate apex cardiogram, it also discloses a device using a plurality of arrays of ultrasonic oscillators, the arrays each being composed of 192 ultrasonic oscillators. In Patent Document 1, it is difficult to obtain a desired apex cardiogram unless many ultrasonic oscillators are thus used. Therefore, a product unavoidably becomes very expensive and is not practical. Further, to create a clear apex cardiogram, care has to be taken to the detection of the rib position, a transmission frequency is limited, the positions of receiving elements are limited, and the device needs to satisfy these conditions, which also leads to a complicated structure and a higher manufacturing cost.
Further, Patent Document 1 states as a merit that the device requires only one person for the measurement while conventional mechanocardiogram recording devices require a plurality of people for handling However, the scanning with the ultrasonic oscillators requires the operation by a doctor or an engineer, and thus the measurement takes a lot of time and trouble, and only data during the scanning can be obtained. That is, the device is not suitable for capturing cardiac apex beat of a measurement subject over long hours.
On the other hand, the sensor that captures the vibration generated on the body surface in a non-constraining manner as is disclosed in Patent Documents 2 to 5 is composed of a three-dimensional knitted fabric, a film surrounding the three-dimensional knitted fabric, a microphone, and so on, and is capable of obtaining biological signals through the body surface only by being in contact with the person. That is, the device can obtain biological signal data only by being attached to the body of the person, without any operation by a doctor or the like. However, in Patent Documents 2 to 5, an autonomic nervous function and heart rate variability are mainly analyzed, and an analysis regarding cardiac apex beat which is the movement of the apex of the left ventricle is not done.
Further, the sensors used in Patent Documents 2 to 5 are suitable for passively capturing acoustic vibration of 20 to 100 Hz propagated from the dorsal body surface of a person but are not suitable for extracting cardiac apex beat which is in-vivo vibration of 10 Hz or lower.
Diagnosing using the apex cardiogram of Non-patent Document 3 requires professional knowledge for reading the waveform of the apex cardiogram, and the determination result may differ depending on each determiner, and thus there is still room for improvement in terms of determination accuracy and promptness. Therefore, the diagnosis using the apex cardiogram is only for auxiliary use. Moreover, such an apex cardiogram is used for diagnosis as information relating to cardiac apex beat, though auxiliary, but the inference of various health conditions (the presence/absence of a disease such as heart disease, the identification of a disease, whether the physical condition is good or not, and so on) of a person using an index relating to cardiac apex beat, other than the apex cardiogram, has not been done.
The present invention was made in consideration of the above, and has an object to provide a biological signal analysis device that can be manufactured easily and at a low cost, is capable of long-time measurement, is capable of extracting information regarding cardiac apex beat by analyzing biological signals collected by a biological signal detection sensor, and further is capable of inferring various health conditions of a person, and to provide a computer program and a recording medium.
To solve the aforesaid problems, the present inventors have developed a biological signal detection sensor that includes a three-dimensional knitted fabric, a microphone, and so on and is capable of obtaining biological signal data of 0.5 to 80 Hz. The biological signal detection sensor only passively captures acoustic vibration propagated from the body surface, and the captured biological signal data contains various biological sounds and in-vivo vibrations. As a result of widening a collectable frequency band to a low-frequency band, the biological signal data contains data regarding cardiac apex beat, but we faced a new problem that the data of the cardiac apex beat is hidden behind data of heart sound and thus the extraction of only the data of the cardiac apex beat is difficult. Under such circumstances, as a result of intensive studies for solving such a problem, the present inventors have completed the present invention.
Specifically, the biological signal analysis device of the present invention includes:
a frequency analyzing means that frequency-analyzes biological signal data obtained by a biological signal detection sensor through a body surface; and
a boundary frequency identifying means that finds a boundary frequency between vibration generated by cardiac apex beat and vibration generated by heart sound in the biological signal data, from a result of the frequency analysis of the biological signal data which result is obtained from the frequency analyzing means.
Preferably, the boundary frequency identifying means includes a means that finds, in the result of the frequency analysis, a power spectrum sudden changing point which is a boundary between harmonic vibration and random vibration, and identifies the boundary frequency based on the sudden changing point.
Preferably, the boundary frequency identifying means includes a means that finds the power spectrum sudden changing point, in consideration of a result of a frequency analysis of heart sound data measured simultaneously.
Preferably, the boundary frequency identifying means includes:
a log-log plot means that represents, in log-log axes, a waveform resulting from addition averaging of the results of the frequency analyses of the biological signal data and the heart sound data, using a log difference method; and
a sudden changing point identifying means that finds a fluctuation changing point from the waveform represented in log-log axes, and identifies the fluctuation changing point as the power spectrum sudden changing point.
Preferably, as the frequency analyzing means, short time Fourier transform is employed, and
the boundary frequency identifying means includes a means that finds the power spectrum sudden changing point from an analysis result of the short time Fourier transform.
Preferably, the frequency analyzing means includes a means that outputs the analysis result of the short time Fourier transform as image data showing time, frequency, and a degree of power spectrum variation, and
the boundary frequency identifying means includes a means that finds the power spectrum sudden changing point from the image data.
Preferably, the biological signal analysis device further includes:
a correlation data storage unit in which correlation data relating to the boundary frequency is stored; and
a measurement-time condition inferring means that infers a measurement-time health condition of a measurement subject, by referring to the correlation data.
Preferably, the correlation data is correlation data of the boundary frequency and heart rate, and
the measurement-time condition inferring means includes a means that infers a measurement-time boundary frequency of the measurement subject by collating a measurement-time heart rate of the measurement subject with the correlation data of the boundary frequency and the heart rate.
Preferably, the correlation data is correlation data of the boundary frequency and a fluctuation characteristic of heart rate variability, and
the measurement-time condition inferring means includes a means that infers a measurement-time fluctuation characteristic of heart rate variability of the measurement subject by collating a measurement-time boundary frequency of the measurement subject with the correlation data of the boundary frequency and the fluctuation characteristic of the heart rate variability.
Preferably, the biological signal analysis device further includes:
a cardiac apex beat waveform extracting means that filters the biological signal data, with an upper limit value being set to the boundary frequency identified by the boundary frequency identifying means, to find a waveform of the vibration generated by the cardiac apex beat; and
a waveform classifying means that classifies the waveform of the vibration generated by the cardiac apex beat, and
the measurement-time condition inferring means includes a means that causes the waveform classifying means to classify a measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, which waveform is obtained by the cardiac apex beat waveform extracting means, and infers the measurement-time health condition of the measurement subject based on data of a result of the classification of the measurement-time waveform.
Preferably, correlation data of a result of the waveform classification and health condition is stored in the correlation data storage unit in advance, and
the measurement-time condition inferring means includes a means that causes the waveform classifying means to classify the measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, which waveform is obtained by the cardiac apex beat waveform extracting means, and infers the measurement-time health condition of the measurement subject by collating the result of the classification of the measurement-time waveform with the pre-stored correlation data of the waveform classification result and the health condition.
Preferably, the biological signal analysis device includes a model creating means that uses information on waveforms and the health conditions as training data, the waveforms being output by the cardiac apex beat waveform extracting means to be classified by the waveform classifying means, to create, by machine learning, an inference model for inferring the health condition from the waveform information, and
the measurement-time condition inferring means includes a means that receives the information of the measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, and outputs a value of the measurement-time health condition based on the obtained information of the measurement-time waveform, using the inference model created by the model creating means.
Preferably, the waveform classifying means uses the following means (A) and (B):
(A) a means that classifies the waveform by a mathematical approach using Fourier series expansion; and
(B) a means that classifies the waveform based on a combination of one or two or more of the following data (1) to (4) obtained by engineering approaches:
(1) data indicating a shape of a time waveform in a predetermined time range of the waveform of the vibration generated by the cardiac apex beat;
(2) data of a graph of the time waveform, in which frequency and power spectrum are taken on a horizontal axis and a vertical axis, the data being obtained by the frequency analyzing means;
(3) image data in a predetermined time range that is the analysis result of the short time Fourier transform of the time waveform and shows the time, the frequency, and the degree of the power spectrum variation; and
(4) data regarding a type and a period of a waveform, in a predetermined time range, which data is extracted from a correlogram,
to classify waveforms of vibrations generated by the cardiac apex beat of healthy people into five, and
infers whether the measurement-time health condition of the measurement subject is a healthy condition or not based on whether or not the result of the classification of the measurement-time waveform of the measurement subject corresponds to any one of the five results of the waveform classification. Preferably, the biological signal detection sensor includes:
a three-dimensional knitted fabric;
a housing film that sealingly covers a periphery of the three-dimensional knitted fabric;
a microphone disposed on an outer side of the housing film;
a case covering the microphone; and
an external disturbance mixture preventing member that functions to prevent external disturbance from mixing into the microphone in the case.
Preferably, the external disturbance mixture preventing member is gel.
Further, the present invention provides a computer program causing a computer to process biological signal data that is obtained by a biological signal detection sensor through a body surface, to function as a biological signal analysis device, the computer program causing the computer to execute:
a procedure that frequency-analyzes the biological signal data; and
a procedure that identifies a boundary frequency between vibration generated by cardiac apex beat and vibration generated by heart sound in the biological signal data, from a result of the frequency analysis.
Preferably, the procedure that identifies the boundary frequency finds, in the result of the frequency analysis, a power spectrum sudden changing point which is a boundary between harmonic vibration and random vibration, and identifies the boundary frequency based on the sudden changing point.
Preferably, the procedure that identifies the boundary frequency finds the power spectrum sudden changing point, in consideration of a result of a frequency analysis of heart sound data measured simultaneously.
Preferably, the procedure that identifies the boundary frequency represents, in log-log axes, a waveform resulting from addition averaging of the results of the frequency analyses of the biological signal data and the heart sound data, using a log difference method, finds a fluctuation changing point from the waveform represented in log-log axes, and identifies the fluctuation changing point as the power spectrum sudden changing point.
Preferably, in the frequency analyzing procedure, short time Fourier transform is employed, and
the procedure that identifies the boundary frequency finds the power spectrum sudden changing point from an analysis result of the short time Fourier transform.
Preferably, the frequency analyzing procedure outputs the analysis result of the short time Fourier transform as image data showing time, frequency, and a degree of power spectrum variation, and
the procedure that identifies the boundary frequency finds the power spectrum sudden changing point from the image data.
Preferably, the computer is caused to execute a procedure that infers a measurement-time health condition of a measurement subject, by referring to correlation data that relates to the boundary frequency and is stored in a correlation data storage unit.
Preferably, the correlation data is correlation data of the boundary frequency and heart rate, and
the procedure that infers the measurement-time condition infers a measurement-time boundary frequency of the measurement subject by collating a measurement-time heart rate of the measurement subject with the correlation data of the boundary frequency and the heart rate.
Preferably, the correlation data is correlation data of the boundary frequency and a fluctuation characteristic of heart rate variability, and
the procedure that infers the measurement-time condition infers a measurement-time fluctuation characteristic of heart rate variability of the measurement subject by collating a measurement-time boundary frequency of the measurement subject with the correlation data of the boundary frequency and the fluctuation characteristic of the heart rate variability.
Preferably, the computer is further caused to execute:
a procedure that filters the biological signal data, with an upper limit value being set to the boundary frequency identified by the boundary frequency identifying means, to find a waveform of the vibration generated by the cardiac apex beat; and
a procedure that classifies the waveform of the vibration generated by the cardiac apex beat, and
the procedure that infers the measurement-time condition causes the procedure that classifies the waveform to classify a measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, and infers the measurement-time health condition of the measurement subject based on a result of the classification of the measurement-time waveform.
Preferably, correlation data of a result of the waveform classification and health condition is stored in the correlation data storage unit, and
the procedure that infers the measurement-time condition causes the procedure that classifies the waveform to classify the measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, and infers the measurement-time health condition of the measurement subject by collating the result of the classification of the measurement-time waveform with the correlation data of the waveform classification result and the health condition.
Preferably, information on waveforms and the health conditions are used as training data, the waveforms being output by the execution of the procedure that extracts the cardiac apex beat waveform to be classified by the procedure that classifies the waveform, and an inference model for inferring the health condition is created from the waveform information by machine learning, and
the procedure that infers the measurement-time health condition receives the information of the measurement-time waveform of the vibration generated by the cardiac apex beat of the measurement subject, and outputs a value of the measurement-time health condition based on the obtained information of the measurement-time waveform, using the inference model.
Preferably, the procedure that classifies the waveform uses the following procedures (A) and (B):
(A) a procedure that classifies the waveform by a mathematical approach using Fourier series expansion; and
(B) a procedure that classifies the waveform based on a combination of one or two or more of the following data (1) to (4) obtained by engineering approaches:
(1) data indicating a shape of a time waveform in a predetermined time range of the waveform of the vibration generated by the cardiac apex beat;
(2) data of a graph of the time waveform, in which frequency and power spectrum are taken on a horizontal axis and a vertical axis, the data being obtained by the frequency analyzing means;
(3) image data in a predetermined time range which is the analysis result of the short time Fourier transform of the time waveform and shows the time, the frequency, and the degree of the power spectrum variation; and
(4) data regarding a type and a period of a waveform, in a predetermined time range, which data is extracted from a correlogram,
to classify waveforms of vibrations generated by the cardiac apex beat of healthy people into five, and
infers whether the measurement-time health condition of the measurement subject is a healthy condition or not based on whether or not the result of the classification of the measurement-time waveform of the measurement subject corresponds to any one of the five results of the waveform classification.
Further, the present invention provides a recording medium in which the above-described computer program is recorded.
According to the present invention, it is possible to find a boundary frequency between cardiac apex beat and heart sound. The boundary frequency which is a new index regarding the cardiac apex beat correlates with the fluctuation characteristic of heart rate variability, and the use of the boundary frequency enables the inference of health conditions (the presence/absence of a disease such as a heart disease, the identification of a disease, whether or not the physical condition is good or not, and so on). Further, according to the present invention, since the boundary frequency can be represented by a function of heart rate, the use of correlation data based on this makes it possible to find the boundary frequency in real time.
Further, in the present invention, the waveforms of vibrations generated by cardiac apex beat found by using the aforesaid boundary frequency are classified, and the waveforms are associated with health conditions of people (the presence/absence of a disease such as a heart disease, the identification of a disease, whether or not the physical condition is good or poor, and so on), which enables the inference of a health condition. This makes it possible to automate the evaluation of a waveform of an apex cardiogram which has been done by a professional such as a doctor or the like. In particular, using the waveform information of the cardiac apex beat and the health conditions as training data to create, by machine learning, an inference model for the inference of a health condition from the waveform information makes it possible to improve inference accuracy of the health condition.
The present invention will be hereinafter described in more detail based on embodiments of the present invention illustrated in the drawings.
First, the structure of a biological signal detection sensor 1 used in this embodiment will be described based on
The three-dimensional knitted fabric 10 is formed of a pair of ground knitted fabrics disposed apart from each other and connecting yarns connecting the ground knitted fabrics. For example, the ground knitted fabrics each can be formed to have a flat knitted fabric structure (fine meshes) continuous both in a wale direction and a course direction using yarns of twisted fibers or to have a knitted fabric structure having honeycomb (hexagonal) meshes. The connecting yarns impart certain rigidity to the three-dimensional knitted fabric so that one of the ground knitted fabrics and the other ground knitted fabric are kept at a predetermined interval. Therefore, applying tension in the planar direction makes it possible to cause string vibration of the yarns of the facing ground knitted fabrics forming the three-dimensional knitted fabric or of the connecting yarns connecting the facing ground knitted fabrics. Accordingly, cardiovascular sound/vibration being a biological signal causes the string vibration and is propagated in the planar direction of the three-dimensional knitted fabric.
As the material of the yarns forming the ground knitted fabrics or of the connecting yarns of the three-dimensional knitted fabric, various materials are usable, and examples thereof include synthetic fibers and regenerated fibers such as polypropylene, polyester, polyamide, polyacrylonitrile, and rayon, and natural fibers such as wool, silk, and cotton. These materials each may be used alone or any combination of these may be used. Preferably, the material is a polyester-based fiber represented by polyethylene terephthalate (PET), polybutylene terephthalate (PBT), and the like, a polyamide-based fiber represented by nylon 6, nylon 66, and the like, a polyolefin-based fiber represented by polyethylene, polypropylene, and the like, or a combination of two kinds or more of these fibers. Further, the shapes of the ground yarns and the connecting yarns are not limited either, and round cross-section yarns, modified cross-section yarns, hollow yarns, or the like may be used. Carbon yarns, metallic yarns, and so on are also usable.
The following products are usable as the three-dimensional knitted fabric, for instance.
(a) product number: 49013D (manufactured by Sumie Textile Co., Ltd.), 10 mm thickness material:
front-side ground knitted fabric . . . twisted yarn of two polyethylene terephthalate fiber false twisted yarns with 450 decitexes/108 f
rear-side ground knitted fabric . . . twisted yarn of two polyethylene terephthalate fiber false twisted yarns with 450 decitexes/108 f
connecting yarn . . . polytrimethylene terephthalate monofilament with 350 decitexes/1 f
(b) product number: AKE70042 (manufactured by Asahi Kasei Corporation), 7 mm thickness
(c) product number: T28019C8G (manufactured by Asahi Kasei Corporation), 7 mm thickness
The three-dimensional knitted fabric 10 is covered with the housing film 20. In this embodiment, the housing film 20 is composed of two films 21, 22 made of a synthetic resin, they are arranged to cover the front surface and the rear surface of the three-dimensional knitted fabric 10, and peripheral edge portions of the films 21, 22 are bonded by welding. Consequently, the three-dimensional knitted fabric 10 is sealingly housed in the housing film 20. In bonding the peripheral edge portions of the films 21, 22, the peripheral edge portions are preferably bonded such that the three-dimensional knitted fabric 10 is slightly pressed in the thickness direction by the films 21, 22. This increases the tension of the three-dimensional knitted fabric 10 to more facilitate the generation of the string vibration of the yarns forming the three-dimensional knitted fabric 10.
The case 40 is attached to the outer side of the housing film 20, and the microphone 30 is disposed in the case 40. The gel 50 as an external disturbance mixture preventing member is filled around the microphone 30 in the case 40. The case 40 is made of a synthetic resin and has a function of preventing acoustic vibration propagated to the microphone 30 from spreading out, and the gel 50 inhibits external vibration from being captured by the microphone 30. A cord 30a that carries detected acoustic vibration data is connected to the microphone 30.
The biological signal detection sensor 1 is set in contact with any of various parts of a person, for example, the back, the chest, or the lumbar when used. Vibration of the body surface is propagated to the housing film 20 and the three-dimensional knitted fabric 10 to be captured by the microphone 30, but the biological signal detection sensor 1 may be attached to the surface of clothes when used, not limited to being directly pasted on the skin surface.
Here, the detection performance of the novel biological signal detection sensor 1 employed in this embodiment will be described as compared with the conventional signal detection sensor disclosed in Patent Document 2 and so on by the present inventors. Note that in the below, the biological signal detection sensor 1 of this embodiment will be called 4SR (Sound Sensing System using Stochastic Resonance), and the conventional biological signal detection sensor 1000 will be called 3SR (Sound Sensing System using Resonance).
(Evaluation of Sensing Performance of 4SR in Frequency Band of Heart Sound)
First, the basic principle of amplitude amplification of 4SR will be described using the photograph (a-4) of the cross section of the 3D net. In the 3D net used in this experiment, the upper and lower surfaces are base fabric layers formed of multifilament yarns, and piles formed of a monofilament yarn in the middle are knitted to the upper and lower base fabric layers to be friction-connected to the monofilament yarns. When the weight of a subject is applied to the upper surface of the 3D net, the piles knitted in an X-shape (hereinafter called X-piles) of the 3D net disposed in the air pack bend, and body movement causes the piles to rub against one another, so that sound and vibration are generated. The sound and vibration cause a stochastic resonance phenomenon. Further, tension is generated when the X-piles bend, so that a natural oscillator is formed. The natural frequency of the X-piles under tension is around 20 Hz. Therefore, owing to the two resonance mechanisms of the stochastic resonance and the string vibration, 4SR which is a sealed air pack is thought to amplify the vibration of 0.5 to 80 Hz.
Next, the performance of amplifying heart sound will be studied. As illustrated in
f0=0.161√{square root over (s0/m0 )} (1)
s0 represents the strength of a spring, and m0 represents the weight of the vibration plate.
The mechanically reproduced sound (hereinafter, called Reproduction PCG) output from the speaker passes through a 3D net placed on the upper surface of the accelerometer, and thereafter is measured by the microphones of 3SR and 4SR. The microphone 1010 of 3SR together with a 3D net 1020 is wrapped with an elastomer film 1030 and beads 1040. Damping is caused when the air in a slot flows to make the air around the microphone flow due to spring force of the internal 3D net. The magnitude of the damping is determined by a slit width in a peripheral portion of a cord 1010a attached to the microphone as illustrated in the photograph (a-7). Since 3SR has the microphone 1010 arranged in the elastomer film 1030, the cord 1010a has to be drawn out of the elastomer film 1030, and a gap is unavoidably formed in this portion. In 4SR, on the other hand, since the three-dimensional knitted fabric (3D net) 10 is structured such that the air pack 1A housed in the housing film 20 is separated from the gel pack 1B where the microphone 30 is disposed, the cord 30a of the microphone 30 is drawn out from the gel pack 1B, which does not obstruct the sealability of the air pack 1A. That is, in 4SR, the whole periphery of the three-dimensional knitted fabric (3D net) 10 is sealed with the elastomer film (housing film 20), and when a load is applied to 4SR, a uniform air pressure is applied to the elastomer film, so that tension is generated in the elastomer film and an environment where acoustic vibration is easily propagated is formed (see
(Evaluation of Sensing Performance of 4SR in Frequency Band of Cardiac Apex Beat)
In ACG of a healthy person, cardiac apex beat is composed of short-duration pulsation recognized only in early systole (a systolic wave is composed of a short-duration, sharp positive wave and a subsequent dicrotic limb, and an A-wave originating in left atrial contraction is recognized as a trace) and is described as tapping, and its good recording in ACG is difficult. Therefore, a subject (heart rate: 62/min) who was in a sitting posture so that cardiac apex beat described as tapping could be sensed by palpation was selected, and F-APW (biological information measured from the anterior chest surface: Front Acoustic Pulse Wave) was measured using 3SR and 4SR, was recorded in a data logger, and filtered, and the following analysis was conducted.
The input waveform recorded in the data logger is frequency-analyzed, a ratio of the power spectra of 4SR and 3SR (PSD4SR/PSD3SR) at each frequency is defined as a gain, and sensing performance is evaluated based on the gain. On the anterior chest surface at a position corresponding to the left fifth intercostal space, 4SR is placed 6 cm left of the chest midline, and 3SR is placed 10 cm left thereof (corresponding to the vicinity of V4 induction of an electrocardiogram). Since the installation positions of the microphones for PCG and the 4SR sensing systems are places where vesicular breathing sound is likely to mix, F-APW is measured for ten seconds while the breath is held. Note that the procedure of the measurement experiment was to continue breathing at rest for five minutes, and hold the breath after two-second inspiration. As a subject, a healthy male in his forties without any underlying disease in the cardiocirculatory system was selected.
(APW Measurement)
Next, APW is measured. All the measurements were conducted in a sitting posture. 4SR for capturing F-APW considered as including cardiac apex beat is placed on the anterior chest surface at a position corresponding to the fifth left intercostal space on the left midclavicular line (near V4) and 10 cm left of the chest midline, and 4SR for R-APW (biological information measured from the posterior chest: Rear Acoustic Pulse Wave) is placed on the posterior chest surface at a position that is the same height as that of 4SR for F-APW which is on the anterior chest surface, and is 6 cm left of the dorsal midline. 4SR for capturing L-APW of the lumbar (biological signal measured from the lumbar: Lumbar Acoustic Pulse Wave) is placed on the third to fourth lumbar vertebrae median which is directly behind the umbilical region. The microphone for PCG is placed at a position corresponding to the cardiac apex. Further, in ECG (Electrocardiogram), II induction is obtained.
The position of the cardiac apex microphone sensor on the anterior chest is an optimal position for recording first sound, third sound, and fourth sound, but at this position, it is difficult to record second sound well. It is said that the second sound is generated by the hyperextension of the semilunar valve caused by arterial blood trying to flow back to the ventricle in late systole and subsequent reaction of the arterial wall (cardiohemic system), and the second sound becomes weak at the cardiac apex.
Since the installation positions of the microphones for PCG and the 4SR sensing system are places where vesicular breathing sound is likely to mix, the result of the ten-second measurement during which the breath is held is used for the calculation. Note that the APW measurement experiments at rest are conducted under the following four conditions: a) a state where natural breathing is continued for five minutes; b) a state where the breath is held for thirty seconds after two-second inspiration; c) a state of thirty-second natural breathing; and d) a state where the breath is held for thirty seconds after two-second inspiration.
Subjects were 70 healthy males and females in a group with a relatively thin chest wall and a group with a relatively thick chest wall and thus have different features.
As for fifty healthy people having a relatively thin chest wall, their ages were distributed in twenties to sixties, and a ratio of the males was 78%. The thin chest wall facilitates the observation of a sine wave or a triangular wave which is the feature of cardiac apex beat. Using these subjects, later-described five waveform types were defined, and BF (Boundary Frequency) between vibration ascribable to cardiac apex beat and vibration ascribable to heart sound was decided.
Using twenty subjects in their twenties who had a relatively thick chest wall and were muscular, the compatibility of the analysis method was evaluated and its extensibility was studied, based on the previously decided five waveform types.
Clinical profiles of the fifty healthy subjects with a relatively thin chest wall are as follows:
average age 38.2 (±10.3) years old, 39 males, 11 females, average BMI value 23 (±3.4), out of them, four have hypertension and ten are smokers.
The clinical profiles of the twenty subjects with a relatively thick chest wall are as follows: average age 23.4 (±3.6) years old, 20 males, average BMI value 20.8 (±2.4), out of them, one has hypertension and five are smokers.
As data to be analyzed, used was data during ten seconds from a part where RR variation was within 15%, in the first thirty-second period during which the breath was held, and as heart rate, an average value of data in this ten-second period was used.
(Performance Evaluation Result and APW Measurement Result of 4SR Sensing System)
In
Next, the damping performance of a mechanical vibration system constituting the sensing system will be studied.
W=πcωα2 (2)
[C is damping coefficient, ω is angular frequency, and α is amplitude]
Comparison between the dynamic characteristics of 4SR and 3SR shows that 3SR has a substantially equal spring characteristic to that of 4SR, and in addition, is noticeably influenced by damping due to the air that flows in/out.
Further, amplitude magnification Z/Y used in the evaluation of vibration displacement is given by the following, where Z is relative displacement and Y is displacement of cardiac apex beat.
Here, the natural angular frequencies of the mechanical vibration systems of the sensors 4SR and 3SR are represented by the following,
and, a damping ratio is represented by the following.
The following values were substituted.
4SR=485N/m, =24 kg , ζ4SR=0.20,
3SR=9367N/m,
=24 kg, ζ3SR=0.53
It is seen that, since the natural frequency of the sensing system of 4SR is set low, 4SR is capable of the measurement of a frequency band of 0.45 Hz or higher in the case of the measurement in the range of Z/Y≈1, and 3SR is capable of the measurement of a frequency band of 2.10 Hz or higher. Note that, in the case of the damping ratio=0.7, 4SR is capable of the measurement of Hz or higher, and 3SR is capable of the measurement of 3.15 Hz or higher. By thus changing the damping ratio, it is possible to adjust the characteristic of a mechanical filter.
Cardiac auscultation is to capture the characteristic of heart sound in a high frequency range of around 50 to 100 Hz and interpret and confirm it, and it is important to record acoustic vibration of high frequencies. On the other hand, cardiac apex beat is short-duration pulsation that is within a frequency range lower than that of heart sound and is recognized only in early systole, and it is difficult to record a good apex cardiogram if emphasis is put on high-frequency measurement. Since the required specifications of the sensing system thus differ depending on which of heart sound and cardiac apex beat is to be measured, quantitative evaluations of the sensing systems were conducted separately for these.
Incidentally, the air damping of 3SR had a small influence on random vibration related to heart sound, but as illustrated in
This agrees with the simulation results in
A vibration transmitting mechanism will be discussed. Since tension is generated also in the elastomer film (the housing film denoted by reference sign 20 in 4SR and the elastomer film denoted by reference sign 1030 in 3SR) itself, low-frequency vibration of around 1 Hz is transmitted to the internal 3D net in 4SR through tension variation. As a result, vibration of 13 Hz or higher was measured by 3SR and 4SR as having substantially equal power spectra between 13 to 30 Hz illustrated in
Conventionally, the frequency to be analyzed by 3SR was set to 10 to 30 Hz, but as a result of the above experiment result, the vibration frequency of F-APW to be measured by 4SR was set to 0.5 to 80Hz in consideration of a performance ensured range of the microphone, which is 0.1 Hz or higher, and the frequency to be analyzed of heart sound was set to 25 to 45 Hz or 40 to 80 Hz instead of 100 Hz or higher because most of the heart sound gather in a low audible spectrum. Further, the sampling frequency was set to 1000 Hz in 4SR though set to 200 Hz in 3SR.
Here, a comparative experiment was conducted in which acoustic vibration was measured using the biological signal detection sensor 1 (4SR) of this embodiment and a structure having only the air pack 1A with the gel pack 1B removed from 4SR (4SR without gel pack).
As is apparent from these charts, 4SR has higher detection sensitivity to acoustic vibration information than the structure only with the air pack 1A (4SR without gel pack). In particular, when the sound source was a person, 4SR was capable of clearly capturing the acoustic vibration information because an amplitude change of about ±0.2 [V] occurred, while it was difficult for the structure only with the air pack 1A to capture an amplitude change. Further, as illustrated in
Further, the microphone 30 is arranged on a disk-shaped support member (member denoted by reference sign 30b in
Further, a time waveform of MAV (micro acoustic vibration) output from the trunk and measured by a microphone attached directly to the chest and a time waveform of F-APW (acoustic vibration information influenced by the response of the mechanical vibration system of the aforesaid members composing 4SR and further influenced by the resonance of the natural oscillator and the stochastic resonance) measured from 4SR were frequency-analyzed, and a gain of F-APW/MAV was found through the calculation of their amplitude magnifications. The gain of F-APW/MAV was 42 dB in the frequency band of the fundamental harmonic of cardiac apex beat, was slightly more than dB in a higher harmonic frequency band, and was 20 dB in the lowest frequency band of heart sound of around 20 to 30 Hz, which is a frequency band of the natural oscillator. The effects of the response of the mechanical vibration system and the stochastic resonance appeared as over 40 dB at the maximum and as 10 dB at the minimum. An amplitude difference was slightly less than 200 times. In F-APW, a waveform in a 1 to 1.5 Hz band was observed, and information on heart sound and cardiac apex beat was contained. A time lag of MAV and F-APW was 0.12 seconds.
Biological Signal Analysis Device
Next, a biological signal analysis device 100 having a computer function, in which a computer program for processing data obtained from the biological signal detection sensor 1 of this embodiment, is set will be described based on
The biological signal analysis device 100 processes time-series data of biological signals obtained by the biological signal detection sensor 1 to obtain a waveform of cardiac apex beat and further infers biological conditions, that is, various health conditions (the presence/absence of a disease such as a heart disease, the identification of a disease, whether a physical condition is good or not, and so on) of a person. The biological signal analysis device 100 is constituted by a computer (including a personal computer, a microcomputer incorporated in a device, and so on) and receives the time-series data of the biological signals transmitted from the microphone 30 of the biological signal detection sensor 1. It includes: a frequency analyzing means 110 that performs predetermined processing using the received time-series data; and a boundary frequency identifying means 120.
As a storage unit, it has a correlation data storage unit 150 in which correlation data of a boundary frequency between vibration ascribable to cardiac apex beat and vibration ascribable to heart sound with heart rate is stored, the correlation data being created through advance measurement In more detail, in the biological signal analysis device 100, a computer program causing the execution of procedures functioning as the frequency analyzing means 110 and the boundary frequency identifying means 120 is stored in the storage unit (including not only a recording medium such as an internal hard disk as the computer (biological signal analysis device 100) but also various types of removable recording media and a recording medium of another computer connected through communication means). Note that the biological signal analysis device 100 can also be actualized using an electronic circuit having one recording circuit or more in which the computer program implementing the frequency analyzing means 110 and the boundary frequency identifying means 120 is incorporated.
The computer program can be stored in a recording medium to be provided. The recording medium storing the computer program may be a non-transitory recording medium. The non-transitory recording medium is not limited, and examples thereof include recording media such as a flexible disk, a hard disk, CD-ROM, MO (magneto-optical disk), DVD-ROM, and a memory card. It is also possible to transfer the computer program to the computer through a communication line to install it therein.
Here, since the biological signal detection sensor 1 is capable of capturing heart sound as described above, the frequency analyzing means 110, upon receiving the obtained biological signal data, frequency-analyzes biological signal data obtained through the body surface by the biological signal detection sensor 1.
The boundary frequency identifying means 120 finds BF (Boundary Frequency) between the vibration ascribable to cardiac apex beat and the vibration ascribable to heart sound in the biological signal data, from the frequency analysis result of the biological signal data obtained from the frequency analyzing means 110. The boundary frequency identifying means 120 includes a means that finds a power spectrum sudden changing point which is the boundary between harmonic vibration and random vibration, in the frequency analysis result obtained by the frequency analyzing means 110 and identifies the boundary frequency based on the sudden changing point.
According to this embodiment, finding the boundary frequency makes it possible to distinguish the waveform of cardiac apex beat from the waveform of heart sound to obtain the waveform of cardiac apex beat from the biological signal.
Further, the biological signal analysis device 100 of this embodiment uses the boundary frequency itself as a determination criterion. Specifically, as is apparent from a later-described experiment, the boundary frequency is represented by a quadratic function with heart rate, and from correlation data between these, the boundary frequency can be found based on the heart rate. From the correlation between the boundary frequency and the heart rate, it is seen that the heart rate of the heart in a highly efficient condition is about 70 to 80/min. This is a state where an exercise load of the heart is small, and the boundary frequency in this case is smallest. Therefore, from a measurement-time boundary frequency obtained by measuring the heart rate and using the correlation data, it is possible to infer the state of the heart at this instant.
Preferably, the biological signal analysis device 100 of this embodiment has: the correlation data storage unit 150, indicated by the imaginary line in
The measurement-time condition inferring means 130 requires information on the heart rate when inferring the condition. As for the heart rate information, a heart rate calculating means 140 indicated by the imaginary line in
The correlation data stored in the correlation data storage unit 150 is correlation data of the boundary frequency between vibration ascribable to cardiac apex beat and the vibration ascribable to heart sound with the heart rate, which are measured in advance. This correlation data is configured by a quadratic function indicating the correlation between the boundary frequency and the heart rate obtained when the boundary frequency is found. Therefore, collating the measured heart rate with the correlation data stored in the storage unit of the correlation data storage unit 150 makes it possible to find the boundary frequency.
The boundary frequencies stored in the correlation data storage unit 150 are found through an analysis of a lot of biological signal data in advance. The boundary frequency is found as the aforesaid sudden changing point identified by the boundary frequency identifying means 120.
The sudden changing point can be found as follows. First, the biological signal data obtained through the body surface by the biological signal detection sensor 1 is frequency-analyzed, and in consideration of this frequency analysis result and the frequency analysis result of heart sound data measured simultaneously, the boundary frequency between the vibration ascribable to cardiac apex beat and the vibration ascribable to heart sound, in the biological signal data, is extracted.
Specifically, the frequency analysis results of the biological signal data found by the frequency analyzing means 110 are subjected to addition averaging, the obtained waveform is represented in log-log axes using frequency and power spectrum, and the frequency analysis results of the heart sound data are subjected to addition averaging, the obtained waveform is represented in log-log axes using frequency and power spectrum, a fluctuation changing point is found from a waveform of a log difference of the two waveforms represented in log-log axes, this fluctuation changing point is considered as a practical disappearing point where higher harmonic components of the cardiac apex beat become very small, and a frequency at the disappearing point is considered as corresponding to the aforesaid sudden changing point, whereby the boundary frequency can be found.
As the frequency analyzing means 110, a STFT (Short Time Fourier Transform) method can also be employed. In this case, the boundary frequency identifying means 120 finds the aforesaid power spectrum sudden changing point from the analysis result of the short time Fourier transform. Preferably, the analysis result of the short time Fourier transform is output as image data showing time, frequency, and a power spectrum variation degree. This enables the boundary frequency identifying means 120 to find the power spectrum sudden changing point from the image data of the short time Fourier transform. A specific means for identifying the sudden changing point, that is, a means for finding the boundary frequency will be described later.
Here, the biological signal detection sensor 1 is attached to the back, the chest, the lumbar, or the like of the body of a person to capture sound and vibration in the body. The obtained biological signal is a collection of various biological sounds and in-vivo vibrations. On the other hand, cardiac apex beat is vibration hidden behind a heart sound waveform, and it is difficult to separate them. In this embodiment, however, the aforesaid BF (Boundary Frequency) for distinguishing the cardiac apex beat waveform and the heart sound waveform is found, which enables the inference of the biological condition relating to the cardiac apex beat based on the heart rate.
Hereinafter, methods of the frequency analysis by the frequency analyzing means 110 and the identification of the boundary frequency by the boundary frequency identifying means 120 will be described in detail. As described above, this embodiment uses two methods: the method of using the waveform represented in log-log axes to find, in the waveform, the fluctuation changing point based on the frequency at which the harmonic vibration of higher harmonics practically disappears, and considering this fluctuation changing point as the power spectrum sudden changing point; and the method of finding the power spectrum changing point from the analysis result of the short time Fourier transform.
[Method of Extracting Boundary Frequency (BF)]
(BF extraction method based on disappearing frequency of harmonic vibration of higher harmonics of APW)
Heart sound is random vibration as described in Non-patent Document 1, while the later-described experiment result has made it obvious that a cardiac apex beat wave is mainly composed of harmonic vibration composed of a fundamental harmonic and higher harmonics. Therefore, we have focused on frequency-analyzing the waveform of APW obtained from the 4SR sensing system attached to the anterior chest, performing the addition averaging, and finding the power spectrum changing points of the harmonic vibration of the higher harmonics and the random vibration, that is, the power spectrum sudden changing point. The higher harmonic components of CAB (Cardiac Apex Beat) of the cardiac apex beat wave become smaller in power spectrum as the frequency becomes higher. On the other hand, the power spectrum of the random vibration of CAS (Cardiac Acoustic Sound) of heart sound changes in magnitude without depending on frequency. Here, a frequency at which the power spectrum of the higher harmonic components of CAB becomes small and a variation behavior of CAS changes is called boundary frequency, which will be hereinafter called BF. Note that the occurrence of heart rate variability or blood pressure variability leads to instantaneous variation of the waveform of APW, and one of the possible reasons of the instantaneous variation of the APW waveform is variation of the dominant frequency of PCG.
Because of individual differences of subjects, impedance regarding the transmission characteristic of acoustic vibration of the trunk also differs, which may make it difficult to find BF. Therefore, the frequency at which the higher harmonic components of the cardiac apex beat wave disappear (the frequency at which the power spectrum of the higher harmonic components becomes very small and is practically on a negligible level) is found using the log difference method. In the addition averaging of APW, Fourier transform is performed with a time window of 8.2 seconds and 90% overlapping. The waveform resulting from the addition averaging is represented in log-log axes. This is because the addition theorem is usable for the power spectrum waveform represented in log-log axes. Note that the PCG waveform is used for heart sound. When the log difference method is used, an absolute value of the PCG power spectrum in a frequency band of 0.5 to 20 Hz becomes large and the power spectrum waveform can be reversed in phase. Consequently, the power spectra of the higher harmonic components of the cardiac apex beat and the random vibration of the heart sound are decreased to facilitate finding BE
(BF Extraction Method Using Short Time Fourier Transform)
The power spectrum STFT (Short Time Fourier Transform) is applied to F-APW, PCG, APW×PCG, and APW×PCG−1 to visualize a fluctuation changing point. A formula of the short time Fourier transform used in this analysis is as follows.
[Formula 4]
STFT(τ, ω)=∫−∞∞x(t)W(t−τ)e−jωtdt (4)
x(t) is measured signal data, and W(t — r) is a window function. To obtain the frequency resolution of 0.5 Hz or lower, the number of points of the window function is set to 212=4096. Under the sampling frequency of 1000 Hz, the 4096 points correspond to 4.096 seconds, and 4.096 seconds correspond to 0.244 Hz. As a result, the frequency resolution can be set equal to or less than Hz, specifically, set to 0.244 Hz. Then, a shift width of the window function was set to 50 points. t represents measurement time, and r represents time-series data based on the moving time of the window function. Multiplying the measured signal data by the window function makes it possible to create frequency information that changes with time. Note that the frequency information is generated every 50/1000=0.05 seconds.
Next, using BF identified based on the aforesaid two BF extraction methods, a cardiac apex beat waveform (Front CAB) and heart sound (Front CAS) are drawn from F-APW, and similarly Rear CAB, Rear CAS, Lumbar CAB, and Lumbar CAS are drawn from R-APW and L-APW. Further, since it has been found out in a preliminary experiment and a preliminary study that HR correlates with BF, a correlation chart of BF and HR is created using data of fifty subjects for the purpose of studying the correlation between BF and HR.
The BF extraction method based on the frequency at which the higher harmonic components disappear is an approach to decide BF from the frequency, while the BF extraction method using the short time Fourier transform is an approach to decide BF based on the power spectrum threshold. Specifically, the former is a method to decide, as BF, the frequency at which a decreasing change of the power spectrum height shifts to an increasing change as the frequency is changed, and the latter is a method to set 70% of an average value of a power spectrum group present around BF as the threshold to decide BF.
Here,
In the case of the subject with HR: 58/min, a systolic wave measured from Front CAB presents a late systolic bulge of ACG, but from ECG and PCG (Front CAS), it is less likely that the subject has dilated cardiomyopathy. It is thought that the cardiac apex beat is influenced by gravitational force to resemble a sustained pattern.
In the case of the subject with HR: 71/min, in ECG, a P-wave is not clear and is too small to be seen. In Front CAB, an A-wave originating in left atrial contraction (atrial bulging wave) as is confirmed in ACG is confirmed, and in Front CAS as well, a low-frequency wave is present earlier than first sound of PCG, which does not correspond to ectopic rhythm.
In the case of the subject with HR: 79 min, a waveform measured from Front CAB presents a sustained pattern, but an A-wave, a RFW (Rapid Filling Wave) waveform, and Front CAS indicate that this does not correspond to a case of left ventricular hypertrophy and leakage from the left ventricle in systole (mitral valve regurgitation or the like).
In the case of the subject with HR: 93/min, ECG presents a low voltage and presents a high heart rate in a resting state, but Front CAB having low-frequency components and high-frequency components and Front CAS in which first sound does not rise nor weaken rule out the possibility of an underlying heart disease.
(Correlation of Heart Rate Variability Fluctuation Characteristic with Heart Rate and BF)
A power spectrum slope on a log-log plot is called a fluctuation coefficient here, and a steep changing point of the fluctuation coefficient will be hereinafter called a fluctuation changing point. At the fluctuation changing point, complexity is lost and a breakpoint occurs. The fluctuation coefficient is represented by power spectrum P(f) and frequency (f), which are log P(f) and log(f) on the log-log plot, and in the case where log P(f) and log(f) have a linear relationship, the following formula holds.
[Formula 5]
log P(f)=κ−n log(f)=κ+log(1/fn) (5)
[κ is Constant and n is Fluctuation Coefficient]
To evaluate the fluctuation characteristic of heart rate variability found from ECG, experimental data measured from an experimental protocol of APW measurement is used. The fluctuation characteristic is calculated using Formula (5). An obtained value of n is a minus numerical value (about −0.5 to −1.5) but is expressed by an absolute value because Non-patent Document 3 defines it as 1/f fluctuation in a cardiac cycle.
To evaluate the fluctuation of the heart rate variability, data during no breathing, natural breathing, and no breathing each having thirty-second duration, that is, data during totally ninety seconds are used. As the heart rate here, an average value during the ninety seconds is used. [Results of BF Extraction by the Extraction Methods]
(Result of BF extraction based on disappearing frequency of harmonic vibration of higher harmonics of F-APW)
(1) F-APW is found from 4SR data without filtering processing and is frequency-analyzed. Hereafter, the result is represented by F-APW.
(2) F-APW is subjected to addition averaging (8.2-second time window with 90% overlapping).
(3) A waveform is extracted from PCG and is frequency-analyzed. Hereinafter, the result will be represented by PCG.
(4) The same addition averaging (8.2-second time window with 90% overlapping) as that in (2) is performed.
(5) F-APW×PCG−1 is created using the power spectra F-APW and PCG of F-APW and PCG found in (2) and (4), and a breakpoint is found.
(6) From the breakpoint on a waveform of F-APW×PCG−1, the disappearing frequency of higher harmonics of a cardiac apex beat wave having fluctuation and an appearing frequency of a random vibration component of heart sound which becomes white noise are identified.
(7) A BF line is drawn from the disappearing point of the higher harmonics on F-APW×PCG−1, and an intersection of the BF line and F-APW is identified as BF.
(8) 0.5 Hz to BF is a CAB band (hereinafter, CAB (0.5-BF)), and BF to 50 Hz is a CAS band (hereinafter, CAS (BF-50)).
(9) A 0.5 Hz to BF bandpass filter is applied to F-APW and it is confirmed that the higher harmonic components of CAB decrease near BE
(10) A BF to 50 Hz bandpass filter is applied to F-APW and it is confirmed that the power spectrum changes near BF without depending on frequency.
Next,
The study of the obtained results shows that Front CAB (5-BF) is formed of the combination of a plurality of 5 to 10 Hz waveforms. On the other hand, Rear CAB (0.5-BF) and Lumbar CAB (0.5-BF) had harmonic vibration waveforms mainly composed of a waveform with a higher frequency than that of the waveform illustrated in Front CAB (0.5-BF), and had waveforms different from that of Front CAB (0.5-BF). Further, heart sound was random vibration formed of a high-frequency waveform in which multiple vibration components and multiple vibration waveforms are combined. (Result of BF Extraction Using Short Time Fourier Transform)
In the STFT charts, the vertical axis represents frequency and the horizontal axis represents the time when STFT is output. The color intensity indicates a degree of power spectrum variation, and by changing a range of the color intensity, it is also possible to express the degree of the fluctuation coefficient represented by Formula (5). The degree of the power spectrum variation is indicated by the intensities of Red, Green, and Blue (RGB). In a place where the power spectrum variation is small and a part where the fluctuation coefficient does not suddenly change, and in random vibration with a small amplitude, planar Red, Green, and Blue appear. Conversely, in a place where the power spectrum variation is large and a part where the fluctuation coefficient suddenly changes, Red, Green, and Blue appear in a line shape.
The aforesaid log difference (F-APW×PCG−1) is illustrated in
In
From the two linearly drawn lines of APW×PCG−1 and APW×PCG, variation behaviors of the higher harmonic components of the cardiac apex beat and the random vibration of the heart sound were detected, and it was possible to identify the disappearing frequency of the higher harmonic components of the cardiac apex beat, the heart sound appearing frequency, and BF.
Next, the focused points for the studies of the charts in
In
In APW×PCG in
BFs in
Further, BF did not correlate with age, SBP, DBP, or BMI. Further, in the study of BF in a control group of twenty young people, a high correlation was observed between BF found by STFT and BF decided from a BF-HR curve (R2=0.85). That is, since MAV (micro acoustic vibration) output from the trunk is collected using the biological signal detection sensor 1 (4SR) of this embodiment, by the analysis by the biological signal analysis device 100, it is possible to identify BF and the cardiac apex beat wave and classify a waveform (a sine wave or a triangular wave) even from a subject in “a thick chest wall group”, the identification from which was conventionally difficult.
Therefore, by storing in advance the correlation charts illustrated in
(Relationship of Heart Rate Variability Fluctuation Characteristic with Heart Rate and BF)
A correlation between the fluctuation characteristic of heart rate variability and the heart rate found from ECG in the APW measurement experiment was studied in fifty subjects. As illustrated in
The fluctuation of the heart rate variability is 1/f when the heart rate is 70 to 80/min in the resting state. When this is collated with BF, the fluctuation of the heart rate variability is 1/f when the aforesaid BF is around 8 to 12 Hz, and when BF has a value near this range, it can be inferred that the heart is moving efficiently. Therefore, by storing the correlation chart illustrated in
Next, a second embodiment of the present invention will be described based on
The cardiac apex beat waveform extracting means 160 filters biological signal data with the upper limit value being set to the boundary frequency (BF) identified by the boundary frequency identifying means 120, and finds a waveform of vibration generated by cardiac apex beat. In this embodiment, for the filtering, a 0.5 Hz to BF bandpass filter is used as described above. The cardiac apex beat waveforms (CAB) and the heart sound waveforms (CAS) are drawn in
The waveform classifying means 170 classifies a waveform of vibration generated by cardiac apex beat (cardiac apex beat waveform) output from the cardiac apex beat waveform extracting means 160. In this embodiment, for the classification, used is a combination of (A) a means that classifies the waveform by a mathematical approach using the Fourier series expansion, and (B) a means that classifies the waveform based on data obtained by an engineering approach. In this embodiment, cardiac apex beat waveforms of healthy people are classified into five using the two approaches, which will be described in detail later.
The means that classifies the waveform based on the data obtained by the engineering approach classifies the waveform using at least one of the following four determination elements. For higher classification accuracy, a combination of two or more of the determination elements is preferably used for the classification.
(1) Data indicating the shape of a time waveform in a predetermined time range of the aforesaid waveform of the vibration generated by the cardiac apex beat (the charts (b) in
(2) Data of a graph of the time waveform, in which frequency and power spectrum are taken on the horizontal axis and the vertical axis, the data being obtained by the frequency analyzing means (the charts (c) in
(3) Image data in a predetermined time range which is the analysis result of the short time Fourier transform of the time waveform and shows the time, the frequency, and the degree of the power spectrum variation (the charts (d) in
(4) Data regarding the type and period of a waveform, in a predetermined time range, which data is extracted from a correlogram (the charts (e) in
Out of the above, it is first required that the time waveforms of the healthy people are classified into five waveforms by the mathematical approach using the Fourier series expansion, and as its verification, amplitude, frequency, and time phase are analyzed by the engineering approach, and data corresponding to the five waveforms are found from the shape data of the time waveform of Front CAB (0.5-BF) (the charts (b) in
Information on the classification results of the waveforms classified by the waveform classifying means 170 is stored in the correlation data storage unit 150. Here, the information of the waveform classification result is stored in association with information regarding the health condition (the presence/absence of a disease such as a heart disease, the identification of a disease, whether the physical condition is good or not, and so on). To create this data, cardiac apex beat waveforms of as many people as possible are measured and they are associated with the health conditions at the time of the measurement.
The measurement-time condition inferring means 130 causes the waveform classifying means 170 to classify a measurement-time waveform of vibration generated by cardiac apex beat of a measurement subject which waveform is found by the cardiac apex beat waveform extracting means 160, and collates the measurement-time waveform classification result with the information of the waveform classification result stored in advance by accessing the correlation data storage unit 150. Consequently, the health condition corresponding to the measurement-time waveform classification result is output.
As illustrated in
In this case, the measurement-time condition inferring means 130 receives the measurement-time waveform information of the vibration generated by the cardiac apex beat of the measurement subject, and by using the inference model created by the model creating means 180, outputs a value of the measurement-time health condition based on the obtained measurement-time waveform information. Such a configuration to combine the machine learning improves the accuracy of inferring the health condition from the waveform information. Note that learning-target feature quantities of the waveform information are amplitude, frequency, time phase, and so on constituting each of the waveforms illustrated in
Note that the other configuration including the biological signal detection sensor 1, of the second embodiment is the same as that of the first embodiment.
Hereinafter, the methods of the waveform classification and evaluation will be more specifically described.
[Evaluation Method of Time Waveform of Front CAB]
(Analysis of Time Waveform of Front CAB)
Next, the evaluation method of the Front CAB waveform will be described.
The aforesaid apex cardiogram (ACG) is known as one of the physical diagnostic tools serving as a criterion of the non-invasive evaluation and as an index of disease condition and treatment, regarding the movement of the cardiovascular system, but a supine posture is not an absolute requirement as a posture during examination, and a sitting posture or a left lateral recumbent posture is sometimes taken. Therefore, if information regarding a cardiac apex beat normal pattern is available for each of the measurement postures including all of a supine posture, a left lateral recumbent posture, and a sitting posture, it is possible to study the relationship with the waveform analysis and the time phase analysis of PCG and ECG.
For this reason, Front CAB (0.5-BF), which is a waveform in a 0.5 to BF frequency band extracted from the frequency analysis results of the F-APW waveform measured by a microphone placed on the chest wall, is generated, and this is compared with a clinically known cardiac apex beat waveform. As focusing points of the waveform analysis, attention is given to “tapping” of an E-wave which is short-duration pulsation recognized only in early systole even by ocular inspection and palpation, and to a long-duration, strong “sustained”. Here, for the waveform analysis, the shape of the waveform of Front CAB (0.5-BF) and the time phase when the peak of ACG appears are used.
The evaluation of the amplitudes of the time waveforms in Front CAB (0.5-BF) is performed in the time phase of systole. An A-wave originating in left atrial contraction and a trace of the A-wave are confirmed in Front CAB (0.5-BF). A bottom value after the A-wave is considered as a C point of ACG, and the total amplitude up to the peak value corresponding to the E-wave of ACG is measured and is defined as the amplitude of the systolic wave. This total amplitude from the C point up to the E-wave in ACG is calculated also in the time waveforms of Front CAB (5-BF), Rear CAB (0.5-BF), and Lumbar CAB (0.5-BF). When the amplitudes of the time waveforms measured at the respective positions are compared, that of Front CAB (0.5-BF) is used as a reference value. An amplitude ratio is a ratio in each waveform. The bottom value of the amplitude of the time waveform is present between a P-wave and an R-wave in an electrocardiogram. The peak value of the amplitude of the time waveform is present in a range after the R-wave up to the starting point of second sound of PCG. Next, the time waveforms of Front CAB (0.5-BF) of fifty subjects are classified.
(Higher Harmonic Component Analysis Method Using Short Time Fourier Transform)
Next, a wave steepness analysis method and a higher harmonic component analysis method using the short time Fourier transform will be described.
Targets of the study are higher harmonics forming the time waveform of Front CAB (0.5-BF), and it was studied at which position of the time waveform the higher harmonic components are superimposed, using Formula (4). In the analysis, the time waveform is divided into those in frequency bands of 1 to 4 Hz and 4 to 16 Hz. The number of window points for the 1 to 4 Hz waveform analysis was 29=512. The 512 points correspond to 0.512 seconds because the sampling frequency is 1000 Hz, and according to the calculation, a component of the power spectrum of 1.95 Hz appears. A shift width of a Hanning window function is set to 51 points. At this time, frequency information is created every 51/1000 =0.051 seconds.
The number of window points for capturing an E-wave formed of higher harmonics was set to 27=128. Because the sampling frequency is 1000 Hz, 128 points correspond to 0.128 seconds, and a component of the power spectrum appears at 7.8 Hz. A reason why the resolution was set to seconds (7.8 Hz) is that the duration of systole in the case of 90/min heart rate is assumed to be seconds, and the required resolution of a detection frequency for detecting the peak in the duration of the systole is set to 0.13 seconds. In this case, the shift width of the Hanning window function was set to 12 points. At this time, the frequency information is created every 12/1000=0.012 seconds.
(Vibration Discrimination Using Autocorrelation Function)
A correlogram is drawn using ACF (Autocorrelation function), and waveform distinction is performed. In the correlogram, the vertical axis represents ACF and the horizontal axis represents time. ACF in each time interval T is calculated, and the type and period of the waveform are captured from the correlogram. Note that ACF, that is, Rx(τ) is represented by the following.
τ represents time interval, that is, lag, and T represents total interval.
A normalized autocorrelation function takes the maximum value 1 when τ=0, and this value is a mean square value. When the ACF value is 0.2 or less, it is in an interval of 95% or higher confidence of the absence of autocorrelation. In the vibration type analysis, ten-second data is used for determination. [Evaluation Result of Time Waveforms of Front CAB].
Using Front CAB (0.5-BF) measured from the fifty subjects, normal patterns in a sitting posture were classified and verified.
The method of the classification into the five waveform patterns is based on the Fourier series expansion. However, the mathematical expression is not suitable for intuitive understanding. Therefore, in studying the waveforms of Front CAB, the characteristics of the waveforms of Front CAB were classified, using PSD (power spectral density) when focus is on amplitude, using STFT when focus is on amplitude and phase, and using a correlogram when the presence/absence of periodicity is studied.
In this specification, the five waveforms resulting from the classification are named a pseudo sine wave (S-wave), a pseudo triangular wave group (T-wave (triangular wave), C-wave (cosine wave), G-wave (Gauss wave)), and a pseudo full-wave rectified wave (R-wave). Out of these, the S wave and the G-wave correspond to “sustained” and the T-wave, the C-wave, and the R-wave correspond to “tapping”, as is said in ocular inspection/palpation, in view of the duration of a systolic wave. Note that the prefix “pseudo” is appended because they have similar waveform shapes to those of a sine wave, a triangular wave, and a full-wave rectified wave though not complete sine wave, triangular wave, and full-wave rectified wave.
What are necessary for the waveform analysis are a waveform of a fundamental harmonic, and amplitudes, frequencies, and phases of waveforms near the bottom and the peak.
In (c) in
In more detail, based on PSD, the classification into the following waveform patterns is possible.
First, a systolic E-wave indicated by the S-wave of Front CAB resembles a sine wave, with less contribution of higher harmonic components. That is, only a fundamental frequency component is dominant (the “X” pattern illustrated in the charts).
Next, a systolic E-wave indicated by the T-wave, the C-wave, and the G-wave of Front CAB is in the process of converging to a triangular wave and is classified as a triangular wave. That is, odd-order frequency components are dominant, and their amplitudes become smaller in proportion to a reciprocal of the square of the order as the order becomes higher (the “Y” pattern illustrated in the charts).
In the R-wave of Front CAB, because the amplitude of a fundamental frequency component is small and has a small difference from the amplitude of a component of a frequency that is an integral multiple of the frequency of the fundamental frequency component, a triangular wave part is replaced by a full-wave rectified waveform of a sine wave. In corresponding PSD, an even-order frequency component is dominant, and the amplitude becomes smaller as the order becomes higher (the “Z” pattern illustrated in the charts).
STFT is suitable for studying the influence of phase information, and combining PSD and STFT improves the accuracy of the waveform classification. Using the correlogram with PSD and STFT to study periodicity between the waveforms of the A-wave and the E-wave makes it possible to classify a given waveform as a variant of a normal waveform from the viewpoint of amplitude, phase, and periodicity.
The STFT charts (d) in
In the waveform I (S-wave), there are three time zones where the fundamental harmonic appears, there are two time zones where 5 to 10 Hz components appear, and an appearance interval of an atrial bulging wave and a systolic wave is short, which indicates a square wave. In the waveform II (T-wave), the appearance of 5 Hz or lower components concentrates in a 2.3 to 2.5 second period in STFT, 5 to 10 Hz components concentrate in a 2.5 to 2.7 second period, and these components are roughly within a systolic time phase, which indicates a triangular wave. Similarly, the waveform III (C-wave) is also within the systolic time phase, and an interval between an atrial bulging wave and a systolic wave is short, which indicates a triangular wave. In the waveform IV (G-wave), a peak-to-peak interval of the wave in 5 to 10 Hz is close to that of the waveform I (S-wave), its shape is a combined shape of the S-wave and the T-wave, and it as a whole presents the feature of a triangular wave. In the waveform V (R-wave), a peak value in STFT in 5 Hz or lower presents a declining contour line, which indicates that the power spectra of the fundamental harmonic up to higher harmonics are equal in magnitude, further there are three or more peak values in not lower than 5 Hz nor higher than 10 Hz, and the appearance time phase is the same as that of a component of Hz or lower, which indicates the feature of a full-wave rectified waveform.
Using the correlograms (e) in
[Theoretical Ground and Objective Classification Criteria of Waveform Classification]
Here, the mathematical approaches to a cardiac apex beat wave using the Fourier series expansion will be described in detail. It has been shown that cardiac apex beat waves of healthy people are classified into five by these mathematical approaches, and PSD, STFT, and the correlogram described above serve as engineering approach-based supports for the classification into five.
First, let a given periodic function be f(t), then f(t) can be expanded to the following Fourier series.
Note that T is period, and ω=2π/T.
[Formula 7]
f(t)=α0Σn=1∞(αn cos nωt+bn sin nωt) (7)
The S-wave corresponds to a state where a first-order component is dominant in Formula (7) because its waveform shape resembles that of a sine wave. However, since it also includes some higher harmonic components, it can be interpreted from the STFT analysis result that its main component is a fundamental frequency accompanied by the high-order frequency components.
This agrees with the observation result of the short duration “Sustained” in systole.
It is thought that the T-wave can be expanded to the Fourier series as follows because it resembles a triangular wave.
Formula (8) shows that an expansion coefficient sharply decreases as the order of n becomes higher.
It also shows that the order that appears as a spectrum is an odd order.
This can explain the shape of PSD.
The influence of waveforms of higher harmonics other than the sine wave accounts for the result that the autocorrelation coefficient shown by the correlogram has a peak value when t=0 and the absolute value of the autocorrelation coefficient is smaller than 1 at the other times.
This means that a short-duration, sharp waveform originating in ventricular contraction in early systole is superimposed on the sine wave, which agrees with the observation result of short-duration “Tapping” in early systole.
Of the C-wave, a waveform originating in ventricular contraction can also be discussed in the same manner as is discussed about the T-wave. Its difference from the T-wave is that in the Fourier series expansion of Formula (8), a first-order component is more dominant in the C-wave than in the T-wave.
It is seen from the result of the STFT analysis that a waveform originating in ventricular contraction is formed mainly of a fundamental frequency and a second harmonic and of frequency components up to the seventh harmonic, and the correlogram agrees with the observation result of short-duration “Tapping” in early systole similarly to the T-wave. From the STFT analysis, the G-wave is considered as a combination of the S-wave and the T-wave, and a first-order component is dominant therein. Since it more resembles the S-wave rather than the T-wave, its duration in systole is long, which agrees with the observation result of “Sustained”.
In the R-wave, as for a waveform originating in ventricular contraction, the waveform itself is not unimodal but multimodal, but under the hypothesis that it resembles a full-wave rectified waveform, it can be expanded to the following Fourier series.
From Formula (9), it is considered that the R-wave starts with the second-order component with an even-order harmonic being superimposed thereon. An expansion coefficient becomes smaller as n becomes larger. It is thought that the R-wave has the shape illustrated in the PSD chart because the R-wave also includes the first-order component.
The analysis of the correlogram shows a reason why components other than the first-order component are dominant as in T-wave. This agrees with the observation result of short-duration “tapping” in early systole.
As described above, from the correlograms, the waveforms are classified into the S-wave and the G-wave which have the primitive period of the sine wave, and classified into the T-wave, the C-wave, and the R-wave in which the higher harmonics are prominent. Further, when attention is paid to the difference in the autocorrelation function with respect to lag, the S-wave and the G-wave can be distinguished on the+side of the correlogram, and the T-wave and the C-wave can be distinguished on the−side of the correlogram, and the R-wave can be distinguished by the notch indicating multimodality.
That is, in the classification using the correlogram, attention is paid to sinusoidal periodicity for the S-wave and the G-wave, while attention is paid to distorted wave periodicity for the T-wave, the C-wave, and the R-wave. The periodicities of systolic waves are classified by the autocorrelation function on the+side of the correlogram, and the periodicities of the atrial bulging waves and other waveforms are classified by the autocorrelation function on the−side of the correlogram. It was possible to sort the S-wave and the G-wave according to a difference in autocorrelation with respect to lag, and sort the T-wave, the C-wave, and the R-wave according to the value (0, −0.2, and −0.5 in the above example) of the autocorrelation coefficient on the−side. Then, for the irregularities of the atrial bulging wave and the other waveforms, the drawn shape (notch) of the autocorrelation coefficient was used. Since the R-wave is multimodal, it had two notches or more.
Using the charts (b) to (d) in
Therefore, in this embodiment, when the classification result of the cardiac apex beat wave of a measurement subject obtained by the waveform classifying means 170 corresponds to any of the aforesaid five patterns which are waveform patterns of healthy people, the measurement-time condition inferring means 130 determines that the subject is healthy and outputs this result, and when the classification result does not correspond to any of these, it determines that there is an abnormality finding regarding a heart disease and outputs this result.
Lumbar CAB (0.5-BF) to the amplitude of Front CAB (0.5-BF). Here, the amplitude ratio of Front CAB (5-BF) to that of the Front CAB (0.5-BF) was 0.2 because it does not include a 0.5 Hz to 5 Hz frequency component. The amplitude ratio of Rear CAB (0.5-BF) was 1.1, and the amplitude ratio of Lumbar CAB (0.5-BF) was 2.6. The mean value and the standard deviation were calculated using the amplitude ratios of the fifty subjects.
Next, the power spectra in
From the above, according to the present invention, the following specific effects were obtained.
The conventional 3SR is a sensing system that is developed for on-vehicle use and damps external vibration by the air damping effect of the air pack to extract first sound and second sound of heart sound including 20 Hz or higher frequency components and periodic components of heart rate variability. On the other hand, 4SR of the present invention is developed for use in rooms, and it has been found out that 4SR is suitable as a sensor for capturing the period and the amplitude of cardiac apex beat of around 1 to 2 Hz.
The waveform of F-APW was similar to that of a cardiac apex beat wave. Then, from F-APW, it was possible to generate waveforms of a cardiac apex beat wave and heart sound. It is considered that the waveforms of R-APW and L-APW became different due to a resonance phenomenon caused by pressure pulsation. However, from R-APW and L-APW whose waveforms were different from that of F-APW, it is also possible to extract vibration information (CAB) regarding the movement of the heart including cardiac apex beat and heart sound (CAS) by using the vibration analysis, and the possibility of clinical significance is suggested.
By combining the BF extraction method based on the disappearing frequency of the higher harmonics in F-APW and the BF extraction method using the STFT method, the BF-heart rate correlation chart is created, BF can be represented by the quadratic function of heart rate according to the actual measurement results, and it is possible to decide BF based on the heart rate in real time only with data of F-APW, without using PCG. Further, a possibility is suggested that the use of the visualization method by the STFT method enables not only the vibration analysis of the vibration of the heart and pressure pulsation of the aorta but also diagnosis using a change degree.
The experiment result of the fluctuation of heart rate variability showed that when the heart rate was 70 to 80/min in the resting state, the fluctuation was 1/f and the movement of the heart is efficient. The experiment result of the fluctuation of heart rate variability showed that when BF was around 10 Hz, the fluctuation of the heart rate variability was 1/f and the movement of the heart is efficient. The use of these characteristics has a possibility that a change in the circulatory organ function caused by drug can be quantified.
As described above, BF is decided based on the heart rate (for ten seconds in this case). The tendency of the value of the fluctuation (for ninety seconds) of the heart rate variability obtained from BF agreed with the tendency of the value of the fluctuation (for ninety seconds) of heart rate variability directly obtained from the heart rate (for ninety seconds). The result found when their relational formula was substituted supported the validity of BF.
Further, by applying the 0.5 Hz to BF band-pass filtering to APW, it was possible to draw a waveform resembling that of ACG, and it was possible to classify Front CAB (0.5-BF) representing the cardiac apex beat roughly into the five waveforms, by the frequency analysis, the time-series analysis method using the STFT method, and the waveform analysis method using the correlogram.
When heart rate variability occurs under respiratory depression in the resting state, BF shifts on the regression curve of the BF-HR correlation chart (see
Further, in a supine posture and left and right lateral recumbent postures, the same BF-HR correlation charts were obtained as those obtained in the sitting posture (see
According to the present invention, signals belonging to two frequency bands are separated through the identification of BF.
Specifically, since signals corresponding to cardiac apex beat and heart sound can be separated in APW, it is possible to quantitatively recognize the results of different types of operations of cardiac activity by observing one physiological signal, which makes it possible to obtain useful complementary information regarding the cardiac condition. Further, since the signal separation increases waveform resolution, there is a possibility that, for example, a frequency as the boundary between two frequency bands, higher harmonic components of Front CAB near the boundary, and a dominant frequency component of Front CAS serve as new parameters of the cardiocirculatory system.
Further, Front CAB has pressure information of a systolic wave in an apex cardiogram, and Front CAS has time phase information relating to the cardiac cycle, and the combination of the pressure and the time phase information is thought to be usable for diagnosis from the complementary information.
For example, first heart sound can be identified from Front CAS, and a notch corresponding to the first sound of an anacrotic limb of a systolic wave can be observed both in a wide range and a narrow range in STFT illustrated in (b), (c) in
Measuring data along the aorta makes it possible to separate the data into those originating in cardiac apex beat and heart sound in hemodynamics of the aorta, and may make it possible to obtain new findings regarding the hemodynamics of the aorta.
As is apparent from the above, it is seen that by extracting APW (Acoustic Pulse Wave) of 0.5 to 80 Hz and applying predetermined filtering, it is possible to generate an ACG-equivalent waveform (Cardiac Apex Beat: CAB) and a heart sound-equivalent waveform (Cardiac Acoustic Sound: CAS).
Further, CAB and CAS are different in the generation source and transmission route, and vibration resulting from the superposition of these two vibrations is APW. Considering that systolic performance is a characteristic unique to the cardiac muscles and is independent of afterload and preload, the waveform of CAB in systole has information regarding the influence of chemical stimulation and hormone at the time of contraction. Considering that CAS has movement information of blood ejected from the ventricle in systole, APW has information on preload and afterload and on systolic performance.
Therefore, APW, CAB, and CAS have a lot of information regarding sound, vibration, and pressure, and observing them together with a polarization phenomenon that ECG has enables diagnosis of cardiac failure.
It has been also found out that the real-time measurement of APW, CAB, and CAS additionally provides clinically significant information on acoustic vibration and pressure, and the engineering approach having been used for diagnosis of a failure of a machine is applicable to health diagnosis. Further, the 4SR system of the present invention that measures APW serves as a physical diagnostic tool capable of repeatedly obtaining biological information non-invasively over a long period, and is expected to have wide applications including AI analysis including the aforesaid mechanical learning.
1 biological signal detection sensor (4SR)
10 three-dimensional knitted fabric
20 housing film
30 microphone
40 case
50 gel
100 biological signal analysis device
110 frequency analyzing means
120 boundary frequency identifying means
130 measurement-time condition inferring means
140 heart rate calculating means
150 correlation data storage unit
160 cardiac apex beat waveform extracting means
170 waveform classifying means
180 model creating means
1000 biological signal detection sensor (3SR)
| Number | Date | Country | Kind |
|---|---|---|---|
| 2020-180963 | Oct 2020 | JP | national |
| 2021-107109 | Jun 2021 | JP | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/JP2021/039924 | 10/28/2021 | WO |
