CN113384264A - 基于雷达的呼吸频率检测方法及睡眠监测设备 - Google Patents
基于雷达的呼吸频率检测方法及睡眠监测设备 Download PDFInfo
- Publication number
- CN113384264A CN113384264A CN202110656598.XA CN202110656598A CN113384264A CN 113384264 A CN113384264 A CN 113384264A CN 202110656598 A CN202110656598 A CN 202110656598A CN 113384264 A CN113384264 A CN 113384264A
- Authority
- CN
- China
- Prior art keywords
- target
- radar
- echo signal
- frame
- phase
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
- A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- 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/4815—Sleep quality
-
- 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/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/62—Sense-of-movement determination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Pulmonology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明适用于人体生理特征检测技术领域,提供了一种基于雷达的呼吸频率检测方法及睡眠监测设备,上述方法包括:获取雷达的回波信号;确定雷达探测范围内至少一个目标相对于雷达的方位角;针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。本发明通过方位角区分不同的目标,分别根据各个目标的方位角及雷达回波信号,分别确定各个目标的呼吸频率,可同时检测多个目标的呼吸频率,应用场景更广泛。
Description
技术领域
本发明属于人体生理特征检测技术领域,尤其涉及一种基于雷达的呼吸频率检测方法及睡眠监测设备。
背景技术
睡眠监测设备主要用于监测人体睡眠时的呼吸深度、频率等,从而对人体的健康情况进行评估。现有睡眠监测设备主要分为接触式和非接触式两种。非接触式主要应用微波雷达传感器,通过测量人体呼吸时的相位变化分析呼吸频率。非接触式睡眠监测设备相对于接触式睡眠监测设备不依赖对用户的接触效果,在睡眠监测领域得到了广泛的应用。
现有技术中,非接触式睡眠监测设备只能检测单一用户的呼吸频率,无法同时监测两个或多个用户的呼吸频率,应用场景受限。
发明内容
有鉴于此,本发明实施例提供了一种基于雷达的呼吸频率检测方法及睡眠监测设备,以解决现有技术中非接触式睡眠监测设备无法检测多用户的呼吸频率的问题。
本发明实施例的第一方面提供了一种基于雷达的呼吸频率检测方法,包括:
获取雷达的回波信号;
确定雷达探测范围内至少一个目标相对于雷达的方位角;
针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。
本发明实施例的第二方面提供了一种睡眠监测设备,包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序,处理器执行计算机程序时实现如本发明实施例第一方面提供的基于雷达的呼吸频率检测方法的步骤。
本发明实施例的第三方面提供了一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现如本发明实施例第一方面提供的基于雷达的呼吸频率检测方法的步骤。
本发明实施例提供了一种基于雷达的呼吸频率检测方法及睡眠监测设备,上述方法包括:获取雷达的回波信号;确定雷达探测范围内至少一个目标相对于雷达的方位角;针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。本发明实施例通过方位角区分不同的目标,分别根据各个目标的方位角及雷达回波信号,分别确定各个目标的呼吸频率,可同时检测多个目标的呼吸频率,应用场景更广泛。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于雷达的呼吸频率检测方法的实现流程示意图;
图2是本发明实施例提供的雷达的发射天线及接收天线的示意图;
图3是本发明实施例提供的雷达发射信号的示意图;
图4是本发明实施例提供的雷达架设示意图;
图5是本发明实施例提供的基于雷达的呼吸频率检测装置的示意图;
图6是本发明实施例提供的睡眠监测设备的示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
为了说明本发明的技术方案,下面通过具体实施例来进行说明。
随着人们越来越注意自身健康状态,睡眠质量监测成为人们关注的重点需求之一。睡眠监测设备主要用于监测人体睡眠时的呼吸深度、频率等参数,根据对呼吸深度及频率的分析,得到睡眠深度、睡眠时长、身体健康状态等高阶信息。
现有睡眠质量监测设备主要分为接触式和非接触式两种;
接触式包括手环、睡眠监测枕、睡眠监测床垫等。这类设备对用户使用习惯和接触效果依赖程度较高,一旦用户没有佩戴、或是接触效果不好,其监测结果就会大打折扣,甚至无法获得监测结果。
非接触式主要指微波雷达传感器。现有的雷达传感器通常放置在床头或者天花板上,通过测量人体呼吸时的相位变化,分析呼吸频率。但现有技术中无法同时测量床上的两个人的呼吸频率,使得使用场景大大受限。
基于以上问题,参考图1,本发明实施例提供了一种基于雷达的呼吸频率检测方法,包括:
S101:获取雷达的回波信号;
S102:确定雷达探测范围内至少一个目标相对于雷达的方位角;
S103:针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率
雷达的发射天线按帧发射信号,该发射信号可以为FMCW(Frequency ModulatedContinuous Wave,调频连续波)信号,参考图3。其中,T为调频时长,B为调频带宽,Tframe为两帧信号之间的间隔。Nr个接收天线分别接收回波信号,并进行混频、滤波、AD采样等处理,得到Nr个接收天线的数字回波信号;对Nr个接收天线的数字回波信号采样,将一帧内Nr个接收通道的数字回波信号记为S(i_NS,i_Nr),其中,i_Nr=1,2,…,Nr,i_NS=1,2,…,N,N为采样点数,N=fs×T,fs为采样频率。其中,雷达的回波信号可以为预设数量帧的回波信号。雷达的回波信号包含目标的呼吸频率特征,因此首先获取雷达的回波信号。其中,雷达可以为毫米波雷达。
导向矢量是阵列天线的所有阵元对具有单位能量窄带信源的响应,本质是描述空间相位差的,反应目标的方位信息。本发明实施例中确定目标的方位角,根据目标的方位角确定该目标对应的导向矢量。
例如,参考图4,雷达架设在床体正上方(例如,天花板上),目标为躺在床上的两人(目标1和目标2),两人相对于雷达的位置不同,分别确定两人相对于雷达的方位角,进而分别确定两人的导向矢量,再根据导向矢量及雷达的回波信号确定携带该目标特征信息的目标频谱。频谱是频率谱密度的简称,频谱峰值处对应的频率的信号强度最大,因此,该目标频谱峰值处的频率即可确定为目标的呼吸频率。
本发明实施例中,通过方位角区分不同的目标,分别根据各个目标的方位角及雷达回波信号,得到各个携带目标特征信息的目标频谱,进而根据频谱特性确定各个目标的呼吸频率,可同时检测多个目标的呼吸频率,应用场景更广泛。
一些实施例中,雷达的回波信号包括:预设数量帧回波信号;S103可以包括:
S1031:针对每帧回波信号,根据该目标对应的导向矢量对该帧回波信号进行波束形成,得到该帧回波信号的目标波束;根据该帧回波信号的目标波束确定该帧回波信号的目标相位;
S1032:将各帧回波信号的目标相位按时间顺序排列得到该目标对应的目标序列,并对该目标对应的目标序列进行傅里叶变换,得到该目标对应的目标频谱。
波束形成的目的是使多阵元构成的基阵经适当处理得到在预定方向的指向性。也即使回波信号具有一定的指向性,获取回波信号中携带的该目标的信息。进而根据目标波束确定目标相位,经过傅里叶变换得到目标频谱。其中,预设数量可根据实际应用需求设置。
基于以上,雷达的回波信号包括:预设数量帧回波信号,可每预设数量帧根据上述方法确定一次目标的呼吸频率。
其中,该帧回波信号的目标波束Sz的计算公式可以为:
Sz=S·A'
其中,A'该目标对应的导向矢量A的转置。
一些实施例中,S1031可以包括:
步骤1:对该帧回波信号的目标波束进行傅里叶变换,得到该帧回波信号的目标波束的一维距离像;
步骤2:获取该目标与雷达之间的垂直距离,并根据该目标与雷达之间的垂直距离确定该目标所在的距离单元;
步骤3:提取该帧回波信号的目标波束的一维距离像中该目标所在的距离单元处的相位,作为该帧回波信号的单元相位;
步骤4:根据该帧回波信号的单元相位,确定该帧回波信号的目标相位。
本发明实施例中采用傅里叶变换结合该目标所在的距离单元提取得到该帧回波信号的目标相位。
该帧回波信号的目标波束的一维距离像的计算公式可以为:Sfft=FFT(Sz);
该帧回波信号的单元相位的计算公式可以为:p=phase(Sfft(ind));其中,ind为该目标所在的距离单元。
一些实施例中,根据该帧回波信号的单元相位,确定该帧回波信号的目标相位,可以包括:
步骤1:获取该帧回波信号的上一帧回波信号的单元相位;
步骤2:将该帧回波信号的单元相位与该帧回波信号的上一帧回波信号的单元相位进行解缠,得到该帧回波信号对应的解缠后的相位;
步骤3:将该帧回波信号对应的解缠后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位;
步骤4:将该帧回波信号对应的解缠后的相位作为新的该帧回波信号的单元相位。
分别对各帧回波信号进行处理,其中,对于第一帧信号,其上一帧回波信号的单元相位可以为0,根据第一帧回波信号的单元相位和0相位计算得到第一帧回波信号的目标相位。同时,将解缠后的相位作为新的第一帧回波信号的单元相位,参与到第二帧回波信号的目标相位的计算中。
该帧回波信号对应的解缠后的相位p_unwrap的计算公式可以为:p_unwrap=unwrap(p,p_last);其中,p_last为该帧回波信号的上一帧回波信号的单元相位。
该帧回波信号的目标相位p_diff的计算公式可以为:
p_diff=p_unwrap-p_last。
一些实施例中,在将该帧回波信号的单元相位与该帧回波信号的上一帧回波信号的单元相位进行解缠,得到该帧回波信号对应的解缠后的相位之后,根据该帧回波信号的单元相位,确定该帧回波信号的目标相位还可以包括:
步骤5:对该帧回波信号对应的解缠后的相位进行滤波,得到滤波后的相位;
上述将该帧回波信号对应的解缠后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位,可以包括:
将该帧回波信号对应的滤波后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位。
本发明实施例中,为提高计算的准确度,还可对该帧回波信号对应的解缠后的相位滤波,采用滤波后的相位确定该帧回波信号的目标相位,提高了计算的准确性。其中,可采用IIR滤波器((Infinite Impulse Response,递归滤波器)进行滤波,通带可设置为0.1Hz~1Hz。
一些实施例中,该目标所在的距离单元ind的计算公式可以为:
其中,B为雷达的调频带宽,H为该目标与雷达之间的垂直距离,c为光速。
其中,该目标与雷达之间的垂直距离为目标所在水平面与雷达所在水平面之间的距离。
一些实施例中,雷达包括多个接收天线,且各个接收处于同一直线上;S102可以包括:
S1021:针对每个目标,获取该目标与雷达之间的垂直距离和水平距离,确定该目标与雷达之间的水平距离在雷达的各个接收天线所在竖直平面上的投影距离;并根据投影距离及该目标与雷达之间的垂直距离,确定该目标相对于雷达的方位角;
以雷达的轴线为0°方向;
当该目标位于雷达的左侧时,该目标相对于雷达的方位角theta1的计算公式为:
当该目标位于雷达的右侧时,该目标相对于雷达的方位角theta1的计算公式为:
其中,W为该目标与雷达之间的水平距离在雷达的各个接收天线所在竖直平面上的投影距离,H为该目标与雷达之间的垂直距离。
例如,参考图2,雷达的各个接收天线均处于同一直线上,且各个接收天线等间距排列,间隔为λ/2。目标相对于雷达的方位角为目标投影到接收天线所在竖直平面上后相对于雷达的角度。其中,目标与雷达之间的水平距离为雷达与目标投影到同一水平面上后两者之间的直线距离。为简化计算,可定义左侧目标的方位角为负,右侧目标的方位角为正,简单的区分左右目标。
一些实施例中,雷达部署于双人床中心轴线的上方,雷达的探测范围覆盖双人床;雷达包括多个接收天线,且各个接收天线处于与双人床中心轴线垂直的同一直线上;
至少一个目标为两个目标,两个目标包括第一目标和第二目标,且第一目标和第二目标分别位于双人床中心轴线的两侧;
S102可以包括:
S1022:根据雷达与双人床的垂直距离以及双人床的宽度,确定第一目标相对于雷达的方位角,以及第二目标相对于雷达的方位角。
一些实施例中,第一目标相对于雷达的方位角θ1的计算公式为:
第二目标相对于雷达的方位角θ2的计算公式为:
其中,w为双人床的宽度,h为雷达与双人床的垂直距离。
一些实施例中,该目标对应的导向矢量A的计算公式可以为:
A=[1,πsin(theta1),2πsin(theta1),…(Nr-1)πsin(theta1)]
其中,theta1为该目标相对于雷达的方位角,Nr为雷达的通道数。
应理解,上述实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
参考图5,本发明实施例还提供了一种基于雷达的呼吸频率检测装置,包括:
信号获取模块21,用于获取雷达的回波信号;
方位角确定模块22,用于确定雷达探测范围内至少一个目标相对于雷达的方位角;
呼吸频率确定模块23,用于针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。
一些实施例中,雷达的回波信号包括预设数量帧回波信号;呼吸频率确定模块23可以包括:
目标相位确定子模块231,用于针对每帧回波信号,根据该目标对应的导向矢量对该帧回波信号进行波束形成,得到该帧回波信号的目标波束;根据该帧回波信号的目标波束确定该帧回波信号的目标相位;
目标频谱输出子模块232,用于将各帧回波信号的目标相位按时间顺序排列得到该目标对应的目标序列,并对该目标对应的目标序列进行傅里叶变换,得到该目标对应的目标频谱。
一些实施例中,目标相位确定子模块231可以包括:
傅里叶变换单元2311,用于对该帧回波信号的目标波束进行傅里叶变换,得到该帧回波信号的目标波束的一维距离像;
距离单元确定单元2312,用于获取该目标与雷达之间的垂直距离,并根据该目标与雷达之间的垂直距离确定该目标所在的距离单元;
单元相位确定单元2313,用于提取该帧回波信号的目标波束的一维距离像中该目标所在的距离单元处的相位,作为该帧回波信号的单元相位;
目标相位确定单元2314,用于根据该帧回波信号的单元相位,确定该帧回波信号的目标相位。
一些实施例中,目标相位确定单元2314可以具体用于:
获取该帧回波信号的上一帧回波信号的单元相位;
将该帧回波信号的单元相位与该帧回波信号的上一帧回波信号的单元相位进行解缠,得到该帧回波信号对应的解缠后的相位;
将该帧回波信号对应的解缠后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位;
将该帧回波信号对应的解缠后的相位作为新的该帧回波信号的单元相位。
一些实施例中,目标相位确定单元2214还具体用于:
对该帧回波信号对应的解缠后的相位进行滤波,得到滤波后的相位;
将该帧回波信号对应的解缠后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位,包括:
将该帧回波信号对应的滤波后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位。
一些实施例中,该目标所在的距离单元ind的计算公式为:
其中,B为雷达的调频带宽,H为该目标与雷达之间的垂直距离,c为光速。
一些实施例中,雷达包括多个接收天线,且各个接收处于同一直线上;方位角确定模块22可以包括:
第一方位角确定子模块221,用于针对每个目标,获取该目标与雷达之间的垂直距离和水平距离,确定该目标与雷达之间的水平距离在雷达的各个接收天线所在竖直平面上的投影距离;并根据投影距离及该目标与雷达之间的垂直距离,确定该目标相对于雷达的方位角;
以雷达的轴线为0°方向;
当该目标位于雷达的左侧时,该目标相对于雷达的方位角theta1的计算公式为:
当该目标位于雷达的右侧时,该目标相对于雷达的方位角theta1的计算公式为:
其中,W为该目标与雷达之间的水平距离在雷达的各个接收天线所在竖直平面上的投影距离,H为该目标与雷达之间的垂直距离。
一些实施例中,雷达部署于双人床中心轴线的上方,雷达的探测范围覆盖双人床;雷达包括多个接收天线,且各个接收天线处于与双人床中心轴线垂直的同一直线上;
至少一个目标为两个目标,两个目标包括第一目标和第二目标,且第一目标和第二目标分别位于双人床中心轴线的两侧;
方位角确定模块22可以包括:
第二方位角确定子模块222,用于根据雷达与双人床的垂直距离以及双人床的宽度,确定第一目标相对于雷达的方位角,以及第二目标相对于雷达的方位角。
一些实施例中,第一目标相对于雷达的方位角θ1的计算公式为:
第二目标相对于雷达的方位角θ2的计算公式为:
其中,w为双人床的宽度,h为雷达与双人床的垂直距离。
一些实施例中,该目标对应的导向矢量A的计算公式为:
A=[1,πsin(theta1),2πsin(theta1),…(Nr-1)πsin(theta1)]
其中,theta1为该目标相对于雷达的方位角,Nr为雷达的通道数。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将睡眠监测设备的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述装置中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
图6是本发明一实施例提供的睡眠监测设备的示意框图。如图6所示,该实施例的睡眠监测设备4包括:一个或多个处理器40、存储器41以及存储在存储器41中并可在处理器40上运行的计算机程序42。处理器40执行计算机程序42时实现上述各个基于雷达的呼吸频率检测方法实施例中的步骤,例如执行图1所示的步骤S101至S103。或者,处理器40执行计算机程序42时实现上述基于雷达的呼吸频率检测装置实施例中各模块/单元的功能,例如实现图5所示模块21至23的功能。
示例性地,计算机程序42可以被分割成一个或多个模块/单元,一个或者多个模块/单元被存储在存储器41中,并由处理器40执行,以完成本申请。一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述计算机程序42在睡眠监测设备4中的执行过程。例如,计算机程序42可以被分割成信号获取模块21、方位角确定模块22及呼吸频率确定模块23。
信号获取模块21,用于获取雷达的回波信号;
方位角确定模块22,用于确定雷达探测范围内至少一个目标相对于雷达的方位角;
呼吸频率确定模块23,用于针对每个目标,根据该目标相对于雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。
其它模块或者单元在此不再赘述。
睡眠监测设备4包括但不仅限于处理器40、存储器41。本领域技术人员可以理解,图6仅仅是睡眠监测设备的一个示例,并不构成对睡眠监测设备4的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如睡眠监测设备4还可以包括输入设备、输出设备、网络接入设备、总线等。
处理器40可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器41可以是睡眠监测设备的内部存储单元,例如睡眠监测设备的硬盘或内存。存储器41也可以是睡眠监测设备的外部存储设备,例如睡眠监测设备上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,存储器41还可以既包括睡眠监测设备的内部存储单元也包括外部存储设备。存储器41用于存储计算机程序42以及睡眠监测设备所需的其他程序和数据。存储器41还可以用于暂时地存储已经输出或者将要输出的数据。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
在本申请所提供的实施例中,应该理解到,所揭露的睡眠监测设备和方法,可以通过其它的方式实现。例如,以上所描述的睡眠监测设备实施例仅仅是示意性的,例如,模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,计算机程序包括计算机程序代码,计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括:能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、电载波信号、电信信号以及软件分发介质等。需要说明的是,计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括是电载波信号和电信信号。
以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围,均应包含在本申请的保护范围之内。
Claims (10)
1.一种基于雷达的呼吸频率检测方法,其特征在于,包括:
获取雷达的回波信号;
确定所述雷达探测范围内至少一个目标相对于所述雷达的方位角;
针对每个目标,根据该目标相对于所述雷达的方位角确定该目标对应的导向矢量;根据该目标对应的导向矢量及所述雷达的回波信号,确定该目标对应的目标频谱;将该目标对应的目标频谱峰值处的频率作为该目标的呼吸频率。
2.如权利要求1所述的基于雷达的呼吸频率检测方法,其特征在于,所述雷达的回波信号包括预设数量帧回波信号;所述根据该目标对应的导向矢量及所述雷达的回波信号,确定该目标对应的目标频谱,包括:
针对每帧回波信号,根据该目标对应的导向矢量对该帧回波信号进行波束形成,得到该帧回波信号的目标波束;根据该帧回波信号的目标波束确定该帧回波信号的目标相位;
将各帧回波信号的目标相位按时间顺序排列得到该目标对应的目标序列,并对该目标对应的目标序列进行傅里叶变换,得到该目标对应的目标频谱。
3.如权利要求2所述的基于雷达的呼吸频率检测方法,其特征在于,所述根据该帧回波信号的目标波束确定该帧回波信号的目标相位,包括:
对该帧回波信号的目标波束进行傅里叶变换,得到该帧回波信号的目标波束的一维距离像;
获取该目标与所述雷达之间的垂直距离,并根据该目标与所述雷达之间的垂直距离确定该目标所在的距离单元;
提取该帧回波信号的目标波束的一维距离像中该目标所在的距离单元处的相位,作为该帧回波信号的单元相位;
根据该帧回波信号的单元相位,确定该帧回波信号的目标相位。
5.如权利要求3所述的基于雷达的呼吸频率检测方法,其特征在于,所述根据该帧回波信号的单元相位,确定该帧回波信号的目标相位,包括:
获取该帧回波信号的上一帧回波信号的单元相位;
将该帧回波信号的单元相位与该帧回波信号的上一帧回波信号的单元相位进行解缠,得到该帧回波信号对应的解缠后的相位;
将该帧回波信号对应的解缠后的相位减去该帧回波信号的上一帧回波信号的单元相位,得到该帧回波信号的目标相位;
将该帧回波信号对应的解缠后的相位作为新的该帧回波信号的单元相位。
6.如权利要求1所述的基于雷达的呼吸频率检测方法,其特征在于,所述雷达包括多个接收天线,且各个接收处于同一直线上;
所述确定所述雷达探测范围内至少一个目标相对于所述雷达的方位角,包括:
针对每个目标,获取该目标与所述雷达之间的垂直距离和水平距离,确定该目标与所述雷达之间的水平距离在所述雷达的各个接收天线所在竖直平面上的投影距离;并根据所述投影距离及该目标与所述雷达之间的垂直距离,确定该目标相对于所述雷达的方位角;
以所述雷达的轴线为0°方向;
当该目标位于所述雷达的左侧时,所述该目标相对于所述雷达的方位角theta1的计算公式为:
当该目标位于所述雷达的右侧时,所述该目标相对于所述雷达的方位角theta1的计算公式为:
其中,W为该目标与所述雷达之间的水平距离在所述雷达的各个接收天线所在竖直平面上的投影距离,H为该目标与所述雷达之间的垂直距离。
7.如权利要求1所述的基于雷达的呼吸频率检测方法,其特征在于,所述雷达部署于双人床中心轴线的上方,所述雷达的探测范围覆盖所述双人床;所述雷达包括多个接收天线,且各个接收天线处于与所述双人床中心轴线垂直的同一直线上;
所述至少一个目标为两个目标,所述两个目标包括第一目标和第二目标,且所述第一目标和所述第二目标分别位于所述双人床中心轴线的两侧;
所述确定所述雷达探测范围内至少一个目标相对于所述雷达的方位角,包括:
根据所述雷达与所述双人床的垂直距离以及所述双人床的宽度,确定所述第一目标相对于所述雷达的方位角,以及所述第二目标相对于所述雷达的方位角。
9.一种睡眠监测设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至8任一项所述的基于雷达的呼吸频率检测方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至8任一项所述的基于雷达的呼吸频率检测方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110656598.XA CN113384264B (zh) | 2021-06-11 | 2021-06-11 | 基于雷达的呼吸频率检测方法及睡眠监测设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110656598.XA CN113384264B (zh) | 2021-06-11 | 2021-06-11 | 基于雷达的呼吸频率检测方法及睡眠监测设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113384264A true CN113384264A (zh) | 2021-09-14 |
CN113384264B CN113384264B (zh) | 2023-06-06 |
Family
ID=77620872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110656598.XA Active CN113384264B (zh) | 2021-06-11 | 2021-06-11 | 基于雷达的呼吸频率检测方法及睡眠监测设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113384264B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113729674A (zh) * | 2021-09-24 | 2021-12-03 | 青岛海信日立空调系统有限公司 | 一种呼吸检测方法及装置 |
CN113729676A (zh) * | 2021-10-12 | 2021-12-03 | 九州云合(山东)智能科技有限公司 | 一种智能生命体征监护设备 |
CN114403848A (zh) * | 2022-01-05 | 2022-04-29 | 河南牧原智能科技有限公司 | 一种动物群体呼吸频率检测的方法、系统及存储介质 |
CN116035558A (zh) * | 2023-03-02 | 2023-05-02 | 中国科学技术大学 | 基于波束形成的抗干扰呼吸检测方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2517621A1 (en) * | 2011-04-29 | 2012-10-31 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radar apparatus for detecting multiple life-signs of a subject, a method and a computer program product |
CN108388912A (zh) * | 2018-02-08 | 2018-08-10 | 南京理工大学 | 基于多传感器特征优化算法的睡眠分期方法 |
US20180279884A1 (en) * | 2017-03-29 | 2018-10-04 | Texas Instruments Incorporated | Multi-person vital signs monitoring using millimeter wave (mm-wave) signals |
CN108652601A (zh) * | 2018-05-24 | 2018-10-16 | 厦门精益远达智能科技有限公司 | 一种基于调频连续波毫米波雷达的睡眠监控方法、装置和雷达系统 |
CN110346790A (zh) * | 2019-07-09 | 2019-10-18 | 长沙莫之比智能科技有限公司 | 一种基于毫米波雷达的非接触式生命体征监测方法、装置及系统 |
CN110464320A (zh) * | 2019-08-22 | 2019-11-19 | 四川长虹电器股份有限公司 | 多目标人体心率和呼吸频率测量系统及方法 |
KR102091974B1 (ko) * | 2019-07-25 | 2020-03-24 | 정기섭 | Fmcw 레이더를 이용한 위치 추적 및 호흡(심박) 패턴 측정 방법 |
CN111481184A (zh) * | 2020-04-24 | 2020-08-04 | 华侨大学 | 基于毫米波雷达技术的多目标呼吸心率监测方法与系统 |
CN111856455A (zh) * | 2020-09-04 | 2020-10-30 | 四川长虹电器股份有限公司 | 匹配不同雷达带宽的多目标心率呼吸测量方法及系统 |
CN111990968A (zh) * | 2020-08-07 | 2020-11-27 | 北京航空航天大学杭州创新研究院 | 一种信号监测方法、设备、存储介质及系统 |
CN112686094A (zh) * | 2020-12-03 | 2021-04-20 | 华中师范大学 | 一种基于毫米波雷达的非接触式身份识别方法及系统 |
CN112716462A (zh) * | 2020-12-09 | 2021-04-30 | 北京航空航天大学 | 一种可控照射指向的窄波束毫米波人体心跳/呼吸体征监测装置 |
CN112882016A (zh) * | 2019-11-29 | 2021-06-01 | 南京理工大学 | 基于改进的稳健自适应波束形成的多人生命体征检测方法 |
-
2021
- 2021-06-11 CN CN202110656598.XA patent/CN113384264B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2517621A1 (en) * | 2011-04-29 | 2012-10-31 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | A radar apparatus for detecting multiple life-signs of a subject, a method and a computer program product |
US20180279884A1 (en) * | 2017-03-29 | 2018-10-04 | Texas Instruments Incorporated | Multi-person vital signs monitoring using millimeter wave (mm-wave) signals |
CN108388912A (zh) * | 2018-02-08 | 2018-08-10 | 南京理工大学 | 基于多传感器特征优化算法的睡眠分期方法 |
CN108652601A (zh) * | 2018-05-24 | 2018-10-16 | 厦门精益远达智能科技有限公司 | 一种基于调频连续波毫米波雷达的睡眠监控方法、装置和雷达系统 |
CN110346790A (zh) * | 2019-07-09 | 2019-10-18 | 长沙莫之比智能科技有限公司 | 一种基于毫米波雷达的非接触式生命体征监测方法、装置及系统 |
KR102091974B1 (ko) * | 2019-07-25 | 2020-03-24 | 정기섭 | Fmcw 레이더를 이용한 위치 추적 및 호흡(심박) 패턴 측정 방법 |
CN110464320A (zh) * | 2019-08-22 | 2019-11-19 | 四川长虹电器股份有限公司 | 多目标人体心率和呼吸频率测量系统及方法 |
CN112882016A (zh) * | 2019-11-29 | 2021-06-01 | 南京理工大学 | 基于改进的稳健自适应波束形成的多人生命体征检测方法 |
CN111481184A (zh) * | 2020-04-24 | 2020-08-04 | 华侨大学 | 基于毫米波雷达技术的多目标呼吸心率监测方法与系统 |
CN111990968A (zh) * | 2020-08-07 | 2020-11-27 | 北京航空航天大学杭州创新研究院 | 一种信号监测方法、设备、存储介质及系统 |
CN111856455A (zh) * | 2020-09-04 | 2020-10-30 | 四川长虹电器股份有限公司 | 匹配不同雷达带宽的多目标心率呼吸测量方法及系统 |
CN112686094A (zh) * | 2020-12-03 | 2021-04-20 | 华中师范大学 | 一种基于毫米波雷达的非接触式身份识别方法及系统 |
CN112716462A (zh) * | 2020-12-09 | 2021-04-30 | 北京航空航天大学 | 一种可控照射指向的窄波束毫米波人体心跳/呼吸体征监测装置 |
Non-Patent Citations (1)
Title |
---|
梁福来,等: "UWB MIMO生物雷达多静止人体目标成像方法研究", 雷达学报 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113729674A (zh) * | 2021-09-24 | 2021-12-03 | 青岛海信日立空调系统有限公司 | 一种呼吸检测方法及装置 |
CN113729676A (zh) * | 2021-10-12 | 2021-12-03 | 九州云合(山东)智能科技有限公司 | 一种智能生命体征监护设备 |
CN114403848A (zh) * | 2022-01-05 | 2022-04-29 | 河南牧原智能科技有限公司 | 一种动物群体呼吸频率检测的方法、系统及存储介质 |
CN116035558A (zh) * | 2023-03-02 | 2023-05-02 | 中国科学技术大学 | 基于波束形成的抗干扰呼吸检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113384264B (zh) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113384264B (zh) | 基于雷达的呼吸频率检测方法及睡眠监测设备 | |
Wang et al. | mmHRV: Contactless heart rate variability monitoring using millimeter-wave radio | |
US10548503B2 (en) | Health related diagnostics employing spectroscopy in radio / microwave frequency band | |
Islam et al. | Non-contact vital signs monitoring for multiple subjects using a millimeter wave FMCW automotive radar | |
CN112816960B (zh) | 车内生命探测方法、装置、设备和存储介质 | |
CN110584631A (zh) | 一种基于fmcw雷达的静态人体心跳和呼吸信号提取方法 | |
CN112415485B (zh) | 毫米波雷达的角度超分辨方法、装置及终端设备 | |
CN111142102B (zh) | 一种呼吸数据计算方法以及相关设备 | |
Cho et al. | Detection of heart rate through a wall using UWB impulse radar | |
EP3157420A1 (en) | Monitoring the body using microwaves | |
EP2858566B1 (en) | Multi-spectral imaging system and method for remote biometric measurement of human physiological parameters | |
Yakovleva et al. | EEG analysis in structural focal epilepsy using the methods of nonlinear dynamics (Lyapunov exponents, Lempel–Ziv complexity, and multiscale entropy) | |
CN108177685A (zh) | 一种基于健康监测的新型交通工具 | |
Ding et al. | A novel real-time human heart rate estimation method for noncontact vital sign radar detection | |
US20190320940A1 (en) | Apparatus and method for determining bio-information of target using impulse radar | |
KR20170019785A (ko) | 원형편파 헬릭스 안테나를 이용한 고감도 비침습 생체신호 검출 장치 및 방법 | |
CN113854990A (zh) | 一种心跳检测方法及装置 | |
CN114966656A (zh) | 一种基于毫米波设备的定位方法及装置 | |
Hunter et al. | Distribution of independent components of binocular natural images | |
Chang et al. | Dl-aided nomp: a deep learning-based vital sign estimating scheme using fmcw radar | |
CN116893415B (zh) | 基于毫米波稀疏阵列系统设计的快速成像方法及系统 | |
CN116269297A (zh) | 基于雷达的生命体征检测方法、装置、雷达及存储介质 | |
CN113925481A (zh) | 测量生理信号方法及生理信号测量装置 | |
Walid et al. | Accuracy assessment and improvement of FMCW radar-based vital signs monitoring under Practical Scenarios | |
KR101814290B1 (ko) | 임펄스 레이더를 이용한 타겟의 위치 정보 결정 장치 및 방법 |
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 |