講義資料4-1 - 奈良先端科学技術大学院大学
Download
Report
Transcript 講義資料4-1 - 奈良先端科学技術大学院大学
音響信号処理特論
音響信号処理の基礎２
（その１）
2013年6月21日
奈良先端科学技術大学院大学
猿渡 洋
本日の予定
イントロダクション
音環境の制御
空間音響（聴空間の知覚）
人間はどのようにして音源位置を把握するのか
聴覚特性を利用した仮想音源の合成
音場再現（聴空間の合成）
音場制御とは？
逆フィルタ制御による音場再現
空間内の音を操る
従来の応用分野（音の前方，左右の定位感制御）
オーディオシステム
PA（Public Address）システム
新しい応用分野（真のステレオフォニクス）
空間上の「音」を消す
• 制御点での音圧を抑制（騒音の低減）
空間上に「音」を創る
• 厳密な音の再現（すべての音響情報の保存と再現）
• いわゆる「音のVirtual Reality」
音場制御：
空間の音の状態を物理的に制御
市販のオーディオシステム
ステレオからマルチチャンネルサラウンドへ
従来のステレオ再生
5.1chシステム
音空間の再現
これまでの再生手法
音の到来方向（音圧，到来時間差）を制御
さらに厳密に再現するために
元の音空間の持つ反射・残響などの空間的な情報を
再現
再生する室内の反射・残響の影響を除去
必要な知見
元の空間での音の伝播・音の知覚
再生する空間における音の振る舞い
空間知覚
音源（音像）の位置を知る能力
方向の知覚（方向定位）
距離の知覚（距離定位）
音源
距離
方向
音源の方向知覚（１）
左右の耳に到達する音波の
わずかな時間の違いと強さ
の違いを手がかり
到達距離に差
両耳間時間差
両耳間音圧差
Amp.
音圧差
time
時間差
Amp.
time
音源の方向知覚（２）
混同の円錐
2点間の音圧差と時間差だけでは，音源の位置を
一意的に明らかにすることはできない
音圧差と時間差の情報により，
音像が知覚される円錐底面の周
音源の方向知覚（３）
両耳間の時間差と音圧差以外の情報
胴体，頭部，耳介などにおける音波の回折によるスペ
クトルの変化（HRTF）
頭部の運動によるスペクトルの変化
視覚との相互作用によるスペクトルパターンの
学習
広帯域の音源は定位しやすい
高域成分が音像定位に寄与
純音では定位精度は悪くなる
音源の距離知覚
自由空間の場合
慣れ親しんだ音の場合は音圧レベルが主な手がかり
音源が多数ある場合には，音源の音圧レベルを相対的
に知覚
ある距離を越えると，空気吸音減衰が有効
有響下の空間の場合
直接音と反射音の強さの比
直接音と反射音間の時間遅れ
反射音は距離感だけでなく，拡がり感や方向知覚（定
位感）に影響を与える
先行音効果（１）
最初に到来する音源の方向が音像の定位に支配
的に影響すること
2つのスピーカの間に音像を知覚
右側のスピーカからの音しか
知覚できない
先行音効果（２）
先行音効果の成立する条件
コヒーレントな信号
信号の到来時間差や音圧差が影響
音圧
先行音効果の成立する領域
時間
頭部伝達関数（１）
Head Related Transfer Function: HRTF
自由空間のある１点から鼓膜面上までの伝達関数
• 室内の空間情報は一切含まれていない
• 純粋に胴体や，頭部，耳介の影響のみが含まれる
実際には外耳道入口での音圧が測定される
• マイクロホンを鼓膜面の直前まで入れるのは，被験者に
対する危険を伴う
• 外耳道伝達関数の個人差が大きい
音源の方向や距離によって変化する
頭部伝達関数（２）
音源が受聴者の右90度に存在するときの
左耳のHRTFの周波数振幅特性
HRTFを用いた仮想音の合成
ドライソースと任意の音源位置からの頭部
伝達関数を畳み込み
HRTF
convolution
空間内の音を操る
従来の応用分野（音の前方，左右の定位感制御）
オーディオシステム
PA（Public Address）システム
新しい応用分野（真のステレオフォニクス）
空間上の「音」を消す
• 制御点での音圧を抑制（騒音の低減）
空間上に「音」を創る
• 厳密な音の再現（すべての音響情報の保存と再現）
• いわゆる「音のVirtual Reality」
音場制御：
空間の音の状態を物理的に制御
音を消す
直感的に分かりやすく言うと
騒音から発した音波を他のスピーカから発した音波で
干渉させることにより，騒音を低減する方法
科学的に厳密に言うと
騒音の低減を目的として処理された信号により，ス
ピーカを駆動する方法
騒音源
スピーカ
制御点
Controller
能動騒音制御(Active Noise Control)
ANCの応用分野
自動車騒音（室内騒音，排気騒音，タイヤ騒音）
航空機の室内騒音
建築音響分野（塀，壁，吸音材，配管系の振動制御）
電気製品の騒音（冷蔵庫，掃除機）
発電所，工場における騒音
「音空間」の再現
再びマルチチャンネルサラウンド
従来のステレオ再生
音の移動感や包囲感が向上
5.1chシステム
場としての
臨場感は？？
音場再現
ダミーヘッド＋マイクロホン
室内音響伝達特性
逆システム
原音場
再現音場
音の「移動感」「距離感」「包囲感」の完全な再現
コンサートホール（原音場）における受聴者周辺の音空間（音場）を
リスニングルームにて再現
全ての音響情報の保存と再現
再現規模による分類
空間上の点を制御
バイノーラルシステム
トランスオーラルシステム
耳元周辺の制御
空間上のフィールドを制御
境界音場制御
多点入力・多点出力のトランス
オーラルシステム
（多チャネル音場再現システム）
頭部周辺の制御
ハードウェア的な分類
没入型
デバイスを身に付けることにより，仮想世界の合成
情報以外の感覚情報を遮断する
ヘッドホン，HMD
開放型
特定のデバイスを身に付けずに，日常生活において
違和感のない仮想世界を実現する
スピーカ
バイノーラルシステム
Inverse
System
Head And Torso Simulator (HATS)
原音場
聴取者
ダミーヘッド（もしくはHATS）を用いて収音，ヘッド
フォンで再生
システムが簡易
聴取者の頭部回転や移動に弱く，音像が頭内定位する
トランスオーラルシステム
Inverse
System
原音場
再現音場
ダミーヘッドで収音，2個以上のスピーカで再生
再現音場の伝達特性を考慮するため，逆フィルタ
を使用（残響やクロストークの抑圧）
境界音場制御に基づく音場再現
Control
area
Inverse
System
Kirchhoff-Helmholtz積分方程式と逆システム理論をベー
ス
ある境界上の音圧，法線方向粒子速度を完全に制御するこ
とにより，その境界内部の音場を完全に制御
音場再現システムの比較
名称
基本理論
制御領域
空間的な再現能力
モノフォニック
波動
なし
なし
ステレオ
波動
なし
拡がり感
マルチスピーカ
波動
スピーカに囲まれ
た内部
方向，拡がり感
バイノーラル
システム
両耳の2点
方向，距離，拡がり感
トランスオーラル
システム
両耳の2点
方向，距離，拡がり感
境界音場制御
波動＋システム
制御領域を自由に
決定可能
方向，距離，拡がり感
トランスオーラルシステムの作り方
H11
H21
G11
G12
G21
G31
H31
H12
H22
H32
G22
G32
原音を収録
再現する系の室内伝達関数を測定
伝達関数のインパルス応答から逆フィルタを設計
収録した原音に逆フィルタを畳み込む
室内における音の伝播
Loudspeaker
二次反射音
直接音
一次反射音
Microphone
壁面・物体での多重反射により複雑な音場を形成
室内伝達関数
音源からパルスを発生させたときの受音信号
単位インパルスに対する応答＝インパルス応答
室内の反射状況をインパルス信号列で表現
時間領域のインパルス応答
周波数振幅特性
室内伝達関数の測定方法
時間的に幅の狭いパルスをスピーカから出力し，
その応答を測定
全帯域で大きなエネルギーをもつ信号を用いて，
クロススペクトル法で計算
TSP（time stretched pulse）
その他（楽音信号，音声，白色雑音，その他）
Time Stretched Pulse (TSP)を用いた測定
インパルス応答の位相を周波数の2乗に比例して変化さ
せることにより，時間軸を引き伸ばした信号
時間的に瞬時周波数が均等に遷移
信号長が約0.005秒の場合
信号長が約11秒の場合
逆フィルタ
室内伝達特性の逆特性を持つフィルタ
入力信号に逆フィルタを畳み込むことにより，音源－
制御点間の反射や残響の影響を除去できる
逆フィルタ
1
H ()
G()
X ()
Loudspeaker
H ()
G()
D()
Microphone
多点逆システム
（MINT: Multi-input/output INverse Theorem）
X ()
H1 ()
H2 ()
G1()
X ()
D()
G2 ()
安定な逆システムを構成するために複数のスピーカ，逆
フィルタを用いる
一方の伝達系に０成分があっても，もう一方でカバーすること
が可能
多チャネル-多点制御系と逆フィルタ
X1
X2
S1
GJ1
H1J
HIJ
X1
C1
G1I
C2
G2I
GJI
CJ
X2
…
…
…
XJ
D1
D2
…
HI2
SI
G21
Control
point
…
HI１
…
XJ
H12
G11
…
X1
X2
H11
…
XJ
Inverse filter
Loudspeaker
DJ
I: 二次音源数， J: 制御点数
逆フィルタの基本的な設計法
時間領域での設計
二次音源数≦制御点数
二次音源数＞制御点数
最小二乗解（近似解）
MINT法（厳密解）
• 正確な逆フィルタを求めることが出来るが，計算量大
周波数領域での設計
二次音源数＜制御点数
二次音源数＝制御点数
二次音源数＞制御点数
最小二乗解（近似解）
唯一の解
SVDによる一般逆行列
• 少ない計算量で求めることが出来るが，必ずしも正確でない
特異値分解(SVD)による逆フィルタ設計
H11
H21
G11
G12
G21
G31
H31
H12
H22
G22
G32
H32
H11() H12()
G11() G21() G31()
H () H 21() H 22() G()
G
(
)
G
(
)
G
(
)
22
32
12
H31() H32()
特異値分解(SVD)による逆フィルタ設計
伝達関数 G n
を周波数ごとに特異値分解
I: 2 次音源数、 J: 制御点数 （I > J）
μj(ωn)： 特異値
I
行I
列
J
行
I
列
J
行
J
列
J
行
I
列
1 n
O
H
U
V
n
O
J n
伝達関数 直交行列 特異値が並ぶ行列 直交行列
G
特異値分解(SVD)による逆フィルタ設計
I
行J
列
1
I
行I
列
行J
列
J
行J
列
I
1
n
H
1
n
J n
逆フィルタ
H
O
V
U
O
特異値の逆数が並ぶ行列
特異値の逆数を 0 に置き換えることで逆特性の緩和が可能
0 に置き換える数を打ち切り数と呼ぶ
打ち切り数を増やすことでより緩和される
打ち切り数 0 の場合、最小ノルム解となる
多チャネル音場再現システム
Equalized Area
Inverse
System
Primary field
Secondary field
音場再現の応用例：
多ユーザ化
各制御点を独立に制御できるので，制御点の配置を工夫
することによって複数の受聴者に異なる音の提示が可能
Screen
こんにちは！
Visual data
Multimedia
contents
Hello!
Audio data
Signal
processing
高品質・高臨場感な音場再現をめざして
スピーカ数
受音点数
逆フィルタ
音場再現システム
制御領域・帯域
再生環境
室内伝達特性の変動の要因
物理的（物体的）な配置の変化
スピーカや受音点の位置
新たな物体の出現
室内の自然環境変化
温度（音速が変化）
湿度（周波数特性が変化）
風（見かけの音速変化や屈折による伝達経路の変化，
音源の指向特性の変化）
逆システムの安定性
伝達系の互いの経路が似て
いることを意味する
誤差の拡大（雑音感が強く
なる）
条件数
逆システムの安定性に関す
る指標
大きいほど逆行列が不安定
Condition number
逆システムが不安定だと
Frequency [Hz]
2次音源数16，制御点数6における条件数
スピーカ数と安定性
878.9 Hz
Condition number
Condition number
219.7 Hz
Number of secondary sources
Number of secondary sources
スピーカ数を増やすことにより，安定性が向上
制御点数と安定性
スピーカ数を16としたときの条件数の変化
制御点数を増やすことにより，伝達経路の独立性が低下
制御領域を広くしたり，制御周波数帯域を広くするため
には，多くの制御点が必要
室内伝達特性の時変性とその対処
逆フィルタのオフライン設計
特別なハードウェアを必要としない
伝達特性の変化には追従できない
逆フィルタのオンライン処理（適応フィルタ）
伝達特性の変化に追従可能
制御点には信号を逐次観測するためのセンサとなる
マイクロホンが必要
適応的な音場再現システム
ダミーヘッド
＋マイクロホン
③ 推定した変動に基づき，
逆フィルタを更新
スピーカ
ユーザー周りの
音空間を制御
逆フィルタ
エラーセンサ
② 観測した音と原
音から室内環境
の変動を推定
環境変動の推定
原音場
（コンサートホールなど）
① エラーセンサで
現在再生されて
いる音を収録
再現音場
（リスニングルーム）
アプローチ
伝達系の変動を推定し，それに応じた厳密な
逆フィルタを再設計
伝達系の変動に応じて逆特性を緩和した逆
フィルタを設計
伝達系の変動に応じてシステムの規模を変更
室内インパルス応答の温度変化の補正
線形伸縮処理によるインパルス応答補正
インパルス応答の時間軸伸縮処理
室温以外の環境条件が不変であれば，インパルス応答
を適当な伸縮率で伸縮すれば室温変化が吸収可能（引
地，板倉， 1994）
室内温度変化に対して，逆フィルタ設計
に用いる室内インパルス応答を伸縮し，
逆フィルタを再設計
再現精度の比較(温度変化：約１℃)
提案法は温度計に基づく伸縮処理よりも良好な結果
提案法では信号観測位置による差はない
逆フィルタの緩和処理
on-line adaptive-TSVD: OATSVD
適応的な逆フィルタの緩和に基づく環境変化に
ロバストな音場再現システム
逆フィルタの緩和
逆フィルタの精度を落とすことによって、音質劣化を防止
打ち切り特異値分解（TSVD）により逆フィルタを緩和
利点
全ての制御点に信号観測用センサを配置する必要がない
アルゴリズムも簡素
試聴実験による再現音の比較
再生音の音質
音像定位
音質は従来法に比べ向上
定位精度は従来法と同等
音場制御の未来図
交通情報
DVD鑑賞
（吹き替え音声）
ロードノイズ
エンジンノイズ
音声認識カーナビ
の使用
DVD鑑賞
（低音強調）