非合作跳频信号参数的盲压缩感知估计

2016-12-03郭树旭

雷达学报 2016年5期

关键词：压缩比参数估计先验

陈莹钟菲郭树旭*

①(吉林大学电子科学与工程学院长春 130012)

②(长春工程学院长春 130012)

非合作跳频信号参数的盲压缩感知估计

陈莹①钟菲②郭树旭*①

①(吉林大学电子科学与工程学院长春 130012)

②(长春工程学院长春 130012)

针对非合作跳频通信系统采样速率高，先验信息少等问题，论文提出基于盲压缩感知重构理论参数估计算法。利用稀疏编码与正交基变换交替迭代的思想实现信号精确重构，并根据重构结果直接对跳频信号进行参数估计。与传统的压缩感知理论相比，盲压缩感知理论避免了对信号先验信息的需求，有效解决了非合作通信系统中先验信息少的问题。首先，建立信号模型，然后利用正交块对角盲压缩感知算法(Orthonormal Block Diagonal Blind Compressed Sensing, OBD-BCS)实现信号的重构，并估算出跳变频率及跳变周期。通过实验分析，该方法可以在低信噪比环境下恢复信号原始结构及信息，完成参数估计。

盲压缩感知；跳频信号；非合作；参数估计

1 引言

信息安全在军事通信及雷达传输领域至关重要，为解决无线通信领域中存在的窃听，电子干扰及信道混叠等问题，非合作跳频通信被广泛应用于其中。非合作跳频通信系统发射端与接收端的载波频率在既定集合内同步发生跳变，具有良好的抗截获、抗干扰及组网能力[1]。在非合作跳频通信系统中，由于发射机与接收机之间具有非匹配性，使得在接收端接收到的跳频信号往往只包含有限信息，信号结构信息及相关调制参数均未知。而且，系统中同时含有其他干扰信号，如宽频带噪声、多频带噪声等。因此，如何在缺少先验信息，并且存在其他干扰条件下实现对信号的参数估计，已经成为非合作通讯侦查领域一个重要研究课题。

目前，跳频信号参数估计的主流方法是时频分析法，1997年，Barbarossa等人[2]首次提出运用Wigner-Ville分布(Wigner-Ville Distribution, WVD)，在时域分析跳频信号参数特性，这种方法计算简单，但是由于存在交叉干扰，在一定程度上影响了实验结果，降低了估计精度。为解决这一问题，文献[3,4]提出基于平滑伪 Wigner-Ville 分布(Smoothed Pseudo Wigner-Ville Distribution, SPWVD)和基于重排列平滑伪 Wigner-Ville 分布(Rearrangement Smoothed Pseudo Wigner-Ville Distribution, RSPWVD)的跳频信号参数估计方法，这两种方法利用信号分量能量重新分布实现参数估计，由于算法具有时频不变性及能量守恒的特点，因此能够改善频谱重叠现象，增强时频聚合性。为解决时频分析法中信噪比阀值的约束性，文献[5]利用多物种粒子群实现参数估计，这种方法通过搜索过完备字典，得到跳频分量最优匹配原子，再利用粒子群参数求解估计出相应跳变参数。但其对粒子距离测度阀值有严格要求，过大或过小都会影响估计精度。文献[6]引入压缩感知(Compressive Sensingm, CS)参数估计算法，利用少量不相干观测值实现跳频信号高速采样及参数估计，这种算法有效地降低了采样频率，但是在构造观测矩阵的过程中仍然需要部分先验信息，因此并不能为非合作跳频通信系统中信号的解调与破译提供足够的支持。文献[7]引入局部特征尺度分解算法，实现跳变参数盲估计，有效克服了时频不确定性及交叉干扰等影响，但是在实现过程中其他参数估计结果依赖于跳速估计值，当跳速估计错误时，会产生较大的误差。

近几年来，跳频信号的频率范围有逐渐增大的趋势，同时考虑到非合作通信系统中信号先验信息少，并且存在其他干扰信号，本文提出基于盲压缩感知(Blind Compressed Sensing, BCS)理论参数估计法，首先利用OBD-BCS算法实现非合作跳频信号重构，然后通过对重构结果进行频域分析，估计出相应跳变参数。理论分析及仿真实验表明，该方法能够在信噪比低，先验信息未知的条件下准确估计出信号的跳变频率及跳变周期。

2 盲压缩感知理论

压缩感知理论是一种基于稀疏表示[8–10]，不相干观测及最优逼近理论提出的采样及重构方法，目前已被广泛应用于雷达信号的参数估计等领域[11]。对于信号根据其稀疏特性，构造一个N×N的正交变换基则信号可表示为是x的加权系数。另构造一个与正交基不相干的平稳测量矩阵即可得到压缩后观测向量利用l0范数的最优化求解，可以得到信号逼近信号即

CS理论在构造矩阵Ψ时，信号的先验信息是至关重要的，但在实际应用过程中，信号信息往往不能被事先获知，如本文涉及的非合作跳频信号。为解决这一问题，BCS理论被提出[12]，它将CS理论和字典学习(Dictionary Learning, DL)结合起来[13]，一方面满足CS理论可以从少量数据中重构出高维信号的特性。另一方面，与字典学习过程一样，不需要获取信号稀疏基的先验信息。BCS理论省略了信号稀疏表示的过程，在采样和重构的过程中避免了对信号先验信息的需求，这使得BCS理论适用于非合作通信系统中信号的处理。

然后利用稀疏编码与正交基变换的不断交替迭代，直到满足

由于省略了稀疏表示的过程，信号观测值的维度较高，这使得BCS理论的解具有不确定性，为了确保解的唯一性，正交基P需要满足以下限制条件之一[14,15]：

(1) P是某个有限已知基字典Λ的正交子集。

(2) P可以被基字典Λ稀疏表示。

(3) P是一个块正交对角矩阵。

条件(1)可以被理解为CS理论重构过程，条件(2)可以被视作字典学习过程，考虑到跳频通信系统的跳变特性，本文主要是在正交基P满足条件(3)的前提下实现信号的重构，即有：

相对于条件(1)和条件(2)，当P满足条件(3)时，在对P进行变换过程中，可以对每一块单独变换而不需要改变整体序列。此时，中含有K项非零列，为正交观测矩阵，并有

3 信号模型及重构过程

3.1 跳频信号模型

跳频信号(Frequency Hopping, FH)载波频率按照一定变化规律迅速发生跳变，其间隔通常为数据调制带宽，因此，信号对于随机性，结构一致性以及调制带宽有严格要求。其信号模型为：

式中，S表示信号功率，K表示总的跳变次数，Th表示每次跳变持续时间，T0表示起跳时刻，fi表示第i次的跳变频率，所对应时刻为T0+(i+1)Th, n(t)表示系统中所含有的高斯白噪声，窗函数

图1给出了根据定义式(4)仿真得到的FH信号时频曲线图，从图中可以观察出跳变频率与跳变周期的对应关系。

图1 FH信号时频曲线图Fig. 1 Time and frequency relation of frequency hopping signals

3.2 跳频信号重构

跳频信号具有稀疏特性，这满足BCS理论前提条件。对于跳频信号而言，每一次跳变都可以表示为一段持续时间为Th的正弦信号，在全部采样周期内，信号可以看成是K个正弦信号的线性组合。根据跳频信号这一特点，可以将其表示成由于在BCS理论中，正交基满足式(3)，因此可以推断出信号Xi在Pi下分别满足稀疏特性。

BCS理论是根据字典学习过程进一步推导得出的，字典学习的目的是利用基字典D求解：

为了使式(6)达到最小，我们对正交基P进行按列迭代。

具体的重构方法为：

步骤4 固定正交基P，根据B=DS，再利用OMP算法对S进行重构，可以得到稀疏矩阵

其中，Tr[.]表示矩阵的求迹运算，式(8)表示矩阵BTAPS主对角线各元素之和最大，根据奇异值分解理论(Singular Value Decomposition, SVD)，可得：

步骤8 重复步骤4–步骤8，直到满足式(2) 。

我们将信号变换到频域范围内，可以得到如图2所示的信号原始结构和重构结果对比图，从图中可以观察出，原始信号和重构信号基本重合，且重构误差接近10–9，因此BCS理论能够实现非合作跳频信号重构。

图2 FH信号重构频谱图Fig. 2 Reconstruction of frequency hopping signals

4 跳变参数估计

4.1 跳变频率估计

跳频信号的自相关函数及功率谱具有瞬时变化的特点，其全部能量包含在一定跳变周期内的窄带宽中，这类非平稳信号难以在时域环境下分析其特性，因此我们可以利用傅里叶变换，将信号变换到频域，以获得其全部频域信息，进而实现跳变频率估计。具体实现方法为：首先对重构跳频信号进行快速傅里叶变换(Fast Fourier Transform, FFT)，可以得到如图2所示的频谱图；然后对频谱图进行谱峰搜索，找到峰值所对应频率点，即为估计频率。由于信号中存在噪声，为提高估计精度，我们取信号峰值附近相邻M个点的频率平均值作为估计结果，即有

而峰值个数即为跳变次数，即k=K。

4.2 跳变周期估计

假设，实验中跳频信号是同步获得，则基于BCS理论可知：只有当信号长度w等于采样点个数m时，信号才能实现精确重构，式(2)的解才具有唯一性，则信号可表示为其中，为K阶稀疏矩阵，此时满足：

的值达到最大。

因此，以式(11)的计算结果作为判定条件，当Q值达到最大时，利用跳变次数k及采样点个数m求得信号的跳变周期

4.3 参数估计流程

假设实验中所用到的非合作跳频信号含有高斯白噪声，则具体参数估计步骤可表述为：

步骤4 观测峰值个数k。

步骤5 利用观测到的峰值个数k，及实验获得的采样点个数m，根据可以计算出信号的跳变周期

5 仿真结果及性能分析

为验证非合作跳频通信系统中，基于BCS理论的重构及跳变参数估计的有效性，现根据定义式(4)仿真得到一段含有高斯白噪声的跳频信号，信号包含1280个采样点，共发生10次跳变，跳变频率fK={40,15,45,20,5,35,10,25,50,30} kHz，采样频率fs=100 kHz，跳变周期Th=1.28 ms，起跳时刻T0=0时刻。实验所采用的评价指标为：

其中，式(12)为信号的重构误差，式(13)为归一化平均误差(Normalized Mean Square Error, NMSE),为第i次的参数估计值，为第i次的实际参数值，Nobs为总的估计次数。

5.1 信号重构实验

实验中，我们对比了取不同观测值时OBDBCS, IRLS及OMP算法的重构精度，并针对OBDBCS算法取信噪比(Signal Noise Ratio, SNR)为–5～5 dB，压缩比为L=2, L=4及L=8 3种情况进行具体分析，根据式(12)可以得到如图3，图4所示非合作跳频信号重构误差拟合曲线。

图3 非合作跳频信号重构误差Fig. 3 Reconstruction error of non-cooperative frequency hopping signal

图4 不同压缩比的跳频信号重构误差Fig. 4 Reconstruction error under different compression ratio

通过图3中拟合曲线可以观察出，对于传统CS理论而言，由于在非合作通信系统中未知信号先验信息，导致正交基的选取存在一定偏差，因此IRLS及OMP算法的重构误差较大，而BCS理论避免了对信号先验信息的需求，因此对于非合作信号而言，BCS理论重构精度优于传统CS理论。

分析图4中的拟合曲线可知，OBD-BCS算法的重构精度与压缩比及SNR有关。当SNR相同时，信号压缩比降低的同时，采样点数增多，重构精度逐渐增强；压缩比一定时，SNR减小会导致噪声干扰增强，重构误差增大。

5.2 跳变频率估计实验

实验中，SNR取0～5 dB，在不同SNR下，我们针对压缩比为L=2, L=4及L=8 3种情况分别进行多次仿真实验。为了减小估计误差，我们以Auf–Aui≤ 10 dB为判定条件，其中，Auf为峰值幅度，Aui为第i个点的幅值。选取峰值附近M个点的频率平均值作为估计结果，并得到频率估计归一化平均误差曲线如图5所示。

图5 不同压缩比的频率估计Fig. 5 NMSE of frequency estimation under different compression ratio

实验结果表明，SNR降低或压缩比增大，均会导致信号频率估计误差增大。由于我们求取了平均值作为频率估计值，因此在不同SNR及压缩比条件下，信号频率估计误差较为接近。

5.3 跳变周期估计实验

为了验证周期估计方法的有效性，我们比较了不同压缩比及SNR条件下的实验结果，在实验中，当连续多个采样点的幅值达到Au ≥ 70 dB时，则视为一个峰值，即发生一次跳变；如果有不连续的采样点的幅值达到70 dB，则可视为是噪声干扰。图6给出了压缩比为L=2, L=4及L=8, SNR为–8～5 dB时的拟合曲线。

实验数据表明，当L=2时，跳变周期估计误差在SNR=–5 dB时出现转折，当SNR ≥ –5 dB时，误差开始接近为0；同理当L=4及L=8时，周期估计误差分别在SNR=–4 dB及SNR=–2 dB时出现转折。通过分析实验过程及数据，可以发现这是由于跳变周期估计值受峰值个数影响所导致，当SNR较低时，受噪声干扰，信号的峰值较为杂乱，不能精准地确定峰值个数，因此会影响跳变周期的估计结果。随着信噪比增大，噪声干扰减弱，可以精确搜索出峰值个数，因此重构误差接近为0。仿真实验结果证明了盲压缩感知理论用于参数估计的可行性及有效性，适用于非合作通信系统中信号的处理。

图6 不同压缩比的周期估计Fig. 6 Hop period estimation under different compression ratio

6 结论

本文提出利用BCS理论对非合作跳频信号进行参数估计，BCS理论兼具CS理论及DL理论特性，能够有效克服稀疏表示过程中缺少信号先验信息的问题，因此能够被应用于收发双方不具备匹配协议的非合作跳频通信系统中。为了确保非合作跳频信号处理过程中BCS理论解的唯一性，通过正交基更新过程构造正交块对角矩阵，使其满足BCS理论的限制条件，并通过OBD-BCS算法实现非合作环境中信号的精确重构及参数估计。通过理论分析及仿真实验表明，BCS理论能够有效克服系统中存在的干扰，从少量的采样值中重构出原始信号结构，且不需要获得信号的先验信息。实验结果显示，随着非合作跳频信号压缩比减小及SNR的增大，噪声干扰逐渐减弱，重构及估计误差逐渐减小。

综上所述，BCS理论能够在信噪比低、先验信息少的非合作跳频通信系统中，利用部分采样值实现对信号的精确重构及参数的准确估计。这为非合作跳频通信过程中信号的处理提供了一个全新方法。

[1]Wang L, Zhang B, and Zhao Y. Compressive sampling and rapid reconstruction of broadband frequency hopping signals with interference[J]. Circuits, Systems, and Signal Processing, 2015, 34(5): 1535–1547.

[2]Barbarossa S and Scaglione A. Parameter estimation of spread spectrum frequency-hopping signals using timefrequency distributions[C]. First IEEE Signal Processing Working on Signal Processing Advances in Wireless Communications, Paris, 1997: 213–216.

[3]赵俊, 张朝阳, 赖利峰, 等. 一种基于时频分析的跳频信号参数盲估计方法[J]. 电路与系统学报, 2003, 8(3): 46–50. Zhao Jun, Zhang Chao-yang, Lai Li-feng, et al.. Blind parameter estimation of frequency-hopping signals based on time-frequency analysis[J]. Journal of Circuits and Systems, 2003, 8(3): 46–50.

[4]Lv J F, Han Y, and Liang X P. The parameter estimation of non-cooperative frequency hopping signals based on the algorithm of RSPWVD[J]. Applied Mechanics and Materials, 2014, 556/562: 4779–4783.

[5]郭建涛, 王宏远, 余本海. 基于粒子群算法的跳频信号参数估计[J]. 计算机应用研究, 2010, 27(2): 512–514. Guo Jian-tao, Wang Hong-yuan, and Yu Ben-hai. Parameter estimation of frequency hopping signal based on particle swarm optimization[J]. Application Research of Computers, 2010, 27(2): 512–514.

[6]Li B, Li Y, and Zhu Y. Compressive frequency estimation for frequency hopping signal[C]. TENCON 2013–2013 IEEE Region 10 Conference (31194), 2013: 1–4.

[7]吕晨杰, 王斌, 唐涛. 采用局部特征尺度分解的跳频信号参数盲估计算法[J]. 信号处理, 2015, 31(3): 308–313. Lv Chen-jie, Wang Bin, and Tang Tao. Blind parameter estimation of frequency hopping signal using local characteristic-scale decomposition[J]. Journal of Signal Processing, 2015, 31(3): 308–313.

[8]Donoho D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289–1306.

[9]Candès E J, Romberg J, and Tao T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information[J]. IEEE Transactions on Information Theory, 2006, 52(2): 489–509.

[10]Candè E J and Wakin M B. An introduction to compressive sampling[J]. IEEE Magazine Signal Processing, 2008, 25(2): 21–30.

[11]朱莹, 张弓, 张劲东. 基于 DCS 的统计 MIMO 雷达信号模型及参数估计[J]. 雷达学报, 2012, 1(2): 143–148. Zhu Ying, Zhang Gong, and Zhang Jin-dong. Signal model and parameters estimation of statistical MIMO radar based on distributed compressed sensing[J]. Journal of Radars, 2012, 1(2): 143–148.

[12]Gleichman S and Eldar Y C. Blind compressed sensing[J]. IEEE Transactions on Information Theory, 2011, 57(10): 6958–6975.

[13]Kreutz-Delgado K, Murray J F, Rao B D, et al.. Dictionary learning algorithms for sparse representation[J]. Neural Computation, 2003, 15(2): 349–396.

[14]Fang B, Huang G, and Gao J. Sub-nyquist sampling and reconstruction model of LFM signals based on blind compressed sensing in FRFT domain[J]. Circuits, Systems, and Signal Processing, 2015, 34(2): 419–439.

[15]方标, 黄高明, 高俊. LFM 宽带雷达信号的盲压缩感知模型[J].航空学报, 2014, 35(8): 2261–2270. Fang Biao, Huang Gao-ming, and Gao Jun. A blind compressed sensing, model for linear frequency modulated wideband radar signals[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(8): 2261–2270.

Blind Compressed Sensing Parameter Estimation of Non-cooperative Frequency Hopping Signal

Chen Ying①Zhong Fei②Guo Shuxu①

①(College of Electronic Science & Engineering, Jilin University, Changchun 130012, China)

②(Changchun Institute of Technology, Changchun 130012, China)

To overcome the disadvantages of a non-cooperative frequency hopping communication system, such as a high sampling rate and inadequate prior information, parameter estimation based on Blind Compressed Sensing (BCS) is proposed. The signal is precisely reconstructed by the alternating iteration of sparse coding and basis updating, and the hopping frequencies are directly estimated based on the results. Compared with conventional compressive sensing, blind compressed sensing does not require prior information of the frequency hopping signals; hence, it offers an effective solution to the inadequate prior information problem. In the proposed method, the signal is first modeled and then reconstructed by Orthonormal Block Diagonal Blind Compressed Sensing (OBD-BCS), and the hopping frequencies and hop period are finally estimated. The simulation results suggest that the proposed method can reconstruct and estimate the parameters of noncooperative frequency hopping signals with a low signal-to-noise ratio.

Blind Compressed Sensing (BCS); Frequency hopping signal; Non-cooperative; Parameter estimation

TN911

2095-283X(2016)05-0531-07

10.12000/JR15106

陈莹, 钟菲, 郭树旭. 非合作跳频信号参数的盲压缩感知估计[J]. 雷达学报, 2016, 5(5): 531–537.

10.12000/JR15106.

Reference format: Chen Ying, Zhong Fei, and Guo Shuxu. Blind compressed sensing parameter estimation of non-cooperative frequency hopping signal[J]. Journal of Radars, 2016, 5(5): 531–537. DOI: 10.12000/JR15106.