CN103067931A - Method and device capable of positioning interference source in wireless network - Google Patents
Method and device capable of positioning interference source in wireless network Download PDFInfo
- Publication number
- CN103067931A CN103067931A CN201110321864XA CN201110321864A CN103067931A CN 103067931 A CN103067931 A CN 103067931A CN 201110321864X A CN201110321864X A CN 201110321864XA CN 201110321864 A CN201110321864 A CN 201110321864A CN 103067931 A CN103067931 A CN 103067931A
- Authority
- CN
- China
- Prior art keywords
- grid
- interference
- drive test
- base station
- residential quarter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
The invention provides a method and device capable of positioning an interference source in a wireless network. The method includes the steps of gaining a plurality of communities which are most interfered, defining an interference area according to the plurality of communities, dividing the interference area into a plurality of lattice cells, respectively adopting a ray tracing propagation model on each community of the plurality of communities, confirming gathering of the lattice cells which produce interference on the community in the interference area, getting intersection of the gathering of the lattice cells corresponding to each community, and gaining the gathering of target lattice cell where the interference source is placed. The method and device capable of positioning the interference source is capable of improving positioning accuracy of the interference source.
Description
Technical field
The present invention relates to wireless communication technology field, particularly method and the device of interference source in a kind of positioned radio network.
Background technology
At present, following two kinds of methods are mainly adopted in the wireless network interference analysis:
Method one
Utilize the network measure data, which residential quarter exists quality problems, possible up-downgoing interference and interference type in the phase-split network.
This method generally can only be positioned at interference source the residential quarter one-level, can't further analyze the disturbed condition of the inner zones of different in residential quarter.
Even in conjunction with mobile positioning technique, for example determine mobile terminal locations according to measurement report, and then infer the disturbed condition at diverse location place in the network that because the restriction of the precision of location technology itself, the precision of interference analysis also can only reach about 200m * 200m.
Method two
According to network configuration data, adopt wireless signal propagation model, the environmental geography information that combining geographic information system (GIS, Geographic Information System) provides, by simulation means, interference and the signal distributions situation of diverse location in calculating, the phase-split network; Simultaneously, in order to improve the precision of propagation model, adopt drive test data that propagation model is proofreaied and correct.
In the method, geography information generally is two-dimentional but not three-dimensional, and propagation model then adopts traditional empirical model, such as OH (Okumara-Hata) model, cost231 model etc.Under the dense city environment, building blocks and causes wireless signal reflection, refraction, diffraction, cause all kinds of non line of sight (NLOS, Non Line of Sight) propagate, though the experience propagation model of proofreading and correct through drive test data also can't the accurate description wireless signal strength and propagation distance between relation.Therefore, the positioning accuracy of utilizing the simulation analysis means of experience propagation model to arrive is also limited, preferably also can only reach about 300m * 300m.
As seen, the positioning accuracy of prior art interference analysis is limited, and particularly for outer disturb (interference unit, rival's base station and the repeater etc.) of net, solve its orientation problem is a difficult problem always.
Summary of the invention
Technical problem to be solved by this invention provides method and the device of interference source in a kind of positioned radio network, to improve the positioning accuracy to interference source.
For solving the problems of the technologies described above, it is as follows to the invention provides technical scheme:
The method of interference source in a kind of positioned radio network comprises:
Obtain the strongest a plurality of residential quarters that are interfered, determine an interference region according to described a plurality of residential quarters;
Described interference region is divided into a plurality of grids;
For each residential quarter in described a plurality of residential quarters, adopt respectively the ray trace propagation model, determine in the described interference region this residential quarter to be caused the grid set of interference;
Common factor is got in the grid set that each residential quarter is corresponding, obtained the target grid set at interference source place.
Above-mentioned method, wherein, described employing ray trace propagation model determines in the described interference region this residential quarter to be caused the grid set of interference, specifically comprises:
Adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, and received power is greater than the number of the grid of default power threshold in the described interference region;
Adjust the transmitting power of described base station, until described number when equaling preset value, obtains described grid set.
Above-mentioned method, wherein, the transmitting power of the described base station of described adjustment is:
When described number is 0, increase step by step the transmitting power of described base station;
When described number during greater than described preset value, reduce step by step the transmitting power of described base station.
Above-mentioned method wherein, also comprises:
When the grid number in the set of described target grid greater than 1 the time, utilize drive test frequency sweep data to determine the target grid at interference source place.
Above-mentioned method, wherein, the described target grid that utilizes drive test frequency sweep data to determine the interference source place specifically comprises:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively ray tracing algorithm, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
The device of interference source in a kind of positioned radio network comprises:
The interference region determination module is used for obtaining the strongest a plurality of residential quarters that are interfered, and determines an interference region according to described a plurality of residential quarters;
The grid division module is used for described interference region is divided into a plurality of grids;
The first locating module is used for each residential quarter for described a plurality of residential quarters, adopts respectively the ray trace propagation model, determines in the described interference region this residential quarter to be caused the grid set of interference;
The module that seeks common ground is used for common factor is got in grid set corresponding to each residential quarter, obtains the target grid set at interference source place.
Above-mentioned device, wherein, described the first locating module specifically is used for:
Adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, and received power is greater than the number of the grid of default power threshold in the described interference region;
Adjust the transmitting power of described base station, until described number when equaling preset value, obtains described grid set.
Above-mentioned device, wherein, described the first locating module specifically is used for:
When adjusting the transmitting power of described base station, when described number is 0, increase step by step the transmitting power of described base station; When described number during greater than described preset value, reduce step by step the transmitting power of described base station.
Above-mentioned device wherein, also comprises:
The second locating module is used for the grid number gathered when described target grid greater than 1 the time, utilizes drive test frequency sweep data to determine the target grid at interference source place.
Above-mentioned device, wherein, described the second locating module specifically is used for:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively the ray trace propagation model, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
Compared with prior art, the present invention has adopted the ray trace propagation model to come the interference source in the wireless network is positioned, because the ray trace propagation model has been considered the wireless signal non-line-of-sight propagation, computational accuracy is high, is applicable to the complex environments such as dense city.
Description of drawings
Fig. 1 is the method flow diagram according to interference source in the positioned radio network of the embodiment of the invention;
Fig. 2 is the schematic diagram of grid division in the embodiment of the invention;
Fig. 3 is the structure drawing of device according to interference source in the positioned radio network of the embodiment of the invention.
Embodiment
Under the complex dense urban area circumstance, wireless signal has four kinds of circulation ways: direct projection, reflection, scattering and diffraction.The rule-of-thumb approach such as OH, COST231 radio transmission model because can't accurate analysis, calculate the non-line-of-sight propagation features such as reflection, diffraction of signal, computational accuracy is not high, predicated error is larger.Specifically, generally about 10dB, its positioning accuracy can only reach about 300m * 300m the received signal strength mean square deviation of this propagation model prediction.
The ray trace propagation model has been considered the wireless signal non-line-of-sight propagation, and computational accuracy is high, is applicable to the complex environments such as dense city.Specifically, the residential quarter received signal strength and the actual test received signal strength mean square deviation that adopt the ray trace model to obtain can accomplish that about 3dB, its positioning accuracy can reach about 30m * 30m.
Based on this, the embodiment of the invention adopts the ray trace propagation model to come the interference source in the wireless network is positioned, and raising that can be larger is to the positioning accuracy of interference source.
For ease of better understanding embodiments of the invention, at first the ray trace propagation model is simply introduced here.
The ray trace propagation model will be regarded ray/beam as from the wireless signal that emission source (such as the base station) gives off based on " locality " of wireless signal radio-frequency field, and ray energy is independently being propagated in the ray tube separately.The signal communication process is reduced to direct projection, reflection and diffraction; Utilize 3D building information in the communication environments that electronic chart provides, all possible signal was broadcast the path between search was followed the trail of from the launch point to the acceptance point, calculate signal propagation losses and acceptance point field intensity on every paths, the field intensity in all paths that same acceptance point place is arrived is carried out vector superposed, obtains the total received field strength in acceptance point place and signal strength signal intensity.
Ray trace is divided into the forward ray trace and return radiation is followed the tracks of.In forward ray trace process, the emission sources such as base station towards periphery space uniform send some rays, from source point, follow the tracks of respectively the propagation path of every ray, judge that at the acceptance point place beam is to the contribution of acceptance point place field intensity.
Return radiation is followed the tracks of and to be set out by acceptance point/field point, and each bar of backward tracing arrives a possible path of point from source point, and very little and very little to the location precision propagation path of field intensity when ignoring those and arriving fields and put.In general negligible through the ray behind twice diffraction.Return radiation is followed the tracks of and at first will be sought visible of source point and show and split table, and table is shown and split to visible an of point, and source point and common visible an of point show and split table, then follows the trail of reflection, the diffraction situation of ray, obtains final result of calculation.
Based on above-mentioned Radio Frequency Tracking propagation model, the embodiment of the invention provides method and the device of interference source in the following positioned radio network
With reference to Fig. 1, the method for interference source can comprise the steps: in the positioned radio network of the embodiment of the invention
Step 101: obtain the strongest a plurality of residential quarters that are interfered, determine an interference region according to described a plurality of residential quarters;
Step 102: described interference region is divided into a plurality of grids;
Step 103: for each residential quarter in described a plurality of residential quarters, adopt respectively the ray trace propagation model, determine in the described interference region this residential quarter to be caused the grid set of interference;
Step 104: common factor is got in the grid set that each residential quarter is corresponding, obtained the target grid set at interference source place;
Step 105: when the grid number in the set of described target grid greater than 1 the time, utilize drive test frequency sweep data to determine the target grid at interference source place.
Wherein, step 105 is optional step.Because after completing steps 104, just obtained the target grid set at interference source place, if the grid number in the set of described target grid equals 1, this grid is the target grid at interference source place so; If the grid number in the described target grid set is greater than 1, generally speaking, its grid number can be not too many, takes so the artificial mode of getting rid of to obtain easily the target grid at interference source place yet.
In step 101, can according to the network measure data, obtain the strongest a plurality of residential quarters that are interfered, then, just can determine an interference region roughly according to described a plurality of residential quarters, for example, with the directional diagram phase of described a plurality of residential quarters with, obtain described interference region.
In addition, can also be according to the network measure data, and in conjunction with conventional wireless propagation models such as OH, COST231, come the described interference region of Primary Location.For example, by adjacent cell support (Neighboring Cell Support, NCS) data, drive test frequency sweep data, in conjunction with measurement report record (Measurement Report Recording, MRR), frequency is distributed support (Frequency Allocation Support, FAS), idle channel is measured (Idle Channel Measurements, ICM) data, which the residential quarter that can tentatively obtain being interfered has, by show the position at these places, residential quarter at map, and in conjunction with traditional propagation model, just can access the interference source of these residential quarters, thereby locate described interference region.
The positioning accuracy of present this method mainly depends on Timing Advance (Timing Advance, TA), TA is per 550 meters yardsticks, 0～550 meter of TA=0 correspondence, 550～1100 meters of TA=1 correspondences, ..., so the positioning accuracy of this method is difficult to improve, the scope of the interference region of Primary Location is roughly about 500m * 500m.
In the present embodiment, uplink interference MRR, FAS, the ICM data of receiving according to 3-5 base station cell (CELL1, CELL2, CELL3, CELL4, CELL5), the extraction the strongest 3-4 residential quarter that is interfered, and with the interference region about this Primary Location 500m * 500m.
Below above-mentioned several measurement data are introduced respectively.
(1)NCS
Can access the disturbed coefficient of residential quarter according to the NCS statistics, mainly be to use its relative signal intensity measurement function to extract co-channel interference factor C/I＜9, the adjacent frequently data of interference factor C/A＜-9, between the residential quarter, form interference matrix, according to the interference coefficient of described interference matrix calculation plot frequency.Interference coefficient can be divided into disturbed coefficient and interference source coefficient.Disturbed coefficient is exactly to be subject to same, the adjacent traffic ratio of frequently disturbing in adjacent area, and unit is percentage; The interference source coefficient be exactly this residential quarter as the adjacent area, to same, the adjacent traffic ratio of frequently disturbing that other adjacent areas produce, unit is percentage.
After getting the above-mentioned NCS measurement result file of making of instruction, with related tool destination file is processed, obtained calculating the NCS statistics that interference factor etc. needs.Cooperate again the frequency situation of main plot, adjacent cell namely can calculate disturbed coefficient and the interference source coefficient of main adjacent cell.
Wherein, the computing formula of the C/I of main plot, C/A is:
The main plot is become by the frequency array of the coefficient of an adjacent area interference by C/I interference factor, C/A interference factor and main plot, and formula is as follows:
The disturbed overall coefficient of main plot is all adjacent area interference coefficient sums, and formula is as follows:
In the above-mentioned formula:
The adjacent area refers to the actual Broadcast Control Channel that can measure (Broadcast Control Channel, BCCH) frequency point cell, contains definition (DEF) and undefined (UNDEF) two parts;
RepArfcn is the measurement report sum of adjacent area, derives from the statistics of above-mentioned NCS;
TimesRelSS and TimesRelSS2 distinguish the relative signal intensity counter of corresponding C/I and C/A, derive from the statistics of above-mentioned NCS;
With again and again count, the neighbour counts again and again, the total frequency number in adjacent area derives from residential quarter definition data (Cell Design Data, CDD) medium frequency data, the time of CDD will be consistent with the time that NCS measures.
According to the disturbed coefficient data after the NCS analysis, get front 5～10 main plot of disturbed overall coefficient rank, main plot, front 2 adjacent cell corresponding to disturbed coefficient coefficient of rank then got in interference adjacent area corresponding to main plot, and with these residential quarters as the needed serious residential quarter that is interfered.
(2)MRR
1. filter out the residential quarter of up average signal strength＞=-85 and uplink speech quality＜=-90.
2. filter out the residential quarter of descending average signal strength＞=-85 and downlink voice quality＜=-90.
Get 1,2 both residential quarters and sort from small to large (the less interference of value is larger) according to the uplink and downlink voice quality, obtain the needed serious residential quarter that is interfered.
(3)FAS
Getting and measuring frequency is that the residential quarter is at frequency and avMedian-110＞=-95 (avMedian: the residential quarter mean value of all interference level samples in measuring period) of usefulness, sort from big to small (the larger interference of value is larger) according to avMedian, obtain the needed serious residential quarter that is interfered.
(4)ICM
The residential quarter of getting (ICM1+ICM2+ICM3)/(ICM1+ICM2+ICM3+ICM4+ICM5)＞=50% screens, and from big to small ordering (it is larger to be worth the outer interference of larger net) obtains the needed serious residential quarter that is interfered.The implication of ICM is as shown in the table:
| The interference fringe grade | The level value scope |
| ICM1 | -110～-100dBm |
| ICM2 | -100～-95dBm |
| ICM3 | -95～-90dBm |
| ICM4 | -90～-85dBm |
| ICM5 | -85～-47dBm |
The above-mentioned result who obtains respectively according to NCS, MRR, FAS, ICM data is carried out combined analysis, can extract the strongest 3-4 residential quarter that is interfered.The merging method for example will appear at the strongest residential quarter that is interfered the most, residential quarter among above-mentioned 4 results simultaneously; Perhaps, the strongest residential quarter that is interfered the most, the more residential quarter of number of times will be appeared among above-mentioned 4 results.
In the present embodiment, suppose that 4 the strongest residential quarters that are interfered that obtain are: CELL1, CELL2, CELL3, CELL4.
In step 102, the interference region of Primary Location is divided into a plurality of grid G (i, j), i is longitude, j is latitude.As previously mentioned, the positioning accuracy of the ray trace propagation model of embodiment of the invention employing is higher than existing traditional experience propagation model.
In addition, in the ray trace propagation model, can also improve tracking accuracy by increasing radiographic density, operand can increase like this.A ray line (θ, Φ) is sent in the base station, and θ and Φ are respectively level angle, the vertical angle of line, and this ray is through direct projection, reflection, diffraction, arrives certain grid G (i, j) in this base station certain limit.By adjusting the density of θ and Φ variable, it is big or small for behind the geographical grid about 30m * 30m one by one to guarantee that interference region is divided into, and one and above signal estimation point are arranged in each grid.
Can be with reference to the grid division to the overlay area of base station to the grid division of interference region, with reference to Fig. 2, specific as follows:
According to covering sideband signal strength threshold value, as-100dbm, determine that the residential quarter covers border, for example radius R
0=1200m; The angle of coverage scope Φ of appointment take the azimuth, residential quarter as symmetry axis
0, such as 90 degree.
According to parameters R
0And Φ
0, the coverage of base station cell is divided into the grid that a plurality of length are L * L, according to the required precision of interference analysis, can get L=30m herein.
Work as R0=1200m, when Φ 0=90 spent, the grid number had approximately:
In step 103, for each residential quarter, adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, received power is greater than the number of the grid of default power threshold in the described interference region, and constantly adjust described transmitting power and calculate, until described number equals preset value.Because the transmitting power of interference source is unknown, can be first with 0.5w or 1w as transmitting power, then adopt a convergence algorithm, guarantee that described number equals preset value, for example 5.Particularly:
When described number is 0, increases step by step the transmitting power of described base station, and calculate received power in the described interference region greater than the number of the grid of default power threshold, until described number stops when being increased to preset value;
When described number during greater than described preset value, reduce step by step the transmitting power of described base station and calculate in the described interference region received power greater than the number of the grid of default power threshold, until described number stops when being reduced to preset value.
For example, when the interference grid number that falls into interference region when reckoning equaled 0, it was P (t)=P (t)+Δ P that transmit signal strength P (t) increases Δ P, Δ P is 1dB, 2dB, 3dB, 4dB...... is until there be 3-5 to disturb grid to fall in the interference region;
The interference grid number that falls into interference region when reckoning is greater than 5 the time, and it is P (t)=P (t)-Δ P that transmit signal strength P (t) reduces Δ P, and Δ P is 1dB, 2dB, and 3dB, to disturb grid to fall into interference region interior until there be 3-5 for 4dB.......
Grid set (namely falling into the interference grid of interference region) corresponding to each residential quarter used as shown in the formula subrepresentation:
CELL1=∑ G (ia, ja) is a=1 wherein, and 2 ..., 5
CELL2=∑ G (ib, jb) is b=1 wherein, and 2 ..., 5
CELL3=∑ G (ic, jc) is c=1 wherein, and 2 ... 5
CELL4=∑ G (id, jd) is d=1 wherein, and 2 ... 5
Below introduce the method according to the received power of ray trace propagation model computation grid that adopts in this step.
Suppose certain base station cell cellID, send a ray line (θ from this base station, Φ), θ and Φ are respectively level angle, the vertical angle of line, and this ray is through direct projection, reflection, diffraction, arrive certain grid G (i in this base station certain limit, j), calculate in accordance with the following steps the signal received power that comes from line (θ, Φ) that G (i, j) receives:
Step S1: from the base station cell table, find the effective isotropic radiated power (Effective Isotropic Radiated Power, EiRP) of cellID:
The transmitting power P of EiRP=base station cellID
T0Gain on the maximum direction of-feeder line/mixer loss+antenna, its unit is dbm.
Step S2: the horizontal azimuth Dir and vertical angle of declination Tilt that from the table of residential quarter, find cellID, according to ray line (θ, horizontal gain figure Φ) and vertical gain figure, calculate the horizontal azimuth Dir and the loss that vertical angle of declination Tilt causes that depart from antenna owing to the direction θ of line and Φ, the transmitting power (dbm of unit) when obtaining line and sending from antenna end:
dbmP
t＝EiRP-HLoss(Dir，θ)-VLoss(Tilt，Φ)
Illustrate:
Each antenna in cell model, antenna horizontal gain figure, vertical gain figure need to provide as configuration data, are known;
To the antenna damping that horizontal direction produces, calculate the θ of ray line and the difference between antenna in cell horizontal azimuth Dir | Dir-θ |, substitution horizontal gain graph model obtains level decay gain HLoss (Dir, θ);
Similarly, according to | Tilt-Φ |, can calculate vertical decay gain VLoss (Tilt, Φ).
Step S3: with the antenna end transmitting power dbmP take dbm as unit
tBe converted to the antenna end transmitting power wP take W as unit
t
dbmP
t＝10lg(wP
t×10
3)
＝30+10lg(wP
t)
lg(wP
t)＝(dbmP
t-30)/10
Obtain wP by following formula
tOccurrence:
Step S4: follow the tracks of the path of ray line (θ, Φ) from the antenna end of base station cellID to grid G (i, j), the received power that comes from line (θ, Φ) that G (i, j) receives is:
Wherein:
Wavelength X, the m of unit;
Computational methods λ=c/f, c=3 * 10
5Km/s=3 * 10
8M/s; F is frequency, such as 900M, and the Hz of unit.
When f=900MHz=900 * 10
6Hz=9 * 10
8During Hz,
Its wavelength X
900Be λ
900=3 * 10
8M/s ÷ (9 * 108Hz)=1/3m=0.333m,
For GSM 1800 networks, when f=1800, λ
1800=λ
900/ 2=0.1665m
∏: circumference ratio, 3.1415926
dis(cellid，G(i，j))：
The path total length of ray line (θ, Φ) from the antenna end of base station cellID to grid G (i, j), unit is m.This path is comprised of mulitpaths such as direct projection, reflection, diffraction
R
u: on the propagation path, the reflection coefficient when the u secondary reflection occurs.Suppose to occur altogether the m secondary reflection, reflection coefficient is relevant with reflection angle.
T
v: on the propagation path, the diffraction coefficient when the v time diffraction occurs.Suppose to occur altogether n time diffraction.Diffraction coefficient is relevant with the diffraction angle.
Suppose total n bar propagation path of (i, j) from cellid to G.By the way, calculate every the propagation path line that comes from cellid that G (i, j) receives
k(θ
k, Φ
k) received power P
r(line
k) (w of unit), with its addition, obtain the total received power P that G (i, j) receives from cellid
r(cellid, G (i, j)):
By following formula, be the P of W with unit
r(cellid, G (i, the j)) unit of being converted to is the received power of dbm:
DbmPr (cellid, G (i, j))=10lg{Pr (cellid, G (i, j)) * 10
3, the dbm of unit.
In step 104, common factor is got in the grid set that each residential quarter is corresponding, obtain the target grid set DIM at interference source place:
DIM＝∑G(ia，ja)∩∑G(ib，jb)∩∑G(ic，jc)∩∑G(id，jd)
When the grid number among the set DIM was 1, this grid was the target grid at interference source place; When the grid number of set among the DIM greater than 1 the time, for example, it comprises that 4 are disturbed grid DIM1, DIM2, DIM3, DIM4 then can continue execution in step 105, further locating which grid is most probable interference source.
In the step 105, the described target grid that utilizes drive test frequency sweep data to determine the interference source place comprises:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively ray tracing algorithm, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
Be exemplified below:
At first, from the frequency sweep drive test route of interference region, select 4 drive test point interference signal intensity { INT1, INT2, INT3, INT4}
Then, utilize ray tracing method, respectively from DIM1, DIM2, DIM3, DIM4 calculates the signal level of 4 drive test points:
DIM1{rxleve1，rxleve2，rxleve3，rxleve4}
DIM2{rxleve1，rxleve2，rxleve3，rxleve4}
DIM3{rxleve1，rxleve2，rxleve3，rxleve4}
DIM4{rxleve1，rxleve2，rxleve3，rxleve4}
At last, from DIM1, DIM2, DIM3, the signal level of 4 drive test points that DIM4 calculates is with actual interference signal strength signal intensity { INT1, the INT2 of frequency sweep drive test route, INT3, INT4} compares, take the frequency sweep data with calculate that the DIM of mean square deviation minimum value between the data is as the target grid at interference source place.
Wherein, the mean square deviation between frequency sweep data and the reckoning data is:
SDIM1=((rxleve1-INT1) ∧ 2+ (rxleve2-INT2) ∧ 2+ (rxleve3-INT3) ∧ 2+ (rxleve4-INT4) ∧ 2)/4) square root.
The computational methods of SDIM2, SDIM3, SDIM4 are similar.
4 some frequency sweep data and calculate mean square deviation sum minimum value=MIN{SDIM1 between the data then, SDIM2, SDIM3, the DIM that SDIM4} is corresponding are the target grid at interference source place.
Corresponding to the method for interference source in the above-mentioned positioned radio network, the embodiment of the invention also provides the device of interference source in a kind of positioned radio network.
With reference to Fig. 3, described device can comprise such as lower module:
Interference region determination module 10 is used for obtaining the strongest a plurality of residential quarters that are interfered, and determines an interference region according to described a plurality of residential quarters;
The first locating module 30 is used for each residential quarter for described a plurality of residential quarters, adopts respectively the ray trace propagation model, determines in the described interference region this residential quarter to be caused the grid set of interference;
The module 40 that seeks common ground is used for common factor is got in grid set corresponding to each residential quarter, obtains the target grid set at interference source place.
Wherein, described the first locating module specifically is used for:
Adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, and received power is greater than the number of the grid of default power threshold in the described interference region;
Adjust the transmitting power of described base station, until described number when equaling preset value, obtains described grid set.
Wherein, the transmitting power of the described base station of described adjustment is:
When described number is 0, increase step by step the transmitting power of described base station;
When described number during greater than described preset value, reduce step by step the transmitting power of described base station.
Preferably, described device also comprises:
The second locating module (not shown) is used for the grid number gathered when described target grid greater than 1 the time, utilizes drive test frequency sweep data to determine the target grid at interference source place.
Wherein, described the second locating module specifically is used for:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively the ray trace propagation model, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
In sum, the embodiment of the invention has adopted the ray trace propagation model to come the interference source in the wireless network is positioned, because the ray trace propagation model has been considered the wireless signal non-line-of-sight propagation, computational accuracy is high, be applicable to the complex environments such as dense city, its positioning accuracy can reach about 30m * 30m.
Should be noted that at last, above embodiment is only unrestricted in order to technical scheme of the present invention to be described, those of ordinary skill in the art is to be understood that, can make amendment or be equal to replacement technical scheme of the present invention, and not breaking away from the spiritual scope of technical solution of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.
Claims (10)
1. the method for interference source in the positioned radio network is characterized in that, comprising:
Obtain the strongest a plurality of residential quarters that are interfered, determine an interference region according to described a plurality of residential quarters;
Described interference region is divided into a plurality of grids;
For each residential quarter in described a plurality of residential quarters, adopt respectively the ray trace propagation model, determine in the described interference region this residential quarter to be caused the grid set of interference;
Common factor is got in the grid set that each residential quarter is corresponding, obtained the target grid set at interference source place.
2. the method for claim 1 is characterized in that, described employing ray trace propagation model determines in the described interference region this residential quarter to be caused the grid set of interference, specifically comprises:
Adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, and received power is greater than the number of the grid of default power threshold in the described interference region;
Adjust the transmitting power of described base station, until described number when equaling preset value, obtains described grid set.
3. method as claimed in claim 2 is characterized in that, the transmitting power of the described base station of described adjustment is:
When described number is 0, increase step by step the transmitting power of described base station;
When described number during greater than described preset value, reduce step by step the transmitting power of described base station.
4. such as claim 1,2 or 3 described methods, it is characterized in that, also comprise:
When the grid number in the set of described target grid greater than 1 the time, utilize drive test frequency sweep data to determine the target grid at interference source place.
5. method as claimed in claim 4 is characterized in that, the described target grid that utilizes drive test frequency sweep data to determine the interference source place specifically comprises:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively ray tracing algorithm, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
6. the device of interference source in the positioned radio network is characterized in that, comprising:
The interference region determination module is used for obtaining the strongest a plurality of residential quarters that are interfered, and determines an interference region according to described a plurality of residential quarters;
The grid division module is used for described interference region is divided into a plurality of grids;
The first locating module is used for each residential quarter for described a plurality of residential quarters, adopts respectively the ray trace propagation model, determines in the described interference region this residential quarter to be caused the grid set of interference;
The module that seeks common ground is used for common factor is got in grid set corresponding to each residential quarter, obtains the target grid set at interference source place.
7. device as claimed in claim 6 is characterized in that, described the first locating module specifically is used for:
Adopt the ray trace propagation model, the base station of calculating this residential quarter is during with default transmitting power emission, and received power is greater than the number of the grid of default power threshold in the described interference region;
Adjust the transmitting power of described base station, until described number when equaling preset value, obtains described grid set.
8. device as claimed in claim 7 is characterized in that, described the first locating module specifically is used for: when adjusting the transmitting power of described base station, when described number is 0, increase step by step the transmitting power of described base station; When described number during greater than described preset value, reduce step by step the transmitting power of described base station.
9. such as claim 6,7 or 8 described devices, it is characterized in that, also comprise:
The second locating module is used for the grid number gathered when described target grid greater than 1 the time, utilizes drive test frequency sweep data to determine the target grid at interference source place.
10. device as claimed in claim 9 is characterized in that, described the second locating module specifically is used for:
From the drive test frequency sweep track of described interference region, choose a predetermined number drive test point, and obtain interference signal intensity corresponding to each described drive test point;
For each grid in the described target grid set, adopt respectively the ray trace propagation model, when calculating transmits from this grid, the received signal strength of each described drive test point, and calculate the described interference signal intensity of all drive test points and the mean square deviation between the described received signal strength;
Determine that grid corresponding to the minimum value of described mean square deviation is described target grid in the described target grid set.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110321864.XA CN103067931B (en) | 2011-10-20 | 2011-10-20 | The method of interference source and device in a kind of positioned radio network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110321864.XA CN103067931B (en) | 2011-10-20 | 2011-10-20 | The method of interference source and device in a kind of positioned radio network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103067931A true CN103067931A (en) | 2013-04-24 |
| CN103067931B CN103067931B (en) | 2015-10-07 |
Family
ID=48110356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110321864.XA Active CN103067931B (en) | 2011-10-20 | 2011-10-20 | The method of interference source and device in a kind of positioned radio network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103067931B (en) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103986536A (en) * | 2014-05-29 | 2014-08-13 | 中国人民解放军军械工程学院 | Method for generating equally-segmented ray beams for positive ray tracing |
| CN105188102A (en) * | 2015-09-24 | 2015-12-23 | 上海斐讯数据通信技术有限公司 | Wireless channel automatic selecting method and system |
| CN105530701A (en) * | 2015-12-11 | 2016-04-27 | 北京北方烽火科技有限公司 | Interference source positioning method and device |
| CN105722135A (en) * | 2016-03-08 | 2016-06-29 | 福建泰克通信有限公司 | Border base station road test system |
| CN106230525A (en) * | 2016-07-27 | 2016-12-14 | 汪瑞 | A kind of radio interference source localization method |
| WO2017000498A1 (en) * | 2015-06-30 | 2017-01-05 | 中兴通讯股份有限公司 | Method and device for locating interference source |
| CN106332271A (en) * | 2015-06-30 | 2017-01-11 | 中国移动通信集团公司 | Interference source location method and device |
| CN106714227A (en) * | 2016-08-31 | 2017-05-24 | 成都中星世通电子科技有限公司 | Automatic positioning method for mobile monitoring station |
| CN107889049A (en) * | 2016-09-30 | 2018-04-06 | 中国移动通信有限公司研究院 | A kind of interference positioning method and positioner |
| CN108260146A (en) * | 2018-01-24 | 2018-07-06 | 深圳市科虹通信有限公司 | A kind of method and its system based on the optimization of interference matrix double-flow design |
| 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 |
| CN108540252A (en) * | 2017-03-02 | 2018-09-14 | 中国移动通信集团广东有限公司 | A kind of neighborhood disturbance treatment method and device |
| CN109314868A (en) * | 2017-10-24 | 2019-02-05 | 深圳市大疆创新科技有限公司 | Interference source localization method, unmanned plane and the computer readable storage medium of base station |
| CN109996239A (en) * | 2017-12-29 | 2019-07-09 | 中国移动通信集团设计院有限公司 | A kind of method and apparatus of LTE network structural optimization analysis |
| CN111246564A (en) * | 2018-11-28 | 2020-06-05 | 中国移动通信集团浙江有限公司 | External interference positioning method and device based on MR data |
| CN111278040A (en) * | 2018-12-05 | 2020-06-12 | 中国移动通信集团四川有限公司 | Interference source positioning method, device, equipment and computer storage medium |
| CN111294920A (en) * | 2018-12-10 | 2020-06-16 | 中国移动通信集团四川有限公司 | Method, device, equipment and medium for identifying interference in LTE network |
| CN111372275A (en) * | 2018-12-25 | 2020-07-03 | 中国移动通信集团有限公司 | Method, device and equipment for positioning interference source |
| CN112731361A (en) * | 2020-12-21 | 2021-04-30 | 成都瀚德科技有限公司 | Bicycle positioning method and device of radio monitoring vehicle |
| CN113015200A (en) * | 2019-12-18 | 2021-06-22 | 中国移动通信集团四川有限公司 | Method, system and device for positioning interference source |
| CN114286375A (en) * | 2021-12-16 | 2022-04-05 | 北京邮电大学 | Mobile communication network interference positioning method |
| CN114449649A (en) * | 2020-11-06 | 2022-05-06 | 中国移动通信集团设计院有限公司 | Interference source positioning method and device based on MRO data |
| WO2023072050A1 (en) * | 2021-10-28 | 2023-05-04 | 中兴通讯股份有限公司 | Interference source positioning method and apparatus, and computer-readable medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101692734A (en) * | 2009-09-29 | 2010-04-07 | 深圳市优网科技有限公司 | Method and device for predicting cell interference |
| WO2010106556A2 (en) * | 2009-03-20 | 2010-09-23 | Centre Of Excellence In Wireless | Cognitive interference management in wireless networks with relays, macro cells, micro cells, pico cells and femto cells |
-
2011
- 2011-10-20 CN CN201110321864.XA patent/CN103067931B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010106556A2 (en) * | 2009-03-20 | 2010-09-23 | Centre Of Excellence In Wireless | Cognitive interference management in wireless networks with relays, macro cells, micro cells, pico cells and femto cells |
| CN101692734A (en) * | 2009-09-29 | 2010-04-07 | 深圳市优网科技有限公司 | Method and device for predicting cell interference |
Non-Patent Citations (1)
| Title |
|---|
| 成艳: "基于射线跟踪模型的GSM网络干扰分析预测研究", 《中国优秀硕士学位论文全文数据库，信息科技辑》 * |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103986536B (en) * | 2014-05-29 | 2016-04-13 | 中国人民解放军军械工程学院 | A kind of forward ray tracings etc. separate the generation method of beam |
| CN103986536A (en) * | 2014-05-29 | 2014-08-13 | 中国人民解放军军械工程学院 | Method for generating equally-segmented ray beams for positive ray tracing |
| WO2017000498A1 (en) * | 2015-06-30 | 2017-01-05 | 中兴通讯股份有限公司 | Method and device for locating interference source |
| CN106332271A (en) * | 2015-06-30 | 2017-01-11 | 中国移动通信集团公司 | Interference source location method and device |
| CN106332154A (en) * | 2015-06-30 | 2017-01-11 | 中兴通讯股份有限公司 | Interference source positioning method and apparatus |
| CN106332271B (en) * | 2015-06-30 | 2020-03-17 | 中国移动通信集团公司 | Interference source positioning method and device |
| CN105188102A (en) * | 2015-09-24 | 2015-12-23 | 上海斐讯数据通信技术有限公司 | Wireless channel automatic selecting method and system |
| CN105188102B (en) * | 2015-09-24 | 2019-07-12 | 上海斐讯数据通信技术有限公司 | A kind of wireless channel automatic selecting method and system |
| CN105530701B (en) * | 2015-12-11 | 2018-12-18 | 北京北方烽火科技有限公司 | A kind of interference source localization method and device |
| CN105530701A (en) * | 2015-12-11 | 2016-04-27 | 北京北方烽火科技有限公司 | Interference source positioning method and device |
| CN105722135A (en) * | 2016-03-08 | 2016-06-29 | 福建泰克通信有限公司 | Border base station road test system |
| CN106230525A (en) * | 2016-07-27 | 2016-12-14 | 汪瑞 | A kind of radio interference source localization method |
| CN106714227B (en) * | 2016-08-31 | 2019-09-10 | 成都中星世通电子科技有限公司 | Mobile monitoring station automatic positioning method |
| CN106714227A (en) * | 2016-08-31 | 2017-05-24 | 成都中星世通电子科技有限公司 | Automatic positioning method for mobile monitoring station |
| CN107889049A (en) * | 2016-09-30 | 2018-04-06 | 中国移动通信有限公司研究院 | A kind of interference positioning method and positioner |
| CN108540252A (en) * | 2017-03-02 | 2018-09-14 | 中国移动通信集团广东有限公司 | A kind of neighborhood disturbance treatment method and device |
| CN108540252B (en) * | 2017-03-02 | 2020-01-17 | 中国移动通信集团广东有限公司 | Cell interference processing method and device |
| CN109314868B (en) * | 2017-10-24 | 2022-03-22 | 深圳市大疆创新科技有限公司 | Interference source positioning method of base station, unmanned aerial vehicle and computer readable storage medium |
| CN109314868A (en) * | 2017-10-24 | 2019-02-05 | 深圳市大疆创新科技有限公司 | Interference source localization method, unmanned plane and the computer readable storage medium of base station |
| CN109996239A (en) * | 2017-12-29 | 2019-07-09 | 中国移动通信集团设计院有限公司 | A kind of method and apparatus of LTE network structural optimization analysis |
| CN109996239B (en) * | 2017-12-29 | 2022-06-03 | 中国移动通信集团设计院有限公司 | Method and device for optimizing and analyzing LTE (Long term evolution) network structure |
| CN108260146A (en) * | 2018-01-24 | 2018-07-06 | 深圳市科虹通信有限公司 | A kind of method and its system based on the optimization of interference matrix double-flow design |
| 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 |
| CN111246564A (en) * | 2018-11-28 | 2020-06-05 | 中国移动通信集团浙江有限公司 | External interference positioning method and device based on MR data |
| CN111278040A (en) * | 2018-12-05 | 2020-06-12 | 中国移动通信集团四川有限公司 | Interference source positioning method, device, equipment and computer storage medium |
| CN111294920B (en) * | 2018-12-10 | 2021-05-18 | 中国移动通信集团四川有限公司 | Method, device, equipment and medium for identifying interference in LTE network |
| CN111294920A (en) * | 2018-12-10 | 2020-06-16 | 中国移动通信集团四川有限公司 | Method, device, equipment and medium for identifying interference in LTE network |
| CN111372275A (en) * | 2018-12-25 | 2020-07-03 | 中国移动通信集团有限公司 | Method, device and equipment for positioning interference source |
| CN113015200A (en) * | 2019-12-18 | 2021-06-22 | 中国移动通信集团四川有限公司 | Method, system and device for positioning interference source |
| CN114449649A (en) * | 2020-11-06 | 2022-05-06 | 中国移动通信集团设计院有限公司 | Interference source positioning method and device based on MRO data |
| CN114449649B (en) * | 2020-11-06 | 2023-09-05 | 中国移动通信集团设计院有限公司 | Interference source positioning method and device based on MRO data |
| CN112731361A (en) * | 2020-12-21 | 2021-04-30 | 成都瀚德科技有限公司 | Bicycle positioning method and device of radio monitoring vehicle |
| CN112731361B (en) * | 2020-12-21 | 2024-04-16 | 成都瀚德科技有限公司 | Bicycle positioning method and device for radio monitoring bicycle |
| WO2023072050A1 (en) * | 2021-10-28 | 2023-05-04 | 中兴通讯股份有限公司 | Interference source positioning method and apparatus, and computer-readable medium |
| CN114286375A (en) * | 2021-12-16 | 2022-04-05 | 北京邮电大学 | Mobile communication network interference positioning method |
| CN114286375B (en) * | 2021-12-16 | 2023-08-18 | 北京邮电大学 | Mobile communication network interference positioning method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103067931B (en) | 2015-10-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103067931B (en) | The method of interference source and device in a kind of positioned radio network | |
| CN101631349B (en) | Method, device and wireless operation maintenance center for positioning terminal | |
| CN100403052C (en) | Method and system for estimating position of mobile device | |
| KR100635319B1 (en) | Method and system for determining the speed and distance of a mobile unit | |
| US8131312B2 (en) | Method and system for construction of radio environment model | |
| US20170338901A1 (en) | Method for predicting indoor three-dimensional space signal field strength using an outdoor-to-indoor propagation model | |
| US11800373B2 (en) | System and method for 3D propagation modelling for planning of a radio network | |
| CN102869038B (en) | Measurement compensation method and device in positioning of base station, server and system | |
| CN101090301B (en) | Radio wave route loss simulation measuring method | |
| CN109862505A (en) | Car tracing positioning system and tracking positioning method in tunnel based on bus or train route collaboration | |
| CN101777958B (en) | Method for forecasting group delays within certain range near receiving point | |
| Li et al. | Cramer-rao lower bound analysis of data fusion for fingerprinting localization in non-line-of-sight environments | |
| CN102438251B (en) | Method for calculating declination angle of antenna of mobile communication base station | |
| EP1302783A1 (en) | Location of a mobile terminal | |
| Hein et al. | Synergies between satellite navigation and location services of terrestrial mobile communication | |
| CN110118953A (en) | A kind of localization method, apparatus and system | |
| Wang et al. | RSSI-based localization in cellular networks | |
| RU2657237C1 (en) | One-way method of the radio frequency sources location | |
| RU2444705C1 (en) | System for real-time location of mobile objects | |
| Corre et al. | A new approach for radio propagation modeling in urban environment: Knife-edge diffraction combined with 2D ray-tracing | |
| PILEGGI | 5G positioning: preliminary analysis of early data sets | |
| US7231217B2 (en) | Determining neighbour lists | |
| KR101495503B1 (en) | Method for Estimating Azimuth, Apparatus And Computer-Readable Recording Medium with Program Therefor | |
| KR20160090722A (en) | Informing Method of Location in Indoor Space | |
| CN104936198B (en) | A kind of method and apparatus of determining neighboring BS relationship |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |