JPH0296665A - Apparatus for estimating frequency and phase - Google Patents
Apparatus for estimating frequency and phaseInfo
- Publication number
- JPH0296665A JPH0296665A JP63247324A JP24732488A JPH0296665A JP H0296665 A JPH0296665 A JP H0296665A JP 63247324 A JP63247324 A JP 63247324A JP 24732488 A JP24732488 A JP 24732488A JP H0296665 A JPH0296665 A JP H0296665A
- Authority
- JP
- Japan
- Prior art keywords
- sine wave
- phase
- frequency
- value
- time
- 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
Links
- 238000000034 method Methods 0.000 abstract description 9
- 239000011159 matrix material Substances 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measuring Phase Differences (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、雑音の重畳した正弦波から周波数と位相を推
定する周波数・位相推定装置に関するものである。DETAILED DESCRIPTION OF THE INVENTION [Field of Industrial Application] The present invention relates to a frequency/phase estimating device that estimates frequency and phase from a sine wave on which noise is superimposed.
この種の周波数・位相推定装置としては、昭和５９年１
２月２７日に出願された特願昭５９−２７７３６５号「
周波数・位相推定装置」に開示されたものがある。この
従来の周波数・位相推定装置では、雑音の重畳した信号
において、位相方向に加わっている雑音成分にのみ着目
し、入力された正弦波の初期位相と周波数に対する最適
解を求めていた。それによると、時刻ｉにおける信号Ｓ
、が、Ｓｉ ＝Ａ（ｅｘｐ ｊ （θ。＋ｉω）但し
、Ａ！ ＝Ａ＋ｅｍｉ ｉ ｅａ＆は振幅方向の雑音
と表される正弦波において、その位相θｉが、θ ＝
θ。＋ｉω
であり、これに雑音（ｅｉ）が重畳し等価的に、θ１＝
θ。＋ｉω＋ｅ。As this type of frequency/phase estimation device,
Patent Application No. 59-277365 filed on February 27th “
There is one disclosed in "Frequency/Phase Estimation Device". This conventional frequency/phase estimation device focuses only on the noise component added in the phase direction in a signal with superimposed noise, and finds the optimal solution for the initial phase and frequency of the input sine wave. According to it, the signal S at time i
, but Si = A(exp j (θ.+iω) However, A! = A+emi i ea& is a sine wave represented as noise in the amplitude direction, and its phase θi is θ =
θ. +iω, and noise (ei) is superimposed on this, equivalently, θ1=
θ. +iω+e.
なる位相として観測されるとき、雑音（６ｉ）による影
響をできるだけ受けないように初期位相θ。と周波数ω
を推定する場合に、時刻１からｋまでの観測値を用い、
推定誤差の分散を最小にする方法として最小自乗法を用
いている。最小自乗法を用いる場合には、評価関数Ｊを
、
と定め、その逐次回帰法の解として、
ｙｔはＹの転置を表す
を時刻ｋにおける最適解として求めていた。When observed as a phase, the initial phase θ is set so as to be as unaffected by noise (6i) as possible. and frequency ω
When estimating, using the observed values from time 1 to k,
The method of least squares is used to minimize the variance of estimation errors. When using the least squares method, the evaluation function J is determined as follows, and as a solution of the successive regression method, yt represents the transposition of Y is obtained as the optimal solution at time k.
この従来技術では、振幅方向に含まれている雑音成分に
つ、いては無視し、振幅を一定と考え、位相方向に含ま
れている雑音成分のみにより推定を行っているため、振
幅成分を無視しただけの劣化が生じている。従って、正
確に周波数・位相を推定できないという問題点がある。This conventional technology ignores the noise component included in the amplitude direction, considers the amplitude to be constant, and performs estimation based only on the noise component included in the phase direction, so the amplitude component is ignored. A certain amount of deterioration has occurred. Therefore, there is a problem that the frequency and phase cannot be accurately estimated.
本発明の目的は、この問題点を解決した周波数・位相推
定装置を提供することにある。An object of the present invention is to provide a frequency/phase estimation device that solves this problem.
本発明は、正弦波あるいは雑音の重畳した正弦波を入力
とし、入力正弦波の初期位相と周波数を推定する周波数
・位相推定装置において、入力正弦波の初期位相の推定
値と入力正弦波の周波数の推定値から次時刻における正
弦波の値を発生する手段と、
この正弦波と前記入力正弦波との位相誤差を検出する手
段と、
検出した位相誤差に前記入力正弦波の振幅による重み付
けをする手段と、
前記重み付けを考慮し推定のための第１及び第２の時変
係数を発生する係数発生手段と、前記重み付けされた位
相誤差に第１の時変係数を乗じた値により前記初期位相
の推定値を逐次的に更新する初期位相推定手段と、
前記重み付けされた位相誤差に第２の時変係数を乗じた
値により前記周波数の推定値を逐次的に更新する周波数
推定手段とを′有することを特徴としている。The present invention provides a frequency/phase estimation device that receives a sine wave or a sine wave on which noise is superimposed as input and estimates the initial phase and frequency of the input sine wave. means for generating a value of a sine wave at the next time from an estimated value of , means for detecting a phase error between this sine wave and the input sine wave, and weighting the detected phase error by the amplitude of the input sine wave. means for generating first and second time-varying coefficients for estimation in consideration of the weighting; and determining the initial phase by a value obtained by multiplying the weighted phase error by the first time-varying coefficient. an initial phase estimating means for sequentially updating an estimated value of the frequency; and a frequency estimating means for sequentially updating an estimated value of the frequency by a value obtained by multiplying the weighted phase error by a second time-varying coefficient. It is characterized by having
本発明が対象とするのは、時刻ｉにおける信号ＳＩが、
Ｓ五 ＝Ａ門 ｅｘｐｊ（θ。十 ｉ ω）と表される
正弦波において、その位相θ直が、θ、＝θ。＋ｉω
であり、これに雑音（ｅｉ）が重畳し等価的に、θ、；
θ。＋ｉω＋ｅｊ
なる位相として観測されるとき、初期位相θ。と周波数
ωが未知パラメータであり、雑音（ｅ、）による影響を
できるだけ受けないように観測位相値θ、から初期位相
θ。と周波数ωを推定する問題である。The object of the present invention is that the signal SI at time i is a sine wave expressed as S5 = A gate expj (θ. 10 i ω), and its phase θ is θ, = θ. +iω, and noise (ei) is superimposed on this, equivalently, θ,;
θ. The initial phase θ is observed as +iω+ej. and frequency ω are unknown parameters, and the initial phase θ is calculated from the observed phase value θ, so as to be as unaffected by noise (e, ) as possible. The problem is to estimate the frequency ω.
時刻１からｋまでの観測値を用い、推定誤差の分散を最
小にする方法として、最小自乗法が良く知られている。The least squares method is well known as a method of minimizing the variance of estimation errors using observed values from time 1 to time k.
最小自乗法を上記の問題に当てはめた場合には、評価関
数Ｊを従来技術で述べたように、
と定義するが、ここで振幅成分を考慮し、雑音の重畳し
た観測振幅値Ａ１によって重み付けをし評価関数Ｊを、
と定義する。評価関数Ｊの初期位相θ。及び周波数ωに
関する偏微分係数をそれぞれ０とおくと、以下の方程式
が得られる。When applying the least squares method to the above problem, the evaluation function J is defined as described in the prior art, but here, considering the amplitude component, weighting is performed using the observed amplitude value A1 on which noise is superimposed. The evaluation function J is defined as follows. Initial phase θ of evaluation function J. If the partial differential coefficients with respect to ω and frequency ω are respectively set to 0, the following equation is obtained.
または
初期位相θ。と周波数ωの最小分散推定値は＋１１式の
解として、
と与えられる。この際、最適解を毎時剤逐次的に求める
手段として逐次回帰法が知られている。それによれば時
刻ｉにおける初期位相θ。の推定値を１゜１、周波数ω
の推定値を５１とすると時刻ｋにおいて、
と修正することになる。更に、
ここで（３）式右辺第２項の（）内の値を位相誤差と呼
ぶと、初期位相推定値１゜鯉は位相誤差に振幅Ａ、を乗
じ、更に時刻にと振幅Ａ１に依存する時変係数を乗じ一
時刻前の推定値１゜ト■から減することにより修正し、
また周波数推定価δ１は同じく位相誤差に振幅Ａ、を乗
じ、更に時刻にと振幅Ａ、に依存する時変係数を乗じ一
時刻前の推定値δト１から減することにより修正すると
いう構成をとることにより（３）式が実現される。この
方法によれば、毎時剤その時点までの観測値を用いた最
良の推定が方程式を解く形で与えられるため、極めて短
時間に最適な推定値が得られる。or the initial phase θ. The minimum variance estimate of the frequency ω and the frequency ω is given as the solution to equation +11. At this time, the successive regression method is known as a means for successively finding the optimal solution every hour. According to this, the initial phase θ at time i. The estimated value of 1°1, the frequency ω
If the estimated value of is set to 51, then at time k, it will be corrected as follows. Furthermore, if we call the value in parentheses in the second term on the right side of equation (3) the phase error, then the initial phase estimate of 1° is obtained by multiplying the phase error by the amplitude A, and then depending on the time and the amplitude A1. Corrected by multiplying by the time-varying coefficient and subtracting it from the estimated value one hour ago by 1°,
In addition, the frequency estimation value δ1 is corrected by multiplying the phase error by the amplitude A, further multiplying the time by a time-varying coefficient that depends on the amplitude A, and subtracting it from the estimated value δto1 one time ago. By taking, equation (3) is realized. According to this method, the best estimate using the observed values up to that point in time is given in the form of solving an equation, so the best estimate can be obtained in an extremely short time.
以下本発明の実施例について図面を参照して詳細に説明
する。Embodiments of the present invention will be described in detail below with reference to the drawings.
第１図は本発明の実施例を示すブロック図である。この
周波数・位相推定装置は、入力正弦波の初期位相の推定
値と入力正弦波の周波数の推定値から次時刻における正
弦波の値を発生する正弦波発生器４と、この正弦波発生
器４からの正弦波と入力正弦波との位相誤差を検出する
ディジタル位相検波器３と、検出した位相誤差に入力正
弦波の振幅による重み付けをする重み付は回路５と、重
み付けを考慮し推定のための第１及び第２の時変係数を
発生する係数発生器８と、重み付けされた位相誤差に第
１の時変係数を乗じた値により初期位相の推定値を逐次
的に更新する初期位相推定器６と、重み付けされた位相
誤差に第２の時変係数を乗じた値により周波数の推定値
を逐次的に更新する周波数推定器７とを有している。な
お、１はサンプルホールド回路、２はＡ／Ｄ変換器であ
る。FIG. 1 is a block diagram showing an embodiment of the present invention. This frequency/phase estimation device includes a sine wave generator 4 that generates a sine wave value at the next time from an estimated value of the initial phase of an input sine wave and an estimated value of the frequency of the input sine wave; A digital phase detector 3 detects the phase error between the input sine wave and the input sine wave, and a weighting circuit 5 weights the detected phase error by the amplitude of the input sine wave. a coefficient generator 8 that generates first and second time-varying coefficients, and an initial phase estimation that sequentially updates an estimated value of the initial phase by a value obtained by multiplying the weighted phase error by the first time-varying coefficient. 6, and a frequency estimator 7 that sequentially updates the estimated frequency using a value obtained by multiplying the weighted phase error by a second time-varying coefficient. Note that 1 is a sample and hold circuit, and 2 is an A/D converter.
サンプルホールド回路１は、雑音の重畳された正弦波信
号をサンプルホールドする。A sample and hold circuit 1 samples and holds a sine wave signal on which noise is superimposed.
Ａ／Ｄ変換器２は、サンプルホールドされた正弦波をデ
ィジタル化する。The A/D converter 2 digitizes the sampled and held sine wave.
ディジタル位相検波器３は、Ａ／Ｄ変換器２の出力と正
弦波発生器４の発生する正弦波とを比較し、位相誤差を
出力する。この位相検波器の実現手段は、例えば文献（
１１，Ｃ，Ｌｉｎｄｓｅｙ ’Ａ ５ｕｒｖｅｙ ｏｒ
Ｄｉｇｉｔａｌ Ｐｈａｓｅ−Ｌｏｃｋｅｄ Ｌｏｏｐ
″Ｐｒｏｃ、 ＩＥＥＥ、 ｖｏｌ、６９、 Ｎｏ、
４． Ａｐｒｆｌ １９８１）に詳しく示されてい
る。Digital phase detector 3 compares the output of A/D converter 2 and the sine wave generated by sine wave generator 4, and outputs a phase error. The means for realizing this phase detector can be found, for example, in the literature (
11, C, Lindsey'A 5urvey or
Digital Phase-Locked Loop
"Proc, IEEE, vol. 69, No.
4. Aprfl 1981).
正弦波発生器４は、初期位相推定器６及び周波数推定器
７により推定された初期位相と周波数とに合致した正弦
波を発生するものであり、位相を入力としその値に対応
した正弦波の値を出力とする読み出し専用メモリ （Ｒ
ＯＭ）４４、加算器４３、乗算器４２、時刻を発生する
カウンタ４１により構成することができ、位相推定値と
周波数推定値をアドレスとして正弦波の標本値を出力す
る。The sine wave generator 4 generates a sine wave that matches the initial phase and frequency estimated by the initial phase estimator 6 and frequency estimator 7, and receives the phase as input and generates a sine wave corresponding to the value. Read-only memory (R
OM) 44, an adder 43, a multiplier 42, and a counter 41 that generates time, and outputs a sample value of a sine wave using the phase estimated value and frequency estimated value as addresses.
重み付は回路５は、係数発生器５１と乗算器５２゜５３
とで構成され、信号振幅の大きさにより重み付けの係数
を係数発生器５１で発生させ、前記位相誤差に乗算器５
２．５３で重み付けする回路である。The weighting circuit 5 includes a coefficient generator 51 and multipliers 52 and 53.
A coefficient generator 51 generates a weighting coefficient according to the magnitude of the signal amplitude, and a multiplier 5
This is a circuit that weights by 2.53.
係数発生器８は、（３）式のＧｋ−Ｉｙｋを計算し各々
の係数を発生する回路であり、逆行列計算器８１とカウ
ンタ８２とから構成され、カウンタ８２によって得られ
る時刻と各時刻に得られる重み付は係数とから逆行列計
算器８１によって逆行列を計算し、初期位相推定と周波
数推定のためのそれぞれ第１及び第２の時変係数を発生
する。係数発生器８は毎時側逆行列の計算を行うのでマ
イクロプロセッサ等により実現される。The coefficient generator 8 is a circuit that calculates Gk-Iyk in equation (3) and generates each coefficient, and is composed of an inverse matrix calculator 81 and a counter 82, and is An inverse matrix is calculated from the obtained weighting coefficients by an inverse matrix calculator 81 to generate first and second time-varying coefficients for initial phase estimation and frequency estimation, respectively. The coefficient generator 8 calculates the inverse matrix on the hourly side, so it is realized by a microprocessor or the like.
初期位相推定器６は、レジスタ６１と乗算器６２と加算
器６３とから構成され、重み付けされた位相誤差と第１
の時変係数を乗算器６２にて掛は合わせ、その値を初期
位相推定値を保持するレジスタ６１の内容と加算器６３
において加算し、初期位相推定値の更新値として再びレ
ジスタに格納すると共に、正弦波発生器４に入力する。The initial phase estimator 6 is composed of a register 61, a multiplier 62, and an adder 63, and is configured to calculate the weighted phase error and the first
A multiplier 62 multiplies the time-varying coefficients of
is added and stored in the register again as an updated value of the initial phase estimation value, and is also input to the sine wave generator 4.
周波数推定器７は、レジスタ７１と乗算器７２と加算器
７３とから構成され、重み付けされた位相誤差と第２の
時変係数を乗算器７２にて掛は合わせ、その値を周波数
推定値を保持するレジスタ７１の内容と加算器７３にお
いて加算し、周波数推定値の更新値として再びレジスタ
７１に格納すると共に、正弦波発生器４に入力する。The frequency estimator 7 is composed of a register 71, a multiplier 72, and an adder 73, and multiplies the weighted phase error and a second time-varying coefficient together in the multiplier 72, and uses the resulting value as a frequency estimation value. It is added to the contents of the held register 71 in an adder 73, and stored in the register 71 again as an updated value of the frequency estimate, and is input to the sine wave generator 4.
次に、本実施例の動作を説明する。Next, the operation of this embodiment will be explained.
雑音の重畳された正弦波信号が周波数・位相推定装置に
入力されると、サンプルホールド回路１でサンプルされ
、Ａ／Ｄ変換器２でディジタル化される。ディジタル化
された正弦波は、ディジタル位相検波器３において、初
期位相推定器６及び周波数推定器７により推定された初
期位相と周波数とに合致した正弦波を発生する正弦波発
生器４０発生する正弦波と比較され、位相誤差が出力さ
れる。この位相誤差は、重み付は回路５の乗算器５２、
５３に入力される。When a noise-superimposed sine wave signal is input to the frequency/phase estimation device, it is sampled by a sample and hold circuit 1 and digitized by an A/D converter 2. The digitized sine wave is passed to the digital phase detector 3 by a sine wave generator 40 that generates a sine wave that matches the initial phase and frequency estimated by the initial phase estimator 6 and frequency estimator 7. The phase error is output. This phase error is weighted by the multiplier 52 of the circuit 5;
53.
一方、Ａ／Ｄ変換器２の出力は、重み付は回路５の係数
発生器５１に入力され、係数発生器は信号振幅の大きさ
により重み付けの係数を発生し、この重み付は係数は乗
算器５２．５３に入力される。乗算器５２．５３は、前
記位相誤差に重み付は係数を掛は合わせることによって
、位相誤差に重み付けを行う。乗算器５２の出力は初期
位相推定器６０乗算器６２に入力され、乗算器５３の出
力は周波数推定器７の乗算器７２に入力される。On the other hand, the output of the A/D converter 2 is input to the coefficient generator 51 of the weighting circuit 5, and the coefficient generator generates a weighting coefficient depending on the magnitude of the signal amplitude. 52 and 53. Multipliers 52 and 53 weight the phase errors by multiplying them by weighting coefficients. The output of the multiplier 52 is input to the initial phase estimator 60 multiplier 62, and the output of the multiplier 53 is input to the multiplier 72 of the frequency estimator 7.
一方、重み付は回路５の係数発生器５１の出力である重
み付は係数は、係数発生器８の逆行列計算器８１に入力
される。逆行列計算器８１では、カウンタ８２によって
得られる時刻と各時刻に得られる重み付は係数とから逆
行列を計算し、初期位相推定と周波数推定のためのそれ
ぞれ第１及び第２の時変係数を発生する。初期位相推定
のための第１の時変係数は初期位相推定器６の乗算器６
２に入力され、周波数推定のための第２の時変係数は周
波数推定器７の乗算器７２に入力される。On the other hand, the weighted coefficients output from the coefficient generator 51 of the weighted circuit 5 are input to the inverse matrix calculator 81 of the coefficient generator 8 . The inverse matrix calculator 81 calculates an inverse matrix from the time obtained by the counter 82 and the weighting coefficient obtained at each time, and calculates the first and second time-varying coefficients for initial phase estimation and frequency estimation, respectively. occurs. The first time-varying coefficient for initial phase estimation is the multiplier 6 of the initial phase estimator 6.
2 and the second time-varying coefficient for frequency estimation is input to the multiplier 72 of the frequency estimator 7.
初期位相推定器６では、重み付けされた位相誤差と第１
の時変係数を乗算器６２にて掛は合わせ、その値を初期
位相推定値を保持するレジスタ６１の内容と加算器６３
において加算し、初期位相推定値の更新値として再びレ
ジスタ６１に格納する。The initial phase estimator 6 uses the weighted phase error and the first
A multiplier 62 multiplies the time-varying coefficients of
is added and stored in the register 61 again as an updated value of the initial phase estimate.
周波数推定器７では、重み付けされた位相誤差と第２の
時変係数を乗算器７２にて掛は合わせ、その値を周波数
推定値を保持するレジスタ７１の内容と加算器７３にお
いて加算し、周波数推定値の更新値として再びレジスタ
７１に格納する。更新された周波数推定値は正弦波発生
器４においてカウンタ４１の内容と乗算器４２で乗算さ
れ、更に初期位相推定値と加算器４３で加算され、時刻
ｋにおける最適位相推定値として読み出し専用メモリ４
４に入力され正弦波の値が出力される。なお、この読み
出し専用メモリ４４のアドレスの最大値を２πに対応さ
せておくことにより、位相推定値はモジュロ２πで行わ
れる。In the frequency estimator 7, a multiplier 72 multiplies the weighted phase error and the second time-varying coefficient, and an adder 73 adds the resulting value to the contents of a register 71 that holds the frequency estimate. The estimated value is stored in the register 71 again as an updated value. The updated frequency estimate is multiplied by the content of the counter 41 in the sine wave generator 4 in the multiplier 42, and further added to the initial phase estimate in the adder 43, and stored in the read-only memory 4 as the optimal phase estimate at time k.
4, and the value of the sine wave is output. Note that by making the maximum value of the address of this read-only memory 44 correspond to 2π, the phase estimation value is performed modulo 2π.
なお、カウンタ４１．８２は本装置の動作開始時点から
毎時側カウントアツプしており、これらカウンタの値は
等しいため共用することも可能である。Note that the counters 41 and 82 have been incrementing every hour since the start of the operation of this device, and since the values of these counters are equal, they can also be used in common.
以上により、入力される雑音の重畳した正弦波から周波
数と位相を逐次的に最小分散推定をすることができる。As described above, it is possible to successively estimate the minimum variance of the frequency and phase from the input sine wave on which noise is superimposed.
推定値はそれぞれレジスタ６１．７１の内容として与え
られる。なお、プログラマブルな信号処理プロセッサで
上記と等価な動作を容易に実現することができることは
明らかである。The estimated values are given as the contents of registers 61, 71, respectively. Note that it is clear that operations equivalent to the above can be easily realized using a programmable signal processing processor.
本発明によれば、従来技術よりも特性良くしかも正確に
周波数・位相を推定できるという効果がある。According to the present invention, there is an effect that the frequency and phase can be estimated more accurately and with better characteristics than the conventional technology.
第１図は本発明の一実施例を示すブロック図である。 １・・・・・サンプルホールド回路 ２・・・・・Ａ／Ｄ変換器 ３・・・・・ディジタル位相検波器 ４・・・・・正弦波発生器 ５・・・・・重み付は回路 ６・・・・・初期位相推定器 ７・・・・・周波数推定器 ８・・・・・係数発生器 FIG. 1 is a block diagram showing one embodiment of the present invention. 1...Sample hold circuit 2...A/D converter 3...Digital phase detector 4...Sine wave generator 5...Weighting is a circuit 6...Initial phase estimator 7...Frequency estimator 8...Coefficient generator
Claims (1)
、入力正弦波の初期位相と周波数を推定する周波数・位
相推定装置において、 入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の値を発生する手段と
、 この正弦波と前記入力正弦波との位相誤差を検出する手
段と、 検出した位相誤差に前記入力正弦波の振幅による重み付
けをする手段と、 前記重み付けを考慮し推定のための第１及び第２の時変
係数を発生する係数発生手段と、前記重み付けされた位
相誤差に第１の時変係数を乗じた値により前記初期位相
の推定値を逐次的に更新する初期位相推定手段と、 前記重み付けされた位相誤差に第２の時変係数を乗じた
値により前記周波数の推定値を逐次的に更新する周波数
推定手段とを有することを特徴とする周波数・位相推定
装置。(1) In a frequency/phase estimation device that receives a sine wave or a sine wave with superimposed noise as input and estimates the initial phase and frequency of the input sine wave, the estimated value of the initial phase of the input sine wave and the frequency of the input sine wave are calculated. means for generating a value of a sine wave at the next time from an estimated value; means for detecting a phase error between this sine wave and the input sine wave; and means for weighting the detected phase error by the amplitude of the input sine wave. and coefficient generation means for generating first and second time-varying coefficients for estimation in consideration of the weighting, and calculating the initial phase by a value obtained by multiplying the weighted phase error by the first time-varying coefficient. comprising: initial phase estimating means for sequentially updating the estimated value; and frequency estimating means for sequentially updating the estimated value of the frequency by a value obtained by multiplying the weighted phase error by a second time-varying coefficient. A frequency/phase estimation device characterized by:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63247324A JP2527008B2 (en) | 1988-10-03 | 1988-10-03 | Frequency / phase estimation device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63247324A JP2527008B2 (en) | 1988-10-03 | 1988-10-03 | Frequency / phase estimation device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0296665A true JPH0296665A (en) | 1990-04-09 |
| JP2527008B2 JP2527008B2 (en) | 1996-08-21 |
Family
ID=17161704
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63247324A Expired - Lifetime JP2527008B2 (en) | 1988-10-03 | 1988-10-03 | Frequency / phase estimation device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2527008B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008136443A1 (en) * | 2007-04-26 | 2008-11-13 | The University Of Tokyo | Sine wave parameter estimation method |
| CN102879642A (en) * | 2012-10-31 | 2013-01-16 | 南京大学 | Frequency estimation method for sine signal |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61154317A (en) * | 1984-12-27 | 1986-07-14 | Nec Corp | Frequency and phase estimating device |
| JPS63252014A (en) * | 1987-04-08 | 1988-10-19 | Kokusai Denshin Denwa Co Ltd <Kdd> | Phase locked loop system |
-
1988
- 1988-10-03 JP JP63247324A patent/JP2527008B2/en not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61154317A (en) * | 1984-12-27 | 1986-07-14 | Nec Corp | Frequency and phase estimating device |
| JPS63252014A (en) * | 1987-04-08 | 1988-10-19 | Kokusai Denshin Denwa Co Ltd <Kdd> | Phase locked loop system |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008136443A1 (en) * | 2007-04-26 | 2008-11-13 | The University Of Tokyo | Sine wave parameter estimation method |
| JP5553334B2 (en) * | 2007-04-26 | 2014-07-16 | 国立大学法人 東京大学 | Sinusoidal parameter estimation method |
| CN102879642A (en) * | 2012-10-31 | 2013-01-16 | 南京大学 | Frequency estimation method for sine signal |
| CN102879642B (en) * | 2012-10-31 | 2014-12-17 | 南京大学 | Frequency estimation method for sine signal |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2527008B2 (en) | 1996-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20040085096A1 (en) | Efficient digital method of and system for determining the instantaneous phase and amplitude of a vibratory accelerometer and other sensors | |
| JPH0829229A (en) | Coriolis mass flowmeter | |
| JPH0783964A (en) | Method for estimating amplitude and frequency of sine-wave signal, method for obtaining display of signal characteristic and system for monitoring plurality of analog-signal states in electronic circuit | |
| Hayes et al. | Sinusoidal frequency estimation by gradient descent | |
| JP2002252991A (en) | Motor control apparatus | |
| JP2009038885A (en) | Signal extracting device and reactive power compensator containing the same | |
| JPH0290082A (en) | Measuring device for detecting distance "h" up to body | |
| WO2002103377A1 (en) | Jitter measuring apparatus and jitter measuring method | |
| JPH0296665A (en) | Apparatus for estimating frequency and phase | |
| JP3887192B2 (en) | Independent component analysis method and apparatus, independent component analysis program, and recording medium recording the program | |
| JPH10213613A (en) | Frequency measuring apparatus | |
| JP2527011B2 (en) | Frequency / phase estimation device | |
| JPS61154317A (en) | Frequency and phase estimating device | |
| JPH0296666A (en) | Apparatus for estimating frequency and phase | |
| JP4344356B2 (en) | Detector, method, program, recording medium | |
| JP2003047299A (en) | Method of controlling vector of pw motor | |
| JPH1183914A (en) | Period detecting equipment and period detecting method | |
| JPS5876891A (en) | Voice pitch extraction | |
| JPH0415404B2 (en) | ||
| JPH0770955B2 (en) | Filter coefficient calculation setting method | |
| Petrović | Computational Effective Modified Newton–Raphson Algorithm for Power Harmonics Parameters Estimation | |
| JP3641469B2 (en) | Modulated carrier frequency detection method and apparatus | |
| JP4320949B2 (en) | Signal incident angle estimation method | |
| JP3239162B2 (en) | Ghost measurement method by waveform matching | |
| Wilkes et al. | Spectral line tracking for nonstationary random processes |