Disclosure of Invention
In order to solve the defects in the prior art, the application provides a high dynamic frequency acquisition tracking method of a low-orbit satellite, which adopts a chirp modulation mode, accurately determines the frequency deviation and timing deviation of a current chirp signal by processing the chirp signal, tracks the chirp signal and dynamically compensates the frequency deviation caused by the Doppler frequency change rate.
The technical effect to be achieved by the application is realized through the following scheme:
according to a first aspect of the present invention, there is provided a method for capturing and tracking a high dynamic frequency of a low-orbit satellite, comprising the steps of:
step 1: pre-compensating the input signal received by the receiving end;
step 2: performing mediation calculation on the compensated input signal and a local chirp signal and performing fast Fourier transform calculation;
step 3: calculating integer frequency offset according to the peak value of the fast Fourier transform, and calculating decimal frequency offset by using an interpolation method;
step 4: calculating frequency deviation and timing deviation according to the integer frequency deviation and the decimal frequency deviation;
step 5: calculating Doppler frequency deviation and Doppler change rate by utilizing the frequency deviation and the timing deviation, and feeding the Doppler frequency deviation and the Doppler change rate back to an input signal for compensation; the compensated input signal continues to step 2 to continue the capturing and tracking of the signal with high dynamics.
Preferably, in step 1, the mode of precompensation is: and calibrating the temperature compensation crystal oscillator or the external crystal oscillator according to the accurate clock to obtain the initial frequency offset of the crystal oscillator of the receiving end, and pre-compensating the input signal according to the initial frequency offset.
Preferably, in step 2, a local chirp signal is generated from the frame structure of the input signal, with a local up chirp signal or a local down chirp signal corresponding to the input signal.
Preferably, in the frame structure of the local chirp signal:
defining a local upchirp signal:
,
;
defining a local downchirp signal:
,
;
where BW is the bandwidth of the signal, Ts Is the time of one symbol, e is the euler number, j is the imaginary unit.
Preferably, the manner of mediating calculation and performing the fast fourier transform calculation in step 2 is:
multiplying the received upchirp signal with a local downchirp signal to obtain a direct current signal, and setting the position of the maximum value after N points of fast Fourier transformation at the position of the 1 st point;
according to the formula
To obtain the received signal r
0 (t);
Wherein s is
0 (t) is the original signal of the system,
e is Euler number, j is imaginary number unit;
the sampled and demodulated signal is:
,
wherein: b is the signal sampling bandwidth;
after fast fourier transformation, it is obtained that:
preferably, in step 3, the formula after the fast fourier transform is obtained:
the integer frequency offset is:
;
fractional frequency offset is calculated by using a spectral line difference method according to the formula:
for the position where the maximum value is in the N-point fast Fourier transform is k
0 Corresponding amplitude is X
k0 Taking k
0 The position on the right side is k
0+1 Amplitude is X
k0+1 Taking k
0 The position on the left side is k
0-1 Amplitude is X
k0-1 Due to |X
k0 | >| X
k0+a I, therefore->
Between->
And
wherein a= ±1 represents a direction; demodulating the downchirp signal to obtain +.>
ThenFrequency deviation->
Timing deviation->
。
Preferably, in step 5, when the signal is affected by the doppler rate of change, the input signal is expressed according to the following formula:
wherein:
for Doppler frequency offset +.>
For Doppler rate of change, +.>
The sampling period of the signal is represented by θ, and the phase is represented by θ;
multiplication of the input signal with the local chirp signal yields:
wherein the method comprises the steps of
For noise Gaussian random signal, the maximum likelihood criterion is used to estimate +.>
Of 4 typical phases; after estimating the minimum mean square error, the solution of the equation can be much more based on the minimum mean square error criterionThe rate of change of the puler.
Preferably, the 4 typical phases are:
the estimated minimum mean square error is:
wherein:
is n
i The values are N/8, 3N/8, 5N/8 and the phases of the chirp signals corresponding to the positions of 7N/8;
the estimated value of the Doppler change rate is:
the estimated value of the Doppler change rate is:
。
according to a second aspect of the present invention, there is provided a server comprising: a memory and at least one processor;
the memory stores a computer program, and the at least one processor executes the computer program stored by the memory to implement the low-orbit satellite high-dynamic frequency acquisition tracking method described in any one of the above.
According to a third aspect of the present invention, there is provided a computer readable storage medium having stored therein a computer program which when executed implements the low orbit satellite high dynamic frequency acquisition tracking method of any of the above.
According to one embodiment of the invention, the beneficial effects of the invention are: the method provides a channel framing method, which is used for frequency synchronization of burst channels; the method for estimating and tracking the frequency offset can be used for carrying out accurate frequency offset calculation on the signals in a high dynamic scene; in the communication process, firstly, a pre-calibration mode is used for estimating the frequency offset brought by a terminal due to a transmitting channel, a receiving channel and a crystal oscillator, and then the frequency offset introduced by the terminal is calibrated; in the signal receiving process, the information carried in the signal frame synchronizing signal is facilitated, and the initial frequency offset and the frequency change rate are estimated and compensated. The frequency change rate in the communication frame is tracked by the pilot signal to compensate for the frequency error due to the Doppler change rate. The method and the device fundamentally solve the influence of Doppler frequency offset and Doppler change rate on the reception in low-orbit satellite communication, greatly improve the communication length limit of burst signals, and can transmit more information in one communication process.
Detailed Description
For the purposes, technical solutions and advantages of the present application, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
In an embodiment of the present application, a method for capturing and tracking a high dynamic frequency of a low-orbit satellite as shown in fig. 1 includes the following steps:
step 1: pre-compensating the input signal received by the receiving end;
in this step, the mode of precompensation is: and (3) calibrating the temperature compensation crystal oscillator or the external crystal oscillator according to the accurate clock to obtain the initial frequency offset of the crystal oscillator of the receiving end, and pre-compensating the input signal according to the initial frequency offset to prevent the influence caused by the time drift of the crystal oscillator.
Step 2: performing mediation calculation on the compensated input signal and a local chirp signal and performing fast Fourier transform calculation;
the local chirp signal in this step is generated according to the frame structure of the input signal, and the frequency synchronization frame is as shown in fig. 2, and has a local up chirp signal or a local down chirp signal corresponding to the input signal; the frame head comprises 99 upchirp signals and 1 downchirp signal, the synchronization head comprises 1 upchirp signal and 1 downchirp signal, and the intervals between the frame head and the synchronization head and between the synchronization heads are 10ms.
In the frame structure of the local chirp signal:
defining a local upchirp signal:
,
;
defining a local downchirp signal:
,
;
where BW is the bandwidth of the signal,
is the time of one symbol, e is the euler number, j is the imaginary unit.
The mode of mediating calculation and performing fast Fourier transform calculation is as follows:
multiplying the received upchirp signal with a local downchirp signal to obtain a direct current signal, and setting the position of the maximum value after N points of fast Fourier transformation at the position of the 1 st point;
according to the formula
To obtain the received signal r
0 (t);
Wherein s is
0 (t) is the original signal of the system,
for frequency offset, e is the Euler number and j is the imaginary unit.
The sampled and demodulated signal is:
,
wherein: b is the signal sampling bandwidth;
after fast fourier transformation, it is obtained that:
step 3: calculating integer frequency offset according to the peak value of the fast Fourier transform, and calculating decimal frequency offset by using an interpolation method;
specifically, the method is obtained according to the formula after the fast Fourier transform:
the maximum point occurs at:
representing taking the nearest integer as the integer frequency offset;
the resolution due to the N-point FFT (fast fourier transform) is:
for more accurate determination of frequency offset, for N-point FFT, the position where the maximum value is located is k
0 Corresponding FFT amplitude is X
k0 Taking k
0 The position on the right side is k
0+1 Amplitude is X
k0+1 Taking k
0 The position on the left side is k
0-1 Amplitude is X
k0-1 Calculation by adopting spectral line interpolation method
Wherein->
= ±1, indicating direction, if |x
k0+1 | >| X
k0-1 I, then->
=1, otherwise, ++>
= -1 due to |x
k0 | >| X
k0+a I, therefore->
Between->
And->
Between them; and the estimation of the integer frequency offset and the decimal frequency offset is completed.
Step 4: calculating frequency deviation and timing deviation according to the integer frequency deviation and the decimal frequency deviation;
according to the upper partThe integer frequency offset and the decimal frequency offset in the steps are used for demodulating the downlink chirp signal, and the same way is adopted to obtain
Frequency deviation->
Timing deviation->
The timing deviation parameter represents the deviation between the sampling time of the receiving end and the sending end, and the calculation compensation can be performed through the deviation so as to ensure the accuracy of the subsequent input signal compensation.
Step 5: calculating frequency deviation and timing deviation to obtain frequency deviation and frequency change rate, and feeding the frequency deviation and the frequency change rate back to an input signal for compensation;
when the signal is affected by the doppler rate of change, the input signal is expressed according to the following formula:
wherein:
for Doppler frequency offset +.>
For Doppler rate of change, +.>
The sampling period of the signal is represented by θ, and the phase is represented by θ;
multiplication of the input signal with the local chirp signal yields:
wherein the method comprises the steps of
For noise Gaussian random signal, the maximum likelihood criterion is used to estimate +.>
The 4 typical phases of (2) are:
the estimated minimum mean square error is:
wherein:
is n
i The values are N/8, 3N/8, 5N/8 and the phases of the chirp signals corresponding to the positions of 7N/8;
according to the minimum mean square error criterion, solving the equation to obtain the estimated value of the Doppler change rate is as follows:
and (3) locally making corresponding compensation signals for the estimated Doppler frequency offset and Doppler change rate, and feeding the compensation signals back to the input signals for compensation.
The compensated input signal continues to step 2 to continue the capturing and tracking of the input signal, so that the capturing and tracking process of the signal under high dynamic state can be completed, as shown in fig. 3.
According to a second aspect of the present invention, as shown in fig. 4, there is provided a server comprising: amemory 401 and at least oneprocessor 402;
thememory 401 stores a computer program, and the at least oneprocessor 402 executes the computer program stored in thememory 401 to implement the low-orbit satellite high dynamic frequency acquisition tracking method described above.
According to a third aspect of the present invention, there is provided a computer readable storage medium having stored therein a computer program which when executed implements the low-orbit satellite high dynamic frequency acquisition tracking method described above.
According to one embodiment of the invention, the beneficial effects of the invention are: the method provides a channel framing method, which is used for frequency synchronization of burst channels; the method for estimating and tracking the frequency offset can be used for carrying out accurate frequency offset calculation on the signals in a high dynamic scene; in the communication process, firstly, a pre-calibration mode is used for estimating the frequency offset brought by a terminal due to a transmitting channel, a receiving channel and a crystal oscillator, and then the frequency offset introduced by the terminal is calibrated; in the signal receiving process, the information carried in the signal frame synchronizing signal is facilitated, and the initial frequency offset and the frequency change rate are estimated and compensated. The frequency change rate in the communication frame is tracked by the pilot signal to compensate for the frequency error due to the Doppler change rate. The method and the device fundamentally solve the influence of Doppler frequency offset and Doppler change rate on the reception in low-orbit satellite communication, greatly improve the communication length limit of burst signals, and can transmit more information in one communication process.
It should be noted that the foregoing detailed description is exemplary and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments in accordance with the present application. As used herein, the singular is intended to include the plural unless the context clearly indicates otherwise. Furthermore, it will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, steps, operations, devices, components, and/or groups thereof.
It should be noted that the terms "first," "second," and the like in the description and claims of the present application and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are capable of operation in sequences other than those illustrated or otherwise described herein.
Furthermore, the terms "comprise" and "have," as well as any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those elements but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways, such as rotated 90 degrees or at other orientations, and the spatially relative descriptors used herein interpreted accordingly.
In the above detailed description, reference is made to the accompanying drawings, which form a part hereof. In the drawings, like numerals typically identify like components unless context indicates otherwise. The illustrated embodiments described in the detailed description, drawings, and claims are not meant to be limiting. Other embodiments may be utilized, and other changes may be made, without departing from the spirit or scope of the subject matter presented herein.
The above is only a preferred embodiment of the present invention, and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.