一种穿墙雷达墙体参数估计方法技术

一种穿墙雷达墙体参数估计方法，主要解决现有方法时延估计结果精度较差而导致墙体参数估计结果偏差大的问题。通过七个步骤利用正交匹配追踪稀疏重建算法估计墙体前表面和后表面反射波的时延差，将M个观测位置的时延差估计值和时延差理论值构造成目标函数，通过最小化目标函数对墙体的厚度和相对介电常数进行准确估计，在结合单基地收发同置下墙体后表面和前表面反射波幅度比对墙体的电导率进行准确估计。在墙体参数估计过程中采用正交匹配追踪稀疏重建算法对墙体前表面和后表面的双程传输时延差进行估计，提高了低信噪比下时延差估计结果的分辨率和准确度，从而保证了墙体参数估计的准确度。本发明专利技术特别适用于低信噪比情况下对薄层墙体的参数估计。

A Wall Parameter Estimation Method for Wall-Crossing Radar

A wall parameter estimation method for through-wall radar mainly solves the problem of large deviation of wall parameter estimation results caused by poor accuracy of time delay estimation results of existing methods. Seven steps are used to estimate the time delay difference of reflected wave from the front and back surfaces of the wall by using the orthogonal matching pursuit sparse reconstruction algorithm. The estimated time delay difference and theoretical value of time delay difference of M observation positions are constructed as objective functions. The thickness and relative dielectric constant of the wall are accurately estimated by minimizing the objective function. The back and front surfaces of the wall are combined with the single-base receiving and receiving at the same location. The amplitude ratio of reflected wave is used to estimate the conductivity of the wall accurately. In the process of wall parameter estimation, orthogonal matching pursuit sparse reconstruction algorithm is used to estimate the two-way transmission delay difference between the front and back surfaces of the wall, which improves the resolution and accuracy of the estimation result under low signal-to-noise ratio, thus ensuring the accuracy of wall parameter estimation. The invention is especially suitable for parameter estimation of thin wall under low signal-to-noise ratio.

【技术实现步骤摘要】
一种穿墙雷达墙体参数估计方法
本专利技术涉及穿墙雷达
，尤其涉及一种穿墙雷达墙体参数估计方法，该方法主要应用于城市执法、灾害救援和军事行动等领域，特别适用于低信噪比情况下对薄层墙体的参数估计。
技术介绍
穿墙雷达是一种利用电磁波的低频穿透特性对墙后隐蔽目标进行探测的透视成像技术，在穿墙雷达探测过程中，如果墙体参数(介电常数、电导率、墙体厚度)已知,很多成像算法(例如后向投影算法等)容易消除由墙体引起的位置偏移效应。但在实际应用中，墙体参数是无法提前知晓的，墙体参数的估计精度会造成成像质量下降、目标位置定位的偏差以及出现虚假目标等问题。因此，如何对墙体参数进行有效准确估计，是目前穿墙雷达面临的难题之一。针对穿墙雷达的墙体参数估计，已有若干学者对其进行研究并提出了一系列方法。从回波利用情况来看分为两种，一种是从墙体前后表面回波中提取包含墙体参数的信息，通过相关函数公式推算确定墙体参数；另外一种是采集墙后目标回波数据并对其进行成像定位，通过多次图像质量评估、定位修正寻求最优墙体参数。目前，通过测量墙体前后表面反射回波时延差的方法可以估计墙体参数，现有的时延估计方法通常采用基于快速傅里叶变换和子空间超分辨率方法，但这些方法受到墙体的厚度大小和测量数据信噪比的制约，墙体前表面和后表面反射回波的时延估计结果精度较差，从而导致墙体参数估计结果具有很大偏差。
技术实现思路
本专利技术要解决的技术问题是针对上述现有技术的不足，提供一种穿墙雷达墙体参数估计方法，实现在穿墙雷达探测过程中对墙体参数进行准确估计。利用正交匹配追踪稀疏重建算法估计墙体前表面和后表面反射波的时延差，基于M个观测位置的时延差估计值和时延差理论值构造目标函数，通过最小化目标函数实现对墙体的厚度和相对介电常数进行准确估计，在结合单基地收发同置下墙体后表面和前表面反射波幅度比对墙体的电导率进行准确估计。为实现上述目的，本专利技术采用如下技术方案：一种穿墙雷达墙体参数估计方法，包括以下步骤：步骤1：穿墙雷达发射天线和接收天线到墙体前表面的距离均为r，发射天线保持不动，接收天线沿着水平测线方向按固定步长移动M次，共得到M个观测位置，在每个观测位置记录N个均匀频率点的测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维向量Tm＝[Tm(f0),Tm(f1),…,Tm(fN-1)]T，fn＝f0+nΔf为第n(n＝0,1,…,N-1)个频率点的频率，f0为穿墙雷达系统的起始频率，Δf为频率步进间隔；步骤2：将收发天线置于自由空间中，对应每个和步骤1相同的收发天线间距，记录N个均匀频率点的天线直达波测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维向量bm＝[bm(f0),bm(f1),…,bm(fN-1)]T；步骤3：采用背景对消法去除步骤1得到的测量数据中的天线直达波信号，得到墙体回波测量数据。将第m(m＝0,1,…,M-1)个观测位置墙体回波测量数据表示为N×1维向量ym＝Tm-bm。设置最大双程传输时延τmax，将最大双程传输时延τmax均匀划分为Q个时延网格，则可以得到Q×1维双程传输时延向量τ＝[τ0,τ1,…,τQ-1]T。则第m(m＝0,1,…,M-1)个观测位置墙体回波测量数据表示为矩阵向量的形式，如式(1)所示：ym＝Asm+nm(1)其中，ym＝[ym(f0),ym(f1),…,ym(fN-1)]T为第m个观测位置对应的N×1维墙体回波测量数据向量，sm＝[sm(0),sm(1),…,sm(Q-1)]T为Q×1维幅度向量，nm为N×1维测量噪声向量，A＝[a0,a1,…,aQ-1]为N×Q维字典矩阵，其第q(q＝0,1,…,Q-1)列如式(2)所示：步骤4：在第m(m＝0,1,…,M-1)个观测位置，采用正交匹配追踪稀疏重建算法对各个收发天线间距下的墙体回波测量数据进行时延估计，得到墙体前表面和后表面反射波的时延差估计值具体步骤为：①初始化残差rm0＝ym，支撑集Ω0为空集，迭代次数k＝0；②计算残差rmk与字典矩阵的列向量内积中的最大值所对应的索引集，即Λk＝argmaxq{um(q)}，其中相关系数um(q)＝|<rmk,aq>|,q＝0,1,…,Q-1；③更新支撑集Ωk+1＝Ωk∪Λk，计算④更新残差⑤迭代次数k加1，当k＜2时返回步骤②，否则停止迭代；⑥得到墙体前表面和后表面反射波的双程传输时延差估计值，记为步骤5：通过几何模型计算墙体前表面和后表面反射波的理论时延差Δtm(d,εr)，在第m个观测位置，墙体前表面和后表面反射波的理论时延差表示如下：其中2Lm是第m个观测位置发射天线和接收天线的距离，c为电磁波在真空中传播速度，d为墙体厚度，εr为墙体的相对介电常数。xm表示第m个观测位置对应的折射点P的位置，可表示为步骤6：构造目标函数f(d,εr)，得到墙体厚度d和相对介电常数εr的估计值；利用步骤4得到的M个观测位置的时延差估计值和步骤5得到的M个观测位置的时延差理论值Δtm(d,εr)构造目标函数如下：通过求解式(5)所示目标函数的最小值得到墙体厚度d和相对介电常数εr的估计值。步骤7：利用步骤6的求解结果估计墙体的电导率σ，具体方法如下：将收发同置天线放置于距离墙体前表面r处，得到墙体前表面和后表面的反射波幅值R1和R2，因此，墙体后表面和前表面的反射波的幅度比为求解式(6)得到墙体损耗衰减率表达式为将步骤6估计的墙体厚度d和相对介电常数εr带入式(8)，求得墙体损耗衰减率α。针对电磁波损耗较低的墙体，利用下式可以准确计算出墙体的电导率σ其中自由空间波阻抗η0＝120π。采用上述技术方案所产生的有益效果在于：在墙体参数估计过程中采用正交匹配追踪稀疏重建算法对墙体前表面和后表面反射波的双程传输时延差进行估计，显著提高了低信噪比下时延差估计结果的分辨率和准确度，从而保证了墙体参数估计的准确度。本专利技术提供的穿墙雷达墙体参数估计方法，特别适用于低信噪比情况下对薄层墙体的参数估计。附图说明图1为本专利技术实施例提供的一种穿墙雷达墙体参数估计方法的流程图；图2为本专利技术实施例提供的墙体参数反演场景示意图；图3为本专利技术实施例提供的采用本专利技术的穿墙雷达墙体参数估计方法进行墙体参数估计的结果图。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。实施例如图1所示，一种穿墙雷达墙体参数估计方法，该方法是通过下述步骤实现的：步骤1：穿墙雷达发射天线和接收天线到墙体前表面的距离均为r，发射天线保持不动，接收天线沿着水平测线方向按固定步长移动M次，共得到M个观测位置，在每个观测位置记录N个均匀频率点的测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维向量Tm＝[Tm(f0),Tm(f1),…,Tm(fN-1)]T，fn＝f0+nΔf为第n(n＝0,1,…,N-1)个频率点的频率，f0为穿墙雷达系统的起始频率，Δf为频率步进间隔；步骤2：将收发天线置于自由空间中，对应每个和步骤1相同的收发天线间距，记录N个均匀频率点的天线直达波测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维本文档来自技高网...

【技术保护点】
1.一种穿墙雷达墙体参数估计方法，其特征在于：该方法是利用正交匹配追踪稀疏重建算法估计墙体前表面和后表面反射波的时延差，基于M个观测位置的时延差估计值和时延差理论值构造目标函数，通过最小化目标函数实现对墙体的厚度和相对介电常数进行准确估计，在结合单基地收发同置下墙体后表面和前表面反射波幅度比对墙体的电导率进行准确估计并通过具体步骤实现。

【技术特征摘要】
1.一种穿墙雷达墙体参数估计方法，其特征在于：该方法是利用正交匹配追踪稀疏重建算法估计墙体前表面和后表面反射波的时延差，基于M个观测位置的时延差估计值和时延差理论值构造目标函数，通过最小化目标函数实现对墙体的厚度和相对介电常数进行准确估计，在结合单基地收发同置下墙体后表面和前表面反射波幅度比对墙体的电导率进行准确估计并通过具体步骤实现。2.如权利要求1所述的一种穿墙雷达墙体参数估计方法，其特征在于该方法是通过下述步骤实现的：步骤1：穿墙雷达发射天线和接收天线到墙体前表面的距离均为r，发射天线保持不动，接收天线沿着水平测线方向按固定步长移动M次，共得到M个观测位置，在每个观测位置记录N个均匀频率点的测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维向量Tm＝[Tm(f0),Tm(f1),…,Tm(fN-1)]T，fn＝f0+nΔf为第n(n＝0,1,…,N-1)个频率点的频率，f0为穿墙雷达系统的起始频率，Δf为频率步进间隔；步骤2：将收发天线置于自由空间中，对应每个和步骤1相同的收发天线间距，记录N个均匀频率点的天线直达波测量数据，将第m(m＝0,1,…,M-1)个观测位置测量数据表示为N×1维向量bm＝[bm(f0),bm(f1),…,bm(fN-1)]T；步骤3：采用背景对消法去除步骤1得到的测量数据中的天线直达波信号，得到墙体回波测量数据，将第m(m＝0,1,…,M-1)个观测位置墙体回波测量数据表示为N×1维向量ym＝Tm-bm，设置最大双程传输时延τmax，将最大双程传输时延τmax均匀划分为Q个时延网格，则可以得到Q×1维双程传输时延向量τ＝[τ0,τ1,…,τQ-1]T，则第m(m＝0,1,…,M-1)个观测位置墙体回波测量数据表示为矩阵向量的形式，如式(1)所示：ym＝Asm+nm(1)其中，ym＝[ym(f0),ym(f1),…,ym(fN-1)]T为第m个观测位置对应的N×1维墙体回波测量数据向量，fn＝f0+nΔf为第n(n＝0,1,…,N-1)个频率点的频率，f0为穿墙雷达系统的起始频率，Δf为步进频率，sm＝[sm(0),sm(1),…,sm(Q-1)]T为Q×1维幅度向量，nm为N×1维测量噪声向量，A＝[a0,a1,…,aQ-1]为N×Q维字典矩阵，其第q(q＝0,1,…,Q-1)列如式(2)所示：步骤4：在第m(m＝0,1,…,M-1)个观测位置，采用正交匹配追踪稀疏重建算法对各个收发天线间距下的墙体回波测量数据进行时延估计，得到墙体前表面和后表面反射波的时延差估计值具体步骤为：①初始化残差rm0＝ym，支撑集Ω0为空集，迭代次数k＝0；②计算残差rmk与字...

【专利技术属性】
技术研发人员：屈乐乐，方中里，杨天虹，孙延鹏，
申请(专利权)人：沈阳航空航天大学，
类型：发明
国别省市：辽宁,21