(57)【要約】【課題】本発明は、高調波信号の基本周波数を推定するための確実な方法を提供する。【解決手段】本発明によると、高調波信号の基本周波数を推定する方法は、基本周波数仮定値（ｆ０´）を形成するステップと、基本周波数仮定値に基づいて、くし型フィルタを提供するステップと、くし型フィルタを用いて、与えられた高調波信号のフィルタリングを行うステップと、くし型フィルタにおける各歯毎に、基本周波数仮定値のテストを行うステップと、を備える。【選択図】図３
2008242431.tif 000002
【特許請求の範囲】
【請求項１】
高調波信号の基本周波数を推定する方法であって、
基本周波数仮定値（ｆ０´）を形成するステップと、
前記基本周波数仮定値に基づいて、くし型フィルタを提供するステップと、
前記くし型フィルタを用いて、与えられた高調波信号のフィルタリングを行うステップと、
前記くし型フィルタにおける各歯毎に、前記基本周波数仮定値のテストを行うステップと、
前記テストに基づいて、前記与えられた高調波信号の推定された基本周波数を表わす信号を出力するステップと、
を備える方法。
【請求項２】
前記基本周波数仮定値（ｆ０´）は、前記信号のサンプリング解像度に基づいて形成される、請求項１に記載の方法。
【請求項３】
前記くし型フィルタは、前記基本周波数仮定値（ｆ０´）と、その可能な高調波とを含む、請求項１に記載の方法。
【請求項４】
前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯で発見される第１の値と前記基本周波数仮定値から期待される第２の値との差を、所定の閾値と比較するステップを備える、請求項１に記載の方法。
【請求項５】
前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯における信号のゼロ交差点の間の間隔と前記基本周波数仮定値から期待される信号のゼロ交差点の間の間隔との対応する次数の差異を、所定の閾値と比較するステップを備える、請求項１に記載の方法。
【請求項６】
前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯における信号の自己相関のピークの位置と前記基本周波数仮定値から期待される信号の自己相関のピークの位置との間の差異を、所定の閾値と比較するステップを備える、請求項１に記載の方法。
【請求項７】
前記閾値は、前記信号に存在する干渉によって適応的に設定される、請求項４、５、または６のいずれか１項の記載の方法。
【請求項８】
高調波及び分数調波に対する前記くし型フィルタの歯のプロトタイプ割り当てパターンに基づいて、現在の前記基本周波数仮定値に重み付け値を指定するステップをさらに備える、請求項１から７のいずれか１項に記載の方法。
【請求項９】
前記正しい割り当ては、非線形方式で増幅される、請求項８に記載の方法。
【請求項１０】
前記重み付け値は、前記くし型フィルタの歯における信号のエネルギーにさらに依存する、請求項８または９に記載の方法。
【請求項１１】
前記計算された重み付け値のヒストグラムは、順時的に構成される、請求項１から１０のいずれか１項に記載の方法。
【請求項１２】
請求項１から１１のいずれか１項に記載の方法を用いて、高調波信号の基本周波数の高調波または分数調波をキャンセルする方法。
【請求項１３】
コンピュータデバイス上で作動する場合、請求項１から１２のいずれか１項に記載の方法を具現する、コンピュータソフトウェアプログラム製品。
【請求項１４】
高調波信号の基本周波数を推定するシステムであって、
基本周波数仮定値（ｆ０´）を形成する手段と、
前記基本周波数仮定値に基づいて、くし型フィルタを提供する手段と、
前記くし型フィルタを用いて、与えられた高調波信号のフィルタリングを行う手段と、
前記くし型フィルタにおける各歯毎に、前記基本周波数仮定値のテストを行う手段と、
前記テストに基づいて、推定された基本周波数を表わす信号を出力する手段と、
を備えるシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号の処理に関し、特に、高調波信号の基本周波数を発見する技術に関する。このような技法は、例えば、基底の基本周波数に基づいたモノラル録音における音源の分離、有声音/無声音の決定、または、基本周波数に基づいた性別の検出などの分野において用いられる。しかしながら、本発明は、音響の分野に限定されず、圧力センサから発せられる信号などのような他の信号にも適用することができる。
【背景技術】
【０００２】
音声信号は、多くの高調波成分を含んでいる。これらの高調波成分の基本周波数に関する知識は、数多くの方法で展開することができる。一つの非常に重要な例には、音源の分離がある。音響記録を行う場合、時々、同時に多数の音源が存在したりすることがある。これらは、異なる音声信号、雑音（例えば、ファンの雑音）、あるいは、類似の信号であり得る。信号をさらに分析するためには、先ず、これらの干渉信号を分離する必要がある。通常の適用例としては、音声認識または音響場面分析などがある。
【０００３】
高調波信号の基本周波数を判定するために、異なる従来技術の手法が公知となっている。最も一般的なものとして、自己相関関数を用いたものがある（非特許文献１を参照）。ここで、信号は、一連の帯域通過フィルタを介して周波数帯域に分離され、それぞれの周波数帯域毎に自己相関が判定されて、高調波関係にある周波数は、遅延領域においてピークを共有する。それによって、実際の遅延の倍数及び部分に相当する遅延においてピークがまた発生する。このような付加的ピークは、基本周波数の判定において主ピークと干渉する。
【０００４】
その内容の全体が本明細書に参照として取り入れられている、同一の発明者による欧州特許出願第０５００４０６６号（特許文献１）では、同様に、遅延／間隔領域においてピークを共有する個々の周波数チャネル内の一部の次数のゼロ交差点（ｚｅｒｏ  ｃｒｏｓｓｉｎｇ）の間の間隔を計算することによって、自己相関の使用に取って代わる方法を提案している。つまり、ゼロ交差の間隔の計算を介してチャネルの基本周波数が推定される。同一の基本周波数から高調波が発生すると、高調波は、ゼロ交差間隔を共有する。
【０００５】
例えば、基本周波数に属するチャネル内の２つのゼロ交差点の間の間隔は、第１の高調波における３つのゼロ交差点の間の間隔、及び第２の高調波における４つのゼロ交差点の間の間隔として再び発見される（より詳細には、以下の特許文献１および非特許文献２の著述を参照されたい）。
【０００６】
このような３つまたは４つのゼロ交差点の間の間隔は、さらに、それぞれ２次及び３次などのより高次元のゼロ交差間隔として呼ばれる。しかしながら、このような場合にも、擬似サイドピークが現われる。
【０００７】
非特許文献３では、異なる経路が伴われる。ここでは、「ハーモニック・シーブ（ｈａｒｍｏｎｉｃ  ｓｉｅｖｅ）」とも呼ばれる、くし型フィルタ（ｃｏｍｂ  ｆｉｌｔｅｒ）が、その高調波及び基本周波数における歯（ｔｅｅｔｈ）で設定される。それぞれの歯において発見されたエネルギーは、異なる基本周波数仮定値に対して合算される。仮定値と真の基本周波数とが一致する場合、くし型フィルタ内のすべての歯は、最大をもたらす高いエネルギーを有する。前述した方法においても、真の基本周波数の高調波と分数調波とにおいてサイドピークが生じる。
【特許文献１】欧州特許出願第０５００４０６６号
【非特許文献１】Ｇ．Ｈｕ及びＤ．Ｗａｎｇ：ピッチ・トラッキング及び振幅に基づいたモノラル音声分離（Ｍｏｎａｕｒａｌ  ｓｐｅｅｃｈ  ｓｅｇｒｅｇａｔｉｏｎ  ｂａｓｅｄ  ｏｎ  ｐｉｔｃｈ  ｔｒａｃｋｉｎｇ  ａｎｄ  ａｍｐｌｉｔｕｄｅ）．ＩＥＥＥ  Ｔｒａｎｓ．Ｏｎ  Ｎｅｕｒａｌ  Ｎｅｔｗｏｒｋｓ（２００４）
【非特許文献２】Ｍａｒｔｉｎ  Ｈｅｃｋｍａｎｎ及びＦｒａｎｋ  Ｊｏｕｂｌｉｎによる論文「ピッチに基づいたロボットに対する音源分離（Ｓｏｕｎｄ  Ｓｏｕｒｃｅ  Ｓｅｐａｒａｔｉｏｎ  ｆｏｒ  ａ  Ｒｏｂｏｔ  Ｂａｓｅｄ  ｏｎ  Ｐｉｔｃｈ）」，Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｉｎｔｅｌｌｉｇｅｎｔ  Ｒｏｂｏｔｓ  ａｎｄ  Ｓｙｓｔｅｍｓ（ＩＲＯＳ）、Ｅｄｍｏｎｔｏｎ、Ｃａｎａｄａ、Ａｕｇｕｓｔ  ２００５、ｐｐ．２０３−２０８）
【非特許文献３】Ｈ．Ｄｕｉｆｈｕｉｓらによる論文「Ｈ．Ｄｕｉｆｈｕｉｓ、Ｌ．Ｗｉｌｌｅｍｓ及びＲ．Ｓｌｕｙｔｅｒ：音声におけるピッチの測定：ピッチ認識のゴールドスタイン理論の具現（Ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  ｐｉｔｃｈ  ｉｎ  ｓｐｅｅｃｈ：Ａｎ  ｉｍｐｌｅｍｅｎｔａｔｉｏｎ  ｏｆ  Ｇｏｌｄｓｔｅｉｎ’ｓ  ｔｈｅｏｒｙ  ｏｆ  ｐｉｔｃｈ  ｐｅｒｃｅｐｔｉｏｎ．）」，Ｊ．Ａｃｏｕｓｔ．Ｓｏｃ．Ａｍ．Ｐｐ．１５６８−１５８０（１９８２）
【発明の開示】
【発明が解決しようとする課題】
【０００８】
したがって、本発明の目的は、高調波信号の基本周波数を推定するための確実な方法を提供することである。
【課題を解決するための手段】
【０００９】
このような目的は、独立請求項の特徴により、本発明によって達成される。好ましい実施形態は、従属請求項において定義されている。
【００１０】
本発明の第１態様によると、高調波信号の基本周波数を推定する方法は、基本周波数仮定値（ｆ０´）を形成するステップと、前記基本周波数仮定値に基づいて、くし型フィルタを提供するステップと、前記くし型フィルタを用いて、与えられた高調波信号のフィルタリングを行うステップと、前記くし型フィルタにおける各歯毎に、前記基本周波数仮定値のテストを行うステップと、を備える。該方法はさらに、前記テストに基づいて、前記与えられた高調波信号の推定された基本周波数を表わす信号を出力するステップを備える。
【００１１】
前記基本周波数仮定値（ｆ０´）は、前記信号のサンプリング解像度に基づいて形成することができる。前記くし型フィルタは、前記基本周波数仮定値（ｆ０´）と、その可能な高調波とを含むことができる。
【００１２】
さらに、前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯で発見される第１の値と前記基本周波数仮定値から期待される第２の値との差異を、所定の閾値と比較するステップを備えることができる。
【００１３】
別の態様によると、前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯における信号のゼロ交差点の間の間隔と前記基本周波数仮定値から期待される信号のゼロ交差点の間の間隔との差を、所定の閾値と比較するステップを備えることができる。あるいは、前記基本周波数仮定値のテストを行うステップは、前記くし型フィルタの歯における信号の自己相関のピークの位置と前記基本周波数仮定値から期待される信号の自己相関のピークの位置との間の差異を、所定の閾値と比較するステップを備えることができる。すべての場合において、前記閾値は、前記信号に存在する干渉（ｄｉｓｔｕｒｂａｎｃｅ）によって、適応的に設定されることができる。前記方法は、さらに、高調波及び分数調波に対する前記くし型フィルタの歯のプロトタイプ割り当てパターンに基づいて、現在の基本周波数仮定値に重み付け値を指定するステップを備えることができる。なお、正しい割り当ては、非線形方式で増幅することができる。前記重み付け値は、前記くし型フィルタの歯における信号のエネルギーにさらに依存することができる。
【００１４】
本発明のさらに別の態様によると、前記計算された重み付け値のヒストグラムは、順時的に構成することができる。
【００１５】
上記方法は、高調波信号において、基本周波数の高調波または分数調波をキャンセルするために用いることができる。
【００１６】
本発明は、高調波信号の基本周波数の抽出において、結果を向上するために活用することができる。とりわけ、真の基本周波数の高調波及び分数調波における擬似サイドピークの問題点が、提案された方法によって著しく軽減される。
【００１７】
上述したような、及び更なる本発明の態様及び長所は、添付の図面を参照した、以下の本発明の詳細な説明に基づいて、より一層明らかになるだろう。
【発明を実施するための最良の形態】
【００１８】
図１は、本発明の第１実施形態による高調波信号の基本周波数を推定する方法１００のフローチャートである。
【００１９】
ステップ１１０において、与えられた高調波信号の基本周波数に対する仮定値が形成される。ステップ１２０において、ステップ１１０で形成された基本周波数仮定値に基づいて、くし型フィルタが提供あるいは設定される。当業者によく知られているように、くし型フィルタの伝達関数は、髪の毛のくしと類似している。スペクトル領域に多数の「歯」を有し、ここに情報が保持される。このような歯の外部の情報は、除去される。
【００２０】
ここで、くし型フィルタは、調査された基本周波数と、それの可能な高調波とを含むように設定される。すなわち、くし型フィルタは、くしの「歯」が、調査された基本周波数及びそれの可能な高調波において発生するように設定される。
【００２１】
ステップ１３０において、高調波信号が、くし型フィルタを用いてフィルタリングを受ける。次いで、ステップ１４０において、くし型フィルタのそれぞれの歯に対し、基本周波数仮定値がテストされる。このテストの間に、基本周波数仮定値から期待される値は、くし型フィルタの歯で発見される値と比較され、発見された偏差に基づいて、当該歯が仮定値に属するか否かが検討される。それについて用いられる閾値は、期待値に対して相対的に、あるいは絶対的に設定することができる。
【００２２】
現在調査している基本周波数が、信号の真の基本周波数と一致する場合、くし型フィルタのすべての歯に、高調波によって信号が生じる。一部の歯が空いていると、その基底のチャネルに、現在調査している基本周波数の高調波ではない周波数によって信号が生じたことを意味し、これは、現在調査している基本周波数は、信号の真の基本周波数ではなく、高調波または分数調波であることを暗示する。
【００２３】
真の基本周波数を推定するために、すべての可能な基本周波数は、上述の方法でテストされる。
【００２４】
図２は、より確実に高調波信号の基本周波数の時間的経過を発見する方法のフローチャートを示した図であり、ここで、本発明のさらに別の実施形態に関する高調波信号の基本周波数の推定方法が採用される。特に、特許文献１の上記ゼロ交差に基づいたアルゴリズムと提案された方法との組み合わせについて説明する。しかし、提案された方法は、例えば、非特許文献１に提案されているように、基本周波数の判定のための他の技法と組み合わせることもできる。
【００２５】
その準備として、ステップ２１０において、アナログからデジタルへ信号が変換され、ステップ２２０において、一連の帯域通過フィルタまたはフィルタ・バンクを介して、周波数領域に変換される。フィルタ・バンクを用いた周波数領域での変換の結果として、信号は、フィルタ帯域幅に与えられた解像度でその周波数成分に分離する一方、帯域通過信号である、このような周波数成分のそれぞれに対する時間情報は保持する。次いで、それぞれの帯域通過信号毎に、現在の基本周波数仮定値に対する関係にかかる情報が収集され得る。
【００２６】
以下において、ゼロ交差間隔が用いられる場合、現在の基本周波数仮定値に対する相違する帯域通過信号の関係の算定が行われる方法について説明する。
【００２７】
真の基本周波数を探し出すために、すべての可能な基本周波数がスキャンされ、基本周波数仮定値として使用されなければならない。ゼロ交差点の間の間隔が基本周波数の推定における基礎となる場合、基本周波数に対する妥当な離散化（ｄｉｓｃｒｅｔｉｚａｔｉｏｎ）は、サンプリング解像度である。サンプリング率を１６ｋＨｚ、最小基本周波数を１００Ｈｚと仮定する。これは、１６０個のサンプルに相当するゼロ交差点の間の間隔に該当し、第１の基本周波数仮定値として用いることができる。第２の基本周波数仮定値として用いることのできる、次の可能な基本周波数は、１５９個のサンプルに相当する間隔を有し、したがって、１００．３Ｈｚの周波数を有する。可能な基本周波数の範囲は、自由に決定することができ、信号のサンプリング・レイトのみによって制限される。
【００２８】
それぞれの帯域通過信号毎に、ステップ２３０において、ゼロ交差点の判定を行うことができる。また、連続したゼロ交差点の間の間隔が計算されてもよい。これは、調査中にある帯域通過信号において、支配的な周波数または基本周波数の非常に正確な推定を可能とする。さらに、３つのゼロ交差点の間の間隔が計算され、２次ゼロ交差間隔として称することができる。このように、与えられた次数までゼロ交差間隔が計算され得る。実際的な最大次数の値は、７である。
【００２９】
ステップ２４０において、間隔ヒストグラムが構成される。先ず、ステップ２４１において、スキャンされた基本周波数仮定値毎に、該当のくし型フィルタが設定される。くし型フィルタは、帯域通過信号に基づき、周波数領域において設計される。通過帯域がくし型フィルタの歯に該当する周波数の一つを含む帯域通過信号がフィルタを通過し、他の信号は拒否される。くし型フィルタを設定する場合、どの次数までゼロ交差間隔を計算すべきであるのかを考慮しなければならない。また、この次数まで歯が設定される。現在の基本周波数ｆ０´は１００Ｈｚとし、最大ゼロ交差間隔の次数を５とすると、くしは、１００、２００、３００、４００、及び５００Ｈｚの周波数に相当するチャネルを構成するようになる（図３ａと比較）。
【００３０】
ここで、図３ａは、基本周波数仮定値が１００Ｈｚである場合における５つの歯を有するくし型フィルタを示した図である。
【００３１】
図３ｂは、信号の真の基本周波数と基本周波数仮定値とが一致する（両方とも１００Ｈｚ）場合におけるくし型フィルタの割り当てを示した図である。
【００３２】
図３ｃは、基本周波数仮定値が真の基本周波数の２倍である場合（ｆ０´=２００Ｈｚ及びｆ０=１００Ｈｚ）におけるくし型フィルタの割り当てを示した図である。
【００３３】
図３ｄは、基本周波数仮定値が真の基本周波数の半分である場合（ｆ０´=５０Ｈｚ及びｆ０=１００Ｈｚ）におけるくし型フィルタの割り当てを示した図である。この場合、基本周波数仮定値の第１の分数調波（１／２）の倍数の歯がくしに含まれる。
【００３４】
図３ｅは、基本周波数仮定値と信号の真の基本周波数とが一致する場合（両方とも１００Ｈｚ）における基本周波数仮定値（１．ｄ参照）の第１の分数調波（１／２）の倍数の歯で拡張されるくし型フィルタの割り当てを示した図である。
【００３５】
ステップ２４２において、くし型フィルタにおけるチャネルのゼロ交差間隔が、現在の基本周波数のゼロ交差間隔と比較される。それによって、くしの歯に対して仮定されたチャネルの次数が考慮され得る（例えば、１００Ｈｚチャネルは、第１の次数と比較され、２００Ｈｚチャネルは、第２の次数と比較される。以下、同様）。チャネルを現在の基本周波数と比較する代わりに、平均または中間値として、平均値を用いることができる。
【００３６】
本発明の一実施形態において、くし型フィルタの歯は、くし型フィルタチャネルにおいて測定される実際の周波数値と、現在調査中の基本周波数とに基づいて、現在の基本周波数の高調波である周波数によって信号を生じるものまたは信号を生じないものとしてラベリングされてもよい。つまり、比較値（すなわち、現在の基本周波数）から、それぞれの歯の偏差によって、現在の基本周波数に属するものとして、または属さないものとしてラベリングされてもよい。このような比較において、許容偏差の閾値が導入され得る。
【００３７】
現在の基本周波数ｆ０｀が、信号の真の基本周波数ｆ０と一致する場合、くしのすべての歯がラベリングされたり、あるいは設定されてもよい（図３ｂと比較）。現在の基本周波数ｆ０｀が真の基本周波数（第１の高調波）の２倍である場合、くしにおける１つおきの歯のみが、ラベリングされ、あるいは設定される（図３ｃと比較）。最後に、現在の基本周波数が真の基本周波数の半分である場合（第１の分数調波）、くしのすべての歯がラベリングされるか、設定され、これに加え、現在の基本周波数の半分の倍数の歯がラベリングされたり、設定されてもよい（図３ｄと比較）。後者を検出するために、現在の基本周波数の半分の倍数の周波数が、くし型フィルタに含まれ得る。現在の基本周波数が真の基本周波数と同一である場合における、第１の分数調波の倍数によって拡張されたくし型フィルタの割り当てが、図３ｅに示されている。
【００３８】
次のステップ２４３において、現在の基本周波数が真の基本周波数の高調波または分数調波である場合に発見される典型的な割り当てパターンと比較することにより、発見されたくし型フィルタの割り当てパターンに対する重み付け値が判定される。
【００３９】
図３に示されているようなくし型フィルタに対する、これらの予め定義されたプロトタイプ割り当てパターンに基づいて誤ったパターンを区分し、正しいパターンを強調する規則を定めることができる。一つの戦略として、非線形的な方式で、正しい割り当てパターンを増幅し、誤った割り当てパターンを抑制することもあり得る。正しい割り当てが、最大の重み付け値を獲得し、選択された高調波及び分数調波の割り当てが、重み付け値０の結果となるように、歯の割り当てを組み合わせる別の手法があり得る。
【００４０】
すなわち、割り当てパターンに基づいて、真の基本周波数のこのような高調波及び分数調波を抑制する方法を開発することができる。つまり、テストされる基本周波数が真の基本周波数である場合、くしの歯の割り当てパターンの情報を用い、テストされる基本周波数が高調波または分数調波である場合、典型的な割り当てパターンを用いて、テストされる基本周波数のヒストグラムにおいて高調波及び分数調波のピークを抑制する方法が適用されてもよい。
【００４１】
ステップ２４４において、２次元のヒストグラムが形成される。ヒストグラムは、ｘ軸に時間を示し、ｙ軸に、異なる基本周波数仮定値のゼロ交差間隔を示す。ヒストグラムに表示された値は、それらの累積発生を示すものである。このような累積発生を計算するためには、ステップ２４３において判定された重み付け値が、ヒストグラムに加えられる。
【００４２】
次に、ステップ２５０において、基本周波数ｆ０をトラッキングする方法が続けられる。
【００４３】
図４（ａ及びｂ）は、特許文献１または非特許文献２に開示されたように、計算されたゼロ交差間隔のヒストグラムに基づいて基本周波数を判定した結果（ａ）を、本発明に関連して提案された方法を追加的に用いた場合の結果（ｂ）と比較する。
【００４４】
第１の高調波と、第１及び第２の分数調波がキャンセルされるように、割り当てが組み合わせられる。ｘ軸には、時間が秒単位で与えられ、ｙ軸には、ゼロ交差点の間の間隔がミリ秒単位で与えられている。すなわち、ヒストグラムは２次元であって、ｘ軸に時間を、ｙ軸に異なる基本周波数仮定値のゼロ交差間隔を示す。ヒストグラムに表示される値は、それらの累積発生である。基本周波数に関する情報を抽出するために用いられる方法によって、ｙ軸は、さらに、自己相関のピークの遅延または基本周波数の一部の類似した周波数表示を示すことができる。示された間隔の値は、周波数に直接変換されることが可能である。
【００４５】
図４ｂより、ヒストグラムにおける高調波及び分数調波の大幅な減少が明らかに見られる。
【００４６】
基本周波数の抽出のために、くし型フィルタを活用する最新の手法において、信号を周波数帯域に分離するために採用される、先行する帯域通過フィルタの周波数選択度によって、くし型フィルタの精度が決定される（例えば、非特許文献３）。選択度とフィルタの立ち上がり時間との間には、トレードオフの関係があるはずである。他の効果を無視すると、立ち上がり時間の増大は、取得可能な選択度を制限する。支配的な周波数の推定のために帯域通過信号のゼロ交差間隔を追加に用いた場合、立ち上がり時間を増大しなくても、選択度の向上が可能となる。帯域通過フィルタによって与えられるものよりも高い精度で歯を基本周波数とラベリングするステップは、提案された方法を、このようなラベリングが行われないために、その後の抑制が不可能である従来の技術の方法と明らかに区分するものである。
【００４７】
実際の適用例において、本発明は、処理対象の音響信号を表現する信号が供給され、推定された基本周波数を表わす信号を出力するコンピュータシステムとして具現することができる。このような出力信号は、例えば、音声認識及び補聴器等で有用な、例えば、音源の分離等のような他の適用例のために用いることができる。
【図面の簡単な説明】
【００４８】
【図１】本発明の第一実施形態による高調波信号の基本周波数を推定する方法のフローチャートである。
【図２】本発明の別の実施形態による高調波信号の基本周波数を推定する方法のフローチャートである。
【図３】くし型フィルタと周波数との関係を示した図である。
【図４】ゼロ交差間隔のヒストグラムが計算された場合の基本周波数の推定の結果を比較した図である。
【図１】
2008242431.tif 000003
【図２】
2008242431.tif 000004
【図３】
2008242431.tif 000005
【図４】
2008242431.tif 000006
【外国語明細書】
１．TITLE OF INVENTION
2008242431.tif 000007
２．DETAILED DESCRIPTION OF INVENTION
2008242431.tif 000008
2008242431.tif 000009
2008242431.tif 000010
2008242431.tif 000011
2008242431.tif 000012
2008242431.tif 000013
2008242431.tif 000014
2008242431.tif 000015
2008242431.tif 000016
2008242431.tif 000017
2008242431.tif 000018
2008242431.tif 000019
2008242431.tif 000020
2008242431.tif 000021
３．BRIEF DESCRIPTION OF DRAWINGS
2008242431.tif 000022
CLAIMS
2008242431.tif 000023
2008242431.tif 000024
2008242431.tif 000025
ABSTRACT
2008242431.tif 000026
REPRESENTATIVE OF DRAWING
Fig. 3
Figure 1
2008242431.tif 000027
Figure 2
2008242431.tif 000028
Figure 3
2008242431.tif 000029
Figure 4
2008242431.tif 000030