水声OFDM系统中卷积码译码设计及其DSP实现

水声信道可用频带窄，具有强多途、高噪声，同时伴有明显的随机包络起伏。这些特性使水声信道具有典型的频率选择性衰落和时间选择性衰落的特点。信号经过复杂的水声信道会出现随机错误和突发错误。因此需要加入纠错编码与交织相结合的时频编码方案进行错误纠正[1]。在水声通信中，纠错编码是不可或缺的技术之一。
目前通常采用的编码方案有卷积码、Turbo码和LDPC码等。后两者具有接近香农限的性质，在水声OFDM系统中都表现出很好的纠错特性。在编码增益上，卷积码比Turbo码和LDPC码分别低3.8 dB和4 dB[2]，但其实现简单，实时性好。因此从水声OFDM系统硬件实现的角度综合考虑，本文采用卷积码作为系统的差错控制方案。

为提高系统性能和实时性，本文在方案中采用了卷积交织和Viterbi软译码等方法。经过大量仿真及水池实验，确定了交织和编译码的最佳参数。在TMS320DM642上实现时，通过采用蝶形运算等优化方法减少了译码复杂度，提高了运算速度。最后通过海洋实验验证了其性能。

1 水声OFDM通信系统原理及编译码方案确定

图1为水声OFDM基带通信系统图[3]。由于水声信道存在严重的频率选择性衰落，采用卷积编码实现各个子载波上的纠错。为了抵抗由于时变带来的突发错误，采用交织技术将突发错误打散为随机错误。对卷积交织后的串行数据流进行DQPSK映射、IFFT、插入循环前缀以及并串转换后形成OFDM符号，加入同步前导符号后经水声信道发送出去。在接收端提取同步信息后，经过一系列与发送端相反的逆过程实现OFDM解调，解调后的数据流再经过解交织和Viterbi译码后还原数据信息。

卷积码码型的选择主要根据系统的相干带宽和相干时间来确定。为了提高水声信道的频带利用率，通常采用1/2码率的卷积码。约束长度的选择以不同信噪比条件下卷积码性能为参考。仿真和水池实验结果表明：信噪比较高时，约束长度越大卷积码性能越好；信噪比较低时，约束长度越小卷积码性能越好。由于水声信道的信噪比较低，本文最终确定使用生成矩阵为[7，5]的(2，1，3)卷积码，从复杂度与性能上折中考虑，译码回溯深度定为16，Viterbi软译码量化级数定为3bit均匀量化。交织长度和深度的确定要综合考虑实时性和译码性能，仿真结果表明交织长度为4个OFDM符号时，最佳交织深度为96。

该系统在Matlab平台上进行水池实验，实验结果如图2所示。实验中对采用三个不同方案的OFDM系统进行了性能比较，以示波器观测到的接收端信号幅度作为SNR大小的参照。实验中噪声幅度约为150 mV。从图中可以明显看出Viterbi软译码性能最佳。

2 DSP优化实现

2.1 交织的优化

硬件实现时需考虑所使用的交织器类型。常见的交织器主要有块交织器、卷积交织器和随机交织器。在相同的交织长度和交织深度下，卷积交织器与块交织相比延时小，与随机交织相比硬件复杂度低，是较实用的类型。
根据卷积交织原理[4]，在DSP上实现时，如果采用一般算法，则需要用到二维数组和多重嵌套循环，DSP处理将耗费大量时间和存储空间。交织的实质是改变原始序列中数据的顺序，一旦交织深度和宽度确定，其数据的顺序改变规律即可相应确定。因此为了节省处理时间和片内存储空间，本文采用查表的方法实现卷积交织。

2.2 Viterbi译码器的优化

实现Viterbi译码可以分为初始化、度量值更新和回溯跟踪三部分，其中最复杂的是第二部分。本文从两个方面进行了优化：一是对计算量最多的蝶形运算采用宏定义的方式，并在记录度量长度时采用双数组，能够加快译码计算速度；二是定义了一种软判决量度，减少计算复杂度。

2.2.1 蝶形运算

在Viterbi译码程序设计中，使用对称的蝶形运算实现。(2，1，3)卷积码的蝶形结构如图3所示。每一个蝶形包括当前状态为i(i=0，1)和i+2的两个节点的加_比较_选择运算，它们的0和1分支在篱笆图的下一个节点合并。为了加快运算速度，程序设计上对蝶形运算采用宏定义的方式，定义两个蝶形运算的宏[5]。一个设当前度量为old_pm[ ]， 经过分支度量计算后得到下一状态度量new_pm[ ]。另一个则从下一状态new_pm[ ]开始，经过分支度量计算后得到old_pm[ ]。

通过宏定义，大量的蝶形运算在编译时间里进行宏展开代入表达式，这样运算速度会比子程序设计快， 而且不需要大量的寄存器变量。程序中还定义了两个记录度量长度的数组，循环进行数组更新，不需要再定义中间变量或指针进行数组的更新，可以提升DSP运算速度。

2.2.2 软判决量度

计算接收到的比特与理论输出比特的距离度量公式为：bm[i0 i1]=metric[i0 r0]+metric[i1 r1]。其中i0、i1为理论输出数