頭部伝達関数の第 2 ピークが正中面上方の音像定位に及ぼす影響

一般社団法人 電子情報通信学会
THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS
信学技報
IEICE Technical Report
頭部伝達関数の第 2 ピークが正中面上方の音像定位に及ぼす影響
飯田 一博†
石井 要次‡
{†先進工学部, ‡大学院工学研究科 } 千葉工業大学
275-0016 千葉県習志野市津田沼 2-17-1
E-mail:
†kazuhiro.iida@it-chiba.ac.jp,
‡s1179501rq@s.chibakoudai.jp
あらまし 4kHz 付近の最も周波数の低いピーク(P1)と P1 より高い周波数で最も周波数の低い 2 つのノッチ(N1,
N2)だけで再構成されたパラメトリック頭部伝達関数(HRTF)は,前方と後方において実測 HRTF と同等の音像定位
精度を有することが報告されている.しかし,上方については,被験者によっては音像定位精度が低下する.本論
文では,上半球正中面において 2 つの音像定位実験を行った.実験結果は以下のとおりである.(1) N1N2P1 に P2
を加えることにより,90°および 120°で音像定位精度が向上した.(2) N1N2P1+P2 と実測 HRTF の音像定位誤差の差
は,目標とする 7 方向すべてにおいて 10°未満であった.(3) P2 だけを再現しても音像は上方に知覚されなかった.
また,HRTF を観察したところ,90°での N1 のレベルはノッチの検知限を超えないが,P2 を加えることにより N1
の相対レベルは検知可能となることが明らかになった.以上より,P2 それ自身はスペクトラルキューではないが,
P2 は N1 を強調することにより,正中面上方の音像定位に重要な役割を果たしていると言える.
Keyword 頭部伝達関数,音像定位,正中面,スペクトラルキュー,ピーク,ノッチ
Effects of the second lowest frequency spectral peak in head-related transfer
functions on the sound localization for the upper direction in the median plane
Kazuhiro IIDA† and
Yohji ISHII‡
{†Faculty of Advanced Engineering, ‡Graduate School } Chiba Institute of Technology
2-17-1 Tsudanuma, Narashino, Chiba, 275-0016 Japan
E-mail:
†kazuhiro.iida@it-chiba.ac.jp,
‡s1179501rq@s.chibakoudai.jp
Abstract It is reported that the parametric head-related transfer function (HRTF) recomposed of only the lowest frequency
peak around 4 kHz (P1) and the two lowest frequency notches (N1 and N2) above P1 provided approximately the same
localization performance as the measured HRTFs for the front and rear directions. However, for the upper directions, the
localization performance for some of the subjects decreased. In the present study, two localization tests were carried out in the
upper median plane. The results revealed the following: (1) the localization performance was improved by adding P2 to
N1N2P1 at the vertical angles of 90° and 120°; (2) the difference in localization error between N1N2P1+P2 and the measured
HRTF was less than 10° at all seven target vertical angles; and (3) a sound image was not perceived in the upward direction by
reproducing only P2. An observation of the HRTFs revealed that the level of N1 doesn’t exceed the threshold for the notch
detection at a vertical angle of 90°. However, by adding P2, the relative level of N1 becomes detectable. Therefore, P2 is not a
spectral cue, but, by enhancing N1, plays an important role in the localization for the upper directions in the median plane.
Keyword head-related transfer function, sound localization, the median plane, spectral cue, peak, notch
1. は じ め に
は周波数応答曲線における明確なノッチは,音源が前
ヒ ト の 頭 部 伝 達 関 数 (HRTF: head-related transfer
半球正中面の上方から下方に移動すると,低い周波数
function)の 5 kHz 以 上 の ノ ッ チ や ピ ー ク が 音 像 の 上 昇
に 移 動 す る こ と を 示 し た .Raykar et al. (2005)は 耳 介 に
角知覚に顕著に貢献していることは広く知られている.
よる深いノッチは仰角知覚に重要な特徴であり,ノッ
Shaw and Teranishi (1968)は 音 源 の 仰 角 が –45°か ら 45°
チ周波数は耳介の物理的な寸法や形状に関連すると報
に 変 化 す る と ノ ッ チ 周 波 数 が 6 kHz か ら 10 kHz に 変 化
告している.
す る こ と を 報 告 し て い る . Butler and Belendiuk (1977)
ノッチやピークの生成メカニズムも研究されてい
This article is a technical report without peer review, and its polished and/or extended version may be published elsewhere.
Copyright ©20●● by IEICE
る . Gardner and Gardner (1973)は 耳 介 の 窪 み を 順 次 埋
像定位性能に影響しないこと,つまり正中面定位に利
めて正中面音像定位実験を行い,耳介の閉塞を進める
用 さ れ て い る 特 性 は 3.8 か ら 16 kHz で あ る こ と を 示 し
と音像定位性能が低下すること,および低域成分と比
た.さらに彼らは,フィルタ処理なしと比較して,あ
較して高域成分がより重要であることを示した.
る 方 向 に 1.5 倍 以 上 の 回 答 が 得 ら れ た フ ィ ル タ を ス ペ
Musicant and Butler (1984)は , 4 kHz 低 域 通 過 ノ イ ズ と
クトラルキューと定義し,前方のキューは低域側のカ
比 較 し て , 4 kHz 高 域 通 過 ノ イ ズ は 有 意 に 重 要 で あ る
ッ ト オ フ 周 波 数 が 4 kHz か ら 8 kHz の 間 に あ る 1 オ ク
こ と , お よ び 4 kHz 高 域 通 過 ノ イ ズ に お い て は 耳 介 の
タ ー ブ ノ ッ チ と 13 kHz 以 上 の エ ネ ル ギ ー の 上 昇 ,上 方
閉塞により多くの前後誤判定が発生することを報告し
の キ ュ ー は 7 kHz か ら 9 kHz の 間 の 1/4 オ ク タ ー ブ の
た . Iida et al. (1998)は 耳 介 の 主 要 な 3 つ の 窪 み , す な
ピ ー ク , そ し て 後 方 の キ ュ ー は 10 kHz か ら 12 kHz の
わ ち 舟 状 窩 (scapha), 三 角 窩 (fossa), 耳 甲 介 腔 (concha)
間 に あ る 小 さ な ピ ー ク (こ の ピ ー ク の 高 域 側 お よ び 低
を 粘 土 で 閉 塞 し て 音 像 定 位 実 験 と HRTF 測 定 を 行 っ た
域 側 で の エ ネ ル ギ ー の 減 衰 を 伴 う )で あ る と 報 告 し た .
(粘 土 に は 外 耳 道 と 同 じ 直 径 の 穴 を 開 け ,外 耳 道 を 外 と
Iida et al. (2007)は 上 昇 角 音 像 定 位 を シ ミ ュ レ ー ト す
繋 げ て い る ). そ の 結 果 , 3 kHz 以 下 で は 全 て の 窪 み を
る パ ラ メ ト リ ッ ク HRTF モ デ ル を 提 案 し た . パ ラ メ ト
閉 塞 し て も ,そ の HRTF は 閉 塞 し て い な い 耳 介 の HRTF
リ ッ ク HRTF は , 音 源 の 上 昇 角 に 依 存 し な い 4 kHz の
と 同 等 で あ る が , 耳 甲 介 腔 を 閉 塞 す る だ け で 4 kHz 付
ピ ー ク (Shaw and Teranishi, 1968)を 下 限 周 波 数 と し て ,
近 の ピ ー ク と 4 kHz か ら 10 kHz の 間 の 鋭 い ノ ッ チ が 消
受 聴 者 の 実 測 HRTF か ら 抽 出 し た ノ ッ チ や ピ ー ク か ら
滅 し た . ま た , 耳 甲 介 腔 を 閉 塞 し た HRTF の 音 像 定 位
再構成される.ノッチとピークは周波数の順にラベル
精 度 は 閉 塞 し て い な い HRTF と 比 較 し て 有 意 に 低 下 し
が つ け ら れ る (P1, N1, P2, N2, ・ ・ ・ の よ う に ). ノ ッ
た.これらより,彼らは正中面音像定位のスペクトラ
チとピークは中心周波数,レベル,先鋭度でパラメト
ル キ ュ ー は 耳 甲 介 腔 の 伝 達 関 数 の 4 kHz 以 上 に 存 在 す
リックに表現される.彼らは上半球正中面において音
る と 結 論 付 け た . Takemoto et al. (2012) は FDTD
像 定 位 実 験 を 行 い 以 下 の こ と を 示 し た . (1) す べ て の
(finite-difference time-domain)法 を 用 い て 頭 部 形 状 か ら
ノ ッ チ と ピ ー ク で 再 構 成 さ れ た パ ラ メ ト リ ッ ク HRTF
HRTF を 計 算 し , 耳 介 が HRTF の 基 本 的 な ピ ー ク ・ ノ
は 被 験 者 本 人 の 実 測 HRTF と 同 等 の 定 位 性 能 を 有 す る .
ッチパターンを決定していると報告した.ピークは耳
(2) 4 kHz 付 近 の 第 1 ピ ー ク (P1)と ,P1 周 波 数 よ り 高 く
介のノーマルモードと一致する.ノッチ周波数では耳
最 も 周 波 数 の 低 い 2 つ の ノ ッ チ (N1, N2)だ け で 再 構 成
甲 介 舟 (cymba)と 三 角 窩 に 腹 が 生 じ ,耳 甲 介 腔 に 節 が 生
さ れ た パ ラ メ ト リ ッ ク HRTF は , 前 方 と 後 方 に お い て
じる.
は 被 験 者 本 人 の 実 測 HRTF と 同 等 の 定 位 性 能 を 有 す る .
ノッチとピークについては,微細構造よりむしろ概
(3) し か し 上 方 に お い て は , N1, N2, P1 で 再 構 成 さ れ
形 が 重 要 で あ る こ と が 報 告 さ れ て い る . Asano et al.
た パ ラ メ ト リ ッ ク HRTF の 音 像 定 位 性 能 は , 被 験 者 に
(1990)は ARMA モ デ ル で 平 滑 化 し た HRTF を ヘ ッ ド ホ
よ っ て は 本 人 の 実 測 HRTF と 比 較 し て 低 下 す る .(4) N1
ンで再生する正中面音像定位実験を行い,仰角判断の
と N2 の 周 波 数 は 音 源 の 上 昇 角 に 強 く 依 存 す る が , P1
主 要 な 手 掛 か り は 5 kHz 以 上 の 高 域 に あ り , ピ ー ク や
周波数はほぼ一定で上昇角に依存しない.これらの結
ノッチの微細部分ではなく巨視的なパターンの情報が
果 よ り ,彼 ら は N1 と N2 は ,少 な く と も 前 方 と 後 方 で
利 用 さ れ て い る こ と を 示 し た . Kulkarni and Colburn
は,スペクトラルキューとして重要な役割を果たして
(1998) は チ ュ ー ブ ホ ン (tube-phones)を 用 い て , 水 平 面
いると結論付けた.さらに彼らは,ヒトの聴覚システ
内 (方 位 角 : –45°, 0°, 45°, 180°)の 被 験 者 本 人 の HRTF
ム が 耳 入 力 信 号 か ら N1 と N2 を 分 析 す る た め の 参 照 情
の振幅スペクトルを 7 段階で平滑化して音像定位実験
報 と し て P1 を 用 い て い る 可 能 性 を 指 摘 し た .
を行い,微細なスペクトル構造は,ノッチやピークの
先 に 述 べ た よ う に ,Hebrank and Wright (1974)は 上 方
概形と比較すると,相対的に重要ではないことを示し
の ス ペ ク ト ラ ル キ ュ ー は 7 kHz か ら 9 kHz の 間 の 1/4
た.ただし,ノッチの 1 つが消滅するような過度な平
オクターブのピークであると主張した.このピークは
滑化を施すと,音像の仰角は上昇した.
P2 と 一 致 し , さ ら に Blauert (1969/70)が 提 案 し た 上 方
上昇角知覚のスペクトラルキューについても研究
の方向決定帯域とも一致する.しかし,狭帯域信号を
が 進 め ら れ て い る .Hebrank and Wright (1974)は 上 半 球
用いて得られた方向決定帯域が広帯域信号に対するス
正 中 面 (–30°か ら +210°ま で ,30°間 隔 )に 設 置 さ れ た ス ピ
ペクトラルキューの役割を果たすか否かは不明であ
ーカを用いて無響室で音像定位実験を行った.刺激は
る.
帯 域 幅 が low-pass, high-pass, band-pass, band-stop
本 論 文 は 2 つ の 目 的 を 有 す る . 1 つ は , N1,N2,P1 に
filters で 制 限 さ れ た 白 色 雑 音 で あ る .彼 ら は 16 kHz 以
P2 を 加 え , N1,N2,P1,P2 で 再 構 成 さ れ た パ ラ メ ト リ ッ
上 も し く は 3.8 kHz 以 下 の ス ペ ク ト ル エ ネ ル ギ ー は 音
ク HRTF が 上 方 に お い て 実 測 HRTF と 同 等 の 音 像 定 位
性能を有するか否かを検証することである.もう 1 つ
30°か ら 150°で は , 被 験 者 に よ っ て は N1,N2 は 浅 く
の 目 的 は , 正 中 面 上 方 の 音 像 定 位 に お け る P2 の 役 割
て不明瞭な場合があるので,以下のアルゴリズムによ
を解明することである.
り N1,N2 を 抽 出 し た .
(1) FFT に よ り 512 サ ン プ ル の HRIR の 振 幅 ス ペ ク ト
本論文では,方位角と仰角ではなく,側方角と上昇
角 で 表 現 し た 耳 軸 座 標 系 (Morimoto and Aokata, 1984)
ルを求め,差分法で極小値をすべて検出する.
(2) Iida and Ishii (2011)が 報 告 し た 以 下 の 回 帰 式 を 用
を用いる.上昇角は,音源を通る矢状面において,音
源と耳軸を結ぶ直線と水平線がなす角と定義する.
い て N1,N2 周 波 数 を 算 出 す る . 彼 ら は , 正 面 方 向 の
N1,N2 周 波 数 は 被 験 者 に 強 く 依 存 す る が , 上 昇 角 の 関
2. パ ラ メ ト リ ッ ク HRTF の 再 構 成 方 法
数 と し て の N1,N2 周 波 数 の 振 舞 い は 被 験 者 間 で 共 通 で
2.1. HRTF の 測 定
あるとみなせると報告している.
無 響 室 に お い て ,22–24 歳 の 3 人 の 男 性 被 験 者 (MKI,
OIS, OTK)お よ び 1 人 の 女 性 被 験 者 (YSD)の 上 半 球 正 中
面 7 方 向 (30°間 隔 )の HRTF を 測 定 し た .測 定 用 信 号 は
swept sine wave で ,サ ン プ リ ン グ 周 波 数 は 48 kHz で あ
る . 測 定 用 信 号 は 上 半 球 正 中 面 に 30°間 隔 で 設 置 し た
直 径 80mm の ス ピ ー カ (FOSTEX FE83E)か ら 提 示 し た .
ス ピ ー カ か ら 被 験 者 の 頭 部 中 心 ま で の 距 離 は 1.2m で
ある.被験者の外耳道入口で測定用信号を収録するた
め ,耳 栓 型 マ イ ク ロ ホ ン (Iida et al., 2014)を 用 い た .マ
10
𝑓𝑓/0
60
𝛽𝛽 = 1.001×1067 ×𝛽𝛽 8 − 6.431×106= ×𝛽𝛽 = + 8.686×
×𝛽𝛽 @ − 3.265×1060 ×𝛽𝛽 + 𝑓𝑓/0 0
[Hz]
(2)
𝑓𝑓/@ 𝛽𝛽 = 1.310×1067 ×𝛽𝛽 8 − 5.154×106= ×𝛽𝛽 = + 5.020×
1060 ×𝛽𝛽 @ + 2.563×10×𝛽𝛽 + 𝑓𝑓/@ 0
[Hz]
(3)
こ こ で 𝑓𝑓/0 と 𝑓𝑓/@ は そ れ ぞ れ N1,N2 周 波 数 を 表 し ,𝛽𝛽は 上
昇角である.
イクロホンの振動板は被験者の外耳道入口に位置し,
(3) (1)で 検 出 し た 極 小 値 の う ち , (2)で 算 出 し た N1,
blocked-entrances condition (Shaw and Teranishi, 1968)と
N2 周 波 数 か ら 0.2 オ ク タ ー ブ 以 内 で 最 も 深 い 極 小 値 を
み な せ る . HRTF は 以 下 の よ う に 求 め た .
𝐻𝐻𝐻𝐻𝐻𝐻𝐻𝐻%,' 𝜔𝜔 = 𝐺𝐺%,' 𝜔𝜔
𝐹𝐹(𝜔𝜔)
そ れ ぞ れ N1,N2 と す る .そ の 理 由 は ,上 昇 角 知 覚 に お
け る N1,N2 の 丁 度 可 知 差 (JND)は 0.1 か ら 0.2 オ ク タ ー
(1)
ブ と 考 え ら れ る (Iida and Ishii, 2011)か ら で あ る .
P1,P2 周 波 数 は 上 昇 角 に 依 存 し な い の で , 30° か ら
150°の P1,P2 は 0°と 同 様 と し た .
こ こ で , 𝐹𝐹(𝜔𝜔)は 被 験 者 の い な い 自 由 空 間 に お い て 被 験
2.3. パ ラ メ ト リ ッ ク HRTF の 生 成
て被験者の外耳道入口で耳栓型マイクロホンにより測
成したノッチとピークを重ね合わせることにより生成
定したインパルス応答のフーリエ変換である.
し た . 実 測 HRTF, N1,N2,P1 で 再 構 成 し た パ ラ メ ト リ
者の頭部中心に相当する位置で測定したインパルス応
答 𝑓𝑓(𝑡𝑡)の フ ー リ エ 変 換 で あ り ,𝐺𝐺%,' 𝜔𝜔 は 自 由 空 間 に お い
2.2. ノ ッ チ と ピ ー ク の 抽 出
N1, N2, P1, P2 は 耳 介 で 生 成 さ れ て い る た め (Shaw
パ ラ メ ト リ ッ ク HRTF は , 2 次 の IIR フ ィ ル タ で 作
ッ ク HRTF(以 降 N1N2P1 と 呼 ぶ ), N1,N2,P1,P2 で 再 構
成 し た パ ラ メ ト リ ッ ク HRTF(以 降 N1N2P1+P2 と 呼 ぶ )
and Teranishi, 1968; Lopez-Poveda and Meddis, 1996;
の 例 を FIG. 1 に 示 す .パ ラ メ ト リ ッ ク HRTF は 再 構 成 し
Takemoto et al., 2012),0°と 180°に つ い て は 被 験 者 の 頭
たノッチとピークを精度よく再現し,その他の帯域の
部 イ ン パ ル ス 応 答 (HRIR)の 初 期 部 分 か ら 抽 出 し た .そ
スペクトルは平坦である.
の ア ル ゴ リ ズ ム (Iida et al., 2014) を 以 下 に 示 す .
20
(1) HRIR の 絶 対 値 が 最 大 と な る サ ン プ ル を 検 出 す る .
(3) 全 て の 要 素 を ゼ ロ と し た 512 ポ イ ン ト の 配 列 を
用 意 し , (2) で 切 り 出 し た HRIR を 上 書 き す る . た だ
し , 最 大 サ ン プ ル を 配 列 の 257 番 目 に 合 わ せ る
Amplitude [dB]
(2) 検 出 し た サ ン プ ル を 中 心 と し て 4 次 96 サ ン プ ル
の Blackman-Harris 窓 で HRIR を 切 り 出 す .
Measured
P1
10
N1N2P1
P2
0
-10
N2
(4) FFT に よ り 512 ポ イ ン ト の 配 列 の 振 幅 ス ペ ク ト
ルを求め,差分法で極大値および極小値をすべて検出
する.
N1N2P1+P2
-20
N1
0
4
8
12
Frequency [kHz]
16
20
(5) 3 kHz 以 上 で 最 も 周 波 数 の 低 い 2 つ の 極 大 値 を そ
れ ぞ れ P1,P2 と し ,P1 よ り 高 い 帯 域 で 最 も 周 波 数 の 低
い 2 つ の 極 小 値 を N1,N2 と す る .
FIG. 1. Examples of a measured HRTF and parametric HRTFs
(N1N2P1 and N1N2P1+P2) for a vertical angle of 90°.
ホ ン に 覆 わ れ て い な い (FIG. 2(a))). ヘ ッ ド ホ ン か ら 耳
3. 音 像 定 位 実 験 1
N1N2P1+P2 が 上 方 に お い て 実 測 HRTF と 同 等 の 音 像
栓 型 マ イ ク ロ ホ ン ま で の 200 Hz か ら 17 kHz ま で の 代
定位精度を有するか否かを検証するため, 上半球正中
表 的 な peak-to-peak レ ン ジ は 約 20 dB で あ っ た が , 補
面において音像定位実験を行った.
正 フ ィ ル タ G に よ り 約 3 dB ま で 低 減 し た (FIG. 3).
DT990 で は ,耳 介 が ヘ ッ ド ホ ン に 覆 わ れ て い る た め
3.1. 実 験 方 法
被 験 者 は 2.1 節 で HRTF の 測 定 を 行 っ た 4 人 (MKI,
OIS, OTK, YSD)で あ る . HRTF を 以 下 の 3 種 類 を 用 い
(FIG. 2(b)), 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃を 測 定 す る に は ヘ ッ ド ホ ン の 着 脱
が 必 要 と な る .𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃は ヘ ッ ド ホ ン の 装 着 位 置 に よ り
た . (1) 被 験 者 本 人 の 実 測 HRTF, (2) 被 験 者 本 人 の
変化するので,それを測定してもヘッドホンの着脱を
N1N2P1, (3) 被 験 者 本 人 の N1N2P1+P2.
行 う と 正 確 な 補 正 は 保 証 さ れ な い .以 上 の 理 由 に よ り ,
音 像 定 位 実 験 は 防 音 室 で 行 っ た . 室 の 内 法 は 4.6m
(W) × 5.8m (D) × 2.8m (H)で , 暗 騒 音 は 19.5 dB(A)
で あ っ た . 実 験 シ ス テ ム の 構 成 は , ノ ー ト PC (DELL
XPS M1330),オ ー デ ィ オ イ ン タ フ ェ イ ス (RME Fireface
DT990 で は 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の 補 正 は 行 わ な か っ た .参 考 の た め ,
FIG. 4 に 4 人 の 被 験 者 の DT990 の 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の 測 定 例 を 示
す . peak-to-peak レ ン ジ は , Møller et al. (1995)の 報 告
と 同 様 30 dB を 超 え る .
400) , ア ン プ (Marantz PM4001) , A/D コ ン バ ー タ
(Roland M-10MX), オ ー プ ン タ イ プ ヘ ッ ド ホ ン , 耳 栓
(a)
(b)
型マイクロホンである.
外耳道入口を開放した状態における鼓膜での音圧 P
は,閉塞した状態の外耳道入口の音圧に以下に示す補
正フィルタ G を施してヘッドホンから提示することに
FIG. 2. Headphones used in the experiments. (a): K1000 (AKG),
よ り 得 ら れ る (Møller et al., 1995).
𝐺𝐺 =
≜
0
E×FGH
0
E×FGH
×
IJKLNKOKP QIRJKSTRUOJ
IJKLNKOKP QILKSVKWVUO
×𝑃𝑃𝑃𝑃𝑃𝑃,
(4)
10
(5)
0
(6)
ここで,S は音源,M は耳栓型マイクロホンの伝達関
Relativelevel[dB]
𝑃𝑃 = 𝑆𝑆×𝐻𝐻𝐻𝐻𝐻𝐻𝐻𝐻×𝐺𝐺,
(b): DT990 (Beyerdynamic).
-10
-20
数 ,𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃は 閉 塞 し た 外 耳 道 入 口 で 測 定 し た ヘ ッ ド ホ
-30
ン の 伝 達 関 数 , 𝑍𝑍\]'^]_]` と 𝑍𝑍a\]bcad_\ は そ れ ぞ れ 外 耳 道
-40
お よ び ヘ ッ ド ホ ン の イ ン ピ ー ダ ン ス , 𝑍𝑍']be]fed_ は 外 耳
道から外側をみた放射インピーダンスである.G の第
2 項 は PDR (pressure division ratio)と 呼 ば れ , Møller et
al. は PDR が 1 と な る ヘ ッ ド ホ ン を FEC(free-air
0
4000
8000
12000
16000
20000
24000
Frequency [Hz]
FIG. 3. An example of transfer functions between K1000 and the
earplug-type microphone with a compensation filter, G.
equivalent coupling to the ear)ヘ ッ ド ホ ン と 定 義 し た .
0
2 種 類 の ヘ ッ ド ホ ン , K1000 (AKG)お よ び DT990 PRO
-10
(Beyerdynamic)(以 降 DT990 と 記 述 )を 用 い た .
-20
K1000 で は , 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の 補 正 は 200 Hz か ら 17 kHz の
範 囲 で 行 っ た . 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃 は 以 下 の 手 順 で 測 定 し た . 被
験者は防音室の中央に座り,耳栓型マイクロホンを外
耳 道 に 装 着 し た . マ イ ク ロ ホ ン の 振 動 板 は 2.1 節 の
HRTF 測 定 と 同 様 に 外 耳 道 入 口 に 位 置 し た . 被 験 者 に
ヘ ッ ド ホ ン を 装 着 し て M 系 列 信 号 (48kHz サ ン プ リ ン
グ , 12 次 , 繰 返 し な し )を 提 示 し , 耳 栓 型 マ イ ク ロ ホ
ン ま で の 伝 達 関 数 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃を 測 定 し た . さ ら に , ヘ ッ ド
ホンを装着したまま,その位置を変えないようにして
耳 栓 型 マ イ ク ロ ホ ン だ け を 取 り 外 し た (耳 介 は ヘ ッ ド
MPTF[dB]
本 実 験 で は ,彼 ら が FEC ヘ ッ ド ホ ン と み な せ る と し た
MKI
OIS
OTK
YSD
-30
-40
-50
-60
0
4000
8000
12000
16000
20000
24000
Frequency [Hz]
FIG. 4. Transfer functions between DT990 and the earplug-type
microphone, 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃.
K1000 と DT990 の 再 生 条 件 の 違 い を 表 現 す る た め に ,
実 測 HRTF と 比 べ て や や 後 方 に 分 布 す る 傾 向 が あ っ た .
以 降 ,K1000 を 用 い て 式 (6)が 1 と み な せ る 条 件 を FEC
N1N2P1+P2 の 回 答 は ,0°, 30°, 60°, 180°で は 実 測 HRTF
条 件 ,DT990 を 用 い て 1/(𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃)が 残 る 条 件 を 準 FEC
お よ び N1N2P1 と 同 様 で あ っ た . 90°, 120°で は 回 答 分
条 件 と 表 す .Møller et al. (1995)が 定 義 し た FEC ヘ ッ ド
布 は N1N2P1 よ り 実 測 HRTF に 近 づ い た .し か し ,150°
ホ ン と の 違 い に 注 意 さ れ た い .準 FEC 条 件 は 理 論 的 厳
では上方寄りに知覚した.
密さに欠けるが,伝達関数の測定が不要であるため,
被 験 者 OTK で は ,実 測 HRTF に 対 し て は 0°, 30°, 90°,
一般の受聴者が使いやすいという利点を有している.
180°で は 目 標 方 向 に 回 答 し た .し か し ,60°, 120°, 150°
つ ま り , 準 FEC 条 件 は 実 用 性 の 観 点 か ら 検 証 し た .
で は 90°付 近 に 回 答 し た . N1N2P1 の 回 答 は 0°, 150°,
音 源 信 号 は 200Hz か ら 17kHz の 広 帯 域 ホ ワ イ ト ノ
180°で は 実 測 HRTF と 同 様 で あ っ た .し か し ,30°で は
イ ズ で あ る . 刺 激 は 被 験 者 の 左 右 の 外 耳 道 入 口 で 63
前 方 に 寄 る 傾 向 が み ら れ , 60°, 90°, 120°で は 0°と 90°
dB SPL で 提 示 し た . 刺 激 の 提 示 時 間 は 1.2 秒 (立 上 が
付 近 の 両 方 に 回 答 が 分 布 し た . N1N2P1+P2 の 回 答 は ,
り ,立 下 が り そ れ ぞ れ 0.1 秒 を 含 む )で あ る .目 標 方 向
す べ て の 目 標 方 向 に お い て 実 測 HRTF と 同 様 で あ っ た .
は 上 半 球 正 中 面 に お け る 30°間 隔 の 7 方 向 で あ る . 被
N1N2P1 で 観 測 さ れ た 回 答 の 分 離 は み ら れ な か っ た .
験者が目標方向を推測することを避けるため,回答方
被 験 者 YSD で は , 実 測 HRTF に 対 す る 回 答 は 逆 S
法 に は マ ッ ピ ン グ 法 を 用 い た . ノ ー ト PC の 画 面 上 に
字 カ ー ブ を 描 き , 30°, 60°, 120°で や や 上 方 に 寄 る 傾 向
正中面を示す円周と矢印を表示した.被験者のタスク
が あ っ た . N1N2P1 の 回 答 は , 90°, 120°, 180°で は 実 測
は,知覚した上昇角を円周上にスタイラスペンでクリ
HRTF と 同 様 で あ っ た . し か し , 0°, 30°, 60°で は 実 測
ックすることである.頭内定位した場合はチェックボ
HRTF と 比 べ て や や 上 方 に 分 布 す る 傾 向 が あ っ た .150°
ックスに印を入れるよう指示した.
で は 目 標 方 向 付 近 に 回 答 し た . N1N2P1+P2 の 回 答 は
実 測 HRTF と 2 種 類 の パ ラ メ ト リ ッ ク HRTF の 音 像
定 位 実 験 は 別 々 に 行 っ た .Møller et al. (1995)は 2 つ の
60°を 除 き N1N2P1 と 同 様 で あ っ た . 60°で は 回 答 の 分
散 は N1N2P1 よ り 大 き い .
条 件 ,す な わ ち (1) 1 人 の 被 験 者 の HRTF を ラ ン ダ ム 化
し た 条 件 と (2) 数 人 の HRTF を ラ ン ダ ム 化 し た 条 件 を
比較し,同じ刺激群を分離して実験しても混合して実
Subj:
MKI
験しても有意な差はないことを示している.
1 回 の 試 行 で 35 刺 激 (7 方 向 ×5 回 )を ラ ン ダ ム 化 し て
被 験 者 に 提 示 し た .1 試 行 の 所 要 時 間 は 約 7 分 で あ る .
3 種 類 の HRTF に つ い て 2 試 行 ず つ 実 験 を 行 っ た . つ
OIS
ま り , 被 験 者 は 各 刺 激 に 対 し て 10 回 ず つ 回 答 し た .
3.2. 実 験 結 果
3.2.1. FEC 条 件
A. 回 答 分 布
FIG. 5 に 実 測 HRTF,N1N2P1,N1N2P1+P2 に 対 す る
OTK
4 人の被験者の回答を示す.縦軸は知覚した上昇角,
横軸は目標上昇角である.円の直径は回答の頻度に比
被 験 者 MKI で は ,実 測 HRTF に 対 す る 回 答 は ほ ぼ 対
角 線 上 に 分 布 し た .N1N2P1 の 回 答 は 0°, 30°, 150°, 180°
で は 実 測 HRTF と 同 様 で あ っ た . し か し , 他 の 3 方 向
(60°, 90°, 120°)で は 回 答 の 分 散 は 実 測 HRTF よ り 大 き
い . N1N2P1+P2 の 回 答 は , 0°, 30°, 150°, 180°で は 実 測
HRTF お よ び N1N2P1 と 同 様 で あ っ た . 90°, 120°で は
YSD
Perceivedrisingangle[deg.]
例 し , 回 答 方 向 は 5°間 隔 で 丸 め て い る .
Targetrisingangle[deg.]
Measured HRTF
N1N2P1
N1N2P1+P2
分 散 は 実 測 HRTF よ り や や 大 き く , N1N2P1 よ り や や
小 さ い . し か し , 60°で は 分 散 は N1N2P1 よ り 大 き い .
被 験 者 OIS で は ,実 測 HRTF に 対 す る 回 答 は 逆 S 字
FIG. 5. Responses to the measured HRTF, N1N2P1, and
N1N2P1+P2 under the FEC condition. The ordinate represents
カ ー ブ を 描 き , 30°, 60°, 120°で や や 上 方 寄 り に 回 答 す
the perceived vertical angle, and the abscissa represents the
る 傾 向 が あ っ た . N1N2P1 の 回 答 は , 0°, 30°, 60°, 180°
target vertical angle. The diameter of each circle is proportional
で は 実 測 HRTF と 同 様 で あ っ た が ,90°, 120°, 150°で は
to the number of responses with a resolution of 5°.
B. 平 均 上 昇 角 定 位 誤 差
C. 前 後 誤 判 定 率
各 HRTF と 各 目 標 方 向 の 平 均 上 昇 角 定 位 誤 差 (𝐸𝐸𝐸𝐸)を
TABLE I に 前 後 誤 判 定 率 を 示 す . こ こ で , 前 後 誤 判
求 め た (FIG. 6).こ こ で ,平 均 上 昇 角 定 位 誤 差 は 目 標 の
定率は上半球正中面において目標方向と異なる四分円
上昇角と回答した上昇角の差の絶対値の平均と定義す
に音像を知覚した率と定義する.
る . さ ら に パ ラ メ ト リ ッ ク HRTF と 実 測 HRTF の 𝐸𝐸𝐸𝐸の
0°と 180°で は , い ず れ の HRTF に お い て も 前 後 誤 判
差 𝛥𝛥𝛥𝛥𝛥𝛥を 求 め た (FIG. 7).
定 率 は 0%で あ っ た .し か し ,他 の 4 方 向 に お い て は ,
(2010)の 報 告 と 同 様 , 水 平 面 上 の 0°, 180°で は 小 さ く
そ れ よ り 高 く , N1N2P1+P2 の 前 後 誤 判 定 率 は N1N2P1
(3.1°, 1.8°),水 平 面 か ら 上 昇 し た 30°か ら 150°で は 大 き
よりやや低い.
実 測 HRTF の 𝐸𝐸𝐸𝐸は ,Carlile et al. (1997)や Majdak et al.
N1N2P1 と N1N2P1+P2 の 前 後 誤 判 定 率 は 実 測 HRTF の
い (18.1-29.1°)傾 向 が あ る .N1N2P1, N1N2P1+P2 で も 同
HRTF 間 の 前 後 誤 判 定 率 に 統 計 的 有 意 な 差 が あ る か
様 の 角 度 依 存 性 が み ら れ た .N1N2P1, N1N2P1+P2 の 𝐸𝐸𝐸𝐸
否かを明らかにするためにカイ 2 乗検定を行った.そ
は 実 測 HRTF よ り 大 き い が ,N1N2P1+P2 の 𝐸𝐸𝐸𝐸は 0°, 30°,
の 結 果 , い ず れ の 目 標 方 向 に お い て も HRTF 間 に 有 意
れ は 被 験 者 OIS が 上 方 寄 り に 回 答 し た こ と に 依 る が ,
TABLE I Ratio of front-back confusion for each HRTF and target
90°, 120°で N1N2P1 よ り 小 さ い . 特 に 90°で は 10.8°減
な差は認められなかった.
少 し た . た だ し , 150°で は N1N2P1 よ り 増 加 し た . こ
vertical angle under the FEC condition.
上方に寄った原因は不明である.
N1N2P1 と 実 測 HRTF の 𝐸𝐸𝐸𝐸の 差 (𝛥𝛥𝛥𝛥𝛥𝛥)は 150°,180°で は
小 さ く (0.6°, 0.4°) , 他 の 上 昇 角 に お い て は 8.1° か ら
14.0°で あ っ た . N1N2P1+P2 の 𝛥𝛥𝛥𝛥𝛥𝛥は 90°, 120°, 180°に
お い て 小 さ く (–0.4°, 0.4°, 0.4°),そ れ 以 外 の 上 昇 角 で も
2������������������������
.12-
�
��
��
��
���
���
���
�.12-
����
����
����
�
����
����
����
������
����
����
����
�
����
����
����
���������
����
����
����
�
����
����
����
5.9°か ら 8.9°で あ り ,す べ て の 上 昇 角 で 10°未 満 で あ っ
D. 頭 内 定 位 率
た.
被 験 者 YSD が 60°の 実 測 HRTF に 対 し て 1 度 だ け 頭
Measured
Meanvertical
localizationerror,Ev [deg.]
60
N1N2P1
内 に 音 像 を 知 覚 し た こ と を 除 き , い ず れ の HRTF に お
N1N2P1+P2
50
いても頭内定位は生じなかった.
40
3.2.2. 準 FEC 条 件
A. 回 答 分 布
30
FIG. 8 に 実 測 HRTF,N1N2P1,N1N2P1+P2 に 対 す る
20
4 人の被験者の回答を示す.
10
0
被 験 者 MKI で は ,実 測 HRTF に 対 す る 回 答 は ほ ぼ 対
0
30
60
90
120
150
180
Ave.
Targetverticalangle[deg.]
角 線 上 に 分 布 し た .N1N2P1 の 回 答 は ,90°と 120°以 外
で は 実 測 HRTF と 同 様 で あ っ た . 90°と 120°で は 回 答
FIG. 6. Mean vertical localization error, 𝐸𝐸𝐸𝐸, for each HRTF and
target vertical angle under the FEC condition.
の 分 散 は 実 測 HRTF よ り 大 き い .N1N2P1+P2 の 回 答 は
120°を 除 き 実 測 HRTF と 同 様 で あ っ た . 120°で は 上 方
寄 り に 回 答 す る 傾 向 が あ っ た .90°で は ,回 答 の 分 散 は
N1N2P1 よ り 小 さ く , 実 測 HRTF と 同 等 で あ っ た .
Differenceinmeanvertical
localizationerror,Δ Ev [deg.]
30
N1N2P1
被 験 者 OIS で は ,実 測 HRTF に 対 す る 回 答 は 逆 S 字
N1N2P1+P2
カ ー ブ を 描 き , 30°, 60°で は や や 上 方 に , 90°で は や や
20
後 方 に 回 答 す る 傾 向 が あ っ た . N1N2P1 の 回 答 は , 0°,
180°で は 実 測 HRTF と 同 様 で あ っ た .し か し , 90°, 120°
10
0
-10
で は や や 後 方 に 回 答 し , 30°, 150°で は 回 答 の 分 散 は 大
き い .N1N2P1+P2 の 回 答 は ,90°, 150°を 除 き 実 測 HRTF
0
30
60
90
120
150
180
Ave.
Targetverticalangle[deg.]
FIG. 7. Difference in 𝐸𝐸𝐸𝐸 between the parametric HRTF and the
measured HRTF, 𝛥𝛥𝛥𝛥𝛥𝛥, the under FEC condition.
と 同 様 で あ っ た . 90°, 150°で は N1N2P1 と 同 様 の 分 布
であった.
被 験 者 OTK で は , 実 測 HRTF に 対 し て は 0°, 180°で
は 目 標 方 向 に 回 答 し た が , 30°で は 前 方 か ら 後 方 , 60°
か ら 150° で は 上 方 か ら 後 方 の 広 い 範 囲 に 回 答 が 分 布
し た . N1N2P1 の 回 答 は , 0°, 180°で は 実 測 HRTF と 同
様 で あ っ た . し か し 30°で は 回 答 の 分 散 は 実 測 HRTF
特 に 60°, 90°, 120°で は , N1N2P1 と 比 較 し て そ れ ぞ れ
と 比 べ て 大 き い . 60°か ら 150°で は 後 方 に 回 答 し た .
13.4°, 17.5°, 11.8°減 少 し た .
N1N2P1+P2 の 回 答 は 30°を 除 き 実 測 と 同 様 で あ っ た .
30°で は 回 答 の 分 散 は 実 測 HRTF よ り も 大 き い .
被 験 者 YSD で は , 実 測 HRTF に 対 す る 回 答 は 逆 S
字 カ ー ブ を 描 き , 30°, 60°で は や や 上 方 に , 90°で は や
や 後 方 に 回 答 す る 傾 向 が あ っ た .N1N2P1 の 回 答 は 60°,
N1N2P1 と 実 測 HRTF の 𝐸𝐸𝐸𝐸の 差 (𝛥𝛥𝛥𝛥𝛥𝛥)は 0°, 150°, 180°
で は 小 さ く (3.2°, 2.7°, 2.2°) , 他 の 上 昇 角 に お い て は
9.9°か ら 23.5°で あ っ た . N1N2P1+P2 の 𝛥𝛥𝛥𝛥𝛥𝛥は 30°で は
19.6°で あ っ た が , 他 の 方 向 で は 10°未 満 で あ っ た .
90°を 除 き 実 測 HRTF と 同 様 で あ っ た . 60°, 90°で は 回
は す べ て の 目 標 方 向 で N1N2P1 と 同 様 で あ っ た . た だ
し ,90°の 回 答 の 分 散 は や や 狭 ま り ,150°の 回 答 の 分 散
はやや広がった.
Measured
60
Meanvertical
localizationerror,Ev [deg.]
答 は 上 方 か ら 後 方 に 広 く 分 布 し た .N1N2P1+P2 の 回 答
Subj:
MKI
N1N2P1
N1N2P1+P2
50
40
30
20
10
0
0
30
60
90
120
150
180
Ave.
Targetverticalangle[deg.]
FIG. 9. Mean vertical localization error, 𝐸𝐸𝐸𝐸, for each HRTF and
target vertical angle under the quasi-FEC condition.
OIS
Differenceinmeanvertical
localizationerror, Δ Ev [deg.]
30
OTK
N1N2P1
20
10
0
0
30
Perceivedrisingangle[deg.]
-10
YSD
N1N2P1+P2
60
90
120
150
180
Ave.
Targetverticalangle[deg.]
FIG. 10. Difference in 𝐸𝐸𝐸𝐸 between the parametric HRTF and the
measured HRTF, 𝛥𝛥𝛥𝛥𝛥𝛥, under the quasi-FEC condition.
C. 前 後 誤 判 定 率
Targetrisingangle[deg.]
Measured HRTF
N1N2P1
N1N2P1+P2
FIG. 8. Responses to the measured HRTF, N1N2P1, and
N1N2P1+P2 under the quasi-FEC condition.
TABLE II に 前 後 誤 判 定 率 を 示 す .0°と 180°で は ,い
ず れ の HRTF に お い て も 前 後 誤 判 定 率 は 0%で あ っ た .
し か し ,他 の 4 方 向 に お い て は ,N1N2P1 と N1N2P1+P2
の 前 後 誤 判 定 率 は 実 測 HRTF の そ れ よ り 高 い .
TABLE III に カ イ 2 乗 検 定 の 結 果 を 示 す . 120°に お
B. 平 均 上 昇 角 定 位 誤 差
各 HRTF と 各 目 標 方 向 の 平 均 上 昇 角 定 位 誤 差 (𝐸𝐸𝐸𝐸)お
よ び , パ ラ メ ト リ ッ ク HRTF と 実 測 HRTF の 𝐸𝐸𝐸𝐸の 差
い て 実 測 HRTF と N1N2P1 の 間 に , 30°に お い て 実 測
HRTF と N1N2P1+P2 の 間 に 統 計 的 有 意 な 差 (p < 0.05)
が認められた.
(𝛥𝛥𝛥𝛥𝛥𝛥)を そ れ ぞ れ 求 め た (FIGs. 9, 10).
実 測 HRTF の 𝐸𝐸𝐸𝐸は ,FEC 条 件 と 同 様 に 水 平 面 上 の 0°,
180°で は 小 さ く (4.7°, 2.1°), 水 平 面 か ら 上 昇 し た 30°
TABLE II Ratio of front-back confusion for each subject, HRTF,
and target vertical angle under the quasi-FEC condition.
か ら 150°で は 大 き い (16.2°-37.6°)傾 向 が あ る .N1N2P1,
N1N2P1+P2 で も 同 様 の 角 度 依 存 性 が み ら れ た .し か し ,
こ れ ら の 𝐸𝐸𝐸𝐸は 実 測 HRTF よ り 大 き い .N1N2P1+P2 の 𝐸𝐸𝐸𝐸
は 30°を 除 く す べ て の 目 標 方 向 で N1N2P1 よ り 小 さ い .
2������������������������
.12-
�
��
��
��
���
���
���
�.12-
����
����
����
�
����
����
����
������
����
����
����
�
����
����
����
���������
����
����
����
�
����
����
����
TABLE III Results of chi-square tests for the ratio of front-back
confusion
under the FEC condition. *: p < 0.05, **: p < 0.01.
"$ "&
!"%
!# !# 同 様 に 分 布 し た .す な わ ち ,MKI と OIS は 前 方 か ら 後
方 ,OTK は 上 方 か ら 後 方 ,YSD は 斜 め 後 ろ か ら 後 方 に
広く分布した.いずれの被験者においても,上方に回
一 方 , P1, P2, P1P2 に 対 す る 回 答 は 被 験 者 内 で ほ ぼ
答が集中することはなかった.
(a)MKI
D. 頭 内 定 位 率
3 人 の 被 験 者 (MKI, OIS, OTK)で は ,い ず れ の HRTF
においても頭内定位は生じなかった.しかし被験者
YSD で は ,N1N2P1 (30°, 90°で 1 回 ,60°で 3 回 )お よ び
N1N2P1+P2 (30°, 60°, 90°, 120°で 1 回 )で 頭 内 定 位 が 生
じた.
P1
4. 音 像 定 位 実 験 2
P2Measured N1N2P1+P2
P1P2
N1N2P1
(b)OIS
音 像 定 位 実 験 1 で は ,N1N2P1 に P2 を 加 え る こ と に
より正中面上方の音像定位精度が向上することを示し
た .音 像 定 位 実 験 2 の 目 的 は ,P2 が 正 中 面 上 方 の 音 像
定位で担う役割を解明することである.
4.1. 実 験 方 法
実験方法は音像定位実験 1 と同様である.ただし,
被 験 者 本 人 の 90°の (1) P1, (2) P2, (3) P1P2 で 再 構 成 し
P1
た 3 種 類 の パ ラ メ ト リ ッ ク HRTF を 用 い た .
P2Measured N1N2P1+P2
P1P2
N1N2P1
(c) OTK
FIG. 11 に 実 測 HRTF と P1, P2, P1P2 で 再 構 成 し た パ
ラ メ ト リ ッ ク HRTF の 例 を 示 す .パ ラ メ ト リ ッ ク HRTF
は再構成したピークを精度よく再現し,その他の帯域
のスペクトルは平坦である.
20
Amplitude [dB]
Measured
P1
P2
P1P2
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
10
(d)YSD
0
-10
-20
0
4
8
12
Frequency [kHz]
16
20
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
FIG. 11. Examples of a measured HRTF and parametric HRTFs
FIG. 12. Responses to P1, P2, and P1P2 under the FEC condition.
(P1, P2, and P1P2) for a vertical angle of 90°.
For comparison, the responses to the measured HRTF, N1N2P1,
and N1N2P1+P2 for a target vertical angle of 90° obtained in
4.2. 実 験 結 果
4.2.1. 回 答 分 布
FEC 条 件 で の P1, P2, P1P2 に 対 す る 回 答 分 布 を FIG.
12 に 示 す .音 像 定 位 実 験 1 で 求 め た 目 標 方 向 90°に お
け る 実 測 HRTF, N1N2P1, N1N2P1+P2 の 回 答 分 布 も 併
せ て 示 す . N1N2P1+P2 の 回 答 分 布 は , 3.2.1.項 で 示 し
た よ う に MKI で や や 分 散 が 大 き く な る も の の , 実 測
HRTF の そ れ と 同 様 で あ る .
localization test 1 are also shown.
準 FEC 条 件 で の 回 答 分 布 を FIG. 13 に 示 す . 各 被 験
者 の 回 答 分 布 を Fig.13 に 示 す .N1N2P1+P2 で は ,3.2.2.
項 で 示 し た よ う に , OIS で や や 分 散 が 大 き く な る も の
の , 実 測 HRTF と 同 様 の 回 答 分 布 で あ っ た .
一 方 , P1, P2, P1P2 に 対 す る 回 答 は , FEC 条 件 と 同
様 に ,被 験 者 内 で ほ ぼ 同 様 に 分 布 し た .MKI は 前 方 か
ら 後 方 ,OIS は 前 方 か ら 斜 め 後 ろ に 広 く 分 布 し た .OTK
と YSD は 斜 め 後 ろ か ら 後 方 に 分 布 し た .い ず れ の 被 験
者 に お い て も ,上 方 に 回 答 が 集 中 す る こ と は な か っ た .
𝐸𝐸𝐸𝐸(52.9°, 51.6°, 54.6°)は 実 測 HRTF と 比 較 し て 顕 著 に
大きい値であった.
Meanrisingangle
localizationerror,Ev [deg.]
60
(a)MKI
(a)FECcondition
50
40
30
20
10
0
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
P1
P2
P1+2
60
Meanrisingangle
localizationerror,Ev [deg.]
(b)OIS
Measured N1N2P1 N1N2P1+2
(b)Quasi-FEC condition
50
40
30
20
10
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
0
P1
P2
P1+2
Measured N1N2P1 N1N2P1+2
FIG. 14. Mean vertical localization error, Ev, under (a) the FEC
(c) OTK
condition and (b) the quasi-FEC condition at a target vertical
angle of 90°.
4.2.3. 頭 内 定 位 率
P1, P2, P1P2 の い ず れ に お い て も , ど の 被 験 者 も 頭
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
(d)YSD
内定位は生じなかった.
5. 考 察
5.1. 従 来 の 研 究 と の 比 較
先 に 述 べ た よ う に ,Hebrank and Wright (1974)は フ ィ
ル タ 処 理 な し と 比 較 し て , あ る 方 向 に 1.5 倍 以 上 の 回
答が得られたフィルタをスペクトラルキューと定義し
た . 彼 ら は , 上 方 の キ ュ ー は 7 kHz か ら 9 kHz の 間 の
P1
P2Measured N1N2P1+P2
P1P2
N1N2P1
FIG. 13. Responses to P1, P2, and P1P2 under the quasi-FEC
condition.
4.2.2. 平 均 上 昇 角 定 位 誤 差
平 均 上 昇 角 定 位 誤 差 (𝐸𝐸𝐸𝐸)を Fig.14 に 示 す .
FEC 条 件 (a)に お い て は ,N1N2P1+P2 の 𝐸𝐸𝐸𝐸(17.7°)は 実
測 HRTF(18.1°)と 同 等 で あ っ た . し か し , P1, P2, P1P2
の 𝐸𝐸𝐸𝐸(56.7°, 48.8°, 42.1°)は 実 測 HRTF と 比 較 し て 顕 著
に大きい値であった.
準 FEC 条 件 (b)に お い て は ,実 測 HRTF と N1N2P1+P2
の 𝐸𝐸𝐸𝐸は そ れ ぞ れ 27.9°, 34.4°で あ っ た .P1, P2, P1+P2 の
1/4 オ ク タ ー ブ の ピ ー ク で あ る と 主 張 し た . こ の ピ ー
ク は P2 と 一 致 し ,さ ら に Blauert (1969/70)が 提 案 し た
上方の方向決定位帯域とも一致する.
し か し , 第 3,4 章 の 結 果 は , 広 帯 域 信 号 に 対 し て は
P2 を 再 現 し た だ け で は 上 方 に 音 像 定 位 し な い こ と を
明 確 に 示 し て い る . つ ま り , Hebrank and Wright の 主
張するピークは上方のスペクトラルキューではない.
5.2. 正 中 面 上 方 に お け る P2 の 役 割
FIG. 15 は , FEC 条 件 の 目 標 方 向 90° に お い て ,
N1N2P1 で 上 方 と 前 方 に 回 答 し , N1N2P1+P2 で は 上 方
に 回 答 し た 被 験 者 OTK の N1, N2, P1, P2 の 周 波 数 と レ
ベ ル の 関 係 を 示 し た も の で あ る .白 丸 は 0°,黒 丸 は 90°
の N1, N2, P1, P2 を 示 す . 2 本 の 破 線 は , Moore et al.
(1989)が 中 心 周 波 数 8 kHz,帯 域 幅 が 中 心 周 波 数 の 25%
のノッチを用いて求めた 3 人の被験者のノッチ検知閾
の最大値と最小値を示している.すなわち,ノッチレ
で は FEC 条 件 よ り も 準 FEC 条 件 の 方 が 𝐸𝐸𝐸𝐸は や や 小 さ
い . こ れ は , 準 FEC 条 件 で は FEC 条 件 と 比 較 し て 後
方 に 知 覚 す る 傾 向 が あ る こ と に 起 因 す る (FIGs. 5,8 参
ベ ル が –9dB 以 上 で は す べ て の 被 験 者 が ノ ッ チ を 検 知
照 ).後 方 に 音 像 を 知 覚 す る 本 質 的 な 理 由 は 不 明 で あ る
で き ず , –20dB 以 下 で は す べ て の 被 験 者 が 検 知 し た .
が , Fig.4 に 示 し た 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の peak-to-peak レ ン ジ の 広
0°の N1, N2 の レ ベ ル は と も に 検 知 閾 を 超 え て い る .
さ が 根 底 に あ る 理 由 の 1 つ か も し れ な い .準 FEC 条 件
し か し , 90°の N1 は 検 知 で き な い レ ベ ル で あ り , N2
において上方から後方への回答が頻繁にみられた被験
は被験者により検知の可否が分かれるレベルである.
者 OTK の 𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の peak-to-peak レ ン ジ は 他 の 被 験 者 と
周 波 数 軸 上 で み る と ,90°で は ノ ッ チ と ピ ー ク は 周 波
数 の 低 い 方 か ら P1, P2, N1, N2 の 順 で 並 ん で い る . P2
比べて大きく,定位誤差の大きさと対応している.
N1N2P1(b)で も 0°, 180°で は FEC 条 件 と 準 FEC の 𝐸𝐸𝐸𝐸
を 再 現 し な い 場 合 , P1 は N1 か ら 離 れ た 周 波 数 に あ る
の 差 は 小 さ い (–1.5°, 2.1°).し か し ,30°, 60°, 90°で は 準
た め 対 比 効 果 は 期 待 で き な い . し か し P2 を 再 現 す る
FEC 条 件 の 𝐸𝐸𝐸𝐸は FEC 条 件 と 比 較 し て 顕 著 に 大 き く ,そ
と , P2 か ら み た N1 の 相 対 的 な レ ベ ル は –14.7dB と な
る.これは被験者によってはノッチの検知が可能なレ
ベルである.
以 上 よ り ,P2 そ れ 自 体 は ス ペ ク ト ラ ル キ ュ ー で は な
い が ,N1 を 強 調 す る こ と に よ り ,正 中 面 上 方 の 音 像 定
位精度を向上させる重要な役割を果たしている.
条 件 の 𝐸𝐸𝐸𝐸は FEC 条 件 と 比 較 し て 顕 著 に 大 き く ,そ の 差
は そ れ ぞ れ 27.9°, 16.3°で あ っ た .
(a)MeasuredHRTF
Meanvertical
localizationerror,Ev [deg.]
FEC
50
0
Notch
undetectable
-10
30
P2
20
90°N1
partly
detectable
-20
detectable
-30
10
0
90°N2
4000
0
30
60
90
120
150
180
Ave.
Targetverticalangle[deg.]
0°N1
8000
12000
Frequency[Hz]
(b)N1N2P1
60
0°N2
0
Quasi-FEC
40
P1
16000
20000
FIG. 15. Relationship between the frequencies and the levels of
N1, N2, P1, and P2 for subject OTK. Open and closed circles
denote N1, N2, P1, and P2 for 0° and 90°, respectively. Broken
Meanvertical
localizationerror,Ev [deg.]
10
Level [dB]
N1N2P1+P2(c) で は , 30°, 90° を 除 き FEC 条 件 と 準
FEC の 𝐸𝐸𝐸𝐸の 差 は 小 さ い . し か し , 30°, 90°で は 準 FEC
60
20
-40
の 差 は そ れ ぞ れ 16.5°, 14.9°, 23.0°で あ っ た .
50
10
0
0
30
60
90
120
(c)N1N2P1+P2
60
Meanvertical
localizationerror,Ev [deg.]
実 測 HRTF(a)で は ,0°, 180°に お い て は ,FEC 条 件 と
準 FEC の 𝐸𝐸𝐸𝐸の 差 は 小 さ い (1.5°, 0.3°). し た が っ て , 前
後 方 向 の 音 像 定 位 に お い て は ,𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の 補 正 は 必 要 で
Ave.
Targetverticalangle[deg.]
50
比 較 す る . FIG. 16 に 各 HRTF と 各 条 件 の 𝐸𝐸𝐸𝐸を 示 す .
180
20
which is 8 kHz and the bandwidth of which is 25 % of the center
FEC 条 件 と 準 FEC 条 件 の 平 均 上 昇 角 定 位 誤 差 (𝐸𝐸𝐸𝐸)を
150
30
subjects for the detection of the notch, the center frequency of
5.3. FEC 条 件 と 準 FEC 条 件 の 比 較
Quasi-FEC
40
lines indicate the maximum and minimum thresholds of three
frequency, as reported by Moore et al. (1989).
FEC
FEC
Quasi-FEC
150
180
40
30
20
10
0
0
30
60
90
120
Ave.
Targetverticalangle[deg.]
は な い . し か し , 30°, 60°, 90°に お い て は , 準 FEC 条
FIG. 16. Mean vertical localization error, Ev, for the FEC and
件 の 𝐸𝐸𝐸𝐸は FEC 条 件 と 比 較 し て 顕 著 に 大 き く ,そ の 差 は
quasi-FEC conditions for (a) the measured HRTF, (b) N1N2P1,
そ れ ぞ れ 16.4°, 13.0°, 9.8°で あ っ た . 一 方 , 120°, 150°
and (c) N1N2P1+P2.
以 上 を ま と め る と ,実 測 HRTF,N1N2P1,N1N2P1+P2
の い ず れ に お い て も , 0°お よ び 120°か ら 180°で は 準
FEC 条 件 は FEC 条 件 と 同 等 の 音 像 定 位 性 能 を 有 し て い
る .し か し 30°か ら 90°に お い て は ,準 FEC 条 件 の 𝐸𝐸𝐸𝐸は
FEC 条 件 と 比 較 し て 顕 著 に 大 き く ,𝑀𝑀×𝑃𝑃𝑃𝑃𝑃𝑃の 補 正 が 必
要である.
Blauert, J. (1969/70). “Sound localization in the median plane,”
Acustica, 22, 205–213.
Butler, A., and Belendiuk, K. (1977). “Spectral cues utilizes in
the localization of sound in the median sagittal plane,” J. Acoust.
Soc. Am. 61, 1264–1269.
Carlile, S., Leong, P., and Hyams, S. (1997). “The nature and
distribution of errors in sound localization by human listeners,”
6. 結 論
本 論 文 で は , 音 像 定 位 実 験 1 に よ り , N1, N2, P1 に
P2 を 加 え て 再 構 成 し た パ ラ メ ト リ ッ ク HRTF
(N1N2P1+P2)に よ り 正 中 面 上 方 で 実 測 HRTF と 同 等 の
音像定位精度が得られるか否かを検証した.さらに,
音像定位実験 2 により,正中面上方の音像定位におけ
Hear. Res. 114, 179–196.
Gardner, B., and Gardner, S. (1973). “Problem of localization in
the median plane: effect of pinna cavity occlusion,” J. Acoust.
Soc. Am. 53, 400–408.
Hebrank, J., and Wright, D. (1974). “Spectral cues used in the
る P2 の 役 割 を 考 察 し た .
FEC 条 件 に お け る 結 果 を 以 下 に 示 す .
localization of sound sources on the median plane,” J. Acoust.
(1) N1N2P1 に P2 を 加 え る こ と に よ り ,正 中 面 の 90°
Soc. Am. 56, 1829–1834.
お よ び 120°で 音 像 定 位 性 能 が 向 上 し た .
(2) N1N2P1+P2 と 実 測 HRTF の 間 の 平 均 上 昇 角 定 位
誤 差 の 差 は 上 半 球 正 中 面 の 全 7 方 向 に お い て 10°未 満
であった.
Iida, K., Yairi, M., and Morimoto, M. (1998). “Role of pinna
cavities in median plane localization,” Proceedings of 16th
International Congress on Acoustics, 845–846.
(3) P1, P2, P1P2 を 再 現 し て も 音 像 は 上 方 に 知 覚 さ れ
Iida, K., Itoh, M., Itagaki, A., and Morimoto, M. (2007).
なかった.言い換えれば,これらのスペクトラルピー
“Median plane localization using parametric model of the
クは上方のスペクトラルキューではない.
head-related transfer function based on spectral cues,” Appl.
(4) 0°の N1 と N2 の レ ベ ル は と も に ノ ッ チ の 検 知 閾
を 超 え る が , 90°の N1 は 検 知 で き な い . し か し , P2
を 加 え る こ と に よ り , P2 か ら み た N1 の 相 対 レ ベ ル は
検知可能となる.
(5) 以 上 よ り , P2 は N1 を 強 調 す る こ と に よ り , 正
中面上方の音像定位精度を向上する重要な役割を果た
していると考えられる.
さ ら に , 準 FEC 条 件 に お け る 結 果 を 以 下 に 示 す .
(6) N1N2P1 に P2 を 加 え る こ と に よ り ,正 中 面 の 60°
か ら 120°で 音 像 定 位 性 能 が 向 上 し た .
(7) 0°お よ び 120°か ら 180°で は , FEC 条 件 と 同 等 の
音 像 定 位 性 能 が 得 ら れ た . し か し , 30°か ら 90°で は 音
像 定 位 性 能 は FEC 条 件 と 比 較 し て 低 下 し た .
謝
辞
Acoust. 68, 835–850.
Iida,
K.,
and
Ishii,
Y.
(2011).
“Individualization
of
the
head-related transfer functions in the basis of the spectral cues
for sound localization,” in Principles and Applications of Spatial
Hearing, edited by Y. Suzuki, D. Brungard, Y. Iwaya, K. Iida, D.
Cabrera, and H. Kato (World Scientific, Singapore), pp. 159–
178.
Iida, K., Nishioka, S., and Ishii, Y. (2014). “Personalization of
head-related transfer functions in the median plane based on the
anthropometry of the listener ’s pinnae,” J. Acoust. Soc. Am. 136,
317–333.
Kulkarni, A., and Colburn, H. S. (1998). “Role of spectral detail
in sound- source localization,” Nature 396, 747–749.
本 研 究 の 一 部 は 科 研 費 ( 基 盤 研 究 (A) 15H01790) に
Lopez-Poveda, E. A., and Meddis, R. (1996). “A physical model
より実施した.また,神戸大学の森本政之名誉教授に
of sound diffraction and reflections in the human concha,” J.
は有益な議論をしていただいた.記して感謝いたしま
Acoust. Soc. Am. 100, 3248–3259.
す.
Majdak, P., Goupell, M. J., and Laback, B. (2010). “3-D
文
献
Asano, F., Suzuki, Y., and Sone, T. (1990). “Role of spectral
cues in median plane localization,” J Acoust. Soc. Am. 88, 159–
168.
localization
of
virtual
sound
sources:
Effects
of
visual
environment, pointing method, and training,” Atten. Percept.
Psychophys. 72, 454–469.
Møller, H., Hammershøi, D., Jensen, C. J., and Sørensen, M. F.
(1995). “Transfer characteristics of headphones measured on
human ears,” J. Audio Eng. Soc. 43, 203–217.
Moore, B. C. J., Oldfield. R., and Dooley, G. J. (1989).
“Detection and discrimination of peaks and notches at 1 and 8
kHz,” J. Acoust. Soc. Am. 85, 820–836.
Morimoto, M., and Aokata, H. (1984). “Localization cues of
sound sources in the upper hemisphere,” J. Acoust. Soc. Jpn. (E)
5, 165–173.
Musicant,
A.,
and
Butler,
R.
(1984).
“The
influence
of
pinnae-based spectral cues on sound localization,” J. Acoust.
Soc. Am. 75, 1195–1200.
Raykar, V. C., Duraiswami, R., and Yegnanarayana, B. (2005).
“Extracting the frequencies of the pinna spectral notches in
measured head related impulse responses,” J. Acoust. Soc. Am.
118, 364–374.
Shaw, E. A. G., and Teranishi, R. (1968). “Sound pressure
generated in an external-ear replica and real human ears by a
nearby point source,” J. Acoust. Soc. Am. 44, 240–249.
Takemoto, H., Mokhtari, P., Kato, H., Nishimura, R., and Iida, K.
(2012). “Mechanism for generating peaks and notches of
head-related transfer functions in the median plane,” J. Acoust.
Soc. Am. 132, 3832–3841.