WO2020184246A1 - 無線端末収容判定装置、無線端末収容判定方法、及びプログラム - Google Patents

無線端末収容判定装置、無線端末収容判定方法、及びプログラム Download PDF

Info

Publication number
WO2020184246A1
WO2020184246A1 PCT/JP2020/008607 JP2020008607W WO2020184246A1 WO 2020184246 A1 WO2020184246 A1 WO 2020184246A1 JP 2020008607 W JP2020008607 W JP 2020008607W WO 2020184246 A1 WO2020184246 A1 WO 2020184246A1
Authority
WO
WIPO (PCT)
Prior art keywords
density function
probability density
base station
received power
interference signal
Prior art date
Application number
PCT/JP2020/008607
Other languages
English (en)
French (fr)
Inventor
一光 坂元
鈴木　賢司
陽平 片山
洋輔 藤野
浩之 福本
Original Assignee
日本電信電話株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to US17/435,093 priority Critical patent/US11785475B2/en
Publication of WO2020184246A1 publication Critical patent/WO2020184246A1/ja

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W28/00—Network traffic management; Network resource management
    • H04W28/02—Traffic management, e.g. flow control or congestion control
    • H04W28/0231—Traffic management, e.g. flow control or congestion control based on communication conditions
    • H04W28/0236—Traffic management, e.g. flow control or congestion control based on communication conditions radio quality, e.g. interference, losses or delay
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/18—Network planning tools
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/20—Monitoring; Testing of receivers
    • H04B17/26—Monitoring; Testing of receivers using historical data, averaging values or statistics
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H04B17/318—Received signal strength
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H04B17/345—Interference values
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/373—Predicting channel quality or other radio frequency [RF] parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/22—Traffic simulation tools or models
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/02—Arrangements for optimising operational condition
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/06—Testing, supervising or monitoring using simulated traffic
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/08—Testing, supervising or monitoring using real traffic
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W28/00—Network traffic management; Network resource management
    • H04W28/02—Traffic management, e.g. flow control or congestion control
    • H04W28/0247—Traffic management, e.g. flow control or congestion control based on conditions of the access network or the infrastructure network
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W28/00—Network traffic management; Network resource management
    • H04W28/16—Central resource management; Negotiation of resources or communication parameters, e.g. negotiating bandwidth or QoS [Quality of Service]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W52/00—Power management, e.g. TPC [Transmission Power Control], power saving or power classes
    • H04W52/04—TPC
    • H04W52/18—TPC being performed according to specific parameters
    • H04W52/22—TPC being performed according to specific parameters taking into account previous information or commands
    • H04W52/225—Calculation of statistics, e.g. average, variance
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • H04W64/003—Locating users or terminals or network equipment for network management purposes, e.g. mobility management locating network equipment
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W84/00—Network topologies
    • H04W84/02—Hierarchically pre-organised networks, e.g. paging networks, cellular networks, WLAN [Wireless Local Area Network] or WLL [Wireless Local Loop]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W88/00—Devices specially adapted for wireless communication networks, e.g. terminals, base stations or access point devices
    • H04W88/02—Terminal devices
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W88/00—Devices specially adapted for wireless communication networks, e.g. terminals, base stations or access point devices
    • H04W88/08—Access point devices
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W88/00—Devices specially adapted for wireless communication networks, e.g. terminals, base stations or access point devices
    • H04W88/18—Service support devices; Network management devices
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W92/00—Interfaces specially adapted for wireless communication networks
    • H04W92/02—Inter-networking arrangements
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W92/00—Interfaces specially adapted for wireless communication networks
    • H04W92/04—Interfaces between hierarchically different network devices
    • H04W92/10—Interfaces between hierarchically different network devices between terminal device and access point, i.e. wireless air interface
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W92/00—Interfaces specially adapted for wireless communication networks
    • H04W92/16—Interfaces between hierarchically similar devices
    • H04W92/20—Interfaces between hierarchically similar devices between access points
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Definitions

  • the present invention relates to a wireless terminal accommodation determination device, a wireless terminal accommodation determination method, and a program.
  • IoT devices Internet of Things
  • terminals things with various sensors and communication functions
  • IoT devices Remote control of IoT is becoming more and more popular every year.
  • LPWA Low Power Wide Area network
  • LPWA includes, for example, LoRaWAN and Sigfox that use an unlicensed band (frequency band that does not require a radio station license), and LTE-M (Long Term Evolution for Machines) that uses a license band (frequency band that does not require a radio station license).
  • NB Narrow Band
  • IoT devices are often fixedly installed and used mainly in specific locations.
  • radio wave propagation is highly location-dependent.
  • the IoT device is installed in a dead zone where the received power of radio waves from the base station is low, communication may be disabled for a long period of time. Therefore, it is important to design the station so that sufficient received power can be secured at the planned installation location of the IoT device.
  • Patent Document 1 As a conventional technique for estimating the communication success rate, for example, there are techniques described in Patent Document 1 and Patent Document 2.
  • the technique described in Patent Document 1 is a technique for measuring a communication success rate by temporarily installing a radio and actually trying communication.
  • Patent Document 2 is a technique for estimating the throughput in a multi-cell environment composed of a plurality of base stations by using Monte Carlo simulation.
  • Patent Document 1 has a problem that a large amount of human cost is required for measurement preparation and measurement work because a terminal is temporarily installed at a planned installation location and the communication success rate is measured in an actual environment.
  • a wide area wireless communication system such as LPWA
  • the human cost required for measuring the communication success rate in a real environment becomes enormous.
  • Patent Document 2 estimates the throughput by Monte Carlo simulation, so that measurement in a real environment is not necessary.
  • the technique described in Patent Document 2 repeatedly executes a process of probabilistically generating traffic of a large number of terminals using random numbers a huge number of times (for example, tens of thousands of times or more), thereby causing peripheral terminals.
  • the propagation loss between the base station and the terminal is calculated in consideration of the distance attenuation (long section variation) and shadowing (short section variation) based on the propagation model.
  • the process of calculating the received power instantaneous value of the desired signal and the interference signal by applying the influence of the instantaneous fluctuation using a random number is enormous. It needs to be repeated a number of times (for example, tens of thousands or more). Therefore, a huge amount of calculation time is required.
  • the present invention has been made in consideration of such circumstances, and is a technique for estimating the communication success rate at the planned installation location of a new wireless terminal at a higher speed without performing measurement in a real environment or Monte Carlo simulation.
  • the purpose is to provide.
  • One aspect of the present invention is an average received power calculation unit that calculates the average received power for each of the minute meshes constituting the area for each of the plurality of base stations, and a terminal belonging to each base station based on the average received power. Is determined, and for each base station, the probability density function of the interference signal power between the base stations is shown by substituting the average received power of the interference signal from the base station into the probability density function of the instantaneous value of the received power. A second probability density function showing the probability density function of the interference signal power from each terminal by creating a density function and substituting the average received power of the interference signal from each terminal into the probability density function of the instantaneous value of the received power.
  • the interference signal power CDF creation unit that creates the cumulative distribution function of the total interference signal power and the average received power of the desired signal from the desired terminal for each base station are substituted into the probability density function of the received power instantaneous value.
  • the desired signal power PDF creation unit that creates the fourth probability density function indicating the probability density function of the desired signal power, and the desired terminal of the desired terminal based on the cumulative distribution function and the fourth probability density function for each base station. It is a wireless terminal accommodation determination device including a communication success rate estimation unit for calculating a communication success rate.
  • one aspect of the present invention is the above-mentioned wireless terminal accommodation determination device, in which the interference signal power CDF creation unit obtains an instantaneous value of received power of a signal arriving at a base station from all terminals including the desired terminal.
  • the interference signal power CDF creation unit obtains an instantaneous value of received power of a signal arriving at a base station from all terminals including the desired terminal.
  • one aspect of the present invention is the above-mentioned wireless terminal accommodation determination device, in which the average received power calculation unit further calculates the average interference signal power of inter-base station interference from the positional relationship of each base station, and the above-mentioned
  • the interference signal power CDF creation unit is based on the instantaneous value of the received power created in consideration of the influence of the instantaneous fluctuation on the average received power, the number of terminals, the traffic pattern, and the average interference signal power.
  • the 1st probability density function and the 2nd probability density function are created.
  • one aspect of the present invention is the above-mentioned wireless terminal accommodation determination device, further including a repeat determination unit, and the desired signal power PDF creation unit sorts the base stations in descending order of the calculated average received power. Then, for each iterative process, the number of base stations considered in the communication success rate estimation is increased in descending order of power, a desired signal power PDF is created for the base station to be calculated, and the communication success rate estimation unit is the cumulative.
  • the communication success rate is calculated based on the distribution function, the fourth probability density function of the base station to be calculated, and the thermal noise power, and the repetition determination unit determines that the calculated communication success rate is the first predetermined value.
  • one aspect of the present invention belongs to each base station based on the average received power calculation step of calculating the average received power for each of the minute meshes constituting the area and the average received power for each of the plurality of base stations.
  • the probability density function of the interference signal power between the base stations is shown by substituting the average received power of the interference signal from the base station into the probability density function of the instantaneous value of the received power. 1 Probability density function is created, and the average received power of the interference signal from each terminal is substituted into the probability density function of the instantaneous value of the received power to show the probability density function of the interference signal power from each terminal.
  • a density function is created, and by convoluting the first probability density function and the second probability density function, a third probability density function showing the probability density function of the total interference signal power is created, and the third probability density function is created.
  • Interference signal power CDF creation step to create a cumulative distribution function of total interference signal power based on, and for each base station, the average received power of the desired signal from the desired terminal is substituted into the probability density function of the received power instantaneous value.
  • the desired signal power PDF creation step for creating a fourth probability density function indicating the probability density function of the desired signal power, and the desired signal power PDF creation step for each base station based on the cumulative distribution function and the fourth probability density function. It is a wireless terminal accommodation determination method having a communication success rate estimation step for calculating a communication success rate of a terminal.
  • one aspect of the present invention is a program for operating a computer as the above-mentioned wireless terminal accommodation determination device.
  • the present invention it is possible to estimate the communication success rate at the planned installation location of the new wireless terminal at a higher speed without performing measurement in a real environment or Monte Carlo simulation.
  • the uplink communication from the terminal (IoT device) to the base station is performed by the Class A protocol.
  • the downlink communication from the base station to the terminal is performed by the Class A protocol when transmitting Ac (Acknowledge) for the uplink communication, and the data transferred from the application server (not shown) to the base station.
  • Ac Acknowledge
  • the data transferred from the application server (not shown) to the base station When transmitting, it is performed by the Class B protocol.
  • the signal transmitted from the terminal is received and demodulated by the plurality of base stations. Then, the successfully received signal is sent to a higher-level network server (not shown). Then, among the signals successfully received by the network server, the signal having the highest received power is selected (site diversity).
  • site diversity among the signals successfully received by the network server, the signal having the highest received power is selected (site diversity).
  • the base station has a plurality of antennas
  • the signal transmitted from the terminal is received by the plurality of antennas and demodulated. Then, among the signals successfully received by the base station, the signal having the highest received power is selected (antenna diversity).
  • the terminal for which the communication success rate is estimated is referred to as a “desired terminal”, and the other terminals are referred to as “interference terminals”.
  • a base station with which a desired terminal communicates is referred to as a “desired base station”, and other base stations are referred to as an “interference base station”.
  • FIG. 1 is a block diagram showing a functional configuration of the wireless terminal accommodation determination device 1a according to the first embodiment of the present invention.
  • the wireless terminal accommodation determination device 1a includes an average received power calculation unit 11, a radio equipment database 12, a desired signal power PDF creation unit 13a, an interference signal power CDF creation unit 14, and a communication success rate. It includes an estimation unit 15a.
  • the average received power calculation unit 11 corresponds to a general radio wave propagation simulator.
  • the average received power calculation unit 11 imports map data (for example, data such as terrain height, building height, land use classification, etc.) from an external device or the like.
  • the average received power calculation unit 11 divides the map data of the area designated by the operator into minute meshes (for example, meshes at intervals of 5 m). Then, the average received power calculation unit 11 calculates the average received power of uplink communication between the base station and each mesh (that is, a virtual terminal assumed to be installed at a position corresponding to each mesh). Further, the average received power calculation unit 11 calculates the average interference signal power of the interference between the base stations generated between the base stations based on the position of each base station.
  • the radio equipment database 12 holds parameters related to radio equipment such as the number of antennas and thermal noise power of each base station. In addition, the radio equipment database 12 holds parameters related to radio equipment such as the number of antennas, thermal noise power, and traffic pattern (communication frequency and communication time per communication) of each terminal.
  • the communication frequency may be the value obtained by dividing the communication frequency of the terminal by the number of channels.
  • the desired signal power PDF creation unit 13a considers the influence of instantaneous fluctuations on the average received power between the mesh on which the desired terminal is located and the base station, and considers the influence of instantaneous fluctuations on the probability density function (PDF) of the instantaneous value of the received power. ) Is created for each base station.
  • PDF probability density function
  • the interference signal power CDF creation unit 14 considers the influence of instantaneous fluctuations and the traffic pattern on the average received power between the mesh on which the interfering terminal is located and the base station, and considers the influence of the instantaneous fluctuation and the traffic pattern, and the probability density function (PDF) of the instantaneous value of the received power. : Probability Density Function) is created. Further, the interference signal power CDF creation unit 14 creates a PDF of the interference signal power based on the average received power between the desired base station and the interference base station and the traffic pattern of the interference terminal under the interference base station. The interference signal power CDF creation unit 14 totals the created PDFs to create a CDF. Hereinafter, the CDF created by totaling the PDFs of the interference signals is referred to as "interference signal power CDF".
  • the communication success rate estimation unit 15a is composed of the desired signal power PDF created by the desired signal power PDF creation unit 13a, the interference signal power CDF created by the interference signal power CDF creation unit 14, and the thermal noise power of the base station. , The probability that the desired signal power to interference and noise power ratio (SINR; Signal-to-Interference plus Noise power Ratio) is equal to or higher than the required SINR is calculated. Then, the communication success rate estimation unit 15a calculates the communication success rate based on the calculated probability in consideration of the effects of the site diversity and the antenna diversity.
  • the required SINR indicates a value of SINR required for successful communication in the wireless communication of the modulation method used.
  • FIG. 2 is a diagram for explaining a communication success rate estimation process by the wireless terminal accommodation determination device 1a according to the first embodiment of the present invention.
  • a case of estimating the communication success rate in the channel used only for the communication of Class A among the plurality of channels of LoRaWAN will be described.
  • Class A and Class B when the interference signal power CDF creation unit 14 creates the PDF of the interference signal power from the interference base station, Class The communication occurrence time of B may also be taken into consideration.
  • terminals # 1 and terminal # 6 are newly installed in an environment in which a wireless area is constructed by three base stations, communication between terminals # 1 and terminal # 6 The case of estimating the success rate will be described.
  • Each base station has one antenna.
  • the average received power calculation unit 11 corresponds to a general radio wave propagation simulator.
  • the average received power calculation unit 11 divides the area designated by the operator into minute meshes. Then, the average received power calculation unit 11 calculates the average received power of uplink communication between the base station and each mesh (that is, a virtual terminal assumed to be installed at a position corresponding to each mesh).
  • the average received power calculation unit 11 acquires the value of the average received power in the mesh in which the terminal #j is located, and the average received power Rup (i, j) for uplink communication between the base station #i and the terminal #j. And. Further, the average received power calculation unit 11 calculates the average interference signal power of the interference between the base stations generated between the base stations based on the position of each base station.
  • the average interference signal power of interference between base stations from base station #i to base station #k is expressed as RICI (i, k).
  • the interference signal power CDF creation unit 14 first creates a PDF of the interference signal power from the terminal (second probability density function) and a PDF of the interference signal power between the base stations (first probability density function).
  • the interference signal power CDF creation unit 14 creates a PDF of the instantaneous received power value for each terminal in consideration of the influence of instantaneous fluctuations and the traffic pattern on the average received power between the desired base station and the terminal of the desired terminal. By convolving a plurality of created PDFs, a PDF of interference signal power from a terminal (PDF of interference signals from all terminals at a desired base station of a desired terminal) is obtained.
  • the instantaneous value of received power is calculated by randomly determining and synthesizing the phases of radio waves arriving at base station #i by a plurality of routes using random numbers. Then, by repeating this, a PDF of the instantaneous value of the received power can be created.
  • this iterative calculation is generally performed an enormous number of times (about tens of thousands of times). Therefore, there is a problem that a huge amount of calculation time is required.
  • the instantaneous fluctuation due to fading is expressed by the following analysis formula.
  • the desired signal power PDF is obtained without performing Monte Carlo simulation.
  • the PDF of the amplitude a that fluctuates instantaneously due to Rayleigh fading can be expressed by the following equation (1) using Rup (i, j) (see Non-Patent Document 1). ).
  • the number of terminals is totaled assuming that each terminal belongs to the base station #i having the maximum Rup (i, j). Then, for each base station, the communication occurrence time of Ac (downlink communication by Class A) to the terminal is calculated based on the traffic pattern of the uplink communication of the terminal to which the base station belongs, and the calculated communication occurrence time out of the total time is calculated. During that time, the PDF is created assuming that interference between base stations of the interference power RICI (i, k) occurs.
  • base stations are generally installed at high positions such as on the rooftops of high-rise buildings, and it is common for base stations to be in an environment where there are no reflectors or obstacles in the vicinity and there is visibility. is there. Therefore, for the interference between base stations, the average interference signal power calculated based on the positional relationship between the base stations is used.
  • rice fading is assumed and the average interference signal calculated from the positional relationship between the base stations.
  • the interference signal power CDF creation unit 14 uses the interference base station # 2 and the interference base station # 3 to create the base station # 1.
  • the two PDFs of the interference signal power to the base station # 1 and the PDF of the interference signal power from the terminal to the base station # 1 are convoluted.
  • the interference signal power CDF creation unit 14 creates a PDF (defined as PDF I (r)) (third probability density function) of the total interference signal power, and the interference signal power CDF is created by the following equation (3).
  • PDF I (r) third probability density function
  • the interference signal power CDF creation unit 14 uses the interference base station # 1 and the interference base station # 3 as the base station. The two PDFs of the interference signal power to # 2 and the PDF of the interference signal power from the terminal to the base station # 2 are convoluted. As a result, the interference signal power CDF creation unit 14 creates a PDF (third probability density function) of the total interference signal power, and creates an interference signal power CDF (CDF I_2 (y)) (cumulative distribution function).
  • the interference signal power CDF creation unit 14 uses the interference base station # 1 and the interference base station # 2 as the base station. The two PDFs of the interference signal power to # 3 and the PDF of the interference signal power from the terminal to the base station # 3 are convoluted. As a result, the interference signal power CDF creation unit 14 creates a PDF (third probability density function) of the total interference signal power, and creates an interference signal power CDF (CDF I_3 (y)) (cumulative distribution function).
  • the desired signal power PDF creation unit 13a has the average received power Rup (1,1), Rup (2,1) of the desired signal from the terminal # 1 to the base station # 1, the base station # 2, and the base station # 3. ), and R Stay up-(a 3,1), in consideration of the influence of instantaneous variation, the desired signal power PDF (PDF S_1 (y), PDF S_2 (y), and creates a PDF S_3 (y)).
  • the communication success rate estimation unit 15a is composed of the desired signal power PDF created by the desired signal power PDF creation unit 13a, the interference signal power CDF created by the interference signal power CDF creation unit 14, and the thermal noise power of the base station. , Calculate the probability that the SINR will be greater than or equal to the required SINR.
  • the probability in base station #i is represented by P- up (i).
  • the communication success rate P'up is calculated by the following equation (4).
  • the communication success rate of the terminal # 6 is subsequently estimated by the same procedure.
  • the processing up to the creation of the interference signal power CDF for each base station, which is carried out by the interference signal power CDF creation unit 14, is common to all terminals (that is, the created interference signal power CDF is the all terminals. Commonly available). Therefore, the communication success rate of the terminal # 6 can be estimated using the same result.
  • the desired signal power PDF creation unit 13a creates the desired signal power PDF from the terminal # 6 to each base station
  • the communication success rate estimation unit 15a calculates P- up (i)
  • the communication success rate estimation unit 15a need only calculate the uplink communication success rate P 'Stay up-by the formula (4).
  • FIG. 3 is a flowchart showing the operation of the wireless terminal accommodation determination device 1a according to the first embodiment of the present invention.
  • the desired terminal (the terminal for which the communication success rate is estimated) is referred to as terminal # j, and the desired number of terminals is referred to as J.
  • the average received power calculation unit 11 divides the area designated by the operator into minute meshes.
  • the average received power calculation unit 11 calculates the average received power of uplink communication from each divided mesh to each base station. Further, the average received power calculation unit 11 calculates the average interference signal power of interference between base stations based on the position of each base station (step S001).
  • the interference signal power CDF creation unit 14 aggregates the number of terminals belonging to each base station, assuming that each terminal belongs to the base station having the maximum average received power (step S002).
  • the interference signal power CDF creation unit 14 creates a PDF of the interference signal power between base stations for each base station based on the traffic pattern of the uplink communication of the terminal to which the base station belongs. Further, the interference signal power CDF creation unit 14 creates a PDF of the instantaneous received power value for each terminal in consideration of the influence of the instantaneous fluctuation and the traffic pattern on the average received power between the base station and the terminal. Then, the interference signal power CDF creation unit 14 creates a PDF of the interference signal power from the terminal for each base station by convolving the created PDF (step S003).
  • a value of 1 is assigned to the variable i that identifies the base station (step S004).
  • the interference signal power CDF creation unit 14 determines the PDF of the interference signal power between the base stations from the interference base station to the desired base station when the base station #i is set as the desired base station, and the PDF of the interference signal power from the terminal to the desired base station. By convolving the interference signal power with the PDF, the PDF of the total interference signal power is created, and the interference signal power CDF is created (step S005).
  • step S006 It is determined whether or not the variable i is equal to the number of base stations I (step S006). If the variable i is not equal to the number of base stations I (step S006-No), the value 1 is added to the variable i (step S007), and the process proceeds to step S005. On the other hand, when the variable i is equal to the number of base stations I (step S006-Yes), the process proceeds to step S008.
  • the value 1 is assigned to the variable j that identifies the desired terminal (step S008).
  • the desired signal power PDF creation unit 13a creates a desired signal power PDF for each base station in consideration of the influence of instantaneous fluctuations on the average received power of the desired signal from the terminal #j to each base station (step S009). ..
  • the communication success rate estimation unit 15a calculates the probability that the SINR becomes equal to or higher than the required SINR for each base station based on the desired signal power PDF, the interference signal power CDF, and the thermal noise power of the base station. Then, the communication success rate estimation unit 15a calculates the uplink communication success rate by the above equation (4) based on the calculated probability (step S010).
  • step S011 It is determined whether or not the variable j is equal to the desired number of terminals J (step S011). If the variable j is not equal to the desired number of terminals J (step S011-No), the value 1 is added to the variable j (step S012), and the process proceeds to step S005. On the other hand, when the variable j is equal to the desired number of terminals J (step S011-Yes), the processing of the flowchart shown in FIG. 3 ends.
  • the wireless terminal accommodation determination device 1a calculates the average received power for each mesh constituting the area for each of the plurality of base stations.
  • the calculation unit 11 and the terminal belonging to each base station are determined based on the average received power, and for each base station, the average received power of the interference signal from the base station is calculated as the PDF (probability density function) of the instantaneous value of the received power.
  • the PDF probability density function
  • a PDF (second probability density function) showing the probability density function of the interference signal power from each terminal is created, and the first probability density function and the second probability density function are convoluted.
  • the interference signal power CDF By substituting the interference signal power CDF creation unit 14 that creates the cumulative distribution function) and the average received power of the desired signal from the desired terminal into the PDF (probability density function) of the received power instantaneous value for each base station.
  • the desired signal power PDF creation unit 13a that creates a PDF (fourth probability density function) showing the probability density function of the desired signal power, the interference signal power CDF (cumulative distribution function), and the desired signal power PDF for each base station.
  • a communication success rate estimation unit 15a for calculating the communication success rate of the desired terminal based on (fourth probability density function) is provided.
  • the wireless terminal accommodation determination device 1a communicates at the planned installation location of the new wireless terminal without performing measurement in the actual environment or Monte Carlo simulation.
  • the success rate can be estimated faster.
  • FIG. 4 is a block diagram showing a functional configuration of the wireless terminal accommodation determination device 1b according to the second embodiment of the present invention.
  • the wireless terminal accommodation determination device 1b includes an average received power calculation unit 11, a radio equipment database 12, a desired signal power PDF creation unit 13b, an interference signal power CDF creation unit 14, and a communication success rate. It includes an estimation unit 15b and a repeat determination unit 16. The point that the repetitive determination unit 16 is provided is different from the functional configuration of the wireless terminal accommodation determination device 1a in the first embodiment.
  • the desired signal power PDF creation unit 13b sorts the base stations in descending order of the average received power based on the average received power Rup (i, j) calculated by the average received power calculation unit 11.
  • the desired signal power PDF creation unit 13b determines how many base stations are to be calculated in the communication success rate estimation in descending order of the average received power, depending on the number of times the processing is currently performed in the iterative processing described later. To do. Then, the desired signal power PDF creation unit 13b creates the desired signal power PDF from the Rup (i, j) of the base station to be calculated.
  • the communication success rate estimation unit 15b includes a desired signal power PDF created by the desired signal power PDF creation unit 13b, an interference signal power CDF created by the interference signal power CDF creation unit 14, and a base station each time in the iterative processing.
  • the probability that the SINR becomes equal to or greater than the required SINR is calculated from the thermal noise power of.
  • the communication success rate estimation unit 15b calculates the communication success rate based on the calculated probability in consideration of the effects of the site diversity and the antenna diversity.
  • the communication success rate estimation unit 15b ends the repetition process when the repetition determination unit 16 determines the end of the repetition process, and outputs the communication success rate calculated at that time as a result.
  • the repetition determination unit 16 determines the amount of change from the communication success rate calculated when the calculated communication success rate becomes equal to or higher than the target value (first predetermined value) in the design work, or in the previous iterative process. Is within a predetermined value (second predetermined value), the end of the repeated process is determined.
  • the target value P target communication success rate in the design work e.g., 0.9, etc.
  • the convergence target value of the change in the communication success rate in the iterative process ⁇ P target e.g 0.02 etc.
  • the average received power calculation unit 11 calculates the average received power Rup (i, 1) for uplink communication from terminal # 1 to base station #i.
  • i is any of 1 to I
  • I represents the number of base stations.
  • the average received power calculation unit 11 calculates the average interference signal power of the interference between the base stations generated between the base stations.
  • the desired signal power PDF creation unit 13b sorts the base stations in descending order of Rup (i, 1).
  • the average received power is calculated for each antenna, and the base stations are sorted by using the highest value among these average received powers as the representative value of the base station.
  • the base stations may be sorted by the average value of these average received powers.
  • the base stations for which the communication success rate is calculated in the nth processing are the top n base stations from the one with the highest average received power.
  • the number of units to be increased for each iterative process is not limited to one, and may be increased by two or exponentially increased from one to two to four.
  • the desired signal power PDF creation unit 13b sets the average received power of the desired signal from the terminal # 1 to the above n upper base stations in the same manner as in the first embodiment described above.
  • the desired signal power PDF is created in consideration of the influence of the instantaneous fluctuation.
  • the communication success rate estimation unit 15b interferes with the desired signal power PDF created by the desired signal power PDF creation unit 13b and the interference signal power CDF creation unit 14 in the same manner as in the first embodiment described above. From the signal power CDF and the thermal noise power of the base station, the communication success rate P'up (n) in the nth process is calculated.
  • the repetition determination unit 16 determines that the communication success rate P'(n) calculated in the nth process satisfies the following equation (5), or the communication success rate P'calculated in the n-1th process. When the amount of change from up (n-1) satisfies the following equation (6), the end of the iterative process is determined. Then, the communication success rate estimation unit 15b outputs the communication success rate P'up (n) calculated at that time to an external device or the like as the estimation result of the communication success rate of the terminal # 1.
  • the repetition determination unit 16 determines that the repetition processing is continued, and shifts to the n + 1th processing. That is, the number of base stations to be calculated is increased by one, the desired signal power PDF creation unit 13b creates a desired signal power PDF for the upper n + 1 base stations, and the communication success rate estimation unit 15b P'up (n + 1). ) Is calculated.
  • the communication success rate of the terminal # 6 is subsequently estimated by the same procedure, but it is carried out by the interference signal power CDF creation unit 14 as in the first embodiment.
  • the processing up to the creation of the interference signal power CDF for each base station is common to all terminals. Therefore, the communication success rate of the terminal # 6 can be estimated using the same result. That is, the desired signal power PDF creation unit 13b creates the desired signal power PDF from the terminal # 6 to each base station, and the communication success rate estimation unit 15b increases the uplink communication success rate P'up (n) by the above equation (4). ) Just be calculated.
  • FIG. 5 is a flowchart showing the operation of the wireless terminal accommodation determination device 1b according to the second embodiment of the present invention.
  • the operator determines the target value P target of the communication success rate in the design work and the convergence target value ⁇ P target of the change in the communication success rate in the iterative processing.
  • the repetition determination unit 16 acquires and sets these determined values (step S101).
  • the average received power calculation unit 11 divides the area designated by the operator into minute meshes. The average received power calculation unit 11 calculates the average received power of uplink communication from each mesh to each base station. Further, the average received power calculation unit 11 calculates the average interference signal power of interference between base stations from the positional relationship of each base station (step S102).
  • the interference signal power CDF creation unit 14 aggregates the number of terminals belonging to each base station, assuming that each terminal belongs to the base station having the maximum average received power (step S103).
  • the interference signal power CDF creation unit 14 creates a PDF of the interference signal power between base stations for each base station based on the traffic pattern of the uplink communication of the terminal to which the base station belongs. Further, the interference signal power CDF creation unit 14 creates a PDF of the instantaneous value of the received power for each terminal in consideration of the influence of the instantaneous fluctuation and the traffic pattern on the average received power between the base station and the terminal. Then, the interference signal power CDF creation unit 14 creates a PDF of the interference signal power from the terminal for each base station by convolving the created PDF (step S104).
  • the interference signal power CDF creation unit 14 sets each base station as a desired base station in order by the same procedure as in steps S004 to S007 shown in FIG. 3, and the interference signal power between the base stations from the interference base station to the desired base station. By convolving the PDF of the above and the PDF of the interference signal power from the terminal to the desired base station, the PDF of the total interference signal power is created, and the interference signal power CDF is created (step S105).
  • the value 1 is assigned to the variable j that identifies the desired terminal (step S106).
  • the desired signal power PDF creation unit 13b sorts the base stations in descending order of the average received power (step S107).
  • the value 1 is assigned to the variable n that identifies the number of times of processing in the iterative processing (step S8).
  • the desired signal power PDF creation unit 13b creates a desired signal power PDF for each base station in consideration of the influence of instantaneous fluctuations on the average received power of the desired signal from the terminal #j to the upper n base stations ( Step S109).
  • the communication success rate estimation unit 15b calculates the probability that the SINR is equal to or higher than the required SINR for each base station from the desired signal power PDF, the interference signal power CDF, and the thermal noise power of the base station. Then, the communication success rate estimation unit 15b calculates the uplink communication success rate by the above equation (4) based on the calculated probability (step S110). Then, the communication success rate estimation unit 15b determines whether or not at least one of the above equations (5) and (6) is satisfied (step S111). If neither is satisfied (step S111-No), the value 1 is added to the variable n (step S112), and the process proceeds to step S109.
  • step S111-Yes when at least one of the above equations (5) and (6) is satisfied (step S111-Yes), the communication success rate estimation unit 15b finishes the calculation of the communication success rate of the terminal #j, and in step S113. Transition.
  • step S113 It is determined whether or not the variable j is equal to the desired number of terminals J (step S113). If the variable j is not equal to the desired number of terminals J (step S113-No), the value 1 is added to the variable j (step S14), and the process proceeds to step S108. On the other hand, when the variable j is equal to the desired number of terminals J (step S113-Yes), the processing of the flowchart shown in FIG. 5 ends.
  • the wireless terminal accommodation determination device 1b according to the second embodiment calculates the communication success rate only for the base stations having a high contribution to the communication success rate (high average received power). .. As a result, the wireless terminal accommodation determination device 1b according to the second embodiment can further shorten the calculation time required for estimating the communication success rate.
  • a part or all of the wireless terminal accommodation determination device 1a and the wireless terminal accommodation determination device 1b in the above-described embodiment may be realized by a computer.
  • the program for realizing this function may be recorded on a computer-readable recording medium, and the program recorded on the recording medium may be read by the computer system and executed.
  • the term "computer system” as used herein includes hardware such as an OS and peripheral devices.
  • the "computer-readable recording medium” refers to a flexible disk, a magneto-optical disk, a portable medium such as a ROM or a CD-ROM, or a storage device such as a hard disk built in a computer system.
  • a "computer-readable recording medium” is a communication line for transmitting a program via a network such as the Internet or a communication line such as a telephone line, and dynamically holds the program for a short period of time. It may also include a program that holds a program for a certain period of time, such as a volatile memory inside a computer system that serves as a server or a client in that case. Further, the above program may be for realizing a part of the above-mentioned functions, and may be further realized for realizing the above-mentioned functions in combination with a program already recorded in the computer system. It may be realized by using hardware such as PLD (Programmable Logic Device) or FPGA (Field Programmable Gate Array).
  • PLD Programmable Logic Device
  • FPGA Field Programmable Gate Array

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Quality & Reliability (AREA)
  • Probability & Statistics with Applications (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

無線端末収容判定装置は、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部と、全干渉信号電力の累積分布関数を作成する干渉信号電力ＣＤＦ作成部と、所望信号電力の確率密度関数を示す第４確率密度関数を作成する所望信号電力ＰＤＦ作成部と、基地局ごとに前記累積分布関数と前記第４確率密度関数とに基づいて前記所望端末の通信成功率を算出する通信成功率推定部とを備える。

Description

無線端末収容判定装置、無線端末収容判定方法、及びプログラム
　本発明は、無線端末収容判定装置、無線端末収容判定方法、及びプログラムに関する。
　ＩｏＴ（Internet of Things；モノのインターネット）技術の発展により、各種センサ及び通信機能を備えたモノ（以下、「ＩｏＴ機器」又は「端末」という。）がネットワークに接続され、遠隔データ収集及びＩｏＴ機器の遠隔制御等を行うことが年々盛んになっている。また、近年、ＩｏＴ機器向けの無線通信方式として、ＬＰＷＡ（Low Power Wide Area network）が注目されている。ＬＰＷＡには、例えば、アンライセンスバンド（無線局免許を必要としない周波数帯）を用いるＬｏＲａＷＡＮ及びＳｉｇｆｏｘ、ライセンスバンド（無線局免許を必要とする周波数帯）を用いるＬＴＥ－Ｍ（Long Term Evolution for Machines）及びＮＢ（Narrow Band）－ＩｏＴ等の無線通信方式がある。
　ＩｏＴ機器は、主に特定の場所に固定して設置され、利用されることが多い。一方、電波伝搬は、場所依存性が高い。これにより、基地局からの電波の受信電力が低い不感地帯にＩｏＴ機器が設置された場合、長期的に通信不可の状態となりうる。そのため、ＩｏＴ機器の設置予定場所において十分な受信電力を確保できるように置局設計を行うことが重要となる。
　また、基地局に収容される端末数が多すぎる場合、あるいは、周辺の基地局等の干渉源からの電波干渉が多い場合、新規にＩｏＴ機器を設置しても十分な通信品質が確保できないことがある。この場合、通信の失敗が頻発することになる。そのため、上述したような設置予定場所における受信電力の評価だけでなく、基地局の無線リソースの使用状況、及び周辺の電波干渉の発生状況を考慮して、新規ＩｏＴ機器を収容可能か否かを事前に判定することが重要となる。そして、当該判定を行うためには、通信成功率を推定する手段が必要となる。
　通信成功率を推定する従来技術として、例えば、特許文献１及び特許文献２に記載の技術がある。特許文献１に記載の技術は、無線機を仮設置して実際に通信を試行することによって通信成功率を測定する技術である。また、特許文献２に記載の技術は、複数の基地局によって構成されるマルチセル環境におけるスループットを、モンテカルロシミュレーションを用いて推定する技術である。
特開２０１４－７２８５３号公報 特許第５０７７３４７号公報
岩井誠人著、「移動通信における電波伝搬　－無線通信シミュレーションのための基礎知識－」、株式会社コロナ社、ｐ．１６－P．１９、２０１２年１１月２２日
　特許文献１に記載の技術は、設置予定場所に端末を仮設置して実環境で通信成功率を測定するため、測定準備及び測定作業に多くの人的コストが必要になるという課題がある。特に、ＬＰＷＡのような広域の無線通信システムでは、実環境での通信成功率の測定に要する人的コストは膨大になる。
　特許文献２に記載の技術は、モンテカルロシミュレーションによりスループットを推定するため、実環境での測定は不要である。しかしながら、特許文献２に記載の技術は、乱数を用いて多数の端末のトラフィックを確率的に発生させる処理を、膨大な回数（例えば、数万回以上）、繰り返し実行することによって、周辺端末との干渉を考慮したスループットの推定を行う。そのため、膨大な計算時間が必要になるという課題がある。
　また、特許文献２では、基地局と端末との間の伝搬損失を、伝搬モデルに基づく距離減衰（長区間変動）とシャドウイング（短区間変動）とを考慮して計算する。しかしながら、特に、端末が固定して設置されることが多いＬＰＷＡのような無線通信システムにおいては、フェージングによる瞬時変動も考慮して通信成功率を推定することが重要になる。フェージングによる瞬時変動を考慮して精度高く通信成功率を推定するためには、例えば乱数を用いて瞬時変動の影響を与えて所望信号及び干渉信号の受信電力瞬時値を計算する処理を、膨大な回数（例えば、数万回以上）、繰り返し実行する必要がある。そのため、更に膨大な計算時間を必要になる。
　本発明は、このような事情を考慮してなされたものであり、実環境での測定やモンテカルロシミュレーションを行うことなく、新規無線端末の設置予定場所における通信成功率をより高速に推定する技術を提供することを目的とする。
　本発明の一態様は、複数の基地局の各々に関して、エリアを構成する微小メッシュごとに平均受信電力を算出する平均受信電力算出部と、前記平均受信電力に基づいて各基地局に所属する端末を決定し、各基地局について、前記基地局からの干渉信号の平均受信電力を受信電力瞬時値の確率密度関数に代入することによって基地局間の干渉信号電力の確率密度関数を示す第１確率密度関数を作成し、各端末からの干渉信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって各端末からの干渉信号電力の確率密度関数を示す第２確率密度関数を作成し、前記第１確率密度関数及び前記第２確率密度関数の畳み込みを行うことによって全干渉信号電力の確率密度関数を示す第３確率密度関数を作成し、前記第３確率密度関数に基づいて全干渉信号電力の累積分布関数を作成する干渉信号電力ＣＤＦ作成部と、各基地局について、所望端末からの所望信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって所望信号電力の確率密度関数を示す第４確率密度関数を作成する所望信号電力ＰＤＦ作成部と、基地局ごとに、前記累積分布関数と前記第４確率密度関数とに基づいて前記所望端末の通信成功率を算出する通信成功率推定部と、を備える無線端末収容判定装置である。
　また、本発明の一態様は、上記の無線端末収容判定装置であって、前記干渉信号電力ＣＤＦ作成部は、前記所望端末を含めた全端末から基地局に到来する信号の受信電力瞬時値を用いて前記第２確率密度関数を作成することによって全端末共通に利用可能な前記累積分布関数を作成する。
　また、本発明の一態様は、上記の無線端末収容判定装置であって、前記平均受信電力算出部は、各基地局の位置関係から基地局間干渉の平均干渉信号電力をさらに算出し、前記干渉信号電力ＣＤＦ作成部は、前記平均受信電力に瞬時変動の影響を考慮して作成された受信電力瞬時値と、端末数と、トラフィックパターンと、前記平均干渉信号電力とに基づいて、前記第１確率密度関数と前記第２確率密度関数とを作成する。
　また、本発明の一態様は、上記の無線端末収容判定装置であって、繰り返し判定部をさらに備え、前記所望信号電力ＰＤＦ作成部は、算出された平均受信電力の高い順に前記基地局をソートし、繰り返し処理ごとに、通信成功率推定において考慮される基地局数を電力の高い順に増加させ、計算対象の基地局について所望信号電力ＰＤＦを作成し、前記通信成功率推定部は、前記累積分布関数と、計算対象の基地局の前記第４確率密度関数と、熱雑音電力とに基づいて通信成功率を算出し、前記繰り返し判定部は、算出された前記通信成功率が第１所定値以上になった場合、又は、１つ前の前記繰り返し処理において算出された前記通信成功率からの値の変化が第２所定値以内になった場合に、前記繰り返し処理を終了する。
　また、本発明の一態様は、複数の基地局の各々に関して、エリアを構成する微小メッシュごとに平均受信電力を算出する平均受信電力算出ステップと、前記平均受信電力に基づいて各基地局に所属する端末を決定し、各基地局について、前記基地局からの干渉信号の平均受信電力を受信電力瞬時値の確率密度関数に代入することによって基地局間の干渉信号電力の確率密度関数を示す第１確率密度関数を作成し、各端末からの干渉信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって各端末からの干渉信号電力の確率密度関数を示す第２確率密度関数を作成し、前記第１確率密度関数及び前記第２確率密度関数の畳み込みを行うことによって全干渉信号電力の確率密度関数を示す第３確率密度関数を作成し、前記第３確率密度関数に基づいて全干渉信号電力の累積分布関数を作成する干渉信号電力ＣＤＦ作成ステップと、各基地局について、所望端末からの所望信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって所望信号電力の確率密度関数を示す第４確率密度関数を作成する所望信号電力ＰＤＦ作成ステップと、基地局ごとに、前記累積分布関数と前記第４確率密度関数とに基づいて前記所望端末の通信成功率を算出する通信成功率推定ステップと、を有する無線端末収容判定方法である。
　また、本発明の一態様は、上記の無線端末収容判定装置としてコンピュータを機能させるためのプログラムである。
　本発明によれば、実環境での測定やモンテカルロシミュレーションを行うことなく、新規無線端末の設置予定場所における通信成功率をより高速に推定することができる。
本発明の第１の実施形態に係る無線端末収容判定装置１ａの機能構成を示すブロック図である。 本発明の第１の実施形態に係る無線端末収容判定装置１ａによる通信成功率推定処理を説明するための図である。 本発明の第１の実施形態に係る無線端末収容判定装置１ａの動作を示すフローチャートである。 本発明の第２の実施形態における無線端末収容判定装置１ｂの機能構成を示すブロック図である。 本発明の第２の実施形態に係る無線端末収容判定装置１ｂの動作を示すフローチャートである。
　以下、図面を参照して、本発明の実施形態について詳しく説明する。なお、以下の実施形態では、ＬｏＲａＷＡＮの、Ｃｌａｓｓ　Ａ及びＣｌａｓｓ　Ｂのプロトコルを用いる無線通信システムを例に挙げて説明する。
　また、以下の実施形態では、端末（ＩｏＴ機器）から基地局への上り通信は、Ｃｌａｓｓ　Ａのプロトコルにより行われる。また、基地局から端末への下り通信は、上記上り通信に対するＡｃｋ（Acknowledge）を送信する場合にはＣｌａｓｓ　Ａのプロトコルにより行われ、アプリケーションサーバ（図示せず）から基地局に転送されてきたデータ送信する場合にはＣｌａｓｓ　Ｂのプロトコルにより行われる。
　また、以下の実施形態では、周辺に複数の基地局がある場合、端末から送信された信号は複数の基地局において受信され、復調される。そして、受信に成功した信号は、上位のネットワークサーバ（図示せず）へ送られる。そして、ネットワークサーバにおいて受信に成功した信号のうち、最も受信電力が高い信号が選択される（サイトダイバーシチ）。また、基地局が複数のアンテナを備える場合、端末から送信された信号は複数のアンテナで受信され、復調される。そして、基地局において受信に成功した信号のうち、最も受信電力が高い信号が選択される（アンテナダイバーシチ）。
　また、以下の実施形態では、下り通信では、直前の上り通信の際に選択された基地局及びアンテナにより下り通信が行われることで、高い通信品質が確保される。
　なお、以下の説明において、通信成功率を推定する対象となる端末を「所望端末」といい、その他の端末を「干渉端末」という。また、所望端末が通信を行う基地局を「所望基地局」といい、その他の基地局を「干渉基地局」という。
＜第１の実施形態＞
　以下、本発明の第１の実施形態について説明する。
［無線端末収容判定装置の機能構成］
　図１は、本発明の第１の実施形態に係る無線端末収容判定装置１ａの機能構成を示すブロック図である。図１に示すように、無線端末収容判定装置１ａは、平均受信電力算出部１１と、無線設備データベース１２と、所望信号電力ＰＤＦ作成部１３ａと、干渉信号電力ＣＤＦ作成部１４と、通信成功率推定部１５ａと、を具備する。
　平均受信電力算出部１１は、一般的な電波伝搬シミュレータに相当するものである。平均受信電力算出部１１は、地図データ（例えば、地形の高さ、建物の高さ、土地の利用分類等のデータ）を外部の機器等からインポートする。平均受信電力算出部１１は、作業者によって指定されたエリアの地図データを、微小なメッシュ（例えば、５ｍ間隔のメッシュ）に分割する。そして、平均受信電力算出部１１は、基地局と各メッシュ（すなわち、各メッシュに相当する位置に設置されているものと仮定した仮想端末）との間の上り通信の平均受信電力を算出する。また、平均受信電力算出部１１は、各基地局の位置に基づいて、基地局間で発生する基地局間干渉の平均干渉信号電力を算出する。
　なお、以下の説明は全て、作業者によって指定されたエリア内に位置する基地局と端末とを判定の対象とする場合についての説明である。
　無線設備データベース１２は、各基地局の、アンテナ本数、及び熱雑音電力等の無線設備に関わるパラメータを保持する。また、無線設備データベース１２は、各端末の、アンテナ本数、熱雑音電力、及びトラフィックパターン（通信頻度、及び１回あたりの通信時間）等の無線設備に関わるパラメータを保持する。なお、ＬｏＲａＷＡＮのように、チャネルが複数あり、使用されるチャネルが通信の度にランダムに選択される場合、端末の通信頻度をチャネル数で割った値を通信頻度とすればよい。
　所望信号電力ＰＤＦ作成部１３ａは、所望端末が位置するメッシュと基地局との間の平均受信電力に、瞬時変動の影響を考慮して、受信電力瞬時値の確率密度関数（ＰＤＦ：Probability Density Function）を基地局ごとに作成する。以下、所望信号の受信電力瞬時値のＰＤＦを、「所望信号電力ＰＤＦ」という。
　干渉信号電力ＣＤＦ作成部１４は、干渉端末が位置するメッシュと基地局との間の平均受信電力に、瞬時変動の影響とトラフィックパターンとを考慮して、受信電力瞬時値の確率密度関数（ＰＤＦ：Probability Density Function）を作成する。また、干渉信号電力ＣＤＦ作成部１４は、所望基地局と干渉基地局との間の平均受信電力と干渉基地局配下の干渉端末のトラフィックパターンとに基づいて、干渉信号電力のＰＤＦを作成する。干渉信号電力ＣＤＦ作成部１４は、作成したＰＤＦを集計してＣＤＦを作成する。以下、干渉信号のＰＤＦを集計することによって作成されたＣＤＦを、「干渉信号電力ＣＤＦ」という。
　通信成功率推定部１５ａは、所望信号電力ＰＤＦ作成部１３ａで作成された所望信号電力ＰＤＦと、干渉信号電力ＣＤＦ作成部１４で作成された干渉信号電力ＣＤＦと、基地局の熱雑音電力とから、所望信号電力対干渉及び雑音電力比（ＳＩＮＲ；Signal-to-Interference plus Noise power Ratio）が所要ＳＩＮＲ以上となる確率を算出する。そして、通信成功率推定部１５ａは、算出された確率に基づき、サイトダイバーシチ及びアンテナダイバーシチの効果を考慮して、通信成功率を算出する。ここで、所要ＳＩＮＲとは、用いられる変調方式の無線通信において、通信が成功するために必要となるＳＩＮＲの値を示す。
［通信成功率推定］
　以下、通信成功率の推定処理について説明する。
　図２は、本発明の第１の実施形態に係る無線端末収容判定装置１ａによる通信成功率推定処理を説明するための図である。なお、以下では、ＬｏＲａＷＡＮの複数のチャネルのうち、Ｃｌａｓｓ　Ａの通信のみで使用されるチャネルにおける通信成功率を推定する場合について説明する。なお、Ｃｌａｓｓ　Ａ及びＣｌａｓｓ　Ｂの通信で使用されるチャネルにおける通信成功率を推定する場合には、干渉信号電力ＣＤＦ作成部１４で干渉基地局からの干渉信号電力のＰＤＦを作成する際に、Ｃｌａｓｓ　Ｂの通信発生時間も考慮すればよい。
　以下、図２に示すように、基地局３台で面的に無線エリアが構築されている環境に端末＃１及び端末＃６を新規に設置する際に、端末＃１及び端末＃６の通信成功率を推定する場合について説明する。各基地局は、それぞれ１本のアンテナを備える。基地局＃ｉ（ｉ＝１，２，３）における、端末＃ｊ（ｊ＝１，２，・・・，１５）からの上り通信の平均受信電力をＲｕｐ（ｉ，ｊ）と表す。
　上述したように、平均受信電力算出部１１は、一般的な電波伝搬シミュレータに相当するものである。平均受信電力算出部１１は、作業者によって指定されたエリアを、微小なメッシュに分割する。そして、平均受信電力算出部１１は、基地局と各メッシュ（すなわち、各メッシュに相当する位置に設置されているものと仮定した仮想端末）との間の上り通信の平均受信電力を算出する。平均受信電力算出部１１は、端末＃ｊが位置するメッシュにおける平均受信電力の値を取得し、基地局＃ｉと端末＃ｊとの間の上り通信の平均受信電力Ｒｕｐ（ｉ，ｊ）とする。また、平均受信電力算出部１１は、各基地局の位置に基づて、基地局間で発生する基地局間干渉の平均干渉信号電力を算出する。基地局＃ｉから基地局＃ｋへの基地局間干渉の平均干渉信号電力を、ＲＩＣＩ（ｉ，ｋ）と表す。
　干渉信号電力ＣＤＦ作成部１４は、まず、端末からの干渉信号電力のＰＤＦ（第２確率密度関数）と基地局間の干渉信号電力のＰＤＦ（第１確率密度関数）とを作成する。
　干渉信号電力ＣＤＦ作成部１４は、所望端末の所望基地局と端末との間の平均受信電力に、瞬時変動の影響とトラフィックパターンを考慮して、受信電力瞬時値のＰＤＦを端末ごとに作成し、作成された複数のＰＤＦの畳み込みを行うことによって、端末からの干渉信号電力のＰＤＦ（所望端末の所望基地局における全端末からの干渉信号のＰＤＦ）を得る。
　なお、受信電力瞬時値を算出する単純な方法として、モンテカルロシミュレーション法がある。モンテカルロシミュレーション法では、基地局＃ｉに複数の経路で到来する電波の位相を乱数を用いてランダムに決定して合成することで、受信電力瞬時値を算出する。そして、これを繰り返し行うことによって、受信電力瞬時値のＰＤＦを作成することができる。しかしながら、この繰り返し計算は、一般的に膨大な回数（数万回程度）行われる。そのため、膨大な計算時間が必要になるという課題がある。
　そこで、本実施形態においては、フェージングによる瞬時変動を、下記の解析式で表す。これにより、本実施形態では、モンテカルロシミュレーションを行うことなく、所望信号電力ＰＤＦを得る。
　受信信号の振幅をａとした場合、レイリーフェージングにより瞬時変動する振幅ａのＰＤＦは、Ｒｕｐ（ｉ，ｊ）を用いて、以下の式（１）で表すことができる（非特許文献１参照）。
Figure JPOXMLDOC01-appb-M000001
　これを電力ｒ＝ａ２のＰＤＦに変換することで、受信電力瞬時値のＰＤＦは、以下の式（２）で表すことができる。
Figure JPOXMLDOC01-appb-M000002
　一方、基地局間の干渉信号電力のＰＤＦを作成するためには、各端末がどの基地局に所属するか（すなわち、端末への下り通信がどの基地局から送信されるか）を決定し、各基地局に所属する端末の端末数を集計する必要がある。ここでは、各端末は、Ｒｕｐ（ｉ，ｊ）が最大となる基地局＃ｉに所属するものとして、端末数が集計される。そして、基地局ごとに、所属する端末の上り通信のトラフィックパターンに基づいて、端末へのＡｃｋ（Ｃｌａｓｓ　Ａによる下り通信）の通信発生時間が算出され、全体の時間のうち前記算出した通信発生時間の間、干渉電力ＲＩＣＩ（ｉ，ｋ）の基地局間干渉が発生するものとしてＰＤＦが作成される。
　なお、基地局は、一般的に高層ビル屋上の等の高い位置に設置されることが多く、基地局同士は、周辺に反射物や遮蔽物がなく見通しがある環境であることが一般的である。そのため、基地局間干渉については、基地局同士の位置関係に基づいて算出される平均干渉信号電力を用いるものとしている。但し、基地局間の見通し経路を到来する直接波の他に、電力の低い反射波が到来するような場合においては、ライスフェージングを仮定し、基地局同士の位置関係から算出される平均干渉信号電力に、ライスフェージングによる瞬時変動の影響を考慮することによって、基地局間の干渉信号電力瞬時値のＰＤＦを作成すればよい。
　そして、干渉信号電力ＣＤＦ作成部１４は、基地局＃１が所望基地局である場合の干渉信号電力ＣＤＦを作成するにあたっては、干渉基地局＃２と干渉基地局＃３とから基地局＃１への干渉信号電力のＰＤＦ２つと、基地局＃１への端末からの干渉信号電力のＰＤＦとの畳み込みを行う。これにより、干渉信号電力ＣＤＦ作成部１４は、全干渉信号電力のＰＤＦ（ＰＤＦＩ（ｒ）と定義する）（第３確率密度関数）を作成し、以下の式（３）によって干渉信号電力ＣＤＦ（ＣＤＦＩ＿１（ｙ））（累積分布関数）を作成する。
Figure JPOXMLDOC01-appb-M000003
　同様にして、干渉信号電力ＣＤＦ作成部１４は、基地局＃２が所望基地局である場合の干渉信号電力ＣＤＦを作成するにあたっては、干渉基地局＃１と干渉基地局＃３とから基地局＃２への干渉信号電力のＰＤＦ２つと、基地局＃２への端末からの干渉信号電力のＰＤＦとの畳み込みを行う。これにより、干渉信号電力ＣＤＦ作成部１４は、全干渉信号電力のＰＤＦ（第３確率密度関数）を作成し、干渉信号電力ＣＤＦ（ＣＤＦＩ＿２（ｙ））（累積分布関数）を作成する。
　同様にして、干渉信号電力ＣＤＦ作成部１４は、基地局＃３が所望基地局である場合の干渉信号電力ＣＤＦを作成するにあたっては、干渉基地局＃１と干渉基地局＃２とから基地局＃３への干渉信号電力のＰＤＦ２つと、基地局＃３への端末からの干渉信号電力のＰＤＦとの畳み込みを行う。これにより、干渉信号電力ＣＤＦ作成部１４は、全干渉信号電力のＰＤＦ（第３確率密度関数）を作成し、干渉信号電力ＣＤＦ（ＣＤＦＩ＿３（ｙ））（累積分布関数）を作成する。
　所望信号電力ＰＤＦ作成部１３ａは、端末＃１から基地局＃１、基地局＃２、及び基地局＃３への所望信号の平均受信電力Ｒｕｐ（１，１）、Ｒｕｐ（２，１）、及びＲｕｐ（３，１）に、瞬時変動の影響を考慮して、所望信号電力ＰＤＦ（ＰＤＦＳ＿１（ｙ）、ＰＤＦＳ＿２（ｙ）、及びＰＤＦＳ＿３（ｙ））を作成する。
　通信成功率推定部１５ａは、所望信号電力ＰＤＦ作成部１３ａで作成された所望信号電力ＰＤＦと、干渉信号電力ＣＤＦ作成部１４で作成された干渉信号電力ＣＤＦと、基地局の熱雑音電力とから、ＳＩＮＲが所要ＳＩＮＲ以上となる確率を算出する。ここで、基地局＃ｉにおける確率をＰｕｐ（ｉ）で表す。そして、サイトダイバーシチ効果を考慮すると、通信成功率Ｐ’ｕｐは以下の式（４）によって算出される。
Figure JPOXMLDOC01-appb-M000004
　上記により端末＃１の通信成功率の推定が完了すると、続いて同手順で端末＃６の通信成功率を推定する。ここで、干渉信号電力ＣＤＦ作成部１４で実施される、基地局ごとの干渉信号電力ＣＤＦの作成までの処理は、全端末で共通である（すなわち、作成された干渉信号電力ＣＤＦは、全端末共通に利用可能である）。そのため、同じ結果を用いて端末＃６の通信成功率推定を行うことができる。すなわち、所望信号電力ＰＤＦ作成部１３ａが端末＃６から各基地局への所望信号電力ＰＤＦを作成し、通信成功率推定部１５ａがＰｕｐ（ｉ）を算出し、通信成功率推定部１５ａが、上記式（４）により上り通信成功率Ｐ’ｕｐを算出するだけでよい。
［無線端末収容判定装置の動作］
　以下、通信成功率推定処理における無線端末収容判定装置１ａの動作について説明する。
　図３は、本発明の第１の実施形態に係る無線端末収容判定装置１ａの動作を示すフローチャートである。なお、以下、所望端末（通信成功率を推定する対象の端末）を端末＃ｊとし、所望端末数をＪとする。
　平均受信電力算出部１１は、作業者によって指定されたエリアを微小なメッシュに分割する。平均受信電力算出部１１は、分割された各メッシュから各基地局への上り通信の平均受信電力を算出する。また、平均受信電力算出部１１は、各基地局の位置に基づいて基地局間干渉の平均干渉信号電力を算出する（ステップＳ００１）。
　干渉信号電力ＣＤＦ作成部１４は、平均受信電力が最大となる基地局に各端末が所属するものとして、各基地局の所属端末数を集計する（ステップＳ００２）。
　次に、干渉信号電力ＣＤＦ作成部１４は、所属端末の上り通信のトラフィックパターンに基づいて、基地局間の干渉信号電力のＰＤＦを基地局ごとに作成する。また、干渉信号電力ＣＤＦ作成部１４は、基地局と端末との間の平均受信電力に、瞬時変動の影響とトラフィックパターンとを考慮して、受信電力瞬時値のＰＤＦを端末ごとに作成する。そして、干渉信号電力ＣＤＦ作成部１４は、作成されたＰＤＦの畳み込みを行うことで、端末からの干渉信号電力のＰＤＦを基地局ごとに作成する（ステップＳ００３）。
　基地局を識別する変数ｉに値１が代入される（ステップＳ００４）。
　干渉信号電力ＣＤＦ作成部１４は、基地局＃ｉを所望基地局としたときの、干渉基地局から所望基地局への基地局間の干渉信号電力のＰＤＦと、所望基地局への端末からの干渉信号電力のＰＤＦとの畳み込みを行うことで全干渉信号電力のＰＤＦを作成し、干渉信号電力ＣＤＦを作成する（ステップＳ００５）。
　変数ｉが基地局数Ｉと等しいか否かが判定される（ステップＳ００６）。変数ｉが基地局数Ｉと等しくない場合（ステップＳ００６－Ｎｏ）、変数ｉに値１が加算され（ステップＳ００７）、ステップＳ００５に移行する。一方、変数ｉが基地局数Ｉと等しい場合（ステップＳ００６－Ｙｅｓ）、ステップＳ００８に移行する。
　所望端末を識別する変数ｊに値１が代入される（ステップＳ００８）。
　所望信号電力ＰＤＦ作成部１３ａは、端末＃ｊから各基地局への所望信号の平均受信電力に、瞬時変動の影響を考慮して、所望信号電力ＰＤＦを基地局ごとに作成する（ステップＳ００９）。
　通信成功率推定部１５ａは、所望信号電力ＰＤＦと、干渉信号電力ＣＤＦと、基地局の熱雑音電力とに基づいて、ＳＩＮＲが所要ＳＩＮＲ以上となる確率を基地局ごとに算出する。そして、通信成功率推定部１５ａは、算出された確率に基づき、上記式（４）によって上り通信成功率を算出する（ステップＳ０１０）。
　変数ｊが所望端末数Ｊと等しいか否かが判定される（ステップＳ０１１）。変数ｊが所望端末数Ｊと等しくない場合（ステップＳ０１１－Ｎｏ）、変数ｊに値１が加算され（ステップＳ０１２）、ステップＳ００５に移行する。一方、変数ｊが所望端末数Ｊと等しい場合（ステップＳ０１１－Ｙｅｓ）、図３に示すフローチャートの処理が終了する。
　なお、以上の説明では、Ｃｌａｓｓ　Ａの通信で用いられるチャネルが複数あり、端末からのＣｌａｓｓ　Ａの上り通信に対するＡｃｋ送信（Ｃｌａｓｓ　Ａによる下り通信）が前記複数チャネルのいずれかのチャネルで行われることを前提に説明してきた。但し、Ａｃｋ送信用のチャネルが用意され、Ｃｌａｓｓ　Ａの上り通信と下り通信とが干渉しない場合には、干渉信号電力ＣＤＦを作成するにあたって基地局間の干渉信号電力のＰＤＦを作成する必要はなく、端末からの干渉信号電力のＰＤＦのみを用いればよい。
　以上説明したように、本発明の第１の実施の実施形態に係る無線端末収容判定装置１ａは、複数の基地局の各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部１１と、前記平均受信電力に基づいて各基地局に所属する端末を決定し、各基地局について、前記基地局からの干渉信号の平均受信電力を受信電力瞬時値のＰＤＦ（確率密度関数）に代入することによって基地局間の干渉信号電力の確率密度関数を示すＰＤＦ（第１確率密度関数）を作成し、各端末からの干渉信号の平均受信電力を前記受信電力瞬時値の前記ＰＤＦ（確率密度関数）に代入することによって各端末からの干渉信号電力の確率密度関数を示すＰＤＦ（第２確率密度関数）を作成し、前記第１確率密度関数及び前記第２確率密度関数の畳み込みを行うことによって全干渉信号電力の確率密度関数を示すＰＤＦ（第３確率密度関数）を作成し、前記第３確率密度関数に基づいて全干渉信号電力の累積分布関数を示す干渉信号電力ＣＤＦ（累積分布関数）を作成する干渉信号電力ＣＤＦ作成部１４と、各基地局について、所望端末からの所望信号の平均受信電力を前記受信電力瞬時値の前記ＰＤＦ（確率密度関数）に代入することによって所望信号電力の確率密度関数を示すＰＤＦ（第４確率密度関数）を作成する所望信号電力ＰＤＦ作成部１３ａと、基地局ごとに、前記干渉信号電力ＣＤＦ（累積分布関数）と前記所望信号電力ＰＤＦ（第４確率密度関数）とに基づいて前記所望端末の通信成功率を算出する通信成功率推定部１５ａと、を備える。
　上記の構成を備えることにより、本発明の第１の実施の実施形態に係る無線端末収容判定装置１ａは、実環境での測定やモンテカルロシミュレーションを行うことなく、新規無線端末の設置予定場所における通信成功率をより高速に推定することができる。
＜第２の実施形態＞
　第１の実施形態では、全基地局によるサイトダイバーシチ効果を考慮して通信成功率を推定する場合について説明した。一方、以下に説明する第２の実施形態では、通信成功率への寄与度の高い基地局のみを計算対象として通信成功率を推定する場合について説明する。なお、以下の説明では、第１の実施形態との相違点についてのみ説明する。
［無線端末収容判定装置の機能構成］
　図４は、本発明の第２の実施形態における無線端末収容判定装置１ｂの機能構成を示すブロック図である。図４に示すように、無線端末収容判定装置１ｂは、平均受信電力算出部１１と、無線設備データベース１２と、所望信号電力ＰＤＦ作成部１３ｂと、干渉信号電力ＣＤＦ作成部１４と、通信成功率推定部１５ｂと、繰り返し判定部１６と、を具備する。繰り返し判定部１６を具備する点が、第１の実施形態における無線端末収容判定装置１ａの機能構成との相違点である。
　所望信号電力ＰＤＦ作成部１３ｂは、平均受信電力算出部１１で算出された平均受信電力Ｒｕｐ（ｉ，ｊ）に基づいて、平均受信電力の高い順に基地局をソートする。所望信号電力ＰＤＦ作成部１３ｂは、後述する繰り返し処理において、現在何回目の処理なのかに応じて、平均受信電力の高い順に何台の基地局を通信成功率推定における計算対象とするかを決定する。そして、所望信号電力ＰＤＦ作成部１３ｂは、計算対象の基地局のＲｕｐ（ｉ，ｊ）から所望信号電力ＰＤＦを作成する。
　通信成功率推定部１５ｂは、繰り返し処理における各回において、所望信号電力ＰＤＦ作成部１３ｂで作成された所望信号電力ＰＤＦと、干渉信号電力ＣＤＦ作成部１４で作成された干渉信号電力ＣＤＦと、基地局の熱雑音電力とから、ＳＩＮＲが所要ＳＩＮＲ以上となる確率を算出する。そして、通信成功率推定部１５ｂは、算出された確率に基づいて、サイトダイバーシチやアンテナダイバーシチの効果を考慮して、通信成功率を算出する。そして、通信成功率推定部１５ｂは、繰り返し判定部１６にて繰り返し処理終了の判定がなされると繰り返し処理を終了し、その時点で算出された通信成功率を結果として出力する。
　繰り返し判定部１６は、算出された通信成功率が当該設計作業における目標値（第１所定値）以上となった場合、又は、１つ前の繰り返し処理で算出された通信成功率からの変化量が所定値（第２所定値）以内になった場合に、繰り返し処理終了の判定を行う。
［通信成功率推定］
　以下、通信成功率の推定処理について説明する。第１の実施形態と同様に、端末＃１と端末＃６とを新規に設置する際に、端末＃１と端末＃６の通信成功率を推定する場合について説明する。
　まず端末＃１の通信成功率の推定において、作業者により、当該設計作業における通信成功率の目標値Ｐｔａｒｇｅｔ（例えば、０．９等）と、繰り返し処理における通信成功率の変化の収束目標値ΔＰｔａｒｇｅｔ（例えば、０．０２等）とが決定される。そして、繰り返し判定部１６に対してそれらの値が設定される。
　平均受信電力算出部１１は、基地局＃ｉへの端末＃１からの上り通信の平均受信電力Ｒｕｐ（ｉ，１）を算出する。なお、ｉは、１～Ｉのいずれかであり、Ｉは基地局数を表す。また、第１の実施形態を同様に、平均受信電力算出部１１は、基地局間で発生する基地局間干渉の平均干渉信号電力を算出する。
　無線設備データベース１２、及び干渉信号電力ＣＤＦ作成部１４の構成は、第１の実施形態と同様であるため、説明を省略する。
　所望信号電力ＰＤＦ作成部１３ｂは、Ｒｕｐ（ｉ，１）の高い順に基地局をソートする。なお、基地局が複数のアンテナを備える場合、各アンテナについてそれぞれ平均受信電力が算出されるが、これらの平均受信電力のうち一番高い値を当該基地局の代表値として基地局をソートしてもよいし、これらの平均受信電力を平均した値を用いて基地局をソートしてもよい。
　以下、繰り返し処理ごとに、通信成功率の計算対象とする基地局を１台ずつ増加させる場合について説明する。すなわち、ｎ回目の処理において通信成功率の計算対象となる基地局は、平均受信電力が高いほうから上位ｎ台の基地局である。
　なお、繰り返し処理ごとに増やす台数は１台に限られるものはなく、２台ずつ増やしたり、１台→２台→４台と指数関数的に増やしたりしてもよい。所望信号電力ＰＤＦ作成部１３ｂは、ｎ回目の処理において、上述した第１の実施形態と同様の方法で、端末＃１から上記の上位ｎ台の基地局への所望信号の平均受信電力に、瞬時変動の影響を考慮して、所望信号電力ＰＤＦを作成する。
　通信成功率推定部１５ｂは、上述した第１の実施形態と同様の方法で、所望信号電力ＰＤＦ作成部１３ｂで作成された所望信号電力ＰＤＦと、干渉信号電力ＣＤＦ作成部１４で作成された干渉信号電力ＣＤＦと、基地局の熱雑音電力とから、ｎ回目の処理における通信成功率Ｐ’ｕｐ（ｎ）を算出する。
　繰り返し判定部１６は、ｎ回目の処理で算出された通信成功率Ｐ’（ｎ）が以下の式（５）を満たす場合、又は、ｎ－１回目の処理で算出された通信成功率Ｐ’ｕｐ（ｎ－１）からの変化量が以下の式（６）を満たす場合に、繰り返し処理終了の判定を行う。そして、通信成功率推定部１５ｂは、その回において算出された通信成功率Ｐ’ｕｐ（ｎ）を、端末＃１の通信成功率の推定結果として外部の装置等へ出力する。
Figure JPOXMLDOC01-appb-M000005
Figure JPOXMLDOC01-appb-M000006
　上記条件を満たさない場合、繰り返し判定部１６は、繰り返し処理継続の判定を行い、ｎ＋１回目の処理に移行する。すなわち、計算対象の基地局が１台増やされ、所望信号電力ＰＤＦ作成部１３ｂは、上位ｎ＋１台の基地局について所望信号電力ＰＤＦを作成し、通信成功率推定部１５ｂは、Ｐ’ｕｐ（ｎ＋１）の値を算出する。
　上記により端末＃１の通信成功率の推定が完了すると、続いて同様の手順で端末＃６の通信成功率を推定するが、第１の実施形態と同様、干渉信号電力ＣＤＦ作成部１４で実施される、基地局ごとの干渉信号電力ＣＤＦの作成までの処理は、全端末で共通である。そのため、同じ結果を用いて端末＃６の通信成功率推定を行うことができる。すなわち、所望信号電力ＰＤＦ作成部１３ｂが端末＃６から各基地局への所望信号電力ＰＤＦを作成し、通信成功率推定部１５ｂが、上記式（４）により上り通信成功率Ｐ’ｕｐ（ｎ）を算出するだけでよい。
［無線端末収容判定装置の動作］
　以下、通信成功率推定処理における無線端末収容判定装置１ｂの動作について説明する。
　図５は、本発明の第２の実施形態に係る無線端末収容判定装置１ｂの動作を示すフローチャートである。
　まず、作業者により、当該設計作業における通信成功率の目標値Ｐｔａｒｇｅｔと、繰り返し処理における通信成功率の変化の収束目標値ΔＰｔａｒｇｅｔとが決定される。繰り返し判定部１６は、決定されたこれらの値を取得し設定する（ステップＳ１０１）。
　平均受信電力算出部１１は、作業者によって指定されたエリアを微小なメッシュに分割する。平均受信電力算出部１１は、各メッシュから各基地局への上り通信の平均受信電力を算出する。また、平均受信電力算出部１１は、各基地局の位置関係から基地局間干渉の平均干渉信号電力を算出する（ステップＳ１０２）。
　干渉信号電力ＣＤＦ作成部１４は、各端末が平均受信電力最大の基地局に所属するものとして、各基地局の所属端末数を集計する（ステップＳ１０３）。
　次に、干渉信号電力ＣＤＦ作成部１４は、所属端末の上り通信のトラフィックパターンに基づいて、基地局間の干渉信号電力のＰＤＦを基地局ごとに作成する。また、干渉信号電力ＣＤＦ作成部１４は、基地局と端末との間の平均受信電力に、瞬時変動の影響とトラフィックパターンを考慮して、受信電力瞬時値のＰＤＦを端末ごとに作成する。そして、干渉信号電力ＣＤＦ作成部１４は、作成されたＰＤＦの畳み込みを行うことで、端末からの干渉信号電力のＰＤＦを基地局ごとに作成する（ステップＳ１０４）。
　干渉信号電力ＣＤＦ作成部１４は、図３に示したステップＳ００４～Ｓ００７と同様の手順により、各基地局を順に所望基地局として、干渉基地局から所望基地局への基地局間の干渉信号電力のＰＤＦと、所望基地局への端末からの干渉信号電力のＰＤＦとの畳み込みを行うことで、全干渉信号電力のＰＤＦを作成し、干渉信号電力ＣＤＦを作成する（ステップＳ１０５）。
　所望端末を識別する変数ｊに値１が代入される（ステップＳ１０６）。
　所望信号電力ＰＤＦ作成部１３ｂは、平均受信電力の高い順に基地局をソートする（ステップＳ１０７）。
　繰り返し処理における何回目の処理かを識別する変数ｎに値１が代入される（ステップＳ８）。所望信号電力ＰＤＦ作成部１３ｂは、端末＃ｊから上位ｎ台の基地局への所望信号の平均受信電力に、瞬時変動の影響を考慮して、所望信号電力ＰＤＦを基地局ごとに作成する（ステップＳ１０９）。
　通信成功率推定部１５ｂは、所望信号電力ＰＤＦと、干渉信号電力ＣＤＦと、基地局の熱雑音電力とから、ＳＩＮＲが所要ＳＩＮＲ以上となる確率を基地局ごとに算出する。そして、通信成功率推定部１５ｂは、算出された確率に基づいて、上記式（４）により上り通信成功率を算出する（ステップＳ１１０）。そして、通信成功率推定部１５ｂは、上記式（５）及び式（６）のうち少なくとも一方を満たすか否かを判定する（ステップＳ１１１）。どちらも満たさない場合（ステップＳ１１１－Ｎｏ）、変数ｎに値１が加算され（ステップＳ１１２）、ステップＳ１０９に移行する。
　一方、上記式（５）及び式（６）のうち少なくとも一方を満たす場合（ステップＳ１１１－Ｙｅｓ）、通信成功率推定部１５ｂは、端末＃ｊの通信成功率の算出を終了し、ステップＳ１１３に移行する。
　変数ｊが所望端末数Ｊと等しいか否かが判定される（ステップＳ１１３）。変数ｊが所望端末数Ｊと等しくない場合（ステップＳ１１３－Ｎｏ）、変数ｊに値１が加算され（ステップＳ１４）、ステップＳ１０８に移行する。一方、変数ｊが所望端末数Ｊと等しい場合（ステップＳ１１３－Ｙｅｓ）、図５に示すフローチャートの処理が終了する。
　以上説明したように、第２の実施形態に係る無線端末収容判定装置１ｂは、通信成功率への寄与度の高い（平均受信電力の高い）基地局のみを計算対象として通信成功率を算出する。これにより、第２の実施形態に係る無線端末収容判定装置１ｂは、通信成功率の推定に要する計算時間を、より短縮させることができる。
　上述した実施形態における無線端末収容判定装置１ａ及び無線端末収容判定装置１ｂの一部又は全部を、コンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。更に「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、上述した機能の一部を実現するためのものであっても良く、更に上述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＰＬＤ（Programmable Logic Device）やＦＰＧＡ（Field Programmable Gate Array）等のハードウェアを用いて実現されるものであってもよい。
　以上、図面を参照して本発明の実施形態を説明してきたが、上記実施形態は本発明の例示に過ぎず、本発明が上記実施形態に限定されるものではないことは明らかである。したがって、本発明の技術思想及び要旨を逸脱しない範囲で構成要素の追加、省略、置換、及びその他の変更を行ってもよい。
１ａ、１ｂ…無線端末収容判定装置、１１…平均受信電力算出部、１２…無線設備データベース、１３ａ、１３ｂ…所望信号電力ＰＤＦ作成部、１４…干渉信号電力ＣＤＦ作成部、１５ａ、１５ｂ…通信成功率推定部、１６…判定部

Claims (6)

  1. 　複数の基地局の各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部と、
    　前記平均受信電力に基づいて各基地局に所属する端末を決定し、各基地局について、前記基地局からの干渉信号の平均受信電力を受信電力瞬時値の確率密度関数に代入することによって基地局間の干渉信号電力の確率密度関数を示す第１確率密度関数を作成し、各端末からの干渉信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって各端末からの干渉信号電力の確率密度関数を示す第２確率密度関数を作成し、前記第１確率密度関数及び前記第２確率密度関数の畳み込みを行うことによって全干渉信号電力の確率密度関数を示す第３確率密度関数を作成し、前記第３確率密度関数に基づいて全干渉信号電力の累積分布関数を作成する干渉信号電力ＣＤＦ作成部と、
    　各基地局について、所望端末からの所望信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって所望信号電力の確率密度関数を示す第４確率密度関数を作成する所望信号電力ＰＤＦ作成部と、
    　基地局ごとに、前記累積分布関数と前記第４確率密度関数とに基づいて前記所望端末の通信成功率を算出する通信成功率推定部と、
    　を備える無線端末収容判定装置。
  2. 　前記干渉信号電力ＣＤＦ作成部は、
    　前記所望端末を含めた全端末から基地局に到来する信号の受信電力瞬時値を用いて前記第２確率密度関数を作成することによって全端末共通に利用可能な前記累積分布関数を作成する
    　請求項１に記載の無線端末収容判定装置。
  3. 　前記平均受信電力算出部は、
    　各基地局の位置関係から基地局間干渉の平均干渉信号電力をさらに算出し、
    　前記干渉信号電力ＣＤＦ作成部は、
    　前記平均受信電力に瞬時変動の影響を考慮して作成された受信電力瞬時値と、端末数と、トラフィックパターンと、前記平均干渉信号電力とに基づいて、前記第１確率密度関数と前記第２確率密度関数とを作成する
    　請求項１または請求項２に記載の無線端末収容判定装置。
  4. 　繰り返し判定部
    　をさらに備え、
    　前記所望信号電力ＰＤＦ作成部は、
    　算出された平均受信電力の高い順に前記基地局をソートし、繰り返し処理ごとに、通信成功率推定で考慮される基地局数を電力の高い順に増加させ、計算対象の基地局について所望信号電力ＰＤＦを作成し、
    　前記通信成功率推定部は、
    　前記累積分布関数と、計算対象の基地局の前記第４確率密度関数と、熱雑音電力とに基づいて通信成功率を算出し、
    　前記繰り返し判定部は、
    　算出された前記通信成功率が第１所定値以上になった場合、又は、１つ前の前記繰り返し処理において算出された前記通信成功率からの値の変化が第２所定値以内になった場合に、前記繰り返し処理を終了する
    　請求項１から請求項３のうちいずれか一項に記載の無線端末収容判定装置。
  5. 　複数の基地局の各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出ステップと、
    　前記平均受信電力に基づいて各基地局に所属する端末を決定し、各基地局について、前記基地局からの干渉信号の平均受信電力を受信電力瞬時値の確率密度関数に代入することによって基地局間の干渉信号電力の確率密度関数を示す第１確率密度関数を作成し、各端末からの干渉信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって各端末からの干渉信号電力の確率密度関数を示す第２確率密度関数を作成し、前記第１確率密度関数及び前記第２確率密度関数の畳み込みを行うことによって全干渉信号電力の確率密度関数を示す第３確率密度関数を作成し、前記第３確率密度関数に基づいて全干渉信号電力の累積分布関数を作成する干渉信号電力ＣＤＦ作成ステップと、
    　各基地局について、所望端末からの所望信号の平均受信電力を前記受信電力瞬時値の前記確率密度関数に代入することによって所望信号電力の確率密度関数を示す第４確率密度関数を作成する所望信号電力ＰＤＦ作成ステップと、
    　基地局ごとに、前記累積分布関数と前記第４確率密度関数とに基づいて前記所望端末の通信成功率を算出する通信成功率推定ステップと、
    　を有する無線端末収容判定方法。
  6. 　請求項１から請求項４のうちいずれか一項に記載の無線端末収容判定装置としてコンピュータを機能させるためのプログラム。
PCT/JP2020/008607 2019-03-08 2020-03-02 無線端末収容判定装置、無線端末収容判定方法、及びプログラム WO2020184246A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US17/435,093 US11785475B2 (en) 2019-03-08 2020-03-02 Wireless terminal accommodation determination apparatus, wireless terminal accommodation determination method and program

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019-042476 2019-03-08
JP2019042476A JP7096499B2 (ja) 2019-03-08 2019-03-08 無線端末収容判定装置、無線端末収容判定方法、及びプログラム

Publications (1)

Publication Number Publication Date
WO2020184246A1 true WO2020184246A1 (ja) 2020-09-17

Family

ID=72354609

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/008607 WO2020184246A1 (ja) 2019-03-08 2020-03-02 無線端末収容判定装置、無線端末収容判定方法、及びプログラム

Country Status (3)

Country Link
US (1) US11785475B2 (ja)
JP (1) JP7096499B2 (ja)
WO (1) WO2020184246A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001320761A (ja) * 2000-03-01 2001-11-16 Yrp Mobile Telecommunications Key Tech Res Lab Co Ltd チャネル割当方法および装置
JP2004215075A (ja) * 2003-01-07 2004-07-29 Ntt Docomo Inc Ｃｄｍａ適応アンテナアレイ基地局設計方法および装置、ｃｄｍａ適応アンテナアレイ基地局設計プログラムおよびそのプログラムを記録した記録媒体
JP2012244563A (ja) * 2011-05-24 2012-12-10 Nippon Telegr & Teleph Corp <Ntt> 無線通信方法
JP2015050502A (ja) * 2013-08-30 2015-03-16 株式会社国際電気通信基礎技術研究所 無線通信装置、収容予測装置、境界情報取得装置、収容予測方法、境界情報取得方法、及びプログラム

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008143026A1 (ja) 2007-05-24 2008-11-27 Nec Corporation スループット推定方法及びシステム
US8098590B2 (en) 2008-06-13 2012-01-17 Qualcomm Incorporated Apparatus and method for generating performance measurements in wireless networks
JP5938832B2 (ja) 2012-10-01 2016-06-22 株式会社日立産機システム 無線通信システム、電力監視システム、および、無線通信システムの親機選択方法
JP6955144B2 (ja) 2017-06-29 2021-10-27 富士通株式会社 無線機器の設置位置決定装置、無線機器の設置位置決定方法及び無線機器の設置位置決定プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001320761A (ja) * 2000-03-01 2001-11-16 Yrp Mobile Telecommunications Key Tech Res Lab Co Ltd チャネル割当方法および装置
JP2004215075A (ja) * 2003-01-07 2004-07-29 Ntt Docomo Inc Ｃｄｍａ適応アンテナアレイ基地局設計方法および装置、ｃｄｍａ適応アンテナアレイ基地局設計プログラムおよびそのプログラムを記録した記録媒体
JP2012244563A (ja) * 2011-05-24 2012-12-10 Nippon Telegr & Teleph Corp <Ntt> 無線通信方法
JP2015050502A (ja) * 2013-08-30 2015-03-16 株式会社国際電気通信基礎技術研究所 無線通信装置、収容予測装置、境界情報取得装置、収容予測方法、境界情報取得方法、及びプログラム

Also Published As

Publication number Publication date
US20220132357A1 (en) 2022-04-28
US11785475B2 (en) 2023-10-10
JP7096499B2 (ja) 2022-07-06
JP2020145638A (ja) 2020-09-10

Similar Documents

Publication Publication Date Title
Wang et al. MmWave vehicle-to-infrastructure communication: Analysis of urban microcellular networks
Bai et al. Coverage and rate analysis for millimeter-wave cellular networks
Weiler et al. Quasi-deterministic millimeter-wave channel models in MiWEBA
Pätzold et al. A space‐time channel simulator for MIMO channels based on the geometrical one‐ring scattering model
JP4819303B2 (ja) 移動通信システムにおける基地局設置設計方法及び基地局設置設計装置並びにプログラム
Weng et al. Indoor massive MIMO channel modelling using ray-launching simulation
Ford et al. Markov channel-based performance analysis for millimeter wave mobile networks
Oehmann et al. SINR model with best server association for high availability studies of wireless networks
CN112584420B (zh) 一种信号测量方法、终端及网络侧设备
Zeman et al. Accuracy comparison of propagation models for mmWave communication in NS-3
CN114710747A (zh) 通导融合系统的室分节点布局方法及装置
Andrade et al. Radio channel spatial propagation model for mobile 3G in smart antenna systems
WO2020184246A1 (ja) 無線端末収容判定装置、無線端末収容判定方法、及びプログラム
Burghal et al. Band assignment in dual band systems: A learning-based approach
JP4071660B2 (ja) 伝搬環境模擬装置およびプログラム
WO2020184253A1 (ja) 到達率推定装置、到達率推定方法、及びプログラム
Xue et al. A probabilistic propagation graph modeling method for channel characterization in dense urban environments
Lecci et al. Accuracy vs. complexity for mmWave ray-tracing: A full stack perspective
Risi et al. Development and comparative analysis of path loss models using hybrid wavelet-genetic algorithm approach
JP7161136B2 (ja) 置局設計装置、置局設計方法及びプログラム
Gupta Physical Channel Models for Emerging Wireless Communication Systems
Comisso et al. Shannon capacity evaluation for 5G communications using the 3D random waypoint mobility model
Yamada et al. Plane-wave and vector-rotation approximation technique for reducing computational complexity to simulate MIMO propagation channel using ray-tracing
CN113346970B (zh) 面向无线三维信道的用户级信道空域特征建模方法
Joseph et al. Path Loss and Models: A Survey and Future Perspective for Wireless Communication Networks

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20770795

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20770795

Country of ref document: EP

Kind code of ref document: A1