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 PDF

Info

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
Application number
CN201810883492.1A
Other languages
Chinese (zh)
Other versions
CN108828509B (en
Inventor
扶彩霞
李万春
马俊虎
陈锐滨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NORTH AUTOMATIC CONTROL TECHNOLOGY INSTITUTE
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810883492.1A priority Critical patent/CN108828509B/en
Publication of CN108828509A publication Critical patent/CN108828509A/en
Application granted granted Critical
Publication of CN108828509B publication Critical patent/CN108828509B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

A kind of multi-platform multi radiation sources bearing relation determination method
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.
CN201810883492.1A 2018-08-06 2018-08-06 Multi-platform multi-radiation source bearing relation judgment method Active CN108828509B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
HUANG WEIPING等: "An intermediary method used in direction-finding location system", 《2010 2ND INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING SYSTEMS》 *
周龙健等: "运动单站对运动目标的快速高精度测向定位研究", 《电子信息对抗技术》 *
黄剑伟等: "一种改进的测向交叉定位方法", 《航天电子对抗》 *

Cited By (4)

* Cited by examiner, † Cited by third party
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