CN110448290A - 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及*** - Google Patents

一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及*** Download PDF

Info

Publication number
CN110448290A
CN110448290A CN201910764062.2A CN201910764062A CN110448290A CN 110448290 A CN110448290 A CN 110448290A CN 201910764062 A CN201910764062 A CN 201910764062A CN 110448290 A CN110448290 A CN 110448290A
Authority
CN
China
Prior art keywords
signal
echo
heart rate
terahertz
rate detection
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
Application number
CN201910764062.2A
Other languages
English (en)
Other versions
CN110448290B (zh
Inventor
梁晓林
邓建钦
张超群
刘朝阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Electronics Technology Instruments Co Ltd CETI
Original Assignee
China Electronics Technology Instruments Co Ltd CETI
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Electronics Technology Instruments Co Ltd CETI filed Critical China Electronics Technology Instruments Co Ltd CETI
Priority to CN201910764062.2A priority Critical patent/CN110448290B/zh
Publication of CN110448290A publication Critical patent/CN110448290A/zh
Application granted granted Critical
Publication of CN110448290B publication Critical patent/CN110448290B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 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
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02444—Details of sensor
    • 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Physiology (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开了一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***，该方法包括：接收心率检测实时的I/Q两路回波信号；分别抑制所述回波信号中的静态杂波和非静态杂波；对杂波抑制后的回波信号进行复信号解调；采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；对回波信号中心跳信号信噪比进行增强；采用小波变换对信噪比增强后的回波信号进行分解和重构；进行基于变时窗选择的信号频谱估计；抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。

Description

一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置 及***
技术领域
本发明属于心率检测的技术领域，尤其涉及一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***。
背景技术
无论是在医学中还在自然灾害救援中，精确估计人员的呼吸频率、心率均有广泛的应用前景。如在医学中可以用来诊断或预防呼吸类和心血管类疾病，在自然灾害中可以为救援人员提供可靠的依据，降低人员的生命财产的损失。在呼吸频率与心率检测过程中，由于心跳幅度处于毫米级，相比于呼吸幅度更小，通常情况下在远距离且穿透障碍物(如废墟等)的情况下很难准确估计。
目前，常用的心率检测手段包括接触式与非接触式两种。其中，接触式方法主要是指利用传统的心电图、光电容积描记波信号等手段通过在人员身体上粘贴电极片、佩戴传感器等设备实现心率的检测。然而，这会给患者造成一定的困扰及心理负担，极易产生排斥情绪。鉴于上述手段的局限性，随着雷达在医疗、工程等方面的广泛应用，非接触检测手段引发广泛关注。
利用雷达实现心率的提取是一种被广泛应用的非接触信号检测方式，具有全天候、高精度等优势。然而，发明人在研究过程中发现，目前非接触式心率提取手段多是基于微波与低频毫米波波段，鉴于心跳幅度通常小于4mm，在利用连续波雷达检测时，由心跳导致的微多普勒效应在微波与低频毫米波频段不敏感，很难实现精准检测。尤其是在远距离穿墙情况下，微波与低频毫米波雷达更难实现心率的提取。
发明内容
为克服上述现有技术的不足，本发明提供了一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***，利用太赫兹雷达估计心率，实现远距离穿墙情况下心率的精确检测。
根据本公开的一个或多个实施例的一个方面，提供一种用于远距离人员心率检测的太赫兹穿墙雷达。
一种用于远距离人员心率检测的太赫兹穿墙雷达，包括：射频输出端和射频接收端，将倍频、射频上变频、再耦合后的太赫兹波作为本振信号，同时通过第一极化器自所述射频输出端输出，所述射频接收端通过第二极化器后实时接收I/Q两路回波信号。
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
一种基于太赫兹穿墙雷达的远距离人员心率检测方法，该方法包括：
接收心率检测实时的I/Q两路回波信号；
分别抑制所述回波信号中的静态杂波和非静态杂波；
对杂波抑制后的回波信号进行复信号解调；
采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
对回波信号中心跳信号信噪比进行增强；
采用小波变换对信噪比增强后的回波信号进行分解和重构；
进行基于变时窗选择的信号频谱估计；
抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
进一步地，该方法中，所述静态杂波根据背景滤波算法进行估计，并采用实时接收的I/Q两路回波信号分别减去估计的两路静态杂波进行静态杂波的抑制。
进一步地，该方法中，对杂波抑制后的回波信号根据雷达与生命体之间的距离、心跳幅度、心率、波长和相位噪声进行复信号解调，并对复信号调节后的信号进行对数处理。
进一步地，该方法中，采用5阶巴特沃斯滤波器对回波信号中的高频和低频噪声干扰进行抑制。
进一步地，该方法中，采用滑动滤波器对回波信号中心跳信号信噪比进行增强。
进一步地，该方法中，采用小波变换对信噪比增强后的回波信号进行分解和重构中采用的小波函数为Morlet小波。
进一步地，该方法中，所述进行基于变时窗选择的信号频谱估计的具体步骤包括：
对于具有固定时长的信号，频率分辨率为：i＝1,2,…,q，其中，fs为频率，wi表示第i个窗口的长度，q表示窗口数目；
对每一个时窗内的信号进行傅里叶变换，同时对频谱进行加窗处理，窗口宽度为0.8-2.4Hz；
对于得到的q组频谱进行累计运算。
根据本公开的一个或多个实施例的一个方面，提供一种计算机可读存储介质。
一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
根据本公开的一个或多个实施例的一个方面，提供一种终端设备。
一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测装置。
一种基于太赫兹穿墙雷达的远距离人员心率检测装置，基于所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，包括：
信号接收模块，被配置为接收心率检测实时的I/Q两路回波信号；
杂波抑制模块，被配置为分别抑制所述回波信号中的静态杂波和非静态杂波；
复信号解调模块，被配置为对杂波抑制后的回波信号进行复信号解调；
噪声干扰抑制模块，被配置为采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
信噪比增强模块，被配置为对回波信号中心跳信号信噪比进行增强；
分解重构模块，被配置为采用小波变换对信噪比增强后的回波信号进行分解和重构；
频谱估计模块，被配置为进行基于变时窗选择的信号频谱估计；
谐波干扰抑制模块，被配置为抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测***。
一种基于太赫兹穿墙雷达的远距离人员心率检测***，包括依次连接的所述太赫兹穿墙雷达、多路采集设备和所述的一种基于太赫兹穿墙雷达的远距离人员心率检测装置。
以上一个或多个技术方案存在以下有益效果：
本发明公开的一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***，分别利用背景滤波算法与LTS算法抑制回波信号中的杂波干扰；对I/Q两路信号进行复信号解调，及进行对数处理；利用带通滤波器抑制高、低频噪声；利用滑动滤波器改善信号信噪比；利用小波变换对心跳信号进行重构；设计变时窗策略，循环利用有限的时域信号；对不同长度的时域信号频谱进行累积；通过频域加窗估计心率，有效实现远距离穿墙情况下心率的精确检测。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明一个或多个实施例的一种超大规模数据库缺失数据填充方法及装置总体流程图；
图2为本发明一个或多个实施例的射频输出端示意图；
图3为本发明一个或多个实施例的射频接收端示意图；
图4为本发明一个或多个实施例的雷达采集的I/Q信号示意图；
图5为本发明一个或多个实施例的滤除杂波干扰后的I/Q信号示意图；
图6为本发明一个或多个实施例的回波信号处理后的信号示意图；
图7为本发明一个或多个实施例的经过滤波后的回波信号示意图；
图8为本发明一个或多个实施例的滑动滤波后的信号示意图；
图9为本发明一个或多个实施例的重构后的心跳信号示意图；
图10为本发明一个或多个实施例的变时窗选择示意图；
图11为本发明一个或多个实施例的基于变时窗技术的信号频谱示意图；
图12为本发明一个或多个实施例的基最终的心率信号示意图。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
需要注意的是，附图中的流程图和框图示出了根据本公开的各种实施例的方法和***的可能实现的体系架构、功能和操作。应当注意，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，所述模块、程序段、或代码的一部分可以包括一个或多个用于实现各个实施例中所规定的逻辑功能的可执行指令。也应当注意，在有些作为备选的实现中，方框中所标注的功能也可以按照不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，或者它们有时也可以按照相反的顺序执行，这取决于所涉及的功能。同样应当注意的是，流程图和/或框图中的每个方框、以及流程图和/或框图中的方框的组合，可以使用执行规定的功能或操作的专用的基于硬件的***来实现，或者可以使用专用硬件与计算机指令的组合来实现。
在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
实施例一
根据本公开的一个或多个实施例的一个方面，提供一种用于远距离人员心率检测的太赫兹穿墙雷达。
如图2和图3所示，一种用于远距离人员心率检测的太赫兹穿墙雷达，包括：射频输出端和射频接收端，将倍频、射频上变频、再耦合后的太赫兹波作为本振信号，同时通过第一极化器自所述射频输出端输出，所述射频接收端通过第二极化器后实时接收I/Q两路回波信号。
本实施例中的THz雷达接收机(零中频)能够实时接收I/Q两路回波信号。本实施例的THz雷达中，为了有效提高参考本振信号与发射信号的相位一致性，提供***识别精度，将倍频、射频上变频、再耦合后的THz波作为本振信号，此外通过线性的频谱搬移降低信号的非线性相位误差，降低***数据补偿难度，如图2所示。在输出端与接收端增加极化器，能够有效消除环境杂波，提高接收机信噪比，获取目标特性，增加目标识别维度，接收端如图3所示。
为了有效提高参考本振信号与发射信号的相位一致性，提供***识别精度，将倍频、射频上变频、再耦合后的THz波作为本振信号，此外通过线性的频谱搬移降低信号的非线性相位误差，降低***数据补偿难度。在输出端与接收端增加极化器，能够有效消除环境杂波，提高接收机信噪比，获取目标特性，增加目标识别维度。
实施例二
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
一种基于太赫兹穿墙雷达的远距离人员心率检测方法，该方法包括：
步骤1：接收心率检测实时的I/Q两路回波信号；
步骤2：分别抑制所述回波信号中的静态杂波和非静态杂波；
步骤3：对杂波抑制后的回波信号进行复信号解调；
步骤4：采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
步骤5：对回波信号中心跳信号信噪比进行增强；
步骤6：采用小波变换对信噪比增强后的回波信号进行分解和重构；
步骤7：进行基于变时窗选择的信号频谱估计；
步骤8：抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
在本实施例步骤1中采集THz雷达采集的I/Q信号如图4所示。
在本实施例步骤2杂波抑制中，步骤2主要是抑制回波信号中杂波干扰，主要包括静态杂波与非静态杂波。其中，静态杂波主要来自于测试环境中的静态目标，可以认为其是时不变信号。本发明中利用背景滤波算法对其进行估计，可表示为
其中，n＝1,…,N，表示数字化的回波信号的样本数据。为了抑制静态杂波，将数字化的回波I/Q的每一个样本分别减去式(1)与式(2)的结果，可表示为
抑制回波中的非静态杂波，可得到：
W1＝IΤ-X(XΤX)-1XΤIT (5)
W2＝QΤ-X(XΤX)-1XΤQT (6)
其中，X＝[x1,x2],x1＝[0,1,…,N-1]Τ,滤除杂波干扰后的I/Q信号如图5所示。
在本实施例步骤3信号解调中，
对I/Q信号进行复信号解调(complex signal demodulation,CSD)，可表示为：
其中，d0表示雷达与生命体之间的距离，mh心跳幅度，fh表示心率，λ表示波长，表示相位噪声。
为了抑制解调后信号中存在的常数项干扰，对信号进行对数处理，可表示为：
由此可知，心跳信号主要存在于上式的虚部，因此，对式(8)进行虚部处理，可得到：
回波信号处理后的信号如图6所示。
在本实施例步骤4高低频杂波抑制中，
经过步骤3的处理后，接下来主要是利用5阶巴特沃斯滤波器抑制回波中存在的高频与低频噪声干扰，可表示为：
其中，κi andχi表示滤波系数。经过滤波后的回波信号如图7所示。
在本实施例步骤5信噪比增强中，
本发明利用滑动滤波器实现心跳信号信噪比的增强，可表示为：
其中，表示小于N/7的最大整数。滑动滤波后的信号图8所示。
在本实施例步骤6信号重构中，
为了抑制带外噪声对心跳信号的影响，利用小波变换(Wavelet Transform，WT)实现信号的分解与重构。在对信号进行分解过程中，小波变换可表示为：
其中，ψ((n-b)/a)表示变换过程中采用的小波函数，a表示尺度因子，b表示扩展因子。
本发明中采用的小波函数为Morlet小波，可表示为：
利用分解后的扩展因子与尺度因子，信号可重构为
其中，
重构后的心跳信号如图9所示。
在本实施例步骤7基于变时窗选择的信号频谱估计中，
为了有效估计心率，步骤7中提出了一种变时长的信号选择策略。通常，对于具有固定时长的信号，频率分辨率为：
Δf＝1/Tw (17)
为了准确估计心率，必须满足：
Δf＜＜fh (18)
即
fr＝ρ×Δf (19)
其中，ρ表示正整数。
为了满足上式要求，同时提高估计精度，本发明提出如图10所示的变时窗选择策略，此种情况下频率分辨率可表示为：
其中，ωi表示第i个窗口的长度，q表示窗口数目。
为了估计心率，对每一个时窗内的信号进行傅里叶变换，同时对频谱进行加窗处理，窗口宽度为0.8-2.4Hz。此外，对于得到的q组频谱进行累计运算，可表示为：
基于变时窗技术的信号频谱如图11所示。
在本实施例步骤8带内谐波抑制中，
由步骤7得到的信号频谱除了心率，还存在谐波干扰。因此，需要对谐波进行抑制。本发明中对带内谐波的抑制可表示为：
H[i]＝L[i]+jL[i] (22)
其中，
κ表示频谱分量。
由此心率可估计为：
fr＝H[μr] (24)
其中，μr表示累积频谱的峰值。最终的心率信号如图12所示。
实施例三
根据本公开的一个或多个实施例的一个方面，提供一种计算机可读存储介质。
一种计算机可读存储介质，其中存储有多条指令，所述指令适于由终端设备的处理器加载并执行所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
实施例四
根据本公开的一个或多个实施例的一个方面，提供一种终端设备。
一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，所述指令适于由处理器加载并执行所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
这些计算机可执行指令在设备中运行时使得该设备执行根据本公开中的各个实施例所描述的方法或过程。
在本实施例中，计算机程序产品可以包括计算机可读存储介质，其上载有用于执行本公开的各个方面的计算机可读程序指令。计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
本公开所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开内容操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开内容的各个方面。
实施例五
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测装置。
一种基于太赫兹穿墙雷达的远距离人员心率检测装置，基于所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，包括：
信号接收模块，被配置为接收心率检测实时的I/Q两路回波信号；
杂波抑制模块，被配置为分别抑制所述回波信号中的静态杂波和非静态杂波；
复信号解调模块，被配置为对杂波抑制后的回波信号进行复信号解调；
噪声干扰抑制模块，被配置为采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
信噪比增强模块，被配置为对回波信号中心跳信号信噪比进行增强；
分解重构模块，被配置为采用小波变换对信噪比增强后的回波信号进行分解和重构；
频谱估计模块，被配置为进行基于变时窗选择的信号频谱估计；
谐波干扰抑制模块，被配置为抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
以上实施例三、四和五的装置中涉及的各步骤与方法实施例二相对应，具体实施方式可参见实施例一的相关说明部分。术语“计算机可读存储介质”应该理解为包括一个或多个指令集的单个介质或多个介质；还应当被理解为包括任何介质，所述任何介质能够存储、编码或承载用于由处理器执行的指令集并使处理器执行本发明中的任一方法。
实施例六
根据本公开的一个或多个实施例的一个方面，提供一种基于太赫兹穿墙雷达的远距离人员心率检测***。
一种基于太赫兹穿墙雷达的远距离人员心率检测***，包括依次连接的所述太赫兹穿墙雷达、多路采集设备和所述的一种基于太赫兹穿墙雷达的远距离人员心率检测装置。
在本实施例中利用多路采集设备实现上述信号的精确采集，数字化的回波信号分别记为I与Q，储存在电脑中，以便后期数据处理应用。
以上一个或多个技术方案存在以下有益效果：
本发明公开的一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***，分别利用背景滤波算法与LTS算法抑制回波信号中的杂波干扰；对I/Q两路信号进行复信号解调，及进行对数处理；利用带通滤波器抑制高、低频噪声；利用滑动滤波器改善信号信噪比；利用小波变换对心跳信号进行重构；设计变时窗策略，循环利用有限的时域信号；对不同长度的时域信号频谱进行累积；通过频域加窗估计心率，有效实现远距离穿墙情况下心率的精确检测。
本领域技术人员应该明白，上述本发明的各模块或各步骤可以用通用的计算机装置来实现，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。本发明不限制于任何特定的硬件和软件的结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (10)

1.一种用于远距离人员心率检测的太赫兹穿墙雷达，其特征在于，包括：射频输出端和射频接收端，将倍频、射频上变频、再耦合后的太赫兹波作为本振信号，同时通过第一极化器自所述射频输出端输出，所述射频接收端通过第二极化器后实时接收I/Q两路回波信号。
2.一种基于太赫兹穿墙雷达的远距离人员心率检测方法，其特征在于，该方法包括：
接收心率检测实时的I/Q两路回波信号；
分别抑制所述回波信号中的静态杂波和非静态杂波；
对杂波抑制后的回波信号进行复信号解调；
采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
对回波信号中心跳信号信噪比进行增强；
采用小波变换对信噪比增强后的回波信号进行分解和重构；
进行基于变时窗选择的信号频谱估计；
抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
3.如权利要求2所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，其特征在于，所述静态杂波根据背景滤波算法进行估计，并采用实时接收的I/Q两路回波信号分别减去估计的两路静态杂波进行静态杂波的抑制。
4.如权利要求2所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，其特征在于，对杂波抑制后的回波信号根据雷达与生命体之间的距离、心跳幅度、心率、波长和相位噪声进行复信号解调，并对复信号调节后的信号进行对数处理。
5.如权利要求2所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，其特征在于，该方法中，采用5阶巴特沃斯滤波器对回波信号中的高频和低频噪声干扰进行抑制；
采用滑动滤波器对回波信号中心跳信号信噪比进行增强。
采用小波变换对信噪比增强后的回波信号进行分解和重构中采用的小波函数为Morlet小波。
6.如权利要求2所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，其特征在于，该方法中，所述进行基于变时窗选择的信号频谱估计的具体步骤包括：
对于具有固定时长的信号，频率分辨率为：其中，fs为频率，wi表示第i个窗口的长度，q表示窗口数目；
对每一个时窗内的信号进行傅里叶变换，同时对频谱进行加窗处理，窗口宽度为0.8–2.4Hz；
对于得到的q组频谱进行累计运算。
7.一种计算机可读存储介质，其中存储有多条指令，其特征在于，所述指令适于由终端设备的处理器加载并执行如权利要求2-6任一项所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
8.一种终端设备，其包括处理器和计算机可读存储介质，处理器用于实现各指令；计算机可读存储介质用于存储多条指令，其特征在于，所述指令适于由处理器加载并执行如权利要求2-6任一项所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法。
9.一种基于太赫兹穿墙雷达的远距离人员心率检测装置，其特征在于，基于如权利要求2-6任一项所述的一种基于太赫兹穿墙雷达的远距离人员心率检测方法，包括：
信号接收模块，被配置为接收心率检测实时的I/Q两路回波信号；
杂波抑制模块，被配置为分别抑制所述回波信号中的静态杂波和非静态杂波；
复信号解调模块，被配置为对杂波抑制后的回波信号进行复信号解调；
噪声干扰抑制模块，被配置为采用滤波器对回波信号中的高频和低频噪声干扰进行抑制；
信噪比增强模块，被配置为对回波信号中心跳信号信噪比进行增强；
分解重构模块，被配置为采用小波变换对信噪比增强后的回波信号进行分解和重构；
频谱估计模块，被配置为进行基于变时窗选择的信号频谱估计；
谐波干扰抑制模块，被配置为抑制信号频谱中除心率的谐波干扰，得到最终的心率信号。
10.一种基于太赫兹穿墙雷达的远距离人员心率检测***，其特征在于，包括依次连接的如权利要求1所述的太赫兹穿墙雷达、多路采集设备和如权利要求2-6任一项所述的一种基于太赫兹穿墙雷达的远距离人员心率检测装置。
CN201910764062.2A 2019-08-19 2019-08-19 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及*** Active CN110448290B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910764062.2A CN110448290B (zh) 2019-08-19 2019-08-19 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910764062.2A CN110448290B (zh) 2019-08-19 2019-08-19 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***

Publications (2)

Publication Number Publication Date
CN110448290A true CN110448290A (zh) 2019-11-15
CN110448290B CN110448290B (zh) 2022-05-03

Family

ID=68487571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910764062.2A Active CN110448290B (zh) 2019-08-19 2019-08-19 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***

Country Status (1)

Country Link
CN (1) CN110448290B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142084A (zh) * 2019-12-11 2020-05-12 中国电子科技集团公司第四十一研究所 一种微动太赫兹谱识别检测算法
CN111982856A (zh) * 2020-08-27 2020-11-24 中电科仪器仪表有限公司 一种基于太赫兹波的物质无标志检测识别方法
CN114983373A (zh) * 2022-06-02 2022-09-02 谢俊 一种检测人体心率的方法
CN116643247A (zh) * 2023-07-25 2023-08-25 南京航空航天大学 一种舱内活体检测杂波抑制方法和装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5443073A (en) * 1991-09-12 1995-08-22 Drexel University System and method of impedance cardiography monitoring
CN102046076A (zh) * 2008-04-03 2011-05-04 Kai医药公司 非接触式生理运动传感器及其使用方法
CN103674001A (zh) * 2013-11-19 2014-03-26 南京航空航天大学 一种基于增强自适应时频峰值滤波的光纤陀螺去噪方法
CN104434059A (zh) * 2014-10-29 2015-03-25 上海无线电设备研究所 一种利用太赫兹波进行非接触式探测生命体征的探测方法
CN106175723A (zh) * 2016-06-27 2016-12-07 中国人民解放军第三军医大学第附属医院 一种基于fmcw宽带雷达的多生命监护***
KR20170055307A (ko) * 2015-11-11 2017-05-19 (주)비스토스 초음파 도플러 태아감시 장치의 심박 검출 신호처리 방법
CN106805940A (zh) * 2015-12-02 2017-06-09 由国峰 一种连续波生物雷达体征检测装置
CN109375204A (zh) * 2018-10-26 2019-02-22 中电科仪器仪表有限公司 基于雷达的目标检测方法、***、设备及介质
CN109805931A (zh) * 2019-02-18 2019-05-28 中电科仪器仪表有限公司 基于太赫兹多普勒雷达的远距离生命微动信号检测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5443073A (en) * 1991-09-12 1995-08-22 Drexel University System and method of impedance cardiography monitoring
CN102046076A (zh) * 2008-04-03 2011-05-04 Kai医药公司 非接触式生理运动传感器及其使用方法
CN103674001A (zh) * 2013-11-19 2014-03-26 南京航空航天大学 一种基于增强自适应时频峰值滤波的光纤陀螺去噪方法
CN104434059A (zh) * 2014-10-29 2015-03-25 上海无线电设备研究所 一种利用太赫兹波进行非接触式探测生命体征的探测方法
KR20170055307A (ko) * 2015-11-11 2017-05-19 (주)비스토스 초음파 도플러 태아감시 장치의 심박 검출 신호처리 방법
CN106805940A (zh) * 2015-12-02 2017-06-09 由国峰 一种连续波生物雷达体征检测装置
CN106175723A (zh) * 2016-06-27 2016-12-07 中国人民解放军第三军医大学第附属医院 一种基于fmcw宽带雷达的多生命监护***
CN109375204A (zh) * 2018-10-26 2019-02-22 中电科仪器仪表有限公司 基于雷达的目标检测方法、***、设备及介质
CN109805931A (zh) * 2019-02-18 2019-05-28 中电科仪器仪表有限公司 基于太赫兹多普勒雷达的远距离生命微动信号检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋鹏涛: "基于时变窗长时频峰值滤波的地震信号去噪研究", 《中国优秀硕士学位论文全文数据库（电子期刊） 基础科学辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111142084A (zh) * 2019-12-11 2020-05-12 中国电子科技集团公司第四十一研究所 一种微动太赫兹谱识别检测算法
CN111142084B (zh) * 2019-12-11 2023-04-07 中国电子科技集团公司第四十一研究所 一种微动太赫兹谱识别检测算法
CN111982856A (zh) * 2020-08-27 2020-11-24 中电科仪器仪表有限公司 一种基于太赫兹波的物质无标志检测识别方法
CN111982856B (zh) * 2020-08-27 2023-08-01 中电科思仪科技股份有限公司 一种基于太赫兹波的物质无标志检测识别方法
CN114983373A (zh) * 2022-06-02 2022-09-02 谢俊 一种检测人体心率的方法
CN116643247A (zh) * 2023-07-25 2023-08-25 南京航空航天大学 一种舱内活体检测杂波抑制方法和装置
CN116643247B (zh) * 2023-07-25 2023-10-27 南京航空航天大学 一种舱内活体检测杂波抑制方法和装置

Also Published As

Publication number Publication date
CN110448290B (zh) 2022-05-03

Similar Documents

Publication Publication Date Title
CN110448290A (zh) 一种基于太赫兹穿墙雷达的远距离人员心率检测方法、装置及***
Varanini et al. An efficient unsupervised fetal QRS complex detection from abdominal maternal ECG
Mert ECG feature extraction based on the bandwidth properties of variational mode decomposition
CN108056770A (zh) 一种基于人工智能的心率检测方法
US11253204B2 (en) Method for assessing electrocardiogram signal quality
Fonseca et al. Lomb-scargle periodogram applied to heart rate variability study
CN109805931A (zh) 基于太赫兹多普勒雷达的远距离生命微动信号检测方法
CN108338784A (zh) 基于eemd的小波熵阈值的心电信号去噪方法
CN106175731A (zh) 非接触式生命体征监测的信号处理***
CN106889981B (zh) 一种用于提取胎儿心率的智能终端
Satija et al. Noise‐aware dictionary‐learning‐based sparse representation framework for detection and removal of single and combined noises from ECG signal
US20190142289A1 (en) Robust real-time heart rate monitoring method based on heartbeat harmonics using small-scale radar
CN110327029A (zh) 一种基于微波感知的心率与心率变异性监测方法
Mir et al. ECG denoising and feature extraction techniques–a review
Biswal ECG signal analysis using modified S‐transform
Saini et al. Effective automated method for detection and suppression of muscle artefacts from single‐channel EEG signal
Pangerc et al. Robust detection of heart beats in multimodal records using slope-and peak-sensitive band-pass filters
Rodrigues Fetal beat detection in abdominal ECG recordings: global and time adaptive approaches
US11540763B2 (en) Control method and system for filtering power line interferences
CN107890347B (zh) 滤除工频干扰信号的控制方法及***
Figuera et al. Spectrally adapted Mercer kernels for support vector nonuniform interpolation
CN106073800A (zh) 基于绝对差值和提取的动态光谱数据处理方法及其装置
Yamamoto et al. MUSIC-based non-contact heart rate estimation with adaptive window size setting
Yamamoto et al. Non-contact heartbeat detection by MUSIC with discrete cosine transform-based parameter adjustment
Li et al. Robust suppression of nonstationary power-line interference in electrocardiogram signals

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 266555 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Applicant after: CLP kesiyi Technology Co.,Ltd.

Address before: 266555 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province

Applicant before: CHINA ELECTRONICS TECHNOLOGY INSTRUMENTS Co.,Ltd.

GR01 Patent grant
GR01 Patent grant