一种基于频域TOEPLITZ矩阵重构的跳频信号快速测向方法与流程

一种基于频域toeplitz矩阵重构的跳频信号快速测向方法
技术领域
1.本发明涉及阵列信号处理技术领域，主要涉及一种基于频域toeplitz矩阵重构的跳频信号快速测向方法。


背景技术：

2.无人机起源于军事领域，随着近年来相关技术的进步，民用无人机快速普及，已经在医疗、娱乐、运输、教育等众多领域得到了广泛的应用。但无人机在为人民生活带来了诸多便利的同时，也带来一些新的社会问题。“黑飞”，是指民用无人机未经报备在禁飞区飞行的现象。随着民用无人机普及，国际上无人机“黑飞”相关的安全事故层出不穷，其中部分事件带来了相当严重的后果，肇事“黑飞手”却不知所踪。
3.无线电信号传输距离远，受不良环境影响小，民用无人机遥控器发出的遥控信号中包含无人机型号和“黑飞手”位置信息，信号波达方向(direction of arrival，doa)是其重要的特征参数，对其进行截获处理，有望实现“黑飞手”定位，为“黑飞”现象治理提供有力支撑。但民用无人机遥控器发出的无线电遥控信号属于典型的跳频信号，该信号起源于军事用途，专门为反敌方截获侦察设计，因此对民用无人机遥控信号进行截获和信号处理的挑战很大。主要有以下几方面困难：
4.1)使用开放频段通信，与众多信号混叠：民用无人机遥控信号使用ism开放频段通信，该频段允许任何人发射一定功率信号，常见的wifi、蓝牙和无人机图传信号也都属于该频段，监测“黑飞”无人机遥控信号时面对的干扰信号很多；
5.2)信号微弱，截获信号信噪比低：受开放频段协议限制，民用无人机遥控信号功率非常低，接收信号的信噪很低；
6.3)信号载频跳变速度快，难以压缩采样：民用无人机遥控信号属于典型的跳频信号，信号载频跳变速度很快，平稳时间仅有1-3ms，无法使用压缩采样技术；
7.4)信号载频跳变范围广，数据采样存储压力大：民用无人机遥控信号载频跳变范围达到80mhz，无法压缩采样情况下要求设备具有很高的采样速率，也需要能够相适应的数据存储速率。
8.民用无人机遥控信号功率非常小，而且通信频段电磁环境非常复杂，采集信号中包含无人机图传信号、wifi信号等干扰信号。遥控信号是个跳频信号，载频跳变速度非常快，1-3ms跳变一次，跳变范围非常广，达到80mhz。信道在时域、频域、空域都和wifi 信号、无人机图传信号存在混叠。基于无线电技术监测“黑飞”民用无人机虽然具有很多优势，但是也存在很大挑战。


技术实现要素：

9.发明目的：本发明提供了一种基于频域toeplitz矩阵重构的跳频信号快速测向方法，使用天线阵列截获“黑飞”无人机遥控器发出的遥控信号，然后运用阵列信号处理技术，对数据进行空时频多域处理，进而对跳频信号进行测向，找到“黑飞手”，为有关部门治理“黑飞”问题提供有力的技术支撑。在保证对跳频信号测向精度的同时，降低计算复杂度，减少硬件成本。
10.技术方案：为实现上述目的，本发明采用的技术方案为：
11.一种基于频域toeplitz矩阵重构的跳频信号快速测向方法，包括以下步骤：
12.步骤s1、建立均匀线阵以及对应的接收信号数据模型；
13.步骤s2、通过快速傅里叶变换fft获取频域数据模型，并构建谱峰搜索函数，找到 k个峰值并提取对应的横坐标，作为信号的精准瞬时频率估计值；
14.步骤s3、基于提取的所述峰值对应的k组频域数据，利用协方差矩阵的toeplitz 矩阵特性，构建基于频域数据的协方差矩阵；
15.步骤s4、基于所述协方差矩阵的特征分解，计算噪声子空间e
nk
，并根据root
‑ꢀ
music原理构造多项式求解，获取最接近单位圆的根
16.步骤s5、找到k个最接近于单位圆上的根作为最终解，并基于每个根对应的信号频率解出最终的方向角。
17.进一步地，所述步骤s1中建立的均匀线阵包括m个排列成一条直线的天线阵元，阵元间距d满足d0＝λ/2，其中λ为入射信号的最大载波波长；k个入射信号的入射角满足：θk∈[-π/2,π/2],k＝1,2,
…
,k，频率依次为f1,f2,
…
,fk；则所述均匀线阵在t时刻接收k个跳频信号的接收信号数据模型如下：
[0018]
x(t)＝as(t) n(t)
[0019]
其中x为m
×
1的矩阵，a＝[a(θ1),a(θ2),
…
,a(θk)]，代表方向矩阵，且满足：
[0020][0021]
j代表虚数符号；d0表示均匀线阵中两个阵元之间的间距；c代表光速；代表信号向量，代表噪声向量。
[0022]
进一步地，所述步骤s2中首先将每个天线阵元接收到的信号数据重新排列获得：
[0023][0024]
将时域数据x进行快速傅里叶变换fft，获取当前各个接收信号的实时频率值，即频域数据：
[0025][0026]
q＝1,2,
…
,j,m＝1,2,
…
,m
[0027]
其中，xm(tj)为xm的第j个元素，构造谱峰搜索函数如下：
[0028][0029]
由于不同的信号不共享同一瞬时频率，可以找到k个峰值的横坐标f1,f2,
…
,fk作为信号的精确瞬时频率估计值。
[0030]
进一步地，步骤s3中所述协方差矩阵由m个频域数据利用矩阵的toeplitz性质构造而成；由于其中是固定的相位差； |vk|是固定的幅值；qk,k＝1,2,
…
,k，代表各谱峰值的横坐标索引，由toeplitz性质构建协方差矩阵如下：
[0031][0032]
其中qk是的相位。
[0033]
进一步地，步骤s4中，首先对协方差矩阵进行特征分解如下：
[0034][0035]
得到信号的噪声子空间e
nk
，其中d
sk
表示的最大特征值，d
nk
表示包含其他特征值的对角矩阵；e
sk
代表信号子空间对应的最大特征值的特征向量，噪声子空间e
nk
由对应的其他特征向量组成；
[0036]
然后根据root-music原理构造多项式如下：
[0037][0038]
其中a(z)＝[1,z1,
…
,z
m-1
]为方向向量；获取最接近单位圆的根：
[0039][0040]
即为所需根
[0041]
进一步地，，步骤s5中最终方向角的估计结果如下：
[0042][0043]
有益效果：本发明具备以下优点：
[0044]
(1)、具备强大的抗噪声抗干扰能力；本发明运用阵列信号处理技术，充分利用了天线阵列的空间增益，显著提高了接收信号信噪比。针对干扰问题，本作品使用时-空
‑ꢀ
频多域特征提取，首先确定各个信号的瞬时频率，再对这些频率上的频域数据进行处理，大大提高了数据的可靠性。
[0045]
(2)、低复杂度和硬件成本下的高精度测向能力；本发明所需阵列天线数量和快拍数少，大大降低了硬件成本，另外算法中不需要谱峰搜索等大计算量的步骤，使用频域数据和toeplitz特性直接构造出协方差矩阵，可以满足工程中的实时测向需求。
附图说明
[0046]
图1是本发明电磁矢量l阵结构示意图。
[0047]
图2是本发明方法50次独立估计得到doa参数的散点图。
[0048]
图3是在相同条件下，本发明方法与其他四种方法的doa估计均方根误差随着快拍数变化的性能对比图。
[0049]
图4是在相同条件下，本发明方法与其他四种方法的doa估计均方根误差随着信噪比变化的性能对比图。
[0050]
图5a是本发明实施例中采用的软件无线电平台示意图；
[0051]
图5b是本发明实施例中采用的可以发射跳频信号的无人机遥控器示意图；
[0052]
图6是本发明实施例中软件显示界面示意图。
具体实施方式
[0053]
下面结合附图对本发明作更进一步的说明。所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
[0054]
本发明中用(
·
)
t
表示转置，(
·
)
*
表示共轭，(
·
)h表示共轭转置，大写字母如x表示矩阵，小写字母如x(
·
)表示矢量，表示kronecker积，
⊙
表示khatri-rao积，angle(
·
)表示取复数的相角。
[0055]
本发明中涉及的电磁矢量l阵结构如图1所示，均匀线阵由m个天线排列成一条直线组成，阵元间距均为d0＝λ/2，λ为入射信号的最大载波波长。由于阵元间距没有大于最大波长信号的一半，因此不会产生相位模糊的问题，可以直接通过超分辨率的算法求出结果。然而由于跳频信号持续时间短，频率跳变快，普通的算法并不能很好地计算出估计结果，因为这些算法普遍来说需要很多的快拍数，且运算时间长，在频率高速跳变的场景下并不适用。本发明均匀线阵中基于频域toeplitz矩阵重构的跳频信号快速测向方法的具体实现如下：
[0056]
步骤s1、建立均匀线阵以及对应的接收信号数据模型。
[0057]
假设有k,(k《m)个信号入射到该均匀线阵上，入射角为θk∈[-π/2,π/2],k＝ 1,2,
…
,k，频率依次为f1,f2,
…
,fk；该线阵在t时刻接收k个跳频信号的接收数据模型为：
[0058]
x(t)＝as(t) n(t)
[0059]
其中x为m
×
1矩阵，a＝[a(θ1),a(θ2),
…
,a(θk)]，代表方向矩阵，且满足：
[0060][0061]
j代表虚数符号；d0表示均匀线阵中两个阵元之间的间距；c代表光速；代表信号向量，代表噪声向量。
[0062]
考虑一组由j个快拍组成的连续时间离散信号，其总的接收信号数据模型可以合并为：
[0063]
x＝as n
[0064]
其中x＝[x(t1),x(t2),
…
,x(tj)]，s＝[s(t1),s(t2),
…
,s(tj)]，n＝[n(t1),n(t2),
…
,n(tj)]。
[0065]
步骤s2、通过快速傅里叶变换fft获取频域数据模型，并构建谱峰搜索函数，找到k个峰值并提取对应的横坐标，作为信号的精准瞬时频率估计值。
[0066]
首先对时域数据按每个阵元接收到的数据重新排列得到进行快速傅里叶变化获得当前各个信号的实时频率值
[0067]
[0068]
q＝1,2,
…
,j,m＝1,2,
…
,m
[0069]
其中，xm(tj)为xm的第j个元素。因此，相同元素个数的频域数据可由如下公式表示：
[0070][0071]
即为信号的离散频谱。对每一列的幅度求和以获得频谱，其定义为
[0072][0073]
p即为谱峰搜索函数。假设不同的信号不共享相同的瞬时频率，我们可以找到横坐标f1,f2,
…
,fk为信号的精确频率估计值的k个峰值。
[0074]
步骤s3、基于提取的所述峰值对应的k组频域数据，利用协方差矩阵的toeplitz 矩阵特性，构建基于频域数据的协方差矩阵。
[0075]
由于方向矩阵a在时间上是独立的，因此fft变换公式可以改写为
[0076][0077]
q＝1,2,
…
,j,m＝1,2,
…
,m
[0078]
其中是a的第m行。由上式可知，每列中的元素幅值相同，相位差固定，且满足
[0079][0080]
其中qk,k＝1,2,
…
,k表示对应于f1,f2,l,fk的列索引。
[0081]
第k个信号的协方差矩阵rk,k＝1,2,
…
,k如下表示：
[0082][0083]
其中是信号的功率，sk,k＝1,2,
…
,k表示信号矩阵s的第k行，r
noise
是噪声协方差矩阵。信号的协方差矩阵近似为toeplitz矩阵(在无噪声情况下)，所以协方差矩阵表达式可以重写为
[0084][0085]
由于协方差矩阵中相邻行元素之间的相位差与式
具有相同的关系，将协方差矩阵的第一行替换为的第qk列，表示为：
[0086][0087]
其中qk是的相位。因此，可以根据和toeplitz特性恢复近似协方差矩阵即
[0088][0089]
步骤s4、基于所述协方差矩阵的特征分解，计算噪声子空间e
nk
，并根据root-music 原理构造多项式求解，获取最接近单位圆的根
[0090]
首先对协方差矩阵进行特征分解如下：
[0091][0092]
得到信号的噪声子空间e
nk
，其中d
sk
表示的最大特征值，d
nk
表示包含其他特征值的对角矩阵；e
sk
代表信号子空间对应的最大特征值的特征向量，噪声子空间e
nk
由对应的其他特征向量组成；
[0093]
然后根据root-music原理构造多项式如下：
[0094][0095]
其中a(z)＝[1,z1,
…
,z
m-1
]为方向向量；获取最接近单位圆的根：
[0096][0097]
即为所需根
[0098]
步骤s5、找到k个最接近于单位圆上的根作为最终解，并基于每个根对应的信号频率解出最终的方向角。
[0099]
最终的精确doa估计可以通过以下方式计算
[0100][0101]
对应于其他频率的其他doa估计也可以由其他协方差矩阵k＝1,2,
…
,k独立地估计，其中数据源来自的第qk列。
[0102]
下面提供一份具体仿真实施例对本发明进行进一步说明。
[0103]
设有2个不相关的跳频信号入射到8个阵元的均匀线阵中。其中这两个不相关信号的doa分别为θ1＝20
°
,θ2＝40
°
，跳频频率为2.40ghz-2.48ghz。采用均方根误差 (rmse)作为性能度量来证明所提出的频域根music算法的有效性，定义如下：
[0104][0105]
其中i是蒙特卡罗实验的数量。表示第i次实验中第k个信号的doa估计。
[0106]
图2为利用本发明的方法对其doa进行估计的散点图，共进行了50次独立估计。仿真条件为：阵元数m＝8；快拍数j＝512；信噪比snr＝0db。很明显可以看出本发明的方法可以有效地估计doa参数，估计值非常接近真实值。
[0107]
图3为利用本发明方法和其他四种算法pm、root-music、capon和esprit算法的估计性能对比图。估计算法性能的标准为均方根误差(rmse)，仿真条件为：阵元数 m＝8；快拍数j＝[200,1000]；信噪比固定为snr＝-5db。可以明显看出，本发明的方法性能优于其他三种算法。
[0108]
图4为利用本发明方法和其他四种算法pm、root-music、capon和esprit算法的估计性能对比图。估计算法性能的标准为均方根误差(rmse)，仿真条件为：阵元数 m＝8；快拍数j＝512；信噪比(snr)从-10db变化到14db。由图可知，本发明的方法性能优于其他三种算法。
[0109]
图5为本发明运用实际测试的硬件，图5a为usrp2943r软件无线电平台，此硬件是national instruments的软件无线电平台，高速链路连接到主机，基于主机的labview 软件用它来控制usrp硬件和发送/接收数据。usrp 2943r具4个可同步正交采样通道，宽达1.2ghz mhz至6ghz的工作频率范围，0db到95db的增益范围，2.5ppm 的频率精度，以及80mhz的最大瞬时实时带宽。图5b是可以发射跳频信号的无人机遥控器，可以发射2.40ghz-2.48ghz的跳频信号，用于实际测试测向精度。
[0110]
图6为本发明使用的软件显示界面示意图，其中展示了单次测向结果，频率计算结果，最大监测频率，最小监测频率等。还用散点图直观显示了角度变化的倾向和历史记录，非常方便观察。
[0111]
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。