CN110677865B - Method for positioning external interference source of mobile communication network - Google Patents
Method for positioning external interference source of mobile communication network Download PDFInfo
- Publication number
- CN110677865B CN110677865B CN201910909518.XA CN201910909518A CN110677865B CN 110677865 B CN110677865 B CN 110677865B CN 201910909518 A CN201910909518 A CN 201910909518A CN 110677865 B CN110677865 B CN 110677865B
- Authority
- CN
- China
- Prior art keywords
- ray
- grid
- interference
- point
- interference source
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/08—Testing, supervising or monitoring using real traffic
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/345—Interference values
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/10—Scheduling measurement reports ; Arrangements for measurement reports
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention discloses a mobile communication network external trunkA disturbance source positioning method belongs to the field of mobile communication. Firstly, confirming a detection region R, rasterizing a three-dimensional space, recording the spatial distribution of terrain and ground objects in the region R, and realizing scene modeling. Then reasonably planning a drive test track and manually carrying out drive test according to the measurement radius of the drive test sweep frequency measurement terminal and the distribution condition of interfered points in the region R to obtain an initial measurement result set DTSet, selecting N measurement points from the DTSet to form a back tracking starting point set MPSet, reversely sending a primary signal ray to each measurement point in the MPSet according to the signal coverage range of the measurement point to generate possible interference signal propagation paths PathSet (MPSet), carrying out three-dimensional modeling on each propagation path in the MPSet, and constructing a grid-path index IgtAnd screening a grid set CanGSet serving as an interference source candidate position from the initial grid set, and selecting a grid corresponding to the optimal target value as an interference source positioning grid. The invention improves the positioning accuracy and reduces the manpower and material resource investment.
Description
Technical Field
The invention belongs to the field of mobile communication, and particularly relates to a method for positioning an external interference source of a mobile communication network.
Background
In a complex urban environment, as the wireless propagation environment of the network is increasingly complex, multiple networks are deployed in an overlapping manner, network interference seriously affects the information transmission quality and user experience.
The intra-network interference is mainly same frequency and adjacent frequency interference, and is generated by other network element equipment with the same or similar parameter (such as frequency in L TE network and Physical Cell Identification (PCI)) setting in the network, the position and parameter configuration of the network element equipment are known, the extra-network interference is stray, blocking or intermodulation interference and the like generated by network elements of other networks outside the network, such as interference generated by W L AN network to L TE network, or independent interference equipment outside the network, such as various signal shields.
In the interference positioning optimization work, firstly, the interference is found and the interference type is analyzed, then, the position of an interference source is positioned, and finally, the interference source is processed to eliminate or weaken the interference. Currently, the intra-network interference positioning optimization can achieve targets by analyzing and adjusting network configuration parameters, such as optimization and adjustment of frequency and PCI parameters, and a mature and practical processing method is available at present.
For the optimization of the positioning of the interference outside the network, the interference existing in the network can be found through means such as DT (Driving test) drive test, frequency sweep and the like or according to the performance index analysis of a network management system. However, since the deployment position of the network element or device serving as the off-network interference source is unknown, or the deployment position is known but the configuration parameters are unknown, how to quickly and accurately locate the position of the interference source is always a big difficulty.
The currently and practically adopted method for positioning the external interference of the mobile communication network mainly comprises the steps of collecting interference signals by using a frequency spectrograph, a directional antenna and other equipment in an area with network interference, analyzing the frequency spectrum characteristics of the interference signals and judging the interference types and directions by means of drive test. Then, combining the strength and direction of the interference signal, adopting a three-point positioning method or other geometric positioning methods to gradually reduce the search range and gradually realize the positioning of the interference source along the direction of the interference source.
The interference positioning method has the following problems in practical application: 1) in a complex urban environment, a large amount of multipath non-line-of-sight propagation exists in interference signal propagation, and a simple triangulation method assumes that the interference signal is propagated linearly and is inconsistent with the actual situation, so that a large positioning error exists; 2) although the positioning method of continuous measurement and gradual search can achieve certain positioning accuracy, the consumed manpower and material resources are large.
Disclosure of Invention
The invention provides an external interference source positioning method of a mobile communication network based on interference signal propagation path analysis, which aims at mobile communication networks which are deployed and operated in dense urban areas, including but not limited to L TE cellular networks, wireless private networks and the like, and considers the influence of terrain and ground objects in the environment on wireless signal propagation, so as to reduce the investment of road test and frequency sweeping manpower and material resources while ensuring that the ideal positioning accuracy is achieved.
The method comprises the following steps:
step one, aiming at an area R covered by a mobile communication network0According to the interference information counted by the base station, the drive test frequency sweep analysis result and the user complaint information, confirming a detection area R of the interference source;
and secondly, carrying out three-dimensional space rasterization on the detection region R of the interference source, and recording the spatial distribution of terrain and ground objects in the search region R of the interference source by adopting a uniform grid acceleration structure to realize scene modeling.
And step three, reasonably planning a drive test track and manually implementing drive test according to the measurement radius of the drive test sweep frequency measurement terminal and the distribution condition of interfered points in the region R to obtain an initial measurement result set DTset.
An initial measurement result set DTSet ═ { mp };
the information recorded by each measuring point mp comprises a measuring time mp.time, a measuring point longitude mp.lat, a measuring point latitude mp.lon, an interference signal strength mp.str and a frequency band mp.freq received at the measuring point.
And step four, analyzing the initial measurement result set DTSet, selecting N measurement points from the initial measurement result set DTSet, calculating the signal angle of each measurement point, using the signal angle as a starting point of the back ray tracking, and forming a back tracking starting point set MPSet.
Step five, sequentially selecting each measuring point from the back tracking starting point set MPSet, and reversely sending rays according to the signal coverage angle range of the measuring point to generate a possible interference signal propagation reverse ray path PathSet (MPSet);
the specific single signal back ray tracing steps are as follows:
step 501, based on the current emission starting point mpiEmitting a primary ray r0And adding the information nodeinfo of the transmission starting point into a data structure ray L ist0；
Step 502, judging the currently tracked ray rkWhether k is directly into the top surface of the region R orIf the ray exceeds the area range, the current ray tracing is ended, the end point is the intersection point of the ray and the boundary surface, and the end point information nodeinfo is recordedeAnd added to ray L ist to obtain a complete slave mpiThe signal propagation path starting at the end point and ending at the end point, and the flow proceeds to step 508; otherwise, go to step 503;
step 503, judging the currently tracked ray rkIf the number of reflections and diffractions occurring in the previous propagation process exceeds the predetermined threshold RDmax, the current ray tracing process is ended and the process proceeds to step 508. Otherwise, go to step 504;
step 504, ray r is performed by using a uniform grid acceleration structurekAnd fast collision detection with the topographic asperities in the interference source search region R.
Step 505, judge ray rkWith the current uniform grid UG in the forward directionuWhether a collision is detected by the intersection of the inner building bid and the triangle tid; if so, go to step 506, otherwise, the current ray rkContinue to propagate to the next uniform grid in the forward direction, go back to step 502;
step 506, find ray rkThe collision point with the terrain bump is taken as a ray rkEnd point r ofkEp, r iskEp information nodeinfokAdding the data into a data structure ray L ist, judging whether the collision point is on a terrain TIN surface, if so, terminating the current ray tracing, and entering step 508, otherwise, entering step 507 when the collision point is on the building surface;
step 507, judging whether the reflection or diffraction is performed according to the position of the intersection point, and generating a corresponding secondary reflection line or a corresponding diffraction line; returning to step 502 to continue tracking and judging a new reflection line rk+1Or around a ray rk+1Whether the area is out of range.
The method specifically comprises the following steps: if the intersecting collision point is the side surface or the top surface of the building, the next-stage reflection line r is generatedk+1Setting a ray rk+1Is of the type reflection, ray origin rk+1Sp is rkEp and according to the primary ray rkDirection calculation ofOut of the reflected ray rk+1Direction along which the collision-generated reflected ray r continues to be tracedk+1Returning to step 502 to judge whether the area is out of range again; otherwise, the position of the crossed collision point is positioned on the edge of the building to generate the next-level ray rk+1Setting a new section of ray rk+1Is a diffraction, ray origin rk+1Sp is rkEp and according to the primary ray rkThe direction of (a) is calculated around the ray rk+1Direction along which the ray-around r generated by the collision is continuously tracedk+1Then, go back to step 502 to determine again whether the area is out of range.
Step 508, aiming at nodeinfo information in the data structure ray L ist, constructing mp from the starting pointiStarting from a plurality of rays r0,r1,...,rRDmaxForming an interference signal propagation path, and calculating the propagation parameters of the path;
the propagation parameters include reflection/diffraction loss coefficient, each segment rkThe number of ray orders, ray rkAnd k is more than or equal to 1 and less than or equal to RDmax at the starting point and the end point.
One propagation path is denoted as traj ═ r0,r1,...,rRDmax＞。
Tracking all interference measurement points mp in the set of origin MPSets for the back-raysi∈ MPSet, a propagation path traj (mp) sequentially constructed from a data structure ray L isti) The set of possible propagation traces constituting the interfering signal is:
PathSet(MPSet)＝{traj(mpi),mpi∈MPSet}。
step six, modeling the three-dimensional rays in the possible propagation path set PathSet (MPSet) of the interference signal, and constructing a grid-path index Igt；
The trellis-path index is represented as: i isgt＝{Gr|trInf1,trInf2,...,trInfNt}；
IgtRecording individual grids G within an interferer search region RrTrajectory information of all paths/trajectories covered over which interference signals may propagate; trInfiRepresenting a propagation pathtrack information of traj;
for a certain propagation path traj, traversing the ray r in sequence0,r1,...,rRDmaxExperienced individual grids GrAnalyzing the recording track information trInfiGenerating a raster-path index IgtIndex record of < Gr|trInfi＞；
Step seven, indexing I by analyzing the grid-pathgtAnd screening out a grid set CanGSet serving as a candidate position of the interference source, calculating a plurality of target values of grids of each candidate position, and selecting a grid corresponding to the optimal target value as an interference source positioning grid.
The screening and filtration mode is as follows:
1) filtering and screening out backward ray tracing starting point mpi∈ MPSet.
2) Filtering out signal starting point emission point mp from which covered rays comeiThe number of the grids is less than a set threshold From _ num;
3) and filtering out grids with covered strong signal ray ratios smaller than a set threshold S _ rate.
The process of selecting the grid corresponding to the optimal target value as the interference source positioning grid specifically comprises the following steps:
first step, combine trellis-path index IgtComputing each grid G of the candidate grids CanGSetrNumber of covered strong signal ray paths GrPathnum, integrated interference signal strength difference GrPathloss, number of traceback starting points G from which the grid covered rays comer.fromnum；
Second, respectively solving the number G of covered strong signal ray pathsrMaximum Path of pathnummaxIntegrated signal strength difference GrMinimum value of pathloss P L ossminAnd the number of covering rays from the back-tracking starting point GrMaximum From of Frommax。
Thirdly, solving the median Path of each sub-target according to the scales of the three target valuesmid，PLossmid，FrommidFor reducing the result of weighting the dimensionsAnd (3) modifying each sub-target function into:
wherein f is1(Gr) Is a grid GrThe strong signal ray constraint value of (f) is calculated2(Gr) Is a grid GrThe integrated signal interference difference value constraint value of f3(Gr) Is a grid GrThe signal source constraint value computation function.
Fourthly, setting a weighting coefficient w according to the importance degree of each sub-targetiAnd converting the multi-objective optimization problem into a single-objective optimization problem, finishing the evaluation and positioning of the interference source, and finally obtaining the single-objective optimization problem by the following formula:
and fifthly, calculating an evaluation value of each grid through a single objective function in the grid set CanGSet, wherein the grid with the minimum evaluation value is the interference source positioning grid.
The invention has the advantages that:
1) the method is oriented to dense urban areas, and provides a method for positioning an unknown interference source outside the network based on a reverse ray tracing path aiming at an actual TD-L TE network which is deployed and operated and considering the influence of the environment on signal propagation, so that the ideal positioning accuracy is achieved and the investment of manpower and material resources is reduced.
2) The method for positioning the interference source outside the mobile communication network comprises the steps that protruding terrain objects such as slopes, buildings and the like exist in an interference search area, and interference signals are directly transmitted, reflected and diffracted in the process of propagation, so that multiple paths from the interference source can be received by one position point. The invention starts from the path itself, proposes from the strong signal path measuring points with wide distribution and less shielding, simulates the propagation of signals, generates the reverse propagation path of the signals, and realizes the positioning of the interference source by analyzing the distribution of the paths and the corresponding propagation loss.
3) The invention provides a method for positioning an interference source outside a mobile communication network, and aims to avoid positioning the position of the interference source according to the number of covered ray paths.
Drawings
Fig. 1 is a flow chart of a method for locating an external interference source of a mobile communication network according to the present invention;
FIG. 2 is a schematic diagram of the present invention for determining the detection range of an interference source;
FIG. 3 is a flow chart of the present invention for generating possible interfering signal propagation back ray paths;
FIG. 4 is a schematic diagram of a two-dimensional x-y planar ray map of the present invention;
FIG. 5 is a flow chart of a grid filtering algorithm of the present invention;
fig. 6 is a schematic diagram illustrating the principle of selecting a grid corresponding to an optimal target value as an interference source positioning grid according to the present invention.
Detailed Description
The following describes embodiments of the present invention in detail and clearly with reference to the examples and the accompanying drawings.
A large number of ground objects such as buildings and the like exist in a complex dense urban area, and a slope and other convex terrains exist, so that the interference signal propagation has Non-line-of-Sight (N L OS, Non-L ine of Sight) propagation such as reflection and diffraction, and meanwhile, the interference signal propagation has a multipath phenomenon, namely the interference signal reaches the same path measuring point through a plurality of propagation paths.
The method positions the position of the interference source according to the signal intensity of the interference signal obtained by drive test sweep frequency, provides a reverse ray tracking model through the signal in consideration of non-line-of-sight and multipath propagation in the process of the propagation of the interference signal, simulates the possible propagation paths of the interference signal, analyzes and positions the position of the off-network interference source according to the propagation paths, and is specifically completed by the following three stages:
the first stage is as follows: selecting a starting point set for carrying out reverse ray tracking on the signal from strong signal sweep frequency drive test measuring points;
the interference signal propagation path according to which the position of the interference source is positioned is generated by the backward ray tracing, and the backward ray tracing starting point is the selected sweep frequency drive test measuring point.
And a second stage: possible propagation paths for interfering signals are generated by back ray tracing.
In the stage, an acceleration structure such as a uniform grid and a KD-tree is adopted to model objects such as convex terrains such as slopes and buildings in an interference source search area, and the spatial distribution of the terrains and the objects is recorded. Then, each sweep road test measurement point mp in the set of ray tracing origin obtained from the first stagei∈ MPSet, reversely sending out signals according to the signal coverage angle range of the measuring point, generating interference signals by reverse ray tracing from possible interference source position to the point mpiPossible signal propagation paths trajiThe propagation path set PathSet is formed as { traj ═ traji}；
And a third stage: and positioning the position of the interference source step by step according to multi-target constraint.
Using interference signal propagation paths traj obtained in the second stagei∈ PathSet, screening out a plurality of interference source candidate positions loc according to 3 target constraint conditions of interference source positioningkForming a set of candidate positions CanGSet ═ { loc ═ loc of interference sourcesk}; then, each possible candidate location loc is evaluatedkAnd picking out the candidate position loc ∈ CanGSet with the optimal score as the positioning result of the off-network interference source.
The TD-L TE network is taken as an example to illustrate the specific content of the invention, but the positioning method of the invention is also applicable to other mobile cellular networks, such as FDD-L TE, WCDMA, and wireless private networks such as power wireless private network, trunking system and the like.
As shown in fig. 1, the method comprises the following steps:
step one, aiming at an area R covered by a mobile communication network0According to the interference information counted by the base station, the drive test sweep frequency analysis result and the user complaint information confirmation area R0Defining the detection range R of the interference source according to the position of the interference;
the method comprises the following specific steps:
as shown in fig. 2, firstly, according to the base station interference statistical information provided by the network management system, the interfered degree of the interfered cell/base station is evaluated, and the interfered degree exceeding the base station interference threshold thld is screenedeintThe cell/base station as the first kind of interfered point is marked as Ib1，Ib2Etc.;
the base station interference statistic information comprises: background noise, PRB interference, etc.;
then, the sweep frequency drive test data obtained in advance are analyzed, and interference signal intensity exceeding the sweep frequency drive test interference threshold value thld is screened outdtThe sweep frequency drive test measuring point is taken as a second type interfered point and is marked as Id1，Id2Etc.;
continuously analyzing the position of the complaint point and the interference severity degree according to the complaint data of the user, thereby determining the position point which is seriously interfered in the network as a third type interfered point which is marked as Iu1，Iu2Etc.;
finally, constructing a rectangular bounding box R according to the positions of the three types of interference points1Will enclose the box R1And expanding the set distance He to the periphery to obtain a new rectangle R, namely the interference source searching area.
And secondly, carrying out three-dimensional space rasterization on the detection search region R of the interference source, and recording the spatial distribution of terrain and ground objects in the detection search region R of the interference source by adopting a uniform grid acceleration structure to realize scene modeling.
In order to calculate an interference signal propagation track by utilizing reverse ray tracking and shield and collision of a wireless signal and a terrain feature in a propagation process, three-dimensional space rasterization is carried out on an interference source searching area R, and a uniform grid acceleration structure is adopted to record the spatial distribution of the terrain feature in the interference source searching area R, so that modeling of the terrain feature, namely modeling of the terrain feature, is realized.
The method comprises the following specific steps:
first, the minimum longitude minlon and the minimum latitude minlat in the interference source search region R are used as the origin of coordinates, and the three-dimensional space in which the interference source search region R is located is divided into Lw×Lw×LhA rectangular parallelepiped grid;
one-dimensional number of grid is GrR ∈ N is the one-dimensional number of the grid, N is the natural number set, and the search region R is regarded as the grid GrComposition, expressed as R ═ GrGrid G, r ∈ NrThe x-y-z axis three-dimensional discrete coding coordinate representation is adopted, namely G (gx, gy, gz), gx, gy, gz ∈ N, and the conversion relation between the one-dimensional coordinate r and the three-dimensional coordinate (gx, gy, gz) is rho: r → (gx, gy, gz), mu (gx, gy, gz) → r.
Correspondingly, the ground two-dimensional plane in which the interference source search region R is located is divided into the side lengths of LwA square grid of (a);
then, downloading grid type discrete digital high-rise model DEM data which represent terrain height in the interference source searching region R, and converting the DEM data into vectorized irregular triangular network TIN;
the TIN is formed by a set of triangles trian representing the topographic relief within the region RtComposition is carried out; denoted as tin (r) ═ triant,t∈N}；t＝0,1,2,...。
And continuing to complete building modeling according to the building information in the region R provided by the three-dimensional electronic map.
Buildings are identified by bid, and a building is represented by the set of vertices { bvretex } of the building floor outline polygon and its building height bhight.
Finally, an acceleration structure represented in the form of a uniform grid is constructed for increasing the speed of the back ray tracing collision detection.
The identification bid for the building satisfies the following condition: the bottom surface polygon of identity bid contains a uniform grid UGuOr marking the sides of the bottom polygon of bid and the uniform gridGrid UGuIntersecting;
the identification tid for the TIN triangle satisfies the following condition: the projected triangle of the identification tid on the two-dimensional ground plane contains a uniform grid UGuOr identifying the sides of the projected triangle of tid on the two-dimensional ground plane and the uniform grid UGuIntersecting;
the uniform grid acceleration structure I is formed by the identification points meeting the conditionsug＝{UGu|bid,tid}。
The scene modeling is used for the backward ray tracing part in one part and spatial rasterization in the other part for positioning the interference source part.
Step three, in the interference source detection region R, measuring the radius R of the measuring terminal according to the sweep frequency drive testmeaAnd the distribution condition of interfered points in the region R, reasonably planning a drive test track and manually implementing the drive test to obtain an initial measurement result set DTset.
The interfered points comprise: a cell/base station interfered point, a drive test interfered point and a user complaint interfered point.
Firstly, in an area R where an interference source is located, a plurality of planning measurement points ap are arranged along a road on a three-dimensional map to form a planning measurement point set DPSet ═ ap };
the setting of each planned measurement point ap should guarantee: for any point p in the region R, at least N is found in the set DPSetmeaMore than or equal to 3 measurement points are planned, such that points p and NmeaThe distances between the planned measuring points do not exceed the measuring radius rmeaI.e. p is in NmeaWithin the measurement range of each planned measurement point;
then, according to the DPSet set of the planning measurement points, a sweep frequency drive test line is planned, and the following conditions are ensured: 1) the sweep frequency drive test line passes through all the planning measurement points ap in the set DPSet; 2) the path length of the sweep frequency drive test line is as small as possible;
finally, performing drive test on the planned sweep frequency drive test line by using measurement equipment such as a sweep frequency drive test instrument, and forming an initial measurement result set DTset (mp);
the information recorded by each measuring point mp comprises a measuring time mp.time, a measuring point longitude mp.lat, a measuring point latitude mp.lon, an interference signal strength mp.str and a frequency band mp.freq received at the measuring point.
Step four, analyzing the initial measurement result set DTSet, selecting N measurement points from the initial measurement result set DTSet, analyzing the signal approximate direction of the measurement points according to the road measurement point signals to obtain the measurement point signal center angle, using the measurement point signal center angle as the starting point of the back ray tracking, and forming a back tracking starting point set MPSet ═ { mp ═ mp { (mp)i}。
In order to ensure that the reverse ray tracing generates as many propagation paths passing through the real position of the interference source as possible so as to improve the interference positioning precision of the subsequent stage, according to the following conditions, part of sweep frequency path measuring points mp are selected from a sweep frequency path measuring point set DTset ═ mp } points obtained by sweep frequency path measuringiAs the starting point of the back ray tracing and form a back ray tracing starting point set
The following conditions are satisfied: 1) measuring point mpiThe received interference signal is stronger; 2) measuring point mpiLittle signal obstruction brought by surrounding buildings and sloping fields, i.e. mpiShielding less at the periphery; 3) the positions of the measurement points in the MPSet are distributed as widely as possible.
The method comprises the following specific steps:
step 401, extracting the interference signal strength threshold value thld with the signal strength higher than the set signal strength threshold value from the drive test measurement point set DTSet ═ { mp }strMeasuring point mpiForming a candidate set of traceback start points
mpi.str＞thldstrIllustrates the arrival of an interference signal from an interference source at point mpiHas small signal loss and point mpiThe distance between the interference source and the interference source is short;
step 402, according to the spatial distribution of each measurement point in the candidate back tracking starting point set COPset1, some measurement points which are seriously shielded by the terrain such as surrounding buildings or the convex terrain such as sloping fields and are not suitable for being used as the back ray tracking starting points are removed from the candidate back tracking starting point set COPset 2.
Less signal shielding is brought by convex terrains such as buildings, ground objects and sloping fields around the measuring point in the set COPSEt2, and the generation of more non-line-of-sight N L OS propagation paths in the process of back ray tracking is facilitated, so that a larger positioning error is caused;
step 403, according to the spatial distribution of the measurement points in the COPset2, selecting N measurement points according to the following 2 conditions to form a back tracking start point set
1) For arbitrary 2 measuring points mpi,mpj∈ MPSet, i ≠ j, distance dist (mp) between these two pointsi,mpj) Greater than a measurement point distance threshold thlddisI ≠ j, i.e. dist (mp)i,mpj)＞thlddis；
2) Set MPSet ═ mpiThe positions of the measuring points in the device are distributed as widely as possible, and the measuring points are prevented from being distributed in a narrow and long area, such as a straight road.
Specifically, a rectangular bounding box rec (MPSet) is required to be made for all measurement points in the MPSet, and the aspect Ratio of the bounding box rec (MPSet) needs to meet an upper and lower limit threshold value Ratio L W ∈ [ thld [, ]RL,thldRH]；
Wherein thldRLLower threshold of aspect ratio, thldRHIs the aspect ratio upper threshold.
Step 404, analyzing the set of back tracking start points MPSet ═ mpiEach measuring point mpiAnd the signal size of the peripheral road measuring points, calculating mpiThe center azimuth angle of the back ray tracing signal emission of (1);
step five, sequentially selecting each measuring point from the back tracking starting point set MPSet, and reversely emitting rays to an azimuth angle range [ angle-delta, angle + delta ] according to the signal angle of the measuring point to generate a possible interference signal propagation back ray path PathSet (MPSet);
where Δ is the signal coverage angle.
In dense urban areas, due to the shielding, reflection and diffraction of wireless signals by topographic features such as sloping fields, buildings and the like, multipath and non-line-of-sight propagation exist in interference signal propagation. Thus, considering the analog interference signal propagation process, from the selected strong interference measurement point mpi∈ MPSet starts, reversely sends out signals, simulates the interference signal propagation process, and positions the specific position of the interference source by analyzing the interference signal propagation path.
The method realizes the generation of possible interference signal propagation paths based on a ray tracking model, sets the whole interference source search area R as a cuboid to interfere with a measurement point mpi∈ MPSet is used as a starting point, and a back-tracing primary ray is emitted to the surface of a building in the region R, the boundary surface and the top surface of the region R, and then the traces of the primary ray and the subsequent secondary ray are traced, so that the propagation path of the interference signal is obtained.
From a single disturbance measurement point mpiEach propagation path obtained by ray tracing in the back tracing process of an emitted primary ray is recorded in a data structure ray L ist < Nodeinfo >iAn ordered set of reflected or diffracted collision points, and ray end point information; and each propagation path consists of a plurality of sections of ordered rays, and each section of rays is a directed line segment with a three-dimensional space R starting point of r.sp and an end point of r.ep.
As shown in fig. 3, the specific ray tracing steps are as follows:
step 501, based on the current emission starting point mpiEmitting a primary ray r0And adding the information nodeinfo of the transmission starting point into a data structure ray L ist0；
r0Starting point r of0Sp is mpiI.e. r0.sp＝mpi；mpiStr is taken as the initial signal strength of the signal source;
step 502, judging the currently tracked ray rkIf the k is directly on the top surface of the region R or exceeds the region range, the current ray tracing is ended and the process proceeds to step 508; otherwise, go to step 503;
when the current ray tracing is finished, the end point is the intersection point of the ray and the boundary surface, and end point information nodeinfo is recordedeAdding to ray L ist to obtain a complete slave mpiA signal propagation path starting at the end point and ending at the end point.
Step 503, judging the currently tracked ray rkIf the reflection and diffraction times in the previous propagation process exceed the specified threshold value RDmax, the current ray tracking process is ended and the process goes to step 508; otherwise, go to step 504;
exceeding the specified threshold value RDmax indicates that the interference signal has been attenuated too much during the propagation along the ray trajectory/path, and it makes no sense to continue tracking the ray generation propagation path.
Step 504 of accelerating the structure I with a uniform gridug＝{UGu| bid, tid } making ray rkAnd fast collision detection with the topographic asperities in the interference source search region R.
Step 505, judge ray rkCurrent uniform grid UG in forward directionuWhether or not to interact with the ray rkDetecting a collision; if so, go to step 506, otherwise, the current ray rkContinue to propagate to the next uniform grid in the forward direction, go back to step 502;
the collision being, in particular, with rays rkAnd a uniform grid UGuIntersection test of inner building bid and triangle tid.
Step 506, find ray rkThe collision point with the terrain bump is taken as a ray rkEnd point r ofkEp, r iskEp information nodeinfokAdding the data into a data structure ray L ist, judging whether the collision point is on a terrain TIN surface, if so, terminating the current ray tracing, and entering step 508, otherwise, entering step 507 when the collision point is on the building surface;
step 507, judging whether the reflection or diffraction is performed according to the position of the intersection point, and generating a corresponding secondary reflection line or a corresponding diffraction line; returning to step 502 to judge the new reflection ray r againk+1Or around a ray rk+1Whether it exceeds the range of the areaAnd (5) enclosing.
The method specifically comprises the following steps: if the intersecting collision point is the side surface or the top surface of the building, the next-stage reflection line r is generatedk+1Setting a ray rk+1Is of the type reflection, ray origin rk+1Sp is rkEp and according to the primary ray rkThe direction of (d) is calculated to obtain a reflection ray rk+1Direction along which the collision-generated reflected ray r continues to be tracedk+1Returning to step 502 to judge whether the area is exceeded; otherwise, the position of the crossed collision point is positioned on the edge of the building to generate the next-level ray rk+1Setting a new section of ray rk+1Is a diffraction, ray origin rk+1Sp is rkEp and according to the primary ray rkThe direction of (a) is calculated around the ray rk+1Direction along which the ray-around r generated by the collision is continuously tracedk+1Returning to step 502 to judge whether the area is out of range again.
Step 508, aiming at nodeinfo information in the data structure ray L ist, constructing mp from the starting pointiStarting from a plurality of rays r0,r1,...,rRDmaxForming an interference signal propagation path, and calculating the propagation parameters of the path;
the propagation parameters include reflection/diffraction loss coefficient, each segment rkThe number of ray orders, ray rkAnd k is more than or equal to 1 and less than or equal to RDmax at the starting point and the end point.
A start point mpiThe path of travel from is denoted traj ═ r0,r1,...,rRDmax> (ii). Wherein the path starting point traj.sp ═ r0.sp＝mpi。
Tracking all interference measurement points mp in the set of origin MPSets for the back-raysi∈ MPSet, a propagation path traj (mp) sequentially constructed from a data structure ray L isti) The set of possible propagation trajectories for the interfering signal for the MPSet is composed as follows:
PathSet(MPSet)＝{traj(mpi),mpi∈MPSet}。
modeling three-dimensional rays in a possible propagation path set PathSet (MPSet) of the interference signal, and constructing a grid-path index to support subsequent screening of the position of the interference source based on path analysis;
the trellis-path index is represented as: i isgt＝{Gr|trInf1,trInf2,...,trInfNt}；
IgtRecording individual grids G within an interferer search region RrPaths/traces traj ═ r that cover all possible interfering signals propagation0,r1,...,rRDmaxTrInf track information for > ∈ PathSet (MPSet)i；
trInfiThe method comprises the following steps: 1) path identifier trajiID ═ i; 2) sp r of the path start point0.sp＝mpm(ii) a 3) Interference signal slave mpmDeparture arrival grid GrPropagation distance dist (G)r|traji) And propagation loss att (G)r|traji)。
For a certain propagation path traj, traversing the ray r in sequence0,r1,...,rRDmaxExperienced individual grids GrAnalyzing the recording track information trInfiGenerating a raster-path index IgtIndex record of < Gr|trInfi＞；
As shown in fig. 4, the specific steps of ray traversal are as follows:
for propagation path trajiThe k-th ray r in (1)kFirst, according to the starting point rkSp and end point rkEp coordinates, calculating ray rkLength r ofkLen; according to the starting point rkSp coordinate O (x)0,y0,z0) Completion ray initiation point rkSp to grid mapping, resulting in the starting point rkCoordinate G (gx) of the grid on which sp is located0,gy0,gz0) Using the grid as ray rkThe current latest traversed raster;
then, let ray rkIs G (gx, gy, gz), respectively calculated from the starting point rkSp along the ray rkAdvancing to the next x-plane, y-plane, z-plane of the lattice passgridA distance l ofx，ly，lz(ii) a According to lx，lyAnd lzCalculating the next grid closest to the traversed position of the current ray; judging that the minimum value is lx，lyOr lz；
When the minimum value is lxTime, ray rkThe relative starting point run length is lxThe next grid coordinate reached is G (gx + dir)xGy, gz), calculate rkInformation on arrival at the grid < trinfi,Gr＞,r＝μ(gx+dirxGy, gz), generate index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update lx:＝lx+dxI.e. the distance from the ray origin r the next time it reaches the x-planekA distance sp; update PassdGrid: ═ G (gx + dir)x,gy,gz)；
When the minimum value is lyTime, ray rkThe relative starting point run length is lyThe next grid coordinate reached is G (gx, gy + dir)yGz), calculate rkInformation on arrival at the grid < trinfi,Gr＞,r＝μ(gx,gy+diryGz), generate index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update ly:＝ly+dyI.e. the distance from the ray origin r the next time y-plane is reachedkA distance sp; update PassdGrid:G (gx, gy + dir)y,gz)；
When the minimum value is lzTime, ray rkThe relative starting point run length is lzIs reached with the next grid coordinate G (gx, gy, gz + dir)z) Calculating rkInformation on arrival at the grid < trinfi,Gr＞,r＝μ(gx,gy,gz+dirz) Generating an index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update lz:＝lz+dzI.e. the distance from the ray origin r the next time it reaches zkA distance sp; update PassdGrid:G (gx, gy, gz + dir)z)；
Wherein d isxIs a ray rkLength of line segment sectioned by adjacent x-planes; dyIs a ray rkThe length of the line segment sectioned by the adjacent y-plane; dzIs a ray rkThe length of the line segment sectioned by the adjacent z-plane; dirx，diryAnd dirzAre respectively rays rkUnit vector uI in x, y and z directionsx，uIy，uIzLength of (d). When unit vector uIxWhen aligned with the x-axis, i.e. unit vector uIxIn the positive direction, dirxThe value is +1, otherwise the unit vector uIxIn the negative direction, dirxA value of-1; define dir in the same wayyAnd dirzThe value of (1).
Finally, when lx，ly，lzAre all longer than the total length r of the ray segmentkLen, then the ray traversal has been completed, terminating the computation.
Step seven, indexing I by analyzing the grid-pathgtAnd screening out grids serving as candidate positions of the interference source.
The method for positioning the interference source based on the signal propagation path provided by the invention has the core that:
simulating an interference signal propagation process through back ray tracking to generate a possible interference signal propagation path; building a grid-path index I by three-dimensional ray modelinggt(ii) a Then analysis IgtAnd finding out grids passed by an interference signal propagation path in the interference source searching region R, screening grids serving as interference source candidate positions from the grids possibly serving as positions of the interference source, calculating a plurality of target values of each candidate position grid, selecting a candidate position grid with an optimal target value, and finishing the grid positioning of the interference source.
The method specifically comprises the following steps:
firstly, the grids covering the propagation path of the interference source constitute a grid set OriGSet, each grid in the set OriGSet is likely to be the position of the interference source, but the probability is different, and the probability of becoming the interference source needs to be calculated through subsequent screening and evaluation operations. Considering that the set OriGSet is large in scale, screening is performed by adopting three filtering modes to obtain a grid set CanGSet of the interference source candidate position:
1) filtering and screening out backward ray tracing starting point mpi∈ MPSet.
Because the grids around the back tracking start point are dense, but most rays of the grids come from the same back tracking start point, the probability of containing the position of an interference point is extremely small;
2) filtering out signal starting point emission point mp from which covered rays comeiThe number of the grids is less than a set threshold From _ num;
3) and filtering out grids with covered strong signal ray ratios smaller than a set threshold S _ rate.
Specifically, the set OriGSet is preliminarily filtered through the following iterative filtering process to generate the interference source candidate position grid set CanGSet with the expected size, as shown in fig. 5, and the filtering algorithm steps are as follows:
initializing parameters of an iterative filtering process;
the parameters comprise the number of starting points From _ num From which the rays covered by the initial grid at least come, the ratio of strong signal paths to the minimum value S _ rate in the propagation paths covered by the grid, the size interval [ NGmin, NGmax ] of a target candidate grid set, the current iteration frequency iteraCount and the total iteration frequency count;
the second step updates the set OriGSet: and (3) screening out the peripheral grids of the back tracking starting point, and initializing: CanGSet ═ OriGSet.
The third step is to update CanGSet: deleting grids with the number of emission points From which the covered rays come being less than From _ num;
step four, updating CanGSet: deleting the grids with the strong signal ray ratio smaller than the S _ rate;
fifthly, judging whether the total number of grids in the candidate grid set CanGSet is in an interval [ NGmin, NGmax ] or whether the current iteration number iterCount exceeds the total iteration number count, if so, terminating iteration and finishing filtering; otherwise, entering the sixth step;
judging whether | CanGSet | < NGmin, if yes, indicating that the filtering condition is too strict, and properly reducing From _ num and S _ rate values; otherwise, CanGSet > NGmax, which indicates that the filtering strength needs to be increased, and the From _ num and S _ rate values are properly increased; and returning to the second step.
Finally, a candidate grid set CanGSet ═ G, which is much smaller than OriGSet in size, is obtainedrAnd the calculation time for evaluating and positioning the interference source by adopting multi-target constraint can be effectively reduced.
Step eight, calculating a plurality of target values of each candidate position grid in a candidate grid set CanGSet, and selecting a grid corresponding to the optimal target value as an interference source positioning grid;
first, a grid G satisfying the following three constraint targets is selected from an interference source candidate position grid set CanGSetrAs candidate interference source positions:
1) grid GrThe strong interference signal propagation paths from the set of backward ray tracing starting points MPSet are received as many as possible. In particular, grid G is traversedrAll ray paths covered, satisfying str (G)r|traji)＞thldstrOf the interference signal propagation path traji∈ PathSet (MPSet) as many as possible;
2) grid GrThe received total interference difference from any 2 propagation paths is as small as possible, i.e. G is assumedrCovers two strong interference signal propagation paths traji,trajj∈ PathSet (MPSet), requires | str (G)r|traji)-str(Gr|trajj) | is as small as possible;
3) grid GrThe covered strong interference signal propagation paths come from as many interference measurement points mp as possible as the starting points of the back ray tracingi∈MPSet；
Then, respectively calculating a plurality of target values of each candidate position grid, and selecting a grid corresponding to the optimal target value as an interference source positioning grid;
the method converts multiple targets into a single target and solves the multi-target optimization problem. Firstly, the optimal solution of each sub-target is obtained by a single-target optimization method, then a plurality of sub-targets are normalized, the distance relative value of each sub-target function from each optimal solution is calculated, and a single target function is obtained after linear weighting combination.
The specific process is as follows:
first step combining trellis-path index IgtComputing each grid G of the candidate grids CanGSetrNumber of covered strong signal ray paths GrPathnum, integrated interference signal strength difference GrPathloss, number of traceback starting points G from which the grid covered rays comer.fromnum；
The integrated interference signal difference of the grid means: the grid receives signals sent from a plurality of different route measuring points, and the sum of the grid transmitting power difference values obtained by deducing the signals sent from two route measuring points from a signal back propagation path is counted. Based on the integrated interference signal difference definition, the pseudo code for calculating the integrated interference signal difference of each grid point in CanGSet can be as follows, and can pass through G after the calculation is finishedrPathloss access grid GrThe integrated interference signal difference value:
grid GrThe method for calculating the number of the transmitting starting points of the covered strong interference signal propagation path comprises the following steps: traverse GrThe number of the transmission starting points of all the covered strong interference signal propagation paths is counted, and the strong interference signal propagation paths and the transmission starting points can be calculated simultaneously. After the calculation is completed, can pass through GrFromnum access grid GrNumber of strong interfering signal paths covered.
The second step respectively calculates respective optimal solutions according to three objective functions, namely the maximum value Path of the number of covered strong signal ray pathsmaxMinimum value of integrated signal strength difference P L ossminAnd the maximum From of the number of covering rays From the back-tracking originmax。
Thirdly, solving the median Path of each sub-target according to the scales of the three target valuesmid，PLossmid，FrommidTo reduceAnd modifying each sub-target function into the following functions according to the influence of the dimension on the weighting result:
wherein f is1(Gr)，f2(Gr)，f3(Gr) Are respectively grids GrThe strong signal ray constraint value calculation function, the comprehensive signal interference difference value constraint value calculation function and the signal source constraint value calculation function.
Fourthly, setting a weighting coefficient w according to the importance degree of each sub-targetiAnd converting the multi-objective optimization problem into a single-objective optimization problem, finishing the evaluation and positioning of the interference source, and finally obtaining the single-objective optimization problem by the following formula:
in CanGSet, an evaluation value of each grid is calculated through a single objective function, and the grid with the minimum evaluation value is a positioning result.
As shown in fig. 6, mpi(i ═ 1,2,3,4) is the selected back ray tracing starting point, bi(i-1, 2,3,4) is a building, and the dotted line is the point mp from which the back ray trace startsiEmitted backward tracing ray generated by simulating interference signal propagation and forming interference signal propagation path, diamond I1As the actual source of interference, the square loci(i ═ 1,2,3) are the candidate grid positions of the interference source obtained after screening, and constitute the set CanGSet ═ loc of the candidate positions of the interference source1,loc2,loc3}。
Passing through 3 interferer candidate locations l in the graphociThe rays are all denser, but the signal strength and path loss information of the rays covered by the 3 grids are different. Positioning a target interference source grid loc from a candidate grid set CanGSet by analyzing ray information covered by the grid and using a multi-objective optimization methodaimBy locaimAs the interference source location. Due to loc1The multi-objective constraint is best satisfied among the three candidate grids, so the location result grid is loc1Center point distance I1The closer the distance, the more ideal the positioning result.
The invention divides the interference source searching area into L sizew×Lw×LhThe three-dimensional grid converts the problem of positioning the interference source into the problem of positioning the grid where the interference source is located, and all subsequent related interference positioning is completed based on the grid. Ray traj generated when the back ray tracesi∈ PathSet passes through 1 grid, the grid is said to cover ray traji。
In the process of positioning the off-network interference source based on the path analysis, if the core idea is not changed, the method selects the drive test points to generate the ray path, screens the grid based on multi-target constraint and the like, but completes the three-dimensional modeling of the ray path by adopting other modes, and still belongs to the protection scope of the patent.
The key points of the technology of the invention are as follows:
(1) off-network interference source positioning based on propagation path analysis
And simulating a reverse ray tracking model, selecting a drive test point set with strong signals, wide distribution and less shielding, carrying out ray tracking to generate a ray path, and carrying out three-dimensional modeling on the ray path to analyze and process the obtained grid-ray record so as to position the grid where the interference source is located.
(2) Positioning method based on multi-target constraint
In the positioning of the off-network interference source based on the path analysis, multi-target constraint is provided, namely conditions which should be met by the interference source are reversely deduced according to the signal propagation condition of a real scene, and the grid position of the candidate interference source is evaluated according to the multi-target constraint.
Claims (5)
1. A method for positioning an external interference source of a mobile communication network is characterized by comprising the following steps:
step one, aiming at an area R covered by a mobile communication network0According to the interference information counted by the base station, the drive test frequency sweep analysis result and the user complaint information, confirming an interference source search area R;
step two, carrying out three-dimensional space rasterization on the interference source searching region R, and recording the spatial distribution of terrain and ground objects in the interference source searching region R by adopting a uniform grid acceleration structure to realize scene modeling;
thirdly, reasonably planning a drive test track and manually implementing drive test according to the measurement radius of the drive test sweep frequency measurement terminal and the distribution condition of interfered points in the region R to obtain an initial measurement result set DTset;
an initial measurement result set DTSet ═ { mp };
the information recorded by each measuring point mp comprises measuring time mp.time, measuring point longitude mp.lat, measuring point latitude mp.lon, interference signal strength mp.str and frequency band mp.freq received at the measuring point;
analyzing the initial measurement result set DTSet, selecting N measurement points from the initial measurement result set DTSet, calculating the signal angle of each measurement point as the starting point of the back ray tracking, and forming a back tracking starting point set MPSet;
step five, sequentially selecting each measuring point from the back tracking starting point set MPSet, and reversely sending rays according to the signal coverage angle range of the measuring point to generate a possible interference signal propagation reverse ray path PathSet (MPSet);
the specific single signal back ray tracing steps are as follows:
step 501, based on the current emission starting point mpiEmitting a primary ray r0And adding the information nodeinfo of the transmission starting point into a data structure ray L ist0；
Step 502, judging the currently tracked ray rkIf the current ray tracing is finished, the end point is the ray and boundary surfaceIntersection, recording end point information nodeinfoeAnd added to ray L ist to obtain a complete slave mpiThe signal propagation path starting at the end point and ending at the end point, and the flow proceeds to step 508; otherwise, go to step 503;
step 503, judging the currently tracked ray rkIf the reflection and diffraction times in the previous propagation process exceed the specified threshold value RDmax, the current ray tracking process is ended and the process goes to step 508; otherwise, go to step 504;
step 504, ray r is performed by using a uniform grid acceleration structurekFast collision detection with the topographic asperities in the interference source search region R;
step 505, judge ray rkWith the current uniform grid UG in the forward directionuWhether a collision is detected by the intersection of the inner building bid and the triangular identifier tid; if so, go to step 506, otherwise, the current ray rkContinue to propagate to the next uniform grid in the forward direction, go back to step 502;
the identification tid for the TIN triangle satisfies the following condition: the projected triangle of the identification tid on the two-dimensional ground plane contains a uniform grid UGuOr identifying the sides of the projected triangle of tid on the two-dimensional ground plane and the uniform grid UGuIntersecting;
step 506, find ray rkThe collision point with the terrain bump is taken as a ray rkEnd point r ofkEp, r iskEp information nodeinfokAdding the data into a data structure ray L ist, judging whether the collision point is on a terrain TIN surface, if so, terminating the current ray tracing, and entering step 508, otherwise, entering step 507 when the collision point is on the building surface;
step 507, judging whether the reflection or diffraction is performed according to the position of the intersection point, and generating a corresponding secondary reflection line or a corresponding diffraction line; returning to step 502 to judge whether the new reflection ray or the new winding ray exceeds the area range;
the method specifically comprises the following steps: if the intersecting collision point is the side surface or the top surface of the building, the next-stage reflection line r is generatedk+1Setting a ray rk+1Is of the type reflection, ray origin rk+1Sp is rkEp and according to the primary ray rkThe direction of (d) is calculated to obtain a reflection ray rk+1Direction along which the collision-generated reflected ray r continues to be tracedk+1Returning to step 502 to judge whether the area is out of range again; otherwise, the position of the crossed collision point is positioned on the edge of the building to generate the next-level ray rk+1Setting a new section of ray rk+1Is a diffraction, ray origin rk+1Sp is rkEp and according to the primary ray rkThe direction of (a) is calculated around the ray rk+1Direction along which the ray-around r generated by the collision is continuously tracedk+1Returning to step 502 to judge whether the area is out of range again;
step 508, aiming at nodeinfo information in the data structure ray L ist, constructing mp from the starting pointiStarting from a plurality of rays r0,r1,...,rRDmaxForming an interference signal propagation path, and calculating the propagation parameters of the path;
the propagation parameters include reflection/diffraction loss coefficient, each segment rkThe number of ray orders, ray rkK is more than or equal to 1 and is less than or equal to RDmax;
one propagation path is denoted as traj ═ r0,r1,...,rRDmax＞；
Tracking all interference measurement points mp in the set of origin MPSets for the back-raysi∈ MPSet, a propagation path traj (mp) sequentially constructed from a data structure ray L isti) The set of possible propagation traces constituting the interfering signal is:
PathSet(MPSet)＝{traj(mpi),mpi∈MPSet}；
step six, modeling the three-dimensional rays in the possible propagation path set PathSet (MPSet) of the interference signal, and constructing a grid-path index Igt；
The trellis-path index is represented as: i isgt＝{Gr|trInf1,trInf2,...,trInfNt}；
IgtRecording individual grids G within an interferer search region RrTrajectory information of all paths/trajectories covered over which interference signals may propagate; trInfiTrajectory information indicating a propagation path traj;
for a certain propagation path traj, traversing the ray r in sequence0,r1,...,rRDmaxExperienced individual grids GrAnalyzing the recording track information trInfiGenerating a raster-path index IgtIndex record of<Gr|trInfi>；
Step seven, indexing I by analyzing the grid-pathgtAnd screening out a grid set CanGSet serving as a candidate position of the interference source, calculating a plurality of target values of grids of each candidate position, and selecting a grid corresponding to the optimal target value as an interference source positioning grid.
2. The method as claimed in claim 1, wherein the step one comprises the following steps:
firstly, according to the base station interference statistical information provided by the network management system, evaluating the interfered degree of the interfered cell/base station, and screening the interference degree exceeding the base station interference threshold value thldeintThe cell/base station of (1) is used as a first type interfered point;
then, the sweep frequency drive test data obtained in advance are analyzed, and interference signal intensity exceeding the sweep frequency drive test interference threshold value thld is screened outdtThe sweep frequency drive test measuring point is used as a second type interfered point;
continuously analyzing the position of the complaint point and the interference severity degree according to the complaint data of the user, thereby determining the position point which is seriously interfered in the network as a third type interfered point;
finally, constructing a rectangular bounding box R according to the positions of the three types of interference points1Will enclose the box R1And expanding the set distance He to the periphery to obtain a new rectangle R, namely the interference source searching area.
3. The method for locating the interference source outside the mobile communication network as claimed in claim 1, wherein in the sixth step, the specific steps of ray traversal are as follows:
for propagation path trajiThe k-th ray r in (1)kFirst, according to the starting point rkSp and end point rkEp coordinates, calculating ray rkLength r ofkLen; according to the starting point rkSp coordinate O (x)0,y0,z0) Completing the ray starting point rkSp to grid mapping, resulting in the starting point rkCoordinates of the grid on which sp is located
Using the grid as ray rkThe current latest traversed raster;
then, let ray rkIs G (gx, gy, gz), respectively calculated from the starting point rkSp along the ray rkDistance l from the advancing direction to the next x-plane, y-plane, z-plane of the lattice passgridx，ly，lz(ii) a According to lx，lyAnd lzCalculating the next grid closest to the traversed position of the current ray; judging that the minimum value is lx，lyOr lz；
When the minimum value is lxTime, ray rkThe relative starting point run length is lxThe next grid coordinate reached is G (gx + dir)xGy, gz), calculate rkInformation on arrival at the grid < trInfi,Gr＞,r＝μ(gx+dirxGy, gz), generate index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update lx:＝lx+dxI.e. the distance from the ray origin r the next time it reaches the x-planekA distance sp; update PassdGrid: ═ G (gx + dir)x,gy,gz)；
When the minimum value is lyTime, ray rkThe relative starting point run length is lyThe next grid coordinate reached is G (gx, gy + dir)yGz), calculate rkInformation on arrival at the grid < trInfi,Gr＞,r＝μ(gx,gy+diryGz), generate index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update ly:＝ly+dyI.e. the distance from the ray origin r the next time y-plane is reachedkA distance sp; update PassdGrid:G (gx, gy + dir)y,gz)；
When the minimum value is lzTime, ray rkThe relative starting point run length is lzIs reached with the next grid coordinate G (gx, gy, gz + dir)z) Calculating rkInformation on arrival at the grid < trInfi,Gr＞,r＝μ(gx,gy,gz+dirz) Generating an index entry < Gr|trInfi> and add the index entry to index IgtPerforming the following steps; update lz:＝lz+dzI.e. the distance from the ray origin r the next time it reaches zkA distance sp; update PassdGrid:G (gx, gy, gz + dir)z)；
Wherein d isxIs a ray rkLength of line segment sectioned by adjacent x-planes; dyIs a ray rkThe length of the line segment sectioned by the adjacent y-plane; dzIs a ray rkThe length of the line segment sectioned by the adjacent z-plane; dirx，diryAnd dirzAre respectively rays rkUnit vector uI in x, y and z directionsx，uIy，uIzLength of (d); when unit vector uIxWhen aligned with the x-axis, i.e. unit vector uIxIn the positive direction, dirxThe value is +1, otherwise the unit vector uIxIn the negative direction, dirxA value of-1; define dir in the same wayyAnd dirzR ∈ N is a one-dimensional number of the grid, N is a natural number set, and the conversion relation between the one-dimensional coordinate r and the three-dimensional coordinate (gx, gy, gz) is rho: r → (gx, gy, gz), μ (gx, gy, gz) → r;
finally, when lx，ly，lzAre all longer than the total length r of the ray segmentkLen, then the ray traversal has been completed, terminating the computation.
4. The method as claimed in claim 1, wherein the screening algorithm in step seven is as follows:
firstly, initializing parameters of an iterative filtering process;
the parameters comprise the number of starting points From _ num From which the rays covered by the initial grid at least come, the ratio of strong signal paths to the minimum value S _ rate in the propagation paths covered by the grid, the size interval [ NGmin, NGmax ] of a target candidate grid set, the current iteration frequency iteraCount and the total iteration frequency count;
and secondly, updating a grid set OriGSet covering the propagation path of the interference source: and (3) screening out the peripheral grids of the back tracking starting point, and initializing: CanGSet ═ OriGSet;
and step three, updating the CanGSet: deleting grids with the number of emission points From which the covered rays come being less than From _ num;
step four, updating CanGSet: deleting the grids with the strong signal ray ratio smaller than the S _ rate;
fifthly, judging whether the total number of grids in the candidate grid set CanGSet is in an interval [ NGmin, NGmax ] or whether the current iteration number iterCount exceeds the total iteration number count, if so, terminating iteration and finishing filtering; otherwise, entering the sixth step;
sixthly, judging whether the CanGSet is less than NGmin or not, if so, indicating that the filtering condition is too strict, and properly reducing the From _ num and S _ rate values; otherwise, CanGSet > NGmax, which indicates that the filtering strength needs to be increased, and the From _ num and S _ rate values are properly increased; and returning to the second step.
5. The method according to claim 1, wherein the step seven of selecting the grid corresponding to the optimal target value as the grid for positioning the interference source specifically comprises:
first step, combine trellis-path index IgtComputing each grid G of the candidate grids CanGSetrNumber of covered strong signal ray paths GrPathnum, integrated interference signal strength difference GrPathloss, grid-covered rayNumber of back tracking start points from Gr.fromnum；
Second, respectively solving the number G of covered strong signal ray pathsrMaximum Path of pathnummaxThe difference G of the integrated interference signal intensityrMinimum value of pathloss P L ossminAnd the number of covering rays from the back-tracking starting point GrMaximum From of Frommax；
Thirdly, solving the median Path of each sub-target according to the scales of the three target valuesmid，PLossmid，FrommidTo reduce the influence of the dimension on the weighting result, each sub-targeting function is modified as follows:
wherein f is1(Gr) Is a grid GrThe strong signal ray constraint value of (f) is calculated2(Gr) Is a grid GrThe integrated signal interference difference value constraint value of f3(Gr) Is a grid GrA signal source constraint value calculation function of;
fourthly, setting a weighting coefficient w according to the importance degree of each sub-targetiAnd converting the multi-objective optimization problem into a single-objective optimization problem, finishing the evaluation and positioning of the interference source, and finally obtaining the single-objective optimization problem by the following formula:
and fifthly, calculating an evaluation value of each grid through a single objective function in the grid set CanGSet, wherein the grid with the minimum evaluation value is the interference source positioning grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910909518.XA CN110677865B (en) | 2019-09-25 | 2019-09-25 | Method for positioning external interference source of mobile communication network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910909518.XA CN110677865B (en) | 2019-09-25 | 2019-09-25 | Method for positioning external interference source of mobile communication network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110677865A CN110677865A (en) | 2020-01-10 |
| CN110677865B true CN110677865B (en) | 2020-08-04 |
Family
ID=69079186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910909518.XA Active CN110677865B (en) | 2019-09-25 | 2019-09-25 | Method for positioning external interference source of mobile communication network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110677865B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111257653B (en) * | 2020-01-16 | 2021-06-22 | 湘潭大学 | Electromagnetic radiation prediction method under underground pedestrian passageway scene |
| CN116056204A (en) * | 2021-10-28 | 2023-05-02 | 中兴通讯股份有限公司 | Interference source positioning method, device and computer readable medium |
| CN114286375B (en) * | 2021-12-16 | 2023-08-18 | 北京邮电大学 | Mobile communication network interference positioning method |
| CN114629761B (en) * | 2022-03-22 | 2023-10-27 | 吉林省广播电视研究所(吉林省广播电视局科技信息中心) | Frequency modulation signal inverse triangular frequency measurement demodulation method |
| CN117761678B (en) * | 2024-02-22 | 2024-04-26 | 成都鹰谷米特科技有限公司 | Complex environment target detection method and chip based on V frequency band |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102694607A (en) * | 2011-03-25 | 2012-09-26 | 北京海曼无限信息技术有限公司 | Networked software radio signal monitoring method aiming at mobile interference source and system thereof |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102083103B (en) * | 2009-11-27 | 2012-12-05 | 中国移动通信集团北京有限公司 | Method and apparatus for determining positioning reference point in interference source positioning system |
| CN103796219B (en) * | 2012-10-26 | 2017-02-08 | 同济大学 | Long distance cofrequency interference source detection and positioning method for TD-LTE system |
| CN107144812B (en) * | 2017-06-30 | 2020-04-17 | 北京德辰科技股份有限公司 | Single-station mobile automatic positioning method |
| CN108391279A (en) * | 2018-02-08 | 2018-08-10 | 南瑞集团有限公司 | A kind of outer interference source failure method for tracking and positioning of LTE electric power wireless private network net |
-
2019
- 2019-09-25 CN CN201910909518.XA patent/CN110677865B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102694607A (en) * | 2011-03-25 | 2012-09-26 | 北京海曼无限信息技术有限公司 | Networked software radio signal monitoring method aiming at mobile interference source and system thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110677865A (en) | 2020-01-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110677865B (en) | Method for positioning external interference source of mobile communication network | |
| CN109963287B (en) | Antenna direction angle optimization method, device, equipment and medium | |
| Rautiainen et al. | Verifying path loss and delay spread predictions of a 3D ray tracing propagation model in urban environment | |
| CN109188355B (en) | Method for optimizing and optimally distributing receiving antennas of multi-point positioning system | |
| CN102546039B (en) | Radio wave propagation prediction method and device | |
| Wölfle et al. | A fast and enhanced ray optical propagation model for indoor and urban scenarios, based on an intelligent preprocessing of the database | |
| CN105005039B (en) | Satellite-signal localization method and system based on 3D modeling scene dynamics fingerprint | |
| US20070093212A1 (en) | Radio wave propagation characteristic estimation system, and its method and program | |
| CN110907901B (en) | Sea warfare passive interference calculation and evaluation method based on TDSBR | |
| CN104093205A (en) | Method for deploying anchor nodes of wireless positioning system based on received signal strength indication | |
| Schmitz et al. | Efficient rasterization for outdoor radio wave propagation | |
| CN106921989B (en) | A kind of communication network field strength distribution determines method and device | |
| US20190320326A1 (en) | Method and system for radio communication network planning | |
| CN115205454A (en) | Space segmentation acceleration method and system suitable for ray tracing wireless channel modeling | |
| JP3823917B2 (en) | Ray reception determination method and system, and radio wave propagation characteristic estimation method using the same | |
| CN114286375B (en) | Mobile communication network interference positioning method | |
| CN106528956B (en) | Data interpolations field intensity prediction method based on ray tracing models | |
| CN113301503B (en) | Weak coverage area detection method and device | |
| KR100927962B1 (en) | Ray Tracing Device and Method Using Path Preprocessing | |
| CN114584992A (en) | Measurement and control station alternative station address acquisition method and measurement and control station layout planning method | |
| Schmitz et al. | Simulation of radio wave propagation by beam tracing. | |
| Schmitz et al. | Beam tracing for multipath propagation in urban environments | |
| KR20010011922A (en) | Prediction and analysis method of cell planning for wireless communication | |
| CN114363909A (en) | Azimuth angle determining method and device, electronic equipment and storage medium | |
| CN114900841B (en) | Site selection method of communication site and communication device |
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 |