CN108828509A - A kind of multi-platform multi radiation sources bearing relation determination method - Google Patents
A kind of multi-platform multi radiation sources bearing relation determination method Download PDFInfo
- Publication number
- CN108828509A CN108828509A CN201810883492.1A CN201810883492A CN108828509A CN 108828509 A CN108828509 A CN 108828509A CN 201810883492 A CN201810883492 A CN 201810883492A CN 108828509 A CN108828509 A CN 108828509A
- Authority
- CN
- China
- Prior art keywords
- platform
- radiation source
- matrix
- error
- observation station
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to electronics investigative technique fields, are related to a kind of multi-platform multi radiation sources bearing relation determination method.The present invention provides a kind of methods using direction cross positioning, obtain the azimuth information of multiple radiation sources, in conjunction with the position of platform information for having position error, carry out multi-platform and multi radiation sources bearing relation based on the principle of least square and determine.Beneficial effects of the present invention are that the present invention can accurately determine the bearing relation between multi-platform multi radiation sources, and method is simple, works well.
Description
Technical field
The invention belongs to electronics investigative technique fields, are related to a kind of multi-platform multi radiation sources bearing relation determination method.
Background technique
Future battlefield electromagnetic environment will become extremely complex with the development of information technology, and countermeasure search is in complexity
Electromagnetic environment in will become particularly difficult, for example the difficulties such as target location difficulty, the accuracy reliability of surveillance information can be faced
Topic.Realize accurate electromagnetism operation, electronic intelligence reconnaissance is still indispensable important link, in electronic reconnaissance technology,
The position for accurately scouting out platform and electromagnetism threat radiation sources, is the precondition and guarantee for realizing accurate electromagnetism operation, in practical war
In, radiation source is often associated together with platform, is possible to equip a variety of electromagnetism threat radiation sources on a platform, then
Electromagnetism threat radiation sources and platform bearer relationship, which are decrypted, just at essential electronic reconnaissance task is carried out on high-level
Battle field situation and threat estimating key and main foundation, and to the important guarantee that realization of goal is precisely hit.
In practical battlefield, since the attachment relationship of the platforms such as electromagnetism threat radiation sources and large ship is generally not a pair of
One relationship is all in most cases one-to-many situation, i.e. often there are multiple radiation sources on a platform, for such
Situation, it is contemplated that in passive location system, there is position error always for the positioning of platform.
Summary of the invention
The purpose of the present invention is in view of the above-mentioned problems, providing a kind of method using direction cross positioning, and acquisition is multiple
The azimuth information of radiation source is carried out multi-platform in conjunction with the position of platform information for having position error based on the principle of least square
Determine with multi radiation sources bearing relation.
The technical scheme is that:
A kind of space signal power estimation method based on oblique projection, each signal and when the orientation of interference in known spatial
By oblique projection operator realize do not influenced by other interference signal residual powers accurate signal power estimation, following step only with
To angle, θ0For the signal power estimation at place, which is characterized in that include the following steps:
S1, intersect the azimuth that principle obtains multiple radiation sources using direction finding, model is as shown in Figure 1 in direction finding interior extrapolation method.
If observing station coordinates (xm,ym), m=1,2 ... M, M be observation station quantity, a total of K radiation source, by data correlation it
Afterwards, if the observation matrix θ of the bearing data composition of each observation station observationM×KFor:
Wherein, θm,k(m=1,2 ... M, k=1,2 ... K) physical significance indicate the kth that observes of m-th of observation station
The bearing data of a radiation source.For convenience, by θM×KIt is expressed as vector form:θ=[θ1 θ2 ... θk], wherein θk
=[θ1,k θ2,k ... θM,k]T.It is in view of measurement error obtains azimuth determination value:Wherein
n1For observation noise;
S2, known platform position coordinates be P=[xn,yn;xn+1,yn+1;...], n=1,2...N, N are platform
Number adds random error variable in the true coordinate of platform, obtains the platform coordinate matrix containing errorUsing containing wrong
The platform coordinate of differenceBearing data matrix α of the computing platform relative to observation stationM×N, survey station coordinate (xm,ym), m=1,
2 ... M, M are that observation station number calculation formula is as follows:
Wherein n2For random noise, αm,n(m=1,2 ... M, n=1,2 ... N) physical significance indicate n-th of platform
Bearing data relative to m-th of observation station;For convenience, by αM×NIt is expressed as vector form:α=[α1 α2 ... αn],
Wherein, αn=[α1,n α2,n ... αM,n]T。
S3, it is directed to radiation source bearing data matrix θM×KWith platform bearing data matrix αM×N, calculate k-th of radiation source
Relative to the error matrix Δ between the azimuth and platform azimuth of each observation stationk M×N:
Wherein em,n=[| θm,k-αm,1| |θm,k-αm,2| ... |θm,k-αm,n|], θk M×1Indicate that k-th of radiation source is opposite
In the bearing data vector of M radiation source, that is, matrix θM×KKth column.I is the vector that element is all 1;
S4, the corresponding error matrix Δ of k-th of radiation source is calculatedk M×NThe corresponding error sum of squares of each column difference, then
Obtain the error sum of squares matrix E of k-th of radiation sourcek=[E1 E2 ... En], n=1,2...N, if decision threshold is η, if
Have:
En≤ η, then k-th of radiation source belongs to n-th of platform.
Beneficial effects of the present invention are that the present invention can accurately determine the bearing relation between multi-platform multi radiation sources, side
Method is simple, works well.
Detailed description of the invention
Fig. 1 direction finding intersects schematic diagram
Fig. 2 emulates schematic diagram
Fig. 3 determines accuracy with decision threshold change curve
Fig. 4 determines accuracy with angle error change curve
Fig. 5 determines accuracy with platform spacing change curve
Fig. 6 determines accuracy with platform position error change curve
Specific embodiment
Below in conjunction with attached drawing and emulation, technical solution of the present invention is further described.
This example will carry out simulating, verifying to mentioned method using matlab, for simplicity, make to algorithm model following false
If:
1. all observation stations and target are all in X/Y plane;
2. all observation station direction finding precisions having the same;
3. all Engineering Errors are all added in angle measurement error;
If in the square region that observation station, platform and radiation source are 50km × 50km in region, it is known that 3 platforms of distribution,
Their coordinate is [10.6,35.5;29.2 30.8;46.4,35.3], unit km.The radiation source number being distributed on three platforms
Amount is [2,2,1], is observed to obtain the azimuth information of radiation source, 3 sights to above-mentioned radiation source using 3 fixed observer stations
The coordinate of survey station be respectively (0.5,0.5), (25,5), (50,1), unit km, emulation schematic diagram it is as shown in Figure 2.Observation station
Angle error obeys the Gaussian Profile that mean value is zero, and mutually indepedent between the angle error of 3 observation stations.
Simulation result is as shown in the figure:
Simulation scenarios 1:
0.5 degree of angle error, platform spacing 5km, platform position error are 0.5km, threshold variation range is:1~5 is (single
Position:km);
Decision threshold is bigger it can be seen from Fig. 3 result, and platform determines that accuracy is higher with radiation source bearing relation, when
Decision threshold determines that accuracy growth rate is very fast, after thresholding is greater than 4km, the variation of accuracy becomes in 1km~4km
Gesture is slower, but accuracy reaches 99%, illustrate when the setting of accuracy decision threshold in 4km or more, using this method into
Row platform determines that effect is rather good with radiation source bearing relation.
Simulation scenarios 2:
Angle error variation range:0.5~5 (degree), platform spacing 5km, platform position error are 0.5km, thresholding is
4km;
Angle error is bigger it can be seen from Fig. 4 result, and platform determines that accuracy is lower with radiation source bearing relation, when
Angle error determines accuracy 80% or more, after thresholding is greater than 2 degree, accuracy is gradually decreased at 1~2 degree
50% or so, illustrate after observation station angle error is greater than 2 degree, carries out platform using this method and sentence with radiation source bearing relation
It is very poor to determine effect.
Simulation scenarios 3:
0.5 degree of angle error, platform spacing variation range are:1~5km), platform position error be 0.5km, thresholding is
4km;
Platform spacing is bigger it can be seen from Fig. 5 result, and platform determines that accuracy is higher with radiation source bearing relation, when
Decision threshold determines that accuracy growth rate is very fast, after thresholding is greater than 3km, the variation of accuracy becomes in 1km~3km
Gesture relatively slowly even almost no longer changes, and accuracy only maintains between 80%~90% as seen from the figure, it can be seen that
When within 3km, capable of seeing apparent effect using this method, but when spacing is in 3km or more, utilize this method between platform
It carries out platform and determines that accuracy variation effect is unobvious with radiation source bearing relation.
Simulation scenarios 4:
0.5 degree of angle error, platform spacing 5km, platform position error variation range are:0.1~1km, thresholding 4km;
Platform positioning is bigger it can be seen from Fig. 6 result, and platform determines that accuracy is lower with radiation source bearing relation, but
Accuracy maintains 90% or more always, and variation effect is unobvious, it can be seen that when within 1km, utilizing the party between platform
Method carries out platform and determines that accuracy variation effect is all good with radiation source bearing relation.
Claims (1)
1. a kind of multi-platform multi radiation sources bearing relation determination method, which is characterized in that include the following steps:
S1, intersect the azimuth that former method obtains multiple radiation sources using direction finding:If observation station coordinates is (xm,ym), m=1,2 ... M,
M is observation station quantity, a total of K radiation source, after data correlation, what the bearing data of each observation station observation formed
Observation matrix θM×KFor:
Wherein, θm,kPhysical significance indicate the bearing data of k-th of radiation source that m-th of observation station observes, m=1,
2 ... M, k=1,2 ... K;θM×KBe at vector form:θ=[θ1 θ2...θk], wherein θk=[θ1,k θ2,k...θM,k
]T;It is in view of measurement error obtains azimuth determination value:Wherein n1For observation noise;
S2, known platform position coordinates be P=[xn,yn;xn+1,yn+1;...], n=1,2...N, N are platform number,
Random error variable is added in the true coordinate of platform, obtains the platform coordinate matrix containing errorUsing containing error
Platform coordinateBearing data matrix α of the computing platform relative to observation stationM×N:
Wherein n2For random noise, αm,nPhysical significance indicate bearing data of n-th of platform relative to m-th of observation station, m=
1,2 ... M, n=1,2 ... N;αM×NVector form:α=[α1 α2...αn], wherein αn=[α1,n α2,n...αM,n]T;
S3, it is directed to radiation source bearing data matrix θM×KWith platform bearing data matrix αM×N, calculate k-th of radiation source relative to
Error matrix Δ between the azimuth and platform azimuth of each observation stationk M×N:
Wherein em,n=[| θm,k-αm,1| |θm,k-αm,2|...|θm,k-αm,n|], θk M×1Indicate k-th of radiation source relative to M spoke
Penetrate the bearing data vector in source, that is, matrix θM×KKth column, I is the vector that element is all 1;
S4, the corresponding error matrix Δ of k-th of radiation source is calculatedk M×NThe corresponding error sum of squares of each column difference, obtains k-th
The error sum of squares matrix E of radiation sourcek=[E1 E2...En], n=1,2...N, if decision threshold is η, if having:
En≤ η, then k-th of radiation source belongs to n-th of platform.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883492.1A CN108828509B (en) | 2018-08-06 | 2018-08-06 | Multi-platform multi-radiation source bearing relation judgment method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883492.1A CN108828509B (en) | 2018-08-06 | 2018-08-06 | Multi-platform multi-radiation source bearing relation judgment method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108828509A true CN108828509A (en) | 2018-11-16 |
CN108828509B CN108828509B (en) | 2022-01-11 |
Family
ID=64153630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810883492.1A Active CN108828509B (en) | 2018-08-06 | 2018-08-06 | Multi-platform multi-radiation source bearing relation judgment method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108828509B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110673090A (en) * | 2019-10-14 | 2020-01-10 | 电子科技大学 | Passive multi-station multi-target positioning method based on DBSCAN |
CN111999696A (en) * | 2020-08-02 | 2020-11-27 | 中国人民解放军海军工程大学 | Multi-platform direction-finding cross-positioning optimization method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105601A (en) * | 2013-01-07 | 2013-05-15 | 哈尔滨工业大学 | Maximum posterior principle radiation source position method based on grid search |
US20140079248A1 (en) * | 2012-05-04 | 2014-03-20 | Kaonyx Labs LLC | Systems and Methods for Source Signal Separation |
CN106654564A (en) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学(威海) | Phase interferometer based on broadband conformal antenna array and parameter estimation method thereof |
CN107861096A (en) * | 2017-11-03 | 2018-03-30 | 中国人民解放军陆军炮兵防空兵学院 | Least square direction-finding method based on voice signal reaching time-difference |
-
2018
- 2018-08-06 CN CN201810883492.1A patent/CN108828509B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140079248A1 (en) * | 2012-05-04 | 2014-03-20 | Kaonyx Labs LLC | Systems and Methods for Source Signal Separation |
CN103105601A (en) * | 2013-01-07 | 2013-05-15 | 哈尔滨工业大学 | Maximum posterior principle radiation source position method based on grid search |
CN106654564A (en) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学(威海) | Phase interferometer based on broadband conformal antenna array and parameter estimation method thereof |
CN107861096A (en) * | 2017-11-03 | 2018-03-30 | 中国人民解放军陆军炮兵防空兵学院 | Least square direction-finding method based on voice signal reaching time-difference |
Non-Patent Citations (3)
Title |
---|
HUANG WEIPING等: "An intermediary method used in direction-finding location system", 《2010 2ND INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING SYSTEMS》 * |
周龙健等: "运动单站对运动目标的快速高精度测向定位研究", 《电子信息对抗技术》 * |
黄剑伟等: "一种改进的测向交叉定位方法", 《航天电子对抗》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110673090A (en) * | 2019-10-14 | 2020-01-10 | 电子科技大学 | Passive multi-station multi-target positioning method based on DBSCAN |
CN110673090B (en) * | 2019-10-14 | 2022-08-05 | 电子科技大学 | Passive multi-station multi-target positioning method based on DBSCAN |
CN111999696A (en) * | 2020-08-02 | 2020-11-27 | 中国人民解放军海军工程大学 | Multi-platform direction-finding cross-positioning optimization method |
CN111999696B (en) * | 2020-08-02 | 2023-07-04 | 中国人民解放军海军工程大学 | Multi-platform direction-finding cross positioning optimization method |
Also Published As
Publication number | Publication date |
---|---|
CN108828509B (en) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106407677B (en) | A kind of multi-object tracking method in the case of missing measurement | |
CN109212475B (en) | Dual-machine passive positioning method based on azimuth angle and pitch angle information | |
CN103941233B (en) | The radar interval alternately radiation control method of tracking worked in coordination with by a kind of multi-platform main passive sensor | |
CN108061877A (en) | A kind of passive track-corelation direction cross positioning method based on angle information | |
CN108089148A (en) | A kind of passive track-corelation direction cross positioning method based on time difference information | |
CN107015199A (en) | A kind of double unmanned plane direction finding time difference positioning methods for considering UAV Attitude angle | |
CN109145474A (en) | A kind of infrared imaging device recon effect analogy method suitable for emulation of fighting | |
CN108828509A (en) | A kind of multi-platform multi radiation sources bearing relation determination method | |
CN105676170B (en) | A kind of False Intersection Points minimizing technology positioned for Multi-Target Passive with tracking | |
CN108562945A (en) | A kind of localization method and its device of more magnetic movement objective bodies | |
CN108152789A (en) | Utilize the passive track-corelation data correlation and localization method of RSS information | |
CN103926569B (en) | Three-dimensional radar net is based on the associated centralization of cross bearing point and interferes discrimination method with distributed compacting | |
CN106646395A (en) | Radar echo deduction method for flight target | |
CN105425206A (en) | Steady least square positioning method in nonsynchronous wireless network | |
CN106643297A (en) | Estimation and correction method for vector miss distance parameters of motion platform | |
CN103499809B (en) | A kind of Pure orientation double computer cooperation target following location path planing method | |
CN101308206A (en) | Circumferential track mobile target tracking method under white noise background | |
CN104977559B (en) | Target positioning method in interference environment | |
US10295662B2 (en) | Producing data describing target measurements | |
Hanson | Using multiple MEMS IMUs to form a distributed inertial measurement unit | |
CN109188417A (en) | The method that single station Passive Positioning is carried out to scan-type radiation source using unmanned aerial vehicle platform | |
CN104391272A (en) | Method and system for target positioning by using direction finding data | |
CN109884582A (en) | The method of target three-dimensional coordinate is quickly determined using one-dimensional direction finding | |
CN111999696B (en) | Multi-platform direction-finding cross positioning optimization method | |
KR101310199B1 (en) | Passive ranging module |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240313 Address after: No. 351 Tiyu Road, Xiaodian District, Taiyuan City, Shanxi Province 030000 Patentee after: NORTH AUTOMATIC CONTROL TECHNOLOGY INSTITUTE Country or region after: China Address before: 611731, No. 2006, West Avenue, hi tech West District, Sichuan, Chengdu Patentee before: University of Electronic Science and Technology of China Country or region before: China |