特許 7601231 音響特性計算装置、音響特性計算方法、プログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-12-09

(45)【発行日】2024-12-17

(54)【発明の名称】音響特性計算装置、音響特性計算方法、プログラム

(51)【国際特許分類】

   H04R 29/00 20060101AFI20241210BHJP

   G10K 15/00 20060101ALI20241210BHJP

【ＦＩ】

H04R29/00 310

G10K15/00 L

G10K15/00 K

【請求項の数】 4

(21)【出願番号】P 2023536253

(86)(22)【出願日】2021-07-20

(86)【国際出願番号】 JP2021027075

(87)【国際公開番号】W WO2023002554

(87)【国際公開日】2023-01-26

【審査請求日】2023-12-20

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(74)【代理人】

【識別番号】100121706

【弁理士】

【氏名又は名称】中尾 直樹

(74)【代理人】

【識別番号】100128705

【弁理士】

【氏名又は名称】中村 幸雄

(74)【代理人】

【識別番号】100147773

【弁理士】

【氏名又は名称】義村 宗洋

(72)【発明者】

【氏名】石川 憲治

(72)【発明者】

【氏名】白木 善史

(72)【発明者】

【氏名】守谷 健弘

【審査官】▲徳▼田 賢二

(56)【参考文献】

【文献】特開２０１１－２１１３１１（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｈ０４Ｒ ２９／００

Ｇ１０Ｋ １５／００

(57)【特許請求の範囲】

【請求項1】

Lを略円形の音源であるトランスデューサを用いて生じさせたパラメトリックアレイによって発生した可聴音（以下、復調音という）Sによる光の位相変化を測定する音場測定装置における光路、Cを光路Lとトランスデューサを通り復調音Sの伝搬方向と平行な直線とが交わる点、Xを点Cからの距離がxである光路L上の点、zをトランスデューサから光路Lまでの距離、aをトランスデューサの半径、f_d=|f₁-f₂|（ただし、f₁は搬送波の周波数、f₂は側帯波の周波数）を復調音Sの周波数、cを音速とし、
復調音Sによる光位相の変化量φ_sから、周波数f_dにおける光位相の複素振幅dを計算する第１計算部と、
q_diff(ξ’,η’)をガウシアンビーム展開法を用いて定義される関数とし、
点Xにおける関数値q_diff(ξ,η)（ただし、ξ=x/a, η=2z/ka², k₁=f₁/c, k₂=f₂/c, k=(k₁+k₂)/2）と、光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’とを計算する第２計算部と、
周波数f_dにおける光位相の複素振幅dと点Xにおける関数値q_diff(ξ,η)と光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’を用いて、点Xにおける周波数f_dの復調音の複素振幅pを計算する第３計算部と、
を含む音響特性計算装置。

【請求項2】

請求項１に記載の音響特性計算装置であって、
関数q_diff(ξ’,η’)は、次式で定義される関数であり、

【数16】

【数17】

（ただし、r₁ ^(mm’)=(k_aB_{k_bm’}+k_bB_{k_am})+i(k_a+k_b)η’B_{k_am}B_{k_bm’}, r₂ ^(mm’)=(k_aB_{k_am}+k_bB_{k_bm’})η’-i(k_a+k_b), s₁ ^(mm’)=(k_a+k_b)²B_{k_am}B_{k_bm’}, s₂ ^(mm’)=-i(k_a+k_b)k_ak_b(B_{k_am}-B_{k_bm’})², B_{k_am}=B_m ⁽¹⁾/k_a+i/δ, B_{k_bm’}=B_m’ ⁽²⁾/k_b+i/δ, k_a=k₁/k, k_b=-k₂/kであり、A_m, A_m’, B_m ⁽¹⁾, B_m’ ⁽²⁾, δは所定の定数）
前記第３計算部は、次式により、複素振幅pを計算する

【数18】

（ただし、C_aoは音場測定時の物理的条件から定まる値）
ことを特徴とする音響特性計算装置。

【請求項3】

Lを略円形の音源であるトランスデューサを用いて生じさせたパラメトリックアレイによって発生した可聴音（以下、復調音という）Sによる光の位相変化を測定する音場測定装置における光路、Cを光路Lとトランスデューサを通り復調音Sの伝搬方向と平行な直線とが交わる点、Xを点Cからの距離がxである光路L上の点、zをトランスデューサから光路Lまでの距離、aをトランスデューサの半径、f_d=|f₁-f₂|（ただし、f₁は搬送波の周波数、f₂は側帯波の周波数）を復調音Sの周波数、cを音速とし、
音響特性計算装置が、復調音Sによる光位相の変化量φ_sから、周波数f_dにおける光位相の複素振幅dを計算する第１計算ステップと、
q_diff(ξ’,η’)をガウシアンビーム展開法を用いて定義される関数とし、
前記音響特性計算装置が、点Xにおける関数値q_diff(ξ,η)（ただし、ξ=x/a, η=2z/ka², k₁=f₁/c, k₂=f₂/c, k=(k₁+k₂)/2）と、光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’とを計算する第２計算ステップと、
前記音響特性計算装置が、周波数f_dにおける光位相の複素振幅dと点Xにおける関数値q_diff(ξ,η)と光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’を用いて、点Xにおける周波数f_dの復調音の複素振幅pを計算する第３計算ステップと、
を含む音響特性計算方法。

【請求項4】

請求項１または２に記載の音響特性計算装置としてコンピュータを機能させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、パラメトリックアレイの音響特性を測定する技術に関する。

【背景技術】

【0002】

パラメトリックアレイは、媒質の非線形性によって空間中に音が生じる非線形音響現象である（非特許文献１参照）。より詳しくは、パラメトリックアレイとは、２つの異なる周波数を有する音が伝搬するとき、当該周波数の差で表される周波数を有する音が非線形的に生じる現象のことである。この現象を利用して可聴音を再生するパラメトリックアレイスピーカーは、限られた空間への音の再生や空間音響表現技術、騒音制御などに用いられる。

【0003】

マイクロホンを用いて、パラメトリックアレイにより生じた可聴音（以下、復調音という）を測定する場合、擬音と呼ばれる雑音が生じることが知られている。擬音は、パラメトリックアレイ中の大振幅超音波によって生じるマイクロホンを含む受音系の非線形雑音である。原理的に擬音の周波数が復調音のそれと一致することから、測定対象である復調音と擬音が重畳され両者は区別できない形で検出される。

【先行技術文献】

【非特許文献】

【0004】

【文献】Woon-Seng Gan, Jun Yang, and Tomoo Kamakura, “A review of parametric acoustic array in air,” Applied Acoustics, Vol.73, Issue 12, pp.1211-1219, 2012.

【発明の概要】

【発明が解決しようとする課題】

【0005】

復調音の音響特性を測定するため、擬音の原因となるマイクロホンに入射する超音波成分を空間中で減衰させ復調音成分を透過させる音響フィルタを用いて擬音を除去する方法が考えられる。しかし、音響フィルタによる擬音の除去効果は限定的であり、完全に擬音を除去することはできない。つまり、復調音の音響特性を正確に測定することができない。

【0006】

そこで本発明では、正確にパラメトリックアレイの音響特性を測定する技術を提供することを目的とする。

【課題を解決するための手段】

【0007】

本発明の一態様は、Lを略円形の音源であるトランスデューサを用いて生じさせたパラメトリックアレイによって発生した可聴音（以下、復調音という）Sによる光の位相変化を測定する音場測定装置における光路、Cを光路Lとトランスデューサを通り復調音Sの伝搬方向と平行な直線とが交わる点、Xを点Cからの距離がxである光路L上の点、zをトランスデューサから光路Lまでの距離、aをトランスデューサの半径、f_d=|f₁-f₂|（ただし、f₁は搬送波の周波数、f₂は側帯波の周波数）を復調音Sの周波数、cを音速とし、復調音Sによる光位相の変化量φ_sから、周波数f_dにおける光位相の複素振幅dを計算する第１計算部と、q_diff(ξ’,η’)をガウシアンビーム展開法を用いて定義される関数とし、点Xにおける関数値q_diff(ξ,η)（ただし、ξ=x/a, η=2z/ka², k₁=f₁/c, k₂=f₂/c, k=(k₁+k₂)/2）と、光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’とを計算する第２計算部と、周波数f_dにおける光位相の複素振幅dと点Xにおける関数値q_diff(ξ,η)と光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’を用いて、点Xにおける周波数f_dの復調音の複素振幅pを計算する第３計算部と、を含む。

【発明の効果】

【0008】

本発明によれば、正確にパラメトリックアレイの音響特性を測定することが可能となる。

【図面の簡単な説明】

【0009】

【図1】音場測定装置８００の構成を示すブロック図である。

【図2】音響特性計算装置１００の構成を示すブロック図である。

【図3】音響特性計算装置１００の動作を示すフローチャートである。

【図4】点Xの位置を示す図である。

【図5】本発明の実施形態における各装置を実現するコンピュータの機能構成の一例を示す図である。

【発明を実施するための形態】

【0010】

以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。

【0011】

各実施形態の説明に先立って、この明細書における表記方法について説明する。

【0012】

^（キャレット）は上付き添字を表す。例えば、x^{y^z}はy^zがxに対する上付き添字であり、x_y^zはy^zがxに対する下付き添字であることを表す。また、_（アンダースコア）は下付き添字を表す。例えば、x^y_zはy_zがxに対する上付き添字であり、x_{y_z}はy_zがxに対する下付き添字であることを表す。

【0013】

ある文字xに対する^xや~xのような上付き添え字の”^”や”~”は、本来”x”の真上に記載されるべきであるが、明細書の記載表記の制約上、^xや~xと記載しているものである。
＜技術的背景＞
本発明の実施形態では、光による音計測技術（参考非特許文献１参照）を用いて、非接触にパラメトリックアレイを測定することで、原理的に擬音の存在することのない測定を実現する。［背景技術］で説明したように、擬音はマイクロホンの振動膜の非線形応答に起因するものである。光による音計測技術では振動膜に相当するものが存在しないため、原理的に擬音が発生することはない。
（参考非特許文献１：Kenji Ishikawa, Kohei Yatabe, Nachanant Chitanont, Yusuke Ikeda, Yasuhiro Oikawa, Takashi Onuma, Hayato Niwa, and Minoru Yoshii, “High-speed imaging of sound using parallel phase-shifting interferometry,” Optics Express, Vol.24, Issue 12, pp.12922-12932, 2016.）
まず、音響光学効果と呼ばれる音による媒質の屈折率変化を利用した音場測定法について説明する。音響光学効果によると、音による光位相の変化量φ_sは、次式により表される。

【0014】

【数1】

【0015】

ただし、k_lは光の波数、n₀は定常状態の空気屈折率、P₀は定常状態の大気圧、γは空気の比熱比である。式(1)の積分は光の伝搬経路（以下、光路と呼ぶ）Lに沿った音圧pの線積分を表す。ここでは、x軸を光路Lと平行となるように定義している。また、音圧pはある時刻における音圧である。なお、k_l, n₀, P₀, γは測定時の物理的条件から定まる定数である。

【0016】

式(1)からわかるように、音による光位相の変化量φ_sは、空間内の一点における音圧を表す量ではなく、光路Lに沿った音圧の線積分値に比例する量である。

【0017】

本発明の実施形態では、ガウシアンビーム展開法（GBE: Gaussian Beam Expansion)を用いることにより、測定値である光位相の変化量φ_sから空間内の一点における復調音の複素振幅pを計算する。具体的には、復調音を生じさせる原因である搬送波と側帯波をそれぞれガウシアンビームの和で近似する。ここで、ガウシアンビームとは、音圧を一定時間区間でフーリエ変換することにより得られる複素振幅がガウス関数で表現される波のことである。また、ピストン運動をする音源から生じる波については、ガウシアンビームの和として近似することができる。

【0018】

以下、パラメトリックアレイを発生させる音源として、円形に複数のマイクロホンを並べて構成されたパラメトリックアレイスピーカーのような、円形の音源を仮定する。音源は、完全な円でなく、略円形の音源であってもよい。本発明の実施形態で扱うピストン運動をする音源から生じる波は、パラメトリックアレイスピーカーを構成するマイクロホンがすべて同じ音を出すことで生じるものである。

【0019】

周波数f₁の搬送波、周波数f₂の側帯波の相互作用により生じる周波数f_d=|f₁-f₂|の復調音の複素振幅をpとする。参考非特許文献２によれば、円形の音源によって生成される復調音の複素振幅pは、ガウシアンビーム展開法を用いて次式のように表すことができる。

【0020】

【数2】

【0021】

ただし、ξ=x/a, η=2z/ka², k=(k₁+k₂)/2, k_d=k₁-k₂であり、aは音源の半径、k₁は搬送波の波数、k₂は側帯波の波数、p₀ ⁽¹⁾は搬送波の振幅、p₀ ⁽²⁾は側帯波の振幅、βは空気の非線形係数、ρは空気の密度、cは音速、iは虚数単位を表す。また、z軸を復調音の伝搬方向と平行となるように定義している。なお、波数k₁, k₂はそれぞれk₁=f₁/c, k₂=f₂/cと表すことができる。したがって、k_d=f_d/cである。

【0022】

また、q_diff(ξ, η)は復調音のビーム形状を表す項であり、

【0023】

【数3】

【0024】

【数4】

【0025】

と定義される。ここで、mは搬送波を近似するために用いるガウシアンビームを表すパラメータ、m’は側帯波を近似するために用いるガウシアンビームを表すパラメータである。式(3)からわかるように、搬送波、側帯波をそれぞれ10のガウシアンビームを用いて近似している。これは、10のガウシアンビームを用いると十分な近似精度が得られることが経験的に知られていることに基づくものである。また、A_m, A_m’ (m=1, …, 10, m’=1, …, 10)はそれぞれガウシアンビームの係数であり、例えば、参考非特許文献３の表１のA_nの値を用いることができる。また、E₁(x)は指数積分関数である。r₁ ^(mm’), r₂ ^(mm’), s₁ ^(mm’), s₂ ^(mm’)(m=1, …, 10, m’=1, …, 10)は、それぞれ次式で定義される（参考非特許文献２参照）。

【0026】

【数5】

【0027】

【数6】

【0028】

【数7】

【0029】

【数8】

【0030】

【数9】

【0031】

【数10】

【0032】

ここで、k_a=k₁/k, k_b=-k₂/kである。また、B_m ⁽¹⁾, B_m’ ⁽²⁾ (m=1, …, 10, m’=1, …, 10)は所定の定数であり、例えば、参考非特許文献３の表１のB_nの値を用いることができる。δは音源の正規化された曲率半径を表す。音源が平面状であればδは無限大となり、項i/δは0となる。

【0033】

つまり、式(2)は、搬送波と側帯波をガウシアンビームの和で近似した場合における搬送波と側帯波の相互作用により生じる復調音の複素振幅を表す式である。
（参考非特許文献２：Desheng Ding, “A simplified algorithm for the second-order sound fields,” The Journal of Acoustical Society of America, Vol.108, Issue 6, pp.2759-2764, 2000.）
（参考非特許文献３：J. J. Wen and M. A. Breazeale, “A diffraction beam field expressed as the superposition of Gaussian beams,” The Journal of Acoustical Society of America, Vol.83, Issue5, pp.1752-1756, 1988.）
式(2)を光路Lに沿って線積分することにより、次式の通り、周波数f_dにおける光位相の複素振幅dが得られる。

【0034】

【数11】

【0035】

ただし、C_ao=k_l(n₀-1)/γP₀である。

【0036】

式(4)の複素振幅dは、光位相の変化量φ_sの周波数スペクトルであり、光位相の変化量φ_sを一定の時間区間でフーリエ変換して、予め定めた周波数f_dの成分を取り出したものに相当する。

【0037】

式(2)、式(4)より、周波数f_dの復調音の複素振幅pは、次式により計算することができる。

【0038】

【数12】

【0039】

式(5)中の複素振幅dは測定結果として得られる光位相の変化量φ_sから求まる値、C_aoは測定時の物理的条件から定まる値である。また、関数q_diffの値とその積分値は、搬送波と側帯波をそれぞれ近似するために用いるガウシアンビームの数を定め、音源の半径aと曲率半径δを設定することにより、計算することができる。

【0040】

以上からわかるように、搬送波と側帯波がガウスビームの和で表すことができると仮定することにより、復調音の複素振幅pが位置と波数の関数（式(2)参照）として定義でき、空間の一点における復調音の複素振幅pが、測定結果として得られる光位相の変化量φ_sから近似的に求められるようになる。

【0041】

光による音計測技術により原理的に擬音の存在しない測定を実現することで、パラメトリックアレイの音響特性を正確に測定することが可能となる。擬音の影響は特にトランスデューサの近傍で顕著であるため、これまでトランスデューサの近傍におけるパラメトリックアレイの音響特性を測定することは困難であった。上記方法を用いて音場を測定することにより音源放射特性の詳細な計測が可能となるため、上記方法はトランスデューサの検査や高性能化開発に役立つ。また、上記方法によれば、設置場所の制約によってパラメトリックスピーカの近傍で聴取する状況において音圧レベルや再生音の品質を正確に測定する方法を実現することもできる。
＜第１実施形態＞
音響特性計算装置１００は、光を用いて音場を測定する音場測定装置８００の出力である音による光位相の変化量を入力とし、空間内の一点における復調音の複素振幅を出力する。そこで、まず、図１を参照して音場測定装置８００を説明する。

【0042】

図１は、音場測定装置８００の構成の一例を示す図である。図１に示すように音場測定装置８００は、トランスデューサ８１０と、光源８２０と、位相変化測定器８３０とを含む。ここで、トランスデューサ８１０とは、パラメトリックアレイを生じさせることにより復調音を発生させる略円形の音源のことであり、例えば、円形に複数のマイクロホンを並べて構成されたパラメトリックアレイスピーカーとして構成することができる。また、光源８２０と位相変化測定器８３０とを含む構成部を光測定装置という。光測定装置は、パラメトリックアレイの音響特性を求めるために用いる復調音による光位相の変化量を測定する。光測定装置としては、光位相の変化量を測定可能な装置であればどのようなものでもよく、例えば、参考非特許文献１や参考非特許文献４に記載のものを用いることができる。なお、光路Lは、図１に示すように、トランスデューサ８１０から放射される超音波（放射音）の伝播方向と直交するように設定する。
（参考非特許文献４：A. Torras-Rosell, S. Barrera-Figueroa, and F. Jacobsen,“Sound field reconstruction using acousto optic tomography,” The Journal of Acoustical Society of America, Vol.131, Issue 5, pp.3786-3793, 2012.）
以下、音場測定装置８００の動作について説明する。まず、トランスデューサ８１０は、復調音を発生させ、音場を生成する。次に、光源８２０は光を放射する。光源８２０から放射された光は、音場を透過することにより、音による位相変調を受ける。音による位相変調を受けた光は、位相変化測定器８３０に入力され、位相変化測定器８３０により当該位相変調の量に依存した光量変化を生じ、位相変化測定器８３０は、この変化した光量の分布、すなわち、復調音による光位相の変化量を出力する。

【0043】

以下、図２～図３を参照して音響特性計算装置１００を説明する。図２は、音響特性計算装置１００の構成を示すブロック図である。図３は、音響特性計算装置１００の動作を示すフローチャートである。図２に示すように音響特性計算装置１００は、第１計算部１１０と、第２計算部１２０と、第３計算部１３０と、記録部１９０を含む。記録部１９０は、音響特性計算装置１００の処理に必要な情報を適宜記録する構成部である。なお、音場測定装置８００と音響特性計算装置１００を含む装置のことを音響特性測定装置１０００という。

【0044】

以下では、Lを略円形の音源であるトランスデューサ８１０を用いて生じさせたパラメトリックアレイによって発生した復調音Sによる光の位相変化を測定する音場測定装置８００における光路、Cを光路Lとトランスデューサ８１０を通り復調音Sの伝搬方向と平行な直線とが交わる点、Xを点Cからの距離がxである光路L上の点、zをトランスデューサ８１０から光路Lまでの距離、aをトランスデューサ８１０の半径、f_d=|f₁-f₂|（ただし、f₁は搬送波の周波数、f₂は側帯波の周波数）を復調音Sの周波数、cを音速とする。図４は、点Xの位置を示す図である。

【0045】

図３に従い音響特性計算装置１００の動作について説明する。

【0046】

Ｓ１１０において、第１計算部１１０は、音場測定装置８００の出力である復調音Sによる光位相の変化量φ_sを入力とし、光位相の変化量φ_sから周波数f_dにおける光位相の複素振幅dを計算する。具体的には、第１計算部１１０は、変化量φ_sを一定時間区間でフーリエ変換して、周波数f_dの成分を取り出すことすることにより、複素振幅dを得る。

【0047】

Ｓ１２０において、第２計算部１２０は、次式で定義される関数q_diff(ξ’,η’)を用いて、点Xにおける関数値q_diff(ξ,η)（ただし、ξ=x/a, η=2z/ka², k₁=f₁/c, k₂=f₂/c, k=(k₁+k₂)/2）と、光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’とを計算する。

【0048】

【数13】

【0049】

【数14】

【0050】

（ただし、r₁ ^(mm’)=(k_aB_{k_bm’}+k_bB_{k_am})+i(k_a+k_b)η’B_{k_am}B_{k_bm’}, r₂ ^(mm’)=(k_aB_{k_am}+k_bB_{k_bm’})η’-i(k_a+k_b), s₁ ^(mm’)=(k_a+k_b)²B_{k_am}B_{k_bm’}, s₂ ^(mm’)=-i(k_a+k_b)k_ak_b(B_{k_am}-B_{k_bm’})², B_{k_am}=B_m ⁽¹⁾/k_a+i/δ, B_{k_bm’}=B_m’ ⁽²⁾/k_b+i/δ, k_a=k₁/k, k_b=-k₂/kであり、A_m, A_m’, B_m ⁽¹⁾, B_m’ ⁽²⁾, δは所定の定数）
つまり、関数q_diff(ξ’,η’)はガウシアンビーム展開法を用いて定義される関数である。

【0051】

なお、A_m, A_m’, B_m ⁽¹⁾, B_m’ ⁽²⁾ (m=1, …, 10, m’=1, …, 10), δは予め記録部１９０に記録しておくとよい。

【0052】

Ｓ１３０において、第３計算部１３０は、Ｓ１１０で計算した周波数f_dにおける光位相の複素振幅d、Ｓ１２０で計算した点Xにおける関数値q_diff(ξ,η)と光路Lに沿った関数q_diff(ξ’,η)の線積分値∫_Lq_diff(ξ’,η)dξ’とを入力とし、次式により、点Xにおける周波数f_dの復調音の複素振幅pを計算する。

【0053】

【数15】

【0054】

（ただし、C_aoは音場測定時の物理的条件から定まる値）
なお、C_ao=k_l(n₀-1)/γP₀であり、定数C_aoの計算に用いる光の波数k_l、定常状態の空気屈折率n₀、定常状態の大気圧P₀、空気の比熱比γは予め記録部１９０に記録しておくとよい。

【0055】

本発明の実施形態によれば、光による音計測技術を用いて得た復調音による光位相の変化量を用いることにより、正確にパラメトリックアレイの音響特性である復調音の複素振幅を求めることが可能となる。

【0056】

以下、本実施形態を応用した例について説明する。
（応用例１）
本実施形態を用いて、パラメトリックスピーカの周波数特性を測定する。測定対象とするパラメトリックスピーカを所望の位置に設置する。復調音の周波数を要素とするベクトル（以下、復調音周波数ベクトルという）から周波数を１つ取り出し、本実施形態によりある点での当該周波数における複素振幅を計算する。復調音周波数ベクトルのすべての要素に対して当該計算を繰り返すことにより、ある点での各周波数の復調音の複素振幅を要素とするベクトル（以下、復調音複素振幅ベクトルという）を生成する。これにより、パラメトリックスピーカの周波数特性を表す復調音周波数ベクトルと復調音複素振幅ベクトルの組が得られる。
（応用例２）
本実施形態を用いて、パラメトリックアレイの時間信号の音圧を測定する。音場測定装置８００によって測定された離散時間信号である光位相の変化量φ_sに対して離散フーリエ変換を行い、離散周波数スペクトルを得る。離散周波数スペクトルのすべての周波数における複素振幅に対して本実施形態によりある点における当該周波数の復調音の複素振幅を計算し、ある点における離散周波数スペクトルを生成する。生成されたある点における離散周波数スペクトルに対して逆フーリエ変換を行い、当該点における音圧の時間波形、すなわち、パラメトリックアレイの時間信号の音圧を得る。

【0057】

なお、離散時間信号φ_sに対して、予め設定されたサンプル数が得られ次第、離散フーリエ変換を行うようにしてもよい。当該サンプル数の離散時間信号に対して、上記処理を実行することにより、ある点における音圧の時間波形を得る。このとき、計算に時間がかかる第２計算部１２０での計算を予め実行しておくことで、ある点における音圧の時間波形をリアルタイムで得ることができるようになる。測定終了まで上記所定のサンプル数ごとの処理を繰り返すことにより、パラメトリックアレイの時間信号の音圧をリアルタイムに測定することができる。リアルタイムに得られる復調音の音圧の時間波形は、例えば、能動音場制御装置、オンライン診断装置、設置条件・駆動信号の最適化装置への入力として利用することができる。
＜補記＞
図５は、上述の各装置を実現するコンピュータ２０００の機能構成の一例を示す図である。上述の各装置における処理は、記録部２０２０に、コンピュータ２０００を上述の各装置として機能させるためのプログラムを読み込ませ、制御部２０１０、入力部２０３０、出力部２０４０などに動作させることで実施できる。

【0058】

本発明の装置は、例えば単一のハードウェアエンティティとして、キーボードなどが接続可能な入力部、液晶ディスプレイなどが接続可能な出力部、ハードウェアエンティティの外部に通信可能な通信装置（例えば通信ケーブル）が接続可能な通信部、ＣＰＵ（Central Processing Unit、キャッシュメモリやレジスタなどを備えていてもよい）、メモリであるＲＡＭやＲＯＭ、ハードディスクである外部記憶装置並びにこれらの入力部、出力部、通信部、ＣＰＵ、ＲＡＭ、ＲＯＭ、外部記憶装置の間のデータのやり取りが可能なように接続するバスを有している。また必要に応じて、ハードウェアエンティティに、ＣＤ－ＲＯＭなどの記録媒体を読み書きできる装置（ドライブ）などを設けることとしてもよい。このようなハードウェア資源を備えた物理的実体としては、汎用コンピュータなどがある。

【0059】

ハードウェアエンティティの外部記憶装置には、上述の機能を実現するために必要となるプログラムおよびこのプログラムの処理において必要となるデータなどが記憶されている（外部記憶装置に限らず、例えばプログラムを読み出し専用記憶装置であるＲＯＭに記憶させておくこととしてもよい）。また、これらのプログラムの処理によって得られるデータなどは、ＲＡＭや外部記憶装置などに適宜に記憶される。

【0060】

ハードウェアエンティティでは、外部記憶装置（あるいはＲＯＭなど）に記憶された各プログラムとこの各プログラムの処理に必要なデータが必要に応じてメモリに読み込まれて、適宜にＣＰＵで解釈実行・処理される。その結果、ＣＰＵが所定の機能（上記、…部、…手段などと表した各構成部）を実現する。

【0061】

本発明は上述の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。また、上記実施形態において説明した処理は、記載の順に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されるとしてもよい。

【0062】

既述のように、上記実施形態において説明したハードウェアエンティティ（本発明の装置）における処理機能をコンピュータによって実現する場合、ハードウェアエンティティが有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記ハードウェアエンティティにおける処理機能がコンピュータ上で実現される。

【0063】

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。具体的には、例えば、磁気記録装置として、ハードディスク装置、フレキシブルディスク、磁気テープ等を、光ディスクとして、ＤＶＤ（Digital Versatile Disc）、ＤＶＤ－ＲＡＭ（Random Access Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）、ＣＤ－Ｒ（Recordable）／ＲＷ（ReWritable）等を、光磁気記録媒体として、ＭＯ（Magneto-Optical disc）等を、半導体メモリとしてＥＥＰ－ＲＯＭ（Electronically Erasable and Programmable-Read Only Memory）等を用いることができる。

【0064】

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ－ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。

【0065】

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記憶装置に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。

【0066】

また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、ハードウェアエンティティを構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

【0067】

上述の本発明の実施形態の記載は、例証と記載の目的で提示されたものである。網羅的であるという意思はなく、開示された厳密な形式に発明を限定する意思もない。変形やバリエーションは上述の教示から可能である。実施形態は、本発明の原理の最も良い例証を提供するために、そして、この分野の当業者が、熟考された実際の使用に適するように本発明を色々な実施形態で、また、色々な変形を付加して利用できるようにするために、選ばれて表現されたものである。すべてのそのような変形やバリエーションは、公正に合法的に公平に与えられる幅にしたがって解釈された添付の請求項によって定められた本発明のスコープ内である。

【図1】

【図2】

【図3】

【図4】

【図5】