CN110830077A - A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals - Google Patents
A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals Download PDFInfo
- Publication number
- CN110830077A CN110830077A CN201910913971.8A CN201910913971A CN110830077A CN 110830077 A CN110830077 A CN 110830077A CN 201910913971 A CN201910913971 A CN 201910913971A CN 110830077 A CN110830077 A CN 110830077A
- Authority
- CN
- China
- Prior art keywords
- code
- pseudo
- data
- zero
- block
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 238000001514 detection method Methods 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 11
- 230000018199 S phase Effects 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000003780 insertion Methods 0.000 claims description 2
- 230000037431 insertion Effects 0.000 claims description 2
- 238000004891 communication Methods 0.000 abstract description 24
- 238000013461 design Methods 0.000 description 7
- 238000013480 data collection Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012966 insertion method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7075—Synchronisation aspects with code phase acquisition
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/16—Circuits
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
一种改善多路突发信号接收性能的快速捕获方法,用于解决报文通信用户使用特殊伪码进行信号调制的捕获问题。该方法首先进行零中频数据与伪码的采集,然后将零中频数据与伪码进行划分块处理,利用双块拓展方法完成零中频数据块与伪码块的拓展;根据零中频数据与主扩频伪码相关计算结果的排列规律,完成相位并行搜索与频率并行搜索,进而完成所有伪码相位的捕获。本发明捕获方法极大地缩短了捕获的搜索时间,同时捕获性能优异,使这种特殊的伪码结构能够应用于实际工程。
A fast acquisition method for improving the receiving performance of multi-channel burst signals is used to solve the acquisition problem that a message communication user uses a special pseudo code for signal modulation. The method first collects the zero-IF data and pseudo-code, then divides the zero-IF data and pseudo-code into blocks, and uses the double-block expansion method to complete the expansion of the zero-IF data block and the pseudo-code block; According to the arrangement rule of frequency pseudo-code correlation calculation results, phase parallel search and frequency parallel search are completed, and then the acquisition of all pseudo-code phases is completed. The acquisition method of the invention greatly shortens the acquisition search time, and at the same time has excellent acquisition performance, so that this special pseudo-code structure can be applied to practical engineering.
Description
技术领域technical field
本发明涉及一种改善多路突发信号接收性能的快速捕获方法,属于扩频通信技术领域。The invention relates to a fast acquisition method for improving multi-channel burst signal receiving performance, and belongs to the technical field of spread spectrum communication.
背景技术Background technique
报文通信接收机的信号格式为短突发模式,信号持续时间长短可变。该信号由同步头,勤务段与数据段构成。同步头由周期短码组成,不进行数据调制,用于报文通信接收机捕获信号。报文通信接收机接收各个终端发送的报文通信信号,其中各个用户同步头部分均由相同的周期短码组成。当多个报文通信用户同时发送入站信号时,存在用户碰撞的问题,同时用户间多址抑制能力有限。针对该问题,设计了一种新的伪码。该伪码设计基于两个周期较短的扩频码实现,分别称为主扩频码和副扩频码。The signal format of the message communication receiver is short burst mode, and the signal duration is variable. The signal consists of a sync header, a service segment and a data segment. The synchronization header is composed of periodic short codes without data modulation, and is used to capture signals by the message communication receiver. The message communication receiver receives message communication signals sent by each terminal, wherein each user synchronization header part is composed of the same period short code. When multiple message communication users send inbound signals at the same time, there is a problem of user collision, and the ability to suppress multiple access between users is limited. To solve this problem, a new pseudocode is designed. The pseudo-code design is based on two short-period spreading codes, which are called the main spreading code and the secondary spreading code respectively.
周期性重复主扩频码,通过在每个主扩频码周期后插入不等长的副扩频码,构成特殊的伪码结构。与恒周期重复码不同,在每个主扩频码周期后插入了若干个码片的副扩频码序列,通过合理设计每个周期后插入的扩频码码片数,可以有效地破坏同步头的周期性,从而降低不同用户同时入站碰撞的概率,具体插入方式见图1所示。The main spreading code is repeated periodically, and a special pseudo-code structure is formed by inserting sub-spreading codes of unequal lengths after each period of the main spreading code. Different from the constant cycle repeating code, a sub-spreading code sequence of several chips is inserted after each main spreading code cycle. By rationally designing the number of spreading code chips inserted after each cycle, the synchronization can be effectively destroyed. The periodicity of the header reduces the probability of simultaneous inbound collisions by different users. The specific insertion method is shown in Figure 1.
在报文通信信号特殊伪码结构设计中,共包含X1个周期的主扩频伪码,主扩频码间插入副扩频码的个数满足等差数列的关系,公差数为d。这种伪码结构设计可以解决多用户同时入站的碰撞问题,同时改善多用户间干扰抑制能力。但是,对于报文通信接收机来说,新的伪码结构破坏了周期伪码的特性,给捕获系统的设计提出新的要求。现有的捕获方法不适用特殊伪码结构的捕获。In the special pseudo-code structure design of the message communication signal, a total of X1 cycles of the main spreading pseudo-code are included. This pseudo-code structure design can solve the collision problem of multiple users entering the station at the same time, and at the same time improve the interference suppression ability between multiple users. However, for the message communication receiver, the new pseudo-code structure destroys the characteristics of periodic pseudo-code, which puts forward new requirements for the design of the acquisition system. Existing capture methods are not suitable for capture of special pseudo-code structures.
发明内容SUMMARY OF THE INVENTION
本发明的技术解决问题是:克服现有技术的不足,提供一种改善多路突发信号接收性能的快速捕获方法,实现特殊伪码结构的调制信号快速捕获,任何时刻入站信号的捕获概率均大于99%,且捕获时间较短。The technical solutions of the present invention are: overcoming the deficiencies of the prior art, providing a fast acquisition method for improving the receiving performance of multi-channel burst signals, realizing fast acquisition of modulated signals with special pseudo-code structure, and acquiring probability of inbound signals at any time. All are greater than 99%, and the capture time is short.
本发明的技术解决方案是:The technical solution of the present invention is:
一种改善多路突发信号接收性能的快速捕获方法,包括如下步骤:A fast acquisition method for improving multi-channel burst signal receiving performance, comprising the following steps:
(1)零中频数据与伪码采集(1) Zero-IF data and pseudo-code acquisition
根据捕获开始标识产生采集起始时刻,以2倍的伪码速率完成对零中频数据与伪码的采集,采集结束后,生成采集结束标识;所述伪码为主扩频码,采集的零中频数据存储在数据区,采集的伪码存储在伪码区,采集起始时刻对应零中频数据第S个主扩频码;The acquisition start time is generated according to the acquisition start flag, and the acquisition of the zero-IF data and the pseudocode is completed at twice the pseudocode rate. After the acquisition is completed, the acquisition end flag is generated; The intermediate frequency data is stored in the data area, the collected pseudo code is stored in the pseudo code area, and the acquisition start time corresponds to the S-th main spreading code of the zero intermediate frequency data;
(2)双块拓展(2) Double block expansion
采集结束后,根据零中频数据在数据区存储的位置以及伪码在伪码区存储的位置,将零中频数据与伪码进行划分块处理,根据双块拓展方法,完成零中频数据块与伪码块的拓展;其中零中频数据与伪码划分的块数相同;After the acquisition, according to the location where the zero-IF data is stored in the data area and the location where the pseudo-code is stored in the pseudo-code area, the zero-IF data and the pseudo-code are divided into blocks. The expansion of the code block; wherein the zero-IF data is the same as the number of blocks divided by the pseudo code;
(3)相位并行搜索(3) Phase parallel search
对拓展后的每块零中频数据和每块伪码进行FFT运算,将第k块零中频数据与第k块伪码的FFT运算结果进行共轭相乘运算,然后对运算结果进行IFFT运算,得到零中频数据与伪码的相关计算结果,存储在搜索区中;k=1,2,……,K,K为步骤(2)中零中频数据或伪码的划分块数;Perform FFT operation on each block of zero-IF data after expansion and each block of pseudocode, perform conjugate multiplication operation on the kth block of zero-IF data and the FFT operation result of the kth block of pseudocode, and then perform IFFT operation on the operation result, Obtain the correlation calculation result of the zero-IF data and the pseudo-code, and store in the search area; k=1, 2, .
(4)频率并行搜索(4) Frequency parallel search
根据零中频数据与伪码的相关计算结果在搜索区中的排列规律,获取相同伪码相位计算的相关结果进行FFT运算,完成每一个伪码相位对应的频率搜索;According to the arrangement rule of the correlation calculation result between the zero intermediate frequency data and the pseudo code in the search area, obtain the correlation result of the same pseudo code phase calculation and perform FFT operation to complete the frequency search corresponding to each pseudo code phase;
(5)捕获判决(5) Capture judgment
根据步骤(4)中FFT运算输出的结果,进行能量计算,同时计算噪声的功率,根据捕获检测概率与虚警概率,设定检测门限,完成所有伪码相位的捕获判决。According to the result of the FFT operation output in step (4), the energy calculation is performed, and the power of the noise is calculated at the same time.
所述步骤(1)中,设零中频数据采集时间为T毫秒,则:In the described step (1), set the zero-IF data acquisition time to be T milliseconds, then:
设原始的零中频数据共包含X1个主扩频码,T毫秒共采集到X2个主扩频码,则X2<X1;Assuming that the original zero-IF data contains a total of X1 main spreading codes, and a total of X2 main spreading codes are collected in T milliseconds, then X2<X1;
伪码采集时间根据伪码的周期决定,如果伪码周期长度大于T毫秒,伪码采集时间等于T毫秒;如果伪码周期长度小于等于T毫秒,伪码采集时间等于一个主扩频码周期长度。The pseudo code collection time is determined according to the pseudo code period. If the pseudo code period length is greater than T milliseconds, the pseudo code collection time is equal to T milliseconds; if the pseudo code period length is less than or equal to T milliseconds, the pseudo code collection time is equal to the length of a main spreading code period. .
所述步骤(2)中,对伪码进行划分块处理的方法如下:In the step (2), the method for dividing the pseudocode into blocks is as follows:
设主扩频码用s表示,一个主扩频码周期包含Q个码片,以2倍的伪码速率完成采集后共得到2Q个半码片,组成一个伪码块;Suppose the main spreading code is represented by s, one main spreading code period contains Q chips, and after the acquisition is completed at twice the pseudo-code rate, a total of 2Q half-chips are obtained to form a pseudo-code block;
各个伪码块的取值相同,令sk(i)表示第k个伪码块内第i个半码片,则:The values of each pseudo-code block are the same, let s k (i) represent the i-th half-chip in the k-th pseudo-code block, then:
sk(i)=s(i)1≤i≤2Qs k (i)=s(i)1≤i≤2Q
其中s(i)表示任意一个伪码块的第i个半码片。where s(i) represents the i-th half-chip of any pseudo-code block.
所述步骤(2)中,对零中频数据进行划分块处理的方法如下:In the step (2), the method for dividing the zero-IF data into blocks is as follows:
一个零中频数据块的长度与一个伪码块的长度相同,包含2Q个采样数据,在T毫秒采集时间内,采集得到的零中频数据用r表示,d表示主扩频码插入副扩频码的公差数,rk(i)表示第k个零中频数据块内第i个数据,则The length of a zero-IF data block is the same as that of a pseudo-code block, including 2Q sampling data. During the acquisition time of T milliseconds, the collected zero-IF data is represented by r, and d represents the insertion of the main spreading code into the sub-spreading code. The tolerance number of , r k (i) represents the i-th data in the k-th zero-IF data block, then
rk(i)=r(2Qk-2Q-2d+kd+i)1≤i≤2Q。r k (i)=r(2Qk−2Q−2d+kd+i)1≤i≤2Q.
所述步骤(2)中,按照如下方法完成伪码块的拓展:In the described step (2), the expansion of the pseudocode block is completed according to the following method:
伪码块拓展后的长度为8Q,表示sk(i)拓展后的伪码块,拓展方法如下所示:The extended length of the pseudo-code block is 8Q, represents the extended pseudo-code block of s k (i), and the extension method is as follows:
sk(i)为第k个伪码块内第i个半码片。s k (i) is the i-th half-chip in the k-th pseudo-code block.
所述步骤(2)中,按照如下方法完成零中频数据块的拓展:In the described step (2), the expansion of the zero-IF data block is completed according to the following method:
零中频数据块拓展后的长度为8Q,表示rk(i)拓展后的零中频数据块,拓展方法如下所示:The extended length of the zero-IF data block is 8Q, Represents the zero-IF data block after expansion of r k (i). The expansion method is as follows:
rk(i)为第k个零中频数据块内的第i个数据。r k (i) is the i-th data in the k-th zero-IF data block.
所述步骤(4)中,对每一个伪码相位对应的频率进行搜索的实现方法如下:In the described step (4), the implementation method of searching the frequency corresponding to each pseudo code phase is as follows:
采集的零中频数据共包含X2个主扩频码,由于采集起始时刻对应零中频数据第S个主扩频码,根据零中频数据与伪码的相关计算结果在搜索区中的排列规律,相同伪码相位计算得到的相关值依次向后移动2*d*S个相位;The collected zero-IF data contains a total of X2 main spreading codes. Since the acquisition start time corresponds to the S-th main spreading code of the zero-IF data, according to the arrangement rule of the correlation calculation result between the zero-IF data and the pseudo code in the search area, The correlation values calculated by the same pseudo-code phase are shifted backward by 2*d*S phases in sequence;
对采集时间对应的原始数据任何一个主扩频周期进行遍历搜索,挑选出相同伪码相位计算得到的相关值,然后进行FFT计算,完成每一个伪码相位对应的频率搜索;其中相同伪码相位计算得到的相关值共有X2个,与主扩频码个数相同。Perform a traversal search on any main spread spectrum cycle of the original data corresponding to the acquisition time, select the correlation value calculated by the same pseudo-code phase, and then perform FFT calculation to complete the frequency search corresponding to each pseudo-code phase; among them, the same pseudo-code phase There are X2 correlation values in total, which are the same as the number of main spreading codes.
FFT计算点数最好满足2的n次方,n为整数,满足2n-X2≥0且取最小值,如果X2小于2n,在相同伪码相位计算得到的相关值后,添加2n-X2个零进行FFT计算;如果X2等于2n,直接进行FFT计算。The number of FFT calculation points is preferably 2 to the nth power, n is an integer, 2 n -X2≥0 and takes the minimum value, if X2 is less than 2 n , after the correlation value calculated by the same pseudo code phase, add 2 n - X2 zeros are used for FFT calculation; if X2 is equal to 2 n , FFT calculation is performed directly.
所述步骤(5)中根据FFT运算输出的结果行能量计算的实现方式如下:In the described step (5), the implementation mode of calculating the row energy according to the result of the FFT operation output is as follows:
根据FFT运算输出的同向分量与正交分量,计算捕获检测量I:According to the in-direction component and quadrature component output by the FFT operation, the captured detection quantity I is calculated:
其中,IPi为FFT运算输出的第i个同向分量,QPi为FFT运算输出的第i个正交分量;Wherein, IP i is the i-th in-direction component output by the FFT operation, and QP i is the i-th quadrature component output by the FFT operation;
在仅有噪声存在下,根据捕获检测量I计算得到噪声方差则仅有噪声时,捕获检测量I符合标准的中心分布;若信号存在,则捕获检测量I服从非中心分布,非中心分布的参数λ=2X*snr,X等于1,snr为相干累加后的信噪比。In the presence of only noise, the noise variance is calculated according to the captured detection quantity I Then when there is only noise, the capture detection quantity I meets the center of the standard distribution; if the signal exists, the capture detection quantity I obeys the non-centrality distribution, non-central The distribution parameter λ=2X*snr, X is equal to 1, and snr is the signal-to-noise ratio after coherent accumulation.
所述步骤(5)中,检测门限Vt计算如下:In the step (5), the detection threshold V t is calculated as follows:
Pfa为虚警概率。 Pfa is the false alarm probability.
本发明与现有技术相比的优点在于:The advantages of the present invention compared with the prior art are:
(1)针对特殊的伪码结构,采用主扩频信号进行能量的累加,丢弃副扩频码的信号能量;根据在主扩频码插入副扩频码的规律,对零中频数据完成块的划分,该划分利用主扩频码插入的副扩频码的特征,且采样速率是伪码速率的两倍,划分的零中频数据块的起始时刻顺次向后进行偏移,这种划分方法为后续FFT快速计算提供了条件。(1) For the special pseudo-code structure, the main spread spectrum signal is used to accumulate energy, and the signal energy of the sub-spread spectrum code is discarded; The division uses the characteristics of the sub-spreading code inserted by the main spreading code, and the sampling rate is twice the rate of the pseudo code. The method provides conditions for the fast calculation of the subsequent FFT.
(2)对零中频数据块与伪码块进行双块拓展运算,这种拓展运算方式保证了零中频数据与伪码进行FFT快速计算的正确性,可以迅速完成伪码相位的并行搜索。(2) Double-block expansion operation is performed on the zero-IF data block and the pseudo-code block. This expansion operation method ensures the correctness of the fast FFT calculation between the zero-IF data and the pseudo-code, and can quickly complete the parallel search of the pseudo-code phase.
(3)针对特殊伪码结构的捕获方法,每次计算的相同码相位对应的相关值的位置是变化,且这种变换满足一定的规律。捕获数据采集的起始时刻位于第S个主扩频伪码区间,计算零中频数据与相同伪码相位的相关值依次向后移动2*d*S的相位。根据这种排列规律,读取相同伪码相位计算的相关结果进行FFT运算,完成每一个伪码相位对应的频率搜索。频率的并行搜索大量节约了串行搜索的捕获计算时间。(3) For the acquisition method of the special pseudo-code structure, the position of the correlation value corresponding to the same code phase calculated each time is changed, and this transformation satisfies a certain law. The start time of the captured data collection is located in the S-th main spread spectrum pseudo-code interval, and the correlation value between the zero-IF data and the same pseudo-code phase is calculated and the phase is shifted backward by 2*d*S in turn. According to this arrangement rule, the correlation results of the same pseudo-code phase calculation are read and the FFT operation is performed to complete the frequency search corresponding to each pseudo-code phase. The parallel search of frequency saves the acquisition computation time of the serial search by a large amount.
附图说明Description of drawings
图1为特殊结构伪码结构设计;Figure 1 is a special structure pseudo code structure design;
图2为特殊结构伪码结构构成;Fig. 2 is a special structure pseudo-code structure composition;
图3为本发明流程图;Fig. 3 is the flow chart of the present invention;
图4报文通信信号特殊结构伪码序列;Fig. 4 message communication signal special structure pseudo code sequence;
图5零中频数据块的划分;Figure 5 is divided into zero-IF data blocks;
图6进行扩展后的零中频数据块;Fig. 6 carries out the zero intermediate frequency data block after extension;
图7进行扩展后的伪码块。Fig. 7 Pseudo-code block after extension.
具体实施方式Detailed ways
多个报文通信用户同时发送入站信号时,存在用户碰撞的问题,且用户间多址抑制能力有限,使用图2所示的扩频码可以解决上述问题。该伪码设计基于两个周期较短的扩频码实现,分别称为主扩频码和副扩频码。周期性重复主扩频码,通过在每个主扩频码周期后插入不等长的副扩频码,构成特殊的伪码结构。每个报文通信用户用这种特殊的伪码调制同步头信息进行入站,该同步头包含X1个主扩频码,一个主扩频码周期包含Q个码片,d表示主扩频码插入副扩频码的公差数。报文通信接收机接收报文通信用户的入站信息,通过同步头进行信号捕获。报文通信接收机接收的中频信号通过AD采样后输入至数字下变频模块,进行数字下变频,得到零中频数据送至低通滤波器,该零中频数据仅包含多普勒信息。When multiple message communication users send inbound signals at the same time, there is the problem of user collision, and the multiple access suppression capability between users is limited. The spreading code shown in FIG. 2 can solve the above problem. The pseudo-code design is based on two short-period spreading codes, which are called the main spreading code and the secondary spreading code respectively. The main spreading code is repeated periodically, and a special pseudo-code structure is formed by inserting sub-spreading codes of unequal lengths after each period of the main spreading code. Each message communication user uses this special pseudo code to modulate the synchronization header information for inbound, the synchronization header contains X1 main spreading codes, one main spreading code cycle contains Q chips, and d represents the main spreading code Tolerance number for inserting subspreading codes. The message communication receiver receives the inbound information of the message communication user, and performs signal acquisition through the synchronization header. The IF signal received by the message communication receiver is sampled by AD and input to the digital down-conversion module for digital down-conversion to obtain zero-IF data and send it to the low-pass filter. The zero-IF data only contains Doppler information.
本发明即是对通过低通滤波处理后的零中频数据进行捕获计算,捕获计算核心要解决同步头伪码插入副扩频码破坏伪码周期性的问题,设计新的捕获方法进行报文通信信号的捕获,从而改善报文通信接收机接收多路突发信号的性能。The invention is to capture and calculate the zero-intermediate frequency data processed by the low-pass filtering. The core of the capture and calculation is to solve the problem that the pseudo code of the synchronization header is inserted into the sub-spread spectrum code to destroy the periodicity of the pseudo code, and a new capture method is designed for message communication. Signal acquisition, thereby improving the performance of the message communication receiver to receive multiple burst signals.
报文通信接收机接收采样的单路中频信号为:The single-channel IF signal received and sampled by the message communication receiver is:
上式中,tj为第j个采样时刻,本地接收信号的采样速率为伪码速率的两倍;τ为信号到达报文通信接收机的时延;A为信号幅度,fL1为射频载波频率,fd为多普勒频移,fIF为中频载波频率,为初始相位,C(t)为伪码,n(tj)为第j个采样时刻的加性高斯白噪声。In the above formula, t j is the jth sampling time, and the sampling rate of the locally received signal is twice the pseudo code rate; τ is the time delay of the signal reaching the message communication receiver; A is the signal amplitude, and f L1 is the radio frequency carrier frequency, f d is the Doppler shift, f IF is the intermediate frequency carrier frequency, is the initial phase, C(t) is the pseudo code, and n(t j ) is the additive white Gaussian noise at the jth sampling time.
将中频信号进行数字下变频,得到零中频数据送至低通滤波器,最终得到可以用于捕获的零中频数据,原始的零中频数据共包含X1个主扩频码,采集后的零中频数据用r(x)表示,该零中频数据仅包含多普勒信息,x表示采集的第x个数据。The IF signal is digitally down-converted, and the zero-IF data is sent to the low-pass filter, and finally the zero-IF data that can be used for capture is obtained. The original zero-IF data contains X1 main spreading codes in total. The collected zero-IF data It is represented by r(x), the zero-IF data only contains Doppler information, and x represents the xth data collected.
针对这种特殊伪码结构,设计的快速捕获方法如图3所示,步骤如下:For this special pseudo-code structure, the designed fast capture method is shown in Figure 3, and the steps are as follows:
(1)零中频数据与伪码采集:根据捕获开始标识产生采样起始时刻,以2倍的伪码速率完成对零中频数据与伪码的采集,此处的伪码为主扩频码,用于接收信号的捕获计算;采集起始时刻对应零中频数据第S个主扩频码。(1) Acquisition of zero-IF data and pseudo-code: The sampling start time is generated according to the capture start flag, and the collection of zero-IF data and pseudo-code is completed at twice the pseudo-code rate. The pseudo-code here is the main spreading code, It is used for the acquisition calculation of the received signal; the acquisition start time corresponds to the S-th main spreading code of the zero-IF data.
根据捕获开始标识产生采集起始时刻,没有时间基准,产生的时刻是随机的,零中频数据共包含X1个主扩频码,采集结果存储分别存储在数据区与伪码区中,每个区存储地址对应一个采集结果,设零中频数据采集时间为T毫秒,共包含X2个主扩频码,X2<X1;伪码采集时间根据伪码的周期决定,如果伪码周期长度大于T毫秒,伪码采集时间等于T毫秒;如果伪码周期长度小于等于T毫秒,伪码采集时间等于一个主扩频码周期长度;零中频数据与伪码采集结束后,生成采集结束标识。The acquisition start time is generated according to the acquisition start flag. There is no time reference, and the generated time is random. The zero-IF data contains X1 main spreading codes in total. The acquisition results are stored in the data area and the pseudo code area respectively. The storage address corresponds to a collection result, and the zero-IF data collection time is set to T milliseconds, including X2 main spreading codes, X2<X1; the pseudo code collection time is determined according to the pseudo code period, if the pseudo code period length is greater than T milliseconds, The pseudo code collection time is equal to T milliseconds; if the pseudo code period length is less than or equal to T milliseconds, the pseudo code collection time is equal to the period length of one main spreading code; after the zero-IF data and pseudo code collection is completed, the collection end flag is generated.
(2)双块拓展:根据零中频数据在数据区存储的位置以及伪码在伪码区存储的位置,将零中频数据与伪码进行划分块处理,其中零中频数据与伪码划分的块数相同,零中频数据块与伪码块组成双块;根据双块拓展的方法,完成零中频数据块与伪码块的划分与拓展;(2) Double-block expansion: According to the location where the zero-IF data is stored in the data area and the location where the pseudo-code is stored in the pseudo-code area, the zero-IF data and the pseudo-code are divided into blocks, where the zero-IF data and the pseudo-code are divided into blocks. The number is the same, the zero-IF data block and the pseudo-code block form a double block; according to the method of double-block expansion, the division and expansion of the zero-IF data block and the pseudo-code block are completed;
主扩频码用s表示,一个主扩频码周期为Q个码片,以2倍的伪码速率完成采集后共2Q个半码片,组成一个伪码块,s(i)表示一个伪码块的第i的半码片,k表示伪码块的块号,sk(i)表示第k个伪码块内第i个半码片,所有的伪码块取值相同,用如下公式表示:The main spreading code is represented by s. A main spreading code period is Q chips. After the acquisition is completed at twice the pseudo code rate, a total of 2Q half chips form a pseudo code block, and s(i) represents a pseudo code block. The i-th half-chip of the code block, k represents the block number of the pseudo-code block, s k (i) represents the i-th half-chip in the k-th pseudo-code block, and all pseudo-code blocks have the same value, which is used as follows The formula says:
sk(i)=s(i)1≤i≤2Qs k (i)=s(i)1≤i≤2Q
一个零中频数据块的长度与一个伪码块的长度相同,包含2Q个采样数据,在T毫秒时间内,采集得到的零中频数据用r表示,r(x)表示第x个零中频数据;k表示零中频数据块的块号;d表示主扩频码插入副扩频码的公差数,rk(i)表示第k个零中频数据块内第i个数据,用2Qk-2Q-2d+kd+i计算得到零中频数据的位置,可用如下公式表示:The length of a zero-IF data block is the same as that of a pseudo-code block, and contains 2Q sample data. Within T milliseconds, the collected zero-IF data is represented by r, and r(x) represents the xth zero-IF data; k represents the block number of the zero-IF data block; d represents the tolerance number of the main spreading code inserted into the sub-spreading code, r k (i) represents the i-th data in the k-th zero-IF data block, using 2Qk-2Q-2d +kd+i calculates the position of zero intermediate frequency data, which can be expressed by the following formula:
rk(i)=r(2Qk-2Q-2d+kd+i)1≤i≤2Qr k (i)=r(2Qk-2Q-2d+kd+i)1≤i≤2Q
k=1,2,……,K,K为零中频数据或伪码的划分块总数。k=1, 2, ..., K, K is the total number of divided blocks of zero-IF data or pseudocode.
零中频数据块与伪码块进行双块拓展后的长度为8Q,表示sk(i)拓展后的伪码块,拓展方法如公式所示:The length of the zero-IF data block and the pseudo-code block after double-block expansion is 8Q, represents the expanded pseudo-code block of s k (i), and the expansion method is shown in the formula:
表示rk(i)拓展后的零中频数块,拓展方法如公式所示 Represents the zero-IF block after expansion of r k (i), the expansion method is shown in the formula
(3)相位并行搜索:对拓展后的每块零中频数据和每块伪码分别进行FFT运算,将第k块零中频数据与第k块伪码的FFT运算结果进行共轭相乘运算,然后对运算结果进行IFFT运算,得到零中频数据与伪码的相关计算结果存储在搜索区中,为主扩频码每个相位的位置。(3) Phase parallel search: perform FFT operation on each block of zero-IF data after expansion and each block of pseudocode, and perform conjugate multiplication operation on the kth block of zero-IF data and the FFT operation result of the kth block of pseudocode, Then perform IFFT operation on the operation result to obtain the correlation calculation result of zero-IF data and pseudocode stored in the search area, is the position of each phase of the main spreading code.
按照上述步骤与方法,完成(2)中所有零中频数据块与伪码块的计算与结果存储。According to the above steps and methods, the calculation and result storage of all zero-IF data blocks and pseudo-code blocks in (2) are completed.
(4)频率并行搜索:捕获开始标识产生时刻未知,零中频数据与伪码采集起始时刻位于主扩频码的任何位置,需要对采集起始时刻位于主扩频码任何位置进行考虑。根据零中频数据与伪码相关计算在区中的排列规律,读取相同伪码相位计算的相关结果进行FFT运算,完成每一个伪码相位对应的频率搜索;(4) Frequency parallel search: The acquisition start time is unknown, and the acquisition start time of zero intermediate frequency data and pseudo code is located at any position of the main spreading code. It is necessary to consider that the acquisition start time is located at any position of the main spreading code. According to the arrangement rule of zero intermediate frequency data and pseudo code correlation calculation in the area, read the correlation result of the same pseudo code phase calculation and perform FFT operation to complete the frequency search corresponding to each pseudo code phase;
采集存储的零中频数据共包含X2个主扩频码周期,采样起始时刻是随机的,采样起始时刻可以位于主扩频码中任何一个周期。零中频数据起始采样时刻在第S个主扩频,则相同伪码相位计算得到的相关值依次向后移动2*d*S个相位。对采集时间位于任何一个主扩频周期进行遍历,挑选出相同伪码相位计算得到的相关值,然后进行FFT计算;相同伪码相位计算得到的相关值共有X2个,与主扩频伪码周期数相同;FFT计算点数最好满足2的n次方(n为整数,满足2n-X2≥0且取最小值),如果X2小于2n,相同伪码相位计算得到的相关值后输入2n-X2个零进行FFT计算,如果X2等于2n,直接进行FFT计算。The zero-IF data collected and stored contains a total of X2 periods of the main spreading code. The sampling start time is random, and the sampling start time can be located in any period of the main spreading code. The initial sampling time of the zero-IF data is at the S-th main spread spectrum, and the correlation value calculated by the same pseudo-code phase is shifted backward by 2*d*S phases in sequence. Traverse the acquisition time in any main spread spectrum period, select the correlation value calculated by the same pseudo code phase, and then perform FFT calculation; the correlation value calculated by the same pseudo code phase has a total of X2, which is the same as the main spread spectrum pseudo code period. The number is the same; the number of FFT calculation points should preferably satisfy the nth power of 2 (n is an integer, 2 n -X2≥0 and take the minimum value), if X2 is less than 2 n , enter 2 after the correlation value calculated by the same pseudo code phase n -X2 zeros are used for FFT calculation. If X2 is equal to 2 n , FFT calculation is performed directly.
(5)捕获判决:根据FFT运算输出的结果,进行能量计算,同时计算噪声的功率;根据捕获检测概率与虚警概率,设定检测门限,完成所有伪码相位的捕获判决。(5) Acquisition judgment: According to the result of the FFT operation, the energy calculation is performed, and the power of the noise is calculated at the same time; according to the acquisition detection probability and the false alarm probability, the detection threshold is set to complete the acquisition judgment of all pseudo code phases.
根据FFT运算输出的同向分量与正交分量,计算捕获检测量I:According to the in-direction component and quadrature component output by the FFT operation, the captured detection quantity I is calculated:
其中,IPi为FFT运算输出的第i个同向分量,QPi为FFT运算输出的第i个正交分量;I服从自由度为2的分布。Among them, IP i is the i-th in-direction component output by the FFT operation, and QP i is the i-th quadrature component output by the FFT operation; distributed.
在仅有噪声存在下,根据捕获检测量I计算得到噪声方差则仅有噪声时I符合标准的中心分布,若信号存在,则I服从非中心分布,非中心分布参数为λ=2X*snr(进行了的归一化,X等于1,snr为相干累加后的信噪比),可以得到自由度为2的中心分布与非中心分布的概率分布,根据检测估值理论与上述计算关系,捕获门限Vt计算如下,虚警概率Pfa为已知量。In the presence of only noise, the noise variance is calculated according to the captured detection quantity I Then I only meet the center of the standard when there is only noise distribution, if the signal exists, then I obey the non-centrality distribution, the non-central distribution parameter is λ=2X*snr (performed The normalization of , X is equal to 1, and snr is the signal-to-noise ratio after coherent accumulation), the center with 2 degrees of freedom can be obtained Distribution and non-centrality The probability distribution of the distribution, according to the detection evaluation theory and the above calculation relationship, the acquisition threshold V t is calculated as follows, and the false alarm probability P fa is a known quantity.
为使本发明的目的、技术方案及优点更加清楚明白,用实际仿真数据对本发明进一步详细说明。In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention is further described in detail by using actual simulation data.
在报文通信项目中,对上述捕获方法进行性能仿真验证,测试条件如下:In the message communication project, the performance simulation verification of the above capture method is carried out, and the test conditions are as follows:
信号中心频点:85.9MHzSignal center frequency: 85.9MHz
主扩频码伪码周期数N2:50个The number of pseudo-code cycles of the main spreading code N2: 50
主扩频码一个周期码片的个数N1:1023The number of chips in one period of the main spreading code N1: 1023
伪码速率:1.023MHzPseudo code rate: 1.023MHz
副扩频码长度:1275个Sub-spreading code length: 1275
插入的副扩频码长度公差d:1Inserted sub-spreading code length tolerance d: 1
生成伪码速率X MHz:1.023MHzGenerate pseudocode rate X MHz: 1.023MHz
伪码周期:1msPseudocode period: 1ms
多普勒频率为300HzDoppler frequency is 300Hz
数据与伪码的采样频率:2.046MHzSampling frequency of data and pseudocode: 2.046MHz
信号入站频度:1s/次Signal inbound frequency: 1s/time
双拓展数据与伪码进行FFT-IFFT的点数:8192点Number of points for FFT-IFFT with double extended data and pseudocode: 8192 points
频率搜索FFT点数:32点Frequency search FFT points: 32 points
报文通信用户的入站方式为随机入站,到达报文通信接收机的时间满足泊松随机分布,报文通信接收机入口载噪比C/N0为36dBHz。同步头伪码采用上述特殊的伪码结构,主扩频采用周期长度为1023的m序列,副扩频码采用周期长度为2046的m序列,插入副扩频码长度公差d等于1,插入后的伪码结构如图4所示,组成报文通信信号特殊结构伪码序列,产生10万个入站信号完成上述捕获方法的仿真验证。The inbound mode of the message communication user is random inbound, the time to arrive at the message communication receiver satisfies the Poisson random distribution, and the ingress carrier-to-noise ratio C/N0 of the message communication receiver is 36dBHz. The synchronization header pseudo code adopts the above-mentioned special pseudo code structure, the main spread spectrum adopts the m sequence with the period length of 1023, the sub spread spectrum code adopts the m sequence with the period length of 2046, and the length tolerance d of the inserted sub spread spectrum code is equal to 1. The pseudo-code structure of Fig. 4 is shown in Figure 4, which constitutes a special-structure pseudo-code sequence of the message communication signal, and generates 100,000 inbound signals to complete the simulation verification of the above capture method.
数据的采集时刻通过随机数产生,matlab随机生成一个小于24的整数作为数据采集的起始时刻。零中频数据采集时间为26ms,包含25个完成的主扩频码;伪码采集时间根据伪码的周期决定,采集时间为1ms;以2倍的伪码速率完成对零中频数据与伪码的采集,采集的速率为2.046MHz;数据与伪码共划分25块,划分的零中频数据块如图5所示,进行拓展后的零中频数据块如图6所示,进行拓展后的伪码数据块如图7所示,扩展后的数据块的长度为8184。The data collection time is generated by random numbers, and matlab randomly generates an integer less than 24 as the start time of data collection. The zero-IF data collection time is 26ms, including 25 completed main spreading codes; the pseudo-code collection time is determined according to the pseudo-code cycle, and the collection time is 1ms; the zero-IF data and the pseudo-code are completed at twice the pseudo-code rate. Acquisition, the acquisition rate is 2.046MHz; the data and pseudocode are divided into 25 blocks, the divided zero-IF data block is shown in Figure 5, and the expanded zero-IF data block is shown in Figure 6, the expanded pseudocode The data block is shown in Figure 7, and the length of the extended data block is 8184.
对10万个入站信号进行仿真验证,用户成功接入99920个,接收成功率为99.92%,捕获时间约22ms,可以快速的完成多路随机接入信号的捕获。采用这种伪码结构设计,可以提高用户容量,同时多址抑制能力得到增强,使得改善多路突发信号的接收性能,使得这种特殊的伪码结构能够应用于实际工程。同时该方法可以根据项目的具体需求,可以进行参数化设计,对数据采集时间、FFT计算的点数,均根据实际需求可以变化,具有较高的扩展性。100,000 inbound signals are simulated and verified, and users successfully access 99,920 signals, the receiving success rate is 99.92%, and the acquisition time is about 22ms, which can quickly complete the acquisition of multi-channel random access signals. The use of this pseudo-code structure design can increase user capacity, and at the same time, the multiple access suppression capability is enhanced, which improves the receiving performance of multi-channel burst signals, and enables this special pseudo-code structure to be applied to practical engineering. At the same time, the method can be parameterized according to the specific needs of the project, and the data acquisition time and the number of FFT calculation points can be changed according to the actual needs, and it has high scalability.
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。The above is only the best specific embodiment of the present invention, but the protection scope of the present invention is not limited to this. Substitutions should be covered within the protection scope of the present invention.
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。Contents that are not described in detail in the specification of the present invention belong to the well-known technology of those skilled in the art.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910913971.8A CN110830077B (en) | 2019-09-25 | 2019-09-25 | A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910913971.8A CN110830077B (en) | 2019-09-25 | 2019-09-25 | A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110830077A true CN110830077A (en) | 2020-02-21 |
| CN110830077B CN110830077B (en) | 2021-07-13 |
Family
ID=69548289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910913971.8A Active CN110830077B (en) | 2019-09-25 | 2019-09-25 | A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110830077B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113225104A (en) * | 2021-04-12 | 2021-08-06 | 中国电子科技集团公司第三十八研究所 | Method and system for capturing multi-user burst spread spectrum signal in asynchronous communication system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080180321A1 (en) * | 2007-01-30 | 2008-07-31 | Chi-Shin Wang | Methods and systems for acquisition, reacquisition and tracking of weak navigational signals |
| CN105553507A (en) * | 2015-12-22 | 2016-05-04 | 北京理工大学 | Full coherence accumulation time-frequency domain parallel capturing method based on FFT |
| CN109617569A (en) * | 2018-11-28 | 2019-04-12 | 西安空间无线电技术研究所 | A Pseudocode Design Method for Improving the Receiving Performance of Multi-channel Burst Signals |
-
2019
- 2019-09-25 CN CN201910913971.8A patent/CN110830077B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080180321A1 (en) * | 2007-01-30 | 2008-07-31 | Chi-Shin Wang | Methods and systems for acquisition, reacquisition and tracking of weak navigational signals |
| CN105553507A (en) * | 2015-12-22 | 2016-05-04 | 北京理工大学 | Full coherence accumulation time-frequency domain parallel capturing method based on FFT |
| CN109617569A (en) * | 2018-11-28 | 2019-04-12 | 西安空间无线电技术研究所 | A Pseudocode Design Method for Improving the Receiving Performance of Multi-channel Burst Signals |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113225104A (en) * | 2021-04-12 | 2021-08-06 | 中国电子科技集团公司第三十八研究所 | Method and system for capturing multi-user burst spread spectrum signal in asynchronous communication system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110830077B (en) | 2021-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101777933B (en) | Generation and capture system of encrypted frame hopping spread spectrum signal of air fleet link | |
| CN106646546B (en) | A kind of satellite-signal multidimensional quick capturing method and system | |
| CN103199887B (en) | General capture method applied to direct spread spectrum signal | |
| CN101072044A (en) | Long-period spreading code frequency folding time frequency parallel searching method | |
| CN104765052B (en) | GEO navigation satellite high-sensitivity carrier tracking method | |
| CN101206255A (en) | Method for capturing GPS weak signal | |
| CN113972929B (en) | Method for capturing spread spectrum signal under high dynamic Doppler | |
| CN103809193B (en) | Improve the capture systems of Weak GNSS signal processing gain | |
| CN104483684A (en) | Method for rapidly capturing weak signals of Beidou D1 satellite navigation system | |
| CN111147102A (en) | Low signal-to-noise ratio code capturing method based on frequency modulation signal | |
| CN107888230B (en) | Multiphase coherent integration acquisition method applied to weak signal environment | |
| CN108196274A (en) | Be applicable in BOC (n, n) signal without fuzziness catching method and device | |
| CN103885071A (en) | Satellite signal capturing method and equipment | |
| CN115567163A (en) | A method, system and device for blind detection of frequency hopping signals | |
| CN107911328B (en) | Frame synchronization judging device and judging method | |
| CN117938199A (en) | Spread spectrum signal three-level capturing method, system, equipment and medium | |
| EP1109325B1 (en) | Apparatuses and method for acquisition of asynchronous wideband CDMA signals | |
| JP2003517241A (en) | Multi-bit spread spectrum signaling | |
| CN113225104A (en) | Method and system for capturing multi-user burst spread spectrum signal in asynchronous communication system | |
| CN110830077A (en) | A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals | |
| Wang et al. | Research on acquisition method of dsss system based on pmf-fft | |
| JP3639839B2 (en) | Spread spectrum communication device and high-speed synchronization establishment method thereof | |
| CN116719063A (en) | A method and system for rapid acquisition of ground-based navigation signals | |
| CN109150234B (en) | Direct sequence spread spectrum signal transmission method based on direct differential coherent accumulation | |
| JP4335913B2 (en) | Method and system for capturing a received impulse radio signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |



























