JPWO2013118681A1 - Main power area estimation device and main power area estimation method - Google Patents

Abstract

The GPS position information acquisition unit 601 and the delay position information acquisition unit 602 acquire delay information, PRACH-PD position information, and GPS position information. The main power area estimation unit 613 calculates the distribution of the mobile devices 100 based on the delay information, the PRACH-PD position information, or the GPS position information. Then, the main power area estimation unit 613 estimates the radius of the main power area where the estimation target sector is the main power based on the calculated distribution of the mobile devices 100 and the position of the antenna 201.

Description

  The present invention relates to a main power area estimation device and a main power area estimation method for estimating a main power area of a sector.

  Conventionally, it is desired to accurately obtain a main power area in which a sector formed by a base station is a main power. This sector is determined by the radio wave radiated from the antenna, and it is difficult to accurately distinguish the main power area. Therefore, for example, it is possible to actually perform radio wave measurement at each point to obtain the main power area of the sector. As another method, for example, it is conceivable to grasp the main power area by Voronoi division. For example, Patent Document 1 (Japanese Patent Laid-Open No. 2005-252548) describes that an area is divided by performing Voronoi division using the position of a base station.

JP 2005-252548 A

  However, the method of actually performing radio wave measurement at each point requires a very large number of man-hours for measurement itself, and it is difficult to measure the main power area for each sector in a wide range such as the whole of Japan. Further, in the method using Voronoi division, it is obtained based on the position of the base station, and information on the actual main power area is not taken into consideration, so it is difficult to accurately determine the main power area of the sector. It is.

  Therefore, an object of the present invention is to provide a main power area estimation apparatus and a main power area estimation method that can easily and accurately estimate a main power area of a sector.

  The main power area estimation apparatus according to the present invention obtains correspondence information that associates positioning information of a mobile device obtained by predetermined positioning with a sector identifier of a sector where the mobile device is located. And the base station forming the sector corresponding to the predetermined sector identifier, and the distribution of mobile stations obtained from the positioning position information associated with the predetermined sector identifier among the correspondence information acquired by the correspondence information acquisition unit Main power area estimation means for estimating the radius of the main power area where the sector corresponding to the predetermined sector identifier is the main power based on the position of the antenna of the station.

  Further, the main power area estimation method according to the present invention is a method for acquiring correspondence information in which positioning information of a mobile device obtained by predetermined positioning and a sector identifier of a sector where the mobile device is located are associated with each other. Out of the correspondence information acquired in the information acquisition step and the correspondence information acquisition step, the mobile station distribution obtained from the positioning position information associated with the predetermined sector identifier and the sector corresponding to the predetermined sector identifier are formed. And a main power area estimation step of estimating a radius of a main power area in which a sector corresponding to a predetermined sector identifier is a main power based on the position of the antenna of the base station to be operated.

  In these inventions, the sector corresponding to the predetermined sector identifier becomes the main force based on the distribution of mobile devices obtained from the positioning position information associated with the predetermined sector identifier and the position of the antenna. The radius of the main power area is estimated. In other words, since the radius of the main power area of the sector can be obtained based on the actual distribution of mobile devices, the main power area of the sector can be obtained accurately. In addition, it is not necessary to perform radio wave measurement at each point, and it is not necessary to incorporate a special function for determining the sector's main power area in the mobile device itself. Can be sought.

  The correspondence information acquisition means is a first positioning based on a propagation time of radio waves transmitted and received between the mobile device and the base station, and the positioning position information is obtained by the first positioning. First correspondence information that is first positioning position information and that obtains the first correspondence information when the correspondence information is first correspondence information in which the first positioning position information and the sector identifier are associated with each other. The information acquisition means, and the predetermined positioning is a second positioning different from the first positioning, the positioning position information is the second positioning position information obtained by the second positioning, and the correspondence information is Including at least one of second correspondence information acquisition means for acquiring the second correspondence information when the second positioning information and the sector identifier are associated with each other. Is preferred. In this case, the first correspondence information or the second correspondence information can be obtained by the first correspondence information acquisition means or the second correspondence information acquisition means.

  Further, the first positioning position information includes position information calculated based on propagation time information related to the propagation time, and the main power area estimation means obtains a distance between the antenna and the mobile device, and determines the obtained antenna and mobile device. It is preferable to obtain the distribution of the mobile devices based on the distance. In this case, the distribution of mobile devices is obtained using the propagation time information, and the main power area of the sector can be obtained based on this distribution.

  Further, the first positioning position information or the second positioning position information preferably includes coordinate information indicating the position of the mobile device, and the main power area estimation means preferably obtains the distribution of the mobile device based on the coordinate information. In this case, the distribution of the mobile device is obtained using the coordinate information, and the main power area of the sector can be obtained based on this distribution.

  Moreover, it is preferable to further include an antenna position calculation unit that calculates the position of the antenna based on the distribution of mobile devices obtained based on the first correspondence information. In this case, the position of the antenna can be calculated based on the distribution of mobile devices. Thereby, in addition to the main power area of the sector, the position of the antenna can be grasped. Therefore, it is possible to grasp the position of the antenna based on the distribution of the mobile devices without using the information regarding the position of the antenna included in the equipment information of the base station input by the worker or the like.

  Further, among the first correspondence information, a sector that determines the type of the sector corresponding to the predetermined sector identifier based on the distribution of mobile devices obtained from the first positioning position information associated with the predetermined sector identifier It is preferable to further include a type determination unit. In this case, the sector type can be determined based on the distribution of mobile devices. Thereby, in addition to the main power area of a sector, the type of sector can be grasped. Therefore, it is possible to grasp the sector type based on the distribution of the mobile devices without using information on the sector type included in the base station facility information input by the worker or the like.

  In addition, the radiation direction for calculating the radiation direction of the radio wave radiated from the antenna based on the distribution of the mobile devices obtained from the first positioning position information associated with the predetermined sector identifier among the first correspondence information It is preferable to further include a calculation means. In this case, the radio wave radiation direction can be calculated based on the distribution of mobile devices. Thereby, in addition to the main power area of the sector, the radiation direction of the radio wave can be grasped. Therefore, it is possible to grasp the radiation direction of radio waves based on the distribution of mobile devices without using information on the radiation direction of radio waves included in the equipment information of the base station input by an operator or the like.

  Further, based on the distribution of mobile devices obtained from the second positioning position information associated with the predetermined sector identifier among the second correspondence information, the base station forming the sector corresponding to the predetermined sector identifier It is preferable to further include a radiation width calculating means for calculating the radiation width of the radio wave radiated from the antenna. In this case, the radiation width of the radio wave radiated from the antenna can be accurately calculated based on the distribution of the mobile devices obtained from the second positioning position information associated with the predetermined sector identifier.

  According to the present invention, the main power area of a sector can be estimated easily and accurately.

It is a block diagram which shows the function structure of the mobile communication system to which the main influence area estimation apparatus is applied. It is a figure which shows the relationship between BTS and a sector. It is a figure which shows the principle of PRACH-PD positioning. It is a figure which shows delay information. It is a figure which shows the positioning point of the mobile device in a directional sector. It is a figure which shows PRACH-PD position information. It is a figure which shows PRACH-PD position information. It is a figure which shows GPS position information. It is a block diagram which shows the detail of the main power area map preparation apparatus. It is a figure which shows the positioning point of the mobile device in a directional sector. It is a figure which shows the positioning point of the mobile device in an omnidirectional sector. It is a figure which shows the result of having totaled PRACH-PD position information in a directional sector for every positioning point. It is a figure which shows the result of having totaled the PRACH-PD position information in an omnidirectional sector for every positioning point. It is a flowchart which shows the flow of the process regarding the 3rd method in a sector classification judgment part. It is a figure explaining the 1st method of calculating the position of an antenna. It is a figure explaining the 2nd method of calculating the position of an antenna. It is a figure explaining the 3rd method of calculating the position of an antenna. It is a figure explaining the 3rd method of calculating the position of an antenna. It is a figure explaining the 4th method of calculating the position of an antenna. It is the figure explaining the 5th and 6th method of calculating the position of an antenna. It is a figure explaining the 7th method of calculating the position of an antenna. It is a figure which shows the number of signals for every latitude and longitude. It is a figure which shows an intermediate | middle table. It is a figure which shows signal density distribution according to distance. It is a figure which shows the flow of a process in the main influence area estimation part. It is a figure which shows the number of signals for every propagation time in delay information. It is a figure which shows an intermediate | middle table. It is a figure which shows the flow of a process in the main influence area estimation part. It is the figure explaining the method of estimating the radius of the main influence area in the main influence area estimation part. It is a figure explaining the method of calculating a radiation direction. It is a figure which shows the number of signals for every positioning point. It is a figure explaining the method of calculating a radiation width. It is a figure which shows the graph which simplified density distribution. It is a figure which shows the main power area map. It is a figure which shows the modification of the method of producing a main influence area map.

  Next, a preferred embodiment of a mobile communication system to which a main power area estimation device and a main power area estimation method according to the present invention are applied will be described. First, the functional block configuration of the mobile communication system 10 will be described. FIG. 1 is a block diagram showing a functional configuration of a mobile communication system. As shown in FIG. 1, a mobile communication system 10 includes a mobile device 100, a BTS (base station) 200, an RNC (radio network control device) 300, an exchange 400, a GPS location management unit 501, and a location collection unit. 502 and the main power area map creation apparatus 600 are comprised.

  The mobile device 100 includes a GPS function unit 101. The GPS function unit 101 uses GPS (second positioning) to acquire detailed position information (second positioning position information) (hereinafter referred to as “GPS coordinate information”) indicating the location of the mobile device 100. Is. The GPS coordinate information is acquired when a predetermined application (an application using GPS coordinate information) installed in the mobile device 100 is executed. When the GPS function unit 101 acquires the GPS coordinate information, the GPS function unit 101 outputs the acquired GPS coordinate information to the GPS position management unit 501.

  The BTS 200 radiates radio waves for communication with the mobile device 100 from the antenna 201, and forms a sector that is a communication area. There are two types of sectors: a directional sector and an omnidirectional sector. A directional sector is a sector formed using a directional antenna. For example, as shown in FIG. 2A, when the antenna 201 includes six antennas having directivity and each antenna radiates radio waves in different directions by 60 degrees, six directional sectors C1 to C6 are provided. Is formed. An omnidirectional sector is a sector formed by, for example, radiating radio waves in one direction from one antenna. For example, as shown in FIG. 2B, the omnidirectional sector C10 is formed by radiating radio waves from the antenna 201 in all directions.

  The RNC 300 includes a communication control unit 301 and a position specifying unit 302. The communication control unit 301 is a part that performs communication connection with the mobile device 100 via the BTS 200, for example, a portion that performs communication connection processing based on communication connection processing based on outgoing call processing or incoming call processing from the mobile device 100. It is.

  The position specifying unit 302 specifies the position in the sector of the mobile device 100 using the PRACH-PD positioning method (first positioning). The location of the mobile device 100 in the sector using the PRACH-PD positioning method is executed when the mobile device 100 performs communication connection (for example, connection processing at the time of outgoing call, incoming call, or handover). The Specifically, as illustrated in FIG. 3, the position specifying unit 302 transmits a propagation time (first positioning position information) until the PRACH signal transmitted to the mobile device 100 is returned by the mobile device 100 and reaches the BTS 200. , Propagation time information). The measurement result of the propagation time is every predetermined unit time obtained based on the phase period of the radio wave.

  As shown in FIG. 4, the position specifying unit 302 obtains the positioning time when the position of the mobile device 100 is measured using the PRACH-PD positioning method, the sector identifier of the sector where the mobile device 100 is located, and the propagation time. The associated information is calculated as delay information (first correspondence information). This delay information can be calculated regardless of whether the sector is a directional sector or an omnidirectional sector. Note that the distance between the antenna 201 and the mobile device 100 can be obtained from the propagation time, and it can be said that the propagation time represents the position of the mobile device 100. That is, it can be said that the delay information is position information indicating the position of the mobile device 100.

  In addition to calculating the delay information using the PRACH-PD positioning method, the position specifying unit 302 also uses PRACH-PD position information (using the coordinates of the mobile device 100 as information for specifying the location of the mobile device 100). First correspondence information) is calculated. Specifically, when the sector is a directional sector, as illustrated in FIG. 5, the position specifying unit 302 assumes that the mobile device 100 is located on a straight line L along the directivity direction of the radio wave. Then, the position P1 on the straight line L obtained based on the separation distance is calculated as the position of the mobile device 100 (first positioning position information). Then, as shown in FIG. 6, the position specifying unit 302 is on a straight line L obtained based on the positioning time when the position of the mobile device 100 is measured, the sector identifier of the sector where the mobile device 100 is located, and the separation distance. PRACH-PD position information that associates the coordinates (latitude / longitude) of the mobile device 100 is calculated.

  When the sector is an omni-directional sector, the position specifying unit 302 assumes that the mobile device 100 is located at the position P10 of the antenna 201 as illustrated in FIG. Then, as shown in FIG. 7, the position specifying unit 302 determines the positioning time when the position of the mobile device 100 is positioned, the sector identifier of the sector where the mobile device 100 is located, and the coordinates (latitude / longitude) of the mobile device 100. The PRACH-PD position information in which (the position of the antenna 201) is associated is calculated. When the sector is an omni-directional sector, almost all the positions of the mobile device 100 that are positioned using the PRACH-PD positioning method coincide with the position of the antenna 201.

  The position specifying unit 302 outputs the calculated delay information and PRACH-PD position information to the exchange 400 through the communication control unit 301.

  The exchange 400 includes a communication control unit 401 and a location information management unit 402. Similar to the communication control unit 301 of the RNC 300, the communication control unit 401 is a part that performs communication connection processing.

  The position information management unit 402 stores the delay information and the PRACH-PD position information transmitted from the position specifying unit 302. As shown in FIGS. 4, 6, and 7, delay information (see FIG. 4) and PRACH-PD position information (see FIGS. 6 and 7) are shown in separate tables. The time and the coordinate information of the PRACH-PD position information may be held in the same table.

  As shown in FIG. 8, the GPS location management unit 501 associates the time when GPS coordinate information is acquired by the mobile device 100, the sector identifier of the sector where the mobile device 100 is located, and the latitude and longitude. The associated result is stored as GPS position information.

  The position collection unit 502 collects GPS position information stored in the GPS position management unit 501, delay information stored in the position information management unit 402, and PRACH-PD position information. Then, the position collection unit 502 outputs the collected information to the main power area map creating apparatus 600 at a predetermined timing or in response to a request from the main power area map creating apparatus 600.

  The main power area map creation device 600 creates a main power area map related to a sector. As shown in FIG. 9, the main power area map creation device 600 includes a GPS position information acquisition unit (second correspondence information acquisition unit) 601, a delay position information acquisition unit (first correspondence information acquisition unit) 602, and an operating sector determination. Unit 604, facility change detection unit 605, facility information storage unit 606, all-sector main power area map creation unit 607, active sector main power area map creation unit 608, main power area map storage unit 609, main power area creation parameter estimation unit (main power) Area estimation device) 610 and main power area map output unit 620.

  The GPS position information acquisition unit 601 acquires GPS position information (see FIG. 8) output from the position collection unit 502. The delay position information acquisition unit 602 acquires delay information and PRACH-PD position information (see FIGS. 4, 6, and 7) output from the position collection unit 502.

  The main power area creation parameter estimation unit 610 is based on at least one of GPS position information, delay information, and PRACH-PD position information acquired by the GPS position information acquisition unit 601 or the delay position information acquisition unit 602. The sector type and the main power area are calculated. In detail, as shown in FIG. 9, the main influence area creation parameter estimation section 610 includes a sector type determination section (sector type determination means) 611, an antenna position calculation section (antenna position calculation means) 612, a main influence area estimation section. (Main power area estimation means) 613, a radiation direction calculation section (radiation direction calculation means) 615, and a radiation width calculation section (radiation width calculation means) 616 are configured.

  The sector type determination unit 611 determines whether the sector formed by the BTS 200 is a directional sector or an omnidirectional sector based on the PRACH-PD position information acquired by the delay position information acquisition unit 602. To do. Hereinafter, three methods for determining the sector type will be described.

  Here, details of the PRACH-PD position information of the mobile device 100 calculated using the PRACH-PD positioning method will be described. First, a case where a sector is a directional sector and a plurality of mobile devices 100 are located in the sector will be described. In this case, since the measurement result of the propagation time is every predetermined unit time, as shown in FIG. 10, the coordinate position of each mobile device 100 in the PRACH-PD position information is the distance to the antenna 201 in each mobile device 100. Based on the above, it has a property to be any one of positions P1, P2, P3, P4, and P5 on the straight line L along the directivity direction of the radio wave.

  However, if PRACH-PD positioning is performed in a state where the mobile device 100 exists at the boundary between a plurality of sectors, the positioning point of the mobile device 100 may be approximated to the center position between the sectors. In the vicinity of the antenna 201, a plurality of sectors are adjacent at close intervals. In particular, in the vicinity of the antenna 201, as a positioning point of the mobile device 100, a position other than the straight line L (for example, the position P11, FIG. P12, P13, etc.) may be calculated. However, since the frequency with which the mobile device 100 is located at the sector boundary is low, the positioning point of the mobile device 100 is rarely calculated as a position other than the straight line L. For example, it is assumed that several hundred mobile devices 100 are measured at positions P1 to P5, respectively. In addition, it is assumed that several mobile devices 100 are measured at P11 to P13, respectively.

  Next, a case where the sector is an omnidirectional sector and a plurality of mobile devices 100 are located in the sector will be described. In this case, as shown in FIG. 11, the positioning points of each mobile device 100 in the PRACH-PD position information have a property of gathering at the position P <b> 10 of the antenna 201. However, as described above, when PRACH-PD positioning is performed in a state where the mobile device 100 exists at the boundaries of a plurality of sectors, the positioning point of the mobile device 100 may be approximated to the center position between the sectors. For this reason, although the frequency is low, positions other than the position P10 of the antenna 201 (for example, positions P14, P15, P16, and P17 in FIG. 11) may be calculated as positioning points of the mobile device 100. For example, it is assumed that the measurement is performed assuming that several hundred mobile devices 100 exist at the position P10. It is assumed that several mobile devices 100 are measured at P14 to P17, respectively.

  First, a first method of sector type determination by the sector type determination unit 611 will be described. The sector type determination unit 611 adds up PRACH-PD position information having sector identifiers for which sector types are determined for each positioning point. When the positions of the positioning points are concentrated on one point, the sector type determination unit 611 determines that the target sector is an omnidirectional sector.

  Next, a second method of sector type determination by the sector type determination unit 611 will be described. The sector type determination unit 611 adds up PRACH-PD position information having sector identifiers for which sector types are determined for each positioning point. Then, the sector type determination unit 611 counts the number of terminals (number of signals) of the mobile device 100 for each positioning point. Then, the sector type determination unit 611 obtains a straight line connecting the positioning point with the largest number of signals and the positioning point with the second largest number of signals, and the positioning points are equal to or greater than a predetermined threshold (for example, three points) on the obtained straight line. Etc.) Determine whether it exists. In addition to the positioning point on the straight line, a positioning point existing within a predetermined distance (for example, 20 m) from the straight line may be used. When it is determined that the positioning point on the straight line is equal to or greater than the predetermined threshold, the sector type determination unit 611 determines that the target sector is a directional sector. In the case of a directional sector, this utilizes the property that positioning points are arranged on a straight line L along the direction of radio wave as shown in FIG.

  Next, a third method of sector type determination by the sector type determination unit 611 will be described. The sector type determination unit 611 adds up PRACH-PD position information having sector identifiers for which sector types are determined for each positioning point. For example, when the PRACH-PD position information for the directional sector C11 shown in FIG. 10 is tabulated, as shown in FIG. 12, each position P1 to P5, which is a positioning point where the mobile device 100 is positioned, is displayed. For each of P11 to P13, the number of terminals (number of signals) of the mobile device 100 is tabulated. The latitude of the position P2 is Y2, the longitude is X2, and the number of signals is 400. The latitude of the position P5 is Y5, the longitude is X5, and the number of signals is 200.

  Further, for example, when the PRACH-PD position information for the omnidirectional sector C12 shown in FIG. 11 is tabulated, as shown in FIG. 13, each position P10, which is a positioning point where the mobile device 100 is positioned, For each of P14 to P17, the number of terminals (number of signals) of the mobile device 100 is tabulated. The latitude of the position P10 is Y10, the longitude is X10, and the number of signals is 400. The latitude of the position P15 is Y15, the longitude is X15, and the number of signals is 2.

  Next, the sector type determination unit 611 extracts two positioning points with the larger number of signals in order from the number of signals for each positioning point in the aggregated PRACH-PD position information. For example, in the total result of the PRACH-PD position information for the directional sector C11 shown in FIG. 12, the position P2 (X2, Y2) and the position P5 (X5, Y5) are extracted as two positioning points with a large number of signals. Is done. Further, for example, in the total result of the PRACH-PD position information for the omnidirectional sector C12 shown in FIG. 13, the position P10 (X10, Y10) and the position P15 (X15, Y15) are two positioning points with a large number of signals. Is extracted.

  Next, the sector type determination unit 611 compares the number of signals of the extracted two positioning points. As a result of the comparison, the sector type determination unit 611 determines that the sector is an omnidirectional sector when the difference in the number of signals at the two positioning points is larger than a predetermined value, or when the ratio is smaller than the predetermined value, If the difference in the number of signals is smaller than the predetermined value or the ratio is larger than the predetermined value, the sector is determined to be a directional sector. This is based on the fact that in the case of the omnidirectional sector, the positioning points of almost all the mobile devices 100 are the positions of the antenna 201 and are rarely positioned at positions other than the antenna 201. Further, in the case of the directional sector, the fact that the mobile device 100 is positioned relatively evenly at each position on the straight line in the direction of the radio wave is utilized.

  For example, the number of signals of two positioning points extracted from the PRACH-PD position information for the directional sector C11 shown in FIG. 12 is the number of signals 400 at the position P2 and the number of signals 200 at the position P5. . Further, the number of signals at the two positioning points extracted from the PRACH-PD position information for the omnidirectional sector C12 shown in FIG. 13 is the number of signals 400 at the position P10 and the number of signals 2 at the position P15. . Therefore, the sector type determination unit 611 determines the sector C12 (see FIG. 11) having a large difference in the number of signals as a non-directional sector as shown in FIG. 13, and the difference in the number of signals as shown in FIG. A small sector C11 (see FIG. 10) is determined as a directional sector.

  Next, the flow of processing related to the third method in the sector type determination unit 611 will be described using the flowchart shown in FIG. First, the sector type determination unit 611 acquires PRACH-PD position information having a sector identifier for which the sector type is determined, and totals the acquired PRACH-PD position information for each positioning point (step S101). Next, the sector type determination unit 611 extracts two positioning points with the larger number of signals in order from the number of signals for each positioning point in the aggregated PRACH-PD position information (step S102). When two positioning points cannot be extracted (step S102: NO), the sector type determination unit 611 determines that the sector to be determined is an omnidirectional sector (step S105).

  On the other hand, when two positioning points can be extracted (step S102: YES), the sector type determination unit 611 determines whether or not the difference in the number of signals at the two extracted positioning points is equal to or greater than a predetermined value (step S102). S103). In this determination, the difference in the number of signals is greater than or equal to a predetermined value when the smaller signal number of the two signals is less than a predetermined threshold value (for example, 5%) with respect to the larger signal number. It can also be judged as being.

  When the difference in the number of signals is equal to or greater than the predetermined value (step S103: YES), the sector type determination unit 611 determines that the sector to be determined is an omnidirectional sector (step S105). On the other hand, if the difference in the number of signals is not greater than or equal to the predetermined value (step S103: NO), the sector type determination unit 611 determines that the sector to be determined is a directional sector (step S104).

  As described above, the sector type determination unit 611 determines the sector type for each sector using the PRACH-PD position information. Further, among the above-described three methods for determining the sector type, a predetermined two or more methods may be used in combination. That is, the sector type determination unit 611 determines whether or not the positions of the respective positioning points of the PRACH-PD position information are concentrated at one point as in the above-described first method of sector type determination. The target sector is determined to be an omnidirectional sector. When the position of the positioning point is not concentrated on one point, the sector type determination unit 611 determines the positioning point with the largest number of signals and the positioning point with the second largest number of signals as in the second method of sector type determination described above. If a straight line to be connected is obtained, it is determined whether or not the positioning point is determined to be greater than or equal to a predetermined threshold (for example, 3 or more) on the obtained straight line, Is determined to be an omnidirectional sector. When it is determined that the positioning point is equal to or greater than the predetermined threshold, the sector type determination unit 611 uses the second number of signals at the positioning point with the largest number of signals as in the third method of sector type determination described above. If the difference from the number of signals at a large number of positioning points is larger than a predetermined value, or if the ratio is smaller than a predetermined value, it is determined that the target sector is an omnidirectional sector. On the other hand, when the difference between the number of signals at the positioning point with the largest number of signals and the number of signals at the positioning point with the second largest number of signals is smaller than a predetermined value, or when the ratio is larger than the predetermined value, the sector type determination unit 611 The target sector is determined to be an omnidirectional sector. In this case, the accuracy of determining the sector type can be improved. The sector type determination by the sector type determination unit 611 is not essential, and the sector type can be stored in the facility information storage unit 606 in advance.

  The antenna position calculation unit 612 calculates the position of the antenna 201 based on the distribution of the mobile device 100 obtained from the PRACH-PD position information. Note that the antenna position calculation unit 612 uses a different method to determine the position of the antenna 201 based on whether the sector formed by the BTS 200 for which the position of the antenna 201 is to be calculated is a directional sector or an omnidirectional sector. Estimate the position. As the sector type, the sector type determined by the sector type determining unit 611 can be used. Hereinafter, as the first to seventh methods, a method for calculating the position of the antenna 201 when the sector formed by the BTS 200 is a directional sector will be described, and as the eighth and ninth methods, omnidirectionality will be described. A method for calculating the position of the antenna 201 in the case of a sector will be described.

  First, a first method for calculating the position of the antenna 201 will be described. When the sector formed by the BTS 200 for which the antenna 201 is to be calculated is a directional sector, the antenna position calculation unit 612 delays the PRACH-PD position information for the BTS 200 for which the antenna 201 is to be calculated. Obtained from the position information obtaining unit 602. Note that identification information for identifying the BTS is added to the PRACH-PD position information in addition to the information shown in FIG. 6, for example.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain BTS for each positioning point, and maps it as shown in FIG. Here, it is assumed that the antenna 201 of the BTS 200 to be calculated includes three antennas having directivity, and three sectors C21, C22, and C23 are formed. Here, the position of the antenna 201 of the BTS B1 is obtained.

  When the antenna 201 has directivity, when PRACH-PD positioning is performed, the positioning points of the mobile device 100 are arranged along the direction of radio wave for each sector as described above. FIG. 15 (a) shows the positioning points of each mobile device 100 associated with sector identifiers C21 to C23 for BTS B1.

  Next, the antenna position calculation unit 612 obtains straight lines L21, L22, and L23 that pass through the mapped positioning points for each sector, as shown in FIG. For example, when the straight line L21 is obtained, it can be obtained based on two positioning points having a large number of signals among the positioning points of the mobile device 100 for the sector C21. And the antenna position calculation part 612 calculates the position of intersection A1 of the straight lines L21-L23 as the position of the antenna 201 of BTS B1.

  Next, a second method for calculating the position of the antenna 201 will be described. Here, it is assumed that the BTS 200 for which the position of the antenna 201 is to be calculated forms six sectors as in FIG. When the sector formed by the BTS 200 for which the antenna 201 is to be calculated is a directional sector, the antenna position calculation unit 612 delays the PRACH-PD position information for the BTS 200 for which the antenna 201 is to be calculated. Obtained from the position information obtaining unit 602.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain BTS for each positioning point, and maps it as shown in FIG. Here, it is assumed that the antenna 201 of the BTS 200 to be calculated includes six antennas having directivity, and six sectors C31, C32, C33, C34, C35, and C36 are formed. Here, the position of the antenna 201 of the BTS B2 is obtained. FIG. 16 (a) shows the positioning points of each mobile device 100 associated with sector identifiers C31 to C36 for BTS B2.

  Next, as shown in FIG. 16B, the antenna position calculation unit 612 obtains straight lines L31, L32, L33, L34, L35, and L36 passing through the mapped positioning points for each sector. When this straight line is obtained, it can also be obtained based on two positioning points having a large number of signals, as in the first method described above.

  And the antenna position calculation part 612 calculates | requires the intersection in the predetermined two combination of the straight lines L31-L36, respectively. Specifically, for a sector whose radio wave radiation direction is about 180 degrees (for example, about 180 degrees ± 5 degrees), the intersection point of the straight line passing through the positioning point is calculated. The case where it deviates from the position of an antenna is considered. Therefore, for a sector whose radio wave radiation direction differs by about 180 degrees, the intersection of straight lines passing through the positioning point is not obtained. Therefore, the intersection of the straight line L31 and the straight line L32, the intersection of the straight line L31 and the straight line L33, the intersection of the straight line L31 and the straight line L35, the intersection of the straight line L31 and the straight line L36, the intersection of the straight line L32 and the straight line L33, and the straight line L32 Intersection with straight line L34, Intersection with straight line L32 and straight line L36, Intersection with straight line L33 and straight line L34, Intersection with straight line L33 and straight line L35, Intersection with straight line L34 and straight line L35, Intersection with straight line L34 and straight line L36 , And the intersection of the straight line L35 and the straight line L36, respectively.

  Then, the antenna position calculation unit 612 calculates the average position A2 of the obtained intersections as the position of the antenna 201 of BTS B2. In addition, when the distance between the average position of each intersection and the position of each intersection is within a predetermined distance (for example, about several meters), the obtained average position of each intersection can be calculated as the antenna position. .

  Next, a third method for calculating the position of the antenna 201 will be described. When calculating the position of the antenna 201 of the BTS 200 that forms a certain directional sector, the antenna position calculation unit 612 acquires PRACH-PD position information having a sector identifier that is a calculation target of the position of the antenna 201, and acquires delayed position information. From the unit 602.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain sector for each positioning point, and calculates the number of signals for each positioning point as in FIG. When the aggregation results are mapped, as shown in FIG. 17, the positioning points are arranged in the direction of the radio wave.

  And the antenna position calculation part 612 calculates | requires the relationship between the distance from an origin, and the number of signals in each positioning point by making one end point of the alignment direction of a positioning point into the origin among each positioning point shown in FIG. This is represented by a graph in FIG. In FIG. 18, the horizontal axis represents the distance from the origin, and the vertical axis represents the number of signals.

  In the example shown in FIG. 18, a portion with a large number of signals is biased toward the origin position side. Therefore, the antenna position calculation unit 612 calculates the position of the positioning point as the origin among the positioning points as the position of the antenna 201. This uses a wireless characteristic that the closer to the BTS 200, the more easily the mobile device 100 belongs to the BTS 200 (the number of signals increases).

  In FIG. 18, when the portion with a large number of signals is biased to the side opposite to the origin (the side far from the origin), the antenna position calculation unit 612 determines the position of the other end point in the positioning point arrangement direction. Is calculated as the position of the antenna 201.

  Next, a fourth method for calculating the position of the antenna 201 will be described. When calculating the position of the antenna 201 of the BTS 200 that forms a certain directional sector, the antenna position calculation unit 612 acquires PRACH-PD position information having a sector identifier that is a calculation target of the position of the antenna 201, and acquires delayed position information. From the unit 602.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain sector for each positioning point, and calculates the number of signals for each positioning point as in FIG. When the aggregation results are mapped, as shown in FIG. 19A, the positioning points are arranged in the direction of the radio wave. Here, it is assumed that there are 11 positioning points (positioning points T1 to T11).

  And the antenna position calculation part 612 extracts two positioning points with many signals sequentially from the positioning points shown in FIG. For example, it is assumed that the number of signals at the positioning point T2 is 400, the number of signals at the positioning point T5 is 200, and the positioning points T2 and T5 are extracted. As shown in FIG. 19B, the antenna position calculation unit 612 obtains a straight line L40 passing through the positioning point T2 and the positioning point T5.

  As shown in FIGS. 19B and 19C, the antenna position calculation unit 612 has positioning points T1 positioned at both ends of the straight line L40 among the positioning points within a predetermined distance from the obtained straight line L40. , T11. Thereby, as described above, even when the mobile device 100 is positioned at a position different from the directivity direction of the radio wave when performing PRACH-PD positioning, these positioning points are determined from the antenna position candidates. This can be excluded, and the calculation accuracy of the antenna position can be improved.

  Next, the antenna position calculation unit 612 calculates the number of signals at each positioning point, as shown in FIG. Then, the antenna position calculation unit 612 calculates, as the position of the antenna 201, the position of the positioning point on the side where the portion with a large number of signals is biased among the positioning points T1 and T11. In the example shown in FIG. 19D, it is assumed that a portion with a large number of signals is biased toward the positioning point T1. Therefore, the antenna position calculation unit 612 calculates the position of the positioning point T1 as the position of the antenna 201.

  In addition, the position of the antenna 201 can be calculated by the following method in addition to calculating the position of the positioning point on the side where the portion having a large number of signals is biased as the position of the antenna 201. Specifically, for example, when the positioning points T1 to T11 are obtained as shown in FIG. 19C, the antenna position calculation unit 612 calculates the total number of signals at the positioning points T1 to T5 (hereinafter referred to as “S ( T1 to T5) ”) and the total number of signals at the positioning points T7 to T11 (hereinafter referred to as“ S (T7 to T11) ”). That is, it is only necessary that the number of positioning points in the group including the positioning point T1 is the same as the number of positioning points in the group including the positioning point T11 (for example, the positioning points shown in FIG. If it is up to T10, it can be divided into positioning points T1 to T5 and positioning points T6 to T10.) Then, the antenna position calculation unit 612 determines the positioning points T1, T11 when the ratio of the total value of the number of signals of one group and the total value of the number of signals of the other group is equal to or greater than a predetermined threshold. Among them, the positioning point belonging to the group having the larger total number of signals can be calculated as the position of the antenna 201. Here, the ratio is equal to or greater than a predetermined threshold, for example, when the total value of the number of signals in one group is greater than twice the total value of the number of signals in the other group.

  Next, a fifth method for calculating the position of the antenna 201 will be described. When calculating the position of the antenna 201 that forms a certain directional sector, the antenna position calculation unit 612 obtains the PRACH-PD position information having the sector identifier for which the position of the antenna 201 is calculated, as the delay position information acquisition unit 602. Get from.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain sector for each positioning point, and calculates the number of signals for each positioning point as in FIG. When the aggregation results are mapped, as shown in FIG. 20, the positioning points are arranged in the direction of the radio wave. Furthermore, as described above, it is assumed that the positioning points T20 to T25 are measured at a place other than the straight line L41 that is the direction of the radio wave by the PRACH-PD positioning.

  Positioning points T <b> 20 to T <b> 25 that are measured at locations other than on the straight line L <b> 41 are particularly frequently positioned near the installation position of the antenna 201. Therefore, the antenna position calculation unit 612 calculates the position of the antenna 201 based on the distribution of the positioning points T20 to T25. In the example shown in FIG. 20, for example, at a predetermined position (for example, an end point of a positioning point aligned on the straight line L41) in an area R1 surrounded by a circle with a predetermined radius (for example, 1 km) centering on the end point on the straight line L41. Calculation is made assuming that the antenna 201 exists.

  Next, a sixth method for calculating the position of the antenna 201 will be described. When calculating the position of the antenna 201 that forms a certain directional sector, the antenna position calculation unit 612 obtains the PRACH-PD position information having the sector identifier for which the position of the antenna 201 is calculated, as the delay position information acquisition unit 602. Get from. Here, each positioning point indicated by the acquired PRACH-PD position information includes positioning points arranged on a straight line L41 and positioning points T20 to T25 located on other than the straight line L41 as shown in FIG. . Further, out of the positioning points arranged on the straight line L41, the positioning points at both ends are set as positioning points T26 and T27, respectively.

  The antenna position calculation unit 612 extracts positioning points T20 to T25 that are measured at places other than the straight line L41. And the antenna position calculation part 612 calculates | requires the sum total of the distance of the positioning point T26 located in the edge part of the positioning point on a straight line L41, and each extracted positioning point T20-T25. Similarly, the antenna position calculation unit 612 calculates the sum of the distances between the positioning points T27 located at the end portions of the positioning points arranged on the straight line L41 and the extracted positioning points T20 to T25. The antenna position calculation unit 612 has a positioning point having a smaller total sum of distances from the positioning points T20 to T25 among the positioning points T26 and T27 at both ends of the positioning point arranged on the straight line L41 (in the example of FIG. 20, the positioning point). T26) is calculated as the position of the antenna 201.

  Next, a seventh method for calculating the position of the antenna 201 will be described. When calculating the position of the antenna 201 that forms a certain directional sector, the antenna position calculation unit 612 obtains the PRACH-PD position information having the sector identifier for which the position of the antenna 201 is calculated, as the delay position information acquisition unit 602. Get from.

  The antenna position calculation unit 612 aggregates the acquired PRACH-PD position information for a certain sector for each positioning point, and calculates the number of signals for each positioning point as in FIG. When the aggregation results are mapped, as shown in FIG. 21A, the positioning points are arranged in the directivity direction of the radio wave. Furthermore, it is assumed that positioning points have been measured at locations other than on a straight line along the direction of radio wave by PRACH-PD positioning.

  And the antenna position calculation part 612 extracts two positioning points with many signals sequentially from the positioning points shown in FIG. For example, it is assumed that the number of signals at the positioning point T32 is 400, the number of signals at the positioning point T35 is 200, and the positioning points T32 and T35 are extracted. As shown in FIG. 21B, the antenna position calculation unit 612 obtains a straight line L42 passing through the positioning point T32 and the positioning point T35.

  As shown in FIGS. 21 (b) and 21 (c), the antenna position calculation unit 612 has positioning points T31 positioned at both ends of the straight line L40 among positioning points within a predetermined distance from the obtained straight line L42. , T39. Thereby, as described above, even when the mobile device 100 is positioned at a position different from the directivity direction of the radio wave when performing PRACH-PD positioning, these positioning points are determined from the antenna position candidates. This can be excluded, and the calculation accuracy of the antenna position can be improved.

  Next, as shown in FIG. 21D, the antenna position calculation unit 612 calculates the dispersion of the positioning points around the positioning fixed point T31 and the dispersion of the positioning points around the positioning fixed point T39. Then, the antenna position calculation unit 612 calculates the position of the positioning point T31 located on the side with large dispersion as the position of the antenna 201. This is based on the fact that positioning points measured at a place other than on the straight line L42 are gathered particularly in the vicinity of the installation position of the antenna 201.

  Of the first to seventh methods for calculating the position of the antenna 201 described above, the position of the antenna 201 can be calculated by combining two or more predetermined methods. In this case, the calculation accuracy of the position of the antenna 201 can be improved.

  Next, an eighth method for calculating the position of the antenna 201 will be described. The antenna position calculation unit 612, when the sector formed by the BTS 200 that is the calculation target of the antenna 201 is an omni-directional sector, stores the PRACH-PD position information including the sector identifier that is the calculation target of the antenna 201, Obtained from the delay position information obtaining unit 602. In this case, almost all the positions of the mobile device 100 coincide with the position of the antenna 201. Therefore, the antenna position calculation unit 612 calculates the positioning point with the highest signal density as the position of the antenna 201.

  Next, a ninth method for calculating the position of the antenna 201 will be described. The antenna position calculation unit 612, when the sector formed by the BTS 200 that is the calculation target of the antenna 201 is an omni-directional sector, stores the PRACH-PD position information including the sector identifier that is the calculation target of the antenna 201, Obtained from the delay position information obtaining unit 602.

  The antenna position calculation unit 612 has a signal number in which the ratio between the positioning point with the highest signal density (hereinafter referred to as “first high-density point”) and the number of signals at the first high-density point is 5% or more. The positioning point is calculated as the position of the antenna 201. Here, there are sectors that have the same sector identifier and are operated by different antennas. In such a case, by using the eighth method, a plurality of antenna positions are estimated for one omnidirectional sector. For example, when drawing a power area, the power area is centered on each antenna position. Draw.

  As described above, the antenna position calculation unit 612 calculates the position of the antenna 201. Note that the calculation of the position of the antenna 201 by the antenna position calculation unit 612 is not essential, and the position of the antenna 201 can be stored in the facility information storage unit 606 in advance.

  The main power area estimation unit 613 calculates the distribution of the mobile device 100 based on the delay information, the PRACH-PD position information, or the GPS position information, and determines the estimation target based on the calculated distribution and the position of the antenna 201. The radius of the main power area in which the sector to be the main power is estimated. The main influence area estimation unit 613 estimates the radius of the main influence area by the following method according to the acquired information (delay information, PRACH-PD position information, or GPS position information).

  First, the case where the radius of the main power area where the estimation target sector is the main power is estimated using the PRACH-PD position information will be described. Here, it is assumed that the sector to be estimated is a directional sector. When estimating the radius of the main power area in which a certain sector is the main power, the main power area estimation unit 613 determines the PRACH-PD position information (for example, FIG. 6) having the sector identifier of the sector for which the main power area is estimated. Is acquired from the delay position information acquisition unit 602. Here, a case where the radius of the main power area is calculated for the sector with the sector identifier C8 will be described.

  Then, as shown in FIG. 22, the main power area estimation unit 613 adds up the acquired PRACH-PD position information for each latitude and longitude, and calculates the number of signals for each latitude and longitude. The PRACH-PD position information used here preferably uses only information about positioning points arranged on a straight line along the radio wave radiation direction. Further, as a straight line along the radio wave radiation direction, a straight line connecting a positioning point with the largest number of signals and a positioning point with the second largest number of signals can be used. Next, the main power area estimation unit 613 acquires position information of the antenna 201 that forms a sector that is an estimation target of the main power area. Here, the position information of the antenna 201 calculated by the antenna position calculation unit 612 is used, or when the position information of the antenna 201 is stored in the facility information storage unit 606 in advance, the stored position information of the antenna 201 is used. Can be used.

  Next, the main power area estimation unit 613 determines the distance between the antenna 201 and the mobile device 100 based on the latitude / longitude (latitude / longitude shown in FIG. 22) of the positioning result of the mobile device 100 and the position of the antenna. Calculate the distance. Then, as shown in FIG. 23, the main power area estimation unit 613 calculates the sector identifier, the distance between the antenna 201 and the mobile device 100, and the number of signals from the PRACH-PD position information collected for each latitude and longitude. The associated intermediate table is calculated.

  Next, as shown in FIG. 24, the main power area estimation unit 613 obtains a signal density distribution by distance, with the horizontal axis representing the distance and the vertical axis representing the number of signals, as shown in FIG. Then, the main influence area estimation unit 613 obtains the distance D from the antenna position (origin position) to the portion where the accumulated density of the number of signals is 90% in the obtained signal density distribution by distance. The main power area estimation unit 613 estimates the obtained distance D as the radius of the main power area where the sector C8 (sector identifier C8 sector, hereinafter the same) is the main power.

  Next, the flow of processing when the radius of the main power area in which the sector to be estimated is the main power is estimated using the PRACH-PD position information will be described with reference to FIG. First, the main power area estimation unit 613 acquires, for example, the PRACH-PD position information illustrated in FIG. 6 from the delay position information acquisition unit 602 (step S201: correspondence information acquisition step). Next, the main power area estimation unit 613 totals the acquired PRACH-PD position information for each latitude and longitude, and calculates the number of signals for each latitude and longitude (see FIG. 22) (step S202).

  Next, the main power area estimation unit 613 acquires position information of the antenna 201 that forms a sector that is an estimation target of the main power area (step S203). Then, the main power area estimation unit 613 calculates the distance between the antenna 201 and the mobile device 100 (step S204) and associates the sector identifier, the distance, and the number of signals (see FIG. 23). Is calculated (step S205).

  Next, the main power area estimation unit 613 calculates a signal density distribution by distance (see FIG. 24) from the intermediate table (step S206). Then, the main power area estimation unit 613 calculates the distance D to the portion where the cumulative density of signals is 90%, and estimates the calculated distance D as the radius of the main power area where the sector C8 is the main power ( Step S207: main power area estimation step).

  As described above, when the acquired PRACH-PD position information is information on the directional sector, the main power area estimation unit 613 has a main power area in which the sector is the main power based on the PRACH-PD position information. Can be obtained.

  Next, a case where the radius of the main power area in which the estimation target sector is the main power is estimated using the delay information will be described. When estimating the radius of the main power area in which a certain sector is the main power, the main power area estimation unit 613 has delay information (for example, the information shown in FIG. 4) having the sector identifier of the sector for which the main power area is estimated. ) Is acquired from the delay position information acquisition unit 602. Here, a case where the radius of the main power area is calculated for the sector with the sector identifier C7 will be described. The sector here may be either a directional sector or an omnidirectional sector.

  Then, as shown in FIG. 26, the main power area estimation unit 613 aggregates the acquired delay information for each propagation time, and calculates the number of signals for each propagation time. Next, the main power area estimation unit 613 calculates the distance between the antenna 201 and the mobile device 100 based on the propagation time, and uses the sector identifier as shown in FIG. And an intermediate table in which the distance between the antenna 201 and the mobile device 100 is associated with the number of signals.

  Next, as shown in FIG. 24, the main power area estimation unit 613 obtains a signal density distribution by distance, with the horizontal axis representing the distance and the vertical axis representing the number of signals, as shown in FIG. Then, the main influence area estimation unit 613 obtains the distance D from the antenna position (origin position) to the portion where the accumulated density of the number of signals is 90% in the obtained signal density distribution by distance. However, a value of a predetermined cumulative density can be used in addition to the 90% value used as the cumulative density. The main power area estimation unit 613 estimates the obtained distance D as the radius of the main power area where the sector C7 is the main power.

  Next, the flow of processing when estimating the radius of the main power area in which the sector to be estimated is the main power using the delay information will be described with reference to FIG. First, the main power area estimation unit 613 acquires, for example, the delay information illustrated in FIG. 4 from the delay position information acquisition unit 602 (step S301: correspondence information acquisition step). Next, the main power area estimation unit 613 aggregates the acquired delay information for each propagation time, and calculates the number of signals for each propagation time (see FIG. 26) (step S302).

  Next, the main power area estimation unit 613 calculates an intermediate table (see FIG. 27) in which the sector identifier, the distance, and the number of signals are associated (step S303). Then, the main influence area estimation unit 613 calculates a signal density distribution by distance (see FIG. 24) from the intermediate table (step S304).

  Next, the main power area estimation unit 613 calculates the distance D to the portion where the cumulative density of signals is 90%, and estimates the calculated distance D as the radius of the main power area where the sector C7 is the main power. (Step S305: main power area estimation step).

  As described above, the main power area estimation unit 613 can obtain the radius of the main power area in which the directional sector and the non-directional sector are the main power based on the acquired delay information.

  Next, the case where the radius of the main influence area where the sector to be estimated is the main influence is estimated using the GPS position information will be described. When estimating the radius of the main power area in which a certain sector is the main power, the main power area estimation unit 613 has GPS position information (for example, as shown in FIG. 8) having the sector identifier of the sector for which the main power area is estimated. Information) is acquired from the GPS position information acquisition unit 601. Here, a case where the radius of the main power area is calculated for the sector with the sector identifier C20 will be described.

  The main power area estimation unit 613 acquires position information of the antenna 201 that forms a sector that is an estimation target of the main power area. Here, the position information of the antenna 201 calculated by the antenna position calculation unit 612 is used, or when the position information of the antenna 201 is stored in the facility information storage unit 606 in advance, the stored position information of the antenna 201 is used. Can be used.

  The main power area estimation unit 613 calculates the distance between the antenna 201 and the mobile device 100 based on the latitude / longitude (latitude / longitude shown in FIG. 8) of the GPS position information and the position of the antenna. Then, as shown in FIG. 24, the main influence area estimation unit 613 obtains a signal density distribution by distance with the horizontal axis representing the distance and the vertical axis representing the number of signals, as shown in FIG.

  Then, the main influence area estimation unit 613 obtains the distance D from the antenna position (origin position) to the portion where the accumulated density of the number of signals is 90% in the obtained signal density distribution by distance. However, a value of a predetermined cumulative density can be used in addition to the 90% value used as the cumulative density. The main power area estimation unit 613 estimates the obtained distance D as the radius of the main power area where the sector C20 is the main power.

  For example, when the sector C20 is a directional sector, as shown in FIG. 29A, a sector area R2 having a radius D from the position of the antenna 201 is a main power area where the sector C20 is the main power. Presumed. Further, when the sector C20 is an omnidirectional sector, as shown in FIG. 29B, an area R3 in a circle having a radius D from the position of the antenna 201 is a main power area where the sector C20 is the main power area. Is estimated as

  The radiation direction calculation unit 615 obtains the distribution of the mobile device 100 based on the PRACH-PD position information, and calculates the radiation direction of the radio wave radiated from the antenna 201 based on the obtained distribution. Note that the radiation direction calculation unit 615 does not need to perform a process of calculating the radiation direction of the radio wave for the omnidirectional sector. Whether or not the sector is an omnidirectional sector can be determined by using the determination result of the sector type determination unit 611 and the sector type information stored in the facility information storage unit 606. Here, the radiation direction of the radio wave radiated from the antenna 201 forming the directional sector is calculated.

  Specifically, when the radiation direction calculation unit 615 calculates the radiation direction of the radio wave of the antenna 201 that forms a certain sector, the PRACH-PD position information having the sector identifier of the sector formed by the antenna 201 that is the calculation target (For example, information shown in FIG. 6) is acquired from the delay position information acquisition unit 602. Here, a case will be described in which the radiation direction of the radio wave radiated from the antenna 201 forming the sector with the sector identifier C8 is calculated.

  In addition, the radiation direction calculation unit 615 acquires position information of the antenna that forms the sector C8. Here, the position information of the antenna 201 calculated by the antenna position calculation unit 612 is used, or when the position information of the antenna 201 is stored in the facility information storage unit 606 in advance, the stored position information of the antenna 201 is used. Can be used. When the position of the mobile device 100 and the position of the antenna 201 in the acquired PRACH-PD position information are mapped, as shown in FIG. 30A, the positioning points are arranged in the direction of the radio wave from the position of the antenna 201. Have

  The radiation direction calculation unit 615 aggregates the acquired PRACH-PD position information of the sector identifier C8 for each positioning point (for each latitude and longitude), and calculates the number of signals for each positioning point as shown in FIG.

  And the radiation | emission direction calculation part 615 extracts two in order from a positioning point with many signals among each positioning point shown in FIG. 31 in order. For example, as shown in FIG. 30B, the number of signals at the positioning point T42 (X2, Y2) is 400, and the number of signals at the positioning point T45 (X5, Y5) is 200. The positioning points T42, T45. Is extracted. Further, the radiation direction calculation unit 615 obtains a straight line L43 passing through the positioning point T42 and the positioning point T45 as shown in FIG. 30 (b).

  The radiation direction calculation unit 615 obtains the radiation direction of the radio wave radiated from the antenna 201 forming the sector C8 from the position of the antenna 201 and the straight line L43. For example, as shown in FIG. 30C, the north direction is 0 degree, and the radiation direction of the radio wave can be expressed by the angle of the straight line L43 with respect to the north direction. Note that the radiation direction calculation process by the radiation direction calculation unit 615 is not essential, and the main power area map creation process or the like can be performed using data stored in the facility information accumulation unit 606 in advance.

  The radiation width calculation unit 616 calculates the distribution of the mobile device 100 based on the GPS position information, and calculates the radiation width of the radio wave radiated from the antenna forming the predetermined sector based on the calculated distribution. Note that the emission width calculation unit 616 does not need to perform processing for calculating the emission width of radio waves for non-directional sectors. Whether or not the sector is an omnidirectional sector can be determined by using the determination result of the sector type determination unit 611 and the sector type information stored in the facility information storage unit 606. Here, the case where the radiation width of the radio wave radiated from the antenna 201 forming the directional sector is calculated will be described.

  Specifically, the radiation width calculation unit 616 obtains GPS position information (for example, information shown in FIG. 8) having the sector identifier of the sector formed by the radio wave for which the radiation width is to be calculated, as a GPS position information acquisition unit 601. Get from. Here, a case where the radiation width of the radio wave forming the sector of the sector identifier C20 is calculated will be described.

  The radiation width calculation unit 616 acquires position information of the antenna 201 that forms the sector C20 formed by the radio wave that is the target of radiation width calculation. Here, the position information of the antenna 201 calculated by the antenna position calculation unit 612 can be used, or when the position information of the antenna 201 is stored in advance, the stored position information of the antenna 201 can be used. .

  When the position of the mobile device 100 and the position of the antenna 201 in the acquired GPS position information are mapped, as shown in FIG. 32 (a), the mobile device has a predetermined spread from the position of the antenna 201 to the radio wave radiation direction side. 100 positioning points are scattered.

  Further, the radiation width calculation unit 616 creates a plurality of radial sector regions S centered on the position of the antenna 201 as shown in FIG. As the radius of the fan-shaped region S, a predetermined value set in advance or the radius of the main power area estimated by the main power area estimation unit 613 can be used.

  Next, the radiation width calculation unit 616 includes a positioning point whose distance to the antenna 201 is less than a predetermined distance among positioning points included in the extracted GPS position information, and a positioning whose distance to the antenna 201 is equal to or greater than a predetermined distance. A positioning point other than a point is extracted. Specifically, as shown in FIG. 32B, the positioning points in the area U1 where the distance to the antenna 201 is less than the predetermined distance and the area U2 where the distance to the antenna 201 is equal to or larger than the predetermined distance are excluded. Then, the remaining positioning points are extracted. Thereby, for example, a positioning point measured by a radio wave or the like that circulates in the vicinity of the antenna 201 in a direction opposite to the radiation direction of the radio wave, which is caused by the characteristics of the antenna 201, can be excluded.

  Based on each extracted positioning point, the radiation width calculation unit 616 uses the angle with respect to the radio wave radiation direction as a reference, and the number of positioning points (signals) in each sector area S, as shown in FIG. The density distribution of the number of signals is created with the number) as the vertical axis. The radio wave radiation direction used here is a value calculated by the radiation direction calculation unit 615, or when the radio wave radiation direction is stored in the facility information storage unit 606 in advance, the stored radio wave radiation direction is used. Direction can be used.

  The emission width calculation unit 616 calculates the emission width of radio waves based on the created density distribution of the number of signals. First, a procedure for obtaining an angle range on the positive angle side with respect to 0 ° (radiation direction of radio waves) as the radio wave emission width will be described. The radiation width calculation unit 616 obtains the total value of the number of signals for each angle from the larger angle side toward the smaller angle side. Then, the radiation width calculation unit 616 obtains an angle at which the total value of the number of signals is 5% of the total number of signals. In the example shown in FIG. 32C, it is assumed that an angle of 30 ° is obtained. The radiation width calculation unit 616 sets an angle at which the total value of the number of signals thus obtained is 5% of the total number of signals as an angle boundary on the positive angle side in the radio wave radiation width. Next, a procedure for obtaining an angle range on the negative angle side with respect to 0 ° as the radio wave emission width will be described. Similarly to the above, the radiation width calculation unit 616 calculates the total value of the number of signals for each angle from the smaller angle side to the larger side. Then, the radiation width calculation unit 616 obtains an angle at which the total value of the number of signals is 5% of the total number of signals. In the example shown in FIG. 32C, it is assumed that an angle of −30 ° is obtained. The emission width calculation unit 616 sets the angle at which the total value of the number of signals obtained in this way is 5% of the total number of signals as an angle boundary on the minus angle side in the emission width of the radio wave. Then, the radiation width calculation unit 616 calculates 60 °, which is an angle range between the angle boundary 30 ° on the plus angle side and the angle boundary −30 ° on the minus angle side, as the radiation width of the radio wave.

  The radial width calculator 616 gradually approaches the positive angle side boundary and the negative angle side angle boundary toward the angle 0 degree while keeping the absolute values of the respective angles the same. The total value of the number of signals for each angle sandwiched between the angle boundary and the angle boundary on the negative angle side is 90% of the total number of signals (this 90% value is an example, and other values may be used). Each of the angle boundaries at the time can be an angle boundary on the plus angle side and an angle boundary on the minus angle side in the radiation width of the radio wave.

  In the above description, an angle that is 5% of the total number of signals is used when obtaining the total value of the number of signals. However, the angle is not limited to 5%, and other values may be used. For example, when the population calculated using the traffic information of the sector is allocated to the main power area created in the present embodiment, it is preferable to obtain the radiation width by the following method. Here, the density distribution graph shown in FIG. 32C will be described using a simplified density distribution graph (see FIG. 33). When calculating the population using traffic information in the sector, it is conceivable to perform statistical processing on the assumption that the population is uniformly distributed in the main power area. That is, in each sector area S shown in FIG. 32A, the sector area S having a large number of positioning points (having a large population) is subjected to statistical processing, etc., assuming that the number of mobile devices 100 is smaller than the actual number. The fan-shaped area S having a small number of positioning points (small population) is subjected to statistical processing and the like assuming that the number of mobile devices 100 is larger than the actual number.

  As described above, when the mobile devices 100 are uniformly distributed within a predetermined radiation width, the actual number of mobile devices 100 is increased or decreased for each sector area S. May cause errors. When calculating the radiation width so as to reduce this error, a rectangle obtained based on the total number of signals in the triangle St (see FIG. 33B for details) and the radiation width shown in FIG. It is preferable that the total number of signals in Ss (refer to FIG. 33C for details) is the same.

The rectangle Ss is formed by a horizontal axis and a vertical axis passing through the origin, a straight line Ls extending in the horizontal axis direction, and a straight line Lsk extending in the vertical axis direction. The angle indicated by the straight line Lsk (the angle indicated by the horizontal axis at the intersection of the straight line Lsk and the horizontal axis) is the boundary of the angle on one side when the radio wave radiation width is determined. Hereinafter, the angle indicated by the horizontal axis at the intersection of the straight line Lsk and the horizontal axis is referred to as an angle k. Note that the total number of signals in the rectangle Ss can be obtained by integrating the straight line Ls. The straight line Ls can be expressed by the following formulas (1) and (2).
Ls = h (d ≦ k, h = total number of signals in triangle St / k) (1)
Ls = 0 (d> k) (2)

  The triangle St is formed by a horizontal axis and a vertical axis passing through the origin, and an inclination Lt of the number of signals.

Therefore, in order to make the total number of signals in the triangle St equal to the total number of signals in the rectangle Ss, an angle k satisfying the following expression (3) is obtained.

  The angle k obtained using the equation (3) is a boundary of one side angle (here, an angle on the plus side with respect to the radio wave radiation direction) when the radiation width is determined. Similarly, using equation (3), the boundary of the angle on the other side when determining the radiation width (here, the angle on the minus side with respect to the radio wave radiation direction) is obtained. In this way, the emission width calculation unit 616 can calculate the emission width of the radio wave particularly suitable for calculating the population in the sector using the equation (3). Note that the radiation width calculation process by the radiation width calculation unit 616 is not essential, and a main power area map creation process or the like can also be performed using data stored in the facility information storage unit 606 in advance.

  The equipment change detection unit 605 includes the sector type determined by the sector type determination unit 611, the position information of the antenna 201 calculated by the antenna position calculation unit 612, the radius of the main power area estimated by the main power area estimation unit 613, The radio wave radiation direction calculated by the radiation direction calculation unit 615 and the radio wave radiation width calculated by the radiation width calculation unit 616 are acquired. Then, the equipment change detection unit 605 compares the acquired information with the information stored in the equipment information storage unit 606. If there is a change, the equipment change detection unit 605 uses the information stored in the equipment information storage unit 606. Update. Note that the equipment change detection unit 605 is effective for reducing the amount of calculation processing, but without the equipment change detection unit 605, the sector type calculation result or the like may be stored in the equipment information storage unit 606. .

  In this manner, the facility change detection unit 605 updates various information about the antenna and the sector. Therefore, various actual information about the antenna and the sector can be obtained without waiting for the operator to input the facility information (various information about the antenna and the sector) of the BTS 200.

  The facility information storage unit 606 stores information related to the facilities of the BTS 200. Specifically, the sector information, the position information of the antenna 201, the radius of the main power area whose main power is the target sector, and the antenna 201 The radiation direction and the radiation width of the radio wave radiated from are stored.

  Based on various information stored in the facility information storage unit 606, the all-sector main power area map creation unit 607 is the main power of these sectors for all sectors whose information is stored in the facility information storage unit 606. Create a power area map.

  For example, when creating the main power area map of the directional sector, the all-sector main power area map creation unit 607 corresponds to the antenna 201 corresponding to the sector to be created among the pieces of information stored in the facility information storage unit 606. , Information on the radius of the main power area, the radiation direction of the radio wave, and the radio wave emission width are extracted. Whether or not the sector is a directional sector can be determined based on the sector information stored in the facility information storage unit 606. Based on the extracted information, the all-sector main power area map creation unit 607, as shown in FIG. 34A, the antenna position is the center of the sector, the radius of the main power area is the sector radius, A fan-shaped main influence area map is created with the width of the fan as the divergence angle and the direction of the radio wave as the direction of the fan. Similarly, when creating the main power area map of the omnidirectional sector, the all-sector main power area map creation unit 607 is the antenna corresponding to the sector to be created among the pieces of information stored in the facility information storage unit 606. Information on the position 201 and the radius of the main power area is extracted. Based on the extracted information, the all-sector main power area map creation unit 607, as shown in FIG. 34 (b), a circular main power area map centered on the antenna position and having a radius of the main power area as a radius. Create The all-sector main power area map creation unit 607 creates the main power area map as a circle when the sector type is an omnidirectional sector.

  The operating sector determination unit 604 acquires PRACH-PD position information from the delay position information acquisition unit 602, and determines an operating sector. For example, first, the active sector determination unit 604 arranges PRACH-PD position information for a certain sector in time series. Then, the operating sector determination unit 604 determines whether or not the positioning interval of the PRACH-PD position information is open for a predetermined time or more, and if the positioning interval is open for a predetermined time or more, it is assumed that the sector is not operating. to decide. Further, the operating sector determination unit 604 determines whether the sector is operating using delay information, location registration information generated when the mobile device 100 performs location registration, a communication traffic log, and the like. You can also.

  The operating sector main power area map creation unit 608 acquires information on the sector that is determined to be operating by the operating sector determination unit 604. The operating sector determination unit 604 acquires the main power area map for the sector determined to be operating from the all-sector main power area map creation unit 607. The operating sector main power area map creation unit 608 geographically integrates the acquired main power area map and creates a main power area map related to the operating sector. Here, for example, a main power area map of each operating sector in Japan is created. The operating sector main power area map creation unit 608 stores the created main power area map in the main power area map storage unit 609.

  The main power area map storage unit 609 stores the main power area map created by the operating sector main power area map creation unit 608.

  The main power area map output unit 620 outputs the main power area map stored in the main power area map storage unit 609 to other devices as necessary. The main power area map output from the main power area map output unit 620 is used for facility design planning and the like.

  The present embodiment is configured as described above, and the main power area estimation unit 613 calculates the distribution of the mobile device 100 based on the delay information, the PRACH-PD position information, or the GPS position information. Based on the position 201, the radius of the main power area where the sector to be estimated is the main power is estimated. That is, since the radius of the main power area can be obtained based on the actual distribution of the mobile devices 100, the main power area of the sector can be obtained accurately. In addition, it is not necessary to perform radio wave measurement at each point, and it is not necessary to incorporate a special function for determining the sector's main power area in the mobile device itself. Can be sought.

  Further, when the delay information is acquired, the main power area estimation unit 613 can calculate the distribution of the mobile device 100 based on the delay information and can estimate the main power area of the sector based on this distribution.

  Further, when the PRACH-PD position information is acquired, the main power area estimation unit 613 calculates the distribution of the mobile devices 100 based on the PRACH-PD position information, and based on this distribution, determines the main power area of the sector. Can be estimated.

  Further, the sector type determination unit 611 determines the sector type based on the distribution of the mobile devices 100. Thereby, in addition to the main power area of a sector, the type of sector can be grasped. Therefore, the sector type can be grasped based on the distribution of the mobile devices 100 without using the information regarding the sector type included in the equipment information of the BTS 200 input by the operator or the like.

  Further, the antenna position calculation unit 612 calculates the position of the antenna based on the distribution of the mobile device 100. Thereby, in addition to the main power area of the sector, the position of the antenna can be grasped. Therefore, it is possible to grasp the position of the antenna based on the distribution of the mobile device 100 without using information regarding the position of the antenna included in the equipment information of the BTS 200 input by the worker or the like.

  Further, the radiation direction calculation unit 615 calculates the radiation direction of the radio wave based on the distribution of the mobile device 100. Thereby, in addition to the main power area of the sector, the radiation direction of the radio wave can be grasped. Accordingly, it is possible to grasp the radiation direction of radio waves based on the distribution of the mobile devices 100 without using information regarding the radio wave radiation direction included in the equipment information of the BTS 200 input by an operator or the like.

  Further, the radiation width calculation unit 616 calculates the radiation width of the radio wave radiated from the antenna 201 based on the distribution of the mobile device 100 obtained using the GPS position information. In this way, the radio wave emission width can be accurately calculated using the GPS position information.

  Next, a modified example when the all-sector main power area map creation unit 607 creates a main power area map will be described. First, the all-sector main power area map creation unit 607 obtains the position of the antenna 201 and the radius of the main power area corresponding to the sector to be created among the pieces of information stored in the facility information storage unit 606. Further, the all-sector main power area map creation unit 607 acquires the sector radio wave reach of the sector for which the main power area map is to be created. This sector radio wave reachable range is a theoretical radio wave reachable range, and can be obtained by an existing method such as the Okumura-Kashiwa method. This sector radio wave reachable range can be stored in advance in the facility information accumulation unit 606, and the sector radio wave reachable range may be acquired from the facility information accumulation unit 606.

  Next, as shown in FIG. 35, the all-sector main power area map creation unit 607 creates a circular area Ci defined by the radius ri of the main power area centered on the antenna position of the base station, and the sector radio wave reach range Ri. (A hatched portion in FIG. 35) is obtained as a main power area map MRi in which the sector is the main power. Thus, by using the sector radio wave reachable range, the main power area map can be created with higher accuracy.

  Furthermore, in the above-described embodiment, the delay position information acquisition unit 602 acquires the delay information and the PRACH-PD position information output from the position collection unit 502. However, the present invention is not limited to this. For example, 3GPP (Third You may acquire the positioning position result of the moving apparatus 100 obtained by Timing Advance Type2 positioning defined by the standard of Generation Partnership Project. This Timing Advance Type 2 positioning is a technique for positioning the position of the mobile device 100 based on the propagation time (delay amount) of radio waves, as in the PRACH-PD positioning method. Based on the location information of the mobile device 100 obtained by this Timing Advance Type 2 positioning, in the functional units in the main power area creation parameter estimation unit 610, as in the case of using the delay information or the PRACH-PD location information described above. The radius of the main power area can also be obtained.

  Furthermore, in the above-described embodiment, the GPS position information acquisition unit 601 acquires the GPS position information output from the position collection unit 502. However, the present invention is not limited to this. For example, the standard of LTE (Long Term Evolution) You may acquire the positioning position result of the moving apparatus 100 obtained by the OTDOA (Observed time difference arrival) positioning defined by the standard. This OTDOA positioning is a technique for estimating the position of the mobile device 100 by a method such as so-called three-point surveying based on the propagation time (delay amount) of radio waves to the antennas of a plurality of adjacent base stations. Based on the position information of the mobile device 100 obtained by this OTDOA positioning, the main power area estimation unit 613 can also determine the radius of the main power area, as in the case of using the GPS position information described above.

  601 ... GPS position information acquisition unit (second correspondence information acquisition unit), 602 ... Delay position information acquisition unit (first correspondence information acquisition unit), 610 ... Main power area creation parameter estimation unit (main power area estimation device) 611 ... Sector type determination unit (sector type determination unit), 612 ... Antenna position calculation unit (antenna position calculation unit), 613 ... Main power area estimation unit (main power area estimation unit), 615 ... Radiation direction calculation unit (radiation) Direction calculation means), 616... Radiation width calculation section (radiation width calculation means).

Claims (11)

Correspondence information acquisition means for acquiring correspondence information in which positioning position information of a mobile device obtained by predetermined positioning and a sector identifier of a sector where the mobile device is located are associated with each other;
Of the correspondence information acquired by the correspondence information acquisition means, the mobile station distribution obtained by the positioning position information associated with the predetermined sector identifier, and the sector corresponding to the predetermined sector identifier Main power area estimation means for estimating a radius of a main power area where the sector corresponding to the predetermined sector identifier is the main power based on the position of the antenna of the base station to be formed;
A main power area estimation device. The correspondence information acquisition means includes
The predetermined positioning is a first positioning based on a propagation time of a radio wave transmitted and received between the mobile device and the base station, and the positioning position information is obtained by the first positioning. First positioning information is obtained when the correspondence information is first correspondence information in which the first positioning position information and the sector identifier are associated with each other. The main power area estimation apparatus according to claim 1, which is correspondence information acquisition means. The correspondence information acquisition means includes
In the predetermined positioning, the positioning position information is second positioning position information obtained by the second positioning, and the correspondence information is associated with the second positioning position information and the sector identifier. 2. The main power area estimation apparatus according to claim 1, which is a second correspondence information acquisition unit that acquires the second correspondence information when the second correspondence information is provided. The predetermined positioning includes a first positioning based on a propagation time of radio waves transmitted and received between the mobile device and the base station, and a second positioning different from the first positioning,
The correspondence information acquisition means includes
The predetermined positioning is the first positioning, the positioning position information is first positioning position information obtained by the first positioning, and the correspondence information is the first positioning position information. And the first correspondence information acquisition means for acquiring the first correspondence information when the sector identifier is associated with the first correspondence information;
as well as,
The predetermined positioning is the second positioning, the positioning position information is second positioning position information obtained by the second positioning, and the correspondence information is the second positioning position information. Second correspondence information acquisition means for acquiring the second correspondence information when the second correspondence information is associated with the sector identifier,
The main power area estimation apparatus according to claim 1, comprising: The first positioning position information includes position information calculated based on propagation time information related to the propagation time,
The main power area estimation means obtains a distance between the antenna and the mobile device, and obtains a distribution of the mobile device based on the obtained distance between the antenna and the mobile device.
The main influence area estimation apparatus according to claim 2 or 4.   The main position according to any one of claims 2, 4 and 5, further comprising antenna position calculation means for calculating the position of the antenna based on the distribution of the mobile devices obtained based on the first correspondence information. Power area estimation device.   Of the first correspondence information, based on the distribution of the mobile station obtained by the first positioning position information associated with the predetermined sector identifier, the sector type corresponding to the predetermined sector identifier is determined. The main power area estimation apparatus according to any one of claims 2 and 4 to 6, further comprising sector type determination means for determining.   Based on the distribution of the mobile device obtained from the first positioning position information associated with the predetermined sector identifier among the first correspondence information, the radiation direction of the radio wave radiated from the antenna is calculated. The main influence area estimation apparatus as described in any one of Claims 2 and 4-7 further equipped with the radiation direction calculation means to do.   Of the second correspondence information, a sector corresponding to the predetermined sector identifier is formed based on the distribution of the mobile devices obtained by the second positioning position information associated with the predetermined sector identifier. 5. The main power area estimation apparatus according to claim 3, further comprising radiation width calculating means for calculating a radiation width of a radio wave radiated from an antenna of the base station. The first positioning position information or the second positioning position information includes coordinate information indicating the position of the mobile device,
The main power area estimation means obtains the distribution of the mobile devices based on the coordinate information.
The main influence area estimation apparatus as described in any one of Claims 2-9. Correspondence information acquisition step of acquiring correspondence information in which positioning position information of a mobile device obtained by predetermined positioning and a sector identifier of a sector where the mobile device is located are associated with each other;
Among the correspondence information acquired in the correspondence information acquisition step, the distribution of the mobile device obtained by the positioning position information associated with the predetermined sector identifier, and the sector corresponding to the predetermined sector identifier A main power area estimation step for estimating a radius of a main power area in which a sector corresponding to the predetermined sector identifier is a main power based on the position of the antenna of the base station to be formed;
Main power area estimation method including

Images