JP7174249B2 - 到達率推定装置、到達率推定方法、及びプログラム - Google Patents

到達率推定装置、到達率推定方法、及びプログラム Download PDF

Info

Publication number
JP7174249B2
JP7174249B2 JP2019044959A JP2019044959A JP7174249B2 JP 7174249 B2 JP7174249 B2 JP 7174249B2 JP 2019044959 A JP2019044959 A JP 2019044959A JP 2019044959 A JP2019044959 A JP 2019044959A JP 7174249 B2 JP7174249 B2 JP 7174249B2
Authority
JP
Japan
Prior art keywords
arrival rate
received power
average received
calculated
base stations
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2019044959A
Other languages
English (en)
Other versions
JP2020150366A (ja
Inventor
一光 坂元
賢司 鈴木
陽平 片山
洋輔 藤野
浩之 福本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2019044959A priority Critical patent/JP7174249B2/ja
Priority to PCT/JP2020/008678 priority patent/WO2020184253A1/ja
Priority to US17/437,196 priority patent/US11991545B2/en
Publication of JP2020150366A publication Critical patent/JP2020150366A/ja
Application granted granted Critical
Publication of JP7174249B2 publication Critical patent/JP7174249B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/022—Site diversity; Macro-diversity
    • 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
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0802—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using antenna selection
    • H04B7/0805—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using antenna selection with single receiver and antenna switching
    • H04B7/0814—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using antenna selection with single receiver and antenna switching based on current reception conditions, e.g. switching to different antenna when signal level is below threshold
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L43/00—Arrangements for monitoring or testing data switching networks
    • H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
    • H04L43/0876—Network utilisation, e.g. volume of load or congestion level
    • H04L43/0894—Packet rate
    • 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
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/04—Arrangements for maintaining operational condition
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Mobile Radio Communication Systems (AREA)

Description

本発明は、到達率推定装置、到達率推定方法、及びプログラムに関する。
ＩｏＴ（Internet of Things；モノのインターネット）技術の発展により、各種センサ及び通信機能を備えたモノ（以下、「ＩｏＴ機器」又は「端末」という。）がネットワークに接続され、遠隔データ収集及びＩｏＴ機器の遠隔制御等を行うことが年々盛んになっている。また、近年、ＩｏＴ機器向けの無線通信方式として、ＬＰＷＡ（Low Power Wide Area network）が注目されている。ＬＰＷＡには、例えば、アンライセンスバンド（無線局免許を必要としない周波数帯）を用いるＬｏＲａＷＡＮ及びＳｉｇｆｏｘ、ライセンスバンド（無線局免許を必要とする周波数帯）を用いるＬＴＥ－Ｍ（Long Term Evolution for Machines）及びＮＢ（Narrow Band）－ＩｏＴ等の無線通信方式がある。
ＩｏＴ機器は、主に特定の場所に固定して設置され、利用されることが多い。一方、電波伝搬は、場所依存性が高い。これにより、基地局からの電波の受信電力が低い不感地帯にＩｏＴ機器が設置された場合、長期的に通信不可の状態となりうる。そのため、ＩｏＴ機器の設置予定場所において十分な受信電力を確保できるように置局設計を行うことが重要となる。
一般に、置局設計において基地局を配置したときの通信可能エリアを推定するために、電波伝搬シミュレータが用いられる。屋外評価用の電波伝搬シミュレータでは、３Ｄ（three-dimensional；３次元の）地図データ（地形高さ、建物高さ、及び土地利用分類等のデータ）が用いられ、例えば奥村－秦モデル等の電波伝搬モデルに基づいて平均受信電力が算出される（例えば、非特許文献１参照）。一般に、平均受信電力は、基地局からの距離に対する減衰（長区間変動）と、建物等の遮蔽物による場所的な受信電力変動（短区間変動）と、を考慮して算出される。また、平均受信電力は、基地局周辺エリアの微小なメッシュ（例えば、５ｍ間隔のメッシュ）ごとに算出される。そして、置局設計者は、電波伝搬シミュレータによって算出された平均受信電力の空間的分布を用いて、置局場所を選定する。
ところで、実際の受信点においては、建物等による反射、回折、及び透過を経て、複数の経路を介して到来する電波が合成されて受信される。そのため、周辺の車両や人間等の移動、及び樹木の揺れ等によって、到来する各電波の位相関係が変化する。これにより、実際の受信点においては、電波の受信電力は瞬時変動する。この瞬時変動は、レイリーフェージングと呼ばれる。
長区間変動、短区間変動、及び瞬時変動を考慮して平均受信電力を算出する方法として、レイトレース法による電波伝搬解析法がある。レイトレース法による電波伝搬解析法では、送信点から複数の経路（パス）を介して受信点へそれぞれ到来した複数の信号が算出され、当該複数の信号の位相関係に基づいて複素数として合成されることで、受信電力瞬時値が算出される。しかしながら、レイトレース法によって算出される、上記パスが異なる複数の信号の位相について信頼性が担保される場合とは、建物情報及び送受信点の位置の精度が、およそ波長の１／１０以下である場合に限られる。そのため、レイトレース法を用いる電波伝搬シミュレータは、一般に、複素数での合成を行わず、電力次元で各パスの信号を合成することで、瞬時変動を平均した平均受信電力を算出する。なお、一般に、無線通信では移動端末が前提であり、特定の位置における厳密な受信電力瞬時値の推定が不要であるため、置局設計では平均受信電力が用いられる。
"Atoll 3.3.0 Technical Reference Guide for Radio Networks", Forsk, AT330_TRR_E1, p.57-p.121, March 2015 岩井誠人著、「移動通信における電波伝搬 －無線通信シミュレーションのための基礎知識－」、株式会社コロナ社、ｐ．１６－ｐ．１９、２０１２年１１月２２日
上述したように、ＩｏＴ機器は、主に特定の場所に固定して設置される。そのため、設置予定場所への電波の到達性を確保しつつ、効率的な（例えば、コストが最小になる）置局設計を行うためには、平均受信電力の推定だけでは不十分であり、瞬時変動の影響を考慮した場合の電波の到達率を推定する必要がある。ここでいう電波の到達率とは、希望波の受信電力が受信局の受信感度以上となる確率である。
また、ＬＰＷＡでは、一般に、端末から基地局への上り通信において、端末から送信された信号は複数の基地局において受信され、復調される。受信に成功した信号は、上位のネットワークサーバへ送信される。そして、受信に成功した信号のうち、最も受信電力が高い信号がネットワークサーバによって選択される（サイトダイバーシチ）。また、基地局が複数のアンテナを備える場合、端末から送信された信号は複数のアンテナによって受信され、復調される。受信に成功した信号は、上位のネットワークサーバへ送られる。そして、受信に成功した信号のうち、最も受信電力が高い信号がネットワークサーバによって選択される（アンテナダイバーシチ）。基地局から端末への下り通信においては、直前の上り通信において選択された基地局及びアンテナから下り通信を行うことによって、高い通信品質が確保される。したがって、これらのダイバーシチ効果を考慮して到達率を推定する必要がある。
本発明は、このような事情を考慮してなされたものであり、電波の到達率をより高い精度で推定する技術を提供することを目的とする。
本発明の一態様は、複数の基地局及び前記基地局が備える複数のアンテナの各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部と、前記平均受信電力を、レイリーフェージングにより瞬時変動する振幅の確率密度関数から得られる受信電力瞬時値の累積分布関数に代入することで通信到達率を算出し、算出された通信到達率に基づいてサイトダイバーシチ及びアンテナダイバーシチ効果が考慮された通信到達率を算出する到達率推定部と、を備える到達率推定装置である。
また、本発明の一態様は、上記の到達率推定装置であって、前記到達率推定部は、受信電力の瞬時変動の影響を解析式によって算出する。
また、本発明の一態様は、上記の到達率推定装置であって、前記平均受信電力算出部は、上り通信の平均受信電力と下り通信の平均受信電力とを算出し、前記到達率推定部は、上り通信の到達率と下り通信の到達率の結合確率を算出する。
また、本発明の一態様は、上記の到達率推定装置であって、繰り返し判定部をさらに備え、前記平均受信電力算出部は、算出された平均受信電力の値が高い順に前記基地局をソートし、前記到達率推定部は、考慮する基地局数を電力の高い順に徐々に増やしながら到達率を繰り返し算出し、前記繰り返し判定部は、算出された到達率が目標値以上になった場合、又は１つ前の繰り返し処理で算出された到達率からの変化が収束目標値以内になった場合に前記繰り返し処理を終了させる。
また、本発明の一態様は、上記の到達率推定装置であって、厳密計算実施判定テーブル作成部と、厳密計算実施判定部と、をさらに備え、前記平均受信電力算出部は、算出された平均受信電力の値が高い順に前記基地局をソートし、前記厳密計算実施判定テーブル作成部は、平均受信電力と基地局数とから厳密な到達率計算の実施要否を判定するための厳密計算実施判定テーブルを作成し、前記厳密計算実施判定部は、前記厳密計算実施判定テーブルを参照して厳密計算の実施要否を判定する。
また、本発明の一態様は、上記の到達率推定装置であって、グルーピング部をさらに備え、前記平均受信電力算出部は、算出された平均受信電力の値が高い順に前記基地局をソートし、前記グルーピング部は、隣接メッシュと上位数基地局の平均受信電力の差分が小さい場合に前記メッシュのグルーピングを行い、前記到達率推定部は、グループ内の１つのメッシュについて瞬時変動の影響とダイバーシチ効果を考慮して到達率を算出し、算出された前記到達率を前記グループ内の他のメッシュについての到達率と見なす。
また、本発明の一態様は、複数の基地局及び前記基地局が備える複数のアンテナの各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出ステップと、前記平均受信電力を、レイリーフェージングにより瞬時変動する振幅の確率密度関数から得られる受信電力瞬時値の累積分布関数に代入することで通信到達率を算出し、算出された通信到達率に基づいてサイトダイバーシチ及びアンテナダイバーシチ効果が考慮された通信到達率を算出する到達率推定ステップと、を有する到達率推定方法である。
また、本発明の一態様は、上記の到達率推定装置としてコンピュータを機能させるためのプログラムである。
本発明によれば、電波の到達率をより高い精度で推定するができる。
本発明の第１の実施形態に係る到達率推定装置１ａの機能構成を示すブロック図である。 本発明の第１の実施形態に係る到達率推定装置１ａによる到達率推定処理を説明するための図である。 本発明の第１の実施形態に係る到達率推定装置１ａの平均受信電力算出部１１によって算出される平均受信電力を示す図である。 本発明の第１の実施形態に係る到達率推定装置１ａの動作を示すフローチャートである。 本発明の第２の実施形態に係る到達率推定装置１ｂの機能構成を示すブロック図である。 本発明の第２の実施形態に係る到達率推定装置１ｂによる到達率推定処理を説明するための図である。 本発明の第２の実施形態に係る到達率推定装置１ｂの平均受信電力算出部１１によって算出される平均受信電力を示す図である。 本発明の第２の実施形態に係る到達率推定装置１ｂの動作を示すフローチャートである。 本発明の第３の実施形態に係る到達率推定装置１ｃの機能構成を示すブロック図である。 本発明の第３の実施形態に係る到達率推定装置１ｃが用いる厳密計算実施判定テーブルの一例を示す図である。 本発明の第３の実施形態に係る到達率推定装置１ｃの動作を示すフローチャートである。 本発明の第４の実施形態に係る到達率推定装置１ｄの機能構成を示すブロック図である。 本発明の第４の実施形態に係る到達率推定装置１ｄの動作を示すフローチャートである。
以下、図面を参照して、本発明の実施形態について詳しく説明する。
＜第１の実施形態＞
以下、本発明の第１の実施形態について説明する。
［到達率推定装置の機能構成］
以下、本実施形態に係る到達率推定装置１ａの機能構成について説明する。
図１は、本発明の第１の実施形態に係る到達率推定装置１ａの機能構成を示すブロック図である。図１に示すように、到達率推定装置１ａは、平均受信電力算出部１１と、無線設備データベース１２と、到達率推定部１３ａと、を具備する。
平均受信電力算出部１１は、一般的な電波伝搬シミュレータに相当するものである。平均受信電力算出部１１は、地図データ（例えば、地形の高さ、建物の高さ、土地の利用分類等のデータ）を外部の機器等からインポートする。平均受信電力算出部１１は、設計者によって設定されたエリアの地図データを、微小なメッシュ（例えば、５ｍ間隔のメッシュ）に分割する。そして、平均受信電力算出部１１は、基地局から送出される下り通信の電波の、各メッシュ（すなわち、各メッシュに相当する位置に設置されているものと仮定した仮想端末）における平均受信電力を算出する。同様に、平均受信電力算出部１１は、各メッシュから基地局への上り通信の電波の平均受信電力を算出する。
無線設備データベース１２は、各基地局の、アンテナ本数、及び受信感度等の無線設備に関わるパラメータを保持する。
到達率推定部１３ａは、平均受信電力算出部１１によって算出された、各基地局と各メッシュとの間の上り通信及び下り通信の平均受信電力に対して、瞬時変動及びダイバーシ
チ効果の影響を考慮して、電波の到達率を算出する。
［到達率推定］
以下、電波の到達率の推定処理について説明する。
図２は、本発明の第１の実施形態に係る到達率推定装置１ａによる到達率推定処理を説明するための図である。
ここでは、図２に示すメッシュ＃１の位置における、電波の到達率の推定処理を一例として説明する。また、図２に示すエリアに存在する基地局の数は、Ｉ局であるものとする。なお、図２においては、基地局を４局のみ図示しており、その他の基地局の記載は省略している。また、図２に示す各基地局が備えるアンテナの数は、いずれも２基であるものとする。
ここで、基地局＃ｉ（ｉ＝１～Ｉ）のアンテナ＃ｊ（ｊ＝１～２）における上り通信の平均受信電力をＲｕｐ（ｉ，ｊ）と表す。また、基地局＃ｉのアンテナ＃ｊからの下り通信の平均受信電力をＲｄｏｗｎ（ｉ，ｊ）と表す。
図３は、本発明の第１の実施形態に係る到達率推定装置１ａの平均受信電力算出部１１によって算出される平均受信電力を示す図である。平均受信電力算出部１１は、図３（Ａ）に示す、基地局＃ｉとメッシュ＃１との間の上り通信の平均受信電力であるＲｕｐ（ｉ，ｊ）を、基地局ごとに算出する。また、平均受信電力算出部１１は、図３（Ｂ）に示す、基地局＃ｉとメッシュ＃１との間の下り通信の平均受信電力であるＲｄｏｗｎ（ｉ，ｊ）を、基地局ごとに算出する。
以下、到達率推定部１３ａによる到達率の推定処理について説明する。まず、上り通信について説明する。到達率推定部１３ａは、基地局＃ｉのアンテナ数、及び受信感度を示す情報を、無線設備データベース１２から取得する。そして、到達率推定部１３ａは、上り通信の平均受信電力であるＲｕｐ（ｉ，ｊ）に対してレイリーフェージングによる瞬時変動の影響を考慮した場合における、受信電力瞬時値が基地局＃ｉの受信感度Ｔｕｐ（ｉ）以上となる確率を、基地局＃ｉのアンテナ＃ｊの上り通信到達率Ｐｕｐ（ｉ，ｊ）として算出する。
ここで、受信電力瞬時値を算出する単純な方法として、モンテカルロシミュレーション法がある。モンテカルロシミュレーション法では、基地局＃ｉのアンテナ＃ｊに複数の経路で到来する電波の位相をランダムに変化させて合成することで、受信電力瞬時値を算出する。そして、これを繰り返し行うことによって、受信電力瞬時値の累積分布関数（ＣＤＦ；Cumulative Distribution Function）を作成する。しかしながら、この繰り返し計算は、一般的に膨大な回数（例えば数万回以上）行われる。そのため、膨大な計算時間が必要になるという課題がある。そこで、本実施形態では、レイリーフェージングの瞬時変動を下記の解析式で表すことによって、モンテカルロシミュレーションを行うことなく、受信電力瞬時値のＣＤＦを得る。
受信信号の振幅をａとした場合、レイリーフェージングにより瞬時変動する振幅ａの確率密度関数（ＰＤＦ；Probability Density Function）は、平均受信電力Ｒｕｐ（ｉ，ｊ）を用いて、以下の式（１）で表すことができる（非特許文献２参照）。
Figure 0007174249000001
これを電力ｒ＝ａ２のＰＤＦに変換することで、受信電力瞬時値のＰＤＦは、以下の式（２）で表すことができる。
Figure 0007174249000002
上記のＰＤＦからＣＤＦを算出すると、ＣＤＦは、以下の式（３）で表すことができる。
Figure 0007174249000003
基地局＃ｉのアンテナ＃ｊの上り通信到達率Ｐｕｐ（ｉ，ｊ）は、式（３）のｙに、基地局＃ｉの受信感度Ｔｕｐ（ｉ）を入力することによって、算出することができる。基地局＃ｉのアンテナ＃ｊの上り通信到達率Ｐｕｐ（ｉ，ｊ）は、以下の式（４）で表すことができる。
Figure 0007174249000004
以上により、レイリーフェージングの瞬時変動を解析式で表すことによって、モンテカルロシミュレーションを行うことなく、基地局＃ｉのアンテナ＃ｊの上り通信到達率Ｐｕｐ（ｉ，ｊ）を短時間で算出することができる。
次に、各アンテナで受信成功した信号のうち、最も受信電力が高い信号を選択する際のアンテナダイバーシチ効果を考慮すると、基地局＃ｉの上り通信到達率Ｐ'ｕｐ（ｉ）は以下の式（５）式で表すことができる。
Figure 0007174249000005
さらに、複数の基地局で受信成功した信号のうち、最も受信電力が高い信号を選択する際のサイトダイバーシチ効果を考慮すると、上り通信到達率Ｐ''ｕｐは、以下の式（６）で表すことができる。
Figure 0007174249000006
以下、上り通信と下り通信の一連の通信シーケンスを想定する場合の到達率の算出について説明する。ここでは、上り通信と下り通信とがそれぞれ１回ずつ行われる場合を想定する。すなわち、上り通信及び下り通信の一連の通信において高い通信品質が確保されるように基地局及びアンテナが選択される仕組みを有する場合において、到達率を推定する場合について説明する。上り通信と下り通信とがそれぞれ１回ずつ行われる場合とは、例えば、ＩｏＴ機器が取得したセンサデータをアップロードして、基地局がＩｏＴ機器へＡｃｋ（Acknowledge;肯定応答）信号を返送する場合等である。
到達率推定部１３ａは、基地局＃ｉのアンテナ数と端末の受信感度とを、無線設備データベース１２から取得する。到達率推定部１３ａは、下り通信の平均受信電力であるＲｄｏｗｎ（ｉ，ｊ）に対してレイリーフェージングによる瞬時変動の影響を考慮した場合における、受信電力瞬時値が端末の受信感度Ｔｄｏｗｎ以上となる確率を、基地局＃ｉのアンテナ＃ｊの下り通信到達率Ｐｄｏｗｎ（ｉ，ｊ）として算出する。上り通信の到達率推定と同様に、レイリーフェージングの瞬時変動を解析式で表すことで、モンテカルロシミュレーションを行うことなく、基地局＃ｉのアンテナ＃ｊの下り通信到達率Ｐｄｏｗｎ（ｉ，ｊ）を短時間で算出することができる。
レイリーフェージングにより瞬時変動する下り通信の電力のＣＤＦは、平均受信電力Ｒｄｏｗｎ（ｉ，ｊ）を用いて、以下の式（７）で表すことができる。
Figure 0007174249000007
基地局＃ｉのアンテナ＃ｊの下り通信到達率Ｐｄｏｗｎ（ｉ，ｊ）は、式（７）のｙに、端末の受信感度Ｔｄｏｗｎを入力することによって算出することができる。基地局＃ｉのアンテナ＃ｊの下り通信到達率Ｐｄｏｗｎ（ｉ，ｊ）は、以下の式（８）で表すことができる。
Figure 0007174249000008
アンテナダイバーシチ効果を考慮すると、基地局＃ｉの下り通信到達率Ｐ'ｄｏｗｎ（ｉ）は、以下の式（９）で表すことができる。
Figure 0007174249000009
基地局＃ｉの、上り通信の到達率と下り通信の到達率との結合確率Ｐｊｏｉｎｔ（ｉ）は、以下の式（１０）で表すことができる。
Figure 0007174249000010
さらに、サイトダイバーシチ効果を考慮すると、上り通信の到達率と下り通信の到達率との結合確率Ｐ'ｊｏｉｎｔは、以下の式（１１）で表すことができる。
Figure 0007174249000011
［到達率推定装置の動作］
以下、到達率推定処理における到達率推定装置１ａの動作について説明する。
図４は、本発明の第１の実施形態に係る到達率推定装置１ａの動作を示すフローチャートである。なお、ここでは上り通信の到達率推定について説明する。
平均受信電力算出部１１は、設計者によって指定されたエリアの地図データを微小なメッシュに分割する。ここで、メッシュ数をＭとする。平均受信電力算出部１１は、分割された各メッシュから各基地局への上り通信の平均受信電力を算出する（ステップＳ００１）。
到達率推定部１３ａは、無線設備データベース１２から、各基地局のアンテナ数及び受信感度を示す情報を取得する（ステップＳ００２）。
メッシュを識別する変数ｍに１が代入される（ステップＳ００３）。
到達率推定部１３ａは、メッシュ＃ｍの平均受信電力を用いて、レイリーフェージングによる瞬時変動の影響を考慮した際の、当該メッシュ＃ｍについての受信電力瞬時値のＣＤＦを作成する（ステップＳ００４）。
到達率推定部１３ａは、作成されたＣＤＦに基づいて、基地局の受信感度以上となる確率を算出する。到達率推定部１３ａは、アンテナダイバーシチ及びサイトダイバーシチの効果を考慮して、メッシュ＃ｍの上り通信の到達率を算出する（ステップＳ００５）。
変数ｍがメッシュ数Ｍと等しいか否かについての判定がなされる（ステップＳ００６）。変数ｍがメッシュ数Ｍと等しくない場合（ステップＳ００６－Ｎｏ）、変数ｍに１が加算され（ステップＳ００７）、ステップＳ００４に移行する。一方、変数ｍがメッシュ数Ｍと等しい場合（ステップＳ００６－Ｙｅｓ）、図４のフローチャートが示す到達率推定装置１ａの動作が終了する。
以上説明したように、本発明の第１の実施形態に係る到達率推定装置１ａは、複数の基地局及び前記基地局が備える複数のアンテナの各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部１１と、前記平均受信電力を、レイリーフェージングにより瞬時変動する振幅の確率密度関数から得られる受信電力瞬時値の累積分布関数に代入することで通信到達率を算出し、算出された通信到達率に基づいてサイトダイバーシチ及びアンテナダイバーシチ効果が考慮された通信到達率を算出する到達率推定部１３ａとを備える。
このような構成を備えることによって、本発明の第１の実施形態に係る到達率推定装置１ａは、電波の到達率をより高い精度で推定することができる。
＜第２の実施形態＞
以下、本発明の第２の実施形態について説明する。
上述した第１の実施形態では、全ての基地局によるサイトダイバーシチ効果を考慮して到達率を算出する場合について説明した。これに対し、以下、第２の実施形態では、到達率への寄与度が高い基地局のみを計算対象として到達率を算出する場合について説明する。また、以下、第１の実施形態との相違点についてのみ説明する。なお、ここでは上り通信の到達率推定について説明する。
［到達率推定装置の構成］
以下、本実施形態に係る到達率推定装置１ｂの機能構成について説明する。
図５は、本発明の第２の実施形態に係る到達率推定装置１ｂの機能構成を示すブロック図である。図５に示すように、到達率推定装置１ｂは、平均受信電力算出部１１と、無線設備データベース１２と、到達率推定部１３ｂと、繰り返し判定部１４と、を具備する。このように、繰り返し判定部１４を具備する点が、第１の実施形態に係る到達率推定装置１ａの機能構成との相違点である。
到達率推定部１３ｂは、平均受信電力算出部１１によって算出された、平均受信電力Ｒｕｐ（ｉ，ｊ）に基づいて、平均受信電力の値が高い順に基地局をソートする。そして、到達率推定部１３ｂは、到達率推定において計算対象とする基地局を、平均受信電力が高い順に徐々に増やしながら、到達率を繰り返し算出する。そして、到達率推定部１３ｂは、繰り返し判定部１４によって繰り返し処理終了の判定がなされた場合、繰り返し処理を終了する。そして、到達率推定部１３ｂは、その時点において算出された到達率を、計算対象のメッシュの位置における到達率として出力する。
繰り返し判定部１４は、算出された到達率が置局設計における目標値以上となった場合、又は、１つ前の繰り返し処理で算出された到達率からの変化量が所定値以内になった場合に、繰り返し処理終了の判定を行う。
［到達率推定］
以下、電波の到達率の推定処理について説明する。
図６は、本発明の第２の実施形態に係る到達率推定装置１ｂによる到達率推定処理を説明するための図である。
ここでは、図６に示すメッシュ＃１の位置における電波の到達率を推定を推定する場合を一例として説明する。また、図６に示すエリアに存在する基地局の数は、Ｉ局であるものとする。なお、図６においては、基地局を５局のみ図示しており、残りの基地局の記載は省略している。また、図６に示す各基地局が備えるアンテナの数は、いずれも２基であるものとする。
まず、設計者によって、置局設計における到達率の目標値Ｐｔａｒｇｅｔ（例えば、０．９等の値）と繰り返し処理における到達率の変化の収束目標値ΔPｔａｒｇｅｔ（例えば、０．０２等の値）とが決定される。そして、これら決定されたＰｔａｒｇｅｔの値及びΔPｔａｒｇｅｔの値が、繰り返し判定部１４に設定される。
平均受信電力算出部１１は、基地局＃ｉ（ｉ＝１～Ｉ）のアンテナ＃ｊ（ｊ＝１～２）における平均受信電力Ｒｕｐ（ｉ，ｊ）を算出する。
到達率推定部１３ｂは、基地局＃ｉのアンテナ数及び受信感度を示す情報を、無線設備データベース１２から取得する。次に、到達率推定部１３ｂは、Ｒｕｐ（ｉ，ｊ）の値の高い順に基地局をソートする。ここでは、一例として、図７に示すリストの上から順に、平均受信電力が高い基地局であるものとする。なお、ここでは、基地局が備える２本のアンテナについてそれぞれ平均受信電力が算出されるが、平均受信電力の高い方を当該基地局の代表値として基地局をソートしてもよいし、２本のアンテナの平均受信電力を平均した値によって基地局をソートしてもよい。
以下、到達率の計算の対象とする基地局を、繰り返し処理ごとに１局ずつ増やしていく場合について説明する。すなわち、ｋ回目の繰り返し処理において到達率の計算の対象とする基地局は、平均受信電力の値が高いほうから上位ｋ局の基地局（上位数基地局）である。なお、繰り返し処理ごとに増やす局数は１局に限られるものではなく、２局ずつ増やしたり、例えば、１局、２局、４局、・・・、のように、指数関数的に増やしたりしてもよい。
上述した第１の実施形態と同様の方法により、到達率推定部１３ｂは、ｋ回目の繰り返し処理における到達率Ｐ''ｕｐ（ｋ）を算出する。到達率推定部１３ｂは、繰り返し判定部１４によって繰り返し処理終了の判定がなされるまで、計算対象の基地局を１局ずつ増やしながら到達率の算出を繰り返す。繰り返し判定部１４によって繰り返し処理終了の判定がなされた場合、到達率推定部１３ｂは、繰り返し処理を終了し、その時点において算出された到達率Ｐ''ｕｐ（ｋ）を、計算対象のメッシュの位置における到達率として出力する。
繰り返し判定部１４は、ｋ回目の繰り返し処理で算出された到達率Ｐ''ｕｐ（ｋ）が、以下の式（１２）を満たすか否かを判定する。
Figure 0007174249000012
式（１２）を満たす場合、もしくはｋ－１回目の（１つ前の）繰り返し処理において算出された到達率Ｐｕｐ''（ｋ－１）からの変化量が以下の式（１３）を満たす場合に、繰り返し判定部１４は、繰り返し処理終了の判定を行う。
Figure 0007174249000013
繰り返し判定部１４によって繰り返し処理終了の判定がなされると、到達率推定部１３ｂは、その回において算出された到達率Ｐ''ｕｐ（ｋ）を、計算対象のメッシュの位置における到達率として出力する。
［到達率推定装置の動作］
以下、到達率推定処理における到達率推定装置１ｂの動作について説明する。
図８は、本発明の第２の実施形態に係る到達率推定装置１ｂの動作を示すフローチャートである。
まず、設計者によって、置局設計における到達率の目標値Ｐｔａｒｇｅｔと、繰り返し処理における到達率の変化の収束目標値ΔＰｔａｒｇｅｔと、が決定される。そして、これら決定されたＰｔａｒｇｅｔの値及びΔＰｔａｒｇｅｔの値が、繰り返し判定部１４に設定される（ステップＳ１０１）。
平均受信電力算出部１１は、設計者によって指定されたエリアの地図データを微小なメッシュに分割する。ここで、メッシュ数をＭとする。平均受信電力算出部１１は、分割された各メッシュについて、各基地局における平均受信電力を算出する（ステップＳ１０２）。
到達率推定部１３ｂは、無線設備データベース１２から、各基地局のアンテナ数及び受信感度を示す情報を取得する（ステップＳ１０３）。
メッシュを識別する変数ｍに１が代入される（ステップＳ１０４）。
到達率推定部１３ｂは、メッシュ＃ｍについて、平均受信電力の値が高い順に基地局をソートする（ステップＳ１０５）。
繰り返し処理の実行回数を識別する変数ｋに、１が代入される（ステップＳ１０６）。
到達率推定部１３ｂは、メッシュ＃ｍの平均受信電力を用いて、レイリーフェージングによる瞬時変動の影響を考慮した際の、当該メッシュ＃ｍについての受信電力瞬時値のＣＤＦを作成する。
到達率推定部１３ｂは、ＣＤＦから、基地局の受信感度以上となる確率を算出する。到達率推定部１３ｂは、アンテナダイバーシチ、及び、平均受信電力の高いほうから上位ｋ局の基地局によるサイトダイバーシチの効果を考慮して、到達率を算出する（ステップＳ１０７）。
到達率推定部１３ｂは、上記の式（１１）及び式（１２）のうち少なくとも一つを満たすか否かを判定する（ステップＳ１０８）。式（１１）及び式（１２）のいずれも満たさない場合（ステップＳ１０８－Ｎｏ）、変数ｋに１が加算され（ステップＳ１０９）、ステップＳ１０７に移行する。一方、式（１１）及び式（１２）のうち少なくとも一方を満たす場合（ステップＳ１０８－Ｙｅｓ）、当該メッシュにおける到達率の算出を終了し、ステップＳ１１０に移行する。
変数ｍがメッシュ数Ｍと等しいか否かについての判定がなされる（ステップＳ１１０）。変数ｍがメッシュ数Ｍと等しくない場合（ステップＳ１１０－Ｎｏ）、変数ｍに１が加算され（ステップＳ１１１）、ステップＳ１０５に移行する。一方、変数ｍがメッシュ数Ｍと等しい場合（ステップＳ１１０－Ｙｅｓ）、図８のフローチャートが示す到達率推定装置１ｂの動作が終了する。
以上説明したように、第２の実施形態に係る到達率推定装置１ｂは、到達率への寄与度の高い（すなわち、平均受信電力の高い）基地局のみを計算対象として到達率を算出するため、到達率推定に要する計算時間をより短縮させることができる。
＜第３の実施形態＞
以下、本発明の第３の実施形態について説明する。
上述した第１の実施形態及び第２の実施形態では、全てのメッシュについて、レイリーフェージングによる瞬時変動及びダイバーシチ効果を考慮した厳密な到達率計算を行う場合について説明した。これに対し、以下、第３の実施形態では、事前処理として、各メッシュについて到達率を厳密に算出するか否かを判定し、厳密計算が必要なメッシュについてのみ、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率計算を行う場合について説明する。また、以下、第１の実施形態との相違点についてのみ説明する。なお、ここでは上り通信の到達率推定について説明する。
［到達率推定装置の構成］
以下、本実施形態に係る到達率推定装置１ｃの機能構成について説明する。
図９は、本発明の第３の実施形態に係る到達率推定装置１ｃの機能構成を示すブロック図である。図９に示すように、到達率推定装置１ｃは、平均受信電力算出部１１と、無線設備データベース１２と、到達率推定部１３ｃと、厳密計算実施判定テーブル作成部１５と、厳密計算実施判定テーブル記憶部１６と、厳密計算実施判定部１７と、を具備する。このように、厳密計算実施判定テーブル作成部１５と、厳密計算実施判定テーブル記憶部１６と、厳密計算実施判定部１７と、を具備する点が、第１の実施形態における到達率推定装置１ａの機能構成との相違点である。
厳密計算実施判定テーブル作成部１５は、事前処理として、平均受信電力が受信感度からｘ［ｄＢ］のマージンがある基地局がｙ局ある場合の到達率を算出する。そして、厳密計算実施判定テーブル作成部１５は、上記の厳密計算が必要であるか否か（厳密計算の実施要否）の判定に用いられるテーブルである、厳密計算実施判定テーブルを作成する。厳密計算実施判定テーブル作成部１５は、算出された到達率が置局設計における到達率の目標値Ｐｔａｒｇｅｔ以上である場合には（上記の厳密計算が不要であることを表す）「不要」を示す情報を、及び、算出された到達率が置局設計における到達率の目標値Ｐｔａｒｇｅｔ未満である場合には（上記の厳密計算が必要であることを表す）「必要」を示す情報を、厳密計算実施判定テーブルに設定する。厳密計算実施判定テーブル作成部１５は、作成した厳密計算実施判定テーブルを、厳密計算実施判定テーブル記憶部１６に記憶させる。
厳密計算実施判定部１７は、平均受信電力算出部１１によって算出された平均受信電力Ｒｕｐ（ｉ，ｊ）を、値が高い順にソートする。そして、厳密計算実施判定部１７は、厳密計算実施判定テーブルに設定された値と、平均受信電力Ｒｕｐ（ｉ，ｊ）の値が上位の基地局の平均受信電力と、に基づいて、上記の厳密計算が必要であるか否かの判定結果（「不要」又は「必要」）を示す情報を到達率推定部１３ｃに出力する。
到達率推定部１３ｃは、厳密計算実施判定部１７からの出力が「不要」を示す情報である場合には、目標値Ｐｔａｒｇｅｔを到達率の推定結果として出力する。また、到達率推定部１３ｃは、厳密計算実施判定部１７からの出力が「必要」を示す情報である場合には、上述した第１の実施形態又は第２の実施形態と同様の方法を用いて、瞬時変動やダイバーシチ効果を考慮した厳密な到達率計算を行う。
［到達率推定］
以下、電波の到達率の推定処理について説明する。
図１０は、本発明の第３の実施形態に係る到達率推定装置１ｃが用いる厳密計算実施判定テーブルの一例を示す図である。
まず、設計者によって、置局設計における到達率の目標値Ｐｔａｒｇｅｔ（例えば、０．９等の値）が決定される。決定された到達率の目標値Ｐｔａｒｇｅｔは、厳密計算実施判定テーブル作成部１５に設定される。
厳密計算実施判定テーブル作成部１５は、事前処理として、平均受信電力が受信感度からｘ［ｄＢ］のマージンがある基地局がｙ局存在する場合の到達率を算出する。図１０は、ｘの値が５ｄＢ間隔で、基地局数ｙが１～４局である場合の例を示す。
例えば、ｘ＝１５［ｄＢ］、及びｙ＝３である場合には、基地局の受信感度＋１５［ｄＢ］の平均受信電力となる基地局が３局あるという想定のもとで、上述した第１の実施形態と同様の方法により、瞬時変動及びダイバーシチ効果が考慮された到達率が算出される。ここでは、算出された到達率は目標値Ｐｔａｒｇｅｔ以上であるものとする。この場合、厳密計算実施判定テーブル作成部１５は、厳密計算実施判定テーブルに対して「不要」を示す情報を設定する。
また、例えば、ｘ＝５［ｄＢ］、及びｙ＝２である場合には、基地局の受信感度＋５［ｄＢ］の平均受信電力となる基地局が２局存在するという想定のもとで、上述した第１の実施形態と同様の方法により、瞬時変動及びダイバーシチ効果が考慮された到達率が算出される。ここでは、算出された到達率は目標値Ｐｔａｒｇｅｔ未満であるものとする。この場合、厳密計算実施判定テーブル作成部１５は、厳密計算実施判定テーブルに対して「必要」を示す情報を設定する。以上の手順により、厳密計算実施判定テーブル作成部１５は、厳密計算実施判定テーブルを作成する。なお、図１０は、基地局が４局である場合において作成された厳密計算実施判定テーブルの一例である。
厳密計算実施判定部１７は、平均受信電力算出部１１によって算出された平均受信電力Ｒｕｐ（ｉ，ｊ）を、値の高い順にソートする。そして、厳密計算実施判定部１７は、それぞれの平均受信電力について、基地局の受信感度に対して何［ｄＢ］のマージンがあるかについての計算を行う。
例えば、平均受信電力Ｒｕｐ（ｉ，ｊ）の値が高い順に、マージンが１８［ｄＢ］、１６［ｄＢ］、４［ｄＢ］、及び３［ｄＢ］である場合、平均受信電力Ｒｕｐ（ｉ，ｊ）の値が１５［ｄＢ］以上２０［ｄＢ］未満の基地局は２局である。そのため、厳密計算実施判定部１７は、厳密計算実施判定テーブルにおけるｘ＝１５［ｄＢ］かつ基地局２局である場合に相当するセルに設定された情報（「不要」）を読み取って、当該情報を到達率推定部１３ｃへ出力する。到達率推定部１３ｃは、厳密計算実施判定部から「不要」を示す情報が出力されると、目標値Ｐｔａｒｇｅｔを到達率の推定結果として外部の機器等へ出力する。
また、例えば、平均受信電力Ｒｕｐ（ｉ，ｊ）の値が高い順にマージンが１８［ｄＢ］、６［ｄＢ］、４［ｄＢ］、３［ｄＢ］である場合、平均受信電力Ｒｕｐ（ｉ，ｊ）の値が１５［ｄＢ］以上２０［ｄＢ］未満の基地局は１局である。そのため、厳密計算実施判定部１７は、厳密計算実施判定テーブルにおけるｘ＝１５［ｄＢ］かつ基地局１局である場合に相当するセルに設定された情報（「必要」）を読み取って、当該情報を到達率推定部１３ｃへ出力する。到達率推定部１３ｃは、厳密計算実施判定部から「必要」を示す情報が出力されると、上述した第１の実施形態又は第２の実施形態と同様の方法により、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率の計算を行う。
なお、図１０に示した厳密計算実施判定テーブルにおいて、ｘ＝０ｄＢの欄に着目した場合、０［ｄＢ］以上の基地局は４局であるため「不要」と判定される。このように、厳密計算実施判定テーブルにおいて「必要」と「不要」とが混在する場合も発生し得る。しかしながら、この場合、「必要」として扱うのか「不要」として扱うのかについて、予め設計指針に応じて定義しておけばよい。
例えば、多少の計算時間を犠牲にして、曖昧なメッシュについては全て厳密に計算すべきであるという設計指針である場合には、厳密計算実施判定部１７は、例えば、当てはまる複数のパターンのうち１つでも「必要」を示す情報があれば、「必要」を示す情報を到達率推定部１３ｃに出力するようにすればよい。逆に、計算時間の短さを優先すべきであるという設計指針である場合には、厳密計算実施判定部１７は、例えば、当てはまる複数のパターンのうち１つでも「不要」を示す情報があれば、「不要」を示す情報を到達率推定部１３ｃ出力し、全てが「必要」を示す情報であるときにのみ、「必要」を示す情報を到達率推定部１３ｃに出力するようにすればよい。
［到達率推定装置の動作］
以下、到達率推定処理における到達率推定装置１ｃの動作について説明する。
図１１は、本発明の第３の実施形態に係る到達率推定装置１ｃの動作を示すフローチャートである。
まず、設計者によって、置局設計における到達率の目標値Ｐｔａｒｇｅｔが決定される。そして、決定された到達率の目標値Ｐｔａｒｇｅｔが、厳密計算実施判定テーブル作成部１５に設定される（ステップＳ２０１）。
厳密計算実施判定テーブル作成部１５は、上述した手順によって厳密計算実施判定テーブルを作成する（ステップＳ２０２）。
平均受信電力算出部１１は、設計者によって指定されたエリアの地図データを微小なメッシュに分割する。ここで、メッシュ数をＭとする。平均受信電力算出部１１は、分割された各メッシュについて、各基地局における平均受信電力を算出する（ステップＳ２０３）。
厳密計算実施判定部１７及び到達率推定部１３ｃは、無線設備データベース１２から、各基地局のアンテナ数及び受信感度を示す情報を取得する（ステップＳ２０４）。
メッシュを識別する変数ｍに１が代入される（ステップＳ２０５）。
厳密計算実施判定部１７は、メッシュ＃ｍについて、平均受信電力の値が高い順に基地局をソートする。そして、厳密計算実施判定部１７は、それぞれの平均受信電力について、基地局の受信感度に対して何［ｄＢ］のマージンがあるかについての計算を行う（ステップＳ２０６）。そして、厳密計算実施判定部１７は、厳密計算実施判定テーブルを参照し、算出されたマージンに基づいて、上述した厳密計算が必要であるか否かの判定結果を示す情報（すなわち、「必要」を示す情報又は「不要」を示す情報）を得る（ステップＳ２０７）。厳密計算実施判定部１７は、当該判定結果を示す情報を、到達率推定部１３ｃへ出力する。
厳密計算が必要であるか否かの判定結果が「不要」を示す情報である場合（ステップＳ２０８－Ｙｅｓ）、到達率推定部１３ｃは、目標値Ｐｔａｒｇｅｔを到達率の推定結果として出力する（ステップＳ２０９）。一方、厳密計算が必要であるか否かの判定結果が「必要」を示す情報である場合（ステップＳ２０８－Ｎｏ）、到達率推定部１３ｃは、上述した第１の実施形態又は第２の実施形態と同様の方法により、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率の計算を行う（ステップＳ２１０）。
変数ｍがメッシュ数Ｍと等しいか否かについての判定がなされる（ステップＳ２１１）。変数ｍがメッシュ数Ｍと等しくない場合（ステップＳ２１１－Ｎｏ）、変数ｍに１が加算され（ステップＳ２１２）、ステップＳ２０６に移行する。一方、変数ｍがメッシュ数Ｍと等しい場合（ステップＳ２１１－Ｙｅｓ）、図１１のフローチャートが示す到達率推定装置１ｃの動作が終了する。
以上説明したように、第３の実施形態に係る到達率推定装置１ｃは、事前処理として、各メッシュについて、到達率を厳密に算出する必要があるか否かについての判定をそれぞれ行う。そして、到達率推定装置１ｃは、到達率を厳密に算出する必要があると判定されたメッシュについてのみ、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率の計算を行う。そのため、第３の実施形態に係る到達率推定装置１ｃは、到達率推定に要する計算時間をより短縮させることができる。
＜第４の実施形態＞
以下、本発明の第４の実施形態について説明する。
上述した第１の実施形態及び第２の実施形態では、全メッシュについて、レイリーフェージングによる瞬時変動及びダイバーシチ効果を考慮した厳密な到達率計算を行う場合について説明した。これに対し、以下、第４の実施形態では、隣接メッシュ同士で上位数局の平均受信電力の値を比較し、差分が小さい場合、これら隣接メッシュをグルーピングし、１つの代表メッシュについてのみ、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率計算を行う場合について説明する。また、以下、第１の実施形態との相違点についてのみ説明する。なお、ここでは上り通信の到達率推定について説明する。
［到達率推定装置の構成］
以下、本実施形態に係る到達率推定装置１ｄの機能構成について説明する。
図１２は、本発明の第４の実施形態に係る到達率推定装置１ｄの機能構成を示すブロック図である。図１２に示すように、到達率推定装置１ｄは、平均受信電力算出部１１と、無線設備データベース１２と、到達率推定部１３ｄと、グルーピング部１８と、を具備する。このように、グルーピング部１８を具備する点が、第１の実施形態における到達率推定装置１ａの機能構成との相違点である。
グルーピング部１８は、平均受信電力算出部１１によって算出された平均受信電力の値が高い順に基地局をソートする。グルーピング部１８は、上位所定数局の平均受信電力と、隣接メッシュにおける同一の基地局の平均受信電力とを比較し、差分が小さい場合、これら隣接メッシュをグルーピングする。
到達率推定部１３ｄは、グルーピングされた複数のメッシュのうち、１つの代表メッシュについて、上述した第１の実施形態又は第２の実施形態と同様の方法を用いて、瞬時変動やダイバーシチ効果を考慮した厳密な到達率計算を行う。そして、到達率推定部１３ｄは、グルーピングされたその他のメッシュに対して、代表メッシュに対する計算結果を適用する（すなわち、グルーピングされたその他のメッシュについての到達率は、１つの代表メッシュについて算出された到達率と同一であるものと見なされる）。
なお、代表メッシュは、対象範囲においてなるべく中心に位置するメッシュであることが好ましい。例えば横３×縦３の９メッシュをグルーピングする場合、中心に位置するメッシュを代表メッシュとすることが好ましい。
［到達率推定装置の動作］
以下、到達率推定処理における到達率推定装置１ｄの動作について説明する。
図１３は、本発明の第４の実施形態に係る到達率推定装置１ｄの動作を示すフローチャートである。
平均受信電力算出部１１は、設計者によって指定されたエリアの地図データを微小なメッシュに分割する。ここで、メッシュ数をＭとする。平均受信電力算出部１１は、分割された各メッシュについて、各基地局における平均受信電力を算出する（ステップＳ３０１）。
到達率推定部１３ｄは、無線設備データベース１２から、各基地局のアンテナ数及び受信感度を示す情報を取得する（ステップＳ３０２）。
メッシュを識別する変数ｍに１が代入される（ステップＳ３０３）。
メッシュ＃ｍの到達率が推定済みであるか否か（すなわち、メッシュ＃ｍをグルーピング済みであるか否か）についての判定がなされる（ステップＳ３０４）。メッシュ＃ｍの到達率が推定済みである場合（ステップＳ３０４－Ｙｅｓ）、ステップＳ３０７に移行する。
一方、メッシュ＃ｍの到達率が推定済みでない場合（ステップＳ３０４－Ｎｏ）、グルーピング部１８は、メッシュ＃ｍのグループを生成する。そして、グルーピング部１８は、メッシュ＃ｍ及び当該グループに隣接する隣接メッシュのうち、まだグルーピングされていないメッシュについて、平均受信電力算出部１１によって算出された平均受信電力の値が高い順に基地局をソートし、上位所定数局の平均受信電力の値を取得する。
そして、グルーピング部１８は、隣接メッシュにおける同一の基地局の平均受信電力と、グループ内の各メッシュにおける同一の基地局の平均受信電力と、の差分が小さい場合に、当該隣接メッシュをグループに含める処理を繰り返すことによってグルーピングを行う（ステップＳ３０５）。なお、隣接メッシュにおける同一の基地局の平均受信電力と、グループ内の各メッシュにおける同一の基地局の平均受信電力と、の差分が小さい場合とは、例えば、対象とする各基地局の平均受信電力について、隣接メッシュにおける平均受信電力と、グループ内の各メッシュにおける平均受信電力のうち最大値及び最小値と、の差分がいずれも所定値以内である場合である。
到達率推定部１３ｄは、グルーピングされた複数のメッシュのうち、１つの代表メッシュについて、上述した第１の実施形態又は第２の実施形態と同様の方法を用いて、瞬時変動やダイバーシチ効果を考慮した厳密な到達率計算を行う（ステップＳ３０６）。
変数ｍがメッシュ数Ｍと等しいか否かについての判定がなされる（ステップＳ３０７）。変数ｍがメッシュ数Ｍと等しくない場合（ステップＳ３０７－Ｎｏ）、変数ｍに１が加算され（ステップＳ３０８）、ステップＳ３０４に移行する。一方、変数ｍがメッシュ数Ｍと等しい場合（ステップＳ３０７－Ｙｅｓ）、図１２のフローチャートが示す到達率推定装置１ｄの動作が終了する。
以上説明したように、第４の実施形態に係る到達率推定装置１ｄは、上位数局の、平均受信電力が近似する複数のメッシュをグルーピングする。そして、到達率推定装置１ｄは、グルーピングされた複数のメッシュのうち、１つの代表メッシュについてのみ、瞬時変動及びダイバーシチ効果を考慮した厳密な到達率の計算を行う。そして、到達率推定装置１ｄは、グルーピングされたその他のメッシュに対しては、上記代表メッシュについて同算出された計算結果を適用する。これにより、第４の実施形態に係る到達率推定装置１ｄは、到達率推定に要する計算時間をより短縮させることができる。
上述した実施形態における到達率推定装置１ａ～１ｄの一部又は全部を、コンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。更に「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、上述した機能の一部を実現するためのものであっても良く、更に上述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＰＬＤ（Programmable Logic Device）やＦＰＧＡ（Field Programmable Gate Array）等のハードウェアを用いて実現されるものであってもよい。
以上、図面を参照して本発明の実施形態を説明してきたが、上記実施形態は本発明の例示に過ぎず、本発明が上記実施形態に限定されるものではないことは明らかである。したがって、本発明の技術思想及び要旨を逸脱しない範囲で構成要素の追加、省略、置換、及びその他の変更を行ってもよい。
１ａ～１ｄ…到達率推定装置、１１…平均受信電力算出部、１２…無線設備データベース、１３ａ～１３ｄ…到達率推定部、１４…判定部、１５…厳密計算実施判定テーブル作成部、１６…厳密計算実施判定テーブル記憶部、１７…厳密計算実施判定部、１８…グルーピング部

Claims (8)

  1. 複数の基地局及び前記基地局が備える複数のアンテナの各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出部と、
    前記平均受信電力を、レイリーフェージングにより瞬時変動する振幅の確率密度関数から得られる受信電力瞬時値の累積分布関数に代入することで通信到達率を算出し、算出された通信到達率に基づいてサイトダイバーシチ及びアンテナダイバーシチ効果が考慮された通信到達率を算出する到達率推定部と、
    を備える到達率推定装置。
  2. 前記到達率推定部は、
    受信電力の瞬時変動の影響を解析式によって算出する
    請求項１に記載の到達率推定装置。
  3. 前記平均受信電力算出部は、
    上り通信の平均受信電力と下り通信の平均受信電力とを算出し、
    前記到達率推定部は、
    上り通信の到達率と下り通信の到達率の結合確率を算出する
    請求項１又は請求項２に記載の到達率推定装置。
  4. 繰り返し判定部
    をさらに備え、
    前記平均受信電力算出部は、
    算出された平均受信電力の値が高い順に前記基地局をソートし、
    前記到達率推定部は、
    考慮する基地局数を電力の高い順に徐々に増やしながら到達率を繰り返し算出し、
    前記繰り返し判定部は、
    算出された到達率が目標値以上になった場合、又は１つ前の繰り返し処理で算出された到達率からの変化が収束目標値以内になった場合に前記繰り返し処理を終了させる
    請求項１又は請求項２に記載の到達率推定装置。
  5. 厳密計算実施判定テーブル作成部と、
    厳密計算実施判定部と、
    をさらに備え、
    前記平均受信電力算出部は、
    算出された平均受信電力の値が高い順に前記基地局をソートし、
    前記厳密計算実施判定テーブル作成部は、
    平均受信電力と基地局数とから厳密な到達率計算の実施要否を判定するための厳密計算実施判定テーブルを作成し、
    前記厳密計算実施判定部は、
    前記厳密計算実施判定テーブルを参照して厳密計算の実施要否を判定する
    請求項１又は請求項２に記載の到達率推定装置。
  6. グルーピング部
    をさらに備え、
    前記平均受信電力算出部は、
    算出された平均受信電力の値が高い順に前記基地局をソートし、
    前記グルーピング部は、
    隣接メッシュと上位数基地局の平均受信電力の差分が小さい場合に前記メッシュのグルーピングを行い、
    前記到達率推定部は、
    グループ内の１つのメッシュについて瞬時変動の影響とダイバーシチ効果を考慮して到達率を算出し、算出された前記到達率を前記グループ内の他のメッシュについての到達率と見なす
    請求項１又は請求項２に記載の到達率推定装置。
  7. 複数の基地局及び前記基地局が備える複数のアンテナの各々に関して、エリアを構成するメッシュごとに平均受信電力を算出する平均受信電力算出ステップと、
    前記平均受信電力を、レイリーフェージングにより瞬時変動する振幅の確率密度関数から得られる受信電力瞬時値の累積分布関数に代入することで通信到達率を算出し、算出された通信到達率に基づいてサイトダイバーシチ及びアンテナダイバーシチ効果が考慮された通信到達率を算出する到達率推定ステップと、
    を有する到達率推定方法。
  8. 請求項１から請求項６のうちいずれか一項に記載の到達率推定装置としてコンピュータを機能させるためのプログラム。
JP2019044959A 2019-03-12 2019-03-12 到達率推定装置、到達率推定方法、及びプログラム Active JP7174249B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2019044959A JP7174249B2 (ja) 2019-03-12 2019-03-12 到達率推定装置、到達率推定方法、及びプログラム
PCT/JP2020/008678 WO2020184253A1 (ja) 2019-03-12 2020-03-02 到達率推定装置、到達率推定方法、及びプログラム
US17/437,196 US11991545B2 (en) 2019-03-12 2020-03-02 Arrival rate estimation apparatus, arrival rate estimation method and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019044959A JP7174249B2 (ja) 2019-03-12 2019-03-12 到達率推定装置、到達率推定方法、及びプログラム

Publications (2)

Publication Number Publication Date
JP2020150366A JP2020150366A (ja) 2020-09-17
JP7174249B2 true JP7174249B2 (ja) 2022-11-17

Family

ID=72427931

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019044959A Active JP7174249B2 (ja) 2019-03-12 2019-03-12 到達率推定装置、到達率推定方法、及びプログラム

Country Status (3)

Country Link
US (1) US11991545B2 (ja)
JP (1) JP7174249B2 (ja)
WO (1) WO2020184253A1 (ja)

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL131700A0 (en) * 1999-03-08 2001-03-19 Mintz Yosef Method and system for mapping traffic congestion
JP3439396B2 (ja) * 1999-09-24 2003-08-25 日本電信電話株式会社 電磁環境設計方法および設計プログラムを記録した記録媒体
US7522657B2 (en) * 2003-10-20 2009-04-21 William Marsh Rice University Throughput maximization in wireless communication systems
US7373168B1 (en) * 2005-01-12 2008-05-13 The Aerospace Corporation Power controlled fading communication channel system
JP2007259445A (ja) * 2006-03-20 2007-10-04 Fujitsu Ltd Ｏｆｄｍ通信システムにおける送信装置及び方法
JP5228951B2 (ja) * 2009-01-30 2013-07-03 沖電気工業株式会社 パケット中継システム及び無線ノード
EP2723133A1 (en) * 2012-10-18 2014-04-23 NTT DoCoMo, Inc. Method for controlling operation within a cell of a wireless cellular network, base station and wireless cellular network
JP6097409B2 (ja) * 2013-01-16 2017-03-15 エヌイーシー（チャイナ）カンパニー， リミテッドＮＥＣ（Ｃｈｉｎａ）Ｃｏ．，Ｌｔｄ． Ｔｄｄシステムにおけるｄｌ／ｕｌリソース設定方法及び装置
US10848355B2 (en) * 2014-06-03 2020-11-24 Valorbec Societe En Commandite Methods and systems for cognitive radio spectrum monitoring
US11606667B2 (en) * 2017-03-17 2023-03-14 WangsLabs, LLC Precise positioning system and method of using the same
JP6955144B2 (ja) * 2017-06-29 2021-10-27 富士通株式会社 無線機器の設置位置決定装置、無線機器の設置位置決定方法及び無線機器の設置位置決定プログラム
TWI666947B (zh) * 2018-05-16 2019-07-21 Yuan Ze University 多用戶隨機存取訊號之分析方法

Also Published As

Publication number Publication date
US11991545B2 (en) 2024-05-21
JP2020150366A (ja) 2020-09-17
US20220174520A1 (en) 2022-06-02
WO2020184253A1 (ja) 2020-09-17

Similar Documents

Publication Publication Date Title
Wang et al. MmWave vehicle-to-infrastructure communication: Analysis of urban microcellular networks
Weiler et al. Quasi-deterministic millimeter-wave channel models in MiWEBA
Karttunen et al. Spatially consistent street-by-street path loss model for 28-GHz channels in micro cell urban environments
Bor-Yaliniz et al. Environment-aware drone-base-station placements in modern metropolitans
JP6696859B2 (ja) 品質推定装置及び品質推定方法
Filiposka et al. Terrain-aware three-dimensional radio-propagation model extension for NS-2
Kim et al. Geometric optics-based propagation prediction model in urban street canyon environments
Wang et al. An extension of spatial channel model with spatial consistency
Eisenblätter et al. UMTS radio network evaluation and optimization beyond snapshots
Fernandes et al. Comparison of artificial intelligence and semi-empirical methodologies for estimation of coverage in mobile networks
JP7161136B2 (ja) 置局設計装置、置局設計方法及びプログラム
JP7174249B2 (ja) 到達率推定装置、到達率推定方法、及びプログラム
Andrade et al. Radio channel spatial propagation model for mobile 3G in smart antenna systems
CN112929907B (zh) 天线参数的确定方法及装置
JP7096499B2 (ja) 無線端末収容判定装置、無線端末収容判定方法、及びプログラム
Xue et al. A probabilistic propagation graph modeling method for channel characterization in dense urban environments
JP2004304302A (ja) 伝搬環境模擬装置およびプログラム
Yang et al. Coverage analysis of heterogeneous cellular networks in urban areas
Comisso et al. Shannon capacity evaluation for 5G communications using the 3D random waypoint mobility model
Kim et al. Achievable rate estimation of IEEE 802.11 ad visual big-data uplink access in cloud-enabled surveillance applications
Meesa-ard et al. Analyzing the impact of mobility over α-μ generalized fading channels in wireless communication
JP3854973B2 (ja) フェージングシミュレータ
JP3817558B2 (ja) フェージングシミュレータ
WO2020008181A1 (en) Evaluating the wireless performance of a building
JP7280654B1 (ja) 制御装置、コードブック生成装置、制御方法、コードブック生成方法及びプログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210702

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20221004

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221017

R150 Certificate of patent or registration of utility model

Ref document number: 7174249

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150