JP2011010088A - 最大ドップラー周波数推定装置および最大ドップラー周波数推定方法 - Google Patents

Abstract

【課題】パイロットシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合において、フェージング変動による最大ドップラー周波数の推定精度を向上させること。
【解決手段】変移検出部１０６は、受信信号の時間的な位相変移を検出するＣＰＩＣＨシンボルの受信パワー測定部１０４は、位相変移に対するフェージング変動以外の要因による影響の大きさを推定する。変移回数計数部１０７は、位相変移に対するフェージング変動以外の要因による影響の大きさの推定の結果に基づいて、位相変移の回数をカウントする。ｆｄ推定値検出部１０８は、変移回数計数部１０７によりカウントした位相変移の回数より、最大ドップラー周波数の推定値を検出する。
【選択図】図２

Description

本発明は、最大ドップラー周波数推定装置および最大ドップラー周波数推定方法に関し、特にフェージング変動の速さに関連する最大ドップラー周波数を精度良く検出する最大ドップラー周波数推定装置および最大ドップラー周波数推定方法に関する。

移動無線通信においては、フェージングによって受信信号の包絡線および位相が変動する。移動無線通信システムによっては、送信機および受信機において既知のパイロットシンボルを送信および受信し、受信機はそのパイロットシンボルの変動を測定する。また、受信機は、パイロットシンボルの変動の測定値に基づいて、フェージング変動が受信信号に与える影響を推定する。この方式は、チャネル推定と呼ばれる。チャネル推定では、雑音の影響を抑圧し、フェージング変動の推定精度を高めるために平均化処理を行う場合が多い。この際、最適な平均区間長は、フェージング変動の速さ、すなわち送受信機の移動速度差に依存することが知られている。そのため、フェージング変動の速さに関連する最大ドップラー周波数を推定する技術が開発されており、無線通信装置単独で推定する方法として、特許文献１の方法が知られている。

図１は、従来のパイロットシンボルの受信信号点のＩ−Ｑ平面における配置を示す図である。

図１に示すように、現在（時刻ｉ）のパイロットシンボルの受信信号点の座標＃１を（Ｉ_ｉ、Ｑ_ｉ）とし、その直前（時刻ｉ−１）のパイロットシンボルの受信信号点の座標＃２を（Ｉ_ｉ−１、Ｑ_ｉ−１）とする。直前の受信信号点の座標＃２から現在の受信信号点の座標＃１への変移は、Ｑ軸およびＩ軸を超えて、第１象限から第３象限へと変移する。この場合、変移の回数は２回と検出される。ここで、受信信号は、フェージング変動およびフェージング変動以外の要因、例えば白色雑音の影響を受けているため、検出された変移の回数は、フェージング変動およびフェージング変動以外の要因による両方の影響を受けた変移の回数である。図１に示すように、フェージング変動のみの影響を受けた場合の現在（時刻ｉ）のパイロットシンボルの信号点の座標＃３が（Ｉ’_ｉ、Ｑ’_ｉ）のとき、直前の受信信号点の座標＃２からフェージング変動のみの影響を受けた受信信号点の座標＃３への変移は、Ｑ軸のみを超えて、第１象限から第２象限へと変移する。この場合、変移回数は１回と検出される。

特開２００４−８６７５号公報

しかしながら、従来の装置においては、図１に示すように、フェージング変動の影響およびフェージング変動以外の要因による影響を受けた場合と、フェージング変動のみの影響を受けた場合とで、パイロットシンボルの信号点の変移する象限が異なり、検出される変移の回数が異なる。従って、パイロットシンボルの受信信号点の変移の回数から検出した最大ドップラー周波数は、フェージング変動のみの影響を受けたときに検出される所望の最大ドップラー周波数と異なり、最大ドップラー周波数を精度良く検出することができないという問題がある。

本発明はかかる点に鑑みてなされたものであり、パイロットシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合において、フェージング変動による最大ドップラー周波数の推定精度を向上させることができる最大ドップラー周波数推定装置および最大ドップラー周波数推定方法を提供することを目的とする。

本発明の最大ドップラー周波数推定装置は、受信信号の時間的な位相変移を検出する位相変移検出手段と、前記位相変移に対するフェージング変動以外の要因による影響の大きさを推定する推定手段と、前記推定の結果に基づいて前記位相変移の回数をカウントする計数手段と、カウントした前記回数より最大ドップラー周波数の推定値を検出する推定値検出手段と、を具備する構成を採る。

本発明の最大ドップラー周波数推定方法は、最大ドップラー周波数を推定する最大ドップラー周波数推定装置における最大ドップラー周波数推定方法であって、受信信号間の時間的な位相変移を検出するステップと、前記位相変移に対するフェージング変動以外の要因による影響の大きさを推定するステップと、前記推定の結果に基づいて前記位相変移の回数をカウントするステップと、カウントした前記回数より最大ドップラー周波数の推定値を検出するステップと、を具備するようにした。

本発明によれば、パイロットシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合において、フェージング変動による最大ドップラー周波数の推定精度を向上させることができる。

従来のパイロットシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態１に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態１に係るＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態１に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図 本発明の実施の形態１に係る閾値の設定方法を示す図 本発明の実施の形態１に係る閾値の設定方法の他の例を示す図 本発明の実施の形態１に係る閾値の設定時の試験結果の一例を示す図 本発明の実施の形態２に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態２に係るＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態２に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図 本発明の実施の形態３に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態３に係るＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態３に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図 本発明の実施の形態４に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態４に係るＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態４に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図 本発明の実施の形態５に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態５に係るＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図 本発明の実施の形態５に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図 本発明の実施の形態５に係る閾値を示す図 本発明の実施の形態６に係る最大ドップラー周波数推定装置の構成を示すブロック図 本発明の実施の形態６に係る最大ドップラー周波数推定装置の動作の一部を示すフロー図

以下、本発明の実施の形態について、図面を参照して詳細に説明する。

（実施の形態１）
図２は、本発明の実施の形態１に係る最大ドップラー周波数推定装置１００の構成を示すブロック図である。最大ドップラー周波数推定装置１００は、例えば、無線受信装置に搭載される。

最大ドップラー周波数推定装置１００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、ＣＰＩＣＨシンボルの受信パワー測定部１０４と、計数判定部１０５と、変移検出部１０６と、変移回数計数部１０７と、最大ドップラー周波数推定値検出部１０８（以下「ｆｄ推定値検出部１０８」と記載する）とから主に構成される。以下に、各構成について詳細に説明する。

アンテナ１０１は、最大ドップラー周波数推定装置１００に到来した信号を受信して無線受信部１０２へ出力する。

無線受信部１０２は、アンテナ１０１から入力した受信信号に対して増幅処理および周波数変換処理を行い、受信直交ベースバンド信号の同相成分（以下、「Ｉ成分」と記載する）および直交成分（以下、「Ｑ成分」と記載する）を生成して取り出す。また、無線受信部１０２は、取り出したＩ成分およびＱ成分をＣＰＩＣＨシンボルのデータ抽出部１０３へ出力する。

ＣＰＩＣＨシンボルのデータ抽出部１０３は、無線受信部１０２から入力したＩ成分およびＱ成分に対し、逆拡散を施しＣＰＩＣＨシンボルのＩ成分およびＱ成分を抽出する。また、ＣＰＩＣＨシンボルのデータ抽出部１０３は、抽出したＣＰＩＣＨシンボルのＩ成分およびＱ成分をＣＰＩＣＨシンボルの受信パワー測定部１０４および変移検出部１０６へ出力する。

ＣＰＩＣＨシンボルの受信パワー測定部１０４は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルの受信電力を測定する。そして、ＣＰＩＣＨシンボルの受信パワー測定部１０４は、測定値を計数判定部１０５へ出力する。

計数判定部１０５は、ＣＰＩＣＨシンボルの受信パワー測定部１０４から入力した受信電力の測定値に基づいて、受信信号に含まれるフェージング変動以外の要因による影響の可能性を推定する。具体的には、計数判定部１０５は、ＣＰＩＣＨシンボルの受信電力が閾値以上の場合にはフェージング変動以外の要因による影響の可能性が小さいと判定し、ＣＰＩＣＨシンボルの受信電力が閾値未満の場合にはフェージング変動以外の要因による影響の可能性が大きいと判定する。また、計数判定部１０５は、フェージング変動以外の要因による影響の可能性が小さいと判定した場合には、変移を計数するように変移回数計数部１０７を制御し、フェージング変動以外の要因による影響の可能性が大きいと判定した場合には、変移を計数しないように変移回数計数部１０７を制御する。ここで、フェージング変動以外の要因とは、雑音等である。

変移検出部１０６は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分に基づいて、ＣＰＩＣＨシンボルの受信信号点が位置するＩ−Ｑ平面上の位相領域の変移を検出する。この変移は、受信信号点が、複数の区切られた位相領域において、異なる位相領域に変移したか否かにより検出する。例えば、この位相領域を、Ｉ−Ｑ平面におけるＩ軸とＱ軸とによって区切られた第１象限から第４象限の４つの象限とする。また、変移検出部１０６は、変移の検出結果を変移回数計数部１０７へ出力する。ここで、ＣＰＩＣＨシンボルの受信信号点は、フェージング変動およびフェージング変動以外の要因の影響によって変移する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部１０５の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。

ｆｄ推定値検出部１０８は、観測時間内のＣＰＩＣＨシンボルの全数に対する変移回数計数部１０７から入力した計数結果の計数値の比率と、予め設定された最大ドップラー周波数ｆｄとの関係式から、最大ドップラー周波数の推定値を検出する。そして、ｆｄ推定値検出部１０８は、検出した最大ドップラー周波数推定値を最大ドップラー周波数推定データ（ｆｄ推定データ）として出力する。以上で、最大ドップラー周波数推定装置１００の構成の説明を終える。

次に、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置について、図３を用いて説明する。図３は、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を示す図である。

図３において、現在（時刻ｉ）におけるＣＰＩＣＨシンボルの受信信号点の座標＃３０１を（Ｉ_ｉ、Ｑ_ｉ）とすると、原点＃３０２とＣＰＩＣＨシンボルの受信信号点の座標＃３０１との距離Ｒ（Ｒ＝√I_i^2+Q_i^2）は、現在（時刻ｉ）のＣＰＩＣＨシンボルの受信電力の平方根により求められる。

ここで、原点＃３０２と受信信号点の座標＃３０１との距離、すなわち受信電力の平方根Ｒが予め設定した閾値Ｐ以上の場合は、フェージング変動以外の要因による影響の可能性が小さく、誤検出の可能性が小さいと推定できる。逆に、受信電力の平方根Ｒが閾値Ｐ未満の場合は、フェージング変動以外の要因による影響の可能性が大きく、誤検出の可能性が大きいと推定できる。

従って、計数判定部１０５は、ＣＰＩＣＨシンボルの受信信号点が横縞の網掛けの範囲内のとき（Ｒ＜Ｐのとき）は、そのＣＰＩＣＨシンボルの受信信号点を変移の回数の計数対象外とし、ＣＰＩＣＨシンボルの受信信号点が横縞の網掛けの範囲外のとき（Ｒ≧Ｐのとき）は、そのＣＰＩＣＨシンボルの受信信号点を変移の回数の計数対象とする。図３においては、計数判定部１０５は、ＣＰＩＣＨシンボルの現在の受信信号＃３０１は、変移の回数の計数対象外と判定する。

次に、最大ドップラー周波数推定装置１００の動作について、図４を用いて説明する。図４は、最大ドップラー周波数推定装置１００の動作の一部を示すフロー図である。

最初に、計数判定部１０５は、ＣＰＩＣＨシンボルの受信電力の平方根（原点とＣＰＩＣＨシンボルの受信信号点との距離）を求め、求めた平方根が予め設定された閾値Ｐ以上か否かを判定する（ステップＳＴ４０１）。

次に、ＣＰＩＣＨシンボルの受信電力の平方根が閾値Ｐ以上の場合には（ステップＳＴ４０１：ＹＥＳ）、計数判定部１０５は、誤検出の可能性が小さいと推定し、変移回数計数時間に「１」を加算する（変移回数計数時間＋１）（ステップＳＴ４０２）。

次に、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ４０３）。

Ｉ成分の符号変化がある場合には（ステップＳＴ４０３：ＹＥＳ）、変移回数計数部１０７は、Ｑ軸を超えて変移したことを検出し、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ４０４）。

一方、Ｉ成分の符号変化がない場合には（ステップＳＴ４０３：ＮＯ）、変移回数計数部１０７は、変移回数に「１」を加算せずに（ステップＳＴ４０４の処理を行わずに）、ステップＳＴ４０５に処理を進める。

また、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ４０５）。

Ｑ成分の符号変化がある場合には（ステップＳＴ４０５：ＹＥＳ）、変移検出部１０６は、Ｉ軸を超えて変移したことを検出し、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ４０６）。

一方、Ｑ成分の符号変化がない場合には（ステップＳＴ４０５：ＮＯ）、変移回数計数部１０７は、変移回数に「１」を加算せずに（ステップＳＴ４０６の処理を行わずに）、処理を終了する。

また、ステップＳＴ４０１において、ＣＰＩＣＨシンボルの受信電力の平方根が閾値Ｐ未満の場合には（ステップＳＴ４０１：ＮＯ）、計数判定部１０５は、誤検出の可能性が大きいと推定し、変移回数および変移回数計数時間を更新しないように制御して、計数判定の処理を終了する。以上で、最大ドップラー周波数推定装置１００の動作の説明を終える。

ここで、実際に最大ドップラー周波数推定装置１００が運用される電波伝搬環境で予め試験を行い、最大ドップラー周波数推定装置１００単独で検出した最大ドップラー周波数と、ＧＰＳなどの位置情報を用いて求めた移動速度若しくは車速を用いて求めた最大ドップラー周波数などの、外部入力装置を用いて求めた最大ドップラー周波数との差が小さくなるように閾値Ｐを算出することができる。

図５は、閾値Ｐの設定方法を示す図である。図５は、試験により得られた閾値と、上記の最大ドップラー周波数の差分との関係の一例を示す図である。図５において、閾値は、ＣＰＩＣＨシンボルの受信電力と比較する閾値を示し、検出ｆｄの差分は、最大ドップラー周波数推定装置１００単独で検出した最大ドップラー周波数と、外部入力装置を用いて求めた最大ドップラー周波数との差を示す。図５の場合、最終的に予め設定される閾値Ｐの値は、−２２となる（矢印＃５０１）。

また、図６は、閾値Ｐの設定方法の他の例を示す図である。図６は、最大ドップラー周波数推定装置１００において検出したＣＰＩＣＨシンボルの受信電力と受信する遅延波の数とに応じて閾値を設定するものである。

無線通信装置が受信する遅延波の数を検出できるときに、検出した受信電力が０ｄＢ以上かつ−２０ｄＢ以下の場合において、検出した遅延波の数が１の場合は、閾値として−２２ｄＢを選択する。また、検出した遅延波数が２の場合は、閾値として−２２ｄＢを選択し、検出した遅延波数が３の場合は、閾値として−２０ｄＢを選択し、検出した遅延波化数が４以上の場合は、閾値として−１８ｄＢを選択する。

これにより、ＣＰＩＣＨシンボルの受信電力に加え、受信する遅延波の数が最大ドップラー周波数の検出精度に影響する場合でも、検出した受信電力および遅延波の数に応じて、フェージング変動以外の要因による影響の可能性を判定する閾値を設定することができる。これにより、最大ドップラー周波数の検出精度をさらに高めることができる。

図６に示すＣＰＩＣＨシンボルの受信電力と遅延波の数に応じた閾値の算出は以下の方法により算出することができる。まず、実際に本発明が運用される電波伝搬環境で試験を行い、ある閾値設定時において、検出したＣＰＩＣＨシンボルの受信電力、検出した遅延波の数、および検出した最大ドップラー周波数と外部入力装置を用いて求めた最大ドップラー周波数との差分を記録する。そして、記録した結果から、ＣＰＩＣＨシンボルの受信電力および遅延波の数毎に、最大ドップラー周波数推定装置１００単独で検出した最大ドップラー周波数と外部入力装置を用いて求めた最大ドップラー周波数数との差が小さくなるように、閾値を算出する。

図７は、閾値Ｐの設定時の試験結果の一例を示す図である。図７において、閾値は、ＣＰＩＣＨシンボルの受信電力と比較する閾値を示す。また、検出ｆｄの差分は、最大ドップラー周波数推定装置１００単独で検出した最大ドップラー周波数と外部入力装置を用いて求めた最大ドップラー周波数との差を示す。また、検出した遅延波数は、ＣＰＩＣＨシンボルの受信信号を求める際に受信した遅延波の数を示す。

図７において、検出した受信電力が−２０ｄＢ以上、かつ検出した遅延波数が１つであるケースに着目する。着目するケースに該当する試験結果は、「閾値−１８ｄＢを設定したときに、検出した受信電力が−２２ｄＢであり、検出ｆｄの差分が５０Ｈｚ」（矢印＃７０１）の場合と、「閾値−２０ｄＢを設定したときに検出した受信電力が−２２ｄＢであり、検出ｆｄの差分が３０Ｈｚ」（矢印＃７０２）の場合の二つである。この該当する二つ試験結果の中で、検出ｆｄの差分が小さい方を閾値とするため、図７の場合は、閾値として−２０ｄＢを選択する。遅延波数が２以上である場合も、同様の方法により閾値を選択する。

次に、図７において、検出した受信電力が０ｄＢ以上かつ−２０ｄＢ未満であり、かつ検出した遅延波数が１であるケースに着目する。着目するケースに該当する試験結果は、「閾値−２２ｄＢを設定したときに、検出した受信電力が−１５ｄＢ、検出ｆｄの差分が５０Ｈｚ」（矢印＃７０３）の場合のみである。よって、検出した受信電力が０ｄＢ以上かつ−２０ｄＢ未満であり、かつ検出し遅延波数が１であるケースの場合は、閾値として−２２ｄＢを選択する。遅延波数が２以上である場合も、同様の方法により閾値を選択する。

なお、受信電力の指標として、ＲＳＣＰ(Received Signal Code Power)、またはＩＳＣＰ（Interference Signal Code Power）の逆数、その他の指標であってもよい。また、上記の指標を２つ以上組み合わせてもよい。これにより、フェージング以外の影響の推定精度をさらに向上させ、最大ドップラー周波数の検出精度を高めることができる。

また、受信ＳＩＲ（Signal to Interference Ratio）に応じて閾値を設定してもよい。受信ＳＩＲが高い場合は、閾値ａを選択し、低い場合は閾値ｂを選択する。これにより、受信ＳＩＲが異なる場合でも、フェージング変動以外の要因による影響の可能性の判定する閾値を変えることで、フェージング変動以外の要因による影響の可能性の判定精度を高めることができる。なお、上記の閾値は、図５を用いて説明した閾値の設定方法と同様に、検出ｆｄの差分が小さくなるように決めてもよい。

このように、本実施の形態によれば、ＣＰＩＣＨシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合に、フェージング変動による最大ドップラー周波数の推定精度を向上させることができる。

（実施の形態２）
図８は、本発明の実施の形態２に係る最大ドップラー周波数推定装置８００の構成を示すブロック図である。最大ドップラー周波数推定装置８００は、例えば、無線受信装置に搭載される。

図８に示す最大ドップラー周波数推定装置８００は、図２に示す実施の形態１に係る最大ドップラー周波数推定装置１００に対して、ＣＰＩＣＨシンボルの受信パワー測定部１０４を除き、ＣＰＩＣＨシンボルのＩＱ成分測定部８０１を追加し、計数判定部１０５の代わりに計数判定部８０２を有する。なお、図８において、図２と同一構成である部分には同一の符号を付してその説明を省略する。

最大ドップラー周波数推定装置８００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、変移検出部１０６と、変移回数計数部１０７と、ｆｄ推定値検出部１０８と、ＣＰＩＣＨシンボルのＩＱ成分測定部８０１と、計数判定部８０２とから主に構成される。以下に、本実施の形態における上記の実施の形態１とは異なる各構成について詳細に説明する。

ＣＰＩＣＨシンボルのデータ抽出部１０３は、無線受信部１０２から入力したＩ成分およびＱ成分に対し、逆拡散を施しＣＰＩＣＨシンボルのＩ成分およびＱ成分を抽出する。また、ＣＰＩＣＨシンボルのデータ抽出部１０３は、抽出したＣＰＩＣＨシンボルのＩ成分およびＱ成分をＣＰＩＣＨシンボルのＩＱ成分測定部８０１および変移検出部１０６へ出力する。

ＣＰＩＣＨシンボルのＩＱ成分測定部８０１は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分から、Ｉ−Ｑ平面上の受信信号点の座標を求める。そして、ＣＰＩＣＨシンボルのＩＱ成分測定部８０１は、求めた座標を計数判定部８０２へ出力する。

計数判定部８０２は、ＣＰＩＣＨシンボルのＩＱ成分測定部８０１から入力した座標と閾値とを比較してフェージング変動以外の要因による影響の可能性を推定する。具体的には、計数判定部８０２は、受信信号点の座標とＩ軸およびＱ軸との距離を求める。また、計数判定部８０２は、求めた各々の距離が所定の範囲内である場合にはフェージング変動以外の要因による影響の可能性が大きいと判定し、求めた各々の距離が所定の範囲外である場合にはフェージング変動以外の要因による影響の可能性は小さいと判定する。また、計数判定部８０２は、フェージング変動以外の要因による影響の可能性が小さいと判定した場合には、変移を計数するように変移回数計数部１０７を制御し、フェージング変動以外の要因による影響の可能性が大きいと判定した場合には、変移を計数しないように変移回数計数部１０７を制御する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部８０２の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。以上で、最大ドップラー周波数推定装置８００の構成の説明を終える。

次に、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置について、図９を用いて説明する。図９は、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を説明する図である。

図９において、ＣＰＩＣＨシンボルの現在の受信信号点Ｘ（Ｉ、Ｑ）のＩ成分の値は、Ｑ軸からの距離を示し、Ｉ成分の絶対値により求めることができる。また、ＣＰＩＣＨシンボルの現在の受信信号点Ｘ（Ｉ、Ｑ）のＱ成分の値は、Ｉ軸からの距離を示し、Ｑ成分の絶対値により求めることができる。

ここで、受信信号のＩ成分の絶対値が予め設定した閾値ａ以上の場合には、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が小さく、誤検出の可能性は小さいと推定できる。逆に、Ｉ成分の絶対値が予め設定した閾値ａ未満の場合には、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が大きく、誤検出の可能性は大きいと推定できる。Ｑ成分の絶対値も同様に予め設定した閾値ｂとの大小を判定し、フェージング変動以外の要因による影響の可能性を推定する。

これを用いて、ＣＰＩＣＨシンボルの受信信号のＩ成分の値が横縞の網掛け範囲内の場合（｜Ｉ｜＜ａ）は、Ｉ成分の変移検出結果を変移回数の計数対象外とし、横縞の網掛け範囲外の場合（｜Ｉ｜≧ａ）は、Ｉ成分の変移検出結果を変移回数の計数対象とする。

Ｑ成分についても同様に、ＣＰＩＣＨシンボルの受信信号のＱ成分の値が縦縞の網掛けの範囲内の場合（｜Ｑ｜＜ｂ）は、Ｑ成分の変移検出結果を変移回数の計数対象外とし、縦縞の網掛け範囲外の場合（｜Ｑ｜≧ｂ）は、Ｑ成分の変移検出結果を変移回数の計数対象とする。

図９の場合、計数判定部８０２は、現在の受信信号点Ｘ（Ｉ、Ｑ）のＩ成分の変移検出結果を変移回数の計数対象とし、現在の受信信号点Ｘ（Ｉ、Ｑ）のＱ成分の変移検出結果は変移回数の計数対象外とするように、変移回数計数部１０７を制御する。

次に、最大ドップラー周波数推定装置８００の動作について、図１０を用いて説明する。図１０は、最大ドップラー周波数推定装置８００の動作の一部を示すフロー図である。

最初に、計数判定部８０２は、受信信号点のＩ成分の絶対値を求め、求めたＩ成分の絶対値が予め設定した閾値ａ以上か否かを判定する（ステップＳＴ１００１）。

受信信号点のＩ成分の絶対値が予め設定した閾値ａ未満の場合には（ステップＳＴ１００１：ＮＯ）、計数判定部８０２は、受信信号点のＱ成分の絶対値を求め、求めたＱ成分の絶対値が予め設定した閾値ｂ以上か否かを判定する（ステップＳＴ１００２）。

受信信号点のＱ成分の絶対値が予め設定した閾値ｂ未満である場合には（ステップＳＴ１００２：ＮＯ）、処理を終了する。

一方、ステップＳＴ１００１において、受信信号点のＩ成分の絶対値が予め設定した閾値ａ以上の場合には（ステップＳＴ１００１：ＹＥＳ）、ステップＳＴ１００２の処理を行わずにステップＳＴ１００３に処理を進める。

また、ステップＳＴ１００２において、受信信号点のＱ成分の絶対値が予め設定した閾値ｂ以上の場合には（ステップＳＴ１００２：ＹＥＳ）、計数判定部８０２は、変移回数計数時間に１を加算する（変移回数計数時間＋１）（ステップＳＴ１００３）。

次に、計数判定部８０２は、受信信号点のＩ成分の絶対値が予め設定した閾値ａ以上か否かを判定する（ステップＳＴ１００４）。

受信信号点のＩ成分の絶対値が予め設定した閾値ａ未満の場合には（ステップＳＴ１００４：ＮＯ）、ステップＳＴ１００５及びステップＳＴ１００６の処理を行わずに、後述するステップＳＴ１００７に処理を進める。

一方、受信信号点のＩ成分の絶対値が予め設定した閾値ａ以上の場合には（ステップＳＴ１００４：ＹＥＳ）、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ１００５）。

Ｉ成分の符号変化がない場合には（ステップＳＴ１００５：ＮＯ）、ステップＳＴ１００６の処理を行わずに、後述するステップＳＴ１００７に処理を進める。

一方、Ｉ成分の符号変化がある場合には（ステップＳＴ１００５：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ１００６）。

次に、計数判定部８０２は、受信信号点のＱ成分の絶対値が予め設定した閾値ｂ以上か否かを判定する（ステップＳＴ１００７）。

受信信号点のＱ成分の絶対値が予め設定した閾値ｂ未満の場合には（ステップＳＴ１００７：ＮＯ）、ステップＳＴ１００８及びステップＳＴ１００９の処理を行わずに処理を終了する。

一方、受信信号点のＱ成分の絶対値が予め設定した閾値ｂ以上の場合には（ステップＳＴ１００７：ＹＥＳ）、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ１００８）。

Ｑ成分の符号変化がない場合には（ステップＳＴ１００８：ＮＯ）、ステップＳＴ１００９の処理を行わずに処理を終了する。

一方、Ｑ成分の符号変化がある場合には（ステップＳＴ１００８：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算し（変移回数＋１）（ステップＳＴ１００９）、処理を終了する。

このように、本実施の形態によれば、ＣＰＩＣＨシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合に、フェージング変動による最大ドップラー周波数の推定精度を向上させることができる。

（実施の形態３）
図１１は、本発明の実施の形態３に係る最大ドップラー周波数推定装置１１００の構成を示すブロック図である。最大ドップラー周波数推定装置１１００は、例えば、無線受信装置に搭載される。

図１１に示す最大ドップラー周波数推定装置１１００は、図２に示す実施の形態１に係る最大ドップラー周波数推定装置１００に対して、ＣＰＩＣＨシンボルの受信パワー測定部１０４を除き、ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１およびメモリ・変移量算出部１１０２を追加し、計数判定部１０５の代わりに計数判定部１１０３を有する。なお、図１１において、図２と同一構成である部分には同一の符号を付してその説明を省略する。

最大ドップラー周波数推定装置１１００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、変移検出部１０６と、変移回数計数部１０７と、ｆｄ推定値検出部１０８と、ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１と、メモリ・変移量算出部１１０２と、計数判定部１１０３とから主に構成される。以下に、本実施の形態における上記の実施の形態１とは異なる各構成について詳細に説明する。

ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分から、Ｉ−Ｑ平面上の受信信号点の座標を求める。そして、ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１は、求めた座標をメモリ・変移量算出部１１０２へ出力する。

メモリ・変移量算出部１１０２は、ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１から入力した、座標のＩ成分およびＱ成分の値を保存する。また、メモリ・変移量算出部１１０２は、ＣＰＩＣＨシンボルのＩＱ成分測定部１１０１から入力した、現在のＣＰＩＣＨシンボルの受信信号点のＩ成分およびＱ成分と、保存している直前のＣＰＩＣＨシンボルの受信信号点のＩ成分およびＱ成分とに基づいて、Ｉ成分およびＱ成分のそれぞれの変移の大きさ（以下、「変移量」と記載する）を検出する。そして、メモリ・変移量算出部１１０２は、検出した変移量を計数判定部１１０３へ出力する。

計数判定部１１０３は、メモリ・変移量算出部１１０２から入力した変移量と閾値とを比較してフェージング変動以外の要因による影響の可能性を推定する。具体的には、計数判定部１１０３は、変移量が閾値以上の場合にはフェージング変動以外の要因による影響の可能性が小さいと判定し、変移量が閾値未満の場合にはフェージング変動以外の要因による影響の可能性が大きいと判定する。また、計数判定部１１０３は、フェージング変動以外の要因による影響の可能性が小さいと判定した場合には、変移を計数するように変移回数計数部１０７を制御し、フェージング変動以外の要因による影響の可能性が大きいと判定した場合には、変移を計数しないように変移回数計数部１０７を制御する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部１１０３の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。以上で、最大ドップラー周波数推定装置１１００の構成の説明を終える。

次に、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置について、図１２を用いて説明する。図１２は、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を説明する図である。

図１２において、Ｉ_ｉは、ＣＰＩＣＨシンボルの現在（時刻ｉ）の受信信号点ＸのＩ成分の値を示し、Ｑ_ｉは、Ｑ成分の値を示す。また、図１２において、Ｉ_ｉ−１は、直前（時刻ｉ−１）のＣＰＩＣＨシンボルの受信信号点ＹのＩ成分の値を示し、Ｑ_ｉ−１は、Ｑ成分の値を示す。この場合、Ｉ成分の変移量はＩ_ｉ−Ｉ_ｉ−１の絶対値で示され、Ｑ成分の変移量はＱ_ｉ−Ｑ_ｉ−１の絶対値で示される。

ここで、受信信号点のＩ成分の変移量が、予め設定した閾値ａ以上の場合は、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が小さく、誤検出の可能性が小さいと推定できる。逆に、受信信号点のＩ成分の変移量が予め設定した閾値ａ未満の場合は、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が大きく、誤検出の可能性が大きいと推定できる。Ｑ成分の変移量も同様に予め設定した閾値ｂとの大小を判定し、フェージング変動以外の要因による影響の可能性を推定する。

これを用いて、ＣＰＩＣＨシンボルの受信信号のＩ成分の変移量が横縞の網掛け範囲内の場合（Ｉ成分の変移量（Ｉ_ｉ−Ｉ_ｉ−１）＜ａ）は、Ｉ成分の変移検出結果を計数対象外とし、ＣＰＩＣＨシンボルの受信信号のＩ成分の変移量が横縞の網掛け範囲外の場合（Ｉ成分の変移量（Ｉ_ｉ−Ｉ_ｉ−１）≧ａ）は、Ｉ成分の変移検出結果を変移回数の計数対象とする。Ｑ成分についても同様に、Ｑ成分の変移量が縦縞の網掛けの範囲内の場合（Ｑ成分の変移量（Ｑ_ｉ−Ｑ_ｉ−１）＜ｂ）は、Ｑ成分の変移検出結果を変移回数の計数対象外とし、Ｑ成分の変移量が縦縞の網掛け範囲外のとき（Ｑ成分の変移量（Ｑ_ｉ−Ｑ_ｉ−１）≧ｂ）は、Ｑ成分の変移検出結果を変移回数の計数対象とする。

図１２の場合、計数判定部１１０３は、現在の受信信号点Ｘ（Ｉ_ｉ、Ｑ_ｉ）のＩ成分の変移検出結果を変移回数の計数対象とし、現在の受信信号点Ｘ（Ｉ_ｉ、Ｑ_ｉ）のＱ成分の変移検出結果は変移回数の計数対象外とするように、変移回数計数部１０７を制御する。

次に、最大ドップラー周波数推定装置１１００の動作について、図１３を用いて説明する。図１３は、最大ドップラー周波数推定装置１１００の動作の一部を示すフロー図である。

最初に、計数判定部１１０３は、Ｉ成分の変移量を求め、求めたＩ成分の変移量が閾値ａ以上であるか否かを判定する（ステップＳＴ１３０１）。

Ｉ成分の変移量が閾値ａ以上である場合には（ステップＳＴ１３０１：ＹＥＳ）、ステップＳＴ１３０２の処理を行わずに、ステップＳＴ１３０３に処理を進める。

一方、Ｉ成分の変移量が閾値ａ未満である場合には（ステップＳＴ１３０１：ＮＯ）、計数判定部１１０３は、Ｑ成分の変移量を求め、求めたＱ成分の変移量が閾値ｂ以上であるか否かを判定する（ステップＳＴ１３０２）。

Ｑ成分の変移量が閾値ｂ未満である場合には（ステップＳＴ１３０２：ＮＯ）、ステップＳＴ１３０３〜ステップＳＴ１３０９の処理を行わずに、処理を終了する。

一方、Ｑ成分の変移量が閾値ｂ以上である場合には（ステップＳＴ１３０２：ＹＥＳ）、計数判定部１１０３は、変移回数計数時間に１を加算する（変移回数計数時間＋１）（ステップＳＴ１３０３）。

次に、計数判定部１１０３は、Ｉ成分の変移量が閾値ａ以上であるか否かを判定する（ステップＳＴ１３０４）。

Ｉ成分の変移量が閾値ａ未満の場合には（ステップＳＴ１３０４：ＮＯ）、ステップＳＴ１３０５及びステップＳＴ１３０６の処理を行わずに、後述するステップＳＴ１３０７に処理を進める。

一方、Ｉ成分の変移量が閾値ａ以上の場合には（ステップＳＴ１３０４：ＹＥＳ）、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ１３０５）。

Ｉ成分の符号変化がない場合には（ステップＳＴ１３０５：ＮＯ）、ステップＳＴ１３０６の処理を行わずに、後述するステップＳＴ１３０７に処理を進める。

一方、Ｉ成分の符号変化がある場合には（ステップＳＴ１３０５：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ１３０６）。

次に、計数判定部１１０３は、Ｑ成分の変移量が閾値ｂ以上であるか否かを判定する（ステップＳＴ１３０７）。

Ｑ成分の変移量が閾値ｂ未満である場合には（ステップＳＴ１３０７：ＮＯ）、ステップＳＴ１３０８及びステップＳＴ１３０９の処理を行わずに、処理を終了する。

一方、Ｑ成分の変移量が閾値ｂ以上である場合には（ステップＳＴ１３０７：ＹＥＳ）、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ１３０８）。

Ｑ成分の符号変化がない場合には（ステップＳＴ１３０８：ＮＯ）、ステップＳＴ１３０９の処理を行わずに、処理を終了する。

一方、Ｑ成分の符号変化がある場合には（ステップＳＴ１３０８：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算し（変移回数＋１）（ステップＳＴ１３０９）、処理を終了する。

このように、本実施の形態によれば、ＣＰＩＣＨシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合に、フェージング変動の影響のみを受けた最大ドップラー周波数の推定精度を向上させることができる。

（実施の形態４）
図１４は、本発明の実施の形態４に係る最大ドップラー周波数推定装置１４００の構成を示すブロック図である。最大ドップラー周波数推定装置１４００は、例えば、無線受信装置に搭載される。

図１４に示す最大ドップラー周波数推定装置１４００は、図２に示す実施の形態１に係る最大ドップラー周波数推定装置１００に対して、ＣＰＩＣＨシンボルの受信パワー測定部１０４を除き、ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１及びメモリ・シンボルの平均ＩＱ成分算出部１４０２を追加し、計数判定部１０５の代わりに計数判定部１４０３を有する。なお、図１４において、図２と同一構成である部分には同一の符号を付してその説明を省略する。

最大ドップラー周波数推定装置１４００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、変移検出部１０６と、変移回数計数部１０７と、ｆｄ推定値検出部１０８と、ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１と、メモリ・シンボルの平均ＩＱ成分算出部１４０２と、計数判定部１４０３とから主に構成される。以下に、本実施の形態における上記の実施の形態１とは異なる各構成について詳細に説明する。

ＣＰＩＣＨシンボルのデータ抽出部１０３は、無線受信部１０２から入力したＩ成分およびＱ成分に対し、逆拡散を施しＣＰＩＣＨシンボルのＩ成分およびＱ成分を抽出する。また、ＣＰＩＣＨシンボルのデータ抽出部１０３は、抽出したＣＰＩＣＨシンボルのＩ成分およびＱ成分をＣＰＩＣＨシンボルのＩＱ成分測定部１４０１および変移検出部１０６へ出力する。

ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分から、Ｉ−Ｑ平面上の受信信号点の座標を求める。そして、ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１は、求めた座標をメモリ・シンボルの平均ＩＱ成分算出部１４０２へ出力する。

メモリ・シンボルの平均ＩＱ成分算出部１４０２は、ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１から入力した、直前（時刻ｉ−１）のＣＰＩＣＨシンボルの受信信号点の座標のＩ成分およびＱ成分を保存する。また、メモリ・シンボルの平均ＩＱ成分算出部１４０２は、現在（時刻ｉ）において、ＣＰＩＣＨシンボルのＩＱ成分測定部１４０１から入力した、現在（時刻ｉ）のＣＰＩＣＨシンボルの受信信号点の座標のＩ成分およびＱ成分に、保存している直前（時刻ｉ−１）のＣＰＩＣＨシンボルの受信信号点の座標のＩ成分およびＱ成分を同相加算し、加算結果のＩ成分およびＱ成分に０．５を乗算する。そして、メモリ・シンボルの平均ＩＱ成分算出部１４０２は、乗算結果を計数判定部１４０３へ出力する。なお、２シンボルの平均Ｉ成分およびＱ成分を算出するために０．５を乗算したが、０．５を乗算する代わりに、平均Ｉ成分および平均Ｑ成分に対する閾値を２倍にしてもよい。

計数判定部１４０３は、メモリ・シンボルの平均ＩＱ成分算出部１４０２から入力した乗算結果と閾値とを比較してフェージング変動以外の要因による影響の可能性を推定する。具体的には、計数判定部１４０３は、Ｉ成分またはＱ成分の乗算結果が閾値以上の場合にはフェージング変動以外の要因による影響の可能性が小さいと判定し、Ｉ成分またはＱ成分の乗算結果が閾値未満の場合にはフェージング変動以外の要因による影響の可能性が大きいと判定する。また、計数判定部１４０３は、フェージング変動以外の要因による影響の可能性が小さいと判定した場合には、変移を計数するように変移回数計数部１０７を制御し、フェージング変動以外の要因による影響の可能性が大きいと判定した場合には、変移を計数しないように変移回数計数部１０７を制御する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部１４０３の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。以上で、最大ドップラー周波数推定装置１４００の構成の説明を終える。

次に、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置について、図１５を用いて説明する。図１５は、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を説明する図である。

ＣＰＩＣＨシンボルの現在（時刻ｉ）の受信信号点の座標Ｘを（Ｉ_ｉ、Ｑ_ｉ）、直前（時刻ｉ−１）の受信信号点の座標Ｙを（Ｉ_ｉ−１、Ｑ_ｉ−１）とすると、現在の受信信号点と直前の受信信号点との平均信号点の座標Ｚにおいて、Ｉ成分（平均Ｉ成分）は「Ｉ＝（Ｉ_ｉ−１＋Ｉ_ｉ）／２」により求めることができ、Ｑ成分（平均Ｑ成分）は「Ｑ＝（Ｑ_ｉ−１＋Ｑ_ｉ）／２」により求めることができる。

ここで、平均Ｉ成分が予め設定した閾値ａ以上の場合は、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が小さく、誤検出の可能性が小さいと推定できる。逆に、平均Ｉ成分が予め設定した閾値ａ未満の場合は、Ｉ成分におけるフェージング変動以外の要因による影響の可能性が大きく、誤検出の可能性が大きいと推定できる。平均Ｑ成分も同様に予め設定した閾値ｂとの大小を判定し、フェージング変動以外の要因による影響の可能性を推定する。

これを用いて、ＣＰＩＣＨシンボルの平均Ｉ成分の値Ｉが横縞の網掛け範囲内の場合（｜Ｉ｜＜ａ）は、Ｉ成分の変移検出結果を変移回数の計数対象外とし、横縞の網掛け範囲外の場合（｜Ｉ｜≧ａ）は、Ｉ成分の変移検出結果を変移回数の計数対象とする。Ｑ成分についても同様に、平均Ｑ成分の値Ｑが縦縞の網掛けの範囲内の場合（｜Ｑ｜＜ｂ）は、Ｑ成分の変移検出結果を変移回数の計数対象外とし、縦縞の網掛け範囲外の場合（｜Ｑ｜≧ｂ）はＱ成分の変移検出結果を変移回数の計数対象とする。

図１５の場合、計数判定部１４０３は、現在の受信信号点Ｘ（Ｉ_ｉ、Ｑ_ｉ）のＩ成分およびＱ成分の変移検出結果は変移回数の計数対象外とするように、変移回数計数部１０７を制御する。

なお、本実施の形態では、現在（時刻ｉ）のＣＰＩＣＨシンボルと直前（時刻ｉ−１）のＣＰＩＣＨシンボルとの各１シンボルを用い、受信信号の平均Ｉ成分および平均Ｑ成分を求める区間の長さである平均区間長を２シンボルとしたが、２シンボルでなくてもよい。

直前（時刻ｉ−１）と、現在（時刻ｉ）と、直後（時刻ｉ＋１）との３シンボルの受信信号の平均Ｉ成分および平均Ｑ成分を算出し、算出した平均Ｉ成分及び平均Ｑ成分を用いて、現在（時刻ｉ）検出した変移検出結果を変移回数の計数対象とするか否かを判定してもよい。現在を含まない、過去と未来、即ち直前（時刻ｉ−１）と直後（時刻ｉ＋１）とに受信したＣＰＩＣＨシンボルの平均値を用いて同様に判定してもよい。また、現在の１シンボルと過去の２シンボルと未来の２シンボル（時刻ｉ−２、ｉ−１、ｉ、ｉ＋１、ｉ＋２の５シンボル全て）を用いて同様の判定をしてもよい。また、時刻ｉ−２、ｉ−１、ｉ、ｉ＋１、ｉ＋２の５シンボルの内からいずれか４シンボル、または３シンボル、または２シンボルを用いて同様の判定をしてもよい。過去と未来の数を必ずしも等しくする必要はない。

次に、最大ドップラー周波数推定装置１４００の動作について、図１６を用いて説明する。図１６は、最大ドップラー周波数推定装置１４００の動作の一部を示すフロー図である。

最初に、計数判定部１４０３は、平均Ｉ成分を求め、求めた平均Ｉ成分の絶対値が閾値ａ以上であるか否かを判定する（ステップＳＴ１６０１）。

平均Ｉ成分の絶対値が閾値ａ以上の場合には（ステップＳＴ１６０１：ＹＥＳ）、ステップＳＴ１６０２の処理を行わずに、後述するステップＳＴ１６０３に処理を進める。

一方、平均Ｉ成分の絶対値が閾値ａ未満の場合には（ステップＳＴ１６０１：ＮＯ）、計数判定部１４０３は、平均Ｑ成分を求め、求めた平均Ｑ成分の絶対値が閾値ｂ以上であるか否かを判定する（ステップＳＴ１６０２）。

平均Ｑ成分の絶対値が閾値ｂ未満の場合には（ステップＳＴ１６０２：ＮＯ）、ステップＳＴ１６０３〜ステップＳＴ１６０９の処理を行わずに、処理を終了する。

一方、平均Ｑ成分の絶対値が閾値ｂ以上の場合には（ステップＳＴ１６０２：ＹＥＳ）、計数判定部１４０３は、変移回数計数時間に１を加算する（変移回数計数時間＋１）（ステップＳＴ１６０３）。

次に、計数判定部１４０３は、平均Ｉ成分の絶対値が閾値ａ以上であるか否かを判定する（ステップＳＴ１６０４）。

平均Ｉ成分の絶対値が閾値ａ未満の場合には（ステップＳＴ１６０４：ＮＯ）、ステップＳＴ１６０５およびステップＳＴ１６０６の処理を行わずに、後述するステップＳＴ１６０７に処理を進める。

一方、平均Ｉ成分の絶対値が閾値ａ以上の場合には（ステップＳＴ１６０４：ＹＥＳ）、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ１６０５）。

Ｉ成分の符号変化がない場合には（ステップＳＴ１６０５：ＮＯ）、ステップＳＴ１６０６の処理を行わずに、後述するステップＳＴ１６０７に処理を進める。

一方、Ｉ成分の符号変化がある場合には（ステップＳＴ１６０５：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ１６０６）。

次に、計数判定部１４０３は、平均Ｑ成分の絶対値が閾値ｂ以上であるか否かを判定する（ステップＳＴ１６０７）。

平均Ｑ成分の絶対値が閾値ｂ未満の場合には（ステップＳＴ１６０７：ＮＯ）、ステップＳＴ１６０８およびステップＳＴ１６０９の処理を行わずに、処理を終了する。

一方、平均Ｑ成分の絶対値が閾値ｂ以上の場合には（ステップＳＴ１６０７：ＹＥＳ）、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ１６０８）。

Ｑ成分の符号変化がない場合には（ステップＳＴ１６０８：ＮＯ）、ステップＳＴ１６０９の処理を行わずに、処理を終了する。

一方、Ｑ成分の符号変化がある場合には（ステップＳＴ１６０８：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算し（変移回数＋１）（ステップＳＴ１６０９）、処理を終了する。

なお、本実施の形態では、ＣＰＩＣＨシンボルの受信信号の平均Ｉ成分および平均Ｑ成分を用いてフェージング変動以外の要因による影響の可能性を推定する方法を示したが、平均Ｉ成分および平均Ｑ成分から算出した受信電力を用いて、実施の形態１と同様の方法でフェージング変動以外の要因による影響の可能性を推定してもよい。

因みに、上記の実施の形態２の場合、フェージング変動が低速であり、フェージング変動によるＣＰＩＣＨシンボルの受信信号の変移量が小さい場合、受信信号の変移量に対するフェージング変動以外の要因による影響の可能性が大きく、検出した変移回数がフェージング変動以外の要因による影響の計数である可能性が大きい。

一方、本実施の形態の構成によりＣＰＩＣＨシンボルの受信信号を平均することにより、フェージング変動が低速な場合において、フェージング変動以外の要因による影響を平滑化し、その影響を小さくでき、受信信号のばらつきを小さくすることができる。このばらつきが小さい受信信号を用いることで、平均しない受信信号を用いる場合に比べ、計数判定に用いる閾値の設定を容易に行うことができる。

このように、本実施の形態によれば、上記の実施の形態１の効果に加えて、フェージング変動が低速な場合において、計数判定に用いる閾値を容易に設定することができる。

（実施の形態５）
図１７は、本発明の実施の形態５に係る最大ドップラー周波数推定装置１７００の構成を示すブロック図である。最大ドップラー周波数推定装置１７００は、例えば、無線受信装置に搭載される。

図１７に示す最大ドップラー周波数推定装置１７００は、図２に示す実施の形態１に係る最大ドップラー周波数推定装置１００に対して、ＣＰＩＣＨシンボルの受信パワー測定部１０４を除き、ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１と、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２と、第２のメモリ・シンボルの平均ＩＱ成分算出部１７０３とを追加し、計数判定部１０５の代わりに計数判定部１７０４を有する。なお、図１７において、図２と同一構成である部分には同一の符号を付してその説明を省略する。

最大ドップラー周波数推定装置１７００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、変移検出部１０６と、変移回数計数部１０７と、ｆｄ推定値検出部１０８と、ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１と、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２と、第２のメモリ・シンボルの平均ＩＱ成分算出部１７０３と、計数判定部１７０４とから主に構成される。以下に、本実施の形態における上記の実施の形態１とは異なる各構成について詳細に説明する。

ＣＰＩＣＨシンボルのデータ抽出部１０３は、無線受信部１０２から入力したＩ成分およびＱ成分に対しから、逆拡散を施しＣＰＩＣＨシンボルのＩ成分およびＱ成分を抽出する。また、ＣＰＩＣＨシンボルのデータ抽出部１０３は、抽出したＣＰＩＣＨシンボルのＩ成分およびＱ成分をＣＰＩＣＨシンボルのＩＱ成分測定部１７０１および変移検出部１０６へ出力する。

ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分から、Ｉ−Ｑ平面上の受信信号点の座標を求める。そして、ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１は、求めた座標を第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２および第２のメモリ・シンボルの平均ＩＱ成分算出部１７０３へ出力する。

第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２は、ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１から入力した座標より、所定の平均区間長における平均Ｉ成分および平均Ｑ成分を求める。そして、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２は、求めた平均Ｉ成分および平均Ｑ成分を計数判定部１７０４へ出力する。

第２のメモリ・シンボルの平均ＩＱ成分算出部１７０３は、ＣＰＩＣＨシンボルのＩＱ成分測定部１７０１から入力した座標より、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２の平均区間長とは異なる平均区間長における平均Ｉ成分および平均Ｑ成分を求める。そして、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２は、求めた平均Ｉ成分および平均Ｑ成分を計数判定部１７０４へ出力する。

計数判定部１７０４は、第１のメモリ・シンボルの平均ＩＱ成分算出部１７０２から入力した平均Ｉ成分および平均Ｑ成分、および第２のメモリ・シンボルの平均ＩＱ成分算出部１７０３から入力した平均Ｉ成分および平均Ｑ成分と、閾値とを各々比較してフェージング変動以外の要因による影響の可能性を推定する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部１７０４の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。以上で、最大ドップラー周波数推定装置１７００の構成の説明を終える。

次に、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置について、図１８を用いて説明する。図１８は、本実施の形態におけるＣＰＩＣＨシンボルの受信信号点のＩ−Ｑ平面における配置を説明する図である。

図１８の平均信号点１（Ｉ_１、Ｑ_１）と平均信号点２（Ｉ_２、Ｑ_２）は、各々異なる平均区間長で平均した平均信号点である。図１８において、Ｉ_１は、平均信号点１における平均Ｉ成分であり、Ｑ_１は、平均信号点１における平均Ｑ成分である。また、図１８において、Ｉ_２は、平均信号点２における平均Ｉ成分であり、Ｑ_２は、平均信号点２における平均Ｑ成分である。

ここで、平均信号点１における平均Ｉ成分が予め設定した閾値ａ１以上の場合は、平均区間長１の平均Ｉ成分におけるフェージング変動以外の要因による影響の可能性が小さく、誤検出の可能性が小さいと推定できる。逆に、平均信号点１における平均Ｉ成分が予め設定した閾値ａ１未満の場合は、平均区間長１の平均Ｉ成分におけるフェージング変動以外の要因による影響の可能性が大きく、誤検出の可能性が大きいと推定できる。平均信号点１における平均Ｑ成分も同様に、予め設定した閾値ｂ１との大小を判定し、フェージング変動以外の要因による影響の可能性を推定する。同様に、平均信号点２における平均Ｉ成分および平均Ｑ成分について、閾値ａ２および閾値ｂ２と大小判定を行い、フェージング変動以外の要因による影響の可能性を推定する。

これを用いて、ＣＰＩＣＨシンボルの平均信号点１の平均Ｉ成分の値Ｉ_１が横縞の網掛け範囲内のとき（｜Ｉ_１｜＜ａ１）、またはＣＰＩＣＨシンボルの平均信号点２の平均Ｉ成分の値Ｉ_２が右斜線縞の網掛け範囲内のとき（｜Ｉ_２｜＜ａ２）は、Ｉ成分の変移検出結果を変移回数の計数対象外とする。また、ＣＰＩＣＨシンボルの平均信号点１の平均Ｉ成分の値Ｉ_１が横縞の網掛け範囲外のとき（｜Ｉ_１｜≧ａ１）、かつＣＰＩＣＨシンボルの平均信号点２の平均Ｉ成分の値Ｉ_２が右斜線縞の網掛け範囲外のとき（｜Ｉ_２｜≧ａ２）は、Ｉ成分の変移検出結果を変移回数の計数対象とする。Ｑ成分についても同様に、ＣＰＩＣＨシンボルの平均信号点１の平均Ｑ成分の値Ｑ_１が縦縞の網掛けの範囲内のとき（｜Ｑ_１｜＜ｂ１）、またはＣＰＩＣＨシンボルの平均信号点２の平均Ｑ成分の値Ｑ_２が左斜線縞の網掛けの範囲内のとき（｜Ｑ_２｜＜ｂ２）は、Ｑ成分の変移検出結果を変移回数の計数対象外とする。また、ＣＰＩＣＨシンボルの平均信号点１の平均Ｑ成分の値Ｑ_１が縦縞の網掛け範囲外のとき（｜Ｑ_１｜≧ｂ１）、かつＣＰＩＣＨシンボルの平均信号点２の平均Ｑ成分の値Ｑ_２が左斜線縞の網掛け範囲外のとき（｜Ｑ_２｜≧ｂ２）は、Ｑ成分の変移検出結果を変移回数の計数対象とする。

図１８の場合、計数判定部１７０４は、平均信号点１（Ｉ_１、Ｑ_１）および平均信号点２（Ｉ_２、Ｑ_２）のＩ成分およびＱ成分の変移検出結果は変移回数の計数対象外とするように、変移回数計数部１０７を制御する。

次に、最大ドップラー周波数推定装置１７００の動作について、図１９および図２０を用いて説明する。図１９は、最大ドップラー周波数推定装置１７００の動作の一部０を示すフロー図である。また、図２０は、閾値を示す図である。

最初に、計数判定部１７０４は、平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１（図２０において、ａ１＝−１８）以上か否かを判定する（ステップＳＴ１９０１）。

平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１未満の場合には（ステップＳＴ１９０１：ＮＯ）、計数判定部１７０４は、平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１（図２０において、ｂ１＝−１８）以上か否かを判定する（ステップＳＴ１９０２）。

平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１未満の場合には（ステップＳＴ１９０２：ＮＯ）、計数判定部１７０４は、平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２（図２０において、ａ２＝−２５）以上か否かを判定する（ステップＳＴ１９０３）。

平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２未満の場合には（ステップＳＴ１９０３：ＮＯ）、計数判定部１７０４は、平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２（図２０において、ｂ２＝−２５）以上か否かを判定する（ステップＳＴ１９０４）。

平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２未満の場合には（ステップＳＴ１９０４：ＮＯ）、処理を終了する。

一方、ステップＳＴ１９０１において、平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１以上の場合（ステップＳＴ１９０１：ＹＥＳ）、ステップＳＴ１９０２において、平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１以上の場合（ステップＳＴ１９０２：ＹＥＳ）、ステップＳＴ１９０３において、平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２以上の場合（ステップＳＴ１９０３：ＹＥＳ）、およびステップＳＴ１９０４において、平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２以上の場合（ステップＳＴ１９０４：ＹＥＳ）の少なくとも何れか１つに該当する場合には、計数判定部１７０４は、変移回数計数時間に１を加算する（変移回数計数時間＋１）（ステップＳＴ１９０５）。

次に、計数判定部１７０４は、平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１（図２０において、ａ１＝−１８）以上か否かを判定する（ステップＳＴ１９０６）。

平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１以上の場合には（ステップＳＴ１９０６：ＹＥＳ）、ステップＳＴ１９０７の処理を行わずに、ステップＳＴ１９０８に処理を進める。

一方、平均ＩＱ成分１の平均Ｉ成分の絶対値が閾値ａ１未満の場合には（ステップＳＴ１９０６：ＮＯ）、計数判定部１７０４は、平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２（図２０において、ａ２＝−２５）以上か否かを判定する（ステップＳＴ１９０７）。

平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２未満の場合には（ステップＳＴ１９０７：ＮＯ）、ステップＳＴ１９０８およびステップＳＴ１９０９の処理を行わずに、後述するステップＳＴ１９１０に処理を進める。

一方、平均ＩＱ成分２の平均Ｉ成分の絶対値が閾値ａ２以上の場合には（ステップＳＴ１９０７：ＹＥＳ）、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ１９０８）。

Ｉ成分の符号変化がない場合には（ステップＳＴ１９０８：ＮＯ）、ステップＳＴ１９０９の処理を行わずに、後述するステップＳＴ１９１０に処理を進める。

一方、Ｉ成分の符号変化がある場合には（ステップＳＴ１９０８：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ１９０９）。

次に、計数判定部１７０４は、平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１（図２０において、ｂ１＝−１８）以上か否かを判定する（ステップＳＴ１９１０）。

平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１以上の場合には（ステップＳＴ１９１０：ＹＥＳ）、ステップＳＴ１９１１の処理を行わずに、ステップＳＴ１９１２に処理を進める。

一方、平均ＩＱ成分１の平均Ｑ成分の絶対値が閾値ｂ１未満の場合には（ステップＳＴ１９１０：ＮＯ）、計数判定部１７０４は、平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２（図２０において、ｂ２＝−２５）以上か否かを判定する（ステップＳＴ１９１１）。

平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２未満の場合には（ステップＳＴ１９１１：ＮＯ）、ステップＳＴ１９１２およびステップＳＴ１９１３の処理を行わずに、処理を終了する。

一方、平均ＩＱ成分２の平均Ｑ成分の絶対値が閾値ｂ２以上の場合には（ステップＳＴ１９１１：ＹＥＳ）、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ１９１２）。

Ｑ成分の符号変化がない場合には（ステップＳＴ１９１２：ＮＯ）、ステップＳＴ１９１３の処理を行わずに、処理を終了する。

一方、Ｑ成分の符号変化がある場合には（ステップＳＴ１９１２：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算し（変移回数＋１）（ステップＳＴ１９１３）、処理を終了する。

なお、本実施の形態では、いずれか一方の閾値を満たした場合に、変移検出結果を変移回数の計数対象とする例を示したが、本発明はこれに限ったものではなく、全ての閾値を満たした場合に、変移検出結果を変移回数の計数対象としてもよい。

上記により、平均区間長によってフェージング変動以外の要因による影響の可能性が異なる場合でも、最大ドップラー周波数の検出精度を高めることができる。

例えば、平均区間長が２シンボルの１種類のみを有し、フェージング変動が高速であり、かつ受信信号２シンボルの同相成分の符合が各々反対で大きさがほぼ等しい場合、この２シンボルを平均区間長で平均した平均Ｉ成分および平均Ｑ成分のＩ成分およびＱ成分は小さくなる。この結果、フェージング変動が高速であるにも関わらず、最大ドップラー周波数が小さく算出され、最大ドップラー周波数の検出精度が低下する場合がある。

一方、平均区間長として２シンボル（現在（時刻ｉ）の受信信号と直前（時刻ｉ−１）の受信信号）と１シンボル（現在（時刻ｉ）の受信信号）の２種類を有し、フェージング変動が高速であり、かつ受信信号２シンボルの同相成分の符合がそれぞれ反対で大きさがほぼ等しい場合、２シンボルで平均した平均Ｉ成分および平均Ｑ成分のＩ成分およびＱ成分は小さく算出される。しかし、１シンボルで平均した平均Ｉ成分および平均Ｑ成分は、２シンボルで平均した平均Ｉ成分および平均Ｑ成分に比べて大きく算出される。この結果、変移回数を正しく計数することができ、最大ドップラー周波数の検出精度を高めることができる。

このように、本実施の形態によれば、上記の実施の形態１の効果に加えて、フェージング変動が高速である場合において、最大ドップラー周波数の推定精度を向上させることができる。

（実施の形態６）
図２１は、本発明の実施の形態６に係る最大ドップラー周波数推定装置２１００の構成を示すブロック図である。最大ドップラー周波数推定装置２１００は、例えば、無線受信装置に搭載される。

図２１に示す最大ドップラー周波数推定装置２１００は、図２に示す実施の形態１に係る最大ドップラー周波数推定装置１００に対して、ＣＰＩＣＨシンボルの受信パワー測定部１０４を除き、ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１を追加し、計数判定部１０５の代わりに計数判定部２１０２を有する。なお、図２１において、図２と同一構成である部分には同一の符号を付してその説明を省略する。

最大ドップラー周波数推定装置２１００は、アンテナ１０１と、無線受信部１０２と、ＣＰＩＣＨシンボルのデータ抽出部１０３と、変移検出部１０６と、変移回数計数部１０７と、ｆｄ推定値検出部１０８と、ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１と、計数判定部２１０２とから主に構成される。以下に、本実施の形態における上記の実施の形態１とは異なる各構成について詳細に説明する。

ＣＰＩＣＨシンボルのデータ抽出部１０３は、無線受信部１０２から入力したＩ成分およびＱ成分に対し、逆拡散を施しＣＰＩＣＨシンボルのＩ成分およびＱ成分を抽出する。また、ＣＰＩＣＨシンボルのデータ抽出部１０３は、抽出したＣＰＩＣＨシンボルのＩ成分およびＱ成分をＣＰＩＣＨシンボルのＳＩＲ測定部２１０１および変移検出部１０６へ出力する。

ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１は、ＣＰＩＣＨシンボルのデータ抽出部１０３から入力したＣＰＩＣＨシンボルのＩ成分およびＱ成分からＩＳＣＰおよびＲＳＣＰを算出し、算出したＩＳＣＰおよびＲＳＣＰから受信ＳＩＲを算出する。そして、ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１は、算出した受信ＳＩＲを計数判定部２１０２へ出力する。

計数判定部２１０２は、ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１から入力した受信ＳＩＲに基づいて、受信信号に含まれるフェージング変動以外の要因による影響の可能性を推定する。具体的には、計数判定部２１０２は、受信ＳＩＲが閾値以上の場合にはフェージング変動以外の要因による影響の可能性が小さいと判定し、受信ＳＩＲが閾値未満の場合にはフェージング変動以外の要因による影響の可能性が大きいと判定する。また、計数判定部２１０２は、フェージング変動以外の要因による影響の可能性が小さいと判定した場合には、変移を計数するように変移回数計数部１０７を制御し、フェージング変動以外の要因による影響の可能性が大きいと判定した場合には、変移を計数しないように変移回数計数部１０７を制御する。

変移回数計数部１０７は、予め設定された観測時間内において、計数判定部２１０２の制御に従って、変移検出部１０６から入力した変移の検出結果をカウントし、変移したＣＰＩＣＨシンボルの個数を計数する。また、変移回数計数部１０７は、計数結果をｆｄ推定値検出部１０８へ出力する。以上で、最大ドップラー周波数推定装置２１００の構成の説明を終える。

次に、最大ドップラー周波数推定装置２１００の動作について、図２２を用いて説明する。図２２は、最大ドップラー周波数推定装置２１００の動作の一部を示すフロー図である。

最初に、計数判定部２１０２は、ＣＰＩＣＨシンボルのＳＩＲ測定部２１０１により算出した受信ＳＩＲが閾値Ｐ以上であるか否かを判定する（ステップＳＴ２２０１）。

受信ＳＩＲが閾値Ｐ未満の場合には（ステップＳＴ２２０１：ＮＯ）、ステップＳＴ２２０２〜ステップＳＴ２２０６の処理を行わずに、処理を終了する。

一方、受信ＳＩＲが閾値Ｐ以上の場合には（ステップＳＴ２２０１：ＹＥＳ）、計数判定部２１０２は、誤検出の可能性が小さいと推定し、変移回数計数時間に「１」を加算する（変移回数計数時間＋１）（ステップＳＴ２２０２）。

次に、変移回数計数部１０７は、検出された変移を計数する。この際、まず、変移検出部１０６は、Ｉ成分の符号変化があるか否かを判定する（ステップＳＴ２２０３）。

Ｉ成分の符号変化がない場合には（ステップＳＴ２２０３：ＮＯ）、ステップＳＴ２２０４の処理を行わずに、ステップＳＴ２２０５に処理を進める。

Ｉ成分の符号変化がある場合には（ステップＳＴ２２０３：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算する（変移回数＋１）（ステップＳＴ２２０４）。

次に、変移検出部１０６は、Ｑ成分の符号変化があるか否かを判定する（ステップＳＴ２２０５）。

Ｑ成分の符号変化がない場合には（ステップＳＴ２２０５：ＮＯ）、ステップＳＴ２２０６の処理を行わずに、処理を終了する。

一方、Ｑ成分の符号変化がある場合には（ステップＳＴ２２０５：ＹＥＳ）、変移回数計数部１０７は、変移回数に「１」を加算し（変移回数＋１）（ステップＳＴ２２０６）、処理を終了する。

このように、本実施の形態によれば、ＣＰＩＣＨシンボルの受信信号がフェージング変動以外の要因による影響を受ける場合に、フェージング変動の影響のみを受けた最大ドップラー周波数の推定精度を向上させることができる。

なお、本実施の形態において、現在の受信信号から受信ＳＩＲを算出したが、本発明はこれに限らず、直前を含む過去の受信信号を用いて算出した受信ＳＩＲを用いてもよい。

また、上記の実施の形態１〜実施の形態６において、Ｗ−ＣＤＭＡシステムを対象として説明しているが、本発明はこれに限らず、受信信号の変移によりフェージング変動の最大ドップラー周波数を推定する、Ｗ−ＣＤＭＡシステム以外の無線通信システムに適用可能である。また、上記の実施の形態１〜実施の形態６において、ＣＰＩＣＨシンボルを用いたが、本発明はこれに限らず、パイロットチャネルであれば、任意のチャネルのシンボルを用いることができる。

本発明にかかる最大ドップラー周波数推定装置および最大ドップラー周波数推定方法は、特にフェージング変動の速さに関連する最大ドップラー周波数を精度良く検出するのに好適である。

１００ 最大ドップラー周波数推定装置
１０１ アンテナ
１０２ 無線受信部
１０３ ＣＰＩＣＨシンボルのデータ抽出部
１０４ ＣＰＩＣＨシンボルの受信パワー測定部
１０５ 計数判定部
１０６ 変移検出部
１０７ 変移回数計数部
１０８ ｆｄ推定値検出部

Claims (8)

1. 受信信号の時間的な位相変移を検出する位相変移検出手段と、
   前記位相変移に対するフェージング変動以外の要因による影響の大きさを推定する推定手段と、
   前記推定の結果に基づいて前記位相変移の回数をカウントする計数手段と、
   カウントした前記回数より最大ドップラー周波数の推定値を検出する推定値検出手段と、
   を具備する最大ドップラー周波数推定装置。
2. 前記推定手段は、受信信号の受信パワーと閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記受信パワーが前記閾値以上の場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
3. 受信信号を周波数変換してＩ成分およびＱ成分を生成する生成手段をさらに具備し、
   前記推定手段は、前記Ｉ成分の値と第１閾値との比較および前記Ｑ成分の値と第２閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記Ｉ成分の値が前記第１閾値以上の場合または前記Ｑ成分の値が前記第２閾値以上の場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
4. 受信信号を周波数変換してＩ成分およびＱ成分を生成する生成手段と、
   前記生成手段により生成した現在のＩ成分およびＱ成分と直前のＩ成分およびＱ成分との演算により、Ｉ成分の変移量およびＱ成分の変移量を求める変移量取得手段とをさらに具備し、
   前記推定手段は、前記変移量取得手段により求めた前記Ｉ成分の変移量と第１閾値との比較、および前記変移量取得手段により求めた前記Ｑ成分の変移量と第２閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記Ｉ成分の変移量が前記第１閾値以上の場合または前記Ｑ成分の変移量が前記第２閾値以上の場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
5. 受信信号を周波数変換してＩ成分およびＱ成分を生成する生成手段をさらに具備し、
   前記推定手段は、所定の平均化時間における前記Ｉ成分の平均値と第１閾値との比較および前記Ｑ成分の平均値と第２閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記Ｉ成分の平均値が前記第１閾値以上の場合または前記Ｑ成分の平均値が前記第２閾値以上の場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
6. 受信信号を周波数変換してＩ成分およびＱ成分を生成する生成手段をさらに具備し、
   前記推定手段は、異なる平均化時間における前記Ｉ成分の各平均値と第１閾値との比較、および前記異なる平均化時間における前記Ｑ成分の各平均値と第２閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記Ｉ成分の各平均値および前記Ｑ成分の各平均値のうち、前記第１閾値以上または前記第２閾値以上の平均値が少なくとも１つある場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
7. 前記推定手段は、受信信号の受信ＳＩＲと閾値との比較により前記フェージング変動以外の要因による影響の大きさを推定し、
   前記計数手段は、前記受信ＳＩＲが前記閾値以上の場合に前記回数をカウントする請求項１記載の最大ドップラー周波数推定装置。
8. 最大ドップラー周波数を推定する最大ドップラー周波数推定装置における最大ドップラー周波数推定方法であって、
   受信信号間の時間的な位相変移を検出するステップと、
   前記位相変移に対するフェージング変動以外の要因による影響の大きさを推定するステップと、
   前記推定の結果に基づいて前記位相変移の回数をカウントするステップと、
   カウントした前記回数より最大ドップラー周波数の推定値を検出するステップと、
   を具備する最大ドップラー周波数推定方法。

Images