JP2001513954A - 特異値分解を用いたチャンネルインパルス応答の推定 - Google Patents

Abstract

(57)【要約】 本発明は、ＣＤＭＡ受信器におけるＣＩＲ（チャンネルインパルス応答）及びＳＩＮＲ（信号対干渉及び雑音比）の推定に係る。より詳細には、本発明は、通信システムのチャンネルインパルス応答(ＣＩＲ)、例えば、デジタル移動無線ネットワーク（ＧＳＭネットワーク）の無線チャンネルのＣＩＲを決定する方法に係る。特に、本発明は、既知のトレーニングシーケンスを受け取るのに基づいてＣＩＲを決定することに係る。又、本発明は、干渉打消しへの適用及びＣＤＭＡ受信器への使用にも係る。本発明は、ネットワーク分析により得られたマトリクスを「対角マトリクス」へと最初に変換し、そしてそれにより得られたマトリクスを反転することにより、公知技術に関連した問題を軽減できるという実現から達成される。更に、本発明は、サイドローブが除去されるまで、又はそれらがほとんど影響しなくなるように、反転マトリクスにおいて固有値の逆数に減少する数値を乗算する。

Description

【発明の詳細な説明】 特異値分解を用いたチャンネルインパルス応答の推定発明の分野 本発明は、ＣＤＭＡ受信器におけるＣＩＲ（チャンネルインパルス応答）及び ＳＩＮＲ（信号対干渉及び雑音比）の推定に係る。より詳細には、本発明は、通 信システムのチャンネルインパルス応答(ＣＩＲ)、例えば、デジタル移動無線ネ ットワーク（ＧＳＭネットワーク）の無線チャンネルのＣＩＲを決定する方法に 係る。特に、本発明は、既知のトレーニングシーケンスを受け取るのに基づいて ＣＩＲを決定することに係る。又、本発明は、干渉打消しへの適用及びＣＤＭＡ 受信器への使用にも係る。先行技術の説明 ＣＩＲを決定するためには、送信信号の一部分が分からねばならない。ＧＳＭ ネットワークの場合に、同期バースト（ＳＢ）が信号の有用な部分である。ＳＢ は、各ベースステーションから少なくとも１つのチャンネルを経て送信され、そ してそれらは規則的なパターンで送信される。ＧＳＭプロトコルをデコードする 必要はない。ＳＢにおけるデータと、それらが発生するパターンは、両方とも、 固定であり、そして全てのベースステーションに対して実質的に同一である。 ＳＢを用いてＣＩＲを決定する利点は、それらが比較的長いノイズ状の所定の 送信信号を表わすことである。典型的に、６４ビットが２３７μ秒の周期にわた って送信される。従って、ＣＩＲを決定する場合には、ＳＢを抽出してＣＩＲを 決定するように処理できるに充分なほど、受信信号のバーストに同期することが 必要となる。 この処理は、推定技術を使用することにより行なわれる。ＣＩＲを推定するた めに、既知のトレーニングシーケンスが送信されＳ_tx(ｔ)、これは通信チャンネ ルによって歪まされ、受信信号Ｓ_rx(ｔ)を発生する。ＣＩＲを推定する際の問題 は、ＦＩＲフィルタを通過した後の既知のＳ_tx(ｔ)を受信信号Ｓ_rx(ｔ)にできる だけ近づけるようにＦＩＲフィルタのタップ重み｛α｝(これはＣＩＲを近似す るのに用いられる)を決定することである。 送信信号Ｓ_tx(ｔ)及び受信信号Ｓ_rx(ｔ)が分かると、ＣＩＲを次のものから推 定することができる。 １．送信信号の既知のサンプルＴ_k＝Ｓ_tx(ｔ₀＋ｋτ)、−Ｎ_c≦ｋ＜Ｎ(Ｎ＋Ｎ_c サンプルは、全トレーニングシーケンスを表わし、番号付けは、トレーニング シーケンスの開始からＮ_cサンプルの後に発生する受信信号の既知のクリーンサ ンプルから開始するよう構成されることに注意されたい)；及び ２．受信信号の測定サンプルＲ_k＝Ｓ_rx(ｔ₀＋ｋτ)、０≦ｋ＜Ｎ(トレーニン グシーケンスの最初のＮ_c個のサンプルが歪んでいると仮定し、無視する)。 タップ重み｛α_k｝は、次のような簡単な相関により決定される。 このアルゴリズムは、送信信号のノイズ状特性に依存し、従って、その自己相関 関数は、低い時間サイドローブを有していなければならない。しかしながら、ト レーニングシーケンスの「クリーン」な部分の相関特性と、その後の重み付けに 対して部分相関のみが行われることとにより、このアルゴリズムのサイドローブ 性能が制限されることが分かった。 米国特許第５，４７３，６３２号は、チャンネルが実際に出力するものと、推 定予想モデルが出力するものとの間のエラーの平方を最小にすることにより複素 パルス応答（ＣＩＲ）を予想する「最小自乗手法」を使用することによりＣＩＲ を推定する。この手法は、米国特許第５，４７３，６３２号の第６カラムの３５ 行目に記載された式を与え、これは種々の文献から知られているものである。 米国特許第５，４７３，６３２号の手法には問題があると考えられる。基本的 に、米国特許第５，４７３，６３２号は、マトリクスを形成する。問題は、実際 に、ＣＩＲを計算するためにこのマトリクスに対して何が行なわれるかに関連し ている。 米国特許第５，４７３，６３２号の第５カラムの２７ないし６１行目及び請求 項２を参照すれば、マトリクスが反転される前に「係数α」が追加され、そして「 係数α」の追加後に、マトリクスが反転される。実際には、マトリクスの反転は 、ＣＩＲを計算する上で問題を招く。発明の要旨 本発明は、公知技術でＣＩＲ及びＳＩＮＲを決定する際に遭遇した問題を軽減 することに向けられる。 本質的に、本発明は、ネットワーク分析により得られたマトリクスを「対角マ トリクス」へと最初に変換し、そしてそれにより得られたマトリクスを反転する ことにより、公知技術に関連した問題を軽減できるという実現から達成される。 従って、「係数α」は追加されない。 対角マトリクスを用いる効果は、ある固有値（例えば、問題のある小さな固有 値の幾つか）を除去又は無視し、これにより、米国特許第５，４７３，６３２号 においてマトリクスが反転されたときのエラー及びノイズ増加を回避できること である。 対角マトリクスを用いる別の効果は、対角マトリクスの特性に起因するもので 、即ち、対角エレメントのみを反転するだけでマトリクスを反転できることであ る。これは、小さな固有値により生じる問題を露呈し、反転において非常に大き なエレメントを招く。 しかしながら、幾つかの固有値を無視又は除去することにより、図１に示すよ うに、大きなサイドローブの発生という更に別の問題を生じる。 従って、本発明は、大きなサイドローブの発生の問題を軽減することに向けら れた更に別の技術も開示する。この更に別の技術は、上記技術の改良と考えられ 、小さな固有値を破棄する（上記の手法である）のではなく、むしろ、サイドロ ーブが除去されるまで、又はそれらがほとんど影響しなくなるように、反転マト リクスにおいて固有値の逆数に減少する数値を乗算することに基づいている。図面の簡単な説明 以下、添付図面を参照して本発明の好ましい実施形態を詳細に説明する。 図１は、ε＝２００、Ｓ／Ｎ＝０ｄＢを用いる場合に米国特許第５，４７３， ６３２号に開示されたアルゴリズムを用いて推定されるＣＩＲを示す図である。 図２は、ε＝４００を使用し、ノイズがない状態で米国特許第５，４７３，６ ３２号に開示されたアルゴリズムを用いて推定されるＣＩＲを示す図である。 図３は、ＳＶＤを用いて推定されるＣＩＲであって、１０個の最大固有値を含 むＣＩＲを示す図である。 図４は、ＳＶＤを用いて推定されるＣＩＲであって、１６個の最大固有値を含 むＣＩＲを示す図である。 図５は、ＳＶＤを用いて推定されるＣＩＲであって、２５個の最大固有値を含 むＣＩＲを示す図である。 図６は、ＳＶＤを用いて推定されるＣＩＲであって、３６個の最大固有値を含 むＣＩＲを示す図である。 図７は、ＳＶＤを使用すると共に、ｎ₁＝４、ｎ₂＝２５及び追加ノイズなしの 状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図８は、ＳＶＤを使用すると共に、ｎ₁＝４、ｎ₂＝２５及びＳ／Ｎ−０ｄＢの 状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図９は、ＳＶＤを使用すると共に、ｎ₁＝４、ｎ₂＝２０及びノイズなしの状態 で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１０は、ＳＶＤを使用すると共に、ｎ₁＝４、ｎ₂＝２０及びＳ／Ｎ＝０ｄＢ の状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１１は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５及びノイズなしの状 態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１２は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５及びＳ／Ｎ＝０ｄＢ の状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１３は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝６４及びノイズなしの状 態で傾斜関数を使用して推定されるＣＩＲを示す図である(−６０ｄＢまでのｄ Ｂ目盛であることに注意されたい)。 図１４は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝６４及びＳ／Ｎ−４０ｄ Ｂの状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１５は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５及びＳ／Ｎ＝４０ｄ Ｂの状態で傾斜関数を使用して推定されるＣＩＲを示す図である。 図１６は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５、経路１：６．９μ ｓ、０ｄＢ、経路２：１８．５μｓ、−２０ｄＢ、経路３：２９．５μｓ、−１ ０ｄＢ、ノイズなしの状態で傾斜関数を使用して推定されるＣＩＲを示す図 である。 図１７は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５、経路１：６．９μ ｓ、−１０ｄＢＮ経路２：１８．５μｓ、０ｄＢ、経路３：２９．５μｓ、０ｄ Ｂ、経路４：４０．４μｓ、−６ｄＢ、Ｓ／Ｎ−２０ｄＢの状態で傾斜関数を使 用して推定されるＣＩＲを示す図である。 図１８は、ＳＶＤを使用すると共に、ｎ₁＝１、ｎ₂＝２５、経路１：６．９μ ｓ、０ｄＢ、経路２：１２．２μｓ、０ｄＢ、ノイズなしの状態で傾斜関数を使 用して推定されるＣＩＲを示す図である。好ましい実施形態の詳細な説明 タップ重みを推定する１つの方法は、測定信号に最も厳密に一致する受信信号 を予想する重みを決定することである。 既知のサンプルＴ_k＝Ｓ_tx(ｔ₀＋ｋτ)、−Ｍ≦ｋ＜Ｎ（Ｍは、仮定したＣＩＲ の範囲から「ダーティ」サンプルの数を引いたものをもつ送信信号Ｓ_tx(ｔ)と、測 定サンプルＲ_k＝Ｓ_rx(ｔ₀＋ｋτ)、０≦ｋ＜Ｎをもつ受信信号Ｓ_rx(ｔ)とが与え られる。 所望のＦＩＲタップ重み（α_k）は、測定された受信信号と、受信信号の推定 値（即ちチャンネルを近似するＦＩＲを通過した既知の送信信号）との間の累積 平方エラーを最小にすることを必要とすることにより推定され、即ち次のものを 最小にする。 これは、同時方程式の次の系へ変換することができる。 又はマトリクス表示法を使用する。 Ａｃ＝ＴＲ （２．４） 但し、マトリクスＡは、次の式で与えられるＭ×Ｍ相関マトリクスである。 Ｔは、共役の遅延信号のＭ×Ｎマトリクスである。 ｔ_i,j＝Ｔ_i,j ０≦ｉ≦Ｍ、０≦ｊ≦Ｎ （２．６） ベクトルｃは、必要とされるＣＩＲ(α_k)である。 簡単な相関は、次の式に対応することに注意されたい。 ｃ＝ＴＲ （２．７） マトリクス式（２．４）は、次のような慣例的な解を有する。 ｃ＝(Ａ^-1Ｔ)Ｒ （２．８） これは、単純なマトリクス乗算を用いて受信信号のサンプルからＣＩＲを決定で きることを示す。 米国特許第５，４７３，６３２号を用いた最小自乗推定 式（２．８）は、米国特許第５，４７３，６３２号の第６カラム、３５行目に 開示された式と同じである。 この式を使用する場合には、式が一般的に不完全な条件にあるという点で式を 直接解くときに複雑なものとなる。実際に、ＳＢにおけるＧＳＭトレーニングシ ーケンスについては、マトリクスの固有値が１０¹⁰以上の範囲にわたって変化す ることが分かっている。これは、マトリクスの反転が丸めエラーの影響を非常に 受け易く、そしてその結果が測定信号のノイズの影響を非常に受け易くすると考 えられる。 従来の技術を適用してこのマトリクスを反転するか又は式を解くことにより、 ＩＥＥＥ倍精度演算を用いたときでも、有用な結果が得られないと考えられる。 米国特許第５，４７３，６３２号は、マトリクスの反転を安定化しそしてノイ ズに対して適度に弾力性のある結果を発生するアルゴリズムを開示している。米 国特許第５，４７３，６３２号は、次の式によりＣＩＲを推定する。 ｃ＝((Ａ＋εＩ)^-1Ｔ)Ｒ （２．９） 但し、Ｉは、認識マトリクスであり、そしてεは、定数（ノイズ項と称する）で ある。この式は、米国特許第５，４７３，６３２号の第６カラムの４４行目に記 載された式に等しい。このアルゴリズムを用いて得られる結果は、次の通りであ る。 このアルゴリズムの性能は、変数εの性能に厳密に依存する。大きさ１９２ の対角エレメントを有するマトリクスＡについては、異なる値のεを用いて、信 号経路より成るＣＩＲを決定する米国特許第５，４７３，６３２号のアルゴリズ ムの性能が次の図に示されている。 例えば、図１及び２を参照すれば、分解能（メインローブの巾）とノイズ性能 との間に妥協があることが明らかである。約−１８ｄＢのサイドローブレベルが このアルゴリズムの特性であると考えられる。 ε≡１００を使用することにより最良の比較し得る性能が得られ、これに基づ いて、米国特許第５，４７３，６３２号との比較を行う。 ＳＶＤを用いた最小自乗推定 本発明においては、特異値分解（ＳＶＤ）を不完全な条件のマトリクスに適用 して最小自乗推定を実行するように選択される。このように、対角マトリクスヘ と反転されるべきマトリクスを変換し、次いで、その新たに形成された対角マト リクスを反転する。 ＳＶＤを使用し、マトリクスＡ（正の有限）を次の式で表すことができる。 Ａ＝ＵＳＶ^t （２．１０） 但し、Ｕ及びＶは直交であり、そしてＳは対角である。実際に、Ｓの対角エレメ ントは、Ａの固有値（全て正）であり、これらは、減少する順に配置される。 ｓ_1.1≧ｓ_2.2≧…≧ｓ_n.n （２．１１） マトリクスが特異なものである場合には、これらの幾つかがゼロである。従って 、Ａの反転は、次のようになる。 Ａ^-1＝ＶＳ^-1Ｕ^t （２．１２） 但し、Ｓ^-1は、Ｓの（対角）エレメントを反転することにより形成される。不完 全な条件のマトリクスの作用はここでは明確であり、（ほぼゼロの）固有値が反 転され、そして反転に対して非常に大きな貢献を与える。これらはノイズを増幅 する。 従来のＳＶＤ技術は、固有値の幾つかだけを反転し、あるポイントで停止する 。１／λ₁をＳの反転に入れるのではなく、本発明はゼロを使用する。これら技 術及び他の技術の性能は、以下で評価する。 式（２．８）のマトリクスの係数は、既知のトレーニングシーケンスの純粋な 関数であり、従って、それらは前もって計算することができる。このマトリクス は、一度反転することができ、従って、解はいつでもマトリクス乗算を必要とす る。従って、１組の信号サンプルからＣＩＲ推定値を得るためには、Ｍ＊Ｎの複 素乗算と、Ｍ＊（Ｎ−１）の複素加算が必要となり、全部で６＊Ｍ＊Ｎ＋２＊Ｍ ＊（Ｎ−１）の演算が必要となる。 断固とした最小自乗問題にＳＶＤを適用するには、マトリクスの擬似反転にお いてＮ個の最大固有値からの貢献を含ませるだけでよい。ここにモデリングされ るケースでは、マトリクスが次元６４を有し、従って、Ｎは１ないし６４の範囲 である。以下のグラフは、Ｎ＝１０、１６、２５、３６の場合に対するＳＶＤア ルゴリズムの性能を示す。固有値の数を増加する作用は、推定装置の分解能を高 めることであり、即ちメインローブの巾を減少することである。 図示されていないが、シミユレーションを行い、１０ないし２０個の固有値の 範囲内でノイズ性能が著しく低下することが示された。 このＳＶＤ解決策は、擬似反転のランクが最小であるので、優れたノイズ性能 を形成する。結果について最も厄介なことは、約−１３ｄＢの持続的な時間サイ ドローブがあることである。フーリエ理論と同様に、裁断された時間波形が−１ ３ｄＢの周波数サイドローブを与える場合には、ある数の後に単にカットオフす るのではなく、次々の固有値からの貢献が下方傾斜することが、時間サイドロー ブに対する制御を与える上で非常に有用であることが分かった。 バートレット、ブラックマン・ハリス、ポイゾン、レイツ、ハミング、カイザ ー・ベッセル、ランゾス、タキー及び他の良く知られたウインドウ関数のような 多数の傾斜関数を使用することができる。この実施形態で選択される関数は、ハ ミングウインドウから導出されるが、以下に述べる累乗コサイン傾斜が最も有用 であることが分かっている。 但し、ｗ_iは、ｉ番目の固有値からの貢献に適用される重みである。例えば、図 ７は、ｎ₁＝４、ｎ₂＝２５の作用を示している。 −１３ｄＢのサイドローブは、約−２２ｄＢに減少されたことに注意されたい (上記米国特許第５，４７３，６３２号の性能より著しく良好である)。メインロ ーブの巾は、低いサイドローブをもつε＝１００の米国特許第５，４７３，６３ ２号のアルゴリズムより狭い。Ｓ／Ｎ比が０ｄＢの状態の性能が図８に示されて いる。 この性能は、ε＝１００の米国特許第５，４７３，６３２号の解決策で得られ るものに類似しているが、メインローブはより狭く、これがこのノイズ性能を低 下することに注意されたい。 ノイズ性能は、ｎ₁＝４、ｎ₂＝２０の場合に図９及び１０に示すように、ロー ブの巾を犠牲として改善することができる。 満足な解決策を見出すために著しい試みが行なわれた。その目標は、次の通り であった。 − 米国特許第５，４７３，６３２号の性能に勝る顕著な縁を与えるための２ ０ｄＢ以上のダイナミックレンジ(好ましくは２５ｄＢ)、及び − Ｓ／Ｎ−０ｄＢにおける相関と同様のノイズ性能(これは、ノイズの貢献 を不必要に増大する解決策を望まないことを意味する)。 以上のことから決定された結果は、ｎ₁＝１、ｎ₂＝２５で得られる重み係数を使 用し、そしてその性能は、図１１及び１２に示されている。 更に複雑なＣＩＲに対するこのアルゴリズムの性能を簡単に説明する。図１３ の例に示すように、傾斜アルゴリズムがどれほど大きなダイナミックレンジを与 えることができるか調べることが重要である。 この例は、ノイズの影響を非常に受け易い。図１４は、Ｓ／Ｎ＝４０ｄＢでの 性能を示す。この僅かな量のノイズでも著しく増大されることに注意されたい。 この同じ量のノイズは、上記のように決定されそして図１５に示された好ましい 結果に対して目立ったノイズを発生しない。 図１６及び１７は、更に複雑なＣＩＲに対する結果を示す。 この実施形態は、アルゴリズムがピーク間を弁別できるようにすることが分か っており、そして図１８は、５μＳ以上の時間的分離を伴う密接したピーク分 離を示している。 又、上述したマトリクス反転は、干渉打消しに関連してチャンネルインパルス 応答ベクトルｈを推定するときにも必要となる。 Φｈ＝Ψ （２．１４） 但し、Φは、局部的に発生される信号間のクロス相関マトリクスであり、そして Ψは、局部的に発生された信号と受信信号との間のクロス相関マトリクスである 。この状態においても、上述した原理及び方法が等しく適用される。 同様に、干渉制限された環境においては、通常、多数のセンサ（アンテナアレ ー）から受け取られた信号の最適な組合せを用いてＣＤＭＡが適用される。最適 な組合わせは、信号対干渉及び雑音比（ＳＩＮＲ）を最大にする。ＳＩＮＲを最 大にするためのアンテナアレーの重みベクトルは、例えば、次の式から得られる 。 ｗ＝αＲ^-1Ｕ_d ^* （２．１５） 但し、ｗはアンテナ素子の重みベクトルであり、αは定数であり、Ｒは受け取っ た干渉及びノイズ相関マトリクスであり、そしてＵ_d ^*は所望の信号ベクトルの共 役ベクトルである。重みベクトルを得る効率的な方法は、ＣＩＲを最初に計算す ることであることが分かった。ある最適な組合わせ技術は、最初にＣＩＲ(及び ノイズ共変数マトリクスＱ)を計算し、これを多次元ＭＬＳＥに使用して、アン テナアレーからの信号が合成される。通常、移動無線チャンネルにおいてはマト リクスの反転が問題である。上記の実施形態は、マトリクス反転に直接適用する ことができる。重みベクトルを得るか又はＱの反転を使用するために、信号が最 適に合成される。 要約すれば、本発明により種々の変更や修正が包含されるが、ここに示した好 ましい実施形態は、式（２．８）に規定されたマトリクス乗算を用いてＣＩＲを 発生し、そしてＳＶＤを使用して、マトリクスＡの反転を決定するものである。 反転に含まれる固有値は、ｎ₁＝１、ｎ₂＝２５の状態で式（２．１３）に基づい て傾斜化（重み付け）される。

───────────────────────────────────────────────────── フロントページの続き (81)指定国 ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ， ＤＫ，ＥＳ，ＦＩ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，Ｌ Ｕ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＯＡ(ＢＦ，ＢＪ，ＣＦ ，ＣＧ，ＣＩ，ＣＭ，ＧＡ，ＧＮ，ＭＬ，ＭＲ，ＮＥ， ＳＮ，ＴＤ，ＴＧ)，ＡＰ(ＧＨ，ＧＭ，ＫＥ，ＬＳ，Ｍ Ｗ，ＳＤ，ＳＺ，ＵＧ，ＺＷ)，ＥＡ(ＡＭ，ＡＺ，ＢＹ ，ＫＧ，ＫＺ，ＭＤ，ＲＵ，ＴＪ，ＴＭ)，ＡＬ，ＡＭ ，ＡＴ，ＡＵ，ＡＺ，ＢＡ，ＢＢ，ＢＧ，ＢＲ，ＢＹ， ＣＡ，ＣＨ，ＣＮ，ＣＵ，ＣＺ，ＤＥ，ＤＫ，ＥＥ，Ｅ Ｓ，ＦＩ，ＧＢ，ＧＥ，ＧＨ，ＧＭ，ＧＷ，ＨＵ，ＩＤ ，ＩＬ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ， ＬＣ，ＬＫ，ＬＲ，ＬＳ，ＬＴ，ＬＵ，ＬＶ，ＭＤ，Ｍ Ｇ，ＭＫ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，ＮＺ，ＰＬ，ＰＴ ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＳＬ， ＴＪ，ＴＭ，ＴＲ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，Ｖ Ｎ，ＹＵ，ＺＷ

Claims (1)

1. 【特許請求の範囲】 １．通信システムの推定予想モデルにより使用するために通信システムのチャン ネルインパルス応答ＣＩＲを推定する方法において、 測定された受信信号と、その受信信号の所定の推定値との間のエラーの平方を 最小にすることによりＣＩＲを予想し、上記推定値は、無線チャンネルの近似 に基づくものであり、 上記ＣＩＲの第１のマトリクス式を形成し、 上記第１マトリクスを対角マトリクスに変換し、そして 上記対角マトリクスを反転する、 という段階を含むことを特徴とする方法。 ２．上記対角マトリクスは、特異値分解ＳＶＤを第１マトリクスに適用すること により得られる請求項１に記載の方法。 ３．上記反転段階に続いて、その反転された対角マトリクスに生じる過剰な固有 値の貢献を下方に傾斜化する段階を更に含む請求項１又は２に記載の方法。 ４．ウインドウ関数を使用して、過剰な固有値を下方傾斜化する請求項３に記載 の方法。 ５．バートレット、ブラックマン、ハリス、ガウス、ポイゾン、レイツ、ハミン グ(累乗コサイン)、カイザー・ベッセル、ランゾス、タキー、又は他の同様の 既知のウインドウ関数のいずれかを使用して、上記過剰な固有値を下方傾斜化 させる請求項４に記載の方法。 ６．ＧＳＭネットワーク、干渉打消し又はＣＤＭＡ受信器に使用される請求項１ ないし５のいずれかに記載の方法。