(57)【要約】【課題】マルチキャリア、ＭＩＭＯ方式およびビームフォーミング処理を組合せたシステムに対する試験装置を、小さな回路規模、少ない消費電力で実現できるようにする。【解決手段】レイヤ周波数領域信号生成部３１が出力するＲレイヤ分の各キャリアの変調信号に対し、窓関数演算部３２が窓関数の周波数特性の畳み込み演算を行なう。一方、フェージング設定部５１が出力する各パスの伝搬路特性に対して、ビームフォーミング等価演算部５２がビームフォーミング処理と等価の演算処理を行い、その演算結果に対して、フーリエ変換部５３がフーリエ変換を行う。演算部５４は両演算結果Ｈt 、Ｆsymの乗算を行い、各受信アンテナで受信される信号のスペクトラム情報を求め、これを時間領域信号生成部３３によるフーリエ逆変換処理により時間領域の信号に変換し、シフト加算部３４によって窓関数の長さ分ずらして加算することで、各受信アンテナの受信信号を生成する。【選択図】図４
2018101841.tif 000002
【特許請求の範囲】
【請求項１】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で求めたＮ×Ｍ×Ｕパス分の伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２）と、
前記ビームフォーミング等価演算部によって得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３）と、
前記フーリエ変換部によって得られた周波数領域における伝搬路特性と前記窓関数演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えたことを特徴とするＭＩＭＯ方式システムの試験装置。
【請求項２】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３′）と、
前記フーリエ変換部で求めた周波数領域における伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２′）と、
前記ビームフォーミング等価演算部の演算結果と前記窓関数演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４′）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えたことを特徴とするＭＩＭＯ方式システムの試験装置。
【請求項３】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記窓関数演算部の演算結果に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２″）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３′）と、
前記フーリエ変換部によって求められた周波数領域における伝搬路特性とビームフォーミング等価演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４″）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えたことを特徴とするＭＩＭＯ方式システムの試験装置。
【請求項４】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについて求めた伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記ビームフォーミング処理と等価の演算処理で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記周波数領域における伝搬路特性と前記窓関数の周波数特性の畳み込み演算の結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含むことを特徴するＭＩＭＯ方式システムの試験方法。
【請求項５】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについて求めた伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記フーリエ変換によって得られた伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記ビームフォーミング処理と等価の演算処理で得られた伝搬路特性と前記窓関数の周波数特性の畳み込み演算の結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含むことを特徴するＭＩＭＯ方式システムの試験方法。
【請求項６】
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記窓関数の畳み込み演算の演算結果に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについての伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記周波数領域における伝搬路特性と、前記ビームフォーミング処理と等価の演算処理の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含むことを特徴するＭＩＭＯ方式システムの試験方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、基地局から移動体端末へのダウンリンク信号を、基地局側アンテナ数Ｎ、端末側アンテナ数Ｍで伝達するＭＩＭＯ（Multi Input Multi Output）方式に対応した端末あるいはその端末に内蔵される回路基板や集積回路等を試験対象とし、送受信アンテナ間に想定されるＮ×Ｍチャネルの伝搬路に対してフェージング処理を行なう機能を有する試験装置の回路規模を小さくするための技術に関する。
【背景技術】
【０００２】
ＭＩＭＯ方式は、図９に示すように、端末側へのダウンリンク信号Ｓtx1 〜ＳtxN を、Ｎ本(この例ではＮ＝４とする)の基地局側アンテナ（以下、送信アンテナ）Ａtx1 〜ＡtxN から送信し、Ｍ本(この例ではＭ＝２とする)の端末側アンテナ（以下、受信アンテナ）Ａrx1 〜ＡrxM で受信させる。
【０００３】
したがって、各送信アンテナと各受信アンテナの間にＮ×Ｍの伝搬路（チャネル）が想定され、また各チャネルについて異なる複数Ｕ（例えばＵ＝４）のパスが想定される。パスを含めた各チャネルの伝搬特性をＨ(1,1,1〜U)〜Ｈ(N,M,1〜U)とすると、ＭＩＭＯ方式に対応した移動体端末やそれに用いる回路等を試験する場合には、ダウンリンク信号に対し、各チャネルの伝搬特性およびパスについての損失、遅延、ドップラシフト等の特性を加味した演算処理を行い、最終的にＭ本の受信アンテナＡrx1 〜ＡrxM から出力される受信信号Ｓrx1 〜ＳrxM を生成して試験対象１に与える必要がある。
【０００４】
一方、近年では変調方式として、ＯＦＤＭ（Orthogonal Frequency Division Multiplexing）、ＵＦＭＣ（Universal Filtered Multicarrier）、ＧＦＤＭ（Generalized Frequency Division Multiplexing）、ＦＢＭＣ（Filtered Bank Multi-Carrier）等のマルチキャリア変調方式による高速な信号伝達が実現されており、このマルチキャリア変調方式とＭＩＭＯ方式との組合せにより、さらに高速な情報通信が可能なＭＩＭＯ方式システムが実現され、そのシステムを試験する装置が必要となる。
【０００５】
また、次の世代（第５世代）の通信方式では、より高い周波数帯を用いることが提案されている。このように通信に使用する周波数帯が高くなると、アンテナ個々のサイズを小さく構成できるため、多数のアンテナ素子を縦横に配列したアレーアンテナ構造を採用して、それらのアンテナ素子に与えるダウンリンク信号の位相制御により、通信対象の移動体端末の存在する方向に効率的に電波を放射させる、所謂ビームフォーミングが可能となる。したがって、このような次世代の移動体端末を試験対象とする試験装置においては、アレー化された多数のアンテナについてのビームフォーミングの演算処理が必要となる。
【０００６】
図１０は、マルチキャリア変調方式、ＭＩＭＯ方式およびアレーアンテナによるビームフォーミング処理を組合せたシステムを試験するための試験装置の構成例を示している。
【０００７】
この試験装置１０は、マルチキャリア変調方式のうち、複数Ｋのサブキャリアを用いて端末との通信を行なうＯＦＤＭに対応したものであり、レイヤ周波数領域信号生成部１１において、試験対象に伝達しようとするＲ系列の伝送データ（レイヤあるいストリームと呼ばれる）に対してそれぞれ複数Ｋのサブキャリア毎の変調信号（コンスタレーションデータ）Ｓsym(1,1)〜Ｓsym(1,K)、Ｓsym(2,1)〜Ｓsym(2,K)、……、Ｓsym(R,1)〜Ｓsym(R,K)を生成出力する。この変調信号Ｓsym は、ＯＦＤＭシンボル毎に、周波数軸上にＫ個のコンスタレーションデータを並べたデータをＲ系列分含む周波数領域の信号である。
【０００８】
これらの変調信号Ｓsym はビームフォーミング処理部１２に入力され、Ｎ本の送信アンテナから出射される電波のビーム特性が所望特性となるように演算処理され、１送信アンテナ当り複数Ｋのサブキャリア毎のビームフォーミング処理信号Ｓbf(1,1)〜Ｓbf(1,K)、Ｓbf(2,1)〜Ｓbf(2,K)、……、Ｓbf(N,1)〜Ｓbf(N,K)に変換される。なお、以下の説明では、図面も含めて、ｊ個一組の信号Ｓx(i,1)〜Ｓx(i,j)をＳx(i,1〜j)と略記する場合がある。
【０００９】
これらのビームフォーミング処理信号Ｓbf は、Ｎ組の時間領域信号生成部１３(1)〜１３(N)に入力される。各時間領域信号生成部１３(i)（ｉ＝１〜Ｎ）は、Ｋ個１組のビームフォーミング処理信号Ｓbf(i,1〜K)に対して、フーリエ逆変換（ＩＦＦＴ）処理、サイクリックプレフィクス（ＣＰ）付加処理、帯域制限処理等を行い、ＯＦＤＭ方式で規定された時間軸上の信号に変換する。
【００１０】
これによって、各時間領域信号生成部１３(1)〜１３(N)からは、Ｎ本の送信アンテナＡtx1 〜ＡtxN に与えるための送信信号(ダウンリンク信号)Ｓtx1 〜ＳtxN が出力されることになる。
【００１１】
これらの送信信号Ｓtx1 〜ＳtxN は、Ｎ×Ｍチャネルの伝搬路の特性を模擬する伝搬路シミュレータ２０に入力される。
【００１２】
伝搬路シミュレータ２０は、Ｎ本の送信アンテナとＭ本の受信アンテナの間に形成されるＮ×Ｍのチャネルと、それら各チャネルについてそれぞれ複数Ｕのパスとを想定し、これらＮ×Ｍ×Ｕの各パスに所望の遅延とフェージングを付加し、Ｍ本の受信アンテナがそれぞれ受ける受信信号を仮想的に生成するものである。
【００１３】
この伝搬路シミュレータ２０は、無線通信における受信レベル変動の分布を表すレイリーフェージングを与えるものであり、Ｎ系列の送信信号Ｓtx1 〜ＳtxN それぞれに設定される複数Ｕのパスに対し、所定の遅延を与えて出力する遅延設定部２１と、ドップラシフト、ＭＩＭＯ相関の情報が付与されたレイリー分布の伝搬路の特性を求めるフェージング設定部２２と、遅延設定部２１から出力される全パス分の遅延処理信号Ｓtx(1,1,1〜U)、Ｓtx(2,1,1〜U)、……、Ｓtx(N,M,1〜U)と、フェージング設定部２２で得られた伝搬特性Ｈ(1,1,1〜U）、Ｈ(2,1,1〜U）、……、Ｈ(N,M,1〜U）とを用いた積和演算（行列の乗算）により、Ｎ×Ｍ×Ｕの仮想的な伝搬路を経由してＭ本の受信アンテナで受信される受信信号Ｓrx1 〜ＳrxM を生成する演算部２３とを有している。
【００１４】
ここで、遅延設定部２１は、例えばメモリによる１サンプル単位の遅延とリサンプルフィルタによる１サンプル以下の遅延の組み合わせで、各パスに所望の遅延を付与している。
【００１５】
また、演算部２３の演算処理は、例えば、
Ｓrx1＝ΣＨ(1,1,i)・Ｓtx(1,1,i)＋ΣＨ(2,1,i)・Ｓtx(2,1,i)＋……
＋ΣＨ(N,1,i)・Ｓtx(N,1,i)
Ｓrx2＝ΣＨ(1,2,i)・Ｓtx(1,2,i)＋ΣＨ(2,2,i)・Ｓtx(2,2,i)＋……
＋ΣＨ(N,2,i)・Ｓtx(N,2,i)
……
ＳrxM＝ΣＨ(1,M,i)・Ｓtx(1,M,i)＋ΣＨ(2,M,i)・Ｓtx(2,M,i)＋……
＋ΣＨ(N,M,i)・Ｓtx(N,M,i)
となる。ただし、記号Σは、ｉ＝１〜Ｕまでの総和を表す。
【００１６】
このようにして得られた受信信号Ｓrx1 〜ＳrxM を試験対象１に与えることで、試験装置側で設定した送受信アンテナ間の伝搬路の状態に対する試験対象１の動作を試験することができる。
【００１７】
なお、伝搬路シミュレータは含まれていないが、上記したようにマルチキャリア変調方式とＭＩＭＯ方式とを組合せたシステムを試験するため試験装置は、例えば次の特許文献１に開示されている。
【先行技術文献】
【特許文献】
【００１８】
【特許文献１】特開２０１４−９３７５８号公報
【発明の概要】
【発明が解決しようとする課題】
【００１９】
上記構成の試験装置のように、ビームフォーミング処理を行なうシステムではアレー化された送信アンテナの数Ｎが例えば１２８のように非常に大きくなるため、それに伴って、フーリエ逆変換処理をＮ系列分並列的に行なう時間領域信号生成部１３が１２８組必要となり、回路規模が非常に大きくなってしまう。
【００２０】
また、伝搬路シミュレータ２０の遅延設定部２１は、前記したように、メモリとリサンプルフィルタの組合せにより任意の遅延を付与するハードウエア構成を要するため、上記のように１２８系列の信号に対してそれぞれ設定した複数Ｕのパスに任意の遅延を付与するためには、やはりその回路規模が非常に大きくなり、装置が大型化し、製造コストおよび消費電力が大きくなってしまう。
【００２１】
本発明は、上記問題を解決し、マルチキャリア変調方式、ＭＩＭＯ方式およびビームフォーミング処理を組合せたシステムで送信アンテナ数が多い場合でも、小さな回路規模、少ない消費電力で実現できる試験装置および試験方法を提供することを目的としている。
【課題を解決するための手段】
【００２２】
前記目的を達成するために、本発明の請求項１のＭＩＭＯ方式システムの試験装置は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で求めたＮ×Ｍ×Ｕパス分の伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２）と、
前記ビームフォーミング等価演算部によって得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３）と、
前記フーリエ変換部によって得られた周波数領域における伝搬路特性と前記窓関数演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えている。
【００２３】
また、本発明の請求項２のＭＩＭＯ方式システムの試験装置は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３′）と、
前記フーリエ変換部で求めた周波数領域における伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２′）と、
前記ビームフォーミング等価演算部の演算結果と前記窓関数演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４′）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えている。
【００２４】
また、本発明の請求項３のＭＩＭＯ方式システムの試験装置は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記レイヤ周波数領域信号生成部が出力するＲ×Ｋ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう窓関数演算部（３２）と、
前記窓関数演算部の演算結果に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なうビームフォーミング等価演算部（５２″）と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求めるフェージング設定部（５１）と、
前記フェージング設定部で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求めるフーリエ変換部（５３′）と、
前記フーリエ変換部によって求められた周波数領域における伝搬路特性とビームフォーミング等価演算部の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める演算部（５４″）と、
前記演算部の演算結果に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（３３）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成するシフト加算部（３４）とを備えている。
【００２５】
また、本発明の請求項４のＭＩＭＯ方式システムの試験方法は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについて求めた伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記ビームフォーミング処理と等価の演算処理で得られた全てのパスの伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記周波数領域における伝搬路特性と前記窓関数の周波数特性の畳み込み演算の結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含んでいる。
【００２６】
また、本発明の請求項５のＭＩＭＯ方式システムの試験方法は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについて求めた伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記フーリエ変換によって得られた伝搬路特性に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記ビームフォーミング処理と等価の演算処理で得られた伝搬路特性と前記窓関数の周波数特性の畳み込み演算の結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含んでいる。
【００２７】
また、本発明の請求項６のＭＩＭＯ方式システムの試験方法は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式および前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、該伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号をＲ×Ｋ系列分生成する段階と、
前記Ｒ×Ｎ系列分の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行なう段階と、
前記窓関数の畳み込み演算の演算結果に対し、前記アンテナ数Ｎの送信アンテナによる放射ビーム特性を所望特性にするためのビームフォーミング処理と等価の演算処理を行なう段階と、
前記送信アンテナと前記受信アンテナの間に想定される全てのパスについての伝搬路特性を求める段階と、
前記全てのパスについての伝搬路特性に対し、パス毎の遅延を加味したフーリエ変換を行ない、周波数領域における伝搬路特性を求める段階と、
前記周波数領域における伝搬路特性と、前記ビームフォーミング処理と等価の演算処理の演算結果との乗算により、前記各受信アンテナでそれぞれ受信される信号のスペクトラム情報を求める段階と、
前記スペクトラム情報に対してフーリエ逆変換処理を行ない、前記各受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記各受信アンテナでそれぞれ受信される連続した受信信号を生成する段階とを含んでいる。
【発明の効果】
【００２８】
このように、本発明の請求項１、４では、Ｒレイヤのデータ信号列についてそれぞれキャリア毎の変調信号を生成し、この変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、窓関数の周波数特性の畳み込み演算を行なうとともに、全てのパスについての伝搬路特性に対して、複数の送信アンテナによる放射ビーム特性を決めるビームフォーミング処理と同等の演算処理を行い、その結果に対して、パス毎の遅延を加味したフーリエ変換を行ない、ビームフォーミングを加味した周波数領域における伝搬路特性を求め、この周波数領域における伝搬路特性と窓関数の周波数特性の畳み込み演算の結果との乗算により、各受信アンテナでそれぞれ受信される信号の周波数領域の情報（スペクトラム情報）を求め、これをフーリエ逆変換処理して時間領域の信号を生成し、これを窓関数の長さ分ずらして加算することで、各受信アンテナでそれぞれ受信される連続した受信信号を生成している。
【００２９】
このように、本発明では、伝搬路の特性と変調信号との乗算演算を周波数領域において行ない、その演算結果から時間領域信号を生成しているので、従来方式のように、送信アンテナ毎の周波数領域の信号をフーリエ逆変換して時間領域の信号に変換してから伝搬路特性を付与する場合に比べて、フーリエ逆変換を行なう回路および伝搬路特性を生成する回路の規模を格段に小さくすることができる。
【００３０】
例えば、Ｒ＝２、Ｎ＝１２８、Ｍ＝８、キャリア数Ｋの場合、従来方式では、時間領域信号生成のために、Ｋ個一組の信号を１２８（＝Ｎ）組分並列的にフーリエ逆変換する必要があるが、本発明では、最小の場合として、Ｋ個一組の信号を８（＝Ｍ）組分並列的にフーリエ逆変換すればよく、回路規模をＭ／Ｎに縮小できる。
【００３１】
ただし、本発明では、全パスの伝搬路特性を周波数領域に変換するためのフーリエ変換処理が必要となるが、このフーリエ変換処理では、時間領域における各パスの遅延量は、周波数領域では各パスの周波数成分の回転速度に対応するため、従来時間領域で行なっていたメモリとリサンプルフィルタの組合せにより各パスに遅延を付与するハードウエアは、フーリエ変換における回転処理に置き換えられることになり、両者のハードウエアの規模を比較して本発明の方が格段に有利となる。
【００３２】
しかも、ビームフォーミング等価演算部の演算処理をフェージング設定部の出力に対して行ない、その結果に対してフーリエ変換を行なうようにしているので、フーリエ変換処理はＲ×Ｍ×Ｕ系列分行なえばよく、その結果を用いた演算部の演算も、最小でＲ×Ｋ系列の入力に対しておこなえばよいので、フーリエ変換処理の増加よりも演算部の演算処理の減少効果の方が格段に大きく、これによる回路規模のさらなる縮小が可能である。
【００３３】
また、本発明の請求項２、４では、ビームフォーミング処理と等価な演算処理を、伝搬路の特性に対するフーリエ変換によって得られた結果に対して行なう点で、請求項１、３の発明と異なるが、伝搬路の特性と変調信号との乗算演算を周波数領域において行なっている点では、請求項１、３の発明と同じであり、時間領域信号生成のための演算処理を行なう回路規模を従来装置に比べてＭ／Ｎに縮小でき、Ｍに対してＮが格段に大きいシステムを小規模に実現できる。
【００３４】
また、本発明の請求項３、６では、ビームフォーミング処理と等価な演算処理を、窓関数の畳み込み演算の結果に対して行なう点で、請求項１、２、４、５の発明と異なるが、伝搬路の特性と変調信号との乗算演算を周波数領域において行なっている点では、請求項１、２、４、５の発明と同じであり、時間領域信号生成のための演算処理を行なう回路規模を従来装置に比べてＭ／Ｎに縮小でき、Ｍに対してＮが格段に大きいシステムを小規模に実現できる。
【図面の簡単な説明】
【００３５】
【図１】本発明の原理を説明するためのタイミング図
【図２】時間領域の窓関数の一例を示す図
【図３】時間領域の窓関数の別の例を示す図
【図４】本発明の実施形態の構成を示す図
【図５】本発明の実施形態の要部の構成図
【図６】本発明の実施形態の要部の構成図
【図７】本発明の他の実施形態の構成を示す図
【図８】本発明の他の実施形態の構成を示す図
【図９】マルチパスＭＩＭＯの伝搬路の一例を示す図
【図１０】従来装置の構成図
【発明を実施するための形態】
【００３６】
以下、図面に基づいて本発明の実施の形態を説明するが、具体的な構成を説明する前に、本発明の試験装置の原理について説明する。
【００３７】
本発明は、前記したＯＦＤＭ、ＵＦＭＣ、ＧＦＤＭ、ＦＢＭＣなどのマルチキャリア変調方式で、Ｎ×Ｍ  ＭＩＭＯ（Ｎ＞Ｍ）を実施する場合の伝搬路シミュレータとして適用できるものであり、３Ｄ−ＭＩＭＯ／Ｍａｓｓｉｖｅ−ＭＩＭＯのように、送信アンテナ数が受信アンテナ数に比べて非常に多い場合に特に有効である。以降では、変調方式としては、主にＯＦＤＭを念頭に置いて説明する。
【００３８】
本発明では、次式（１）のように、ＭＩＭＯ伝搬路の特性の時間変化が無視できる程度の時間スパン毎（Ｔｃ毎）に、その時間スパン内ではＭＩＭＯ伝搬路の特性が一定であるとして、周波数領域でＭＩＭＯ伝搬路処理を施すというものである。
【００３９】
Ｔｃ≪１／ｆｄ    （ｆｄ：ドップラ周波数）        ……（１）
【００４０】
例えば、ＯＦＤＭの場合には、次式（２）のように、１ＯＦＤＭシンボル長Ｔsym（＝有効データ長＋サイクリックプレフィクス長）をＰ分割した長さＴｃが式（１）を満たすようにする（Ｐ分割は必ずしも等分割でなくても良い）。
【００４１】
Ｔｃ＝Ｔsym／Ｐ    （Ｐ＝１，２，３，……）      ……（２）
【００４２】
図１は時間領域におけるＰ＝２の例を示すものであり、図１の（ａ）に示すＯＦＤＭの信号列に対して、図１の（ｂ１）〜（ｂ４）のように、１シンボル長Ｔsym をＴｃ＝Ｔsym／２で２分割するような矩形窓関数を乗算して切り取った波形に対し、マルチパス伝搬処理、フィルタ処理などを施した波形を図１の（ｃ１）〜（ｃ４）のようにＴｃ分ずれた状態で得て、加算処理することで最終的な送信信号を得る。
【００４３】
実際に用いるローカライゼーション（信号切出し）用の窓関数の長さＴｃ′は、対応する周波数特性の広がりを抑えるために端部を丸めることによって、Ｔｃより若干大きくしてもよく、この窓関数のタイミングを順次Ｔｃずつずらしながら乗算した波形のそれぞれに、ＭＩＭＯチャネルのマルチパス伝搬路処理、フィルタ処理などを施す。また、分割された一つ一つの波形の時間長は、マルチパスの遅延時間分およびフィルタ処理による広がり分ＴｄだけＴｃ′よりも長くなる。それらをＴｃずつずらしながら加算した波形を処理結果とする考え方である。処理結果の波形は、Ｎ×Ｍ  ＭＩＭＯ伝搬路の場合には、Ｍ系統分計算される。
【００４４】
図２は、ローカライゼーション用の窓関数（区間長：Ｔｃ′）の詳細を示すもので、例えば、１ＯＦＤＭシンボルを複数個に分割する場合などに利用可能である。ナイキスト基準を満たすような特性になっており、この窓関数をＴｃずつずらした区間が連続的につながるような特性である。この時間領域における窓関数のロールオフが大きいほど、周波数領域での広がりが抑えられて、後述する窓関数演算部３２におけるフィルタのタップ数を抑えることができる。
【００４５】
また、変調方式がＦＢＭＣで、１シンボル情報が広がった時間（Ｖ・Ｔｃ）（Ｖはオーバラッピングファクタ）において、ＭＩＭＯ伝搬路の特性が一定である（Ｖ・Ｔｃ≪１／ｆｄ）と見なすことができる場合には、図３に示しているように、そのＦＢＭＣ用のローカライゼーション用窓関数そのものを使うこともできる。
【００４６】
上記処理は時間軸上で想定したものであるが、本発明は、最後の加算処理以外を周波数領域で等価な処理を実施するとともに、従来ではレイヤ信号に対して行なっているビームフォーミング処理と等価な処理を他の信号に対して行なうことにより、試験装置全体の回路規模を縮小している。
【００４７】
次に、本発明を適用した試験装置の実施形態を説明する。
図４は、本発明の実施形態の試験装置３０の構成を示している。
【００４８】
試験装置３０は、１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式で且つ送信アンテナ数Ｎ、受信アンテナ数ＭのＭＩＭＯ方式およびアンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、送信アンテナと受信アンテナの間にＮ×Ｍのチャネルと各チャネルにそれぞれ複数Ｕのパスを有する擬似的な伝搬路を想定し、その伝搬路を経由してＭ個の受信アンテナで受信される受信信号を生成して試験対象に与えるＭＩＭＯ方式システムの試験装置である。なお、以下の説明は、マルチキャリア変調方式がＯＦＤＭの場合とし、ＯＦＤＭでは端末との通信に用いる複数のキャリアを「サブキャリア」と呼んでいるので、以下の説明でもこの「サブキャリア」と記す。
【００４９】
この試験装置３０は、レイヤ周波数領域信号生成部３１、窓関数演算部３２、時間領域信号生成部３３、シフト加算部３４、伝搬路シミュレータ５０を有している。
【００５０】
レイヤ周波数領域信号生成部３１は、試験対象に伝達しようとするＲ系列の伝送データ（レイヤあるいストリームと呼ばれる）に対してそれぞれ複数Ｋのサブキャリア毎の変調信号（コンスタレーションデータ）Ｓsym(1,1)〜Ｓsym(1,K)、Ｓsym(2,1)〜Ｓsym(2,K)、……、Ｓsym(R,1)〜Ｓsym(R,K)を生成出力する。この変調信号は、ＯＦＤＭシンボル毎に、周波数軸上にＫ個のコンスタレーションデータを並べたデータをＲ系列分含む周波数領域の信号である。なお、レイヤ数Ｒは、原理上、試験対象の受信アンテナ数Ｍ以下の値である。
【００５１】
ここで、そのコンスタレーションデータを、次のようにＳsym,r,kという記号の複素数で表すことにする。
【００５２】
Ｓsym,r,k      ……（３）
sym ：ＯＦＤＭシンボル番号
ｒ＝｛１，２，３，……，Ｒ｝：送信レイヤ番号インデックス
ｋ＝｛１，２，３，……，Ｋ｝：サブキャリア番号インデックス
【００５３】
なお、図４では、信号Ｓsym,r,kのインデックスが分かりやすいように、Ｓsym（r,k）の形で表す（他の信号についても同様）。
【００５４】
また、コンスタレーションが配置される間隔（サブキャリア間隔）をｆscとする。ｆscはＯＦＤＭシンボル長Ｔsym（＝有効データ長＋サイクリックプレフィクス長）と次の関係にある。
有効データ長＝１／ｆsc    ……（４）
Ｔsym＝（１／ｆsc）＋サイクリックプレフィクス長    ……（５）
【００５５】
また、窓関数演算部３２は、レイヤ周波数領域信号生成部３１から出力されるＫ×Ｒ系列分の変調信号に対し、時間領域でのローカライゼーション用の窓関数の周波数特性の畳み込み演算を行ない、時間領域でローカライゼーション用の窓関数の掛け算に相当する結果を得る（ただし、窓関数の時間長は、伝搬路特性の変化が無視できる程度の時間長とする必要がある。）
【００５６】
より具体的に説明すると、窓関数演算部３２は次の処理を実施する。
時間領域で区間長Ｔｃ′のローカライゼーション用の窓関数ｆwτ（τは時間軸方向のインデックス）を乗算するのと等価な処理として、ローカライゼーション用の窓関数のフーリエ変換（Ｃoep,i）［ｉは周波数方向の係数インデックス、ｐは１ＯＦＤＭシンボル中の窓関数の番号］を周波数領域で畳み込み処理する。
【００５７】
ここでは、長さＴsym の１ＯＦＤＭシンボルをＰ分割したｐ番目（ｐ＝１，２，３，……,Ｐ）の区間長Ｔｃ′のローカライゼーション用の窓関数の乗算計算について以下のように数式化する。ただし、（Ｔｃ′＋Ｔｄ）と１／ｆsc（レイヤ周波数領域信号生成部３１の出力のサブキャリア間隔ｆscでＩＦＦＴを実施するときの時間領域のサイクル）の大小関係に応じて、処理後の周波数領域上でのサンプリング間隔が決まる。
【００５８】
（ａ）  （Ｔｃ′＋Ｔｄ）＞１／ｆscの場合
周波数領域でのサンプリング間隔を細かくすることで、時間長（Ｔｃ′＋Ｔｄ）の波形が時間領域でのエイリアシング（オーバーラップ）を生じないようにする。つまり、次式（６）、（７）に示しているように、補間係数をＤscとし、（Ｔｃ′＋Ｔｄ）＜１／（ｆsc／Ｄsc）を満たすように、畳み込み処理結果の周波数領域でのサンプリング間隔が１／Ｄsc倍となるように補間処理をする。（元のサブキャリア同士の間にＤsc−１個のゼロを置いてから畳み込みによるフィルタ処理を実施することに相当する処理である。）
【００５９】
【数１】
2018101841.tif 000003
【００６０】
ここで、式（６）の記号Σは、ｉ＝−TapNum／２〜TapNum／２までの総和を表す。また、window(ｉ)はタップ長がＤsc・（TapNum ＋１）の窓関数であるとする。ＤＦＴ(ｆwτ)は、時間スパン１／（ｆsc／Ｄsc）に渡る離散的フーリエ変換である。ｆwτはその中心が時間０に位置しており、Ｔｃ・（１／２＋ｐ−１）だけ波形を遅らせるとｐに対応する位置に移動するような波形であるとする。
【００６１】
また、式（７）のｅＺは、時間領域でＴｃ・（１／２＋ｐ−１）だけ波形を遅らせるのに相当する周波数領域上での回転を与える項である。
【００６２】
また、式（６）において、ｋ’は補間処理後の周波数インデックスを表しているが、
ｋ’＝｛Ｄsc・（０−<Ｋ／２>），Ｄsc・（１−<Ｋ／２>），…，Ｄsc・（Ｋ−１−
<Ｋ／２>）｝
の位置にそれぞれｋ＝｛１，２，３，……Ｋ｝の変調波Ｓsym,r,kが対応する関係にある。なお、式（６）の括弧記号<Ａ>は、Ａを越えない最大の整数を表す（以下、同様）。
【００６３】
また、式（６）の記号％は剰余演算子であり、ｇは、ｋ’をＤscで割った時の余りである。ただし、式（６）は、
−Ｄsc・（<Ｋ／２>＋TapNum／２）≦ｋ’≦Ｄsc・（<Ｋ／２>＋TapNum／２）
の範囲で計算する必要があり、ｉ＜０およびｉ＞Ｋの範囲でＳsym,r,i＝０とする。
【００６４】
（ｂ）  （Ｔｃ′＋Ｔｄ）＜１／ｆscの場合
次式（８）、（９）に示すように、畳み込み処理結果のサンプリング間隔は、周波数領域信号生成部３１の出力と変わらないように畳み込み処理する（補間無し）。
【００６５】
【数２】
2018101841.tif 000004
【００６６】
ここで、式（８）の記号Σは、ｉ＝−TapNum／２〜TapNum／２までの総和を表す。また、window(ｉ)はタップ長が（TapNum＋１）の窓関数であるとする。ＤＦＴ(ｆwτ)は、時間スパン１／ｆscに渡る離散的フーリエ変換である。ｆwτはその中心が時間０に位置しており、Ｔｃ・（１／２＋ｐ−１）だけ波形を遅らせるとｐに対応する位置に移動するような波形であるとする。
【００６７】
また、式（９）のｅＺ’は、時間領域でＴｃ・（１／２＋ｐ−１）だけ波形を遅らせるのに相当する周波数領域上での回転を与える項である。
【００６８】
また、式（８）において、ｋ’は周波数インデックスを表しているが、
ｋ’＝｛−<Ｋ／２>，１−<Ｋ／２>，…，Ｋ−１−<Ｋ／２>｝
の位置にそれぞれｋ＝｛１，２，３，……Ｋ｝の変調波Ｓsym,r,kが対応する関係にある。ただし、式（８）は、
−（<Ｋ／２>＋TapNum／２）≦ｋ’≦（<Ｋ／２>＋TapNum／２）
の範囲で計算する必要があり、ｉ＜０およびｉ＞Ｋの範囲でＳsym,r,i＝０とする。
【００６９】
伝搬路シミュレータ５０は、無線通信における受信レベル変動の分布を表すレイリーフェージングを与えるものであるが、本実施形態では、周波数領域でフェージングを付与する処理を行なっている。
【００７０】
この伝搬路シミュレータ５０は、フェージング設定部５１、ビームフォーミング等価演算部５２、フーリエ変換部５３および演算部５４を有している。
【００７１】
フェージング設定部５１は、伝搬路を形成する全てのパスについての振幅特性と位相特性にレイリー分布に従う変動を与えるとともに、試験対象の移動によるドップラシフト等を含む伝搬路特性Ａｔを求める。
【００７２】
このフェージング設定部５１の構成は任意であるが、例えば図５に示しているように、加法性白色ガウス雑音発生器（ＡＷＧＮ）５１ａで生成したレイリー分布の基になる白色ガウス雑音信号Ｇn(1,1,1〜U)〜Ｇn(N,M,1〜U)をドップラフィルタ５１ｂに入力してドップラスペクトルを付与し、その出力Ｄp(1,1,1〜U)〜Ｄp(N,M,1〜U)をＭＩＭＯ相関設定部５１ｃに入力してＭＩＭＯ相関を与え、その出力Ｍc(1,1,1〜U)〜Ｍc(N,M,1〜U)を補間器５１ｄに入力してその出力をパス毎の伝搬路特性Ａt(1,1,1〜U)〜Ａt(N,M,1〜U)として出力する。なお、補間器５１ｄは、窓関数演算部３２によるＴｃ間隔での信号切り出しとそれに乗算する伝搬路の特性とレートを合わせるための補間処理を行なうものである。この構成の他に、多数の素波の組合せにより伝搬路特性を求める方式も採用できる。
【００７３】
ビームフォーミング等価演算部５２は、フェージング設定部５１で求めたＮ×Ｍ×Ｕパス分の伝搬路特性Ａｔに対し、アンテナ数Ｎの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式と等価な演算処理を行ない、送信アンテナの放射ビーム特性が加味された全てのパスの伝搬路特性を求める。
【００７４】
このビームフォーミング処理は、入力される伝搬路特性Ａｔに対し、所望の放射ビーム特性が得られるように予め指定されたビームフォーミング行列Ｗを乗算するものである。
【００７５】
フーリエ変換部５３は、ビームフォーミング等価演算部５２によって得られた全てのパスの伝搬路特性ＡＷに遅延を与えてフーリエ変換することで、ＭＩＭＯ伝搬路の特性を表す周波数領域の信号Ｈｔを生成する。
【００７６】
具体的には、時間長Ｔｃ毎に、Ｒ×Ｍ個の、チャネルのインパルス応答×ビームフォーミング行列のフーリエ変換を実施する。各チャネルのパス数をＵ個とすると、時間ｔにおける、ｒ番目の送信レイヤとｍ番目の受信アンテナの間の伝搬路特性Ａt,ｒ.m.uのインパルス応答は次の式で表現されるとする。
ｈt,r,m＝ΣＡt,r,m,u・δ(ｔ−τｕ)    ……（１０）
ただし、記号Σは、ｕ＝１〜Ｕまでの総和を表す。
【００７７】
このフーリエ変換は次の式（１１）で表現できる。
【００７８】
【数３】
2018101841.tif 000005
【００７９】
ただし、記号Σはｕ＝１〜Ｕまでの総和を表す。ｋ’は周波数軸上のインデックスであり、
−Ｄsc・（<Ｋ／２>＋TapNum／２）≦ｋ’≦Ｄsc・（<Ｋ／２>＋TapNum／２）
の範囲をとる。また、Δｆはサブキャリア間隔を示している。
【００８０】
上記の式（１１）で示されるように、フーリエ変換部５３は、回転と累積加算の演算ブロックのみで構成され、回転の情報に遅延情報が含まれるため、従来装置のように、メモリとリサンプルフィルタの組合せにより各パスに遅延を付与する構成に比べて回路規模を格段に小さくできる。
【００８１】
演算部５４は、窓関数演算部３２の演算結果に対し、フーリエ変換部５３の出力を乗算することで、時間領域でのＭＩＭＯ伝搬路特性を周波数領域で与え、各受信アンテナで受信される信号の周波数領域の情報（スペクトラム情報）を求める。
【００８２】
この処理は、次式（１２）のように、周波数インデックスｋ’毎に、伝搬路行列を乗算することによって、ｍ番目の受信アンテナの周波数領域での受信信号を算出する。
【００８３】
Ｓｆsym,p,m,k'＝ΣＨt,r,m,k'・Ｆsym,p,r,k'    ……（１２）
ただし、記号Σは、ｒ＝１〜Ｒまでの総和を表し、Ｈt,r,m,k'の時間インデックスｔは、（sym,p）に対応した時間であるとする。前記同様に、周波数軸上のインデックスｋ'は、
−Ｄsc・（<Ｋ／２>＋TapNum／２）≦ｋ'≦Ｄsc・（<Ｋ／２>＋TapNum／２）
の範囲をとるものとする。
【００８４】
このようにして、伝搬路シミュレータ５０によって周波数領域におけるＭＩＭＯ伝搬路特性が与えられた信号Ｓｆsym,p,m,k'は、時間領域信号生成部３３に入力される。時間領域信号生成部３３は、図６に示すように、帯域制限フィルタ３３ａおよびフーリエ逆変換部３３ｂを有している。
【００８５】
帯域制限フィルタ３３ａは、次式（１３）に示すように、入力信号Ｓfsym,p,m,k'に対し、帯域制限フィルタの特性（BandFilk）の周波数領域における乗算を実施して、帯域制限を行なう。なお、この帯域制限処理は省略することも可能である。
Ｓbsym,p,m,k'＝Ｓfsym,p,m,k'・BandFilk'    ……（１３）
【００８６】
フーリエ逆変換部３３ｂは、次式（１４）のように、帯域制限された周波数領域の信号Ｓbsym,p,m,k'（あるいは伝搬路シミュレータ５０の出力信号Ｓfsym,p,m,k'）に対して、高速フーリエ逆変換ＩＦＦＴを行なうことで時間領域の信号Ｓtsym,p,m,τに変換する。
【００８７】
Ｓtsym,p,m,τ＝ＩＦＦＴ（Ｓbsym,p,m,k'）    ……（１４）
ただし、τ＝｛１，２，３，……，Ｎfft｝は時間のインデックスであるとする。ＮfftはＦＦＴポイント数とする。
【００８８】
さらに、ｋ'が、
Ｄsc・（<Ｋ／２>＋TapNum／２）＜ｋ'＜Ｎfft−Ｄsc・（<Ｋ／２>＋TapNum／２）
にある場合、Ｓbsym,p,m,k'＝０であり、Ｓbsym,p,m,k'は、Ｎfftを周期として周期的とする。即ち、整数ｉについて、Ｓbsym,p,m,k'＝Ｓbsym,p,m,(k'+i・Ｎfft)が成り立つものとする。
【００８９】
この時間領域に変換された信号Ｓtsym,p,m,τは、シフト加算部３４により、次式（１５）に示すように、前記窓関数の時間領域における長さ分ずらされて順次加算され、連続性を維持した受信信号が生成されることになる。つまり、上式（１４）の処理結果を、図１に示したように時間Ｔｃずつずらしながら加算することで１系列分の連続した受信信号を得る。これをＭ系列分並列的に行なうことで、Ｍ系列分の連続した受信信号を生成できる。
【００９０】
【数４】
2018101841.tif 000006
【００９１】
式（１５）で、ｆs は時間領域におけるサンプリング周波数であるとする。
【００９２】
なお、ビームフォーミング処理は、通常はレイヤ周波数領域信号生成部３１の出力に対して行なうものであるが、この実施形態のように、フェージング設定部５１の出力に対して行なっても、最終的な計算結果は同じになる。これについて以下に示す。
【００９３】
元のＲレイヤの変調信号の行列[Ｓs]、窓関数処理された信号の行列[Ｓs′]、ビームフォーミング行列[Ｗ]、送信アンテナと受信アンテナ間のチャネルの伝搬路の特性の行列［ａ０］、そのインパルス応答のフーリエ変換した行列［Ｈ０］、演算部５４の出力の行列［Ｓｆ］を、それぞれ次のように表す。ただし、記号Fourie{Ｂ｝はＢのフーリエ変換を表す。
【００９４】
【数５】
2018101841.tif 000007
【００９５】
そして、仮に、ビームフォーミング処理をＲレイヤの変調信号列に対して行なった場合、次の式（１６）が成り立つ。記号INTP{Ｃ｝は、Ｃに対する窓関数処理を表す。
【００９６】
［Ｓｆ］＝［Ｈ０］INTP｛[Ｗ][Ｓs]｝
＝［Ｈ０］[Ｗ]INTP｛[Ｓs]｝
＝［Ｈ０］[Ｗ][Ｓs′]          ……（１６）
【００９７】
上記式（１６）において、フーリエ変換後の行列［Ｈ０］を元の特性の行列［ａ０］で表して、展開すると以下のように表すことができる。
【００９８】
［Ｓｆ］＝Σ｛［ａ０］ｅＺ ｝[Ｗ][Ｓs′]          ……（１７）
＝Σ｛［ａ０］[Ｗ]ｅＺ ｝[Ｓs′]          ……（１８）
Ｚ＝−２π・τｕ・Δｆ・ｋ′
ただし、記号Σは、ｕ＝１〜Ｕまでの総和を表す。
【００９９】
上記式（１８）のΣ｛［ａ０］[Ｗ]ｅＺ ｝の項は、伝搬路のフェージング特性［ａ０］にビームフォーミング行列を乗じた結果（ＡＷ）をフーリエ変換したものであり、そのフーリエ変換の演算結果に窓関数処理で得られた信号[Ｓs′]を乗じているので、上記実施形態の構成による処理が実現されていることになる。つまり、実施形態のように、ビームフォーミング処理をフェージング設定部５１の出力に対して行なっても最終的な演算結果は変わらない。
【０１００】
このように、実施形態の試験装置３０は、Ｒレイヤのサブキャリア毎の変調信号に対し、時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、窓関数の周波数特性の畳み込み演算を行なうとともに、全てのパスについての伝搬路特性に対してビームフォーミング処理と等価な演算処理を行なってフーリエ変換を行ない、周波数領域における伝搬路特性を求め、この周波数領域における伝搬路特性と窓関数の周波数特性の畳み込み演算の結果との乗算により、各受信アンテナでそれぞれ受信される信号の周波数領域の情報（スペクトラム情報）を求め、これをフーリエ逆変換処理して時間領域の信号を生成し、これを窓関数の長さ分ずらして加算することで、各受信アンテナでそれぞれ受信される連続した受信信号を生成している。
【０１０１】
このため、従来方式のように、送信アンテナ毎の周波数領域の信号をフーリエ逆変換して時間領域の信号に変換してから伝搬路特性を付与する場合に比べて、フーリエ逆変換を行なう回路の規模を格段に小さくすることができる。
【０１０２】
例えば、Ｎ＝１２８、Ｍ＝８、レイヤ数Ｒ＝２、サブキャリア数Ｋの場合、従来方式では、時間領域信号生成のために、Ｋ個一組の信号を１２８（＝Ｎ）組分並列的にフーリエ逆変換する必要があるが、本実施形態では、Ｄsc・Ｋ個一組の信号を８（＝Ｍ）組分並列的にフーリエ逆変換すればよい。ここで、補間係数Ｄscが１（補間無しの場合）であれば乗算回数を、Ｍ／Ｎに縮小できる。に縮小できる。
【０１０３】
また、補間する場合には、乗算回数を、
（Ｍ・２Ｌ′・log２２Ｌ′）／（Ｎ・２Ｌ・log２２Ｌ）
に縮小することができ、Ｄsc・Ｍ＜Ｎであれば、従来回路より少ない乗算回数で実現できる。ただし、Ｌは、（２Ｌ−１）＜Ｋ≦２Ｌを満たす整数、Ｌ′は、（２Ｌ′−１）＜Ｄsc・Ｋ≦２Ｌ′を満たす整数である。
【０１０４】
しかも、ビームフォーミング処理を、フェージング設定部５１が出力する伝搬路の特性に対して行なう構成としたことで、窓関数演算部３２は、Ｒ×Ｋ系列の信号列に対する演算処理で、最小（Ｄsc＝１）でＲ×Ｋ系列の信号列を出力する構成で済む。また、演算部５４も、最小でＲ×Ｋ系列の信号列に対する演算処理を行なう構成で済む。
【０１０５】
したがって、例えば送信アンテナ数Ｎ＝１２８、受信アンテナ数Ｍ＝４の場合、レイヤ数Ｒを２とすれば、従来のＮ×Ｋ系列の信号列の演算処理を時間軸上で行なう場合に比べて、概略計算でＲ／Ｎ＝１／６４の回路規模で演算処理が可能であり、レイヤ数ＲをＭに等しくしても、Ｒ／Ｎ＝１／３２の回路規模で演算処理が可能となる。
【０１０６】
また、本実施形態の場合、フェージングの情報を周波数領域に変換するためのフーリエ変換処理が必要となるが、このフーリエ変換処理では、時間領域における各パスの遅延量は、周波数領域では各パスの周波数成分の回転速度に対応するため、従来時間領域で行なっていたメモリとリサンプルフィルタの組合せにより各パスに遅延を付与するハードウエアは、フーリエ変換における回転処理に置き換えられることになり、両者のハードウエアの規模を比較して本実施形態の方が格段に有利である。しかも、フーリエ変換の前に、ビームフォーミング処理を行なうことで、フーリエ変換処理する信号の系列数をＲ／Ｎに減らせるので、回路規模の増大は僅少で済む。
【０１０７】
前記実施形態では、ビームフォーミング処理と等価の演算処理をフェージング設定部５１の出力に対して行なっていたが、前記式（１６）は、次のように変形できる。
【０１０８】
［Ｓｆ］＝［Ｈ０］［Ｗ］［Ｓs′］
＝｛［Ｈ０］［Ｗ］｝×［Ｓs′］    ……（１６ａ）
＝［Ｈ０］×｛［Ｗ］［Ｓs′］｝    ……（１６ｂ）
【０１０９】
上記式（１６ａ）は、｛［Ｈ０］［Ｗ］｝の演算で得られた結果に対して［Ｓs′］を乗算するものであり、図７に示す試験装置３０′のように、フェージング設定部５１の出力に対するフーリエ変換処理をフーリエ変換部５３′で行い、その演算結果Ｈｔ（［Ｈ０］に相当）に対して、ビームフォーミング等価演算部５２′でビームフォーミング行列［Ｗ］の乗算処理を行なった結果ＨＷと、窓関数演算部３２の演算結果Ｆsym（［Ｓs′］に相当）とを演算部５４′で乗算処理する構成となる。
【０１１０】
また、上記式（１６ｂ）は、｛［Ｗ］［Ｓs′］｝の演算で得られた結果に対して［Ｈ０］を乗算するものであり、図８に示す試験装置３０″のように、フェージング設定部５１の出力に対するフーリエ変換処理をフーリエ変換部５３′で行った結果Ｈｔ（［Ｈ０］に相当）と、窓関数演算部３２の演算結果Ｆsym（［Ｓs′］に相当）にビームフォーミング等価演算部５２″でビームフォーミング行列［Ｗ］の乗算処理を行なって得られた結果ＦＷとを演算部５４″で乗算処理する構成となる。
【０１１１】
これらの試験装置３０′、３０″の場合も、伝搬路の特性と変調信号との乗算演算を周波数領域において行なっている点で前記実施形態と同じであり、時間領域信号生成のための演算処理を行なう回路規模を、従来装置に比べてＭ／Ｎに縮小でき、Ｍに対してＮが格段に大きいシステムを小規模に実現できる。
【０１１２】
なお、上記実施形態では、時間領域信号生成部３３に帯域制限フィルタ３３ａが設けられているが、これを省略して演算部５４の出力をフーリエ逆変換部３３ｂに直接入力してもよい。
【０１１３】
また、帯域制限フィルタの処理をフーリエ逆変換処理の後に時間領域で行なうことも可能であるが、その場合、フーリエ逆変換処理で得られた時間領域の信号に対して畳み込み演算処理を行なう必要がある。これに対し、本実施形態のように、帯域制限フィルタをフーリエ逆変換処理の前段に設けておけば、フィルタ処理を周波数領域での乗算処理で済ませることができ、畳み込み演算に比べて格段に少ない演算量で処理を実行でき、帯域制限フィルタを設ける場合であっても高速処理できる。
【０１１４】
また、上記説明は、マルチキャリア変調方式がＯＦＤＭの場合で説明したが、他のマルチキャリア変調方式のＵＦＭＣ、ＧＦＤＭ、ＦＢＭＣ等を用いたＭＩＭＯシステムについても本発明を同様に適用できる。
【０１１５】
特に、第４世代Ｅｖｏｌｕｔｉｏｎおよび第５世代の携帯電話方式で利用されることが期待される３Ｄ−ＭＩＭＯ／Ｍａｓｓｉｖｅ−ＭＩＭＯにおいては、基地局の送信アンテナ数の方が、移動機の受信アンテナ数よりも圧倒的に多い状況となっており、本発明が非常に有効である。
【符号の説明】
【０１１６】
１……試験対象、３０、３０′、３０″……ＭＩＭＯ方式システムの試験装置、３１……レイヤ周波数領域信号生成部、３２……窓関数演算部、３３……時間領域信号生成部、３３ａ……帯域制限フィルタ、３３ｂ……フーリエ逆変換部、３４……シフト加算部、５０……伝搬路シミュレータ、５１……フェージング設定部、５１ａ……ＡＷＧＮ、５１ｂ……ドップラフィルタ、５１ｃ……ＭＩＭＯ相関設定部、５１ｄ……補間器、５２、５２′、５２″……ビームフォーミング等価演算部、５３、５３′……フーリエ変換部、５４、５４′、５４″……演算部
【図１】
2018101841.tif 000008
【図２】
2018101841.tif 000009
【図３】
2018101841.tif 000010
【図４】
2018101841.tif 000011
【図５】
2018101841.tif 000012
【図６】
2018101841.tif 000013
【図７】
2018101841.tif 000014
【図８】
2018101841.tif 000015
【図９】
2018101841.tif 000016
【図１０】
2018101841.tif 000017