CN104753561B - Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication - Google Patents
Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication Download PDFInfo
- Publication number
- CN104753561B CN104753561B CN201310739788.3A CN201310739788A CN104753561B CN 104753561 B CN104753561 B CN 104753561B CN 201310739788 A CN201310739788 A CN 201310739788A CN 104753561 B CN104753561 B CN 104753561B
- Authority
- CN
- China
- Prior art keywords
- sequence
- spread spectrum
- pseudo
- tau
- sequences
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种抑制水声通信中多途干扰的直接序列扩频调制方法,该直接序列扩频调制方法将互补序列对作为直接序列扩频调制的伪随机序列,引入水声通信,利用其构造的零相关窗内理想相关函数,实现比普通的有限长随机序列具有更好的多途抑制能力,能够完全抑制水声扩频通信中的多途干扰;根据实际需要,可选择两种“互补序列”对的组合方法构造直接序列扩频调制的伪随机序列:采用时分正交组合方式能够降低发射占空比,从而起到冷却功率放大器的作用,这种信号采用正交解调时对载波相位不敏感,不需要精确的载波相位同步;采用载波相位正交组合方式的信息发射效率同普通的m序列、Gold序列等伪随机序列扩频调制一致,在信噪比较高时,该方式的解调误码率较好。
The invention provides a direct sequence spread spectrum modulation method for suppressing multi-channel interference in underwater acoustic communication. The direct sequence spread spectrum modulation method uses a complementary sequence pair as a pseudo-random sequence of direct sequence spread spectrum modulation, introduces underwater acoustic communication, and utilizes The ideal correlation function in the zero-correlation window constructed by it has better multi-path suppression ability than ordinary finite-length random sequences, and can completely suppress multi-path interference in underwater acoustic spread spectrum communication; according to actual needs, two kinds of The combination method of "complementary sequence" pairs constructs the pseudo-random sequence of direct sequence spread spectrum modulation: the time-division orthogonal combination method can reduce the transmission duty cycle, thereby cooling the power amplifier. It is not sensitive to the carrier phase and does not require precise carrier phase synchronization; the information transmission efficiency of the carrier phase orthogonal combination method is consistent with that of ordinary m-sequence, Gold sequence and other pseudo-random sequence spread spectrum modulation. When the signal-to-noise ratio is high, The demodulation error rate of this mode is better.
Description
技术领域technical field
本发明属于水声扩频通信技术领域,尤其涉及一种抑制水声通信中多途干扰的直接序列扩频调制方法。The invention belongs to the technical field of underwater acoustic spread spectrum communication, in particular to a direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication.
背景技术Background technique
在水声通信中,多途效应引起的符号间干扰(ISI)是限制信息传输速率的主要原因。对此,通过利用直接序列扩频调制不仅能提高信号的传输距离,还可根据PRN码尖锐的自相关特性,对延时大于一个扩频码片的多途信号有很强的抑制作用。因此,基于DSSS调制体制的水声通信受到了广泛的研究。在DSSS调制技术中,常用的扩频码有m序列、Gold码、混沌码等。但是在工程使用中,只能将有限长度的扩频码用于信号的调制和解调。然而,由文献“Welch L R.Lower bounds on the maximum cross-correlation of signals[J].IEEETransactions Information Theory,1974,20:397-399.”中给出的Welch界可知,对于长度有限的随机序列,不能获得理想的自相关或者互相关函数。因此,非理想的自相关和互相关使得扩频系统不能完全抑制符号间干扰(ISI)和多址干扰(MAI),改善随机序列的相关特性,成为抑制多途干扰的有效手段之一。In underwater acoustic communication, the inter-symbol interference (ISI) caused by the multipath effect is the main reason for limiting the information transmission rate. In this regard, the use of direct sequence spread spectrum modulation can not only increase the transmission distance of the signal, but also have a strong inhibitory effect on multi-path signals with a delay greater than one spread spectrum chip according to the sharp autocorrelation characteristics of the PRN code. Therefore, underwater acoustic communication based on DSSS modulation scheme has been extensively studied. In DSSS modulation technology, commonly used spread spectrum codes include m-sequence, Gold code, chaos code and so on. However, in engineering applications, only limited-length spreading codes can be used for signal modulation and demodulation. However, from the Welch bound given in the document "Welch L R. Lower bounds on the maximum cross-correlation of signals [J]. IEEETransactions Information Theory, 1974, 20:397-399." , the ideal autocorrelation or cross-correlation function cannot be obtained. Therefore, non-ideal auto-correlation and cross-correlation make the spread spectrum system unable to completely suppress inter-symbol interference (ISI) and multiple access interference (MAI), improving the correlation characteristics of random sequences has become one of the effective means to suppress multi-path interference.
为了减小甚至消除符号间干扰(ISI)和多址干扰(MAI),在文献“Sta′nczakS.Boche H.Haardt M.Are LAS-codes a miracle?[J].Global TelecommunicationsConference,2001.GLOBECOM'01.IEEE,2001,1:589-593.”中给出了用于LAS-CDMA系统的LS码。该LS码是一种具有零相关窗(IFW)的“互补序列”(Complementary Series),在零相关窗内,“互补序列”对具有理想的自相关和互相关函数。这种扩频码的提出,主要应用于移动无线通信网络,用以减小CDMA系统的符号间干扰(ISI)和多址干扰(MAI),提高蜂窝网络的容量。In order to reduce or even eliminate intersymbol interference (ISI) and multiple access interference (MAI), in the document "Sta′nczak S. Boche H. Haardt M. Are LAS-codes a miracle? [J]. Global Telecommunications Conference, 2001. GLOBECOM' 01. IEEE, 2001, 1:589-593." provides the LS codes for the LAS-CDMA system. The LS code is a "complementary series" with a zero-correlation window (IFW). Within the zero-correlation window, the pair of "complementary sequences" has ideal autocorrelation and cross-correlation functions. The proposal of this spreading code is mainly used in mobile wireless communication networks to reduce the intersymbol interference (ISI) and multiple access interference (MAI) of the CDMA system and improve the capacity of the cellular network.
发明内容Contents of the invention
本发明的目的在于,为解决现有的直接序列扩频调制方法不能完全抑制水声通信中多途干扰的技术问题,提供一种抑制水声通信中多途干扰的直接序列扩频调制方法,该直接序列扩频调制方法将互补序列对作为直接序列扩频调制的伪随机序列,引入水声通信,利用其构造的零相关窗内理想相关函数,实现比普通的有限长随机序列具有更好的多途抑制能力。The purpose of the present invention is to provide a direct sequence spread spectrum modulation method for suppressing multi-path interference in underwater acoustic communication in order to solve the technical problem that the existing direct sequence spread spectrum modulation method cannot completely suppress multi-path interference in underwater acoustic communication, The direct sequence spread spectrum modulation method uses the complementary sequence pair as the pseudo-random sequence of the direct sequence spread spectrum modulation, introduces underwater acoustic communication, and uses the ideal correlation function in the zero correlation window constructed by it to achieve better performance than ordinary finite-length random sequences. multipath suppression capability.
为了实现上述目的,本发明提供一种抑制水声通信中多途干扰的直接序列扩频调制方法,该直接序列扩频调制方法包括:In order to achieve the above object, the present invention provides a direct sequence spread spectrum modulation method for suppressing multi-channel interference in underwater acoustic communication, the direct sequence spread spectrum modulation method comprising:
步骤1)构造两个满足“互补序列”对关系的序列;Step 1) Construct two sequences that satisfy the "complementary sequence" pair relationship;
步骤2)将步骤1)中获得的两个序列进行正交组合生成伪随机序列,以该伪随机序列作为扩频码;Step 2) Perform orthogonal combination of the two sequences obtained in step 1) to generate a pseudo-random sequence, and use the pseudo-random sequence as a spreading code;
步骤3)根据扩频信号的传输带宽,计算出单个扩频码码片的持续时间,根据信息传输速率,计算出单个信息码片的持续时间,然后计算出单个信息码片所对应的扩频码的码片个数;Step 3) Calculate the duration of a single spread code chip according to the transmission bandwidth of the spread spectrum signal, calculate the duration of a single information chip according to the information transmission rate, and then calculate the spread spectrum corresponding to a single information chip The number of code chips;
步骤4)根据步骤3)中计算的结果,将发射的信息码片按照对应个数关系分别与步骤2)中得到的伪随机序列做模2相乘运算后,组合生成复基带扩频信号;Step 4) According to the result calculated in step 3), the transmitted information chips are respectively multiplied by modulo 2 with the pseudo-random sequence obtained in step 2) according to the corresponding number relationship, and combined to generate a complex baseband spread spectrum signal;
步骤5)将步骤4)中得到的复基带扩频信号和信号载波相乘并取实部,得到扩频发射信号。Step 5) Multiply the complex baseband spread spectrum signal obtained in step 4) by the signal carrier and take the real part to obtain the spread spectrum transmission signal.
作为上述技术方案的进一步改进,步骤1)中生成的每个序列在连续时间域上表示为:As a further improvement of the above technical solution, each sequence generated in step 1) is represented in the continuous time domain as:
其中,序列的长度为N,t表示时间,cn表示第n个伪随机序列的序列值,rect(·)为方波函数TC表示伪随机序列单个码片的持续时间;Among them, the length of the sequence is N, t represents time, c n represents the sequence value of the nth pseudo-random sequence, rect( ) is a square wave function T C represents the duration of a single chip of the pseudo-random sequence;
所述的“互补序列”对的相关函数定义为:The correlation function of the "complementary sequence" pair is defined as:
其中,τ表示序列Ai和Aj的延时,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。Among them, τ represents the delay of sequences A i and A j , A i , A j , B i , B j are all pseudo-random sequences satisfying the relationship of "complementary sequence". When i=j, the above formula describes the autocorrelation characteristic ; i≠j describes the cross-correlation characteristics.
作为上述技术方案的进一步改进,所述步骤4)中发射的信息码片表示为:As a further improvement of the above technical solution, the information chips transmitted in step 4) are expressed as:
其中,dn表示第n个信息序列值,t表示时间,Td表示单个信息序列持续的时间。Among them, d n represents the value of the nth information sequence, t represents time, and T d represents the duration of a single information sequence.
作为上述技术方案的进一步改进,所述步骤5)中的信号载波表示为:As a further improvement of the above technical solution, the signal carrier in step 5) is expressed as:
其中,表示虚数单位,ω为载波角频率,为载波初始相位,t表示时间。in, Indicates the imaginary unit, ω is the carrier angular frequency, is the initial phase of the carrier, and t represents the time.
作为上述技术方案的进一步改进,所述步骤2)中正交组合生成伪随机序列采用时分方式,将两个序列按时间的先后顺序排列后,在两个序列之间加入一定长度的值为0的序列,得到的伪随机序列表示为:As a further improvement of the above technical solution, in the step 2) the orthogonal combination generates a pseudo-random sequence in a time-division manner. After the two sequences are arranged in chronological order, a certain length of value 0 is added between the two sequences sequence, the obtained pseudo-random sequence is expressed as:
其中,An和Bn表示序列,W表示在两个序列之间加入的0的个数,S表示伪随机序列,其长度为2(N+W);Among them, A n and B n represent sequences, W represents the number of 0s added between the two sequences, S represents a pseudo-random sequence, and its length is 2(N+W);
根据相关函数的定义,公式(5)的相关函数表示为:According to the definition of the correlation function, the correlation function of formula (5) is expressed as:
其中,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。Among them, A i , A j , B i , and B j are all pseudo-random sequences satisfying the "complementary sequence" relationship. When i=j, the above formula describes the autocorrelation characteristic; i≠j describes the cross-correlation characteristic.
作为上述技术方案的进一步改进,所述步骤5)中得到扩频发射信号表示为:As a further improvement of the above technical solution, the spread spectrum transmission signal obtained in step 5) is expressed as:
其中,TC表示伪随机序列单个码片的持续时间,t表示时间,表示信息码片,和表示序列,ω为载波角频率为载波初始相位,Td表示单个信息序列持续的时间。Among them, T C represents the duration of a single chip of the pseudo-random sequence, t represents the time, Indicates the information chip, with Indicates the sequence, ω is the carrier angular frequency is the initial phase of the carrier, and T d represents the duration of a single information sequence.
作为上述技术方案的进一步改进,所述步骤2)中正交组合生成伪随机序列采用载波相位正交方式,将两个序列分别调至在相位相差为90度的正交载波上,得到的伪随机序列表述为:As a further improvement of the above technical solution, in the step 2) the orthogonal combination generates the pseudo-random sequence using the carrier phase quadrature method, and the two sequences are respectively adjusted to the orthogonal carrier with a phase difference of 90 degrees, and the pseudo-random sequence obtained is A random sequence is expressed as:
CS=An+jBn (8)CS=A n + jB n (8)
其中,An和Bn表示序列,表示虚数单位,CS表示组合后的复数形式的伪随机序列,其长度为N;Among them, A n and B n represent sequences, Represents the imaginary number unit, CS represents the pseudo-random sequence in complex number form after combination, and its length is N;
根据相关函数的定义,公式(8)的相关函数为:According to the definition of correlation function, the correlation function of formula (8) is:
所述相关函数的实部表示为:The real part of the correlation function Expressed as:
所述相关函数的虚部表示为:The imaginary part of the correlation function Expressed as:
其中,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。Among them, A i , A j , B i , and B j are all pseudo-random sequences satisfying the "complementary sequence" relationship. When i=j, the above formula describes the autocorrelation characteristic; i≠j describes the cross-correlation characteristic.
作为上述技术方案的进一步改进,所述步骤5)中得到扩频发射信号表示为:As a further improvement of the above technical solution, the spread spectrum transmission signal obtained in step 5) is expressed as:
其中,TC表示伪随机序列单个码片的持续时间,t表示时间,表示信息码片,和表示序列,表示组合后的复数形式的伪随机序列,ω为载波角频率,为载波初始相位,Td表示单个信息序列持续的时间。Among them, T C represents the duration of a single chip of the pseudo-random sequence, t represents the time, Indicates the information chip, with represents a sequence, Indicates the combined pseudo-random sequence in the form of complex numbers, ω is the carrier angular frequency, is the initial phase of the carrier, and T d represents the duration of a single information sequence.
作为上述技术方案的进一步改进,所述的“互补序列”对采用LS码。As a further improvement of the above technical solution, the "complementary sequence" pair adopts LS code.
本发明的一种抑制水声通信中多途干扰的直接序列扩频调制方法的优点在于:The advantages of a direct sequence spread spectrum modulation method for suppressing multi-path interference in underwater acoustic communication of the present invention are:
本发明的采用“互补序列”对作为水声扩频通信中的伪随机序列,能构造出理想的自相关和互相关函数,在理论上能够完全抑制水声扩频通信中的多途干扰;提供两种“互补序列”对的组合方法,使用时可根据具体情况选择:1)时分正交组合方式:互补序列之间添加的值为0的序列可以降低发射占空比,从而起到冷却功率放大器的作用;这种信号采用正交解调时对载波相位不敏感,不需要精确的载波相位同步;在信噪比较低时,该方法的解调误码率优于方式2);2)载波相位正交组合方式:该方法的信息发射效率同普通的m序列、Gold序列等伪随机序列扩频调制一致,比方式1)效率更高;在信噪比较高时,该方法的解调误码率优于方式1)。The present invention adopts the "complementary sequence" pair as a pseudo-random sequence in underwater acoustic spread spectrum communication, which can construct ideal autocorrelation and cross-correlation functions, and can completely suppress multi-path interference in underwater acoustic spread spectrum communication in theory; Provide two combination methods of "complementary sequence" pairs, which can be selected according to the specific situation: 1) Time-division orthogonal combination method: the sequence with a value of 0 added between complementary sequences can reduce the emission duty cycle, thereby cooling The role of the power amplifier; this kind of signal is not sensitive to the carrier phase when using quadrature demodulation, and does not require precise carrier phase synchronization; when the signal-to-noise ratio is low, the demodulation error rate of this method is better than that of mode 2); 2) Carrier phase quadrature combination method: The information transmission efficiency of this method is consistent with that of ordinary m-sequence, Gold sequence and other pseudo-random sequence spread spectrum modulation, which is more efficient than method 1); when the signal-to-noise ratio is high, this method The demodulation bit error rate is better than mode 1).
附图说明Description of drawings
图1是本发明的一种抑制水声通信中多途干扰的直接序列扩频调制方法的流程图。Fig. 1 is a flowchart of a direct sequence spread spectrum modulation method for suppressing multi-path interference in underwater acoustic communication according to the present invention.
图2是本发明中“互补序列”对采用时分方式正交组合生成伪随机序列的结构示意图。Fig. 2 is a schematic diagram of the structure of the "complementary sequence" pair in the present invention, which is orthogonally combined to generate a pseudo-random sequence in a time-division manner.
图3是本发明中采用时分方式正交组合生成伪随机序列的自相关函数。Fig. 3 is the autocorrelation function of pseudo-random sequence generated by orthogonal combination in time-division mode in the present invention.
图4是本发明中采用时分方式正交组合生成伪随机序列的互相关函数。Fig. 4 is the cross-correlation function of pseudo-random sequence generated by orthogonal combination in time-division mode in the present invention.
图5是本发明中采用载波相位正交方式正交组合生成伪随机序列的自相关函数。Fig. 5 is the autocorrelation function of the pseudo-random sequence generated by the quadrature combination of carrier phase quadrature in the present invention.
图6是本发明中采用载波相位正交方式正交组合生成伪随机序列的互相关函数。Fig. 6 is the cross-correlation function of the pseudo-random sequence generated by the quadrature combination of carrier phase quadrature in the present invention.
图7是本发明中基于时分方式得到的扩频发射信号的湖试CZT时频搜索三维图。FIG. 7 is a three-dimensional diagram of the CZT time-frequency search of the lake test CZT time-frequency search based on the time-division method in the present invention.
图8是本发明中基于载波相位正交方式得到的扩频发射信号的湖试CZT时频搜索三维图。FIG. 8 is a three-dimensional diagram of the CZT time-frequency search of the lake test CZT time-frequency search of the spread spectrum transmission signal obtained based on the carrier phase quadrature method in the present invention.
具体实施方式detailed description
下面结合附图和实施例对本发明所述一种抑制水声通信中多途干扰的直接序列扩频调制方法进行详细说明。A direct sequence spread spectrum modulation method for suppressing multi-path interference in underwater acoustic communication according to the present invention will be described in detail below in conjunction with the accompanying drawings and embodiments.
如图1所示,本发明的一种抑制水声通信中多途干扰的直接序列扩频调制方法,所述直接序列扩频调制方法包括:As shown in Fig. 1, a kind of direct sequence spread spectrum modulation method of the present invention suppresses multi-path interference in underwater acoustic communication, described direct sequence spread spectrum modulation method comprises:
步骤1)构造两个满足“互补序列”对关系的序列A码和B码,即两个序列互为正交序列,且A码的非周期自相关或者互相关函数与B码的自相关或者互相关函数在零延时处相等,在其他地方绝对值相等,符号相反;Step 1) Construct two sequences A code and B code that satisfy the "complementary sequence" pair relationship, that is, the two sequences are mutually orthogonal sequences, and the aperiodic autocorrelation or cross-correlation function of A code and the autocorrelation or The cross-correlation functions are equal at zero delay, and are equal in absolute value and opposite in sign elsewhere;
步骤2)将步骤1)中获得的两个序列进行正交组合生成伪随机序列,以该伪随机序列作为扩频码;Step 2) Perform orthogonal combination of the two sequences obtained in step 1) to generate a pseudo-random sequence, and use the pseudo-random sequence as a spreading code;
步骤3)根据扩频信号的传输带宽,计算出单个扩频码码片的持续时间,根据信息传输速率,计算出单个信息码片的持续时间,然后计算出单个信息码片所对应的扩频码的码片个数;Step 3) Calculate the duration of a single spread code chip according to the transmission bandwidth of the spread spectrum signal, calculate the duration of a single information chip according to the information transmission rate, and then calculate the spread spectrum corresponding to a single information chip The number of code chips;
步骤4)根据步骤3)中计算的结果,将发射的信息码片按照对应个数关系分别与步骤2)中得到的伪随机序列做模2相乘运算后,组合生成复基带扩频信号;Step 4) According to the result calculated in step 3), the transmitted information chips are respectively multiplied by modulo 2 with the pseudo-random sequence obtained in step 2) according to the corresponding number relationship, and combined to generate a complex baseband spread spectrum signal;
步骤5)将步骤4)中得到的复基带扩频信号和信号载波相乘并取实部,得到扩频发射信号。Step 5) Multiply the complex baseband spread spectrum signal obtained in step 4) by the signal carrier and take the real part to obtain the spread spectrum transmission signal.
基于上述直接序列扩频调制方法,假设序列的长度都为N,则满足“互补序列”对关系的两个序列在连续时间域上可以表示为:Based on the above direct sequence spread spectrum modulation method, assuming that the length of the sequence is N, the two sequences satisfying the "complementary sequence" pair relationship can be expressed in the continuous time domain as:
其中,t表示时间,cn表示第n个伪随机序列的序列值,rect(·)为方波函数TC表示伪随机序列单个码片的持续时间。Among them, t represents time, c n represents the sequence value of the nth pseudo-random sequence, rect( ) is a square wave function T C represents the duration of a single chip of the pseudo-random sequence.
所述的“互补序列”对的相关函数可定义为:The correlation function of the "complementary sequence" pair can be defined as:
其中,τ表示序列Ai和Aj的延时,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。由公式(2)可以看出,“互补序列”对构造出的相关函数的自相关和互相关特性都是理想的,即自相关在0延时处有个脉冲;在0延时外,自相关和互相关都为0。Among them, τ represents the delay of sequences A i and A j , A i , A j , B i , B j are all pseudo-random sequences satisfying the relationship of "complementary sequence". When i=j, the above formula describes the autocorrelation characteristic ; i≠j describes the cross-correlation characteristics. It can be seen from the formula (2) that the autocorrelation and cross-correlation characteristics of the constructed correlation function of the "complementary sequence" are ideal, that is, the autocorrelation has a pulse at 0 delay; Both correlation and cross-correlation are 0.
另外,所述步骤4)中发射的信息码片可表示为:In addition, the information chips transmitted in step 4) can be expressed as:
其中,dn表示第n个信息序列值,t表示时间,Td表示单个信息序列持续的时间。Among them, d n represents the value of the nth information sequence, t represents time, and T d represents the duration of a single information sequence.
所述步骤5)中的信号载波可表示为:The signal carrier in the step 5) can be expressed as:
其中,表示虚数单位,ω为载波角频率,为载波初始相位,t表示时间。在本发明中,可作为伪随机序列的“互补序列”对及其构造方法有很多。在本实施例中,以LS码为例,作为水声扩频通信的“互补序列”对扩频码。假设(cN,sN)和(c′N,s′N)是同组互补序列对,当N=2时,基本的互补序列对码组可表示为:in, Indicates the imaginary unit, ω is the carrier angular frequency, is the initial phase of the carrier, and t represents the time. In the present invention, there are many pairs of "complementary sequences" that can be used as pseudo-random sequences and their construction methods. In this embodiment, the LS code is taken as an example as a "complementary sequence" pair spreading code for underwater acoustic spread spectrum communication. Assuming (c N ,s N ) and (c′ N ,s′ N ) are the same complementary sequence pair, when N=2, the basic complementary sequence pair code group can be expressed as:
按照如下的树形迭代,每一次迭代,可以获得长度增加一倍的互补序列对码组。According to the following tree iteration, each iteration can obtain a complementary sequence pair code group whose length is doubled.
经过n次扩展,可以获得2n对编码长度为2nN的“互补序列”对。After n times of expansion, 2 n pairs of "complementary sequence" pairs with a coding length of 2 n N can be obtained.
由于C码的非周期自相关函数与S码的自相关函数在零延时处相等,在其他地方绝对值相等,符号相反。因此,将C码和S码通过某种方式组合起来,在C码和S码的正交区间内,相关函数叠加,形成“零相关窗”,由此形成LS码。Since the non-periodic autocorrelation function of C code and the autocorrelation function of S code are equal at zero delay, the absolute value is equal in other places, and the sign is opposite. Therefore, the C code and the S code are combined in a certain way, and in the orthogonal interval of the C code and the S code, the correlation functions are superimposed to form a "zero correlation window", thereby forming the LS code.
为了将“互补序列”对中的序列A码和B码(A码和B码的序列值只能取+1或者-1)在零相关窗内构造理想相关函数,需要保证序列A码和B码不能出现相关运算,即序列A码和B码在某种意义上达到正交。为了在调制时达到正交的目的,本发明给出了两种可行的方法:In order to construct the ideal correlation function of the sequence A code and B code (the sequence value of A code and B code can only take +1 or -1) in the "complementary sequence" pair in the zero correlation window, it is necessary to ensure that the sequence A code and B code Correlation operations cannot occur, that is, sequence A code and B code are orthogonal in a certain sense. In order to achieve the purpose of orthogonality during modulation, the present invention provides two feasible methods:
1)时分方式:将“互补序列”对中的序列A码和B码按时间的先后顺序排列,在序列A码和B码之间加入一定长度的值为0的序列,由此实现两个序列A码和B码的正交组合生成伪随机序列。如图1所示的伪随机序列的结构可表示为:1) Time-division method: Arrange the sequence A code and B code in the "complementary sequence" pair in chronological order, and add a sequence with a certain length of 0 between the sequence A code and B code, thereby realizing two The orthogonal combination of sequence A code and B code generates a pseudo-random sequence. The structure of the pseudo-random sequence shown in Figure 1 can be expressed as:
其中,An和Bn表示序列,W表示在两个序列之间加入的0的个数,S表示伪随机序列,其长度为2(N+W);Among them, A n and B n represent sequences, W represents the number of 0s added between the two sequences, S represents a pseudo-random sequence, and its length is 2(N+W);
根据相关函数的定义,公式(7)的相关函数表示为:According to the definition of the correlation function, the correlation function of formula (7) is expressed as:
其中,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。Among them, A i , A j , B i , and B j are all pseudo-random sequences satisfying the "complementary sequence" relationship. When i=j, the above formula describes the autocorrelation characteristic; i≠j describes the cross-correlation characteristic.
显然,当延时|τ|<W时,采用时分方式正交组合的“互补序列”对形成的伪随机序列都能得到理想的相关函数。Obviously, when the time delay |τ|<W, the ideal correlation function can be obtained for the pseudo-random sequence formed by the "complementary sequence" pair orthogonally combined in a time-division manner.
2)载波相位正交方式:将“互补序列”对中的序列A码和B码分别调至在相位相差为90度的两个正交载波上,由此实现两个序列的正交组合生成伪随机序列。得到的伪随机序列可表示为:2) Carrier Phase Orthogonal Mode: The sequence A code and B code in the "complementary sequence" pair are respectively adjusted to two orthogonal carriers with a phase difference of 90 degrees, thereby realizing the orthogonal combination generation of the two sequences pseudorandom sequence. The obtained pseudo-random sequence can be expressed as:
CS=An+jBn (9)CS=A n + jB n (9)
其中,An和Bn表示序列,表示虚数单位,CS表示组合后的复数形式的伪随机序列,其长度为N;Among them, A n and B n represent sequences, Represents the imaginary number unit, CS represents the pseudo-random sequence in complex number form after combination, and its length is N;
根据相关函数的定义,公式(9)的相关函数表示为:According to the definition of the correlation function, the correlation function of formula (9) is expressed as:
相关函数的实部表示为:real part of the correlation function Expressed as:
相关函数的虚部表示为:imaginary part of the correlation function Expressed as:
其中,Ai,Aj,Bi,Bj均为满足“互补序列”关系的伪随机序列,当i=j时,上式描述自相关特性;i≠j则描述互相关特性。Among them, A i , A j , B i , and B j are all pseudo-random sequences satisfying the "complementary sequence" relationship. When i=j, the above formula describes the autocorrelation characteristic; i≠j describes the cross-correlation characteristic.
显然,从公式(11)的实部相关函数可以看出,采用载波相位正交方式正交组合的“互补序列”对形成的伪随机序列,其相关函数的实部,在任意延时上有理想的自相关和互相关函数。Obviously, from the correlation function of the real part of formula (11), it can be seen that the real part of the correlation function of the pseudo-random sequence formed by the "complementary sequence" pair orthogonally combined in the carrier phase quadrature method has Ideal autocorrelation and cross-correlation functions.
从公式(12)中可以看出,采用载波相位正交方式正交组合的扩频码的相关函数的虚部表示“互补序列”对中序列A码和B码的互相关作用的结果,虽然它的值在任意延时区域内并不恒等于0,但由于序列A码和B码正交,故虚部的幅度很小,对整个相关函数的影响也很小。为了充分利用组合码的理想相关特性,对于这种扩频码的相关解调,只需要取出相关结果的实部作为有效输出结果。It can be seen from the formula (12) that the imaginary part of the correlation function of the spread spectrum codes orthogonally combined in the carrier phase quadrature mode represents the result of the cross-correlation between the sequence A code and B code in the "complementary sequence", although Its value is not always equal to 0 in any delay area, but because the sequence A code and B code are orthogonal, the magnitude of the imaginary part is very small, and the impact on the entire correlation function is also small. In order to make full use of the ideal correlation characteristics of the combination code, for the correlation demodulation of this spreading code, only the real part of the correlation result needs to be taken out as an effective output result.
根据扩频信号的直接序列扩频调制方式,当采用时分方式时,结合公式(1)、(3)和(7)得到的扩频发射信号表示为:According to the direct sequence spread spectrum modulation method of the spread spectrum signal, when the time division method is adopted, the spread spectrum transmitted signal obtained by combining the formulas (1), (3) and (7) is expressed as:
当采用载波相位正交方式时,结合公式(1)、(3)和(9)得到的扩频发射信号表示为:When the carrier phase quadrature method is adopted, the spread spectrum transmission signal obtained by combining formulas (1), (3) and (9) is expressed as:
根据扩频通信的基本理论可知,假设1个信息码片内有至少一个完整的伪随机扩频序列,则由公式(13)和公式(14)计算得到的扩频发射信号在基带上的相关特性,分别由公式(3)和公式(4)的相关函数确定。According to the basic theory of spread spectrum communication, assuming that there is at least one complete pseudo-random spread spectrum sequence in one information chip, the correlation of the spread spectrum transmitted signal on the baseband calculated by formula (13) and formula (14) characteristics, which are determined by the correlation functions of Equation (3) and Equation (4), respectively.
下面将结合某次湖试的实际发射信号,对本发明的具体实施方式做进一步的详细描述。The specific embodiment of the present invention will be further described in detail below in combination with the actual transmission signal of a certain lake test.
实验发射的基于LS码的扩频发射信号参数如下:The parameters of the spread spectrum transmitted signal based on the LS code transmitted in the experiment are as follows:
信号带宽:4KHz~8KHz,中心频率:6KHz。Signal bandwidth: 4KHz ~ 8KHz, center frequency: 6KHz.
步骤201):生成“互补序列”对的序列A码和B码,该步骤可以按照公式(5)和(6)的迭代方式,生成所需要长度的“互补序列”对C码和S码。Step 201): Generate the sequence A code and B code of the "complementary sequence" pair. This step can generate the required length of the "complementary sequence" pair C code and S code according to the iterative method of formulas (5) and (6).
步骤202):根据步骤201)中的“互补序列”对C码和S码进行正交组合生成伪随机序列,获得两种伪随机扩频序列:1)生成512个码片长度的C码和S码,在C码和S码间插入512个码片的0,按照时分正交形式组合成长度为2048个码片的扩频码;2)生成1024个码片长度的C码和S码,按照载波相位正交形式组合成长度为1024个码片的复扩频码。Step 202): According to the "complementary sequence" in step 201), the C code and the S code are orthogonally combined to generate a pseudo-random sequence, and two pseudo-random spread spectrum sequences are obtained: 1) Generate a C code with a length of 512 chips and S code, 512 chips of 0 are inserted between the C code and the S code, and combined into a spreading code with a length of 2048 chips according to the time-division orthogonal form; 2) Generate a C code and an S code with a length of 1024 chips , combined into a complex spreading code with a length of 1024 chips according to the carrier phase orthogonal form.
步骤203):根据扩频信号的传输带宽,计算出单个扩频码码片的持续时间TC;根据信息传输速率,计算出单个信息码片的持续时间Td;然后计算出单个信息码片对应的扩频序列的码片个数。Step 203): Calculate the duration T C of a single spread spectrum code chip according to the transmission bandwidth of the spread spectrum signal; calculate the duration T d of a single information chip according to the information transmission rate; then calculate a single information chip The number of chips of the corresponding spreading sequence.
步骤204):按照步骤203)计算的结果,将发射的信息码片按照对应个数关系分别与步骤202)中得到的两种伪随机扩频序列做模2相乘运算,并根据公式(7)和(9)给出的组合关系,组合生成复基带扩频信号。Step 204): According to the result calculated in step 203), the transmitted information chips are multiplied by modulo 2 with the two pseudo-random spread spectrum sequences obtained in step 202) according to the corresponding number relationship, and according to the formula (7 ) and (9) to combine to generate a complex baseband spread spectrum signal.
步骤205:根据公式(13)或者(14)进行直接序列扩频调制,将复基带扩频信号和信号载波相乘并取实部,得到扩频发射信号。Step 205: Perform direct sequence spread spectrum modulation according to formula (13) or (14), multiply the complex baseband spread spectrum signal by the signal carrier and take the real part to obtain the spread spectrum transmit signal.
如图3和图4示出了根据上述湖试实验,采用时分方式正交组合生成伪随机序列的自相关和互相关结果,它由长度为512个码片的C码和S码按照时分正交的形式组合,在C码和S码之间插入了512个码片长度的0。从图中可知,形成的零相关窗(IFW)的区间为(-512,512),在零相关窗内,组合码具有理想的自相关和互相关函数。Figure 3 and Figure 4 show the autocorrelation and cross-correlation results of the pseudo-random sequence generated by the time-division orthogonal combination according to the above-mentioned lake test experiment, which consists of the C code and the S code with a length of 512 chips according to the time-division regularization Combination in the form of intersection, 512 chip-length 0s are inserted between the C code and the S code. It can be seen from the figure that the interval of the formed zero-correlation window (IFW) is (-512,512), and within the zero-correlation window, the combined code has ideal autocorrelation and cross-correlation functions.
如图5和图6示出了根据上述湖试实验,采用载波相位正交方式正交组合生成伪随机序列的自相关和互相关结果,它由长度为1024个码片的C码和S码按照复数实部和虚部的形式组合。从图中可知,在任意延时内,这种组合码相关函数的实部具有理想的自相关和互相关值,其与公式(11)所描述的理论结果一致。Figure 5 and Figure 6 show the autocorrelation and cross-correlation results of the pseudo-random sequence generated by the orthogonal combination of the carrier phase orthogonal method according to the above-mentioned test experiment, which consists of the C code and the S code with a length of 1024 chips Combined in the form of real and imaginary parts of complex numbers. It can be seen from the figure that within any time delay, the real part of the correlation function of this combined code has ideal autocorrelation and cross-correlation values, which is consistent with the theoretical results described in formula (11).
如图7和图8示出了根据上述湖试实验,两种“互补序列”对组合码的扩频信号在做CZT时频搜索得到的三维结果。可以看出,在搜索结果的主峰外的区域非常平坦,这意味着着两种扩频信号能具有很好的相关特性,能达到很好的多途信号抑制能力。Figure 7 and Figure 8 show the three-dimensional results obtained by CZT time-frequency search of two "complementary sequences" on the spread spectrum signal of the combined code according to the above-mentioned lake test. It can be seen that the area outside the main peak of the search results is very flat, which means that the two spread spectrum signals can have good correlation characteristics and can achieve good multi-path signal suppression capabilities.
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than limit them. Although the present invention has been described in detail with reference to the embodiments, those skilled in the art should understand that modifications or equivalent replacements to the technical solutions of the present invention do not depart from the spirit and scope of the technical solutions of the present invention, and all of them should be included in the scope of the present invention. within the scope of the claims.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310739788.3A CN104753561B (en) | 2013-12-26 | 2013-12-26 | Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310739788.3A CN104753561B (en) | 2013-12-26 | 2013-12-26 | Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104753561A CN104753561A (en) | 2015-07-01 |
CN104753561B true CN104753561B (en) | 2017-05-17 |
Family
ID=53592747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310739788.3A Active CN104753561B (en) | 2013-12-26 | 2013-12-26 | Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104753561B (en) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106470046A (en) * | 2015-08-14 | 2017-03-01 | 宜春市等比科技有限公司 | A kind of utilization 4b/6b technology and the method with double pseudo-random code spread-spectrum |
CN105407066B (en) * | 2015-12-22 | 2018-07-20 | 中国船舶重工集团公司第七一五研究所 | A method of selection optimal pilot sequence inhibits underwater sound ofdm signal peak-to-average ratio |
CN105656620A (en) * | 2016-02-19 | 2016-06-08 | 上海果路交通科技有限公司 | Pseudorandom sequence generating method |
CN107276928B (en) * | 2016-04-08 | 2021-09-03 | 深圳光启合众科技有限公司 | Channel estimation method and device |
US10003377B1 (en) | 2016-12-19 | 2018-06-19 | Cisco Technology, Inc. | Spread spectrum acoustic communication techniques |
US10141973B1 (en) | 2017-06-23 | 2018-11-27 | Cisco Technology, Inc. | Endpoint proximity pairing using acoustic spread spectrum token exchange and ranging information |
CN108429591B (en) * | 2017-11-13 | 2021-05-18 | 西北工业大学 | Multi-carrier underwater acoustic communication method suitable for deep sea channel |
CN110266622B (en) * | 2018-03-12 | 2020-12-01 | 中国科学院声学研究所 | An Orthogonal Multi-Carrier M-Element Chaotic Phase Modulation Spread Spectrum Underwater Acoustic Communication Method |
US10404319B1 (en) | 2018-07-09 | 2019-09-03 | Cisco Technology, Inc. | Fast correlation of prometheus orthonormal sets (PONS) for communications |
US10396846B1 (en) | 2018-10-12 | 2019-08-27 | Cisco Technology, Inc. | Adaptive decoding of spread spectrum signals using multiple correlator peaks |
US10601459B1 (en) | 2018-11-02 | 2020-03-24 | Cisco Technology, Inc. | Efficient handling of clock offset in spread spectrum decoders |
CN111756470B (en) * | 2020-06-23 | 2024-05-14 | 中国人民解放军63892部队 | Noise modulation method based on pseudo-random sequence |
CN114705883B (en) * | 2022-02-21 | 2024-07-05 | 天津大学 | Layered flow measurement method based on complementary coding |
CN118921129B (en) * | 2024-10-11 | 2024-12-27 | 广东海洋大学 | An underwater acoustic wave data communication system |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079674A (en) * | 2007-04-25 | 2007-11-28 | 北京交通大学 | A water sound communication method of shallow sea remote horizontal channel |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006157643A (en) * | 2004-11-30 | 2006-06-15 | Naoki Suehiro | Radio communications system, radio communication method and communications equipment |
-
2013
- 2013-12-26 CN CN201310739788.3A patent/CN104753561B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079674A (en) * | 2007-04-25 | 2007-11-28 | 北京交通大学 | A water sound communication method of shallow sea remote horizontal channel |
Non-Patent Citations (2)
Title |
---|
"Analysis of Channel Effects on Direct-Sequence and Frequency-Hopped Spread-Spectrum Acoustic Communication";Lee Freitag, Milica Stojanovic, Sandipa Singh, and Mark Johnson;《IEEE JOURNAL OF OCEANIC ENGINEERING》;20011031;第26卷(第4期);586-593 * |
"LS 码在水声扩频中的应用及估计";张帅,高勇;《通信技术》;20131130;第46卷(第11期);20-24 * |
Also Published As
Publication number | Publication date |
---|---|
CN104753561A (en) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104753561B (en) | Direct sequence spread spectrum modulation method for suppressing multipath interference in underwater acoustic communication | |
CA2818775C (en) | Zero correlation zone sequences for communication system | |
CA2818766C (en) | Zero correlation zone sequences for communication system | |
CN103501201B (en) | A kind of frequency hopping arteries and veins position based on linear FM signal coding underwater acoustic communication method | |
CN101567727A (en) | Differential cyclic shift spread-spectrum underwater sound communication method | |
CN101471687B (en) | User terminal for wireless sensing network based on m sequence self correlation | |
CN101277125B (en) | Method and apparatus for generating zero correlation section code, transmitting and receiving spread-spectrum code | |
CN102170314A (en) | Hyperbolic frequency-modulation spread spectrum acoustic communication method | |
CN103684521A (en) | Fast and accurate synchronization method for spread spectrum underwater acoustic communication | |
CN102833008B (en) | Spread spectrum underwater acoustic communication method | |
CN101534278A (en) | Time-frequency expansion Orthogonal Frequency Division Multiplexing transmitting and receiving device, method and system | |
CN105915480A (en) | Efficient chaotic communication scheme based on orthogonal chaotic generator | |
CN103905085A (en) | Burst hybrid spread spectrum underwater sound covert communication method | |
CN105933104A (en) | Multi-user difference chaotic communication system based on Walsh codes | |
CN101355374B (en) | Method for generating signal of non-interference quasi-synchronous CDMA communication system | |
CN113242198B (en) | PI-DCSK modulation and demodulation method and system based on Walsh codes | |
CN100433600C (en) | CDMA spectrum-spreading method, despreading method and receiver | |
CN107786296A (en) | A kind of reconstructing method for GOLD sequences | |
CN102104576A (en) | A Mixed Carrier Direct Sequence Spread Spectrum CDMA Transmission System | |
CN104779993B (en) | Deep space TT&C system and method based on frequency domain equalization | |
CN103095335B (en) | Parallel combined spread spectrum communication method of variable combination | |
CN100452789C (en) | Short-distance radio network data transmitting method | |
Zhang et al. | An efficient spread spectrum pulse position modulation scheme for point-to-point underwater acoustic communication | |
Cotae et al. | A cyclostationary blind Doppler estimation method for underwater acoustic communications using direct-sequence spread spectrum signals | |
CN203872203U (en) | Code division multiplexing quadrature frequency division multiple access communication system signal emitting device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |