基于平滑伪魏格纳分布的二相编码信号参数估计方法技术

本发明专利技术涉及一种基于平滑伪魏格纳分布的二相编码信号参数估计方法，主要解决现有方法不能兼顾低信噪比与小计算量的问题。其实现步骤是：1)产生13位巴克码序列编码的二相编码信号；2)对二相编码信号作FFT变换求出其功率谱；3)对功率谱进行平滑处理，并根据平滑后功率谱的3dB带宽功率谱重心来估计信号载频；4)对二相编码信号进行平滑伪魏格纳分布变换，得到载频处的幅度序列；5)搜索幅度序列最小负尖峰间距得到信号码元宽度，并计算相邻负尖峰所包含的子码个数；6)根据负尖峰位置以及子码个数，得到码元序列的估计。本发明专利技术在低信噪比情况下提高了对二相编码信号的参数的估计精度，减小了计算量，可用于对目标的准确侦察。

【技术实现步骤摘要】

本专利技术属于数据处理
，特别涉及一种平滑伪魏格纳分布SPWVD的二相编码BPSK信号参数估计方法，可用于低截获概率雷达中对目标的准确侦察。
技术介绍
为了减小平均发射功率，降低信号被截获的概率，低截获概率雷达通常采用具有大时宽带积的脉冲压缩信号。相位编码信号由于具有较高的距离分辨率和较强的抗侦察、抗干扰能力而成为一种常用的脉冲压缩信号。由于雷达接收到的回波信号往往湮没在大的噪声背景中，因此，在采用数字接收时，如何在低信噪比下快速地估计出相位编码信号的相关参数具有重要的现实意义。针对BPSK信号的参数估计问题,吴江标,万方,郁春来提出了“基于小波变换法的相位编码信号脉内特征提取”见航天电子对抗,2005(3):38-40，文中运用小波变换提取小波脊线的方法来识别相位编码信号，在低信噪比的条件下小波脊线的提取受到噪声的影响很大，影响信号参数估计精度，而且该方法运算量大，不适合工程应用；殷吉昊,汪玲,陈天麒提出了“一种多个二相编码信号载频、码率与编码序列估计方法”见信号处理，2006,22(5):639-543，文中运用时频分析方法-短时傅里叶变换STFT对BPSK信号的码元信息及编码序列进行估计，但是STFT的抗噪性能差，在低信噪比的环境下对信号参数估计精度低；吕明久,周铭,韩俊,邓静提出了“一种低信噪比二相编码信号时频分析方法”见雷达科学与技术,2010(6):549-552.文中提出对信号STFT变换后的结果再进行双阈值突变点检测，该方法虽然实现了在低信噪比的条件下对BPSK信号的参数估计，但是估计精度低，而且在估计过程需要对脉冲进行累积，所需的处理时间长，计算量大。
技术实现思路
本专利技术的目的在于针对上述现有方法的不足，提出一种基于平滑伪魏格纳分布的二相编码信号参数估计方法，以提高在低信噪比条件下对二相编码BPSK信号载频、码宽和编码序列的估计精度，同时减少计算量。为了实现上述目的，本专利技术的技术方案包括如下：(1)产生由13位巴克码序列[11111-1-111-11-11]编码的二相编码信号s(t)；(2)对二相编码信号s(t)进行FFT变换，得到其功率谱Y(m)；(3)对功率谱Y(m)进行平滑处理，得到平滑后的功率谱Ys(k)；(4)根据平滑后功率谱Ys(k)的3dB带宽功率谱重心，估计出信号载频(5)获得幅度序列5a)对二相编码信号s(t)进行平滑伪魏格纳分布SPWVD变换，得到二维幅频矩阵SPWVD(t,f)，5b)对二维幅频矩阵SPWVD(t,f)在频率维进行搜索,得到信号载频处的幅度序列(6)计算二相编码信号s(t)的码元宽度6a)搜索幅度序列得到各个负尖峰位置，计算出相邻负尖峰间隔Δ；6b)取最小相邻负尖峰间隔Δ记为Δmin，则信号码元宽度的估计值为：其中,Ts为信号采样周期；(7)计算由负尖峰分隔开的每段数据所包含的子码个数l：l＝Δ/Δmin；(8)计算二相编码信号的码元序列8a)根据负尖峰位置以及由负尖峰分隔开的每段数据所包含的子码个数l，粗略估计出码元序列8b)对粗略估计出的码元序列按照如下公式器得到码元序列的精确估计其中，||·||是取范数运算，是由估计出的信号载频码元宽度码元序列以及时间t重新构造出的二相编码信号，argmin是取最优解。本专利技术与现有技术相比具有如下优点：1、载频估计精度高本专利技术采用改进的功率谱平滑公式对二相编码BPSK信号的功率谱进行滤波处理，并根据平滑后功率谱的3dB带宽重心来估计信号载频，相比于现有技术对时频矩阵最大值序列取平均值来估计载频的方法，提高了在低信噪比条件下对二相编码BPSK信号的载频估计精度。2、码元宽度的估计精度高本专利技术采用平滑伪魏格纳分布SPWVD变换对二相编码信号的码元宽度进行估计，相比于现有技术采用ZAM变换的方法，SPWVD变换具有更好的时频聚集性，提高了在低信噪比条件下对二相编码信号的码元宽度估计精度。3、计算量小本专利技术采用平滑伪魏格纳分布SPWVD变换将二相编码信号的相位信息映射为信号载频处的幅度信息，相比于现有技术采用短时傅里叶变换的方法，对二相编码信号的参数估计不需要对数据进行多次累积以及多个阈值判断，减小了计算量。下面结合附图和具体实施方式对本专利技术做进一步详细说明。附图说明图1是本专利技术的实现流程图；图2是用本专利技术得到的二相编码信号的载频估计精度随信噪比变化曲线；图3是用本专利技术与现有方法得到的二相编码信号的码元宽度估计精度随信噪比变化曲线对比图。具体实施方式参照图1，本专利技术的实现步骤如下：步骤1:获得二相编码信号s(t)。按如下公式产生由13位巴克码序列[11111-1-111-11-11]编码的二相编码信号s(t)：其中,SNR为信噪比，fc为载频，t为时间，τ为回波延时,TP为码元宽度，n为巴克码序列长度，φ为初始相位，φi为每个子脉冲的相位，对于二相编码信号，φi∈{0,π本文档来自技高网...

【技术保护点】
基于平滑伪魏格纳分布的二相编码信号参数估计方法，包括：(1)产生由13位巴克码序列[1 1 1 1 1 ‑1 ‑1 1 1 ‑1 1 ‑1 1]编码的二相编码信号s(t)；(2)对二相编码信号s(t)进行FFT变换，得到其功率谱Y(m)；(3)对功率谱Y(m)进行平滑处理，得到平滑后的功率谱Ys(k)；(4)根据平滑后功率谱Ys(k)的3dB带宽功率谱重心，估计出信号载频(5)获得幅度序列5a)对二相编码信号s(t)进行平滑伪魏格纳分布SPWVD变换，得到二维幅频矩阵SPWVD(t,f)，5b)对二维幅频矩阵SPWVD(t,f)在频率维进行搜索,得到信号载频处的幅度序列(6)计算二相编码信号s(t)的码元宽度6a)搜索幅度序列得到各个负尖峰位置，计算出相邻负尖峰间隔Δ；6b)取最小相邻负尖峰间隔Δ记为Δmin，则信号码元宽度的估计值为：T^b=Δmin·Ts,]]>其中,Ts为信号采样周期；(7)计算由负尖峰分隔开的每段数据所包含的子码个数l：l＝Δ/Δmin；(8)计算二相编码信号的码元序列8a)根据负尖峰位置以及由负尖峰分隔开的每段数据所包含的子码个数l，粗略估计出码元序列8b)对粗略估计出的码元序列按照如下公式器得到码元序列的精确估计P^m=argmin{||s^(f^0,T^b,P~m,t)-s(t)||},]]>其中，||·||是取范数运算，是由估计出的信号载频码元宽度码元序列以及时间t重新构造出的二相编码信号，arg min是取最优解。...

【技术特征摘要】
1.基于平滑伪魏格纳分布的二相编码信号参数估计方法，包括：(1)产生由13位巴克码序列[11111-1-111-11-11]编码的二相编码信号s(t)；(2)对二相编码信号s(t)进行FFT变换，得到其功率谱Y(m)；(3)对功率谱Y(m)进行平滑处理，得到平滑后的功率谱Ys(k)；(4)根据平滑后功率谱Ys(k)的3dB带宽功率谱重心，估计出信号载频(5)获得幅度序列5a)对二相编码信号s(t)进行平滑伪魏格纳分布SPWVD变换，得到二维幅频矩阵SPWVD(t,f)，5b)对二维幅频矩阵SPWVD(t,f)在频率维进行搜索,得到信号载频处的幅度序列(6)计算二相编...

【专利技术属性】
技术研发人员：张林让，蒙妍，许建华，彭炜杰，赵永志，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61