JPH0296666A - 周波数・位相推定装置 - Google Patents
周波数・位相推定装置Info
- Publication number
- JPH0296666A JPH0296666A JP63247330A JP24733088A JPH0296666A JP H0296666 A JPH0296666 A JP H0296666A JP 63247330 A JP63247330 A JP 63247330A JP 24733088 A JP24733088 A JP 24733088A JP H0296666 A JPH0296666 A JP H0296666A
- Authority
- JP
- Japan
- Prior art keywords
- value
- frequency
- sine wave
- phase
- estimated
- 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
- 238000005070 sampling Methods 0.000 claims description 5
- 230000002265 prevention Effects 0.000 description 9
- 238000001228 spectrum Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function 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
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
本発明は、雑音の重畳した正弦波から周波数と位相を推
定する周波数・位相推定装置に関するものである。
定する周波数・位相推定装置に関するものである。
この種の周波数・位相推定装置としては、昭和５９年１
２月２７日に出願された特願昭５９−２７７３６５号「
周波数・位相推定装置」に開示されたものがある。この
従来の周波数・位相推定装置では、雑音の重畳した信号
において、位相方向に加わっている雑音成分にのみ着目
し、入力された正弦波の初期位相と周波数に対する最適
解を求めていた。それによると、時刻ｉにおける信号Ｓ
、が、Ｓ、！ＡｉｅＸｐ ｊ （θ、＋ ｉ ω）と
表される正弦波において、その位相θ８が、θ、；θ。
２月２７日に出願された特願昭５９−２７７３６５号「
周波数・位相推定装置」に開示されたものがある。この
従来の周波数・位相推定装置では、雑音の重畳した信号
において、位相方向に加わっている雑音成分にのみ着目
し、入力された正弦波の初期位相と周波数に対する最適
解を求めていた。それによると、時刻ｉにおける信号Ｓ
、が、Ｓ、！ＡｉｅＸｐ ｊ （θ、＋ ｉ ω）と
表される正弦波において、その位相θ８が、θ、；θ。
＋ｉω
であり、これに雑音（ｅ！）が重畳し等測的に、θ、＝
θ。＋ｉω＋８ｉ なる位相として観測されるとき、雑音（ｅ、）による影
響をできるだけ受けないように初期位相θ。と周波数ω
を推定する場合に、時刻１からｋまでの観測値を用い、
推定誤差の分散を最小にする方法として最小自乗法を用
いている。最小自乗法を用いる場合には、評価関数Ｊを
、 と定め、その逐次回帰法の解として、 ・ ・ ・（１） ＹＬはＹの転置を表す を時刻ｋにおける最適解として求めていた。
θ。＋ｉω＋８ｉ なる位相として観測されるとき、雑音（ｅ、）による影
響をできるだけ受けないように初期位相θ。と周波数ω
を推定する場合に、時刻１からｋまでの観測値を用い、
推定誤差の分散を最小にする方法として最小自乗法を用
いている。最小自乗法を用いる場合には、評価関数Ｊを
、 と定め、その逐次回帰法の解として、 ・ ・ ・（１） ＹＬはＹの転置を表す を時刻ｋにおける最適解として求めていた。
この従来技術では、サンプリング定理のスペクトラムの
繰り返しに起因して複数のキャリアスペクトルが発生し
、周波数の推定値ωは複数のうちどれか一つを推定して
いた。これは従来技術が信号を周期Ｔでサンプルしてい
るため、１／Ｔ間隔の繰り返しで±１／２Ｔ内の帯域の
スペクトラムが繰り返し発生されるためである。従来技
術では周波数の推定値が、この複数のキャリアスペクト
ルの内どれか一つのスペクトラムを推定していた。
繰り返しに起因して複数のキャリアスペクトルが発生し
、周波数の推定値ωは複数のうちどれか一つを推定して
いた。これは従来技術が信号を周期Ｔでサンプルしてい
るため、１／Ｔ間隔の繰り返しで±１／２Ｔ内の帯域の
スペクトラムが繰り返し発生されるためである。従来技
術では周波数の推定値が、この複数のキャリアスペクト
ルの内どれか一つのスペクトラムを推定していた。
本来、±１／２Ｔ帯域内の信号以外はこの帯域の繰り返
しであるため、たとえ他のスペクトルを推定したとして
も、±１／２Ｔ帯域内の信号を推定したとみなすべきで
、従来技術では±１／２Ｔ帯域外の推定値でもそのまま
の値として捕らえられるという問題点があった。
しであるため、たとえ他のスペクトルを推定したとして
も、±１／２Ｔ帯域内の信号を推定したとみなすべきで
、従来技術では±１／２Ｔ帯域外の推定値でもそのまま
の値として捕らえられるという問題点があった。
本発明の目的は、この問題点を解決した周波数・位相推
定装置を提供することにある。
定装置を提供することにある。
本発明は、正弦波あるいは雑音の重畳した正弦波を入力
とし、入力正弦波の初期位相と周波数を推定する周波数
・位相推定装置において、入力正弦波をある周期Ｔでサ
ンプルし正弦波サンプル値を出力するサンプル手段と、 入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の推定値を発生ずる手
段と、 前記正弦波サンプル値と前記正弦波の推定値との位相誤
差を検出する手段と、 前記位相誤差に第１の時変係数を乗じた値に従って前記
初期位相の推定値を逐次的に更新する初期位相推定手段
と、 前記位相誤差に第２の時変係数を乗じた値により新しい
周波数の推定値を求める周波数推定手段と、 前記新しい周波数の推定値が±π／Ｔの範囲の値にする
ために２πｉ／Ｔ（ｉ＝・・、−ｉ、ｏ。
とし、入力正弦波の初期位相と周波数を推定する周波数
・位相推定装置において、入力正弦波をある周期Ｔでサ
ンプルし正弦波サンプル値を出力するサンプル手段と、 入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の推定値を発生ずる手
段と、 前記正弦波サンプル値と前記正弦波の推定値との位相誤
差を検出する手段と、 前記位相誤差に第１の時変係数を乗じた値に従って前記
初期位相の推定値を逐次的に更新する初期位相推定手段
と、 前記位相誤差に第２の時変係数を乗じた値により新しい
周波数の推定値を求める周波数推定手段と、 前記新しい周波数の推定値が±π／Ｔの範囲の値にする
ために２πｉ／Ｔ（ｉ＝・・、−ｉ、ｏ。
１、・・）を減算し、前記周波数の推定値を更新する手
段とを有することを特徴としている。
段とを有することを特徴としている。
本発明によれば、従来技術で推定されていた周波数の値
δｉに対して以下の操作によりスペクトルの繰り返しを
防ぎ、問題点を解決している。弐（１）によって得られ
た周波数推定値δは真の値δ。
δｉに対して以下の操作によりスペクトルの繰り返しを
防ぎ、問題点を解決している。弐（１）によって得られ
た周波数推定値δは真の値δ。
に対して
δ＝５、±ｎ／Ｔ （ｎは整数）
のいずれかの値をとる。この値δを±１／２Ｔ帯域内の
真値のスペクトルにするため、次式による操作または同
様な操作を行う。
真値のスペクトルにするため、次式による操作または同
様な操作を行う。
’２Ｂ’ ＝ （ｍｏｄ （２；’Ｔ＋ （ｉｎｔ （
２ａｂｓ（δ１Ｔ）／π）２＋３）π／４．π／２）−
π／４）Ｔ ・・・（２）但しａｂ
ｓ （ｘ）；ｘの絶対値 ｉｎｔ （ｘ）；ｘを越えない最大の整数ｍｏｄ （ｘ
、ｚ）＝ｘ−ｉｎｔ （ｘ／ｚ）ｚ〔実施例〕 以下本発明の実施例について図面を参照して詳細に説明
する。
２ａｂｓ（δ１Ｔ）／π）２＋３）π／４．π／２）−
π／４）Ｔ ・・・（２）但しａｂ
ｓ （ｘ）；ｘの絶対値 ｉｎｔ （ｘ）；ｘを越えない最大の整数ｍｏｄ （ｘ
、ｚ）＝ｘ−ｉｎｔ （ｘ／ｚ）ｚ〔実施例〕 以下本発明の実施例について図面を参照して詳細に説明
する。
第１図は本発明の実施例を示すブロック図である。この
周波数・位相推定装置は、入力正弦波をある周期Ｔでサ
ンプルし正弦波サンプル値を出力するサンプル手段を構
成するサンプルホールド回路１及びＡ／Ｄ変換器２と、
入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の推定値を発生する正
弦波発生器４と、正弦波サンプル値と正弦波の推定値と
の位相誤差を検出するディジタル位相検波器３と、推定
のための第１及び第２の時変係数を発生する係数発生器
８と、位相誤差に第１の時変係数を乗じた値に従って初
期位相の推定値を逐次的に更新する初期位相推定器６と
、位相誤差に第２の時変係数を乗じた値により新しい周
波数の推定値を求める周波数推定器７と、新しい周波数
の推定値が±π／Ｔの範囲の値にするために２πｉ／Ｔ
（ｉ−・・、−１，０，１，・・）を減算し周波数の推
定値を更新する繰り返し防止回路９とを有している。
周波数・位相推定装置は、入力正弦波をある周期Ｔでサ
ンプルし正弦波サンプル値を出力するサンプル手段を構
成するサンプルホールド回路１及びＡ／Ｄ変換器２と、
入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の推定値を発生する正
弦波発生器４と、正弦波サンプル値と正弦波の推定値と
の位相誤差を検出するディジタル位相検波器３と、推定
のための第１及び第２の時変係数を発生する係数発生器
８と、位相誤差に第１の時変係数を乗じた値に従って初
期位相の推定値を逐次的に更新する初期位相推定器６と
、位相誤差に第２の時変係数を乗じた値により新しい周
波数の推定値を求める周波数推定器７と、新しい周波数
の推定値が±π／Ｔの範囲の値にするために２πｉ／Ｔ
（ｉ−・・、−１，０，１，・・）を減算し周波数の推
定値を更新する繰り返し防止回路９とを有している。
サンプルホールド回路１は、雑音の重畳された正弦波信
号をサンプルホールドする。
号をサンプルホールドする。
Ａ／Ｄ変換器２は、サンプルホールドされた正弦波をデ
ィジタル化する。
ィジタル化する。
ディジタル位相検波器３は、Ａ／Ｄ変換器２の出力と正
弦波発生器４の発生する正弦波とを比較し、位相誤差を
出力する。この位相検波器の実現手段は、例えば文献（
Ｗ、Ｃ，Ｌｉｎｄｓｅｙ ”Ａ ５ｕｒｖｅｙ ｏｆＤ
ｊｇｉｔａｌ Ｐｈａｓｅ−Ｌｏｃｋｅｄ Ｌｏｏｐ’
Ｐｒｏｃ、 ＩＰＥＥ、 ｖｏｌ。
弦波発生器４の発生する正弦波とを比較し、位相誤差を
出力する。この位相検波器の実現手段は、例えば文献（
Ｗ、Ｃ，Ｌｉｎｄｓｅｙ ”Ａ ５ｕｒｖｅｙ ｏｆＤ
ｊｇｉｔａｌ Ｐｈａｓｅ−Ｌｏｃｋｅｄ Ｌｏｏｐ’
Ｐｒｏｃ、 ＩＰＥＥ、 ｖｏｌ。
６９、 Ｎｏ、 ４． Ａｐｒｉｌ １９８１）に
詳しく示されている。
詳しく示されている。
正弦波発生器４は、初期位相推定器６及び周波数推定器
７と繰り返し防止回路９により推定された初期位相と周
波数とに合致した正弦波を発生するものであり、位相を
入力としその値に対応した正弦波の値を出力とする読み
出し専用メモリ （ＲＯＭ）４４、加算器４３、乗算器
４２、時刻を発生するカウンタ４１により構成すること
ができ、位相推定値と周波数推定値をアドレスとして正
弦波の標本値を出力する。
７と繰り返し防止回路９により推定された初期位相と周
波数とに合致した正弦波を発生するものであり、位相を
入力としその値に対応した正弦波の値を出力とする読み
出し専用メモリ （ＲＯＭ）４４、加算器４３、乗算器
４２、時刻を発生するカウンタ４１により構成すること
ができ、位相推定値と周波数推定値をアドレスとして正
弦波の標本値を出力する。
係数発生器８は、＋１１式の２／にと６／ｋ（ｋ＋１）
を計算し各々の係数を発生する回路であり、逆行列計算
器８１とカウンタ８２とから構成され、カウンタ８２に
よって得られる時刻から逆行列計算器８１によって逆行
列を計算し、初期位相推定と周波数推定のためのそれぞ
れ第１及び第２の時変係数を発生する。係数発生器８は
毎時側逆行列の計算を行うのでマイクロプロセッサ等に
より実現される。
を計算し各々の係数を発生する回路であり、逆行列計算
器８１とカウンタ８２とから構成され、カウンタ８２に
よって得られる時刻から逆行列計算器８１によって逆行
列を計算し、初期位相推定と周波数推定のためのそれぞ
れ第１及び第２の時変係数を発生する。係数発生器８は
毎時側逆行列の計算を行うのでマイクロプロセッサ等に
より実現される。
初期位相推定器６は、レジスタ６１と乗算器６２と加算
器６３とから構成され、位相誤差と第１の時変係数を乗
算器６２にて掛は合わせ、その値を初期位相推定値を保
持するレジスタ６１の内容と加算器６３において加算し
、初期位相推定値の更新値として再びレジスタに格納す
ると共に、正弦波発生２Ｓ４に入力する。
器６３とから構成され、位相誤差と第１の時変係数を乗
算器６２にて掛は合わせ、その値を初期位相推定値を保
持するレジスタ６１の内容と加算器６３において加算し
、初期位相推定値の更新値として再びレジスタに格納す
ると共に、正弦波発生２Ｓ４に入力する。
周波数推定器７は、レジスタ７１と乗算２Ｓ７２と加算
器７３とから構成され、位相誤差と第２の時変係数を乗
算器７２にて掛は合わせ、その値を周波数推定値を保持
するレジスタ７１の内容と加算器７３において減算し、
時刻ｋまでの周波数推定値を繰り返し防止回路９に一時
出力する。
器７３とから構成され、位相誤差と第２の時変係数を乗
算器７２にて掛は合わせ、その値を周波数推定値を保持
するレジスタ７１の内容と加算器７３において減算し、
時刻ｋまでの周波数推定値を繰り返し防止回路９に一時
出力する。
繰り返し防止回路９は、式（２）の演算を行い、その結
果を周波数推定値の更新値として再びレジスタ７１に格
納すると共に、正弦波発生器４に入力する。この繰り返
し防止回路９は、マイクロプロセッサ等により容易に実
現される。
果を周波数推定値の更新値として再びレジスタ７１に格
納すると共に、正弦波発生器４に入力する。この繰り返
し防止回路９は、マイクロプロセッサ等により容易に実
現される。
次に、本実施例の動作を説明する。
雑音の重畳された正弦波信号が周波数・位相推定装置に
入力されると、サンプルホールド回路１でサンプルされ
、Ａ／Ｄ変換器２でディジタル化される。ディジタル化
された正弦波は、ディジタル位相検波器３において、初
期位相推定器６及び周波数推定器７と繰り返し防止回路
９により推定された初期位相と周波数とに合致した正弦
波を発生する正弦波発生器４の発生する正弦波と比較さ
れ、位相誤差が出力される。この位相誤差は、初期位相
推定器６の乗算器６２及び周波数推定器７の乗算器７２
にそれぞれ入力される。
入力されると、サンプルホールド回路１でサンプルされ
、Ａ／Ｄ変換器２でディジタル化される。ディジタル化
された正弦波は、ディジタル位相検波器３において、初
期位相推定器６及び周波数推定器７と繰り返し防止回路
９により推定された初期位相と周波数とに合致した正弦
波を発生する正弦波発生器４の発生する正弦波と比較さ
れ、位相誤差が出力される。この位相誤差は、初期位相
推定器６の乗算器６２及び周波数推定器７の乗算器７２
にそれぞれ入力される。
一方、係数発生器８では、カウンタ８２によって得られ
る時刻から逆行列計算器８１によって逆行列を計算し、
初期位相推定と周波数推定のためのそれぞれ第１及び第
２の時変係数を発生する。初期位相推定のための第１の
時変係数は初期位相推定器６の乗算器６２に入力され、
周波数推定のための第２の時変係数は周波数推定器７の
乗算器７２に入力される。
る時刻から逆行列計算器８１によって逆行列を計算し、
初期位相推定と周波数推定のためのそれぞれ第１及び第
２の時変係数を発生する。初期位相推定のための第１の
時変係数は初期位相推定器６の乗算器６２に入力され、
周波数推定のための第２の時変係数は周波数推定器７の
乗算器７２に入力される。
初期位相推定器６では、位相誤差と第１の時変係数を乗
算器６２にて掛は合わせ、その値を初期位相推定値を保
持するレジスタ６１の内容と加算器６３において加算し
、初期位相推定値の更新値として再びレジスタ６１に格
納する。
算器６２にて掛は合わせ、その値を初期位相推定値を保
持するレジスタ６１の内容と加算器６３において加算し
、初期位相推定値の更新値として再びレジスタ６１に格
納する。
周波数推定器７では、位相誤差と第２の時変係数を乗算
器７２にて掛は合わせ、その値を周波数推定値を保持す
るレジスタ７１の内容と加算器７３において減算し、時
刻ｋまでの周波数推定値を繰り返し防止回路９に一時出
力する。
器７２にて掛は合わせ、その値を周波数推定値を保持す
るレジスタ７１の内容と加算器７３において減算し、時
刻ｋまでの周波数推定値を繰り返し防止回路９に一時出
力する。
繰り返し防止回路９により式（２）の演算を行い、その
結果を周波数推定値の更新値として再びレジスタ７１に
格納する。更新された周波数推定値は正弦波発生器４に
おいてカウンタ４１の内容と乗算器４２で乗算され、更
に初期位相推定値と加算器４３で加算され、時刻ｋにお
ける最適位相推定値として読み出し専用メモリ４４に入
力され正弦波の値が出力される。なお、この読み出し専
用メモリ４４のアドレスの最大値を２πに対応させてお
くことにより、位相推定値はモジュロ２πで行われる。
結果を周波数推定値の更新値として再びレジスタ７１に
格納する。更新された周波数推定値は正弦波発生器４に
おいてカウンタ４１の内容と乗算器４２で乗算され、更
に初期位相推定値と加算器４３で加算され、時刻ｋにお
ける最適位相推定値として読み出し専用メモリ４４に入
力され正弦波の値が出力される。なお、この読み出し専
用メモリ４４のアドレスの最大値を２πに対応させてお
くことにより、位相推定値はモジュロ２πで行われる。
なお、カウンタ４１．８２は本装置の動作開始時点から
毎時刻カウントアツプしており、これらカウンタの値は
等しいため共用することも可能である。
毎時刻カウントアツプしており、これらカウンタの値は
等しいため共用することも可能である。
以上により、入力される雑音の重畳した正弦波から周波
数と位相を逐次的に最小分散推定をすることができる。
数と位相を逐次的に最小分散推定をすることができる。
、推定値はそれぞれレジスタ６１．７１の内容として与
えられる。なお、プログラマブルな信号処理プロセッサ
で上記と等価な動作を容易に実現することができること
は明らかである。
えられる。なお、プログラマブルな信号処理プロセッサ
で上記と等価な動作を容易に実現することができること
は明らかである。
本発明によれば、従来技術では誤って推定していた周波
数イーコ、±ｎ／Ｔをナイキスト帯域内の真値δ、に戻
すことにより、正確に推定できるという効果がある。
数イーコ、±ｎ／Ｔをナイキスト帯域内の真値δ、に戻
すことにより、正確に推定できるという効果がある。
第１図は本発明の一実施例を示すプロ・ツク図である。
ｌ・・・・・サンプルホールド回路
２・・・・・Ａ／Ｄ変換器
３・・・・・ディジタル位相検波器
正弦波発生器
初期位相推定器
周波数推定器
係数発生器
繰り返し防止回路
Claims (1)
- （１）正弦波あるいは雑音の重畳した正弦波を入力とし
、入力正弦波の初期位相と周波数を推定する周波数・位
相装置において、 入力正弦波をある周期Ｔでサンプルし正弦波サンプル値
を出力するサンプル手段と、 入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の推定値を発生する手
段と、 前記正弦波サンプル値と前記正弦波の推定値との位相誤
差を検出する手段と、 前記位相誤差に第１の時変係数を乗じた値に従って前記
初期位相の推定値を逐次的に更新する初期位相推定手段
と、 前記位相誤差に第２の時変係数を乗じた値により新しい
周波数の推定値を求める周波数推定手段と、 前記新しい周波数の推定値が±π／Ｔの範囲の値にする
ために２πｉ／Ｔ（ｉ＝・・、−１、０、１、・・）を
減算し、前記周波数の推定値を更新する手段とを有する
ことを特徴とする周波数・位相推定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63247330A JP2527009B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63247330A JP2527009B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0296666A true JPH0296666A (ja) | 1990-04-09 |
| JP2527009B2 JP2527009B2 (ja) | 1996-08-21 |
Family
ID=17161798
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63247330A Expired - Lifetime JP2527009B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2527009B2 (ja) |
-
1988
- 1988-10-03 JP JP63247330A patent/JP2527009B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2527009B2 (ja) | 1996-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4498184B2 (ja) | 直線性補償回路 | |
| CN107209028B (zh) | 分析器装置 | |
| NO180809B (no) | Fremgangsmåte til korrigering av fasefeil og forsterkningsfeil | |
| JPH0783964A (ja) | 正弦波信号の振幅および周波数を推定する方法、信号の特性の表示を得る方法、および電子回路において複数個のアナログ信号の状態をモニタするためのシステム | |
| US6819279B2 (en) | Method and apparatus for the recovery of signals acquired by an interleaved system of digitizers with mismatching frequency response characteristics | |
| JPH11326407A (ja) | 周波数分析方法及びこの方法を用いた掃引型スペクトラム・アナライザ | |
| US5872538A (en) | Frequency domain correction of I/Q imbalance | |
| JP3610343B2 (ja) | 測定用受信器のスプリアス信号除去方法 | |
| US6794857B2 (en) | Apparatus and method for measuring a phase delay characteristic | |
| JPH0296666A (ja) | 周波数・位相推定装置 | |
| JPH10213613A (ja) | 周波数測定装置 | |
| JP2527008B2 (ja) | 周波数・位相推定装置 | |
| JP2527011B2 (ja) | 周波数・位相推定装置 | |
| Chakraborty et al. | New adaptive algorithm for delay estimation of sinusoidal signals | |
| KR100706217B1 (ko) | 검파 장치, 방법 및 기록 매체 | |
| JPS61154317A (ja) | 周波数・位相推定装置 | |
| JP2624696B2 (ja) | スペクトル推定装置 | |
| JP2005214932A (ja) | 信号処理装置、この信号処理装置を用いた電圧測定装置及び電流測定装置 | |
| JPH1183914A (ja) | 周期検出装置および周期検出方法 | |
| US7705609B2 (en) | Phase frequency distortion measurement system | |
| JP3641469B2 (ja) | 変調搬送波周波数検出方法および装置 | |
| JP3914818B2 (ja) | 回転角度検出装置 | |
| JP2006352390A (ja) | デジタルフィルタ | |
| JP2015004634A (ja) | 回転角度検出装置、画像処理装置及び回転角度検出方法 | |
| JPH1078460A (ja) | 直流成分除去装置および素子定数測定装置 |