一种基于Bernoulli混沌的四相编码OFDM雷达信号设计方法

霍 凯赵晶晶

（国防科学技术大学电子科学与工程学院 长沙 410073）

一种基于Bernoulli混沌的四相编码OFDM雷达信号设计方法

霍 凯*赵晶晶

（国防科学技术大学电子科学与工程学院 长沙 410073）

相位编码正交频分复用（OFDM）雷达具有许多优良的性能，近年来受到雷达界的广泛关注。但在实际应用中，相位编码OFDM雷达信号存在优良编码数量不多、长度受限、峰均比过高等问题。该文针对这些问题，提出了一种基于Bernoulli混沌的四相编码OFDM雷达信号设计方法，可以构建任意数量和长度的编码，具有更大的设计自由度，并且通过初相加权，得到包络峰均比小于2的混沌四相编码OFDM信号，该信号在高分辨、多普勒等方面具有较多优良特性。

正交频分复用雷达；相位编码；混沌；包络峰均比；模糊函数

引用格式：霍凯， 赵晶晶.一种基于Bernoulli混沌的四相编码OFDM雷达信号设计方法［J］.雷达学报， 2016， 5（4）: 361-372.DOI: 10.12000/JR16050.

Reference format: Huo Kai and Zhao Jingjing.A design method of four-phase-coded OFDM radar signal based on Bernoulli chaos［J］.Journal of Radars， 2016， 5（4）: 361-372.DOI: 10.12000/JR16050.

1 引言

正交频分复用［1］（Orthogonal Frequency-Division Multiplexing， OFDM）技术首先成熟应用于通信领域，Jankiraman等人［2，3］于1998年将OFDM技术引入到雷达系统中，他们设计了名叫“PANDORA”的雷达，采用了OFDM-LFM信号体制。近几年来，OFDM新体制雷达引起了包括美、德等多国在内的世界范围内研究者的兴趣，在国际期刊和会议上涌现出不少高水平的相关论文，从不同角度和背景研究了相关的技术和理论，推动着OFDM雷达向前发展。

OFDM雷达同时发射多个相互正交的子载频来获得大的时宽带宽积，具有距离高分辨能力和多普勒分辨能力。与传统的雷达相比，OFDM雷达具有波形设计灵活、易于数字化、强抗干扰能力和环境适应能力等优点［4-6］。但它在应用上还存在一些技术难点，例如发射波形的包络峰均比（Peak to Mean Envelope Power Ratio， PMEPR）较高［7-10］、多普勒频偏导致的邻道干扰［11，12］等。

相位编码［13，14］是OFDM雷达研究最为广泛的一种信号形式，因为它具有低截获、设计灵活等特点，并且为降低PMEPR提供了较好的条件。然而相位编码OFDM也有其自身的问题，其中最主要的是如何设计具有良好性质的编码。现有已知的具有良好自相关性质的编码数目有限，例如理想的巴克码、Huffman码等自相关旁瓣较低，但是其长度受限且序列固定，容易被截获，安全性较差；互补编码集［13］自相关和互相关性能都很好，但其数量有限，寻码过程复杂。因此，设计出一种PMEPR满足发射需求、编码数目和长度灵活多样、自相关性质优良、生成方案简单的信号就成为制约相位编码OFDM雷达走向实际应用的瓶颈问题之一。通过定义以PMEPR和模糊函数性能为指标的优化函数，并通过优化算法求解编码集，Sebt等人［15］，Alimosaymer等人［16］，Guo等人［17］提出了相应的相位编码OFDM信号优化算法。

本文针对相位编码OFDM雷达的上述问题，提出了一种基于Bernoulli混沌的相位编码OFDM信号设计方法，综合考虑编码数目和长度灵活度、PMEPR、自相关性质等因素，形成一套简单易行的编码生成方案，得到的编码方式灵活，长度任意且具有初值敏感性，具有很好的低截获和抗干扰性能，并且PMEPR小于2，自相关性质优良。目前，已有较多文献研究将混沌理论应用到单载频相位编码信号中［18］，但确切研究混沌应用于相位编码OFDM雷达的文献较少，中国科学院的丁凯［19］在博士论文中研究了混沌码OFDM-SAR系统，在OFDMSAR系统成像、抗干扰性能等方面开展了工作，西北工业大学的Huang等人［20］提出了一种混沌相位编码多载频步进频脉冲串信号，并研究了这种信号的特性。但在混沌相位编码波形设计、编码生成方案等关键问题上未开展深入研究。

本文首先简单介绍了相位编码OFDM雷达信号及其特点，然后提出基于Bernoulli混沌的四相编码OFDM雷达信号生成方案，重点分析了信号的各种特性，最后对这种新的信号体制的应用前景进行了展望。

2 相位编码OFDM雷达信号及其特点

OFDM雷达实质上是一种采用了多载频发射信号体制的雷达，它同时发射多个相互正交的子载频，每个子载频上可采用各种调制方式获得大的时宽带宽积。

设雷达发射信号为：

其中，f0为发射中心载频，M是脉冲个数，Tr是脉冲重复周期，u（t）为OFDM多载频信号的复包络：

其中，N是子载频个数，wn是第n个子载频的频率加权系数，rect（·）为单位矩形窗函数，Tp是脉冲宽度，Δf是相邻子载频之间的频率间隔，满足OFDM条件：

此处，第n个子载频与第m个子载频之间正交是指：

若Tr=Tp，则OFDM雷达为连续波体制。

OFDM雷达的子载频上可以进行任何形式的调制，最常见的有线性调频调制（Linear Frequency Modulation， LFM）和相位编码（Phase Coded， PC）调制，分别称为LFM-OFDM和PC-OFDM，其与基本OFDM信号的区别如图1所示。图1（a）是基本的OFDM信号，每个子载频是单载频，子载频之间的频率间隔为Δf；图1（b）是LFM-OFDM信号，每个子载频采用线性调频，带宽均为Δf；图1（c）是PC-OFDM信号，每个子载频上采用相位编码调制，且子载频频率间隔Δf和码元宽度tb满足Δf= 1/tb。这种PC-OFDM信号脉内、脉间子载频上编码可以相同也可以不同，设计非常灵活，而且还具备一些独特优点，如抗干扰能力强、易于数字化、通信雷达一体化等，因此这种信号是OFDM雷达研究的热点。

PC-OFDM雷达信号的复包络可表示为［7］：

其中，P为发射脉冲数，N为子载频数，K为编码长度，wn，m为第m个脉冲内第n个子载频的加权系数，包括幅度及初相加权，an，k，m为第m个脉冲内第n个子载频的第k个相位编码。tb为码元宽度，Δf= 1/tb为子载频间隔，则脉冲宽度Tp=Ktb，而带宽B=NΔf=N/tb，因此PC-OFDM信号的时间-带宽积表示为BT=KN，为基本OFDM信号时间-带宽积K倍，且与传统的线性调频信号（Linear Frequency Modulation， LFM）信号时间-带宽积受限于调频率不同，PC-OFDM信号时间-带宽积可以通过改变子载频数N和编码长度K灵活设计。

图1 3种OFDM雷达信号时频结构示意图Fig.1 Time-frequency structure of 3 kinds of OFDM radar signal

PC-OFDM雷达信号在各类OFDM雷达信号中具有独特的优势：首先，它继承了相位编码信号的低截获的特点，信号被侦收复制的概率小，安全性能高；其次，它的设计最灵活，由于编码种类繁多、长度各异，因此其各种参数设计自由度大；再次，PC-OFDM信号体制支持通信功能的开发，有利于实现雷达、通信一体化。

PC-OFDM雷达信号也存在一些不足：首先是发射波形的包络峰均比（Peak-to-Mean Envelope Power Ratio， PMEPR）问题，由于A/D采样、放大器等硬件有一定的动态范围，限制了发射信号的PMEPR不能过高，文献［13］表明采用合适的相位编码可以使OFDM雷达发射波形的PMEPR降低到实用要求的2以内，因此如何找到满足条件的相位编码方式就成为一个必须面对的问题；其次现有已知的具有良好自相关性质的编码数目有限，无法满足OFDM雷达实际应用时对编码数目和长度的要求，因此如何设计编码数目和长度灵活多样、具有良好自相关性质的编码是另外一个重要的问题；另外，PC-OFDM雷达信号的处理比较复杂，需要针对不同应用寻求不同的处理方法。

3 基于Bernoulli混沌的四相编码OFDM雷达信号

混沌是指可通过确定性系统产生的，貌似随机的不规则运动［21］。时域混沌信号为对初始值敏感的非周期随机信号，变换到频域为类噪声的宽带连续功率谱，具有很好的信号隐身和保密性能。相比于传统的伪随机信号，混沌信号具有数量众多、可选择性强、长度不受限、且不存在周期性的优点；相比于随机噪声信号，混沌信号可控、易于产生、存储和恢复。因此，选择混沌序列作为相位编码的基础，可以保证编码信号的随机性、安全性，而且编码数据和长度要求可以满足。

用混沌相位编码调制OFDM信号，在多载频技术的基础上结合了相位编码技术，它的编码类型和长度将更加灵活，除了宽带所带来的高的距离分辨能力，还具有较低的自相关函数旁瓣，高的频谱利用率，图钉型模糊函数。根据混沌信号的初值敏感性可生成大量的随机噪声信号，同时又可根据其初始条件和映射规则确定再生同样的随机噪声信号。因此在对雷达回波信号进行压缩处理时，可根据混沌初始条件重生混沌相位编码OFDM信号，可省去体积庞大、电路复杂的延时线部件。由于混沌相位编码OFDM信号的每一个发射脉冲信号都根据不同的初始条件和映射规则产生混沌序列，所以不同脉冲信号具有正交性，因此，将该信号用于成像雷达，可获得低的截获概率，同时在抗干扰、抗衰落、抑制多径与杂波干扰等方面有优越的性能。

（1） f 周期点周期无上界，

这意味着，存在一个不可数集合，此集合中只含有混沌轨道，任意两个轨道之间既不趋近也不趋远，且任意轨道不趋于周期轨道。这样的不可数区间称为混沌吸引域。混沌一般分为两种：一是连续混沌，通过微分或积分获得；一是离散混沌，通过差分迭代获得。考虑到混沌产生、处理的简便，以及混沌相位编码的性质，本文选取Bernoulli混沌作为产生相位编码的基础，它是离散混沌序列的一种，可由下式产生：

其中，xk为过程变量，a， b为混沌映射参数，ε为一接近于0的极小正值，混沌吸引域为［-0.5， 0.5）。b值越接近2，混沌序列值在混沌吸引域［-0.5， 0.5）的遍历性越好。

通过设置x0和a， b后，式（8）可以产生Bernoulli混沌序列。然后对Bernoulli混沌序列进行量化编码，得到四相编码序列，量化规则为：

其中，xk为截取的Bernoulli混沌序列，Np为相位编码数，这里取Np=4。实际上按照本文的设计方法，Np取不同值可以获得二相及多相相位编码。由于相位编码数Np越大，信号安全性越好，自相关旁瓣越低，但信号的复杂度越大，综合考虑安全性能和实现复杂度，优选四相编码，即Np=4；然后，用得到的一个四相编码序列对第m个OFDM雷达脉冲信号各子载频进行相位调制，若每个子载频上编码都相同，则an，k，m=ak，如图1（c）中第1个脉冲所示；若每个子载频上编码不同，则重复生成四相编码序列的过程N次，得到an，k，并令an，k，m=an，k，如图1（c）中第2个脉冲所示。考虑到实现简单，并且后续需要进行初相加权以控制PMEPR，本文只研究每个子载频上用的混沌相位编码序列相同，即an，k，m=ak。最后对每个子载频进行初相加权等手段控制PMEPR，即可获得混沌相位编码OFDM雷达信号。具体步骤如下图2所示。

图2 基于Bernoulli混沌的四相编码OFDM信号生成过程Fig.2 Generation procedure of 4-phase coded OFDM signal based on Bernoulli chaos

步骤1 设定混沌初值x0以及混沌映射参数b，根据1维离散Bernoulli迭代映射方程式（8）产生混沌序列x0， x1， ··， xK-1， ··，可取， a=0.5， 1.5＜b＜2，这里b的取值应接近2，但不能等于2，最好取1.9＜b＜2，这样混沌的遍历性较好；

步骤2 剔除前C个序列点，C的取值范围是C ≥ 1，取K个Bernoulli混沌序列点进行量化编码，量化规则如式（9），这里剔掉前C个序列点是为了使选取的混沌序列性质更好，为了获得随机性更好的混沌序列，还可以采用其它选取方式，如采用固定采样率抽取等，不会影响选取序列的性质；

步骤3 对N个子载频进行wn加权。令wn=，其中为幅度加权，θn为初相加权，若取，则表示只进行初相加权，若取θn=0，则表示只进行幅度加权；令an，k，m=ak，得到各行均相等的相位编码矩阵；用该矩阵对加权后的各子载频进行相位编码调制，得到PC-OFDM信号第m个脉冲上的复包络u（t）。对于脉冲串，可重复上述操作获得每个脉冲互不相同的脉冲串，也可以重复发射这一个脉冲，获得相同脉冲的脉冲串。最终得到的混沌四相编码OFDM雷达脉冲的包络可表示为：

上述过程中，产生基于Bernoulli混沌的多相编码序列并用于PC-OFDM雷达信号的调制是本文方法的关键，相对于常用的编码序列和随机编码序列，这种方法下产生的编码序列有以下主要优势，第一，基于Bernoulli混沌的四相编码序列长度任意，例如，相关性能良好的二相Barker码，比K=13长的二相Barker码全部未知，而多相Barker码序列长度可以更长，但寻码过程复杂，而基于Bernoulli混沌的四相编码序列长度基本不受限制，一般而言，随着编码序列长度增加，混沌特性更加明显，相关性能更好；第二，基于Bernoulli混沌的四相编码序列值更加灵活且易于保存、复现，与Barker码、Frank码、P3、P4码序列值相对固定不同，由于Bernoulli混沌所具有的类随机性，其序列值灵活多变近似于随机编码，但是与随机编码存储较难、基本无法复现的缺点相比，基于Bernoulli混沌的四相编码序列，仅混沌初值是随机选取的，因此只须保存混沌初值，经过相同的混沌映射及量化编码就可以复现性质良好的编码序列。当然，必须看到基于Bernoulli混沌的四相编码序列的相关性能很大程度上受到混沌初值选择的影响，不同混沌初值所产生的混沌序列量化编码序列在相关性能上的差异是很大的，因此需要进一步研究混沌初值选择算法，来获得满足应用需要的编码序列，一种可行的思路是预先设置编码序列峰值旁瓣比（Peak-to-Side-Lobe Ratio， PSLR）的阈值，再在混沌吸引域中随机产生混度初值，并按式（8）、式（9）得到编码序列，计算PSLR，若满足预先设置的阈值条件，则输出序列和对应的混沌初值，否则重复上述步骤。这样的方法保证了编码序列的随机性，但是同时增加了寻码时间，尽管产生一个混沌编码序列的过程算法复杂度很低，但是当阈值条件设置较高时，迭代过程平均耗时势必增加。因此，需要根据实际需要，权衡寻码时间和编码序列PSLR，进行设置。其他混沌初值选择算法是进一步优化该种编码序列的方向之一。

另外，对于PC-OFDM雷达信号而言，子载频数N和编码长度K是其较为重要的两个参数。在码元宽度tb保持不变（相应地，子载频间隔Δf也保持不变）的情况下，子载频数N越大，信号带宽越大，距离分辨率越高，且包络的起伏越小；编码长度K与信号的PMEPR无关，但编码长度越长，信号脉宽越长，速度分辨率越高，且信号的自相关函数旁瓣相对越低。基于Bernoulli混沌的四相编码，由于其编码长度K没有限制，可以自由选择，因此用于PC-OFDM雷达信号的调制可以增加信号设计的灵活性，优化其时间-带宽积，使得信号同时具备高的测距、测速精度，降低脉冲压缩旁瓣，且信号抗干扰、抗截获性能有所提高。

4 信号特性分析

下面从信号PMEPR、模糊函数两方面对基于Bernoulli混沌的四相编码OFDM雷达信号特性进行分析，得出信号综合优化方案。

4.1 发射信号PMEPR

信号u（t）的PMEPR可定义为［7］：

其中，T为一个信号周期。

对于式（10）定义的混沌四相编码OFDM雷达脉冲的包络，由，且采用均一化载频加权系数，即有，包络功率均值为：

因此，降低信号峰值旁瓣比的问题变为降低信号峰值功率，且信号各码元包络是相同的，因此降低整个混沌四相编码OFDM信号峰值的问题可以简化为降低一个码元内信号包络的峰值，且仅与子载频数N与子载频加权系数的选择有关。对于相位编码OFDM信号，常用的降低PMEPR的载频加权方式有3种：一种是幅度加权，采用各种窗函数调制各子载频幅度，即频域加窗处理；二是相位加权，即各子载频加上不同的初相；三是幅度加权和相位加权同时使用。针对混沌四相编码OFDM雷达信号，我们研究了不同途径的效果。

只对各子载频幅度进行加权，首先令θn=0，所用窗函数如表1所示，并对表1中的窗函数进行归一化处理，使得。

表1 幅度加权的窗函数Tab.1 Window functions used for amplitude weighting

只考虑单码元信号包络，图3所示为不同窗函数幅度加权条件下，32个子载频合成OFDM单码元信号包络，过采样率为2， Kaiser窗形状参数β=75。由图可得，当各子载频采用矩形窗函数加权时，包络峰值达到，此时PMEPR为最大值N；对于广义余弦窗Hanning窗而言，包络主瓣峰值及旁瓣略有降低，主瓣展宽并不明显，包络形状没有明显改变，PMEPR降低不明显；对于Kaiser窗而言，由于选取了较大的形状参数值β，主瓣明显展宽，峰值降低明显，能量分布更加均匀。

进一步地，设置子载频数N为不同数值，分别计算信号PMEPR，可以得到不同窗函数的PMEPR随N变化的曲线如图4所示。从图中看出，PMEPR随子载频个数N的增大近乎直线增加，尽管经过Kaiser窗加权，PMEPR大大降低，但仍未达到技术指标要求，说明单一采用幅度加权不能控制信号的PMEPR。

（2） 初相加权

图3 不同窗函数加权单码元OFDM信号包络Fig.3 Single bit of OFDM signal envelope for different window functions

图4 不同窗函数加权OFDM信号PMEPR随子载频个数N的变化Fig.4 PMEPR varies with N for different window functions

表2 初相加权的初始相位Tab.2 Initial phase terms used for phase weighting

32个子载频合成单码元OFDM信号包络见图5。由图可得，初相加权完全改变了信号包络形状，且大大降低了包络峰值，且Newman初相加权和Schroeder初相加权信号包络形状完全一样，只是平移了一个采样点。对不同子载频个数OFDM单码元信号包络，分别计算PMEPR，可以得到不同初相加权信号PMEPR随N变化的曲线如图6所示。从图中看出，表2中的两种初相单独加权是等效的，均可以将PMEPR控制在2左右，并且随着N增大，PMEPR趋于稳定在1.8左右，若要得到更低的PMEPR，则需要将初相加权与幅度加权相结合。

图5 不同初相加权单码元OFDM信号包络Fig.5 Single bit of OFDM signal envelope for different window functions

图6 不同初相加权OFDM信号PMEPR随子载频个数N的变化Fig.6 PMEPR varies with N for different window functions

（3） 幅度和初相加权

选取|wn|为表1中的不同窗函数，θn采用表2中的不同初相。

图7为32个子载频OFDM单码元信号包络。图7（a）中，对于Newman初相，经过比较，优选Kaiser窗形状参数选为β=1.2，此时Kaiser窗近似于矩形窗，但是旁瓣水平相对于矩形窗较低，组合上Newman初相可进一步降低信号包络峰值水平，能量分布更加平均，而组合Hanning窗，信号包络峰值反而有所升高。对于Schroeder初相，其值与加权幅度之间存在耦合，图7（b）中仍选择β=75的Kaiser窗，可以看到，相对于矩形窗情况，无论是Hanning窗还是Kaiser窗，其包络形状大致不变，而波形进一步平滑，幅值降低。

由图8得到两种具有较好效果的组合加权方式：β=1.2的Kaiser窗组合Newman初相加权，Hanning窗组合Schroeder初相加权。且随着N增加，PMEPR值略有降低，且渐趋于平稳。

信号包络峰均比（PMEPR）仅与子载频个数及加权方式有关。以上3种子载频加权方式中，单独使用幅度加窗，信号包络形状未见明显变化，降低PMEPR的效果有限；单独使用初相加权，可以有效降低信号PMEPR；幅度组合初相加权，当选择合适的窗函数与初相进行加权时，可进一步降低信号PMEPR，其中，两种较为有效的方式为形状参数β=1.2的Kaiser窗函数组合Newman初相加权和Hanning窗组合Schroeder初相加权。这两种加权方式下，随着子载频数的增加，信号PMEPR略有下降，渐趋于平稳，均在1.5以下，达到了雷达发射机的技术指标。

图7 不同初相组合幅度加权单码元OFDM信号包络Fig.7 Single bit of OFDM signal envelope for different combination of amplitude and initial phase weighting

图8 不同窗函数、不同初相加权OFDM信号PMEPR随子载频个数的变化Fig.8 PMEPR varies with N for different combination of amplitude and initial phase weighting

4.2 模糊函数及特性分析

模糊函数是分析、比较雷达信号特性的有力工具，它可以用来衡量波形的分辨力、模糊度、测量精度以及抑制杂波能力。其定义为［23］：

PC-OFDM单脉冲信号的模糊函数表达式：

当k＜0或k＞K-1时，ak=0。因此，信号模糊函数性能与子载频加权及编码序列的自相关函数有关。

实验采用上一节中讨论的，包络峰均比控制较好的两种加权方式对N=32个子载频进行加权，并用同一个K=32长的Bernoulli混沌量化四相编码序列对信号各载频进行调制，在混沌吸引域［-0.5， 0.5）中随机产生Bernoulli混沌初值为0.3147。信号码元宽度tb=0.2 ms，则脉冲宽度Tp=6.4 ms，带宽B= 160 MHz，发射基带频率f0=10 GHz。求取两种信号模糊函数如图9所示。两种加权方式下，信号的模糊函数均呈图钉状，具有单一的峰值，其余能量较为均匀地分布于距离-多普勒平面，较为狭窄的中心峰值意味着该种信号具有距离、多普勒高分辨，因而同时具有高精度测速和宽带成像的潜力。

（1） 距离分辨率

其中Ra（i）为相位编码序列非周期自相关函数。特别地，令η=0，整倍码元延迟处，自相关函数表达式进一步写为：

因此，整码元延迟处的旁瓣水平仅由相位编码编码序列的旁瓣水平决定，而子载频加权、子载频数仅影响码元间的旁瓣水平。

图9 不同载频加权方式信号模糊函数图Fig.9 Ambiguity function of signal for different sub-carrier weighting

图10中给出了子载频数分别为N=32， N=64的相位编码OFDM信号与单频相位编码信号归一化自相关函数值对比。由图可见，在码元的整倍延迟处，相位编码OFDM信号自相关函数旁瓣值与所用相位编码序列的自相关函数旁瓣值是一致的，这与式（17）的结论是一致的。雷达信号的距离分辨率与其带宽成反比，而相位编码OFDM雷达信号通过增加子载频数N，可以增大带宽，从而提高距离分辨率，从图中可得，随着子载频数N的增加，主瓣更窄，且整倍码元之间的旁瓣宽度和高度变低。

图10 相位编码OFDM信号模糊函数自相关函数Fig.10 Auto-correlation function of PC-OFDM signal

为了进一步分析子载频加权方式的影响，图11比较了图9中两种加权信号的自相关函数。为了分析两种加权方式对信号距离分辨的影响，截取延迟为一个码元长度内的曲线，可以明显看出，Kaiser窗组合Newman初相加权，主瓣更窄，码元内部旁瓣更低。综合信号PMEPR、距离分辨考虑，Kaiser窗、Newman初相加权是一种较为理想的子载频加权方式。

最后，考虑编码序列长度对相位编码序列自相关函数旁瓣水平的影响。以峰值旁瓣比（Peak-to-Side Lobe Ratio， PSLR）作为指标衡量信号的相关函数旁瓣水平，以相同的混沌初值产生不同长度的Bernoulli混沌量化编码，产生不同相位编码数的子载频数N=32的相位编码OFDM信号，并求取自相关函数PSLR值。重复上述步骤并求最终的平均值，以消除混沌初值对结果的影响。由图12可知，随着编码长度K的增加，信号自相关函数PSLR不断降低，相关性能更加优良。这是由于Bernoulli量化相位编码序列自身的相关性能随着序列长度的增加趋于优良所决定的。然而，相位编码序列长度也不能无限增加，在脉冲宽度的一定的情况下，码元宽度随着序列长度增加而减小，而在现有硬件环境下，相位跳变的速率是有限的，序列长度过长，势必导致误差的产生，且后续处理较为复杂。

图11 不同载频加权方式信号自相关函数Fig.11 Auto-correlation function of signal for different sub-carrier weighting

图12 相位编码OFDM信号自相关函数峰值旁瓣比（PSLR）随编码长度变化曲线Fig.12 PSLR curve varies with K of auto-correlation function of PC-OFDM

图13为本文方法与现有的经典算法［14，15］在距离分辨率上的比较。其中图13（a）为子载频数N=15，相位编码长度K=128，随机产生混沌初值x0=-0.332746454505278， β=1.2 Kaiser窗、Newman初相加权四相编码OFDM信号，图13（b）为采用15位多相Barker码调制，同样为β=1.2 Kaiser窗、Newman初相加权的IS PC-OFDM信号，图13（c）则为最小二乘准则下同时最优化模糊函数、PMEPR 的PC-OFDM信号。其中，图13（a），图13（b）中信号PMEPR均等于1.7107，而图13（c）中信号PMEPR等于1.8949。经过比较，我们可以看到，本文方法通过增长编码序列的长度使得自相关函数旁瓣水平与其他两种方法相当。信号时间-带宽积更大，但是相应地信号处理复杂度有所提高，且旁瓣水平略高于另两种方法，并依赖于初值选择和编码长度。然而，相对于图13（b）中采用相对固定的序列进行编码，本文方法抗截获、抗干扰性能更好，且参数设置更为灵活；图13（c）中的最小二乘准则下的优化算法原则上可以优化设计任意子载频数和编码长度的PC-OFDM信号，而其迭代次数也随着这两个参数的增大而增大，相对而言，本文方法算法复杂度较低，且信号PMEPR控制得也更好。

图13 本文方法与现有经典算法在距离分辨上的比较Fig.13 Comparison with classic algorithm in range resolution

（2） 多普勒分辨率

令τ=0，得到模糊函数零时延截线，即多普勒模糊函数

图14所示为两种加权信号模糊函数零时延截线，与式（18）中的推导一致，加权方式对零时延截线没有影响，主瓣宽度等于1/Tp，且随着多普勒频移 υ的增大，幅度下降明显。

图14 Kaiser窗，Newman初相加权信号模糊函数零时延截线Fig.14 Zero time delay curve of ambiguity function of signal weighted by Kaiser window and Newman initial phase

图15 Doppler频偏对脉冲压缩结果的影响Fig.15 Effects of Doppler offset to results of pulse compression

由本节的分析，基于Bernoulli混沌，通过量化编码产生大量序列长度任意、相关性能较为优良、易于存储与复现的相位编码序列，使得雷达波形设计更加灵活，且具有保密性、抗干扰性，优势明显；经过窗函数和初相组合加权，解决了相位编码OFDM信号包络起伏较大的问题，PMEPR降至1.5左右，且进一步地，压缩了信号自相关函数码元内部的旁瓣值。相位编码OFDM雷达信号模糊函数近似图钉型，同时具有距离、多普勒高分辨，应用前景广阔；而另一方面，相位编码OFDM雷达信号对于多普勒频偏较为敏感，需要在后续处理时加以补偿。综合信号PMEPR和模糊函数特性及后续处理算法的复杂度，一般优选N＞32， 64≤K≤128，Kaiser窗混合Newman初相加权。

5 结束语

相位编码OFDM雷达符合雷达多功能、智能化、数字化、网络化的发展潮流，具有广泛的应用前景，值得我们展开深入研究。本文介绍了相位编码OFDM雷达信号及其特点，提出了基于Bernoulli混沌的四相编码OFDM雷达信号生成方案，并分析了信号的各种特性，得出了一些初步的结论。相对于目前常用的相位编码OFDM信号波形优化算法，基于Bernoulli混沌的四相编码OFDM雷达信号波形参数设计更加灵活，包络起伏控制良好，易于产生，存储、保密性、抗干扰性优良，但是相对而言自相关旁瓣水平较高，依赖于混度初始值选取，且编码长度要求较长，仍需进一步研究改进。

［1］ 佟学俭， 罗涛.OFDM移动通信技术原理与应用［M］.北京: 人民邮电出版社， 2003.Tong Xuejian and Luo Tao.The Principle and Application of OFDM Mobile Communication Technology［M］.Beijing: The People's Posts and Telecommunications Press， 2003.

［2］Jankiraman M， Wessels B J， and Van Genderen P.Design of a multifrequency FMCW radar［C］.The 28th European Microwave Conference， Amsterdam， The Netherland， 1998: 548-589.

［3］Jankiraman M， Wessels B J， and Van Genderen P.Pandora multifrequency FMCW/SFCW radar［C］.IEEE National Radar Conference， Alexandria， VA， USA， 2000: 750-757.

［4］Alimosaymer M and Mohseni R.Systematic approach in designing wavelet packet modulation- orthogonal frequency division multiplexing radar signal by applying the criterion of least-squares［J］.IET Signal Processing， 2014， 8（5）: 475-482.

［5］Sen S.Adaptive OFDM radar waveform design for improved micro-Doppler estimation［J］.IEEE Sensors Journal， 2014， 14（10）: 3548-3556.

［6］Sen S.PAPR-constrained pareto-optimal waveform design for OFDM-STAP radar［J］.IEEE Transactions on Geoscience and Remote Sensing， 2014， 52（6）: 3658-3669.

［7］Levanon N and Mozeson E.Radar Signals［M］.New York: John Wiley & Sons， Inc.， Hoboken， 2004: 327-372.

［8］ 李自琦， 梅进杰， 胡登鹏， 等.限幅法降低OFDM雷达通信一体化系统PAPR研究［J］.雷达科学与技术， 2014， 12（4）: 406-410.Li Ziqi， Mei Jinjie， Hu Dengpeng， et al..Research on deliberate clipping for PAPR reduction of integrated radar and communication systems based on OFDM signals［J］.Radar Science and Technology， 2014， 12（4）: 406-410.

［9］Huang Tianyao and Zhao Tong.Low PMEPR OFDM radar waveform design using the iterative least squares algorithm［J］.IEEE Signal Processing Letters， 2015， 22（11）: 1975-1979.

［10］Wang Wenqin， So H C， Huang Longting， et al..Low peakto-average ratio OFDM chirp waveform diversity design［C］.IEEE International Conference on Acoustic， Speech and Signal Processing （ICASSP）， Florence， Italy， 2014: 8351-8354.

［11］霍凯， 姜卫东， 黎湘， 等.一种新的OFDM相位编码频率步进雷达信号及其特性［J］.电子与信息学报， 2011， 33（3）: 677-683.Huo Kai， Jiang Weidong， Li Xiang， et al..A new OFDM phase-coded stepped-frequency radar signal and its characteristic［J］.Journal of Electronics & Information Technology， 2011， 33（3）: 677-683.

［12］Huo Kai， Deng Bin， Liu Yongxiang， et al..Hign resolution range profile analysis based on multicarrier phase-coded waveforms of OFDM radar［J］.Chinese Journal of Systems Engineering and Electronics， 2011， 22（3）: 421-427.

［13］Levanon N.Multifrequency complementary phase-coded radar signal［J］.IEE Proceedings-Radar， Sonar and Navigation， 2000， 147（6）: 276-284.

［14］Levanon N.Multicarrier radar signals with low peak-tomean envelope power ratio［J］.IEE Proceedings-Radar，Sonar and Navigation， 2003， 150（2）: 71-77.

［15］Sebt M A， Sheikhi A， and Nayebi M M.Orthogonal frequency-division multiplexing radar signal design with optimised ambiguity function and low peak-to-average power ratio［J］.IET Radar， Sonar & Navigation， 2009， 3（2）: 122-132.

［16］Alimosaymer Mostafa and Mohseni Reze.Systematic approach in designing wavelet packet modulationorthogonal frequency division multiplexing radar signal by applying the criterion of least squares［J］.IET Signal Processing， 2014， 8（5）: 475-482.

［17］Guo T and Qiu R.OFDM waveform design compromising spectral nulling， side-lobe suppression and range resolution［C］.IEEE Radar Conference， Cincinnati， USA，2014: 1424-1429.

［18］Harman S and Qineti Q.The Diversity of Chaotic Waveforms in Use and Characteristics［M］.UK: IET Michael Faraday House， 2009: 33-39.

［19］丁凯.混沌信号合成孔径雷达研究［D］.［博士论文］， 中国科学院电子学研究所， 2009.Ding Kai.Researches on chaotic signal synthetic aperture radar［D］.［Ph.D.dissertation］， Institute of Electronics，Chinese Academy of Sciences， 2009.

［20］Huang Qiongdan， Li Yong， Cheng Wei， et al..A new multicarrier chaotic phase coded stepped-frequency pulse train radar signal and its characteristic analysis［C］.IEEE 10th Conference on Industrial Electronics and Applications，Auckland， New Zealand， 2015: 444-448.

［21］郑维敏.正反馈［M］.北京: 清华大学出版社， 1998.Zheng Weimin.Positive Feedback［M］.Beijing: Tsinghua University Press， 1998.

［22］Li T Y and Yorke J A.Period three implies chaos［J］.The American Mathematical Monthly， 1975， 82（10）: 985-992.

［23］Ruggiano M and Van Genderen P.Radar and communication waveform: wideband ambiguity function and narrowband approximation［C］.IET Conference Publications， Edinburgh，UK， 2007: 1-5.

霍 凯（1983-），男，湖北黄冈人，博士，毕业于国防科学技术大学，现任国防科学技术大学电子科学与工程学院讲师，主要研究方向为雷达波形设计、雷达信号处理。

E-mail: huokai2001@163.com

赵晶晶（1990-），女，籍贯江苏省扬州市，现为国防科学技术大学电子科学与工程学院在读博士研究生，主要研究方向为雷达波形设计。

E-mail: xiaosatianyu@163.com

A Design Method of Four-phase-coded OFDM Radar Signal Based on Bernoulli Chaos

Huo Kai Zhao Jingjing
（College of Electronic Science and Engineering， NUDT， Changsha 410073， China）

Orthogonal Frequency-Division Multiplexing （OFDM） radar is receiving increasing attention in the radar field in recent years and is showing excellent performance.However， for practical applications， there are several problems with phase-coded OFDM radar， such as the existence of few good codes， limited length capability， and a high Peak-to-Mean-Envelope Power Ratio （PMEPR）.To address those problems， in this paper， we propose a design method for a four-phase-coded OFDM radar signal based on Bernoulli chaos， which can construct codes of arbitrary amounts and lengths and demonstrate more agility and flexibility.By adopting original phase weighting， this method can obtain a chaotic four-phase-coded OFDM signal with a PMEPR less than two.This signal has excellent performance with respect to high resolution and Doppler radar application.Key words: Orthogonal Frequency Division Multiplexing （OFDM） radar； Phase-coded； Chaos； Peak-to-Mean Envelope Power Ratio （PMEPR）； Ambiguity function

TN911.72

A

2095-283X（2016）04-0361-12

10.12000/JR16050

2016-03-02；改回日期：2016-07-03；网络出版：2016-07-22

霍凯 huokai2001@163.com

自然科学基金青年科学基金（61501481），中国航天科技集团公司航天科技创新基金

Foundation Items: This Research is Funded in Part by the Natural Science Foundation of China （61501481）， The Foundation of China Aerospace Science and Technology Corporation