(57)【要約】【課題】スピーカ・マイクロホンアレイ間の音響伝達特性の測定に要する時間を短縮することができる音響伝達特性測定装置を提供する。【解決手段】スピーカアレイと、マイクロホンアレイと、ＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を計算する時空間周波数領域ＳＴＳＰ信号計算部と、時空間周波数領域ＳＴＳＰ信号に基づいて、ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算してスピーカアレイから出力する時空間領域ＳＴＳＰ信号計算部と、時空間周波数領域逆ＳＴＳＰ信号を計算する時空間周波数領域逆ＳＴＳＰ信号計算部と、時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する時空間領域逆ＳＴＳＰ信号計算部と、マイクロホンアレイに収音される信号である入力信号と、時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出部を含む。【選択図】図１
2015142248.tif 000002
【特許請求の範囲】
【請求項１】
複数のスピーカを直線状に配置してなるスピーカアレイと、
複数のマイクロホンを直線状に配置してなるマイクロホンアレイと、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力するＳＴＳＰ信号計算部と、
前記時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する逆ＳＴＳＰ信号計算部と、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出部と、
を含む音響伝達特性測定装置。
【請求項２】
請求項１に記載の音響伝達特性測定装置であって、
Δωを時間周波数の増加量を調整するパラメータとし、Δkxを空間周波数の増加量を調整するパラメータとし、ω0を時間周波数の初期値を調整するパラメータとし、kx0を空間周波数の初期値を調整するパラメータとし、i,jは空間方向（スピーカ方向）xiおよび時間方向tjのインデックスを示すものとし、Nxをスピーカ・マイクロホン数とし、Ntを時間方向の信号長とし、
前記ＳＴＳＰ信号計算部が、時空間領域ＳＴＳＰ信号を
【数１２】
2015142248.tif 000014
として計算し、
前記逆ＳＴＳＰ信号計算部が、時空間領域逆ＳＴＳＰ信号を
【数１３】
2015142248.tif 000015
として計算する
音響伝達特性測定装置。
【請求項３】
複数のスピーカを直線状に配置してなるスピーカアレイと、
複数のマイクロホンを直線状に配置してなるマイクロホンアレイと、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する時空間周波数領域ＳＴＳＰ信号計算部と、
前記時空間周波数領域ＳＴＳＰ信号に基づいて、前記ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力する時空間領域ＳＴＳＰ信号計算部と、
前記時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する時空間周波数領域逆ＳＴＳＰ信号計算部と、
前記時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する時空間領域逆ＳＴＳＰ信号計算部と、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出部と、
を含む音響伝達特性測定装置。
【請求項４】
請求項３に記載の音響伝達特性測定装置であって、
前記時空間周波数領域ＳＴＳＰ信号計算部が、
kx,nを空間周波数とし、ωmを時間周波数とし、nを空間周波数kx,nのインデクスとし、mを時間周波数ωmのインデクスとし、Δωを時間周波数の増加量を調整するパラメータとし、Δkxを空間周波数の増加量を調整するパラメータとし、ω0を時間周波数の初期値を調整するパラメータとし、kx0を空間周波数の初期値を調整するパラメータとし、前記時空間周波数領域ＳＴＳＰ信号を、
【数１４】
2015142248.tif 000016
として計算する
音響伝達特性測定装置。
【請求項５】
請求項１から４の何れかに記載の音響伝達特性測定装置であって、
前記伝達特性推定値算出部が、
２次元円状畳み込みを用いる方法、２次元線状畳み込みを用いる方法、円状畳み込みと線状畳み込みを組み合わせる方法、のいずれかにより時空間音響伝達特性の推定値を算出する
音響伝達特性測定装置。
【請求項６】
複数のスピーカを直線状に配置してなるスピーカアレイと、複数のマイクロホンを直線状に配置してなるマイクロホンアレイとを用いる音響伝達特性測定方法であって、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力するＳＴＳＰ信号計算ステップと、
前記時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する逆ＳＴＳＰ信号計算ステップと、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出ステップと、
を含む音響伝達特性測定方法。
【請求項７】
複数のスピーカを直線状に配置してなるスピーカアレイと、複数のマイクロホンを直線状に配置してなるマイクロホンアレイとを用いる音響伝達特性測定方法であって、
合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する時空間周波数領域ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域ＳＴＳＰ信号に基づいて、前記ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算して前記スピーカアレイから出力する時空間領域ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する時空間周波数領域逆ＳＴＳＰ信号計算ステップと、
前記時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する時空間領域逆ＳＴＳＰ信号計算ステップと、
前記マイクロホンアレイに収音される信号である入力信号と、前記時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する伝達特性推定値算出ステップと、
を含む音響伝達特性測定方法。
【請求項８】
コンピュータを、請求項１から５の何れかに記載の音響伝達特性測定装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、スピーカアレイとマイクロホンアレイを用いて、それらの間の音響伝達特性を測定する技術に関する。
【背景技術】
【０００２】
音響伝達特性を測定する従来技術にTime Stretched Pulse（ＴＳＰ）信号を用いるものがある（非特許文献１）。ＴＳＰ信号とは、パルス信号を時間方向に伸長した信号である。非特許文献１の音響伝達特性の測定技術は、ＴＳＰ信号を１個のスピーカから出力し、１個のマイクロホンで収音することで、単一のスピーカと単一のマイクロホンの間の伝達特性を、高精度かつ頑健に測定することを可能にする技術である。この技術を多数のスピーカ・マイクロホン間で適用する場合には、各スピーカからＴＳＰ信号を出力し、マイクロホンアレイで観測した信号に対し、逆ＴＳＰ信号を畳み込むことにより、音響伝達特性を算出する。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Y. Suzuki, F. Asano, H. Kim, and T. Sone. “An optimal computer-generated pulse signal suitable for the the measurement of very long impulse responses,” J. Acoust. Soc. Am., vol. 97, no. 2, pp.1119-1123, 1995.
【発明の概要】
【発明が解決しようとする課題】
【０００４】
従来技術では、多数のスピーカ、マイクロホン間の音響伝達特性を測定する場合において、各スピーカからＴＳＰ信号を出力する必要があるため、全スピーカからの特性を測定するためには長い時間を要するという課題があった。そこで本発明では、スピーカ・マイクロホンアレイ間の音響伝達特性の測定に要する時間を短縮することができる音響伝達特性測定装置を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の音響伝達特性測定装置は、スピーカアレイと、マイクロホンアレイと、ＳＴＳＰ信号計算部と、逆ＳＴＳＰ信号計算部と、伝達特性推定値算出部を含む。
【０００６】
スピーカアレイは、複数のスピーカを直線状に配置してなる。マイクロホンアレイは、複数のマイクロホンを直線状に配置してなる。ＳＴＳＰ信号計算部は、合成時の平面波の角度と周波数とが時間方向に変化する信号であるＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算してスピーカアレイから出力する。逆ＳＴＳＰ信号計算部は、時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する。伝達特性推定値算出部は、マイクロホンアレイに収音される信号である入力信号と、時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する。
【発明の効果】
【０００７】
本発明の音響伝達特性測定装置によれば、スピーカ・マイクロホンアレイ間の音響伝達特性の測定に要する時間を短縮することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の実施例１の音響伝達特性測定装置の構成を示すブロック図。
【図２】本発明の実施例１の音響伝達特性測定装置の動作を示すフローチャート。
【図３】本発明の実施例２の音響伝達特性測定装置の構成を示すブロック図。
【図４】本発明の実施例２の音響伝達特性測定装置の動作を示すフローチャート。
【発明を実施するための形態】
【０００９】
以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。
【実施例１】
【００１０】
＜Space-Time Stretched Pulse（ＳＴＳＰ）信号について＞
以下、本発明で用いられるSpace-Time Stretched Pulse（ＳＴＳＰ）信号について定義する。ＳＴＳＰ信号とは、簡単にはＴＳＰ信号を空間方向に伸長させた信号のことである。従って、ＳＴＳＰ信号は、パルス信号を時間方向および空間方向に伸長させた信号であり、合成時の平面波の角度と周波数とが時間方向に変化する信号であるとも表現できる。
【００１１】
以下、図１、図２を参照して実施例１の音響伝達特性測定装置について説明する。図１は、本実施例の音響伝達特性測定装置１の構成を示すブロック図である。図２は、本実施例の音響伝達特性測定装置１の動作を示すフローチャートである。
【００１２】
図１に示すように、本実施例の音響伝達特性測定装置１は、ＳＴＳＰ信号計算部１１と、スピーカアレイ１２と、逆ＳＴＳＰ信号計算部１３と、マイクロホンアレイ１４と、伝達特性推定値算出部１５を含む。ＳＴＳＰ信号計算部１１は、時空間周波数領域ＳＴＳＰ信号計算部１１１と、時空間領域ＳＴＳＰ信号計算部１１２を含む。逆ＳＴＳＰ信号計算部１３は、時空間周波数領域逆ＳＴＳＰ信号計算部１３１と、時空間領域逆ＳＴＳＰ信号計算部１３２を含む。
【００１３】
スピーカアレイ１２は、複数のスピーカ１２−１，１２−２，…，１２−Ｎｘを直線状に配置してなる。マイクロホンアレイ１４は、複数のマイクロホン１４−１，１４−２，…，１４−Ｎｘを直線状に配置してなる。ただし、Nxをスピーカ・マイクロホン数とし、Nxは２以上の任意の整数とする。
【００１４】
＜ＳＴＳＰ信号計算部１１＞
以下、ＳＴＳＰ信号計算部１１について説明する。ＳＴＳＰ信号計算部１１の時空間周波数領域ＳＴＳＰ信号計算部１１１は、前述のＳＴＳＰ信号の時空間周波数領域の表現である時空間周波数領域ＳＴＳＰ信号を、時間周波数と空間周波数とが線形変化する条件で計算する（Ｓ１１１）。より詳細には、時空間周波数領域ＳＴＳＰ信号計算部１１１は、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)を入力とし、時空間周波数領域ＳＴＳＰ信号F(kx,n,ωm)を、例えば以下の式（１）で計算する。
【００１５】
【数１】
2015142248.tif 000003
【００１６】
ここで、kx,nは空間周波数、ωmは時間周波数、Δωは時間周波数の増加量を調整するパラメータ、Δkxは空間周波数の増加量を調整するパラメータ、ω0は時間周波数の初期値を調整するパラメータ、kx0は空間周波数の初期値を調整するパラメータ、n,mはそれぞれ空間周波数kx,n，時間周波数ωmのインデックスを示しており、1≦n≦Nkx,1≦m≦Nωである。NkxとNωはそれぞれ、空間周波数ビン数、時間周波数ビン数である。空間周波数kx,n，時間周波数ωmは例えば以下のように定義される。
【００１７】
【数２】
2015142248.tif 000004
【００１８】
ここで、Δx,Δtはそれぞれ空間（スピーカ・マイクロホン）および時間サンプリング間隔である。次に、ＳＴＳＰ信号計算部１１の時空間領域ＳＴＳＰ信号計算部１１２は、時空間周波数領域ＳＴＳＰ信号に基づいて、ＳＴＳＰ信号の時空間領域の表現である時空間領域ＳＴＳＰ信号を計算してスピーカアレイ１２から出力する（Ｓ１１２）。より詳細には、時空間領域ＳＴＳＰ信号計算部１１２は、前述した時空間周波数領域ＳＴＳＰ信号F(kx,n,ωm)に対して、例えば２次元逆離散フーリエ変換を適用することにより、時空間領域ＳＴＳＰ信号f(xi,tj)を計算する。時空間領域ＳＴＳＰ信号計算部１１２は、例えば以下の式（２）で定義される２次元逆離散フーリエ変換を実行して、時空間領域ＳＴＳＰ信号f(xi,tj)を計算する。
【００１９】
【数３】
2015142248.tif 000005
【００２０】
ただし、i,jは空間方向（スピーカ方向）xiおよび時間方向tjのインデックスを示す。前述したように、Nxをスピーカ・マイクロホン数とし、Ntを時間方向の信号長とすると、インデクスi,jは1≦i≦Nx，1≦j≦Ntである。時空間領域ＳＴＳＰ信号計算部１１２は、この時空間領域ＳＴＳＰ信号f(xi,tj)をスピーカアレイ１２を構成する各スピーカ１２−１，１２−２，…，１２−Ｎｘから出力する。上述した、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)、及びその他のパラメータの具体的な値として、例えばΔω=400, ω0=0, Δkx=4000, kx0=0, Nx=512, Nt=512, Δx=0.06m, Δt=0.625ms（サンプリング周波数16kHz）を用いることができる。
【００２１】
＜逆ＳＴＳＰ信号計算部１３＞
以下、逆ＳＴＳＰ信号計算部１３について説明する。逆ＳＴＳＰ信号計算部１３の時空間周波数領域逆ＳＴＳＰ信号計算部１３１は、時空間周波数領域ＳＴＳＰ信号に畳み込むことでパルス信号を生成する信号である時空間周波数領域逆ＳＴＳＰ信号を計算する（Ｓ１３１）。より詳細には、時空間周波数領域逆ＳＴＳＰ信号計算部１３１は、ＳＴＳＰ信号設計パラメータ(Δω,ω0,Δkx,kx0)を入力とし、例えば、以下の式（３）に基づき、時空間周波数領域逆ＳＴＳＰ信号G(kx,n,ωm)を計算する。
【００２２】
【数４】
2015142248.tif 000006
【００２３】
逆ＳＴＳＰ信号計算部１３の時空間領域逆ＳＴＳＰ信号計算部１３２は、時空間周波数領域逆ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を計算する（Ｓ１３２）。より詳細には、時空間領域逆ＳＴＳＰ信号計算部１３２は、時空間周波数領域逆ＳＴＳＰ信号G(kx,n,ωm)に対して、例えば２次元逆離散フーリエ変換を適用することにより、時空間領域逆ＳＴＳＰ信号g(xi,tj)を計算する。時空間領域逆ＳＴＳＰ信号計算部１３２は、例えば以下の式（４）で定義される２次元逆離散フーリエ変換を実行して、時空間領域逆ＳＴＳＰ信号g(xi,tj)を計算する。
【００２４】
【数５】
2015142248.tif 000007
【００２５】
＜伝達特性推定値算出部１５＞
以下、伝達特性推定値算出部１５について説明する。伝達特性推定値算出部１５は、マイクロホンアレイに収音される信号である入力信号と、時空間領域逆ＳＴＳＰ信号に基づいて時空間音響伝達特性の推定値を算出する（Ｓ１５）。より詳細には、マイクロホンアレイ１４に収音される信号である入力信号（マイクロホンアレイ入力信号ともいう）であるf~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)を入力とし、時空間音響伝達特性の推定値h~(xi,tj)（または、時空間周波数領域の表現であるH~(kx,n,ωm)）を出力とする。伝達特性推定値算出部１５が出力する時空間音響伝達特性の推定値は、時空間領域の表現であってもよいし、時空間周波数領域の表現であってもよい。例えば、伝達特性推定値算出部１５は、時空間音響伝達特性の推定値をマイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)の畳み込みによって算出する。伝達特性推定値算出部１５は、例えば二次元の円状畳み込みを用いる方法、または二次元の線状畳み込みを用いる方法、または円状畳み込みと線状畳み込みを組み合わせる方法によって、時空間音響伝達特性の推定値を算出できる。以下に、これら３つの算出方法の一例を示す。
【００２６】
＜二次元円状畳み込みを用いる方法＞
二次元円状畳み込みを用いる場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)の２次元離散フーリエ変換であるF~(kx,n,ωm)と、時空間領域逆ＳＴＳＰ信号g(xi,tj)の２次元離散フーリエ変換であるG(kx,n,ωm)を計算する。伝達特性推定値算出部１５は、例えば以下の式（５）で定義される２次元離散フーリエ変換を実行してF~(kx,n,ωm)を計算する。
【００２７】
【数６】
2015142248.tif 000008
【００２８】
これらを乗算した信号が、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)となる。従って、伝達特性推定値算出部１５は、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)を以下の式（６）で計算する。
【００２９】
【数７】
2015142248.tif 000009
【００３０】
伝達特性推定値算出部１５は、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を、H~(kx,n,ωm)の２次元逆離散フーリエ変換として算出できる。なお、伝達特性推定値算出部１５は、H~(kx,n,ωm)をそのまま出力しても構わない。
【００３１】
＜２次元線状畳み込みを用いる方法＞
２次元線状畳み込みを用いる場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)に対し、以下のように２次元線状畳み込みを行うことで、直接、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を算出する。従って、伝達特性推定値算出部１５は、以下の式（７）により、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)を算出する。
【００３２】
【数８】
2015142248.tif 000010
【００３３】
χ、τはそれぞれ、空間方向、時間方向の遅延を表す。
【００３４】
＜円状畳み込みと線状畳み込みの組み合わせ＞
ここでは、例えば空間方向xiで線状畳み込みを行い、時間方向tjで円状畳み込みを行う場合を考える。この場合、伝達特性推定値算出部１５は、マイクロホンアレイ入力信号f~(xi,tj)と時空間領域逆ＳＴＳＰ信号g(xi,tj)に対して、xi方向の両端にNx/2サンプルだけゼロづめをした信号を、それぞれf−(xi,tj),g−(xi,tj)と設定する。ただしここでは、1≦i≦2Nx,1≦j≦Ntである。伝達特性推定値算出部１５は、f−(xi,tj),g−(xi,tj)の２次元離散フーリエ変換であるF−(kx,n,ωm),G−(kx,n,ωm)を計算する。これらを乗算した信号が、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)となる。従って、伝達特性推定値算出部１５は、以下の式（８）により、時空間周波数領域での時空間音響伝達特性の推定値H~(kx,n,ωm)を計算する。
【００３５】
【数９】
2015142248.tif 000011
【００３６】
伝達特性推定値算出部１５は、時空間領域での時空間音響伝達特性の推定値h~(xi,tj)は、H~(kx,n,ωm)の２次元逆離散フーリエ変換として算出できる。なお、伝達特性推定値算出部１５は、H~(kx,n,ωm)をそのまま出力しても構わない。ただし、この方法ではゼロづめした部分に相当する信号は、測定値としては扱わないものとする。
【００３７】
なお、上述では、周波数領域−時間領域の変換の具体的方法として離散フーリエ変換（および逆変換）を用いた。例えばステップＳ１１２，Ｓ１３２，Ｓ１５である。しかしながら、本発明において周波数領域−時間領域の変換は、離散フーリエ変換（逆変換）には限定されず、ステップＳ１１２，Ｓ１３２，Ｓ１５では、他の方法により、周波数領域−時間領域の変換を実行してもよい。
【００３８】
本実施例の音響伝達特性測定装置１によれば、ＳＴＳＰ信号計算部１１がパルス信号を時間方向だけでなく空間方向にも伸長させたＳＴＳＰ信号を直線状に配置されたスピーカアレイ１２から出力し、この信号を直線状に配置されたマイクロホンアレイ１４で収音し、伝達特性推定値算出部１５が時空間音響伝達特性の推定値をマイクロホンアレイ入力信号と時空間領域逆ＳＴＳＰ信号に基づいて算出するため、一度の測定で伝達特性を求めることが可能となり、時空間音響伝達特性の測定に要する時間を大幅に短縮することができる。
【００３９】
なお、本実施例の音響伝達特性測定装置１は、波面合成法を用い、マイクロホンアレイで収音された音を遠隔地のスピーカアレイで忠実に再生することによって臨場感を実現する音声コミュニケーションにシステムに利用することができる。本実施例の音響伝達特性測定装置１は、波面合成法によるフィルタを生成する際に必要となる伝達特性を求めるのに有用である。
【実施例２】
【００４０】
以下、図３、図４を参照して実施例２の音響伝達特性測定装置について説明する。図３は、本実施例の音響伝達特性測定装置２の構成を示すブロック図である。図４は、本実施例の音響伝達特性測定装置２の動作を示すフローチャートである。
【００４１】
図３に示すように、本実施例の音響伝達特性測定装置２は、ＳＴＳＰ信号計算部２１と、スピーカアレイ１２と、逆ＳＴＳＰ信号計算部２３と、マイクロホンアレイ１４と、伝達特性推定値算出部１５を含む。
【００４２】
実施例１のＳＴＳＰ信号計算部１１、逆ＳＴＳＰ信号計算部１３は、本実施例においてＳＴＳＰ信号計算部２１、逆ＳＴＳＰ信号計算部２３に変更されている。このほかの構成は、実施例１と実施例２において共通するため説明を割愛する。
【００４３】
＜ＳＴＳＰ信号計算部２１＞
本実施例のＳＴＳＰ信号計算部２１は、時空間周波数領域での計算を行わずに、時空間領域ＳＴＳＰ信号を直接計算してスピーカアレイから出力する（Ｓ２１）。具体的には、ＳＴＳＰ信号計算部２１は、以下の式（９）により、時空間領域ＳＴＳＰ信号を直接計算する。
【００４４】
【数１０】
2015142248.tif 000012
【００４５】
＜逆ＳＴＳＰ信号計算部２３＞
本実施例の逆ＳＴＳＰ信号計算部２３は、時空間周波数領域での計算を行わずに、時空間領域ＳＴＳＰ信号に基づいて時空間領域逆ＳＴＳＰ信号を直接計算する（Ｓ２３）。より詳細には、逆ＳＴＳＰ信号計算部２３は、式（９）により計算した時空間領域ＳＴＳＰ信号f(xi,tj)を時空間方向に反転したものを時空間領域の逆ＳＴＳＰ信号とする。具体的には、逆ＳＴＳＰ信号計算部２３は、以下の式（１０）により、時空間領域逆ＳＴＳＰ信号を直接計算する。
【００４６】
【数１１】
2015142248.tif 000013
【００４７】
本実施例の音響伝達特性測定装置２によれば、時空間周波数領域の計算を行わずに、実施例１と同様の効果を得ることができる。  上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。
【００４８】
また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。
【００４９】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。
【００５０】
また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。
【００５１】
このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。
【００５２】
また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【図１】
2015142248.tif 000017
【図２】
2015142248.tif 000018
【図３】
2015142248.tif 000019
【図４】
2015142248.tif 000020