(57)【特許請求の範囲】
【請求項１】
周波数毎の係数値で構成されてステレオ形式の第１音響信号および第２音響信号における特定位置の定位成分を強調または抑圧する第１数値列を生成する第１生成手段と、
前記第１音響信号と前記第２音響信号との周波数毎の位相差に応じた係数値で構成される第２数値列を、前記第１音響信号の振幅と前記第２音響信号の振幅と両者の和成分の振幅とから生成する第２生成手段と、
周波数毎の係数値で構成されて前記定位成分を強調または抑圧する第１処理係数列を前記第１数値列と前記第２数値列とから生成する合成手段と、
前記第１音響信号および前記第２音響信号の各々の各周波数成分に前記第１処理係数列の各係数値を作用させる信号処理手段と
を具備する音響処理装置。
【請求項２】
前記第２生成手段は、
前記第１音響信号および前記第２音響信号の各々の振幅を加算した振幅和と前記和成分の振幅とを周波数毎に算定する振幅算定手段と、
前記振幅和と前記和成分の振幅との相違に応じて前記第２数値列の各係数値を設定する設定手段と
を含む請求項１の音響処理装置。
【請求項３】
前記第１音響信号と前記第２音響信号との加算で和成分を生成する和成分生成手段と、
前記定位成分を抑圧した差成分を前記第１音響信号と前記第２音響信号との間の減算で生成する差成分生成手段とを具備し、
前記第１生成手段は、前記和成分および前記差成分から前記第１数値列を生成し、
前記第２生成手段は、前記和成分生成手段が生成した前記和成分の振幅を前記第２数値列の生成に適用する
請求項１または請求項２の音響処理装置。
【請求項４】
前記第１生成手段は、前記和成分および前記差成分の一方のスペクトルから他方のスペクトルを減算した結果に応じて前記第１数値列を生成する
請求項３の音響処理装置。
【請求項５】
前記定位成分の強調および抑圧の一方に対応する前記第１処理係数列の各係数値を所定値から減算することで前記定位成分の強調および抑圧の他方に対応する第２処理係数列を生成する第３生成手段
を具備する請求項１から請求項４の何れかの音響処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ステレオ形式の音響信号のうち目的の位置（方向）に音像が定位する成分（例えば、右チャネルおよび左チャネルの各々の信号に含まれる成分である。以下では「定位成分」という）を強調または抑圧する技術に関する。
【背景技術】
【０００２】
ＣＤ等の記録媒体に収録された音響信号から所定の定位成分を抑圧する技術が従来から提案されている。例えば特許文献１には、ステレオ形式の音響信号の一方から他方を減算（逆相加算）することで、音像が中央方向に定位する成分（典型的には歌唱音）を抑圧する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−０７９４１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、特許文献１の技術では、所定の成分の抑圧後の音響信号がモノラル形式になるという問題がある。同様の問題は、定位成分を強調する場合にも発生し得る。以上の事情を考慮して、本発明は、定位成分を強調または抑圧したステレオ形式の音響信号を生成することを目的とする。
【課題を解決するための手段】
【０００５】
以上の課題を解決するために、本発明の音響処理装置は、周波数毎の係数値で構成されてステレオ形式の第１音響信号および第２音響信号における特定位置の定位成分を強調または抑圧する第１数値列を生成する第１生成手段と、第１音響信号と第２音響信号との周波数毎の位相差に応じた係数値で構成される第２数値列を、第１音響信号の振幅と第２音響信号の振幅と両者の和成分の振幅とから生成する第２生成手段と、周波数毎の係数値で構成されて定位成分を強調または抑圧する第１処理係数列を第１数値列と第２数値列とから生成する合成手段と、第１音響信号および第２音響信号の各々の各周波数成分に第１処理係数列の各係数値を作用させる信号処理手段とを具備する。以上の構成においては、第１音響信号および第２音響信号の各々の各周波数成分に第１処理係数列を作用させるから、定位成分が強調または抑圧されたステレオ形式の音響信号を生成することが可能である。また、第１音響信号と第２音響信号との位相差が第１処理係数列に反映されるから、第１音響信号と第２音響信号との間に顕著な位相差がある場合でも、定位成分の強調または抑圧の効果を充分に維持できるという利点がある。
【０００６】
本発明の好適な態様において、第２生成手段は、第１音響信号および第２音響信号の各々の振幅を加算した振幅和と和成分の振幅とを周波数毎に算定する振幅算定手段と、振幅和と和成分の振幅との相違に応じて第２数値列の各係数値を設定する設定手段とを含んで構成される。以上の態様においては、第１音響信号と第２音響信号との位相差に応じて第１音響信号および第２音響信号の振幅和と和成分の振幅との相違が変化するという傾向を利用して、第１音響信号と第２音響信号との位相差を高精度に反映した第２数値列を生成することが可能である。
【０００７】
本発明の好適な態様に係る音響処理装置は、第１音響信号と第２音響信号との加算で和成分を生成する和成分生成手段と、定位成分を抑圧した差成分を第１音響信号と第２音響信号との間の減算で生成する差成分生成手段とを具備し、第１生成手段は、和成分および差成分から第１数値列を生成し、第２生成手段は、和成分生成手段が生成した和成分の振幅を第２数値列の生成に適用する。以上の態様においては、第１数値列の生成に利用される和成分が第２数値列の生成にも流用されるから、第２数値列の生成の負荷が軽減されるという利点がある。
【０００８】
本発明の好適な態様において、第１生成手段は、和成分および差成分の一方のスペクトルから他方のスペクトルを減算した結果に応じて第１数値列を生成する。以上の態様においては、和成分および差成分の一方のスペクトルから他方のスペクトルを減算することで定位成分を高精度に強調または抑圧したスペクトルが生成されるから、定位成分を有効に強調または抑圧し得る第１処理係数列を生成できるという利点がある。
【０００９】
差成分生成手段を具備する態様の好適例に係る音響処理装置は、特定位置を示す定位変数を可変に設定する変数設定手段を具備し、差成分生成手段は、変数設定手段が設定した定位変数に応じた比率で第１音響信号および第２音響信号の一方から他方を減算する。以上の態様においては、差成分生成手段による減算時の第１音響信号および第２音響信号の比率が定位変数に応じて可変に設定される。したがって、強調または抑圧の対象となる定位成分の位置（方向）を可変に制御できる（中央方向に限定されない）という利点がある。
【００１０】
本発明の好適な態様に係る音響処理装置は、定位成分の強調および抑圧の一方に対応する第１処理係数列の各係数値を所定値から減算することで定位成分の強調および抑圧の他方に対応する第２処理係数列を生成する第３生成手段を具備する。以上の態様においては、定位成分の強調用および抑圧用の双方の処理係数列を生成することが可能である。また、第１処理係数列の各係数値を所定値から減算することで第２処理係数列が生成されるから、第１処理係数列および第２処理係数列の各々の個別に生成する場合と比較して、第１処理係数列および第２処理係数列の生成の負荷が軽減されるという利点がある。
【００１１】
以上の各態様に係る音響処理装置は、音響信号の処理に専用されるＤＳＰ（Digital Signal Processor）などのハードウェア（電子回路）によって実現されるほか、ＣＰＵ（Central Processing Unit）などの汎用の演算処理装置とプログラム（ソフトウェア）との協働によっても実現される。本発明のプログラムは、周波数毎の係数値で構成されてステレオ形式の第１音響信号および第２音響信号における特定位置の定位成分を強調または抑圧する第１数値列を生成する第１生成処理と、第１音響信号と第２音響信号との周波数毎の位相差に応じた係数値で構成される第２数値列を、第１音響信号の振幅と第２音響信号の振幅と両者の和成分の振幅とから生成する第２生成処理と、周波数毎の係数値で構成されて定位成分を強調または抑圧する第１処理係数列を第１数値列と第２数値列とから生成する合成処理と、第１音響信号および第２音響信号の各々の各周波数成分に第１処理係数列の各係数値を作用させる信号処理とをコンピュータに実行させる。以上のプログラムによれば、本発明の音響処理装置と同様の作用および効果が実現される。本発明のプログラムは、コンピュータが読取可能な記録媒体に格納された形態で利用者に提供されてコンピュータにインストールされるほか、通信網を介した配信の形態でサーバ装置から提供されてコンピュータにインストールされる。
【図面の簡単な説明】
【００１２】
【図１】第１実施形態に係る音響処理装置のブロック図である。
【図２】係数設定部のブロック図である。
【図３】各周波数成分の関係を示す模式図である。
【図４】第２生成部のブロック図である。
【図５】各周波数成分の関係を示す模式図である。
【図６】変形例における第３生成部のブロック図である。
【図７】定位変数と係数値ＣA[k]との関係を示す模式図である。
【図８】定位変数と係数値ＣD[k]との関係を示す模式図である。
【図９】定位変数と位相差との関係を示す模式図である。
【発明を実施するための形態】
【００１３】
＜Ａ：第１実施形態＞
図１は、本発明の第１実施形態に係る音響処理装置１００のブロック図である。音響処理装置１００には信号供給装置１２と放音装置１４と入力装置１６とが接続される。信号供給装置１２は、音響（音声や楽音）の波形を表す時間領域の音響信号ＳIN（ＳIN_L，ＳIN_R）を音響処理装置１００に供給する。左チャネルの音響信号ＳIN_Lおよび右チャネルの音響信号ＳIN_Rは、音響を発生する複数の音源の音像が相異なる位置に定位する（すなわち、音響の振幅や位相が各音源の位置に応じて相違する）ように収音または加工されたステレオ形式の信号である。周囲の音響を収音して音響信号ＳINを生成する収音機器（ステレオマイク）や、可搬型または内蔵型の記録媒体から音響信号ＳINを取得して音響処理装置１００に出力する再生装置や、通信網から音響信号ＳINを受信して音響処理装置１００に出力する通信装置が信号供給装置１２として採用され得る。
【００１４】
音響処理装置１００は、信号供給装置１２が供給する音響信号ＳINから音響信号ＳOUT（ＳOUT_L，ＳOUT_R）を生成する。左チャネルの音響信号ＳOUT_Lおよび右チャネルの音響信号ＳOUT_Rは、音響信号ＳINが表す音響のうち目的の位置に音像が定位する定位成分を強調または抑圧したステレオ形式の信号である。放音装置１４（例えばステレオスピーカやステレオヘッドホン）は、音響処理装置１００が生成した音響信号ＳOUT（ＳOUT_L，ＳOUT_R）に応じた音波を放射する。
【００１５】
入力装置１６は、音響処理装置１００に対する指示を利用者が入力するための機器（例えばマウスやキーボード）である。利用者は、入力装置１６を適宜に操作することで、定位成分の位置（方向）と定位成分の強調／抑圧の何れかの処理とを音響処理装置１００に対して任意に指示することが可能である。
【００１６】
図１に示すように、音響処理装置１００は、演算処理装置２２と記憶装置２４とを具備するコンピュータシステムで実現される。記憶装置２４は、演算処理装置２２が実行するプログラムＰGや演算処理装置２２が使用するデータを記憶する。半導体記録媒体や磁気記録媒体などの公知の記録媒体や複数種の記録媒体の組合せが記憶装置２４として任意に採用される。音響信号ＳIN（ＳIN_L，ＳIN_R）を記憶装置２４に記憶した構成（したがって信号供給装置１２は省略され得る）も好適である。
【００１７】
演算処理装置２２は、記憶装置２４に格納されたプログラムＰGを実行することで、音響信号ＳINから音響信号ＳOUTを生成するための複数の機能（変数設定部３２，周波数分析部３４，係数設定部３６，信号処理部３８，波形合成部４０）を実現する。なお、演算処理装置２２の各機能を複数の集積回路に分散した構成や、専用の電子回路（ＤＳＰ）が各機能を実現する構成も採用され得る。
【００１８】
変数設定部３２は、定位変数αを可変に設定する。定位変数αは、定位成分の音像の位置（方向）を指定する変数である。本形態の変数設定部３２は、入力装置１６に対する利用者からの指示に応じて定位変数α（０≦α≦１）を可変に設定する。例えば、定位変数αが０．５（中央値）である場合には中央（正面）方向が指示される。また、定位変数αが１に近いほど右寄りの方向が指示され、定位変数αが０に近いほど左寄りの方向が指示される。利用者は、放音装置１４からの再生音を聴取しながら随時に入力装置１６の操作で定位変数αを変更することが可能である。
【００１９】
周波数分析部３４は、音響信号ＳIN_Lの周波数スペクトル（複素スペクトル）ＬAと音響信号ＳIN_Rの周波数スペクトル（複素スペクトル）ＲAとを時間軸上の単位区間（フレーム）毎に順次に生成する。周波数スペクトルＬAは、Ｋ個の周波数（周波数帯域）ｆ1〜ｆKの各々に対応する周波数成分ＬAk(ejω)の系列（ＬA1(ejω)〜ＬAK(ejω)）である（ｋ＝１〜Ｋ）。記号ｊは虚数単位を意味する。同様に、周波数スペクトルＲAは、周波数ｆ1〜ｆKの各々に対応する周波数成分ＲAk(ejω)の系列（ＲA1(ejω)〜ＲAK(ejω)）である。周波数スペクトルＬAおよび周波数スペクトルＲAの生成には、短時間フーリエ変換などの公知の周波数分析が任意に採用され得る。なお、通過帯域が相異なるＫ個の帯域通過フィルタで構成されるフィルタバンクも周波数分析部３４として採用され得る。
【００２０】
図１の係数設定部３６は、音響信号ＳIN（ＳIN_L，ＳIN_R）の定位成分を強調（enhance）するための処理係数列Ｇeと音響信号ＳINの定位成分を抑圧（suppress）するための処理係数列Ｇsとを生成する。処理係数列Ｇeおよび処理係数列Ｇsは、周波数スペクトルＬAと周波数スペクトルＲAとを利用して単位区間毎に順次に生成される。処理係数列Ｇeは、周波数ｆ1〜ｆKの各々に対応する係数値Ｇe[k]の系列（Ｇe[1]〜Ｇe[K]）であり、処理係数列Ｇsは、周波数ｆ1〜ｆKの各々に対応する係数値Ｇs[k]の系列（Ｇs[1]〜Ｇs[K]）である。
【００２１】
係数値Ｇe[k]および係数値Ｇs[k]は、音響信号ＳIN_Lの周波数成分ＬAk(ejω)や音響信号ＳIN_Rの周波数成分ＲAk(ejω)に対するゲイン（スペクトルゲイン）に相当し、音響信号ＳIN（ＳIN_L，ＳIN_R）の特性に応じて０以上かつ１以下の範囲内で可変に設定される（０≦Ｇe[k]≦１，０≦Ｇs[k]≦１）。具体的には、定位成分の強調用の処理係数列Ｇeの係数値Ｇe[1]〜Ｇe[K]は、定位成分のパワー（振幅）が大きい周波数ｆkの係数値Ｇe[k]ほど１に近い数値に設定される。他方、定位成分の抑圧用の処理係数列Ｇsの係数値Ｇs[1]〜Ｇs[K]は、定位成分のパワーが大きい周波数ｆkの係数値Ｇs[k]ほど０に近い数値に設定される。
【００２２】
図１の信号処理部３８は、係数設定部３６が生成した処理係数列Ｇeまたは処理係数列Ｇsを周波数スペクトルＬAおよび周波数スペクトルＲAの各々に個別に作用させる（典型的には乗算する）ことで周波数スペクトルＬBと周波数スペクトルＲBとを単位区間毎に順次に生成する。周波数スペクトルＬBは、周波数ｆ1〜ｆKの各々に対応する周波数成分ＬBk(ejω)の系列（ＬB1(ejω)〜ＬBK(ejω)）であり、周波数スペクトルＲBは、周波数ｆ1〜ｆKの各々に対応する周波数成分ＲBk(ejω)の系列（ＲB1(ejω)〜ＲBK(ejω)）である。各単位区間の周波数スペクトルＬAおよび周波数スペクトルＲAの処理には、その単位区間について係数設定部３６が生成した処理係数列Ｇeまたは処理係数列Ｇsが適用される。
【００２３】
入力装置１６に対して利用者から定位成分の強調が指示された場合、信号処理部３８は、以下の数式(1a)および数式(1b)に示すように、周波数スペクトルＬAおよび周波数スペクトルＲAに対する処理係数列Ｇeの乗算で周波数スペクトルＬBおよび周波数スペクトルＲBを生成する。すなわち、周波数スペクトルＬBの各周波数ｆkの周波数成分ＬBk(ejω)は、その周波数ｆkの周波数成分ＬAk(ejω)と係数値Ｇe[k]との乗算値に設定され（数式(1a)）、周波数スペクトルＲBの各周波数ｆkの周波数成分ＲBk(ejω)は、その周波数ｆkの周波数成分ＲAk(ejω)と係数値Ｇe[k]との乗算値に設定される（数式(1b)）。したがって、音響信号ＳIN（ＳIN_L，ＳIN_R）の定位成分を強調した周波数スペクトルＬBおよび周波数スペクトルＲBが生成される。
【数１】
0005454157.tif 000002
【００２４】
他方、入力装置１６に対して利用者から定位成分の抑圧が指示された場合、信号処理部３８は、以下の数式(2a)および数式(2b)に示すように、周波数スペクトルＬBおよび周波数スペクトルＲBの生成に定位成分の抑圧用の処理係数列Ｇsを適用する。すなわち、周波数スペクトルＬAの各周波数成分ＬAk(ejω)と処理係数列Ｇsの各係数値Ｇs[k]との乗算で周波数スペクトルＬBの各周波数成分ＬBk(ejω)を算定し（数式(2a)）、周波数スペクトルＲAの各周波数成分ＲAk(ejω)と各係数値Ｇs[k]との乗算で周波数スペクトルＲBの各周波数成分ＲBk(ejω)を算定する（数式(2b)）。したがって、音響信号ＳIN（ＳIN_L，ＳIN_R）の定位成分を抑圧した周波数スペクトルＬBおよび周波数スペクトルＲBが生成される。
【数２】
0005454157.tif 000003
【００２５】
図１の波形合成部４０は、信号処理部３８による処理後の周波数スペクトルＬBおよび周波数スペクトルＲBからステレオ形式の音響信号ＳOUT_Lおよび音響信号ＳOUT_Rを生成する。具体的には、波形合成部４０は、単位区間毎の周波数スペクトルＬBを逆フーリエ変換で時間領域の信号に変換するとともに前後の単位区間について相互に連結することで音響信号ＳOUT_Lを生成する。同様に、波形合成部４０は、各周波数スペクトルＲBから音響信号ＳOUT_Rを生成する。波形合成部４０が生成した音響信号ＳOUT（ＳOUT_L，ＳOUT_R）が放音装置１４に供給されて音波として再生される。
【００２６】
以上に説明したように、第１実施形態では、音響信号ＳIN_L（周波数スペクトルＬA）および音響信号ＳIN_R（周波数スペクトルＲA）の各々に対して個別に処理係数列Ｇeまたは処理係数列Ｇsが乗算されるから、定位成分を強調または抑圧した音響信号ＳOUT（ＳOUT_L，ＳOUT_R）をステレオ形式のまま生成することが可能である。
【００２７】
次に、係数設定部３６の詳細を説明する。図２は、係数設定部３６のブロック図である。図２に示すように、係数設定部３６は、和成分生成部５２と差成分生成部５４と係数列生成部６０とを含んで構成される。和成分生成部５２は、音響信号ＳIN_Lの周波数成分ＬA1(ejω)〜ＬAK(ejω)と音響信号ＳIN_Rの周波数成分ＲA1(ejω)〜ＲAK(ejω)との加算で単位区間毎に順次に和成分（複素スペクトル）Ｍを生成する。和成分Ｍは、周波数ｆ1〜ｆKの各々に対応する周波数成分Ｍk(ejω)の系列（Ｍ1(ejω)〜ＭK(ejω)）である。数式(3)に示すように、和成分Ｍの周波数ｆkの周波数成分Ｍk(ejω)は、その周波数ｆkの周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との加算（複素数）に相当する。したがって、和成分Ｍは、全部の音源からの音響を混合したモノラル形式の信号に相当する。なお、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との加重和や平均を和成分Ｍとして算定する構成も採用され得る。
【数３】
0005454157.tif 000004
【００２８】
図２の差成分生成部５４は、音響信号ＳIN_Lの周波数成分ＬA1(ejω)〜ＬAK(ejω)と音響信号ＳIN_Rの周波数成分ＲA1(ejω)〜ＲAK(ejω)との間の減算で単位区間毎に順次に差成分（複素スペクトル）Ｓを生成する。差成分Ｓは、周波数ｆ1〜ｆKの各々に対応する周波数成分Ｓk(ejω)の系列（Ｓ1(ejω)〜ＳK(ejω)）である。差成分生成部５４は、変数設定部３２が設定した定位変数αを利用した数式(4)の演算（加重減算）で周波数成分Ｓ1(ejω)〜ＳK(ejω)を算定する。
【数４】
0005454157.tif 000005
【００２９】
数式(4)から理解されるように、定位変数αに応じた可変の比率（重み値）にて周波数成分ＬAk(ejω)から周波数成分ＲAk(ejω)を減算（逆相加算）することで差成分Ｓの各周波数成分Ｓk(ejω)が生成される。したがって、差成分Ｓは、音響信号ＳIN（ＳIN_L，ＳIN_R）のうち定位変数αに応じた位置（方向）の定位成分を他の成分に対して相対的に抑圧した信号（すなわち、定位成分以外の成分を相対的に強調した信号）となる。例えば、定位変数αが０．５（中央値）である場合には、中央方向の定位成分（すなわち、振幅および位相が略同等の成分）を抑圧した差成分Ｓが生成される。また、定位変数αが０．５を上回るほど、中央方向に対して右寄りの定位成分が差成分Ｓでは抑圧され、定位変数αが０．５を下回るほど、中央方向に対して左寄りの定位成分が差成分Ｓでは抑圧される。
【００３０】
数式(4)の記号ｍａｘ（α，１−α）は、定位変数αまたは変数（１−α）のうちの最大値を意味する。定位変数αは１以下の数値に設定されるから、数式(4)の分子の演算のみでは周波数成分Ｓk(ejω)のパワー（振幅）が不足する可能性がある。数式(4)のように最大値ｍａｘ（α，１−α）で除算するのは、周波数成分Ｓk(ejω)のパワーを周波数成分ＬAk(ejω)や周波数成分ＲAk(ejω)と同等に維持するためである。
【００３１】
図２の係数列生成部６０は、和成分生成部５２が生成した和成分Ｍと差成分生成部５４が生成した差成分Ｓとを利用して定位成分の強調用の処理係数列Ｇe（Ｇe[1]〜Ｇe[K]）と定位成分の抑圧用の処理係数列Ｇs（Ｇs[1]〜Ｇs[K]）とを生成する。図２に示すように、係数列生成部６０は、第１生成部６２と第２生成部６４と第３生成部６６と合成部６８とを含んで構成される。
【００３２】
第１生成部６２は、定位成分の強調用の数値列ＣAを和成分Ｍと差成分Ｓとから生成する。数値列ＣAは、周波数ｆ1〜ｆKの各々に対応する係数値（ゲイン）ＣA[k]の系列（ＣA[1]〜ＣA[K]）である。具体的には、第１生成部６２は、以下の数式(5)の演算で数値列ＣAの各係数値ＣA[k]を算定する。
【数５】
0005454157.tif 000006
【００３３】
数式(5)の記号Ｐ[k]は、定位成分を強調したパワースペクトルＰのうち周波数ｆkでのパワーを意味する。パワーＰ[k]は、例えば以下の数式(6a)および数式(6b)で算定される。
【数６】
0005454157.tif 000007
【００３４】
数式(6a)から理解されるように、周波数成分Ｍk(ejω)のパワー|Ｍk(ejω)|２が周波数成分Ｓk(ejω)のパワー|Ｓk(ejω)|２を上回る周波数ｆkのパワーＰ[k]は、和成分Ｍのパワー|Ｍk(ejω)|２から差成分Ｓのパワー|Ｓk(ejω)|２を減算した数値に設定される。すなわち、パワースペクトルＰは、和成分Ｍのパワースペクトルと差成分Ｓのパワースペクトルとの間の減算（スペクトル減算）で生成される。他方、パワー|Ｍk(ejω)|２がパワー|Ｓk(ejω)|２以下となる周波数ｆkのパワーＰ[k]は、和成分Ｍのパワー|Ｍk(ejω)|２と所定の係数（フロアリング係数）βとの乗算値に設定される。以上の説明から理解されるように、数式(6a)および数式(6b)の演算は、和成分Ｍを信号成分（目的音）と仮定するとともに差成分Ｓを雑音成分と仮定した場合に雑音成分を抑圧するためのスペクトル減算（ＳＳ：Spectral Subtraction）に相当する。なお、数式(5)にて振幅Ｐ[k]1/2を周波数成分Ｍk(ejω)の振幅|Ｍk(ejω)|で除算するのは、係数値ＣA[k]を１以下の数値（０≦ＣA[k]≦１）に正規化するためである。
【００３５】
周波数成分Ｓk(ejω)は定位成分を抑圧した成分であるから、数式(6a)および数式(6b)で算定されるパワーＰ[1]〜Ｐ[K]の系列は、音響信号ＳIN（ＳIN_L，ＳIN_R）の定位成分を強調した成分のパワースペクトルＰとなる。したがって、数式(5)で算定される数値列ＣAの係数値ＣA[1]〜ＣA[K]は、定位変数αに応じた位置の定位成分のパワー（振幅）が大きい周波数ｆkの係数値ＣA[k]ほど１に近い数値となり、定位成分のパワーが小さい周波数ｆkの係数値ＣA[k]ほど０に近い数値となる。
【００３６】
以上の説明から理解されるように、第１生成部６２が生成した数値列ＣAを処理係数列Ｇeとして（係数値ＣA[k]を数式(1a)および数式(1b)の係数値Ｇe[k]として）、音響信号ＳIN_Lの周波数成分ＬAk(ejω)および音響信号ＳIN_Rの周波数成分ＲAk(ejω)の各々に乗算する構成（以下「対比例」という）でも、定位成分が強調されたステレオ形式の音響信号ＳOUT（ＳOUT_L，ＳOUT_R）を生成することが可能である。
【００３７】
しかし、音響信号ＳIN_L（周波数成分ＬAk(ejω)）と音響信号ＳIN_R（周波数成分ＲAk(ejω)）とで位相が相違する場合、数式(4)の演算では定位成分が充分に抑圧されず、差成分Ｓに有意なパワーで定位成分が残留する。以上のように差成分Ｓに残留した定位成分は数式(6a)の演算で和成分Ｍから減算される（すなわち、パワースペクトルＰにおける定位成分のパワーが低減される）。したがって、音響信号ＳIN_Lの周波数成分ＬAk(ejω)と音響信号ＳIN_Rの周波数成分ＲAk(ejω)との位相差が大きいほど、第１生成部６２が生成する数値列ＣAの係数値ＣA[k]は小さい数値になるという傾向がある。つまり、数値列ＣAを処理係数列Ｇeとして利用する対比例の構成では、定位成分の強調の効果が不足する結果となる。相互に離れた位置に設置された２個の収音機器を利用した収音で音響空間内での反射後の音響（残響音）の音響信号ＳIN（ＳIN_L，ＳIN_R）を生成した場合や、音響信号ＳINに残響効果を付与した場合には、音響信号ＳIN_Lと音響信号ＳIN_Rとの位相差が顕著となるから、対比例における定位成分の強調の不足は特に顕在化する。
【００３８】
そこで、第２生成部６４は、音響信号ＳIN_L（周波数成分ＬAk(ejω)）と音響信号ＳIN_R（周波数成分ＲAk(ejω)）との位相差に応じた数値列ＣBを生成する。数値列ＣBは、周波数ｆ1〜ｆKの各々に対応する係数値ＣB[k]の系列（ＣB[1]〜ＣB[K]）である。周波数ｆkの周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど、数値列ＣBのうち周波数ｆkの係数値ＣB[k]は大きい数値に設定される。
【００３９】
図３は、音響信号ＳIN_Lの周波数成分ＬAk(ejω)と音響信号ＳIN_Rの周波数成分ＲAk(ejω)と和成分Ｍの周波数成分Ｍk(ejω)とを複素平面内のベクトル（ただし、各成分の符号はＬA，ＲA，Ｍと略記されている）として表記した模式図である。周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との角度θが両者の位相差に相当する。なお、図３では、周波数成分ＬAk(ejω)の振幅|ＬAk(ejω)|と周波数成分ＲAk(ejω)の振幅|ＲAk(ejω)|とが相等しい場合を便宜的に例示した。
【００４０】
図３の部分(A)のように周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)とが同位相である場合（θ＝０）、和成分Ｍの振幅|Ｍk(ejω)|は、周波数成分ＬAk(ejω)の振幅|ＬAk(ejω)|と周波数成分ＲAk(ejω)の振幅|ＲAk(ejω)|との加算値（以下「振幅和」という）Ａkに一致する（Ａk＝|ＬAk(ejω)|＋|ＲAk(ejω)|）。他方、図３の部分(B)および部分(C)に示すように、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差（角度θ）が大きいほど、和成分Ｍの振幅|Ｍk(ejω)|は振幅和Ａkと比較して小さい数値となる。すなわち、振幅和Ａkと和成分Ｍの振幅|Ｍk(ejω)|との相違が大きいほど（振幅|Ｍk(ejω)|が振幅和Ａkと比較して小さいほど）、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいという傾向がある。以上の傾向を利用して、第１実施形態の第２生成部６４は、周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)の振幅和Ａkと和成分Ｍの振幅|Ｍk(ejω)|とを比較した結果に応じて数値列ＣBの各係数値ＣB[k]（ＣB[1]〜ＣB[K]）を決定する。
【００４１】
図４に示すように、第２生成部６４は、振幅算定部６４２と設定部６４４とを含んで構成される。振幅算定部６４２は、周波数成分ＬA1(ejω)〜ＬAK(ejω)の各々の振幅|ＬAk(ejω)｜と、周波数成分ＲA1(ejω)〜ＲAK(ejω)の各々の振幅|ＲAk(ejω)｜と、和成分生成部５２が生成した和成分Ｍの周波数成分Ｍ1(ejω)〜ＭK(ejω)の各々の振幅|Ｍk(ejω)|とを単位区間毎に順次に算定する。設定部６４４は、振幅算定部６４２が算定した各振幅を利用して数値列ＣBの各係数値ＣB[k]を１以上の範囲内で可変に設定する（ＣB[k]≧１）。
【００４２】
具体的には、設定部６４４は、周波数ｆ1〜ｆKの各々について、振幅|ＬAk(ejω)｜および振幅|ＲAk(ejω)｜の振幅和Ａkと和成分Ｍの振幅|Ｍk(ejω)|との差分ΔＡk（ΔＡk＝|ＬAk(ejω)｜＋|ＲAk(ejω)｜−|Ｍk(ejω)|）を算定し、差分ΔＡkが大きいほど係数値ＣB[k]を大きい数値に設定する。差分ΔＡkが０である場合、設定部６４４は係数値ＣB[k]を例えば１に設定する。設定部６４４が差分ΔＡkから係数値ＣB[k]を決定する具体的な方法は任意であるが、例えば、差分ΔＡkの各数値と係数値ＣB[k]の各数値とを対応させたテーブルから係数値ＣB[k]を検索する方法や、差分ΔＡkを変数として係数値ＣB[k]を定義する演算式（例えば差分ΔＡkに対する増加関数）の演算で係数値ＣB[k]を算定する方法が好適である。以上の説明のように、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きい（差分ΔＡkが大きい）ほど、各数値列ＣBの係数値ＣB[k]は大きい数値に設定される。
【００４３】
図２の合成部６８は、第１生成部６２が生成した数値列ＣAと第２生成部６４が生成した数値列ＣBとから処理係数列Ｇe（Ｇe[1]〜Ｇe[K]）を生成する。処理係数列Ｇeの周波数ｆkの係数値Ｇe[k]は、例えば、数値列ＣAの周波数ｆkの係数値ＣA[k]と数値列ＣBの周波数ｆkの係数値ＣB[k]との乗算値（Ｇe[k]＝ＣA[k]×ＣB[k]）に設定される。ただし、係数値ＣA[k]と係数値ＣB[k]との乗算値が１を上回る場合、合成部６８は係数値Ｇe[k]を１に設定する。
【００４４】
以上の説明から理解されるように、定位変数αに応じた位置の定位成分において周波数ｆkの成分が優勢となるほど、または、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど、処理係数列Ｇeの係数値Ｇe[k]は大きい数値に設定される。したがって、周波数スペクトルＬAおよび周波数スペクトルＲAの各々に信号処理部３８が処理係数列Ｇeを乗算することで、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きい場合でも、定位成分を充分に強調したステレオ形式の音響信号ＳOUT（ＳOUT_L，ＳOUT_R）を生成することが可能である。
【００４５】
図２の第３生成部６６は、合成部６８が生成した処理係数列Ｇeを利用して定位成分の抑圧用の処理係数列Ｇs（係数値Ｇs[1]〜Ｇs[K]の系列）を生成する。具体的には、第３生成部６６は、数式(7)に示すように、処理係数列Ｇeの各係数値Ｇe[k]を所定値（本形態では１）から減算することで処理係数列Ｇsの各係数値Ｇs[k]を算定する。
【数７】
0005454157.tif 000008
【００４６】
数式(7)から理解されるように、第３生成部６６が生成する処理係数列Ｇsの係数値Ｇs[1]〜Ｇs[K]は、定位変数αに応じた位置の定位成分のパワー（振幅）が大きい周波数ｆkの係数値Ｇs[k]ほど０に近い数値となり、定位成分のパワーが小さい周波数ｆkの係数値Ｇs[k]ほど１に近い数値となる。したがって、周波数スペクトルＬAおよび周波数スペクトルＲAの各々に信号処理部３８が処理係数列Ｇsを乗算することで、前述の通り、定位成分を抑圧したステレオ形式の音響信号ＳOUT（ＳOUT_L，ＳOUT_R）が生成される。
【００４７】
ところで、前述の対比例では、第１生成部６２が生成した数値列ＣAを数式(7)の処理係数列Ｇeとして抑圧用の処理係数列Ｇsを生成する構成が採用され得る（Ｇs[k]＝１−ＣA[k]）。しかし、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど数値列ＣAの係数値ＣA[k]は小さい数値になるから、対比例の構成では、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど処理係数列Ｇsの係数値Ｇs[k]は大きい数値となる。したがって、定位成分の抑圧が不足するという問題が発生する。他方、第１実施形態では、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど係数値Ｇe[k]は大きい数値となるから、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど処理係数列Ｇsの係数値Ｇs[k]は小さい数値となる。したがって、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きい場合でも、定位成分を充分に抑圧したステレオ形式の音響信号ＳOUT（ＳOUT_L，ＳOUT_R）を生成することが可能である。
【００４８】
以上に説明したように、第１実施形態においては、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が処理係数列Ｇeおよび処理係数列Ｇsに反映されるから、周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)の間に位相差がある場合でも、定位成分の強調または抑圧の効果を充分に維持できるという利点がある。
【００４９】
ところで、特許第３６７０５６２号公報には、ステレオ信号のチャネル間の振幅比に応じた減衰係数ｇi(k)と位相差に応じた減衰係数ｇp(k)との何れかを選択して各チャネルの音響信号に乗算する構成が開示されている。しかし、以上の技術では、振幅比に応じた減衰係数ｇi(k)と位相差に応じた減衰係数ｇp(k)とが択一的に適用されるから、振幅および位相の双方がチャネル間で相違する場合には、定位成分の適切な強調または抑圧が困難であるという問題がある。他方、第１実施形態においては、音響信号ＳIN_L（周波数スペクトルＬA）および音響信号ＳIN_R（周波数スペクトルＲA）から複素スペクトルとして算定された和成分Ｍおよび差成分Ｓが処理係数列Ｇeおよび処理係数列Ｇsの生成に利用されるから、音響信号ＳIN_Lと音響信号ＳIN_Rとの間の振幅差および位相差の双方を反映した処理係数列Ｇeおよび処理係数列Ｇsが生成される。したがって、音響信号ＳIN_Lと音響信号ＳIN_Rとの間で振幅および位相の一方のみが相違する場合に加えて、振幅および位相の双方が相違する場合にも、定位成分を有効に強調または抑圧できるという格別の効果が実現される。
【００５０】
＜Ｂ：第２実施形態＞
本発明の第２実施形態を説明する。第１実施形態では、周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)の振幅和Ａkと和成分Ｍの周波数成分Ｍk(ejω)の振幅|Ｍk(ejω)|との比較の結果に応じて第２生成部６４が数値列ＣBの各係数値ＣB[k]を設定した。第２実施形態では、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との角度θ（すなわち両者の位相差）を算定して係数値ＣB[k]を設定する。なお、以下の各例示において作用や機能が第１実施形態と同等である要素については、第１実施形態と同様の符号を付して各々の詳細な説明を適宜に省略する。
【００５１】
図５は、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)と和成分Ｍの周波数成分Ｍk(ejω)との関係を示す模式図である。図５の三角形ＴMLは、周波数成分ＬAk(ejω)および周波数成分Ｍk(ejω)を２辺とする三角形であり、三角形ＴMRは、周波数成分ＲAk(ejω)および周波数成分Ｍk(ejω)を２辺とする三角形である。三角形ＴMLおよび三角形ＴMRの面積Ｓは相等しく、以下の数式(8)で表現される。数式(8)の角度θ1は、周波数成分ＬAk(ejω)と周波数成分Ｍk(ejω)との角度であり、角度θ2は、周波数成分ＲAk(ejω)と周波数成分Ｍk(ejω)との角度である。
【数８】
0005454157.tif 000009
【００５２】
数式(8)から以下の数式(9a)および数式(9b)が導出される。
【数９】
0005454157.tif 000010
【００５３】
他方、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)と周波数成分Ｍk(ejω)とを３辺とする三角形の面積Ｓは、三角形ＴMLや三角形ＴMRの面積と相等しく、以下の数式(10)で表現される（ヘロンの公式）。
【数１０】
0005454157.tif 000011
【００５４】
第２生成部６４は、数式(10)の演算で算定される面積Ｓを数式(9a)および数式(9b)に適用することで角度θ1および角度θ2を算定する。そして、第２生成部６４は、角度θ1および角度θ2を利用した以下の数式(11)の演算で数値列ＣBの各係数値ＣB[k]を算定する（θ1＋θ2＜π/２）。
【数１１】
0005454157.tif 000012
【００５５】
数式(11)から理解されるように、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差（角度(θ1＋θ2)が大きいほど数値列ＣBの係数値ＣB[k]は大きい数値に設定される。なお、数値列ＣAおよび数値列ＣBを利用して処理係数列Ｇeや処理係数列Ｇsを生成する方法は第１実施形態と同様である。
【００５６】
第２実施形態でも第１実施形態と同様の作用および効果が実現される。第１実施形態や第２実施形態の例示から理解されるように、第２生成部６４は、周波数成分ＬAk(ejω)の振幅|ＬAk(ejω)|と周波数成分ＲAk(ejω)の振幅|ＲAk(ejω)|と周波数成分Ｍk(ejω)の振幅|Ｍk(ejω)|とから数値列ＣBの各係数値ＣB[k]を生成する要素として包括され、各振幅から係数値ＣB[k]を算定する具体的な方法は本発明において任意である。
【００５７】
＜Ｃ：変形例＞
以上の各形態は多様に変形され得る。具体的な変形の態様を以下に例示する。以下の例示から任意に選択された２以上の態様は適宜に併合され得る。
【００５８】
（１）変形例１  
数値列ＣBの各係数値ＣB[k]を算定する方法は以上の例示に限定されない。例えば、図３に例示した関係から理解されるように、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差が大きいほど、和成分Ｍの振幅|Ｍk(ejω)|に対する振幅和Ａkの相対比Ｑk（Ｑk＝Ａk／|Ｍk(ejω)|）は大きい数値となるから、相対比Ｑkを数値列ＣBの各係数値ＣB[k]として算定する構成が採用され得る。また、差分ΔＡk（ΔＡk＝Ａk−|Ｍk(ejω)|）や相対比Ｑkを変数とする所定の関数の演算（写像）で各係数値ＣB[k]を算定する構成も好適である。
【００５９】
（２）変形例２  
以上の各形態では、数値列ＣAと数値列ＣBとを合成した強調用の処理係数列Ｇeから第３生成部６６が抑圧用の処理係数列Ｇsを生成したが、第３生成部６６の代わりに図６の第３生成部６７が処理係数列Ｇsを生成する構成も採用され得る。
【００６０】
図６の第３生成部６７は、第１演算部６７２と第２演算部６７４とを含んで構成される。第１演算部６７２は、第１生成部６２が生成した数値列ＣAから定位成分の抑圧用の数値列ＣD（係数値ＣD[1]〜ＣD[K]の系列）を算定する。例えば、第１演算部６７２は、前述の数式(7)と同様に、数値列ＣAの係数値ＣA[k]を所定位置（例えば１）から減算することで数値列ＣDの各係数値ＣD[k]（ＣD[k]＝１−ＣA[k]）を算定する。したがって、抑圧用の処理係数列Ｇsと同様に、定位変数αに応じた位置の定位成分のパワーが大きい周波数ｆkの係数値ＣD[k]ほど０に近い数値となる。
【００６１】
図６の第２演算部６７４は、周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)の位相差に応じた数値列ＣB'を、第１演算部６７２が生成した数値列ＣDに作用させることで抑圧用の処理係数列Ｇsを生成する。数値列ＣB'は、以上の各形態の数値列ＣBとは逆に、周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)の位相差が大きいほど小さい数値となるように設定された係数値ＣB'[1]〜ＣB'[K]の系列である。例えば、振幅和Ａkに対する和成分Ｍの振幅|Ｍk(ejω)|の相対比（|Ｍk(ejω)|／Ａk）が係数値ＣB'[k]として算定され得る。第２演算部６７４は、例えば、数値列ＣDの係数値ＣD[k]と数値列ＣB'の係数値ＣB'[k]との乗算値を処理係数列Ｇsの係数値Ｇs[k]として算定する．図６の構成でも第１実施形態と同様の効果が実現される。
【００６２】
（３）変形例３  
以上の各形態では、数値列ＣBの全部の係数値ＣB[1]〜ＣB[K]を係数列ＣAの全部の係数値ＣA[1]〜ＣA[K]に作用させて処理係数列Ｇeを生成したが、所定の周波数帯域内の成分に限定して定位成分の強調／抑圧を実行する構成も好適である。例えば、音響信号ＳINのうち人間の音声のパワーが集中する周波数帯域（例えば１００ｋＨｚ〜８ｋＨｚ）についてのみ以上の各形態の処理が実行される（他の帯域については処理せずに再生する）。
【００６３】
また、以下に例示するように、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差に起因した定位成分の強調または抑圧の不足が問題となる係数値ＣA[k]に限定して係数値ＣB[k]を作用させる構成も採用され得る。
【００６４】
図７は、第１生成部６２が生成する数値列ＣAのうち中央方向（α＝０．５）に音像が定位する成分の周波数ｆkに対応する係数値ＣA[k]（縦軸）と定位変数α（横軸）との関係を示す模式図である。図７に示すように、利用者から指示される定位変数αが０．５に近いほど数値列ＣAの各係数列ＣA[k]は１に近い数値に設定される。以上の関係から理解されるように、係数値ＣA[k]が所定の閾値ＣA_THを上回る周波数ｆkの成分は、強調の対象となる定位成分である可能性が高い。
【００６５】
そこで、数値列ＣAの係数値ＣA[1]〜ＣA[K]のうち閾値ＣA_THを上回る係数値ＣA[k]に対して限定的に合成部６８が係数値ＣB[k]を乗算して係数値Ｇe[k]を算定する構成が採用され得る。閾値ＣA_THを下回る係数値ＣA[k]は、係数値ＣB[k]が乗算されることなく係数値Ｇe[k]として採択される。以上の構成によれば、Ｋ個の係数値ＣA[1]〜ＣA[K]にＫ個の係数値ＣB[1]〜ＣB[K]を作用させる第１実施形態の構成と比較して、合成部６８による処理の負荷（乗算の回数）が軽減されるという利点がある。
【００６６】
また、係数値ＣD[k]（ＣD[k]＝１−ＣA[k]）に係数値ＣB'[k]を作用させて抑圧用の処理係数列Ｇsを生成する図６の構成でも同様に、特定の係数値ＣD[k]に対して限定的に係数値ＣB'[k]を作用させる構成が採用され得る。
【００６７】
図８は、図６の構成のもとで第３生成部６７の第１演算部６７２が生成する係数列ＣDのうち中央方向（α＝０．５）に音像が定位する成分の周波数ｆkに対応する係数値ＣD[k]（縦軸）と定位変数α（横軸）との関係を示す模式図である。図８に示すように、利用者から指示される定位変数αが０．５に近いほど数値列ＣDの各係数値ＣD[k]は０に近い数値に設定される。以上の関係から理解されるように、係数値ＣD[k]が所定の閾値ＣD_THを下回る周波数ｆkの成分は、抑圧の対象となる定位成分である可能性が高い。
【００６８】
そこで、数値列ＣDの係数値ＣD[1]〜ＣD[K]のうち閾値ＣD_THを下回る係数値ＣD[k]に対して限定的に図６の第２演算部６７４が係数値ＣB'[k]を乗算して係数値Ｇs[k]を算定する構成が採用され得る。閾値ＣD_THを上回る係数値ＣD[k]は、係数値ＣB'[k]が乗算されることなく係数値Ｇs[k]として採択される。以上の構成によれば、Ｋ個の係数値ＣD[1]〜ＣD[K]にＫ個の係数値ＣB'[1]〜ＣB'[K]を作用させる構成と比較して、第２演算部６７４による処理の負荷（乗算の回数）が軽減されるという利点がある。
【００６９】
（４）変形例４  
定位成分の周波数ｆkに対応する周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差φkは、定位変数αに応じた定位成分の方向に応じて変化する。例えば、中央方向（α＝０．５）の定位成分の周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差φkは０に近い数値となり、定位成分の方向が中央方向から離れる（定位変数αが中央値（０．５）から離れる）ほど位相差φkは増加する。したがって、図９に示すように、定位変数αに対応する位相差φaを含む所定幅の範囲Ａを想定すると、周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差φkが範囲Ａ内にある周波数ｆkの成分は定位成分に該当する可能性が高い。
【００７０】
そこで、係数列ＣAの係数値ＣA[1]〜ＣA[K]のうち周波数成分ＬAk(ejω)と周波数成分ＲAk(ejω)との位相差φkが範囲Ａ内にある各周波数ｆk（すなわち定位成分の周波数ｆk）の各係数値ＣA[k]に限定して、合成部６８が数値列ＣBの係数値ＣB[k]を作用させる構成も採用され得る。位相差φkが範囲Ａ外にある周波数ｆkの各係数値ＣA[k]は、係数値ＣB[k]が乗算されることなく係数値Ｇe[k]として採択される。図６の構成でも同様であり、位相差φkが範囲Ａ内にある周波数kの各係数値ＣD[k]に対して限定的に係数値ＣB'[k]を乗算する構成が採用され得る。
【００７１】
（５）変形例５  
数値列ＣAと数値列ＣBとから処理係数列Ｇeや処理係数列Ｇsを生成する方法は適宜に変更され得る。例えば、数値列ＣAの各係数値ＣA[k]と数値列ＣBの各係数値ＣB[k]との加算値を処理係数列Ｇeの係数値Ｇe[k]として算定する構成が採用され得る。図６の構成でも同様に、各係数値ＣD[k]と各係数値ＣB'[k]との加算で係数値Ｇs[k]を算定する構成が採用され得る。
【００７２】
（６）変形例６  
以上の各形態では、数式(4)から理解されるように、音響信号ＳIN_L（周波数成分ＬAk(ejω)）と音響信号ＳIN_R（周波数成分ＲAk(ejω)）との振幅の比率が定位変数αに応じて線形に変化するように差成分Ｓを算定したが、以下の各態様にて例示するように、和成分Ｍおよび差成分Ｓの算定の方法や定位変数αとの関係は適宜に変更される。
【００７３】
＜第１態様＞
音響信号ＳIN_Lのパワー|ＬAk(ejω)|２および音響信号ＳIN_Rのパワー|ＲAk(ejω)|２から和成分Ｍおよび差成分Ｓを算定する構成が採用され得る。具体的には、和成分生成部５２は、以下の数式(12a)の演算で和成分Ｍ（周波数成分Ｍ1(ejω)〜ＭK(ejω)で構成される複素スペクトル）を生成し、差成分生成部５４は、以下の数式(12b)の演算で差成分Ｓ（周波数成分Ｓ1(ejω)〜ＳK(ejω)で構成される複素スペクトル）を生成する。したがって、差成分Ｓにおいては、音響信号ＳIN_Lおよび音響信号ＳIN_Rの各々のパワーの比率が定位変数αに応じて線形に変化する。なお、数式(12a)および数式(12b)における記号ｅjLは周波数スペクトルＬAの位相スペクトルを意味し、記号ｅjRは周波数スペクトルＲAの位相スペクトルを意味する。
【数１２】
0005454157.tif 000013
【００７４】
以上の方法で和成分Ｍおよび差成分Ｓが生成されると、係数列生成部６０の第１生成部６２は、定位成分を強調したパワースペクトルＰ（パワーＰ[k]）を以下の数式(13a)および数式(13b)の演算で生成し、パワースペクトルＰを利用した数式(14)の演算で処理係数列Ｇeの各係数値Ｇe[k]を算定する。第３生成部６６が処理係数列Ｇeから処理係数列Ｇsを生成する方法は第１実施形態（数式(7)）と同様である。以上の構成でも第１実施形態と同様の効果が実現される。
【数１３】
0005454157.tif 000014
【００７５】
＜第２態様＞
定位成分の位置（方向）を定位変数αの関数ｆ(α)に応じて制御する構成が採用され得る。具体的には、和成分生成部５２は、以下の数式(15a)の演算で和成分Ｍ（周波数成分Ｍ1(ejω)〜ＭK(ejω)）を生成し、差成分生成部５４は、以下の数式(15b)の演算で差成分Ｓ（周波数成分Ｓ1(ejω)〜ＳK(ejω)）を生成する。したがって、差成分Ｓにおける音響信号ＳIN_Lおよび音響信号ＳIN_Rの各々の振幅の比率が定位変数αの関数ｆ(α)に応じて変化する。係数列生成部６０は、以上の各形態と同様の方法で、数式(15a)の和成分Ｍと数式(15b)の差成分Ｓとから処理係数列Ｇeおよび処理係数列Ｇsを算定する。
【数１４】
0005454157.tif 000015
【００７６】
なお、数式(6a)では和成分Ｍのパワー|Ｍk(ejω)|２と差成分Ｓのパワー|Ｓk(ejω)|２との差分をパワーＰ[k]として算定したが、以下の数式(6c)に示すように、和成分Ｍの振幅|Ｍk(ejω)|と差成分Ｓの振幅|Ｓk(ejω)|との差分を振幅Ｐ[k]として算定する構成も採用され得る。数式(6c)で振幅Ｐ[k]を算定する構成では、第１生成部６２が数式(5a)の演算で処理係数列Ｇeの係数値Ｇe[k]を算定する。
【数１５】
0005454157.tif 000016
【００７７】
（７）変形例７  
以上の各形態では、パワースペクトルの減算（数式(6a)，数式(13a)）や振幅スペクトルの減算（数式(6c)）で処理係数列Ｇeを算定したが、処理係数列Ｇeを算定する方法は任意である。例えば、和成分Ｍに含まれる定位成分を信号成分と仮定するとともに差成分Ｓを雑音成分と仮定すると、和成分Ｍおよび差成分Ｓを利用した処理係数列Ｇeの生成には、雑音成分（差成分Ｓ）を抑圧して信号成分（定位成分）を強調するための数値列（処理係数列Ｇe）を生成する公知の音声強調の技術を同様に適用することが可能である。
【００７８】
処理係数列Ｇeの生成に適用され得る技術としては、ウィナーフィルタ（Wiener filter）を利用した音声強調や、ＭＭＳＥ-ＳＴＳＡ法またはＭＡＰ（maximum a posteriori estimation）推定法を利用した音声強調の技術が例示され得る。ＭＭＳＥ-ＳＴＳＡ法については、Y. Ephraim and D. Malah, "Speech enhancement using a minimum mean-square error short-time spectral amplitude estimator", IEEE ASSP, vol.ASSP-32, no.6, p.1109-1121, Dec. 1984に開示され、ＭＡＰ推定法については、T. Lotter and P. Vary, "Speech enhancement by MAP spectral amplitude estimation using a Super-Gaussian speech model", EURASIP Journal on Applied Signal Processing, vol.2005, no,7, p.1110-1126, July 2005に開示されている。
【００７９】
（８）変形例８  
以上の各形態では、処理係数列Ｇeを利用した演算（数式(7)）で処理係数列Ｇsを生成したが、係数列生成部６０が処理係数列Ｇeの生成と同様の方法で和成分Ｍおよび差成分Ｓから直接的に処理係数列Ｇsを生成する構成も採用され得る。具体的には、処理係数列Ｇeの生成について以上の例示した各数式における和成分Ｍと差成分Ｓとを相互に置換すれば、定位成分の抑圧用の処理係数列Ｇsを生成することが可能である。すなわち、以上の各態様における処理係数列Ｇeと処理係数列Ｇsとは相互に置換され得る。
【００８０】
ただし、処理係数列Ｇeを利用して処理係数列Ｇsを算定する第１実施形態の構成によれば、和成分Ｍや差成分Ｓから直接的に処理係数列Ｇsを算定する処理が不要である。したがって、係数列生成部６０による処理の負荷が軽減されるという利点がある。なお、和成分Ｍおよび差成分Ｓから直接的に処理係数列Ｇsを生成し、処理係数列Ｇsの各係数値Ｇs[k]を数式(7)と同様に所定値（例えば１）から減算することで処理係数列Ｇe（係数値Ｇe[k]）を生成する構成も採用され得る。
【００８１】
（９）変形例９  
以上の各形態では、処理係数列Ｇ（Ｇe，Ｇs）を単位区間毎に生成したが、処理係数列Ｇの生成の周期は任意である。例えば、所定個の単位区間の集合を周期として係数列生成部６０が順次に生成した処理係数列Ｇを、信号処理部３８が当該周期内の複数の単位区間について適用する構成も採用され得る。また、以上の各形態では、各単位区間について生成された処理係数列Ｇをその単位区間の周波数スペクトル（ＬA，ＲA）の処理に適用したが、各単位区間の処理係数列Ｇをその単位区間の経過後の各単位区間の周波数スペクトル（ＬA，ＲA）の処理に適用する構成も採用され得る。
【００８２】
（１０）変形例１０  
以上の各形態では、処理係数列Ｇ（Ｇe，Ｇs）を音響信号ＳINの周波数スペクトル（ＬA，ＲA）に乗算したが、信号処理部３８による処理の内容は適宜に変更される。例えば周波数成分ＬAk(ejω)および周波数成分ＲAk(ejω)から処理係数列Ｇ（Ｇe，Ｇs）を減算することで音響信号ＳOUT（ＳOUT_L，ＳOUT_R）を生成する構成が採用され得る。以上の構成における強調用の処理係数列Ｇeは、第１実施形態の処理係数列Ｇeとは逆に、定位成分のパワーが大きい周波数ｆkの係数値Ｇe[k]ほど０に近い数値に設定される。抑圧用の処理係数列Ｇsについても同様であり、第１実施形態の処理係数列Ｇsとは逆に、定位成分のパワーが大きい周波数ｆkの係数値Ｇs[k]ほど１に近い数値に設定される。
【００８３】
（１１）変形例１１  
以上の各形態においては定位成分の強調用の処理係数列Ｇeおよび抑圧用の処理係数列Ｇsの双方を係数列生成部６０が生成したが、係数列生成部６０が処理係数列Ｇeおよび処理係数列Ｇsの一方のみを生成する構成も採用され得る。したがって、例えば、以上の各形態の第３生成部６６は省略され得る。
【００８４】
（１２）変形例１２  
以上の各形態における係数β（フロアリング係数）を、例えば入力装置１６に対する利用者からの指示に応じて可変に設定する構成も好適である。係数βが小さい（０に近い）ほど定位成分の抑圧の度合が強化されるとともに定位成分の強調の度合が低下する。なお、係数βが大きい場合には、抑圧用の処理係数列Gsのうち定位成分以外の大部分の周波数ｆkの係数値Ｇs[k]が充分に小さい数値となるから、音響信号ＳOUT（ＳOUT_L，ＳOUT_R）の音量が不足する可能性がある。そこで、係数βを大きい数値に設定して定位成分の抑圧を実行する場合、音響信号ＳOUTの音量を増加させる構成が好適である。同様に、係数βが小さい場合には、強調用の処理係数列Ｇeの大部分の周波数ｆkの係数値Ｇe[k]が充分に小さい数値となり得るから、係数βを小さい数値に設定して定位成分の強調を実行する場合、音響信号ＳOUTの音量を増加させる構成が好適に採用される。
【符号の説明】
【００８５】
１００……音響処理装置、１２……信号供給装置、１４……放音装置、１６……入力装置、２２……演算処理装置、２４……記憶装置、３２……変数設定部、３４……周波数分析部、３６……係数設定部、３８……信号処理部、４０……波形合成部、５２……和成分生成部、５４……差成分生成部、６０……係数列生成部、６２……第１生成部、６４……第２生成部、６６……第３生成部、６８……合成部。
【図１】
0005454157.tif 000017
【図２】
0005454157.tif 000018
【図３】
0005454157.tif 000019
【図４】
0005454157.tif 000020
【図５】
0005454157.tif 000021
【図６】
0005454157.tif 000022
【図７】
0005454157.tif 000023
【図８】
0005454157.tif 000024
【図９】
0005454157.tif 000025