CN111818555B - Radio monitoring station coverage area evaluation and analysis method based on virtual station building - Google Patents
Radio monitoring station coverage area evaluation and analysis method based on virtual station building Download PDFInfo
- Publication number
- CN111818555B CN111818555B CN202010646709.4A CN202010646709A CN111818555B CN 111818555 B CN111818555 B CN 111818555B CN 202010646709 A CN202010646709 A CN 202010646709A CN 111818555 B CN111818555 B CN 111818555B
- Authority
- CN
- China
- Prior art keywords
- monitoring
- coverage
- data
- station
- module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 209
- 238000004458 analytical method Methods 0.000 title claims description 28
- 238000011156 evaluation Methods 0.000 title claims description 24
- 238000000034 method Methods 0.000 claims abstract description 50
- 238000012360 testing method Methods 0.000 claims description 83
- 238000004364 calculation method Methods 0.000 claims description 34
- 238000007726 management method Methods 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims description 24
- 230000035945 sensitivity Effects 0.000 claims description 23
- 238000013439 planning Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 19
- 238000013213 extrapolation Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 9
- 230000006855 networking Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000012821 model calculation Methods 0.000 claims description 8
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims description 7
- 238000012986 modification Methods 0.000 claims description 7
- 230000004048 modification Effects 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000013523 data management Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 238000013144 data compression Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 238000004891 communication Methods 0.000 abstract description 9
- 230000008569 process Effects 0.000 description 7
- 230000001360 synchronised effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002050 diffraction method Methods 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/02—Arrangements for optimising operational condition
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/06—Testing, supervising or monitoring using simulated traffic
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention relates to the technical field of wireless communication, in particular to a method for evaluating and analyzing the coverage area of a radio monitoring station based on virtual station building.
Description
The application is named as: a radio monitoring station coverage assessment analysis system and method, application date are: 30 days 6 and 7 in 2017, and the application numbers are as follows: 201710519077.3 patent application for inventions.
Technical Field
The invention relates to the technical field of wireless communication, in particular to a radio monitoring station coverage area evaluation and analysis method based on virtual station building.
Background
The radio monitoring station is a basic technical means of radio management work, provides technical support in the aspects of planning assignment of frequency, electromagnetic environment monitoring, radio interference investigation and the like, and determines key factors of the construction scale of the monitoring station by the monitoring coverage range. Meanwhile, related planning files of the industrial and informatization department clearly require monitoring coverage capacity of various radio monitoring sites, so a set of scientific and objective monitoring coverage condition prediction and evaluation method system needs to be established to ensure the benign development of radio monitoring management work, plan and perfect the construction acceptance standards of various radio monitoring sites and unify a monitoring coverage evaluation system method, thereby improving the coverage efficiency of a radio monitoring network.
Factors affecting the coverage of a radio monitoring station: system parameters, radio wave propagation loss and radio transmission source parameters are monitored. The monitoring system parameters mainly comprise fixed parameter monitoring system sensitivity, variable parameters, antenna height, antenna gain and the like; the radio transmission source parameters affecting the coverage area of the monitoring station include the transmission power of the transmission source, the antenna gain, the antenna height and the like. The estimated monitoring range can be regarded as an inverse process of estimating the coverage radius of radio transmission, the monitoring range of the radio fixed monitoring station is predicted, and the propagation loss of radio waves needs to be calculated. At present, electromagnetic wave propagation models are roughly divided into three types, namely theoretical models based on analytical equations, numerical models based on numerical calculations, and statistical models based on measured data. The theoretical model is generally a model established according to a classical electromagnetic wave propagation equation Maxwell equation set or classical geometric propagation reflection, diffuse reflection, diffraction and the like, describes the basic physical law of electromagnetic wave propagation, and has a certain difference from practical application. The numerical model based on numerical calculation is to carry out numerical solution on the analytical model, so that the computer processing is convenient, and the result has no substantive difference compared with the analytical model. The statistical model based on the measured data is an electromagnetic wave propagation model which is established by giving a calculation formula of propagation characteristics through statistical analysis and processing on the basis of actual measurement of a large number of electromagnetic signals. However, wireless propagation environments are different, and the coverage is predicted and evaluated by a model alone, so that the effect and accuracy are relatively poor, so that a standard transmitting system needs to be introduced to carry out drive test, the model is corrected by combining drive test data, and finally, an actual wireless propagation environment model which accords with the current area is obtained, and a prepared prediction result can be obtained.
At present, the coverage area evaluation work of the radio monitoring station in China is less, the acceptance modes of the newly-built station are not uniform, and the coverage area results are different due to the fact that differences among a radio wave propagation prediction model, a road test scheme and final data correction are large.
Disclosure of Invention
The invention aims to provide a radio monitoring station coverage area evaluation and analysis method based on virtual station building, which solves the problem of monitoring coverage area analysis of a radio monitoring station, is used for coverage evaluation of a built station or a proposed station, and is beneficial to station maintenance and scientific station building based on a data mode.
The invention is realized by the following steps:
a radio monitoring site coverage assessment analysis system, comprising:
the coverage analysis subsystem is used for acquiring monitoring system parameters, radio wave propagation loss and radio emission source parameters, and modeling and analyzing according to the acquired parameters;
the GIS geographic information management subsystem comprises a two-dimensional/three-dimensional electronic map module and a map information calling module, and is used for establishing a propagation module, planning a road test route, determining a transmitting point, displaying field intensity area distribution, managing a monitoring station, managing a radio station, inquiring a function and managing GPS track data;
the standard transmitting subsystem comprises a signal source, a power amplifier, an omnidirectional transmitting antenna and a GPS positioning device, and the standard transmitting subsystem generates EIRP power of 1W/3W/5W/10W power grade by covering any frequency point signal of 20MHz-6GHz and amplifying the signal by a power amplifier so as to provide a standard transmitting source for the road test subsystem;
the system comprises a road test subsystem, a monitoring subsystem and a monitoring subsystem, wherein the road test subsystem comprises a standard transmitting module, a radio detection module, an RMTP monitoring network protocol module, a special monitoring transmission network and a GIS three-dimensional electronic map module, and acquires the signal field intensity of the standard transmitting subsystem through the monitoring function of a radio monitoring system, and finally obtains the actual coverage information of a monitoring station by combining the map information and a GPS track;
the station management subsystem comprises a radio monitoring station management module, a radio station management module, a station database calling module and a virtual station management module, and is used for rapidly checking and removing the position of a legal station when interference occurs and reducing the time for determining the interference position;
and the data management subsystem is used for data storage, data calling, data compression, data encryption and data attribute definition.
Further, in the most preferred embodiment of the present invention, the coverage analysis subsystem further includes an ITU-r.p1546 propagation model calculation module, a virtual station building function module, an isoline field strength prediction calculation module, a field strength coverage calculation function module, a coverage model modification module based on drive test data, a monitoring networking area coverage analysis module, and a typical service management module;
the ITU-R.P1546 propagation model calculation module is used for providing a field intensity prediction proposal of 30MHz-3000MHz frequency band ground service, the proposal comprises a correction formula for carrying out interpolation/extrapolation on field intensity under different transmitting power, frequency, antenna height and transmission distance, and any field intensity value of transmitting frequency in the range of 30MHz-3000MHz and the propagation distance of 1 km-1000 km, which cannot be directly read from a provided curve, is obtained by using the correction formula;
the virtual station building functional module is used for realizing automatic calculation of any frequency by various radio wave propagation models and three-dimensional electronic maps and drawing field intensity coverage situation maps under different power levels;
the isoline field intensity prediction calculation module is used for automatically calculating and drawing isolines of different service frequencies under different transmitting power levels to form a boundary area;
the field intensity coverage calculation function module is used for displaying a field intensity coverage map drawn by the model in a map by combining a two-dimensional electronic map, and realizing the functions of monitoring coverage capability display of a monitoring network of a specific radio station, calculating the coverage area and the area coverage rate of a monitoring station/monitoring network, displaying or hiding the coverage map of the monitoring station at will, increasing or deleting the number of the monitoring stations at the same time, automatically comparing the coverage maps under different parameters, calling a local station database, checking the monitoring condition of a certain type of radio station, counting the monitoring coverage rate of various stations, establishing a key monitoring area by using the map and calculating the fixed point field intensity of any position on the map;
the coverage model modification module based on the drive test data is used for optimizing by adopting a standard macro cell model commonly applied in network planning;
the monitoring networking area coverage analysis module obtains the networking coverage condition of the current whole area by utilizing the coverage analysis map of each radio monitoring site, and obtains the monitoring site corresponding to the maximum field intensity of the area by utilizing the radio monitoring site management function of the two-dimensional electronic map to obtain the coverage depth index of the monitoring network;
the typical service management module and other functional modules are used for a key monitoring application to provide rapid monitoring coverage statistics.
A radio monitoring site coverage assessment analysis method comprises the following steps:
s1: marking sites needing simulation of building the sites on the loaded three-dimensional electronic map, wherein the sites comprise site names, site heights, site position information, antenna types, antenna gains and nominal monitoring sensitivity;
s2: setting site service parameters to be simulated, wherein the site service parameters comprise service frequency bands/points, transmitting end height information, geographic information, EIRP power levels and transmitting bandwidths;
s3: drawing rays in a three-dimensional map by taking a monitoring station as a center according to a preset angle, and virtually establishing a plurality of emission sources on each ray;
s4: calculating the monitoring value of the emission source at the monitoring station position by using an ITU-R.P1546 algorithm, and continuously calculating the outside on the ray;
s5: when the calculated value is smaller than the monitoring receiving sensitivity, finding the end point of the sensitivity contour line by combining the last value, namely the ray calculation is completed;
s6: repeating S3 to S6 to complete contour points of all rays and connecting all the points to form a power contour curve;
s7: adjusting other service parameters, wherein the other service parameters comprise frequency and EIRP power, and drawing other power contour curves under different service conditions;
s8: other ray subdivision points in each power contour line are connected to form a plurality of boundary areas, and the field intensity is distinguished through color gradient to form a field intensity coverage map;
s9: and displaying the simulated field intensity coverage range of the monitoring station based on the two-dimensional electronic map.
Further, in the most preferred embodiment of the present invention, when the field intensity coverage evaluation is performed by using actual drive test data, the function of the drive test system is used, which specifically includes the following steps:
s101: preparation before drive test: the method comprises the steps of planning and setting parameters of a transmitting subsystem, a GPS antenna, a network transmission tool, an RMTP monitoring network protocol and monitoring site information;
s102: planning a driving route: selecting at least one main road according to a three-dimensional electronic map, formulating a driving speed and a driving distance according to the Lee's law, carrying out reverse thrust according to the EIRP power of a transmitting system and the sensitivity of a monitoring system, and determining a sensitivity power contour line area;
s103: actual measurement data acquisition: starting a transmitting system on the road test vehicle, starting a GPS positioning device, traveling on a planned road according to a preset vehicle speed, and simultaneously controlling the transmitting system and the monitoring system to record monitoring data, transmitting time and acquiring time according to a certain time interval until all roads are acquired;
s104: and synchronously finishing data acquisition: selecting and controlling a transmitting system to transmit signal requirements of different services, completing data required by a plurality of services by one drive test, and independently repeating the drive test of other services;
s105: and (3) drive test data processing: removing repeated data in the drive test data, removing data without latitude and longitude or drift of latitude and longitude, combining the drive test data and GPS data and importing the data into a system to form a data format which can be displayed on a map;
s106: and (3) actual measurement data display: the processed data can be displayed on a map according to a preset rule, and the data of each road measuring point is drawn on the map;
s107: the method comprises the steps of utilizing measured data to draw a coverage range of a monitoring station, setting a threshold, deducing a plurality of points which are close to or equal to receiving sensitivity according to the drive test data, connecting the points, drawing a measured coverage range, and comparing the measured coverage range with a calculated coverage range.
Further, in the most preferred embodiment of the present invention, if the prediction calculation model of ITU-r.p1546 needs to be corrected based on the drive test data, the method for correcting the prediction model is used, which specifically includes the following steps:
s201: analyzing the GPS data and the monitoring station data to obtain a received signal average field strength value Prx of the monitoring station at each moment, effective radiation power Ptx of a transmitter, a distance d between a transmitting system and the monitoring station, a height Hms on the ground where the transmitting system is located, an effective height Heff of an antenna of the monitoring station, diffraction Loss diffn and a ground object Loss factor Clutter _ Loss;
s202, reading the test data in the step S201: the average field strength value Prx of the received signal, the effective radiation power Ptx of the transmitter, the distance d between the transmitting system and the monitoring station, the height Hms of the transmitting system on the ground, the effective height Heff of the antenna of the monitoring station, the diffraction Loss diffn and the ground feature Loss factor Clutter _ Loss, and the correction calculation is carried out.
Further, in a most preferred embodiment of the present invention, the method for predicting ITU-r.p1546-5 field strength comprises:
s301: determining the first propagation type, the second propagation type and the third propagation type as the types of terrestrial, cold-sea or warm-sea propagation paths, respectively, and if the propagation types are mixed paths, determining a combination of two path types that can be considered as the first propagation type and the second propagation type, and if the paths can be represented by a single type, considering it as the first propagation type;
s302: determining two nominal time percentages;
s303: determining two nominal frequencies;
s304: determining nominal distances of a low end and a high end;
s305: for the first transmission type, proceed as S306 to S311;
s306: for the low end nominal percentage of time, proceed as S307 to S310;
s307: for the low end nominal frequency, proceed as S308 and S309;
s308: for the required distance d and the transmitting/base station antenna height h1, a field strength exceeding 50% of the location point is obtained for the receiving/mobile station antenna at the representative scattering height R above the earth:
s309: if the desired frequency does not match the low-end nominal frequency, repeating S308 for the high-frequency nominal frequency and interpolating or extrapolating the two field strengths;
s310: if the percentage of time required does not match the lower nominal percentage of time, repeating S307 to S309 for the higher nominal percentage of time and interpolating the two field strengths;
s311: if the hybrid path is predicted, the path requirements for each propagation type are from S306 to S310;
s312: correcting the field strength of the terrain clearance angle at the receiver/mobile unit if information in terms of the terrain clearance angle at the receiver/mobile station antenna adjacent the land is available;
s313: calculating an estimated value of field intensity caused by scattering, and taking the maximum value of E and Ets;
s314: correcting the field strength of the antenna height h2 of the receiving/moving station;
s315: if clutter exists around the transmitting/base station terminal, correcting the effect;
s316: correcting the inclined diameter;
s317: if a path less than 1km is given, S301 to S316 for d ═ 1km need to be followed;
s318: if the field strength of the position percentage of the receiving/mobile station antenna adjacent to the land is required to be not 50 percent, correcting the field strength of the required position percentage;
s319: limiting the obtained field strength to a preset maximum value, if the field strength calculation of the hybrid path has been made for less than 50% of the time percentage, the maximum field strength value needs to be calculated by linear interpolation between all terrestrial field strength values and all sea field strength values, which is given by the following formula:
Emax＝Efs+dsEse/dtotaldB(V/m)
s320: the field strength is converted into an equivalent fundamental transmission loss of the path.
Further, in the most preferred embodiment of the present invention, in S302, the method comprises:
the percentage of time required is >1 and <10, where the nominal percentage of time for the low and high ends is 1% and 10%, respectively;
the percentage of time required is >10 and <50, where the nominal percentage of time for the low and high ends is 10% and 50%, respectively;
if the percentage of time required is equal to 1%, 10%, or 50%, then the value may be considered the low end nominal percentage of time without requiring the interpolation process of S310.
Further, in the most preferred embodiment of the present invention, in S303, the method further includes:
when the required frequency is less than 600MHz, the low-end frequency and the high-end frequency are respectively 100MHz and 600 MHz;
when the required frequency is less than 600MHz, the low-end frequency and the high-end frequency are 600MHz and 2000MHz respectively;
if the desired frequency is equal to 100 or 600 or 2000MHz, this value may be considered the low end nominal frequency without the need for interpolation/extrapolation processing in S309.
Further, in the most preferred embodiment of the present invention, in S308, the method comprises:
for a transmit/base station antenna height h1 of less than 10m, the field strength at the required height and distance is determined.
Further, in the most preferred embodiment of the present invention, it specifically further comprises:
s381: determining the values of the low and high nominal heights h1, which should be considered as the low nominal height value h1 if the value of h1 matches any of the nominal values 10, 20, 37.5, 75, 150, 300, 600 or 1200m, without interpolation;
s382: for the low-end nominal value h1, proceeding as S383 to S385;
s383: for the low end nominal distance value, proceed as per S384;
s384: for the required distance d and the transmitting/base station antenna height h1, a field strength exceeding 50% of the location point is obtained for the receiving/mobile station antenna at the representative spread height R;
s385: if the desired distance does not match the lower nominal distance, repeating S384 for the upper nominal distance and interpolating the two field strengths over that distance;
s386: if the desired transmit/base station antenna height h1 does not meet any of the nominal values, S383 through S385 are repeated and h1 is interpolated/extrapolated.
The invention has the beneficial effects that:
(1) the invention can carry out coverage evaluation on the established site or the proposed site so as to realize site maintenance and scientific site establishment based on a data mode;
(2) the method can provide typical barrier heights of different terrains, provide a more accurate correction method for the height of the receiving antenna, and obtain the radio wave propagation loss under any path by utilizing interpolation and extrapolation methods;
(3) the system has EIRP power output meeting multiple power levels in a 60GHz frequency range, and also has a GPS function, so that any output information can be ensured to be marked with GPS time;
(4) the electromagnetic spectrum situation display can be reflected more truly based on the radio monitoring application of the three-dimensional electronic map;
(5) the system can synchronously transmit data of a subsystem and a monitoring system based on GPS time information, supports synchronous and asynchronous data processing, and can acquire field intensity regions at various levels by using a mathematical model.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 is a block diagram of a system for evaluating and analyzing the coverage of a radio monitoring station according to the present invention;
fig. 2 is a flow chart of the radio monitoring station coverage area virtual station building of the present invention;
FIG. 3 is a flow chart of the present invention for synchronous data acquisition;
FIG. 4 is a flow chart of asynchronous data acquisition in accordance with the present invention;
FIG. 5 is a flow chart of the model optimization based on drive test data according to the present invention.
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict.
It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention, and the components related to the present invention are only shown in the drawings rather than drawn according to the number, shape and size of the components in actual implementation, and the type, quantity and proportion of the components in actual implementation may be changed freely, and the layout of the components may be more complicated.
In an embodiment, please refer to fig. 1, a system for evaluating and analyzing coverage of a radio monitoring station includes: the system comprises a coverage analysis subsystem 10, a GIS geographic information management subsystem 20, a standard transmitting subsystem 30, a drive test subsystem 40, a station management subsystem 50 and a data management subsystem 60.
Specifically, the coverage analysis subsystem 10 includes: the system comprises an ITU-R.P1546 propagation model calculation module, a virtual station building function module, an isoline field intensity prediction calculation module, a field intensity coverage calculation function module, a coverage model correction module based on drive test data, a monitoring networking area coverage analysis module, a typical service management module and other function modules.
The ITU-R.P1546 propagation model calculation module comprises an ITU-R P.1546 recommendation which gives a field intensity prediction method of 30MHz-3000MHz frequency band ground service, and the recommendation gives radio wave propagation curves of a transmission source with effective transmission power of 1kW at 100MHz, 600MHz and 2000MHz frequency points under the conditions of certain transmitting antenna heights, receiving antenna heights, propagation path types, time probabilities and place probabilities on the basis of measured data, and 24 curves are calculated. The proposal also provides a correction formula for carrying out interpolation/extrapolation on the field intensity under different transmitting powers, frequencies, antenna heights and transmission distances, and any field intensity value with the transmitting frequency of 30MHz-3000MHz and the propagation distance of 1 km-1000 km, which can not be directly read from the provided curve, can be obtained by using the correction formula. The additional information and method applied by the prediction method mainly comprise the maximum field intensity value, the height h1 of the transmitting antenna is determined, the field intensity is calculated on the basis of an interpolation method on a distance axis, the field intensity is calculated on the basis of an interpolation/extrapolation method on a frequency domain, the field intensity is calculated on the basis of a basic interpolation method on a time domain, the field intensity of a mixed path is calculated, the heights of the reference receiving antenna and other height revisions, the ground clearance angle revisions, the distribution of the predicted points of the land area, the suggested operation steps and the like;
the virtual station building function module mainly aims at automatically calculating any frequency based on various radio wave propagation models and three-dimensional electronic maps and drawing field intensity coverage situation maps under different power levels, and can realize basic evaluation of a radio monitoring station in the early stage of building the station by utilizing the function and guide actual station building. The established virtual station can realize the same coverage analysis of all established stations, and any new monitoring station can be planned, established, deleted and modified on a map at will by utilizing position information or in a mouse click mode; the established virtual station can use the equipment information of a certain existing monitoring station as a template, and can calculate the coverage of the virtual station on a certain radio station by using a propagation model; circular station building simulation is supported, and a proper number of stations are automatically calculated for filling; the method supports polygon station building simulation, can realize the station building simulation of the coverage of any polygon area, and can count the area coverage rate, station coverage rate and coverage rate to a heavy spot area of the whole monitoring network after a new monitoring station is added; for coverage maps generated for virtual stations, the virtual stations may be moved after being mapped, helping to determine the best established position of the virtual stations;
the isoline field intensity prediction calculation module can automatically calculate and draw isolines of different service frequencies under different transmitting power levels to form a boundary area, and the specific operation is as follows: drawing a plurality of rays around the radio monitoring station according to angles by taking the radio monitoring station as a center; secondly, taking out one ray, establishing a transmission source of a typical service at a fixed distance from a monitoring station, calculating a monitoring value of the transmission source at the monitoring station position by using an ITU-R P.1546 algorithm, and recording the value; thirdly, continuously calculating outwards, when the calculated value is smaller than the receiving sensitivity, combining the previous value, finding a point where the receiving sensitivity is equal to, recording the point, and finishing the calculation of the ray; fourthly, calculating points where the ray is equal to the receiving sensitivity according to all the rays, and connecting the points to form a transmitting power contour line; fifthly, connecting points of adjacent rays according to ray subdivision points calculated by the contour line to form a small polygon, and performing color filling according to a field intensity value to form a field intensity coverage map;
the field intensity coverage calculation function module is used for displaying a field intensity coverage map drawn by the model in a map by combining with a two-dimensional electronic map, and can further realize the following functions: monitoring coverage capability display of a monitoring network of a specific radio station; calculating the coverage area and the area coverage rate of a monitoring site/monitoring network; displaying or hiding the coverage map of the monitoring sites randomly, and increasing or deleting the number of the monitoring sites simultaneously; automatically comparing the coverage maps under different parameters; calling a local station database, checking the monitoring condition of a certain type of radio station, and counting the monitoring coverage rate of various stations; establishing key monitoring areas such as airports and borders by using a map; calculating the fixed point field intensity of any position on the map;
the coverage model modification module based on drive test data is optimized by using a standard macro cell model commonly used in network planning, and referring to fig. 5, a model formula is as follows: prx ═ Ptx-Ploss;
wherein:
Ploss＝
K1+K2×log(d)+K3×Hms+K4×log(Hms)+K5×log(Heff)+K6×log(Heff)×log(d)+K7×diffn+Clutter_Loss
prx: the received signal field strength (dBm) of the receiver;
ptx: effective radiated power (dBm) of the transmitter;
ploss: path loss (dB);
d: distance (km) between base station and mobile station;
hms: height (m) above the ground of the mobile station;
heff: the effective height (m) of the base station antenna;
diffn: diffraction losses calculated using the Epstein Peterson, Deygout, or butlington equivalent knife edge diffraction method;
k1& K2: intercept and slope, these factors correspond to a fixed offset and a multiplicating Factor of the logarithm of the distance between the base station and the mobile station;
k3: height factor of the mobile antenna, which parameter should typically be the effective height of the mobile station antenna;
k4: hms Multiplying Factor of Okumura-Hata;
k5: effective antenna height gain Factor, which is the effective antenna height log value's Multiply Factor;
k6: log (heff) x log (d) Factor, which is an Okumura-Hata type multiplyng Factor of log (heff) x log (d) value;
k7: diffraction coefficient, which is the Multiplying Factor value calculated for diffraction;
clutter _ Loss: a ground object loss factor;
the detailed correction calibration method is as follows: firstly, analyzing the GPS data and the monitoring station data to obtain the information of the average field strength value Prx of the received signal, the effective radiation power Ptx of the transmitter, the distance d between the transmitting system and the monitoring station, the height Hms of the transmitting system on the ground, the effective height Heff of the antenna of the monitoring station, the diffraction Loss diffn, the ground object Loss factor Clutter _ Loss and the like at each moment of the monitoring station. And then starting to carry out correction calculation, wherein a specific algorithm is as follows:
reading a group of test data: ptx, Prx, d, Hms, Heff, diffn, Clutter _ Loss;
step two, let yn be Ptx-Prx; xn ═ log (d), Hms, log (hms), log (heff)) ], log (heff)) × log (d), diffn, Clutter _ Loss;
substituting xn and yn into a formula to perform recursive calculation to obtain beta n;
reading the next group of test data, if no test data exists, turning to the fifth step, otherwise, turning to the second step;
step five, calculating the mean square error by combining actual test data by using a model formula and the coefficient of beta n (K1, K2, K3, K4, K5, K6 and K7) obtained by recursion calculation; if the mean square error is less than 8, finishing the algorithm, and replacing each K value in the model formula by beta n to obtain a formula which is the corrected model; if the mean square error is larger than or equal to 8, other actual test data are needed to be used for continuous correction, namely, the step one is carried out.
The monitoring networking area coverage analysis module can obtain the networking coverage condition of the current whole area by utilizing the coverage analysis chart of each radio monitoring station, mainly utilizes the radio monitoring station management function of a two-dimensional electronic map, and the system automatically counts and compares repeated coverage areas among stations to obtain the monitoring station corresponding to the maximum field intensity of the area and obtain the coverage depth index of the monitoring network, wherein the monitoring coverage depth is as follows: the monitoring method refers to the number of overlapping coverage layers of a monitoring network to a certain area, and if K stations in the monitoring network can cover the area, the monitoring coverage depth of the area is K. Meanwhile, the emission sources under various EIRP power levels can be subjected to full-network coverage analysis according to the requirements of the Ministry of industry and information, and a monitoring coverage index is obtained, wherein the monitoring coverage index is as follows: within a particular geospatial area, the ratio of the area of the coverage area to the total area of the area is monitored. To facilitate the actual test evaluation, referring to the concept of "communication probability" in mobile communication, it can be approximated by the following formula: the number of sampling test points/the number of total sampling test points which meet the requirements of effective field intensity measurement and effective spectrum analysis is multiplied by 100 percent;
the service management module supports the planning configuration function of typical services based on the module, such as conventional trunking private network communication service, civil aviation communication service, public network communication service, radar communication service and satellite communication service. Meanwhile, the system also has the management function of an actual station under corresponding services, mainly comprises station GPS position information, station setting basic information, station basic parameters and the like, and is convenient for key monitoring application to provide rapid monitoring coverage statistics.
The GIS geographic information management subsystem 20 mainly includes a basic operation function of a two-dimensional/three-dimensional electronic map and a map information calling function. Basic operation functions such as zoom-in, zoom-out, pan-shift, positioning, and drawing functions of a map; the map information calls functions such as acquisition of information of altitude, building height, terrain and landform and the like in a coverage area, and information of map joints and the like required by field intensity drawing. The invention adopts a three-dimensional electronic map to establish a propagation model, plan a road test route and determine a transmitting point, and the two-dimensional electronic map mainly carries out field intensity area distribution display, monitoring station management, radio station management, query function, GPS track data management and the like.
The standard transmitting subsystem 30 comprises a signal source, a power amplifier, an omnidirectional transmitting antenna and a GPS, wherein a standard transmitting source is established to cover any frequency point signal generation of 20MHz-6GHz, EIRP power of 1W/3W/5W/10W power grade is generated after amplification through a power amplifier, a standard source is provided for the drive test subsystem 40, and the GPS ensures that the signal level acquired by the monitoring system is synchronous with the transmitting source in time, so that data calculation at the same time of transmitting and receiving is facilitated. In order to effectively implement the drive test process, the standard transmitting system is integrated in the vehicle, and a power meter is used for tracing and calibrating.
The drive test subsystem 40 is completed by a standard transmitting module, a radio monitoring module, an RMTP monitoring network protocol module, a special monitoring transmission network module, a GPS positioning device and a GIS three-dimensional electronic map together, and mainly obtains the signal field intensity of the standard transmitting subsystem 30 through the monitoring function of the radio monitoring module, and finally obtains the actual coverage information of the monitoring station by combining the map information and the GPS track. The detailed process is as follows: planning the position of a maximum edge point, acquiring a drive test frequency point and the field intensity of an isoline under a power grade by using a three-dimensional electronic map and technical parameters of a radio monitoring station and adopting an ITU-R.P1546 propagation model, and planning the position area of the drive test edge; the second step is that: the method comprises the following steps of (1) planning a test route and a moving vehicle speed, wherein the route planning has automatic and manual planning functions, automatic planning automatically calculates a main road test path in a current edge position area, and the calculation method is determined by combining topographic and geomorphic information acquired by a three-dimensional electronic map; determining information such as vehicle moving speed, acquisition time, path length and the like according to the Lee's law step; the third step: the method comprises the steps of acquiring and processing drive test data, defining two data acquisition schemes of synchronous and asynchronous, referring to the attached drawings 3 and 4, acquiring data of a standard transmitting module and a monitoring system simultaneously by an RMTP monitoring network protocol module and a special monitoring transmission network module which are based on GPS time information and control a radio monitoring module to work, and automatically screening the data by software based on the GPS time synchronous information. Asynchronous acquisition does not need RMTP monitoring network protocol relative to synchronous acquisition, a standard transmitting module is responsible for receiving GPS geographic data and controlling signal source parameters, a monitoring system control center is responsible for controlling monitoring system parameters and receiving monitoring signal levels, communication exchange is carried out through a special monitoring transmission network module, namely the monitoring system control center issues a transmitting instruction, the standard transmitting module transmits corresponding information back to the control center after finishing instruction requirements, and the control center processes and connects corresponding data in series. The collected data needs to be processed as follows: removing repeated data in the drive test data; removing the drift data without latitude and longitude or the latitude and longitude; combining and connecting the drive test data and the GPS data in series to form a data format which can be displayed in a map; the fourth step: displaying the actually measured data in a two-dimensional electronic map in different colors according to the frequency points and the field intensity values; the fifth step: determining a coverage range according to the measured data, mainly comprising importing, maintaining, sorting, analyzing and combining the drive test data to finally form data convenient to be analyzed and processed, managing, analyzing and counting the drive test data, mainly counting the matching degree of a propagation model calculation result and the drive test data, and simulating the coverage range under the measured condition according to a sensitivity threshold of a monitoring system.
The station management subsystem 50: the system mainly has the functions of managing reported equipment in a radio monitoring station and a station database, managing a key transmitting source, simulating a newly-built monitoring station and the like, and comprises GPS (global positioning system) position information, station construction information and specific technical parameters of the station, wherein the key transmitting source management can be used for determining the position of a legal station when interference occurs, rapidly checking and removing the position of the legal station, and reducing the time for determining the interference position.
The data management subsystem 60: based on the functions of data storage, operation, data compression, data encryption and the like of the database, data operation under multiple dimensions of different services, different EIRP power levels and the like is considered, quick response and correctness screening of the data operation become important, and the key point is the definition of each data attribute. All data are connected in series based on the GPS time information, and the uniqueness of the data can be ensured.
It should be noted that, in the embodiment of the present application, referring to fig. 2, a method for evaluating and analyzing coverage of a radio monitoring station includes the following steps:
s1: marking sites needing simulation of building the sites on the loaded three-dimensional electronic map, wherein the sites comprise site names, site heights, site position information, antenna types, antenna gains and nominal monitoring sensitivity;
s2: setting site service parameters to be simulated, wherein the site service parameters comprise service frequency bands/points, transmitting end height information, geographic information, EIRP power levels and transmitting bandwidths;
s3: drawing rays in a three-dimensional map by taking a monitoring station as a center according to a preset angle, and virtually establishing a plurality of emission sources on each ray;
s4: calculating the monitoring value of the emission source at the monitoring station position by using an ITU-R.P1546 algorithm, and continuously calculating the outside on the ray;
s5: when the calculated value is smaller than the monitoring receiving sensitivity, finding the end point of the sensitivity contour line by combining the last value, namely the ray calculation is completed;
s6: repeating S3 to S6 to complete contour points of all rays and connecting all the points to form a power contour curve;
s7: adjusting other service parameters, wherein the other service parameters comprise frequency and EIRP power, and drawing other power contour curves under different service conditions;
s8: other ray subdivision points in each power contour line are connected to form a plurality of boundary areas, and the field intensity is distinguished through color gradient to form a field intensity coverage map;
s9: and displaying the simulated field intensity coverage range of the monitoring station based on the two-dimensional electronic map.
Further, in the most preferred embodiment of the present invention, when the field intensity coverage evaluation is performed by using actual drive test data, the function of the drive test system is used, which specifically includes the following steps:
s101: preparation before drive test: the method comprises the steps of planning and setting parameters of a transmitting subsystem, a GPS antenna, a network transmission tool, an RMTP monitoring network protocol and monitoring site information;
s102: planning a driving route: selecting at least one main road according to a three-dimensional electronic map, formulating a driving speed and a driving distance according to the Lee's law, carrying out reverse thrust according to the EIRP power of a transmitting system and the sensitivity of a monitoring system, and determining a sensitivity power contour line area;
s103: actual measurement data acquisition: starting a transmitting system on the road test vehicle, starting a GPS positioning device, traveling on a planned road according to a preset vehicle speed, and simultaneously controlling the transmitting system and the monitoring system to record monitoring data, transmitting time and acquiring time according to a certain time interval until all roads are acquired;
s104: and synchronously finishing data acquisition: selecting and controlling a transmitting system to transmit signal requirements of different services, completing data required by a plurality of services by one drive test, and independently repeating the drive test of other services;
s105: and (3) drive test data processing: removing repeated data in the drive test data, removing data without latitude and longitude or drift of latitude and longitude, combining the drive test data and GPS data and importing the data into a system to form a data format which can be displayed on a map;
s106: and (3) actual measurement data display: the processed data can be displayed on a map according to a preset rule, and the data of each road measuring point is drawn on the map;
s107: the method comprises the steps of utilizing measured data to draw a coverage range of a monitoring station, setting a threshold, deducing a plurality of points which are close to or equal to receiving sensitivity according to the drive test data, connecting the points, drawing a measured coverage range, and comparing the measured coverage range with a calculated coverage range.
Further, in the most preferred embodiment of the present invention, if the prediction calculation model of ITU-r.p1546 needs to be corrected based on the drive test data, the method for correcting the prediction model is used, which specifically includes the following steps:
s201: analyzing the GPS data and the monitoring station data to obtain a received signal average field strength value Prx of the monitoring station at each moment, effective radiation power Ptx of a transmitter, a distance d between a transmitting system and the monitoring station, a height Hms on the ground where the transmitting system is located, an effective height Heff of an antenna of the monitoring station, diffraction Loss diffn and a ground object Loss factor Clutter _ Loss;
s202, reading the test data in the step S201: the average field strength value Prx of the received signal, the effective radiation power Ptx of the transmitter, the distance d between the transmitting system and the monitoring station, the height Hms of the transmitting system on the ground, the effective height Heff of the antenna of the monitoring station, the diffraction Loss diffn and the ground feature Loss factor Clutter _ Loss, and the correction calculation is carried out.
Further, in a most preferred embodiment of the present invention, the method for predicting ITU-r.p1546-5 field strength comprises:
s301: determining the first propagation type, the second propagation type and the third propagation type as the types of terrestrial, cold-sea or warm-sea propagation paths, respectively, and if the propagation types are mixed paths, determining a combination of two path types that can be considered as the first propagation type and the second propagation type, and if the paths can be represented by a single type, considering it as the first propagation type;
s302: determining two nominal time percentages;
s303: determining two nominal frequencies;
s304: determining nominal distances of a low end and a high end;
s305: for the first transmission type, proceed as S306 to S311;
s306: for the low end nominal percentage of time, proceed as S307 to S310;
s307: for the low end nominal frequency, proceed as S308 and S309;
s308: for the required distance d and the transmitting/base station antenna height h1, a field strength exceeding 50% of the location point is obtained for the receiving/mobile station antenna at the representative scattering height R above the earth:
s309: if the desired frequency does not match the low-end nominal frequency, repeating S308 for the high-frequency nominal frequency and interpolating or extrapolating the two field strengths;
s310: if the percentage of time required does not match the lower nominal percentage of time, repeating S307 to S309 for the higher nominal percentage of time and interpolating the two field strengths;
s311: if the hybrid path is predicted, the path requirements for each propagation type are from S306 to S310;
s312: correcting the field strength of the terrain clearance angle at the receiver/mobile unit if information in terms of the terrain clearance angle at the receiver/mobile station antenna adjacent the land is available;
s313: calculating an estimated value of field intensity caused by scattering, and taking the maximum value of E and Ets;
s314: correcting the field strength of the antenna height h2 of the receiving/moving station;
s315: if clutter exists around the transmitting/base station terminal, correcting the effect;
s316: correcting the inclined diameter;
s317: if a path less than 1km is given, S301 to S316 for d ═ 1km need to be followed;
s318: if the field strength of the position percentage of the receiving/mobile station antenna adjacent to the land is required to be not 50 percent, correcting the field strength of the required position percentage;
s319: limiting the obtained field strength to a preset maximum value, if the field strength calculation of the hybrid path has been made for less than 50% of the time percentage, the maximum field strength value needs to be calculated by linear interpolation between all terrestrial field strength values and all sea field strength values, which is given by the following formula:
Emax＝Efs+dsEse/dtotaldB(V/m)
s320: the field strength is converted into an equivalent fundamental transmission loss of the path.
Further, in the most preferred embodiment of the present invention, in S302, the method comprises:
the percentage of time required is >1 and <10, where the nominal percentage of time for the low and high ends is 1% and 10%, respectively;
the percentage of time required is >10 and <50, where the nominal percentage of time for the low and high ends is 10% and 50%, respectively;
if the percentage of time required is equal to 1%, 10%, or 50%, then the value may be considered the low end nominal percentage of time without requiring the interpolation process of S310.
Further, in the most preferred embodiment of the present invention, in S303, the method further includes:
when the required frequency is less than 600MHz, the low-end frequency and the high-end frequency are respectively 100MHz and 600 MHz;
when the required frequency is less than 600MHz, the low-end frequency and the high-end frequency are 600MHz and 2000MHz respectively;
if the desired frequency is equal to 100 or 600 or 2000MHz, this value may be considered the low end nominal frequency without the need for interpolation/extrapolation processing in S309.
Further, in the most preferred embodiment of the present invention, in S308, the method comprises:
for a transmit/base station antenna height h1 of less than 10m, the field strength at the required height and distance is determined.
Further, in the most preferred embodiment of the present invention, it specifically further comprises:
s381: determining the values of the low and high nominal heights h1, which should be considered as the low nominal height value h1 if the value of h1 matches any of the nominal values 10, 20, 37.5, 75, 150, 300, 600 or 1200m, without interpolation;
s382: for the low-end nominal value h1, proceeding as S383 to S385;
s383: for the low end nominal distance value, proceed as per S384;
s384: for the required distance d and the transmitting/base station antenna height h1, a field strength exceeding 50% of the location point is obtained for the receiving/mobile station antenna at the representative spread height R;
s385: if the desired distance does not match the lower nominal distance, repeating S384 for the upper nominal distance and interpolating the two field strengths over that distance;
s386: if the desired transmit/base station antenna height h1 does not meet any of the nominal values, S383 through S385 are repeated and h1 is interpolated/extrapolated.
The present embodiment is directed to radio wave propagation model selection: for example, during the evaluation of broadcast television coverage network planning and base station coverage network planning, the models are established by adopting the common Okumura-Hata Onchun model and COST-231Hata model, but the models are relatively bound in the frequency of 100MHz-2GHz and the application environment of urban, suburban or open areas, and cannot adapt to all the geographic environments and the frequency requirements of monitoring systems. In the embodiment, an ITU-R.P1546 model is adopted for modeling, and the ITU-R P.1546 recommendation provides a method for predicting the electric wave propagation field intensity from point to surface in a frequency band of 30MHz to 3000MHz, wherein the method is characterized in that interpolation and extrapolation are carried out on the basis of actual measurement data of 100MHz, 600MHz and 2000MHz to obtain a field intensity curve in the electric wave propagation process, and the curve is a function of distance, antenna height, frequency and time percentage. The method also provides typical barrier heights of different terrains, and provides a more accurate correction method for the height of the receiving antenna. The proposal also gives the relationship of field intensity as a function curve of specific frequency, specific distance and specific receiving and transmitting antenna height, and obtains the radio wave propagation loss under any path by using interpolation and extrapolation methods. The model can obtain a completely consistent prediction result under the condition required by the Okumura-Hata Oncun model, so that the model is better in terms of geographic factors, frequency ranges and the like compared with the traditional propagation model.
Further, in the present embodiment, in terms of a standard transmitting module used for the drive test: the traditional standard transmitting module basically ignores the function of self-tracing calibration, the output EIRP power is fixed, and the embodiment has the function of self-calibration of a transmitting system and can correct the connecting end of a transmitting antenna. Meanwhile, the 50W broadband power amplifier configured in the embodiment ensures that EIRP power output of a plurality of power levels is met in the frequency range of 6 GHz. In addition, the standard transmitting module is additionally provided with a GPS function, so that any output information and GPS time can be ensured to be marked one by one.
Compared with base station coverage evaluation, the method and the device can only carry out drive test verification after model evaluation, adjust the height or the downward inclination angle of the base station antenna according to the drive test result, and improve the coverage requirement. The embodiment includes three coverage evaluation functions, which are coverage evaluation based on the ITU-r.p1546 propagation model, coverage evaluation based on the drive test, and a correction propagation model based on the drive test data, and the model coverage evaluation can be performed again based on the model data after the model correction is completed, thereby reducing the complex workload of the repetitive drive test.
The embodiment applies the three-dimensional electronic map to the field of radio monitoring for the first time, and electromagnetic spectrum situation display is reflected more truly.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (1)
1. A radio monitoring site coverage area evaluation and analysis method applied to a radio monitoring site coverage area evaluation and analysis system is characterized by comprising the following steps:
the radio monitoring site coverage evaluation and analysis system comprises:
the coverage analysis subsystem is used for acquiring monitoring system parameters, radio wave propagation loss and radio emission source parameters, and modeling and analyzing according to the acquired parameters;
the GIS geographic information management subsystem comprises a two-dimensional/three-dimensional electronic map module and a map information calling module, and is used for establishing a propagation module, planning a road test route, determining a transmitting point, displaying field intensity area distribution, managing a monitoring station, managing a radio station, inquiring a function and managing GPS track data;
the standard transmitting subsystem comprises a signal source, a power amplifier, an omnidirectional transmitting antenna and a GPS positioning device, and the standard transmitting subsystem generates EIRP power of 1W/3W/5W/10W power grade by covering any frequency point signal of 20MHz-6GHz and amplifying the signal by a power amplifier so as to provide a standard transmitting source for the road test subsystem;
the system comprises a road test subsystem, a monitoring subsystem and a monitoring subsystem, wherein the road test subsystem comprises a standard transmitting module, a radio detection module, an RMTP monitoring network protocol module, a special monitoring transmission network and a GIS three-dimensional electronic map module, and acquires the signal field intensity of the standard transmitting subsystem through the monitoring function of a radio monitoring system, and finally obtains the actual coverage information of a monitoring station by combining the map information and a GPS track;
the station management subsystem comprises a radio monitoring station management module, a radio station management module, a station database calling module and a virtual station management module, and is used for rapidly checking and removing the position of a legal station when interference occurs and reducing the time for determining the interference position;
the data management subsystem is used for data storage, data calling, data compression, data encryption and data attribute definition;
the coverage analysis subsystem also comprises a propagation model calculation module, a virtual station building function module, an isoline field intensity prediction calculation module, a field intensity coverage calculation function module, a coverage model correction module based on drive test data, a monitoring networking area coverage analysis module and a typical service management module;
the propagation model calculation module is used for providing a field intensity prediction proposal of 30MHz-3000MHz frequency band ground service, the proposal comprises the correction of interpolation/extrapolation to the field intensity under different transmitting power, frequency, antenna height and transmission distance, and any field intensity value of which the transmitting frequency is 30MHz-3000MHz and the propagation distance is 1 km-1000 km, which cannot be directly read from the provided curve, is obtained by using a correction formula;
the virtual station building functional module is used for realizing automatic calculation of any frequency based on various radio wave propagation models and three-dimensional electronic maps and drawing field intensity coverage situation maps at different power levels;
the isoline field intensity prediction calculation module is used for automatically calculating and drawing isolines of different service frequencies under different transmitting power levels to form a boundary area;
the field intensity coverage calculation function module is used for displaying a field intensity coverage map drawn by the model in a map by combining a two-dimensional electronic map, and realizing the functions of monitoring coverage capability display of a monitoring network of a specific radio station, calculating the coverage area and the area coverage rate of a monitoring station/monitoring network, displaying or hiding the coverage map of the monitoring station at will, increasing or deleting the number of the monitoring stations at the same time, automatically comparing the coverage maps under different parameters, calling a local station database, checking the monitoring condition of a certain type of radio station, counting the monitoring coverage rate of various stations, establishing a key monitoring area by using the map and calculating the fixed point field intensity of any position on the map;
the coverage model modification module based on the drive test data is used for optimizing by adopting a standard macro cell model commonly applied in network planning;
the monitoring networking area coverage analysis module obtains the networking coverage condition of the current whole area by utilizing the coverage analysis map of each radio monitoring site, and obtains the monitoring site corresponding to the maximum field intensity of the area by utilizing the radio monitoring site management function of the two-dimensional electronic map to obtain the coverage depth index of the monitoring network;
the typical service management module is used for a key monitoring application to provide rapid monitoring coverage statistics;
the propagation model calculation module adopts ITU-R.P1546 standard;
the method for evaluating and analyzing the coverage area of the radio monitoring station comprises the following steps:
s1: marking sites needing simulation of building the sites on the loaded three-dimensional electronic map, wherein the sites comprise site names, site heights, site position information, antenna types, antenna gains and nominal monitoring sensitivity;
s2: setting site service parameters to be simulated, wherein the site service parameters comprise service frequency bands/points, transmitting end height information, geographic information, EIRP power levels and transmitting bandwidths;
s3: drawing rays in a three-dimensional map by taking a monitoring station as a center according to a preset angle, and virtually establishing a plurality of emission sources on each ray;
s4: calculating the monitoring value of the emission source at the monitoring station position by using an ITU-R.P1546 algorithm, and continuously calculating the outside on the ray;
s5: when the calculated value is smaller than the monitoring receiving sensitivity, finding the end point of the sensitivity contour line by combining the last value, namely the ray calculation is completed;
s6: repeating S3 to S6 to complete contour points of all rays and connecting all the points to form a power contour curve;
s7: adjusting other service parameters, wherein the other service parameters comprise frequency and EIRP power, and drawing other power contour curves under different service conditions;
s8: other ray subdivision points in each power contour line are connected to form a plurality of boundary areas, and the field intensity is distinguished through color gradient to form a field intensity coverage map;
s9: displaying the coverage range of the field intensity of the simulated monitoring station based on a two-dimensional electronic map;
when field intensity coverage evaluation is carried out on actual drive test data, the drive test system function is used, and the method specifically comprises the following steps:
s101: preparation before drive test: the method comprises the steps of planning and setting parameters of a transmitting subsystem, a GPS antenna, a network transmission tool, an RMTP monitoring network protocol and monitoring site information;
s102: planning a driving route: selecting at least one main road according to a three-dimensional electronic map, formulating a driving speed and a driving distance according to the Lee's law, carrying out reverse thrust according to the EIRP power of a transmitting system and the sensitivity of a monitoring system, and determining a sensitivity power contour line area;
s103: actual measurement data acquisition: starting a transmitting system on the road test vehicle, starting a GPS positioning device, traveling on a planned road according to a preset vehicle speed, and simultaneously controlling the transmitting system and the monitoring system to record monitoring data, transmitting time and acquiring time according to a certain time interval until all roads are acquired;
s104: and synchronously finishing data acquisition: selecting and controlling a transmitting system to transmit signal requirements of different services, completing data required by a plurality of services by one-time drive test, or independently repeating drive test on other services;
s105: and (3) drive test data processing: removing repeated data in the drive test data, removing data without latitude and longitude or drift of latitude and longitude, combining the drive test data and GPS data and importing the data into a system to form a data format which can be displayed on a map;
s106: and (3) actual measurement data display: the processed data can be displayed on a map according to a preset rule, and the data of each road measuring point is drawn on the map;
s107: drawing a monitoring station coverage range by utilizing the measured data, setting a threshold, deducing a plurality of points close to or equal to receiving sensitivity according to the drive test data, connecting the points, drawing a measured coverage range, and comparing the measured coverage range with the calculated coverage range;
if the prediction calculation model of ITU-R.P1546 needs to be corrected based on the drive test data, the method for correcting the prediction model is used, and the method specifically comprises the following steps:
s201: analyzing the GPS data and the monitoring station data to obtain a received signal average field strength value Prx of the monitoring station at each moment, effective radiation power Ptx of a transmitter, a distance d between a transmitting system and the monitoring station, a height Hms on the ground where the transmitting system is located, an effective height Heff of an antenna of the monitoring station, diffraction Loss diffn and a ground object Loss factor Clutter _ Loss;
s202, reading the test data in the step S201: the average field strength value Prx of the received signal, the effective radiation power Ptx of the transmitter, the distance d between the transmitting system and the monitoring station, the height Hms of the transmitting system on the ground, the effective height Heff of the antenna of the monitoring station, the diffraction Loss diffn and the ground feature Loss factor Clutter _ Loss, and the correction calculation is carried out.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010646709.4A CN111818555B (en) | 2017-06-30 | 2017-06-30 | Radio monitoring station coverage area evaluation and analysis method based on virtual station building |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010646709.4A CN111818555B (en) | 2017-06-30 | 2017-06-30 | Radio monitoring station coverage area evaluation and analysis method based on virtual station building |
| CN201710519077.3A CN107333284B (en) | 2017-06-30 | 2017-06-30 | A system and method for evaluating and analyzing coverage of radio monitoring site |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710519077.3A Division CN107333284B (en) | 2017-06-30 | 2017-06-30 | A system and method for evaluating and analyzing coverage of radio monitoring site |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111818555A CN111818555A (en) | 2020-10-23 |
| CN111818555B true CN111818555B (en) | 2022-03-22 |
Family
ID=60198348
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010646709.4A Active CN111818555B (en) | 2017-06-30 | 2017-06-30 | Radio monitoring station coverage area evaluation and analysis method based on virtual station building |
| CN201710519077.3A Active CN107333284B (en) | 2017-06-30 | 2017-06-30 | A system and method for evaluating and analyzing coverage of radio monitoring site |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710519077.3A Active CN107333284B (en) | 2017-06-30 | 2017-06-30 | A system and method for evaluating and analyzing coverage of radio monitoring site |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN111818555B (en) |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107846689B (en) * | 2017-11-16 | 2021-09-14 | 上海微波技术研究所（中国电子科技集团公司第五十研究所） | Method for realizing wireless communication network frequency planning based on GIS system |
| CN109495201B (en) * | 2018-12-04 | 2020-05-19 | 苏州恩巨网络有限公司 | Station identification method, device and equipment |
| CN109862572B (en) * | 2019-02-25 | 2022-06-28 | 全球能源互联网研究院有限公司 | Deployment method, evaluation method and electronic device of intelligent interference analysis system |
| CN112020091B (en) * | 2019-05-31 | 2023-04-07 | 中国移动通信集团浙江有限公司 | Weak coverage processing method and device based on digitized geographic information |
| DE102019210659A1 (en) * | 2019-07-18 | 2021-01-21 | Robert Bosch Gmbh | Method for generating a three-dimensional model of the environment using GNSS measurements |
| CN110944338B (en) * | 2019-11-29 | 2022-06-24 | 中建三局智能技术有限公司 | Design method for digital cluster type wireless talkback signal full coverage |
| CN111491321A (en) * | 2020-04-22 | 2020-08-04 | 上海创远仪器技术股份有限公司 | Method for realizing test processing aiming at coverage area of radio monitoring station and test system thereof |
| CN113839729B (en) * | 2020-06-08 | 2024-05-03 | 海鹰航空通用装备有限责任公司 | Frequency modulation broadcast signal coverage prediction evaluation device and method |
| CN113325241A (en) * | 2021-06-03 | 2021-08-31 | 南京嘉安网络技术有限公司 | Method for measuring and processing radio signal field intensity based on position approaching direction |
| CN113917391A (en) * | 2021-09-14 | 2022-01-11 | 成都德辰博睿科技有限公司 | A mobile station positioning subsystem |
| CN114520978B (en) * | 2021-12-18 | 2022-11-18 | 哈尔滨翰才科技有限公司 | Method and system for automatically arranging base stations in network planning simulation |
| CN114423015A (en) * | 2021-12-28 | 2022-04-29 | 上海市无线电监测站 | Vehicle-mounted system for rapid test and comparison of in-use station and implementation method |
| CN114679743B (en) * | 2022-05-30 | 2022-08-16 | 广州天越电子科技有限公司 | Coverage effect evaluation system supporting multiple signal frequency bands |
| CN115102645A (en) * | 2022-06-01 | 2022-09-23 | 苏州恩巨网络有限公司 | Method, apparatus, apparatus and medium for predicting signal strength of radio wave propagation |
| CN115767560A (en) * | 2022-09-23 | 2023-03-07 | 代书胜 | Data monitoring system based on block chain |
| CN115616616B (en) * | 2022-12-19 | 2023-03-14 | 涟漪位置(广州)科技有限公司 | Method, device, equipment and medium for determining reliability of satellite positioning system |
| CN116248207A (en) * | 2023-02-13 | 2023-06-09 | 北京吉亚斯科技有限公司 | Gias Microwave Interference Prediction Model |
| CN116680348A (en) * | 2023-02-13 | 2023-09-01 | 北京吉亚斯科技有限公司 | Electromagnetic Simulation Diffraction Model |
| WO2024176256A1 (en) * | 2023-02-24 | 2024-08-29 | Jio Platforms Limited | System and method for optimum site selection and deployment |
| CN116452246B (en) * | 2023-06-19 | 2023-09-19 | 航天宏图信息技术股份有限公司 | Electromagnetic equipment performance and address selection evaluation method and device and electronic equipment |
| CN116707689B (en) * | 2023-07-05 | 2023-11-21 | 哈尔滨广播器材有限责任公司 | Method and system for distributing monitoring equipment of medium wave broadcast transmitter |
| CN117499906B (en) * | 2023-12-29 | 2024-04-16 | 四川安信科创科技有限公司 | An emergency tactical wireless communication system for high mountain valleys and its design method |
| CN117906481A (en) * | 2024-01-22 | 2024-04-19 | 深圳市华怡丰科技有限公司 | Photoelectric sensor signal linearization processing method and related equipment |
| CN118425881B (en) * | 2024-07-05 | 2024-10-22 | 云南电力试验研究院(集团)有限公司 | Method for calculating electromagnetic wave propagation path and selecting passive positioning measuring station |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101349718A (en) * | 2007-07-20 | 2009-01-21 | 深圳市家国天下科技有限公司 | Method, system and apparatus for generating electromagnetic field three-dimensional model |
| CN102413479A (en) * | 2010-09-26 | 2012-04-11 | 北京迪特卡得通信设备有限公司 | Dynamic evaluation method for cell coverage in communication network |
| CN102625326A (en) * | 2012-03-22 | 2012-08-01 | 浙江大学 | Method and device for multicell joint optimization under coverage of cellular mobile communication network |
| WO2016146788A1 (en) * | 2015-03-17 | 2016-09-22 | Cyprus University Of Technology | System and method for predicting solar power generation |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE60232505D1 (en) * | 2001-04-03 | 2009-07-16 | At & T Mobility Ii Llc | METHOD AND DEVICE FOR ESTIMATING THE MOBILE STATION POSITION |
| CN101363892A (en) * | 2008-09-10 | 2009-02-11 | 俞惟铨 | Radio monitoring and electronic interference integration system |
| CN101835192B (en) * | 2010-04-16 | 2013-08-07 | 华为技术有限公司 | Method for measuring network signal coverage, network controller and communication system |
| US9088552B2 (en) * | 2011-11-30 | 2015-07-21 | Motorola Solutions, Inc. | Method and apparatus for key distribution using near-field communication |
| CN105722135A (en) * | 2016-03-08 | 2016-06-29 | 福建泰克通信有限公司 | Border base station road test system |
| CN106131867A (en) * | 2016-08-31 | 2016-11-16 | 成都九华圆通科技发展有限公司 | A kind of radio cloud monitoring system |
-
2017
- 2017-06-30 CN CN202010646709.4A patent/CN111818555B/en active Active
- 2017-06-30 CN CN201710519077.3A patent/CN107333284B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101349718A (en) * | 2007-07-20 | 2009-01-21 | 深圳市家国天下科技有限公司 | Method, system and apparatus for generating electromagnetic field three-dimensional model |
| CN102413479A (en) * | 2010-09-26 | 2012-04-11 | 北京迪特卡得通信设备有限公司 | Dynamic evaluation method for cell coverage in communication network |
| CN102625326A (en) * | 2012-03-22 | 2012-08-01 | 浙江大学 | Method and device for multicell joint optimization under coverage of cellular mobile communication network |
| WO2016146788A1 (en) * | 2015-03-17 | 2016-09-22 | Cyprus University Of Technology | System and method for predicting solar power generation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111818555A (en) | 2020-10-23 |
| CN107333284A (en) | 2017-11-07 |
| CN107333284B (en) | 2020-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111818555B (en) | Radio monitoring station coverage area evaluation and analysis method based on virtual station building | |
| CN101060689B (en) | A method and equipment for planning the communication system network | |
| EP1191804B1 (en) | Path loss data normalization for growth management of a cellular system | |
| US9473950B2 (en) | Method for the prediction of coverage areas of a cellular network | |
| CA2394508C (en) | Method and apparatus for network planning | |
| CN102883338A (en) | Correction method for propagation model in TD-LTE system | |
| WO2023214176A1 (en) | A method of fast path loss calculation considering environmental factors | |
| CN108259097B (en) | Method and device for correcting wireless propagation model based on MR data | |
| Hanci et al. | Mobile radio propagation measurements and tuning the path loss model in urban areas at GSM-900 band in Istanbul-Turkey | |
| Zhou et al. | AMPLE: An adaptive multiple path loss exponent radio propagation model considering environmental factors | |
| US11533633B2 (en) | Frequency planning for a communication system | |
| CN211453787U (en) | Simulation field intensity test system | |
| CN107222874A (en) | A kind of LTE propagation model revision methods based on variable critical distance | |
| CN101795163A (en) | System for measuring, calculating and evaluating signal covering quality of digital single frequency network | |
| Anamonye et al. | Evaluation and analysis of gsm signals in warri | |
| KR101355363B1 (en) | High-Speed Radio Propagation Analysis Apparatus and Method | |
| KR102119121B1 (en) | Mrthod for calibrating loss of radio-propagation using clutter-morphology calibrating factor based on measurement data | |
| CN106714111A (en) | Positioning method and device for measuring report | |
| Qing | GIS aided radio wave propagation modeling and analysis | |
| Storsaeter et al. | A GIS-based tool for optimizing C-ITS communication infrastructure | |
| Wang et al. | Spectrum environment map dataset based on ray-tracer under 3D dynamic built-up scenarios | |
| Ai et al. | Design of power wireless private network coverage prediction system | |
| JP7804236B2 (en) | Received power estimation device, received power estimation method, and program | |
| Chanaga et al. | Radio Propagation Model Adjusted to Mountainous Topography for TVWS Coverage Study | |
| Rojas | A geographic information system applied to coverage maps of 3g cellular communications networks |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |