(57)【特許請求の範囲】
【請求項１】
複数の音源それぞれに対応する音源の方向又は位置を特定する特定部と、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換する変換処理と、前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記特定された複数の音源それぞれに対応する音源の方向又は位置の中から最適な２つの音源の方向又は位置を選択する選択処理と、前記選択された最適な２つの音源の方向又は位置に基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記特定された各音源の方向又は位置毎でかつ周波数帯域毎の信号に分解する分解処理と、を実行する処理部と、
を備える音源信号変換装置。
【請求項２】
前記選択処理は、
各周波数帯域の各フレームに属する周波数の信号について、前記特定された複数の音源それぞれに対応する音源の方向又は位置の中から任意に選択された２つの音源の方向又は位置の組み合わせの夫々がどの程度確かかを示す第１確率と、時間軸方向に隣接するフレーム間で２つの音源の方向又は位置の組み合わせが遷移する確率である第２確率と、から求まる尤度を最大にする２つの音源の方向又は位置の組み合わせを推定し、各フレームに対して前記推定された２つの音源の方向又は位置を選択する、請求項１に記載の音源信号変換装置。
【請求項３】
前記処理部は、
前記特定された各音源の方向又は位置の組み合わせそれぞれに対して付与されるラベルを生成するラベル生成処理と、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記特定された各音源の方向又は位置毎に生成する生成処理と、
をさらに実行し、
前記選択処理は、前記周波数帯域毎に、
前記フレームのそれぞれに、前記ラベル生成処理にて生成される全てのラベルをラベル候補として割り振る処理と、
前記割り振られたラベル候補毎に、前記ラベル候補に対応する周波数帯域毎のフレームに属する周波数の信号が、前記ラベル候補に対応する前記音源の方向又は位置の組み合わせの夫々がどの程度確かかを示すラベル確率を計算する処理と、
前記時間軸方向に隣接するフレームに割り振られたラベル候補の全ての組合せ間のラベル遷移確率を計算する処理と、
前記ラベル確率及び前記ラベル遷移確率から求まる尤度を最大にするラベル候補の組合せを推定して前記フレーム毎に最適なラベルを付与する処理と、
を実行する、請求項１または２に記載の音源信号変換装置。
【請求項４】
前記処理部は、
前記音源の方向又は位置が第１の既定時間内に第１の既定数以上変化しないように、及び前記分解処理での分解結果の振幅が第２の既定時間内に第２の既定値以上変化しないように、前記ラベル遷移確率を計算する処理を実行する、請求項３に記載の音源信号変換装置。
【請求項５】
前記処理部は、
同時に変化する前記音源の方向又は位置が１つまでに制限されるように、前記ラベル遷移確率を計算する処理を実行する、請求項３又は４に記載の音源信号変換装置。
【請求項６】
前記処理部は、
前記ラベル確率及び前記ラベル遷移確率から求まる尤度を最大にする前記ラベル候補の組合せを、確立伝搬法のアルゴリズムに従って推定する処理を実行する、請求項３乃至５のいずれかに記載の音源信号変換装置。
【請求項７】
前記特定部は、外部より設定された複数の音源それぞれに対応する音源の方向又は位置を特定する請求項１乃至６のいずれかに記載の音源信号変換装置。
【請求項８】
前記処理部はさらに、
前記変換処理にて得られた左右各チャネルそれぞれの複数種の周波数領域毎の信号に対して、周波数帯毎および前記フレーム毎に、音源の方向又は位置及び前記音源の方向又は位置が単一か否かを表す評価値を算出する処理と、
前記算出された音源の方向又は位置とその評価値とに基づき、前記音源の方向又は位置毎に評価値を積算することによりヒストグラムを作成する処理と、
前記作成されたヒストグラムをディスプレイに提示する処理と、
を実行する請求項７に記載の音源信号変換装置。
【請求項９】
複数の音源それぞれに対応する音源の方向又は位置を特定する特定部と、処理部とを有する音源信号変換装置に用いられる音源信号変換方法であって、前記処理部が、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換し、
前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記特定された複数の音源それぞれに対応する音源の方向又は位置の中から最適な２つの音源の方向又は位置を選択し、
前記選択された最適な２つの音源の方向又は位置に基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記特定された各音源の方向又は位置毎でかつ周波数帯域毎の信号に分解し、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記定位毎に生成する、音源信号変換方法。
【請求項１０】
複数の音源それぞれに対応する音源の方向又は位置を特定する特定部と、処理部とを有する音源信号変換装置として用いられるコンピュータに、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換するステップと、
前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記特定された複数の音源それぞれに対応する音源の方向又は位置の中から最適な２つの音源の方向又は位置を選択するステップと、
前記選択された最適な２つの音源の方向又は位置に基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記特定された各音源の方向又は位置毎でかつ周波数帯域毎の信号に分解するステップと、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記定位毎に生成するステップと、
を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音源からのステレオ信号を定位毎の信号に変換する音源信号変換装置、方法、及びプログラムに関する。
【背景技術】
【０００２】
今日、コンパクトディスクなどで提供されている音楽は、スタジオなどで楽器やボーカルなどパート毎に録音を行い、これらの録音した音源をミキシングにより左右２チャネルのステレオ音楽信号を合成して作成される。この際、各音源には異なる定位を与えることが一般的である。
【０００３】
このミキシングされたステレオ信号から、元の音源信号を推定、分離する需要がある。例えば、ボーカルとピアノ音を消去した音を再生し、そこに自分のピアノ演奏や歌声を重ねて楽しむ用途などがある。また音源を分離できれば採譜がより容易になる。
【０００４】
ステレオ信号から３つ以上の音源を分離することは原理的にはできないので、近似的に分離する従来技術が知られている。一般的には、ステレオ信号に対して短時間フーリエ変換などの時間周波数変換を実行し、周波数帯、フレーム毎に振幅比から定位を求める。そして、この定位が予め定めた範囲に入っている信号のみを抽出したり削除したり（バイナリマスク処理）した後に、逆短時間フーリエ変換などの逆時間周波数変換を行い、時間領域の波形に戻すことで音源を分離する。
【０００５】
また、バイナリマスク処理による信号の不連続性の影響を軽減する従来技術も知られている（例えば特許文献１）。この従来技術では、定位情報をクラスタリングして代表値を求める。この代表値は、各音源の定位の推定値である。代表値毎に定位までの距離が近いほど大きな重み付け係数を生成し、この重み付け係数を掛け合わせて音源の分離を行う。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許４７６７２４７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述した従来技術では、各周波数帯、各フレームで鳴っている音源は１つであることを仮定している。しかし、この仮定は正しくない場合がある。例えば、ボーカル、鍵盤楽器（ピアノ等）、ギター、ベース、ドラム等で演奏される典型的なポピュラー音楽では、メロディを担うボーカルの高調波構造の一部周波数帯の信号成分（特に基本周波数成分）は、鍵盤楽器、ギターの高調波周波数信号成分やベースの高次の高調波周波数信号成分、スネアドラム等の音に含まれるノイズ成分などと頻繁に重複する。また短時間フーリエ変換の結果算出される周波数信号成分間の信号漏洩によっても、或る周波数信号成分において重複が発生する。
【０００８】
この場合、従来技術では、ある２つの定位の信号が本来とは異なる別の定位として認識され、音源を正しく分離することができないという課題があった。例えばステレオ信号の左側視聴位置に近い場所でピアノ音が鳴っており、右側視聴位置に近い位置でギター音が同じくらいの強さで鳴っている場合には、従来技術では、両者の音の定位が中央位置付近にあるように計算されてしまう。この結果、従来技術で音源の分離を行うと、本来中央位置付近で鳴っている音（例えばボーカルやベースなど）に関係ない周波数信号成分が混入してしまい、一方、左側視聴位置に近い位置で鳴っているピアノ音と、右側視聴位置に近い位置で鳴っているギター音において、本来あるべき周波数信号成分が欠落してしまうという課題があった。
【０００９】
そこで、本発明は、最大で２つまでの音源の適切な組合せを推測することにより、同一周波数帯でほぼ同一時刻に２つの定位で発音がある場合でも音源を適切に分離可能とすることを目的とする。
【課題を解決するための手段】
【００１０】
態様の一例では、複数の音源それぞれに対応する音源の方向又は位置を特定する特定部と、複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換する変換処理と、前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記特定された複数の音源それぞれに対応する音源の方向又は位置の中から最適な２つの音源の方向又は位置を選択する選択処理と、前記選択された最適な２つの音源の方向又は位置に基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記特定された各音源の方向又は位置毎でかつ周波数帯域毎の信号に分解する分解処理と、を実行する処理部と、を備える。
【発明の効果】
【００１１】
本発明によれば、複数の音源の適切な組合せを推測することにより、同一周波数帯でほぼ同一時刻に複数の定位で発音がある場合でも音源を適切に分離することが可能となる。
【図面の簡単な説明】
【００１２】
【図１】音源信号変換装置の実施形態のブロック図である。
【図２】定位の説明図である。
【図３】ミキシングの説明図である。
【図４】音源信号変換装置の機能を実現可能なコンピュータのハードウェア構成例を示す図である。
【図５】音源信号変換処理の例を示すフローチャートである。
【図６】定位設定処理の例を示すフローチャートである。
【図７】ヒント情報の表示例を示す図である。
【図８】定位の組合せとラベルの対応関係の例を示す図である。
【図９】本実施形態におけるデータ構成例を示す図である。
【図１０】本実施形態における確率モデルの説明図である。
【図１１】選択処理の例を示すフローチャートである。
【図１２】確率伝播法による処理の例を示すフローチャートである。
【発明を実施するための形態】
【００１３】
以下、本発明を実施するための形態について図面と数式を参照しながら詳細に説明する。図１は、音源信号変換装置１００の実施形態のブロック図である。音源信号変換装置１００は、変換部１０１、定位設定部１０２、ラベル生成部１０３、分解部１０４、選択部１０５、及び生成部１０６を備える。
【００１４】
変換部１０１は、ステレオ信号の左チャネル信号Ｌ（ｔ）及び右チャネル信号Ｒ（ｔ）を、それぞれ時間領域の信号から時間周波数領域の信号に変換する。変換処理アルゴリズムとしては、短時間フーリエ変換又はウェーブレット変換を採用することができる。以降の説明では、短時間フーリエ変換を前提に説明を行う。この変換処理の結果、左チャネル信号Ｌ（ｔ）に対応する複素スペクトログラムと、右チャネル信号Ｒ（ｔ）に対応する複素スペクトログラムが、それぞれ生成される。ここで、複素スペクトログラムとは、フレーム（時間区間）毎における周波数帯毎の周波数信号成分を複素数で表現したものである。
【００１５】
定位設定部１０２は、後述する入力装置によりユーザに定位を設定させる。本実施形態では、ユーザが設定した定位をもとに、従来よりも適切な音源信号変換を可能とするものである。
【００１６】
ラベル生成部１０３は、定位設定部１０２で設定した定位から選択される最大で２つまでの定位の組合せ毎にそれぞれ付与されるラベルを生成する。これはステレオ信号に対して定位が２つの場合は、逆問題を一意に解くことができるからである。
【００１７】
選択部１０５は、変換部１０１にて得られるステレオ２チャネルの複素スペクトログラムに基づいて、周波数帯毎およびフレーム毎に、ラベル生成部１０３で生成したラベルのうち最適なラベルを付与する。
【００１８】
分解部１０４は、選択部１０５で付与した最適なラベルに基づいて、周波数帯毎およびフレーム毎に、各チャネルの複素スペクトログラムを定位毎の複素スペクトログラムに分解して抽出する。
【００１９】
生成部１０６は、分解部１０４で抽出した定位毎の複素スペクトログラムをその定位毎の時間領域の信号に変換する。
【００２０】
以上の構成を有する音源信号変換装置１００の動作に説明以下に順次説明する。図２は、定位の説明図である。定位とは、視聴者に対する楽器の配置に相当し、一般的には、図２に示すような視聴者の前方の半円上に各音源が配置されるモデルが想定されることが多い。図２において、黒丸はそれぞれの楽器の定位を表している。また、ＬとＲは左耳と右耳の視聴位置に相当する。各音源の信号がミキシングされることにより、ＬとＲの位置でステレオ試聴を行うと、図２に示すように各音源が定位しているように知覚される。
【００２１】
図３は、ミキシングの説明図である。図３において、Ｓｉはｉ番目の信号源、ＯはＬとＲの中点で半円の中心、θｉは定位を表す角度、ＤＬｉは信号源ＳｉからＬまでの距離、ＤＲｉは信号源ＳｉからＲまでの距離を表す。音の振幅は音源からの距離の逆数に比例する。ここで、ＤＬｉの逆数に相当する乗数をＬｉ、ＤＲｉの逆数に相当する乗数をＲｉとする。半円の半径は特に定められていないので、下記数１式を満たすように正規化する。
【００２２】
【数１】
0006613737.tif 000002
【００２３】
また、θｉとＬｉ、Ｒｉの関係は、下記数２式及び数３式に示すようになる。
【００２４】
【数２】
0006613737.tif 000003
【００２５】
【数３】
0006613737.tif 000004
【００２６】
時刻ｔにおける音源Ｓｉの信号をＳｉ（ｔ）とすると、時刻ｔにおける左耳視聴位置Ｌにおける左チャネル信号Ｌ（ｔ）は下記数４式を、右耳視聴位置Ｒにおける右チャネル信号Ｒ（ｔ）は下記数５式を計算することで、ミキシングが行われる。
【００２７】
【数４】
0006613737.tif 000005
【００２８】
【数５】
0006613737.tif 000006
【００２９】
この左チャネル信号Ｌ（ｔ）と右チャネル信号Ｒ（ｔ）のステレオ信号をヘッドホンで視聴すると、図２に示すように各音源が定位しているように知覚される。
【００３０】
図４は、図１の音源信号変換装置１００の機能をソフトウェア処理として実現可能なコンピュータのハードウェア構成の一例を示す図である。図４に示されるコンピュータは、ＣＰＵ（中央演算処理装置）４０１、メモリ４０２、入力装置４０３、出力装置４０４、外部記憶装置４０５、可搬記録媒体４０９が挿入される可搬記録媒体駆動装置４０６、及び通信インターフェース４０７を有し、これらがバス４０８によって相互に接続された構成を有する。同図に示される構成は図１の音源信号変換装置１００の機能を実現可能なコンピュータの一例であり、そのようなコンピュータはこの構成に限定されるものではない。
【００３１】
ＣＰＵ４０１は、当該コンピュータ全体の制御を行う。メモリ４０２は、プログラムの実行、データ更新等の際に、外部記憶装置４０５（或いは可搬記録媒体４０９）に記憶されているプログラム又はデータを一時的に格納するＲＡＭ等のメモリである。ＣＰＵ４０１は、プログラムをメモリ４０２に読み出して実行することにより、全体の制御を行う。
【００３２】
入力装置４０３は、ユーザによるキーボードやマウス等による入力操作を検出し、その検出結果をＣＰＵ４０１に通知する。また、この入力装置４０３によりユーザは、複数種の定位の組み合わせを入力する。
【００３３】
出力装置４０４は、ＣＰＵ４０１の制御によって送られてくるデータを表示装置や印刷装置に出力する。
【００３４】
外部記憶装置４０５は、例えばハードディスク記憶装置である。主に各種データやプログラムの保存に用いられる。
【００３５】
可搬記録媒体駆動装置４０６は、ＳＤカード、コンパクトフラッシュや、ＣＤ−ＲＯＭ、ＤＶＤ、光ディスク等の可搬記録媒体４０９を収容するもので、外部記憶装置４０５の補助の役割を有する。
【００３６】
通信インターフェース４０７は、例えばＬＡＮ（ローカルエリアネットワーク）又はＷＡＮ（ワイドエリアネットワーク）の通信回線を接続するための装置である。
【００３７】
本実施形態による図１の音源信号変換装置１００の機能は、図５、図６、図１１、及び図１２のフローチャート等で実現される機能を搭載したプログラムをＣＰＵ４０１が実行することで実現される。そのプログラムは、例えば外部記憶装置４０５や可搬記録媒体４０９に記録して配布してもよく、或いは通信インターフェース４０７によりネットワークから取得できるようにしてもよい。
【００３８】
図５は、図４のＣＰＵ４０１が実行する音源信号変換処理の例を示すフローチャートである。この処理は、図４において、ＣＰＵ４０１が例えば外部記憶装置４０５からメモリ４０２に、音源信号変換処理プログラムを読み出して実行する動作である。
【００３９】
ＣＰＵ４０１はまず、図１の変換部１０１の機能に対応する時間周波数変換処理を実行する（ステップＳ１０１）。この処理では、ステレオ信号の左チャネル信号Ｌ（ｔ）及び右チャネル信号Ｒ（ｔ）が、それぞれ時間領域の信号から時間周波数領域の信号である複素スペクトログラムに変換される。この結果得られる左チャネルの複素スペクトログラムＳＰＬ（ｍ，ｂ）を下記数６式に、右チャネルの複素スペクトログラムＳＰＲ（ｍ，ｂ）を下記数７式に示す。
【００４０】
【数６】
0006613737.tif 000007
【００４１】
【数７】
0006613737.tif 000008
【００４２】
ここで、ｍはフレーム番号を、ｂは周波数帯のビン番号を、ＡＬ（ｍ，ｂ）、ＡＲ（ｍ，ｂ）は左チャネル及び右チャネルの各振幅を、ΦＬ（ｍ，ｂ）、ΦＲ（ｍ，ｂ）は左チャネル及び右チャネルの各位相を、それぞれ表す。
【００４３】
次に、ＣＰＵ４０１は、図１の定位設定部１０２の機能に対応する定位設定処理を実行する（ステップＳ１０２）。図６は、図５のステップＳ１０２の定位設定処理の詳細例を示すフローチャートである。
【００４４】
図６のフローチャートにおいて、ＣＰＵ４０１はまず、図５のステップＳ１０１で生成した１組の複素スペクトログラムの全フレーム、全周波数帯に対して、定位と評価値の計算処理を実行する（ステップＳ２０１）。複素スペクトログラムは線形性が保たれているので、前述した数４式と数５式の関係も維持され、下記の数８式と数９式の関係も成り立つ。
【００４５】
【数８】
0006613737.tif 000009
【００４６】
【数９】
0006613737.tif 000010
【００４７】
ここで、ＳＰｉ（ｍ，ｂ）はＳｉ（ｔ）を時間周波数変換したものである。あるフレームｍ、ある周波数帯ｂで単一の定位の音しかない場合は、位相が等しく、かつ左チャネル及び右チャネルの各振幅の比が、ミキシング時の重み付け係数と等しくなる。よって、下記数１０式により、定位を表す角度θ（ｍ，ｂ）（図３のθｉに対応する）が計算できる。
【００４８】
【数１０】
0006613737.tif 000011
【００４９】
またこのとき、単一の定位か否かを表す評価値も計算できる。下記数１１式は、評価値の計算方法の一例である。
【００５０】
【数１１】
0006613737.tif 000012
【００５１】
ここで、Ｋｐ（＞０）は調整用の係数で、値を小さくすると位相が等しいか否かの判定が緩くなる。また、Ｄは位相差を計算する関数である。位相は±πの範囲に折りたたまれているので、この関数は折りたたみを考慮した差分値を計算する。なお、評価値は、振幅値が所定の閾値より小さい場合はゼロにするなどすると、なおよい。
【００５２】
次に、ＣＰＵ４０１は、図６のステップＳ２０１で計算した定位とその評価値を用いて、ヒストグラム作成処理を実行する（ステップＳ２０２）。ヒストグラムは、定位毎（例えば１度刻み）に評価値を積算して、求める。このようにして生成されるヒストグラムにおいては、本来の定位では、その定位での評価値の積算値が大きくなる。
【００５３】
次に、ＣＰＵ４０１は、ステップＳ２０１で計算した定位とその評価値を用いて、クラスタリングを実行する（ステップＳ２０３）。クラスタリングのアルゴリズムとしては、例えばｋ−ｍｅａｎｓ法を適用することができる。
【００５４】
ＣＰＵ４０１は更に、ステップＳ２０２で作成したヒストグラムと、ステップＳ２０３でのクラスタリングにより算出されたクラスタの中心とを、出力装置４０４のディスプレイ等に、ヒント情報として提示する（ステップＳ２０４）。図７は、ＣＰＵ４０１が出力装置４０４のディスプレイ等に表示する、ヒント情報の表示例を示す図である。図７において、太線の半円はヒストグラムを表し、その太線の濃淡や色などでヒストグラムの大小が表現される。また、黒丸はクラスタリングの中心を表す。このようにして、自動的に推定された定位を、ユーザに提示することができる。
【００５５】
ＣＰＵ４０１は、ステップＳ２０４でディスプレイ等に提示したヒント情報を参考に、入力装置２０３のマウス等によりユーザに定位を設定させる（ステップＳ２０５）。ユーザによる指定方法の一例としては、図６に例示される表示上で黒丸の移動、あるいは追加・削除することにより行う。その後、ＣＰＵ４０１は、図６のフローチャートで例示される図５のステップＳ１０２の定位設定処理を終了する。
【００５６】
図５の説明に戻り、ＣＰＵ４０１は、図１のラベル生成部１０３の機能に対応するラベル生成処理を実行する（ステップＳ１０３）。この処理は、ステップＳ１０２の定位設定処理で設定した定位情報から選択される、無音、任意の単一の定位、任意の２つの定位の組合せ毎にそれぞれ付与されるラベルを生成する処理である。本実施形態では、或る周波数帯及び或るフレーム番号において、同時に鳴っている音源の定位を最大で２つまでに制約する（音源が３つ以上同時に鳴っている場合も主要な２つとする）ことで、従来よりも適切な分離処理を行うものである。ラベルの数は、定位数をＮとすると、下記数１２式で示される。
【００５７】
【数１２】
0006613737.tif 000013
【００５８】
すなわち、Ｎ個の定位から０個を選択するコンビネーション（「Ｃ」はコンビネーションの演算を示す）（無音に対応）、Ｎ個の定位から１個を選択するコンビネーション（単一の定位に対応）、およびＮ個の定位から２個を選択するコンビネーション（２つの定位に対応）の総和として、ラベルの数が算出される。
【００５９】
図８は、ユーザが指定した定位の数Ｎが例えば図７に示されるように５である場合における定位の組合せとラベルの対応関係の例を示す図である。ＬＡＢＥＬ０〜１５までの１６個のラベルが生成される。
【００６０】
次に、ＣＰＵ４０１は、図１の選択部１０５の機能に対応する選択処理を実行する（ステップＳ１０４）。この処理では、フレーム毎および周波数帯毎の複素スペクトログラムのそれぞれに、ステップＳ１０３のラベル生成処理で生成されたラベルのうち最適なラベルが付与される。この処理の詳細については、後述する。
【００６１】
さらにＣＰＵ４０１は、図１の分解部１０４の機能に対応する信号分解・抽出処理を実行する（ステップＳ１０５）。この処理では、ステップＳ１０４の選択処理で付与されたラベルに基づき、時間周波数領域の信号である左右チャネルの複素スペクトログラムの信号が定位毎の信号に分解されて抽出される。
【００６２】
最後に、ＣＰＵ４０１は、図１の生成部１０６の機能に対応する生成処理を実行する（ステップＳ１０６）。この処理では、ステップＳ１０５の信号分解・抽出処理で抽出された定位毎の信号が、それぞれを時間周波数領域の信号である複素スペクトログラムから時間領域の信号に変換される。
【００６３】
上述したステップＳ１０１およびＳ１０４〜Ｓ１０６における複素スペクトログラムとラベルと定位との関係を、本実施形態のデータ構成例を示す図９を用いて説明する。図９（ａ）は、図５のステップＳ１０１の時間周波数変換処理（図１の変換部１０１）の出力である複素スペクトログラム、図９（ｂ）は、図５のステップＳ１０４の選択処理（図１の選択部１０５）の出力である最適なラベル、図９（ｃ）は図５の分解部１０４の出力である定位毎に分離した複素スペクトログラムである。
【００６４】
１フレーム分の時間領域の左チャネル信号Ｌ（ｔ）および右チャネル信号Ｒ（ｔ）に対する時間周波数変換処理によって、左チャネル（Ｌ）および右チャネル（Ｒ）毎に、図９（ａ）の１フレーム（ｆｒａｍｅ）分の複数周波数（ｆｒｅｑｕｅｎｃｙ）分、すなわち図９（ａ）の縦１列分の複素スペクトログラムが得られる。図９（ａ）の複素スペクトログラムの各四角には、時間周波数領域の複素信号が格納されているものとする。
【００６５】
本実施形態では、周波数帯毎に、図５のステップＳ１０４からＳ１０６の処理が実行される。今例えば、図９中の網掛けされた部分の任意の周波数帯に着目する。本実施形態では、図９（ａ）の網掛けされた部分の１つの周波数帯の全フレーム分の左チャネル（Ｌ）および右チャネル（Ｒ）の各複素スペクトログラムに基づいて、図５のステップＳ１０４の選択処理において、その周波数帯の全フレームに対して図９（ｂ）の網掛けされた四角部分に示されるように、図８で例示したＬＡＢＥＬ０〜１５の何れかのラベルのうちから最適なラベルが推定されて付与される。
【００６６】
そして、図５のステップＳ１０５の信号分解・抽出処理において、図９（ｂ）の網掛けされた部分の周波数帯に対して付与された全フレーム分の最適なラベルに基づいて、図９（ａ）の網掛けされた部分の１つの周波数帯の全フレーム分の左チャネル（Ｌ）および右チャネル（Ｒ）の各複素スペクトログラムが、図９（ｃ）の網掛けされた部分の１つの周波数帯における全フレーム分の定位毎（定位１〜定位５）の複素スペクトログラムに分解されて抽出される。
【００６７】
以上のステップＳ１０４とＳ１０５の処理が全ての周波数帯に対して実行されることにより、定位毎に、全フレーム分および全周波数帯分の複素スペクトログラムが得られる。その後、１つの定位に対して、フレーム毎に、全周波数帯分の複素スペクトログラムについてステップＳ１０６の生成処理が実行されることにより、各フレームの時間領域の信号が得られる。この処理が、全ての定位に対して実行されることにより、各定位の時間領域の信号が得られる。
【００６８】
図５のステップＳ１０４の選択処理の詳細について、以下に説明する。
【００６９】
まず、本実施形態では、音源の調波構造を仮定しないものとする。この仮定のもとで、ＣＰＵ４０１は、周波数帯毎に、フレーム番号順に連結した直線構造のグラフをそれぞれ作成する。以降の処理はどの周波数帯でも同様なので、或る１つの周波数帯での説明とする。前述した図９の説明では、図９（ａ）の縦方向（ｆｒｅｑｕｅｎｃｙ方向）の関係は考えずに、横方向（ｆｒａｍｅ方向）の網掛けされた部分が、上述した１つの周波数帯に対応する。図１０は、本実施形態における確率モデルの説明図である。破線枠１０００で囲まれた部分が、１つの周波数帯についてフレームを連結して作成されたグラフである。このグラフにおいて、１００１（＃ｉ）、１００１（＃ｉ＋１）などとして示されている四角で示される頂点１００１の部分が、１つの周波数帯における各フレームに対応し、１００２（＃ｉ）などとして示されている辺１００２の部分が、フレーム間の関係に対応する。
【００７０】
本実施形態では、ＣＰＵ４０１は、周波数帯毎の各フレーム、すなわち、図１０の１０００として例示されるグラフの各頂点１００１（図中１００１（＃ｉ）、１００１（＃ｉ＋１））に、それぞれ図５のステップＳ１０３で生成された全てのラベルを、各ラベル候補群１００３（図中１００３（＃ｉ）、１００３（＃ｉ＋１））として割り振る。図５のステップＳ１０３で生成されるラベルは、図８に例示したようにＬＡＢＥＬ０〜ＬＡＢＥＬ１５の１６種類がある。従って、グラフの各頂点１００１に割り振られる各ラベル候補群１００３も、それぞれ図中の０〜１５として示されるように、１６種類ずつのラベル候補を含む。
【００７１】
続いて、ＣＰＵ４０１は、グラフの頂点１００１に割り振られるラベル候補群１００３中の０〜１５の各ラベル候補に関して、各ラベル候補が含まれる頂点１００１に対応する１つの周波数帯の１つのフレームにおいて各ラベル候補に対応する最大で２つまでの定位がどの程度確からしいかを示す各ラベル確率を計算し、各ラベル候補に対応付けて図４のメモリ４０２に記憶する。
【００７２】
さらに、ＣＰＵ４０１は、グラフの辺１００２（図中１００２（＃ｉ））に対応する関係として、グラフの１つの頂点１００１（＃ｉ）に対応するラベル候補群１００３（＃ｉ）中の０〜１５のラベル候補のそれぞれが、時間が増加する方向に隣接する次の頂点１００１（＃ｉ＋１）に対応するラベル候補群１００３（＃ｉ＋１）中の０〜１５のラベル候補のそれぞれに遷移する場合の、各確率値が格納されたラベル遷移確率行列１００４（図中１００４（＃ｉ））を計算する。ラベル遷移確率行列１００４（＃ｉ）中の要素「００」には、頂点１００１（＃ｉ）のラベル候補群１００３（＃ｉ）内の０番目のラベル候補が、次の頂点１００１（＃ｉ＋１）のラベル候補群１００３（＃ｉ＋１）内の０番目のラベル候補に遷移する確率値が格納される。要素「１０」〜「Ｆ０」には、頂点１００１（＃ｉ）のラベル候補群１００３（＃ｉ）内の１〜１５番目の各ラベル候補が、次の頂点１００１（＃ｉ＋１）のラベル候補１００３（＃ｉ＋１）内の０番目のラベル候補に遷移する確率値が格納される。同様に、要素「０１」〜「Ｆ１」には、頂点１００１（＃ｉ）のラベル候補群１００３（＃ｉ）内の０〜１５番目の各ラベル候補が、次の頂点１００１（＃ｉ＋１）のラベル候補群１００３（＃ｉ＋１）内の１番目のラベル候補に遷移する確率値が格納される。以下同様にして、１６×１６のラベル遷移確率行列１００４（＃ｉ）中の各要素に、グラフの１つの頂点１００１（＃ｉ）に対応するラベル候補群１００３（＃ｉ）中の０〜１５のラベル候補のそれぞれが、次の頂点１００１（＃ｉ＋１）に対応するラベル候補群１００３（＃ｉ＋１）中の０〜１５のラベル候補のそれぞれに遷移する場合の、全ての組合せ間の各ラベル遷移確率値が計算されて格納される。
【００７３】
ＣＰＵ４０１は、グラフ１０００を構成する全ての頂点１００１のラベル候補群１００３内の各ラベル候補に関する各ラベル確率と、全ての頂点１００１間の全ての辺１００２に関するラベル遷移確率行列１００４を計算する。
【００７４】
その後、ＣＰＵ４０１は、各頂点１００１のラベル候補群１００３内の各ラベル候補に関する各ラベル確率、及び全ての辺１００２に関するラベル遷移確率行列１００４内の各ラベル遷移確率から求まる尤度を最大にするグラフの全ての頂点１００１におけるラベル候補の組合せを推定して、グラフの頂点１００１毎すなわちフレーム毎に最適なラベルを付与する。これが、前述した図９（ｂ）の網がけされた四角部分となる。
【００７５】
図１１は、以上の動作を実行する図５のステップＳ１０４の選択処理の詳細例を示すフローチャートである。この処理は、各周波数帯毎に実行されるが、以下の説明は簡単のために、１つの周波数帯の処理についてのみ説明する。
【００７６】
まず、ＣＰＵ４０１は、グラフの作成処理を実行する（ステップＳ３０１）。これにより、図１０に例示されるグラフ１０００のデータが、例えば図４のメモリ４０２上に作成される。
【００７７】
次に、ＣＰＵ４０１は、ステップＳ３０１で作成したグラフから未処理の頂点を１つずつ選択しながら（ステップＳ３０２）、ステップＳ３０８で未処理の頂点がなくなったと判定するまで、ステップＳ３０３からＳ３０７の一連の処理を頂点毎に繰り返し実行する。
【００７８】
頂点毎の繰返しにおいて、ＣＰＵ４０１はまず、ステップＳ３０２で選択した頂点に、ラベル候補群を割り振る（ステップＳ３０３）。この処理は、図１０で頂点１００１にラベル候補群１００３を割り振る動作であり、例えば図８に示される１６種類のラベルがラベル候補として割り振られる。
【００７９】
次に、ＣＰＵ４０１は、ステップＳ３０３で割り振ったラベル候補群から未処理のラベル候補を１つずつ選択しながら（ステップＳ３０４）、ステップＳ３０７で未処理のラベル候補がなくなったと判定するまで、ステップＳ３０５とＳ３０６の処理を１６個のラベル候補毎に繰り返し実行する。
【００８０】
ラベル候補毎の繰返しにおいて、ＣＰＵ４０１はまず、信号分解処理を実行する（ステップＳ３０５）。ラベル候補が単一の定位のみの場合（図８のＬＡＢＥＬ１〜ＬＡＢＥＬ５の場合）は、下記数１３式が成り立つのでこれを解いた下記数１４式にて、１つの音源の複素スペクトログラムＳＰｘを計算する。
【００８１】
【数１３】
0006613737.tif 000014
【００８２】
【数１４】
0006613737.tif 000015
【００８３】
上記数１３式および数１４式において、ＳＰＬおよびＳＰＲは、現在処理している周波数帯およびグラフの頂点に対応するフレームに対して図５のステップＳ１０１で算出されている左チャネルおよび右チャネルの各複素スペクトログラムであり、前述した数６式および数７式として計算されている。また、ＬｘおよびＲｘは、ステップＳ３０４で現在選択されているラベル候補に対応する定位を表す角度θｘに基づいて、前述した数２式および数３式（ｘをｉと置き換えたもの）により算出される値であり、図３のＤＬｉおよびＤＲｉの各逆数に対応する。
【００８４】
ラベル候補が２つの定位の場合（図８のＬＡＢＥＬ６〜ＬＡＢＥＬ１５の場合）は、下記数１５式が成り立つのでこれを解いた下記数１６式にて、２つの音源の各複素スペクトログラムＳＰｘ、ＳＰｙを計算する。
【００８５】
【数１５】
0006613737.tif 000016
【００８６】
【数１６】
0006613737.tif 000017
【００８７】
上記数１５式および数１６式において、ＳＰＬおよびＳＰＲは、数１３式および数１４式の場合と同様であり、前述した数６式および数７式として計算されている。また、Ｌｘ，ＲｘおよびＬｙ，Ｒｙはそれぞれ、ステップＳ３０４で現在選択されているラベル候補に対応する２つの定位のそれぞれを表す角度θｘ、θｙに基づき、前述の数２式および数３式（ｘ，ｙをそれぞれｉと置き換えたもの）により算出される値である。
【００８８】
ＣＰＵ４０１は、以上のようにして算出した音源の複素スペクトログラムＳＰｘまたはＳＰｘとＳＰｙを、ステップＳ３０２およびＳ３０４で現在選択しているグラフの頂点（フレームに対応）およびラベル候補に対応させて、図４のメモリ４０２に記憶する。
【００８９】
次に、ＣＰＵ４０１は、ステップＳ３０４で現在選択されているラベル候補のラベル確率を計算し、その結果をステップＳ３０２およびＳ３０４で現在選択しているグラフの頂点（フレームに対応）およびラベル候補に対応させて、図４のメモリ４０２に記憶する（ステップＳ３０６）。確率計算の一例を以下に記述する。
【００９０】
（１）無音のラベルでは、左右の振幅が共に所定の閾値より小さい場合は、無音と判定しＫ０を与え、そうでない場合は（１−Ｋ０）を与える。ここでＫ０は調整用のパラメータで、無音の誤判定を考慮した場合は１よりも若干小さい値を与える。誤判定を考慮しない場合は１を与える。
【００９１】
（２）単一の定位のラベルでは、確率ｐｉは下記数１７式のように計算する。ここでｐＳＩＬＥＮＴは（１）で求めた無音の確率、Ｋ１は調整パラメータである。数１７式は、左右の位相が一致するほど、定位がラベルの定位と一致するほど高い確率を与えるようにしている。Ｄは、数１１式の説明で前述した位相差を計算する関数である。
【００９２】
【数１７】
0006613737.tif 000018
【００９３】
（３）２つの定位からなるラベルでは、確率ｐｃは下記数１８式のように計算する。これは無音の確率ｐＳＩＬＥＮＴと単一の定位の最大確率ｍａｘ（ｐｉ）の合計の相補「１−ｐＳＩＬＥＮＴ−ｍａｘ（ｐｉ）」に相当する。なお、定位の２つの組み合わせではどれも信号分解可能なので、組み合わせによらず全て等しい確率とする。
【００９４】
【数１８】
0006613737.tif 000019
【００９５】
ステップＳ３０６の後、ＣＰＵ４０１は、他に未処理のラベル候補があるか否かを判定し（ステップＳ３０７）、この判定がＹＥＳならば、ステップＳ３０４に戻って他のラベル候補を選択し、そのラベル候補に対して上述したステップＳ３０５とＳ３０６の処理を再度実行する。
【００９６】
ステップＳ３０２で選択しているグラフの頂点に対してステップＳ３０３で割り振られたラベル候補群内の全てのラベル候補に対する処理が終了してステップＳ３０７の判定がＮＯになると、ＣＰＵ４０１は、グラフ上に他に未処理の頂点があるか否かを判定し（ステップＳ３０８）、この判定がＹＥＳならば、ステップＳ３０２に戻って他の未処理の頂点を選択し、その頂点に対して上述したステップＳ３０３からＳ３０７までの一連の処理を再度実行する。
【００９７】
グラフの全ての頂点に対する処理が終了してステップＳ３０８の判定がＮＯになると、ＣＰＵ４０１は、グラフの全ての頂点間の全ての辺（図１０の１００２に対応）に関するラベル遷移確率行列（図１０の１００４に対応）を計算する（ステップＳ３０９）。隣接フレーム間で定位や音源の振幅が急激に変化することは余りないのでこれを考慮する。以下のルールに則って確率を与える。
【００９８】
（１）定位の変化は１つまでとする。つまり２つ変化する場合の確率をゼロとする。
【００９９】
（２）無音を保持する確率を１とする。
【０１００】
（３）取って代わられる定位の振幅をＡ０、取って代わる定位の振幅をＡ１として、数１９式を計算する。ここでｌｎは対数、αは調整係数である。αは振幅の増加と減少で異なる値を設定しても良い。音楽では振幅が急激に増加することはままあるが、急激に減少することはまれであるので、これを考慮して増加時には小さめの値を、減少時には大きめの値にするとよい。なお、定位の振幅はステップＳ３０５の信号分解処理で算出された音源を表す複素スペクトログラムの振幅である。また、無音の振幅は、左右の振幅の大きい方とする。単一の定位と２つの定位の場合は、単一の定位が２つあるものとして考える。
【０１０１】
【数１９】
0006613737.tif 000020
【０１０２】
（４）単一の定位を保持する確率は、同じ定位に取って代わられるものとして計算する。
【０１０３】
（５）２つの定位を保持する確率も（４）と同様に考える。それぞれの定位を保持する確率の相乗平均を計算する。
【０１０４】
次に、ＣＰＵ４０１は、評価関数の計算処理を実行する（ステップＳ３１０）。ステップＳ３０１で作成されたグラフの各頂点に対してそれぞれ１つのラベルを与えた場合のグラフ全体の尤度（確率）は、ステップＳ３０６で算出されたラベル確率と、ステップＳ３０９で算出されたラベル遷移確率をすべて掛けあわせたものになる。ここで、グラフの頂点のラベル確率を分割し確率をグラフの辺（図１０の１００２参照）に集約した評価関数を求めれば、尤度は評価関数をすべて掛けあわせたものとなる。頂点のラベル確率の分割は、「（接続する辺の数）乗根」を計算する。本実施形態では、グラフは図１０の１０００として示されるように直線グラフであるため、グラフの両端の頂点ではその頂点に接続される辺１００２の数は１本であるからそのままの値を採用し、両端以外の頂点ではその頂点に接続される辺１００２の数は２本であるから平方根を求める。
【０１０５】
続いて、ＣＰＵ４０１は、確率伝播法のアルゴリズムを実行する（ステップＳ３１１）。グラフ全体の尤度が最大になるラベルの組み合わせが求めるべき最適なラベルである。しかし、図１０に示されるように、各フレームに対応するグラフの頂点１００１毎にラベル候補群１００３が割り振られるため、ラベル遷移確率行列１００４を見てわかるように、頂点間のラベル候補の組合せは膨大となるため、直接解を求めることが事実上不可能である。このため、本実施形態では、よく知られている確率伝播法のアルゴリズムにより、各頂点のラベル候補の周辺確率を計算することで解を求めるようにする。
【０１０６】
確率伝播法のアルゴリズムとしては、例えば文献「確率モデルによる画像処理技術入門」（田中和之著、森北出版株式会社）に掲載されている手法が採用できる。
【０１０７】
図１２は、確率伝播法による処理の例を示すフローチャートである。確率伝播法では、グラフの辺の頂点の一方から他方へと向かうメッセージと呼ばれる値を求める。そして、ある頂点の周辺確率はその頂点に向かうメッセージを掛け合わせることで計算することができる。メッセージの総数はグラフの辺の数の倍あり、各メッセージの要素数は向かう頂点の状態数に等しい。
【０１０８】
確率伝播法による処理において、ＣＰＵ４０１はまず、メッセージの初期化を行う（ステップＳ４０１）。例えば、各メッセージの要素に状態数の逆数を設定する。
【０１０９】
次に、ＣＰＵ４０１は、更新前のメッセージと図１１のステップＳ３１０で計算した評価関数の行列を使ってメッセージの更新を行う（ステップＳ４０２）。各々のメッセージの更新は独立しているので、処理の並列化が容易であり、処理の高速化がしやすい。これは確率伝播法の利点である。
【０１１０】
ＣＰＵ４０１は、メッセージが収束したかの判定（ステップＳ４０３）と、所定回数メッセージを更新したかの判定（ステップＳ４０４）により、メッセージ更新の終了の判定を行う。一般的に、メッセージの更新を繰り返すと収束するが、極まれに一部のメッセージが収束しない場合があるので更新の上限回数を設けている。
【０１１１】
メッセージ更新が終了しない場合（ステップＳ４０３の判定がＮＯ、かつステップＳ４０４の判定がＮＯ）、ステップＳ４０２のメッセージ更新を繰り返す。
【０１１２】
メッセージが収束した場合（ステップＳ４０３の判定がＹＥＳ）や所定回数メッセージを更新した場合（ステップＳ４０４の判定がＹＥＳ）は、メッセージの更新を終了して、各頂点の周辺確率、つまりラベル候補の確率計算を行う（ステップＳ４０５）。
【０１１３】
その後、ＣＰＵ４０１は、図１２のフローチャートで例示される図１１のステップＳ３１１の確率伝搬法による処理を終了する。
【０１１４】
最後に、ＣＰＵ４０１は、ステップＳ３１１で算出した各頂点のラベル候補の周辺確率が最も高いものを最適なラベルとして決定する（ステップＳ３１２）。
【０１１５】
以上の図１１のフローチャートで例示される図５のステップＳ１０４の選択処理が終了すると、前述したように、ＣＰＵ４０１は、図５のステップＳ１０５で、信号分解・抽出処理を実行する。この処理では、ステップＳ１０４の選択処理で付与された最適なラベルに基づき左右チャネルの複素スペクトログラムの信号が定位毎の信号に分解されて抽出される。この処理は、前述した図１１のステップＳ３０５の処理と同じであるため、選択されたラベルに対応してステップＳ３０５で既に計算され図４のメモリ４０２に保持されていた音源の複素スペクトログラムＳＰｘ（単一の定位の場合）またはＳＰｘとＳＰｙ（２つの定位の場合）を読み出せばよい。
【０１１６】
このようにして読み出された最適なラベルに対応する音源の複素スペクトログラムＳＰｘ（単一の定位の場合）またはＳＰｘとＳＰｙ（２つの定位の場合）が、それぞれに対応する周波数帯およびフレームに配置され、前述したように、定位毎に、その複素スペクトログラムに対して生成処理が実行され（図５のステップＳ１０６）、定位毎に分解された波形信号が生成される。なお、定位毎に、その定位に対応する複素スペクトログラムが割り当てられていない周波数帯およびフレームについては、無音に対応する信号値をパディングして処理する。
【０１１７】
以上説明したように、本実施形態では、最大で２つまで音源（定位）の適切な組み合わせを推測し、これに基づき音源の分離を実行できるようにした。これにより、従来技術で課題となっていた、２つの音源とは異なる音源と推測され正しく音源が分離されない問題が大幅に軽減され、より適切な音源信号変換を実現することが可能となる。
【０１１８】
以上の実施形態に関して、更に以下の付記を開示する。
（付記１）
複数種の定位の組み合わせを入力する入力部と、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換する変換処理と、前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記入力された定位の組み合わせの中から最適な定位の組み合わせを選択する選択処理と、前記選択された最適な定位の組み合わせに基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記設定された複数種の定位毎でかつ周波数帯域毎の信号に分解する分解処理と、を実行する処理部と、
を備える音源信号変換装置。
（付記２）
前記処理部は、
前記設定された複数種の定位及び前記定位の組み合わせそれぞれに対して付与されるラベルを生成するラベル生成処理と、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記定位毎に生成する生成処理と、
をさらに実行し、
前記選択処理は、前記周波数領域毎に、
前記フレームのそれぞれに、前記ラベル生成処理にて生成される全てのラベルをラベル候補として割り振る処理と、
前記割り振られたラベル候補毎に、前記ラベル候補に対応する周波数領域毎のフレームに属する周波数の信号が、前記ラベル候補に対応する前記定位及び前記定位の組み合わせの夫々がどの程度確かかを示すラベル確率を計算する処理と、
前記時間軸方向に隣接するフレームに割り振られたラベル候補の全ての組合せ間のラベル遷移確率を計算する処理と、
前記ラベル確率及び前記ラベル遷移確率から求まる尤度を最大にするラベル候補の組合せを推定して前記フレーム毎に最適なラベルを付与する処理と、
を実行する、付記１に記載の音源信号変換装置。
（付記３）
前記処理部は、
前記定位が第１の既定時間内に第１の既定数以上変化しないように、及び前記分解処理での分解結果の振幅が第２の既定時間内に第２の既定値以上変化しないように、前記ラベル遷移確率を計算する処理を実行する、付記２に記載の音源信号変換装置。
（付記４）
前記処理部は、
同時に変化する前記定位が１つまでに制限されるように、前記ラベル遷移確率を計算する処理を実行する、付記２又は３に記載の音源信号変換装置。
（付記５）
前記処理部は、
前記ラベル確率及び前記ラベル遷移確率から求まる尤度を最大にする前記ラベル候補の組合せを、確立伝搬法のアルゴリズムに従って推定する処理を実行する、付記２ないし４のいずれかに記載の音源信号変換装置。
（付記６）
前記設定部は、外部より複数種の定位及び前記複数種の定位の組み合わせを設定する付記１に記載の音源信号変換装置。
（付記７）
前記処理部はさらに、
前記変換処理にて得られた左右各チャネルそれぞれの複数種の周波数領域毎の信号に対して、周波数帯毎および前記フレーム毎に、定位及び前記定位が単一か否かを表す評価値を算出する処理と、
前記算出された定位とその評価値とに基づき、前記定位毎に評価値を積算することによりヒストグラムを作成する処理と、
前記作成されたヒストグラムをディスプレイに提示する処理と、
を実行する付記６に記載の音源信号変換装置。
（付記８）
複数種の定位の組み合わせを入力する入力部と、処理部とを有する音源信号変換装置に用いられる音源信号変換方法であって、前記処理部が、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換し、
前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記設定された定位の組み合わせの中から最適な定位の組み合わせを選択し、
前記選択された最適な定位の組み合わせに基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記設定された複数種の定位毎でかつ周波数帯域毎の信号に分解し、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記定位毎に生成する、音源信号変換方法。
（付記９）
複数種の定位の組み合わせを入力する入力部と、処理部とを有する音源信号変換装置として用いられるコンピュータに、
複数種の周波数を含むステレオ信号の左右チャネル信号それぞれを、所定時間間隔で分割されたフレーム毎でかつ、複数種の周波数帯域毎の信号に変換するステップと、
前記変換された左右各チャネルそれぞれの複数種の周波数帯域毎の信号に対して、周波数帯域毎および前記フレーム毎に、前記設定された定位の組み合わせの中から最適な定位の組み合わせを選択するステップと、
前記選択された最適な定位の組み合わせに基づいて、前記左右チャネルそれぞれの周波数帯域毎の信号を前記フレーム毎に、前記設定された複数種の定位毎でかつ周波数帯域毎の信号に分解するステップと、
前記分解された周波数帯域毎の信号夫々から、前記周波数帯域毎の信号を含む信号を前記定位毎に生成するステップと、
を実行させるプログラム。
【符号の説明】
【０１１９】
１００  音源信号変換装置
１０１  変換部
１０２  定位設定部
１０３  ラベル生成部
１０４  分解部
１０５  選択部
１０６  生成部
４０１  ＣＰＵ
４０２  メモリ
４０３  入力装置
４０４  出力装置
４０５  外部記憶装置
４０６  可搬記録媒体駆動装置
４０７  通信インターフェース
４０８  バス
４０９  可搬記録媒体
１０００  グラフ
１００１  頂点
１００２  辺
１００３  ラベル候補群
１００４  ラベル遷移確率行列
【図１】
0006613737.tif 000021
【図２】
0006613737.tif 000022
【図３】
0006613737.tif 000023
【図４】
0006613737.tif 000024
【図５】
0006613737.tif 000025
【図６】
0006613737.tif 000026
【図７】
0006613737.tif 000027
【図８】
0006613737.tif 000028
【図９】
0006613737.tif 000029
【図１０】
0006613737.tif 000030
【図１１】
0006613737.tif 000031
【図１２】
0006613737.tif 000032