This article is a technical report without peer review, and its polished and/or extended version may be published elsewhere. Copyright ©20●● by IEICE
実環境雑音下におけるインパルス応答測定波形の
最適切り出し方法の検討
渋澤 功
†金田 豊
††東京電機大学大学院工学研究科 〒120-8551 東京都足立区千住旭町 5 番
E-mail: †[11kmc19@ms, kaneda@c].dendai.ac.jp
あらまし インパルス応答はスピーカや室内の音響特性についての情報を数多く含んでいる。そのためインパル ス応答の測定は大変重要である。しかし、測定環境において室内騒音などが付加されて SN 比が低下するという問 題がある。従来は通常,インパルス応答を時間軸上で切り出していたが、SN 比の良い周波数成分を切り捨ててしま うという問題があった。本報告では Wiener Filter と帯域別切り出しを利用して従来法の問題点を改善した手法を提 案した。そして、提案法を用いた雑音下における最適な切り出し法について検討した。その結果、1)従来法に比べ 提案法の方が真値との平均二乗誤差が低減し提案法の有効性を確認できた、2)実環境雑音下において非定常雑音が 発生しうる環境下では TSP 法よりも M 系列法のほうが雑音の影響を受けにくい測定が行えるということがわかっ た。 キーワード インパルス応答,切り出し,Wiener Filter,TSP 法、M 系列法A study on impulse response waveform extraction under environmental noise
Isao SHIBUSAWA
†Yutaka KANEDA
††Graduate school of Engineering, Tokyo Denki University,
5 Senju-Asahi-cho, Adachi-ku, Tokyo 120-8551, Japan
E-mail: †[11kmc19@ms, kaneda@c].dendai.ac.jp
Abstract The impulse response contains a large amount of information about the acoustic characteristics of a loudspeaker
or a room. Therefore, the measurement of the impulse response is very important. However, the SN ratio is reduced when enviromental noise is added. To increase the SN ratio, the impulse response can be extracted on the time axis. However, this discauses loss of the frequency components with a high SN ratio. In this paper, we propose a new method that mitigates the problems of the conventional method using a frequency-band-wise extraction and Wiener filter, and we study the extraction under real environmental noise. As a result, 1) the effectiveness of the proposed method is confirmed from the reduced mean squared error, 2) under a non-stationary noise environment, the MLS method is less affected by noise than the TSP method.
Keywords impulse response,response extraction,Wiener filter,TSP method, MLS method
1. は じ め に
イ ン パ ル ス 応 答 (IR:Impulse Response) は ス ピ ー カ や 室 内 の 音 響 特 性 に つ い て の 情 報 を 数 多 く 含 ん で い る 。 そ の た め 、IR の 測 定 は 大 変 重 要 で あ る [1]。し か し 、測 定 環 境 に お い て 室 内 騒 音 な ど が 付 加 さ れ て SN 比 が 低 下 す る と い う 問 題 が あ る 。 こ の 問 題 に 対 し て こ れ ま で 様 々 な 方 向 か ら の 改 善 が 試 み ら れ て き た [2]。本 稿 で は 実 環 境 雑 音 下 に お い て 計 測 し た IR 波 形 を 時 間 軸 上 で 切 り 出 す こ と で SN 比 を 向 上 す る 手 法 の 検 討 を 行 っ た 。 一 般 に IR の 測 定 は 十 分 に 長 い デ ー タ 長 で 行 わ れ る の で 、 測 定 結 果 の 後 半 は 雑 音 区 間 (も し く は IR が 雑 音 よ り 十 分 に 小 さ い 区 間 )と な っ て い る 。こ の 測 定 結 果 を そ の ま ま 使 う と 不 要 な 雑 音 が 多 く 含 ま れ る こ と で SN 比 が 悪 く な る 。 そ の た め 、 雑 音 部 分 を 切 り 捨 て て 必 要 な 部 分 の み を 切 り 出 す 操 作 が 必 要 と な る 。 図 1 は 定 常 雑 音 が 付 加 さ れ た IR 波 形 の モ デ ル 図 で あ る 。縦 ハ ッ チ 部 分 が IR、斜 ハ ッ チ 部 分 が 雑 音 を 表 し て い る 。 従 来 の 多 く の 方 法 で は イ ン パ ル ス 応 答 の 短 時 間 パ ワ ー が 雑 音 パ ワ ー と 等 し く な る 時 刻 (図 1 の 破 線 ) 付 近 で 切 り 出 し が 行 わ れ て き た [3][4]。こ の 時 間 波 形 に お け る 切 り 出 し (以 降「 波 形 切 り 出 し 」 と 呼 ぶ )は IR に 含 ま れ る 全 周 波 数 成 分 を 同 一 時 刻 で 切 り 出 す も の で あ る 。 し か し 、 実 際 に は 周 波 数 に よ っ て IR の 長 さ は 異 な り 、 ま た 雑 音 の パ ワ ー も 異 な る 。 図 2 は 、 実 測 し た 室 内 イ ン パ ル ス 応 答 (室 容 積120𝑚3、 残 響 時 間 0.6s)に 、 標 準 的 な 室 内 騒 音 で あ る Hoth 雑 音 [5]を 付 加 し た 場 合 の ス ペ ク ト ロ グ ラ ム で 、 横 軸 に 時 間 、 縦 軸 に 周 波 数 と し て 表 し た も の で あ る 。 雑 音 の 低 域 成 分 の パ ワ ー が 大 き い た め 、 図 1 の モ デ ル 図 に 示 す よ う に 約 30ms で IR の 短 時 間 と 雑 音 の 短 時 間 パ ワ ー は 等 し く な り 、 0~30ms の 範 囲 で IR は 切 り 出 さ れ る 。 こ の 波 形 切 り 出 し は 、 図 2 の 破 線 に 示 す よ う に 全 帯 域 同 一 時 刻 で 切 り 出 し が 行 わ れ る こ と に な る 。し か し 、 高 周 波 域 で は IR も 長 時 間 継 続 し て お り 、ま た 雑 音 成 分 も 少 な い 。 全 帯 域 で の 切 り 出 し を 行 っ た 場 合 、 図 2 に 丸 印 で 囲 っ た 、SN の 高 い 周 波 数 成 分 の 応 答 が 切 り 捨 て ら れ る た め 、 誤 差 要 因 と な る 。 こ の 問 題 を 解 決 す る た め 、 本 報 告 で は Wiener Filter を 用 い て 帯 域 別 の 雑 音 抑 圧 を 行 い 、 そ の 後 、 帯 域 別 の 切 り 出 し 時 刻 で IR の 切 り 出 し を 行 う 方 法 を 提 案 す る 。
2. Wiener Filtering
信 号 の SN 比 を 向 上 さ せ る 代 表 的 な 手 法 と し て Wiener Filtering(以 下 WF)が あ る [6]。こ の 手 法 は 測 定 信 号 の 平 均 二 乗 誤 差 を 最 小 に す る フ ィ ル タ を 用 い て 雑 音 を 抑 圧 す る 手 法 で あ る 。 WF は 測 定 信 号 の 短 時 間 フ ー リ エ 変 換 に 作 用 す る 時 変 フ ィ ル タ で 次 式 の 時 間 -周 波 数 特 性G(𝑛, 𝑘)を 持 つ 。 G(𝑛, 𝑘) = 𝑃𝑠(𝑛, 𝑘) (𝑃⁄ 𝑠(𝑛, 𝑘)+ 𝑃𝑁(𝑘)) (1) た だ し 、n は 離 散 時 間 を 表 し 、k は 離 散 周 波 数 を 表 す 。 ま た 、 雑 音 は 定 常 雑 音 を 仮 定 し𝑃𝑁(𝑘) は 雑 音 の 短 時 間 パ ワ ー ス ペ ク ト ル の 平 均 値 、𝑃𝑠(𝑛, 𝑘) は 信 号 (真 値 )の 短 時 間 パ ワ ー ス ペ ク ト ル を 表 す 。式 (1)よ り 、あ る 時 間 ― 周 波 数 区 間 に お い て 、𝑃𝑠(𝑛, 𝑘) ≫ 𝑃𝑁(𝑘) で あ れ ば G(𝑛, 𝑘) ≈ 1と な っ て 信 号 は そ の ま ま 出 力 さ れ る 。 一 方 、 𝑃𝑠(𝑛, 𝑘) ≪ 𝑃𝑁(𝑘) す な わ ち 、ほ ぼ 雑 音 の み の 区 間 で あ れ ばG(𝑛, 𝑘) ≈ 0 と な っ て 雑 音 は 抑 圧 さ れ る 。 し か し 、現 実 的 に は 真 値 は 知 る こ と が で き な い の で 雑 音 を 含 ん だ 測 定 結 果 の パ ワ ー ス ペ ク ト ル𝑃𝑆𝑁(𝑛, 𝑘) を 用 い て G(𝑛, 𝑘) = (𝑃𝑆𝑁(𝑛, 𝑘) − 𝑃𝑁(𝑘)) 𝑃⁄ 𝑆𝑁(𝑛, 𝑘) (2) と 近 似 計 算 し た 特 性 を 用 い る 。 図 1 雑 音 が 付 加 さ れ た IR と 波 形 切 り 出 し の モ デ ル 図 図 2 Hoth 雑 音 が 加 わ っ た イ ン パ ル ス 応 答 (IR)の ス ペ ク ト ロ グ ラ ム と 時 間 波 形 上 で 切 り 出 し (全 帯 域 で の 切 り 出 し )の 欠 点3. 提 案 す る 手 法
本 稿 で 提 案 す る 手 法 の 処 理 手 順 を 図 3 に 従 っ て 説 明 す る 。 ① 雑 音 が 加 わ っ た IR に 短 時 間 フ ー リ エ 変 換 を 行 っ て サ ブ バ ン ド 信 号 と す る 。 ② こ の サ ブ バ ン ド 信 号 に WF を 行 う 。た だ し WF に お け る 雑 音 の パ ワ ー ス ペ ク ト ル𝑃𝑁(𝑘) は 測 定 さ れ た 信 号 に お い て 雑 音 の み が 含 ま れ る と 考 え ら れ る 時 間 区 間 を 用 い て 計 算 し て お く 。ま た 、式 (2)の 分 子 項 が 負 値 に な る 場 合 はG(𝑛, 𝑘) = 0 と し た 。 ③ WF 出 力 に 対 し て 帯 域 別 切 り 出 し を 行 っ て イ ン パ ル ス 応 答 収 束 時 刻 以 降 の 雑 音 成 分 を 抑 圧 す る 。(詳 細 は 4 節 ) ④ サ ブ バ ン ド 信 号 を 逆 DFT し 、 overlap-add 法 に よ っ て IR を 再 合 成 す る 。図 3 提 案 法 の 処 理 手 順
4. 帯 域 別 切 り 出 し の 検 討
WF は 平 均 的 な 雑 音 パ ワ ー を 用 い て 雑 音 を 抑 圧 す る 。 従 っ て IR が 収 束 し て 、 雑 音 の み と な っ た 時 間 -周 波 数 区 間 (𝑃𝑆𝑁(𝑛, 𝑘)に 雑 音 の み し か 含 ま れ て い な い 区 間 ) で あ っ て も 、 雑 音 が 平 均 値 よ り も 大 き い 場 合 に は 式 (2) の 分 子 は 0 と は な ら ず 雑 音 は 残 っ て し ま い 誤 差 要 因 と な る 。 そ こ で 、 信 号 に WF 法 を 適 用 し た あ と に 帯 域 別 切 り 出 し を 行 っ て 誤 差 の 抑 圧 を 行 う 。 こ こ で 、式 (2)の 分 子 項 が 負 値 に な る 場 合 に つ い て 考 え て み る と 𝑃𝑆𝑁(𝑛, 𝑘) − 𝑃𝑁(𝑘) = 𝑃𝑆(𝑛, 𝑘) + 𝑃̃𝑁(𝑛, 𝑘) − 𝑃(𝑘) ≤ 0 (3) す な わ ち 𝑃𝑆(𝑛, 𝑘) ≤ P(k) − 𝑃̃𝑁(𝑛, 𝑘) (4) た だ し 、𝑃𝑆(𝑛, 𝑘)は IR の 真 値 、 𝑃̃𝑁(𝑛, 𝑘)は 雑 音 パ ワ ー の 時 間 -周 波 数 瞬 時 値 で あ る 。 式 (4)は IR の 真 値 が 雑 音 の 平 均 値 か ら の 瞬 時 変 動 よ り 小 さ い 。言 い 換 え る と IR が 雑 音 の 変 動 以 下 に 収 束 し て い る こ と を 表 す 。よ っ て 本 手 法 で は 式 (2)の 分 子 項 が 最 初 に 負 値 と な っ た 時 刻 を IR の 切 り 出 し 時 刻 と 定 め る も の と す る 。 次 に こ の 事 の 妥 当 性 を 確 認 す る た め に 行 っ た シ ミ ュ レ ー シ ョ ン に つ い て 述 べ る 。 ま ず 、 高 い SN 比 で 室 内 イ ン パ ル ス 応 答 ( 残 響 時 間 0.6 秒 ) を 測 定 し 、 こ れ を 「 真 値 」 と 考 え る 。 次 に 、 こ れ に Hoth 雑 音 を 付 加 し た も の を 処 理 対 象 と し た 。イ ン パ ル ス 応 答 は 窓 長 32 の ハ ニ ン グ 窓 を 用 い シ フ ト 長 1 で 短 時 間 フ ー リ エ 変 換 を 行 っ た 。図 4,5,6 に そ れ ぞ れ 、 真 値 波 形 、 WF を 行 う 前 お よ び 後 の 波 形 を 示 す 。 た だ し 、波 形 は 750-2250Hz 帯 域 の サ ブ バ ン ド 信 号 を 時 間 波 形 に 変 換 し て 示 し た 。図 よ り 、WF 処 理 前 (図 5)に 見 ら 図 4 IR の 真 値 波 形 (750-2250Hz) 図 5 WF を 適 用 す る 前 の 750-2250Hz 帯 域 の IR 図 6 WF を 適 用 し た 後 の 750-2250Hz 帯 域 の IR れ た 雑 音 成 分 は 処 理 後 (図 6)に お い て 大 幅 に 抑 圧 さ れ て い る こ と が わ か る が 、消 し 残 り 雑 音 波 形 が 見 ら れ る 。 次 に 図 6 の 波 形 を さ ま ざ ま な 時 刻 で 切 り 出 し た 場 合0
200
400
600
800
-10
-5
0
5
時間 (ms)
相対振幅
0
200
400
600
800
-10
-5
0
5
時間 (ms)
相対振幅
0
200
400
600
800
-10
-5
0
5
時間 (ms)
相対振幅
消 し 残 り 雑 音の 真 値 と の 誤 差 を 計 算 し た 。図 7 は 、図 6 の IR の 切 り 出 し 時 刻 τ を 横 軸 に 、 縦 軸 を そ の 切 り 出 し 時 刻 に お け る IR 波 形 誤 差E(τ) E(τ) = ∑ (|ℎ(𝑛)| − |ℎ̂(𝜏, 𝑛)|)2 𝑛 ∑ |ℎ(𝑛)| 2 𝑛 ⁄ を 表 し て い る 。 た だ し 、ℎ(𝑛), ℎ̂(𝜏, 𝑛)は そ れ ぞ れ 真 値 お よ び WF 適 用 後 に 切 り 出 し た IR 波 形 (750-2250Hz)で あ る 。 図 よ り 最 も 誤 差 が 小 さ か っ た 時 刻 は 式 (2) の 分 子 項 が 最 初 に 負 値 と な る 時 刻 に 近 い 結 果 と な っ た 。 ま た 、 ほ か の 周 波 数 帯 域 に つ い て も 同 じ よ う な 結 果 と な り 、 提 案 す る 切 り 出 し の 時 刻 の 妥 当 性 が 確 認 さ れ た 。 最 後 に 時 間 波 形 で 切 り 出 し を 行 っ た 場 合 (波 形 切 り 出 し )と の 性 能 比 較 を 行 っ た 。シ ミ ュ レ ー シ ョ ン 条 件 は 同 じ も の と し た 。 波 形 切 り 出 し を し た 結 果 は 図 2 の ス ペ ク ト ロ グ ラ ム の 切 り 出 し 時 刻 の 左 側 の 部 分 で あ る 。 一 方 、図 8 に 提 案 法 を 適 用 し た IR の ス ペ ク ト ロ グ ラ ム を 示 す 。波 形 切 り 出 し 法 で は 切 り 捨 て ら れ て い た 30ms 以 降 の 部 分 が 提 案 法 で は 保 存 さ れ て い る 。 次 に 波 形 切 り 出 し お よ び 提 案 法 で 得 ら れ た イ ン パ ル ス 応 答 を フ ー リ エ 変 換 し て 周 波 数 振 幅 特 性 を 求 め 、 真 値 の 周 波 数 振 幅 特 性 と の 誤 差 を 周 波 数 帯 域 ご と に 計 算 し 、 そ の 結 果 を 図 9 の 太 線 と 破 線 に 示 す 。 図 よ り 提 案 法 の 方 が 3k~ 12kHz に お い て 10dB 程 度 誤 差 が 減 少 し て い る こ と が わ か る 。 こ の こ と よ り 提 案 法 の 有 効 性 を 確 認 す る こ と が で き た 。
5. 実 環 境 測 定 の 問 題 点 と 対 策
提 案 法 は 定 常 雑 音 を 前 提 と す る も の で あ る 。 し か し 、 実 環 境 に は ド ア の 開 閉 音 や 歩 行 音 、 衝 撃 音 な ど の 非 定 常 雑 音 が 発 生 す る 可 能 性 が あ る 。 以 下 で は こ の 非 定 常 雑 音 が 本 手 法 に 及 ぼ す 影 響 に つ い て 検 討 し た 。 IR 測 定 信 号 の 代 表 的 な も の と し て 掃 引 正 弦 波 の 一 つ で あ る TSP 信 号 と 擬 似 雑 音 の 一 つ で あ る M 系 列 法 が 知 ら れ て い る 。 こ れ ら 2 つ の 信 号 を 用 い た 場 合 の 非 定 常 雑 音 の 影 響 に つ い て シ ミ ュ レ ー シ ョ ン を 行 っ た 。 ま ず 最 初 に TSP 法 を 用 い た シ ミ ュ レ ー シ ョ ン 結 果 に つ い て 述 べ る 。IR は こ れ ま で と 同 じ 室 内 IR を 用 い た 。 こ の IR を TSP 信 号 と た た み こ ん で TSP 応 答 を 計 算 し 、 こ れ に 非 定 常 雑 音 を 加 算 し た 。 非 定 常 雑 音 と し て は 、 別 途 録 音 し た 机 を 叩 く 音 を 用 い 、 そ の ス ペ ク ト ロ グ ラ ム を 図 10 に 示 す 。非 定 常 雑 音 に は 室 内 の 低 周 波 定 常 雑 音 が 付 加 さ れ た も の と な っ て い る 。 そ し て 、 TSP 応 答 に 非 定 常 雑 音 を 付 加 し た 時 の ス ペ ク ト ロ グ ラ ム を 図 11 に 示 す 。図 に お い て 時 刻 0 か ら 直 線 的 に 周 波 数 が 上 昇 す る TSP 応 答 に 、300ms 付 近 で 非 定 常 雑 音 が 加 算 さ 図 7 切 り 出 し 時 刻τと IR 波 形 の 誤 差 E(τ) 図 8 提 案 法 を 適 用 し た IR の ス ペ ク ト ロ グ ラ ム 図 9 周 波 数 振 幅 特 性 の 周 波 数 帯 域 ご と の 誤 差 れ て い る 。 次 に TSP 応 答 に 逆 TSP 信 号 を 畳 み 込 み IR を 求 め た 。 求 め た IR の ス ペ ク ト ロ グ ラ ム を 図 12 に 示 す 。 図 よ り TSP 法 で 求 め た IR は 非 定 常 雑 音 の 影 響 を200
400
600
800
-12.1
-12
-11.9
-11.8
-11.7
切り出し時刻τ (ms)
誤差
(d
B
)
時間 (ms)
周波数
(k
H
z)
0
30
100
200
300
400
0
5
10
15
20
-80
-60
-40
-20
0
(dB)
0
5
10
15
-30
-25
-20
-15
-10
-5
0
周波数帯域 (kHz)
誤差
(d
B
)
波形切出
提案法
図 10 今 回 用 い た 非 定 常 雑 音 強 く 受 け て い る と い う こ と が わ か る 。 次 に 図 12 の IR に 対 し て 提 案 法 を 適 用 し た IR の ス ペ ク ト ロ グ ラ ム を 図 13 に 示 す 。図 よ り 非 定 常 雑 音 下 で は 提 案 法 を 用 い て も 雑 音 が 抑 圧 で き ず 大 き な 誤 差 が 発 生 し て い る こ と が 分 か る 。ま た 、WF で は IR の 存 在 し て い な い 時 間 区 間 (た と え ば 図 12 の 400ms 以 降 )に お い て 定 常 雑 音 の パ ワ ー ス ペ ク ト ル𝑃𝑁(𝑘)を 推 定 す る が 、こ の 区 間 に 非 定 常 雑 音 が 存 在 す る と𝑃𝑁(𝑘)が 誤 推 定 さ れ 、 大 き な 誤 差 要 因 と な る 。 次 に M 系 列 を 用 い た シ ミ ュ レ ー シ ョ ン 結 果 に つ い て 述 べ る 。シ ミ ュ レ ー ト し た M 系 列 応 答 に 非 定 常 雑 音 を 付 加 し た 場 合 の ス ペ ク ト ロ グ ラ ム を 図 14 に 示 す 。図 よ り M 系 列 応 答 の ス ペ ク ト ロ グ ラ ム は 測 定 時 間 全 体 に ラ ン ダ ム に 分 布 し て い る こ と 、 お よ び 、 300ms 付 近 に 非 定 常 雑 音 が 存 在 す る こ と が わ か る 。次 に M 系 列 応 答 に 逆 M 系 列 信 号 を 畳 み 込 ん で 求 め た IR の ス ペ ク ト ロ グ ラ ム を 図 15 に 示 す 。図 よ り M 系 列 法 で 求 め た IR に お い て は 、非 定 常 雑 音 が 測 定 時 間 全 体 的 に 拡 散 し て 、 定 常 雑 音 化 し て い る こ と が わ か る 。 次 に 図 15 の IR に 対 し て 提 案 法 を 適 用 し た IR の ス ペ ク ト ロ グ ラ ム を 図 16 に 示 す 。非 定 常 雑 音 が 定 常 化 さ れ て い る の で 提 案 法 が う ま く 機 能 し 雑 音 成 分 の 抑 圧 が で き て い る こ と が 分 か る 。 図 17 に 図 13 お よ び 図 16 の 場 合 の 周 波 数 帯 域 ご と の 誤 差 を 示 す 。 図 よ り 、 非 定 常 雑 音 の 影 響 を 受 け て い な い 高 周 波 域 で は TSP の 誤 差 が 小 さ い が 全 体 と し て は M 系 列 法 が 有 利 で あ る こ と が 分 か る 。 以 上 の 結 果 か ら 実 環 境 に お い て 非 定 常 雑 音 が 混 入 す る 可 能 性 の あ る 部 屋 な ど で 測 定 を す る 際 に は M 系 列 法 を 用 い て 測 定 す れ ば 提 案 法 を 有 効 に 使 用 で き る と い う こ と が わ か っ た 。 図 11 非 定 常 雑 音 を 付 加 し た TSP 応 答 図 12 TSP 法 で 求 め た IR 図 13 図 12 の IR に 提 案 法 を 適 用 し た IR
時間 (ms)
周波数
(k
H
z)
0
200
400
600
0
5
10
15
20
-80
-60
-40
-20
0
(dB)
時間 (ms)
周波数
(k
H
z)
0
200
400
600
0
5
10
15
20
-80
-60
-40
-20
0
(dB)
時間 (ms)
周波数
(k
H
z)
0
200
400
0
5
10
15
20
-80
-60
-40
-20
0
(dB)
時間 (ms)
周波数
(k
H
z)
0
200
400
0
5
10
15
20
-80
-60
-40
-20
0
(dB)
図 14 非 定 常 雑 音 を 付 加 し た M 系 列 応 答 図 15 M 系 列 法 で 求 め た IR 図 16 図 15 の IR に 提 案 法 を 適 用 し た IR 図 17 TSP 法 と M 系 列 法 で の 周 波 数 振 幅 特 性 の 周 波 数 帯 域 ご と の 誤 差
6. む す び
本 稿 で は 、 計 測 し た イ ン パ ル ス 応 答 (IR:Impulse Response) に 含 ま れ る 雑 音 成 分 の 低 減 方 法 と し て Wiener Filtering(WF)と 帯 域 別 切 り 出 し を 利 用 し た 手 法 を 提 案 し た 。 WF は 通 常 、 残 留 雑 音 の 問 題 が あ る こ と が 知 ら れ て い る 。 そ こ で 、 提 案 法 で は 、 帯 域 分 割 を 行 っ た イ ン パ ル ス 応 答 に WF を 適 用 し た 後 、 出 力 が 最 初 に 零 と な る 時 刻 (測 定 結 果 の 短 時 間 パ ワ ー が 平 均 的 雑 音 パ ワ ー よ り 小 さ く な る 時 刻 )以 前 の 帯 域 信 号 を 切 り 出 す こ と で 残 留 雑 音 の 低 減 を 行 う 。室 内 IR を 用 い た シ ミ ュ レ ー シ ョ ン に よ り 、 こ の 切 り 出 し 時 刻 の 妥 当 性 、 お よ び 時 間 波 形 そ の も の を あ る 時 刻 で 切 り 出 す 方 法 に 対 す る 有 効 性 を 示 し た 。 今 回 提 案 し た 方 法 は 定 常 雑 音 を 前 提 と し て い る が 実 環 境 で は 非 定 常 雑 音 の 影 響 を 受 け る 場 合 も 多 い 。 イ ン パ ル ス 応 答 測 定 信 号 で あ る TSP と M 系 列 を 比 較 し た 結 果 、 非 定 常 雑 音 環 境 下 で は 、 非 定 常 雑 音 を 定 常 雑 音 化 す る M 系 列 法 が 本 手 法 に は 適 し て い る こ と を 示 し た 。文 献
[1] 橘 秀 樹 , ‘’ 室 内 音 響 測 定 の 現 状 と 今 後 の 課 題 , ’’ 日 本 音 響 学 会 誌 49(2), pp. 97 -102, (1993). [2] 落 合 裕 一 ,金 田 豊 , ‘’ 全 帯 域 で S N 比 を 一 定 と す る イ ン パ ル ス 応 答 測 定 信 号 (C SN-TSP) の 実 環 境 に お け る 有 効 性 の 検 討 ,’’ 音 響 学 会 秋 季 講 演 論 文 集 , 3-10-11, pp.621-622, (2010/09).[3] L. Fa iget, C . Le gros , and R . R ui z, ’’ Opt imizat io n o f t he I mp uls e R es po ns e Le ngt h: App licat io n to No is y a nd Highly R eve rbera nt R oo ms , ’’ J. Aud io E ng. Soc , 46, 9, (1998 ). [4] ISO 3382 - 1 :2009 (E) pp. 8 . [5] IE EE P2 69/D9 , pp. 80 -82, J an. 200 2. [6] 浅 野 太 , ‘’ 音 の ア レ イ 信 号 処 理 ’’ pp. 48 -54, 日 本 音 響 学 会 (2011 ).