(57)【要約】ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができるピッチ周波数推定装置を提供する。この装置において、スペクトル抽出部（１０４）は、音声スペクトルからピッチ調波スペクトルを抽出する。スペクトル平均値計算部（１０６）は、スペクトル抽出部（１０４）によって抽出されたピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する。推定部は、スペクトル平均値計算部（１０６）によって計算された平均値を用いてピッチ周波数を推定する。【特許請求の範囲】
【請求項１】
音声スペクトルからピッチ調波スペクトルを抽出する抽出手段と、
前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算手段と、
前記平均値を用いてピッチ周波数を推定する推定手段と、
を有するピッチ周波数推定装置。
【請求項２】
前記ピッチ調波スペクトルのパワの加算値を、前記複数のピッチ周波数候補の各々に対応づけて計算する加算値計算手段をさらに有し、
前記推定手段は、
前記加算値を用いてピッチ周波数を推定する、
請求項１記載のピッチ周波数推定装置。
【請求項３】
前記推定手段は、
前記平均値および前記加算値を、前記複数のピッチ周波数候補の各々に対応づけて互いに乗算する乗算手段と、
前記複数のピッチ周波数候補のうち、前記乗算手段による乗算の結果の最大値に対応するピッチ周波数候補を、推定ピッチ周波数として決定する決定手段と、
を有する請求項２記載のピッチ周波数推定装置。
【請求項４】
前記平均値計算手段は、
前記音声スペクトルにおけるパワの最大値に対応する周波数成分を基準周波数として用いて、前記平均値の計算を行う、
請求項２記載のピッチ周波数推定装置。
【請求項５】
前記加算値計算手段は、
前記音声スペクトルにおけるパワの最大値に対応する周波数成分を基準周波数として用いて、前記加算値の計算を行う、
請求項２記載のピッチ周波数推定装置。
【請求項６】
前記加算値のべき乗を計算するべき乗計算手段をさらに有し、
前記乗算手段は、
前記べき乗計算手段による計算の結果を前記平均値に乗算し、
前記べき乗計算手段は、
前記べき乗の計算に用いられる乗数を可変に設定する、
請求項３記載のピッチ周波数推定装置。
【請求項７】
前記平均値計算手段は、
前記加算値を用いて、前記平均値の計算を行う、
請求項２記載のピッチ周波数推定装置。
【請求項８】
前記ピッチ調波スペクトルの振幅を制限する振幅制限手段をさらに有する、
請求項２記載のピッチ周波数推定装置。
【請求項９】
前記音声スペクトルの有声性を判定する判定手段をさらに有し、
前記抽出手段は、
前記判定手段による判定の結果、前記音声スペクトルの有声性が所定レベル以下の場合は、前記ピッチ調波スペクトルの抽出を回避する、
請求項２記載のピッチ周波数推定装置。
【請求項１０】
音声スペクトルからピッチ調波スペクトルを抽出する抽出ステップと、
前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算ステップと、
前記平均値を用いてピッチ周波数を推定する推定ステップと、
を有するピッチ周波数推定方法。
【請求項１１】
音声スペクトルからピッチ調波スペクトルを抽出する抽出ステップと、
前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算ステップと、
前記平均値を用いてピッチ周波数を推定する推定ステップと、
をコンピュータに実現させるためのピッチ周波数推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、ピッチ周波数推定装置およびピッチ周波数推定方法に関し、特に、周波数領域でピッチ周波数推定を行うピッチ周波数推定装置およびピッチ周波数推定方法に関する。
【背景技術】
【０００２】  一般に、時間領域または周波数領域において音声のピッチ周波数を推定する方法としては、音声波形の自己相関関数による自己相関法や、ＬＰＣ（Ｌｉｎｅａｒ  Ｐｒｅｄｉｃｔｉｖｅ  Ｃｏｄｉｎｇ）分析の残差信号の自己相関関数による変形相関法などが知られている。【０００３】  また、雑音抑圧や音声符号化などの音声処理を周波数領域において行う場合は、周波数領域でピッチ周波数を推定すると整合性が良くなることがある。周波数領域でのピッチ周波数推定方法としては、周波数スペクトルに対する自己相関関数の最大化によりピッチ周波数を算出する方法があり、その一般式は次の式（１）によって表される。この式において、自己相関関数Ｒ（ｉ）を最大にするピッチ周波数候補ｉを推定ピッチ周波数とする。
【数１】
2006006366.tif 000002
ここで、ｋは離散周波数成分であり、Ｐ（ｋ）はピッチ調波スペクトルのパワであり、ＰＭＩＮおよびＰＭＡＸはそれぞれピッチ周波数候補ｉの最小値および最大値である。【０００４】  ところで、周波数領域での自己相関を用いたピッチ周波数推定方法では、音声信号のホルマントの影響により誤って倍のピッチ周波数が算出されてしまうこと（倍ピッチ周波数誤り）がある。【０００５】  ホルマントの影響を低減しつつピッチ周波数推定を行う従来の方法としては、例えば、非特許文献１に開示されたものがある。この方法では、スペクトル包絡の情報でスペクトルを平坦化した後のスペクトルが用いられる。
【非特許文献１】″Ａ  ｓｐｅｃｔｒａｌ  ａｕｔｏｃｏｒｒｅｌａｔｉｏｎ  ｍｅｔｈｏｄ  ｆｏｒ  ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  ｔｈｅ  ｆｕｎｄａｍｅｎｔａｌ  ｆｒｅｑｕｅｎｃｙ  ｏｆ  ｎｏｉｓｅ−ｃｏｒｒｕｐｔｅｄ  ｓｐｅｅｃｈ″，Ｍ．Ｌａｈａｔ，ＩＥＥＥ  Ｔｒａｎｓ．ｏｎ  Ａｃｏｕｓｔｉｃｓ，Ｓｐｅｅｃｈ，ａｎｄ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，ｖｏｌ．ＡＳＳＰ−３５，ｎｏ．６，ｐｐ．７４１−７５０，１９８７
【発明の開示】
【発明が解決しようとする課題】
【０００６】  しかしながら、上記従来のピッチ周波数推定方法では、スペクトルの平坦化処理を伴うため、ピッチ周波数推定に要する演算量が増大するという問題がある。【０００７】  本発明の目的は、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができるピッチ周波数推定装置およびピッチ周波数推定方法を提供することである。
【課題を解決するための手段】
【０００８】  本発明のピッチ周波数推定装置は、音声スペクトルからピッチ調波スペクトルを抽出する抽出手段と、前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算手段と、前記平均値を用いてピッチ周波数を推定する推定手段と、を有する構成を採る。【０００９】  本発明のピッチ周波数推定方法は、音声スペクトルからピッチ調波スペクトルを抽出する抽出ステップと、前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算ステップと、前記平均値を用いてピッチ周波数を推定する推定ステップと、を有するようにした。【００１０】  本発明のピッチ周波数推定プログラムは、音声スペクトルからピッチ調波スペクトルを抽出する抽出ステップと、前記ピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する平均値計算ステップと、前記平均値を用いてピッチ周波数を推定する推定ステップと、をコンピュータに実現させるようにした。
【発明の効果】
【００１１】  本発明によれば、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施の形態に係るピッチ周波数推定装置の構成を示すブロック図
【図２Ａ】本発明の一実施の形態において、抽出された音声パワスペクトルの例を示す図
【図２Ｂ】本発明の一実施の形態において、乗数をある値に設定した条件の下で平均値および加算値を乗算した結果を示す図
【図２Ｃ】本発明の一実施の形態において、乗数を他の値に設定した条件の下で平均値および加算値を乗算した結果を示す図
【発明を実施するための最良の形態】
【００１３】  以下、本発明の実施の形態について、図面を用いて詳細に説明する。【００１４】  図１は、本発明の一実施の形態に係るピッチ周波数推定装置の構成を示すブロック図である。ピッチ周波数推定装置１００は、ハニング窓部１０１、ＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）部１０２、有声性判定部１０３、スペクトル抽出部１０４、スペクトル振幅制限部１０５、スペクトル平均値計算部１０６、スペクトル加算部１０７、べき乗計算部１０８、乗算部１０９および最大値抽出部１１０を有する。【００１５】  ハニング窓部１０１は、所定時間単位のフレーム単位に分割された入力音声信号に対して、ハニングウィンドウなどを利用した窓掛け処理を施してＦＦＴ部１０２に出力する。【００１６】  ＦＦＴ部１０２は、ハニング窓部１０１から入力されたフレーム、つまりフレーム単位に分割された音声信号に対してＦＦＴを行って音声信号を周波数領域に変換する。これにより、音声パワスペクトルを取得する。よって、フレーム単位の音声信号は、所定の周波数帯域を有する音声パワスペクトルとなる。このようにして生成された音声パワスペクトルは、有声性判定部１０３、スペクトル抽出部１０４およびスペクトル振幅制限部１０５に出力される。【００１７】  有声性判定部１０３は、ＦＦＴ部１０２から音声パワスペクトルの有声性、つまり元の音声信号が有声であるか無声であるかを判定する。判定結果は、スペクトル抽出部１０４に出力される。【００１８】  スペクトル抽出部１０４は、有声性判定部１０３によって音声パワスペクトルが有声性なしと判定された場合、ピッチ調波スペクトルの抽出を回避する。これによって、スペクトル抽出部１０４の演算量、ひいてはピッチ周波数推定装置１００の全体の演算量を低減することができる。【００１９】  一方、音声パワスペクトルが有声性ありと判定された場合、スペクトル抽出部１０４はピッチ調波スペクトルの抽出を行う。より具体的には、音声パワスペクトルにおけるピークを抽出することにより、ピッチ調波スペクトルの抽出を行う。【００２０】  また、スペクトル抽出部１０４は、スペクトル振幅制限部１０５による音声パワスペクトルの振幅制限が行われた場合、抽出されたピッチ調波スペクトルにその振幅制限の結果を反映させることにより、ピッチ調波スペクトルの振幅を制限する。このようにして、ピッチ周波数推定の精度に与えられ得るホルマントの影響を低減することができる。ピッチ調波スペクトルは、スペクトル平均値計算部１０６およびスペクトル加算部１０７に出力される。【００２１】  スペクトル振幅制限部１０５は、ＦＦＴ部１０２によって取得された音声パワスペクトルの振幅が所定の閾値を超過しないように制限する。音声パワスペクトルの振幅制限の結果は、スペクトル抽出部１０４に出力される。【００２２】  スペクトル平均値計算部１０６は、スペクトル抽出部１０４からのピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する。すなわち、ピッチ調波スペクトルにおいて、ピッチ周波数候補を所定の最小値から所定の最大値までシフトさせながら、ピッチ周波数候補の整数倍に相当する周波数成分のパワの平均値を計算する。計算された平均値は、乗算部１０９に出力される。【００２３】  また、スペクトル平均値計算部１０６は、平均値の計算を行うとき、パワの最大値に対応する周波数成分を、平均値計算対象の周波数帯域における基準周波数として用いる。【００２４】  具体的には、基準周波数からピッチ周波数候補の整数倍に相当する周波数を減算して得られる周波数におけるパワと、基準周波数からピッチ周波数候補の整数倍に相当する周波数を加算して得られる周波数におけるパワと、を用いて、平均値の計算を行う。これにより、音声の準周期特性および雑音の影響ならびにピッチ周波数推定誤差により生じるピッチ高調波における誤差の累積を低減することができ、より正確にピッチ周波数の推定を行うことができる。【００２５】  なお、ピッチ調波スペクトルのパワの平均値は、後述するピッチ調波スペクトルのパワの加算値を特定の値で除して得られる値である。よって、スペクトル平均値計算部１０６は、スペクトル加算部１０７によって計算された加算値を取得し、これを用いて平均値の算出を行っても良い。【００２６】  スペクトル加算部１０７は、スペクトル抽出部１０４からのピッチ調波スペクトルのパワの加算値を、複数のピッチ周波数候補の各々に対応づけて計算する。すなわち、ピッチ調波スペクトルにおいて、ピッチ周波数候補を所定の最小値から所定の最大値までシフトさせながら、ピッチ周波数候補の整数倍に相当する周波数成分のパワを加算する。パワの加算によって得られた加算値はべき乗計算部１０８に出力される。【００２７】  また、スペクトル加算部１０７は、パワの加算を行うとき、パワの最大値に対応する周波数成分を、加算値計算対象の周波数帯域における基準周波数として用いる。【００２８】  具体的には、基準周波数からピッチ周波数候補の整数倍に相当する周波数を減算して得られる周波数におけるパワと、基準周波数からピッチ周波数候補の整数倍に相当する周波数を加算して得られる周波数におけるパワと、を用いて、加算値の計算を行う。これにより、音声の準周期特性および雑音の影響ならびにピッチ周波数推定誤差により生じるピッチ高調波における誤差の累積を低減することができ、より正確にピッチ周波数の推定を行うことができる。【００２９】  べき乗計算部１０８は、スペクトル加算部１０７によって算出された加算値のべき乗の値を計算する。算出されたべき乗の値は乗算部１０９に出力される。また、べき乗計算部１０８は、べき乗の計算に用いられる乗数を可変に設定する。乗数の可変設定つまり乗数の調整については後述する。【００３０】  乗算部１０９および最大値抽出部１１０の組み合わせは、複数のピッチ周波数候補の各々に対応づけて計算された平均値を用いてピッチ周波数を推定する推定部を構成する。【００３１】  推定部において、乗算部１０９は、ピッチ調波スペクトルのパワの平均値とピッチ調波スペクトルのパワの加算値とを、複数のピッチ周波数候補の各々に対応づけて乗算する。より具体的には、加算値のべき乗計算結果を平均値に乗算する。乗算結果は、最大値抽出部１１０に出力される。【００３２】  最大値抽出部１１０は、乗算部１０９で計算された乗算結果の最大値を抽出する。また、所定の最小値から所定の���大値までの複数のピッチ周波数候補のうち、乗算結果が最大となるときのピッチ周波数候補を推定ピッチ周波数として決定し、図示されない後段の処理部に出力する。【００３３】  次いで、上記構成を有するピッチ周波数推定装置１００におけるピッチ周波数推定動作について説明する。【００３４】  まず、ＦＦＴ部１０２では、次の式（２）で表される音声パワスペクトルＳＦ２（ｋ）を取得する。ここで、ｋは離散周波数成分を示す。ＨＦは、ピッチ周波数推定用の上限周波数成分であり、例えばＨＦ＝１［ｋＨｚ］である。Ｒｅ｛ＤＦ（ｋ）｝およびＩｍ｛ＤＦ（ｋ）｝は、それぞれＦＦＴ変換後の入力音声スペクトルＤＦ（ｋ）の実数部および虚数部を示す。
【数２】
2006006366.tif 000003
【００３５】  なお、式（２）では、スペクトルのパワ値を用いているが、パワ値の代わりに、平方根をとったスペクトル振幅値を用いても良い。【００３６】  また、有声性判定部１０３では、音声パワスペクトルＳＦ２（ｋ）の有声性を判定する。【００３７】  より具体的には、第１に、フレームｍの音声パワスペクトルＳＦ２（ｋ）の和Ｓ２（ｍ）と、推定雑音スペクトルパワの移動平均値Ｎ２（ｍ）と、を次の式（３）および（４）を用いてそれぞれ計算する。ここで、αは移動平均係数であり、ΘＮは、音声か雑音かを判定するための閾値である。
【数３】
2006006366.tif 000004
【数４】
2006006366.tif 000005
【００３８】  そして、第２に、音声と雑音との比ＳＮＲを式（５）を用いて計算し、その計算結果に基づいて有声性判定を行う。例えば式（６）に示すように、比ＳＮＲが閾値ΘＶよりも大きい場合は有声性ありと判定し、比ＳＮＲが閾値ΘＶ以下の場合は有声性なしと判定する。なお、ここでは有声性ありと判定された場合を例にとり、ピッチ周波数推定動作の説明を続ける。
【数５】
2006006366.tif 000006
【数６】
2006006366.tif 000007
【００３９】  そして、スペクトル抽出部１０４では、式（７）を用いて音声パワスペクトルＳＦ２（ｋ）のピークを抽出することにより、ピッチ調波スペクトルＰＦ（ｋ）の抽出を行う。
【数７】
2006006366.tif 000008
【００４０】  このとき、音声の準周期特性および雑音の影響により生じ得るピッチ調波スペクトルの位置ずれを考慮して、抽出されたピークの近傍にある音声パワスペクトルＳＦ２（ｋ−１）およびＳＦ２（ｋ＋１）を、ピッチ調波スペクトルＰＦ（ｋ−１）およびＰＦ（ｋ＋１）として一緒に抽出し、これら以外の周波数成分における音声パワスペクトルをゼロとみなす。【００４１】  また、スペクトル振幅制限部１０５で音声パワスペクトルの振幅制限が行われた場合、スペクトル抽出部１０４では、ピッチ調波スペクトルＰＦ（ｋ）にその振幅制限の結果を反映させることにより、ピッチ調波スペクトルＰＦ（ｋ）の振幅を制限する。【００４２】  すなわち、抽出されたピッチ調波スペクトルＰＦ（ｋ）を所定値と比較する。所定値は、周波数帯域ＨＦにおける音声パワスペクトルＳＦ２（ｋ）の平均値と乗算係数δとの積であり、式（８）によって求められる。そして、ピッチ調波スペクトルＰＦ（ｋ）が所定値を超過する場合には、式（９）を用いてピッチ調波スペクトルＰＦ（ｋ）の振幅に減衰係数を乗算することにより、ピッチ調波スペクトルＰＦ（ｋ）の振幅を制限する。減衰係数は式（１０）によって求められる。
【数８】
2006006366.tif 000009
【数９】
2006006366.tif 000010
【数１０】
2006006366.tif 000011
【００４３】  また、抽出されたピッチ調波スペクトルＰＦ（ｋ−１）およびＰＦ（ｋ＋１）に対しても同様に、式（１１）および（１２）を用いて振幅の制限を行う。
【数１１】
2006006366.tif 000012
【数１２】
2006006366.tif 000013
【００４４】  そして、スペクトル平均値計算部１０６では、式（１３）を用いて、ピッチ調波スペクトルＰＦ（ｋ）のパワの平均値ＰＡ（ｉ）を計算する。
【数１３】
2006006366.tif 000014
【００４５】  ここで、Ｎ（ｉ）＝ＮＦ／ｉであり、ＮＬ（ｉ）＝ｊ／ｉであり、ＮＨ（ｉ）＝（ＨＦ−ｊ）／ｉである。また、ｉはピッチ周波数候補であり、ＰＭＩＮおよびＰＭＡＸはそれぞれピッチ周波数候補の最小値および最大値である。また、ｊは、周波数帯域ＨＦにおける音声パワスペクトルＳＦ２（ｋ）の最大値に対応する周波数成分であり、ｎは、ピッチ周波数の整数倍の係数である。【００４６】  そして、スペクトル加算部１０７では、式（１４）を用いて、ピッチ調波スペクトルＰＦ（ｋ）のパワの加算値ＰＢ（ｉ）を計算する。
【数１４】
2006006366.tif 000015
【００４７】  ここで、式（１３）および（１４）を比較して分かるように、平均値ＰＡ（ｉ）および加算値ＰＢ（ｉ）の間には式（１５）で表される関係がある。したがって、スペクトル加算部１０７で式（１４）を用いて加算値ＰＢ（ｉ）を計算してから、スペクトル平均値計算部１０６で式（１３）の代わりに式（１５）を用いて平均値ＰＡ（ｉ）を計算した場合は、ピッチ周波数推定における演算量をさらに低減することができる。
【数１５】
2006006366.tif 000016
【００４８】  そして、べき乗計算部１０８では、例えば式（１６）を用いて、加算値ＰＢ（ｉ）のべき乗を計算する。
【数１６】
2006006366.tif 000017
【００４９】  そして、乗算部１０９では、式（１７）を用いて、べき乗計算結果ＰＣ（ｉ）を平均値ＰＡ（ｉ）に乗算する。
【数１７】
2006006366.tif 000018
【００５０】  そして、最大値抽出部１１０では、乗算結果ＰＤ（ｉ）の最大値ＰＤ−ｍａｘを抽出し、そのときのピッチ周波数候補ｐを推定ピッチ周波数として決定する。このようにしてピッチ周波数推定動作が行われる。【００５１】  続いて、半ピッチ周波数誤りおよび倍ピッチ周波数誤りの発生を防止するための条件（以下「防止条件」と言う）について説明する。ここでは、ピッチ調波スペクトルのパワの平均値のみを用いてピッチ周波数推定を行った場合（以下「第１のケース」と言う）と、ピッチ調波スペクトルのパワの平均値および加算値を用いてピッチ周波数推定を行った場合（以下「第２のケース」と言う）と、を例にとって説明する。【００５２】  まず、第１のケースにおける防止条件を定量的に求める。【００５３】  正しく推定されたピッチ周波数ｐに対する平均値ＰＡ（ｐ）を式（１８）で表した場合、半ピッチ周波数ｐ／２に対する平均値ＰＡ（ｐ／２）は式（１９）によって求められる。
【数１８】
2006006366.tif 000019
【数１９】
2006006366.tif 000020
【００５４】  ここで、ｘは、半ピッチ周波数ｐ／２を推定したときの、ピッチ周波数ｐに対する加算値ＰＢ（ｐ）の増加倍率を示す係数である。平均値ＰＡのみの最大化によりピッチ周波数を推定する場合、式（１８）および（１９）を比較して分かるように、ＰＡ（ｐ）＞ＰＡ（ｐ／２）つまりｘ＜１の条件を満たすときに、半ピッチ周波数誤りの発生を防止することができる。すなわち、加算値ＰＢの増加量がＰＢ（ｐ）未満のときに、半ピッチ周波数誤りの発生を防止することができる。【００５５】  また、倍ピッチ周波数２ｐに対する平均値ＰＡ（２ｐ）は式（２０）によって求められる。
【数２０】
2006006366.tif 000021
【００５６】  ここで、ｙは、倍ピッチ周波数２ｐを推定したときの、ピッチ周波数ｐに対する加算値ＰＢ（ｐ）の減少倍率を示す係数である。平均値ＰＡのみの最大化によりピッチ周波数を推定する場合、式（１８）および（２０）を比較して分かるように、ＰＡ（ｐ）＞ＰＡ（２ｐ）つまりｙ＞０．５の条件を満たすときに、倍ピッチ周波数誤りの発生を防止することができる。すなわち、加算値ＰＢの減少量が０．５ＰＢ（ｐ）より大きいときに、倍ピッチ周波数誤りの発生を防止することができる。【００５７】  次いで、第２のケースにおける防止条件を定量的に求める。【００５８】  前述の式（１７）で表される乗算結果ＰＤ（ｉ）を、半ピッチ周波数ｐ／２および倍ピッチ周波数２ｐに対してそれぞれ求めると、式（２１）および（２２）に示すとおりとなる。
【数２１】
2006006366.tif 000022
【数２２】
2006006366.tif 000023
【００５９】  式（１７）で表される乗算結果ＰＤ（ｉ）の最大化によってピッチ周波数を推定する場合、ＰＤ（ｐ）＞ＰＤ（ｐ／２）の条件を満たすときに、半ピッチ周波数誤りの発生を防止することができる。また、ＰＤ（ｐ）＞ＰＤ（２ｐ）の条件を満たすときに、倍ピッチ周波数誤りの発生を防止することができる。【００６０】  ここで、スペクトル抽出部１０４で抽出された音声パワスペクトルＳＦ２（ｋ）の例を図２Ａに示す。この例において、Ｐ２、Ｐ４、Ｐ５およびＰ６で示されるピークによりピッチ調波スペクトルが構成されると仮定する。【００６１】  また、図２Ｂに、加算値ＰＢ（ｉ）のべき乗の乗数を１に設定した条件の下で、平均値ＰＡ（ｉ）および加算値ＰＢ（ｉ）を互いに乗算した結果の例を示し、図２Ｃに、加算値ＰＢ（ｉ）のべき乗の乗数を３に設定した条件の下で、平均値ＰＡ（ｉ）および加算値ＰＢ（ｉ）を互いに乗算した結果の例を示す。【００６２】  そして、式（２１）を用いて半ピッチ周波数誤りの防止条件ＰＤ（ｐ）＞ＰＤ（ｐ／２）を変換すると、乗数が１の場合はｘ＜０．４１４となり、乗数が３の場合はｘ＜０．１８９となる。また、式（２２）を用いて倍ピッチ周波数誤りの防止条件ＰＤ（ｐ）＞ＰＤ（２ｐ）を変換すると、乗数が１の場合はｙ＞０．２９３となり、乗数が３の場合はｙ＞０．１５９となる。すなわち、乗数が１の場合は加算値ＰＢの増加量が０．４１４ＰＢ（ｐ）未満のときに、または、乗数が３の場合は加算値ＰＢの増加量が０．１８９ＰＢ（ｐ）未満のときに、半ピッチ周波数誤りの発生を防止することができる。また、乗数が１の場合は加算値ＰＢの減少量が０．２９３ＰＢ（ｐ）より大きいときに、または、乗数が３の場合は加算値ＰＢの減少量が０．１５９ＰＢ（ｐ）より大きいときに、倍ピッチ周波数誤りの発生を防止することができる。【００６３】  さらに、第１のケースにおける防止条件と第２のケースにおける防止条件とを比較する。この比較の結果として、倍ピッチ周波数誤りの防止条件については、第１のケースに比べて第２のケースの方が緩和されていることが分かる。すなわち、倍ピッチ周波数誤り発生の主因はホルマントによるピッチ調波スペクトル振幅値の変動であるが、この変動によって倍ピッチ周波数誤りの防止条件を満たさなくなる確率が、第１のケースよりも第２のケースの方が低くなる。したがって、ピッチ調波スペクトルのパワの平均値および加算値を用いてピッチ周波数推定を行うことにより、ホルマントの影響を低減することができ、ピッチ周波数推定の精度を向上することができる。【００６４】  さらに、べき乗の乗数を調整することによって、半ピッチ周波数誤りの発生率または倍ピッチ周波数誤りの発生率を自在に調整することができる。例えば、前述のとおり、乗数が１の場合と比べて乗数が３の場合は、半ピッチ周波数誤りが生じやすくなるが、倍ピッチ周波数誤りが生じにくくなる。逆に言えば、乗数が３の場合に比べて乗数が１の場合は、倍ピッチ周波数誤りが生じやすくなるが、半ピッチ周波数誤りが生じにくくなる。したがって、実際の場合は、音声や雑音の状態に応じて乗数を選択することによって、より正確にピッチ周波数を推定することができる。例えば、雑音の多い環境下でピッチ周波数推定が行われる場合は、乗数をより小さい値とすることによって、半ピッチ周波数誤りの発生率を低減することができる。一方、乗数をより大きい値とすることによって、ホルマントの影響による倍ピッチ周波数誤りの発生を低減することができる。【００６５】  ここで、同じ条件下で且つ同じピッチ調波スペクトルを用いてシミュレーションを行うことにより、式（１）で示される自己相関法に基づくピッチ周波数推定と本実施の形態に係るピッチ周波数推定との各推定誤り率を算出する。シミュレーションの諸条件は次のとおりである。ハニング窓長が３２０であり、ＦＦＴ変換長は５１２であり、移動平均係数αは０．０２であり、閾値ΘＶは２であり、乗算係数δは６であり、ピッチ周波数候補の最小値ＰＭＩＮは６２．５Ｈｚであり、ピッチ周波数候補の最大値ＰＭＡＸは３９０Ｈｚである。また、乗数βは３とした。下記の表は、算出された推定誤り率の一覧である。この表から分かるように、適切な乗数を選択することによって、本実施の形態に係るピッチ周波数推定は自己相関法に基づくものに比べて推定誤り率を低減することができる。
【表１】
2006006366.tif 000024
【００６６】  このように、本実施の形態によれば、ピッチ調波スペクトルのパワの平均値であって、複数のピッチ周波数候補の各々に対応づけて計算された平均値を用いて、ピッチ周波数を推定する、すなわち、周波数スペクトル上での自己相関を用いることなくピッチ周波数推定を行うため、ホルマントの影響を低減するためのスペクトル平坦化処理を不要とすることができるとともに、例えば、ピッチ調波スペクトルのパワに関する所定の定量的な条件が満たされる場合に半ピッチ周波数誤りや倍ピッチ周波数誤りの発生を防止することができ、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができる。【００６７】  また、本実施の形態によれば、ピッチ調波スペクトルのパワの平均値および加算値であって、複数のピッチ周波数候補の各々に対応づけて計算された平均値および加算値を、複数のピッチ周波数候補の各々に対応づけて互いに乗算し、乗算結果の最大値に対応するピッチ周波数候補を推定ピッチ周波数として決定する、すなわち、平均値および加算値の乗算値を関数としてピッチ周波数の推定を行うため、スペクトル平坦化処理を行うことなくホルマントの影響を低減することができ、ピッチ周波数推定の精度を向上することができる。【００６８】  なお、本実施の形態のピッチ周波数推定装置およびピッチ周波数推定方法は、音声符号化や音声強調などの音声信号処理を行う音声信号処理装置および音声信号処理方法に適用することができる。【００６９】  また、本発明は様々な実施の形態を採ることが可能であり、本実施の形態で説明したもののみに限定されない。例えば、上記のピッチ周波数推定方法をソフトウェアとしてコンピュータに実行させるようにしても良い。すなわち、上記の実施の形態で説明したピッチ周波数推定方法を実行するプログラムを予め例えばＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）等の記録媒体に記録しておき、そのプログラムをＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｏｒ  Ｕｎｉｔ）によって動作させることで、本発明のピッチ周波数推定方法を実行することができる。【００７０】  なお、上記各実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されても良いし、一部又は全てを含むように１チップ化されても良い。【００７１】  ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。【００７２】  また、集積回路化の手法はＬＳＩに限るものではなく、専用回路又は汎用プロセッサで実現しても良い。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサーを利用しても良い。【００７３】  さらには、半導体技術の進歩又は派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行っても良い。バイオ技術の適応等が可能性としてありえる。【００７４】  本明細書は、２００４年７月１３日出願の特願２００４−２０６３８７に基���く。この内容はすべてここに含めておく。
【産業上の利用可能性】
【００７５】  本発明のピッチ周波数推定装置およびピッチ周波数推定方法は、音声符号化や音声強調などの音声信号処理を行う装置および方法に適用することができる。
【図１】
2006006366.tif 000025
【図２Ａ】
2006006366.tif 000026
【図２Ｂ】
2006006366.tif 000027
【図２Ｃ】
2006006366.tif 000028
【手続補正書】
【提出日】平成18年5月12日(2006.5.12)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【書類名】請求の範囲
【請求項１】（削除）
【請求項２】（削除）
【請求項３】（削除）
【請求項４】（削除）
【請求項５】（削除）
【請求項６】（削除）
【請求項７】（削除）
【請求項８】（削除）
【請求項９】（削除）
【請求項１０】（削除）
【請求項１１】（削除）
【請求項１２】音声パワスペクトルからピッチ調波スペクトルを抽出する抽出手段と、
前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算する平均値計算手段と、
前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算する加算値計算手段と、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するべき乗計算手段と、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する決定手段と、
を具備するピッチ周波数推定装置。
【請求項１３】前記平均値計算手段は、周波数帯域中で最大の音声パワスペクトルを与える周波数を中心として、ピッチ周波数候補の整数倍を加減算した周波数における前記平均値を前記複数のピッチ周波数候補毎に計算する、
請求項１２記載のピッチ周波数推定装置。
【請求項１４】前記加算値計算手段は、周波数帯域中で最大の音声パワスペクトルを与える周波数を中心として、ピッチ周波数候補の整数倍を加減算した周波数における前記加算値を前記複数のピッチ周波数候補毎に計算する、
請求項１２記載のピッチ周波数推定装置。
【請求項１５】前記乗数が可変である、
請求項１２記載のピッチ周波数推定装置。
【請求項１６】前記ピッチ調波スペクトルの振幅を制限する振幅制限手段をさらに具備する、
請求項１２記載のピッチ周波数推定装置。
【請求項１７】前記音声パワスペクトルの有声性の有無を判定する判定手段をさらに具備し、
前記抽出手段は、前記有声性が有る場合は前記ピッチ調波スペクトルの抽出を行い、前記有声性が無い場合は前記ピッチ調波スペクトルの抽出を回避する、
請求項１２記載のピッチ周波数推定装置。
【請求項１８】前記判定手段は、音声と雑音との比であるＳＮＲと閾値との大小関係に従って前記音声パワスペクトルの有声性の有無を判定する、
請求項１７記載のピッチ周波数推定装置。
【請求項１９】音声パワスペクトルからピッチ調波スペクトルを抽出し、
前記ピッチ調波スペクトルのパワの平均値および加算値を複数のピッチ周波数候補毎に計算し、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算し、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、
前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する、
ピッチ周波数推定方法。
【請求項２０】音声パワスペクトルからピッチ調波スペクトルを抽出するステップと、
前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算するステップと、
前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算するステップと、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するステップと、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算するステップと、
前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定するステップと、
をコンピュータに実現させるためのピッチ周波数推定プログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００２
【補正方法】変更
【補正の内容】
【０００２】
ａｌ  ｆｒｅｑｕｅｎｃｙ  ｏｆ  ｎｏｉｓｅ−ｃｏｒｒｕｐｔｅｄ  ｓｐｅｅｃｈ”，Ｍ．Ｌａｈａｔ，ＩＥＥＥ  Ｔｒａｎｓ．ｏｎ  Ａｃｏｕｓｔｉｃｓ，Ｓｐｅｅｃｈ，ａｎｄ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ，ｖｏｌ．ＡＳＳＰ−３５，ｎｏ．６，ｐｐ．７４１−７５０，１９８７
【発明の開示】
【発明が解決しようとする課題】
［０００６］  しかしながら、上記従来のピッチ周波数推定方法では、スペクトルの平坦化処理を伴うため、ピッチ周波数推定に要する演算量が増大するという問題がある。
［０００７］  本発明の目的は、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができるピッチ周波数推定装置およびピッチ周波数推定方法を提供することである。
【課題を解決するための手段】
［０００８］（補正後）本発明のピッチ周波数推定装置は、音声パワスペクトルからピッチ調波スペクトルを抽出する抽出手段と、前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算する平均値計算手段と、前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算する加算値計算手段と、前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するべき乗計算手段と、前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する決定手段と、を具備する構成を採る。
［０００９］（削除）
［００１０］（削除）
【発明の効果】
［００１１］  本発明によれば、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができる。
【図面の簡単な説明】
２
【手続補正書】
【提出日】平成19年1月10日(2007.1.10)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
音声パワスペクトルからピッチ調波スペクトルを抽出する抽出手段と、
前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算する平均値計算手段と、
前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算する加算値計算手段と、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するべき乗計算手段と、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する決定手段と、
を具備するピッチ周波数推定装置。
【請求項２】
前記平均値計算手段は、周波数帯域中で最大の音声パワスペクトルを与える周波数を中心として、ピッチ周波数候補の整数倍を加減算した周波数における前記平均値を前記複数のピッチ周波数候補毎に計算する、
請求項１記載のピッチ周波数推定装置。
【請求項３】
前記加算値計算手段は、周波数帯域中で最大の音声パワスペクトルを与える周波数を中心として、ピッチ周波数候補の整数倍を加減算した周波数における前記加算値を前記複数のピッチ周波数候補毎に計算する、
請求項１記載のピッチ周波数推定装置。
【請求項４】
前記乗数が可変である、
請求項１記載のピッチ周波数推定装置。
【請求項５】
前記ピッチ調波スペクトルの振幅を制限する振幅制限手段をさらに具備する、
請求項１記載のピッチ周波数推定装置。
【請求項６】
前記音声パワスペクトルの有声性の有無を判定する判定手段をさらに具備し、
前記抽出手段は、前記有声性が有る場合は前記ピッチ調波スペクトルの抽出を行い、前記有声性が無い場合は前記ピッチ調波スペクトルの抽出を回避する、
請求項１記載のピッチ周波数推定装置。
【請求項７】
前記判定手段は、音声と雑音との比であるＳＮＲと閾値との大小関係に従って前記音声パワスペクトルの有声性の有無を判定する、
請求項６記載のピッチ周波数推定装置。
【請求項８】
音声パワスペクトルからピッチ調波スペクトルを抽出し、
前記ピッチ調波スペクトルのパワの平均値および加算値を複数のピッチ周波数候補毎に計算し、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算し、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、
前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する、
ピッチ周波数推定方法。
【請求項９】
音声パワスペクトルからピッチ調波スペクトルを抽出するステップと、
前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算するステップと、
前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算するステップと、
前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するステップと、
前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算するステップと、
前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定するステップと、
をコンピュータに実現させるためのピッチ周波数推定プログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ピッチ周波数推定装置およびピッチ周波数推定方法に関し、特に、周波数領域でピッチ周波数推定を行うピッチ周波数推定装置およびピッチ周波数推定方法に関する。
【背景技術】
【０００２】
一般に、時間領域または周波数領域において音声のピッチ周波数を推定する方法としては、音声波形の自己相関関数による自己相関法や、ＬＰＣ（Linear Predictive Coding）分析の残差信号の自己相関関数による変形相関法などが知られている。
【０００３】
また、雑音抑圧や音声符号化などの音声処理を周波数領域において行う場合は、周波数領域でピッチ周波数を推定すると整合性が良くなることがある。周波数領域でのピッチ周波数推定方法としては、周波数スペクトルに対する自己相関関数の最大化によりピッチ周波数を算出する方法があり、その一般式は次の式（１）によって表される。この式において、自己相関関数R(i)を最大にするピッチ周波数候補iを推定ピッチ周波数とする。
【数１】
2006006366.tif 000029
ここで、kは離散周波数成分であり、P(k)はピッチ調波スペクトルのパワであり、PMINおよびPMAXはそれぞれピッチ周波数候補iの最小値および最大値である。
【０００４】
ところで、周波数領域での自己相関を用いたピッチ周波数推定方法では、音声信号のホルマントの影響により誤って倍のピッチ周波数が算出されてしまうこと（倍ピッチ周波数誤り）がある。
【０００５】
ホルマントの影響を低減しつつピッチ周波数推定を行う従来の方法としては、例えば、非特許文献１に開示されたものがある。この方法では、スペクトル包絡の情報でスペクトルを平坦化した後のスペクトルが用いられる。
【非特許文献１】"A spectral autocorrelation method for measurement of the fundamental frequency of noise-corrupted speech", M. Lahat, IEEE Trans. on Acoustics, Speech, and Signal Processing, vol. ASSP-35, no. 6, pp. 741-750, 1987
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記従来のピッチ周波数推定方法では、スペクトルの平坦化処理を伴うため、ピッチ周波数推定に要する演算量が増大するという問題がある。
【０００７】
本発明の目的は、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができるピッチ周波数推定装置およびピッチ周波数推定方法を提供することである。
【課題を解決するための手段】
【０００８】
本発明のピッチ周波数推定装置は、音声パワスペクトルからピッチ調波スペクトルを抽出する抽出手段と、前記ピッチ調波スペクトルのパワの平均値を複数のピッチ周波数候補毎に計算する平均値計算手段と、前記ピッチ調波スペクトルのパワの加算値を前記複数のピッチ周波数候補毎に計算する加算値計算手段と、前記複数のピッチ周波数候補毎に、１以上のいずれかの乗数を用いて前記加算値のべき乗の値を計算するべき乗計算手段と、前記複数のピッチ周波数候補毎に前記べき乗の値を前記平均値に乗算し、前記複数のピッチ周波数候補のうち最大の乗算結果を与えるピッチ周波数候補を推定ピッチ周波数として決定する決定手段と、を有する構成を採る。
【発明の効果】
【０００９】
本発明によれば、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができる。
【発明を実施するための最良の形態】
【００１０】
以下、本発明の実施の形態について、図面を用いて詳細に説明する。
【００１１】
図１は、本発明の一実施の形態に係るピッチ周波数推定装置の構成を示すブロック図である。ピッチ周波数推定装置１００は、ハニング窓部１０１、ＦＦＴ（Fast Fourier Transform）部１０２、有声性判定部１０３、スペクトル抽出部１０４、スペクトル振幅制限部１０５、スペクトル平均値計算部１０６、スペクトル加算部１０７、べき乗計算部１０８、乗算部１０９および最大値抽出部１１０を有する。
【００１２】
ハニング窓部１０１は、所定時間単位のフレーム単位に分割された入力音声信号に対して、ハニングウィンドウなどを利用した窓掛け処理を施してＦＦＴ部１０２に出力する。
【００１３】
ＦＦＴ部１０２は、ハニング窓部１０１から入力されたフレーム、つまりフレーム単位に分割された音声信号に対してＦＦＴを行って音声信号を周波数領域に変換する。これにより、音声パワスペクトルを取得する。よって、フレーム単位の音声信号は、所定の周波数帯域を有する音声パワスペクトルとなる。このようにして生成された音声パワスペクトルは、有声性判定部１０３、スペクトル抽出部１０４およびスペクトル振幅制限部１０５に出力される。
【００１４】
有声性判定部１０３は、ＦＦＴ部１０２から音声パワスペクトルの有声性、つまり元の音声信号が有声であるか無声であるかを判定する。判定結果は、スペクトル抽出部１０４に出力される。
【００１５】
スペクトル抽出部１０４は、有声性判定部１０３によって音声パワスペクトルが有声性なしと判定された場合、ピッチ調波スペクトルの抽出を回避する。これによって、スペクトル抽出部１０４の演算量、ひいてはピッチ周波数推定装置１００の全体の演算量を低減することができる。
【００１６】
一方、音声パワスペクトルが有声性ありと判定された場合、スペクトル抽出部１０４はピッチ調波スペクトルの抽出を行う。より具体的には、音声パワスペクトルにおけるピークを抽出することにより、ピッチ調波スペクトルの抽出を行う。
【００１７】
また、スペクトル抽出部１０４は、スペクトル振幅制限部１０５による音声パワスペクトルの振幅制限が行われた場合、抽出されたピッチ調波スペクトルにその振幅制限の結果を反映させることにより、ピッチ調波スペクトルの振幅を制限する。このようにして、ピッチ周波数推定の精度に与えられ得るホルマントの影響を低減することができる。ピッチ調波スペクトルは、スペクトル平均値計算部１０６およびスペクトル加算部１０７に出力される。
【００１８】
スペクトル振幅制限部１０５は、ＦＦＴ部１０２によって取得された音声パワスペクトルの振幅が所定の閾値を超過しないように制限する。音声パワスペクトルの振幅制限の結果は、スペクトル抽出部１０４に出力される。
【００１９】
スペクトル平均値計算部１０６は、スペクトル抽出部１０４からのピッチ調波スペクトルのパワの平均値を、複数のピッチ周波数候補の各々に対応づけて計算する。すなわち、ピッチ調波スペクトルにおいて、ピッチ周波数候補を所定の最小値から所定の最大値までシフトさせながら、ピッチ周波数候補の整数倍に相当する周波数成分のパワの平均値を計算する。計算された平均値は、乗算部１０９に出力される。
【００２０】
また、スペクトル平均値計算部１０６は、平均値の計算を行うとき、パワの最大値に対応する周波数成分を、平均値計算対象の周波数帯域における基準周波数として用いる。
【００２１】
具体的には、基準周波数からピッチ周波数候補の整数倍に相当する周波数を減算して得られる周波数におけるパワと、基準周波数からピッチ周波数候補の整数倍に相当する周波数を加算して得られる周波数におけるパワと、を用いて、平均値の計算を行う。これにより、音声の準周期特性および雑音の影響ならびにピッチ周波数推定誤差により生じるピッチ高調波における誤差の累積を低減することができ、より正確にピッチ周波数の推定を行うことができる。
【００２２】
なお、ピッチ調波スペクトルのパワの平均値は、後述するピッチ調波スペクトルのパワの加算値を特定の値で除して得られる値である。よって、スペクトル平均値計算部１０６は、スペクトル加算部１０７によって計算された加算値を取得し、これを用いて平均値の算出を行っても良い。
【００２３】
スペクトル加算部１０７は、スペクトル抽出部１０４からのピッチ調波スペクトルのパワの加算値を、複数のピッチ周波数候補の各々に対応づけて計算する。すなわち、ピッチ調波スペクトルにおいて、ピッチ周波数候補を所定の最小値から所定の最大値までシフトさせながら、ピッチ周波数候補の整数倍に相当する周波数成分のパワを加算する。パワの加算によって得られた加算値はべき乗計算部１０８に出力される。
【００２４】
また、スペクトル加算部１０７は、パワの加算を行うとき、パワの最大値に対応する周波数成分を、加算値計算対象の周波数帯域における基準周波数として用いる。
【００２５】
具体的には、基準周波数からピッチ周波数候補の整数倍に相当する周波数を減算して得られる周波数におけるパワと、基準周波数からピッチ周波数候補の整数倍に相当する周波数を加算して得られる周波数におけるパワと、を用いて、加算値の計算を行う。これにより、音声の準周期特性および雑音の影響ならびにピッチ周波数推定誤差により生じるピッチ高調波における誤差の累積を低減することができ、より正確にピッチ周波数の推定を行うことができる。
【００２６】
べき乗計算部１０８は、スペクトル加算部１０７によって算出された加算値のべき乗の値を計算する。算出されたべき乗の値は乗算部１０９に出力される。また、べき乗計算部１０８は、べき乗の計算に用いられる乗数を可変に設定する。乗数の可変設定つまり乗数の調整については後述する。
【００２７】
乗算部１０９および最大値抽出部１１０の組み合わせは、複数のピッチ周波数候補の各々に対応づけて計算された平均値を用いてピッチ周波数を推定する推定部を構成する。
【００２８】
推定部において、乗算部１０９は、ピッチ調波スペクトルのパワの平均値とピッチ調波スペクトルのパワの加算値とを、複数のピッチ周波数候補の各々に対応づけて乗算する。より具体的には、加算値のべき乗計算結果を平均値に乗算する。乗算結果は、最大値抽出部１１０に出力される。
【００２９】
最大値抽出部１１０は、乗算部１０９で計算された乗算結果の最大値を抽出する。また、所定の最小値から所定の最大値までの複数のピッチ周波数候補のうち、乗算結果が最大となるときのピッチ周波数候補を推定ピッチ周波数として決定し、図示されない後段の処理部に出力する。
【００３０】
次いで、上記構成を有するピッチ周波数推定装置１００におけるピッチ周波数推定動作について説明する。
【００３１】
まず、ＦＦＴ部１０２では、次の式（２）で表される音声パワスペクトルSF2(k)を取得する。ここで、kは離散周波数成分を示す。HFは、ピッチ周波数推定用の上限周波数成分であり、例えばHF=1[kHz]である。Re{DF(k)}およびIm{DF(k)}は、それぞれＦＦＴ変換後の入力音声スペクトルDF(k)の実数部および虚数部を示す。
【数２】
2006006366.tif 000030
【００３２】
なお、式（２）では、スペクトルのパワ値を用いているが、パワ値の代わりに、平方根をとったスペクトル振幅値を用いても良い。
【００３３】
また、有声性判定部１０３では、音声パワスペクトルSF2(k)の有声性を判定する。
【００３４】
より具体的には、第１に、フレームmの音声パワスペクトルSF2(k)の和S2(m)と、推定雑音スペクトルパワの移動平均値N2(m)と、を次の式（３）および（４）を用いてそれぞれ計算する。ここで、αは移動平均係数であり、ΘNは、音声か雑音かを判定するための閾値である。
【数３】
2006006366.tif 000031
【数４】
2006006366.tif 000032
【００３５】
そして、第２に、音声と雑音との比SNRを式（５）を用いて計算し、その計算結果に基づいて有声性判定を行う。例えば式（６）に示すように、比SNRが閾値ΘVよりも大きい場合は有声性ありと判定し、比SNRが閾値ΘV以下の場合は有声性なしと判定する。なお、ここでは有声性ありと判定された場合を例にとり、ピッチ周波数推定動作の説明を続ける。
【数５】
2006006366.tif 000033
【数６】
2006006366.tif 000034
【００３６】
そして、スペクトル抽出部１０４では、式（７）を用いて音声パワスペクトルSF2(k)のピークを抽出することにより、ピッチ調波スペクトルPF(k)の抽出を行う。
【数７】
2006006366.tif 000035
【００３７】
このとき、音声の準周期特性および雑音の影響により生じ得るピッチ調波スペクトルの位置ずれを考慮して、抽出されたピークの近傍にある音声パワスペクトルSF2(k-1)およびSF2(k+1)を、ピッチ調波スペクトルPF(k-1)およびPF(k+1)として一緒に抽出し、これら以外の周波数成分における音声パワスペクトルをゼロとみなす。
【００３８】
また、スペクトル振幅制限部１０５で音声パワスペクトルの振幅制限が行われた場合、スペクトル抽出部１０４では、ピッチ調波スペクトルPF(k)にその振幅制限の結果を反映させることにより、ピッチ調波スペクトルPF(k)の振幅を制限する。
【００３９】
すなわち、抽出されたピッチ調波スペクトルPF(k)を所定値と比較する。所定値は、周波数帯域HFにおける音声パワスペクトルSF2(k)の平均値と乗算係数δとの積であり、式（８）によって求められる。そして、ピッチ調波スペクトルPF(k)が所定値を超過する場合には、式（９）を用いてピッチ調波スペクトルPF(k)の振幅に減衰係数を乗算することにより、ピッチ調波スペクトルPF(k)の振幅を制限する。減衰係数は式（１０）によって求められる。
【数８】
2006006366.tif 000036
【数９】
2006006366.tif 000037
【数１０】
2006006366.tif 000038
【００４０】
また、抽出されたピッチ調波スペクトルPF(k-1)およびPF(k+1)に対しても同様に、式（１１）および（１２）を用いて振幅の制限を行う。
【数１１】
2006006366.tif 000039
【数１２】
2006006366.tif 000040
【００４１】
そして、スペクトル平均値計算部１０６では、式（１３）を用いて、ピッチ調波スペクトルPF(k)のパワの平均値PA(i)を計算する。
【数１３】
2006006366.tif 000041
【００４２】
ここで、N(i)=NF/iであり、NL(i)=j/iであり、NH(i)=(HF-j)/iである。また、iはピッチ周波数候補であり、PMINおよびPMAXはそれぞれピッチ周波数候補の最小値および最大値である。また、jは、周波数帯域HFにおける音声パワスペクトルSF2(k)の最大値に対応する周波数成分であり、nは、ピッチ周波数の整数倍の係数である。
【００４３】
そして、スペクトル加算部１０７では、式（１４）を用いて、ピッチ調波スペクトルPF(k)のパワの加算値PB(i)を計算する。
【数１４】
2006006366.tif 000042
【００４４】
ここで、式（１３）および（１４）を比較して分かるように、平均値PA(i)および加算値PB(i)の間には式（１５）で表される関係がある。したがって、スペクトル加算部１０７で式（１４）を用いて加算値PB(i)を計算してから、スペクトル平均値計算部１０６で式（１３）の代わりに式（１５）を用いて平均値PA(i)を計算した場合は、ピッチ周波数推定における演算量をさらに低減することができる。
【数１５】
2006006366.tif 000043
【００４５】
そして、べき乗計算部１０８では、例えば式（１６）を用いて、加算値PB(i)のべき乗を計算する。
【数１６】
2006006366.tif 000044
【００４６】
そして、乗算部１０９では、式（１７）を用いて、べき乗計算結果PC(i)を平均値PA(i)に乗算する。
【数１７】
2006006366.tif 000045
【００４７】
そして、最大値抽出部１１０では、乗算結果PD(i)の最大値PD_maxを抽出し、そのときのピッチ周波数候補pを推定ピッチ周波数として決定する。このようにしてピッチ周波数推定動作が行われる。
【００４８】
続いて、半ピッチ周波数誤りおよび倍ピッチ周波数誤りの発生を防止するための条件（以下「防止条件」と言う）について説明する。ここでは、ピッチ調波スペクトルのパワの平均値のみを用いてピッチ周波数推定を行った場合（以下「第１のケース」と言う）と、ピッチ調波スペクトルのパワの平均値および加算値を用いてピッチ周波数推定を行った場合（以下「第２のケース」と言う）と、を例にとって説明する。
【００４９】
まず、第１のケースにおける防止条件を定量的に求める。
【００５０】
正しく推定されたピッチ周波数pに対する平均値PA(p)を式（１８）で表した場合、半ピッチ周波数p/2に対する平均値PA(p/2)は式（１９）によって求められる。
【数１８】
2006006366.tif 000046
【数１９】
2006006366.tif 000047
【００５１】
ここで、xは、半ピッチ周波数p/2を推定したときの、ピッチ周波数pに対する加算値PB(p)の増加倍率を示す係数である。平均値PAのみの最大化によりピッチ周波数を推定する場合、式（１８）および（１９）を比較して分かるように、PA(p)>PA(p/2)つまりx<1の条件を満たすときに、半ピッチ周波数誤りの発生を防止することができる。すなわち、加算値PBの増加量がPB(p)未満のときに、半ピッチ周波数誤りの発生を防止することができる。
【００５２】
また、倍ピッチ周波数2pに対する平均値PA(2p)は式（２０）によって求められる。
【数２０】
2006006366.tif 000048
【００５３】
ここで、yは、倍ピッチ周波数2pを推定したときの、ピッチ周波数pに対する加算値PB(p)の減少倍率を示す係数である。平均値PAのみの最大化によりピッチ周波数を推定する場合、式（１８）および（２０）を比較して分かるように、PA(p)>PA(2p)つまりy>0.5の条件を満たすときに、倍ピッチ周波数誤りの発生を防止することができる。すなわち、加算値PBの減少量が0.5PB(p)より大きいときに、倍ピッチ周波数誤りの発生を防止することができる。
【００５４】
次いで、第２のケースにおける防止条件を定量的に求める。
【００５５】
前述の式（１７）で表される乗算結果PD(i)を、半ピッチ周波数p/2および倍ピッチ周波数2pに対してそれぞれ求めると、式（２１）および（２２）に示すとおりとなる。
【数２１】
2006006366.tif 000049
【数２２】
2006006366.tif 000050
【００５６】
式（１７）で表される乗算結果PD(i)の最大化によってピッチ周波数を推定する場合、PD(p)> PD(p/2)の条件を満たすときに、半ピッチ周波数誤りの発生を防止することができる。また、PD(p)> PD(2p)の条件を満たすときに、倍ピッチ周波数誤りの発生を防止することができる。
【００５７】
ここで、スペクトル抽出部１０４で抽出された音声パワスペクトルSF2(k)の例を図２Ａに示す。この例において、P2、P4、P5およびP6で示されるピークによりピッチ調波スペクトルが構成されると仮定する。
【００５８】
また、図２Ｂに、加算値PB(i)のべき乗の乗数を１に設定した条件の下で、平均値PA(i)および加算値PB(i)を互いに乗算した結果の例を示し、図２Ｃに、加算値PB(i)のべき乗の乗数を３に設定した条件の下で、平均値PA(i)および加算値PB(i)を互いに乗算した結果の例を示す。
【００５９】
そして、式（２１）を用いて半ピッチ周波数誤りの防止条件PD(p)> PD(p/2)を変換すると、乗数が１の場合はx<0.414となり、乗数が３の場合はx<0.189となる。また、式（２２）を用いて倍ピッチ周波数誤りの防止条件PD(p)> PD(2p)を変換すると、乗数が１の場合はy>0.293となり、乗数が３の場合はy>0.159となる。すなわち、乗数が１の場合は加算値PBの増加量が0.414PB(p)未満のときに、または、乗数が３の場合は加算値PBの増加量が0.189PB(p)未満のときに、半ピッチ周波数誤りの発生を防止することができる。また、乗数が１の場合は加算値PBの減少量が0.293PB(p)より大きいときに、または、乗数が３の場合は加算値PBの減少量が0.159PB(p)より大きいときに、倍ピッチ周波数誤りの発生を防止することができる。
【００６０】
さらに、第１のケースにおける防止条件と第２のケースにおける防止条件とを比較する。この比較の結果として、倍ピッチ周波数誤りの防止条件については、第１のケースに比べて第２のケースの方が緩和されていることが分かる。すなわち、倍ピッチ周波数誤り発生の主因はホルマントによるピッチ調波スペクトル振幅値の変動であるが、この変動によって倍ピッチ周波数誤りの防止条件を満たさなくなる確率が、第１のケースよりも第２のケースの方が低くなる。したがって、ピッチ調波スペクトルのパワの平均値および加算値を用いてピッチ周波数推定を行うことにより、ホルマントの影響を低減することができ、ピッチ周波数推定の精度を向上することができる。
【００６１】
さらに、べき乗の乗数を調整することによって、半ピッチ周波数誤りの発生率または倍ピッチ周波数誤りの発生率を自在に調整することができる。例えば、前述のとおり、乗数が１の場合と比べて乗数が３の場合は、半ピッチ周波数誤りが生じやすくなるが、倍ピッチ周波数誤りが生じにくくなる。逆に言えば、乗数が３の場合に比べて乗数が１の場合は、倍ピッチ周波数誤りが生じやすくなるが、半ピッチ周波数誤りが生じにくくなる。したがって、実際の場合は、音声や雑音の状態に応じて乗数を選択することによって、より正確にピッチ周波数を推定することができる。例えば、雑音の多い環境下でピッチ周波数推定が行われる場合は、乗数をより小さい値とすることによって、半ピッチ周波数誤りの発生率を低減することができる。一方、乗数をより大きい値とすることによって、ホルマントの影響による倍ピッチ周波数誤りの発生を低減することができる。
【００６２】
ここで、同じ条件下で且つ同じピッチ調波スペクトルを用いてシミュレーションを行うことにより、式（１）で示される自己相関法に基づくピッチ周波数推定と本実施の形態に係るピッチ周波数推定との各推定誤り率を算出する。シミュレーションの諸条件は次のとおりである。ハニング窓長が３２０であり、ＦＦＴ変換長は５１２であり、移動平均係数αは０．０２であり、閾値ΘVは２であり、乗算係数δは６であり、ピッチ周波数候補の最小値PMINは６２．５Ｈｚであり、ピッチ周波数候補の最大値PMAXは３９０Ｈｚである。また、乗数βは３とした。下記の表は、算出された推定誤り率の一覧である。この表から分かるように、適切な乗数を選択することによって、本実施の形態に係るピッチ周波数推定は自己相関法に基づくものに比べて推定誤り率を低減することができる。
【表１】
2006006366.tif 000051
【００６３】
このように、本実施の形態によれば、ピッチ調波スペクトルのパワの平均値であって、複数のピッチ周波数候補の各々に対応づけて計算された平均値を用いて、ピッチ周波数を推定する、すなわち、周波数スペクトル上での自己相関を用いることなくピッチ周波数推定を行うため、ホルマントの影響を低減するためのスペクトル平坦化処理を不要とすることができるとともに、例えば、ピッチ調波スペクトルのパワに関する所定の定量的な条件が満たされる場合に半ピッチ周波数誤りや倍ピッチ周波数誤りの発生を防止することができ、ピッチ周波数推定に要する演算量を低減しつつ、ピッチ周波数を正確に推定することができる。
【００６４】
また、本実施の形態によれば、ピッチ調波スペクトルのパワの平均値および加算値であって、複数のピッチ周波数候補の各々に対応づけて計算された平均値および加算値を、複数のピッチ周波数候補の各々に対応づけて互いに乗算し、乗算結果の最大値に対応するピッチ周波数候補を推定ピッチ周波数として決定する、すなわち、平均値および加算値の乗算値を関数としてピッチ周波数の推定を行うため、スペクトル平坦化処理を行うことなくホルマントの影響を低減することができ、ピッチ周波数推定の精度を向上することができる。
【００６５】
なお、本実施の形態のピッチ周波数推定装置およびピッチ周波数推定方法は、音声符号化や音声強調などの音声信号処理を行う音声信号処理装置および音声信号処理方法に適用することができる。
【００６６】
また、本発明は様々な実施の形態を採ることが可能であり、本実施の形態で説明したもののみに限定されない。例えば、上記のピッチ周波数推定方法をソフトウェアとしてコンピュータに実行させるようにしても良い。すなわち、上記の実施の形態で説明したピッチ周波数推定方法を実行するプログラムを予め例えばＲＯＭ（Read Only Memory）等の記録媒体に記録しておき、そのプログラムをＣＰＵ（Central Processor Unit）によって動作させることで、本発明のピッチ周波数推定方法を実行することができる。
【００６７】
なお、上記各実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されても良いし、一部又は全てを含むように１チップ化されても良い。
【００６８】
ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
【００６９】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路又は汎用プロセッサで実現しても良い。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサーを利用しても良い。
【００７０】
さらには、半導体技術の進歩又は派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行っても良い。バイオ技術の適応等が可能性としてありえる。
【００７１】
本明細書は、２００４年７月１３日出願の特願２００４−２０６３８７に基づく。この内容はすべてここに含めておく。
【産業上の利用可能性】
【００７２】
本発明のピッチ周波数推定装置およびピッチ周波数推定方法は、音声符号化や音声強調などの音声信号処理を行う装置および方法に適用することができる。
【図面の簡単な説明】
【００７３】
【図１】本発明の一実施の形態に係るピッチ周波数推定装置の構成を示すブロック図
【図２Ａ】本発明の一実施の形態において、抽出された音声パワスペクトルの例を示す図
【図２Ｂ】本発明の一実施の形態において、乗数をある値に設定した条件の下で平均値および加算値を乗算した結果を示す図
【図２Ｃ】本発明の一実施の形態において、乗数を他の値に設定した条件の下で平均値および加算値を乗算した結果を示す図
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】全図
【補正方法】変更
【補正の内容】
【図１】
2006006366.tif 000052
【図２Ａ】
2006006366.tif 000053
【図２Ｂ】
2006006366.tif 000054
【図２Ｃ】
2006006366.tif 000055
【国際調査報告】
2006006366.tif 000056
2006006366.tif 000057
2006006366.tif 000058
2006006366.tif 000059
2006006366.tif 000060