专利名称：生命探测雷达中的呼吸与心跳信号的提取方法及装置的制作方法
生命探测雷达中的呼吸与心跳信号的提取方法及装置技术领域：
本发明涉及信号处理技术领域，尤其涉及一种生命探测雷达中的呼吸与心跳信 号的提取方法及装置。背景技术：
穿墙生命探测雷达融合了雷达和生物医学工程技术，可穿透非金属介质(如砖 墙、废墟等)，非接触式、远距离检测到人体的呼吸和心脏跳动等信息，广泛应用于震后 救援、特殊病人监护以及反恐斗争等领域。
穿墙生命探测雷达的基本原理是发射机通过振荡器产生连续波信号 COS(2 JIfot),由发射天线向目标辐射，连续波信号穿透非金属障碍物，遇到目标后返回。 回波信号经过接收天线进入低噪声放大器，混频解调后再进行预处理，射频干扰和杂波 干扰则被过滤掉，得到基带信号，基带信号经过数模转换后送入信号处理单元，信号处 理单元提取出呼吸和心跳信号。如果目标是人体，则由人的呼吸和心跳共同引起的胸腔 周期性的运动将导致回波信号产生一个相移，相当于对回波进行了调制。其中回波信号 的表达式为
R (t) = cos[2 π f0t-4 π d0/ λ -4 π χ (t) / λ + Φ (t_2d0/c)]
其中，fO为载波频率，λ为载波波长，c为光速，Φ (t)为发射机的相位噪声， d0为雷达-目标距离；χ⑴=Δ lCos (ω !t) + Δ 2cos (ω 2t),代表胸腔运动，A1A2Co1Co2分 别为呼吸和心跳信号的幅度与频率。经过混频解调，得到基带信号为
B (t) = cos[ θ +4 π [ Δ lCos (ω jt) + Δ 2cos (ω 2t) ]/ λ + Δ Φ (t)]
式中θ是距离礼产生的恒定相移，Δ Φ (t)是总的剩余相位噪声。基带信号携 带呼吸和心跳信息。对B(t)进行贝塞尔展开，可知B (t)中包含In1 ω Jm2Co2(mi，m2为 任意整数)组成的谐波分量。经过进一步的滤波和信号处理便可以得到代表呼吸与心跳 信息的ω^ ω2分量。
雷达接收机采集的人体回波复合信号具有非平稳特征，随机性强，复合信号中 心跳信号收到了呼吸运动和体表微动回波信号的压制和干扰，使检测系统无法检测出规 则的心跳信号，因此，必须才去有效的信号处理方法将呼吸和心跳信号分离，提取出比 呼吸信号能力小得多的心跳信号。
目前的呼吸和心跳信号提取方法通常基于数字滤波技术或快速傅立叶变换 (FFT)从频域上入手，因而呼吸信号谐波干扰和目标信号非平稳特性的影像未得到很好 解决。通常需采用自适应处理的方法，将提取的呼吸信号作为参考输入，与同时携带呼 吸、心跳频率成分的信号进行自适应处理，从而实现消除呼吸谐波干扰，检测出心跳波 形的目的。但该方法的前提是必须对呼吸信号基频做精确估计，才能达到理想的结果。 也有运用小波变换理论对检测信号进行分析，以求从混合信号中提取规律性较强的心跳 信号，但如何选择恰当的小波基函数是难点，并且基函数选定后，不能随输入信号的变 化而改变，不具备自适应性。
在实际应用中，呼吸与心跳信号究竟是以何种关系形成是我们选择信号处理关 键所在。线性非时变系统比较容易分析，具有精确的数学解析表达式，同时还可以设计 出各种线性非时变系统，用以实现各种信号处理功能。线性非时变系统满足叠加原理， 在分离加法合成的信号时特别有用。对于占据不同频段的两个信号分量，可用一个线性 滤波器在时域上实现分离。然而，雷达接收信号受系统非线性影像，会经常遇到不属于 相加性组合的信号，若单处的依靠线性滤波来分离和处理这些信号分量，无法达到预期 效果。
发明内容
基于此，有必要提供一种更精确的生命探测雷达中的呼吸与心跳信号的提取方法。
一种生命探测雷达中的呼吸与心跳信号的提取方法，包括以下步骤对雷达接 收信号进行预处理；对预处理后的信号进行经验模态分解，得到有限个固有模态函数； 对分解出的每个固有模态函数进行快速傅立叶变换；分别计算每个固有模态函数对应的 总能量、呼吸频带能量和心跳频带能量；根据所述总能量与呼吸频带能量、心跳频带能 量之间的比值关系选取对应的固有模态函数分别重构呼吸和心跳信号。
其中，所述对雷霆接收信号进行预处理的步骤为对雷达接收信号低通滤波， 滤除高频噪声。
其中，按照如下公式对分解出的每个回有模态函数进行快速傅立叶变换
权利要求
1.一种生命探测雷达中的呼吸与心跳信号的提取方法，包括以下步骤 对雷达接收信号进行预处理；对预处理后的信号进行经验模态分解，得到有限个固有模态函数； 对分解出的每个固有模态函数进行快速傅立叶变换； 分别计算每个固有模态函数对应的总能量、呼吸频带能量和心跳频带能量； 根据所述总能量与呼吸频带能量、心跳频带能量之间的比值关系选取对应的固有模 态函数分别重构呼吸和心跳信号。
2.如权利要求1所述的生命探测雷达中的呼吸与心跳信号的提取方法，其特征在于， 所述对雷达接收信号进行预处理的步骤为对雷达接收信号低通滤波，滤除高频噪声。
3.如权利要求1所述的生命探测雷达中的呼吸与心跳信号的提取方法，其特征在于， 按照如下公式对分解出的每个固有模态函数进行快速傅立叶变换
4.如权利要求3所述的生命探测雷达中的呼吸与心跳信号的提取方法，其特征在于， 在频域上按照如下公式计算每个固有模态函数对应的总能量
5.如权利要求4所述的生命探测雷达中的呼吸与心跳信号的提取方法，其特征在于， 计算所述总能量与呼吸频带能量、心跳频带能量之间的比值关系为
6.—种生命探测雷达中的呼吸与心跳信号的提取装置，其特征在于，包括 预处理模块，用于对雷达接收信号进行预处理；经验模态分解模块，用于对预处理后的信号进行经验模态分解，得到有限个固有模 态函数；傅立叶变换模块，用于对分解出的每个固有模态函数进行快速傅立叶变换； 计算模块，用于分别计算每个固有模态函数对应的总能量、呼吸频带能量和心跳频 带言糧；重构模块，根据所述总能量与呼吸频带能量、心跳频带能量之间的比值关系选取对 应的固有模态函数分别重构呼吸和心跳信号。
7.如权利要求6所述的生命探测雷达中的呼吸与心跳信号的提取装置，其特征在于， 所述预处理模块为低通滤波器，用于对雷达接收信号进行低通滤波，滤除高频噪声。
8.如权利要求6所述的生命探测雷达中的呼吸与心跳信号的提取装置，其特征在于， 所述傅立叶变换模块按照如下公式对分解出的每个有模态函数进行快速傅立叶变换
9.如权利要求8所述的生命探测雷达中的呼吸与心跳信号的提取方法，其特征在于， 所述计算模块包括总能量计算模块，用于按照如下公式在频域上计算每个固有模态函数对应的总能量
10.如权利要求9所述的生命探测雷达中的呼吸与心跳信号的提取装置，其特征在 于，所述重构模块包括比值计算模块，按照如下公式计算所述总能量与呼吸频带能量、心跳频带能量之间 的比值关系
全文摘要
本发明提供了一种生命探测雷达中的呼吸与心跳信号的提取方法及装置，所述方法包括以下步骤对雷达接收信号进行预处理；对预处理后的信号进行经验模态分解，得到有限个固有模态函数；对分解出的每个固有模态函数进行快速傅立叶变换；分别计算每个固有模态函数对应的总能量、呼吸频带能量和心跳频带能量；根据所述总能量、呼吸频带能量和心跳频带能量之间的比值关系选取对应的固有模态函数分别重构呼吸和心跳信号。采用本发明提供的方法及装置，能更精确的提取呼吸与心跳信号，且具有自适应的特点，能够适应环境的变化。
文档编号A61B5/0205GK102018503SQ201010517158
公开日2011年4月20日 申请日期2010年10月21日 优先权日2010年10月21日
发明者吴晶, 李烨, 潘水洋, 王强 申请人:中国科学院深圳先进技术研究院