(57)【特許請求の範囲】
【請求項１】
１つの第１無線局装置と複数の第２無線局装置により構成され、前記第１無線局装置は複数本のアンテナで構成される第１アンテナ群を備え、前記第２無線局装置は複数本のアンテナで構成される第２アンテナ群を備え、前記第１無線局装置の前記第１アンテナ群および前記複数の第２無線局装置の全てまたはその一部の前記第２アンテナ群により構成されるＭＩＭＯ（Multiple Input Multiple Output）チャネルを介して複数の信号系統を同一周波数チャネルおよび同一時刻に空間多重してＭＩＭＯ通信する無線通信システムであって、
前記第１無線局装置は、
前記第１アンテナ群の各アンテナから送信する際に付与する複数の送信ウエイトを成分とする複数種類の送信ウエイトベクトルを記憶する手段と、
該送信ウエイトベクトルを該第１無線局装置および前記第２無線局装置の双方で既知であるパターンのチャネル推定用信号に乗算した信号を生成する手段と、
複数種類の送信ウエイトベクトルに対応した複数の前記信号を前記第１アンテナ群の各アンテナを介して送信する手段と、
その信号を受信した前記第２無線局装置が返送する推定ＳＮＲに関する情報を受信する手段と、
その受信した情報をもとに複数の前記第２無線局装置との間で同時に空間多重して通信することができる当該第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列を決定する手段と、
該決定した、同時に空間多重して通信することができる前記第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列の対応をデータベース化して記録する手段と、
該データベース化された情報から、送信すべきデータが存在する第２無線局装置の間で空間多重ができる該第２無線局装置の組み合わせを検索する手段と、
空間多重ができる第２無線局装置の組み合わせが検索された場合には、その組み合わせに対応した送信ウエイト行列を用い送信する信号を要素とする送信信号ベクトルに該送信ウエイト行列を乗算する手段と、
該乗算結果を前記第1アンテナ群の各アンテナから送信する手段と、を備え、
前記第２無線局装置は、
前記第１無線局装置より既知であるパターンのチャネル推定用信号に複数の送信ウエイトベクトルを乗算した信号を受信する手段と、
該受信信号に対しチャネル推定する手段と、
前記第２のアンテナ群の各アンテナ毎の伝達関数を要素とするベクトルすなわち伝達関数ベクトルに対し、該ベクトルの各成分の絶対値のべき乗値ないしはその近似値を全アンテナに対して総和をとる手段と、
該総和の値が大きい方からひとつまたは複数の伝達関数ベクトルを選択する手段と、
選択された該伝達関数ベクトル毎に、該ベクトル以外の他の伝達関数ベクトルとの内積ないしはその近似値を演算する手段と、
該内積値が所定の閾値以下の伝達関数ベクトルの中から、ひとつまたは複数個を内積値の小さい方から選択する手段と、
その結果を第１無線局装置に対して通知する手段と、を備える
ことを特徴とする無線通信システム。
【請求項２】
１つの第１無線局装置と複数の第２無線局装置により構成され、前記第１無線局装置は複数本のアンテナで構成される第１アンテナ群を備え、前記第２無線局装置は複数本のアンテナで構成される第２アンテナ群を備え、前記第１無線局装置の前記第１アンテナ群および前記複数の第２無線局装置の全てまたはその一部の前記第２アンテナ群により構成されるＭＩＭＯ（Multiple Input Multiple Output）チャネルを介して複数の信号系統を同一周波数チャネルおよび同一時刻に空間多重してＭＩＭＯ通信する無線通信システムにおける無線通信方法であって、
前記第１無線局装置は、
前記第1アンテナ群の各アンテナから送信する際に付与する複数の送信ウエイトを成分とする複数種類の送信ウエイトベクトルを該第１無線局装置および前記第２無線局装置の双方で既知であるパターンのチャネル推定用信号に乗算した信号を生成する処理と、
複数種類の送信ウエイトベクトルに対応した複数の前記信号を前記第１アンテナ群の各アンテナを介して送信する処理と、
その信号を受信した前記第２無線局装置が返送する推定ＳＮＲに関する情報を受信する処理と、
その受信した情報をもとに複数の前記第２無線局装置との間で同時に空間多重して通信することができる当該第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列を決定する処理と、
該決定した、同時に空間多重して通信することができる前記第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列の対応をデータベース化して記録する処理と、
該データベース化された情報から、送信すべきデータが存在する第２無線局装置の間で空間多重ができる該第２無線局装置の組み合わせを検索する処理と、
空間多重ができる第２無線局装置の組み合わせが検索された場合には、その組み合わせに対応した送信ウエイト行列を用い送信する信号を要素とする送信信号ベクトルに該送信ウエイト行列を乗算する処理と、
該乗算結果を前記第1アンテナ群の各アンテナから送信する処理とを実施し、
前記第２無線局装置は、
前記第１無線局装置より既知であるパターンのチャネル推定用信号に複数の送信ウエイトベクトルを乗算した信号を受信する処理と、
該受信信号に対しチャネル推定する処理と、
前記第２のアンテナ群の各アンテナ毎の伝達関数を要素とするベクトルすなわち伝達関数ベクトルに対し、該ベクトルの各成分の絶対値のべき乗値ないしはその近似値を全アンテナに対して総和をとる処理と、
該総和の値が大きい方からひとつまたは複数の伝達関数ベクトルを選択する処理と、
選択された該伝達関数ベクトル毎に、該ベクトル以外の他の伝達関数ベクトルとの内積ないしはその近似値を演算する処理と、
該内積値が所定の閾値以下の伝達関数ベクトルの中から、ひとつまたは複数個を内積値の小さい方から選択する処理と、
その結果を第１無線局装置に対して通知する処理とを実施する
ことを特徴とする無線通信方法。 
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ひとつの基地局装置が複数の端末局装置と同時に空間多重して通信を行うマルチユーザＭＩＭＯ（Multiple Input Multiple Output）通信の手法を利用した無線通信システムおよび無線通信方法に関する。
【背景技術】
【０００２】
近年、２．４ＧＨｚ帯または５ＧＨｚ帯を用いた高速無線アクセスシステムとして、ＩＥＥＥ８０２．１１ｇ規格、ＩＥＥＥ８０２．１１ａ規格などの普及が目覚しい。これらのシステムでは、マルチパスフェージング環境での特性を安定化させるための技術である直交周波数分割多重（ＯＦＤＭ：Orthogonal Frequency Division Multiplexing）変調方式を用い、最大で５４Ｍｂｐｓの伝送速度を実現している。ただし、ここでの伝送速度とは物理レイヤ上での伝送速度であり、実際にはＭＡＣ（Medium Access Control）レイヤでの伝送効率が５０〜７０％程度であるため、実際のスループットの上限値は３０Ｍｂｐｓ程度である。一方で、有線ＬＡＮの世界ではEthernetの１００Base−Tインタフェースをはじめ、各家庭にも光ファイバを用いたＦＴＴＨ（Fiber to the home）が普及し、また１００Ｍｂｐｓの高速回線の提供が普及しており、無線ＬＡＮの世界においても更なる伝送速度の高速化が求められている。そのための技術としては、ＭＩＭＯ技術が有力である。このＭＩＭＯ技術とは、送信局装置側において複数の送信アンテナから同一チャネル上で異なる独立な信号を空間多重して送信し、受信局装置側において同じく複数のアンテナを用いて信号を受信し、各送信アンテナ／受信アンテナ間の伝達関数行列を求め、この行列を用いて送信局装置側の各アンテナから送信した独立な信号を受信側にて推定し、データを再生するものである。
【０００３】
このＭＩＭＯ通信の究極の形態として、ひとつの基地局装置が複数の端末局装置と同時に空間多重して通信を行うマルチユーザＭＩＭＯ通信がある。従来のＭＩＭＯ通信技術は、基地局装置と通信する端末局装置は時分割で棲み分けられ、常に１対１の通信であった。通常、これをシングルユーザＭＩＭＯ通信と呼ぶ。このシングルユーザＭＩＭＯ通信では、端末局装置は比較的少数のアンテナを小さな筐体の中に備え、各アンテナ間の相関は非常に強い状況である。この様な状態では、時空間資源を十分に有効活用することはできない。一般には、アンテナ間の距離を離し、アンテナ間の相関を低くすることで高次の空間多重を実現し、且つ、各サブチャネル（ＭＩＭＯ伝送において各信号系列の伝送に用いられる仮想チャネル）のＳＮＲを高く保つことが可能である。
【０００４】
マルチユーザＭＩＭＯ通信は、ひとつひとつの端末においては空間多重するＭＩＭＯチャネル数を少なくする一方、複数の異なる端末と同時に同一周波数チャネルで空間多重を行う。図１に示すように、マルチユーザＭＩＭＯシステムが、基地局装置１０１、端末局装置＃１〜＃３（１０２〜１０４）を備えているとする。このとき、実際にひとつの基地局装置が収容する端末局装置数は多数であるが、その中の数局装置を選び出し（端末局装置＃１〜＃３:１０２〜１０４）通信を行う。各端末局装置は基地局装置と比較して送受信アンテナ数が一般的に少ない。例えば基地局装置から端末局装置方向への通信（ダウンリンク）を行う場合を考える。基地局装置１０１は多数のアンテナを用いて、複数の指向性ビームを形成する。例えば、各端末局装置１０２〜１０４に対してそれぞれ３つのＭＩＭＯチャネルを割り当て、全体としては９系統の信号系列を送信する場合を考える。その際、端末局装置１０２に対して送信する信号は、端末局装置１０３および端末局装置１０４方向には指向性利得が極端に低くなるように調整し、この結果として端末局装置１０３および端末局装置１０４への干渉を抑制する。同様に、端末局装置１０３に対して送信する信号は、端末局装置１０２および端末局装置１０４方向には指向性利得が極端に低くなるように調整する。同様の処理を端末局装置１０４にも施す。この様に指向性制御を行う理由は、例えば端末局装置１０２においては、端末局装置１０３および端末局装置１０４で受信した信号の情報を知る術がないため、端末間での協調的な受信処理ができない。つまり、アンテナが３本しかない端末局装置１０２のみの受信処理において、９系統の全ての信号系列を信号分離することは非常に厳しい。そこで、各端末局装置には他の端末局装置宛の信号が受信されないように、送信側で干渉分離を事前に行う。
【０００５】
次に、指向性ビームの形成方法について、以下に説明を加える。例えば、基地局装置１０１の第ｊアンテナと端末局装置１０２の第１受信アンテナとの間の伝達関数をｈ１ｊと表記することにする。基地局装置１０１のｊ＝１〜９の全てのアンテナに関する伝達関数を用い、行ベクトルｈ１を（ｈ１１，ｈ１２，ｈ１３，…，ｈ１８，ｈ１９）と表記する。同様に基地局装置１０１と端末局装置１０２の第２受信アンテナ、第３受信アンテナの伝達関数をｈ２ｊおよびｈ３ｊとし、対応する行ベクトルｈ２およびｈ３を（ｈ２１，ｈ２２，ｈ２３，…，ｈ２８，ｈ２９）、（ｈ３１，ｈ３２，ｈ３３，…，ｈ３８，ｈ３９）とする。端末局装置１０３、端末局装置１０４の受信アンテナも同様に、行ベクトルｈ４〜ｈ９を（ｈ４１，ｈ４２，ｈ４３，…，ｈ４８，ｈ４９）〜（ｈ９１，ｈ９２，ｈ９３，…，ｈ９８，ｈ９９）とする。加えて、基地局装置１０１が送信する９系統の信号をｔ１〜ｔ９と表記し、これを成分とする列ベクトルをＴｘ［ａｌｌ］＝(ｔ１，ｔ２，ｔ３，…，ｔ８，ｔ９)Ｔと表記する。ここで、右肩のＴの文字はベクトル、行列の転置を表す。また同様に、端末局装置１０２〜１０４の９本のアンテナでの受信信号をｒ１〜ｒ９と表記し、これを成分とする列ベクトルをＲｘ［ａｌｌ］＝(ｒ１，ｒ２，ｒ３，…，ｒ８，ｒ９)Ｔと表記する。最後に、行ベクトルｈ１〜ｈ９を第１から第９行成分とする行列を、全体伝達関数行列Ｈ［ａｌｌ］と表記する。この様に表記した場合、システム全体としては以下の関係式が成り立つ。
【０００６】
【数１】
0004708205.tif 000002
【０００７】
ここで、受信側において効率的に信号分離を行えるように、送信局装置側では送信指向性制御を行う。ここでは９行９列の送信ウエイト行列Ｗを導入し、（式１）を以下のように書き換える。
【０００８】
【数２】
0004708205.tif 000003
【０００９】
さらに、送信ウエイト行列Ｗを列ベクトルｗ１〜ｗ９に分解し、Ｗ＝（ｗ１，ｗ２，ｗ３，…，ｗ８，ｗ９）と表記すると、以下のように表すことができる。
【００１０】
【数３】
0004708205.tif 000004
【００１１】
ここで、例えば６つの行ベクトルｈ４〜ｈ９と３つの列ベクトルｗ１〜ｗ３の乗算（各成分同士を乗算したものの総和、複素ベクトルの場合は内積とは異なる）が全てゼロになるようにｗ１〜ｗ３を選ぶことを考える。同様に、行ベクトルｈ１〜ｈ３およびｈ７〜ｈ９と列ベクトルｗ４〜ｗ６の積、行ベクトルｈ１〜ｈ６と列ベクトルｗ７〜ｗ９の積の全てがゼロになるように選ぶことにする。すると、（式３）に示す９行９列の行列は、３行３列の９個の部分行列を用いて表記すると以下のように表すことができる。
【００１２】
【数４】
0004708205.tif 000005
【００１３】
ここで、部分行列Ｈ［１］、Ｈ［２］、Ｈ［３］は３行３列の行列であり、またこの数式において「０」と示した部分行列は成分が全てゼロの３行３列の行列である。この様な条件を満たす変換行列Ｗを選択することで、（式４）は以下の３つの関係式に分解できる。
【００１４】
【数５】
0004708205.tif 000006
【００１５】
【数６】
0004708205.tif 000007
【００１６】
【数７】
0004708205.tif 000008
【００１７】
ここで、Ｔｘ［１］＝(ｔ１，ｔ２，ｔ３)Ｔ、Ｔｘ［２］＝(ｔ４，ｔ５，ｔ６)Ｔ、Ｔｘ［３］＝(ｔ７，ｔ８，ｔ９)Ｔ、Ｒｘ［１］＝(ｒ１，ｒ２，ｒ３)Ｔ、Ｒｘ［２］＝(ｒ４，ｒ５，ｒ６)Ｔ、Ｒｘ［３］＝(ｒ７，ｒ８，ｒ９)Ｔとした。この様にして、３つのシングルユーザＭＩＭＯ通信とみなすことができるようになる。以上がマルチユーザＭＩＭＯシステムの概要である。
【００１８】
マルチユーザＭＩＭＯ通信においては、（式３）に示した送信ウエイト行列の算出が必須である。一般の端末局装置の場合には、端末局装置の場所移動や周囲の動きにより、伝搬チャネルは時間的に変動する。シングルユーザＭＩＭＯ通信の場合にも固有モード伝送のために送信指向性制御を行うことがあるが、１台の送信局装置と１台の受信局装置間のＭＩＭＯチャネルのみを考慮して最適化を行えばよかった。例えばシングルユーザＭＩＭＯ通信では、送信ウエイト行列としてはＭＩＭＯチャネルに関する伝達関数行列を対角化するユニタリー行列などが用いられるが、その伝達関数行列のサイズ（行列の行数、列数）は比較的小さく、行列演算も比較的演算量が少なくて済んだ。
【００１９】
しかし、マルチユーザＭＩＭＯ通信においては、同時に空間多重される端末数の増加に伴い行列のサイズが増大し、行列の対角化をはじめとする各種演算処理の計算量が指数関数的に増大する。一方で、ある端末局装置が静止している場合でも、他の端末局装置に動きがあると、再度、全体での送信ウエイト行列の最適化を行う必要がある。これは、ある端末局装置に対する信号送信が他の端末局装置に対して干渉とならないような指向性制御を行う場合、ひとつの端末局装置に動きがあると、他の端末局装置の送信指向性を全て調整しなおす必要があるからである。
【００２０】
また、多数の端末局装置がひとつの基地局装置の配下にある場合、同時に空間多重できる端末はそれらの全体の端末の中の一部であり、少数である。通信時には、多数の端末から同時に空間多重するのに適した端末の組み合わせを選択する必要があるが、ある端末局装置への送信信号が他の端末局装置に対してどの程度干渉となるかがシステム全体の特性を左右するため、同時に空間多重する端末局装置の組み合わせの選び方が重要である。しかし、その組み合わせの最適解はチャネルの時変動と共に変わるため、逐次、組み合わせの最適化処理を行わなければならない。例えば全体でN局装置の端末局装置があり、その中からＭ局の装置を選択して空間多重をする場合、その組み合わせはＮＣＭ通りとなる。この膨大な組み合わせに対し、それぞれ送信指向性制御の送信ウエイト行列を求め、初めて特性の優劣が評価できるのであるが、前述の様に送信ウエイト行列の演算は容易ではない。したがって、この様に逐次環境が激しく変動する環境では、マルチユーザＭＩＭＯにおいて求められる送信指向性制御の負荷は非常に大きくなる。この様に、マルチユーザＭＩＭＯ通信の運用中における処理・演算量の削減は必須の課題となっている。このための対策として、非特許文献１に示す方式が提案されている。
【非特許文献１】Jaehak Chung, Chan-Soo Hwang, Kiho Kim and Young Kyun Kim“A Random Beamforming Technique in ＭＩＭＯ Systems Exploiting Multiuser Diversity”IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS, VOL. 21, NO. 5, JUNE 2003
【発明の開示】
【発明が解決しようとする課題】
【００２１】
ここで、上述の従来技術において用いられるランダムビームフォーミングの原理は以下の通りである。まず、基地局装置は自局装置のもつアンテナで所定の数の信号系列を空間多重して送信する際の送信ウエイト行列を何らかの手段で定める。例えば、ランダムに選んだ直交ベクトルであったり、ランダムに仮定した伝達関数行列に対して特異値分解を行って得たユニタリー行列でも構わない。この送信ウエイト行列を用いて、（式２）と同様に既知のパターンのチャネル推定用信号を送信する。各端末局装置は、受信した信号をもとにチャネル推定を行い、基地局装置との間の伝達関数行列を取得する。この伝達関数行列をもとに受信指向性形成を行った場合、空間多重された各信号系列を受信した場合の推定ＳＮＲ（Signal to Noise Ratio：信号対雑音比）を求め、これを基地局装置に対して返送する。基地局装置はこの返送された信号系列毎の推定ＳＮＲをもとに、各端末が所望のＳＮＲで信号を受信できるように、空間多重する端末の組み合わせを選択する。この組み合わせに対し、空間多重数や伝送モード（変調方式と誤り訂正の符号化率の組み合わせ）の最適化を行う。
【００２２】
具体的には、例えば基地局装置が４つの信号系列を空間多重して送信するケースを例に取って説明する。ランダムな送信ウエイト行列を用いた際に端末局装置１０２においては信号系列＃１〜信号系列＃２の推定ＳＮＲは高いが信号系列＃３〜信号系列＃４の推定ＳＮＲは非常に低く、逆に端末局装置１０３においては信号系列＃３〜信号系列＃４の推定ＳＮＲは高いが信号系列＃１〜信号系列＃２の推定ＳＮＲは非常に低いとする。この場合、端末局装置１０２宛ての信号を信号系列＃１〜信号系列＃２に割り当て、端末局装置１０３宛ての信号を信号系列＃３〜信号系列＃４に割り当てることで同時に信号を空間多重して送信しても、事前に端末局装置が学習した受信指向性形成のための受信ウエイト行列を用いて信号受信すれば、他局装置宛の信号は抑制できる。実際の基地局装置の判断は、単なる推定ＳＮＲではなく、干渉信号の推定ＳＮＲを考慮した信号対干渉雑音比ＳＮＩＲ（Signal to Noise Interference Ratio）をもとに推定品質を判断し、決定する。
【００２３】
図７は従来技術における基地局装置の送信処理フローを示す図である。
基地局装置が信号送信にあたり、ランダムな送信ウエイト行列を適宜生成し（Ｓ１０１）、複数の既知のプリアンブル信号を要素とする送信信号ベクトルに送信ウエイト行列を乗算し（Ｓ１０２）、空間多重しながらこれらの信号を送信する（Ｓ１０３）。その後、端末局装置からの推定ＳＮＲ情報が収容された返信信号を受信し（Ｓ１０４）、これらの報告された推定ＳＮＲ情報から空間多重可能な端末の組み合わせを検索する（Ｓ１０５）。この結果として空間多重が可能と判断された端末局装置に対し、先の送信ウエイト行列を用いて空間多重しながら信号を送信する（Ｓ１０６）。
【００２４】
図８は従来技術における端末局装置の受信処理フロー（事前処理）を示す図である。
事前処理フローとは、実際の信号を受信する前の処理であり、図７におけるステップＳ１０３とステップＳ１０４の間に実施される処理に相当する。図８より、まず、端末局装置がランダムな送信ウエイト行列を乗算したプリアンブル信号を受信すると（Ｓ１１１）、チャネル推定により伝達関数行列を取得する（Ｓ１１２）。その後、伝達関数行列の特異値分解を実施する（Ｓ１１３）。この特異値分解によって得られた固有値及び固有ベクトルを用いて推定ＳＮＲを算出する（Ｓ１１４）。その後、この推定ＳＮＲ情報を基地局装置に対して通知する（Ｓ１１５）。
【００２５】
図９は従来技術における端末局装置の受信処理フロー（データ受信処理）を示す図である。
基地局装置では、図７のステップＳ１０６によりデータを送信する。端末局装置は、先の図８の事前処理に続けて、基地局装置から送信されるデータの受信処理を図９のように行う。まず信号を受信すると（Ｓ１２１）、図８の事前処理のステップＳ１１３にて求めた伝達関数行列を対角化するユニタリー行列を受信信号ベクトルに乗算し（Ｓ１２２）、他局装置宛の信号の抑圧を行う。その後、自局装置宛の信号のチャネル推定を行い（Ｓ１２３）、ここで求めた伝達関数行列を用いて後続するデータの信号検出処理を実施し（Ｓ１２４）、信号検出後の信号を合成してデータを再生し、出力する（Ｓ１２５）。以上が従来技術であるランダムビームフォーミングの概要である。
【００２６】
しかし、このランダムビームフォーミング法では、ランダムに行ったひとつの送信指向性制御に適した端末局装置の組み合わせが存在する場合に利用可能であり、且つ、それらの端末宛の送信データが送信バッファ内に蓄積されていて初めて意味がある。しかし、実際には送信データのある限定された端末局装置の中で、たった一つのランダムな送信ウエイト行列に適した端末局装置が存在する可能性は低い。つまり、このランダムビームフォーミング法は膨大な数の端末局装置が存在し、且つそれぞれに対し送信すべきデータが定常的に存在する場合に有効であるが、現実にはその様な条件は存在しない。したがって、新たなマルチユーザＭＩＭＯ通信における処理・演算量の削減方法の開発が課題となっている。
【００２７】
そこでこの発明は、マルチユーザＭＩＭＯ通信の運用中における処理・演算量の削減することができる無線通信システムおよび無線通信方法を提供することを目的としている。
【課題を解決するための手段】
【００２８】
上記目的を達成するために、本発明は、１つの第１無線局装置と複数の第２無線局装置により構成され、前記第１無線局装置は複数本のアンテナで構成される第１アンテナ群を備え、前記第２無線局装置は複数本のアンテナで構成される第２アンテナ群を備え、前記第１無線局装置の前記第１アンテナ群および前記複数の第２無線局装置の全てまたはその一部の前記第２アンテナ群により構成されるＭＩＭＯ（Multiple Input Multiple Output）チャネルを介して複数の信号系統を同一周波数チャネルおよび同一時刻に空間多重してＭＩＭＯ通信する無線通信システムであって、前記第１無線局装置は、前記第１アンテナ群の各アンテナから送信する際に付与する複数の送信ウエイトを成分とする複数種類の送信ウエイトベクトルを記憶する手段と、該送信ウエイトベクトルを該第１無線局装置および前記第２無線局装置の双方で既知であるパターンのチャネル推定用信号に乗算した信号を生成する手段と、複数種類の送信ウエイトベクトルに対応した複数の前記信号を前記第１アンテナ群の各アンテナを介して送信する手段と、その信号を受信した前記第２無線局装置が返送する推定ＳＮＲに関する情報を受信する手段と、その受信した情報をもとに複数の前記第２無線局装置との間で同時に空間多重して通信することができる当該第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列を決定する手段と、該決定した、同時に空間多重して通信することができる前記第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列の対応をデータベース化して記録する手段と、該データベース化された情報から、送信すべきデータが存在する第２無線局装置の間で空間多重ができる該第２無線局装置の組み合わせを検索する手段と、空間多重ができる第２無線局装置の組み合わせが検索された場合には、その組み合わせに対応した送信ウエイト行列を用い送信する信号を要素とする送信信号ベクトルに該送信ウエイト行列を乗算する手段と、該乗算結果を前記第1アンテナ群の各アンテナから送信する手段と、を備え、前記第２無線局装置は、前記第１無線局装置より既知であるパターンのチャネル推定用信号に複数の送信ウエイトベクトルを乗算した信号を受信する手段と、該受信信号に対しチャネル推定する手段と、前記第２のアンテナ群の各アンテナ毎の伝達関数を要素とするベクトルすなわち伝達関数ベクトルに対し、該ベクトルの各成分の絶対値のべき乗値ないしはその近似値を全アンテナに対して総和をとる手段と、該総和の値が大きい方からひとつまたは複数の伝達関数ベクトルを選択する手段と、選択された該伝達関数ベクトル毎に、該ベクトル以外の他の伝達関数ベクトルとの内積ないしはその近似値を演算する手段と、該内積値が所定の閾値以下の伝達関数ベクトルの中から、ひとつまたは複数個を内積値の小さい方から選択する手段と、その結果を第１無線局装置に対して通知する手段と、を備えることを特徴とする無線通信システムである。
【００２９】
また本発明は、１つの第１無線局装置と複数の第２無線局装置により構成され、前記第１無線局装置は複数本のアンテナで構成される第１アンテナ群を備え、前記第２無線局装置は複数本のアンテナで構成される第２アンテナ群を備え、前記第１無線局装置の前記第１アンテナ群および前記複数の第２無線局装置の全てまたはその一部の前記第２アンテナ群により構成されるＭＩＭＯ（Multiple Input Multiple Output）チャネルを介して複数の信号系統を同一周波数チャネルおよび同一時刻に空間多重してＭＩＭＯ通信する無線通信システムにおける無線通信方法であって、前記第１無線局装置は、前記第1アンテナ群の各アンテナから送信する際に付与する複数の送信ウエイトを成分とする複数種類の送信ウエイトベクトルを該第１無線局装置および前記第２無線局装置の双方で既知であるパターンのチャネル推定用信号に乗算した信号を生成する処理と、複数種類の送信ウエイトベクトルに対応した複数の前記信号を前記第１アンテナ群の各アンテナを介して送信する処理と、その信号を受信した前記第２無線局装置が返送する推定ＳＮＲに関する情報を受信する処理と、その受信した情報をもとに複数の前記第２無線局装置との間で同時に空間多重して通信することができる当該第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列を決定する処理と、該決定した、同時に空間多重して通信することができる前記第２無線局装置の組み合わせおよびその際に用いるべき送信ウエイト行列の対応をデータベース化して記録する処理と、該データベース化された情報から、送信すべきデータが存在する第２無線局装置の間で空間多重ができる該第２無線局装置の組み合わせを検索する処理と、空間多重ができる第２無線局装置の組み合わせが検索された場合には、その組み合わせに対応した送信ウエイト行列を用い送信する信号を要素とする送信信号ベクトルに該送信ウエイト行列を乗算する処理と、該乗算結果を前記第1アンテナ群の各アンテナから送信する処理とを実施し、前記第２無線局装置は、前記第１無線局装置より既知であるパターンのチャネル推定用信号に複数の送信ウエイトベクトルを乗算した信号を受信する処理と、該受信信号に対しチャネル推定する処理と、前記第２のアンテナ群の各アンテナ毎の伝達関数を要素とするベクトルすなわち伝達関数ベクトルに対し、該ベクトルの各成分の絶対値のべき乗値ないしはその近似値を全アンテナに対して総和をとる処理と、該総和の値が大きい方からひとつまたは複数の伝達関数ベクトルを選択する処理と、選択された該伝達関数ベクトル毎に、該ベクトル以外の他の伝達関数ベクトルとの内積ないしはその近似値を演算する処理と、該内積値が所定の閾値以下の伝達関数ベクトルの中から、ひとつまたは複数個を内積値の小さい方から選択する処理と、その結果を第１無線局装置に対して通知する処理とを実施することを特徴とする無線通信方法である。
【発明の効果】
【００３０】
本発明によれば、マルチユーザＭＩＭＯ通信を行う際の送信指向性制御を行うための送信ウエイト行列を、高次の行列の特異値分解や直交化などの膨大な線形演算を行うことなしに簡易に求めることが可能となる。また、この演算の簡易化技術として提案されていた従来のランダムビームフォーミング技術では、送信ウエイトベクトルの組み合わせである送信ウエイト行列を一方的に基地局装置が決定し、この送信ウエイト行列に適合した端末局装置を検索してマルチユーザＭＩＭＯ通信を行なうとしていたために、偶然にもその条件に一致している端末局装置が存在し、且つそれらの端末局装置宛てのデータが送信バッファ内に存在する場合にのみ有効であった。しかし、本発明では、ランダムに生成する送信ウエイトベクトルを、端末局装置から返信される情報に応じて自由に組み合わせて利用することが可能となるため、本技術が有効に適用可能な状況の頻度が飛躍的に増大する。その結果として、マルチユーザＭＩＭＯ通信を安定して有効に利用することが可能となり、周波数資源の効率的な運用が可能となる。
【発明を実施するための最良の形態】
【００３１】
以下、本発明の一実施形態による無線通信システムを図面を参照して説明する。
図１はマルチユーザＭＩＭＯによる無線通信システムの構成を示す図である。
この図において、１０１は基地局装置、１０２〜１０４は端末局装置＃１〜＃３を示す。なお無線通信システムにおいて、実際は１つの基地局装置が収容する端末局装置数は多数であるが、その中の数局装置を選び出し（図では端末局装置＃１〜＃３:１０２〜１０４）通信を行う。各端末局装置は基地局装置と比較して送受信アンテナ数が一般的に少ない。
【００３２】
まず基地局装置１０１は、多数の異なる送信ウエイトベクトルをデータテーブル内に記憶している。また基地局装置１０１は、チャネル推定用の所定のプリアンブル信号に、これらの送信ウエイトベクトルを乗算した変換プリアンブル信号も合わせてデータテーブル内に記憶している。例えば、ＯＦＤＭ変調方式を用いたシステムであれば、１シンボルあたりの信号パターンを記録しておく。これらの多数の変換プリアンブル信号は連続的に複数パターンを一括して送信する。例えば、１００個の送信ウエイトベクトルを用意している場合には、１００ＯＦＤＭシンボルに渡り一括して信号送信を行う。
【００３３】
各端末局装置１０２〜１０４ではこれらの変換プリアンブル信号を全て受信する。変換前のプリアンブル信号のパターンは端末局装置１０２〜１０４においても既知であり、既知のパターンを用いて受信アンテナ毎に相関をとり、伝達関数情報を取得する。端末局装置１０２〜１０４は、この様にして求めた伝達関数情報を各アンテナに対応した成分とする伝達関数ベクトルを、各送信ウエイトベクトルに対する情報として記憶する。同一の送信ウエイトベクトルに対する変換プリアンブル信号の送信は、通信の間に適宜繰り返し送信される。その都度、伝達関数ベクトルを求めるが、この伝達関数ベクトルは１回の信号受信毎に更新しても良いし、複数回の受信に対して平均化した値としてもよい。また均等な平均化ではなく、忘却係数μ（０＜μ＜１）を導入し、過去のベクトルＶｏｌｄ、新規に取得したベクトルＶｎｅｗに対し、更新後の伝達関数ベクトルをμＶｎｅｗ＋（１−μ）Ｖｏｌｄの様に与えても構わない。なお、基地局装置１０１が用いる送信ウエイトベクトルは、送信する順番やその他の制御情報等により、基地局装置１０１と端末局装置間１０２〜１０４で対応をとることが可能なようにしておく。なお、以上の説明では異なる送信ウエイトベクトルを乗算した変換プリアンブル信号は別々のタイミングで送信するとしていたが、お互いが直交関係にあるパターンであるならば、同時に送信を行っても個別に伝達関数ベクトルを取得可能である。
【００３４】
各端末局装置１０２〜１０４は、この様にして求めた複数の伝達関数ベクトルに対し、それぞれの各成分の絶対値ないしは絶対値のべき乗（例えば２乗）の全成分の総和（以降、伝達関数ベクトル絶対値と呼ぶ）を求め、この値が大きくなる方から所定の数だけ伝達関数ベクトルを求める。ここで、この伝達関数ベクトル絶対値は、各伝達関数ベクトルに対する推定ＳＮＲに比例する物理量である。また各端末局装置１０２〜１０４は、この様にして求めた伝達関数ベクトルに対し、他の伝達関数ベクトルとのベクトル積（内積）をとり、この値が小さくなる伝達関数ベクトルを検索する。ここでは、所定の数だけ小さい方から伝達関数ベクトルを検索する他、所定の閾値以下となる伝達関数ベクトルを検索することとしても良い。なお、閾値は単一である必要は無く、複数段階の閾値を設けても良い。
【００３５】
以上処理の結果を各端末局装置１０２〜１０４は基地局装置１０１に通知する。通知する内容は、単に伝達関数ベクトル絶対値が大きい伝達関数ベクトルの識別番号（例えば、その信号を受信した順番等）と、その伝達関数ベクトルに対して内積値が小さい伝達関数ベクトルの識別番号の組であっても良い。さらには、伝達関数ベクトル絶対値が大きい伝達関数ベクトルの識別番号に加えてその伝達関数ベクトル絶対値を通知したり、各伝達関数ベクトルに対して内積値が小さい伝達関数ベクトルの識別番号に加えてその内積値を通知しても良い。その他、伝達関数ベクトル絶対値を内積値で除算した値（またはその対数値等）であっても良い。以上のそれぞれの値は、値そのものを通知する他、ある程度の刻み幅で階級化しておき、その値が所属する階級値を報告に用いることも出来る。この様に情報の加工を端末局装置１０２〜１０４側で行い、この結果を基地局装置に集約することで基地局装置１０１の処理量を軽減することが出来る。
【００３６】
基地局装置１０１においては、この結果をもとに効率的な空間多重の組み合わせを求める。例えば、送信ウエイトベクトルを１００種類設定したとする。端末局装置１０１において、１００種類の中で１および２番の伝達関数ベクトルの絶対値が大きいとする。１番の伝達関数ベクトルに対し内積値が小さいベクトルが例えば２０〜３０であったとする。また同様に、２番の伝達関数ベクトルに対し内積値が小さいベクトルが例えば２５〜３５、また３番の伝達関数ベクトルに対し内積値が小さいベクトルが例えば４０〜５０、４番の伝達関数ベクトルに対し内積値が小さいベクトルが例えば６０〜７０であったとする。この場合、例えば基地局装置１０１がある端末局装置に対して２５番から３０番のいずれかの送信ウエイトベクトルを用いて信号送信したとしても、それに重畳する形で１番と２番の送信ウエイトベクトルを用いて端末局装置１０２に対してふたつの信号系列を空間多重することが可能である。
【００３７】
例えば、端末局装置１０３において、２６番と２７番の伝達関数ベクトルの伝達関数ベクトル絶対値が大きく、２６番と２７番の伝達関数ベクトルが１番および２番の伝達関数ベクトルとの内積値が共に小さな値であるならば、１番、２番、２６番、２７番の送信ウエイトベクトルからなる送信ウエイト行列を用いて4つの信号系列を空間多重し、端末局装置１０２および端末局装置１０３に対して同時にマルチユーザＭＩＭＯ通信を行うことが可能である。補足であるが、例えば端末局装置１０４において、伝達関数ベクトルの４５番及び６５番の絶対値が大きな値であるとする。しかし、伝達関数ベクトルの４５番ないしは６５番と１番および２番との内積が大きければ端末基地局装置１０４は同時に空間多重するのには適さない。ただし、同一の端末局装置間においては、空間多重する際の伝達関数ベクトルの内積が必ずしも小さくなる必要はない。また、実際の通信時においてはチャネルの時変動により先の送信ウエイトベクトルに対して補正が必要な場合がある。しかし、ＭＩＭＯ通信固有のパスまでは大きく変動することは少ないので、端末局装置の良好な組み合わせパターンまでは変動は小さいものと考えられる。さらには、基地局装置１０１において、基地局装置１０１の送信アンテナと端末局装置の受信アンテナの間の伝達関数行列に関する情報が取得可能な場合、送信ウエイト行列の最適解ないしは準最適解は例えば「Quentin H. Spencer, A. Lee Swindlehurst and Martin Haardt, “Zero-Forcing Methods for Downlink Spatial Multiplexing in Multiuser ＭＩＭＯ Channels”, IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 52, NO. 2, FEBRUARY 2004」に示す様な方法を用いることで求めることが出来る。
【００３８】
これらの方法は、初期値を出発点として、最小二乗誤差を小さくするような繰り返し処理（イタラティブアプローチ）を用いるのであるが、その初期値として先に求めた送信ウエイトベクトルを利用すると収束速度が速くなり、演算量を減らすことが可能である。
【００３９】
以下、本発明の一実施形態による各種処理フローを順を追って説明する。
図２は本発明の一実施形態による基地局装置の送信信号記憶フローを示す図である。
図２の処理フローにおいては、実際に基地局装置１０１内部で実施しなくても、全く同様の情報を別途作成し、基地局装置１０１内のメモリに保存することでも代用が出来るが、従来方式との対比の関係であえて示しておく。まず基地局装置１０１は、適宜、ランダムな送信ウエイトベクトルを生成する（ステップＳ１）。ここでの送信ウエイトベクトルは、各成分を全くランダムに与えるものであっても構わないし、その結果生成される指向性がランダムになるように生成されたものであっても構わない。後者の場合には、基地局装置１０１のアンテナ配列に対して物理的な仮想指向性を仮定し、この指向性を形成するウエイトをアレーアンテナ技術により定めることになる。この送信ウエイトベクトルをチャネル推定用の既知のプリアンブル信号に乗算し、各アンテナから送信すべき信号を生成する（ステップＳ２）。この生成した信号をメモリに保存する（ステップＳ３）。以上が事前処理であるが、基地局装置１０１の製造時に別途作成した信号をメモリに記録しておけばそれでもよい。
【００４０】
図３は本発明の一実施形態による基地局装置の事前処理フローを示す図である。
次に基地局装置１０１は定期的に、または特に送受信するデータが無く無線帯域に空きがある場合、ないしは新しい端末局装置を新規に収容する際に、メモリ内に記録されている送信ウエイトベクトルが乗算された信号を連続的に送信する（ステップＳ１１）。例えば、１００種類の送信ウエイトベクトルのパターンがあるのであれば、１００シンボル以上をかけて送信する。もちろん、複数回に分割してこれらを送信しても構わない。この様に送信した信号をもとに、端末局装置が各自でチャネル推定を行い推定ＳＮＲを推定した結果を算出し、この情報を基地局装置１０１に返送した信号を受信する（ステップＳ１２）。この報告された情報をもとに、各端末局装置間の希望送信ウエイトベクトルと、その際にペアにしても干渉とならない送信ウエイトベクトルの組み合わせとして、複数端末局装置間で全てが両立する組み合わせを検索する（ステップＳ１３）。選択された組み合わせに対し、その送信ウエイトベクトルにより構成される（各行ベクトルが送信ウエイトベクトルである）送信ウエイト行列を生成する（ステップＳ１４）。これらの組み合わせは一つの端末局装置につき一つの組み合わせだけではなく、複数の組み合わせに所属するようにデータベース化して記憶する（ステップＳ１５）。実際のデータの送信にあたっては、このデータベース化された情報を参照し、空間多重する端末局装置の組み合わせと送信ウエイト行列を決定する。
【００４１】
図４は本発明の一実施形態による基地局装置の送信処理フローを示す図である。
実際にデータを送信する際には、基地局装置１０１は送信すべきデータが存在する端末局装置の中で、空間多重することが可能な組み合わせがあるか否かを検索する（ステップＳ２１）。そして基地局装置１０１は、組み合わせが存在するか否かを条件判断し（ステップＳ２２）、組み合わせがある場合にはその組み合わせ及びその組み合わせに対応した送信ウエイト行列を用い、送信信号ベクトルにその送信ウエイト行列を乗算して信号を送信し（ステップＳ２３）、マルチユーザＭＩＭＯ通信を実現する。また基地局装置１０１はステップＳ２２の判断において組み合わせがない場合には、単独の端末局装置との間でのシングルユーザＭＩＭＯ通信となるが、その端末局装置にとって推定ＳＮＲが高い送信ウエイトベクトルを用いて送信信号に乗算し、信号を送信する（ステップＳ２４）。以上が基地局装置側における処理である。端末局装置側においては、図３で示した基地局装置１０１のステップＳ１１とステップＳ１２との間に行われる事前処理と、図４で示した基地局装置１０１のステップＳ２３（またはステップＳ２４）の処理の後に対応するデータ受信処理を実施する。
【００４２】
図５は本発明の一実施形態による端末局装置の受信処理フロー（事前処理）を示す図である。
端末局装置１０２〜１０４は、基地局装置１０１の上記ステップＳ１１の処理により信号が送信されると、そのランダムな送信ウエイトベクトルが乗算された多数のプリアンブル信号を受信し（ステップＳ３１）、それぞれの信号毎に個別にチャネル推定により伝達関数ベクトルを取得する（ステップＳ３２）。そして伝達関数ベクトルから推定ＳＮＲを算出する。ここでの推定ＳＮＲとは、例えば伝達関数ベクトルの各成分の絶対値のべき乗値の総和、ないしはその近似値によって与えられる。受信した全ての送信ウエイトベクトル毎の推定ＳＮＲ値を求め（ステップＳ３３）、この中から推定ＳＮＲ値が大きい伝達関数ベクトルを検索する（ステップＳ３４）。ここでは、所定の数だけ推定ＳＮＲ値が大きい方から選択したり、閾値以上の伝達関数ベクトルを選択したりする。
【００４３】
この様にして選択された伝達関数ベクトルに対し、他の伝達関数ベクトルとの内積を求め（ステップＳ３５）、この内積値が小さくなる伝達関数ベクトルの組み合わせを検索する（ステップＳ３６）。ここでも、所定の数だけ内積値が小さい方から選択したり、閾値以下の伝達関数ベクトルを選択したりする。この様な組み合わせ結果の情報を、基地局に対して通知する（ステップＳ３７）。ステップＳ３７で通知する内容の詳細は、推定ＳＮＲ値や内積値の全ての値を含んでいたり、これらの値を階級化し、どの階級に所属するかの値であっても構わない。また、具体的な推定ＳＮＲ値や内積値の情報を省略し、単なる組み合わせ情報として、例えばステップＳ３４で選択された伝達関数ベクトルと、ステップＳ３６で選択された伝達関数ベクトルとの内積値が小さくなると判断された伝達関数ベクトルとの組み合わせ情報（実際には、各伝達関数ベクトルの識別番号の組み合わせ情報等）だけであっても構わない。また、それらの中間として、組み合わせ情報に付加する形で、具体的な推定ＳＮＲ値や内積値の全体の情報の一部のみを収容しても構わない。
【００４４】
なお、報告される組み合わせにおいて、良好な推定ＳＮＲを示す伝達関数ベクトルに対しては内積値が大きくなる一方、この伝達関数ベクトルと組合される（内積値が小さい）伝達関数ベクトルとの内積は小さいような受信ウエイトベクトルを算出しておき、これをデータの受信処理の中で利用することも可能である。
【００４５】
図６は、本発明の一実施形態による端末局装置の受信処理フロー（データ受信処理）を示す図である。
これは、図５で示した事前処理の結果を受けて、基地局装置１０１が図３及び図４で示した処理によりユーザデータを含む信号を送信した際の端末局装置１０２〜１０４側の受信処理である。
【００４６】
まず端末局装置１０２〜１０４はそれぞれ信号を受信すると（ステップＳ４１）、各信号系列のＭＩＭＯチャネル推定を行う（ステップＳ４２）。ここでのチャネル推定は、少なくとも自端末局装置宛の信号系列については全て行う。一方、他端末局装置宛に関しては、可能であればチャネル推定を行い、他端末局装置宛の信号の干渉抑制のための受信ウエイト行列を算出する（ステップＳ４３）。基地局装置１０１から空間多重の際に用いた送信ウエイト行列における送信ウエイトベクトルの情報が通知されている場合には、この情報をもとに受信ウエイト行列を選択しても構わない。この様な受信ウエイト行列を受信信号ベクトルに乗算し（ステップＳ４４）、その結果の信号をもとに信号検出処理を行う（ステップＳ４５）。その後、信号検出された信号を合成し、データを再生の後出力する（ステップＳ４６）。
【００４７】
以上の実施形態を説明するため、各種パラメータを特定の条件に仮定して説明を行ったが、当然ながらその他の一般的なパラメータによって実施可能である。すなわち、以上述べた実施形態は全て本発明を例示的に示すものであって限定的に示すものではなく、本発明は他の種々の変形態様及び変更態様で実施することが出来る。従って本発明の範囲は特許請求の範囲及びその均等範囲によってのみ規定されるものである。
【００４８】
以上本発明の一実施形態について説明したが、本発明によれば、マルチユーザＭＩＭＯ通信を行う際の送信指向性制御を行うための送信ウエイト行列を、高次の行列の特異値分解や送信ウエイトベクトルの直交化などの膨大な線形演算を行うことなしに簡易に求めることが可能となる。つまり、マルチユーザＭＩＭＯ通信における処理・演算量の削減を行うことができる。また、この演算の簡易化技術として提案されていた従来のランダムビームフォーミング技術では、送信ウエイトベクトルの組み合わせである送信ウエイト行列を一方的に基地局装置が決定し、この送信ウエイト行列に適合した端末局装置を検索してマルチユーザＭＩＭＯ通信を行なうとしていたために、偶然にもその条件に一致している端末局装置が存在し、且つそれらの端末局装置宛てのデータが送信バッファ内に存在する場合にのみ有効であった。しかし、本発明では、ランダムに生成する送信ウエイトベクトルを、端末局装置から返信される情報に応じて自由に組み合わせて利用することが可能となるため、本技術が有効に適用可能な状況の頻度が飛躍的に増大する。その結果として、マルチユーザＭＩＭＯ通信を安定して有効に利用することが可能となり、周波数資源の効率的な運用が可能となる。
【００４９】
なお上述の各装置は内部に、コンピュータシステムを有していても構わない。そして、上述した処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われることも可能である。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【００５０】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【図面の簡単な説明】
【００５１】
【図１】マルチユーザＭＩＭＯにようる無線通信システムの構成を示す図である。
【図２】基地局装置の送信信号記憶フローを示すである。
【図３】基地局装置の事前処理フローを示す図である。
【図４】基地局装置の送信処理フローを示す図である。
【図５】端末局装置の受信処理フロー（事前処理）を示す図である。
【図６】端末局装置の受信処理フロー（データ受信処理）を示す図である。
【図７】従来技術における基地局装置の送信処理フローを示す図である。
【図８】従来技術における端末局装置の受信処理フロー（事前処理）を示す図である。
【図９】従来技術における端末局装置の受信処理フロー（データ受信処理）を示す図である。
【符号の説明】
【００５２】
１０１・・・基地局装置
１０２〜１０４・・・端末局装置
【図１】
0004708205.tif 000009
【図２】
0004708205.tif 000010
【図３】
0004708205.tif 000011
【図４】
0004708205.tif 000012
【図５】
0004708205.tif 000013
【図６】
0004708205.tif 000014
【図７】
0004708205.tif 000015
【図８】
0004708205.tif 000016
【図９】
0004708205.tif 000017