CN109245781A - Anti-interference method when a kind of annulus conformal array sky - Google Patents
Anti-interference method when a kind of annulus conformal array sky Download PDFInfo
- Publication number
- CN109245781A CN109245781A CN201811088759.4A CN201811088759A CN109245781A CN 109245781 A CN109245781 A CN 109245781A CN 201811088759 A CN201811088759 A CN 201811088759A CN 109245781 A CN109245781 A CN 109245781A
- Authority
- CN
- China
- Prior art keywords
- interference
- array element
- array
- signal
- annulus
- 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
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/12—Neutralising, balancing, or compensation arrangements
- H04B1/123—Neutralising, balancing, or compensation arrangements using adaptive balancing or compensation means
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
The invention discloses anti-interference methods when a kind of sky of annulus conformal array, improve the robustness and rapidity of anti-interference antenna.To conformal when the aerial array of carrier surface does anti-interference process, since carrier has certain occlusion effect, " dark space " can be generated to the signal of certain incident angles, cause part array element that can not receive, adaptive array element selects the very good solution problem, when interfering more, array element and delay unit can be increase accordingly, autocorrelation matrix order increased dramatically, matrix inversion becomes not realizing in engineering, Anti-Jamming Technique avoids autocorrelation matrix and directly inverts when sky based on Eigenvalues Decomposition, reduces difficulty in computation.The present invention is based on anti-interference methods when the sky of the selection of adaptive array element and Eigenvalues Decomposition, realize interference free performance of the antenna in complex environment, flexibility is good, and calculation amount is small, it is easy to accomplish.
Description
Technical field
The present invention relates to anti-interference method when the sky of the selection of adaptive array element and Eigenvalues Decomposition, specifically a kind of annulus is total
Anti-interference method when shape battle array sky.
Background technique
In recent years, satellite navigation flourishes, and has either all obtained in military or civil field, satellite navigation system
It is widely applied, importance increasingly highlights.But since satellite distribution is in tens of thousands of kilometers from the ground of high-altitude, satellite-signal is being propagated
There is greatly decaying on the way, the power for reaching ground is typically well below ambient noise, therefore is highly prone to various interference signals
Influence, cause the navigation neceiver performance of day line back end sharply to decline, under adverse circumstances, or even be unable to complete positioning calculation.
Therefore it need to increase anti-interference process to front end antenna and effectively inhibit interference signal.
According to the situation of structuring the formation of actual antennas array element, linear array, face battle array and conformal array can be divided into.Most conformal arrays by
In blocking for carrier, only part array element can effectively receive signal, if still participating in calculating with whole array elements, it will result is caused to go out
Existing deviation, null may not be able to be directed at interference radiating way, therefore must carry out adaptive array element selection.In addition, antenna applications environment is multiple
Miscellaneous, interference source number is more, especially introduces time delay unit again in the case where more array elements, and difficulty in computation increased dramatically,
Autocorrelation matrix in treatment process is inverted extremely difficult in practical projects, it is therefore desirable in the feelings for not influencing performance as far as possible
Conversion process is carried out under condition, reduces difficulty in computation, is convenient for Project Realization.This patent is based on the selection of adaptive array element and characteristic value point
The anti-interference method when sky of solution, the effective solution above problem also achieve good result in engineering in practice, and antenna is anti-dry
The performance and robustness and rapidity disturbed all greatly improve.
It is rarer for the adaptive array element selection of annulus conformal array antenna and in existing scheme, certain angles are entered
Penetrating signal may cause interference free performance to decline in the presence of blocking；Autocorrelation matrix, which generallys use, directly inverts or dimension reduction method,
Direct inversion calculation amount is big, is not easy to realize；Dimension reduction method inevitably loses precision, will lead to interference free performance decline；
A large amount of operation is to hardware requirement height, and the increase of calculation amount will definitely increase the power consumption of anti-interference antenna, and small at this stage
Type, the development trend of low-power consumption are not inconsistent.
Summary of the invention
The purpose of the present invention is to provide anti-interference methods when a kind of annulus conformal array sky, to solve in above-mentioned background technique
The problem of proposition.
To achieve the above object, the invention provides the following technical scheme:
Anti-interference method when a kind of annulus conformal array sky, part array element caused by being blocked based on carrier cannot effectively receive letter
Number, so as to cause interference free performance decline, by the adaptively selected of array element, effective array element is determined, in later period Anti-interference algorithm
In processing, using the method for Eigenvalues Decomposition, the difficulty of inverting of higher dimensional matrix is reduced.
As a further solution of the present invention: by the azimuth information of incoming signal, the effective array element of adaptive selection is mentioned
High anti-interfering performance.
As a further solution of the present invention: the inverse matrix of autocorrelation matrix is not obtained by the method directly inverted, and
It is the method based on Eigenvalues Decomposition, improves computational efficiency.
As a further solution of the present invention: not changing original Anti-interference algorithm framework, need to only increase and increase certainly in front end
Adapt to array element selecting module, compatible original Anti-interference algorithm.
As a further solution of the present invention: the plane projection relationship based on array element and signal, array element 1 are reference array element,
When signal is with 0 ° < P ° < 90 ° of azimuth incidence, the only the 1st, 2,3,4 array elements can receive signal, remaining array element is due to being hidden
It keeps off and can not effectively receive, the matrix inversion based on Eigenvalues Decomposition is embodied in Anti-interference algorithm, Anti-interference algorithm process master
It comprises the steps of:
1. being selected as algorithm by adaptive array element provides effective information；
2. Signal Pretreatment: by received signal shift frequency to base band, realizing that IQ is decomposed, and by digital low-pass filtering to band
External signal is inhibited；
3. calculating autocorrelation matrix: being calculated by the accumulation of certain sample points；
4. autocorrelation matrix is inverted: can effectively avoid directly inverting for autocorrelation matrix, and precision by Eigenvalues Decomposition
Height is close with directly inverting；
5. automatically adjusting the amplitude and phase of each signal according to weight, make to form the zero point got deeply stuck in interference radiating way, in the phase
It hopes in sense and forms gain.
Compared with prior art, the beneficial effects of the present invention are: it is dry the invention proposes resisting when a kind of annulus conformal array sky
Method is disturbed, is selected by adaptive array element, can effectively solve the problem that the problem of interference free performance caused by being blocked due to carrier is declined,
System robustness is improved, is provided a strong guarantee for rear end satellite navigation receiver normal work；Pass through the feature of autocorrelation matrix
Value is decomposed, and directly inverting for higher dimensional matrix is avoided, and in the range of loss of significance allows, is reduced difficulty in computation, is promoted operation
Speed；And anti-interference method when based on the sky of the selection of adaptive array element and Eigenvalues Decomposition, it need to only increase in the front end of original frame
Add the adaptively selected module of array element, changes smaller, good compatibility；Design difficulty is low in engineer application, small power consumption, wants to hardware
It asks not high, meets the development trend of current miniaturization, low-power consumption.
Detailed description of the invention
The array element distribution and incoming signal plane projection structural representation of anti-interference method when Fig. 1 is a kind of annulus conformal array sky
Figure.
Anti-interference algorithm process flow diagram when the annulus conformal array sky of anti-interference method when Fig. 2 is a kind of annulus conformal array sky
Structural schematic diagram.
The adaptive array element choosing principles structural schematic diagram of anti-interference method when Fig. 3 is a kind of annulus conformal array sky.
Anti-interference method structure when the sky based on Eigenvalues Decomposition of anti-interference method when Fig. 4 is a kind of annulus conformal array sky
Schematic diagram.
Specific embodiment
Below in conjunction with the embodiment of the present invention, technical scheme in the embodiment of the invention is clearly and completely described,
Obviously, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based in the present invention
Embodiment, every other embodiment obtained by those of ordinary skill in the art without making creative efforts, all
Belong to the scope of protection of the invention.
In the embodiment of the present invention, anti-interference method when a kind of annulus conformal array sky, the plane projection based on array element and signal
Relationship, array element 1 is reference array element, when signal is with 0 ° < P ° < 90 ° of azimuth incidence, the only the 1st, 2,3,4 array elements can receive
Signal, remaining array element can not be received effectively due to being blocked, and the matrix inversion based on Eigenvalues Decomposition is embodied in anti-interference calculation
In method, Anti-interference algorithm process is mainly comprised the steps of:
(1) algorithm is selected as by adaptive array element and effective information is provided；
(2) Signal Pretreatment: by received signal shift frequency to base band, realize that IQ is decomposed, and pass through digital low-pass filtering pair
Out of band signal is inhibited；
(3) it calculates autocorrelation matrix: being calculated by the accumulation of certain sample points；
(4) autocorrelation matrix is inverted: can effectively avoid directly inverting for autocorrelation matrix, and precision by Eigenvalues Decomposition
Height is close with directly inverting；
(5) amplitude and phase that each signal is automatically adjusted according to weight made to form the zero point got deeply stuck in interference radiating way, in the phase
It hopes in sense and forms gain.
Complete technical solution:
For without loss of generality, it is assumed that the deflection of incoming signal is that be straight line AB in the plane vertical with incoming signal, then
Straight line AB includes that the array element of incoming signal side is effective array element, i.e. with quantitative mathematical relationship expression: the folder with reference array element
Angular region array element be receive the effective array element of signal；
Anti-interference method is to carry out Eigenvalues Decomposition to autocorrelation matrix when sky based on Eigenvalues Decomposition, utilizes auto-correlation
The Hermitian matrix properties of matrix obtain the inverse matrix of autocorrelation matrix by characteristic value and feature vector, so as to avoid
High level matrix is directly inverted；
Autocorrelation matrix RxFor non-singular matrix, Eigenvalues Decomposition is carried out to it:
Wherein, λiIndicate RxIth feature value, qiFor RxIth feature vector, Λ is characterized value λiWhat is constituted is diagonal
Matrix, Q are individual features vector qiThe unitary matrice of composition.Due to RxIt is Hermitian matrix, according to the spy of Hermitian matrix
Property it is available:
According to the characteristic of Hermitian matrix, JACOBI method is can be used in the Eigenvalues Decomposition of autocorrelation matrix.JACOBI
Method asks Hermitian matrix to obtain characteristic value and feature vector using unitary transformation, is by not during computer solving
Disconnected iteration and the process solved, when the element on Λ off-diagonal is close to 0, the element on Λ diagonal line is approximately RxSpy
Value indicative, therefore suitable threshold value is set and can guarantee precision, i.e., it is acquired by (2) formulaFor RxTrue inverse matrix.
In complex field, planar rotation matrix are as follows:
Wherein, sin θ > 0, cos θ > 0, θ are rotation angle, θ1+θ4=θ2+θ3.It is computedThat is QijIt is the tenth of the twelve Earthly Branches
Matrix, andIt is also the spin matrix of complex plane.
If R ∈ Cn×nIt is Hermitian matrix, remembers R0=R makees a series of rotation similarity transformations to R, i.e.,
It is found that RkIt is still Hermitian matrix.QkRk-1Only change Rk-1The i-th row, jth row element,
Only change Rk-1The i-th row, j row, i column, j column element.
It can be to the R after rotation transformation according to matrix multiplication principlek, RkThe expression-form of each element:
According to matrix theory, the diagonal entry of complex conjugate symmetry matrix is real number, and characteristic value is also real number, so Rk
Diagonal entryAlso it should be real number, thus,
θ4-θ3=θ2-θ1 (7)
According to the thought of JACOBI algorithm, off diagonal elementHave:
Simultaneous (7), (8) can obtain:
θ2-θ4=θ1-θ3=angle (rij) (9)
rijAnd rjiIt is a pair of of conjugate complex number, can respectively indicates are as follows:
Joint type (8), (9), (10), abbreviation can obtain:
sinθcosθ(rjj-rii)+|rij|cos2θ-|rij|sin2θ=0 (11)
It further calculates, has:
Formula (12) can be according to rij,rii,rjjUniquely determine the angle θ, the θ in formula (9)1,θ2,θ3,θ4Four unknown numbers cannot obtain
Unique solution out takes θ for convenience of calculating2=θ4=-θ1=-θ3, i.e. θ1=∠ rij/ 2, complex plane spin matrix QijIt can be rewritten as:
To which formula (5), (6) become:
It can be seen that from above derive when meeting formula (16),Zero can be melted into.
θ, θ are calculated according to formula (16)1, then obtain complex plane spin matrix Qij, choose suitable iteration order and thresholding
Value, is not difficult to obtain more actually RxCharacteristic value and feature vector.
It is obvious to a person skilled in the art that invention is not limited to the details of the above exemplary embodiments, Er Qie
In the case where without departing substantially from spirit or essential attributes of the invention, the present invention can be realized in other specific forms.Therefore, no matter
From the point of view of which point, the present embodiments are to be considered as illustrative and not restrictive, and the scope of the present invention is by appended power
Benefit requires rather than above description limits, it is intended that all by what is fallen within the meaning and scope of the equivalent elements of the claims
Variation is included within the present invention.
In addition, it should be understood that although this specification is described in terms of embodiments, but not each embodiment is only wrapped
Containing an independent technical solution, this description of the specification is merely for the sake of clarity, and those skilled in the art should
It considers the specification as a whole, the technical solutions in the various embodiments may also be suitably combined, forms those skilled in the art
The other embodiments being understood that.
Claims (5)
1. anti-interference method when a kind of annulus conformal array sky, which is characterized in that the method blocks caused part based on carrier
Array element cannot effectively receive signal, so as to cause interference free performance decline, by the adaptively selected of array element, determine effective battle array
Member, using the method for Eigenvalues Decomposition, reduces the difficulty of inverting of higher dimensional matrix in the processing of later period Anti-interference algorithm.
2. anti-interference method when a kind of annulus conformal array sky according to claim 1, which is characterized in that the method is led to
The azimuth information of incoming signal is crossed, the effective array element of adaptive selection improves interference free performance.
3. anti-interference method when a kind of annulus conformal array sky according to claim 1, which is characterized in that the method, no
The inverse matrix of autocorrelation matrix is obtained by the method directly inverted, but the method based on Eigenvalues Decomposition, raising calculate effect
Rate.
4. anti-interference method when a kind of annulus conformal array sky according to claim 1, which is characterized in that the method, no
Change original Anti-interference algorithm framework, need to only increase and increase adaptive array element selecting module, compatible original anti-interference calculation in front end
Method.
5. anti-interference method when a kind of annulus conformal array sky according to claim 1, which is characterized in that the method, base
In the plane projection relationship of array element and signal, array element 1 is reference array element, when signal is with 0 ° < P ° < 90 ° of azimuth incidence, only
1st, 2,3,4 array elements can receive signal, remaining array element can not be received effectively due to being blocked, the square based on Eigenvalues Decomposition
Battle array, which is inverted, to be embodied in Anti-interference algorithm, and Anti-interference algorithm process mainly comprises the steps of:
(1) algorithm is selected as by adaptive array element and effective information is provided；
(2) Signal Pretreatment: by received signal shift frequency to base band, realize that IQ is decomposed, and outer to band by digital low-pass filtering
Signal is inhibited；
(3) it calculates autocorrelation matrix: being calculated by the accumulation of certain sample points；
(4) autocorrelation matrix is inverted: directly inverting for autocorrelation matrix can effectively avoid by Eigenvalues Decomposition, and precision is high,
It is close with directly inverting；
(5) amplitude and phase that each signal is automatically adjusted according to weight make to form the zero point got deeply stuck in interference radiating way, believe in expectation
Number side is upwardly formed gain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811088759.4A CN109245781A (en) | 2018-09-18 | 2018-09-18 | Anti-interference method when a kind of annulus conformal array sky |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811088759.4A CN109245781A (en) | 2018-09-18 | 2018-09-18 | Anti-interference method when a kind of annulus conformal array sky |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109245781A true CN109245781A (en) | 2019-01-18 |
Family
ID=65058251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811088759.4A Pending CN109245781A (en) | 2018-09-18 | 2018-09-18 | Anti-interference method when a kind of annulus conformal array sky |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109245781A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104880714A (en) * | 2015-06-15 | 2015-09-02 | 中国人民解放军国防科学技术大学 | GNSS antenna array anti-interference method based on two-stage filtering structure |
| CN106940689A (en) * | 2017-03-07 | 2017-07-11 | 电子科技大学 | High-precision Eigenvalue Decomposition implementation method based on Jacobi iterative algorithms |
| US20170201354A1 (en) * | 2010-05-28 | 2017-07-13 | Cohere Technologies, Inc. | Modulation and equalization in an orthonormal time-frequency shifting communications system |
| CN107544076A (en) * | 2017-08-10 | 2018-01-05 | 北京航空航天大学 | A kind of navigation neceiver jamproof system suitable for rotating carrier |
-
2018
- 2018-09-18 CN CN201811088759.4A patent/CN109245781A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170201354A1 (en) * | 2010-05-28 | 2017-07-13 | Cohere Technologies, Inc. | Modulation and equalization in an orthonormal time-frequency shifting communications system |
| CN104880714A (en) * | 2015-06-15 | 2015-09-02 | 中国人民解放军国防科学技术大学 | GNSS antenna array anti-interference method based on two-stage filtering structure |
| CN106940689A (en) * | 2017-03-07 | 2017-07-11 | 电子科技大学 | High-precision Eigenvalue Decomposition implementation method based on Jacobi iterative algorithms |
| CN107544076A (en) * | 2017-08-10 | 2018-01-05 | 北京航空航天大学 | A kind of navigation neceiver jamproof system suitable for rotating carrier |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108462521B (en) | Anti-interference realization method of self-adaptive array antenna | |
| Shahbazpanahi et al. | Robust adaptive beamforming for general-rank signal models | |
| CN106972882B (en) | Self-adaptive beam forming method of co-prime array based on virtual domain space power spectrum estimation | |
| CN106483493B (en) | A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method | |
| CN107104720B (en) | Mutual-prime array self-adaptive beam forming method based on covariance matrix virtual domain discretization reconstruction | |
| CN105335615B (en) | A kind of two dimension angular and polarization parameter combined estimation method of low complex degree | |
| CN108594166B (en) | Two-dimensional direction of arrival estimation method and device | |
| CN106872934B (en) | L-type Electromagnetic Vector Sensor Array decorrelation LMS ESPRIT method for parameter estimation | |
| CN106788655B (en) | Interference coherent robust beam forming method for unknown mutual coupling information under mutual coupling condition | |
| CN109782218A (en) | A kind of non-circular signal DOA estimation method of relevant distribution based on double parallel antenna array | |
| CN108872930B (en) | Extended aperture two-dimensional joint diagonalization DOA estimation method | |
| CN109270483B (en) | Three-dimensional array virtual extended coherent source two-dimensional direction of arrival estimation method | |
| CN106291453A (en) | A kind of two-dimentional ripple based on beam space conversion reaches angle assay method | |
| CN112255629A (en) | Sequential ESPRIT two-dimensional incoherent distribution source parameter estimation method based on combined UCA array | |
| CN113835063B (en) | Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method | |
| CN113075610B (en) | DOA estimation method for differential array interpolation based on co-prime polarization array | |
| CN112327292B (en) | DOA estimation method for two-dimensional sparse array | |
| CN114779236A (en) | Improved meter-wave radar low-elevation height measurement method based on spatial smoothing MUSIC | |
| CN112733333A (en) | Two-dimensional direction finding estimation method based on polynomial root finding in co-prime area array | |
| CN112016037A (en) | Two-dimensional direction finding estimation method based on dimensionality reduction Capon root finding in co-prime area array | |
| CN109245781A (en) | Anti-interference method when a kind of annulus conformal array sky | |
| Yang et al. | Robust adaptive beamforming against array calibration errors | |
| CN107247274A (en) | A kind of four quick self-adapted methods of array element anti-interference antenna passage | |
| Chen et al. | A computationally efficient iterative MLE for GPS AOA estimation | |
| CN110208830B (en) | Navigation anti-interference method based on space-time two-dimensional sparse array |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190118 |
|
| RJ01 | Rejection of invention patent application after publication |