CN107167802A - 一种基于超宽带雷达的呼吸信号检测算法 - Google Patents
一种基于超宽带雷达的呼吸信号检测算法 Download PDFInfo
- Publication number
- CN107167802A CN107167802A CN201710372374.XA CN201710372374A CN107167802A CN 107167802 A CN107167802 A CN 107167802A CN 201710372374 A CN201710372374 A CN 201710372374A CN 107167802 A CN107167802 A CN 107167802A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- breath signal
- frequency
- noise
- 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.)
- Pending
Links
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
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
-
- 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
- G01S13/08—Systems for measuring distance only
- G01S13/10—Systems for measuring distance only using transmission of interrupted, pulse modulated waves
- G01S13/106—Systems for measuring distance only using transmission of interrupted, pulse modulated waves using transmission of pulses having some particular characteristics
-
- 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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
-
- 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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于超宽带雷达的呼吸信号检测算法。该方法包括Motion Filter滤波、回波数据分块、呼吸信号增强、呼吸信号相关叠加和多重系数衰减的高斯噪声擦除等多个步骤。当雷达观测环境中有人体对象时,对回波数据沿慢时间方向分块,各块数据都包含有呼吸信号,但是各块间噪声不相关。频域中各块数据相关叠加,呼吸信号得到增强,同时噪声稳态随机分布于整个频率观测空间上。单块距离—频率数据呼吸信号增强有效保护弱呼吸信号,叠加和距离—频率数据呼吸信号增强有效增强呼吸信号同时抑制整个频率观测空间中的噪声。多重系数衰减的高斯噪声擦除过程逐步衰减噪声并期望噪声均值衰减为0。该方法能够有效去除白噪声和色噪声成分,检测弱呼吸信号。
Description
技术领域
本发明涉及电磁场与微波技术领域,具体涉及一种基于超宽带雷达的呼吸信号检测算法。
背景技术
人类面临诸多自然灾害和人为灾害,灾害过程中或灾后被困生命受到严重的安全威胁,及时有效的生命定位和搜索成为一个至关重要的问题。老年人的实时健康状态监控是需要关注的问题。在某些犯罪环境和反恐、巷战等军事斗争环境中,犯罪分子或者敌对分子隐藏在各种掩体中,此时需要对他们进行及时准确的判断和定位。对生命体的及时准确的存在性判断和准确的位置判断,对生命呼吸状态的实时监控,以及其他潜在的应用,都会产生巨大价值。因此,非接触式的生命信号检测技术受到越来越多的关注和研究。
当前已经有多种应用超宽带雷达脉冲信号进行呼吸信号检测的算法,包括应用快速傅里叶变换(FFT)的方法、应用小波变换的方法、应用奇异值分解的方法、应用等幅追踪的方法等。但应用FFT方法的检测算法存在高次谐波,应用小波变换方法的检测算法无法获得位置信息,应用奇异值分解方法的算法计算复杂度高,应用等幅追踪方法的检测算法要求源信号信噪比高。本发明提出一种新的实用的基于超宽带雷达的呼吸信号检测算法。本算法针对呼吸信号微弱的问题,提出频域呼吸信号增强的方法,针对噪声种类多且干扰强的问题,提出分块策略和多重系数衰减的高斯噪声擦除方法,并且算法计算复杂度低,检测正确率高,鲁棒性好,具有工程应用价值。
发明内容
本发明提出一种基于超宽带雷达的呼吸信号检测算法,有效抑制多种噪声干扰,有效检测弱呼吸信号并获得一维位置信息,方法步骤如下:
步骤1:超宽带雷达发射冲击脉冲并采样量化回波数据,形成二维慢时间—快时间矩阵r(m,n);
步骤2:对二维矩阵沿慢时间方向应用Motion Filter滤波去除背景杂波获得结果rmf(m,n);
步骤3:对步骤2的结果沿慢时间方向分块为rmf i(m,n),优先分成20~30块,块数设为K,i=1,2,...K;
步骤4:获取各分块频域幅度谱并进行呼吸信号增强,结果为Rmf i(f,l):
对各块二维慢时间—快时间数据沿慢时间方向进行快速傅里叶变换(FFT),并获取二维距离—频率频域幅度谱,对各块频域幅度谱进行呼吸信号增强,步骤如下:
1)设置频率搜索范围F,优先设置F=0.15~0.45Hz,该范围F必须包含呼吸频率f0;
2)设置阈值y和增强比k;
3)扫描每个距离上F范围内的频率幅度值的变化,计算峰值区间内的峰值和PKs并获得峰值区间旁的落差ΔR;
4)找到最大的落差,依据该落差与阈值y的差别,判断该峰是否是呼吸频率,如果是,则进行k倍增强;如果不是,则不对该频点的频率幅度值进行处理;
步骤5:对Rmf i(f,l)进行呼吸信号相关叠加,获得叠加和Rmfs(f,l):
步骤6:对Rmfs(f,l)通过步骤4中的呼吸增强步骤进行呼吸信号增强,获得结果Rs(f,l);
步骤7:对Rs(f,l)去除噪声均值并进行多重系数衰减的高斯噪声擦除过程,步骤如下:
1)获取步骤6后的距离—频域叠加和结果Rs(f,l);
2)由式获得噪声均值Rs(f,l)减去该均值并取绝对值abs为取绝对值函数;
3)获取当前擦除系数a和每个距离上的高斯噪声均值,将每个距离上的频率幅度减去擦除系数与高斯噪声均值的乘积,取绝对值为频率幅度值;
4)迭代20~30次的步骤3),最终获得检测结果。
在步骤1中,应用超宽带雷达采集回波数据时,慢时间采样数据应足量。
在步骤2中,Motion Filter有多种滤波器参数,其中有FIR2、FIR3、FIR4等,可选择其中一种滤波器进行滤波。滤波器如下:
FIR2是一个二阶差分滤波器,使用以下公式:
y(n)=w(1)·x(n)+w(2)·x(n-1)=x(n)-x(n-1) w=[1,-1]
FIR3是一个三阶差分滤波器,使用以下公式:
y(n)=w(1)·x(n)+w(2)·x(n-1)+w(3)·x(n-2) w=[1,-0.8,-0.2]
FIR4是一个四阶差分滤波器,使用以下公式:
y(n)=w(1)·x(n)+w(2)·x(n-1)+w(3)·x(n-2)+w(4)·x(n-3)
w=[1,-0.6,-0.3,-0.1]
在步骤3中,对二维回波数据沿慢时间方向分成多块,相邻两块的距离定义为块间距,块间距优先设置为40~70,块长L根据频率分辨率要求设置,PRF是脉冲重复频率,Δf是频率分辨率。当存在呼吸信号时,每一块数据都会包含呼吸频率信号,而各块间噪声不相关。
在步骤4和步骤6中,应用呼吸信号增强方法进行单块距离—频率数据呼吸信号增强和叠加和距离—频率数据呼吸信号增强。优先设置频率搜索范围F=0.15~0.45Hz,在此范围内搜索呼吸信号;峰值搜索方法中,根据频率分辨率设置峰值区间长度并获得峰值区间内的频率采样点的和作为当前峰值和PKs;落差ΔR是峰值和与紧邻峰值区间的两个频点频率幅度值的和DPKs求差,即,ΔR=PKs-DPKs。
在步骤5中,对多块距离—频率数据进行呼吸信号相关叠加,某个距离上的慢时间回波信号r(t)的频域信号R(f)为:
R(f)=a·δ(f-f0)+N(f)
其中,f0为呼吸频率,δ(f)为冲击信号,N(f)为频域噪声,K块数据Rk(f),k=1,2,……,K
Rk(f)=ak·δ(f-f0)+Nk(f)
则K个距离-频域的Rk(f)叠加和Rs(f)为
其中当某个距离上没有呼吸信号时ak=0,则
当这个距离上有呼吸信号时ak≠0,则
Nk(f)是频域随机噪声,由于多块的叠加效应,Ns(f)是稳态随机分布,则有
n(f)是高斯噪声,呼吸信号相关叠加有效增强呼吸信号,而在单块呼吸信号增强过程中误放大的随机噪声不会产生叠加干扰,整个频率观测空间上的随机噪声的叠加使得其呈现稳态随机分布。
在步骤7中,擦除系数a随擦除次数增多而减小。多重系数衰减的高斯噪声擦除过程逐步减弱整个频率观测空间中的噪声,期望噪声均值衰减为0。设Rs(f,l)为步骤7的2)后的结果,设Rs(f,l,j)为第j次擦除后的距离—频率结果,j=0时为原结果,即
Rs(f,l,0)=Rs(f,l)
在第j(j>0)次擦除时,擦除系数为aj,沿频率向的高斯噪声均值为
则有
如果某个距离上存在呼吸信号,则该距离上只会保留下呼吸信号。擦除系数a随擦除次数x的增多递减,递减函数模型可选择指数模型、对数模型、余弦函数模型,优先选择余弦函数模型a=k·cos(c·j)+b,其中j是擦除次数,j∈[0,J],J优先设置为20~30,a是擦除系数,递减范围优先设置为1.6~0,k、c、b是函数参数。
本发明的优点:
本发明对回波数据分块,利用各块数据间呼吸信号的频率相关性和噪声的频率不相关性,进行呼吸信号相关叠加,呼吸信号得到有效增强而随机分布的噪声信号不会得到增强;本发明应用呼吸信号增强的方法对单块呼吸信号增强,显著保护弱呼吸信号;对叠加和呼吸信号增强,则有效抑制整个频率观测空间上的多种噪声;多重系数衰减的高斯噪声擦除过程逐步减弱整个频率观测空间上的噪声均值,期望噪声均值衰减为0。该算法各个步骤都有效提高信噪比,并且计算效率高。
附图说明
图1本发明方法流程图;
图2是回波数据分块示意图;
图3是呼吸增强示意图;
图4是呼吸信号相关叠加示意图;
图5是多重系数衰减的高斯噪声擦除过程示意图;
图6为对实施样例的二维回波数据应用本发明进行检测的结果。
具体实施方式
通过实施样例,进一步阐述本发明。本发明方法流程图如图1所示。
(1)本样例中,观测环境中设置人体目标,使用超宽带雷达PulseON440以8Hz的脉冲重复频率对观测环境发射脉冲并采样量化回波数据,形成二维慢时间—快时间数据r(m,n)。
(2)对回波数据r(m,n)在慢时间方向应用四阶差分Motion Filter FIR4滤波去除背景噪声,获得结果rmf(m,n)。
(3)对rmf(m,n)沿慢时间方向分块为rmf i(m,n)。分块示意图如图2所示。图2中,虚线阵表示回波慢时间—快时间矩阵。三个实线矩形方框1、2、3各表示一个分块。块间距、块长的定义如图所示。块长以频率分辨率为依据,分辨率要求越高则块长越长;块间距最小为1,最大为块长。本样例中块长设置为256,块间距设置为64,块数为30块。
(4)获取各分块慢时间方向频域幅度谱,对各块频域幅度谱进行呼吸信号增强,结果为Rmf i(f,l)。呼吸增强方法示意图如图3所示,当搜索到呼吸频率处有峰值时,对该峰值进行扩大。由于单块呼吸信号弱,阈值应设置为较小值,增强比也应设置较小。呼吸增强过程如下:
1)设置频率检测范围F=0.15~0.45Hz;
2)设置阈值为y=0.2,设置增强比为k=2;
3)扫描每个距离上F范围内的频率幅度值的变化,计算峰值区间内的峰值和PKs并获得峰值区间旁的落差ΔR:
设置峰值区间长度为3,获得峰值区间内的频率采样点的和作为当前峰值和,落差是峰值和与紧邻峰值区间的两个频点频率幅度值的和DPKs求差,即,ΔR=PKs-DPKs。
4)找到最大的落差,依据该落差与阈值y的差别,判断该峰是否是呼吸频率,如果是,则进行k倍增强;如果不是,则不对该频点的频率幅度值进行处理;
(5)对Rmf i(f,l)进行呼吸信号相关叠加获得叠加和Rmfs(f,l)。呼吸信号相关叠加示意图如图4所示,图4(a)是单块慢时间—快时间数据变换到频域中某个距离上的频率分布示意图,多块距离—频率数据叠加并归一化后形成图4(b),呼吸信号增强,频域噪声稳态随机分布于整个频域观测空间上。
(6)设置合理阈值和增强比对叠加和进行呼吸信号增强,获得结果Rs(f,l)。由于叠加和呼吸信号较强,阈值应设置为较大值,数据归一化后设置为y=0.8,增强比也应设置较大,设置为k=6;步骤如(4)中的1)~4);
(7)对Rs(f,l)去除噪声均值并进行多重系数衰减的高斯噪声擦除,该方法示意图如图5所示。图5(a)是步骤(5)后的结果,图5(b)中虚线是噪声均值示意图,图5(c)是去除一次高斯噪声均值的结果。该方法步骤如下:
1)获取距离—频域叠加和结果Rs(f,l);
2)由式获得噪声均值Rs(f,l)减去该均值并取绝对值abs为取绝对值函数;
3)擦除系数衰减函数应用余弦函数模型,(j是擦除次数,设置j∈[0 20]),获取当前擦除系数a和高斯噪声均值。将每个距离上的频率幅度减去擦除系数与高斯噪声均值的乘积,取绝对值为频率幅度值如下:
设Rs(f,l)为步骤2)后的结果,设Rs(f,l,j)为第j次擦除后的距离—频率结果,j=0时为原结果,即
Rs(f,l,0)=Rs(f,l)
在第j(j>0)次擦除时,擦除系数为aj,沿频率向的高斯噪声均值为
则有
4)迭代步骤3)20次,最终获得检测结果,如图6所示,人体目标的呼吸频率为0.31Hz,人体目标对雷达的参考距离为3.34m,检测结果正确。
Claims (7)
1.一种基于超宽带雷达的呼吸信号检测算法,其特征在于,所述方法包括以下步骤:
步骤1:超宽带雷达发射冲击脉冲并采样量化回波数据,形成二维慢时间—快时间矩阵r(m,n);
步骤2:对二维矩阵沿慢时间方向应用Motion Filter滤波去除背景杂波获得结果rmf(m,n);
步骤3:对步骤2的结果沿慢时间方向分块为rmf i(m,n),优先分成20~30块,块数设为K,i=1,2,...K;
步骤4:获取各分块频域幅度谱并进行呼吸信号增强,结果为Rmf i(f,l):
对各块二维慢时间—快时间数据沿慢时间方向进行快速傅里叶变换(FFT),并获取二维距离—频率频域幅度谱,对各块频域幅度谱进行呼吸信号增强,步骤如下:
1)设置频率搜索范围F,优先设置F=0.15~0.45Hz,该范围F必须包含呼吸频率f0;
2)设置阈值y和增强比k;
3)扫描每个距离上F范围内的频率幅度值的变化,计算峰值区间内的峰值和PKs并获得峰值区间旁的落差ΔR;
4)找到最大的落差,依据该落差与阈值y的差别,判断该峰是否是呼吸频率,如果是,则进行k倍增强;如果不是,则不对该频点的频率幅度值进行处理;
步骤5:对Rmf i(f,l)进行呼吸信号相关叠加,获得叠加和Rmfs(f,l):
步骤6:对Rmfs(f,l)通过步骤4中的呼吸增强步骤进行呼吸信号增强,获得结果Rs(f,l);
步骤7:对Rs(f,l)去除噪声均值并进行多重系数衰减的高斯噪声擦除过程,步骤如下:
1)获取步骤6后的距离—频域叠加和结果Rs(f,l);
2)由式获得噪声均值Rs(f,l)减去该均值并取绝对值abs()为取绝对值函数;
3)获取当前擦除系数a和每个距离上的高斯噪声均值,将每个距离上的频率幅度减去擦除系数与高斯噪声均值的乘积,取绝对值为频率幅度值;
4)迭代20~30次的步骤3),最终获得检测结果。
2.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤3中,对二维回波数据沿慢时间方向分成多块,相邻两块的距离优先设置为40~70,块长L根据频率分辨率要求设置,PRF是脉冲重复频率,Δf是频率分辨率。
3.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤4和步骤6中,应用呼吸信号增强方法进行单块距离—频率数据呼吸信号增强和叠加和距离—频率数据呼吸信号增强。
4.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤4和步骤6中,优先设置频率搜索范围F=0.15~0.45Hz,在此范围内搜索呼吸信号;峰值搜索方法中,根据频率分辨率设置峰值区间长度并获得峰值区间内的频率采样点的和作为当前峰值和PKs;落差ΔR是峰值和与紧邻峰值区间的两个频点频率幅度值的和DPKs求差,即,ΔR=PKs-DPKs。
5.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤5中,对多块距离—频率数据进行呼吸信号相关叠加,某个距离上的慢时间回波信号r(t)的频域信号R(f)为:
R(f)=A·δ(f-f0)+N(f)
其中,A为呼吸频率信号幅度值,f0为呼吸频率,δ(f)为冲击信号,N(f)为频域噪声,K块数据Rk(f),k=1,2,......,K
Rk(f)=Ak·δ(f-f0)+Nk(f)
则K个距离—频域的Rk(f)叠加和Rs(f)为
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>R</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中当某个距离上没有呼吸信号时Ak=0,则
<mrow>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
当这个距离上有呼吸信号时Ak≠0,则
<mrow>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
Nk(f)是频域随机噪声,由于多块的叠加效应,Ns(f)是稳态随机分布,则有
<mrow>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mrow>
n(f)是高斯噪声,呼吸信号相关叠加有效增强呼吸信号,而在单块呼吸信号增强过程中误放大的随机噪声不会产生叠加干扰。
6.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤7中,多重系数衰减的高斯噪声擦除过程逐步减弱整个频率观测空间中的噪声,期望噪声均值衰减为0,设Rs(f,l)为权利要求1中步骤7的2)后的结果,设Rs(f,l,j)为第j次擦除后的距离—频率结果,j=0时为原结果,即
Rs(f,l,0)=Rs(f,l)
在第j(j>0)次擦除时,擦除系数为aj,沿频率向的高斯噪声均值为
<mrow>
<mover>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>f</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
则有
<mrow>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>a</mi>
<mi>b</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>f</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>(</mo>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
如果某个距离上存在呼吸信号,则该距离上只会保留下呼吸信号。
7.根据权利要求1所述的基于超宽带雷达的呼吸信号检测算法,其特征在于,在步骤7中,擦除系数a随擦除次数x的增多递减,递减函数模型可选择指数模型、对数模型、余弦函数模型,优先选择余弦函数模型a=k·cos(c·j)+b,其中j是擦除次数,j∈[0,J],J优先设置为20~30,a是擦除系数,递减范围优先设置为1.6~0,k、c、b是函数参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710372374.XA CN107167802A (zh) | 2017-05-24 | 2017-05-24 | 一种基于超宽带雷达的呼吸信号检测算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710372374.XA CN107167802A (zh) | 2017-05-24 | 2017-05-24 | 一种基于超宽带雷达的呼吸信号检测算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107167802A true CN107167802A (zh) | 2017-09-15 |
Family
ID=59821977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710372374.XA Pending CN107167802A (zh) | 2017-05-24 | 2017-05-24 | 一种基于超宽带雷达的呼吸信号检测算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107167802A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108888249A (zh) * | 2018-06-07 | 2018-11-27 | 北京邮电大学 | 一种非接触式车内多人生命体征监测的方法及装置 |
CN109009124A (zh) * | 2018-06-05 | 2018-12-18 | 南通大学 | 基于超宽带雷达的呼吸频率测量及目标定位方法 |
CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
CN113720882A (zh) * | 2021-08-30 | 2021-11-30 | 天津大学 | 一种基于mq3传感器和uwb激光雷达的酒驾检测系统及方法 |
CN113786177A (zh) * | 2021-08-24 | 2021-12-14 | 煤炭科学研究总院 | 生命体征信息提取方法、装置及电子设备 |
CN113786175A (zh) * | 2021-08-16 | 2021-12-14 | 煤炭科学研究总院 | 一种基于n次差时强体征信息提取的超宽带雷达回波信号处理方法和装置 |
CN114166334A (zh) * | 2021-11-23 | 2022-03-11 | 中国直升机设计研究所 | 一种非消声风洞旋翼噪声测点的声衰减系数校准方法 |
US11554665B1 (en) | 2022-08-29 | 2023-01-17 | Tianjin University | Method of detecting of driving under influence of alcohol based on MQ3 sensor and ultra wide band radar |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073042A (zh) * | 2010-10-27 | 2011-05-25 | 中国人民解放军第四军医大学 | 一种多通道基于uwb雷达式生命探测仪的多个人体目标二维定位方法 |
CN103110422A (zh) * | 2012-12-18 | 2013-05-22 | 中国人民解放军第四军医大学 | 基于生物雷达检测的呼吸和心跳实时分离方法 |
CN103245976A (zh) * | 2013-05-23 | 2013-08-14 | 中国人民解放军第四军医大学 | 基于uwb生物雷达的人体目标与周围环境结构兼容探测方法 |
CN104133199A (zh) * | 2014-07-08 | 2014-11-05 | 中国科学院电子学研究所 | 用于生命探测雷达的微弱呼吸信号增强方法 |
-
2017
- 2017-05-24 CN CN201710372374.XA patent/CN107167802A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073042A (zh) * | 2010-10-27 | 2011-05-25 | 中国人民解放军第四军医大学 | 一种多通道基于uwb雷达式生命探测仪的多个人体目标二维定位方法 |
CN103110422A (zh) * | 2012-12-18 | 2013-05-22 | 中国人民解放军第四军医大学 | 基于生物雷达检测的呼吸和心跳实时分离方法 |
CN103245976A (zh) * | 2013-05-23 | 2013-08-14 | 中国人民解放军第四军医大学 | 基于uwb生物雷达的人体目标与周围环境结构兼容探测方法 |
CN104133199A (zh) * | 2014-07-08 | 2014-11-05 | 中国科学院电子学研究所 | 用于生命探测雷达的微弱呼吸信号增强方法 |
Non-Patent Citations (2)
Title |
---|
YI KE WANG ETAL: "Supervised Automatic Detection of UWB Ground-Penetrating Radar Targets Using the Regression SSIM Measure", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
李鑫: "UWB脉冲雷达呼吸信号检测算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109009124A (zh) * | 2018-06-05 | 2018-12-18 | 南通大学 | 基于超宽带雷达的呼吸频率测量及目标定位方法 |
CN109009124B (zh) * | 2018-06-05 | 2021-08-06 | 南通大学 | 基于超宽带雷达的呼吸频率测量及目标定位方法 |
CN108888249A (zh) * | 2018-06-07 | 2018-11-27 | 北京邮电大学 | 一种非接触式车内多人生命体征监测的方法及装置 |
CN110507293A (zh) * | 2019-07-26 | 2019-11-29 | 中国电子科技集团公司第三十八研究所 | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 |
CN110554381A (zh) * | 2019-08-30 | 2019-12-10 | 湖南正申科技有限公司 | 一种用于冲激脉冲式穿墙雷达的人体静止目标加速检测方法 |
CN113786175A (zh) * | 2021-08-16 | 2021-12-14 | 煤炭科学研究总院 | 一种基于n次差时强体征信息提取的超宽带雷达回波信号处理方法和装置 |
CN113786177A (zh) * | 2021-08-24 | 2021-12-14 | 煤炭科学研究总院 | 生命体征信息提取方法、装置及电子设备 |
CN113786177B (zh) * | 2021-08-24 | 2023-03-03 | 煤炭科学研究总院有限公司 | 生命体征信息提取方法、装置及电子设备 |
CN113720882A (zh) * | 2021-08-30 | 2021-11-30 | 天津大学 | 一种基于mq3传感器和uwb激光雷达的酒驾检测系统及方法 |
CN114166334A (zh) * | 2021-11-23 | 2022-03-11 | 中国直升机设计研究所 | 一种非消声风洞旋翼噪声测点的声衰减系数校准方法 |
US11554665B1 (en) | 2022-08-29 | 2023-01-17 | Tianjin University | Method of detecting of driving under influence of alcohol based on MQ3 sensor and ultra wide band radar |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167802A (zh) | 一种基于超宽带雷达的呼吸信号检测算法 | |
CN110507293B (zh) | 一种超宽带穿墙雷达人体呼吸及心跳检测方法及系统 | |
CN106874833B (zh) | 一种振动事件的模式识别方法 | |
Guan et al. | Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter | |
CN103760543B (zh) | 一种基于多模式的恒虚警目标检测方法 | |
CN106019254B (zh) | 一种uwb冲击生物雷达多人体目标距离向分离辨识方法 | |
CN106842164B (zh) | 基于小波去噪的非合作脉冲压缩雷达微弱目标检测方法 | |
CN107144829A (zh) | 一种高效的激光雷达回波信号去噪方法 | |
CN103245976B (zh) | 基于uwb生物雷达的人体目标与周围环境结构兼容探测方法 | |
CN106093868B (zh) | 一种基于双源ir-uwb生物雷达的强反射杂波消除方法 | |
CN109031229B (zh) | 一种杂波环境下目标跟踪的概率假设密度方法 | |
CN104215959B (zh) | 一种多机动目标径向初速度和径向加速度的估计方法 | |
CN104793253A (zh) | 基于数学形态学的航空电磁数据去噪方法 | |
CN109633633A (zh) | 一种基于分段分类增强处理的生命信号增强方法 | |
CN103116740B (zh) | 一种水下目标识别方法及其装置 | |
CN103169449B (zh) | 呼吸信号检测方法和装置 | |
CN106199539A (zh) | 基于白化滤波器的地杂波抑制方法 | |
CN104849701B (zh) | 一种车载防撞雷达系统的杂波抑制方法 | |
CN113640768A (zh) | 一种基于小波变换的低分辨雷达目标识别方法 | |
CN106569190B (zh) | 一种高海况下海面目标检测装置及其检测方法 | |
CN103885044B (zh) | 一种基于clean算法的窄带雷达回波杂噪抑制方法 | |
CN107121705A (zh) | 一种基于自动反相校正和峰度值比较的探地雷达回波信号去噪算法 | |
Bagwe et al. | Automated Radar Signal Analysis Based on Deep Learning | |
CN105699971A (zh) | 一种sar雷达运动目标成像方法 | |
CN115345216A (zh) | 一种融合先验信息的fmcw雷达干扰消除方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170915 |
|
WD01 | Invention patent application deemed withdrawn after publication |