基于压缩传输的TDOA估计方法技术

本发明专利技术公开了一种基于压缩传输的TDOA估计方法，主要解决现有定位方法依赖于原信号重构且运算效率低，在带宽有限和低信噪比的条件下稳健性差和估计效率低的问题。其方案是：通过分组压缩的思想，设计压缩矩阵，对接收设备获取到的目标信号进行信号压缩；对压缩后的信号进行频域变换，得到频域信号和转换矩阵；由频域信号和转换矩阵，计算DFT系数算子和延迟矩阵；通过DFT系数算子和延迟矩阵，构建最小化差异函数，并求解该函数得到目标信号延迟的估计值。本发明专利技术能避免重构原信号的冗余操作，减少了传输带宽占用率，提高了估计效率和稳健性。可用于空中运动中的飞机、卫星、无人机等目标进行无源定位。

TDOA estimation method based on compressed transmission

The invention discloses a TDOA estimation method based on compressed transmission, which mainly solves the problems of poor robustness and low estimation efficiency under the condition of limited bandwidth and low signal-to-noise ratio. The scheme is as follows: compression matrix is designed to compress the target signal acquired by the receiving equipment through the idea of group compression; frequency domain transformation is carried out for the compressed signal to obtain the frequency domain signal and conversion matrix; DFT coefficient operator and delay matrix are calculated from the frequency domain signal and conversion matrix; DFT coefficient operator and delay matrix are constructed to minimize the difference function, and The estimation of target signal delay is obtained by solving the function. The invention can avoid the redundant operation of reconstructing the original signal, reduce the occupancy rate of transmission bandwidth, and improve the estimation efficiency and robustness. It can be used for passive positioning of aircraft, satellite, UAV and other targets in air motion.

【技术实现步骤摘要】
基于压缩传输的TDOA估计方法
本专利技术属于通信
，更进一步涉及一种TDOA估计方法，可用于空中运动中的飞机、卫星、无人机等目标进行无源定位。
技术介绍
在多站无源定位过程中，多站的天线会截获地面设置的雷达或者通信设备等目标发出的电磁信号。到达时间差TDOA定位是一种利用时间差进行多站定位的方法。通过测量信号到达监测站的时间差，可以确定目标的距离。传统TDOA定位方法通常需要将奈奎斯特采样信号完整地传输回中心站，在带宽有限的条件下，将会占用较大部分的带宽。为了减少带宽占用，引入信号压缩的思想，压缩传输过程中的信号。由于压缩传输不满足时间等间距性，所以传统的方法在计算目标信号时延差前一般需要重构压缩前的原始信号，这将大大增加了计算的复杂程度。陈墩金等在2016年申请的专利“一种基于压缩感知和近邻差值比较的室内定位方法”中提出，首先构建室内网格地图,通过锚节点采集未知节点信号强度,将采集到的信号强度值构成压缩采样向量,再通过基追踪重构算法准确恢复出原信号,最后使用RSSI差值比较法,获取未知节点的精确定位。田子建等在2017年申请的专利“基于压缩感知子空间重构的TOA矿井目标定位方法”中，提出通过在基于TOA矿井目标定位方法中引入压缩感知子空间重构,使得感知节点采样数据量减少的情况下,也能够进行目标的精确定位。上述两种方法在定位前都需要通过重构算法恢复出原信号，即引入了信号重构的步骤，增加了系统的计算量和复杂度。
技术实现思路
本专利技术的目的在于针对上述现有技术的不足，提出一种基于压缩传输的TDOA估计方法，以降低带宽占用率，减少TDOA定位的计算量和复杂度，提高了对目标进行定位的效率。本专利技术的技术方案是，通过设计压缩矩阵，对观测平台获取到的信号，进行压缩处理，然后再通过数据链路传输到中心站；中心站通过构建压缩信号的频域离散傅里叶变换模型和求解时延的最小差异函数进行数据处理；最后通过波谷搜索，找到匹配的时延信息。其实现步骤包括如下：(1)设计压缩矩阵，对接收设备获取到的目标信号进行信号压缩：(1a)利用分组压缩的思想，推导出压缩矩阵A，该压缩矩阵A的第m行第k列元素值Am,k表示如下：式中，qm＝rand(Q)-1，M表示压缩后信号序列长度，N表示压缩前接收机接收信号长度，定义压缩倍数Q＝N/M，rand(Q)是随机函数，能够生成1～Q之间的随机数字；(1b)通过压缩矩阵A对目标信号si进行时域压缩，得到压缩信号yi；(2)对压缩信号yi进行频域离散傅里叶变换DFT，得到频域信号Yi和转换矩阵H；(3)根据频域信号Yi和转换矩阵H估计接收信号si的时延τ：(3a)根据步骤(2)得出的Yi和H，得到DFT系数DPC算子U(p)和延迟矩阵g(p,τ)：式中，Y1(p)是Yi的第p个元素，hpk是H第p行第k列的元素，||表示求模函数，*表示取共轭操作，τ表示目标信号时延；(3b)由步骤(3a)的结果构建两者的最小化差异函数F(τ)：式中，为信号si方差的估计值，||||表示向量求模函数，argmin表示函数取值最小时的自变量取值；3c)对最小化差异函数F(τ)进行求解，得出时延估计值本专利技术与现有技术相比具有以下优点：第一，本专利技术在传感器接收到获取目标信号后，利用信号压缩原理，将目标信号进行压缩传输，从而降低带宽的占用率和计算的复杂度。第二，本专利技术在TDOA估计过程中，采用PDC非重构估计方法，避免了现有技术在对压缩后信号进行TDOA估计时的重构步骤，降低了计算复杂度，提高了TDOA的估计效率。附图说明图1是本专利技术的使用场景图；图2是本专利技术的实现流程图；图3是本专利技术在信噪比SNR＝10db的条件下的时延估计仿真结果图；图4是本专利技术在信噪比SNR＝20db的条件下的时延估计仿真结果图；图5是本专利技术时延估计的信噪比SNR与均方误差MSE的关系图。具体实施方式下面结合附图对本专利技术做进一步的描述。参照图1，对本专利技术的使用场景做进一步的说明。如图1所示，该使用场景大致分成三个部分：目标反射基站信号、接收信号压缩传输和中心站处理信号。本专利技术针对其中两个接收设备，利用一对接收设备接收的目标信号，进行数据压缩，经过数据链路传输到中心站进行数据处理和时延估计，并通过测量信号到达监测站的时延差，确定目标的距离。参照图2，对本专利技术的具体实现步骤如下。步骤1，设计压缩矩阵，对接收设备获取到的目标信号进行信号压缩。通常压缩矩阵可以设置哈达码矩阵，但是在实际应用中，对信号处理比较复杂。本专利技术中，引入了分组压缩的思想将压缩矩阵进行重新分配。(1a)利用分组压缩的思想，推导出压缩矩阵A，该压缩矩阵A的第m行第k列元素值Am,k表示如下：式中，qm＝rand(Q)-1，M表示压缩后信号序列长度，N表示压缩前接收机接收信号长度，定义压缩倍数Q＝N/M，rand(Q)是随机函数，能够生成1～Q之间的随机数字；(1b)接收设备获取原始信号后，分别对第i个接收机接收到的目标信号si以压缩矩阵A进行时域压缩,即两者相乘，得到序列长为M的接收信号yi：yi＝Asi+wi，其中i＝1,2，i表示接收机的序号，yi表示压缩后的信号，A表示压缩矩阵，si表示一对接收机接收到的信号，其信号之间的关系满足s2(n)＝s1(n+τ)，s1(n)表示目标信号s1的第n个元素，s2(n+τ)表示目标信号s2的第n+τ个元素，τ表示信号时延，wi表示在压缩过程中引入的噪声信号。步骤2，对压缩后信号yi进行频域变换，得到频域信号Yi和转换矩阵H。频域变换的方法包含连续傅里叶变换，离散傅里叶变换，快速傅里叶变换和傅里叶级数等变换方法。对于离散的信号，本专利技术中，采用离散傅里叶变换中的乘子法获取其频域形式，其实现步骤如下：(2a)根据离散傅里叶变换乘子公式定义,计算出目标信号和压缩信号的DFT矩阵乘子F和G；其中F的第n行第k列的元素值Fn,k和G的第m行第p列的元素值Gm,p分别表示如下：式中，e表示以自然常数为底的指数函数,j表示虚部符号。(2b)将离散傅里叶变换矩阵与压缩信号相乘后，得到频域信号。Yi＝Gyi＝GAsi+Gwi，＝GAF-1Si+Gwi，式中，Yi表示频域中的第i个信号矢量，G表示傅里叶变换矩阵，yi表示时域中的第i个压缩信号，A表示压缩矩阵，si表示一对接收机接收到的信号，wi表示在压缩过程中引入的噪声信号。定义转换矩阵H:由H得到Yi的第p列Yi(p)：式中，hpk是H的第p行第k列的元素。考虑时延影响，由时域s2(n)＝s1(n+τ)的关系式，推导出两个频域信号之间的关系式:步骤3，计算DFT系数算子和延迟矩阵。根据步骤(2)得出的Yi和H，得到DFT系数算子U(p)和延迟矩阵g(p,τ)：式中，Y1(p)是Yi的第p个元素，*表示取共轭操作，||表示求模函数，hpk是H的第p行第k列元素，UR(p)和UI(p)分别是U(p)实部和虚部，gR(p,τ)和gI(p,τ)分别是g(p,τ)实部和虚部。步骤4，通过离散傅里叶变换乘子法，得到最小化差异函数，并求解函数估计接收信号的时延。(4a)由步骤(3)的结果构建最小化差异模型F(τ)：式中，||||表示向量求模函数，argmin表示函数取值最小时的自变量取值；为信号si的方差估计，其表示如下：(本文档来自技高网...

【技术保护点】
1.一种基于压缩传输的TDOA估计方法，其特征在于，包括：(1)设计压缩矩阵，对接收设备获取到的目标信号进行信号压缩：(1a)利用分组压缩的思想，推导出压缩矩阵A，该压缩矩阵A的第m行第k列元素值Am,k表示如下：

【技术特征摘要】
1.一种基于压缩传输的TDOA估计方法，其特征在于，包括：(1)设计压缩矩阵，对接收设备获取到的目标信号进行信号压缩：(1a)利用分组压缩的思想，推导出压缩矩阵A，该压缩矩阵A的第m行第k列元素值Am,k表示如下：式中，qm＝rand(Q)-1，M表示压缩后信号序列长度，N表示压缩前接收机接收信号长度，定义压缩倍数Q＝N/M，rand(Q)是随机函数，能够生成1～Q之间的随机数字；(1b)通过压缩矩阵A对目标信号si进行时域压缩，得到压缩信号yi；(2)对压缩信号yi进行频域离散傅里叶变换DFT，得到频域信号Yi和转换矩阵H；(3)根据频域信号Yi和转换矩阵H构建DFT系数乘积(PDC)模型，并求解得到估计接收信号si的时延τ：(3a)根据步骤(2)得出的Yi和H，得到DFT系数算子U(p)和延迟矩阵g(p,τ)：式中，Y1(p)是Yi的第p个元素，hpk是H第p行第k列的元素，||表示求模函数，*表示取共轭操作，τ表示目标信号时延；(3b)由步骤(3a)的结果构建两者的最小化差异函数F(τ)：式中，σ2为信号si的方差，||||表示向量求模函数，argmin表示函数取值最小时的自变量取值；3c)对最小化差异函数F(τ)进行求解，得出时延估计值2.根据权利要求1所述的方法，其中(1b)中得到的压缩信号yi表示如下...

【专利技术属性】
技术研发人员：鲍丹，蔡晶晶，钱奔奔，武斌，秦国栋，刘高高，李鹏，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61