"参照信号が挿入された受信信号を、前記参照信号のレプリカで複素除算することによりチャネル応答推定値を算出して伝搬路の状態を推定するチャネル推定装置であって、前記複素除算により算出されたチャネル応答推定値を、IDFT処理により時間領域のチャネル応答推定値に変換するIDFT部と、前記IDFT処理後の所望波信号のチャネル応答推定値に窓関数を乗算する窓関数乗算部と、前記窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換するDFT部と、前記複素除算により算出されたチャネル応答推定値から置き換えデータを算出するサンプル補正部と、前記DFT部から出力されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、前記サンプル補正部が算出した前記置き換えデータで置き換えるサンプル置換部と、を備えることを特徴とするチャネル推定装置。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータを前記置き換えデータとすることを特徴とする請求項1に記載のチャネル推定装置。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータの平均値を前記置き換えデータとすることを特徴とする請求項1に記載のチャネル推定装置。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータの平均値を算出し、算出した平均値に忘却係数を乗算したデータを前記置き換えデータとすることを特徴とする請求項1に記載のチャネル推定装置。参照信号が挿入された受信信号を、前記参照信号のレプリカで複素除算することによりチャネル応答推定値を算出して伝搬路の状態を推定するチャネル推定方法であって、前記複素除算により算出されたチャネル応答推定値を、IDFT処理により時間領域のチャネル応答推定値に変換するステップと、前記IDFT処理後の所望波信号のチャネル応答推定値に窓関数を乗算するステップと、前記窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換するステップと、前記複素除算により算出されたチャネル応答推定値から置き換えデータを算出するステップと、DFT処理後のチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、前記置き換えデータで置き換えるステップと、を含むことを特徴とするチャネル推定方法。しかしながら、実際のシステムにおいて、送受信に使用するIFFT/FFT処理のポイント数と、データ送受信に使用するサブキャリア数とが等しくないことから、図5に示すように、窓関数を乗算する際、実際には少なからず、所望波信号成分有効な遅延波成分をも除去することになる。そのため、上述したIDFT/DFT処理を用いるチャネル推定では、図6に示すように、DFT処理後においてチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分に歪みが生じるため、チャネル応答の推定精度が低下する。また、上述したIDFT/DFT処理を用いるチャネル推定では、図7に示すように、窓関数対象時間内に雑音がある場合は、その影響により、チャネル応答の推定精度が低下する。特に、図8のように、遅延波を含む受信信号で遅延が0の所望波の信号レベルが高いとき遅延の影響で所望波が窓関数対象時間からはみ出したとき、IFFT/FFT処理のポイント数と、データ送受信に使用するサブキャリア数が等しくないことから発生する、インパルス応答の折り返し図中右端の情報は、窓関数を乗算することで完全に消えてしまうため、チャネル応答の推定精度が低下する。本発明は、このような問題点に鑑みてなされたものであり、本発明の目的は、チャネル応答の推定精度を向上させることのできるチャネル推定装置およびチャネル推定方法を提供することにある。上記目的を達成するため、本発明のチャネル推定装置は、参照信号が挿入された受信信号を、前記参照信号のレプリカで複素除算することによりチャネル応答推定値を算出して伝搬路の状態を推定するチャネル推定装置であって、前記複素除算により算出されたチャネル応答推定値を、IDFT処理により時間領域のチャネル応答推定値に変換するIDFT部と、前記IDFT処理後の所望波信号のチャネル応答推定値に窓関数を乗算する窓関数乗算部と、前記窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換するDFT部と、前記複素除算により算出されたチャネル応答推定値から置き換えデータを算出するサンプル補正部と、前記DFT部から出力されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、前記サンプル補正部が算出した前記置き換えデータで置き換えるサンプル置換部とを備えることを特徴とする。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータを前記置き換えデータとすることが好ましい。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータの平均値を前記置き換えデータとすることが好ましい。前記サンプル補正部は、前記複素除算により算出されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、抽出したデータの平均値を算出し、算出した平均値に忘却係数を乗算したデータを前記置き換えデータとすることが好ましい。また、本発明のチャネル推定方法は、参照信号が挿入された受信信号を、前記参照信号のレプリカで複素除算することによりチャネル応答推定値を算出して伝搬路の状態を推定するチャネル推定方法であって、前記複素除算により算出されたチャネル応答推定値を、IDFT処理により時間領域のチャネル応答推定値に変換するステップと、前記IDFT処理後の所望波信号のチャネル応答推定値に窓関数を乗算するステップと、前記窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換するステップと、前記複素除算により算出されたチャネル応答推定値から置き換えデータを算出するステップと、DFT処理後のチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、前記置き換えデータで置き換えるステップとを含むことを特徴とする。本発明は、チャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分に発生する歪みを解消することができるので、チャネル応答の推定精度を向上させることができる。本発明の実施の形態について、図面を参照して説明する。図9は、本発明のチャネル推定装置が用いられるマルチキャリア伝送システム例えばOFDMにおける基地局受信機の構成を示す図であり、2本アンテナの例を示している。図9に示す基地局受信機では、無線受信部11-1、11-2にて信号を受信し、A/D変換部12-1、12-2にて受信信号をアナログ/デジタル変換する。CP除去部13-1、13-2は、デジタル変換された受信信号からCPサイクリックプレフィックス:CyclicPrefixを除去し、FFT部14-1、14-2は、CP除去後の受信信号を高速フーリエ変換FFTする。FFT部14-1、14-2の出力する受信信号は2分配され、一方がチャネル推定部15-1、15-2へ、他方がチャネル等化部16へ入力される。チャネル推定部チャネル推定装置15-1、15-2は、参照信号のレプリカを用いて、FFT部14-1、14-2の出力する受信信号からチャネル応答を推定する。チャネル等化部16は、チャネル推定部15-1、15-2で推定されたチャネル応答を用いて、FFT部14-1、14-2から出力される受信信号に対してチャネル等化を行う。サブキャリアデマッピング部17は、チャネル等化後の信号をデマッピングし、復調部18は、デマッピング後の信号を復調して受信データを算出する。図10は、本発明の第1の実施形態に係るチャネル推定装置の構成図である。第1の実施形態に係るチャネル推定装置は、図6に示すチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分における歪みを解消するために、帯域両端の高域部分と低域部分の値を複素除算で算出されたチャネル応答推定値で置き換えるものである。第1の実施形態に係るチャネル推定装置は、参照信号レプリカ生成部1と、複素除算部2と、IDFT部3と、窓関数乗算部4と、DFT部5と、サンプル補正部6と、サンプル置換部7とを備えている。参照信号レプリカ生成部1は、送信機側で送信信号に挿入された参照信号のレプリカを生成し、複素除算部2は、受信信号を参照信号のレプリカで複素除算することによりチャネル応答推定値を算出する。IDFT部3は、複素除算部2で算出した周波数領域のチャネル応答推定値を、IDFT処理により時間領域のチャネル応答推定値に変換する。窓関数乗算部4は、IDFT処理後の所望波信号のチャネル応答推定値に窓関数矩形窓でもその他の窓でもよいを乗算する。DFT部5は、窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換する。サンプル補正部6は、複素除算部2で算出されたチャネル応答推定値から置き換えデータを算出し、サンプル置換部7は、DFT部5から出力されたチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、サンプル補正部6が算出した置き換えデータで置き換える。図11は、複素除算部から出力されたチャネル応答推定値を示す図である。サンプル補正部6は、複素除算部2から出力されたチャネル応答推定値のうち、データ割当帯域両端の高域部分と低域部分から置き換えポイント数分のチャネル応答推定値を抽出して置き換えデータとし、サンプル置換部7は、図6に示すデータ割当帯域両端の高域部分と低域部分の歪みのあるデータを、この置き換えデータで置き換える。図12は、サンプル置換部から出力されるチャネル応答推定値を示す図である。サンプル置換部7からは、データ割当帯域両端の高域部分と低域部分において歪みのないチャネル応答推定値を出力することができる。サンプル補正部6は、複素除算部2で算出されたチャネル応答推定値から、以下の方法で置き換えデータを算出する。なお、ポイント数は、システムの帯域と割当データの帯域から可変とする。第1の方法は、データ割当帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、そのデータをそのまま置き換えデータとする方法であり、第2の方法は、データ割当帯域両端の高域部分と低域部分から置き換えポイント数分のデータを抽出し、それぞれ平均値を求めて、そのデータを置き換えデータとする方法であり、第3の方法は、第2の方法で求めた平均値に忘却係数を乗算して求めたデータを置き換えデータとする方法である。例えば、ポイント数を3、サンプル置換部7の出力をy’,y’,y’、複素除算部2の出力をy,y,y、忘却係数を0<α≦1とすると、第1の方法で置き換えデータを算出した場合、サンプル置換部7の出力は、となり、第2の方法で置き換えデータを算出した場合、サンプル置換部7の出力は、となり、第3の方法で置き換えデータを算出した場合、サンプル置換部7の出力はとなる。上述したように、本発明の第1の実施形態に係るチャネル推定装置は、チャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分に発生する歪みを解消することができるので、チャネル応答の推定精度を向上させることができる。図13は、本発明の第2の実施形態に係るチャネル推定装置の構成図である。第2の実施形態に係るチャネル推定装置は、図7に示すように、窓関数対象時間内に雑音がある場合に、その影響により、チャネル応答の推定精度が低下し、図8に示すように、インパルス応答の折り返し図中右端の情報が、窓関数を乗算することで消失してしまうことにより、チャネル応答の推定精度が低下するという問題に対応するために、受信信号のうち、一番レベルの強いインパルス成分所望波遅延波を特定し、その前後の数サンプルを窓関数にて乗算するものである。第2の実施形態に係るチャネル推定装置は、参照信号レプリカ生成部1と、複素除算部2と、IDFT部3と、窓関数乗算部4と、DFT部5と、ピークサーチ部8と、窓スライド部9とを備えている。参照信号レプリカ生成部1は、送信機側で送信信号に挿入された参照信号のレプリカを生成し、複素除算部2は、受信信号を参照信号のレプリカで複素除算することによりチャネル応答推定値を算出する。IDFT部3は、複素除算部2で算出された周波数領域のチャネル応答推定値をIDFT処理により時間領域のチャネル応答推定値に変換する。ピークサーチ部8は、IDFT部3の出力であるチャネル応答推定値のうち、所定の範囲内で最大レベルのインパルス成分所望遅延波をサーチする。ここで、所定の範囲とは、可変であるが、一般的には、CP長であるとする。窓スライド部9は、ピークサーチ部8で特定した最大レベルのインパルス成分のサンプル点から、前後数サンプルの窓関数矩形窓でもその他の窓でも良いを導出する。このときの前後サンプル数は、マルチキャリア伝送システムの帯域幅や、CP長、送信機と受信機の距離等から可変とする。窓関数乗算部4は、窓スライド部9で導出した窓関数をIDFT部3から出力されたチャネル応答推定値に乗算する。DFT部5は、窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換する。図14は、窓関数対象時間内に雑音がある場合に、ピークサーチ部8により所定の範囲内で最大レベルのインパルス成分をサーチし、窓スライド部9により、ピークサーチ部8で特定した最大レベルのインパルス成分のサンプル点から、前後数サンプルの窓関数を導出窓関数対象時間を図14のようにやや右側にスライドし、算出した窓関数を、窓関数乗算部4により、IDFT処理後のチャネル応答推定値に適用する場合を示す図である。図15は、インパルス応答の折り返し図中右端の情報がある場合に、ピークサーチ部8により所定の範囲内で最大レベルのインパルス成分をサーチし、窓スライド部9により、ピークサーチ部8で特定した最大レベルのインパルス成分のサンプル点から、前後数サンプルの窓関数を導出窓関数対象時間を図15のように図中左端、図中右端に設定し、算出した窓関数を、窓関数乗算部4により、IDFT処理後のチャネル応答推定値に適用する場合を示す図である。上述したように、本発明の第2の実施形態に係るチャネル推定装置は、窓関数対象時間内に不要な雑音成分が混在しないように窓関数を導出することができるので、チャネル応答の推定精度を向上させることができる。図16は、本発明の第3の実施形態に係るチャネル推定装置の構成図である。第3の実施形態に係るチャネル推定装置は、第1の実施形態のチャネル推定装置と第2の実施形態のチャネル推定装置とを組み合わせたものである。第3の実施形態に係るチャネル推定装置は、参照信号レプリカ生成部1と、複素除算部2と、IDFT部3と、窓関数乗算部4と、DFT部5と、サンプル補正部6と、サンプル置換部7と、ピークサーチ部8と、窓スライド部9とを備えている。参照信号レプリカ生成部1は、送信機側で送信信号に挿入された参照信号のレプリカを生成し、複素除算部2は、受信信号を参照信号のレプリカで複素除算することによりチャネル応答推定値を算出する。IDFT部3は、複素除算部2で算出した周波数領域のチャネル応答推定値をIDFT処理により時間領域のデータに変換する。ピークサーチ部8は、IDFT部3の出力であるチャネル応答推定値のうち、所定の範囲内で最大レベルのインパルス成分所望遅延波をサーチする。窓スライド部9は、ピークサーチ部8で特定した最大レベルのインパルス成分のサンプル点から、前後数サンプルの窓関数矩形窓でもその他の窓でも良いを導出する。窓関数乗算部4は、窓スライド部9で導出した窓関数をIDFT部3から出力されたチャネル応答推定値に乗算する。DFT部5は、窓関数が乗算されたチャネル応答推定値を、DFT処理により周波数領域のチャネル応答推定値に変換する。サンプル補正部6は、複素除算部2で算出されたチャネル応答推定値から置き換えデータを算出し、サンプル置換部7は、DFT部5から出力されるチャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分のデータを、サンプル補正部6が算出した置き換えデータで置き換える。本発明の第3の実施形態に係るチャネル推定装置は、チャネル応答推定値のデータが割り当てられている周波数帯域両端の高域部分と低域部分に発生する歪みを解消することができるので、チャネル応答の推定精度を向上させることができる。また、本発明の第3の実施形態に係るチャネル推定装置は、窓関数対象時間内に不要な雑音成分が混在しないように窓関数を導出することができるので、チャネル応答の推定精度を向上させることができる。"