JP3263091B2 - 周波数分析装置 - Google Patents
周波数分析装置Info
- Publication number
- JP3263091B2 JP3263091B2 JP40447190A JP40447190A JP3263091B2 JP 3263091 B2 JP3263091 B2 JP 3263091B2 JP 40447190 A JP40447190 A JP 40447190A JP 40447190 A JP40447190 A JP 40447190A JP 3263091 B2 JP3263091 B2 JP 3263091B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- signal
- instantaneous
- order correlation
- output
- 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 - Fee Related
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Description
【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、音声信号や音響信号等
の時系列信号の認識を行う前処理等において、その時系
列信号の周波数を抽出する周波数分析装置に関するもの
である。
の時系列信号の認識を行う前処理等において、その時系
列信号の周波数を抽出する周波数分析装置に関するもの
である。
【０００２】
【従来の技術】従来、この種の分野の技術としては、安
居院・中嶋著「コンピュータ音声処理」（１９８６−１
−１０）秋葉出版、Ｐ．１１４−１２７に記載されるも
のがあった。以下、その構成を図を用いて説明する。
居院・中嶋著「コンピュータ音声処理」（１９８６−１
−１０）秋葉出版、Ｐ．１１４−１２７に記載されるも
のがあった。以下、その構成を図を用いて説明する。
【０００３】図２は従来の周波数分析装置の一構成例を
示す構成ブロック図、及び図３及び図４は図２の信号波
形図である。なお、図３は入力時系列信号ｘ（ｔ）、及
び図４は入力時系列信号ｘ（ｔ）の周波数特性ｘ（ｔ）
を示す図である。
示す構成ブロック図、及び図３及び図４は図２の信号波
形図である。なお、図３は入力時系列信号ｘ（ｔ）、及
び図４は入力時系列信号ｘ（ｔ）の周波数特性ｘ（ｔ）
を示す図である。
【０００４】図３に示す入力時系列信号ｘ（ｔ）の周波
数ｆo を分析して抽出する場合、該入力時系列信号ｘ
（ｔ）をアナログ／ディジタル変換器（以下、Ａ／Ｄ変
換器という）１によりディジタル信号に変換した後、フ
ーリエ変換器２へ与える。フーリエ変換器２では、Ａ／
Ｄ変換器１の出力に対して高速フーリエ変換（以下、Ｆ
ＦＴという）等のアルゴリリズムを適用し、図４に示す
ような入力時系列信号ｘ（ｔ）に対する周波数特性Ｘ
（ｆ）を求め、ピーク検出器３へ送る。ピーク検出器３
では、周波数特性Ｘ（ｆ）のピークを検出することによ
り、そのピークに対応する周波数ｆo を入力時系列信号
ｘ（ｔ）の周波数として出力する。
数ｆo を分析して抽出する場合、該入力時系列信号ｘ
（ｔ）をアナログ／ディジタル変換器（以下、Ａ／Ｄ変
換器という）１によりディジタル信号に変換した後、フ
ーリエ変換器２へ与える。フーリエ変換器２では、Ａ／
Ｄ変換器１の出力に対して高速フーリエ変換（以下、Ｆ
ＦＴという）等のアルゴリリズムを適用し、図４に示す
ような入力時系列信号ｘ（ｔ）に対する周波数特性Ｘ
（ｆ）を求め、ピーク検出器３へ送る。ピーク検出器３
では、周波数特性Ｘ（ｆ）のピークを検出することによ
り、そのピークに対応する周波数ｆo を入力時系列信号
ｘ（ｔ）の周波数として出力する。
【０００５】
【発明が解決しようとする課題】しかしながら、上記構
成の周波数分析装置では、次のような課題があった。従
来の装置では、まず、フーリエ変換器２によって入力時
系列信号Ｘ（ｔ）を求めた後、ピーク検出器３で周波数
ｆo を求めている。ここで、フーリエ変換器２では、離
散的な信号値からフーリエ係数を求める離散的フーリエ
変換（以下、ＤＦＴという）を用いて周波数特性Ｘを求
めても良い。しかし、ＤＦＴを用いると、データ数Ｎ＝
２L （但し、Ｌ；パラメータ）の場合、Ｎ2 回の乗算回
数を必要とし、演算時間が長くなる。そこで、演算時間
を大幅に減少させるために、ＦＦＴを用いると、（Ｎ／
２）・Ｌの乗算回数でスペクトルが求まる。
成の周波数分析装置では、次のような課題があった。従
来の装置では、まず、フーリエ変換器２によって入力時
系列信号Ｘ（ｔ）を求めた後、ピーク検出器３で周波数
ｆo を求めている。ここで、フーリエ変換器２では、離
散的な信号値からフーリエ係数を求める離散的フーリエ
変換（以下、ＤＦＴという）を用いて周波数特性Ｘを求
めても良い。しかし、ＤＦＴを用いると、データ数Ｎ＝
２L （但し、Ｌ；パラメータ）の場合、Ｎ2 回の乗算回
数を必要とし、演算時間が長くなる。そこで、演算時間
を大幅に減少させるために、ＦＦＴを用いると、（Ｎ／
２）・Ｌの乗算回数でスペクトルが求まる。
【０００６】ところが、この様なＦＦＴを用いた周波数
分析装置において、高精度に周波数ｆo を求めるために
は、周波数特性Ｘ（ｆ）をＦＦＴを用いて高精度に求め
なければならない。しかし、周波数特性Ｘ（ｆ）を高精
度に求めるためには、パラメータＬを大きくしなければ
ならないが、それによって演算量が大きくなり、処理時
間が長くなるという問題があり、それを解決することが
困難であった。本発明は前記従来技術の持っていた課題
として、精度の高い周波数を求めるためには演算量が大
きくなって処理時間が長くなるという点について解決し
た周波数分析装置を提供するものである。
分析装置において、高精度に周波数ｆo を求めるために
は、周波数特性Ｘ（ｆ）をＦＦＴを用いて高精度に求め
なければならない。しかし、周波数特性Ｘ（ｆ）を高精
度に求めるためには、パラメータＬを大きくしなければ
ならないが、それによって演算量が大きくなり、処理時
間が長くなるという問題があり、それを解決することが
困難であった。本発明は前記従来技術の持っていた課題
として、精度の高い周波数を求めるためには演算量が大
きくなって処理時間が長くなるという点について解決し
た周波数分析装置を提供するものである。
【０００７】
【課題を解決するための手段】前記課題を解決するため
に、本発明のうちの第１の発明は、所定の時間間隔Ｔｓ
で離散化された時系列信号である入力信号列Ｘ n （但
し、ｎは時刻、ｎ＝０，１，２，…）の周波数ｆ n を抽
出する周波数分析装置において、前記入力信号列Ｘ n の
振幅値を２乗して瞬時パワーＰ n （＝Ｘ n 2 ）を求める
瞬時パワー演算手段と、前記入力信号列Ｘ n を前記１時
間間隔Ｔｓだけ遅延させ、この遅延した信号Ｘ n-1 の振
幅値と該入力信号列Ｘ n の振幅値とを乗算して瞬時１次
相関関数Ｃ n （＝Ｘ n ・Ｘ n-1 ）を求める瞬時１次相関
演算手段と、前記瞬時パワーＰ n の高周波成分を除去し
て除去後の信号Ｐ na を出力する第１の高周波除去手段
と、前記瞬時１次相関関数Ｃ n の高周波成分を除去して
除去後の信号Ｃ na を出力する第２の高周波除去手段と、
前記信号Ｃ na を前記信号Ｐ na で除算することにより正規
化して１次相関係数γ n （＝Ｃ na ／Ｐ na ）を求める正規
化手段と、前記１次相関係数γ n の逆余弦を演算して前
記入力信号列Ｘ n の周波数ｆ n （＝（１／２πＴｓ）・
arccos・γ n 、但しarccosは逆余弦関数）を求める周波
数演算手段と、を有している。
に、本発明のうちの第１の発明は、所定の時間間隔Ｔｓ
で離散化された時系列信号である入力信号列Ｘ n （但
し、ｎは時刻、ｎ＝０，１，２，…）の周波数ｆ n を抽
出する周波数分析装置において、前記入力信号列Ｘ n の
振幅値を２乗して瞬時パワーＰ n （＝Ｘ n 2 ）を求める
瞬時パワー演算手段と、前記入力信号列Ｘ n を前記１時
間間隔Ｔｓだけ遅延させ、この遅延した信号Ｘ n-1 の振
幅値と該入力信号列Ｘ n の振幅値とを乗算して瞬時１次
相関関数Ｃ n （＝Ｘ n ・Ｘ n-1 ）を求める瞬時１次相関
演算手段と、前記瞬時パワーＰ n の高周波成分を除去し
て除去後の信号Ｐ na を出力する第１の高周波除去手段
と、前記瞬時１次相関関数Ｃ n の高周波成分を除去して
除去後の信号Ｃ na を出力する第２の高周波除去手段と、
前記信号Ｃ na を前記信号Ｐ na で除算することにより正規
化して１次相関係数γ n （＝Ｃ na ／Ｐ na ）を求める正規
化手段と、前記１次相関係数γ n の逆余弦を演算して前
記入力信号列Ｘ n の周波数ｆ n （＝（１／２πＴｓ）・
arccos・γ n 、但しarccosは逆余弦関数）を求める周波
数演算手段と、を有している。
【０００８】第２の発明は、第１の発明の周波数分析装
置において、前記第１の高周波除去手段は、インパルス
応答Ｗ n により前記瞬時パワーＰ n を畳み込み演算して
前記 信号Ｐ na を出力し、前記第２の高周波除去手段は、
インパルス応答Ｗ n により前記瞬時１次相関関数Ｃ n を
畳み込み演算して前記信号Ｃ na を出力する構成にしてい
る。
置において、前記第１の高周波除去手段は、インパルス
応答Ｗ n により前記瞬時パワーＰ n を畳み込み演算して
前記 信号Ｐ na を出力し、前記第２の高周波除去手段は、
インパルス応答Ｗ n により前記瞬時１次相関関数Ｃ n を
畳み込み演算して前記信号Ｃ na を出力する構成にしてい
る。
【０００９】
【作用】第１及び第２の発明は、以上のように周波数分
析装置を構成したので、時系列信号である入力信号列Ｘ
n が入力されると、瞬時パワー演算手段及び瞬時１次相
関演算手段により、その入力信号列Ｘ n から瞬時パワー
Ｐ n 及び瞬時１次相関関数Ｃ n が求められる。さらに、
その瞬時パワーＰ n 及び瞬時１次相関関数Ｃ n の高周波
成分が第１及び第２の高周波除去手段によってそれぞれ
除去され、この第１及び第２の高周波除去手段の各出力
信号Ｐ na ，Ｃ na から、正規化手段によって１次相関係数
γ n が求められる。その後、１次相関係数γ n から、周
波数演算手段によって入力信号列Ｘ n の周波数ｆ n が求
められる。これにより、従来装置のように、フーリエ変
換を用いずに周波数ｆ n を抽出でき、少ない演算量で高
精度に入力信号列Ｘ n の周波数ｆ n を抽出することがで
きる。従って、前記課題を解決できるのである。
析装置を構成したので、時系列信号である入力信号列Ｘ
n が入力されると、瞬時パワー演算手段及び瞬時１次相
関演算手段により、その入力信号列Ｘ n から瞬時パワー
Ｐ n 及び瞬時１次相関関数Ｃ n が求められる。さらに、
その瞬時パワーＰ n 及び瞬時１次相関関数Ｃ n の高周波
成分が第１及び第２の高周波除去手段によってそれぞれ
除去され、この第１及び第２の高周波除去手段の各出力
信号Ｐ na ，Ｃ na から、正規化手段によって１次相関係数
γ n が求められる。その後、１次相関係数γ n から、周
波数演算手段によって入力信号列Ｘ n の周波数ｆ n が求
められる。これにより、従来装置のように、フーリエ変
換を用いずに周波数ｆ n を抽出でき、少ない演算量で高
精度に入力信号列Ｘ n の周波数ｆ n を抽出することがで
きる。従って、前記課題を解決できるのである。
【００１０】
【実施例】図１は、本発明の実施例を示す周波数分析装
置の構成ブロック図である。この周波数分析装置は、音
声信号や音響信号等のアナログ信号である時系列信号ｘ
（ｔ）を入力するための入力端子５１を有し、その入力
端子５１がアナログ／ディジタル変換用のＡ／Ｄ変換器
５２に接続されている。さらに、Ａ／Ｄ変換器５２の出
力側には、瞬時パワー演算手段である瞬時パワー演算回
路５３と、瞬時１次相関演算手段である瞬時１次相関演
算回路５４が接続されている。
置の構成ブロック図である。この周波数分析装置は、音
声信号や音響信号等のアナログ信号である時系列信号ｘ
（ｔ）を入力するための入力端子５１を有し、その入力
端子５１がアナログ／ディジタル変換用のＡ／Ｄ変換器
５２に接続されている。さらに、Ａ／Ｄ変換器５２の出
力側には、瞬時パワー演算手段である瞬時パワー演算回
路５３と、瞬時１次相関演算手段である瞬時１次相関演
算回路５４が接続されている。
【００１１】瞬時パワー演算回路５３は、Ａ／Ｄ変換器
５２の出力から瞬時パワーを求める回路であり、Ａ／Ｄ
変換器５２の出力の２乗を演算して出力する乗算器５３
ａで構成されている。瞬時１次相関演算回路５４は、Ａ
／Ｄ変換器５２の出力から瞬時１次相関関数を求める回
路であり、Ａ／Ｄ変換器５２の出力を１サンプル遅延さ
せる遅延器５４ａと、Ａ／Ｄ変換器５２の出力及び遅延
器５４ａの出力の乗算を演算する乗算器５４ｂとで構成
されている。乗算器５３ａ，５４ｂの出力側には、同一
のフィルタ特性を有する第１の高周波除去手段であるロ
ーパスフィルタ５５と第２の高周波除去手段であるロー
パスフィルタ５６がそれぞれ接続されている。
５２の出力から瞬時パワーを求める回路であり、Ａ／Ｄ
変換器５２の出力の２乗を演算して出力する乗算器５３
ａで構成されている。瞬時１次相関演算回路５４は、Ａ
／Ｄ変換器５２の出力から瞬時１次相関関数を求める回
路であり、Ａ／Ｄ変換器５２の出力を１サンプル遅延さ
せる遅延器５４ａと、Ａ／Ｄ変換器５２の出力及び遅延
器５４ａの出力の乗算を演算する乗算器５４ｂとで構成
されている。乗算器５３ａ，５４ｂの出力側には、同一
のフィルタ特性を有する第１の高周波除去手段であるロ
ーパスフィルタ５５と第２の高周波除去手段であるロー
パスフィルタ５６がそれぞれ接続されている。
【００１２】ローパスフィルタ５５，５６は、それぞれ
乗算器５３ａ，５４ｂの出力の高周波成分を除去して出
力する機能を有し、その出力側には正規化手段である正
規化回路５７と周波数演算手段である周波数演算回路５
８が接続されている。正規化回路５７は、除算器等で構
成され、ローパスフィルタ５５，５６の出力から１次相
関係数を演算する機能を有し、周波数演算回路５８は、
逆余弦関数演算回路等で構成され、正規化回路５７の出
力から入力信号列Ｘ n の周波数ｆ n を演算する回路であ
る。
乗算器５３ａ，５４ｂの出力の高周波成分を除去して出
力する機能を有し、その出力側には正規化手段である正
規化回路５７と周波数演算手段である周波数演算回路５
８が接続されている。正規化回路５７は、除算器等で構
成され、ローパスフィルタ５５，５６の出力から１次相
関係数を演算する機能を有し、周波数演算回路５８は、
逆余弦関数演算回路等で構成され、正規化回路５７の出
力から入力信号列Ｘ n の周波数ｆ n を演算する回路であ
る。
【００１３】以上のように構成される周波数分析装置の
動作について説明する。入力端子５１から入力された時
系列信号ｘ（ｔ）は、Ａ／Ｄ変換器５２により予め定め
られたサンプリング時間間隔Ｔｓ毎に離散化され、入力
信号列Ｘn （ｎ＝０，１，２，……）に変換される。乗
算器５３ａは、Ａ／Ｄ変換器５２の出力Ｘn により、時
刻ｎの瞬時パワーＰn を、 Ｐn ＝Ｘn 2 ・・・（１） により演算する。
動作について説明する。入力端子５１から入力された時
系列信号ｘ（ｔ）は、Ａ／Ｄ変換器５２により予め定め
られたサンプリング時間間隔Ｔｓ毎に離散化され、入力
信号列Ｘn （ｎ＝０，１，２，……）に変換される。乗
算器５３ａは、Ａ／Ｄ変換器５２の出力Ｘn により、時
刻ｎの瞬時パワーＰn を、 Ｐn ＝Ｘn 2 ・・・（１） により演算する。
【００１４】遅延器５４ａは、Ａ／Ｄ変換器５２の出力
Ｘn を入力し、１サンプル遅延させて出力Ｘn-1 を出力
する。乗算器５４ｂは、Ａ／Ｄ変換器５２の出力Ｘn 及
び遅延器５４ａの出力Ｘn-1 から、瞬時１次相関関数Ｃ
n を、 Ｃn ＝Ｘn ・Ｘn-1 ・・・（２） により演算する。
Ｘn を入力し、１サンプル遅延させて出力Ｘn-1 を出力
する。乗算器５４ｂは、Ａ／Ｄ変換器５２の出力Ｘn 及
び遅延器５４ａの出力Ｘn-1 から、瞬時１次相関関数Ｃ
n を、 Ｃn ＝Ｘn ・Ｘn-1 ・・・（２） により演算する。
【００１５】ローパスフィルタ５５，５６は、各乗算器
５３ａ，５４ｂの出力である瞬時パワーＰn 及び瞬時１
次相関関数Ｃn をそれぞれ入力し、瞬時パワーＰn 及び
瞬時１次相関関数Ｃn の高周波成分を除去する。即ち、
ローパスフィルタ５５，５６のインパルス応答をＷn と
すると、ローパスフィルタ５５，５６は、次式（３）の
演算により、Ｐna，Ｃnaを演算して出力する。
５３ａ，５４ｂの出力である瞬時パワーＰn 及び瞬時１
次相関関数Ｃn をそれぞれ入力し、瞬時パワーＰn 及び
瞬時１次相関関数Ｃn の高周波成分を除去する。即ち、
ローパスフィルタ５５，５６のインパルス応答をＷn と
すると、ローパスフィルタ５５，５６は、次式（３）の
演算により、Ｐna，Ｃnaを演算して出力する。
【００１６】
【数１】 正規化回路５７は、ローパスフィルタ５５，５６の各出
力であるＰna，Ｃnaから１次相関係数γn を次式（４）
により演算して出力する。 γn ＝Ｃna／Ｐna ・・・（４） １次相関係数γn は出力Ｃnaを出力Ｐnaで正規化したも
のに相当する。周波数演算回路５８は、正規化回路５７
で演算された１次相関係数γn により、周波数ｆn を次
式（５）により演算する。
力であるＰna，Ｃnaから１次相関係数γn を次式（４）
により演算して出力する。 γn ＝Ｃna／Ｐna ・・・（４） １次相関係数γn は出力Ｃnaを出力Ｐnaで正規化したも
のに相当する。周波数演算回路５８は、正規化回路５７
で演算された１次相関係数γn により、周波数ｆn を次
式（５）により演算する。
【００１７】
【数２】 arccosは逆余弦関数であり、ＴｓはＡ／Ｄ変換器５２に
おける予め定めたサンプリング時間間隔である。周波数
演算回路５８で演算周波数ｆn を時刻ｎの周波数として
出力する。
おける予め定めたサンプリング時間間隔である。周波数
演算回路５８で演算周波数ｆn を時刻ｎの周波数として
出力する。
【００１８】本実施例は、次のような利点を有してい
る。図５は、本実施例による周波数分析結果の一例であ
る。信号として、１００Ｈｚから３００Ｈｚまで、直線
的に周波数が変化する信号を入力し、サンプリング時間
間隔１／８００秒で離散化した後、上述の方法により求
めた周波数をプロットしたものである。横軸はサンプル
番号であり、サンプリング時間間隔１／８００秒をかけ
ることにより、実際の時刻（秒）となる。縦軸は周波数
（Ｈｚ）である。同図より、本実施例によれば、高精度
に周波数を求められることが分かる。なお、本発明は、
図示の実施例に限定されず、種々の変形が可能である。
例えば、図１の周波数分析装置をプロセッサを用いたプ
ログラム制御等により実行するようにすることも可能で
ある。
る。図５は、本実施例による周波数分析結果の一例であ
る。信号として、１００Ｈｚから３００Ｈｚまで、直線
的に周波数が変化する信号を入力し、サンプリング時間
間隔１／８００秒で離散化した後、上述の方法により求
めた周波数をプロットしたものである。横軸はサンプル
番号であり、サンプリング時間間隔１／８００秒をかけ
ることにより、実際の時刻（秒）となる。縦軸は周波数
（Ｈｚ）である。同図より、本実施例によれば、高精度
に周波数を求められることが分かる。なお、本発明は、
図示の実施例に限定されず、種々の変形が可能である。
例えば、図１の周波数分析装置をプロセッサを用いたプ
ログラム制御等により実行するようにすることも可能で
ある。
【００１９】
【発明の効果】以上詳細に説明したように、第１及び第
２の発明によれば、入力信号列Ｘ n から瞬時パワーＰ n
及び瞬時１次相関関数Ｃ n を求め、その瞬時パワーＰ n
と瞬時１次相関関数Ｃ n との高周波成分を除去して１次
相関係数γ n を求めた後、入力信号列Ｘ n の周波数ｆ n
を求めるようにしたので、フーリエ変換を用いずに入力
信号列Ｘ n の周波数ｆ n を抽出できる。これにより、少
ない演算量で高速処理が可能となると共に、単純な演算
方法にもかかわらず、高精度に入力信号列Ｘ n の周波数
ｆ n を抽出することができる。
２の発明によれば、入力信号列Ｘ n から瞬時パワーＰ n
及び瞬時１次相関関数Ｃ n を求め、その瞬時パワーＰ n
と瞬時１次相関関数Ｃ n との高周波成分を除去して１次
相関係数γ n を求めた後、入力信号列Ｘ n の周波数ｆ n
を求めるようにしたので、フーリエ変換を用いずに入力
信号列Ｘ n の周波数ｆ n を抽出できる。これにより、少
ない演算量で高速処理が可能となると共に、単純な演算
方法にもかかわらず、高精度に入力信号列Ｘ n の周波数
ｆ n を抽出することができる。
【図面の簡単な説明】
【図１】本発明の実施例の周波数分析装置の構成ブロッ
ク図である。
ク図である。
【図２】従来の周波数分析装置の構成ブロック図であ
る。
る。
【図３】図２の信号波形図である。
【図４】図２の信号波形図である。
【図５】本実施例による周波数分析結果の一例を示す図
である。
である。
５３ 瞬時パワー演算回路 ５４ 瞬時１次相関演算回路 ５５ 第１のローパスフィルタ ５６ 第２のローパスフィルタ ５７ 正規化回路 ５８ 周波数演算回路 ｘ（ｔ） 時系列信号
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭58−143394（ＪＰ，Ａ) 特開 昭58−220200（ＪＰ，Ａ) 特開 昭59−212898（ＪＰ，Ａ) 特開 昭61−32097（ＪＰ，Ａ) 特開 昭61−32098（ＪＰ，Ａ) 特開 昭61−32099（ＪＰ，Ａ) 特開 昭62−180398（ＪＰ，Ａ) 特開 昭63−261180（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) G10L 11/00 - 11/04
Claims (2)
- 【請求項１】 所定の時間間隔Ｔｓで離散化された時系
列信号である入力信号列Ｘ n （但し、ｎは時刻、ｎ＝
０，１，２，…）の周波数ｆ n を抽出する周波数分析装
置において、 前記入力信号列Ｘ n の振幅値を２乗して瞬時パワーＰ n
（＝Ｘ n 2 ）を求める瞬時パワー演算手段と、 前記入力信号列Ｘ n を前記１時間間隔Ｔｓだけ遅延さ
せ、この遅延した信号Ｘ n-1 の振幅値と該入力信号列Ｘ
n の振幅値とを乗算して瞬時１次相関関数Ｃ n （＝Ｘ n
・Ｘ n-1 ）を求める瞬時１次相関演算手段と、 前記瞬時パワーＰ n の高周波成分を除去して除去後の信
号Ｐ na を出力する第１の高周波除去手段と、 前記瞬時１次相関関数Ｃ n の高周波成分を除去して除去
後の信号Ｃ na を出力する第２の高周波除去手段と、 前記信号Ｃ na を前記信号Ｐ na で除算することにより正規
化して１次相関係数γ n （＝Ｃ na ／Ｐ na ）を求める正規
化手段と、 前記１次相関係数γ n の逆余弦を演算して前記入力信号
列Ｘ n の周波数ｆ n （＝（１／２πＴｓ）・arccos・γ
n 、但しarccosは逆余弦関数）を求める周波数演算手段
と、 を有することを特徴とする周波数分析装置。 - 【請求項２】 前記第１の高周波除去手段は、インパル
ス応答Ｗ n により前記瞬時パワーＰ n を畳み込み演算し
て前記信号Ｐ na を出力し、前記第２の高周波除去手段
は、インパルス応答Ｗ n により前記瞬時１次相関関数Ｃ
n を畳み込み演算して前記信号Ｃ na を出力する構成にし
たことを特徴とする請求項１記載の周波数分析装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP40447190A JP3263091B2 (ja) | 1990-12-20 | 1990-12-20 | 周波数分析装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP40447190A JP3263091B2 (ja) | 1990-12-20 | 1990-12-20 | 周波数分析装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH04219800A JPH04219800A (ja) | 1992-08-10 |
| JP3263091B2 true JP3263091B2 (ja) | 2002-03-04 |
Family
ID=18514141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP40447190A Expired - Fee Related JP3263091B2 (ja) | 1990-12-20 | 1990-12-20 | 周波数分析装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3263091B2 (ja) |
-
1990
- 1990-12-20 JP JP40447190A patent/JP3263091B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH04219800A (ja) | 1992-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US3973112A (en) | System for determining transfer function | |
| Boashash et al. | An efficient real-time implementation of the Wigner-Ville distribution | |
| US4901244A (en) | Apparatus for, and method of, analyzing signals | |
| Nguyen | Parametric harmonic analysis | |
| Papoulis et al. | Detection of hidden periodicities by adaptive extrapolation | |
| US20050094344A1 (en) | Method of compensating for distorted secondary current of current transformer | |
| CA1228673A (en) | Apparatus for and method of, analysing signals | |
| JPH0389174A (ja) | 電子測定装置及び周波数推定方法 | |
| JP3263091B2 (ja) | 周波数分析装置 | |
| CA2142509C (en) | Complex cepstrum analyzer for speech signals | |
| JP3263067B2 (ja) | 周波数分析方法 | |
| US9268745B2 (en) | Method for fast wavelet transform of a signal | |
| JP4344356B2 (ja) | 検波装置、方法、プログラム、記録媒体 | |
| Bittanti et al. | Poles identification of an analog filter for nuclear spectroscopy via subspace-based techniques | |
| JP2957572B1 (ja) | 地震応答スペクトル演算装置 | |
| JP3234365B2 (ja) | 応答特性計測装置 | |
| JPH0464059A (ja) | 分析データ処理装置 | |
| US20020101228A1 (en) | Method and apparatus for determining the energy of a signal | |
| JP2770046B2 (ja) | 分解能偏り誤差の推定方法およびその装置 | |
| JP2705064B2 (ja) | 線形予測分析装置 | |
| GB2212619A (en) | Analyzing signals | |
| JPH06300888A (ja) | 原子炉出力解析装置 | |
| Linggard et al. | High-speed digital filter bank | |
| JP2888916B2 (ja) | スペクトル分析方法及びその装置 | |
| JPS5858625B2 (ja) | 信号処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20011211 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081221 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091221 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |