CN117547228B - Data quality analysis method of biological pulse diagnosis signal - Google Patents
Data quality analysis method of biological pulse diagnosis signal Download PDFInfo
- Publication number
- CN117547228B CN117547228B CN202410028486.3A CN202410028486A CN117547228B CN 117547228 B CN117547228 B CN 117547228B CN 202410028486 A CN202410028486 A CN 202410028486A CN 117547228 B CN117547228 B CN 117547228B
- Authority
- CN
- China
- Prior art keywords
- zero crossing
- pulse
- data
- waveform
- acquisition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 19
- 238000004458 analytical method Methods 0.000 title claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 14
- 238000004891 communication Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000241 respiratory effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000013186 photoplethysmography Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 230000005534 acoustic noise Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229910052744 lithium Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000002688 persistence Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7221—Determining signal validity, reliability or quality
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
Abstract
The application relates to the technical field of pulse diagnosis signal acquisition, in particular to the technical field of continuity and signal quality analysis of acquired pulse diagnosis signals, and particularly discloses a data quality analysis method of biological pulse diagnosis signals, which comprises the steps of data acquisition, data filtering, continuity judgment, zero crossing frequency judgment, special pulse condition characteristic matching and secondary zero crossing frequency K judgment 2 And if the pulse signal quality falls into the threshold range Q of the preset qualified signal quality, repeating the steps to finally obtain an analysis conclusion whether the current pulse signal quality is credible or not, and prompting whether the current pulse signal quality has acquisition conditions or not. The pulse signal continuity and clutter mixing condition is judged by gradually calculating the zero crossing times of the pulse signals step by step and matching the zero crossing times, and signals which are discontinuous, impure or not in accordance with the biological pulse characteristic range are taken as low-quality or unqualified pulse signals to be removed, so that the reliability of the adopted pulse signals is ensured.
Description
Technical Field
The invention relates to the technical field of pulse diagnosis signal acquisition, in particular to the technical field of signal quality analysis and continuity of acquired pulse diagnosis signals, and particularly relates to a data quality analysis method of biological pulse diagnosis signals.
Background
The pulse signal is a basic characteristic signal of the living being, and the pulse signal can change along with different environments and different states of the living being, so that the change of the pulse signal can objectively reflect the change of different activities or states of the living being in different time periods. In the prior art, the wearable device for carrying out pulse diagnosis or statistics recording of exercise quantity based on pulse acquisition has various mature products, such as a portable fingertip detection pulse diagnosis instrument disclosed in publication No. CN216907945U, and relates to the technical field of pulse diagnosis instruments. For another example, publication number CN108065922a discloses that the computer-intelligent pulse-taking fixing device can realize accurate positioning and digital positioning of pulse-taking, and can be matched with advanced pulse-taking probes and intelligent devices to realize accurate intelligent pulse-taking and networking big data pulse-taking.
Because of pulse diagnosis modes, rules are described in Chinese medical theory, the existing pulse diagnosis technology is mainly aimed at more portable pulse acquisition and intelligent application, and the fresh technology is aimed at analyzing the quality of the acquired pulse signals, that is, at present, the field of pulse diagnosis technology is almost not subjected to credibility research on whether the pulse signals can truly and objectively reflect biological states. That is, in the prior art, the biological pulse signal is directly collected in a sensor manner, and the pulse signal collected by default is completely trusted and is directly subjected to subsequent analysis or application. In fact, although the biological pulse signals are easily obtained by the prior art, various non-biological pulse signals are easily mixed in the process of acquisition, and the non-biological pulse signals can be collectively called noise signals, and mainly comprise acquisition devices, such as a sensor and a biological body surface, which have a variable gap, and vibration or acoustic noise generated in the environment can confuse the pulse signals. Therefore, no matter what the end use based on pulse acquisition is, the method and the device aim to solve the problem that the acquired pulse information is distorted in credibility due to the fact that other clutter noise is mixed, namely the method and the device are used for analyzing the credibility of the quality of the acquired pulse signals, and the problem that the real pulse information cannot be objectively reflected due to the fact that noise or special pulse signals are mixed due to the fact that the used pulse signals are distorted is avoided, so that whether the current acquisition environment can meet acquisition conditions is determined.
Disclosure of Invention
In order to solve the problem that the reliability is distorted due to the fact that other clutter noise is mixed in pulse information acquired through a sensor in the prior art, the application provides a data quality analysis method of biological pulse diagnosis signals, which is used for analyzing the quality reliability of the acquired pulse signals, so that pulse signals which are not reliable or do not meet the subsequent use conditions are discarded or selected not to be acquired, only the pulse information with the reliability is used, and therefore whether the current pulse acquisition environment can meet the acquisition conditions is determined. The invention only analyzes the pulse signal data to evaluate whether the pulse signal acquired under the current environment/state is truly credible or usable, and does not prejudge the state of the biological body, such as any other aspect of human or animal.
In order to achieve the above purpose, the technical scheme adopted in the application is as follows:
the invention provides a data quality analysis method of biological pulse diagnosis signals, which comprises the following steps:
step STP100, data acquisition, namely acquiring biological pulse original data through an acquisition terminal attached to a biological body surface, wherein the duration time period of the acquired data is T 1 =2s, obtain analog data C 0 ；
Step STP200, data filtering, and filtering the analog data C 0 The respiratory signal and the power frequency noise signal are filtered and temporarily stored in the acquisition terminal;
step STP300, continuity judgment, judgment of analog data C 0 Whether the signal duration of (2) is up to 2 seconds; if not, executing step STP100, and if not, executing step STP400;
step STP400, zero crossing number judgment, judgment of analog data C 0 Number of initial zero crossings K of waveform 1 And according to the initial zero crossing number K 1 Searching whether the zero crossing times and the initial zero crossing times K exist 1 Equal special pulse characteristics;
if the special pulse condition characteristics exist, executing a zero crossing compensation step;
if no special pulse condition exists, the secondary zero crossing times K are carried out 2 Counting;
step STP500, judging the number of times K of the secondary zero crossing 2 Falls within a threshold range Q of preset qualified signal quality, if K 2 E Q, then determine the current analog data C 0 Is qualified in signal quality; if K 2 ∉ Q, then determine the current analog data C 0 Is unqualified;
step STP600, repeating steps STP 100-STP 500 three times to obtain three times of double zero crossing K 2 ，K 3 And K 4 ；
If K 2 ，K 3 And K 4 If any two of the gesture is qualified, prompting a user to keep current gesture acquisition until the current gesture acquisition is finished through a visual intelligent terminal which establishes wireless communication connection with the acquisition terminal;
if K 2 ，K 3 And K 4 If any two or more of the wireless communication terminals are unqualified, the visual intelligent terminal which establishes wireless communication connection with the acquisition terminal prompts the user that the current environment or gesture does not accord with the acquisition condition.
In order to increase the calculation speed, preferably, in the step STP400 zero crossing number judgment, the analog data C is judged 0 Number of initial zero crossings K of waveform 1 Also previously include the analog data C 0 By analogue to digital conversion into digital data C 1 And band-pass filtering to obtain digital data C 2 Step (C) of judging the digital data C 2 The zero crossing times method of (2) is as follows:
and counting the times of F (T) F (t+1) < 0, wherein F (T) is the value of a sampling point, T is the sampling times in a time period T1, and the sampling frequency fs=100 Hz, and t=1, 2 and 3 … 200.
Further preferably, the step STP400 searches for the presence of zero crossings and the simulationData C 0 Number of initial zero crossings K 1 The step of the same special pulse condition characteristics specifically comprises the following steps:
step STP401, determining analog data C 0 A waveform period T of (2);
step STP402, calculating waveform period T and simulating data C 0 The number k of waveform zero crossings;
step STP403, searching for special pulse waveform preset in the acquisition terminal and calculating the number of zero crossings k in the waveform unit period n And analog data C 0 The number x of special pulse waveform matched with the number k of the zero crossing times of the waveform n ；
Way1: when x is n When=1, the compensation value B of the currently matched special pulse waveform is executed n Wherein the number of zero crossings after compensation is k+B n ；
Way2: when x is n If not less than 2, counting the upper zero crossing times m and the lower zero crossing times n respectively and performing second matching, wherein the second zero crossing matching process is to make m, n and k respectively n Matching the upper zero crossing times and the lower zero crossing times contained in the first set, wherein k=m+n;
if the number x of the waveforms of the special pulse condition is matched for the second time n When=1, then Way1 is performed;
if the number x of the waveforms of the special pulse condition is matched for the second time n And if not less than 2, executing Way3:
wherein x is n N=1, 2,3 for the number of special pulse waveforms for the nth match;
way3: matching m, n and k n The order in which the upper zero crossing number and the lower zero crossing number are present is included,
if the number x of the special pulse waveforms is matched for the third time n When=1, then Way1 is performed; otherwise, step STP500 is performed.
The beneficial effects are that:
the pulse signal continuity and clutter mixing condition is judged by gradually calculating the zero crossing times of the pulse signals step by step and matching the zero crossing times, and signals which are discontinuous, impure or not in accordance with the biological pulse characteristic range are taken as low-quality or unqualified pulse signals to be removed, so that the reliability of the adopted pulse signals is ensured.
The invention adopts the highest three zero crossing times to match, can carry out compensation processing on the signals with special pulse characteristics, and avoids the problem of result distortion caused by the elimination of special pulse signals as clutter signals or the evaluation based on normal pulse conditions.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive faculty for a person skilled in the art.
Fig. 1 is a flow chart of the pulse signal quality analysis and judgment according to the present invention.
Fig. 2 is a flow chart diagram of zero crossing compensation.
FIG. 3 is a waveform diagram of an embodiment.
FIG. 4 is another pulse waveform diagram of an embodiment.
FIG. 5 is a plot of sampling points for a pulse condition according to an embodiment.
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
In the description of the present application, it should be noted that, if the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," "outer," and the like indicate an azimuth or a positional relationship based on that shown in the drawings, or an azimuth or a positional relationship that a product of the application conventionally puts in use, it is merely for convenience of describing the present application and simplifying the description, and does not indicate or imply that the device or element to be referred to must have a specific azimuth, be configured and operated in a specific azimuth, and thus should not be construed as limiting the present application. Furthermore, the terms "first," "second," and the like in the description of the present application, if any, are used for distinguishing between the descriptions and not necessarily for indicating or implying a relative importance.
Furthermore, the terms "horizontal," "vertical," and the like in the description of the present application, if any, do not denote a requirement that the component be absolutely horizontal or overhang, but rather may be slightly inclined. As "horizontal" merely means that its direction is more horizontal than "vertical", and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.
In the description of the present application, it should also be noted that, unless explicitly stated and limited otherwise, the terms "disposed," "mounted," "connected," and "connected" should be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the terms in this application will be understood by those of ordinary skill in the art in a specific context.
Example 1:
the embodiment provides a data quality analysis method of biological pulse diagnosis signals, which is shown in fig. 1, and includes the following steps:
step STP100, data acquisition, namely acquiring biological pulse original data through an acquisition terminal attached to a biological body surface, wherein the duration time period of the acquired data is T 1 =2s, obtain analog data C 0 The method comprises the steps of carrying out a first treatment on the surface of the The acquisition terminal can adopt a wearable device to acquire human body pulse based on the existing photoplethysmography (PPG) sensor or pressure wave sensor.
Step STP200, data filtering, and filtering the analog data C 0 The respiratory signal and the power frequency noise signal are filtered and temporarily stored in the acquisition terminal; wherein, the setting of the filter can be customized according to the actual requirement to more finely and completely filter the analog data C 0 The technology is very mature, in the prior art, filters with different model parameters are already available, and different frequency bandwidths can be set according to requirements for filtering, which is not an improvement of the technology of the present invention, and belongs to the mature prior art, and not described in detail herein.
Step STP300, continuity judgment, judgment of analog data C 0 Whether the signal duration of (2) is up to 2 seconds; if not, executing step STP100, and if not, executing step STP400; the step is used for eliminating discrete pulse signals caused by fluctuation of a data acquisition environment, for example, in the acquisition process of an acquisition terminal, the surface of a living body cannot be in good contact with a sensor of the acquisition terminal, gaps exist, so that the discrete signals are caused, the bad pulse signals with the persistence of less than 2 seconds can be abandoned, the subsequent analysis can not be carried out, and the unreal feedback in the subsequent practical application caused by the distorted discrete signals is avoided. When the acquired pulse signals can be continuously acquired and have continuity reaching 2 seconds, the current acquisition environment is considered to be in accordance with the condition, and the acquired pulse signals are in accordance with the basic requirement under the time domain condition and can be used for subsequent quality analysis.
Step STP400, zero crossing number judgment, judgment of analog data C 0 Number of initial zero crossings K of waveform 1 And according to the initial passZero times K 1 Searching whether the zero crossing times and the initial zero crossing times K exist 1 Equal special pulse characteristics;
if the special pulse condition characteristics exist, executing a zero crossing compensation step;
if no special pulse condition exists, the secondary zero crossing times K are carried out 2 Counting;
the step is an improvement of the key technology of the invention, and the zero crossing times mentioned in the step refer to the acquired analog data C 0 The number of times the waveform formed in the time domain crosses the time axis is shown in figures 3-5. When the signal waveform intersects with the time axis, an intersection point is generated, the intersection point is marked as a zero crossing, the basic frequency of the pulse signal can be determined by counting the zero crossing times of the pulse signal in unit time, and in theory, the credibility of the current pulse signal can be analyzed by judging whether the zero crossing times are in the range of the preset zero crossing times or not; however, judging the pulse signal quality based on the continuous pulse signal statistics zero crossing times is not objectionable, because even though the pulse signal quality is subjected to fine filtering treatment, the zero crossing times of the pulse signal can be obviously beyond the normal value range because the living beings are in different physiological states, if the pulse signal quality is judged and eliminated only according to the zero crossing times, the special pulse signal can be eliminated as unqualified clutter, thus leading to real analysis distortion, and leading to the false elimination of the high-quality pulse signal which is mistakenly regarded as clutter. To solve this technical problem, in this embodiment, the currently acquired analog data C will be used in the zero-crossing determination process 0 Number of initial zero crossings K 1 Matching with special pulse characteristics preset in the acquisition terminal or an intelligent terminal capable of carrying out real-time information interaction with the acquisition terminal, and if special pulse matching conditions can be met, executing zero crossing compensation according to the matched special pulse, thereby avoiding the problem that the current pulse signal is misjudged and removed because of the fact that the current pulse signal belongs to the special pulse and the zero crossing times exceeds the limiting value. Specifically, assume that the current analog data C 0 Number of initial zero crossings K 1 After matching the special pulse condition, the zero crossing frequency of the same special pulse condition is 18, theThe compensation value of the special pulse condition is-12, then K is calculated after zero-crossing compensation is performed 1 =18+ (-12) =6, that is, the current K will be 1 =6 is included as the actual zero crossing number in the subsequent determination, and not actually acquired K 1 The determination of 18 is incorporated, so that the problem of misdetermination of a specific pulse condition caused by a specific physiological phase of the living being is avoided.
Step STP500, judging the number of times K of the secondary zero crossing 2 Falls within a threshold range Q of preset qualified signal quality, if K 2 E Q, then determine the current analog data C 0 Is qualified in signal quality; if K 2 ∉ Q, then determine the current analog data C 0 Is unqualified; as described above, the special pulse condition after zero crossing compensation, namely K 1 After=6 as the object of quality judgment, if the zero crossing frequency range still does not belong to the qualified quality, the currently acquired analog data C is recorded 0 And if the sample is unqualified, carrying out the next acquisition.
Step STP600, repeating steps STP 100-STP 500 three times to obtain three times of double zero crossing K 2 ，K 3 And K 4 ；
If K 2 ，K 3 And K 4 If any two of the gesture is qualified, prompting a user to keep current gesture acquisition until the current gesture acquisition is finished through a visual intelligent terminal which establishes wireless communication connection with the acquisition terminal;
if K 2 ，K 3 And K 4 If any two or more of the wireless communication terminals are unqualified, the visual intelligent terminal which establishes wireless communication connection with the acquisition terminal prompts the user that the current environment or gesture does not accord with the acquisition condition. The purpose of adopting the statistics of the number of times of zero crossing is to improve the feasibility and redundancy of measurement, thereby avoiding the difficulty in carrying out the statistics due to the too severe acquisition conditions or the requirement of the pulse signal quality, which is unfavorable for the implementation of the existing pulse acquisition-based application, and further reducing the practicability of the invention.
Example 2:
in this embodiment, in order to increase the calculation speed based on embodiment 1, the steps of this embodiment are further optimizedIn particular, in the step STP400 zero crossing number judgment, the analog data C is judged 0 Number of initial zero crossings K of waveform 1 Also previously include the analog data C 0 By analogue to digital conversion into digital data C 1 And band-pass filtering to obtain digital data C 2 Step (C) of judging the digital data C 2 The zero crossing times method of (2) is as follows:
counting the times of f (T) f (t+1) less than or equal to 0, wherein f (T) is the value of a sampling point, and T is the time period T 1 The number of samples in the same, the sampling frequency fs=100 hz, t=1, 2,3 … 200, and the number of times of two adjacent equal to 0 is calculated 1 time.
Compared with embodiment 1, the embodiment has lower hardware calculation requirement and lower energy consumption, and the matching essential technology in the technical scheme is to calculate the zero crossing times and to simulate data C 0 In other words, it is a waveform based on time domain, which causes calculation of analog data C 0 Is g (t) 0 ) When another g (t) 0 ) When=0, it corresponds to t 0 The number is zero crossing times, though the principle is simpler, the singlechip with simple functions cannot be completed in actual calculation, and the singlechip is required to be realized by a more complex functional circuit, however, the singlechip with simple functions does not have good practicability for the wearable device as an acquisition terminal after the hardware power consumption is increased, and the built-in lithium battery and the circuit board are required to be considered in consideration of the energy supply problem, so that the problem of occupation of physical space is required to be considered, and therefore, the power consumption of the circuit for acquiring and supplying energy can be realized by adopting small-size wearable devices such as a bracelet and a ring. Using digital data C 2 After the number of zero crossings due to digital data C 2 As shown in fig. 5, any zero crossing point is an intersection point between a connecting line of two adjacent sampling points and a time axis, so that positive and negative fractions are necessarily present in values of the two adjacent sampling points, and in special cases, a sampling point is 0 may be present, so long as the product of the two adjacent sampling points is less than or equal to zero, and the current digital data C can be obtained as long as the number of times that the statistical product is less than or equal to zero 2 Is not zero crossing number of adjacentThe values of the two sampling points are both positive numbers or both negative numbers, and no positive-negative condition exists, so that statistics of zero crossing times cannot be affected by the number of the other sampling points. The calculation basis is only simple operation and positive and negative values, and the comparison of the value 0 can be completed by a simple comparison circuit and a singlechip, so that the quick operation result can be realized by extremely low consumption of power consumption, and the requirement of small-size wearable equipment on an internal circuit can be met.
Example 3:
this embodiment is further refined based on any of the above embodiments, as shown in fig. 2, in step STP400, whether the zero crossing number and the analog data C exist is searched for 0 Number of initial zero crossings K 1 The step of the same special pulse condition characteristics specifically comprises the following steps:
step STP401, determining analog data C 0 A waveform period T of (2); the waveform period in this step refers to waveform units forming a pulse waveform, for example, in a continuous sine wave, the average number of zero crossings of a single waveform unit is 2, the number of zero crossings of the whole sine wave is 2n+1, and n is 360 °;
step STP402, calculating waveform period T and simulating data C 0 The number k of waveform zero crossings;
step STP403, searching for special pulse waveform preset in the acquisition terminal and calculating the number of zero crossings k in the waveform unit period n And analog data C 0 Digital data C 2 The number x of special pulse waveform matched with the number k of the zero crossing times of the waveform n The method comprises the steps of carrying out a first treatment on the surface of the For example, referring to fig. 3, the current zero crossing number k=3, it is searched whether k exists n A special pulse waveform of =3 exists; if there is no matched special pulse condition, it indicates x n If =0, directly executing step STP500 to perform the second zero crossing statistics, if x n If not 1, the following method is adopted according to the situation:
way1: when x is n When=1, the compensation value B of the currently matched special pulse waveform is executed n Wherein the number of zero crossings after compensation is k+B n The method comprises the steps of carrying out a first treatment on the surface of the For example, current k=3, matchCompensation value B of matched special pulse condition n -2, then the number of zero crossings after compensation is k=3+ (-2) =1, then the number of zero crossings valid is k=1;
way2: when x is n If not less than 2, counting the upper zero crossing times m and the lower zero crossing times n respectively and performing second matching, wherein the second zero crossing matching process is to make m, n and k respectively n Matching the upper zero crossing times and the lower zero crossing times contained in the first set, wherein k=m+n;
if the number x of the waveforms of the special pulse condition is matched for the second time n When=1, then Way1 is performed;
if the number x of the waveforms of the special pulse condition is matched for the second time n And if not less than 2, executing Way3:
wherein x is n N=1, 2,3 for the number of special pulse waveforms for the nth match;
way3: matching m, n and k n The same applies to the order of occurrence of the upper zero crossing times and the lower zero crossing times, as shown in fig. 3, taking k=3 as an example, which are the upper zero crossing, the lower zero crossing and the upper zero crossing in sequence according to the time sequence, then m=2, n=1, if the number of special pulse waveforms x is matched for the third time n When=1, then Way1 is performed; otherwise, step STP500 is performed.
The beneficial effects are that: in this embodiment, the successive step-by-step matching is adopted, and most pulse conditions do not belong to special pulse conditions, so that in a normal case, after the first initial zero crossing frequency calculation, statistics of the second zero crossing rate or zero crossing frequency is rapidly performed, and whether the current pulse data signal is a trusted signal is judged according to whether the counted zero crossing frequency belongs to a preset zero crossing frequency range, so that a complete data quality analysis is completed. Because the special pulse condition is not matched with a proper object when the special pulse condition is matched for the first time, the calculation aiming at zero crossing compensation is not carried out under most conditions, so that the calculation amount is greatly reduced, the power consumption of the calculation is reduced, and the invention has wide application space in the aspects of miniaturization design and application of the acquisition terminal and can be even accommodated in a conventional ring carrier. On the other hand, although the biological objects of the special pulse condition are few but not absent, such as a virtual pulse, a floating pulse, a short pulse, a flood pulse, a slippery pulse and the like, the invention creates the problem that the zero crossing times including but not limited to the special pulse condition can be preset in advance, including the zero crossing up-down sequence and zero crossing compensation, so that the unnecessary zero crossing times in the unit waveform period are brought into the follow-up application and are inverted to data distortion due to the specificity of the pulse condition after the invention is adopted for collection and analysis; the problem that the acquired pulse data is invalid data due to the fact that the overrun of the zero crossing times is used as clutter is solved, and meanwhile, the high efficiency of the conventional pulse condition calculation and the targeted and comprehensive inclusion of special pulse conditions are considered. After the acquired pulse signals are analyzed, the invention can embody the real and objective pulse frequency data of organisms and solve the problem of data distortion.
The foregoing description is only of the preferred embodiments of the present application and is not intended to limit the same, but rather, various modifications and variations may be made by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims (2)
1. The data quality analysis method of the biological pulse diagnosis signal is characterized by comprising the following steps:
step STP100, data acquisition, namely acquiring biological pulse original data through an acquisition terminal attached to a biological body surface, wherein the duration time period of the acquired data is T 1 =2s, obtain analog data C 0 ；
Step STP200, data filtering, and filtering the analog data C 0 The respiratory signal and the power frequency noise signal are filtered and temporarily stored in the acquisition terminal;
step STP300, continuity judgment, judgment of analog data C 0 Whether the signal duration of (2) is up to 2 seconds; if not, executing step STP100, and if not, executing step STP400;
step STP400, zero crossing number judgment, judgment of analog data C 0 Number of initial zero crossings K of waveform 1 And according to the initial zero crossing number K 1 Searching whether the zero crossing times and the initial zero crossing times K exist 1 Equal special pulse characteristics;
if the special pulse condition characteristics exist, executing a zero crossing compensation step;
if no special pulse condition exists, the secondary zero crossing times K are carried out 2 Counting;
step STP500, judging the number of times K of the secondary zero crossing 2 Falls within a threshold range Q of preset qualified signal quality, if K 2 E Q, then determine the current analog data C 0 Is qualified in signal quality; if K 2 ∉ Q, then determine the current analog data C 0 Is unqualified;
step STP600, repeating steps STP 100-STP 500 three times to obtain three times of double zero crossing K 2 ，K 3 And K 4 ；
If K 2 ，K 3 And K 4 If any two of the gesture is qualified, prompting a user to keep current gesture acquisition until the current gesture acquisition is finished through a visual intelligent terminal which establishes wireless communication connection with the acquisition terminal;
if K 2 ，K 3 And K 4 If any two or more of the wireless communication terminals are unqualified, prompting that the current environment or posture of the user does not accord with the acquisition conditions through a visual intelligent terminal which establishes wireless communication connection with the acquisition terminal;
searching in the STP400 whether there is zero crossing number and the analog data C 0 Number of initial zero crossings K 1 The step of the same special pulse condition characteristics specifically comprises the following steps:
step STP401, determining analog data C 0 A waveform period T of (2);
step STP402, calculating waveform period T and simulating data C 0 The number k of waveform zero crossings;
step STP403, searching for special pulse waveform preset in the acquisition terminal and calculating the number of zero crossings k in the waveform unit period n And analog data C 0 The number x of special pulse waveform matched with the number k of the zero crossing times of the waveform n ；
Way1: when x is n When=1, then perform current matchingCompensation value B of special pulse waveform n Wherein the number of zero crossings after compensation is k+B n ；
Way2: when x is n If not less than 2, counting the upper zero crossing times m and the lower zero crossing times n respectively and performing second matching, wherein the second zero crossing matching process is to make m, n and k respectively n Matching the upper zero crossing times and the lower zero crossing times contained in the first set, wherein k=m+n;
if the number x of the waveforms of the special pulse condition is matched for the second time n When=1, then Way1 is performed;
if the number x of the waveforms of the special pulse condition is matched for the second time n And if not less than 2, executing Way3:
wherein x is n N=1, 2,3 for the number of special pulse waveforms for the nth match;
way3: matching m, n and k n The order in which the upper zero crossing number and the lower zero crossing number are present is included,
if the number x of the special pulse waveforms is matched for the third time n When=1, then Way1 is performed; otherwise, step STP500 is performed.
2. The method for analyzing data quality of biological pulse diagnosis signal according to claim 1, wherein in the step STP400 zero crossing number judgment, the analog data C is judged 0 Number of initial zero crossings K of waveform 1 Also previously include the analog data C 0 By analogue to digital conversion into digital data C 1 And band-pass filtering to obtain digital data C 2 Step (C) of judging the digital data C 2 The zero crossing times method of (2) is as follows:
counting the times of F (T) F (t+1) < 0, wherein F (T) is the value of the sampling point, and T is the time period T 1 The sampling frequency fs=100 hz, t=1, 2,3 … 200.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410028486.3A CN117547228B (en) | 2024-01-09 | 2024-01-09 | Data quality analysis method of biological pulse diagnosis signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410028486.3A CN117547228B (en) | 2024-01-09 | 2024-01-09 | Data quality analysis method of biological pulse diagnosis signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117547228A CN117547228A (en) | 2024-02-13 |
| CN117547228B true CN117547228B (en) | 2024-03-19 |
Family
ID=89813146
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410028486.3A Active CN117547228B (en) | 2024-01-09 | 2024-01-09 | Data quality analysis method of biological pulse diagnosis signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117547228B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07148126A (en) * | 1993-11-30 | 1995-06-13 | Seiko Epson Corp | Pulse wave analyzing apparatus |
| JPH10248819A (en) * | 1997-03-13 | 1998-09-22 | Seiko Epson Corp | Pulse-wave diagnostic device |
| JPH11113863A (en) * | 1997-10-09 | 1999-04-27 | Seiko Epson Corp | Cardiac function diagnostic device |
| CN1813628A (en) * | 2006-01-18 | 2006-08-09 | 哈尔滨工业大学 | Pulse condition detecting device capable of automatically regulating pulse-feeling pressure and pulse condition characteristic extracting method |
| CN107822608A (en) * | 2017-10-26 | 2018-03-23 | 中国民航大学 | Pulse wave feature extracting method based on gauss hybrid models |
| CN115089133A (en) * | 2022-06-20 | 2022-09-23 | 中国农业银行股份有限公司 | Pulse condition signal identification method and device |
| CN217793021U (en) * | 2022-06-07 | 2022-11-15 | 常州金姆健康科技有限公司 | Pulse signal detects draws pulse diagnosis appearance |
| CN116502124A (en) * | 2023-04-27 | 2023-07-28 | 芯海科技（深圳）股份有限公司 | Pulse wave signal quality evaluation method and device and computer readable storage medium |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8491487B2 (en) * | 2009-02-11 | 2013-07-23 | Edwards Lifesciences Corporation | Detection of parameters in cardiac output related waveforms |
| US10765331B2 (en) * | 2015-04-02 | 2020-09-08 | Microsoft Technology Licensing, Llc | Wearable pulse sensing device signal quality estimation |
| CN110123270A (en) * | 2019-04-02 | 2019-08-16 | 博脉有限公司 | Pulse condition analyser output signal modification method and system |
-
2024
- 2024-01-09 CN CN202410028486.3A patent/CN117547228B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07148126A (en) * | 1993-11-30 | 1995-06-13 | Seiko Epson Corp | Pulse wave analyzing apparatus |
| JPH10248819A (en) * | 1997-03-13 | 1998-09-22 | Seiko Epson Corp | Pulse-wave diagnostic device |
| JPH11113863A (en) * | 1997-10-09 | 1999-04-27 | Seiko Epson Corp | Cardiac function diagnostic device |
| CN1813628A (en) * | 2006-01-18 | 2006-08-09 | 哈尔滨工业大学 | Pulse condition detecting device capable of automatically regulating pulse-feeling pressure and pulse condition characteristic extracting method |
| CN107822608A (en) * | 2017-10-26 | 2018-03-23 | 中国民航大学 | Pulse wave feature extracting method based on gauss hybrid models |
| CN217793021U (en) * | 2022-06-07 | 2022-11-15 | 常州金姆健康科技有限公司 | Pulse signal detects draws pulse diagnosis appearance |
| CN115089133A (en) * | 2022-06-20 | 2022-09-23 | 中国农业银行股份有限公司 | Pulse condition signal identification method and device |
| CN116502124A (en) * | 2023-04-27 | 2023-07-28 | 芯海科技（深圳）股份有限公司 | Pulse wave signal quality evaluation method and device and computer readable storage medium |
Non-Patent Citations (2)
| Title |
|---|
| Song Zha ; Zhaofeng Wu ; Shiqi Xing ; Jingyi Xie.A High-Performance Energy Selective Surface Design Method based on Series Resonant.《2022 IEEE MTT-S International Microwave Workshop Series on Advanced Materials and Processes for RF and THz Applications (IMWS-AMP)》.2022,第1-3页. * |
| 基于光电容积脉搏波的呼吸信号和呼吸频率提取算法研究;肖胜朗;《中国优秀硕士学位论文全文数据库 (基础科学辑)》;20220415(第4期期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117547228A (en) | 2024-02-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101449973B (en) | Judgment index generation method and device for cardiac interference signal identification | |
| US20180116597A1 (en) | Physiological sign information acquisition method and system | |
| CN109009124A (en) | Respiratory rate measurement and object localization method based on ULTRA-WIDEBAND RADAR | |
| CN112998690B (en) | Pulse wave multi-feature fusion-based respiration rate extraction method | |
| CN111973151A (en) | Infectious disease monitoring system and method based on wearable intelligent bandage | |
| CN113848544A (en) | Human body existence induction detection method and device based on Doppler radar and storage medium | |
| CN110960211A (en) | Embedded-based real-time electrocardio monitoring system | |
| CN101382837B (en) | Computer mouse control device of compound motion mode | |
| CN117547228B (en) | Data quality analysis method of biological pulse diagnosis signal | |
| CN103315732B (en) | Specific acupoint electrograph meter for selecting human body acupoint bio-electricity signal | |
| Wan et al. | Combining parallel adaptive filtering and wavelet threshold denoising for photoplethysmography-based pulse rate monitoring during intensive physical exercise | |
| Xu et al. | Digital filter design for peak detection of surface EMG | |
| CN108392176A (en) | A kind of Sleep architecture detection method based on the acquisition of heart impact signal | |
| CN1249478A (en) | Aid system for doctors in same community | |
| Liu et al. | Ultrasound-based 3-D gesture recognition: Signal optimization, trajectory, and feature classification | |
| Karilainen et al. | Mobile patient monitoring based on impedance-loaded SAW-sensors | |
| CN107811619B (en) | Portable pulse diagnosis instrument and analysis method thereof | |
| CN105919568A (en) | Gabor transformation based method and device for extracting and analyzing breathing and heartbeat signals | |
| CN105411568A (en) | Automatic medical monitoring method | |
| CN105105728B (en) | Pulse wave measuring method and device | |
| CN205849445U (en) | A kind of human body physiological characteristics information collecting device merging multisensor | |
| Zhang et al. | Research on gesture definition and electrode placement in pattern recognition of hand gesture action SEMG | |
| CN203354541U (en) | Acupoint electronic graph meter used for non-injury human body visceral organ disease diagnosing and treatment | |
| Xu et al. | Wearable psychological stress monitoring equipment and data analysis based on a wireless sensor | |
| CN111000540A (en) | Student physique health detection system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |