CN102879790A - Anti-interference system and method based on digital beam forming and space-time zeroing cascade - Google Patents

Anti-interference system and method based on digital beam forming and space-time zeroing cascade Download PDF

Info

Publication number
CN102879790A
CN102879790A CN2011101957153A CN201110195715A CN102879790A CN 102879790 A CN102879790 A CN 102879790A CN 2011101957153 A CN2011101957153 A CN 2011101957153A CN 201110195715 A CN201110195715 A CN 201110195715A CN 102879790 A CN102879790 A CN 102879790A
Authority
CN
China
Prior art keywords
signal
interference
vector
satellite
space
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.)
Pending
Application number
CN2011101957153A
Other languages
Chinese (zh)
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.)
BEIJING TELLHOW SCI-TECH Co Ltd
Original Assignee
BEIJING TELLHOW SCI-TECH Co Ltd
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 BEIJING TELLHOW SCI-TECH Co Ltd filed Critical BEIJING TELLHOW SCI-TECH Co Ltd
Priority to CN2011101957153A priority Critical patent/CN102879790A/en
Publication of CN102879790A publication Critical patent/CN102879790A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention provides an anti-interference system and an anti-interference method based on digital beam forming and space-time zeroing cascade. The method is to cascade the digital beam forming and space-time zeroing, thus, the interference resistance of the whole system can be improved. The anti-interference system comprises an antenna array, a radio frequency module and an anti-interference processing module which are connected in sequence, wherein the anti-interference processing module comprises a digital beam forming module and a space-time zeroing processing module; the digital beam forming module is used for weighting and synthesizing various antenna array element signals into a reference beam signal pointed to a satellite, and synchronously obtaining an orthogonal vector signal orthogonal to a steering vector; output of the digital beam forming module is cascaded to the space-time zeroing processing module; and the space-time zeroing processing module is used for carrying out self-adaptive zeroing interference resistance processing to an interference signal in the reference beam signal, removing the interference signal in the reference beam signal, and outputting an effective satellite signal.

Description

Anti-interference system and method based on digital beam forming and space-time zeroing cascade
Technical Field
The invention relates to the technical field of satellite navigation, in particular to a satellite navigation anti-interference system and a method based on digital beam forming and space-time zeroing cascade.
Background
Currently, the Global Navigation Satellite System (GNSS) mainly includes GPS in the united states, GLONASS in russia, beidou in china, and GALILEO in europe. The signals of the Global Navigation Satellite System (GNSS) are generally interfered when being received. Taking the GPS in the united states as an example, if the interference power exceeds the GPS signal reception power by 24dB, the C/a code receiver of the commercial GPS cannot keep track of the signal. Experiments show that the jammer with the power of 1W can enable a C/A code receiver within 85 kilometers to be incapable of working.
Many contexts require high interference rejection of GNSS receivers. In a differential network base station of the GNSS, a receiver needs to have anti-interference performance; the GNSS time of the communication network base station and the mobile phone network base station is provided by a receiver with better anti-interference performance; ocean fleet activities are global, and it cannot be guaranteed that the ocean fleet does not work in unfriendly or even malicious environments. Near an aviation airport where airplanes approach each other, the management and scheduling of airports requires GNSS signals to be of higher accuracy and not threatened by interference. The GNSS receiver is required to have certain anti-interference capability in city control networks, geodetic surveying, precision engineering surveying, crustal motion monitoring, resource exploration and the like and special application related to safety departments.
Aiming at navigation signal interference, the technical system is divided into compression interference, deception interference and distributed three-dimensional interference. The pressing type interference is to press the signal at the front end of the GNSS receiver by transmitting an interference signal with a certain level through a transmitter, so that the receiver cannot receive satellite signals; the deceptive jamming is that a transmitter transmits a fake signal which has the same parameters as the GNSS signal but different information codes, so that a receiver receives wrong positioning information; the distributed three-dimensional interference is the omnibearing three-dimensional interference which is carried out by applying a plurality of various interference machines on the ground and in the air. The specific interference pattern is various, such as wideband gaussian noise, continuous wave, swept continuous wave, pulsed continuous wave, am continuous wave, narrowband/wideband fm signal, etc.
The development of the anti-jamming technology of the GNSS receiver mainly includes the following three technologies: 1. time domain and frequency domain filtering technology: this technique is implemented at a digital intermediate frequency. It uses DFT (discrete fourier transform) technique in the processing of digital intermediate frequency signals. Taking simple frequency domain amplitude processing as an example, if there is no radio frequency interference in the signal, the thermal noise power spectrum is fairly uniform in the frequency domain; if there is a narrow-band interference in the signal, it will have abnormal spectral lines in the frequency domain, and such abnormal spectral lines are adaptively filtered out before the DFT inversion. Mitre corporation developed a GPS narrowband interference frequency domain removal chip. 2. The spatial filtering technology comprises the following steps: an adaptive array nulling technique is applied to a GNSS receiver. The self-adaptive antenna array comprises a plurality of antenna array elements, each array element is connected with a microwave network, the microwave network is connected with a processor, the processor processes signals from the microwave network and then feeds back and adjusts the microwave network, the gain and the phase of each array element are controlled to change, and a zero point facing to an interference direction is generated in an antenna directional diagram, so that interference is counteracted. The number of the nulls is determined by the number of the antennas, and generally M array elements can control M-1 null points. Under an ideal condition, the adaptive antenna can improve the anti-interference performance of the GNSS receiver by 40-50 dB. 3. Space-time adaptive (STAP) techniques: compared with the simple time domain and frequency domain technology, the spatial domain filtering technology has obvious advantages, simple realization and small calculated amount, but has two defects: if the number of array units is M, the maximum number of nulls which can be generated by the array is M-1, and the maximum number of interferences which can be eliminated by the array is the same. This statement does not exclude the possibility that a null could eliminate both interferers, since the aperture of the array is of the order of the carrier wavelength, and the actual angular separation of the two interferers, although relatively large, is still small compared to the beam width. In practical application, the number of array elements is limited in consideration of constraints such as cost, size and power, so that the anti-interference performance of the adaptive array is reduced. Secondly, if the angle interval of one interference source to a certain GNSS satellite is very close, the signal of the GNSS satellite is attenuated to be unusable aiming at the null of the interference space. In practice, although the angular separation of the interference from the satellite is already relatively large, it is still relatively small compared to the beamwidth of the array. In order to eliminate the influence of different error sources, signals of as many GNSS satellites as possible are required. Aiming at the defects, a space-time adaptive (STAP) technology is provided on the basis of a space domain filtering technology, time domain delay is added on the space domain filtering, a space-time two-dimensional filter is constructed by utilizing a space domain array and the time domain delay, each stage of delay forms FIR filtering when each array element channel is seen in isolation, and interference can be removed in the time domain; from the same time delay node, different array elements form self-adaptive filtering of a space domain, a space interference source can be distinguished, and space domain null is formed to inhibit space domain interference; the STAP can be considered that the space-time domain filtering is popularized to the space-time domain or the time domain filtering is popularized to the space-time two-dimensional filtering, and the STAP has the capability of eliminating interference on a space-frequency two-dimensional plane. On the premise of not increasing array elements, the method greatly increases the degree of freedom of the array and qualitatively improves the anti-interference capability of narrow-band interference. However, the method has the disadvantages of large operation amount and difficult realization.
The Digital Beam Forming (DBF) technology in navigation signal processing is also a new technology which is developed currently, and the basic idea of the technology is to utilize the direction function product theorem of an antenna array, control the direction function of the antenna array by weighting on antenna array elements, compensate the echo phase difference caused by the propagation path difference caused by different sensors in spatial positions, realize in-phase superposition, and achieve the purpose of controlling an antenna array directional diagram to dynamically generate high-gain narrow beams in the useful signal direction. To form multiple beams, a plurality of beamformers may be used that direct steering vectors in different directions.
Disclosure of Invention
The invention provides a satellite navigation anti-interference system and a method based on digital beam forming and space-time zero-setting cascade connection, which cascade digital beam forming and space-time zero-setting processing and improve the anti-interference performance of the whole system.
The invention can also form a plurality of independent beams to carry out anti-interference processing, the working state of the multi-beam realizes flexible scanning of the beams and quick change of the waveform through continuous change of the weight, and each beam can deal with a plurality of interference directions by combining self-adaptive interference zero setting.
The technical scheme of the invention is as follows:
an anti-interference system based on digital beam forming and space-time zeroing cascade connection comprises an antenna array, a radio frequency module and an anti-interference processing module which are sequentially connected, and is characterized in that the anti-interference processing module comprises a digital beam forming module and a space-time zeroing processing module; the digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector obtained by the satellite pointing control information, so that the reference beam signal forms the maximum gain in the satellite direction, and meanwhile, an orthogonal vector signal orthogonal to the steering vector is obtained, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal.
The digital beam forming module comprises a guiding vector calculation module and an orthogonal vector calculation module, wherein the guiding vector calculation module calculates satellite position information according to the output of the inertial navigation system and ephemeris forecast, and simultaneously calculates the pitch angle and the azimuth angle of the satellite in a view field to obtain an N-dimensional guiding vector pointing to the satellite, so that the SINR of the designated satellite in the view field is maximized; the orthogonal vector calculation module extracts one dimension of the N-dimensional input signal as a reference beam signal, and simultaneously multiplies the N-dimensional input signal by an N-1-dimensional orthogonal matrix of the N-dimensional steering vector to obtain an N-1-dimensional orthogonal vector signal orthogonal to the steering vector.
The space-time zero-setting processing module comprises an adaptive zero-setting anti-interference weight calculation module, and the adaptive zero-setting anti-interference weight calculation module calculates an adaptive zero-setting weight through a random gradient LMS adaptive algorithm on an input N-1 dimensional orthogonal vector signal under a minimum mean square error criterion.
A joint anti-interference system comprising a plurality of independent digital beam forming and space-time zero-setting based cascades comprises an antenna array, a radio frequency module and an anti-interference processing module which are connected in sequence, and is characterized in that the anti-interference system comprises a plurality of independent anti-interference processing subsystems, and each anti-interference processing subsystem comprises the radio frequency module, the digital beam forming module and the space-time zero-setting processing module which are connected in sequence respectively; the digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector provided by the satellite pointing control information, so that the reference beam signal forms maximum gain in the satellite direction, an orthogonal vector signal orthogonal to the steering vector is obtained at the same time, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal.
An anti-interference method based on digital beam forming and space-time zero-setting cascade is characterized in that two technologies of digital beam forming and space-time zero-setting anti-interference are combined, through the digital beam forming technology, according to a steering vector obtained by satellite directional control information, each antenna array element signal is weighted and synthesized into a reference beam signal pointing to a satellite, so that the reference beam signal forms maximum gain in the satellite direction, an orthogonal vector signal orthogonal to the steering vector is obtained at the same time, and an interference signal in the reference beam signal is predicted; and then, the interference signals in the predicted reference beam signals are subjected to adaptive zero-adjusting anti-interference processing by a space-time zero-adjusting anti-interference technology, so that the interference signals in the reference beam signals are eliminated, and useful satellite signals are output.
The steering vector obtained according to the pointing control information of the satellite refers to: according to the output of the inertial navigation system and ephemeris forecast, calculating satellite position information, and simultaneously calculating a pitch angle and an azimuth angle of the satellite in the view field to obtain a guide vector pointing to the satellite, so that the SINR of the designated satellite in the view field is maximized; the steering vector is obtained by the following method:
assuming that the array element positions are written in vector form as
Figure BSA00000536923900041
The signal direction vector is
Figure BSA00000536923900042
If the modulus of the signal direction vector is 1, the spherical coordinate of the signal direction vector is
Figure BSA00000536923900043
Wherein,
Figure BSA00000536923900044
when the angle between the signal direction vector and the positive direction of the z axis is theta, the x axis rotates anticlockwise to reach the angle rotated by the projection of the signal direction vector on the xoy plane when viewed from the positive direction of the z axis, and then the rectangular coordinate of the signal direction vector is expressed as follows:
the wave path difference between the arrival of the signal at each array element and the origin is expressed as
Figure BSA00000536923900046
Dot multiplication of an array element vector and a signal vector; the expression from which the steering vector is derived is:
Figure BSA00000536923900047
the weighting and synthesizing of the antenna array element signals into a reference beam signal pointing to a satellite means that: extracting one dimension of the N-dimensional input signal as a reference beam signal according to a steering vector pointing to the satellite, wherein the one-dimension reference beam signal is obtained by the following method:
if the system receives the N-dimensional array signal:
X(θ)=(x0,x1,…,xN-1)
the maximum composite signal after digital beamforming is:
d=X(θ)AH(θ)
wherein A (theta) is a guide vector, AH(θ) represents the conjugate transpose of steering vector a (θ); d is the one-dimensional reference beam signal, and d includes all signals of a certain satellite;
the orthogonal vector signal orthogonal to the steering vector acquired at the same time is the remaining N-1 dimensional signal:
Y=X(θ)P(AH(θ))=X(θ)(I-AH(θ)A(θ)/(A(θ)AH(θ))) representing said N-1 orthogonal vector signals; wherein, PRepresents an orthogonal projection, and I represents an identity matrix; y does not contain any navigation signals.
The space-time zero-setting anti-interference processing comprises adaptive zero-setting anti-interference weight calculation, wherein the adaptive zero-setting anti-interference weight calculation refers to that an N-1-dimensional orthogonal vector signal is sent to an adaptive filter for processing, and the adaptive filter comprises two parts: the filtering process: calculating a response of the linear filter output to the input signal; an estimation error is then generated by comparing the output result field expected response. The self-adaptive process comprises the following steps: automatically adjusting filter parameters according to the estimate; these two processes work together to form a feedback loop; in the filtering process, the expected response d (n) is processed together with the input vector x (n), and given an input, the filter produces an outputAn estimate as an expected response d (n); defining an estimation error e (n) as the difference between the expected response and the actual filter output, adding the estimation error e (n) and an input vector x (n) to the adaptive control part, wherein a feedback loop surrounding the weight vector is closed loop, and the closed loop algorithm is the following according to a stochastic gradient LMS adaptive algorithm proposed by Widrow and Hoff:
Figure BSA00000536923900052
in the above formula, mu is the step size factor, w (n) is the real-time weight, d*Representing a conjugate operation.
The steering vector includes corrections and compensations for the amplitude and phase of the antenna, X (θ) — (X) for the received array signal0,x1,…,xN-1) After compensation correction, it is expressed as:
Xcorrect(θ)=(x0,x1/M1(θ,φ),…,xN-1/MN-1(θ,φ))
the above formula uniformly represents the correction and compensation of amplitude and phase; wherein M isiAnd (theta, phi) represents the amplitude and phase difference between each antenna element and the reference element obtained by measurement.
The antenna array comprises a non-uniform array or an nonlinear array, and the phase difference of the non-uniform array or the nonlinear array caused by the wave path difference is as follows: if the unit vector for a certain incoming wave traveling direction in space is
Figure BSA00000536923900053
I.e. the connecting line direction of the carrier and the satellite, and the wave path difference is the vector and the displacement vector r from each array element to the reference array elementiThe inner product of (2) is then the compensation parameter of each array element is
Way I: sin ( 2 π · s → · r i → λ )
and a path Q: cos ( 2 π · s → · r i → λ )
where λ is the wavelength.
The invention has the technical effects that:
the invention designs an anti-interference system based on digital beam forming and space-time zero-tuning cascade aiming at high requirements of anti-interference performance of a GNSS (Global Navigation Satellite System) receiver in many contexts, and the anti-interference system comprises an antenna array, a radio frequency module and an anti-interference processing module, can replace an antenna and a radio frequency part of a traditional GNSS receiver to be connected with a baseband part of a Navigation receiver, and can also replace the antenna part to be directly connected with a general GNSS Navigation receiver. The system cascades the digital beam forming and the space-time zero-setting treatment, combines the digital beam forming and the space-time zero-setting anti-interference, and improves the anti-interference performance of the whole system because the array gain of the digital beam forming is approximately equal to the number of array elements and the anti-interference capacity of N array elements can be improved by about N times; the invention also designs a combined anti-interference system comprising a plurality of independent anti-interference subsystems based on digital beam forming and space-time zero-setting cascade connection, which comprises a plurality of independent anti-interference subsystems based on digital beam forming and space-time zero-setting cascade connection, can simultaneously form a plurality of independent beams for anti-interference processing, each subsystem can carry out independent pointing control and self-adaptive anti-interference processing aiming at a certain satellite and can coordinate to work, the working state of the multi-beam realizes flexible scanning and rapid agility of the waveform of the beam through continuous weight value conversion, and meanwhile, the combined anti-interference zero-setting method can deal with a plurality of interference directions, has stronger protection capability on satellite signals in a complex interference scene, stronger satellite pointing capability and stronger robustness, and has correction capability on amplitude and phase inconsistency among receiving channels.
Drawings
Fig. 1 is a schematic overall structure diagram of an anti-interference system according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of an overall structure of an anti-interference system with an up-conversion universal interface according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a cascade relationship between a digital beam forming module and a space-time zeroing processing module of the anti-interference system of the present invention.
Fig. 4 is a schematic diagram of a plurality of independent joint anti-interference systems based on digital beam forming and space-time nulling cascade connection according to the present invention.
Fig. 5 is a diagram of an actual antenna gain pattern.
FIG. 6 is a schematic diagram of non-uniform matrix wave path difference calculation.
Fig. 7 is a schematic diagram of a space-time STAP filter structure.
FIG. 8 is a schematic diagram of an adaptive nulling LMS algorithm according to the present invention.
Detailed Description
The embodiments of the present invention will be described in further detail below with reference to the accompanying drawings.
As shown in fig. 1, which is a schematic structural diagram of an anti-interference system according to an embodiment of the present invention, the anti-interference system includes three main components, namely an antenna array 3 (note: the number of array elements 31 is variable, and a quaternary array is shown in the schematic diagram), a multi-channel down-conversion radio frequency module 4, and an anti-interference processing module 5; the output of the anti-jamming system is connected to the GNSS satellite navigation receiver baseband section 6. The multi-channel down-conversion radio frequency module 4 is responsible for completing the amplification and frequency conversion work of signals, converting the frequency of an input GNSS carrier frequency signal into a low intermediate frequency for digital signal processing, and amplifying the signals to a level suitable for sampling; the output digital intermediate frequency signal enters an anti-interference processing module 5 for anti-interference processing, and the output signal processed by the anti-interference processing module 5 is subjected to baseband processing and navigation positioning calculation. Another general connection method is shown in fig. 2, the rf section 4 includes a multi-channel down-conversion rf module 41 and an up-conversion rf module 42, which is responsible for re-converting the processed if signal to a carrier frequency to ensure the versatility of the interface with the GNSS receiver. Therefore, the anti-interference system designed by the invention can replace the connection between the antenna and the radio frequency part of the traditional GNSS receiver and the baseband part of the navigation receiver, and can also replace the direct connection between the antenna part and the general GNSS navigation receiver.
The anti-interference processing module 5 of the present invention cascades the digital beam forming module and the space-time zeroing processing module, as shown in fig. 3, a schematic diagram of the cascade relation between the digital beam forming module and the space-time zeroing processing module in the anti-interference processing module 5 of the anti-interference system of the present invention is shown. The anti-interference processing module comprises a digital beam forming module and a space-time zero-setting processing module; the digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector obtained by the satellite pointing control information, so that the reference beam signal forms the maximum gain in the satellite direction, and meanwhile, an orthogonal vector signal orthogonal to the steering vector is obtained, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal. The digital beam forming module comprises a guiding vector calculation module and an orthogonal vector calculation module, wherein the guiding vector calculation module calculates satellite position information according to the output of the inertial navigation system and ephemeris forecast, and simultaneously calculates the pitch angle and the azimuth angle of the satellite in a view field to obtain an N-dimensional guiding vector pointing to the satellite, so that the SINR of the designated satellite in the view field is maximized; the orthogonal vector calculation module extracts one dimension of the N-dimensional input signal as a reference beam signal, and simultaneously multiplies the N-dimensional input signal by an N-1-dimensional orthogonal matrix of the N-dimensional steering vector to obtain an N-1-dimensional orthogonal vector signal orthogonal to the steering vector. The space-time zero-setting processing module comprises an adaptive zero-setting anti-interference weight calculation module, and the adaptive zero-setting anti-interference weight calculation module calculates an adaptive zero-setting weight through a random gradient LMS adaptive algorithm on the input N-1 dimensional orthogonal vector signal under the minimum mean square error criterion.
As shown in fig. 4, an anti-interference system including multiple independent digital beam forming and space-time zeroing cascade connections includes an antenna array, a radio frequency module, and an anti-interference processing module, which are connected in sequence, where the anti-interference system includes multiple independent anti-interference processing subsystems, and each anti-interference processing subsystem includes a radio frequency module, a digital beam forming module, and a space-time zeroing processing module, which are connected in sequence; each subsystem can perform independent pointing control and self-adaptive anti-interference processing on a certain satellite, and can work coordinately to ensure that the system can effectively cover satellite signals. Fig. 4 is a schematic diagram showing a cascade relationship of a plurality of independent anti-interference subsystems of the anti-interference system of the present invention. Each anti-interference processing subsystem respectively comprises a multi-channel down-conversion radio frequency module DDC, a digital beam forming module DBF-1 and a space-time zero-setting processing module STAP-1, a multi-channel down-conversion radio frequency module DDC, a digital beam forming module DBF-2 and a space-time zero-setting processing module STAP-2 which are connected in sequence, a multi-channel down-conversion radio frequency module DDC, a digital beam forming module DBF-N and a space-time zero-setting processing module STAP-N, and intermediate frequency signals IF _1, IF _2, IF. The digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector provided by the satellite pointing control information, so that the reference beam signal forms the maximum gain in the satellite direction, and meanwhile, an orthogonal vector signal orthogonal to the steering vector is obtained, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal.
An anti-interference method based on digital beam and space-time zero-setting cascade is characterized in that two technologies of digital beam forming and space-time zero-setting anti-interference are combined, through the digital beam forming technology, according to a steering vector obtained by satellite directional control information, antenna array element signals are weighted and synthesized into a reference beam signal pointing to a satellite, so that the reference beam signal forms maximum gain in the satellite direction, an orthogonal vector signal orthogonal to the steering vector is obtained at the same time, and an interference signal in the reference beam signal is predicted; and then, performing further anti-interference processing on the interference signal in the predicted reference beam signal by a space-time zero-setting anti-interference technology, eliminating the interference signal in the reference beam signal and outputting a useful satellite signal.
The difference between digital beam forming and traditional adaptive nulling is that a phased array technology is adopted on a reference signal, a guide vector (namely an array flow pattern) of a satellite is utilized to enable the reference signal beam to point to the satellite, and due to practical physical limitations such as array scale, the reference beam also has side lobes in a non-satellite direction, and interference signals can enter from the side lobes, so that useful satellite signals and interference signals exist in the reference signal, and at the moment, the output of beam forming is cascaded with traditional adaptive nulling anti-interference so as to further improve the anti-interference performance. As shown in fig. 3, the N-dimensional input signal outputs a one-dimensional reference signal according to a steering vector provided by the pointing control information, and at the same time, the N-1-dimensional orthogonal projection can be obtained by the N-1-dimensional orthogonal matrix of the N-dimensional input signal and the steering vector, because of the orthogonality with the steering vector, the N-1-dimensional output signal does not include a satellite signal but includes an interference signal, so that cancellation processing of an interference signal in a reference beam can be performed by an adaptive nulling algorithm, and the output after the interference cancellation is a satellite signal in the beam pointing direction.
The N-dimensional steering vector obtained from the pointing control information is: generally, according to Inertial Navigation System (INS) output and ephemeris forecast, satellite position information is calculated, and a pitch angle and an azimuth angle of a satellite in a field of view are calculated, so that a steering vector pointing to the satellite can be obtained, and a beam forming module maximizes SINR of a specified satellite in the field of view according to the steering vector.
Digital beamforming requires the synthesis of the best signal in a certain direction, and this N-dimensional vector is the steering vector. In array signal processing, array flow patterns are generally represented by N-dimensional vectors
A ( α ) = ( 1 , e j φ 1 , e j φ 2 , . . , e j φ N - 1 )
The task of resisting disturbance is to find and
Figure BSA00000536923900082
and an orthogonal N-dimensional weight value vector W is combined with a pointing control module in the navigation processor to calculate satellite coordinates by using ephemeris and time, and the satellite coordinates are converted into a steering vector and output to a beam forming module and a preset zero anti-interference weight value calculation module.
Assuming that the array element positions are written in vector form as
Figure BSA00000536923900091
The signal direction vector is
Figure BSA00000536923900092
The signal direction vector has a modulus of 1, and the spherical coordinate is
Figure BSA00000536923900093
The included angle between the signal direction vector and the positive direction of the z axis is shown, and theta is the angle rotated by the x axis anticlockwise when the signal direction vector is viewed from the positive direction of the z axis to the angle rotated by the projection of the signal direction vector on the xoy plane. The rectangular coordinates of the signal bearing vector can be expressed as:
Figure BSA00000536923900094
the difference in the path length of the signal between the arrival of each array element and the origin can be expressed as
Figure BSA00000536923900095
Dot multiplication of an array element vector and a signal vector; the expression from which the steering vector is derived is:
Figure BSA00000536923900096
according to the output of an Inertial Navigation System (INS) and ephemeris forecast, satellite position information is calculated, the pitch angle and the azimuth angle of the satellite in the view field are calculated at the same time, and the pitch angle and the azimuth angle are substituted into the formula, so that a guide vector pointing to the satellite can be obtained, and the SINR of the specified satellite in the view field is maximized.
According to the above expression of the steering vector, if the one-dimensional linear array is taken as an example for analysis, the steering vector can be obtained as follows:
A ( θ ) = ( 1 , e - j 2 π · d · sin θ λ , e - j 2 π · d · sin θ λ * 2 , . . . , e - j 2 π · d · sin θ λ * ( N - 1 ) )
the system receives the N-dimensional array signals as follows:
X(θ)=(x0,x1,…,xN-1)
the maximum composite signal after digital beamforming is:
d=X(θ)AH(θ)
wherein A isH(θ) represents the conjugate transpose of the steering vector a (θ). d is the one-dimensional reference signal, and according to the beam forming principle, d includes all signals of a certain satellite.
The remaining N-1 dimensional signals are: y ═ X (θ) P(AH(θ))=X(θ)(I-AH(θ)A(θ)/(A(θ)AH(θ))) in any of the N-1 dimensions, where P isFor orthogonal projection, I represents an identity matrix, Y represents the orthogonal vector signal of the N-1 dimension, and since it is orthogonal projection, Y does not contain any navigation signal, d is interference with any component that can be predicted by Y,interference in d is achieved by eliminating the predicted portion of Y in d.
The space-time zero-setting anti-interference processing comprises adaptive zero-setting anti-interference weight calculation and adaptive zero-setting anti-interference weight calculation, and a space-time adaptive two-dimensional processing STAP technology is generally adopted. Space-time adaptive (STAP) technology utilizes a space domain array and time domain delay to construct a space-time two-dimensional filter. Each array element channel is isolated, each stage of time delay forms FIR filtering, and interference can be removed in a time domain; from the same time delay node, different array elements form self-adaptive filtering of a space domain, a space interference source can be distinguished, and space domain null suppression space domain interference is formed. The STAP can be considered that the space-time domain filtering is popularized to the space-time domain or the time domain filtering is popularized to the space-time two-dimensional filtering, and the STAP has the capability of eliminating interference on a space-frequency two-dimensional plane.
The number of array elements is set as M, and the number of time delay units is set as N. The time delay of each time delay unit is Ts. The weight w is an MN × 1-dimensional weight vector, and the two-dimensional filter structure is shown in fig. 7.
For the constrained optimization problem, the value obtained by calculating the weight vector through an LMS (Least-Mean-Square) algorithm under the minimum Mean Square error criterion represents an estimation, and when the iteration number tends to be infinite, the estimated expected value approaches a wiener solution for the generalized stationary process. The spatial nulling technique obtained based on the principle adopts a self-adaptive algorithm to correct the weighted value of the antenna, so that the gain and the phase of each array element are changed, zero points facing to the direction of an interference source are generated in a directional diagram of the antenna array, and the strong interference level can be effectively reduced to a low noise level, thereby achieving the purpose of resisting interference.
The LMS algorithm is a member of the family of stochastic gradient algorithms. The algorithm uses deterministic gradients in the stochastic input wiener filter recursive computation. One significant feature of the LMS algorithm is its simplicity. In addition, the method does not need to calculate related correlation functions and matrix inversion operation, and is easy to implement.
In general, the LMS includes two basic processes:
■ Filtering Process: firstly, calculating the response of the output of the linear filter to an input signal; an estimation error is then generated by comparing the output result field expected response.
■ adaptive procedure: the filter parameters are automatically adjusted according to the estimate.
The space-time zero-setting anti-interference processing is to send an N-1-dimensional orthogonal vector signal into an adaptive filter to perform adaptive zero-setting anti-interference weight calculation, and the adaptive filter comprises two parts: the filtering process: calculating a response of the linear filter output to the input signal; an estimation error is then generated by comparing the output result field expected response. The self-adaptive process comprises the following steps: automatically adjusting filter parameters according to the estimate; these two processes work together to form a feedback loop, as shown in fig. 8.
During the filtering process, the expected response d (n) participates in the processing together with the input vector y (n). In this case, given an input, the filter produces an output
Figure BSA00000536923900101
As an estimate of the expected response d (n). Thus, the estimation error e (n) can be defined as the difference between the desired response and the actual filter output. The estimation error e (n) and the input vector y (n) are both applied to the adaptive control part, so that the feedback loop around the weight vectors is closed loop. According to the random gradient LMS adaptive algorithm proposed by Widrow and Hoff, the closed-loop algorithm is as follows:
Figure BSA00000536923900102
in the above formula, mu is the step size factor, w (n) is the real-time weight, d*Representing a conjugate operation.
For a DBF system, calibration and compensation of the antenna amplitude and phase are key technologies for implementation. Since neither the actual rf nor the antenna is ideal, the rf may introduce a fixed phase shift, and the antenna elements have different gains and phases in different directions, an actual antenna gain pattern may be as shown in fig. 5. Besides the pointing information, the main control part of the system needs to provide compensation coefficients caused by channel and radio frequency imperfections and apply the compensation coefficients to the compensation correction of the N-dimensional steering vectors.
In order to correct the effect of antenna inconsistency, it is necessary to measure the amplitude and phase of the antenna, and to measure the gain and phase difference M between each antenna element and the reference elementi(theta, phi). The specific measurement method can be completed by a darkroom method or a specially designed receiver algorithm in an open place. For a received array signal X (theta) — (X)0,x1,…,xN-1) After compensation correction, it is expressed as:
Xcorrect(θ)=(x0,x1/M1(θ,φ),…,xN-1/MN-1(θ,φ))
the above equation represents the correction and compensation of amplitude and phase.
Because the array cannot be a linear array or even a uniform array due to the limitation of the conditions of a carrier and the like, the phase calculation of the array flow pattern is not necessarily a multiple relation, and the phase difference caused by the wave path difference needs to be obtained.
As shown in FIG. 6, the unit vector of the traveling direction of a certain incoming wave in space is
Figure BSA00000536923900111
(i.e. the direction of the connection between the carrier and the satellite), the wave path difference is the vector and the displacement vector r from each array element to the reference array elementiInner product of, compensation parameters of each array element are
Way I: sin ( 2 π · s → · r i → λ )
and a path Q: cos ( 2 π · s → · r i → λ )
where λ is the wavelength.
It should be noted that the above-mentioned embodiments enable a person skilled in the art to more fully understand the invention, without restricting it in any way. Therefore, although the present invention has been described in detail in this specification and examples, it will be understood by those skilled in the art that the present invention may be modified and equivalents may be substituted; all technical solutions and modifications thereof which do not depart from the spirit and scope of the present invention are intended to be covered by the scope of the present invention.

Claims (10)

1. An anti-interference system based on digital beam forming and space-time zeroing cascade connection comprises an antenna array, a radio frequency module and an anti-interference processing module which are sequentially connected, and is characterized in that the anti-interference processing module comprises a digital beam forming module and a space-time zeroing processing module; the digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector provided by the satellite pointing control information, so that the reference beam signal forms maximum gain in the satellite direction, an orthogonal vector signal orthogonal to the steering vector is obtained at the same time, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal.
2. The anti-jamming system based on digital beam forming and space-time zeroing cascade of claim 1, wherein the digital beam forming module comprises a steering vector calculation module and an orthogonal vector calculation module, the steering vector calculation module calculates satellite position information according to inertial navigation system output and ephemeris forecast, and simultaneously calculates a pitch angle and an azimuth angle of a satellite in a field of view to obtain a steering vector pointing to the satellite, so that SINR of the designated satellite in the field of view is maximized; the orthogonal vector calculation module extracts one dimension of the N-dimensional input signal as a reference beam signal, and simultaneously multiplies the N-dimensional input signal by an N-1-dimensional orthogonal matrix of the N-dimensional steering vector to obtain an N-1-dimensional orthogonal vector signal orthogonal to the steering vector.
3. The anti-jamming system based on digital beam forming and space-time zero-ing cascade of claim 2, wherein the space-time zero-ing processing module comprises an adaptive nulling anti-jamming weight calculation module, and the adaptive nulling anti-jamming weight calculation module calculates an adaptive nulling weight from an input N-1-dimensional orthogonal vector signal under a minimum mean square error criterion through a stochastic gradient LMS adaptive algorithm.
4. A joint anti-interference system comprising a plurality of independent digital beam forming and space-time zero-setting based cascades comprises an antenna array, a radio frequency module and an anti-interference processing module which are connected in sequence, and is characterized in that the anti-interference system comprises a plurality of independent anti-interference processing subsystems, and each anti-interference processing subsystem comprises the radio frequency module, the digital beam forming module and the space-time zero-setting processing module which are connected in sequence respectively; the digital beam forming module is used for weighting and synthesizing each antenna array element signal into a reference beam signal pointing to the satellite according to a steering vector provided by the satellite pointing control information, so that the reference beam signal forms maximum gain in the satellite direction, an orthogonal vector signal orthogonal to the steering vector is obtained at the same time, and an interference signal in the reference beam signal is predicted; the output of the digital beam forming module is cascaded with the space-time zero-setting processing module, and the space-time zero-setting processing module is used for receiving the reference beam signal output by the digital beam forming module and the orthogonal vector signal orthogonal to the guide vector, performing adaptive zero-setting anti-interference processing on the interference signal in the reference beam signal, eliminating the interference signal in the reference beam signal, and outputting a useful satellite signal.
5. An anti-interference method based on digital beam forming and space-time zero-setting cascade is characterized in that two technologies of digital beam forming and space-time zero-setting anti-interference are combined, through the digital beam forming technology, according to a guide vector provided by satellite direction control information, antenna array element signals are weighted and synthesized into a reference beam signal pointing to a satellite, the reference beam signal forms maximum gain in the satellite direction, meanwhile, an orthogonal vector signal orthogonal to the guide vector is obtained, and an interference signal in the reference beam signal is predicted; and then, the interference signals in the predicted reference beam signals are subjected to adaptive zero-adjusting anti-interference processing by a space-time zero-adjusting anti-interference technology, so that the interference signals in the reference beam signals are eliminated, and useful satellite signals are output.
6. The anti-jamming method based on digital beam forming and space-time zeroing cascade of claim 5, wherein the steering vector provided according to the pointing control information of the satellite is: according to the output of the inertial navigation system and ephemeris forecast, calculating satellite position information, and simultaneously calculating a pitch angle and an azimuth angle of the satellite in the view field to obtain a guide vector pointing to the satellite, so that the SINR of the designated satellite in the view field is maximized; the steering vector is obtained by the following method:
assuming that the array element positions are written in vector form as
Figure FSA00000536923800021
The signal direction vector is
Figure FSA00000536923800022
If the modulus of the signal direction vector is 1, the spherical coordinate of the signal direction vector is
Figure FSA00000536923800023
Wherein,when the angle between the signal direction vector and the positive direction of the z axis is theta, the x axis rotates anticlockwise to reach the angle rotated by the projection of the signal direction vector on the xoy plane when viewed from the positive direction of the z axis, and then the rectangular coordinate of the signal direction vector is expressed as follows:
Figure FSA00000536923800025
the wave path difference between the arrival of the signal at each array element and the origin is expressed as
Figure FSA00000536923800026
Dot multiplication of an array element vector and a signal vector; the expression from which the steering vector is derived is:
Figure FSA00000536923800027
7. the method of claim 6, wherein the directing the reference beam signal to the satellite is: extracting one dimension of the N-dimensional input signal as a reference beam signal according to a steering vector pointing to the satellite, wherein the one-dimension reference beam signal is obtained by the following method:
if the system receives the N-dimensional array signal:
X(θ)=(x0,x1,…,xN-1)
the maximum composite signal after digital beamforming is:
d=X(θ)AH(θ)
wherein A (theta) is a guide vector, AH(θ) represents the conjugate transpose of steering vector a (θ); d is the one-dimensional reference beam signal, and d includes all signals of a certain satellite;
the orthogonal vector signal orthogonal to the steering vector acquired at the same time is the remaining N-1 dimensional signal:
Y=X(θ)P(AH(θ))=X(θ)(I-AH(θ)A(θ)/(A(θ)AH(θ))) representing said N-1 orthogonal vector signals; wherein, P ^ represents orthogonal projection, I represents unit matrix; y does not contain any navigation signals.
8. The anti-interference method based on digital beam forming and space-time zeroing cascade of claim 7, wherein the space-time zeroing anti-interference process comprises adaptive zeroing anti-interference weight calculation, the adaptive zeroing anti-interference weight calculation refers to sending an N-1-dimensional orthogonal vector signal into an adaptive filter for processing, and the adaptive filter comprises two parts: the filtering process: calculating a response of the linear filter output to the input signal; an estimation error is then generated by comparing the output result field expected response. The self-adaptive process comprises the following steps: automatically adjusting filter parameters according to the estimate; these two processes work together to form a feedback loop; in the filtering process, the expected response d (n) is processed together with the input vector x (n), and given an input, the filter produces an output
Figure FSA00000536923800031
An estimate as an expected response d (n); defining the estimation error e (n) as desiredAdding an estimation error e (n) and an input vector x (n) to the adaptive control portion in response to the difference between the actual filter output and the estimated error e (n), a feedback loop around the weight vector is closed-loop according to a stochastic gradient LMS adaptive algorithm proposed by Widrow and Hoff, the closed-loop algorithm is as follows:
Figure FSA00000536923800032
in the above formula, mu is the step size factor, w (n) is the real-time weight, d*Representing a conjugate operation.
9. The method according to one of claims 1 to 8, wherein the steering vector comprises corrections and compensations of the antenna amplitude and phase, and X (θ) — (X) is used for the received array signal0,x1,…,xN-1) After compensation correction, it is expressed as:
Xcorrect(θ)=(x0,x1/M1(θ,φ),…,xN-1/MN-1(θ,φ))
the above formula uniformly represents the correction and compensation of amplitude and phase; wherein M isiAnd (theta, phi) represents the amplitude and phase difference between each antenna element and the reference element obtained by measurement.
10. The anti-jamming method based on digital beam forming and space-time nulling cascade of claims 1 to 8, wherein the antenna array comprises a non-uniform array or an amorphous array, and the phase difference of the non-uniform array or the amorphous array due to the wave path difference is: if the unit vector for a certain incoming wave traveling direction in space is
Figure FSA00000536923800033
I.e. the direction of the connection between the carrier and the satellite, the wave path difference is the vector and eachDisplacement vector r from array element to reference array elementiThe inner product of (2) is then the compensation parameter of each array element is
Way I: sin ( 2 π · s → · r i → λ )
and a path Q: cos ( 2 π · s → · r i → λ )
where λ is the wavelength.
CN2011101957153A 2011-07-13 2011-07-13 Anti-interference system and method based on digital beam forming and space-time zeroing cascade Pending CN102879790A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101957153A CN102879790A (en) 2011-07-13 2011-07-13 Anti-interference system and method based on digital beam forming and space-time zeroing cascade

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101957153A CN102879790A (en) 2011-07-13 2011-07-13 Anti-interference system and method based on digital beam forming and space-time zeroing cascade

Publications (1)

Publication Number Publication Date
CN102879790A true CN102879790A (en) 2013-01-16

Family

ID=47481184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101957153A Pending CN102879790A (en) 2011-07-13 2011-07-13 Anti-interference system and method based on digital beam forming and space-time zeroing cascade

Country Status (1)

Country Link
CN (1) CN102879790A (en)

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103105614A (en) * 2013-01-17 2013-05-15 陕西北斗恒通信息科技有限公司 Space and time domain joint anti-jamming method based on inertial navigation assisting
CN103323861A (en) * 2013-06-19 2013-09-25 电子科技大学 Method for improving steady-state performance of adaptive algorithm
CN103514638A (en) * 2012-06-15 2014-01-15 深圳市金溢科技有限公司 Positioning device based on DSRC technology, RSU, and application system
CN103630911A (en) * 2013-12-11 2014-03-12 北京北斗星通导航技术股份有限公司 Navigation information processing method and apparatus
CN103760575A (en) * 2014-01-15 2014-04-30 北京北斗星通导航技术股份有限公司 Anti-interference Beidou satellite navigation receiver board card and receiver terminal thereof
CN104880714A (en) * 2015-06-15 2015-09-02 中国人民解放军国防科学技术大学 GNSS antenna array anti-interference method based on two-stage filtering structure
CN105044660A (en) * 2015-06-18 2015-11-11 许昌学院 Rapid correction method of array channel inconsistency error
CN105158773A (en) * 2015-07-06 2015-12-16 北京航天长征飞行器研究所 Multi-beam directional anti-interference processing method for satellite navigation reception system
CN106329123A (en) * 2016-08-19 2017-01-11 西安空间无线电技术研究所 Nulling method of closed-loop nulling antenna for suppressing different forms of interference
CN106483536A (en) * 2015-08-26 2017-03-08 罗克韦尔柯林斯公司 Satellite-signal collection using the antenna possessing beamforming capabilities
CN106842263A (en) * 2017-01-21 2017-06-13 中国人民解放军信息工程大学 Based on the satellite navigation method that navigation signal space time processing is combined with vector tracking
CN106933121A (en) * 2015-12-30 2017-07-07 北京空间飞行器总体设计部 The spot beam anternma control method of task based access control spatial and temporal distributions characteristic
CN108196223A (en) * 2017-12-22 2018-06-22 电子科技大学 A kind of array df method based on STAP processing
CN109286421A (en) * 2018-09-13 2019-01-29 垣纬多媒体卫星通信(上海)有限公司 A kind of synthesis anti-interference method and its system handling satellite communication signals
CN109507698A (en) * 2018-09-28 2019-03-22 西南电子技术研究所(中国电子科技集团公司第十研究所) The anti-interference steering vector automatic correction system of satellite navigation
CN109633692A (en) * 2018-11-26 2019-04-16 西南电子技术研究所(中国电子科技集团公司第十研究所) GNSS navigation satellite signal anti-interference processing method
CN109725335A (en) * 2018-12-11 2019-05-07 上海无线电设备研究所 More star formation of the digital multiple beam methods in satellite navigation system
CN110108298A (en) * 2019-04-22 2019-08-09 北京控制工程研究所 A kind of front and back resolves fault-tolerance combined navigation method parallel
CN110366825A (en) * 2017-03-03 2019-10-22 英特尔Ip公司 It is eliminated using the burr that inverse burr injects
CN110429996A (en) * 2019-07-23 2019-11-08 熊军 A kind of joint multiple antennas calibrating installation for group delay and non-linear distortion
CN111147103A (en) * 2019-12-12 2020-05-12 南京理工大学 Airspace real-time self-adaptive anti-interference system aiming at short-time pulse interference
CN111211826A (en) * 2020-01-10 2020-05-29 中国人民解放军战略支援部队航天工程大学 Recursive structure beam forming method and device
CN111510171A (en) * 2020-03-27 2020-08-07 中国人民解放军海军工程大学 Anti-interference device and anti-interference method for time slot switching working mode
CN111540371A (en) * 2020-04-22 2020-08-14 深圳市友杰智新科技有限公司 Method and device for beamforming microphone array and computer equipment
CN111934728A (en) * 2020-06-16 2020-11-13 和芯星通科技(北京)有限公司 Digital multi-beam antenna processing method, device and equipment
CN111929646A (en) * 2020-08-14 2020-11-13 中国地质大学(北京) Beam scanning ground penetrating radar system and intelligent beam scanning detection method
CN112532308A (en) * 2020-12-09 2021-03-19 中国电子科技集团公司第五十四研究所 Anti-interference zero setting system
CN112737661A (en) * 2020-12-16 2021-04-30 中电科航空电子有限公司 Interference coordination method and system for airborne Beidou equipment and other aircraft iridium equipment
CN112769469A (en) * 2021-01-23 2021-05-07 成都振芯科技股份有限公司 Method and device for controlling operation array element number based on beam forming
CN113300087A (en) * 2021-04-21 2021-08-24 北京爱科迪通信技术股份有限公司 Multimode self-adaptive phased array satellite antenna
CN113852383A (en) * 2021-09-23 2021-12-28 上海航天电子通讯设备研究所 Anti-interference processing system for VHF frequency band burst signals
CN114460605A (en) * 2020-11-10 2022-05-10 西安开阳微电子有限公司 Anti-interference method for navigation receiver
CN114978257A (en) * 2022-05-24 2022-08-30 中国人民解放军国防科技大学 Adaptive beamforming satellite MIMO efficient system
CN115336185A (en) * 2022-06-30 2022-11-11 上海移远通信技术股份有限公司 Communication method and terminal device
CN116106939A (en) * 2023-01-19 2023-05-12 湖南大佳数据科技有限公司 Satellite navigation anti-interference test method, system, equipment and medium
CN116224406A (en) * 2023-04-23 2023-06-06 东南大学 Low-cost multi-antenna GNSS/INS ultra-tight combined multistage anti-interference system
CN116449398A (en) * 2023-04-10 2023-07-18 中国矿业大学 Self-adaptive anti-interference method for satellite navigation receiver in antenna array element mutual coupling environment
CN117452442A (en) * 2023-12-21 2024-01-26 中国人民解放军国防科技大学 High-precision satellite navigation anti-interference method based on space-frequency self-adaptive processing

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998016077A2 (en) * 1996-10-10 1998-04-16 Teratech Corporation Communication system using geographic position data
CN1897486A (en) * 2006-05-12 2007-01-17 西安电子科技大学 Satellite figured antenna anti-intrusion based on high-resolution algorithm
CN101949702A (en) * 2010-07-28 2011-01-19 北京泰豪联星技术有限公司 Quick self-testing method for GNSS PVT quality by using MEMS accelerometer

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998016077A2 (en) * 1996-10-10 1998-04-16 Teratech Corporation Communication system using geographic position data
CN1897486A (en) * 2006-05-12 2007-01-17 西安电子科技大学 Satellite figured antenna anti-intrusion based on high-resolution algorithm
CN101949702A (en) * 2010-07-28 2011-01-19 北京泰豪联星技术有限公司 Quick self-testing method for GNSS PVT quality by using MEMS accelerometer

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
段亚峰: "数字多波束GNSS抗干扰处理机实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 05, 15 May 2011 (2011-05-15), pages 61 - 68 *
项建弘: "基于空时自适应处理的GPS调零技术应用研究", 《中国期刊网全文数据库》, 14 March 2011 (2011-03-14) *

Cited By (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103514638B (en) * 2012-06-15 2016-08-03 深圳市金溢科技股份有限公司 A kind of positioner based on DSRC technology, RSU and application system
CN103514638A (en) * 2012-06-15 2014-01-15 深圳市金溢科技有限公司 Positioning device based on DSRC technology, RSU, and application system
CN103105614A (en) * 2013-01-17 2013-05-15 陕西北斗恒通信息科技有限公司 Space and time domain joint anti-jamming method based on inertial navigation assisting
CN103323861A (en) * 2013-06-19 2013-09-25 电子科技大学 Method for improving steady-state performance of adaptive algorithm
CN103323861B (en) * 2013-06-19 2015-06-10 电子科技大学 Method for improving steady-state performance of adaptive algorithm
CN103630911A (en) * 2013-12-11 2014-03-12 北京北斗星通导航技术股份有限公司 Navigation information processing method and apparatus
CN103760575A (en) * 2014-01-15 2014-04-30 北京北斗星通导航技术股份有限公司 Anti-interference Beidou satellite navigation receiver board card and receiver terminal thereof
CN104880714A (en) * 2015-06-15 2015-09-02 中国人民解放军国防科学技术大学 GNSS antenna array anti-interference method based on two-stage filtering structure
CN105044660A (en) * 2015-06-18 2015-11-11 许昌学院 Rapid correction method of array channel inconsistency error
CN105158773A (en) * 2015-07-06 2015-12-16 北京航天长征飞行器研究所 Multi-beam directional anti-interference processing method for satellite navigation reception system
CN105158773B (en) * 2015-07-06 2018-03-09 北京航天长征飞行器研究所 A kind of multi-beam for satellite navigation reception system points to anti-interference processing method
CN106483536A (en) * 2015-08-26 2017-03-08 罗克韦尔柯林斯公司 Satellite-signal collection using the antenna possessing beamforming capabilities
CN106483536B (en) * 2015-08-26 2019-04-02 罗克韦尔柯林斯公司 It is acquired using the satellite-signal for the antenna for having beamforming capabilities
CN106933121A (en) * 2015-12-30 2017-07-07 北京空间飞行器总体设计部 The spot beam anternma control method of task based access control spatial and temporal distributions characteristic
CN106933121B (en) * 2015-12-30 2019-04-30 北京空间飞行器总体设计部 The spot beam anternma control method of task based access control spatial and temporal distributions characteristic
CN106329123A (en) * 2016-08-19 2017-01-11 西安空间无线电技术研究所 Nulling method of closed-loop nulling antenna for suppressing different forms of interference
CN106329123B (en) * 2016-08-19 2019-03-22 西安空间无线电技术研究所 A kind of adjusting zero method of the closed loop Nulling antenna for inhibiting different form to interfere
CN106842263A (en) * 2017-01-21 2017-06-13 中国人民解放军信息工程大学 Based on the satellite navigation method that navigation signal space time processing is combined with vector tracking
CN110366825A (en) * 2017-03-03 2019-10-22 英特尔Ip公司 It is eliminated using the burr that inverse burr injects
CN110366825B (en) * 2017-03-03 2021-11-19 苹果公司 Glitch cancellation using inverse glitch implantation
CN108196223A (en) * 2017-12-22 2018-06-22 电子科技大学 A kind of array df method based on STAP processing
CN109286421A (en) * 2018-09-13 2019-01-29 垣纬多媒体卫星通信(上海)有限公司 A kind of synthesis anti-interference method and its system handling satellite communication signals
CN109286421B (en) * 2018-09-13 2021-03-26 上海垣信卫星科技有限公司 Comprehensive anti-interference method and system for processing satellite communication signals
CN109507698A (en) * 2018-09-28 2019-03-22 西南电子技术研究所(中国电子科技集团公司第十研究所) The anti-interference steering vector automatic correction system of satellite navigation
CN109507698B (en) * 2018-09-28 2022-07-08 西南电子技术研究所(中国电子科技集团公司第十研究所) Automatic correction system for anti-interference guide vector of satellite navigation
CN109633692A (en) * 2018-11-26 2019-04-16 西南电子技术研究所(中国电子科技集团公司第十研究所) GNSS navigation satellite signal anti-interference processing method
CN109633692B (en) * 2018-11-26 2022-07-08 西南电子技术研究所(中国电子科技集团公司第十研究所) GNSS navigation satellite signal anti-interference processing method
CN109725335A (en) * 2018-12-11 2019-05-07 上海无线电设备研究所 More star formation of the digital multiple beam methods in satellite navigation system
CN110108298A (en) * 2019-04-22 2019-08-09 北京控制工程研究所 A kind of front and back resolves fault-tolerance combined navigation method parallel
CN110429996A (en) * 2019-07-23 2019-11-08 熊军 A kind of joint multiple antennas calibrating installation for group delay and non-linear distortion
CN111147103A (en) * 2019-12-12 2020-05-12 南京理工大学 Airspace real-time self-adaptive anti-interference system aiming at short-time pulse interference
CN111147103B (en) * 2019-12-12 2021-09-21 南京理工大学 Airspace real-time self-adaptive anti-interference system aiming at short-time pulse interference
CN111211826A (en) * 2020-01-10 2020-05-29 中国人民解放军战略支援部队航天工程大学 Recursive structure beam forming method and device
CN111211826B (en) * 2020-01-10 2023-08-04 中国人民解放军战略支援部队航天工程大学 Recursive structure beam forming method and device
CN111510171A (en) * 2020-03-27 2020-08-07 中国人民解放军海军工程大学 Anti-interference device and anti-interference method for time slot switching working mode
CN111510171B (en) * 2020-03-27 2021-10-26 中国人民解放军海军工程大学 Anti-interference device and anti-interference method for time slot switching working mode
CN111540371B (en) * 2020-04-22 2020-11-03 深圳市友杰智新科技有限公司 Method and device for beamforming microphone array and computer equipment
CN111540371A (en) * 2020-04-22 2020-08-14 深圳市友杰智新科技有限公司 Method and device for beamforming microphone array and computer equipment
CN111934728B (en) * 2020-06-16 2022-02-11 北京北斗星通导航技术股份有限公司 Digital multi-beam antenna processing method, device and equipment
CN111934728A (en) * 2020-06-16 2020-11-13 和芯星通科技(北京)有限公司 Digital multi-beam antenna processing method, device and equipment
CN111929646A (en) * 2020-08-14 2020-11-13 中国地质大学(北京) Beam scanning ground penetrating radar system and intelligent beam scanning detection method
CN114460605A (en) * 2020-11-10 2022-05-10 西安开阳微电子有限公司 Anti-interference method for navigation receiver
CN112532308A (en) * 2020-12-09 2021-03-19 中国电子科技集团公司第五十四研究所 Anti-interference zero setting system
CN112737661A (en) * 2020-12-16 2021-04-30 中电科航空电子有限公司 Interference coordination method and system for airborne Beidou equipment and other aircraft iridium equipment
CN112769469B (en) * 2021-01-23 2023-02-24 成都振芯科技股份有限公司 Method and device for controlling and operating array element number based on beam forming
CN112769469A (en) * 2021-01-23 2021-05-07 成都振芯科技股份有限公司 Method and device for controlling operation array element number based on beam forming
CN113300087A (en) * 2021-04-21 2021-08-24 北京爱科迪通信技术股份有限公司 Multimode self-adaptive phased array satellite antenna
CN113852383A (en) * 2021-09-23 2021-12-28 上海航天电子通讯设备研究所 Anti-interference processing system for VHF frequency band burst signals
CN113852383B (en) * 2021-09-23 2022-08-09 上海航天电子通讯设备研究所 Anti-interference processing system for VHF frequency band burst signals
CN114978257A (en) * 2022-05-24 2022-08-30 中国人民解放军国防科技大学 Adaptive beamforming satellite MIMO efficient system
CN115336185A (en) * 2022-06-30 2022-11-11 上海移远通信技术股份有限公司 Communication method and terminal device
CN115336185B (en) * 2022-06-30 2024-04-23 上海移远通信技术股份有限公司 Communication method and terminal equipment
CN116106939A (en) * 2023-01-19 2023-05-12 湖南大佳数据科技有限公司 Satellite navigation anti-interference test method, system, equipment and medium
CN116449398A (en) * 2023-04-10 2023-07-18 中国矿业大学 Self-adaptive anti-interference method for satellite navigation receiver in antenna array element mutual coupling environment
CN116449398B (en) * 2023-04-10 2023-11-03 中国矿业大学 Self-adaptive anti-interference method for satellite navigation receiver in antenna array element mutual coupling environment
CN116224406A (en) * 2023-04-23 2023-06-06 东南大学 Low-cost multi-antenna GNSS/INS ultra-tight combined multistage anti-interference system
CN116224406B (en) * 2023-04-23 2024-04-26 东南大学 Low-cost multi-antenna GNSS/INS ultra-tight combined multistage anti-interference system
CN117452442A (en) * 2023-12-21 2024-01-26 中国人民解放军国防科技大学 High-precision satellite navigation anti-interference method based on space-frequency self-adaptive processing
CN117452442B (en) * 2023-12-21 2024-04-05 中国人民解放军国防科技大学 High-precision satellite navigation anti-interference method based on space-frequency self-adaptive processing

Similar Documents

Publication Publication Date Title
CN102879790A (en) Anti-interference system and method based on digital beam forming and space-time zeroing cascade
Fernández-Prades et al. Robust GNSS receivers by array signal processing: Theory and implementation
CN110412620B (en) Anti-interference antenna signal processing device
KR101708381B1 (en) Anti-jamming signal processing apparatus and method using a dual-polarized antenna array
Vagle et al. Performance analysis of GNSS multipath mitigation using antenna arrays
CN101900819B (en) Blind adaptive space-time array anti-interference method for navigation receiver
Li et al. A robust anti-jamming navigation receiver with antenna array and GPS/SINS
CN110850445B (en) Pulse interference suppression method based on space-time sampling covariance inversion
US10705176B2 (en) Signal direction processing for an antenna array
CN102033227A (en) Weak target detection method for passive radar taking global positioning system (GPS) navigation satellite as external radiation source
Chiba et al. Digital beam forming (DBF) antenna system for mobile communications
JP2006270847A (en) Antenna device
Daneshmand GNSS interference mitigation using antenna array processing
Wang et al. Design of optimum sparse array for robust MVDR beamforming against DOA mismatch
Jia et al. Adaptive blind anti-jamming algorithm using acquisition information to reduce the carrier phase bias
Azimi-Sadjadi et al. Wideband DOA estimation algorithms for multiple moving sources using unattended acoustic sensors
CN104391305A (en) Satellite navigation deception interference suppression method based on deception interference DOA estimation
Bao et al. Spoofing mitigation in Global Positioning System based on C/A code self-coherence with array signal processing
Kalyanaraman et al. GPS adaptive array phase compensation using a software radio architecture
CN106338742B (en) Dimensionality reduction self-adaptive multiple-beam gps signal anti-interference method based on cross-spectrum criterion
CN109116387A (en) Navigational satellite reflected signal intermediate frequency Enhancement Method
Zheng Adaptive antenna array processing for GPS receivers.
CN105353386A (en) Anti-interference method and device for navigation receiver through employing inertial navigation equipment
Daneshmand et al. GNSS interference and multipath suppression using array antenna
Cui et al. Adaptive antenna array based non-cooperative co-channel interference cancellation system for ultrashort wave radios and non-ideality analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 100085, Beijing, Haidian District on the West Road, No. 8 hospital (on the floor of science and technology building), building 4, East 701 room

Applicant after: BEIJING TELLHOW SCI-TECH CO., LTD.

Address before: 100083 Beijing city Haidian District Wangzhuang Road No. 1, Tsinghua Tongfang Technology Plaza B block, room 905

Applicant before: Beijing Tellhow Sci-tech Co., Ltd.

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: BEIJING TELLHOW LIANXING SCI-TECH CO., LTD. TO: BEIJING BEIFANG LIANXING TECHNOLOGY CO., LTD.

C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130116