JP2009524812A - 信号解析器 - Google Patents

信号解析器 Download PDF

Info

Publication number
JP2009524812A
JP2009524812A JP2008551879A JP2008551879A JP2009524812A JP 2009524812 A JP2009524812 A JP 2009524812A JP 2008551879 A JP2008551879 A JP 2008551879A JP 2008551879 A JP2008551879 A JP 2008551879A JP 2009524812 A JP2009524812 A JP 2009524812A
Authority
JP
Japan
Prior art keywords
signal
filter
data
curve
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.)
Pending
Application number
JP2008551879A
Other languages
English (en)
Inventor
シンプソン,マーティン,ポール
Original Assignee
ガッディディース リミテッド
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 by ガッディディース リミテッド filed Critical ガッディディース リミテッド
Publication of JP2009524812A publication Critical patent/JP2009524812A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • G01R23/165Spectrum analysis; Fourier analysis using filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0007Frequency selective voltage or current level measuring
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

信号を処理するための方法が、信号を受信するステップと、各バンドパスフィルタが少なくとも1/√2のQファクタを有する複数のバンドパスフィルタで信号をフィルタリングするステップと、フィルタの出力からデータを抽出するステップと、を具える。主要なデータ抽出が、フィルタの出力のローカルピークを特定するステップを具えており、それぞれのローカルピークが時間及び周波数方向の双方のローカルピークである。さらに本方法は、カーブフィッティングの手順を実行するステップを有利に具えており、当該手順の出力が、抽出データの一部に少なくとも近似する曲線を規定する曲線データを具えている。さらに、本方法は、信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更するステップを具えている。
【選択図】図5

Description

本発明は、信号を処理するための方法及びシステムに関する。
全ての信号が少なくとも1の正弦波信号を具えていることが知られている。このような信号は、任意の周波数範囲(mHzからEHzの範囲の周波数帯域の中の例えば、亜音速、可聴周波、ビデオ、無線通信、レーダー、可視光線、x線、ガンマ線)にあり、例えば信号の特性を特定及び/又はデータを圧縮するように、多くの場合、解析、測定、及び/又は特定を必要とする。
フーリエ解析及びFFT(高速フーリエ変換)といった従来の信号解析法は、時間窓の中の入力信号の周波数の測定値を提供するが、この方法では、入力信号の中の正弦波を微細な時間分解能に対して特定及び測定できない。さらに、過渡信号又は衝撃信号の性質により、時間内に局在していない窓の中の広帯域の周波数情報をもたらし正弦波情報を不鮮明にする。
信号を処理するために複数のバンドパスフィルタの使用することが知られている。例えば、GB 2034056がスペクトル表示装置を開示している。この装置では、複数のバンドパスフィルタが周波数帯域に入力信号を分割して、それぞれの帯域が、制御信号に応じて対応する周波数帯域における信号レベルを保持するための各保持回路に与えられる。保持回路から受信する対応する周波数帯域の信号レベルを表示するためのセクションに表示がなされる。所定の時間周期で制御信号を与える動作手段のための制御スイッチが設けられている。保持回路は通常は非作動で、ディスプレイは周波数帯域における信号の瞬間的なレベルを表示する。制御スイッチを作動させると、保持回路が周波数帯域におけるレベルを保持して、ディスプレイが所定の時間周期これらの保持レベルを表示し、その後、この装置は瞬間的なレベルの表示に自動的に復帰する。
この開示に関する装置は、受信信号の中の複数の周波数帯域の振幅レベルを表示する。これは、受信信号の中のごく一部分の成分の比較的簡単な解析であり、入力信号の非常に単純な解析以外に使用することができない。
このため、本発明の目的は従来技術を改良することである。
本発明に係る第1の態様によれば、信号を受信するステップと、複数のバンドパスフィルタで信号をフィルタリングするステップと、フィルタの出力からデータを抽出するステップと、を具えた信号の処理方法であって、データを抽出するステップが、フィルタの出力の中の局所ピークであって、それぞれの局所ピークが時間及び周波数の双方の方向の局所ピークを特定するステップを具える方法が提供されている。
本発明に係る第2の態様によれば、信号を受信するための受信部と、信号をフィルタリングするための複数のバンドパスフィルタと、フィルタの出力からデータを抽出するよう構成された処理装置と、を具えた信号を処理するためのシステムであって、データ抽出が、フィルタ出力の中の局所ピークであって、それぞれの局所ピークが時間及び周波数の双方の方向の局所ピークを特定することを具えているシステムが提供されている。
本発明によれば、微細、正確且つ詳細に信号の成分を表す信号を処理する方法を提供することが可能となる。本発明は、非常に高い時間分解能で正弦波信号の存在、変化を判別する能力を与える。さらに、本発明は、高度に正確な周波数の分解能で正弦波信号の周波数を判別する能力を与える。さらに、本発明は、高い精度で複数の正弦波信号を識別する能力を与える。また、受信信号の中の正弦波信号に関する位相情報を取得可能で有る。本方法を使用して、信号の中の周期的成分を検出することが可能であり、さらに、信号の中の衝撃成分を検出することができる。例えば、音声信号では、手をたたくことが衝撃成分であり、バンドパスフィルタの出力から局所ピークを抽出する方法を使用してこれを特定できる。
本発明の第1の利点は、周波数及び時間情報の双方が同じ分解能であることである。例えば、高速フーリエ変換により、大量の周波数情報を有しない(低い周波数を制限した小さい窓)ことで時間分解能が増加し、又は低い周波数情報が増加して高周波の時間情報がスメアされる(立てば、直前の1024個のサンプルの間に、平均の実数の大きさR及び平均の虚数の大きさから成るxHzの実成分を有する)。本発明に係るシステムは、時間Tで、大きさM(及び位相P)を具えたxHzの正弦波を与え、周波数xHzでのxHzはフィルタの出力で、xはナイキスト周波数に一致する値である。2以上のフィルタ出力のストリームを有して、この出力を平行して進む他のフィルタ出力と比較することによって、信号の1又はそれ以上の成分の正弦波/非正弦波を判別することが可能となる。
データ抽出は、フィルタの出力の中に局所ピークを具えており、各局所ピークは時間及び周波数の双方の方向の局所ピークである。特定されるピークは、フィルタ出力の中の振幅のピークであり、これは隣接するフィルタのピークの振幅よりも大きい。フィルタの出力が三次元トポロジーとして表される場合、特定したピークはこのようなとトポロジーの頂上である。信号処理では、フィルタリングに続いて、フィルタの出力からデータが抽出される。各フィルタの出力は、特定の周波数(フィルタの中心周波数)で振幅を変える信号である。データの抽出方法は、信号処理の最終的な目的(圧縮又は解析等)に依存する。
従来技術のGB2034056では、フィルタの出力の中の局所ピークは特定されず、各フィルタの現状のエネルギレベルを単に表示するように、装置が動作する。本発明のシステムでは、フィルタの出力の中の局所ピークを特定する。例えば、フィルタfによって戻される値が92、94、92、94、96及び94である場合、2つの局所ピークがこのような一連の値の中に存在し、「94」及び「96」であるが、これらは単に時間方向の局所ピークである。本発明のシステムは、時間及び周波数の双方の方向で局所ピークであるこれらのピークを特定する。これは、時間方向における局所ピークのいかなる値も、それらが双方の次元における局所ピークである場合、フィルタfにそれぞれ隣接するフィルタに関して戻されるピークの値よりも大きくなければならないことを意味する。
効果的に、フィルタの出力が、時間、周波数及び振幅(又はエネルギ)の軸を具えた3次元のトポロジーに変換される。それぞれの個別のフィルタの出力は2次元曲線であり、全てのフィルタの出力は3次元面である。フィルタによって帰還した値で規定される面で、時間及び周波数の双方の領域における局所ピークは丘の頂上として表され、フィルタ出力の処理において特定されるポイントが存在する。重要なことは、非常に複雑な信号が比較的少数の特定のピークに減少することであるが、それにも係わらず、このようなピークは入力信号の中に含まれている重要な情報を含んでいる。
有益なことに、バンドパスフィルタは、カスケード接続された複数のバンドパスフィルタであり、フィルタのそれぞれの配列が同じ中心周波数及び1/√2よりも大きいQを有している。隣接するフィルタ及びカスケード接続したフィルタは、同じQの値を有する必要がない。
有益なことに、各カスケート段階の出力において局所ピークを検出する。フィルタによる曲線出力の振幅のピークから抽出されるデータにより、元の信号を構成する正弦波の再構成が可能となる。受信した複雑な信号から曲線によって規定される面のピークの振幅及び時間を規定する座標への一連の曲線への信号の変換が、効果的に発生する。
フィルタの出力からのデータ抽出は、代替的に、各フィルタによって発生する曲線の下側の領域であり、又は、波の中の絶対的又は相対的なレベルを通過する信号上のポイントである。信号のターニング・ポイントを抽出するデータとすることが可能である。フィルタの出力によって規定される面の中に存在する形態を規定するよう使用する突起部の発見といった方法を具えたトポロジーレベルでデータ抽出を実行できる。
有益なことに、本プロセスは、さらに、カーブフィッティング処理の実行を具えており、その出力は、抽出データの一部に少なくとも近似する曲線を規定する曲線データを具えている。カーブフィッティングは、元の信号のさらなる変換である。このようなケースでは、抽出データが面の局所ピークである場合、局所ピークを通過する(又はこれらの局所ピークに近似する)曲線を規定する曲線データを発生し得る。
例示を目的とする例として、信号の処理が10個のフィルタを使用する場合、時間とともに振幅を変える信号は各フィルタの出力であり、10個の別々の信号となる。共に見られるこれらの信号は面を規定する。抽出データはこれらの面のピークであり、カーブフィッティングが、規定した面においてピークにそれぞれ適合する分離した3次元曲線を規定する曲線データを出力する。
カーブフィッティングは、カーブフィッティングで規定される曲線が周波数領域を越えて進むように、実際には様々なフィルタ出力に由来するポイントを結合する。カーブフィッティングは、抽出データの部分を形成する複数のポイントを結合する面を事実上規定する。既知のフィルタのスタートアップ時及び終了時のエラーを補償するために、時間内でこれらの曲線のバックグラウンドを反復するのもまた有益である。
データ抽出処理の反転を実行して、(エラー修正又は圧縮といった)中間処理又は中間処理無しに、元の信号に少なくとも近似する(同一でもよい)信号を再構成可能である。再同期するために、ピーク間の距離の間続く正弦波を発生させてもよい。例えば、正のピークのみを検出する場合、正弦波は1つの全周期存続する。正及び負のピークを検出する場合、正弦波は1/2周期続く。様々な方法によって正弦波を発生させてもよい。マイクロプロセッサでは、正弦関数を呼び出すことによって、又は参照テーブルを使用することによって、これを実現する。波を開始/終了させることによるエラーが次の共振周期に伝わるように、共振器/発信器によって正弦波を有利に発生させてもよい。信号が周波数が変化する信号であると認められる場合、代わりに共振器/発信器を新たな周波数に再調整する。有益なことに、カーブフィッティング処理が局所ピーク上で実行されている場合、正弦波の周波数をピーク間の距離よりも小さな間隔で更新してもよい。
さらに、本処理は、信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更するステップを具えている。これは、様々な正弦波のより優れた識別(拡大)をサポートし、及び/又は静的な閾値レベルの比較器、例えば固定した精密/機械的装置が可能となる。これは、ある状況では、処理に電力の有利性を与える。
さらに、バンドパスフィルタの数を可変にすることが可能であり、本方法はさらに、信号をフィルタリングする際にバンドパスフィルタの数を変更するステップを具えている。複雑な信号が様々な周波数で多数の様々な成分を具えている。信号の中の成分の数の信号を処理するよう使用される数のフィルタを採用するのがより効果的である。比較的広い周波数帯域をカバーするようフィルタを設定することが可能であり、このような帯域の中で活動が検出される場合、複数の狭帯域フィルタを使用してこのような帯域の中の様々な成分をフィルタリングすることが可能である。信号を処理する際にこのような処理を連続的及びフレキシブルにすることが可能である。これは、ある状況では、処理に電力の有利性を与える。
ある実施例では、複数のバンドパスフィルタが、受信信号の成分毎にほぼ3の整数倍を具えており、処理がさらに、3つのバンドパスフィルタのセットを用いた受信信号の各成分の追跡を具えている。複雑な受信信号が中に20成分を有している場合、60のフィルタを使用して成分を追跡可能である。いずれかの側に(周波数の大きさの)広帯域フィルタを具えた狭帯域フィルタで信号の各成分をフィルタリング可能であり、成分毎に3つのフィルタを与える。成分が周波数を変える場合、変化の方向を検出する中央のフィルタのいずれかの側に広帯域のフィルタを具えた3つのフィルタのセットによって、これを追跡する。成分が周波数を変える際にフィルタが動作する周波数を自動的及び動的に調整する。これは、ある状況では、処理に電力の有利性を与える。
本処理に使用されるフィルタの数を、受信信号の中で活性なオクターブの数に基づいて識別可能である。複数のバンドパスフィルタは、受信信号の周波数のオクターブ毎に約1つ具えることができるが、より詳細な信号の解釈では、複数のバンドパスフィルタが、受信信号の周波数のオクターブ毎にほぼ12の整数倍を具えている。使用されるフィルタの数が増えると、それらの間の周波数の間隔が減少して、システムの分解能が大きくなる。
好適には、本処理がさらに、調波関係判別段階を実行するよう構成されており、当該処理の出力が、異なる正弦波間の調波関係を特定する調波関係データを具えている。抽出したデータは、受信信号の中に(音声信号のケースでは)音楽及び/又は音声の調波が有ることを示す。これを自動的に検出可能であり、調波のデータを生成及び記憶することで、例えば元の信号を再構成する際に必要とする抽出データの量を減らす。
本処理は、さらに、信号をフィルタリングする際に1又はそれ以上のフィルタのゲインを変更することを具えている。これの主要な利点は、例えば、固定した精密/機械式比較器といった静的な閾値レベルの比較器が可能なことである。
本システムの一実施例では、本処理が、さらに、抽出データの再構成の表示を具えている。このような再構成は、x軸の時間及びy軸の周波数及びz軸の振幅を具えたデータの3次元面の図とすることが可能である。
本発明に係る第3の態様によれば、信号を受信するステップと、信号の正弦波分離を実行するステップと、分離した正弦波からデータを抽出するステップと、カーブフィッティングの処理を実行するステップと、を具えた信号を処理するための方法であって、当該処理による出力が、抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えている処理するための方法が提供されている。
本発明に係る第4の態様によれば、信号を受信するための受信部と、信号の正弦波分離を実行するためのフィルタリング装置と、分離した正弦波からデータを抽出して、カーブフィッティングの処理を実行するよう構成された処理装置と、を具えた信号を処理するためのシステムであって、処理の出力が、抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えている処理するためのシステムが提供されている。
本発明に係るこれらの態様によれば、一連の振幅/時間/周波数曲線として信号の解析/処理が可能である。解析は、圧縮/認識を可能にする。処理は、再同期が後に続く周波数/振幅調整を有してもよい。
本発明に係る第5の態様によれば、信号を受信するステップと、少なくとも1のバンドパスフィルタで信号をフィルタリングするステップと、信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更するステップと、フィルタ又は各フィルタの出力からデータを抽出するステップと、を具える信号を処理するための方法が提供されている。
本発明に係る第6の態様によれば、信号を受信するための受信部と、信号をフィルタリングするための少なくとも1のバンドパスフィルタと、信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更して、フィルタ又は各フィルタの出力からデータを抽出するよう構成された処理装置と、を具える信号を処理するためのシステムが提供されている。
本発明のこれらの態様によれば、様々な正弦波のより良好な判別(拡大)、及び/又は静的な閾値レベルの比較器、例えば固定した精密/機械的装置が可能となる。
本発明係る第7の態様によれば、信号を受信するステップと、複数のバンドパスフィルタで信号をフィルタリングするステップと、信号をフィルタリングする際にバンドパスフィルタの数を変更するステップと、フィルタの出力からデータを抽出するステップと、を具える信号を処理するための方法が提供されている。
本発明係る第8の態様によれば、信号を受信するための受信部と、信号をフィルタリングするための複数のバンドパスフィルタと、信号をフィルタリングする際にバンドパスフィルタの数を変更して、フィルタの出力からデータを抽出するよう構成された処理装置と、を具える信号を処理するためのシステムが提供されている。
本発明のこれらの態様によれば、処理能力の最適化の改善が可能となる。
本発明に係る第9の態様によれば、信号を受信するステップと、少なくとも300のバンドパスフィルタで信号をフィルタリングするステップと、フィルタの出力からデータを抽出するステップと、を具える信号を処理するための方法が提供されている。
本発明に係る第10の態様によれば、信号を受信するための受信部と、信号をフィルタリングするための少なくとも300のバンドパスフィルタと、フィルタの出力からデータを抽出するよう構成された処理装置と、を具える信号を処理するためのシステムが提供されている。
誤解を避けるために言えば、本発明を、(携帯型又は独立型のいずれかの)独立型の専用装置、又は例えば、携帯型通信機器、PDA、パーソナルコンピュータ、ラップトップ、ゲーム機又は自動販売機等といった他の機能を具える大きな装置又はホスト装置/システムの中に具えて実施してもよい。このような装置、システム又は器具は、単一の集積回路を具えており、又は代替的に別の集積回路の別の構成部品で実施することによって異なる機能を与えてもよい。
当業者が多くの他のシステム、装置及び方法を本発明に組み込むよう有利に構成し得ることに留意されたい。
本書に記載した本発明の態様の多くの実施例は、電子部品、回路及び/又はブロック図に関連して図示及び説明されており、当業者は、プロセッサの制御の下で動作するアルゴリズム又はソフトウェアが、プロセッサ、アルゴリズム及び/又はソフトウェアが実施例の部分を形成する場合には、このような部品、回路及び/又はブロック図を有利にシミュレートするものと認識するであろう。本発明の実施例の任意の態様の部分を形成する1又はそれ以上の部品、回路及び/又はブロック図を、個別部品、組立部品又はシミュレートした部品を組み合わせて有利に形成してもよいことが、当業者によって理解されるであろう。組立部品は、例えば、集積回路プロセス及び/又はマイクロマシン・プロセスで作製された部品を具えてよい。シミュレートした部品を使用する1つの利点は、値が一定であるが又は可変であることを要しているかに係わらず、高精度の成分値の制御を与えることである。シミュレートした部品を使用する他の利点は、隣接するフィルタが互いに干渉せず、フィルタ出力が(四捨五入による誤差は別として)数学的に理想なことである。しかしながら、これは、高い耐性を有する部品が将来的に有用となる可能性があるため、シミュレートしていない実施例を除外するものではない。
任意の部品がシミュレートされる場合、及び/又はなんらかの処理がデジタル領域で実行される場合、当業者は、このような処理の前に、信号をデジタル形式で出力するか、又は1又はそれ以上のアナログ・デジタル変換器(添付図面には図示せず)を使用してデジタル化した又はサンプリング形式の信号を発生させる必要があることを認識するであろう。このようなサンプリングの速さは、解析される信号の最も高い周波数に依存する。
本書に記載の本発明に係る多くの実施例を、バンドパスフィルタ(BPF)及び/又はノッチフィルタ(NF)に関連して図示及び説明するが、このようなフィルタを図示及び/又説明する場合に、フィルタが高次のフィルタ、ローパスフィルタ(LPF)、ハイパスフィルタ(HPF)、及び/又はアンプといった能動部品を有利に具えていることを当業者は認識するであろう。好適な実施例におけるバンドパスフィルタは、有限インパルス応答フィルタ(IIR)を具えている。ある実施例は、事象発生器の入力部又は他のデータ簡約段階に与えられるBPF出力データを与える代替的な方法を有利に具えており、このような代替的な事象発生器の入力データは、例えば、高速フーリエ変換(FFT)処理が後に続く指数窓関数からのデータを具えていて、入力信号がこのような指数窓関数の入力部に与えられ、速いデータ演算速度という利点を有する。
1又はそれ以上の部品、回路、及び/又はブロック図の機能を、1又はそれ以上のプロセッサによって、順次処理、逐次処理及び/又は並列の組み合わせで実行してもよい。このようなプロセッサは、デジタル・シグナル・プロセッサ、並列プロセッサ、状態機械、マイクロプロセッサ、マイクロコントローラ、パーソナルコンピュータ、メインフレームコンピュータ等を具えている。
解析に要する信号及び本発明によって解析される信号が、任意のソースから出力され、このようなソースの例が、アナログ又はデジタル形式で予め取り込まれ又は発生し及び記憶された信号を含んでおり、センサ、トランスデューサ及び/又はシナリオシミュレータからリアルタイムで受信した信号を含んでいることを当業者は認識するであろう。このような信号発生源の特定の例として:ハイドロフォン、マイクロフォン、無線通信機又はラジオ受信機、TV又はビデオ機器、レーダー機器、電波望遠鏡、(MRI,CAT,IR,超音波及びX線といった)医用画像機器、知的医学解析システム、医用又は科学シミュレーション機器を有している。
本発明のさらなる態様及び利点が、所定の実施例のみによって添付図面を参照して本発明に係る好適な実施例に関する以下の説明から明らかとなろう。
様々な異なる信号処理法を図1から図4に示す。これらの方法を実施するためのシステム及び装置の例を図5に示す。一般に、信号を処理するための装置は、信号を受信するための受信機と、信号をフィルタリングするための1又はそれ以上のフィルタと、データ抽出/簡約を処理を実行するための処理装置とを具えている。この処理装置は、フィルタの動作をうまく動的に制御する。
本発明の態様による装置は、1又はそれ以上の処理段階の組み合わせを具えている。この処理段階は、正弦波分離段階と、データ抽出/簡約段階と、カーブフィッティング段階と、調波関係判別段階と、データ処理及び操作段階と、及び/又は逆変換段階とを有している。
正弦波分離段階は、1又はそれ以上のBPF及び/又はフーリエ変換を具えている。このような正弦波分離段階は、入力信号に関する振幅、周波数、時間及び/又はQ情報を出力する。当業者は、大文字Qを使用してフィルタ又は同調回路の線質係数として知られている特性を指定することを認識するであろう。
データ抽出/簡約段階を使用して分離した正弦波から、サイン波の振幅ピーク、又はこの波によって規定される領域といったデータを抽出する。データ簡約は、正弦波分離段階からデータ出力の量を減らし、このようなデータ簡約は特定の適用に適しており、特定の適用に関する十分な情報を保持しつつデータ簡約を最大源にする。
信号解析器の第1の実施例を図1のフローチャートを参照して説明する。このフローチャートは、装置がどのように入力信号を変換するのかを示す。ステップS1001で入力信号が装置に与えられる。1又はそれ以上のBPF及び/又はフーリエ変換を具えるステップS1002の正弦波分離段階で、装置が、入力信号の中の正弦波に関する情報を判別して、入力信号に関する振幅、周波数、時間及び/又はQ情報を含む情報を出力する。
このような情報全ては、ステップS1003のデータ抽出/簡約段階に供給与えられる(本例のデータ簡約法は、以下に詳細に説明されている)。信号解析法の最も簡単な形式では、データ簡約が行われない。この段階S1003で、フィルタ出力からのデータ抽出が、フィルタ出力の局所ピークの特定を具えており、各局所ピークは時間及び周波数方向の双方における局所ピークである。そして、このようなデータ(局所ピーク)は、信号解析法が使用する適用に応じて処理される。
フィルタ出力から他のデータを抽出できる。例えば、このようなデータは各フィルタによって2つのゼロ交差点間で規定される曲線の下の領域とすることができる。抽出が行われる前に各フィルタの出力を処理して、例えばフィルタ出力の曲線が正弦波であるかどうかを判別できる。そして、これらの曲線が正弦波である場合にのみ、データ抽出段階で使用される。
図1の実施例では、ステップS1003からデータ出力がステップS1004に与えられる。ステップS1004のカーブフィッティング段階で、装置が局所ピークを結合させることによって形成されるラインの特性を判別するが、各局所ピークは時間方向及び周波数方向の双方における局所ピークである。当該ラインの特性のこのような判別は、ライン全体又は部分の曲率及び/又は形状特性を記述する係数の判別を有してよい。ステップS1004からの情報出力がステップS1005及びステップS1006の一方又は双方に与えられる。
装置は、ステップS1005で調波関係判別段階を有してもよく、ここでは、入力信号の中の分離した正弦波に関する情報を使用して1又はそれ以上の調波系の存在を特定する。各調波系は、例えば、基本周波数、2次の調波、3次の調波、4次の調波、5次の調波及び他の調波のうちの1又はそれ以上を具える。ステップS1005からの出力が、存在する場合、ステップS1006に与えられる。
ステップS1006で、ステップS1004及び/又はステップS1005において判別された情報が、特定の適用の要請に従って装置によって使用される。このような情報の使用は、例えば、記憶、さらなる処理、さらなる解析、信号の再構成、及び/又は例えば、画像表示装置及び/又は音声出力装置又はトランスデューサ/アクチュエータへの出力を含んでいる。本装置は、正弦波分離段階が関連する線質係数(Q)特性を具える1又はそれ以上のBPFを具える場合、動的なQ制御を有してもよい。本装置は、1又はそれ以上のBPFの動的な中心周波数制御を有してもよい。
このような動的なQ及び/又は動的な中心周波数制御は、ステップS1007で表されており、その制御はステップS1002からS1005のうちの1又はそれ以上から入力データを抽出し、ステップS1007の出力がステップS1002に与えられ、ステップS1007からの出力情報を使用して、正弦波分離段階の中でQ及び/又は1又はそれ以上のBPFの中心周波数を制御する(このような動的なQ及び中心周波数制御を以下に詳細に説明する)。
信号解析器の第2の実施例を図2のフローチャートを参照して説明する。ステップS2001で入力信号が装置に与えられる。複数のBPFを具えるステップS2002の正弦波分離段階で、装置が、入力信号の中の正弦波に関する情報を判別して、入力信号に関する振幅、周波数、時間及び/又はQ情報を含む情報を出力する。
ステップS2002で、各BPFが、IIRフィルタ、又は有限インパルス応答(FIR)フィルタ又はIIR及びFIRフィルタの組み合わせを具えている。このようなFIRフィルタは、後にFFTが続く指数窓関数の使用によって実現されるIIRインパルス応答を具えている。ステップS2002からのこのような情報全ては、ステップS2003のデータ簡約段階に与えられる(データ簡約法の例を以下に詳細に説明する)。ステップS2003からのデータ出力が、ステップS2004,S2005,S2006のうちの1又はそれ以上に供給される。
本装置は、ステップS2004のカーブフィッティング段階を有しており、本装置が局所ピークを結合することによって形成されるラインの特性を判別するが、各局所ピークは時間及び周波数方向の双方における局所ピークである。ラインの特性のこのような判別は、曲率及び/又はライン全体又は部分の形状特性を記述する係数の判別を有する。ステップS2004からの情報の出力は、存在する場合には、ステップS2005及びステップS2006のうちの一方又は双方に与えられる。
また、本装置は、ステップS2005で調波関係判別段階を有しており、ここでは入力信号の中の分離した正弦波に関する情報を使用して1又はそれ以上の調波系の存在を特定する。各調波系は、例えば、基本周波数、2次の調波、3次の調波、4次の調波、5次の調波及び他の調波のうちの1又はそれ以上を具える。ステップS2005からの情報の出力は、存在する場合には、ステップS2006に与えられる。
ステップS2006で、ステップS2003及び/又はステップS2004及び/又はステップS2005において判別された情報が、特定の適用の要請に従って装置によって使用される。このような情報の使用は、例えば、記憶、さらなる処理、さらなる解析、信号の再構成、及び/又は例えば、画像表示装置及び/又は音声出力装置及び又はトランスデューサ/アクチュエータへの出力を含んでいる。
本装置は、正弦波分離段階が関連する線質係数(Q)特性を具える1又はそれ以上のBPFを具える場合、動的なQ制御を有してもよい。本装置は、1又はそれ以上のBPFの動的な中心周波数制御を有している。このような動的なQ及び/又は動的な中心周波数制御は、ステップS2007で表されており、ステップS2002からS2005のうちの1又はそれ以上から入力データを抽出し、ステップS2007の出力がステップS2002に与えられ、ステップS2007からの出力情報を使用して、正弦波分離段階の中でQ及び/又は1又はそれ以上のBPFの中心周波数を制御する(このような動的なQ及び中心周波数制御を以下に詳細に説明する)。
信号解析器の第3の実施例を図3のフローチャートを参照して説明する。ステップS3001で入力信号が装置に与えられる。単一のIIR BPFを具えるステップS3002の正弦波分離段階で、本装置が、入力信号の中の正弦波に関する情報を判別し、入力信号に関する振幅、周波数、時間及びQ情報を含む情報を出力する。ステップS3002からのこれらの情報が、ステップS3003のデータ簡約段階に与えられる(データ簡約法の例を以下に詳細に説明する)。
ステップS3003からのデータ出力が、ステップS3004,S3005又はS3006のうちの1又はそれ以上に与えられる。本装置は、ステップS3004のカーブフィッティング段階を有しており、本装置が局所ピークを結合することによって形成されるラインの特性を判別するが、各局所ピークは時間及び周波数方向の双方における局所ピークである。このようなラインの特性の判別は、ライン全体又は部分の曲率及び/又は形状特性を記述する係数の判別を有する。ステップS3004からの情報の出力は、存在する場合には、ステップS3006に与えられる。
また、本装置は、振幅及び/又は位相情報の解析を実行するステップS3005の処理段階を有している。ステップS3005からの出力が、存在する場合には、ステップS3006に与えられる。ステップS3006で、ステップS3003及び/又はステップS3004及び/又はステップS3005において判別された情報が、特定の適用の要請に従って本装置によって使用される。このような情報の使用は、例えば、記憶、さらなる処理、さらなる解析、信号の再構成、及び/又は例えば、画像表示装置及び/又は音声出力装置及び又はトランスデューサ/アクチュエータ装置への出力を含んでいる。
本装置は、動的なQ制御を有している。このような動的なQ制御は、ステップS3007で表されており、ステップS3002からS3005のうちの1又はそれ以上から入力データを取得し、ステップS3007の出力がステップS3002に与えられ、ステップS3007からの出力情報を使用して、正弦波分離段階の中でBPFのQを制御する。
本装置は、1又はそれ以上のBPFの動的な中心周波数制御を有している。このような動的な中心周波数の制御は、ステップS3008で示されており、ステップ3002からステップS3005のうちの1又それ以上から入力データを取得し、ステップS3008からの出力情報を使用して正弦波分離段階の中でBPFの中心周波数を制御する(このような動的な中心周波数制御を以下に詳細に説明する)。
信号解析器の第4の実施例を図4のフローチャートを参照して説明する。ステップS4001で入力信号が装置に供給される。ステップS4001で、入力信号が本装置に与えられる。単一のBPFを具えるステップS4002の正弦波分離段階で、本装置が、入力信号の中の正弦波に関する情報を判別し、入力信号に関する振幅、周波数、時間及び/又はQ情報を含む情報を出力する。
ステップS4002で、BPFがIIRインパルス応答(後にFFTが続く指数窓関数を使用することによってこれを実現してもよい)を具えた有限インパルス応答(FIR)フィルタを具えている。ステップS4002からのこのような情報全ては、ステップS4003のデータ簡約段階に与えられる(データ簡約法の例を以下に詳細に説明する)。
テップS4003からのデータ出力が、ステップS4004,S4005又はS4006のうちの1又はそれ以上に供給される。本装置は、ステップS4004のカーブフィッティング段階を任意に有しており、本装置が局所振幅ピークを結合することによって形成されるラインの特性を判別するが、各局所ピークは時間及び周波数方向の双方における局所ピークである。このようなラインの特性の判別は、ライン全体又は部分の曲率及び/又は形状特性を記述する係数の判別を有する。
ステップS4004からの情報の出力は、存在する場合には、ステップS4006に与えられる。本装置は、振幅及び又は位相情報の解析を実行するステップS4005の処理段階を任意に有している。ステップS4005からの出力は、存在する場合には、ステップS4006に供給される。ステップS4006で、ステップS4003及び/又はステップS4004及び/又はステップS4005において判別された情報が、特定の適用の要請に従って本装置によって使用される。このような情報の使用は、例えば、記憶、さらなる処理、さらなる解析、信号の再構成、及び/又は例えば、画像表示装置及び/又は音声出力装置及び又は他のトランスデューサ又はアクチュエータ装置を含んでいる。
本装置は、動的なQ制御を任意に有してよい。本装置は、1又はそれ以上のBPFの動的な中心周波数制御を有してよい。このような動的なQ及び/又は動的な中心周波数の制御は、ステップS4007で表されており、ステップS4002からS4005のうちの1又はそれ以上から入力データを取得し、ステップS4007の出力がステップS4002に与えられ、ステップS4007からの出力情報を使用して正弦波分離段階の中でQ及び/又はBPFの中心周波数を制御する(このような動的なQ及び中心周波数の制御を以下に詳細に説明する)。
図5は、信号を処理するためのシステムの好適な実施例を示す。本システムは、信号400を受信するための受信部を有している。さらに、本システムは、正弦波分離段階451と、データ抽出/簡約段階452と、信号事象解析器413とを有している。正弦波分離段階451は、複数のバンドパスフィルタ(BPF)を具えている。データ簡約段階452は、複数の事象発生器を具えている。信号事象解析器413は、カーブフィッティング段階、調波関係判別段階、データ処理及び操作段階及び/又は逆変換段階のうちの1又はそれ以上を実行できる。一連の事象発生器421及び信号事象解析器413は、信号処理器として実行し得る。
入力信号400が、各BPF(401−406)の入力部に与えられ、各BPFの出力がそれぞれの事象発生器(421−426)の入力部に与えられる。より高次のバンドパスフィルタの配列を使用でき、例えば、それぞれが直列にカスケード接続された2及び3のフィルタを具えた4次及び6次のフィルタ配列を使用できる。このようなケースでは、フィルタリングの中間段階でデータ抽出を実行でき、抽出が実行されて第1のフィルタ(等)の出力及び最終出力となる。
各事象発生器の出力が信号事象解析器413の入力部に与えられる。信号事象解析器からのデータ出力414は、入力信号400の記述及び解析を与えており、事象情報に関する曲線の係数及び/又は情報を具えている。このような曲線の係数は曲線を記述しており、各曲線がある最大振幅事象から別の最大振幅事象に移動するラインとして定義される。通信インターフェイス450は、より大きな装置又はシステム(図示せず)とやり取りを行うデータを制御し易くする。
ここで、単一のBPF、単一の事象発生器及び単一の事象解析器を示す図6を参照して、振幅が変化する単一の正弦波を具える信号の例を示す図7を参照し、事象ストリームの例を示す図8を参照し、事象ストリームの代替的な例及び曲線の例を示す図10を参照し、どのように動的なQ及び/又は中心周波数の制御を使用するのかを示す図11を参照して、図5の装置の動作説明をすることとする。
図7は、振幅が変化する1つの正弦波である信号100を示す。この信号100は、図5のBPFのうちの1つからの出力である。信号100を事象ストリームによって正確に表し得ることが判別され、各事象が信号の中のピークに関する情報を具えている。このような事象情報は2つの部分を具えており、第1の部分は、例えばエネルギ、電流又は電圧を表すピーク振幅に比例する値を具え、第2の部分はピークが生じる時間の記録を具えている。事象発生器204は、このような事象情報がBPF209からのデータ出力よりもあまりデータを具えていないため、データ簡約を実行する。
矢印104,105及び106によって表される事象は、事象の例であり、信号100の事象ストリームの部分を形成するが、ピーク101,102及び103でそれぞれ発生する。代替的に、事象ストリームが、負のピーク107及び108でそれぞれ生じる点線矢印109及び110で示す負のピーク事象を具えてもよい。さらなる実施例では、事象ストリームが、正及び負のピークの双方に関する事象を具えてもよい。
図6は、1/√2のQファクタを有するBPF209と、事象発生器204と、信号事象解析器220とを具える装置の実施例を示す。BPF209は、図5のフィルタ401−406のうちの1を具えている。事象発生器204は、図5の事象発生器421−426のうちの1を具えている。信号事象解析器220は、図5の信号事象解析器413の一部又は全部を具えている。
入力信号206は、バンドパスフィルタ209への入力部であるノード200に接続されており、バンドパスフィルタはIIRフィルタを具えている。バンドパスフィルタの出力である信号207が、事象発生器204の入力部に与えられる。事象発生器204は、図7に示す正の1つのピークに関する事象が発生するように、アルゴリズムを動かす回路及び/又は処理装置を具えている。図6の事象発生器204の出力を信号ノード205として示しており、ノード207で信号が存在する間、連続的に事象ストリーム208を与える。このような事象ストリームの部分を矢印104,105及び106として図8に示す(104,105及び106は図7で同様にラベルされた事象を示す)。図6を再び参照すると、事象発生器204が、例えば値の変化、値の差と組み合わせて様々な方法でピークを認識する。変化値法は、例えば、サンプル値が減少する場合にピークを認識するように、1又はそれ以上の前の振幅のサンプル値を記憶することによって機能する。微分値法は、例えば、微分した信号がゼロに達する場合に事象が発生するように、信号207を微分する回路及び/又は処理装置を実行するアルゴリズムを比較することによって機能する。
入力信号206は、バンドパスフィルタの中心周波数又はこれに近い周波数の1つの正弦波を具え、又は2以上の正弦波であって、それらのうちの一方がバンドパスフィルタの中心周波数又はこれに近い周波数を具えている。
事象208は、事象が発生して信号事象解析器220によって受信される時点でのノード205の出力である。信号事象解析器は事象情報を使用して、バンドパスフィルタの中心周波数又はその近くでの入力信号における正弦波の存在を、単一の事象(又は半周期)に対して非常に高い時間分解能で判別する。このような細かい時間分解能の判別は、応答時間といった影響が補償されるように、信号事象解析器の中のフィルタQ(線質係数)及び/又はフィルタの次数に関するデータ及び/又フィルタの特性に関するデータの値の使用を必要とする。信号事象解析器は、入力信号200に関連する出力データ221を出力する。
信号事象解析器220は、事象が発生する際に事象ストリーム情報を使用し、又は代替的に信号事象解析器は後で使用するための事象情報を記憶する。
ある実施例では、事象発生器204の中に記憶手段を含めることに利点があり、このようなケースでは、事象発生器がフィルタの出力サンプルデータ、ノード207を記憶し、事象情報を取得するためにこのようなデータを後で処理する。
好適な実施例では、信号事象解析器220が、1又はそれ以上のアルゴリズムを有しておりラインの部分の曲率を判別するが、このようなラインは、このようなアルゴリズムがカーブフィッティング段階を形成する事象の上部を結合することによって形成される。事象607で開始して事象608で終わる事象ストリームが図10に示す時間の入力信号(図6の200)を表す例を図10に示す。入力信号は、曲線604,605,及び606によって表されている。曲線604は、このケースでは直線であり、事象振幅600で開始して事象振幅601で終了する。曲線605は事象振幅601で開始して事象振幅602で終了する。曲線606は事象振幅602で開始して事象振幅603で終了する。各ラインの曲率を3つの係数で定義し、この3つの係数が事象情報から導かれる。
当業者は、曲線の係数のこのような導出を、例えば既知のカーブフィティング・アルゴリズムで取得してもよいことを認識するであろう。入力信号のバリエーションの多くの例で、ある曲線による係数が後に続く曲線の1又はそれ以上の係数を形成する。代替的に、ある実施例では、例えば、図10のポイント600からポイント603までの曲線全体といった、より複雑な曲線を記述する係数を導き出すのが有益である。入力信号を一連の1又はそれ以上の曲線を表す一連の係数で記述することで、信号を解析且つ処理する新たな方法を提供する。
図6を再び参照すると、ある実施例では、(図1から図4のステップS1007,S2007,S3007,S4007でそれぞれ示し、図5の1又はそれ以上のフィルタ401から406のQを制御するよう使用する)動的なQ制御を含めるのが有益である。このような動的なQ制御は、ノード205での事象振幅及び/又はノード207でのフィルタ出力の振幅に関するフィルタ209のQ制御を与える。フィルタQのこのようなフィードバック制御は、例えば、ノード205とやり取りする信号事象解析器220といった他のソースからの信号によって付加的又は代替的に影響を及ぼされる。図11は、Qフィードバック制御250(符号200,204,205,206,207,209,220及び221を付けた項目は図6において同様に符号を付けた項目と同じ機能を有する)を介してフィルタ209のQを制御するよう使用される信号251,252,253及び254の例を示す。
ノード207及び/又はノード205での信号の振幅の増加に関連して増加するようフィルタQが制御される場合には、ノード207及び/又はノード205での信号振幅の制限の制御を要して動揺/不安定性を防止しており、このような制限は、例えば、フィードバック信号の付加的な制御ループ、フィードバック信号のスルーレート制御、事象発生器204の中のゲイン/減衰の制御、及び/又はエネルギ測定値及び/又は1又はそれ以上の前のサンプルの信号レベルに応じたフィードバック信号、の組み合わせによって実現される。
ある実施例では、ノード205での事象振幅及び/又はノード207でのフィルタ出力振幅に対する事象発生器204の中のゲイン/減衰を制御するのが有益である。
実施例は、フィルタQフィードバック制御を有しており、入力信号の周波数弁別を増加させる。このような弁別は、フィルタの中心周波数と異なる周波数の入力信号がフィルタQが増加するにつれてより大きく減衰する(フィルタの応答、例えば、図9の符号503がより狭くなる)ために生じ、フィルタの中心周波数と等しい周波数の入力信号はあまり減衰しない(フィルタ応答曲線が狭くなるとフィルタ応答曲線の中心周波数での振幅が増大する)。
実施例は、フィルタQ及び/又はゲイン/減衰フィードバック制御を有しており、事象発生器204及び/又は信号事象解析器220におけるデータ処理の量を減らし、例えば、データの値が所定の範囲の中に保持され整数演算を使用して処理を実行可能である。
実施例は、フィルタQ及び/又はゲイン/減衰フィードバック制御を有しており、フィルタ209の応答時間を増やす。例えば、フィルタQが低いレベルで開始して信号を検出すると増加する場合、このような応答時間を増やしてよい。
ここで再び図6を参照すると、ある実施例では、(図1から図4のステップS1007,S2007,S3008,S4007でそれぞれ示し、図5の1又はそれ以上のフィルタ401から406の中心周波数を制御するよう使用する)動的な中心周波数制御を有するのが有益である。このような動的な中心周波数制御は、ノード205での事象振幅及び/又はノード207でのフィルタ出力振幅に対するフィルタ209の中心周波数の制御を与える。フィルタ中心周波数のこのようなフィードバック制御は、他のソース、例えば、ノード205とやり取りする信号事象解析器220からの信号によって付加的又は代替的に影響を及ぼされる。図11は、中心周波数フィードバック制御250(符号200,204,205,206,207,209,220及び221を付けた項目は図6において同様に符号を付けた項目と同じ機能を有する)を介してフィルタ209の中心周波数を制御するよう使用される信号251,252,253及び254の例を示す。
実施例は、フィルタ中心周波数フィードバック制御を有してもよく、入力信号の周波数が最新のBPFの中心周波数に近いが等しくない場合、例えば入力信号を追尾又は追跡する。
図6及び/又は図11の装置の適用は、ある信号の存在の非常に速い検出の要求及び/又はある信号への変化の迅速な検出の要求を含めてもよい。一例は、信号全体の強度の変化によるゆっくりとした振幅の変化を迅速な変調の変化と区別する必要がある場合、振幅の復調である。
第2の例は、変調のエッジが瞬時の位相反転を具える位相の復調であり、このような信号では、既知の周波数で発生する事象ストリームが半周期だけ一時的に突然にシフトするが、このような時間シフトは第1又は第2のシフトした事象で検出される。第3の例は、副搬送波変調の小さなバーストの検出であり、このようなバーストは、数周期、例えば4周期の副搬送波変調を具えている。第4の例は、位相ロックループの迅速な安定化を要する場合である。
再び図5を参照すると、(解析するのを要する)入力信号400が少なくとも1の正弦波信号を具えており、少なくとも1の正弦波信号が周波数の範囲内にある。図6に関して本書に記載の任意のバリエーションを図5の装置に含めてもよい。
各BPF401−406におけるフィルタの中心周波数は、あるBPFから次のBPFに向かって増加し、このような方式では、解析すべき周波数範囲がカバーされる。中心周波数及び各BPFのフィルタの帯域幅は、隣接するBPFの応答曲線が重ならないようになっていて、この一例を図9に示し、図9では中心周波数500−502をそれぞれ具えるフィルタの応答503−505がBPF402−404に対応する(図5)。明りょうのために、3つのフィルタ応答曲線のみを図9に示す。
図5を再び参照すると、各事象発生器が事象ストリームを出力し、事象発生器421−426からの事象ストリーム407−412がそれぞれ信号事象解析器413に入力される。
6つのBPFを図5に示すが、特定の装置で解析するのに要する入力信号の既知の特性に応じて特定の装置において任意の数を(各フィルタから事象を取得するための手段とともに)有利に含めてもよい。例えば、適用は、2個、10個、100個又は1000個のBPSを含めることによって利益を得て、例えば音声解析の分野において、オクターブ当たり300よりも多いBPF以上を含めることは有益であり、中心周波数が対数的に増加する。直線的又は代替的に非線形的に中心周波数が増加する場合、他のアプリケーションがBPFを要する。
信号事象解析器413は、出力データ414及びインターフェイス450を介してより大きなシステム(図示せず)の部分と通信するが、このようなより大きなシステムは、例えば、メモリ、追加のプロセッサ、入力インターフェイス、出力インターフェイス、ディスプレイ、アクチュエータ、有線又は無線通信機器のうちの1又はそれ以上を具えている。
データ簡約段階452の出力は、事象ストリームを具えており入力信号を記述するが、解析を用いて元の信号を再構成又は同期できる。信号事象解析器は、入力信号の中の各正弦波の存在を示す情報を出力する。複数の事象ストリームの使用により、1つの事象に関して非常に高い時間分解能での入力信号の周波数成分の解析がし易くなる。
信号事象解析器は、特定の装置の信号解析の要求に従って入力信号400が解析及び/又は処理される方法で、少なくとも1の事象発生器からの事象ストリームの事象からの事象情報を処理するが、このような解析及び/又は処理は、既知又は予期される入力信号400の特性に従っている。このような特性は、1つの信号のタイプから別の信号のタイプに変化し、下記の信号のタイプを説明する(以下の例では、BPF401−405の中心周波数をそれぞれf−fにする)。
第1の例では、入力信号400が周波数fで1つの正弦波を具えている場合、事象発生器423が周波数fで事象を出力する。また、事象発生器422及び424は、図9を参照すると応答曲線503および505が周波数f501がポイント506で示すバンドパスフィルタ503及び505の中にあることを示すため、(振幅が小さい)周波数fで事象を出力する。信号事象解析器413は、各事象ストリーム(各事象発生器)からの事象率を判別し、信号事象解析器が各事象発生器に関する各BPFの中心周波数を識別するため、信号事象解析器は事象発生器422及び424からの事象率がBPF402及び404の対応する中心周波数とは異なるものと判別する。
信号事象解析器は、全ての事象率及び振幅情報を使用して入力信号の解析を提供し得る。このような解析は、事象率が事象発生器のフィルタ処理した中心周波数と異なる場合に事象ストリーム情報の廃棄を有しており、減少した振幅情報とともにこのような異なる情報を使用する。このような例では、fの事象ストリーム409が入力信号の解析を記述又は提供する。
入力信号の第2の例では、入力信号400がほぼ等しい振幅の周波数f及びfの2つの正弦波を具えている場合、事象発生器422は、周波数fとfとの関係に従って事象の振幅がどれほど変化しても(この振幅の変化は、2つの近接周波数のうなり周波数に関する)周波数fで事象を出力する。事象発生器423からの事象をさらに解析処理の部分として使用してもよい。
入力信号の第3の例では、入力信号400がインパルスを具えている場合、ピークの出力レベルが全ての事象ストリームで生じ、事象ストリームの振幅が周期毎に同じ率で減衰する。
ある適用は、例えば、事象率がフィルタに対応する率と異なる場合にこれらの事象が抑制又は削除されることで信号事象解析器に出力されないように、各事象発生器の中の追加の情報によって利益を得る。
1又はそれ以上の処理装置がBPF及び/又は事象発生器をシミュレートする場合、例えば、BPFの群(このような群の前に偽信号防止のフィルタを必要とする)のサンプルレートを減らすことによって、処理の量を減らすことが好都合である。
1又はそれ以上の処理装置がBPF及び/又は事象発生器をシミュレートする場合、例えば、プレフィルタによって正弦波を検出する場合、BPF及び事象発生器の群の処理を実行するだけで処理の量を減らすのが好都合である(プレフィルタの帯域幅がBPF全体の群の周波数範囲をカバーする場合)。
ある実施例では、3次元(3D)のマップ又はランドスケープ(landscape)として信号/事象情報を表すことによって、信号解析を改善してもよい。3次元のグラフが多くの場合、2次元のx/y面から延びる高さを表現する追加の次元zを具えており、従来の2次元のx/yグラフと関連して表現又は記載されるものと当業者は認識するであろう。事象情報のこのような3次元マップの例は、x軸が時間を表し、y軸が周波数を表し、z軸が振幅を表すものである。事象情報を、少なくとも2つの事象発生器から又は少なくとも2つのフィルタから直接所得したデータから所得してもよい。全ての事象をプロットするか、又は事象データ全体の一部をプロットしてもよく、例えば事象は、隣接する事象に関連する所定の振幅を具えている。
一例として、一定の周波数及び振幅の単一の正弦波を具える入力信号を考える。この例では、3次元マップが、z方向に一定の振幅を有するx方向に向いた水平な直線を具えている。
さらに、このような3次元マッピングを解析の目的で処理してもよく、又は例えば視覚化解析のために表示及び/又はプリントしてもよい。
好適な実施例では、カーブフィッティング手段への入力データとして少なくとも1の事象ストリームからの事象情報を使用することによって、追加の信号解析及び/又は信号データ圧縮を実現する。このようなカーブフィッティングを、振幅及び時間からなる2次元の図10を参照して本書で説明するが、このようなカーブフィッティングを(第3の次元として周波数を追加して)都合良く3次元に拡張してもよく、このようなケースでは、ピーク事象の振幅値を結合することによってそれぞれの3次元曲線を形成し(ピーク事象はその近くよりも高い事象値である)、ピークの振幅は事象ストリームに対して近接又は隣接している。このようなラインは3次元空間を通る(本書では3次元マッピングとして説明する)。
入力信号は、1又はそれ以上のこのような3次元ラインを具えてもよい。図15は、事象1503の上部と1507の上部とを結合する曲線を形成するこのような3次元のラインの例を示す(図は、3次元のうちの2次元を表しており、3次元が時間、振幅及び周波数であり軸1500、1502及び1503によってそれぞれ表され、実際には互いに直交する)。第1の事象1503は時間t及び周波数tで発生し、第2の事象1504はより小さい振幅(t,f)で、第3の事象1505はさらに小さい振幅(t,f)で、第4の事象は第3の事象よりも大きい振幅であり、第5の事象1507(t,f)はさらに大きい振幅である。明りょうのために、曲線を形成するのに使用しない事象ストリームからの事象を図示していない。
このようなカーブフィッティングは、長時間にわたる入力信号の挙動を記述する1又はそれ以上の組の曲線係数を与える。このような曲線係数は、入力信号のデータの圧縮表現を与える。
このようなカーブフィッティングの例は、入力信号がある周波数の正弦波を具えており、周波数が時間とともに増加して振幅が時間とともに減少するように、正弦波が変化する。このような例では3次元のラインが、入力信号が周波数及び振幅を変化させるまで、z方向に一定の振幅を有してx方向に向いた水平な直線を具えている。入力信号が変化する場合、ラインの周波数(y軸)が斜めに上昇し、振幅(z軸)が下降する。
一般に、カーブフィッティングは、1次、2次、3次又はそれ以上の次元で実行してよい。多次元の例は、図5の複数の装置を具える装置を具えており、図5のそれぞれの装置は別々の入力信号を有しており、図5の複数の装置は配列を形成し、複数のセンサ又はセンサ配列から別々の入力信号を所得する。このような別々の入力信号をベクター又は配列で形成又は記述してもよい。このようなセンサの配列の例は、ビデオカメラセンサ、目(目を模したもの)、多チャンネル録音(バイノーラル録音を含む)、ステレオヒアリングのシミュレーションを含んでいる。図5の複数の任意の装置の中の様々な処理段階からのデータを組み合わせ、出力データを形成するよう処理及び/又は比較してもよい。このような出力データを、本書に記載の使用又は適用のために使用してもよいが、これに限定されない。さらに、図5の装置の中の様々な処理段階からのデータを、図5の他の装置の任意の処理段階の入力部に与えてもよい。当業者は、本発明に係る装置の代替的な多次元の実施例から利益を得ることができる他の適用が有るものと認識するであろう。
調波関係の判別段階は、入力信号の中の別々になった正弦波に関する情報を使用して、1又はそれ以上の調波系の存在を特定するアルゴリズムを具えている。各調波系は、基本周波数を、2次の調波、3次の調波、4次の調波、5次の調波及び他の調波のうちの1又はそれ以上を具えてもよい。図12を参照すると、このような特定は、特定の時間、例えば時間1200で生じてもよく、例としての調波系は、周波数1210の基本周波数1206と、周波数1211の2次の調波1207と、周波数1212の3次の調波1208と、周波数1213の4次の調波1209とを具えている(明りょうのため、より高次の調波を図示していない)。同じ調波系の特定は、後の時間(例えば、時間1201から1204)で発生してもよく、このようなケースでは、このような後の特定が例えば声の高さの上昇を表示してもよい。
調波関係を判別する信号解析の多くの適用が有り、その例として:音声認識、生体認証解析、楽器の特定、化学系統の認識、データ圧縮、消音を含んでいる。
ある実施例では、追加のデータ処理及び操作段階を含めるのが都合がよく、このような追加段階への入力データは、データ簡約段階、カーブフィッティング段階、調波関係判別段階のうちの1又はそれ以上(本書に記載の実施例の組み合わせによるこのような段階から)から取得される。このような追加の処理段階は、特定の適用のために必要とされ、例えば、これを使用して逆変換段階に情報を提供する。このような逆変換段階は、正弦波を形成して出力信号(このような出力信号は振幅、周波数及び位相情報を具えている)を形成するようそれらを追加するパラレルシンセサイザの形式をとってもよい。データ処理及び操作段階及び逆変換段階を用いた例としての適用は、楽器の音の高さの移動、発声の高さの移動、信号範囲全体の高さの移動、消音、(部屋の周波数応答を解析しその後で演奏している音声の特定の周波数の振幅を変えて補償することによる)室内音響補償を有している。
ある実施例では、代替的なデータ簡約段階を使用して事象を発生させるのが有益である。代替的な例は、単一の一定閾値及び複数の一定閾値を有している。
単一の一定閾値の代替的なデータ簡約段階を図13を参照して説明し、関連するフィルタ周波数に関する正弦波のyes/noの存在情報を与える。事象1303が、例としての正弦波1300が閾値1301を越えて上昇する時間1302で発生する。そのようにして発生した事象は、信号1300の存在を示している。ある適用は、信号が閾値を下回る場合に事象の発生の利益を得てもよく、このような事象の例を矢印付き点線1305で示す。
複数の一定閾値の代替的なデータ簡約段階を図14を参照して説明する。複数の一定閾値の方法は、単一の一定閾値の方法と同じような方法で機能するが、2又はそれ以上の閾値トリガ回路を使用する。3つの例としての閾値を、符号1401、1402及び1403で示しており、(時間1404で)事象1405、(時間1406で)1407及び(時間1408で)1409をそれぞれ発生させる。特定の閾値に関して発生した事象は、閾値及び/又は事象が発生した時間に関する情報を具えている。複数の一定閾値の方法では、より多くの閾値レベルを使用するため、事象の間が縮まってデータ簡約があまり発生せず、極限ではレベルが事実上連続的になってデータ簡約が生じないようになる。
当業者は、閾値のトリガに履歴を含めることが有益であることを認識するであろう。
アクティブ又は動的な閾値の変更を、単一又は複数の一定閾値とともに使用してもよい。このような動的な閾値の変更は、閾値トリガ回路への1又はそれ以上の信号入力を動的に変えることによって、及び/又は関連するフィルタのQを変えることによって実現する。閾値トリガ回路への入力の例は、信号振幅及び閾値振幅のレベルを有している。このような単一の閾値の動的な変更を代替的な方法として使用して、例えば図5から図9に関する本書に記載の実施例に関する事象を取得してもよい。
ある実施例では、どのようなデータ簡約段階を使用しても、所定の時間周期で事象を覆うのが有益である。例としての事象マスキング処理は、マスクの時間周期の間に発生する後に続く事象が後に続く処理段階に達しないように、事象がマスクの時間カウンタをトリガする(例えば、マスクの時間周期は関連するBPFの時間周期の整数倍でよい)。
図16は、3次元面を示す。この面は、576のバンドパスフィルタの0から2000の時間分にわたる出力に対応する。各フィルタは6次のフィルタであり、576の異なる周波数をそれぞれ用いるために3つのフィルタを直列にカスケード接続する。面のトポロジーは、各カスケードにおける最終的なフィルタの出力に対応する。最終的なフィルタは曲線から成る出力を有しており、これを用いて面を発生させる。このような面を発生させるのに使用した入力信号は、異なる周波数の3つの正弦波を具えていた。多くの適用で、入力信号が様々な周波数の非常に多くの様々な正弦波を具えていることに留意されたい。それにも係わらず、入力信号が何を含んでいても、面を生成する。
入力信号を処理する方法における主要なステップは、データ抽出である。このようなデータ抽出は、フィルタ出力の中の局所ピークを特定するステップを具えており、それぞれの局所ピークは時間及び周波数方向の双方の局所ピークである。図16のトポロジーの視覚化から分かるように、非常に多くのフィルタがそれ自身の局所ピークを有しているが、図16の面の解析が、周波数及び時間方向の双方における局所ピークのピークの特定に関係している。このようなピークの1つは符号160であるが、複数の双方向局所ピークが図16に示すように面に存在する。
図17は図16のトポロジーに関するデータ抽出の結果を示す。周波数方向及び時間方向の双方における図16の面における局所ピークである各ピークが、空間の中の点で表示されている(図面パッケージの機能として、それらを面から上に延びるスパイクとして示す)。一例が符号162でマークされており、これは図16のピーク160に対応する。この図から、カーブフィッティングをどのように用いて特定したピークを表す座標を結ぶのかが分かる。このような方法では、元の信号を全体として新たな方法で処理して表示する。
抽出された実際のデータ(局所ピーク)を、信号の解析方法が採用される適用に応じて、様々な方法で使用できる。例えば、医療検査(例えば、心拍数の監視)及び機器検査(例えば、振動)といった検査環境において、図17の解析データの特定の特性を神経回路網等といった方法を用いて探査してもよい。
図18は、2の値であるQ及び441.0Hzの入力信号に応じた441.191Hzの中心周波数を具えた6次のバンドパスフィルタの入力、出力及びピーク検出の一例を示す。入力及び出力は、入力及び出力信号(符号164)の位相及び振幅がほぼ一致する際に、互いに上部に現れる。示される差は、入力信号とフィルタの中心周波数との間の0.191Hzの差による。スパイク166は、単一の時間方向におけるフィルタの出力に関するピーク検出(その値は1.00であり、裸眼で見付けるのが難しい)である。
図19は、2の値であるQ及び441.0Hzの同じ入力信号に応じた448.385Hzの中心周波数を具えた6次のバンドパスフィルタの入力、出力及びピーク検出の一例を示す。この場合、入力信号164及び出力信号165は位相が明らかに一致せず、ピーク166の値は1よりも小さい(この値は、約0.95であり、裸眼で見付けるのが難しい)。この差は、入力信号とフィルタの中心周波数との間の7.385Hzの差によるものである。時間及び周波数方向の双方における局所ピークの検出については、図19のピーク166は、図18に示すようにピーク166が隣接するフィルタの振幅よりも小さい振幅を有するため、双方の方向における局所ピークではない。
図20は、400Hz及び600Hzの正弦波に過渡的衝突成分を加えたものから成る信号を示しており、12次の低域共鳴フィルタを用いてフィルタ処理される。上側の部分168は、振幅対時間を具えた元の信号であり、符号170で示す(ハンドクラップ(hand clap)といった)衝撃が現れている。その下は、データ抽出の結果であり、それぞれのドットが検出した局所ピークの存在を示す。受信した信号を衝撃モードで解析するが、このモードは信号の中の周期成分を解析するために使用されるモードとはわずかに異なる。図20では、特定される局所ピークは、負の振幅(すなわち、ゼロよりも低いピーク値を有する)を有するピークに関する局所ピークを含んでいる。特定したピークは、周囲に対するピークであり、正の振幅を有している。周期的な解析では、正の振幅を有する隣接するピークを有する局所ピークのみを特定する。図において、衝撃170が周波数の中で低下する一連のラインとして示されている。衝撃モードでは、下部Qフィルタを正弦波モードに対して使用する。
図21は、図20の信号168の正弦波的な再構成を示しており、6次の高共鳴フィルタを用いて、400Hz及び600Hzのサイン波に過渡的衝突成分を加えたものから成る。過渡的衝突成分170はサイン波172の正弦波的な再構成に影響を及ぼさないことに留意されたい。
図22は、元の信号174のリード(lead)のうちの2つを示す12誘導ECG(心電図)信号の解析を示している。この解析は2つの信号に対して個別に実行され、解析が重なる。図20に示すように、それぞれのドットが検出される局所ピークの存在を示す。ECG信号を作り出す心臓の中で発生した電子的衝撃信号の性質を示す周波数及び時間方向に進む曲線が見受けられた。この図に見られるように、本発明の実際の適用では、上記のようなデータ抽出技術を用いて元の信号の中の非常に多くの情報を特定可能である。
一般に、上述した1又はそれ以上の装置及びその組み合わせによる入力信号の解析及び/又は処理の例は、解析−認識−再構成、圧縮−再構成、処理−再構成、シミュレーション−モデル化−予測−再構成を有しているがこれらに限定されず、これらの例を以下に示す。
解析/認識/再構成の例は:レーダー、ソナー、音声、文章、楽器、楽譜、翻訳、ヒアリングの訂正、機械故障の予測、侵入者検知、SETI(地球外知的生命体探査)、生体信号を含む。
圧縮/再構成の例は:音響機器、テレビ、衛星中継、ビデオ、電気通信を含む単一/複数の正弦波信号を含む。
処理/再構成の例は:フィルタ処理、音の高さの移動、暗号化/復号化、消音、信号除去、ノイズ除去、信号補正(例えば、多くの精密なハイファイ装置)を含む。
シミュレーション/モデル化/予測/再構成の例は:
これらの配列を組み合わせる(一方の出力部から別の入力部に供給等する)ことによって、(システムの「メモリ容量」及び「処理速度」次第で)システムのモデル化が可能となるであろう。
シミュレーションにより、より良好な予測が可能となるため:物体の誘導、天気の予測、地理的な分析、地質学的な探査、プランニング、環境災害の予測、人工知能、衛生/遺伝子工学、食料、動物、野菜、(ウイルス性/細菌性の)病気、薬学研究、(例えば、ガンの増殖の)シミュレーション、治療法に利点を有する。
図1は、第1の信号解析器を示すフローチャートである。 図2は、第2の信号解析器を示すフローチャートである。 図3は、第3の信号解析器を示すフローチャートである。 図4は、第4の信号解析器を示すフローチャートである。 図5は、本発明の一態様に係る装置の好適な実施例を示す。 図6は、本発明の一態様に係る装置の一例を示す。 図7は、振幅を変える単一の正弦波を具えた信号の一例及び信号のピークで発生する事象の例を示す。 図8は、事象ストリームの例を示す。 図9は、フィルタ応答の例を示す。 図10は、事象ストリームの例及び事象の上部を結合する曲線の例を示す。 図11は、BPFの線質係数(Q)及び/又は中心周波数フィードバック制御を与える装置の例を示す。 図12は、調波関係の判別段階によって特定される調波系の例を示す。 図13は、事象を発生させるよう使用する第1の代替的なデータ簡約法の例を示す。 図14は、事象を発生させるよう使用する第2の代替的なデータ簡約法の例を示す。 図15は、3次元曲線の例を示す。 図16は、複数のバンドパスフィルタの出力に対応する3次元面を示す。 図17は、図16の面について実行されるデータ抽出の結果を示す。 図18は、異なる2つの6次のバンドパスフィルタの入力、出力及びピーク検出を示す。 図19は、異なる2つの6次のバンドパスフィルタの入力、出力及びピーク検出を示す。 図20は、信号及び当該信号のデータ抽出の出力の概要図を示す。 図21は、図20の信号の正弦波的な再構成を示す。 図22は、信号及び当該信号のデータ抽出の出力の図20と同じような概要図を示す。

Claims (48)

  1. 信号を受信するステップと、
    複数のバンドパスフィルタで前記信号をフィルタリングするステップと、
    前記フィルタの出力からデータを抽出するステップと、を具えた信号の処理方法であって、
    前記データを抽出するステップが、前記フィルタの出力の中の局所ピークであって、それぞれの局所ピークが時間及び周波数の双方の方向の局所ピークとなる局所ピークを特定するステップを具えることを特徴とする方法。
  2. それぞれのバンドパスフィルタが少なくとも1/√2のQファクタ(factor)を有することを特徴とする請求項1に記載の方法。
  3. さらに、カーブフィティング処理を実行するステップを具えており、
    当該処理による出力が、前記抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えることを特徴とする請求項1又は2に記載の方法。
  4. さらに、前記信号をフィルタリングする際に少なくとも1のバンドパスフィルタの前記Qファクタを変更するステップを具えることを特徴とする請求項1から3のいずれか1項に記載の方法。
  5. 前記バンドパスフィルタの数が可変であり、
    前記方法が、さらに、前記信号をフィルタリングする際に前記バンドパスフィルタの数を変更するステップを具えることを特徴とする請求項1から4のいずれか1項に記載の方法。
  6. 前記複数のバンドパスフィルタが、少なくとも300のフィルタを具えることを特徴とする請求項1から5のいずれか1項に記載の方法。
  7. 前記複数のバンドパスフィルタが、前記受信した信号の周波数のオクターブ毎にほぼ1つ具えていることを特徴とする請求項1から5のいずれか1項に記載の方法。
  8. 前記複数のバンドパスフィルタが、前記受信した信号の周波数のオクターブ毎にほぼ12の整数倍具えていることを特徴とする請求項1から5のいずれか1項に記載の方法。
  9. 前記複数のバンドパスフィルタが、前記受信した信号の成分毎にほぼ3の整数倍具えていることを特徴とする請求項1から8のいずれか1項に記載の方法。
  10. さらに、3つのバンドパスフィルタのセットを用いて前記受信した信号の各成分をトラッキングするステップを具えることを特徴とする請求項9に記載の方法。
  11. さらに、調波関係の判別処理を実行するステップを具えており、
    当該処理の出力が、異なる正弦波間の調波関係を特定する調波関係のデータを具えていることを特徴とする請求項1から10のいずれか1項に記載の方法。
  12. さらに、前記信号をフィルタリングする際に、1又はそれ以上のフィルタのゲインを変更するステップを具えていることを特徴とする請求項1から11のいずれか1項に記載の方法。
  13. さらに、前記抽出したデータの説明を表示するステップを具えることを特徴とする請求項1から12のいずれか1項に記載の方法。
  14. さらに、前記抽出したデータから複雑な信号を再構成するステップを具えることを特徴とする請求項1から13のいずれか1項に記載の方法。
  15. 信号を受信するための受信部と、
    前記信号をフィルタリングするための複数のバンドパスフィルタと、
    前記フィルタの出力からデータを抽出するよう構成された処理装置と、を具えた信号を処理するためのシステムであって、
    前記データ抽出が、前記フィルタ出力の中の局所ピークであって、それぞれの局所ピークが時間及び周波数の双方の方向の局所ピークとなる局所ピークを特定するステップを具えていることを特徴とするシステム。
  16. それぞれのバンドパスフィルタが少なくとも1/√2のQファクタを有することを特徴とする請求項15に記載のシステム。
  17. 前記処理装置が、さらに、カーブフィティングの処理を実行するよう構成されており、
    当該処理による出力が、前記抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えることを特徴とする請求項15又は16に記載のシステム。
  18. 前記処理装置が、さらに、前記信号をフィルタリングする際に少なくとも1のバンドパスフィルタの前記Qファクタを変更するよう構成されていることを特徴とする請求項15から17のいずれか1項に記載のシステム。
  19. 前記処理装置が、さらに、前記信号をフィルタリングする際に前記バンドパスフィルタの数を変更するよう構成されていることを特徴とする請求項15から18のいずれか1項に記載のシステム。
  20. 前記複数のバンドパスフィルタが、少なくとも300のフィルタを具えることを特徴とする請求項15から19のいずれか1項に記載のシステム。
  21. 前記複数のバンドパスフィルタが、前記受信信号の周波数のオクターブ毎にほぼ1つ具えていることを特徴とする請求項15から19のいずれか1項に記載のシステム。
  22. 前記複数のバンドパスフィルタが、前記受信信号の周波数のオクターブ毎にほぼ12の整数倍具えていることを特徴とする請求項15から19のいずれか1項に記載のシステム。
  23. 前記複数のバンドパスフィルタが、前記受信信号の成分毎にほぼ3の整数倍具えていることを特徴とする請求項15から22のいずれか1項に記載のシステム。
  24. 前記処理装置が、さらに、3つのバンドパスフィルタのセットを用いて前記受信信号の各成分をトラッキングするよう構成されていることを特徴とする請求項23に記載のシステム。
  25. 前記処理装置が、さらに、調波関係の判別処理を実行するよう構成されており、
    前記処理の出力が、異なる正弦波間の調波関係を特定する調波関係データを具えていることを特徴とする請求項15から24のいずれか1項に記載のシステム。
  26. 前記処理装置が、さらに、前記信号をフィルタリングする際に、1又はそれ以上のフィルタのゲインを変更するよう構成されていることを特徴とする請求項15から25のいずれか1項に記載のシステム。
  27. さらに、前記抽出したデータの説明を表示するための表示装置を具えていることを特徴とする請求項15から26のいずれか1項に記載のシステム。
  28. 前記処理装置が、さらに、前記抽出したデータから複雑な信号を再構成するよう構成されていることを特徴とする請求項15から27のいずれか1項に記載のシステム。
  29. 信号を受信するステップと、
    前記信号の正弦波分離を実行するステップと、
    前記分離した正弦波からデータを抽出するステップと、
    カーブフィッティングの処理を実行するステップと、を具えた信号を処理するための方法であって、
    当該処理による出力が、前記抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えていることを特徴とする処理するための方法。
  30. 前記抽出したデータが、フィルタの正弦波出力の振幅のピークを規定することを特徴とする請求項29に記載の方法。
  31. 前記曲線データが複数の2次元曲線を規定することを特徴とする請求項29又は30に記載の方法。
  32. 前記曲線データのそれぞれの曲線が、振幅のピークに少なくとも近似する曲線を規定することを特徴とする請求項30又は31に記載の方法。
  33. 前記曲線データのそれぞれの曲線が、それぞれのフィルタの前記振幅のピークに少なくとも近似する曲線を規定することを特徴とする請求項32に記載の方法。
  34. 前記曲線データが少なくとも1の3次元表面を規定することを特徴とする請求項29又は30に記載の方法。
  35. 前記曲線データの面又は各面が、振幅のピークに少なくとも近似する面を規定することを特徴とする請求項30又は34に記載の方法。
  36. 信号を受信するための受信部と、
    前記信号の正弦波分離を実行するためのフィルタリング装置と、
    前記分離した正弦波からデータを抽出して、カーブフィッティングの処理を実行するよう構成された処理装置と、を具えた信号を処理するためのシステムであって、
    前記処理の出力が、前記抽出したデータの一部に少なくとも近似する曲線を規定する曲線データを具えていることを特徴とするシステム。
  37. 前記抽出したデータが、前記フィルタの正弦波出力の振幅のピークを規定することを特徴とする請求項36に記載のシステム。
  38. 前記曲線データが複数の2次元曲線を規定することを特徴とする請求項36又は37に記載のシステム。
  39. 前記曲線データのそれぞれの曲線が、振幅のピークに少なくとも近似する曲線を規定することを特徴とする請求項36又は37に記載のシステム。
  40. 前記曲線データのそれぞれの曲線が、それぞれのフィルタの前記振幅のピークに少なくとも近似する曲線を規定することを特徴とする請求項39に記載のシステム。
  41. 前記曲線データが少なくとも1の3次元表面を規定することを特徴とする請求項36又は37に記載のシステム。
  42. 前記曲線データの面又は各面が、振幅のピークに少なくとも近似する面を規定することを特徴とする請求項37又は41に記載のシステム。
  43. 信号を受信するステップと、
    少なくとも1のバンドパスフィルタで前記信号をフィルタリングするステップと、
    前記信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更するステップと、
    前記フィルタ又は各フィルタの出力からデータを抽出するステップと、
    を具えることを特徴とする信号を処理するための方法。
  44. 信号を受信するための受信部と、
    前記信号をフィルタリングするための少なくとも1のバンドパスフィルタと、
    前記信号をフィルタリングする際に少なくとも1のバンドパスフィルタのQファクタを変更して、前記フィルタ又は各フィルタの出力からデータを抽出するよう構成された処理装置と、
    を具えることを特徴とする信号を処理するためのシステム。
  45. 信号を受信するステップと、
    複数のバンドパスフィルタで前記信号をフィルタリングするステップと、
    前記信号をフィルタリングする際に前記バンドパスフィルタの数を変更するステップと、
    前記フィルタの出力からデータを抽出するステップと、
    を具えることを特徴とする信号を処理するための方法。
  46. 信号を受信するための受信部と、
    前記信号をフィルタリングするための複数のバンドパスフィルタと、
    前記信号をフィルタリングする際に前記バンドパスフィルタの数を変更して、前記フィルタの出力からデータを抽出するよう構成された処理装置と、
    を具えることを特徴とする信号を処理するためのシステム。
  47. 信号を受信するステップと、
    少なくとも300のバンドパスフィルタで前記信号をフィルタリングするステップと、
    前記フィルタの出力からデータを抽出するステップと、
    を具えることを特徴とする信号を処理するための方法。
  48. 信号を受信するための受信部と、
    前記信号をフィルタリングするための少なくとも300のバンドパスフィルタと、
    前記フィルタの出力からデータを抽出するよう構成された処理装置と、
    を具えることを特徴とする信号を処理するためのシステム。
JP2008551879A 2006-01-28 2007-01-26 信号解析器 Pending JP2009524812A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB0601761A GB2434649A (en) 2006-01-28 2006-01-28 Signal analyser
PCT/GB2007/000281 WO2007085853A2 (en) 2006-01-28 2007-01-26 Signal analyser

Publications (1)

Publication Number Publication Date
JP2009524812A true JP2009524812A (ja) 2009-07-02

Family

ID=36061060

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008551879A Pending JP2009524812A (ja) 2006-01-28 2007-01-26 信号解析器

Country Status (5)

Country Link
US (1) US20090030628A1 (ja)
EP (1) EP1974225B1 (ja)
JP (1) JP2009524812A (ja)
GB (1) GB2434649A (ja)
WO (1) WO2007085853A2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013053883A (ja) * 2011-09-01 2013-03-21 Kyowa Electron Instr Co Ltd リアルタイム周波数解析方法

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007077524A2 (en) * 2006-01-03 2007-07-12 Ati Technologies, Inc Harmonic edge synthesizer, enhancer and methods
GB2467724B (en) * 2008-10-28 2013-06-05 Psyanna Ltd Signal processing
JP5164807B2 (ja) * 2008-11-14 2013-03-21 アンリツ株式会社 信号解析装置
US20120004911A1 (en) * 2010-06-30 2012-01-05 Rovi Technologies Corporation Method and Apparatus for Identifying Video Program Material or Content via Nonlinear Transformations
US8527268B2 (en) 2010-06-30 2013-09-03 Rovi Technologies Corporation Method and apparatus for improving speech recognition and identifying video program material or content
US8761545B2 (en) 2010-11-19 2014-06-24 Rovi Technologies Corporation Method and apparatus for identifying video program material or content via differential signals
CN102201037B (zh) * 2011-06-14 2013-07-10 中国农业大学 农业灾害预测方法
WO2015047065A1 (en) * 2013-09-27 2015-04-02 Latvijas Universitates Agentura "Latvijas Universitates Polimeru Mehanikas Instituts" Method and device for estimating the real and imaginary parts of the frequency response
EP3058043B1 (en) * 2013-10-18 2021-07-14 ExxonMobil Chemical Patents Inc. High polymer load polyolefin adhesive compositions
GB2536226B (en) * 2015-03-09 2019-11-27 Crfs Ltd Frequency discriminator
CN106227843A (zh) * 2016-07-28 2016-12-14 云南大学 一种地质灾害野外调查信息采集系统及方法
CN110796310A (zh) * 2019-10-30 2020-02-14 黄淮学院 一种区域性地质灾害的易发性预测方法和系统
CN114694938B (zh) * 2022-03-30 2023-03-24 安徽省昌盛电子有限公司 一种电流杂波过滤系统
CN116054402B (zh) * 2023-01-09 2023-09-22 上海正泰智能科技有限公司 用于拓扑识别的电流解调方法、调制方法、设备及介质

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3731188A (en) * 1971-04-19 1973-05-01 Tracor Signal analysis of multiplicatively related frequency components in a complex signal
JPS5843699B2 (ja) * 1978-09-26 1983-09-28 日本ビクター株式会社 スペクトル表示装置
FR2534718A1 (fr) * 1982-10-19 1984-04-20 Efcis Analyseur spectral, notamment pour la reconnaissance de parole
JPS6129770A (ja) * 1984-07-20 1986-02-10 Yokogawa Medical Syst Ltd ドプラ信号分析装置
US5224170A (en) * 1991-04-15 1993-06-29 Hewlett-Packard Company Time domain compensation for transducer mismatch
JPH08248070A (ja) * 1995-03-08 1996-09-27 Anritsu Corp 周波数スペクトル分析装置
EP0977384A1 (en) * 1998-07-31 2000-02-02 Hewlett-Packard Company Method and apparatus for measuring spectral energy of mobile handset
US6448909B1 (en) * 2000-01-28 2002-09-10 The United States Of America As Represented By The Secretary Of The Navy Analog continuous wavelet transform circuit
JP2001312289A (ja) * 2000-02-21 2001-11-09 Matsushita Electric Ind Co Ltd 帯域分割用フィルタ回路ならびにそれを用いた信号分析装置および信号加工装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013053883A (ja) * 2011-09-01 2013-03-21 Kyowa Electron Instr Co Ltd リアルタイム周波数解析方法

Also Published As

Publication number Publication date
GB0601761D0 (en) 2006-03-08
WO2007085853A3 (en) 2007-10-18
WO2007085853A2 (en) 2007-08-02
US20090030628A1 (en) 2009-01-29
EP1974225A2 (en) 2008-10-01
GB2434649A (en) 2007-08-01
EP1974225B1 (en) 2019-07-03

Similar Documents

Publication Publication Date Title
JP2009524812A (ja) 信号解析器
KR101238362B1 (ko) 음원 거리에 따라 음원 신호를 여과하는 방법 및 장치
Evans et al. Method for calculating self-noise spectra and operating ranges for seismographic inertial sensors and recorders
US20020039425A1 (en) Method and apparatus for removing noise from electronic signals
KR102191736B1 (ko) 인공신경망을 이용한 음성향상방법 및 장치
Moukadem et al. Time-Frequency domain for segmentation and classification of Non-stationary signals: The Stockwell transform applied on Bio-signals and electric signals
Tenoudji Analog and digital signal analysis
Singh et al. Countermeasures to replay attacks: A review
Melgoza et al. Comparing radar receiver pulse deinterleaving performance of differing window functions for bandpass FIR filter design
Hosseini et al. Time difference of arrival estimation of sound source using cross correlation and modified maximum likelihood weighting function
JP5077847B2 (ja) 残響時間推定装置及び残響時間推定方法
WO2020250797A1 (ja) 情報処理装置、情報処理方法、及びプログラム
JP2006209123A (ja) 高調波信号の基本周波数を求める方法
JP2000296118A (ja) 生体信号解析方法および装置
Costa et al. Sparse time-frequency representations for polyphonic audio based on combined efficient fan-chirp transforms
MacDonald et al. Ultrasonic vocalizations, their recording, and bioacoustic analysis
JP6933303B2 (ja) 波源方向推定装置、波源方向推定方法、およびプログラム
Zhang et al. A matlab-based signal processing toolbox for the characterization and analysis of musical vibrato
Candy et al. Vibrational processing of a dynamic structural flight system: A multichannel spectral estimation approach
JP7461192B2 (ja) 基本周波数推定装置、アクティブノイズコントロール装置、基本周波数の推定方法及び基本周波数の推定プログラム
Chithra et al. A Comprehensive Study of Time-Frequency Analysis of Musical Signals
Kaneko Time-Frequency Analysis in Health Informatics
Varanis et al. Signal Analysis in Chaotic Systems: A Comprehensive Assessment through Time-Frequency Analysis
WO2023141565A1 (en) Data augmentation system and method for multi-microphone systems
WO2021044470A1 (ja) 波源方向推定装置、波源方向推定方法、およびプログラム記録媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100108

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20100227

A072 Dismissal of procedure [no reply to invitation to correct request for examination]

Free format text: JAPANESE INTERMEDIATE CODE: A073

Effective date: 20110531