JPH0340872B2 - - Google Patents
Info
- Publication number
- JPH0340872B2 JPH0340872B2 JP5904282A JP5904282A JPH0340872B2 JP H0340872 B2 JPH0340872 B2 JP H0340872B2 JP 5904282 A JP5904282 A JP 5904282A JP 5904282 A JP5904282 A JP 5904282A JP H0340872 B2 JPH0340872 B2 JP H0340872B2
- Authority
- JP
- Japan
- Prior art keywords
- circuit
- power spectrum
- reference power
- correlation coefficient
- maximum value
- 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
【発明の詳細な説明】
〔発明の技術分野〕
この発明は電力スペクトルの相互相関係数を計
算することにより、入力信号の中から基準となる
電力スペクトルと同等のスペクトルを持つ信号を
検出する相関処理装置に関する。[Detailed Description of the Invention] [Technical Field of the Invention] The present invention is a correlation method for detecting a signal having a spectrum equivalent to a reference power spectrum from an input signal by calculating a cross-correlation coefficient of a power spectrum. It relates to a processing device.
〔発明の技術的背景とその問題点〕
第１図は従来の相関処理装置を示すものであ
る。入力信号iはアナログ／デイジタル（Ａ／
Ｄ）変換回路１１に供給され、この変換回路１１
において所定周期のサンプリングパルス信号sに
より、デイジタル時系列信号に変換される。この
サンプルデータは第１の記憶回路１２に供給され
記憶される。この記憶されるサンプル数は2N個
とする。この記憶回路１２より読出されたサンプ
ルデータは順次電力スペクトル計算回路１３に供
給され、例えばFFT（Fast Fourier Transform）
方式によつてＮ個の電力スペクトルデータxi（ｉ
＝０、１、…、Ｎ−１）が求められる。この電力
スペクトルデータxiは第２の記憶回路１４に供給
され記憶される。[Technical background of the invention and its problems] FIG. 1 shows a conventional correlation processing device. Input signal i is analog/digital (A/
D) is supplied to the conversion circuit 11, and this conversion circuit 11
At the sampling pulse signal s of a predetermined period, the signal is converted into a digital time series signal. This sample data is supplied to the first storage circuit 12 and stored therein. The number of samples stored is 2N. The sample data read out from the storage circuit 12 is sequentially supplied to the power spectrum calculation circuit 13, and is processed using, for example, FFT (Fast Fourier Transform).
N power spectrum data x i (i
=0, 1,...,N-1) are obtained. This power spectrum data x i is supplied to the second storage circuit 14 and stored therein.
一方、第３の記憶回路１５には予じめ検出した
い信号に対応する基準の電力スペクトルデータyi
（ｉ＝０、１、…、Ｍ−１；Ｍ＜Ｎ）が記憶され
ている。この記憶回路１５および１４に記憶され
たデータはそれぞれ順次読出され、相互相関計算
回路１６に供給される。この計算回路１６では次
の計算が行なわれる。 On the other hand, the third storage circuit 15 stores reference power spectrum data y i corresponding to the signal to be detected in advance.
(i=0, 1, . . . , M-1; M<N) are stored. The data stored in storage circuits 15 and 14 are read out in sequence and supplied to cross-correlation calculation circuit 16. This calculation circuit 16 performs the following calculation.
但し、
ｋ＝０、１、…、Ｎ−１
＝１／Ｍk+M+1
〓i=k
xi
＝１／ＭM-1
〓i=0
yi
この計算によつて得られた相関係数ｒ(k)は判定
回路１７に供給され、この判定回路１７において
相関係数ｒ(k)の最大値ｒ(k)naxが判定される。そ
して、このｋの値により基準のスペクトルパター
ンとその電力スペクトルパターンが一致する信号
の中心周波数が求められる。第２図ａ，ｂ，ｃは
それぞれxi、yi、ｒ(k)の一例を示すものである。 However, k=0, 1,..., N-1 = 1/M k+M+1 〓 i=k x i = 1/M M-1 〓 i=0 y iThe phase obtained by this calculation The correlation coefficient r(k) is supplied to a determination circuit 17, which determines the maximum value r(k) nax of the correlation coefficient r(k). Then, from this value of k, the center frequency of the signal whose power spectrum pattern matches the reference spectrum pattern is determined. Figures 2a, b, and c show examples of x i , y i , and r(k), respectively.
ところで、入力信号の電力スペクトルxiは周波
数分解能Δ（但し、Δ＝s／2N）毎の離散的なデ
ー
タであるため、第３図ａ，ｂに示す如く信号の中
心周波数0の値によつて本来は同一スペクトルの
筈が、かなり異なつたスペクトルデータとなる可
能性を有している。尚、第３図ａ，ｂ中点線は連
続スペクトルエンベロープを示し、
第３図ａは0＝（ｍ＋１／２）×Δの場合
第３図ｂは0＝ｍ×Δの場合
（但し、ｍは正整数）
であり、0の値によつてこの両者を両極端として
その中間のスペクトルパターンをとることにな
る。 By the way, since the power spectrum x i of the input signal is discrete data for each frequency resolution Δ (Δ = s / 2N), it is determined by the value of the center frequency 0 of the signal as shown in Figure 3 a and b. Therefore, what should originally be the same spectrum has the possibility of becoming quite different spectrum data. The dotted line between a and b in Figure 3 shows the continuous spectral envelope. (a positive integer), and a value of 0 makes these two extremes and a spectral pattern in between.
このことは、即ち基準のスペクトルパターンが
例えば第３図ａの形で記憶されているとすると、
第３図に示すａとｂの中間パターンの場合にはそ
れだけ相関係数の最大値ｒ(k)naxが変動すること
になる。このように、従来の装置は信号の中心周
波数と周波数分解能との相対関係によつて相関係
数が大きく変動するという欠点を有している。 This means that if the reference spectral pattern is stored, for example, in the form shown in Figure 3a,
In the case of an intermediate pattern between a and b shown in FIG. 3, the maximum value r(k) nax of the correlation coefficient varies accordingly. As described above, the conventional apparatus has the drawback that the correlation coefficient varies greatly depending on the relative relationship between the center frequency of the signal and the frequency resolution.
この発明は上記事情に基づいてなされたもの
で、その目的とするところは入力信号の中心周波
数が如何なる場合であつても所定の相関係数を得
ることが可能な相関処理装置を提供しようとする
ものである。
The present invention has been made based on the above circumstances, and its purpose is to provide a correlation processing device that can obtain a predetermined correlation coefficient no matter what the center frequency of the input signal. It is something.
この発明は電力スペクトル計算回路によつて決
まる周波数分解能Δの範囲内にわたりΔ／ｎず
つサンプリング位置が変化されたｎ個の基準電力
スペクトルを予じめ記憶しておき、この基準電力
スペクトルそれぞれと入力信号より求めた電力ス
ペクトルとの相関係数の最大値を求め、この最大
値群よりさらに最大値を求めることにより変動の
少ない相関係数を得るものである。
This invention stores in advance n reference power spectra whose sampling positions are changed by Δ/n over a range of frequency resolution Δ determined by a power spectrum calculation circuit, and inputs each of these reference power spectra. The maximum value of the correlation coefficient with the power spectrum determined from the signal is determined, and the maximum value is determined from this group of maximum values to obtain a correlation coefficient with little fluctuation.
以下、この発明の一実施例について図面を参照
して説明する。尚、第４図中第１図と同一部分に
は同一符号を付す。
An embodiment of the present invention will be described below with reference to the drawings. Note that the same parts in FIG. 4 as in FIG. 1 are given the same reference numerals.
第４図において、入力信号iの電力スペクトル
xiは第１図と同様の動作によつて得られる。一
方、第３の記憶回路４１1〜４１oには周波数領域
におけるサンプリング位置が電力スペクトル計算
回路１３によつて決まる周波数分解能Δの範囲
内においてΔ／ｎずつ変化された基準電力スペ
クトルy1i〜yoiがそれぞれ記憶されている。この
基準電力スペクトルy1i〜yoiは切換回路４２によ
つて順次読出され、第２の記憶回路１４より読出
された電力スペクトルxiとともに相互相関計算回
路１６に供給される。即ち、最初に記憶回路４１
１より基準電力スペクトルy1iが読出され、このス
ペクトルy1iとスペクトルxiとの相関係数が求めら
れる。第１の判定回路１７ではこの中から最大値
r1（k1）が判定される。次に、記憶回路４１2より
スペクトルy2iが読出され、これとスペクトルxiと
の相関係数の最大値r2（k2）が同様にして求めら
れる。以下、順次第３の記憶回路より基準電力ス
ペクトルが読出されて相関係数の最大値が求めら
れ、r1（k1）、r2（k2）、…、ro（ko）のｎ個の最大値
群が得られる。この結果、サンプリング位置を
Δ／ｎずつ変化させたｎ個の基準スペクトルパ
ターンと入力信号とのｎ個の相関係数の最大値が
得られたことになる。したがつて、このｎ個の最
大値群の中から第２の判定回路４３によつて最大
値ｒ(k)naxを判定すれば、基準の周波数パターン
に最も近いときの相関係数を得ることができる。 In Figure 4, the power spectrum of input signal i
x i can be obtained by the same operation as in FIG. On the other hand, the third storage circuits 41 1 to 41 o store reference power spectra y 1i to y whose sampling positions in the frequency domain are varied by Δ/n within the range of frequency resolution Δ determined by the power spectrum calculation circuit 13 . oi is memorized respectively. The reference power spectra y 1i to y oi are sequentially read out by the switching circuit 42 and supplied to the cross-correlation calculation circuit 16 together with the power spectrum x i read out from the second storage circuit 14 . That is, first the memory circuit 41
1 , the reference power spectrum y 1i is read out, and the correlation coefficient between this spectrum y 1i and the spectrum x i is determined. The first judgment circuit 17 selects the maximum value from among these.
r 1 (k 1 ) is determined. Next, the spectrum y 2i is read out from the storage circuit 41 2 , and the maximum value r 2 (k 2 ) of the correlation coefficient between this and the spectrum x i is found in the same manner. Thereafter, the reference power spectrum is sequentially read out from the storage circuit 3 and the maximum value of the correlation coefficient is determined, and n of r 1 (k 1 ), r 2 (k 2 ), ..., r o (k o ) is calculated. A group of maximum values is obtained. As a result, the maximum value of the n correlation coefficients between the input signal and the n reference spectrum patterns whose sampling positions are changed by Δ/n is obtained. Therefore, if the second determination circuit 43 determines the maximum value r(k) nax from among these n maximum value groups, it is possible to obtain the correlation coefficient when it is closest to the reference frequency pattern. Can be done.
このような構成とすることにより入力信号iの
中心周波数0が周波数分解能Δの中のどの位置
にあつても常に変動の少ない相関係数ｒ(k)naxを
得ることが可能となる。 With such a configuration, it is possible to always obtain a correlation coefficient r(k) nax with little variation no matter where the center frequency 0 of the input signal i is located within the frequency resolution Δ.
以上、詳述したようにこの発明によれば、入力
信号の中心周波数が如何なる場合であつても所定
の相関係数を得ることが可能な相関処理装置を提
供できる。
As described in detail above, according to the present invention, it is possible to provide a correlation processing device that can obtain a predetermined correlation coefficient regardless of the center frequency of the input signal.
第１図は従来の相関処理装置の一例を示す構成
図、第２図ａ，ｂ，ｃおよび第３図ａ，ｂはそれ
ぞれ従来およびこの発明の動作を説明するために
示す図、第４図はこの発明に係わる相関処理装置
の一実施例を示す構成図である。
１１……Ａ／Ｄ変換回路、１３……電力スペク
トル計算回路、１６……相互相関計算回路、１７
……第１の判定回路、４１1〜４１o……第３の記
憶回路、４２……切換回路、４３……第２の判定
回路。
FIG. 1 is a block diagram showing an example of a conventional correlation processing device, FIG. FIG. 1 is a configuration diagram showing an embodiment of a correlation processing device according to the present invention. 11...A/D conversion circuit, 13...Power spectrum calculation circuit, 16...Cross correlation calculation circuit, 17
...first determination circuit, 41 1 to 41 o ... third storage circuit, 42 ... switching circuit, 43 ... second judgment circuit.
Claims (1)
Ａ／Ｄ変換する回路と、このサンプリングされた
信号より電力スペクトルを求める回路と、この回
路の周波数分解能Δの範囲内にわたりΔ／ｎず
つサンプリング位置が変化されたｎ個の基準電力
スペクトルを記憶する回路と、この基準電力スペ
クトルを順次読出す回路と、この読出された基準
電力スペクトルおよび前記求められた電力スペク
トルより相互相関係数を求める回路と、この求め
られた相関係数より最大値を求める第１の判定回
路と、前記基準電力スペクトルを変えて得られた
ｎ個の最大値群の中からさらに最大値を求める第
２の判定回路とを具備したことを特徴とする相関
処理装置。1. A circuit that A/D converts an input signal using a predetermined sampling pulse signal, a circuit that obtains a power spectrum from this sampled signal, and a sampling position that is changed by Δ/n within the frequency resolution Δ range of this circuit. a circuit for storing n reference power spectra; a circuit for sequentially reading out the reference power spectra; a circuit for calculating a cross-correlation coefficient from the read reference power spectrum and the determined power spectrum; a first determination circuit that determines a maximum value from the correlation coefficient obtained by changing the reference power spectrum; and a second determination circuit that further determines a maximum value from a group of n maximum values obtained by changing the reference power spectrum. A correlation processing device characterized by:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5904282A JPS58176774A (en) | 1982-04-09 | 1982-04-09 | Correlation processor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5904282A JPS58176774A (en) | 1982-04-09 | 1982-04-09 | Correlation processor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS58176774A JPS58176774A (en) | 1983-10-17 |
| JPH0340872B2 true JPH0340872B2 (en) | 1991-06-20 |
Family
ID=13101847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5904282A Granted JPS58176774A (en) | 1982-04-09 | 1982-04-09 | Correlation processor |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS58176774A (en) |
-
1982
- 1982-04-09 JP JP5904282A patent/JPS58176774A/en active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS58176774A (en) | 1983-10-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0447330B2 (en) | ||
| CN1480736A (en) | Method for testing electronic component and its instrument | |
| EP0177557A1 (en) | COUNTING DEVICE AND METHOD FOR SAMPLE FREQUENCIES. | |
| JPH06506057A (en) | Method for determining the transmission characteristics of electrical conductors | |
| EP0072706B1 (en) | Sound signal processing apparatus | |
| IE43171L (en) | Digital device | |
| US4833717A (en) | Voice spectrum analyzing system and method | |
| US3715509A (en) | Method and means for providing resolution level selection in a spectrum analyzer | |
| JPH0340872B2 (en) | ||
| JP2867769B2 (en) | Sound measurement method and device | |
| JPH0340876B2 (en) | ||
| JPH0340875B2 (en) | ||
| JPH0340874B2 (en) | ||
| US6968353B2 (en) | Interpolator and interpolation method when interpolating a digital input signal sequence at interpolation instants prescribed by a control signal | |
| JPH0340873B2 (en) | ||
| JP3139803B2 (en) | Impulse response measurement device | |
| JPH10160507A (en) | Peak detector | |
| JP2000284008A (en) | Frequency measuring method and frequency measuring device | |
| RU2099720C1 (en) | Digital spectrum analyzer | |
| JPH07312551A (en) | A / D conversion method, conversion device thereof, and digital arithmetic processing device thereof | |
| RU2133041C1 (en) | Method determining spectrum of electric signals | |
| JPS5816838B2 (en) | Multi-frequency signal reception method | |
| JP3438376B2 (en) | Periodic signal processing method and apparatus | |
| JP2589864B2 (en) | Signal comparison device | |
| JPH029360B2 (en) |