JP2002014600A - 相関関数測定方法及び装置 - Google Patents

相関関数測定方法及び装置

Info

Publication number
JP2002014600A
JP2002014600A JP2000296812A JP2000296812A JP2002014600A JP 2002014600 A JP2002014600 A JP 2002014600A JP 2000296812 A JP2000296812 A JP 2000296812A JP 2000296812 A JP2000296812 A JP 2000296812A JP 2002014600 A JP2002014600 A JP 2002014600A
Authority
JP
Japan
Prior art keywords
calibration
correlation function
calibration data
frequency
signal processing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2000296812A
Other languages
English (en)
Other versions
JP4038009B2 (ja
Inventor
Makoto Kurosawa
誠 黒澤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Advantest Corp
Original Assignee
Advantest Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Advantest Corp filed Critical Advantest Corp
Priority to JP2000296812A priority Critical patent/JP4038009B2/ja
Priority to US09/771,527 priority patent/US6594605B2/en
Priority to DE10165055A priority patent/DE10165055B4/de
Priority to DE10103822A priority patent/DE10103822B4/de
Priority to GB0102749A priority patent/GB2366030B/en
Publication of JP2002014600A publication Critical patent/JP2002014600A/ja
Application granted granted Critical
Publication of JP4038009B2 publication Critical patent/JP4038009B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Complex Calculations (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Holo Graphy (AREA)

Abstract

(57)【要約】 【課題】 高精度に相関ベクトルを校正することができ
る相関関数測定方法及び装置を提供する。 【解決手段】 第１の観測信号を処理して第１のスペク
トルを求める第１の信号処理手段１８ａ〜２２ａと、第
２の観測信号を処理して第２のスペクトルを求める第２
の信号処理手段１８ｂ〜２２ｂと、第１の信号処理手段
及び第２の信号処理手段に校正信号を入力し、第１の信
号処理手段により求められたスペクトルと、第２の信号
処理手段により求められたスペクトルとから、各周波数
に応じた校正値を求める校正値演算手段２６と、校正値
演算手段により求められた校正値を用いつつ、第１のス
ペクトルと第２のスペクトルとの相関関数を求める相関
関数演算手段３０とを有している。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、相関関数測定方法
及び装置に係り、特に、複数の入力信号の相関ベクトル
を測定する相関関数測定方法及び装置に関する。
【０００２】
【従来の技術】近年、電波ホログラフィ技術を用いて波
源像を可視化する技術が実用化され、例えば、携帯電話
の基地局からの電波を観測するために用いられている。
電波ホログラムを解析することにより、基地局からの電
波を視覚的に観測することができ、新たに基地局を設置
すべきか否かを判断する有効な材料を得ることができ
る。
【０００３】電波ホログラフィ技術を用いた波源像可視
化方法及び装置については、本願出願人が、例えば特願
平７−２８９８４８号明細書で提案している。
【０００４】本願出願人が提案した波源像可視化方法及
び装置を図８を用いて説明する。図８は、提案されてい
る波源像可視化装置を示すブロック図である。
【０００５】図８に示すように、固定アンテナ（図示せ
ず）により受信された観測信号ａは、前変換部２１８
ａ、データメモリ２２０ａ、フーリエ変換部２２２ａに
より所定の信号処理が為されて、相関ベクトル演算部２
３０に入力されるようになっており、走査アンテナ（図
示せず）により受信された観測信号ｂは、前変換部２１
８ｂ、データメモリ２２０ｂ、フーリエ変換部２２２ｂ
により所定の信号処理が為されて、相関ベクトル演算部
２３０に入力されるようになっている。前置換部２１８
ａ、２１８ｂは、観測信号ａ、ｂを中心周波数ｆ0、帯
域幅ｂｗで帯域制限し、ＩＦ信号に変換して出力するも
のであり、データメモリ２２０ａ、２２０ｂは、ＩＦ信
号をＡ／Ｄ変換して記憶するものであり、フーリエ変換
部２２２ａ、２２２ｂは、フーリエ変換を行ってスペク
トルＳA（ｆ）、ＳB（ｆ）を出力するものである。
【０００６】そして、相関ベクトル演算部２３０では、
以下の式
【０００７】
【数１】
【０００８】に基づいて、相関ベクトルが求められる。
なお、＊は、複素共役を示している。
【０００９】このような波源像可視化装置では、前変換
部２１８ａ、２１８ｂにおいて帯域制限や周波数変換等
の測定条件を変更すると、前変換部２１８ａ、２１８ｂ
の周波数特性が変化する。しかも、この周波数特性の変
化は、前変換部２１８ａと前置換部２１８ｂとの相互間
でもばらつく。
【００１０】このため、予め校正データを取得しておい
て、その校正データで相関ベクトルを校正することが必
要である。そこで、図８に示す波源像可視化装置では、
以下のようにして校正データρ（ｃ）を求める。
【００１１】即ち、校正データρ（ｃ）を求める際に
は、スイッチ２１６ａ、２１６ｂ、２２４を校正側ｃに
設定し、信号発生器２１２から校正信号を出力して、相
関ベクトル演算部２３０で相関ベクトルを求める。そし
て、求められた相関ベクトルの値を、校正データρ
（ｃ）として、校正データメモリ２２８に記憶する。
【００１２】そして、電波監視を行う際には、スイッチ
２１６ａ、２１６ｂ、２２４を測定側ｍに設定し、相関
ベクトル演算部２３０で相関ベクトルを求める。求めら
れた相関ベクトルは、校正部２２９において校正データ
ρ（ｃ）を用いて校正され、波源像再生処理部（図示せ
ず）に出力される。
【００１３】波源像再生処理部は、所定の画像処理を行
い、処理された画像は、表示部（図示せず）のディスプ
レイに表示される。
【００１４】
【発明が解決しようとする課題】しかしながら、提案さ
れている波源像可視化装置では、校正信号を入力した際
の帯域周波数平均の相関ベクトルを校正データρ（ｃ）
として用いていた。このため、前変換部２１８ａ、２１
８ｂの特性が周波数によって変化する場合には、周波数
によって校正精度がばらついてしまっていた。
【００１５】また、前変換部２１８ａ、２１８ｂの特性
は、測定条件の変更、周囲温度の変化、経時変化等によ
りばらつくものである。提案されている波源像可視化装
置では、これらの要因により校正精度が低下してしまっ
ていた。
【００１６】本発明の目的は、高精度に相関ベクトルを
校正することができる相関関数測定方法及び装置を提供
することにある。
【００１７】
【課題を解決するための手段】上記目的は、第１の観測
信号を第１の信号処理手段により処理して第１のスペク
トルを求め、第２の観測信号を第２の信号処理手段によ
り処理して第２のスペクトルを求め、前記第１のスペク
トルと前記第２のスペクトルとの相関関数を測定する相
関関数測定方法であって、相関関数の測定に先立って、
前記第１の信号処理手段及び前記第２の信号処理手段に
校正信号を入力し、前記第１の信号処理手段により求め
られたスペクトルと、前記第２の信号処理手段により求
められたスペクトルとから、前記校正信号の周波数に応
じた校正値を求め、相関関数の測定の際には、周波数に
応じた前記校正値を用いつつ、前記第１のスペクトルと
前記第２のスペクトルとの相関関数を求めることを特徴
とする相関関数測定方法により達成される。これによ
り、相関関数の測定に先立って、各周波数毎に校正デー
タを取得しておき、周波数に応じた校正データを用いて
相関関数を求めるので、相関関数を高い精度で求めるこ
とができる。従って、高精度に相関関数を測定すること
ができる相関関数測定方法を提供することができる。
【００１８】また、上記の相関関数測定方法において、
前記校正値を求める際には、前記校正信号の周波数を掃
引することにより、各周波数に応じた前記校正値を求め
ることが望ましい。
【００１９】また、上記の相関関数測定方法において、
前記校正値から回帰式を求め、前記校正値を更新する際
には、少なくとも前記回帰式を決定するのに必要な校正
値を再度求めることにより前記回帰式を更新し、更新さ
れた前記回帰式に基づいて前記校正値を更新することが
望ましい。
【００２０】また、上記の相関関数測定方法において、
前記校正値を求める際には、前記第１の信号処理手段及
び前記第２の信号処理手段に変調波を入力することによ
り、各周波数に応じた校正値を求めることが望ましい。
【００２１】また、上記の相関関数測定方法において、
前記校正値を求める際には、前記第１の信号処理手段及
び前記第２の信号処理手段に、第１の校正信号を周波数
を掃引しつつ入力し、前記第１の信号処理手段により求
められたスペクトルと、前記第２の信号処理手段により
求められたスペクトルとから、各周波数に応じた振幅校
正データを求め、前記第１の信号処理手段及び前記第２
の信号処理手段に第２の校正信号を入力し、前記第１の
信号処理手段により求められたスペクトルと、前記第２
の信号処理手段により求められたスペクトルとから、各
周波数に応じた位相校正データを求め、前記振幅校正デ
ータと前記位相校正データとから各周波数に応じた前記
校正値を求めることが望ましい。これにより、頻繁に取
得する必要のない振幅校正データを第１の校正信号を用
いて取得し、頻繁に取得する必要のある位相校正データ
を第２の校正信号を用いて取得し、これらのデータを用
いて各周波数に応じた校正値を求めるので、校正精度を
劣化することなく、迅速に校正を行うことができる。
【００２２】また、上記の相関関数測定方法において、
前記第２の校正信号は、変調波より成ることが望まし
い。
【００２３】また、上記の相関関数測定方法において、
前記振幅校正データを求める際には、前記第１の校正信
号の周波数を所定のステップで掃引し、補間を行うこと
により前記振幅校正データを求めることが望ましい。
【００２４】また、上記の相関関数測定方法において、
前記位相校正データを、前記振幅校正データより高い頻
度で更新することが望ましい。
【００２５】また、上記の相関関数測定方法において、
前記振幅校正データを更新した際の前記位相校正データ
に対して、更新された前記位相校正データが所定値以上
変化している場合に、前記振幅校正データを更新するこ
とが望ましい。
【００２６】また、上記目的は、第１の観測信号を処理
して第１のスペクトルを求める第１の信号処理手段と、
第２の観測信号を処理して第２のスペクトルを求める第
２の信号処理手段と、前記第１の信号処理手段及び前記
第２の信号処理手段に校正信号を入力し、前記第１の信
号処理手段により求められたスペクトルと、前記第２の
信号処理手段により求められたスペクトルとから、各周
波数に応じた校正値を求める校正値演算手段と、前記校
正値演算手段により求められた前記校正値を用いつつ、
前記第１のスペクトルと前記第２のスペクトルとの相関
関数を求める相関関数演算手段とを有することを特徴と
する相関関数測定装置により達成される。これにより、
相関関数の測定に先立って、各周波数毎に校正データを
取得しておき、周波数に応じた校正データを用いて相関
関数を求めるので、相関関数を高い精度で求めることが
できる。従って、高精度に相関関数を測定することがで
きる相関関数測定装置を提供することができる。
【００２７】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記校正信号の周波数を掃引す
ることにより、各周波数に応じた前記校正値を求めるこ
とが望ましい。
【００２８】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記校正値から回帰式を求め、
少なくとも前記回帰式を決定するのに必要な校正値を再
度求めることにより前記回帰式を更新し、更新された前
記回帰式に基づいて前記校正値を更新することが望まし
い。
【００２９】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記第１の信号処理手段及び前
記第２の信号処理手段に変調波を入力することにより、
各周波数に応じた前記校正値を求めることが望ましい。
【００３０】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記第１の信号処理手段及び前
記第２の信号処理手段に、第１の校正信号を周波数を掃
引しつつ入力し、前記第１の信号処理手段により求めら
れたスペクトルと、前記第２の信号処理手段により求め
られたスペクトルとから、各周波数に応じた振幅校正デ
ータを求め、前記第１の信号処理手段及び前記第２の信
号処理手段に第２の校正信号を入力し、前記第１の信号
処理手段により求められたスペクトルと、前記第２の信
号処理手段により求められたスペクトルとから、各周波
数に応じた位相校正データを求め、前記振幅校正データ
と前記位相校正データとから各周波数に応じた前記校正
値を求めることが望ましい。これにより、頻繁に取得す
る必要のない振幅校正データを第１の校正信号を用いて
取得し、頻繁に取得する必要のある位相校正データを第
２の校正信号を用いて取得し、これらのデータを用いて
各周波数に応じた校正値を求めるので、校正精度を劣化
することなく、迅速に校正を行うことができる。
【００３１】また、上記の相関関数測定装置において、
前記第２の校正信号は、変調波より成ることが望まし
い。
【００３２】また、上記の相関関数測定装置において、
校正値演算手段は、前記第１の校正信号の周波数を所定
のステップで掃引し、補間を行うことにより前記振幅校
正データを求めることが望ましい。
【００３３】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記位相校正データを、前記振
幅校正データより高い頻度で更新することが望ましい。
【００３４】また、上記の相関関数測定装置において、
前記校正値演算手段は、前記振幅校正データを更新した
際の前記位相校正データに対して、更新された前記位相
校正データが所定値以上変化している場合に、前記振幅
校正データを更新することが望ましい。
【００３５】
【発明の実施の形態】［第１実施形態］本発明の第１実
施形態による相関関数測定方法及び装置を図１を用いて
説明する。図１は、本実施形態による相関関数測定装置
を示すブロック図である。
【００３６】なお、本実施形態では、電波監視を行う際
に固定アンテナ及び走査アンテナから入力される観測信
号の相関関数を測定する場合を例に説明するが、本発明
の原理は、電波監視を行う場合のみならず、複数の観測
信号の相関関数を測定する際に広く適用することが可能
である。
【００３７】（相関関数測定装置）図１に示すように、
本実施形態による相関関数測定装置は、校正信号の周波
数を制御する校正周波数制御部１０と、一定の振幅で校
正信号を発生する信号発生器１２と、校正信号を分配す
るパワースプリッタ１４と、校正時と測定時とで回路を
切り換えるスイッチ１６ａ、１６ｂと、所定の周波数成
分のみを帯域制限し、中間周波数ＩＦ（Intermediate-F
requency）に変換して出力する前変換部１８ａ、１８ｂ
と、前変換部１８ａ、１８ｂから出力されるＩＦ信号を
サンプリングして記憶するデータメモリ２０ａ、２０ｂ
と、データメモリ２０ａ、２０ｂに記憶されたデータを
フーリエ変換してスペクトルＳA（ｆ）、ＳB（ｆ）を求
めるフーリエ変換部２２ａ、２２ｂと、校正時と測定時
とで回路を切り換えるスイッチ２４ａ、２４ｂと、各校
正周波数に応じたスペクトルから校正データＹ（ｆC）
を求める校正データ演算部２６と、各校正周波数毎の校
正データを記憶する校正データメモリ２８と、校正デー
タを参照しつつ観測信号のスペクトルの相関ベクトルを
求める相関ベクトル演算部３０とを有している。
【００３８】スイッチ１６ａの校正側ｃはパワースプリ
ッタ１４に接続されており、スイッチ１６ａの測定側ｍ
は固定アンテナに接続されている。また、スイッチ１６
ｂの校正側ｃはパワースプリッタ１４に接続されてお
り、スイッチ１６ｂの測定側ｍは走査アンテナに接続さ
れている。また、スイッチ２４ａ、２４ｂの校正側ｃは
校正データ演算部２６に接続されており、スイッチ２４
ａ、２４ｂの測定側ｍは相関ベクトル演算部３０に接続
されている。このため、スイッチ１６ａ、１６ｂ、２４
ａ、２４ｂを校正側ｃに設定すると、校正データを取得
する回路が構成され、スイッチ１６ａ、１６ｂ、２４
ａ、２４ｂを測定側ｍに設定すると、相関関数を測定す
る回路が構成される。
【００３９】本実施形態による相関関数測定装置は、相
関関数の測定に先立って、各周波数毎の校正データを取
得しておき、相関関数の測定を行う際に、各周波数毎に
校正データを参照しつつ相関ベクトルを求めることに主
な特徴がある。
【００４０】図８に示す提案されている波源像可視化装
置では、一定の周波数帯域に対して、帯域内周波数平均
の校正データρ（ｃ）を１つだけ求め、この校正データ
ρ（ｃ）を用いて校正を行っていたため、当該周波数帯
域において不均一なスペクトル分布が存在する場合に
は、正確な相関ベクトルを求めることができなかった。
【００４１】これに対し、本実施形態では、相関関数の
測定に先立って、各周波数毎に予め校正データを求めて
おき、各周波数毎の校正データを参照しつつ相関ベクト
ルを求めるので、高精度に相関ベクトルを求めることが
できる。
【００４２】従って、本実施形態によれば、高精度に相
関関数の測定を行うことができる相関関数測定装置を提
供することができる。
【００４３】以下、本実施形態による相関関数測定装置
の各構成部分について個々に詳述する。
【００４４】（ａ） 校正周波数制御部１０ 校正周波数制御部１０は、信号発生器１２から出力され
る校正信号の周波数を、所定のステップで掃引するよう
に制御するものである。
【００４５】具体的には、ｆ0±ｂｗ／２の周波数帯域
において、ｆS／Ｍのステップで校正周波数ｆCを掃引す
る。ここで、ｆ0は中心周波数であり、ｂｗは帯域幅で
あり、ｆSは発振器２１からデータメモリ２０ａ、２０
ｂに供給されるサンプリング周波数であり、Ｍは任意の
整数である。
【００４６】校正周波数ｆCは、信号発生器１２のみな
らず、校正データ演算部２６にも入力されるようになっ
ている。
【００４７】（ｂ） 信号発生器１２ 信号発生器１２は、校正周波数制御部１０から供給され
る校正周波数ｆCに応じて、校正信号を発生するもので
ある。
【００４８】信号発生器１２は、周波数に依存すること
なく、一定の振幅の校正信号を発生する。
【００４９】（ｃ） パワースプリッタ１４ パワースプリッタ１４は、信号発生器１２から出力され
る校正信号を２つの経路に分配するものである。
【００５０】パワースプリッタ１４の出力は、スイッチ
１６ａ、１６ｂの校正側ｃに接続されている。
【００５１】（ｄ） スイッチ１６ａ、１６ｂ スイッチ１６ａ、１６ｂは、校正時と測定時とで入力信
号を切り換えるものである。スイッチ１６ａ、１６ｂに
入力される信号は高周波であるため、高周波用のスイッ
チが用いられる。
【００５２】スイッチ１６ａの校正側ｃには、パワース
プリッタ１４が接続されており、スイッチ１６ａの測定
側ｍには、例えば固定アンテナ（図示せず）が接続され
るようになっている。
【００５３】スイッチ１６ａを校正側ｃに設定すると、
信号発生器１２からの校正信号がスイッチ１６ａを介し
て前変換部１８ａに入力されるようになり、スイッチ１
６ａを測定側ｍに設定すると、固定アンテナからの観測
信号ａがスイッチ１６ａを介して前変換部１８ａに入力
されるようになる。
【００５４】一方、スイッチ１６ｂの校正側ｃには、パ
ワースプリッタ１４が接続されており、スイッチ１６ｂ
の測定側ｍには、例えば走査アンテナ（図示せず）が接
続されるようになっている。スイッチ１６ｂを校正側ｃ
に設定すると、信号発生器１２からの校正信号がスイッ
チ１６ｂを介して前変換部１８ｂに入力されるようにな
り、スイッチ１６ｂを測定側ｍに設定すると、走査アン
テナからの観測信号ｂがスイッチ１６ｂを介して前変換
部１８ｂに入力されるようになる。
【００５５】即ち、スイッチ１６ａ、１６ｂを校正側ｃ
に設定すると、校正データを取得する回路が構成され、
スイッチ１６ａ、１６ｂを測定側ｍに設定すると、相関
関数の測定を行う回路が構成される。
【００５６】（ｅ） 前変換部１８ａ、１８ｂ スイッチ１６ａ、１６ｂの出力側には、それぞれ前変換
部１８ａ、１８ｂが設けられている。
【００５７】前変換部１８ａ、１８ｂは、所定の周波数
帯域の成分のみを通過させ、中間周波数ＩＦに変換して
出力するものである。
【００５８】前変換部１８ａ、１８ｂは、例えば、ＲＦ
スペクトラムアナライザを用いて構成することができ
る。ＲＦスペクトラムアナライザを用いて前置換部１８
ａ、１８ｂを構成する場合には、ゼロスパンモードに設
定し、基準周波数ｆREF（図示せず）で位相をロックす
る。
【００５９】（ｆ） データメモリ２０ａ、２０ｂ 前変換部１８ａ、１８ｂの出力側には、データメモリ２
０ａ、２０ｂが設けられており、データメモリ２０ａ、
２０ｂには、それぞれサンプリング周波数ｆSが供給さ
れるようになっている。
【００６０】データメモリ２０ａ、２０ｂは、前変換部
１８ａ、１８ｂから出力されるＩＦ信号をサンプリング
周波数ｆSに基づいてサンプリングし、Ａ／Ｄ変換した
後、メモリに記憶する。データメモリ２０ａ、２０ｂに
は、Ｍ個のデータが記憶される。なお、サンプリング周
波数ｆSは、ｂｗ＝ｆS／２となるように設定されてい
る。
【００６１】（ｇ） フーリエ変換部２２ａ、２２ｂ データメモリ２０ａ、２０ｂの出力側には、フーリエ変
換部２２ａ、２２ｂが設けられている。
【００６２】フーリエ変換部２２ａ、２２ｂは、データ
メモリ２０ａ、２０ｂから出力されるデータをフーリエ
変換して、スペクトルＳA（ｆ）、ＳB（ｆ）を出力する
ものである。
【００６３】（ｈ） スイッチ２４ａ、２４ｂ フーリエ変換部２２ａ、２２ｂの出力側には、スイッチ
２４ａ、２４ｂが設けられている。
【００６４】スイッチ２４ａ、２４ｂは、スイッチ１６
ａ、１６ｂと同様、校正時と測定時とで回路を切り換え
るものである。
【００６５】スイッチ２４ａの校正側ｃは、校正データ
演算部２６に接続されており、スイッチ２４ａの測定側
ｍは、相関ベクトル演算部３０に接続されている。ま
た、スイッチ２４ｂの校正側ｃは、校正データ演算部２
６に接続されており、スイッチ２４ｂの測定側ｍは、相
関ベクトル演算部３０に接続されている。
【００６６】スイッチ２４ａ、２４ｂを校正側ｃに設定
すると、フーリエ変換部２２ａ、２２ｂから出力される
スペクトルＳA（ｆC）、ＳB（ｆC）が校正データ演算部
２６に入力されるようになり、スイッチ２４ａ、２４ｂ
を測定側ｍに設定すると、フーリエ変換部２２ａ、２２
ｂから出力されるスペクトルＳA（ｆ）、ＳB（ｆ）が相
関ベクトル演算部３０に入力されるようになる。
【００６７】（ｉ） 校正データ演算部２６ 校正データ演算部２６は、各周波数毎に校正データを演
算するものである。
【００６８】校正データ演算部２６には、校正周波数制
御部１０から校正周波数ｆCが入力されるようになって
おり、スイッチ２４ａ、２４ｂの校正側ｃからスペクト
ルＳ A（ｆC）、ＳB（ｆC）が入力されるようになってい
る。
【００６９】校正データ演算部２６は、以下の式から、
校正周波数ｆC毎に校正データＹ（ｆC）を求め、校正デ
ータメモリ２８に出力する。なお、＊は、複素共役を示
している。
【００７０】
【数２】
【００７１】（ｊ） 校正データメモリ２８ 校正データメモリ２８は、校正データ演算部２６により
求められた各周波数毎の校正データＹ（ｆC）を記憶す
るものである。
【００７２】即ち、ｆ0±ｂｗ／２の周波数帯域におい
て、ｆS／Ｍのステップで校正データＹ（ｆC）が求めら
れ、これら校正データＹ（ｆC）が校正データメモリ２
８に記憶される。
【００７３】（ｋ） 相関ベクトル演算部３０ 相関ベクトル演算部３０は、相関ベクトルを演算するも
のである。
【００７４】相関ベクトル演算部３０には、フーリエ変
換部２２ａ、２２ｂから出力されたスペクトルＳ
A（ｆ）、ＳB（ｆ）が、スイッチ２４ａ、２４ｂを介し
て入力されるようになっており、相関ベクトル演算部３
０は観測信号の周波数ｆに応じて校正データメモリ２８
から校正データＹ（ｆC）を適宜参照し、以下の式に基
づいて相関ベクトルを演算する。
【００７５】
【数３】
【００７６】相関ベクトル演算部３０により求められた
相関ベクトルは、波形像再生制御部（図示せず）等に出
力され、表示部（図示せず）に表示される。
【００７７】なお、本実施形態による相関関数測定装置
は、本願出願人による特願平７−２８９８４８号明細書
に記載された波源像可視化装置等に適用可能である。波
源像可視化装置に適用する場合の構成等については、当
該明細書を参照されたい。
【００７８】（相関関数測定方法）次に、本実施形態に
よる相関関数測定方法を図１を用いて説明する。
【００７９】まず、相関関数の測定に先立って、校正デ
ータＹ（ｆC）を取得する。校正データＹ（ｆC）を取得
する際には、各スイッチ１６ａ、１６ｂ、２４ａ、２４
ｂをそれぞれ校正側ｃに設定する。これにより、信号発
生器１２から出力された校正信号が前変換部１８ａ、１
８ｂに入力されるようになり、また、フーリエ変換され
たスペクトルＳA（ｆC）、ＳB（ｆC）が校正データ演算
部２６に入力されるようになる。
【００８０】次に、校正周波数ｆCを掃引して、各校正
周波数ｆC毎に校正データＹ（ｆC）を求める。具体的に
は、校正周波数制御部１０から信号発生器１２に供給す
る校正周波数ｆCを、ｆS／Ｍのステップ、ｆ0±ｂｗ／
２の周波数帯域で順次掃引していき、各校正周波数ｆC
毎に校正データＹ（ｆC）を求めていく。こうして求め
られた校正データＹ（ｆC）は、校正データメモリ２８
に順次記憶される。こうして、校正データの取得が完了
する。
【００８１】相関関数の測定を行う際には、各スイッチ
１６ａ、１６ｂ、２４ａ、２４ｂを測定側ｍに設定す
る。これにより、固定アンテナから観測信号ａが入力さ
れるようになり、走査アンテナから観測信号ｂが入力さ
れるようになる。また、フーリエ変換部２２ａ、２２ｂ
から出力されるスペクトルＳA（ｆ）、ＳB（ｆ）が、相
関ベクトル演算部３０に入力されるようになる。
【００８２】相関ベクトルを演算する際には、校正デー
タメモリ２８に記憶された校正データＹ（ｆC）を、観
測信号の周波数に応じて適宜参照し、相関ベクトルを演
算する。こうして、観測信号ａ、ｂの周波数成分に応じ
て各周波数毎に校正が行われ、高い精度で相関ベクトル
が求められる。なお、相関関数測定方法の詳細は、本願
出願人による特願平７−２８９８４８号明細書を参照さ
れたい。
【００８３】測定条件を変更した場合、周囲温度が変化
した場合、経時変化が生じた場合など、校正データの信
憑性が低くなった場合には、再度上記と同様にして校正
データＹ（ｆC）を取得し、校正データを更新する。
【００８４】そして、新たな校正データＹ（ｆC）を参
照しつつ、相関ベクトルを演算すればよい。
【００８５】このように本実施形態によれば、相関関数
の測定に先立って、各周波数毎に校正データを取得して
おき、周波数に応じた校正データを用いて相関ベクトル
を求めるので、相関ベクトルを高い精度で求めることが
できる。また、校正データの信憑性が低くなった場合に
は、適宜校正データの更新を行うので、相関ベクトルの
校正精度を維持することができる。従って、本実施形態
によれば、高精度に相関関数の測定を行うことができる
相関関数測定方法及び装置を提供することができる。
【００８６】（変形例（その１））次に、本実施形態に
よる相関関数測定方法及び装置の変形例（その１）を図
２及び図３を用いて説明する。図２は、本変形例による
相関関数測定装置を示すブロック図である。図３は、周
波数変調と位相変調とを示すタイムチャートである。
【００８７】本変形例は、校正信号として変調波を用い
ることに主な特徴がある。変調波としては、例えば上記
の帯域幅ｂｗより広い周波数帯域で変調されたＦＭ波等
を用いることができる。
【００８８】図１に示す相関関数測定装置では、校正デ
ータを取得する際に校正信号の周波数を掃引していたた
め、校正データの取得に長時間を要していたが、本変形
例では、多様な周波数成分を含む変調波を校正信号に用
いるので、短時間で校正データを取得することができ
る。従って、本変形例によれば、校正データの取得をよ
り迅速化することができる。
【００８９】但し、変調波は、振幅が周波数に対して均
一でないため、以下のように、理論振幅分布で補正する
必要がある。
【００９０】理論振幅分布による補正について、方形波
により周波数変調されたＦＭ変調波を用いる場合を例に
説明する。図３（ａ）は周波数変調のタイムチャートで
あり、図３（ｂ）は位相変調のタイムチャートである。
【００９１】周波数ωSで搬送波の周波数をω±ωmに変
化させる場合には、次のような式が成立する。
【００９２】
【数４】
【００９３】
【数５】
【００９４】そして、ｇ（ｔ）をフーリエ展開すると、
【００９５】
【数６】
【００９６】と表すことができる。この式から、ω0を
中心として、ωSごとに変調波のスペクトルが存在する
ことがわかる。
【００９７】νmodのスペクトル分布を一般の場合につ
いて計算するのは困難であるため、ω0＝ｎωS、ωm＝
ｍωS（ｎ≫１、ｎ≫ｍ）として、スペクトル分布を計
算すると、νmodは容易にフーリエ展開することがで
き、１／ｎを１に対して無視すると、
【００９８】
【数７】
【００９９】となる。従って、方形波による周波数変調
の場合には、このような理論振幅分布で補正を行えばよ
い。
【０１００】このように本変形例によれば、変調波を用
いて校正データを取得するので、短時間で校正データを
取得することができる。
【０１０１】（変形例（その２））次に、本実施形態に
よる相関関数測定方法及び装置の変形例（その２）を図
４を用いて説明する。図４は、本変形例による相関関数
測定装置を示すブロック図である。
【０１０２】上述したように、前変換部１８ａ、１８ｂ
の周波数特性は、周囲温度の変化や経時変化によって変
動する。しかし、変動量の主たるものは、ベクトルのオ
フセットである。従って、周囲温度が変化したり、経時
変化が生じたとしても、必ずしもすべての校正データを
更新する必要はない。そこで、本変形例では、校正デー
タメモリ２８に記憶された校正データの一部のみを更新
することにより、データ更新の迅速化を図っている。
【０１０３】即ち、図４に示すように、校正データメモ
リ２８には、回帰式メモリ３２が接続されており、回帰
式メモリ３２には校正データ演算部２６により求められ
た校正データが入力されるようになっている。
【０１０４】回帰式メモリ３２は、校正データメモリ２
８に記憶された校正データに基づいて回帰式を作成し、
回帰式を決定するのに必要な校正データのみを取得して
回帰式を更新する。そして、更新された回帰式に基づい
て、校正データメモリ２８の校正データを更新する。
【０１０５】図１に示す第１実施形態による相関関数測
定装置では、ｆS／Ｍのステップですべての校正データ
を更新するようにしていたため、データ更新に長時間を
要していたが、本変形例では、一部の校正データのみを
更新するため、データ更新の迅速化を図ることができ
る。
【０１０６】次に、本変形例による相関関数測定方法を
図４を用いて説明する。
【０１０７】まず、図１に示す第１実施形態による相関
関数測定方法と同様にして、信号発生器１２から校正信
号を発生し、校正データを取得する。具体的には、ｆ0
±ｂｗ／２の周波数帯域において、ｆS／Ｍのステップ
で校正周波数ｆCを掃引し、Ｍ個の校正データＹ（ｆC）
を取得する。取得された校正データＹ（ｆC）は、校正
データメモリ２８に記憶される。
【０１０８】次に、取得したＭ個の校正データＹ
（ｆC）を、最小二乗法により、一定の曲線に回帰させ
る。例えば、取得したＭ個の校正データから、ｎ次の多
項式
【０１０９】
【数８】
【０１１０】に回帰させる。こうして求められた回帰式
は、回帰式メモリ３２に記憶される。
【０１１１】そして、図１に示す第１実施形態による相
関関数測定方法と同様にして、電波の監視を行う。
【０１１２】校正データの更新を行う際には、回帰曲線
を決定するのに必要なポイントについてのみ校正データ
を取得し、回帰式メモリ３２に記憶された回帰式を更新
する。例えば、変極点をｎ＋１個以上測定することによ
り回帰曲線を決定しうる場合には、当該変極点のデータ
のみを取得し、回帰式を更新する。そして、更新された
回帰式から導き出される校正データを、校正データメモ
リ２８に入力することにより、校正データの更新を完了
する。
【０１１３】なお、回帰式を更新した際に標準偏差を計
算し、その標準偏差の値が一定値以上の場合に、Ｍ個の
校正データのすべてを更新するようにしてもよい。
【０１１４】このように、本変形例では、必要なポイン
トについてのみ校正データを取得すればよいので、より
迅速に校正データを更新することができる。
【０１１５】［第２実施形態］本発明の第２実施形態に
よる相関関数測定方法及び装置を図５を用いて説明す
る。図５は、本実施形態による相関関数測定装置を示す
ブロック図である。図１乃至図４に示す第１実施形態に
よる相関関数測定方法及び装置と同一の構成要素には、
同一の符号を付して説明を省略または簡潔にする。
【０１１６】図２に示す第１実施形態の変形例（その
１）による相関関数測定方法及び装置では、校正データ
の取得の迅速化を図るべく、校正信号としてＦＭ波等の
変調波を用い、理論振幅分布で補正することにより、校
正データＹ（ｆC）を取得している。
【０１１７】変調波は多様な周波数成分を含むため、Ｃ
Ｗ波を用いた場合のように校正周波数を掃引する必要は
なく、短時間で校正データＹ（ｆC）を取得することが
できる。
【０１１８】しかしながら、信号発生器から出力される
変調波は、単純な関数で変調されているわけではないた
め、理論振幅分布を求めることは非常に困難である。こ
のため、変調波を校正信号として用いた場合には、正確
な理論振幅分布で補正することができず、校正精度の劣
化を招いてしまう。
【０１１９】本発明の第２実施形態による相関関数測定
方法及び装置は、校正精度を劣化することなく、迅速に
校正を行うことに主な特徴がある。
【０１２０】（相関関数測定装置）まず、本実施形態に
よる相関関数測定装置について説明する。
【０１２１】図５に示すように、本実施形態による相関
関数測定装置は、校正信号の周波数を制御する校正信号
制御部１１０と、校正信号を発生する信号発生器１１２
と、校正信号を分配するパワースプリッタ１１４と、校
正時と測定時とで回路を切り換えるスイッチ１１６ａ、
１１６ｂと、所定の周波数成分のみを帯域制限し、中間
周波数ＩＦ（Intermediate-Frequency）に変換して出力
する前変換部１１８ａ、１１８ｂと、前変換部１１８
ａ、１１８ｂから出力されるＩＦ信号をサンプリングし
て記憶するデータメモリ１２０ａ、１２０ｂと、データ
メモリ１２０ａ、１２０ｂに記憶されたデータをフーリ
エ変換してスペクトルＳA（ｆ）、ＳB（ｆ）を求めるフ
ーリエ変換部１２２ａ、１２２ｂと、校正時と測定時と
で回路を切り換えるスイッチ１２４ａ、１２４ｂと、各
校正周波数に応じたスペクトルから校正データＹ（ｆ）
を求める演算部１２６と、各校正周波数毎の校正データ
Ｙ（ｆ）を記憶するメモリ１２８と、振幅校正データＭ
ａｇ（ｆ）と位相校正データＰｈａｓｅ（ｆ）とをそれ
ぞれ求める演算部１３０と、振幅校正データＭａｇ
（ｆ）と位相校正データＰｈａｓｅ（ｆ）とを記憶する
メモリ１３２と、振幅校正データＭａｇ（ｆ）と位相校
正データＰｈａｓｅ（ｆ）とから校正データＹ′（ｆ）
を求める演算部１３４と、校正データＹ′（ｆ）を記憶
するメモリ１３６と、校正データＹ′（ｆ）を参照しつ
つ観測信号のスペクトルの相関ベクトルを求める相関ベ
クトル演算部１３８とを有している。
【０１２２】スイッチ１１６ａの校正側ｃはパワースプ
リッタ１１４に接続されており、スイッチ１１６ａの測
定側ｍは固定アンテナに接続されている。また、スイッ
チ１１６ｂの校正側ｃはパワースプリッタ１１４に接続
されており、スイッチ１１６ｂの測定側ｍは走査アンテ
ナに接続されている。また、スイッチ１２４ａ、１２４
ｂの校正側ｃは演算部１２６に接続されており、スイッ
チ１２４ａ、１２４ｂの測定側ｍは相関ベクトル演算部
１３８に接続されている。スイッチ１１６ａ、１１６
ｂ、１２４ａ、１２４ｂを校正側ｃに設定すると、校正
データを取得する回路が構成され、スイッチ１１６ａ、
１１６ｂ、１２４ａ、１２４ｂを測定側ｍに設定する
と、相関関数の測定を行う回路が構成されるようになっ
ている。
【０１２３】本実施形態による相関関数測定装置は、変
調波を用いて位相校正データＰｈａｓｅ（ｆ）を取得
し、ＣＷ波（Continuous wave、連続波）を用いて振幅
校正データＭａｇ（ｆ）を取得し、これらのデータを用
いて所定の演算を行うことにより校正データＹ′（ｆ）
を求めることに主な特徴がある。
【０１２４】前変換部１１８ａ、１１８ｂの周波数応答
のドリフトは、位相については変動幅が比較的大きい
が、振幅については変動幅が小さいという特徴を有して
いる。このため、位相については頻繁に校正する必要が
あるが、振幅については必ずしも頻繁に校正する必要は
ない。
【０１２５】一方、変調波を用いて校正データを取得す
る場合には、振幅については正確に測定できないもの
の、位相については正確に測定できる。しかも、変調波
は多様な周波数成分を含んでいるため、極めて短時間で
校正データを取得することが可能である。
【０１２６】そこで、本実施形態では、頻繁に更新する
必要のある位相校正データＰｈａｓｅ（ｆ）を変調波を
用いて高速に取得し、振幅校正データＭａｇ（ｆ）をＣ
Ｗ波を用いて正確に取得し、これらのデータを用いて演
算することにより校正データＹ′（ｆ）を求めている。
従って、本実施形態によれば、校正精度を劣化すること
なく、全体として迅速に校正を行うことができる。
【０１２７】また、本実施形態では、ＣＷ波より成る校
正信号を所定のステップで掃引し、補間を行うことによ
り振幅校正データＭａｇ（ｆ）を求める。従って、本実
施形態によれば、更に校正の迅速化を図ることができ
る。
【０１２８】以下、本実施形態による相関関数測定装置
の各構成部分について個々に詳述する。
【０１２９】（ａ） 校正信号制御部１１０ 校正信号制御部１１０は、信号発生器１１２から出力す
る校正信号の種類や周波数を適宜制御するものである。
【０１３０】振幅校正データＭａｇ（ｆ）を取得する際
には、信号発生器１１２からＣＷ波が出力されるよう
に、信号発生器１１２を制御する。そして、信号発生器
１１２に校正周波数ｆCを入力し、信号発生器１１２か
ら出力される校正信号の周波数を掃引する。具体的に
は、ｆ0±ｂｗ／２の周波数帯域において、ｆS／Ｍのス
テップで校正周波数ｆCを掃引する。
【０１３１】ここで、ｆ0は中心周波数であり、ｂｗは
帯域幅であり、ｆSは発振器１２１からデータメモリ１
２０ａ、１２０ｂに供給されるサンプリング周波数であ
り、Ｍは任意の整数である。
【０１３２】なお、校正周波数ｆCは、信号発生器１１
２のみならず、演算部１２６にも入力されるようになっ
ている。
【０１３３】位相校正データＰｈａｓｅ（ｆ）を取得す
る際には、信号発生器１１２から変調波が出力されるよ
うに、信号発生器１１２を制御する。変調波には多様な
周波数成分が含まれているので、短時間で位相校正デー
タＰｈａｓｅ（ｆ）を取得することが可能となる。
【０１３４】（ｂ） 信号発生器１１２ 信号発生器１１２は、校正信号制御部１１０からの指示
に応じて、ＣＷ波や変調波を適宜出力する。
【０１３５】ＣＷ波を出力する際には、制御部１１０か
ら供給される校正周波数ｆCに応じた周波数の校正信号
を発生する。信号発生器１１２は、周波数に依存するこ
となく、一定の振幅の校正信号を発生する。
【０１３６】変調波を出力する際には、例えば上記の帯
域幅ｂｗより広い周波数帯域で変調されたＦＭ波等を出
力する。
【０１３７】（ｃ） パワースプリッタ１１４ パワースプリッタ１１４は、信号発生器１１２から出力
される校正信号を２つの経路に分配するものである。
【０１３８】パワースプリッタ１１４の出力は、スイッ
チ１１６ａ、１１６ｂの校正側ｃに接続されている。
【０１３９】（ｄ） スイッチ１１６ａ、１１６ｂ スイッチ１１６ａ、１１６ｂは、校正時と測定時とで入
力信号を切り換えるものである。スイッチ１１６ａ、１
１６ｂに入力される信号は高周波であるため、高周波用
のスイッチが用いられる。
【０１４０】スイッチ１１６ａの校正側ｃには、パワー
スプリッタ１１４が接続されており、スイッチ１１６ａ
の測定側ｍには、例えば固定アンテナ（図示せず）が接
続されるようになっている。
【０１４１】スイッチ１１６ａを校正側ｃに設定する
と、信号発生器１１２からの校正信号がスイッチ１１６
ａを介して前変換部１１８ａに入力されるようになり、
スイッチ１１６ａを測定側ｍに設定すると、固定アンテ
ナからの観測信号ａがスイッチ１１６ａを介して前変換
部１１８ａに入力されるようになる。
【０１４２】一方、スイッチ１１６ｂの校正側ｃには、
パワースプリッタ１１４が接続されており、スイッチ１
１６ｂの測定側ｍには、例えば走査アンテナ（図示せ
ず）が接続されるようになっている。スイッチ１１６ｂ
を校正側ｃに設定すると、信号発生器１１２からの校正
信号がスイッチ１１６ｂを介して前変換部１１８ｂに入
力されるようになり、スイッチ１１６ｂを測定側ｍに設
定すると、走査アンテナからの観測信号ｂがスイッチ１
１６ｂを介して前変換部１１８ｂに入力されるようにな
る。
【０１４３】即ち、スイッチ１１６ａ、１１６ｂを校正
側ｃに設定すると、校正データを取得する回路が構成さ
れ、スイッチ１１６ａ、１１６ｂを測定側ｍに設定する
と、相関関数の測定を行う回路が構成される。
【０１４４】（ｅ） 前変換部１１８ａ、１１８ｂ スイッチ１１６ａ、１１６ｂの出力側には、それぞれ前
変換部１１８ａ、１１８ｂが設けられている。
【０１４５】前変換部１１８ａ、１１８ｂは、所定の周
波数帯域の成分のみを通過させ、中間周波数ＩＦに変換
して出力するものである。
【０１４６】前変換部１１８ａ、１１８ｂは、例えば、
ＲＦスペクトラムアナライザを用いて構成することがで
きる。ＲＦスペクトラムアナライザを用いて前変換部１
１８ａ、１１８ｂを構成する場合には、ゼロスパンモー
ドに設定し、基準周波数ｆRE F（図示せず）で位相をロ
ックする。
【０１４７】前変換部１１８ａ、１１８ｂには、周波数
応答のドリフトが生ずるが、このドリフトは、振幅の変
動幅は小さく、位相の変動幅は比較的大きい。これは、
高周波において熱的外因が回路に影響した場合、位相、
即ち遅延時間に対して影響が大きいためと考えられる。
【０１４８】（ｆ） データメモリ１２０ａ、１２０ｂ 前変換部１１８ａ、１１８ｂの出力側には、データメモ
リ１２０ａ、１２０ｂが設けられており、データメモリ
１２０ａ、１２０ｂには、それぞれサンプリング周波数
ｆSが供給されるようになっている。
【０１４９】データメモリ１２０ａ、１２０ｂは、前変
換部１１８ａ、１１８ｂから出力されるＩＦ信号をサン
プリング周波数ｆSに基づいてサンプリングし、Ａ／Ｄ
変換した後、メモリに記憶する。データメモリ１２０
ａ、１２０ｂには、Ｍ個のデータが記憶される。なお、
サンプリング周波数ｆSは、ｂｗ＝ｆS／２となるように
設定されている。
【０１５０】（ｇ） フーリエ変換部１２２ａ、１２２
ｂ データメモリ１２０ａ、１２０ｂの出力側には、フーリ
エ変換部１２２ａ、１２２ｂが設けられている。
【０１５１】フーリエ変換部１２２ａ、１２２ｂは、デ
ータメモリ１２０ａ、１２０ｂから出力されるデータを
フーリエ変換して、スペクトルＳA（ｆ）、ＳB（ｆ）を
出力するものである。
【０１５２】（ｈ） スイッチ１２４ａ、１２４ｂ フーリエ変換部１２２ａ、１２２ｂの出力側には、スイ
ッチ１２４ａ、１２４ｂが設けられている。
【０１５３】スイッチ１２４ａ、１２４ｂは、スイッチ
１１６ａ、１１６ｂと同様、校正時と測定時とで回路を
切り換えるものである。
【０１５４】スイッチ１２４ａの校正側ｃは、演算部１
２６に接続されており、スイッチ１２４ａの測定側ｍ
は、相関ベクトル演算部１３８に接続されている。ま
た、スイッチ１２４ｂの校正側ｃは、演算部１２６に接
続されており、スイッチ１２４ｂの測定側ｍは、相関ベ
クトル演算部１３８に接続されている。
【０１５５】スイッチ１２４ａ、１２４ｂを校正側ｃに
設定すると、フーリエ変換部１２２ａ、１２２ｂから出
力されるスペクトルＳA（ｆC）、ＳB（ｆC）が演算部１
２６に入力されるようになり、スイッチ１２４ａ、１２
４ｂを測定側ｍに設定すると、フーリエ変換部１２２
ａ、１２２ｂから出力されるスペクトルＳA（ｆ）、ＳB
（ｆ）が相関ベクトル演算部１３８に入力されるように
なる。
【０１５６】（ｉ） 演算部１２６ 演算部１２６は、各周波数毎の校正データＹ（ｆC）を
演算するものである。
【０１５７】演算部１２６には、校正信号制御部１１０
から校正周波数ｆCが入力されるようになっており、ス
イッチ１２４ａ、１２４ｂの校正側ｃからスペクトルＳ
A（ｆC）、ＳB（ｆC）が入力されるようになっている。
【０１５８】演算部１２６は、以下の式から、校正周波
数ｆC毎に校正データＹ（ｆC）を求め、校正データメモ
リ１２８に出力する。なお、＊は、複素共役を示してい
る。
【０１５９】
【数９】
【０１６０】（ｊ） メモリ１２８ 演算部１２６の出力側には、メモリ１２８が設けられて
いる。
【０１６１】メモリ１２８は、演算部１２６により求め
られた各周波数毎の校正データＹ（ｆ）を記憶するもの
である。
【０１６２】振幅校正データを求める際には、ｆ0±ｂ
ｗ／２の周波数帯域において、ｆS／Ｍのステップで校
正データＹ（ｆ）が求められ、これら校正データＹ
（ｆ）がメモリ１２８に記憶される。
【０１６３】位相校正データを求める際にも、各周波数
毎の校正データＹ（ｆ）が記憶される。
【０１６４】（ｋ） 演算部１３０ メモリ１２８の出力側には、演算部１３０が設けられて
いる。
【０１６５】演算部１３０は、振幅校正データＭａｇ
（ｆ）及び位相校正データＰｈａｓｅ（ｆ）を演算する
ものである。
【０１６６】振幅校正データＭａｇ（ｆ）を求める際に
は、ＣＷ波より成る校正信号をｆS／Ｍのステップで掃
引しつつ、以下の式により、データＳ（ｆ）を求める。
なお、必ずしもｆS／Ｍのステップで校正信号を掃引す
る必要はなく、ｆS／Ｍより小さいステップで校正信号
を掃引するようにしてもよい。
【０１６７】
【数１０】
【０１６８】そして、データＳ（ｆ）について補間を行
い、これにより振幅校正データＭａｇ（ｆ）を求める。
なお、補間の際には、例えばスプライン補間を用いるこ
とができる。
【０１６９】補間を行うことにより振幅校正データＭａ
ｇ（ｆ）を求めるので、ＣＷ波を用いているにもかかわ
らず、振幅校正データＭａｇ（ｆ）を迅速に求めること
ができる。
【０１７０】位相校正データＰｈａｓｅ（ｆ）を求める
際には、以下の式により、位相校正データＰｈａｓｅ
（ｆ）を求める。
【０１７１】
【数１１】
【０１７２】（ｌ） メモリ１３２ 演算部１３０の出力側には、メモリ１３２が設けられて
いる。
【０１７３】メモリ１３２は、演算部１３０により求め
られた振幅校正用データＭａｇ（ｆ）と位相校正用デー
タＰｈａｓｅ（ｆ）とを、それぞれ記憶しておくもので
ある。
【０１７４】（ｍ） 演算部１３４ メモリ１３２の出力側には、演算部１３４が設けられて
いる。
【０１７５】演算部１３４は、振幅校正用データＭａｇ
（ｆ）と位相校正用データＰｈａｓｅ（ｆ）を用いて、
校正データＹ′（ｆ）を演算するものである。
【０１７６】校正データＹ′（ｆ）は、以下の式により
求められる。
【０１７７】
【数１２】
【０１７８】（ｎ） メモリ１３６ 演算部１３４の出力側には、メモリ１３６が設けられて
いる。
【０１７９】メモリ１３６は、演算部１３４で求められ
た校正データＹ′（ｆ）を記憶しておくものである。
【０１８０】（ｏ） 相関ベクトル演算部１３８ 相関ベクトル演算部１３８は、相関ベクトルを演算する
ものである。
【０１８１】相関ベクトル演算部１３８には、フーリエ
変換部１２２ａ、１２２ｂから出力されたスペクトルＳ
A（ｆ）、ＳB（ｆ）が、スイッチ１２４ａ、１２４ｂを
介して入力されるようになっている。
【０１８２】相関ベクトル演算部１３８は、観測信号の
周波数に応じて校正データメモリ１３６から校正データ
Ｙ′（ｆ）を適宜参照し、以下の式に基づいて相関ベク
トルを演算する。
【０１８３】
【数１３】
【０１８４】相関ベクトル演算部１３８により求められ
た相関ベクトルは、波形像再生制御部（図示せず）等に
出力され、表示部（図示せず）に表示される。
【０１８５】なお、本実施形態による相関関数測定装置
は、本願出願人による特願平７−２８９８４８号明細書
に記載された波源像可視化装置等に適用可能である。波
源像可視化装置に適用する場合の構成等については、当
該明細書を参照されたい。
【０１８６】（データの精度）次に、本実施形態で取得
される振幅校正データＭａｇ（ｆ）の精度について図６
を用いて説明する。
【０１８７】図６は、ＣＷ波を用いて取得した位相校正
データと振幅校正データとを示すグラフであり、いずれ
も７．５ＭＨｚ〜１２．５ＭＨｚの範囲で測定したもの
である。
【０１８８】図６（ａ）の実施例１は、スプライン補間
を用いて取得した位相校正データを示しており、実施例
２は、スプライン補間を用いて取得した振幅校正データ
を示している。測定点は９点であり、測定の所要時間は
約５秒であった。なお、図中の丸は測定点を示してい
る。
【０１８９】図６（ｂ）の比較例１は、補間を行うこと
なく取得した位相校正データを示しており、比較例２
は、補間を行うことなく取得した振幅校正データを示し
ている。測定点は１２９点であり、測定の所要時間は約
１分であった。
【０１９０】実施例２と比較例２とを比較して分かるよ
うに、両者の振幅の差は約０．１ｄＢ以下に抑えられて
いる。このことから分かるように、スプライン補間を用
いる場合であっても、振幅校正データＭａｇ（ｆ）を高
精度に求めることが可能である。
【０１９１】次に、本実施形態で取得される位相校正デ
ータＰｈａｓｅ（ｆ）の精度について図７を用いて説明
する。
【０１９２】図７（ａ）は、変調波を用いて取得した位
相校正データと振幅校正データとを示すグラフである。
図７（ａ）の実施例３は、変調波を用いて取得された位
相校正データを示しており、実施例４は、変調波を用い
て取得された振幅校正データを示している。
【０１９３】図７（ｂ）は、図７（ａ）に示すデータを
取得する際に用いた変調波のスペクトル示すグラフであ
り、中心周波数２ＧＨｚ、スパン５０ＭＨｚで表したも
のである。なお、図７（ｂ）に示す変調波の２ＧＨｚ±
２．５ＭＨｚを校正信号として用いて、図７（ａ）に示
すデータが得られている。
【０１９４】図７（ａ）に示す実施例３と、図６（ｂ）
に示す比較例１と比較して分かるように、両者の位相差
は約０．２度以下に抑えられている。このことから分か
るように、変調波を用いた場合であっても、位相校正デ
ータＰｈａｓｅ（ｆ）を高精度に求めることが可能であ
る。
【０１９５】（相関関数測定方法）次に、本実施形態に
よる相関関数測定方法を図５を用いて説明する。
【０１９６】まず、相関関数の測定に先立って、校正デ
ータＹ′（ｆ）を取得する。校正データＹ′（ｆ）を取
得する際には、各スイッチ１１６ａ、１１６ｂ、１２４
ａ、１２４ｂをそれぞれ校正側ｃに設定する。これによ
り、信号発生器１１２から出力された校正信号が前変換
部１１８ａ、１１８ｂに入力されるようになり、また、
フーリエ変換されたスペクトルＳA（ｆC）、ＳB（ｆC）
が演算部１２６に入力されるようになる。
【０１９７】次に、振幅校正データＭａｇ（ｆ）を取得
すべく、ＣＷ波を出力するように信号発生器１１２を制
御する。
【０１９８】そして、校正周波数ｆCを掃引して、各校
正周波数ｆC毎に校正データＹ（ｆC）を求める。具体的
には、校正信号制御部１１０から信号発生器１１２に供
給する校正周波数ｆCを、ｆS／Ｍのステップ、ｆ0±ｂ
ｗ／２の周波数帯域で順次掃引していき、演算部１２６
において、各校正周波数ｆC毎に校正データＹ（ｆC）を
求めていく。こうして求められた校正データＹ（ｆC）
は、メモリ１２８に順次記憶される。
【０１９９】更に、演算部１３０において、上述したよ
うな演算を行ない、これにより、振幅校正データＭａｇ
（ｆ）を求める。こうして求められた振幅校正データＭ
ａｇ（ｆ）は、メモリ１３２に記憶される。例えば、帯
域幅ｂｗが１０ＭＨｚ、測定点が９点の場合には、振幅
校正データＭａｇ（ｆ）は、約６秒で取得することが可
能である。
【０２００】次に、位相校正データＰｈａｓｅ（ｆ）を
取得すべく、変調波を出力するように信号発生器１１２
を制御する。
【０２０１】そして、演算部１２６において、各校正周
波数ｆC毎の校正データＹ（ｆC）を求める。例えば、帯
域幅ｂｗが１０ＭＨｚの場合には、位相校正データＰｈ
ａｓｅ（ｆ）は、約０．５秒で取得することが可能であ
る。
【０２０２】更に、演算部１３０において、上述したよ
うな演算を行い、これにより、位相校正データＰｈａｓ
ｅ（ｆ）を求める。こうして求められた位相校正データ
Ｐｈａｓｅ（ｆ）は、メモリ１３２に記憶される。
【０２０３】次に、演算部１３４において、振幅校正デ
ータＭａｇ（ｆ）と位相校正データＰｈａｓｅ（ｆ）と
を用いて上述したような演算を行うことにより、校正デ
ータＹ′（ｆ）を求める。こうして求められたＹ′
（ｆ）は、メモリ１３６に記憶される。
【０２０４】こうして、校正データＹ′（ｆ）の取得が
完了する。
【０２０５】相関関数の測定を行う際には、各スイッチ
１１６ａ、１１６ｂ、１２４ａ、１２４ｂを測定側ｍに
設定する。これにより、固定アンテナから観測信号ａが
入力されるようになり、走査アンテナから観測信号ｂが
入力されるようになる。また、フーリエ変換部１２２
ａ、１２２ｂから出力されるスペクトルＳA（ｆ）、ＳB
（ｆ）が、相関ベクトル演算部１３８に入力されるよう
になる。
【０２０６】相関ベクトルを演算する際には、メモリ１
３６に記憶された校正データＹ′（ｆ）を、観測信号の
周波数に応じて適宜参照し、相関ベクトルを演算する。
こうして、観測信号ａ、ｂの周波数成分に応じて各周波
数毎に校正が行われ、高い精度で相関ベクトルが求めら
れる。なお、相関関数測定方法の詳細は、本願出願人に
よる特願平７−２８９８４８号明細書を参照されたい。
【０２０７】周囲温度が変化した場合、経時変化が生じ
た場合など、校正データＹ′（ｆ）の信憑性が低くなっ
た場合には、上記と同様にして位相校正データＰｈａｓ
ｅ（ｆ）のみを取得し、所定の演算を行うことにより校
正データＹ′（ｆ）を更新する。
【０２０８】そして、新たな校正データＹ′（ｆ）を参
照しつつ、相関ベクトルを演算すればよい。
【０２０９】なお、振幅校正データＭａｇ（ｆ）の信憑
性も低くなったと考えられる場合には、位相校正データ
Ｐｈａｓｅ（ｆ）のみならず、振幅校正データＭａｇ
（ｆ）をも取得し、これらのデータを用いて演算するこ
とにより校正データＹ′（ｆ）を求めればよい。例え
ば、振幅校正データＭａｇ（ｆ）を取得した際における
位相校正用データＰｈａｓｅ（ｆ）を記憶しておき、こ
の位相校正データＰｈａｓｅ（ｆ）に対して所定値以上
位相校正データＰｈａｓｅ（ｆ）の値が変化した場合に
は、振幅校正用データＭａｇ（ｆ）の信憑性も低くなっ
ている蓋然性が高いため、振幅校正用データＭａｇ
（ｆ）を更新するようにしてもよい。位相校正用データ
Ｐｈａｓｅ（ｆ）の値が所定以上変化したということ
は、例えば標準偏差等により求めることが可能である。
【０２１０】このように本実施形態によれば、頻繁に取
得する必要のある位相校正データＰｈａｓｅ（ｆ）を変
調波を用いて取得し、頻繁に取得する必要のない振幅校
正データＭａｇ（ｆ）はＣＷ波を用いて取得し、これら
のデータを用いて各周波数に応じた校正データＹ′
（ｆ）を求めるので、校正精度を劣化することなく、迅
速に校正を行うことができる。
【０２１１】また、本実施形態では、ＣＷ波より成る校
正信号を所定のステップで掃引し、補間を行うことによ
り振幅校正データＭａｇ（ｆ）を求めるので、更に校正
を迅速に行うことができる。
【０２１２】［変形実施形態］本発明は上記実施形態に
限らず種々の変形が可能である。
【０２１３】例えば、第１及び第２実施形態では、電波
監視を行う際に固定アンテナ及び走査アンテナから入力
される観測信号の相関関数を測定する場合を例に説明し
たが、電波監視を行う場合のみならず、複数の観測信号
の相関関数を測定する際に広く適用することができる。
【０２１４】また、第１及び第２実施形態では、２つの
観測信号の相関ベクトルを演算する場合に適用したが、
３つ以上の相関ベクトルを演算する場合にも応用するこ
とができる。
【０２１５】また、第２実施形態では、変調波を用いて
位相校正データを取得したが、必ずしも変調波を用いて
位相校正データを取得しなくてもよい。例えば、迅速な
校正が要求されない場合には、連続波を用いて位相校正
データを取得してもよい。
【０２１６】
【発明の効果】以上の通り、本発明によれば、相関関数
の測定に先立って、各周波数毎に校正データを取得して
おき、周波数に応じた校正データを用いて相関ベクトル
を求めるので、相関ベクトルを高い精度で求めることが
できる。また、校正データの信憑性が低くなった場合に
は、適宜校正データの更新を行うので、相関ベクトルの
校正精度を維持することができる。従って、本発明によ
れば、高精度に相関関数の測定を行うことができる相関
関数測定方法及び装置を提供することができる。
【０２１７】また、本発明によれば、頻繁に取得する必
要のある位相校正データＰｈａｓｅ（ｆ）を変調波を用
いて取得し、頻繁に取得する必要のない振幅校正データ
Ｍａｇ（ｆ）はＣＷ波を用いて取得し、これらのデータ
を用いて各周波数に応じた校正データＹ′（ｆ）を求め
るので、校正精度を劣化することなく、迅速に校正を行
うことができる。
【０２１８】また、本発明によれば、ＣＷ波より成る校
正信号を所定のステップで掃引し、補間を行うことによ
り振幅校正データＭａｇ（ｆ）を求めるので、更に校正
を迅速に行うことができる。
【図面の簡単な説明】
【図１】本発明の第１実施形態による相関関数測定装置
を示すブロック図である。
【図２】本発明の第１実施形態の変形例（その１）によ
る相関関数測定装置を示すブロック図である。
【図３】周波数変調と位相変調とを示すタイムチャート
である。
【図４】本発明の第１実施形態の変形例（その２）によ
る相関関数測定装置を示すブロック図である。
【図５】本発明の第２実施形態による相関関数測定装置
を示すブロック図である。
【図６】ＣＷ波を用いて取得した位相校正データと振幅
校正データとを示すグラフである。
【図７】変調波を用いて取得した位相校正データと振幅
校正データとを示すグラフである。
【図８】提案されている波源像可視化装置を示すブロッ
ク図である。
【符号の説明】
１０…校正周波数制御部 １２…信号発生器 １４…パワースプリッタ １６ａ、１６ｂ…スイッチ １８ａ、１８ｂ…前変換部 ２０ａ、２０ｂ…データメモリ ２１…発振器 ２２ａ、２２ｂ…フーリエ変換部 ２４ａ、２４ｂ…スイッチ ２６…校正データ演算部 ２８…校正データメモリ ３０…相関ベクトル演算部 ３２…回帰式メモリ １１０…校正信号制御部 １１２…信号発生器 １１４…パワースプリッタ １１６ａ、１１６ｂ…スイッチ １１８ａ、１１８ｂ…前変換部 １２０ａ、１２０ｂ…データメモリ １２１…発振器 １２２ａ、１２２ｂ…フーリエ変換部 １２４ａ、１２４ｂ…スイッチ １２６…演算部 １２８…メモリ １３０…演算部 １３２…メモリ １３４…演算部 １３６…メモリ １３８…相関ベクトル演算部 ２１２…信号発生器 ２１４…パワースプリッタ ２１６ａ、２１６ｂ…スイッチ ２１８ａ、２１８ｂ…前変換部 ２２０ａ、２２０ｂ…データメモリ ２２２ａ、２２２ｂ…フーリエ変換部 ２２４…スイッチ ２２８…校正データメモリ ２２９…校正部 ２３０…相関ベクトル演算部
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） Ｈ０４Ｂ 17/00 Ｈ０４Ｂ 17/00 Ｃ

Claims (18)

    【特許請求の範囲】
  1. 【請求項１】 第１の観測信号を第１の信号処理手段に
    より処理して第１のスペクトルを求め、第２の観測信号
    を第２の信号処理手段により処理して第２のスペクトル
    を求め、前記第１のスペクトルと前記第２のスペクトル
    との相関関数を測定する相関関数測定方法であって、 相関関数の測定に先立って、前記第１の信号処理手段及
    び前記第２の信号処理手段に校正信号を入力し、前記第
    １の信号処理手段により求められたスペクトルと、前記
    第２の信号処理手段により求められたスペクトルとか
    ら、前記校正信号の周波数に応じた校正値を求め、 相関関数の測定の際には、周波数に応じた前記校正値を
    用いつつ、前記第１のスペクトルと前記第２のスペクト
    ルとの相関関数を求めることを特徴とする相関関数測定
    方法。
  2. 【請求項２】 請求項１記載の相関関数測定方法におい
    て、 前記校正値を求める際には、前記校正信号の周波数を掃
    引することにより、各周波数に応じた前記校正値を求め
    ることを特徴とする相関関数測定方法。
  3. 【請求項３】 請求項１又は２記載の相関関数測定方法
    において、 前記校正値から回帰式を求め、 前記校正値を更新する際には、少なくとも前記回帰式を
    決定するのに必要な校正値を再度求めることにより前記
    回帰式を更新し、更新された前記回帰式に基づいて前記
    校正値を更新することを特徴とする相関関数測定方法。
  4. 【請求項４】 請求項１記載の相関関数測定方法におい
    て、 前記校正値を求める際には、前記第１の信号処理手段及
    び前記第２の信号処理手段に変調波を入力することによ
    り、各周波数に応じた校正値を求めることを特徴とする
    相関関数測定方法。
  5. 【請求項５】 請求項１記載の相関関数測定方法におい
    て、 前記校正値を求める際には、前記第１の信号処理手段及
    び前記第２の信号処理手段に、第１の校正信号を周波数
    を掃引しつつ入力し、前記第１の信号処理手段により求
    められたスペクトルと、前記第２の信号処理手段により
    求められたスペクトルとから、各周波数に応じた振幅校
    正データを求め、前記第１の信号処理手段及び前記第２
    の信号処理手段に第２の校正信号を入力し、前記第１の
    信号処理手段により求められたスペクトルと、前記第２
    の信号処理手段により求められたスペクトルとから、各
    周波数に応じた位相校正データを求め、前記振幅校正デ
    ータと前記位相校正データとから各周波数に応じた前記
    校正値を求めることを特徴とする相関関数測定方法。
  6. 【請求項６】 請求項５記載の相関関数測定方法におい
    て、 前記第２の校正信号は、変調波より成ることを特徴とす
    る相関関数測定方法。
  7. 【請求項７】 請求項５又は６記載の相関関数測定方法
    において、 前記振幅校正データを求める際には、前記第１の校正信
    号の周波数を所定のステップで掃引し、補間を行うこと
    により前記振幅校正データを求めることを特徴とする相
    関関数測定方法。
  8. 【請求項８】 請求項５乃至７のいずれか１項記載の相
    関関数測定方法において、 前記位相校正データを、前記振幅校正データより高い頻
    度で更新することを特徴とする相関関数測定方法。
  9. 【請求項９】 請求項８記載の相関関数測定方法におい
    て、 前記振幅校正データを更新した際の前記位相校正データ
    に対して、更新された前記位相校正データが所定値以上
    変化している場合に、前記振幅校正データを更新するこ
    とを特徴とする相関関数測定方法。
  10. 【請求項１０】 第１の観測信号を処理して第１のスペ
    クトルを求める第１の信号処理手段と、 第２の観測信号を処理して第２のスペクトルを求める第
    ２の信号処理手段と、 前記第１の信号処理手段及び前記第２の信号処理手段に
    校正信号を入力し、前記第１の信号処理手段により求め
    られたスペクトルと、前記第２の信号処理手段により求
    められたスペクトルとから、各周波数に応じた校正値を
    求める校正値演算手段と、 前記校正値演算手段により求められた前記校正値を用い
    つつ、前記第１のスペクトルと前記第２のスペクトルと
    の相関関数を求める相関関数演算手段とを有することを
    特徴とする相関関数測定装置。
  11. 【請求項１１】 請求項１０記載の相関関数測定装置に
    おいて、 前記校正値演算手段は、前記校正信号の周波数を掃引す
    ることにより、各周波数に応じた前記校正値を求めるこ
    とを特徴とする相関関数測定装置。
  12. 【請求項１２】 請求項１０又は１１記載の相関関数測
    定装置において、 前記校正値演算手段は、前記校正値から回帰式を求め、
    少なくとも前記回帰式を決定するのに必要な校正値を再
    度求めることにより前記回帰式を更新し、更新された前
    記回帰式に基づいて前記校正値を更新することを特徴と
    する相関関数測定装置。
  13. 【請求項１３】 請求項１０記載の相関関数測定装置に
    おいて、 前記校正値演算手段は、前記第１の信号処理手段及び前
    記第２の信号処理手段に変調波を入力することにより、
    各周波数に応じた前記校正値を求めることを特徴とする
    相関関数測定装置。
  14. 【請求項１４】 請求項１０記載の相関関数測定装置に
    おいて、 前記校正値演算手段は、前記第１の信号処理手段及び前
    記第２の信号処理手段に、第１の校正信号を周波数を掃
    引しつつ入力し、前記第１の信号処理手段により求めら
    れたスペクトルと、前記第２の信号処理手段により求め
    られたスペクトルとから、各周波数に応じた振幅校正デ
    ータを求め、前記第１の信号処理手段及び前記第２の信
    号処理手段に第２の校正信号を入力し、前記第１の信号
    処理手段により求められたスペクトルと、前記第２の信
    号処理手段により求められたスペクトルとから、各周波
    数に応じた位相校正データを求め、前記振幅校正データ
    と前記位相校正データとから各周波数に応じた前記校正
    値を求めることを特徴とする相関関数測定装置。
  15. 【請求項１５】 請求項１４記載の相関関数測定装置に
    おいて、 前記第２の校正信号は、変調波より成ることを特徴とす
    る相関関数測定装置。
  16. 【請求項１６】 請求項１４又は１５記載の相関関数測
    定装置において、 校正値演算手段は、前記第１の校正信号の周波数を所定
    のステップで掃引し、補間を行うことにより前記振幅校
    正データを求めることを特徴とする相関関数測定装置。
  17. 【請求項１７】 請求項１４乃至１６のいずれか１項記
    載の相関関数測定装置において、 前記校正値演算手段は、前記位相校正データを、前記振
    幅校正データより高い頻度で更新することを特徴とする
    相関関数測定装置。
  18. 【請求項１８】 請求項１７記載の相関関数測定装置に
    おいて、 前記校正値演算手段は、前記振幅校正データを更新した
    際の前記位相校正データに対して、更新された前記位相
    校正データが所定値以上変化している場合に、前記振幅
    校正データを更新することを特徴とする相関関数測定装
    置。
JP2000296812A 2000-02-03 2000-09-28 相関関数測定方法及び装置 Expired - Fee Related JP4038009B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2000296812A JP4038009B2 (ja) 2000-02-03 2000-09-28 相関関数測定方法及び装置
US09/771,527 US6594605B2 (en) 2000-02-03 2001-01-29 Correlation function measuring method and apparatus
DE10165055A DE10165055B4 (de) 2000-02-03 2001-01-29 Korrelationsfunktions-Messverfahren und -Vorrichtung
DE10103822A DE10103822B4 (de) 2000-02-03 2001-01-29 Korrelationsfunktions-Messverfahren und -Vorrichtung
GB0102749A GB2366030B (en) 2000-02-03 2001-02-05 Correlation function measuring method and apparatus

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
JP2000026863 2000-02-03
JP2000127080 2000-04-27
JP2000-127080 2000-04-27
JP2000-26863 2000-04-27
JP2000296812A JP4038009B2 (ja) 2000-02-03 2000-09-28 相関関数測定方法及び装置

Publications (2)

Publication Number Publication Date
JP2002014600A true JP2002014600A (ja) 2002-01-18
JP4038009B2 JP4038009B2 (ja) 2008-01-23

Family

ID=27342239

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000296812A Expired - Fee Related JP4038009B2 (ja) 2000-02-03 2000-09-28 相関関数測定方法及び装置

Country Status (4)

Country Link
US (1) US6594605B2 (ja)
JP (1) JP4038009B2 (ja)
DE (2) DE10103822B4 (ja)
GB (1) GB2366030B (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008128667A (ja) * 2006-11-16 2008-06-05 Toshiba Corp 電波発射源可視化方法および装置

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6979993B2 (en) * 2001-03-14 2005-12-27 Advantest Corporation Frequency analyzing method, frequency analyzing apparatus, and spectrum analyzer
GB0113627D0 (en) 2001-06-05 2001-07-25 Univ Stirling Controller and method of controlling an apparatus
US8214682B2 (en) * 2008-04-11 2012-07-03 Oracle America, Inc. Synchronizing signals related to the operation of a computer system
US10271233B2 (en) 2013-03-15 2019-04-23 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection with temporal feature extraction within a spectrum
US10231206B2 (en) 2013-03-15 2019-03-12 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying signal-emitting devices
US10237770B2 (en) 2013-03-15 2019-03-19 DGS Global Systems, Inc. Systems, methods, and devices having databases and automated reports for electronic spectrum management
US8750156B1 (en) 2013-03-15 2014-06-10 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying open space
US10122479B2 (en) 2017-01-23 2018-11-06 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection with temporal feature extraction within a spectrum
US10257729B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems, methods, and devices having databases for electronic spectrum management
US10219163B2 (en) 2013-03-15 2019-02-26 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US9622041B2 (en) 2013-03-15 2017-04-11 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10244504B2 (en) 2013-03-15 2019-03-26 DGS Global Systems, Inc. Systems, methods, and devices for geolocation with deployable large scale arrays
US10257728B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US10299149B2 (en) 2013-03-15 2019-05-21 DGS Global Systems, Inc. Systems, methods, and devices for electronic spectrum management
US11646918B2 (en) 2013-03-15 2023-05-09 Digital Global Systems, Inc. Systems, methods, and devices for electronic spectrum management for identifying open space
US10257727B2 (en) 2013-03-15 2019-04-09 DGS Global Systems, Inc. Systems methods, and devices having databases and automated reports for electronic spectrum management
US10101049B2 (en) 2015-11-12 2018-10-16 Oracle International Corporation Determining parameters of air-cooling mechanisms
US10700794B2 (en) 2017-01-23 2020-06-30 Digital Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time within an electromagnetic spectrum
US10459020B2 (en) 2017-01-23 2019-10-29 DGS Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time within a spectrum
US10498951B2 (en) 2017-01-23 2019-12-03 Digital Global Systems, Inc. Systems, methods, and devices for unmanned vehicle detection
US10529241B2 (en) 2017-01-23 2020-01-07 Digital Global Systems, Inc. Unmanned vehicle recognition and threat management
US10943461B2 (en) 2018-08-24 2021-03-09 Digital Global Systems, Inc. Systems, methods, and devices for automatic signal detection based on power distribution by frequency over time

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4169245A (en) * 1972-07-26 1979-09-25 E-Systems, Inc. Spectral correlation
JPS59157575A (ja) * 1983-02-27 1984-09-06 Anritsu Corp スペクトラムアナライザ
US4947176A (en) * 1988-06-10 1990-08-07 Mitsubishi Denki Kabushiki Kaisha Multiple-beam antenna system
US5146471A (en) * 1989-03-23 1992-09-08 Echelon Systems Corporation Correlator for spread spectrum communications systems
JP3642844B2 (ja) 1995-11-08 2005-04-27 株式会社アドバンテスト 波源像可視化方法及び装置
US5748314A (en) * 1995-11-08 1998-05-05 Advantest Corporation Correlation function measurement method and apparatus, and wave source image visualization method and apparatus based on correlation function measurement
US6917644B2 (en) * 1996-04-25 2005-07-12 Sirf Technology, Inc. Spread spectrum receiver with multi-path correction
US5901183A (en) * 1996-09-25 1999-05-04 Magellan Corporation Signal correlation technique for a receiver of a spread spectrum signal including a pseudo-random noise code that reduces errors when a multipath signal is present
US6295442B1 (en) * 1998-12-07 2001-09-25 Ericsson Inc. Amplitude modulation to phase modulation cancellation method in an RF amplifier

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008128667A (ja) * 2006-11-16 2008-06-05 Toshiba Corp 電波発射源可視化方法および装置

Also Published As

Publication number Publication date
JP4038009B2 (ja) 2008-01-23
GB2366030B (en) 2005-01-12
US20010020220A1 (en) 2001-09-06
DE10103822B4 (de) 2009-06-04
DE10103822A1 (de) 2001-11-29
GB0102749D0 (en) 2001-03-21
GB2366030A (en) 2002-02-27
US6594605B2 (en) 2003-07-15
DE10165055B4 (de) 2009-10-15

Similar Documents

Publication Publication Date Title
JP2002014600A (ja) 相関関数測定方法及び装置
JP4618445B2 (ja) 確率密度関数で重み付けした積分処理を用いた線形性評価方法と、それを用いた回路シミュレータ、評価装置、通信回路、およびプログラム
US8364098B2 (en) PLL calibration
JP2001268150A (ja) リニアライザ
JP2007135214A (ja) フィルタ等化方法
JP3852938B2 (ja) 広帯域変調ｐｌｌおよびその変調度調整方法
JP2019191344A (ja) 位相特性校正装置及び位相特性校正方法
US6041083A (en) Method and system for tuning resonance modules
JP2022524455A (ja) 側波帯抑圧方法、装置、コンピュータ機器及びコンピュータプログラム
US9270311B2 (en) Methods and systems for calibrating an analog filter
JP2020193848A (ja) 位相特性校正装置および位相特性校正方法
JP4923276B2 (ja) 周波数変換装置の周波数応答の特徴付け方法
US6393258B1 (en) Process and system for adjusting the level of the spurious lines of the output frequency spectrum of a single-sideband frequency transposition device, in particular incorporated into a mobile telephone
US6995571B1 (en) Vector network analyzer mixer calibration using the unknown thru calibration
US5646541A (en) Apparatus and method for measuring circuit network
CN114285466A (zh) 基于双光频梳及受激布里渊散射的微波频率测量方法
JP3561184B2 (ja) Ｉｑスプリッタ装置
JP3719675B2 (ja) Ｆｍｃｗレーダ装置およびｆｍｃｗレーダ用信号処理方法
JP4376751B2 (ja) トラッキング発生器及び試験信号発生方法
CN114499705B (zh) 频响平坦度校准方法、装置、电子设备及存储介质
CN114024627B (zh) 针对宽带调制信号实现快速全频段频响补偿处理的方法、装置、处理器及其存储介质
CN115792770B (zh) 矢量网络分析仪通道间固有相参校准数据获取方法及系统
JP6568608B2 (ja) 信号測定装置及び信号測定方法
JP2003315395A (ja) ベクトル・ネットワーク・アナライザおよびその位相測定方法
JP2003110640A (ja) 直交検波回路

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040318

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070619

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070808

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20071030

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20071102

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101109

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111109

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121109

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121109

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131109

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees