(57)【要約】【課題】抑圧対象の音響成分のうちの残存成分を聴感的に知覚され難くする。【解決手段】分離フィルタ生成部３２は、音響信号の特定成分が抑圧されるように各周波数の係数値ｇ(k,τ)が設定された分離フィルタＧ(τ)を、時間軸上で相互に重複する単位期間毎に生成する。移相フィルタ生成部３４は、各周波数における位相を、分離フィルタＧ(τ)のうち当該周波数の係数値ｇ(k,τ)と単位期間毎に変動する変動値ｎ(k,τ)とに応じた移相量θ(k,τ)だけ変化させる移相フィルタＰ(τ)を単位期間毎に生成する。【選択図】図２
2014178578.tif 000002
【特許請求の範囲】
【請求項１】
音響信号の特定成分が抑圧されるように複数の周波数の各々の係数値が設定された分離フィルタを、時間軸上で相互に重複する単位期間毎に生成する分離フィルタ生成手段と、
前記複数の周波数の各々における位相を、前記分離フィルタのうち当該周波数の係数値と単位期間毎に変動する変動値とに応じた移相量だけ変化させる移相フィルタを、単位期間毎に生成する移相フィルタ生成手段と
を具備する音響処理装置。
【請求項２】
前記移相量は、各周波数の前記係数値が前記特定成分を抑圧する度合が大きいほど大きい数値となる制御値と前記変動値との乗算に応じて設定される
請求項１の音響処理装置。
【請求項３】
前記変動値は、単位期間毎に生成された乱数である
請求項１または請求項２の音響処理装置。
【請求項４】
前記分離フィルタと前記移相フィルタとを前記音響信号に作用させる信号処理手段
を具備する
請求項１から請求項３の何れかの音響処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の音響成分が混合された音響信号から特定の音響成分を分離（抽出または抑圧）する技術に関する。
【背景技術】
【０００２】
音響信号の目的音成分と非目的音成分とが混在した音響信号について非目的音成分を抑圧（各成分を分離）するための各種の技術が従来から提案されている。例えば特許文献１には、音響信号の単位期間（フレーム）毎のスペクトルに抑圧係数列（例えばウィナーフィルタ）を作用させることで音響信号のうち目的音成分以外の非目的音成分を抑圧（目的音成分を強調）する技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１１３１９０公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、非目的音成分を確実に抑圧し得る完全な処理係数列を推定することは現実には困難であるから、処理係数列を適用した処理後の音響信号には非目的音成分の一部の成分（以下「残存成分」という）が残存し得る。以上の事情を背景として、本発明は、抑圧対象の音響成分のうちの残存成分を聴感的に知覚され難くすることを目的とする。
【課題を解決するための手段】
【０００５】
以上の課題を解決するために、本発明の音響処理装置は、音響信号の特定成分が抑圧されるように複数の周波数の各々の係数値（例えば係数値ｇ(k,τ)）が設定された分離フィルタを、時間軸上で相互に重複する単位期間毎に生成する分離フィルタ生成手段と、複数の周波数の各々における位相を、分離フィルタのうち当該周波数の係数値と単位期間毎に変動する変動値（例えば変動値ｎ(k,τ)）とに応じた移相量だけ変化させる移相フィルタを、単位期間毎に生成する移相フィルタ生成手段とを具備する。例えば、移相量は、各周波数の係数値が特定成分を抑圧する度合が大きいほど大きい数値となる制御値（例えば制御値Ｆ(k,τ)）と変動値との乗算に応じて設定される。
【０００６】
以上の構成では、音響信号の複数の周波数の各々の移相量が、音響信号の特定成分を抑圧する分離フィルタのうち当該周波数の係数値と単位期間毎に変動する変動値とに応じて設定される。すなわち、音響信号の各周波数の移相量に対する変動値の変動の影響（移相量の変動の度合）が分離フィルタの係数値に応じて周波数毎に制御される。したがって、例えば、音響信号のうち特定成分が優勢な周波数（例えば分離フィルタの係数値が小さい周波数）に対する移相量を、音響信号のうち特定成分以外の音響成分が優勢な周波数に対する移相量と比較して大きく変動させることが可能である。各単位期間は時間軸上で相互に重複する（すなわち、音響信号の各周波数の成分が複数の単位期間にわたる）から、音響信号のうち単位期間毎の移相量の変動が大きい周波数の成分ほど時間軸上の広範囲に分散される。特定成分のエネルギーの最大値は、各周波数成分が時間軸上の広範囲に分散されるほど小さくなる。したがって、以上の構成によれば、分離フィルタの作用の前後で各周波数の位相が維持される構成と比較して、特定成分のうち分離フィルタの係数値では完全には抑圧されない残存成分を聴感的に知覚され難くすることが可能である。
【０００７】
本発明の好適な態様において、変動値は、複数の周波数の各々について個別に設定される。以上の構成によれば、複数の周波数にわたる残存成分を効果的に知覚され難くすることが可能である。変動値の好適例は、単位期間毎に生成された乱数（例えば一様乱数）である。以上の構成によれば、乱数を発生させる簡便な構成により、残存成分を時間軸上に効果的に分散させることが可能である。
【０００８】
本発明の好適な態様に係る音響処理装置は、分離フィルタと移相フィルタとを音響信号に作用させる信号処理手段を具備する。以上の構成によれば、分離フィルタと移相フィルタとを作用させた音響信号を生成することが可能である。
【０００９】
以上の各態様に係る音響処理装置は、音響信号の処理に専用されるＤＳＰ（Digital Signal Processor）などのハードウェア（電子回路）によって実現されるほか、ＣＰＵ（Central Processing Unit）等の汎用の演算処理装置とプログラムとの協働によっても実現される。具体的には、本発明のプログラムは、音響信号の特定成分が抑圧されるように複数の周波数の各々の係数値が設定された分離フィルタを、時間軸上で相互に重複する単位期間毎に生成する分離フィルタ生成処理と、複数の周波数の各々における位相を、分離フィルタのうち当該周波数の係数値と単位期間毎に変動する変動値とに応じた移相量だけ変化させる移相フィルタを、単位期間毎に生成する移相フィルタ生成処理とをコンピュータに実行させる。本発明のプログラムは、コンピュータが読取可能な記録媒体に格納された形態で提供されてコンピュータにインストールされ得る。記録媒体は、例えば非一過性（non-transitory）の記録媒体であり、ＣＤ-ＲＯＭ等の光学式記録媒体（光ディスク）が好例であるが、半導体記録媒体や磁気記録媒体等の公知の任意の形式の記録媒体を包含し得る。また、例えば、本発明のプログラムは、通信網を介した配信の形態で提供されてコンピュータにインストールされ得る。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１実施形態に係る音響処理装置のブロック図である。
【図２】フィルタ生成部の具体的な構成のブロック図である。
【発明を実施するための形態】
【００１１】
＜第１実施形態＞
図１は、本発明の第１実施形態に係る音響処理装置１００のブロック図である。図１に示すように、音響処理装置１００には信号供給装置１２と放音装置１４とが接続される。信号供給装置１２は、音響信号ＳXを音響処理装置１００に供給する。音響信号ＳXは、複数の音響成分（例えば楽音や音声）の混合音の波形を示す時間領域信号である。例えば、複数種の楽器の演奏音（歌唱音等の音声も含意する）の混合音を示す音響信号ＳXが音響処理装置１００に供給される。周囲の音響を収音して音響信号ＳXを生成する収音機器や、可搬型または内蔵型の記録媒体から音響信号ＳXを取得して音響処理装置１００に供給する再生装置や、通信網から音響信号ＳXを受信して音響処理装置１００に供給する通信装置が信号供給装置１２として採用され得る。
【００１２】
第１実施形態の音響処理装置１００は、信号供給装置１２から供給される音響信号ＳXに対する音響処理で音響信号ＳYを生成する信号処理装置である。音響信号ＳYは、音響信号ＳXに包含される複数の音響成分のうち１以上の特定の音響成分（以下「第１音響成分」という）を抽出した音響の波形を示す時間領域信号である。すなわち、音響処理装置１００は、音響信号ＳXのうち第１音響成分以外の１以上の音響成分（以下「第２音響成分」という）を抑圧することで音響信号ＳYを生成する。第１音響成分は、音響信号ＳXのうち強調対象となる目的音成分とも換言され、第２音響成分は、音響信号ＳXのうち抑圧対象となる非目的音成分とも換言される。放音装置１４（例えばスピーカやヘッドホン）は、音響処理装置１００から供給される音響信号ＳYに応じた音波を放射する。なお、音響信号ＳYをデジタルからアナログに変換するＤ/Ａ変換器の図示は便宜的に省略した。
【００１３】
図１に示すように、音響処理装置１００は、演算処理装置２２と記憶装置２４とを具備するコンピュータシステムで実現される。記憶装置２４は、演算処理装置２２が実行するプログラムや演算処理装置２２が使用する各種のデータを記憶する。半導体記録媒体や磁気記録媒体等の公知の記録媒体や複数種の記録媒体の組合せが記憶装置２４として任意に採用され得る。音響信号ＳXを記憶装置２４に記憶した構成（したがって信号供給装置１２は省略され得る）も好適である。
【００１４】
図１の演算処理装置２２は、記憶装置２４に記憶されたプログラムを実行することで、音響信号ＳXから音響信号ＳYを生成するための複数の機能（周波数分析部２０，フィルタ生成部３０，信号処理部４０，波形生成部５０）を実現する。なお、演算処理装置２２の各機能を複数の集積回路に分散した構成や、専用の電子回路（例えばＤＳＰ）が一部の機能を実現する構成も採用され得る。
【００１５】
図１の周波数分析部２０は、音響信号ＳXを時間軸上で区分した複数の単位期間（フレーム）の各々についてスペクトル（複素スペクトル）ＱX(τ)を算定する。任意の１個の単位期間のスペクトルＱX(τ)は、周波数軸上の相異なる周波数に対応する複数の周波数成分（複素数）Ｘ(k,τ)の系列である。記号kは、周波数軸上の任意の１個の周波数（周波数帯域）を指示する変数であり、記号τは、時間軸上の複数の単位期間の何れかを指示する変数である。相前後する各単位期間は時間軸上で相互に部分的に重複する。スペクトルＱX(τ)の生成には例えば短時間フーリエ変換等の公知の周波数分析が任意に利用される。なお、通過帯域が相違する複数の帯域通過フィルタで構成されるフィルタバンクを周波数分析部２０として採用することも可能である。また、音響信号ＳXを記憶装置２４に記憶した構成では、記憶装置２４から周波数分析部２０に音響信号ＳXが供給される。
【００１６】
図１のフィルタ生成部３０は、処理フィルタＭ(τ)を単位期間毎に順次に生成する。本実施形態の処理フィルタＭ(τ)は、音響信号ＳX（スペクトルＱX(τ)）の第２音響成分を抑圧するためのフィルタであり、相異なる周波数（または周波数帯域）に対応する複数の係数値ｍ(k,τ)の系列である。
【００１７】
信号処理部４０は、音響信号ＳXの各単位期間のスペクトルＱX(τ)に当該単位期間の処理フィルタＭ(τ)を作用させることで、音響信号ＳYのスペクトル（複素スペクトル）ＱY(τ)を単位期間毎に順次に生成する。具体的には、信号処理部４０は、以下の数式(1)で表現される通り、音響信号ＳXの各単位期間の周波数成分Ｘ(k,τ)に対し、当該単位期間についてフィルタ生成部３０が生成した処理フィルタＭ(τ)のうち当該周波数成分Ｘ(k,τ)と共通の周波数に対応する係数値ｍ(k,τ)を乗算することで、音響信号ＳYのスペクトルＱY(τ)の周波数成分Ｙ(k,τ)を算定する。
【数１】
2014178578.tif 000003
【００１８】
波形生成部５０は、信号処理部４０が単位期間毎に生成するスペクトルＱY(τ)から時間領域の音響信号ＳYを生成する。具体的には、波形生成部５０は、各単位期間のスペクトルＱY(τ)を例えば短時間逆フーリエ変換で時間領域の信号に変換し、相前後する各単位期間の信号を時間軸上で相互に重複させた状態で加算することにより音響信号ＳYを生成する。波形生成部５０が生成した音響信号ＳYが放音装置１４に供給されることで音響として放射される。
【００１９】
第１実施形態のフィルタ生成部３０が生成する処理フィルタＭ(τ)は、分離フィルタＧ(τ)と移相フィルタＰ(τ)とで構成される。分離フィルタＧ(τ)は、音響信号ＳXの各周波数成分Ｘ(k,τ)の強度を変調するフィルタであり、相異なる周波数に対応する複数の係数値ｇ(k,τ)の系列である。また、移相フィルタＰ(τ)は、音響信号ＳXの各周波数成分Ｘ(k,τ)の位相を変調するフィルタであり、相異なる周波数に対応する複数の係数値ｐ(k,τ)の系列である。処理フィルタＭ(τ)の各係数値ｍ(k,τ)は、以下の数式(2)で表現される通り、分離フィルタＧ(τ)の係数値ｇ(k,τ)と移相フィルタＰ(τ)の係数値ｐ(k,τ)との乗算値に相当する。
【数２】
2014178578.tif 000004
【００２０】
＜フィルタ生成部３０＞
図２は、フィルタ生成部３０の具体的な構成のブロック図である。図２に示すように、フィルタ生成部３０は、分離フィルタ生成部３２と移相フィルタ生成部３４と変動値生成部３６とを含んで構成される。分離フィルタ生成部３２は分離フィルタＧ(τ)を生成し、移相フィルタ生成部３４は移相フィルタＰ(τ)を生成する。数式(2)を参照して前述した通り、分離フィルタ生成部３２が生成した分離フィルタＧ(τ)と移相フィルタ生成部３４が生成した移相フィルタＰ(τ)とに応じた処理フィルタＭ(τ)が信号処理部４０による第２音響成分の抑圧（第１音響成分の強調）に適用される。
【００２１】
分離フィルタ生成部３２は、分離フィルタＧ(τ)を単位期間毎に順次に生成する。本実施形態の分離フィルタ生成部３２は、音響信号ＳXの第２音響成分が抑圧されるように分離フィルタＧ(τ)の各係数値ｇ(k,τ)を設定する。具体的には、分離フィルタＧ(τ)は、各係数値ｇ(k,τ)が以下の数式(3)で表現されるウィナー（Wiener）フィルタである。
【数３】
2014178578.tif 000005
数式(3)の記号|Ｘ(k,τ)|2は周波数成分Ｘ(k,τ)のパワーであり、記号|Ｅ(k,τ)|2は、音響信号ＳXから推定される第１音響成分のパワーである。第１音響成分のパワー|Ｅ(k,τ)|2の推定には公知の技術が任意に採用される。数式(1)および数式(2)から理解される通り、係数値ｇ(k,τ)は、周波数成分Ｘ(k,τ)の強度に作用し、周波数成分Ｘ(k,τ)の利得（スペクトルゲイン）を意味する。数式(3)から理解される通り、第２音響成分が第１音響成分に対して優勢な周波数（例えば第２音響成分の強度が第１音響成分の強度を上回る周波数）の係数値ｇ(k,τ)ほど小さい数値となるように、各係数値ｇ(k,τ)は０以上かつ１以下の範囲内で音響信号ＳXの音響特性に応じて可変に設定される。
【００２２】
ところで、分離フィルタＧ(τ)は以上のように生成されるから、音響信号ＳXのスペクトルＱX(τ)に分離フィルタＧ(τ)を作用させるだけの構成（以下「対比例」という）でも、音響信号ＳXから第２音響成分を抑圧することは可能である（Ｙ(k,τ)＝ｇ(k,τ)Ｘ(k,τ)）。対比例では、音響信号ＳXのスペクトルＱX(τ)の位相が分離フィルタＧ(τ)の作用後にも維持される。しかし、第１音響成分のパワー|Ｅ(k,τ)|2を音響信号ＳXから完全に推定することは現実には困難であるから、対比例のもとで音響信号ＳXに分離フィルタＧ(τ)を作用させた音響には、第２音響成分の一部（残存成分）が残存し得る。第１実施形態の移相フィルタＰ(τ)は、分離フィルタＧ(τ)では抑圧し切れない残存成分を聴感的に知覚され難くするためのフィルタである。
【００２３】
図２の移相フィルタ生成部３４は、音響信号ＳXに作用させる移相フィルタＰ(τ)を単位期間毎に順次に生成する。移相フィルタＰ(τ)は、前述の通り、相異なる周波数（または周波数帯域）に対応する複数の係数値ｐ(k,τ)の系列である。具体的には、移相フィルタＰ(τ)の各係数値ｐ(k,τ)は、以下の数式(4)で表現される通り、スペクトルＱX(τ)における各周波数成分Ｘ(k,τ)の位相を移相量θ(k,τ)だけ変化させるように設定される。
【数４】
2014178578.tif 000006
【００２４】
数式(4)の記号ｊは虚数単位を意味する。第１実施形態の移相量θ(k,τ)は、制御値Ｆ(k,τ)と変動値ｎ(k,τ)とに応じて可変に設定される。具体的には、移相フィルタ生成部３４は、以下の数式(5)で表現される通り、制御値Ｆ(k,τ)と変動値ｎ(k,τ)とを乗算することで移相量θ(k,τ)を算定する。
【数５】
2014178578.tif 000007
【００２５】
図２の変動値生成部３６は、数式(5)の変動値ｎ(k,τ)を複数の周波数の各々について単位期間毎に算定する。変動値ｎ(k,τ)は、単位期間毎に刻々と変動する数値であり、例えば単位期間毎に発生する乱数である。例えば、変動値生成部３６は、所定範囲(-π≦ｎ(k,τ)≦π)内の各数値が略一様な確率で発生する一様乱数を変動値ｎ(k,τ)として算定する。したがって、変動値ｎ(k,τ)は正数または負数に設定され得る。
【００２６】
数式(5)の制御値Ｆ(k,τ)は、複数の周波数の各々について単位期間毎に設定される。具体的には、任意の１個の単位期間における任意の１個の周波数の制御値Ｆ(k,τ)は、当該単位期間の分離フィルタＧ(τ)のうち当該周波数の係数値ｇ(k,τ)に応じて設定される。具体的には、移相フィルタ生成部３４は、係数値ｇ(k,τ)により音響信号ＳXが抑圧される度合が大きい（係数値ｇ(k,τ)が小さい）ほど制御値Ｆ(k,τ)を大きい数値に設定する。例えば、前掲の数式(6)から理解される通り、所定値（数式(6)の例示では１）から係数値ｇ(k,τ)を減算した数値が制御値Ｆ(k,τ)として好適である。
【００２７】
以上の説明から理解される通り、本実施形態では、各周波数の移相量θ(k,τ)に対する変動値ｎ(k,τ)の変動の影響（移相量θ(k,τ)の変動の度合）が、分離フィルタＧ(τ)の係数値ｇ(k,τ)（制御値Ｆ(k,τ)）に応じて単位期間毎および周波数毎に制御される。係数値ｇ(k,τ)が小さい（第２音響成分が優勢である）ほど制御値Ｆ(k,τ)は大きな数値に設定されるから、音響信号ＳXのうち第２音響成分が優勢な周波数の移相量θ(k,τ)は、第１音響成分が優勢な周波数の移相量θ(k,τ)と比較して大きく変動するという傾向がある。
【００２８】
音響信号ＳXのうち時間軸上の特定の時点に存在する特定の周波数の残存成分に着目する。各単位期間は時間軸上で相互に重複する（すなわち、音響信号の各周波数の成分が複数の単位期間にわたる）から、残存成分は、複数の単位期間に重複して存在する。各単位期間の残存成分は、当該単位期間の移相フィルタＰ(τ)の係数値ｐ(k,τ)の作用で移相量θ(k,τ)だけ位相が変化する。移相量θ(k,τ)は、単位期間毎に変動する変動値ｎ(k,τ)に応じて設定されるから、各単位期間の残存成分の位相は単位期間毎に相異なる移相量θ(k,τ)だけ変動する。すなわち、各単位期間の残存成分は、時間軸上の相異なる時点に移動する。波形生成部５０では、各単位期間が相互に重複した状態で加算されるから、移相フィルタＰ(τ)の作用前に時間軸上の１点に存在していた残存成分は、移相フィルタＰ(τ)の作用後には時間軸上の複数の時点に分散される。すなわち、音響信号ＳXの残存成分の強度（エネルギー）が時間軸上で分散される。移相フィルタＰ(τ)は、音響信号ＳXのうち第２音響成分の残存成分に残響を付与する要素とも換言される。
【００２９】
各周波数成分Ｘ(k,τ)の移相量θ(k,τ)は、第２音響成分が優勢な周波数ほど大きいから、音響信号ＳXに移相フィルタＰ(τ)を作用させることで、音響信号ＳXの複数の周波数成分Ｘ(k,τ)のうち第２音響成分が優勢な周波数の周波数成分Ｘ(k,τ)ほど時間軸上の広範囲に分散される。したがって、第１実施形態によれば、音響信号ＳXの位相が分離フィルタＧ(τ)の作用後にも維持される対比例と比較して、第２音響成分の残存成分を聴感的に知覚され難くすることができる。
【００３０】
＜第２実施形態＞
本発明の第２実施形態を以下に説明する。なお、以下に例示する各構成において作用や機能が第１実施形態と同等である要素については、第１実施形態の説明で参照した符号を流用して各々の詳細な説明を適宜に省略する。
【００３１】
第１実施形態では、音響信号ＳXの第１音響成分を抽出した音響信号ＳYを生成した。第２実施形態の音響処理装置１００は、音響信号ＳY1および音響信号ＳY2を音響信号ＳXから選択的に生成する。音響信号ＳY1は、第１実施形態の音響信号ＳYと同様に、第１音響成分を目的音成分として音響信号ＳXから抽出（第２音響成分を非目的音成分として抑圧）した時間領域信号であり、音響信号ＳY2は、第２音響成分を目的音成分として音響信号ＳXから抽出（第１音響成分を非目的音成分として抑圧）した時間領域信号である。音響信号ＳY1または音響信号ＳY2の何れを生成するかは、例えば利用者からの指示に応じて選択される。
【００３２】
第２実施形態のフィルタ生成部３０は、音響信号ＳXの第１音響成分を抽出するための処理フィルタＭ1(τ)と音響信号ＳXの第２音響成分を抽出するための処理フィルタＭ2(τ)とを選択的に生成する。処理フィルタＭi(τ)（ｉ＝１,２）は、相異なる周波数に対応する複数の係数値ｍi(k,τ)の系列であり、分離フィルタＧi(τ)と移相フィルタＰi(τ)とを含んで構成される。分離フィルタＧi(τ)は、音響信号ＳXの各周波数成分Ｘ(k,τ)の強度を変調する複数の係数値ｇi(k,τ)で構成され、移相フィルタＰi(τ)は、各周波数成分Ｘ(k,τ)の位相を変調する複数の係数値ｐi(k,τ)で構成される。
【００３３】
信号処理部４０は、音響信号ＳXの各単位期間のスペクトルＱX(τ)に当該単位期間の処理フィルタＭi(τ)を作用させることで、音響信号ＳYiのスペクトルＱYi(τ)を単位期間毎に順次に生成する。具体的には、信号処理部４０は、以下の数式(7)の演算で音響信号ＳYiの周波数成分Ｙi(k,τ)を算定する。
【数６】
2014178578.tif 000008
【００３４】
処理フィルタＭi(τ)の各係数値ｍi(k,τ)は、以下の数式(8)で表現される通り、係数値ｇi(k,τ)と係数値ｐi(k,τ)との乗算値に相当する。
【数７】
2014178578.tif 000009
【００３５】
分離フィルタＧ1(τ)の各係数値ｇ1(k,τ)は、第１実施形態の係数値ｇ(k,τ)と同様に前掲の数式(3)で算定され、周波数成分Ｘ(k,τ)にて第２音響成分が優勢であるほど小さい数値となる。したがって、数式(7)および数式(8)から理解される通り、処理フィルタＭ1(k,τ)は、第１実施形態の処理フィルタＭ(k,τ)と同様に、音響信号ＳXの第２音響成分を抑圧するとともに、第２音響成分の残存成分（分離フィルタＧ1(τ)では抑圧し切れない成分）が移相により聴感的に知覚され難くなるように作用する。
【００３６】
他方、分離フィルタＧ2(τ)の各係数値ｇ2(k,τ)は、周波数成分Ｘ(k,τ)において第１音響成分が優勢であるほど小さい数値に設定される。例えば、係数値ｇ2(k,τ)は、以下の数式(9)の演算で係数値ｇ1(k,τ)に応じて算定される。
【数８】
2014178578.tif 000010
以上の説明から理解される通り、処理フィルタＭ2(k,τ)は、音響信号ＳXの第１音響成分を抑圧するとともに、第１音響成分の残存成分（分離フィルタＧ2(τ)では抑圧し切れない成分）が移相により聴感的に知覚され難くなるように作用する。。
【００３７】
第２実施形態においても第１実施形態と同様の効果が実現される。また、第２実施形態では、音響信号ＳXから第１音響成分を抽出した音響信号ＳY1に加え、音響信号ＳXから第２音響成分を抽出した音響信号ＳY2を生成することができる。
【００３８】
＜変形例＞
以上に例示した各形態には様々な変形が加えられる。具体的な変形の態様を例示すれば以下の通りである。なお、以下の例示から２以上の態様を任意に選択して組合せてもよい。
【００３９】
（１）分離フィルタＧ(τ)の形式は前述の各形態での例示に限定されない。例えば、各係数値ｇ(k,τ)が以下の数式(10)で表現されるウィナーフィルタを分離フィルタＧ(τ)として利用することも可能である。
【数９】
2014178578.tif 000011
数式(10)の記号max{ }は、括弧内の最大値を選択する演算（すなわち差分|Ｘ(k,τ)|2−|Ｉ(k,τ)|2を非負値に制限する演算）を意味する。記号|Ｉ(k,τ)|2は、音響信号ＳXから推定される第２音響成分のパワーである。第２音響成分のパワー|Ｉ(k,τ)|2の推定には公知の技術が任意に採用される。
【００４０】
また、第２実施形態では、分離フィルタＧ1(τ)の係数値ｇ1(k,τ)から分離フィルタＧ2(τ)の各係数値ｇ2(k,τ)を算定したが、第１音響成分を抑圧するための分離フィルタＧ2(τ)の各係数値ｇ2(k,τ)を、以下に例示される数式(11)または(12)の演算で算定することも可能である。
【数１０】
2014178578.tif 000012
以上の説明から理解される通り、分離フィルタＧ(τ)は、音響信号ＳXの特定成分（第１音響成分または第２音響成分）が抑圧されるように複数の周波数の各々の係数値ｇ(k,τ)が設定されたフィルタとして包括的に表現される。
【００４１】
（２）変動値生成部３６が生成する変動値ｎ(k,τ)は一様乱数に限定されない。例えば各数値の発生確率が正規分布に従う正規乱数を変動値ｎ(k,τ)として利用することも可能である。もっとも、変動値ｎ(k,τ)は乱数に限定されない。例えば、単位期間毎に所定量ずつ増加または減少する数値を変動値ｎ(k,τ)として利用することも可能である。また、前述の各形態では、周波数毎に個別に変動値ｎ(k,τ)を設定したが、複数の周波数にわたり共通の変動値ｎ(τ)を単位期間毎に生成することも可能である。以上の説明から理解される通り、前述の各形態の変動値ｎ(k,τ)は、単位期間毎に変動する数値として包括的に表現される。
【００４２】
（３）第１実施形態では、所定値（１）から係数値ｇ(k,τ)を減算した数値を制御値Ｆ(k,τ)として算定した。以上の構成では、係数値ｇ(k,τ)が１に近い場合（第１音響成分が優勢である場合）にも制御値Ｆ(k,τ)が有意な正数に設定され、結果的に周波数成分Ｘ(k,τ)が移相される。以上の傾向を考慮すると、第１音響成分が優勢である周波数成分Ｘ(k,τ)については位相の変動を積極的に抑制する構成が好適である。例えば、係数値ｇ(k,τ)が所定の閾値ＴHを上回る程度に周波数成分Ｘ(k,τ)にて第１音響成分が優勢である場合には、移相フィルタ生成部３４は、制御値Ｆ(k,τ)を強制的に０に設定する。前掲の数式(5)から理解される通り、制御値Ｆ(k,τ)が０の周波数成分については移相量θ(k,τ)も０となる。すなわち、係数値ｇ(k,τ)が所定の閾値ＴHを上回る周波数成分Ｘ(k,τ)については位相が変更されない。したがって、以上の構成によれば、音響信号ＳXの第１音響成分の波形歪を抑制できるという利点がある。なお、第２実施形態においても同様の構成が採用される。具体的には、係数値ｇi(k,τ)が所定の閾値ＴHi（閾値ＴH1と閾値ＴH2との異同は不問）を上回る周波数の制御値Ｆi(k,τ)（移相量θi(k,τ)）が０に設定される。
【００４３】
（４）第２実施形態においては、音響信号ＳY1および音響信号ＳY2を選択的に生成したが、音響信号ＳY1および音響信号ＳY2を並列に生成することも可能である。以上の構成によれば、例えば音響信号ＳY1および音響信号ＳY2の各々に別個の処理(例えば残響付与等の各種の効果付与)を実行したうえで両者を混合するといった特徴的な音響処理が実現される。
【００４４】
（５）携帯電話機等の端末装置と通信するサーバ装置で音響処理装置１００を実現することも可能である。例えば、音響処理装置１００は、端末装置から受信した音響信号ＳXから音響信号ＳYを生成して端末装置に送信する。なお、音響信号ＳXのスペクトルＱX(τ)を端末装置から受信する構成（例えば端末装置が周波数分析部２０を具備する構成）では音響処理装置１００から周波数分析部２０が省略され、スペクトルＱY(τ)を端末装置に送信する構成（例えば端末装置が波形生成部５０を具備する構成）では音響処理装置１００から波形生成部５０が省略される。また、音響処理装置１００が生成した処理フィルタＭ(τ)を端末装置等の外部装置に転送して音響信号ＳXの処理に適用する構成では、信号処理部４０も音響処理装置１００から省略され得る。
【符号の説明】
【００４５】
１００……音響処理装置、１２……信号供給装置、１４……放音装置、２０……周波数分析部、２２……演算処理装置、２４……記憶装置、３０……フィルタ生成部、３２……分離フィルタ生成部、３４……移相フィルタ生成部、３６……変動値生成部、４０……信号処理部、５０……波形生成部
【図１】
2014178578.tif 000013
【図２】
2014178578.tif 000014