Costas序列在雷达信号设计中的应用研究

1、Costas序列在雷达信号设计中的应用研究Costas序列在雷达信号设计中的应用研究类别：通信网络0 引言 雷达的距离分辨率和速度(径向)分辨率取决于雷达选择的信号形式，雷达信号在频域上占据的频带越宽，则其距离分辨率越好；雷达信号在时域上持续宽度越大，则其速度分辨率越佳。从提高雷达分辨率的角度设计雷达信号，就要求信号模糊函数的主峰高而尖锐，副峰低而平坦。常用的雷达信号，例如：LFM(线性调频)信号对速度不敏感，不能用来测速，NLFM(非线性调频信号)自相关函数的旁瓣电平有所改善，但在模糊函数高多普勒频率截面上仍存在较大的距离旁瓣，大目标或杂波的旁瓣将掩盖旁瓣附近小目标的主瓣，在多目标环境中，多

2、个目标响应旁瓣的合成，甚至可能掩盖较强目标响应的主瓣。又由于线性调频信号存在着多普勒频移与距离的耦合，因此当目标回波的多普勒频移较大时将产生较大的测距误差。 本文介绍一种FHSS(跳频扩频)信号，雷达信号如果按Costas序列的构成方法进行跳频，将得到理想的模糊函数性能。 1 Costas序列的概念 设P为n阶置换矩阵，若序列P的(离散)自相关函数R(，d)副瓣的最大值不大于1，则称置换矩阵P为n阶Costas矩阵，序列P称为Costas序列。 为线性调频信号与Costas信号及其离散自相关函数的图形。其中：图(a)为正斜率的线性调频信号，图(e)为其离散自相关函数；图(b)为负斜率的线性调频

3、信号，图(f)为其离散自相关函数；图(c)为6阶Costas信号，6阶Costas序列可用有限域理论生成，其本原元为=3，放置函数为：y(k)=k(mod7)，k=1，2，6，图(g)为其离散自相关函数；图(d)为7阶Costas信号，图(h)为其离散自相关函数。 a) 线性调频信号的离散自相关函数具有较高的副瓣，而Costas信号离散自相关函数副瓣的最大值为1。 b) 在图(e)(或图(f)中，副瓣的总和为：2(5+4+3+2+1)=30，在图(g)中，副瓣的总和(即1的个数)为：6(6-1)=30。还可以发现，在图(e)(或图(f)中每一行(或列)副瓣的数值与图(g)中对应的行(或列)副瓣

4、1的和相等。由模糊体积不变性原理，可以知道任何调制都不能改变模糊曲面下的总容积。该容积只取决于信号的能量，而与信号的形式无关，但可选择适当的信号形式，也就是选择不同形状的模糊曲面，使其与特定的目标环境图相匹配，以实现在所需要分辨目标的区域，使模糊图的体积分布小些，从而提高分辨力。例如：可以采用非线性调频信号来改善线性调频信号的副瓣性能和提高多普勒灵敏度，从这种意义上说，可以认为Costas信号是一种特殊形式的非线性调频信号。其模糊函数图形具有像相位编码信号一样的图钉状特性。在时域宽度一定的情况下，信号按照Costas序列构成方法进行FHSS，利用Costas序列特殊的序列结构，以及信号频谱的扩

5、展换取了信号理想的模糊函数性能。Costas序列特殊的序列结构可使模糊函数的副峰低而平坦，信号频谱的扩展可使模糊函数的主峰高而尖锐。 判断n阶置换序列P是否为Costas序列，需要计算序列P的自相关函数，当n较大时，计算工作量是很大的。利用有限域理论，能够快速构造Costas序列。 2 有限域简介 具有有限个元素的域称为有限域，有限域又称为伽罗瓦(Galois)域，将q阶有限域记作GF(q)。 任何两个元素个数相同的有限域是同构的。两个同构的域，如果不管它们的实际背景而只考虑它们的代数性质，可以将它们等同起来看做一个域。 GF(q)(q25时，其Costas序列的数量将显著减少。特别地，当n趋

6、于无穷大时，Costas序列数趋于0。 d) 我们还不能回答，当n满足什么条件时，不存在Costas序列，但当n=32、33或43时，不存在Costas序列。 4 信号的模糊函数 设u(t)为复包络，令 x(，fd)称为信号u(t)的模糊函数。 x(，fd)沿着fd=0的截面为： x()称为信号的距离模糊函数。 由式(4)可知，x()为信号的自相关函数，对能量型信号，其傅里叶变换存在，信号自相关函数的傅里叶变换等于信号傅里叶变换幅值的平方，即信号的能谱密度。对功率型信号，当信号的自相关函数绝对可积时，信号自相关函数的傅里叶变换等于信号的功率谱密度。雷达信号在频域上占据的频带越宽，则其距离分辨力

7、越好。 x(，fd)沿着=0的截面为： x(fd)称为信号的速度模糊函数。 由式(5)可知，s(fd)为信号复包络模(实包络)平方的傅里叶变换，换句话说，速度模糊函数只与信号的幅度有关，而与信号的相位和频率调制无关。雷达信号在时域上持续宽度越大，则其速度分辨力越好。 当u(t)为离散数字信号时，u(t)=u(nTs)=u(n)(n=0，1，M-1)，Ts表示采样周期。对于给定的=m(-M+1mM-1)，设fd=fk= ，式中fs=1Ts表示采样频率，代入模糊函数定义式(3)中，得 5 FHSS信号与Costas序列信号 单位能量FHSS信号的复包络可表达为： 式中：y(k)称为跳频算子；B表示

8、信号占据的频带宽度；频隙Fb=BN，Tb表示时隙，信号占据的时域宽度T=NTb。第k个时隙发送信号的中心频率为fc+fk。为防止子脉冲频谱发生交叠，频隙应不小于子脉冲的频带宽度。 当y(k)为N阶Costas序列的放置函数时，上式表示的FHSS信号即为按照Costas序列构成方法进行FHSS的信号，这种信号具有理想的模糊函数性能。对子脉冲为恒载频的Costas序列跳频编码信号，k(t)=0；对子脉冲为LFM的Costas序列跳频编码信号， 6 多目标散射和多路径传输环境的数学模型与雷达信号设计 设雷达发射信号为u(t)，假设散射目标为多个点目标，信道为AWGN(加性白高斯噪声)信道，则接收信号

9、为： 式中：z(t)为零均值加性复高斯白噪声过程； i(t)、i(t)和wi(t)分别为第i条路径的衰减因子(以t为变量的复随机过程)、传播时延和多普勒角频率。 若散射目标不是点目标，则可将v(t)看做由连续多径分量组成的，此时：式中：c(，t)是信道时变冲激响应h(，t)的复包络(等效低通时变冲激响应)；h(，t)表示在t-时刻施加的单位冲激在t时刻的信道响应。当存在大量路径时，应用中心极限定理，c(，t)可建模为一个以t为变量的复高斯随机过程。 若信道为非时变的，则 如果将目标散射电磁波看成是信道的特性，那么对目标的探测就可等同于对信道特性的检测。因为信道是随机时变的，所以要研究能表征其统

10、计特征的物理量-散射函数。散射函数与信道等效低通冲激响应(信道冲激响应的复包络)的频移、时移自相关函数具有双傅里叶变换的关系。散射函数是时延和多普勒频率的函数，它表征了信道平均输出功率的量度，在多目标散射和多径传输的信道中，对目标的检测可归结为对信道的散射函数的检测。文献9给出了匹配滤波器输出的平均信号干扰比的数学表达式，如果信号的模糊函数具有二维函数的特性，那么匹配滤波器的输出即为信道的散射函数。所以为了正确检测出目标，要求信号的模糊函数具有二维函数的针状特性：主峰尖锐，且模糊体积很小。但事实上，这样的信号是不存在的，由模糊体积不变性知道，任何调制都不能改变模糊曲面下的总容积，该容积只决定于

11、信号的能量，而与信号的形式无关。Costas序列信号的模糊函数具有图钉状的特性：主峰尖锐，副峰很小，由图1可知，Costas序列信号的模糊体积不会减小。为了解决这个问题，可用2个信号来代替1个信号，即在同一个重复周期内，发射2个不同的Costas跳频编码序列信号，因为2个信号模糊函数的副峰不会都出现在同一个位置，所以在接收端可通过综合2个信号的信息使2个信号模糊函数的边瓣相互抵消，从而获得主峰尖锐且模糊体积很小的针状特性。 在同时使用多个同阶Costas跳频编码序列信号的多用户雷达系统中，需要减小各信号之间的互相关，多用户系统中雷达信号设计问题，限于篇幅，本文不再讨论，有关这方面的内容参见文献

12、3。 7 计算机仿真结果 有限域GF(11)上共有4个本原元，分别为2、6、7和8。其中，2和6为一对互逆的本原元，7和8为另一对互逆的本原元。 本文给出利用GF(11)上本原元=2和=6生成的2个Welch Costas序列信号的自模糊函数、自相关函数、离散傅里叶变换、互模糊函数以及采用双信号发、收的计算机仿真结果。在图2图4中：为时延，fd为多普勒频率，Tb为时隙，T为信号时域长度。子脉冲的带宽时宽乘积BT=5.76。为了满足奈奎斯特(Nyquist)取样定理，采用N=128点，用10段调频序列，调频范围分别是0.0450.09，0.1350.09，0.1350.18，0.2250.18，

13、0.2250.27，0.3150.27，0.3150.36，0.4050.36，0.4050.45，0.4950.45(归一化频率)，对这10段调频序列进行Costas跳频编码。采用N=128点，用10段恒载频序列，恒载频序列的频率分别是0.008，0.016，0.024，0.032，0.04，0.048，0.056，0.064，0.072，0.08(归一化频率)，对这10段恒载频序列进行Costas跳频编码。用双信号s1(t)和s2(t)发射、接收的仿真图。图(a)和图(b)的s1(t)采用N=128点，调频范围为00.5(归一化频率)的正斜率LFM信号；s2(t)采用N=128点，调频范围

14、为0.50(归一化频率)的负斜率LFM信号。图(a)表示：用s1(t)和s2(t)共同发射，在接收端对同一个目标(时延0=2，多普勒频率0=0.002 5)联合接收的情况。图(b)表示：用s1(t)和s2(t)共同发射，在接收端对两个不同的目标(s1(t)所对目标时延丁0=2，多普勒频率0=0.002 5，s2(t)所对目标时延1=6，多普勒频率1=0.04)联合接收的情况。图(c)和图(d)的s1(t)为用本原元=2生成的Welch Costas序列信号，含有10个子脉冲，子脉冲为LFM矩形包络的脉冲信号，每个子脉冲含有128点，调频斜率正负相间，每段序列的调频范围与图2的仿真条件相同；s2

15、(t)为用本原元=6生成的WelchCostas序列信号，含有10个子脉冲，子脉冲为LFM矩形包络的脉冲信号，每个子脉冲含有128点，调频斜率正负相间，每段序列的调频范围与图2的仿真条件相同。图(c)表示：用s1(t)和s2(t)共同发射，在接收端对同一个目标(时延0=2，多普勒频率0=0.0045)联合接收的情况。图(d)表示：用s1(t)和s2(t)共同发射，在接收端对两个不同的目标(s，(t)所对目标时延0=2，多普勒频率0=0.004 5，s2(t)所对目标时延1=8，多普勒频率1=0.036)联合接收的情况。 在信号时域长度相同的情况下，子脉冲为LFM矩形包络的Costas信号的有效

16、频谱比子脉冲为恒载频矩形包络的Costas信号的有效频谱有所扩展(扩展的倍数约为5.76倍)，距离分辨率主瓣前者比后者变得尖窄。另外，在仿真中发现，对子脉冲为恒载频的Costas信号，改变频隙，对信号的模糊函数性能影响较大，图中仿真结果对应的频隙等于时隙的倒数，即归一化频隙为11280.008。由于速度模糊函数为信号复包络模(实包络)平方的傅里叶变换，与信号的相位和频率调制无关，而信号包络为矩形脉冲(门函数)，因此速度模糊函数的辐频特性具有取样函数的模|sa(x)| =|sin xx|的形状。 采用正斜率和负斜率的两种LFM信号，对同一个目标这两种信号的边瓣将相互抵消，从而获得主峰尖锐、且模糊体积很小的针状特性。但对两个不同的目标，边瓣将不能相互抵消，引起虚假信号。采用两种Costas序列信号，对同一个目标这两种Costas序列信号的边瓣将相互抵消，从而获得主峰尖锐、且模糊体积很小的针状特性，并且对两个不同的目标，边瓣也能相互抵消，不会引起虚假信号。 8 结束语 利用Co