(57)【要約】  本発明は、無線ネットワーク用の受信機に関する。受信機は、受信ユニットと、チャネル推定ユニット（１８）とを備える。受信ユニットは、各々がサイクリックプレフィックスを有する複数のシンボルを形成するいくつかの時間領域信号を受信するように構成される。チャネル推定ユニット（１８）は、受信ユニットによって受信された信号内の既知のパイロットに基づいて、ユーザ（ｕ）に対して時間領域内のチャネル推定値（＾ｈｏｓ）を決定し、ユーザ（ｕ）に関連付けられたサブキャリアの数（Ｎｃ）に関連する基本期間を有するＤＦＴ行列（Ｖ）を使用してユーザ（ｕ）に対してチャネル推定値（＾ｈｏｓ）を決定するように構成され、チャネル推定値（＾ｈｏｓ）は、Ｌ個のタップ（ただし、Ｌはユーザに関連付けられたサブキャリアの数（Ｎｃ）より小さい）の長さを有する。また本発明は、無線チャネルのチャネル推定値（＾ｈｏｓ）を決定する方法に関する。
2010526511.tif 000002
【特許請求の範囲】
【請求項１】
無線ネットワーク用の受信機（３０）であって、当該受信機は、
複数のシンボルを形成するいくつかの時間領域信号を受信するように構成された受信ユニット（３０）であって、各々のシンボルがサイクリックプレフィックスを有する、受信ユニット（３０）と、
前記受信ユニットによって受信された前記信号内の既知のパイロットに基づいてユーザ（ｕ）に対して前記時間領域内のチャネル推定値（＾ｈｏｓ）を決定するように構成されたチャネル推定ユニット（１８）と、
を備え、
前記チャネル推定ユニット（１８）が、ユーザ（ｕ）に関連付けられたサブキャリアの数（Ｎｃ）に関連する基本期間を有するＤＦＴ行列（Ｖ）を使用して前記ユーザ（ｕ）に対してチャネル推定値（＾ｈｏｓ）を決定するように構成され、前記チャネル推定値（＾ｈｏｓ）が、Ｌ個のタップの長さを有し、Ｌは、前記ユーザに関連付けられたサブキャリアの数（Ｎｃ）より少ないことを特徴とする受信機。
【請求項２】
前記チャネル推定ユニット（１８）が、
前記チャネル推定値を
【数１】
2010526511.tif 000045
として推定するように構成されたチャネル推定部（２１）を備え、Ｖ＊Ｎｃ・ｋが、基本期間Ｎｃ・ｋを有するＤＦＴ行列の共役であり、Ｄ＊ｕが、前記ユーザ（ｕ）に対応する周波数領域内のパイロットシンボルを含む行列であるＤｕの共役であり、Ｒｕが、前記ユーザ（ｕ）の周波数領域内の受信サンプルを示し、ｋが定数であることを特徴とする、請求項１に記載の受信機。
【請求項３】
前記基本期間（Ｎｃ・ｋ）が、前記ユーザに関連付けられたサブキャリアの数Ｎｃをオーバーサンプリング係数ｋで乗算した値に等しいことを特徴とする、請求項１又は２に記載の受信機。
【請求項４】
前記オーバーサンプリング係数ｋが、１．０〜２．５の範囲内にあることを特徴とする、請求項３に記載の受信機。
【請求項５】
前記オーバーサンプリング係数ｋが、２であることを特徴とする、請求項４に記載の受信機。
【請求項６】
前記ＤＦＴ行列（Ｖ）の行数が、検索ウィンドウの長さＬに等しいことを特徴とする、前記請求項のいずれか１項に記載の受信機。
【請求項７】
前記チャネル推定ユニット（１８）が、チャネルタップ低減推定値＾ｈｏｓ＿ｒｅｄを形成するように、前記チャネル推定値
【数２】
2010526511.tif 000046
内の行のサブセットを選択するように構成されたチャネルタップ低減部（２２）を備えることを特徴とする、請求項６に記載の受信機。
【請求項８】
前記チャネルタップ低減部（２２）が、前記チャネル推定値（＾ｈｏｓ）に対して前記検索ウィンドウ内で第１の終点（＾ｍＲ）と第２の終点（＾ｍＬ）とを決定し、前記チャネルタップ低減推定値（＾ｈｏｓ＿ｒｅｄ）が、前記第１及び第２の終点の間の行のみを含む行のサブセットを選択するように構成されることを特徴とする、請求項７に記載の受信機。
【請求項９】
前記チャネルタップ低減部（２２）が、
第１の方向の前記チャネル推定値（＾ｈｏｓ）のいくつかの可能な長さの各々の決定値を決定し、最も好適な決定値に関連付けられた長さを選択することにより、前記第１の終点（＾ｍＲ）を決定するように構成された第１の適応チャネル推定器（２５）と、
第２の方向の前記チャネル推定値（＾ｈｏｓ）のいくつかの可能な長さの各々の決定値を決定し、最も好適な決定値に関連付けられた長さを選択することにより、前記第２の終点（＾ｍＬ）を決定するように構成された第２の適応チャネル推定器（２６）と、
を備えることを特徴とする、請求項８に記載の受信機。
【請求項１０】
前記第１の適応チャネル推定器（２５）が、前記チャネル推定値（＾ｈｏｓ）の予測されたタイミングに関連付けられた開始値（Ｎｓｔａｒｔ）に可能な長さを関連付けるように構成されることを特徴とする、請求項９に記載の受信機。
【請求項１１】
前記第２の適応チャネル推定器（２６）が、前記第１の終点（＾ｍＲ）に可能な長さを関連付けるように構成されることを特徴とする、請求項１０に記載の受信機。
【請求項１２】
前記チャネル推定値（＾ｈｏｓ）の決定されたタイミングと前記開始値（Ｎｓｔａｒｔ）によって定義される前記チャネル推定値（＾ｈｏｓ）の前記予測されたタイミングとの間のタイミングミスアラインメントを決定するように構成されたタイミングミスアラインメント決定ユニット（２９）を備えることを特徴とする、請求項１０に記載の受信機。
【請求項１３】
前記タイミングミスアラインメント決定ユニット（２９）が、前記開始サンプル値（Ｎｓｔａｒｔ）と前記第１及び／又は前記第２の終点との間の予測された関係に基づいて前記タイミングミスアラインメントを決定するように構成されることを特徴とする、請求項１２に記載の受信機。
【請求項１４】
ＵＭＴＳ規格のリリース８に準拠して動作するように構成されることを特徴とする前記請求項のいずれか１項に記載の受信機。
【請求項１５】
前記請求項のいずれか１項に記載の受信機を備えることを特徴とする基地局（２）。
【請求項１６】
前記タイミングミスアラインメントユニット（２９）に動作可能に接続され、前記ミスアラインメントに関連する情報を前記関連するユーザ（ｕ）へ送信するように構成される送信ユニット（３１）を備えることを特徴とする、請求項１２又は１３を参照する請求項１５に記載の基地局（２）。
【請求項１７】
請求項１６に記載の基地局を備えることを特徴とする無線ネットワーク。
【請求項１８】
請求項１から１４のいずれか１項に記載の受信機を備えることを特徴とする端末（１）。
【請求項１９】
無線チャネルのチャネル推定値（＾ｈ）を決定する方法であって、前記方法が、
複数のシンボルを形成するいくつかの時間領域信号（ｙｎ）を受信することであって、各々のシンボルがサイクリックプレフィックスを有する、ことと、
前記受信信号内の既知のパイロットに基づいて、ユーザ（ｕ）に対して前記時間領域内の前記チャネル推定値（＾ｈ）を決定することと、
を含み、
ユーザ（ｕ）に対して前記チャネル推定値（＾ｈ）を決定するステップが、前記ユーザ（ｕ）に関連付けられたサブキャリアの数（Ｎｃ）に関連する基本期間を有するＤＦＴ行列（Ｖ）を使用することを含むことを特徴とする方法。
【請求項２０】
ユーザ（ｕ）に対して前記チャネル推定値（＾ｈ）を決定することが、
基本期間Ｎｃ・ｋを有する前記ＤＦＴ行列（Ｖ）の共役Ｖ＊Ｎｃ・ｋを決定するステップと、
前記ユーザ（ｕ）に対応する周波数領域内のパイロットシンボルを含むパイロット行列Ｄｕの共役Ｄ＊ｕを決定するステップと、
前記ユーザ（ｕ）に対して前記プレフィックスを除去し、周波数領域内の受信サンプルのベクトル（Ｒｕ）を形成するように前記受信信号（ｙｎ）を処理するステップと、
前記チャネル推定値を
【数３】
2010526511.tif 000047
として推定するステップと、
を含むことを特徴とする、請求項１９に記載の方法。
【請求項２１】
チャネルタップ低減チャネル推定値（＾ｈｏｓ＿ｒｅｄ）を形成するように、前記チャネル推定値
【数４】
2010526511.tif 000048
行列の行のサブセットが選択されることを特徴とする、請求項１９又は２０に記載の方法。
【請求項２２】
行のサブセットの前記選択ステップが、前記チャネル推定値＾ｈｏｓに対して第１の終点（＾ｍＲ）及び第２の終点（＾ｍＬ）を決定するステップと、前記チャネルタップ低減チャネル推定値（＾ｈｏｓ＿ｒｅｄ）が、前記第１及び第２の終点の間の行しか含まないように、行のサブセットを選択するステップとを含むことを特徴とする、請求項２１に記載の方法。
【請求項２３】
前記第１の終点（＾ｍＲ）の決定が、開始サンプル値（Ｎｓｔａｒｔ）から始まる第１の方向の前記チャネル推定値＾ｈｏｓのいくつかの可能な長さの各々に対して決定値を決定することと、最も好適な決定値に関連付けられた長さを選択することとを含み、
前記第２の終点（＾ｍＬ）の決定が、前記第１の方向とは逆の第２の方向の前記チャネル推定値＾ｈｏｓのいくつかの可能な長さの各々に対して決定値を決定することと、最も好適な決定値に関連付けられた長さを選択することとを含む、
ことを特徴とする、請求項２２に記載の方法。
【請求項２４】
前記チャネルタップ低減チャネル推定値の予測されたタイミングと前記決定されたチャネルタップ低減チャネル推定値（＾ｈｏｓ＿ｒｅｄ）との間のタイミングミスアラインメントが決定されることを特徴とする、請求項２２に記載の方法。
【請求項２５】
前記タイミングミスアラインメントに関連する情報が、前記ユーザ（ｕ）へ送信されることを特徴とする、請求項２４に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
技術分野
本発明は、受信ユニットを備える無線ネットワーク用の受信機及びユーザに対してチャネル推定値を決定するように構成されたチャネル推定ユニットに関する。
【０００２】
本発明は、また、受信信号内の既知のパイロットに基づいてユーザに対して無線チャネルのチャネル推定値を決定する方法に関する。
【背景技術】
【０００３】
背景
現代のセルラー無線システムでは、無線ネットワークは、ネットワーク内の端末の挙動を厳重に制御している。アップリンク送信パラメータ、すなわち、制御対象の端末（ＷＣＤＭＡ用語ではＵＥとも呼ばれる）から基地局（ＷＣＤＭＡ用語ではノードＢとも呼ばれる）への通信チャネルを介した送信のためのパラメータは、例えば、周波数、タイミング及び電力である。アップリンク送信パラメータは、基地局から端末へのダウンリンク制御信号方式を介して調整される。
【０００４】
端末の制御を可能にするためには、アップリンク信号の測定が必要である。タイミングミスアラインメント（timing misalignment）の決定は、基地局が受信信号のタイミングを基地局内の基準クロックからの基準クロック信号と比較することを必要とする。タイミングミスアラインメントは、基地局内の基準クロックと端末内のクロックとの間の予測できない伝搬遅延及び相互ドリフトによって引き起こされる。同様に、送信電力の決定は、所定のしきい値に対する基地局内で受信された信号の信号電力の決定を必要とする。受信電力は、端末の送信電力だけでなく伝搬時の信号の減衰にも左右される。信号の減衰は、経時的に変化する。
【０００５】
端末と基地局との間の通信チャネル上のアップリンク無線通信は、マルチパス伝搬、フェージング、周波数エラー、往復伝搬時間などの問題を抱える。このアップリンク通信チャネルは、多くの場合、エアインタフェースと呼ばれる。上記エアインタフェースを介して送信される情報内にビットエラー及びブロックエラーが発生することは珍しくない。基地局のアップリンク受信機は、ビットエラー及びブロックエラー率を低減するように構成されている。基地局受信機の１つの不可欠な部分は、アップリンク通信チャネルを推定するように構成されたチャネル推定ユニットである。
【０００６】
ＤＦＴ（ＦＦＴとして実施される）が時間領域信号上で基地局受信機によって使用される時はいつでも、上記時間領域信号は、周期的で周期時間がＤＦＴの長さに等しいと仮定する。このことを分散チャネルで送信される信号について確保するために、サイクリックプレフィックス（cyclic prefix）を使用することができる。サイクリックプレフィックスは、アップリンク通信チャネルの遅延拡散より大きいかそれに等しい長さを有していなければならない。また、いかなるタイミングミスアラインメントもサイクリックプレフィックスより低い解像度を有していなければならない。それ故、タイミングミスアラインメントの決定アルゴリズムは、好ましくはこのサイクリックプレフィックスの長さより低い解像度を有する。
【発明の概要】
【０００７】
概要
先行技術のチャネル推定ユニットの１つの問題は、その中で使用されるアルゴリズムが高い計算の複雑さを有するか、又は大量の事前計算された値を必要とするということである。本発明の１つの目的はこの問題を低減することである。
【０００８】
本発明の一実施形態によれば、これは、無線ネットワーク用の受信機によって達成されている。この受信機は、受信ユニットと、チャネル推定ユニットとを備える。受信ユニットは、複数のシンボルを形成するいくつかの時間領域信号を受信するように構成され、各々のシンボルはサイクリックプレフィックスを有している。チャネル推定ユニットは、受信ユニットによって受信された信号内の既知のパイロットに基づいてユーザに対して時間領域内のチャネル推定値を決定し、ユーザに関連付けられたサブキャリアの数に関連する基本期間を備えたＤＦＴ行列を使用して上記ユーザに対してチャネル推定値を決定するように構成され、上記チャネル推定値は、Ｌ個のタップの長さを有し、Ｌは上記ユーザに関連付けられたサブキャリアの数より小さい。
【０００９】
本発明の利点の１つは、チャネル推定ユニットは、実施が容易であるように、すなわち、計算回数及び予め計算された値の数の両方が小さく維持されるように設計されているという点である。
【００１０】
一実施例では、受信機は基地局内に提供されるが、端末又はその他の装置内に実施することもできる。
【００１１】
本発明の一実施形態によれば、チャネル推定ユニットは、チャネル推定値を
【数１】
2010526511.tif 000003
として推定するように構成されたチャネル推定部を備える。一実施例では、基本期間（Ｎｃ・ｋ）は、ユーザに関連付けられたサブキャリアの数にオーバーサンプリング係数ｋを乗算して得た値に等しい。オーバーサンプリング係数ｋが＞１の場合、オーバーサンプリングされたチャネルモデルが提供され、これは、さらなる性能の向上を意味する。
【００１２】
本発明の一実施形態によれば、ＤＦＴ行列（Ｖ）内の行数は、検索ウィンドウの長さＬに等しい。しかし、この数は、チャネルタップ低減推定値＾ｈｏｓ＿ｒｅｄ（以下、文字の前に記号＾が記載されている場合は、本来当該文字の上に記載される記号＾を、便宜上当該文字の前にずらして記載したものとする。）を形成するようにチャネル推定値
【数２】
2010526511.tif 000004
内の行のサブセットを選択するように構成されたチャネルタップ低減部によってさらに低減することができる。一実施例では、チャネルタップ低減部は、チャネル推定値＾ｈｏｓに対して検索ウィンドウ内の第１の終点＾ｍＲ（以下、＾ｍＲと示す）及び第２の終点＾ｍＬ（以下、＾ｍＬと示す）を決定し、チャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄが第１及び第２の終点の間の行しか含まないように行のサブセットを選択するように構成されている。従って、チャネルタップ低減部は、一実施例では、第１の方向のチャネル推定値＾ｈｏｓのいくつかの可能な長さの各々の決定値を決定し、最も好適な決定値に関連付けられた長さを選択することで第１の終点＾ｍＲを決定するように構成された第１の適応チャネル推定器と、第２の方向のチャネル推定値のいくつかの可能な長さの各々の決定値を決定し、最も好適な決定値に関連付けられた長さを選択することで第２の終点＾ｍＬを決定するように構成された第２の適応チャネル推定器とを備える。
【００１３】
上記チャネルタップ低減ユニットによって、タイミングミスアラインメントの決定が可能になる。従って、受信機は、チャネル推定値＾ｈｏｓの決定されたタイミングと開始値（Ｎｓｔａｒｔ）によって定義されるチャネル推定値＾ｈｏｓの予測されたタイミングとの間のタイミングミスアラインメントを決定するように構成されたタイミングミスアラインメント決定ユニットを備える。タイミングミスアラインメント決定ユニットは、一実施例では開始サンプル値（Ｎｓｔａｒｔ）と第１及び／又は第２の終点との間の予測された関係に基づいてタイミングミスアラインメントを決定するように構成されている。上記タイミングミスアラインメント決定ユニットを使用することで、タイミングミスアラインメントは高精度に決定することができる。
【００１４】
一実施形態では、受信機は、ＵＭＴＳ規格のリリース８に準拠して動作するように構成されている。
【００１５】
本発明は、添付の特許請求の範囲のいずれかによる無線チャネルのチャネル推定値を決定する方法をさらに備える。本発明は、また、特許請求の範囲に規定されたチャネルタップ除去手順に基づくタイミングミスアラインメント決定の方法に関する。
【００１６】
以下に添付の図面を参照しながら本発明を例示的に説明する。
【図面の簡単な説明】
【００１７】
図面の簡単な説明
【図１】図１は、無線ネットワーク内の端末及び基地局の一実施例を概略的に示す。
【図２】図２は、図１の端末の送信部の一実施例のブロック図である。
【図３】図３は、図１の基地局の受信機の一実施例のブロック図である。
【図４】図４は、図３の受信機の受信部の一実施例のブロック図である。
【図５】図５は、図３の受信機内のチャネル推定ユニットの一実施例のブロック図である。
【図６】図６は、本発明の第１の例によるチャネル推定値の周波数応答を示す概略グラフである。
【図７】図７は、本発明の第２の例によるチャネル推定値の周波数応答を示す概略グラフである。
【図８】図８は、図５のチャネル推定ユニット内のチャネルタップ低減部が用いる原理の一実施例を概略的に示す。
【図９】図９は、図５のチャネル推定ユニット内のチャネルタップ低減部のブロック図である。
【図１０】図１０は、無線チャネルのチャネル推定値を決定する方法の一実施例の流れ図である。
【００１８】
詳細な説明
図１を参照すると、セルラー無線システムは、いくつかの端末１（図１には１つだけ示す）と基地局２とを備える。無線ネットワークは、ネットワーク内の端末１の挙動を厳重に制御している。端末１は、ＷＣＤＭＡ用語ではユーザ装置（ＵＥ）と呼ばれる。基地局は、ＷＣＤＭＡ用語ではノードＢと呼ばれる。端末は、端末１から基地局２へ通信チャネル３上で送信可能なように信号を処理するように構成された送信部１２を備える。送信部について以下に詳細に説明する。アップリンク送信パラメータ、すなわち、無線ネットワークによって制御される端末１から基地局２への通信チャネル３を介した送信のためのパラメータは、例えば、周波数、タイミング及び電力である。
【００１９】
基地局は、アップリンク送信を受信するように構成された受信機３０とダウンリンク送信を送信するように構成された送信機３１とを備える。アップリンク送信パラメータは、ダウンリンク通信チャネル４を介した基地局から端末へのダウンリンク制御信号方式によって調整される。端末は、上記ダウンリンク送信を受信するように構成された受信部１３をさらに備える。
【００２０】
端末１の制御を可能にするためには、アップリンク信号の測定が必要である。タイミングミスアラインメントの決定は、基地局が受信信号のタイミングを基地局内の基準クロック５からの基準クロック信号と比較することを必要とする。タイミングミスアラインメントは、基地局２内の基準クロック５と端末１内の内部クロック６との間の予測できない伝搬遅延及び相互ドリフトによって引き起こされる。同様に、送信電力の決定は、所定のしきい値に対する基地局２内で受信された信号の信号電力の決定を必要とする。受信電力は、端末１の送信電力だけでなく伝搬時の信号の減衰にも左右される。信号の減衰は、経時的に変化する。
【００２１】
図２を参照すると、端末１の送信部１２は、時間領域内の信号ｃ１（ｔ）を離散周波数領域信号Ｃ１に変換するように構成されたＤＦＴ変換器７を備える。添え字「１」は、信号が端末に関連付けられたユーザ１に関連することを示す。周波数領域信号は、いくつかのサブキャリア信号Ｃ１（ｆ１），Ｃ１（ｆ２），Ｃ１（ｆ３），．．．，Ｃ１（ｆｎ）に分割される。ＤＦＴ変換器は、特に周波数信号を１２〜１２００のサブキャリアに分割するように構成されている。使用されるサブキャリアの数は、用途と端末の容量によって変わる。サブキャリアの数の選択については本明細書では説明しない。一実施例によれば、ＤＦＴ変換器７は高速フーリエ変換（ＦＦＴ）として実施される。
【００２２】
離散周波数領域信号Ｃ１は、マッピングユニット８に供給される。マッピングユニット８は、１ユーザから離散周波数領域信号Ｃ１を受信するように構成されている。離散周波数領域信号Ｃ１（ｆ１），Ｃ１（ｆ２），Ｃ１（ｆ３），．．．，Ｃ１（ｆｎ）は、上記ユーザに関連付けられたいくつかのサブキャリアに関連付けられている。基地局２は、どのサブキャリアが上記ユーザに関連付けられているかについての情報を有する。本明細書に記載するユーザの送信部１２では、他のユーザ（２〜ｍ）の他のすべてのサブキャリアがゼロにセットされる。サブキャリアのこのマッピングは、ユーザ間でサブキャリアの重なりがないように実行される。マッピングユニット８は、マッピング後の離散周波数領域信号をＩＤＦＴ変換器９へ供給するように構成されている。ＩＤＦＴ変換器９は、受信したマッピング後の周波数領域信号に対して逆離散変換動作を実行し、時間領域信号ｃｎ（ｔ）を出力するように構成されている。一実施例によれば、ＩＤＦＴ変換器９は逆高速フーリエ変換（ＩＦＦＴ）として実施される。
【００２３】
ＩＤＦＴ変換器９から受信した時間領域信号ｃｎ（ｔ）は、サイクリックプレフィックスユニット１０に供給される。サイクリックプレフィックスユニット１０は、時間領域信号ｃｎ（ｔ）にサイクリックプレフィックスを追加するように構成されている。言い換えれば、サイクリックプレフィックスユニットは、サイクリックプレフィックスｓｎ（ｔ）を有する出力信号を提供すべく、時間領域信号ｃｎ（ｔ）の先頭に時間領域信号ｃｎ（ｔ）の末尾から得た部分を追加するように構成されている。サイクリックプレフィックスの長さは、アップリンク通信チャネル３の遅延拡散より大きいか、又はそれに等しくなるように選択されるという特徴がある。サイクリックプレフィックスユニットから出力されたサイクリックプレフィックスｓｎ（ｔ）を有する信号は、無線送信機１１に供給されて送信される。
【００２４】
図３を参照すると、基地局２の受信機は、受信部とチャネル推定ユニット１８とを備える。受信部１４は、ユーザごとに受信信号ｙｎ（ｔ）内のデータを抽出するように構成されている。端末１と基地局２との間の通信チャネル３上のアップリンク無線通信は、マルチパス伝搬、フェージング、周波数エラー、往復伝搬時間などの問題を抱える。アップリンク通信チャネル３は、多くの場合、エアインタフェースと呼ばれる。上記エアインタフェースを介して送信される情報内にビットエラー及びブロックエラーが発生することは珍しくない。基地局２のアップリンク受信機は、ビットエラー及びブロックエラー率を低減するように構成されている。基地局受信機の１つの不可欠な部分は、アップリンク通信チャネル３を推定するように構成されたチャネル推定ユニット１８である。
【００２５】
図４を参照すると、受信機３０の受信部１４は、アップリンク通信チャネル３上で送信された信号ｙｎ（ｔ）を受信するように構成されている。サイクリックプレフィックスは、上記サイクリックプレフィックスを除去するように構成されたユニット１５内で受信信号ｙｎ（ｔ）から除去される。サイクリックプレフィックス除去ユニット１５は、時間領域内の受信信号を離散周波数領域信号に変換するように構成されたＤＦＴ変換器１６に接続されている。一実施例によれば、ＤＦＴ変換器７は、高速フーリエ変換（ＦＦＴ）として実施される。ＤＦＴ変換器１６は、抽出ユニット１７に接続されている。抽出ユニット１７は、ユーザごとに、上記ユーザに関連付けられた周波数範囲内にある信号の部分を抽出するように構成されている。抽出されたデータは、所与のユーザｕについてＲｕと呼ばれる。
【００２６】
図５を参照すると、チャネル推定ユニット１８は、既知のパイロットを含むサンプルＲｕを受信するように構成されている。この実施例では、受信したサンプルＲｕは、Ｎｃ×Ｎｒのサイズの行列であるＮｃ個のサブキャリアとＮｒ個の受信機アンテナとを備えた現在のユーザの１つの基準ＳＣ−ＦＤＭＡシンボルに対応する。チャネル推定ユニットは、既知のパイロットを含む受信サンプルＲｕに基づいてアップリンク通信チャネル３のチャネル推定値を決定するように構成されている。
【００２７】
はじめに、本明細書で提案するチャネル推定の理論を説明する。受信サンプルは、
Ｒｕ＝ＤｕＶＮｈ＋Ｗ
とモデル化できる。ここで、Ｗは、加法的ガウス雑音を含むＮｃ×Ｎｒのサイズの行列である。
ｈは、時間領域チャネルインパルス応答のモデルを含むＬ×Ｎｒのサイズの行列、Ｌは、チャネルのタップ数で、関連するチャネルタップ数を含む時間領域チャネルインパルス応答のモデルについては以下に詳細に説明する。
ＶＮは、基本期間Ｎを備えたＤＦＴを含むＮｃ×Ｌのサイズの行列、添え字Ｎは、ＤＦＴの基本期間を示し、Ｎの適した値については以下に詳細に説明する。
Ｄｕは、所定のユーザｕの主対角の既知のパイロットを含むＮｃ×Ｎｃのサイズの周波数領域内の対角行列である。
【００２８】
ＤＦＴ行列は、
【数３】
2010526511.tif 000005
に等しい。このＤＦＴ行列は、以下の説明から分かるように、チャネルモデルｈの推定値の解像度、位置及び長さを設計するために使用される。
【００２９】
まず、Ｎの値について説明する。ＤＦＴ行列ＶＮを使用する時の重なるフォールディングを回避するため、基本期間Ｎは、割り当てられたサブキャリアの数Ｎｃより大きいか、又はそれに等しくなければならない。すなわち、Ｎ≧Ｎｃである。
【００３０】
しかし、基本期間Ｎの上限は、最小二乗（ＬＳ）チャネル推定値を考慮することで導出することができる。チャネルｈの最小二乗（ＬＳ）推定値は、
【数４】
2010526511.tif 000006
に等しく、＊は、複素共役及び移項（エルミート演算子）を表す。パイロットシンボルは、周波数領域内で大きさの１単位を有するものと仮定する。すなわち、Ｄ＊ｕＤｕ＝Ｉで、Ｉは恒等行列である。次に、チャネルｈの推定値は、次のように記述することができる。
【数５】
2010526511.tif 000007
【００３１】
最小二乗チャネル推定値に含まれるＬ×Ｌ行列Ｖ＊ＮＶＮを、Λ＝Ｖ＊ＮＶＮで示す。ただし、Λの行ｒと列ｑ内の要素は
【数６】
2010526511.tif 000008
に等しい。Λのオフ対角要素はゼロである。すなわち、Λ（ｒ，ｑ）＝０（Ｎｃ（ｒ−ｑ）／Ｎ＝Ｋ（ここで、Ｋは整数）の場合、ｒ≠ｑ）である。従って、行列Λは、Ｎ＝Ｎｃ／Ｋならば、スカラー定数と置換することができる。これは、チャネルタップ数Ｌがいくつであってもあてはまる。
【００３２】
Ｎ≧Ｎｃ（重なるフォールディングを回避するため）及びＮ＝Ｎｃ／Ｋ（Λをスカラー定数と置換するため）の２つの要件を組み合わせると、Ｎ＝Ｎｃとなる。
【００３３】
Ｎ＝Ｎｃを選択すると、
【数７】
2010526511.tif 000009
となる。従って、チャネルｈの推定値は、
【数８】
2010526511.tif 000010
と記述することができる。
【００３４】
しかし、このチャネル推定値
【数９】
2010526511.tif 000011
を使用することは、以下に示す１つの欠点を有する。
【００３５】
【数１０】
2010526511.tif 000012
になるように、＾ＨＭＦを周波数領域内の一致フィルタチャネル推定値を含む行列Ｎｃ×Ｎｒとして表す。
ここで、行ｋ（周波数インデックスに対応する）及び列ｒ（受信アンテナに対応する）内の要素は、
【数１１】
2010526511.tif 000013
で表される。
【００３６】
チャネルタップ数を低減しない時間領域チャネル推定値は、周波数領域内の一致したフィルタにＮｃ×Ｎｃのサイズの逆ＤＦＴ行列Ｖ＊Ｎを適用することで達成される。時間領域内のこの一致フィルタチャネル推定値は、
【数１２】
2010526511.tif 000014
で表されるＮｃ×Ｎｒ行列であり、行ｍ（時間インデックスに対応する）及び列ｒ（受信アンテナに対応する）内の要素は、
【数１３】
2010526511.tif 000015
に等しい（ｍ＝０，Ｋ，Ｎｃ−１、及びｒ＝０，Ｋ，Ｎｒ−１の場合）。このチャネル推定値は、Ｎｃ個のチャネルタップを有することに留意されたい。すなわち、チャネルタップ数の低減は行われていない。時間領域内でＬ個のチャネルタップを選択する処理は、一致フィルタチャネル推定値に矩形ウィンドウを適用することで
【数１４】
2010526511.tif 000016
と表現することができる。
ここで、
【数１５】
2010526511.tif 000017
である。
【００３７】
上記のように、Ｌ個のタップの選択は、Ｎｃ×Ｎｃの代わりにＬ×Ｎｃのサイズの逆ＤＦＴ行列Ｖ＊Ｎを適用することで実行される。一実施例では、Ｌ個のタップの選択は、上記のように、一致フィルタ＾ｈＭＦにウィンドウを適用することで実行される。
【００３８】
図６にウィンドウ処理の周波数領域解釈を示す。この図で、Ｎｃ＝１２０、Ｎ＝Ｎｃ、及びＬ＝３、６及び１２である。図から分かるように、サブキャリアインデックスの数が大きいチャネルは、サブキャリアインデックスゼロのチャネル推定値に影響する。ここで、時間領域内でチャネルタップ数が低減した推定チャネルは、最大のサブキャリアインデックスの値が最小のサブキャリアインデックスの値にほぼ等しい周波数応答を有する。これは、明らかに推定チャネルに対する不要な制約である。１つの例外は、平坦な周波数応答を有するタップが１つのチャネルであり、その場合、円形畳み込みが当然制約を課さない。しかし、その他のすべてのチャネルでは、この円形畳み込みは望ましくない。
【００３９】
この欠点を克服するために、周波数領域への変換後に、時間領域チャネル推定値を補間してからチャネルタップ数を低減する。上記補間の導入によって円形畳み込みの影響が低減する。時間領域補間を提供する１つの方法は、周波数領域へインパルス応答を変換し、ゼロパディング（zero padding）を追加し、次に、時間領域への変換を行ってからチャネルタップ数を低減する方法である。詳細は、行列表記で表されるこの補間チャネル推定値は、
【数１６】
2010526511.tif 000018
で表される。ここで、
【数１７】
2010526511.tif 000019
は、基本期間Ｎｃを含むＤＦＴ行列で、
【数１８】
2010526511.tif 000020
である。従って、Ｖ＊ＮｃＫｏｓは、基本期間ＮｃＫｏｓを含む逆ＤＦＴ行列（ただし、Ｋｏｓは、オーバーサンプリング係数）である。前述の最小二乗チャネル推定値を挿入すると、
【数１９】
2010526511.tif 000021
が得られる。
【００４０】
図７にこの補間チャネル推定値＾ｈｏｓの周波数領域解釈を示す。この図では、Ｎ＝Ｋｏｓ・Ｎｃ、Ｎｃ＝１２０、Ｋｏｓ＝２、Ｌ＝３、６及び１２である。図から分かるように、大きいインデックスを備えたサブキャリアの小さいインデックスのサブキャリアでの補間チャネル推定値への影響はきわめて小さい。オーバーサンプリング係数Ｋｏｓは、一実施例では、１．５〜２．５の範囲内である。
【００４１】
上記から分かるように、オーバーサンプリング係数Ｋｏｓの導入は、チャネル推定値に近似が導入されるということを意味する。ここで、
【数２０】
2010526511.tif 000022
と仮定するが、これは、Ｎ＝Ｎｃの場合にのみ真である。しかし、この近似で比較的小さいエラーが導入される。
【００４２】
従って、以下のチャネル推定値が提案される。
【数２１】
2010526511.tif 000023
ここで、検索ウィンドウの事前選択された長さＬからのチャネルタップ数の低減が＾ｈｏｓからいくつかの行を選択するだけで実行することができる。チャネルタップ数は、使用したチャネル推定値の最大遅延拡散に関連して選択することができる。対角行列Ｄ＊ｕを用いた乗算は、Ｎｃ回の複素乗算だけで実行され、Ｖ＊ＫｏｓＮｃを用いた乗算は、Ｎｃ個の入力値がゼロでない基本期間ＫｏｓＮｃを含む逆高速フーリエ変換（ＩＦＦＴ）を用いて実施することができることに留意されたい。
【００４３】
それ故、図５のチャネル推定ユニット１８は、ＤＴＦ行列計算器１９を備える。ＤＴＦ行列計算器１９は、基本期間ＮｃＫｏｓを備えるものとして以下のようにＤＦＴ行列を計算するように構成されている。
【数２２】
2010526511.tif 000024
チャネル推定ユニット１８は、パイロット行列Ｄｕを形成するように構成されたパイロット行列計算器２０をさらに備える。パイロット行列Ｄｕは、特定のユーザｕの主対角上に既知のパイロットを含むサイズがＮｃ×Ｎｃの周波数領域内の対角行列である。チャネル推定部２１は、チャネル推定値＾ｈｏｓを提供するように構成されている。従って、チャネル推定部２１は、ＤＴＦ行列計算器１９によって提供されるＤＦＴ行列に対して複素共役及び移項（エルミート演算子）Ｖ＊ＫｏｓＮｃを形成するように構成されている。チャネル推定部２１は、パイロット行列計算器２０によって提供されるパイロット行列Ｄｕに対して複素共役及び移項（エルミート演算子）を形成するようにさらに構成されている。チャネル推定部２１は、その後、チャネル推定値を
【数２３】
2010526511.tif 000025
として計算するように構成されている。一実施例では、チャネル推定部２１は、チャネルタップ低減部２２にチャネル推定値＾ｈｏｓを供給するように構成されている。チャネルタップ推定器２２は、本明細書に記載するチャネル推定に限らず、いかなる方法で推定されたチャネルにも適用可能である。
【００４４】
チャネルタップ低減部２２によって使用される、チャネル推定値＾ｈｏｓを用いて推定されるチャネル応答の最大遅延拡散に基づいてチャネルタップ数を低減する手順について以下に説明する。
【００４５】
図８で、チャネル推定値＾ｈｏｓの位置及び長さがチャネルタップ位置を適合させる２つのステージで決定される。第１に、開始位置Ｎｓｔａｒｔを一定に保ちながらインパルス応答の開始位置Ｎｓｔａｒｔから右側にチャネル長を調整することで、チャネルの終点＾ｍＲが推定される。次に、第２のステージで、終点＾ｍＲを一定に保ちながらインパルス応答の左側のチャネルの始点＾ｍＬが調整される。
【００４６】
別の実施形態（図示せず）では、位置Ｎｓｔａｒｔを一定に保ちながらインパルス応答の開始位置Ｎｓｔａｒｔから左側にチャネル長を調整することで、始点＾ｍＬが最初に決定される。次に、第２のステージで、始点＾ｍＬを一定に保ちながらインパルス応答の右側のチャネルの終点＾ｍＲが調整される。
【００４７】
始点＾ｍＬ及び終点＾ｍＲは、一実施例では、インパルス応答内のエネルギーの所定量が始点＾ｍＬと終点＾ｍＲとの間の範囲に収まるように選択される。別の実施例では、始点＾ｍＬ及び終点＾ｍＲは、始点＾ｍＬと終点＾ｍＲとの間の範囲外のインパルス応答＾ｈｏｓの振幅が所定の値を下回るように選択される。
【００４８】
図９は、始点＾ｍＬ及び終点＾ｍＲを決定する方法を実施する一実施例がチャネルタップ低減部２２で実施される様子を示す。一致フィルタ推定器２３は、時間領域内の一致フィルタチャネル推定値を計算するように構成されている。一致フィルタチャネル推定値は、
【数２４】
2010526511.tif 000026
として計算される。この一致フィルタチャネル推定値は、上記チャネル推定値＾ｈｏｓに対応する。
【００４９】
共分散行列計算器２４は、受信信号＾ｈｏｓの共分散行列ＰＹを計算するように構成されている。受信信号の共分散行列ＰＹは、ＰＹ（ｓ）＝Ｒ＊ｕＲｕとして計算される。共分散行列ＰＹは、エルミート行列であり、従って、オフ対角要素の半分だけを計算すればよい。
【００５０】
右方への適応チャネル推定器２５は、一致フィルタ推定器２３からの一致フィルタチャネル推定値＾ｈＭＦと共分散行列計算器２４からの共分散行列ＰＹとに基づいて終点＾ｍＲを決定するように構成されている。右方への適応チャネル推定器２５は、開始インデックスＮｓｔａｒｔを一定に保ちながら終了インデックスｍＲを変更するように構成されている。開始インデックスＮｓｔａｒｔは、設計パラメータであり、以下に詳細に説明する。評価される終了インデックスｍＲの値はすべて、ＮｓｔａｒｔとＬ−１との間のＮΔ刻みの整数で、ＮΔは、終点＾ｍＲを決定する際に用いる解像度を表す。
【００５１】
右方への適応チャネル推定器２５は、
【数２５】
2010526511.tif 000027
として、又は行列表記
【数２６】
2010526511.tif 000028
ここで、
【数２７】
2010526511.tif 000029
を用いて残留雑音共分散行列ＱＲを計算するように構成されている。
【００５２】
一実施例では、残留雑音共分散行列ＱＲは反復して計算される。この実施例によれば、
【数２８】
2010526511.tif 000030
であり、
＾ｍＲ＞Ｎｓｔａｒｔの場合、
【数２９】
2010526511.tif 000031
である。
【００５３】
従って、残留雑音行列ＱＲを計算するのに１回の複素乗算しか必要でない。さらに、残留雑音共分散行列ＱＲは、エルミート行列であり、オフ対角要素の半分だけを計算すればよい。
【００５４】
右方への適応チャネル推定器２５は、その後、残留雑音共分散行列ＱＲに基づいてＮｓｔａｒｔとＬ−１との間のＮΔ刻みの可能な終了インデックスｍＲの各々の決定値βｒｉｇｈｔを計算するように構成されている。決定値βｒｉｇｈｔは、一実施例では、
【数３０】
2010526511.tif 000032
として計算され、｜ＱＲ（ｍＲ）｜は、残留雑音共分散行列ＱＲ（ｍＲ）の行列式を表し、
【数３１】
2010526511.tif 000033
は、設計変数ＫＡＩＣを含む赤池様チャネル長処罰係数（Akaike like channel length punishment factor）である。一実施例では、設計変数は、ＫＡＩＣ＝２に設定され、それによって古典的赤池情報量基準が適用される。しかし、この設計変数ＫＡＩＣを増やすことで、より保守的なチャネル長選択を達成することができる。
【００５５】
２つの受信機アンテナを含む一実施例では、共分散行列の行列式は、
【数３２】
2010526511.tif 000034
に等しい。
【００５６】
行列式は、負ではない実数値である。
【００５７】
各終了インデックスｍＲについて決定値βｒｉｇｈｔが計算されると、決定値βｒｉｇｈｔの最小値に対応する終了インデックスｍＲの値が終点
【数３３】
2010526511.tif 000035
として決定される。
【００５８】
左方への適応チャネル推定器２６は、一致フィルタ推定器２３からの一致フィルタチャネル推定値＾ｈＭＦと、共分散行列計算器２４からの共分散行列ＰＹと、右方への適応チャネル推定器２５から入力される終点＾ｍＲとに基づいて始点＾ｍＬを決定するように構成されている。左方への適応チャネル推定器２６は、終点＾ｍＲを一定に保ちながら開始インデックスｍＬを変更するように構成されている。それ故、評価される開始インデックスの値はすべて、０と終点＾ｍＲとの間のＮΔ刻みの整数で、ＮΔは、始点＾ｍＬを決定する際に用いる解像度を表す。
【００５９】
左方への適応チャネル推定器２６は、
【数３４】
2010526511.tif 000036
として、又は行列表記ＱＬ（ｍＬ）＝ＫＳｃａｌｅＰＹ−ｈ＊０ｈ０、
ここで、
【数３５】
2010526511.tif 000037
を用いて残留雑音共分散行列ＱＬを計算するように構成されている。
【００６０】
一実施例では、残留雑音行列ＱＬは、反復して計算される。従って、残留雑音行列ＱＬは、
【数３６】
2010526511.tif 000038
と計算され、
ｍＬ＜＾ｍＲの場合、
【数３７】
2010526511.tif 000039
である。
【００６１】
従って、残留雑音行列ＱＬを計算するのに１回の複素乗算しか必要でない。さらに、残留雑音共分散行列ＱＬは、エルミート行列であり、オフ対角要素の半分だけを計算すればよい。
【００６２】
左方への適応チャネル推定器２６は、その後、残留雑音共分散行列ＱＬに基づいて０と終点＾ｍＲとの間のＮΔ刻みの可能な開始インデックスｍＬの各々の決定値βｌｅｆｔを計算するように構成されている。決定値βｌｅｆｔは、一実施例では、
【数３８】
2010526511.tif 000040
として計算される。
【００６３】
各開始インデックスｍＬに対して決定値βｌｅｆｔが計算されると、決定値βｌｅｆｔの最小値に対応する開始インデックスｍＬの値が始点
【数３９】
2010526511.tif 000041
として決定される。
【００６４】
チャネルタップ低減器２７は、チャネル推定ユニット１８のチャネル推定部２１によって提供されるチャネル推定値＾ｈｏｓと、左方への適応チャネル推定器２６によって推定される始点＾ｍＬと、右方への適応チャネル推定器２５によって推定される終点＾ｍＲとに基づいてチャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄを提供するように構成されている。チャネルタップ低減器２７は、始点＾ｍＬ及び終点＾ｍＲによって与えられる範囲外のチャネル推定値内の行を削除することでチャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄを提供するように構成されている。
【００６５】
一実施例では、受信機は、ユーザデータシンボルの周波数領域等化器を備える。これは、ユーザデータシンボルが例えばＳＣ−ＦＤＭＡシンボルの場合にあてはまる。この実施例によれば、ＤＴＦ変換器２８は、チャネルタップ低減器２７から入力されるチャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄに基づいて周波数領域内のチャネル推定値＾Ｈを計算するように構成されている。それ故、ＤＦＴ変換器２８は、周波数インデックスｋ＝０，Ｋ，Ｎｃ−１及び受信アンテナｒ＝０，Ｋ，Ｎｒ−１に対して時間領域チャネル推定値を周波数領域
【数４０】
2010526511.tif 000042
に変換するように構成されている。
【００６６】
タイミングミスアラインメント計算器２９は、予測されたチャネルタップ低減チャネル推定値（＾ｈｏｓ＿ｒｅｄ，ｅｘｐ）と現時点で計算されたチャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄとの間のタイミングミスアラインメントを決定するように構成されている。
【００６７】
タイミングミスアラインメント計算を詳細に説明する前に、右方への適応チャネル推定器２５によって使用される開始インデックスＮｓｔａｒｔについて説明する。一実施例では、開始インデックスＮｓｔａｒｔは、チャネル応答の受信の予測されたタイミングを表す。従って、開始インデックスＮｓｔａｒｔは、推定始点＾ｍＬ，ｅｘｐに関連して選択される。例えば、開始インデックスＮｓｔａｒｔは、Ｎｓｔａｒｔ＝＾ｍＬ，ｅｘｐとして選択することができる。別の実施例では、開始インデックスＮｓｔａｒｔは、Ｎｓｔａｒｔ＝＾ｍＬ，ｅｘｐ＋ｋｏｎｓｔ（ｋｏｎｓｔは、チャネルタップ数を低減したチャネル推定値の長さに基づく事前設定した定数又は定数のセットを表す）として選択される。
【００６８】
上記開始インデックスの選択に関する一実施例では、タイミングミスアラインメント計算器２９は、ほぼ予測始点＾ｍＬ，ｅｘｐと決定された始点＾ｍＬとの間の時間差に基づいてタイミングミスアラインメントを決定するように構成されている。タイミングミスアラインメント情報は、上記情報をチャネル推定値が計算された対象のユーザへ送信するように構成された送信ユニット３１に供給される。上記ユーザに関連付けられた端末１の受信ユニット１３は、タイミングミスアラインメント情報を受信するように構成されている。端末１は、上記タイミングミスアラインメント情報に基づいてその送信タイミングを調整するようにさらに構成されている。
【００６９】
図１０に、以下のステップを含む無線チャネルのチャネル推定値＾ｈｏｓを決定する方法を示す。第１のステップ３２で、各々がサイクリックプレフィックスを有する複数のシンボルを形成するいくつかの時間領域信号（ｙｎ）が受信される。第２のステップ３３で、プレフィックスを除去し受信サンプルのベクトル（Ｒｕ）を形成するように受信信号（ｙｎ）が処理される。第３のステップ３４で、基本期間Ｎｃ・ｋを有するＤＦＴ行列（Ｖ）の共役Ｖ＊Ｎｃが決定される。第４のステップ３５で、上記ユーザ（ｕ）に対応する周波数領域内のパイロットシンボルを含むパイロット行列Ｄｕの共役Ｄ＊ｕが決定される。第３のステップ３４と第４のステップ３５は、任意の順序で実行することができる。第５のステップで、チャネル推定値が
【数４１】
2010526511.tif 000043
として決定される。ここで、計算回数が少なく事前計算される値が小さいチャネル推定値を提供する必要がある。
【００７０】
しかし、第６のステップ３７、第７のステップ３８及び第８のステップ３９に記載するように、チャネル推定値のサイズは、チャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄを形成するように、チャネル推定値
【数４２】
2010526511.tif 000044
の行列から行のサブセットを選択することでさらに低減することができる。第６のステップ３７で、チャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄの終点＾ｍＲが推定される。この推定は、一実施例では、開始サンプル値Ｎｓｔａｒｔから始まる右方向のチャネル推定値＾ｈｏｓのいくつかの可能な長さの各々に対して決定値を決定するステップと、最も好適な決定値に関連付けられた長さを選択するステップとを含む。第７のステップ３７では、チャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄの始点（＾ｍＬ）が推定される。一実施例では、この推定は、左方向のチャネル推定値＾ｈｏｓのいくつかの可能な長さの各々に対して決定値を決定するステップと、最も好適な決定値に関連付けられた長さを選択するステップとを含む。これらの長さは、一実施例では、推定終点＾ｍＲに関連する。第８のステップ３９では、チャネル推定値が始点と終点との間の行のみを含むように低減される。
【００７１】
一実施例では、第９のステップ４０で、チャネルタップ低減チャネル推定値の予測されたタイミングと決定されたチャネルタップ低減チャネル推定値＾ｈｏｓ＿ｒｅｄとの間のタイミングミスアラインメントが決定される。一実施例では、第１０のステップ４１で、タイミングミスアラインメントに関連する情報がユーザ（ｕ）へ送信される。
【図１】
2010526511.tif 000049
【図２】
2010526511.tif 000050
【図３】
2010526511.tif 000051
【図４】
2010526511.tif 000052
【図５】
2010526511.tif 000053
【図６】
2010526511.tif 000054
【図７】
2010526511.tif 000055
【図９】
2010526511.tif 000056
【図１０】
2010526511.tif 000057
【図８】
2010526511.tif 000058
【国際調査報告】
2010526511.tif 000059
2010526511.tif 000060
2010526511.tif 000061
2010526511.tif 000062