CN107192982A

CN107192982A - 一种多用户水声超短基线定位系统及其信号形式设计方法

Info

Publication number: CN107192982A
Application number: CN201710356360.9A
Authority: CN
Inventors: 陈洲; 王熙赢
Original assignee: Jiangsu Hi-Target Ocean Information Technology Co Ltd
Current assignee: Jiangsu Hi-Target Ocean Information Technology Co Ltd
Priority date: 2017-05-19
Filing date: 2017-05-19
Publication date: 2017-09-22

Abstract

本发明公开了一种基于多用户水声超短基线定位系统及其信号形式设计方法；水声超短基线定位系统由安装在船体上的声头和水声应答器两部分组成。本发明利用迹函数在二元域上构造出一类基于GMW序列的低相关区序列，得到在低相关区内互相关函数边峰取值为‑1的低相关区序列集，当水声超短基线系统在同时定位多个水声应答器时，若将同步接收的时间点控制在码分多址的低相关区之内的话可以提升码分多址直接序列扩频系统的容量；进行一个由正余弦型符号函数组成的副载波调制：BOC调制，BOC调制将信号功率调制到载波频率两边的旁瓣上，提高了频带利用率和增加Gabor带宽，同时在临近一个码片范围内的相关函数主瓣更尖锐，进一步保证多用户信号检测和分离。

Description

一种多用户水声超短基线定位系统及其信号形式设计方法

技术领域

本发明属于水声工程领域，涉及一种基于多用户水声超短基线定位系统及其信号形式设计方法。

背景技术

水声超短基线系统是一种常见的水声定位技术，与其它基线相比，超短基线定位系统基线基阵尺寸小，易于安装，轻便灵活，操作方便等优点。传统的水声超短基线系统只支持单一用户或者多个窄带用户的同时定位跟踪，且当进行多用户同时定位跟踪时，用户之间的干扰是多用户数量受限的主要原因。在解决填海造地、水运交通时的铺排作业过程中，一般采用超短基线和长基线组合导航的技术，一方面指导铺排施工，另一方面对完成的排体进行高精度水下定位检测，此时需要对个水声应答器的高精度定位信息，因此水声超短基线系统要满足多个宽带用户的精确定位跟踪需求。

发明内容

本发明所要解决的技术问题是针对背景技术的不足提供了一种基于多用户水声超短基线定位系统及其信号形式设计方法。

本发明为解决上述技术问题采用以下技术方案

一种多用户水声超短基线定位系统，包含声头和水声应答器，所述水声应答器包含入水检测电路模块、收发合置换能器模块、MCU模块、DSP模块、发射电路模块、接收电路模块、压力传感器模块、电源模块；

入水检测电路模块，用于检测水声应答器是否入水，入水即给系统供电，离开水即断电；

收发合置换能器模块，用于接收声头发射的声信号，并发射定位信号给声头；

MCU模块，与DSP模块连接，用于检测声头发射的唤醒信号，以及控制DSP模块的工作；

DSP模块，接受MCU模块的控制，用于解算声头发射的定位信号；

发射电路模块，用于对应答声信号进行数模转换和功率放大后发送给收发合置换能器；

接收电路模块，用于对收发合置换能器接收的声信号进行功率放大、滤波处理和模数转换；

压力传感器模块，用于检测水声应答器的入水深度；

电源模块，用于负责整个水声应答器的电源供给；

所述声头包含换能器基阵模块、接收电路模块、发射电路模块、FPGA模块、ARM模块、甲板处理单元模块、声基阵位置和姿态校正模块；

换能器基阵模块，用于发射声信号给水声应答器，并接收水声应答器发射的定位信号；

接收电路模块，用于对水声应答器发射的定位信号进行前期预处理；

发射电路模块，用于对发射声信号进行前期预处理；

FPGA模块，用于接收电路模块之后的数据采集，并驱动DAC模块和同步声学数据、姿态数据、GPS数据；

ARM模块，用于与上位机通信；

甲板处理单元模块，用于处理声学数据，并发送控制命令；

声基阵位置和姿态校正模块，用于提供声头的精确位置和姿态信息；

一种基于多用户水声超短基线定位系统的信号形式设计方法，具体包含如下步骤：

步骤1，构造GMW序列；

步骤2，根据GMW序列构造低相关序列集；

步骤3，进行BOC调制。

作为本发明一种基于多用户水声超短基线定位系统的信号形式设计方法的进一步优选方案，GMW序列如下：

其中，α∈GF(pⁿ)且是本原元，是从扩域GF(pⁿ)映射到基域GF(p^m)的函数，I是不为1的陪集首，且0≤k＜2ⁿ-2，0≤I＜2ⁿ-1，gcd(I，2^m-1)＝1。

作为本发明一种基于多用户水声超短基线定位系统的信号形式设计方法的进一步优选方案，所述步骤2具体步骤如下：

同一本原产生的两个GMW序列，定义如下：

若m、n、r₁、r₂为整数，m|n,gcd(r₁,2ⁿ-1)＝gcd(r₂,2^m-1)＝1，T＝(2ⁿ-1)/(2^m-1)，α是二元域的本原元，则两个GMW序列定义为：

其中，x、y/{x_i}、{y_i}是两个GMW序列，r₁、r₂属于陪集，α∈GF(pⁿ)且是本原元。互相关函数定义为：

其中，R_x,y(τ)是互相关函数，R_x',y'(τ/T)是m序列的互相关函数，x',y'是m序列，T是序列长度，τ是序列相位。且有：

其中，x'、y'/{x_i'}、{y_i'}是m序列，T是序列长度，α∈GF(pⁿ)且是本原元，r₁、r₂属于陪集。

由同一本原元产生的GMW序列具有良好的互相关特性，满足低相关序列要求，但是由同一本原元生成的序列数量有限；

构造二元低相关序列集，具体步骤如下：

步骤2.1，同一本原元产生的两个GMW序列如下：

其中，a、b是GMW序列，m|n，gcd(r₁,2ⁿ-1)＝gcd(r₂,2^m-1)＝1，T＝(2ⁿ-1)/(2^m-1)，α是二元域的本原元。a'、b'是a、b分别对应着GF(2^m)上的m序列，即：

a'＝{a_i'}＝{Tr₁ ^m(α^iTr)}

b'＝{b_i'}＝{Tr₁ ^m(α^iTs)}

其中，a'、b'/{a_i'}、{b_i'}是m序列，T是序列长度，r、s属于陪集。

步骤2.2，计算m序列a'和b'的互相关函数，若有M-1个时延满足：

R_a',b'(l₁)＝R_a',b'(l₂)＝…R_a',b'(l_M-1)＝-1

其中，R_a',b'(l)是m序列a'和b'的互相关函数，l是m序列的不同相位。

步骤2.3，基于GMW序列a、b构成的复合序列集是：

其中，A是复合序列集，S是左移算子，a、b是两个GMW序列，T是序列长度，在序列集中任意两个序列的互相关函数满足：

其中，R_x,y(τ)是复合序列集其中任两个序列的互相关函数，R_x',y'(τ/T)是构成复合序列集的两个GMW序列的互相关函数，T是序列长度，τ是序列相位。此时，构造出了序列数目是M，低相关值等于-1，低相关区长度等于T的一个LCZ序列集A。

作为本发明一种基于多用户水声超短基线定位系统的信号形式设计方法的进一步优选方案，所述步骤3具体如下：

在BPSK调制基础上，进行一个由正余弦型符号函数组成的副载波调制就是BOC调制，BOC调制偶数阶信号的功率谱密度公式如下：

其中，f、f_s、f_c分别是载波频率、副载波频率和码元速率。

BPSK信号调制是基带脉冲控制载波相位的一种数字调制形式，它的功率谱密度公式：

其中，A是幅度，f、f_c分别是载波频率和码元速率。

本发明采用以上技术方案与现有技术相比，具有以下技术效果：

1、本发明利用迹函数在二元域上构造出一类基于GMW序列的低相关区序列，得到在低相关区内互相关函数边峰取值为-1的低相关区序列集，当水声超短基线系统在同时定位多个水声应答器时，若将同步接收的时间点控制在码分多址的低相关区之内的话可以提升码分多址直接序列扩频系统的容量。

2、本发明利用在BPSK调制基础上，进行一个由正余弦型符号函数组成的副载波调制：BOC调制，BOC调制将信号功率调制到载波频率两边的旁瓣上，提高了频带利用率和增加Gabor带宽，同时在临近一个码片范围内的相关函数主瓣更尖锐，进一步保证多用户信号检测和分离。

附图说明

图1为水声应答器硬件详细设计框图；

图2为声头硬件详细设计框图；

图3为低相关序列集中序列的周期自相关和周期互相关图；

图4为水声超短基线系统三个用户同时工作一个用户同步接收时低相关性能比较图；

图5为水声超短基线系统三个用户同时工作一个用户同步接收时低相关性能比较的局部放大图；

图6为BOC(1,1)和BPSK的基带功率谱密度曲线；

图7为给定初始条件的BOC(1,1)调制和BPSK调制的自相关图；

图8为基于低相关序列的BOC(1,1)和BPSK性能比较图；

图9为基于低相关序列的BOC(1,1)和BPSK性能比较的局部放大图；

图10为本发明的系统结构图。

图中标号具体如下：1-声头，2-水声应答器。

具体实施方式

下面结合附图对本发明的技术方案做进一步的详细说明：

一种多用户水声超短基线定位系统，如图10所示，包含声头和水声应答器，如图1所示，所述水声应答器包含入水检测电路、收发合置换能器、MCU模块、DSP模块、接收电路、压力传感器、电源模块；

如图2所示，所述声头包含换能器基阵模块、接收电路模块、发射电路模块、FPGA模块、ARM模块、甲板处理单元模块、声基阵位置和姿态校正模块:换能器基阵模块用于发射声信号给水声应答器，并接收水声应答器发射的定位信号；接收电路模块(前放、滤波、AGC、ADC)对水声应答器发射的定位信号进行前期预处理；FPGA模块用于接收电路模块之后的数据采集，并驱动DAC模块和同步声学数据、姿态数据、GPS数据；ARM模块用于与上位机通信；甲板处理单元模块用于处理声学数据，并发送控制命令；声基阵位置和姿态校正模块用于提供声头的精确位置和姿态信息。

低互相关序列集的构造。所述一种基于扩频通信的多用户水声超短基线定位系统信号形式设计算法利用GMW序列构造低相关序列集，在序列集内的任意两个序列具有一定长度的低互相关窗，当某个序列同步接收时可以有效减小其它序列的多址干扰。

迹函数实质是从扩域GF(pⁿ)映射到基域GF(p^m)的函数，定义为：

GMW序列，定义式为：

其中，α∈GF(pⁿ)且是本原元，0≤k＜2ⁿ-2，0≤I＜2ⁿ-1，gcd(I，2^m-1)＝1；

基于GMW序列的低相关区序列集的构造如下，同一本原产生的两个GMW序列，定义如下：

若m、n、r₁、r₂为整数，m|n,gcd(r₁,2ⁿ-1)＝gcd(r₂,2^m-1)＝1，T＝(2ⁿ-1)/(2^m-1)，α是二元域的本原元。则两个GMW序列定义为：

互相关函数定义为：

其中，x',y'是m序列，且有：

由同一本原元产生的GMW序列具有良好的互相关特性，满足低相关序列要求。但是由同一本原元生成的序列数量有限，下面介绍Long等人基于GMW序列复合形式，构造二元低相关序列集，步骤如下：

第1步同一本原元产生的两个GMW序列如下：

其中，m|n，gcd(r₁,2ⁿ-1)＝gcd(r₂,2^m-1)＝1，T＝(2ⁿ-1)/(2^m-1)，α是二元域的本原元。a'、b'是a、b分别对应着GF(2^m)上的m序列，即：

a'＝{a_i'}＝{Tr₁ ^m(α^iTr)}

b'＝{b_i'}＝{Tr₁ ^m(α^iTs)}

第2步计算m序列a'和b'的互相关函数，若有M-1个时延满足：

R_a',b'(l₁)＝R_a',b'(l₂)＝…R_a',b'(l_M-1)＝-1

第3步基于GMW序列a、b构成的复合序列集是：

其中，S是左移算子。在序列集中任意两个序列的互相关函数满足：

此时，构造出了序列数目是M，低相关值等于-1，低相关区长度等于T的一个LCZ序列集A。

低相关序列集中任一个序列拥有和m序列一样完美的周期自相关函数，任两个序列的周期互相关函数在零时延附近有一段长度是T＝(2ⁿ-1)/(2^m-1)，数值是-1/1023低相关窗。在水声超短基线多用户定位跟踪模式下，单用户对声头发射的信号进行同步接收时，其它用户的多址干扰决定了系统多用户的数量，如果把同步时间点放在低相关区域内将有效地减小多址的干扰。定位信号持续时间超过水声信道时延扩展的情况下，并不是所有移位相关结果都影响多径带来的干扰，对于超出水声信道时延扩展的相关结果，对性能没有影响，如果把多径的位置控制在低相关区域内将有效地减小多径的干扰。低相关序列集中序列的周期自相关和周期互相关图如图3，上子图是低相关序列集中的任一个序列的循环自相关函数，下子图是低相关序列集内任意两个序列的循环互相关函数；

当水声超短基线系统三个水声应答(三个用户)同时工作时，声头将低相关序列集中的三通道序列叠加发射，每个水声应答器对应ID的低相关序列同步接收，此时其它水声应答器的低相关序列对当前水声应答器来说是多址干扰，此时的同步时间点必然是在低互相关区域内，其它水声应答器对当前同步接收水声应答器的多址干扰较小。水声超短基线系统三个用户同时工作一个用户同步接收时低相关性能比较图如图4，上子图发射信号形式是低相关序列集中的三个序列三通道叠加，中子图发射信号形式是不具有优良互相关特性的三个低相关序列三通道叠加，下子图发射信号形式是一个低相关序列与Gold序列和m序列三通道叠加，且水声超短基线系统三个用户同时工作一个用户同步接收时低相关性能比较的局部放大图如图5；

BOC调制。二进制偏移载波键控(BOC)技术是Galileo导航系统的主要调制体制，是针对原有信号频带与GPS增发的信号频带共用的目的提出来的。传统的BPSK调制方式是将信号的能量集中在载频附近。BOC调制方式具有频谱分裂的特点，将信号的能量推向载频两边。这既能使两种调制方式并存提高频带利用率，又能增加Gabor带宽，同时在临近一个码片范围内的相关函数主瓣更尖锐，进一步保证多用户信号检测和分离。

在BPSK调制基础上，进行一个由正余弦型符号函数组成的副载波调制就是BOC调制。BOC(fs，fc)调制偶数阶信号的功率谱密度公式如下：

其中，f_s、f_c分别是副载波频率和码率与基准码元速率的比值。

BPSK信号调制，也称二进制相移键控，是基带脉冲控制载波相位的一种数字调制形式，它的功率谱密度公式：

其中，A是幅度，f_c是码元速率。当幅度A为1，基准码元速率为500cps，可得BOC(1,1)和BPSK的基带功率谱密度曲线如图6，由图可以看出，BOC(1,1)调制在频带宽度不变的条件下将频谱分裂成左右对称结构，f_s、f_c这两个参数决定了两边主瓣之间主、旁瓣的个数。因此推算BOC(1,1)主瓣之间主、旁瓣的个数为两个，即证明主瓣间旁瓣的个数为零；

由维纳-辛钦定理可知，在功率谱密度已知条件下可推导出自相关函数。假设是频率带限的接收系统，它的归一化自相关函数与功率谱密度的关系如下：

其中，B是水听器频带宽度，给出初始条件和功率谱密度就得到BOC(1,1)调制和BPSK调制的自相关函数如图7，由图可以看出，在相同的频带宽度时，BPSK调制方式的自相关曲线没有BOC(1,1)调制的自相关曲线尖锐，这对时延精度的提高大有帮助；

(3.2)基于BOC调制的低相关窗水声同步系统。下面分别对基于低相关序列的不同调制方式的同步性能进行比较研究，调制方式采用传统的BPSK调制方式和改进的BOC调制方式，同样采用水声超短基线系统中三个水声应答同时工作时，声头将低相关序列集中的三通道序列叠加发射，每个水声应答器对应ID的低相关序列同步接收，基于低相关序列的BOC(1,1)和BPSK性能比较如图8，基于低相关序列的BOC(1,1)和BPSK性能比较的局部放大图如图9，由图8和图9可以看出，在取得低相关区的基础上，基于BOC调制方式呈现明显的频谱分裂特性，且基于BOC调制方式在正负一个码片范围内的归一化相关峰最大边峰的值是0.6158，主峰宽度约为0.3ms。而基于BPSK调制方式在正负一个码片范围内的相关峰最大边峰的值是0.8747，主峰宽度约为1ms。可知次大峰的值得到明显的改善，且主峰宽度明显减小。这可以大大减小多径干扰的影响，提高在正负码片范围内的定位精度。

Claims

1.一种多用户水声超短基线定位系统，其特征在于：包含声头和水声应答器，所述水声应答器包含入水检测电路模块、收发合置换能器模块、MCU模块、DSP模块、发射电路模块、接收电路模块、压力传感器模块、电源模块；

压力传感器模块，用于检测水声应答器的入水深度；

电源模块，用于负责整个水声应答器的电源供给；

发射电路模块，用于对发射声信号进行前期预处理；

ARM模块，用于与上位机通信；

甲板处理单元模块，用于处理声学数据，并发送控制命令；

声基阵位置和姿态校正模块，用于提供声头的精确位置和姿态信息。

2.一种基于多用户水声超短基线定位系统的信号形式设计方法，其特征在于：具体包含如下步骤：

步骤1，构造GMW序列；

步骤2，根据GMW序列构造低相关序列集；

步骤3，进行BOC调制。

3.根据权利要求2所述的一种基于多用户水声超短基线定位系统的信号形式设计方法，

其特征在于：GMW序列如下：

<mrow> <mi>u</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Tr</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>(</mo> <msup> <mi>&alpha;</mi> <mi>k</mi> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mi>I</mi> </msup> <mo>&rsqb;</mo> </mrow>

4.根据权利要求2所述的一种基于多用户水声超短基线定位系统的信号形式设计方法，其特征在于：所述步骤2具体步骤如下：

同一本原产生的两个GMW序列，定义如下：

<mrow> <mi>x</mi> <mo>=</mo> <mo>{</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Tr</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>(</mo> <msup> <mi>&alpha;</mi> <mi>i</mi> </msup> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> </msup> <mo>&rsqb;</mo> <mo>}</mo> </mrow>

<mrow> <mi>y</mi> <mo>=</mo> <mo>{</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Tr</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>(</mo> <msup> <mi>&alpha;</mi> <mi>i</mi> </msup> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mn>2</mn> </msub> </msup> <mo>&rsqb;</mo> <mo>}</mo> </mrow>

<mrow> <msup> <mi>x</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>x</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mrow> <mo>(</mo> <msup> <mi>&alpha;</mi> <mrow> <msub> <mi>iTr</mi> <mn>1</mn> </msub> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>

<mrow> <msup> <mi>y</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>y</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mrow> <mo>(</mo> <msup> <mi>&alpha;</mi> <mrow> <msub> <mi>iTr</mi> <mn>1</mn> </msub> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>

其中，x'、y'/{x′_i}、{y′_i}是m序列，T是序列长度，α∈GF(pⁿ)且是本原元，r₁、r₂属于陪集。

构造二元低相关序列集，具体步骤如下：

步骤2.1，同一本原元产生的两个GMW序列如下：

<mrow> <mi>a</mi> <mo>=</mo> <mo>{</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Tr</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>(</mo> <msup> <mi>&alpha;</mi> <mi>i</mi> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mi>r</mi> </msup> <mo>&rsqb;</mo> <mo>}</mo> </mrow>

<mrow> <mi>b</mi> <mo>=</mo> <mo>{</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mo>&lsqb;</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Tr</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>(</mo> <msup> <mi>&alpha;</mi> <mi>i</mi> </msup> <mo>)</mo> <mo>)</mo> </mrow> <mi>s</mi> </msup> <mo>&rsqb;</mo> <mo>}</mo> </mrow>

<mrow> <msup> <mi>a</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>a</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mrow> <mo>(</mo> <msup> <mi>&alpha;</mi> <mrow> <mi>i</mi> <mi>T</mi> <mi>r</mi> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>

<mrow> <msup> <mi>b</mi> <mo>&prime;</mo> </msup> <mo>=</mo> <mo>{</mo> <msubsup> <mi>b</mi> <mi>i</mi> <mo>&prime;</mo> </msubsup> <mo>}</mo> <mo>=</mo> <mo>{</mo> <msubsup> <mi>Tr</mi> <mn>1</mn> <mi>m</mi> </msubsup> <mrow> <mo>(</mo> <msup> <mi>&alpha;</mi> <mrow> <mi>i</mi> <mi>T</mi> <mi>s</mi> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>

其中，a'、b'/{a′_i}、{b′_i}是m序列，T是序列长度，r、s属于陪集。

R_a',b'(l₁)＝R_a',b'(l₂)＝…R_a',b'(l_M-1)＝-1

步骤2.3，基于GMW序列a、b构成的复合序列集是：

5.根据权利要求2所述的一种基于多用户水声超短基线定位系统的信号形式设计方法法，其特征在于：所述步骤3具体如下：

<mrow> <msub> <mi>G</mi> <mrow> <mi>B</mi> <mi>O</mi> <mi>C</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>f</mi> <mi>c</mi> </msub> <mo>&CenterDot;</mo> <msup> <mrow> <mo>&lsqb;</mo> <mfrac> <mrow> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>&pi;</mi> <mi>f</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>f</mi> <mi>s</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>&pi;</mi> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>c</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&pi;</mi> <mi>f</mi> </mrow> </mfrac> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow>

其中，f、f_s、f_c分别是载波频率、副载波频率和码元速率。

<mrow> <msub> <mi>G</mi> <mrow> <mi>B</mi> <mi>P</mi> <mi>S</mi> <mi>K</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msup> <mi>A</mi> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>c</mi> </msub> </mrow> </mfrac> <mo>&CenterDot;</mo> <msup> <mrow> <mo>&lsqb;</mo> <mfrac> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&pi;</mi> <mi>f</mi> <mfrac> <mn>1</mn> <msub> <mi>f</mi> <mi>c</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&pi;</mi> <mi>f</mi> <mfrac> <mn>1</mn> <msub> <mi>f</mi> <mi>c</mi> </msub> </mfrac> </mrow> </mfrac> <mo>&rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow>

其中，A是幅度，f、f_c分别是载波频率和码元速率。