(57)【要約】  複数の実施形態の１つによる信号処理の方法は、受信信号の決定論的成分を推定することを含む。この推定は、伝送路の推定応答に基づく。推定された決定論的成分に基づいて、受信信号の非決定論的成分が推定される。推定された非決定論的成分のなまった部分に基づいて、雑音推定が得られて、受信信号が、雑音推定に基づいて補正される。別の実施形態による方法は、なまった位置における受信サンプルを算出されたモデルからの値に置換することを含む。
2009505511.tif 000002
【特許請求の範囲】
【請求項１】
伝送路を介して受信されるマルチキャリア信号に基づいて、前記伝送路の応答を推定するステップと、
受信された前記マルチキャリア信号により伝えられる複数のデータ値を推定するステップであって、前記推定するステップが、
（Ａ）前記伝送路の前記推定応答、および、（Ｂ）受信された前記マルチキャリア信号により伝えられる既知のパターンのデータ値に基づいて基準信号を算出するステップと、
前記基準信号と前記複数のデータ値を伝える受信された前記マルチキャリア信号の一部分との時間領域の差異に基づいて雑音推定を算出するステップと、
周波数領域において前記雑音推定に従い前記複数のデータ値を伝える受信された前記マルチキャリア信号の前記一部分を補正するステップと、
を含むステップと
を含むことを特徴とする信号処理の方法。
【請求項２】
前記基準信号を算出するステップが、（Ａ）前記伝送路の前記推定応答と、（Ｂ）前記受信されたマルチキャリア信号のパイロットキャリアにより伝えられる既知のパターンのデータ値とに基づくことを特徴とする請求項１に記載の信号処理の方法。
【請求項３】
前記雑音推定を算出するステップが、時間領域の差異のなまった値を特定するステップを含むことを特徴とする請求項１に記載の信号処理の方法。
【請求項４】
前記なまった値を特定するステップが、前記時間領域の差異の各値に対してしきい値との比較を実行するステップを含むことを特徴とする請求項３に記載の信号処理の方法。
【請求項５】
前記なまった値を特定するステップが、前記時間領域の差異の値から複数の測定値を得るステップと、前記複数の測定値のそれぞれをしきい値と比較するステップとを含むことを特徴とする請求項３に記載の信号処理の方法。
【請求項６】
前記補正するステップが、前記複数のデータ値を伝える前記受信されたマルチキャリア信号の前記一部分から前記雑音推定を取り去るステップを含むことを特徴とする請求項１に記載の記号処理の方法。
【請求項７】
伝送路を介したマルチキャリア信号により伝えられるシンボルを推定するステップと、
前記推定シンボル、基準値の集合、および、前記伝送路の推定応答に基づいて、前記シンボルを伝える前記受信されたマルチキャリア信号の一部分のモデルの算出するステップと、
時間領域において、前記推定シンボルのなまった値を特定するステップと、
前記なまった値の位置に基づいて、前記算出されたモデルの値に従って前記受信されたマルチキャリア信号を補正するステップと、
を含むことを特徴とする信号処理の方法。
【請求項８】
前記伝送路の前記推定応答が、前記受信されたマルチキャリア信号に基づくことを特徴とする請求項７に記載の信号処理の方法。
【請求項９】
前記モデルを算出するステップが、前記推定シンボルの各周波数領域値に対して前記基準値の集合から１つを選択するステップを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１０】
前記モデルを算出するステップが、前記推定シンボルの各値に対してコンステレーション点の集合から１つを選択するステップを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１１】
前記モデルを算出するステップが、前記推定シンボルの各値に対して前記コンステレーション点の集合から最も近い点を選択するステップを含むことを特徴とする請求項１０に記載の信号処理の方法。
【請求項１２】
前記なまった値を特定するステップが、前記推定シンボルの各値に対してしきい値との比較を実行するステップを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１３】
前記なまった値を特定するステップが、前記推定シンボルの前記値から複数の測定値を得るステップと、前記複数の測定値のそれぞれをしきい値と比較するステップとを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１４】
前記補正するステップが、前記受信されたマルチキャリア信号を前記算出されたモデルの対応する値と置換するステップを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１５】
前記伝送路の推定応答に従って前記補正された受信マルチキャリア信号を等化するステップを含むことを特徴とする請求項７に記載の信号処理の方法。
【請求項１６】
伝送路の応答を推定するステップと、
前記伝送路を介して受信されたマルチキャリア信号の少なくとも１つのなまった成分を特定するステップと、
各特定されたなまった成分に対して、前記なまった成分に対応する周波数における前記伝送路の前記推定応答を修正するステップと、
前記受信されたマルチキャリア信号により伝えられる複数のデータ値を推定するステップであって、前記推定ステップが、
（Ａ）前記伝送路の前記修正された推定応答、および、（Ｂ）受信された前記マルチキャリア信号により伝えられる既知のパターンのデータ値に基づいて基準信号を算出するステップと、
前記基準信号と前記既知のパターンを伝える受信された前記マルチキャリア信号の一部分との差異に基づいて雑音推定を算出するステップと、
前記雑音推定に従って、前記複数のデータ値を伝える受信された前記マルチキャリア信号の前記一部分を補正するステップと
を含むステップと、
を含むことを特徴とする信号処理の方法。
【請求項１７】
前記伝送路の前記推定応答が、前記受信信号に基づくことを特徴とする請求項１６に記載の信号処理の方法。
【請求項１８】
前記なまった成分に対応する周波数において前記伝送路の推定応答を修正する前記ステップが、前記なまった成分に対応する前記周波数における前記推定応答に対しての置換の値を算出するステップを含むことを特徴とする請求項１６に記載の信号処理の方法。
【請求項１９】
前記置換の値を算出するステップが、別の周波数の前記推定応答の値に基づいて前記置換の値を補間するステップを含むことを特徴とする請求項１８に記載の信号処理の方法。
【請求項２０】
前記基準信号を算出するステップが、前記受信された信号のパイロットキャリアによって伝えられる既知のパターンのデータ値に基づくことを特徴とする請求項１６に記載の信号処理の方法。
【請求項２１】
前記雑音推定を算出するステップが、前記基準信号と複数のパイロットキャリアを含む前記受信信号の一部分との差異に基づくことを特徴とする請求項１６に記載の信号処理の方法。
【請求項２２】
前記雑音推定を算出するステップが、前記パイロットキャリアの周波数以外の周波数における前記雑音推定の値を補間するステップを含むことを特徴とする請求項２１に記載の信号処理の方法。
【請求項２３】
前記補正するステップが、前記複数のデータ値を伝える前記受信されたマルチキャリア信号の前記一部分から前記雑音推定を取り去るステップを含むことを特徴とする請求項１６に記載の信号処理の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、無線通信に関する。
【背景技術】
【０００２】
マルチキャリア技術は、地上での使用から衛星通信まで、さらには、ホームネットワークから広範囲のテレビジョン放送まで、通信の多くの異なる分野で適用または試験されている。この技術の潜在的な利点には、処理能力の向上およびより優れた帯域幅利用が含まれる。マルチキャリア技術−直交周波数分割多重（ＯＦＤＭ（ｏｒｔｈｏｇｏｎａｌ  ｆｒｅｑｕｅｎｃｙ  ｄｉｖｉｓｉｏｎ  ｍｕｌｔｉｐｌｅｘｉｎｇ））−のスペクトル的に効率的な１つ形態は、デジタル音声放送（ＤＡＢ（ｄｉｇｉｔａｌ  ａｕｄｉｏ  ｂｒｏａｄｃａｓｔｉｎｇ））、デジタル地上テレビジョン放送（ＤＶＢ−ＴまたはＩＳＤＢ−Ｔ（digital terrestrial television broadcasting））、ハンドヘルドデバイス用のデジタルテレビジョン放送（ＤＶＢ−Ｈ（digital television broadcasting for handheld devices））、地上およびハンドヘルド用のデジタルマルチメディア放送（ＤＭＢ−Ｔ／Ｈ（terrestrial and handheld digital multimedia broadcasting））、ツイストペアチャネル上での高ビットレートデジタルサブスクライバサービス用のデジタルサブスクライバライン（ＤＳＬ（ｄｉｇｉｔａｌ  ｓｕｂｓｃｒｉｂｅｒ  ｌｉｎｅ））、ならびに、ＩＥＥＥ８０２．１１および８０２．１６などの規格下でのワイヤレスネットワークなどの異なる用途で既に一般に使用されている。
【０００３】
マルチキャリアのシステム設計は、複数の異なる構成から選択されることができる。例えば、単入力単出力（ＳＩＳＯ（ａ  ｓｉｎｇｌｅ−ｉｎｐｕｔ，ｓｉｎｇｌｅ−ｏｕｔｐｕｔ））システムは、１つのマルチキャリア信号を伝送し、一方、多重入力多重出力（ＭＩＭＯ（ｍｕｌｔｉｐｌｅ−ｉｎｐｕｔ，ｍｕｌｔｉｐｌｅ−ｏｕｔｐｕｔ））システムは、同一の帯域で複数のマルチキャリア信号を伝送し、送信機および受信機においてアンテナアレイを含む。適用される変調方式も、多くの異なる設計オプションから選択されることができる。典型的には、横軸ＰＳＫ（ＱＰＳＫ）、横軸振幅変調（ＱＡＭ（ｑｕａｄｒａｔｕｒｅ  ａｍｐｌｉｔｕｄｅ  ｍｏｄｕｌａｔｉｏｎ））、または、オフセットＱＡＭなどの直角位相相偏移キーイング（ＰＳＫ（ｐｈａｓｅ−ｓｈｉｆｔ  ｋｅｙｉｎｇ））方式が使用される。典型的なＱＡＭ方式には、４−、１６−、６４−、２５６−および１０２４−ＱＡＭが含まれる。
【０００４】
ＯＦＤＭの適用の１つの例である地上デジタル映像放送すなわち「ＤＶＢ−Ｔ」は、欧州電気通信標準化機構の非特許文献１に規定されている。ＤＶＢ−Ｔ規格は、５、６、７または８ＭＨｚの帯域幅で伝送されるマルチキャリア信号を用いた、２Ｋモード（約４．５ｋＨｚ間隔の１７０５本のサブキャリアを用いる）および８Ｋモード（約１．１ｋＨｚ間隔の６８１７本のサブキャリアを用いる）を規定する。また、関連するハンドヘルド規格（ＤＶＢ−Ｈ）は、４ｋモード（３４０９本のサブキャリアを用いる）を規格する。ＤＶＢシステムで使用される変調方式は、ＱＰＳＫ、１６−ＱＡＭおよび／または６４−ＱＡＭを含んでよい。
【０００５】
ＯＦＤＭシステムは、一般に、シングルキャリアシステムよりインパルス雑音に耐性がある。この耐性は、より長い期間のＯＦＤＭシンボルに加えてサブキャリア間でのインパルス雑音のエネルギーの分布に起因する。一般に、インパルス雑音に対するロバストネスは、サブキャリアの数と共に増大する。また、より小さなコンステレーションおよびより強いコードにより、インパルス雑音に対する耐性が増大する。
【０００６】
【非特許文献１】ＥＮ  ３００  ７４４
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、近くの自動車のイグニションシステムにより、ヒーター、冷蔵庫、蛍光灯または調光器などの電気器具の切り替えにより、あるいは、ヘアドライヤまたは電気掃除機などの電気器具の作動により引き起こされる干渉が、ＯＦＤＭ信号を介して受信されるテレビジョン画像をフリーズさせるのに十分な干渉の原因となることがある。他の電子装置からのおよび／または受信機自体からの放射線も、１つまたは複数のキャリアの信号の受信に干渉することがある。１つまたは複数のこのような原因による干渉を軽減することが望ましい。
【課題を解決するための手段】
【０００８】
一実施形態による信号処理の方法は、伝送路を介して受信されるマルチキャリア信号に基づく伝送路の応答の推定と、受信されたマルチキャリア信号により伝えられた複数のデータ値の推定とを含む。複数のデータ値の推定は、基準信号の計算、雑音推定の計算、および、複数のデータ値を伝える受信されたマルチキャリア信号の一部分の補正を含む。基準信号の計算は、（Ａ）伝送路の推定応答、および、（Ｂ）受信されたマルチキャリア信号により伝えられる既知のパターンのデータ値に基づく。雑音推定の計算は、基準信号と複数のデータ値を伝える受信されたマルチキャリア信号の一部分との間の時間領域の差異に基づく。受信されたマルチキャリア信号の一部分の補正は、周波数領域において雑音推定に従って実行される。
【０００９】
別の実施形態による信号処理の方法は、伝送路を介して受信されるマルチキャリア信号による伝えられるシンボルの推定と、シンボルを伝える受信されたマルチキャリア信号の一部分のモデルの計算を含む。モデルの計算は、推定されるシンボル、基準値の集合、および、伝送路の推定応答に基づく。この方法は、また、推定されたシンボルのなまった値の時間領域においての特定と、なまった値の位置に基づいた、計算されたモデルの値に従う受信されたマルチキャリア信号の補正とを含む。
【００１０】
別の実施形態による信号処理の方法は、伝送路の応答の推定と、伝送路を介して受信されたマルチキャリア信号の少なくとも１つのなまった成分の特定とを含む。この方法は、また、それぞれの特定されたなまった成分に対しての、なまった成分の対応する周波数における伝送路の推定応答の修正を含む。この方法は、また、受信されたマルチキャリア信号により伝えられる複数のデータ値の推定を含む。複数のデータ値の推定は、（Ａ）伝送路の修正された推定応答と、（Ｂ）受信されたマルチキャリア信号により伝えられる既知のパターンのデータ値とに基づく基準信号の計算を含む。複数のデータ値の推定は、また、基準信号と既知のパターンを伝える受信されたマルチキャリア信号の一部分との差異に基づいた雑音推定の計算と、雑音推定に従う、複数のデータ値を伝える受信されたマルチキャリア信号の一部分の補正とを含む。
【００１１】
別の実施形態は、伝送路の推定応答が、受信された信号のそれぞれの複数の異なるキャリア周波数における伝送路の推定応答を含む；（Ａ）複数のキャリア周波数の１つにおける伝送路の推定応答と（Ｂ）複数のキャリア周波数の別の１つにおける伝送路の推定応答との比が経時的に変動する；複数のデータ値のそれぞれが、信号の異なるキャリア上で受信される；この方法が、補正された信号に基づいた伝送路の応答の再推定を含む；および／または、この方法が、再推定されたチャネル応答に基づいた雑音推定またはモデルの再計算を含む、上記のいずれかに従う方法を含む。別の実施形態は、上記のいずれかに従う１つまたは複数の方法の機械可読の命令を有するデータ記憶媒体を含む。
【００１２】
一実施形態による装置は、受信されたマルチキャリア信号により伝えられる複数のデータ値を推定するように構成される。装置は、（Ａ）伝送路の推定応答と（Ｂ）受信されたマルチキャリア信号により伝えられる既知のパターンのデータ値とに基づいて基準信号を計算するように構成される基準信号発生器を含む。装置は、また、基準信号と複数のデータ値を伝える受信されたマルチキャリア信号の一部分との間の時間領域での差異に基づいた雑音推定を計算するように構成される雑音推定器を含む。装置は、また、周波数領域において雑音推定に従って、複数のデータ値を伝える受信されたマルチキャリア信号の一部分を補正するように構成される信号補正器を含む。基準信号発生器は、受信されたマルチキャリア信号に基づいて伝送路の推定応答を計算するように構成されるチャネル応答推定器を含む。
【００１３】
別の実施形態による信号処理のための装置は、伝送路を介して受信されたマルチキャリア信号により伝えられるシンボルを推定するように構成されるシンボル推定器を含む。装置は、また、シンボルを伝える受信されるマルチキャリア信号の一部分のモデルを計算するように構成されるモデル計算機を含む。モデル計算機は、推定されるシンボル、基準値の集合、および、伝送路の推定応答に基づいてモデルを計算するように構成される。装置は、また、時間領域において、推定されるシンボルのなまった値を特定するように構成される雑音検出器と；なまった値の位置に基づき、計算されたモデルの値に従って、受信されたマルチキャリア信号を補正するように構成される信号補正器とを含む。
【００１４】
別の実施形態による装置は、伝送路を介して受信されるマルチキャリア信号により伝えられる複数のデータ値を推定するように構成される。装置は、受信されたマルチキャリア信号の少なくとも１つのなまった成分を特定するように構成される基準信号発生器を含む。基準信号発生器は、また、特定されたなまった値のそれぞれに対して、なまった成分に対応する周波数における伝送路の推定応答を修正するように構成される。基準信号発生器は、また、（Ａ）伝送路の修正された推定応答と（Ｂ）受信されたマルチキャリア信号により伝えられる既知のパターンのデータ値とに基づいて基準信号を計算するように構成される。装置は、また、基準信号と既知のパターンを伝える受信されたマルチキャリア信号の一部分との差異に基づいて雑音推定を計算するように構成される雑音推定器と；雑音推定に従って、複数のデータ値を伝える受信されたマルチキャリア信号の一部分を補正するように構成される信号補正器とを含む。
【発明を実施するための最良の形態】
【００１５】
以下の記述では、記号ｊは、２つの異なる形で使用される。いくつかの例では、記号ｊは、虚数の平方根−１を示す（ｅjφ(i)に見られるように）。別の例では、記号ｊは、行列の列などの指数を示すのに使用される（Ｗijに見られるように）。両方の使用例とも、当技術分野で一般的であり、当業者は、記号ｊが現れるそれぞれの例の文脈から、これらの２つのうちのどちらが意図されているかを理解するであろう。
【００１６】
以下の記述では、下線はベクトルを示すのに使用され（ｈに見られるように）、ボールド体は、行列を示すのに使用される（Ｗに見られるように）。時間領域のベクトルおよび行列は小文字で表示され（ｈに見られるように）、周波数領域のベクトルおよび行列は大文字で表示される（Ｗに見られるように）。
【００１７】
その文脈により明白に限定されない限り、「得る」という用語は、計算、算出、測定、推定、受信（例えば、外部装置から）、および、検索（例えば、記憶素子から）を含む任意のその通常の意味を示すのに使用される。その文脈により明白に限定されない限り、「算出する」という用語は、その通常の意味に加えて、計算または他の算出の結果に従って値を選択する（例えば、リストまたは表から）ことを含み、「算出された値」という用語は、その通常の意味に加えて、計算または他の算出の結果に従って選択された（例えば、リストまたは表から）値を含む。
【００１８】
図１に、代表的なＯＦＤＭシステムのブロック図を示す。送信機は、それぞれが対応する１つのＭ個のサブキャリア上に変調されたデータ値を表すＭ個の複素数値を有するシンボルＸを受信する。データ値は、テレビジョンまたは他の高帯域幅信号などの１つの直列データストリームから導出されることができる。別法として、データ値は、多重接続方式などでの、２つ以上（最大Ｍ）の異なるデータストリームから導出されることができる。ユーザトラフィックに対応するデータ値に加えて、データ値は、１つまたは複数のパイロット信号も含むことができ、このパイロット信号の値は、（例えば、既知の擬似ランダム系列に基づいた）時間および／または周波数の確定関数である。シンボルＸは、また、データ値間に点在するならびに／またはデータ値から離れて集まるナル値を含むことができる（例えば、周波数領域の最低端および／または最高端のキャリアに対応する値の）。
【００１９】
サブキャリアは、ｍ−ａｒｙＰＳＫ（例えば、ＱＰＳＫ）またはｍ−ａｒｙＱＡＭ（例えば、１６−ＱＡＭまたは６４−ＱＡＭ）などの方式に従って変調されることができる。変調方式は、サブキャリアごとに変化してよく、このような変化により、ある制約条件下で容量を増加させ、またはビット誤り率を減少させることができる。誤り検出、エラー訂正、ならびに／あるいは、エンコード、パンクチュアリングおよび／またはインターリビングなどの過剰コード化処理は、変調の前にデータ値上で前もって実行されてよい。
【００２０】
図２Ａに示すように、送信機は、逆高速フーリエ変換すなわちＦＦＴ（少なくと、Ｍ個の点において）をシンボルＸに適用し、対応する複素時間領域シークエンスｘを得る。送信機は、シンボル間にガードインターバルを挿入して、シンボル間干渉を減少させることができる。典型的には、送信機は、シークエンスｘの巡回拡大によりこのインターバルを挿入する。例えば、送信機は、シークエンスｘの最後のＮビットのコピーをそのシークエンスの最初の部分にサイクリックプレフィクス（ｃｙｃｌｉｃ  ｐｒｅｆｉｘ）として付加することができ、サイクリックプレフィクスの長さＮは、チャネルインパルス応答の予想長さＬと少なくとも同じ長さになるように選択される。ＤＶＢ送信機は、シークエンスの長さの１／３２、１／１６、１／８または１／４のガードインターバルを適用するように構成されることができる。ガードインターバルが挿入される前あるいは後のどちらかで、時間領域シークエンスは、伝送のために連続化される。２つの以上のこのような送信機が、ＭＩＭＯシステム内に存在してよい。
【００２１】
Ｍ個のサブキャリアのいくつかは、パイロット信号を伝えるためにリザーブされてよい。典型的には、パイロット信号は、受信機において既知のアプリオリであってよいデータ値のパターンを形成する。パイロットキャリアの割り当ては、周波数領域の全体に配分されてよく、チャネル応答が、伝送の帯域幅の全体において調べられるようにする。キャリアの割り当て方式は、固定されるあるいは経時的に変化してよい。ＤＶＢ規格は、各シンボル間のパイロット信号を用いてキャリアが変調される連続パイロットと、いくつかのシンボル間のパイロット信号を用いてならびに別のシンボル間のトラフィック信号を用いてキャリアが変調される分散パイロットを含む方式を規定する。図３に、このような割り当て方式の一例を示しており、白丸がトラフィックチャネルを示しており、黒丸がパイロットチャネルを示している。この実施例では、分散パイロットは、周波数において１２キャリアの間隔で発生し、３キャリアごとにあるシンボルから次のシンボルへずれており時間において４シンボルの期間を有するパターンに従う。
【００２２】
伝送路の効果は、ランダムノイズ（例えば、付加白色ガウス雑音）およびパルチパスフェージング（ｍｕｌｔｉｐａｔｈ  ｆａｄｉｎｇ）を含む。チャネル効果は、時間の関数ｈおよび／または周波数の関数Ｈとしてモデル化される。受信機では、信号が並行処理され、巡回拡大が除去される。得られる時間領域信号ｙは、図２Ｂに示すような受信される周波数領域シンボルＹに（例えば、ＦＦＴを使用して）変換される。シンボルＹが等化されて、伝送されるシンボルＸの推定が得られる。さらなる処理は、復調と復号化を含んでよい。
【００２３】
シングルキャリアシステムと比較して、ＯＦＤＭシステムなどのマルチキャリアシステムは、符号間干渉（ＩＳＩ）の減少およびマルチパス効果に対するロバストネスの向上を示すことができる。経時的な連続シンボルの直交性は、巡回拡大などのガードインターバル（例えば、プレフィクス）を挿入することにより維持される。ＯＦＤＭシステムは、また、それぞれのナローサブチャネルの全体で比較的均一の周波数応答を有する傾向がある。
【００２４】
低移動度の携帯装置から車載用の車両用機器までの範囲の移動性の用途においてマルチキャリアまたはＯＦＤＭ信号の受信を補助することが望ましい場合がある。例えば、可動の環境下でＤＶＢ信号の受信を補助することが望ましい場合がある。チャネル応答は、チャネル推定がより重要ならびに困難にもなるように、移動性の用途において経時的に変化することが期待される場合がある。
【００２５】
インパルス雑音がないとき、１つのシンボル区間を介して受信されるマルチキャリア信号は、伝送されるシンボル、チャネル応答およびランダム雑音プロセスの関数として時間領域内でモデル化されることができる。
【００２６】
【数１】
2009505511.tif 000003
【００２７】
上式で、ｙ（ｋ）は時間ｉで受信される信号を示し、ｈ（ｌ）は遅延ｌでのチャネルインパルス応答を示し、ｘ（ｋ）は、時間ｋでの時間領域伝送信号を示し、ｎ（ｋ）は時間ｋでのランダムノイズプロセスの値を示し、Ｍはサブバンドの数を示し、Ｌはチャネルインパルス応答の長さを示す。ランダムノイズプロセスｎは、ゼロ平均および分散σ2のガウス統計を有すると想定されてよい。
【００２８】
周波数領域では、式（１．１）のモデルは、
Ｙ＝ＨＸ＋Ｎ          （１．２）
として表されることができ、上式では、Ｙはｙの周波数変換を示し；列ベクトルＸはｘの周波数変換を示し；列ベクトルＮはｎの周波数変換を示し；Ｈは、対角がベクトルＷLｈである対角行列であり、（Ｍ×Ｌ）行列ＷLは、周波数変換行列である。使用される周波数変換は、典型的には、離散フーリエ変換（ＤＦＴ）であり、この離散フーリエ変換は、種々の高速フーリエ変換（ＦＦＴ）方法の任意のものを使用して実行されることができる。
【００２９】
一例では、行列ＷLは、（Ｍ×Ｍ）ＤＦＴ行列Ｗの最初のＬ列であり、その要素は以下のように定義される。
【００３０】
【数２】
2009505511.tif 000004
【００３１】
上式で、指数ｉおよびｊ、０＜ｉ，ｊ＜Ｍ−１、は、それぞれ行および列を示す。変換処理は、
【００３２】
【数３】
2009505511.tif 000005
【００３３】
のスケーリングファクタを含んでよい。別法として、この係数は、代わりに逆変換処理で適用されてよく、あるいは、２つの処理に分割されてもよい（例えば、変換処理および逆変換処理のそれぞれが
【００３４】
【数４】
2009505511.tif 000006
【００３５】
によって縮小される）。ベクトルＷLｈは、また、Ｍ個の要素の列ベクトルＨとして表されることができる。（文脈から理解できるように、記号Ｈは、本明細書では、Ｍ個未満の要素のチャネル推定値ベクトルを示すのにも使用される。）  
インパルス雑音の存在下では、１つのシンボル区間を介して受信されるマルチキャリア信号は、以下の式に従って時間領域内でモデル化されることができる。
【００３６】
【数５】
2009505511.tif 000007
【００３７】
上式では、ｙ(0)（ｋ）は、時間ｋで受信される信号を示し、ｉ（ｋ）は時間ｋでのインパルス雑音事象の状態を示す。
【００３８】
周波数領域では、式（２．１）のモデルは、
Ｙ(0)＝ＨＸ＋Ｉ＋Ｎ          （２．２）
として表されることができ、上式では、列ベクトルＹ(0)はｙ(0)の周波数変換；列ベクトルＩはｉの周波数変換を示す。
【００３９】
インパルス雑音事象は、ポアソン分布に従って経時的に発生することが想定される。各インパルス雑音事象は、シンボル区間に比べて短い期間を有し、各事象のエネルギーが、シンボル内のいくつかの、多くの、あるいは、全てのキャリア上に分配されるようにする。受信される信号内でのインパルス雑音の影響を軽減することが望まれる場合がある。
【００４０】
図４Ａに、一実施形態による方法Ｍ１００のためのフローチャートを示す。タスクＴ１００は、伝送路の推定応答を得る。既知の値（パイロット信号など）が、サブキャリアのいくつかのセットＰを介して伝送される用途において、Ｘのこれらの要素の情報は、少なくともそれらのサブキャリア上の推定されるチャネル応答
【００４１】
【数６】
2009505511.tif 000008
【００４２】
を得るのに使用されることができる。式（２．１）または（２．２）のモデル内の種々のパラメータのうち、受信信号は既知であり、推定されるチャネル応答も少なくとも１つの受信シンボルに基づくことができるようにする。一実装形態では、タスクＴ１００は、
【００４３】
【数７】
2009505511.tif 000009
【００４４】
に従って伝送路の応答の推定
【００４５】
【数８】
2009505511.tif 000010
【００４６】
を得る。ここで、ＸPは既知である伝送されるシンボルのこれらの要素を含む列ベクトルである。チャネル推定法のための他の方法および構造は、本明細書に記載されており、タスクＴ１００は、任意のこのような方法を実施するために実行されてよく、あるいはそうでない場合、任意のこのような方法により算出される推定チャネル応答を得る（例えば、他のストレージまたは他の装置から）ために実行されることができる。
【００４７】
上述したように、パイロット信号は、時間および／または周波数におけるサブキャリアの範囲の全体に亘って配分されることができる。いくつかの用途では、ヌルキャリアが、パイロット信号と共にあるいはパイロット信号の代わりに使用されることができる。この記述では、既知の伝送値に関連するパイロット信号の例を主として使用するが、他のデータ構造（パケットヘッダなど）または時間ベースの活動（同期系列など）に起因する伝送値の事前情報も、チャネル応答の推定に適用されることができる。
【００４８】
タスクＴ２００は、基準信号を算出する。基準信号は、伝送される信号Ｘの既知の要素などの、受信信号により伝えられるデータの既知の部分に基づくことができる。基準信号は、また、推定チャネル応答などの、受信信号から得られる情報に基づいて計算されることができる。タスクＴ２００の少なくともいくつかの実装形態では、基準信号が、受信信号の決定論的成分を推定する。例えば、タスクＴ２００は、伝送路を介して受信される伝送信号の既知の位置を推定することができる。
【００４９】
タスクＴ２００は、タスクＴ１００で得られた推定チャネル応答に基づいて基準信号を算出するように構成されることができる。一実装形態では、タスクＴ２００は、推定チャネル応答の結果および１つまたは複数のパイロット信号などの伝送シンボルの既知の位置として基準信号を算出する。例えば、タスクＴ２００は、周波数領域の式
【００５０】
【数９】
2009505511.tif 000011
【００５１】
に従って基準信号を算出するために実施されることができる。別の実装形態では、タスクＴ２００は、式
【００５２】
【数１０】
2009505511.tif 000012
【００５３】
に従って、時間領域のコンボリューションとして基準信号を算出する。
【００５４】
タスクＴ３００は、基準信号に基づいて雑音推定を算出する。図４Ｂに、タスクＴ３００の実装形態Ｔ３１０を含む、方法Ｍ１００の実装形態Ｍ１１０のためのフローチャートを示す。雑音推定算出のタスクＴ３１０は、受信信号の少なくとも一部分から基準信号を除去するサブタスクを実行する。基準信号の除去は、受信信号の余り内のインパルス雑音の検出および／または特徴づけを容易にすることができる。
【００５５】
タスクＴ３１０は、受信信号の少なくとも一部分から基準信号を取り去るために実施されることができる。例えば、タスクＴ３１０は、受信されたパイロット信号から基準信号を取り去るために実施されることができる。少なくともいくつかの実装形態では、タスクＴ３００は、受信信号の推定される決定論的成分（例えば、基準信号）と受信される信号の少なくとも一部分との距離に基づいて、受信信号の推定される非決定論的成分を算出する。
【００５６】
タスクＴ４００は、基準信号に基づいて補正信号を算出する。例えば、タスクＴ４００は、雑音推定に基づいて基準信号を補正するために（すなわち、受信信号内の雑音の影響を低下させるために）実施されることができる。図５に、方法Ｍ１００の実装形態Ｍ１２０のためのフローチャートを示す。方法Ｍ１２０は、補正信号を得るために、受信信号から雑音推定を取り去ることにより、雑音推定に従って受信信号を補正するタスクＴ４００の実装形態Ｔ４１０を含む。
【００５７】
タスクＴ４１０は、周波数領域の式Ｙ(1)＝Ｙ(0)−Ｉ(1)に従って補正信号Ｙ(1)を算出するために実施されることができる。タスクＴ４１０は、また、受信信号からインパルス雑音推定を取り去る前に時間領域のインパルス雑音推定ｉ(1)を周波数領域に変換するために実行されることができる。図５Ｂに、タスクＴ３１０を含む方法Ｍ１１０の実装形態Ｍ１３０のためのフローチャートを示す。
【００５８】
図６Ａに、一実施形態による装置１００のブロック図を示しており、これは、情報信号Ｓ１０を受信するようにならびに対応する補正信号Ｓ４０を生成するように構成される。基準信号発生器１１０は、基準信号Ｓ２０を算出するように構成される。例えば、基準信号発生器は、本明細書に記載するタスクＴ２００の実施を行うことができる。雑音推定計算機１２０は、基準信号Ｓ２０に基づいて、雑音推定Ｓ３０を算出するように構成される。例えば、雑音推定計算機１２０は、本明細書に記載するタスクＴ３００の実施を行うことができる。信号補正器１３０は、受信信号Ｓ１０および雑音推定Ｓ３０に基づいて補正信号Ｓ４０を生成するように構成される。例えば、信号補正器１３０は、本明細書に記載されるタスクＴ４００の実施を行うように構成されることができる。
【００５９】
装置１００の種々の素子は、例えば、同一のチップまたは１つのチップセット内の２つ以上のチップに存在する電子デバイスおよび／または光学デバイスとして実装されることができるが、そのような限定のない他の配置も考慮される。装置１００の１つまたは複数の素子は、マイクロプロセッサ、埋め込みプロセッサ、ＩＰコア、デジタル信号プロセッサ、ＦＰＧＡ（フィールドプログラマブルゲートアレイ）（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙｓ））、ＡＳＳＰ（アプリケーションスペシフィックスタンダードプロダクト）（ａｐｐｌｉｃａｔｉｏｎ−ｓｐｅｃｉｆｉｃ  ｓｔａｎｄａｒｄ  ｐｒｏｄｕｃｔｓ））、ＡＳＩＣ（アプリケーションスペシフィックインテグレーテッドサーキット）（ａｐｐｌｉｃａｔｉｏｎ−ｓｐｅｃｉｆｉｃ  ｉｎｔｅｇｒａｔｅｄ  ｃｉｒｃｕｉｔｓ））などの、論理素子（例えば、トランジスタ、ゲート）の固定アレイまたはプログラム可能な１つまたは複数のアレイ上で実行される１つまたは複数の集合の命令として、全体としてあるいは部分的に、実装されることができる。また、１つまたは複数のこのような素子は、共通の構造を有することができる（例えば、異なる時間に異なる素子に対応するコードの一部を実行するのに使用されるプロセッサ、または、異なる時間に異なる素子の代わりに処理を実行する電子および／または光学デバイスの配置）。
【００６０】
図６Ｂに、装置１００の実装形態１０２のブロック図を示す。装置１０２では、信号補正器１３０が、受信信号Ｓ１０から雑音推定Ｓ３０を取り去るように構成される、可算器などのコンバイナ１３２として実装される。
【００６１】
図７Ａに、方法Ｍ１１０の実装形態Ｍ２００のためのフローチャートを示す。方法Ｍ２００は、受信信号の推定される非決定論的成分に基づいて雑音推定を算出するタスクＴ３００の実装形態Ｔ３２０を含む。
【００６２】
タスクＴ３２０は、受信信号から基準信号を除去するサブタスクを含む。例えば、タスクＴ３２０は、時間領域において受信信号と基準信号との差異を算出するように構成される。１つのこのような実装形態において、タスクＴ３２０は、式
【００６３】
【数１１】
2009505511.tif 000013
【００６４】
ここで、
【００６５】
【数１２】
2009505511.tif 000014
【００６６】
は受信信号の推定される非決定論的成分を表す、に従って受信信号から基準信号を取り去る。この場合、タスクＴ３２０は、受信信号ｙ(0)から基準信号を取り去る前に、周波数領域の基準信号
【００６７】
【数１３】
2009505511.tif 000015
【００６８】
を時間領域に変換することができる。  
タスクＴ３２０は、また、推定される非決定論的成分の１つまたは複数のなまった値を特定するサブタスクを含む。このサブタスクは、推定される非決定論的成分の値をしきい値Ｓ６０と比較することを含んでよい。しきい値Ｓ６０との比較は、推定される非決定論的成分の各値に対して実行されることができる。一例において、しきい値Ｓ６０は、推定される非決定論的成分
【００６９】
【数１４】
2009505511.tif 000016
【００７０】
の各時間領域サンプルの出力（例えば、大きさの２乗）と比較される。しきい値Ｓ６０より大きい値（あるいは、別法として、しきい値Ｓ６０以上）は、なまっていると特定される。別の実施例では、この比較は、出力以外（値の対数など）の時間領域のそれぞれの値の測定により実行されることができる。
【００７１】
別法として、しきい値Ｓ６０は、ある時間枠での値の大きさの２乗の合計などの、ある時間枠内の１組の値に関連する測定値と比較されることができる。時間枠の幅（これは、シンボルの幅より大きい、と等しい、あるいは、より小さい）は、インパルス雑音事象の推定期間に従って選択されることができる。１つのこのような実施例では、しきい値Ｓ６０は、その時間枠の幅に亘っての推定される非決定論的成分の平均信号出力と比較される。別のこのような実地例では、しきい値Ｓ６０は、その時間枠の幅に亘っての推定される非決定論的成分の合計信号出力と比較される。時間枠の測定値がしきい値Ｓ６０より大きい（あるいは、別法として、しきい値Ｓ６０以上の）場合、その時間枠内の値はなまっていると特定される。別法として、値が、なまっていると特定された少なくとも一定数の時間枠内にある場合、その値はなまっていると特定されることができる。値が２つ以上の時間枠内にあるようにするために、隣接する時間枠は重複してよく、あるいは、各値が１つの時間枠内のみにあるようにするために重複しないでよい。
【００７２】
１つまたは複数の信号条件に従ってしきい値Ｓ６０の値を得ることが望ましい場合がある。しきい値Ｓ６０は、受信信号Ｓ１０の出力の測定値から、あるいは、信号Ｓ１０の一部（例えば、受信パイロット信号）の出力の測定値から、得られることができる。例えば、しきい値Ｓ６０は、サンプルごとの平均チャネル出力の推定
【００７３】
【数１５】
2009505511.tif 000017
【００７４】
に基づいてよく、これは、以下のような式に従って算出されることができる：
【００７５】
【数１６】
2009505511.tif 000018
【００７６】
平均は、Ｈ(1)の値が利用できるシンボルの全てのサブキャリアに亘って、あるいは、シンボルのサブキャリアの集合Ｐ（例えば、パイロット信号の集合）に亘って得られることができる。式（３．１）での平均化は、平均値または中央値を使用して算出されることができる。別法として、しきい値Ｓ６０は、別の統計フィルタを使用して算出されるチャネル出力の測定値に基づいてよい。  
【００７７】
しきい値Ｓ６０は、荷重係数を含んでよい。一例では、しきい値Ｓ６０は値
【００７８】
【数１７】
2009505511.tif 000019
【００７９】
を有し、ここで、ｗ1は、高出力の信号のための間違い警報の確率に対する所望の検出率に平衡を保たせるように選択されることができる荷重係数である。一例では、ｗ1の値は２５である。
【００８０】
所定の値の集合の中からの選択に従いしきい値Ｓ６０の値を変化させることが望ましい場合がある。例えば、ｗ1の値は、電流信号特性とモデルのセットの１つとの間の適合に従って選択されることができる。別法として、しきい値Ｓ６０またはその係数は、例えば、本明細書に記載する選択および／または算出に従って得られた値に固定されてもよい。
【００８１】
タスクＴ３２０は、なまっていると特定された推定される非決定論的成分の値から雑音推定ｉ(1)を得るように構成される。雑音として分類される成分の各値に対して、その値の一部または全てが雑音に起因するとみなすことができる。雑音としてのサンプルの分類が正確である場合、サンプル出力の全てが雑音に起因と考えるほうが通常よい。というのは、なまったサンプルの雑音成分は、信号成分より大きい傾向があるからである。この場合、非直線性のブランキング（すなわち、出力ゼロ化）が、雑音推定の算出に使用されることができる。一例では、雑音推定ｉ(1)は、以下の式に従って時間領域で計算される。
【００８２】
【数１８】
2009505511.tif 000020
【００８３】
しかし、雑音としてのサンプルの分類が間違っている場合、サンプル出力の全てが雑音に起因すると考えることは、さらなる信号歪を引き起こす可能性がある。この場合、非直線性のリミッティング（または、クリッピング）が、代わりに、雑音推定の計算に使用されることができる。例えば、信号は、しきい値Ｓ６０を超えるサンプル出力の一部分のみが雑音に起因すると考えることにより、以下のような式に従ってクリッピングされることができる。
【００８４】
【数１９】
2009505511.tif 000021
【００８５】
図７Ｂに、方法Ｍ２００の実装形態Ｍ２１０のためのフローチャートを示す。方法Ｍ２１０は、受信信号から雑音推定を取り去ることにより、雑音推定に従って補正信号を算出するタスクＴ４００の実装形態Ｔ４１０を含む。
【００８６】
図８Ａに、タスクＴ１００のインスタンスに適用されることができる、パイロット信号ベクトルＸPの大きさ対周波数のプロットを示す。図８Ｂに、タスクＴ２００のインスタンスによって算出されることができる、対応する基準信号
【００８７】
【数２０】
2009505511.tif 000022
【００８８】
の大きさ対周波数のプロットを示す。図８Ｃに、タスクＴ３２０のインスタンスで算出されることができる非決定論的成分
【００８９】
【数２１】
2009505511.tif 000023
【００９０】
の大きさの２乗対時間のプロットを示す。また、図８Ｃに、雑音推定ｉ(1)を算出するために成分
【００９１】
【数２２】
2009505511.tif 000024
【００９２】
に適用されるしきい値Ｓ６０を示す。図示を容易にするために、図８Ａおよび８Ｂは、実数値ベクトルのプロットを示すが、実際には、少なくともベクトル
【００９３】
【数２３】
2009505511.tif 000025
【００９４】
は、典型的には、複素数値を有する。
【００９５】
図９に、装置１００の実装形態２００のブロック図を示す。基準信号発生器１１０の実装形態２１０は、周波数領域の値に作用するように構成され、雑音推定計算機１２０の実装形態２２０は、時間領域の値に作用するように構成される。また、図９に示す配置は、装置１０２の一部として構成されることができる変換ブロックＸ１０を含む。別法として、変換ブロックＸ１０は、分離素子またはプロセスとして構成されることができ、装置１００のインスタンスを含むシステムが、たとえ装置１００が作動されていない場合でも、情報信号Ｓ１０の周波数変換を実行できるようにする。変換ブロックＸ１０は、既知のあるいは展開することができるアーキテクチャまたはプロセスに従って、離散フーリエ変換（例えば、高速フーリエ変換）、または、用途に適した他の周波数変換を実行することができる。
【００９６】
基準信号発生器２１０は、基準信号Ｓ２０のインスタンスＳ２２０を算出する（例えば、受信信号Ｓ１０の周波数変換に基づいて）ように構成される。雑音推定計算機２２０は、基準信号Ｓ２２０および受信信号Ｓ１０に基づいて、雑音推定Ｓ３０のインスタンスＳ２３０を算出するように構成される。信号補正器１３０は、雑音推定Ｓ２３０および受信信号Ｓ１０の周波数変換に基づいて、補正信号Ｓ４０のインスタンスＳ２４０を算出するように構成される。
【００９７】
図１０に、装置２００の実装形態２０２のブロック図を示す。パターンＳ５０は、パイロット信号パターンＸPなどの、伝送信号の既知の部分である。基準信号発生器２１２および雑音推定計算機２２２は、それぞれ、基準信号発生器２１０および雑音推定計算機２２０の実装形態である。
【００９８】
基準信号発生器２１２は、伝送路の推定応答を算出する（例えば、受信信号Ｓ１０に基づいて）ように構成されるチャネル応答推定器１０５を含む。コンバイナ２３５は、その入力の両方に基づく出力を生成するように構成される。この実施例では、コンバイナ２３５は、推定チャネル応答にパターンＳ５０を重ね合わせて周波数領域で基準信号を得るように構成される。パターンＳ５０は、記憶され、生成され、ならびに／または、他の装置から受信されることができる。コンバイナ２３５は、基準信号
【００９９】
【数２４】
2009505511.tif 000026
【０１００】
を得るために周波数領域ベクトル
【０１０１】
【数２５】
2009505511.tif 000027
【０１０２】
およびＸPをかけるように構成される乗算器として実装されることができる。
【０１０３】
逆変換ブロックＸ２０は、周波数領域内で基準信号に逆周波数変換を実行して時間領域内で基準信号Ｓ２２２を得るように構成される。逆離散フーリエ変換（例えば、逆高速フーリエ変換）を実行するように構成されることができる逆変換ブロックＸ２０は、既知のあるいは展開することができる任意のアーキテクチャまたはプロセスに従って構築されることができる。逆変換ブロックＸ２０は、また、分離素子またはプロセスとして構成されることができ、装置１００のインスタンスを含むシステムが、たとえ装置１００が起動されていない場合でも、ブロックＸ２０を使用して逆周波数変換を実行できるようにする。基準信号Ｓ２２２は、受信信号の決定論的成分を推定することができる。  
別の実装形態では、基準信号発生器２１０は、推定チャネル応答
【０１０４】
【数２６】
2009505511.tif 000028
【０１０５】
と既知のパターンｘPとをコンボリューションして時間領域内で基準信号
【０１０６】
【数２７】
2009505511.tif 000029
【０１０７】
を得るように構成される。この場合、基準信号発生器２１０は、コンボルバとして実装されるコンバイナを含むことができる。
【０１０８】
雑音推定器２２２は、受信信号Ｓ１０から基準信号Ｓ２２２を取り去るように構成されるコンバイナ２３０を含む。出力計算機２５０は、重ね合わせされた信号の大きさの２乗を計算するように構成される。別の実装形態では、出力計算機２５０は、本明細書に記載するように、重ね合わせされた信号に時間枠を適用するように構成されることができる。比較器２４０は、出力値をしきい値Ｓ６０と比較して時間領域内で雑音推定を得るように構成される。例えば、比較器２４０は、式（３．２）または（３．３）に記載する処理を実行するように構成されることができる。
【０１０９】
変換ブロックＸ３０は、時間領域の雑音推定に周波数変換を実行して周波数領域内に雑音推定Ｓ２３２を得るように構成される。離散フーリエ変換（例えば、高速フーリエ変換）を実行するように構成されることができる変換ブロックＸ３０は、既知のあるいは展開することができる任意のアーキテクチャまたはプロセスに従って構築されることができる。変換ブロックＸ１０およびＸ３０は、異なる時間に異なる信号に作用する同一の構造またはプロセスとして実装されることができる。
【０１１０】
信号補正器１３０は、周波数領域での雑音推定Ｓ２３２と受信信号Ｓ１０との差異として補正信号Ｓ２４２を生成するように構成される。装置２０２の別の実装形態では、信号補正器１３０は、時間領域内で受信信号Ｓ１０から雑音推定を取り去るように構成される。
【０１１１】
上述のように、伝送信号が既知である信号Ｓ１０のサブキャリアに対して、タスクＴ１００は、対応する既知の伝送信号により受信信号を分割することによりチャネル応答の推定を得るように実装されることができる。
【０１１２】
【数２８】
2009505511.tif 000030
【０１１３】
タスクＴ１００は、また、集合Ｐにないサブキャリアに対してチャネル応答推定を得るために、（時間領域および／または周波数領域において）補間処理を含むように構成されることができる。推定Ｈ(0)は、推定チャネル応答Ｈ(1)として使用されることができ、あるいは、Ｈ(0)からＨ(1)を得るためにさらなる処理が実行されることができる。例えば、推定チャネル応答に、特に、チャネル応答の長さがサブキャリアの数よりもかなり少ないような用途に対しては、フィルタリングまたはカーブフィッティング処理を実行することが望ましい場合がある。フィルタリングされたチャネル応答推定を得るタスクＴ１００の実装形態は、本明細書に記載されており、推定チャネル応答は、最小平均２乗（ｌｅａｓｔ−ｍｅａｎ−ｓｑｕａｒｅｓ（ＬＭＳ））アルゴリズムなどの予測プロセスを例えば使用することにより、先の受信シンボルから得られる。
【０１１４】
タスクＴ１００は、初期推定Ｈ(0)に基づいてフィルタリングされたチャネル反応推定Ｈ(1)を得るように実装されることができる。図１１Ａに、このようなタスクを実行するのに使用されることができるチャネル推定器１０５の実装形態１１２のブロック図を示す。チャネル推定器１１２は、キャリアｐ∈Ｐに対して推定Ｈ0（ｐ）からフィルタリングされた推定Ｈ1（ｐ）を生成するように構成されるフィルタ４００を含む。このような処理は、Ｐの他の要素に対してフィルタリングされた推定を得るために連続してならびに／または並行して実行されることができる。
【０１１５】
タスクＴ１００の実装形態Ｔ１０２は、初期チャネル応答推定を異なる領域に変換する、予想されるチャネルインパルス応答長さに従ってそれの端を切り捨てる、ならびに、切られた信号に逆変換を実行することによりチャネル応答推定を再算出する。このようなフィルタリングは、応答が経時的に比較的素早く変化する（例えば、あるシンボルから次のシンボルへ）チャネルに対して特に適する場合がある。一例では、タスクＴ１０２は、その逆離散フーリエ変換を利用することによりＨ(0)を時間領域シークエンスｈ(0)に変換する：
ｈ(0)＝ＩＤＦＴ（Ｈ(0)）＝ＷHＨ(0)          （４．２）
上式では、上付き文字Ｈは、共役転置処理を示す。
【０１１６】
典型的には、伝送シンボルＸのＭ個の値のいくつかのみが、既知のアプリオリである（例えば、パイロット信号の集合Ｐ）。この場合、初期推定ベクトルＨ(0)は、Ｘのそれらの値に対応する要素のみを含むことができ、変換行列Ｗは、これらの要素に対応するＷの列のみを含む減らした行列ＷPにより置き換えられることができる。
【０１１７】
Ｈ(0)に逆変換を適用することにより、Ｍ個の要素のベクトルｈ(0)が生成される。しかし、チャネルインパルス応答は、既知であってよく、あるいは、Ｍより小さい長さＬを有すると予想される場合がある。その結果、ｈ(0)の最初のＬ個の値のみが有効であると想定され、さらに、残りの値が雑音であると想定されることができる。タスクＴ１０２は、ｈ(0)の端を切り捨て（例えば、長さＬに）、Ｈ(0)の領域に戻す切られたベクトルｈ1の変換を適用してチャネル応答推定Ｈ(1)を得る。一実施形態では、ｈ(0)の切り捨ては、Ｗの最初のＬ個の列のみを含む減少された変換行列ＷLを適用することにより実行される：
Ｈ(1)＝ＤＦＴ（ＩＤＦＴ（Ｈ(0)）＝ＷLＷHＨ(0)          （４．３）
一実施形態によるシステムまたは装置は、このような処理を実行するように構成されるフィルタ４００の実装形態を含むことができる。当業者なら、変換行列（例えば、ＷL、ＷH）が、実行時間の前に算出されて記憶されることができることを理解するであろう。
【０１１８】
タスクＴ１００の別の実装形態Ｔ１０４は、１つまたは複数の先行するシンボル区間の間で受信される情報に基づいて、チャネル応答推定Ｈ(1)を得る。このようなフィルタリングは、応答が経時的に比較的ゆっくりと変化する（例えば、あるシンボルから次のシンボルへのランダム雑音プロセス未満）チャネルに対して特に適する場合がある。図１１Ｂに、チャネル推定器１０５の実装形態１１４のブロック図を示す。推定器１１４は、遅延素子、および、キャリアｐ∈Ｐに対してシンボル区間ｌの間でチャネル応答のフィルタリングされた推定
【０１１９】
【数２９】
2009505511.tif 000031
【０１２０】
を得る予測フィルタ４１０を含んでおり、
【０１２１】
【数３０】
2009505511.tif 000032
【０１２２】
は、ｌより前のシンボル区間に対応する１つまたは複数の推定Ｈ(0)（ｐ）（ならびに、場合によっては現行の推定
【０１２３】
【数３１】
2009505511.tif 000033
【０１２４】
）に基づく。
【０１２５】
タスクＴ１０４は、線形予測を実行するように構成されることができる。例えば、タスクＴ１０４は、先行するシンボル区間の間に受信された信号からの加重値の合計に基づいてチャネル推定を得ることができる。いくつかの例では（連続パイロットを伝えるキャリア上など）、加重値は、連続するシンボル区間に対応することができる。別の例では（分散パイロットを伝えるキャリア上など）、加重値は、連続しないシンボル区間（例えば、４つのシンボル区間ごと）に対応することができる。図１２Ａに、先行するシンボル区間の間に受信された信号のＫ個の加重値の合計を算出する予測フィルタ４１０の実装形態４１２を示す。Ｋの値は、固定されてよく、あるいは、チャネル応答の経時的に予想される変化率に基づいて選択されることができる。一実施形態では、Ｋの値は６である。
【０１２６】
いくつかのフィルタ加重は適応型（例えば、固定されるのではなく、定期的にあるいはある事象ごとに更新されるで）であることが望ましい場合がある。例えば、加重の値は、Ｈ(0)（Ｐ）またはＨ(1)（Ｐ）の１つまたは複数の現行値、Ｈ(0)（Ｐ）またはＨ(1)（Ｐ）の１つまたは複数の過去の値、ならびに／または、別のキャリアからのチャネル応答推定の特性に基づいてよい。いくつかの場合、加重の値は、現行の推定
【０１２７】
【数３２】
2009505511.tif 000034
【０１２８】
と
【０１２９】
【数３３】
2009505511.tif 000035
【０１３０】
との差異などの誤差信号に基づいてよい。  
タスクＴ１０４は、最小平均２乗（ＬＭＳ）あるいは以下の式などの最急勾配法の適用に従って１組のフィルタ加重を更新するように構成されることができる：
【０１３１】
【数３４】
2009505511.tif 000036
【０１３２】
上式では、Ｗl（ｐ）はシンボル区間ｌに対する加重値のベクトル［Ｗ0，Ｗ1，．．．，ＷM-1］を示し；
【０１３３】
【数３５】
2009505511.tif 000037
【０１３４】
は加重に対する入力値のベクトルに対応する複素共役を示し；ｅl（ｐ）は、値
【０１３５】
【数３６】
2009505511.tif 000038
【０１３６】
を有する誤差信号を示し；μは、固定され、選択され、または、適合型であってよいステップサイズパラメータを示し；ｒは、連続パイロットキャリアに対しては１を有して分散パイロットキャリアに対しては４を有する。別の実施形態では、タスクＴ１０４は、正規化最小平均２乗または再帰的最小２乗アルゴリズムなどの関連する方法を適用するように構成されることができる。図１２Ｂに、チャネル推定器１１４の実装形態１１６のブロック図を示しており、予測フィルタ４１２の実装形態４１４内での１つまたは複数の加重の値は、このようなアルゴリズムに基づいて更新されることができる。
【０１３７】
経時的に変化するチャネル応答の一部分に予測作用を集中させることが望ましい場合がある。例えば、平均推定からの現行の推定の距離を予測することが望ましい場合がある。図１３Ａに、アベレージャ５１０および予測フィルタ４１０の実装形態４１６を含むチャネル推定器１０５の実装形態１１８のブロック図を示す。アベレージャ５１０は、２つ以上の雑音推定に基づく有限インパルス応答（ｆｉｎｉｔｅ  ｉｍｐｕｌｓｅ  ｒｅｓｐｏｎｓｅ（ＦＩＲ））を有する移動式アベレージャとして実装されることができる。分散パイロットキャリアに対して、アベレージャ５１０は、以下の式に従って平均値を生成することができる：
【０１３８】
【数３７】
2009505511.tif 000039
【０１３９】
上式で、αは、固定され、選択され、または、適応されてよい加重係数（０＜α＜１）である。別法として、アベレージャ５１０は、１つまたは複数の先行状態に基づく無限インパルス応答（ｉｎｆｉｎｉｔｅ  ｉｍｐｕｌｓｅ  ｒｅｓｐｏｎｓｅ（ＩＩＲ））を有する減衰積分器として実装されることができる：
【０１４０】
【数３８】
2009505511.tif 000040
【０１４１】
上式で、βは、固定され、選択され、または、適応されてよい加重係数（０＜β＜１）である。一実施形態では、係数βは、６３／６４の値を有する。
【０１４２】
フィルタ４１６は、平均推定からの現行の雑音推定の距離を示す値Ｇを受信する。図１３Ａの一例では、Ｇl（ｐ）の値は、
【０１４３】
【数３９】
2009505511.tif 000041
【０１４４】
として算出される。分散パイロットキャリアに対して、フィルタ４１６の加重値は、以下の正規化ＬＭＳの式に従って更新されることができる：
【０１４５】
【数４０】
2009505511.tif 000042
【０１４６】
上式で、Ｇl（ｐ）は、キャリアｐおよびシンボル区間ｌに対しての加重に対する入力値のベクトル［Ｇ0，Ｇ1，．．．，ＧM-1］を示し、ｒは、連続パイロットキャリアに対しては１の値を有し、分散パイロットキャリアに対しては４の値を有する。
【０１４７】
いつかの用途では（例えば、非常に雑音のある環境では）、他のキャリア（隣接するキャリアなど）からの情報に基づいて加重の値を更新することが好適である場合がある。１組の（あるいは、全ての）キャリアを介する加重値補正の平均化は、より堅調なチャネル応答推定をもたらすことができる。分散パイロットの一例では、Ｐにおける分散パイロット（ＳＰ）に対しての補正値は平均化され、これらの全てのキャリアに対して同一の加重値が使用される。
【０１４８】
【数４１】
2009505511.tif 000043
【０１４９】
同様の処理が、Ｐにおける連続パイロットキャリアに対して共通の加重値ベクトルを得るために適用されてよい。分散パイロットに対する加重は、この実施例では、シンボル区間ごとに更新されるが、キャリア上に新しいパイロット信号が出現するのに対応して、４つのシンボル区間ごとのみにおいてこれらのキャリアのそれぞれに対して新しいフィルタリングされた推定を算出することが望ましい場合があることに留意されたい。当業者には、タスクＴ１０４も、１つまたは複数の後続のシンボル区間の間（１つまたは複数の先行するシンボル区間の間に受信される情報に加えてあるいはそれの代わりとしてのどちらかで）受信される情報に基づいてチャネル応答推定Ｈ(1)を得るように実装されることができることを理解するであろう。
【０１５０】
本明細書に記載するタスクＴ１００の実装形態は、パイロット信号が伝送される間のシンボル区間に対してのチャンネル応答推定を得るために分散パイロットキャリアに適用されることができる。このようなタスクは、また、分散パイロット信号の間のシンボル区間におけるキャリアに対してのチャネル応答を得るために経時的に補間を実行することができる。
【０１５１】
図１３Ｂに、現行推定の区間に先行する３つのシンボル区間に対しての推定チャネル応答を算出する補間器５２０のブロック図を示す。この実施例では、この線形補間器に対する加重νは、以下のとおりである。  
［ν0,0，ν1,0］＝［０，１］
［ν0,1，ν1,1］＝［１／４，３／４］
［ν0,2，ν1,2］＝［１／２，１／２］
［ν0,3，ν1,3］＝［３／４，１／４］
補間器５２０の他の実装形態は、非線形係数（例えば、２次または３次などの高次多項式、または、他のいくつかの関数または統計）を適用することができ、および／または、３つ以上の推定に基づいて補間を実行することができる。
【０１５２】
タスクＴ１００の一実装形態は、データキャリアに対してのチャネル推定を得るために周波数軸に沿って補間をさらに行うことができる。このような補間の可能性のある形態は、２つの最も近いパイロットチャネル推定の間の線形補完から、そのシンボル区間に対してのより多くの（場合によっては、全ての）利用可能なチャネル推定に基づくＦＩＲまたはＩＩＲ補間にまで及ぶ。中間のパイロットおよび／またはデータキャリアに対してのチャネル応答推定の補間は、時間および周波数領域においての２次元（例えば、サーフェスフィッティング）処理として実施されることもできる。
【０１５３】
時間領域においてのインパルス雑音に対しての推定の算出は、計算上において集中的である場合がある。例えば、時間領域の方法は、各受信信号に対しての１つまたは複数の追加の変換処理を含む場合がある。６４−キャリアシステム（例えば、ワイヤレスＬＡＮ用）では、このような複雑さも許容される場合があるが、このような処理を補助するための追加チップ領域のコストが、数千のキャリア（例えば、ＤＶＢ−Ｔ／Ｈ）を有するシステムに対しては非常に高価になる可能性がある。しかし、いくつかの実装形態では、チャネル応答を推定する、他のタイプの雑音を補正する、あるいは、コンステレーション点の集合に受信信号を写像する手順などの受信信号上での他の手順の間に実行される１つまたは複数の変換（および／または逆変換）処理の結果を再利用するために、本明細書に記載する時間領域の方法を構成することが可能である。同様に、異なる時間に異なる信号上でこのようなタスクを実行するように構成されることができる構造および／またはプロセスを含むシステム内でこのような方法を実行するための装置を構成することが可能である。
【０１５４】
受信信号での周波数領域インパルス雑音（すなわち、「狭帯域干渉」）の影響を軽減することが望ましい場合がある。このような雑音事象のエネルギーは、典型的には、時間領域インパルス雑音事象より狭い帯域幅に亘って配分され、雑音事象は、１つまたは複数のシンボルに亘る１つまたは複数のキャリアを回復不能になまらせる。
【０１５５】
狭帯域干渉は、他の通信伝送のスペクトルが所期の伝送のスペクトルと重複した場合に発生する場合がある。いくつかのデジタルテレビジョンシステムでは、例えば、アナログテレビジョン信号は、同一のスペクトル内で放送され、アナログルマ（ｌｕｍａ）および／またはクロマテレビジョン（ｃｈｒｏｍａ  ｔｅｌｅｖｉｓｉｏｎ）信号（および／または、これらの信号の調波）用のキャリアが、デジタルテレビジョン受信機に対する狭帯域の同一チャネル干渉を示すことができるようにする。
【０１５６】
図１４Ａに、方法Ｍ１００の実装形態Ｍ３００のためのフローチャートを示す。方法Ｍ３００は、推定チャネル応答を得るタスクＴ１００の実装形態Ｔ１１０を含む。タスクＴ１１０は、受信信号のなまったキャリア（例えば、狭帯域干渉によってなまったキャリア）を特定するサブタスクを含む。いくつかの実装形態では、タスクＴ１１０は、受信信号内のなまったキャリアを特定するために、受信信号Ｓ１０にしきい値を適用する。例えば、タスクＴ１１０は、Ｓ１０の各パイロット信号をしきい値と比較することができる。タスクＴ１１０は、また、修正推定Ｈ(1)を生成するため、特定されたキャリアに従って初期推定チャネル応答Ｈ(0)を修正する。
【０１５７】
タスクＴ１１０は、例えば式（４．１）に従って、対応する既知の伝送信号により受信信号Ｓ１０を分割することにより、チャネル応答の初期推定Ｈ(0)を得るように構成されることができる。このような推定は、少なくとも最初の反復に対しては、十分に正確である可能性があるが、任意の方法または本明細書に記載するチャネル応答の推定方法が、最初のまたは後続の反復に使用されることができる。例えば、タスクＴ１１０は、ＩＦＦＴフィルタリングまたはカーブフィッティングなどの追加の処理、あるいは、最小平均２乗プロセスなどの先行するシンボルに基づく予測プロセスを使用して初期推定を得るように構成されることができる。
【０１５８】
タスクＴ１１０は、受信信号Ｓ１０の各パイロット値出力（例えば、大きさの２乗）としきい値Ｓ７０とを比較することにより、なまったパイロット信号を特定するように構成されることができ、しきい値Ｓ７０を超える（あるいは、別法として、しきい値Ｓ７０以上の）出力を有するキャリアがなまっていると特定されるようにする。別法として、しきい値Ｓ７０は、時間枠内の値の大きさの２乗の合計などの、ある時間枠内の１組の値に関連する測定値と比較されることができる。周波数におけるこの時間枠の幅は、インパルス雑音事象の推定される帯域幅に従って選択されることができる。そのような一例では、しきい値Ｓ７０は、時間枠の幅の全体に亘る受信信号Ｓ１０の平均パイロット信号出力と比較される。別のそのような実施例では、しきい値Ｓ７０は、時間枠の幅の全体に亘る受信信号Ｓ１０の合計パイロット信号出力と比較される。時間枠の測定値がしきい値Ｓ７０より大きい（あるいは、別法として、しきい値Ｓ７０以上である）場合、時間枠内の値はなまっていると特定される。別法として、値は、なまっていると特定される少なくとも一定数の時間枠内にある場合に、なまっていると特定されることができる。値が２つ以上の時間枠内にあるようにするために、隣接する時間枠は周波数において重複してよく、あるいは、各値が１つの時間枠内のみにあるようにするために重複しないでよい。
【０１５９】
１つまたは複数の信号条件に従ってしきい値Ｓ７０の値を得ることが望ましい場合がある。しきい値Ｓ７０は、受信信号Ｓ１０の出力の測定値から、あるいは、信号Ｓ１０の一部（例えば、受信パイロット信号）の出力の測定値から、得られることができる。例えば、しきい値Ｓ７０は、サンプルごとの平均チャネル出力の推定
【０１６０】
【数４２】
2009505511.tif 000044
【０１６１】
に基づいてよく、これは、以下の式などの式に従って算出されることができる。
【０１６２】
【数４３】
2009505511.tif 000045
【０１６３】
平均は、Ｈ(0)の値が利用できるシンボルの全てのサブキャリアに亘って、あるいは、シンボルのサブキャリアの集合Ｐ（例えば、パイロット信号の集合）に亘って得られることができる。式（５．１）での平均化は、平均値または中央値を使用して算出されることができる。別法として、しきい値Ｓ７０は、別の統計フィルタを使用して算出されるチャネル出力の測定値に基づいてよい。  
しきい値Ｓ７０は、荷重係数を含んでよい。一例では、しきい値Ｓ７０は値
【０１６４】
【数４４】
2009505511.tif 000046
【０１６５】
を有し、ここで、ｗ2は、高出力の信号のための間違い警報の確率に対する所望の検出率に平衡を保たせるように選択されることができる荷重係数である。一例では、ｗ2の値は９である。
【０１６６】
１組の所定の値の中からの選択に従ってしきい値Ｓ７０の値を変化させることが望ましい場合がある。例えば、ｗ2の値は、電流信号特性とモデルのセットの１つとの間の適合に従って選択されることができる。別法として、しきい値Ｓ７０またはその係数は、例えば、本明細書に記載する選択および／または算出に従って得られた値に固定されてもよい。
【０１６７】
図１４Ｂに、方法Ｍ３００の実装形態Ｍ３１０のためのフローチャートを示す。方法Ｍ３１０は、初期チャネル応答推定Ｈ(0)に基づいて、修正された推定チャネル応答Ｈ(1)を得るタスクＴ１００の実装形態Ｍ１２０を含む。タスクＴ１２０は、１つまたは複数の初期推定Ｈ(0)に基づいて、なまった周波数に対しての新しいチャネル応答値を得るサブタスクを含む。タスクＴ１２０は、なまったと特定されるキャリアのそれぞれに対しての新しいチャネル応答値を得るように構成されることができる。
【０１６８】
タスクＴ１２０は、なまった周波数に隣接するキャリアに対応するＨ(0)の値（例えば、なまった周波数に隣接する集合Ｐ内のキャリアに対応する
【０１６９】
【数４５】
2009505511.tif 000047
【０１７０】
の値）から新しい値を得るように実施されることができる。例えば、タスクＴ１２０は、なまったキャリアの両側の１つまたは複数のＨ(0)の値に基づいて線形（または、高次）補間を実行するように実装されることができる。推定応答Ｈ(0)の形態に応じて、新しい値は、パイロットおよび／またはトラフィック周波数での推定チャネル応答から得られることができる。タスクＴ１２０は、また、２つ以上のなまったキャリアのギャップに跨って補間するように実施されることができる。  
図１５Ａに、受信信号
【０１７１】
【数４６】
2009505511.tif 000048
【０１７２】
の大きさの２乗対周波数のプロットを示す。また、図１５Ａに、タスクＴ１１０の実装形態に従った、しきい値Ｓ７０と
【０１７３】
【数４７】
2009505511.tif 000049
【０１７４】
の出力値との比較を示す。図示を容易にするために、１つのなまったキャリアｃのみを示すが、タスクＴ１１０は、信号
【０１７５】
【数４８】
2009505511.tif 000050
【０１７６】
のシンボル内の２つ以上のなまったキャリアを特定することができる。例えば、タスクＴ１１０は、信号
【０１７７】
【数４９】
2009505511.tif 000051
【０１７８】
のなまっている２つ以上の隣接するキャリアを特定することができる。  
図１５Ｂおよび１５Ｃに、推定チャネル応答ベクトルＨPの大きさ対周波数のプロットを示しており、黒丸は、初期推定応答ベクトル
【０１７９】
【数５０】
2009505511.tif 000052
【０１８０】
を示す。図１５Ｂに、タスクＴ１２０の実装形態に従った、修正された推定応答ベクトル
【０１８１】
【数５１】
2009505511.tif 000053
【０１８２】
（白丸）に対してのなまったキャリアにおいての新しい値を得るための線形補間処理を示す。別のキャリア（例えば、なまっていない）においては、
【０１８３】
【数５２】
2009505511.tif 000054
【０１８４】
および
【０１８５】
【数５３】
2009505511.tif 000055
【０１８６】
の値は同一である。  
【０１８７】
図１５Ｃに、タスクＴ１２０の別の実装形態に従った、修正された推定応答ベクトル
【０１８８】
【数５４】
2009505511.tif 000056
【０１８９】
（白丸）に対してのなまったキャリアにおいての新しい値を得るためのカーブフィッティング処理（例えば、高次補間処理）を示す。図１５Ｂおよび１５Ｃは、単一のなまったキャリアｃに対しての修正値の補間を示しているが、いくつかの場合では、タスクＴ１２０は、２つ以上のなまったキャリアのギャップに跨って補間することが望ましい場合がある。図示を容易にするために、図１５Ｂおよび１５Ｃは、実数値ベクトルのプロットを示すが、実際には、推定チャネル応答ベクトルＨは、典型的には、複素数値を有する。当業者なら、異なる周波数における伝送路の応答は、シンボルごとに異なって変化する場合があり、その結果、集合Ｐ内の２つの異なる周波数におけるチャネルの応答の間の比は、典型的には、経時的に変化することを理解するであろう。
【０１９０】
図１６Ａに、方法Ｍ１００の実装形態Ｍ３２０のためのフローチャートを示す。方法Ｍ３２０は、雑音推定算出のタスクＴ３００の実装形態Ｔ３６０を含む。タスクＴ３６０は、対応する伝送された値が既知であり、それらの値からシンボル内の他のキャリアに対しての雑音推定値を得るキャリアの１つの集合Ｐに対しての雑音推定値を算出するサブタスクを含む。例えば、タスクＴ３６０は、１組のパイロットキャリアに対して算出された雑音推定値
【０１９１】
【数５５】
2009505511.tif 000057
【０１９２】
に基づいた雑音推定Ｉ(1)内のトラフィックキャリアに対しての雑音推定値を補間するように実装される。  
【０１９３】
図１７Ａおよび図１７Ｂに、２種類の雑音推定ベクトルＩ(1)の大きさ対周波数のプロットを示してあり、黒丸は、やはりベクトル
【０１９４】
【数５６】
2009505511.tif 000058
【０１９５】
である値を示す。図１７Ａに、タスクＴ３６０の一実装形態に従った、
【０１９６】
【数５７】
2009505511.tif 000059
【０１９７】
で表されない周波数（白丸）におけるＩ(1)の値を得るための線形補間処理の結果を示す。図１７Ｂに、タスクＴ３６０の他の実装形態に従った、
【０１９８】
【数５８】
2009505511.tif 000060
【０１９９】
で表されない周波数（白丸）におけるＩ(1)の値を得るためのカーブフィッティング処理（例えば、高次補間処理）を示す。このような場合、
【０２００】
【数５９】
2009505511.tif 000061
【０２０１】
のゼロ値の間のキャリアに対しての値は、図１７Ｂに示すように、フィットする曲線に従って補間されることができ、
【０２０２】
【数６０】
2009505511.tif 000062
【０２０３】
のゼロ値の間のキャリアに対してゼロ値を前提にしておくことが望ましい場合がある。タスクＴ３６０は、また、補間された雑音推定をフィルタリングまたはスムージングする処理を含むことができる。
【０２０４】
図１６Ｂに、方法Ｍ３２０の実装形態Ｍ３３０のためのフローチャートを示す。方法Ｍ３３０は、雑音推定算出のタスクＴ３６０の実装形態Ｔ３７０を含む。タスクＴ３７０は、対応する伝送された値が既知であるキャリアの集合Ｐなどの、受信シンボルの一部分から基準信号を取り去ることにより、雑音推定ベクトル
【０２０５】
【数６１】
2009505511.tif 000063
【０２０６】
の値を算出するサブタスクを含む。例えば、タスクＴ３７０は、
【０２０７】
【数６２】
2009505511.tif 000064
【０２０８】
などの式に従って、周波数領域内の受信パイロット信号から基準信号を取り去るように実施されることができる。
【０２０９】
図１８に、タスクＴ１２０、Ｔ３７０およびＴ４１０を含む方法Ｍ３１０の実装形態Ｍ３４０のフローチャートを示す。図１９に、このような方法の一実施形態の段階を示す一連のプロットを示す。図１９の全てのプロットで、水平軸は、シンボル帯域幅の同一の位置を示しているが、種々のプロットの垂直軸は、全てが同じスケールではない。
【０２１０】
図１９のプロット（Ａ）は、既知の伝送ベクトルＸPの一部を示しており、プロット（Ｂ）は、受信された信号
【０２１１】
【数６３】
2009505511.tif 000065
【０２１２】
の対応する部分を示す。プロット（Ｃ）は、対応する初期推定チャネル応答ベクトル
【０２１３】
【数６４】
2009505511.tif 000066
【０２１４】
の一部を示す。この実施形態では、ベクトル
【０２１５】
【数６５】
2009505511.tif 000067
【０２１６】
は、（式４．１）に従って算出される。  
図１９のプロット（ｄ）は、しきい値Ｓ７０に対する、受信信号
【０２１７】
【数６６】
2009505511.tif 000068
【０２１８】
の出力値を示す。この場合では、
【０２１９】
【数６７】
2009505511.tif 000069
【０２２０】
の第２および第３の値が、しきい値を超えており、なまっていると特定されることが分かるであろう。プロット（Ｅ）は、新しい値（白丸）が初期チャネル応答推定の隣接する値に基づいてなまった位置において補間されている、修正された推定チャネル応答ベクトル
【０２２１】
【数６８】
2009505511.tif 000070
【０２２２】
の対応する部分を示す。（プロット（Ｃ）および（Ｅ）は、シンボルの同一の部分に対応しているが、２つのプロットの垂直軸のスケールは異なっている。）  
図１９のプロット（Ｆ）は、基準信号
【０２２３】
【数６９】
2009505511.tif 000071
【０２２４】
の一部分を示す。この実施例では、ベクトル
【０２２５】
【数７０】
2009505511.tif 000072
【０２２６】
は、周波数領域の式
【０２２７】
【数７１】
2009505511.tif 000073
【０２２８】
に従って、プロット（Ａ）および（Ｅ）に示す部分から算出される。プロット（Ｇ）は、式
【０２２９】
【数７２】
2009505511.tif 000074
【０２３０】
に従って算出される雑音推定ベクトル
【０２３１】
【数７３】
2009505511.tif 000075
【０２３２】
の対応する部分を示す。雑音推定ベクトルＩ(1)の他の要素に対しての値は、本明細書に記載されており、例えば図１７Ａおよび１７Ｂに示されるように、
【０２３３】
【数７４】
2009505511.tif 000076
【０２３４】
の値から補間されることができる。次いで、雑音推定に基づく補正信号Ｙ(1)は、例えば周波数領域の式Ｙ(1)＝Ｙ(0)−Ｉ(1)に従って得られることができる。
【０２３５】
図２０に、周波数領域内の値に作用するように構成される装置１００の実装形態３００のブロック図を示す。変換ブロックＸ１０は、装置３００の一部としてあるいは分離素子またはプロセスとして構成されることができる。基準信号発生器１１０の実装形態３１０は、（例えば受信信号Ｓ１０に基づいて）基準信号Ｓ２０のインスタンスＳ３２０を算出するように構成される。雑音推定計算機１２０の実装形態３２０は、基準信号Ｓ３２０および受信信号Ｓ１０の少なくとも一部分に基づいて、雑音推定Ｓ３０のインスタンスＳ３３０を算出するように構成される。信号補正器１３０は、雑音推定Ｓ３３０および受信信号Ｓ１０に基づいて、補正信号Ｓ４０のインスタンスＳ３４０を算出するように構成される。
【０２３６】
図２１に、装置３００の実装形態３０２のブロック図を示す。基準信号発生器３１２および雑音推定計算機３２２は、それぞれ、基準信号発生器３１０および雑音推定計算機３２０の実装形態である。
【０２３７】
基準信号発生器３１２は、（例えば、受信信号に基づいてよい）伝送路の推定応答を算出するように構成されるチャネル応答推定器１１０を含む。例えば、チャネル応答推定器１０５は、式（４．１）または本明細書に記載する別の処理に従って、初期推定チャネル応答Ｈ(0)を算出するように構成されることができる。
【０２３８】
出力計算機２５５は、受信信号の大きさの２乗またはその一部分（例えば、集合Ｐ内のキャリアに対応する値）を算出するように構成される。別の実装形態では、出力計算機２５５は、本明細書に記載するように、受信信号に時間枠を適用するように構成されることができる。比較器２６０は、なまったキャリアを特定するために、本明細書に記載するように、出力値をしきい値Ｓ７０と比較するように構成される。例えば、比較器２６０は、シンボルのどのキャリアがなまっているかの論理（例えば、バイナリ）命令を出力するように構成されることができる。
【０２３９】
補間器２７０は、比較器に従ってなまっていると特定されたキャリアに対しての新しい推定チャネル応答値を補間し、修正された推定チャネル応答ベクトル
【０２４０】
【数７５】
2009505511.tif 000077
【０２４１】
を生成するように構成されることができる。このような補間は、線形（例えば、各側の集合Ｐ内で最も近いなまっていないキャリア）であってよく、あるいは、高次多項式または別の関数に従ってよい。
【０２４２】
コンバイナ２３５は、その入力の両方に基づいた出力を生成するように構成される。この実施例では、コンバイナ２３５は、修正チャネル応答
【０２４３】
【数７６】
2009505511.tif 000078
【０２４４】
にパターンＳ５０を重ね合わせて周波数領域で基準信号を得るように構成される。パイロット信号パターンＸPなどの伝送信号の既知の部分であるパターンＳ５０は、記憶され、生成され、および／または、他の装置から受信されることができる。コンバイナ２３５は、基準信号
【０２４５】
【数７７】
2009505511.tif 000079
【０２４６】
を得るために周波数領域ベクトル
【０２４７】
【数７８】
2009505511.tif 000080
【０２４８】
とＸPとの要素ごとの乗算を実行するように構成される乗算器として実装されることができる。  
雑音推定器３２２は、受信信号Ｓ１０の対応する部分（例えば、ＹP）から基準信号Ｓ３２２を取り去って雑音推定ベクトル
【０２４９】
【数７９】
2009505511.tif 000081
【０２５０】
を生成するように構成されるコンバイナ２８０を含む。補間器２９０は、
【０２５１】
【数８０】
2009505511.tif 000082
【０２５２】
内の値から
【０２５３】
【数８１】
2009505511.tif 000083
【０２５４】
に表されないキャリア（例えば、トラフィックキャリア）に対しての値を補間することにより、雑音推定ベクトルＩ(1)を生成するように構成される。信号補正器１３０は、Ｙ(1)＝Ｙ(0)−Ｉ(1)などの式に従って補正信号Ｙ(1)を算出するように構成される。
【０２５５】
雑音推定に基づく受信信号によって伝えられるデータ値を推定することが望ましい場合がある。図２２Ａに、そのような実施形態による方法Ｍ４００のフローチャートを示す。方法Ｍ４００は、本明細書に記載する方法Ｍ１００の一実装形態および等化タスクＴ６００を含む。タスクＴ６００は、推定チャネル応答Ｈ(1)に従って補正信号Ｙ(1)を等化して、推定伝送信号Ｘ(1)を得る。一例では、タスクＴ６００は、Ｘ(1)＝（Ｈ(1)）-1Ｙ(1)、ここで、Ｈ(1)＝ｄｉａｇ（Ｈ(1)）、などの式に従って、ワンタップ等化処理を介して伝送シンボルＸの推定を得る。タスクＴ６００は、また、受信信号のなまったトラフィックキャリア（例えば、タスクＴ１１０でなまっていると特定されたパイロットの間または近傍のキャリア）によって伝えられるデータ値を推定するために、等化信号内の誤り訂正符号を復号化することを含むように実装されることができる。
【０２５６】
方法Ｍ１００の特定の実装形態に応じて、タスクＴ６００は、集合Ｐ内にないキャリアに対しての推定チャネル応答Ｈ(1)の値を得るために、補間処理（例えば、線形補間処理、または、多項式、スプラインなどのカーブフィッティング処理を介する）の実行を含む場合がある。図２３Ａに、推定チャネル応答ベクトルＨの大きさ対周波数のプロットを示す。黒丸は、集合Ｐ内のキャリアに対応する値を示しており、白丸は、線形補間を介して算出された他のキャリアに対しての値を示す。図２３Ｂに、白丸がカーブフィッティングなどの高次補間処理を介して算出された他のキャリアに対しての値を示している同様のプロットを示す。図示を容易にするために、図２３Ａおよび２３Ｂは、実数値ベクトルのプロットを示しているが、実際には、推定チャネル応答ベクトルＨは、典型的には、複素数値を有する。
【０２５７】
図２２Ｂに、方法Ｍ４００の実装形態Ｍ４１０のフローチャートを示す。方法Ｍ４１０は、補正信号に従って、チャネル応答を再推定するタスクＴ５００を含む。例えば、タスクＴ５００は、
【０２５８】
【数８２】
2009505511.tif 000084
【０２５９】
などの式に従って、補正信号および既知の伝送された値に基づいた再推定チャネル応答を算出するように構成されることができる。別法として、タスクＴ５００は、本明細書に開示する任意の他のチャネル応答の推定手法を適用するように実装されることができ、タスクＴ１００およびＴ５００が同一の手法を使用する必要はない。  
一例では、タスクＴ１００は、式
【０２６０】
【数８３】
2009505511.tif 000085
【０２６１】
に従って推定
【０２６２】
【数８４】
2009505511.tif 000086
【０２６３】
を得て、タスクＴ５００は、式
【０２６４】
【数８５】
2009505511.tif 000087
【０２６５】
に従ってチャネル応答を再推定する。別の実施例では、タスクＴ１００は、本明細書に記載する逆変換または予測手法を使用して推定
【０２６６】
【数８６】
2009505511.tif 000088
【０２６７】
を得て、タスクＴ５００は、式
【０２６８】
【数８７】
2009505511.tif 000089
【０２６９】
に従ってチャネル応答を再推定する。タスクＴ５００は、集合Ｐ内にないキャリア（例えば、トラフィックチャネル）に対してのＨ(2)の値を得るためのサブタスクを含むことができる。このサブタスクは、補間処理を含むことができ、Ｈ(2)の値も、利用可能な場合、Ｈ(1)からのトラフィックチャネル値に基づくことができる。
【０２７０】
方法Ｍ４１０は、また、推定伝送信号Ｘ(1)を得るために再推定チャネル応答Ｈ(2)に従って補正信号Ｙ(1)を等化する、本明細書に記載する等化タスクＴ６００の実装形態Ｔ６１０を含む。一例では、タスクＴ６１０は、Ｘ(1)＝（Ｈ(2)）-1Ｙ(1)、ここで、Ｈ(2)＝ｄｉａｇ（Ｈ(2)）、などの式に従って、ワンタップ等化処理を介して伝送信号Ｘの推定を得る。タスクＴ６１０は、集合Ｐ内にないキャリア（例えば、トラフィックキャリア）に対してのＨ(2)の値を得るサブタスクを含むことができる。このサブタスクは、補間処理を含むことができ、Ｈ(2)の値も、利用可能な場合、Ｈ(1)からのトラフィックチャネル値に基づくことができる。
【０２７１】
図２２Ｃに、一実施形態による方法Ｍ５００のフローチャートを示す。方法Ｍ５００では、タスクＴ２００、Ｔ３００およびＴ４００が、推定チャネル応答として
【０２７２】
【数８８】
2009505511.tif 000090
【０２７３】
（または、Ｈ(2)）を使用して繰り返される。再推定および反復のこのようなプロセスは、現行の推定チャネル応答を使用して等化タスクＴ６１０が実行されるまで、所望の回数だけ繰り返されてよい。
【０２７４】
図２４に、チャネル応答が最後の反復の後に再推定されないという点で方法Ｍ５００と異なる実施形態による方法Ｍ６００のフローチャートを示す。本明細書に記載する装置１００は、また、例えば、等化器（ワンタップ等化器など）、ならびに、反復制御およびデータ信号方向のための適切な決定論理を含むことにより、方法Ｍ４００、Ｍ５００、またはＭ６００の実装形態を実行するように構成されることができる。
【０２７５】
チャネル応答の再推定タスクＴ５００の異なる反復では、異なるチャネル推定手法を使用してよい。また、方法Ｍ１００の異なる実装形態が、異なる反復において使用されてよい。例えば、１つの反復において、方法Ｍ２００の一実装形態（例えば、時間領域内で雑音推定を算出する）を含んでよく、一方、先行するならびに／または後続の反復においては、方法Ｍ３００の一実装形態（例えば、周波数領域内で雑音推定を算出する）を含んでよい。
【０２７６】
上述する実施形態は、雑音推定に従って受信信号を補正するように構成される方法および装置を含む。さらなる実施形態による方法Ｍ７００は、受信信号の算出モデルの値に従って受信信号Ｓ１０を補正する。
【０２７７】
図２５Ａに、方法Ｍ７００の実装形態のためのフローチャートを示す。タスクＴ１００は、本明細書に記載する任意のチャネル応答推定手順に従って、伝送路の推定応答Ｈaを得る。タスクＴ７１０は、受信信号Ｓ１０によって伝えられるシンボルＸを推定する。タスクＴ７２０は、推定シンボルの１つまたは複数のなまった値を特定する。タスクＴ７３０は、推定シンボルと、伝送路の推定応答と、基準信号の集合とに基づく受信信号のモデルｓmを算出する。タスクＴ７４０は、なまった値の位置に基づいてならびに算出されたモデルｓmの値に従って、受信信号Ｓ１０を補正する。
【０２７８】
タスクＴ７１０は、伝送路の推定応答を受信信号Ｙに適用することによりシンボルＸを推定するように構成されることができる。例えば、タスクＴ７１０は、受信信号に等化処理を実行するように構成されることができる。このような１つ実施例では、タスクＴ７１０は、Ｓ＝（Ｈa）-1Ｙ、ここでＨa＝ｄｉａｇ（Ｈa）、などの式に従って、ワンタップ等化処理を介して伝送シンボルＸの推定Ｓを算出する。
【０２７９】
タスクＴ７２０は、推定シンボルＳの時間領域の値ｓkをしきい値Ｓ８０と比較することにより、推定シンボルの１つまたは複数のなまった値を特定するように構成されることができる。しきい値Ｓ８０との比較は、推定シンボルの各値ｓkに対して実行されてよい。一例では、しきい値Ｓ８０は、推定シンボルｓの各時間領域の値の出力（例えば、大きさの２乗）と比較される。しきい値Ｓ８０より大きい（あるいは、別法として、しきい値Ｓ８０以上の）値は、なまっていると特定される。別の実施例では、この比較は、出力以外（値の対数など）の時間領域の各値の測定によって実行されることができる。
【０２８０】
別法として、しきい値Ｓ８０は、時間枠内の値の大きさの２乗の合計などの、ある時間枠内の１組の値ｓkに関連する測定値と比較されることができる。時間枠の幅（これは、シンボルの幅より大きい、と等しい、あるいは、より小さい）は、インパルス雑音事象の推定期間に従って選択されることができる。１つのこのような実施例では、しきい値Ｓ８０は、その時間枠の幅に亘っての推定される非決定論的成分の平均信号出力と比較される。別のこのような実施例では、しきい値Ｓ８０は、その時間枠の幅に亘っての推定される非決定論的成分の合計信号出力と比較される。時間枠の測定値がしきい値Ｓ８０より大きい（あるいは、別法として、しきい値Ｓ８０以上の）場合、その時間枠内の値はなまっていると特定される。別法として、値が、なまっていると特定された少なくとも一定数の時間枠内にある場合、その値はなまっていると特定されることができる。値が２つ以上の時間枠内にあるようにするために、隣接する時間枠は重複してよく、あるいは、各値が１つの時間枠内のみにあるようにするために重複しないでよい。
【０２８１】
１つまたは複数の信号条件に従ってしきい値Ｓ８０の値を得ることが望ましい場合がある。しきい値Ｓ８０は、受信信号Ｓ１０の出力の測定値から、あるいは、信号Ｓ１０の一部（例えば、受信パイロット信号）の出力の測定値から、得られることができる。例えば、しきい値Ｓ８０は、サンプルごとの平均チャネル出力の推定
【０２８２】
【数８９】
2009505511.tif 000091
【０２８３】
に基づいてよく、これは、以下の式などの式に従って算出されることができる。
【０２８４】
【数９０】
2009505511.tif 000092
【０２８５】
平均は、Ｈaの値が利用できるシンボルの全てのサブキャリアに亘って、あるいは、シンボルのサブキャリアの集合Ｐ（例えば、パイロット信号の集合）に亘って得られることができる。式（５．１）での平均化は、平均値または中央値を使用して算出されることができる。別法として、しきい値Ｓ８０は、別の統計フィルタを使用して算出されるチャネル出力の測定値に基づいてよい、式（５．１）は、精製された形のＨa（フィルタリングならびに／または補間された形など）、先行モデルのＨa、または、Ｈaとは別に算出された推定などの、Ｈa以外のチャネル応答推定を使用して算出されることもできる。  
しきい値Ｓ８０は、荷重係数を含んでよい。一例では、しきい値Ｓ８０は値
【０２８６】
【数９１】
2009505511.tif 000093
【０２８７】
を有し、ここで、ｗ1は、高出力の信号のための間違い警報の確率に対する所望の検出率に平衡を保たせるように選択されることができる荷重係数である。一例では、ｗ1の値は２５である。
【０２８８】
１組の所定の値の中からの選択に従ってしきい値Ｓ８０の値を変化させることが望ましい場合がある。例えば、ｗ1の値は、電流信号特性とモデルのセットの１つとの間の適合に従って選択されることができる。別法として、しきい値Ｓ８０またはその係数は、例えば、本明細書に記載する選択および／または算出に従って得られた値に固定されてもよい。
【０２８９】
タスクＴ７２０は、なまっていると特定された推定シンボル内の時間領域位置の命令を出力するように構成されることができる。一例では、タスクＴ７２０は、ｋ∈｛０，１，．．．，Ｍ−１｝に対する以下の式に従って、２値ベクトルまたはマスクｑを算出する。
【０２９０】
【数９２】
2009505511.tif 000094
【０２９１】
推定シンボル内のなまった値の位置に関する情報は、受信信号Ｓ１０の対応するサンプルの値を修正するのに使用されることができる。１つの単純な手法では、ゼロの値を、推定シンボルのなまった位置に対応する受信信号Ｓ１０のサンプルに割り当てる。しかし、受信信号Ｓ１０の単純な出力ゼロ化は、特に、なまっているサンプル位置の分類が不正確である場合に、歪を招く場合がある。
【０２９２】
タスクＴ７３０は、なまった位置におけるＳ１０のサンプルを修正するのに使用されることができる受信信号Ｓ１０のモデルを算出する。タスクＴ７３０は、基準信号の集合に基づいて、１組の周波数領域の値を選択する（あるいはそうでない場合は、算出する）サブタスクを含むように実装されることができる。例えば、タスクＴ７３０は、推定シンボルＳの各値を基準信号の集合の１つの写像するように構成されることができる。このような場合、タスクＴ７３０は、推定シンボルＳの各値に対して、ある距離尺度に従って基準信号の集合の中の最も近い基準値を選択することにより、写像を行うように構成されることができる（例えば、ユークリッド距離）。
【０２９３】
ベースバンドデジタル通信信号は、典型的には、１つの有限集合の状態を有しながら伝送される。基準信号の集合が、この集合の可能性のある伝送状態を表すことが望ましい場合がある。例えば、基準信号の集合は、信号を伝送するのに使用される変調方式により決定されるコンステレーション点の１つの集合であってよい。図２６に、１６−ＱＡＭ変調方式のためのＩ−Ｑ平面内のコンステレーション点の集合の一例を示す。また、図２６に、１６の対応する領域へのＩ−Ｑ平面の区分を示しており、平面内の任意の点が、各軸に沿ったその領域の境界値とその点のＩおよびＱ値とを比較することにより、素早く適したコンステレーション点に写像されることができるようにする。このような形では、タスクＴ７３０は、推定シンボルＳの各値に対して、コンステレーション点の集合から周波数領域値を選択するように構成されることができる。
【０２９４】
タスクＴ７３０は、また、伝送路の応答の推定Ｈbを、選択された（あるいはそうでない場合は、算出された）周波数領域値に適用するサブタスクを含むことができる。推定Ｈbは、チャネル応答推定Ｈaと同一であってよく、あるいは、推定Ｈaの精製された形または先行モデルであってよく、あるいは、推定Ｈaとは別に算出されてよい。一例では、推定Ｈbは、フィルタリング（ＩＦＦＴ／ＦＦＴフィルタリングなど）および／または補間（例えば、パイロットキャリアにおける値から）によって推定Ｈaから得られる。図２５Ｂに、２つのサブタスクを有するタスクＴ７３０の実装形態Ｔ７３２を含む方法Ｍ７００の実装形態Ｍ７１０のためのフローチャートを示す。
【０２９５】
受信信号Ｓ１０に基づいて、タスクＴ７４０は、算出されたモデルの値に従って補正信号を算出する。タスクＴ７４０は、また、推定シンボルのなまった値の位置に基づいて受信信号Ｓ１０を補正するように実装されることができる。タスクＴ７４０の実装形態Ｔ７４２は、なまった値に対応する受信信号のサンプルをモデルの対応する値に入れ替えることにより受信信号Ｓ１０を補正する。１つのこのような実施例では、タスクＴ７４２は、以下の式などに従って補正信号ｙ(1)を算出する。
【０２９６】
【数９３】
2009505511.tif 000095
【０２９７】
方法Ｍ１００に関連して図２２Ａ、Ｂ、Ｃおよび２４に示すように、方法Ｍ７００も、補正信号の等化および／または反復などのさらなる処理を含むように実装されることができる。例えば、方法Ｍ７００の１つまたは複数のタスクは、補正信号Ｓ７４０で反復されてよく（例えば、改良されたチャネル推定および／またはシンボル推定を適用する）、この反復は、等化処理の前に実行されることができる。方法Ｍ７００は、また、受信信号Ｓ１０の同一シンボルおよび／または連続するシンボルに対してのどちらでも、本明細書に記載する方法Ｍ１００、Ｍ４００、Ｍ５００またはＭ６００の実装形態とあいまって（例えば、前にあるいは後で）実行されてよい。別の実施例では、雑音推定のタスクＴ３２０は、タスクＴ７３０で算出されるモデルの対応する値に基づくしきい値を適用するように構成される。
【０２９８】
図２７に、情報信号Ｓ１０を受信して対応する補正信号Ｓ７４０を生成するように構成される実施形態による装置７００のブロック図を示す。チャネル応答推定器１０５は、本明細書に開示するチャネル応答推定手順に従って伝送路の推定応答を算出するように構成される。シンボル推定器７１０は、推定シンボルＳ７１０を算出するように構成される。例えば、シンボル推定器７１０は、本明細書に開示するタスクＴ７１０の実装形態を実行することができる。雑音検出器７２０は、推定シンボルＳ７１０の１つまたは複数のなまった値の位置Ｓ７２０を示すように構成される。例えば、雑音検出器７２０は、本明細書に開示するタスクＴ７２０の実装形態を実行することができる。モデル計算機７３０は、受信信号のモデルＳ７３０を算出するように構成される。例えば、モデル計算機７３０は、本明細書に開示するタスクＴ７３０の実装形態を実行することができる。信号補正器７４０は、受信信号Ｓ１０、示された位置Ｓ７２０、および、算出されたモデルＳ７３０に基づいて補正信号Ｓ４０を生成するように構成される。例えば、信号補正器７４０は、本明細書に開示するタスクＴ７４０の実装形態を実行するように構成されることができる。
【０２９９】
装置７００の種々の素子は、例えば、同一のチップまたは１つのチップセット内の２つ以上のチップに存在する電子デバイスおよび／または光学デバイスとして実装されることができるが、そのような制限のない他の配置も考慮される。装置７００の１つまたは複数の素子は、マイクロプロセッサ、埋め込みプロセッサ、ＩＰコア、デジタル信号プロセッサ、ＦＰＧＡ（フィールドプログラマブルゲートアレイ）（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙｓ））、ＡＳＳＰ（アプリケーションスペシフィックスタンダードプロダクト）（ａｐｐｌｉｃａｔｉｏｎ−ｓｐｅｃｉｆｉｃ  ｓｔａｎｄａｒｄ  ｐｒｏｄｕｃｔｓ））、ＡＳＩＣ（アプリケーションスペシフィックインテグレーテッドサーキット）（ａｐｐｌｉｃａｔｉｏｎ−ｓｐｅｃｉｆｉｃ  ｉｎｔｅｇｒａｔｅｄ  ｃｉｒｃｕｉｔｓ））などの、論理素子（例えば、トランジスタ、ゲート）の固定アレイまたはプログラム可能な１つまたは複数のアレイ上で実行される１つまたは複数の集合の命令として、全体としてあるいは部分的に、実装されることができる。また、１つまたは複数のこのような素子は、共通の構造を有することができる（例えば、異なる時間に異なる素子に対応するコードの一部を実行するのに使用されるプロセッサ、または、異なる時間に異なる素子に対して処理を実行する電子および／または光学デバイスの配置）。
【０３００】
図２８に、装置７００の実装形態７０２のブロック図を示す。シンボル推定器７１２、雑音検出器７２２、モデル計算機７３２、および信号補正器７４２は、それぞれ、シンボル推定器７１０、雑音検出器７２０、モデル計算機７３０、および信号補正器７４０の実装形態である。
【０３０１】
シンボル推定器７１２は、推定チャネル応答に従って受信信号Ｓ１０を等化するように構成されるコンバイナを含む。この実施例では、コンバイナ７１２は、チャネル応答推定Ｈaに受信信号Ｓ１０を重ね合わせて周波数領域内で推定シンボルＳ７１２を得るように構成される。コンバイナ７１２は、チャネル推定Ｈaにより周波数領域ベクトルＹの要素ごとの区分を（連続してならびに／または並行して）実行して推定シンボルＳを得るように構成される乗算器を含むように実装されることができる。
【０３０２】
（上述した）逆変換ブロックＸ２０は、周波数領域内で推定シンボルＳに逆周波数変換を実行して時間領域内で値ｓを得るように構成される。逆変換ブロックＸ２０は、また、分離素子またはプロセスとして構成されることができ、装置７００のインスタンスを含むシステムが、たとえ装置７００が起動されていない場合でも、ブロックＸ２０を使用して逆周波数変換を実行できるようにする。
【０３０３】
出力計算機２５０は、それぞれの時間領域推定シンボル値ｓkの大きさの２乗を算出するように構成される。別の実装形態では、出力計算機２５０は、本明細書に開示する推定シンボル値に時間枠を適用するように構成されることができる。比較器７２４は、出力値をしきい値Ｓ８０と比較してなまった値を特定するように構成される。例えば、比較器７２４は、式（５．２）に記載されるような処理を実行するように構成されることができる。比較器７２４は、連続してならびに／または並行してシンボル値に作用するように構成されることができる。
【０３０４】
モデル計算機７３２は、バリューセレクタ７３４およびコンバイナ７３６を含む。バリューセレクタ７３４は、基準信号の１つ集合から、コンステレーション点の１つの集合などの周波数領域値を選択するように構成される。例えば、バリューセレクタ７３４は、基準値と推定シンボルの値との距離を算出するように構成されることができる。別法として、バリューセレクタ７３４は、図２６の実施例に示す領域境界の集合に従って推定シンボルの値を分類するように構成されることができる。
【０３０５】
コンバイナ７３６は、選択された周波数領域値に推定チャネル応答を適用するように構成される。コンバイナ７３６は、算出されたモデルのリプレゼンテーション（ｒｅｐｒｅｓｅｎｔａｔｉｏｎ）Ｓmを得るために周波数領域値と推定チャネル応答との要素ごとの乗算を実行する（連続してならびに／または並行して）ように構成される乗算器として実装されることができる。逆変換ブロックＸ４０は、１つのブロックＸ２０あるいは複数のブロックＸ２０の別のインスタンスとして実装されることができ、Ｘ４０は、異なる時間に異なる信号に作用する同一の構造またはプロセスとして実装されることができる。
【０３０６】
いくつかの実装形態は、シンボル推定器７１２によりならびにコンバイナ７３６により使用される推定チャネル応答が同一であるあるいは別々に得られるように、構成されるが、装置７００は、これら２つの要素が異なるチャネル推定に適用されるように実装されることもできる。信号補正器７４２は、例えば上述の式（５．３）に従って受信信号Ｓ１０に選択的置換処理を実行するように実装される。  
【０３０７】
方法Ｍ２００の一例での処理を説明する式は、以下のように反復実行のｎ回目の繰り返しに一般化されることができる：
【０３０８】
【数９４】
2009505511.tif 000096
【０３０９】
上式で、関数ｆは、タスクＴ３２０に関連して本明細書に記載するなまった値を特定する処理を示す。  
【０３１０】
方法Ｍ３００の一例の処理を説明する式は、以下のように反復実行のｎ回目の繰り返しに一般化されることができる：
【０３１１】
【数９５】
2009505511.tif 000097
【０３１２】
上式で、関数ｇ1およびｇ2は、それぞれ、タスクＴ１１０に関連して本明細書に記載するなまったキャリアを特定する処理、および、タスクＴ３１０に関連して本明細書に記載する雑音推定の微分処理を示す。  
【０３１３】
方法Ｍ７００の一例の処理を説明する式は、以下のように反復実行のｎ回目の繰り返しに一般化されることができる：
【０３１４】
【数９６】
2009505511.tif 000098
【０３１５】
上式で、関数ｒ１、ｒ２およびｒ３は、それぞれ、タスクＴ７２０に関連して本明細書に記載するなまった値を特定する処理、タスクＴ７３０に関連して本明細書に記載するモデルを算出する処理、および、タスクＴ７４０に関連して本明細書に記載する置換処理を示す。
【０３１６】
伝送シンボルＸに対しての推定が得られると、最良のコンステレーション点を選択するために（例えば、最小距離評価に従って）、推定される伝送シンボルの各要素に対して決定処理が実行されることができる（例えば、サブキャリアに対して使用される変調方式に従って）。
【０３１７】
本明細書に記載する方法および装置は、対応する受信シンボルに基づいて補正信号を算出するのに使用されることができる。このようなシンボルごとの処理は、シンボルのそれぞれの受信ストリームに対して繰り返されることができる。チャネル応答推定の場合、本明細書に記載するように、予め受信されたシンボルからの情報も使用されることができる。本明細書に記載する方法および装置の別の実装形態では、先行するシンボルからの雑音推定に関する情報も使用されることができる。経時的に繰り返すパルス様雑音事象の場合、例えば、これらの方法または装置は、そのような事象の区間を特定して探知するように構成されることができる。これらの情報は、特定のシンボルに対しての方法Ｍ１００の適切な実装形態を予測するのに適用されることができ、処理サイクルをより適切に経時的に配分することにより電力消費を低下させるのを助けることができる。
【０３１８】
狭帯域干渉事象は、経時的に相互に関連付けられることもできる。例えば、このような事象は、隣接するシンボル内のほぼ同一のキャリアの集合に影響を与える場合があり、および／または、ある期間に従ってこれらの周波数上で繰り返す場合がある。このような場合、１つまたは複数の以下のシンボルに対する処理において雑音推定情報を使用することが望ましい場合がある。例えば、本明細書に記載する方法または装置のこのような実装形態は、シンボルに対しての初期および／または修正チャネル推定を算出する前に後続の受信シンボルから雑音推定ＩまたはＩPを取り去るように構成されることができる。別の実施例では、本明細書に記載する方法または装置の実装形態は、しきい値Ｓ６０、Ｓ７０および／またはＳ８０の算出により先行するシンボル内でなまっていると特定されたパイロットを除外するように構成されることができる。このような手法は、特に、チャネル応答がそうでなければ経時的に比較的一定である場合には、なまった周波数の改良された復号化を補助することができる。
【０３１９】
一実施形態による方法または装置が適用される特定の方式に応じて、パイロット成分の伝送出力は、トラフィック成分の伝送出力と異なってよい。ＤＶＢシステムでは、例えば、パイロットキャリアの出力は、トラフィックキャリアの出力の４／３に等しい。パイロットキャリア値から、雑音推定および／または推定チャネル応答値などの、トラフィックキャリア値の補間におけるこれらのファクタを明らかにすることが望ましい場合がある。
【０３２０】
別の実施形態では、方法Ｍ２００の実装形態が実行され、ここでは、基準信号または決定論的成分がゼロであると想定される。例えば、タスクＴ１００およびＴ２００が省略されてよく、雑音推定が受信信号から直接算出されるようにする（例えば、時間領域内で受信信号をしきい値化することによって）。その結果、この方法で得られるチャネル推定は、方法Ｍ１００、Ｍ４００、Ｍ５００、Ｍ６００、またはＭ７００のインスタンスにおいて初期チャネル推定として適用される。
【０３２１】
方法Ｍ５００、またはＭ６００の実装形態は、反復において等化処理を含むことができ、第２または以降のパスでの伝送信号Ｘのベクトルが、トラフィック信号さらにはパイロット信号からの成分を含むことができるようにする。また、先行する反復からのＸの推定に対して誤り訂正符号化およびそれに続くリエンコード（例えば、誤り検出、エラー訂正、および／または、エンコード、パンクチュアリングおよび／またはインターリビング）を実行することによりＸに対しての新しい推定を得ることが望ましい場合がある。
【０３２２】
本明細書に記載するインパルス雑音推定の方法が実行されたかどうかならびに／またはどのように実行されたかは、受信信号の１つまたは複数の特性および／または利用可能な供給出力またはエネルギーなどの、処理環境の他のファクタに応じて選択されることができる。例えば、別の実施形態は、受信信号の検出されるドップラー周波数に基づいて１つまたは複数のタスクおよび／または素子がアクティブ化または非アクティブ化される（例えば、方法Ｍ５００またはＭ６００の実装形態におけるさらなる反復）ような、ならびに／または、受信信号の検出されるドップラー周波数に基づいて１つまたは複数のタスクまたは素子の複雑さ（例えば、補間処理の複雑さまたはチャネル応答推定タスクの性質）が変化するようなシステム、方法および装置を含む。このような実施形態は、検出されるドップラー周波数に基づいてインパルス雑音推定の方法Ｍ１００またはＭ７００をアクティブ化または非アクティブ化することができる。少なくともいくつかの実施形態では、このような選択は、電力節約に利点をもたらすことができ、これは、典型的には、電池式（例えば、ハンドヘルド）でさらに熱伝達が制限される用途において重要な問題である。
【０３２３】
装置１００または７００の１つまたは複数の素子は、トランジスタおよび／またはゲートなどの論理素子のアレイとして実装されることができる。素子は、一体に１つのチップにあるいは２つ以上のチップに跨って実装されることができる。１つまたは複数のこれらの素子を含むチップは、ビタビ復号、リード−ソロモンの前進型誤信号訂正などの他の処理機能を受信信号に実行するように構成されるアレイを含むこともできる。
【０３２４】
装置１００または７００の実装形態は、１つまたは複数の構造を有する素子を使用して、ベクトル、行列、ならびに、他の計算および決定を補助するために、対応する加算、引き算、乗算、割り算、比較、加重、および、成分の他の数学的および／または論理的処理を実行するように構成されることができる。このような構造は、ハードウェア、ソフトウェアおよび／またはファームウェアに実装されることができる。このような構造または素子のいくつかまたは全ては、１つのデバイスまたはシステム内で他の方法および／または装置と共に共用されることもできる。
【０３２５】
記載する実施形態の先の提示は、任意の当業者が本発明を製作または使用することを可能にするために与えられた。これらの実施形態の種々の修正形態が可能であり、本明細書に示される包括的な原理も、他の実施形態に適用されることができる。例えば、本発明は、部分的にあるいは全体として、ハードワイヤード回路として、特定用途向けＩＣの中で作られる回路構成として、あるいは、非揮発性記憶装置へロードされるファームウェアプログラム、または、データ記憶媒体からもしくはデータ記憶媒体へロードされる機械可読コード（このコードは、マイクロプロセッサまたは他のデジタル信号処理ユニットなどの論理素子のアレイにより実行可能な命令である）などのソフトウェアプログラムとして実装されることができる。したがって、本発明は、上記に示す実施形態に限定されることを意図しておらず、むしろ、本明細書に任意の形式で開示する原理および新規の特徴に従う最も広い範囲と一致する。
【図面の簡単な説明】
【０３２６】
【図１】ＯＦＤＭシステムの実施例を示すブロック図である。
【図２Ａ】逆変換ＦＦＴ処理の入力および出力を示す図である。
【図２Ｂ】ＦＦＴ処理の入力および出力を示す図である。
【図３】連続パイロットおよび分散パイロット（ｓｃａｔｔｅｒｅｄ  ｐｉｌｏｔｓ）を含むパイロット方式の一部分を示す図である。
【図４Ａ】一実施形態による方法Ｍ１００のためのフローチャートである。
【図４Ｂ】方法Ｍ１００の実装形態Ｍ１１０のためのフローチャートである。
【図５Ａ】方法Ｍ１００の実装形態Ｍ１２０のためのフローチャートである。
【図５Ｂ】方法Ｍ１１０の実装形態Ｍ１３０のためのフローチャートである。
【図６Ａ】一実施形態による装置１００のブロック図である。
【図６Ｂ】装置１００の実装形態１０２のブロック図である。
【図７Ａ】方法Ｍ１００の実装形態Ｍ２００のためのフローチャートである。
【図７Ｂ】方法Ｍ２００の実装形態Ｍ２１０のためのフローチャートである。
【図８Ａ】パイロットベクトルＸPの振幅対周波数のプロットを示す図である。
【図８Ｂ】基準信号ＹRの振幅対周波数のプロットを示す図である。
【図８Ｃ】非決定論的成分ｙNの大きさの２乗対時間のプロットを示す図である。
【図９】装置１００の実装形態２００のブロック図である。
【図１０】装置２００の実装形態２０２のブロック図である。
【図１１Ａ】チャネル応答推定器１１０の実装形態１１２のブロック図である。
【図１１Ｂ】チャネル応答推定器１１０の実装形態１１４のブロック図である。
【図１２Ａ】予測フィルタ４１０の実装形態４１２のブロック図である。
【図１２Ｂ】チャネル応答推定器１１０の実装形態１１６のブロック図である。
【図１３Ａ】チャネル応答推定器１１０の実装形態１１８のブロック図である。
【図１３Ｂ】補間器５２０のブロック図である。
【図１４Ａ】方法Ｍ１００の実装形態Ｍ３００のためのフローチャートである。
【図１４Ｂ】方法Ｍ３００の実装形態Ｍ３１０のためのフローチャートである。
【図１５Ａ】受信信号部分ＹPの大きさの２乗対時間のプロットを示す図である。
【図１５Ｂ】推定されるチャネル応答ベクトルＨPの大きさ対周波数のプロットおよび線形補間処理を示す図である。
【図１５Ｃ】推定されるチャネル応答ベクトルＨPの大きさ対周波数のプロットおよびカーブフィッティング処理を示す図である。
【図１６Ａ】方法Ｍ１００の実装形態Ｍ３２０のフローチャートである。
【図１６Ｂ】方法Ｍ１００の実装形態Ｍ３３０のフローチャートである。
【図１７Ａ】雑音推定ベクトルＩPの大きさ対周波数のプロットおよび線形補間によりＩPから得られるベクトルＩの付加的な値を示す図である。
【図１７Ｂ】雑音推定ベクトルＩPの大きさ対周波数のプロットおよびカーブフィッティング処理によりＩPから得られるベクトルＩの付加的な値を示す図である。
【図１８】方法Ｍ３１０の実装形態Ｍ３４０のフローチャートである。
【図１９Ａ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｂ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｃ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｄ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｅ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｆ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図１９Ｇ】方法Ｍ３４０の一例の段階を示す一連のプロットを示す図である。
【図２０】装置１００の実装形態３００のブロック図である。
【図２１】装置３００の実装形態３０２のブロック図である。
【図２２Ａ】一実施形態による方法Ｍ４００のフローチャートである。
【図２２Ｂ】方法Ｍ４００の実装形態Ｍ４１０のフローチャートである。
【図２２Ｃ】一実施形態による方法Ｍ５００のフローチャートである。
【図２３Ａ】推定されるチャネル応答ベクトルＨPの大きさ対周波数のプロットおよび線形補間によりＨPから得られるベクトルＨの付加的な値を示す図である。
【図２３Ｂ】推定されるチャネル応答ベクトルＨPの大きさ対周波数のプロットおよびカーブフィッティング処理によりＨPから得られるベクトルＨの付加的な値を示す図である。
【図２４】一実施形態による方法Ｍ６００のフローチャートである。
【図２５Ａ】一実施形態による方法Ｍ７００のフローチャートである。
【図２５Ｂ】方法Ｍ７００の実装形態Ｍ７１０のフローチャートである。
【図２６】コンステレーション点の集合および対応する領域内のＩ−Ｑ平面の区分を示す図である。
【図２７】一実施形態による装置７００のブロック図である。
【図２８】装置７００の実装形態７０２のブロック図である。
【図１】
2009505511.tif 000099
【図２Ａ】
2009505511.tif 000100
【図２Ｂ】
2009505511.tif 000101
【図３】
2009505511.tif 000102
【図４Ａ】
2009505511.tif 000103
【図４Ｂ】
2009505511.tif 000104
【図５Ａ】
2009505511.tif 000105
【図５Ｂ】
2009505511.tif 000106
【図６Ａ】
2009505511.tif 000107
【図６Ｂ】
2009505511.tif 000108
【図７Ａ】
2009505511.tif 000109
【図７Ｂ】
2009505511.tif 000110
【図８Ａ】
2009505511.tif 000111
【図８Ｂ】
2009505511.tif 000112
【図８Ｃ】
2009505511.tif 000113
【図９】
2009505511.tif 000114
【図１０】
2009505511.tif 000115
【図１１Ａ】
2009505511.tif 000116
【図１１Ｂ】
2009505511.tif 000117
【図１２Ａ】
2009505511.tif 000118
【図１２Ｂ】
2009505511.tif 000119
【図１３Ａ】
2009505511.tif 000120
【図１３Ｂ】
2009505511.tif 000121
【図１４Ａ】
2009505511.tif 000122
【図１４Ｂ】
2009505511.tif 000123
【図１５Ａ】
2009505511.tif 000124
【図１５Ｂ】
2009505511.tif 000125
【図１５Ｃ】
2009505511.tif 000126
【図１６Ａ】
2009505511.tif 000127
【図１６Ｂ】
2009505511.tif 000128
【図１７Ａ】
2009505511.tif 000129
【図１７Ｂ】
2009505511.tif 000130
【図１８】
2009505511.tif 000131
【図１９Ａ】
2009505511.tif 000132
【図１９Ｂ】
2009505511.tif 000133
【図１９Ｃ】
2009505511.tif 000134
【図１９Ｄ】
2009505511.tif 000135
【図１９Ｅ】
2009505511.tif 000136
【図１９Ｆ】
2009505511.tif 000137
【図１９Ｇ】
2009505511.tif 000138
【図２０】
2009505511.tif 000139
【図２１】
2009505511.tif 000140
【図２２Ａ】
2009505511.tif 000141
【図２２Ｂ】
2009505511.tif 000142
【図２２Ｃ】
2009505511.tif 000143
【図２３Ａ】
2009505511.tif 000144
【図２３Ｂ】
2009505511.tif 000145
【図２４】
2009505511.tif 000146
【図２５Ａ】
2009505511.tif 000147
【図２５Ｂ】
2009505511.tif 000148
【図２６】
2009505511.tif 000149
【図２７】
2009505511.tif 000150
【図２８】
2009505511.tif 000151
【国際調査報告】
2009505511.tif 000152