(57)【要約】【課題】オープンループＳＣＷ−ＭＩＭＯにおいて、より正確なランク数とレートの予測を行うことを可能とする。【解決手段】マッピング部３２は、ストリーム数を１から送信側無線通信装置が備える物理アンテナの数まで仮定し、送信側無線通信装置が備える物理アンテナより送信されるデータの通信状況を示す通信状況情報に基づいて、ストリームと物理アンテナとの対応関係を、仮定したストリーム数毎に決定する。ランク・レート予測部３３は、マッピング部３２で決定したストリームを送信する物理アンテナの通信状況情報に基づいて、送信側無線端末装置が送信するデータのキャパシティを仮定したストリーム数毎に予測する。また、ランク・レート予測部３３は、予測した仮定したストリーム数毎のキャパシティに基づいて、最適なストリーム数を決定する。【選択図】図１
2009152839.tif 000002
【特許請求の範囲】
【請求項１】
オープンループＳＣＷ−ＭＩＭＯシステムに含まれる無線通信装置であって、
ストリーム数を１から送信側無線通信装置が備える物理アンテナの数まで仮定し、送信側無線通信装置が備える前記物理アンテナより送信されるデータの通信状況を示す通信状況情報に基づいて、ストリームと前記物理アンテナとの対応関係を、前記仮定したストリーム数毎に決定するマッピング部と、
前記マッピング部で決定した前記ストリームを送信する前記物理アンテナの前記通信状況情報に基づいて、前記送信側無線端末装置が送信するデータのキャパシティを前記仮定したストリーム数毎に予測する予測部と、
前記予測部が予測した前記仮定したストリーム数毎の前記キャパシティに基づいて、最適な前記ストリーム数を決定するストリーム数決定部と、
を備えたことを特徴とする無線通信装置。
【請求項２】
前記マッピング部は、同一の前記ストリームを送信する前記物理アンテナで構成される集合内の各前記物理アンテナの前記通信状況情報の総和を算出し、各前記集合について算出した前記通信状況情報の前記総和の分散が最も少なくなるように前記物理アンテナで構成される前記集合を決定する
ことを特徴とする請求項１に記載の無線通信装置。
【請求項３】
オープンループＳＣＷ−ＭＩＭＯシステムに含まれる無線通信装置の無線通信方法であって、
ストリーム数を１から送信側無線通信装置が備える物理アンテナの数まで仮定し、送信側無線通信装置が備える前記物理アンテナより送信されるデータの通信状況を示す通信状況情報に基づいて、ストリームと前記物理アンテナとの対応関係を、前記仮定したストリーム数毎に決定するマッピングステップと、
前記マッピングステップで決定した前記ストリームを送信する前記物理アンテナの前記通信状況情報に基づいて、前記送信側無線端末装置が送信するデータのキャパシティを前記仮定したストリーム数毎に予測する予測ステップと、
前記予測ステップで予測した前記仮定したストリーム数毎の前記キャパシティに基づいて、最適な前記ストリーム数を決定するストリーム数決定ステップと、
を含むことを特徴とする無線通信方法。
【請求項４】
前記マッピングステップは、同一の前記ストリームを送信する前記物理アンテナで構成される集合内の各前記物理アンテナの前記通信状況情報の総和を算出し、各前記集合について算出した前記通信状況情報の前記総和の分散が最も少なくなるように前記物理アンテナで構成される前記集合を決定する
ことを特徴とする請求項３に記載の無線通信方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、オープンループＳＣＷ−ＭＩＭＯ技術を使用する無線通信装置および無線通信方法に関する。
【背景技術】
【０００２】
ＭＩＭＯ（Ｍｕｌｔｉｐｌｅ−Ｉｎｐｕｔ  Ｍｕｌｔｉｐｌｅ−Ｏｕｔｐｕｔ、多入力・多出力）技術は、送受信局に複数のアンテナを装備し、同一周波数干渉環境下でも異なる信号系列の同時空間多重化伝送を可能とする技術である（例えば、非特許文献１参照）。これにより、伝送路を空間的に多重化して伝送レートを向上させたり、あるいは時空間符号化による送受信ダイバーシチで誤り率を低減させたりといった効果が得られる。そのため、ＭＩＭＯ技術は様々な高速無線通信システムに採用されている。
【０００３】
ＭＩＭＯは、送受信局が同時に符号／復号、変調／復調するコードワードの数によってＳｉｎｇｌｅ  Ｃｏｄｅ  Ｗｏｒｄ（シングルコードワード、ＳＣＷ）ＭＩＭＯとＭｕｌｔｉ−Ｃｏｄｅ  Ｗｏｒｄ（マルチコードワード、ＭＣＷ）ＭＩＭＯとに分かれる。ＳＣＷ−ＭＩＭＯでは、送受信局が符号／復号、変調／復調するデータは１つの（同一の）コードワードである。
【０００４】
図３は、ＳＣＷ−ＭＩＭＯを使用する送信局の構成例を示した図である。送信局は、ターボエンコーダ部１１と、ＱＡＭ（Ｑｕａｄｒａｔｕｒｅ  Ａｍｐｌｉｔｕｄｅ  Ｍｏｄｕｌａｔｉｏｎ、直交振幅変調）部１２と、シリアル／パラレル変換部１３と、マッピング部１４と、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ  Ｆｒｅｑｕｅｎｃｙ  Ｄｉｖｉｓｉｏｎ  Ｍｕｌｔｉｐｌｅｘｉｎｇ、直交波周波数分割多重）変調部１５と物理アンテナ１６と、受信部１７と、レート予測部１８とを備える。なお、送信局は、ＯＦＤＭ変調部１５と物理アンテナ１６とをそれぞれＭＴ個含む。
【０００５】
図３を参照して、送信局がデータを送信する手順について説明する。初めに、ターボエンコーダ部１１は、送信対象のデータを誤り訂正符号化する。続いて、ＱＡＭ部１２は、誤り訂正符号化されたデータより変調シンボルを作成する。続いて、受信部１７が受信局より受信した予測ランク情報に基づいて、シリアル／パラレル変換部１３は、変調シンボルをシリアル／パラレル変換し、独立に送信するＭ個のストリームを生成する。ＳＣＷ−ＭＩＭＯにおいて、送信局は、送信するデータを分割し、受信局に対して並列に送信する。この送信するデータを分割する数をランク数と呼ぶ。また、受信局に対して送信する並列数をストリーム数と呼ぶ。ランク数とストリーム数は同じ数である。なお、予測ランク情報には、受信局が予測した一番キャパシティ（通信容量）が良いと想定されるランク数が含まれている。
【０００６】
続いて、マッピング部１４は、ストリーム毎に、ストリームをどのＯＦＤＭ変調部１５および物理アンテナ１６に関連付けるのか決定する。この関連付けをマッピングと呼ぶ。続いてＯＦＤＭ変調部１５は、マッピングされたストリームをＯＦＤＭ変調する。物理アンテナ１６は、ＯＦＤＭ変調されたストリームを受信局に対して送信する。
【０００７】
また、ＭＩＭＯシステムは、Ｐｒｅ−ｃｏｄｅｒ（プリコーダ）を使用するかどうかによって、大きくＣｌｏｓｅｄ  Ｌｏｏｐ  ＭＩＭＯ（クローズド・ループ・マイモ）システムとＯｐｅｎ  Ｌｏｏｐ  ＭＩＭＯ（オープン・ループ・マイモ）システムに分かれる。
【０００８】
Ｃｌｏｓｅｄ  Ｌｏｏｐ  ＭＩＭＯシステムは、予め用意したＰｒｅ−ｃｏｄｅｒによってＢｅａｍ  Ｆｏｒｍｉｎｇ（ビームフォーミング）を作成し、通信を行う。例えば、物理アンテナ１６に｛０，１｝の要素からなる送信アレー重みを乗じることにより、ストリームを送信する物理アンテナ１６を選択する。
【０００９】
Ｏｐｅｎ  Ｌｏｏｐ  ＭＩＭＯシステムは、Ｍ個のストリームのうち一部のストリームをコピーして、物理アンテナ１６にマッピングし、通信を行う。これにより、Ｏｐｅｎ  Ｌｏｏｐ  ＭＩＭＯシステムでは送信ダイバシチゲインを得られる。
【００１０】
Ｏｐｅｎ  Ｌｏｏｐ  ＭＩＭＯシステムでは、最大のキャパシティを得るために、受信局でランクとレートの予測を行い、送信局にフィードバックしている。なお、ランクは独立に送信するストリームを意味する。また、レートはいわゆるＭＣＳ（エム・シー・エス、Ｍｏｄｕｌａｔｉｏｎ  ａｎｄ  Ｃｏｄｉｎｇ  Ｓｃｈｅｍｅ）であり、送信コードワードの変調方式および符号化率を意味する。また、受信局が行うランクとレートの予測は、受信局で推定した空間チャネル特性Ｈに基づいて行う。
【００１１】
図４は、ＳＣＷ−ＭＩＭＯを使用する受信局の構成例を示した図である。受信局は、物理アンテナ２１と、ＯＦＤＭ復調部２２と、データ分離部２３と、パラレル／シリアル変換部２４と、対数尤度比（ＬＬＲ、Ｌｏｇ−Ｌｉｋｅｌｉｈｏｏｄ  Ｒａｔｉｏ）算出部２５と、ターボデコーダ部２６と、ランク・ＣＱＩ（Ｃｈａｎｎｅｌ  Ｑｕａｌｉｔｙ  Ｉｎｄｉｃａｔｉｏｎ）予測部２７とを備える。なお、受信局は、物理アンテナ２１とＯＦＤＭ復調部２２とをそれぞれＭＲ個含む。
【００１２】
図４を参照して、受信局がデータを受信する手順について説明する。初めに、物理アンテナ２１は、送信局より送信されたデータを受信する。続いて、ＯＦＤＭ復調部２２は、物理アンテナ２１が受信したデータをＯＦＤＭ復調する。続いて、データ分離部２３は、ＭＭＳＥ（Ｍｉｎｉｍｕｍ−Ｍｅａｎ  Ｓｑｕａｒｅ−Ｅｒｒｏｒ）空間フィルタにより、各ストリームのデータを分離する。続いて、パラレル／シリアル変換部２４は、データ分離部２３が分離したデータをパラレル／シリアル変換する。続いて、対数尤度比算出部２５は、シリアル変換されたデータの誤り訂正を行う。続いて、ターボデコーダ部２６は、誤り訂正を行った信号をオリジナルの信号に復号する。
【００１３】
また、データ分離部２３が出力するＭ個の出力は、送信局より送信されたストリームデータの推定値を表す。なお、送信局の物理アンテナ１６毎に用意されるパイロット信号の送受信により、送信局の物理アンテナ１６と受信局の物理アンテナ２１との間のチャネル数の推定ができる。図３および図４の例では、Ｍ＝ＭＴである。
【００１４】
また、各ストリームデータの推定値と、無線通信に影響を与える雑音（熱雑音と他のストリームからの干渉）の電力比をＰｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲ（Ｓ／Ｎ比、Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  ｒａｔｉｏ）として定義する。ランク・ＣＱＩ予測部２７は、Ｐｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲに基づいて、ランクとレートの予測を行う。
【００１５】
以下、Ｐｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲの求め方を説明する。なお、以下の数式ではマトリックスのサイズを下付きの中括弧｛行×列｝で表記する。送信局の物理アンテナ１６の数ＭＴと受信局の物理アンテナ２１の数ＭＲ間のチャネルインパルス行列は（１）式と仮定する。
【００１６】
【数１】
2009152839.tif 000003
【００１７】
また、送信信号ベクトルを（２）式、受信信号ベクトルを（３）式、データ分離部２３からの復号信号ベクトルを（４）式、送信信号総電力をρ、雑音電力を（５）式、ＯＦＤＭシンボルＩＤをｎ、サブキャリアＩＤをｋ、受信アンテナ応答の単位行列を（６）式とすれば、（７）式および（８）式となる。ここで、（８）式において上付きのＨはエルミート転置を表す。また、（９）式はＭＭＳＥ空間フィルタのウェートを表す。
【００１８】
【数２】
2009152839.tif 000004
【００１９】
【数３】
2009152839.tif 000005
【００２０】
【数４】
2009152839.tif 000006
【００２１】
【数５】
2009152839.tif 000007
【００２２】
【数６】
2009152839.tif 000008
【００２３】
【数７】
2009152839.tif 000009
【００２４】
【数８】
2009152839.tif 000010
【００２５】
【数９】
2009152839.tif 000011
【００２６】
また、受信信号ベクトルｒは（１０）式となる。（１０）式において、Ｎは雑音を表す。
【００２７】
【数１０】
2009152839.tif 000012
【００２８】
また、ＭＭＳＥ空間フィルタの要素を（１１）式とし、（１２）式で示すチャネルインパルス行列の要素を（１３）式とすると、（８）式を（１４）式と書き換えることができる。このとき、ｑ＝１，２，・・・，ＭＴであり、送信局の物理アンテナ１６のｑ番目から送信された送信信号ベクトルはＳｑであり、物理アンテナ１６のｑ番目から送信され、受信局で受信した復号信号ベクトルの推定値は（１５）式である。
【００２９】
【数１１】
2009152839.tif 000013
【００３０】
【数１２】
2009152839.tif 000014
【００３１】
【数１３】
2009152839.tif 000015
【００３２】
【数１４】
2009152839.tif 000016
【００３３】
【数１５】
2009152839.tif 000017
【００３４】
送信局の物理アンテナ１６のｑ番目から送信された送信信号のＳＩＮＲ（Ｓｉｇｎａｌ  ｔｏ  Ｉｎｔｅｒｆｅｒｅｎｃｅ  ｐｌｕｓ  Ｎｏｉｓｅ  ｐｏｗｅｒ  Ｒａｔｉｏ、受信信号電力対干渉および雑音電力比）値であるＳＩＮＲｑは、（１４）式に基づいて計算すると（１６）式となる。このとき、ｑ＝１，・・・，ＭＴである。
【００３５】
【数１６】
2009152839.tif 000018
【００３６】
ランク（ストリーム）の予測方法としては、例えば、６４ＱＡＭの限界性能カーブＣ６４ＱＡＭに基づいて、ストリーム毎にキャパシティ（Ｃａｐ）と、有効ＳＩＮＲ（ＳＩＮＲＥｆｆ、Ｅｆｆｅｃｔｉｖｅ  ＳＩＮＲ）値を算出し、算出した値に基づいてランク数を予測する方法が考えられる。
【００３７】
以下、ランク数を予測する方法の具体例を説明する。Ｃａｐは（１７）式となる。また、ＳＩＮＲＥｆｆは（１８）式となる。このとき、ｑ＝１，・・・，ＭＴであり、Ｋは総サブキャリア数であり、Ｎは総ＯＦＤＭシンボル数である。また、Ｃ６４ＱＡＭの逆関数は（１９）式となる。
【００３８】
【数１７】
2009152839.tif 000019
【００３９】
【数１８】
2009152839.tif 000020
【００４０】
【数１９】
2009152839.tif 000021
【００４１】
図５は６４ＱＡＭの限界性能カーブＣ６４ＱＡＭを示した図である。図の縦軸はキャパシティ（ｂｉｔ／Ｈｚ／ｓ）の値を示し、横軸はＳＮＲ値を示している。図示する例では、Ｃ６４ＱＡＭはＳＮＲ値が大きくなるにつれ、キャパシティの値も大きくなっている。
【００４２】
（１６）式、（１７）式、（１８）式に基づいて、ランクが１からＭＴまでの各ランクにおけるキャパシティ（Ｆｉｎａｌ＿Ｃａｐａｃｉｔｙ）の予測を行う。
【００４３】
（ランク＝１と仮定した場合のキャパシティ予測）
ランク＝１の場合、ＭＴ個のストリームデータの内１つのストリームデータを取得することで送信データを復号できるため、最大のＳＩＮＲ値に基づいて予測キャパシティ値を算出する。なお、予測キャパシティ値の算出に使用するＳＩＮＲＥｆｆをＦｉｎａｌ＿ＳＩＮＲＥｆｆとする。
【００４４】
（ステップＳ１０１）Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝１）を１番良いＳＩＮＲＥｆｆｑとする。ｑ＝１，２，・・・，ＭＴである。
【００４５】
（ステップＳ１０２）ステップＳ１０１で決定したＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝１）と、予め用意されているＭＣＳ（エム・シー・エス、Ｍｏｄｕｌａｔｉｏｎ  ａｎｄ  Ｃｏｄｉｎｇ  Ｓｃｈｅｍｅ）のＰＥＲ（Ｐａｃｋｅｔ  Ｅｒｒｏｒ  Ｒａｔｅ、パケットの誤り率）カーブとに基づいて、要求するパケットの誤り率を満たすことができる最大のレートを選択する。
【００４６】
（ステップＳ１０３）ステップＳ１０２で決定したレートの変調方式における限界性能カーブと、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝１）とに基づいて、Ｃａｐａｃｉｔｙ（ランク＝１）を決定する。なお、限界性能カーブは図５に示した通りである。
【００４７】
（ステップＳ１０４）ステップＳ１０３で決定したＣａｐａｃｉｔｙ（ランク＝１）をＦｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝１）とする。上述した手順により、ランク＝１と仮定した場合のキャパシティの予測結果を算出することができる。
【００４８】
（ランク＝Ｍ、１＜Ｍ＜ＭＴと仮定した場合のキャパシティ予測）
ランク＝Ｍ、１＜Ｍ＜ＭＴと仮定した場合、ＭＴ個のストリームデータの内、Ｍ番目に良いＳＩＮＲ値に基づいて予測キャパシティ値を算出する。なお、予測キャパシティ値の算出に使用するＳＩＮＲＥｆｆをＦｉｎａｌ＿ＳＩＮＲＥｆｆとする。
【００４９】
（ステップＳ２０１）Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝Ｍ）をＭ番目に良いＳＩＮＲＥｆｆｑとする。ｑ＝１，２，・・・，ＭＴである。
【００５０】
（ステップＳ２０２）ステップＳ２０１で決定したＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝Ｍ）と、予め用意されているＭＣＳのＰＥＲカーブとに基づいて、要求するパケットの誤り率を満たすことができる最大のレートを選択する。
【００５１】
（ステップＳ２０３）ステップＳ２０２で決定したレートの変調方式における限界性能カーブと、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝Ｍ）とに基づいて、Ｃａｐａｃｉｔｙ（ランク＝Ｍ）を決定する。なお、限界性能カーブは図５に示した通りである。
【００５２】
（ステップＳ２０４）ランクがＭの場合、データをＭ並列で送受信するため、ステップＳ２０３で決定したＣａｐａｃｉｔｙ（ランク＝Ｍ）に基づいてＦｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝Ｍ）＝Ｍ×Ｃａｐａｃｉｔｙ（ランク＝Ｍ）とする。上述した手順により、ランク＝Ｍと仮定した場合のキャパシティの予測結果を算出することができる。
【００５３】
（ランク＝ＭＴと仮定した場合のキャパシティ予測）
ランク＝ＭＴの場合、ＭＴ個のストリームデータの内１つでも欠けると送信データを復号することが出来ないため、最小のＳＩＮＲ値に基づいて予測キャパシティ値を算出する。なお、予測キャパシティ値の算出に使用するＳＩＮＲＥｆｆをＦｉｎａｌ＿ＳＩＮＲＥｆｆとする。
【００５４】
（ステップＳ３０１）Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝ＭＴ）を１番低いＳＩＮＲＥｆｆｑとする。ｑ＝１，２，・・・，ＭＴである。
【００５５】
（ステップＳ３０２）ステップＳ３０１で決定したＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝ＭＴ）と、予め用意されているＭＣＳのＰＥＲカーブとに基づいて、要求するパケットの誤り率を満たすことができる最大のレートを選択する。
【００５６】
（ステップＳ３０３）ステップＳ３０２で決定したレートの変調方式における限界性能カーブと、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝ＭＴ）とに基づいて、Ｃａｐａｃｉｔｙ（ランク＝ＭＴ）を決定する。なお、限界性能カーブは図５に示した通りである。
【００５７】
（ステップＳ３０４）ランクがＭＴの場合、データをＭＴ並列で送受信するため、ステップＳ３０３で決定したＣａｐａｃｉｔｙ（ランク＝ＭＴ）よりＦｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝ＭＴ）＝ＭＴ×Ｃａｐａｃｉｔｙ（ランク＝ＭＴ）とする。上述した手順により、ランク＝Ｍと仮定した場合のキャパシティの予測結果を算出することができる。
【００５８】
上述した手順により求めたＦｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝１）、Ｆｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝２）、・・・、Ｆｉｎａｌ＿Ｃａｐａｃｉｔｙ（ランク＝ＭＴ）のうち、最大の値となるランク数をランク予測の結果とする。また、決定したランク数に対応するレートをレート予測の結果（ＣＱＩ）として、送信局にフィードバックする。
【００５９】
【非特許文献１】ＩＥＴＦ  ＲＦＣ３６２６：Ｏｐｔｉｍｉｚｅｄ  Ｌｉｎｋ  Ｓｔａｔｅ  Ｒｏｕｔｉｎｇ  Ｐｒｏｔｏｃｏｌ
【発明の開示】
【発明が解決しようとする課題】
【００６０】
上述した通り、オープンループＳＣＷ−ＭＩＭＯにおいて、ランク数と送信局の物理アンテナ数とが一致しない場合もありえる。しかしながら、上述したランク数の決定方法では、ランクと送信局の物理アンテナとのマッピングについては考慮していない。ランクと物理アンテナとの組み合わせによってはレートが異なるため、正確なランク数とレートの予測を行うことが出来ないという問題があった。
【００６１】
本発明は上記の課題を解決するためになされたものであり、オープンループＳＣＷ−ＭＩＭＯにおいて、より正確なランク数とレートの予測を行うことが可能な無線通信装置を提供することを目的とする。
【課題を解決するための手段】
【００６２】
本発明は、オープンループＳＣＷ−ＭＩＭＯシステムに含まれる無線通信装置であって、ストリーム数を１から送信側無線通信装置が備える物理アンテナの数まで仮定し、送信側無線通信装置が備える前記物理アンテナより送信されるデータの通信状況を示す通信状況情報に基づいて、ストリームと前記物理アンテナとの対応関係を、前記仮定したストリーム数毎に決定するマッピング部と、前記マッピング部で決定した前記ストリームを送信する前記物理アンテナの前記通信状況情報に基づいて、前記送信側無線端末装置が送信するデータのキャパシティを前記仮定したストリーム数毎に予測する予測部と、前記予測部が予測した前記仮定したストリーム数毎の前記キャパシティに基づいて、最適な前記ストリーム数を決定するストリーム数決定部と、を備えたことを特徴とする無線通信装置である。
【００６３】
また、本発明の無線通信装置において、前記マッピング部は、同一の前記ストリームを送信する前記物理アンテナで構成される集合内の各前記物理アンテナの前記通信状況情報の総和を算出し、各前記集合について算出した前記通信状況情報の前記総和の分散が最も少なくなるように前記物理アンテナで構成される前記集合を決定することを特徴とする。
【００６４】
また、本発明は、オープンループＳＣＷ−ＭＩＭＯシステムに含まれる無線通信装置の無線通信方法であって、ストリーム数を１から送信側無線通信装置が備える物理アンテナの数まで仮定し、送信側無線通信装置が備える前記物理アンテナより送信されるデータの通信状況を示す通信状況情報に基づいて、ストリームと前記物理アンテナとの対応関係を、前記仮定したストリーム数毎に決定するマッピングステップと、前記マッピングステップで決定した前記ストリームを送信する前記物理アンテナの前記通信状況情報に基づいて、前記送信側無線端末装置が送信するデータのキャパシティを前記仮定したストリーム数毎に予測する予測ステップと、前記予測ステップで予測した前記仮定したストリーム数毎の前記キャパシティに基づいて、最適な前記ストリーム数を決定するストリーム数決定ステップと、を含むことを特徴とする無線通信方法である。
【００６５】
また、本発明の無線通信方法において、前記マッピングステップは、同一の前記ストリームを送信する前記物理アンテナで構成される集合内の各前記物理アンテナの前記通信状況情報の総和を算出し、各前記集合について算出した前記通信状況情報の前記総和の分散が最も少なくなるように前記物理アンテナで構成される前記集合を決定することを特徴とする。
【発明の効果】
【００６６】
本発明によれば、オープンループＳＣＷ−ＭＩＭＯにおいて、より正確なランク数とレートの予測を行うことができる。
【発明を実施するための最良の形態】
【００６７】
以下、本発明の一実施形態について説明する。本実施形態では、オープンループＳＣＷ−ＭＩＭＯ技術を使用する無線通信システムを想定する。図１は本発明の一実施形態における無線通信装置の受信部分に関する構成を示した図である。
【００６８】
無線通信装置の受信部分は、物理アンテナ２１と、ＯＦＤＭ復調部２２と、データ分離部２３と、パラレル／シリアル変換部２４と、対数尤度比（ＬＬＲ、Ｌｏｇ−Ｌｉｋｅｌｉｈｏｏｄ  Ｒａｔｉｏ）算出部２５と、ターボデコーダ部２６と、ＳＮＲ（Ｓ／Ｎ比、Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  ｒａｔｉｏ）計算部３１と、マッピング部３２と、ランク・レート予測部３３とを備える。なお、無線通信装置の受信部分は、物理アンテナ２１とＯＦＤＭ復調部２２とをそれぞれＭＲ個含む。また、無線通信装置の送信部分は、従来技術と同様の構成であり、図３に示した構成と同様である。
【００６９】
図１を参照して、無線通信装置の受信部分がテータを受信する手順について説明する。初めに、物理アンテナ２１は、他の無線通信装置より送信されたデータを受信する。続いて、ＯＦＤＭ復調部２２は、物理アンテナ２１が受信したデータをＯＦＤＭ復調する。続いて、データ分離部２３は、ＭＭＳＥ（Ｍｉｎｉｍｕｍ−Ｍｅａｎ  Ｓｑｕａｒｅ−Ｅｒｒｏｒ）空間フィルタにより、各ストリームのデータを分離する。続いて、パラレル／シリアル変換部２４は分離した各ストリームのデータをパラレル／シリアル変換する。続いて、対数尤度比算出部２５は、シリアル変換されたデータの誤り訂正を行う。続いて、ターボデコーダ部２６は、誤り訂正を行った信号をオリジナルの信号に復号する。
【００７０】
また、データ分離部２３が出力するＭ個の出力は、送信局より送信されたストリームデータの推定値を表す。なお、送信局の物理アンテナ１６毎に用意されるパイロット信号の送受信により、送信局の物理アンテナ１６と受信局の物理アンテナ２１との間のチャネル数の推定ができる。本実施形態では、Ｍ＝ＭＴである。
【００７１】
また、各ストリームデータの推定値と、無線通信に影響を与える雑音（熱雑音と他のストリームからの干渉）の電力比をＰｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲ（Ｓ／Ｎ比、Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  ｒａｔｉｏ）として定義する。ＳＮＲ計算部３１は、Ｐｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲを算出する。ランク・レート予測部３３は、ＳＮＲ計算部３１が算出したＰｏｓｔ−ｐｒｏｃｅｓｓｉｎｇ  ＳＮＲに基づいて、ランクとレートの予測を行う。
【００７２】
次に、図２を参照して、本実施形態におけるランクおよびレートの予測方法について説明する。図２は、本実施形態におけるランクおよびレートの測定手順を示すフローチャートである。なお、以下の数式ではマトリックスのサイズを下付きの中括弧｛行×列｝で表記する。
【００７３】
（ステップＳ１１０１）従来技術と同様に、（８）式および（１６）式に基づいて、ＳＮＲ計算部３１はＳＩＮＲｑ（ｎ，ｋ）を算出する。その後、ステップＳ１１０２に進む。
【００７４】
（ステップＳ１１０２）従来技術と同様に、（１７）式および（１８）式に基づいて、ＳＮＲ計算部３１はＳＩＮＲＥｆｆｑ（ｎ，ｋ）を算出する。なお、ｑ＝１，２，・・・，ＭＴである。また、ＳＮＲ計算部３１は算出した値をマッピング部３２に入力する。その後、ステップＳ１１０３に進む。
【００７５】
（ステップＳ１１０３）ランク＝１と仮定した場合のキャパシティ予測を行うためｉｎｉｔ＝１とする。その後、ステップＳ１１０４に進む。なお、本実施形態ではｉｎｉｔ値＝ランク数である。
【００７６】
ステップＳ１１０４〜ステップＳ１１１４の処理は、ランクＭが、Ｍ＝１の場合と、１＜Ｍ＜ＭＴの場合と、Ｍ＝ＭＴの場合とで異なる。以下、ランクＭがＭ＝１の場合と、１＜Ｍ＜ＭＴの場合と、Ｍ＝ＭＴの場合との処理の手順について説明する。
【００７７】
（ランクＭ＝１の場合）
（ステップＳ１１０４）マッピング部３２は、ストリームと物理アンテナ１６とのマッピングを行う。ランクＭ＝１の場合は、全ての物理アンテナ１６は同一のストリームデータを送信するため、マッピングは１通りとなる。その後、ステップＳ１１０５に進む。
【００７８】
（ステップＳ１１０５）ステップＳ１１０４で決定したマッピングにおいて、データ分離部２３はストリームと物理アンテナ１６間のチャネル行列Ｈ´を算出する。その後、ステップＳ１１０６に進む。なお、チャネル行列Ｈ´は（２０）式より算出することができる。
【００７９】
【数２０】
2009152839.tif 000022
【００８０】
（ステップＳ１１０６）データ分離部２３は、ＭＭＳＥウェートω´を算出する。続いて、ステップＳ１１０７に進む。なお、ＭＭＳＥウェートω´は（２１）式、（２２）式より算出することができる。
【００８１】
【数２１】
2009152839.tif 000023
【００８２】
【数２２】
2009152839.tif 000024
【００８３】
（ステップＳ１１０７）ＳＮＲ計算部３１は、ランクＭ＝1の場合のＳＩＮＲ値を算出する。続いて、ステップＳ１１０８に進む。なお、チャネルインパルス行列を（２３）式、ＭＭＳＥ空間フィルタの要素を（２４）式、チャネルインパルス行列の要素を（２５）式とすると、ＳＩＮＲ値は（２６）式より算出することができる。また、ランク数は１であるためｑ＝１である。
【００８４】
【数２３】
2009152839.tif 000025
【００８５】
【数２４】
2009152839.tif 000026
【００８６】
【数２５】
2009152839.tif 000027
【００８７】
【数２６】
2009152839.tif 000028
【００８８】
（ステップＳ１１０８）ステップＳ１１０７で算出したＳＩＮＲに基づいて、ランク・レート予測部３３はキャパシティＣａｐとＳＩＮＲＥｆｆとを算出する。続いて、ステップＳ１１０９に進む。なお、キャパシティは（２７）式、ＳＩＮＲＥｆｆは（２８）式より算出することができる。また、ランク数は１であるためｑ＝１である。
【００８９】
【数２７】
2009152839.tif 000029
【００９０】
【数２８】
2009152839.tif 000030
【００９１】
（ステップＳ１１０９）ランク・レート予測部３３はｉｎｉｔが１であるか否か判断する。ランクＭ＝1の場合はｉｎｉｔ＝１であるため、ステップＳ１１１０に進む。
【００９２】
（ステップＳ１１１０）ランク・レート予測部３３はランク＝１の場合のＦｉｎａｌ＿Ｃａｐ（ランク＝１）およびＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝１）を求める。その後、ステップＳ１１１３に進む。なお、（２９）式に示す通り、Ｆｉｎａｌ＿Ｃａｐ（ランク＝１）はステップＳ１１０８で算出したＣａｐである。また、（３０）式に示す通り、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝１）はステップＳ１１０８で算出したＳＩＮＲＥｆｆである。
【００９３】
【数２９】
2009152839.tif 000031
【００９４】
【数３０】
2009152839.tif 000032
【００９５】
（ステップＳ１１１３）ランク・レート予測部３３はｉｎｉｔがＭＴ以上か否か判断する。ｉｎｉｔがＭＴ以上と判断した場合ステップＳ１１１５に進み、それ以外はステップＳ１１１４に進む。
（ステップＳ１１１４）ｉｎｉｔの値に１を加算する。その後ステップＳ１１０４に戻る。
【００９６】
ランクＭが１の場合における、ステップＳ１１０４〜ステップＳ１１１４の処理の手順は以上の通りである。
【００９７】
（ランクＭが１＜Ｍ＜１の場合）
（ステップＳ１１０４）マッピング部３２は、ストリームと物理アンテナ１６とのマッピングを行う。ランクＭが１＜Ｍ＜１の場合、ステップＳ１１０３で入力されたＳＩＮＲＥｆｆｑに基づいて、ＭＴ個の物理アンテナ１６をＭ個の集合とする。このときの条件として、各集合に属する物理アンテナ１６のＳＩＮＲＥｆｆの総和がなるべく等しくなるよう（総和の分散がなるべく小さくなるよう）に集合を決定する。なお、ＳＩＮＲＥｆｆｑは真数とする。また、集合Ｍに属する物理アンテナ１６に、集合Ｍ毎に異なるストリームをマッピングする。その後、ステップＳ１１０５に進む。
【００９８】
ＳＣＷ−ＭＩＭＯでは、送信側無線通信装置が、送信するデータを分割し、並列に送信するため、分割した分割データ全てが届かなければ受信側無線通信装置は送信されたデータを復元できない。よって、分割されたうちの１つの分割データのみの送信が遅れた場合においても、その遅れた分割データの影響でデータ転送のキャパシティは悪くなってしまう。このことにより、本実施形態では各集合に属する物理アンテナ１６のＳＩＮＲＥｆｆの総和の分散がなるべく小さくなるようにし、分割した全ての分割データがより速く受信側無線通信装置で揃うようにストリームと物理アンテナ１６とのマッピングを行う。
【００９９】
ステップＳ１１０４でのマッピング方法について具体例を用いて説明する。例えば、ＭＴ＝４すなわち物理アンテナが４個であり、Ｍ＝２すなわちストリームが２つであり、物理アンテナ１のＳＩＮＲＥｆｆ１＜物理アンテナ２のＳＩＮＲＥｆｆ２＜物理アンテナ３のＳＩＮＲＥｆｆ３＜物理アンテナ４のＳＩＮＲＥｆｆ４であり、ＳＩＮＲＥｆｆ１＋ＳＩＮＲＥｆｆ４≒ＳＩＮＲＥｆｆ２＋ＳＩＮＲＥｆｆ３である場合、マッピング部３２は、ストリームの一方を物理アンテナ１と物理アンテナ４にマッピングし、もう一方のストリームを物理アンテナ２と物理アンテナ３にマッピングする。
【０１００】
（ステップＳ１１０５）ステップＳ１１０４で決定したマッピングにおいて、データ分離部２３はストリームと物理アンテナ２１間のチャネル行列Ｈ´を算出する。その後、ステップＳ１１０６に進む。なお、上記の具体例におけるチャネル行列Ｈ´は（３１）式より算出することができる。
【０１０１】
【数３１】
2009152839.tif 000033
【０１０２】
（ステップＳ１１０６）データ分離部２３は、ＭＭＳＥウェートω´を算出する。続いて、ステップＳ１１０７に進む。なお、ＭＭＳＥウェートω´は（３２）式、（３３）式より算出することができる。
【０１０３】
【数３２】
2009152839.tif 000034
【０１０４】
【数３３】
2009152839.tif 000035
【０１０５】
（ステップＳ１１０７）ＳＮＲ計算部３１は、ランクＭが１＜Ｍ＜１の場合のＳＩＮＲ値を算出する。続いて、ステップＳ１１０８に進む。なお、チャネルインパルス行列を（３４）式、ＭＭＳＥ空間フィルタの要素を（３５）式、チャネルインパルス行列の要素を（３６）式とすると、ＳＩＮＲ値は（３７）式より算出することができる。また、ｑ＝１，２，・・・，ＭＴである。
【０１０６】
【数３４】
2009152839.tif 000036
【０１０７】
【数３５】
2009152839.tif 000037
【０１０８】
【数３６】
2009152839.tif 000038
【０１０９】
【数３７】
2009152839.tif 000039
【０１１０】
（ステップＳ１１０８）ステップＳ１１０７で算出したＳＩＮＲに基づいて、ランク・レート予測部３３はキャパシティＣａｐとＳＩＮＲＥｆｆとを算出する。続いて、ステップＳ１１０９に進む。なお、キャパシティは（３８）式、ＳＩＮＲＥｆｆは（３９）式より算出することができる。また、ｑ＝１，２，・・・，ＭＴである。
【０１１１】
【数３８】
2009152839.tif 000040
【０１１２】
【数３９】
2009152839.tif 000041
【０１１３】
（ステップＳ１１０９）ランク・レート予測部３３はｉｎｉｔが１であるか否か判断する。１であると判断した場合ステップＳ１１１０に進み、それ以外はステップＳ１１１１に進む。現時点ではｉｎｉｔ＞１であるため、ステップＳ１１１１に進む。
【０１１４】
（ステップＳ１１１１）ステップＳ１１０４で各集合に属する物理アンテナ１６のＳＩＮＲＥｆｆの総和がなるべく等しくなるようにし、ストリームと物理アンテナのマッピングを行ったため、各ストリームはＳＩＮＲＥｆｆの平均値で受信することが出来ると考えられる。よって、ステップＳ１１０８で算出したキャパシティＣａｐとＳＩＮＲＥｆｆに基づいて、ランク・レート予測部３３は、ＳＩＮＲＥｆｆの平均値ＳＩＮＲＥｆｆａｖｇを算出する。また、ランク・レート予測部３３は、算出したＳＩＮＲＥｆｆａｖｇに基づいて、キャパシティＣａｐａｖｇを算出する。なお、ＳＩＮＲＥｆｆａｖｇは（４０）式、Ｃａｐａｖｇは（４１）式より算出することができる。また、ｑ＝１，２，・・・，ＭＴである。
【０１１５】
【数４０】
2009152839.tif 000042
【０１１６】
【数４１】
2009152839.tif 000043
【０１１７】
（ステップＳ１１１２）ランク・レート予測部３３はランクＭが１＜Ｍ＜１の場合のＦｉｎａｌ＿Ｃａｐ（ランク＝Ｍ）およびＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝Ｍ）を求める。その後、ステップＳ１１１３に進む。なお、（４２）式に示すとおり、Ｆｉｎａｌ＿Ｃａｐ（ランク＝Ｍ）はＣａｐａｖｇのＭ倍の値となる。また、（４３）式に示すとおり、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝Ｍ）は、ＳＩＮＲＥｆｆａｖｇとなる。
【０１１８】
【数４２】
2009152839.tif 000044
【０１１９】
【数４３】
2009152839.tif 000045
【０１２０】
（ステップＳ１１１３）ランク・レート予測部３３はｉｎｉｔがＭＴ以上か否か判断する。ｉｎｉｔがＭＴ以上と判断した場合ステップＳ１１１５に進み、それ以外はステップＳ１１１４に進む。
（ステップＳ１１１４）ｉｎｉｔの値に１を加算する。その後ステップＳ１１０４に戻る。
【０１２１】
ランクＭが１＜Ｍ＜１の場合における、ステップＳ１１０４〜ステップＳ１１１４の処理の手順は以上の通りである。
【０１２２】
（ランクＭ＝ＭＴの場合）
（ステップＳ１１０４）マッピング部３２は、ストリームと物理アンテナ１６とのマッピングを行う。ランクＭ＝ＭＴの場合、ストリームと物理アンテナ１６とは１対１の関係となるため、１つのストリームを１つの物理アンテナ１６にマッピングする。その後、ステップＳ１１０５に進む。
【０１２３】
（ステップＳ１１０５）ステップＳ１１０４で決定したマッピングにおいて、データ分離部２３はストリームと物理アンテナ２１間のチャネル行列Ｈ´を算出する。なお、ストリームと物理アンテナ１６とは１対１の関係であるため、ストリームと物理アンテナ２１間のチャネル行列Ｈ´は、従来技術における送信局の物理アンテナ１６の数ＭＴと受信局の物理アンテナ２１の数ＭＲ間のチャネルインパルス行列Ｈと同様に求めることが出来る。その後、ステップＳ１１０６に進む。
【０１２４】
（ステップＳ１１０６）データ分離部２３は、ＭＭＳＥウェートω´を算出する。続いて、ステップＳ１１０７に進む。なお、ストリームと物理アンテナ１６とは１対１の関係であるため、従来技術におけるＭＭＳＥウェートωと同様にＭＭＳＥウェートω´を求めることが出来る。その後、ステップＳ１１０７に進む。
【０１２５】
（ステップＳ１１０７）ＳＮＲ計算部３１は、ランクＭ＝ＭＴの場合のＳＩＮＲ値を算出する。続いて、ステップＳ１１０８に進む。なお、ストリームと物理アンテナ１６とは１対１の関係であるため、従来技術におけるＳＩＮＲ値と同様にランクＭ＝ＭＴの場合のＳＩＮＲ値を求めることが出来る。
【０１２６】
（ステップＳ１１０８）ステップＳ１１０７で算出したＳＩＮＲに基づいて、ランク・レート予測部３３はキャパシティＣａｐとＳＩＮＲＥｆｆとを算出する。続いて、ステップＳ１１０９に進む。なお、キャパシティは（４４）式、ＳＩＮＲＥｆｆは（４５）式より算出することができる。
【０１２７】
【数４４】
2009152839.tif 000046
【０１２８】
【数４５】
2009152839.tif 000047
【０１２９】
（ステップＳ１１０９）ランク・レート予測部３３はｉｎｉｔが１であるか否か判断する。ランクＭ＝ＭＴの場合はｉｎｉｔ＞１であるため、ステップＳ１１１１に進む。
【０１３０】
（ステップＳ１１１１）ストリームと物理アンテナ１６とは１対１の関係であるため、各ストリームはＳＩＮＲＥｆｆの平均値で受信することが出来ると考えられる。よって、ステップＳ１１０８で算出したキャパシティＣａｐとＳＩＮＲＥｆｆに基づいて、ランク・レート予測部３３は、ＳＩＮＲＥｆｆの平均値ＳＩＮＲＥｆｆａｖｇを算出する。また、ランク・レート予測部３３は、算出したＳＩＮＲＥｆｆａｖｇに基づいて、キャパシティＣａｐａｖｇを算出する。なお、ＳＩＮＲＥｆｆａｖｇは（４６）式、Ｃａｐａｖｇは（４７）式より算出することができる。
【０１３１】
【数４６】
2009152839.tif 000048
【０１３２】
【数４７】
2009152839.tif 000049
【０１３３】
（ステップＳ１１１２）ランク・レート予測部３３はランクＭ＝ＭＴの場合のＦｉｎａｌ＿Ｃａｐ（ランク＝ＭＴ）およびＦｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝ＭＴ）を求める。その後、ステップＳ１１１３に進む。なお、（４８）式に示すとおり、Ｆｉｎａｌ＿Ｃａｐ（ランク＝ＭＴ）はＣａｐａｖｇのＭＴ倍の値となる。また、（４９）式に示すとおり、Ｆｉｎａｌ＿ＳＩＮＲＥｆｆ（ランク＝ＭＴ）は、ＳＩＮＲＥｆｆａｖｇとなる。
【０１３４】
【数４８】
2009152839.tif 000050
【０１３５】
【数４９】
2009152839.tif 000051
【０１３６】
（ステップＳ１１１３）ランク・レート予測部３３はｉｎｉｔがＭＴ以上か否か判断する。ランクＭ＝ＭＴの場合、ｉｎｉｔ＝ＭＴであるためステップＳ１１１５に進む。
【０１３７】
ランクＭ＝ＭＴの場合における、ステップＳ１１０４〜ステップＳ１１１４の処理の手順は以上の通りである。
【０１３８】
以下、ランクＭが１≦Ｍ≦ＭＴの各値の場合において算出したＦｉｎａｌ＿Ｃａｐに基づいて、ランク・レート予測部３３が予測する予測ランク数および予測レートの決定方法を説明する。
【０１３９】
（ステップＳ１１１５）ランクＭが１≦Ｍ≦ＭＴの各値の場合においてステップＳ１１１２で算出したＦｉｎａｌ＿Ｃａｐを比較し、ランク・レート予測部３３はＦｉｎａｌ＿Ｃａｐの値が最大となるランク数を算出する。このランク数を予測ランク数とする。その後、ステップＳ１１１６に進む。
【０１４０】
（ステップＳ１１１６）ランク・レート予測部３３は、ステップＳ１１１５で算出した予測ランク数のときのＦｉｎａｌ＿ＳＩＮＲを決定し、決定したＦｉｎａｌ＿ＳＩＮＲに基づいてレート（ＭＣＳ）を決定する。これを予測レートする。その後、ステップＳ１１１７に進む。
【０１４１】
（ステップＳ１１１７）ステップＳ１１１５およびステップＳ１１１６で決定した予測ランク数および予測レートと、ステップＳ１１０４で決定した予測ランク数の場合でのストリームと物理アンテナ１６とのマッピング情報を送信局に送信する。その後、処理を終了する。
【０１４２】
なお、送信局となる無線通信装置は、受信局となる無線通信装置より受信した予測ランクおよびマッピング情報に基づいてストリームを送信する物理アンテナ１６を決定し、受信局となる無線通信装置に対してデータの送信を行う。
【０１４３】
上述したとおり、ランク数が１からＭＴ（ＭＴは送信局の物理アンテナ数）の場合を想定し、各ランク数の場合におけるストリームと物理アンテナとのマッピングを決定し、決定したマッピングの場合のキャパシティをランク数が１からＭＴの場合毎に算出する。その後、ランク数毎に算出したキャパシティを比較し、一番良いキャパシティの場合のランク数を予測ランク数と決定し、決定した予測ランク数に基づいたレートを予測レートとする。このことにより、オープンループＳＣＷ−ＭＩＭＯにおいて、より正確なランク数とレートの予測を行うことが可能となる。
【０１４４】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【図面の簡単な説明】
【０１４５】
【図１】本発明の一実施形態における無線通信装置の受信部分に関する構成を示した図である。
【図２】本実施形態におけるランクおよびレートの測定手順を示すフローチャートである。
【図３】ＳＣＷ−ＭＩＭＯを使用する送信局の構成例を示した図である。
【図４】ＳＣＷ−ＭＩＭＯを使用する受信局の構成例を示した図である。
【図５】６４ＱＡＭの限界性能カーブＣ６４ＱＡＭを示した図である。
【符号の説明】
【０１４６】
１１・・・ターボエンコーダ部、１２・・・ＱＡＭ部、１３・・・シリアル／パラレル変換部、１４・・・マッピング部、１５・・・ＯＦＤＭ変調部、１６，２１・・・物理アンテナ、１７・・・受信部、１８・・・レート予測部、２１・・・物理アンテナ、２２・・・ＯＦＤＭ復調部、２３・・・データ分離部、２４・・・パラレル／シリアル変換部、２５・・・対数尤度比算出部、２６・・・ターボデコーダ部、２７・・・ランク・ＣＱＩ予測部、３１・・・ＳＮＲ計算部、３２・・・マッピング部、３３・・・ランク・レート予測部
【図１】
2009152839.tif 000052
【図２】
2009152839.tif 000053
【図３】
2009152839.tif 000054
【図４】
2009152839.tif 000055
【図５】
2009152839.tif 000056