(57)【特許請求の範囲】
【請求項１】
エミュレータ（４１８）の疑似無線チャネルを介して、試験対象である電気機器（４００）と通信する方法であって、
複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）の重みを決定するために、アンテナ素子で取得された論理空間相互相関及び空間相関の費用関数を最適化するステップと、
これら重みに基づいて、少なくとも２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で、疑似無線チャネルの少なくとも１つのビーム（５００、７００）を形成するステップと、
無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器（４００）と通信するために、エミュレータ（４１８）と接続された複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に対して疑似無線チャネルをシフトするステップ（３０００）と
を備えていることを特徴とする方法。
【請求項２】
請求項１記載の方法において、該方法は、少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の少なくとも１つのビーム（５００、７００）を、様々な方向で試験対象機器（４００）に向けることにより、疑似無線チャネルをシフトするステップ（２７００）を備えていることを特徴とする方法。
【請求項３】
請求項１記載の方法において、該方法は、２つの連続する通信時点の間に、電気機器（４００）に対する２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の角度（Δθ）よりも小さい分又は大きい分だけビーム（５００、７００）の方向を変更するために、２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の間の距離（ｄｓ）より短い分又は長い分だけ、疑似無線チャネルの電力角度スペクトルをシフトするステップを備えていることを特徴とする方法。
【請求項４】
請求項１の方法において、該方法は、２つの連続する通信時点の間に、電気機器（４００）に対する２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の角度（Δθ）分だけビーム（５００、７００）の方向を変更するために、２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の間の距離（ｄｓ）分だけ、疑似無線チャネルの電力角度スペクトルをシフトするステップを特徴とする方法。
【請求項５】
請求項１記載の方法において、該方法は、
第１の時点において、少なくとも１つのパスの信号の少なくとも１つのビーム（５００、７００）を、少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で形成する（２３０２）ために、疑似無線チャネルを複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に与えるステップと、
疑似無線チャネルを、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に対して新しい位置に設定するステップと、
第２の時点において、少なくとも１つのビーム（５００、７００）を、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）の少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で形成するために、新しい位置に設定された無線チャネルを、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に与えるステップと
を備えていることを特徴とする方法。
【請求項６】
エミュレータ（４１８）の疑似無線チャネルを介して、試験対象である電気機器（４００）と通信する試験システムであって、該試験システムは、
複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）の重みを決定するために、アンテナ素子で取得された論理空間相互相関及び空間相関の費用関数を最適化し、
これら重みに基づいて、少なくとも２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で、疑似無線チャネルの少なくとも１つのビーム（５００、７００）を形成し、
無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器（４００）と通信するために、エミュレータ（４１８）と接続された複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に対して疑似無線チャネルをシフトする
よう構成されていることを特徴とする試験システム。
【請求項７】
請求項６記載の試験システムにおいて、該試験システムは、少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の少なくとも１つのビーム（５００、７００）を、複数の異なる方向で試験対象機器（４００）に向けることにより、疑似無線チャネルをシフトするよう構成されていることを特徴とする試験システム。
【請求項８】
請求項６記載の試験システムにおいて、該試験システムは、２つの連続する通信時点の間で、電気機器（４００）に対する２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の角度（Δθ）よりも小さい分又は大きい分だけビーム（５００、７００）の方向を変更するために、２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）間の距離（ｄｓ）より短い分又は長い分だけ、疑似無線チャネルの電力角度スペクトルをシフトするよう構成されていることを特徴とする試験システム。
【請求項９】
請求項６記載の試験システムにおいて、該試験システムは、２つの連続する通信時点の間で、電気機器（４００）に対する２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）の角度（Δθ）分だけビーム（５００、７００）の方向を変更するために、２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）間の距離（ｄｓ）分だけ、疑似無線チャネルの電力角度スペクトルをシフトするよう構成されていることを特徴とする試験システム。
【請求項１０】
請求項６記載の試験システムにおいて、試験システムは、
第１の時点において、少なくとも１つのパスの信号の少なくとも１つのビーム（５００、７００）を、少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で形成する（２３０２）ために、疑似無線チャネルを、エミュレータ（４１８）と接続された複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に与え、
疑似無線チャネルを、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に対して新しい位置に設定し、
第２の時点において、少なくとも１つのビーム（５００、７００）を、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）の少なくとも１つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で形成するために、新しい位置に設定された無線チャネルを、複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に与える
よう構成されていることを特徴とする試験システム。
【請求項１１】
疑似無線チャネルを介して、試験対象である電気機器（４００）と通信するエミュレータであって、該エミュレータは、
複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）の重みを決定するために、アンテナ素子で取得された論理空間相互相関及び空間相関の費用関数を最適化し、
これら重みに基づいて、少なくとも２つのアンテナ素子（４０２〜４１６、６００２〜６０１０）で、疑似無線チャネルの少なくとも１つのビーム（５００、７００）を形成し、
無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器（４００）と通信するために、エミュレータ（４１８）と接続された複数のアンテナ素子（４０２〜４１６、６００２〜６０１０）に対して疑似無線チャネルをシフトする
よう構成されていることを特徴とするエミュレータ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、無響チャンバ内での機器のオーバーザエアー（ＯＴＡ）試験システムに関する。
【背景技術】
【０００２】
無線周波数信号が送信機から受信機に送信されると、信号は、さまざまな到来角、信号遅延、偏波、及び電力を有する１又は複数のパスに沿って無線チャネル内で伝搬し、これにより、受信信号において様々な持続時間及び強度のフェージングが生じる。加えて、他の送信機によるノイズ及び干渉が、無線接続に干渉する。
【０００３】
送信機及び受信機について、現実の状況をエミュレートする無線チャネルエミュレータを用いて試験することが可能である。デジタル無線チャネルエミュレータでは、チャネルはＦＩＲフィルタでモデル化されるのが通常であり、これはさまざまな遅延によって遅延された信号を、チャネル係数すなわちタップ係数によって重み付けすることと、重み付けされた信号成分を合算することとにより、チャネルモデル及び入力信号間の畳み込みを形成する。チャネル係数は、実際のチャネルの時間的な振る舞いに対応する時間の関数である。無線チャネルのエミュレーション試験は、導電性の接続を介して実行されて、送信機及び受信機はケーブルで接続されるのが通常である。
【０００４】
加入者端末及び無線システムの基地局間の通信を、実際の加入者端末が無響チャンバ（anechoic chamber）内のエミュレータの複数のアンテナ素子に囲まれているようなＯＴＡ（Over The Air：オーバーザエアー）試験を用いて、テストすることができる。エミュレータは、チャネルモデルに従って、加入者端末及び基地局間の通信経路（パス）をエミュレートする基地局と接続されているか、又はそのような基地局として動作する。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
通信が幾つかの方向から生じる場合の加入者端末をテストするために、加入者端末は機械的に回転される。しかしながら、加入者端末及びエミュレータ間のケーブルは充分に回転することができないために、回転によってケーブリングに問題が生じる。さらに、モータ及びモータで駆動するターニングテーブル等の回転メカニズムの利用及びその制御により、試験システムの複雑さが増大する。したがって、より優れたＯＴＡ試験ソリューションが必要とされている。
本発明の目的は、ソリューションの向上を提供するものである。
【課題を解決するための手段】
【０００６】
本発明の一態様によると、エミュレータの疑似無線チャネル（simulated radio channel）を介して、試験対象電気機器と通信する方法が提供される。この方法は、無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器と通信するために、エミュレータと接続された複数のアンテナ素子に対して疑似無線チャネルをシフトすることを特徴とする。
本発明の別の態様によると、エミュレータの疑似無線チャネルを介して、試験対象電気機器と通信する試験システムが提供される。試験システムは、無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器と通信するために、エミュレータと接続された複数のアンテナ素子に対して疑似無線チャネルをシフトするよう構成されている。
本発明の別の態様によると、疑似無線チャネルを介して試験対象電気機器と通信するエ
ミュレータが提供される。エミュレータは、無響チャンバにおいて疑似無線チャネルに対して様々な方向を用いて試験対象機器と通信するために、エミュレータと接続された複数のアンテナ素子に対して疑似無線チャネルをシフトするよう構成されている。
【０００７】
本発明は、幾つかの作用効果を奏することができる。疑似無線チャネルをＤＵＴ（試験対象機器）に対して電気的にシフトすることにより、ケーブル又は複雑性の問題が生じることなく、ＤＵＴを異なる方向から試験することができる。
【図面の簡単な説明】
【０００８】
【図１】無線信号の伝播を示す図である。
【図２】受信ビームの電力アジマススペクトルを示す図である。
【図３】送信ビームの電力アジマススペクトルを示す図である。
【図４】ＯＴＡ試験チャンバ内の測定構成を示す図である。
【図５】アンテナ素子によってモデル化されるビームを示す図である。
【図６】アンテナ素子群及び関連するアンテナ群スイッチングネットワークを示す図である。
【図７】アンテナ素子群に囲まれているＤＵＴを示す図である。
【図８】ＭＩＭＯ構成におけるアンテナの遅延の制御を表す図である。
【図９】ＯＴＡチャンバにおけるアンテナの遅延の制御を表す図である。
【図１０】ＯＴＡチャンバにおけるＡｏＡを表す図である。
【図１１】図１０のアンテナ素子のアンテナ重み付けを表す図である。
【図１２】３つのアンテナ素子との空間相関を表す図である。
【図１３】アンテナ素子の重み付け及び結果として生じたＰＡＳを表す図である。
【図１４】論理的相関及び理想的な空間相関を表す図である。
【図１５】６つのクラスタのＰＤＰを表す図である。
【図１６】８つのチャネルの遅延タップマッピングを表す図である。
【図１７】ＤＵＴの空間分解能がアンテナ素子の間隔よりも狭い状況を表す図である。
【図１８】図１７の状況におけるＰＡＳを表す図である。
【図１９】ＤＵＴの空間分解能がアンテナ素子の間隔よりも広い状況を表す図である。
【図２０】図１９の状況におけるＰＡＳを表す図である。
【図２１】３つ及び５つのアンテナ素子を有するＰＡＳを表す図である。
【図２２】アンテナ素子の偏波を表す図である。
【図２３】一状態にある、スイッチを有するエミュレータを示す図である。
【図２４】別の状態にある、スイッチ有するエミュレータを示す図である。
【図２５】ＦＩＲフィルタを有するエミュレータを示す図である。
【図２６】ＦＩＲフィルタを示す図である。
【図２７】一位置にある２つのビームを有するＯＴＡチャンバを示す図である。
【図２８】別の位置にシフトされた２つのビームを有するＯＴＡチャンバを示す図である。
【図２９】減衰器及び移相器を有するエミュレータを示す図である。
【図３０】方法のフロー図を示す図である。
【発明を実施するための形態】
【０００９】
ＤＵＴに対する正確な相関、到来角及び偏波特性を可能にするために、ＯＴＡにおいてチャネルインパルス応答及びアンテナ重み付けの最適化が生成される。上記のソリューションで使用するのは、例えば送信機、多次元無線チャネルエミュレータ(multidimensional radio channel emulator)、無響チャンバ、無響チャンバ内の別々の無線チャネルと接続されたアンテナ素子、及び無響チャンバ内のＤＵＴである。疑似無線チャネルは、ＤＵＴ上の到来角に対してシフトされ、様々な時点の様々な角度の通信に、同一の電力角度スペクトルを使用できるようにする。
【００１０】
図１は、送信機及び受信機間の無線信号の伝搬を示している。送信機１００は、少なくとも１つのアンテナ素子１０４〜１１０からなるアンテナ１０２を備えている。アンテナは、例えば、アンテナ素子間の間隔が、例えば無線信号の波長の半分等のように一定である、ＵＬＡ（均一線形アレイ）アンテナである。この例では、送信機１００は無線システムの基地局である。同様にして、受信機１１２は、少なくとも１つのアンテナ素子１１６〜１２２からなるアンテナ１１４を備えている。この例では、受信機１１２は無線システムの加入者端末である。送信機１００が無線信号を送信すると、送信ビーム１２４は角φ１に向けられ、そしてその角度広がりであるδφは、ｘδφstdである。ここでｘはゼロよりも大きな実数であり、δφstdは角度φ１の基準偏差である。送信ビーム１２４は、放射線を反射及び／又は散乱させる少なくとも１つのクラスタ１２６、１２８に当たる。クラスタ１２６及び１２８はそれぞれ、クラスタ１２６、１２８内で主に反射及び／散乱する幾つかのアクティブな領域１２６０〜１２６４、１２８０〜１２８４を有している。クラスタ１２６、１２８は固定されているか又は運動しており、かつクラスタ１２６、１２８は、建造物、列車、山等、自然の物体又は人工の物体である。アクティブな領域は、オブジェクト上のいくつかのより微細な構造的特徴である。
【００１１】
反射及び／又は散乱されたビームは、受信機１１２のアンテナ１１４に対して向けられる。アンテナ１１４は受信角ψ1を有しており、その角度広がりであるδψは、ｙδψstdである。ここでｙはゼロより大きな実数であり、δψstdは角ψ１の基準偏差である。そして、クラスタ１２６から反射及び／又は散乱されたビーム１３０が受信される。同様に、アンテナ１１４は受信角ψ２からのビームも有しており、その角度広がりはδψ２である。少なくとも１つのクラスタ１２６、１２８を介した送信機１００から受信機１１２への伝搬は、見通し線に沿ってまっすぐ進む信号に対して、信号に追加的な遅延を生じさせる。
【００１２】
無線チャネル内のクラスタ１２６、１２８がマルチパス伝搬を生じさせている。概して、パスとクラスタ１２６、１２８は、１つの受信パスが１つのクラスタから到来するという対応を有している。したがって、無線チャネルは、クラスタの電力、遅延、公称のＡｏＡ（到来角）及びＡｏＤ（発射角）、並びに到来及び発射の両端におけるクラスタの角度広がりによって説明されてもよい。さらに加えて、送信機のアンテナアレイの情報が必要である。この情報は、アンテナアレイの形状及びアンテナのフィールドパターン（ビーム）のパラメータ値を含む。またさらに、加入者端末の速度ベクトル及び／又はクラスタのドップラ周波数成分も必要とされる場合もある。
【００１３】
表１は、市街地環境における、無線チャネルのクラスタ遅延線モデル(clustered delay
line model)の例を表している。クラスタ１及び３は、異なる遅延及び電力を有する３つのアクティブな領域を有している。
【表１】
0005358021.tif 000002
【００１４】
ＡＳＤ（発射角広がり）は全てのクラスタで一定であると仮定し、この例ではＡＳＤ＝２°とする。同様に、ＡＳＡ（到来角広がり）も全てのクラスタで一定であると仮定し、この例ではＡＳＡ＝１５°である。加えて、ＸＰＲ（Cross Polarization Power ratio:交差偏波電力比）も全てのクラスタで一定であると仮定し、この例ではＸＰＲ＝７ｄＢである。これらはクラスタによって異なっていてもよい。
【００１５】
無線チャネルのインパルス応答推定Ｈｕ，ｓ，ｎ（ｔ，τ）は、次の数学的形式で表される。
【数１】
0005358021.tif 000003
ここでＦｔｘ，ｓは送信アンテナのフィールドパターン、Ｆｒｘ，ｕは受信アンテナのフィールドパターン、ｄｓは、ＵＬＡ送信アンテナ内の連続する２つのアンテナ素子間の距離、ｄｕは、ＵＬＡ受信アンテナ内の連続する２つのアンテナ素子間の距離、ｋは波数（ｋ＝２π／λ０で、λ０は無線信号の波長）であり、Ｐｎはクラスタ電力、Ｍはクラスタ内のアクティブな領域の数を意味し、ｍはアクティブな領域のインデクス、ｎはクラスタのインデクス、Φｎ，ｍは散乱ｎ，ｍの一定位相項（a constant phase term）、νｎ
，ｍはインデクスｎ，ｍを有するアクティブな領域のドップラ周波数であり、τは遅延である。
【００１６】
インデクスｎ，ｍを有するアクティブな領域のドップラ周波数は、次のように表現される。
【数２】
0005358021.tif 000004
ここでｖ（ただし、ｖの上部に−）は速度ベクトルであり、‖ｖ‖（ｖの上部に―）はアクティブな領域及び受信機間の相対速度である。
【００１７】
式（１）におけるインパルス応答推定は、受信アンテナを全方向性であると仮定すると、次のような形式に単純化される。
【数３】
0005358021.tif 000005
【００１８】
インパルス応答推定は無線チャネル推定とも称され、それに従って無線チャネルがその中を伝搬する信号を歪ませるモデルである。ＳＣＭ（空間チャネルモデル）、ＳＣＭＥ（SCM Extended）、ＷＩＮＮＥＲ（Wireless World Initiative New Radio）、及びＩＭＴ−Ａｄｖａｎｃｅｄ（International Mobile Telecommunications）のようなチャネルモデルは、双方向性クラスタを含んだ形状モデルである。送信機端及び受信機端の電力アジマススペクトルは、図２及び図３の例のようになる。チャネル係数の生成は、式（１）において放射線（別々の方向）を合算することによって実行されるが、クラスタは上述のパラメータによって定義される。言い換えると、鏡面散乱は本モデルに不可欠な部分ではなく、チャネル係数を生成するための道具に過ぎない。
ドップラ及び可能な基地局のアンテナ相関、並びにチャネル電力遅延プロファイル等のフェージングはチャネル係数に含まれている。
【００１９】
ＤＵＴのアンテナ相関及び他のＤＵＴのアンテナ効果（影響）のみが、ＯＴＡチャンバ内の実際の無線送信に関連していない。
【数４】
0005358021.tif 000006
ドップラ周波数νｎ，ｍは、ＡｏＡ角に基づいて決定される。計算結果は、ベクトル係数Ｈｎ（ｔ，τ）を有する離散インパルス応答である。Ｈｎ（ｔ，τ）の次元は、１×Ｓであり、ここでＳは基地局アンテナの数である。
【００２０】
クラスタｎのマッピングは、クラスタの公称の方向及びクラスタの角度広がりに依存して、適宜のエミュレータチャネル及びＯＴＡアンテナに対して実行される。
２つのＯＴＡアンテナによるクラスタの到来角度広がりの概算は、不正確さの元になる。このことは、まばらなＯＴＡアンテナ配置及び狭いクラスタ(narrow clusters)である場合、すなわちΔθ＞＞ＡｏＡ角度広がりの場合に、特に不正確となる。ＤＵＴ側でのクラスタの角度広がり値は、例えば、ＳＣＭでは３５°、ＷＩＮＮＥＲでは３°〜２２°、ＩＭＴ−Ａｄｖａｎｃｅｄでは３°〜２２°、及びＴＧｎでは１４°〜５５°であり、モデルのシナリオに依存する。
可能なＤＵＴアンテナ間の非相関を生成するために、単一のクラスタを少なくとも２つ
のＯＴＡアンテナに分割することが必要な場合もある。信号が単一のＯＴＡアンテナから送信された場合は、角度広がりがなくＤＵＴと完全に相関している鏡面反射の場合に等しい。
【００２１】
図２は、５つのクラスタからの受信ビームの電力アジマススペクトルを示している。図２では、Ｘ軸は度で表示された角度であり、Ｙ軸がデシベルで表示された電力である。５つのビーム２００、２０２、２０４、２０６、及び２０８は、様々な到来角で受信される。ビーム２００、２０２、２０４、２０６、及び２０８は、様々な時点で受信される。すなわち、それらの少なくとも１つは、残りのビームとは異なる遅延を有している。
【００２２】
図３は、図２の例に係るものと同一の５つのクラスタに対する送信ビームの、電力アジマススペクトルを示している。図３では、Ｘ軸は度で表示された角度であり、Ｙ軸がデシベルで表示された電力である。５つのビーム３００、３０２、３０４、３０６、及び３０８は、わずかに異なる発射角で送信されるが、それは反射及び／又は散乱しているクラスタが、角度において極僅かしか分散していないためである。
【００２３】
図４は、ＯＴＡ試験チャンバを表している。ＤＵＴ４００が中心にあり、チャンバアンテナ素子４０２、４０４、４０６、４０８、４１０、４１２、４１４、及び４１６は、ＤＵＴ４００の周りに円形状に均一間隔（例えば、８つの素子であれば４５°）で存在している。Ｋ個のＯＴＡアンテナの方向を、θｋ、ここでｋ＝1,…,Ｋ、かつ角度ドメインにおけるアンテナの間隔をΔθとする。角度Δθは、電気機器４００に対する２つのアンテナ素子４０２〜４１６の間隔の測定値を表している。アンテナ素子のそれぞれが単一のエミュレータ出力ポートと接続されている。単一アンテナ素子であれば、エミュレータ構成は１×８ＳＩＭＯであり、２つのアンテナ素子であれば２×８ＭＩＭＯを有する構成等となる。
ＭＳ（ＤＵＴ）アンテナ特性は、未知であると仮定する。言い換えると、この情報は、このＯＴＡモデリングでは使用されない。
【００２４】
テストチャンバは無響室(anechoic room)である。加入者端末等のＤＵＴ４００は、アンテナ素子４０２、４０４、４０６、４０８、４１０、４１２、４１４、及び４１６に囲まれており、これらのアンテナ素子はエミュレータ４１８に接続されており、エミュレータは例えばＥＢ（Ｅｌｅｋｔｒｏｂｉｔ）Ｐｒｏｐｓｉｍ（登録商標）Ｃ８である。エミュレータ４１８は、プロセッサ、メモリ、及び適宜のコンピュータプログラムを備えている。この例では、８つのアンテナ素子が、一定の角度４５°で分離されて円形状に存在する。通常は、少なくとも２つのアンテナ素子４０２〜４１６が存在し、分離角度Δθで互いから分離している。少なくとも３つのアンテナ素子４０２〜４１６が存在する場合、任意の連続する２つのアンテナ素子４０２〜４１６に関し、分離角度Δθは同一であるか又は相違していてもよい。アンテナ素子４０２〜４１６は、ＤＵＴ４００から同一距離又は異なる距離にあり、かつアンテナ素子４０２〜４１６は、全角度又は全立体角度で配置される代わりに、セクターのみに配置されてもよい。ＤＵＴ４００も、アンテナ内に１又は複数の素子を備えている。
【００２５】
ＤＵＴ４００と空中無線で通信することにより、試験にパス方向が自由に含まれるような方法で、アンテナ設計、偏波及び配置効果の試験をすることができる。これは、エミュレータ４１８及びＤＵＴ４００間でケーブル接続が用いられる場合には不可能である。
エミュレータ４１８は、試験のためのチャネルモデルを有している。チャネルモデルは、試験を行う人間によって選択される。さらに、干渉及びノイズは、希望の方法及び程度で試験に入力される。使用されるチャネルモデルは、現実の無線システムから記録されたチャネルに基づくプレイバックモデル、又は人工的に生成されたモデル、又はプレイバックモデル及び人工的に生成されたモデルの組み合わせである。
【００２６】
エミュレータ４１８が無線システムの基地局と接続されている又は基地局として動作していると仮定し、アンテナ素子４０２〜４１６はＤＵＴ４００に送信し、ＤＵＴ４００は無線チャネルの受信側の加入者端末として動作するものとする。ＤＵＴのアンテナ特性は未知であり、かつその情報は下記の例では無視されるものとする。ＯＴＡアンテナ素子４０２〜４１６はＤＵＴから角度θｋの方向にあるものとし、ここでｋは１，．．．，Ｋであり、Ｋはアンテナ素子の数である。アンテナ素子４０２〜４１６の角度広がりは一定で、θｋ＋１−θｋ＝Δθである。
【００２７】
エミュレータ４１８内にある形状モデルが、ＯＴＡアンテナ素子４０２〜４１６にマッピングされる。エミュレータ４１８は、基地局からの送信放射線がクラスタに当たる状況をシミュレートする。エミュレータ４１８はさらに、各クラスタからの反射及び／又は散乱ビームを形成し、かつ、クラスタの発射電力及び遅延を、少なくとも１つのアンテナ素子４０２〜４１６に適するよう分割する。したがって、アンテナ素子４０２〜４１６は、クラスタの反射及び／又は散乱ビームを再生するよう制御される。
クラスタからの反射及び／又は散乱ビームを表しているビームの角度が、アンテナ素子４０２〜４１６の角度θｋと、例えば１°等の閾値を超えて異なっている場合がよくある。その場合、ビームは少なくとも２つのアンテナ素子４０２〜４１６を用いて送信される。
【００２８】
一実施形態では、シミュレートされたクラスタの電力は、アンテナ角度θｋ及びクラスタ角度ψｎに基づいて、２つのアンテナ素子間で分割される。クラスタ角度ψｎに最も近接しているアンテナ素子ｋの角度θｋは、次の数式によって得られる。
【数５】
0005358021.tif 000007
ここで、ｍｉｎはθｊの全ての値の中での式の最小値を意味し、ｉｎｔは分割の整数値（０を含む）を意味する。ｋの値は、
【数６】
0005358021.tif 000008
である。すると、第２のアンテナ素子ｋ＋１は、角度θｋ＋Δθ＝θｋ＋１を有するものになる。したがって、選択されたアンテナ素子は、それらの間で少なくとも主として１つのクラスタから反射及び／又は散乱されたビームがＤＵＴ４００に向けられているアンテナ素子である。
【００２９】
クラスタｎのためのＯＴＡアンテナの選択は、ＡｏＡの公称の角度ψｎに最も近い２つのθｋの値を選択することによって行われる。クラスタｎの電力は、θｋ及びψｎ間の距離によって、２つのＯＴＡアンテナ間で分割される。例えば、ψｎがθｋ及びθｋ＋１のちょうど真ん中にあれば、電力は５０パーセントずつになるよう分割される。
各アンテナ素子４０２〜４１６の重みＷｎ，ｋは、次のように算出される。
【数７】
0005358021.tif 000009
上記において、ｉは１及び２のいずれかであり、ｋはクラスタｎの角度ψｎに最も近いアンテナ素子のインデクスである。クラスタｎのアンテナ素子ｋへの電力Ｐｎは、重みＷｎ，ｋと掛け合わされ、Ｐｋ＋Ｐｋ＋１＝Ｐｎになるようにする。
【００３０】
図１０は、ＯＴＡチャンバ内のＡｏＡを示している。線１０００はＡｏＡベクトル、円はＤＵＴ４００の周囲のＯＴＡアンテナ素子である。
図１１は、図１０のアンテナ素子のアンテナ重み付けを示している。曲線１１００は、受信機から見たアンテナ素子のビームを表す。アンテナ素子の２つの重み１１０２及び１１０４はノンゼロであるが、残りのものはゼロである。
【００３１】
ＤＵＴの周りに円形状にあるアンテナ素子は８つであり、すなわちＫ＝８かつΔθ＝４５°であり、単一の基地局、単一のクラスタ、クラスタ電力２、ＡｏＡψｎ＝３７°であると仮定する。アンテナ素子４０２（アンテナｋ）のための電力Ｐｋ、及びアンテナ素子４０４（アンテナｋ＋１）のための電力は、以下のように表される
【数８】
0005358021.tif 000010
【００３２】
図５は、算出された電力分配を有するアンテナ素子４０２、４０４によって形成されたビーム５００を示している。異なるアンテナ素子へ与えられた信号もまた、互いに対して位相シフトされて、方向電力スペクトルが修正されるようにする。位相シフトは、ベースバンド信号を、信号の電力及び関連する遅延を設定する適宜の複素係数(complex coefficients)で重み付けすることにより形成される。位相シフトは、互いに対して無線周波数信号を遅延させることによっても実行される。例えば、所望の遅延は、デジタル遅延のバンク（例えば、デジタル有限インパルス応答フィルタ構造）から選択される。疑似無線チャネルの様々なパスの様々なビームが、様々な時点で形成される。疑似無線チャネルの１つのパスの１つのビームが、様々な時点で形成されてもよい。疑似無線チャネルの様々なパスの複数の様々なビームが、１つの時点で形成されてもよい。
【００３３】
図６は、アンテナ素子群６００を示している。一実施形態では、アンテナは、少なくとも１群６００のアンテナ素子６００２、６００４、６００６、６００８、６０１０からなる。このように、例えばアンテナ素子４０２の代わりに、１つのアンテナ素子だけでなく、幾つかの素子６００２、６００４、６００６、６００８、６０１０が存在していてもよい。アンテナ素子群は、例えば５つの素子を備えている。通常は、１つのアンテナ素子４０２〜４１６の代わりに、一群６００の少なくとも２つのアンテナ素子６００２、６００４、６００６、６００８、６０１０が存在する。
【００３４】
ＯＴＡアンテナ素子へのマッピングは、単一のＯＴＡアンテナ素子が一群６００のアンテナ素子６００２、６００４、６００６、６００８、６０１０に置き換えられた場合は、より単純かつ正確である。１つの群がＮ個のアンテナ素子６００２、６００４、６００６、６００８、６０１０を備えていると仮定する。
各アンテナ群６００に与えられる素子６００２、６００４、６００６、６００８、６０１０の数は、チャネルモデルの（クラスタ毎の）到来アジマス広がりに基づいて選択される。各群は単一のエミュレータ出力ポートによって供給され、かつ各群のアンテナ素子６００２、６００４、６００６、６００８、６０１０は、少なくとも１つのスプリッタ、コンバイナ、減衰器及び／又は移相器を備えたスイッチングネットワーク６２０によってエミュレータと接続されている。一実施形態では、スイッチング（すなわち、アンテナ素子の選択）は、全ての群で同様であり、一測定につき一回のみ実行される。
【００３５】
エミュレータからの信号に基づいて、ビームコントローラ６２２がビームに対して１つの群に幾つのアンテナ素子が必要であるかを制御する。通常は、最大数までの任意の正数のアンテナ素子が使用される。
一実施形態では、奇数の素子が使用される。例えばＮ＝５であれば、選択肢は１、３又は５であり、これはチャネルモデルのシナリオに依存する。チャネルモデルにあるのが狭いクラスタであれば、ビームに対して３つの素子があれば充分である。クラスタが幅広い場合は、ビームに対して最大数の素子が使用される。
【００３６】
一つの群内のアンテナ素子の選択は、次の数式で表される。
【数９】
0005358021.tif 000011
上記式に置いて、ここでＺ＝Ｎ−２ｊであってｊは０，．．．，（Ｎ−３）／２を意味し、roundは除算の値を最も近い整数値まで丸めること（最小値は１）を意味している。
【００３７】
チャネルモデルのＯＴＡアンテナへのマッピングは、次のルールを適用することによって実行される。クラスタの公称の方向によって、クラスタのそれぞれを、適宜のエミュレータチャネル及びＯＴＡアンテナ素子に設定する。クラスタｎに対するＯＴＡアンテナ素子の選択は、最も近いＯＴＡアンテナ群の中心θｋをクラスタの公称のＡｏＡψｎとすることによって達成される。スイッチ６２２によって一群の中のアンテナ素子の数、例えばＺ’、を選択する。
【００３８】
図７は、アンテナ素子群６００〜６１４に囲まれているＤＵＴ４００を示している。この例では、６００〜６１４の各群は、３つのアンテナ素子を有している。ビーム７００はグループ６０２を用いて形成される。群が８つで各群の素子が５つであれば、円の全体が均一に配置されたアンテナ素子で覆われる。クラスタが極度に幅広くて、例えばΔθよりも広いような非常に大きなビームを必要とする場合は、クラスタは２以上のアンテナ群にマッピングされる。
【００３９】
１つのビームを形成するために幾つかの群が使用されることもある。群は、２つのアンテナ素子を選択するための方程式（４）及び（５）に関連して説明されたものと同じやり方で適用される。そして１つのビームに対して、２つのアンテナ素子を選択する代わりに、２つのアンテナ素子群が選択される。図７では、ビーム７００が群６００及び６０２を使用して形成される。
一実施形態では、固定重みがアンテナ素子にインプリメントされて、ガウス又はラプラシアン形のクラスタ電力アジマススペクトルが複製されるようにする。
【００４０】
少なくとも２つのアンテナ素子を用いた受信は、同様の方式で実行される。したがって、この方法はアップリンク及びダウンリンクの両方に適用される。アンテナ素子４０２〜４１６がＤＵＴ４００からの信号を受信していると仮定する。少なくとも２つのアンテナ
素子４０２〜４１６によって受信された信号は、エミュレータ４１８で結合されて、１つの疑似無線チャネルのパスの信号の受信ビームを形成する。この結合は、式（４）及び（５）で計算される重みＷｎｋ＋１を用いて、２つのアンテナ素子又はアンテナ素子群からの電力を重み付けすることを含んでいる。さらに、ビームの形及び方向は、複素係数又は別の種類の位相シフトを用いて重み付けされる。
【００４１】
このような実施形態は、３ＧＰＰ（Third Generation Partnership Project：スリージーピーピー）のＬＴＥ（Long Term Evolution：エルティーイー）、ＷｉＭＡＸ（Worldwide Interoperability for Microwave Access：ワイマックス）、Ｗｉ−Ｆｉ（ワイファイ）、及び／又はＷＣＤＭＡ（Wide-band Code Division Multiple Access：ダブリュシーディーエムエー）に適用することができる。ＭＩＭＯ（Multiple In Multiple Out：マイモ）もまた適用可能であるが、ＭＩＭＯでは本実施形態とは異なる方式で、信号がアンテナ素子に分配される。図８は、２つの送信アンテナ素子８００、８０２、及び２つの受信アンテナ素子８０４、８０６を有するＭＩＭＯ構成を示している。エミュレータ８１２の遅延素子８１４〜８２０には、異なるパスを示している２つの遅延タップ８０８、８１０がある。各送信アンテナ８００、８０２からの信号は、同一の遅延（タップ８０８、８１０）で信号を遅延させる遅延素子８１４〜８２０に与えられる。２つの遅延（タップ８０８、８１０）で遅延している遅延素子８１４及び８２０の出力は、結合されてアンテナ素子８０６に与えられる。同様に、両方の遅延（遅延タップ８０８、８１０）で遅延している遅延素子８１６及び８１８の出力も、結合されてアンテナ素子８０４に与えられる。
【００４２】
図９は、本発明の実施形態の一例を示している。この例でも、２つの送信アンテナ素子９００、９０２、及びＯＴＡ試験の無響チャンバ９２２内にある複数のアンテナ素子の内の２つの受信アンテナ素子９０４、９０６がある。２つの遅延タップ９０８、９１０はエミュレータ９１２の遅延素子９１４〜９２０における異なるパスを示している。送信アンテナ９００からの信号は、遅延素子９１４、９１６に与えられる。遅延素子９１４は、遅延タップ９０８に相当する遅延で信号を遅延させ、かつ遅延素子９１６は、遅延タップ９１０に相当する遅延で信号を遅延させる。
【００４３】
送信アンテナ９０２からの信号は、遅延素子９１８、９２０に与えられる。遅延素子９１８は、遅延タップ９１０に相当する遅延で信号を遅延させ、かつ遅延素子９２０は、遅延タップ９０８に相当する遅延で信号を遅延させる。同一の遅延（遅延タップ９０８）で遅延している遅延素子９１４及び９２０の出力は、結合されてアンテナ素子９０６に与えられる。同様に、同一の遅延（遅延タップ９１０）で遅延している遅延素子９１６及び９１８の出力は、結合されてアンテナ素子９０４に与えられる。このように、異なる遅延タップが異なるＡｏＡを表している場合は、異なる遅延タップが異なるアンテナ素子９０４、９０６に供給される。
【００４４】
ＯＴＡチャンバ内部の空間効果の生成は、正弦波の和をベースにしたチャネルモデリングに対応している。時空間チャネルモデルのためのパラメータを算出する技術は、Ｌｐノルム法（Lｐ-norm method）と称され、ＯＴＡチャネルモデリングのために精密化されている。正確な空間相関モデリングのためには、Ｌ２−ｎｏｒｍ  Ｅρ（ｇ1，ｇ2，．．．，ｇｋ）等の費用関数が最適化される。
【数１０】
0005358021.tif 000012
上記式において、ρ（Δｍ，ψ０，σψ）は、アンテナ素子のアンテナ素子間隔Δｍに関する論理空間相互相関であり、ψ０は公称のＡｏＡ、σψは角度広がりであり、ρ（Δｍ）（ただし、ρの上部に〜が付いている）は、実際のＯＴＡアンテナ素子で取得された
空間相互相関である。そして、項Δｍはｄｓ／λ０と定義され、λ０は無線信号の電磁波の波長である。Ｌｐノルム最適化の代わりに、何らかの他のノルムが使用されてもよく、最適化は何らかの他の費用関数に基づいていてもよい。
【００４５】
最適化の目的は、利得係数(gain factor)Ｇに対する上述の費用関数を最小化すること、又は直接に重みｗに対する費用関数を最小化することによって、重みｗのＯＴＡアンテナ素子を決定することである。あるいは、勾配法、半空間法（half space method）等によって、最適化が実行されてもよい。
相互相関ρは、利得係数ｇの関数として表すことができるのが通常である。
ｆ（ｇ1，ｇ2，．．．，ｇｋ）＝ρ                  （９）
ここでｆは関数である。同様に、項Ｇは相互相関ρの関数として表される。
Ｇ＝（ｇ1，ｇ2，．．．，ｇｋ）＝ｆ−１（ρ）      （１０）
なお、ｆ−１はｆの逆関数である。
【００４６】
ラプラシアン形のＰＡＳ（電力角度スペクトル）の論理相互相関関数ρ（Δｍ，ψ０，σψ）は、次のように定義される。
【数１１】
0005358021.tif 000013
実際には、トランケートされたラプラシアンＰＡＳのために又は離散近似によって、算出される。ＯＴＡアンテナ素子で取得された空間相関は、次のように定義される。
【数１２】
0005358021.tif 000014
【００４７】
上記式において、項目Ｋ’は、ＯＴＡチャンバ内の複数のアンテナ素子の中のアクティブなアンテナ素子の数を表している。アクティブなアンテナ素子が希望のビームを形成する。空間アンテナの８つのＯＴＡ素子のソリューションでは、Ｋ’＝３、θｋ∈｛０°，４５°，９０°，１３５°，１８０°，２２５°，２７０°，３１５°，３６０°｝であり、ｇｋはｇｋ⊂[０，１]であるように制限される。Ｍの実用的な値は約５０である。最適化は、Ｋ’次元空間に二分探索を適用することによって数値的に実行されるが、それは方程式（８）が凸関数であるからである。二分探索を用いれば、必要になるのは、約ｌｏｇ２ＬＫ’＝Ｋ’ｌｏｇ２Ｌの反復（すなわち、方程式（８）の計算）のみとなり、ここでＬはｇｋ⊂[０，１]のポイントの数である。例えばＬ＝１０００、Ｋ’＝３であれば、３０反復のみが必要とされる。これらのパラメータを用いた場合、ブルートフォース法（力ずくの方法）では、方程式（８）に対して１０００３＝１０９のソリューションが必要になる。
【００４８】
式（８）は、（１１）及び（１２）を適用して、勾配法及び半空間法等の数値最適化方法を用いることによって計算される。
表記を単純化するために、重みをベクトルＧで以下のように表す。
Ｇ＝（ｇ1，ｇ2，．．．，ｇｋ）                  （１３）
そして、一連の位相の項をベクトルＡｍとする。
Ａｍ＝（aｍ１，aｍ２，．．．，aｍＫ’）
＝(exp（-j2πΔｍsinθ１）, exp（-j2πΔｍsinθ２），
．．．，
exp（-j2πΔｍsinθＫ’））                  （１４）
さらに、論理相互相関をスカラーρｍとする。
ρｍ＝ρ（Δｍ，ψ０，σψ）                      （１５）
【００４９】
Ｅρは勾配＝ゼロを解くことによって最小化される。
【数１３】
0005358021.tif 000015
ここでＵｋは第ｋ番目のユニットベースベクトルである。上記の勾配式は、一組のＫ’の方程式に加工されてもよく、重みｇｋに関して解かれる。
【数１４】
0005358021.tif 000016
式（１７）は、一組の解析的な式、すなわち勾配の式（１６）をゼロにするものを表している。
【００５０】
正確な相関モデリングには、下記のパラメータによるＬ２−ｎｏｒｍの最適化が用いられる。まず第１に、公称の到来角ψ０に最も近いＫ’のアンテナ素子が、第１のＫ’のインデクスｋｉを次のシーケンスから取り出すことによって探し出される。
|exp(jψ０)-exp(jθｋ１)|<|exp（jψ０）-exp(jθｋ2)|<  ・・・・・・
<|exp(jψ０)-exp(jθｋK)|
アンテナ素子ｋ１が、公称の到来角ψ０に最も近い。そして第２に、公称の到来角ψ０が、式（８）及び（９）でゼロに設定されて、空間相互相関ρ（Δｍ）（ただし、ρの上に〜が付いている）が形成される。
【数１５】
0005358021.tif 000017
【００５１】
換言すると、正確な相関モデリングのためには、式（８）及び（１１）における公称の到来角ψ０は、最も近いＯＴＡアンテナ素子の方向θｋｉに丸められる。それから、Ｋ’個のアクティブなＯＴＡアンテナ素子（奇数個のアンテナ）がアンテナｋｉの周囲に対称的に選択される。大体の目安として、数Ｋ’はＫ’Δθ＜１８０°であるべきである。例えば、Ｋ’＝３、Δθ＝４５°であれば、式（１２）にはアンテナ素子角度θｋ１＝−４５°、θｋ１＝０°及びθｋ１＝４５°が設定される。ここで係数｛ｇｋ−１，ｇｋ，ｇｋ＋１｝は式（８）を最小化することによって決定される。他の係数ｇｋはゼロと等しい。例えば、図１０では、アンテナのインデクスは３、２、及び４であり、対応する方向θＫは９０°、４５°及び１３５°である。
【００５２】
最後に、クラスタｎそして同時にアンテナ素子ｋの重みであるＷｎ，ｋが以下のように表される。
Ｗｎ，ｋ＝√ｇｋ                  （１９）
なお、係数ｇｋは各クラスタｎに対して個別に決定される。
【００５３】
図１２は、３つのアンテナ素子の論理空間相関１２００、及び８つのＯＴＡアンテナ素子を４５°の間隔で有している３５°ラプラシアンＰＡＳの理想的な空間相関１２０２の例を示している。
ＯＴＡチャンバでは、アンテナ素子の位置は固定されている。任意の到来角（ＡｏＡ）をモデリングする場合は、ＯＴＡアンテナ間の方向は、補間される必要がある。これは、丸めをすることなく実際の公称の到来角ψ０を用いることで、式（６）のノルムを最小化することによって行われる。そうでなければ、手順は上述の通りである。
【００５４】
図１３及び１４の例では、アンテナ素子の重みは同時最適化によって決定された。図１３は、アンテナ素子の重み１３０２、１３０４、１３０６、１３０８及び１３１０、並びに結果として生じるＰＡＳ１３００を表している。図１４は、３つのアンテナ素子の理論空間相関１４００、及び理想的な空間相関１４０２を表している。この例は、８つのアンテナを４５°の間隔で有しているＲＸ上の２−素子ＵＬＡを示しており、ＡｏＡ＝１００°である。ターゲットＡｏＡは１００°であり、波長分離０．５を有するターゲット相互相関ρ（Δｍ，ψ０，ρψ）は、｜ρ｜＝０．２４７６になる。結果として生じる相関マトリクスＲｒｘ＿ａｂｓが以下に記されており、結果として生じた図１３におけるＰＡＳの最大値は１０１°である。
【数１６】
0005358021.tif 000018
【００５５】
無線チャネルのエミュレーションでは、チャネルインパルス応答がエミュレータに与えられて、かつ送信信号で畳み込みされる。畳み込みのエミュレーションでは、異なるＭＩＭＯチャネル（Ｔｘ／Ｒｘアンテナ対）のインパルス応答は、等しい電力遅延プロフィル及びタップ数を有している。ＯＴＡ環境のためのモデリングは相違する。チャネルインパルス応答はディスアセンブルされて、クラスタ（タップ）のＡｏＡ情報に基づいて、異なるＯＴＡアンテナ素子に対してリアセンブルされる。ＳＣＭモデルの実現例の６つのクラスタ１５００、１５０２、１５０４、１５０６、１５０８、及び１５１０の元々のＰＤＰ（電力遅延プロファイル）を図１５に示す。ＯＴＡエミュレーションのケースにおける、８つのチャネルの遅延タップのマッピングを図１６に示す。６つのクラスタがあり、それぞれが異なる遅延を有している。
【００５６】
電力重み付けでＯＴＡアンテナへタップをマッピングすることに加えて、元々のフェージング信号も、ドップラシフティングによって修正される必要がある。これは希望の相関及びＡｏＡ効果を得るために必要なことである。各形状チャネルモデルでは、移動中のモバイル端末が想定されている。端末の動きは、移動角θｖの特定方向を有する速度ベクトルによって表される。
平面波が、方向ψｎの代わりにアンテナ素子ｋの方向θｋを有する場合は、式（２）は以下のようになる。
【数１７】
0005358021.tif 000019
【００５７】
そして、ＯＴＡアンテナ素子ｋ及びクラスタｎに対するドップラ補正項Ｃの結果は、以下のようになる。
Ｃｋ，ｎ＝νｋ−νｎ              （２１）
最後に、重みＷｎ, ｋによる電力重み付けに加えて、ＯＴＡアンテナ素子ｋによって送
信された、クラスタｎの疑似無線チャネルＨｓ，ｎ（ｔ，τ）は、例えば次のような方法でドップラ補正項Ｃによってドップラシフトされる。
【数１８】
0005358021.tif 000020
上記式に置いて、Ｈｓ，ｎ（ｔ，τ）は、式（１）、（３）及び（４）に表されたインパルス応答Ｈｕ，ｓ，ｎ（ｔ，τ）及びＨｓ，ｎ（ｔ，τ）に基づいた、最初の疑似無線チャネルである。疑似無線チャネル
【数１９】
0005358021.tif 000021
が今度は、重みＷｎ，ｋを用いてアンテナ素子４０２〜４１６及びＤＵＴ４００に対してシフトされる予定の又は既にシフトされた、アンテナ素子ｋの疑似無線チャネルを表す。さらに、疑似無線チャネル
【数２０】
0005358021.tif 000022
も、項ｅｘｐ（ｊ２πＣｋ，ｎｔ）によってドップラシフトされた可能性のあるアンテナ素子ｋの疑似無線チャネルを表す。
【００５８】
ＯＴＡチャネルモデリングの正確さも考慮に入れられる。図１７は、ＤＵＴの空間分解能が２４°であり、ＯＴＡアンテナ素子の数は８であり、かつアンテナ素子の間隔は４５°である状況を示している。ＯＴＡアンテナ素子は円でマーキングされている。基準数１７００はＡｏＡベクトルを意味し、基準数１７０２は速度ベクトルを意味している。
図１８は、図１７に示されている状況において受信機から見たＰＡＳ１８００を表している。円１８０２はアンテナ素子の相対電力を示している。ＰＡＳ１８００は２つのピークを有しているが、これは望ましくない。
【００５９】
図１９は、ＤＵＴの空間分解能が２４°であり、ＯＴＡアンテナ素子の数は１６であり、アンテナ素子の間隔は２２．５°である状況を示している。ＯＴＡアンテナ素子は円でマーキングされている。基準数１９００はＡｏＡベクトルを意味しており、基準数１９０２は速度ベクトルを意味している。
図２０は、図１９に示されている状況において受信機から見たＰＡＳ２０００を表している。円２００２は、アンテナ素子の相対電力を示している。ＰＡＳ２０００は１つのピークのみを有しており、かつそれが望ましい。ＤＵＴのアンテナアレイサイズが空間分解能を決定する。大体の目安として、λ／２ＵＬＡの分解能は９６°／Ｍであり、ここでＭはＤＵＴのアンテナの数である。例えば、２アンテナＵＬＡ、すなわちＭ＝２であれば４８°ＡｏＡという結果になり、４アンテナＵＬＡは２４°ＡｏＡという結果になる。したがって、ＯＴＡアンテナ素子間の間隔は、ＤＵＴの空間分解能より小さいことが望ましい。
【００６０】
ラプラシアン形ＰＡＳ及び３５°ｒｍｓのアジマス広がりを用いて、１波長サイズであってＯＴＡチャンバの8つのアンテナ素子を有するアレイ、並びに２波長サイズであってＯＴＡチャンバの１６のアンテナ素子を有するアレイを制御することが可能である。
クラスタＰＡＳモデリングに用いられるＯＴＡアンテナの数が、どのくらいの大きさのＤＵＴのアレイサイズが正確な相関を有することができるかを決定する。ＤＵＴのサイズは小型の方がよいが、ＯＴＡ内により多くのアンテナ素子があればそれだけ、ＤＵＴはより大きなサイズを有することができる。
【００６１】
図２１は、５つの送信アンテナ素子を有するＰＡＳ２１００、及び３つの送信アンテナ
素子を有するＰＡＳ２１０２を表している。
ドップラ及びＯＴＡのアンテナ素子の可能な相関、さらにチャネル遅延プロファイルを含んだフェージングが、チャネル係数には含まれている。
チャネル係数は式（１８）の修正版である以下の式によって生成される。
【数２１】
0005358021.tif 000023
ここでｋは、ｋ＝２π／λ０と定義される波数である。
【００６２】
ＯＴＡチャンバがデュアル偏波アンテナ素子を有する場合は、チャネル係数の式は、Ｖ偏波及びＨ偏波で別々に表される。
【数２２】
0005358021.tif 000024
上記式において、Ｆｔｘ，ｓＶ、Ｆｔｘ，ｓＨ  はそれぞれ、アンテナ素子のＶ（垂直）偏波及びＨ（水平）偏波のフィールドパターンである。位相項Φｎ，ｍｖｖ等は、ランダム初期位相∈［０，２π］であり、κｎ，ｍはｋ交差偏波電力比（ＸＰＲ）である。
【００６３】
しかし、ドップラ周波数νｎ，ｍはやはりＡｏＡ角に基づいて決定される。生じるものは、マトリクス係数Ｈｎ（ｔ，τ）を有する離散インパルス応答である。Ｈｎ（ｔ，τ）の次元は、単一偏波の場合は１×Ｓであり、デュアル偏波の場合は２×Ｓであり、Ｓは基地局内のアンテナ素子の数である。このステップは、例えばＳＣＭＥ又はＷＩＮＮＥＲモデル等の形状チャネルモデルのマトラボ実装によって行われる。
【００６４】
次に、クラスタの公称の方向及びクラスタの角度広がりに依存して、適宜のエミュレータチャネル及びＯＴＡアンテナ素子に対して、クラスタｎのマッピングが実行される。選択された方法は、正確なＡｏＡ（方程式（５）を参照）、正確な空間相関（方程式（１８）を参照）、及び両方のバランスのとれた組み合わせ（図１４及びその説明を参照）のいずれが重視されるかに依存している。単一偏波の場合の方法が説明されているが、デュアル偏波の場合にも適用することが可能である。唯一の違いは、デュアル偏波の場合には、方程式（２４）及び（２５）によるＶ（垂直）及びＨ（水平）に偏波されたチャネルインパルス応答が、Ｖ偏波及びＨ偏波されたＯＴＡアンテナ素子に対して別々にマッピングされることである。
【００６５】
図２２は、８つの均一な間隔にあるデュアル偏波されたチャンバアンテナ素子２２０２，２２０４，２２０６，２２０８，２２１０，２２１２，２２１４，及び２２１６を有するＯＴＡチャンバアンテナの設定を表している。図２２では、Ｖ偏波されたＯＴＡアンテナ素子は紙面（アジマス面）と実際は垂直である。デュアル偏波されたＯＴＡアンテナ素子が図２２のように用いられる場合は、エミュレータ構成は、１つの基地局出力信号１×
１６ＳＩＭＯを有するもの、２つの基地局出力信号２×１６ＳＩＭＯするもの等である。例えば、アンテナＡ１Ｖは、第１のＯＴＡアンテナ位置にある垂直（Ｖ）偏波された素子であることを意味し、Ａ８Ｈは、第８のＯＴＡアンテナ位置にある水平（Ｈ）偏波された素子を意味している。
上記で説明されたことは、ＯＴＡチャンバ内のアンテナ素子を有する疑似無線チャネルの、少なくとも１つのアンテナパターンを形成する方法に関している。次に、例えば、ＤＵＴ及びアンテナ素子に対してどのようにして各アンテナパターンを回転するかを検証する。
【００６６】
図２３及び２４は、単純化された方法による、疑似無線チャネルのシフトを示している。図２３は、エミュレータ４１８が１つの送信機又は受信機のみをシミュレーションし、１つのビーム５００のみが形成される一実施形態を示している。この例では、エミュレータ４１８は、スイッチ２３００、コントローラ２３０２、及びメモリ２３０４を有している。基地局からの信号を表す入力信号が、スイッチ２３００に供給される。エミュレータ４１８のコントローラ２３０２は、スイッチ２３００の状態を制御することにより、推定無線チャネルのシフトを制御する。重み付けは、振幅及び位相のシフトを含んでいる。スイッチ２３００を適切に制御することにより、アンテナ素子４０２〜４１６に供給された疑似無線チャネルの信号によって形成されたビーム５００が、ＤＵＴ５００の周囲をアンテナ素子４０２〜４１６の周囲に沿って電気的に回転される。図２３の例では、スイッチＡ及びＨは、アンテナ素子４１０及び４１２にＲＦ信号を供給するために閉じた状態（オン状態）に切り替えられており、スイッチＢ〜Ｇは開いた状態（オフ状態）に切り替えられて通信のある時点でのビーム形成にアクティブに参加しないようにされる。アンテナ素子４１０及び４１２に供給されたＲＦ信号は、ＤＵＴ４００に向けられたビームを形成するよう互いに干渉する。
【００６７】
図２４は、別の時点のエミュレータ４１８及びスイッチ２３００を示しており、その時点は、通信の時点の時系列において、図２３に示されている時点のすぐ次の時点である。図２３では、アンテナパターンを決定する閉じられた状態のスイッチは、スイッチＡ及びＧである。図２４における閉じられた状態のスイッチは、１インクリメントの回転により、スイッチＧ及びＨである。図２３及び２４に示された時点の間におけるスイッチＡ〜Ｈの状態の変化は、疑似無線チャネルを、アンテナ素子４０２〜４１６に対して形状的にシフトする。図２３及び２４に示された時点の間におけるスイッチＡ〜Ｈの状態の変化はさらに、２つのアンテナ素子４０２〜４１６間の距離ｄｓに対応する離散ステップによって、ビーム５００のシフトを生じさせる。したがって、ビーム５００の方向２４００は、２つのアンテナ素子４０２〜４１６の角度Δθと同一である角度Δαだけ離散的に回転する。ビーム５００の角度Δαの回転は、非常に頻繁に生じるので、ビーム５００はＤＵＴ４００の周囲を少なくとも一度は完全に回転する。
【００６８】
開いた状態のスイッチは、ＲＦ信号をアンテナ素子４０２〜４１６へ又はアンテナ素子４０２〜４１６からＲＦ信号へ、ゼロの重み付けをすることに相当している。そして閉じられた状態のスイッチは、ＲＦ信号をアンテナ素子４０２〜４１６へ又はアンテナ素子４０２〜４１６からＲＦ信号へ、実際の重みの重み付けをすることに相当している。したがって、スイッチＡ〜Ｈを全く必要とせずに、各スイッチＡ〜Ｈが、ＲＦ信号を実際の重みで乗算するマルチプライアに置き換えられる場合もある。次にコントローラ２３０２が、時間に応じてアンテナ素子４０２〜４１６の１つのアンテナ素子から次のアンテナ素子へ重みをシフトして、ＤＵＴ４００の試験を１つのアンテナパターンで２以上の方向から実行できるようにする。つまり、疑似無線チャネルはアンテナ素子４０２〜４１６に対して新しい位置に設定され、それは各重みを１つのアンテナ素子から別のアンテナ素子へとシフトすることによって行われる。
【００６９】
図２５は、エミュレータ４１８が２つの送信機又は受信機をシミュレーションする実施形態を示している。第１の信号Ｓ１及び第２の信号Ｓ２がエミュレータ４１８に供給される。信号Ｓ１及びＳ２は、例えば２つの異なる基地局の信号である。あるいは、エミュレータが信号Ｓ１及びＳ２を出力し、その後信号は、例えば２つの異なる基地局によって受信される。ＤＵＴ４００に向けて送信する場合、信号Ｓ１及びＳ２のそれぞれがＦＩＲフィルタ２５００〜２５３０に供給され、これらがチャネルモデルに従って信号を歪める。ＤＵＴから受信する場合は、信号Ｓ１及びＳ２のそれぞれが、チャネルモデルに従って信号を歪めたＦＩＲフィルタ２５００〜２５３０から受信される。信号Ｓ１、Ｓ２の数は３以上であることが通常である。異なる基地局の信号Ｓ１、Ｓ２は、別々のＦＩＲフィルタに供給されて、信号が独特の方法で歪められるようにする。少なくとも１つのビーム５００のアンテナ素子４０２〜４１６に対するシフトは、例えば、ＦＩＲフィルタ２５００〜２５３０における重み係数を変更することによって実行される。
【００７０】
図２６は、アナログ／デジタルコンバータ２６００、重み付け素子２６０２、シフトレジスタとして配置されている遅延素子２６０４、マルチプライア２６０６、サマー２６０８、ドップラ素子２６１０及びデジタル−アナログコンバータ２６１２を備えているＦＩＲフィルタのブロック図を示している。アナログ／デジタルコンバータ２６００は、アナログ信号Ｓ1又はＳ２を受信する。重み付け素子２６０２及びドップラ素子２６１０を有していないＦＩＲフィルタの基本的な機能は、次の通りである。アナログ／デジタルコンバータ２６００からのデジタル入力信号Ｘ（ｎ）は、各遅延素子２６０４において遅延され、これらの遅延は同一又は異なる長さの時間を有しており、遅延された信号はマルチプライア２６０６で希望のチャネル係数ｈｊ（ｉ）によって乗算される。なお、ｉ＝０，．．．，Ｎであり、ｊ＝1，．．．，Ｋは、図２５のＦＩＲフィルタ２５００〜２５３０を意味している。チャネルモデルはチャネル係数ｈｊ＝ｈ（０），．．．，ｈ（Ｎ）によって定義され、これらのチャネル係数は、無線チャネルのチャネル推定又はＦＩＲフィルタのタップ係数とも称される。信号成分は加算器２６０８で合算されて、かつ合算された信号はデジタル／アナログコンバータ２６１２においてアナログ形式に変換される。
【００７１】
一実施形態では、各ＦＩＲフィルタ２５００〜２５３０は、重み付け素子２６０２を備えている。重み付け素子２６０２は、全ての遅延信号成分が遅延前又は遅延後に重み付けされている間は、アナログ／デジタルコンバータ２６００及びデジタル／アナログコンバータ２６１２の間のどこかに位置している。重み付け素子２６０２は、疑似無線チャネルＨｎ，ｋを重みＷｎ，ｋで乗算して、積Ｗｎ，ｋＨｎ，ｋを与えるためのマルチプライアである。
【００７２】
一実施形態では、各ＦＩＲフィルタはさらに、ドップラ素子２６１０を備えている。ドップラ素子は、重み付けされた無線チャネルＷｎ，ｋＨｎ，ｋを、ドップラシフトｅｘｐ（ｊ２πＣｎ，ｋｔ）で乗算して、積Ｗｎ，ｋＨｎ，ｋｅｘｐ（ｊ２πＣｎ，ｋｔ）を与えるためのマルチプライアである。ドップラ素子２６１０は、全ての遅延信号成分が遅延前又は遅延後にドップラシフトされている間は、アナログ／デジタルコンバータ２６００及びデジタル／アナログコンバータ２６１２の間のいずれかの位置にある。
【００７３】
重み付け係数ｗは、式（２）、（６）、（８）、及び（１９）を用いて計算され、加算器２６０８によるチャネル信号出力は、重みＷｎ，ｋによって乗算される。ただし、重みはチャネル係数とも結合されて、遅延信号がＦＩＲフィルタにおいて積Ｗｎ，ｋｈｊによって乗算されるようにする。同様に、ドップラシフトｅｘｐ（ｊ２πＣｎ，ｋｔ）は、マルチプライア２６０６においてその乗算結果と結合されて、Ｗｎ，ｋｈｊｅｘｐ（ｊ２πＣｎ，ｋｔ）を形成する。項Ｗｎ，ｋｈｊｅｘｐ（ｊ２πＣｎ，ｋｔ）が合算されると、所望の出力Ｗｎ，ｋＨｎ，ｋｅｘｐ（ｊ２πＣｎ，ｋｔ）が得られる。ＦＩＲフィルタ２５００〜２５３０における重みＷｎ，ｋは変更されて、少なくとも１つのビーム５００が
生成されてＤＵＴ４００の周囲をシフトするようにする。疑似無線チャネルの回転は、式（２）、（６）、（８）、又は（１９）での最適化において角度ψｎ，ｍ、ψ０又はψｎを変えることにより、重みの計算に含まれる。そして、重み係数Ｈ＝［ｈ１，．．．，ｈｋ］が、変更のために実際の無線チャネルの特性を測定する方法と同様の方法で変えられる。
【００７４】
一般に、無線チャネル係数ｈは、実係数又は複素係数である。シミュレータは直交変調を用いており、ここでは信号が２つの部分に分割されるため、複素チャネル係数が必要とされることが多い。実信号部分Ｉ（同位相）は位相シフトされていないキャリアによって乗算され、虚信号部分Ｑ（直交位相）は、位相シフトされたキャリアによって乗算される。このように、信号ｘは、ｘ＝Ｉ＋ｊＱの形で表わされ、ここでＩは実信号部分、Ｑは虚信号部分、そしてｊは虚数単位である。
【００７５】
数式では、ＦＩＲフィルタの出力信号ｙ（ｎ）は、遅延信号及びチャネル係数の積の和の畳み込みとして表される。
【数２３】
0005358021.tif 000025
なお、＊は畳み込み演算を意味し、ｎは信号素子のインデクスを意味している。信号ｘ及びｙ、並びにチャネルインパルス応答推定ｈは、スカラー形式、ベクトル形式、又はマトリクス形式で処理される。
【００７６】
図２７は、アンテナ素子４０２〜４１６のアンテナパターンにおいて２つのビーム５００、２７００が存在する通信の１つの時点を表している。この例では、第１の時点において、ビーム２７００はアンテナ素子４１０によって形成され、ビーム５００はアンテナ素子４１４及び４１６によって形成されるものとする。
【００７７】
図２８は、図２７における時点に関する通信のすぐ次の時点を示している。同じビーム５００、２７００が、ＤＵＴ４００及びアンテナ素子４０２〜４１６に対して、エミュレータ４１８によってその次の位置に回転されている。ビームの動きは、疑似無線チャネルの角度広がりのシフトを表している。すなわち、疑似無線チャンネル自体が変更されたのではなく、ＤＵＴ４００及びアンテナ素子４０２〜４１６に対して回転されたことを意味する。アンテナパターン又は電力角度スペクトルは、アンテナ素子４０２〜４１６に対してシフトされてもよい。ビーム２５００は今やアンテナ素子４１０及び４１２の間にあり、かつそれらによって形成されている。ビーム５００も、同じ分量だけ同じ方向に移動されている。ビーム５００及び２７００のそれぞれを、２つのアンテナ素子４０２〜４１６の角度Δθのインクリメントによってシフトする代わりに、エミュレータ４１８がビーム５００、２５００のそれぞれを角度Δθより少ない又は多い分量だけシフトしてもよい。推定無線チャネルが、新しい１つのクラスタ及び新しい複数のパスを有する新しい無線チャネルに変更されると、その新しい推定無線チャネルは、同様の方法で、アンテナ素子４０２〜４１６及びＤＵＴ４００に対してシフトされる。
【００７８】
通常は、ビーム５００、２７００のそれぞれを、無響チャンバにおいて通信の異なる時点で異なる方向からＤＵＴ４００に向けるために、エミュレータ４１８は疑似無線チャネルを、複数のアンテナ素子４０２〜４１６に対してシフトさせる。この例では、疑似無線チャネルは、各ビームをアンテナ素子４０２〜４１６に対して所定の分量だけシフトし、これによってＤＵＴ４００に対してもシフトすることにより、アンテナ素子４０２〜４１６に対して新しい位置に設定される。設定及びシフトの数は２以上である。各ビームがＤＵＴ４００に向けられる順序は、時間的順序である必要はない。
【００７９】
エミュレータ４１８は、各時点で個別にアンテナ素子の重みを形成する。あるいは、エミュレータ４１８は、あらかじめ重みが記憶されているメモリ２３０４を有しており、かつエミュレータ４１８のコントローラ２３０２が、アンテナ素子４０２〜４１６に対して無線チャネルをシフトするために、エミュレータ４１８及びＤＵＴ４００が通信している間に全ての時点において新しい重みをメモリ２３０４から抽出してもよい。希望の重みは、エミュレータ４１８又は何らかの別のコンピュータで事前に算出されてメモリ２３０４に記憶される。
【００８０】
重みは、式（２）、（６）、（８）及び（１９）を用いて計算される。疑似無線チャネルの回転は、式（２）、（６）、（８）、又は（１９）との最適化において角度ψｎ，ｍ、ψ０又はψｎを変えることによって、重みの計算に含まれる。角度ψｎ，ｍ、ψ０又はψｎは、疑似無線チャネルの各アンテナパターンの新しい方向のための重みを計算するために、例えば２０°インクリメントされる。ここで疑似無線チャネルのアンテナパターンが、初期到来角１０°、１３０°、２５６°及び３００°を有する４つのビームを有するものとする。表２は、ＤＵＴ４００の周りを２０°のインクリメントで一周する回転を示している。
【表２】
0005358021.tif 000026
【００８１】
図２３では、メモリ２３０４は、通信の各時点に関するスイッチＡ〜Ｈの状態の制御のためのデータを記憶している。エミュレータ４１８及びＤＵＴ４００の間の通信中に、スイッチＡ〜Ｈを異なる時点で所望の状態に設定して疑似無線チャネルを回転するために、エミュレータ４１８のコントローラ２３０２は、メモリからデータを抽出する。所望のデータは、事前に計算されてメモリ２３０４に記憶されている。
【００８２】
図２９は、減衰器２９００〜２９３０及び移相器２９３２〜２９６２を介して、信号がアンテナ素子４０２〜４１６に分配される一実施形態を示している。疑似無線チャネルは、減衰器２９００〜２９３０のコントローラ２３０２から適宜の重み付けコマンドを受信し、そして、移相器２９３２〜２９６２のコントローラ２３０２から適宜の重み付けコマンドを受信することによって、形成かつ回転される。増幅器であってもよいが、減衰器２９００〜２９３０への重み付けコマンドは、減衰器の振幅を増大又は低減させる。移相器２９３２〜２９６２への重み付けコマンドは、移相器においてＲＦ信号の位相を進行又は遅延させる。移相器は、制御可能な遅延を有する部品であってもよい。アナログ移相器の動作は、非線形現象に基づいている。アナログ移相器は、制御電圧に応じてキャパシタンスを変更するもので、例えばバラクタダイオードである。デジタル的には、減衰器及び移相器は、１又は２つのマルチプライア、及び、実部が減衰（増幅）を決定して虚部が位相シフトを決定するような複素重みを用いて実行される。アンテナ素子４０２〜４１６への信号は、様々な重みを用いて位相シフトかつ減衰されて、ＤＵＴ４００のための各推定無線チャネルを形成かつシフトする。
【００８３】
正弦波の加算に基づくエミュレータでは、アンテナ素子４０２〜４１６のためのＲＦ信
号が、減衰器２９００〜２９３０及び移相器２９３２〜２９６２において形成され、ＲＦ信号は、所望のアンテナパターンを形成し、そのアンテナパターンをアンテナ素子４０２〜４１６に対してシフトするために、少なくとも２つのアンテナ素子４０２〜４１６を介して送信又は受信される。アンテナ素子４０２〜４１６は、ＤＵＴ４００のための少なくとも１つのビーム５００、２５００を形成する。この種の通信は、干渉によりＤＵＴ４００において空気中でフェージングをもたらす。アンテナ素子４０２〜４１６の全てが、少なくとも１つのビームを形成する時にアクティブである。アンテナ素子４０２〜４１６の一部のみがアクティブであることも可能である。先行技術における正弦波の加算のシミュレータのソリューションは、Patzold,M.らによる”A Study of Stochastic and Deterministic Procedures for the Design of Simulation Models for Spatial Channels”,１９２４〜１９３１頁、パーソナル、インドア、及び移動無線通信に関するProc.第１３回ＩＥＥＥ国際シンポジウム、ＩＥＥＥ  ＰＩＭＲＣ２００２、ボルトガル共和国リスボン、２００２年９月１５〜１８。このソリューションには、ＤＵＴ４００に対する疑似無線チャネルのシフトを電気的に実現することが欠けている。
【００８４】
上記では、疑似無線チャネルのシフトが２次元的に説明された。しかし一実施形態では、ＯＴＡチャンバ内の平面上に配置されていないアンテナ素子を利用して、疑似無線チャネルのシフトが３次元的に実行される。そして、少なくとも１つのビームを有する角度スペクトルの方向が立体角で決定される。
【００８５】
一実施形態では、ＤＵＴ４００は、少なくとも１つの対応するクラスタからの反射及び／又は散乱に関連した少なくとも１つのビーム内にある、アンテナ素子４０２〜４１６からの信号を受信する可能性がある。あるいは、アンテナ素子４０２〜４１６が、少なくとも１つの対応するクラスタからの反射及び／又は散乱に関連した少なくとも１つのビーム内にある、ＤＵＴ４００からの信号を受信する可能性がある。電力角度スペクトルは、ＤＵＴ４００に対してシフトするが、アンテナ素子４０２〜４１６及びＤＵＴ４００の位置関係はそのまま変わらない。
【００８６】
図３０は、エミュレータの疑似無線チャネルを介して、試験対象電気機器と通信する方法のフロー図を示している。ステップ３０００では、無響チャンバ内において疑似無線チャネルに対して様々な方法を用いて試験対象電気機器と通信するために、エミュレータと接続された複数のアンテナ素子に対して疑似無線チャネルをシフトする。
【００８７】
本実施形態は、例えば、ＡＳＩＣ又はＶＬＳＩ回路（特定用途向け集積回路、超大規模集積回路）で実施される。あるいは又はさらに、方法ステップの実施形態は、エミュレータの疑似無線チャネルを介して、試験対象の電気機器と通信するためのコンピュータ処理を実行させる命令を含むコンピュータプログラムとして実現される。エミュレータは、電気回路及び／又はコンピュータプログラムに基づいて、アンテナ素子の使用、ビームの形成、及び無響チャンバ内の各ビームの少なくとも１つのシフトを制御する。
【００８８】
コンピュータプログラムは、コンピュータ又はプロセッサによって読み取り可能なコンピュータプログラム媒体上に記憶される。コンピュータプログラム媒体は、例えば電気的、磁気的、光学的、赤外線又は半導体の、システム、デバイス又は伝送媒体であるが、これに限定するものではない。コンピュータプログラム媒体は、コンピュータ可読媒体、プログラム記憶媒体、記録媒体、コンピュータ可読メモリ、ランダムアクセスメモリ、消去可能プログラマブル読み取り専用メモリ、コンピュータ可読ソフトウェア配布パッケージ、コンピュータ可読信号、コンピュータ可読遠隔通信信号、コンピュータ可読印刷物、コンピュータ可読圧縮ソフトウェアパッケージ、の少なくとも１つを含んでいる。
【００８９】
本発明の一例を添付図面を参照して説明したが、本発明がこれに限定されるものではな
いことは明らかであり、添付の特許請求の範囲内で幾つかの方法で修正することが可能である。
【図１】
0005358021.tif 000027
【図２】
0005358021.tif 000028
【図３】
0005358021.tif 000029
【図４】
0005358021.tif 000030
【図５】
0005358021.tif 000031
【図６】
0005358021.tif 000032
【図７】
0005358021.tif 000033
【図８】
0005358021.tif 000034
【図９】
0005358021.tif 000035
【図１０】
0005358021.tif 000036
【図１１】
0005358021.tif 000037
【図１２】
0005358021.tif 000038
【図１３】
0005358021.tif 000039
【図１４】
0005358021.tif 000040
【図１５】
0005358021.tif 000041
【図１６】
0005358021.tif 000042
【図１７】
0005358021.tif 000043
【図１８】
0005358021.tif 000044
【図１９】
0005358021.tif 000045
【図２０】
0005358021.tif 000046
【図２１】
0005358021.tif 000047
【図２２】
0005358021.tif 000048
【図２３】
0005358021.tif 000049
【図２４】
0005358021.tif 000050
【図２５】
0005358021.tif 000051
【図２６】
0005358021.tif 000052
【図２７】
0005358021.tif 000053
【図２８】
0005358021.tif 000054
【図２９】
0005358021.tif 000055
【図３０】
0005358021.tif 000056