JP6753329B2 - フィルタ生成装置、及びフィルタ生成方法 - Google Patents

Description

本発明は、フィルタ生成装置、及びフィルタ生成方法に関する。

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

頭外定位再生においては、２チャンネル（以下、ｃｈと記載）のスピーカから発した測定信号（インパルス音等）を聴取者本人の耳に設置したマイクロフォン（以下、マイクとする）で録音する。そして、インパルス応答で得られた収音信号に基づいて、処理装置がフィルタを作成する。作成したフィルタを２ｃｈのオーディオ信号に畳み込むことにより、頭外定位再生を実現することができる。

特許文献１には、個人化された室内インパルス応答のセットを取得する方法が開示されている。特許文献１では、聴取者の各耳の近くにマイクを設置している。そして、スピーカを駆動した時のインパルス音を、左右のマイクが録音する。

特表２００８−５１２０１５号公報 特許４１８４４２０号公報

インパルス応答測定時には、暗騒音や電源ノイズといった外乱がある。このため、インパルス応答測定では、同じ条件でのインパルス応答測定を複数回行い、マイクで収音された収音信号を同期加算している（特許文献２）。これにより、外乱の影響を排除することができ、Ｓ／Ｎ比を向上させることができる。同期加算を行う場合、同期加算回数が多ければ多い程、外乱の影響が少なくなる。しかしながら、測定時にはユーザは動かずに静止していなければならず、また測定音を何度も聴き続けることは、ユーザにとって負担となる。

本発明は上記の点に鑑みなされたもので、ユーザの負担を軽減しつつ、伝達特性に応じたフィルタを適切に生成することができるフィルタ生成装置、及びフィルタ生成方法を提供することを目的とする。

本発明にかかるフィルタ生成装置は、測定信号を出力する音源から出力された前記測定信号を収音して、収音信号を取得するマイクと、前記収音信号に基づいて、前記音源から前記マイクまでの伝達特性に応じたフィルタを生成するフィルタ生成部と、を備え、前記フィルタ生成部は、受聴者にマイクが装着された状態で取得された前記収音信号を第１の同期加算回数で同期加算することで、第１の同期加算信号を生成する第１の同期加算部と、前記受聴者以外にマイクが装着された状態で取得された前記収音信号を、前記第１の同期加算回数よりも多い第２の同期加算回数で同期加算することで、第２の同期加算信号を生成する第２の同期加算部と、前記第１の同期加算信号に対応する第１のスペクトルと、前記第２の同期加算信号に対応する第２のスペクトルを取得するように、前記第１及び第２の同期加算信号を周波数領域に変換する変換部と、所定の周波数以下の帯域における第２のスペクトルを用いて、前記第１のスペクトルを補正して、第３のスペクトルを生成する補正部と、前記第３のスペクトルを時間領域に逆変換する逆変換部と、を備えたものである。

本発明にかかるフィルタ生成方法は、音源から出力された測定信号をマイクで収音することで伝達特性に応じたフィルタを生成するフィルタ生成方法であって、受聴者にマイクが装着された状態で取得された収音信号を第１の同期加算回数で同期加算することで、第１の同期加算信号を生成するステップと、前記受聴者以外にマイクが装着された状態で取得された収音信号を、前記第１の同期加算回数よりも多い第２の同期加算回数で同期加算することで、第２の同期加算信号を生成するステップと、前記第１の同期加算信号に対応する第１のスペクトルと、前記第２の同期加算信号に対応する第２のスペクトルを取得するように、前記第１及び第２の同期加算信号を周波数領域に変換するステップと、所定の周波数以下の帯域における第２のスペクトルを用いて、前記第１のスペクトルを補正して、第３のスペクトルを生成するステップと、前記第３のスペクトルを時間領域のデータに逆変換するステップと、を備えたものである。

本発明によれば、ユーザの負担を軽減しつつ、伝達特性に応じたフィルタを適切に生成することができるフィルタ生成装置、及びフィルタ生成方法を提供することができる。

本実施の形態に係る頭外定位処理装置を示すブロック図である。 フィルタを生成するフィルタ生成装置の構成を示す図である。 ダミーヘッドを用いて、１６回の同期加算回数で測定された同期加算信号の対数パワースペクトルのＲＡＷデータである。 ダミーヘッドを用いて、６４回の同期加算回数で測定された同期加算信号の対数パワースペクトルのＲＡＷデータである。 ダミーヘッドを用いて、１６回の同期加算回数で測定された同期加算信号に対して処理を行った対数パワースペクトルである。 ダミーヘッドを用いて、６４回の同期加算回数で測定された同期加算信号の対数パワースペクトルである。 同期加算回数が１６回と６４回との定常波減衰率を示すグラフである。 個人測定において、６４回の同期加算回数で測定された同期加算信号の対数パワースペクトルである。 フィルタ生成方法の概要を示すフローチャートである。 補正前の対数パワースペクトルを示すグラフである。 補正後の対数パワースペクトルを示すグラフである。 フィルタ生成装置の構成を示す制御ブロック図である。 フィルタ生成方法を示すフローチャートである。 個人測定の対数パワースペクトルと補正後の対数パワースペクトルを示す例１である。 個人測定の対数パワースペクトルと補正後の対数パワースペクトルを示す例２である。 個人測定の対数パワースペクトルと補正後の対数パワースペクトルを示す例３である。 個人測定の対数パワースペクトルと補正後の対数パワースペクトルを示す例４である。 個人測定の対数パワースペクトルと補正後の対数パワースペクトルを示す例５である。

本実施の形態では、スピーカからマイクまでの伝達特性を測定している。そして、測定された伝達特性に基づいて、フィルタ生成装置がフィルタを生成している。

本実施の形態にかかるフィルタ生成装置で生成したフィルタを用いた音像定位処理の概要について説明する。ここでは、音像定位処理装置の一例である頭外定位処理について説明する。本実施形態にかかる頭外定位処理は、個人の空間音響伝達特性（空間音響伝達関数ともいう）と外耳道伝達特性（外耳道伝達関数ともいう）を用いて頭外定位処理を行うものである。外耳道伝達特性は、外耳道入口から鼓膜までの伝達特性である。本実施形態では、スピーカから聴取者の耳までの空間音響伝達特性、及びヘッドホンを装着した状態での外耳道伝達特性の逆特性を用いて頭外定位処理を実現している。

本実施の形態にかかる頭外定位処理装置は、パーソナルコンピュータ、スマートホン、タブレットＰＣなどの情報処理装置であり、プロセッサ等の処理手段、メモリやハードディスクなどの記憶手段、液晶モニタ等の表示手段、タッチパネル、ボタン、キーボード、マウスなどの入力手段、ヘッドホン又はイヤホンを有する出力手段を備えている。

実施の形態１．
本実施の形態にかかる音場再生装置の一例である頭外定位処理装置１００を図１に示す。図１は、頭外定位処理装置のブロック図である。頭外定位処理装置１００は、ヘッドホン４３を装着するユーザＵに対して音場を再生する。そのため、頭外定位処理装置１００は、ＬｃｈとＲｃｈのステレオ入力信号ＸＬ、ＸＲについて、音像定位処理を行う。ＬｃｈとＲｃｈのステレオ入力信号ＸＬ、ＸＲは、ＣＤ（Compact Disc）プレイヤーなどから出力されるアナログのオーディオ再生信号、又は、mp3(MPEG Audio Layer-3)等のデジタルオーディオデータである。である。なお、頭外定位処理装置１００は、物理的に単一な装置に限られるものではなく、一部の処理が異なる装置で行われてもよい。例えば、一部の処理がパソコンなどにより行われ、残りの処理がヘッドホン４３に内蔵されたＤＳＰ(Digital Signal Processor)などにより行われてもよい。

頭外定位処理装置１００は、頭外定位処理部１０と、フィルタ部４１、フィルタ部４２、及びヘッドホン４３を備えている。

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

空間音響伝達特性は、４つの伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに応じたフィルタを有している。４つの伝達特性に応じたフィルタは、後述するフィルタ生成装置を用いて求めることができる。

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

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

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

フィルタ部４１は、補正されたＬｃｈ信号をヘッドホン４３の左ユニット４３Ｌに出力する。フィルタ部４２は、補正されたＲｃｈ信号をヘッドホン４３の右ユニット４３Ｒに出力する。ユーザＵは、ヘッドホン４３を装着している。ヘッドホン４３は、Ｌｃｈ信号とＲｃｈ信号をユーザＵに向けて出力する。これにより、ユーザＵの頭外に定位された音像を再生することができる。

（フィルタ生成装置）
図２を用いて、空間音響伝達特性（以下、伝達特性とする）を測定して、フィルタを生成するフィルタ生成装置について説明する。図２は、フィルタ生成装置２００の測定構成を模式的に示す図である。なお、フィルタ生成装置２００は、図１に示す頭外定位処理装置１００と共通の装置であってもよい。あるいは、フィルタ生成装置２００の一部又は全部が頭外定位処理装置１００と異なる装置となっていてもよい。

図２に示すように、フィルタ生成装置２００は、ステレオスピーカ５とステレオマイク２を有している。ステレオスピーカ５が測定環境に設置されている。測定環境は、ユーザＵの自宅の部屋やオーディオシステムの販売店舗やショールーム等でもよい。

本実施の形態では、フィルタ生成装置２００の処理装置（図２では不図示）が、伝達特性に応じたフィルタを適切に生成するための演算処理を行っている。処理装置は、例えば、ＭＰ３（ＭＰＥＧ−１ Ａｕｄｉｏ Ｌａｙｅｒ−３）プレイヤー、ＣＤプレイヤー等の音楽プレイヤーなどを有している。処理装置は、パーソナルコンピュータ（ＰＣ）、タブレット端末、スマートホン等であってもよい。

ステレオスピーカ５は、左スピーカ５Ｌと右スピーカ５Ｒを備えている。例えば、受聴者１の前方に左スピーカ５Ｌと右スピーカ５Ｒが設置されている。左スピーカ５Ｌと右スピーカ５Ｒは、インパルス応答測定を行うためのインパルス音等を出力する。

以下、本実施の形態では、音源となるスピーカの数を２（ステレオスピーカ）として説明するが、測定に用いる音源の数は２に限らず、１以上であればよい。すなわち、1chのモノラル、または、5.1ch、7.1ch等の、いわゆるマルチチャンネル環境においても同様に、本実施の形態を適用することができる。

ステレオマイク２は、左のマイク２Ｌと右のマイク２Ｒを有している。左のマイク２Ｌは、受聴者１の左耳９Ｌに設置され、右のマイク２Ｒは、受聴者１の右耳９Ｒに設置されている。具体的には、左耳９Ｌ、右耳９Ｒの外耳道入口又は鼓膜位置にマイク２Ｌ、２Ｒを設置することが好ましい。マイク２Ｌ、２Ｒは、ステレオスピーカ５から出力された測定信号を収音して、収音信号を取得する。例えば、測定信号はインパルス信号やＴＳＰ（Ｔｉｍｅ Ｓｔｒｅｃｈｅｄ Ｐｕｌｅ）信号等でもよい。マイク２Ｌ、２Ｒは収音信号を後述するフィルタ生成装置２００に出力する。受聴者１は、人でもよく、ダミーヘッドでもよい。すなわち、本実施形態において、受聴者１は人だけでなく、ダミーヘッドを含む概念である。

上記のように、左右のスピーカ５Ｌ、５Ｒで出力されたインパルス音をマイク２Ｌ、２Ｒで測定することでインパルス応答が測定される。フィルタ生成装置２００は、インパルス応答測定に基づいて取得した収音信号をメモリなどに記憶する。これにより、左スピーカ５Ｌと左マイク２Ｌとの間の伝達特性Ｈｌｓ、左スピーカ５Ｌと右マイク２Ｒとの間の伝達特性Ｈｌｏ、右スピーカ５Ｌと左マイク２Ｌとの間の伝達特性Ｈｒｏ、右スピーカ５Ｒと右マイク２Ｒとの間の伝達特性Ｈｒｓが測定される。すなわち、左スピーカ５Ｌから出力された測定信号を左マイク２Ｌが収音することで、伝達特性Ｈｌｓが取得される。左スピーカ５Ｌから出力された測定信号を右マイク２Ｒが収音することで、伝達特性Ｈｌｏが取得される。右スピーカ５Ｒから出力された測定信号を左マイク２Ｌが収音することで、伝達特性Ｈｒｏが取得される。右スピーカ５Ｒから出力された測定信号を右マイク２Ｒが収音することで、伝達特性Ｈｒｓが取得される。

そして、フィルタ生成装置２００は、収音信号に基づいて、左右のスピーカ５Ｌ、５Ｒから左右のマイク２Ｌ、２Ｒまでの伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに応じたフィルタを生成する。具体的には、フィルタ生成装置２００は、伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓを所定のフィルタ長で切り出して、演算処理を行う。このようにすることで、フィルタ生成装置２００は、頭外定位処理装置１００の畳み込み演算に用いられるフィルタとして生成する。図１で示したように、頭外定位処理装置１００が、左右のスピーカ５Ｌ、５Ｒと左右のマイク２Ｌ、２Ｒとの間の伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに応じたフィルタを用いて頭外定位処理を行う。すなわち、伝達特性に応じたフィルタをオーディオ再生信号に畳み込むことにより、頭外定位処理を行う。

ここで、電源ノイズや空調等によるいわゆる暗騒音（定在波、定常波）と近い周波数帯域である、低周波数帯域に対して、測定により得られる特性の精度をさらに向上させるための検討について説明する。ここでは、ダミーヘッドによる詳細な測定と、その測定によって得られる特性を用いた各個人の特性データの補正を検討する。

上述した暗騒音等または突発的なノイズ等の外乱による影響を軽減するため、フィルタ生成装置２００は、同期加算を行っている。左のスピーカ５Ｌ又は右のスピーカ５Ｒが同じ測定信号を一定の時間間隔で繰り返し出力する。そして、左マイク２Ｌ、及び右マイク２Ｒが複数の測定信号を収音して、それぞれの測定信号に対応する収音信号を同期して加算している。例えば、同期加算回数が１６回の場合、左のスピーカ５Ｌ又は右のスピーカ５Ｒが測定信号を１６回出力する。そして、左マイク２Ｌ、及び右マイク２Ｒが１６個の収音信号を同期して加算している。このようにすることで、暗騒音等または突発的なノイズ等の外乱による影響を軽減することができ、適切なフィルタを生成することができる。

左のスピーカ５Ｌ又は右のスピーカ５Ｒは、前の測定信号の残響などが無い状態で、次の測定信号を出力する必要がある。よって、測定信号を出力する時間間隔をある程度長くしなければならない。そのため、同期加算回数が多くなると、全体の測定時間が長くなってしまう。受聴者１は、測定中、動かずに静止していなければならない。受聴者１がユーザＵ個人の場合、測定時間を長くすることはユーザＵにとって負担になる。そのため、本実施の形態では、ユーザ個人の測定では、同期加算回数を少なくしている。

一方、同期加算回数を多くすることで、外乱の影響をより少なくすることができる。そのため、ダミーヘッドを用いた測定では、同期加算回数を多くしても、ユーザＵの放胆とはならない。よって、本実施の形態では、ダミーヘッドを用いた測定と、ユーザＵ個人の測定で同期加算回数を変えている。

例えば、受聴者１としてのダミーヘッドにステレオマイク２を装着した状態では、同期加算回数が６４回の測定が行われる。一方、実際のユーザＵにマイク２を装着した状態では、同期加算回数が１６回の測定が行われる。ここで、ダミーヘッドにステレオマイク２を装着した状態で得られた測定をコンフィギュレーション測定とし、コンフィギュレーション測定に基づくデータをコンフィギュレーションデータとする。実際に頭外定位受聴を行うユーザＵにマイク２を装着した状態の測定を個人測定とし、個人測定に基づくデータを個人測定データとする。フィルタ生成装置２００は、個人測定データをコンフィギュレーションデータで補正する。

具体的には、補正上限周波数よりも低い低周波数帯域（補正帯域とも称する）については、個人測定データをコンフィギュレーションデータで補正する。例えば、低周波数帯域については、個人測定データの値（例えば、パワー、又は振幅）をコンフィギュレーションデータの値（例えば、パワー、又は振幅）で置き換える。補正上限周波数よりも高い高周波数帯域については、個人測定データの値をそのまま用いる。このように、フィルタ生成装置２００は、コンフィギュレーションデータと個人測定データを合成することで、伝達特性に応じたフィルタを生成する。本実施の形態では、位相スペクトルは補正せず、パワースペクトルのみを補正している。

個人測定の同期加算回数をコンフィギュレーション測定の同期加算回数よりも少なくすることで、ユーザの負担を軽減することができる。すなわち、個人測定の同期加算回数を少なくすることで、ユーザＵが実際に測定信号を受聴する測定時間を短くすることができる。これにより、ユーザ負担を軽減することができる。また、コンフィギュレーション測定での同期加算回数を多くすることで、フィルタの低周波数帯域を適切に設定することができる。

ここで、同期加算回数による測定データの違いについて、説明する。図３は、同期加算回数を１６回とした測定データを示し、図４は、同期加算回数を６４回とした測定データを示す。図３、図４は、同期加算した同期加算信号を高速フーリエ変換（ＦＦＴ）で解析した対数パワースペクトルを示している。また、図３、図４とも受聴者１としてダミーヘッドを用いた場合の測定データを示している。本実施の形態の測定では、サンプリング周波数は４８ｋＨｚ、測定フレーム長は８１９２サンプルとしている。図３、図４は、８１９２サンプルのデータ（以下、ＲＡＷデータ）の対数パワースペクトルを示している。

図３、図４は、４つの伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓの対数パワースペクトルを示している。図３では、１６回の同期加算回数を１セットとして、５セットの測定を行った結果を示し、図４は６４回の同期加算回数を１セットとして、５セットの測定を行った結果を示している。したがって、図３、図４のそれぞれには、伝達測定Ｈｌｓについて、５つの対数パワースペクトルが示されている。同様に、伝達特性、Ｈｌｏ、Ｈｒｏ、Ｈｒｓについても、それぞれ５つの対数パワースペクトルが示されている。図３、図４には、それぞれ２０個の対数パワースペクトルが示されている。

図３、図４中の丸枠に囲まれた箇所から見て取れるように、約４０Ｈｚ〜２００Ｈｚの周波数帯域において、１６回の同期加算回数よりも６４回の同期加算回数のほうが、伝達特性が安定しており、精度が高いことがわかる。すなわち、同期加算回数が１６回の場合、図３に示すように、約４０Ｈｚ〜２００Ｈｚの周波数帯域において、セット毎のばらつきが大きくなる。

図５、図６は、マイク特性の補正、４０９６サンプル長へのフィルタの切り出し、及び窓掛けの処理が行われた同期加算信号の対数パワースペクトルを示している。図５は、同期加算回数が１６回の測定データ、すなわち、図３に対応するＲＡＷデータに処理を行った対数パワースペクトルを示している。図６は、同期加算回数が６４回の測定データ、すなわち、図４に対応するＲＡＷデータに処理を行った対数パワースペクトルを示している。

この場合も、図５、図６中の丸枠に囲まれた箇所から見て取れるように、約４０Ｈｚ〜２００Ｈｚの周波数帯域において、１６回の同期加算回数よりも６４回の同期加算回数のほうが、伝達特性が安定しており、精度が高いことがわかる。すなわち、同期加算回数が１６回の場合、図５に示すように、約４０Ｈｚ〜２００Ｈｚの周波数帯域において、セット毎のばらつきが大きくなる。

図７に、同期加算による定常波減衰率を示す。図７は、サンプリング周波数48kHz、同期フレーム内サンプル数8192、の場合における、純音１Ｈｚから２００Ｈｚまでの１Ｈｚ毎の定常波減衰率を示している。また、図７では、同期加算回数が１６回と６４回の場合の定常波減衰率を示している。同期加算回数が６４回の場合、概ね−２０ｄＢ以上の減衰率が得られることが分かる。すなわち、同期加算回数が６４回の場合、外乱による定常波が十分に減衰していることが分かる。さらに、同期加算回数が１６回と比較すると、６４回の場合、全体的に数十ｄＢの改善が得られていることが分かる。したがって、２００Ｈｚ以下の低周波数帯域において、同期加算回数を６４回とすることで、外乱の影響を十分低減することができる。

暗騒音の周波数帯域と近い低周波数帯域の測定制度を向上させるためには、同期加算回数を多くすることが好ましい。そこで、本実施の形態では、低周波数帯域については、ダミーヘッドを用いたコンフィギュレーション測定を行うことで、同期加算回数を多くしている。すなわち、ダミーヘッドにステレオマイク２を装着した状態で伝達特性の測定を行うことで、同期加算回数を多くした場合でも、ユーザの負担を軽減することができる。そして、フィルタ生成装置２００は、個人測定データをコンフィギュレーションデータで補正する。

ここで、個人測定データの一例を図８に示す。図８は、受聴者１をユーザＵとした場合の測定結果を示すグラフである。図８では、図６と同様に、マイク特性の補正、４０９６サンプル長へのフィルタの切り出し、窓掛けを行ったデータに対して、ＦＦＴで解析した対数パワースペクトルを示している。図８は、同期加算回数を６４回にした時の個人測定データを示している。

図６と図８とを比較してわかるように、低周波数帯域の対数パワースペクトルの形状は、コンフィギュレーションデータと個人測定データとの間で同等となっている。理論的にも、低周波数帯域の頭部伝達関数は、各個人でほとんど差がないことが分かっている。すなわち、低周波数帯域の対数パワースペクトルの形状は、ユーザＵによる個人差がほとんどない。したがって、低周波数帯域の個人測定データを、コンフィギュレーションデータで補正することが可能である。

ここで、図６、図８等に示した対数パワースペクトルにおいて、同期加算信号の時間波形におけるサンプル値の二乗総和（＝セグメンタルパワー）について、伝達特性Ｈｌｓ，Ｈｒｓのいずれか大きい方が１となるように、データを正規化している。すなわち、４つの伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓに同じ係数を乗じることで、正規化が行われている。しかしながら、正規化を行ったとしても、図６、図８の丸枠に示すように、低周波数帯域のレベルに違いが生じていることが分かる。

そこで、本実施の形態では、調整帯域におけるコンフィギュレーションデータと個人測定データに応じて、レベル調整を行うことが好ましい。調整帯域は、補正上限周波数よりも高い周波数を含んでいる。調整帯域は、例えば、２００Ｈｚ〜５００Ｈｚとなっている。すなわち、このレベル調整の詳細については後述する。

次に、本実施の形態にかかるフィルタ生成方法について、図９を用いて説明する。図９は、フィルタ生成方法の概要を示すフローチャートである。

まず、コンフィギュレーション測定を行うため、フィルタ生成装置２００が、ダミーヘッドを用いて、同期加算回数６４回での測定を行う（Ｓ１１）。すなわち、図２に示した測定環境において、受聴位置にダミーヘッドを設置して、ダミーヘッドにステレオマイク２を装着する。そして、ステレオスピーカ５が、同じ測定信号を６４回出力する。ステレオマイク２が収音した６４個の収音信号を同期して加算する。これにより、伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓのそれぞれに対応する同期加算信号が取得される。

次に、フィルタの切り出しを行う（Ｓ１２）。例えば、Ｓ１１で取得した同期加算信号に対して、４０９６サンプル長へのフィルタの切り出し処理を前処理として行う。同期加算信号は部屋の残響などを考慮し、十分に長い時間のデータであるため、フィルタ生成装置２００は、必要なサンプル数のデータ長に切り出しを行う。なお、フィルタ生成装置２００が、切り出されたフィルタに対して、ＤＣ成分のカット、マイク特性の補正、及び窓掛け等の処理を前処理として行ってもよい。

そして、フィルタ生成装置２００が、前処理されたデータをコンフィギュレーションデータとして保存する（Ｓ１３）。具体的には、フィルタ生成装置２００が前処理されたコンフィギュレーションデータを、周波数領域のデータに変換する。フィルタ生成装置２００が周波数領域のデータをコンフィギュレーションデータとして保存する。例えば、フィルタ生成装置２００は、ＦＦＴを行うことで、対数パワースペクトルと位相スペクトルを算出する。対数パワースペクトルと位相スペクトルとがコンフィギュレーションデータとしてメモリなどに保存される。

次に、個人測定データを取得するため、ユーザＵにステレオマイク２を装着して、同期加算回数１６回での測定を行う（Ｓ２１）。すなわち、ユーザＵが図２で示した測定環境の受聴位置に座り、ステレオマイク２を装着する。そして、ステレオスピーカ５が、同じ測定信号を１６回出力する。ステレオマイク２が収音した１６個の収音信号を同期して加算する。これにより、伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓのそれぞれに対応する同期加算信号が取得される。

次に、フィルタの切り出しを行う（Ｓ２２）。例えば、Ｓ２１で取得した同期加算信号に対して、４０９６サンプル長へのフィルタの切り出し処理を前処理として行う。同期加算信号は部屋の残響などを考慮し、十分に長い時間のデータであるため、フィルタ生成装置２００は、必要なサンプル数のデータ長に切り出しを行う。なお、フィルタ生成装置２００が、切り出されたフィルタに対して、ＤＣ成分のカット、マイク特性の補正、及び窓掛け等の処理を前処理として行ってもよい。

次に、フィルタ生成装置２００が、個人測定データに対して、コンフィギュレーションデータを用いた補正を行う（Ｓ２３）。そのため、まず、フィルタ生成装置２００は、Ｓ２２で前処理された個人測定データを周波数領域のデータに変換する。例えば、フィルタ生成装置２００は、ＦＦＴを行うことで、対数パワースペクトルと位相スペクトルを算出する。

そして、個人測定データの対数パワースペクトルをコンフィギュレーションデータの対数パワースペクトルで補正する。具体的には、フィルタ生成装置２００は、補正上限周波数よりも低い低周波数帯域では、個人測定データのパワー値をコンフィギュレーションデータのパワー値で置き換える。フィルタ生成装置２００は、補正上限周波数よりも高い高周波数帯域では、個人測定データのパワー値をそのまま用いる。すなわち、フィルタ生成装置２００は、低周波数帯域のコンフィギュレーションデータのパワー値と、高周波数帯域の個人測定データのパワー値とを組み合わせることで、補正データを生成している。

なお、フィルタ生成装置２００は、補正を行う際において、個人測定データとコンフィギュレーションデータとのレベルを調整してもよい。具体的には、調整帯域の個人測定データとコンフィギュレーションデータとの対数パワースペクトルに基づいて、コンフィギュレーションデータの対数パワースペクトルのレベル調整を行っている。調整帯域は、第１の周波数と第２の周波数との間の帯域である。第１の周波数は第２の周波数より高く、かつ、上記した、補正上限周波数よりも高い。また、ここでは、第２の周波数を補正上限周波数よりも高くしているが、第１の周波数は補正上限周波数よりも低い周波数であってもよい。

図１０、図１１に補正前の対数パワースペクトルと補正後の対数パワースペクトルの一例を示す。図１０では、補正前の個人測定データが実線で示され、コンフィギュレーションデータが破線で示されている。図１１では、補正後のデータが実線で示され、コンフィギュレーションデータが破線で示されている。低周波数帯域において、補正後の対数パワースペクトルとコンフィギュレーションデータは一致する。

具体的な一例では、補正上限周波数は１５０Ｈｚ、第１の周波数は５００Ｈｚ、第２の周波数は２００Ｈｚである。すなわち、調整帯域は、２００Ｈｚ〜５００Ｈｚとなる。フィルタ生成装置２００は、個人測定データにおける１５０Ｈｚ以下のパワー値を、コンフィギュレーションンデータで置換する。個人測定データを補正する低周波数帯域は、最低周波数から１５０Ｈｚの帯域である。個人測定データを補正しない高周波数帯域は補正上限周波数よりも高い帯域である。補正上限周波数は１００Ｈｚ以上、２００Ｈｚ以下とすることが好ましい。

次に、フィルタ生成装置２００の処理装置と、その処理について詳細に説明する。図１２は、フィルタ生成装置２００の処理装置２１０を示す制御ブロック図である。図１３は、処理装置２１０における処理を示すフローチャートである。

処理装置２１０は、フィルタ生成装置（フィルタ生成部）として機能する。処理装置２１０は、測定信号生成部２１１、収音信号取得部２１２、第１の同期加算部２１３、第２の同期加算部２１４、波形切り出し部２１５、ＤＣカット部２１６、第１の窓掛部２１７、正規化部２１８、位相合わせ部２１９、第１の変換部２２０、レベル調整部２２１、第１の補正部２２２、第１の逆変換部２２３、第２の窓掛部２２４、第２の変換部２２５、第２の補正部２２６、第２の逆変換部２２７、及び第３の窓掛部２２８を備えている。

例えば、処理装置２１０は、パーソナルコンピュータ、スマートホン、タブレット端末などの情報処理装置であり、音声入力インターフェース（ＩＦ）と音声出力インターフェースを備えている。すなわち、処理装置２１０は、ステレオマイク２、及びステレオスピーカ５に接続される入出力端子を有する音響デバイスである。

測定信号生成部２１１は、Ｄ／Ａ変換器やアンプなどを備えており、測定信号を生成する。測定信号生成部２１１は、生成した測定信号をステレオスピーカ５にそれぞれ出力する。左スピーカ５Ｌと右スピーカ５Ｒがそれぞれ伝達特性を測定するための測定信号を出力する。左スピーカ５Ｌによるインパルス応答測定と、右スピーカ５Ｒによるインパルス応答測定がそれぞれ行われる。測定信号はインパルス音等の測定音を含んでいる。

ステレオマイク２の左マイク２Ｌ、右マイク２Ｒがそれぞれ測定信号を収音し、収音信号を処理装置２１０に出力する。収音信号取得部２１２は、左マイク２Ｌ、右マイク２Ｒからの収音信号を取得する。なお、収音信号取得部２１２は、Ａ／Ｄ変換器、及びアンプなどを有しており、左マイク２Ｌ、右マイク２Ｒからの収音信号をＡ／Ｄ変換、増幅などしてもよい。収音信号取得部２１２は、取得した収音信号を第１の同期加算部２１３又は第２の同期加算部２１４に出力する。

個人測定の場合、測定信号生成部２１１は、１６回の測定信号を繰り返し、左スピーカ５Ｌ又は右スピーカ５Ｒに出力する。そして測定信号生成部２１１は、１６回の測定信号に対応する収音信号を第１の同期加算部２１３に出力する。第１の同期加算部２１３は、１６回の収音信号を同期加算することで、第１の同期加算信号を生成する。第１の同期加算部２１３は、それぞれ伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓ毎に同期加算信号を生成する。

コンフィギュレーション測定の場合、測定信号生成部２１１は、６４回の測定信号を繰り返し、左スピーカ５Ｌ又は右スピーカ５Ｒに出力する。そして測定信号生成部２１１は、６４回の測定信号に対応する収音信号を第２の同期加算部２１４に出力する。第２の同期加算部２１４は、６４回の収音信号を同期加算することで、第２の同期加算信号を生成する。第２の同期加算部２１４は、それぞれ伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓ毎に同期加算信号を生成する。

第１の同期加算信号は、個人測定データとなり、第２の同期加算信号はコンフィギュレーションデータとなる。

次に、波形切り出し部２１５が、第１及び第２の同期加算信号から必要なデータサンプル長の波形を切り出す（Ｓ３１）。具体的には、８１９２サンプル長の第１及び第２の同期加算信号から４０９６サンプル長のデータを切り出す。

ＤＣカット部２１６は、切り出し後の第１及び第２の同期加算信号のＤＣ成分（直流成分）をカットする（Ｓ３２）。これにより、第１及び第２の同期加算信号のＤＣノイズ成分が除去される。

第１の窓掛部２１７は、ＤＣ成分カット後の第１及び第２の同期加算信号に対して、第１の窓掛けを行う（Ｓ３３）。窓関数は、同期加算信号の絶対最大値を基準に前後窓長の異なる窓関数の半分を掛ける。例えば、窓関数はハニング窓でもよいし、ハミング窓でもよい。また、全体に窓関数をかけずに、両端の一部にのみ窓関数をかけてもよい。第１の窓掛部２１７で用いられる窓関数は特に限定されるものではない。

なお、Ｓ３１からＳ３３の処理は、第１の同期加算信号及び第２の同期加算信号に対して同じとなっている。すなわち、切り出すサンプル長と窓関数は、第１の同期加算信号と第２の同期加算信号の間で同じとなっている。また、第１の同期加算信号と第２の同期加算信号との処理順は特に限定されるものではない。第２の同期加算信号に対してＳ３１〜Ｓ３３の前処理を行った後、第１の同期加算信号に対してＳ３１〜Ｓ３３の前処理を行ってもよい。あるいは、第１の同期加算信号に対してＳ３１〜Ｓ３３の前処理を行った後、第２の同期加算信号に対してＳ３１〜Ｓ３３の前処理を行ってもよい。すなわち、第１の同期加算信号に対して第２の同期加算信号よりも先にＳ３１〜Ｓ３３の前処理を行ってもよく、第２の同期加算信号に対して第１の同期加算信号よりも先にＳ３１〜Ｓ３３の前処理を行ってもよい。

次に、正規化部２１８は、窓掛処理後の同期加算信号に対して正規化を行う（Ｓ３４）。具体的には、正規化部２１８は、伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓの４つの同期加算信号についてデータの二乗総和を求める。正規化部２１８は、４つの二乗総和の中の最大値が１となるような係数を求める。正規化部２１８は、伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓの４つの同期加算信号に対して、その係数を掛ける。例えば、第１の同期加算信号において、伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓの係数Ｋ１は同じ値である。第２の同期加算信号において、伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓの係数Ｋ２は同じ値である。

位相合わせ部２１９は、正規化後の第１の同期加算信号と第２の同期加算信号の位相合わせを行う（Ｓ３５）。具体的には、位相合わせ部２１９は、伝達特性Ｈｌｓ、Ｈｌｏ，Ｈｒｏ、Ｈｒｓのそれぞれに対して、絶対最大値をもつサンプル位置を求める。そして、第１の同期加算信号と第２の同期加算信号とにおいて、絶対最大値を持つサンプル位置が同じになるように、第２の同期加算信号をシフトする。

例えば、伝達特性Ｈｌｓの第１の同期加算信号と、伝達特性Ｈｌｓの第２の同期加算信号の位相合わせを行う場合を説明する。伝達特性Ｈｌｓの第１の同期加算信号の絶対最大値がサンプル位置Ｎ１であり、伝達特性Ｈｌｓの第２の同期加算信号の絶対最大値がサンプル位置Ｎ２であるとする。この場合、第１の同期加算信号と第２の同期加算信号の絶対最大値がサンプル位置Ｎ１で一致するよう、第２の同期加算信号を（Ｎ１−Ｎ２）だけシフトする。

同様に伝達特性Ｈｌｏについても、第１の同期加算信号と第２の同期加算信号の絶対最大値が一致するように、第２の同期加算信号をシフトする。伝達特性Ｈｒｏについても、第１の同期加算信号と第２の同期加算信号の絶対最大値が一致するように、第２の同期加算信号をシフトする。伝達特性Ｈｒｓについても、第１の同期加算信号と第２の同期加算信号の絶対最大値が一致するように、第２の同期加算信号をシフトする。なお、位相合わせの方法は上記の手法に限らず、第１の同期加算信号と第２の同期加算信号の相関などを用いてもよい。

次に、第１の変換部２２０は、位相合わせ後の第１及び第２の同期加算信号を周波数領域のデータに変換する（Ｓ３６）。第１の変換部２２０は、ＦＦＴを用いて、第１の同期加算信号を第１の対数パワースペクトル及び第１の位相スペクトルを生成する。同様に第１の変換部２２０は、ＦＦＴを用いて、第２の同期加算信号を第２の対数パワースペクトル及び第２の位相スペクトルを生成する。

第１の対数パワースペクトル及び第１の位相スペクトルは、個人測定データであり、第２の対数パワースペクトル及び第２の位相スペクトルは、コンフィギュレーションデータである。なお、第１の変換部２２０は、対数パワースペクトルの代わりに振幅スペクトルを生成してもよい。また、第１の変換部２２０は、離散フーリエ変換や離散コサイン変換により、同期加算信号を周波数領域のデータに変換してもよい。

レベル調整部２２１は、対数パワースペクトルの基準値に基づいて、コンフィギュレーションデータのレベル調整を行う（Ｓ３７）。具体的には、レベル調整部２２１は、第１の対数パワースペクトルと第２の対数パワースペクトルの基準値を求める。基準値は、例えば、所定の周波数範囲における対数パワースペクトルの平均値である。なお、レベル調整部２２１は、一定値以上の外れ値を、除外してもよい。あるいは、レベル調整部２２１は、一定値以上の外れ値を一定値に制限してもよい。なお、基準値の算出方法は、これに限られるものではない。例えば、ケプストラムスムージング、移動平均、直線近似等によるスムージングや変換を施したデータの平均値を基準値として用いることも可能であり、または、それらの中央値を基準値として用いることができる。

レベル調整部２２１は、第１の対数パワースペクトルの基準値を第１の基準値として算出し、第２の対数パワースペクトルの基準値を第２の基準値として算出する。そして、レベル調整部２２１は、第１の基準値及び第２の基準値に基づいて、第２の対数パワースペクトルのレベル調整を行う。具体的には、第２の基準値が、第１の基準値と一致するように第２の対数パワースペクトルのパワー値を調整する。例えば、第１の基準値と第２の基準値の比に応じた係数Ｋ３を第２の対数パワースペクトルに加算または減算する。なお、対数パワースペクトルに代えて振幅スペクトルを用いる場合は、係数Ｋ３を乗算することで振幅値を調整する。係数Ｋ３には、周波数によらない一定値を用いることができる。このようにして、レベル調整部２２１は、第１の対数パワースペクトルに基づいて、第２の対数パワースペクトルのレベル調整を行う。

第１の補正部２２２は、レベル調整後に対数パワースペクトルを用いて、第１の対数パワースペクトルを補正する（Ｓ３８）。具体的には、第１の対数パワースペクトルの低周波数帯域のパワー値を第２の対数パワースペクトルのパワー値に置き換える。これにより、図１０に示した対数パワースペクトルが図１１に示す対数パワースペクトルに補正される。なお、低周波数帯域とは、上記の通り、補正上限周波数以下の帯域である。例えば、補正上限周波数は１５０Ｈｚであるため、低周波数帯域は最低周波数〜１５０Ｈｚとなる。補正上限周波数よりも高い高周波数帯域では、第１の補正部２２２が第１の対数パワースペクトルのパワー値を補正せずに、そのまま用いる。なお、第１の補正部２２２により補正された対数パワースペクトルを第１の補正データ、又は第３の対数パワースペクトルとも称する。

第１の逆変換部２２３が第３の対数パワースペクトルを時間領域に逆変換する（Ｓ３９）。具体的には、第１の逆変換部２２３が逆高速フーリエ変換（ＩＦＦＴ）を用いて、第１の補正データを時間領域に逆変換する。例えば、第１の逆変換部２２３が第３の対数パワースペクトルと第１の位相スペクトルとに逆離散フーリエ変換を施すことで、第１の補正データが時間領域のデータとなる。第１の逆変換部２２３は、逆離散フーリエ変換ではなく、逆離散コサイン変換等により、逆変換を行ってもよい。

第２の窓掛部２２４は、逆変換された第１の補正データに第２の窓掛けを施す（Ｓ４０）。第２の窓掛けの処理は、Ｓ３３の第１の窓掛けの処理と同じ処理であるため、説明を省略する。第２の窓掛けに用いる窓関数は、第１の窓掛けに用いる窓関数と同じものでもよく、異なるものでもよい。

第２の変換部２２５は、第２の窓掛け後の第１の補正データを周波数領域に変換する（Ｓ４１）。第２の変換部２２５は、第１の変換部２２０と同様に、ＦＦＴを用いて、時間領域における第２の窓掛け後の第１の補正データを周波数領域の第１の補正データに変換する。第２の変換部２２５が算出した対数パワースペクトル及び位相スペクトルを第４の対数パワースペクトル、及び第４の位相スペクトルとする。第４の対数パワースペクトル、及び第４の位相スペクトルは、第２の窓掛け後の対数パワースペクトル、及び第４の位相スペクトルである。

次に、第２の補正部２２６は、第２の窓掛けによる減衰率から、第３の対数パワースペクトルを補正する（Ｓ４２）。具体的には、第２の補正部２２６は、Ｓ３８で算出された第３の対数パワースペクトルと、Ｓ４１で算出された第４の対数パワースペクトルとのパワーの減衰率を求める。第２の補正部２２６は、第２の窓掛け前後の第１の補正データを比較して、所定の周波数帯域におけるパワーの減数率を算出する。そして、第２の補正部２２６は、減衰率に応じて、第３の対数パワースペクトルに対する第２の補正を行う。なお、第２の補正部２２６によって補正された対数パワースペクトルを第５の対数パワースペクトル、又は第２の補正データとする。

ここで、減衰率を算出するための周波数帯域を算出用帯域とする。算出用帯域は、対数パワースペクトルの一部の帯域である。算出用帯域は、同期加算信号のサンプル数やサンプリングレートを用いて求めることができる。算出用帯域は、所定の周波数よりも低い周波数の帯域である。算出用帯域は、低周波数帯域と異なる帯域となっていてもよく、同じ帯域であってもよい。

第２の補正部２２６は、算出用帯域における第３の対数パワースペクトルのパワー値と第４の対数パワースペクトルのパワー値を比較することで、第２の窓掛け処理による減衰率を求める。そして、算出用帯域に第３の対数パワースペクトルのパワー値を減衰率に応じて底上げする。例えば、算出用帯域における第３の対数パワースペクトルのパワー値に減衰率に応じた値を加えたり、乗じたりすることで、算出用帯域における第３の対数パワースペクトルのパワー値が底上げされる。具体的には、第４の対数パワースペクトルと第５の対数パワースペクトルとの減衰率が１となるように、第２の補正部２２６は第３の対数パワースペクトルを補正する。

そして、第２の逆変換部２２７は、第５の対数パワースペクトルを時間領域に逆変換する（Ｓ４３）。第２の逆変換部２２７は、Ｓ３９と同様に逆離散フーリエ変換等を行って、第２の補正データを時間領域に変換する。例えば、第２の逆変換部２２７が第５の対数パワースペクトルと第１の位相スペクトルとに逆離散フーリエ変換を施すことで、第２の補正データが時間領域のデータとなる。第２の逆変換部２２７は、逆離散フーリエ変換ではなく、逆離散コサイン変換により、逆変換を行ってもよい。

そして、第３の窓掛部２２８は、時間領域の第２の補正データに対して、窓掛けを行う（Ｓ４４）。第３の窓掛部２２８は、Ｓ４０の窓掛けと同じ窓関数を用いて、窓掛けを行う。これにより、処理が終了する。

上記の処理を行うことで、処理装置２１０が伝達特性に応じたフィルタを生成することができる。低周波数帯域の特性は、周波数帯域が近い、電源ノイズや空調等によるいわゆる暗騒音（定在波、定常波）の影響を排除することが難しい。また、低周波数帯域の特性は、個人差が小さい。よって、低周波数帯域については、コンフィギュレーションデータで個人測定データを置き換えている。これより、伝達特性に応じたフィルタを適切に生成することができる。処理装置２１０は、伝達特性Ｈｌｓ、Ｈｌｏ、Ｈｒｏ、Ｈｒｓ毎にフィルタを生成する。そして、処理装置２１０が生成したフィルタを、図１の畳み込み演算部１１、１２、２１、２２に設定する。このようにすることで、適切に頭外定位処理することができる。

頭外定位処理装置１００のユーザＵは短時間の簡易な測定のみでよいため、ユーザＵの負担を軽減することができる。上記のフィルタを用いた結果、頭外定位された再生音の音質を向上することができる。これにより、聴感上、（１）耳の周囲に残る低周波数帯域の音像が明瞭となる、（２）左右の偏りが補正され違和感が減少する、（３）中低域の音圧バランスが良くなる、等の効果が得られる。

個人測定データの対数パワースペクトルと補正後の対数パワースペクトルを図１４〜図１８に示す。図１４〜図１８は異なる５人のユーザＵに対して測定された個人測定データの対数パワースペクトルと、補正後の対数パワースペクトルを示している。図１４〜図１８において、太線が補正後の対数パワースペクトルであり、細線が補正前の個人測定スペクトルである。また、図１４〜図１８では、同じコンフィギュレーションデータが用いられている。図１４〜図１８から、低周波数帯域の特性のばらつきが、補正処理により安定化していることが分かる。

なお、第１の補正部２２２が、低周波数帯域のパワー値を置き換えることで、第１の補正を行っていたが、補正する方法は、特に限定されるものではない。補正上限周波数の近傍に境界周波数帯域を設定して、境界周波数帯域において、指数関数的あるいは、線形的にパワー値を漸近させてもよい。

例えば、補正上限周波数を２００Ｈｚとし、２００Ｈｚ〜１ｋＨｚを境界周波数帯域とすることができる。２００Ｈｚ以下の低周波数帯域では、第１の対数パワースペクトルのパワー値を第２の対数パワースペクトルのパワー値で置換する。１ｋＨｚ以上では、第１の対数パワースペクトルのパワー値をそのまま用いる。境界周波数帯（２００Ｈｚ〜１ｋＨｚ）では、２００Ｈｚのパワー値と１ｋＨｚのパワー値を漸近的につなぎ合わせる関数に基づいて、パワー値を設定する。この関数は、例えば、指数関数や線形関数とすることができる。

さらには、個人測定に応じて、補正上限周波数を可変とすることも可能である。例えば、一定の周波数幅を指定して、その周波数幅の範囲内で、第１の対数パワースペクトルと第２の対数パワースペクトルの差異が最小となる周波数点を探索する。探索された周波数点を補正上限周波数とすることができる。例えば、周波数幅が５０Ｈｚとして探索した場合において、８０Ｈｚ〜１３０Ｈｚの周波数幅で、第１の対数パワースペクトルと第２の対数パワースペクトルの差異が最小となったとする。この場合、補正上限周波数を１３０Ｈｚとすることができる。

コンフィギュレーション測定での同期加算回数を６４回、個人測定での同期加算回数を１６回としたが、それぞれの同期加算回数はこれに限られるものではない。すなわち、コンフィギュレーション測定での同期加算回数が個人測定での同期加算回数よりも多ければよい。個人測定での同期加算回数は２回以上であればよい。

個人測定での同期加算回数をコンフィギュレーション測定での同期加算回数よりも少なくすることで、個人測定時間を短縮することができる。よって、ユーザＵの負担を軽減することができる。

ダミーヘッドを用いることで、同期加算回数を多くすることができるので、外乱などの影響を低減することができる。なお、ダミーヘッドを用いてコンフィギュレーション測定を行うことで、ユーザＵの負担を軽減することができるが、コンフィギュレーション測定は、個人測定を行った個人（ユーザＵ）と異なる個人であってもよい。すなわち、１人のコンフィギュレーションデータを複数のユーザＵに使用するようにしてもよい。このようにしても、ユーザＵの負担を軽減することができる。

処理装置２１０において実施された全ての処理は必須ではない。例えば、Ｓ３１〜Ｓ３４の処理、及びＳ３５などの処理の一部又は全部は省略することが可能である。また、レベル調整部２２１によるＳ３７を行うことで、適切にフィルタを生成することができるが、適宜省略することも可能である。Ｓ４０〜Ｓ４４等の処理の一部又は全部を省略することも可能である。

なお、処理装置２１０は、物理的な単一な装置に限られるものではない。すなわち、処理装置２１０の一部の処理を他の装置で行うことも可能である。例えば、他の装置で測定したコンフィギュレーションデータを用意しておく。そして、処理装置２１０は、コンフィギュレーションデータの第２の対数パワースペクトルをメモリなどに格納しておく。コンフィギュレーションデータを予めメモリに格納しておくことで、複数のユーザＵの個人測定データの補正に用いることができる。

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

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

Ｕ ユーザ
１ 受聴者
２Ｌ 左マイク
２Ｒ 右マイク
５Ｌ 左スピーカ
５Ｒ 右スピーカ
９Ｌ 左耳
９Ｒ 右耳
１０ 頭外定位処理部
１１ 畳み込み演算部
１２ 畳み込み演算部
２１ 畳み込み演算部
２２ 畳み込み演算部
２４ 加算器
２５ 加算器
４１ フィルタ部
４２ フィルタ部
４３ ヘッドホン
１００ 頭外定位処理装置
２００ フィルタ生成装置
２１０ 処理装置
２１１ 測定信号生成部
２１２ 収音信号取得部
２１３ 第１の同期加算部
２１４ 第２の同期加算部
２１５ 波形切り出し部
２１６ ＤＣカット部
２１７ 第１の窓掛部
２１８ 正規化部
２１９ 位相合わせ部
２２０ 第１の変換部
２２１ レベル調整部
２２２ 第１の補正部
２２３ 第１の逆変換部
２２４ 第２の窓掛部
２２５ 第２の変換部
２２６ 第２の補正部
２２７ 第２の逆変換部
２２８ 第３の窓掛部

Claims (8)

1. 測定信号を出力する音源から出力された前記測定信号を収音して、収音信号を取得するマイクと、
   前記収音信号に基づいて、前記音源から前記マイクまでの伝達特性に応じたフィルタを生成するフィルタ生成部と、を備え、
   前記フィルタ生成部は、
   受聴者にマイクが装着された状態で取得された前記収音信号を第１の同期加算回数で同期加算することで、第１の同期加算信号を生成する第１の同期加算部と、
   前記受聴者以外にマイクが装着された状態で取得された前記収音信号を、前記第１の同期加算回数よりも多い第２の同期加算回数で同期加算することで、第２の同期加算信号を生成する第２の同期加算部と、
   前記第１の同期加算信号に対応する第１のスペクトルと、前記第２の同期加算信号に対応する第２のスペクトルを取得するように、前記第１及び第２の同期加算信号を周波数領域に変換する変換部と、
   所定の周波数以下の帯域における第２のスペクトルを用いて、前記第１のスペクトルを補正して、第３のスペクトルを生成する補正部と、
   前記第３のスペクトルを時間領域に逆変換する逆変換部と、を備えたフィルタ生成装置。
2. 前記所定の周波数よりも高い第１の周波数と、前記第１の周波数よりも低い第２の周波数との間の帯域に含まれるデータを用いて、前記第２のスペクトルのレベルを調整する調整部をさらに備え、
   前記補正部が前記調整部で調整された前記第２のスペクトルのデータを用いて、前記所定の周波数以下の帯域における前記第１のスペクトルのデータを補正する請求項１に記載のフィルタ生成装置。
3. 前記補正部が、前記所定の周波数以下の帯域における前記第１のスペクトルのデータを、前記第２のスペクトルのデータに置換する請求項１、又は２に記載のフィルタ生成装置。
4. 前記第２の同期加算部が、ダミーヘッドに前記マイクを装着した状態で取得された前記収音信号を加算することで、前記第２の同期加算信号を生成する請求項１から３のいずれか１項に記載のフィルタ生成装置。
5. 音源から出力された測定信号をマイクで収音することで伝達特性に応じたフィルタを生成するフィルタ生成方法であって、
   受聴者にマイクが装着された状態で取得された収音信号を第１の同期加算回数で同期加算することで、第１の同期加算信号を生成するステップと、
   前記受聴者以外にマイクが装着された状態で取得された収音信号を、前記第１の同期加算回数よりも多い第２の同期加算回数で同期加算することで、第２の同期加算信号を生成するステップと、
   前記第１の同期加算信号に対応する第１のスペクトルと、前記第２の同期加算信号に対応する第２のスペクトルを取得するように、前記第１及び第２の同期加算信号を周波数領域に変換するステップと、
   所定の周波数以下の帯域における第２のスペクトルを用いて、前記第１のスペクトルを補正して、第３のスペクトルを生成するステップと、
   前記第３のスペクトルを時間領域のデータに逆変換するステップと、を備えたフィルタ生成方法。
6. 前記所定の周波数よりも高い第１の周波数と、前記第１の周波数よりも低い第２の周波数との間の帯域に含まれるデータを用いて、前記第２のスペクトルのレベルを調整するステップをさらに備え、
   前記補正するステップでは、レベルが調整された前記第２のスペクトルを用いて、前記所定の周波数以下の帯域における前記第１のスペクトルを補正する請求項５に記載のフィルタ生成方法。
7. 前記補正するステップでは、前記所定の周波数以下の帯域における前記第１のスペクトルのデータを、前記第２のスペクトルのデータに置換する請求項５、又は６に記載のフィルタ生成方法。
8. ダミーヘッドに前記マイクを装着した状態で取得された前記収音信号を加算することで、前記第２の同期加算信号を生成する請求項５〜７のいずれか１項に記載のフィルタ生成方法。