一种分段线性调频信号实时生成方法

1.本发明属于信号处理技术领域，具体涉及分段线性调频(piecewise linear frequency modulation,plfm)信号实时生成方法。


背景技术：

2.合成孔径雷达(synthetic aperture radar,sar)作为遥感领域重要的技术手段之一，以全天时、全天候的优势广泛应用于国防、环境监测、自然资源普查和灾害监测等领域。目前，lfm信号是sar系统中最常用的发射波形，但是lfm信号自相关函数的峰值旁瓣比(peak side lobe ratio，pslr)很高，为-13分贝(decibel，db)，会降低sar系统成像性能。为了抑制旁瓣高度，通常采用加权窗函数、滤波或优化等方法，但是这些方法会降低输出信噪比(signal to noise ratio,snr)。非线性调频信号(non-linear frequency modulation,nlfm)自相关函数具有很低的峰值旁瓣比(peak side lobe ratio，pslr)，并且不会降低输出snr，但是目前的技术手段仅能得到nlfm信号的离散值，难以进行工程推广。
3.与nlfm信号相比，plfm信号可以通过调节线性调频成分同样获得很低的pslr，并且这一过程也没有损失输出snr。通过相关实验，plfm信号可以避免snr的1～2db损失，相当于节省天线发射功率25％，对于能量紧缺的雷达系统中，采用plfm信号作为发射波形可以提升系统性能。
4.目前常用的技术手段是通过地面处理计算得到信号的离散值，再把该信号的离散值保存到现场可编程门阵列(field－programmable gate array，fpga)的只读存储器(read-only memory，rom)。在使用时读取rom保存的离散数据，进行数模转换，生成plfm信号。而这种方法需要先采集信号并保存后，才能调用并转换生成plfm信号，受制于rom有限的存储空间，这种存储模式导致plfm信号只能用于科学实验，无法进行工程推广。
5.通过调节线性调频成分可以改变plfm信号自相关函数的峰值旁瓣比，但难以快速地获得低pslr的plfm信号。因此，如何缩短生成plfm信号的时间，实现plfm信号的实时生成是目前亟待解决的问题。


技术实现要素：

6.为解决上述技术问题，本发明提供一种分段线性调频信号实时生成方法。
7.本发明采用的技术方案如下：
8.一种分段线性调频信号实时生成方法，包括如下步骤：
9.步骤101、根据系统参数确定分段线性调频信号参数，包括信号带宽、信号脉宽、分段数量和峰值旁瓣比pslr；
10.步骤102、根据峰值旁瓣比确定分段线性调频信号的功率谱函数；
11.步骤103、对功率谱函数进行积分得到群时延函数；
12.步骤104、根据信号参数和群时延函数计算分段线性调频信号的分段瞬时频率函
数；
13.步骤105、对分段瞬时频率函数积分得到相位函数，并根据相位函数确定分段线性调频信号的时域表达式；
14.步骤106、根据时域表达式生成数字信号，并将数字信号进行数模转换，得到分段线性调频信号。
15.进一步地，所述步骤102包括：
16.以峰值旁瓣比pslr为索引，寻找具有相同峰值旁瓣比pslr的窗函数，然后以该窗函数作为分段线性调频信号的近似功率谱函数；
17.其中，窗函数为汉明窗、升余弦窗或泰勒窗。
18.进一步地，当窗函数为泰勒窗时，将泰勒窗作为分段线性调频信号的近似功率谱函数，其表示为：
[0019][0020]
其中，f为频率轴，表示泰勒系数，b和t表示分段线性信号的带宽和脉宽，η和分别表示旁瓣水平和等高旁瓣数，用以控制峰值旁瓣比pslr，m为计算过程量。
[0021]
进一步地，所述步骤103包括：
[0022]
对窗函数对应的近似功率谱函数进行积分，得到分段线性信号的群时延函数g(f)，表示为：
[0023][0024]
其中，c为常数，表示为：
[0025][0026]
进一步地，当所述窗函数为泰勒窗时，其对应的群时延函数表示为：
[0027][0028]
进一步地，所述步骤104包括：
[0029]
将信号的带宽b均匀分为n段，其中第n(n＝1,2,
…
,n)段的起始点为bn，终点为b
n 1
，每段长度为b/n，且：
[0030][0031]
由此得到带宽向量为：
[0032][0033]
然后根据式(2)中的时延函数计算出频率分段向量对应的脉宽向量
[0034][0035]
其中：
[0036]
tn＝g(bn)
ꢀꢀ
(8)
[0037]
根据脉宽向量和带宽向量计算分段线性调频信号的分段线性频率函数，包括：
[0038]
首先计算每段线性频率函数的调频率kn为：
[0039][0040]
其中，n＝1,2,
…
,n；
[0041]
然后得到分段线性频率函数f(t)为：
[0042][0043]
其中，t为时间轴，ωn＝[tn,t
n 1
),n＝1,2,
…
,n。
[0044]
进一步地，所述步骤105包括：
[0045]
分段线性调频信号的相位函数可以通过下式求解：
[0046][0047]
根据式(10)和(11)，相位函数表示为：
[0048][0049]
其中：
[0050]
c1＝0,c
n 1
＝cn π(b
n 1
bn)(t
n 1-tn),n＝1,2,
…
,n-1
ꢀꢀ
(13)
[0051]
最后，分段线性调频信号的时域表达式为：
[0052][0053]
其中，j为虚数单位，rect(x)表示矩形函数，定义为：
[0054]
[0055]
进一步地，所述步骤106包括：
[0056]
根据系统参数确定信号的信号参数，所述信号参数包括信号带宽、信号脉宽、分段数量n，峰值旁瓣比pslr，系统采样率f；
[0057]
将信号带宽均匀分为n段，得到带宽向量
[0058]
根据峰值旁瓣比pslr确定功率谱密度函数；
[0059]
根据功率谱密度函数确定群时延函数；
[0060]
根据群时延函数计算脉宽向量
[0061]
根据信号脉宽和系统采样率计算时间轴向量，为：
[0062][0063]
其中，其中，表示向下取整，ti＝-t/2 (i-1)t/i,i＝1,2
…
,i；
[0064]
根据分段瞬时频率函数计算分段瞬时频率向量，为：
[0065][0066]
其中：
[0067][0068]
对分段瞬时频率向量进行积分，得到相位向量；
[0069]
根据相位向量确定信号时域离散向量；
[0070]
最后将信号时域离散向量进行数模转换，得到分段线性调频信号。
[0071]
有益效果：
[0072]
传统方法在生成信号时需要插值，使得计算复杂度上升。通过本发明方法，可以避免插值过程，并且只涉及到正余弦、线性函数等基本函数的计算，因此能够降低计算复杂度，有利于信号的实时生成。本发明大大降低信号设计的复杂度，是一种适合fpga实时生成的分段线性调频信号，进一步的，相对于lfm信号能节省发射功率，降低snr损失，从而提高系统性能。
附图说明
[0073]
图1为本发明实施例的一种分段线性调频信号实时生成方法流程示意图；
[0074]
图2为本发明实施例的一种分段线性调频信号实时生成装置结构示意图；
[0075]
图3为本发明实施例的一种基于fpga硬件实现的装置结构示意图；
[0076]
图4为本发明实施例的一种基于sar系统实现的生成plfm信号的时序结构示意图；
[0077]
图5为本发明实施例的plfm信号的时域信号实部图像；
[0078]
图6为本发明实施例进行脉冲压缩后的时域信号图像。
具体实施方式
[0079]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0080]
本发明提供的一种分段线性调频信号实时生成方法，根据需求初始化信号参数；根据设定的功率谱密度函数确定群时延函数；根据分段数量和带宽确定带宽向量；根据群时延函数和带宽向量确定脉宽宽向量；根据带宽向量和脉宽向量确定瞬时频率函数；根据瞬时频率函数计算瞬时频率向量；对瞬时频率向量积分得到信号相位；根据信号相位确定时域信号的离散值。
[0081]
如图1，图2所示，本发明的一种分段线性调频信号实时生成方法具体包括如下步骤：
[0082]
步骤101、根据系统参数确定分段线性调频信号的信号参数，包括信号带宽、信号脉宽、分段数量和峰值旁瓣比pslr。
[0083]
步骤102、根据峰值旁瓣比确定分段线性调频信号的功率谱函数，具体地，以峰值旁瓣比pslr为索引，寻找具有相同pslr的窗函数，然后以该窗函数作为分段线性调频信号的近似功率谱函数。
[0084]
其中，窗函数可以包括：汉明窗(hamming)、升余弦窗(raised cosine)及泰勒窗(taylor)。以泰勒窗为例，将泰勒窗作为分段线性调频信号的近似功率谱函数，其表示为：
[0085][0086]
其中，f为频率轴，表示泰勒系数，b和t表示分段线性信号的带宽和脉宽，η和分别表示旁瓣水平和等高旁瓣数，用以控制峰值旁瓣比pslr，m为计算过程量。
[0087]
步骤103、对窗函数对应的近似功率谱函数进行积分，可以得到分段线性信号的群时延函数g(f)，表示为：
[0088][0089]
其中c为常数，表示为：
[0090][0091]
以泰勒窗为例，其对应的群时延函数可以表示为：
[0092][0093]
步骤104、根据信号参数和群时延函数计算分段线性调频信号的分段瞬时频率函数，包括：
[0094]
首先将信号的带宽b均匀分为n段，其中第n(n＝1,2,
…
,n)段的起始点为bn，终点
为b
n 1
，每段长度为b/n，且：
[0095][0096]
由此得到带宽向量为：
[0097][0098]
然后根据式(2)中的时延函数计算出频率分段向量对应的脉宽向量
[0099][0100]
其中：
[0101]
tn＝g(bn)
ꢀꢀ
(8)
[0102]
进一步，根据脉宽向量和带宽向量计算分段线性调频信号的分段线性频率函数。具体为，首先计算每段线性频率函数的调频率kn为：
[0103][0104]
然后可以得到分段线性频率函数f(t)为：
[0105][0106]
其中，t为时间轴。
[0107]
步骤105、对分段瞬时频率函数积分得到相位函数，并根据相位函数确定分段线性调频信号的时域表达式，包括：
[0108]
分段线性调频信号的相位函数可以通过下式求解：
[0109][0110]
进一步，根据式(10)和(11)，相位函数表示为：
[0111][0112]
其中：
[0113]
[0114]
最后，分段线性调频信号的时域表达式为：
[0115][0116]
其中，j为虚数单位，rect(x)表示矩形函数，定义为：
[0117][0118]
步骤106、根据时域表达式生成数字信号，并将数字信号进行数模转换，得到分段线性调频信号，包括：
[0119]
根据系统参数确定信号的信号参数：包括信号带宽、信号脉宽、分段数量n，峰值旁瓣比pslr，采样率f；将信号带宽均匀分为n段，得到带宽向量
[0120]
根据峰值旁瓣比pslr确定功率谱密度函数；
[0121]
根据功率谱密度函数确定群时延函数；
[0122]
根据群时延函数计算脉宽向量
[0123]
根据信号脉宽和系统采样率计算时间轴向量，为：
[0124][0125]
其中，其中，表示向下取整，ti＝-t/2 (i-1)t/i,i＝1,2
…
,i；
[0126]
根据分段瞬时频率函数计算分段瞬时频率向量，为：
[0127][0128]
其中：
[0129][0130]
对分段瞬时频率向量进行积分，得到相位向量；
[0131]
根据相位向量确定信号时域离散向量；
[0132]
最后将信号时域离散向量进行数模转换，得到分段线性调频信号。
[0133]
下面结合实施例对本发明再作进一步详细的描述。
[0134]
如图3所示，根据输入指令，初始化信号带宽、脉宽等参数。接着利用群时延函数计算每段线性函数。进一步，将线性函数合成瞬时频率函数。然后通过积分得到信号的相位值，并确定信号时域离散值。最后，进行数模转换得到信号输出。
[0135]
本发明实施例提供一种基于sar系统实现的生成plfm信号的时序结构，如图4所示，包括：
[0136]
在某时刻接收到上位机发射的输入指令后，在下一个脉冲重复时间pri开始前先
初始化信号参数；为了不冲掉存储在fpga的ram中的下pri发射的雷达波形数据，以下个pri信号发射窗口结束为开始计算时刻，开始执行计算过程；计算信号，在计算完成对应的plfm信号后，等待输出在下下个pri的信号发射窗口对应模拟plfm信号输出。
[0137]
选择初始化参数为：信号带宽300mhz,信号脉宽20us，分段数量为100。选定-40db泰勒窗函数作为初始化的功率谱密度函数，且系统的采样率为360mhz。生成的plfm信号的时域信号实部如图5所示，进行脉冲压缩结果如图6所示。plfm信号与lfm信号的性能指标结果如下表所示：
[0138][0139]
从上面的描述中可以看出，采用本发明提供的方法，能大大降低信号设计的复杂度，是一种适合fpga实时生成的plfm信号，进一步地，相对于lfm信号能够节省发射功率，降低信噪比损失，从而提高系统性能。
[0140]
以上所述，仅为本发明的实施例而已，并非用于限定本发明的保护范围。本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。