(57)【要約】      （修正有）【課題】伝搬パスパラメータの時間変動特性を各パス単位で生成する時変マルチパス生成装置、該装置を用いた多重波フェージングシミュレータ及びマルチパス生成方法を提供する。【解決手段】マルチパス変動を模擬する時変マルチパス生成装置が、複数の伝搬パスの生成条件を制御するパラメータ制御部；パス生成パラメータ及びデータファイルを記憶するデータ記憶部；乱数を生成して出力する乱数生成部；及び複数の伝搬パス分の複数の時変振幅関数及び時変位相関数を生成し、それらを時系列的に並べて時変遮蔽振幅関数及び時変遮蔽位相関数を伝搬パスの数だけ生成し、乱数生成部から出力された乱数を用いて初期振幅、初期位相、初期時間遅延及び初期到来方向をパスパラメータとして生成し、該複数の伝搬パスの各々の初期振幅及び初期位相に対してそれぞれ時変遮蔽振幅関数及び時変遮蔽位相関数を個別に重畳して時変伝搬パスを計算生成するパス生成部からなる。【選択図】        図３
2004254250.tif 000002
【特許請求の範囲】
【請求項１】
無線通信におけるマルチパス変動を模擬する時変マルチパス生成装置であって：
複数の伝搬パスの生成条件を制御するパラメータ制御部；
複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルを記憶するデータ記憶部；
前記パラメータ制御部から指定された乱数パラメータに基づき乱数を生成して出力する乱数生成部；及び
前記データ記憶部から得たデータファイル及びパス生成パラメータ並びに前記乱数生成部から得た複数の乱数に基づいて複数の伝搬パス分の複数の時変振幅関数及び時変位相関数を生成し、それらを時系列的に並べて時変遮蔽振幅関数及び時変遮蔽位相関数を伝搬パスの数だけ生成し、前記データ記憶部から得た初期値生成パラメータに基づき前記乱数生成部から出力された乱数を用いて初期振幅、初期位相、初期時間遅延及び初期到来方向を前記伝搬パスの各々のパスパラメータとして生成し、該複数の伝搬パスの各々の初期振幅及び初期位相に対してそれぞれ前記時変遮蔽振幅関数及び時変遮蔽位相関数を個別に重畳して複数の時変伝搬パスを計算生成するパス生成部；
からなることを特徴とする時変マルチパス生成装置。
【請求項２】
無線通信におけるマルチパス変動を模擬する時変マルチパス生成装置であって：
複数の伝搬パスの生成条件を制御するパラメータ制御部；
複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルを記憶するデータ記憶部；
前記パラメータ制御部から指定された乱数パラメータに基づき乱数を生成して出力する乱数生成部；
前記データ記憶部から得た遮蔽物パラメータ及び前記乱数生成部から得た乱数に基づき、複数の伝搬パスへの遮蔽特性としての複数の時変振幅関数及び時変位相関数を生成して出力する時変関数生成部；及び
該時変関数生成部から得た複数の伝搬パス分の複数の時変振幅関数及び時変位相関数を生成し、それらを時系列的に並べて時変遮蔽振幅関数及び時変遮蔽位相関数を伝搬パスの数だけ生成し、前記データ記憶部から得た初期値生成パラメータに基づき前記乱数生成部から出力された乱数を用いて初期振幅、初期位相、初期時間遅延及び初期到来方向を前記伝搬パスの各々のパスパラメータとして生成し、該複数の伝搬パスの各々の初期振幅及び初期位相に対してそれぞれ前記時変遮蔽振幅関数及び時変遮蔽位相関数を個別に重畳して複数の時変伝搬パスを計算生成するパス生成部；
からなることを特徴とする時変マルチパス生成装置。
【請求項３】
請求項１又は２のいずれかに記載の時変マルチパス生成装置であって：
前記パス生成部において、前記複数の伝搬パスに対して前記データ記憶部で記憶された到来波角度差に対するパラメータ間相関特性に応じた相関を有する乱数を前記乱数発生部で発生させ、遮蔽時間間隔、遮蔽振幅及び遮蔽発生時間間隔を有相関乱数値として生成して前記時変関数の定数パラメータとして与えることにより時変振幅関数及び時変位相関数を生成する；
ことを特徴とする時変マルチパス生成装置。
【請求項４】
請求項１，２又は３のいずれかに記載の時変マルチパス生成装置であって：
前記パス生成部が、生成した複数の時変伝搬パスの位相に対してパラメータ制御部から得た移動局の移動速度と各伝搬波に与えた初期到来方向からドップラー変動に対する時変位相回転を各伝搬パスに対して計算し、複数の時変伝搬パスの各時変位相特性に対して前記時変位相回転を加える；
ことを特徴とする時変マルチパス生成装置。
【請求項５】
請求項１乃至４のいずれかに記載された時変マルチパス生成装置であって：
前記データ記憶部が、前記パラメータ制御部から得た移動局に適用する複数のアンテナの指向性利得パターンファイルを記憶し；
前記パス生成部が、移動局の進行方向に合わせて少なくとも一つの前記指向性利得パターンファイルを読み出し、複数のＮ個の伝搬パスの初期到来波方向に対する指向性利得を各伝搬パスの初期振幅に対して乗じてＮ個の時変伝搬パスを計算生成し出力し、かつ複数のアンテナパターン毎に前記計算生成を行ってアンテナパターンの数Ｍ個に対応する総計Ｎ×Ｍ個の時変伝搬パスを計算生成する；
ことを特徴とする時変マルチパス生成装置。
【請求項６】
請求項１乃至５のいずれかに記載された時変マルチパス生成装置；
該マルチパス生成装置から得た複数の時変伝搬パスの複素振幅を実部と虚部に分けてアナログ出力するパス出力部；及び
少なくとも１つのデジタル信号入力端子と，少なくとも１つのデジタル信号出力端子と，前記デジタル信号入力端子から入力される複数の信号から直交入力信号（Ｑ成分）を作成し，前記複数の信号の入力信号（Ｉ成分）と該直交入力信号（Ｑ成分）とを、それぞれＮ個の複数の縦続接続した遅延素子と乗算器からなるトランスバーサル回路に入力し、前記パス出力部から出力される各伝搬パスの初期時間遅延が前記トランスバーサル回路の遅延素子通過後の各段の合計遅延と等しくなる乗算器にそれぞれ入力し、前記遅延素子で遅延させた各デジタル信号にそれぞれのパス変動を乗算し、各遅延時間毎の乗算結果を加算して生成されるＩ成分及びＱ成分の信号を再び合成する信号合成部を複数個（Ｍ個）有して前記パス出力部から出力されるアンテナ指向性利得パターン毎の時変伝搬パス出力（Ｎ×Ｍ個）に対して前記入力信号との信号合成を行い，該デジタル信号出力端子に出力する信号合成部；
から構成されることを特徴とする多重波フェージングシミュレータ。
【請求項７】
無線通信におけるマルチパス変動を模擬するための時変マルチパスを生成する方法であって：
複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルをデータ記憶部に記憶する段階（Ｓ１）；
パラメータ制御部において、前記データ記憶部からパス生成パラメータファイルを読み出す段階（Ｓ１）；
パラメータ制御部において、Ｍ個のアンテナ指向性利得パターンファイル等のアンテナに関するデータファイルをデータ記憶部から読み出す段階（Ｓ２）；
乱数生成部において、生成波のパスパラメータ初期値生成条件を満たす乱数を生成する段階（Ｓ３）；
パス生成部において、前記乱数に基づいてＮ個のパスパラメータの初期振幅値等の初期条件を設定する段階（Ｓ３）；
アンテナ指向性利得パターンファイルから得た到来波方向に対するアンテナ指向性利得を、Ｎ個の初期振幅値に乗算し、Ｎ×Ｍ個の伝搬パスパラメータを設定する段階（Ｓ３）；
パラメータ制御部において、データ記憶部から時変関数の定数パラメータ生成用条件ファイルを読み出す段階（Ｓ５）；
パス生成部５において、乱数生成部で生成した乱数に基づき時変関数の定数パラメータを決定する段階（Ｓ６）；
パス生成部５において、決定した時変関数の定数パラメータに基づき各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性を生成し、それらを各波の振幅及び位相に対して演算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する段階（Ｓ７）；及び
生成した遮蔽変動特性に基づいて、時変マルチパスを生成する段階
から構成される方法。
【請求項８】
無線通信におけるマルチパス変動を模擬するための時変マルチパスを生成する方法であって：
複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルをデータ記憶部に記憶する段階（Ｓ１）；
パラメータ制御部において、前記データ記憶部からパス生成パラメータファイルを読み出す段階（Ｓ１）；
パラメータ制御部において、Ｍ個のアンテナ指向性利得パターンファイル等のアンテナに関するデータファイルをデータ記憶部から読み出す段階（Ｓ２）；
乱数生成部において、生成波のパスパラメータ初期値生成条件を満たす乱数を生成する段階（Ｓ３）；
パス生成部において、前記乱数に基づいてＮ個のパスパラメータの初期振幅値等の初期条件を設定する段階（Ｓ３）；
アンテナ指向性利得パターンファイルから得た到来波方向に対するアンテナ指向性利得を、Ｎ個の初期振幅値に乗算し、Ｎ×Ｍ個の伝搬パスパラメータを設定する段階（Ｓ３）；
パラメータ制御部において、パス生成パラメータファイルに基づき計算モデルを指定する段階（Ｓ１１）；
乱数生成部において、遮蔽物のパラメータに対する各乱数を生成する段階（Ｓ１２）；
前記の乱数に基づき、時変関数生成部において、指定された計算モデルのＮ個の時変関数を生成する段階（Ｓ１２）；
パス生成部５において、生成した時変関数に基づき各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性を生成し、それらを各波の振幅及び位相に対して演算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する段階（Ｓ１３）；及び
生成した遮蔽変動特性に基づいて、時変マルチパスを生成する段階
から構成される方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、一般に移動通信システム等の無線通信システムにおいて、マルチパス変動を模擬する時変マルチパス生成装置、該装置を用いた多重波フェージングシミュレータ及びマルチパス生成方法に関し、特に入力パラメータに対応する各伝搬パスのパラメータ（到来方向特性、位相特性、遅延特性、レベル特性）に時変パス変動を付加することにより時変マルチパスを生成する技法を用いた時変マルチパス生成装置、該装置を用いた多重波フェージングシミュレータ及びマルチパス生成方法に関する。
【０００２】
【従来の技術】
図１に、基地局アンテナから電波を送信し、建物等で反射や回折を経て到達する電波を移動局が受信した場合の概念図を示す。図１の右側には、遅延プロファイルの例を示す。図１に示した遅延プロファイルにおいて、横軸は移動局に到来する電波の伝搬遅延時間（以下、時間遅延と呼ぶ）であり、縦軸は受信電力である。時間遅延、受信電力は絶対的なものであっても、相対的なものであっても良い。
【０００３】
遅延プロファイルの▲１▼，▲２▼，▲３▼，▲４▼，▲５▼の電波は平面波（パス）と見なせる。時間遅延の最も小さいパスは基地局から最短距離で到達したパス▲１▼であり、時間遅延が▲１▼より大きなパスは遠方の建物や山岳等で反射、回折を経て到達したパスである。時間遅延が大きなパスほど伝搬距離が大きくなるために伝搬損失が増加し、反射、回折回数も増加するため受信電力は小さくなる傾向がある。
図２は、移動局が道路に沿って移動する度に測定した遅延プロファイルの測定例を示すグラフである（非特許文献１）。移動局の走行に伴ってパスの時間遅延と受信電力が刻々と変化している様子が分かる。
【０００４】
移動通信システムにおいてある方式を設計する際にその方式がパス変動に耐えられるかどうかを判断するため、或いはその方式の中のパラメータの最適値を決定するため等に、基地局と移動局との間の空間における伝送品質をシミュレートするシミュレータが必要である。移動通信における伝送品質をシミュレーション装置で評価する場合等には、フェージングシミュレータを用いて、送信機から入力された変調信号にフェージングをかけて受信機に対して出力する。その際、フェージングシミュレータは、例えばＩＴＵ（Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｔｅｌｅｃｏｍｍｕｎｉｃａｔｉｏｎ  Ｕｎｉｏｎ）等で推奨されているプロファイルモデル（Ｖｅｈｉｃｕｌａｒ−Ｂモデル（非特許文献２）に従いユーザが入力した時間遅延を有する伝搬波に変調信号を重畳した後、仲上−ライスフェージングあるいはレイリーフェージングをかけるか、又はレイトレースにより遅延プロファイルを計算した上で各時間遅延の伝搬波に変調信号を重畳した後、同様のフェージングをかけていた。
【０００５】
伝送速度の低い従来の移動通信システムでは、受信レベル特性（例：レイリーフェージングモデル；非特許文献３）や遅延特性（例：等価伝送路モデル；非特許文献４）を考慮すれば十分であったため、実環境を正確に模擬したパスの伝搬特性を明らかにする必要はなかった。従って、複数のパスをベクトル合成したパスの変動として明らかにされている仲上−ライス分布やレイリー分布に従うレベル変動を、パス間での相関無しに与えるフェージングシミュレータが考えられていた。
【０００６】
【非特許文献１】
進士昌明著  “無線通信の電波伝搬”  １９９２年  第２１０頁
【０００７】
【非特許文献２】
ＩＴＵ−Ｒ  Ｍ．１２２５，  ”Ｇｕｉｄｅｌｉｎｅｓ  ｆｏｒ  ｅｖａｌｕａｔｉｏｎｓ  ｏｆ  ｒａｄｉｏ  ｔｒａｎｓｍｉｓｓｉｏｎ  ｔｅｃｈｎｏｌｏｇｉｅｓ  ｆｏｒ  ＩＭＴ−２０００，”  １９９７
【０００８】
【非特許文献３】
Ｊａｋｅｓ  Ｗ．  Ｃ．  Ｊｒ．，  ｅｄ．  ：“Ｍｉｃｒｏｗａｖｅ  ｍｏｂｉｌｅ  ｃｏｍｍｕｎｉｃａｔｉｏｎｓ”，  Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ，  Ｉｎｃ．，  Ｎｅｗ  Ｙｏｒｋ  ，  １９７４
【０００９】
【非特許文献４】
Ｙ．  Ｋａｒａｓａｗａ，  Ｔ．  Ｋｕｒｏｄａ，  ａｎｄ  Ｈ．  Ｉｗａｉ：  “Ｔｈｅ  ｅｑｕｉｖａｌｅｎｔ  ｔｒａｎｓｍｉｓｓｉｏｎ−ｐａｔｈ  ｍｏｄｅｌ，  −Ａ  ｔｏｏｌ  ｆｏｒ  ａｎａｌｙｚｉｎｇ  ｅｒｒｏｒ  ｆｌｏｏｒ  ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ  ｄｕｅ  ｔｏ  ｉｎｔｅｒ−ｓｙｍｂｏｌ  ｉｎｔｅｒｆｅｒｅｎｃｅ  ｉｎ  Ｎａｋａｇａｍｉ−Ｒｉｃｅ  ｆａｄｉｎｇ  ｅｎｖｉｒｏｎｍｅｎｔｓ，”  ＩＥＥＥ  Ｔｒａｎｓ．  Ｖｅｈ．  Ｔｅｃｈｎｏｌ．  ４６  ［１］  ｐｐ．１９４−２０２，  １９９７
【００１０】
【非特許文献５】
：山田，富里，松本：「実伝搬データを用いた時空等化器のシステム評価」，信学技報Ａ・Ｐ２０００−９６  （２０００−１０）
【００１１】
【非特許文献６】
Ｓ．  Ｉｃｈｉｔｓｕｂｏ，  Ｋ．  Ｔｓｕｎｅｋａｗａ，  ａｎｄ  Ｙ．  Ｅｂｉｎｅ：  “Ｍｕｌｔｉｐａｔｈ  Ｐｒｏｐａｇａｔｉｏｎ  Ｍｏｄｅｌ  ｏｆ  Ｓｐａｔｉｏ−Ｔｅｍｐｏｒａｌ  Ｄｉｓｐｅｒｓｉｏｎ  Ｏｂｓｅｒｖｅｄ  ａｔ  Ｂａｓｅ  Ｓｔａｔｉｏｎ  ｉｎ  Ｕｒｂａｎ  Ａｒｅａｓ”，  ＩＥＥＥ  Ｊｏｕｒｎａｌ  ｏｎ  ｓｅｌｅｃｔｅｄ  ａｒｅａｓ  ｉｎ  ｃｏｍｍｕｎｉｃａｔｉｏｎｓ，  Ｖｏｌ．２０，  Ｎｏ．６，  Ａｕｇｕｓｔ  ２００２
【００１２】
【発明が解決しようとする課題】
例えば第４世代移動通信方式（広帯域伝送方式）においては、遅延時間軸上で不要なパスを除去する適応等化技術や空間軸（電波到来方向）上で不要なパスを分離するアダプティブアレーアンテナ（ＡＡＡ）が必要不可欠であると考えられており、その実現に向けた研究開発が精力的に進められている（非特許文献５）。また、伝送速度を高速化（ビットレートを大きく）すればそれに応じて到来波を複数の伝搬遅延波（単一のパス、又は複数のパスを合成したもの）に分離することができるため、伝送速度の低い従来の伝送方式で問題であったフェージングの影響を軽減でき、高品質な通信が実現できる。そのため、高精度な伝送特性評価を行うためには送信局から送信された無線波が実伝搬路内を伝搬して受信局まで達するまでの複数の伝搬パスのそれぞれについて伝搬特性を再現できることが重要となる。
【００１３】
パスの伝搬特性には、例えば伝搬遅延特性や到来方向特性があり、その形状の分析やモデル化の検討が盛んに進められている（非特許文献６）。例えば伝搬遅延特性に関しては、時間遅延に対して指数関数的に減衰する指数関数モデルでおおよそ近似できることが知られている（非特許文献３）。しかしながら、そのようなモデルに基づいた伝搬パスに対する時間変動特性の与え方は、単なる仲上−ライスフェージングやレイリーフェージング変動を重畳するという方法でしかなく、これらは複数のパスの合成結果として現れる変動特性を与えているに過ぎなかった。
【００１４】
さらに、適応等化器やＡＡＡを組み合わせた技術を用いる通信方式の評価を行うためには、パスの到来方向特性、遅延時間特性、及び受信レベル特性の全てを考慮したパスモデルが不可欠である。
【００１５】
本来の各パス単位のレベル変動の要因として、周囲に存在する自動車、人体といった地物の移動によるパスの遮蔽を挙げる事ができる。さらに、移動局が移動する場合には、建物等によるパス遮蔽がそれに加わると考えることができる。これらの遮蔽によるレベル変動を各パスの変動として考慮すべきであるが、従来のフェージングシミュレータではそのような評価装置とはなっていなかった。
【００１６】
また、受信局への到来方向がほぼ同じである伝搬波については、同一の遮蔽物によるパス遮蔽変動を受けている可能性が高く、従ってその変動特性には強い相関があるはずであるが、そのような空間的なパスの相関特性を考慮できるフェージングシミュレータは未だ提供されていない。
【００１７】
そこで、本発明の一目的は、入力されたシミュレーション条件に対応する伝搬パスパラメータ（レベル特性、位相特性、時間遅延特性、到来方向特性）の時間変動特性（時変マルチパス変動）を、各パス単位で生成し、かつ空間的なパスの相関特性を考慮した時変パス特性として生成する時変マルチパス生成装置、該装置を用いた多重波フェージングシミュレータ及びマルチパス生成方法を提供することにある
【００１８】
【課題を解決するための手段】
上記目的を達成するための本発明の一特徴に従った、無線通信におけるマルチパス変動を模擬する時変マルチパス生成装置は、  複数の伝搬パスの生成条件を制御するパラメータ制御部；  複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルを記憶するデータ記憶部；  パラメータ制御部から指定された乱数パラメータに基づき乱数を生成して出力する乱数生成部；及び  データ記憶部から得たデータファイル及びパス生成パラメータ並びに乱数生成部から得た複数の乱数に基づいて複数の伝搬パス分の複数の時変振幅関数及び時変位相関数を生成し、それらを時系列的に並べて時変遮蔽振幅関数及び時変遮蔽位相関数を伝搬パスの数だけ生成し、データ記憶部から得た初期値生成パラメータに基づき乱数生成部から出力された乱数を用いて初期振幅、初期位相、初期時間遅延及び初期到来方向を伝搬パスの各々のパスパラメータとして生成し、該複数の伝搬パスの各々の初期振幅及び初期位相に対してそれぞれ時変遮蔽振幅関数及び時変遮蔽位相関数を個別に重畳して複数の時変伝搬パスを計算生成するパス生成部；
からなることを特徴とする。
【００１９】
本発明の他の特徴に従った、無線通信におけるマルチパス変動を模擬する時変マルチパス生成装置は、  複数の伝搬パスの生成条件を制御するパラメータ制御部；  複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルを記憶するデータ記憶部；  パラメータ制御部から指定された乱数パラメータに基づき乱数を生成して出力する乱数生成部；  データ記憶部から得た遮蔽物パラメータ及び乱数生成部から得た乱数に基づき、複数の伝搬パスへの遮蔽特性としての複数の時変振幅関数及び時変位相関数を生成して出力する時変関数生成部；及び  該時変関数生成部から得た複数の伝搬パス分の複数の時変振幅関数及び時変位相関数を生成し、それらを時系列的に並べて時変遮蔽振幅関数及び時変遮蔽位相関数を伝搬パスの数だけ生成し、データ記憶部から得た初期値生成パラメータに基づき乱数生成部から出力された乱数を用いて初期振幅、初期位相、初期時間遅延及び初期到来方向を伝搬パスの各々のパスパラメータとして生成し、該複数の伝搬パスの各々の初期振幅及び初期位相に対してそれぞれ時変遮蔽振幅関数及び時変遮蔽位相関数を個別に重畳して複数の時変伝搬パスを計算生成するパス生成部；  からなることを特徴とする。
【００２０】
本発明の他の特徴に従った、多重波フェージングシミュレータは、上記のような時変マルチパス生成装置；  該マルチパス生成装置から得た複数の時変伝搬パスの複素振幅を実部と虚部に分けてアナログ出力するパス出力部；及び  少なくとも１つのデジタル信号入力端子と，少なくとも１つのデジタル信号出力端子と，デジタル信号入力端子から入力される複数の信号から直交入力信号（Ｑ成分）を作成し，複数の信号の入力信号（Ｉ成分）と該直交入力信号（Ｑ成分）とを、それぞれＮ個の複数の縦続接続した遅延素子と乗算器からなるトランスバーサル回路に入力し、パス出力部から出力される各伝搬パスの初期時間遅延がトランスバーサル回路の遅延素子通過後の各段の合計遅延と等しくなる乗算器にそれぞれ入力し、遅延素子で遅延させた各デジタル信号にそれぞれのパス変動を乗算し、各遅延時間毎の乗算結果を加算して生成されるＩ成分及びＱ成分の信号を再び合成する信号合成部を複数個（Ｍ個）有してパス出力部から出力されるアンテナ指向性利得パターン毎の時変伝搬パス出力（Ｎ×Ｍ個）に対して入力信号との信号合成を行い，該デジタル信号出力端子に出力する信号合成部；
から構成されることを特徴とする。
【００２１】
本発明の他の特徴に従った、無線通信におけるマルチパス変動を模擬するための時変マルチパスを生成する方法は、  複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルをデータ記憶部に記憶する段階（Ｓ１）；
パラメータ制御部において、データ記憶部からパス生成パラメータファイルを読み出す段階（Ｓ１）；  パラメータ制御部において、Ｍ個のアンテナ指向性利得パターンファイル等のアンテナに関するデータファイルをデータ記憶部から読み出す段階（Ｓ２）；  乱数生成部において、生成波のパスパラメータ初期値生成条件を満たす乱数を生成する段階（Ｓ３）；  パス生成部において、乱数に基づいてＮ個のパスパラメータの初期振幅値等の初期条件を設定する段階（Ｓ３）；
アンテナ指向性利得パターンファイルから得た到来波方向に対するアンテナ指向性利得を、Ｎ個の初期振幅値に乗算し、Ｎ×Ｍ個の伝搬パスパラメータを設定する段階（Ｓ３）；  パラメータ制御部において、データ記憶部から時変関数の定数パラメータ生成用条件ファイルを読み出す段階（Ｓ５）；  パス生成部５において、乱数生成部で生成した乱数に基づき時変関数の定数パラメータを決定する段階（Ｓ６）；  パス生成部５において、決定した時変関数の定数パラメータに基づき各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性を生成し、それらを各波の振幅及び位相に対して演算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する段階（Ｓ７）；  及び生成した遮蔽変動特性に基づいて、時変マルチパスを生成する段階から構成される。
【００２２】
本発明のさらに他の特徴に従った、無線通信におけるマルチパス変動を模擬するための時変マルチパスを生成する方法は、  複数の伝搬パスの生成に必要な各種パス生成パラメータ及びデータファイルをデータ記憶部に記憶する段階（Ｓ１）；  パラメータ制御部において、データ記憶部からパス生成パラメータファイルを読み出す段階（Ｓ１）；  パラメータ制御部において、Ｍ個のアンテナ指向性利得パターンファイル等のアンテナに関するデータファイルをデータ記憶部から読み出す段階（Ｓ２）；  乱数生成部において、生成波のパスパラメータ初期値生成条件を満たす乱数を生成する段階（Ｓ３）；  パス生成部において、乱数に基づいてＮ個のパスパラメータの初期振幅値等の初期条件を設定する段階（Ｓ３）；  アンテナ指向性利得パターンファイルから得た到来波方向に対するアンテナ指向性利得を、Ｎ個の初期振幅値に乗算し、Ｎ×Ｍ個の伝搬パスパラメータを設定する段階（Ｓ３）；  パラメータ制御部において、パス生成パラメータファイルに基づき計算モデルを指定する段階（Ｓ１１）；  乱数生成部において、遮蔽物のパラメータに対する各乱数を生成する段階（Ｓ１２）；  該乱数に基づき、時変関数生成部において、指定された計算モデルのＮ個の時変関数を生成する段階（Ｓ１２）；  パス生成部５において、生成した時変関数に基づき各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性を生成し、それらを各波の振幅及び位相に対して演算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する段階（Ｓ１３）；  及び生成した遮蔽変動特性に基づいて、時変マルチパスを生成する段階  から構成される。
【００２３】
【作用】
本発明は、各パスの合成結果としての変動特性である仲上−ライスフェージングやレイリーフェージングなどの変動を各パスの変動として模擬するのではなく、現実の伝搬路内で実際にパス変動の要因と考えられる車両や人体などの移動物体及び建物や構造物などの地物による遮蔽によって生じるパス変動を各パスの変動として与えるもので、そのようなパス変動を考慮することから更に空間的なパスの変動相関をも考慮できるパス変動を生成して多重波フェージングを模擬することができる。
【００２４】
【発明の実施の形態】
以下に図面を参照しながら、本発明の実施の形態について説明する。
【００２５】
図３は、本発明の実施形態に係る時変マルチパス生成装置の構成例を示す図である。本明細書において「時変」とは、時間的変動の意味であり、例えばあるパスの初期値として与えた方向、位相、遅延等に対して時間的変動を与えるということを意味する。本実施形態に係る時変マルチパス生成装置９は、以下の処理部から構成される。すなわち、入力デバイス１、パラメータ入力・制御部２、データ記憶部３、乱数生成部４、パス生成部５、時変関数生成部６、出力データのバッファメモリ部７、時変伝搬パスの出力インターフェース部８から構成される。また、図４に本発明の多重波シミュレータの実施の形態を示す。多重波シミュレータ１４は、時変伝搬パスの出力インターフェース部８、時変マルチパス生成装置９、パス出力部１０、信号合成部１１、デジタル信号入力端子１２、デジタル信号出力端子１３から構成される。以下に各部の説明を行う。
【００２６】
入力デバイス１は、キーボード、マウス、フロッピー（登録商標）ディスク、メモリーカード、ＣＤ−ＲＯＭ、ＤＶＤ、ＭＯなどから成り、パラメータ入力・制御部２において表示されるパス生成に必要な各種パラメータの入力、データ記憶部３に格納すべきデータファイル類の入出力を行うデバイスである。
【００２７】
パラメータ入力・制御部２は、パス生成パラメータ、各種データファイルに関する情報を表示すると共に、入力デバイス１を用いて必要なパラメータの入力を行う。また、パス生成に必要な時変関数（時間を変数とした関数）やアンテナ指向性利得パターンなどのデータファイルを必要に応じて予めパラメータ入力・制御部２を介して入力され、データ記憶部３に格納する。データ記憶部３に予め必要なデータが記憶されている場合には、パラメータ入力・制御部３を介してデータを入力する必要はない。その場合には、単にパラメータ制御部として機能しても良い。更にパス生成のためのパラメータ指定等を行い、データ記憶部３、パス生成部５及び時変関数生成部６に対して必要な情報のやり取りを指示しながらパス生成の開始から終了までの制御を実行する。
【００２８】
データ記憶部３には、パス生成パラメータファイル、時変関数ファイル、アンテナの指向性利得パターンファイル、到来波方向密度関数ファイル、分布関数ファイル等が保存される。
【００２９】
パス生成パラメータファイルは、パラメータ入力・制御部２から入力された伝搬パスの周波数Ｆ［ＧＨｚ］、生成波数Ｎ、出力アンテナブランチ数Ｍ、パスパラメータ初期値生成条件、遅延プロファイル関数指定、平均遅延スプレッドＳＤ、時変関数タイプ指定値、定数パラメータ生成用条件ファイル、遮蔽時間間隔ＴＡ［ｓｅｃ］の平均値＜ＴＡ＞及び標準偏差σＡ、遮蔽振幅ＬＢ［ｄＢ］の平均値＜ＬＢ＞及び標準偏差σＢ、遮蔽発生時間間隔ＴＣ［ｓｅｃ］の平均値＜ＴＣ＞及び標準偏差σＣ、落込み係数Ｃｄの平均値＜Ｃｄ＞及び標準偏差σＣｄ、立ち上がり係数Ｃｒの平均値＜Ｃｒ＞及び標準偏差σＣｒ、遮蔽時標準偏差σＳ（ＬＢ）、到来波角度差に対するパラメータ間相関特性ファイルρＡ（Ω）、ρＢ（Ω）、ρＣ（Ω）、遮蔽変動計算モデルファイル、遮蔽物寸法（Ｗｈ［ｍ］、Ｄｈ［ｍ］、Ｈｈ［ｍ］）の確率密度関数ファイル、伝搬パスへの遮蔽物体の進入角度θｓ［ｄｅｇｒｅｅ］  の確率密度関数ファイル、送受信点間距離（最短パスの伝搬距離）Ｌ０、遮蔽物の進入位置Ｌｄ［ｍ］  の確率密度関数ファイル、遮蔽物の移動速度Ｖｓ［ｍ／ｓ］  の確率密度関数ファイル、移動局の速度Ｖ［ｍ／ｓ］、データサンプリング間隔Ｄｔ［ｓｅｃ］等の情報から成る。
【００３０】
時変関数ファイルには、パラメータ入力部２から入力されるパスの振幅及び位相に対する複数の遮蔽変動の関数Ａｊ（ｔ）、φｊ（ｔ）  が格納されている。
【００３１】
アンテナの指向性利得パターンファイルには、シミュレーションにおいて使用する各種アンテナの指向性利得パターン（Ｖ偏波パターン及びＨ偏波パターン）が複数個格納され、それらは計算により作成されたものであっても測定により得られたものでもよく、同じフォーマットで記述されていないものであれば前記パラメータ入力・制御部２を介して同一フォーマットに変換したものが格納される。
【００３２】
到来波方向密度関数ファイルは、到来波方向に対するパスの生起確率を規定した確率密度関数が複数個格納されており、それらはパス生成時に指定される環境条件に対してそれぞれ異なる確率密度関数が用いられる。
【００３３】
分布関数ファイルには、各種環境条件に対する乱数生起用確率密度関数（指数分布関数、べき乗分布関数、一様分布関数等）が格納されている。
【００３４】
乱数生成部４は、前記パラメータ入力・制御部２から指定された乱数生成パラメータに基づき各種乱数を要求された個数だけ生成し、パス生成部５及び時変関数生成部６に必要な乱数生成結果をわたす。前記パラメータ入力・制御部２からの生成要求が継続する場合には、前記パラメータ入力・制御部２から停止要求が出されるまで前記乱数の生成を繰り返し行い、パス生成部５及び時変関数生成部６に、生成した必要な乱数を繰り返しわたす。
【００３５】
パス生成部５は、前記パラメータ入力・制御部２からの生成要求により、パス生成の開始及び停止を行う。パス生成部５はまず、前記パラメータ入力・制御部２から指定されるパス生成パラメータファイルの記述内容に基づき、伝搬パス数（Ｎ個とする）分のパスパラメータ（振幅、位相、時間遅延、到来方向）の初期値を前記乱数生成部４からの乱数生成結果に基づき設定する。パス生成部５で生成する時変遮蔽変動を前記データ記憶部３に格納した時変振幅関数及び時変位相関数を適用して生成するか、もしくは前記時変関数生成部６において新たに計算・生成される時変振幅関数及び時変位相関数を適用して生成するかは、前記パス生成パラメータファイルの記述内容に基づき選択される。
【００３６】
データ記憶部３に格納した時変関数を適用する場合、時変関数は時変関数ファイルを読み込んで決まる。該時変関数が有する複数の定数パラメータは前記パス生成パラメータファイルの記述内容に基づき前記乱数生成部４に要求した乱数出力結果より生成し、Ｎ個の伝搬パス分の異なる時変関数を生成する。
【００３７】
時変関数生成部６において新たに計算・生成される時変振幅関数及び時変位相関数を適用する場合、パス生成パラメータファイルの記述内容に基づく遮蔽計算モデルに従って時変関数生成部６が計算・出力する遮蔽変動特性（振幅、位相）の計算結果を時変関数として適用すると共に、前記パス生成パラメータファイルの記述内容に基づき前記乱数生成部４に要求した別の乱数出力結果より遮蔽条件を変えてＮ個の伝搬パス分の異なる時変関数を生成する。
【００３８】
また、パス生成部５は、パラメータ入力・制御部２からの停止要求を受けるまで、前記手順を繰り返して時変関数をＮ個分ずつ生成し、それぞれの伝搬パスに対応する時変関数を時系列的に継続接続して連続する時変関数を生成する。このように生成したＮ個の連続時変振幅関数及びＮ個の連続時変位相関数を、それぞれＮ個の伝搬パスの初期振幅から減算しかつ初期位相に加算することにより、Ｎ個の伝搬パスの時変振幅特性及び時変位相特性の生成が実現される。
【００３９】
図５に示すように、パス生成部５は、移動局の速度及び初期到来波方向（移動局の進行方向を基準に規定される）の情報から、各波の時間的位相回転Δφｉ（ｔ）を次式で計算し、これをそれぞれのパスに対する初期位相に加算することにより、Ｎ個の伝搬パスの時変位相特性にドップラー位相回転の効果を与えることができる。
【００４０】
【数１】
2004254250.tif 000003
但し、Ｖ  ［ｍ／ｓｅｃ］は移動局の速度、ｔ［ｓｅｃ］は時間、Ωｉ＝（θｉ，φｉ）はｉ番目の到来波の初期到来方向、λは波長である。
【００４１】
更に、パス生成部５は、パス生成パラメータファイルの記述内容に基づき、指定されたアンテナの指向性利得パターンファイルのデータをパラメータ入力・制御部２を介してデータ記憶部３から読み込み、各パスの初期到来方向に応じた指向性利得を初期振幅に乗じて前記パス生成処理を行い、Ｍ個（Ｍ≧１）の複数のアンテナパターン毎にパス生成処理を実行してＭ個のグループに分けたＮ個の時変伝搬パスを出力することもできる。
【００４２】
時変関数生成部６は、前記パラメータ入力・制御部２が指定したパス生成パラメータファイルの記述内容に基づき、指定された計算モデルとモデル計算に用いる各種計算条件を前記乱数生成部４の出力結果を用いて変更しながらＮ個の伝搬パス分の時変遮蔽特性（振幅，位相）を計算し、パス生成部５にその結果をわたす。パス生成部５あるいはパラメータ入力・制御部２の要求毎に必要な個数の時変遮蔽関数を計算して出力をつづけ、計算停止要求が来るまで実行を続ける。
【００４３】
出力データのバッファメモリ部７は、パス生成部が生成した時変伝搬パスの出力データを一時的に格納し、パラメータ入力・制御部２からの停止要求が来ない場合は、該格納したデータが予め設定されたデータ容量の閾値を超えた時点で出力インターフェース部８にデータを転送し、パラメータ入力・制御部２からの停止要求が来た時点で格納されているデータ及びパス出力部５から転送されてくるデータを順次、全て出力インターフェース部８に転送する。
【００４４】
図４を参照すると、時変伝搬パスの出力インターフェース部８は、パス生成部５で生成したＮ個の時変伝搬パス特性（時変振幅、時変位相、初期時間遅延、初期到来方向）の結果を図４に示すパス出力部９のような外部装置に出力するインターフェース部であり、生成データは一旦出力データのバッファメモリ部７に格納された後に本出力インターフェース部８に転送される。該時変特性は、Ｎ×Ｍ個の時系列的変動値として本インターフェース部８を介して外部装置に出力される。
【００４５】
パス出力部１０は、前記パス生成部５において計算生成され、前記出力インターフェース部８を介して得られるＮ個の複数の時変伝搬パスを時変振幅及び時変位相から複素振幅の実部と虚部の変動へと変換すると共に、デジタル／アナログ変換を行って２Ｎ×Ｍ個のアナログ信号を出力する。
【００４６】
信号合成部１１は、前記パス出力部１０からのＮ個の複素振幅の実部時変信号ｗＲ（１１）  〜ｗＲ（ＭＮ）とＮ個の複素振幅の虚部時変信号ｗＩ（１１）〜ｗＩ（ＭＮ）を一つのグループとするＭ個のグループの時変伝搬パス信号を受信入力し、デジタル入力信号端子１２において入力されたＩＦ又はＲＦ信号ｕ（ｔ）をＭ個に分配し、分配したそれぞれのデジタル信号毎に位相をπ／２だけ変化させた直交信号（Ｑ成分）を作成し、前記デジタル入力信号（Ｉ成分：ｕＩ（ｔ））とＱ成分：ｕＱ（ｔ）をそれぞれＮ個の複数の遅延素子を縦続接続したトランスバーサル型回路（図１６参照）に入力し、パス出力部１０からのそれぞれの伝搬パスの時間遅延に対応するＩ成分のＮ個の遅延信号にＮ個の複素振幅の実部時変信号をそれぞれ乗算し、同様にＱ成分のＮ個の遅延信号にＮ個の複素振幅の虚部時変信号をそれぞれ乗算し、それらを加算合成することにより、Ｍ個の合成信号ｙ（Ｍ）（ｔ）を生成し、デジタル信号出力端子１３から出力する。入力信号と出力信号の関係を次式に示す。
【００４７】
【数２】
2004254250.tif 000004
ここに上式においては、キャリア成分表記を省略している。
【００４８】
図６は、本実施形態に係る時変マルチパス生成装置９の動作を示すフローチャートである。以下にこのフローチャートの各ステップＳ１−Ｓ１５に沿って本発明の動作を説明する。尚、本発明の時変マルチパス波の生成を実行する前に、パス生成パラメータファイル、時変関数ファイル、アンテナの指向性利得パターンファイル、到来波方向密度関数ファイル等のデータファイルがオペレータによって事前にデータ記憶部３に格納・保存されているものとする。
【００４９】
Ｓ１：パラメータ入力・制御部２は、パス生成パラメータファイルをデータ記憶部３から読み出し、生成波数Ｎ、出力アンテナブランチ数Ｍに対応する生成データ格納用のメモリ領域を出力データ用バッファメモリ部７に確保する。パラメータ入力・制御２は、パス生成パラメータファイルをパス生成部５と共有しあるいはパス生成部５へも供給して、ステップ２（Ｓ２はステップ２を表す。）へと進む。
【００５０】
Ｓ２：パラメータ入力・制御部２は、出力アンテナブランチ数Ｍに対応するＭ個のアンテナ指向性利得パターンファイルと到来波方向密度関数ファイルをデータ記憶部３から読み出して、ステップ３へと進む。
【００５１】
Ｓ３：パス生成部５が、生成波のパスパラメータ初期値生成条件を満たす乱数生成・出力を乱数生成部４に対して要求し、乱数生成部４により生成した乱数を得る。この乱数によりＮ個のパスパラメータの初期条件（初期振幅、初期位相、初期時間遅延、初期到来方向）を設定する。
【００５２】
初期振幅及び初期時間遅延は、例えばパス生成パラメータファイル内で指定される遅延プロファイル関数指定が指数分布関数（係数は指定値）であれば、生成波の遅延時間間隔をΔＴ、生成波数Ｎ、平均遅延スプレッド値ＳＤによりＮ波の初期振幅及び初期時間遅延がそれぞれ決定される。また初期位相は、例えば［０，２π］の範囲で一様分布であるとすれば、一様乱数をＮ個生成することにより決定され初期到来方向は到来波方向密度関数ファイルに基づき、この到来波方向密度関数に従う乱数をＮ個生成することにより決定される。これらの初期値を与えるための乱数生成をどのように生成するかは、パス生成パラメータファイルで指定されるパスパラメータ初期値生成条件によって規定される。
【００５３】
Ｎ個のパスパラメータの初期値を決定した後、それぞれの到来波方向に対するアンテナ指向性利得をデータ記憶部３内のアンテナ指向性利得パターンファイルから読み出して初期振幅値に乗算し、アンテナブランチ数Ｍ個に対するＮ×Ｍ個の伝搬パスパラメータを設定して、判断ステップ４へと進む。
【００５４】
図７は、パス数：１０、遅延プロファイルの形状：矩形型、遅延スプレッド：１［μｓｅｃ］、到来方向特性：［０，  ２π］の一様分布、位相特性：［０，  ２π］の一様分布及び移動局の速度：０［ｍ／ｓ］としたパスパラメータ初期値生成条件により生成したパスパラメータ初期値を示した３種のグラフである。図７（ａ）は初期時間遅延と初期振幅を表す遅延プロファイル、図７（ｂ）は初期位相、図７（ｃ）は初期到来方向を示している。例えば、上記のような条件により、生成波のパスパラメータ初期値生成条件を規定する。
【００５５】
Ｓ４：パラメータ入力・制御部２は、データ記憶部３から得たパス生成パラメータファイルの時変関数タイプ指定値に従って、時変関数をデータ記憶部３に格納した関数を用いるのか、それとも時変空間生成部６におけるモデル計算によって新たに生成するのかを選択する。通常、この選択はパス生成パラメータファイルの時変関数タイプ指定値に従ってなされるため、本判断ステップ４において、それ以降のステップ５又はステップ１１のいずれに進むのかを選択することになる（ステップ５又はステップ１１に進む）。
【００５６】
Ｓ５：パラメータ入力・制御部２は、データ記憶部３から時変関数ファイルを読み込み、それをパス生成部５と共有しあるいはパス生成部５へと供給する。該時変関数ファイルに基づき、時変関数のタイプを決定すると共に該関数の定数パラメータ生成用条件ファイルをデータ記憶部３内のパス生成パラメータファイルから読み込む。また、データサンプリング間隔及び到来波角度差に対するパラメータ間相関特性ファイルも読み込んで、ステップ６へと進む。
【００５７】
Ｓ６：パス生成部５は、パス生成パラメータファイル内の遮蔽時間間隔ＴＡ、遮蔽振幅ＬＢ、遮蔽発生時間間隔ＴＣに対する確率密度関数に従う乱数の生成・出力を乱数生成部４に対して要求してその乱数を得る。このとき、これらの乱数に相関を与える必要がある場合には、到来波角度差に対するパラメータ間相関特性ファイルで規定される初期到来方向の差に応じた相関を与えて乱数発生を要求してその乱数を得る。このようにして得られた乱数により選択された時変関数の定数パラメータを決定して、ステップ７へと進む。
【００５８】
ここで人などによる遮蔽を表した図８のグラフを参照して、時変関数ファイルとして与える時変振幅関数の例を説明する。図８（ａ）に示す関数ｆ（ｔ）においてＴＡ［ｓｅｃ］は遮蔽時間間隔、ＬＢ［ｄＢ］は遮蔽振幅、ＴＣ［ｓｅｃ］は遮蔽発生時間間隔であり、これらの値により図８（ｂ）に示すような関数形に変形できるので、関数形を決める定数パラメータとなる。そこで、これらの定数パラメータを決める条件として、例えば、遮蔽時間間隔ＴＡ［ｓｅｃ］が平均値＜ＴＡ＞［ｓｅｃ］及び標準偏差σＡ［ｓｅｃ］の正規分布に従い、遮蔽振幅ＬＢ［ｄＢ］が平均値＜ＬＢ＞［ｄＢ］及び標準偏差σＢ［ｄＢ］の正規分布に従い、かつ遮蔽発生時間間隔ＴＣ［ｓｅｃ］が平均値＜ＴＣ＞［ｓｅｃ］及び標準偏差σＣ［ｓｅｃ］の正規分布に従うと規定してあれば、それぞれの分布に従うＮ個の正規乱数を別々に乱数生成部４に要求して得ることができる。
【００５９】
図９は、トラックやバス等の大きな物による遮蔽を表した、時変関数ファイルとして与える時変振幅関数の別の例を示す図である。図９（ａ）に示す関数ｇ（ｔ）においてＴＡ［ｓｅｃ］は遮蔽時間間隔、ＬＢ［ｄＢ］は遮蔽振幅、ＴＣ［ｓｅｃ］は遮蔽発生時間間隔、Ｃｄ［ｄＢ／ｓｅｃ］は落込み係数、Ｃｒ［ｄＢ／ｓｅｃ］は立ち上がり係数、σＳ［ｄＢ］は遮蔽時標準偏差である。これらの値に依存して図９（ｂ）に示すような関数形に変形できるので、関数形を決める定数パラメータとなる。そこで、これらの定数パラメータを決める条件として、例えば、遮蔽時間間隔ＴＡ［ｓｅｃ］が平均値＜ＴＡ＞［ｓｅｃ］及び標準偏差σＡ［ｓｅｃ］の正規分布に従い、遮蔽振幅ＬＢ［ｄＢ］が平均値＜ＬＢ＞［ｄＢ］及び標準偏差σＢ［ｄＢ］の正規分布に従い、遮蔽発生時間間隔ＴＣ［ｓｅｃ］が平均値＜ＴＣ＞［ｓｅｃ］及び標準偏差σＣ［ｓｅｃ］の正規分布に従い、落込み係数Ｃｄ［ｄＢ／ｓｅｃ］が平均値＜Ｃｄ＞［ｄＢ／ｓｅｃ］及び標準偏差σＣｄ［ｄＢ／ｓｅｃ］  の正規分布に従い、立ち上がり係数Ｃｒ［ｄＢ／ｓｅｃ］が平均値＜Ｃｒ＞［ｄＢ／ｓｅｃ］及び標準偏差σＣｒ［ｄＢ／ｓｅｃ］  の正規分布に従い、かつ遮蔽時変動が平均値０及び標準偏差σＳ［ｄＢ］の正規分布に従うと規定してあれば、それぞれの分布に従うＮ個の正規乱数を別々に乱数生成部４に要求して得ることができる。図８（人による遮蔽）及び図９（トラック、バス等の大きな物による遮蔽）で示した時変関数として、予め想定した環境条件に応じたものを用意しておく。従って、本発明ではそれらを時変関数ファイルとして必要なだけデータ記憶部３に格納することにより任意に実現できる。
【００６０】
Ｓ７：パス生成部５が、ステップ６で決定した時変関数の定数パラメータに基づき、各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性を生成し、各波の振幅及び位相に対して時系列的に該時変遮蔽振幅特性を減算及び該時変遮蔽位相特性を加算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する。それと同時に、パス生成部５は、移動局の速度及び初期到来波方向（移動局の進行方向を基準に規定される）の情報から、各波の時間的位相回転Δφｉ（ｔ）を式（１）で計算し、これをそれぞれのパスについての初期位相に加算することにより、Ｎ個の伝搬パスの時変位相特性にドップラー位相回転の効果を与えることができる。そしてステップ８へと進む。
【００６１】
Ｓ８：パス生成部５が、次回の乱数発生によって継続加算される時変データが現在の結果に影響を及ばさないデータであることを確認し、出力データとしてバッファメモリ７への転送を開始する（ステップ９に進む）。
【００６２】
Ｓ９：パス生成部５が、生成した各パスの振幅及び位相に対する遮蔽変動特性の生成データを順次繰り返し、時系列的に継続接続・加算して出力データを生成し続ける。すなわち、パス生成部５は、パラメータ入力・制御部２からのパス生成停止要求が有るか否かをチェックし、停止要求が無い場合にはステップ６に戻って各パスの振幅及び位相に対する遮蔽変動特性の生成を繰り返す。この際、各パスの振幅及び位相に加減算する時変遮蔽振幅特性及び時変遮蔽位相特性は遮蔽発生時間間隔毎に継続接続して生成する。
【００６３】
従って、遮蔽時間間隔ＴＡ［ｓｅｃ］が遮蔽発生時間間隔ＴＣ［ｓｅｃ］よりも大きい値になる場合には、図１０に示すように遮蔽時間が重複する（図１０（ｂ）の上のグラフ参照）。この場合にはこれら２つの遮蔽特性を単純加算したものを時変遮蔽振幅特性及び時変遮蔽位相特性とする（図１０（ｂ）の下のグラフ参照）。
【００６４】
一方、パラメータ入力・制御部２からパス生成停止要求が有る場合には、停止要求が生じるまでの生成データを出力データのバッファメモリ７に転送して終了し、ステップ１０へと進む。
【００６５】
Ｓ１０：出力データのバッファメモリ部７が、パス生成部５から転送されてきた生成データを格納すると共に、予め設定しておいた規定データ量を超えた時点から該格納データを出力インターフェース部８に対して出力する。このとき、パラメータ入力・制御部２からの実行停止要求が有れば、バッファメモリ部７のメモリ使用容量が前記規定データ量を下回っても全てのデータを吐き出すまで該格納データを出力インターフェース部８に対して出力して終了する（パス生成終了）。
【００６６】
Ｓ１１：パラメータ入力・制御部２が、先ず、パス生成パラメータファイルの遮蔽変動計算モデルファイル、伝搬パスの周波数Ｆ［ＧＨｚ］及び送受信点間距離（最短パスの伝搬距離）Ｌ０［ｍ］に基づき、時変関数生成部４に対して計算モデルを指定して、ステップ１２へと進む。
【００６７】
Ｓ１２：パス生成部５が、遮蔽物寸法（Ｗｈ［ｍ］，Ｄｈ［ｍ］，Ｈｈ［ｍ］）の確率密度関数、伝搬パスへの遮蔽物体の進入角度θｓ［ｄｅｇｒｅｅ］  の確率密度関数、遮蔽物の進入位置Ｌｄ［ｍ］  の確率密度関数、遮蔽物の移動速度Ｖｓ［ｍ／ｓ］  の確率密度関数にそれぞれ従うような、遮蔽物寸法（Ｗｈ［ｍ］，Ｄｈ［ｍ］，Ｈｈ［ｍ］）、伝搬パスへの遮蔽物体の進入角度θｓ［ｄｅｇｒｅｅ］、遮蔽物の進入位置Ｌｄ［ｍ］、遮蔽物の移動速度Ｖｓ［ｍ／ｓ］に対する各乱数を、乱数生成部４に対して要求して生成させて各乱数を得る。こうして得た乱数パラメータを時変関数生成部６へと転送して、計算モデルの各種パラメータを決定し、時変関数生成部６にそのモデルでのＮ個の時変関数の生成を要求する。このとき、これらの乱数生成に相関を与える必要がある場合には、到来波角度差に対するパラメータ間相関特性ファイルで規定される初期到来方向の差に応じた相関を与えて乱数発生を要求して乱数を得る。その後ステップ１３へと進む。
【００６８】
図１１に、計算モデルによって時変振幅関数及び時変位相関数を生成する場合の計算手順のモデルのフローチャートを示す。
【００６９】
図１２に計算モデルの概要を図示する。図１２に示す実施例で示す計算モデルにおいては、送受信点間距離（最短パスの伝搬距離）Ｌ０［ｍ］、伝搬パスの周波数Ｆ［ＧＨｚ］、移動局からの遮蔽物の進入位置Ｌｄ［ｍ］、遮蔽物の進入角度θ［ｄｅｇｒｅｅ］を与え、また遮蔽物の大きさを短半径ａ［ｍ］、長半径ｂ［ｍ］の楕円柱として与えている。遮蔽物が完全に見通し内のパス上に存在するとき、受信電界Ｅは以下のように表される。
【００７０】
【数３】
2004254250.tif 000005
ここで、Ｅ１，Ｅ２はそれぞれナイフエッジ１，２にて回折し到来する電波の電界強度である。
さらに、Ｅ１，Ｅ２は式（３）の形で与えられる。
【００７１】
【数４】
2004254250.tif 000006
ここで、Ｅｆは遮蔽物がない場合の受信点における自由空間電界であり、
【００７２】
【数５】
2004254250.tif 000007
ただしＰは送信アンテナの放射電力である。Ｃ、Ｓはフレネル積分であり、次のように表される。
【００７３】
【数６】
2004254250.tif 000008
ここでωはクリアランス関数と呼ばれ、
【００７４】
【数７】
2004254250.tif 000009
と表される。Ｌ０は送受信間距離［ｍ］、λは波長［ｍ］、（ｘ１，  ｙ１），  （ｘ２，  ｙ２）  は図１２で示すナイフエッジ１，２の座標を意味する。
【００７５】
図１１を参照しながら、本計算モデルによる時変関数の計算手順を説明する。先ずステップＳ４０１において、送受信点間距離（最短パスの伝搬距離）Ｌ０［ｍ］及び移動局からの遮蔽物の進入位置Ｌｄ［ｍ］、遮蔽物の進入角度θ［ｄｅｇｒｅｅ］、及び遮蔽物寸法（本モデルでは、短半径ａ［ｍ］、長半径ｂ［ｍ］の楕円柱を仮定）のパラメータを用いて遮蔽物の中心が見通し線上に位置する場合のナイフエッジ１，２の位置を決め、ナイフエッジ１及びナイフエッジ２での回折パスの電界Ｅ１，  Ｅ２  を計算し、式（２）にてベクトル合成した電界Ｅ  を求める。
【００７６】
ステップＳ４０２において、遮蔽物の位置を見通し線からデータサンプリング間隔×遮蔽物の速度［ｍ］で決まるステップ距離だけ戻し（図１２において下方向）、ステップＳ４０１と同様にナイフエッジ１，２の位置を決定し、回折パス１，２の電界を計算し、ベクトル合成した電界Ｅを求める。受信電界Ｅが遮蔽物が無い状態における電界（直接波の電界）Ｅｆを上回るまで遮蔽物の位置をステップ距離だけ繰り返し戻してこの手順を続ける。
【００７７】
ステップＳ４０３において、遮蔽物の位置を見通し線からデータサンプリング間隔×遮蔽物の速度［ｍ］  で決まるステップ距離だけ進めた（図１２において上方向）位置とし、ステップＳ４０１と同様にナイフエッジ１，２の位置を決定し、回折パス１，２の電界を計算し、ベクトル合成した電界Ｅを求める。この手順を、電界Ｅが遮蔽物が無い状態におけるレベルを上回るまで遮蔽物の位置をステップ距離だけ繰り返し進めて続ける。
【００７８】
ステップＳ４０４において、ステップＳ４０１、ステップＳ４０２及びＳ４０３で求めた電界Ｅの振幅及び位相を遮蔽物の移動時間に対応させ、時変関数としてメモリ領域に保存する。
【００７９】
図１３はＬ０＝１３．５［ｍ］、Ｆ＝７［ＧＨｚ］、Ｌｄ＝３［ｍ］、θ＝４５［ｄｅｇ］及びａ＝０．４９［ｍ］、ｂ＝０．２０［ｍ］の計算モデルパラメータを仮定した場合の受信電力変動の計算例を図示（実線）しており、受信電力Ｐと電界Ｅとの関係は次式で与えられる。
【００８０】
【数８】
2004254250.tif 000010
図１３には、同条件で行った測定結果も併せて示してある（点線）。尚、測定では遮蔽物を人体としており、横軸は見通し線上に人体が存在している場合の位置を基準としている。図１３より、受信電力変動特性、即ち時変振幅関数は本計算モデルで良く近似できていることが分かる。
【００８１】
Ｓ１３：時変関数生成部６で生成されるＮ個の時変関数（時変振幅関数及び時変位相関数）が各Ｎ波のパスの時変遮蔽振幅特性及び時変遮蔽位相特性であることから、パス生成部５がこれらを各波の振幅及び位相に対して時系列的に該時変遮蔽振幅特性を減算しかつ該時変遮蔽位相特性を加算することにより、各パスの振幅及び位相に対する遮蔽変動特性を生成する。それと同時に、パス生成部５は、移動局の速度及び初期到来波方向（移動局の進行方向を基準に規定される）の情報から、各波の時間的位相回転Δφｉ（ｔ）を式（１）で計算し、これをそれぞれのパスに対する初期位相に加算することにより、Ｎ個の伝搬パスの時変位相特性にドップラー位相回転の効果を与えることができる。その後ステップ１４に進む。
【００８２】
Ｓ１４：ステップ１２で決定した各乱数に基づき、ステップ１３で生成した各パスの振幅及び位相に対する遮蔽変動特性の生成データを、次のステップ１５において順次繰り返し図１０に示すように時系列的に継続接続・加算して出力データを生成し続けるが、本ステップ１４において、パス生成部５が、次の乱数発生によって継続加算される時変データが現在の結果に影響を及ばさないデータであることを確認して、出力データとしてバッファメモリ７への転送を開始する（ステップ１５に進む）。
【００８３】
Ｓ１５：パス生成部５は、パラメータ入力・制御部２からのパス生成停止要求が有るか否かをチェックし、停止要求が無い場合にはステップ１２に戻り各パスの振幅及び位相に対する遮蔽変動特性の生成を繰り返す。この際、各パスの振幅及び位相に加減算する時変遮蔽振幅特性及び時変遮蔽位相特性は遮蔽発生時間間隔毎に継続接続して生成する。従って、ステップ９で記述したように、図１０に示すような遮蔽時間が重複する時変特性を生成する場合がある。この場合にはこれら２つの遮蔽特性を単純加算したものを時変遮蔽振幅特性及び時変遮蔽位相特性とする。
【００８４】
一方、パラメータ入力・制御部２からのパス生成停止要求が有る場合には、停止要求が生じるまでの生成データを出力データのバッファメモリ７に転送して終了し、ステップ１０へと進む。
【００８５】
図１４は、各パスにレベル変動を与え終わった状態の例を示すグラフである。遮蔽発生間隔は平均値０．５［ｓｅｃ］の指数分布とし、移動局の速度は０［ｍ／ｓ］、各遮蔽に対し全く同一の進入位置、進入角度、遮蔽物の大きさ、及び遮蔽物の移動速度を与えた場合の各パスの時変振幅変動の一例である。本例では各遮蔽の遮蔽振幅ＬＢ［ｄＢ］、遮蔽時間間隔ＴＡは１［ｓｅｃ］となっており、また、図１０に示すような遮蔽の重複がない場合の例である。これらの遮蔽変動を与える条件は、ここに表記した条件に留まらず、模擬しようとする環境条件に応じた計算条件に代えられるものである。
【００８６】
図１５に、本実施形態に係る多重波シミュレータ１３（図４）のパス出力部１０の構成例を示す。パス出力部１０は、実部／虚部分離変換部１５、Ｄ／Ａ変換器１６から構成される。実部／虚部分離変換部１５は、時変マルチパス生成装置から入力されるＭ×Ｎ個の伝搬パス変動の時変複素振幅を実部及び虚部の信号に分離し、Ｄ／Ａ変換器１６に出力する。Ｍ×２Ｎ個のＤ／Ａ変換器１６は、デジタル／アナログ変換を行い、時変パス信号をアナログ信号として出力する。
【００８７】
図１６に、本実施形態に係る多重波シミュレータ１４の信号合成部１１の構成例を示す。信号合成部１１は、直交信号生成器１７、トランスバーサル型回路１８、遅延素子１９、乗算器２０、加算器２１から構成される。各乗算器２０には乗算ウェイトとしてパス出力部１０のアナログ出力を与える。直交信号生成器１７はデジタル信号入力端子１２から入力された信号（Ｉ成分）から直交信号（Ｑ成分）を生成する。Ｉ成分及びＱ成分はそれぞれ第１信号から第Ｍ信号までＭ分割される。ここで、Ｎはパスの数であり、Ｍはアンテナの数である。Ｉ成分の第１信号は第Ｉ−１番、第２信号は第Ｉ−２番、…第Ｍ信号は第Ｉ−Ｍ番のトランスバーサル型回路１８に入力される。同様にＱ成分の第１信号は第Ｑ−１番、第２信号は第Ｑ−２番、…第Ｍ信号は第Ｑ−Ｍ番のトランスバーサル型回路１８に入力される。
【００８８】
このような構成となっているので、パス出力部１０の出力は、アンテナの個数Ｍに対応するＭ個のグループ化されたＮ個のパス信号であり、それぞれ実部信号と虚部信号がある。即ち、ｗＲ（１１），  ｗＲ（１２），…ｗＲ（１Ｎ），ｗＲ（２１），  ｗＲ（２２），…ｗＲ（２Ｎ），…，ｗＲ（Ｍ１），  ｗＲ（Ｍ２），…ｗＲ（ＭＮ），及びｗＩ（１１），  ｗＩ（１２），…ｗＩ（１Ｎ），ｗＩ（２１），  ｗＩ（２２），…ｗＩ（２Ｎ），…，ｗＩ（Ｍ１），  ｗＩ（Ｍ２），…ｗＩ（ＭＮ）の信号が出力される。
【００８９】
また、各乗算器２０に与えるウェイトは第Ｉ−１番のトランスバーサル型回路１８の場合は前段から後段に向ってそれぞれｗＲ（１１），  ｗＲ（１２），…ｗＲ（１Ｎ）のパス変動信号を与え、第Ｉ−Ｍ番のトランスバーサル型回路１８の場合は前段から後段に向ってそれぞれｗＲ（Ｍ１），  ｗＲ（Ｍ２），…ｗＲ（ＭＮ）のパス変動信号を与える。同様にして、第Ｑ−１番のトランスバーサル型回路１８の場合は前段から後段に向ってそれぞれｗＩ（１１），  ｗＩ（１２），…ｗＩ（１Ｎ）のパス変動信号を与え、第Ｑ−Ｍ番のトランスバーサル型回路１８の場合は前段から後段に向ってそれぞれｗＩ（Ｍ１），  ｗＩ（Ｍ２），…ｗＩ（ＭＮ）のパス変動信号を与える。各トランスバーサル型回路１８の出力は、アンテナに対応するグループ毎にＩ成分及びＱ成分を加算し、Ｍ個の出力信号を得る。これによって、アンテナ毎に指向性パターンの影響を加味した出力を入力信号に乗算することができる。
【００９０】
図１７において、本発明の実施形態に係るフェージングシミュレータが出力した結果の例を細い線で示す。Ｋ＝−∞，  ３，  ６，  １０，  １５［ｄＢ］とした場合のそれぞれの仲上ライス分布を太い線で示す。
【００９１】
また、図１８において、屋内環境にて受信レベル特性について行った測定結果の一例を細い線で示す。Ｋ＝−∞，  ６，  １０［ｄＢ］とした場合のそれぞれの仲上ライス分布を太い線で示す。図面に示されるように、本実施形態におけるフェージングシミュレータの出力結果や測定結果の一例は共にＫファクタを変更した仲上ライス分布と良く一致している。
【００９２】
【発明の効果】
以上説明したように、本発明の実施形態に従えば、再現しようとするシミュレーション条件に対応する伝搬パスパラメータ（レベル特性、位相特性、時間遅延特性、到来方向特性）の時間変動特性（時変マルチパス変動）を、従来の単純なレイリーフェージング変動や仲上−ライスフェージング変動などのように全体として近似的に与えるのではなく、各伝搬パスをパス単位で生成し、かつ空間的なパスの相関特性を考慮した時変パス特性として生成することができる。また、本発明の実施形態で生成した複数のパスをベクトル合成することにより、移動局が静止した状態でもレイリーフェージング変動や仲上−ライスフェージング変動などが再現できることが確認でき、本来の伝搬メカニズムで生起するパス変動を模擬することのできる多重波フェージングシミュレータを作ることができ、実際の伝搬環境の特性を再現する上で極めて有効かつ重要なものである。
【図面の簡単な説明】
【図１】基地局から送信された電波が反射や回折を経て移動局に到達する様子の概念図及び各波の遅延プロファイルを示すグラフである。
【図２】移動局が移動している場合の遅延プロファイルの測定例を示すグラフである。
【図３】本発明の実施形態に係る時変マルチパス生成装置の構成例を示すブロック図である。
【図４】本発明の実施形態に係る多重波シミュレータの構成例を示すブロック図である。
【図５】本発明の実施形態に係る受信局の移動による位相の変化を示す図である。
【図６】本発明の実施形態に係る時変マルチパス生成装置の動作を示すフローチャートである。
【図７】本発明の実施形態に係るパスパラメータ初期値の生成結果の例を示すチャートである。
【図８】本発明の実施形態に係る時変関数の例を示すグラフである。
【図９】本発明の実施形態に係る他の時変関数の例を示すグラフである。
【図１０】本発明の実施形態に係る遮蔽時間の重複が生じるケースでのパス変動出力を説明するグラフである。
【図１１】本発明の実施形態に係る時変関数生成部の動作を示すフローチャートである。
【図１２】遮蔽によるレベル変動計算の例を示す模式図である。
【図１３】遮蔽によるレベル変動計算結果と測定値を比較するためのグラフである。
【図１４】各パスに遮蔽によるレベル変動を与え終わった状態の例を示すグラフである。
【図１５】本発明の実施形態に係るパス出力部の構成例を示すブロック図である。
【図１６】本発明の実施形態に係る信号合成部の構成例を示すブロック図である。
【図１７】本発明の実施形態に係るフェージングシミュレータが出力した結果の例を仲上ライス分布とともに示すグラフである。
【図１８】合成した受信レベルの測定結果の一例を仲上ライス分布とともに示すグラフである。
【符号の説明】
１    入力デバイス
２    パラメータ入力・制御部
３    データ記憶部
４    乱数生成部
５    パス生成部
６    時変関数生成部
７    出力データバッファメモリ部
８    出力インターフェース部
９    時変マルチパス生成装置
１０    パス出力部
１１    信号合成部
１２    デジタル信号入力端子
１３    デジタル信号出力端子
１４    多重波シミュレータ
１５    実部／虚部分離変換部
１６    Ｄ／Ａ変換器
１７    直交信号生成器
１８    トランスバーサルアレイ回路
１９    遅延素子
２０    乗算器
２１    加算器
【図１】
2004254250.tif 000011
【図２】
2004254250.tif 000012
【図３】
2004254250.tif 000013
【図４】
2004254250.tif 000014
【図５】
2004254250.tif 000015
【図６】
2004254250.tif 000016
【図７】
2004254250.tif 000017
【図８】
2004254250.tif 000018
【図９】
2004254250.tif 000019
【図１０】
2004254250.tif 000020
【図１１】
2004254250.tif 000021
【図１２】
2004254250.tif 000022
【図１３】
2004254250.tif 000023
【図１４】
2004254250.tif 000024
【図１５】
2004254250.tif 000025
【図１６】
2004254250.tif 000026
【図１７】
2004254250.tif 000027
【図１８】
2004254250.tif 000028
【手続補正書】
【提出日】平成15年2月25日(2003.2.25)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００８５
【補正方法】変更
【補正の内容】
【００８５】
図１４は、各パスにレベル変動を与え終わった状態の例を示すグラフである。遮蔽発生間隔は平均値０．５［ｓｅｃ］の指数分布とし、移動局の速度は０［ｍ／ｓ］、各遮蔽に対し全く同一の進入位置、進入角度、遮蔽物の大きさ、及び遮蔽物の移動速度を与えた場合の各パスの時変振幅変動の一例である。本例では各遮蔽の遮蔽振幅ＬＢは１０［ｄＢ］、遮蔽時間間隔ＴＡは１［ｓｅｃ］となっており、また、図１０に示すような遮蔽の重複がない場合の例である。これらの遮蔽変動を与える条件は、ここに表記した条件に留まらず、模擬しようとする環境条件に応じた計算条件に代えられるものである。