特許 6586129 ＭＩＭＯ方式システムの試験装置および試験方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6586129

(24)【登録日】2019年9月13日

(45)【発行日】2019年10月2日

(54)【発明の名称】ＭＩＭＯ方式システムの試験装置および試験方法

(51)【国際特許分類】

   H04B 17/29 20150101AFI20190919BHJP

   H04L 27/26 20060101ALI20190919BHJP

   H04B 7/0413 20170101ALI20190919BHJP

【ＦＩ】

   H04B17/29

   H04L27/26 100

   H04B7/0413

【請求項の数】6

【全頁数】31

(21)【出願番号】特願2017-120619(P2017-120619)

(22)【出願日】2017年6月20日

(65)【公開番号】特開2019-9493(P2019-9493A)

(43)【公開日】2019年1月17日

【審査請求日】2018年8月24日

(73)【特許権者】

【識別番号】000000572

【氏名又は名称】アンリツ株式会社

(74)【代理人】

【識別番号】100072604

【弁理士】

【氏名又は名称】有我 軍一郎

(74)【代理人】

【識別番号】100140501

【弁理士】

【氏名又は名称】有我 栄一郎

(72)【発明者】

【氏名】小林 武史

【審査官】 太田 龍一

(56)【参考文献】

【文献】 特開２０１７−１０８２３６（ＪＰ，Ａ）

【文献】 特開２０１４−１１６７４７（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１３／０００６６０１（ＵＳ，Ａ１）

【文献】 Study on channel model for frequencies from 0.5 to 100 GHz (Release 14)，3GPP TR38.901 V14.0.0，２０１７年 ３月，p.36 (7.5-22)，ＵＲＬ，https://www.3gpp.org/ftp/Specs/archive/38_series/38.901

【文献】 上橋 俊介，小川 恭孝，西村 寿彦，大鐘 武雄，時変動マルチユーザＭＩＭＯ環境における2段階の圧縮センシングを用いたチャネル予測，電子情報通信学会技術研究報告，２０１６年 ２月２４日，Vol.115 No.472，pp.357-362

【文献】 奥村 香菜子，小川 恭孝，西村 寿彦，大鐘 武雄，インパルス応答推定に基づくチャネル予測手法に関する検討，電子情報通信学会技術研究報告，２０１６年 ６月１５日，Vol.116,No.110，pp.81-86

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｈ０４Ｂ １７／２９

Ｈ０４Ｂ ７／０４１３

Ｈ０４Ｌ ２７／２６

(57)【特許請求の範囲】

【請求項1】

１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記複数Ｕの受信アンテナを通して受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める第１の伝搬路演算部（４１）と、
前記第１の伝搬路演算部の演算結果に対し、前記複数Ｕの受信アンテナが前記試験対象とともに移動することで生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める第２の伝搬路演算部（４２）と、
前記第２の伝搬路演算部の演算結果に対し、前記複数Ｕの受信アンテナの受信特性を表す受信アンテナ特性を乗じて、前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性を求める第３の伝搬路演算部（４３）と、
前記第３の伝搬路演算部の演算結果のうち、前記複数Ｕの受信アンテナ毎にドップラ周波数シフトが共通とみなせる前記素波を複数まとめて１単位とするとともに、それぞれドップラ周波数が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求めるフーリエ変換部（４４）と、
前記フーリエ変換部の演算結果に対して、前記レイヤ周波数領域信号生成部で生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記複数Ｕの受信アンテナまでの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成する伝搬信号演算部（５１）と、
前記伝搬信号演算部で得られる前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号の抽出処理を行なう窓関数演算部（５２）と、
前記窓関数演算部の演算結果に対し、前記複数Ｕの受信アンテナ毎にＬ単位の信号の加算処理を行い、前記複数Ｋのキャリア当りＵ系列の伝搬信号を生成するパス信号加算部（５３）と、
前記パス信号加算部の演算結果に対してフーリエ逆変換処理を行ない、前記複数Ｕの受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（５４）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記複数Ｕの受信アンテナを通して受信される連続した受信信号を生成するシフト加算部（５５）とを備えたことを特徴とするＭＩＭＯ方式システムの試験装置。

【請求項2】

１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記試験対象が有する前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナにそれぞれ入射するための入射波を生成して、電波暗室内でＵ'本のプローブアンテナを通して前記試験対象に前記入射波を与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める第１の伝搬路演算部（４１）と、
前記第１の伝搬路演算部の演算結果に対し、前記試験対象が移動することによって生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める第２の伝搬路演算部（４２）と、
前記第２の伝搬路演算部の演算結果のうち、前記試験対象からみた電波の到来方向が同一とみなせる前記素波を複数まとめて１単位とするとともに、それぞれ到来方向が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求めるフーリエ変換部（４４′）と、
前記フーリエ変換部の演算結果に対して、前記レイヤ周波数領域信号生成部で生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記試験対象が有する前記複数Ｕの受信アンテナの入射経路までの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当りＬ系列で生成する伝搬信号演算部（５１′）と、
前記伝搬信号演算部で得られたＬ系列の前記伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当りＬ系列の伝搬信号の抽出処理を行なう窓関数演算部（５２′）と、
前記窓関数演算部の演算結果に対し、前記Ｕ'本のプローブアンテナ毎にＬ単位の信号の重みづけ合成処理を行い、前記複数Ｋのキャリア当りＵ'系列の伝搬信号を生成する重み付け演算部（５３′）と、
前記重み付け演算部の演算結果に対してフーリエ逆変換処理を行ない、前記Ｕ'本のプローブアンテナから出力するための時間領域の信号を生成する時間領域信号生成部（５４′）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記Ｕ'本のプローブアンテナからそれぞれ出力されて前記試験対象に向けて入射する連続した入射波を生成するシフト加算部（５５′）とを備えたことを特徴とするＭＩＭＯ方式システムの試験装置。

【請求項3】

前記フーリエ変換部が、
前記Ｌ単位分の伝搬路の特性に対するフーリエ変換演算処理を行なうフーリエ変換演算手段（４４ａ）と、
前記フーリエ変換演算手段の演算結果に対する周波数軸方向の補間処理を行う補間手段（４４ｃ）とを含むことを特徴とする請求項１または請求項２記載のＭＩＭＯ方式システムの試験装置。

【請求項4】

１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナを通して受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成する段階と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める段階と、
前記ビームフォーミング特性が付与された伝搬路の特性に対し、前記複数Ｕの受信アンテナが前記試験対象とともに移動することで生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める段階と、
前記試験対象の移動状況を加味した伝搬路の特性に対し、前記複数Ｕの受信アンテナの受信特性を表す受信アンテナ特性を乗じて、前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性を求める段階と、
前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性の演算結果のうち、前記複数Ｕの受信アンテナ毎にドップラ周波数シフトが共通とみなせる前記素波を複数まとめて１単位とするとともに、それぞれドップラ周波数が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階と、
前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性に対して、前記生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記複数Ｕの受信アンテナまでの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成する段階と、
前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成された周波数領域の伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号の抽出処理を行なう段階と、
前記窓関数の周波数特性の畳み込み演算により抽出された前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号に対し、前記複数Ｕの受信アンテナ毎にＬ単位の信号の加算処理を行い、前記複数Ｋのキャリア当りＵ系列の伝搬信号を生成する段階と、
前記加算処理により生成された前記複数Ｋのキャリア当りＵ系列の伝搬信号に対してフーリエ逆変換処理を行ない、前記複数Ｕの受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記複数Ｕの受信アンテナを通して受信される連続した受信信号を生成する段階とを含むことを特徴とするＭＩＭＯ方式システムの試験方法。

【請求項5】

１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記試験対象が有する前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナにそれぞれ入射するための入射波を生成して、電波暗室内でＵ'本のプローブアンテナを通して前記試験対象に前記入射波を与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成する段階と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める段階と、
前記ビームフォーミング特性が付与された伝搬路の特性に対し、前記試験対象が移動することにより生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める段階と、
前記試験対象の移動状況を加味した伝搬路の特性のうち、前記試験対象からみた電波の到来方向が同一とみなせる前記素波を複数まとめて１単位とするとともに、それぞれ到来方向が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階と、
前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性に対して、前記生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記試験対象が有する前記複数Ｕの受信アンテナの入射経路までの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当りＬ系列で生成する段階と、
前記複数Ｋのキャリア当りＬ系列で生成された伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当りＬ系列の伝搬信号の抽出処理を行なう段階と、
前記窓関数の周波数特性の畳み込み演算により抽出された伝搬信号に対し、前記Ｕ'本のプローブアンテナ毎にＬ単位の信号の重みづけ合成処理を行い、前記複数Ｋのキャリア当りＵ'系列の伝搬信号を生成する段階と、
前記重み付け処理の演算結果に対してフーリエ逆変換処理を行ない、前記Ｕ'本のプローブアンテナから出力するための時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記Ｕ'本のプローブアンテナからそれぞれ出力されて前記試験対象に向けて入射する連続した入射波を生成する段階とを含むことを特徴とするＭＩＭＯ方式システムの試験方法。

【請求項6】

前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階は、
前記Ｌ単位分の伝搬路の特性に対するフーリエ変換演算処理を行なう段階と、
該フーリエ変換演算処理の演算結果に対する周波数軸方向の補間処理を行う段階とを含むことを特徴とする請求項４または請求項５記載のＭＩＭＯ方式システムの試験方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、基地局から移動体端末へのダウンリンク信号を、基地局側アンテナ数Ｓ、端末側アンテナ数Ｕで伝達するＭＩＭＯ（Multi Input Multi Output）方式に対応した端末あるいはその端末に内蔵される回路基板や集積回路等を試験対象とし、送受信アンテナ間に想定されるＳ×Ｕチャネルの伝搬路に対してフェージング処理を行なう機能を有する試験装置の回路規模を小さくするための技術に関する。

【背景技術】

【0002】

ＭＩＭＯ方式は、図１０に示すように、端末側へのダウンリンク信号Ｓtx1 〜ＳtxS を、Ｓ本(この例ではＳ＝４とする)の基地局側アンテナ（以下、送信アンテナ）Ａtx1 〜ＡtxS から送信し、Ｕ本(この例ではＵ＝２とする)の端末側アンテナ（以下、受信アンテナ）Ａrx1 〜ＡrxU で受信させる。

【0003】

したがって、各送信アンテナと各受信アンテナの間にＳ×Ｕの伝搬路（チャネル）が想定され、また各チャネルについて異なる複数Ｌ（例えばＬ＝４）のパスが想定される。パスを含めた各チャネルの伝搬特性をＨ(1,1,1〜L)〜Ｈ(S,U,1〜L)とすると、ＭＩＭＯ方式に対応した移動体端末やそれに用いる回路等を試験する場合には、ダウンリンク信号に対し、各チャネルの伝搬特性およびパスについての損失、遅延、ドップラシフト等の特性を加味した演算処理を行い、最終的にＵ本の受信アンテナＡrx1 〜ＡrxU から出力される受信信号Ｓrx1 〜ＳrxU を生成して試験対象１に与える必要がある。

【0004】

一方、近年では変調方式として、ＯＦＤＭ（Orthogonal Frequency Division Multiplexing）、ＵＦＭＣ（Universal Filtered Multicarrier）、ＧＦＤＭ（Generalized Frequency Division Multiplexing）、ＦＢＭＣ（Filtered Bank Multi-Carrier）等のマルチキャリア変調方式による高速な信号伝達が実現されており、このマルチキャリア変調方式とＭＩＭＯ方式との組合せにより、さらに高速な情報通信が可能なＭＩＭＯ方式システムが実現され、そのシステムを試験する装置が必要となる。

【0005】

また、次の世代（第５世代）の通信方式では、より高い周波数帯を用いることが提案されている。このように通信に使用する周波数帯が高くなると、アンテナ個々のサイズを小さく構成できるため、多数のアンテナ素子を縦横に配列したアレーアンテナ構造を採用して、それらのアンテナ素子に与えるダウンリンク信号の位相制御により、通信対象の移動体端末の存在する方向に効率的に電波を放射させる、所謂ビームフォーミングが可能となる。したがって、このような次世代の移動体端末を試験対象とする試験装置においては、アレー化された多数のアンテナについてのビームフォーミングの演算処理が必要となる。

【0006】

図１１は、マルチキャリア変調方式、ＭＩＭＯ方式およびアレーアンテナによるビームフォーミング処理を組合せたシステムを試験するための試験装置の構成例を示している。

【0007】

この試験装置１０は、マルチキャリア変調方式のうち、複数Ｋのサブキャリアを用いて端末との通信を行なうＯＦＤＭに対応したものであり、レイヤ周波数領域信号生成部１１において、試験対象に伝達しようとするＲ系列の伝送データ（レイヤあるいはストリームと呼ばれる）に対してそれぞれ複数Ｋのサブキャリア毎の変調信号（コンスタレーションデータ）Ｓsym(1,1)〜Ｓsym(1,K)、Ｓsym(2,1)〜Ｓsym(2,K)、……、Ｓsym(R,1)〜Ｓsym(R,K)を生成出力する。この変調信号Ｓsym は、ＯＦＤＭシンボル毎に、周波数軸上にＫ個のコンスタレーションデータを並べたデータをＲ系列分含む周波数領域の信号である。

【0008】

これらの変調信号Ｓsym はビームフォーミング処理部１２に入力され、Ｓ本の送信アンテナから出射される電波のビーム特性が所望特性となるように演算処理され、１送信アンテナ当り複数Ｋのサブキャリア毎のビームフォーミング処理信号Ｓbf(1,1)〜Ｓbf(1,K)、Ｓbf(2,1)〜Ｓbf(2,K)、……、Ｓbf(S,1)〜Ｓbf(S,K)に変換される。なお、以下の説明では、図面も含めて、ｊ個一組の信号Ｓx(i,1)〜Ｓx(i,j)をＳx(i,1〜j)と略記する場合がある。

【0009】

これらのビームフォーミング処理信号Ｓbf は、Ｓ組の時間領域信号生成部１３(1)〜１３(S)に入力される。各時間領域信号生成部１３(i)（ｉ＝１〜Ｓ）は、Ｋ個１組のビームフォーミング処理信号Ｓbf(i,1〜K)に対して、フーリエ逆変換（ＩＦＦＴ）処理、サイクリックプレフィクス（ＣＰ）付加処理、帯域制限処理等を行い、ＯＦＤＭ方式で規定された時間軸上の信号に変換する。

【0010】

これによって、各時間領域信号生成部１３(1)〜１３(S)からは、Ｓ本の送信アンテナＡtx1 〜ＡtxS に与えるための送信信号(ダウンリンク信号)Ｓtx1 〜ＳtxS が出力されることになる。

【0011】

これらの送信信号Ｓtx1 〜ＳtxS は、Ｓ×Ｕチャネルの伝搬路の特性を模擬する伝搬路シミュレータ２０に入力される。

【0012】

伝搬路シミュレータ２０は、Ｓ本の送信アンテナとＵ本の受信アンテナの間に形成されるＳ×Ｕのチャネルと、それら各チャネルについてそれぞれ複数Ｌのパスとを想定し、これらＳ×Ｕ×Ｌの各パスに所望の遅延とフェージングを付加し、Ｕ本の受信アンテナがそれぞれ受ける受信信号を仮想的に生成するものである。

【0013】

この伝搬路シミュレータ２０は、無線通信における受信レベル変動の分布を表すレイリーフェージングを与えるものであり、Ｓ系列の送信信号Ｓtx1 〜ＳtxS それぞれに設定される複数Ｌのパスに対し、所定の遅延を与えて出力する遅延設定部２１と、ドップラシフト、ＭＩＭＯ相関の情報が付与されたレイリー分布の伝搬路の特性を求めるフェージング設定部２２と、遅延設定部２１から出力される全パス分の遅延処理信号Ｓtx(1,1,1〜L)、Ｓtx(2,1,1〜L)、……、Ｓtx(S,U,1〜L)と、フェージング設定部２２で得られた伝搬特性Ｈ(1,1,1〜L）、Ｈ(2,1,1〜L）、……、Ｈ(S,U,1〜L）とを用いた積和演算（行列の乗算）により、Ｓ×Ｕ×Ｌの仮想的な伝搬路を経由してＵ本の受信アンテナで受信される受信信号Ｓrx1 〜ＳrxU を生成する演算部２３とを有している。

【0014】

ここで、遅延設定部２１は、例えばメモリによる１サンプル単位の遅延とリサンプルフィルタによる１サンプル以下の遅延の組み合わせで、各パスに所望の遅延を付与している。

【0015】

また、演算部２３の演算処理は、例えば、
Ｓrx1＝ΣＨ(1,1,i)・Ｓtx(1,1,i)＋ΣＨ(2,1,i)・Ｓtx(2,1,i)＋……
＋ΣＨ(S,1,i)・Ｓtx(S,1,i)
Ｓrx2＝ΣＨ(1,2,i)・Ｓtx(1,2,i)＋ΣＨ(2,2,i)・Ｓtx(2,2,i)＋……
＋ΣＨ(S,2,i)・Ｓtx(S,2,i)
……
ＳrxU＝ΣＨ(1,U,i)・Ｓtx(1,U,i)＋ΣＨ(2,U,i)・Ｓtx(2,U,i)＋……
＋ΣＨ(S,U,i)・Ｓtx(S,U,i)
となる。ただし、記号Σは、ｉ＝１〜Ｌまでの総和を表す。

【0016】

このようにして得られた受信信号Ｓrx1 〜ＳrxU を試験対象１に与えることで、試験装置側で設定した送受信アンテナ間の伝搬路の状態に対する試験対象１の動作を試験することができる。

【0017】

なお、伝搬路シミュレータは含まれていないが、上記したようにマルチキャリア変調方式とＭＩＭＯ方式とを組合せたシステムを試験するため試験装置は、例えば次の特許文献１に開示されている。

【先行技術文献】

【特許文献】

【0018】

【特許文献1】特開２０１４−９３７５８号公報

【発明の概要】

【発明が解決しようとする課題】

【0019】

上記構成の試験装置のように、ビームフォーミング処理を行なうシステムではアレー化された送信アンテナの数Ｓが例えば１２８のように非常に大きくなるため、それに伴って、フーリエ逆変換処理をＳ系列分並列的に行なう時間領域信号生成部１３が１２８組必要となり、回路規模が非常に大きくなってしまう。

【0020】

また、伝搬路シミュレータ２０の遅延設定部２１は、前記したように、メモリとリサンプルフィルタの組合せにより任意の遅延を付与するハードウエア構成を要するため、上記のように１２８系列の信号に対してそれぞれ設定した複数Ｌのパスに任意の遅延を付与するためには、やはりその回路規模が非常に大きくなり、装置が大型化し、製造コストおよび消費電力が大きくなってしまう。

【0021】

本発明は、上記問題を解決し、マルチキャリア変調方式、ＭＩＭＯ方式およびビームフォーミング処理を組合せたシステムで送信アンテナ数が多い場合でも、小さな回路規模、少ない消費電力で実現できる試験装置および試験方法を提供することを目的としている。

【課題を解決するための手段】

【0022】

前記目的を達成するために、本発明の請求項１のＭＩＭＯ方式システムの試験装置は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記複数Ｕの受信アンテナを通して受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める第１の伝搬路演算部（４１）と、
前記第１の伝搬路演算部の演算結果に対し、前記複数Ｕの受信アンテナが前記試験対象とともに移動することで生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める第２の伝搬路演算部（４２）と、
前記第２の伝搬路演算部の演算結果に対し、前記複数Ｕの受信アンテナの受信特性を表す受信アンテナ特性を乗じて、前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性を求める第３の伝搬路演算部（４３）と、
前記第３の伝搬路演算部の演算結果のうち、前記複数Ｕの受信アンテナ毎にドップラ周波数シフトが共通とみなせる前記素波を複数まとめて１単位とするとともに、それぞれドップラ周波数が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求めるフーリエ変換部（４４）と、
前記フーリエ変換部の演算結果に対して、前記レイヤ周波数領域信号生成部で生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記複数Ｕの受信アンテナまでの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成する伝搬信号演算部（５１）と、
前記伝搬信号演算部で得られる前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号の抽出処理を行なう窓関数演算部（５２）と、
前記窓関数演算部の演算結果に対し、前記複数Ｕの受信アンテナ毎にＬ単位の信号の加算処理を行い、前記複数Ｋのキャリア当りＵ系列の伝搬信号を生成するパス信号加算部（５３）と、
前記パス信号加算部の演算結果に対してフーリエ逆変換処理を行ない、前記複数Ｕの受信アンテナでそれぞれ受信される時間領域の信号を生成する時間領域信号生成部（５４）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記複数Ｕの受信アンテナを通して受信される連続した受信信号を生成するシフト加算部（５５）とを備えたことを特徴とする。

【0023】

また、本発明の請求項２のＭＩＭＯ方式システムの試験装置は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記試験対象が有する前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナにそれぞれ入射するための入射波を生成して、電波暗室内でＵ'本のプローブアンテナを通して前記試験対象に前記入射波を与えるＭＩＭＯ方式システムの試験装置において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成するレイヤ周波数領域信号生成部（３１）と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める第１の伝搬路演算部（４１）と、
前記第１の伝搬路演算部の演算結果に対し、前記試験対象が移動することによって生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める第２の伝搬路演算部（４２）と、
前記第２の伝搬路演算部の演算結果のうち、前記試験対象からみた電波の到来方向が同一とみなせる前記素波を複数まとめて１単位とするとともに、それぞれ到来方向が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求めるフーリエ変換部（４４′）と、
前記フーリエ変換部の演算結果に対して、前記レイヤ周波数領域信号生成部で生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記試験対象が有する前記複数Ｕの受信アンテナの入射経路までの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当りＬ系列で生成する伝搬信号演算部（５１′）と、
前記伝搬信号演算部で得られたＬ系列の前記伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当りＬ系列の伝搬信号の抽出処理を行なう窓関数演算部（５２′）と、
前記窓関数演算部の演算結果に対し、前記Ｕ'本のプローブアンテナ毎にＬ単位の信号の重みづけ合成処理を行い、前記複数Ｋのキャリア当りＵ'系列の伝搬信号を生成する重み付け演算部（５３′）と、
前記重み付け演算部の演算結果に対してフーリエ逆変換処理を行ない、前記Ｕ'本のプローブアンテナから出力するための時間領域の信号を生成する時間領域信号生成部（５４′）と、
前記時間領域信号生成部が生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記Ｕ'本のプローブアンテナからそれぞれ出力されて前記試験対象に向けて入射する連続した入射波を生成するシフト加算部（５５′）とを備えたことを特徴とする。

【0024】

また、本発明の請求項３のＭＩＭＯ方式システムの試験装置は、請求項１または請求項２記載のＭＩＭＯ方式システムの試験装置において、
前記フーリエ変換部が、
前記Ｌ単位分の伝搬路の特性に対するフーリエ変換演算処理を行なうフーリエ変換演算手段（４４ａ）と、
前記フーリエ変換演算手段の演算結果に対する周波数軸方向の補間処理を行う補間手段（４４ｃ）とを含むことを特徴とする。

【0025】

また、本発明の請求項４のＭＩＭＯ方式システムの試験方法は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナを通して受信される受信信号を生成して前記試験対象に与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成する段階と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める段階と、
前記ビームフォーミング特性が付与された伝搬路の特性に対し、前記複数Ｕの受信アンテナが前記試験対象とともに移動することで生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める段階と、
前記試験対象の移動状況を加味した伝搬路の特性に対し、前記複数Ｕの受信アンテナの受信特性を表す受信アンテナ特性を乗じて、前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性を求める段階と、
前記複数Ｕの受信アンテナの受信特性を加味した伝搬路の特性の演算結果のうち、前記複数Ｕの受信アンテナ毎にドップラ周波数シフトが共通とみなせる前記素波を複数まとめて１単位とするとともに、それぞれドップラ周波数が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階と、
前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性に対して、前記生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記複数Ｕの受信アンテナまでの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成する段階と、
前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列で生成された周波数領域の伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号の抽出処理を行なう段階と、
前記窓関数の周波数特性の畳み込み演算により抽出された前記複数Ｋのキャリア当り（Ｕ×Ｌ）系列の伝搬信号に対し、前記複数Ｕの受信アンテナ毎にＬ単位の信号の加算処理を行い、前記複数Ｋのキャリア当りＵ系列の伝搬信号を生成する段階と、
前記加算処理により生成された前記複数Ｋのキャリア当りＵ系列の伝搬信号に対してフーリエ逆変換処理を行ない、前記複数Ｕの受信アンテナでそれぞれ受信される時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記複数Ｕの受信アンテナを通して受信される連続した受信信号を生成する段階とを含むことを特徴とする。

【0026】

また、本発明の請求項５のＭＩＭＯ方式システムの試験方法は、
１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式、複数Ｓの送信アンテナ、複数Ｕの受信アンテナを用いるＭＩＭＯ方式および前記複数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、前記複数Ｓの送信アンテナからＮ個の散乱体を含む伝搬環境を経て前記試験対象が有する前記複数Ｕの受信アンテナに至る擬似的な伝搬路を想定し、該伝搬路を経由して前記複数Ｕの受信アンテナにそれぞれ入射するための入射波を生成して、電波暗室内でＵ'本のプローブアンテナを通して前記試験対象に前記入射波を与えるＭＩＭＯ方式システムの試験方法において、
前記試験対象に伝達するためのＲレイヤ分のデータ信号列に対して、それぞれ前記複数Ｋのキャリア毎の周波数領域の変調信号を（Ｒ×Ｋ）系列分生成する段階と、
前記複数Ｓの送信アンテナと前記データ信号列のレイヤの数で特定される（Ｓ×Ｒ）系列分のビームフォーミング特性に、前記複数Ｓの送信アンテナが送出する信号が該複数Ｓの送信アンテナの特性に従って散乱体１つ当りＭ本の素波として出力されて対応する散乱体に到達するまでの伝搬路を模擬する特性を乗じて、ビームフォーミング特性が付与された伝搬路の特性を求める段階と、
前記ビームフォーミング特性が付与された伝搬路の特性に対し、前記試験対象が移動することにより生じるドップラ周波数シフトを付与するための位相特性を乗じて、前記試験対象の移動状況を加味した伝搬路の特性を求める段階と、
前記試験対象の移動状況を加味した伝搬路の特性のうち、前記試験対象からみた電波の到来方向が同一とみなせる前記素波を複数まとめて１単位とするとともに、それぞれ到来方向が異なるＬ単位分の伝搬路の特性に対してフーリエ変換処理を行い、前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階と、
前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性に対して、前記生成された（Ｒ×Ｋ）系列分の変調信号を乗算して、前記複数Ｓの送信アンテナから前記試験対象が有する前記複数Ｕの受信アンテナの入射経路までの擬似的な伝搬路を経由した周波数領域の伝搬信号を前記複数Ｋのキャリア当りＬ系列で生成する段階と、
前記複数Ｋのキャリア当りＬ系列で生成された伝搬信号に対し、該各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、前記複数Ｋのキャリア当りＬ系列の伝搬信号の抽出処理を行なう段階と、
前記窓関数の周波数特性の畳み込み演算により抽出された伝搬信号に対し、前記Ｕ'本のプローブアンテナ毎にＬ単位の信号の重みづけ合成処理を行い、前記複数Ｋのキャリア当りＵ'系列の伝搬信号を生成する段階と、
前記重み付け処理の演算結果に対してフーリエ逆変換処理を行ない、前記Ｕ'本のプローブアンテナから出力するための時間領域の信号を生成する段階と、
前記生成した時間領域の信号を、前記窓関数の長さ分ずらして加算して、前記Ｕ'本のプローブアンテナからそれぞれ出力されて前記試験対象に向けて入射する連続した入射波を生成する段階とを含むことを特徴とする。

【0027】

また、本発明の請求項６のＭＩＭＯ方式システムの試験方法は、請求項４または請求項５記載のＭＩＭＯ方式システムの試験方法において、
前記複数Ｋのキャリア間隔毎の周波数領域での伝搬路の特性を求める段階は、
前記Ｌ単位分の伝搬路の特性に対するフーリエ変換演算処理を行なう段階と、
該フーリエ変換演算処理の演算結果に対する周波数軸方向の補間処理を行う段階とを含むことを特徴とする。

【発明の効果】

【0028】

このように、本発明の請求項１、４では、Ｒレイヤのデータ信号列についてそれぞれキャリア毎の変調信号を生成する一方で、ビームフォーミング特性と送信アンテナ特性とで決まる散乱体を含めた伝搬路の特性に、試験対象の移動に伴うドップラ周波数シフトの特性と受信アンテナの特性とを加味して、送信アンテナから各受信アンテナまでの伝搬路の特性を求め、これをドップラ周波数シフトが共通とみなせる伝搬路ごとにまとめてレイリー分布に従うフェージング特性を付与し、フーリエ変換処理により、キャリア間隔毎の周波数領域の伝搬路の特性を求める。そして、この伝搬路の特性とキャリア毎の変調信号との演算により、ドップラ周波数シフトが共通とみなせる伝搬路ごとの所望の伝搬路特性とビームフォーミング特性が付与された周波数領域での伝搬信号を生成し、この伝搬信号に対し、各系列に対応するドップラ周波数で回転する時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、窓関数の周波数特性の畳み込み演算を行ない、それぞれドップラ周波数が異なる系列を受信アンテナごとに加算して、これをフーリエ逆変換処理して時間領域の信号を生成し、これを窓関数の長さ分ずらして加算することで、各受信アンテナでそれぞれ受信される連続した受信信号を生成している。

【0029】

また、本発明の請求項２、５では、Ｒレイヤのデータ信号列についてそれぞれキャリア毎の変調信号を生成する一方で、ビームフォーミング特性と送信アンテナ特性とで決まる散乱体を含めた伝搬路の特性に、試験対象の移動に伴うドップラ周波数シフトの特性を加味して、送信アンテナから散乱体（クラスタ）までの伝搬路の特性を求め、これを到来方向が同一とみなせる伝搬路ごとにまとめてレイリー分布に従うフェージング特性を付与し、フーリエ変換処理により、到来方向が同一とみなせる伝搬路ごとの周波数領域の伝搬路の特性を求める。そして、この伝搬路の特性とキャリア毎の変調信号との演算により、到来方向が同一とみなせる伝搬経路ごとの所望の伝搬路特性とビームフォーミング特性が付与された周波数領域での伝搬信号を生成し、この伝搬信号に対し、各系列に対応するドップラ周波数で回転する時間領域での窓関数の乗算による信号切出しに相当する周波数領域での処理として、窓関数の周波数特性の畳み込み演算を行ない、それぞれの到来方向が異なる系列に重み付けを施し、各プローブアンテナから送信する信号に変換して、これをフーリエ逆変換処理して時間領域の信号を生成し、これを窓関数の長さ分ずらして加算することで、電波暗室内の各プローブアンテナからそれぞれ出力されて試験対象が有する各受信アンテナに向けて入射するための連続した入射波を生成している。

【0030】

このように、本発明では、伝搬路の特性と変調信号との乗算処理を周波数領域において行ない、その演算結果から時間領域信号を生成しているので、従来方式のように、送信アンテナ毎の周波数領域の信号をフーリエ逆変換して時間領域の信号に変換してから伝搬路特性を付与する場合に比べて、フーリエ逆変換を行なう回路および伝搬路特性を生成する回路の規模を格段に小さくすることができる。

【0031】

特に、伝搬路の特性を求める演算において、ビームフォーミング特性と散乱体を考慮した送信アンテナ特性との演算を優先的に行なっているので、それ以降のドップラ周波数シフトの演算や受信アンテナ特性の演算の処理では、送信アンテナ数に依存した演算処理がなくなるため、次の世代（第５世代）の通信方式で提案されているように、基地局側（送信側）がアレーアンテナのような多数のアンテナ素子を用いるシステムを試験する場合、後続の演算処理の規模を格段に小さくでき、極めて有効である。

【0032】

ただし、本発明では、散乱体を考慮し、ドップラ周波シフト（到来方向）が共通とみなせるパスの伝搬路特性を周波数領域に変換するためのフーリエ変換処理が必要となるが、このフーリエ変換処理では、時間領域における各パスの遅延量は、周波数領域では各パスの周波数成分の回転速度に対応するため、従来時間領域で行なっていたメモリとリサンプルフィルタの組合せにより各パスに遅延を付与するハードウエアは、フーリエ変換における回転処理に置き換えられることになり、両者のハードウエアの規模を比較して本発明の方が格段に有利となる。

【0033】

また、請求項３、６では、フーリエ変換によりキャリア間隔毎の周波数領域での伝搬路の特性を求めるための構成として、Ｌ単位分の伝搬路の特性に対するフーリエ変換演算処理と、その演算結果に対する周波数軸方向の補間処理との２段階処理を用いているので、フーリエ変換の周波数領域でのサンプリング間隔を、例えば、時間領域の伝搬路特性の最大遅延量で決まる必要最低限の間隔とし、その後で周波数領域での伝搬路特性を補間する２段階のパイプライン処理とすることが可能となり、各処理にかけられる時間が延び、共通回路リソースを複数回使い回すことで、フーリエ変換における回路規模を小さくすることができる。

【0034】

また、回路規模を抑えつつも、従来の時間領域で処理する伝搬路シミュレータの処理結果を精度よく再現できる構成となっている。

【0035】

即ち、ＩＦＦＴ処理間隔といった比較的粗い時間間隔ごとの伝搬路特性を時間領域で計算して、それをフーリエ変換することでＩＦＦＴ時間間隔ごとの伝搬路の周波数特性を求める。その特性をデータに乗算するだけでは、ＩＦＦＴ時間間隔の中での伝搬路特性を一定値で近似することになるが、本発明では、ドップラ周波数シフト（あるいは到来方向）が共通とみなすことができる経路ごとのデータに周波数領域でドップラシフト分の周波数シフトを与えているため、ＩＦＦＴ時間間隔の中での伝搬路特性を時間変化するものにすることができる。従来の時間領域で処理する伝搬路シミュレータと比較した再現性の精度は、ドップラ周波数シフト（あるいは到来方向）を共通とみなす近似の度合いに応じて調整することができる。

【図面の簡単な説明】

【0036】

【図1】本発明の原理を説明するためのタイミング図

【図2】時間領域の窓関数の一例を示す図

【図3】時間領域の窓関数の別の例を示す図

【図4】ドップラ周波数シフトを付与する場合の周波数シフト分で回転する時間領域の窓関数の周波数特性の例を示す図

【図5】本発明の実施形態の構成を示す図

【図6】本発明の実施形態の要部の構成図

【図7】試験環境の一例を示す図

【図8】試験環境の別の例を示す図

【図9】本発明の他の実施形態の構成を示す図

【図10】マルチパスＭＩＭＯの伝搬路の一例を示す図

【図11】従来装置の構成図

【発明を実施するための形態】

【0037】

以下、図面に基づいて本発明の実施の形態を説明するが、具体的な構成を説明する前に、本発明の試験装置の原理について説明する。

【0038】

本発明は、前記したＯＦＤＭ、ＵＦＭＣ、ＧＦＤＭ、ＦＢＭＣなどのマルチキャリア変調方式で、Ｓ×Ｕ ＭＩＭＯ（Ｓ＞Ｕ）を実施する場合の伝搬路シミュレータとして適用できるものであり、３Ｄ−ＭＩＭＯ／Ｍａｓｓｉｖｅ−ＭＩＭＯのように、送信アンテナ数が受信アンテナ数に比べて非常に多い場合に特に有効である。以降では、変調方式としては、主にＯＦＤＭを念頭に置いて説明する。

【0039】

本発明では、次式（１）のように、時間スパンＴｃの中で差が無視できる程度のドップラ周波数の差Δｆ_ｄを同一とみなして周波数領域でＭＩＭＯ伝搬路処理を施すというものである。

【0040】

Ｔｃ≪１／Δｆ_ｄ ……（１）
ただし、Δｆ_ｄは、同一とみなすことができるドップラ周波数の上限と下限の差である。

【0041】

例えば、ＯＦＤＭの場合には、次式（２）のように、１ＯＦＤＭシンボル長Ｔsym（＝有効データ長＋サイクリックプレフィクス長）をＰ分割した長さＴｃが式（１）を満たすようにする（Ｐ分割は必ずしも等分割でなくても良い）。

【0042】

Ｔｃ＝Ｔsym／Ｐ （Ｐ＝１，２，３，……） ……（２）

【0043】

図１は時間領域におけるＰ＝２の例を示すものであり、図１の（ａ）に示すＯＦＤＭの信号列に対して、図１の（ｂ１）〜（ｂ４）のように、１シンボル長Ｔsym をＴｃ＝Ｔsym／２で２分割するような矩形窓関数を乗算して切り取った波形に対し、マルチパス伝搬処理、フィルタ処理などを施した波形を図１の（ｃ１）〜（ｃ４）のようにＴｃ分ずれた状態で得て、加算処理することで最終的な送信信号を得る。

【0044】

実際に用いるローカライゼーション（信号切出し）用の窓関数の長さＴｃ′は、対応する周波数特性の広がりを抑えるために端部を丸めることによって、Ｔｃより若干大きくしてもよく、この窓関数のタイミングを順次Ｔｃずつずらしながら乗算した波形のそれぞれに、ＭＩＭＯチャネルのマルチパス伝搬路処理、フィルタ処理などを施す。また、分割された一つ一つの波形の時間長は、マルチパスの遅延時間分およびフィルタ処理による広がり分ＴｄだけＴｃ′よりも長くなる。それらをＴｃずつずらしながら加算した波形を処理結果とする考え方である。処理結果の波形は、Ｓ×Ｕ ＭＩＭＯ伝搬路の場合には、Ｕ系統分計算される。

【0045】

図２は、ローカライゼーション用の窓関数（区間長：Ｔｃ′）の詳細を示すもので、例えば、１ＯＦＤＭシンボルを複数個に分割する場合などに利用可能である。ナイキスト基準を満たすような特性になっており、この窓関数をＴｃずつずらした区間が連続的につながるような特性である。この時間領域における窓関数のロールオフが大きいほど、周波数領域での広がりが抑えられて、後述する窓関数演算部５２におけるフィルタのタップ数を抑えることができる。

【0046】

また、変調方式がＦＢＭＣで、１シンボル情報が広がった図３のような窓関数の時間長（Ｖ・Ｔｃ）（Ｖはオーバラッピングファクタ）において、ドップラ周波数差が十分小さいとみなせる範囲の経路のドップラ周波数シフト分を同一として処理して、本発明を適用することができる。

【0047】

窓関数演算により変調信号に対するドップラ周波数シフトを付与する場合には、ＣＰ＋有効データの長さ分（またはそれを分割した長さ分）の窓関数の周波数特性のサンプル位置を周波数方向にシフトした関数を畳み込むことで、変調信号に対するドップラ周波数シフトを実現できる（図４は、その窓関数の周波数特性であり、その周波数軸上でのサンプル位置をドップラ周波数シフト分だけずらした図）。

【0048】

上記処理は時間軸上で想定したものであるが、本発明は、最後の加算処理以外を周波数領域で等価な処理を実施するとともに、従来ではレイヤ信号に対して行なっているビームフォーミング処理と等価な処理を伝搬路の特性に対して行なうことにより、試験装置全体の回路規模を縮小している。

【0049】

次に、本発明の試験装置が行なう演算処理について説明する。
テクニカルレポート３ＧＰＰ ＴＲ３８．９００には、第５世代のシステムシミュレーションで使用されている、Ｕ個中の受信アンテナ番号ｕ、Ｓ個中の送信アンテナ番号ｓ、Ｎ個中のクラスタ（散乱体）番号ｎ、各クラスタで散乱するＭ個中の素波（Ｒａｙと表記される場合もある）番号ｍの伝搬路モデル(クラスタモデル)の非見通し波伝搬路係数Ｈ^NLOSは、次式（５）で定義されている。

【0050】

【数1】

【0051】

上式（５）において、Ｐｎはｎ番目のクラスタから出射される電波のパワー、Ｆ_rx,u,θ(…)、Ｆ_rx,u,φ(…)は、受信アンテナの垂直偏波成分および水平偏波成分に対する特性、ｅｘｐ（ｊΦ^θθ）、ｅｘｐ（ｊΦ^φφ）は偏波が変化しない成分の位相係数、ｅｘｐ（ｊΦ^θφ）、ｅｘｐ（ｊΦ^φθ）は偏波が変化する成分の位相係数であり、√ｋ^−１はその割合を示す。また、Ｆ_tx,u,θ(…)、Ｆ_tx,u,φ(…)は、送信アンテナの垂直偏波成分および水平偏波成分に対する特性、ｅｘｐ{ｊ2π(ｒ^T_rx,n,m・d_rx,u ）/λ₀}、ｅｘｐ{ｊ2π(ｒ^T_tx,n,m・d_tx,s ）/λ₀}は、受信アンテナおよび送信アンテナの位置に依存する位相のずれを示し、ｅｘｐ{ｊ2π(ｒ^T_rx,n,m・ｖ）ｔ/λ₀}は、ドップラ効果による周波数のずれを示している。

【0052】

なお、上記式で、上に記号＾が付されたｒ^T_rx,n,m は、受信機側のｎ番目のクラスタで散乱したｍ番目の素波の受信アンテナへの到来方向の単位ベクトルを示し、上に記号^が付されたr^T_tx,n,m は送信機側のｎ番目のクラスタで散乱されるｍ番目の素波の送信アンテナからの到来方向の単位ベクトルを示し、上に記号−が付されたd_rx,u は受信アンテナ素子ｕの位置ベクトルを示し、上に記号−が付されたd_tx,s は送信アンテナ素子ｓの位置ベクトルを示し、上に記号−が付されたｖは、移動速度ベクトルを示す。

【0053】

ここで、上記式（６）を構成する要素を次のように定義する。
（ａ）送信アンテナ特性と、送信アンテナから放射されてｎ番目のクラスタで散乱されるｍ番目の素波の伝搬経路による位相回転とを組合せた要素ｆ_rx,s,n,m を次のように定義する。

【0054】

【数2】

【0055】

（ｂ）ｎ番目のクラスタで散乱されるｍ番目の素波の位相と、その素波に与えられるドップラ周波数シフトとを組合せた要素Ｈ_n,m を次のように定義する。

【0056】

【数3】

【0057】

（ｃ）受信アンテナ特性と、ｎ番目のクラスタで散乱されるｍ番目の素波が受信アンテナで受信されるまでの位相回転とを組合せた要素ｆ_rx,u,n,m を次のように定義する。

【0058】

【数4】

【0059】

さらに、次のように定義する。

【0060】

【数5】

【0061】

ここで、ξ_tx,s,n の縦にＭ個並んだ行列要素のｆ_rx,s,n,m には、式（６）のように、偏波（水平偏波と垂直偏波）の項が縦に２つ並んだ行列が含まれるため、ξ_tx,s,n 全体としては２Ｍ×１行列となる。また、Ψ_n の対角行列要素のＨ_n,m には、式（７）のように、偏波の項が縦横に２ずつ並んだ行列が含まれるため、Ψ_n 全体としては２Ｍ×２Ｍ行列となる。また、Ξ_rx,u,n の対角行列要素のｆ_rx,u,n,m には、式（８）のように、偏波の項が縦に２つ並んだ行列の転置が含まれるため、ｆ_rx,u,n,m 全体としてはＭ×２Ｍ行列となる。。

【0062】

上記のように定義した上で、送信アンテナから試験対象の受信アンテナあるいは受信アンテナとフロントエンドまでを含めた伝搬路を模擬し、試験対象のアンテナ端子あるいは中間周波信号入力端子に対してケーブルを介して試験用信号を与える試験環境を想定した場合、ｕ番目の受信アンテナのサブキャリアデータｋに対して乗算すべき周波数領域における伝搬路の特性を表すＬ×Ｒ行列の式は、次式（１２）で表される。

【0063】

【数6】

【0064】

上記式（１２）において、ｗ_1，1 〜ｗ_Ｓ,Ｒ を要素とするＳ×Ｒ行列は、ビームフォーミング特性を示す行列であり、式（１２）を構成する５つの行列を簡略化し、
Ω≡［ζ］［Ξ］［Ψ］［ξ］［ｗ］
と表すと、そのうち４つの行列［Ξ］、［Ψ］、［ξ］、［ｗ］の演算により、送信アンテナから受信アンテナに至る伝搬経路の時間領域の特性が得られ、これをドップラ周波数が共通とみなせるもの同士をまとめて、フーリエ変換処理を行なうための行列［ζ］の演算を加えることにより、レイリー分布特性が付与された伝搬路の周波数領域での特性が得られる。

【0065】

なお、[ζ_ijk]は、Ｌ×（Ｎ・Ｍ）の行列であり、
ｉ＝｛１，２，…，Ｌ｝
ｊ＝｛１，２，…，Ｎ・Ｍ｝＝（全クラスタに含まれる素波のインデックス）

【0066】

ζ_ijk は、以下のように表される。

【0067】

【数7】

【0068】

ここで、式（１３）の｛ａｎｇｌｅ……ａｎｇｌｅ(i)｝は、試験対象からｎ番目のクラスタのｍ番目の素波を見た方向をＬ通りの方向に量子化するとｉ番目の方向であるときに「１」、それ以外のとき「０」となる比較演算であるとする。また、τ_n,m は、ｎ番目のクラスタのｍ番目の素波から試験対象までの伝搬遅延、Δｆ_０ は、伝搬路の特性を周波数領域に変換するためのフーリエ変換処理結果の周波数方向のサンプル間隔（時間領域の遅延の広がりからエリアシングが生じないと判断される必要十分な間隔）であり、ｅ^{−ｊ２π…} は、周波数インデックスｋ_０によって回転することで、時間領域におけるτ_n,m 分の遅延を与えるものである。

【0069】

また、Ｋ0 は、フーリエ変換処理後の補間処理において、Ｋ個のサブキャリアで構成されるＯＦＤＭ変調帯域をＢ分割したときの各帯域毎のサブキャリア数Ｋ／Ｂそれぞれの周波数における特性を計算する際に必要となる範囲をカバーするのに必要な数とする。つまり、Ｋ／Ｂの周波数幅に対応するΔf_０ 間隔のサンプル数よりも、前後に補間処理のインパルス応答の半分ずつ程度長い範囲をカバーする数とする。つまり、｛ｎ，ｍ｝方向がｉ番目のドップラ周波数と見なせる場合には、その成分に対応する要素に遅延に相当する回転を与えてｉ行目に足しこむことを表している。ここで、Ｂは、ＯＦＤＭ変調帯域をＢ分割した帯域ごとにビームフォーミング特性を別にすることを考えたものである。全サブキャリアＫに渡って全て共通のビームフォーミングを施す場合にはＢ＝１とする。

【0070】

本発明は、上記式（１２）に基づいて試験に必要な信号を生成するものである。以下、その実施形態を説明する。図５は、本発明の実施形態の試験装置３０の構成を示している。

【0071】

試験装置３０は、１つの移動体端末に対する通信に複数Ｋのキャリアを用いるマルチキャリア変調方式で且つ送信アンテナ数Ｓ、受信アンテナ数ＵのＭＩＭＯ方式およびアンテナ数Ｓの送信アンテナによる放射ビーム特性を設定するビームフォーミング処理方式を採用するシステムを試験対象とし、送信アンテナと受信アンテナの間にＳ×Ｕのチャネルと各チャネルにそれぞれ複数Ｌのパス（送受信のアンテナ間にＮ個の散乱体が存在する伝搬環境で、各散乱体から散乱されたＭ個ずつの素波が入射し、それをＬ通りにまとめたものと仮定）を含む擬似的な伝搬路を想定し、その伝搬路を経由してＵ個の受信アンテナで受信される受信信号を生成して試験対象に与えるＭＩＭＯ方式システムの試験装置である。なお、以下の説明は、マルチキャリア変調方式がＯＦＤＭの場合とし、ＯＦＤＭでは端末との通信に用いる複数のキャリアを「サブキャリア」と呼んでいるので、以下の説明でもこの「サブキャリア」と記す。

【0072】

この試験装置３０は、レイヤ周波数領域信号生成部３１、伝搬路特性演算部４０、伝搬信号演算部５１、窓関数演算部５２、パス信号加算部５３、時間領域信号生成部５４、シフト加算部５５を有している。

【0073】

レイヤ周波数領域信号生成部３１は、試験対象に伝達しようとするＲ系列の伝送データ（レイヤあるいはストリームと呼ばれる）に対してそれぞれ複数Ｋのサブキャリア毎の変調信号（コンスタレーションデータ）Ｓsym(1,1)〜Ｓsym(1,K)、Ｓsym(2,1)〜Ｓsym(2,K)、……、Ｓsym(R,1)〜Ｓsym(R,K)を生成出力する。この変調信号は、ＯＦＤＭシンボル毎に、周波数軸上にＫ個のコンスタレーションデータを並べたデータをＲ系列分含む周波数領域の信号である。なお、レイヤ数Ｒは、原理上、試験対象の受信アンテナ数Ｕ以下の値である。

【0074】

ここで、そのコンスタレーションデータを、次のようにＳsym,r,kという記号の複素数で表すことにする。

【0075】

Ｓsym,r,k ……（１４）
sym ：ＯＦＤＭシンボル番号
ｒ＝｛１，２，３，……，Ｒ｝：送信レイヤ番号インデックス
ｋ＝｛１，２，３，……，Ｋ｝：サブキャリア番号インデックス

【0076】

なお、図５では、信号Ｓsym,r,kのインデックスが分かりやすいように、Ｓsym（r,k）の形で表す（他の信号についても同様）。

【0077】

また、コンスタレーションデータが配置される間隔（サブキャリア間隔）をｆscとする。ｆscはＯＦＤＭシンボル長Ｔsym（＝有効データ長＋サイクリックプレフィクス長）と次の関係にある。
有効データ長＝１／ｆsc ……（１５）
Ｔsym＝（１／ｆsc）＋サイクリックプレフィクス長 ……（１６）

【0078】

一方、送信アンテナから散乱体を含む伝搬環境を経て受信アンテナに至る伝搬路についての周波数領域での特性を求める伝搬路特性演算部４０は、各受信アンテナに到来する電波のドップラ周波数（Ｌ種類に量子化）が同一とみなすことができる複数の素波の集まりとしてのレイリーノイズ生成処理およびビームフォーミング等価処理を、ＦＦＴ単位の時間間隔ごとに実施するものであり、第１の伝搬路演算部４１〜第３の伝搬路演算部４３により構成される。これら各演算部の処理はＦＦＴ単位の時間ごとのパイプライン処理として実装できる。

【0079】

第１の伝搬路演算部４１は、レイヤ毎に複数の送信アンテナから放射されるビームの特性を表すビームフォーミングのデータと、送信アンテナから送信アンテナ特性に従って出力されて、散乱体１つ当りＭ本の素波が散乱体に到達するまでの伝搬路を模擬する伝搬路特性のデータとに基づいて、送信アンテナから散乱体に至る伝搬路にビームフォーミング特性が付与された特性を求める。なお、この実施例では、送信アンテナから出力された信号を、互いに直交する二つの偏波成分（水平偏波と垂直偏波）に区別して説明するが、単一の直線偏波であってもよい。

【0080】

この第１の伝搬路演算部４１は、具体的には、図５に示しているように、ビームフォーミング設定手段４１ａ、送信アンテナ特性設定手段４１ｂおよび演算手段４１ｃにより構成される。ビームフォーミング設定手段４１ａは、ＯＦＤＭ変調波帯域をＢ分割して、それぞれの分割帯域毎にビームフォーミング特性が異なる可能性を考慮し、送信アンテナ数Ｓとレイヤ数Ｒの行列要素からなるビームフォーミング行列をＢ組設定できるようになっている。ＯＦＤＭ変調波帯域を分割しない場合、Ｂ＝１とする。

【0081】

ビームフォーミング設定手段４１ａおよび送信アンテナ特性設定手段４１ｂに設定されるデータは、前記式（１２）で示したＳ×Ｒ行列の［ｗ］と、前記したように２つの偏波成分を含む２Ｍ×１行列の要素ξ_tx,s,m がＮ×Ｓ行列並んで、全体として２・Ｎ・Ｍ×Ｓ行列で構成される［ξ］の要素であり、それらの行列に対する演算手段４１ｃの行列乗算（積和演算）により、二つの偏波成分についてそれぞれ（Ｎ・Ｍ）×Ｒ行列のデータがＢ通りずつ生成されることになる。これを図５では、（２ＮＭ＊Ｒ）＊Ｂで示している。

【0082】

したがって、この第１の伝搬路演算部４１の演算結果［Ａ１］には、送信アンテナ数Ｓに依存した系列数が無くなり、第５世代として提案されているように送信アンテナ数が極めて多い（例えば１２８等）システムであっても、後続の演算処理を小規模な構成で短期間に行なうことができる。なお、送信アンテナ数Ｓ以外の要素の数値例について言えば、受信アンテナ数Ｕ＝２、レイヤー数Ｒ＝２、散乱体数Ｎ＝２４、素波数Ｍ＝２０等である。

【0083】

なお、通常の送信機では、ビームフォーミング処理は送信データに対して施すが、フェージングを含むシミュレータの処理としては、伝搬路特性側との演算を先にしてＳ個の経路を持つ演算部の数を減らす。つまり、この部分の演算量は送信アンテナ数Ｓに比例したものになるが、この実施形態のように、ビームフォーミング特性と伝搬路特性側との演算を最初に実行することで、Ｓ個の経路を持つ演算部をこの演算部のみとすることができる。この演算処理は、ＦＦＴを実施する時間間隔という比較的長い時間で演算すれば良く、乗算器を時分割で何度も使うことを考えると、実装上の乗算器数は小さくて済む。

【0084】

第２の伝搬路演算部４２は、第１の伝搬路演算部４１の演算結果［Ａ１］に対し、受信アンテナを含めて試験対象が移動することにより生じるドップラ周波数シフトを付与した伝搬路の特性を求める。この第２の伝搬路演算部４２は、ドップラ周波数シフトにより回転する位相のデータを要素とする前記行列［Ψ］を与えるドップラシフト特性設定手段４２ａと、第１の伝搬路演算部４１の演算結果［Ａ１］に対して、ドップラシフト特性設定手段４２ａにより設定された行列［Ψ］の積和演算を行なう演算手段４２ｂにより構成されている。

【0085】

ここで、前記したように、第１の伝搬路演算部４１の演算結果［Ａ１］＝［ξ］［ｗ］は、二つの偏波成分を含めてＢ通りの（２・Ｎ・Ｍ）×Ｒ行列であり、ドップラシフト特性設定手段４２ａにより設定される行列［Ψ］は、前記したように、２Ｍ×２Ｍ行列の要素Ｈ_n,m がＮ×Ｎ行列並んで、全体として（２・Ｎ・Ｍ）×（２・Ｎ・Ｍ）の対角行列であるので、その演算結果［Ａ２］は、Ｂ通りの（２・Ｎ・Ｍ）×Ｒ行列となる。

【0086】

また、第３の伝搬路演算部４３は、第２の伝搬路演算部４２で得られた演算結果［Ａ２］に対し、二つの偏波成分に対する受信アンテナの受信特性を表す受信アンテナ特性を付与して受信アンテナを含めた伝搬路の特性を求める。この第３の伝搬路演算部４３は、受信アンテナ特性のデータを要素とする前記行列［Ξ］を与える受信アンテナ特性設定手段４３ａと、第２の伝搬路演算部４２の演算結果［Ａ２］に対して、受信アンテナ特性設定手段４３ａにより設定された行列［Ξ］の積和演算を行なう演算手段４３ｂにより構成されている。

【0087】

ここで、受信アンテナ特性設定手段４３ａにより設定された行列［Ξ］は、前記したように、Ｍ×２Ｍ行列の要素ｆ_rx,u,n,mがＮ×Ｎ行列並んでいるが、それが全受信アンテナ分Ｕ個あり、全体として（Ｕ・Ｎ・Ｍ）×（２・Ｎ・Ｍ）の対角行列であるので、演算手段４３ｂによる行列［Ξ］、［Ａ２］の演算結果［Ａ３］は、２つの偏波成分が受信アンテナで合成されたことにより、Ｂ通りのＵ・Ｎ・Ｍ×Ｒ行列で示される。

【0088】

このようにして得られた伝搬路の特性［Ａ３］＝［Ξ］［Ψ］［ξ］［ｗ］は、各送信アンテナから所望のビームフォーミング特性が付与されて放射され、Ｎ個全ての散乱体からそれぞれＭ個ずつの経路で所望のドップラ周波数シフトが付与された状態で各受信アンテナに到来して受信されるまでの伝搬路の特性を示している。

【0089】

この伝搬路の特性は、フーリエ変換部４４に与えられる。フーリエ変換部４４は、第３の伝搬路演算部４３で得られた伝搬路の特性のうち、各受信アンテナ毎にドップラ周波数シフトが共通とみなせる経路をまとめて１単位とすることでレイリー分布に従うフェージング特性を付与するとともに、全体をＬ単位にグループ化した伝搬路の特性に対して、フーリエ変換処理を行ない、サブキャリア間隔毎の周波数領域での伝搬路の特性に変換する。

【0090】

これを実現するために、フーリエ変換部４４には、ドップラ周波数シフトが共通とみなせるもの同士をまとめてフーリエ変換処理を行なうフーリエ変換演算手段４４ａと、その演算に用いるパラメータを設定するための係数設定手段４４ｂと、フーリエ変換演算手段４４ａの演算結果に対して補間処理を行なう補間手段４４ｃとを有している。

【0091】

フーリエ変換演算手段４４ａは、係数設定手段４４ｂに設定されたパラメータを用いて、演算結果［Ａ３］に対する前記［ζ］の演算処理を行なう。前記したように、［Ａ３］はＢ通りのＵ・Ｎ・Ｍ×Ｒ行列であり、［ζ］は、Ｕ個の受信アンテナそれぞれにＬ×（Ｎ・Ｍ）行列であるから、その演算結果［Ａ４］は、Ｕ・Ｌ×Ｒ×Ｋ０行列のＢ通りとなる。

【0092】

なお、係数設定手段４４ｂには、上記フーリエ変換処理に必要な係数、Ｐｎ、Ｍ、τ_n,m 、ＡoＡ_n,m（Angle of Arrival：クラスタｎの素波ｍが到来する角度方向の情報）等が予め設定される。

【0093】

上記したように、フーリエ変換部４４は、演算結果［Ａ３］をＵ本の各受信アンテナから見たときのドップラ周波数（Ｌ種類に量子化）が同一とみなすことができる素波（Ｒａｙ）同士をまとめて、伝搬路の時間拡がりから決まる（時間領域の遅延拡がりからエリアシングが生じないと判断される）必要十分なサンプリング間隔周波数領域信号に変換する。つまり、Ｕアンテナ毎に、時間領域での素波（Ｎ・Ｍ個）の遅延に対応した回転を与え、同一のドップラ周波数とみなせる素波同士を加算するものであり、その処理結果は、Ｕ・Ｌ×Ｒ行列で、Ｋ０通りとなる（Ｋ０は周波数方向のサンプル数）。ここで、ＬおよびＫ０を小さくすることで、現実的な回路リソースに圧縮することができる。遅延を周波数領域における回転として実施することで、従来より行なわれている「サンプル遅延＋リサンプルフィルタ」によって直接ＩＱデータに遅延を施す方式よりも、回路規模を小さくすることができる。

【0094】

補間手段４４ｃは、フーリエ変換演算手段４４ａで得られたフーリエ変換結果［Ａ４］を、周波数軸方向に補間して、サンプル間隔がＯＦＤＭ信号のサブキャリア間隔となるようにする。この補間処理により、周波数領域におけるＵ・Ｌ×Ｒ×Ｋ行列の伝搬路特性［Ａ５］が得られる。

【0095】

上記のように、このフーリエ変換部４４は、フーリエ変換演算手段４４ａと補間手段４４ｃを有しているので、フーリエ変換の処理を実施する際の周波数領域でのサンプリング間隔を、例えば、時間領域の伝搬路特性の最大遅延量で決まる必要最低限の間隔とし、その演算結果に対し、補間手段４４ｃにより周波数領域での伝搬路特性を補間するという２段階のパイプライン処理が可能となり、各処理にかけられる時間が延び、共通回路リソースを複数回使い回すことで、フーリエ変換における回路規模を小さくすることができる。なお、フーリエ変換演算手段４４ａの周波数領域でのサンプリング間隔をより小さくして補間手段４４ｃによる補間処理を省略する構成も可能である。

【0096】

伝搬信号演算部５１は、フーリエ変換部４４で求められた伝搬路特性［Ａ５］を、レイヤ周波数領域信号生成部３１で生成されたＲ×Ｋ系列分の変調信号Ｓsymに乗算して、サブキャリア当りＵ×Ｌ系列の伝搬信号ＳＡsym(u,l,k)を生成する。なお、伝搬信号ＳＡsym(u,l,k)および後述する各伝搬信号ＳＢsym(…)〜ＳＥsym(…)について図面ではＳＡ(…)〜ＳＥ(…)と表記する。

【0097】

窓関数演算部５２は、伝搬信号演算部５１により生成されたキャリア当りＵ×Ｌ系列の伝搬信号ＳＡsym(u,l,k)に対し、各系列に対応するドップラ周波数で回転する窓関数の乗算による信号切出しに相当する周波数領域での処理として、前記窓関数の周波数特性の畳み込み演算を行ない、キャリア当りＵ×Ｌ系列の伝搬信号ＳＢの抽出処理を行なう。

【0098】

つまり、（Ｕ・Ｌ）個の経路ごとに、各経路に対応するドップラ周波数で回転する窓関数を時間領域で掛け算することに対応する、周波数領域での補間畳み込み演算を行なう。ただし、窓関数の時間長は、各経路内のドップラ周波数が同一とみなせる複数の素波（Ｒａｙ）間の位相関係の変化が無視できる程度の時間長とする必要がある。窓関数の時間長はＯＦＤＭのＣＰ長＋有効データ長とするのが最もシンプルである。この演算の処理結果は、Ｄscを補間レートとすると、Ｕ・Ｌ×Ｋ・Ｄscとなる。このように窓関数演算部５２によりＯＦＤＭのサブキャリアデータそのものにドップラ周波数シフトを付加することで、フェージングモデルを忠実に再現することができる。

【0099】

構造的には、変調波そのもののドップラ効果による周波数ずれは、窓関数を時間領域で掛け算することに対応する周波数領域での畳み込みで使用されるフィルタ係数として、ドップラ周波数分だけずれたポイントがサンプルされたものを利用することで変調波信号の周波数シフト（周波数領域でのリサンプル）を実現する。

【0100】

ここで、窓関数演算部５２の基本的な処理について説明すると、時間領域で区間長Ｔｃ′のローカライゼーション用の窓関数ｆw_τ（τは時間軸方向のインデックス）を乗算するのと等価な処理として、ローカライゼーション用の窓関数のフーリエ変換（Ｃoe_p,i）［ｉは周波数方向の係数インデックス、ｐは１ＯＦＤＭシンボル中の窓関数の番号］を周波数領域で畳み込み処理する。

【0101】

ここでは、長さＴsym の１ＯＦＤＭシンボルをＰ分割したｐ番目（ｐ＝１，２，３，……,Ｐ）の区間長Ｔｃ′のローカライゼーション用の窓関数の乗算計算について以下のように数式化する。ただし、（Ｔｃ′＋Ｔｄ）と１／ｆsc（レイヤ周波数領域信号生成部３１の出力のサブキャリア間隔ｆscでＩＦＦＴを実施するときの時間領域のサイクル）の大小関係に応じて、処理後の周波数領域上でのサンプリング間隔が決まる。

【0102】

（ａ） （Ｔｃ′＋Ｔｄ）＞１／ｆscの場合
周波数領域でのサンプリング間隔を細かくすることで、時間長（Ｔｃ′＋Ｔｄ）の波形が時間領域でのエイリアシング（オーバーラップ）を生じないようにする。つまり、次式（１８）、（１９）に示しているように、（Ｔｃ′＋Ｔｄ）＜１／（ｆsc／Ｄsc）を満たすように、畳み込み処理結果の周波数領域でのサンプリング間隔が１／Ｄsc倍となるように補間処理をする。（元のサブキャリア同士の間にＤsc−１個のゼロを置いてから畳み込みによるフィルタ処理を実施することに相当する処理である。）

【0103】

【数8】

【0104】

ここで、式（１８）の記号Σは、ｉ＝−TapNum／２〜TapNum／２までの総和を表す。また、window(ｉ)はタップ長がＤsc・（TapNum ＋１）の窓関数でＤＦＴ(fw_τ)の周波数領域での広がりを制限するものであるとする。ＤＦＴ(ｆw_τ)は、時間スパン１／（ｆsc／Ｄsc）に渡る離散的フーリエ変換である。ｆw_τは例えば図２に示すような窓関数であり、その中心が時間０に位置しており、Ｔｃ・（１／２＋ｐ−１）だけ波形を遅らせるとｐに対応する位置に移動するような波形であるとする。

【0105】

また、式（１９）のｅ^Ｚは、時間領域でＴｃ・（１／２＋ｐ−１）だけ波形を遅らせるのに相当する周波数領域上での回転を与える項である。また、式（１９）の、Ｒｅｓａｍｐｌｅ()は、窓関数のＤＦＴ結果を(u,l)に対応するドップラ周波数分だけシフトするようにリサンプルし、全体の位相を前後のＩＤＦＴ結果が連続波形となるように(sym,p)に従って回転させる演算子であるとする。

【0106】

また、式（１８）において、ｋ’は補間処理後の周波数インデックスを表しているが、
ｋ’＝｛Ｄsc・（０−<Ｋ／２>），Ｄsc・（１−<Ｋ／２>），…，Ｄsc・（Ｋ−１−
<Ｋ／２>）｝
の位置にそれぞれｋ＝｛１，２，３，……Ｋ｝の伝搬信号ＳＡsym(u,l,k)が対応する関係にある。なお、式（１８）の括弧記号<Ａ>は、Ａを越えない最大の整数を表す（以下、同様）。

【0107】

また、式（１８）の記号％は剰余演算子であり、ｇは、ｋ’をＤscで割った時の余りである。ただし、式（１８）は、
−Ｄsc・（<Ｋ／２>＋TapNum／２）≦ｋ’≦Ｄsc・（<Ｋ／２>＋TapNum／２）
の範囲で計算する必要があり、ｉ＜０およびｉ＞Ｋの範囲でＳＡsym(u,l,i)＝０とする。

【0108】

（ｂ） （Ｔｃ′＋Ｔｄ）＜１／ｆscの場合
次式（２０）、（２１）に示すように、畳み込み処理結果のサンプリング間隔は、伝搬信号演算部５１の出力と変わらないように畳み込み処理する（補間無し）。

【0109】

【数9】

【0110】

ここで、式（２０）の記号Σは、ｉ＝−TapNum／２〜TapNum／２までの総和を表す。また、window(ｉ)はタップ長が（TapNum＋１）の窓関数でＤＦＴ(ｆw_τ)の周波数領域での広がりを制限するものであるとする。ＤＦＴ(ｆw_τ)は、時間スパン１／ｆscに渡る離散的フーリエ変換である。ｆw_τは例えば図２に示す窓関数であり、その中心が時間０に位置しており、Ｔｃ・（１／２＋ｐ−１）だけ波形を遅らせるとｐに対応する位置に移動するような波形であるとする。

【0111】

また、式（２１）のｅ^Ｚ’は、時間領域でＴｃ・（１／２＋ｐ−１）だけ波形を遅らせるのに相当する周波数領域上での回転を与える項である。また、式（２０）の、Ｒｅｓａｍｐｌｅ()は、窓関数のＤＦＴ結果を(u,l)に対応するドップラ周波数分だけシフトするようにリサンプルし、全体の位相を前後のＩＤＦＴ結果が連続波形となるように(sym,p)に従って回転させる演算子であるとする。

【0112】

また、式（２０）において、ｋ’は周波数インデックスを表しているが、
ｋ’＝｛−<Ｋ／２>，１−<Ｋ／２>，…，Ｋ−１−<Ｋ／２>｝
の位置にそれぞれｋ＝｛１，２，３，……Ｋ｝の伝搬信号ＳＡsym(u,l,k)が対応する関係にある。ただし、式（２０）は、
−（<Ｋ／２>＋TapNum／２）≦ｋ’≦（<Ｋ／２>＋TapNum／２）
の範囲で計算する必要があり、ｉ＜０およびｉ＞Ｋの範囲でＳＡsym（u,l,i)＝０とする。

【0113】

上記説明では、ドップラ周波数シフトについて述べていないが、前記したように、ドップラ周波数シフトを考慮した場合、変調波そのもののドップラ効果による周波数ずれは、窓関数を時間領域で掛け算することに対応する周波数領域での畳み込みで使用されるフィルタ係数として、図４の窓関数周波数特性として示したように、ドップラ周波数分だけずれたポイントがサンプルされたものを利用することで変調波信号の周波数シフト（周波数領域でのリサンプル）を実現している。

【0114】

パス信号加算部５３は、窓関数演算部５２で抽出されたＵ×Ｌ系列の伝搬信号ＳＢsym(u,l,k')を用いて、Ｕ個の受信アンテナ毎に異なるドップラ周波数シフトを施したＬ経路の信号を加算する演算を行う。

【0115】

パス信号加算部５３の演算結果ＳＣsym(u,l,k')は、時間領域信号生成部５４に入力される。時間領域信号生成部５４は、図６に示すように、帯域制限フィルタ５４ａおよびフーリエ逆変換演算手段５４ｂを有しており、パス信号加算部５３で得られた演算結果に対してフーリエ逆変換処理を行ない、各受信アンテナに与える時間領域の受信信号を生成する。

【0116】

帯域制限フィルタ５４ａは、次式（２２）に示すように、入力信号ＳＣsym(u,l,k')に対し、帯域制限フィルタの特性（BandFil_k'）の周波数領域における乗算を実施して、帯域制限を行なう。なお、この帯域制限処理は省略することも可能である。
ＳＤsym(u,p,k')＝ＳＣsym(u,l,k')・BandFil_k' ……（２２）

【0117】

フーリエ逆変換演算手段５４ｂは、次式（２３）のように、帯域制限された周波数領域の信号ＳＤsym(u,p,k')あるいはパス信号加算部５３の出力信号ＳＣsym(u,l,k')に対して、高速フーリエ逆変換ＩＦＦＴを行なうことで時間領域の信号ＳＥsym(u,p,τ)に変換する。

【0118】

ＳＥsym（u,p)τ＝ＩＦＦＴ{ＳＤsym(u,p,k')} ……（２３）
ただし、τ＝｛１，２，３，……，Ｎfft｝は時間のインデックスであるとする。ＮfftはＦＦＴポイント数とする。

【0119】

さらに、ｋ'が、
Ｄsc・（<Ｋ／２>＋TapNum／２）＜ｋ'＜Ｎfft−Ｄsc・（<Ｋ／２>＋TapNum／２）
にある場合、ＳＤsym(u,p,k'）＝０であり、ＳＤsym(u,p,k'）は、Ｎfftを周期として周期的とする。即ち、整数ｉについて、ＳＤsym(u,p,k'）＝ＳＤsym(u,p,k'+i・Ｎfft)が成り立つものとする。

【0120】

なお、時間領域信号生成部５４において、帯域制限フィルタ５４ａの処理をフーリエ逆変換演算手段５４ｂの処理の後に時間領域で行なうことも可能であるが、その場合、フーリエ逆変換の処理で得られた時間領域の信号に対して畳み込み演算処理を行なう必要がある。これに対し、本実施形態のように、帯域制限フィルタ５４ａをフーリエ逆変換演算手段５４ｂの前段に設けておけば、フィルタ処理を周波数領域での乗算処理で済ませることができ、畳み込み演算に比べて格段に少ない演算量で処理を実行でき、帯域制限フィルタを設ける場合であっても高速処理できる。

【0121】

シフト加算部５５は、時間領域信号生成部５４が生成した時間領域の受信信号ＳＥsym(u,p)を、前記窓関数の長さ分ずらして加算して、各受信アンテナで受信された連続する受信信号Ｓout を生成する。つまり、上式（２３）の処理結果を、図１に示したように時間Ｔｃずつずらしながら加算することで１系列分の連続した受信信号を得る。これをＵ系列分並列的に行なうことで、次式（２４）で表されるＵ系列分の連続した受信信号Ｓout を生成できる。ただし、式（２４）で、ｆs は時間領域におけるサンプリング周波数である。

【0122】

【数10】

【0123】

このように、実施形態の試験装置３０は、伝搬路の特性と変調信号との乗算処理を周波数領域において行ない、その演算結果にドップラ周波数シフトを考慮した周波数領域での窓関数演算によるリサンプル処理を行い、その演算結果から時間領域信号を生成しているので、従来方式のように、送信アンテナ毎の周波数領域の信号をフーリエ逆変換して時間領域の信号に変換してから伝搬路特性を付与する場合に比べて、フーリエ逆変換を行なう回路および伝搬路特性を生成する回路の規模を格段に小さくすることができる。

【0124】

特に、伝搬路の特性を求める演算において、ビームフォーミング特性と散乱体を考慮した送信アンテナ特性との演算を優先的に行なっているので、それ以降のドップラ周波数シフト、受信アンテナ特性およびレイリー分布を付与する演算処理では、送信アンテナ数に依存した演算処理がなくなるため、次の世代（第５世代）の通信方式で提案されているように、基地局側（送信側）がアレーアンテナのような多数のアンテナを用いるシステムを試験する場合、後続の演算処理の規模を格段に小さくでき、極めて有効である。

【0125】

例えば、Ｓ＝１２８、Ｕ＝２、レイヤ数Ｒ＝２、サブキャリア数Ｋの場合、従来方式では、時間領域信号生成のために、Ｋ個一組の信号を１２８（＝Ｓ）組分並列的にフーリエ逆変換する必要があるが、本実施形態では、Ｄsc・Ｋ個一組の信号を２（＝Ｕ）組分並列的にフーリエ逆変換すればよい。ここで、補間レートＤscが１（補間無しの場合）であれば乗算回数を、Ｕ／Ｓに縮小できる。

【0126】

また、補間する場合には、乗算回数を、
（Ｓ・２^Ｑ′・log_２２^Ｑ′）／（Ｕ・２^Ｑ・log_２２^Ｑ）
に縮小することができ、Ｄsc・Ｓ＜Ｕであれば、従来回路より少ない乗算回数で実現できる。ただし、Ｑは、（２^Ｑ−１）＜Ｋ≦２^Ｑを満たす整数、Ｑ′は、（２^Ｑ′−１）＜Ｄsc・Ｋ≦２^Ｑ′を満たす整数であり、ＩＦＦＴポイント数である。

【0127】

また、本実施形態の場合、伝搬路の特性を周波数領域に変換するためのフーリエ変換処理が必要となるが、このフーリエ変換処理では、時間領域における各パスの遅延量は、周波数領域では各パスの周波数成分の回転速度に対応するため、従来時間領域で行なっていたメモリとリサンプルフィルタの組合せにより各パスに遅延を付与するハードウエアは、フーリエ変換における回転処理に置き換えられることになり、両者のハードウエアの規模を比較して本実施形態の方が格段に有利である。

【0128】

前記実施形態は、送信アンテナから受信アンテナまでを含む伝搬路を想定して、図７のように、受信アンテナの出力（ＲＦ信号）を受けるためのＵ個のポートＰort1 〜ＰortU にケーブル２を介して試験用信号（ＲＦ信号）Ｓrx1 〜ＳrxU をそれぞれ与えて試験対象１の動作試験を行なうケーブル接続環境の例であるが、試験装置３０の実際の構成として、前記したシフト加算までの処理をアンテナで送受信する電波の無線周波数帯で行なうことは困難なので、前記したシフト加算までの処理をベースバンドで行い、得られたベースバンド信号を図示しない周波数変換部により無線周波数帯の信号（ＲＦ信号）に変換して、ケーブル２を介して試験対象１に与える構成となる。

【0129】

また、ケーブル接続環境の別の例として、試験対象１のポートＰort1 〜ＰortU が、フロントエンドから出力される中間周波数帯の信号（ＩＦ信号）を受けるポートである場合には、前記したように、シフト加算までの処理で得られたベースバンド信号を図示しない周波数変換部により中間周波数帯の信号（ＩＦ信号）に変換して、ケーブル２を介して各ポートＰort１〜ＰortＵに与えればよい。この場合、受信アンテナの特性にフロントエンドの周波数特性等も含めることが可能となり、フロントエンドの特性も含めた伝搬路を模擬した試験が行なえる。なお、試験用信号を無線周波数帯で出力する場合も中間周波数帯で出力する場合も、周波数帯がシフトしているだけで基本的には同等であり、受信アンテナで受信される受信信号とみなせる。

【0130】

また、上記のようなケーブル接続環境の試験だけでなく、例えば、図８のように、内部に受信アンテナＡrx1 〜ＡrxU を有する移動体端末本体のような試験対象１に対して電波による試験を行なう、所謂ＯＴＡ（Ｏｖｅｒ−Ｔｈｅ−Ａｉｒ）環境試験の場合には、電波暗室５内に、試験対象１とＵ’個のプローブアンテナＰrobe1 〜ＰrobeU' を配置し、試験用信号Ｓrx1 〜ＳrxU' が与えられたプローブアンテナＰrobe1 〜ＰrobeU' と試験対象１との間で電波を送受信して試験を行なう。

【0131】

このようなＯＴＡ環境試験の場合には、図９に示す構成の試験装置３０′で対応することができる。この試験装置３０′の伝搬路特性演算部４０′は、受信アンテナを含めた試験対象を試験することを考慮して、前記実施形態の試験装置３０の伝搬路特性演算部４０に対して、二つの偏波成分に対する受信アンテナの受信特性を付与するための第３の伝搬路演算部４３が省かれ、パス信号加算部５３の代わりに重み付け演算部５３′が用いられているが、その他の構成については前記実施形態と基本的に同等である。

【0132】

この構成例では、第２の伝搬路演算部４２の演算結果［Ａ２］には、二つの偏波成分それぞれについての演算結果が含まれており、これを直接受けるフーリエ変換部４４′では、第２の伝搬路演算部４２の演算結果［Ａ２］を、試験対象から見たときの到来方向（Ｌ方向に量子化）が同一とみなすことができる入射波を水平、垂直の偏波面（２面）ごとにそれぞれまとめて、伝搬路インパルス応答の時間拡がりから決まる必要十分なサンプリング間隔の周波数領域の特性に変換する。

【0133】

つまり、偏波面毎に、時間領域での素波（Ｎ・Ｍ個）の遅延に対応した回転を与え、同一の到来方向とみなせる入射波同士を加えるものであり、その処理結果は、２・Ｌ×Ｒ行列×Ｋ０（Ｋ０は周波数方向のサンプル数）となり、前記同様に、ＬおよびＫ０を小さくすることで、現実的な回路リソースに圧縮することができる。

【0134】

フーリエ変換部４４′のフーリエ変換演算手段４４ａの演算結果［Ａ４′］を受ける補間手段４４ｃにおいても、二つの偏波成分毎に周波数方向の補間処理を行なうため、その処理結果［Ａ５′］は、２・Ｌ×Ｒ行列×Ｋとなる。前記実施形態で説明したように、フーリエ変換演算手段４４ａのサンプリング間隔を大きくして補間処理を行なう方式と、フーリエ変換演算手段４４ａのサンプリング間隔を小さくして補間処理を省く方式のいずれの構成も可能であるが、前者の方式の方が回路規模を小さくする観点で有利である。

【0135】

また、伝搬信号演算部５１′の演算では、サブキャリア当りＳ×Ｒ系列の変調信号Ｓsym に対して、フーリエ変換部４４′の演算結果［Ａ５′］の行列を乗算するので、その演算で得られる伝搬信号ＳＡsym は、２つの偏波成分についてそれぞれサブキャリア当りＬ系列となり、窓関数演算部５２′の演算で得られる伝搬信号ＳＢsym も、２つの偏波成分についてそれぞれサブキャリア当りＬ系列（補間レートＤsc＝１と仮定）となる。

【0136】

この窓関数演算部５２′は、前記実施形態と同様に、周波数領域での補間畳み込み演算を行なうが、窓関数の時間長は、各経路内の到来方向が同一とみなせる複数の素波間の位相関係の変化が無視できる程度の時間長とする必要がある。

【0137】

また、重み付け演算部５３′では、各プローブアンテナについてＬ経路の信号に重みづけ加算演算を行う、即ち、窓関数演算部５２′の演算結果に対して、Ｌ個の到来方向の信号をＵ'個（水平・垂直偏波合わせて２・Ｕ'個）のＯＴＡ用プローブアンテナから出力する際の重みの乗算（偏波面ごとのＵ'×Ｌ行列乗算）を行なう。したがって、その演算結果ＳＣsym は、２つの偏波成分についてそれぞれキャリア当りＵ'系列となる。

【0138】

なお、この重み付けの演算処理は、シフト加算部５５′の処理結果に対して実施することもできる。一般に、Ｕ'＞Ｌであるため、その場合はＩＦＦＴ数は減るが、行列演算を実施するのに使うことができる乗算時間は時間領域のサンプル間隔となるため、行列乗算に必要な乗算器数がリソース使いまわし時間減少により増える効果と、ＩＦＦＴ数が減る効果のバランスで有利な方で実装すればよい。また、この重み付け演算処理をデジタル処理で行なわずに、最後にアナログ的に実施することもできる。

【0139】

さらに、時間領域信号生成部５４′の演算結果ＳＥsym およびシフト加算部５５′の演算結果Ｓout についても、２つの偏波成分についてそれぞれＵ'系列となる。

【0140】

この実施形態の試験装置３０′の場合、シフト加算部５５′で生成した２つの偏波成分についてＵ'系列の入射波を、実際の無線通信で使用するアナログの周波数帯に変換して、図８に示したように、電波暗室５内で、ＯＴＡ環境試験に用いられるプローブアンテナＰrobe1 〜ＰrobeU' を介して二つの偏波成分を移動体端末本体などの受信アンテナを有する試験対象１に送信して試験を行なうことになる。

【0141】

この試験装置３０′の場合も、前記実施形態と同様に、伝搬路の特性と変調信号との乗算処理を周波数領域において行ない、その演算結果にドップラ周波数シフトを考慮した周波数領域での窓関数演算によるリサンプル処理を行い、その演算結果から時間領域信号を生成しているので、従来方式のように、送信アンテナ毎の周波数領域の信号をフーリエ逆変換して時間領域の信号に変換してから伝搬路特性を付与する場合に比べて、フーリエ逆変換を行なう回路および伝搬路特性を生成する回路の規模を格段に小さくすることができる。

【0142】

特に、伝搬路の特性を求める演算において、ビームフォーミング特性と散乱体を考慮した送信アンテナ特性との演算を優先的に行なっているので、それ以降のドップラ周波数シフトおよびレイリー分布を付与する演算処理では、送信アンテナ数に依存した演算処理がなくなるため、次の世代（第５世代）の通信方式で提案されているように、基地局側（送信側）がアレーアンテナのような多数のアンテナを用いるシステムを試験する場合、後続の演算処理の規模を格段に小さくでき、極めて有効である。

【0143】

前記実施形態では、送信アンテナから出射される電波に二つの直交する偏波成分が含まれる場合について説明したが、送信アンテナから出射される電波を、単一の直線偏波成分としてもよい。この場合、各伝搬路演算部の演算は単一の直線偏波成分に対して行なえばよい。

【0144】

上記説明は、マルチキャリア変調方式がＯＦＤＭの場合で説明したが、他のマルチキャリア変調方式のＵＦＭＣ、ＧＦＤＭ、ＦＢＭＣ等を用いたＭＩＭＯシステムについても本発明を同様に適用できる。

【0145】

特に、第４世代Ｅｖｏｌｕｔｉｏｎおよび第５世代の携帯電話方式で利用されることが期待される３Ｄ−ＭＩＭＯ／Ｍａｓｓｉｖｅ−ＭＩＭＯにおいては、基地局の送信アンテナ数の方が、移動体端末の受信アンテナ数よりも圧倒的に多い状況となっており、本発明が非常に有効である。

【符号の説明】

【0146】

１……試験対象、３０、３０′……ＭＩＭＯ方式システムの試験装置、３１……レイヤ周波数領域信号生成部、４０、４０′……伝搬路特性演算部、４１……第１の伝搬路演算部、４２……第２の伝搬路演算部、４３……第３の伝搬路演算部、４４、４４′……フーリエ変換部、５１、５１′……伝搬信号演算部、５２、５２′……窓関数演算部、５３……パス信号加算部、５３′……重み付け演算部、５４、５４′……時間領域信号生成部、５４ａ……帯域制限フィルタ、５４ｂ……フーリエ逆変換演算手段、５５、５５′……シフト加算部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】