CN114371492A - Interference situation analysis method based on terminal fusion positioning - Google Patents

Interference situation analysis method based on terminal fusion positioning Download PDF

Info

Publication number
CN114371492A
CN114371492A CN202111559192.6A CN202111559192A CN114371492A CN 114371492 A CN114371492 A CN 114371492A CN 202111559192 A CN202111559192 A CN 202111559192A CN 114371492 A CN114371492 A CN 114371492A
Authority
CN
China
Prior art keywords
interference
interference source
terminal
observation
power
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
Application number
CN202111559192.6A
Other languages
Chinese (zh)
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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202111559192.6A priority Critical patent/CN114371492A/en
Publication of CN114371492A publication Critical patent/CN114371492A/en
Pending legal-status Critical Current

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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses an interference situation analysis method based on terminal fusion positioning, and belongs to the field of satellite navigation. In the invention, the terminal receives the carrier-to-noise ratio information before and after interference to determine the position of the interference source; monitoring signal information fusion of the interference source by combining direction-finding monitoring equipment to obtain an accurate position of the interference source; and predicting and deducing the power of the interference source by the path loss model according to the position of the terminal and the received power to obtain the situation analysis of the interference source. The method and the system can integrate the terminal interference information and the interference information received by the monitoring station in a complex environment, and adjust the prediction of the position and the power of the interference source in real time, so that the situation presentation effect of the interference source is more excellent.

Description

Interference situation analysis method based on terminal fusion positioning
Technical Field
The invention belongs to the field of satellite navigation, and particularly relates to an interference situation analysis method based on terminal fusion positioning.
Background
In real-life production applications, the role of Global Navigation Satellite System (GNSS) is becoming more and more important. In recent years, the development of navigation technology presents the development situation of multi-system multi-level fusion. And if the navigation signal is interfered or can not be used, the navigation signal can bring great influence to the life of people and even cause serious loss. It is therefore very necessary and meaningful to study the navigation disturbance situation.
The navigation interference situation refers to that whether an interference source exists in the motion area range is judged by knowing surrounding geographic information and each received signal, the overall interference situation is known, and meanwhile, prior knowledge is provided for the inversion of the subsequent navigation signal environment, and convenience is provided for the subsequent navigation. The navigation interference situation supplements the traditional navigation situation which mainly identifies the target, and the accuracy of situation analysis is improved through learning of the situation level algorithm, so that the navigation interference situation is presented accurately. At present, most of the research fields of situation awareness are network systems, the navigation interference field is not provided with system research articles temporarily, most of the interference source positioning methods are interference source position positioning methods such as angle, arrival time and frequency and the like through observation quantity, the methods only pay attention to the position information of an interference source, interference source power information is not solved, and positioning accuracy still has a space for improvement.
Disclosure of Invention
In view of this, the invention provides an interference situation analysis method based on terminal fusion positioning, which has high precision, can solve the interference source position and power information, and can obtain the situation of the interference source by combining a signal propagation attenuation model.
In order to achieve the purpose, the invention adopts the technical scheme that:
an interference situation analysis method based on terminal fusion positioning comprises the following steps:
(1) determining a first presumed position of an interference source according to carrier-to-noise ratio change conditions of received signals of a plurality of observation terminals before and after interference;
(2) monitoring by direction-finding monitoring equipment to obtain interference source signal angle information, positioning a second presumed position of the interference source through direction-finding cross, and fusing the first presumed position and the second presumed position to obtain an accurate position of the interference source;
(3) according to the position and the received power of the observation terminal, the power of the interference source is estimated by a path loss model;
(4) and (4) obtaining the situation of the interference source according to the accurate position of the interference source obtained in the step (2) and the power of the interference source obtained in the step (3).
Further, the specific mode of the step (1) is as follows:
carrier-to-noise ratio variation u of received signal of observation terminal i before and after interferenceiComprises the following steps:
Figure BDA0003419994590000021
wherein the content of the first and second substances,
Figure BDA0003419994590000022
in order to achieve a carrier-to-noise ratio after interference,
Figure BDA0003419994590000023
is the interference front carrier to noise ratio; according to uiThe positive or negative of the value determines the interference type, uiGreater than 0 is suppression disturbance, ui< 0 is deceptive jamming;
calculating the weight w of the observation terminal ii：
Figure BDA0003419994590000024
n is the total number of the observation terminals;
the position of the observation terminal i is qi＝(xi，yi) And obtaining a first estimated position Pt of the interference source according to the position and the weight of each observation terminal:
Figure BDA0003419994590000025
further, in the step (2), a specific manner of fusing the first presumed position and the second presumed position is as follows:
carrying out weighted summation on the first presumed position and the second presumed position to obtain the accurate position of the interference source; wherein, the weight of the first guess position is α, the weight of the second guess position is 1- α, and the value of the weight α is as follows:
Figure BDA0003419994590000031
d is the average interval of the observation terminals, and the unit is meter.
Further, the specific mode of the step (3) is as follows:
the method comprises the following steps of defining a circular sampling range by taking the accurate position of an interference source as a center of a circle, wherein the radius r of the sampling range is selected by one of the following two modes:
1) taking one fourth of the radius of the distribution range of the observation terminal;
2) sequencing the observation terminals according to the sequence of the carrier-to-noise ratio variation values from large to small, finding the observation terminal positioned at a quarter of the sequence, and taking the distance from the observation terminal to an interference source as a r value;
and respectively averaging the received power and the interference-to-signal ratio of the observation terminal in the sampling range, and summing the two averages to obtain the presumed power of the interference source.
Further, the specific mode of the step (4) is as follows:
and (3) calculating the received power of each stepping point in the area range according to the accurate position of the interference source obtained in the step (2) and the power of the interference source obtained in the step (3) through a surrounding geographical environment and a path loss model to obtain the situation of the interference source in the area range, and displaying the situation and the received power of each observation terminal together in a visual form.
The invention has the beneficial effects that:
(1) the invention adopts a mode of fusing two positioning phases to obtain the accurate position of the interference source, can improve the positioning precision and simultaneously reduces the defect of a single positioning method.
(2) The method estimates the power of the interference source depending on the actual environment and the geographic condition, can invert the interference situation and perform visual display, and realizes the situation analysis function of the interference source.
Detailed Description
The technical solution of the present invention is further described below.
An interference situation analysis method based on terminal fusion positioning comprises the following steps:
(1) the terminal receives the carrier-to-noise ratio information before and after interference, and intelligently positions and determines the position of the interference source;
and comparing the satellite information received by the terminal before and after interference, wherein the carrier-to-noise ratio of the received signal is changed.
Setting the carrier-to-noise ratio variation value of the received signals of the terminal before and after interference as ui，
Figure BDA0003419994590000041
Wherein the content of the first and second substances,
Figure BDA0003419994590000042
in order to achieve a carrier-to-noise ratio after interference,
Figure BDA0003419994590000043
is the interference front carrier to noise ratio. Can be according to uiThe positive or negative of the value determines the interference type, uiGreater than 0 is suppression disturbance, ui< 0 is spoofing interference.
Calculating a weight w for each terminaliThe calculation method can be adjusted according to the specific surrounding environment, and the weighting mode core is to modify the weights by using the relative distance between the interference source and the receiver, for example, the weight of each sample point can be calculated from the carrier-to-noise ratio variation value according to the following formula:
Figure BDA0003419994590000044
each terminal position is qi＝(xi，yi) And weighting and calculating the position of the interference source according to different weights to obtain the position of the interference source as follows:
Figure BDA0003419994590000045
(2) the method comprises the following steps of combining interference source signal information monitored by direction-finding monitoring equipment, and obtaining the accurate position of an interference source by using two positioning fusion modes, wherein the specific modes are as follows:
and (3) measuring the angle information of the interference signal by the monitoring equipment, positioning the position of the interference source through direction finding and crossing, and performing data fusion on the position and the position positioned in the step (1). The fusion standard is judged according to the number of the terminals and the positioning accuracy of the terminals, the terminal positioning effect is the best within a certain range at the interval of the terminals, and the two positioning results use different weighting coefficients to obtain the final position coordinate result. If the average terminal interval is d meters and the terminal positioning weight is alpha, then
Figure BDA0003419994590000046
The positioning weight of the direction-finding monitoring equipment is 1-alpha.
(3) Predicting and deducing the power of an interference source by a path loss model according to the position of the terminal and the received power;
and (3) combining the position information of the interference source obtained in the step (1) and the position information of the interference source and the power intensity of the signal received by the terminal and the monitoring equipment with the geographic information and the surrounding environment, and realizing the power prediction of the interference source by selecting the terminal and utilizing the interfered interference-signal ratio and power information of the terminal.
Terminal average power PSAnd average interference-to-signal ratio JSRSThe calculation is shown in formula (1) and formula (2)
Figure BDA0003419994590000051
Figure BDA0003419994590000052
And taking the value of r as the number of terminal samples in a circle with the positioning result of the interference source as the center of the circle and r as the radius, wherein one method is to take a quarter value of an observation range, the other method is to sort the carrier-to-noise ratio change values, find the terminal at the quarter of the total carrier-to-noise ratio change values from large to small, and take the distance from the terminal at the position to the interference source as the value of r. PjFor post-interference received power of terminal samples, JSRjAfter it interferes withThe interference-to-signal ratio.
The interference source power prediction result PN is shown in formula (3):
PN＝PS+JSRS(3)
(4) and obtaining the situation analysis of the interference source according to the position and the power of the interference source deduced by prediction.
And obtaining the position and power information of the interference source through the first three steps, and calculating the power value of each reference point in a certain area range, namely the received signal strength of each stepping point in the area range through the surrounding geographic environment and the path loss model. And combining the power value of the terminal after the previous interference to perform interface display, and realizing the inversion and visualization of the navigation interference situation.
In an actual navigation environment, a signal observation range is not an ideal environment, and the situation that an interference source exists can be met with a certain probability, so that a signal received by a terminal is interfered, and the navigation effect is not ideal. In order to better determine the interference degree and the navigation positioning precision, the interference situation perception technology is researched and analyzed.
The interference situation analysis method based on the terminal fusion positioning provided by the invention fully utilizes the information obtained by the terminal and the monitoring equipment and the properties of the signal, and can obtain a more accurate interference situation inversion effect by using an intelligent positioning algorithm, fusion positioning and power prediction.
In a word, the method can integrate the terminal interference information and the interference information received by the monitoring station in a complex environment, and adjust the prediction of the position and the power of the interference source in real time, so that the situation presentation effect of the interference source is more excellent.

Claims (5)

1. An interference situation analysis method based on terminal fusion positioning is characterized by comprising the following steps:
(1) determining a first presumed position of an interference source according to carrier-to-noise ratio change conditions of received signals of a plurality of observation terminals before and after interference;
(2) monitoring by direction-finding monitoring equipment to obtain interference source signal angle information, positioning a second presumed position of the interference source through direction-finding cross, and fusing the first presumed position and the second presumed position to obtain an accurate position of the interference source;
(3) according to the position and the received power of the observation terminal, the power of the interference source is estimated by a path loss model;
(4) and (4) obtaining the situation of the interference source according to the accurate position of the interference source obtained in the step (2) and the power of the interference source obtained in the step (3).
2. The interference situation analysis method based on terminal fusion positioning according to claim 1, wherein the specific manner of the step (1) is as follows:
carrier-to-noise ratio variation u of received signal of observation terminal i before and after interferenceiComprises the following steps:
Figure FDA0003419994580000011
wherein the content of the first and second substances,
Figure FDA0003419994580000012
in order to achieve a carrier-to-noise ratio after interference,
Figure FDA0003419994580000013
is the interference front carrier to noise ratio; according to uiThe positive or negative of the value determines the interference type, uiGreater than 0 is suppression disturbance, ui< 0 is deceptive jamming;
calculating the weight w of the observation terminal ii：
Figure FDA0003419994580000014
n is the total number of the observation terminals;
the position of the observation terminal i is qi＝(xi，yi) Obtaining a first estimated position P of the interference source according to the position and weight of each observation terminalt：
Figure FDA0003419994580000015
3. The interference situation analysis method based on terminal fusion positioning according to claim 1, wherein in the step (2), the specific way of fusing the first presumed position and the second presumed position is as follows:
carrying out weighted summation on the first presumed position and the second presumed position to obtain the accurate position of the interference source; wherein, the weight of the first guess position is α, the weight of the second guess position is 1- α, and the value of the weight α is as follows:
Figure FDA0003419994580000021
d is the average interval of the observation terminals, and the unit is meter.
4. The interference situation analysis method based on terminal fusion positioning according to claim 1, wherein the specific manner of step (3) is as follows:
the method comprises the following steps of defining a circular sampling range by taking the accurate position of an interference source as a center of a circle, wherein the radius r of the sampling range is selected by one of the following two modes:
1) taking one fourth of the radius of the distribution range of the observation terminal;
2) sequencing the observation terminals according to the sequence of the carrier-to-noise ratio variation values from large to small, finding the observation terminal positioned at a quarter of the sequence, and taking the distance from the observation terminal to an interference source as a r value;
and respectively averaging the received power and the interference-to-signal ratio of the observation terminal in the sampling range, and summing the two averages to obtain the presumed power of the interference source.
5. The interference situation analysis method based on terminal fusion positioning according to claim 1, wherein the specific manner of step (4) is as follows:
and (3) calculating the received power of each stepping point in the area range according to the accurate position of the interference source obtained in the step (2) and the power of the interference source obtained in the step (3) through a surrounding geographical environment and a path loss model to obtain the situation of the interference source in the area range, and displaying the situation and the received power of each observation terminal together in a visual form.
CN202111559192.6A 2021-12-20 2021-12-20 Interference situation analysis method based on terminal fusion positioning Pending CN114371492A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111559192.6A CN114371492A (en) 2021-12-20 2021-12-20 Interference situation analysis method based on terminal fusion positioning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111559192.6A CN114371492A (en) 2021-12-20 2021-12-20 Interference situation analysis method based on terminal fusion positioning

Publications (1)

Publication Number Publication Date
CN114371492A true CN114371492A (en) 2022-04-19

Family

ID=81139417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111559192.6A Pending CN114371492A (en) 2021-12-20 2021-12-20 Interference situation analysis method based on terminal fusion positioning

Country Status (1)

Country Link
CN (1) CN114371492A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115480475A (en) * 2022-09-02 2022-12-16 北京无线电计量测试研究所 Time service interference detection method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115480475A (en) * 2022-09-02 2022-12-16 北京无线电计量测试研究所 Time service interference detection method
CN115480475B (en) * 2022-09-02 2023-11-28 北京无线电计量测试研究所 Time service interference detection method

Similar Documents

Publication Publication Date Title
Wielgosz et al. Regional ionosphere mapping with kriging and multiquadric methods
CN106908821B (en) Indoor and outdoor seamless positioning switching method
CA2808155C (en) Adaptive method for estimating the electron content of the ionosphere
CN104635203B (en) Radio interference source direction-finding and positioning method based on particle filter algorithm
CA2866496A1 (en) Location agent geofence
CN106226732B (en) The indoor wireless positioning and tracing method filtered based on TOF and iteration without mark
CN103152745A (en) Method of locating mobile node with strong adaptivity
CN114371492A (en) Interference situation analysis method based on terminal fusion positioning
CN104507097A (en) Semi-supervised training method based on WiFi (wireless fidelity) position fingerprints
CN114417580A (en) Method for evaluating influence of observation system on assimilation performance of global ionosphere data
Zhang et al. Ultra wideband positioning technology for accident ships under adverse sea condition
Huang et al. Analysis and improvement of ionospheric thin shell model used in SBAS for China region
CN108650629A (en) A kind of indoor three-dimensional location based on radio communication base station
CN111932017A (en) Short-term forecasting method suitable for single-station ionized layer TEC
CN115524661A (en) Short wave time difference positioning method for joint optimization of height of ionized layer and target position
CN111781570A (en) Radar online precision analysis method based on real-time ADS-B data
CN101777958B (en) Method for forecasting group delays within certain range near receiving point
AU2014235257A1 (en) Directional pruning of transmitters to improve position determination
Groves et al. Intelligent urban positioning, shadow matching and non-line-of-sight signal detection
Yang et al. Geometric dilution of precision for far-distance TDOA location of shortwave
CN115857058A (en) Ionosphere data analysis model construction method and terminal thereof
CN111427076B (en) Vehicle performance evaluation method and system based on dual-mode satellite information deep fusion
CN108738129A (en) A kind of base station positioning method based on drive test data
CN102170317B (en) Method for analyzing spatial randomness of group path adjacent to receiving point
Lv et al. Research on shadow matching algorithm based on consistency probability weighting

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