CN106921989B

CN106921989B - A kind of communication network field strength distribution determines method and device

Info

Publication number: CN106921989B
Application number: CN201510998126.7A
Authority: CN
Inventors: 徐海勇; 高随祥; 孙伟; 毕文通; 姜志鹏; 范政; 王晋龙; 周欣; 宋皓明
Original assignee: China Mobile Group Beijing Co Ltd
Current assignee: China Mobile Group Beijing Co Ltd
Priority date: 2015-12-25
Filing date: 2015-12-25
Publication date: 2019-09-27
Anticipated expiration: 2035-12-25
Also published as: CN106921989A

Abstract

The invention discloses a kind of communication network field strength distributions to determine method and device.The method of the present invention includes: the outdoor grid obtained in target area, wherein outdoor grid refers to that the grid for not falling within interior of building, a grid are the area of space being sized；According to the spatial position of grid where the objective emission antenna in outdoor grid and target area, the path loss of determining outdoor grid to objective emission antenna；According to the transmission power of objective emission antenna and outdoor grid to the path loss of the objective emission antenna, field strength of the objective emission antenna in outdoor grid is determined.Further, field strength of the objective emission antenna indoors in grid can also be determined in the field strength in outdoor grid based on identified objective emission antenna.The present invention, which can establish one kind, can reach higher field intensity prediction precision and avoid excessively huge calculation amount, while also can be suitably used for the field-strength prediction model of multiple network several scenes.

Description

A kind of communication network field strength distribution determines method and device

Technical field

The present invention relates to fields of communication technology more particularly to a kind of communication network field strength distribution to determine method and device.

Background technique

With the continuous growth of service of mobile communication, the demand of network capacity constantly expands, due to having for frequency resource Limit is fundamentally needing to wirelessly communicate to meet the needs of network capacity further expansion and guarantee the quality of network simultaneously System can be realized the reasonable effective network coverage, and rationally the effective network coverage depends on the field intensity prediction of degree of precision. It therefore, is design cordless communication network and mobile communication system engineering design to the field intensity prediction of network signal in communication network Scientific basis.

The field intensity prediction of communication network is usually to utilize typical radio waves propagation model, to the field in communication network coverage area Strong distribution, propagated loss and communication probability are measured and are estimated.In field of communication technology, from the 1950s, Scholar both domestic and external establishes a variety of field datas that are based on successively and statisticallys analyze resulting field-strength prediction model, mainly includes Empirical model, half determine model and deterministic models.Empirical model mainly estimates model parameter by measured data, does not consider Diffraction and reflection, realization is relatively easy, but precision is lower；Half determines model mainly for certain in the area of field intensity prediction Special scene (such as scene shown in FIG. 1) provides the calculation formula of the fining based on ray theory, thus to these areas Domain carries out field intensity prediction, and application scenarios are limited；Deterministic models mainly pass through such as accurate ray-tracing procedure, send out from antenna Each point that end can be influenced to antenna is penetrated, the electromagnetism wave ray emulation refined, calculating is with high costs, also can not be big Scale uses.

Accordingly, with respect to existing field-strength prediction model, higher field intensity prediction precision can be reached by how establishing one kind And excessively huge calculation amount is avoided, while also can be suitably used for the field-strength prediction model of multiple network several scenes, it is industry institute Urgently studies and solve the problems, such as.

Summary of the invention

The embodiment of the present invention provides a kind of communication network field strength distribution and determines method and device, can reach to establish one kind To higher field intensity prediction precision and excessively huge calculation amount is avoided, while also can be suitably used for the field of multiple network several scenes Strong prediction model.

The communication network field strength distribution that one embodiment of the present of invention provides determines method, comprising:

The outdoor grid in target area is obtained, the outdoor grid refers to the grid for not falling within interior of building, one Grid is the area of space being sized；

According to the spatial position of grid where the objective emission antenna in the outdoor grid and the target area, really Path loss of the fixed outdoor grid to the objective emission antenna；

According to the transmission power of the objective emission antenna and the outdoor grid to the road of the objective emission antenna Diameter loss determines field strength of the objective emission antenna in the outdoor grid.

Wherein, according to the space bit of grid where the objective emission antenna in the outdoor grid and the target area Set, determine the outdoor grid to the objective emission antenna path loss, comprising:

By n times iterative process, be calculated each outdoor grid to the objective emission antenna ray path parameter, N is the integer more than or equal to 1；

According to each outdoor grid to the ray path parameter of the objective emission antenna, each outdoor is determined respectively Path loss of the grid to the objective emission antenna；

Wherein, include: in each iterative process

Determine the reference grid during current iteration；Wherein, the reference grid in first time iterative process is the mesh The outdoor grid where transmitting antenna is marked, the reference grid for the second time and in later iterative process is in a preceding iterative process Reference grid adjacent chamber outside grid；

According to where the relatively described objective emission antenna of grid outside the ray path parameter of reference grid and adjacent chamber The spatial deviation amount of outdoor grid calculates the ray path parameter of grid outside adjacent chamber；Wherein, in first time iterative process Reference grid ray path parameter be setting value.

Specifically, according to the relatively described objective emission day of grid outside the ray path parameter of reference grid and adjacent chamber The spatial deviation amount of outdoor grid where line calculates the ray path parameter of grid outside adjacent chamber, comprising:

To grid outside an adjacent chamber, the benchmark grid adjacent in the every one-dimensional square in space with grid outside the adjacent chamber are obtained The ray path parameter of lattice, and outdoor grid where the relatively described objective emission antenna of grid outside the adjacent chamber is obtained in sky Between spatial deviation amount in every one-dimensional square；

The ray path parameter of grid outside an adjacent chamber is calculated according to following formula:

Wherein, C_pathIndicate the ray path parameter of grid outside an adjacent chamber being calculated,Indicate that this is adjacent The ray path parameter of outdoor grid reference grid adjacent on the i-th dimension direction of space, step⁽ⁱ⁾Indicate grid outside the adjacent chamber Spatial deviation amount of the outdoor grid on the i-th dimension direction of space where the relatively described objective emission antenna of lattice, D is space Dimension, 1≤i≤D, i, D are positive integer.

Wherein, which is determined to the ray path parameter of the objective emission antenna according to an outdoor grid To the path loss of the objective emission antenna, comprising:

Ray path parameter, the outdoor grid to the target according to an outdoor grid to the objective emission antenna Azimuth and angle of declination of the distance and the outdoor grid of grid to the objective emission antenna where transmitting antenna, determine Path loss of the outdoor grid to the objective emission antenna.

Specifically, it is arrived according to an outdoor grid to the ray path parameter of the objective emission antenna, the outdoor grid The distance of grid where the objective emission antenna and the outdoor grid to the objective emission antenna azimuth with have a down dip Angle, determine the outdoor grid to the objective emission antenna path loss, comprising:

Determine according to the following formula an outdoor grid to the objective emission antenna path loss:

L=a+b × lg (d)+c × ln (C_path)+c_w×f(α,β)

Wherein, L indicates an outdoor grid to the path loss of the objective emission antenna, and d is the outdoor grid to institute State the linear distance of objective emission antenna, C_pathFor the ray path parameter value of the outdoor grid, α is the outdoor grid described in The azimuth of objective emission antenna, β are angle of declination of the outdoor grid to the objective emission antenna, and f (α, β) is the target The antenna waveform diagram attenuation function of transmitting antenna, a, b, c, c_wFor constant coefficient.

Further, according to the transmission power of the objective emission antenna and the outdoor grid to the objective emission The path loss of antenna determines field strength of the objective emission antenna in the outdoor grid, comprising:

According to the following formula, field strength of the objective emission antenna in the outdoor grid is determined:

P=sendpower-L

Wherein, P is field strength of the objective emission antenna in an outdoor grid, and sendpower is objective emission antenna Transmission power, L are path loss of the objective emission antenna to the outdoor grid.

Preferably, determine the objective emission antenna after the field strength in the outdoor grid, further includes:

According to field strength of the identified objective emission antenna in each outdoor grid, institute is determined according to following formula Field strength of the objective emission antenna in the indoor grille in target area is stated, an indoor grille, which refers to, falls into interior of building Grid:

Wherein, P_inFor field strength of the objective emission antenna in an indoor grille, P_jTo be built with where the indoor grille The field strength of grid, η outside jth room in all J adjacent outdoor grids of object₀Indicate that electromagnetic wave once penetrates the decaying number of wall Value, η₁Indicate that unit distance stops body attenuation value, d_jIndicate will the distance between grid outside the indoor grille and jth room,It indicates to choose all J values being calculated maximum value, j, J are the integer more than or equal to 0.

The communication network field strength distribution determining device that one embodiment of the present of invention provides, comprising:

Module is obtained, for obtaining the outdoor grid in target area, the outdoor grid, which refers to, not to be fallen in building The grid in portion, a grid are the area of space being sized；

First determining module, for according to where the objective emission antenna in the outdoor grid and the target area The spatial position of grid, determine the outdoor grid to the objective emission antenna path loss；

Second determining module, for the transmission power and the outdoor grid according to the objective emission antenna described in The path loss of objective emission antenna determines field strength of the objective emission antenna in the outdoor grid.

Wherein, the first determining module is specifically used for:

Wherein, include: in each iterative process

In some embodiments, the first determining module is specifically used for:

Wherein, the first determining module is specifically used for:

In some embodiments, the first determining module is specifically used for:

L=a+b × lg (d)+c × ln (C_path)+c_w×f(α,β)

Wherein, the second determining module is specifically used for:

P=sendpower-L

Preferably, communication network field strength distribution determining device provided by one embodiment of the present of invention, further includes: third Determining module, for the field strength according to the identified objective emission antenna in each outdoor grid, according to following formula Determine field strength of the objective emission antenna in the indoor grille in target area, an indoor grille, which refers to, falls into building Internal grid:

The embodiment of the invention provides a kind of communication network field strength distributions to determine method and device, the embodiment of the present invention energy A kind of field-strength prediction model is enough provided, the model is based on the rasterizing to target area, according to outdoor grid and target area The spatial position of grid where interior objective emission antenna, the path loss of determining outdoor grid to the objective emission antenna, Transmission power and outdoor grid further according to objective emission antenna determine objective emission to the path loss of objective emission antenna Field strength of the antenna in outdoor grid.It can be seen that determining in a kind of communication network field strength distribution provided in an embodiment of the present invention It can determine the distribution of communication network field strength in target area in method, it is contemplated that building distribution and the influence of topography and geomorphology, Applicable scene is not strict with simultaneously, it is substantially accurate only to need to building and terrain information, and it is usual to calculate cost In feasible range, therefore it can be realized and a kind of can reach higher field intensity prediction precision and avoid excessively huge calculating Amount, while also can be suitably used for the field-strength prediction model of multiple network several scenes.

Detailed description of the invention

To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment Attached drawing is briefly introduced, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this For the those of ordinary skill in field, without any creative labor, it can also be obtained according to these attached drawings His attached drawing.

Fig. 1 is a kind of application scenarios schematic diagram that field intensity prediction half determines model in the prior art；

Fig. 2 is an example of rasterizing map in the embodiment of the present invention；

Fig. 3 is the example that rasterizing map is converted into connected graph in the embodiment of the present invention；

Fig. 4 is an example at public the parameter azimuth and angle of declination in of antenna in the embodiment of the present invention；

Fig. 5 is one of the public parameter of antenna in the embodiment of the present invention based on the decaying of the waveform diagram of azimuth and angle of declination in A example；

Fig. 6 is the flow diagram that the communication network field strength distribution that one embodiment of the present of invention provides determines method；

Fig. 7 determines for the communication network field strength distribution that one embodiment of the present of invention provides calculates each outdoor grid in method The iterative process schematic diagram of the ray path parameter of lattice；

Fig. 8 is a kind of example for iterative process initialization in two-dimensional space that one embodiment of the present of invention provides；

Fig. 9 is the example for the update reference grid that one embodiment of the present of invention provides；

Figure 10 is the example updated to lattice-shaped state value outside adjacent chamber that one embodiment of the present of invention provides；

Figure 11 is the example to raster symbol-base parameter outside adjacent chamber that one embodiment of the present of invention provides；

Figure 12 is the example for confirming the ray path parameter of grid outside adjacent chamber that one embodiment of the present of invention provides；

Figure 13 is the example for the update reference grid that one embodiment of the present of invention provides；

Figure 14 is the example to raster symbol-base parameter outside adjacent chamber that one embodiment of the present of invention provides；

Figure 15 is the example for confirming the ray path parameter of grid outside adjacent chamber that one embodiment of the present of invention provides；

Figure 16 is a kind of showing for the updated reference grid set of iteration several times that one embodiment of the present of invention provides Example；

Figure 17 is an adjacent gate of the updated reference grid of iteration several times that one embodiment of the present of invention provides The example of the ray path parameter of lattice；

Figure 18 is the result example that a kind of iterative process that one embodiment of the present of invention provides is completed；

The example of indoor grille in a kind of two-dimensional scene that Figure 19 provides for one embodiment of the present of invention；

Figure 20 is the example that an indoor field strength is calculated based on outdoor grid field strength that one embodiment of the present of invention provides；

Figure 21 is the structural schematic diagram for the communication network field strength distribution determining device that one embodiment of the present of invention provides.

Specific embodiment

To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with attached drawing to the present invention make into It is described in detail to one step, it is clear that the described embodiments are only some of the embodiments of the present invention, rather than whole implementation Example.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts All other embodiment, shall fall within the protection scope of the present invention.

The research main purpose of field-strength prediction model is the planning of science activities in order to realize communication network, the reasonable siting of station Covering can be effectively provided, eliminate coverage hole, and provide convenience for the network optimization from now on, at present for field intensity prediction, It has studied both at home and abroad for many years, there is many research achievements.Field-strength prediction model is broadly divided into following 3 class:

1, empirical model: this class model is mainly classified by the area to field intensity prediction, for example is divided into city, suburb Area, rural area etc. provide different field intensity prediction formula further according to different classifications.This class model mainly includes Hata mould Type, Egli model, Cost231 model etc..

However although empirical model can quickly calculate the field strength of each point as a result, the still distribution and feature of building Coverage prediction result there is no and have an impact, antenna radiation pattern and day are only mainly considered in coverage prediction result The influence of line height etc..Therefore for many places, forecasting accuracy is very poor.Meanwhile in practical application, for field strength to be carried out The area of prediction, it is generally difficult to which this area is classified as a certain seed type.Therefore, empirical model, which is applicable in scene, certain limitation. Meanwhile empirical model can only mostly carry out the prediction of two-dimentional field strength.

2, half determine model: this class model provides base primarily directed to certain special scenes in the area of field intensity prediction In the calculation formula of the fining of ray theory, field intensity prediction is carried out to these regions.

Semidefiniteness model can carry out Accurate Prediction to certain scenes, but in practical application, building distribution and ground Shape is extremely complex, and most of place is difficult to be classified as a certain scene.Therefore, small part area is only applicable in practical application Field intensity prediction verifying, scene as shown in Figure 1 can be according to the antenna height got, the location parameter of specific region, angle Parameter etc. carries out the field intensity prediction in the specific region.

3, deterministic models: this class model is mainly with accurate ray-tracing procedure, to from antenna transmitting terminal to antenna The each point that can be influenced, the electromagnetism wave ray emulation refined.Ray tracing technique is that optical ray technology exists Application in electromagnetism calculating field, can accurately take into account the various routes of transmission of electromagnetic wave, including direct projection, reflection, diffraction, Transmission etc., and the various factors in view of influencing radio wave propagation, to do accurate prediction for different concrete scenes.It is this kind of Model mainly includes Volcano ray tracing model, WaveSight model and WinProp model etc..

Although deterministic models can accurately predict the field strength of any point as a result, still ray-tracing procedure is practical Calculation amount is very big when use, and operation time is too long, can only calculate small region, is difficult real in large-scale area With.Meanwhile ray-tracing procedure, to building information, the required precision of topography profile is very high, while also needing ground, wall material The information that matter etc. hardly results in.And the software package that ray-tracing algorithm mostly uses offshore company packaged at present, nothing Method obtains specific algorithm and calculation formula.

Aiming at the shortcomings in the prior art, the embodiment provides a kind of communication network field strength distributions to determine method And device.A kind of communication network field strength distribution provided by embodiment through the invention determines method and device, can establish For one kind compared with empirical model, calculated result accuracy is high, it is contemplated that building distribution and the influence of topography and geomorphology, while can be with Carry out the field-strength prediction model of three-dimensional field intensity prediction；And for the model compared with semidefiniteness model, applicability is wide, to any Any scene in area is all agreed to be applicable in；Meanwhile it is feasible to calculate the time compared with ray tracing model for the model, and to building The precision of object and terrain information does not have very high requirement.

A kind of communication network field strength distribution provided by the embodiment of the present invention determines the shifting that method and device can be convenient The multiple network that 2G network, 3G network, 4G network even 5G network etc. may develop in the future is planted, can be subsequent net Network application, which lays the foundation, can reach higher field intensity prediction precision and avoid excessively huge calculation amount, while also can be suitably used for The field-strength prediction model of several scenes.

When should be pointed out that practical application field-strength prediction model in engineering, need first to include basic data preparation. Method and device is determined for a kind of communication network field strength distribution provided by the embodiment of the present invention, it is usually required mainly for prepare following base Plinth data: the high accuracy three-dimensional electronic map of target area；The basic work parameter evidence of antenna is calculated in target area；And The frequency sweep data of target area.

Wherein, the high accuracy three-dimensional electronic map of target area can be by such as providing storage, showing, analytically It manages the softwares such as the GIS (Geographic Information System, GIS-Geographic Information System) of spatial data function to obtain, ground Diagram data is capable of providing the shape of building, height, the shape of road, direction, and width etc. reflects target area landform landforms Geography information.

Based on the electronic map of the target area got, a kind of communication network field strength point provided by the embodiment of the present invention Cloth determine method and device also need in advance to the high accuracy three-dimensional electronic map of the target area got need first into Row three-dimensional map rasterizing, this is that a kind of communication network field strength distribution provided by the embodiment of the present invention is determined method and device Realize basis in a particular application.

Such as in some specific embodiments, based on accessed electronic map, it can be counted first by required The target area of calculation carries out a × a × a rasterizing and divides, i.e., being averagely divided into several sizes is a × a × a grid, Division direction is parallel with solid axes.Wherein, the size of a can according to demand and/or accessed electronic three-dimensional map Precision determine, such as 1 meter, 5 meters or 20 meters different numerical value.

Fig. 2 shows a kind of by the scene of target area progress rasterizing division for 5 meters, with one in target area Building for the building that every layer of size is 200m × 200m, can be divided into 7 layers according to 5 meters one layer by 35 meters of height, therein every One layer of planar grid for being divided into 1600 5m × 5m.Entire building are divided into 11200 5m × 5m × 5m space lattices.

It divides in obtained all grids, has in a kind of grid it can be seen that rasterizing will be carried out to target area Heart point is in the polyhedron where falling in building, and the central point of also a kind of grid is then the polyhedron where falling in building Outside, the communication network field strength distribution provided by the embodiment of the present invention determines in method, and central point is fallen in where building The intracorporal grid of multi-panel is considered indoor grille, and grid central point fallen in outside the polyhedron where building is considered outdoor Grid.

For any one transmitting antenna in the target area, all non-building grids are regarded as a point, each Connect a line between point point corresponding with grid adjacent cells where it.Then the regioinvertions that a transmitting antenna influences are become One connected graph.Fig. 3 shows map rasterizing and is converted into an example of connected graph.Wherein, the grid of Fig. 3 Oxford gray Lattice are indoor grille, and white grid is outdoor grid, and the line per adjacent two between the central point of white grid is used for Indicate that this two white grids are connected.It can be seen that only connected graph, indoor grille can be converted into for outdoor grid Due to being located in building, cannot be connected to outdoor grid.

It should be pointed out that in the description of the present invention, a transmitting antenna indicates that one or more belongs to the same base The transmitting antenna for cell of standing a, it can transmitting antenna is interpreted as a base station cell.

In addition, the work parameter evidence of the antenna of required acquisition, being capable of providing includes position (such as longitude and latitude where antenna Degree), antenna height, transmission power, antenna gain, antenna model, horizontal azimuth, vertical angle of declination, whether be microcellulor day The data informations such as line.The frequency sweep data of required acquisition are capable of providing including the corresponding antenna model of each data, each data Measure the data informations such as place (such as longitude and latitude), measured value, time of measuring.

Accordingly and the map datum that gets according to the work parameter of antenna, for any in the target area that is rasterized into Grid can be determined by the position of the grid central point of this in map datum and the position of transmitting antenna launch point The central point of the grid to the linear distance between the launch point of transmitting antenna, and can also determine the grid central point with Angle between line between the launch point of transmitting antenna and antenna principal direction, including horizontal sextant angle, i.e. azimuth, and it is vertical Angle, i.e. angle of declination.Fig. 4 shows an outdoor grid to the azimuth of transmitting antenna and the example of angle of declination, wherein α shown in (a) in Fig. 4 is the azimuth of the outdoor grid to objective emission antenna in the horizontal direction, (b) in Fig. 4 Shown in β be the angle of declination of the outdoor grid to objective emission antenna in vertical direction.

Further, decayed according to the work parameter of antenna according to antenna waveform diagram can also be obtained by the model of antenna.Fig. 5 Show an example based on the decaying of the waveform diagram of azimuth and angle of declination, wherein (a) in Fig. 5 is azimuth and decaying The functional relation of value, abscissa are azimuth, and ordinate is pad value, and (b) in Fig. 5 is the function pass of angle of declination and pad value System, abscissa is angle of declination, and ordinate is pad value.For a grid, by azimuth and angle of declination and pass through antenna Model obtain antenna waveform diagram decaying, can determine on transmitting antenna to the grid direction by azimuth and angle of declination Caused path attenuation value, the pad value decline equal to the horizontal pad value as caused by azimuth is vertical with as caused by angle of declination The sum of depreciation.

Based on above pretreatment work, Fig. 6 shows a kind of communication network field of one embodiment of the present of invention offer The flow diagram of distribution determination method, the process can be realized by the combination of computer software or software and hardware by force, the stream Journey includes the following steps:

Step 601: obtaining the outdoor grid in target area, the outdoor grid refers to the grid for not falling within interior of building Lattice, a grid are the area of space being sized.

Step 602: according to the spatial position of grid where the objective emission antenna in outdoor grid and target area, really Path loss of the fixed outdoor grid to objective emission antenna.

Step 603: being damaged according to the transmission power of objective emission antenna and outdoor grid to the path of objective emission antenna Consumption, determines field strength of the objective emission antenna in outdoor grid.

Wherein, step 601 can be is realized more based on obtained rasterizing maps in pretreatment work, by locating in advance Accessed rasterizing map, obtains the target area of rasterizing, so as to get in the target area in science and engineering work Outdoor grid and indoor grille, specifically, each grid can be identified by the center position of the grid, such as The longitude and latitude of central point.

Specifically, described in step 602, according to where the objective emission antenna in outdoor grid and target area The spatial position of grid determines the realization of path loss of the outdoor grid to objective emission antenna, can be through n times iteration mistake Each outdoor grid is calculated to the ray path parameter of objective emission antenna in journey first, and N is the integer more than or equal to 1； Further according to the ray path parameter of each outdoor grid to objective emission antenna, to determine each outdoor grid to mesh respectively Mark the path loss of transmitting antenna.

Wherein, include in each iterative process: determining reference grid during current iteration first；Wherein, first Reference grid in secondary iterative process is the outdoor grid where objective emission antenna, for the second time and in later iterative process Reference grid is grid outside the adjacent chamber of the reference grid in a preceding iterative process；Further according to the ray path ginseng of reference grid The spatial deviation amount of outdoor grid outside several and adjacent chamber where grid relative target transmitting antenna calculates adjacent chamber The ray path parameter of outer grid；Wherein, the ray path parameter of the reference grid in first time iterative process is setting value.

Specifically, according to grid relative target transmitting antenna institute outside the ray path parameter of reference grid and adjacent chamber Outdoor grid spatial deviation amount, calculate adjacent chamber outside grid ray path parameter, comprising:

Due to for iterative calculation, the characteristics of according to interative computation, it should be readily understood that can pass through for formula (1) By target area average grid in preprocessing process, and assume for indicating the relatively described objective emission of outdoor grid The step of spatial deviation amount of the outdoor grid on the i-th dimension direction of space where antenna⁽ⁱ⁾Using a step as linear module, One step is size of the grid in the dimension, is calculated to simplify.Formula (1) can be converted to following form to simplify meter It calculates:

Wherein,Indicate grid is adjacent on the i-th dimension direction of space outside the adjacent chamber reference grid in space i-th Tie up the spatial deviation amount on direction.

It can be seen that above-mentioned pretreatment and it is assumed that making the outdoor grid where the relatively described objective emission antenna of each grid Spatial deviation amount of the lattice on the i-th dimension direction of space, can directly by interative computation to adjacent thereto previous Spatial deviation amount of the outdoor grid on the i-th dimension direction of space where the relatively described objective emission antenna of grid is added 1 obtains, i.e., indicates the outdoor grid outside an adjacent chamber where the relatively described objective emission antenna of grid in space in formula (1) The step of spatial deviation amount on i-th dimension direction⁽ⁱ⁾, can by grid outside the adjacent chamber on the i-th dimension direction of space phase Spatial deviation amount of the adjacent reference grid on the i-th dimension direction of space1 is added to obtain.

Such as by taking three-dimensional space as an example, by the way that target area to be averagely divided into multiple a × a × a in preprocessing process The grid and step of size⁽ⁱ⁾Linear module be a step, and the width of a step be a.Had according to formula (1) and formula (2):

Based on above explanation, in order to be described herein it is more clear with conveniently, will use formula (2) as one phase of calculating The formula of the ray path parameter of grid outside adjacent room.

Since the coverage of each transmitting antenna is centered on grid where the transmitting antenna, with certain radius (ratio Such as 3 kms) all grids in range, therefore, in some specific embodiments, to the objective emission day line computation mesh When marking transmitting antenna to field strength on each outdoor grid, it is generally preferred that obtain all in objective emission antenna coverage Outdoor grid.

Each outdoor grid is calculated to objective emission antenna about by n times iterative process to above-mentioned in order to clearer The process of ray path parameter namely a kind of process of rasterizing ray tracing are illustrated, below will be by the way that grid is arranged The mode of several auxiliary parameters specifically describes a kind of realization of above-mentioned rasterizing ray tracing process, can specifically pass through calculating The mode of machine programming realizes the rasterizing ray tracing process.

It, can be by the position in each every one-dimensional square of grid installation space in currently known three-dimensional space state Set parameter, it is assumed that each grid is a × a × a size cube, and being with the central point position coordinate of a grid should The spatial position coordinate of grid, and assume that the width that linear module is a step and a step is a, to a grid installation space position Coordinate are as follows: the step number step on the direction reference axis x^(x), step number step on the direction reference axis y^(y), step on the direction reference axis z Number step^(z)。

In addition, being additionally provided with state value to each outdoor grid: flag can be used to identify each outdoor grid and calculate Each outdoor grid is to the state in the interative computation of the ray path parameter of objective emission antenna.

Meanwhile to each outdoor grid, it is provided with the direction x and updates state value: uflag^(x), the direction y update state value: uflag^(y), the direction z updates state value: uflag^(z).It is respectively intended to mark at no point in the update process, each outdoor grid, if quilt The reference grid in its adjacent direction x, y and z updates.If grid is updated by certain one-dimensional reference grid, remaining phase Adjacent reference grid with dimension is no longer updated.

By taking any one objective emission antenna as an example, using the outdoor grid where objective emission antenna launch point as starting point grid Lattice traverse all grids of transmitting antenna influence, i.e., iterative calculation successively according to following process, to be calculated every Ray path parameter of a outdoor grid to the objective emission antenna.Fig. 7 shows a kind of state value mark room based on setting The iterative process schematic diagram of the ray path parameter of outer grid, it is described that specific step is as follows:

Step 701: initialization: first by the ray path parameter C of all outdoor grids_pathIt is set as 0, by all outdoor grid The space position parameter step of lattice^(x)、step^(y)、step^(z)And state value flag is set as -1, uflag^(x), uflag^(y)With uflag^(z)All it is set as 0.The set G of the i.e. reference grid of grid to be processed is set, and set G is initially empty.

By the ray path parameter C of starting point grid (the outdoor grid i.e. where the objective emission antenna transmitting terminal)_pathIf It is set to 1, by the step of starting point grid^(x)、step^(y)And step^(z)It is set as 0, state value flag is set as 0.

Fig. 8 shows a kind of example of rasterizing ray tracing process initialization in two-dimensional space, it should be understood that It is, in embodiments of the present invention for convenience, by the example to all about rasterizing ray tracing with two-dimensional space In example illustrate, and by that should can see the description of the embodiment of the present invention in three-dimensional even more higher-dimension in two-dimensional space Space in, the concrete operations thinking of process described in the embodiment of the present invention is almost the same.

As shown in figure 8, the grid of Dark grey represents building grid, i.e., non-outdoor grid.In each outdoor grid most The number in bracket is the number of the outdoor grid below, and No. 0 grid is starting point grid.When initial, setting starting point grid C_pathEqual to 1, step^(x)And step^(y)It is equal to 0, state value flag is equal to 0, the C of remaining outdoor grid_pathEqual to 0, step^(x) And step^(y)It is equal to -1, state value flag is equal to -1, uflag^(x)And uflag^(y))All it is set as 0.

Step 702: it updates reference grid set G: determining when the reference grid set in previous iterative process, it is specific real It can be now and first empty set G, all flag are then put into G for 0 outdoor grid, then by all outdoors in G The flag of grid is updated to 1, it can be seen that after above-mentioned initialization, the reference grid in iterative calculation is for the first time For the outdoor grid where objective emission antenna.

Based on initialization shown in Fig. 8, Fig. 9 shows the situation for updating reference grid for the first time, at this point, in benchmark set Element be starting point grid, i.e., number be 0 outdoor grid (0), while the state value flag of the outdoor grid (0) being updated to 1。

Step 703: all grids in processing reference grid set G: for each reference grid in set G, obtaining Grid outside each adjacent chamber adjacent thereto, and updating grid flag outside each adjacent chamber is 0.

If the flag of grid is -1 outside an adjacent chamber, the flag for updating grid outside the adjacent chamber is 0.Such as Figure 10 institute Show, reference grid is outdoor grid (0) at this time, and grid is respectively outdoor grid (1), outdoor grid (2) and room outside adjacent chamber Outer grid (3).The flag of these three outdoor grids is shown in (a) in -1, Figure 10 at this time.Therefore these three outdoor grids Flag be updated to 0, as shown in (b) in Figure 10.

Step 704: grid outside the adjacent chamber for being 0 for flag, according to the ray path parameter of reference grid and adjacent The spatial deviation amount of outdoor grid where the relatively described objective emission antenna of outdoor grid calculates grid outside adjacent chamber Ray path parameter.

Specifically, it to grid outside an adjacent chamber, obtains adjacent in the every one-dimensional square in space with grid outside the adjacent chamber Reference grid ray path parameter, and obtain the outdoor where the relatively described objective emission antenna of grid outside the adjacent chamber Spatial deviation amount of the grid in the every one-dimensional square in space.

Wherein, an outdoor grid gn may be only adjacent on the direction reference axis x with a reference grid, or only with One reference grid is adjacent on the direction reference axis y, or only adjacent on the direction reference axis z with a reference grid, can also Can simultaneously with a reference grid on the direction reference axis x it is adjacent and with another reference grid it is adjacent on the direction reference axis y, Or simultaneously with a reference grid on the direction reference axis y it is adjacent and with another reference grid on the direction reference axis z phase Neighbour, or simultaneously with a reference grid on the direction reference axis x it is adjacent and with another reference grid on the direction reference axis z It is adjacent, or simultaneously with a reference grid on the direction reference axis x it is adjacent and with another reference grid in the direction reference axis y It is upper adjacent and adjacent etc. on the direction reference axis z with third reference grid.

Specifically, if grid is g outside an adjacent chamber_n, g_nCurrent ray path parameter isAccording to formula (2):

If outdoor grid g_nWith a reference grid g₁It is adjacent on the direction solid axes x, and g_nCurrent uflag^(x)=0, by g₁Ray path parameter be expressed asThen:

It enables

g_nStep^(x)=max (g_nCurrent step^(x),g₁Step^(x)+1)

g_nStep^(y)=max (g_nCurrent step^(y),g₁Step^(y))

g_nStep^(z)=max (g_nCurrent step^(z),g₁Step^(z))

uflag^(x)=1

If outdoor grid g_nWith a reference grid g₂It is adjacent on the direction solid axes y, and g_nCurrent uflag^(y)=0, by g₂Ray path parameter be expressed asThen:

It enables

g_nStep^(y)=max (g_nCurrent step^(y),g₂Step^(y)+1)

g_nStep^(x)=max (g_nCurrent step^(x),g₂Step^(x))

g_nStep^(z)=max (g_nCurrent step^(z),g₂Step^(z))

uflag^(y)=1

If outdoor grid g_nWith a reference grid g₃It is adjacent on the direction solid axes z, and g_nCurrent uflag^(z)=0, by g₃Ray path parameter be expressed asThen:

It enables

g_nStep^(z)=max (g_nCurrent step^(z),g₃Step^(z)+1)

g_nStep^(x)=max (g_nCurrent step^(x),g₃Step^(x))

g_nStep^(y)=max (g_nCurrent step^(y),g₃Step^(y))

uflag^(z)=1

Specifically, as follows to the calculating of grid outside the adjacent chamber of reference grid to as shown in figure 11:

Wherein, outdoor grid (1) only uflag adjacent in the x direction with reference grid (0) and outdoor grid (1)^(x)= 0, then for outdoor grid (1):

Outdoor grid (1)

The step of outdoor grid (1)^(x)=max (the current step of outdoor grid (1)^(x), the step of reference grid (0)^(x)+ 1)=max (- 1,0+1)=1；The step of outdoor grid (1)^(y)=max (the current step of outdoor grid (1)^(y), reference grid (0) step^(y))=max (- 1,0)=0, the uflag of outdoor grid (1)^(x)=1.

Likewise, outdoor grid (2) only uflag adjacent in the x direction with reference grid (0) and outdoor grid (2)^(x) =0, then for outdoor grid (2):

Outdoor grid (2)

The step of outdoor grid (2)^(x)=max (the current step of outdoor grid (2)^(x), the step of reference grid (0)^(x)+ 1)=max (- 1,0+1)=1；The step of outdoor grid (2)^(y)=max (the current step of outdoor grid (2)^(y), reference grid (0) step^(y))=max (- 1,0)=0, the uflag of outdoor grid (2)^(x)=1.

Wherein, outdoor grid (3) only uflag adjacent in y-direction with reference grid (0) and outdoor grid (3)^(y)= 0, then for outdoor grid (3):

Outdoor grid (3)

The step of outdoor grid (3)^(y)=max (the current step of outdoor grid (3)^(y), the step of reference grid (0)^(y)+ 1)=max (- 1,0+1)=1；The step of outdoor grid (3)^(x)=max (the current step of outdoor grid (3)^(x), reference grid (0) step^(x))=max (- 1,0)=0, the uflag of outdoor grid (3)^(y)=1.

Further, it based on above counted intermediate parameters, calculates the outdoor grid that all state values are 0: updating The ray path parameter C for the outdoor grid that current all state values are 0_path: it can work as according to above-mentioned formula (2) to calculate Preceding state is the ray path parameter C of each of 0 outdoor grid_path。

Specifically, according to formula (2), as shown in figure 12, outdoor grid (1), outdoor grid (2) and outdoor grid at this time (3) state value is 0, calculates the ray path parameter C of these three grids_path:

Outdoor grid (1):

Outdoor grid (2):

Outdoor grid (3):

Step 705: grid outside the adjacent chamber of the reference grid in this iterative process is determined as next iteration process In reference grid, i.e., based on above-mentioned setting and update the process of state value, obtaining three flag by step 704 is 0 Outdoor grid, therefore the outdoor grid for being 0 using these three state values is as the reference grid of second of iterative process, i.e., second Reference grid in iterative process is grid outside the adjacent chamber of the reference grid in first time iterative process, continues above-mentioned change It include to empty benchmark set G again, then outdoor grid (1), outdoor grid (2) and outdoor grid (3) are put into base for operation In quasi- set G, outside the adjacent chamber of the outdoor grid (1) of acquisition, outdoor grid (2) and outdoor grid (3) the step of grid.Such as Figure 13 Shown, adjacent outdoor grid includes outdoor grid (4), room with outdoor grid (1), outdoor grid (2) and outdoor grid (3) Outer grid (5) and outdoor grid (6) continue according to the step 704 in above-mentioned calculating process, respectively the outer grid (1) of process chamber, Outdoor grid (2) and outdoor grid (3), that is, calculate all ginsengs of outdoor grid (4), outdoor grid (5) and outdoor grid (6) Number, calculated result are as shown in figure 14.Then according to step 704, outdoor grid (4), outdoor grid (5) and outdoor grid are calculated (6) ray path parameter C_path:

Outdoor grid (4):

Outdoor grid (5):

Outdoor grid (6):

It is as shown in figure 15 to update result.

Further, for example, after being updated by iteration several times, the reference grid in reference grid set G is outdoor grid Lattice (7) and outdoor grid (8).As shown in figure 16.

Then outdoor grid (10) is as the neighbours on the direction reference axis y of outdoor grid (8) and outdoor grid (10) uflag^(y)=0, the result that the parameter of outdoor grid (10) calculates in step 704 are as follows:

The step of outdoor grid (10)^(y)=max (the current step of outdoor grid (10)^(y), outdoor grid (8) step^(y)+ 1)=2

The step of outdoor grid (10)^(x)=max (the current step of outdoor grid (10)^(x), the step of reference grid (8)^(x))=max (- 1,2)=2, the uflag of outdoor grid (10)^(y)=1.

The ray path parameter C of outdoor grid (10) is calculated again_path:

Outdoor grid (10):

Calculated result is as shown in figure 17, it can be seen that calculated result has embodied building and blocked to signal.

The iteration meter that all outdoor grids in the objective emission antenna footprint can be carried out by above procedure All outdoor grids that can be traversed in the objective emission antenna footprint are calculated, until that cannot obtain again in reference grid set G Until getting new reference grid, i.e. set G completes for all outdoor grids empty, in the objective emission antenna footprint When calculating, iterative process terminates, and outdoor grids all at this time completes the calculating of ray path parameter, to obtain The ray path parameter C of the objective emission antenna currently calculated outdoor grid each of in its coverage area_path.Figure 18 is shown Final calculation result based on above-mentioned example.

By above-mentioned iterative process, be calculated each outdoor grid to objective emission antenna ray path parameter Afterwards, then it can further realize and determine the outdoor grid to target described in the step 602 in process as shown in FIG. 6 The path loss of transmitting antenna, i.e., it is true respectively according to the ray path parameter of each outdoor grid to the objective emission antenna Path loss of the fixed each outdoor grid to objective emission antenna.

Specifically, ray path parameter, the outdoor in step 602, according to an outdoor grid to objective emission antenna Grid to grid where objective emission antenna distance and the outdoor grid to objective emission antenna azimuth with have a down dip Angle, determine the outdoor grid to objective emission antenna path loss.

Wherein, it is sent out according to an outdoor grid to the ray path parameter of objective emission antenna, the outdoor grid to target Azimuth and angle of declination of the distance and the outdoor grid of grid to objective emission antenna, determine the outdoor where penetrating antenna Grid to objective emission antenna path loss, specifically can include:

Determine according to the following formula an outdoor grid to objective emission antenna path loss:

L=a+b × lg (d)+c × ln (C_path)+c_w×f(α,β)……………………(3)

Wherein, L indicates an outdoor grid to the path loss of objective emission antenna, and d is that the outdoor grid is sent out to target Penetrate the linear distance of antenna, C_pathFor the ray path parameter value of the outdoor grid, α is the outdoor grid to objective emission antenna Azimuth, β be the outdoor grid to objective emission antenna angle of declination, f (α, β) for objective emission antenna antenna waveform diagram Attenuation function, a, b, c, c_wFor constant coefficient.

Wherein, the C of all outdoor grids can be determined by above-mentioned iterative process_path；α, β and f (α, β) can be above The described public parameter according to antenna is accordingly and the position of grid and antenna obtains.

Specifically, for constant coefficient a, b, c, c_w, can be obtained by obtaining frequency sweep data by least square in training To optimal constant coefficient value, detailed process is described as follows:

Obtain the frequency sweep data of objective emission antenna；

For a frequency sweep data of objective emission antenna, according to the reception power of grid outdoor where the frequency sweep data and The transmission power of objective emission antenna determines the path loss of outdoor grid where the frequency sweep data；Determine the frequency sweep data institute In d, C of outdoor grid_path,α,β；

By the path loss and parameter of grid outdoor where identified every frequency sweep data, substitute into above-mentioned for determining one A outdoor grid to the path loss formula (3) of the objective emission antenna, be calculated constant coefficient a in above-mentioned formula, b、c、c_w。

Specifically, matrix A and column vector b can be initially set up, wherein the line number of A can be according to accessed frequency sweep number According to item number it is corresponding, the line number of A and the element number of b are identical；Wherein the columns of A be 4, be path loss Parameters in Formula (a, B, c, c_w) number.

For every frequency sweep data in the coverage area of objective emission antenna (being set as a), construct row vector [1, lg (d), ln(C_path), f (α, β)], wherein d, C_path, α and β be respectively the outdoor grid where a frequency sweep data path loss formula Parameter.Meanwhile pad value is calculated, the transmission power that pad value is specifically equal to antenna a subtracts the received power value of frequency sweep data. Row vector is put into A, pad value is put into and the position in corresponding b.

Such as: certain frequency sweep data belong to antenna a, and corresponding path loss parameters of formula is respectively as follows: d=395.2, ln (C_path)=5.71, α=53.5, β=10.2.The transmission power of antenna a is 39, and the reception power of frequency sweep data is -81.5, then Pad value is 120.5.Then, situation of the information of this frequency sweep data in A, b is as follows:

Wherein, line number of the data in A is identical as position of the pad value in b.All frequency sweeps of all antennas After data operate in above manner, final A and b is obtained, is then calculated: x=(A^T A)^-1 A^Tb

Then path loss equation coefficients are as follows: a=x (1), b=x (2), c=x (3), c_w=x (4).

For example pass through all constant parameter (a, b, c, the c according to the method described above of data provided by a certain common carrier_w) into Row training, the parameter trained are as shown in table 1:

Constant parameter example in 1 path loss formula (3) of table

Network	a	b	c	c_w
					LTE (F-band)	68.2	24	-2.4	0.48
LTE (D frequency range)	65.7	24	-2.9	0.48
					GSM900	58.8	20.6	-1.6	0.51
DCS1800	62.3	20.6	-2.2	0.51

For process as shown in FIG. 6, by determining that outdoor grid is damaged to the path of objective emission antenna in step 602 It, in step 603 can be according to the transmission power of the objective emission antenna and the outdoor grid to the mesh after consumption The path loss for marking transmitting antenna determines field strength of the objective emission antenna in the outdoor grid.For example, can according to Lower formula determines field strength of the objective emission antenna in the outdoor grid:

P=sendpower-L ... ... ... ... ... ... (4)

It can be seen that by the above process, for each transmitting antenna in target area, may be by above Path loss formula (3) and field strength calculation formula (4) calculate each transmitting antenna to grid outdoor each of in its coverage area The field strength of lattice.

For example, with LTE (for the F-band in (Long Term Evolution, long term evolution) network, if certain The transmission power sendpower of a antenna is 54dBm, the outdoor straight line of grid away from the antenna of certain in the antenna coverage away from From being 100 meters, ray path parameter is 0.0005, and horizontal sextant angle is 50 degree, and vertical angle is 8 degree, antenna waveform diagram pad value It is 13, i.e. sendpower=54dBm, d=100m, C_path=0.0005, f (α, β)=13 substitutes into formula above and obtains path damage Consume L=142.24, the field intensity value that antenna is generated in the outdoor grid are as follows: P=-88.24dBm.

Preferably, in order to avoid a large amount of useless calculating, for each transmitting antenna, for the 0th layer of all outdoor grid Lattice all calculate field strength positioned at all outdoor grids of ground floor, and for the 1st layer and with the outdoor grid in upper layer, only It needs to calculate the outdoor grid field strength for being located at and making a circle in building week, that is, calculates the field of the outdoor grid adjacent with building grid By force.

Further, communication network field strength distribution provided by one embodiment of the present of invention determines method by as schemed Process shown in 6 determines that objective emission antenna, can also be according to identified objective emission after the field strength in outdoor grid Field strength of the antenna in each outdoor grid, determines indoor grille of the objective emission antenna in target area according to following formula Interior field strength a, wherein indoor grille refers to the grid for falling into interior of building:

Wherein, P_inFor field strength of the objective emission antenna in an indoor grille, P_jTo be built with where the indoor grille The field strength of grid, η outside jth room in all J adjacent outdoor grids of object₀Indicate that electromagnetic wave once penetrates the decaying number of wall Value, such as concrete wall, average attenuation η₀It is 13；η₁Indicate that unit distance stops body attenuation value, such as usually can be with Take 0.7；d_jIndicate will the distance between grid outside the indoor grille and jth room,It indicates to all J being calculated A value chooses maximum value, and j, J are the integer more than or equal to 0.

Specifically, shown in Figure 19 for two-dimensional scene, wherein Figure 19 Oxford gray grid is indoor grille.For The field strength that objective emission antenna generates on the indoor grille of Dark grey is found out (set institute with the indoor grille to be calculated first The indoor grille to be calculated is identical for the geographical height in i) place, and all outdoors adjacent with building where indoor grille i Grid, for Figure 19, all oblique line grids are all outdoor grids adjacent with the indoor grille of grey.If a shared K A outdoor grid, the outdoor grid of each of this K outdoor grid determine that field strength is respectively as follows: by process as shown in FIG. 6 P_jThe field strength of the outdoor grid of each of this K determined outdoor grid is substituted into formula (5) by (j=1 ..., K) To obtain the field strength of indoor grille i.

Figure 20 shows the example that an indoor field strength is calculated based on outdoor grid field strength, wherein Figure 20 Oxford gray Grid is indoor grille, identical as Figure 19, it is assumed that some transmitting antenna in adjacent with Dark grey grid all outdoor grids, Field strength only is generated in grayish three grids, field intensity value is as shown in the figure.

It is then as follows to the calculating process of the field intensity value of an indoor grille shown in dotted line frame according to formula (5):

Utilize the field strength for the indoor grille that the outdoor field strength that field strength is -89.5 is calculated are as follows:

P₁=-89.5- (13+0.7 (2 × 5))

Utilize the field strength for the indoor grille that the outdoor field strength that field strength is -93.2 is calculated are as follows:

Utilize the field strength for the indoor grille that the outdoor field strength that field strength is -95.7 is calculated are as follows:

P₃=-95.7- (13+0.7 (2 × 5))

Therefore, the field strength of the indoor grille are as follows:

P=max { P₁,P₂,P₃}=- 109.5dBm

Therefore, a kind of communication network provided by the embodiment of the present invention determines that method can not only determine outdoor grid Field strength, can also field strength further based on the outdoor grid determined, determine the field strength of indoor grille.It can be seen that By above-mentioned process, for each transmitting antenna, outdoor grid all in (for example being 3km) and room in coverage area Field strength on outer grid can be determined, correspondingly, each grid, including indoor grille and outdoor grid, by surrounding All transmitting antennas influence of coverage area to the grid can also be determined, it is possible thereby to the field strength distribution of being determined Target area establish the fingerprint base for reflecting each grid Cover Characteristics, for example can receive with each grid most strong Top N antenna antenna id and field intensity value be grid fingerprint characteristic, table 2 shows the Top7 field strength example of three grids, Wherein, antenna id (identification) can be by the public parameter of antenna according to acquisition.

2 grid characteristic fingerprint library example of table

Antenna id	80031-2	72397-3	72397-1	72396-3	72396-1	79272-3	72396-2
								Field strength	-94.5	-81.6	-87.1	-95.1	-100.1	-104.9	-105.1
Antenna id	72397-2	72397-3	72397-1	72396-3	72396-1	79272-3	72396-2
								Field strength	-76.6	-81.6	-87.1	-95.1	-100.1	-104.9	-105.1
Antenna id	88159-1	88159-3	88159-2	91520-2	72788-3	72788-2	88259-1
								Field strength	-81.2	-86.5	-92.1	-95.8	-97.5	-97.9	-98.5

A kind of provided communication network field strength distribution determines the obtained target area of method through the embodiment of the present invention In each grid Cover Characteristics fingerprint base, can further be used for includes covering field strength analysis, weak overlay area is determining, surveys Amount report positioning, user location determination etc..

By above description as can be seen that a kind of communication network field strength distribution determination side provided in embodiments of the present invention In method, the outdoor grid in target area is obtained first, further according to the objective emission antenna in outdoor grid and target area The spatial position of place grid, the path loss of determining outdoor grid to objective emission antenna, finally according to objective emission antenna Transmission power and outdoor grid to the path loss of objective emission antenna, determine objective emission antenna in outdoor grid Field strength.It can be seen that being capable of providing a kind of field-strength prediction model in the embodiment of the present invention, which is based on to target area Rasterizing, the field strength of outdoor grid is calculated and is iterated step by step using grid outdoor where antenna as starting point grid first, I.e. it can be appreciated that a kind of ray tracing of rasterizing, to calculate the path loss of each outdoor grid, and then determines The field strength of each outdoor grid out, at the same provide it is a kind of based on the outdoor grid field strength determined by the above method into one The method for determining indoor field strength is walked, so that the field strength on all grids in target area can be predicted out, it can accurate The field strength distribution situation for determining target area.

Therefore, compared with prior art, a kind of communication network field strength distribution provided in an embodiment of the present invention determines that method is examined Building distribution and the influence of topography and geomorphology are considered, higher calculated result accuracy can be reached compared with empirical model, and can be with Carry out three-dimensional field intensity prediction；Also, a kind of communication network field strength distribution provided in an embodiment of the present invention determines method to application Scene is not strict with, and most of application scenarios can be suitable for, and relatively half determines that model can have wider array of applicability；This Outside, calculating cost needed for a kind of communication network field strength distribution provided in an embodiment of the present invention determines method in feasible region, And building and terrain information are only needed substantially accurately, required precision does not need stringent as determining model.Finally, this A kind of communication network field strength distribution that inventive embodiments provide determines that method is not only easily transplanted to 2G network, 3G network, 4G Network, the even more advanced network of the 5G network for future can adaptability application, can be subsequent network Development It lays the foundation.

Based on the same technical idea, one embodiment of the present of invention also provides a kind of determining dress of communication network field strength distribution It sets, which can be performed above-mentioned communication network field strength distribution and determine embodiment of the method.A kind of communication provided in an embodiment of the present invention Network field strength distribution determining device is as shown in figure 21, comprising:

Module 2101 is obtained, for obtaining the outdoor grid in target area, the outdoor grid, which refers to, does not fall within building Grid inside object, a grid are the area of space being sized；

First determining module 2102, for according to the objective emission antenna in the outdoor grid and the target area The spatial position of place grid, determine the outdoor grid to the objective emission antenna path loss；

Second determining module 2103, for being arrived according to the transmission power and the outdoor grid of the objective emission antenna The path loss of the objective emission antenna determines field strength of the objective emission antenna in the outdoor grid.

Wherein, the first determining module 2102, is specifically used for: by n times iterative process, each outdoor grid is calculated and arrives The ray path parameter of the objective emission antenna, N are the integer more than or equal to 1；

Wherein, include: in each iterative process

Specifically, the first determining module 2102, is used for: to grid outside an adjacent chamber, obtaining and grid outside the adjacent chamber The ray path parameter of adjacent reference grid in the every one-dimensional square in space, and to obtain outside the adjacent chamber grid relatively described Spatial deviation amount of the outdoor grid in the every one-dimensional square in space where objective emission antenna；According in description above Provided in formula (1) calculate an adjacent chamber outside grid ray path parameter.Wherein, formula (1) can be according to iteration The characteristics of operation, by by target area average grid, and assuming for indicating an outdoor grid in preprocessing process The parameter of spatial deviation amount of the outdoor grid on the i-th dimension direction of space where the relatively described objective emission antenna is with one Step is linear module, and a step is size of the grid in the dimension, is reduced to formula (2).

Wherein, the first determining module 2102, is specifically used for: according to an outdoor grid to the ray road of objective emission antenna The distance of grid where diameter parameter, the outdoor grid to objective emission antenna and the outdoor grid are to objective emission antenna Azimuth and angle of declination, determine the outdoor grid to objective emission antenna path loss.

Specifically, the first determining module 2102, is used for: determining one according to provided formula (3) in the above description Path loss of a outdoor grid to the objective emission antenna.

Wherein, the second determining module 2103, is specifically used for: being determined according to provided formula (4) in the above description Field strength of the objective emission antenna in the outdoor grid.

Preferably, the communication network field strength distribution determining device that one embodiment of the present of invention provides can also include: the Three determining modules 2104, for the field strength according to the identified objective emission antenna in each outdoor grid, according to Formula provided in description above (5) determines field of the objective emission antenna in the indoor grille in target area By force, an indoor grille refers to the grid for falling into interior of building.

The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.

These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.

These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.

Although preferred embodiments of the present invention have been described, it is created once a person skilled in the art knows basic Property concept, then additional changes and modifications may be made to these embodiments.So it includes excellent that the following claims are intended to be interpreted as It selects embodiment and falls into all change and modification of the scope of the invention.

Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies Within, then the present invention is also intended to include these modifications and variations.

Claims

1. a kind of communication network field strength distribution determines method characterized by comprising

The outdoor grid in target area is obtained, the outdoor grid refers to the grid for not falling within interior of building, a grid For the area of space being sized；

By n times iterative process, ray path parameter of each outdoor grid to the objective emission antenna, and root is calculated According to the ray path parameter of each outdoor grid to the objective emission antenna, determine each outdoor grid described in respectively The path loss of objective emission antenna, N are the integer more than or equal to 1；

It is damaged according to the transmission power of the objective emission antenna and the outdoor grid to the path of the objective emission antenna Consumption determines field strength of the objective emission antenna in the outdoor grid；

Wherein, include: in each iterative process

Determine the reference grid during current iteration；Wherein, the reference grid in first time iterative process is target hair The outdoor grid where antenna is penetrated, the reference grid for the second time and in later iterative process is the base in a preceding iterative process Grid outside the adjacent chamber of quasi- grid；

Outdoor where the relatively described objective emission antenna of grid outside the ray path parameter of reference grid and adjacent chamber The spatial deviation amount of grid calculates the ray path parameter of grid outside adjacent chamber；Wherein, the base in first time iterative process The ray path parameter of quasi- grid is setting value.

2. the method as described in claim 1, which is characterized in that according to the ray path parameter of reference grid and adjacent chamber The spatial deviation amount of outdoor grid where the relatively described objective emission antenna of outer grid calculates grid outside adjacent chamber and penetrates Thread path parameter, comprising:

To grid outside an adjacent chamber, the reference grid adjacent in the every one-dimensional square in space with grid outside the adjacent chamber is obtained Ray path parameter, and the outdoor grid where obtaining the relatively described objective emission antenna of grid outside the adjacent chamber are every in space Spatial deviation amount in one-dimensional square；

Wherein, C_pathIndicate the ray path parameter of grid outside an adjacent chamber being calculated,It indicates outside the adjacent chamber The ray path parameter of grid reference grid adjacent on the i-th dimension direction of space, step⁽ⁱ⁾Indicate grid phase outside the adjacent chamber To spatial deviation amount of the outdoor grid where the objective emission antenna on the i-th dimension direction of space, D is space dimension Number, 1≤i≤D, i, D are positive integer.

3. the method as described in claim 1, which is characterized in that according to outdoor grid penetrating to the objective emission antenna Thread path parameter, determine the outdoor grid to the objective emission antenna path loss, comprising:

Ray path parameter, the outdoor grid to the objective emission according to an outdoor grid to the objective emission antenna Azimuth and angle of declination of the distance and the outdoor grid of grid to the objective emission antenna, determine the room where antenna Path loss of the outer grid to the objective emission antenna.

4. method as claimed in claim 3, which is characterized in that according to outdoor grid penetrating to the objective emission antenna The distance of grid where thread path parameter, the outdoor grid to the objective emission antenna and the outdoor grid are to the mesh Mark transmitting antenna azimuth and angle of declination, determine the outdoor grid to the objective emission antenna path loss, comprising:

L=a+b × lg (d)+c × ln (C_path)+c_w×f(α,β)

Wherein, L indicates an outdoor grid to the path loss of the objective emission antenna, and d is the outdoor grid to the mesh Mark the linear distance of transmitting antenna, C_pathFor the ray path parameter value of the outdoor grid, α is the outdoor grid to the target The azimuth of transmitting antenna, β are angle of declination of the outdoor grid to the objective emission antenna, and f (α, β) is the objective emission The antenna waveform diagram attenuation function of antenna, a, b, c, c_wFor constant coefficient.

5. method according to any one of claims 1 to 4, which is characterized in that according to the transmitting of the objective emission antenna Power and the outdoor grid determine the objective emission antenna in the room to the path loss of the objective emission antenna Field strength in outer grid, comprising:

P=sendpower-L

Wherein, P is field strength of the objective emission antenna in an outdoor grid, and sendpower is the transmitting of objective emission antenna Power, L are path loss of the objective emission antenna to the outdoor grid.

6. method according to any one of claims 1 to 4, which is characterized in that determine the objective emission antenna described After field strength in outdoor grid, further includes:

According to field strength of the identified objective emission antenna in each outdoor grid, the mesh is determined according to following formula Field strength of the transmitting antenna in the indoor grille in target area is marked, an indoor grille refers to the grid for falling into interior of building Lattice:

Wherein, P_inFor field strength of the objective emission antenna in an indoor grille, P_jFor building phase where with the indoor grille The field strength of grid, η outside jth room in the outdoor grid of adjacent all J₀Indicate that electromagnetic wave once penetrates the attenuation value of wall, η₁ Indicate that unit distance stops body attenuation value, d_jIndicate will the distance between grid outside the indoor grille and jth room, It indicates to choose all J values being calculated maximum value, j, J are the integer more than or equal to 0.

7. a kind of communication network field strength distribution determining device characterized by comprising

Module is obtained, for obtaining the outdoor grid in target area, the outdoor grid, which refers to, does not fall within interior of building Grid, a grid are the area of space being sized；

First determining module, for each outdoor grid to be calculated to the objective emission antenna by n times iterative process Ray path parameter, and according to the ray path parameter of each outdoor grid to the objective emission antenna, it determines respectively each For the outdoor grid to the path loss of the objective emission antenna, N is the integer more than or equal to 1；

Second determining module, for according to the transmission power of the objective emission antenna and the outdoor grid to the target The path loss of transmitting antenna determines field strength of the objective emission antenna in the outdoor grid；

Wherein, include: in each iterative process

8. device as claimed in claim 7, which is characterized in that first determining module is specifically used for:

9. device as claimed in claim 7, which is characterized in that first determining module is specifically used for:

10. device as claimed in claim 9, which is characterized in that first determining module is specifically used for:

L=a+b × lg (d)+c × ln (C_path)+c_w×f(α,β)

11. the device as described in any one of claim 7 to 10, which is characterized in that second determining module is specific to use In:

P=sendpower-L

12. the device as described in any one of claim 7 to 10, which is characterized in that further include:

Third determining module, for the field strength according to the identified objective emission antenna in each outdoor grid, according to Following formula determines field strength of the objective emission antenna in the indoor grille in target area, and an indoor grille, which refers to, to be fallen Enter the grid of interior of building: