JPH06180711A - 離散フーリエ変換値の遅延時間補正装置 - Google Patents

Abstract

(57)【要約】 【目的】 本発明は、ディジタル信号を離散フーリエ変
換したときの遅延時間補正装置に関し、１サンプリング
周期以下の遅延Δｔに対するフーリエ変換値を補正す
る。 【構成】 サンプリング周期をτとし、サンプリング点
をｔ＝−Δｔ，τ−Δｔ，・・・，（Ｎ−１）τ−Δｔ
とした合計Ｎ個のサンプル値ｇ（−Δｔ），ｇ（τ−Δ
ｔ），・・・，ｇ（（Ｎ−１）τ−Δｔ）に対して離散
フーリエ変換して、周波数成分Ｇ’（ｆ）＝Ｇ（ｆ）・
ｅｘｐ（−ｊ２πｆΔｔ）を得るフーリエ変換部と、位
相補正データγ＝ｅｘｐ（ｊ２πｆΔｔ）を算出する位
相補正算出部と、前記Ｇ’（ｆ）と、前記位相補正デー
タγとを乗算する位相補正乗算部とを備えて、Ｇ（ｆ）
＝Ｇ’（ｆ）・ｅｘｐ（ｊ２πｆΔｔ）により、基準サ
ンプルよりΔｔだけ遅延しているサンプルでの周波数成
分Ｇ’（ｆ）を基に、遅延Δｔに対する補正をして、基
準点サンプルでの周波数成分Ｇ（ｆ）を求めるように構
成する。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、ディジタル信号を離散
フーリエ変換したときの遅延時間補正装置に関する。

【０００２】近年、ディジタル計算機の高速化やマイク
ロプロセッサの普及が進み、音声、画像、電波等の解析
処理の分野において、ディジタル信号解析の主流である
フーリエ変換は、その重要性をますます高めている。

【０００３】例えば、電波天文学において、宇宙からの
電波を、専用のパラボラアンテナで受信した時系列のア
ナログ信号を、所定のサンプリング数Ｎでサンプリング
し、それぞれでのサンプリング位置での振幅をディジタ
ル変換した後、離散的フーリエ変換による周波数分析を
行うことにより、該受信した電波の周波数成分を認識す
ることができ、該電波の発信元である恒星等を構成して
いる元素を解析することができる。

【０００４】このとき、発信元の位置解析の分解能を高
める為、長大な距離だけ離れた複数個、例えば、２ケ所
のパラボラアンテナで、同じ発信元である恒星等の電波
を受信するようにして、等価的に直径の大きなパラボラ
アンテナで受信したのと等価な発信元位置の解析精度、
周波数成分の得られる、所謂、電波干渉計が知られてい
る。

【０００５】この電波干渉計等においては、電波の発信
元からパラボラアンテナ迄の距離が異なるため、受信し
た電波にタイムラグ（位相のずれ）が生じ、受信した２
つの生の電波情報を用いたとき、そのままの形では、相
関を得ることができないことになる。

【０００６】計算機等で、この受信信号をディジタル処
理する場合、上記のタイムラグが、サンプリング周期の
整数倍である場合には、そのずれているサンプリング周
期分だけ、クロックをずらす等して、サンプリングデー
タを補正することで、２つの受信電波の位相を合わせる
ことができるが、１サンプリング周期以内のずれである
Δｔ分のずれについての補正は、上記のようなクロック
をずらせる等の補正だけでは困難となる。

【０００７】

【従来の技術】図１３は、２つのパラボラアンテナで受
信した電波の信号処理を説明する図であり、図１３(a)
は、２つのパラボラアンテナ PA1,PA2で同じ発信源から
の電波を受信した場合の位相ずれを説明している図であ
り、図１３(b) は、受信した電波信号のサンプリング点
のずれを説明している。

【０００８】図１３(a) に示されているように、同じ発
信源からの電波を所定の基幹距離だけ離れて設けられて
いる２つのパラボラアンテナ PA1,PA2で受信した場合、
その基幹距離に比例した位相ずれ｛図１３(a) では、遅
延時間Ｄで示す｝を生じる。

【０００９】図１３(b) は、同一電波源より発せられた
電波を、上記２つのパラボラアンテナ PA1,PA2で受信し
た場合の、同一波面の電波信号に時刻信号をとったとき
の番号を示したもので、例えば、図示されているよう
に、２つのパラボラアンテナ PA1,PA2で受信した電波信
号の位相は、2.3 サンプリング分ずれている。

【００１０】この内、２サンプリング周期分のずれは、
サンプリングクロックで位相を合わせることにより、補
正可能であるが、０．3 サンプリングクロック分の位相
のずれを合わせることは、従来のクロックによる補正で
は、極めて困難なものとなる。

【００１１】従来、時系列の２つのサンプル・データ列
の相関を求める場合、相関結果として、できる限り高い
ピーク値（高い相関値）を得るためには、被相関データ
の基準データに対する上記遅延時間をできる限り少な
く、即ち、時間軸で遅延補正してから、相関を求める方
法を採っていた。

【００１２】又、相関処理にフーリエ変換を利用した一
例として、上記時間領域で相関を求める代わりに、２つ
のサンプル・データ列をフーリエ変換した周波数領域で
相関を求める方法が知られている。

【００１３】

【発明が解決しようとする課題】通常、ディジタル計算
機でフーリエ変換する場合、離散的フーリエ変換の方法
が用いられる。特に、演算速度を高速化した高速フーリ
エ変換（ＦＦＴ）法が用いられる。

【００１４】被相関データ（例えば、上記パラボラアン
テナ PA2で受信したデータ) の基準データ (例えば、上
記パラボラアンテナ PA1で受信したデータ) に対する遅
延時間をＤとした場合、次の関係があるとする。 Ｄ＝t0＋Δｔ 但し、 t0＝ｎ・τ（ｎは整数） ０≦Δｔ≦τ τ：サンプリング周期 この場合、t0はサンプリング周期τの整数倍であるの
で、この部分の遅延時間は補正可能であるが、前述のよ
うに、サンプリング周期τ以下のΔｔの補正はできな
い。

【００１５】このようなとき、上記時間領域で相関をと
る方法では、相関結果として高いピーク値を得ることが
できなく、又、フーリエ変換を用いた場合でも、被相関
データをフーリエ変換した周波数領域での位相は、一致
しなくなり、相関結果として高いピーク値（高い相関
値）を得ることができないという問題があった。

【００１６】本発明は上記従来の欠点に鑑み、ディジタ
ル信号を離散フーリエ変換したときの遅延時間補正装置
において、１サンプリング周期以下の遅延Δｔに対する
フーリエ変換値を補正する離散フーリエ変換値の遅延時
間補正装置を提供することを目的とするものである。

【００１７】

【課題を解決するための手段】図１〜図３は、本発明の
原理構成図であり、それぞれ、第１〜第３の原理図を示
している。上記の問題点は下記の如くに構成した離散フ
ーリエ変換値の遅延時間補正装置によって解決される。

【００１８】(1) アナログ信号に対するサンプリング周
期をτとし、サンプリング点をｔ＝−Δｔ，τ−Δｔ，
・・・，（Ｎ−１）τ−Δｔとした合計Ｎ個のサンプル
値ｇ（−Δｔ），ｇ（τ−Δｔ），・・・，ｇ（（Ｎ−
１）τ−Δｔ）に対して離散フーリエ変換して、周波数
成分Ｇ’（ｆ）＝Ｇ（ｆ）・ｅｘｐ（−ｊ２πｆΔｔ）
を得るフーリエ変換部 1と、位相補正データγ＝ｅｘｐ
（ｊ２πｆΔｔ）を算出する位相補正算出部 2と、前記
Ｇ’（ｆ）と、前記位相補正データγとを乗算する位相
補正乗算部 3とを備えるように構成する。｛図１参照｝ (2) 上記位相補正算出部 2として、上記Δｔをτで除し
て遅延残差ΔＷを演算するΔＷ演算手段 4と、周波数イ
ンデックスｋを発生するｋ発生手段 5と、前記ΔＷに、
前記ｋを乗算したΔＷ・ｋを演算するｋ乗算手段 6と、
全サンプル点数Ｎで、前記ΔＷ・ｋを除した値である位
相回転角Δωを演算するＮ除算手段 7と、前記Δωより
γ＝ｅｘｐ（ｊ２πΔω）を得るγ演算手段 8とを備え
るように構成する。｛図２参照｝ (3) 上記位相補正算出部 2において、周波数インデック
スｋを発生するｋ発生手段 5と、前記ΔＷに、前記ｋを
乗算したΔＷ・ｋを演算するｋ乗算手段 6と、全サンプ
ル点数Ｎで、前記ΔＷ・ｋを除した値である位相回転角
Δωを演算するＮ除算手段 7とからなる部分を、上記サ
ンプル点の数Ｎ＝２ⁿ としたときのｍ（ｍ≧ｎ）ビット
からなるカウンタ 21 と、上記カウンタ 21 の全ビット
をリバースするビットリバース部 22 と、上記ビットリ
バース部 22 でビットリバースされたｍビットの最上位
ビットの前に、小数点が存在するものとして、全ビット
をｋ／Ｎの値として抽出し、上記ΔＷを乗算して、ΔＷ
・ｋ／Ｎを得る演算手段 23 とで置き換えて、上記Δω
を得るように構成する。｛図３参照｝ (4) 前記γ演算手段 8として、前記Δωの小数点以下の
値に対応して、前記γ値が記載されたルックアップテー
ブルで構成する。

【００１９】

【作用】一般に、フーリエ変換は、(1) 式で与えられ
る。

【００２０】

【数１】

【００２１】ここで、ｇ（ｔ）をｇ（ｔ−Δｔ）と置き
換え、ｇ（ｔ−Δｔ）のフーリエ変換Ｇ’（ｆ）の演算
を行うと (2)式,(3)式が得られる。

【００２２】

【数２】

【００２３】上記ｇ（ｔ）とｇ（ｔ−Δｔ）との関係を
図４に示してある。図４において、時刻０，１，２，〜
は、サンプリングクロックの番号を示しており、ΔＷは
Δｔに対応し、ΔＷはサンプリング周期τに対するΔｔ
の割合を示している。Δｔを０≦Δｔ≦τとすれば、０
≦ΔＷ≦１である。

【００２４】上記において、(3) 式から明らかな如く、
時間がΔｔだけずれた信号をフーリエ変換したときに得
られるＧ’（ｆ）に、ｅｘｐ（ｊ２πｆΔｔ）（これ
を、位相補正データγと呼ぶ）を乗算すれば、基準時刻
ｔでのフーリエ変換Ｇ（ｆ）が得られることを示してい
る。

【００２５】そこで、本発明においては、サンプリング
点をｔ＝−Δｔ，τ−Δｔ，・・・・，（Ｎ−１）τ−
Δｔとしたときの信号の値ｇ（−Δｔ），ｇ（τ−Δ
ｔ），・・・，ｇ（（Ｎ−１）τ−Δｔ）に対して離散
フーリエ変換して、その周波数成分Ｇ’（ｆ）を、図１
のフーリエ変換部 1で算出する。

【００２６】上記フーリエ変換部 1で算出されたＧ’
（ｆ）と、上記位相補正データγ＝ｅｘｐ（ｊ２πｆΔ
ｔ）とを、位相補正乗算部 3で乗算することにより、上
記遅延時間Δｔ分を補正した、基準時刻ｔでのフーリエ
変換Ｇ（ｆ）を求めることができる。

【００２７】上記位相補正データγ＝ｅｘｐ（ｊ２πｆ
Δｔ）を得る為には、ｆΔｔを得る必要がある。ここ
で、Δωを位相回転角とし、Δω＝ｆΔｔとして、全サ
ンプリル点数をＮとしたとき、全サンプル点数Ｎの内の
周波数インデックスｋの周波数をｆとし、ｆｓ（＝１／
τ）をサンプリング周波数とすると、Δｆ＝ｆｓ／Ｎで
表したとき、上記周波数インデックスｋ、即ち、ｋ番目
の周波数ｆ＝ｋ・Δｆで表せるので 、ｆ＝ｋ・ｆｓ／Ｎ ＝ｋ・１／（τ・Ｎ）───────(4) となる。

【００２８】上記では、周波数インデックスｋについ
て、式の上でのみ表現したが、上記周波数インデックス
ｋの物理的な意味について、具体的に説明する。フーリ
エ変換は、一見不規則に見える信号から、ある規則性を
見出す手法であるスペクトル分析の方法である。

【００２９】このフーリエスペクトル分析は、上記不規
則な信号が、ある正弦波（又は、余弦波）の基本成分Δ
ｆと、その基本周波数Δｆの整数倍の成分をどの程度の
割合で含んでいるかを分析するもので、上記のように、
全サンプリング点数をＮとすると、各周波数成分の周波
数ｆ＝Δｆ，２Δｆ，〜，ｋΔｆ，〜，ＮΔｆで表せ
る。従って、上記周波数インデックスｋは、上記フーリ
エスペクトル分析を行ったときの周波数成分の基本周波
数Δｆのｋ倍目（即ち、前述のｋ番目の周波数）である
ことを意味していることになる。

【００３０】又、前述の図４の説明からも明らかなよう
に、 ΔＷ＝Δｔ／τ ─────────(5) 上記の(4) 式と、(5) 式から、 Δω＝ｆΔｔ＝ｋ・Δｔ／（τ・Ｎ）＝ｋ・ΔＷ・１／
Ｎ＝ΔＷ・ｋ／Ｎ─(6) 上記 (6)式を実現するため、図２のΔＷ演算手段 4でΔ
ｔをτで除してΔＷを求め、ｋ発生手段 5で、上記周波
数インデックスｋを発生して、ｋ乗算手段 6で、ΔＷに
上記周波数インデックスｋを乗算して、ΔＷ・ｋを求
め、Ｎ乗算手段 7で、上記ΔＷ・ｋをＮで除すことで、
上記 (6)式のΔωを求めることができる。このΔωを用
いて、γ演算手段 8により、前述の位相補正データγ＝
ｅｘｐ（ｊ２πΔω）を求めることができると、前述の
(3)式により、遅延時間Δｔ分の補正をしたフーリエ変
換を行うことができる。

【００３１】上記γ演算手段 8で、Δωより、位相補正
データγの値を求める場合、各Δωの小数点以下の値に
対するγの値をテーブルにしたルックアップテーブルを
用いることにより、短時間に、Δωよりγを求めること
ができる。上記Δωは位相回転値で、前述のように０〜
１の値を、０〜２πで表しているので、小数点以下の値
をとればよい。又、小数点以下の値のみでよいため、メ
モリの容量も少なくすることができる。

【００３２】次に、上記周波数インデックスｋを求める
ｋ発生手段 5について、以下に説明する。上記離散的フ
ーリエ変換の演算速度を高速化するのに、通常、高速フ
ーリエ変換(FFT) の手法が用いられる。

【００３３】上記高速フーリエ変換のアルゴリズムにつ
いては、文献「“高速フーリエ変換入門", "高速フーリ
エ変換(FFT) の使い方",安居院猛著, 廣済堂産報出版」
に詳しく記載されているので、ここでは、その詳細は省
略するが、本願発明に関連する部分について、その要点
を以下に示す。

【００３４】高速フーリエ変換(FFT) の手法には、種々
の方法が提案されているが、ここでは、説明の便宜上、
最も基本的な基数２の高速フーリエ変換の手法の内、時
間的間引きによる手法を基に、上記周波数インデックス
ｋの発生方法を説明する。

【００３５】高速フーリエ変換における高速化は、より
少ないデータ数に対する離散的フーリエ変換を繰り返し
提要して、全体の離散的フーリエ変換を求める手法であ
り、上記基数２の高速フーリエ変換は、最も簡単にフー
リエ演算、具体的には、所謂バタフライ演算が行える２
つのデータに対する離散的フーリエ変換を繰り返す手法
である。

【００３６】図５は、基数２の高速フーリエ変換を説明
する図であり、具体的には、サンプル値データの数Ｎ＝
８の場合の高速フーリエ変換での手法を示している。図
示されているように、ここでは、全ての演算が、２項間
の、所謂たすき掛け構成の演算によって成り立ってお
り、このたすき掛け演算は、基数２の高速フーリエ変換
における基本演算，又は、上記バタフライ演算と呼ばれ
ているものである。

【００３７】例えば、第１段目の最上段の基本演算は、
g(0)と、g(4)との演算｛ここで、g(0),g(4) の"0","4"
が、前述の周波数インデックスｋに対応している｝であ
り、第２段目では，図示のP₀とR₀,P₁ とR₁,Q₀ とS₀, お
よび Q₁ とS₁との演算であり、更に、第３段目では、E₀
とH₀,E₁ とH_1, ─などの基本演算から構成されている。

【００３８】このように、高速フーリエ変換(FFT) にお
いては、データの並べ換え操作を行う必要がある。 (図
５参照） 図６は、ビットリバース方法による、上記周波数インデ
ックスｋを求める手法を示した図である。

【００３９】上記図５の例では、サンプル値データg(n)
は、g(0),g(4),(g2),g(6),g(1),g(5),g(3),g(7) の順序
に並び替えられる。このように並び換えておけば、段数
により定まる基本演算対に対する、上記基本演算を施す
ことにより、最終結果は、最上段が直流成分 G₀ に、次
に、基本成分 G₁,そして、周波数の順に、高調波 G
_n(=7) が整列することになる。

【００４０】このように、最終結果を周波数の順に整列
するように高速フーリエ変換を行う為には、入力データ
の並び換えが必要となるが、この並べ換え操作の例を示
したものが、上記図６である。

【００４１】この並べ換え操作は、ｎ番目のデータ g
(n) において、ｎが２進数 b₁,b₂, ─,b_i と表示された
場合、２進数で前後のビット符号を入れ換えた b_i ,b
_i-1,─,b₂,b₁ となる番号へ移動する、所謂ビットリバ
ースの方法で簡単に行うことができる。

【００４２】このビットリバース（ビット入れ換え）の
方法で求めたサンプル値データg(n)の順序が、前述の周
波数インデックスｋである。前述の (6)式から明らかな
ように、Δω＝ｆΔｔ＝ｋ・Δｔ／（τ・Ｎ）＝ｋ・Δ
Ｗ・１／Ｎ＝ΔＷ・ｋ／Ｎである。但し、サンプル値デ
ータＮ＝２ⁿ とするサンプル値データＮは、ｎビットか
らなるデータとなる。

【００４３】ここで、Δω＝ｋ・ΔＷ／２ⁿ とすると、 Δω＝（ｋ・２^14-n）・ΔＷ／２¹⁴──────────(6) と変換される｛但し、ビットリバース用のカウンタとし
て、最高16K 点の高速フーリエ変換(FFT) を想定した場
合、上記ビットリバース演算用のカウンタは、14ビット
カウンタとなる｝ので、14ビットからなるカウンタをビ
ットリバースすることで、ΔＷとの乗算結果を、高速フ
ーリエ変換の点数Ｎに関係なく１４ビット右シフト｛即
ち、(6) 式で２¹⁴で除する) して、Δωを求めることが
できる。

【００４４】ここで、例えば、２５６（２⁸)点の高速フ
ーリエ変換(FFT) の場合を考えると、図３に示したよう
に、14ビットカウンタの上位６ビットは、常に“０”で
ある。そこで、このカウンタを全ビット反転して、上記
周波数インデックスｋを求めると、下位６ビットに
“０”が入るため、周波数インデックスｋを２⁶ 倍 (ｋ
・２^14-8) していることになる。

【００４５】つまり、カウンタ値を、高速フーリエ変換
(FFT) の点数Ｎに関係なく、全ビット反転することで、
Ｋ・２^14-nが得られるため、そのあとの演算は、上記式
(6)に示したように演算が行うこと、即ち、上記１４ビ
ット右シフトすることでΔωを求めてもよいし、上記１
４ビット右シフトすることなく、上記ビットリバースし
た状態で、その最上位ビットに、小数点があるものとし
て、Δωを求めるようにしてもよい。

【００４６】このようにシフト処理を行わない方法を採
用することで、少ないハードウェア量で、Δωを求める
ことができる。このように、本発明によれば、基準信号
と、Δｔの時間遅延をもった信号とを別々にフーリエ変
換した後、周波数領域で位相補正を行う処理となるた
め、サンプリング周期τ以下のような遅延時間の補正で
も精度よく行うことができる効果がある。又、周波数イ
ンデックスｋを高速フーリエ変換(FFT) の点数Ｎで除す
る演算を、最高の高速フーリエ変換点数Ｎ＝２^m を想定
したｍビットカウンタをビットリバースするだけで求め
る手法を用いることにより、少ないハードウェアで、上
記遅延時間の補正を行うことができるようになる。

【００４７】

【実施例】以下本発明の実施例を図面によって詳述す
る。前述の図１〜図３は、本発明の原理構成図であり、
図４は、ｇ（ｔ−Δｔ）とｇ（ｔ）との関係を説明する
図であり、図５は、基数２の高速フーリエ変換を説明す
る図であり、図６は、ビットリバース方法による周波数
インデックスｋを求める手法を示した図であり、図７〜
図９は、本発明の一実施例を示した図であって、図７は
全体の構成例を示し、図８は、遅延補正部の詳細を示
し、図９は、相関データ演算装置の構成例を示してお
り、図１０〜図１２はビットリバース装置の実施例を示
した図であって、図１０はシフタを使用した場合の構成
例を模式的に示し、図１１は、シフタを使用しない場合
の構成例を示し、図１２は、その実施例を示している。

【００４８】本発明においては、ディジタル信号を離散
フーリエ変換したときの遅延時間補正装置において、デ
ィジタル信号に対するサンプリング周期をτとし、サン
プリング点をｔ＝−Δｔ，τ−Δｔ，・・・，（Ｎ−
１）τ−Δｔとした合計Ｎ個のサンプル値ｇ（−Δ
ｔ），ｇ（τ−Δｔ），・・・，ｇ（（Ｎ−１）τ−Δ
ｔ）に対して離散フーリエ変換して、周波数成分Ｇ’
（ｆ）＝Ｇ（ｆ）・ｅｘｐ（−ｊ２πｆΔｔ）を得るフ
ーリエ変換部 1と、位相補正データγ＝ｅｘｐ（ｊ２π
ｆΔｔ）を算出する位相補正算出部 2と、前記Ｇ’
（ｆ）と、前記位相補正データγとを乗算する位相補正
乗算部 3が、本発明を実施するのに必要な手段である。
尚、全図を通して同じ符号は同じ対象物を示している。

【００４９】以下、図１〜図６を参照しながら、図７〜
図１２によって、本発明の離散フーリエ変換値の遅延時
間補正装置の構成と動作を説明する。先ず、図７におい
て、入力されるサンプル値データを、ｇ（ｔ）からＤだ
け進んだｇ（ｔ−Ｄ）とする。これは、ｇ（ｔ−Ｄ）を
基準にすると、ｇ（ｔ）はＤだけ遅れたデータというこ
とになる。

【００５０】そして、Ｄ＝t0＋Δｔとし、t0はサンプル
周期τの整数倍、具体的には、τ,２τ，─，（Ｎ−
１）τとし、Δｔは、１高速フーリエ変換(FFT) の周期
の間変化しない値とする。

【００５１】時間軸での遅延補正部 11 では、時間軸上
で、上記 t0 の補正を行う。これは、τの整数倍の時間
遅延は、この遅延補正部 11 で補正することを意味す
る。具体的には、図８に示した先入れ先だし(FIFO)バッ
ファ 20 を利用して、補正対象となるデータを、上記先
入れ先だし(FIFO)バッファ 20 に書き込み、所定の遅延
時間 (τの整数倍) 後に、順次読み出すことにより、ｇ
（ｔ−Ｄ) （Ｄ＝t0＋Δｔ）のデータよりt0だけ遅延し
たｇ(t−Δｔ) のデータを得ることができる。

【００５２】以後の遅延処理では、上記サンプル周期τ
より短いΔｔの遅延を補正することになる。そこで、図
７のＦＦＴ部 12 では、上記t0の補正が行われたｇ（ｔ
−Δｔ）を高速フーリエ変換法によってフーリエ変換を
行う。

【００５３】前述の図４では、上記ｇ（ｔ−Δｔ）と，
ｇ（ｔ）との関係を示している。サンプリング点は、ｔ
＝−Δｔ，τ−Δｔ，・・・，（Ｎ−１）τ−Δｔの合
計Ｎ個であり、そのサンプル値は、それぞれ、ｇ（−Δ
ｔ），ｇ（τ−Δｔ），・・・，ｇ（（Ｎ−１）τ−Δ
ｔ）となる。

【００５４】ｇ（ｔ−Δｔ）を基準とすると、ｇ（ｔ）
は、図４からも明らかなように、Δｔだけ遅延した位置
の値である。サンプリング点は、時刻（ｔ−Δｔ）であ
るのでｇ（ｔ−Δｔ）の値は得られるが、ｇ（ｔ）の値
は得られない。

【００５５】そこで、本発明では、前述の数式(2),(3)
で示した理論により、ｇ（ｔ−Δｔ）に対してフーリエ
変換を行い、その周波数成分Ｇ’（ｆ）を、上記ＦＦＴ
部 12 で得る。前述のように、図４において、ΔＷはΔ
ｔに対応し、ΔＷはサンプル周期τに対するΔｔの割
合、ΔＷ＝Δｔ／τを表している。従って、Δｔを０≦
Δｔ≦τとすれば、ΔＷは、０≦ΔＷ≦１である。

【００５６】上記のように求めたＧ’（ｆ）は、次の位
相補正乗算部 14 において、位相補正データ計算部 13
で求めた位相補正データγ＝ｅｘｐ（ｊ２πΔω）と乗
算されることで、所望のフーリエ変換後の周波数成分Ｇ
（ｆ）を得ることができる。

【００５７】図７に示した、上記位相補正データ計算部
13 では、入力Δｔを除算器 15 でτで除してΔＷを得
る。このΔＷに、前述の作用欄で詳細に説明したｋ発生
部 5、具体的には、後述するように、例えば、所定のＦ
ＦＴ点数Ｎ＝２ⁿ の場合には、最高ＦＦＴ点数Ｎ＝２^m
に対応したｍビットカウンタをビットリバースして、ｍ
−ｎビットシフトする手段で、周波数インデックスｋを
求めることができるので、上記ΔＷと、乗算器 6で乗算
して、ΔＷ・ｋを得る。

【００５８】次に、上記ΔＷ・ｋをシフタ 7を通して全
サンプル点数Ｎ(=２ⁿ ) で割ることで、前述のΔωを得
る。具体的には、Ｎ＝２ⁿ で表されるようにすると、上
記の除算は、シフタ 7で実現することができる。

【００５９】上記位相補正データ計算部 13 における位
相補正データメモリ 19 は、各Δωの小数点以下の値に
対して、γ＝ｅｘｐ（ｊ２πΔω）の値をルックアップ
テーブルの形で格納されているので、上記Δωの小数点
以下の値に対するγの値を、ルックアップテーブルに対
する読み出しで、位相補正データγを得ることができ
る。

【００６０】前述のように、上記位相補正乗算部 14
は、ＦＦＴ部 12 で求められたＧ’（ｆ）の値に、上記
位相補正データγを乗ずることにより、所望のＧ（ｆ）
を得る。Ｇ（ｆ）は、ｇ（ｔ）に対して離散的にフーリ
エ変換して得られた周波数成分である。

【００６１】図９は、上記の実施例を用いて、遅延時間
が異なる信号の相関を求める場合を示している。基準信
号ｈ（ｔ）に対して、被相関信号ｇ（ｔ−Ｄ）は、Ｄだ
け時間が進んでいる信号である。

【００６２】このｇ（ｔ−Ｄ）は、本実施例の遅延時間
補正装置 30 でフーリエ変換後、位相補正を行いＧ
（ｆ）を得て、上記基準信号ｈ（ｔ）をＦＦＴ部 32 で
フーリエ変換したＨ（ｆ）と、相関部 31 とで相関を取
る。

【００６３】精度よく遅延時間のあったデータの相関を
求めることにより、複数セグメントに渡って、各周波数
領域についての相関を取り累積することで、雑音に埋も
れた信号など、S/N 比の低い微弱信号、例えば、電波信
号などの検出及び解析を精度よく行うことができる。

【００６４】上記の図７で示した位相補正データ計算部
13 においては、ΔＷ・ｋをシフタ7でｎビットシフト
することで、FFT 点数Ｎ(=２ⁿ ）による除算、即ち、Δ
Ｗ・Ｋ／Ｎを実現した例で説明したが、この実施例で
は、ｋ発生部 5と、乗算器 6と、シフタ 7を必要とし、
ハードウェア量が大きくなるという問題が残る。

【００６５】この方法は、図３で説明したように、最高
ＦＦＴ点数Ｎ（＝２¹⁴＝16K 点) のフーリエ変換を求め
る為の、前述の周波数インデックスｋを求める場合、作
用欄で詳細に説明したように、１４(=m)ビットのカウン
タを設けておき、例えば、ＦＦＴ点数Ｎ（＝２⁸ =256
点) の周波数インデックスｋを求めるときには、図１０
に示したように、上記１４(=m)ビットのカウンタをビッ
トリバースして、「m-n=14-8=6」ビットのシフトをする
ことで、上記周波数インデックスｋを求め、、更に、n=
８ビットシフトすることで、ｋ／Ｎを求める手法を利用
している。

【００６６】然しながら、図１１に示したように、上記
１４ビットカウンタをビットリバースして、その儘のデ
ータ形式の最上位ビットの前に、小数点があるものとす
ると、図１０で説明したようなシフト動作をすることな
く、即ち、図７で説明したシフタ 18 を設けることな
く、上記ｋ／Ｎを求めることができることが分かる。

【００６７】図１２は、上記の点に着目して、ΔＷ・ｋ
／Ｎを求める場合の実施例を示したものである。総ＦＦ
Ｔ点数Ｎ（＝２^m ) のｍビットカウンタ 21 を設けて、
例えば、サンプリングクロックでカウントする。

【００６８】このカウンタ 21 を、次のビットリバース
部 22 、具体的には、図１２に示したビットリバース部
22 をワイヤで実現する。今のＦＦＴ点数Ｎ（＝２ⁿ )
であると、全ビットのビットリバースした値の下位ｍ−
ｎビットは、常に、“０”となる。そこで、本発明にお
いては、ビットリバース部 22 で、上記カウンタ 21 の
全ビットをビットリバースした値の最上位ビットの左
に、小数点があるものとして、上記ビットリバース部 2
2 の値を取り出すことで、上記ｋ／Ｎの演算を実現する
こきができる。このビットリバース部22 から求められ
たｋ／Ｎに乗算器 23 でΔＷと乗算することで、上記Δ
ωを求めることができる。

【００６９】従って、本実施例を使用することで，図７
で説明した位相補正データ計算部 13 の、ｋ発生部 5
と、乗算器 6と、シフタ 7の部分を、図１２に示したカ
ウンタ21 と、ワイヤ接続のみで実現したビットリバー
ス部 22 の乗算器 23 に置き換えることができ、ｋ発生
部でのｎ−ｍビットのシフタと、シフタ 7等を削減する
ことができ、ｋ／Ｎを、簡単で、回路規模の小さい回路
で得ることができる効果が得られる。

【００７０】

【発明の効果】以上、詳細に説明したように、本発明の
離散フーリエ変換値の遅延時間補正装置によれば、図７
からも明らかなように、被相関信号ｇ（ｔ−Ｄ）をフー
リエ変換した後、γ＝ｅｘｐ（ｊ２πｆΔｔ）で示され
る位相補正データを乗算するという、周波数領域での位
相補正を行うため、Δｔといったサンプリング周期τ以
下のような遅延時間の補正でも、精度よく遅延時間の補
正を行うことができる効果がある。

【００７１】又、上記γ＝ｅｘｐ（ｊ２πｆΔｔ）＝ｅ
ｘｐ（ｊ２πΔω）で示される位相補正データを作成す
る際、最高ＦＦＴ点数Ｎ（＝２^m ) に対応したｍビット
のカウンタと、該カウンタをワイヤ接続でビットリバー
スするビットリバース部を設けて、そのビットリバース
部の出力値の最上位ビットの前に、小数点があるものと
して、ＦＦＴ点数Ｎ（＝２^m ) に対応するｋ／Ｎの値を
求めるようにすることで、簡単で、回路規模の小さい回
路で、ｋ／Ｎを求めることができる効果がある。

【図面の簡単な説明】

【図１】本発明の原理構成図（その１）

【図２】本発明の原理構成図（その２）

【図３】本発明の原理構成図（その３）

【図４】ｇ（ｔ）とｇ（ｔ−Δｔ）との関係を説明する
図

【図５】基数２の高速フーリエ変換を説明する図

【図６】ビットリバース方法による周波数インデックス
ｋを求める手法を示した図

【図７】本発明の一実施例を示した図（その１）

【図８】本発明の一実施例を示した図（その２）

【図９】本発明の一実施例を示した図（その３）

【図１０】ビットリバース装置の実施例を示した図（そ
の１）

【図１１】ビットリバース装置の実施例を示した図（そ
の２）

【図１２】ビットリバース装置の実施例を示した図（そ
の３）

【図１３】２つのパラボラアンテナで受信した電波の信
号処理を説明する図

【符号の説明】

1 フーリエ変換部 2 位相補正算
出部 3 位相補正乗算部 4 ΔＷ演算手
段 5 ｋ発生手段, ｋ発生部 6 ｋ乗算手段 7 Ｎ除算手段（シフタ） 8 γ演算手段 11 時間軸での遅延補正部 12,32 ＦＦＴ部 13 位相補正デ
ータ計算部 14 位相補正乗算部 15 除算器 23 乗算器 19 位相補正データメモリ 20 先入れ先出し(FIFO)バッファ 21 ｍビットカウンタ，カウンタ 22 ビットリバース部 30 遅延時間補正装置 31 相関部 Ｎ ＦＦＴ点数，サンプル点数 Ｇ（ｆ） 基準データの周波数成分 Ｇ’（ｆ） Δｔ遅延したデータの周波数成分 γ 位相補正データ ｇ（ｔ） サンプル値 ｋ 周波数インデックス Δｔ 遅延時間 ΔＷ 遅延残差 Δω 位相回転角 ｆｓ サンプリング周波数 τ サンプル周期，サンプリング周期 Δｆ 周波数チャネルの間隔（フーリエスペクトルの
基本成分）

───────────────────────────────────────────────────── フロントページの続き (72)発明者 久保 慎一 神奈川県川崎市中原区上小田中1015番地 富士通株式会社内

Claims (4)

【特許請求の範囲】
1. 【請求項１】アナログ信号に対するサンプリング周期を
   τとし、サンプリング点をｔ＝−Δｔ，τ−Δｔ，・・
   ・，（Ｎ−１）τ−Δｔとした合計Ｎ個のサンプル値ｇ
   （−Δｔ），ｇ（τ−Δｔ），・・・，ｇ（（Ｎ−１）
   τ−Δｔ）に対して離散フーリエ変換して、周波数成分
   Ｇ’（ｆ）＝Ｇ（ｆ）・ｅｘｐ（−ｊ２πｆΔｔ）を得
   るフーリエ変換部(1) と、位相補正データγ＝ｅｘｐ
   （ｊ２πｆΔｔ）を算出する位相補正算出部(2) と、前
   記Ｇ’（ｆ）と、前記位相補正データγとを乗算する位
   相補正乗算部(3) とを備えたことを特徴とする離散フー
   リエ変換値の遅延時間補正装置。
2. 【請求項２】上記位相補正算出部(2) として、上記Δｔ
   をτで除して遅延残差ΔＷを演算するΔＷ演算手段(4)
   と、周波数インデックスｋを発生するｋ発生手段(5)
   と、前記ΔＷに、前記ｋを乗算したΔＷ・ｋを演算する
   ｋ乗算手段(6) と、全サンプル点数Ｎで、前記ΔＷ・ｋ
   を除した値である位相回転角Δωを演算するＮ除算手段
   (7) と、前記Δωよりγ＝ｅｘｐ（ｊ２πΔω）を得る
   γ演算手段(8) とを備えたことを特徴とする請求項１に
   記載の離散フーリエ変換値の遅延時間補正装置。
3. 【請求項３】上記位相補正算出部(2) において、周波数
   インデックスｋを発生するｋ発生手段(5) と、前記ΔＷ
   に、前記ｋを乗算したΔＷ・ｋを演算するｋ乗算手段
   (6) と、全サンプル点数Ｎで、前記ΔＷ・ｋを除した値
   である位相回転角Δωを演算するＮ除算手段(7) とから
   なる部分を、 上記サンプル点の数Ｎ＝２ⁿ としたと、のｍ（ｍ≧ｎ）
   ビットからなるカウンタ(21)と、上記カウンタ(21)の全
   ビットをリバースするビットリバース部(22)と、上記ビ
   ットリバース部(22)でビットリバースされたｍビットの
   最上位ビットの前に、小数点が存在するものとして、全
   ビットをｋ／Ｎの値として抽出し、上記ΔＷを乗算し
   て、ΔＷ・ｋ／Ｎを得る演算手段とで置き換えて、上記
   Δωを得ることを特徴とする請求項１に記載の離散フー
   リエ変換値の遅延時間補正装置。
4. 【請求項４】前記γ演算手段(8) として、前記Δωの小
   数点以下の値に対応して、前記γ値が記載されたルック
   アップテーブルよりになることを特徴とする請求項２，
   又は、３に記載の離散フーリエ変換値の遅延時間補正装
   置。