CN104007453B - The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary - Google Patents
The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary Download PDFInfo
- Publication number
- CN104007453B CN104007453B CN201410190833.9A CN201410190833A CN104007453B CN 104007453 B CN104007453 B CN 104007453B CN 201410190833 A CN201410190833 A CN 201410190833A CN 104007453 B CN104007453 B CN 104007453B
- Authority
- CN
- China
- Prior art keywords
- signal
- interference
- road
- carried out
- frequency domain
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Noise Elimination (AREA)
Abstract
The invention discloses the frequency domain spatial domain associating anti-interference method of a kind of probabilistic search auxiliary, implement according to following steps: step 1, Signal Pretreatment, be the most each divided into I road and Q road signal to export；Step 2, frequency domain process, and the I road finally exported step 1 and Q road signal are divided into two branch roads to be respectively processed, and the signal after then postponing is synthesized by adder with the signal obtained；Step 3, spatial processing, the output signal being finally synthesizing step 2 carries out spatial domain power inversion adaptive-filtering, suppresses the broadband interference in signal, obtains final output signal.The inventive method solve the associating of frequency domain spatial domain anti-interference time, arrowband interference is submerged in broadband interference and makes Anti-interference algorithm reduce performance or malfunctioning problem, expands the scope of application of frequency domain spatial domain associating Anti-interference algorithm.
Description
Technical field
The invention belongs to GLONASS Anti-Jamming Technique field, relate to a kind of probability and search
The frequency domain spatial domain associating anti-interference method of rope auxiliary.
Background technology
Along with GLONASS, (Global Navigation Satellite System is called for short
GNSS) technology is day by day ripe, no matter at military field or the equal table of civil area GNSS
Revealing its unrivaled importance, the research and development of GNSS technology and application are by increasingly
The concern of many people.But the signal power of aeronautical satellite is low, it is easy to by various consciously or unconsciously
Interference effect, particularly can be suppressed by strong jamming in wartime, and cisco unity malfunction.
In prior art, when complex electromagnetic environment, arrowband interference and broadband interference the feelings deposited
Under condition, airspace filter, space-time joint (STAP), frequency domain is generally used to cascade spatial domain scheduling algorithm
Carry out anti-interference process.Although simple spatial domain filter algorithms can filter simultaneously arrowband interference and
Broadband interference, but arrowband interference will reduce array antenna degree of freedom, reduce anti-interference number.Empty
Time united algorithm complicated, amount of calculation is excessive, high to hardware requirement.Frequency domain cascade spatial domain is first used
Frequency domain filtering technology filters arrowband interference, then filters more interference through airspace filter, but works as
Arrowband coexists with broadband interference, and when arrowband jamming power is less than broadband interference, carries out frequency domain filter
Ripple inspection does not measures arrowband interference position, and common frequency domain cascade spatial domain Anti-interference algorithm can be caused to lose
Effect.Therefore, existing technology is respectively arranged with pluses and minuses, in addition it is also necessary to develop new aeronautical satellite anti-interference
Method.
Summary of the invention
The frequency domain spatial domain associating that it is an object of the invention to provide a kind of probabilistic search auxiliary is anti-interference
Method, solves prior art in interference RST, and arrowband interference and broadband interference are also deposited,
The especially mid frequency of arrowband interference falls at broadband interference band alleviating distention in middle-JIAO, and power is done less than broadband
Disturb, cause systematic function to be decreased obviously or the problem of cisco unity malfunction.
The technical solution adopted in the present invention is, the frequency domain spatial domain associating of a kind of probabilistic search auxiliary
Anti-interference method, implements according to following steps:
Step 1, Signal Pretreatment；
Step 2, frequency domain process；
Step 3, spatial processing.
The invention has the beneficial effects as follows: disturbed GNSS signal is processed through radio-frequency front-end
The intermediate-freuqncy signal obtained carries out A/D collection, Digital Down Convert, filtering, interpolation, then adopts
With the frequency domain filtering suppression arrowband interference of overlapping windowing FFT, finally signal is passed through again spatial domain
Adaptive-filtering module carries out broadband interference and filters, and completes the purpose of filtering interference signals.Above-mentioned
The frequency domain compact technique suppression arrowband interference of carrying out overlapping windowing FFT time, owing to arrowband is done
Disturb and may be submerged in a period of time in broadband interference, after FFT, can not find arrowband
Interference peak value, applied probability statistic algorithm, arrowband interference is easily detected out and carries out effectively
Filter, also can use the preferable statistic algorithm of real-time as requested.For airspace filter module
Use quick self-adapted filtering algorithm to carry out power inversion filtering, be greatly improved the dynamic of system
Energy.This algorithm improves the associating of frequency domain the spatial domain stability of Anti-interference algorithm, real-time, expands
The scope of application of frequency domain spatial domain associating Anti-interference algorithm improves the capacity of resisting disturbance of system.
Accompanying drawing explanation
Fig. 1 is the system block diagram that the inventive method is relied on；
Fig. 2 is the theory diagram of the Signal Pretreatment part in the inventive method；
Fig. 3 is the theory diagram of the frequency domain filtering part in the inventive method；
Fig. 4 is the theory diagram of the airspace filter part in the inventive method.
Detailed description of the invention
The present invention is described in detail with detailed description of the invention below in conjunction with the accompanying drawings.
With reference to Fig. 1, it is the frequency domain spatial domain associating anti-interference method of probabilistic search of the present invention auxiliary,
Three parts of the system relied on, i.e. Signal Pretreatment part, frequency domain filtering part and spatial domain filter
Ripple part, wherein,
The structure of Signal Pretreatment part is, including multi-channel rf front end, A/D data acquisition unit,
Each road A/D data acquisition unit is with digital down converter, low pass filter and withdrawal device successively
Connect；
The structure of frequency domain filtering part is, including two branch roads, wherein a branch road is, including
Window added device, FFT processing module, threshold value generate and spectral line processing module, IFFT processing module,
50% delay disposal module is connected in sequence；Another branch road is, including 50% delay disposal
Module, window added device, FFT processing module, threshold value generate and at spectral line processing module, IFFT
Reason module is connected in sequence, and I, Q two-way of two branch roads is last to be connected with adder simultaneously,
Output signal after corresponding addition is for next step airspace filter.
Window added device (Blackman window) have employed Sidelobe Suppression effect preferable Blackman window
Function is implemented, and is substantially reduced the leakage of signal spectrum；Probability statistics are then according to maximum in every frame
Occur that the statistical value of position compares with setting thresholding, be then mono-tone interference more than thresholding, enter
Row processes；
The structure of airspace filter part is, by filter construction and adaptive filter algorithm two parts
Composition, adaptive filter algorithm is the gradient variable-step least mean square algorithm of present invention innovation
(being called for short GVSS-NLMS).
The frequency domain spatial domain associating anti-interference method of probabilistic search of the present invention auxiliary, depends on above-mentioned
Structure and principle, implement according to following steps:
Step 1, Signal Pretreatment
With reference to Fig. 2, the intermediate-freuqncy signal obtained after radio frequency front-end processing is carried out A/D data acquisition,
The data signal collected is carried out the process of Digital Down Convert；After again Digital Down Convert being processed
The I road obtained and Q road signal carry out low-pass filtering simultaneously, carry out extraction process the most simultaneously,
Finally I road and Q road signal each it is divided into export, for subsequent treatment；
Step 2, frequency domain process
With reference to Fig. 3, what frequency domain filtering processed comprises the concrete steps that: the I finally exported step 1
Road and Q road signal are divided into two branch roads to be respectively processed,
One branch road is: 2.11) I road and Q road signal carry out windowing process, drawing of window function
Enter the sequence boundaries of blocking obtained so that follow-up fast Fourier transform (FFT) to become
Smooth, therefore, it is possible to alleviate the energy leakage of follow-up FFT；
2.12) signal of step 1 output is carried out FFT, output signal spectrum；
2.13) signal spectrum data are carried out probability statistics process, detect the arrowband in frequency spectrum
Interference, compares the ratio of arrowband interference and broadband interference maximum, the position occurred according to maximum
Putting the most identical, identical, statistics number n adds 1, and difference then statistics number n subtracts 1；
Set upper limit NmaxWith a lower limit NminIf, statistics number beyond the upper limit or
Keep constant less than lower limit then n；Set 2 thresholding α and β, n more than α (β < α < Nmax)
One mono-tone interference is described, has carried out zero setting process；If n is less than β (β > Nmin) then recognize
For not having mono-tone interference, test of many times is needed to obtain generating suitable threshold value；
2.14) according to the threshold value generated, frequency spectrum data is carried out spectral line process；
2.15) data signal after processing spectral line carries out IFFT conversion (inverse FFT),
Signal after being processed by frequency domain carries out IFFT and is transformed to time-domain signal；
2.16) again time-domain signal is carried out 50% delay disposal；
Another branch road is:
2.21) signal of step 1 output is carried out 50% delay；
2.22) I road and Q road signal being carried out windowing process, the introducing of this window function makes soon
It is smoothened that what speed Fourier transformation (FFT) obtained blocks sequence boundaries, therefore, it is possible to
Alleviate the energy leakage of follow-up FFT,
2.23) FFT is carried out, signal spectrum after output windowing；
2.24) signal spectrum data are carried out probability statistics process, detect the arrowband in frequency spectrum
Interference, compares the ratio of arrowband interference and broadband interference maximum, the position occurred according to maximum
Putting the most identical, identical, statistics number n adds 1, and difference then statistics number n subtracts 1；
Set upper limit NmaxWith a lower limit NminIf, statistics number beyond the upper limit or
Keep constant less than lower limit then n；Set 2 thresholding α and β, n more than α (β < α < Nmax)
One mono-tone interference is described, has carried out zero setting process；If n is less than β (β > Nmin) then recognize
For not having mono-tone interference, test of many times is needed to obtain generating suitable threshold value；
2.25) according to the threshold value generated, frequency spectrum data is carried out spectral line process；
2.26) data after processing spectral line carry out IFFT conversion (inverse FFT), will frequency
Signal after the process of territory carries out IFFT and is transformed to time-domain signal；
2.27) by step 2.16) postpone after signal and step 2.26) signal that obtains passes through
Adder synthesizes；
Step 3, airspace filter process
With reference to Fig. 4, it is adaptive that the output signal being finally synthesizing step 2 carries out spatial domain power inversion
Should filter, the broadband interference in signal is suppressed, obtain final output signal,
The renewal of power inversion adaptive filter algorithm weights based on GVSS-LMS algorithm is public
Formula is:
Wherein X (n) is input signal vector, XTN () is the transposition of input signal vector, W (n) is power
Vector, μ (n) is variable step size, s0For constraint vector, W (0) s0=1 for ensureing airspace filter the
One tunnel output is not zero, and β is smoothing parameter constant (levels off to 1), and γ is the constant more than zero,
G (n) is smooth gradient vector, and I is unit matrix.
The frequency domain spatial domain associating anti-interference method of probabilistic search of the present invention auxiliary, its operation principle
It is:
Time-domain signal is mapped to frequency domain and goes to carry out the processing method of AF panel, signal can be made
Frequency characteristic can show more prominent, thus obtains more preferable treatment effect.Due to arrowband
Interference and broadband interference and deposit and arrowband jamming power less than broadband interference time, frequency domain filtering FFT
After can't detect or detect inaccurate arrowband interference frequency.Use the side of probability statistics auxiliary
Arrowband interference is detected the maximum by average energy value and arrowband disturb and compares by method, energy
Enough effectively detect arrowband interference position and filtered, the most again after IFFT will process
Signal transform to time domain and carry out airspace filter.Airspace filter is with based on linear restriction minimum equal
The Power-inversion algorithm of side (LCMV) carries out algorithm iteration, is output as purpose with power minimum.
GVSS-NLMS algorithm used by power inversion utilizes a firstorder filter to gradient vector
Carry out smooth to reduce influence of noise, use with smooth after gradient vector to update step-length, energy
The most effectively estimate the optimal step size that the mean square norm of system errors minimizes, improve and receive
Hold back speed, it is achieved that the Fast Convergent on the premise of not reducing steady-state error, compare existing
New variable step-size LMS has faster convergence rate, greatly strengthen real-time.
Claims (1)
1. the frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary, it is characterised in that according to
Lower step is implemented:
Step 1, Signal Pretreatment
The intermediate-freuqncy signal obtained after radio frequency front-end processing is carried out A/D data acquisition, the data that will collect
Signal carries out the process of Digital Down Convert；The I road obtained after again Digital Down Convert being processed and Q road signal
Carry out low-pass filtering simultaneously, carry out extraction process the most simultaneously, be the most each divided into I road and Q road signal
Export；
Step 2, frequency domain process
The I road finally exported step 1 and Q road signal are divided into two branch roads to be respectively processed,
One of them branch road is:
2.11) I road and Q road signal carry out windowing process；
2.12) signal of step 1 output is carried out FFT, output signal spectrum；
2.13) signal spectrum data are carried out probability statistics process, detect the arrowband interference in frequency spectrum,
Relatively arrowband interference and the ratio of broadband interference maximum, the most identical according to the position that maximum occurs,
Identical, statistics number n adds 1, and difference then statistics number n subtracts 1；Set upper limit NmaxWith one
Lower limit NminIf statistics number is beyond the upper limit or less than lower limit, and n keeps constant；Set 2 thresholdings
α and β, n are more than α, β < α < Nmax, a mono-tone interference has been described, has carried out zero setting process；If
N is less than β, β > Nmin, then it is assumed that there is no mono-tone interference, obtain generating suitable threshold value；
2.14) according to the threshold value generated, frequency spectrum data is carried out spectral line process；
2.15) data signal after processing spectral line carries out IFFT conversion, the signal after being processed by frequency domain
Carry out IFFT and be transformed to time-domain signal；
2.16) again time-domain signal is carried out 50% delay disposal；
Another branch road is:
2.21) signal of step 1 output is carried out 50% delay；
2.22) I road and Q road signal are carried out windowing process,
2.23) FFT is carried out, signal spectrum after output windowing；
2.24) signal spectrum data are carried out probability statistics process, detect the arrowband interference in frequency spectrum,
Relatively arrowband interference and the ratio of broadband interference maximum, the most identical according to the position that maximum occurs,
Identical, statistics number n adds 1, and difference then statistics number n subtracts 1；Set upper limit NmaxWith one
Lower limit NminIf statistics number is beyond the upper limit or less than lower limit, and n keeps constant；Set 2 thresholdings
α and β, n are more than α, β < α < Nmax, a mono-tone interference has been described, has carried out zero setting process；If
N is less than β, β > Nmin, then it is assumed that there is no mono-tone interference, generate suitable threshold value；
2.25) according to the threshold value generated, frequency spectrum data is carried out spectral line process；
2.26) data after processing spectral line carry out IFFT conversion, and the signal after being processed by frequency domain is carried out
IFFT is transformed to time-domain signal；
2.27) by step 2.16) postpone after signal and step 2.26) signal that obtains passes through adder
Synthesize；
Step 3, spatial processing
The output signal being finally synthesizing step 2 carries out spatial domain power inversion adaptive-filtering, in signal
Broadband interference suppress, obtain final output signal,
The more new formula of power inversion adaptive filter algorithm weights based on GVSS-LMS algorithm is:
Wherein, X (n) is input signal vector, XTN () is the transposition of input signal vector, W (n) is
Weight vector, μ (n) is variable step size, s0For constraint vector, W (0) s0=1 for ensureing airspace filter first
Road output is not zero, and β is smoothing parameter constant, and γ is the constant more than zero, g (n) be smooth gradient to
Amount, I is unit matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410190833.9A CN104007453B (en) | 2014-05-07 | 2014-05-07 | The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410190833.9A CN104007453B (en) | 2014-05-07 | 2014-05-07 | The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104007453A CN104007453A (en) | 2014-08-27 |
| CN104007453B true CN104007453B (en) | 2016-12-07 |
Family
ID=51368197
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410190833.9A Expired - Fee Related CN104007453B (en) | 2014-05-07 | 2014-05-07 | The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104007453B (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105375938A (en) * | 2014-10-16 | 2016-03-02 | 航天恒星科技有限公司 | Signal processing method based on inertial navigation auxiliary space frequency, apparatus and receiver |
| CN104393883B (en) * | 2014-11-01 | 2017-03-15 | 中国民航大学 | ADS B pressing type disturbance restraining methods based on improved power inversion method |
| CN105572695A (en) * | 2014-11-07 | 2016-05-11 | 航天恒星科技有限公司 | Anti-interference method and system |
| CN104360354B (en) * | 2014-12-05 | 2018-06-08 | 北京北斗星通导航技术股份有限公司 | Null tone anti-interference processing method and device |
| CN107689814A (en) * | 2016-08-03 | 2018-02-13 | 北京遥感设备研究所 | A kind of MSK band spread receivers narrow-band interference rejection method |
| CN107121685A (en) * | 2017-06-08 | 2017-09-01 | 南京理工大学 | A kind of miniature spaceborne high-dynamic GNSS receiver and its air navigation aid |
| CN107589430A (en) * | 2017-09-07 | 2018-01-16 | 中国民航大学 | ADS B pressing type disturbance restraining methods based on minimum dispersion method |
| CN109541646A (en) * | 2018-12-27 | 2019-03-29 | 中国电子科技集团公司第二十研究所 | A kind of large gain satellite navigation aerial anti-interference processing method |
| CN110504988B (en) * | 2019-04-22 | 2022-04-19 | 长沙翼盾电子科技有限公司 | Time-frequency domain combined narrow-band interference resisting method |
| CN111812682A (en) * | 2020-07-24 | 2020-10-23 | 华力智芯(成都)集成电路有限公司 | Narrow-band interference resistant circuit |
| CN113721273A (en) * | 2021-07-27 | 2021-11-30 | 四创电子股份有限公司 | Double-antenna Beidou portable monitoring terminal device |
| CN115656940B (en) * | 2022-11-23 | 2023-05-30 | 中国人民解放军空军预警学院 | Space-frequency domain narrowband composite interference resisting method and system |
| CN117560259B (en) * | 2023-11-14 | 2024-08-13 | 成都玖锦科技有限公司 | MSK and 2FSK signal identification method and system based on bandwidth estimation |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101561490B (en) * | 2009-06-01 | 2012-06-27 | 中国民航大学 | Blind self-adaptive GPS interference suppression method based on codeword structure |
| CN101718873A (en) * | 2009-11-13 | 2010-06-02 | 西安电子科技大学 | Homing signal space-time anti-interference digital signal processor |
| CN101799551B (en) * | 2010-02-08 | 2014-03-26 | 中国民航大学 | Inhibition method of space-time blind self-adaptive GPS (Global Positioning System) interference based on despreading and respreading technology |
| CN102207549A (en) * | 2010-05-05 | 2011-10-05 | 北京泰豪联星技术有限公司 | Integrated anti-interference satellite navigation receiving system and anti-interference processing method thereof |
| CN101986634B (en) * | 2010-09-26 | 2013-02-06 | 北京大学 | Time-frequency expanding anti-jamming method, equipment and system based on lapped transformation algorithm |
| CN102645659B (en) * | 2012-04-18 | 2014-01-08 | 西安理工大学 | Frequency domain filtering algorithm based on frequency statistics |
-
2014
- 2014-05-07 CN CN201410190833.9A patent/CN104007453B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104007453A (en) | 2014-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104007453B (en) | The frequency domain spatial domain associating anti-interference method of probabilistic search auxiliary | |
| CN103197300B (en) | Real-time processing method for cancellation of direct wave and clutter of external radiation source radar based on graphic processing unit (GPU) | |
| CN102645659B (en) | Frequency domain filtering algorithm based on frequency statistics | |
| CN103728594B (en) | Based on the external illuminators-based radar sea clutter disturbance restraining method of hyperchannel NLMS | |
| CN101937088B (en) | Space and time two-dimensional anti-interference realization method for satellite navigation receiving equipment | |
| CN104808219A (en) | Novel space-time joint anti-interference method | |
| CN103064090A (en) | Anti-interference method | |
| CN103675842B (en) | A kind of empty frequency combines anti-interference realization method | |
| CN104360355B (en) | Anti-interference method and device | |
| CN103700374A (en) | Method for determining system time delay in acoustic echo cancellation and acoustic echo cancellation method | |
| CN103558614B (en) | Passage in a kind of GPS and observation territory joined multi-path suppressing method | |
| CN103353550A (en) | Method for measuring signal frequency and harmonic parameters of electric power system | |
| CN110504988B (en) | Time-frequency domain combined narrow-band interference resisting method | |
| CN107144879A (en) | A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation | |
| CN103312373B (en) | Big Dipper generation time service type receiver suppresses the time-domain device of arrowband interference | |
| CN104898132A (en) | Navigation anti-interference algorithm combining threshold processing and space-frequency adaptive algorithm | |
| Li et al. | High speed maneuvering target detection based on joint keystone transform and CP function | |
| CN106289239A (en) | A kind of method eliminating the interference of wideband time domain in the pulsar data time of advent | |
| CN104502925A (en) | Narrowband interference resisting system and method based on adaptive signal processing | |
| CN104076369A (en) | Frequency domain anti-interference method and device based on adaptive threshold judgment | |
| CN103763229A (en) | CDMA communication anti-interference system based on frequency domain filtering | |
| CN108919227B (en) | Multichannel FBLMS implementation method based on GPU acceleration | |
| CN103344945A (en) | Direct wave and method for inhibiting multi-path interference of direct wave | |
| CN102508264B (en) | Method for inhibiting pulse interference of satellite navigation system | |
| CN103871422A (en) | Time-domain aliasing blind signal separation method based on three-factor iteration joint block diagonalization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161207 |