JP6981330B2 - 頭外定位処理装置、頭外定位処理方法、及びプログラム - Google Patents

Description

本発明は、頭外定位処理装置、頭外定位処理方法、及びプログラムに関する。

例えば、音像定位技術として、ヘッドホンを用いて受聴者の頭部の外側に音像を定位させる頭外定位技術がある。頭外定位技術では、ヘッドホンから耳（鼓膜）までの特性（外耳道伝達特性）をキャンセルし、ステレオスピーカから耳までの４本の特性（空間音響伝達特性）を与えることにより、音像を頭外に定位させている。

頭外定位再生においては、ヘッドホンから発した測定信号（インパルス音等）を聴取者本人の耳に設置したマイクロフォン（以下、マイクとする）で録音する（特許文献１）。インパルス応答で得られた収音信号に基づいて、処理装置がヘッドホン特性を測定する。処理装置がヘッドホン特性に対する逆フィルタを作成する。処理装置が、空間音響伝達特性を畳み込んだ後、逆フィルタをさらに畳み込むことにより、頭外定位再生を実現することができる。

特開２０１４−３８６０８号公報

しかしながら、ヘッドホンの装着状態に応じて、ヘッドホンのスピーカユニットから鼓膜までの外耳道伝達特性（外耳道伝達関数ともいう）が変化してしまうことがある。例えば、ヘッドホンを装着する都度、ヘッドホンの装着位置がずれるおそれがある。この場合、ヘッドホンのスピーカユニットの位置がずれるため、外耳道伝達特性が変化してしまう。つまり、耳にマイクを装着して、外耳道伝達特性を測定する測定時と、耳からマイクを外して、頭外定位受聴を行う受聴時とで、ヘッドホンの装着位置が変わってしまう。この場合、外耳道伝達特性が変化し、適切に頭外定位処理を行うことができなくなるおそれがある。

本実施形態は上記の点に鑑みなされたもので、適切に頭外定位処理を行うことができる頭外定位処理装置、頭外定位処理方法、及びプログラムを提供することを目的とする。

本実施形態にかかる頭外定位処理装置は、参照用マイクを有する音声出力部と、前記参照用マイクが収音した収音信号に基づく参照信号を取得する参照信号取得部と、前記音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記参照信号取得部が取得した第１の参照信号を記憶する第１の記憶部と、前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照信号取得部が取得した第２の参照信号を記憶する第２の記憶部と、前記測定状態において、前記測定用マイクが収音した収音信号に基づいて取得された伝達特性を記憶する第３の記憶部と、前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出する変換関数算出部と、前記変換関数を用いて、前記伝達特性の周波数特性を補正する補正部と、前記補正された周波数特性に基づいて、フィルタを生成するフィルタ生成部と、前記フィルタを用いて、頭外定位処理を行う頭外定位処理部と、を備えたものである。

本実施形態にかかる頭外定位処理方法は、参照用マイクを有する音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記測定用マイクで収音された収音信号に基づく伝達特性を取得するステップと、前記測定状態において、前記参照用マイクで収音された収音信号に基づく第１の参照信号を取得するステップと、前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照用マイクで収音された収音信号に基づく第２の参照信号を取得するステップと、前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出するステップと、前記変換関数を用いて、前記伝達特性の周波数特性を補正するステップと、前記補正された周波数特性に基づいて、フィルタを生成するステップと、前記フィルタを用いて、頭外定位処理を行うステップと、を備えたものである。

本実施形態にかかるプログラムは、コンピュータに、参照用マイクを有する音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記測定用マイクで収音された収音信号に基づく伝達特性を取得するステップと、前記測定状態において、前記参照用マイクで収音された収音信号に基づく第１の参照信号を取得するステップと、前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照用マイクで収音された収音信号に基づく第２の参照信号を取得するステップと、前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出するステップと、前記変換関数を用いて、前記伝達特性の周波数特性を補正するステップと、前記補正された周波数特性に基づいて、フィルタを生成するステップと、前記フィルタを用いて、頭外定位処理を行うステップと、を実行させるものである。

本実施形態によれば、適切に頭外定位処理を行うことができる頭外定位処理装置、頭外定位処理方法、及びプログラムを提供することができる。

本実施の形態に係る頭外定位処理装置を示すブロック図である。 外耳道伝達特性の測定構成を示す図である。 外耳道伝達特性を補正するための構成を示す図である。 外耳道伝達特性を補正するための処理を示すフローチャートである。 変換関数を算出するための処理を示すフローチャートである。 パワースペクトルの極値間のベクトルを示す図である。 ベクトルの変換手法を示す図である。 参照信号と外耳道伝達特性の対数パワースペクトルを示す図である。

フィルタを用いた音像定位処理の概要について説明する。本実施形態にかかる頭外定位処理は、空間音響伝達特性と外耳道伝達特性を用いて頭外定位処理を行うものである。空間音響伝達特性は、スピーカなどの音源から外耳道までの伝達特性である。外耳道伝達特性は、ヘッドホン又はイヤホンのスピーカユニットから鼓膜までの伝達特性であり、ヘッドホン特性ともいう。本実施形態では、ヘッドホン又はイヤホンを装着した状態での外耳道伝達特性を測定し、それらの測定データを用いてフィルタを生成する処理に特徴を有するものである。

本実施の形態にかかる頭外定位処理は、パーソナルコンピュータ、スマートホン、タブレットＰＣなどのユーザ端末で実行される。ユーザ端末は、プロセッサ等の処理手段、メモリやハードディスクなどの記憶手段、液晶モニタ等の表示手段、タッチパネル、ボタン、キーボード、マウスなどの入力手段を有する情報処理装置である。ユーザ端末は、データを送受信する通信機能を有していてもよい。さらに、ユーザ端末には、左右の出力ユニットを有する音声出力部（ヘッドホン又はイヤホン）が接続される。以下、音声出力部をヘッドホンとした場合の構成について例示する。

（頭外定位処理装置）
本実施の形態にかかる音場再生装置の一例である頭外定位処理装置１００を図１に示す。図１は、頭外定位処理装置１００のブロック図である。頭外定位処理装置１００は、ヘッドホン４３を装着するユーザＵに対して音場を再生する。そのため、頭外定位処理装置１００は、ＬｃｈとＲｃｈのステレオ入力信号ＸＬ、ＸＲについて、音像定位処理を行う。ＬｃｈとＲｃｈのステレオ入力信号ＸＬ、ＸＲは、ＣＤ（Compact Disc）プレイヤーなどから出力されるアナログのオーディオ再生信号、又は、mp3(MPEG Audio Layer-3)等のデジタルオーディオデータである。なお、オーディオ再生信号、又はデジタルオーディオデータをまとめて再生信号と称する。すなわち、ＬｃｈとＲｃｈのステレオ入力信号ＸＬ、ＸＲが再生信号となっている。

なお、頭外定位処理装置１００は、物理的に単一な装置に限られるものではなく、一部の処理が異なる装置で行われてもよい。例えば、一部の処理がパソコンなどにより行われ、残りの処理がヘッドホン４３に内蔵されたＤＳＰ(Digital Signal Processor)などにより行われてもよい。

頭外定位処理装置１００は、頭外定位処理部１０、フィルタ部４１、フィルタ部４２、及びヘッドホン４３を備えている。頭外定位処理部１０、フィルタ部４１、及びフィルタ部４２は、具体的にはプロセッサ等により実現可能である。

頭外定位処理部１０は、畳み込み演算部１１〜１２、２１〜２２、及び加算器２４、２５を備えている。畳み込み演算部１１〜１２、２１〜２２は、空間音響伝達特性を用いた畳み込み処理を行う。頭外定位処理部１０には、ＣＤプレイヤーなどからのステレオ入力信号ＸＬ、ＸＲが入力される。頭外定位処理部１０には、空間音響伝達特性が設定されている。頭外定位処理部１０は、各ｃｈのステレオ入力信号ＸＬ、ＸＲに対し、空間音響伝達特性のフィルタ（以下、空間音響フィルタとも称する）を畳み込む。空間音響伝達特性は被測定者の頭部や耳介で測定した頭部伝達関数ＨＲＴＦでもよいし、ダミーヘッドまたは第三者の頭部伝達関数であってもよい。

４つの空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓを１セットとしたものを空間音響伝達関数とする。畳み込み演算部１１、１２、２１、２２で畳み込みに用いられるデータが空間音響フィルタとなる。空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓを所定のフィルタ長で切り出すことで、空間音響フィルタが生成される。

空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓのそれぞれは、インパルス応答測定などにより、事前に取得されている。例えば、ユーザＵが左右の耳にマイクをそれぞれ装着する。ユーザＵの前方に配置された左右のスピーカが、インパルス応答測定を行うための、インパルス音をそれぞれ出力する。そして、スピーカから出力されたインパルス音等の測定信号をマイクで収音する。マイクでの収音信号に基づいて、空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓが取得される。左スピーカと左マイクとの間の空間音響伝達特性Ｈｌｓ、左スピーカと右マイクとの間の空間音響伝達特性Ｈｌｏ、右スピーカと左マイクとの間の空間音響伝達特性Ｈｒｏ、右スピーカと右マイクとの間の空間音響伝達特性Ｈｒｓが測定される。

そして、畳み込み演算部１１は、Ｌｃｈのステレオ入力信号ＸＬに対して空間音響伝達特性Ｈｌｓに応じた空間音響フィルタを畳み込む。畳み込み演算部１１は、畳み込み演算データを加算器２４に出力する。畳み込み演算部２１は、Ｒｃｈのステレオ入力信号ＸＲに対して空間音響伝達特性Ｈｒｏに応じた空間音響フィルタを畳み込む。畳み込み演算部２１は、畳み込み演算データを加算器２４に出力する。加算器２４は２つの畳み込み演算データを加算して、フィルタ部４１に出力する。

畳み込み演算部１２は、Ｌｃｈのステレオ入力信号ＸＬに対して空間音響伝達特性Ｈｌｏに応じた空間音響フィルタを畳み込む。畳み込み演算部１２は、畳み込み演算データを、加算器２５に出力する。畳み込み演算部２２は、Ｒｃｈのステレオ入力信号ＸＲに対して空間音響伝達特性Ｈｒｓに応じた空間音響フィルタを畳み込む。畳み込み演算部２２は、畳み込み演算データを、加算器２５に出力する。加算器２５は２つの畳み込み演算データを加算して、フィルタ部４２に出力する。

フィルタ部４１、４２にはヘッドホン特性（ヘッドホンのスピーカユニットとマイク間の特性）をキャンセルする逆フィルタが設定されている。そして、頭外定位処理部１０での処理が施された再生信号（畳み込み演算信号）に逆フィルタを畳み込む。フィルタ部４１で加算器２４からのＬｃｈ信号に対して、Ｌｃｈ側のヘッドホン特性の逆フィルタを畳み込む。同様に、フィルタ部４２は加算器２５からのＲｃｈ信号に対して、Ｒｃｈ側のヘッドホン特性の逆フィルタを畳み込む。逆フィルタは、ヘッドホン４３を装着した場合に、ヘッドホンユニットからマイクまでのヘッドホン特性をキャンセルする。マイクは、外耳道入口から鼓膜までの間ならばどこに配置してもよい。逆フィルタは、後述するように、ユーザＵ本人の特性の測定結果から算出されている。

フィルタ部４１は、処理されたＬｃｈ信号ＹＬをヘッドホン４３の左ユニット４３Ｌに出力する。フィルタ部４２は、処理されたＲｃｈ信号ＹＲをヘッドホン４３の右ユニット４３Ｒに出力する。ユーザＵは、ヘッドホン４３を装着している。ヘッドホン４３は、Ｌｃｈ信号ＹＬとＲｃｈ信号ＹＲ（以下、Ｌｃｈ信号ＹＬとＲｃｈ信号をまとめてステレオ信号ともいう）をユーザＵに向けて出力する。これにより、ユーザＵの頭外に定位された音像を再生することができる。

このように、頭外定位処理装置１００は、空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに応じた空間音響フィルタと、ヘッドホン特性の逆フィルタを用いて、頭外定位処理を行っている。以下の説明において、空間音響伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに応じた空間音響フィルタと、ヘッドホン特性の逆フィルタとをまとめて頭外定位処理フィルタとする。２ｃｈのステレオ再生信号の場合、頭外定位フィルタは、４つの空間音響フィルタと、２つの逆フィルタとから構成されている。そして、頭外定位処理装置１００は、ステレオ再生信号に対して合計６個の頭外定位フィルタを用いて畳み込み演算処理を行うことで、頭外定位処理を実行する。頭外定位フィルタは、ユーザＵ個人の測定に基づくものであることが好ましい。例えば，ユーザＵの耳に装着されたマイクが収音した収音信号に基づいて、頭外定位フィルタが設定されている。

このように空間音響フィルタと、ヘッドホン特性の逆フィルタはオーディオ信号用のフィルタである。これらのフィルタが再生信号（ステレオ入力信号ＸＬ、ＸＲ）に畳み込まれることで、頭外定位処理装置１００が、頭外定位処理を実行する。

（外耳道伝達特性の測定装置）
次に、逆フィルタを生成するために、外耳道伝達特性を測定する測定装置２００について、図２を用いて説明する。図２は、被測定者１に対して外耳道伝達特性を測定するための構成を示している。測定装置２００は、マイクユニット２と、ヘッドホン４３と、処理装置２０１と、を備えている。なお、ここでは、被測定者１は、図１のユーザＵと同一人物となっている。

処理装置２０１には、マイクユニット２と、ヘッドホン４３と、が接続されている。なお、マイクユニット２は、ヘッドホン４３に脱着可能に取り付けられていてもよい。マイクユニット２は、左マイク２Ｌと、右マイク２Ｒとを備えている。左マイク２Ｌは、被測定者１の左耳９Ｌに装着される。右マイク２Ｒは、被測定者１の右耳９Ｒに装着される。処理装置２０１は、頭外定位処理装置１００と同じ処理装置であってもよく、異なる処理装置であってよい。以下の説明では、頭外定位処理装置１００と処理装置２０１が同じ装置であるとして説明する。

ヘッドホン４３は、ヘッドホンバンド４３Ｂと、左ユニット４３Ｌと、右ユニット４３Ｒとを、有している。左ユニット４３Ｌと、右ユニット４３Ｒとはそれぞれ、左右の耳９Ｌ、９Ｒに対して音を出力する出力ユニットである。ヘッドホンバンド４３Ｂは、左ユニット４３Ｌと右ユニット４３Ｒとを連結する。左ユニット４３Ｌは被測定者１の左耳９Ｌに向かって音を出力する。右ユニット４３Ｒは被測定者１の右耳９Ｒに向かって音を出力する。ヘッドホン４３は密閉型、開放型、半開放型、または半密閉型等であり、ヘッドホンの種類を問わない。マイクユニット２が被測定者１に装着された状態で、被測定者１がヘッドホン４３を装着する。すなわち、左マイク２Ｌ、右マイク２Ｒが装着された左耳９Ｌ、右耳９Ｒにヘッドホン４３の左ユニット４３Ｌ、右ユニット４３Ｒがそれぞれ装着される。ヘッドホンバンド４３Ｂは、左ユニット４３Ｌと右ユニット４３Ｒとをそれぞれ左耳９Ｌ、右耳９Ｒに押し付ける付勢力を発生する。

左マイク２Ｌは、ヘッドホン４３の左ユニット４３Ｌから出力された音を収音する。右マイク２Ｒは、ヘッドホン４３の右ユニット４３Ｒから出力された音を収音する。左マイク２Ｌ、及び右マイク２Ｒのマイク部は、外耳孔近傍の収音位置に配置される。左マイク２Ｌ、及び右マイク２Ｒは、ヘッドホン４３に干渉しないように構成されている。すなわち、左マイク２Ｌ、及び右マイク２Ｒは左耳９Ｌ、右耳９Ｒの適切な位置に配置された状態で、被測定者１がヘッドホン４３を装着することができる。

処理装置２０１は、ヘッドホン４３に対して測定信号を出力する。これにより、ヘッドホン４３はインパルス音などを発生する。具体的には、左ユニット４３Ｌから出力されたインパルス音を左マイク２Ｌで測定する。右ユニット４３Ｒから出力されたインパルス音を右マイク２Ｒで測定する。測定信号の出力時に、マイク２Ｌ、２Ｒが収音信号を取得することで、インパルス応答測定が実施される。

処理装置２０１は、インパルス応答測定に基づく収音信号をメモリなどに記憶する。これにより、左ユニット４３Ｌと左マイク２Ｌとの間の伝達特性（すなわち、左耳の外耳道伝達特性）と、右ユニット４３Ｒと右マイク２Ｒとの間の伝達特性（すなわち、右耳の外耳道伝達特性）が取得される。左マイク２Ｌで取得された左耳の外耳道伝達特性をＬｃｈ（左ｃｈ）の外耳道伝達特性とし、右マイク２Ｒで取得された右耳の外耳道伝達特性をＲｃｈ（右ｃｈ）の外耳道伝達特性とする。処理装置２０１が伝達特性の測定データを所定のフィルタ長で切り出すことで、フィルタ係数が求められる。処理装置２０１は、フィルタ係数から外耳道伝達特性（ヘッドホン特性）を打ち消すような逆フィルタを算出する。

処理装置２０１は、伝達特性の測定データをそれぞれ記憶するメモリなどを有している。なお、処理装置２０１は、外耳道伝達特性を測定するための測定信号として、インパルス信号やＴＳＰ（Ｔｉｍｅ Ｓｔｒｅｔｃｈｅｄ Ｐｕｌｓｅ）信号等を発生する。測定信号はインパルス音等の測定音を含んでいる。

（外耳道伝達特性の補正）
本実施の形態にかかる頭外定位処理装置１００は、ヘッドホン４３の装着状態に応じて、外耳道伝達特性を補正する。そして、頭外定位処理装置１００が、補正後の外耳道伝達特性に基づいて、逆フィルタを算出している。このようにすることで、頭外定位処理装置１００が、ヘッドホン装着状態に適応した逆フィルタを用いて、頭外定位処理することができる。

図３を用いて、外耳道伝達特性の補正処理について説明する。図３は、外耳道伝達特性を補正するための構成を示す図である。図３では、説明の簡略化のため、左ユニット４３Ｌのみを示しているが、右ユニット４３Ｒについても同様の構成となっている。従って、右ユニット４３Ｒに関する説明を適宜省略する。

頭外定位処理装置１００は、測定信号生成部１１１と、Ｄ／Ａコンバータ１１２と、Ａ／Ｄコンバータ１２１と、参照信号取得部１２２と、Ａ／Ｄコンバータ１３１と、ＥＣＴＦ取得部１３２と、メモリ１４０と、変換関数算出部１５１と、補正部１５２と、逆フィルタ生成部１５３と、を備えている。メモリ１４０は、第１の記憶部１４１と、第２の記憶部１４２と、第３の記憶部１４３と、を備えている。

左ユニット４３Ｌは、ハウジング４５と、ヘッドホンスピーカ４６と、内蔵マイク４８とを備えている。ハウジング４５は、耳を覆うイヤーカップを備えている。ハウジング４５には、ヘッドホンスピーカ４６と、内蔵マイク４８とが設けられている。ヘッドホンスピーカ４６は、磁器回路や振動板等を有しており、ユーザＵの左耳に対して音を出力する。

内蔵マイク４８は、左ユニット４３Ｌに内蔵されている。内蔵マイク４８は、ハウジング４５で覆われた空間に配置されている。内蔵マイク４８は、ヘッドホンスピーカ４６から出力された音を収音する。内蔵マイク４８で収音した信号を参照用収音信号とする。つまり、内蔵マイク４８は、参照用収音信号を収音するための参照用マイクである。内蔵マイク４８と、ヘッドホンスピーカ４６は、ハウジング４５に固定されている。よって、ヘッドホンスピーカ４６に対する内蔵マイク４８の位置は一定となる。

なお、ヘッドホン４３はＢｌｕｅｔｏｏｔｈ（登録商標）等の無線通信を用いたワイヤレスタイプであってもよい。さらに、ヘッドホン４３は、一部の処理を実施するＤＳＰを備えていてもよい。ヘッドホン４３のＤＳＰ等が図３に示すブロックの一部、又は全ての処理を行ってもよい。例えば、Ｄ／Ａコンバータ１１２、Ａ／Ｄコンバータ１２１、Ａ／Ｄコンバータ１３１等がヘッドホン４３に内蔵されていてもよい。

さらに、左耳９Ｌ（図３では不図示）には、外耳道伝達特性を測定するための左マイク２Ｌが配置されている。左マイク２Ｌは、図２に示したように、外耳道伝達特性を測定するための測定用マイクである。左マイク２Ｌは左耳９Ｌに対して脱着可能に設けられている。外耳道伝達特性の測定時には、左耳９Ｌに左マイク２Ｌが装着される。また、ヘッドホン４３により音楽を頭外定位受聴する時（以下、単に受聴時とする）、左マイク２Ｌは、左耳９Ｌから取り外される。なお、左マイク２Ｌは、左ユニット４３Ｌに脱着可能に取り付けられていてもよく、左ユニット４３Ｌとは独立した構成となっていてもよい。例えば、左マイク２Ｌは、左耳９Ｌに直接取り付けられていてもよい。

外耳道伝達特性の測定時には、左耳９Ｌに左マイク２Ｌを装着した装着状態となり、この状態を測定状態とする。ヘッドホン４３により音楽を頭外定位受聴する受聴時には、左耳９Ｌから左マイク２Ｌを取り外した状態を非装着状態となり、この状態を受聴状態とする。

測定信号生成部１１１は、測定信号を生成する。測定信号生成部１１１で生成された測定信号は、Ｄ／Ａコンバータ１１２でＤ／Ａ変換されて、ヘッドホンスピーカ４６に出力される。ヘッドホンスピーカ４６が伝達特性を測定するための測定信号を出力する。

測定状態において、マイク２Ｌは、ヘッドホンスピーカ４６からの測定信号を収音する。マイク２Ｌで収音された収音信号（特性用収音信号ともいう）は、Ａ／Ｄコンバータ１３１でＡ／Ｄ変換される。Ａ／Ｄ変換された特性用収音信号は、ＥＣＴＦ取得部１３２に出力される。なお、インパルス応答測定を複数回行って、特性用収音信号を同期加算してもよい。

ＥＣＴＦ取得部１３２は、特性用収音信号に基づいて、外耳道伝達特性（ＥＣＴＦ）を取得する。例えば、ＥＣＴＦ取得部１３２は、ＦＦＴ（高速フーリエ変換）により、時間領域の特性用収音信号から周波数領域の外耳道伝達特性を算出する。これにより、外耳道伝達特性のパワー特性（パワースペクトル）と、位相特性（位相スペクトル）が生成される。なお、パワースペクトルの代わりに振幅スペクトルを生成してもよい。なお、ＥＣＴＦ取得部１３２は、離散フーリエ変換や離散コサイン変換等により、特性用収音信号を周波数領域のデータ（周波数特性）に変換することができる。

測定状態、及び、受聴状態の両方で、内蔵マイク４８は、ヘッドホンスピーカ４６からの測定信号を収音する。内蔵マイク４８で収音された収音信号（参照用収音信号ともいう）は、Ａ／Ｄコンバータ１２１でＡ／Ｄ変換される。Ａ／Ｄ変換された参照用収音信号は、参照信号取得部１２２に出力される。具体的には、受聴者Ｕが左マイク２Ｌを装着した測定状態と、装着していない受聴状態との両方で、インパルス応答測定が実施される。なお、インパルス応答測定を複数回行って、参照用収音信号を同期加算してもよい。

参照信号取得部１２２は、参照用収音信号に基づいて、参照信号を取得する。例えば、参照信号取得部１２２は、ＦＦＴ（高速フーリエ変換）により、時間領域の参照用収音信号から周波数領域の参照信号を算出する。これにより、参照用収音信号のパワー特性（パワースペクトル）と、位相特性（位相スペクトル）が生成される。なお、パワースペクトルの代わりに振幅スペクトルを生成してもよい。なお、参照信号取得部１２２は、離散フーリエ変換や離散コサイン変換等により、参照用収音信号を周波数領域のデータ（周波数特性）に変換することができる。

ここで、測定状態において、内蔵マイク４８で収音された参照用収音信号を第１の参照用収音信号とする。第１の参照用収音信号に基づいて取得された参照信号を第１の参照信号とする。第１の参照用収音信号は、特性用収音信号とは、実質的に同時に収音される。つまり、同じインパルス応答測定での測定された第１の参照用収音信号は、特性用収音信号に基づいて、第１の参照信号と外耳道伝達特性とが取得される。

受聴状態において、内蔵マイク４８で収音された参照用収音信号を第２の参照用収音信号とする。第２の参照用収音信号に基づいて取得された参照信号を第２の参照信号とする。第１の参照信号と第２の参照信号は、ヘッドホン４３の装着状態に応じて変化する信号である。換言すると、第１の参照信号と第２の参照信号との差が、ヘッドホン４３の装着状態の違いに相当する。

メモリ１４０は、第１の記憶部１４１と、第２の記憶部１４２と、第３の記憶部１４３とを備えている。第１の記憶部１４１は、第１の参照信号を記憶する。第２の記憶部１４２は、第２の参照信号を記憶する。第３の記憶部１４３は、外耳道伝達特性を記憶する。第１の記憶部１４１、第２の記憶部１４２、及び第３の記憶部１４３は、それぞれメモリ１４０において静的に確保しておいてもよく、任意の領域を動的に確保してもよい。

具体的には、第３の記憶部１４３は、外耳道伝達特性のパワースペクトル、及び位相スペクトルを記憶する。第１の記憶部１４１は、第１の参照信号のパワースペクトルを記憶する。第２の記憶部１４２は、第２の参照信号のパワースペクトルを記憶する。第１の記憶部１４１、及び第２の記憶部１４２は、それぞれ参照信号の位相スペクトルを記憶していてもよく、記憶していなくてもよい。

第１の記憶部１４１と、第２の記憶部１４２と、第３の記憶部１４３とは、物理的に単一なメモリ装置であってもよく、異なる装置であってもよい。例えば、１つのメモリ１４０が、外耳道伝達特性、第１の参照信号、及び第２の参照信号を記憶していてもよい。あるは、２つ以上のメモリに分けて、外耳道伝達特性、第１の参照信号、及び第２の参照信号が記憶されていてもよい。

変換関数算出部１５１は、第１の記憶部１４１、及び第２の記憶部１４２から第１の参照信号と第２の参照信号を読み出す。そして、変換関数算出部１５１は、第１の参照信号と第２の参照信号とに基づいて、周波数特性の変換関数を算出する。変換関数算出部１５１は、２つのパワースペクトルに基づいて、変換関数を算出している。つまり、変換関数算出部１５１は、２つの参照信号のパワースペクトルを比較することで、変換関数を算出している。

補正部１５２は、第３の記憶部１４３から外耳道伝達特性を読み出す。そして、補正部１５２は、変換関数を用いて、外耳道伝達特性を補正する。補正部１５２は、外耳道伝達特性のパワースペクトルを補正している。ここでは、補正部１５２は、パワースペクトルのみを補正しており、位相スペクトルについては補正していないが、位相スペクトルについて補正してもよい。

逆フィルタ生成部１５３は、補正された外耳道伝達特性（パワースペクトル）を用いて、逆フィルタを算出する。具体的には、逆フィルタ生成部１５３は、逆離散フーリエ変換により、補正後のパワースペクトル（振幅特性）と位相スペクトル（位相特性）とを用いて時間信号を算出する。逆フィルタ生成部１５３は、時間信号に基づいて、所定のフィルタ長の逆フィルタを算出する。

以下、逆フィルタを生成するための処理に基づいて、図４を用いて説明する。図４は、逆フィルタを生成するための処理を示すフローチャートである。

まず、測定状態での測定により、ＥＣＴＦ取得部１３２、及び参照信号取得部１２２が、外耳道伝達特性ＥＣＴＦ、及び第１の参照信号Ｒｅｆ１を取得する（Ｓ１１）。つまり、ユーザＵがマイク２Ｌ、及びヘッドホン４３を装着した装着状態で、頭外定位処理装置１００がインパルス応答測定を実施する。これにより、マイク２Ｌ、及び内蔵マイク４８がそれぞれ収音信号を収音する。そして、特性用収音信号に基づいて、ＥＣＴＦ取得部１３２が外耳道伝達特性を算出する。参照用収音信号に基づいて、参照信号取得部１２２が第１の参照信号Ｒｅｆ１を算出する。

次に、受聴状態での測定により、参照信号取得部１２２が第２の参照信号Ｒｅｆ２を取得する（Ｓ１２）。つまり、ユーザＵがマイク２Ｌを取り外した状態で、ヘッドホン４３を装着する。ユーザＵがヘッドホン４３のみを装着した状態で、インパルス応答測定が実施される。参照用収音信号に基づいて、参照信号取得部１２２が第２の参照信号Ｒｅｆ２を算出する。

次に、変換関数算出部１５１が、第１の参照信号Ｒｅｆ１と第２の参照信号Ｒｅｆ２から、変換関数を算出する（Ｓ１３）。補正部１５２が、外耳道伝達特性ＥＣＴＦに対して、変換関数を適用して、補正後の外耳道伝達特性（以下、補正特性ＡｄＥＣＴＦとする）を算出する（Ｓ１４）。ここでは、対数パワースペクトルに対して、変換関数が適用されている。つまり、補正部１５２が、対数パワースペクトルのみを補正している。逆フィルタ生成部１５３が補正特性ＡｄＥＣＴＦから逆フィルタを算出する（Ｓ１５）。

次に、Ｓ１３〜Ｓ１５の処理について、図５を用いて詳細に説明する。図５は、Ｓ１３〜Ｓ１５の処理の１例を示すフローチャートである。つまり、図５は、変換関数算出部１５１、補正部１５２、及び逆フィルタ生成部１５３における処理を示すフローチャートである。

まず、変換関数算出部１５１が外耳道伝達特性ＥＣＴＦ、第１の参照信号Ｒｅｆ１、第２の参照信号Ｒｅｆ２の対数パワースペクトルを計算し、正規化する（Ｓ２１）。変換関数算出部１５１、補正部１５２、逆フィルタ生成部１５３は、正規化後のデータについて、以下の処理を実施する。

変換関数算出部１５１が第１の参照信号Ｒｅｆ１と第２の参照信号の対数パワースペクトルＰＳＤ１、ＰＳＤ２において、極値ＥＸ１、ＥＸ２を求める（Ｓ２２）。変換関数算出部１５１は、第１の参照信号Ｒｅｆ１の対数パワースペクトルＰＳＤ１の極値ＥＸ１と、第２の参照信号Ｒｅｆ２の対数パワースペクトルＰＳＤ２の極値ＥＸ２を求める。図６に極値ＥＸ１，ＥＸ２を示す。図６は、対数パワースペクトルＰＳＤ１、ＰＳＤ２の一部を示す図であり、横軸が周波数（Ｈｚ）、縦軸がパワー（ｄＢ）となっている。

対数パワースペクトルＰＳＤ１、ＰＳＤ２は通常、複数の極値を有している。ここで、対数パワースペクトルＰＳＤ１の複数の極値ＥＸ１を低周波数側から順にＥＸ１−１、ＥＸ１−２、ＥＸ１−３、・・・・ＥＸ１−Ｎとする。同様に、対数パワースペクトルＰＳＤ２の複数の極値ＥＸ２を低周波数側から順にＥＸ２−１、ＥＸ２−２、ＥＸ２−３、・・・・ＥＸ２−Ｎとする。なお、図６ではＮ＝４となっており、極値ＥＸ１−１〜ＥＸ１−４と極値ＥＸ２−１〜ＥＸ２−４が図示されている。

変換関数算出部１５１は、極値ＥＸ１、ＥＸ２について、対応する極値間の変化ベクトルＶを算出する（Ｓ２３）。つまり、図６に示すように、変換関数算出部１５１は、極値ＥＸ１−１と極値ＥＸ２−１の変化ベクトルを変化ベクトルＶ１として求める。同様に、変換関数算出部１５１は、極値ＥＸ１−２と極値ＥＸ２−２の変化ベクトルを変化ベクトルＶ２として求める。変換関数算出部１５１は、変化ベクトルＶ１〜ＶＮを算出する。図６の例では、Ｎ＝４であるため、４つの変化ベクトルＶ１〜Ｖ４が図示されている。変化ベクトルＶは、周波数と対数パワーとを要素とする２次元ベクトルである。

なお、対数パワースペクトルＰＳＤ１の極値ＥＸ１と、対数パワースペクトルＰＳＤ２の極値ＥＸ２の数が異なる場合、最も近い極値間の変化ベクトルＶを求めればよい。例えば、対数パワースペクトルＰＳＤの極値ＥＸ１の数が、対数パワースペクトルＰＳＤ２の極値ＥＸ２の数よりも小さい場合、対数パワースペクトルＰＳＤ１の極大値に最も近い対数パワースペクトルＰＳＤの極大値をペアとして、変化ベクトルが求められる。同様に、対数パワースペクトルＰＳＤ１の極小値に最も近い対数パワースペクトルＰＳＤの極小値をペアとして、変化ベクトルが求められる。

次に、変換関数算出部１５１は、変化ベクトルＶに基づいて、変換関数を求める（Ｓ２４）。ここでは、変換関数算出部１５１は、ベクトル変換手法を用いて、変換関数を算出している。具体的には、図７に示すように、対数パワースペクトルのグラフ内に格子状の制御点を配置する。そして、複数の変化ベクトルＶ１〜ＶＮに基づいて、制御点のメッシュ構造を変化させる。

例えば、（５，５）にある制御点が、（６，６）に変化した場合、隣接した制御点がその変化に連動して移動するようにする。変換関数算出部１５１は、極値に近い制御点を変化ベクトルＶに基づいて移動させる。変換関数算出部１５１は、複数の変化ベクトルＶ１〜ＶＮに基づいて、対数パワースペクトルのグラフ上における全制御点の移動先を求める。つまり、変換関数算出部１５１は、変換関数（変換用メッシュ）を設定する。本ベクトル変換手法は、例えば、画像の形状変化や３次元データのモーフィング（Ｍｏｒｐｈｉｎｇ）などに用いられるメッシュ構造と同等の手法であるため、詳細な説明は省略する。

補正部１５２は、外耳道伝達特性ＥＣＴＦの対数パワースペクトルＰＳＤに変換関数を適用する（Ｓ２５）。ここでは、説明のため、外耳道伝達特性ＥＣＴＦの対数パワースペクトルをＰＳＤとし、位相スペクトルをＡＳＤとする。補正部１５２は、変換関数を用いて、対数パワースペクトルＰＳＤを補正する。補正後の対数パワースペクトルをＡｄＰＳＤとする。つまり、補正部１５２は、対数パワースペクトルＰＳＤに、変換関数を適用することで、補正後の対数パワースペクトルＡｄＰＳＤを算出する。なお、外耳道伝達特性ＥＣＴＦを補正した補正特性ＡｄＥＣＴＦは、位相スペクトルＡＳＤと、補正後の対数パワースペクトルＡｄＰＳＤとから構成される。

逆フィルタ生成部１５３は、補正後の対数パワースペクトルＡｄＰＳＤに基づいて、逆フィルタを生成する（Ｓ２６）。具体的には、逆フィルタ生成部１５３は、逆離散フーリエ変換又は逆離散コサイン変換等により、補正後の振幅特性（対数パワースペクトルＡｄＰＳＤ）と位相特性（位相スペクトルＡＳＤ）から時間信号を算出する。この時間信号が補正された外耳道伝達特性ＡｄＥＣＴＦ（以下、時間領域の補正特性ＡｄＥＣＴＦとも称する）を示す。逆フィルタ生成部１５３は、時間領域の補正特性ＡｄＥＣＴＦから、振幅と位相が反転する逆フィルタを生成する。逆フィルタの生成方法については公知の手法を用いることができるため説明を省略する。

この逆フィルタが図１で示したフィルタ部４１に設定される。また、右ユニット４３に対して同様の処理を行うことで求められた逆フィルタは、フィルタ部４２に設定される。

図８に、外耳道伝達特性ＥＣＴＦの対数パワースペクトルＰＳＤと補正特性ＡｄＥＣＴＦの対数パワースペクトルＡｄＰＳＤと、を示す。さらに、図８は、第１の参照信号Ｒｅｆ１の対数パワースペクトルＰＳＤ１と、第２の参照信号Ｒｅｆ２の対数パワースペクトルＰＳＤ２と、を示す。このように、対数パワースペクトルＰＳＤ１、ＰＳＤ２に基づく変換関数を、対数パワースペクトルＰＳＤに適用することで、補正後の対数パワースペクトルＡｄＰＳＤを求めることができる。

対数パワースペクトルＰＳＤ１と対数パワースペクトルＰＳＤ２の差は、ヘッドホン４３の装着状態の変化に対応する。補正部１５２が、対数パワースペクトルＰＳＤ１と対数パワースペクトルＰＳＤ２とに基づく変換関数を用いることで、適切に外耳道伝達特性を補正することができる。逆フィルタ生成部１５３が補正後の外耳道伝達特性から逆フィルタを算出している。従って、頭外定位処理装置１００が精度の高い頭外定位処理を行うことができる。

具体的には、頭外定位受聴を行う前の事前測定として、測定状態でのインパルス応答測定を行う。これにより、メモリ１４０に、第１の参照信号Ｒｅｆ１、及び外耳道伝達特性ＥＣＴＦを予め記憶させておくことができる。事前測定が終了したら、ユーザＵがマイクユニット２を取り外す。そして、ユーザＵがヘッドホン４３を装着すると、受聴状態でのインパルス応答測定を実施して、第２の参照信号Ｒｅｆ２を取得する。ヘッドホン４３を装着する毎に、参照信号取得部１２２が、第２の参照信号Ｒｅｆ２を取得することが好ましい。つまり、ユーザＵがヘッドホン４３を装着して、頭外定位受聴を行う前に、参照信号取得部１２２が、第２の参照信号Ｒｅｆ２を取得する。

上記の処理により、頭外定位処理装置１００が、外耳道伝達特性ＥＣＴＦ、第１の参照信号Ｒｅｆ１、及び第２の参照信号Ｒｅｆ２から逆フィルタを算出する。このようにすることで、適応化された逆フィルタを用いて頭外定位処理を行うことができる。つまり、ヘッドホン４３の装着状態が変化しても、頭外定位処理装置１００が、適切な逆フィルタを用いて、頭外定位処理を行うことができる。なお、受聴時にヘッドホン４３を装着したことを自動検知して、受聴状態での測定を自動で行ってもよく。あるいは、ユーザＵが、タッチパネルなど操作部を操作して、受聴状態での測定を行うことを指示してもよい。

上記のように、変換関数算出部１５１は、パワースペクトルに基づいて、変換関数を算出している。したがって、第１の記憶部１４１、及び第２の記憶部１４２は、それぞれパワースペクトルのみを記憶していればよい。すなわち、第１の記憶部１４１と第２の記憶部１４２は、位相スペクトルを記憶していなくてもよい。あるいは、第１の記憶部１４１と第２の記憶部１４２は、時間領域の収音信号をそのまま参照信号として記憶していてもよい。同様に、第３の記憶部１４３は、時間領域の収音信号を外耳道伝達特性として記憶していてもよい。そして、変換関数算出部１５１が変換関数を算出する都度、離散フーリエ変換等を行って、参照信号及び外耳道伝達特性のパワースペクトルを算出するようにしてもよい。

なお、変換関数の算出方法は上記の手法に限定されるものではない。変換関数算出部１５１は、上記したベクトル変換に限らず、様々なパラメータや変換手法によって、変換関数を算出することができる。例えば、対数パワースペクトルＰＳＤ１と対数パワースペクトルＰＳＤ２の差に基づいて、変換関数を求めてもよい。

補正前の外耳道伝達特性ＥＣＴＦの対数パワースペクトルＰＳＤ、第１の参照信号Ｒｅｆ１の対数パワースペクトルＰＳＤ１、第２の参照信号Ｒｅｆ２の対数パワースペクトルＰＳＤ２を用いて、補正後の外耳道伝達特性ＥＣＴＦの対数パワースペクトルＡｄＰＳＤを以下の式（１）を変換関数とすることができる。
ＡｄＰＡＤ＝ＰＳＤ＋（ＰＳＤ２−ＰＳＤ１） ・・・（１）

なお、図２の処理装置２０１は、頭外定位処理装置１００と同じ装置であってもよく、異なる装置であってもよい。頭外定位処理装置１００と、処理装置２０１が異なる装置である場合、処理装置２０１によって取得された外耳道伝達特性ＥＣＴＦ及び第１の参照信号Ｒｅｆ１を、頭外定位処理装置１００が用いればよい。

具体的には、処理装置２０１は、外耳道伝達特性ＥＣＴＦと、第１の参照信号Ｒｅｆ１を取得するために、測定状態でのインパルス応答測定を実施する。そして、処理装置２０１は、外耳道伝達特性ＥＣＴＦと第１の参照信号Ｒｅｆ１を頭外定位処理装置１００に無線又は有線で送信する。あるいは、外耳道伝達特性ＥＣＴＦと第１の参照信号Ｒｅｆ１のデータの一部又は全部は外部に保存されていてもよい。処理装置２０１は、外耳道伝達特性ＥＣＴＦと第１の参照信号Ｒｅｆ１を外部記憶装置やクラウドネットワークなどに記憶させる。頭外定位処理装置１００は、外部記憶装置やクラウドネットワークなどに保存されているデータを読み取る、又は受信することで、外耳道伝達特性ＥＣＴＦと第１の参照信号Ｒｅｆ１を取得する。この場合、メモリ１４０は、外耳道伝達特性ＥＣＴＦと第１の参照信号Ｒｅｆ１と第２の参照信号Ｒｅｆ２とを一時的に記憶するメモリであってもよい。つまり、各データの使用後に、データを消去してもよい。

ヘッドホン４３による受聴を行う前に、頭外定位処理装置１００が受聴状態でのインパルス応答測定を行って、第２の参照信号Ｒｅｆ２を取得する。このように取得された外耳道伝達特性ＥＣＴＦ、第１の参照信号Ｒｅｆ１、第２の参照信号Ｒｅｆ２を用いて、頭外定位処理装置１００が上記の補正処理を実施する。したがって、頭外定位処理装置１００は、ヘッドホン装着状態に適応した逆フィルタを求めることができる。これにより、適切に頭外定位処理を行うことができ、高い精度での音像定位効果を得ることができる。

上記処理のうちの一部又は全部は、コンピュータプログラムによって実行されてもよい。上述したプログラムは、様々なタイプの非一時的なコンピュータ可読媒体（ｎｏｎ−ｔｒａｎｓｉｔｏｒｙ ｃｏｍｐｕｔｅｒ ｒｅａｄａｂｌｅ ｍｅｄｉｕｍ）を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（ｔａｎｇｉｂｌｅ ｓｔｏｒａｇｅ ｍｅｄｉｕｍ）を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば光磁気ディスク）、ＣＤ−ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＣＤ−Ｒ、ＣＤ−Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Ｐｒｏｇｒａｍｍａｂｌｅ ＲＯＭ)、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ ＰＲＯＭ)、フラッシュＲＯＭ、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ））を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（ｔｒａｎｓｉｔｏｒｙ ｃｏｍｐｕｔｅｒ ｒｅａｄａｂｌｅ ｍｅｄｉｕｍ)によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。

以上、本発明者によってなされた発明を実施の形態に基づき具体的に説明したが、本発明は上記実施の形態に限られたものではなく、その要旨を逸脱しない範囲で種々変更可能であることは言うまでもない。

Ｕ ユーザ
１ 被測定者
１０ 頭外定位処理部
１１ 畳み込み演算部
１２ 畳み込み演算部
２１ 畳み込み演算部
２２ 畳み込み演算部
２４ 加算器
２５ 加算器
４１ フィルタ部
４２ フィルタ部
４３ ヘッドホン
４５ ハウジング
４６ ヘッドホンスピーカ
４８ 内蔵マイク
１００ 頭外定位処理装置
１１１ 測定信号生成部
１１２ Ｄ／Ａコンバータ
１２１ Ａ／Ｄコンバータ
１２２ 参照信号取得部
１３１ Ａ／Ｄコンバータ
１３２ ＥＣＴＦ取得部
１４０ メモリ
１４１ 第１の記憶部
１４２ 第２の記憶部
１４３ 第３の記憶部
１５１ 変換関数算出部
１５２ 補正部
１５３ 逆フィルタ生成部

Claims (4)

1. 参照用マイクを有する音声出力部と、
   前記参照用マイクが収音した収音信号に基づく参照信号を取得する参照信号取得部と、
   前記音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記参照信号取得部が取得した第１の参照信号を記憶する第１の記憶部と、
   前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照信号取得部が取得した第２の参照信号を記憶する第２の記憶部と、
   前記測定状態において、前記測定用マイクが収音した収音信号に基づいて取得された伝達特性を記憶する第３の記憶部と、
   前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出する変換関数算出部と、
   前記変換関数を用いて、前記伝達特性の周波数特性を補正する補正部と、
   前記補正された周波数特性に基づいて、フィルタを生成するフィルタ生成部と、
   前記フィルタを用いて、頭外定位処理を行う頭外定位処理部と、を備えた頭外定位処理装置。
2. 前記変換関数算出部は、前記第１の参照信号と前記第２の参照信号とのパワースペクトルの極値間の変化ベクトルを求め、
   前記変化ベクトルに基づいて、前記変換関数を算出する請求項１に記載の頭外定位処理装置。
3. 参照用マイクを有する音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記測定用マイクで収音された収音信号に基づく伝達特性を取得するステップと、
   前記測定状態において、前記参照用マイクで収音された収音信号に基づく第１の参照信号を取得するステップと、
   前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照用マイクで収音された収音信号に基づく第２の参照信号を取得するステップと、
   前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出するステップと、
   前記変換関数を用いて、前記伝達特性の周波数特性を補正するステップと、
   前記補正された周波数特性に基づいて、フィルタを生成するステップと、
   前記フィルタを用いて、頭外定位処理を行うステップと、を備えた頭外定位処理方法。
4. コンピュータに、
   参照用マイクを有する音声出力部を装着したユーザの外耳道に測定用マイクが配置された測定状態において、前記測定用マイクで収音された収音信号に基づく伝達特性を取得するステップと、
   前記測定状態において、前記参照用マイクで収音された収音信号に基づく第１の参照信号を取得するステップと、
   前記測定用マイクが前記外耳道から取り外された受聴状態において、前記参照用マイクで収音された収音信号に基づく第２の参照信号を取得するステップと、
   前記第１の参照信号と前記第２の参照信号とに基づいて、周波数特性の変換関数を算出するステップと、
   前記変換関数を用いて、前記伝達特性の周波数特性を補正するステップと、
   前記補正された周波数特性に基づいて、フィルタを生成するステップと、
   前記フィルタを用いて、頭外定位処理を行うステップと、
   を実行させるプログラム。

Images