CN112578410A - Evolution band-limited Gaussian noise interference algorithm for GPS based on LMS - Google Patents

Evolution band-limited Gaussian noise interference algorithm for GPS based on LMS Download PDF

Info

Publication number
CN112578410A
CN112578410A CN202011468606.XA CN202011468606A CN112578410A CN 112578410 A CN112578410 A CN 112578410A CN 202011468606 A CN202011468606 A CN 202011468606A CN 112578410 A CN112578410 A CN 112578410A
Authority
CN
China
Prior art keywords
algorithm
interference
lms
gaussian noise
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011468606.XA
Other languages
Chinese (zh)
Other versions
CN112578410B (en
Inventor
吴涛
邹秋白
吴吉祥
王宏
侯者非
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Air Force Communication Sergeant School Of Pla
Original Assignee
Air Force Communication Sergeant School Of Pla
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Air Force Communication Sergeant School Of Pla filed Critical Air Force Communication Sergeant School Of Pla
Priority to CN202011468606.XA priority Critical patent/CN112578410B/en
Publication of CN112578410A publication Critical patent/CN112578410A/en
Application granted granted Critical
Publication of CN112578410B publication Critical patent/CN112578410B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/21Interference 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)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention relates to an evolutionary band-limited Gaussian noise interference algorithm for a GPS (global positioning system) based on LMS (least mean square), and an expression function x (nT) of an interference signal is determined; establishing an interference machine and a GPS receiver model; optimizing a flicker function r (nT) in the expression function x (nT) by adopting an evolutionary algorithm; the degree of influence of the interfering signal on the receiver is determined. The evolution band-limited Gaussian noise interference algorithm has more obvious interference effect than band-limited Gaussian noise interference when different interference amounts exist.

Description

Evolution band-limited Gaussian noise interference algorithm for GPS based on LMS
Technical Field
The invention relates to the technical field of GPS signal interference, in particular to an evolutionary band-limited Gaussian noise interference algorithm for a GPS (Least mean square) based on LMS (Least mean square).
Background
A GPS receiver is an instrument that receives global positioning system satellite signals and determines the position of the ground space. The interference signal is a signal that impairs reception of a desired signal. Signal reception by an enemy GPS receiver may be opposed by signal interference.
Suppressive interference, electronic interference that obscures or completely masks the useful signal received by the hostile electronic device from proper operation. The traditional pressing type GPS jammer can not correctly decode a GPS receiver by transmitting an interference signal with higher power than a GPS signal. The zero setting GPS receiver uses an antenna array, and the self-adaptive algorithm is used for nulling the interference direction, so that the obvious anti-interference effect is achieved on the traditional suppression type interference mode.
Currently, most researchers are still working from the perspective of interfering with ordinary GPS receivers. In 58-62 pages of 58 th-62 th of book 15 of university of air force engineering (nature science edition), of mao hu et al in 2014, "gaussian noise interference parameter selection analysis with band limit for GPS military code signals", and in 257 th-267 pages of 4 th of book 24 of computer measurement and control, "wang jiao in 2016", a gaussian white noise interference optimum parameter selection strategy and a method for judging an optimum interference signal are respectively provided.
Prior to 2017, there were fewer interference studies with zeroed GPS. In a text of interference and anti-interference technology research on a satellite navigation receiver adopting an array antenna, Liyawei in 2017 provides a pulse interference method aiming at a zero GPS, relative time delay among a plurality of interferences is periodically adjusted, and a certain performance gain is obtained. However, no research result for improving the interference effect by the optimization algorithm is disclosed.
Disclosure of Invention
The invention aims to provide an evolutionary band-limited Gaussian noise interference algorithm for a GPS based on LMS, which can obviously improve the effect of signal interference.
The technical scheme for solving the technical problems is as follows:
an evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS, comprising the steps of:
step 1: determining an expression function x (nT) of the interference signal;
step 2: establishing an interference machine and a GPS receiver model;
and step 3: optimizing a flicker function r (nT) in the expression function x (nT) by adopting an evolutionary algorithm;
and 4, step 4: the degree of influence of the interfering signal on the receiver is determined.
An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS, as described above, in step 1,
the complex field expression of the single interference signal x (nT) is shown as the following formula (1);
x(nT)=r(nT)[f1(nT)cos(2πfcnT+β)+jf2(nT)sin(2πfcnT+β)] (1)
in the formula (1), fc1575.42MHz, representing the GPS L1 signal carrier frequency; β represents an initial phase; t1/fsRepresents a sampling period; f. ofsRepresents the sampling frequency; f. of1(nT)) and f2(nT) represents white gaussian noise after h (z) filtering; r (nt) represents a flicker function; n denotes a sample number, and j denotes an imaginary unit.
In the above algorithm for the LMS-based GPS based evolutionary band-limited gaussian noise interference, in step 2, the process of establishing the jammer and GPS receiver models includes the following steps:
step 201: determining the arrangement mode of the jammer and the GPS receiver;
step 202: determining an expression function r (nT) of an array element receiving signal of a GPS receiver;
step 203: obtaining the optimal value W of the weight vector W by adopting an LMS algorithmopt
Step 204: calculating directional diagram
Figure BDA0002829828740000021
In the above-mentioned LMS-based GPS evolutionary band-limited gaussian noise interference algorithm, in step 203, the LMS algorithm is adopted to obtain the optimal value W of the weight vector WoptComprises the following steps:
step 2031: calculating weighted average signals y (K) of the signals of the other array elements except the reference array element;
step 2032: calculating an error signal e (k);
step 2033: the LMS algorithm obtains a weight vector W by utilizing a steepest descent method;
step 2034: obtaining optimal weight vector W by using LMS algorithmopt
In the above evolutionary band-limited gaussian noise interference algorithm for the LMS-based GPS, in step 3, the flicker function r (nt) in the expression function x (nt) is optimized by the evolutionary algorithm, and is implemented by optimizing the sequence a (k); the procedure for optimizing the a (k) sequence is as follows:
step 301: firstly, initializing parameters of an optimization algorithm;
step 302: for each population member, calculating a variation function updating gene of the gene mutation;
step 303: calculating a variation function updating gene of the cross variation;
step 304: calculating each member evaluation value scores;
step 305: the maximum value of all the evaluation values scores is calculated.
The scheme adopted by the invention has the beneficial effects that:
the invention relates to an evolutionary band-limited Gaussian noise interference algorithm for a GPS (global positioning system) based on LMS (least mean square), and an expression function x (nT) of an interference signal is determined; establishing an interference machine and a GPS receiver model; optimizing a flicker function r (nT) in the expression function x (nT) by adopting an evolutionary algorithm; the degree of influence of the interfering signal on the receiver is determined.
The invention optimizes the band-limited Gaussian noise interference flicker function through the evolutionary algorithm, thereby disturbing the convergence of the self-adaptive algorithm based on LMS in the GPS receiver, leading the average zero depth of the GPS receiver under the equal power to be lower and realizing better interference effect.
The evolution band-limited Gaussian noise interference algorithm has more obvious interference effect than band-limited Gaussian noise interference when different interference amounts exist.
Drawings
Fig. 1 is a schematic diagram of the distribution of jammers interfering with a GPS receiver according to the present invention.
Fig. 2 is a model of a zeroed GPS receiver of the present invention.
Fig. 3 is a schematic diagram of a seven-element uniform circular array according to the present invention.
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth by way of illustration only and are not intended to limit the scope of the invention.
As shown in fig. 1-3, an evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS includes the following steps:
step 1: determining an expression function x (nT) of the interference signal;
step 2: establishing an interference machine and a GPS receiver model;
and step 3: optimizing a flicker function r (nT) in the expression function x (nT) by adopting an evolutionary algorithm;
and 4, step 4: the degree of influence of the interfering signal on the receiver is determined.
In the step 1, the step of processing the raw material,
the complex field expression of the single interference signal x (nT) is shown as the following formula (1);
x(nT)=r(nT)[f1(nT)cos(2πfcnT+β)+jf2(nT)sin(2πfcnT+β)] (1)
in the formula (1), fc1575.42MHz, representing the GPS L1 signal carrier frequency; β represents an initial phase; t1/fsRepresents a sampling period; f. ofsRepresents the sampling frequency; f. of1(nT) and f2(nT) represents white gaussian noise after h (z) filtering; r (nt) represents a flicker function; n represents the sample point number, j represents the imaginary unit (j is the same as the common imaginary unit i);
fi(nT) (i ═ 1, 2) represents white gaussian noise after h (z) filtering.
The Z domain expression of H (Z) is:
Figure BDA0002829828740000041
the expression of the flicker function r (nt) is the following formula (2):
Figure BDA0002829828740000042
in the formula (2), function
Figure BDA0002829828740000043
Note Fa(0) 0. k is 0, 1, 2.; n denotes a sample point number.
In the step 2, the process of establishing the jammer and the GPS receiver model includes the following steps:
step 201: determining the arrangement mode of the jammer and the GPS receiver;
step 202: determining an expression function R (nT) of an array element receiving signal of a GPS receiver;
step 203: obtaining the optimal value W of the weight vector W by adopting an LMS algorithmopt
Step 204: calculating directional diagram
Figure BDA0002829828740000044
In the evolution band-limited Gaussian noise interference algorithm, the sequence a (k) is optimized through the evolution algorithm, and then the flicker function r (nT) is optimized, so that a GPS receiver obtains smaller null in the interference direction.
To optimize the a (k) sequence, jammer and GPS receiver models need to be built. A typical jammer and receiver layout is shown in figure 1. Jammers are arranged in different directions of the GPS receiver. Fig. 1 shows six jammers arranged at different azimuth angles. Suppose that the jammer has N totalJA GPS receiver receives NJAn interfering signal.
Suppose that the GPS receiver array element receives a signal R (nT) composed of NsIndividual GPS signal, NJAn interference signal and an additive white Gaussian noise N (nT).
The GPS antenna array selects more uniform circular arrays. Fig. 3 is a seven-array uniform circular array. Selecting the central array element as the reference array element, and assuming the reference array element to receive the signal x0(t) the expression is:
Figure BDA0002829828740000051
where A (t) is the signal amplitude,
Figure BDA0002829828740000052
is the initial phase of the signal.
The model structure of the zero-setting GPS receiver is shown in fig. 2, signals received by M antenna array elements obtain an optimal weight W through an LMS-based adaptive algorithm, and the signals obtained through weighting of the optimal weight W are demodulated through a GPS receiver to obtain GPS information. And calculating the directional diagram G through the optimal weight W. The directional null depth from the interference signal in the directional diagram is much larger than in the other directions. The influence degree of the interference signal on the receiver can be evaluated through the null depth of the directional diagram.
Assuming that M array elements except the array element at the center of circle in the uniform circular array of M array elements are respectively expressed as [ x ]1(t),x2(t),...,xM-1(t)]The propagation time delay between M-1 array elements and the reference array element is [ tau [ tau ] ]1,τ2,...,τM-1]. Considering that for the narrowband signal a (t- τ) ≈ a (t), the array element received signal may be denoted as X ═ Ax0(t) of (d). Wherein X ═ X0(t),x1(t),...,xM-1(t)]T
Figure BDA0002829828740000053
A is called a steering vector. Where j is an imaginary unit and ω represents angular frequency.
For the seven-array-element uniform circular array, the radius of the array element is half wavelength, namely r is lambda/2, and the guide vector of the seven-array element can be obtained
Figure BDA0002829828740000054
The respective expressions are as follows (3):
Figure BDA0002829828740000055
in the formula (3), theta is a signal pitch angle,
Figure BDA0002829828740000056
for azimuth of the signal, e.g. diagram3, respectively. j denotes an imaginary unit.
The GPS receiver array element receiving signal R (nT) is expressed as the following formula (4):
Figure BDA0002829828740000061
in the formula (4), the reaction mixture is,
Figure BDA0002829828740000062
a steering vector representing the ith GPS signal, AjA steering vector, S, representing the jth interfering signali(nT) denotes the received signal of the i-th GPS signal at the reference array element, Jj(nT) represents the received signal of the jth interference signal in the reference array element, and N (nT) represents additive white Gaussian noise. Here, j is a number of the interference signal, and is different from j in the imaginary unit expressed by the above equations (1) and (3). Since the units of imaginary numbers are commonly denoted by the letter i or j, they are also commonly used in summing functions
Figure BDA0002829828740000063
The manner of representation of (c). Therefore, in the case where the conventional expression is used for both cases, j represents a different meaning. But the person skilled in the art should easily recognize the representation of the serial number j and the representation of the imaginary unit j.
In step 203, an LMS algorithm is used to obtain an optimal value W of the weight vector WoptComprises the following steps:
step 2031: calculating weighted average signals y (k) of the signals of the other array elements except the reference array element;
step 2032: calculating an error signal e (k);
step 2033: the LMS algorithm obtains a weight vector W by utilizing a steepest descent method;
step 2034: obtaining optimal weight vector W by using LMS algorithmopt
The LMS algorithm minimizes the mean square error of the filtered output signal from the desired signal. In the null GPS, a reference array element output signal is used as an expected signal d (k), that is, d (k) x0(k) In that respect And weighted average y (k) of signals of the rest array elements except the reference array element is used as a filtering output signal.
y(k)=WTX(k) (5)
In formula (5), x (k) ═ x1(k),x2(k),...,xM-1(k)]TThe weight vector W is ═ W1,w2,...,wM-1]T
The error signal e (k) is calculated by the following equation (6).
e(k)=d(k)-y(k) (6)
In equation (6), the desired signal d (k) is the reference array element output signal, and y (k) is the filtered output signal.
The LMS algorithm obtains a weight vector by using the steepest descent method according to the formula (7) below.
Figure BDA0002829828740000064
In equation (7), μ represents an iteration step.
Figure BDA0002829828740000074
Representing the gradient and E the mathematical expectation. k is 0, 1, 2.
The LMS algorithm obtains the optimal weight Wopt=[W1 W2 ... WM-1]TThen, the array element weight can be represented by the following formula (8).
W=[1,-W1,-W2,...,-WM-1]T (8)
In the formula (8), M represents the number of array elements, and the number of array elements is 7 in fig. 3.
Directional diagram
Figure BDA0002829828740000071
Are represented by the following formula (9).
Figure BDA0002829828740000072
In the formula (9), abs is a function for calculating an absolute value.
In the step 3, the flicker function R (nT) in the expression function x (nT) is optimized by adopting an evolutionary algorithm, and the optimization is realized by optimizing a (k) sequence; the procedure for optimizing the a (k) sequence is as follows:
step 301: firstly, initializing parameters of an optimization algorithm;
step 302: for each population member, calculating a variation function updating gene of the gene mutation;
step 303: calculating a variation function updating gene of the cross variation;
step 304: calculating each member evaluation value scores;
step 305: the maximum value of all the evaluation values scores is calculated.
The evolution band-limited Gaussian noise interference algorithm optimizes the sequence a (k) through the evolution algorithm, and further optimizes the flicker function, so that the GPS receiver obtains smaller nulls in the interference direction.
1. Initialization: selecting a GPS signal number NsNumber of interfering signals NJLength of signal L, number of elements M of GPS receiver, initial signal parameter signal carrier frequency fcAnd a sampling frequency fs. Initializing the pitch angle theta and azimuth angle of each GPS signal and interference signal
Figure BDA0002829828740000073
And the weight W is initialized, the population member number N and the gene length N of the evolutionary algorithm parameters are initializedgInterval of genetic variation [ L, E]Maximum number of iterations NloopCrossover rate α, number of evaluation cycles NtInitializing the Gene a of each population Memberi(j,n)(1≤i≤NJ,1≤j≤NJ,1≤n≤Ng) The best assessment max _ score and corresponding gene max _ a, etc. are initialized.
The following steps 2-5 require a cycle of NloopNext, the process is carried out.
2. For each member, calculating a mutation function Mutate (a (j, n)) updating gene of gene mutation;
the gene update program is as follows:
for i=1:N
for j=1:NJ
k1=randi(1,Ng/2)*2-1
k2=randi(1,Ng/2)*2
αi(j,k1)=randi(L,E)
αi(j,k2)=randi(L,E)
end
end
wherein randi (A, B) represents a random integer within the interval [ A, B ].
3. Calculating a variation function Cross (a (j, n), alpha) updating gene of the Cross variation;
the calculation procedure is as follows:
for k=1:αN
i1=randi(1,N)
i2=randi(1,N)
i3=randi(1,Ng/2)
swap(αi1(1:NJ,i3:i3+Ng/2-1),αi2(1:NJ,i3:i3+Ng/2-1))
end
wherein swap (A, B) represents swaps A and B.
4. Calculating each member evaluation value scores ═ Score (a (j, n));
the calculation procedure is as follows:
for i=1:N
for j=1:Nt
calculating the received signal matrix R (nT) (see equation 4 above)
Taking R (nT) as a first action expected signal d (k), and taking the rest actions X (k)
Computing W by LMS Algorithmopt(calculated by the above equation 5, the above equation 6, the above equation 7)
Calculating directional diagram
Figure BDA0002829828740000081
(see equation 8 above, equation 9 above)
Calculating an interference direction average null
Figure BDA0002829828740000082
Wherein
Figure BDA0002829828740000083
(1≤j≤NJ) Are the interference signal pitch and azimuth angles.
end
score(i)=mean(s)
end
5. The scores maximum is calculated and the correspondence a (j, n) is recorded.
The calculation procedure is as follows:
[max,I]=max(scores)
if max_score<max
max_α=αI
end
wherein the Mutate function takes random values within the range of values for the random positions of a (k). The Cross function Cross swaps the members of the α N, swaps half of the data and swaps randomly. The Score function firstly calculates an interference signal, and then calculates an optimal weight W through an LMS algorithmoptThen return to the directional diagram
Figure BDA0002829828740000091
Depth of null in the direction of the interfering signal.
The performance difference between the band-limited Gaussian noise interference algorithm and the traditional band-limited Gaussian noise interference algorithm is compared and evolved through simulation experiments. Number of interference signals N in experimentjAnd respectively testing 1-6 samples, wherein the zero depth is the average value of 100 groups of sample data. The other parameters are: the array element number M is 7; the radius r of the antenna array is lambda/2; sampling rate fs3.2 GHz; the length of single sample data is 10000; the number Ns of GPS signals is 1; the pitch angle of the GPS signal is 70, and the azimuth angle is 10; the interference signal has a pitch angle of 70 and an azimuth angle of 7040, 70, 100, 150, 250, and 320; the S/N ratio is-30 dB; the dry-to-noise ratio J/N is 30 dB. The results of the experiment are shown in table 1.
TABLE 1 nulling depth contrast for bandlimited Gaussian noise interference Algorithm and Algorithm of the present invention
Number of disturbances Band limited gaussian noise Evolutionary algorithm Performance gain
1 48.85dB 43.07dB 5.78dB
2 48.66dB 42.97dB 5.69dB
3 50.11dB 44.52dB 5.59dB
4 47.60dB 41.94dB 5.66dB
5 41.59dB 32.83dB 8.76dB
6 29.81dB 20.04dB 9.77dB
From table 1, it can be seen that the null depth of the conventional band-limited gaussian noise algorithm in the interference direction is kept stable when the interference amount is 1-4, and the null depth is obviously improved when the interference amount is 5-6. The depth of null of the evolutionary band-limited gaussian noise interference algorithm has similar conclusions at different interference quantities. Compared with the traditional band-limited Gaussian noise interference algorithm, the performance gain of the evolved band-limited Gaussian noise interference algorithm is about 5-9 dB, and the performance difference is increased along with the increase of the interference quantity.
The zero setting GPS based on the LMS algorithm provides an interference algorithm for optimizing a flicker function through an evolutionary algorithm, and the performance of the algorithm is analyzed through the depth of interference zero-trap. The algorithm better disturbs the convergence of LMS through an optimized flicker function, obtains better performance gain, and has performance obviously superior to that of a band-limited Gaussian noise interference algorithm.
The invention optimizes the band-limited Gaussian noise interference flicker function through the evolutionary algorithm, thereby disturbing the convergence of the self-adaptive algorithm based on LMS in the GPS receiver, leading the average zero depth of the GPS receiver under the equal power to be lower and realizing better interference effect. Simulation experiment results show that the interference algorithm has more obvious interference effect than band-limited Gaussian noise interference when the interference amount is different.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (5)

1. An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS, comprising the steps of:
step 1: determining an expression function x (nT) of the interference signal;
step 2: establishing an interference machine and a GPS receiver model;
and step 3: optimizing a flicker function r (nT) in the expression function x (nT) by adopting an evolutionary algorithm;
and 4, step 4: the degree of influence of the interfering signal on the receiver is determined.
2. An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS according to claim 1, characterized in that: in the step 1, the step of processing the raw material,
the complex field expression of the single interference signal x (nT) is shown as the following formula (1);
x(nT)=r(nT)[f1(nT)cos(2πfcnT+β)+jf2(nT)sin(2πfcnT+β)] (1)
in the formula (1), fc=1575.42MHz, representing the GPS L1 signal carrier frequency; β represents an initial phase; t1/fsRepresents a sampling period; f. ofsRepresents the sampling frequency; f. of1(nT)) and f2(nT) represents white gaussian noise after h (z) filtering; r (nt) represents a flicker function; n denotes a sample number, and j denotes an imaginary unit.
3. An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS according to claim 1, characterized in that: in the step 2, the process of establishing the jammer and the GPS receiver model includes the following steps:
step 201: determining the arrangement mode of the jammer and the GPS receiver;
step 202: determining an expression function r (nT) of an array element receiving signal of a GPS receiver;
step 203: obtaining the optimal value of the weight vector W by adopting an LMS algorithm
Figure FDA0002829828730000012
Step 204: calculating directional diagram
Figure FDA0002829828730000011
4. An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS according to claim 3, characterized in that: in step 203, an LMS algorithm is used to obtain an optimal value W of the weight vector WoptComprises the following steps:
step 2031: calculating weighted average signals y (k) of the signals of the other array elements except the reference array element;
step 2032: calculating an error signal e (k);
step 2033: the LMS algorithm obtains a weight vector W by utilizing a steepest descent method;
step 2034: obtaining optimal weight vector W by using LMS algorithmopt
5. An evolutionary band-limited gaussian noise interference algorithm for LMS-based GPS according to claim 1, characterized in that: in the step 3, the flicker function r (nT) in the expression function x (nT) is optimized by adopting an evolutionary algorithm, and the optimization is realized by optimizing a (K) sequence; the procedure for optimizing the a (K) sequence is as follows:
step 301: firstly, initializing parameters of an optimization algorithm;
step 302: for each population member, calculating a variation function updating gene of the gene mutation;
step 303: calculating a variation function updating gene of the cross variation;
step 304: calculating each member evaluation value scores;
step 305: the maximum value of all the evaluation values scores is calculated.
CN202011468606.XA 2020-12-11 2020-12-11 GPS evolution band-limited Gaussian noise interference algorithm based on LMS Active CN112578410B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011468606.XA CN112578410B (en) 2020-12-11 2020-12-11 GPS evolution band-limited Gaussian noise interference algorithm based on LMS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011468606.XA CN112578410B (en) 2020-12-11 2020-12-11 GPS evolution band-limited Gaussian noise interference algorithm based on LMS

Publications (2)

Publication Number Publication Date
CN112578410A true CN112578410A (en) 2021-03-30
CN112578410B CN112578410B (en) 2023-09-01

Family

ID=75134956

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011468606.XA Active CN112578410B (en) 2020-12-11 2020-12-11 GPS evolution band-limited Gaussian noise interference algorithm based on LMS

Country Status (1)

Country Link
CN (1) CN112578410B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5940033A (en) * 1998-01-20 1999-08-17 The United States Of America As Represented By The Secretary Of The Army Apparatus, methods and computer program for evaluating multiple null forming antenna processors and jammers
US6590528B1 (en) * 2001-12-05 2003-07-08 Rockwell Collins, Inc. Low cost interference reduction system for GPS receivers
US20070127608A1 (en) * 2005-12-06 2007-06-07 Jacob Scheim Blind interference mitigation in a digital receiver
US20080246655A1 (en) * 2004-12-17 2008-10-09 Jon Olafur Winkel Spreading Codes for a Satellite Navigation System
US20100080268A1 (en) * 2008-09-23 2010-04-01 Catalin Lacatus Crosscorrelation interference mitigating position estimation systems and methods therefor
US20120139791A1 (en) * 2009-05-07 2012-06-07 Edelhaeusser Thorsten Concept for Determining an Estimated Value of a Location of a Receiving Element
US20170003394A1 (en) * 2015-07-03 2017-01-05 Thales Method for locating a jamming source jamming signals of a satellite navigation system and associated system
US20170111069A1 (en) * 2015-10-20 2017-04-20 The Aerospace Corporation Circuits and methods for reducing an interference signal that spectrally overlaps a desired signal
EP3188426A1 (en) * 2015-12-29 2017-07-05 Thales Method for controlling adaptive interference in a multi-channel receiver
CN107462872A (en) * 2017-07-03 2017-12-12 蔡绍滨 A kind of anti-major lobe suppression algorithm
US20190120969A1 (en) * 2016-11-10 2019-04-25 Cable Television Laboratories, Inc Systems and methods for interference detection in shared spectrum channels
CN111783301A (en) * 2020-07-02 2020-10-16 西北大学 Design method of anti-interference scattering communication system antenna array

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5940033A (en) * 1998-01-20 1999-08-17 The United States Of America As Represented By The Secretary Of The Army Apparatus, methods and computer program for evaluating multiple null forming antenna processors and jammers
US6590528B1 (en) * 2001-12-05 2003-07-08 Rockwell Collins, Inc. Low cost interference reduction system for GPS receivers
US20080246655A1 (en) * 2004-12-17 2008-10-09 Jon Olafur Winkel Spreading Codes for a Satellite Navigation System
US20070127608A1 (en) * 2005-12-06 2007-06-07 Jacob Scheim Blind interference mitigation in a digital receiver
US20100080268A1 (en) * 2008-09-23 2010-04-01 Catalin Lacatus Crosscorrelation interference mitigating position estimation systems and methods therefor
US20120139791A1 (en) * 2009-05-07 2012-06-07 Edelhaeusser Thorsten Concept for Determining an Estimated Value of a Location of a Receiving Element
US20170003394A1 (en) * 2015-07-03 2017-01-05 Thales Method for locating a jamming source jamming signals of a satellite navigation system and associated system
US20170111069A1 (en) * 2015-10-20 2017-04-20 The Aerospace Corporation Circuits and methods for reducing an interference signal that spectrally overlaps a desired signal
EP3188426A1 (en) * 2015-12-29 2017-07-05 Thales Method for controlling adaptive interference in a multi-channel receiver
US20190120969A1 (en) * 2016-11-10 2019-04-25 Cable Television Laboratories, Inc Systems and methods for interference detection in shared spectrum channels
CN107462872A (en) * 2017-07-03 2017-12-12 蔡绍滨 A kind of anti-major lobe suppression algorithm
CN111783301A (en) * 2020-07-02 2020-10-16 西北大学 Design method of anti-interference scattering communication system antenna array

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
BETZ J W 等: "Getting to M direct acquisition of the new military signal", GPS WORLD, no. 4, pages 40 - 46 *
周云锋;李宏;梁宏;廖万友;: "GPS接收阵列自适应抗干扰算法性能比较", 计算机仿真, no. 01, pages 69 - 71 *
李亚伟: "针对采用阵列天线的卫星导航接收机干扰与抗干扰技术研究", 《中国优秀硕士学位论文全文数据库-信息科技辑》, no. 4, pages 136 - 2539 *
李枫: "阵列天线的波束形成技术研究", 《中国优秀硕士学位论文全文数据库-信息科技辑》, no. 2, pages 136 - 140 *
毛虎;吴德伟;卢虎;白盟亮;: "GPSM码信号压制干扰样式效能分析", 电子科技大学学报, no. 03, pages 32 - 38 *

Also Published As

Publication number Publication date
CN112578410B (en) 2023-09-01

Similar Documents

Publication Publication Date Title
CN101349742B (en) Method for optimizing space between broad band phased array elements and measuring frequency and direction of frequency domain multiple targets
CN106772260B (en) Radar array and difference beam directional diagram optimization method based on convex optimized algorithm
CN111224706B (en) Surface multi-beam shaping method based on hybrid adaptive particle swarm algorithm
CN105445709B (en) A kind of thinning array near field passive location amplitude and phase error correction method
CN107944133B (en) Annular antenna array sparse method based on multi-target quantum spider swarm evolution mechanism
CN105842666B (en) Radar Subarray partition optimization method based on difference algorithm
CN106291541A (en) MIMO radar array design methodology based on non-homogeneous Subarray partition
CN108710103A (en) Strong and weak multiple target super-resolution direction finding based on thinned array and Sources number estimation method
CN109932696B (en) Digital beam forming method based on radar channel noise power
CN108508424A (en) A kind of Sidelobe based on antenna array error answers weight vector optimization method
CN112100701B (en) Two-dimensional distributed antenna subarray position optimization method based on genetic algorithm
CN110995331B (en) Beam forming method based on multipoint accurate control
CN112100811A (en) Antenna array directional diagram synthesis method based on adaptive wind-driven optimization algorithm
CN109212466B (en) Quantum dragonfly evolution mechanism-based broadband direction finding method
CN111273269A (en) IPSO-BP-based radar target positioning method of frequency diversity array
CN112578410A (en) Evolution band-limited Gaussian noise interference algorithm for GPS based on LMS
CN104346532B (en) MIMO (multiple-input multiple-output) radar dimension reduction self-adaptive wave beam forming method
CN104539331A (en) Array antenna beam forming method based on improved hybrid invasive weed optimization
CN116882149A (en) Antenna array synthesis method based on hybrid differential drosophila optimization algorithm
CN107918112A (en) Combined transceiving Beamforming Method based on the collaboration of active passive radar
CN111830495A (en) Airborne radar self-adaptive beam forming algorithm based on convex optimization learning
CN116933542A (en) Sparse linear array pattern synthesis method and system based on improved invasive weed algorithm
CN111695240B (en) Simulation modeling method and system for non-circular beam antenna pattern
CN114488027A (en) Wave beam zero setting and widening method of conformal array antenna
Zhao et al. An improved mixed-integer multi-objective particle swarm optimization and its application in antenna array design

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