Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH0340875B2 - - Google Patents
[go: Go Back, main page]

JPH0340875B2 - - Google Patents

Info

Publication number
JPH0340875B2
JPH0340875B2 JP5904582A JP5904582A JPH0340875B2 JP H0340875 B2 JPH0340875 B2 JP H0340875B2 JP 5904582 A JP5904582 A JP 5904582A JP 5904582 A JP5904582 A JP 5904582A JP H0340875 B2 JPH0340875 B2 JP H0340875B2
Authority
JP
Japan
Prior art keywords
circuit
frequency
power spectrum
signal
oscillation
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
Application number
JP5904582A
Other languages
Japanese (ja)
Other versions
JPS58176777A (en
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 filed Critical
Priority to JP5904582A priority Critical patent/JPS58176777A/en
Publication of JPS58176777A publication Critical patent/JPS58176777A/en
Publication of JPH0340875B2 publication Critical patent/JPH0340875B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation 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] This invention relates to a correlation process that detects a signal having a spectrum equivalent to a reference power spectrum from an input signal by calculating a correlation coefficient of power spectra. Regarding equipment.

〔発明の技術的背景とその問題点〕[Technical background of the invention and its problems]

第1図は従来の相関処理装置を示すものであ
る。入力信号iはアナログ/デイジタル(A/
D)変換回路11に供給され、この変換回路11
において所定周期のサンプリングパルス信号sに
より、デイジタル時系列信号に変換される。この
サンプルデータは第1の記憶回路12に供給され
記憶される。この記憶されるサンプル数は2N個
とする。この記憶回路12より読出されたサンプ
ルデータは順次電力スペクトル計算回路13に供
給され、例えばFFT(Fast Fourier Transform)
方式によつてN個の電力スペクトルデータxi(i
=0、1、…、N−1)が求められる。この電力
スペクトルデータxiは第2の記憶回路14に供給
され記憶される。
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).
Depending on the method, N pieces of power spectrum data xi(i
=0, 1,...,N-1) are obtained. This power spectrum data xi is supplied to the second storage circuit 14 and stored therein.

一方、第3の記憶回路15には予じめ検出した
い信号に対応する基準の電力スペクトルデータyi
(i=0、1、…、M−1;M<N)が記憶され
ている。この記憶回路15および14に記憶され
たデータはそれぞれ順次読出され、相互相関計算
回路16に供給される。この計算回路16では次
の計算が行なわれる。
On the other hand, the third storage circuit 15 stores reference power spectrum data yi 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=0、1、…、N−1 =1/Mk+M-1i=k xi =1/MM=1i=0 yi この計算によつて得られた相関係数γ(k)は判定
回路17に供給され、この判定回路17において
相関係数γ(k)の最大値γ(k)maxが判定される。
そして、このkの値により基準のスペクトルパタ
ーンとその電力スペクトルパターンが一致する信
号の中心周波数が求められる。第2図a,b,c
はそれぞれxi、yi、γ(k)の一例を示すものであ
る。
However, k=0, 1,..., N-1 = 1/M k+M-1i=k xi = 1/M M=1i=0 yi Correlation coefficient obtained by this calculation γ(k) is supplied to a determination circuit 17, which determines the maximum value γ(k)max of the correlation coefficient γ(k).
Then, from this value of k, the center frequency of the signal whose power spectrum pattern matches the reference spectrum pattern is determined. Figure 2 a, b, c
are examples of xi, yi, and γ(k), respectively.

ところで、入力信号の電力スペクトルxiは周波
数分解能Δ(但し、Δ=s/2N)毎の離散的なデ
ー タであるため、第3図a,bに示す如く信号の中
心周波数0の値によつて本来は同一スペクトルの
筈が、かなり異なつたスペクトルデータとなる可
能性を有している。尚、第3図a,b中点線は連
続スペクトルエンベロープを示し、 第3図aは0=(m+1/2)×Δの場合 第3図bは0=m×Δの場合 (但し、mは正整数) であり、0の値によつてこの両者を両極端として
その中間のスペクトルパターンをとることにな
る。
By the way, since the power spectrum xi of the input signal is discrete data for each frequency resolution Δ (where Δ=s/2N), it is determined by the value of the center frequency 0 of the signal as shown in Figure 3 a and b. 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.

このことは、即ち基準のスペクトルパターンが
例えば第3図aの形で記憶されているとすると、
第3図に示すaとbの中間のパターンの場合には
それだけ相関係数の最大値r(k)maxが変動する
ことになる。このように、従来の装置は信号の中
心周波数と周波数分解能との相対関係によつて相
関係数が大きく変動するという欠点を有してい
る。
This means that if the reference spectral pattern is stored, for example, in the form shown in Figure 3a,
In the case of a pattern between a and b shown in FIG. 3, the maximum value r(k)max 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.

〔発明の目的〕 この発明は上記事情に基づいてなされたもの
で、その目的とするところは入力信号の中心周波
数が如何なる場合であつても所定の相関係数を得
ることが可能な相関処理装置を提供しようとする
ものである。
[Object of the Invention] This invention was made based on the above circumstances, and its object is to provide a correlation processing device that can obtain a predetermined correlation coefficient no matter what the center frequency of the input signal. This is what we are trying to provide.

〔発明の概要〕[Summary of the invention]

この発明は入力信号の電力スペクトルを求めこ
のスペクトルと基準スペクトルの相関係数のうち
最大のものを求めるとともに、前記入力信号を電
力スペクトル計算回路によつて決る周波数分解能
Δの範囲内でn段階変化させ、この変化したそ
れぞれの信号に対して求めたn個の前記相関係数
の最大値群より最も大きなものを求めることによ
り常に変動の少ない相関係数を得るものである。
This invention obtains the power spectrum of an input signal, finds the maximum correlation coefficient between this spectrum and a reference spectrum, and changes the input signal in n steps within a frequency resolution Δ determined by a power spectrum calculation circuit. By determining the largest value from a group of n maximum values of the correlation coefficients determined for each of the changed signals, a correlation coefficient with little variation is always obtained.

〔発明の実施例〕[Embodiments of the invention]

以下、この発明の一実施例について図面を参照
して説明する。尚、第4図中第1図と同一部分に
は同一符号を付す。
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.

第4図において、入力信号iは周波数変換回路
41に供給される。この変換回路41には局部発
振回路42の出力信号が供給され、この信号によ
り前記入力信号iの周波数が変化される。即ち、
局部発振回路42は出力周波数が所定の範囲変化
可能とされており、この信号により入力信号iの
周波数は前記周波数分解能Δの範囲内を偏移さ
れるようになつている。この周波数変換された信
号はそれぞれA/D変換回路11〜相互相関計算
回路16、第1の判定回路17に順に供給され、
第1図と同様の動作により第1の判定回路17か
ら相関係数の最大値群{γn(k(o)}が求められ
る。
In FIG. 4, input signal i is supplied to frequency conversion circuit 41. In FIG. An output signal from a local oscillation circuit 42 is supplied to this conversion circuit 41, and the frequency of the input signal i is changed by this signal. That is,
The local oscillation circuit 42 is capable of changing its output frequency within a predetermined range, and this signal causes the frequency of the input signal i to be shifted within the range of the frequency resolution Δ. The frequency-converted signals are sequentially supplied to the A/D conversion circuit 11 to the cross-correlation calculation circuit 16 and the first determination circuit 17,
The maximum value group {γn(k( o )}) of the correlation coefficients is obtained from the first determination circuit 17 by the same operation as in FIG.

さて、この装置では先ず最初に周波数変換回路
41の出力周波数が入力信号iの周波数と等しく
なるよう局部発振回路42の周波数が選択され、
このときの相関係数の最大値γ1(k1)が計算され
る。次に、局部発振回路42の周波数がΔ/n
だけ変化され、同様の計算によりこの場合の最大
値γ2(k2)が求められる。以下、局部発振回路4
2の周波数がΔ/nずつ変化されて順次相関係
数の最大値が計算され、γ1(k1)、γ2(k2)、…、γ
o
(kn)のn個の最大値群が求められる。この結
果、入力信号iの周波数はΔ/nずつ幅Δだけ
変化されたことになるから、第3図aから同図b
までのスペクトルパターンのn通りのパターンを
とることになる。したがつて、n個の最大値群の
中から最も大なる最大値γ(k)maxを選択すれば、
基準の周波数パターンに最も近いときの相関係数
を得ることができる。第2の判定回路43はこれ
を行なうものであり、第1の判定回路17で得ら
れたn個の最大値群より最も大なる最大値γ(k)
maxが判定される。
Now, in this device, first, the frequency of the local oscillation circuit 42 is selected so that the output frequency of the frequency conversion circuit 41 is equal to the frequency of the input signal i,
The maximum value γ 1 (k 1 ) of the correlation coefficient at this time is calculated. Next, the frequency of the local oscillation circuit 42 is Δ/n
The maximum value γ 2 (k 2 ) in this case is determined by a similar calculation. Below, local oscillation circuit 4
The maximum value of the correlation coefficient is calculated sequentially by changing the frequency of 2 by Δ/n, and γ 1 (k 1 ), γ 2 (k 2 ), ..., γ
o
A group of n maximum values of (kn) is determined. As a result, the frequency of the input signal i is changed by the width Δ by Δ/n, so from Figure 3a to Figure 3b
There are n different spectral patterns up to. Therefore, if we select the largest maximum value γ(k)max from a group of n maximum values, we get
The correlation coefficient closest to the reference frequency pattern can be obtained. The second determination circuit 43 performs this, and determines the maximum value γ(k) that is the largest among the n maximum value groups obtained by the first determination circuit 17.
max is determined.

このような構成とすることにより入力信号iの
中心周波数0が周波数分解能Δの中のどの位置
にあつても常に変動の少ない相関係数γ(k)max
を得ることが可能となる。
With this configuration, no matter where the center frequency 0 of input signal i is within the frequency resolution Δ, the correlation coefficient γ(k)max always has little variation.
It becomes possible to obtain.

〔発明の効果〕〔Effect of the invention〕

以上、詳述したようにこの発明によれば、入力
信号の中心周波数が如何なる場合であつても所定
の相関係数を得ることが可能な相関処理装置を提
供できる。
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.

【図面の簡単な説明】[Brief explanation of drawings]

第1図は従来の相関処理装置の一例を示す構成
図、第2図a,b,cおよび第3図a,bはそれ
ぞれ従来およびこの発明の動作を説明するために
示す図、第4図はこの発明に係わる相関処理装置
の一実施例を示す構成図である。 11……A/D変換回路、13……電力スペク
トル計算回路、16……相互相関計算回路、17
……第1の判定回路、41……周波数変換回路、
42……局部発振回路、43……第2の判定回
路。
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...frequency conversion circuit,
42... Local oscillation circuit, 43... Second determination circuit.

Claims (1)

【特許請求の範囲】[Claims] 1 発振周波数を変化させることが可能な発振回
路と、この発振信号により入力信号の周波数を変
化させる周波数変換回路と、この変換出力信号を
所定のサンプリングパルス信号でA/D変換する
回路と、このサンプリングされた信号より電力ス
ペクトルを求める回路と、この求められた電力ス
ペクトルおよび基準電力スペクトルが供給され相
互相関係数を求める回路と、この求められた相関
係数より最大値を求める第1の判定回路と、前記
電力スペクトルを求める回路の周波数分解能の範
囲内において前記発振回路の発振周波数をn段階
変化させたたときこの第1の判定回路より得られ
るn個の最大値群の中から最大値を求める第2の
判定回路とを具備したことを特徴とする相関処理
装置。
1. An oscillation circuit that can change the oscillation frequency, a frequency conversion circuit that changes the frequency of an input signal using this oscillation signal, a circuit that A/D converts this conversion output signal with a predetermined sampling pulse signal, and A circuit that obtains a power spectrum from a sampled signal, a circuit that is supplied with the obtained power spectrum and a reference power spectrum and obtains a cross-correlation coefficient, and a first judgment that obtains a maximum value from the obtained correlation coefficient. the maximum value from among the group of n maximum values obtained by the first judgment circuit when the oscillation frequency of the oscillation circuit is changed by n steps within the frequency resolution range of the circuit and the circuit for calculating the power spectrum. A correlation processing device characterized by comprising: a second determination circuit for determining .
JP5904582A 1982-04-09 1982-04-09 Correlation processor Granted JPS58176777A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP5904582A JPS58176777A (en) 1982-04-09 1982-04-09 Correlation processor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP5904582A JPS58176777A (en) 1982-04-09 1982-04-09 Correlation processor

Publications (2)

Publication Number Publication Date
JPS58176777A JPS58176777A (en) 1983-10-17
JPH0340875B2 true JPH0340875B2 (en) 1991-06-20

Family

ID=13101942

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5904582A Granted JPS58176777A (en) 1982-04-09 1982-04-09 Correlation processor

Country Status (1)

Country Link
JP (1) JPS58176777A (en)

Also Published As

Publication number Publication date
JPS58176777A (en) 1983-10-17

Similar Documents

Publication Publication Date Title
WO2018188228A1 (en) High-precision frequency measuring system and method
JPH0447330B2 (en)
EP0177557B1 (en) Counting apparatus and method for frequency sampling
JP2004528742A (en) High-speed filter
JPH095160A (en) Fourier spectrometer
JP3918034B2 (en) Method and apparatus for determining mask limits
US3715509A (en) Method and means for providing resolution level selection in a spectrum analyzer
JP2867769B2 (en) Sound measurement method and device
JPH0340875B2 (en)
US4351032A (en) Frequency sensing circuit
JPH10160507A (en) Peak detector
US3717812A (en) Real time analysis of waves
JPH0340873B2 (en)
JPH0340874B2 (en)
JPH0340872B2 (en)
JPH0340876B2 (en)
RU2099720C1 (en) Digital spectrum analyzer
JP3171026B2 (en) Frequency spectrum analyzer
JPS5816838B2 (en) Multi-frequency signal reception method
JP2589864B2 (en) Signal comparison device
JP2556369B2 (en) Overshoot position detection method for analog waveforms
GB2037523A (en) Frequency Sensing Circuit
SU734747A1 (en) Method of logarithmic conversion of voltage into binary code
SU1649566A1 (en) Device of calculation of signal spectrum
JPS6145909B2 (en)