JP2527011B2 - 周波数・位相推定装置 - Google Patents
周波数・位相推定装置Info
- Publication number
- JP2527011B2 JP2527011B2 JP63250554A JP25055488A JP2527011B2 JP 2527011 B2 JP2527011 B2 JP 2527011B2 JP 63250554 A JP63250554 A JP 63250554A JP 25055488 A JP25055488 A JP 25055488A JP 2527011 B2 JP2527011 B2 JP 2527011B2
- Authority
- JP
- Japan
- Prior art keywords
- phase
- sine wave
- 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.)
- Expired - Lifetime
Links
Landscapes
- Measuring Phase Differences (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Description
【発明の詳細な説明】 （産業上の利用分野） 本発明は、雑音の重畳した正弦波から周波数と位相を
推定する周波数位相推定装置に関するものである。
推定する周波数位相推定装置に関するものである。
（従来の技術） 従来技術として知られているものは、Simon Haykin
著，“Introduction to Ad aptive Filter"（Macmillan
publishing Company,New York,1984年刊行）である。
雑音の重畳した信号に於いて、位相方向に加わっている
雑音成分にのみ着目し入力された正弦波の初期位相と周
波数に対する最適解を求める問題は、次の問題に帰着す
る。すなわち、時刻ｉに於ける信号Siが Si＝Aiexp｛ｊ（θ０＋ｉω）｝ 但し、Ai＝Ａ＋eai:eaiは振幅方向の雑音 と表される正弦波に於てその位相θｉ＝θ０＋ｉωがこ
れに雑音｛ei｝が重畳し等価的に θｉ＝θ０＋ｉω＋ei なる位相として観測されるとき、｛ei｝による影響をで
きるだけ受けないようにθ０とωを推定する問題であ
る。時刻１からｋまでの観測値を用い、推定誤差の分散
を最小にする方法として最小自乗法が良く知られてい
る。上記文献によれば最小自乗法を上記の問題に当ては
めた場合には評価関数を と定めていた。更に文献によれば重み付け関数wiとして
は wi＝rk-1（ｒ＜１） がよく知られていると述べている。上式に対して逐次回
帰法による周波数・位相推定（特願昭59−277365号明細
書「周波数・位相推定装置」昭59年12月27日出願）を用
いると、その解として YtはＹの転置を表わす を時刻ｋに於ける最適解として求めることができた。
著，“Introduction to Ad aptive Filter"（Macmillan
publishing Company,New York,1984年刊行）である。
雑音の重畳した信号に於いて、位相方向に加わっている
雑音成分にのみ着目し入力された正弦波の初期位相と周
波数に対する最適解を求める問題は、次の問題に帰着す
る。すなわち、時刻ｉに於ける信号Siが Si＝Aiexp｛ｊ（θ０＋ｉω）｝ 但し、Ai＝Ａ＋eai:eaiは振幅方向の雑音 と表される正弦波に於てその位相θｉ＝θ０＋ｉωがこ
れに雑音｛ei｝が重畳し等価的に θｉ＝θ０＋ｉω＋ei なる位相として観測されるとき、｛ei｝による影響をで
きるだけ受けないようにθ０とωを推定する問題であ
る。時刻１からｋまでの観測値を用い、推定誤差の分散
を最小にする方法として最小自乗法が良く知られてい
る。上記文献によれば最小自乗法を上記の問題に当ては
めた場合には評価関数を と定めていた。更に文献によれば重み付け関数wiとして
は wi＝rk-1（ｒ＜１） がよく知られていると述べている。上式に対して逐次回
帰法による周波数・位相推定（特願昭59−277365号明細
書「周波数・位相推定装置」昭59年12月27日出願）を用
いると、その解として YtはＹの転置を表わす を時刻ｋに於ける最適解として求めることができた。
（発明が解決しようとする問題点） この従来技術では、初期推定時に於て雑音により誤っ
たベクトルで方向に推定を始めるとその後２Пの不確定
性により最終的な推定値が誤ってしまう。この誤推定は
問題であり、基本的には時間的に古い情報を「忘れる」
重み付けをすれば良いことは明白であるが、従来技術の
ような重み付けでは特性の劣化を招いていた。これは上
記文献では指数重み付けを用いているが、指数重み付け
は非定常環境下でデータの統計的変動に追従するもので
あり、最近のデータに対しても情報が指数的に失われる
ためにこの問題には適していない。本発明はこれを時間
関数による固定的な重み付けを用いることにより解決し
ようとしている。
たベクトルで方向に推定を始めるとその後２Пの不確定
性により最終的な推定値が誤ってしまう。この誤推定は
問題であり、基本的には時間的に古い情報を「忘れる」
重み付けをすれば良いことは明白であるが、従来技術の
ような重み付けでは特性の劣化を招いていた。これは上
記文献では指数重み付けを用いているが、指数重み付け
は非定常環境下でデータの統計的変動に追従するもので
あり、最近のデータに対しても情報が指数的に失われる
ためにこの問題には適していない。本発明はこれを時間
関数による固定的な重み付けを用いることにより解決し
ようとしている。
（問題を解決するための手段） 本発明の周波数・位相推定装置は、正弦波あるいは雑
音の重畳した正弦波を入力とし、入力正弦波の初期位相
と周波数を推定する装置において、 入力正弦波の初期位相の推定値と入力正弦波の周波数
の推定値から次時刻における正弦波の値を発生する手段
と、該正弦波と前期入力正弦波との位相誤差を検出する
手段と、推定開始から一定時間だけ重み付け値を発生す
る手段と、検出された位相誤差に前記重み付け値を乗算
する手段と、前記重み付け値を入力し第一及び第二の時
変係数を発生する係数発生手段と、前記重み付けされた
位相誤差に第一の時変係数を乗じた値により前記初期位
相の推定値を逐次的に更新する初期位相推定手段と、前
記重み付けされた位相誤差に第２の時変係数を乗じた値
により前記周波数の推定値を逐次的に更新する周波数推
定手段とを有することを特徴とする。
音の重畳した正弦波を入力とし、入力正弦波の初期位相
と周波数を推定する装置において、 入力正弦波の初期位相の推定値と入力正弦波の周波数
の推定値から次時刻における正弦波の値を発生する手段
と、該正弦波と前期入力正弦波との位相誤差を検出する
手段と、推定開始から一定時間だけ重み付け値を発生す
る手段と、検出された位相誤差に前記重み付け値を乗算
する手段と、前記重み付け値を入力し第一及び第二の時
変係数を発生する係数発生手段と、前記重み付けされた
位相誤差に第一の時変係数を乗じた値により前記初期位
相の推定値を逐次的に更新する初期位相推定手段と、前
記重み付けされた位相誤差に第２の時変係数を乗じた値
により前記周波数の推定値を逐次的に更新する周波数推
定手段とを有することを特徴とする。
（作用） 本発明が対象とするのは、時刻ｉに於ける信号Siが Si＝Aiexp｛ｊ（θ０＋ｉω）｝ と表される正弦波に於てその位相θｉ＝θ０＋ｉωがこ
れに雑音｛ei｝が重畳し等価的に θｉ＝θ０＋ｉω＋ei なる位相として観測されるとき、θ０とωが未知パラメ
ータであり、｛ei｝による影響をできるだけ受けないよ
うに観測値θｉからθ０とωを推定する問題である。時
刻１からｋまでの観測値を用い、推定誤差の分散を最小
にする方法として最小自乗法が良く知られている。最小
自乗法を上記の問題に当てはめた場合には評価関数を従
来技術で述べたように と定義するが、ここで重み関数wiをwi＝rk-1の様な関数
ではなく、最近のデータに対して均等になるような関
数、例えば wi＝１−ri や 但しｓは立ち上がり関数 で定義する。評価関数の に関する偏微分係数をそれぞれ０とおくと以下の方程式
が得られる。
れに雑音｛ei｝が重畳し等価的に θｉ＝θ０＋ｉω＋ei なる位相として観測されるとき、θ０とωが未知パラメ
ータであり、｛ei｝による影響をできるだけ受けないよ
うに観測値θｉからθ０とωを推定する問題である。時
刻１からｋまでの観測値を用い、推定誤差の分散を最小
にする方法として最小自乗法が良く知られている。最小
自乗法を上記の問題に当てはめた場合には評価関数を従
来技術で述べたように と定義するが、ここで重み関数wiをwi＝rk-1の様な関数
ではなく、最近のデータに対して均等になるような関
数、例えば wi＝１−ri や 但しｓは立ち上がり関数 で定義する。評価関数の に関する偏微分係数をそれぞれ０とおくと以下の方程式
が得られる。
又は θ０とωの最小分散推定値は（１）式の解として と与えられる。この際、最適解を毎時刻逐次的に求める
手段として逐次回帰法が知られている。それによれば時
刻ｉに於ける の推定値を とすると時刻ｋにおいて と修正することになる。更に、 ここで（３）式の２行１列の行列Pk -1ykで１行１列の
値が第一の時変係数で２行１列の値が第二の時変係数と
なる。また（３）式右辺第２項の（・）の値を位相誤差
と呼ぶと、初期位相推定値 は位相誤差に重み関数wiを乗じ、更に時刻ｋと重み関数
wiに依存する第一の時変係数を乗じ一時刻前の推定値 から減ずることにより修正し、又周波数推定値 は同じく位相誤差に重み関数wiを乗じ、更に時刻ｋと重
み関数wiに依存する第二の時変係数を乗じ一時刻前の推
定値 から減ずることにより修正をするという構成をとること
により（３）式が実現される。この方法によれば、毎時
刻その時点までの観測値を用いた最良の推定が方程式を
解く形で与えられるため、極めて短時間に最適な推定値
が得られる。
手段として逐次回帰法が知られている。それによれば時
刻ｉに於ける の推定値を とすると時刻ｋにおいて と修正することになる。更に、 ここで（３）式の２行１列の行列Pk -1ykで１行１列の
値が第一の時変係数で２行１列の値が第二の時変係数と
なる。また（３）式右辺第２項の（・）の値を位相誤差
と呼ぶと、初期位相推定値 は位相誤差に重み関数wiを乗じ、更に時刻ｋと重み関数
wiに依存する第一の時変係数を乗じ一時刻前の推定値 から減ずることにより修正し、又周波数推定値 は同じく位相誤差に重み関数wiを乗じ、更に時刻ｋと重
み関数wiに依存する第二の時変係数を乗じ一時刻前の推
定値 から減ずることにより修正をするという構成をとること
により（３）式が実現される。この方法によれば、毎時
刻その時点までの観測値を用いた最良の推定が方程式を
解く形で与えられるため、極めて短時間に最適な推定値
が得られる。
（実施例） 以下本発明の実施例について図面を参照して詳細に説
明する。第１図は本発明の実施例を示すブロック図であ
る。
明する。第１図は本発明の実施例を示すブロック図であ
る。
雑音の重畳された正弦波信号はサンプルホールド１で
サンプルされ、A/D変換器２でディジタル化される。デ
ィジタル化された正弦波は、ディジタル位相検波器３に
於て、後述する手段により推定された周波数と位相に合
致した正弦波を発生する正弦波発生器４の発生する正弦
波と比較され、位相誤差が出力される。位相検波器の実
現手段は例えば文献（W.C.Lindesy“A survey of Digit
al Phase−Locked Loop"Proc.IEEE,vol.69,No.4,April
1981）に詳しく示されている。
サンプルされ、A/D変換器２でディジタル化される。デ
ィジタル化された正弦波は、ディジタル位相検波器３に
於て、後述する手段により推定された周波数と位相に合
致した正弦波を発生する正弦波発生器４の発生する正弦
波と比較され、位相誤差が出力される。位相検波器の実
現手段は例えば文献（W.C.Lindesy“A survey of Digit
al Phase−Locked Loop"Proc.IEEE,vol.69,No.4,April
1981）に詳しく示されている。
正弦波発生器４は位相を入力としその値に対応した正
弦波の値を出力とする読みだし専用メモリ（ROM）44、
加算器43、乗算器42、時刻を発生するカウンタ41により
構成することができ、後述する位相推定値と周波数推定
値をアドレスとして正弦波の標本値を出力する。重み付
け回路５は、カウンタ54の時間経過に従い重み付けの係
数wiを係数発生器51で発生させ前記位相誤差にディジタ
ル位相変換器３の出力に乗算器52、53で重み付けする回
路である。カウンタ41、54、82は本装置の動作開始時点
から毎時刻カウントアップしている。これらのカウンタ
の値は等しいため共用する事も可能である。
弦波の値を出力とする読みだし専用メモリ（ROM）44、
加算器43、乗算器42、時刻を発生するカウンタ41により
構成することができ、後述する位相推定値と周波数推定
値をアドレスとして正弦波の標本値を出力する。重み付
け回路５は、カウンタ54の時間経過に従い重み付けの係
数wiを係数発生器51で発生させ前記位相誤差にディジタ
ル位相変換器３の出力に乗算器52、53で重み付けする回
路である。カウンタ41、54、82は本装置の動作開始時点
から毎時刻カウントアップしている。これらのカウンタ
の値は等しいため共用する事も可能である。
係数発生器８は（３）式のGk -1ykを計算し各々の係数
を発生する回路であり、カウンタ82によって得られる時
刻と各時刻得られる重み付け係数とから逆行列計算器81
によって逆行列を計算し、初期位相推定と周波数推定の
為のそれぞれ第１及び第２の係数を発生する。係数発生
器８は毎時刻逆行列の計算を行なうのでマイクロプロセ
ッサ等により実現される。
を発生する回路であり、カウンタ82によって得られる時
刻と各時刻得られる重み付け係数とから逆行列計算器81
によって逆行列を計算し、初期位相推定と周波数推定の
為のそれぞれ第１及び第２の係数を発生する。係数発生
器８は毎時刻逆行列の計算を行なうのでマイクロプロセ
ッサ等により実現される。
初期位相推定器６は、重み付けされた位相誤差と第１
の係数を乗算器62にて掛け合わせその値を初期位相推定
値を保持するレジスタ61の内容と加算器63に於て加算
し、初期位相推定値の更新値として再びレジスタに格納
する。
の係数を乗算器62にて掛け合わせその値を初期位相推定
値を保持するレジスタ61の内容と加算器63に於て加算
し、初期位相推定値の更新値として再びレジスタに格納
する。
周波数推定器７は、重み付けされた位相誤差と第２の
係数を乗算器72にて掛け合わせ、その値を周波数推定値
を保持するレジスタ71の内容と加算器73に於て加算し、
周波数推定値の更新値として再びレジスタに格納する。
更新された周波数推定値は正弦発生器４に於てカウンタ
41の内容と乗算器42で乗算され、更に初期位相推定値と
加算器43で加算され、時刻ｋに於ける最適位相推定値と
して読みだし専用メモリ44に入力され正弦波の値が出力
される。尚、この読みだし専用メモリのアドレスの最大
値を２Πに対応させておくことにより位相推定値はモジ
ュロ２Πで行なわれる。以上により、入力される雑音の
重畳した正弦波から周波数と位相を逐次的に最小分散推
定をすることができる。推定値はそれぞれレジスタ61、
71の内容として与えられる。
係数を乗算器72にて掛け合わせ、その値を周波数推定値
を保持するレジスタ71の内容と加算器73に於て加算し、
周波数推定値の更新値として再びレジスタに格納する。
更新された周波数推定値は正弦発生器４に於てカウンタ
41の内容と乗算器42で乗算され、更に初期位相推定値と
加算器43で加算され、時刻ｋに於ける最適位相推定値と
して読みだし専用メモリ44に入力され正弦波の値が出力
される。尚、この読みだし専用メモリのアドレスの最大
値を２Πに対応させておくことにより位相推定値はモジ
ュロ２Πで行なわれる。以上により、入力される雑音の
重畳した正弦波から周波数と位相を逐次的に最小分散推
定をすることができる。推定値はそれぞれレジスタ61、
71の内容として与えられる。
尚、プログラマブルな信号処理プロセッサで上記と等
価な動作を容易に実現することができることは明らかで
ある。
価な動作を容易に実現することができることは明らかで
ある。
（発明の効果） 本発明によれば、従来技術よりも特性良くしかも正確
に周波数・位相を推定できるという効果がある。
に周波数・位相を推定できるという効果がある。
第１図は本発明の実施例を示すブロック図である。図
中、 １……サンプルホールド、２……A/D変換器、 ３……位相検波器、４……正弦波発生器、 ５……重み付け回路、６……初期位相推定器、 ７……周波数推定器、８……係数発生器である。
中、 １……サンプルホールド、２……A/D変換器、 ３……位相検波器、４……正弦波発生器、 ５……重み付け回路、６……初期位相推定器、 ７……周波数推定器、８……係数発生器である。
Claims (1)
- 【請求項１】正弦波あるいは雑音の重畳した正弦波を入
力とし、入力正弦波の初期位相と周波数を推定する装置
において、 入力正弦波の初期位相の推定値と入力正弦波の周波数の
推定値から次時刻における正弦波の値を発生する手段
と、該正弦波と前期入力正弦波との位相誤差を検出する
手段と、推定開始から一定時間だけ重み付け値を発生す
る手段と、検出された位相誤差に前記重み付け値を乗算
する手段と、前記重み付け値を入力し第一及び第二の時
変係数を発生する係数発生手段と、前記重み付けされた
位相誤差に第一の時変係数を乗じた値により前記初期位
相の推定値を逐次的に更新する初期位相推定手段と、前
記重み付けされた位相誤差に第２の時変係数を乗じた値
により前記周波数の推定値を逐次的に更新する周波数推
定手段とを有することを特徴とする周波数・位相推定装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63250554A JP2527011B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63250554A JP2527011B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0296664A JPH0296664A (ja) | 1990-04-09 |
| JP2527011B2 true JP2527011B2 (ja) | 1996-08-21 |
Family
ID=17209639
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63250554A Expired - Lifetime JP2527011B2 (ja) | 1988-10-03 | 1988-10-03 | 周波数・位相推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2527011B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008136443A1 (ja) * | 2007-04-26 | 2008-11-13 | The University Of Tokyo | 正弦波パラメータ推定方法 |
| CN113219248B (zh) * | 2021-05-07 | 2022-09-20 | 南京大学 | 基于时域波形比对的信号分量估计方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0716159B2 (ja) * | 1984-12-27 | 1995-02-22 | 日本電気株式会社 | 周波数・位相推定装置 |
-
1988
- 1988-10-03 JP JP63250554A patent/JP2527011B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0296664A (ja) | 1990-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5444639A (en) | Angular rate sensing system and method, with digital synthesizer and variable-frequency oscillator | |
| US6337888B1 (en) | In-phase and quadrature signal regeneration | |
| US20040085096A1 (en) | Efficient digital method of and system for determining the instantaneous phase and amplitude of a vibratory accelerometer and other sensors | |
| EP0490275A2 (en) | System and method for compensation of in-phase and quadrature phase and gain imbalance | |
| EP0638811B1 (en) | Estimator of amplitude and frequency of a noisy biased sinusoid from short bursts of samples | |
| JPH07170152A (ja) | 入力信号をフィルタ処理する方法、およびトラッキングフィルタ | |
| US5576976A (en) | Amplitude detection and automatic gain control of a sparsely sampled sinusoid by adjustment of a notch filter | |
| Franz et al. | Frequency estimation using warped discrete Fourier transform | |
| Adjrad et al. | Estimation of multicomponent polynomial-phase signals impinging on a multisensor array using state–space modeling | |
| JP2527011B2 (ja) | 周波数・位相推定装置 | |
| US5487015A (en) | Self-oscillating driver circuit for a quartz resonator of an angular rate sensor | |
| Maskell et al. | Adaptive subsample delay estimation using a modified quadrature phase detector | |
| JP2527008B2 (ja) | 周波数・位相推定装置 | |
| US5790412A (en) | Recursive frequency aging estimation and prediction devices, methods and computer programs for crystal oscillators | |
| Chakraborty | A new adaptive filter for estimating and tracking the delay and the amplitude of a sinusoid | |
| JPH0716159B2 (ja) | 周波数・位相推定装置 | |
| US7133791B1 (en) | Two-mean level-crossing time interval estimation method | |
| JP2527009B2 (ja) | 周波数・位相推定装置 | |
| Chakraborty et al. | New adaptive algorithm for delay estimation of sinusoidal signals | |
| US20020191688A1 (en) | Adaptive line enhancer | |
| JP4291563B2 (ja) | パラメータの共用推定のための方法 | |
| JP4344356B2 (ja) | 検波装置、方法、プログラム、記録媒体 | |
| JPH0138270B2 (ja) | ||
| JP2789600B2 (ja) | 周波数推定方式 | |
| JP3402375B2 (ja) | Ｄｑｐｓｋ変調信号評価装置 |