JP5797120B2 - Information analysis apparatus and population distribution calculation method - Google Patents

Description

  The present invention relates to an information analysis apparatus and a population distribution method for calculating a population distribution based on information of a mobile communication network.

  Patent Document 1 describes that the number of terminals per sector is estimated by apportioning the number of locations based on the location registration processing of mobile devices managed by the mobile phone network by the number of incoming calls for each sector. Yes.

International Publication No. 2010/116905

  However, the technique described in Patent Document 1 cannot calculate a population distribution for each terminal attribute. That is, the technique described in Patent Document 1 calculates population distribution information based on the number of incoming calls and the number of visited areas in a certain sector, but the number of incoming calls cannot be calculated for each attribute. . In general, the attribute is information such as the sex and age of the subscriber, but the number of incoming calls is information obtained by measuring the incoming signal to the mobile device or the outgoing signal from the mobile device. Therefore, it is difficult to measure for each attribute. Therefore, with the technique described in Patent Document 1, it is extremely difficult to calculate an accurate population distribution for each attribute in the target sector.

  Therefore, an object of the present invention is to provide an information analysis apparatus and a population distribution calculation method that can accurately calculate a population distribution for each attribute.

In order to solve the above-described problem, the information analysis apparatus of the present invention includes, for each sector, a signal number storage unit that stores the number of signals communicated with the terminal, and a specified attribute in the specified sector. Terminal number acquisition means for acquiring the number of terminals and the number of terminals of all attributes in the sector, incoming number acquisition means for acquiring the number of incoming calls in the designated sector from the signal number storage means, and the designated sector A number-of-zones acquisition means for acquiring the number of areas in a cluster composed of a plurality of sectors including the number of calls received by the number-of-calls acquisition means, and the total number of calls received by summing up the number of calls received the number and the located number in the clusters, obtained by the station number acquisition unit, all in the terminal number and the sector of the attributes specified in the designated sector Based on the attribute number of the terminal and, a, and population distribution calculating means for calculating the population distribution in the sector.

The population distribution calculation method of the present invention includes a terminal number acquisition step of acquiring the number of terminals having a specified attribute in a specified sector and the number of terminals having all the attributes in the sector, An incoming number acquisition step of acquiring the number of incoming calls in the designated sector from a signal number storage means for storing the number of signals communicated in the network; and a presence in a cluster comprising a plurality of sectors including the designated sector. The number- of-areas acquisition step for acquiring the number of areas; the number of incoming calls acquired by the number-of-incoming-acquisition step; the total number of incoming calls obtained by counting the number of incoming calls in units of clusters; the number of areas in the cluster; and the number of terminals acquired by the acquisition step, on the number of terminals all attributes in the terminal number and the sector of the attributes specified in the designated sector And Zui, a, and population distribution calculating step of calculating a population distribution in the sector.

According to the present invention, the number of terminals having the designated attribute in the designated sector and the number of terminals having all the attributes in the sector are obtained, the number of incoming calls in the designated sector is obtained, and the designated sector is obtained. The number of locations in a cluster composed of a plurality of sectors including is acquired. Then, the population distribution in the sector is calculated based on the acquired number of incoming calls, the total number of incoming calls obtained for each cluster, the number of in-zones in the cluster, and the acquired number of terminals. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area.

  In the information analysis apparatus according to the present invention, the number-of-areas acquisition unit may acquire the number of areas totaled for each location registration area including the designated sector.

  According to the present invention, it is possible to calculate the population distribution for each sector and attribute by calculating the population distribution based on the number of areas that are counted in units of location registration areas including the designated sector.

  In the information analysis apparatus according to the present invention, the number-of-zones acquisition unit may acquire the number of zones counted in a cluster unit composed of sectors adjacent to the designated sector.

  According to the present invention, by calculating the population distribution based on the number of located areas that are aggregated in units of clusters composed of sectors adjacent to the designated sector, it is possible to more accurately calculate the number of terminals in units of sectors. Population distribution can be calculated for each attribute.

  Further, in the information analysis apparatus of the present invention, the number-of-zones acquisition means is a cluster unit composed of a plurality of sectors including the designated sector, and is counted in a cluster unit across a plurality of location registration areas. The number of locations may be acquired.

  According to the present invention, it is possible to obtain an excessive number in each sector by obtaining the number of located areas that are cluster units including a plurality of sectors including a designated sector and are aggregated in cluster units across a plurality of location registration areas. It is possible to absorb an error when performing aggregation or under-aggregation.

  In the information analysis apparatus of the present invention, the number-of-zones acquisition unit forms a plurality of cluster units including sectors having different combinations with the designated sector, and calculates the number of zones in the cluster unit. The population distribution calculating means may calculate the population distribution based on the average value of the number of located areas.

  According to the present invention, a plurality of cluster units including sectors having different combinations with the designated sector are formed, the number of areas in the cluster unit is calculated, and the population is calculated based on the average value of the number of terminals. By calculating the distribution, it is possible to calculate the population distribution for each sector and attribute with higher accuracy.

  In the information analysis apparatus of the present invention, the number-of-zones acquisition unit forms a plurality of cluster units including sectors having different combinations with the designated sector, and calculates the number of zones in the cluster unit. The population distribution calculating means may calculate the population distribution in each cluster unit using the number of each area, and calculate the average value of the population distribution.

  According to the present invention, a plurality of cluster units including sectors having different combinations with the designated sector are formed, and the number of spheres in the cluster unit is calculated, and each cluster number is used to calculate each cluster. By calculating the population distribution in units and calculating the average value of the population distribution, the population distribution for each sector and attribute can be calculated with higher accuracy.

  Further, in the information analysis apparatus of the present invention, the population distribution calculating means configures a location registration area including the number of in-zones counted in cluster units composed of sectors around the designated sector and the designated sector. The population distribution may be calculated on the basis of the number of located areas aggregated by the set unit of the sector.

  According to the present invention, based on the number of visited areas counted in units of clusters composed of sectors around the designated sector and the number of visited areas counted in the set unit of sectors constituting the location registration area including the designated sector. By calculating the population distribution, the population distribution for each sector and attribute can be calculated with higher accuracy.

  In the information analysis apparatus of the present invention, the number-of-zones acquisition unit acquires the number of terminals in the cluster by acquiring the number of terminals in each sector that constitutes the cluster and totalizing the number of terminals. Also good.

  Thereby, since the number of terminals acquired for each sector is tabulated, the number of located areas with higher accuracy can be acquired.

Further, in the information analysis apparatus of the present invention, the terminal number acquisition unit is configured to position by crossing a location registration area from the number of terminals having the specified attribute in the specified sector and the set of terminals having the specified attribute. The number of terminals after removal may be calculated by removing a set of registered terminals, and a numerical value that is prorated by a predetermined ratio of the number of terminals and the number of terminals after removal may be used as the number of terminals.

  According to the present invention, the number of terminals having the specified attribute in the specified sector and the number of terminals after removal by removing the number of terminals registered by straddling the position registration area are calculated, and the number of terminals and the terminals after removal are calculated. By making the number apportioned at a predetermined ratio of the number the number of terminals, it is possible to calculate the population distribution for each sector and attribute with higher accuracy.

Further, in the information analysis layer of the present invention, the population distribution calculation means uses the number of terminals L _ij of the attribute j mobile station in the designated sector i, the number of serving areas in the cluster as N, and the number of incoming calls in the cluster. , C _i may be used to calculate the population distribution using the number of terminals L _ij , the number N of visited areas, and the number of incoming calls C _i .

According to the present invention, the number of terminals L _ij of the mobile station having the attribute j in the sector i is set to N, the number of serving areas in the cluster is set to N, and the number of incoming calls in the cluster is set to C _i. To do. Thereby, the population distribution for each sector and attribute can be calculated with higher accuracy.

  According to the present invention, it is possible to accurately calculate the population distribution for each sector and attribute, which eliminates the influence of straddling the location registration area.

It is a figure which shows the system configuration | structure of the communication system 10 of this embodiment. 1 is a block diagram showing a functional configuration of a communication system 10. FIG. 2 is a block diagram showing a hardware configuration of an information analysis apparatus 600. FIG. It is a figure which shows the specific example of the positional infomation table which memorize | stores the positional infomation on the moving apparatus. It is a conceptual diagram for demonstrating terminal number estimation. It is a conceptual diagram for demonstrating the calculation method of terminal number estimation. It is a flowchart which shows the process of terminal number estimation. It is a flowchart which shows the feature-value calculation process. It is a flowchart which shows the population distribution calculation method of the information analyzer 600 in this embodiment. It is a conceptual diagram which shows the relationship between the sector used as a measurement object, and a cluster. It is a block diagram which shows the function structure of the information analyzer 600a in 2nd embodiment. It is a flowchart which shows the process of the information analyzer 600a in 2nd embodiment. It is a conceptual diagram which shows the relationship between the sector used as a measurement object, and a cluster. It is a flowchart which shows the process of the information analyzer 600a in a modification. It is a block diagram which shows the function structure of the information analyzer 600b in 3rd embodiment. It is a flowchart which shows the process of the information analyzer 600b in 3rd embodiment. It is a block diagram which shows the function structure of the information analyzer 600c in 4th embodiment. It is a flowchart which shows the process of the information analyzer 600c in 4th embodiment. It is the schematic diagram which showed the relationship between the transmission timing of the position registration signal transmitted by the mobile device 100, and the estimated stay period. It is explanatory drawing which shows the specific example in the measuring method of the number of entering.

  Embodiments of the present invention will be described with reference to the accompanying drawings. Where possible, the same parts are denoted by the same reference numerals, and redundant description is omitted.

[First embodiment]
[Configuration of communication system]
FIG. 1 is a diagram illustrating a system configuration of a communication system 10 according to the present embodiment. As shown in FIG. 1, the communication system 10 includes a mobile device 100, a BTS (base station) 200, an RNC (radio network control device) 300, an exchange 400, a management center 500, and an HLR (Home Location Register) 700. It consists of The management center 500 includes a social sensor unit 501, a petamining unit 502, a mobile demography unit 503, and a visualization solution unit 504.

  The RNC 300 receives the RRC Connection Request signal transmitted by the mobile device 100 via the BTS 200. At this time, it is possible to count the number of signals in which voice incoming, voice outgoing, packet incoming, and packet outgoing are set as parameters of the RRC Connection Request signal. According to the standard specification “Radio Resource Control (RRC) Protocol Specification: 3GPP TS 25.331”, the RNC 300 stipulates the number of signals per sector, and this method is based on this. The method for counting the number of signals is not limited to 3GPP, and may be based on LTE (Long Term Evolution).

  The exchange 400 collects the location registration signals transmitted by the mobile device 100 via the BTS 200 and the RNC 300. The exchange 400 manages the mobile devices 100 and the like in units of location registration areas, and collects the location registration signals transmitted by the mobile device 100 to obtain the number of mobile devices 100 registered in the location registration area (the number of registrations). ) It is grasped and stored for each location registration area. The exchange 400 outputs the number of mobile devices 100 registered in the stored location registration area (number of registrations) to the management center 500 at a predetermined timing or in response to a request from the management center 500. Here, in general, the RNC 300 is composed of about 1,000 pieces and is arranged throughout Japan. On the other hand, about 300 exchanges 400 are arranged in Japan.

  As described above, the management center 500 includes the social sensor unit 501, the petamining unit 502, the mobile demography unit 503, and the visualization solution unit 504. In each unit, the location registration transmitted by the mobile device 100 is registered. Performs statistical processing using information on signals and outgoing / incoming calls.

  The social sensor unit 501 is a server device that collects data including the number of mobile devices 100 registered in the location registration area stored in the exchange 400 from each exchange 400. The social sensor unit 501 is configured to receive data periodically output from the exchange 400 and acquire data from the exchange 400 according to a predetermined timing in the social sensor unit 501.

  The petamining unit 502 is a server device that converts data received from the social sensor unit 501 into a predetermined data format. For example, the petamining unit 502 performs the sorting process using the user ID as a key or for each area.

  The mobile demography unit 503 is a server device that performs aggregation processing on the data processed in the petamining unit 502, that is, count processing for each item. For example, the mobile demography unit 503 can count the number of users located in a certain area, or count the distribution of the located areas.

  The visualization solution unit 504 is a server device that processes the data aggregated in the mobile demography unit 503 so as to be visible. For example, the visualization solution unit 504 can map the aggregated data on a map. Data processed by the visualization solution unit 504 is provided to companies, government offices or individuals, and is used for store development, road traffic surveys, disaster countermeasures, environmental countermeasures, and the like. It should be noted that the information statistically processed in this way is processed so that individuals are not specified so as not to infringe privacy.

  The social sensor unit 501, petamining unit 502, mobile demography unit 503, and visualization solution unit 504 are all configured by the server device as described above, and although not shown, the basic configuration of a normal information processing device Needless to say, it includes a CPU, a RAM, a ROM, an input device such as a keyboard and a mouse, a communication device that communicates with the outside, a storage device that stores information, and an output device such as a display and a printer.

  FIG. 2 shows a functional configuration of the communication system 10. As shown in FIG. 2, the communication system 10 includes a plurality of mobile devices 100 located in sectors controlled by each of the plurality of BTSs 200, an RNC 300 that controls the BTSs 200, an exchange 400, and an information analysis device 600. Composed. The information analysis apparatus 600 corresponds to the mobile demography unit 503 and the visualization solution unit 504 shown in FIG. Regarding the functions corresponding to the social sensor unit 501 and the petamining unit 502 in FIG. 1, these notations are omitted in FIG.

  The RNC 300 includes an RNC communication control unit 302, a location registration signal receiving unit 303, and a signal number measuring unit 304. The exchange 400 includes an exchange communication control unit 401, a conversion unit 402, a location registration signal processing unit 404, and a storage unit 403. In addition, the information analysis apparatus 600 includes an information analysis apparatus communication control unit 601, an incoming call number acquisition unit 602 (incoming call number acquisition unit), a visited area acquisition unit 603 (present area number acquisition unit), and a terminal number acquisition unit 604 (number of terminals). An acquisition unit), a position information storage unit 605, an incoming call number storage unit 606 (signal number storage unit), an output unit 610, and a population distribution calculation unit 611 (population distribution calculation unit).

  First, the RNC 300 will be described. The RNC communication control unit 302 is a part that performs communication connection with the mobile device 100 via the BTS 200. For example, the RNC communication control unit 302 performs communication connection processing based on transmission processing from the mobile device 100 and communication connection processing based on a location registration request. . In the present embodiment, the RNC communication control unit 302 can transmit an Initial UE Message used for communication connection processing to the exchange 400. This Initial UE Message includes instruction information (location registration signal) indicating a transmission or location registration request, an ID for uniquely identifying the mobile device 100, and location information. In addition, Initial UE Message can also add information on the number of outgoing and incoming signals. Here, the ID may be, for example, ID information as a temporary ID issued by the exchange 400 when the mobile device 100 is connected to the network.

  Note that the RNC 300 can calculate the GAI (Geographical Area ID) as to which position in the sector the mobile device 100 is located based on the signal delay obtained when the RRC connection request is further processed. it can. The position of the mobile device 100 can also be specified based on the sector identifier and the position in the sector.

  The location registration signal receiving unit 303 is a part that receives the location registration signal transmitted by the mobile device 100 via the RNC communication control unit 302.

  The signal number measuring unit 304 is a part that measures the number of sector unit signals (number of signals), which is information indicating the amount of signals transmitted and received by the mobile device 100 located in the sector for each sector.

  Here, as the amount of signal, for example, the number of incoming voice calls, the number of incoming packet calls, the number of outgoing voice calls, the number of outgoing packet calls, and the like determined by a method using the parameters of the RRC connection request signal described above can be used. However, there is no intention to limit to this.

  The amount of signal is preferably measured by a method of summing at least one of the number of incoming voice calls and the number of incoming packets. When voice transmission and packet transmission are used, it can be considered that these (number of transmissions) are related to the action of a human being as a user. That is, when a user makes a call, it is considered that there are many places such as offices, houses, stations, and the like.

  On the other hand, when it is limited to an incoming call, it is not directly related to the action of a human being who is a user of the called mobile device. Therefore, it is possible to obtain population distribution information that more appropriately reflects the actual population distribution by measuring the amount of signal by a method of summing the number of incoming voice calls and the number of incoming packet packets. However, it goes without saying that even if the number of signals per sector is measured from the number of transmissions, it can be sufficiently put into practical use.

  The signal number measurement unit 304 transmits the measured number of sector unit signals to the information analysis apparatus 600 via the RNC communication control unit 302.

  Next, the exchange 400 will be described. The exchange communication control unit 401 is a part that receives an Initial UE Message transmitted from the RNC 300 and performs a communication connection process using the Initial UE Message.

  The conversion unit 402 is a part that converts an ID such as a temporary ID included in the Initial UE Message received by the exchange communication control unit 401 into a telephone number. In the conversion process, the conversion unit 402 extracts a telephone number associated with an ID such as a temporary ID from the subscriber profile information storage unit that stores the subscriber profile information, and converts the telephone number into the extracted telephone number. The subscriber profile information storage unit is provided in, for example, an HLR (Home Location Register) 700, and here manages and stores an ID such as a temporary ID and a telephone number in association with each other.

  In the present embodiment, the mobile device 100 transmits a location registration signal when the mobile device 100 moves across the location registration area or at predetermined intervals. As a result, the exchange communication control unit 401 can register its location in the HLR 700. The standard specification “Mobile Application Part (MAP) specification: 3GPP TS 29.002” describes a method in which the exchange 400 manages location registration. The processing in the exchange 400 according to this embodiment is based on this method.

  Next, the information analysis apparatus 600 will be described. The information analysis device communication control unit 601 is a part that controls communication between the information analysis device 600 and the exchange 400. The information analysis apparatus 600 is an apparatus of the management center 500, and is realized by a hardware configuration as shown in FIG.

  FIG. 3 is a hardware configuration diagram of the information analysis apparatus 600. As shown in FIG. 3, the information analysis apparatus 600 shown in FIG. 2 physically includes a CPU 11, a RAM 12 and a ROM 13 that are main storage devices, an input device 14 such as a keyboard and a mouse that are input devices, a display, and the like. The computer system includes an output device 15, a communication module 16 that is a data transmission / reception device such as a network card, an auxiliary storage device 17 such as a hard disk or a semiconductor memory, and the like. Each function described in FIG. 2 has the input device 14, the output device 15, and the communication module 16 under the control of the CPU 11 by reading predetermined computer software on the hardware such as the CPU 11 and the RAM 12 shown in FIG. 3. This is realized by reading and writing data in the RAM 12 and the auxiliary storage device 17. Hereinafter, each functional block will be described based on the functional blocks shown in FIG.

  The incoming call number acquisition unit 602 is a part that acquires the number of signals in units of sectors measured by the signal number measurement unit 304 and stored in the incoming call number storage unit 606.

  The in-zone number acquisition unit 603 is a part for acquiring the in-zone number of the mobile device in the location registration area including the designated sector based on the in-zone information stored in the HLR 700, the exchange 400, or the like. That is, the number-of-zones acquisition unit 603 identifies the location registration area to which the sector belongs, and acquires the number of locations of the mobile device 100 located in the location registration area.

  The terminal number acquisition unit 604 is a part that acquires the number of terminals in the sector and attribute designated by the operator or the like based on the subscriber information calculated by the location registration signal processing unit 404 and stored in the HLR 700. The number of terminals, for example, the feature value described later is a numerical value calculated for each position of each mobile device, and at least the identifier of the mobile device (may be non-identified), the feature value is included in the management table. , The sector identifier in the area and the position information are stored in association with each other. Then, the terminal number acquisition unit 604 can add up the number of terminals for each attribute j by specifying the sector i and the attribute j using this management table. Note that the attribute j can be determined from the identifier of the mobile device.

  The position information storage unit 605 is a part that stores a position information table describing position information at each positioning time of each mobile device 100. FIG. 4 is a diagram showing a specific example of a position information table that stores the position information. As shown in FIG. 4, in the position information table, a user identifier and a transition state of the position information are stored, and time information and position information are described in association with each other. The location information table stored here may store location information measured using GPS in each mobile device 100, or store location information based on the location information stored in the HLR 700. Also good. The user identifier is an identifier for distinguishing users from each other, and is converted so that the user cannot be specified. This user identifier has been subjected to a so-called de-identification process. For example, the non-discriminating process can be realized by using an irreversible one-way function such as a hash function and converting it into another completely different identifier.

  The incoming call number storage unit 606 is a part that stores the number of signals measured by the signal number measurement unit 304. Here, the number of outgoing / incoming calls is stored for each sector.

  Based on the number of terminals acquired by the number-of-terminals acquisition unit 604, the number of areas acquired by the number-of-areas acquisition unit 603, and the number of incoming calls acquired by the number of incoming calls acquisition unit 602, the population distribution calculation unit 611 This is the part that calculates population distribution information in units. This population distribution information can be calculated for each attribute of the mobile device 100 (user).

Specifically, when the sector i and the attribute j are designated, the population distribution calculation unit 611 is the number of signals per sector, the number of incoming calls C _i acquired by the number of incoming calls acquisition unit 602, and the number of visited areas acquisition unit 603. The corrected number of terminals S _{ij in} which the error due to each mobile device 100 straddling the location registration area is suppressed based on the number of visited areas Na acquired by the above and the number of terminals L _ij acquired by the terminal number acquiring unit 604. The calculation is performed according to the following formula.

The equation (1), relative to the located number Na, the ratio of the number of terminals L _ij in the terminal number sigma (i) L _ij and attribute j sector i in all attributes, clusters that sector i belongs, namely by the present embodiment multiplies the ratio of the incoming number C _i of the sector i in the total incoming number sigma (i) C _i in the location registration area, the more accurate number of terminals S _ij attribute j in sector i It can be obtained. That is, since the number of incoming calls cannot be obtained for each attribute in the number of incoming calls C _i , the number of terminals L _ij in the attribute j of the sector i is calculated by calculating the ratio of the attribute j in all the attributes of the sector i. The number of terminals using the number can be calculated.

In the following, Σ (i) represents a numerical value obtained by accumulatively adding in the range designated for i for convenience. For example, when indicating the total number of incoming calls in the location registration area LAa, Σ (i ) C _i indicates C ₁ + C ₂ +... C _n (n is the number of the sector constituting the location registration area LAa). Σ (j) L _ij is L _i1 + L _i2 +... L _ij and indicates the number of terminals of all attributes in sector i. The same applies to other mathematical expressions.

ここで、Ｋａは、位置登録エリアごとに定まる比例定数である。 It can be understood that the equation (1) is an equation that can calculate a more correct numerical value as described below. First, it is assumed that the ratio between the number of incoming calls in each sector and the number of terminals in the same location registration area is constant. In that case, the number of incoming calls can be expressed as follows.

Here, Ka is a proportionality constant determined for each location registration area.

ここで、

である。 At this time, since the sum of the number of incoming sectors in the location registration area and the number of locations in the location registration area are also in the same proportional relationship, equation (2) can be modified as follows.

here,

It is.

Here, if the ratio for each attribute is correct in the number of terminals based on the location registration signal, the following equation holds.

In this way, the population distribution of the attribute j in the sector i can be calculated based on the number of incoming calls C _i , the number of visited areas Na, and the number of terminals L _ij .

  Here, in order to calculate a more accurate population distribution, it is also possible to multiply the feature amount by a predetermined expansion coefficient. In other words, statistical information (population by census etc.) of a region based on a certain sector or a wider region including that region is stored in advance, and the ratio between the number of terminals and the actual population in that wide region The population distribution closer to the real number can be calculated by obtaining an expansion coefficient as follows and multiplying this ratio by the number of terminals in the sector.

  More specifically, the enlargement factor can be calculated as follows. As an example of the expansion coefficient, the reciprocal of “the product of the location ratio and the terminal penetration rate (that is, the ratio of the number of locations to the population)” can be used. Here, the “area ratio” means the ratio of the area number to the contracted number, and the “popularity ratio” means the ratio of the contracted number to the population. Such an enlargement factor is desirably derived for each population estimation unit described above, but is not essential.

  Further, the expansion coefficient may be derived using, for example, the number of terminals (the number of existing areas) estimated based on the feature amount and the observation period length as follows. That is, a feature amount is obtained from position data, and the number of terminals for each enlargement coefficient calculation unit is obtained based on the feature amount and the observation period length to obtain user number pyramid data, and statistical data (for example, Basic Resident Register) As above, population pyramid data in the same expansion coefficient calculation unit obtained in advance is acquired. And the acquisition rate (namely, the number of area / population) of the position data for every expansion coefficient calculation unit is calculated in the user number pyramid data and the population pyramid data. The “location data acquisition rate (that is, the number of locations / population)” obtained here corresponds to the “product of the location rate and the terminal penetration rate” described above. The reciprocal of the “position data acquisition rate” obtained in this way can be derived as an expansion coefficient. In addition, as an enlargement factor calculation unit for calculating an enlargement factor, for example, every address prefecture, every 5 or 10 years of age, every gender, every hour, etc. may be adopted. A combination of two or more of these may be employed. For example, if the enlargement coefficient calculation unit is “male in the 20s in Tokyo”, the males in the 20s in Japan who live in Tokyo (that is, the address information in the user attribute is Tokyo). The corresponding position data is extracted and the number of terminals is totaled to obtain user number pyramid data, and population pyramid data relating to a 20-year-old man living in Tokyo is obtained from statistical data. When obtaining the above-mentioned user pyramid data, for the condition of “resident in Tokyo”, the address information in the user attribute is Tokyo instead of extracting only the location data of users residing in Tokyo. Extract location data. Then, the location data acquisition rate (that is, the number of people in the area / population) of the expansion coefficient calculation unit (here, a man in the 20s in Tokyo) is calculated from the user number pyramid data and the population pyramid data. The reciprocal of the “position data acquisition rate” can be derived as an expansion factor. In the present application, the enlargement coefficient calculation unit and the population estimation unit are described as being equal. However, this is merely an example, and the present invention is not limited to this.

In the present embodiment, the number of terminals L _ij is used, and the number of location registration signals may be simply counted, or expressed based on a feature amount indicating the amount of terminals as described below. It is good as what was done.

  Here, a feature amount calculation method by the terminal number acquisition unit 604 will be described. The terminal number acquisition unit 604 calculates a feature amount for each first position data (observation target position data) that is each position information of each mobile device 100. For example, the terminal number acquisition unit 604 calculates the difference between the position acquisition time of the second position data and the position acquisition time of the third position data as a feature amount for the first position data. In addition, when the position acquisition time of the second position data is an abnormal value, the number-of-terminals acquisition unit 604 uses the position acquisition time of the first position data and the position acquisition time of the second position data as an example here. When the difference is larger than a predetermined reference value (for example, 1 hour), the time of the second position data is set to a time that is retroactive from the position acquisition time of the first position data by a predetermined time (for example, 1 hour). The feature amount for the first position data is calculated as the acquisition time. Similarly, when the position acquisition time of the third position data is an abnormal value, the terminal number acquisition unit 604, as an example here, the position acquisition time of the first position data and the position acquisition time of the third position data Is greater than a predetermined reference value (for example, 1 hour), the time advanced from the position acquisition time of the first position data to the future by a predetermined time (for example, 1 hour) is set in the third position data. The feature amount for the first position data is calculated using the position acquisition time. The process when the position acquisition times of the second and third position data are abnormal values is not an indispensable process, but the mobile device 100 is located outside the service area by performing the above process. When the position data acquisition time interval becomes abnormally long due to the power of the mobile device 100 being turned off, etc., the influence of the abnormally long acquisition time interval is prevented from becoming excessive. be able to.

  The number-of-terminals acquisition unit 604 estimates the number of terminals located in the observation area during the observation period based on the feature amount of the observation target position data and the observation period length that is the difference between the observation start time and the observation end time. . Although details will be described later, the terminal number acquisition unit 604 estimates a numerical value obtained by dividing the sum of the feature amounts of the observation target position data by twice the observation period length as the number of terminals.

即ち、各端末ａ_ｉの観測期間内のセクタＳの滞在時間ｔ_ｉの総和を観測期間の長さＴで除した結果を、端末数ｍとして推計する。ただし、端末ａ_ｉの観測期間内のセクタＳの滞在時間ｔ_ｉの真の値は観測不能であるが、各端末ａ_ｉは信号（例えば位置登録信号）を発信し、それらの信号は観測可能である。 [Concept of terminal number estimation and calculation method]
Next, the concept and calculation method of terminal number estimation will be described. As model shown in FIG. 5, during a certain observation period (length T), n pieces of terminal _{a 1, a 2, ...,} a n passes through the sector S, the observation period for each terminal a _i It is assumed that the stay time of the sector S is t _i (0 <t _i ≦ T). At this time, the number m of terminals present in the sector S (actually the average value of the number m of terminals present in the sector S within the observation period) is expressed by the following equation (6).

That is, the result of dividing the sum of the stay times t _i of the sectors S within the observation period of each terminal a _i by the length T of the observation period is estimated as the number of terminals m. However, although the true value of the stay time t _i of the sector S within the observation period of the terminal a _i is not observable, each terminal a _i transmits a signal (for example, a location registration signal), and these signals can be observed. It is.

（ｘ_ｉは、端末ａ_ｉが観測期間内にセクタＳで発信した信号の総数）とすると、端末数の推計とは、観測された信号ｑ_ｉｊ（ｊは１以上ｘ_ｉ以下の整数）からｍの値を推計することに他ならない。 The signals transmitted by the terminal a _i from the sector S during the observation period

If x _i is the total number of signals transmitted by the terminal a _i in the sector S during the observation period, the estimation of the number of terminals is based on the observed signal q _ij (j is an integer from 1 to x _i ). It is none other than estimating the value of m.

Now, a method for calculating the number of terminals will be described with reference to FIG. Let p _i be the density at which the signal q _ij is transmitted from the terminal a _i (ie, the number of signals per unit time). At this time, if the probability that the signal is transmitted is independent of the sector, the expected value E (x _i ) of the total number x _i of signals transmitted from the terminal a _i in the sector S within the observation period is E (x _i) since a = _{t i} × _{p i,} the expected value E _{(t i)} the following expression for the residence time _{t i} of the sector S in the observation period of the terminal _{a i} (7) is satisfied.
E (t _i ) = x _i / p _i (7)
Here, when the transmission time of the signal _{q ij} was _{u ij,} density _{p ij} of the signal _{q ij} is given by the following equation (8).
p _ij = 2 / (u _{i (j + 1} ) −u _{i (j−1)} ) (8)
Here, when the signal q _ij is a signal related to the first position data, the signal q _{i (j−1)} is a signal related to the second position data, and the signal q _{i (j + 1)} is related to the third position data. Corresponds to the signal. In the present embodiment, the transmission time _{u i} of the signal of the second position data _{q i (j-1)} transmission time _{u i} of _(j-1) and the signal _{q i} of the third position data _{(j + 1) (} The difference of _{j + 1)} , that is, (u _{i (j + 1)} −u _{i (j−1)} ) in the above equation (8) is set as the feature quantity w _ij for the first position data. Therefore, the above formula (8) is as follows. That is, the feature quantity w _ij can be calculated in association with the reciprocal of the density p _ij .
p _ij = 2 / (u _{i (j + 1)} −u _{i (j−1)} ) = 2 / w _ij (9)

で与えられるため、端末数ｍの推計値Ｅ（ｍ）は以下の式（１１）で計算することができる。

At this time, the density p _i is

Therefore, the estimated value E (m) of the number of terminals m can be calculated by the following equation (11).

As shown in the example of FIG. 6, the terminal a _i transmits the signals q _i1 , q _i2 , and q _i3 within the observation period and the period in which the terminal a _i stays in the sector S, and the signal q _i1 , The signal q _i0 is transmitted immediately before the signal q _i3 , and the signal q _i4 is transmitted immediately after the signal q _i3 , and the transmission times of the signals q _i0 , q _i1 , q _i2 , q _i3 , q _i4 are set to u _i0 , u _i1 , u, respectively. _{Assuming i2} , u _i3 , u _i4 , the above idea is that the stay time t _i of the sector S in the observation period of the terminal a _i is determined from (the midpoint of u _i0 and u _i1 ) (from u _i3 and u _i4 This is equivalent to estimating the period until the midpoint. The terminal a _i transmits the signal q _i4 while staying in the sector S, although it is not within the observation period. However, in order to maintain the unbiasedness of the estimation of the stay time t _i , here, as an example, a process is described in which the end time of the stay time t _i is not estimated to be the same as the end time of the observation period T. .

[Terminal number estimation process]
Hereinafter, the terminal number estimation process according to the terminal number estimation method of the present invention will be described. Here, as an example, it is assumed that the position information included in the position data of the mobile device 100 is given the sector number of the sector where the mobile device 100 is located.

  As shown in FIG. 7, first, the terminal number acquisition unit 604 acquires and stores location data from the HLR 700, the exchange 400, or the like (step S1 in FIG. 7). Thereby, the position information storage unit 605 stores position data for a plurality of times for a large number of users (mobile devices). In addition, after the process of step S1, you may perform the process after step S2 after time. That is, step S1 may be executed as a preliminary preparation for the processing after step S2. Here, the configuration is such that the position data is acquired from the HLR 700 or the exchange 400. However, the present invention is not limited to this, and the position information measured using the GPS or the like in the mobile device 100 is stored in the position information storage unit 605. May be remembered.

  Next, the number-of-terminals acquisition unit 604 acquires observation period information including a set of observation start time and observation end time, and the number-of-terminals acquisition unit 604 is observation area information associated with one or more pieces of position information. Is acquired (step S2). Here, it is assumed that a set of the observation start time T1 and the observation end time T2 is acquired as the observation period information, and the sector number S is acquired as the observation area information.

Next, the number-of-terminals acquisition unit 604 associates from the location information storage unit 605 with the sector number S that includes location acquisition time information that is after the observation start time T1 and before the observation end time T2 and that is observation area information. One or a plurality of position data including the position information (for example, the position information is the sector number S) is acquired as observation target position data (step S3). That is, the terminal number acquisition unit 604 acquires position data that satisfies the following conditions as observation target position data.
Condition 1: The position acquisition time is after the observation start time T1 and before the observation end time T2. That is, it is included in the observation period.
Condition 2: The position information is sector S.

  Next, the following steps S4 and S5 are executed for each of the acquired observation target position data. In step S4, the number-of-terminals acquisition unit 604 includes the same identification information as the first position data for the position data (first position data) for which the feature amount is obtained from the observation target position data. Of the data, the position acquisition time information of the position data immediately before the first position data (second position data) and the position data immediately after the first position data (third Position data) is acquired. It is not essential for the terminal number acquisition unit 604 to acquire the entire second and third position data, and it is only necessary to acquire the position acquisition time information included in the second and third position data.

  In step S5, the terminal number acquisition unit 604 calculates a feature amount for the first position data. The processing contents will be described with reference to FIG. Here, the position acquisition times of the first, second, and third position data are t1, t2, and t3, respectively. Further, a predetermined reference value (reference value relating to the difference between the position acquisition times of the first and second position data) serving as a reference for determining that the position acquisition time t2 of the second position data is an abnormal value is used as a reference. A value A (for example, 1 hour) is set, and a predetermined reference value (a position acquisition time of the first and third position data is determined as a reference for determining that the position acquisition time t3 of the third position data is an abnormal value) A reference value regarding the difference) is set as a reference value B (for example, 1 hour).

  The number-of-terminals acquisition unit 604 determines the difference between the position acquisition times of the first and second position data (that is, the difference between times t1 and t2) Da and the difference between the position acquisition times of the first and third position data (that is, The difference between the times t1 and t3) Db is calculated (step S11 in FIG. 8). Then, the terminal number acquisition unit 604 determines whether or not the difference Da between the position acquisition times of the first and second position data is larger than a predetermined reference value A (for example, 1 hour) (step S12), If the difference Da is larger than the reference value A, the time acquired from the position acquisition time t1 of the first position data in the past by a predetermined time (for example, one hour) is set as the position acquisition time t2 of the second position data. (Step S13). Next, the feature amount calculation unit 17 determines whether or not the difference Db between the position acquisition times of the first and third position data is larger than a predetermined reference value B (for example, 1 hour) (step S14). If the difference Db is larger than the reference value B, the time that has been advanced in the future by a predetermined time (for example, 1 hour) from the position acquisition time t1 of the first position data is set to the position acquisition time t3 of the third position data. (Step S15). Then, the terminal number acquisition unit 604 calculates a difference between the position acquisition time t2 of the second position data and the position acquisition time t3 of the third position data as a feature amount for the first position data (step S16). ). Thus, the processes in steps S4 and S5 for one piece of observation target position data (first position data) are completed.

  Thereafter, the processes in steps S4 and S5 described above are executed for each of the observation target position data, and when the execution is completed for all the observation target position data (positive determination in step S6), the process proceeds to step S7.

In step S7, the terminal number acquisition unit 604 obtains the numerical value obtained by dividing the sum of the feature values w _ij for the observation target position data by twice the observation period length T, as shown in the above-described equation (11). Is estimated as the number of terminals. As is clear from Equation (11), the terminal number acquisition unit 604 obtains the sum of (feature amount w _ij / 2) by dividing each feature amount w _ij for the observation target position data by 2, A numerical value obtained by dividing the obtained sum by the observation period length T may be estimated as the number of terminals. However, since the calculation method in which the sum of the feature quantities w _ij for the observation target position data is divided by twice the observation period length T as in the present embodiment, the number of divisions can be overwhelmingly smaller. There is an advantage that the load can be reduced.

  According to the above embodiment, when estimating the number of terminals using the position data, the correction using the acquisition time information of the previous and subsequent position data is performed to calibrate the influence of the change in the reception interval. The number of terminals can be estimated with high accuracy.

  In addition, by performing the processing when the position acquisition times of the second and third position data described above are abnormal values in the feature amount calculation processing, the mobile device 100 is located outside the service area. To prevent the influence of the abnormally long acquisition time interval from becoming excessive when the position data acquisition time interval becomes abnormally long due to the power of the mobile device 100 being turned off. Can do.

Next, a modified example related to the feature amount will be described. In the embodiment described above, the time difference (time difference between the second position data and the third position data) before and after the position data (first position data) for which the feature amount is to be obtained is calculated using the first An example of calculating the feature value of the position data has been shown. When this is expressed by an equation, the feature amount can be expressed by the following equation (12). In addition, the following formula | equation (12) is only the deformation | transformation of the formula (9) mentioned above, and is equivalent to a formula (9) (namely, it is not what changed the way of thinking of a formula (9)).
w _ij = u _{i (j + 1)} −u _{i (j−1)} (12)
This modification shows another variation of the feature amount calculation method calculated by the feature amount calculation unit 103.

In this modification, the feature quantity calculation unit 103 obtains the type information about the second position data and the third position data (for example, a position data generation factor described later) when obtaining the feature quantity of the first position data. (Generation timing)). Specifically, the feature amount calculation unit 103 calculates a correction coefficient α corresponding to the type information (generation factor here) of the third position data with respect to the time difference between the third position data and the first position data. A value obtained by multiplying the time difference between the first position data and the second position data by a correction coefficient β corresponding to the type information (generation factor in this case) of the second position data is calculated. calculate. However, in addition to the above, the feature amount calculation unit 103 may determine the correction coefficient α or β according to the type information of the first position data, and may also use the type information of the first and second position data. Accordingly, the correction coefficient β may be determined or the correction coefficient α may be determined according to the type information of the first and third position data. Then, the feature amount calculation unit 103 sets a value obtained by adding the values obtained by these multiplications as the feature amount of the first position data. The feature amount calculation process in the feature amount calculation unit 103 is expressed by the following equation (13).
w _ij = α (u _{i (j + 1)} −u _ij ) + β (u _ij −u _{i (j−1)} ) (13)

  As the type information about the second position data and the third position data, for example, when the position data is position registration information, information on the generation factor of the position registration information can be cited. It is included in the generated location registration information. Factors for generating location registration information include that the terminal has crossed a location registration area (Location Area) boundary, that it has been generated based on location registration that is performed periodically, that the attachment process is performed by turning on the terminal, etc. For example, the detachment process is executed when the power is turned off. The set values of the correction coefficients α and β are determined in advance corresponding to these generation factors. Then, the feature amount calculation unit 103 sets the correction coefficient α for the third position data according to the information about the generation factor of the third position data, and sets the correction coefficient α according to the information about the generation factor of the second position data. The correction coefficient β for the position data 2 may be set. Both the correction coefficients α and β may be set in advance to a value of 0 or more and 2 or less. However, this numerical range is not essential.

  For example, in the case of location registration information in which the location of the terminal and the generation opportunity of location registration information are irrelevant, such as location registration information based on location registration performed periodically, the expectation of the time spent in the current sector The value is considered to be the same before and after the location registration information is generated. On the other hand, in the case of location registration information generated by a terminal straddling a location registration area boundary, it can be determined that the terminal has not stayed in the current sector at least before the location registration information is generated. Therefore, the time that the terminal stayed in the current sector before the location registration information is generated is considered as 0, and the type information (generation factor) of the first location data is “location registration area boundary straddle” For example, the correction coefficient β in the above equation (13) (that is, the correction coefficient β related to the time difference from the previous position data) can be set to zero. As a result, it is possible to calculate a feature amount that is more realistic.

  As described above, when the feature amount calculation unit 103 calculates the feature amount for the target position data (first position data), the second and third positions which are position data before and after the first position data. The time difference between the second position data and the third position data is corrected according to the type information about the data (for example, the generation factor of the position data), and the feature amount is calculated using the corrected time difference. Thereby, the feature amount can be calculated with higher accuracy based on the type information of the position data.

  The above-described feature amount indicates the percentage of mobile devices present in the observation period T. FIG. 19 shows another method for measuring the number of mobile devices (mobile device amount) based on the position registration signal. There is a way. FIG. 19 is a schematic diagram illustrating the relationship between the transmission timing of the location registration signal transmitted by the mobile device 100 and the estimated stay period.

As shown in FIG. 19, a series of position information related to the same mobile device acquired within the extension period is arranged in chronological order in the order in which the position information (for example, position information based on the position registration signal) is acquired, When the in-area position data indicating the inside of the predetermined area with information is indicated by a black circle, and the outside-area position data indicating the outside of the predetermined area adjacent to the in-area position data is indicated by a white circle, respectively, the earliest area in time series and acquisition time t _a of the inner position data, estimates stay the corresponding time t _in the apportioning point (midpoint as an example) of the acquisition time t _x of the area outside the location data adjacent to outermost destination area position data The start time of Similarly, when the acquisition time t _c of the last area position data on the sequence, the proportional division point between the acquisition time t _y of the area outside the location data adjacent to area position data of said last (midpoint as an example) The corresponding time t _out is the end time of the estimated stay period. Accordingly, a period indicated by a rectangle in FIG. 19, that is, a period from time t _in to time t _out is calculated as the estimated stay period of the mobile device.

And the terminal number acquisition part 604 extracts the mobile device with which the estimated estimated stay period and the total time zone overlap. FIG. 20 is an explanatory diagram showing a specific example in the method for measuring the number of entrances. FIG. 20 shows a diagram in which a series of position information related to the mobile devices 100a to 100e acquired within the expansion period is arranged in the horizontal axis direction in time series in order of acquisition time of the position information for each mobile device. In FIG. 20, the in-area position data indicating the inside of the predetermined area is indicated by a black circle, and the outside-area position data indicating the outside of the predetermined area is indicated by a white circle. If attention is paid only to the position data in the area where the acquisition time of the position information is within the total time zone (time t ₀ to time t ₁ ), 3 of the mobile devices 100c to 100e that can detect the black circle in the predetermined area. Are extracted.

  However, as the mobile devices in which the estimated stay period indicated by the rectangle in FIG. 20 and the total time zone (time t0 to time t1) overlap, there are five mobile devices 100a to 100e. The five mobile devices 100a to 100e are extracted. As a result, it is estimated that the mobile devices 100a and 100b that are excluded when focusing only on the in-area position data indicated by the black circles are also staying in a predetermined sector somewhere in the counting time period. Become.

  In addition, the method using the estimated stay period described above is an example, and other methods may be adopted. That is, you may measure the black dot which overlapped as a position registration signal as the number of mobile machines as it is. In the example of FIG. 20, in that case, the number of mobile devices can be measured as nine.

  Next, processing for calculating a population distribution using the number of terminals calculated in this way will be described. FIG. 9 is a flowchart showing a population distribution calculation method of the information analysis apparatus 600 according to this embodiment. First, the operator of the information analysis apparatus 600 designates the sector i and the attribute j to be investigated (S101). This designation method may be performed by an operator operation using an input unit (not shown) or the like, and i and j may be automatically designated sequentially.

The number-of-terminals acquisition unit 604 calculates the number of terminals L _ij of the attribute j in the sector i (S102). The number of terminals L _ij is a numerical value obtained by summing up the number of feature amounts or location registration signals as described above, and the number of feature amounts or location registration signals summed up for each attribute j in the sector i at the time of summing up the number of terminals. It is derived based on. It should be noted that the feature amount of only sector i and attribute j may be calculated in advance at the time of feature amount calculation and used. Further, the terminal number acquisition unit 604 calculates the number of terminals Σ (j) L _ij of all attributes in the sector i (S103).

Next, the number of incoming calls C _i in the sector _i is acquired by the number of incoming calls acquisition unit 602 (S104). Then, the number-of-services acquisition unit 603 acquires information on sectors constituting the position registration area LAa to which the sector i belongs (S105), and acquires the number of areas Na of the mobile devices 100 in the position registration area LAa. (S106).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the location registration area is calculated by the incoming call number acquisition unit 602 (S107). Then, the number of terminals L _ij and Σ (i) L _ij acquired by the terminal number acquisition unit 604, the number of areas Na acquired by the area number acquisition unit 603, and the number of incoming calls C acquired by the number of incoming calls acquisition unit 602 _{Based on i} and its sum Σ (i) C _i , a corrected distribution number S _ij that suppresses an error in the number of terminals generated when the mobile device straddles the location registration area is calculated by the population distribution calculation unit 611 ( S108). In this way, it is possible to calculate a population distribution that suppresses errors across location registration areas in attribute j in sector i.

Next, operational effects of the information analysis apparatus 600 according to the first embodiment will be described. According to the information analysis apparatus 600, the terminal number acquisition unit 604 acquires the terminal number L _ij of the specified attribute j in the sector i specified by the operator or the like, and the incoming call number acquisition unit 602 receives the incoming call in the sector i. The number C _i is acquired, and the in-zone number acquisition unit 603 acquires the in-zone number Na in the location registration area LAa which is a cluster including a plurality of sectors including the sector i. The population distribution calculation unit 611 then obtains the number of incoming calls C _i, the total number of incoming calls Σ (i) C _i obtained in cluster units, the number of visited areas Na in the cluster, and the number of obtained terminals. Based on L _ij , the population distribution is calculated by calculating the corrected number of terminals S _ij in the sector. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area.

Note that the number of terminals L _ij in the present embodiment includes the number of terminals indicating the number of mobile stations, the number of entrances, or a feature amount indicating the existence ratio in a certain observation time zone. ) Indicates the amount of terminals existing, that is, the number of terminals.

[Second Embodiment]
In the first embodiment described above, the location registration area LAa is set as one cluster as a cluster unit for counting the number of located areas. However, the present invention is not limited to this. For example, a set of a plurality of sectors adjacent to a sector or located within a predetermined range may be configured as one cluster.

  In the first embodiment, the number of areas acquired by the area number acquisition unit 603 and the number of terminals in the location registration area acquired by the terminal number acquisition unit 604 (the sum of the numbers of terminals in the sector) are There are some errors. In other words, the number of visited areas is a numerical value based on information managed in the HLR 700 and the exchange 400, but since the position information changes as the mobile station moves, the number of terminals also increases or decreases accordingly. However, when the power is turned off, the battery is exhausted, or the service area is out of service area, it may be determined that the service area is present as it is. That is, in the first embodiment, since the information stored in the HLR 700 or the like is used, it is possible to simplify the acquisition processing of the number of located areas, but there may be an error. In the second embodiment, since this is based on the number of terminals acquired by the terminal number acquisition unit 604, this error is absorbed.

  FIG. 10 is a conceptual diagram showing the relationship between sectors to be measured and clusters. FIG. 10A is a conceptual diagram when the location registration area in the first embodiment is a cluster, and FIG. 10B is a concept showing the cluster in the present embodiment (second embodiment). FIG. In FIG. 10A, the location registration area is composed of sectors A1 to A7, and sector A5 is the target sector.

  In FIG. 10B, regardless of the location registration area, a set of sectors composed of sectors A2 to A7 arranged adjacent to each other centering on the target sector A5 is defined as a cluster. As described above, by setting the sector arranged adjacent to the target sector A5 as a cluster serving as a base for calculating the population distribution, it is possible to more accurately reflect the tendency of the geographical incoming frequency. it can. In addition, depending on the sector, there is a possibility that the data will be over-counted or under-counted, but if it is formed across multiple location registration areas in one cluster, this over-counting and under-counting will be offset. Therefore, a more accurate population distribution can be calculated. In the second embodiment as well, as in the first embodiment, it does not preclude the formation of a cluster with sectors constituting the location registration area. In this case, the number of terminals in the cluster is calculated by adding the number of terminals acquired in units of sectors.

  Hereinafter, the information analysis apparatus 600a of the second embodiment will be described. FIG. 11 is a block diagram illustrating a functional configuration of the information analysis apparatus 600a. The information analysis device 600a includes an information analysis device communication control unit 601, an incoming call number acquisition unit 602, a cluster formation unit 603a, a terminal number acquisition unit 604, a location information storage unit 605, an incoming call number storage unit 606, a population distribution calculation unit 611, And an output unit 610. The difference from the information analysis apparatus 600 in the first embodiment is that a cluster forming unit 603 a is provided instead of the in-zone number acquiring unit 603.

  As shown in FIG. 10B, the cluster forming unit 603a forms a set of sectors arranged adjacent to the target sector as a cluster, and the number of terminals of the sectors in the cluster. This is the part that calculates the sum of There are various methods for forming this cluster, such as sectors including peripherals in addition to the adjacent ones. Here, a set of all sectors adjacent to the sectors is considered.

ここでＮｂは、クラスタｂにおける位置登録信号に基づいた端末数の総和であって、

で算出することができる。すなわち、クラスタが十分に広い範囲であれば、位置登録信号に基づいた端末数の影響を受けづらくなる。そして、セクタによっては、過大に集計されたり、過小に集計されたり、する場合がありえるが、一のクラスタにおいて複数の位置登録エリアにまたがって形成されるなど、クラスタを十分広く取った場合、この過大集計および過小集計は相殺されることになり、実質的にＬ_ｉｊの総数と、Ｓ_ｉｊの総数とはほぼ同じであると考えることができる。 Then, the population distribution calculation unit 611 calculates the population distribution based on the number of sectors in the cluster formed by the cluster formation unit 603a. Specifically, the corrected terminal number S _ij is calculated based on the following equation (14).

Here, Nb is the total number of terminals based on the location registration signal in cluster b,

Can be calculated. That is, if the cluster is in a sufficiently wide range, it is difficult to be affected by the number of terminals based on the location registration signal. And depending on the sector, it may be over-counted or over-calculated, but if the cluster is wide enough, such as formed across multiple location registration areas in one cluster, this The over-total and the under-total are offset, and it can be considered that the total number of L _{ij and} the total number of S _ij are substantially the same.

  Next, processing of the information analysis apparatus 600a will be described. FIG. 12 is a flowchart showing the processing of the information analysis apparatus 600a in the second embodiment.

  First, the operator of the information analysis apparatus 600 designates the sector i and the attribute j to be investigated (S101). This designation method is performed by an operation by an operator using an input unit (not shown) or the like.

The number-of-terminals acquisition unit 604 calculates the number of terminals L _ij of the attribute j in the sector i (S102). The method for calculating the number of terminals L _ij is as described above. Further, the terminal number acquisition unit 604 calculates the number of terminals Σ (j) L _ij of all attributes in the sector i (S103).

Next, the number of incoming calls C _i in the sector _i is acquired by the number of incoming calls acquisition unit 602 (S104). Then, the cluster forming unit 603a acquires all sectors adjacent to the sector i as the cluster LAb (S105), and acquires the number of terminals Nb of the mobile device 100 located in the sector belonging to the cluster LAb (S106). ).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the cluster LAb is calculated by the incoming call number acquisition unit 602 (S107a). Then, the number of terminals L _ij and Σ (i) L _ij acquired by the terminal number acquisition unit 604, the number of terminals Nb acquired by the cluster formation unit 603a, the number of incoming calls C _i acquired by the number of incoming calls acquisition unit 602, and the number thereof Based on the total Σ (i) C _i , the population distribution calculation unit 611 calculates a corrected number of terminals S _ij that suppresses an error in the number of terminals caused by the mobile device straddling the location registration area (S108). In this way, it is possible to calculate a population distribution that suppresses errors across location registration areas in attribute j in sector i. In other words, depending on the sector, there may be cases of over-counting or under-counting, but this over-counting and under-counting is offset when formed across multiple location registration areas in one cluster. Therefore, a more accurate population distribution can be calculated.

  Next, a modification of the second embodiment will be described. In the information analysis apparatus 600a described above, the cluster forming unit 603a uses a cluster including all sectors adjacent to the target sector as one cluster and calculates the population distribution. ,there are various things. As described above, one cluster including all sectors adjacent to the target sector may be set as one cluster, or a predetermined number of sectors on the right side may be set as one cluster based on the target sector. For example, FIG. 13 shows an example. In FIG. 13, when there is a target sector A5, the sectors A1 to A7 existing on the right side are treated as one cluster. Then, when calculating the population distribution, the population distribution according to the cluster formation method is once calculated, and the average may be taken as described later.

  FIG. 14 is a flowchart showing a process of the information analysis apparatus 600a that calculates the population distribution by calculating an average of the population distribution calculated in the plurality of clusters, which is a process in the modification.

The number-of-terminals acquisition unit 604 calculates the number of terminals L _ij of the attribute j in the sector i (S102). The method for calculating the number of terminals L _ij is as described above. Further, the terminal number acquisition unit 604 calculates the number of terminals Σ (j) L _ij of all attributes in the sector i (S103).

Next, the number of incoming calls C _i in the sector _i is acquired by the number of incoming calls acquisition unit 602 (S104). The cluster forming unit 603a first sets n = 1, and acquires the cluster LAbn for the sector i by the cluster forming method n = 1 (S105a). The cluster forming unit 603a acquires all sectors adjacent to the sector i as the cluster LAbn (S105), and acquires the number Nb of the mobile devices 100 located in the sectors belonging to the cluster (S106).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the location registration area is calculated by the incoming call number acquisition unit 602 (S107a). Then, the terminal number L _ij and Σ (i) L _ij acquired by the terminal number acquiring unit 604, the terminal number Nbn acquired by the cluster forming unit 603a, the incoming call number C _i acquired by the incoming call number acquiring unit 602, and its Based on the total Σ (i) C _i , the population distribution calculation unit 611 calculates a corrected terminal number Sn _ij that suppresses an error in the number of terminals that occurs when the mobile device straddles the location registration area (S108).

Then, whether or not n has reached the predetermined number is determined by the cluster forming unit 603a. If it is less than the predetermined number, n = n + 1 is set (S110), and the cluster based on the cluster forming method n = 2 is determined. Obtained (S105a). Hereinafter, a plurality of Sn _ij is calculated until n reaches a predetermined number (S106 to S109). Then, the population distribution calculation unit 611 divides each corrected terminal number Sn _ij by a predetermined number, thereby calculating the average value and calculating the final corrected terminal number S _ij .

In this way, the final corrected terminal number S _ij is calculated based on the plurality of corrected terminal numbers Sn _ij , and the more accurate terminal number, that is, the population distribution can be calculated. .

ｎは、クラスタを示す。 In this modification, the corrected number of terminals S _ij is calculated for each cluster in which the set of sectors is changed, and the average value of the corrected number of terminals S _ij is used as the population distribution information. However, the present invention is not limited to this. Instead, the average value (Nbn / n) of the number of terminals Nbn calculated for each cluster is calculated without calculating the corrected number of terminals S _ij for each cluster, and the number of incoming calls for each cluster (Σ (i) An average value (Σ (n) Σ (i) Cin / n) of C _i ) may be calculated, and the corrected number of terminals S _ij may be calculated using these average values. Thereby, calculation processing is reduced.

n indicates a cluster.

Next, operational effects of the information analysis apparatus 600a of the present embodiment will be described. According to the information analysis apparatus 600a, the terminal number acquisition unit 604 acquires the terminal number L _ij of the specified attribute j in the sector i specified by the operator or the like, and the incoming call number acquisition unit 602 receives the incoming call in the sector i. The number C _i is acquired, and the cluster forming unit 603a acquires the number of terminals Nb in a plurality of clusters LAb including a plurality of sectors located in the vicinity of the sector i. The population distribution calculation unit 611 then obtains the number of incoming calls C _i, the total number of incoming calls Σ (i) C _i obtained in units of clusters, the number of terminals Nb in the cluster, and the number of obtained terminals L Based on _ij , the population distribution in the sector (that is, the corrected number of terminals S _ij ) is calculated. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area.

Furthermore, according to the information analysis apparatus 600a in the modification, the terminal number acquisition unit 604 acquires the terminal number L _ij of the specified attribute j in the sector i specified by the operator or the like, and the incoming call number acquisition unit 602 , it acquires the incoming number _{C i} in a sector i. The cluster forming unit 603a acquires the number of terminals Nbn (n is 1 to x) in the cluster LAb composed of a plurality of sectors located in the vicinity of the sector i. The population distribution calculation unit 611 then obtains the number of incoming calls C _i, the total number of incoming calls Σ (i) C _i obtained in cluster units, and the number of terminals Nbn (n is 1 to x) in the cluster. Based on the obtained terminal number L _ij , the corrected terminal number Sn _ij in each sector is calculated. Then, the population distribution calculation unit 611 calculates the average value (corrected terminal number S _ij ) by dividing the corrected terminal number Sn _ij by n to obtain population distribution information. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area. As in the first embodiment, the number of terminals L _ij is information including the amount of terminals.

Instead of calculating the average value of the corrected terminal number Sn _ij , the average value of the terminal number Nbn may be calculated, and the corrected terminal number S _ij may be calculated based on the average value.

According to the information analysis apparatus 600a according to this modification, the total number of incoming calls for each cluster unit is calculated for each set unit of different sectors, and the population distribution is calculated based on the average value of the total number of incoming calls. Thus, the population distribution for each sector and attribute can be calculated with higher accuracy.

[Third embodiment]
Next, a third embodiment will be described. FIG. 15 is a block diagram showing a functional configuration of the information analysis apparatus 600b in the present embodiment. As shown in FIG. 15, the information analysis device 600 b includes an information analysis device communication control unit 601, an incoming call number acquisition unit 602, an in-zone number acquisition unit 603, a cluster formation unit 603 a, a terminal number acquisition unit 604, and a location information storage unit 605. The incoming call number storage unit 606, the population distribution calculation unit 611, and the output unit 610 are configured. The difference from the information analysis apparatus 600a in the second embodiment is that a location number acquisition unit 603 is further provided.

In this information analysis apparatus 600b, the number-of-areas acquisition unit 603 acquires the number-of-areas Na in the location registration area LAa to which the target sector belongs, and the cluster formation unit 603a acquires the number of areas-based. The number of terminals Nb in the cluster LAb formed in this way is acquired. Then, the population distribution calculation unit 611 calculates the corrected terminal numbers Sa _ij and Sb _ij based on the number of areas Na and Nb and the number of terminals L _ij , and calculates the average value as population distribution information. .

  Next, processing of the information analysis apparatus 600b configured as described above will be described. FIG. 16 is a flowchart showing the processing of the information analysis apparatus 600b.

First, the operator of the information analysis apparatus 600 designates the sector i and the attribute j to be investigated (S101). The number-of-terminals acquisition unit 604 calculates the number of terminals L _ij of the attribute j in the sector i (S102). The method for calculating the number of terminals L _ij is as described above. Further, the terminal number acquisition unit 604 calculates the number of terminals Σ (j) L _ij of all attributes in the sector i (S103).

Next, the number of incoming calls C _i in the sector _i is acquired by the number of incoming calls acquisition unit 602 (S104). Then, the location number acquisition unit 603 acquires the location registration area LAa to which the sector i belongs (S105), and acquires the location number Na of the mobile device 100 in the location registration area LAa (S106).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the location registration area LAa is calculated by the incoming call number acquisition unit 602 (S107). Then, the number of terminals L _ij and Σ (i) L _ij acquired by the terminal number acquisition unit 604, the number of areas Na acquired by the area number acquisition unit 603, and the number of incoming calls C acquired by the number of incoming calls acquisition unit 602 _{Based on i} and its sum Σ (i) C _i , a corrected terminal number Sa _ij that suppresses an error in the number of terminals caused by the mobile device straddling the location registration area is calculated by the population distribution calculation unit 611 ( S108).

  Next, the cluster forming unit 603a acquires all sectors adjacent to the sector i as the cluster LAb (S112), and acquires the number Nb of the mobile devices 100 located in the sectors belonging to the cluster (S113). ).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the cluster LAb is calculated by the incoming call number acquisition unit 602 (S114). Then, the number of terminals L _ij and Σ (i) L _ij acquired by the terminal number acquisition unit 604, the number of terminals Nb acquired by the cluster formation unit 603a, the number of incoming calls C _i acquired by the number of incoming calls acquisition unit 602, and the number thereof Based on the total Σ (i) C _i , the corrected terminal number Sb _ij is calculated by the population distribution calculation unit 611 (S115).

Then, to calculate the number of terminals _{S ij} is corrected number of terminals _{Sa ij} and the average value of the corrected number of terminals _{Sb ij,} which is referred to as population distribution information (S116). Thereby, population distribution information with higher accuracy can be calculated.

In this embodiment, the corrected number of terminals Sa _ij based on the location registration area LAa and the corrected number of terminals Sb _ij based on the cluster LAb composed of a set of sectors adjacent to or adjacent to the target sector. The average value is calculated as the corrected number of terminals S _{ij and is used} as the population distribution information. However, the present invention is not limited to this, and the number of visited areas Na based on the location registration area LAa and the total number of incoming calls Σ (i) C _i and the cluster LAb composed of a set of sectors adjacent to or adjacent to the target sector The average value of the number of terminals Nb and the total number of incoming calls Σ (i) C _i may be calculated, and the corrected number of terminals S _ij may be calculated based on these average values.

Next, operational effects of the information analysis apparatus 600b of the third embodiment will be described. According to the information analysis apparatus 600b, the terminal number acquisition unit 604 acquires the terminal number L _ij of the specified attribute j in the sector i specified by the operator or the like, and the incoming call number acquisition unit 602 receives the incoming call in the sector i. The number acquisition section 603 obtains the number C _i , and the location number acquisition unit 603 includes a location registration area LAa that is a cluster composed of a plurality of sectors including the sector i and a predetermined cluster LAb that includes sectors adjacent to the sector i. The number Na of terminals and the number Nb of terminals are acquired. Then, the population distribution calculation unit 611 has acquired the number of incoming calls C _i , the total number of incoming calls Σ (i) C _{i in} which the number of incoming calls is acquired in units of clusters, and the number of areas Na and Nb in the cluster. Based on the terminal numbers L _ij and L _ij , the corrected terminal numbers Sa _ij and Sb _ij in the sector are calculated, and the average value is calculated as the population distribution. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area.

[Fourth embodiment]
Next, a fourth embodiment will be described. FIG. 17 is a block diagram illustrating a functional configuration of the information analysis apparatus 600c according to the fourth embodiment. As shown in FIG. 17, the information analysis device 600 c includes an information analysis device communication control unit 601, an incoming call number acquisition unit 602, an in-zone number acquisition unit 603, an apportioned processed terminal number acquisition unit 604 a, a location information storage unit 605, and the number of incoming calls. A storage unit 606, a population distribution calculation unit 611, and an output unit 610 are included. The difference from the information analysis apparatus 600 in the first embodiment is that a proportionally processed terminal number acquisition unit 604a is provided instead of the terminal number acquisition unit 604.

  The apportioned processed terminal number acquisition unit 604a is calculated by the location registration signal processing unit 404 and based on the subscriber information / location information stored in the HLR 700 or the exchange 400, the feature amount in the designated sector and attribute When the mobile device 100 outputs a location registration signal by acquiring the location number (first terminal number) expressed by the following, the location registration signal is removed and removed. Based on the location registration signal (that is, based on the location registration signal periodically output, the location registration signal when the power is turned on, or the location registration signal output when returning from outside the service area, etc.) This is a part for acquiring the number of terminals (second terminal number: the number of terminals after removal). And based on these 1st terminal numbers and 2nd terminal numbers, the 3rd terminal number used as the apportioned value is calculated, and this is made into the number of terminals.

  As described above, the HLR 700 and the exchange 400 store the location information of the mobile device 100 based on the location registration signal transmitted from the mobile device 100. In general, the mobile device 100 is configured to transmit a location registration signal when a predetermined condition is satisfied, for example, when the location registration area is crossed, when a predetermined period is reached, or when the power is turned on. Yes. When the location registration area is crossed, the location registration signal is attached with information indicating that the location registration processing has been executed by crossing the location registration area. In the HLR 700 and the exchange 400, Location information is stored together with the information.

  The apportioned terminal count acquisition unit 604a removes the location information indicating that the location registration processing has been performed by straddling the location registration area, and thus based on the location registration signal and the like periodically performed. Position information can be acquired, and the number of position information based on this can be acquired.

The apportioned terminal count acquisition unit 604a can calculate the apportioned terminal count by applying the first terminal count and the second terminal count to the following expression.
Number of terminals L _ij = pM _ij + (1−p) N _ij (18)
Here, M _ij is the first number of terminals with attribute j in sector i, and N _ij is the second number of terminals with attribute j in sector i, excluding location information across the location registration area. The number of terminals. Moreover, p is an arbitrary value and is not limited to 1 or less.

In this way, the number of terminals L _ij calculated by the apportioned terminal number acquisition unit 604a, the number of areas Na acquired by the area number acquisition unit 603, and the number of incoming calls C _i acquired by the number of incoming calls acquisition unit 602 are set. Based on this, the population distribution can be calculated. Note that the number Nb of terminals based on the cluster LAb may be calculated instead of the number Na of areas, and the population distribution may be calculated based on this.

  Next, processing of the information analysis apparatus 600c will be described. FIG. 18 is a flowchart showing the processing of the information analysis apparatus 600c.

  First, the operator of the information analysis apparatus 600 designates the sector i and the attribute j to be investigated (S101).

Next, the apportioned processed terminal number acquisition unit 604a acquires the first terminal number M _ij and the second terminal number N _ij of the attribute j in the sector i (S102a and S102b). Then, the number of terminals L _ij is calculated by calculating the number of terminals L _ij = pM _ij + (1−p) M _ij (S102c). Further, the terminal number acquisition unit 604 calculates the number of terminals Σ (j) L _ij of all attributes in the sector i (S103).

Next, the number of incoming calls C _i in the sector _i is acquired by the number of incoming calls acquisition unit 602 (S104). Then, the number-of-services acquisition unit 603 acquires information on sectors constituting the position registration area LAa to which the sector i belongs (S105), and acquires the number of areas Na of the mobile devices 100 in the position registration area LAa. (S106).

Next, the total number of incoming calls Σ (i) C _{i of} all sectors belonging to the location registration area is calculated by the incoming call number acquisition unit 602 (S107). Then, the number of terminals L _ij and Σ (i) L _ij acquired by the terminal number acquisition unit 604, the number of areas Na acquired by the area number acquisition unit 603, and the number of incoming calls C acquired by the number of incoming calls acquisition unit 602 _{Based on i} and its sum Σ (i) C _i , a corrected distribution number S _ij that suppresses an error in the number of terminals generated when the mobile device straddles the location registration area is calculated by the population distribution calculation unit 611 ( S108). In this way, it is possible to calculate a population distribution that suppresses errors across location registration areas in attribute j in sector i.

In the fourth embodiment, the population distribution is calculated based on the number Na of areas acquired using the area number acquisition unit 603. However, the present invention is not limited to this. The number of terminals L _ij calculated by the above equation (14) can be applied to the information analysis devices 600a and 600b described in the third embodiment. In addition, when applied to the second embodiment, a sector may be selected so that the value of the number of terminals L _ij is relatively large in the cluster range with respect to the total number of terminals. desirable. Further, in order to reduce the influence due to the deviation of the incoming frequency of geographic calls, it is desirable that the cluster area be within the narrowest possible geographical range. Therefore, it is desirable to select sectors adjacent to sector i so as to fit in one cluster.

Next, operational effects of the information analysis apparatus 600c of the present embodiment will be described. According to the information analysis device 600c, the terminal number acquisition unit 604 acquires the number of terminals M _ij of the specified attribute j in the sector i specified by the operator or the like, and straddles the location registration area from the number of terminals M _ij . The terminal number N _ij from which the terminal subjected to the location registration process is removed is acquired. Then, the terminal number Lij is calculated based on these terminal numbers Mij and Nij.

Then, call number obtaining unit 602 obtains the incoming number C _i in a sector i, the located number acquisition unit 603 acquires the located number Na in the location registration area LAa containing the sector i. Then, the population distribution calculation unit 611 acquires the acquired number of incoming calls C _i, the total number of incoming calls Σ (i) C _i obtained for the number of incoming calls in cluster units, and the number of areas Na in the location registration area LAa. Based on the number of terminals L _ij , the population distribution in the sector (that is, the corrected number of terminals S _ij ) is calculated. Thereby, it is possible to accurately calculate the population distribution for each sector and for each attribute, which eliminates the influence of straddling the location registration area.

DESCRIPTION OF SYMBOLS 10 ... Communication system, 14 ... Input device, 15 ... Output device, 16 ... Communication module, 17 ... Auxiliary storage device, 17 ... Feature-value calculation part, 100 ... Mobile device, 302 ... Communication control part, 303 ... Position registration signal reception , 304 ... Signal number measurement unit, 400 ... Exchanger, 401 ... Exchange communication control unit, 402 ... Conversion unit, 403 ... Storage unit, 404 ... Location registration signal processing unit, 500 ... Management center, 501 ... Social sensor unit, 502 ... Petamining unit, 503 ... Mobile demography unit, 504 ... Visualization solution unit, 600, 600a, 600b, 600c ... Information analysis device, 601 ... Information analysis device communication control unit, 602 ... Number of incoming calls acquisition unit, 603 ... Location Number acquisition unit, 603a ... Cluster formation unit, 604 ... Terminal number acquisition unit, 604a ... Apportioned terminal count acquisition , 605 ... position information storage unit, 606 ... incoming number storage section, 610 ... output unit, 611 ... population distribution calculator.

Claims (11)

Signal number storage means for storing the number of signals communicated with the terminal for each sector;
A terminal number acquisition means for acquiring the number of terminals having a specified attribute in a specified sector and the number of terminals having all the attributes in the sector ;
An incoming number acquisition means for acquiring the incoming number in the designated sector from the signal number storage means;
A number-of-services acquisition means for acquiring the number of services in a cluster composed of a plurality of sectors including the designated sector;
The number of incoming calls acquired by the number of incoming calls acquisition means, the total number of incoming calls obtained by summing up the number of incoming calls in units of clusters, the number of areas in the cluster, and the designated number acquired by the number of terminal acquisition means An information analysis apparatus comprising: a population distribution calculating unit that calculates a population distribution in the sector based on the number of terminals having a specified attribute in the sector and the number of terminals having all the attributes in the sector.   2. The information analysis apparatus according to claim 1, wherein the in-zone number acquisition unit acquires the in-zone number totaled for each location registration area including the designated sector.   2. The information analysis apparatus according to claim 1, wherein the in-zone number acquisition unit acquires the in-zone number totaled in a cluster unit composed of sectors adjacent to the designated sector. The in-zone number acquisition means includes
2. The information according to claim 1, wherein the number of located areas is a cluster unit including a plurality of sectors including the designated sector, and the total number of areas is aggregated in a cluster unit across a plurality of location registration areas. Analysis equipment. The in-zone number acquisition means forms a plurality of cluster units including sectors different in combination with the designated sector, and calculates the in-zone numbers of the cluster units,
The information analysis apparatus according to claim 1, wherein the population distribution calculation unit calculates the population distribution based on an average value of the number of the locations. The in-zone number acquisition means forms a plurality of cluster units including sectors different in combination with the designated sector, and calculates the in-zone numbers of the cluster units,
The information analysis apparatus according to claim 1, wherein the population distribution calculating unit calculates a population distribution in each cluster unit by using the number of located areas, and calculates an average value of the population distribution. The population distribution calculating means includes the number of in-zones counted in cluster units including peripheral sectors that are further adjacent to adjacent sectors adjacent to the specified sector and are not directly adjacent to the specified sector, and the specified The information analysis apparatus according to claim 1, wherein the population distribution is calculated based on the number of areas that are aggregated in a set unit of sectors constituting the location registration area including the sector.   The said number of area acquisition means acquires the number of areas in the said cluster by acquiring the number of terminals in each sector which comprises a cluster, and totaling this, The one in any one of Claim 3 to 7 characterized by the above-mentioned. The information analysis apparatus according to one item. The terminal number acquisition means includes
Calculate the number of terminals after removal by removing the number of terminals with the specified attribute in the specified sector and the set of terminals registered by crossing the location registration area from the set of terminals with the specified attribute And
The information analysis apparatus according to any one of claims 1 to 8, wherein a numerical value that is apportioned at a predetermined ratio between the number of terminals and the number of terminals after removal is used as the number of terminals. The population distribution calculating means includes
It is assumed that the number of mobile terminals Lij of the attribute j in the designated sector i is
Let N be the number of areas in the cluster,
If the number of incoming calls in sector i is Ci,

The information analysis apparatus according to claim 1, wherein the population distribution is calculated by calculating A terminal number obtaining step for obtaining the number of terminals having the designated attribute in the designated sector and the number of terminals having all the attributes in the sector ;
For each sector, from the number-of-signals storage means for storing the number of signals communicated with the terminal, the number-of-arrivals obtaining step for obtaining the number of incoming calls in the designated sector;
A number- of- services acquisition step of acquiring a number of service areas in a cluster composed of a plurality of sectors including the designated sector;
The number of incoming calls acquired by the step of acquiring the number of incoming calls, the total number of incoming calls obtained by counting the number of incoming calls in units of clusters, the number of in- zones in the cluster, and the designated sector acquired by the step of acquiring the number of terminals. A population distribution calculating method comprising: a population distribution calculating step of calculating a population distribution in the sector based on the number of terminals with the specified attribute in and the number of terminals with all the attributes in the sector.

Images