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

Abstract

【課題】周辺移動散乱体によるドップラースペクトル成分の発生を検知する移動速度推定方法を提供する。
【解決手段】移動速度の推定方法は、ドップラースペクトルを算出するステップと、最大ドップラー周波数を算出するステップと、最大ドップラー周波数から端末の第１の移動速度を決定するステップと、第１の移動速度と所定時間前に決定された端末の第２の移動速度とを比較し、第１の移動速度が、第２の移動速度に第１の係数を掛けた値以下の場合、第１の移動速度を現在時間の移動速度と決定するステップと、を有する。
【選択図】図８

Description

本発明は、無線通信システムにおけるドップラースペクトルを用いた端末の速度推定方法に関し、特に、現在時刻と前後時間の速度推定値との比較結果から信頼性の確認を行う端末の速度推定方法に関する。

無線通信システムは、ユーザによる無線通信の多様な利用を想定し、低データレート通信から高データレート通信、静止画から動画ストリーミング伝送、静止状態や歩行状態での利用から高速移動状態での利用、あるいは屋内での使用から屋外での使用など様々な形態による無線通信サービスを提供している。

無線通信システムを提供する事業者は、このような利用形態に加え、都市部や農村部などサービスを提供する場所、サービスを利用するユーザの人口分布やデータ量（トラフィック密度）など、様々な条件を考慮したシステム構築を行っている。システム構築の一環として、事業者は、どのような場所や時間においても、ユーザが遅滞なく無線通信システムに接続することができ、満足するサービスの提供を受けられるように、利用形態やトラフィック密度を考慮して基地局（またはセル）の配置設計を行う。

しかし、利用形態の更なる多様化やトラフィック量の増加に伴い、平面的なセル配置では、十分なサービスを提供することは困難となる。このため、複数の種類のセルを組み合わせることにより、多様な利用形態やトラフィック特性に合致したサービスを提供するオーバーレイセル構成の構築が進められている。オーバーレイセル構成の１つの例において、セルサイズの大きいマクロセル内にセルサイズの小さい極小セルを重畳（オーバーレイ）させる。極小セルは、マイクロセル、ピコセル、フェムトセルなどである。このような極小セルは、マクロセル内のトラフィック密度が高い場所、高データレート通信が想定される場所、マクロセル基地局からの電波受信が困難な場所等に設置される。

オーバーレイセル構成において、周波数利用効率を向上させるため、あるいは、ハンドオーバの頻度を抑制するために、端末の移動速度に応じて、マクロセルと極小セルの階層を適切に選択する階層選択制御が必要となる。

移動速度が速い端末は、マクロセルのみに接続し、極小セルには接続されないようにする。一方、移動速度の遅い端末は、極小セル内にある場合は極小セルと接続し、極小セルが重畳されていない場所では、マクロセルと接続する。このような階層選択制御により、高速移動中の端末が極小セルに頻繁にハンドオーバすることを防止して、ハンドオーバの制御シグナリングの増加やハンドオーバ増加に伴う通信断を減らすことができる。また、低速移動の端末が極小セルに接続することで、送信電力が抑制されて、システム内の干渉電力を減らすことができる。

階層選択制御を適切に実行するためには、端末の移動速度を正確に測定し、その移動速度に基づいて、マクロセル又は極小セルに接続させることが重要である。

端末の移動速度を推定する方法の１つとして、ドップラースペクトルを利用する方法がある（非特許文献１）。ここでドップラースペクトルとは、伝搬路応答の時間変動をフーリエ変換した結果得られる周波数軸（ドップラー周波数シフト）上の電力波形である。端末と基地局間の無線伝搬環境が電波の到来方向分布が一様である全周散乱モデルと見なせる場合、ドップラースペクトルは、最大ドップラー周波数においてピークを示す。このため、スペクトルのピークを検出することにより、最大ドップラー周波数を測定することができ、最大ドップラー周波数から端末の移動速度を算出することができる。

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月．

端末が自動車内にあり、その端末搭載自動車が一般道路や高速道路を走行すると、対向車や同方向車などの自動車が端末搭載自動車の近傍を走行する。このような対向車や同方向車は、端末から送信される電波、あるいは端末で受信される電波を散乱させる。このような自動車を周辺移動散乱体とする。周辺移動散乱体の存在により、ドップラースペクトルの算出において、端末単体によるドップラースペクトルとは異なるドップラースペクトル成分が発生する。この異なるドップラースペクトル成分は、端末単体によるドップラースペクトルより低い周波数成分または高いスペクトル周波数成分を有している。また、通常は、端末単体によるドップラースペクトルの電力より小さい電力となる。対向車による速度推定値への影響が確認されている文献として、非特許文献２、非特許文献３を挙げる。

端末のドップラースペクトル検出において、周辺移動散乱体によるドップラースペクトル成分を有効成分として検出してしまうと、端末単体によるドップラースペクトルより広いスペクトルとなり、最大ドップラー周波数が高くなってしまう。このため、周辺移動散乱体が存在すると、端末の移動速度が実際より早い速度として推定される可能性がある。

端末のドップラースペクトルを算出し、算出されたドップラースペクトルから端末の移動速度を推定する方法において、現在時刻で算出された移動速度と現在時刻の前後の時刻において算出された移動速度を比較することにより、現在時刻の移動速度に周辺移動散乱体による影響があるか否かを判定し、より正確な端末の移動速度を推定することができる。

本発明は、ドップラースペクトルを用いた端末の移動速度推定において、周辺移動散乱体によるドップラースペクトル成分が存在する場合においても、周辺移動散乱体によるドップラースペクトル成分の発生と、端末の加速又は減速に伴うドップラースペクトルの変動を区別し、正確な移動速度推定を行う移動速度推定方法を提供することを目的とする。

本発明は、端末の移動速度を推定する方法であって、参照信号を受信するステップと、前記受信した参照信号から、伝搬路の時間変動データを作成するステップと、前記伝搬路の時間変動データにＦＦＴを掛けるステップと、前記ＦＦＴの出力からドップラースペクトルを算出するステップと、前記ドップラースペクトルに対するしきい値を設定し、前記ドップラースペクトルの成分のうち、前記しきい値以上となる有効成分を選択するステップと、前記有効成分のうち、最も高い周波数の成分の第１の周波数を検出するステップと、前記有効成分のうち、最も低い周波数の成分の第２の周波数を検出するステップと、前記第１の周波数と前記第２の周波数の差分からドップラースペクトル幅を算出するステップと、前記ドップラースペクトル幅から最大ドップラー周波数を検出するステップと、前記最大ドップラー周波数から端末の第１の移動速度を決定するステップと、前記第１の移動速度と所定時間前に決定された前記端末の第２の移動速度とを比較し、前記第１の移動速度が、前記第２の移動速度に第１の係数を掛けた値以下の場合、前記第１の移動速度を現在時間の移動速度と決定するステップと、を有する。

本発明に係る移動速度推定方法において、周辺移動散乱体によるドップラースペクトル成分を検出し、端末の加速又は減速に伴うドップラースペクトルの変動と区別することにより、端末の移動速度を正確に推定することができる。

ドップラースペクトルの例を示す図である。 最大ドップラー周波数の検出方法を説明する図である。 ドップラースペクトルへのしきい値の適用を説明する図である。 最大ドップラー周波数の推定方法のフローチャートを示す図である。 周辺移動散乱体によるドップラースペクトル成分を示す図である。 端末搭載車と対向車との関係を示す図である。 周辺移動散乱体によるドップラースペクトル成分の変動を示す図である。 本発明による周辺移動散乱体によるドップラースペクトル成分の検出方法を説明する図である。 本発明による周辺移動散乱体によるドップラースペクトル成分の検出方法のフローチャートを示す図である。 本発明による装置構成の例を示す図である。

以下図面を参照して、ドップラースペクトルを用いた端末の速度推定方法について説明する。しかしながら、本発明が、図面又は以下に記載される実施形態に限定されるものではないことを理解されたい。

図１にドップラースペクトルの例を示す。この例は、無線伝搬環境が全周散乱モデルとみなせる理想的な場合を表している。図において、横軸は周波数を表し、縦軸は相対電力を表す。ただし、図１ではベースバンド帯域に変換した後の処理結果を示しているため、実際には横軸の周波数0 Hzがキャリア周波数に対応することに注意が必要である。以降のドップラースペクトルの図においても横軸の周波数の定義は図１と同様である。この例において、最大ドップラー周波数は、100 Hzであり、ドップラースペクトルの周波数100 Hz及び-100 Hzにおいて、高い電力の成分が生じている。このように、受信信号の伝搬路変動を用いて算出されたドップラースペクトルから、電力の高い成分の周波数を検出することで、最大ドップラー周波数を推定することができる。

図２は、ドップラースペクトルから最大ドップラー周波数を算出する方法を簡単に説明した図である。周波数の負側において、最大電力が得られる周波数をf(1)とし、正側において、最大電力が得られる周波数をf(2)とする。理想的なフェージング環境下においては、f(1)とf(2)の間でU字型のスペクトルを示し、f(1)とf(2)は、同じ周波数絶対値を示す。この周波数絶対値が最大ドップラー周波数f_Dとなる。即ち、f(1)=-f_Dであり、f(2)=f_Dとなる。そして、f(1)とf(2)の間の周波数がドップラースペクトル幅となる。

最大ドップラー周波数f_Dが得られると、端末の移動速度νは、f_D×λで算出することができる。ここで、λは、キャリア周波数の波長である。

ドップラースペクトルから最大ドップラー周波数を検出する方法が幾つか考えられる。以下に、ドップラースペクトルに対して固定のしきい値を設定し、しきい値以上となる有効成分について、そのうちの最大周波数の成分と最小周波数の成分から最大ドップラー周波数を算出する方法を説明する。

図３は、ドップラースペクトルにおけるしきい値の設定の一例を示した図である。ここでは、雑音電力N₀ [dBm]を算出し、更に、ΔN [dB]を加算した値をしきい値とする。図３に示されるように、N₀+ΔN [dBm]をしきい値とすることで、ドップラースペクトルの外側の雑音成分を除外することができる。そして、しきい値以上となる有効成分は、直接波成分など幾つかの電力が高い成分を含む。これらの有効成分の中から、最大の周波数を持つ成分（即ち、図３の「Ａ」で示される100 Hzの成分）と最小周波数を持つ成分（即ち、図３の「Ｂ」で示される-100 Hzの成分）を選択することで、最大ドップラー周波数を決定することができる。

図４は、最大ドップラー周波数算出の方法を示すフローチャートである。S10において、無線接続している基地局または端末から参照信号を受信する。ここで、参照信号としたが、システムにより、パイロット信号やビーコンなどとも呼ばれ、伝搬路推定に使用される既知の信号である。S11において、受信した参照信号をベースバンド信号に変換し、復調した後、受信電力を測定することにより、受信信号の伝搬路の時間変動データを作成する。伝搬路の時間変動データは、一定区間毎にバッファに格納される。

受信信号は、端末の移動速度に比例して変動するため、受信信号の伝搬路変動をフーリエ変換することにより、ドップラースペクトルが得られる。S12において、一定区間毎に格納された伝搬路の時間変動データに高速フーリエ変換（FFT）を掛ける。ここで、一定区間は、FFTのポイント数に相当する。S13において、FFTの出力からドップラースペクトルが算出される。

S14において、スペクトルのしきい値を設定し、しきい値以上となる有効成分を検出する。S15において、有効成分の中から最大周波数を持つ有効成分E_{f_max}を選択し、その周波数f(E_{f_max})を検出する。S16において、有効成分の中から最小周波数を持つ有効成分E_{f_min}を選択し、その周波数f(E_{f_min})を検出する。

S17において、検出された２つの周波数を用いて、最大ドップラー周波数を算出する。ここでは、２つの周波数f(E_{f_max})と周波数f(E_{f_min})からドップラースペクトル幅を算出し、その半値を最大ドップラー周波数f_Dとすることができる。即ち、f_D=(f(E_{f_max})−f(E_{f_min}))/2とする。

このように、しきい値を設定し、しきい値以上となる有効成分を抽出することで、雑音成分によるスペクトルを除外して、受信信号によるスペクトルのみを選択することができる。さらに、適切なしきい値を設けることで、雑音成分による誤検出を抑制し、より正確な最大ドップラー周波数を推定することができる。

上記例において、しきい値を雑音電力N₀にΔNを加算した値に設定したが、これに限定されるものではない。例えば、測定されるSNRなどの無線伝搬状況に応じて、しきい値を可変とすることも可能である。

また、周波数f(E_{f_max})と周波数f(E_{f_min})とのうち、絶対値の大きい方を最大ドップラー周波数f_Dとする方法もある。即ち、f_D=max(|f(E_{f_max})|, |f(E_{f_min})|)とする。

次に、周辺移動散乱体によるドップラースペクトル算出への影響とその検知方法について説明する。端末搭載自動車が、一般道路や高速道路を走行する際、対向車や同方向車が存在する。このように端末搭載自動車周辺を移動する車両体（周辺移動散乱体）は、端末から出力される電波、または端末で受信される電波を散乱させる。

周辺移動散乱体が存在すると、端末のドップラースペクトル検出において、端末単体によるドップラースペクトルとは別に、周辺移動散乱体のよるドップラースペクトル成分が発生する。簡単な例を図５に示す。端末単体によるドップラースペクトルは、端末の移動速度に応じた最大ドップラー周波数f₀と-f₀により示される。この端末単体によるドップラースペクトルより高い周波数位置（f₁）において、周辺移動散乱体によるドップラースペクトル成分が発生している。このドップラースペクトル成分の電力は、通常、端末単体によるドップラースペクトルの電力より小さくなる。

しかし、周辺移動散乱体によるドップラースペクトル成分を端末のドップラースペクトルの有効成分として検出してしまうと、ドップラースペクトル幅は、-f₀〜f₁となり、本来のドップラースペクトル幅より広く検出される。図４に示すフローチャートに従い、最大ドップラー周波数を算出すると、ドップラースペクトル幅の半値は(f₁＋f₀)/2となり、実際の最大ドップラー周波数より大きくなる。また、最大周波数と最小周波数の絶対値の大きい方から最大ドップラー周波数を算出しても最大ドップラー周波数はf₁となる。

図６に端末が搭載される自動車Ａと対向車Ｂの関係を簡単に示す。対向車Ｂが周辺移動散乱体となる。自動車Ａは速度v₀で走行し、その最大ドップラー周波数はf₀とする。一方、対向車Ｂは速度vで走行し、その最大ドップラー周波数はfとする。また、時刻tにおける自動車Ａと対向車Ｂとがなす角度をθ(t)とすると、自動車Ａと対向車Ｂの相対速度v₁は、v₁＝v₀＋v×cosθ(t)となる。このため、対向車Ｂによるドップラースペクトル成分の周波数f₁は、f₁＝f₀＋f×cosθ(t)となる。

対向車Ｂが自動車Ａに近づくと、角度θ(t)が90°に近づくため、対向車Ｂによるドップラースペクトル成分の周波数f₁は、f₀に近づく。即ち、図７に示されるように、対向車Ｂによるドップラースペクトル成分は、自動車Ａ単体によるドップラースペクトルより高い周波数位置に現れ、両者が接近するに従い、自動車Ａ単体によるドップラースペクトルの中心に近づく。

対向車Ｂが自動車Ａの横を通り、離れると、対向車Ｂによるドップラースペクトル成分の周波数f₁は、f₁= -f₀-f×cosθ(t)となり、自動車Ａ単体によるドップラースペクトルより低い周波数位置に移動する。

周辺移動散乱体が同方向車の場合も同様に考えることができる。同方向車の速度が自動車Ａの速度より速く、自動車Ａを追い抜いて行く場合、同方向車によるドップラースペクトル成分は、自動車Ａ単体によるドップラースペクトルより高い周波数位置で発生し、追い抜いた後、低い周波数位置に移動する。

更に、自動車Ａの近傍に対向車や同方向車が複数存在する場合、自動車Ａ単体によるドップラースペクトルの片外側または両外側に複数のドップラースペクトル成分が出現することも考えられる。以下に示す本発明による速度推定方法は、このように複数のドップラースペクトル成分が存在している場合においても正確に速度推定を行うことができる。

上記のように周辺移動散乱体によるドップラースペクトル成分が存在すると、端末のドップラースペクトルの正確な検出を妨げることになる。周辺移動散乱体によるドップラースペクトル成分を検出し、端末のドップラースペクトルと区別することで、端末単体によるドップラースペクトルを算出し、正確な速度推定を行うことができる。

本発明の１つの実施例における周辺移動散乱体によるドップラースペクトル成分の検出方法を説明する。

図８は、端末の移動速度の推定を行う時間と、その時のドップラースペクトルの例を示す図である。移動速度の推定は、時間... t_-3、t_-2、t_-1、t₀、t₁、t₂、t₃ ...において行われるとする。その時間間隔Δt_iはシステムや端末により設定される。

現在時刻t₀において算出されるドップラースペクトルP(t₀)は、端末単体のみによるドップラースペクトルであるか、端末単体のドップラースペクトルと周辺移動散乱体によるドップラースペクトル成分とが混在したスペクトルであるかの区別はできない。このため、ドップラースペクトル幅から算出される現在時刻t₀における速度推定値v(t₀)が、実際の速度を表すものであるかの信頼性を確認する。

現在時刻t₀における速度推定値v(t₀)の信頼性を確認するために、ΔT前の速度推定値v(t_-ΔT)との比較を行う。即ち、時間t_-ΔTとt₀との間で、ドップラースペクトル幅に大きな変動があり、速度推定値に影響が生じたかを判定する。ΔTは任意に設定することができる。速度推定値の比較においては、速度推定時の誤差や実際の移動速度の若干の変化を考慮し、現在時刻t₀における速度推定値v(t₀)が、ΔT前の速度推定値v(t_-ΔT)に係数αで重み付けした速度より大きくなるかを調査する。即ち、v(t₀)>α×v(t_-ΔT)となるかを判定する。αは、1.3程度を想定するが、これに制限されるものではない。

v(t₀)がα×v(t_-ΔT)以下となる場合、時刻t_-ΔT〜t₀において、次の例が想定される。
１．周辺移動散乱体がなく、車の速度が維持されている、又は若干速くなっているが、重み付け係数αの範囲内である。
２．周辺移動散乱体がなく、車が減速中である。
３．時刻t₀において周辺移動散乱体によるドップラースペクトル成分が存在するが、重み付け係数αの範囲内である。
４．時刻t_-ΔTにおいて周辺移動散乱体によるドップラースペクトル成分が存在していたが、時刻t₀において周辺移動散乱体が存在しない。

v(t₀)がα×v(t_-ΔT)以下となる場合、速度推定値v(t₀)の信頼性があると判断され、時刻t₀の速度推定値としてv(t₀)がそのまま使用される。

一方、v(t₀)がα×v(t_-ΔT)を超える場合、時刻t_-ΔT〜t₀において、次の例が想定される。
１．周辺移動散乱体がなく、車が加速中である。
２．時刻t_-ΔTにおいて周辺移動散乱体がなく、時刻t₀において周辺移動散乱体によるドップラースペクトル成分が発生している。

図８の例では、現在時刻t₀におけるドップラースペクトルP(t₀)には、端末単体のドップラースペクトルに加え、周波数の低い側に周辺移動散乱体によるドップラースペクトル成分が存在している。このため、算出されるドップラースペクトル幅は、実際の端末単体のドップラースペクトル幅に比べ広くなり、移動速度の推定値v(t₀)も実際より高くなる。

図８では、時刻t_-3を現在時刻のΔT前としている。時刻t_-3において算出されたドップラースペクトルP(t_-3)は、端末単体のドップラースペクトルのみを含んでいる。このため、現在時刻t₀において検出されたドップラースペクトル幅は、時刻t_-3において検出されたドップラースペクトル幅に比べ広くなり、移動速度の推定値v(t₀)は、v(t_-3)より高くなる。

上記の処理において、v(t₀)がα×v(t_-ΔT)を超えると判断されると、更なる信頼性を確認するため、速度推定値v(t₀)とΔT後の速度推定値v(t_ΔT)との比較を行う。速度推定値の比較において、速度推定値v(t_ΔT)と速度推定値v(t₀)を係数αで重み付けした速度とを比較する。即ち、v(t_ΔT)>α×v(t₀)であるかを判定する。αは、1.3程度を想定するが、これに制限されるものではない。また、ΔT前の速度推定値との比較において用いたαと同じ値とすることも、異なる値とすることもできる。

v(t_ΔT)がα×v(t₀)を超える場合、時刻t₀〜t_ΔTにおいて、t_-ΔT〜t₀に継続して、車が加速中であることが想定される。この場合、時刻t₀における速度推定値v(t₀)の信頼性があると判断され、時刻t₀における速度推定値としてv(t₀)がそのまま使用される。

一方、時刻t₀において周辺移動散乱体によるドップラースペクトル成分が発生したことにより、速度推定値v(t₀)が増加した後、時刻t_ΔTにおいてそのドップラースペクトル成分が消滅した場合、速度推定値は減少する。このため、v(t_ΔT)はα×v(t₀)より小さくなる。この場合、時刻t₀における速度推定において、周辺移動散乱体によるドップラースペクトル成分の影響があったものと判断し、速度推定値v(t₀)をΔT前の速度推定値v(t_-ΔT)に設定する。

図８の例では、現在時刻t₀におけるドップラースペクトルP(t₀)の算出において、端末単体のドップラースペクトルに加えて、周辺移動散乱体によるドップラースペクトル成分が存在している。このため、移動速度の推定値v(t₀)も高くなる。一方、ΔT後である時刻t₃において、周辺移動散乱体によるドップラースペクトル成分は消滅し、端末単体によるドップラースペクトルのみ残っている。このため、速度推定値v(t₃)は、α×v(t₀)より小さくなる。この場合、時刻t₀において、周辺移動散乱体による影響があったと判断し、現在時刻の速度推定値v(t₀)は、v(t_-3)に設定される。

図９に上記処理のフローチャートを示す。これは、図４に示されるフローチャートのS17の後に行われるものである。

図９のS21において、現在時刻t₀におけるドップラースペクトル幅の算出から、移動速度v(t₀)の推定を行う。S22において、移動速度推定値v(t₀)の信頼性を検証するため、ΔT前の移動速度推定値v(t_-ΔT)に係数αを掛けた速度との比較を行う。v(t₀)がα×v(t_-ΔT)以下の場合、時刻t₀において、周辺移動散乱体の影響がなく、時刻t_-ΔTから速度が維持されているか減速中であると判断し、S25において、v(t₀)= v(t₀)と設定する。

S22において、v(t₀)がα×v(t_-ΔT)を超える場合、S23において、ΔT後の移動速度推定値v(t_ΔT)とv(t₀)に係数αを掛けた速度との比較を行う。v(t_ΔT)がα×v(t₀)を超える場合、時刻t_-ΔTからt_ΔTにかけて端末が加速中であり、周辺移動散乱体の影響がないと判断し、S24において、v(t₀)= v(t₀)と設定する。これに対し、v(t_ΔT)がα×v(t₀)以下の場合、時刻t₀において、周辺移動散乱体の影響があったと判断し、S26において、v(t₀)= v(t_-ΔT)と設定する。

図１０は、本発明の装置構成を簡単に説明するブロック図である。このような構成は、基地局側、端末側のどちらにも設置することができ、両局のどちらにおいても、最大ドップラー周波数の測定を可能にする。

図１０において、30は信号受信部、31は受信信号電力測定部、32はFFT、33は信号選択部、34は最大ドップラー周波数算出部、35は端末速度推定部、36はバッファ、37は端末速度比較部を表す。信号受信部30は、無線通信を行う送信側局からの信号を受信し、受信信号の中から参照信号を取り出す。参照信号の挿入方法は、システムにより規定される。受信信号電力測定部31は、受信した参照信号をベースバンド信号に変換し、復調した後、受信電力を測定して、受信信号の伝搬路の時間変動データを生成する。受信信号電力測定部31は、伝搬路推定の一部として実装することも可能である。FFT 32は、伝搬路の時間変動データにFFTを掛けて、ドップラースペクトルを算出する。信号選択部33は、ドップラースペクトルの中より、最大ドップラー周波数を算出するために使用する成分を選択する。図４のフローチャートにおいて、信号選択部33は、しきい値を設定し、有効成分中の最大周波数の成分と最小周波数の成分を選択する。最大ドップラー周波数算出部34は、信号選択部33で選択された成分の周波数から、最大ドップラー周波数を算出する。図４の例において、２つの周波数の成分からドップラースペクトル幅を算出し、その半値を最大ドップラー周波数とする。端末速度推定部35は、検出された最大ドップラー周波数より、端末の移動速度を算出する。バッファ36は算出された端末の移動速度を所定期間保存する。端末速度比較部37は、現在時刻に算出された移動速度とバッファ36に保存されるΔT時間前の移動速度を比較する。あるいは、ΔT時間後において、バッファ36に保存される現在時刻の移動速度とΔT時間後の移動速度を比較する。

30 信号受信部
31 受信信号電力測定部
32 FFT
33 信号選択部
34 最大ドップラー周波数算出部
35 端末速度推定部
36 バッファ
37 端末速度比較部

Claims (2)

1. 端末の移動速度を推定する方法であって、
   参照信号を受信するステップと、
   前記受信した参照信号から、伝搬路の時間変動データを作成するステップと、
   前記伝搬路の時間変動データにＦＦＴを掛けるステップと、
   前記ＦＦＴの出力からドップラースペクトルを算出するステップと、
   前記ドップラースペクトルに対するしきい値を設定し、前記ドップラースペクトルの成分のうち、前記しきい値以上となる有効成分を選択するステップと、
   前記有効成分のうち、最も高い周波数の成分の第１の周波数を検出するステップと、
   前記有効成分のうち、最も低い周波数の成分の第２の周波数を検出するステップと、
   前記第１の周波数と前記第２の周波数の差分からドップラースペクトル幅を算出するステップと、
   前記ドップラースペクトル幅から最大ドップラー周波数を検出するステップと、
   前記最大ドップラー周波数から端末の第１の移動速度を決定するステップと、
   前記第１の移動速度と所定時間前に決定された前記端末の第２の移動速度とを比較し、前記第１の移動速度が、前記第２の移動速度に第１の係数を掛けた値以下の場合、前記第１の移動速度を現在時間の移動速度と決定するステップと、
   を有する、移動速度の推定方法。
2. 前記第１の移動速度が、前記第２の移動速度に前記第１の係数を掛けた値を超える場合、
   所定時間後に決定される前記端末の第３の移動速度が、前記第１の移動速度に第２の係数を掛けた値を超えるか否かを判定するステップと、
   前記判定するステップにおいて超えると判定された場合、前記第１の移動速度を現在時間の移動速度と決定するステップと、
   前記判定するステップにおいて超えないと判定された場合、前記第２の移動速度を現在時間の移動速度と決定するステップと、
   を更に有する、請求項１に記載の移動速度の推定方法。