JP6239670B2 - ドップラースペクトルを用いた端末の移動速度推定方法 - Google Patents

ドップラースペクトルを用いた端末の移動速度推定方法 Download PDF

Info

Publication number
JP6239670B2
JP6239670B2 JP2016060664A JP2016060664A JP6239670B2 JP 6239670 B2 JP6239670 B2 JP 6239670B2 JP 2016060664 A JP2016060664 A JP 2016060664A JP 2016060664 A JP2016060664 A JP 2016060664A JP 6239670 B2 JP6239670 B2 JP 6239670B2
Authority
JP
Japan
Prior art keywords
doppler
spectrum
frequency
width
threshold
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
JP2016060664A
Other languages
English (en)
Other versions
JP2017173171A (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.)
SoftBank Corp
Original Assignee
SoftBank 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 SoftBank Corp filed Critical SoftBank Corp
Priority to JP2016060664A priority Critical patent/JP6239670B2/ja
Priority to PCT/JP2017/010027 priority patent/WO2017163988A1/ja
Publication of JP2017173171A publication Critical patent/JP2017173171A/ja
Application granted granted Critical
Publication of JP6239670B2 publication Critical patent/JP6239670B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S11/00Systems for determining distance or velocity not using reflection or reradiation
    • G01S11/02Systems for determining distance or velocity not using reflection or reradiation using radio waves
    • G01S11/10Systems for determining distance or velocity not using reflection or reradiation using radio waves using Doppler effect

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

本発明は、無線通信システムにおけるドップラースペクトルを用いた端末の移動速度推定方法に関し、特に、複数のしきい値を用いた端末の移動速度推定方法に関する。
無線通信システムは、ユーザによる無線通信の多様な利用を想定し、低データレート通信から高データレート通信、静止画から動画ストリーミング伝送、静止状態や歩行状態での利用から高速移動状態での利用、あるいは屋内での使用から屋外での使用など様々な形態による無線通信サービスを提供している。
無線通信システムを提供する事業者は、このような利用形態に加え、都市部や農村部などサービスを提供する場所、サービスを利用するユーザの人口分布やデータ量(トラフィック密度)など、様々な条件を考慮したシステム構築を行っている。システム構築の一環として、事業者は、どのような場所や時間においても、ユーザが遅滞なく無線通信システムに接続することができ、満足するサービスの提供を受けられるように、利用形態やトラフィック密度を考慮して基地局(またはセル)の配置設計を行う。
しかし、利用形態の更なる多様化やトラフィック量の増加に伴い、平面的なセル配置では、十分なサービスを提供することは困難となる。このため、複数の種類のセルを組み合わせることにより、多様な利用形態やトラフィック特性に合致したサービスを提供するオーバーレイセル構成の構築が進められている。オーバーレイセル構成の1つの例において、セルサイズの大きいマクロセル内にセルサイズの小さい極小セルを重畳(オーバーレイ)させる。極小セルは、マイクロセル、ピコセル、フェムトセルなどである。このような極小セルは、マクロセル内のトラフィック密度が高い場所、高データレート通信が想定される場所、マクロセル基地局からの電波受信が困難な場所等に設置される。
オーバーレイセル構成において、周波数利用効率を向上させるため、あるいは、ハンドオーバの頻度を抑制するために、端末の移動速度に応じて、マクロセルと極小セルの階層を適切に選択する階層選択制御が必要となる。
移動速度が速い端末は、マクロセルのみに接続し、極小セルには接続されないようにする。一方、移動速度の遅い端末は、極小セル内にある場合は極小セルと接続し、極小セルが重畳されていない場所では、マクロセルと接続する。このような階層選択制御により、高速移動中の端末が極小セルに頻繁にハンドオーバすることを防止して、ハンドオーバの制御シグナリングの増加やハンドオーバ増加に伴う通信断を減らすことができる。また、低速移動の端末が極小セルに接続することで、送信電力が抑制されて、システム内の干渉電力を減らすことができる。
階層選択制御を適切に実行するためには、端末の移動速度を正確に測定し、その移動速度に基づいて、マクロセル又は極小セルに接続させることが重要である。
端末の移動速度を推定する方法の1つとして、ドップラースペクトルを利用する方法がある(非特許文献1)。ここでドップラースペクトルとは、伝搬路応答の時間変動をフーリエ変換した結果得られる周波数軸(ドップラー周波数シフト)上の電力波形である。端末と基地局間の無線伝搬環境が電波の到来方向分布が一様である全周散乱モデルと見なせる場合、ドップラースペクトルは、最大ドップラー周波数においてピークを示す。このため、スペクトルのピークを検出することにより、最大ドップラー周波数を測定することができ、最大ドップラー周波数から端末の移動速度を算出することができる。
G-H. Park, D. Hong and C-E. Kang, 「A New Doppler Spread Estimation Using FFT」 IEICE Trans. Commun., vol.E86-B, No.9, pp.2799-2803, Sept. 2003. マイティ、星野、三上、「最大ドップラー周波数検出に基づく端末移動速度推定アルゴリズムの見通し外伝搬環境における屋外実験評価」、2015年電子情報通信学会ソサイエティ大会、2015年9月. マイティ、星野、三上、「最大ドップラー周波数検出に基づく端末高速・低速移動判定アルゴリズムの屋外実験評価」、2016年電子情報通信学会総合大会、2016年3月.
端末が自動車内にあり、その端末搭載自動車が一般道路や高速道路を走行すると、対向車や同方向車などの自動車が端末搭載自動車の近傍を走行する。このような対向車や同方向車は、端末から送信される電波、あるいは端末で受信される電波を散乱させる。このような自動車を周辺移動散乱体とする。周辺移動散乱体の存在により、端末のドップラースペクトルの算出において、端末単体によるドップラースペクトルとは異なるドップラースペクトル成分が発生する。この異なるドップラースペクトル成分は、端末単体によるドップラースペクトルより低い周波数成分または高いスペクトル周波数成分を有している。また、通常は、端末単体によるドップラースペクトルの電力より小さい電力となる。対向車による速度推定値への影響が確認されている文献として、非特許文献2、非特許文献3を挙げる。
端末のドップラースペクトルを算出する際、周辺移動散乱体によるドップラースペクトル成分を有効成分として検出してしまうと、端末単体によるドップラースペクトルより広いスペクトルとなり、最大ドップラー周波数が高くなってしまう。このため、周辺移動散乱体が存在すると、端末の移動速度が実際より早い速度として推定される可能性がある。
端末のドップラースペクトルを算出し、算出されたドップラースペクトルから端末の移動速度を推定する方法において、周辺移動散乱体によるドップラースペクトル成分を検知して除外することで、正確な端末のドップラースペクトルを算出し、より実際の移動速度を反映した端末の移動速度を推定することができる。
本発明は、ドップラースペクトルを用いた端末の移動速度推定において、周辺移動散乱体によるドップラースペクトル成分が存在する場合においても、周辺移動散乱体によるドップラースペクトル成分を検知して除外し、正確な移動速度推定を行う移動速度推定方法を提供することを目的とする。
本発明は、端末の移動速度を推定する方法であって、参照信号を受信するステップと、前記受信した参照信号から、伝搬路の時間変動データを作成するステップと、前記伝搬路の時間変動データにFFTを掛けるステップと、前記FFTの出力からドップラースペクトルを算出するステップと、前記ドップラースペクトルの成分のうち、最も電力が高い成分の電力値をピーク電力と設定するステップと、雑音電力を算出するステップと、前記ピーク電力と前記雑音電力の間に、等間隔で所定個数のしきい値を設定するステップと、前記所定個数のしきい値の各々において、しきい値以上となる前記ドップラースペクトルの成分のうち、最も高い周波数の成分と最も低い周波数の成分の周波数の差分から、スペクトル幅を算出するステップと、前記所定個数のしきい値において前記算出されたスペクトル幅のうち、最も高いしきい値のスペクトル幅を第1のドップラースペクトル幅候補とするステップと、前記最も高いしきい値より低いしきい値のスペクトル幅の各々と、前記第1のドップラースペクトル幅候補との差分を算出し、前記算出された差分が周波数差のしきい値以下となるスペクトル幅を第2のドップラースペクトル幅候補とするステップと、前記第2のドップラースペクトル幅候補がない場合、前記第1のドップラースペクトル幅候補をドップラースペクトル幅と決定するステップと、前記第2のドップラースペクトル幅候補がある場合、最も低いしきい値の前記第2のドップラースペクトル幅候補をドップラースペクトル幅と決定するステップと、前記決定されたドップラースペクトル幅から最大ドップラー周波数を算出するステップと、前記最大ドップラー周波数から端末の移動速度を決定するステップと、を有する。
本発明に係る移動速度推定方法において、周辺移動散乱体によるドップラースペクトル成分を検出し、ドップラースペクトル算出から除外するため、端末の移動速度を正確に推定することができる。
ドップラースペクトルの例を示す図である。 最大ドップラー周波数の検出方法を説明する図である。 ドップラースペクトルへのしきい値の適用を説明する図である。 最大ドップラー周波数の推定方法のフローチャートを示す図である。 周辺移動散乱体によるドップラースペクトル成分を示す図である。 端末搭載車と対向車との関係を示す図である。 周辺移動散乱体によるドップラースペクトル成分の変動を示す図である。 複数のしきい値を用いた最大ドップラー周波数の検出方法の実施例を説明する図である。 複数のしきい値を用いた最大ドップラー周波数の検出方法のフローチャートを示す図である。 複数のしきい値を用いた最大ドップラー周波数のその他の検出方法のフローチャートを示す図である。 一部のピーク電力が突出したドップラースペクトルに対する複数しきい値の設定を示す図である。 突出したピーク電力を考慮した複数しきい値の設定を示す図である。 本発明による装置構成の例を示す図である。
以下図面を参照して、ドップラースペクトルを用いた端末の速度推定方法について説明する。しかしながら、本発明が、図面又は以下に記載される実施形態に限定されるものではないことを理解されたい。
図1にドップラースペクトルの例を示す。この例は、無線伝搬環境が全周散乱モデルとみなせる理想的な場合を表している。図において、横軸は周波数を表し、縦軸は相対電力を表す。ただし、図1ではベースバンド帯域に変換した後の処理結果を示しているため、実際には横軸の周波数0 Hzがキャリア周波数に対応することに注意が必要である。以降のドップラースペクトルの図においても横軸の周波数の定義は図1と同様である。この例において、最大ドップラー周波数は、100 Hzであり、ドップラースペクトルの周波数100 Hz及び-100 Hzにおいて、高い電力の成分が生じている。このように、受信信号の伝搬路変動を用いて算出されたドップラースペクトルから、電力の高い成分の周波数を検出することで、最大ドップラー周波数を推定することができる。
図2は、ドップラースペクトルから最大ドップラー周波数を算出する方法を簡単に説明した図である。周波数の負側において、最大電力が得られる周波数をf(1)とし、正側において、最大電力が得られる周波数をf(2)とする。理想的なフェージング環境下においては、f(1)とf(2)の間でU字型のスペクトルを示し、f(1)とf(2)は、同じ周波数絶対値を示す。この周波数絶対値が最大ドップラー周波数fDとなる。即ち、f(1)=-fDであり、f(2)=fDとなる。そして、f(1)とf(2)の間の周波数がドップラースペクトル幅となる。
最大ドップラー周波数fDが得られると、端末の移動速度νは、fD×λで算出することができる。ここで、λは、キャリア周波数の波長である。
ドップラースペクトルから最大ドップラー周波数を検出する方法が幾つか考えられる。以下に、ドップラースペクトルに対して固定のしきい値を設定し、しきい値以上となる有効成分について、そのうちの最大周波数の成分と最小周波数の成分から最大ドップラー周波数を算出する方法を説明する。
図3は、ドップラースペクトルにおけるしきい値の設定の一例を示した図である。ここでは、雑音電力N0 [dBm]を算出し、更に、ΔN [dB]を加算した値をしきい値とする。図3に示されるように、N0+ΔN [dBm]をしきい値とすることで、ドップラースペクトルの外側の雑音成分を除外することができる。そして、しきい値以上となる有効成分は、直接波成分などの幾つかの電力が高い成分を含む。これらの有効成分の中から、最大の周波数を持つ成分(即ち、図3の「A」で示される100 Hzの成分)と最小周波数を持つ成分(即ち、図3の「B」で示される-100 Hzの成分)を選択することで、最大ドップラー周波数を決定することができる。
図4は、最大ドップラー周波数算出の方法を示すフローチャートである。S10において、無線接続している基地局または端末から参照信号を受信する。ここで、参照信号としたが、システムにより、パイロット信号やビーコンなどとも呼ばれ、伝搬路推定に使用される既知の信号である。S11において、受信した参照信号をベースバンド信号に変換し、復調した後、受信電力を測定することにより、受信信号の伝搬路の時間変動データを作成する。伝搬路の時間変動データは、一定区間毎にバッファに格納される。
受信信号は、端末の移動速度に比例して変動するため、受信信号の伝搬路変動をフーリエ変換することにより、ドップラースペクトルが得られる。S12において、一定区間毎に格納された伝搬路の時間変動データに高速フーリエ変換(FFT)を掛ける。ここで、一定区間は、FFTのポイント数に相当する。S13において、FFTの出力からドップラースペクトルが算出される。
S14において、スペクトルのしきい値を設定し、しきい値以上となる有効成分を検出する。S15において、有効成分の中から最大周波数を持つ有効成分Ef_maxを選択し、その周波数f(Ef_max)を検出する。S16において、有効成分の中から最小周波数を持つ有効成分Ef_minを選択し、その周波数f(Ef_min)を検出する。
S17において、検出された2つの周波数を用いて、最大ドップラー周波数を算出する。ここでは、2つの周波数f(Ef_max)と周波数f(Ef_min)からドップラースペクトル幅を算出し、その半値を最大ドップラー周波数fDとすることができる。即ち、fD=(f(Ef_max)−f(Ef_min))/2とする。
このように、しきい値を設定し、しきい値以上となる有効成分を抽出することで、雑音成分によるスペクトルを除外して、受信信号によるスペクトルのみを選択することができる。さらに、適切なしきい値を設けることで、雑音成分による誤検出を抑制し、より正確な最大ドップラー周波数を推定することができる。
上記例において、しきい値を雑音電力N0にΔNを加算した値に設定したが、これに限定されるものではない。例えば、測定されるSNRなどの無線伝搬状況に応じて、しきい値を可変とすることも可能である。
また、周波数f(Ef_max)と周波数f(Ef_min)とのうち、絶対値の大きい方を最大ドップラー周波数fDとする方法もある。即ち、fD=max(|f(Ef_max)|, |f(Ef_min)|)とする。
次に、周辺移動散乱体によるドップラースペクトル算出への影響とその除外方法について説明する。端末搭載自動車が、一般道路や高速道路を走行する際、対向車や同方向車が存在する。このように端末搭載自動車周辺を移動する車両体(周辺移動散乱体)は、端末から出力される電波、または端末で受信される電波を散乱させる。
周辺移動散乱体が存在すると、端末のドップラースペクトル算出において、端末単体によるドップラースペクトルとは別に、周辺移動散乱体によるドップラースペクトル成分が発生する。簡単な例を図5に示す。端末単体によるドップラースペクトルは、端末の移動速度に応じた最大ドップラー周波数f0と-f0により示される。この端末単体によるドップラースペクトルより高い周波数位置(f1)において、周辺移動散乱体によるドップラースペクトル成分が発生している。このドップラースペクトル成分の電力は、通常、端末単体によるドップラースペクトルの電力より小さくなる。
しかし、周辺移動散乱体によるドップラースペクトル成分を端末のドップラースペクトルの有効成分として検出してしまうと、ドップラースペクトル幅は、-f0〜f1となり、本来のドップラースペクトル幅より広く検出される。図4に示すフローチャートに従い、最大ドップラー周波数を算出すると、ドップラースペクトル幅の半値は(f1+f0)/2となり、実際の最大ドップラー周波数より大きくなる。また、最大周波数と最小周波数の絶対値の大きい方から最大ドップラー周波数を算出しても最大ドップラー周波数はf1となる。
図6に端末が搭載される自動車Aと対向車Bの関係を簡単に示す。対向車Bが周辺移動散乱体となる。自動車Aは速度v0で走行し、その最大ドップラー周波数はf0とする。一方、対向車Bは速度vで走行し、その最大ドップラー周波数はfとする。また、時刻tにおける自動車Aと対向車Bとがなす角度をθ(t)とすると、自動車Aと対向車Bの相対速度v1は、v1=v0+2v×cosθ(t)となる。このため、対向車Bによるドップラースペクトル成分の周波数f1は、f1=f0+2f×cosθ(t)と表せる。
対向車Bが自動車Aに近づくと、角度θ(t)が90°に近づくため、対向車Bによるドップラースペクトル成分の周波数f1は、f0に近づく。即ち、図7に示されるように、対向車Bによるドップラースペクトル成分は、自動車A単体によるドップラースペクトルより高い周波数位置に現れ、両者が接近するに従い、自動車A単体によるドップラースペクトルの中心に近づく。
対向車Bが自動車Aの横を通り、離れると、対向車Bによるドップラースペクトル成分の周波数f1は、f1= -f0-2f×cosθ(t)となり、自動車A単体によるドップラースペクトルより低い周波数位置に移動する。
周辺移動散乱体が同方向車の場合も同様に考えることができる。同方向車の速度が自動車Aの速度より速く、自動車Aを追い抜いて行く場合、同方向車によるドップラースペクトル成分は、自動車A単体によるドップラースペクトルより高い周波数位置で発生し、追い抜いた後、低い周波数位置に移動する。
更に、自動車Aの近傍に対向車や同方向車が複数存在する場合、自動車A単体によるドップラースペクトルの片外側または両外側に複数のドップラースペクトル成分が出現することも考えられる。以下に示す本発明による速度推定方法は、このように複数のドップラースペクトル成分が存在している場合においても正確に速度推定を行うことができる。
上記のように周辺移動散乱体によるドップラースペクトル成分が存在すると、端末のドップラースペクトルの正確な検出を妨げることになる。周辺移動散乱体によるドップラースペクトル成分を検出し、端末のドップラースペクトルと区別して、除外することで、端末単体によるドップラースペクトルのみを算出し、正確な速度推定を行うことができる。
本発明の1つの実施例における周辺移動散乱体によるドップラースペクトル成分の検出方法を説明する。図4に示されるフローチャートのS10〜S13の処理により、ドップラースペクトルを算出する。算出したスペクトルの各周波数成分の電力を検査し、最大電力を示すスペクトルピークPpeakを求める。同時に、雑音電力N0を算出する。
スペクトルピークPpeakと雑音電力N0との間に、M個のしきい値を等間隔で設定する。各しきい値の電力をΔNm(m=1〜M)とする。しきい値の間隔は、T=(Ppeak-N0)/(M+1)となり、各しきい値の電力は、ΔNm=N0+T×mで表せる。
各しきい値におけるスペクトル幅ΔfD(m)を求める。スペクトル幅ΔfD(m)は、しきい値ΔNm以上となる周波数成分を抽出し、その中で、最も周波数が低い周波数成分と最も周波数が高い周波数成分との周波数の差分から算出する。
各しきい値におけるスペクトル幅ΔfD(m)を比較する。周辺移動散乱体によるドップラースペクトル成分が1つ存在するとする。最も低いしきい値からしきい値を上げていくと、スペクトル幅が大きく減少するしきい値が出てくる。このとき、そのしきい値より下のしきい値では、周辺移動散乱体によるドップラースペクトル成分が存在して、スペクトル幅が広くなり、そのしきい値以上のしきい値では、周辺移動散乱体によるドップラースペクトル成分がなく、端末単体によるドップラースペクトルのみが存在することになる。
図8にドップラースペクトルの算出方法の1例を示す。端末単体によるドップラースペクトルが周波数-f0〜f0で示され、周辺移動散乱体によるドップラースペクトル成分が周波数f1で示されている。各周波数成分における電力を検査し、最大電力を示すスペクトルピークPpeakを求める。図8の例では、端末単体によるドップラースペクトルの最大周波数及び最小周波数位置において、ピーク電力が生じている。また、雑音電力N0を算出する。
スペクトルピークPpeakと雑音電力N0との間に、M個のしきい値を等間隔で設定する。各しきい値をΔNm(m=1〜M)とする。図8の例では、M=4である。そして、各しきい値におけるスペクトル幅ΔfD(1)〜ΔfD(M)を算出する。ここで、次を仮定する。
仮定:最も高いしきい値(m=M)において、スペクトル幅ΔfD(M)は、端末単体によるドップラースペクトル幅を示す。即ち、周辺移動散乱体によるドップラースペクトル成分の電力は、端末単体によるドップラースペクトルの電力より小さく、最も高いしきい値まで到達しない。
高速フーリエ変換(FFT)を用いてドップラースペクトルを算出する際、有限なサンプル数による演算の影響により、算出されたスペクトルには歪が生じてしまう。そのため、各しきい値におけるスペクトル幅は、1つ上のしきい値のスペクトル幅と同じか、または広くなる。上記の仮定より、最も高いしきい値におけるスペクトル幅ΔfD(M)は、端末単体によるドップラースペクトル幅を表し、各しきい値の中で最も狭いスペクトル幅となる。簡単には、最も高いしきい値におけるスペクトル幅を端末のドップラースペクトル幅と設定することができる。以下に、より正確に端末のドップラースペクトル幅を算出する方法を説明する。
周辺移動散乱体によるドップラースペクトル成分を検出して除外し、端末のドップラースペクトル幅を算出する方法の1つの実施例を説明する。
上で仮定したように、スペクトル幅ΔfD(M)は、端末単体によるドップラースペクトル幅を示す。このため、しきい値ΔNMより低いしきい値におけるスペクトル幅がΔfD(M)より所定値を超えて広くなる場合、端末単体によるドップラースペクトルの外側に別のドップラースペクトル成分が存在することが想定される。別のドップラースペクトル成分の存在が想定された場合、そのドップラースペクトル成分が検知されたしきい値のスペクトル幅をドップラースペクトル幅の算出の候補から除外する。
まず、スペクトル幅ΔfD(M)と1つ下のしきい値のスペクトル幅ΔfD(M-1)を比較し、その差分が周波数差のしきい値fthreshを超えるか否かを判定する。差分が周波数差のしきい値を超えない場合、即ち、ΔfD(M-1)-ΔfD(M)<=fthreshの場合、しきい値ΔNM-1において、周辺移動散乱体によるドップラースペクトル成分は存在せず、端末単体によるドップラースペクトルのみ存在すると判断する。
一方、差分が周波数差のしきい値を超える場合、即ち、ΔfD(M-1)-ΔfD(M)>fthreshの場合、しきい値ΔNM-1において、端末単体によるドップラースペクトルの外側に周辺移動散乱体によるドップラースペクトル成分が存在すると判断し、m=M-1のスペクトル幅ΔfD(M-1)をドップラースペクトル幅の算出の候補から除外する。この場合、m<M-1のスペクトル幅もΔfD(M-1)以上となるため、候補から除外される。この結果、ΔfD(M)がドップラースペクトル幅として決定され、速度推定に使用される。
高速フーリエ変換(FFT)を用いてドップラースペクトルを算出する際、算出されたスペクトルには歪が生じてしまう。このような歪はスペクトル幅の10%程度と考えられる。このため、ドップラースペクトルからスペクトル幅の10%より離れた周波数位置に別のドップラースペクトル成分が存在すると、検知することができるようになる。従って、周波数差のしきい値fthreshは、想定されるドップラースペクトル幅の10%程度に設定することができる。ただし、これに制限されるものではない。
スペクトル幅ΔfD(M-1)とΔfD(M) との差分が周波数差のしきい値を超えない場合、続いて、スペクトル幅ΔfD(M-2)とΔfD(M) との比較を行う。差分が周波数差のしきい値を超える場合、即ち、ΔfD(M-2)-ΔfD(M)>fthreshの場合、m=M-2のスペクトル幅をドップラースペクトル幅の算出の候補から除外する。この場合、m=MとM-1のスペクトル幅ΔfD(M)とΔfD(M-1)が候補となるが、より広い幅を有する可能性があるΔfD(M-1) をドップラースペクトル幅として決定し、速度推定に使用する。
上記処理を、スペクトル幅ΔfD(M) との差分が周波数差のしきい値を超えるまで、あるいは、m=1まで行う。m=1において、スペクトル幅の差分が周波数差のしきい値を超えない場合、全てのしきい値において周辺移動散乱体によるドップラースペクトル成分が存在しないと判断される。
図8の例では、最も高いしきい値のスペクトル幅ΔfD(4)をドップラースペクトル幅算出の候補とする。次に、1つ下のしきい値のスペクトル幅と比較を行う。ΔfD(3)-ΔfD(4)<=fthreshであるため、m=3において、周辺移動散乱体によるドップラースペクトル成分が存在しないと判断し、ΔfD(3) をドップラースペクトル幅算出の新たな候補とする。更に、1つ下のしきい値のスペクトル幅と比較を行い、ΔfD(2)-ΔfD(4)>fthreshであるため、m=2において、周辺移動散乱体によるドップラースペクトル成分が存在すると判断する。この結果、ΔfD(3) をドップラースペクトル幅として決定し、速度推定に使用する。
図9に上記処理のフローチャートを示す。これは、図4に示されるフローチャートのS14〜S17に相当する。
図9のS21において、雑音電力N0とスペクトルピークPpeakを検出する。S22において、スペクトルピークと雑音電力との間にM個のしきい値を等間隔で設定する。しきい値をΔNm(m=1〜M)とする。S23において、各しきい値におけるスペクトル幅ΔfD(m)(m=1〜M)を計算する。スペクトル幅は、各しきい値において検出される周波数成分のうち、最大周波数成分と最小周波数成分の周波数差に相当する。
S24において、m=M-1と設定する。即ち、最も高いしきい値の1つ下のしきい値とする。S25において、スペクトル幅ΔfD(m)とΔfD(M) との比較を行う。差分が周波数差のしきい値fthreshを超える場合、m以下のしきい値のスペクトル幅を候補から除外する。S28において、m+1のスペクトル幅ΔfD(m+1)がドップラースペクトル幅として決定される。
S25において、スペクトル幅ΔfD(m)とΔfD(M) との差分が周波数差のしきい値fthreshを超えない場合、S26において、mを1つ減じて、S25の処理を行う。差分が周波数差のしきい値を超えない場合、S25とS26の処理を繰り返す。
最も低いしきい値(m=1)のスペクトル幅の比較において、周波数差のしきい値を超えない場合、S27において、m=0となり、S28において、m=1のスペクトル幅ΔfD(1)がドップラースペクトル幅として決定される。
図9のS25において、スペクトル幅ΔfD(m)とΔfD(M) との差分が周波数差のしきい値fthreshを超える場合、m以下のしきい値のスペクトル幅は、ドップラースペクトル幅算出の候補から除外される。そして、m+1のしきい値のスペクトル幅ΔfD(m+1)をドップラースペクトル幅として決定している。更に、m+1以上のしきい値のスペクトル幅を候補として、その平均値をドップラースペクトル幅として決定することもできる。
図3で説明したドップラースペクトルの検出方法では、雑音電力にΔN [dB]を加算した値を固定しきい値としてドップラースペクトルの有効成分の検出を行う。しかし、端末における受信信号のSNRが悪い環境では、ドップラースペクトルの電力が低くなり、ドップラースペクトルのピーク電力が固定しきい値より低くなり、有効成分の検出が出来ないことが考えられる。
これに対し、本発明による有効成分の検出では、ドップラースペクトルのピーク電力を見つけ、ピーク電力と雑音電力の間に複数のしきい値を設定している。このように、ピーク電力を検出し、ピーク電力を考慮したしきい値を設定することにより、受信SNRが低い環境においても、ドップラースペクトルの検出が可能となる。
続いて、周辺移動散乱体によるドップラースペクトル成分を検出して除外し、端末のドップラースペクトル幅を算出する方法の別の実施例を説明する。ここでは、しきい値の低い側から処理を行う。
まず、m=1として、最も低いしきい値のスペクトル幅ΔfD(1)とm+1である1つ上のしきい値のスペクトル幅ΔfD(2)の差分(ΔfD(1)-ΔfD(2))を計算する。差分が周波数差のしきい値fthreshを超える場合、即ち、ΔfD(1)-ΔfD(2)> fthreshの場合、m=1のしきい値とm=2のしきい値で、異なるスペクトル形状が検知されたことを意味する。この場合、m=1のしきい値において、端末単体のドップラースペクトルの外側に、異なるスペクトル成分が存在することになる。このとき、m=2のしきい値においては、端末単体のドップラースペクトルの外側に別のスペクトル成分が検出されなかったか、または、m=1のしきい値で検出された別のスペクトル成分の内側に更に別のスペクトル成分が存在することを示す。このため、m=1のしきい値において周辺移動散乱体によるドップラースペクトル成分があると判断し、スペクトル幅ΔfD(1)をドップラースペクトル幅の算出の候補から除外する。
差分ΔfD(1)-ΔfD(2)が周波数差のしきい値fthreshを超えない場合、更に1つ上の段のしきい値のスペクトル幅ΔfD(3)との差分を計算し、周波数差のしきい値fthreshと比較する。ここで、周波数差がしきい値fthreshを超える場合、m=1のスペクトル幅ΔfD(1)をドップラースペクトル幅の算出候補から除外する。一方、周波数差のしきい値fthreshを超えない場合には、処理を繰り返し、最も高いしきい値のスペクトル幅ΔfD(M)との差分を計算し、周波数差のしきい値fthreshと比較する。ここで、周波数差のしきい値fthreshを超える場合、m=1のスペクトル幅ΔfD(1)をドップラースペクトル幅の算出候補から除外する。スペクトル幅ΔfD(M)との比較において、周波数差のしきい値fthreshを超えない場合、全てのしきい値において、端末単体のドップラースペクトルの外側のスペクトル成分が検出されなかったことになり、スペクトル幅ΔfD(1)〜ΔfD(M)の全てが端末単体のドップラースペクトルを表すことになる。このため、最も広いスペクトル幅となる可能性があるスペクトル幅ΔfD(1)を端末のドップラースペクトル幅と決定する。
上記の処理において、m=1のスペクトル幅ΔfD(1)が除外された場合、m=1の処理を終了し、m=2として、同様の処理を行う。即ち、m=2のスペクトル幅ΔfD(2)とm+1のスペクトル幅ΔfD (3)の差分ΔfD(2)-ΔfD(3)を計算する。差分が周波数差のしきい値fthreshを超える場合、m=2のスペクトル幅ΔfD(2)をドップラースペクトル幅の算出候補から除外する。周波数差のしきい値fthreshを超えない場合、ΔfD(3)〜ΔfD(M)との比較を行う。同様の処理をm=3においても行う。
図8の例では、差分ΔfD(1)-ΔfD(2)が周波数差のしきい値fthreshを超えないため、ΔN1とΔN2とにおいて、同じスペクトル幅を持つと判断する。続いて、ΔfD(1)とΔfD(3)との差分を計算し、周波数差のしきい値fthreshと比較する。差分ΔfD(1)-ΔfD(3)が周波数差のしきい値fthreshを超えるため、ΔN1とΔN3とにおいて、異なるスペクトル幅を持つと判断する。即ち、ΔN1(そして、ΔN2)において、端末のドップラースペクトルの外側に別のスペクトル成分が発生していると判断する。このため、m=1のスペクトル幅ΔfD(1)をドップラースペクトル幅算出の候補から除外する。
次に、m=2として、同様の処理を行う。差分ΔfD(2)-ΔfD(3)と周波数差のしきい値fthreshとの比較により、差分が周波数差のしきい値fthreshを超えるため、ΔN2とΔN3とにおいて、異なるスペクトル幅を持つと判断する。このため、m=2のスペクトル幅ΔfD(2)をドップラースペクトル幅算出の候補から除外する。
スペクトル幅ΔfD(2)が除外された場合、m=2の処理を終了し、m=3として、同様の処理を行う。差分ΔfD(3)-ΔfD(4)は周波数差のしきい値fthreshを超えないため、ΔN3とΔN4とにおいて、同じスペクトル幅を持つと判断する。
上記の処理において、m=1と2のスペクトル幅が除外されている。このため、m=3のスペクトル幅ΔfD(3)が、ドップラースペクトル幅として決定される。
図10に上記の処理のフローチャートを示す。これは、図4に示されるフローチャートのS14〜S17に相当する。
図10のS31〜S33は、図9のS21〜S23と同じである。S34において、まず、m=1と設定し、m=1のしきい値におけるスペクトル幅の検査を行う。S35において、m’=m+1として比較対象のしきい値を決定する。S36において、mのしきい値におけるスペクトル幅とm’のしきい値におけるスペクトル幅の差分が周波数差のしきい値fthreshより大きくなるかを判定する。差分が周波数差のしきい値より大きい場合、m=1において、端末単体のドップラースペクトルの外側にスペクトル成分が存在するとして、m=1のスペクトル幅を候補から除外する。S37において、m=M-1でないとして、m=1の処理を終了する。続いて、S38において、mを1増やして同様の処理を続ける。
S36において、スペクトル幅の差分が周波数差のしきい値より小さい場合、S39において、m’がMに到達しているかを判断する。m’=Mでない場合、S40において、m’を1増やす。すなわち、比較するしきい値を1つ上として、S36において、スペクトル幅の差分が周波数差のしきい値より大きくなるかを判定する。
m=1において、S36、S39、S40の処理により、m’=2〜Mのしきい値のスペクトル幅との比較を行い、全てのしきい値において差分が周波数差のしきい値を超えない場合、S39において、m’=Mとなり、m=1のスペクトル幅ΔfD(1)がドップラースペクトル幅に相当すると判断される。
m=1のスペクトル幅がドップラースペクトル幅に相当すると判断されると、m=2〜Mのスペクトル幅も、m=1のスペクトル幅に対して周波数差のしきい値内にあり、ドップラースペクトル幅に相当すると判断される。スペクトル幅ΔfD(1)〜ΔfD(M)の中で、ΔfD(1)が最も広くなる可能性があるため、S41において、ΔfD(1)をドップラースペクトル幅と決定し、速度の推定に使用される。
S36において、mのスペクトル幅が、m’のスペクトル幅より周波数差のしきい値fthreshを超えて大きいと判定されると、mの処理を終了する。S38において、mを1増やして同様の処理を続ける。
S37において、m=M-1となった場合、m=1〜M-1のしきい値において、周辺移動散乱体によるドップラースペクトル成分があると判断し、S42において、ΔfD(M)をドップラースペクトル幅として決定する。
上記実施例において、Ppeakにおけるスペクトル幅の算出は考慮されない。図8に示すドップラースペクトルのようにドップラースペクトルの両端にピーク電力がある場合、ピーク電力位置で、スペクトル幅を算出することができる。しかし、一般的には、不均等な到来電波方向や雑音などの影響により、ドップラースペクトルの一方だけ電力が高くなる、あるいは、スペクトル両端の間に雑音によるピーク電力が生じることがある。このような場合、Ppeak位置において、スペクトル幅を算出することができない。このため、最も高いしきい値(m=M)におけるスペクトル幅までをドップラースペクトル幅算出の対象としている。
しかし、端末近傍のSNR環境が良く、受信電波がほぼ全方向から到来するような場合、図8に示すように、スペクトル両端に同電力のピーク電力が生じる。このような場合において、Ppeak位置におけるスペクトル幅を比較対象とすることもできる。
本発明のその他の実施例を以下に説明する。図11に示すようなドップラースペクトルを想定する。即ち、ドップラースペクトル内の特定の周波数成分の電力が極端に高くなる。この極端に高い電力に対してピーク電力Ppeakが設定される。設定されたPpeakに基づいて、しきい値ΔN1〜ΔN4を設定すると、最も高いしきい値m=4のスペクトル幅ΔfD(4)は、実際のドップラースペクトル幅より狭く算出される可能性がある。このとき、m=3のスペクトル幅ΔfD(3)は、ΔfD(4)より広くなるため、周辺移動散乱体によるドップラースペクトル成分を含むと判断される。この結果、ΔfD(4)がドップラースペクトル幅として決定され、端末の移動速度は、実際より遅く推定されてしまう。
ドップラースペクトルのピーク電力を実際のピーク電力より低く設定することにより、極端に高いピーク電力の影響を抑えることができる。図12に示すように、実際のピーク電力PpeakよりΔPpeak低い電力(Ppeak-ΔPpeak)を新たなピーク電力に設定し、(Ppeak-ΔPpeak)と雑音電力N0の間にしきい値ΔN1〜ΔN4を設定する。最も高いしきい値m=4のスペクトル幅ΔfD(4)は、ドップラースペクトル幅に相当することになり、正確なドップラースペクトル幅を算出することを可能にする。ΔPpeakは、伝搬環境に応じて、適切に設定することができる。
図13は、本発明の装置構成を簡単に説明するブロック図である。このような構成は、基地局側、端末側のどちらにも設置することができ、両局のどちらにおいても、最大ドップラー周波数の測定を可能にする。
図13の(a)において、50は信号受信部、51は受信信号電力測定部、52はFFT、53はスペクトル幅決定部、54は最大ドップラー周波数算出部、55は端末速度推定部を表す。また、(b)はスペクトル幅決定部53を表し、56はしきい値設定部、57はスペクトル幅算出部、58はスペクトル幅比較部を表す。
信号受信部50は、無線通信を行う送信側局からの信号を受信し、受信信号の中から参照信号を取り出す。参照信号の挿入方法は、システムにより規定される。受信信号電力測定部51は、受信した参照信号をベースバンド信号に変換し、復調した後、受信電力を測定して、受信信号の伝搬路の時間変動データを生成する。受信信号電力測定部51は、伝搬路推定の一部として実装することも可能である。FFT 52は、伝搬路の時間変動データにFFTを掛けて、ドップラースペクトルを算出する。スペクトル幅決定部53は、ドップラースペクトルに対して設定される複数のしきい値において算出されるスペクトル幅からドップラースペクトル幅を決定する。最大ドップラー周波数算出部54は、スペクトル幅決定部53で決定されたドップラースペクトル幅から、最大ドップラー周波数を算出する。端末速度推定部55は、検出された最大ドップラー周波数より、端末の移動速度を算出する。
スペクトル幅決定部53において、しきい値設定部56は、ドップラースペクトルに対して複数のしきい値を設定する。スペクトル幅算出部57は、各しきい値におけるスペクトル幅を算出する。スペクトル幅比較部58は、各しきい値のスペクトル幅を比較し、周辺移動散乱体によるドップラースペクトル成分を含むスペクトル幅を検知して除外すると共に、適切なスペクトル幅を選択する。
50 信号受信部
51 受信信号電力測定部
52 FFT
53 スペクトル幅決定部
54 最大ドップラー周波数算出部
55 端末速度推定部
56 しきい値設定部
57 スペクトル幅算出部
58 スペクトル幅比較部

Claims (1)

  1. 端末の移動速度を推定する方法であって、
    参照信号を受信するステップと、
    前記受信した参照信号から、伝搬路の時間変動データを作成するステップと、
    前記伝搬路の時間変動データにFFTを掛けるステップと、
    前記FFTの出力からドップラースペクトルを算出するステップと、
    前記ドップラースペクトルの成分のうち、最も電力が高い成分の電力値をピーク電力と設定するステップと、
    雑音電力を算出するステップと、
    前記ピーク電力と前記雑音電力の間に、等間隔で所定個数のしきい値を設定するステップと、
    前記所定個数のしきい値の各々において、しきい値以上となる前記ドップラースペクトルの成分のうち、最も高い周波数の成分と最も低い周波数の成分の周波数の差分から、スペクトル幅を算出するステップと、
    前記所定個数のしきい値において前記算出されたスペクトル幅のうち、最も高いしきい値のスペクトル幅を第1のドップラースペクトル幅候補とするステップと、
    前記最も高いしきい値より低いしきい値のスペクトル幅の各々と、前記第1のドップラースペクトル幅候補との差分を算出し、前記算出された差分が周波数差のしきい値以下となるスペクトル幅を第2のドップラースペクトル幅候補とするステップと、
    前記第2のドップラースペクトル幅候補がない場合、前記第1のドップラースペクトル幅候補をドップラースペクトル幅と決定するステップと、
    前記第2のドップラースペクトル幅候補がある場合、最も低いしきい値の前記第2のドップラースペクトル幅候補をドップラースペクトル幅と決定するステップと、
    前記決定されたドップラースペクトル幅から最大ドップラー周波数を算出するステップと、
    前記最大ドップラー周波数から端末の移動速度を決定するステップと、
    を有する、移動速度の推定方法。
JP2016060664A 2016-03-24 2016-03-24 ドップラースペクトルを用いた端末の移動速度推定方法 Active JP6239670B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2016060664A JP6239670B2 (ja) 2016-03-24 2016-03-24 ドップラースペクトルを用いた端末の移動速度推定方法
PCT/JP2017/010027 WO2017163988A1 (ja) 2016-03-24 2017-03-13 ドップラースペクトルを用いた端末の移動速度推定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016060664A JP6239670B2 (ja) 2016-03-24 2016-03-24 ドップラースペクトルを用いた端末の移動速度推定方法

Publications (2)

Publication Number Publication Date
JP2017173171A JP2017173171A (ja) 2017-09-28
JP6239670B2 true JP6239670B2 (ja) 2017-11-29

Family

ID=59899400

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016060664A Active JP6239670B2 (ja) 2016-03-24 2016-03-24 ドップラースペクトルを用いた端末の移動速度推定方法

Country Status (2)

Country Link
JP (1) JP6239670B2 (ja)
WO (1) WO2017163988A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018096882A (ja) * 2016-12-14 2018-06-21 ソフトバンク株式会社 ドップラースペクトルを用いた端末の移動速度推定方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020026266A1 (en) * 2018-07-28 2020-02-06 Wisig Networks Private Limited Method and system for classifying speed of a user equipment
CN115019521B (zh) * 2022-05-19 2023-06-09 河北工业大学 一种确定车辆速度的方法和系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20050113468A (ko) * 2004-05-29 2005-12-02 삼성전자주식회사 속도 추정 장치 및 방법
JP4852723B2 (ja) * 2005-09-27 2012-01-11 エスティー‐エリクソン、ソシエテ、アノニム 最大ドップラー周波数及び発振器周波数オフセットを推定するための方法、プログラム、及びモジュール、そのモジュールを含む受信機
JP5196971B2 (ja) * 2007-11-27 2013-05-15 三菱電機株式会社 目標追尾装置
JP5344984B2 (ja) * 2009-04-28 2013-11-20 三菱電機株式会社 移動速度検出装置及び移動速度検出方法
JP4975085B2 (ja) * 2009-11-09 2012-07-11 三菱電機株式会社 受信装置および方法
JP5976137B2 (ja) * 2015-01-14 2016-08-23 ソフトバンク株式会社 移動通信システムにおけるドップラースペクトルを用いた移動局速度推定方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018096882A (ja) * 2016-12-14 2018-06-21 ソフトバンク株式会社 ドップラースペクトルを用いた端末の移動速度推定方法

Also Published As

Publication number Publication date
JP2017173171A (ja) 2017-09-28
WO2017163988A1 (ja) 2017-09-28

Similar Documents

Publication Publication Date Title
CN105101046B (zh) 无线通信系统中的电子设备和无线通信方法
US7953413B2 (en) Handover for cellular radio systems
US9526059B2 (en) Mobility management based on cell size
JP6239670B2 (ja) ドップラースペクトルを用いた端末の移動速度推定方法
EP3017623B1 (en) Method and apparatus of switching communications from a first channel to a second channel of higher-frequency
JP5976137B2 (ja) 移動通信システムにおけるドップラースペクトルを用いた移動局速度推定方法
KR101682030B1 (ko) 무선통신시스템에서 이동 기지국의 송신전력을 제어하기 위한 방법 및 장치
JP6239671B2 (ja) ドップラースペクトルを用いた端末速度推定方法
Wu et al. Reducing handoff delay of wireless access in vehicular environments by artificial neural network-based geographical fingerprint
JP6434479B2 (ja) ドップラースペクトルを用いた端末の移動速度推定方法
JP6180496B2 (ja) 端末速度推定方法
KR101616759B1 (ko) 이동통신 시스템에서 주파수 할당 선택 방법 및 장치
EP3050356B1 (en) Radio resource management measurements in dense network scenarios
CN105981435B (zh) 用于相邻小区相关性配置的方法和设备
KR101285367B1 (ko) 무선통신 시스템에서 이종 망 간의 핸드오프의 수행을결정하기 위한 장치 및 방법
Pourmina et al. Load balancing algorithm by vertical handover for integrated heterogeneous wireless networks
EP3016458A1 (en) Apparatus, Mobile Device, Base Station Transceiver, Adaptation Server, Method and Computer Program for providing information related to a predicted channel state
Maiti et al. Field experimental evaluation of mobile terminal velocity estimation based on Doppler spread detection for mobility control in heterogeneous cellular networks
KR101436187B1 (ko) 단말기에서 수신신호 세기와 시신속도를 이용한 핸드오버 장치 및 방법
Tseng et al. A predictive movement based handover algorithm for broadband wireless networks
Koudouridis et al. A method for the generation of radio signal coverage maps for dense networks
McGuire et al. A robust fuzzy logic handoff algorithm
Miyim et al. Regressive Prediction Approach to Vertical Handover in Fourth Generation Wireless Networks
Raiyn A Cognitive Radio Scheme for Dynamic Resource Allocation Based on QoE
Aldebes et al. Robust Vertical Handover Decision based on Two Thresholds and Multiple Parameters

Legal Events

Date Code Title Description
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: 20171010

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171101

R150 Certificate of patent or registration of utility model

Ref document number: 6239670

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250