1.基本事項
[1.1] パポリス: 工学のための応用フーリエ積分, pp.330, オーム社, 1967.
[1.2] 大賀寿郎, 山崎芳男, 金田豊, 音響システムとディジタル処理, p.124, 電子情報通信学会, 1995.
インパルス応答の応用
[1.3] 橘秀樹, 矢野博夫, 環境騒音・建築音響の測定, p.155, コロナ社, 2004.
[1.4 ] H. Kuttruff, Room Acoustics. London, Elsevier Science Publishers, 1973, pp. 231-243.
2.解説論文
[2.1] 金田: “ インパルス応答測定の際の留意点 ”.音響学会誌 55巻5号 pp.364-369 (1999). [2.2] 佐藤史明, “はじめてのインパルス応答計測,” 音響学会誌, 67, 4, 155-162 (2011). [2.3 ]佐藤史明, “Swept-Sine 法に基づく音響伝搬測定,” 音響学会誌, 63, 6, 322-327 (2007). [2.4] 金田, "はじめての音響信号処理 -ディジタル録音と補間の話-," 音響学会誌, 65巻10号, pp. 531-536 (2009.10 ).3.測定信号
一般
[3.1] S. Muller and P. Massarani, ‘‘Transfer- function measurement with sweeps,’’ J. Audio Eng. Soc., vol. 49, no. 6, pp. 443-471 (2001 June).
「掃引正弦波 SS 信号を用いた測定を中心とした良質な長編解説」 [3.2] 「有色疑似雑音(PN信号)の測定信号作成」
篠原亮,金田豊:"有色疑似雑音を用いたインパルス応答測定の検討," 音響学会秋季講演論文集, 1-Q-26, pp. (2013/09).
TSP
[3.3] N. Aoshima, “Computer-generated pulse signal applied for sound measurement,” J. Acoust. Soc. Am., vol. 69, no. 5, pp. 1484-1488 (1981 May).
[3.4] 鈴木陽一、浅野太、金学胤、曽根敏夫:“時間引き延ばしパルスの設計法に関する考察,” 信学技報、 vol. EA92-86
(1992.12).
[3.5] Y. Suzuki, F. Asano, H. Kim and T. Sone, “An optimum computer-generated pulse signal suitable for the measurement of very long impulse response,” J. Acoust. Soc. Am., vol. 97, no. 2, pp. 1119-1123 (1995 Feb.). [3.6] http://asano.media-interaction.jp/English/doc/tsp/index.htm
Log-SS
[3.7]伊藤次男, “音響測定方法およびその装置,” 特開平5-118906, (1993). [3.8] 藤本卓也, “低域バンドでのSN比改善を目的としたTSP信号に関する検討,” 1999年秋季音響学会講 演論文集 pp. 433-434 (1999.9).[ 参考文献 ]
M系列
[3.12] J. Borish, “An efficient algorithm for measuring the impulse response using pseudorandom noise,” J. Audio Eng. Soc., vol. 31, no. 7/8, pp. 478-488 (1983 July/Aug.).
[3.13] J. Borish, “Self-contained crosscorrelation program for maximum-length sequences,” J. Audio Eng. Soc., vol. 33, no. 11, pp. 888-891 (1985 Nov. ).
[3.14] D. D. Rife and J. Vanderkooy, “Transfer-function measurement using maximum-length sequences,” J. Audio Eng. Soc., vol. 37, no. 6, pp. 419-444 (1989 June).
[3.15] J. Vanderkooy, “Aspects of MLS measuring systems,” J. Audio Eng. Soc., vol. 42, no. 4, pp. 219-231 (1994 Apr.). [3.16] 次の文献の参考文献を参照:金田豊, “M系列を用いたインパルス応答測定における誤差の実験 的検討", 音響学会誌, 52, 10, 752-759 (1996).
その他
[3.17] 森勢, ほか "暗騒音と高調波ひずみに頑健なインパルス応答測定用信号," 電子情報通信学会論文誌, J89-A(1), pp. 7-14 (2008 ). [3.18] G. B. Stan, J. J. Embrocatesa and D. Archamveau, “Comparison of different impulse response measurement techniques,” J. Audio Eng. Soc., vol. 50, no. 4, pp. 249-262 (2002 Apr.).適応形
雑音白色化
[3.19] S. Weinzierl, A. Giese and A. Lindau, “Generalized multiple sweep measurement,” in 126th AES Convention, 7767 (2009 May).
雑音最小化
[3.20] N. Moriya and Y. Kaneda, “Impulse response measurement that maximizes signal-to-noise ratio against ambient noise,” Acoust. Sci. & Tech., vol. 28, no. 1, pp. 43-45 (2007 Jan.).
[3.21] 守谷直也, 金田豊, "雑音に起因する誤差を最小化するインパルス応答測定信号," 音響学会誌, 64巻12号, pp. 695-701 (2008).
SN比一定
[3.22] 落合裕一,金田豊, "全帯域でSN比を一定とするインパルス応答測定法の検討," 音講論集, pp. 879-880 (2010.3).
[3.23] H. Ochiai and Y. Kaneda, “Impulse response measurement with constant signal-to-noise ratio over a wide frequency range,” Acoust. Sci. & Tech., vol. 32, no. 2, pp. 76-78 (2011 Mar.).
[3.24] 中原優樹,金田豊: "CSN-SS信号による残響時間測定効率化の検討," 音響学会秋季講演論文 集, 1-Q-43, pp. (2014/09).
4.誤差
[4.1] 金田: “M系列を用いたインパルス応答測定における誤差の実験的検討”. 音響学会誌 52巻10号 pp. 752-759 (1996).定常雑音
[4.2] 中重亮太,金田豊, “各種インパルス応答測定信号の雑音低減効果について," 音講論集, pp.821-822 (2013.3). [4.3] 「測定結果の2周期目切り出し時のクロスフェード接続」 中重亮太,金田豊:"周期的信号を用いたインパルス応答測定の問題点について," 音響学会秋季講演論 文集, 1-Q-28, pp. 763-764 (2013.9).非線形
[4.4] 守谷直也, 池田亜希, 金田豊, "インパルス応答計測におけるスピーカの非線形歪みに関する検討," 音講論集, pp. 735-736 (2004.9).[4.5] N. Moriya and Y. Kaneda, “Study of harmonic distortion on impulse response measurement with logarithmic time stretched pulse,” Acoust. Sci. & Tech., vol. 26, no. 5, pp. 462-464 (2005 Sept.). [4.6] A. Torras-Rosell and F. Jacobsen, “A new interpretation of distortion artifacts in sweep measurements,” J. Audio Eng. Soc., vol. 59, no. 5, pp. 283-289 (2011 May).
[4.7] C. Dunn and M. O. Hawksford, “Distortion immunity of MLS-derived impulse response measurements,” J. Audio Eng. Soc., vol. 41, no. 5, pp. 314-335 (1993 May).
[4.8] 「非線形歪による主応答の変化」 佐々木長閑,金田豊:"インパルス応答に及ぼすスピーカーの非線形歪の影響," 音響学会秋季講演論文集, 1-Q-27, pp. 761-762 (2013/09). [4.9] 「Log-SS でなくても高調波歪の分離除去はできる」 葛山亮介, 金田豊:"インパルス応答測定用スウィープ信号の高調波歪分離の検討," 音響学会春季講演 論文集, pp.645-646 (2007/3). [4.10] 「DA変換の逆折り返しに起因した混変調歪」 佐藤憲孝,佐々木長閑,金田豊:“AD/DA変換器の非線形特性がインパルス応答に及ぼす影響の検討,” 音響学会春季講演論文集, 1-P4-2, pp. (2014/03). [4.11] 「電源雑音に起因した混変調歪」 佐藤憲孝,金田豊:“掃引正弦波を用いたインパルス応答測定の際に発生する非線形誤差の検討,” 音響学会秋季講演論文集, 1-Q-42, pp. (2014/09).
時変
[4.12] 佐藤史明: “ 室内音響インパルス応答の測定技術 ”. 音響学会誌 58巻10号 pp.669-676 (2002).6.AIF(オーディオインタフェース) の課題
[6.1] 志賀, ほか, "オーディオ・インタフェースの入出力フィルタ特性の検討 ," 信学技報、vol. EA2009-67, (2009.10). [6.2]小林慶弘, 村澤良太, 金田豊, "オーディオインタフェースの特性評価における直流除去特性の影響に ついて," 音響学会秋季講演論文集, pp.769-770 (2008/9). [6.3] 「DAによるクリッピング」 志賀, ほか, "オーディオ・インタフェースのDA変換器におけるクリッピング歪について," 2010年秋季音響 学会講演論文集 2-P-2 (2010.9). [6.4] 「DA・AD 非同期性」 守谷、金田、阪内 : “PC における AD・DA 同時動作の問題点”, 2003年秋季音響学会講演論文集 pp.609-610 (2003.9). [6.5] 「Win XP のクリッピング例」 郡司龍和, 金田豊, "PCによるA/D変換時に発生するクリッピング歪について," 音響学会秋季講演論文集, 1-P-18, pp.747-748 (2011/09).7.測定信号が利用できない場合
[7.1] 山崎芳男, 金田豊, “音・音場のディジタル処理,” p.154- , コロナ社, 2002. [7.2] S. Haykin, 鈴木博訳 “適応フィルタ理論”, p.22-, 科学技術出版, 2001. [7.3] 城戸健一, “ディジタルフーリエ解析(Ⅱ)”, p.68-, コロナ社, 2007.付録目次
付録1.2-1 時不変線形系と正弦波信号 付録3.3-1 時間波形の振幅とスペクトル
付録5.1-1 測定信号による雑音抑圧効果の計算
1 0 0sin
)
(
n n nn
t
A
A
t
y
正弦波信号は、一般形として、A・sin(ωt+θ) と表される。Aは振幅、ω(=2πft)は角周波数、θは位相である。 時不変線形系に、信号 s(t)=sin(ω0t) (A=1、ω=ω0、θ=0 )を入力したときの出力をy(t)と表す。 時不変線形系 s(t) y(t) これを数式で、L
[ s(t) ] = y(t) (1) と表すことにする。以下では、出力 y(t)も同じ周波数の正弦波となることを示す。 入力正弦波の周期 T0だけs(t)を遅らせた信号 s(t-T0) を入力した場合を考える。 時不変性より、出力もT0だけ遅れた信号となる。すなわち、L
[ s(t-T0) ] = y(t-T0) (2) 正弦波は1周期遅らせると同一信号となる、すなわち、 s(t-T0)= s(t) (3) であるので、式(2)に式(3)を代入し、式(1)の関係より、 y(t-T0)=L
[ s(t-T0) ] =L
[ s(t) ] = y(t) (4) よって、s(t) の出力y(t)は、周期T0の周期関数であることがわかる。 (5) y(t) も s(t) と同じ周期 T0(=2π/ω0) を持つ周期関数であることの証明 Step1 Step2 式(5)右辺の、周波数がmω0の正弦波成分に注目して、その半周期の時間をTm=π/(mω0) と表す。 そして、s(t)+s(t-Tm) を入力した場合の出力を考える。線形和の性質および時不変性の性質を利用して、L
[ s(t)+s(t-Tm) ] =L
[ s(t) ] +L
[ s(t-Tm) ] = y(t)+y(t-Tm) (6) 式(6)右辺に、式(5) の関係を代入すると、 フーリエ級数の知識より、周期T0を持つ信号 y(t) は、基本角周波数ω0(=2π/T0) および、 その倍周波数の正弦波の和として、次式のように表される。 y(t) が、周波数ω0の正弦波となることの証明 付録1.2-1 時不変線形系と正弦波信号 時不変線形系に、ある周波数の正弦波を入力したら、出力は同じ周波数の正弦波であることの証明 付録3.3-1 時間波形の振幅とスペクトル
m n n m n m m m m n n m n m m m m n n n m m mT
t
n
B
T
t
m
B
A
T
t
n
A
T
t
m
A
A
t
n
A
t
m
A
A
T
t
y
t
y
)
2
/
(
sin
)
2
/
(
sin
2
)
(
sin
)
(
sin
sin
sin
)
(
)
(
0 0 0 0 0 0 0 0 0 式(7)において、Bn(n=1,2,3,・・・) は、同じ周波数の2つの正弦波、An・sin(nω0t+θn) と An・sin(nω0(t-Tn)+θn) とを加算した場合の振幅を表しており、L
[ s(t)+s(t-Tm) ] =L
[ B・ sin(ω0(t-Tm/2 )) ] = B・y(t-Tm/2 )
m n m n m mm
t
T
BA
n
t
T
BA
BA
0sin
0(
/
2
)
sin
0(
/
2
)
(7) (12) 式(7) と 式(12) における、周波数が mω0 の正弦波の振幅は等しいとおいて、式(9)より、0
m mB
BA
(13) の関係を得る。式(11)より、m=1 ( T1=π/(ω0) ) の場合を除いて B ≠ 0 であるので、0
mA
を得る。以上の議論は、m=1 以外の正弦波に適用できるので、Am=0 (m>1) となり、式(5) より、 (14)
0 1
1 0sin
)
(
t
A
A
t
y
となる。最後に、T1=π/(ω0) の場合を考えると、s(t)+s(t-T1)=0 であるので、線形系の (15) Bn=2・An・cos(nω0・Tn/2) (8) である(三角関数の和の公式より得られる)。 n=m の場合、 Tm=π/(mω0) の関係を代入すると、Bm =2・Am・cos(mω0・Tm/2) = 2・Am・cos(mω0・(π/(mω0) ) /2) = 2・Am・cos(π /2) = 0 (9)
となる。 (周波数mω0の正弦波は、半周期ずらして加算されているので0 となる)
一方、入力信号は、三角関数の和の公式より
s(t)+s(t-Tm) = sin(ω0t)+ sin(ω0(t-Tm )) = B・sin(ω0t-Tm/2) (10)
ただし、 B=2・cos(ω0・Tm/2) (11)
図のような有限長Nの離散SS(掃引正弦波)信号 x(n) を考える。この信号は区間J以外では0または十分に 小さな値を持つものとする。Jを信号x(n)の実効長と呼び、信号の存在する時間区間を実効区間と呼ぶことにす る。さらに、この信号をN点DFTしたもの(周波数スペクトル)をX(k)と表す。 ただし、本節では離散信号x(n)は n=1~Nの時刻で表され、これをDFTしたX(k)はk=1~Nの離散周波数パ ラメータで表されるものとする。 付録3.3-1 時間波形の振幅とスペクトル J N A x(n)
2
/
)
(
2 1 0 2A
J
n
x
N n
掃引正弦波のエネルギは と、表される。また、パーシバルの等式を用いれば、
1 0 2 2 1 0 2)
(
1
2
/
)
(
N k N nk
X
N
A
J
n
x
TSP の場合、|X(k)|2=1 であるので、TSP 波形の振幅は、J
A
2
DFTの定義式が 時間信号 x(n) と、それをDFTした周波数成分 X(k) との間には、下記のパーシバル(Parseval: パーセ バルと表記される場合もある)の等式が成立する。
N
kn
j
k
X
N
n
x
N
kn
j
n
x
k
X
N k N n/
2
exp
)
(
1
)
(
/
2
exp
)
(
)
(
1 0 1 0
1 0 2 1 0 2)
(
1
)
(
N k N nk
X
N
n
x
パーシバルの等式 の場合(通常の場合) DFTの定義式が
N
kn
j
k
X
N
n
x
N
kn
j
n
x
N
k
X
N k N n/
2
exp
)
(
1
)
(
/
2
exp
)
(
1
)
(
1 0 1 0
1 0 2 1 0 2)
(
)
(
N k N nk
X
n
x
の場合(対称性を保った場合) よって、
1 0 2)
(
1
2
N kk
X
N
J
A
付録5.1-1 測定信号による雑音抑圧効果の計算 これらの関係を用いると、雑音抑圧効果 NRP は、
)
4
(
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
)
(
ˆ
)
(
ˆ
)
(
ˆ
)
(
1
)
(
ˆ
)
(
1
)
(
ˆ
)
(
1
)
(
ˆ
)
(
1
)
(
)
(
)
(
)
(
1
)
(
1
1 0 2 0 1 0 2 0 1 0 0 1 0 2 0 1 0 0 1 0 2 1 0 2 0 1 0 0 1 0 2 1 0 2 1 0 2 0 1 0 0 1 0 0 1 0 0 1 0 2 0 1 0 0 1 0 1 2 1 0 0 2 1 0
N k N S N k N N k N S N k N N k N N n N k N N k N N k N k N k N N k N N k N N k N N k N N k N N k N N k N N Nk
S
k
P
N
p
J
k
S
k
P
k
P
p
J
k
S
k
P
k
P
n
s
k
S
k
P
k
P
k
S
N
k
S
k
S
N
k
P
k
P
N
k
P
k
P
N
k
S
k
P
k
P
k
P
N
k
P
N
p
p
NRP
時間信号x(n)と、それをN点DFTした周波数スペクトルX(k)との間には、次式で示されるパーシバルの関係 (→付録3.3-1)が成立する。左辺は時間信号のエネルギーを表し、右辺は周波数スペクトルのエネルギーを 表す。)
1
(
)
(
1
)
(
1 0 2 1 0 2
N k N nk
X
N
n
x
ただし、 第1の等号は、定義)
3
(
)
(
1
)
(
)
(
ˆ
)
(
1
)
(
)
(
ˆ
1 0 2 2 2 1 0 0 0 0
N k N k N N Nk
S
N
k
S
k
S
k
P
N
k
P
k
P
ここで、定常雑音 ni(n) (i=0,1) のパワーを次式で定義し、式(1)を用いてて計算すると、
(
)
1
(
)
(
2
)
1
)
(
1
)
(
1
1 0 2 1 0 2 2 1 0 2 2 1 0 2
N k Ni N k i N k i N n i NiP
k
N
k
N
E
N
k
N
N
E
n
n
N
E
p
ただし、PNi(k)=E[|Ni(k)|2] である。また、P N0(k) と測定信号のパワースぺlクトル|S(k)|2をそれぞれのエネル ギーで正規化したものを 、P^N0(k)、|S^(k)|2と表し、次式で定義する測定信号のスペクトルは、