CN108828509B - Multi-platform multi-radiation source bearing relation judgment method - Google Patents

Abstract

The invention belongs to the technical field of electronic investigation, and relates to a method for judging bearing relationship of multiple platforms and multiple radiation sources. The invention provides a method for cross positioning by using direction finding, which is used for obtaining azimuth angle information of a plurality of radiation sources, and carrying relation judgment of a plurality of platforms and a plurality of radiation sources is carried out based on the least square principle by combining platform position information with positioning errors. The method has the advantages that the method can accurately judge the bearing relationship among the multiple platforms and the multiple radiation sources, and is simple and good in effect.

Description

Multi-platform multi-radiation source bearing relation judgment method

Technical Field

The invention belongs to the technical field of electronic investigation, and relates to a method for judging bearing relationship of multiple platforms and multiple radiation sources.

Background

The electromagnetic environment of the battlefield in the future becomes extremely complex along with the development of information technology, and electronic countermeasure reconnaissance becomes particularly difficult in the complex electromagnetic environment, such as meeting the difficult problems of target positioning, accuracy and reliability of reconnaissance information and the like. In the electronic reconnaissance technology, the positions of a platform and an electromagnetic threat radiation source are accurately reconnaissance, the electronic information reconnaissance is a prerequisite and guarantee for realizing the accurate electromagnetic combat, in an actual battlefield, the radiation source is often associated with the platform, multiple electromagnetic threat radiation sources are possibly equipped on one platform, so that the electromagnetic threat radiation source and the platform are decrypted to bear the relationship, an indispensable electronic reconnaissance task is formed, the key and main basis for performing situation estimation and threat estimation on a high level is provided, and the important guarantee for realizing the accurate strike on a target is provided.

In an actual battlefield, because the attachment relationship between the electromagnetic threat radiation source and the platforms such as large ships is not one-to-one generally, most of the cases are one-to-many cases, that is, a plurality of radiation sources often exist on one platform, and for the cases, the positioning error always exists in the positioning of the platform in consideration of a passive positioning system.

Disclosure of Invention

The invention aims to solve the problems and provides a method for obtaining azimuth angle information of a plurality of radiation sources by using direction-finding cross positioning, and carrying relation judgment of a plurality of platforms and a plurality of radiation sources is carried out based on the least square principle by combining platform position information with positioning errors.

The technical scheme of the invention is as follows:

a method for estimating space signal power based on oblique projection realizes accurate signal power estimation without being influenced by residual power of other interference signals through an oblique projection operator when the directions of various signals and interference in a known space are known, and the following steps only adopt the diagonal theta₀The method for estimating signal power is characterized by comprising the following steps of:

and S1, obtaining azimuth angles of the plurality of radiation sources by using the principle of direction-finding intersection, wherein the model in the direction-finding intersection method is shown in figure 1. Set observation station coordinates (x)_m,y_m) M is 1,2, M is the number of observation stations, there are K radiation sources in total, and after data association, an observation value matrix theta formed by the azimuth data observed by each observation station is set^M×KComprises the following steps:

wherein, theta_m,kThe physical meaning of (M1, 2,. M, K1, 2,. K) represents the azimuth data of the kth radiation source observed by the mth observation station. For convenience, let θ^M×KExpressed in vector form: theta is ═ theta₁ θ₂ ... θ_k]Wherein, theta_k＝[θ_1,k θ_2,k ... θ_M,k]^T. The azimuthal measurements taken to account for measurement errors are:

wherein n is₁To observe noise;

s2, knowing the position coordinate of the platform as P＝[x_n,y_n；x_n+1,y_n+1；......]N is the number of platforms, and random error variables are added to real coordinates of the platforms to obtain a platform coordinate matrix containing errors

Using platform coordinates containing errors

Computing platform azimuth data matrix alpha relative to observation station^M×NCoordinate of the measuring station (x)_m,y_m) M is 1,2.. M, and M is the calculation formula of the number of observation stations as follows:

wherein n is₂Being random noise, α_m,nThe physical meaning of (M1, 2,. M, N1, 2,. N) represents the orientation data of the nth platform relative to the mth observation station; for convenience, will be alpha^M×NExpressed in vector form: alpha ═ alpha₁ α₂ ... α_n]Wherein α is_n＝[α_1,n α_2,n ... α_M,n]^T。

S3 matrix theta for radiation source orientation data^M×KAnd platform azimuth data matrix alpha^M×NCalculating an error matrix delta between the azimuth angle of the kth radiation source relative to each observation station and the platform azimuth angle_k ^M×N：

Wherein e_m,n＝[|θ_m,k-α_m,1| |θ_m,k-α_m,2| ... |θ_m,k-α_m,n|]，θ_k ^M×1Representing the orientation data vector of the k-th radiation source relative to the M radiation sources, i.e. the matrix theta^M×KThe k-th column of (1). I is a vector with elements all being 1;

s4, calculating an error matrix delta corresponding to the k radiation source_k ^M×NThe error sum of squares corresponding to each column difference is obtained, and then the error sum of squares matrix E of the k radiation source is obtained_k＝[E₁ E₂ ... E_n]N is 1,2.. N, the judgment threshold is set as η, if:

E_nwhen eta is less than or equal to eta, the kth radiation source belongs to the nth platform.

The method has the advantages that the method can accurately judge the bearing relationship among the multiple platforms and the multiple radiation sources, and is simple and good in effect.

Drawings

FIG. 1 is a schematic view of a direction-finding intersection

FIG. 2 is a simulation diagram

FIG. 3 is a graph showing the variation of the accuracy of the decision with the decision threshold

FIG. 4 is a graph showing the variation of the determination accuracy with the angle measurement error

FIG. 5 is a graph showing the variation of the determination accuracy with the stage pitch

FIG. 6 is a graph showing the variation of the accuracy with the positioning error of the platform

Detailed Description

The technical solution of the present invention will be further explained with reference to the drawings and simulations.

In this example, matlab is used to perform simulation verification on the proposed method, and for simplicity, the following assumptions are made for the algorithm model:

1. all stations and targets are in the XY plane;

2. all observation stations have the same direction-finding precision;

3. all engineering errors are superposed into the direction-finding errors;

setting an observation station, a platform and a radiation source in a square area with the area of 50km multiplied by 50km, wherein 3 platforms are known to be distributed, and the coordinates of the platforms are [10.6, 35.5; 29.2, 30.8; 46.4, 35.3], in km. The number of radiation sources distributed on the three platforms is [2,2,1], 3 fixed observation stations are used for observing the radiation sources to obtain azimuth angle information of the radiation sources, coordinates of the 3 observation stations are (0.5 ), (25,5), (50,1) respectively, the unit is km, and a simulation schematic diagram is shown in fig. 2. The angle measurement errors of the observation stations are subjected to Gaussian distribution with the mean value of zero, and the angle measurement errors of the 3 observation stations are independent.

The simulation results are shown in the figure:

simulation case 1:

the angle measurement error is 0.5 degree, the platform distance is 5km, the platform positioning error is 0.5km, and the threshold variation range is as follows: 1 to 5 (unit: km);

it can be seen from the results of fig. 3 that the larger the determination threshold is, the higher the determination accuracy of the bearing relationship between the platform and the radiation source is, when the determination threshold is between 1km and 4km, the increase speed of the determination accuracy is fast, and when the threshold is greater than 4km, the change trend of the accuracy is slow, but the accuracy reaches 99%, which indicates that when the determination threshold of the accuracy is set to be greater than 4km, the determination effect of the bearing relationship between the platform and the radiation source by using the method is very good.

Simulation case 2:

the variation range of the angle measurement error is as follows: 0.5-5 degrees, 5km of platform distance, 0.5km of platform positioning error and 4km of threshold;

the result of fig. 4 shows that the larger the angle measurement error is, the lower the accuracy of judging the bearing relationship between the platform and the radiation source is, when the angle measurement error is 1-2 degrees, the accuracy is more than 80%, and when the threshold is more than 2 degrees, the accuracy gradually decreases to about 50%, which indicates that when the angle measurement error of the observation station is more than 2 degrees, the effect of judging the bearing relationship between the platform and the radiation source by using the method is extremely poor.

Simulation case 3:

the angle measurement error is 0.5 degrees, and the change range of the platform distance is as follows: 1-5 km), the platform positioning error is 0.5km, and the threshold is 4 km;

it can be seen from the results of fig. 5 that the larger the distance between the platforms, the higher the accuracy of the judgment of the bearing relationship between the platform and the radiation source, when the judgment threshold is 1km to 3km, the increase speed of the accuracy of the judgment is faster, when the threshold is greater than 3km, the change trend of the accuracy is slower or even almost no longer changes, and it can be seen from the figure that the accuracy is only maintained between 80% and 90%, so that when the distance between the platforms is within 3km, the obvious effect can be seen by using the method, but when the distance is more than 3km, the effect of the accuracy of the judgment of the bearing relationship between the platform and the radiation source by using the method is not obvious.

Simulation case 4:

the angle measurement error is 0.5 degree, the platform distance is 5km, and the variation range of the platform positioning error is as follows: 0.1-1 km, and a threshold of 4 km;

it can be seen from the results of fig. 6 that the larger the platform is positioned, the lower the accuracy of the judgment of the bearing relationship between the platform and the radiation source is, but the accuracy is always maintained at more than 90%, and the change effect is not obvious, so that when the platform is within 1km, the change effect of the accuracy of the judgment of the bearing relationship between the platform and the radiation source by using the method is good.

Claims (1)

1. A method for determining the bearing relationship of a multi-platform multi-radiation source is characterized by comprising the following steps:

s1, azimuth angles of a plurality of radiation sources by using a direction-finding cross original method are adopted: let the observation station coordinate be (x)_m,y_m) M is 1,2, M is the number of observation stations, there are K radiation sources in total, and after data association, the observation value matrix theta is formed by the azimuth data observed by each observation station^M×KComprises the following steps:

wherein, theta_m,kThe physical meaning of (a) denotes the azimuth data of the kth radiation source observed by the mth observation station, M being 1,2,... M, K being 1,2,... K; theta^M×KThe vector-forming form of (a): theta is ═ theta₁ θ₂...θ_k]Wherein, theta_k＝[θ_1,k θ_2,k...θ_M,k]^T(ii) a The azimuthal measurements taken to account for measurement errors are:

wherein n is₁To observe noise;

s2, knowing the position coordinate of the platform as P ═ x_n,y_n；x_n+1,y_n+1；......]N is the number of platforms, and random error variables are added to real coordinates of the platforms to obtain a platform coordinate matrix containing errors

Using platform coordinates containing errors

Computing platform azimuth data matrix alpha relative to observation station^M×N：

Wherein n is₂Being random noise, α_m,nThe physical meaning of (a) denotes the position data of the nth platform relative to the mth observation station, M being 1,2,. M, N being 1,2,. N; alpha is alpha^M×NVector form of (2): alpha ═ alpha₁ α₂...α_n]Wherein α is_n＝[α_1,n α_2,n...α_M,n]^T；

S3 matrix theta for radiation source orientation data^M×KAnd platform azimuth data matrix alpha^M×NCalculating an error matrix delta between the azimuth angle of the kth radiation source relative to each observation station and the platform azimuth angle_k ^M×N：

Wherein e_m,n＝[|θ_m,k-α_m,1| |θ_m,k-α_m,2|...|θ_m,k-α_m,n|]，θ_k ^M×1Representing the orientation data vector of the k-th radiation source relative to the M radiation sources, i.e. the matrix theta^M×KI is a vector with all elements 1;

s4, calculating an error matrix delta corresponding to the k radiation source_k ^M×NThe error square sum corresponding to each column difference value is obtained to obtain an error square sum matrix E of the kth radiation source_k＝[E₁ E₂...E_n]N is 1,2.. N, the judgment threshold is set as η, if:

E_nwhen eta is less than or equal to eta, the kth radiation source belongs to the nth platform.

Images