Signal capturing method for high dynamic satellite communication system
Technical Field
The invention relates to the field of communication, in particular to a signal capturing method of a high-dynamic satellite communication system.
Background
With the rapid development of low-orbit satellite applications and the communication demands of airborne and missile-borne platforms, the capturing of high-dynamic signals is an important problem to be solved in satellite communication. If the conventional methods of matched filtering search, linear search, code phase search, parallel frequency search and the like are continuously used for capturing the high dynamic signals, in order to adapt to large doppler shift, the methods need to repeatedly call, time-sharing or capture for multiple times at the same time, which consumes more hardware resources and calculation time, prolongs the whole synchronization time of the system and affects the communication performance of the receiving terminal.
Therefore, there is a need for a fast acquisition method for a high dynamic satellite communication system.
Disclosure of Invention
In view of this, the present invention provides a signal acquisition method for a high dynamic satellite communication system. The method not only can realize the capturing of the high dynamic satellite signals, but also can calculate the Doppler frequency offset and the code offset of the satellite, and solves the problem of difficult capturing of the high dynamic signals. Compared with the traditional capturing method, the method has the advantages of wide frequency offset range and high capturing performance, reduces the length of the shift register by using a parallel mode, reduces the problem of more occupied resources caused by FFT, and can well realize the capturing function for most of high-dynamic platforms, such as low-orbit satellites, airborne platforms and missile-borne platforms.
The purpose of the invention is realized in the following way:
a method for signal acquisition in a high dynamic satellite communication system, comprising the steps of:
(1) Receiving a 4-times oversampled signal r (k);
(2) The shunt is performed on r (k), dividing r (k) into 4 parallel paths of parallel signals;
(3) Respectively carrying out sliding storage on 4 paths of parallel signals, wherein the storage length is the frame head length;
(4) Dividing the 4 paths of stored signals into a plurality of equal-length segments, and correspondingly, dividing the local code word into the same number of equal-length segments;
(5) Performing correlation operation on the received signal and the local code word in each segment to obtain a correlation operation result y (k) of each segment, and storing the correlation operation result of each segment;
(6) Zero padding is carried out at the tail part of Y (k), then N-point FFT operation is carried out, Y (k) is output, 4 paths of parallel signals can generate 4 paths of Y (k), and the Y (k) is respectively marked as Y 1(k)、Y2(k)、Y3(k)、Y4 (k);
(7) Respectively calculating the amplitude values of Y 1(k)、Y2(k)、Y3(k)、Y4 (k), and sliding and accumulating a plurality of adjacent amplitude values of each path;
(8) Searching the accumulated maximum amplitude value of each path, comparing the maximum amplitude values of the 4 paths, and taking the maximum amplitude value P max in the 4 paths as a final correlation result;
(9) Comparing the maximum amplitude P max with a threshold, if the maximum amplitude P max is larger than the threshold, the capturing is considered to be successful, if the maximum amplitude P max is not larger than the threshold, the frame head is considered not captured, and the steps (3) - (8) are repeated, and sliding storage and FFT operation are continued until the capturing is successful;
(10) After successful acquisition, the Doppler frequency offset and the code offset of the tracked high-dynamic satellite signals are determined according to the position k of the maximum value.
Further, in the step (3), a shift register is used for storing, the length of the shift register is the length of the frame header, and the content of the memory slides once for each symbol clock, namely, each symbol clock removes the oldest data and enters the latest data.
Further, in the step (4), the length of each segment is greater than 1 and less than or equal to M m Is (1) x, wherein:
where R s is the symbol rate and f d is the Doppler frequency offset to be adapted.
Further, in the step (10), the method for calculating the doppler frequency offset f d of the high dynamic satellite signal is as follows:
wherein N is the number of FFT points, M is the length of each segment in step (4), k is the position of the maximum amplitude in the N-point FFT, and R s is the symbol rate.
Further, in step (10), the code offset δ is calculated according to the doppler frequency offset f d:
Wherein f 0 is the carrier frequency.
Compared with the background technology, the invention has the following advantages:
1. compared with the traditional capturing method, the method has the advantages of wide adaptive frequency offset range and high capturing performance.
2. The invention adopts a parallel mode, reduces the length of the shift register, reduces the complexity of operation and reduces the problem of more occupied resources caused by FFT.
3. The method and the device take the maximum value after the amplitude sliding accumulation, reduce the influence of spectrum leakage on the capturing performance, and improve the adaptability to the condition of low signal to noise ratio.
4. The invention is suitable for most high dynamic platforms, such as low-orbit satellites, airborne platforms, missile-borne platforms and the like, and can well realize the capturing function.
Drawings
FIG. 1 is a block diagram of the components of the method of the present invention.
Fig. 2 is a schematic diagram of the segmentation method of step (4).
FIG. 3 is a graph of capture performance of the method of the present invention.
Detailed Description
The present invention will be further described below.
As shown in fig. 1, a signal capturing method of a high dynamic satellite communication system includes the following steps:
(1) A 4 times oversampled signal r (k) is received, r (k) typically being a matched filtered signal. Taking symbol rate R s as 2Msps as an example, it can be seen that the sampling rate of R (k) is 8MHz;
(2) And (3) branching r (k). Dividing r (k) into 4 parallel signals, wherein each signal is equivalent to 1 time of sampling data, namely the speed is 2M Hz;
(3) The 4 paths of signals are respectively stored in a sliding way according to a certain length, wherein the storage length is the frame head length, and the frame head length is calculated by taking the frame head length as 128 as an example.
A shift register may be used for storage, the shift register length being the frame header length 128. The memory content is slid once per symbol clock, i.e. each symbol clock removes one oldest data while entering one newest.
(4) As shown in fig. 2, the 4 paths of stored signals are respectively divided into a plurality of equal-length segments, and the local code word is correspondingly divided into the same number of equal-length segments;
In practical applications, the length M is determined according to the symbol rate R s and the doppler frequency offset f d.
If the doppler frequency offset to be adapted is f d, the length of each segment must be less than or equal to M m Is (1) x.
Calculated by using Doppler frequency offset f d as-220 kHz, the method can be obtained:
I.e. at most 4 symbols per segment length, for convenience of implementation on a hardware FPGA, 4 symbols per segment length may be selected, and there are 128/4=32 segments in total.
(5) And performing correlation operation on the received signal and the local code word in each segment to obtain a correlation operation result y (k) of each segment, and storing the correlation operation result of each segment, wherein 32 correlation results can be obtained in the example.
(6) And (3) supplementing zeros with a certain length at the tail part of Y (k), then sending the zeros into an FFT processing unit for N-point FFT operation, and outputting Y (k), wherein the length of data after zero supplementation is generally 2 N, and N is an integer. The more the zero padding number is, the estimation accuracy of Doppler frequency offset is generally improved, but the estimation accuracy is not increased infinitely. The problem is that the operand is increased and the occupied resources are increased.
In engineering application, less zeros are added so that the total data length can meet 2≡n, and the zero added number can be determined according to the requirements and conditions of estimation precision, signal-to-noise ratio and the like.
To improve accuracy, a 512-point FFT may be performed, with 480 zeros being padded. The estimation accuracy is 2000/4/512=0.98 kHz.
(7) The 4-path parallel signal generates 4 paths Y (k), which are respectively marked as Y 1(k)、Y2(k)、Y3(k)、Y4 (k);
(8) Calculating the amplitude values of Y 1(k)、Y2(k)、Y3(k)、Y4 (k) respectively, and sliding and accumulating 3 or more adjacent amplitude values of each path;
Typically, the number of accumulated magnitudes is no more than 5. In this example, 3 were taken.
(9) Searching the accumulated maximum amplitude value of each path, comparing the maximum amplitude values of the 4 paths, and taking the maximum amplitude value P max in the 4 paths as a final correlation result;
(10) And if the maximum amplitude P max is not greater than the threshold, no frame header is captured, continuing to repeat the steps 3-9, and continuing to perform sliding storage and FFT operation until the capturing is successful.
The magnitude of the simulated output is shown in figure 3 when the capture is successful. It can be seen from the figure that the position of the maximum is 287 th point, i.e. the position k of the maximum is 287.
(11) After successful acquisition, the Doppler frequency offset and the code offset of the tracked high-dynamic satellite signals are determined according to the position k of the maximum value.
The Doppler frequency offset f d of the high dynamic satellite signal is calculated as follows:
Wherein N is the number of FFT points, M is the length of each segment in step (4), k is the kth position where the maximum amplitude is in the N-point FFT, and R s is the symbol rate.
Converting k=287 to signed number, substitutingObtaining:
Is very close to the true value of the doppler frequency offset f d, only 0.3kHz difference, the frequency offset estimate may be considered correct.
The code bias delta and the Doppler bias f d are in a direct proportion relation, and the estimated f d is substituted intoObtaining Doppler code bias. Wherein f 0 is the carrier frequency.
In a word, the method comprises the steps of parallel branching, partial correlation, zero padding FFT, searching for the maximum amplitude, accumulating a plurality of adjacent amplitude values, judging, calculating Doppler frequency offset and code offset and the like, so that the method can realize the capturing of the high dynamic satellite signals, calculate the Doppler frequency offset and the code offset of the satellites at the same time, and solve the problem of difficult capturing of the high dynamic signals. Compared with the traditional capturing method, the method has the advantages of wide frequency offset range, high capturing performance and lower signal to noise ratio, and reduces the length of the shift register, the operation complexity and the occupation of more resources caused by FFT by using a parallel mode. The method is suitable for most of high-dynamic platforms, such as low-orbit satellites, airborne platforms, missile-borne platforms and the like, can well realize the capturing function, and is particularly suitable for signal detection of a communication system.
The above description is only one embodiment of the present invention, but the scope of the present invention is not limited thereto. Any equivalent substitutions or changes made by those skilled in the art according to the technical scheme and the inventive concept of the present invention should be included in the scope of the present invention.