JPH07297859A - デジタル変調解析装置 - Google Patents

Abstract

(57)【要約】 【目的】 本発明は、変調信号の解析を短時間で解析で
きるようにし、かつ小型安価な測定装置を構成すること
を目的とする。 【構成】 ＤＵＴの信号を中間周波数に変換する周波数
変換部１９０を設け、さらに低い中間周波数に変換しＡ
Ｄ変換してメモリに保存するサンプリング部２５０を設
け、直交検波してＩ成分とＱ成分に変換するヒルベルト
変換部３１０を設け、キャリア中心周波数Ωf0と、スロ
ットの同期アドレスdly0を求める初期同期部７４０を設
け、変調誤差εが最小となる最適パラメータ（ゲイン定
数、初期位相、キャリア周波数、バースト・ドループ・
ファクタ）を求め、全シンボルの復調データ（振幅値Ａ
［ｎ］、位相角θ［ｎ］）を復調し、全シンボルの理想
データ（振幅値Ａｒ［ｎ］、位相角値θｒ［ｎ］）を特
定する全シンボル復調／パラメータ推定部４００を設
け、変調精度εデータや各シンボル毎の変調誤差データ
を求める変調精度演算部９００を設ける手段。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】この発明は、デジタル通信方式の
１つであるデジタル・マルチ・チャンネル・アクセス
（ＤＭＣＡ＝Digital Multi Channel Acces）等の、デ
ジタル方式変調の変調度の測定に関する。

【０００２】

【従来の技術】デジタル通信方式の中の民間業務用移動
体通信システムは、ＤＭＣＡの変調方式であるＭ１６Ｑ
ＡＭ（Multi Subcarrier 16 Quadrature Amplitude Mod
ulation,M=4）が採用されている。このＭ１６ＱＡＭ
は、１６ＱＡＭ信号を４つの副搬送波（サブキャリア）
で同時に送信するマルチ・キャリア方式である。この通
信方式によって、送信データの伝送速度は、６４Ｋｂｐ
ｓ（ビット／秒）のデータを４Ｋシンボル／秒の伝送速
度にしている。これについての文献は、財団法人電波シ
ステム開発センター発行のデジタル方式ＭＣＡシステム
標準規格（ＲＣＲーＳＴＤー３２）があり、この文献中
には、変調方式や、符号化や、同期／パイロット挿入
や、下り・上りスロットの構成図や、測定法や、変調精
度の定義等の詳細が説明されている。

【０００３】情報データは、図１４に示すように、１６
ビットを単位として符号化して送信する。この１６ビッ
トデータは、４ビット４組に分割して、４つのサブキャ
リアｃｈ１〜ｃｈ４に乗せて伝送する。各サブキャリア
は、ロールオフ＝０．２のルート・ナイキスト・フィル
タ（帯域圧縮フィルタ）で帯域を圧縮して、各サブキャ
リアの間隔は４．５ＫＨｚにし、この４チャンネルで１
８ＫＨｚの帯域を使用している。この４ｃｈサブキャリ
アの送信系統図を図１５に参考に示す。

【０００４】４ビット単位データは、図１２に示すよう
に、２ビットのＩ成分と２ビットのＱ成分に分けて、こ
れを直交符号化して１６値符号の組み合わせで４ビット
情報を表現している。以後このシンボルのことを理想情
報シンボル（Reference Information Symbol＝ＲＩＳ）
と称す。送信側では、この情報シンボル以外に図１７に
示すように、各サブチャンネル毎に所定のシンボル番号
位置に同期シンボルとパイロットシンボルが挿入されて
いる。この理想同期／パイロットシンボル（Reference
Pilot Symbol＝ＲＰＳ）は、図１２に示すように振幅値
が一定の√１８の値として√１８の円周上に配置され、
かつ図１６に示すようにシンボル番号毎に規定された位
相角が与えられている。これにより受信側のＩ成分とＱ
成分から同期／パイロットデータと情報データの復調が
可能になっている。通常、受信側で復調した場合には、
情報シンボル（Information Symbol＝ＩＳ）や同期／パ
イロットシンボル（Pilot Symbol＝ＰＳ）の位置は、各
理想シンボル点の位置から外れていて誤差を持ってい
る。

【０００５】スロット構成は、図１７に示すように１５
ｍｓｅｃ時間を基本としていて、下り基本スロットの構
成では、６０シンボルで構成していて、この６０シンボ
ルの中には、３同期シンボルと７パイロットシンボルと
５０情報シンボルとで構成している。また、上り基本ス
ロットの構成では、５３シンボルとＡＧＣプリアンブル
とＲａｍｐで構成していて、この５３シンボルの中に
は、３同期シンボルと７パイロットシンボルと４３情報
シンボルとで構成している。また、上りサブスロットの
構成では、基本スロット２０シンボルとＡＧＣプリアン
ブルとＲａｍｐで構成していて、この２０シンボルの中
には、３同期シンボルと３パイロットシンボルと１４情
報シンボルとで構成している。これらの同期シンボルと
パイロットシンボルは、図１６に示すように、スロット
の先頭シンボルをシンボル番号１から始まる連続番号と
して、個々の同期シンボルとパイロットシンボル番号毎
に位相角度が決められている。またフレーム構成は、９
０ｍｓｅｃを基本にしていて、６スロットを基本単位構
成している。これは、６チャンネル多重の上り用の時分
割多重接続方式、すなわちＴＤＭＡ（Time Division Mu
ltiple Acces）で通信する場合に使われる。

【０００６】

【数１】

【０００７】この数１は、変調誤差εの計算式３１であ
る。本デジタル通信方式における変調精度は、変調誤差
計算式３１の式に示すように定義されている。この計算
式で、ｉ：観測ベクトル番号、Ｎ：総観測ベクトル、ｒ
0：信号空間ダイヤグラムの最大半径（√１８）、Ｖm
i：ｉ番目の観測ベクトル、Ｖri：ｉ番目の観測ベクト
ルの理想ベクトル、Ｖ0：原点オフセット・ベクトル、
α：ゲインパラメータ（スカラ量）、φ：位相パラメー
タ（スカラ量）、としていて、Ｖ0、α、φの各パラメ
ータは変調精度εが最小となるように選択するものとす
る。図１３の変調誤差εiは、１個のｉ番目の観測ベク
トルについての変調誤差を示した図である。

【０００８】従来の実施例について、図１０と図１１を
参照して、ＤＭＣＡの変調精度の測定を説明する。本装
置の構成は、図１０に示すように、デジタル変換部１９
０と、信号発生器１８３と、ワークステーション１８４
と、解析ソフト１７０とで構成している。このデジタル
変換部１９０は、周波数変換器１９２と、ＬＰＦ１９３
と、ＡＤ変換器１９４と、バッファメモリ１９５と、制
御部１９６と、ＧＰＩＢインターフェイス１９７とで構
成している。また、この解析ソフト１７０は、図１１に
示すように、チャンネル分離部１７２と、ＩＱ変換部１
７３ａ〜１７３ｄと、ルート・ナイキスト・フィルタ部
１７４ａ〜１７４ｄと、クロック推定部１７５ａ〜１７
５ｄと、キャリア・オフセット推定部１７６ａ〜１７６
ｄと、変調精度誤差演算部１７７とで構成している。

【０００９】被試験器からの被測定信号１９１は、周波
数変換器１９２の一方に入力し、高確度の信号発生器１
８３の信号を他方の入力に与えて目的の中間周波数、例
えば１００ＫＨｚに変換した後、次のローパス・フィル
タＬＰＦ１９３で低域周波数信号を通過させた後、ＡＤ
変換器１９４に供給する。ＡＤ変換器１９４では、この
信号を受けて、１２ビットのデジタル信号に変換した
後、バッファメモリ１９５に記憶する。このメモリ容量
は１Ｍワードあり、９０ｍｓｅｃであるフレームを複数
記録することができる。このＡＤ変換器１９４によるサ
ンプリングと記録は、ワークステーションのＧＰＩＢイ
ンターフェイス１９７のコマンドを制御部１９６が受け
て実行制御する。また、サンプリングされたデータは、
バッファメモリ１９５からＧＰＩＢインターフェイス１
９７を経由してワークステーション１８４の記憶部に転
送される。この転送されたＤＭＣＡデータに対して以後
解析処理を実行して変調精度を求めている。

【００１０】ワークステーション側では、上記の転送デ
ータを受けて、解析ソフト１７０により解析処理を実行
する。まず、チャンネル分離部１７２では、ＤＭＣＡデ
ータを高速フーリエ変換（ＦＦＴ）して周波数軸上に変
換し、Ｍ１６ＱＡＭ信号を４チャンネルのサブチャンネ
ル信号（１６ＱＡＭ）に分離する。次に、ＩＱ変換部１
７３ａ〜１７３ｄでは各チャンネル毎に直交変換してＩ
Ｑ信号に変換する。ルート・ナイキスト・フィルタ部１
７４ａ〜１７４ｄでは、この信号を受けてロールオフ率
０．２のルート・ナイキスト・フィルタ処理をした後、
クロック推定部１７５ａ〜１７５ｄに供給する。

【００１１】クロック推定部１７５ａ〜１７５ｄでは、
この信号を受けて各サブチャンネルのＩＱ信号に対して
非直線性処理を施すと、変調信号のクロック周波数のス
ペクトルが得られ、このスペクトルからクロックの周波
数及び位相の推定を行う。キャリア・オフセット推定部
１７６ａ〜１７６ｄでは、上記でシンボル点が判明して
いるので、予め定められている同期シンボルやパイロッ
トシンボルの位相とシンボル点での位相差から、キャリ
アオフセットの推定を１バースト毎に行う。

【００１２】変調精度誤差演算部１７７では、上記のシ
ンボル点からベクトル計算して、図１３に示す理想ベク
トル点の位置からの誤差ベクトルから変調精度（Vector
Error）、クロック位相誤差（Clock Alignment）、周
波数誤差（Carrier FrequencyOffset）が求まる。この
各種誤差データから、ベクトル誤差の振幅・位相波形、
ベクトル誤差のＦＦＴ、アイパターン、コンスタレーシ
ョン（Constellation）等を画面に表示する。

【００１３】

【発明が解決しようとする課題】上記説明の構成による
実現手段において、ＦＦＴ演算処理等、汎用のワークス
テーションのＣＰＵによる演算の為に処理速度が遅い点
や、大量のＡＤ変換データをそのままＧＰＩＢで転送す
る為に時間がかかる欠点の為に、１回の処理時間が、分
単位の頻度でしか測定できないという不便があった。ま
た、外部に信号発生器１８３と、ワークステーション１
８４とが必要であり、装置が大型になり高価なシステム
構成となっていて好ましくない。

【００１４】そこで、本発明が解決しようとする課題
は、変調信号の解析を短時間で解析できるようにし、か
つ小型安価な測定装置を構成することを目的とする。

【００１５】

【課題を解決する為の手段】第１図は、本発明による第
１の解決手段を示している。上記課題を解決するため
に、本発明の構成では、ＤＵＴの信号を受けて、一定の
中間周波数に変換する周波数変換部１９０を設け、周波
数変換部１９０からの信号を受けて、さらに低い中間周
波数信号に変換した後ＡＤ変換してメモリに保存するサ
ンプリング部２５０を設け、直交検波してＩ成分とＱ成
分に変換するヒルベルト変換部３１０を設け、Ｉ成分と
Ｑ成分信号を受けて、キャリア中心周波数Ωf0と、スロ
ットの同期アドレスdly0を求める初期同期部７４０を設
ける。そして、初期同期部７４０の信号とＩ成分とＱ成
分信号を受けて、変調誤差εが最小となる最適パラメー
タ（ゲイン定数Ａx、初期位相θx、キャリア周波数Ω
x、バースト・ドループ・ファクタσx）を求め、かつ、
全シンボルの復調データ（シンボル番号ｎの振幅値Ａ
［ｎ］、位相角θ［ｎ］）を復調し、全シンボルの理想
データ（理想シンボルの振幅値Ａｒ［ｎ］、理想シンボ
ルの位相角値θｒ［ｎ］）を特定する全シンボル復調／
パラメータ推定部４００を設け、 全シンボル復調／パ
ラメータ推定部４００の信号を受けて、変調精度εデー
タや各シンボル毎の変調誤差データを求める変調精度演
算部９００を設ける構成手段である。これにより、デジ
タル通信システムでＤＵＴが出力するデジタル方式のＭ
ＣＡ信号の変調精度を測定する変調解析手段としてい
る。

【００１６】第１２図は、本発明による第２の解決手段
を示している。上記課題を解決するために、本発明の構
成では、スペクトラムアナライザ１１０内に、ＤＵＴの
信号を受けて、一定の中間周波数信号に変換する周波数
変換部１１１を設け、測定部２００からの信号を受け
て、変調誤差の画面表示用データに変換する表示制御部
１１５を設け、スペクトラムアナライザ１１０からの周
波数変換信号を受けて、さらに低い中間周波数信号に変
換した後ＡＤ変換してメモリに保存するサンプリング部
２５０を設け、直交検波してＩ成分とＱ成分に変換する
ヒルベルト変換部３１０を設け、Ｉ成分とＱ成分信号を
受けて、キャリア中心周波数Ωf0と、スロットの同期ア
ドレスdly0を求める初期同期部７４０を設ける。そし
て、初期同期部７４０の信号とＩ成分とＱ成分信号を受
けて、変調誤差εが最小となる最適パラメータ（ゲイン
定数Ａx、初期位相θx、キャリア周波数Ωx、バースト
・ドループ・ファクタσx）を求め、かつ、全シンボル
の復調データ（シンボル番号ｎの振幅値Ａ［ｎ］、位相
角θ［ｎ］）を復調し、全シンボルの理想データ（理想
シンボルの振幅値Ａｒ［ｎ］、理想シンボルの位相角値
θｒ［ｎ］）を特定する全シンボル復調／パラメータ推
定部４００を設け、全シンボル復調／パラメータ推定部
４００の信号を受けて、変調精度εデータや各シンボル
毎の変調誤差データを求める変調精度演算部９００を設
ける構成手段である。これにより、デジタル通信システ
ムでＤＵＴが出力するデジタル方式のＭＣＡ信号の変調
精度を測定する変調解析手段としている。

【００１７】

【作用】最適パラメータ推定の計算式は、最適パラメー
タの初期値Ａin、θin、Ωin、Ｂin、σinを使うことに
より、近似計算により高速に計算できる作用がある。初
期同期部７４０のキャリア中心周波数検出部７００で
は、概略キャリア中心周波数Ω値を求めている。これに
より、例えば±１０ＫＨｚ程度の中間周波数のずれがあ
った場合でも、支障なく、以後のシンボル復調や最適パ
ラメータ推定を正常に実施できる作用がある。初期同期
部７４０の周波数エラー検出部７６０は、概略キャリア
中心周波数Ω値を受けて、更に正確なキャリア中心周波
数Ωf0が得られる。

【００１８】

【実施例】

（実施例１）本発明の実施例は、周波数変換部１９０を
用いた場合の例である。これについて、図１と図２
（ａ）、（ｂ）と図３（ａ）、（ｂ）、（ｃ）と図４と
図５と図６と図７（ａ）〜（ｅ）と図８を参照して以下
に説明する。最初にシステム構成を説明する。システム
構成は、図１に示すように、周波数変換部１９０と、測
定部２００とでシステムを構成している。周波数変換部
１９０は、信号発生器１０３と、第１ミキサ１０２とで
構成している。測定部２００の内部構成は、サンプリン
グ部２５０と、ヒルベルト（Hilbert）変換部３１０
と、初期同期部７４０と、全シンボル復調／パラメータ
推定部４００と、変調精度演算部９００と、出力部２３
０と制御部２８０とで構成している。サンプリング部２
５０は、第２ミキサ２０２と、第１発振器２０３と、第
１ＬＰＦ（ローパス・フィルタ）２１０と、ＡＤ変換器
２２０と、クロック２２１と、ＡＤバッファメモリ２４
０とで構成している。

【００１９】次に、変調精度の測定全体の動作概要を説
明する。周波数変換部１９０は、被測定信号１０１の周
波数を所定の周波数ｆ２０２＝２１．４２ＭＨｚにダウ
ンコンバート（Down Convert）するものである。この為
に、信号発生器１０３は、所定の周波数ｆ２０２にする
為の周波数信号を第１ミキサ１０２に供給する。第１ミ
キサ１０２は、被測定信号１０１と、信号発生器１０３
の両信号を受けて、所定の周波数にダウンコンバートし
て、測定部２００の第２ミキサ２０２に供給する。第２
ミキサ２０２は、上記信号と第１発振器２０３の２１．
２ＭＨｚの両信号を受けて、２２０ＫＨｚの信号にダウ
ンコンバートする。この出力信号は、第１ＬＰＦ２１０
でローパス・フィルタした後、ＡＤ変換器２２０に供給
する。

【００２０】ＡＤ変換器２２０は、制御部２８０からの
サンプリング指令により、クロック２２１の１００８Ｋ
Ｈｚのサンプリング・クロックを受けて、上記信号を１
０ビットのワードデータに量子化した後、ＡＤバッファ
メモリ２４０に格納保存する。これにより、このＡＤバ
ッファメモリには、１シンボル時間が４ＫＨｚであるか
ら１００８ＫＨｚ／４ＫＨｚ＝２５２ワード／シンボル
の分解能で格納される。ここで、測定に必要なデータの
格納容量は、ＴＤＭＡ通信信号の場合でも有効な１スロ
ットが確保できるようにする為に、９スロット時間分以
上の格納容量を持たせておく。

【００２１】以上説明の動作は、回路で実現している。
以後の動作は、このメモリに格納されたスロットデータ
を順次読みだしてソフト的に解析処理している。本実施
例では、この解析処理を高速にする為にＤＳＰ（Digita
l Signal Processor）を用いている。無論、ＤＳＰの代
わりに、他の高速演算するＣＰＵ等を使用しても良い。

【００２２】次に、ヒルベルト変換部３１０は、ＡＤバ
ッファメモリ２４０より読みだした２２０ＫＨｚ中間周
波数信号のデジタルデータ２４０ａを直交検波（Quadra
nt Modurate Detector）してＩ成分とＱ成分に変換した
ＩＱデータ３００ａを初期同期部７４０と全シンボル復
調／パラメータ推定部４００とに供給する。初期同期部
７４０は、このＩＱデータ３００ａを受けて、キャリア
中心周波数Ωf0と、スロットの同期アドレスdly0を求め
て次段に供給する。全シンボル復調／パラメータ推定部
４００は、この信号を受けて、内部で複数回推定演算し
て、変調誤差εが最小となる最適パラメータＡx、θx、
Ωx、σxを求め、かつ、全シンボルデータを復調して次
段に供給する。変調精度演算部９００は、このデータを
受けて、変調精度誤差εや各シンボル毎の変調誤差やそ
の他の情報データを求めて次段に出力する。ここで、初
期同期部７４０で求めたスロット先頭アドレスdly0によ
る変調精度誤差ε値は、本当に最小値かを確認したい場
合がある。この場合には、スロット先頭アドレスdly0前
後の位置Δdlyに変えたスロット先頭アドレス値を全シ
ンボル復調／パラメータ推定部４００に供給して変調精
度誤差εを求めて、その最小の変調精度誤差εを出力と
しても良い。出力部２３０では、これら解析結果の情報
から各種出力形式に応じたデータに加工変換した後、外
部に出力している。これらの測定全体の制御は、制御部
２８０が行っている。上記説明概要の全体の測定動作を
１サイクルとして、本測定システムの変調精度の測定
は、被測定信号１０１を間欠的にサンプリングして実施
している。

【００２３】以上で、変調精度の測定全体の動作概要を
説明した。次に、各部分の詳細説明をする。この中で、
周波数変換部１９０と、第２ミキサ２０２と、第１発振
器２０３と、第１ＬＰＦ２１０と、ＡＤ変換器２２０
と、クロック２２１と、ＡＤバッファメモリ２４０とで
構成するサンプリング部２５０と、制御部２８０の説明
は、上記の動作概要で容易に理解できるので省略する。
以下に説明するものは、ヒルベルト変換部３１０と、初
期同期部７４０と、全シンボル復調／パラメータ推定部
４００と、変調精度演算部９００と、出力部２３０であ
る。

【００２４】最初に、ヒルベルト変換部３１０の説明を
する。ヒルベルト変換部３１０の内部構成は、図２
（ａ）に示すように、バースト検出部３０２と、ＢＰＦ
３０３と、ヒルベルト変換器３０４と、ＩＱ保存メモリ
３００とで構成している。

【００２５】バースト検出部３０２は、解析に必要な連
続したスロットのデータ部分を取り出して次段に出力す
るものである。ＡＤバッファメモリ２４０に格納されて
いる２５２ワード／シンボル分解能のＡＤメモリデータ
２４０ａを読み出し、解析に必要な連続した２．４スロ
ット時間以上の有効データとして出力する。ＴＤＭＡ通
信信号の場合には、６スロット時間毎の時分割の間欠信
号である為に、９スロット時間のデータの中から解析に
必要な有効スロットデータ部分を検出して出力する。Ｂ
ＰＦ３０３は、上記有効スロットデータを受けて不用な
周波数成分を除去する為帯域フィルタした後、ヒルベル
ト変換器３０４に供給する。

【００２６】ヒルベルト変換器３０４の内部構成は、図
２（ｂ）に示すように、第１周波数３１１と、９０度位
相シフタ３１２と、第３ミキサ３１３と、第４ミキサ３
１４と、第３ＬＰＦ３１５と、第４ＬＰＦ３１６とで構
成している。ヒルベルト変換器３０４は、Ｉ／Ｑ情報を
取り出す為に直交検波を行う部分であり、実数部と虚数
部がヒルベルト（Hilbert）変換対である複素数の信号
に変換して出力する。この為に、図２（ｂ）に示すよう
に、第３ミキサ３１３は、ＢＰＦ３０３からの２２０Ｋ
Ｈｚの周波数信号３０３ａと、第１周波数３１１の２５
２ＫＨｚ信号ｓｉｎωｔを受けて、混合して直交検波信
号に変換して出力する。第３ＬＰＦ３１５は、この信号
を受けて高域をフィルタしたＩ成分（又はＱ成分の何れ
か）信号３１５ａとして再生出力する。こちら側の信号
３１５ａは、便宜上Ｉ成分データと仮定する。

【００２７】他方、第４ミキサ３１４は、周波数信号３
０３ａと、第１周波数３１１の２５２ＫＨｚ信号を９０
度位相シフタ３１２で９０度の位相をずらした２５２Ｋ
Ｈｚコサイン信号ｃｏｓωｔを受けて、混合して直交検
波信号に変換して出力する。第４ＬＰＦ３１６は、この
信号を受けて高域をフィルタしたＱ成分（又はＩ成分の
何れか）信号３１５ｂとして再生出力する。こちら側の
信号３１５ｂは、便宜上Ｑ成分データと仮定する。これ
らＩ成分／Ｑ成分信号は、２５２ＫＨｚ−２２０ＫＨｚ
＝３２ＫＨｚの中間周波数信号であり、４ｃｈのサブキ
ャリア情報を持つている２５２ワード／シンボルの分解
能のデータである。

【００２８】次に、ＩＱ保存メモリ３００は、上記ヒル
ベルト変換器３０４からの両出力信号を受けてメモリに
格納する。以後、初期同期部７４０と全シンボル復調／
パラメータ推定部４００が、このＩＱデータ３００ａを
読み出して使用する。

【００２９】次に、初期同期部７４０の説明をする。初
期同期部７４０では、上記ＩＱ保存メモリ３００からＩ
Ｑデータ３００ａを読みだして、キャリア中心周波数Ω
f0と、スロット先頭アドレスdly0を求めて出力する。初
期同期部７４０の内部構成は、図３（ａ）に示すよう
に、キャリア中心周波数検出部７００と、ｃｈ１のＩＱ
抽出部７１０と、スロット同期検出部７２０と、周波数
エラー検出部７６０とで構成している。

【００３０】まず、概略キャリア中心周波数Ωを求める
為に、キャリア中心周波数検出部７００では、図３
（ｂ）に示すように、複素ＦＦＴ処理７０２と、トータ
ルパワー計算部７０４と、キャリア中心周波数Ω計算部
７０５とで構成している。概略キャリア中心周波数Ωの
概算値の求め方は、ＦＦＴ−ＯＢＷ法によるキャリア周
波数の推定手段を用いて求める。複素ＦＦＴ処理７０２
は、ＩＱ保存メモリ３００からＩＱデータ３００ａを読
み出して周波数スペクトラムに変換する。次に、トータ
ルパワー計算部７０４は、この周波数成分の累積スペク
トラムを求めてトータルパワーを求める。次に、キャリ
ア中心周波数Ω計算部７０５は、トータルパワーの１／
２の点となる周波数として概略キャリア中心周波数Ωを
求め、これをｃｈ１のＩＱ抽出部７１０に供給する。

【００３１】次に、ｃｈ１のＩＱ抽出部７１０では、上
記概略キャリア中心周波数Ω信号と３２ＫＨｚのＩＱデ
ータ３００ａを受けて、ｃｈ１のＩＱ非同期再生データ
７１４ａを求める。ｃｈ１のＩＱ抽出部７１０の内部構
成は、図３（ｃ）に示すように、第５ミキサ７１２と、
ルート・ナイキスト・フィルタ７１３と、非同期再生メ
モリ７１４とで構成している。ｃｈ１のサブキャリア周
波数は、ｆ１ｃ＝Ω−６．７５ＫＨｚから求まる。これ
から第５ミキサ７１２は、この周波数ｆ１ｃと、３２Ｋ
ＨｚのＩＱデータ３００ａを受けて、ダウンコンバート
してｃｈ１のベースバンド信号に変換した後、次段のル
ート・ナイキスト・フィルタ７１３に供給する。ルート
・ナイキスト・フィルタ７１３は、この信号を受けてフ
ィルタしたｃｈ１の再生データを、非同期再生メモリ７
１４に格納保存する。このｃｈ１非同期再生データ７１
４ａは、スロット同期検出部７２０と周波数エラー検出
部７６０が読み出して使用する。

【００３２】スロット同期検出部７２０では、上記非同
期再生データ７１４ａを読みだして、この中から同期／
パイロットシンボル（ＰＳ）を検出してスロット先頭ア
ドレスdly0を求める。即ち、スロットの同期位置を検出
する。スロット先頭アドレスdly0の検出手段は、第１に
同期／パイロットシンボル（ＲＰＳ）の大きさ（Magnit
ude）が常に√１８値の円周上に配置されている点と、
第２に決められたシンボル番号位置に各同期／パイロッ
トシンボル（ＲＰＳ）が配置されている点に着目して検
出する。

【００３３】この同期／パイロットシンボル（ＰＳ）の
検出の説明は、１スロット構成が６０シンボルの下りの
基本スロット構成で、その中のサブチャンネル１の場合
の例で説明する。この場合の同期／パイロットシンボル
（ＰＳ）番号は、１、２、３、６、１４、２２、３０、
３８、４６、５４である。最初は、１スロット分の同期
／パイロットシンボル点の振幅値の分散を求める。この
為に、図７（ａ）に示すように、シンボル番号１のアド
レス値は、非同期再生メモリ７１４の開始アドレスを起
点アドレスdlyxとする。次のシンボル番号２のアドレス
値は、２５２ワード／シンボル単位であるから、起点ア
ドレスdlyx＋２５２のアドレス位置として求まる。以後
同様にして各同期／パイロットシンボル点のアドレス値
が求まる。

【００３４】次に、上記の各アドレス値のシンボルの振
幅値は、読みだしたＩ成分とＱ成分のデータから計算で
求まる。もし真の基本スロットの先頭シンボルアドレス
dly0であれば、これらの振幅値は、図７（ｂ）に示すベ
クトル値推移図の７２０ａ点のように全てが√１８上の
同一振幅値に存在することになり、このときの分散は、
最小の値を示す。逆に真の基本スロットの先頭シンボル
位置から外れている場合は、図７（ｃ）に示すベクトル
値推移図の７２０ｂ点のように√１８の位置から離散し
て存在し、このときの分散は、大きな値を示すことにな
る。このことから、真のスロットの先頭アドレスdly0を
求めるには、起点アドレス位置を増加しながら分散を求
めて、これらの分散データ群の中から分散が最小となる
アドレス位置を求めれば良いことがわかる。

【００３５】本実施例では、処理時間を軽減する為に、
２回に分けて実施している。すなわち、図７（ｄ）に示
す分散推移図ように、１回目は、アドレスを＋９ピッチ
毎にサンプルデータ間隔を間引いて粗く（Coarce）デー
タを読みだして粗分散を求める。この分散値の中から分
散が最小となるアドレス値７２０ｃが求まる。次に、図
７（ｅ）に示す分散推移図ように、２回目は、この分散
が最小となったアドレス値の前後領域７２０ｄについて
のみ、今度は細かく（Fine）分散を求める。この結果、
分散が最も極小となるアドレス値７２０ｅが得られる。
このアドレス値がスロットの先頭アドレスdly0として求
まる。これを、周波数エラー検出部７６０と、全シンボ
ル復調／パラメータ推定部４００に供給する。また、ス
ロット同期が求まった結果、図１６の標準規格の各同期
／パイロットシンボル（ＲＰＳ）毎の公称位相値から、
今まで便宜上のＩ成分データ３１５ａと便宜上Ｑ成分デ
ータ３１６ａが何れであるかが仮の状態であったが、こ
こで判断されて確定する。

【００３６】周波数エラー検出部７６０では、キャリア
中心周波数の誤差周波数値Δｆを求めた後、真のキャリ
ア中心周波数Ωf0を出力する。この為には、上記スロッ
ト先頭アドレスdly0を受け、非同期再生メモリ７１４か
ら同期／パイロットシンボル（ＰＳ）のＩ成分とＱ成分
データを読みだして、各々のシンボルの位相角を計算す
る。各シンボルの誤差位相データθ［ｉ］値は、上記計
算の位相角と、図１６に示す標準規格の各シンボル（Ｒ
ＰＳ）毎の公称位相角との差から求まる。これをプロッ
トすると、図８に示すプロット図が得られる。ここでθ
［ｉ］とは、ｉ番目のシンボルの位相角で、公称位相値
からの差分の角度とする。これらθ［ｉ］のデータは離
散している為、位相直線７６１は、分散を最小とする計
算（最小二乗法）によって求める。求めた位相直線７６
１からキャリア中心周波数Ωとの差、即ちエラー周波数
Δｆが、Δｆ＝ε［（Δθ／（ｔ［５４］−ｔ
［１］）］として求まる。これから、真のキャリア中心
周波数Ωf0は、Ω＋Δｆで求まる。これを、全シンボル
復調／パラメータ推定部４００に供給する。

【００３７】次に、全シンボル復調／パラメータ推定部
４００を説明する。全シンボル復調／パラメータ推定部
４００の内部構成は、図４に示すように、全ｃｈの同期
ＩＱ初期値復調部４０１と、初期値抽出部４０２と、最
適パラメータ推定部８００と、全シンボル復調部５００
と、理想データ推定部６００とで構成している。全シン
ボル復調／パラメータ推定部４００は、初期同期部７４
０からのスロット先頭アドレスdly0と、キャリア中心周
波数Ωf0を受け、ＩＱ保存メモリ３００からＩＱデータ
３００ａを読みだして、変調精度εが最小となるＡx、
Ωx、Ｂx、σxを求める。また、変調精度εが最小とな
る全ｃｈの全シンボルの復調データＡ［ｎ］、θ［ｎ］
を求める。また、このときの全シンボルの理想データＡ
ｒ［ｎ］、θｒ［ｎ］を特定する。これらのデータを変
調精度測定部に供給する。ここでｎは、全シンボル個数
とする。

【００３８】まず全シンボル復調／パラメータ推定部４
００の各部説明の前に、変調精度εを最小とするパラメ
ータの計算式（１）〜（２３）を以下に示す。この計算
結果の中で、計算式（１４）、（１５）、（２１）、
（２３）は、初期値抽出部４０２又は、最適パラメータ
推定部８００で使用する。

【００３９】

【数２】

【００４０】

【数３】

【００４１】

【数４】

【００４２】

【数５】

【００４３】

【数６】

【００４４】

【数７】 この数２〜数７は、シンボルの最適パラメータの計算式
（１）〜（２３）と近似計算の手順を説明している。

【００４５】まず、全ｃｈの同期ＩＱ初期値復調部４０
１を説明する。全ｃｈの同期ＩＱ初期値復調部４０１の
構成は、図５に示すように、ｃｈ１ミキサ４２１と、ｃ
ｈ２ミキサ４２２と、ｃｈ３ミキサ４２３と、ｃｈ４ミ
キサ４２４と、ルート・ナイキストフィルタ４３１〜４
３４とで構成している。ここでは、スロット先頭アドレ
スdly0と、キャリア中心周波数Ωf0を受け、ＩＱ保存メ
モリ３００からＩＱデータ３００ａを読みだして、４ｃ
ｈの各同期／パイロットシンボル（ＰＳ）のＩ成分とＱ
成分に変換した後、振幅値Ａi［ｎ］と位相角θi［ｎ］
を復調する。

【００４６】この為に、ｃｈ１では、Ωf0−６．７５Ｋ
Ｈｚの周波数と、３２ＫＨｚのＩＱデータ３００ａか
ら、ｃｈ１ミキサ４２１でダウンコンバートしてベース
バンドに変換した後ルート・ナイキストフィルタ４３１
でフィルタする。これにより得たＩ成分とＱ成分データ
から、振幅値Ａi［ｎ］と位相角θi［ｎ］を計算して初
期値抽出部４０２に供給する。同様に、ｃｈ２では、Ω
f0−２．２５ＫＨｚの周波数と、３２ＫＨｚのＩＱデー
タ３００ａから、ｃｈ２ミキサ４２２でダウンコンバー
トしてベースバンドに変換した後ルート・ナイキストフ
ィルタ４３２でフィルタする。これにより得たＩ成分と
Ｑ成分データから、振幅値Ａi［ｎ］と位相角θi［ｎ］
を計算して初期値抽出部４０２に供給する。同様に、ｃ
ｈ３では、Ωf0＋２．２５ＫＨｚの周波数と、３２ＫＨ
ｚのＩＱデータ３００ａから、ｃｈ３ミキサ４２３でダ
ウンコンバートしてベースバンドに変換した後ルート・
ナイキストフィルタ４３３でフィルタする。これにより
得たＩ成分とＱ成分データから、振幅値Ａi［ｎ］と位
相角θi［ｎ］を計算して初期値抽出部４０２に供給す
る。同様に、ｃｈ４では、Ωf0＋６．７５ＫＨｚの周波
数と、３２ＫＨｚのＩＱデータ３００ａから、ｃｈ４ミ
キサ４２４でダウンコンバートしてベースバンドに変換
した後ルート・ナイキストフィルタ４３４でフィルタす
る。これにより得たＩ成分とＱ成分データから、振幅値
Ａi［ｎ］と位相角θi［ｎ］を計算して初期値抽出部４
０２に供給する。

【００４７】次に、初期値抽出部４０２では、上記で求
めた各ｃｈの同期／パイロットシンボル（ＰＳ）の振幅
値Ａi［ｎ］と位相角θi［ｎ］を受けて、計算式（１
５）、（２３）から初期値データＡin、θin、Ωin、σ
inを求める。ここで、Ａinはゲイン定数であり、同期／
パイロットシンボル（ＰＳ）が√１８上に配置されてい
ることから求める。θinは初期位相であり、図１６のシ
ンボル番号１の公称位相値を使用する。Ωinは、初期キ
ャリア周波数であり、初期値は、初期同期部７４０から
のキャリア中心周波数Ωf0とする。σinはバースト・ド
ループ・ファクタ（Burst Droop Factor）で初期値は０
とする。これら初期値データは、変調精度εを最小とす
る計算の初回のパラメータとして最適パラメータ推定部
８００に供給する。

【００４８】次に、全シンボル復調部５００を説明す
る。全シンボル復調部５００の構成は、図６に示すよう
に、ｃｈ１ミキサ５２１と、ｃｈ２ミキサ５２２と、ｃ
ｈ３ミキサ５２３と、ｃｈ４ミキサ５２４と、ルート・
ナイキストフィルタ５３１〜５３４と、ｃｈ１復調メモ
リ７４１と、ｃｈ２復調メモリ７４２と、ｃｈ３復調メ
モリ７４３と、ｃｈ４復調メモリ７４４で構成してい
る。全シンボル復調部５００は、スロット先頭アドレス
dly0と、最適パラメータ推定部８００からの復調周波数
関数Ωf(x)を受け、ＩＱ保存メモリ３００からＩＱデー
タ３００ａを読みだして、４ｃｈの全シンボルの振幅値
Ａ［ｎ］と位相角θ［ｎ］を復調して復調メモリ７４１
〜７４４に格納する。ここで、各ミキサに与える復調周
波数関数Ωf(x)については、初回実行時のみ、初期同期
部７４０からのキャリア中心周波数Ωf0を使用する。２
回目以後は、最適パラメータ推定部８００からの復調周
波数関数Ωf(x)を用いて復調する。

【００４９】この為に、図６に示すように、ｃｈ１で
は、ｃｈ１ミキサ５２１とルート・ナイキストフィルタ
５３１で再生して求めた振幅値Ａ［ｎ］と位相角θ
［ｎ］をｃｈ１復調メモリ７４１に格納する。同様にし
て、ｃｈ２では、ｃｈ２ミキサ５２２とルート・ナイキ
ストフィルタ５３２で復調してｃｈ２復調メモリ７４２
に格納する。ｃｈ３では、ｃｈ３ミキサ５２３とルート
・ナイキストフィルタ５３３で復調してｃｈ３復調メモ
リ７４３に格納する。ｃｈ４では、ｃｈ４ミキサ５２４
とルート・ナイキストフィルタ５３４で復調してｃｈ４
復調メモリ７４４に格納する。そして、これら復調メモ
リの振幅値Ａ［ｎ］と位相角θ［ｎ］データは、読み出
されて最適パラメータ推定部８００と理想データ推定部
６００と変調精度演算部９００で使用される。

【００５０】次に、理想データ推定部６００を説明す
る。理想データ推定部６００では、上記の全シンボル復
調部５００で求めた振幅値Ａ［ｎ］と位相角θ［ｎ］を
使って、基準のシンボル点のどれに対応するかを判定し
又は与える。つまり、情報シンボル（ＩＳ）の場合に
は、図１２に示す１６種の基準のシンボル（ＲＩＳ）点
位置の何れかに最適判定をして判定する。一方、同期／
パイロットシンボル（ＰＳ）の場合には、Ａｒ［ｎ］＝
√１８の公称値を与え、また、θ［ｎ］は、各シンボル
番号に対応した同期／パイロットシンボル（ＲＰＳ）の
公称位相値を与える。これらで得た理想シンボルデータ
Ａｒ［ｎ］、θｒ［ｎ］は、最適パラメータ推定部８０
０と、変調精度演算部９００に供給する。ここで、この
情報シンボル（ＩＳ）の最適判定手段は、多値の変調信
号の復調として知られている最尤（さいゆう）判定方法
を使用して判定する。

【００５１】次に、最適パラメータ推定部８００を説明
する。最適パラメータを推定する為には、複数回推定し
て集束させる。最小の変調精度ε値に集束する為には、
最適パラメータ推定を３回程度繰り返すことで集束す
る。１回目の最適パラメータの推定では、初期値として
初期値抽出部４０２からの初期値データＡin、θin、Ω
in、Ｂin、σinを使う。これと理想データ推定部６００
からの理想データＡr［ｎ］、θr［ｎ］と、全シンボル
復調部５００からの全ｃｈの全シンボルの復調データ７
４１ａ〜７４４ａを受けて、変調精度εが最小となる最
適パラメータＡx、θx、Ωx、Ｂx、σxを計算式（１
４）、（１５）、（２１）から求め、また、この最適パ
ラメータから復調周波数関数Ωf(x)を求めて、全シンボ
ル復調部５００に供給する。ここで、復調周波数関数Ω
f(x)は、ゲイン定数Ａxと、初期位相θxと、キャリア周
波数Ωxと、バースト・ドループ・ファクタσxのパラメ
ータを関数とする３２ＫＨｚ相当の周波数信号である。

【００５２】２回目以後の最適パラメータを推定する場
合は、初期値として前回自身が求めた最適パラメータＡ
x、θx、Ωx、Ｂx、σxを使う。これと、前回の復調周
波数関数Ωf(x)を全シンボル復調部５００に供給して復
調した新たな復調データＡ［ｎ］、θ［ｎ］を受け、ま
た、新たに理想データ推定部６００に供給して特定した
新たな理想データＡr［ｎ］、θr［ｎ］を受けて、前回
と同様にして、変調精度εが最小となる新たな最適パラ
メータＡx、θx、Ωx、Ｂx、σxを計算式（１４）、
（１５）、（２１）から求める。また、この最適パラメ
ータから復調周波数関数Ωf(x)を求めて、全シンボル復
調部５００に再び供給する。

【００５３】この推定を、通常３回繰り返すことで変調
精度εの値は集束する。この結果、変調精度εが最小と
なる最適パラメータＡx、θx、Ωx、Ｂx、σxが最終的
に求まる。これらの最適パラメータを次段の変調精度演
算部９００に供給する。この最適パラメータを求める手
段は、特願平５ー１２３０５のＤＱＰＳＫ変調信号評価
装置の線型回帰法の手段を応用して求めている。

【００５４】次に、変調精度演算部９００を説明する。
変調精度演算部９００は、上記理想データ推定部６００
で特定した理想データＡr［ｎ］、θr［ｎ］と、全シン
ボル復調部５００からの各ｃｈの全シンボルの復調デー
タＡ［ｎ］、θ［ｎ］を受けて、変調誤差計算式３１に
基づき変調精度εを求め、また、各シンボル毎の誤差を
求めて出力部２３０に供給する。ここで、計算式３１の
パラメータの記号は、Ｖ0＝Ｂx、α＝Ａ0、φ＝θxとお
く。

【００５５】また、必要により、図１に示すように、ス
ロット先頭アドレスdly0にΔdlyを加算して前後の±１
アドレス位置にスロット先頭位置を変えて、全シンボル
復調／パラメータ推定部４００に供給して再度復調す
る。そして、変調精度εが最小となるデータを選択して
出力部２３０への出力データにしても良い。

【００５６】最後に、出力部２３０の説明をする。出力
部２３０では、これら解析結果の情報から各種出力形式
に応じたデータに加工変換した後、外部に出力してい
る。これらの測定全体の制御は、制御部２８０が行って
いる。出力部２３０は、変調精度εの値を外部に出力す
る。また、必要により、理想データＡr［ｎ］、θr
［ｎ］や、全シンボルの復調データＡ［ｎ］、θ［ｎ］
の情報データから、各種出力形式に応じたデータに加工
変換した後出力しても良い。例えば、ベクトル誤差の振
幅・位相波形、ベクトル誤差のＦＦＴ、アイパターン、
コンスタレーション等の画面表示形式のデータ等であ
る。

【００５７】（実施例２）本発明の実施例は、周波数変
換部１９０と表示手段をスペクトラムアナライザで実現
した場合である。これについて、図９を参照して以下に
説明する。システム構成は、スペクトラムアナライザ１
１０と、測定部２００と、ＳＰＡーＩ／Ｆ２３５と、外
部インターフェィス部１１３と、表示制御部１１５とで
構成している。接続としては、ＤＵＴ１００の出力端子
をスペクトラムアナライザ１１０の周波数変換部１１１
の入力端子に接続し、周波数変換部１１１の出力端子を
測定部２００の入力端子に接続し、測定部２００の出力
端子をＳＰＡーＩ／Ｆ２３５の入力端子に接続し、ＳＰ
ＡーＩ／Ｆ２３５の出力端子をスペクトラムアナライザ
１１０の外部インターフェィス部１１３の入力端子に接
続し、外部インターフェィス部１１３の出力端を表示制
御部１１５の入力端に接続している。

【００５８】スペクトラムアナライザ１１０は、ＤＵＴ
１００からの信号をスペクトラムアナライザの入力端子
で受けて、スペクトラムアナライザ内蔵の周波数変換部
１１１の機能により所定の周波数ｆ２０２＝２１．４２
ＭＨｚにダウンコンバートした後、この信号を出力端子
から出力して測定部２００に供給する。測定部２００
は、実施例１と同様の構成であり、変調精度εデータ
や、その他の測定データを出力部２３０から出力してＳ
ＰＡーＩ／Ｆ２３５に供給する。ＳＰＡーＩ／Ｆ２３５
は、これを受けて、スペクトラムアナライザ１１０の外
部インターフェィス部１１３に対応する信号に変換した
後外部インターフェィス部１１３に供給する。スペクト
ラムアナライザ１１０側の外部インターフェィス部１１
３は、ＳＰＡーＩ／Ｆ２３５からの信号を受けて、表示
制御部１１５に供給する。表示制御部１１５は、この信
号を受けてＤＭＣＡ変調精度の各種表示用に変換した後
表示画面１１６に各種表示形態に応じて表示する。これ
により、一連の測定と表示を実現している。

【００５９】（応用例）上記実施例１の説明（１２ｃ）
では、スロット同期検出部７２０のスロット同期点を検
出する手段において、高速にスロット同期点を求める為
に、粗分散と細かい分散の２回に分けて実施していた。
これに対して、全データに対して間引かずに直接１回で
分散を求めて、この分散が最も極小となるアドレス値、
即ち、スロットの先頭アドレスdly0として求める手段と
しても良い。

【００６０】

【発明の効果】本発明は、以上説明したように構成され
ているので、下記に記載されるような効果を奏する。変
調精度εに集束する迄の繰り返しの最適パラメータ推定
回数は、予め全ｃｈの同期ＩＱ初期値復調部４０１と初
期値抽出部４０２を最適パラメータ推定部８００に、概
略の初期値として最適パラメータＡin、θin、Ωin、Ｂ
in、σinを供給して計算することで軽減される。この結
果、測定解析のサイクル時間は、短時間で解析できる効
果がある。また、初期同期部７４０で、ｃｈ１の同期／
パイロットシンボルのデータから、より誤差の少ないキ
ャリア中心周波数Ωf0を求めることにより、以後の計算
時間短縮の向上をもたらす。また、最適パラメータ推定
計算の計算時間は、最適パラメータの初期値Ａin、θi
n、Ωin、Ｂin、σinを与えて、計算式（１４）（１
５）（２１）（２３）で近似計算することにより、短縮
できる効果が得られる。ＡＤ変換器２２０でサンプリン
グした大量のデータは、直接測定部２００内部で処理
し、解析した結果を外部に出力すれば良い利点が得られ
る。

【００６１】実施例２の場合は、スペクトラムアナライ
ザ１１０が標準で装備している周波数変換機能を利用
し、これに、表示制御部１１５と、ＳＰＡーＩ／Ｆ２３
５を設けることで容易に小型安価にシステムを構成でき
る利点が得られる。

【００６２】

【図面の簡単な説明】

【図１】本発明の実施例１の測定部２００のデジタル変
調解析装置である。

【図２】（ａ）本発明のヒルベルト変換部３１０の構成
図である。 （ｂ）本発明のヒルベルト変換部３１０のヒルベルト変
換器３０４の構成図である。

【図３】（ａ）本発明の初期同期部７４０の構成図であ
る。 （ｂ）本発明のキャリア中心周波数検出部７００の構成
図である。 （ｃ）本発明のｃｈ１のＩＱ抽出部７１０の構成図であ
る。

【図４】本発明の全シンボル復調／パラメータ推定部４
００の構成図である。

【図５】本発明の全ｃｈの同期ＩＱ初期値復調部４０１
の構成図である。

【図６】本発明の全シンボル復調部５００の構成図であ
る。

【図７】（ａ）非同期再生メモリ７１４から起点アドレ
スと同期／パイロットシンボル（ＰＳ）のアドレスを説
明する図である。 （ｂ）基本スロットに一致した場合のベクトル値推移図
である。 （ｃ）基本スロットから外れた場合のベクトル値推移図
である。 （ｄ）１回目の粗分散を求める場合を説明する分散推移
図である。 （ｅ）２回目の狭い領域を細かく分散を求める場合を説
明する分散推移図である。

【図８】各シンボルの公称位相値を減算した位相データ
をプロットした図である。

【図９】本発明の実施例２で、スペクトラムアナライザ
で構成した場合のデジタル変調解析装置である。

【図１０】従来の測定システムの構成図である。

【図１１】従来の測定システムの解析ソフト１７０の構
成ブロック図である。

【図１２】Ｉ成分とＱ成分による１６のシンボル点と円
周上の同期／パイロットシンボル（ＰＳ）を説明する図
である。

【図１３】観測ベクトルの変調誤差を説明する図であ
る。

【図１４】４ｃｈサブキャリア周波数を説明する図と、
中心周波数Ωと、ＩとＱ成分のビットストリームとを説
明する図である。

【図１５】ＤＭＣＡの４ｃｈサブキャリアの送信系統図
である。

【図１６】標準規格のスロット構成と各サブチャンネル
毎の各同期／パイロットシンボルの公称位相値を示す図
である。

【図１７】標準規格のスロット構成と各サブチャンネル
毎のシンボル配置を示す図である。

【符号の説明】

１００ ＤＵＴ（被試験器） １０１ 被測定信号 １０２ 第１ミキサ １０３ 信号発生器 １１０ スペクトラムアナライザ １１１ 周波数変換部 １１３ 外部インターフェィス部 １１５ 表示制御部 １１６ 表示画面 １７０ 解析ソフト １７２ チャンネル分離部 １７３ａ〜１７３ｄ ＩＱ変換部 １７４ａ〜１７４ｄ ルート・ナイキスト・フィルタ
部 １７５ａ〜１７５ｄ クロック推定部 １７６ａ〜１７６ｄ キャリア・オフセット推定部 １７７ 変調精度誤差演算部 １８３ 信号発生器 １８４ ワークステーション １９０ デジタル変換部 １９１ 被測定信号 １９２ 周波数変換器 １９３、２１０ ＬＰＦ（ローパス・フィルタ） １９４、２２０ ＡＤ変換器 １９５ バッファメモリ １９６ 制御部 １９７ ＧＰＩＢインターフェイス ２００ 測定部 ２０２ 第２ミキサ ２０３ 第１発振器 ２１０ 第１ＬＰＦ ２２１ クロック ２３０ 出力部 ２３５ ＳＰＡーＩ／Ｆ ２４０ ＡＤバッファメモリ ２４０ａ ＡＤメモリデータ ２５０ サンプリング部 ２８０ 制御部 ３００ ＩＱ保存メモリ ３００ａ ＩＱデータ ３０２ バースト検出部 ３０３ ＢＰＦ ３０４、３１０ ヒルベルト変換部 ３１１ 第１周波数 ３１２ ９０度位相シフタ ３１３ 第３ミキサ ３１４ 第４ミキサ ３１５ 第３ＬＰＦ ３１６ 第４ＬＰＦ ４００ 全シンボル復調／パラメータ推定部
４００ ４０１ 全ｃｈの同期ＩＱ初期値復調部４０
１ ４０２ 初期値抽出部 ４２１、５２１ ｃｈ１ミキサ ４２２、５２２ ｃｈ２ミキサ ４２３、５２３ ｃｈ３ミキサ ４２４、５２４ ｃｈ４ミキサ ４３１、４３２、４３３、４３４ ルート・ナイキス
トフィルタ ５３１、５３２、５３３、５３４、７１３ ルート・
ナイキストフィルタ ５００ 全シンボル復調部 ６００ 理想データ推定部 ７００ キャリア中心周波数検出部 ７０２ 複素ＦＦＴ処理 ７０４ トータルパワー計算部 ７０５ キャリア中心周波数Ω計算部 ７１０ ｃｈ１のＩＱ抽出部 ７１２ 第５ミキサ ７１４ 非同期再生メモリ ７２０ スロット同期検出部 ７４０ 初期同期部 ７４１ ｃｈ１復調メモリ ７４２ ｃｈ２復調メモリ ７４３ ｃｈ３復調メモリ ７４４ ｃｈ４復調メモリ ７６０ 周波数エラー検出部 ８００ 最適パラメータ推定部 ９００ 変調精度演算部 Ωf0 キャリア中心周波数 dly0 スロット先頭アドレス Ａin 初期値ゲイン定数 θin 初期値位相角 Ωin 初期値キャリア周波数 Ｂin 初期値原点オフセット σin 初期値バースト・ドループ・ファク
タ Ａi［ｎ］ ｃｈ１の振幅値（ｎ＝シンボル番
号） θi［ｎ］ ｃｈ１の位相角値（ｎ＝シンボル番
号） Ａｒ［ｎ］ 理想シンボルの振幅値（ｎ＝シンボ
ル番号） θｒ［ｎ］ 理想シンボルの位相角値（ｎ＝シン
ボル番号） Ａ［ｎ］ 測定振幅値（ｎ＝シンボル番号） θ［ｎ］ 測定位相角（ｎ＝シンボル番号） Ａx ゲイン定数 θx 初期位相 Ωx キャリア周波数 Ｂx 原点オフセット σx バースト・ドループ・ファクタ

Claims (2)

【特許請求の範囲】
1. 【請求項１】 デジタル通信システムで被試験器（ＤＵ
   Ｔ）が出力するデジタル方式のＭＣＡ（Multi Channel
   Acces）信号の変調精度の測定において、 ＤＵＴの信号を受けて、一定の中間周波数に変換する周
   波数変換部（１９０）を設け、 周波数変換部（１９０）からの信号を受けて、さらに低
   い中間周波数信号に変換した後ＡＤ変換してメモリに保
   存するサンプリング部（２５０）を設け、 直交検波（Quadrant Modurate Detector）してＩ成分と
   Ｑ成分に変換するヒルベルト変換部（３１０）を設け、 Ｉ成分とＱ成分信号を受けて、キャリア中心周波数（Ω
   f0）と、スロットの同期アドレス（dly0）を求める初期
   同期部（７４０）を設け、 初期同期部（７４０）の信号とＩ成分とＱ成分信号を受
   けて、変調誤差εが最小となる最適パラメータ（ゲイン
   定数Ａx、初期位相θx、キャリア周波数Ωx、バースト
   ・ドループ・ファクタσx）を求め、かつ、全シンボル
   の復調データ（シンボル番号ｎのベクトル値Ａ［ｎ］、
   位相角θ［ｎ］）を復調し、全シンボルの理想データ
   （理想シンボルのベクトル値Ａｒ［ｎ］、理想シンボル
   の位相角値θｒ［ｎ］）を特定する全シンボル復調／パ
   ラメータ推定部（４００）を設け、 全シンボル復調／パラメータ推定部（４００）の信号を
   受けて、変調精度εデータや各シンボル毎の変調誤差デ
   ータを求める変調精度演算部（９００）を設け、 以上を具備していることを特徴としたデジタル変調解析
   装置。
2. 【請求項２】 デジタル通信システムで被試験器が出力
   するデジタル方式のＭＣＡ信号の変調精度の測定におい
   て、 スペクトラムアナライザ（１１０）内に、ＤＵＴの信号
   を受けて、一定の中間周波数信号に変換する周波数変換
   部（１１１）を設け、測定部（２００）からの信号を受
   けて、変調誤差の画面表示用データに変換する表示制御
   部（１１５）を設け、 測定部（２００）としては、請求項１記載の、サンプリ
   ング部（２５０）と、ヒルベルト変換部（３１０）と、
   初期同期部（７４０）と、全シンボル復調／パラメータ
   推定部（４００）と、変調精度演算部（９００）とを設
   けて構成し、 以上を具備していることを特徴としたデジタル変調解析
   装置。