JP5670024B2 - チャネル推定装置及びその方法並びにプログラム - Google Patents

Description

本発明はチャネル推定装置及びその方法並びにプログラムに関し、特に無線通信システムにおける空間伝搬路のチャネル推定をなす装置及びその方法並びにプログラムに関するものである。

無線通信の大容量化及び高速化に伴い、有限な周波数資源を有効に利用することが不可欠であり、その一つの技術として、空間領域を利用する手法がある。その手法の一例として、複数のアンテナを利用したＭＩＭＯ（Multi Input Multi Output）技術が開発され、実用化が始まっている。このＭＩＭＯ技術は、例えば、複数のデータを複数のアンテナから同時に送信することにより、空間多重を行い、受信側では、受信データから元の複数のデータを取り出す技術である。同時刻に送信するデータ量を増やすことができるので、通信速度が向上することになる。

また変調方式も通信速度に大きく関係しており、例えば、ＱＰＳＫ変調方式では２ビット、１６ＱＡＭ変調方式では４ビット、６４ＱＡＭ変調方式では６ビットと、いずれも同時に送信可能な情報ビット数を増やすことができる。但し、同時に送信可能なビット数が増えれば、それだけ雑音耐性が低下するので、符号化率や符号化方法と合わせて、伝搬路環境の品質に応じた最適な方式を選択するＡＭＣ（Adaptive Modulation and Coding）という技術が開発され用いられている。

移動体通信においても、現在、３ＧＰＰ（Third Generation Pertnership Project）によりＬＴＥ（Long Term Evolution ）と称される新しい通信方式での標準化が進められており、このＭＩＭＯ技術やＡＭＣ技術が採用されている。

このように、送信機及び受信機が複数のアンテナを備えて、空間多重伝送を行う時、Ｓ／Ｎ（信号対雑音比）が十分大きい環境下においては、アンテナ数に比例した通信容量の拡大が可能となる。かかる空間多重伝送方式として、直交周波数分割多重（ＯＦＤＭ：Orthogonal Frequency Multiplexing ）を用いたマルチキャリア変調方式が用いられる。

このマルチキャリア変調方式は、複数のサブキャリアを用いる伝送方式であり、ＯＦＤＭは、各サブキャリア周波数が直交関係にあり、高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）技術を用いて周波数の異なるサブキャリア信号を一括変換処理するものである。そして、サブキャリア信号を時間軸の信号に変換した後、送信するようになっている。

このようなＯＦＤＭ信号の受信に際しては、マルチパスフェージング環境下におけるチャネル特性の高精度な推定が不可欠である。このチャネル推定方法としては、既知パターンであるパイロットシンボルを用いる方法が一般的に利用されている。例えば、先述した３ＧＰＰ ＬＴＥでは、パイロットシンボルの構成として、図９に一例を示す様に、サブフレーム内でパイロットシンボルを時間軸方向および周波数方向に二次元的に分散配置する、いわゆるスキャッタードマッピング（scattered mapping ）が提案されている。

かかるパイロット信号がＲＳ（Reference Signal）信号として、送信アンテナ毎に送信されるようになっている。なお、送信すべきデータは、ＲＳ信号が配置されていない位置に配置されることになる。

図９において、送信アンテナ＃０のＲＳ信号は斜線、送信アンテナ＃１のＲＳ信号は縦線、送信アンテナ＃２のＲＳ信号は黒塗り、送信アンテナ＃３のＲＳ信号は横線でそれぞれ示されている。なお、図９に示したＲＳ信号の配置は、単に一例を示すものである。

ＲＳ信号が配置されている位置のチャネル推定方法の例としては、例えば、特許文献１に開示されているように、受信データをＲＳ信号のレプリカで除算する方法がある。この方法においては、データが配置されている部分のチャネル推定は、例えば、ＲＳ信号により求めたチャネル推定値を用いて補間処理を行うことによって算出するようになっている。

更に、チャネル推定精度を高める方法として、受信データに含まれる雑音成分を効率良く抑圧する技術が特許文献２に示されている。この技術によれば、周波数方向のデータをＩＦＦＴ（Inverse Fast Fourier Transform）技術により時間方向のデータに戻す。そして、この時間方向で、不要な成分や雑音と推定される成分などを、例えば、ゼロデータに置換して、ＦＦＴ処理により周波数方向のデータに戻すことにより、雑音などの抑圧を図っている。

ここで、周波数方向でみてＲＳ信号間やＲＳ信号端において、周波数方向あるいは時間方向の補間により仮推定値を生成し、これを雑音抑圧処理に含めることにより、雑音抑圧の恩恵をより広く享受する方法がある。例えば、図１０に示す送信アンテナ＃０のＲＳ信号間における仮推定値を求める場合について説明する。ここで、仮推定値とは、ＲＳ信号の位置以外のデータ位置のチャネル推定値を求めるための仮の推定値である。

ＲＳ信号間（図１０の黒塗りで示す部分）における仮推定値をｈ_tmp （０，ｂ，ｎ，０，０）とすると、周波数方向補間の場合は、ｈ_rs（０，ｂ，ｎ，０，０）とｈ_rs（０，ｂ，ｎ，０，１）とを用いてｈ_tmp を生成し、また時間方向補間の場合は、ｈ_rs（０，ｂ，ｎ−１，１，０）とｈ_rs（０，ｂ，ｎ，１，０）とを用いてｈ_tmp を生成する。

ここで、ＲＳ信号が配置されている位置のチャネル推定値をｈ_rs（ａ，ｂ，ｎ，ｒ，ｋ）とし、仮推定値をｈ_tmp （ａ，ｂ，ｎ，ｒ，ｋ）とする。ここに、ａは送信アンテナ番号（０〜３）、ｂは受信アンテナ番号（０〜３）、ｎはサブフレームを構成するスロットの番号、ｒはＲＳ信号の種別を示すＲＳ種（０：１ｓｔＲＳ，１：２ｎｄＲＳ）、ｋはＲＳ番号を、それぞれ示す。なお、ｈ_rsおよびｈ_tmp は複素数である。

周波数方向及び時間方向の各仮推定値ｈ_tmp （０，ｂ，ｎ，０，０）は、以下の式（１）および式（２）でそれぞれ示されることになる。なお、Ｗ_f （ｎ，ｉ）は周波数方向の内挿補間係数を示し、Ｗ_t （ｎ，ｉ）は時間方向の内挿補間係数を示す。ここに、ｉは２点間の補間における各点を特定するための番号である。

ｈ_tmp （０，ｂ，ｎ，０，０）＝Ｗ_f （ｎ，０）＊ｈ_rs（０，ｂ，ｎ，０，０）＋Ｗ_f （ｎ，１）＊ｈ_rs（０，ｂ，ｎ，０，１）……（１）

ｈ_tmp （０，ｂ，ｎ，０，０）＝Ｗ_t （ｎ，０）＊ｈ_rs（０，ｂ，ｎ−１，１，０）＋Ｗ_t （ｎ，１）＊ｈ_rs（０，ｂ，ｎ，１，０）……（２）

また、図１１に送信アンテナ＃２のＲＳ信号間における仮推定値ｈ_tmp （２，ｂ，ｎ，０，０）を求める場合について示している。この場合の周波数方向及び時間方向の各仮推定値は、式（３）および（４）でそれぞれ示される。

ｈ_tmp （２，ｂ，ｎ，０，０）＝Ｗ_f （ｎ，０）＊ｈ_rs（２，ｂ，ｎ，０，０）＋Ｗ_f （ｎ，１）＊ｈ_rs（２，ｂ，ｎ，０，１）……（３）

ｈ_tmp （２，ｂ，ｎ，０，０）＝Ｗ_t （ｎ，０）＊ｈ_rs（２，ｂ，ｎ−１，０，０）＋Ｗ_t （ｎ，１）＊ｈ_rs（２，ｂ，ｎ＋１，０，０）……（４）

このように、チャネル推定では、先ず、パイロット信号であるＲＳ信号の位置におけるチャネル推定値ｈ_rsを算出し、次に、このＲＳ信号の位置におけるチャネル推定値ｈ_rsを用いて、式（１）〜（４）などにより得られたＲＳ信号の位置以外のデータ位置の周波数方向および時間方向の仮推定値ｈ_tmp を算出する。

その後、これら周波数方向および時間方向の仮推定値ｈ_tmp のうち一つを適宜選択して、周波数方向または時間方向の補間処理をなすことにより、ＲＳ信号の位置以外のチャネル推定値を求める。こうしてチャネル推定値が求まると、チャネル推定値とパイロット信号を除いた受信データとを用いて等化処理を行う。その後、この等化信号を復調して復調データを導出するのである。

特開２００８−１６０８２２号公報 特開２００８−１６７０８８号公報

ＲＳ信号の位置以外のデータ位置のチャネル推定値を求めるためには、上述したように、ＲＳ信号の位置におけるチャネル推定値から求めた周波数方向および時間方向の仮推定値のうち一つを選択して、周波数方向で補間するか、時間方向で補間するかの選択をする必要がある。そのためには、図１０や図１１のＢで示す時間まで、受信データをバッファリングしておくことが必要である。よって、バッファリングのためのメモリ容量が増大し、また、仮推定値の生成もＢで示す時間まで遅延することになる。

本発明の目的は、チャネル推定精度を維持しつつ、仮推定値の推定処理時間を短縮してチャネル推定処理完了までの処理時間を短くすると共に、受信データのバッファリングのためのメモリ容量の増大を防止可能なチャネル推定装置及びその方法並びにプログラムを提供することである。

本発明によるチャネル推定装置は、空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定装置であって、前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定手段を有し、前記仮推定手段は、前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間手段と、前記一次推定値に対して過去のチャネル推定値による忘却処理を行って前記仮推定値を算出する忘却処理手段とを有し、前記忘却処理後の仮推定値に基づいて前記データ位置のチャネル推定をなすことを特徴とする。

本発明によるチャネル推定方法は、空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定方法であって、前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定ステップを有し、前記仮推定ステップは、前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間ステップと、前記一次推定値に対して過去のチャネル推定値による忘却処理を行って前記仮推定値を算出する忘却処理ステップとを有し、前記忘却処理後の仮推定値に基づいて前記データ位置のチャネル推定をなすことを特徴とする。

本発明によるプログラムは、空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定方法をコンピュータに実行させるためのプログラムであって、前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定処理を有し、前記仮推定処理は、前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間処理と、前記一次推定値に対して過去のチャネル推定値による忘却処理を行って前記仮推定値を算出する忘却処理とを有し、前記忘却処理後の仮推定値に基づいて前記データ位置のチャネル推定をなすことを特徴とする。

本発明によれば、データ位置のチャネル推定値を算出するための仮の推定値の算出において、先ず、周波数方向に内挿補間し、この補間結果に対して、過去のチャネル推定値による忘却処理を行うようにしたので、周波数方向補間と時間方向補間との選択のために、時間方向で少なくとも一つ後の基準信号を待つ必要がなくなり、遅延時間をゼロとすることができる。また、時間方向の忘却処理を行うことにより、チャネル推定精度の向上も可能となる。

本発明の実施の形態が適用される無線通信システムの概略図である。 図１の受信側における復調部の例を示すブロック図である。 図２の仮推定回路３２の一例を示すブロック図であり、本発明の一実施の形態を示す図である。 図３の仮推定回路３２の動作例を示すフローチャートである。 本発明の一実施の形態の動作を説明するための図である。 本発明に関連する技術の機能ブロック図の例である。 図２の仮推定回路３２の他の例を示すブロック図であり、本発明の他の実施の形態を示す図である。 本発明の他の実施の形態の動作を説明するための図である。 ３ＧＰＰ ＬＴＥにより提案されているスキャッタードマッピングの例を示す図である。 本発明に関連する技術の一例を説明する図である。 本発明に関連する技術の他の例を説明する図である。

以下に図面を参照しつつ本発明の実施の形態について説明する。図１は、本発明の実施の形態が適用される通信システムの全体構成図である。図１に示すように、送信側は、例えば基地局１０であり、符号化部１１、変調部１２、Ｄ／Ａ（デジタル／アナログ）変換部１３および複数のアンテナ１４を有する。受信側は、例えば携帯電話機などのユーザ端末２０であり、複数のアンテナ２４、Ａ／Ｄ（アナログ／デジタル）変換部２３、復調部２２および復号化部２１を有する。

送信側の基地局１０では、例えば、先ず基地局１０の図示せぬＣＰＵ（制御部）が送信したいデータを情報ビットとして符号化部１１へ入力する。符号化部１１は、入力された情報ビットに対しＣＲＣ（Cyclic Redundancy Check ）コードの付加や、畳み込み符号化を施す。変調部１２は、入力された符号データを変調し、変調後のデータである複数の変調データをＤ／Ａ変換部１３に出力する。Ｄ／Ａ変換部１３は、変調部１２が出力した変調データをデジタル信号からアナログ信号に変換する。そして、アナログ信号に変換された変調データは複数のアンテナ１４を介して送信される。

受信側のユーザ端末２０は、複数のアンテナ２４を介して送信側１０のアンテナ１４から送信されたデータを受信する。但し、アンテナ２４が受信したデータは、アンテナ１４から出力された後、空間を伝播する際のノイズの影響を受けていることに留意する必要がある。

アンテナ２４が受信したデータはＡ／Ｄ変換部２３に入力される。Ａ／Ｄ変換部２３は、入力されたデータをアナログ信号からデジタル信号に変換する。Ａ／Ｄ変換部２３は、変換後のデジタル信号を復調部２２に出力する。そして、復調部２２は、Ａ／Ｄ変換部２３が出力したデータを復調する。復調部２２は、復調して得た復調データを復号化部２１に出力する。復号化部２１は、入力された復調データに対して誤り訂正復号化を行う。その結果得られる復号データを使用して、後段の図示せぬＣＰＵなどの処理回路が所定の処理を実施する。

ここで、送信側の基地局１０にある変調データを送信データ、受信側のユーザ端末２０にあるＡ／Ｄ変換部２３が出力したデータを受信データと呼ぶこととする。

図１に示した復調部２２の構成例を図２に示している。復調部２２は、ＲＳ推定回路３１、仮推定回路３２、雑音抑圧回路３３、補間推定回路３４、等化回路３５および復調回路３６により構成されている。

チャネル推定は、復調部２２において実施される。先ず、ＲＳ推定回路３１において、受信データに含まれるパイロット信号であるＲＳ信号と、このＲＳ信号のレプリカとを用いて、ＲＳ信号の位置におけるチャネル推定値を算出する。次に、仮推定回路３２において、ＲＳ信号の位置におけるチャネル推定値から仮推定値を算出する。

続いて、雑音抑圧回路３３において、ＲＳ信号の位置におけるチャネル推定値および仮推定値から雑音の抑圧を行い、補間推定回路３４において、雑音抑圧後の結果を周波数方向および時間方向の二次元補間することにより、ＲＳ信号の位置以外のチャネル推定値を求める。こうして、チャネル推定値が求まると、等化回路３５において、チャネル推定値と、パイロット信号を除いた受信データとを用いて等化処理を行い、更に復調回路３６において等化信号を復調し、復調データを算出する。この復調データは後段の復号化部２１に出力されることになる。

仮推定回路３２の一構成例を図３に示している。この仮推定回路３２は、周波数補間回路５１と忘却処理回路５２とにより構成されている。図４のフローチャートに示すように、先ず、周波数補間回路５１において、ＲＳ信号位置のチャネル推定値から一次仮推定値を算出し（ステップＳ１）、次に忘却処理回路５２において、一次仮推定値に対して過去のチャネル推定値による忘却を行い（ステップＳ２）、仮推定を算出する。

この仮推定値ｈ_tmp の算出方法の例を、図５に示しており、その算出式は次の式（５）である。
ｈ_tmp （０，ｂ，ｎ，０，０）＝（１−α）＊｛Ｗ_f （ｎ，０）＊ｈ_rs（０，ｂ，ｎ，０，０）＋Ｗ_f （ｎ，１）＊ｈ_rs（０，ｂ，ｎ，０，１）｝＋α＊ｈ_rs（０，ｂ，ｎ−１，１，０）……（５）

すなわち、先ず、周波数補間回路５１において、式（５）の右辺第一項の｛｝に示すように、周波数補間による一次仮推定値を求めて、更にその一次仮推定値について、忘却処理回路５２において、過去のパイロット信号による忘却処理を行う。ここで、α（０＜α＜１）は忘却係数である。

ここで、参考までに、図１０および図１１において説明した仮推定をなすための仮推定回路の構成を図６に示している。この場合の仮推定回路３２は、周波数補間回路４１と時間補間回路４２と選択回路４３とにより構成されている。
周波数補間回路４１による周波数方向の補間処理が、式（１）または式（３）で行われる。そして、時間補間回路４２による時間方向の補間処理が、式（２）または式（４）で行われる。しかる後に、選択回路４３において、周波数方向の補間結果か、時間方向の補間結果のいずれかが、適宜選択されるようになっている。

この図６の方式では、前述したように、ＲＳ信号の位置におけるチャネル推定値から求めた周波数方向および時間方向の仮推定値のうち一つを選択して、周波数方向で補間するか、時間方向で補間するかの選択をする必要がある。そのためには、図１０や図１１のＢで示した時間まで、受信データをバッファリングしておくことが必要である。よって、バッファリングのためのメモリ容量が増大し、また、仮推定値の生成もＢで示す時間まで遅延することになる。

これに対して、図３〜図５で示した本発明の一実施の形態では、周波数方向の補間を行った後、過去のパイロット（ＲＳ）信号による忘却処理を行うようにしているので、時間方向での少なくとも１つ次のパイロット（ＲＳ）信号を待つ必要がなくなる。よって、受信データのバッファリングの必要もなく、また処理遅延時間もゼロとすることができる。更には、時間方向の忘却処理を行うので、チャネル推定精度も向上することになる。

図７は本発明の他の実施の形態を説明するための機能ブロックであり、図３と同等部分は同一符号により示している。図７は仮推定回路３２の他の構成例を示すものであり、この仮推定回路３２は、第一の周波数補間回路５１と、忘却処理回路５２と、更には第二の周波数補間回路５３とにより構成されている。

図７の機能ブロックの動作例を、図８に示している。図３〜図５で示した先の実施の形態では、時間方向で一つ前のパイロット信号を用いて忘却処理を行っているが、本実施の形態では、時間方向で二つ前パイロット信号を用いて忘却処理を行うようにしている。

これを数式で示すと、下記の式（６）となる。
ｈ_tmp （０，ｂ，ｎ，０，０）＝（１−α）＊｛Ｗ_f （ｎ，０）＊ｈ_rs（０，ｂ，ｎ，０，０）＋Ｗ_f （ｎ，１）＊ｈ_rs（０，ｂ，ｎ，０，１）｝＋α＊｛Ｗ_f （ｎ−１，０）＊ｈ_rs（０，ｂ，ｎ−１，０，０）＋Ｗ_f （ｎ−１，１）＊ｈ_rs（０，ｂ，ｎ−１，０，１）｝……（６）

すなわち、第一の周波数補間回路５１により、式（６）の右辺第一項の｛｝に示すように、周波数補間による第一の一次仮推定値を求め、更に、第二の周波数補間回路５３により、式（６）の右辺第二項の｛｝に示すように、周波数補間による第二の一次仮推定値を求める。そして、これら第一及び第二の一次仮推定値について、忘却処理回路５２において、過去のパイロット信号による忘却処理を行う。なお、忘却処理は、過去の複数の値を用いて行うこともできることは明白である。

上述した実施の形態においては、送信側を基地局とし、受信側を携帯電話機などのユーザ端末としたが、これに限定されることなく、広く無線通信システムに適用されるものである。また、図２、図３、図４及び図７などの各機能ブロックやその動作は、予めその動作手順をプログラムとして記録媒体に格納しておき、これをコンピュータにより読取らせて実行させるように構成できることは勿論である。

１０ 送信側
１１ 符号化部
１２ 変調部
１３ Ｄ／Ａ変換部
１４，２４ アンテナ
２０ 受信側
２１ 復号化部
２２ 復調部
２３ Ａ／Ｄ変換部
３１ ＲＳ推定回路
３２ 仮推定回路
３３ 雑音抑圧回路
３４ 補間推定回路
３５ 等化回路
３６ 復調回路
５１，５３ 周波数補間回路
５２ 忘却回路

Claims (3)

1. 空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定装置であって、
   前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定手段を有し、
   前記仮推定手段は、
   前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間手段と、
   前記一次推定値に対して、前記一次推定値が算出されたデータと同一の周波数であって前記一次推定値におけるデータを受信する前に受信した基準信号を用いて、前記一次推定値に対して忘却処理を行って前記仮推定値を算出する忘却処理手段とを有し、前記忘却処理後の仮推定値および前記基準信号の位置におけるチャネル推定値から雑音の抑圧を行い、前記雑音の抑圧後の結果を周波数方向および時間方向の二次元補間することにより、前記基準信号の位置以外のデータ位置のチャネル推定をなすことを特徴とするチャネル推定装置。
2. 空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定方法であって、
   前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定ステップを有し、
   前記仮推定ステップは、
   前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間ステップと、
   前記一次推定値に対して前記一次推定値が算出されたデータと同一の周波数であって前記一次推定値におけるデータを受信する前に受信した基準信号を用いて、前記一次推定値に対して忘却処理を行って前記仮推定値を算出する忘却処理ステップとを有し、
   前記忘却処理後の仮推定値および前記基準信号の位置におけるチャネル推定値から雑音の抑圧を行い、前記雑音の抑圧後の結果を周波数方向および時間方向の二次元補間することにより、前記基準信号の位置以外のデータ位置のチャネル推定をなすことを特徴とするチャネル推定方法。
3. 空間伝搬路のチャネル推定用の基準信号が周波数方向及び時間方向の二次元領域に分散配置されて付加された受信データ系列を受信して、前記基準信号に基づいて受信データのチャネル推定をなすチャネル推定方法をコンピュータに実行させるためのプログラムであって、
   前記基準信号の位置におけるチャネル推定値に基づいて、前記基準信号の位置以外のデータ位置のチャネル推定値を算出するための仮推定値を算出する仮推定処理を有し、
   前記仮推定処理は、
   前記基準信号の位置におけるチャネル推定値を用いて周波数方向補間により一次推定値を算出する第一の周波数補間処理と、
   前記一次推定値に対して前記一次推定値が算出されたデータと同一の周波数であって前記一次推定値におけるデータを受信する前に受信した基準信号を用いて、前記一次推定値に対して忘却処理を行って前記仮推定値を算出する忘却処理とを有し、
   前記忘却処理後の仮推定値および前記基準信号の位置におけるチャネル推定値から雑音の抑圧を行い、前記雑音の抑圧後の結果を周波数方向および時間方向の二次元補間することにより、前記基準信号の位置以外のデータ位置のチャネル推定をなすことを特徴とするプログラム。

Images