チャンネルの周波数領域の相関性の算出方法であって、サンプルチャンネルのパス遅延分布を解析し、前記サンプルチャンネルのパスエネルギー分布を解析し、前記サンプルチャンネルの前記パス遅延分布と前記パスエネルギー分布に基づいてチャンネルモデルを構築し、構築された前記チャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する、ことを含む方法。遅延スプレッドに基づいて複数の解析区間を選定することをさらに含み、前記パス遅延分布の解析は、前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パス遅延分布を解析し、前記パスエネルギー分布の解析は、前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パスエネルギー分布を解析する、請求項1に記載の方法。前記複数の解析区間のそれぞれにおいて複数のパスの遅延は、均一分布または等差逓増分布を呈し、前記複数の解析区間のそれぞれにおいて前記複数のパスのエネルギー分布は、線形逓減分布、指数逓減分布或は等エネルギー均一分布を呈する、請求項2に記載の方法。チャンネルの周波数領域の相関性の算出装置であって、サンプルチャンネルのパス遅延分布を解析するパス遅延解析手段と、前記サンプルチャンネルのパスエネルギー分布を解析するパスエネルギー解析手段と、前記サンプルチャンネルの前記パス遅延分布と前記パスエネルギー分布に基づいてチャンネルモデルを構築するチャンネルモデル構築手段と、構築された前記チャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する周波数領域相関性算出手段と、を備える算出装置。遅延スプレッドに基づいて複数の解析区間を選定する解析区間選択手段をさらに備え、前記パス遅延解析手段は前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パス遅延分布を解析し、前記パスエネルギー解析手段は前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パスエネルギー分布を解析する、請求項4に記載の算出装置。前記チャンネルモデル構築手段は、前記複数の解析区間のそれぞれにおいて複数のパスの遅延が均一分布又は等差逓増分布を呈し、且つ前記複数の解析区間のそれぞれにおいて前記複数のパスのエネルギー分布が線形逓減分布、指数逓減分布または等エネルギー均一分布を呈する原則に従って、前記チャンネルモデルを構築する請求項5に記載の算出装置。チャンネル推定方法であって、パイロットサブキャリアの最小二乗チャンネル推定を算出し、請求項1~3の何れか一項に記載のチャンネルの周波数領域の相関性の算出方法により得られた、予め記憶されたチャンネルの周波数領域の相関性を読み取り、受信側の送信側に対する移動速度を推定し、前記受信側の前記送信側に対する前記移動速度に基づいてチャンネルの時間領域の相関性を算出し、算出された前記チャンネルの時間領域の相関性と読み取られた前記チャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出し、前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定する、ことを含む方法。前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用して前記最小平均二乗誤差フィルタ係数を算出することは、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用して前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネル推定を取得することは、前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用して前記パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行い、かつ前記最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいて前記データサブキャリアのチャンネルに対して推定を行って前記データサブキャリアのチャンネルを推定することを含む、請求項7に記載の方法。パイロットサブキャリアの最小二乗チャンネル推定を算出するパイロットチャンネル推定手段と、請求項4~6の何れか一項に記載のチャンネルの周波数領域の相関性の算出装置により得られた、予め記憶されたチャンネルの周波数領域の相関性を読み取る周波数領域相関性読取手段と、受信側の送信側に対する移動速度を推定する移動速度推定手段と、前記受信側の前記送信側に対する前記移動速度に基づいてチャンネルの時間領域の相関性を算出する時間領域相関性算出手段と、算出された前記チャンネルの時間領域の相関性と読み取られた前記チャンネルの周波数領域の相関性を利用して最小平均二乗誤差ファイル係数を算出するフィルタ係数算出手段と、算出された前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定するデータチャンネル推定手段と、を備えるチャンネル推定装置。前記フィルタ係数算出手段は、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用してパイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記データチャンネル推定手段は、前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用して前記パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行うパイロット最小平均二乗誤差フィルタ手段と、前記最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいて前記データサブキャリアのチャンネルに対して推定を行って前記データサブキャリアのチャンネルを推定するデータチャンネル推定器とを備える請求項9に記載のチャンネル推定装置。これに鑑みて、本明細書は、MMSEチャンネル推定アルゴリズムを採用する際のシステムコストが大きく削減され、かつシステム遅延が減少される、非常に簡単な、チャンネル相関性を算出する方法及びその装置、ならびに相応するチャンネル推定方法及びその装置を提供する。上記の目的を実現するため、一つの実施形態によれば、チャンネルの周波数領域の相関性の算出方法が提供される。この方法は、サンプルチャンネルのパス遅延分布を解析し、サンプルチャンネルのパスエネルギー分布を解析し、サンプルチャンネルのパス遅延分布とパスエネルギー分布に基づいてチャンネルモデルを構築し、構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する。このチャンネルの周波数領域の相関性の算出方法は、遅延スプレッドに基づいて複数の解析区間を選定することをさらに含むことが好ましい。この場合に、解析区間毎にサンプルチャンネルのパス遅延分布を解析し、且つ、解析区間毎にサンプルチャンネルのパスエネルギー分布を解析することが好ましい。好ましくは、各解析区間において複数のパス遅延は、均一分布又は等差逓増分布、あるいは他の分布を呈する。また好ましくは、各解析空間において各パスのエネルギー分布は線形逓減分布、指数逓減分布、等エネルギー均一分布、あるいは他の分布を呈する。チャンネルの周波数領域の相関性の算出は、異なる最大遅延スプレッドに対して、下記の式により、構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性をそれぞれ算出することが好ましい。ここで、θΔfkkmはチャンネルの周波数領域の相関性関数であり、absは絶対値の計算操作であり、kkmはチャンネルの周波数領域の相関値を算出する両周波数点間のサブキャリア間隔数であり、iはパスのシーケンス番号であり、σi2はi番目のパスのパスエネルギーであり、τiはi番目のパスのパス遅延である、ここで、サンプルチャンネルは、PAチャンネル、EPAチャンネル、VAチャンネル、EVAチャンネル、PBチャンネル、TUチャンネルとETUチャンネル等のチャンネルのうちの一個又は複数個を含む。他の実施形態によれば、チャンネルの周波数領域の相関性の算出装置が提供される。この装置は、サンプルチャンネルのパス遅延分布を解析するパス遅延解析手段と、サンプルチャンネルのパスエネルギー分布を解析するパスエネルギー解析手段と、サンプルチャンネルのパス遅延分布とパスエネルギー分布に基づいてチャンネルモデルを構築するチャンネルモデル構築手段と、構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する周波数相関性算出手段とを備える。このチャンネルの周波数領域の相関性の算出装置は、遅延スプレッドに基づいて複数の解析区間を選定するように配置された解析区間選択手段をさらに備えることが好ましい。この場合に、パス遅延解析手段は解析区間毎にサンプルチャンネルのパス遅延分布を解析し、パスエネルギー解析手段は解析区間毎にサンプルチャンネルのパスエネルギー分布を解析することが好ましい。好ましくは、チャンネルモデル構築手段は、各解析区間において複数のパスの遅延が均一分布、等差逓増分布、或は他の分布を呈し、各解析空間において各パスのエネルギー分布が線形逓減分布、指数逓減分布、等エネルギー均一分布、或は他の分布を呈する原則に従って、チャンネルモデルを構築する。さらに他の実施形態によれば、チャンネル推定方法が提供される。この推定方法は、パイロットサブキャリアの最小二乗チャンネル推定を算出し、上記の実施形態またはその変形例による各チャンネルの周波数領域の相関性の算出方法により得られ、予め記憶されたチャンネルの周波数領域の相関性を読み取り、受信側の送信側に対する移動速度を推定し、受信側の送信側に対する移動速度に基づいてチャンネルの時間領域の相関性を算出し、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出し、算出された最小平均二乗誤差フィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定する。このチャンネル推定方法は、チャンネルの最大遅延スプレッドを推定し、推定された最大遅延スプレッドに基づいて、それに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取ることをさらに含むことが好ましい。このチャンネル推定方法において、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出することは、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してパイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出することを含むことが好ましい。この場合に、算出された最小平均二乗誤差フィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定することは、算出されたパイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用してパイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行うこと、及び、最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいてデータサブキャリアのチャンネルに対して推定を行って、データサブキャリアのチャンネルを推定することを含むことが好ましい。好ましくは、最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに対して線形補間、2次補間、最小平均二乗誤差フィルタリング、ローパスフィルタリングあるいはスプライン補間を行うことにより、データサブキャリアのチャンネルを推定することが好ましい。このチャンネル推定方法によれば、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出することは、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してデータサブキャリアの最小平均二乗誤差フィルタ係数を算出することを含むことが好ましい。この場合に、算出された最小平均二乗誤差フィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネル推定を取得することは、算出されたデータサブキャリアの最小平均二乗誤差フィルタ係数を利用してパイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行って、データサブキャリアのチャンネルを推定することを含むことが好ましい。さらに他の実施形態によれば、チャンネルの推定装置が提供される。この推定装置は、パイロットサブキャリアの最小二乗チャンネル推定を算出するパイロットチャンネル推定手段と、上記の実施形態またはその変形例による各チャンネルの周波数領域の相関性の算出装置により得られ、かつ予め記憶されたチャンネルの周波数領域の相関性を読み取る周波数領域相関性読取手段と、受信側の送信側に対する移動速度を推定する移動速度推定手段と、受信側の送信側に対する移動速度に基づいてチャンネルの時間領域の相関性を算出する時間領域相関性算出手段と、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出するフィルタ係数算出手段、及び算出された最小平均二乗誤差フィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定するデータチャンネル推定手段とを備える。このチャンネル推定装置は、さらに、チャンネルの最大遅延スプレッドを推定する遅延スプレッド推定手段、及び推定された最大遅延スプレッドに基づいて、それに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取る周波数領域相関性読取手段を備えることが好ましい。このチャンネル推定装置において、フィルタ係数算出手段は、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してパイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出することが好ましい。この場合に、データチャンネル推定手段は、算出されたパイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用してパイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行うパイロット最小平均二乗誤差フィルタ手段、及び最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいてデータサブキャリアのチャンネルに対して推定を行って、データサブキャリアのチャンネルを推定するデータチャンネル推定器を備えることが好ましい。好ましくは、データチャンネル推定器は、最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに対して線形補間、2次補間、最小平均二乗誤差フィルタリング、ローパスフィルタリングあるいはスプライン補間を行うことにより、データサブキャリアのチャンネルを推定する。このチャンネル推定装置において、フィルタ係数算出手段は算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してデータサブキャリアの最小平均二乗誤差フィルタ係数を算出し、データチャンネル推定手段は算出されたデータサブキャリアの最小平均二乗誤差フィルタ係数を利用して、パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行ってデータサブキャリアのチャンネルを推定することが好ましい。さらに他の実施形態によれば、機器読取可能のプログラムコードを含む記憶媒体が提供される。情報処理装置に上記のプログラムコードが読み込まれ、情報処理装置がそのプログラムコードを実行すると、そのプログラムコードは情報処理装置に上記のチャンネルの周波数領域の相関性の算出方法またはチャンネル推定方法の何れかを実行させる。さらに他の実施形態によれば、機器実行可能のコマンドを含むコンピュータプログラムが提供される。情報処理装置に上記のコマンドを実行すると、そのコマンドは、情報処理装置に上記のチャンネルの周波数領域の相関性の算出方法またはチャンネル推定方法の何れかを実行させる。図面とともに実施形態の説明を参照することにより、本発明の以上及びその他の目的、特徴と利点をより容易に理解することができる。図面における構成要素は本発明の原理を説明するために示されるものに過ぎない。図面において、同一又は類似の技術的特徴又は構成要素は、同一又は類似の図面符号で示される。以下に、各実施形態を図面に基づいて記述する。明瞭、簡便のため、明細書には実際の実施形態の全ての特徴が記述されているわけではない。しかし、以下のことが理解されるべきであり、例えば、このような実際の実施形態を開発する過程において、開発者の具体的な目標を実現するためには、実施形態に特定される多くのことを決定しなければならない。例えば、システムと業務に関連する様々な制限条件を満たすことであり、これらの制限条件は相違する実施形態によって変化する可能性がある。なお、開発の仕事は非常に複雑であり、時間を費やす可能性があるのだが、ここに開示している内容を利用する当業者には、このような開発の仕事は型どおりの任務だけである。必要以上に詳述することによって本発明が不明瞭になるのを防ぐために、図面には、実施形態に密接な関係がある装置構成及びまたは処理ステップだけが示されており、実施形態にあまり関係のないその他の詳細は省略されている。また、一つの図面又は一種類の実施形態において説明された要素と特徴は、1以上の他の図面あるいは実施形態において示された要素及び特徴と組み合わせることができる。以下、先ず図1~図3を参照して、一つの実施形態によるチャンネルの周波数領域の相関性の算出方法及びその装置の基本的動作原理を説明する。図1に示すように、一実施形態によるチャンネルの周波数領域の相関性の算出方法は、解析サンプルチャンネル選択ステップS110、解析区間選択ステップS120、パス遅延分布解析ステップS130、パスエネルギー分布解析ステップS140、チャンネルモデル構築ステップS150及びチャンネル周波数相関性算出ステップS160等を含む。先ず、解析サンプルチャンネル選択ステップS110において解析用のサンプルチャンネルを特定する。システムへの適用において発生可能な全てのシーンを考慮したうえ、一実施例では、PedestrianAPAチャンネル、ExtendedPedestrianAEPAチャンネル、VehicularAVAチャンネル、ExtendedVehicularAEVAチャンネル、PedestrianBPBチャンネル、TypicalUrbanTUチャンネル及びExtendedTypicalUrbanETUチャンネルのうちの少なくとも一つを解析用のサンプルチャンネルとして選択する。次に、解析区間選択ステップS120において、時刻0からチャンネルの遅延スプレッドまでの区間を観測区間全体とし、必要に応じて解析区間を特定する。例えば、一実施例では、三つの解析区間0T4、0T2と0Tを選択することができる。ここで、Tは構築待ちチャンネルモデルの遅延スプレッドである。勿論、その他の数の解析区間等、例えば四つの解析区間0T8、0T4、0T2と0Tを選択してもよい。その後、パス遅延分布解析ステップS130において、解析区間毎に、全てのサンプルチャンネルのパス遅延分布の統計値を算出する。一実施例として、上記のような三つの解析区間が選択された場合に対して統計的に得られた各サンプルチャンネルのパス遅延分布を、表1に示す。次に、パスエネルギー分布解析ステップS140において、解析区間毎に、全てのサンプルチャンネルのパスエネルギー分布の統計値を算出する。一実施例として、同様に、統計的に得られた各サンプルチャンネルのパスエネルギー分布を表2に示す。その後、チャンネルモデル構築ステップS150において、パス遅延分布解析ステップS130で得られたサンプルチャンネルのパス遅延分布とパスエネルギー分布解析ステップS140で得られたパスエネルギー分布に基づいてチャンネルモデルを構築する。ここで、表1及び表2に示された各サンプルチャンネルのパス遅延分布及びパスエネルギー分布を例としてわかるように、第1の解析区間0T4のマルチパス数の平均値は3429であり、その平均マルチパスエネルギーは総エネルギーの873834である。また、第2の解析区間0T2のマルチパス数の平均値は5286であり、その平均マルチパスエネルギーは総エネルギーの97347である。さらに、第3の解析区間0Tのマルチパス数の平均値は6714であり、その平均マルチパスエネルギーは総エネルギーの100000である。以上の平均データから、ここで構築されたチャンネルモデルのパス遅延分布及びパスエネルギー分布は、第1の解析区間0T4のマルチパス数が4、そのマルチパスエネルギーが総エネルギーの87となり、第2のの解析区間0T2のマルチパス数が5、そのマルチパスエネルギーが総エネルギーの97となり、第3のの解析区間0Tのマルチパス数が7、そのマルチパスエネルギーが総エネルギーの100となるように選択される。図2は構築されたチャンネルモデルの模式図である。図2に示すように、選択されたパラメータによれば、区間0T4において4個のパス、区間T4T2において1個のパス、区間T2Tにおいて2個のパスがある。一実施例によれば、0T4における4個のパスの遅延は均一分布を呈する。また区間T4T2における1個のパスは均一分布を呈する。また区間T2Tにおいて、2個目のパスはチャンネルの最大遅延スプレッドに位置し、区間T2Tにおいて1個目のパスはその区間の中央に位置すると想定される。そのほか、各パスの遅延が等差逓増分布等を呈するとしてもよい。また、各解析区間において、マルチパスの総エネルギーは表2に示された選択パラメータに従い、各パスのエネルギー分布が線形逓減を呈する。そのほか、各パスのエネルギー分布が指数逓減分布を呈するとしてもよく、また各パスのエネルギーが等エネルギーの均一分布を呈するとしてもよく、あるいは各パスのエネルギーがその他の可能な分布などを呈するとしてもよい。最後に、チャンネル周波数領域相関性算出ステップS160において、チャンネルモデル構築ステップS150で構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する。OFDMシステムにおいては、一般的に、その経験した時間領域におけるマルチパスチャンネルの最大遅延スプレッドがシステムに採用された巡回プレフィックスの長さを超えないとみなしている。そのため、巡回プレフィックスの長さが参考できるパラメータとなる。異なる適用環境において、システムの巡回プレフィックスの長さに対して異なる比例値、例えば巡回プレフィックス長さの18、28、38、48、58、68、78、88がチャンネルモデルの最大遅延スプレッドとして選択される。そして下記の15式に従って異なる遅延スプレッドのチャンネルの周波数領域の相関値を算出する。ここで、θΔfkkmはチャンネルの周波数領域の相関性関数であり、absは絶対値の計算操作である。またkkmはチャンネルの周波数領域の相関値を算出する両周波数点間のサブキャリア間隔数であり、iはパスのシーケンス番号であり、σi2はi番目のパスのパスエネルギーであり、τiはi番目のパスのパス遅延である。ここで、端末固有参照信号UESpecificReferenceSignalを有するLTE時分割多重TimeDivisionDuplex、TDDシステムの具体的な状況を考慮し、仮に推定待ちリソースが1リソースブロックResourceBlock、RBに制限されると想定する。この場合、整数変数kkmの値の取る区間を010に設けることができる。下記の表3は、チャンネルモデルの最大遅延スプレッドが48巡回プレフィックスの長さとして選択され、且つkkmの値の取る区間が010である場合の一つのチャンネルの周波数領域の相関性の具体的な算出結果を示している。ここで、チャンネル周波数領域相関性算出ステップS160の算出結果は、チャンネル推定装置が動作する時の使用に備えてテーブルの形式で保存できる。また、上記の解析サンプルチャンネル選択ステップS110と解析区間選択ステップS120は必須なものではない。決められたサンプルチャンネルに対して必要に応じてパス遅延分布解析ステップS130、パスエネルギー分布解析ステップS140、チャンネルモデル構築ステップS150及びチャンネル周波数領域相関性算出ステップS160等のステップを直接実行してもよい。以下に図3を参照して一つの実施形態によるチャンネルの周波数領域の相関性の算出装置の動作原理を説明する。図3に示すように、この実施形態によるチャンネルの周波数領域の相関性の算出装置は、パス遅延解析手段320と、パスエネルギー解析手段330と、チャンネルモデル構築手段340と、周波数領域相関性算出手段350とを備える。パス遅延解析手段320はサンプルチャンネルのパス遅延分布を解析する。またパスエネルギー解析手段330はサンプルチャンネルのパスエネルギー分布を解析する。チャンネルモデル構築手段340はサンプルチャンネルのパス遅延分布及びパスエネルギー分布に基づいてチャンネルモデルを構築する。そして周波数領域相関性算出手段350は構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する。また、この実施形態によるチャンネルの周波数領域の相関性の算出装置は、遅延スプレッドに基づいて複数の解析区間を選択するように配置される解析区間選択手段310をさらに備えてもよい。この場合に、パス遅延解析手段320は解析区間毎にサンプルチャンネルのパス遅延分布を解析し、パスエネルギー解析手段330は解析区間毎にサンプルチャンネルのパスエネルギー分布を解析する。ここで、サンプルチャンネルは、PAチャンネル、EPAチャンネル、VAチャンネル、EVAチャンネル、PBチャンネル、TUチャンネルとETUチャンネル等のチャンネルのうちの一個又は複数個を含むことができる。チャンネルの周波数領域の相関性の算出装置に含まれる解析区間選択手段310、パス遅延解析手段320、パスエネルギー解析手段330、チャンネルモデル構築手段340、及び周波数領域相関性算出手段350における処理のそれぞれは、上記したチャンネルの周波数領域の相関性の算出方法に含まれる解析区間選択ステップS120、パス遅延分布解析ステップS130、パスエネルギー分布解析ステップS140、チャンネルモデル構築ステップS150、及びチャンネル周波数領域相関性算出ステップS160における処理と同様である。そのため、簡潔のためこれら手段の詳細説明を省略する。以上のように、図1~図3を参照して一つの実施形態によるチャンネルの周波数領域の相関性の算出方法及びチャンネルの周波数領域の相関性の算出装置の動作原理を説明した。以下、図4~図7を参照して上記のチャンネルの周波数領域の相関性の算出方法及びチャンネルの周波数領域の相関性の算出装置を利用したチャンネル推定方法とチャンネル推定装置の動作原理を具体的に説明する。図4は他の実施形態によるチャンネル推定方法における処理のフローチャートを示す。図4に示すように、この実施形態によるチャンネル推定方法は、パイロット最小二乗LeastSquare、LSチャンネル推定ステップS410、最大遅延スプレッド推定ステップS420、チャンネル周波数領域相関性読取ステップS430、移動速度推定ステップS440、チャンネル時間領域相関性算出ステップS450、MMSEフィルタ係数算出ステップS460、パイロットMMSEフィルタリングステップS470、及びデータサブキャリアチャンネル推定ステップS480等を含む。先ず、パイロットLSチャンネル推定ステップS410において、システムのパイロットサブキャリアにおける周波数領域のチャンネル応答に対する推定を行う。仮にパイロットの位置するサブキャリアがklであるとすると、そのパイロットサブキャリアのLSチャンネル推定結果は以下のように示される。ここで、YklはパイロットXklの受信値である。最大遅延スプレッド推定ステップS420において、時間領域におけるマルチパスチャンネルの最大遅延スプレッドを推定する。そして、チャンネル周波数領域相関性読取ステップS430において、最大遅延スプレッド推定ステップS420で推定されたマルチパスチャンネルの最大遅延スプレッドをシーリングし、即ち、それより大きく、かつそれに最も近い値を所定の値として選択する。その後、シーリングして選択された所定の値に基づいて相応するチャンネルの周波数領域の相関値の予め記憶された表を読み出す。ここで、LTETDDシステムにおいて、端末固有参照信号はマルチパスチャンネルの最大遅延スプレッドの推定に用いることができないので、最大遅延スプレッド推定ステップS420は必須なものではない。この場合に、一つのチャンネルの周波数領域の相関性の表しか予め記憶できないので、最大遅延スプレッド推定ステップS420を省略し、チャンネル周波数領域相関性読取ステップS430で予め記憶されたチャンネルの周波数領域の相関性の予め記憶された表を直接に読み出せばよい。しかしながら、その他のシステム、例えば共通参照信号CommonReferenceSignalを利用してチャンネル推定を行うLTEシステムの場合、共通参照信号の推定結果により周波数領域のチャンネルを時間領域へ変換すると共にノイズを除去し、最後に時間領域のチャンネルの推定結果に基づいてチャンネルの最大遅延スプレッド推定を粗略に推定することができる。一般的に、チャンネルの最大遅延スプレッド推定値の誤差が大きくても、良好なチャンネル推定効果を取得可能である。そして、チャンネル周波数領域相関性読取ステップS430において、推定された最大遅延スプレッド推定に基づいてそれに相応するチャンネルの周波数領域の相関性の予め記憶された表を選定して読み出す。また、移動速度推定ステップS440において受信側の送信側に対する移動速度を推定する。そしてチャンネル時間領域相関性算出ステップS450において、移動速度推定ステップS440で推定された受信側の送信側に対する移動速度の推定結果を利用して、上記の13式と14式に従って、チャンネルの時間領域の相関性を算出する。例えば、移動速度の推定について、LTETDDシステムでは、共通参照信号を採用することができ、共通参照信号の経たチャンネルの時間領域の相関性を推定することによって、受信側の移動速度を推定する。なお、先に移動速度推定ステップS440とチャンネル時間領域相関性算出ステップS450を実行し、それから、最大遅延スプレッド推定ステップS420とチャンネル周波数領域相関性読取ステップS430を実行してもよい。チャンネルの周波数領域の相関性と時間領域の相関性を取得できればよい。チャンネル周波数領域相関性読取ステップS430とチャンネル時間領域相関性算出ステップS450でチャンネルの周波数領域の相関性と時間領域の相関性をそれぞれ取得する。その後に、MMSEフィルタ係数算出ステップS460において、上記の10式と11式に従って相互相関ベクトルθ及び自己相関行列Фをそれぞれ算出すると共に、5式に従ってMMSEフィルタ係数を算出する。その具体的な実施例において、演算時間を節約するため、チャンネル時間領域相関性算出ステップS450で算出されたチャンネル時間領域相関性及びチャンネル周波数領域相関性読取ステップS430で読み取られたチャンネル周波数領域相関性のみを利用してパイロットサブキャリアのMMSEフィルタ係数を算出する。次に、パイロットMMSEフィルタリングステップS470において、MMSEフィルタ係数算出ステップS460で算出されたパイロットサブキャリアのMMSEフィルタ係数を利用して、パイロットLSチャンネル推定ステップS410で出力されたパイロットサブキャリアにおける周波数領域のチャンネル応答をフィルタリングする。これにより、これらパイロットサブキャリアにおける周波数領域のチャンネル応答の信頼性が向上する。そして、データサブキャリアチャンネル推定ステップS480において、パイロットMMSEフィルタリングステップS470で出力された、MMSEフィルタリングされたパイロットサブキャリアにおける周波数領域のチャンネル応答の推定値を利用して、推定待ちリソース中のデータサブキャリアにおける周波数領域のチャンネル応答に対して推定を行う。推定の過程において、具体的な推定方法は限定されない。一般的に時間領域と周波数領域のそれぞれに対して線形補間を行うような最も簡単なチャンネル推定方法を採用するだけで非常に良い結果を得ることができる。勿論、2次補間、MMSEフィルタリング、ローパスフィルタリング及びスプライン補間等の方法を採用してもよい。例えば、データサブキャリアチャンネル推定ステップS480において線形補間のチャンネル推定方法を採用する場合に、先ずMMSEフィルタリングされた、パイロットサブキャリアにおける周波数領域のチャンネル応答の推定値に対して解析を行うことができる。チャンネルの時間領域に沿った変化がチャンネルの周波数領域に沿った変化よりも早ければ、先に周波数領域の線形補間を行い、それから、その周波数領域の線形補間結果を利用して時間領域の線形補間を行う。さもなければ、先に、時間領域の線形補間を行い、それから、その時間領域の線形補間結果を利用して周波数領域の線形補間を行う。最後に、データサブキャリアチャンネル推定ステップS480のチャンネル推定結果が出力され、受信機における後続の等化、復調等の処理に用いられる。なお、チャンネル周波数領域相関性読取ステップS430で読み取られた、予め記憶されたチャンネルの周波数領域の相関性は、上記の各実施形態によるチャンネルの周波数領域の相関性の算出方法又はチャンネルの周波数領域の相関性の算出装置により算出されたものである。以下に、図5を参照して一つの実施形態によるチャンネル推定装置を説明する。図5に示すように、この実施形態によるチャンネル推定装置は、パイロットチャンネル推定手段501と、周波数領域相関性読取手段507と、移動速度推定手段504と、時間領域相関性算出手段505と、MMSEフィルタ係数算出手段506と、データチャンネル推定手段510とを備える。パイロットチャンネル推定手段501は、パイロットサブキャリアの最小二乗チャンネル推定を算出する。周波数領域相関性読取手段507は、予め記憶されたチャンネルの周波数領域の相関性を読み取る。なお、チャンネルの周波数領域の相関性は、例えば、チャンネル推定装置が有する半導体メモリ、磁気記録媒体あるいは光記録媒体などのメモリに予め記憶される。移動速度推定手段504は、受信側の送信側に対する移動速度を推定する。時間領域相関性算出手段505は、受信側の送信側に対する移動速度に基づいてチャンネルの時間領域の相関性を算出する。MMSEフィルタ係数算出手段506は、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してMMSEフィルタ係数を算出する。データチャンネル推定手段510は、算出されたMMSEフィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネル推定を取得する。この実施形態によるチャンネル推定装置において、周波数領域相関性読取手段507により読み取られた、上記の予め記憶されたチャンネルの周波数領域の相関性も、上記の各実施形態の何れかによるチャンネルの周波数領域の相関性の算出方法又はチャンネルの周波数領域の相関性の算出装置により算出されたものである。また、この実施形態によるチャンネル推定装置において、演算の複雑度を軽減するため、MMSEフィルタ係数算出手段506は時間領域相関性算出手段505により算出されたチャンネルの時間領域の相関性及び周波数領域相関性読取手段507により読み取られたチャンネルの周波数領域の相関性のみを利用してパイロットサブキャリアのMMSEフィルタ係数を算出する。この場合に、データチャンネル推定手段510は、パイロットMMSEフィルタ手段502とデータチャンネル推定器503を備える。パイロットMMSEフィルタ手段502は、MMSEフィルタ係数算出手段506により算出されたパイロットサブキャリアのMMSEフィルタ係数を利用してパイロットチャンネル推定手段501から出力されたパイロットサブキャリアの最小二乗チャンネル推定に対してMMSEフィルタリングを行う。データチャンネル推定器503は、パイロットMMSEフィルタ手段502から出力された、MMSEフィルタリングされたパイロットサブキャリアに基づいて、データサブキャリアのチャンネルに対して推定を行ってデータサブキャリアのチャンネル推定を取得する。データチャンネル推定器503は、推定を行う際に、時間領域と周波数領域のそれぞれに対して線形補間を行うような最も簡単なチャンネル推定方法を採用するだけで、非常に良い結果を得ることができる。勿論、データチャンネル推定器503は、2次補間、MMSEフィルタリング、ローパスフィルタリング及びスプライン補間等の方法を採用してもよい。線形補間のチャンネル推定方法を採用する場合を例とすれば、先にMMSEフィルタリングされたパイロットサブキャリアにおける周波数領域のチャンネル応答の推定値に対して解析を行うことができる。チャンネルの時間領域に沿った変化がチャンネルの周波数領域に沿った変化よりも早ければ、データチャンネル推定器503は、先に、周波数領域の線形補間を行い、それから、その周波数領域の線形補間結果を利用して時間領域の線形補間を行う。さもなければ、データチャンネル推定器503は、先に、時間領域の線形補間を行い、それから、その時間領域の線形補間結果を利用して周波数領域の線形補間を行う。また、共通参照信号を利用してチャンネル推定を行うLTEシステム等に対して、この実施形態によるチャンネル推定装置は、チャンネルの最大遅延スプレッドを推定する遅延スプレッド推定手段508をさらに備えることができる。この場合に、周波数領域相関性読取手段507は、遅延スプレッド推定手段508により推定された最大遅延スプレッドに基づいてそれに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取る。この実施形態によるチャンネル推定装置に含まれるパイロットチャンネル推定手段501、遅延スプレッド推定手段508、周波数領域相関性読取手段507、移動速度推定手段504、時間領域相関性算出手段505、MMSEフィルタ係数算出手段506、パイロットMMSEフィルタ手段502、及びデータチャンネル推定器503における処理過程のそれぞれは、上記の図4を参照して説明されたチャンネル推定方法に含まれるパイロットLSチャンネル推定ステップS410、最大遅延スプレッド推定ステップS420、チャンネル周波数領域相関性読取ステップS430、移動速度推定ステップS440、チャンネル時間領域相関性算出ステップS450、MMSEフィルタ係数算出ステップS460、パイロットMMSEフィルタリングステップS470、及びデータサブキャリアチャンネル推定ステップS480における処理と同様であるので、簡潔のためその具体的な説明を省略する。次に、図6と図7を参照して他の実施形態によるチャンネル推定方法とチャンネル推定装置の動作原理を説明する。図6は他の実施形態によるチャンネル推定方法における処理のフローチャートを示す。図6に示すように、この実施形態によるチャンネル推定方法は、パイロットLSチャンネル推定ステップS610、最大遅延スプレッド推定ステップS620、チャンネル周波数領域相関性読取ステップS630、移動速度推定ステップS640、チャンネル時間領域相関性算出ステップS650、MMSEフィルタ係数算出ステップS660及びMMSEフィルタリングステップS670等を含む。この実施形態によるチャンネル推定方法に含まれるパイロットLSチャンネル推定ステップS610、最大遅延スプレッド推定ステップS620、チャンネル周波数領域相関性読取ステップS630、移動速度推定ステップS640、及びチャンネル時間領域相関性算出ステップS650における処理過程のそれぞれは、図4を参照して説明されたチャンネル推定方法に含まれるパイロットLSチャンネル推定ステップS410、最大遅延スプレッド推定ステップS420、チャンネル周波数領域相関性読取ステップS430、移動速度推定ステップS440、及びチャンネル時間領域相関性算出ステップS450の処理過程と類似している。そのため、ここで簡潔のためその具体的な説明を省略する。図4を参照して説明されたMMSEフィルタ係数算出ステップS460においてパイロットサブキャリアのMMSEフィルタ係数のみを算出することと異なり、この実施形態によるMMSEフィルタ係数算出ステップS660では、チャンネル時間領域相関性算出ステップS650で算出されたチャンネルの時間領域の相関性とチャンネル周波数領域相関性読取ステップS630で読み取られたチャンネルの周波数領域の相関性を利用してデータサブキャリアのMMSEフィルタ係数を算出する。次に、MMSEフィルタリングステップS670において、MMSEフィルタ係数算出ステップS660で算出されたデータサブキャリアのMMSEフィルタ係数を利用してパイロットLSチャンネル推定ステップS610で出力されたパイロットサブキャリアの最小二乗チャンネル推定に対してMMSEフィルタリングを行ってデータサブキャリアのチャンネル推定を取得する。最後にMMSEフィルタリングステップS670のチャンネル推定結果が出力され、受信機における後続の等化、復調等の処理に用いられる。なお、LTETDDシステムにおいては、端末固有参照信号はマルチパスチャンネルの最大遅延スプレッドの推定に用いることができないので、最大遅延スプレッド推定ステップS620は必須なものではない。この場合に、一つのチャンネルの周波数領域の相関性の表しか予め記憶できない。そこで、最大遅延スプレッド推定ステップS620を省略し、チャンネル周波数領域相関性読取ステップS630で予め記憶されたチャンネルの周波数領域の相関性の予め記憶された表を直接読み出せばよい。また、その他のシステム、例えば共通参照信号を利用してチャンネル推定を行うLTEシステムの場合、共通参照信号の推定結果により周波数領域のチャンネルを時間領域へ変換すると共にノイズを除去する。そして最後に時間領域のチャンネルの推定結果に基づいてチャンネルの最大遅延スプレッド推定を粗略に推定できる。そして、推定された最大遅延スプレッド推定に基づいて、チャンネル周波数領域相関性読取ステップS630において、相応するチャンネルの周波数領域の相関性の予め記憶された表を選定して読み出す。また、この実施形態によるチャンネル推定方法において、先に移動速度推定ステップS640とチャンネル時間領域相関性算出ステップS650を実行し、それから、最大遅延スプレッド推定ステップS620とチャンネル周波数領域相関性読取ステップS630を実行してもよい。チャンネルの周波数領域の相関性と時間領域の相関性を取得できればよい。また、チャンネル周波数領域相関性読取ステップS630で読み取られた、予め記憶されたチャンネルの周波数領域の相関性は、上記の各実施形態の何れかによるチャンネルの周波数領域の相関性の算出方法又はチャンネルの周波数領域の相関性の算出装置により得られたものである。以下に図7を参照して他の実施形態によるチャンネル推定装置を説明する。図7に示すように、この実施形態によるチャンネル推定装置は、パイロットチャンネル推定手段701と、周波数領域相関性読取手段707と、移動速度推定手段704と、時間領域相関性算出手段705と、MMSEフィルタ係数算出手段706と、データチャンネル推定手段703とを備える。パイロットチャンネル推定手段701は、パイロットサブキャリアの最小二乗チャンネル推定を算出する。周波数領域相関性読取手段707は、予め記憶されたチャンネルの周波数領域の相関性を読み取る。この実施形態においても、チャンネルの周波数領域の相関性は、例えば、チャンネル推定装置が有する半導体メモリ、磁気記録媒体あるいは光記録媒体などのメモリに予め記憶される。移動速度推定手段704は、受信側の送信側に対する移動速度を推定する。時間領域相関性算出手段705は、受信側の送信側に対する移動速度に基づいてチャンネルの時間領域の相関性を算出する。MMSEフィルタ係数算出手段706は、算出されたチャンネルの時間領域の相関性と読み取られたチャンネルの周波数領域の相関性を利用してMMSEフィルタ係数を算出する。そしてデータチャンネル推定手段703は、算出されたMMSEフィルタ係数とパイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネル推定を取得する。この実施形態によるチャンネル推定装置において、周波数領域相関性読取手段707により読み取られた上記の予め記憶されたチャンネルの周波数領域の相関性も、上記の各実施形態の何れかによるチャンネルの周波数領域の相関性の算出方法又はチャンネルの周波数領域の相関性の算出装置により算出されたものである。同様に、図5を参照して説明されたMMSEフィルタ係数算出手段506がパイロットサブキャリアのMMSEフィルタ係数のみを算出することと異なり、この実施形態によるMMSEフィルタ係数算出手段706は、時間領域相関性算出手段705により算出されたチャンネルの時間領域の相関性と周波数領域相関性読取手段707により読み取られたチャンネルの周波数領域の相関性を利用してデータサブキャリアのMMSEフィルタ係数を算出する。次に、データチャンネル推定手段703は、MMSEフィルタ係数算出手段706により算出されたデータサブキャリアのMMSEフィルタ係数を利用してパイロットチャンネル推定手段701から出力されたパイロットサブキャリアの最小二乗チャンネル推定に対してMMSEフィルタリングを行ってデータサブキャリアのチャンネル推定を取得する。最後にデータチャンネル推定手段703のチャンネル推定結果が出力され、受信機における後続の等化、復調等の処理に用いられる。同様に、共通参照信号を利用してチャンネル推定を行うLTEシステム等に対して、この実施形態によるチャンネル推定装置は、チャンネルの最大遅延スプレッドを推定する遅延スプレッド推定手段708をさらに備えることができる。この場合に、周波数領域相関性読取手段707は、遅延スプレッド推定手段708により推定された最大遅延スプレッドに基づいて、それに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取る。なお、周波数領域相関性読取手段707により読み取られた、予め記憶されたチャンネルの周波数領域の相関性も、上記の各実施形態の何れかによるチャンネルの周波数領域の相関性の算出方法又はチャンネルの周波数領域の相関性の算出装置により得られたものである。また、この実施形態によるチャンネル推定装置に含まれるパイロットチャンネル推定手段701、遅延スプレッド推定手段708、周波数領域相関性読取手段707、移動速度推定手段704、及び時間領域相関性算出手段705における処理過程のそれぞれは、図5を参照して説明されたチャンネル推定装置に含まれるパイロットチャンネル推定手段501、遅延スプレッド推定手段508、周波数領域相関性読取手段507、移動速度推定手段504、及び時間領域相関性算出手段505の処理過程と同様であるので、ここで簡潔のためその具体的な説明を省略する。図8と図9は、本明細書に開示されたチャンネル推定方法とチャンネル推定装置についてのEPAチャンネルとEVAチャンネルでのシミュレーション結果をそれぞれ示している。このシミュレーションは端末固有参照信号を採用するLTETDDシステムにおいて実行され、仮に推定待ちリソースが1リソースブロックに制限されることを想定している。本アルゴリズムは、チャンネルの周波数領域の相関値を算出する際に、仮にチャンネルの最大遅延スプレッドが48の巡回プレフィックスの長さであるとする。ここで、EPAチャンネルの実際の最大遅延スプレッドは410nsである。48の巡回プレフィックスの長さは約23437nsである。以上の説明からわかるように、チャンネル長さの推定値の偏差が非常に大きいにもかかわらず、チャンネルの推定結果が非常に良い。これにより、本明細書に開示されたチャンネルの周波数領域の相関性の算出方法及びその装置、並びにチャンネル推定方法及びその装置のロバストネスと有効性が十分に実証された。上述の装置における各構成モジュール及び手段は、ソフトウェア、ファームウェア、ハードウェアあるいはそれらの組み合わせによって構成されることができる。その構成に使用可能な具体的な手法又は方式は当業者に熟知されているものであり、ここでは説明しない。ソフトウェア又はファームウェアにより実現される場合には、記憶媒体又はネットワークから専用のハードウェア構成を持つコンピュータ例えば図10に示した汎用コンピュータ1100に対してこのソフトウェアを構成するプログラムをインストールし、各種プログラムがインストールされると、そのコンピュータは各種機能等を実行することができる。なお、上記の各実施形態による装置は、例えば、OFDM方式により無線信号を送信し、あるいは受信する基地局または携帯端末のような移動局に組み込まれる。図10においては、中央処理装置CPU1101が、読取専用メモリROM1102に記憶されたプログラム又はメモリ部1108からランダムアクセスメモリRAM1103にアップロードされたプログラムに従って、各種の処理を実行する。RAM1103においても、必要に応じてCPU1101が各種の処理を実行するとき等に必要なデータが記憶されている。CPU1101、ROM1102とRAM1103同士はバス504を介して接続されている。入力/出力インタフェース1105もパス1104に接続されている。入力部1106キーボード、マウス等を含むと、出力部1107ディスプレイ、例えばブラウン管CRT、液晶ディスプレイLCD等とスピーカ等を含むと、メモリ部1108ハードディスク等を含むと、通信部1109ネットワークインターフェースカード、例えばLANカード、モデム等を含むとは、入力/出力インタフェース1105に接続されている。通信部1109がネットワーク、例えばインターネットを経由して通信処理を実行する。必要に応じて、入力/出力インタフェース1105にはドライブ1110も接続されている、例えば磁気ディスク、光ディスク、光磁気ディスク、半導体メモリ等のような取り外し可能な媒体1111は、必要に応じてドライブ1110に取り付けられており、その中から読み出されたコンピュータプログラムが必要に応じてメモリ部1108にインストールされるようにする。ソフトウェアで上記の一連の処理を実現する場合、ネットワーク例えばインターネット、又は記憶媒体例えば取り外し可能な媒体1111からソフトウェアを構成するプログラムがコンピュータ1100にインストールされてもよい。このような記憶媒体は、図10に示したような、その中にプログラムが記憶されているものであって、デバイスと分離して配送されることでユーザにプログラムを提供する取り外し可能な媒体1111には限定されないことを、当業者は理解すべきである。取り外し可能な媒体1111の例として、磁気ディスクフロッピディスク登録商標含む、光ディスクコンパクトディスクリードオンリーメモリCD-ROMやディジタルヴァーサタイルディスクDVDを含む、光磁気ディスクミニディスク登録商標含む及び半導体メモリを含む。または、記憶媒体はROM1102、メモリ部1108に含まれるハードディスクなどでもよい。その中にプログラムが記憶されており、且つこれらを含むデバイスと一緒にユーザに配送される。本発明は更に、機器で読取可能な命令コードが記憶されたプログラムを提供する。上記の命令コードが機器で読み取られて実行されると、その機器は、上記の実施形態の何れかの方法を実行することができる。それに対して、上記の機器で読み取り可能な命令コードを含むプログラムを記憶する記憶媒体も本明細書の開示に含まれる。この記憶媒体は、フロッピーディスク登録商標、光ディスク、光磁気ディスク、メモリカード、メモリスティック等を含むが、これらに限定されない。以上の具体的な実施形態の説明において、一種の実施形態について記述し及び/又は示した特徴は、同一あるいは類似の形態で一つ又は複数の他の実施形態で使用されたり、他の実施形態における特徴と組み合わせたり、あるいは他の実施形態における特徴を代替することができる。なお、専門用語含む/有するが本文で使用される場合には、特徴、要素、ステップ又はコンポーネントの存在を意味するが、一つ又は複数の他の特徴、要素、ステップ又はコンポーネントの存在又は付加を排除しない。また、本明細書に開示された方法は明細書において記述した時間順に実行することには限られず、その他の時間順序に従って、並行にあるいは個別に実行されてもよい。従って、本明細書で記述した方法の実行順序は本発明の技術的範囲への限定を構成しない。上記に本発明の具体的実施例を記述することにより本発明を開示したが、上記の全ての実施例と例示はいずれも例示的なものであり、限定的なものではないことを理解すべきである。当業者は添付の特許請求の範囲の精神と範囲に、本発明に対する各種の修正、改良或は均等物を設計することができる。これらの修正、改良或は均等物は本発明の保護範囲内に含まれるものであると考えられるべきである。以上説明した実施形態及びその変形例に関し、更に以下の付記を開示する。付記1チャンネルの周波数領域の相関性の算出方法であって、サンプルチャンネルのパス遅延分布を解析し、前記サンプルチャンネルのパスエネルギー分布を解析し、前記サンプルチャンネルの前記パス遅延分布と前記パスエネルギー分布に基づいてチャンネルモデルを構築し、構築された前記チャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する、ことを含む方法。付記2遅延スプレッドに基づいて複数の解析区間を選定することをさらに含み、前記パス遅延分布の解析は、前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パス遅延分布を解析し、前記パスエネルギー分布の解析は、前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パスエネルギー分布を解析する、付記1に記載の方法。付記3前記複数の解析区間のそれぞれにおいて複数のパスの遅延は、均一分布または等差逓増分布を呈し、前記複数の解析区間のそれぞれにおいて前記複数のパスのエネルギー分布は、線形逓減分布、指数逓減分布或は等エネルギー均一分布を呈する、付記2に記載の方法。付記4前記チャンネルの周波数領域の相関性の算出は、異なる最大遅延スプレッドに対して、下記の式により、構築された前記チャンネルモデルに基づいて前記チャンネルの周波数領域の相関性をそれぞれ算出し、ここで、θΔfkkmはチャンネルの周波数領域の相関性関数であり、absは絶対値の計算操作であり、kkmはチャンネルの周波数領域の相関値を算出する両周波数点間のサブキャリア間隔数であり、iはパスのシーケンス番号であり、σi2はi番目のパスのパスエネルギーであり、τiはi番目のパスのパス遅延である、付記3に記載の方法。付記5チャンネルの周波数領域の相関性の算出装置であって、サンプルチャンネルのパス遅延分布を解析するパス遅延解析手段と、前記サンプルチャンネルのパスエネルギー分布を解析するパスエネルギー解析手段と、前記サンプルチャンネルの前記パス遅延分布と前記パスエネルギー分布に基づいてチャンネルモデルを構築するチャンネルモデル構築手段と、構築された前記チャンネルモデルに基づいてチャンネルの周波数領域の相関性を算出する周波数領域相関性算出手段と、を備える算出装置。付記6遅延スプレッドに基づいて複数の解析区間を選定する解析区間選択手段をさらに備え、前記パス遅延解析手段は前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パス遅延分布を解析し、前記パスエネルギー解析手段は前記複数の解析区間のそれぞれごとに前記サンプルチャンネルの前記パスエネルギー分布を解析する、付記5に記載の算出装置。付記7前記チャンネルモデル構築手段は、前記複数の解析区間のそれぞれにおいて複数のパスの遅延が均一分布又は等差逓増分布を呈し、且つ前記複数の解析区間のそれぞれにおいて前記複数のパスのエネルギー分布が線形逓減分布、指数逓減分布または等エネルギー均一分布を呈する原則に従って、前記チャンネルモデルを構築する付記6に記載の算出装置。付記8前記周波数領域相関性算出手段は、異なる最大遅延スプレッドに対して、下記の式により、前記チャンネルモデル構築手段により構築されたチャンネルモデルに基づいてチャンネルの周波数領域の相関性をそれぞれ算出し、ここで、θΔfkkmはチャンネル周波数領域の相関性関数であり、absは絶対値の計算操作であり、kkmはチャンネルの周波数領域の相関値を算出する両周波数点間のサブキャリア間隔数であり、iはパスのシーケンス番号であり、σi2はi番目のパスのパスエネルギーであり、τiはi番目のパスのパス遅延である、付記7に記載の算出装置。付記9チャンネル推定方法であって、パイロットサブキャリアの最小二乗チャンネル推定を算出し、付記1~4の何れか一項に記載のチャンネルの周波数領域の相関性の算出方法により得られた、予め記憶されたチャンネルの周波数領域の相関性を読み取り、受信側の送信側に対する移動速度を推定し、前記受信側の前記送信側に対する前記移動速度に基づいてチャンネルの時間領域の相関性を算出し、算出された前記チャンネルの時間領域の相関性と読み取られた前記チャンネルの周波数領域の相関性を利用して最小平均二乗誤差フィルタ係数を算出し、前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定する、ことを含む方法。付記10チャンネルの最大遅延スプレッドを推定し、かつ推定された前記最大遅延スプレッドに基づいて、それに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取る、ことをさらに含む付記9に記載のチャンネル推定方法。付記11前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用して前記最小平均二乗誤差フィルタ係数を算出することは、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用して前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定することは、前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用して前記パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行い、かつ前記最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいて前記データサブキャリアのチャンネルに対して推定を行って前記データサブキャリアのチャンネルを推定することを含む、付記10に記載の方法。付記12前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用して前記最小平均二乗誤差フィルタ係数を算出することは、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用してデータサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネル推定を取得することは、前記データサブキャリアの最小平均二乗誤差フィルタ係数を利用して、前記パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行って前記データサブキャリアのチャンネルを推定する、付記10に記載のチャンネル推定方法。付記13パイロットサブキャリアの最小二乗チャンネル推定を算出するパイロットチャンネル推定手段と、付記5~8の何れか一項に記載のチャンネルの周波数領域の相関性の算出装置により得られた、予め記憶されたチャンネルの周波数領域の相関性を読み取る周波数領域相関性読取手段と、受信側の送信側に対する移動速度を推定する移動速度推定手段と、前記受信側の前記送信側に対する前記移動速度に基づいてチャンネルの時間領域の相関性を算出する時間領域相関性算出手段と、算出された前記チャンネルの時間領域の相関性と読み取られた前記チャンネルの周波数領域の相関性を利用して最小平均二乗誤差ファイル係数を算出するフィルタ係数算出手段と、算出された前記最小平均二乗誤差フィルタ係数と前記パイロットサブキャリアの最小二乗チャンネル推定に基づいてデータサブキャリアのチャンネルを推定するデータチャンネル推定手段と、を備えるチャンネル推定装置。付記14チャンネルの最大遅延スプレッドを推定する遅延スプレッド推定手段をさらに備え、前記周波数領域相関性読取手段は、推定された前記最大遅延スプレッドに基づいて、それに相応する、予め記憶されたチャンネルの周波数領域の相関性を読み取る、付記13に記載のチャンネル推定装置。付記15前記フィルタ係数算出手段は、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用してパイロットサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記データチャンネル推定手段は、前記パイロットサブキャリアの最小平均二乗誤差フィルタ係数を利用して前記パイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行うパイロット最小平均二乗誤差フィルタ手段と、前記最小平均二乗誤差フィルタリングを行ったパイロットサブキャリアに基づいて前記データサブキャリアのチャンネルに対して推定を行って前記データサブキャリアのチャンネルを推定するデータチャンネル推定器とを備える付記14に記載のチャンネル推定装置。付記16前記フィルタ係数算出手段は、前記チャンネルの時間領域の相関性と前記チャンネルの周波数領域の相関性を利用してデータサブキャリアの最小平均二乗誤差フィルタ係数を算出し、前記データチャンネル推定手段は、算出されたデータサブキャリアの前記最小平均二乗誤差フィルタ係数を利用してパイロットサブキャリアの最小二乗チャンネル推定に対して最小平均二乗誤差フィルタリングを行ってデータサブキャリアのチャンネルを推定する、付記14に記載のチャンネル推定装置。