(57)【要約】【課題】従来よりも分析精度が高い線形予測分析技術を提供することを目的とする。【解決手段】自己相関計算部２１が、入力信号から自己相関RO(i)を計算する。予測係数計算部２３が、係数wO(i)と自己相関RO(i)とが乗算されたものである変形自己相関R'O(i)を用いて線形予測分析を行う。ここで、少なくとも一部の各次数iに対して、各次数iに対応する係数wO(i)が、現在又は過去のフレームにおける入力信号のピッチゲインと正の相関関係にある値の増加とともに単調減少する関係にある場合が含まれているとする。係数は、周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、互いに異なる２以上の係数テーブルから取得される。【選択図】図１
2018028699.tif 000002
【特許請求の範囲】
【請求項１】
入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析方法であって、
少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i)を計算する自己相関計算ステップと、
係数と前記自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R' O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算ステップと、を含み、
係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt1には係数wt1(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて前記係数テーブルt0,t1,t2の中の１個の係数テーブルから係数を取得する係数決定ステップを更に含み、
前記周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に前記係数決定ステップで係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが中程度の場合に前記係数決定ステップで係数が取得される係数テーブルを係数テーブルt1とし、周期性の強さ又はピッチゲインが小さい場合に前記係数決定ステップで係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)<wt1(i)≦wt2(i)であり、i=0以外のiのうちの前記以外のiのうちの少なくとも一部の各iについてwt0(i)≦wt1(i)<wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt1(i)≦wt2(i)である、
線形予測分析方法。
【請求項２】
入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析方法であって、
少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i) を計算する自己相関計算ステップと、
係数と前記自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R'O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算ステップと、を含み、
係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて前記係数テーブルt0,t2の少なくとも何れかの係数テーブルから係数を取得する係数決定ステップを更に含み、
前記周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に前記係数決定ステップで係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが小さい場合に前記係数決定ステップで係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)＜wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt2(i)であり、
前記係数決定ステップは、周期性の強さ又はピッチゲインが中程度の場合には、i=0以外の各iについて、ｗ0(i)=β’×wt0(i)+(1-β’)×wt2(i) （ただし0≦β’≦1）である係数を決定する、
線形予測分析方法。
【請求項３】
入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析装置であって、
少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i)を計算する自己相関計算部と、
係数と前記自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R' O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算部と、を含み、
係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt1には係数wt1(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて前記係数テーブルt0,t1,t2の中の１個の係数テーブルから係数を取得する係数決定部を更に含み、
前記周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に前記係数決定部で係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが中程度の場合に前記係数決定部で係数が取得される係数テーブルを係数テーブルt1とし、周期性の強さ又はピッチゲインが小さい場合に前記係数決定部で係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)<wt1(i)≦wt2(i)であり、i=0以外のiのうちの前記以外のiのうちの少なくとも一部の各iについてwt0(i)≦wt1(i)<wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt1(i)≦wt2(i)である、
線形予測分析装置。
【請求項４】
入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析装置であって、
少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i) を計算する自己相関計算部と、
係数と前記自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R'O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算部と、を含み、
係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて前記係数テーブルt0,t2の少なくとも何れかの係数テーブルから係数を取得する係数決定部を更に含み、
前記周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に前記係数決定部で係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが小さい場合に前記係数決定部で係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)＜wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt2(i)であり、
前記係数決定部は、周期性の強さ又はピッチゲインが中程度の場合には、i=0以外の各iについて、ｗ0(i)=β’×wt0(i)+(1-β’)×wt2(i) （ただし0≦β’≦1）である係数を決定する、
線形予測分析装置。
【請求項５】
請求項１又は２の何れかの線形予測分析方法の各ステップをコンピュータに実行させるためのプログラム。
【請求項６】
請求項１又は２の何れかの線形予測分析方法の各ステップをコンピュータに実行させるためのプログラムが記録されたコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、音声信号、音響信号、心電図、脳波、脳磁図、地震波等のディジタル時系列信号の分析技術に関する。
【背景技術】
【０００２】
音声信号、音響信号の符号化では、入力された音声信号や音響信号を線形予測分析して得た予測係数に基づいて符号化する手法が広く用いられている（例えば、非特許文献１，２参照。）。
【０００３】
非特許文献１から３では、図１１に例示する線形予測分析装置により予測係数が計算されている。線形予測分析装置１は、自己相関計算部１１、係数乗算部１２及び予測係数計算部１３を備えている。
【０００４】
入力された時間領域のディジタル音声信号やディジタル音響信号である入力信号は、Nサンプルのフレーム毎に処理される。現時刻で処理対象とするフレームである現フレームの入力信号をXO(n)(n=0,1,…,N-1)とする。nは入力信号における各サンプルのサンプル番号を表し、Nは所定の正の整数である。ここで、現フレームの１つ前のフレームの入力信号はXO(n)(n=-N,-N+1,…,-1)であり、現フレームの１つ後のフレームの入力信号はXO(n)(n=N,N+1,…,2N-1)である。
【０００５】
[自己相関計算部１１]
線形予測分析装置１の自己相関計算部１１は、入力信号XO(n)から自己相関RO(i) (i=0,1,…,Pmax,Pmaxは予測次数)を式(11)により求めて出力する。Pmaxは、N未満の所定の正の整数である。
【数１】
2018028699.tif 000003
【０００６】
［係数乗算部１２］
次に、係数乗算部１２が、自己相関計算部１１から出力された自己相関RO(i)に予め定めた係数wO(i) (i=0,1,…,Pmax)を同じiごとに乗じることにより、変形自己相関R'O(i) (i=0,1,…,Pmax)を求める。すなわち、変形自己相関関数R' O(i)を式(12)により求める。
【数２】
2018028699.tif 000004
【０００７】
［予測係数計算部１３］
そして、予測係数計算部１３が、係数乗算部１２から出力された変形自己相関R'O(i)を用いて例えばLevinson-Durbin法などにより、1次から予め定めた予測次数であるPmax次までの線形予測係数に変換可能な係数を求める。線形予測係数に変換可能な係数とは、PARCOR係数KO(1),KO(2),…,KO(Pmax)や線形予測係数aO(1),aO(2),…,aO(Pmax)等である。
【０００８】
非特許文献１である国際標準ITU-T G.718や非特許文献２である国際標準ITU-T G.729では、係数wO(i)として予め求めておいた60 Hzのバンド幅の固定の係数を用いている。
【０００９】
具体的には、係数wO(i)は式(13)のように指数関数を用いて定義されており、式(13)の中ではf0=60 Hzという固定値が使われている。fsはサンプリング周波数である。
【数３】
2018028699.tif 000005
【００１０】
非特許文献３には、上述の指数関数以外の関数に基づく係数を用いる例が記載されている。しかし、ここで用いられている関数は、サンプリング周期τ(fsに対応する周期に相当）と所定の定数aとに基づく関数であり、やはり固定値の係数が使われている。
【先行技術文献】
【非特許文献】
【００１１】
【非特許文献１】ITU-T Recommendation G.718, ITU, 2008.
【非特許文献２】ITU-T Recommendation G.729, ITU, 1996
【非特許文献３】Yoh'ichi Tohkura, Fumitada Itakura, Shin'ichiro Hashimoto, "Spectral Smoothing Technique in PARCOR Speech Analysis-Synthesis", IEEE Trans. on Acoustics, Speech, and Signal Processing, Vol. ASSP-26, No.6, 1978
【発明の概要】
【発明が解決しようとする課題】
【００１２】
従来の音声信号、音響信号の符号化で用いられている線形予測分析方法では、自己相関RO(i)に固定の係数wO(i)が乗算して得られる変形自己相関R'O(i)を用いて線形予測係数に変換可能な係数を求めていた。よって、自己相関RO(i)への係数wO(i)の乗算による変形を必要としないような、すなわち、変形自己相関R'O(i)ではなく自己相関RO(i)そのものを用いて線形予測係数に変換可能な係数を求めたとしても、線形予測係数に変換可能な係数に対応するスペクトル包絡においてスペクトルのピークが大きくなりすぎることがないような入力信号の場合には、自己相関RO(i)への係数wO(i)の乗算によって、変形自己相関R'O(i)により求まる線形予測係数に変換可能な係数に対応するスペクトル包絡が、入力信号XO(n)のスペクトル包絡を近似する精度が下がってしまう、すなわち、線形予測分析の精度が下がってしまう可能性、があった。
【００１３】
この発明は、従来よりも分析精度が高い線形予測分析方法、装置、プログラム及び記録媒体を提供することを目的とする。
【課題を解決するための手段】
【００１４】
この発明の一態様による線形予測分析方法は、入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析方法であって、少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i)を計算する自己相関計算ステップと、係数と自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R' O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算ステップと、を含み、係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt1には係数wt1(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて係数テーブルt0,t1,t2の中の１個の係数テーブルから係数を取得する係数決定ステップを更に含み、周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に係数決定ステップで係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが中程度の場合に係数決定ステップで係数が取得される係数テーブルを係数テーブルt1とし、周期性の強さ又はピッチゲインが小さい場合に係数決定ステップで係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)<wt1(i)≦wt2(i)であり、i=0以外のiのうちの少なくとも一部の各iについてwt0(i)≦wt1(i)<wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt1(i)≦wt2(i)である。
【００１５】
この発明の一態様による線形予測分析方法は、入力時系列信号に対応する線形予測係数に変換可能な係数を、所定時間区間であるフレームごとに求める、線形予測分析方法であって、少なくともi=0,1,…,Pmaxのそれぞれについて、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)またはiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i) を計算する自己相関計算ステップと、係数と自己相関RO(i)とが対応するiごとに乗算されたものである変形自己相関R' O(i)を用いて、1次からPmax次までの線形予測係数に変換可能な係数を求める予測係数計算ステップと、を含み、係数テーブルt0には係数wt0(i)が格納されており、係数テーブルt2には係数wt2(i)が格納されているとして、現在又は過去のフレームにおける入力時系列信号の周期性の強さ又は入力時系列信号に基づくピッチゲインと正の相関関係にある値を用いて係数テーブルt0,t2の少なくとも何れかの係数テーブルから係数を取得する係数決定ステップを更に含み、周期性の強さ又はピッチゲインと正の相関関係にある値に応じて、周期性の強さ又はピッチゲインが大きい場合、周期性の強さ又はピッチゲインが中程度の場合、周期性の強さ又はピッチゲインが小さい場合の何れかの場合に分類されるとし、周期性の強さ又はピッチゲインが大きい場合に係数決定ステップで係数が取得される係数テーブルを係数テーブルt0とし、周期性の強さ又はピッチゲインが小さい場合に係数決定ステップで係数が取得される係数テーブルを係数テーブルt2として、i=0以外のiのうちの少なくとも一部のiについてwt0(i)＜wt2(i)であり、i=0以外のiのうちの残りの各iについてwt0(i)≦wt2(i)であり、上記係数決定ステップは、周期性の強さ又はピッチゲインが中程度の場合には、i=0以外の各iについて、ｗ0(i)=β’×wt0(i)+(1-β’)×wt2(i) （ただし0≦β’≦1）である係数を決定する。
【発明の効果】
【００１６】
従来よりも分析精度の高い線形予測を実現することができる。
【図面の簡単な説明】
【００１７】
【図１】第一実施形態及び第二実施形態の線形予測装置の例を説明するためのブロック図。
【図２】線形予測分析方法の例を説明するためのフローチャート。
【図３】第二実施形態の線形予測分析方法の例を説明するためのフローチャート。
【図４】第三実施形態の線形予測装置の例を説明するためのブロック図。
【図５】第三実施形態の線形予測分析方法の例を説明するためのフローチャート。
【図６】第三実施形態の具体例を説明するための図。
【図７】変形例を説明するためのブロック図。
【図８】変形例を説明するためのブロック図。
【図９】変形例を説明するためのフローチャート。
【図１０】第四実施形態の線形予測分析装置の例を説明するためのブロック図。
【図１１】従来の線形予測装置の例を説明するためのブロック図。
【発明を実施するための形態】
【００１８】
以下、図面を参照して、線形予測分析装置及び方法の各実施形態を説明する。
【００１９】
［第一実施形態］
第一実施形態の線形予測分析装置２は、図１に示すように、自己相関計算部２１、係数決定部２４、係数乗算部２２及び予測係数計算部２３を例えば備えている。自己相関計算部２１、係数乗算部２２及び予測係数計算部２３の動作は、従来の線形予測分析装置１の自己相関計算部１１、係数乗算部１２及び予測係数計算部１３における動作とそれぞれ同じである。
【００２０】
線形予測分析装置２には、所定時間区間であるフレームごとの時間領域のディジタル音声信号やディジタル音響信号や心電図、脳波、脳磁図、地震波等のやディジタル信号である入力信号XO(n)が入力される。入力信号は、入力時系列信号である。現フレームの入力信号をXO(n)(n=0,1,…,N-1)とする。nは入力信号における各サンプルのサンプル番号を表し、Nは所定の正の整数である。ここで、現フレームの１つ前のフレームの入力信号はXO(n)(n=-N,-N+1,…, -1)であり、現フレームの１つ後のフレームの入力信号はXO(n)(n=N,N+1,…, 2N-1)である。以下では、入力信号XO(n)がディジタル音声信号やディジタル音響信号である場合について説明する。入力信号XO(n) (n=0,1,…,N-1)は、収音された信号そのものであってもよいし、分析のためにサンプリングレートが変換された信号でもよいし、プリエンファシス処理された信号でもよいし、窓かけされた信号でもよい。
【００２１】
また、線形予測分析装置２には、フレームごとのディジタル音声信号やディジタル音響信号のピッチゲインについての情報も入力される。ピッチゲインについての情報は、線形予測分析装置２外にあるピッチゲイン計算部９５０で求められる。
【００２２】
ピッチゲインは、フレームごとの入力信号の周期性の強さのことである。ピッチゲインは、例えば、入力信号やその線形予測残差信号についてのピッチ周期分だけ時間差がある信号間の正規化された相関である。
【００２３】
[ピッチゲイン計算部９５０]
ピッチゲイン計算部９５０は、現フレームの入力信号XO(n) (n=0, 1, …, N-1)および／または現フレームの近傍のフレームの入力信号の全部または一部からピッチゲインGを求める。ピッチゲイン計算部９５０は、例えば、現フレームの入力信号XO(n) (n=0, 1, …, N-1)の全部または一部を含む信号区間のディジタル音声信号やディジタル音響信号のピッチゲインGを求め、ピッチゲインGを特定可能な情報をピッチゲインについての情報として出力する。ピッチゲインを求める方法としては、様々な公知の方法が存在するので、公知の何れの方法を用いてもよい。また、求めたピッチゲインGを符号化してピッチゲイン符号を得る構成とし、ピッチゲイン符号をピッチゲインについての情報として出力してもよい。さらにピッチゲイン符号に対応するピッチゲインの量子化値^Gを得る構成とし、ピッチゲインの量子化値^Gをピッチゲインについての情報として出力してもよい。以下、ピッチゲイン計算部９５０の具体例について説明する。
【００２４】
＜ピッチゲイン計算部９５０の具体例１＞
ピッチゲイン計算部９５０の具体例１は、現フレームの入力信号XO(n) (n=0, 1, …, N-1)が複数個のサブフレームで構成されている場合、かつ、同一のフレームについては線形予測分析装置２よりも先にピッチゲイン計算部９５０が動作される場合、の例である。ピッチゲイン計算部９５０は、まず、2以上の整数であるM個のサブフレームであるXOs1(n) (n=0, 1, …, N/M-1), …, XOsM(n)(n= (M-1)N/M, (M-1)N/M+1, …, N-1)のそれぞれのピッチゲインであるGs1,…, GsMを求める。NはMで割り切れるとする。ピッチゲイン計算部９５０は、現フレームを構成するM個のサブフレームのピッチゲインであるGs1,…, GsMのうちの最大値max(Gs1,…,GsM)を特定可能な情報をピッチゲインについての情報として出力する。
【００２５】
＜ピッチゲイン計算部９５０の具体例２＞
ピッチゲイン計算部９５０の具体例２は、現フレームの入力信号XO(n) (n=0, 1, …, N-1)と１つ後のフレームの一部の入力信号XO(n) (n=N, N+1, …, N+Nn-1)（ただし、Nnは、Nn<Nという関係を満たす所定の正の整数。）とで、先読み部分を含む信号区間が現フレームの信号区間として構成されている場合であり、かつ、同一のフレームについては線形予測分析装置２よりも後にピッチゲイン計算部９５０が動作される場合、の例である。ピッチゲイン計算部９５０は、現フレームの信号区間について、現フレームの入力信号XO(n) (n=0, 1, …, N-1)と１つ後のフレームの一部の入力信号XO(n) (n=N, N+1, …, N+Nn-1)のそれぞれのピッチゲインであるGnow, Gnextを求め、ピッチゲインGnextをピッチゲイン計算部９５０に記憶する。ピッチゲイン計算部９５０は、また、１つ前のフレームの信号区間について求めてピッチゲイン計算部９５０に記憶されていたピッチゲインGnext、すなわち、１つ前のフレームの信号区間のうちの現フレームの一部の入力信号XO(n) (n=0, 1, …, Nn-1)について求めたピッチゲイン、を特定可能な情報をピッチゲインについての情報として出力する。なお、具体例１と同様に、現フレームについては複数のサブフレームごとのピッチゲインを求めてもよい。
【００２６】
＜ピッチゲイン計算部９５０の具体例３＞
ピッチゲイン計算部９５０の具体例３は、現フレームの入力信号XO(n) (n=0, 1, …, N-1)そのものが現フレームの信号区間として構成されている場合であり、かつ、同一のフレームについては線形予測分析装置２よりも後にピッチゲイン計算部９５０が動作される場合、の例である。ピッチゲイン計算部９５０は、現フレームの信号区間である現フレームの入力信号XO(n) (n=0, 1, …, N-1)のピッチゲインGを求め、ピッチゲインGをピッチゲイン計算部９５０に記憶する。ピッチゲイン計算部９５０は、また、１つ前のフレームの信号区間、すなわち、１つ前のフレームの入力信号XO(n) (n=-N, -N+1, …, -1)について求めてピッチゲイン計算部９５０に記憶されていたピッチゲインGを特定可能な情報をピッチゲインについての情報として出力する。
【００２７】
以下、線形予測分析装置２の動作について説明する。図２は、線形予測分析装置２による線形予測分析方法のフローチャートである。
【００２８】
[自己相関計算部２１］
自己相関計算部２１は、入力されたNサンプルのフレーム毎の時間領域のディジタル音声信号やディジタル音響信号である入力信号XO(n)(n=0,1,…,N-1)から自己相関RO(i) (i=0,1,…,Pmax)を計算する（ステップＳ１）。Pmaxは、予測係数計算部２３が求める線形予測係数に変換可能な係数の最大次数であり、N未満の所定の正の整数である。計算された自己相関RO(i) (i=0,1,…,Pmax)は、係数乗算部２２に提供される。
【００２９】
自己相関計算部２１は、入力信号XO(n)を用いて、例えば式(14A)により自己相関RO(i) (i=0,1,…,Pmax)を計算して出力する。すなわち、現在のフレームの入力時系列信号XO(n)とiサンプルだけ過去の入力時系列信号XO(n-i)との自己相関RO(i)を計算する。
【数４】
2018028699.tif 000006
【００３０】
または、自己相関計算部２１は、入力信号XO(n)を用いて、例えば式(14B)により自己相関RO(i) (i=0,1,…,Pmax)を計算する。すなわち、現在のフレームの入力時系列信号XO(n)とiサンプルだけ未来の入力時系列信号XO(n+i)との自己相関RO(i)を計算する。
【数５】
2018028699.tif 000007
【００３１】
または、自己相関計算部２１は、入力信号XO(n)に対応するパワースペクトルを求めてからWiener-Khinchinの定理に従って自己相関RO(i) (i=0,1,…,Pmax)を計算してもよい。また、何れの方法においても、入力信号XO(n) (n=-Np, -Np+1,…, -1, 0,1,…,N-1, N,…, N-1+Nn)というように前後のフレームの入力信号の一部も用いて自己相関RO(i)を計算してもよい。ここで、Np,Nnはそれぞれ、Np<N, Nn<Nという関係を満たす所定の正の整数である。もしくは、MDCT系列をパワースペクトルの近似として代用し、近似されたパワースペクトルから自己相関を求めてもよい。このように自己相関の算出方法は世の中で使われている公知技術の何れかを用いればよい。
【００３２】
［係数決定部２４］
係数決定部２４は、入力されたピッチゲインについての情報を用いて、係数wO(i) (i=0,1,…,Pmax)を決定する（ステップＳ４）。係数wO(i)は、自己相関RO(i)を変形するための係数である。係数wO(i)は、信号処理の分野においては、ラグ窓wO(i)又はラグ窓係数wO(i)とも呼ばれているものである。係数wO(i)は正の値であるので、係数wO(i)が所定の値よりも大きい／小さいことを、係数wO(i)の大きさが所定の値よりも大きい／小さいと表現することがある。また、wO(i)の大きさとは、そのwO(i)の値を意味するものとする。
【００３３】
係数決定部２４に入力されるピッチゲインについての情報は、現フレームの入力信号および／または現フレームの近傍のフレームの入力信号の全部または一部から求まったピッチゲインを特定する情報である。すなわち、係数wO(i)の決定に用いるピッチゲインは、現フレームの入力信号および／または現フレームの近傍のフレームの入力信号の全部または一部から求まったピッチゲインである。
【００３４】
係数決定部２４は、0次からPmax次の全てまたは一部の次数について、ピッチゲインについての情報に対応するピッチゲインの取り得る範囲のうち全てまたは一部で、ピッチゲインについての情報に対応するピッチゲインが大きいほど小さな値を係数wO(0), wO(1), …, wO(Pmax)として決定する。また、係数決定部２４は、ピッチゲインの代わりにピッチゲインと正の相関関係にある値を用いて、ピッチゲインが大きいほど小さな値を係数wO(0), wO(1), …, wO(Pmax)として決定してもよい。
【００３５】
すなわち、係数wO(i) (i=0,1,…,Pmax)は、少なくとも一部の予測次数iに対して、その次数iに対応する係数wO(i)の大きさが、現フレームの入力信号XO(n)の全部または一部を含む信号区間のピッチゲインと正の相関関係にある値の増加とともに単調減少する関係にある場合が含まれているように決定される。
【００３６】
言い換えれば、後述するように、次数iによっては、係数wO(i)の大きさがピッチゲインと正の相関関係にある値の増加とともに単調減少しなくてもよい。
【００３７】
また、ピッチゲインと正の相関関係にある値の取り得る範囲には、係数wO(i)の大きさがピッチゲインと正の相関関係にある値の増加に関わらず一定の範囲があってもよいが、その他の範囲では係数wO(i)の大きさがピッチゲインと正の相関関係にある値の増加とともに単調減少するものとする。
【００３８】
係数決定部２４は、例えば、入力されたピッチゲインについての情報に対応するピッチゲインについての単調非増加関数を用いて、係数wO(i)を決定する。例えば、０より大きい予め定めた値であるαを用いた、以下の式(2)により係数wO(i)を決定する。式(2)において、Gは入力されたピッチゲインについての情報に対応するピッチゲインを意味する。αは、係数wO(i)をラグ窓としてとらえたときのラグ窓の幅、言い換えればラグ窓の強さを調整するための値である。予め定めるαは、例えば、複数のαの候補値について線形予測分析装置２を含む符号化装置とその符号化装置に対応する復号装置とで音声信号や音響信号を符号化復号して、復号音声信号や復号音響信号の主観品質や客観品質が良好である候補値をαとして選択することにより定めればよい。
【数６】
2018028699.tif 000008
【００３９】
または、ピッチゲインGについての予め定めた関数f(G)を用いた、以下の式(2A)により係数wO(i)を決定してもよい。関数f(G)は、f(G)＝αG＋β（αは正の数、βは任意の数）、f(G)＝αG2＋βG+γ（αは正の数、β、γは任意の数）などの、ピッチゲインGと正の相関関係、ピッチゲインGに対して単調非減少の関係となる関数である。
【数７】
2018028699.tif 000009
【００４０】
また、ピッチゲインGを用いて係数wO(i)を決定する式は、上述の(2), (2A)に限らず、ピッチゲインと正の相関関係にある値の増加に対して単調非増加の関係を記述できるものであれば他の式であってもよい。例えば、係数wO(i)を、以下の(3)から(6)の何れかの式により決定してもよい。以下の(3)から(6)の式において、aをピッチゲインに依存して決まる実数とし、mをピッチゲインに依存して決まる自然数とする。例えば、aをピッチゲインと負の相関関係にある値とし、mをピッチゲインと負の相関関係にある値とする。τはサンプリング周期である。
【数８】
2018028699.tif 000010
式(3)はBartlett windowと呼ばれる形式の窓関数であり、式(4)は二項係数により定義されるBinomial windowと呼ばれる形式の窓関数であり、式(5)はTriangular in frequency domain windowと呼ばれる形式の窓関数であり、式(6)はRectangular in frequency domain windowと呼ばれる形式の窓関数である。
【００４１】
なお、0≦i≦Pmaxの各iではなく、少なくとも一部の次数iについてのみ、係数wO(i)がピッチゲインと正の相関関係にある値の増加とともに単調減少してもよい。言い換えれば、次数iによっては、係数wO(i)の大きさがピッチゲインと正の相関関係にある値の増加とともに単調減少しなくてもよい。
【００４２】
例えば、i=0の場合は、上述の式(2)から式(6)の何れかを用いて係数wO(0)の値を決定してもよいし、ITU-T G.718等でも用いられているようなwO(0)＝1.0001,wO(0)＝1.003といった、ピッチゲインと正の相関関係にある値には依存しない、経験的に得られた固定値を用いてもよい。すなわち、1≦i≦Pmaxの各iについては、係数wO(i)はピッチゲインと正の相関関係にある値が大きいほど小さな値を取るが、i=0の係数についてはこの限りではなく固定値を用いてもよい。
【００４３】
[係数乗算部２２]
係数乗算部２２は、係数決定部２４で決定した係数wO(i) (i=0,1,…,Pmax)と、自己相関計算部２１で求めた自己相関RO(i) (i=0,1,…,Pmax)とを同じiごとに乗じることにより、変形自己相関R'O(i) (i=0,1,…,Pmax)を求める（ステップＳ２）。すなわち、係数乗算部２２は、以下の式(7)により自己相関R'O(i)を計算する。計算された自己相関R'O(i)は、予測係数計算部２３に提供される。
【数９】
2018028699.tif 000011
【００４４】
[予測係数計算部２３]
予測係数計算部２３は、係数乗算部２２から出力された変形自己相関R'O(i)を用いて線形予測係数に変換可能な係数を求める（ステップＳ３）。
【００４５】
例えば、予測係数計算部２３は、係数乗算部２２から出力された変形自己相関R'O(i)を用いて、Levinson-Durbin法などにより、1次から予め定めた最大次数であるPmax次までのPARCOR係数KO(1),KO(2),…,KO(Pmax)や線形予測係数aO(1),aO(2),…,aO(Pmax)を計算して出力する。
【００４６】
第一実施形態の線形予測分析装置２によれば、ピッチゲインと正の相関関係にある値に応じて、少なくとも一部の予測次数iに対して、その次数iに対応する係数wO(i)の大きさが、現フレームの入力信号XO(n)の全部または一部を含む信号区間のピッチゲインと正の相関関係にある値の増加とともに単調減少する関係にある場合が含まれている係数wO(i)を自己相関に乗算して変形自己相関を求めて線形予測係数に変換可能な係数を求めることにより、入力信号のピッチゲインが大きいときであってもピッチ成分に起因するスペクトルのピークの発生を抑えた線形予測係数に変換可能な係数を求めることができ、かつ、入力信号のピッチゲインが小さいときであってもスペクトル包絡を表現可能な線形予測係数に変換可能な係数を求めることができ、従来よりも分析精度の高い線形予測を実現することができる。したがって、第一実施形態の線形予測分析装置２を含む符号化装置とその符号化装置に対応する復号装置とで音声信号や音響信号を符号化復号して得られる復号音声信号や復号音響信号の品質は、従来の線形予測分析装置を含む符号化装置とその符号化装置に対応する復号装置とで音声信号や音響信号を符号化復号して得られる復号音声信号や復号音響信号の品質よりも、良い。
【００４７】
［第二実施形態］
第二実施形態は、現在又は過去のフレームにおける入力信号のピッチゲインと正の相関関係にある値と所定の閾値とを比較し、その比較結果に応じて係数wO(i)を決定するものである。第二実施形態は、係数決定部２４における係数wO(i)の決定方法のみが第一実施形態と異なり、他の点については第一実施形態と同様である。以下、第一実施形態と異なる部分を中心に説明し、第一実施形態と同様の部分については重複説明を省略する。
【００４８】
第二実施形態の線形予測分析装置２の機能構成と線形予測分析装置２による線形予測分析方法のフローチャートは、第一実施形態と同じ図１と図２である。第二実施形態の線形予測分析装置２は、係数決定部２４の処理が異なる部分以外は、第一実施形態の線形予測分析装置２と同じである。
【００４９】
第二実施形態の係数決定部２４の処理の流れの例を図３に示す。第二実施形態の係数決定部２４は、図３の各ステップＳ４１Ａ、ステップＳ４２、ステップＳ４３の処理を例えば行う。
【００５０】
係数決定部２４は、入力されたピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値と所定の閾値とを比較する（ステップＳ４１Ａ）。入力されたピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値とは、例えば、入力されたピッチゲインについての情報に対応するピッチゲインそのものである。
【００５１】
係数決定部２４は、ピッチゲインと正の相関関係にある値が所定の閾値以上である場合、すなわち、ピッチゲインが大きいと判断された場合には、予め定めた規則により係数wh(i)を決定し、この決定された係数wh(i) (i=0,1,…,Pmax)をwO(i) (i=0,1,…,Pmax)とする（ステップＳ４２）。すなわち、wO(i)=wh(i)とする。
【００５２】
係数決定部２４は、ピッチゲインと正の相関関係にある値が所定の閾値以上でない場合、すなわち、ピッチゲインが小さいと判断された場合には、予め定めた規則により係数wl(i)を決定し、この決定された係数wl(i) (i=0,1,…,Pmax)をwO(i) (i=0,1,…,Pmax)とする（ステップＳ４３）。すなわち、wO(i)=wl(i)とする。
【００５３】
ここで、wh(i)とwl(i)は、少なくとも一部の各iについてwh(i)<wl(i)という関係を満たすよう決定する。または、wh(i)とwl(i)は、少なくとも一部の各iについてはwh(i)<wl(i)という関係を満たし、それ以外のｉについてはwh(i)≦wl(i)という関係を満たすよう決定する。ここで、少なくとも一部の各iとは、例えば0以外のi（つまり、1≦i≦Pmax）のことである。例えば、wh(i)とwl(i)は、式(2)でピッチゲインGがG1であるときのwO(i)をwh(i)として求め、式(2)でピッチゲインGがG2（ただしG1>G2）であるときのwO(i)をwl(i)として求めるという予め定めた規則により求める。または、例えば、wh(i)とwl(i)は、式(2)でαがα1であるときのwO(i)をwh(i)として求め、式(2)でαがα2（ただしα1>α2）であるときのwO(i)をwl(i)として求めるという予め定めた規則により求める。この場合は、α1とα2は共に式(2)のαと同様に予め定めておく。なお、これらの何れかの規則により予め求めたwh(i) とwl(i)をテーブルに記憶しておき、ピッチゲインと正の相関関係にある値が所定の閾値以上であるか否かによりwh(i)とwl(i)の何れかをテーブルから選択する構成としてもよい。また、wh(i)とwl(i)のそれぞれは、iが大きくなるにつれてwh(i), wl(i)の値が小さくなるように決定される。なお、i=0の係数wh(i), wl(i)については、 wh(0)≦wl(0)の関係を満たしていることは必須ではなく、wh(0)>wl(0)の関係を満たす値を用いてもよい。
【００５４】
第二実施形態によっても、第一実施形態と同様に、入力信号のピッチゲインが大きいときであってもピッチ成分に起因するスペクトルのピークの発生を抑えた線形予測係数に変換可能な係数を求めることができ、かつ、入力信号のピッチゲインが小さいときであってもスペクトル包絡を表現可能な線形予測係数に変換可能な係数を求めることができ、従来よりも分析精度の高い線形予測を実現することができる。
【００５５】
＜第二実施形態の変形例＞
上述の第二実施形態では１個の閾値を用いて係数wO(i)を決定したが、第二実施形態の変形例は２個以上の閾値を用いて係数wO(i)を決定するものである。以下、２個の閾値th1,th2を用いて係数を決定する方法を例に挙げて説明する。閾値th1,th2は、0<th1<th2という関係を満たすとする。
【００５６】
第二実施形態の変形例の線形予測分析装置２の機能構成は、第二実施形態と同じ図１である。第二実施形態の変形例の線形予測分析装置２は、係数決定部２４の処理が異なる部分以外は、第二実施形態の線形予測分析装置２と同じである。
【００５７】
係数決定部２４は、入力されたピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値と、閾値th1,th2とを比較する。入力されたピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値とは、例えば、入力されたピッチゲインについての情報に対応するピッチゲインそのものである。
【００５８】
係数決定部２４は、ピッチゲインと正の相関関係にある値が閾値th2より大きい場合、すなわち、ピッチゲインが大きいと判断された場合には、予め定めた規則により係数wh(i) (i=0,1,…,Pmax)を決定し、この決定された係数wh(i) (i=0,1,…,Pmax)をwO(i) (i=0,1,…,Pmax)とする。すなわち、wO(i)=wh(i)とする。
【００５９】
係数決定部２４は、ピッチゲインと正の相関関係にある値が閾値th1よりも大きく閾値th2以下の場合、すなわち、ピッチゲインが中程度と判断された場合には、予め定めた規則により係数wm(i) (i=0,1,…,Pmax)を決定し、この決定された係数wm(i) (i=0,1,…,Pmax)をwO(i) (i=0,1,…,Pmax)とする。すなわち、wO(i)=wm(i)とする。
【００６０】
係数決定部２４は、ピッチゲインと正の相関関係にある値が閾値th1以下の場合、すなわち、ピッチゲインが小さいと判断された場合には、予め定めた規則により係数wl(i) (i=0,1,…,Pmax)を決定し、この決定された係数wl(i) (i=0,1,…,Pmax)をwO(i) (i=0,1,…,Pmax)とする。すなわち、wO(i)=wl(i)とする。
【００６１】
ここで、wh(i),wm(i),wl(i)は、少なくとも一部の各iについてwh(i)<wm(i)<wl(i)という関係を満たすよう決定するものとする。ここで、少なくとも一部の各iとは、例えば0以外の各i（つまり、1≦i≦Pmax）のことである。または、wh(i),wm(i),wl(i)は、少なくとも一部の各iについてwh(i)<wm(i)≦wl(i)、それ以外のiのうちの少なくとも一部の各iについてwh(i)≦wm(i)<wl(i)、残り少なくとも一部の各iについてwh(i)≦wm(i)≦wl(i)という関係を満たすよう決定する。例えば、wh(i), wm(i), wl(i)は、式(2)でピッチゲインGがG1であるときのwO(i)をwh(i)として求め、式(2)でピッチゲインGがG2（ただしG1>G2）であるときのwO(i)をwm(i)として求め、式(2)でピッチゲインGがG3（ただしG2>G3）であるときのwO(i)をwl(i)として求めるという予め定めた規則により求める。または、例えば、wh(i), wm(i), wl(i)は、式(2)でαがα1であるときのwO(i)をwh(i)として求め、式(2)でαがα2（ただしα1>α2）であるときのwO(i)をwm(i)として求め、式(2)でαがα3（ただしα2>α3）であるときのwO(i)をwl(i)として求めるという予め定めた規則により求める。この場合は、α1, α2, α3は式(2)のαと同様に予め定めておく。なお、これらの何れかの規則により予め求めたwh(i), wm(i), wl(i)をテーブルに記憶しておき、ピッチゲインと正の相関関係にある値と所定の閾値との比較によりwh(i), wm(i), wl(i)の何れかをテーブルから選択する構成としてもよい。
【００６２】
なお、wh(i)とwl(i)を用いて、その間の係数wm(i)を決定しても良い。すなわち、wm(i)=β'×wh(i)+(1-β')×wl(i)によりwm(i)を決定しても良い。ここでβ'は、0≦β'≦1であり、かつ、ピッチゲインGが小さい値をとるときはβ'の値も小さくなり、ピッチゲインGが大きい値をとるときにβ'の値も大きくなる関数β'=c(G)により、ピッチゲインGから求める値である。このようにwm(i)を求めれば、係数決定部２４にはwh(i) (i=0,1,…,Pmax)を記憶したテーブルとwl(i) (i=0,1,…,Pmax)を記憶したテーブルの２つのテーブルだけを記憶しておくことで、ピッチゲインが中程度の場合のうちのピッチゲインが大きいときにはwh(i)に近い係数を得ることができ、逆にピッチゲインが中程度の場合のうちのピッチゲインが小さいときにはwl(i)に近い係数を得ることができる。また、wh(i), wm(i), wl(i)は、iが大きくなるにつれてそれぞれwh(i), wm(i), wl(i)の値が小さくなるよう決定される。なお、i=0の係数wh(0), wm(0), wl(0)については、wh(0)≦wm(0)≦wl(0)の関係を満たしていることは必須ではなく、wh(0)>wm(0)または／およびwm(0)>wl(0)の関係を満たす値を用いてもよい。
【００６３】
第二実施形態の変形例によっても、第二実施形態と同様に、入力信号のピッチゲインが大きいときであってもピッチ成分に起因するスペクトルのピークの発生を抑えた線形予測係数に変換可能な係数を求めることができ、かつ、入力信号のピッチゲインが小さいときであってもスペクトル包絡を表現可能な線形予測係数に変換可能な係数を求めることができ、従来よりも分析精度の高い線形予測を実現することができる。
【００６４】
［第三実施形態］
第三実施形態は、複数個の係数テーブルを用いて係数wO(i)を決定するものである。第三実施形態は、係数決定部２４における係数wO(i)の決定方法のみが第一実施形態と異なり、他の点については第一実施形態と同様である。以下、第一実施形態と異なる部分を中心に説明し、第一実施形態と同様の部分については重複説明を省略する。
【００６５】
第三実施形態の線形予測分析装置２は、係数決定部２４の処理が異なり、図４に例示するように、係数テーブル記憶部２５を更に備えている部分以外は、第一実施形態の線形予測分析装置２と同じである。係数テーブル記憶部２５には、２個以上の係数テーブルが記憶されている。
【００６６】
第三実施形態の係数決定部２４の処理の流れの例を図５に示す。第三実施形態の係数決定部２４は、図５のステップＳ４４、ステップＳ４５の処理を例えば行う。
【００６７】
まず、係数決定部２４は、入力されたピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値を用いて、係数テーブル記憶部２５に記憶された２個以上の係数テーブルから、そのピッチゲインと正の相関関係にある値に応じた１個の係数テーブルtを選択する（ステップＳ４４）。例えば、ピッチゲインについての情報に対応するピッチゲインと正の相関関係にある値は、ピッチゲインについての情報に対応するピッチゲインである。
【００６８】
例えば、係数テーブル記憶部２５に、異なる２個の係数テーブルt0, t1が記憶されており、係数テーブルt0には係数wt0(i) (i=0,1,…,Pmax)が格納されており、係数テーブルt1には係数wt1(i) (i=0,1,…,Pmax)が格納されているとする。２個の係数テーブルt0, t1のそれぞれには、少なくとも一部の各iについてwt0(i)<wt1(i)であり、残りの各iについてwt0(i)≦wt1(i)となるように定められた係数wt0(i) (i=0,1,…,Pmax)と係数wt1(i) (i=0,1,…,Pmax)が格納されている。
【００６９】
このとき、係数決定部２４は、入力されたピッチゲインについての情報により特定されるピッチゲインと正の相関関係にある値が所定の閾値以上であれば係数テーブルt0を係数テーブルtとして選択し、そうでない場合には係数テーブルt1を係数テーブルtとして選択する。すなわち、ピッチゲインと正の相関関係にある値が所定の閾値以上である場合、すなわち、ピッチゲインが大きいと判断された場合には、各ｉについての係数が小さい方の係数テーブルを選択し、ピッチゲインと正の相関関係にある値が所定の閾値より小さい場合、すなわちピッチゲインが小さいと判断された場合には、各ｉについての係数が大きい方の係数テーブルを選択する。
【００７０】
言い換えれば、係数テーブル記憶部２５に記憶されている２個の係数テーブルの中の、ピッチゲインと正の相関関係にある値が第一値である場合に係数決定部２４により選択される係数テーブルを第一係数テーブルとし、係数テーブル記憶部２５に記憶されている２個の係数テーブルの中の、ピッチゲインと正の相関関係にある値が第一値よりも小さい第二値である場合に係数決定部２４により選択される係数テーブルを第二係数テーブルとして、少なくとも一部の各次数iに対して、第二係数テーブルにおける各次数iに対応する係数の大きさは、第一係数テーブルにおける各次数iに対応する係数の大きさよりも大きい。
【００７１】
なお、係数テーブル記憶部２５に記憶されている係数テーブルt0, t1のi=0の係数wt0(0), wt1(0)については、wt0(0)≦wt1(0)の関係を満たしていることは必須ではなく、wt0(0)>wt1(0)の関係にある値であってもよい。
【００７２】
また、例えば、係数テーブル記憶部２５に、異なる３個の係数テーブルt0, t1, t2が記憶されて、係数テーブルt0には係数wt0(i) (i=0,1,…,Pmax)が格納されており、係数テーブルt1には係数wt1(i) (i=0,1,…,Pmax) 、係数テーブルt2には係数wt2(i) (i=0,1,…,Pmax)が格納されているとする。３個の係数テーブルt0, t1, t2のそれぞれには、少なくとも一部の各iについてwt0(i)<wt1(i)≦wt2(i)であり、それ以外のiのうちの少なくとも一部の各iについてwt0(i)≦wt1(i)<wt2(i)であり、残りの各iについてwt0(i)≦wt1(i)≦wt2(i)となるように定められた係数wt0(i) (i=0,1,…,Pmax)と係数wt1(i) (i=0,1,…,Pmax) と係数wt2(i) (i=0,1,…,Pmax)が格納されている。
【００７３】
ここで、0<th1<th2という関係を満たす２個の閾値th1,th2が定められているとする。このとき、係数決定部２４は、
(1) ピッチゲインと正の相関関係にある値>th2の場合、すなわち、ピッチゲインが大きいと判断された場合には、係数テーブルt0を係数テーブルtとして選択し、
(2) th2≧ ピッチゲインと正の相関関係にある値>th1の場合、すなわち、ピッチゲインが中程度であると判断された場合には、係数テーブルt1を係数テーブルtとして選択し、
(3) th1≧ ピッチゲインと正の相関関係にある値の場合、すなわち、ピッチゲインが小さいと判断された場合には、係数テーブルt2を係数テーブルtとして選択する。
【００７４】
なお、係数テーブル記憶部２５に記憶されている係数テーブルt0, t1, t2のi=0の係数wt0(0), wt1(0), wt2(0)については、wt0(0)≦wt1(0)≦wt2(0)の関係を満たしていることは必須ではなく、wt0(0)>wt1(0)または／およびwt1(0)>wt2(0)の関係にある値であってもよい。
【００７５】
そして、係数決定部２４は、その選択された係数テーブルtに格納された各次数iの係数wt(i)を係数wO(i)とする（ステップＳ４５）。すなわち、wO(i)=wt(i)とする。言い換えれば、係数決定部２４は、選択された係数テーブルtから各次数iに対応する係数wt(i)を取得し、取得された各次数iに対応する係数wt(i)をwO(i)とする。
【００７６】
第三実施形態では、第一実施形態及び第二実施形態とは異なり、ピッチゲインと正の相関関係にある値の式に基づいて係数wO(i)を計算する必要がないため、より少ない演算処理量でwO(i)を決定することができる。
【００７７】
＜第三実施形態の具体例＞
以下、第三実施形態の具体例について説明する。線形予測分析装置２には、ハイパスフィルタを通り、12.8 kHzにサンプリング変換され、プリエンファシス処理をされた１フレームあたりNサンプルのディジタル音響信号である入力信号XO(n) (n=0,1,…,N-1)と、ピッチゲインについての情報として現フレームの一部の入力信号XO(n) (n=0, 1, …, Nn)（ただし、Nnは、Nn<Nという関係を満たす所定の正の整数。）についてピッチゲイン計算部９５０で求めたピッチゲインGとが入力される。現フレームの一部の入力信号XO(n) (n=0, 1, …, Nn)についてのピッチゲインGは、ピッチゲイン計算部９５０において当該入力信号の１つ前のフレームの信号区間として現フレームの一部の入力信号XO(n) (n=0, 1, …, Nn)を含めておき、１つ前のフレームの信号区間に対するピッチゲイン計算部９５０の処理においてXO(n) (n=0, 1, …, Nn)に対して計算し記憶したピッチゲインである。
【００７８】
自己相関計算部２１は、入力信号XO(n)から自己相関RO(i) (i=0,1,…,Pmax)を下記の式(8)で求める。
【数１０】
2018028699.tif 000012
【００７９】
係数決定部２４に、ピッチゲインについての情報であるピッチゲインGが入力される。
【００８０】
係数テーブル記憶部２５には、係数テーブルt0と、係数テーブルt1と、係数テーブルt2とが記憶されているものとする。
【００８１】
係数テーブルt0、式(13)の従来法のf0=60Hzの係数テーブルであり、各次数の係数wtO(i)が次のように定められている。
【００８２】
wt0(i)=[1.0001, 0.999566371, 0.998266613, 0.996104103, 0.993084457, 0.989215493, 0.984507263, 0.978971839, 0.972623467, 0.96547842, 0.957554817, 0.948872864, 0.939454317, 0.929322779, 0.918503404, 0.907022834, 0.894909143]
係数テーブルt1は、式(13)の従来法のf0=40Hzのテーブルであり、各次数の係数wt1(i)が次のように定められている。
【００８３】
wt1(i)=[1.0001, 0.999807253, 0.99922923, 0.99826661, 0.99692050, 0.99519245, 0.99308446, 0.99059895, 0.98773878, 0.98450724, 0.98090803, 0.97694527, 0.97262346, 0.96794752, 0.96292276, 0.95755484, 0.95184981]
係数テーブルt2は、式(13)の従来法のf0=20Hzのテーブルであり、各次数の係数wt2(i)が次のように定められている。
【００８４】
wt2(i)=[1.0001, 0.99995181, 0.99980725, 0.99956637, 0.99922923, 0.99879594, 0.99826661, 0.99764141, 0.99692050, 0.99610410, 0.99519245, 0.99418581, 0.99308446, 0.99188872, 0.99059895, 0.98921550, 0.98773878]
ここで、上述のwtO(i), wt1(i), wt2(i)のリストは、Pmax=16として、i=0,1,2,…,16の順に左からiに対応する係数の大きさを並べたものである。すなわち上述の例では、例えばwt0(0)=1.0001であり、wt0(3)=0.996104103である。
【００８５】
図６に係数テーブルt0,t1,t2の係数wt0(i), wt1(i), wt2(i)の大きさをグラフで表す。図６のグラフの点線は係数テーブルt0の係数wt0(i)の大きさを表し、図６のグラフの一点鎖線は係数テーブルt1の係数wt1(i)の大きさを表し、図６のグラフの実線は係数テーブルt2の係数wt2(i)の大きさを表す。図６のグラフの横軸は次数iを意味し、図６のグラフの縦軸は係数の大きさを表す。このグラフからも分かるように、各係数テーブル内では、iの値が大きくなるにしたがって、係数の大きさが単調減少する関係にある。また、同じiの値に対応する異なる係数テーブルの係数の大きさを比較すると、0を除くi≧１のiに対して、言い換えれば少なくとも一部のiについて、wt0(i)<wt1(i)<wt2(i)の関係を満たしている。係数テーブル記憶部２５に記憶される複数の係数テーブルは、このような関係を持つものであれば、上述の例に限らない。
【００８６】
また、非特許文献１や非特許文献２に記載されているように、i=0の係数だけ特別扱いをして、wt0(0)=wt1(0)=wt2(0)=1.0001やwt0(0)=wt1(0)=wt2(0)=1.003という経験的な値を用いてもよい。なお、i=0についてはwt0(i)<wt1(i)<wt2(i)の関係を満たしている必要はなく、また、wt0(0),wt1(0),wt2(0)が必ずしも同じ値でなくともよい。例えば、wt0(0)=1.0001, wt1(0)=1.0, wt2(0)=1.0のように、i=0に関してのみwt0(0), wt1(0), wt2(0)のうちの2つ以上の値の大小関係がwt0(i)<wt1(i)<wt2(i)の関係を満たさなくてもよい。
【００８７】
上述の係数テーブルt0は式(13)においてf0=60Hz, fs=12.8kHzとした場合、係数テーブルt1は式(13)においてf0=40Hz, fs=12.8kHzとした場合、係数テーブルt2は式(13)においてf0=20Hzとした場合の係数値に相当するが、これらはそれぞれ、式(2A)においてf(G)=60, fs=12.8kHzとした場合の係数値、f(G)=40, fs=12.8kHzとした場合、f(G)=20, fs=12.8kHzとした場合、に相当し、式(2A)における関数f(G)はピッチゲインGと正の相関関係にある関数である。つまり、３つの係数テーブルの係数値を予め定める際に、予め定めた３つのピッチゲインを用いて式(2A)により係数値を求めることに代えて、予め定めた３つのf0を用いて式(13)により係数値を求めてもよい。
【００８８】
係数決定部２４は、入力されたピッチゲインGを所定の閾値th1=0.3及び閾値th2=0.6と比較し、G≦0.3の場合は係数テーブルt2を、0.3<G≦0.6の場合は係数テーブルt1を、0.6<Gの場合は係数テーブルt0を選択する。
【００８９】
そして、係数決定部２４は、その選択された係数テーブルtの各係数wt(i)を係数wO(i)とする。すなわち、wO(i)=wt(i)とする。言い換えれば、係数決定部２４は、選択された係数テーブルtから各次数iに対応する係数wt(i)を取得し、取得された各次数iに対応する係数wt(i)をwO(i)とする。
【００９０】
＜第三実施形態の変形例＞
第三実施形態では複数個の係数テーブルのうち何れか１つのテーブルに記憶された係数を係数wO(i)として決定したが、第三実施形態の変形例はこれに加えて複数個の係数テーブルに記憶された係数に基づく演算処理により係数wO(i)を決定する場合を含む。
【００９１】
第三実施形態の変形例の線形予測分析装置２の機能構成は、第三実施形態と同じ図４である。第三実施形態の変形例の線形予測分析装置２は、係数決定部２４の処理が異なり、係数テーブル記憶部２５に含まれる係数テーブルが異なる部分以外は、第三実施形態の線形予測分析装置２と同じである。
【００９２】
係数テーブル記憶部２５には、係数テーブルt0とt2のみが記憶されており、係数テーブルt0には係数wt0(i) (i=0,1,…,Pmax)が格納されており、係数テーブルt2には係数wt2(i) (i=0,1,…,Pmax)が格納されている。２個の係数テーブルt0, t2のそれぞれには、少なくとも一部の各iについてwt0(i)<wt2(i)であり、残りの各iについてwt0(i)≦wt2(i)となるように定められた係数wt0(i) (i=0,1,…,Pmax)と係数wt2(i) (i=0,1,…,Pmax)が格納されている。
【００９３】
ここで、0<th1<th2という関係を満たす２個の閾値th1,th2が定められているとする。このとき、係数決定部２４は、
(1) ピッチゲインと正の相関関係にある値>th2の場合、すなわち、ピッチゲインが大きいと判断された場合には、係数テーブルt0の各係数wt0(i)を係数wO(i)として選択し、
(2) th2≧ピッチゲインと正の相関関係にある値>th1の場合、すなわち、ピッチゲインが中程度であると判断された場合には、係数テーブルt0の各係数wt0(i)と係数テーブルt2の各係数wt2(i)とを用いて、wO(i)=β'×wt0(i)+(1-β')×wt2(i)により係数wO(i)を決定し、
(3) th1≧ピッチゲインと正の相関関係にある値の場合、すなわち、ピッチゲインが小さいと判断された場合には、係数テーブルt2の各係数wt2(i)を係数wO(i)として選択する。
【００９４】
ここでβ'は、0≦β'≦1であり、ピッチゲインGが小さい値をとるときはβ'の値も小さくなり、ピッチゲインGが大きい値をとるときにβ'の値も大きくなる関数β'=c(G)により、ピッチゲインGから求める値である。この構成とすれば、ピッチゲインが中程度の場合のうちのピッチゲインGが小さい時にはwt2(i)に近い値を係数wO(i)とすることができ、逆にピッチゲインが中程度の場合のうちのピッチゲインGが大きい時にはwt0(i)に近い値を係数wO(i)とすることができるので、２つのテーブルだけで、３個以上の係数wO(i)を得ることができる。
【００９５】
なお、係数テーブル記憶部２５に記憶されている係数テーブルt0, t2のi=0の係数wt0(0), wt2(0)については、wt0(0)≦wt2(0)の関係を満たしていることは必須ではなく、wt0(0)>wt2(0)の関係にある値であってもよい。
【００９６】
［第一実施形態から第三実施形態に共通の変形例］
図７及び図８に示すように、上述の全ての実施形態及び変形例において、係数乗算部２２を含まず、予測係数計算部２３において係数wO(i)と自己相関RO(i)とを用いて線形予測分析を行ってもよい。図７と図８は、それぞれ図１と図４に対応する線形予測分析装置２の構成例である。この場合は、予測係数計算部２３は、図９のステップＳ５において、係数wO(i)と自己相関RO(i)とが乗算されたものである変形自己相関R'O(i)ではなく、係数wO(i)と自己相関RO(i)とを直接用いて線形予測分析を行う（ステップＳ５）。
【００９７】
[第四実施形態］
第四実施形態は、入力信号XO(n)に対して従来の線形予測分析装置を用いて線形予測分析を行い、その線形予測分析の結果を用いてピッチゲイン計算部でピッチゲインを得て、得られたピッチゲインに基づく係数wO(i)を用いて本発明の線形予測分析装置により線形予測係数に変換可能な係数を求めるものである。
【００９８】
第四実施形態の線形予測分析装置３は、図１０に示すように、第一線形予測分析部３１、線形予測残差計算部３２、ピッチゲイン計算部３６、第二線形予測分析部３４を例えば備えている。
【００９９】
[第一線形予測分析部３１］
第一線形予測分析部３１は、従来の線形予測分析装置１と同じ動作をする。すなわち、第一線形予測分析部３１は、入力信号XO(n)から自己相関RO(i) (i=0,1,…,Pmax)を求め、自己相関RO(i) (i=0,1,…,Pmax)と予め定めた係数wO(i) (i=0,1,…,Pmax)とを同じiごとに乗じることにより変形自己相関R' O(i) (i=0,1,…,Pmax)を求め、変形自己相関R' O(i) (i=0,1,…,Pmax)から1次から予め定めた最大次数であるPmax次までの線形予測係数に変換可能な係数を求める。
【０１００】
[線形予測残差計算部３２］
線形予測残差計算部３２は、入力信号XO(n)に対して、1次からPmax次までの線形予測係数に変換可能な係数に基づく線形予測や線形予測と等価なまたは類似したフィルタリング処理を行って線形予測残差信号XR(n)を求める。フィルタリング処理は重み付け処理とも言えるので、線形予測残差信号XR(n)は重み付け入力信号であるともいえる。
【０１０１】
[ピッチゲイン計算部３６]
ピッチゲイン計算部３６は、線形予測残差信号XR(n)のピッチゲインGを求め、ピッチゲインについての情報を出力する。ピッチゲインを求める方法としては、様々な公知の方法が存在するので、公知の何れの方法を用いてもよい。ピッチゲイン計算部３６は、例えば、現フレームの線形予測残差信号XR (n) (n=0, 1, …, N-1)を構成する複数個のサブフレームのそれぞれについてピッチゲインを求める。すなわち、2以上の整数であるM個のサブフレームであるXRs1(n) (n=0, 1, …, N/M-1), …, XRsM(n) (n=(M-1)N/M, (M-1)N/M+1, …, N-1)のそれぞれのピッチゲインであるGs1, …, GsMを求める。NはMで割り切れるとする。ピッチゲイン計算部３６は、次に、現フレームを構成するM個のサブフレームのピッチゲインであるGs1, …, GsMのうちの最大値max(Gs1, …, GsM)を特定可能な情報をピッチゲインについての情報として出力する。
【０１０２】
[第二線形予測分析部３４］
第二線形予測分析部３４は、本発明の第一実施形態から第三実施形態及びこれらの変形例の線形予測分析装置２の何れかと同じ動作をする。すなわち、第二線形予測分析部３４は、入力信号XO(n)から自己相関RO(i) (i=0,1,…,Pmax)を求め、ピッチゲイン計算部３６が出力したピッチゲインについての情報に基づいて係数wO(i) (i=0,1,…,Pmax)を決定し、自己相関RO(i) (i=0,1,…,Pmax)と決定した係数wO(i) (i=0,1,…,Pmax)とを用いて変形自己相関R' O(i) (i=0,1,…,Pmax)から1次から予め定めた最大次数であるPmax次までの線形予測係数に変換可能な係数を求める。
【０１０３】
＜ピッチゲインと正の相関関係にある値について＞
第一実施形態においてピッチゲイン計算部９５０の具体例２として説明した通り、ピッチゲインと正の相関関係にある値として、前のフレームの信号処理においてLook-aheadとも呼ばれる先読みして利用するサンプル部分のうち現フレームのサンプルに対応する部分のピッチゲインを用いてもよい。
【０１０４】
また、ピッチゲインと正の相関関係にある値として、ピッチゲインの推定値を用いてもよい。例えば、過去の複数フレームのピッチゲインから予測される現在のフレームについてのピッチゲインの推定値や、過去の複数フレームについてのピッチゲインの平均値や最小値や最大値や重み付線形和を、ピッチゲインの推定値として用いてもよい。また、複数サブフレームについてのピッチゲインの平均値や最小値や最大値や重み付線形和を、ピッチゲインの推定値として用いてもよい。
【０１０５】
また、ピッチゲインと正の相関関係にある値として、ピッチゲインの量子化値を用いてもよい。すなわち、量子化前のピッチゲインを用いてもよいし、量子化後のピッチゲインを用いてもよい。
【０１０６】
なお、上記の各実施形態及び各変形例のピッチゲインと正の相関関係にある値と閾値との比較においては、ピッチゲインと正の相関関係にある値が閾値と同じ値である場合には、閾値を境として隣接する二つの場合の何れか一方に場合分けされるように設定すればよい。すなわち、ある閾値以上の場合としているところを当該閾値より大きい場合とするとともに、当該閾値より小さい場合としているところを当該閾値以下の場合としてもよい。また、ある閾値より大きい場合としているところを当該閾値以上の場合とするとともに、当該閾値以下の場合としているところを当該閾値より小さい場合としてもよい。
【０１０７】
上記装置及び方法において説明した処理は、記載の順にしたがって時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。
【０１０８】
また、線形予測分析方法における各ステップをコンピュータによって実現する場合、線形予測分析方法が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、その各ステップがコンピュータ上で実現される。
【０１０９】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。
【０１１０】
また、各処理手段は、コンピュータ上で所定のプログラムを実行させることにより構成することにしてもよいし、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【０１１１】
その他、この発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。
【図１】
2018028699.tif 000013
【図２】
2018028699.tif 000014
【図３】
2018028699.tif 000015
【図４】
2018028699.tif 000016
【図５】
2018028699.tif 000017
【図６】
2018028699.tif 000018
【図７】
2018028699.tif 000019
【図８】
2018028699.tif 000020
【図９】
2018028699.tif 000021
【図１０】
2018028699.tif 000022
【図１１】
2018028699.tif 000023