CN108008351B - Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation - Google Patents

Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation Download PDF

Info

Publication number
CN108008351B
CN108008351B CN201711222541.9A CN201711222541A CN108008351B CN 108008351 B CN108008351 B CN 108008351B CN 201711222541 A CN201711222541 A CN 201711222541A CN 108008351 B CN108008351 B CN 108008351B
Authority
CN
China
Prior art keywords
signal
amplitude
emission
coordinate
equal
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.)
Active
Application number
CN201711222541.9A
Other languages
Chinese (zh)
Other versions
CN108008351A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201711222541.9A priority Critical patent/CN108008351B/en
Publication of CN108008351A publication Critical patent/CN108008351A/en
Application granted granted Critical
Publication of CN108008351B publication Critical patent/CN108008351B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0278Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves involving statistical or probabilistic considerations

Landscapes

  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

A signal positioning method based on forward attenuation reverse-thrust emission amplitude aggregation comprises the following steps: 1. setting a checking area, and performing mesh generation on a map of the checking area; 2. mounting signal monitoring stations in any A grids in a map grid; 3. calculating or measuring any coordinate position S in the gridi,jTo the coordinate position Sm,nA signal space attenuation value of; 4. suppose the region is divided by Sm,nAll other grid positions may have suspected signal sources, and Si,jAt a signal source emission amplitude of
Figure DDA0001486659420000011
5. If there are A transmit amplitudes at the same coordinate location, there are b transmit amplitude values at
Figure DDA0001486659420000012
Within the confidence interval of (2), the coordinate is considered to be f0Suspected emission coordinates of the signal; 6. sorting all suspected emission coordinates from large to small according to the b value; the larger b, the greater the likelihood of indicating that the location is a true transmit coordinate.

Description

Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation
Technical Field
The invention relates to a signal positioning method based on forward attenuation backward-pushing emission amplitude aggregation, and belongs to the field of signal monitoring.
Background
As the number of radio transmitting devices used and the density of spectral coverage increases, the electrical spectrum resources become increasingly strained. In the face of possible economic loss, even national security, it is important to accurately investigate radio signal sources or illegal stations and ensure smooth communication.
The modern radio direction finding can be roughly divided into the following four systems according to different direction finding principles: amplitude comparison method direction finding, space spectrum estimation method direction finding, arrival time difference method direction finding and an interferometer. (1) Amplitude comparison and direction finding system: the incoming wave direction is measured by using the directional characteristics of the direction-finding antenna array to receive the signal amplitude of the incoming wave in different directions. Taking the widest Watson-Watt system as an example, the system has the advantages of insensitivity to channel interference, high direction finding speed and easiness in realization, but the system has low direction finding precision and direction finding sensitivity and weak capability of resisting wave front distortion. (2) The direction finding principle of the space spectrum estimation direction finding system is as follows: the received signal of each antenna element is compared with the signals of other antenna elements, and taking a typical MUSIC algorithm as an example, the method is used for estimating the Direction of arrival (DOA), and has the advantages of high Direction finding precision, small Direction finding error, high resolution, capability of simultaneously orienting multiple signals and the like, but still has some defects, for example, when multiple incoming waves are mutually coherent, the MUSIC algorithm fails. (3) Arrival time difference direction finding: three or more than three measuring stations are used for carrying out data acquisition on the arrival time of a certain signal and carrying out positioning plane or space positioning on a target source. The direction finding accuracy is high, the direction finding speed is high, the sensitivity is high, and the method is insensitive to the measuring environment. However, the arrival time direction finding requires that the incoming wave signal is a fast time-varying signal and has a fixed modulation mode, and under the condition that the digital modulation is widely applied at present, the requirement of the incoming wave signal on the modulation mode is difficult to meet, so the application is less. (4) An interferometer: the interferometer is divided into a phase interferometer and a correlation interferometer, the phase interferometer is not provided with a pre-established sample library, and the DOA phase difference is solved; the correlation interferometer needs to establish a sample library by the azimuth angle and the pitch angle according to a certain resolution in advance, and phase difference vectors of incoming waves are compared with the sample library one by one.
With the deterioration of the electromagnetic environment, the positioning method using the above principle causes the accuracy to be reduced, the error to be increased and the positioning algorithm to be more and more complex.
Disclosure of Invention
The invention provides a signal positioning method based on forward attenuation reverse-push emission amplitude aggregation, which is used for checking a spatial unknown signal source (emission frequency f)0) The specific location of (a). For convenience of calculation, the amplitude unit of the signals is unified to dBm; to simplify the model, the signal source is set to omni-directional transmission. The method specifically comprises the following steps:
the method comprises the following steps: setting a checking area, meshing a map of the checking area, setting the meshes to be I multiplied by J, and labeling a position coordinate S for each meshi,j(I is less than or equal to I, J is less than or equal to J), forming a coordinate matrix [ S ] of the area]I,J
Step two: and (4) installing signal monitoring stations at any A grids (namely position coordinates) (A is more than or equal to 3) in the map grids obtained in the step one. The kth monitoring station is denoted as Rk(kA) or less) and the position coordinate is Sm,n. At the monitoring station RkMeasuring the amplitude of the signal to be
Figure BDA0001486659400000021
Step three: calculating or measuring any coordinate position S in the grid by using radio wave propagation methodi,jTo the coordinate position Sm,nA signal space attenuation value of
Figure BDA0001486659400000022
Forming all the positions of the region to Sm,nMatrix of attenuation values [ DRk]I,J
Step four: suppose the region is divided by Sm,nAll other grid positions may have suspected signal sources, and Si,jAt a signal source emission amplitude of
Figure BDA0001486659400000023
Is obtained as a coordinate position Sm,nEmitting amplitude matrix at each position of region as reference point
Figure BDA0001486659400000024
Step five: and repeating the second step to the fourth step to obtain A transmitting amplitude matrixes of all A monitoring stations. Thus, for the same coordinate Si,jA number of transmit amplitude values may be obtained. Average amplitude of emission
Figure BDA0001486659400000025
Setting an amplitude confidence coefficient delta, if A emission amplitudes at the same coordinate position exist b (b is more than or equal to 2 and less than or equal to A) emission amplitude values
Figure BDA0001486659400000026
Within the confidence interval of (2), the coordinate position is considered as the frequency f0The suspected emission coordinates of the signal.
Step six: and D, sorting all suspected emission coordinates in the step five from large to small according to the b value. The larger the value of b, the greater the likelihood of indicating that the location is a true transmit coordinate.
The invention has the advantages that: the method avoids a complex passive positioning algorithm, and can realize the positioning of the suspected signal source by actually measuring the signal amplitude data through the monitoring station and combining with an attenuation value matrix obtained by a radio wave propagation algorithm.
Drawings
Fig. 1 is a schematic diagram of the mesh generation and monitoring station of the area to be investigated and the setting position of a suspected signal source.
Detailed Description
The invention provides a signal positioning method based on forward attenuation reverse-thrust emission amplitude aggregation. The present invention will be described in further detail with reference to the accompanying drawings and examples.
Step one, as shown in fig. 1, mesh division is performed on a set investigation region map, the mesh is set to 11 × 11, and a position coordinate S is marked on each meshi,j(i is less than or equal to 11, j is less than or equal to 11), forming a coordinate matrix [ S ] of the area]11,11
Step two, as shown in fig. 1, 3 grids, i.e. position coordinates S, are selected as an example8,1,S1,7,S5,11Department installation signal monitoring station R1、R2、R3And respectively measuring the signal amplitude of
Figure BDA0001486659400000031
Step three, calculating by a parabolic equation method to obtain each coordinate position Si,jTo the coordinate position S8,1,S1,7,S5,11Forming all positions of the region to S respectively8,1,S1,7,S5,11Matrix of attenuation values [ DR1]11,11,[DR2]11,11,[DR3]11,11
Step four, for each monitoring station, assuming that all points of the area except the position of the monitoring station are possible to have suspected signal sources, each grid Si,jThe transmission amplitude obtained by the signal source respectively for the measured values of the monitoring station is
Figure BDA0001486659400000032
Figure BDA0001486659400000033
Thereby respectively obtaining the coordinate position S8,1,S1,7,S5,11Emitting amplitude matrix at each position of region as reference point
Figure BDA0001486659400000034
Step five, the obtained three emission amplitude matrixes
Figure BDA0001486659400000035
Same coordinate Si,jThere are 3 values of the emission amplitude, thus obtaining Si,jAverage emission amplitude of
Figure BDA0001486659400000036
Setting an amplitude confidence delta if there are 3 transmit amplitudes at the same coordinate position, there are b transmit amplitude values
Figure BDA0001486659400000037
Within the confidence interval of (c). Whereby the position coordinate is considered as a frequency f0The suspected emission coordinates of the signal.
And step six, sorting all the suspected emission coordinates according to the b value from large to small, wherein b is more than or equal to 2 and less than or equal to A, and in the example, A is 3, so that b is 2 or b is 3. The larger the value of b, the more likely it is that the position is the true transmit coordinate, i.e.: has 3 transmission amplitude values
Figure BDA0001486659400000038
Is compared with 2 transmission amplitude values
Figure BDA0001486659400000039
Is more likely to be the signal emission location.
The above embodiments are only used to illustrate the implementation process of the present invention, and the protection scope of the present invention is not limited to the above embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention.

Claims (1)

1. A signal positioning method based on forward attenuation reverse-push emission amplitude aggregation is characterized in that: for investigating spatially unknown signal sources, i.e. transmission frequencies f0The specific location of the signal of (a); for convenience of calculation, the amplitude unit of the signals is unified to dBm; in order to simplify the model, a signal source is set to be omnidirectional emission; the method specifically comprises the following steps:
the method comprises the following steps: setting a checking area, meshing a map of the checking area, setting the meshes to be I multiplied by J, and labeling a position coordinate S for each meshi,jForm the coordinate matrix [ S ] of the area]I,J(ii) a Wherein I is less than or equal to I, and J is less than or equal to J;
step two: mounting a signal monitoring station at any A grids, namely position coordinates, in the map grids obtained in the step one; the kth monitoring station is denoted as RkWith a position coordinate Sm,n(ii) a At the monitoring station RkMeasuring the amplitude of the signal to be
Figure FDA0002946734990000011
Wherein A is more than or equal to 3; k is less than or equal to A;
step three: calculating or actually measuring any coordinate position S in grid by radio wave propagation methodi,jTo the coordinate position Sm,nA spatial attenuation value of the signal of
Figure FDA0002946734990000012
Forming all the positions of the region to Sm,nMatrix of attenuation values [ DRk]I,J
Step four: suppose the region is divided by Sm,nAll other grid positions may have suspected signal sources, and Si,jAt a signal source emission amplitude of
Figure FDA0002946734990000013
Is obtained as a coordinate position Sm,nEmitting amplitude matrix at each position of region as reference point
Figure FDA0002946734990000014
Step five: repeating the second step to the fourth step to obtain A emission amplitude matrixes of all A monitoring stations; thus, for the same coordinate Si,jObtaining A transmitting amplitude values; average amplitude of emission
Figure FDA0002946734990000015
Setting an amplitude confidence degree delta, if A emission amplitudes at the same coordinate position exist b emission amplitude values
Figure FDA0002946734990000016
Within the confidence interval of (2), the coordinate position is considered as the frequency f0Suspected emission coordinates of the signal; wherein b is more than or equal to 2 and less than or equal to A;
step six: sorting all suspected emission coordinates in the fifth step from large to small according to the b value; the larger the value of b, the greater the likelihood of indicating that the location is a true transmit coordinate.
CN201711222541.9A 2017-11-29 2017-11-29 Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation Active CN108008351B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711222541.9A CN108008351B (en) 2017-11-29 2017-11-29 Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711222541.9A CN108008351B (en) 2017-11-29 2017-11-29 Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation

Publications (2)

Publication Number Publication Date
CN108008351A CN108008351A (en) 2018-05-08
CN108008351B true CN108008351B (en) 2021-05-14

Family

ID=62054446

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711222541.9A Active CN108008351B (en) 2017-11-29 2017-11-29 Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation

Country Status (1)

Country Link
CN (1) CN108008351B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116456376B (en) * 2023-03-10 2024-01-23 深圳市华谱通信技术有限公司 Method and system for controlling signal measuring station in radio detection network

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101526609A (en) * 2009-03-27 2009-09-09 电子科技大学 Matching locating method based on wireless channel frequency domain amplitude response
CN101547130A (en) * 2008-03-26 2009-09-30 中国人民解放军信息工程大学 Positioning method and positioning device based on wireless sensor network
FR2957185A1 (en) * 2010-03-04 2011-09-09 Canon Kk Threshold determination method for attenuating noise generated in e.g. room, involves determining threshold to be applied to component of sound signal based on noise level and amplitude, and maintaining noise generated by sound signal
CN104010366A (en) * 2014-05-12 2014-08-27 电子科技大学 Environment self-adaptation signal source positioning method
CN104284418A (en) * 2014-10-08 2015-01-14 海南大学 Method for positioning signals of mobile node user in commercial center
CN106714227A (en) * 2016-08-31 2017-05-24 成都中星世通电子科技有限公司 Automatic positioning method for mobile monitoring station
EP3201645A1 (en) * 2014-09-29 2017-08-09 Telecom Italia S.p.A. Positioning method and system for wireless communication networks

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8520647B2 (en) * 2007-09-26 2013-08-27 Aruba Networks, Inc. Wireless client position estimating system and method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101547130A (en) * 2008-03-26 2009-09-30 中国人民解放军信息工程大学 Positioning method and positioning device based on wireless sensor network
CN101526609A (en) * 2009-03-27 2009-09-09 电子科技大学 Matching locating method based on wireless channel frequency domain amplitude response
FR2957185A1 (en) * 2010-03-04 2011-09-09 Canon Kk Threshold determination method for attenuating noise generated in e.g. room, involves determining threshold to be applied to component of sound signal based on noise level and amplitude, and maintaining noise generated by sound signal
CN104010366A (en) * 2014-05-12 2014-08-27 电子科技大学 Environment self-adaptation signal source positioning method
EP3201645A1 (en) * 2014-09-29 2017-08-09 Telecom Italia S.p.A. Positioning method and system for wireless communication networks
CN104284418A (en) * 2014-10-08 2015-01-14 海南大学 Method for positioning signals of mobile node user in commercial center
CN106714227A (en) * 2016-08-31 2017-05-24 成都中星世通电子科技有限公司 Automatic positioning method for mobile monitoring station

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
关于大规模室内目标信号定位自动校准研究;张志涛;《计算机仿真》;20160131;第33卷(第1期);第450-455页 *

Also Published As

Publication number Publication date
CN108008351A (en) 2018-05-08

Similar Documents

Publication Publication Date Title
AU2022259835B2 (en) Direction of arrival estimation
Zhang et al. Distributed angle estimation for localization in wireless sensor networks
US10145933B2 (en) Angle determining system and method
Wen et al. 5G synchronization, positioning, and mapping from diffuse multipath
WO2023010763A1 (en) Method and apparatus for determining positioning parameter, device, and storage medium
CN107271959B (en) Positioning method based on multi-antenna orthogonal frequency division multiplexing modulation signal
CN107171703B (en) Method capable of simulating propagation characteristics of fading signals in indoor multi-antenna communication system
JP2011014980A (en) Device and method of estimating radio wave propagation parameter
CN111693934A (en) Direction finding method and system based on digital interference system
Tayem et al. Hardware implementation of MUSIC and ESPRIT on NI-PXI platform
Meles et al. Measurement based performance evaluation of drone self-localization using AoA of cellular signals
Tonello et al. Radio positioning based on DoA estimation: An implementation perspective
CN108008351B (en) Signal positioning method based on forward attenuation reverse-push emission amplitude aggregation
Alibi et al. 2D DOA estimation method based on channel state information for uniform circular array
US8102297B2 (en) Method and system for deinterleaving signals
Zheng et al. Localization coverage analysis of THz communication systems with a 3D array
CN108828513B (en) Signal source positioning method based on intersection of electric wave propagation attenuation equal differential lines of multiple monitoring points
Lin et al. A real-time/low-cost DOA estimation system based on Rotman lens array antenna
CN116070082A (en) Polarization sensitive array amplitude-phase correction method based on spatial spectrum estimation
Amar et al. Direct position determination: A single-step emitter localization approach
Ahmed et al. Fast indoor localization using WiFi channel state information
Sun et al. 5g positioning based on the wideband electromagnetic vector antenna
Sedunova et al. Utilizing Beamsteering at Millimeter Waves for Indoor Object and Room Geometry Detection
Wu et al. Indoor single station 3D localization based on L-shaped sparse array
Zhen et al. DOA estimation for mixed signals with gain-phase error 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
GR01 Patent grant
GR01 Patent grant