CN106175739B - Phase measurement and optimization method and system between a kind of pulse R R during sleep - Google Patents
Phase measurement and optimization method and system between a kind of pulse R R during sleep Download PDFInfo
- Publication number
- CN106175739B CN106175739B CN201610630794.9A CN201610630794A CN106175739B CN 106175739 B CN106175739 B CN 106175739B CN 201610630794 A CN201610630794 A CN 201610630794A CN 106175739 B CN106175739 B CN 106175739B
- Authority
- CN
- China
- Prior art keywords
- pulse
- interval
- value
- time
- sleep
- 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
Classifications
-
- 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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/024—Measuring pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Measuring devices for evaluating the respiratory organs
- A61B5/0826—Detecting or evaluating apnoea events
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4809—Sleep detection, i.e. determining whether a subject is asleep or not
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4812—Detecting sleep stages or cycles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
- A61B5/4818—Sleep apnoea
-
- 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/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0247—Pressure sensors
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Pulmonology (AREA)
- Anesthesiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种适用于日常生活环境下的有关心率变异性的R-R间期优化技术,特别涉及一种以非接触性的方式监测睡眠中的呼吸、脉搏等生命体征实现对R-R间期优化的睡眠期间的脉搏R-R间期测量和优化方法及系统。The present invention relates to an R-R interval optimization technology applicable to heart rate variability in daily life environment, in particular to a non-contact way to monitor breathing, pulse and other vital signs during sleep to optimize the R-R interval Pulse R-R interval measurement and optimization method and system during sleep.
背景技术Background technique
据《中国心血管病报告(2015年)》估计目前中国有心血管病患者2.9亿,每年有350万人死于心血管疾病。随着生活水平的提高、生活节奏的加快及全球老年化社会的到来,心脏病的发病率和死亡率将愈发上升,由此可见,心脏病已经成为威胁人类的主要疾病之一。现代医学研究表明,获得连续且准确的R-R间期进行心率变异性的分析能够快速了解心脏的功能状态。According to the "China Cardiovascular Disease Report (2015)", it is estimated that there are currently 290 million cardiovascular disease patients in China, and 3.5 million people die of cardiovascular disease every year. With the improvement of living standards, the acceleration of the pace of life and the arrival of the global aging society, the morbidity and mortality of heart disease will increase. It can be seen that heart disease has become one of the main diseases threatening human beings. Modern medical research has shown that obtaining continuous and accurate R-R intervals for heart rate variability analysis can quickly understand the functional state of the heart.
为得到准确的R-R间期值,公认的最标准的方法是医院使用的睡眠多导图(PSG),但PSG由于流程复杂、费用昂贵、用户体验度差等原因而难以普及应用。因此,研究替代PSG、且方便、准确率高的检测方法将为R-R间期的优化与心率变异性的分析提供有力的技术支撑,具有较高的实用价值。In order to obtain accurate R-R interval values, the most recognized method is polysomnography (PSG) used in hospitals, but PSG is difficult to popularize and apply due to the complexity of the process, high cost, and poor user experience. Therefore, the study of a convenient and accurate detection method to replace PSG will provide strong technical support for the optimization of R-R interval and the analysis of heart rate variability, and has high practical value.
目前获取R-R间期的方法大致有以下几种:At present, there are roughly the following methods for obtaining the R-R interval:
1、专利文献CN104398257A公开了一种心电波形周期性统计分析方法,读取给定的若干个周期的心电波形,判别每个周期心电波QRS波群的位置并计算所有周期心电波形的叠加平均图,从而确定R-R间期值。此类方法需要先用PSG测出心电信号,因此,其涉及的流程较为复杂、费用昂贵、用户体验度差等问题,其次,此类方法只能对一段数据做离线分析,不能做到连续实时检测。1. Patent document CN104398257A discloses a periodic statistical analysis method of electrocardiogram waveforms, which reads a number of given cycles of electrocardiogram waveforms, distinguishes the position of each cycle of electrocardiogram QRS complexes and calculates the position of all cycle electrocardiogram waveforms. Overlay the mean plots to determine the R-R interval value. This type of method needs to use PSG to measure the ECG signal first. Therefore, the process involved is more complicated, the cost is expensive, and the user experience is poor. Secondly, this type of method can only do offline analysis of a piece of data, and cannot achieve continuous Real-time detection.
2、中国海洋大学的一份文献公开了一种通过检测R波从而确定R-R间期的方法,从MIT-BIH心电信号数据库中获得数据并对原始信号进行滤波等预处理,然后使用阈值检测算法检测R波,预测后续R-R间期值。此类方法也只能做离线分析,无法获取实时的ECG信号,对R-R间期值的计算也只能停留在理论方面,实用性差。2. A document from Ocean University of China discloses a method for determining the R-R interval by detecting R waves, obtaining data from the MIT-BIH ECG signal database and performing preprocessing such as filtering on the original signal, and then using threshold detection The algorithm detects R waves and predicts subsequent R-R interval values. This type of method can only be analyzed offline, and cannot obtain real-time ECG signals, and the calculation of the R-R interval value can only stay in theory, and has poor practicability.
3、专利文献CN103479349A提供了一种心电信号数据获取及处理方法和系统,利用R-R间期处理单元使得R-R间期的信号带宽落在希尔伯特变换最佳幅频特性所要求的频率带宽之内,并通过重采样、滤波等解决漏检、误检等问题。此类方法对原始ECG信号以及R-R间期的计算过于干预,人为修正过度,计算得出的R-R间期值的准确性有待考量。3. Patent document CN103479349A provides a method and system for acquiring and processing ECG signal data, using the R-R interval processing unit to make the signal bandwidth of the R-R interval fall within the frequency bandwidth required by the Hilbert transform for the best amplitude-frequency characteristics Within, and through resampling, filtering, etc. to solve missed detection, false detection and other problems. Such methods interfere too much with the calculation of the original ECG signal and the R-R interval, and the artificial correction is excessive, and the accuracy of the calculated R-R interval value needs to be considered.
总之,现有R-R间期值的优化方法,要么需要依赖复杂繁琐且用户体验度极低的PSG,在医院环境下才能实施检测排查,无法应用于家庭或差旅等日常生活环境;要么无法做到实时采集与检测,且计算出的R-R间期值的准确度也不高。In short, the existing optimization methods for R-R interval value either need to rely on PSG, which is complex and cumbersome and has extremely low user experience, and can only be detected and checked in a hospital environment, which cannot be applied to daily life environments such as family or business trips; or cannot be done. to real-time collection and detection, and the accuracy of the calculated R-R interval value is not high.
发明内容Contents of the invention
本发明的目的是为了解决上述现有技术的缺点和不足,提供一种睡眠期间的脉搏R-R间期测量和优化方法及系统,通过非接触性且便携度高的设备及检测准确度高的方法对受测者的睡眠情况进行实时检测,且仅需要通过检测脉搏信号,运用多种算法实现准确的R-R间期值的提取,通过脉搏峰值的协方差匹配对脉搏峰值进行初次筛选后得到初步R-R间期值,再通过AR模型对奇异的初步R-R间期值进行二次排查,达到了现有依靠PSG采集出的R-R间期值的效果,并能够满足用户处于居家或差旅等日常生活方式下进行R-R间期的优化的需求,并减少时间和金钱消耗。The purpose of the present invention is to solve the shortcomings and deficiencies of the above-mentioned prior art, and to provide a method and system for measuring and optimizing the pulse R-R interval during sleep, through non-contact and highly portable equipment and a method with high detection accuracy Real-time detection of the subject's sleep condition, and only need to detect the pulse signal, use a variety of algorithms to achieve accurate R-R interval value extraction, and obtain the preliminary R-R interval value after initial screening of the pulse peak value through the covariance matching of the pulse peak value Interval value, and then use the AR model to conduct a secondary investigation on the singular preliminary R-R interval value, which achieves the effect of the existing R-R interval value collected by PSG, and can meet the needs of users in daily life such as home or travel. Under the demand of R-R interval optimization, and reduce time and money consumption.
为了实现上述目的,本发明采用的技术方案如下:In order to achieve the above object, the technical scheme adopted in the present invention is as follows:
首先,本发明提供了一种睡眠期间的脉搏R-R间期测量和优化方法,其包括以下步骤:At first, the present invention provides a kind of pulse R-R interval measurement and optimization method during sleep, and it comprises the following steps:
步骤1:获取并计算出准确的脉率,定位脉搏峰值点,准确获取脉搏波形;Step 1: Obtain and calculate the accurate pulse rate, locate the pulse peak point, and accurately obtain the pulse waveform;
步骤2:根据脉搏波形计算得到初步R-R间期值;Step 2: Calculate the preliminary R-R interval value according to the pulse waveform;
步骤3:设定一时间间期,根据该时间间期内的初步R-R间期值,建立用于表示该时间间期内的有效峰值信息的AR模型;Step 3: Set a time interval, and establish an AR model used to represent the effective peak information in the time interval according to the preliminary R-R interval value in the time interval;
步骤4:所述AR模型根据所述时间间期内的前N个点的初步R-R间期值计算得到第N+1个点的R-R间期值;Step 4: The AR model calculates the R-R interval value of the N+1th point according to the preliminary R-R interval value of the first N points in the time interval;
步骤5:计算出由所述AR模型计算得到的第N+1个点的R-R间期值与所述时间间期内第N+1个点的实际R-R间期值之间的差值;判断该差值是否属于一阈值范围内,如果是,则记录所述时间间期内第N+1个点的实际R-R间期值;否则,视为失败检测,丢弃所述时间间期内第N+1个点的实际R-R间期值;Step 5: Calculate the difference between the R-R interval value of the N+1th point calculated by the AR model and the actual R-R interval value of the N+1th point in the time interval; judge Whether the difference falls within a threshold range, if yes, record the actual R-R interval value of the N+1 point in the time interval; otherwise, consider the detection failure and discard the Nth point in the time interval +1 point of actual R-R interval value;
步骤6:等时间间隔依次抽取与所述时间间期时长相同的多个时间间期,重复所述步骤3~5;直至所有时间间期处理完毕,实现对睡眠期间R-R间期的优化。Step 6: Sequentially extract a plurality of time intervals with the same duration as the time interval at equal time intervals, and repeat steps 3 to 5; until all time intervals are processed, the R-R interval during sleep is optimized.
作为本发明的进一步改进,所述步骤1中,具体包括以下步骤:As a further improvement of the present invention, the step 1 specifically includes the following steps:
步骤11:实时采集睡眠信号;Step 11: collect sleep signals in real time;
步骤12:对睡眠信号进行处理得到脉率;Step 12: Process the sleep signal to obtain the pulse rate;
步骤13:对所述脉率依次进行形态滤波、峰值检测和协方差匹配处理,计算得到准确的脉率数据,并根据所述脉率数据生成脉搏波形。Step 13: Perform morphological filtering, peak detection, and covariance matching processing on the pulse rate in sequence to obtain accurate pulse rate data, and generate a pulse waveform according to the pulse rate data.
作为本发明的进一步改进,所述步骤11是通过压电传感器实现对睡眠信号的实时采集,采集时该压电传感器放置于受测者枕头底面。As a further improvement of the present invention, the step 11 is to realize the real-time collection of sleep signals through a piezoelectric sensor, which is placed on the bottom of the subject's pillow during collection.
作为本发明的进一步改进,所述步骤11中压电传感器采集的睡眠信号的类型为模拟信号类型,并将该模拟信号发送到一A/D转换模块,由A/D转换模块转换成数字信号类型的睡眠信号。As a further improvement of the present invention, the type of the sleep signal collected by the piezoelectric sensor in the step 11 is an analog signal type, and the analog signal is sent to an A/D conversion module, which is converted into a digital signal by the A/D conversion module Types of sleep signals.
作为本发明的进一步改进,所述步骤2中,具体包括以下步骤:As a further improvement of the present invention, the step 2 specifically includes the following steps:
步骤21:抽取脉搏波形的一段脉搏信号,并于该段脉搏信号中建立一脉搏波形脉冲;Step 21: Extract a segment of the pulse signal from the pulse waveform, and create a pulse waveform pulse in the segment of the pulse signal;
步骤22:将脉搏波形脉冲与该段脉搏信号进行移动匹配,得到一组协方差匹配的相关函数;Step 22: moving and matching the pulse waveform pulse with the pulse signal to obtain a set of covariance matching correlation functions;
步骤23:于所述相关函数中建立一周期为6s的移动窗口,并在所述相关函数中以1s为周期移动所述移动窗口,并选取移动窗口每次移动后其内的最大数值点,并将所有选取出来的最大数值点作为脉搏信号中的脉搏波动脉冲;Step 23: Establishing a moving window with a period of 6s in the correlation function, and moving the moving window with a period of 1s in the correlation function, and selecting the maximum value point in the moving window after each movement, And all the selected maximum value points are used as the pulse fluctuation pulse in the pulse signal;
步骤24:根据心率在短时间内不具备突变的特性,对步骤23中形成的脉搏波动脉冲的脉搏周期进行判断,去除异常的峰值,得到有效峰值数据;Step 24: According to the characteristic that the heart rate does not have a sudden change in a short period of time, judge the pulse period of the pulse fluctuation pulse formed in step 23, remove the abnormal peak value, and obtain effective peak value data;
步骤25:于相关函数中将得到的有效峰值数据标记为脉搏峰值点,并计算得到初步R-R间期值。Step 25: mark the obtained effective peak data as the pulse peak point in the correlation function, and calculate the preliminary R-R interval value.
作为本发明的进一步改进,所述步骤3中,设定6s的时间间期,根据该时间间期内的4~6个初步R-R间期值Xi,建立用于表示该时间间期内的有效峰值信息的AR模型,其中i的值为1~6之间的整数,包括1和6。As a further improvement of the present invention, in the step 3, a time interval of 6s is set, and according to 4 to 6 preliminary R-R interval values Xi within the time interval, an effective AR model of peak information, where the value of i is an integer between 1 and 6, including 1 and 6.
所述步骤4中,所述AR模型根据所述时间间期内的第1~3个点的初步R-R间期值X1、X2和X3估计得到第4个点的R-R间期值X4,X4=a×X3+a2×X2+a3×X1;其中a的值为方程a+a2+a3=1的实数解。In the step 4, the AR model estimates the RR interval value X of the fourth point according to the preliminary RR interval values X 1 , X 2 and X 3 of the first to third points in the time interval 4 , X 4 =a×X 3 +a 2 ×X 2 +a 3 ×X 1 ; where the value of a is the real number solution of the equation a+a 2 +a 3 =1.
所述步骤5中,具体包括以下步骤:In the step 5, the following steps are specifically included:
步骤51:计算出由所述AR模型计算得到的第4个点的R-R间期值X4与所述时间间期内第4个点的实际R-R间期值Y4之间的差值;Step 51: Calculate the difference between the RR interval value X4 of the fourth point calculated by the AR model and the actual RR interval value Y4 of the fourth point within the time interval;
步骤52:判断该差值的绝对值|Y4-X4|是否在X4×20%之内,如果是,则记录所述时间间期内第4个点的实际R-R间期值Y4;否则,视为失败检测,丢弃所述时间间期内第4个点的实际R-R间期值Y4。Step 52: Determine whether the absolute value |Y 4 -X 4 | of the difference is within X 4 ×20%, if yes, record the actual RR interval value Y 4 of the fourth point in the time interval ; Otherwise, the detection is regarded as failure, and the actual RR interval value Y 4 of the fourth point in the time interval is discarded.
作为本发明的进一步改进,所述步骤6中,在所述步骤3中的6s时间间期后,以1s为时间间隔,沿时间轴正向方向依次抽取多个6s的时间间期,重复所述步骤3~5;直至所有时间间期处理完毕,实现对睡眠期间R-R间期的优化。As a further improvement of the present invention, in the step 6, after the 6s time interval in the step 3, take 1s as the time interval, sequentially extract a plurality of 6s time intervals along the forward direction of the time axis, and repeat the Steps 3-5 above; until all time intervals are processed, the optimization of the R-R interval during sleep is realized.
为达到本发明的另一目的,本发明还提供了一种与上述睡眠期间的脉搏R-R间期测量和优化方法对应的睡眠期间的脉搏R-R间期测量和优化系统,该系统包括依次连接的压电传感器、A/D转换模块和处理器;In order to achieve another purpose of the present invention, the present invention also provides a pulse R-R interval measurement and optimization system during sleep corresponding to the above-mentioned pulse R-R interval measurement and optimization method during sleep. Electrical sensors, A/D conversion modules and processors;
所述压电传感器用于实时采集受测者睡眠期间的睡眠模拟信号,并传送至所述A/D转换模块;The piezoelectric sensor is used to collect the sleep analog signal during the sleep of the subject in real time, and transmit it to the A/D conversion module;
所述A/D转换模块用于将所述睡眠模拟信号转换成睡眠数字信号,并传送至所述处理器;The A/D conversion module is used to convert the sleep analog signal into a sleep digital signal and send it to the processor;
所述处理器用于对所述睡眠数字信号进行处理得到脉率波形,并根据脉率波形进行处理得到初步R-R间期值,并对初步R-R间期值进行异点排查得到优化后的R-R间期值,实现对睡眠期间R-R间期的优化。The processor is used to process the sleep digital signal to obtain a pulse rate waveform, and to process the pulse rate waveform to obtain a preliminary R-R interval value, and to check the initial R-R interval value to obtain an optimized R-R interval value to optimize the R-R interval during sleep.
通过上述技术方案,本发明达到了以下有益的技术效果:Through the above technical scheme, the present invention has achieved the following beneficial technical effects:
本发明通过非接触性且便携度高的设备及检测准确度高的方法对受测者的睡眠情况进行实时检测,且仅需要通过检测脉搏信号,运用多种算法实现准确的R-R间期值的提取,通过脉搏峰值的协方差匹配对脉搏峰值进行初次筛选后得到初步R-R间期值,再通过AR模型对奇异的初步R-R间期值进行二次排查,达到了现有依靠PSG采集出的R-R间期值的效果,准确性高,并能够满足用户处于居家或差旅等日常生活方式下进行R-R间期的优化的需求,并减少时间和金钱消耗。The present invention detects the sleep condition of the subject in real time through a non-contact and highly portable device and a method with high detection accuracy, and only needs to detect the pulse signal and use various algorithms to realize the accurate R-R interval value. Extraction, the primary R-R interval value is obtained after the initial screening of the pulse peak value through the covariance matching of the pulse peak value, and then the singular preliminary R-R interval value is checked twice through the AR model, reaching the existing R-R interval value collected by PSG The effect of the interval value has high accuracy, and can meet the needs of users to optimize the R-R interval in daily life such as home or business trips, and reduce time and money consumption.
为了更好地理解和实施,下面结合附图详细说明本发明。For better understanding and implementation, the present invention will be described in detail below in conjunction with the accompanying drawings.
附图说明Description of drawings
图1是本发明睡眠期间的脉搏R-R间期测量和优化方法的简要流程图;Fig. 1 is a brief flow chart of pulse R-R interval measurement and optimization method during sleep of the present invention;
图2是本发明睡眠期间的脉搏R-R间期测量和优化系统的结构图;Fig. 2 is the structural diagram of pulse R-R interval measurement and optimization system during sleep of the present invention;
图3是图2中的睡眠期间的脉搏R-R间期测量和优化系统进一步改进后的结构图。FIG. 3 is a further improved structure diagram of the pulse R-R interval measurement and optimization system during sleep in FIG. 2 .
具体实施方式Detailed ways
请参阅图1,本发明提供了一种睡眠期间的脉搏R-R间期测量和优化方法,Please refer to Fig. 1, the present invention provides a kind of pulse R-R interval measurement and optimization method during sleep,
首先,本发明提供了一种睡眠期间的脉搏R-R间期测量和优化方法,其包括以下步骤:At first, the present invention provides a kind of pulse R-R interval measurement and optimization method during sleep, and it comprises the following steps:
步骤1:获取并计算出准确的脉率,定位脉搏峰值点,准确获取脉搏波形;具体地,所述步骤1包括以下步骤:Step 1: Obtain and calculate an accurate pulse rate, locate the pulse peak point, and accurately obtain the pulse waveform; specifically, the step 1 includes the following steps:
步骤11:实时采集睡眠信号;在该步骤11中,是通过压电传感器实现对睡眠信号的实时采集,采集时该压电传感器放置于受测者枕头底面;且压电传感器采集的睡眠信号的类型为模拟信号类型,其将该模拟信号发送到一A/D转换模块,由A/D转换模块转换成数字信号类型的睡眠信号。Step 11: Real-time collection of sleep signals; in this step 11, the real-time collection of sleep signals is realized by a piezoelectric sensor, which is placed on the bottom of the subject's pillow during collection; and the sleep signal collected by the piezoelectric sensor The type is an analog signal type, which sends the analog signal to an A/D conversion module, and is converted into a digital signal type sleep signal by the A/D conversion module.
步骤12:对睡眠信号进行处理得到脉率;具体地,步骤12是对数字类型的睡眠信号进行处理得到脉率;Step 12: Process the sleep signal to obtain the pulse rate; specifically, step 12 is to process the digital sleep signal to obtain the pulse rate;
步骤13:对所述脉率依次进行形态滤波、峰值检测和协方差匹配处理,计算得到准确的脉率数据,并根据所述脉率数据生成脉搏波形。Step 13: Perform morphological filtering, peak detection, and covariance matching processing on the pulse rate in sequence to obtain accurate pulse rate data, and generate a pulse waveform according to the pulse rate data.
步骤2:根据脉搏波形计算得到初步R-R间期值;具体地,所述步骤2包括以下步骤:Step 2: Calculate the preliminary R-R interval value according to the pulse waveform; specifically, the step 2 includes the following steps:
步骤21:抽取脉搏波形的一段脉搏信号,并于该段脉搏信号中建立一脉搏波形脉冲;Step 21: Extract a segment of the pulse signal from the pulse waveform, and create a pulse waveform pulse in the segment of the pulse signal;
步骤22:将脉搏波形脉冲与该段脉搏信号进行移动匹配,得到一组协方差匹配的相关函数;Step 22: moving and matching the pulse waveform pulse with the pulse signal to obtain a set of covariance matching correlation functions;
步骤23:于所述相关函数中建立一周期为6s的移动窗口,并在所述相关函数中以1s为周期移动所述移动窗口,并选取移动窗口每次移动后其内的最大数值点,并将所有选取出来的最大数值点作为脉搏信号中的脉搏波动脉冲;Step 23: Establishing a moving window with a period of 6s in the correlation function, and moving the moving window with a period of 1s in the correlation function, and selecting the maximum value point in the moving window after each movement, And all the selected maximum value points are used as the pulse fluctuation pulse in the pulse signal;
步骤24:根据心率在短时间内不具备突变的特性,对步骤23中形成的脉搏波动脉冲的脉搏周期进行判断,去除异常的峰值,得到有效峰值数据;在本实施例中,该异常的峰值的选取是以其前后两点为基准,如果其与前后两点的差距大,突变大,则认为该峰值点异常,并将其排除。Step 24: According to the characteristic that the heart rate does not have a sudden change in a short period of time, judge the pulse period of the pulse fluctuation pulse formed in step 23, remove the abnormal peak value, and obtain effective peak value data; in this embodiment, the abnormal peak value The selection of the peak point is based on the two points before and after it. If there is a large gap between it and the two points before and after, and the mutation is large, the peak point is considered abnormal and excluded.
步骤25:于相关函数中将得到的有效峰值数据标记为脉搏峰值点,并计算得到初步R-R间期值。Step 25: mark the obtained effective peak data as the pulse peak point in the correlation function, and calculate the preliminary R-R interval value.
步骤3:设定一时间间期,根据该时间间期内的初步R-R间期值,建立用于表示该时间间期内的有效峰值信息的AR模型;在本实施例中,所述步骤3中,设定6s的时间间期,根据该时间间期内的4~6个初步R-R间期值Xi,建立用于表示该时间间期内的有效峰值信息的AR模型,其中i的值为1~6之间的整数,包括1和6。Step 3: set a time interval, according to the preliminary R-R interval value in the time interval, establish an AR model for representing the effective peak information in the time interval; in this embodiment, the step 3 In , a time interval of 6s is set, and an AR model used to represent the effective peak information in this time interval is established according to 4 to 6 preliminary R-R interval values Xi in this time interval, where the value of i is An integer between 1 and 6, including 1 and 6.
步骤4:所述AR模型根据所述时间间期内的前N个点的初步R-R间期值计算得到第N+1个点的R-R间期值;在本实施例中,所述步骤4中,所述AR模型根据所述时间间期内的第1~3个点的初步R-R间期值X1、X2和X3估计得到第4个点的R-R间期值X4,X4=a×X3+a2×X2+a3×X1;其中a的值为方程a+a2+a3=1的实数解。Step 4: The AR model calculates the RR interval value of the N+1th point according to the preliminary RR interval value of the first N points in the time interval; in this embodiment, in the step 4 , the AR model estimates the RR interval value X 4 of the fourth point based on the preliminary RR interval values X 1 , X 2 and X 3 of the first to third points in the time interval, X 4 = a×X 3 +a 2 ×X 2 +a 3 ×X 1 ; where the value of a is the real number solution of the equation a+a 2 +a 3 =1.
步骤5:计算出由所述AR模型计算得到的第N+1个点的R-R间期值与所述时间间期内第N+1个点的实际R-R间期值之间的差值;判断该差值是否属于一阈值范围内,如果是,则记录所述时间间期内第N+1个点的实际R-R间期值;否则,视为失败检测,丢弃所述时间间期内第N+1个点的实际R-R间期值;在本实施例中,所述步骤5中,具体包括以下步骤:Step 5: Calculate the difference between the R-R interval value of the N+1th point calculated by the AR model and the actual R-R interval value of the N+1th point in the time interval; judge Whether the difference falls within a threshold range, if yes, record the actual R-R interval value of the N+1 point in the time interval; otherwise, consider the detection failure and discard the Nth point in the time interval The actual R-R interval value of +1 point; In the present embodiment, in described step 5, specifically comprise the following steps:
步骤51:计算出由所述AR模型计算得到的第4个点的R-R间期值X4与所述时间间期内第4个点的实际R-R间期值Y4之间的差值;Step 51: Calculate the difference between the RR interval value X4 of the fourth point calculated by the AR model and the actual RR interval value Y4 of the fourth point within the time interval;
步骤52:判断该差值的绝对值|Y4-X4是否在X4×20%之内,如果是,则记录所述时间间期内第4个点的实际R-R间期值Y4;否则,视为失败检测,丢弃所述时间间期内第4个点的实际R-R间期值Y4。Step 52: Judging whether the absolute value |Y 4 -X 4 of the difference is within X 4 ×20%, if yes, record the actual RR interval value Y 4 of the fourth point in the time interval; Otherwise, it is regarded as a failed detection, and the actual RR interval value Y 4 of the fourth point in the time interval is discarded.
步骤6:等时间间隔依次抽取与所述时间间期时长相同的多个时间间期,重复所述步骤3~5;直至所有时间间期处理完毕,实现对睡眠期间R-R间期的优化。在本实施例中,所述步骤6中,在所述步骤3中的6s时间间期后,以1s为时间间隔,沿时间轴正向方向依次抽取多个6s的时间间期,重复所述步骤3~5;直至所有时间间期处理完毕,实现对睡眠期间R-R间期的优化。Step 6: Sequentially extract a plurality of time intervals with the same duration as the time interval at equal time intervals, and repeat steps 3 to 5; until all time intervals are processed, the R-R interval during sleep is optimized. In this embodiment, in the step 6, after the 6s time interval in the step 3, a plurality of 6s time intervals are sequentially extracted along the forward direction of the time axis with 1s as the time interval, and the described Steps 3-5: until all the time intervals are processed, the optimization of the R-R interval during sleep is realized.
请参阅图2,另外,本发明还提供了一种与上述睡眠期间的脉搏R-R间期测量和优化方法对应的睡眠期间的脉搏R-R间期测量和优化系统,该系统包括依次连接的压电传感器10、A/D转换模块20和处理器30。Please refer to Fig. 2, in addition, the present invention also provides a pulse R-R interval measurement and optimization system during sleep corresponding to the above-mentioned pulse R-R interval measurement and optimization method during sleep, the system includes piezoelectric sensors connected in sequence 10. An A/D conversion module 20 and a processor 30 .
所述压电传感器10用于实时采集受测者睡眠期间的睡眠模拟信号,并传送至所述A/D转换模块20。The piezoelectric sensor 10 is used to collect sleep analog signals of the subject during sleep in real time and transmit them to the A/D conversion module 20 .
所述A/D转换模块20用于将所述睡眠模拟信号转换成睡眠数字信号,并传送至所述处理器30。The A/D conversion module 20 is used for converting the sleep analog signal into a sleep digital signal and sending it to the processor 30 .
所述处理器30用于对所述睡眠数字信号进行处理得到脉率波形,并根据脉率波形进行处理得到初步R-R间期值,并对初步R-R间期值进行异点排查得到优化后的R-R间期值,实现对睡眠期间R-R间期的优化。The processor 30 is used to process the sleep digital signal to obtain a pulse rate waveform, and process according to the pulse rate waveform to obtain a preliminary R-R interval value, and perform anomaly investigation on the preliminary R-R interval value to obtain an optimized R-R interval value. Interval value to optimize the R-R interval during sleep.
请参阅图3,为避免压电传感器10受到环境因素影响而发生走位现象,同时延长压电传感器10的使用寿命,作为一种更优的技术方案,本发明睡眠期间呼吸暂停的判断系统还包括一传感器安装板40;所述压电传感器10嵌设安装于所述传感器安装板40的中部,且其检测端外漏设置于所述传感器安装板40的表面,并与传感器安装板40的表面平齐。在本实施例中,传感器安装板40的长宽大小优选为与受测者的枕头的长宽大小一致或相近。Please refer to Fig. 3, in order to prevent the piezoelectric sensor 10 from being affected by environmental factors and to prolong the service life of the piezoelectric sensor 10, as a better technical solution, the system for judging apnea during sleep of the present invention is also It includes a sensor mounting plate 40; the piezoelectric sensor 10 is embedded and installed in the middle of the sensor mounting plate 40, and its detection end is leakingly arranged on the surface of the sensor mounting plate 40, and is connected to the sensor mounting plate 40 The surface is even. In this embodiment, the length and width of the sensor mounting plate 40 are preferably consistent with or similar to the length and width of the subject's pillow.
为得到外部干扰小的睡眠数字信号,作为一种更优的技术方案,所述A/D转换模块20包括依次电连接并集成于同一PCB板上的滤波电路、放大电路和A/D转换电路。所述滤波电路与所述压电传感器10电连接,并对由压电传感器10传送的模拟信号进行滤波处理后输入到放大电路;所述放大电路对经滤波处理后的模拟信号进行放大,并传送至所述A/D转换电路;所述A/D转换电路通过PCB板上的串口与处理器30连接,其将依次经过滤波和放大处理的模拟信号转换成数字信号,并通过串口将数字信号传送至处理器30,由处理器30对所述数字信号进行分析处理。In order to obtain a sleep digital signal with little external interference, as a better technical solution, the A/D conversion module 20 includes a filter circuit, an amplifier circuit and an A/D conversion circuit that are electrically connected in sequence and integrated on the same PCB board . The filtering circuit is electrically connected to the piezoelectric sensor 10, and the analog signal transmitted by the piezoelectric sensor 10 is filtered and then input to the amplifying circuit; the amplifying circuit amplifies the filtered analog signal, and Delivered to the A/D conversion circuit; the A/D conversion circuit is connected with the processor 30 through the serial port on the PCB board, which converts the analog signal through filtering and amplifying successively into a digital signal, and converts the digital signal through the serial port The signal is sent to the processor 30, and the processor 30 analyzes and processes the digital signal.
在需要检测脉搏信号时,将整块传感器安装板40放置在受测者睡眠时需要用到的枕头的底面,或套入到枕头的枕巾内,和枕头一起受压;而A/D转换模块20可以放置在床边或其它方便放置的地方。待受测者睡眠时即可通过本系统检测脉搏信号,处理器30即可在受测者的睡眠期间不断处理数据,完成后续分析、处理和优化的工作,实现对睡眠期间R-R间期的优化。具体地,本系统实现对睡眠期间R-R间期的优化的工作过程和工作原理可结合本发明的睡眠期间的脉搏R-R间期测量和优化方法进行理解,故在此不再赘述。When the pulse signal needs to be detected, the entire sensor mounting plate 40 is placed on the bottom surface of the pillow that the subject needs to use when sleeping, or inserted into the pillow cover of the pillow, and pressed together with the pillow; and the A/D conversion module 20 can be placed on the bedside or other convenient places. When the subject is sleeping, the pulse signal can be detected by this system, and the processor 30 can continuously process the data during the sleep period of the subject, complete the subsequent analysis, processing and optimization work, and realize the optimization of the R-R interval during sleep . Specifically, the working process and working principle of this system to realize the optimization of the R-R interval during sleep can be understood in combination with the method for measuring and optimizing the pulse R-R interval during sleep of the present invention, so it will not be repeated here.
在其它变形实施例中,所述处理器30可以替换为计算机、手机、平板电脑、手表等其他智能终端设备。In other modified embodiments, the processor 30 may be replaced by other intelligent terminal devices such as computers, mobile phones, tablet computers, and watches.
相对于现有技术,本发明睡眠期间的脉搏R-R间期测量和优化方法及系统通过非接触性且便携度高的设备及检测准确度高的方法对受测者的睡眠情况进行实时检测,且仅需要通过检测脉搏信号,运用多种算法实现准确的R-R间期值的提取,通过脉搏峰值的协方差匹配对脉搏峰值进行初次筛选后得到初步R-R间期值,再通过AR模型对奇异的初步R-R间期值进行二次排查,达到了现有依靠PSG采集出的R-R间期值的效果,并能够满足用户处于居家或差旅等日常生活方式下进行R-R间期的优化的需求,并减少时间和金钱消耗。Compared with the prior art, the pulse R-R interval measurement and optimization method and system during sleep of the present invention detect the sleep condition of the subject in real time through non-contact and highly portable equipment and methods with high detection accuracy, and It is only necessary to detect the pulse signal and use a variety of algorithms to achieve accurate R-R interval value extraction. After the initial screening of the pulse peak value through the covariance matching of the pulse peak value, the initial R-R interval value is obtained, and then the singular preliminary R-R interval value is obtained through the AR model. The second investigation of the R-R interval value has achieved the effect of the existing R-R interval value collected by PSG, and can meet the needs of users to optimize the R-R interval in daily life such as home or travel, and reduce Time and money consuming.
本发明并不局限于上述实施方式,如果对本发明的各种改动或变形不脱离本发明的精神和范围,倘若这些改动和变形属于本发明的权利要求和等同技术范围之内,则本发明也意图包含这些改动和变形。The present invention is not limited to the above-mentioned embodiments, if the various changes or deformations of the present invention do not depart from the spirit and scope of the present invention, if these changes and deformations belong to the claims of the present invention and the equivalent technical scope, then the present invention is also It is intended that such modifications and variations are included.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610630794.9A CN106175739B (en) | 2016-08-01 | 2016-08-01 | Phase measurement and optimization method and system between a kind of pulse R R during sleep |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610630794.9A CN106175739B (en) | 2016-08-01 | 2016-08-01 | Phase measurement and optimization method and system between a kind of pulse R R during sleep |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106175739A CN106175739A (en) | 2016-12-07 |
CN106175739B true CN106175739B (en) | 2018-03-23 |
Family
ID=57498547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610630794.9A Active CN106175739B (en) | 2016-08-01 | 2016-08-01 | Phase measurement and optimization method and system between a kind of pulse R R during sleep |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106175739B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107252313A (en) * | 2017-05-25 | 2017-10-17 | 深圳市卡迪赛克科技有限公司 | The monitoring method and system of a kind of safe driving, automobile, readable storage medium storing program for executing |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7025729B2 (en) * | 2001-09-14 | 2006-04-11 | Biancamed Limited | Apparatus for detecting sleep apnea using electrocardiogram signals |
US8233973B2 (en) * | 2008-05-02 | 2012-07-31 | Spacelabs Healthcare, Llc | Methods for detection of cardiac arrhythmias |
DE102013001077A1 (en) * | 2013-01-23 | 2014-07-24 | Sissel International GmbH | pillow |
CN103584854B (en) * | 2013-11-29 | 2015-07-08 | 重庆海睿科技有限公司 | Extraction method of electrocardiosignal R waves |
US9545227B2 (en) * | 2013-12-13 | 2017-01-17 | Vital Connect, Inc. | Sleep apnea syndrome (SAS) screening using wearable devices |
CN104523262A (en) * | 2014-11-18 | 2015-04-22 | 南京丰生永康软件科技有限责任公司 | Sleep quality detection method based on electrocardiosignals |
CN105662375B (en) * | 2016-03-17 | 2019-01-29 | 广州中科新知科技有限公司 | The method and apparatus of non-contact detection vital sign parameter signals |
-
2016
- 2016-08-01 CN CN201610630794.9A patent/CN106175739B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106175739A (en) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106937808B (en) | Data acquisition system of intelligent mattress | |
US20130317377A1 (en) | Generative Model-Driven Resource-Efficient Monitoring in Body Sensor Networks | |
Nabar et al. | GeM-REM: Generative model-driven resource efficient ecg monitoring in body sensor networks | |
Castro et al. | Heart sound segmentation of pediatric auscultations using wavelet analysis | |
CN103006195B (en) | Non-contact vital sign data monitoring system and non-contact vital sign data monitoring method on basis of image processing | |
CN105051799A (en) | Method for detecting falls and a fall detector. | |
Sahu et al. | IoT-enabled cloud-based real-time remote ECG monitoring system | |
CN104188652B (en) | Electrocardiogram data quality real-time control method and system | |
CN109528159B (en) | Human sleep and respiration monitoring system and method based on bed body | |
CN106344034A (en) | Sleep quality evaluation system and method | |
CN105982642A (en) | Sleep detection method and detection system based on body vibration signals | |
CN104473628A (en) | Wearable cardiac sound monitoring device | |
CN108186000A (en) | Real-time blood pressure monitor system and method based on heart impact signal and photosignal | |
WO2016168979A1 (en) | Vital sign analysis method and system | |
CN116369888B (en) | Non-contact heart rate variability data acquisition method and device | |
CN104545863A (en) | BCG heart rate extraction method and system based on fuzzy pattern recognition | |
CN107080545A (en) | A kind of lie detection system based on brain electricity | |
WO2019000338A1 (en) | Physiological information measurement method, and physiological information monitoring apparatus and device | |
CN107348971A (en) | A kind of heart disease screening system based on heart sound detection and machine learning algorithm | |
CN110115583A (en) | The method and apparatus of monitoring of respiration | |
CN106175770A (en) | Apneic determination methods and system during a kind of sleep | |
CN204260739U (en) | Electrocardiographic quality of data real-time control system | |
CN106175739B (en) | Phase measurement and optimization method and system between a kind of pulse R R during sleep | |
CN103876731B (en) | A kind of fetal electrocardiogram signal extraction device and method | |
CN203290920U (en) | Portable emotion analysis meter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |