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 PDF

Info

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
Application number
CN201910913971.8A
Other languages
Chinese (zh)
Other versions
CN110830077B (en
Inventor
蒙艳松
朱向鹏
王延光
王登峰
张中英
踪念科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201910913971.8A priority Critical patent/CN110830077B/en
Publication of CN110830077A publication Critical patent/CN110830077A/en
Application granted granted Critical
Publication of CN110830077B publication Critical patent/CN110830077B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/06Receivers
    • H04B1/16Circuits

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

一种改善多路突发信号接收性能的快速捕获方法,用于解决报文通信用户使用特殊伪码进行信号调制的捕获问题。该方法首先进行零中频数据与伪码的采集,然后将零中频数据与伪码进行划分块处理,利用双块拓展方法完成零中频数据块与伪码块的拓展;根据零中频数据与主扩频伪码相关计算结果的排列规律,完成相位并行搜索与频率并行搜索,进而完成所有伪码相位的捕获。本发明捕获方法极大地缩短了捕获的搜索时间,同时捕获性能优异,使这种特殊的伪码结构能够应用于实际工程。

Figure 201910913971

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.

Figure 201910913971

Description

一种改善多路突发信号接收性能的快速捕获方法A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals

技术领域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:

Figure BDA0002215532220000032
Figure BDA0002215532220000032

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,

Figure BDA0002215532220000033
表示rk(i)拓展后的零中频数据块,拓展方法如下所示:The extended length of the zero-IF data block is 8Q,
Figure BDA0002215532220000033
Represents the zero-IF data block after expansion of r k (i). The expansion method is as follows:

Figure BDA0002215532220000041
Figure BDA0002215532220000041

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计算得到噪声方差

Figure BDA0002215532220000043
则仅有噪声时,捕获检测量I符合标准的中心
Figure BDA0002215532220000044
分布;若信号存在,则捕获检测量I服从非中心
Figure BDA0002215532220000045
分布,非中心
Figure BDA0002215532220000046
分布的参数λ=2X*snr,X等于1,snr为相干累加后的信噪比。In the presence of only noise, the noise variance is calculated according to the captured detection quantity I
Figure BDA0002215532220000043
Then when there is only noise, the capture detection quantity I meets the center of the standard
Figure BDA0002215532220000044
distribution; if the signal exists, the capture detection quantity I obeys the non-centrality
Figure BDA0002215532220000045
distribution, non-central
Figure BDA0002215532220000046
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:

Figure BDA0002215532220000047
Figure BDA0002215532220000047

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:

Figure BDA0002215532220000061
Figure BDA0002215532220000061

上式中,tj为第j个采样时刻,本地接收信号的采样速率为伪码速率的两倍;τ为信号到达报文通信接收机的时延;A为信号幅度,fL1为射频载波频率,fd为多普勒频移,fIF为中频载波频率,

Figure BDA0002215532220000062
为初始相位,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,
Figure BDA0002215532220000062
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,

Figure BDA0002215532220000071
表示sk(i)拓展后的伪码块,拓展方法如公式所示:The length of the zero-IF data block and the pseudo-code block after double-block expansion is 8Q,
Figure BDA0002215532220000071
represents the expanded pseudo-code block of s k (i), and the expansion method is shown in the formula:

Figure BDA0002215532220000082
表示rk(i)拓展后的零中频数块,拓展方法如公式所示
Figure BDA0002215532220000082
Represents the zero-IF block after expansion of r k (i), the expansion method is shown in the formula

Figure BDA0002215532220000083
Figure BDA0002215532220000083

(3)相位并行搜索:对拓展后的每块零中频数据和每块伪码分别进行FFT运算,将第k块零中频数据与第k块伪码的FFT运算结果进行共轭相乘运算,然后对运算结果进行IFFT运算,得到零中频数据与伪码的相关计算结果

Figure BDA0002215532220000084
存储在搜索区中,
Figure BDA0002215532220000085
为主扩频码每个相位的位置。(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
Figure BDA0002215532220000084
stored in the search area,
Figure BDA0002215532220000085
is the position of each phase of the main spreading code.

Figure BDA0002215532220000086
Figure BDA0002215532220000086

按照上述步骤与方法,完成(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的

Figure BDA0002215532220000092
分布。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;
Figure BDA0002215532220000092
distributed.

在仅有噪声存在下,根据捕获检测量I计算得到噪声方差

Figure BDA0002215532220000093
则仅有噪声时I符合标准的中心
Figure BDA0002215532220000094
分布,若信号存在,则I服从非中心
Figure BDA0002215532220000095
分布,非中心分布参数为λ=2X*snr(进行了
Figure BDA0002215532220000096
的归一化,X等于1,snr为相干累加后的信噪比),可以得到自由度为2的中心
Figure BDA0002215532220000097
分布与非中心
Figure BDA0002215532220000098
分布的概率分布,根据检测估值理论与上述计算关系,捕获门限Vt计算如下,虚警概率Pfa为已知量。In the presence of only noise, the noise variance is calculated according to the captured detection quantity I
Figure BDA0002215532220000093
Then I only meet the center of the standard when there is only noise
Figure BDA0002215532220000094
distribution, if the signal exists, then I obey the non-centrality
Figure BDA0002215532220000095
distribution, the non-central distribution parameter is λ=2X*snr (performed
Figure BDA0002215532220000096
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
Figure BDA0002215532220000097
Distribution and non-centrality
Figure BDA0002215532220000098
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)

1.一种改善多路突发信号接收性能的快速捕获方法,其特征在于,包括如下步骤:1. a fast acquisition method for improving multi-channel burst signal receiving performance, is characterized in that, comprises the 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. 2.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(1)中,设零中频数据采集时间为T毫秒,则:2. a kind of fast acquisition method that improves multi-channel burst signal reception performance according to claim 1, is characterized in that, in described step (1), set zero intermediate frequency data acquisition time to be T millisecond, 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. . 3.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(2)中,对伪码进行划分块处理的方法如下:3. a kind of fast acquisition method improving multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (2), the method that pseudocode is divided into blocks is processed 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. 4.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(2)中,对零中频数据进行划分块处理的方法如下:4. a kind of fast acquisition method that improves multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (2), the method that the zero-IF data is divided into blocks and processed 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. 5.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(2)中,按照如下方法完成伪码块的拓展:5. a kind of fast acquisition method improving multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (2), completes the expansion of pseudo code block according to the following method: 伪码块拓展后的长度为8Q,s′k(i)表示sk(i)拓展后的伪码块,拓展方法如下所示:The length of the expanded pseudo-code block is 8Q, and s' k (i) represents the expanded pseudo-code block of s k (i). The expansion method is as follows:
Figure FDA0002215532210000021
Figure FDA0002215532210000021
sk(i)为第k个伪码块内第i个半码片。s k (i) is the i-th half-chip in the k-th pseudo-code block.
6.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(2)中,按照如下方法完成零中频数据块的拓展:6. a kind of fast acquisition method improving multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (2), completes the expansion of zero intermediate frequency data block according to the following method: 零中频数据块拓展后的长度为8Q,r′k(i)表示rk(i)拓展后的零中频数据块,拓展方法如下所示:The extended length of the zero-IF data block is 8Q, and r' k (i) represents the zero-IF data block extended by r k (i). The extension method is as follows:
Figure FDA0002215532210000031
Figure FDA0002215532210000031
rk(i)为第k个零中频数据块内的第i个数据。r k (i) is the i-th data in the k-th zero-IF data block.
7.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(4)中,对每一个伪码相位对应的频率进行搜索的实现方法如下:7. a kind of fast acquisition method for improving multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (4), the realization method of searching the frequency corresponding to each pseudo code phase 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. 8.根据权利要求7所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,FFT计算点数最好满足2的n次方,n为整数,满足2n-X2≥0且取最小值,如果X2小于2n,在相同伪码相位计算得到的相关值后,添加2n-X2个零进行FFT计算;如果X2等于2n,直接进行FFT计算。8. The fast acquisition method for improving multi-channel burst signal receiving performance according to claim 7, wherein, the number of FFT calculation points preferably satisfies the nth power of 2, n is an integer, and 2n-X2≥ 0 and take the minimum value, if X2 is less than 2 n , add 2 n -X2 zeros to perform FFT calculation after the correlation value calculated by the same pseudo code phase; if X2 is equal to 2 n , directly perform FFT calculation. 9.根据权利要求1所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(5)中根据FFT运算输出的结果行能量计算的实现方式如下:9. a kind of fast acquisition method improving multi-channel burst signal receiving performance according to claim 1, is characterized in that, in described step (5), the realization mode of energy calculation according to the result row of 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计算得到噪声方差
Figure FDA0002215532210000041
则仅有噪声时,捕获检测量I符合标准的中心χ2分布;若信号存在,则捕获检测量I服从非中心χ2分布,非中心χ2分布的参数λ=2X*snr,X等于1,snr为相干累加后的信噪比。
In the presence of only noise, the noise variance is calculated according to the captured detection quantity I
Figure FDA0002215532210000041
Then when there is only noise, the captured detection quantity I conforms to the standard central χ 2 distribution; if the signal exists, the captured detection quantity I obeys the non-central χ 2 distribution, and the parameter λ=2X*snr of the non-central χ 2 distribution, X is equal to 1 , snr is the signal-to-noise ratio after coherent accumulation.
10.根据权利要求9所述的一种改善多路突发信号接收性能的快速捕获方法,其特征在于,所述步骤(5)中,检测门限Vt计算如下:10. The fast acquisition method for improving multi-channel burst signal reception performance according to claim 9, characterized in that, in the step (5), the detection threshold V t is calculated as follows:
Figure FDA0002215532210000042
Figure FDA0002215532210000042
Pfa为虚警概率。 Pfa is the false alarm probability.
CN201910913971.8A 2019-09-25 2019-09-25 A Fast Acquisition Method to Improve the Receiving Performance of Multi-channel Burst Signals Active CN110830077B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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