一种叠加导频的M元多相序列扩频水声通信方法

一种叠加导频的m元多相序列扩频水声通信方法
技术领域
1.本发明属于水声通信领域，具体来说，涉及多相扩频码序列、叠加导频进行信道估计、最大比合并频域均衡的一种叠加导频的m元多相序列扩频水声通信方法。

背景技术：

2.随着海洋开发和国防建设的发展，利用水声通信技术传递信息的需求大大增加。水声信道是一种带宽有限信道，由于传播损耗的影响，发射信号在高频段的衰减非常大，在低频段容易受到环境噪声和发射换能器的影响而不能被使用；而且水声信道有复杂的多径结构、受多普勒频移严重影响，这些特点使得可靠的水声通信成为研究难点。
3.扩频通信是一种信息处理传输技术，在发送端采用伪随机码(即pn码)作为扩频码，将待传输的数据信息的频谱扩展调制到一个比信息带宽宽得多的频带上，接收端采用与发射端相同的扩频码对接收信号进行解扩解调。扩频通信通过扩频码实现频谱扩展，普通的二元伪随机序列在扩频序列周期选取上受到一定的限制，序列个数也不多，这些缺点导致二元伪随机序列不能满足实际要求。
4.一般的扩频技术其伪随机码速率很高，射频带宽非常宽，在频带受限的情况下，如水声通信中难以满足系统要求，这时就可以采用m元扩频技术。m元扩频通信因其具有抗多径干扰以及隐蔽性等优点，在水声通信中得到广泛使用。然而，常规的m元扩频技术，接收端需要m个相关器，当m较大时，接收端的计算量较大，不易处理。

技术实现要素：

5.为了克服现有技术的问题，本发明提出一种叠加导频的m元多相序列扩频水声通信方法，运用具有优秀相关特性的多相pecan序列作为扩频码序列，该序列的复杂度高、易于生成，可以提高频带利用率，满足工程上的需求；本发明中，在发射端将作为导频信号与数据符号相叠加，比起传统的内插导频方法，可以有效的节省系统带宽；在m元扩频的基础上，本发明采用循环移位技术(csk)，此时在接收端只需与扩频码数目相同的相关器进行相关操作，利用序列的相关性进行通信，只使用fft/ifft操作，降低了系统的复杂度。
6.本发明的技术方案为：
7.一种叠加导频的m元多相序列扩频水声通信方法，包括以下步骤：
8.步骤1：将经串并转换后的待通信数据信息比特平均分给n个单载波数据块，每个单载波块所含的信息量为(j-1)
×
(log2(m1)+log2(m2)+log2(m3))；其中(j-1)是每个单载波数据块并行发射的扩频序列个数，m1为m元扩频的调制阶数，m2为循环移位的调制阶数，m3为相位调制的阶数；其中第k个单载波数据块uk表示为：
9.uk＝[u
k,1
,u
k,2
,
…
,u
k,j-1
],u
k,j
∈uk,j＝1,2,
…
,j-1
[0010]uk,j
为第k个单载波数据块并行发射的第j个扩频序列；
[0011]
步骤2：对于每个单载波数据块，从长度为l，数目为nc的pecan扩频序列集合p中随机选取其中一个扩频序列pn∈p,n＝0,1,
…
,n
c-1作为导频信号，再剩余的n
c-1个扩频序列
随机选择(j-1)个扩频序列，其中
[0012][0013]
步骤3：对步骤2得到的(j-1)个扩频序列中的每个扩频序列进行循环移位和相位调制，并叠加形成数据块，其中第k个单载波块对应的数据块fk为
[0014][0015]
步骤4：将步骤2选择的导频信号pn，与步骤3生成的数据块fk进行叠加，在叠加处的信号为：
[0016]
sk＝pn+fk；
[0017]
步骤5：将步骤4得到的叠加处的信号sk尾部长为l
cp
的部分复制到信号前部，得到第k个数据块的发送信号其中l
cp
为设定的循环前缀长度；
[0018]
步骤6：对于每个数据块，均重复步骤2到步骤5，得到系统的发射端信号s：
[0019][0020]
对发射端信号s进行并串转换后，通过水声信道，传输到接收端，在接收端进行串并转换，得到n个接收单载波数据块，与步骤1中的n个单载波数据块一一对应；对于每个接收单载波数据块，将信号去除循环前缀后得到时域接收信号y，并采用基于叠加导频的系统频域传输模型表示为：
[0021]
yk＝hksk+w
[0022]
其中，yk、hk、sk分别为第k个接收单载波数据块对应的时域接收信号y、时域信道响应h、时域发射信号s分别作l点fft得到的相应频域值，w为信道的加性高斯白噪声的频域形式；
[0023]
步骤7：将步骤6得到的信号yk和导频信号pn作相关处理，并通过ifft得到
[0024][0025]
其中表示相关估计得到的第k个接收单载波数据块的时域信道响应；fft
*
是傅里叶变换fft的复共轭；
[0026]
步骤8：对步骤7得到的时域信道响应进行信道均衡，使用最大比合并mrc频域均衡器，得到频域接收信号yk经过均衡处理后得到的频域信号为zk＝ckyk；ck为第k个接收单载波数据块的mrc均衡器矩阵；
[0027]
步骤9：对步骤8得到的频域信号zk进行解扩和解调，利用相关fft处理得到扩频码的索引位置以及循环移位位置利用循环移位位置对第k个接收单载波数据块进行相位解调；
[0028]
步骤10：对于n个接收单载波数据块均重复步骤7到步骤9，从而完成接收端的解调过程，对解调后的数据进行并串转换输出。
[0029]
进一步的，步骤2具体过程为：
[0030]
从nc个pecan扩频序列中随机选取其中一个扩频序列pn∈p,n＝0,1,
…
,n
c-1作为
导频信号，剩余的n
c-1个扩频序列分成(j-1)个子集，每个子集包含m1个pecan扩频序列；再从每个子集中随机选择一个pecan扩频序列，共得到(j-1)个扩频序列。
[0031]
进一步的，步骤3具体过程为：
[0032]
对于pecan扩频序列pj，将待通信数据信息比特中的b2＝log2(m2)个信息比特映射到pj上，得到序列该序列表示扩频序列pj以b2为移位距离向右循环移位所得到的序列：
[0033][0034]
再将待通信数据信息比特中的b3＝log2(m3)个psk调制符号映射到序列上，得到序列
[0035]
将步骤2得到的(j-1)个扩频序列对应的序列相叠加形成第k个单载波块对应的数据块fk[0036][0037]
进一步的，步骤7中时域信道估计的结果表示为：
[0038][0039]
其中，γ是设定的阈值。
[0040]
进一步的，均衡器的抽头系数表示为：
[0041][0042]
其中，c
t
为第k个接收单载波数据块的mrc均衡器矩阵ck中的第t个均衡器抽头系数。
[0043]
进一步的，步骤9中，扩频码的索引位置为：
[0044][0045]
为第k个接收单载波数据块对应的(j-1)个pecan扩频序列中的第j个pecan扩频序列。
[0046]
进一步的，步骤9中，循环移位位置为：
[0047][0048]
其中表示根据扩频码的索引位置从第k个接收单载波数据块对应的(j-1)个pecan扩频序列中到的pecan扩频序列。
[0049]
有益效果
[0050]
本发明所提出的一种叠加导频的m元多相序列扩频水声通信方法，使用多相序列(pecan)扩频码序列集，由于序列产生方式具有随机性，且序列本身具有强自相关特性，可
以有效增强系统的隐蔽性，提高系统的抗干扰能力。将数据符号与导频信号相叠加，避免导频符号单独占据带宽资源，从而节省带宽；循环前缀单载波块结构可以消除由于多径效应产生符号间干扰。使用最大比合并mrc频域均衡器，可以有效的对抗信道的衰落，并且比起最小均方误差mmse均衡器在误比特率上有更好的表现形式。在系统接收端处，解扩解调时只采用了fft/ifft操作，能够明显地降低系统的计算量，具有很好的应用前景。
[0051]
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0052]
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
[0053]
图1：水声通信整体框架图。
[0054]
图2：多相序列(pecan)星座图。幅值为1，序列长度l＝1024，序列的长度等于序列对应的相位数，即有1024个相位数。
[0055]
图3：在扩频码序列长度l＝1024的情况下，pecan扩频码序列的时域自相关、时域互相关性以及频域自相关性曲线。(a)pecan时域自相关性，(b)pecan时域互相关性，(c)pecan频域自相关性；从图中可以看出，时/频域自相关函数的主峰均很尖锐，具有强自相关性，使得信号易于检测；互相关函数几乎为0，可以提高系统抗干扰能力。
[0056]
图4：循环前缀单载波块结构图。一帧数据信息中的循环前缀单载波块结构示意图。
[0057]
图5：pecan循环移位扩频解码图。
[0058]
图6：在m1＝4，m3＝4下，并行发射的扩频码序列个数j＝2、3、4、5下的ber与信噪比(snr)的关系曲线仿真结果。从中可以看出mrc均衡器的误码率性能好于mmse均衡器，在使用mrc均衡时，j＝2、3、4、5分别在snr＝-5、-2、0、3db下达到10-4
的误码率，随着j的增加，误码率性能越差，在实验中考虑使用j＝2、3。
[0059]
图7：在j＝3，m1＝4下，多进制相位调制阶数m3＝2、4、8下的ber与信噪比(snr)的关系曲线仿真结果。从中可以看出在误码率为10-4
时，使用mrc均衡器比mmse均衡器的性能提升了约0.8db；以及m3的改变对系统误码率的影响很小，在实验中，考虑使用m3＝4，即qpsk调制。
[0060]
图8：在j＝2，m3＝4下，扩频码索引调制阶数m1＝2、4、8下的ber与信噪比(snr)的关系曲线仿真结果，从中可以看出在误码率为10-4
时，使用mrc均衡器比mmse均衡器的性能提升了约0.8db；以及改变m1的数值对系统误码率的影响很小。在湖上实验考虑使用m1＝4、8。
[0061]
图9：数据信息包结构图。
[0062]
图10：发射信号时域波形图。
[0063]
图11：醴泉湖多径信道结构图。
[0064]
图12：信道的时变冲激响应图。
[0065]
图13：醴泉湖所有通道部分接收信号波形图。
具体实施方式
[0066]
本发明提出一种叠加导频的m元多相序列扩频水声通信方法，首先运用具有理想脉冲自相关特性的多相pecan序列作为扩频码序列，由于序列的产生具有随机性以及不可预测性，可以有效提高系统的抗干扰能力；其次在发射端，将用作导频信号的pecan序列，与作为数据的pecan序列相叠加，可以有效节省带宽；并构造循环前缀单载波块结构，避免多径衰落信道引起的符号间干扰(isi)；同时使用最大比合并(maximal ratio combining,mrc)频域均衡器进行接收端信道均衡，可以有效对抗信道衰落；最后在接收端解扩解码时，只使用fft/ifft操作，可以有效降低计算复杂度。
[0067]
本发明的具体步骤为：
[0068]
步骤1：将经串并转换后的待通信数据信息比特平均分给n个单载波数据块，每个单载波块所含的信息量为(j-1)
×
(log2(m1)+log2(m2)+log2(m3))；其中(j-1)是每个单载波数据块并行发射的扩频序列个数，m1为m元扩频的调制阶数，m2为循环移位的调制阶数，m3为相位调制的阶数；其中第k个单载波数据块uk表示为：
[0069]
uk＝[u
k,1
,u
k,2
,
…
,u
k,j-1
],u
k,j
∈uk,j＝1,2,
…
,j-1
[0070]uk,j
为第k个单载波数据块并行发射的第j个扩频序列；
[0071]
步骤2：对于每个单载波数据块，从长度为l，数目为nc的pecan扩频序列集合p中随机选择(j-1)个扩频序列，其中
[0072][0073]
具体过程为：从nc个pecan扩频序列中随机选取其中一个扩频序列pn∈p,n＝0,1,
…
,n
c-1作为导频信号，剩余的n
c-1个扩频序列分成(j-1)个子集，每个子集包含m1个pecan扩频序列；再从每个子集中随机选择一个pecan扩频序列，共得到(j-1)个扩频序列；
[0074]
步骤3：对步骤2得到的(j-1)个扩频序列中的每个扩频序列进行循环移位和相位调制；
[0075]
具体过程为，对于pecan扩频序列pj，将待通信数据信息比特中的b2＝log2(m2)个信息比特映射到pj上，得到序列该序列表示扩频序列pj以b2为移位距离向右循环移位所得到的序列：
[0076][0077]
再将待通信数据信息比特中的b3＝log2(m3)个psk调制符号映射到序列上，得到序列
[0078]
将步骤2得到的(j-1)个扩频序列对应的序列相叠加形成第k个单载波块对应的数据块fk[0079][0080]
步骤4：将步骤2选择的导频信号pn，与步骤3生成的数据块fk进行叠加，在叠加处的信号为：
[0081]
sk＝pn+fk[0082]
步骤5：将步骤4得到的叠加处的信号sk尾部长为l
cp
的部分复制到信号前部，得到第k个数据块的发送信号其中l
cp
为设定的循环前缀长度；
[0083]
步骤6：对于每个数据块，均重复步骤2到步骤5，得到系统的发射端信号s：
[0084][0085]
对发射端信号s进行并串转换后，通过水声信道，传输到接收端，在接收端进行串并转换，得到n个接收单载波数据块，与步骤1中的n个单载波数据块一一对应；对于每个接收单载波数据块，将信号去除循环前缀后得到时域接收信号y，并采用基于叠加导频的系统频域传输模型表示为：
[0086]
yk＝hksk+w
[0087]
其中，yk、hk、sk分别为第k个接收单载波数据块对应的时域接收信号y、时域信道响应h、时域发射信号s分别作l点fft得到的相应频域值，w为信道的加性高斯白噪声的频域形式；
[0088]
步骤7：将步骤6得到的信号yk和导频信号pn作相关处理，并通过ifft得到
[0089][0090]
其中表示相关估计得到的第k个接收单载波数据块的时域信道响应；fft
*
是傅里叶变换fft的复共轭；并由于水声信道具有稀疏性，时域信道估计的结果可以表示为：
[0091][0092]
其中，γ是设定的阈值；
[0093]
步骤8：对步骤7得到的时域信道响应进行信道均衡，使用最大比合并mrc频域均衡器，均衡器的抽头系数表示为：
[0094][0095]
其中，c
t
为第k个接收单载波数据块的mrc均衡器矩阵ck中的第t个均衡器抽头系数；
[0096]
将频域接收信号yk经过均衡处理后得到的频域信号为zk＝ckyk；
[0097]
步骤9：对步骤8得到的频域信号zk进行解扩和解调，利用相关fft处理得到扩频码的索引位置以及循环移位位置：
[0098]
其中扩频码的索引位置为：
[0099][0100]
为第k个接收单载波数据块对应的(j-1)个pecan扩频序列中的第j个pecan扩频序列；
[0101]
循环移位位置为：
[0102]
[0103]
其中表示根据扩频码的索引位置从第k个接收单载波数据块对应的(j-1)个pecan扩频序列中到的pecan扩频序列；
[0104]
利用循环移位位置对第k个接收单载波数据块进行相位解调；
[0105]
步骤10：对于n个接收单载波数据块均重复步骤7到步骤9，从而完成接收端的解调过程，对解调后的数据进行并串转换输出。
[0106]
下面详细描述本发明的实施例，所述实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
[0107]
步骤1：本发明的系统框图如图1所示。将经串并转换后的待通信数据信息比特平均分给n个单载波数据块，每个单载波块所含的信息量为(j-1)
×
(log2(m1)+log2(m2)+log2(m3))；其中(j-1)是每个单载波块并行发射的扩频序列个数，m1为m元扩频的调制阶数，m2为循环移位的调制阶数，m3为相位调制的阶数；其中第k个单载波块uk表示为：
[0108]
uk＝[u
k,1
,u
k,2
,
…
,u
k,j-1
],u
k,j
∈uk,j＝1,2,
…
,j-1
[0109]uk,j
为第k个单载波数据块并行发射的第j个扩频序列。
[0110]
步骤2：对于每个单载波数据块，从长度为l，数目为nc的pecan扩频序列集合p中随机选择(j-1)个扩频序列，其中
[0111][0112]
pecan多相序列的星座图如图2所示，pecan扩频序列的相关性如图3所示。
[0113]
具体过程为：从nc个pecan扩频序列中随机选取其中一个扩频序列pn∈p,n＝0,1,
…
,n
c-1作为导频信号，剩余的n
c-1个扩频序列分成(j-1)个子集，每个子集包含m1个pecan扩频序列；再从每个子集中随机选择一个pecan扩频序列，共得到(j-1)个扩频序列。
[0114]
步骤3：对步骤2得到的(j-1)个扩频序列中的每个扩频序列进行循环移位和相位调制；
[0115]
具体过程为，对于pecan扩频序列pj，将待通信数据信息比特中的b2＝log2(m2)个信息比特映射到pj上，得到序列该序列表示扩频序列pj以b2为移位距离向右循环移位所得到的序列：
[0116][0117]
再将待通信数据信息比特中的b3＝log2(m3)个psk调制符号映射到序列上，得到序列
[0118]
将步骤2得到的(j-1)个扩频序列对应的序列相叠加形成第k个单载波块对应的数据块fk[0119][0120]
基于m元扩频的csk相位调制的数据率rb可以表示为：
[0121]
[0122]
步骤4：将步骤2选择的导频信号pn，与步骤3生成的数据块fk进行叠加，在叠加处的信号为：
[0123]
sk＝pn+fk。
[0124]
步骤5：将步骤4得到的叠加处的信号sk尾部长为l
cp
的部分复制到信号前部，得到第k个数据块的发送信号循环前缀单载波块结构示意图如图4所示，其中l
cp
为设定的循环前缀长度。
[0125]
步骤6：对于每个数据块，均重复步骤2到步骤5，得到系统的发射端信号s：
[0126][0127]
对发射端信号s进行并串转换后，通过水声信道，传输到接收端，在接收端进行串并转换，得到n个接收单载波数据块，与步骤1中的n个单载波数据块一一对应；对于每个接收单载波数据块，将信号去除循环前缀后得到时域接收信号y，并采用基于叠加导频的系统频域传输模型表示为：
[0128]
yk＝hksk+w
[0129]
其中，yk、hk、sk分别为第k个接收单载波数据块对应的时域接收信号y、时域信道响应h、时域发射信号s分别作l点fft得到的相应频域值，w为信道的加性高斯白噪声的频域形式。
[0130]
步骤7：将步骤6得到的信号yk和导频信号pn作相关处理，并通过ifft得到
[0131][0132]
其中表示相关估计得到的第k个接收单载波数据块的时域信道响应；fft
*
是傅里叶变换fft的复共轭；并由于水声信道具有稀疏性，时域信道估计的结果可以表示为：
[0133][0134]
其中，γ是设定的阈值。
[0135]
步骤8：对步骤7得到的时域信道响应进行信道均衡，使用最大比合并mrc频域均衡器，均衡器的抽头系数表示为：
[0136][0137]
其中，c
t
为第k个接收单载波数据块的mrc均衡器矩阵ck中的第t个均衡器抽头系数；
[0138]
将频域接收信号yk经过均衡处理后得到的频域信号为zk＝ckyk。
[0139]
步骤9：对步骤8得到的频域信号zk进行解扩和解调，利用相关fft处理得到扩频码的索引位置以及循环移位位置：
[0140]
其中扩频码的索引位置为：
[0141]
[0142]
为第k个接收单载波数据块对应的(j-1)个pecan扩频序列中的第j个pecan扩频序列；
[0143]
循环移位位置为：
[0144][0145]
其中表示根据扩频码的索引位置从第k个接收单载波数据块对应的(j-1)个pecan扩频序列中到的pecan扩频序列。
[0146]
在计算机仿真时，l＝1024，b2＝10bit，对pecan序列进行循环相关得到循环移位扩频解码结果，解码结果如图5所示。从图中可以看出相关峰明显，从解码结果可知该序列循环移位57次，对应的信息为0000111001。
[0147]
利用循环移位位置对第k个接收单载波数据块进行相位解调；
[0148]
对于n个接收单载波数据块均重复步骤7到步骤9，从而完成接收端的解调过程，对解调后的数据进行并串转换输出。
[0149]
在计算机仿真中，对所提出的方案进行性能分析测试。将各种参数进行比较，误比特率ber与信噪比snr的关系如图6、图7、图8，仿真参数如表1所示。
[0150]
表1：计算机仿真系统参数
[0151]
并行发射的扩频码个数j2，3，4，5m12，4，8csk调制阶数m21024相位调制阶数m32，4，8码片时长tc0.21ms循环前缀长度l
cp
256单载波块个数n10扩频序列长度l1024
[0152]
为了进一步验证所提方案的性能，在醴泉湖进行湖上实验，由编码效率为1/2的二进制分组码(5,7)编码器对信息比特进行卷积编码、随机交织后，经过一个滚降因子为0.25的平方根升余弦脉冲整形滤波器后，形成发射数据。
[0153]
发射的数据信息包结构如图9所示，一个数据块包括一个500ms的同步信号，100ms的零填充信号，数据信号后跟一个100ms的线性调频信号，最后跟一个100ms的零填充周期。数据信息包由10个已编码的数据块组成。
[0154]
在实验中，连续发送四种参数的信号，每个信号持续时间约70秒，发射信号时域波形图如图10所示。实验中发射信号的参数设置如表2所示。
[0155]
表2：湖上实验系统参数
[0156][0157]
发射船与接收船相距约500m，发射换能器在大约8m的深度进行信号发射，接收端是一个线性八元阵，每个水听器相距0.5m，最底端的水听器放置深度大约在12m处，醴泉湖的水声信道多径结构较为复杂，如图11所示，信道冲激响应(cir)接收信号部分通道示意图如图12所示，按照其中，ps是信号的功率，pn是噪声的功率，估计出的信噪比snr大约在12db左右。实验数据处理的结果如表3所示。
[0158]
表3：湖上实验数据处理结果
[0159][0160]
本发明使用多相序列(pecan)扩频码序列集，可以有效增强系统的隐蔽性，提高系统的抗干扰能力。将数据符号与导频信号相叠加，避免导频符号单独占据带宽资源，从而节省带宽；循环前缀单载波块结构可以消除由于多径效应产生符号间干扰。使用最大比合并mrc频域均衡器，可以有效的对抗信道的衰落，并且比起最小均方误差mmse均衡器在误比特率上有更好的表现形式。在系统接收端处，解扩解调时只采用了fft/ifft操作，能够明显地降低系统的计算量，具有很好的应用前景。
[0161]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

技术特征：

1.一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：包括以下步骤：步骤1：将经串并转换后的待通信数据信息比特平均分给n个单载波数据块，每个单载波块所含的信息量为(j-1)
×
(log2(m1)+log2(m2)+log2(m3))；其中(j-1)是每个单载波数据块并行发射的扩频序列个数，m1为m元扩频的调制阶数，m2为循环移位的调制阶数，m3为相位调制的阶数；其中第k个单载波数据块u
k
表示为：u
k
＝[u
k,1
,u
k,2
,
…
,u
k,j-1
],u
k,j
∈u
k
,j＝1,2,
…
,j-1u
k,j
为第k个单载波数据块并行发射的第j个扩频序列；步骤2：对于每个单载波数据块，从长度为l，数目为n
c
的pecan扩频序列集合p中随机选取其中一个扩频序列p
n
∈p,n＝0,1,
…
,n
c-1作为导频信号，再剩余的n
c-1个扩频序列随机选择(j-1)个扩频序列，其中n
c
＝(j-1)
×
m1+1，p
j
＝{p
j
(i)}|
i＝0,1,
…
,l-1
∈p；步骤3：对步骤2得到的(j-1)个扩频序列中的每个扩频序列进行循环移位和相位调制，并叠加形成数据块，其中第k个单载波块对应的数据块f
k
为步骤4：将步骤2选择的导频信号p
n
，与步骤3生成的数据块f
k
进行叠加，在叠加处的信号为：s
k
＝p
n
+f
k
；步骤5：将步骤4得到的叠加处的信号s
k
尾部长为l
cp
的部分复制到信号前部，得到第k个数据块的发送信号其中l
cp
为设定的循环前缀长度；步骤6：对于每个数据块，均重复步骤2到步骤5，得到系统的发射端信号s：对发射端信号s进行并串转换后，通过水声信道，传输到接收端，在接收端进行串并转换，得到n个接收单载波数据块，与步骤1中的n个单载波数据块一一对应；对于每个接收单载波数据块，将信号去除循环前缀后得到时域接收信号y，并采用基于叠加导频的系统频域传输模型表示为：y
k
＝h
k
s
k
+w其中，y
k
、h
k
、s
k
分别为第k个接收单载波数据块对应的时域接收信号y、时域信道响应h、时域发射信号s分别作l点fft得到的相应频域值，w为信道的加性高斯白噪声的频域形式；步骤7：将步骤6得到的信号y
k
和导频信号p
n
作相关处理，并通过ifft得到作相关处理，并通过ifft得到其中表示相关估计得到的第k个接收单载波数据块的时域信道响应；fft
*
是傅里叶变换fft的复共轭；步骤8：对步骤7得到的时域信道响应进行信道均衡，使用最大比合并mrc频域均衡器，得到频域接收信号y
k
经过均衡处理后得到的频域信号为z
k
＝c
k
y
k
；c
k
为第k个接收单载波
数据块的mrc均衡器矩阵；步骤9：对步骤8得到的频域信号z
k
进行解扩和解调，利用相关fft处理得到扩频码的索引位置以及循环移位位置利用循环移位位置对第k个接收单载波数据块进行相位解调；步骤10：对于n个接收单载波数据块均重复步骤7到步骤9，从而完成接收端的解调过程，对解调后的数据进行并串转换输出。2.根据权利要求1所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤2具体过程为：从n
c
个pecan扩频序列中随机选取其中一个扩频序列p
n
∈p,n＝0,1,
…
,n
c-1作为导频信号，剩余的n
c-1个扩频序列分成(j-1)个子集，每个子集包含m1个pecan扩频序列；再从每个子集中随机选择一个pecan扩频序列，共得到(j-1)个扩频序列。3.根据权利要求1所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤3具体过程为：对于pecan扩频序列p
j
，将待通信数据信息比特中的b2＝log2(m2)个信息比特映射到p
j
上，得到序列该序列表示扩频序列p
j
以b2为移位距离向右循环移位所得到的序列：再将待通信数据信息比特中的b3＝log2(m3)个psk调制符号映射到序列上，得到序列将步骤2得到的(j-1)个扩频序列对应的序列相叠加形成第k个单载波块对应的数据块f
k
4.根据权利要求1所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤7中时域信道估计的结果表示为：其中，γ是设定的阈值。5.根据权利要求1所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤8中，均衡器的抽头系数表示为：其中，c
t
为第k个接收单载波数据块的mrc均衡器矩阵c
k
中的第t个均衡器抽头系数。6.根据权利要求1所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤9中，扩频码的索引位置为：
为第k个接收单载波数据块对应的(j-1)个pecan扩频序列中的第j个pecan扩频序列。7.根据权利要求6所述一种叠加导频的m元多相序列扩频水声通信方法，其特征在于：步骤9中，循环移位位置为：其中表示根据扩频码的索引位置从第k个接收单载波数据块对应的(j-1)个pecan扩频序列中到的pecan扩频序列。

技术总结

发明提出一种叠加导频的M元多相序列扩频水声通信方法，首先运用具有理想脉冲自相关特性的多相PeCan序列作为扩频码序列，由于序列的产生具有随机性以及不可预测性，可以有效提高系统的抗干扰能力；其次在发射端，将用作导频信号的PeCan序列，与作为数据的PeCan序列相叠加，可以有效节省带宽；并构造循环前缀单载波块结构，避免多径衰落信道引起的符号间干扰(ISI)；同时使用最大比合并(Maximal RatioCombining,MRC)频域均衡器进行接收端信道均衡，可以有效对抗信道衰落；最后在接收端解扩解码时，只使用FFT/IFFT操作，可以有效降低计算复杂度。低计算复杂度。低计算复杂度。