巴克码及其应用（VDES）

1 巴克码

1.1 简介

巴克码是R.H.巴克于50年代初提出的一种特殊的二进制码组。它是一个非周期序列，并且具有非常良好的自相关特性，可用于数据帧的帧头。其定义如下：一个$n$位的巴克码序列为 { x 1 , x 2 , x 3 , … , x n } \{ {x_1},{x_2},{x_3}, \ldots ,{x_n}\} {x1​,x2​,x3​,…,xn​}，其中$x_i$的可能取值为+1或者-1，其自相关函数（非周期）$\gamma \left(j\right)$满足式：

$$\gamma \left(j\right)=\sum _{i=1}^{n-j}{x}_i{x}_{i+j}=\{\begin{array}{ll}n& \text{j=0}\\ 0,\pm 1& \text{0<j<n}\\ 0,& \text{j⩾n}\end{array}$$

目前，已知的巴克码序列共有9组，其中序列的最大长度为13位，并且13位的巴克码序列为 { + , + , + , + , + , − , − , + , + , − , + , − , + } \{ +,+,+,+,+,-,-,+,+,-,+,-,+\} {+,+,+,+,+,−,−,+,+,−,+,−,+}。其中，‘+’表示+1、‘-’表示-1。

1.2 特性

13位的巴克码序列自相关结果如图所示（自相关函数为偶函数，具有对称性）：

相应的Matlab代码如下：
barker = [1 1 1 1 1 -1 -1 1 1 -1 1 -1 1];	%13位巴克码序列

[y1, x1] = xcorr(barker);	%计算13位巴克码序列的自相关值
        
figure(1);
plot(x1, abs(y1), '-o');
% axis tight;
axis([-13 13 0 14]);
xlabel('\itj', 'FontName', 'Times New Roman');
ylabel('\it\gamma(j)', 'FontName', 'Times New Roman');

该序列的非周期自相关函数只有4种可能的取值：最大值为$n$、最小值为0或者$\pm 1$，其自相关函数的主峰和旁瓣均为等腰三角形，底边宽度为$2T_b$（定义码元周期为$T_b$），并且主峰峰值是旁瓣峰值的$n$倍。对于13个比特长度的巴克码序列，其主峰峰值则是旁瓣峰值的13倍。

2 VDES中的应用

2.1 简介

VDES（甚高频数据交换系统）可以说是AIS（船舶自动识别系统）的加强和升级版系统，在整合了所有AIS功能的基础上，还增加了ASM（特定应用报文信息）和VDE（甚高频数据交换）功能。

因为巴克码序列本身长度过短，实际工程应用中受到很大限制，所以需要寻求新的方法以增加巴克码序列的长度。VDES中的ASM-TER和VDE-TER数据帧的训练序列均为双巴克码序列，即“1 1111100110101 0000011001010”（其中，‘1’表示+1、‘0’表示-1），该序列是由‘1’+“13位巴克码”+“13位反巴克码”组成的一组双巴克码序列。

2.2 特性

双巴克码序列自相关结果如图所示：

相应的Matlab代码如下：
double_barker = [1 1 1 1 1 1 -1 -1 1 1 -1 1 -1 1 -1 -1 -1 -1 -1 1 1 -1 -1 1 -1 1 -1];	%双巴克码序列

[y1, x1] = xcorr(double_barker);	%计算双巴克码序列的自相关值
        
figure(2);
plot(x1, abs(y1), '-o');
axis tight;
xlabel('\itj', 'FontName', 'Times New Roman');
ylabel('\it\gamma(j)', 'FontName', 'Times New Roman');

双巴克码序列自相关函数的主峰峰值经过27个符号的累积，能量为27。此外，双巴克码序列自相关结果中还包含了两个能量为12的副瓣，副瓣是由双巴克码序列中13位巴克码序列的正码和反码对齐导致，具体的对齐情况如图所示：

其中，状态（1）和（3）所示对齐情况计算所得到的自相关值分别对应于两个副瓣，状态（2）所示对齐情况计算所得到的自相关值对应于主峰。

双巴克码序列和13位巴克码序列互相关结果如图所示：

相应的Matlab代码如下：
barker = [1 1 1 1 1 -1 -1 1 1 -1 1 -1 1];	%13位巴克码序列
double_barker = [1 1 1 1 1 1 -1 -1 1 1 -1 1 -1 1 -1 -1 -1 -1 -1 1 1 -1 -1 1 -1 1 -1];	%双巴克码序列

[y1, x1] = xcorr(barker, double_barker);	%计算双巴克码序列和13位巴克码序列的互相关值
    
figure(3);
plot(x1, abs(y1), '-o');
axis tight;
xlabel('\itj', 'FontName', 'Times New Roman');
ylabel('\it\gamma(j)', 'FontName', 'Times New Roman');