基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法
【专利摘要】本发明公开了一种基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法，包括以下步骤：S01、在得到第n次探测回波信号Pn(M)时，用Pn(M)替换掉上一个滑窗Sa(M,N)内最原始的探测回波信号Pa(M)，得到新的滑窗Sa+1(M,N)，其中滑窗Sa(M,N)包括N次探测回波信号，即滑窗宽度为N，每个探测回波信号包括M个距离单元，且n＞N；S02、对滑窗Sa+1(M,N)内的全部距离单元做N点傅立叶变换，得到滑窗Sa+1(M,N)的慢时间维频谱矩阵Fa+1(M,N)。本发明的基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法具有运算速度快、频谱构建效率高等优点。
【专利说明】
基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法
技术领域
[0001] 本发明主要涉及超宽带雷达探测技术领域，特指一种基于迭代滑窗傅里叶变换的 雷达慢时间维频谱构建方法。
【背景技术】
[0002] 超宽带雷达具有很高的分辨率以及良好的障碍穿透能力，广泛地应用于灾害救 援、反恐侦察等领域中。在使用超宽带雷达对隐蔽的受困人体目标(如废墟掩埋人员、被劫 持人质等)进行探测时，需要积累多次探测回波信号，并在不同的距离单元构建慢时间维频 谱，提取呼吸、心跳、微动等频率，以判断是否存在受困人员目标，确定受困人员位置。而目 前采用传统的基于傅里叶变换获得慢时间维频谱的方法，对于每次探测回波信号均需要进 行大量的计算，计算过程繁琐，费时较长。
【发明内容】
[0003] 本发明要解决的技术问题就在于:针对现有技术存在的技术问题，本发明提供一 种运算速度快、频谱构建效率高的基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方 法。
[0004] 为解决上述技术问题，本发明提出的技术方案为：
[0005] S01、在得到第η次探测回波信号?"(1〇时，用Pn(M)替换掉上一个滑窗S a(M，N)内最 原始的探测回波信号Pa(M)，得到新的滑窗Sa+1(M，N)，其中滑Φ5 3(Μ，Ν)包括N次探测回波信 号，即滑窗宽度为Ν，每个探测回波信号包括Μ个距离单元，且η>Ν;
[0006] S02、对滑窗Sa+1(M，N)内的全部距离单元做Ν点傅立叶变换，得到滑窗S a+1(M，N)的 慢时间维频谱矩阵Fa+1 (Μ，N)。
[0007] 优选地，在步骤S02中，依次对滑窗Sa+1(M，N)的全部距离单元按照公式（1)做N点傅 立叶变换，其中公式(1)为：
[0008] Fa+i(m，i) = [Fa(m，i)-Pa(m)]*a(i)+Pn(m)*P(i) (1)
[0009] 其中i = l，2,…，N;m表示第m个距离单元;Fa(m，i)为Sa(M，N)中第m个距离单元的慢 时间维频谱;P a(m)为最原始的探测回波信号中第m个距离单元;α(Ν)和β(Ν)为矢量更新系 数，
[0010] α(?) =θχρ{2πj(i-1 )/N},
[0011] 0(i)=exp{-2Jij(i-l)(N-l)/N} (2)
[0012] 其中j表示虚数单位。
[0013] 优选地，在步骤SOI中，将Pn(M)按照公式(3)更新滑内的探测回波信号， 得到滑窗S a+1(M，N)，其中公式(3)为：
[0014] Pn(M)=Pn%N(M) (3)
[0015] 其中n>N，且从N+1依次递增，％表示取余运算，N为滑窗宽度。
[0016]优选地，初始滑窗Si(M，N)的慢时间维频谱矩阵Fi(M，N)的产生过程为:对第一次累 积满的滑窗Si(M，N)中的全部距离单元，沿着慢时间维做N点快速傅里叶变换，得到慢时间 维频谱矩阵Fi(M，N)，
[0018] 与现有技术相比，本发明的优点在于：
[0019] 本发明的基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法，用每次新获得 的当前探测回波信号更新滑窗内最原始的探测回波信号，当前对应的慢时间维频谱则可在 原有的慢时间维频谱的基础上快速进行迭代更新，从而使频谱构建效率更高。
【附图说明】
[0020] 图1为本发明的方法流程图。
【具体实施方式】
[0021] 以下结合说明书附图和具体实施例对本发明作进一步描述。
[0022] 如图1所示，本实施例的基于迭代滑窗傅里叶变换的超宽带雷达慢时间维频谱构 建方法，包括以下步骤：
[0023] S01、在得到第η次探测回波信号?"(1〇时，用Pn(M)替换掉上一个滑窗S a(M，N)内最 原始的探测回波信号Pa(M)，得到新的滑窗Sa+1(M，N)，其中滑Φ53(Μ，Ν)包括N次探测回波信 号，即滑窗宽度为Ν，每个探测回波信号包括Μ个距离单元，且η>Ν;
[0024] S02、对滑窗Sa+1(M，N)内的全部距离单元做Ν点傅立叶变换，得到滑窗S a+1(M，N)的 慢时间维频谱矩阵Fa+1 (Μ，N)。
[0025]本实施例中，在步骤S02中，依次对滑窗Sa+1(M，N)的全部距离单元按照公式（1)做N 点傅立叶变换，其中公式(1)为：
[0026] Fa+i(m，i) = [Fa(m，i)-Pa(m)]*a(i)+Pn(m)*0(i) (1)
[0027] 其中i = l，2,…，N;m表示第m个距离单元;Fa(m，i)为Sa(M，N)中第m个距离单元的慢 时间维频谱;Pa(m)为最原始的探测回波信号中第m个距离单元;a(N)和β(Ν)为矢量更新系 数，
[0028] α(?) =θχρ{2πj(i-1 )/N},
[0029] P(i)=exp{-23ij(i-l)(N-l)/N} (2)
[0030] 其中j表示虚数单位。可以看出，当前滑窗的慢时间维频谱得到后，后续新探测回 波信号对应更新后的滑窗，其慢时间维频谱即可根据上式进行快速实时迭代更新，其运算 速度快，频谱构建效率高，对应的探测效率更高。
[0031] 本实施例中，在步骤SO 1中，将Pn(M)按照公式(3)更新滑，N)内的探测回波信 号，得到滑窗Sa+1(M，N)，其中公式(3)为：
[0032] Pn(M)=Pn%N(M) (3)
[0033] 其中n>N，且从N+1依次递增，％表示取余运算，N为滑窗宽度。
[0034]本实施例中，初始滑窗Si (Μ，N)的慢时间维频谱矩阵Fi (Μ，N)的产生过程为:对第一 次累积满的初始滑窗SKM，N)中的全部距离单元，沿着慢时间维做N点快速傅里叶变换，得 到慢时间维频谱矩阵Fi(Μ，N)，
[0036] 在具体应用时，当得到第N次探测回波信号时（设定N为2的整数次幂），滑窗积累 满，用矩阵31(11少)表示第一次累积满的滑窗，包含探测回波信号？ 1少2，？3，一，？~;依次对51 (Μ，N)中全部距离单元，沿着慢时间维做N点快速傅里叶变换，得到慢时间维频谱矩阵Fi (M， N)；
[0037] 对于第N+1次探测回波信号PN+1(M)，用PN+1(M)替换滑Φ5ΚΜ，Ν)中的探测回波信号 Pi(Μ)，得到滑窗S2(Μ，Ν)，包括探测回波信号ΡΝ+1，Ρ2，Ρ 3，…，ΡΝ;
[0038] 对于滑窗32(1，《中第m个距离单元，按照公式(5)，计算得到距离单元m的Ν点傅里 叶变换，构建慢时间维频谱，其中公式(5)如下：
[0039]
[0040] 其中 i = l，2,…，Ν。
[0041]依次对滑窗S2(M，N)的全部距离单元，按照公式(5)，计算每个距离单元的N点傅立 叶变换，即可得到滑窗&(11少)的慢时间维频谱矩阵f2(m，n);
[0042]当得到后续的糾2，糾3，糾4，~次探测回波时，重复上述操作，即可快速更新构建 慢时间维频谱矩阵。
[0043]以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例， 凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的 普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护 范围。
【主权项】
1. 一种基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法，其特征在于，包括以 下步骤： 501、 在得到第η次探测回波信号?"(1〇时，用Pn(M)替换掉上一个滑窗Sa(M，N)内最原始 的探测回波信号?^)，得到新的滑窗S a+1(M，N)，其中滑Φ53(Μ，Ν)包括N次探测回波信号，即 滑窗宽度为Ν，每个探测回波信号包括Μ个距离单元，且η>Ν; 502、 对滑窗Sa+1(M，N)内的全部距离单元做Ν点傅立叶变换，得到滑窗Sa+1(M，N)的慢时 间维频谱矩阵F a+1(M，N)。2. 根据权利要求1所述的基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法，其 特征在于，在步骤S02中，依次对滑窗Sa+1(M，N)的全部距离单元按照公式（1)做N点傅立叶变 换，其中公式(1)为： Fa+i(m，i) = [Fa(m，i)-Pa(m)]*a(i)+Pn(m)*P(i) (1) 其中i = l，2,…，N;m表示第m个距离单元;Fa(m，i)为Sa(M，N)中第m个距离单元的慢时间 维频谱;Pa(m)为最原始的探测回波信号中第m个距离单元;a(N)和β(Ν)为矢量更新系数， a(i) = exp{23ij(i-l)/N}, 0(i) = exp{-2Jij(i-l)(N-l)/N} (2) 其中j表示虚数单位。3. 根据权利要求2所述的基于迭代滑窗傅里叶变换的雷达慢时间维频谱构建方法，其 特征在于，在步骤SOI中，将Pn(M)按照公式(3)更新滑内的探测回波信号，得到滑 窗S a+1(M，N)，其中公式(3)为： Pn(M)=Pn%N(M) (3) 其中n>N，且从N+1依次递增，％表示取余运算，N为滑窗宽度。4. 根据权利要求1至3中任意一项所述的基于迭代滑窗傅里叶变换的雷达慢时间维频 谱构建方法，其特征在于，初始滑窗Si (Μ，N)的慢时间维频谱矩阵Fi (Μ，N)的产生过程为:对 第一次累积满的滑窗Si(M，N)中的全部距离单元，沿着慢时间维做N点快速傅里叶变换，得 到慢时间维频谱矩阵F:(Μ，N)，
【文档编号】G01S7/41GK105866761SQ201610479118
【公开日】2016年8月17日
【申请日】2016年6月27日
【发明人】韩明华, 王生水, 衣晓飞
【申请人】湖南华诺星空电子技术有限公司