CN114924121A - Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal - Google Patents
Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal Download PDFInfo
- Publication number
- CN114924121A CN114924121A CN202210787392.5A CN202210787392A CN114924121A CN 114924121 A CN114924121 A CN 114924121A CN 202210787392 A CN202210787392 A CN 202210787392A CN 114924121 A CN114924121 A CN 114924121A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- sub
- change rate
- frequency change
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Abstract
The invention provides a frequency spectrum detection and estimation method for an extraterrestrial celestial body detector signal, which comprises the following steps of: estimating the frequency range and the frequency change rate range of a detector signal to be acquired by a ground station in advance; configuring signal acquisition parameters of a ground station; continuously dividing the frequency bandwidth of the signal acquisition unit into a plurality of sub-frequency bands at equal intervals, and continuously dividing the frequency change rate range into a plurality of frequency change rate estimation value sub-intervals at equal intervals; continuously dividing the frequency bandwidth of each sub-frequency band into a plurality of sub-frequency intervals at equal intervals; and the signal corresponding to the searched power peak value is the target signal. According to the invention, the terrestrial celestial body acquisition signals received by the ground antenna are divided into a plurality of sub-frequency bands, a plurality of sub-frequency intervals and a plurality of frequency change rate estimation value sub-intervals according to the frequency band range, so that rapid signal analysis processing is carried out, the calculation resources are reduced, and rapid signal analysis processing is realized.
Description
Technical Field
The invention belongs to the technical field of spacecraft measurement and control, and particularly relates to a frequency spectrum detection and estimation method for signals of an extraterrestrial celestial body detector.
Background
The flying distance of the extraterrestrial celestial body detector is longer and longer, and the strength of a downlink signal of the extraterrestrial celestial body detector is weaker and weaker. At present, the farthest distance from a Mars detector to the ground reaches 4 hundred million kilometers, the distances from Planet detectors such as Jupiter, Saturn, Tianwang and the like to the ground are farther, the square of the distance from the detector to the ground is in direct proportion to the attenuation degree of signals, and therefore the strength of downlink signals of the detector is obviously weakened along with the increase of the distance. In addition, for the planetary detection task, the detector will go through the process of planetary capture, part of the detector will also realize landing on the planet, and the speed of the detector changes rapidly, thereby causing the dynamics of the downstream signal to be very large.
In the traditional mode, a ground station performs closed-loop signal estimation on a detector signal by adopting a phase-locked loop technology to complete detection on the detector signal. When the detector signal strength is relatively strong and the signal dynamic is relatively small, the closed-loop detection of the detector signal can be effectively finished. However, for the large dynamic and weak signals of the detector, the closed-loop signal estimation method cannot be well adapted, the detection and estimation difficulty of the detector signal is large, and the probability of signal lock losing is obviously increased.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a frequency spectrum detection and estimation method for the signal of the extraterrestrial celestial body detector, which can effectively solve the problems.
The technical scheme adopted by the invention is as follows:
the invention provides a frequency spectrum detection and estimation method for an extraterrestrial celestial body detector signal, which comprises the following steps of:
step 4, dividing the original acquisition signal into a plurality of sections of acquisition signal units by taking the time length T as a period;
for each section of the acquired signal unit d T All execute step 5-step 8, search and collect the signal unit d T Target signal of (1):
step 5, collecting signal unit d T The frequency bandwidth of (a) is continuously divided into M sub-bands at equal intervals, and the acquisition signal of the M sub-band is represented as: d is a radical of T (m),m＝1，2，3，…，M；
Step 6, continuously dividing the frequency change rate range estimated in the step 1 into N frequency change rate estimation value subintervals at equal intervals;
acquisition signal d for m-th sub-band T (m) respectively calculating the frequency and frequency change rate two-dimensional compensation value of each frequency change rate estimation value subinterval, thereby obtaining N two-dimensional compensated acquisition signals g T (m, N), wherein N is 1,2,3., N;
step 7, continuously dividing the frequency bandwidth of the mth sub-band into Q sub-frequency intervals at equal intervals;
collecting signal g after two-dimensional compensation of nth frequency change rate estimation value subinterval of mth subband T (m, n), respectively calculating the Fourier transform result of each sub-frequency interval, thereby obtaining Q signals after Fourier transform, namely the three-dimensional acquisition signal G T (m,n,q)；
Step 8, therefore, for the acquisition signal unit d T Performing fine search calculation on the M sub-frequency bands, the N frequency change rate estimation value sub-intervals and the Q sub-frequency intervals to obtain M, N and Q three-dimensional acquisition signals G T (m,n,q)；
Comparing M N Q three-dimensional collected signals G T The power of (m, n, q), the sub-band number corresponding to the power peak, the sub-interval number of the estimated frequency change rate and the sub-frequency interval number are the searched targetThe position of the target signal is searched for a target signal;
and 9, estimating a frequency spectrum detection estimation value according to the searched target signal.
Preferably, in step 1, the frequency range of the detector signal is denoted as f min ,f max ]The frequency change rate range is represented by [ f' min ,f′ max ]：
Wherein:
f min the estimated value of the minimum frequency value of the detector signal to be acquired is obtained;
f max the frequency maximum value estimation value of the detector signal to be acquired is obtained;
f′ min the minimum value of the frequency change rate of the detector signal to be acquired is estimated;
f′ max the estimated value of the maximum frequency change rate of the detector signal to be acquired is obtained;
in step 2, the configured signal acquisition parameters of the ground station comprise the frequency bandwidth B of the ground station 0 And a center frequency f B Obtained by the following formula:
B 0 ＝f max -f min
preferably, in step 5, a signal collecting unit d is adopted T Frequency bandwidth of (a) and frequency bandwidth B of the ground station arrangement 0 Are equal.
Preferably, in step 5, the collected signal d of each sub-band T (m) is represented by:
wherein: n is a radical of an alkyl radical 0 S (f (t)) is noise and is a target signal;
the meaning is as follows:
of the M subbands, only one subband has the target signal, assuming that there is the target signal in the z-th subband, where z is an unknown value; the other M-1 sub-bands are noise;
the step 6 specifically comprises the following steps:
the frequency change rate range [ f 'estimated in the step 1' min ,f′ max ]Equally spaced and continuously divided into N frequency rate of change estimate subintervals, expressed as: e (N), N ═ 1,2,3, N; therefore, the width Δ f' ═ e (N)/N of each frequency rate estimation value subinterval;
acquisition signal d for m-th sub-band T (m) performing two-dimensional compensation of the frequency and the frequency change rate in the following manner to obtain a two-dimensionally compensated acquisition signal g T (m,n)：
Wherein:
j represents an imaginary unit;
f c (m) represents a lower boundary frequency compensation value of the mth sub-band;
f′ c (n) a lower boundary frequency change rate compensation value representing an nth frequency change rate estimate sub-interval;
wherein, f c (m) and f' c (n) is obtained using the formula:
f c (m)＝(m-1)Δf+f min
f′ c (n)＝(n-1)Δf′+f′ min
wherein: Δ f ═ B 0 and/M, representing the bandwidth of each sub-band.
Preferably, step 7 specifically comprises:
collecting signal g after two-dimensional compensation of nth frequency change rate estimation value subinterval of mth subband T (m, n) is converted into a three-dimensional acquisition signal G by adopting the following formula T (m,n,q)：
Wherein:
FFT represents fourier transform;
the meaning is as follows:
the frequency bandwidth Δ f of the mth sub-band is continuously divided into Q sub-frequency intervals at equal intervals, which is expressed as: fr (Q), Q ═ 1,2,.., Q; two-dimensionally compensated acquisition signal g for each sub-band T (m, n) performing Fourier transform in each sub-frequency interval to obtain a Fourier transformed signal, namely a three-dimensional acquisition signal G T (m,n,q)。
Preferably, in step 8, the signal G is acquired three-dimensionally T The power of (m, n, q) is represented as P T (m, n, q) obtained by:
wherein:
P[N 0 ]represents N 0 Meaning: power of the noise part after Fourier transform;
p [ S (m, n, q) ] represents the power of S (m, n, q), meaning: the target signal s (f (t)) is partially fourier-transformed.
Preferably, in step 8, the power peak is represented by P T (m T ,n T ,q T )；m T ,n T ,q T The meanings are respectively as follows: the subband number where the power peak is located, the frequency change rate estimation value subinterval number, and the subband interval number are expressed as:
(m T ,n T ,q T )＝max[P T (m,n,q)]| m∈[1,M],n∈[1,N],q∈[1,Q]
the step 9 specifically comprises the following steps:
the frequency spectrum detection estimation value of the target signal comprises the following steps: accurate frequency estimation of target signal
Accurate estimation value of frequency change rate
And accurate estimation value P/N of signal power noise spectral density ratio 0 ；
Obtaining a signal acquisition unit d by adopting the following formula T Accurate frequency estimation of intermediate target signal
And frequency rate of change accurate estimate
Obtaining a signal acquisition unit d by adopting the following formula T Accurate estimation value P/N of signal power noise spectral density ratio of intermediate target signal 0 ：
P/N 0 ＝P T (m T ,n T ,q T )/P[N 0 ]
P[N 0 ]＝N/(Δf/Q)
Obtaining a spectral curve G T (m T ,n T Q), wherein Q ═ 1, 2.., Q, includes amplitude frequency curves and phase frequency curves;
wherein:
spectral curve G T (m T ,n T Q) denotes the abscissa FR (q) and the ordinate G T (m T ,n T Q) the curve formed.
The frequency spectrum detection estimation method of the signal of the extraterrestrial celestial body detector provided by the invention has the following advantages:
according to the invention, the terrestrial celestial body acquisition signals received by the ground antenna are divided into a plurality of sub-frequency bands, a plurality of sub-frequency intervals and a plurality of frequency change rate estimation value sub-intervals according to the frequency band range, so that rapid analysis processing of the signals is carried out, the calculation resources are reduced, and rapid analysis processing of the signals is realized.
Drawings
FIG. 1 is a schematic flow chart of a method for estimating the frequency spectrum detection of a signal of an extraterrestrial celestial object detector provided by the present invention;
FIG. 2 is a schematic diagram of a frequency and frequency rate of change two-dimensional compensation grid provided by the present invention;
fig. 3 is a diagram of a spectrum result of a signal of a mars detector provided by the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantageous effects solved by the present invention more clearly apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The invention provides a frequency spectrum detection and estimation method for signals of an extraterrestrial celestial body detector, which comprises the steps of dividing a received acquisition signal frequency band into a plurality of sub-frequency bands, and dividing each sub-frequency band into a plurality of sub-frequency intervals; meanwhile, the estimated frequency change rate range is divided into a plurality of frequency change rate estimation value sub-intervals, so that rapid signal capture and tracking are carried out in three-dimensional spaces of the sub-frequency bands, the sub-frequency intervals and the frequency change rate estimation value sub-intervals, the frequency spectrum of the target signal is recovered, the estimation of the frequency, the frequency change rate and the signal power noise spectral density ratio of the target signal is completed, and the monitoring of the frequency spectrum of the weak signal of the extraterrestrial celestial body detector is realized. The method for processing the sub-frequency bands, the sub-frequency intervals and the frequency change rate estimation value sub-intervals can realize capture and monitoring of extremely weak signals, and is particularly suitable for monitoring and estimating the frequency spectrum of the extraterrestrial celestial body detector signals.
The invention provides a frequency spectrum detection and estimation method for a signal of an extraterrestrial celestial body detector, which comprises the following steps of:
in this step, the frequency range of the detector signal is denoted as f min ,f max ]The frequency change rate range is represented by [ f' min ,f′ max ]；
Wherein:
f min the estimated value of the minimum frequency value of the detector signal to be acquired is obtained;
f max the frequency maximum value estimation value of the detector signal to be acquired is obtained;
f′ min the minimum value of the frequency change rate of the detector signal to be acquired is estimated;
f′ max the frequency change rate maximum value estimation value of the detector signal to be acquired is obtained.
The position of the extraterrestrial celestial object cannot be accurately known, but initial estimation can be carried out according to prior information, so that the relative motion between the extraterrestrial celestial object (namely a detector) and the ground station can be initially determined, and the initial estimation of the frequency and the frequency change rate of the detector signal required to be acquired by the ground station is realized.
in this step, the two configured signal acquisition parameters of the ground station are respectively the frequency bandwidth B of the ground station 0 And a center frequency f B Obtained by the following formula:
B 0 ＝f max -f min
step 4, dividing the original acquisition signal into a plurality of sections of acquisition signal units by taking the time length T as a period;
for each section of the acquired signal unit d T All execute step 5-step 8, search the signal collecting unit d T Target signal of (1):
step 5, collecting signal unit d T Is continuously divided into M sub-bands at equal intervals, and the acquired signal of the M sub-band is represented as: d T (M), M ═ 1,2,3, …, M; wherein, a signal collecting unit d T Frequency bandwidth of (a) and frequency bandwidth B of the ground station arrangement 0 Are equal.
Acquired signal d for each sub-band T (m) is represented by:
wherein: n is a radical of an alkyl radical 0 S (f (t)) is noise and is a target signal;
the meaning is as follows:
of the M subbands, only one subband has the target signal, assuming that there is the target signal in the z-th subband, where z is an unknown value; the other M-1 sub-bands are noise.
Step 6, continuously dividing the frequency change rate range estimated in the step 1 into N frequency change rate estimation value subintervals at equal intervals;
specifically, the frequency change rate range [ f 'estimated in the step 1' min ,f′ max ]The equal-interval continuous division is carried out on N frequency change rate estimation value subintervals, and the frequency change rate estimation value subintervals are expressed as: e (N), N ═ 1,2,3, N; therefore, the width Δ f' of each frequency change rate estimation value subinterval is e (N)/N.
Acquisition signal d for mth subband T (m) respectively calculating the frequency and frequency change rate two-dimensional compensation value of each frequency change rate estimation value subinterval, thereby obtaining N two-dimensional compensated acquisition signals g T (m, N), wherein N is 1,2,3., N;
as a specific implementation manner, the m-th sub-band is collected as the signal d T (m) performing two-dimensional compensation of the frequency and the frequency change rate in the following manner to obtain a two-dimensionally compensated acquisition signal g T (m,n)：
Wherein:
j represents an imaginary unit;
f c (m) represents a lower boundary frequency compensation value of the mth sub-band;
f′ c (n) a lower boundary frequency change rate compensation value representing an nth frequency change rate estimate sub-interval;
wherein, f c (m) and f' c (n) is obtained using the formula:
f c (m)＝(m-1)Δf+f min
f′ c (n)＝(n-1)Δf′+f′ min
wherein: Δ f ═ B 0 and/M, representing the frequency bandwidth of each sub-band.
As a specific implementation mode, the frequency and frequency change rate two-dimensional compensation grids shown in FIG. 2 can be manufactured, and each grid stores corresponding f in advance c (m) and f' c Specifically, for the mth point (M is 1,2,3.., M) on the horizontal axis and the nth point (N is 1,2,3.., N) on the vertical axis, the two-dimensional compensation matrix of the frequency and the frequency change rate is C (M, N), and C (M, N) is (f) (M is 1,2,3.., M), and c (m),f c ′(n))。
for example, in the C (1,3) grid of FIG. 2, f is stored c (1) And f' c (3) The value of (c). Thus, FIG. 2 can be looked up directly to quickly obtain the desired f c (m) and f' c The value of (n).
Step 7, continuously dividing the frequency bandwidth of the mth sub-band into Q sub-frequency intervals at equal intervals;
collecting signal g after two-dimensional compensation of nth frequency change rate estimation value subinterval of mth subband T (m, n), respectively calculating the Fourier transform result of each sub-frequency interval, thereby obtaining Q signals after Fourier transform, namely the three-dimensional acquisition signal G T (m,n,q)；
In this step, the m-th sub-band is two-dimensionally compensated for the nth frequency change rate estimation value sub-interval T (m, n) is converted into a three-dimensional acquisition signal G by adopting the following formula T (m,n,q)：
Wherein:
FFT represents fourier transform;
the meaning is as follows:
the frequency bandwidth Δ f of the mth sub-band is continuously divided into Q sub-frequency intervals at equal intervals, which is expressed as: fr (Q), Q ═ 1, 2.., Q; two-dimensionally compensated acquisition signal g for each sub-band T (m, n) performing Fourier transform in each sub-frequency interval to obtain a Fourier transformed signal, namely a three-dimensional acquisition signal G T (m,n,q)。
Step 8, therefore, for the acquisition signal unit d T Performing fine search calculation on the M sub-frequency bands, the N frequency change rate estimation value sub-intervals and the Q sub-frequency intervals to obtain M, N and Q three-dimensional acquisition signals G T (m,n,q)；
Comparing M N Q three-dimensional collected signals G T The power of (m, n, q), the sub-band number corresponding to the power peak value, the sub-interval number of the frequency change rate estimation value and the sub-frequency interval number are the positions of the searched target signals, and then the target signals are searched;
in this step, a signal G is three-dimensionally acquired T The power of (m, n, q) is denoted P T (m, n, q) obtained by:
wherein:
P[N 0 ]represents N 0 The meaning of power of (a) is: power of the noise part after Fourier transform;
p [ S (m, n, q) ] represents the power of S (m, n, q), meaning: the target signal s (f (t)) is the power after partial fourier transform.
In this step, the peak power value is represented by P T (m T ,n T ,q T )；m T ,n T ,q T The meanings are respectively as follows: the subband number where the power peak is located, the frequency change rate estimation value subinterval number, and the subband interval number are expressed as:
(m T ,n T ,q T )＝max[P T (m,n,q)]| m∈[1,M],n∈[1,N],q∈[1,Q]
and 9, estimating a frequency spectrum detection estimation value according to the searched target signal.
In this step, the estimation value of the frequency spectrum detection of the target signal includes: accurate frequency estimation of a target signal
Accurate estimation of frequency rate of change
And accurate estimation value P/N of signal power noise spectral density ratio 0 ；
Obtaining a signal acquisition unit d by adopting the following formula T Accurate frequency estimation of intermediate target signal
And frequency rate of change accurate estimate
Obtaining a signal acquisition unit d by adopting the following formula T Accurate estimation value P/N of signal power noise spectral density ratio of intermediate target signal 0 ：
P/N 0 ＝P T (m T ,n T ,q T )/P[N 0 ]
P[N 0 ]＝N/(Δf/Q)
Obtaining a spectral curve G T (m T ,n T Q), wherein Q ═ 1, 2.., Q, includes amplitude frequency curves and phase frequency curves;
wherein:
spectral curve G T (m T ,n T Q) denotes the abscissa FR (q) and the ordinate G T (m T ,n T Q) the curve formed.
The computational complexity of the invention is analyzed as follows:
frequency bandwidth of B configured for ground station 0 The complexity of the computation can be expressed as O 1 ：
O 1 ∝(2B 0 ) 2 ×M×N
After dividing into Q sub-frequency intervals, the computational complexity can be represented as O 2 ：
O 2 ∝[(2B 0 /M) 2 ]×M×N＝(2B 0 ) 2 ×N/M
As can be seen, O 2 ＝O 1 /M 2 The more the number of sub-bands (i.e., M) employed, the lower the computational complexity of the present invention and the effective reduction in time required.
Simulation analysis test
Conditions are as follows:
1. the extraterrestrial celestial object target analog signal is based on the following three practical conditions, and the nominal downlink frequency f is assumed 0 ＝8431MHz：
Operating condition one, frequency f 0 +20kHz, frequency change rate-100 Hz/s, signal power noise spectral density ratio P/N 0 ＝20dBHz。
Operating mode two, frequency f 0 +20kHz, frequency change rate 100Hz/s, signal power noise spectral density ratio P/N 0 ＝20dBHz。
Operating mode three, frequency f 0 +20kHz, frequency change rate-10 Hz/s, signal power noise spectral density ratio P/N 0 ＝20dBHz。
2. And (3) data generation: the extraterrestrial celestial body target analog signal is converted into original data through acquisition, conversion and recording equipment (parameter configuration: acquisition and recording bandwidth is 2MHz, and quantization bit number is 16 bits), and the original data is loaded in a file form.
3. The signal level, Doppler and dynamic forecast file adds errors, frequency errors are +/-5 kHz, frequency change rate errors are +/-20 Hz/s, and signal power noise spectral density ratio errors are +/-5 dBHz on the basis of the step 1.
4. The number of subchannels (i.e., Q) is 20 and the number of frequency and frequency variation rate compensation grids is 50, 10.
By adopting the method of the invention, the frequency change rate and the signal power noise spectral density ratio are estimated, and the result is as follows:
operating condition one, frequency f 0 +20.001kHz, frequency change rate-99.017 Hz/s, P/N 0 ＝18.02dBHz。
Operating mode two, frequency f 0 +20.000kHz, frequency-doppler rate of change 98.68Hz/s, P/N 0 ＝18.10dBHz。
Operating mode three, frequency f 0 +19.999kHz, frequency-le rate of change-9.7 Hz/s, P/N 0 ＝17.73dBHz。
Therefore, the method can effectively monitor and estimate the signal spectrum of the extraterrestrial celestial body detector, and according to the results of the three working conditions, the frequency estimation precision reaches 1Hz, the frequency change rate estimation precision reaches 1Hz/s, and the signal power noise spectral density ratio estimation precision reaches 1.5 dBHz.
And (3) analyzing the computational complexity:
conventional analytical method, O 1 ∝(2×2×10 6 ) 2 ×20×10＝3.2×10 15 Process of the invention, O 2 ＝(2×2×10 6 /20) 2 ×20×10＝8×10 12 It can be seen that the computational complexity of the present invention is much less than that of the conventional method.
Analysis of practical application
Conditions are as follows:
1. mars lander signal theoretical value: frequency f 0 362.1kHz, frequency change rate 1Hz/s, signal power noise spectral density ratio P/N 0 ＝30dBHz。
2. And (3) data generation: the signals of the Mars lander are received by a ground antenna, are converted into original data through acquisition, conversion and recording equipment (parameter configuration: acquisition and recording bandwidth is 2MHz, and quantization bit number is 16bit), and are loaded in a file form.
3. By adopting the method of the invention, the frequency change rate and the signal power noise spectral density ratio are estimated, and the result is as follows:
estimation result of Mars lander signal: frequency f 0 362.103kHz, frequency change rate of 0.95Hz/s, P/N 0 31.29 dBHz. The spectrum of the estimated Mars lander signal is shown in FIG. 3.
Therefore, by applying the method, the spectrum of the Mars lander signal is successfully monitored and estimated, the deviation of the frequency of the signal from a theoretical value is 3Hz, the deviation of the frequency change rate from the theoretical value is 0.05Hz/s, and the deviation of the signal power noise spectral density ratio and a baseband test result is 1.29 dBHz. The accuracy of the estimated values of the present invention is thus demonstrated to be very high.
The frequency spectrum detection and estimation method for the signal of the extraterrestrial celestial body detector provided by the invention has the following advantages:
1. according to the invention, the terrestrial celestial body acquisition signals received by the ground antenna are divided into a plurality of sub-frequency bands, a plurality of sub-frequency intervals and a plurality of frequency change rate estimation value sub-intervals according to the frequency band range, so that rapid signal analysis processing is carried out, the calculation resources are reduced, and rapid signal analysis processing is realized.
2. The invention can realize high-precision estimation of signals, the estimation precision of the target signal frequency reaches 1Hz, the estimation precision of the frequency change rate reaches 1Hz/s, and the estimation precision of the signal noise spectral density ratio reaches 1.5 dBHz.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and improvements can be made without departing from the principle of the present invention, and such modifications and improvements should also be considered within the scope of the present invention.
Claims (7)
1. A method for spectrum sensing estimation of a signal from an extraterrestrial celestial object detector, comprising the steps of:
step 1, estimating a frequency range and a frequency change rate range of a detector signal to be acquired by a ground station in advance;
step 2, configuring signal acquisition parameters of the ground station according to the frequency range and the frequency change range of the detector signal to be acquired of the ground station estimated in advance in the step 1;
step 3, the ground station collects the detector signal by adopting the configured signal collection parameters to obtain an original collection signal;
step 4, dividing the original collected signals into a plurality of sections of collected signal units by taking the time length T as a period;
for each section of the acquired signal unit d T All execute step 5-step 8, search and collect the signal unit d T Target signal of (1):
step 5, collecting signal unit d T Is continuously divided into M sub-bands at equal intervals, the M sub-bandThe acquisition signal of (a) is expressed as: d is a radical of T (m),m＝1，2，3，…，M；
Step 6, continuously dividing the frequency change rate range estimated in the step 1 into N frequency change rate estimation value subintervals at equal intervals;
acquisition signal d for m-th sub-band T (m) respectively calculating the frequency and frequency change rate two-dimensional compensation value of each frequency change rate estimation value subinterval, thereby obtaining N two-dimensional compensated acquisition signals g T (m, N), wherein N is 1,2,3, N;
step 7, continuously dividing the frequency bandwidth of the mth sub-band into Q sub-frequency intervals at equal intervals;
collecting signal g after two-dimensional compensation of the mth frequency sub-band in the nth frequency change rate estimation value sub-interval T (m, n), respectively calculating the Fourier transform result of each sub-frequency interval, thereby obtaining Q signals after Fourier transform, namely the three-dimensional acquisition signal G T (m,n,q)；
Step 8, therefore, for the acquisition signal unit d T Performing fine search calculation on the M sub-frequency bands, the N frequency change rate estimation value sub-intervals and the Q sub-frequency intervals to obtain M, N and Q three-dimensional acquisition signals G T (m,n,q)；
Comparing M N Q three-dimensional collected signals G T The power of (m, n, q), the sub-band number corresponding to the power peak value, the sub-interval number of the frequency change rate estimation value and the sub-frequency interval number are the positions of the searched target signals, and then the target signals are searched;
and 9, estimating a frequency spectrum detection estimation value according to the searched target signal.
2. Method for the spectral detection estimation of a extraterrestrial celestial object detector signal as claimed in claim 1, wherein in step 1, the frequency range of the detector signal is denoted as [ f [ ] min ,f max ]The frequency change rate range is represented by [ f' min ,f′ max ]：
Wherein:
f min for detection of the need to collectA minimum estimate of the frequency of the signal;
f max the frequency maximum value estimation value of the detector signal to be acquired is obtained;
f′ min the minimum value of the frequency change rate of the detector signal to be acquired is estimated;
f′ max the frequency change rate maximum value estimation value of the detector signal to be acquired is obtained;
in step 2, the configured signal acquisition parameters of the ground station comprise the frequency bandwidth B of the ground station 0 And a center frequency f B Obtained by the following formula:
B 0 ＝f max -f min
3. method for spectral detection estimation of a extraterrestrial celestial object detector signal according to claim 2, wherein in step 5, a signal acquisition unit d is provided T Frequency bandwidth of (a) and frequency bandwidth B of the ground station arrangement 0 And are equal.
4. Method for spectral detection estimation of a extraterrestrial celestial object detector signal according to claim 3, wherein in step 5, the collected signal d for each sub-band T (m) is represented by:
wherein: n is a radical of an alkyl radical 0 Is noise, s (f (t)) is the target signal;
the meaning is as follows:
of the M subbands, only one subband has the target signal, assuming that there is the target signal in the z-th subband, where z is an unknown value; the other M-1 sub-bands are noise;
the step 6 specifically comprises the following steps:
the frequency change rate range [ f 'estimated in the step 1' min ,f′ max ]The equal-interval continuous division is carried out on N frequency change rate estimation value subintervals, and the frequency change rate estimation value subintervals are expressed as: e (N), N ═ 1,2,3., N; therefore, the width Δ f' ═ e (N)/N of each frequency rate estimation value subinterval;
acquisition signal d for m-th sub-band T (m) performing two-dimensional compensation of the frequency and the frequency change rate in the following manner to obtain a two-dimensionally compensated acquisition signal g T (m,n)：
Wherein:
j represents an imaginary unit;
f c (m) represents a lower boundary frequency compensation value of the mth sub-band;
f c ' (n) represents a lower boundary frequency change rate compensation value for the nth frequency change rate estimate subinterval;
wherein, f c (m) and f c ' (n) is obtained using the formula:
f c (m)＝(m-1)Δf+f min
f c ′(n)＝(n-1)Δf′+f′ min
wherein: Δ f ═ B 0 and/M, representing the frequency bandwidth of each sub-band.
5. The method for detecting and estimating the frequency spectrum of the extraterrestrial celestial object detector signal according to claim 4, wherein the step 7 is specifically:
collecting signal g after two-dimensional compensation of the mth frequency sub-band in the nth frequency change rate estimation value sub-interval T (m, n) is converted into a three-dimensional acquisition signal G by adopting the following formula T (m,n,q)：
Wherein:
FFT represents fourier transform;
the meaning is as follows:
the frequency bandwidth Δ f of the mth sub-band is continuously divided into Q sub-frequency intervals at equal intervals, which is expressed as: fr (Q), Q ═ 1, 2.., Q; two-dimensionally compensated acquisition signal g for each sub-band T (m, n) performing Fourier transform in each sub-frequency interval to obtain a Fourier transformed signal, namely a three-dimensional acquisition signal G T (m,n,q)。
6. Method for spectral detection estimation of a extraterrestrial celestial object detector signal as claimed in claim 5, wherein in step 8, the signal G is acquired in three dimensions T The power of (m, n, q) is represented as P T (m, n, q) obtained by:
wherein:
P[N 0 ]represents N 0 The meaning of power of (a) is: power after fourier transform of the noise part;
p [ S (m, n, q) ] represents the power of S (m, n, q), meaning: the target signal s (f (t)) is partially fourier-transformed.
7. The method of claim 6, wherein in step 8, the power peak is detectedIs expressed as P T (m T ,n T ,q T )；m T ,n T ,q T The meanings are respectively as follows: the subband number where the power peak is located, the frequency change rate estimation value subinterval number, and the subband interval number are expressed as:
(m T ,n T ,q T )＝max[P T (m,n,q)]| m∈[1,M],n∈[1,N],q∈[1,Q]
the step 9 specifically comprises the following steps:
the frequency spectrum detection estimation value of the target signal comprises the following steps: accurate frequency estimation of target signal
Accurate estimation value of frequency change rate
And accurate estimation value P/N of signal power noise spectral density ratio 0 ；
Obtaining a signal acquisition unit d by adopting the following formula T Accurate frequency estimation of intermediate target signal
And frequency rate of change accurate estimate
Obtaining a signal acquisition unit d by adopting the following formula T Accurate estimation value P/N of signal power noise spectral density ratio of intermediate target signal 0 ：
P/N 0 ＝P T (m T ,n T ,q T )/P[N 0 ]
P[N 0 ]＝N/(Δf/Q)
Obtaining a spectral curve G T (m T ,n T Q), wherein Q ═ 1, 2.., Q, includes amplitude frequency curves and phase frequency curves;
wherein:
spectral curve G T (m T ,n T Q) denotes the abscissa FR (q) and the ordinate G T (m T ,n T Q) the resulting curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210787392.5A CN114924121B (en) | 2022-07-04 | 2022-07-04 | Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210787392.5A CN114924121B (en) | 2022-07-04 | 2022-07-04 | Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114924121A true CN114924121A (en) | 2022-08-19 |
| CN114924121B CN114924121B (en) | 2023-03-10 |
Family
ID=82816142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210787392.5A Active CN114924121B (en) | 2022-07-04 | 2022-07-04 | Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114924121B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080177490A1 (en) * | 2005-01-14 | 2008-07-24 | Rohde & Schwarz Gmbh & Co. Kg | Method and system for the detection and/or removal of sinusoidal interference signals in a noise signal |
| JP2012078187A (en) * | 2010-09-30 | 2012-04-19 | Toshiba Corp | Parameter detector, radar device, guidance system and parameter detection method |
| CN104852875A (en) * | 2015-03-26 | 2015-08-19 | 中国人民解放军理工大学通信工程学院卫星通信军队重点实验室 | Frequency offset estimation method for high-dynamic large-frequency-offset burst signals |
| CN109495410A (en) * | 2018-09-28 | 2019-03-19 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | High dynamic PCM/FM signal(-) carrier frequency precise Estimation Method |
| CN111624402A (en) * | 2020-05-31 | 2020-09-04 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Method for accurately estimating carrier frequency of weak PM signal |
-
2022
- 2022-07-04 CN CN202210787392.5A patent/CN114924121B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080177490A1 (en) * | 2005-01-14 | 2008-07-24 | Rohde & Schwarz Gmbh & Co. Kg | Method and system for the detection and/or removal of sinusoidal interference signals in a noise signal |
| JP2012078187A (en) * | 2010-09-30 | 2012-04-19 | Toshiba Corp | Parameter detector, radar device, guidance system and parameter detection method |
| CN104852875A (en) * | 2015-03-26 | 2015-08-19 | 中国人民解放军理工大学通信工程学院卫星通信军队重点实验室 | Frequency offset estimation method for high-dynamic large-frequency-offset burst signals |
| CN109495410A (en) * | 2018-09-28 | 2019-03-19 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | High dynamic PCM/FM signal(-) carrier frequency precise Estimation Method |
| CN111624402A (en) * | 2020-05-31 | 2020-09-04 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Method for accurately estimating carrier frequency of weak PM signal |
Non-Patent Citations (1)
| Title |
|---|
| 薛海中等: "基于局部频谱连续细化的高精度频率估计算法", 《西安电子科技大学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114924121B (en) | 2023-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102098074B (en) | High-dynamic weak-signal rapid capture method for direct sequence spread spectrum system | |
| Goetz et al. | An all-sky search algorithm for continuous gravitational waves from spinning neutron stars in binary systems | |
| CN107436427B (en) | Spatial target motion track and radiation signal correlation method | |
| CN109669183B (en) | Geosynchronous orbit SAR moving target imaging processing device based on Keystone and time-frequency transformation | |
| CN107450055B (en) | High-speed maneuvering target detection method based on discrete linear frequency modulation Fourier transform | |
| CN106772227A (en) | A kind of unmanned plane direction determining method based on the identification of vocal print multiple-harmonic | |
| CN112199996B (en) | Rolling bearing diagnosis method based on parameter self-adaptive VMD and fast Hoyer spectrogram indexes | |
| CN109459745B (en) | Method for estimating speed of moving sound source by using radiation noise | |
| CN103513244A (en) | Multi frame coherent accumulation target track-before-detect method based on dynamic programming | |
| CN103364783A (en) | Moving target radial velocity non-fuzzy estimation method based on single-channel SAR (synthetic aperture radar) | |
| Altadill et al. | A method for real-time identification and tracking of traveling ionospheric disturbances using ionosonde data: First results | |
| Tyler et al. | Radio wave scattering observations of the solar corona First-order measurements of expansion velocity and turbulence spectrum using Viking and Mariner 10 spacecraft | |
| CN104502898A (en) | Maneuvering target parameter estimation method by combining correction RFT (Radon-Fourier Transform) and MDCFT (Modified Discrete Chirp-Fourier Transform) | |
| JP5483343B2 (en) | Ocean surface water quality measurement method, water quality measurement device, and water quality measurement program | |
| CN114924121B (en) | Frequency spectrum detection estimation method for extraterrestrial celestial body detector signal | |
| CN113296073A (en) | Amplitude-phase joint blade clearance measurement method based on envelope constraint | |
| Huang et al. | A practical fundamental frequency extraction algorithm for motion parameters estimation of moving targets | |
| CN102520404A (en) | SAR (Synthetic Aperture Radar) Doppler fuzzy number estimation method based on optimal image quality | |
| CN114330163A (en) | Modeling method for high-frequency ground wave over-the-horizon radar typhoon-ionosphere disturbance dynamics model | |
| CN101846739A (en) | Mixed domain emulation method of SAR (Synthetic Aperture Radar) extended scene primary data | |
| Ding et al. | Slope synchronous chirplet transform and its application to tacho-less order tracking of rotating machineries | |
| CN113009523A (en) | Doppler frequency estimation and compensation method and system for long-time coherent integration capture | |
| Jiao et al. | Performance evaluations of an equatorial GPS amplitude scintillation detector using a machine learning algorithm | |
| CN112113767A (en) | Bearing fault diagnosis method based on proportional band selection criterion | |
| CN114942419B (en) | Method for extracting three-degree-of-freedom micro-motion characteristics of scattering points of ship under long accumulation time |
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 |