JPH0340876B2 - - Google Patents
Info
- Publication number
- JPH0340876B2 JPH0340876B2 JP5904682A JP5904682A JPH0340876B2 JP H0340876 B2 JPH0340876 B2 JP H0340876B2 JP 5904682 A JP5904682 A JP 5904682A JP 5904682 A JP5904682 A JP 5904682A JP H0340876 B2 JPH0340876 B2 JP H0340876B2
- Authority
- JP
- Japan
- Prior art keywords
- circuit
- correlation coefficient
- maximum value
- spectrum
- power spectrum
- 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.)
- Expired
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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)
- Complex Calculations (AREA)
- Image Analysis (AREA)
Description
【発明の詳細な説明】
〔発明の技術分野〕
この発明は電力スペクトルの相互相関係数を計
算することにより、入力信号の中から基準となる
電力スペクトルと同等のスペクトルを持つ信号を
検出する相関処理装置に関する。
算することにより、入力信号の中から基準となる
電力スペクトルと同等のスペクトルを持つ信号を
検出する相関処理装置に関する。
第１図は従来の相関処理装置を示すものであ
る。入力信号iはアナログ／デイジタル（Ａ／
Ｄ）変換回路１１に供給され、この変換回路１１
において所定周期のサンプリングパルス信号sに
より、デイジタル時系列信号に変換される。この
サンプルデータは第１の記憶回路１２に供給され
記憶される。この記憶されるサンプル数は2N個
とする。この記憶回路１２より読出されたサンプ
ルデータは順次電力スペクトル計算回路１３に供
給され、例えばFFT（Fast Fourier Transform）
方式によつてＮ個の電力スペクトルデータxi（ｉ
＝０、１、…、Ｎ−１）が求められる。この電力
スペクトルデータxiは第２の記憶回路１４に供給
され記憶される。
る。入力信号iはアナログ／デイジタル（Ａ／
Ｄ）変換回路１１に供給され、この変換回路１１
において所定周期のサンプリングパルス信号sに
より、デイジタル時系列信号に変換される。この
サンプルデータは第１の記憶回路１２に供給され
記憶される。この記憶されるサンプル数は2N個
とする。この記憶回路１２より読出されたサンプ
ルデータは順次電力スペクトル計算回路１３に供
給され、例えばFFT（Fast Fourier Transform）
方式によつてＮ個の電力スペクトルデータxi（ｉ
＝０、１、…、Ｎ−１）が求められる。この電力
スペクトルデータxiは第２の記憶回路１４に供給
され記憶される。
一方、第３の記憶回路１５には予じめ検出した
い信号に対応する基準の電力スペクトルデータyi
（ｉ＝０、１、…、Ｍ−１；Ｍ＜Ｎ）が記憶され
ている。この記憶回路１５および１４に記憶され
たデータはそれぞれ順次読出され、相互相関計算
回路１６に供給される。この計算回路１６では次
の計算が行なわれる。
い信号に対応する基準の電力スペクトルデータyi
（ｉ＝０、１、…、Ｍ−１；Ｍ＜Ｎ）が記憶され
ている。この記憶回路１５および１４に記憶され
たデータはそれぞれ順次読出され、相互相関計算
回路１６に供給される。この計算回路１６では次
の計算が行なわれる。
但し、
ｋ＝０、１、…、Ｎ−１
＝１／Ｍk+M-1
〓i=k
xi
＝１／ＭM-1
〓i=0
yi
この計算によつて得られた相関係数γ(k)は判定
回路１７に供給され、この判定回路１７において
相関係数γ(k)の最大値γ(k)maxが判定される。
そして、このｋの値により基準のスペクトルパタ
ーンとその電力スペクトルパターンが一致する信
号の中心周波数が求められる。第２図ａ，ｂ，ｃ
はそれぞれxi、yi、γ(k)の一例を示すものであ
る。
回路１７に供給され、この判定回路１７において
相関係数γ(k)の最大値γ(k)maxが判定される。
そして、このｋの値により基準のスペクトルパタ
ーンとその電力スペクトルパターンが一致する信
号の中心周波数が求められる。第２図ａ，ｂ，ｃ
はそれぞれxi、yi、γ(k)の一例を示すものであ
る。
ところで、入力信号の電力スペクトルxiは周波
数分解能Δ（但し、Δ＝s／2N）毎の離散的なデ
ー タであるため、第３図ａ，ｂに示す如く信号の中
心周波数0の値によつて本来は同一スペクトルの
筈が、かなり異なつたスペクトルデータとなる可
能性を有している。尚、第３図ａ，ｂ中点線は連
続スペクトルエンベロープを示し、 第３図ａは0＝（ｍ＋１／２）×Δの場合 第３図ｂは0＝ｍ×Δの場合 （但し、ｍは正整数） であり、0の値によつてこの両者を両極端として
その中間のスペクトルパターンをとることにな
る。
数分解能Δ（但し、Δ＝s／2N）毎の離散的なデ
ー タであるため、第３図ａ，ｂに示す如く信号の中
心周波数0の値によつて本来は同一スペクトルの
筈が、かなり異なつたスペクトルデータとなる可
能性を有している。尚、第３図ａ，ｂ中点線は連
続スペクトルエンベロープを示し、 第３図ａは0＝（ｍ＋１／２）×Δの場合 第３図ｂは0＝ｍ×Δの場合 （但し、ｍは正整数） であり、0の値によつてこの両者を両極端として
その中間のスペクトルパターンをとることにな
る。
このことは、即ち基準のスペクトルパターンが
例えば第３図ａの形で記憶されているとすると、
第３図に示すａとｂの中間パターンの場合にはそ
れだけ相関係数の最大値γ(k)maxが変動するこ
とになる。このように、従来の装置は信号の中心
周波数と周波数分解能との相対関係によつて相関
係数が大きく変動するという欠点を有している。
例えば第３図ａの形で記憶されているとすると、
第３図に示すａとｂの中間パターンの場合にはそ
れだけ相関係数の最大値γ(k)maxが変動するこ
とになる。このように、従来の装置は信号の中心
周波数と周波数分解能との相対関係によつて相関
係数が大きく変動するという欠点を有している。
この発明は上記事情に基づいてなされたもの
で、その目的とするところは入力信号の中心周波
数が如何なる場合であつても所定の相関係数を得
ることが可能な相関処理装置を提供しようとする
ものである。
で、その目的とするところは入力信号の中心周波
数が如何なる場合であつても所定の相関係数を得
ることが可能な相関処理装置を提供しようとする
ものである。
この発明は電力スペクトル計算回路によつて決
まる周波数分解能Δの範囲内にわたりΔ／ｎず
つサンプリング位置が変化されたｎ個の基準電力
スペクトルを予じめ記憶しておき、この基準電力
スペクトルそれぞれと入力信号より求めた電力ス
ペクトルとの相関係数の最大値を求め、この最大
値により相関係数の計算範囲を制御するととも
に、求められた最大値群よりさらに最大値を求め
ることにより変動の少ない相関係数を得るもので
ある。
まる周波数分解能Δの範囲内にわたりΔ／ｎず
つサンプリング位置が変化されたｎ個の基準電力
スペクトルを予じめ記憶しておき、この基準電力
スペクトルそれぞれと入力信号より求めた電力ス
ペクトルとの相関係数の最大値を求め、この最大
値により相関係数の計算範囲を制御するととも
に、求められた最大値群よりさらに最大値を求め
ることにより変動の少ない相関係数を得るもので
ある。
以下、この発明の一実施例について図面を参照
して説明する。尚、第４図中第１図と同一部分に
は同一符号を付す。
して説明する。尚、第４図中第１図と同一部分に
は同一符号を付す。
第４図において、入力信号iの電力スペクトル
xiは第１図と同様の動作によつて得られる。一
方、第３の記憶回路４１1〜４１oには周波数領域
におけるサンプリング位置が電力スペクトル計算
回路１３によつて決まる周波数分解能Δの範囲
内においてΔ／ｎずつ変化された基準電力スペ
クトルy1i〜yoiがそれぞれ記憶されている。この
基準電力スペクトルy1i〜yoiは切換回路４２によ
つて順次読出され、第２の記憶回路１４より読出
された電力スペクトルxiとともに相互相関計算回
路１６に供給される。即ち、最初に記憶回路４１
１より基準電力スペクトルy1iが読出され、このス
ペクトルy1iとスペクトルxiとの相関係数が求めら
れる。第１の判定回路１７ではこの中から最大値
γ1（k1）が判定される。次に、記憶回路４１2より
スペクトルy2iが読出され、これとスペクトルxiと
の相関係数が求められる。このとき、求める信号
の中心周波数に対応する相関係数はk1の近傍の値
であることは明白である。計算範囲制御回路４３
は記憶回路４１2〜４１oより読出される基準電力
スペクトルy2i〜yoiに対する相関係数の計算範囲
をk1の値に基づいて制御するものである。即ち、
基準電力スペクトルy2i〜yoiについては相互相関
計算回路１６における(1)式の計算をｋ＝０、１、
…、Ｎ−１の全範囲について計算せず、k1を中心
として例えば幅2Mの範囲内で計算するよう制御
する。このようにして、順次記憶回路４１2〜４
１oより基準電力スペクトルy2i〜yoiが読出されて
相関係数の最大値が求められ、γ1（k1）、γ2（k2）、
…、γo（ko）のｎ個の最大値群が求められる。こ
の結果、サンプリング位置をΔ／ｎずつ変化さ
せたｎ個の基準スペクトルパターンと入力信号と
のｎ個の相関係数の最大値が得られたことにな
る。したがつて、このｎ個の最大値群の中から第
２の判定回路４４によつて最大値γ(k)maxを判
定すれば、基準の周波数パターンに最も近いとき
の相関係数を得ることができる。
xiは第１図と同様の動作によつて得られる。一
方、第３の記憶回路４１1〜４１oには周波数領域
におけるサンプリング位置が電力スペクトル計算
回路１３によつて決まる周波数分解能Δの範囲
内においてΔ／ｎずつ変化された基準電力スペ
クトルy1i〜yoiがそれぞれ記憶されている。この
基準電力スペクトルy1i〜yoiは切換回路４２によ
つて順次読出され、第２の記憶回路１４より読出
された電力スペクトルxiとともに相互相関計算回
路１６に供給される。即ち、最初に記憶回路４１
１より基準電力スペクトルy1iが読出され、このス
ペクトルy1iとスペクトルxiとの相関係数が求めら
れる。第１の判定回路１７ではこの中から最大値
γ1（k1）が判定される。次に、記憶回路４１2より
スペクトルy2iが読出され、これとスペクトルxiと
の相関係数が求められる。このとき、求める信号
の中心周波数に対応する相関係数はk1の近傍の値
であることは明白である。計算範囲制御回路４３
は記憶回路４１2〜４１oより読出される基準電力
スペクトルy2i〜yoiに対する相関係数の計算範囲
をk1の値に基づいて制御するものである。即ち、
基準電力スペクトルy2i〜yoiについては相互相関
計算回路１６における(1)式の計算をｋ＝０、１、
…、Ｎ−１の全範囲について計算せず、k1を中心
として例えば幅2Mの範囲内で計算するよう制御
する。このようにして、順次記憶回路４１2〜４
１oより基準電力スペクトルy2i〜yoiが読出されて
相関係数の最大値が求められ、γ1（k1）、γ2（k2）、
…、γo（ko）のｎ個の最大値群が求められる。こ
の結果、サンプリング位置をΔ／ｎずつ変化さ
せたｎ個の基準スペクトルパターンと入力信号と
のｎ個の相関係数の最大値が得られたことにな
る。したがつて、このｎ個の最大値群の中から第
２の判定回路４４によつて最大値γ(k)maxを判
定すれば、基準の周波数パターンに最も近いとき
の相関係数を得ることができる。
このような構成とすることにより入力信号iの
中心周波数0が周波数分解能Δの中のどの位置
にあつても常に変動の少ない相関係数γ(k)max
を得ることが可能となる。
中心周波数0が周波数分解能Δの中のどの位置
にあつても常に変動の少ない相関係数γ(k)max
を得ることが可能となる。
以上、詳述したようにこの発明によれば、入力
信号の中心周波数が如何なる場合であつても所定
の相関係数を得ることが可能な相関処理装置を提
供できる。
信号の中心周波数が如何なる場合であつても所定
の相関係数を得ることが可能な相関処理装置を提
供できる。
第１図は従来の相関処理装置の一例を示す構成
図、第２図ａ，ｂ，ｃおよび第３図ａ，ｂはそれ
ぞれ従来およびこの発明の動作を説明するために
示す図、第４図はこの発明に係わる相関処理装置
の一実施例を示す構成図である。 １１……Ａ／Ｄ変換回路、１３……電力スペク
トル計算回路、１６……相互相関計算回路、１７
……第１の判定回路、４１1〜４１o……第３の記
憶回路、４２……切換回路、４３……計算範囲制
御回路、４４……第２の判定回路。
図、第２図ａ，ｂ，ｃおよび第３図ａ，ｂはそれ
ぞれ従来およびこの発明の動作を説明するために
示す図、第４図はこの発明に係わる相関処理装置
の一実施例を示す構成図である。 １１……Ａ／Ｄ変換回路、１３……電力スペク
トル計算回路、１６……相互相関計算回路、１７
……第１の判定回路、４１1〜４１o……第３の記
憶回路、４２……切換回路、４３……計算範囲制
御回路、４４……第２の判定回路。
Claims (1)
- １ 入力信号を所定のサンプリングパルス信号で
Ａ／Ｄ変換する回路と、このサンプリングされた
信号より電力スペクトルを求める回路と、この回
路の周波数分解能Δの範囲内にわたりΔ／ｎず
つサンプリング位置が変化されたｎ個の基準電力
スペクトルを記憶する回路と、この基準電力スペ
クトルを順次読出す回路と、この読出された基準
電力スペクトルおよび前記求められた電力スペク
トルより相互相関係数を求める回路と、この求め
られた相関係数より最大値を求める第１の判定回
路と、この求められた最大値により前記相互相関
係数の計算範囲を制御する回路と、前記第１の判
定回路より出力される基準電力スペクトルを変え
て得られたｎ個の最大値群の中からさらに最大値
を求める第２の判定回路とを具備したことを特徴
とする相関処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5904682A JPS58176778A (ja) | 1982-04-09 | 1982-04-09 | 相関処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5904682A JPS58176778A (ja) | 1982-04-09 | 1982-04-09 | 相関処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS58176778A JPS58176778A (ja) | 1983-10-17 |
| JPH0340876B2 true JPH0340876B2 (ja) | 1991-06-20 |
Family
ID=13101975
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5904682A Granted JPS58176778A (ja) | 1982-04-09 | 1982-04-09 | 相関処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS58176778A (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011003948A (ja) | 2009-06-16 | 2011-01-06 | Sony Corp | データ処理装置及びデータ処理方法、受信装置及び受信方法、同期検出装置及び同期検出方法、並びにコンピューター・プログラム |
-
1982
- 1982-04-09 JP JP5904682A patent/JPS58176778A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS58176778A (ja) | 1983-10-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH02500225A (ja) | ループバック導体グリッドを有するデジタイザシステム | |
| EP0072706B1 (en) | Sound signal processing apparatus | |
| IE43171L (en) | Digital device | |
| JPH04230867A (ja) | 自己整列サンプリング装置 | |
| JPH0738587B2 (ja) | 多重アナログ信号の増幅・サンプル装置 | |
| JPH0340876B2 (ja) | ||
| JP2867769B2 (ja) | 音響測定方法およびその装置 | |
| US3715509A (en) | Method and means for providing resolution level selection in a spectrum analyzer | |
| US4426624A (en) | Device and method for amplifying and sampling multiplexed signals | |
| JPH0340872B2 (ja) | ||
| US4351032A (en) | Frequency sensing circuit | |
| JPH0340875B2 (ja) | ||
| JP3139803B2 (ja) | インパルス応答測定装置 | |
| US4916646A (en) | System for signal to noise separation for electrochemical instruments | |
| JP3272865B2 (ja) | 電力系統ディジタル保護制御装置のａ／ｄ変換方法 | |
| JPH0340874B2 (ja) | ||
| JPH0340873B2 (ja) | ||
| KR920020220A (ko) | 고저항 지락사고 검출장치 | |
| RU2099720C1 (ru) | Цифровой анализатор спектра | |
| JPS5958912A (ja) | アナログ入力装置 | |
| JP2589864B2 (ja) | 信号比較装置 | |
| RU2133041C1 (ru) | Способ определения спектра электрических сигналов | |
| JPH0690637B2 (ja) | 補間方法 | |
| SU924858A1 (ru) | Аналого-цифровой преобразователь | |
| JPS58211275A (ja) | 集録デ−タの保存方法 |