Costas序列模糊函数仿真

---

前言

本文对 Costas 序列进行讲解，并仿真其模糊函数。

---

一、Costas 序列介绍

任意每行每列有且仅有 1 个元素等于 1、其余元素都为 0 的 nxn 阶矩阵称为置换矩阵。Costas 序列是一类特殊的置换矩阵，它与自身任意方向的平移副本之间都至多有 1 个元素 “1” 重合，如矩阵 A 所示（其序列表示为 [4,1,6,7,5,8,3,2]）。

Costas 序列定义为 $N_c^2$ 阶置换矩阵 $A=a_{ij}$（$a_{i,j}$$\in${0,1}，$1\le i\le N_c$，$1\le j\le N_c$），其中 $N_c$ 表示 Costas 序列的长度。

Costas 信号是一种频率跳变信号，是跳频信号中的一种特殊形式。其中跳频信号由若于不同频率的 CW 信号组成，不同频率的 CW 波时域不重叠，即时分复用，而 Costas 信号是跳频频率由 Costas 序列指定的跳频信号。

Costas 信号由不同频率的 CW 子脉冲拼接而成，其带通形式如下：

其中：$Re[\cdot ]$ 表示取实部，$rect(\cdot )$ 为矩形窗函数（在 0~1 等于 1，其余位置等于 0）；$N$ 为 Costas 编码序列的长度，$T_{sp}$ 为单个 CW 子脉冲的持续时间；$f_n=f_c+C_n\Delta f$ 为经过 Costas 序列编码后的各 CW 子脉冲的频率，$C_n$ 为均值调整到 0 的 Costas 序列，$f_c$ 为载波频率，$\Delta f$为单位跳频量。

二、模糊函数

模糊函数是雷达探测波形分析的重要工具，通过对信号波形的模糊函数分析，可以得到信号波形的距离分辨率、多普勒分辨率及多普勒容限特性。

为了分析在不同的时延和多普勒频移条件下接收信号与匹配滤波器的不匹配程度，定义了模糊度函数。

连续时间信号模糊函数的定义为：
$$\chi (\tau ,f_d)=\frac{1}{E}{\int }_{-\infty }^{\infty }b(t)b^{\ast }(t-\tau )e^{j2\pi f_{d}t}dt$$

• 式中，E为信号的总能量；

离散时间序列的模糊函数表示为：
$$\chi (m,k_d)=\frac{1}{E_c}\sum _n{e}_n{e}_{n-m}^{\ast }{e}^{j\frac{2\pi }{N}{k}_{d}n}$$

• 式中，$m=f_s\times \tau$，$f_s$ 为采样率；
• $k_d=\frac{f_d\times f_s}{N}$，N为采样点数

三、Costas 序列模糊函数仿真

1、MATLAB 核心代码

tl = linspace(0,1,Np);      % 模糊图的时延
fd = linspace(-150,150,Np); % 模糊图的频率

ambi = abs(xcorr2(bsxfun(@times, Costas, exp(1j*2*pi*fd'*tl)), Costas)); % 对信号做共轭相乘互相关（计算模糊图）

2、仿真结果

①、Costas 模糊函数图

②、Costas 距离模糊函数图

③、Costas 速度模糊函数图

四、资源自取

链接：Costas序列模糊函数仿真

---

我的qq：2442391036，欢迎交流！

---