发明名称 用于生命探测雷达的微弱呼吸信号检测方法
摘要 本发明提供了一种用于生命探测雷达的微弱呼吸信号检测方法。该微弱呼吸信号检测方法包括:由雷达回波信号,建立M×N维的距离慢时信号矩阵S(m,n);对于每个距离单元,在S(m,n)的慢时方向上进行自相关处理,得到信号矩阵R<sub>x</sub>(m,τ);分别对R<sub>x</sub>(m,τ)的各个行向量进行FFT变换,得到距离频率矩阵<img file="DDA0000575162490000011.GIF" wi="246" he="79" />在<img file="DDA0000575162490000012.GIF" wi="221" he="82" />的距离方向上进行加权累积,获得呼吸信号的频谱<img file="DDA0000575162490000013.GIF" wi="168" he="79" />以及遍历<img file="DDA0000575162490000014.GIF" wi="150" he="79" />的每一个位置,寻找能量比大于预设门限的位置,如果该位置的数目大于0时,判断地下存在生命体。本发明实现了微弱呼吸信号的增强、检测与提取,满足了灾害救援等场合下掩埋人员微弱生命特征检测的需要。
申请公布号 CN104224188B 申请公布日期 2016.06.01
申请号 CN201410488161.X 申请日期 2014.09.23
申请人 中国科学院电子学研究所 发明人 吴世有;陈洁;孟升卫;方广有
分类号 A61B5/11(2006.01)I 主分类号 A61B5/11(2006.01)I
代理机构 中科专利商标代理有限责任公司 11021 代理人 曹玲柱
主权项 一种用于生命探测雷达的微弱呼吸信号检测方法,其特征在于,包括:步骤A:由雷达回波信号,建立M×N维的距离慢时信号矩阵S(m,n),其中,M为距离单元数,N为慢时时刻,m=0,1,…,M‑1,n=0,1,…,N‑1;步骤B:对于每个距离单元,在距离慢时信号矩阵S(m,n)的慢时方向上进行自相关处理,得到信号矩阵R<sub>x</sub>(m,τ);步骤C:分别对信号矩阵R<sub>x</sub>(m,τ)的各个行向量进行FFT变换,得到距离频率矩阵<img file="FDA0000932233410000011.GIF" wi="247" he="86" />其中,f为一索引值,f=0,1,…,2N‑2;步骤D:在距离频率矩阵<img file="FDA0000932233410000012.GIF" wi="223" he="87" />的距离方向上进行加权累积,获得呼吸信号频谱<img file="FDA0000932233410000013.GIF" wi="175" he="79" />以及步骤E:遍历呼吸信号频谱<img file="FDA0000932233410000014.GIF" wi="153" he="78" />的每一个位置,寻找能量比大于预设门限的呼吸信号频谱<img file="FDA0000932233410000015.GIF" wi="153" he="78" />位置,如果该位置的数目大于0时,判断地下存在生命体;如果等于0,则判断地下不存在生命体;其中,所述步骤B包括:步骤B1:对于距离慢时信号矩阵S(m,n)中第m个距离单元上的慢时信号x<sub>m</sub>(n),即信号矩阵S(m,n)的第m个行向量,计算其自相关函数<img file="FDA0000932233410000019.GIF" wi="742" he="82" />其中τ=n<sub>1</sub>‑n<sub>2</sub>,n<sub>1</sub>=0,1,…,N‑1,n<sub>2</sub>=0,1,…,N‑1,E[·]为求期望值;以及步骤B2:将不同的距离单元对应的自相关函数<img file="FDA00009322334100000110.GIF" wi="178" he="75" />作为行向量,构成M×(2N‑1)维的信号矩阵R<sub>x</sub>(m,τ),其中,2N‑1是慢时卷积长度,τ=0,1,…,2N‑2。
地址 100190 北京市海淀区北四环西路19号