(57)【特許請求の範囲】
【請求項１】
音響中心点の周囲に予め設定したリスニングエリアに含まれる１以上のリスニングポイントにおいて、１以上の音源を有する音響システムの音響インテンシティと到来音エネルギーとを算出して表示する音響情報表示装置であって、
前記音源ごとの音源信号が入力される音源信号入力手段と、
前記音源信号の音圧に粒子速度を乗じることで、前記リスニングポイントごとに前記音響インテンシティを算出する音響インテンシティ算出手段と、
前記音響インテンシティ算出手段が算出した音響インテンシティを示す音響インテンシティベクトルと、当該音響インテンシティベクトルの起点である他のリスニングポイントから前記リスニングポイントに向いた単位ベクトルとの内積によって、前記リスニングポイントごとに前記到来音エネルギーを算出する到来音エネルギー算出手段と、
前記音響インテンシティ算出手段が算出した音響インテンシティを表示する音響インテンシティ表示手段と、
前記到来音エネルギー算出手段が算出した到来音エネルギーを表示する到来音エネルギー表示手段と、
を備えることを特徴とする音響情報表示装置。
【請求項２】
前記音響インテンシティ算出手段は、予め設定された算出時間間隔で前記音響インテンシティを算出し、
前記到来音エネルギー算出手段は、当該音響インテンシティを示す音響インテンシティベクトルと前記単位ベクトルとの内積によって、前記算出時間間隔で前記到来音エネルギーを算出し、
前記音響インテンシティ表示手段は、当該音響インテンシティ算出手段が算出した音響インテンシティをリアルタイム表示し、
前記到来音エネルギー表示手段は、当該到来音エネルギー算出手段が算出した到来音エネルギーをリアルタイム表示することを特徴とする請求項１に記載の音響情報表示装置。
【請求項３】
所定の第１軸と、当該第１軸に直交する第２軸との２次元座標系において、前記リスニングポイントの位置と、前記音源の位置とがパラメータとして入力されるパラメータ入力手段をさらに備え、
前記音響インテンシティ算出手段は、
前記リスニングポイントと前記音源との位置関係から、前記音源ごとの前記音源信号の第１軸音圧成分と第２軸音圧成分とを算出し、前記音源ごとの前記音源信号の第１軸音圧成分と第２軸音圧成分とをそれぞれ累計して、前記リスニングポイントごとに第１軸音圧成分と第２軸音圧成分とを算出すると共に、
前記リスニングポイントごとの第１軸音圧成分と第２軸音圧成分とに基づいて前記粒子速度の第１軸成分および第２軸成分を算出し、前記音源ごとの前記音源信号の第１軸粒子速度成分と第２軸粒子速度成分とをそれぞれ累計して、前記粒子速度の第１軸成分累計値および第２軸成分累計値を算出し、前記粒子速度の第１軸成分累計値および第２軸成分累計値のそれぞれに前記音源信号の音圧を乗じて、前記音響インテンシティの第１軸成分および第２軸成分を算出することを特徴とする請求項１または請求項２に記載の音響情報表示装置。
【請求項４】
所定の第１軸と、当該第１軸に直交する第２軸と、当該第１軸および当該第２軸に直交する第３軸との３次元座標系において、前記リスニングポイントの位置と、前記音源の位置とがパラメータとして入力されるパラメータ入力手段をさらに備え、
前記音響インテンシティ算出手段は、
前記リスニングポイントと前記音源との位置関係から、前記音源ごとの前記音源信号の第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とを算出し、前記音源ごとの前記音源信号の第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とをそれぞれ累計して、前記リスニングポイントごとに第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とを算出すると共に、
前記リスニングポイントごとの第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とに基づいて前記粒子速度の第１軸成分と第２軸成分と第３軸成分とを算出し、前記音源ごとの前記音源信号の第１軸粒子速度成分と第２軸粒子速度成分と第３軸粒子速度成分とをそれぞれ累計して、前記粒子速度の第１軸成分累計値と第２軸成分累計値と第３軸成分累計値とを算出し、前記粒子速度の第１軸成分累計値と第２軸成分累計値と第３軸成分累計値とのそれぞれに前記音源信号の音圧を乗じて、前記音響インテンシティの第１軸成分と第２軸成分と第３軸成分とを算出することを特徴とする請求項１または請求項２に記載の音響情報表示装置。
【請求項５】
音響中心点の周囲に予め設定したリスニングエリアに含まれる１以上のリスニングポイントにおいて、１以上の音源を有する音響システムの音響インテンシティと到来音エネルギーとを算出して表示するために、コンピュータを、
前記音源ごとの音源信号が入力される音源信号入力手段、
前記音源信号の音圧に粒子速度を乗じることで、前記リスニングポイントごとに前記音響インテンシティを算出する音響インテンシティ算出手段、
前記音響インテンシティ算出手段が算出した音響インテンシティを示す音響インテンシティベクトルと、当該音響インテンシティベクトルの起点である他のリスニングポイントから前記リスニングポイントに向いた単位ベクトルとの内積によって、前記リスニングポイントごとに前記到来音エネルギーを算出する到来音エネルギー算出手段、
前記音響インテンシティ算出手段が算出した音響インテンシティを表示する音響インテンシティ表示手段、
前記到来音エネルギー算出手段が算出した到来音エネルギーを表示する到来音エネルギー表示手段、
として機能させるための音響情報表示プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響中心点の周囲に予め設定したリスニングエリアに含まれる１以上のリスニングポイントにおいて、１以上の音源を有する音響システムの音響インテンシティと到来音エネルギーとを算出して表示する音響情報表示装置およびそのプログラムに関する。
【背景技術】
【０００２】
近年、５．１ｃｈサラウンドシステムなどのマルチチャンネル音響システムが身近になりつつある。このため、マルチチャンネル音響システムにおいて、ラウドネス、瞬時音響インテンシティ、到来音エネルギーなどの音響情報を表示するデジタルスコープが提案されている。
【０００３】
例えば、特許文献１には、マイクロホンを用いた実測により音響情報を可視化する発明が提案されている。この特許文献１に記載の発明は、２次元平面に展開された空間の測定領域において、測定点での周波数の大きさの経時変化を動画として表示する。
【０００４】
また、例えば、特許文献２には、音像定位の位置を確認できるデジタルスコープが提案されている。この特許文献２に記載の発明は、音響中心（以下、「スウィートスポット」）に位置する聴取者を囲むように音源を配置した状態で、音響情報を表示する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平４−２９０９３０号公報
【特許文献２】特開平５−３０４７００号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、特許文献１，２に記載の発明には、以下のような問題点が存在していた。
特許文献１，２に記載の発明は、スウィートスポット１点でしか音響情報を表示できず、このスウィートスポット以外の点では、どのような音像が形成されているかわからない。このため、特許文献１，２に記載の発明は、例えば、複数の作業者が関わる放送番組の制作現場において、音声コントロールルームなどのリスニングエリア内の１点でしか音響情報を確認できず、極めて不便である。
【０００７】
また、特許文献１，２に記載の発明は、音響情報を算出して表示するために、実際の音場における実測データ（例えば、スタジオのスピーカで音を再生したときの音場の実測データ）を用いる必要があり、放送番組の制作現場に採用することが現実でない。そこで、放送番組の制作現場から、リスニングエリア内の様々な位置での音響情報、すなわち、再生音場を客観的に把握したいという強い要望がある。
【０００８】
そこで、本発明は、前記した問題を鑑みて、リスニングエリア内の様々な位置において、再生音場を客観的に把握できる音響情報表示装置およびそのプログラムを提供することを課題とする。
【課題を解決するための手段】
【０００９】
前記した課題を解決するため、本願第１発明に係る音響情報表示装置は、音響中心点の周囲に予め設定したリスニングエリアに含まれる１以上のリスニングポイントにおいて、１以上の音源を有する音響システムの音響インテンシティと到来音エネルギーとを算出して表示する音響情報表示装置であって、音源信号入力手段と、音響インテンシティ算出手段と、到来音エネルギー算出手段と、音響インテンシティ表示手段と、到来音エネルギー表示手段と、を備えることを特徴とする。
【００１０】
かかる構成によれば、音響情報表示装置は、音源信号入力手段によって、音源ごとの音源信号が入力される。また、音響情報表示装置は、音響インテンシティ算出手段によって、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに音響インテンシティを算出する。つまり、音響インテンシティ算出手段は、リスニングエリア内の様々な位置（リスニングポイント）で、音響インテンシティを算出する。
【００１１】
また、音響情報表示装置は、到来音エネルギー算出手段によって、音響インテンシティ算出手段が算出した音響インテンシティを示す音響インテンシティベクトルと、音響インテンシティベクトルの起点である他のリスニングポイントから算出対象リスニングポイントに向いた単位ベクトルとの内積によって、算出対象リスニングポイントごとに到来音エネルギーを算出する。つまり、到来音エネルギー算出手段は、実測データを用いることなく、シミュレーションによって到来音エネルギーを算出する。この到来音エネルギーは、算出対象リスニングポイントに到来する音（到来音）のエネルギーを示している。
【００１２】
また、音響情報表示装置は、音響インテンシティ表示手段によって、音響インテンシティ算出手段が算出した音響インテンシティを表示する。そして、音響情報表示装置は、到来音エネルギー表示手段によって、到来音エネルギー算出手段が算出した到来音エネルギーを表示する。
【００１３】
また、本願第２発明に係る音響情報表示装置は、音響インテンシティ算出手段が、予め設定された算出時間間隔で音響インテンシティを算出し、到来音エネルギー算出手段が、音響インテンシティを示す音響インテンシティベクトルと単位ベクトルとの内積によって、算出時間間隔で到来音エネルギーを算出し、音響インテンシティ表示手段が、音響インテンシティ算出手段が算出した音響インテンシティをリアルタイム表示し、到来音エネルギー表示手段が、到来音エネルギー算出手段が算出した到来音エネルギーをリアルタイム表示することが好ましい。
【００１４】
また、本願第３発明に係る音響情報表示装置は、所定の第１軸と、当該第１軸に直交する第２軸との２次元座標系において、リスニングポイントの位置と、音源の位置とがパラメータとして入力されるパラメータ入力手段をさらに備え、音響インテンシティ算出手段が、リスニングポイントと音源との位置関係から、音源ごとの音源信号の第１軸音圧成分と第２軸音圧成分とを算出し、音源ごとの音源信号の第１軸音圧成分と第２軸音圧成分とをそれぞれ累計して、リスニングポイントごとに第１軸音圧成分と第２軸音圧成分とを算出すると共に、リスニングポイントごとの第１軸音圧成分と第２軸音圧成分とに基づいて粒子速度の第１軸成分および第２軸成分を算出し、音源ごとの音源信号の第１軸粒子速度成分と第２軸粒子速度成分とをそれぞれ累計して、粒子速度の第１軸成分累計値および第２軸成分累計値を算出し、粒子速度の第１軸成分累計値および第２軸成分累計値のそれぞれに音源信号の音圧を乗じて、音響インテンシティの第１軸成分および第２軸成分を算出することが好ましい。
【００１５】
また、本願第４発明に係る音響情報表示装置は、所定の第１軸と、当該第１軸に直交する第２軸と、当該第１軸および当該第２軸に直交する第３軸との３次元座標系において、リスニングポイントの位置と、音源の位置とがパラメータとして入力されるパラメータ入力手段をさらに備え、音響インテンシティ算出手段が、リスニングポイントと音源との位置関係から、音源ごとの音源信号の第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とを算出し、音源ごとの音源信号の第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とをそれぞれ累計して、リスニングポイントごとに第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とを算出すると共に、リスニングポイントごとの第１軸音圧成分と第２軸音圧成分と第３軸音圧成分とに基づいて粒子速度の第１軸成分と第２軸成分と第３軸成分とを算出し、音源ごとの音源信号の第１軸粒子速度成分と第２軸粒子速度成分と第３軸粒子速度成分とをそれぞれ累計して、粒子速度の第１軸成分累計値と第２軸成分累計値と第３軸成分累計値とを算出し、粒子速度の第１軸成分累計値と第２軸成分累計値と第３軸成分累計値とのそれぞれに音源信号の音圧を乗じて、音響インテンシティの第１軸成分と第２軸成分と第３軸成分とを算出することが好ましい。
【００１６】
また、前記した課題を解決するため、本願第５発明に係る音響情報表示プログラムは、音響中心点の周囲に予め設定したリスニングエリアに含まれる１以上のリスニングポイントにおいて、１以上の音源を有する音響システムの音響インテンシティと到来音エネルギーとを算出して表示するために、コンピュータを、音源信号入力手段、音響インテンシティ算出手段、到来音エネルギー算出手段、音響インテンシティ表示手段、到来音エネルギー表示手段、として機能させることを特徴とする。
【００１７】
かかる構成によれば、音響情報表示プログラムは、音源信号入力手段によって、音源ごとの音源信号が入力される。また、音響情報表示プログラムは、音響インテンシティ算出手段によって、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに音響インテンシティを算出する。つまり、音響インテンシティ算出手段は、リスニングエリア内の様々な位置（リスニングポイント）で、音響インテンシティを算出する。
【００１８】
また、音響情報表示プログラムは、到来音エネルギー算出手段によって、音響インテンシティ算出手段が算出した音響インテンシティを示す音響インテンシティベクトルと、音響インテンシティベクトルの起点である他のリスニングポイントから算出対象リスニングポイントに向いた単位ベクトルとの内積によって、リスニングポイントごとに到来音エネルギーを算出する。つまり、到来音エネルギー算出手段は、実測データを用いることなく、シミュレーションによって到来音エネルギーを算出する。この到来音エネルギーは、算出対象リスニングポイントに到来する音（到来音）のエネルギーを示している。
【００１９】
また、音響情報表示プログラムは、音響インテンシティ表示手段によって、音響インテンシティ算出手段が算出した音響インテンシティを表示する。そして、音響情報表示プログラムは、到来音エネルギー表示手段によって、到来音エネルギー算出手段が算出した到来音エネルギーを表示する。
【発明の効果】
【００２０】
本発明によれば、以下のような優れた効果を奏する。
本願第１，５発明によれば、リスニングエリア内の様々な位置（リスニングポイント）において、実測データを用いることなく、シミュレーションによって到来音エネルギーを算出して表示する。このため、本願第１，５発明によれば、リスニングエリア内の様々な位置において、再生音場を客観的に把握することができる。
【００２１】
本願第２発明によれば、音響インテンシティと到来音エネルギーとをリアルタイムで表示するため、リスニングエリアにおける再生音場の経時変化を把握することができる。
【００２２】
本願第３発明によれば、２次元のリスニングエリアにおいて、再生音場を客観的に把握することができる。
本願第４発明によれば、３次元のリスニングエリアにおいて、再生音場を客観的に把握することができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の第１，２実施形態に係る音響システム用デジタルスコープの構成を示すブロック図である。
【図２】本発明の第１実施形態において、リスニングエリアおよびリスニングポイントを説明する図であり、（ａ）は正方形のリスニングエリアの第１例であり、（ｂ）は正方形のリスニングエリアの第２例であり、（ｃ）は長方形のリスニングエリアの例であり、（ｄ）は円形のリスニングエリアの例であり、（ｅ）は三角形のリスニングエリアの例である。
【図３】本発明の第１実施形態において、音響インテンシティの表示画面の例を示す図である。
【図４】本発明の第１実施形態において、到来音エネルギーの表示画面の例を示す画像である。
【図５】本発明の第１実施形態において、瞬時音響インテンシティの算出の第１例を説明する図である。
【図６】本発明の第１実施形態において、瞬時音響インテンシティの算出の第２例を説明する図である。
【図７】本発明の第１実施形態において、到来音エネルギーの算出を説明する図である。
【図８】本発明の第１実施形態に係る音響システム用デジタルスコープの動作を示すフローチャートである。
【図９】本発明の第２実施形態において、リスニングエリアを説明する図であり、（ａ）は立方体のリスニングエリアであり、（ｂ）は球体のリスニングエリアである。
【図１０】本発明の第２実施形態において、瞬時音響インテンシティの算出を説明する図である。
【図１１】本発明の第２実施形態において、瞬時音響インテンシティの算出を説明する図である。
【図１２】本発明の第２実施形態において、瞬時音響インテンシティの算出を説明する図であり、図１１をＺ軸方向から見た図である。
【図１３】本発明の第２実施形態に係る音響システム用デジタルスコープの動作を示すフローチャートである。
【発明を実施するための形態】
【００２４】
以下、本発明の各実施形態について、本発明に係る音響情報表示装置を音響システム用デジタルスコープとして用いた例で説明する。なお、各実施形態において、同一の機能を有する手段には同一の符号を付し、説明を省略した。
【００２５】
（第１実施形態：２次元）
［音響システム用デジタルスコープの構成］
以下、図１を参照し、本発明の第１実施形態に係る音響システム用デジタルスコープ１の構成について説明する。
なお、第１実施形態では、第１軸（Ｘ軸）を実空間の横方向（左右）とし、第２軸（Ｙ軸）を実空間の奥行方向（前後）とした２次元のリスニングエリアの例で説明する。
【００２６】
音響システム用デジタルスコープ（音響情報表示装置）１は、後記するリスニングエリアに含まれる１以上のリスニングポイントにおいて、音響情報として、音響システムの瞬時音響インテンシティと到来音エネルギーとを算出して表示する。このため、音響システム用デジタルスコープ１は、図１に示すように、パラメータ設定手段（パラメータ入力手段）１０と、パラメータ記憶手段２０と、音源信号入力手段３０と、音響インテンシティ算出手段４０と、到来音エネルギー算出手段５０と、演算結果記憶手段６０と、表示手段７０とを備える。
【００２７】
パラメータ設定手段１０は、音響システム用デジタルスコープ１に必要なパラメータが入力されると共に、入力されたパラメータをパラメータ記憶手段２０に記憶（設定）する。
【００２８】
ここで、パラメータ設定手段１０に入力されるパラメータには、例えば、スピーカ位置と、スピーカ数と、音源信号ファイル名と、メッシュグリッド数と、算出時間間隔と、動画フレームレートと、リスニングエリア特定情報と、リスニングポイント位置とが含まれる。これらパラメータは、例えば、作業者によって適切な値が予め決定される。
【００２９】
スピーカ位置（Ｘｓｐｉ，Ｙｓｐｉ）は、音響システムが備える各スピーカの位置を、例えば、２次元座標系（ＸＹ座標系）で示すパラメータである。なお、添え字ｉは、スピーカの番号を示している。
スピーカ数Ｎは、音響システムが備えるスピーカの数を示すパラメータである。例えば、音響システムがモノラルの場合、スピーカ数Ｎ＝１となる。また、例えば、音響システムがステレオの場合、スピーカ数Ｎ＝２となる。また、例えば、音響システムが５．１ｃｈサラウンドの場合、スピーカ数Ｎ＝６となる。また、例えば、音響システムが２２．２ｃｈサラウンドの場合、スピーカ数Ｎ＝２４となる。
図１では、これらスピーカ位置（Ｘｓｐｉ，Ｙｓｐｉ）およびスピーカ数Ｎは、スピーカＰ（スピーカパラメータ）として図示した。
【００３０】
音源信号ファイル名ｆｉｌｅは、音響システムが備えるスピーカごとの音源信号を格納した音源信号ファイルの名称を示すパラメータである。ここで、例えば、音源信号ファイル名ｆｉｌｅでは、音源信号と同じ数、音源信号ファイルの名称がＮ個指定される。
メッシュグリッド数ｍは、リスニングエリアをメッシュ分割する数を示すパラメータである。
算出時間間隔ｔは、瞬時音響インテンシティを算出する時間間隔を示すパラメータである。
図１では、これら音源信号ファイル名ｆｉｌｅ、メッシュグリッド数ｍおよび算出時間間隔ｔは、算出Ｐ（算出パラメータ）として図示した。
【００３１】
動画フレームレートｆは、音響情報をリアルタイム表示（例えば、モニタでの動画表示）する際のフレームレートを示すパラメータである。この動画フレームレートｆは、図１では、表示Ｐ（表示パラメータ）として図示した。
【００３２】
リスニングエリア特定情報は、リスニングエリアの形状および大きさを特定するパラメータである。このリスニングエリア特定情報については、その詳細を後記する。
リスニングポイント位置（Ｐｘ，Ｐｙ）は、それぞれのリスニングポイントの位置を、例えば、２次元座標系（ＸＹ座標系）で示すパラメータである。
図１では、これらリスニングエリア特定情報およびリスニングポイント位置（Ｐｘ，Ｐｙ）は、リスニングエリアＰ（リスニングエリアパラメータ）として図示した。
【００３３】
以下、図２を参照し、リスニングエリアＬＡおよびリスニングポイントＬＰを説明した後、パラメータとしてのリスニングエリア特定情報を説明する（適宜図１参照）。
リスニングエリアＬＡは、瞬時音響インテンシティおよび到来音エネルギーなどの音響情報の表示対象となる領域であり、例えば、音声コントロールルームを表す領域である。また、リスニングエリアＬＡは、音響中心点であるスウィートスポットＳＳと、リスニングポイントＬＰとが含まれる。
リスニングポイントＬＰは、リスニングエリアＬＡ内での音響情報を算出および表示する位置（点）であり、例えば、音声コントロールルームで音響情報を表示したい位置となる。このリスニングポイントＬＰは、リスニングエリアＬＡ内において、任意の数で、任意の位置に設定できる。
【００３４】
例えば、図２（ａ）に示すように、リスニングエリアＬＡは、スウィートスポットＳＳを中心とし、正方形に設定することができる。この場合、リスニングポイントＬＰは、正方形のリスニングエリアＬＡにおいて、左前、右前、左後および右後の４点に設定することができる。また、図２（ｂ）に示すように、リスニングポイントＬＰは、正方形のリスニングエリアＬＡにおいて、マトリクス状に多数設定することもできる。
【００３５】
また、リスニングエリアＬＡは、正方形以外の形状で設定することもできる。例えば、図２（ｃ）に示すように、リスニングエリアＬＡは、長方形に設定することができる。また、例えば、図２（ｄ）に示すように、リスニングエリアＬＡは、スウィートスポットＳＳを中心とし、円形に設定することができる。さらに、例えば、図２（ｅ）に示すように、リスニングエリアＬＡは、スウィートスポットＳＳを中心とし、正三角形に設定することができる。
【００３６】
パラメータとしてのリスニングエリア特定情報について説明する。
前記したように、リスニングエリア特定情報は、リスニングエリアＬＡの形状および大きさを特定する情報である。例えば、図２（ａ）および図２（ｂ）のようにリスニングエリアＬＡが正方形の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの左前点の座標と、リスニングエリアＬＡの一辺の長さとが含まれる。また、図２（ｃ）のようにリスニングエリアＬＡが長方形の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの左前点の座標と、リスニングエリアＬＡの長辺（横方向の辺）の長さと、リスニングエリアＬＡの短辺（縦方向の辺）の長さとが含まれる。また、図２（ｄ）のようにリスニングエリアＬＡが円形の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの中心点の座標と、リスニングエリアＬＡの半径とが含まれる。また、図２（ｅ）のようにリスニングエリアＬＡが正三角形の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの中心点の座標と、リスニングエリアＬＡの一辺の長さとが含まれる。
【００３７】
なお、図２（ｃ）〜図２（ｅ）において、図面を見やすくするためリスニングポイントＬＰの図示を省略したが、リスニングポイントＬＰは、図２（ｃ）〜図２（ｅ）のリスニングエリアＬＡにおいても、任意の数で任意の位置に設定できる。
また、リスニングエリアＬＡの形状は、図２（ａ）〜図２（ｅ）の例に限定されないことは言うまでもない。
【００３８】
なお、スウィートスポットＳＳは、スピーカで再生する音を最も良い条件で聴くことができる点である。例えば、５．１ｃｈサラウンドにおいて、全スピーカを同心円上に配置した場合、スウィートスポットＳＳは、その同心円の中心点となる。このことから、図２（ａ）〜図２（ｅ）に示すように、スウィートスポットＳＳは、リスニングエリアＬＡの中心に位置することが好ましい。その一方、リスニングエリアＬＡの形状またはスピーカの配置によっては、スウィートスポットＳＳが、リスニングエリアＬＡの中心に位置することが適切でない場合もある。この場合、スウィートスポットＳＳは、リスニングエリアＬＡに含まれていればよい。
【００３９】
以下、図１に戻り、音響システム用デジタルスコープ１の構成について説明を続ける。
パラメータ記憶手段２０は、前記したパラメータを記憶するメモリ、ＨＤＤなどの記憶装置である。
【００４０】
音源信号入力手段３０は、スピーカごとの音源信号が入力されると共に、入力された音源信号を音響インテンシティ算出手段４０に出力する。例えば、音源信号入力手段３０は、音源信号ファイル名ｆｉｌｅで指定されたＮ個の音源信号ファイルを開く。そして、音源信号入力手段３０は、Ｎ個の音源信号ファイルを読み込んで、音源信号ファイルのそれぞれに格納された音源信号のデータサイズｄを求める。さらに、音源信号入力手段３０は、求めたデータサイズｄだけ、音源信号ファイルのそれぞれから音源信号を読み込んで、Ｎ個の音源信号を音響インテンシティ算出手段４０に出力する。
【００４１】
音源信号入力手段３０に入力される音源信号としては、例えば、楽音、Ｖｎ単音、風景音、ノイズまたはｓｉｎ波（トーン）がある。
【００４２】
音響インテンシティ算出手段４０は、音源信号入力手段３０から音源信号が入力されると共に、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに瞬時音響インテンシティを算出する。
【００４３】
まず、音響インテンシティ算出手段４０は、瞬時音響インテンシティを算出するための前処理を行う。例えば、音響インテンシティ算出手段４０は、パラメータ記憶手段２０から、瞬時音響インテンシティの算出に必要なパラメータを読み出す。また、音響インテンシティ算出手段４０は、Ｘ軸方向に２ｍ×Ｙ軸方向に２ｍのリスニングエリアを、メッシュグリッド数ｍ＝１００個のメッシュに分割して、このリスニングエリアに１０１×１０１点をプロットする。そして、音響インテンシティ算出手段４０は、リスニングポイント（プロットされた各点）までの時間遅れを算出する。さらに、音響インテンシティ算出手段４０は、算出した時間遅れを考慮して、音源信号の音圧（算出時間間隔の中心時刻における音源信号の音圧）を用いて、リスニングポイントごとに音圧分布を算出する。
ここでは、リスニングエリア内に多数のリスニングポイント（音場の各点）が、プロットされた各点と一致するように設定されている（図２（ｂ）参照）。
【００４４】
そして、前処理を終了したら、音響インテンシティ算出手段４０は、算出した音圧分布と粒子速度とに基づいて、リスニングポイントごとに、算出時間間隔ｔ（例えば、０．０５秒単位）で瞬時音響インテンシティを算出する。本実施形態では、音響インテンシティ算出手段４０は、瞬時音響インテンシティのＸ軸成分（第１軸成分）およびＹ軸成分（第２軸成分）をそれぞれ算出する。その後、音響インテンシティ算出手段４０は、算出した瞬時音響インテンシティを到来音エネルギー算出手段５０に出力すると共に、演算結果記憶手段６０に記憶する。
なお、瞬時音響インテンシティの算出については、その詳細を後記する。
【００４５】
到来音エネルギー算出手段５０は、音響インテンシティ算出手段４０から瞬時音響インテンシティが入力されると共に、入力された瞬時音響インテンシティを示す瞬時音響インテンシティベクトルと、瞬時音響インテンシティベクトルの起点である他のリスニングポイントから算出対象リスニングポイントに向いた単位ベクトルとの内積によって、算出対象リスニングポイントごとに到来音エネルギーを算出する。ここで、到来音エネルギー算出手段５０は、算出時間間隔ｔで到来音エネルギーを算出する。そして、到来音エネルギー算出手段５０は、算出した到来音エネルギーを演算結果記憶手段６０に記憶する。
なお、到来音エネルギーの算出については、その詳細を後記する。
【００４６】
演算結果記憶手段６０は、演算結果として、音響インテンシティ算出手段４０が算出した瞬時音響インテンシティ（図１では「音響インテンシティ」）と、到来音エネルギー算出手段５０が算出した到来音エネルギーとを記憶するメモリ、ＨＤＤなどの記憶装置である。
【００４７】
表示手段７０は、演算結果記憶手段６０が記憶する演算結果を、図示を省略したモニタにグラフィック表示するものであり、音響インテンシティ表示手段７１と、到来音エネルギー表示手段７３とを備える。
【００４８】
音響インテンシティ表示手段７１は、演算結果記憶手段６０に記憶された瞬時音響インテンシティをモニタにグラフィック表示する。ここで、音響インテンシティ表示手段７１は、例えば、動画フレームレートｆでモニタの表示画面を更新して、瞬時音響インテンシティをリアルタイムでグラフィック表示する。
【００４９】
以下、図３を参照し、瞬時音響インテンシティの表示画面について説明する（適宜図１参照）。
図３の例では、リスニングエリアＬＡは、正方形であると共に、その中心がスウィートスポットＳＳであるとする。また、リスニングポイントＬＰは、スウィートスポットＳＳを中心として、Ｘ軸方向に５個およびＹ軸方向に５個、合計２５個が一定間隔で並んでいるとする。つまり、スウィートスポットＳＳは、リスニングエリアＬＡの中心に位置するリスニングポイントＬＰに重なっている。また、図３では、Ｘ軸およびＹ軸の数値は、リスニングエリアＬＡの中心位置からの距離を示している。
【００５０】
図３に示すように、瞬時音響インテンシティの表示画面には、２５個のリスニングポイントＬＰについて、それぞれの瞬時音響インテンシティが表示されている。この瞬時音響インテンシティの表示画面では、瞬時音響インテンシティが瞬時音響インテンシティベクトルＩｔとして表示されている。つまり、瞬時音響インテンシティベクトルＩｔは、Ｘ軸方向の長さＩｔｘが、音響インテンシティ算出手段４０によって算出された瞬時音響インテンシティのＸ軸成分を示している。また、瞬時音響インテンシティベクトルＩｔは、Ｙ軸方向の長さＩｔｙが、音響インテンシティ算出手段４０によって算出された瞬時音響インテンシティのＹ軸成分を示している。
【００５１】
このように、瞬時音響インテンシティの表示画面によって、作業者が、各リスニングポイントＬＰにおける瞬時音響インテンシティを容易に確認できる。例えば、作業者は、それぞれの瞬時音響インテンシティベクトルＩｔの向きおよび大きさが近似するほど、リスニングエリアＬＡにおける再生音場が良好であると判断できる。
【００５２】
以下、図１に戻り、音響システム用デジタルスコープ１の構成について説明を続ける。
到来音エネルギー表示手段７３は、演算結果記憶手段６０に記憶された到来音エネルギーをモニタにグラフィック表示する。ここで、到来音エネルギー表示手段７３は、例えば、動画フレームレートｆでモニタの表示画面を更新して、到来音エネルギーをリアルタイムでグラフィック表示する。
【００５３】
以下、図４，図５を参照し、到来音エネルギーの表示画面について説明する（適宜図１参照）。
ここでは、図５に示すように、音響システムがステレオの例で説明する。リスニングエリアＬＡは、Ｘ軸方向に２ｍ×Ｙ軸方向に２ｍの正方形とする。また、スウィートスポットＳＳは、リスニングエリアＬＡの中心に位置する。また、スピーカＳＰ１は、スウィートスポットＳＳから左前方３ｍに位置し、スピーカＳＰ２は、スウィートスポットＳＳから右前方３ｍに位置する。ここで、スピーカＳＰ１，ＳＰ２の距離も３ｍとする。つまり、スウィートスポットＳＳを中心とした半径３ｍの円Ｃｒを描くと、スピーカＳＰ１，ＳＰ２は、この円Ｃｒの円周上に位置する。
【００５４】
ここで、リスニングエリアＬＡ内にリスニングポイントＬＰを多数設定する（図２（ｂ）参照）。これらリスニングポイントＬＰのうち、左前、右前、左後および右後に位置する４個を代表的なリスニングポイントとする。このとき、図４には、代表的なリスニングポイントでの到来音エネルギーの表示画面を示した。つまり、図４上段左側の画像が、左前のリスニングポイントでの到来音エネルギーの表示画面であり、図４上段右側の画像が、右前のリスニングポイントでの到来音エネルギーの表示画面である。また、図４下段左側の画像が、左後のリスニングポイントでの到来音エネルギーの表示画面であり、図４下段右側の画像が、右後のリスニングポイントでの到来音エネルギーの表示画面である。
【００５５】
図４に示すように、到来音エネルギーの表示画面では、代表的なリスニングポイントでの到来音エネルギーがカラートーンで表示されている。このカラートーンは、例えば、到来音エネルギーが大きい方から小さい方までを順に、赤色、黄色、緑色、水色、紺色で表している。また、到来音エネルギーの表示画面において、「ｘ」，「ｙ」は、リスニングポイントのＸ軸方向およびＹ軸方向の位置を示している。つまり、図４の到来音エネルギーの表示画面は、音場各点（図２（ｂ）のように多数設定したリスニングエリアの各点）から１点のリスニングポイントに到来する音のエネルギーについて、その１点のリスニングポイント方向の成分分布を示している。このように、到来音エネルギーの表示画面では、作業者が、色の分布により、各リスニングポイントＬＰでの到来音エネルギーを容易に確認できる。
【００５６】
以下、図５，図６を参照し、瞬時音響インテンシティの算出について、第１例および第２例をあげて詳細に説明する（適宜図１参照）。
【００５７】
＜瞬時音響インテンシティの算出：第１例＞
まず、図５を参照して、瞬時音響インテンシティの算出について、第１例を説明する。ここでは、説明を簡易にするため、１個のリスニングポイントＬＰについて説明する。
【００５８】
ある１点での瞬時音響インテンシティＩは、下記の式（１）に示すように、音圧ｐと粒子速度ｕとの積によって求めることができる。
【００５９】
【数１】
0005480748.tif 000002
【００６０】
音響インテンシティ算出手段４０は、下記の手順でリスニングポイントＬＰの瞬時音響インテンシティを算出する。まず、音響インテンシティ算出手段４０は、パラメータ記憶手段２０から、瞬時音響インテンシティの算出に必要なパラメータ（例えば、スピーカ位置およびリスニングポイント位置）を読み出す。そして、音響インテンシティ算出手段４０は、スピーカＳＰ１について、リスニングポイントＬＰとスピーカＳＰ１との位置関係から、音源信号のＸ軸音圧成分（第１軸音圧成分）とＹ軸音圧成分（第２軸音圧成分）とを算出する。
【００６１】
具体的には、音響インテンシティ算出手段４０は、リスニングポイントＬＰを通過するＹ軸方向の軸線と、スピーカＳＰ１の位置（Ｘｓｐ１，Ｙｓｐ１）からリスニングポイントＬＰの位置（Ｐｘ，Ｐｙ）までの線分とのなす角θ１を算出する。そして、音響インテンシティ算出手段４０は、なす角θ１および三角関数を用いて、スピーカＳＰ１の音響信号の音圧からＸ軸音圧成分とＹ軸音圧成分とを算出する。つまり、音響インテンシティ算出手段４０は、ｓｉｎθ１の値と、スピーカＳＰ１の音響信号の音圧絶対値｜ｐ１（ｔ）｜とを乗じて、Ｘ軸音圧成分ｓｉｎθ１＊｜ｐ１（ｔ）｜を算出する。また、音響インテンシティ算出手段４０は、ｃｏｓθ１の値と、スピーカＳＰ１の音響信号の音圧絶対値｜ｐ１（ｔ）｜とを乗じて、Ｙ軸音圧成分ｃｏｓθ１＊｜ｐ１（ｔ）｜を算出する。
【００６２】
また、音響インテンシティ算出手段４０は、スピーカＳＰ１と同様、スピーカＳＰ２についても、リスニングポイントＬＰの位置（Ｐｘ，Ｐｙ）とスピーカＳＰ２の位置（Ｘｓｐ２，Ｙｓｐ２）との関係から、音源信号のＸ軸音圧成分とＹ軸音圧成分とを算出する。つまり、音響インテンシティ算出手段４０は、スピーカＳＰ２について、Ｘ軸音圧成分ｓｉｎθ２＊｜ｐ２（ｔ）｜と、Ｙ軸音圧成分ｃｏｓθ２＊｜ｐ２（ｔ）｜を算出する。
【００６３】
ここで、リスニングポイントＬＰとスピーカＳＰ１，ＳＰ２との距離が離れるほど、音には遅延が生じる。そこで、音響インテンシティ算出手段４０は、この遅延を瞬時音響インテンシティに反映させるべく、リスニングポイントＬＰとスピーカとの距離に応じて、スピーカＳＰ１，ＳＰ２のｐ１（ｔ），ｐ２（ｔ）を遅延させる。
【００６４】
次に、音響インテンシティ算出手段４０は、スピーカＳＰ１，ＳＰ２のＸ軸音圧成分とＹ軸音圧成分を累計して、リスニングポイントＬＰのＸ軸音圧成分とＹ軸音圧成分とを算出する。ここで、音響システム用デジタルスコープ１では、リスニングポイントＬＰへの到来音を扱う。このため、音響インテンシティ算出手段４０は、Ｘ軸方向では、左方向（図５左側）からリスニングポイントＬＰに到来する音を正（＋）の符号とし、右方向（図５右側）からリスニングポイントＬＰに到来する音を負（−）の符号として累計する。また、音響インテンシティ算出手段４０は、Ｙ軸方向では、後方向（図５下側）からリスニングポイントＬＰに到来する音を正（＋）の符号とし、前方向（図５上側）からリスニングポイントＬＰに到来する音を負（−）の符号として累計する。
【００６５】
この場合、図５に示すように、音響インテンシティ算出手段４０は、リスニングポイントＬＰの左側に位置するスピーカＳＰ１について、そのＸ軸音圧成分の符号を正（＋）とする。また、音響インテンシティ算出手段４０は、リスニングポイントＬＰの右側に位置するスピーカＳＰ２について、そのＸ軸音圧成分の符号を負（−）とする。そして、音響インテンシティ算出手段４０は、下記の式（２）で示すように、スピーカＳＰ１，Ｓ２のＸ軸音圧成分を累計して、リスニングポイントＬＰのＸ軸音圧成分を算出する。なお、式（２）では、正（＋）の符号を省略した。
【００６６】
【数２】
0005480748.tif 000003
【００６７】
また、音響インテンシティ算出手段４０は、スピーカＳＰ１，ＳＰ２の両方がリスニングポイントＬＰの前方向（図５上側）に位置するので、Ｙ軸音圧成分の符号を負（−）とする。そして、音響インテンシティ算出手段４０は、下記の式（３）で示すように、スピーカＳＰ１，Ｓ２のＹ軸音圧成分を累計して、リスニングポイントＬＰのＹ軸音圧成分を算出する。
【００６８】
【数３】
0005480748.tif 000004
【００６９】
次に、音響インテンシティ算出手段４０は、リスニングポイントＬＰのＸ軸音圧成分とＹ軸音圧成分とに基づいて粒子速度のＸ軸成分およびＹ軸成分を算出し、スピーカＳＰ１，ＳＰ２の音源信号のＸ軸粒子速度成分（第１軸粒子速度成分）とＹ軸粒子速度成分（第２軸粒子速度成分）とをそれぞれ累計して、粒子速度のＸ軸成分累計値ｕｘ（ｔ）およびＹ軸成分累計値ｕｙ（ｔ）を算出する。具体的には、音響インテンシティ算出手段４０は、下記の式（４）に示すように、空気の固有抵抗ρｃの逆数にリスニングポイントＬＰのＸ軸音圧成分を乗じて粒子速度のＸ軸成分累計値ｕｘ（ｔ）を算出する。また、音響インテンシティ算出手段４０は、空気の固有抵抗ρｃの逆数にリスニングポイントＬＰのＹ軸音圧成分を乗じて粒子速度のＹ軸成分累計値ｕｙ（ｔ）を算出する。
【００７０】
【数４】
0005480748.tif 000005
【００７１】
その後、音響インテンシティ算出手段４０は、粒子速度のＸ軸成分累計値ｕｘ（ｔ）およびＹ軸成分累計値ｕｙ（ｔ）のそれぞれに音源信号の音圧を乗じて、瞬時音響インテンシティのＸ軸成分およびＹ軸成分を算出する。具体的には、音響インテンシティ算出手段４０は、下記の式（５）に示すように、粒子速度のＸ軸成分累計値ｕｘ（ｔ）と音源信号の音圧ｐ（ｔ）とを乗じて瞬時音響インテンシティのＸ軸成分Ｉｘ（ｔ）を算出する。また、音響インテンシティ算出手段４０は、粒子速度のＹ軸成分累計値ｕｙ（ｔ）と音源信号の音圧ｐ（ｔ）とを乗じて瞬時音響インテンシティのＹ軸成分Ｉｙ（ｔ）を算出する。なお、式（５）において、音圧ｐ（ｔ）は、音響インテンシティ算出手段４０に入力されたＮ個の音源信号の音圧合計値である。つまり、瞬時音響インテンシティは、あるリスニングポイントＬＰにおいて、粒子速度ベクトルの累積値と、そのリスニングポイントＬＰに対して各チャンネルから到来する音の音圧合計値とを乗じたものとなる。
【００７２】
【数５】
0005480748.tif 000006
【００７３】
＜瞬時音響インテンシティの算出：第２例＞
この第２例では、音響システムが５．１ｃｈサラウンドの例で説明する。具体的には、図６に示すように、スピーカＳＰ１は、スウィートスポットＳＳから左前方３ｍに位置し、スピーカＳＰ２は、スウィートスポットＳＳから右前方３ｍに位置する。また、スピーカＳＰ３は、スウィートスポットＳＳから右後方３ｍに位置し、スピーカＳＰ４は、スウィートスポットＳＳから真後３ｍに位置する。そして、スピーカＳＰ５は、スウィートスポットＳＳから左後方３ｍに位置する。つまり、スウィートスポットＳＳを中心とした半径３ｍの円Ｃｒを描くと、スピーカＳＰ１〜ＳＰ５は、この円Ｃｒの円周上に位置する。ここでは、説明を簡易にするため、重低音担当（０．１ｃｈ）のスピーカを省略した。
【００７４】
また、図６に示すように、リスニングエリアＬＡは、Ｘ軸方向に２ｍ×Ｙ軸方向に２ｍの正方形とする。また、リスニングポイントＬＰは、リスニングエリアＬＡの左前側に位置する。ここでは、説明を簡易にするため、１個のリスニングポイントＬＰについて説明する。
【００７５】
音響インテンシティ算出手段４０は、スピーカＳＰ１，ＳＰ２と同様、スピーカＳＰ３についても、リスニングポイントＬＰとスピーカＳＰ３との位置関係から、音源信号のＸ軸音圧成分とＹ軸音圧成分とを算出する。具体的には、音響インテンシティ算出手段４０は、スピーカＳＰ３について、Ｘ軸音圧成分ｓｉｎθ３＊｜ｐ３（ｔ）｜と、Ｙ軸音圧成分ｃｏｓθ３＊｜ｐ３（ｔ）｜を算出する。そして、音響インテンシティ算出手段４０は、スピーカＳＰ４について、Ｘ軸音圧成分ｓｉｎθ４＊｜ｐ４（ｔ）｜と、Ｙ軸音圧成分ｃｏｓθ４＊｜ｐ４（ｔ）｜を算出する。さらに、音響インテンシティ算出手段４０は、スピーカＳＰ５について、Ｘ軸音圧成分ｓｉｎθ５＊｜ｐ５（ｔ）｜と、Ｙ軸音圧成分ｃｏｓθ５＊｜ｐ５（ｔ）｜を算出する。
【００７６】
次に、音響インテンシティ算出手段４０は、スピーカＳＰ１〜Ｓ５のＸ軸音圧成分とＹ軸音圧成分を累計して、リスニングポイントＬＰのＸ軸音圧成分とＹ軸音圧成分とを算出する。このとき、音響インテンシティ算出手段４０は、スピーカＳＰ１，ＳＰ２と同様、スピーカＳＰ３〜ＳＰ５のＸ軸音圧成分およびＹ軸音圧成分について、その符号を求める。
【００７７】
ここで、音響インテンシティ算出手段４０は、リスニングポイントＬＰの右側に位置するスピーカＳＰ３，ＳＰ４について、そのＸ軸音圧成分の符号を負（−）とする。また、音響インテンシティ算出手段４０は、リスニングポイントＬＰの左側に位置するスピーカＳＰ５について、そのＸ軸音圧成分の符号を正（＋）とする。その後、音響インテンシティ算出手段４０は、下記の式（６）で示すように、スピーカＳＰ１〜Ｓ５のＸ軸音圧成分を累計して、リスニングポイントＬＰのＸ軸音圧成分を算出する
【００７８】
【数６】
0005480748.tif 000007
【００７９】
また、音響インテンシティ算出手段４０は、スピーカＳＰ３〜ＳＰ５がリスニングポイントＬＰの後方向（図５下側）に位置するので、Ｙ軸音圧成分の符号を正（＋）とする。その後、音響インテンシティ算出手段４０は、下記の式（７）で示すように、スピーカＳＰ１〜Ｓ５のＹ軸音圧成分を累計して、リスニングポイントＬＰのＹ軸音圧成分を算出する。
【００８０】
【数７】
0005480748.tif 000008
【００８１】
次に、音響インテンシティ算出手段４０は、第１例と同様、リスニングポイントＬＰのＸ軸音圧成分とＹ軸音圧成分とに基づいて粒子速度のＸ軸成分およびＹ軸成分を算出し、スピーカＳＰ１〜ＳＰ５の音源信号のＸ軸粒子速度成分とＹ軸粒子速度成分とをそれぞれ累計して、粒子速度のＸ軸成分累計値ｕｘ（ｔ）およびＹ軸成分累計値ｕｙ（ｔ）を算出する。具体的には、音響インテンシティ算出手段４０は、下記の式（８）に示すように、空気の固有抵抗ρｃの逆数にリスニングポイントＬＰのＸ軸音圧成分を乗じて粒子速度のＸ軸成分累計値ｕｘ（ｔ）を算出する。
【００８２】
【数８】
0005480748.tif 000009
【００８３】
また、音響インテンシティ算出手段４０は、下記の式（９）に示すように、空気の固有抵抗ρｃの逆数にリスニングポイントＬＰのＹ軸音圧成分を乗じてＹ軸成分累計値ｕｙ（ｔ）を算出する。
【００８４】
【数９】
0005480748.tif 000010
【００８５】
その後、音響インテンシティ算出手段４０は、前記した式（５）を用いて、粒子速度のＸ軸成分累計値ｕｘ（ｔ）およびＹ軸成分累計値ｕｙ（ｔ）のそれぞれに音源信号の音圧ｐ（ｔ）を乗じて、瞬時音響インテンシティのＸ軸成分Ｉｘ（ｔ）およびＹ軸成分Ｉｙ（ｔ）を算出する。
【００８６】
なお、第１例および第２例では、リスニングポイントＬＰが１個の例で説明したが、音響インテンシティ算出手段４０は、リスニングエリアＬＡ内の全てのリスニングポイントＬＰで瞬時音響インテンシティを算出することは言うまでもない。
また、スピーカの位置および数が、第１例および第２例に限定さないことは言うまでもない。
【００８７】
＜到来音エネルギーの算出＞
以下、図７を参照し、到来音エネルギーの算出について詳細に説明する（適宜図１，図３参照）。
図７では、算出対象リスニングポイントＬＰ１は、到来音エネルギーの算出対象となっているリスニングポイントである。また、他のリスニングポイントＬＰ２は、算出対象リスニングポイントＬＰ１以外のリスニングポイント、つまり、到来音エネルギーの算出対象外のリスニングポイントである。
【００８８】
例えば、図３のように、２５個のリスニングポイントについて、瞬時音響インテンシティを算出したとする。この場合、到来音エネルギー算出手段５０は、２５個のリスニングポイントの中で何れか１個を、算出対象リスニングポイントＬＰ１として選択する。つまり、他のリスニングポイントＬＰ２は、２５個のリスニングポイントのうち、算出対象リスニングポイントＬＰ１として選択されなかった２４個のリスニングポイントである。そして、到来音エネルギー算出手段５０は、選択した算出対象リスニングポイントＬＰ１について、到来音エネルギーを算出する。その後、到来音エネルギー算出手段５０は、２５個のリスニングポイントの中で、まだ到来音エネルギーを算出していないリスニングポイントを、新たな算出対象リスニングポイントＬＰ１として選択する。そして、到来音エネルギー算出手段５０は、新たな算出対象リスニングポイントＬＰ１について、到来音エネルギーを算出する。このように、到来音エネルギー算出手段５０は、リスニングエリア内の全リスニングポイントで到来音エネルギーを算出するまで、算出対象リスニングポイントＬＰ１の選択と、到来音エネルギーの算出とを繰り返す。
【００８９】
説明を簡易にするため、図７に示すように、他のリスニングポイントＬＰ２が１個の場合を考える。ここで、他のリスニングポイントＬＰ２について、瞬時音響インテンシティは、他のリスニングポイントＬＰ２を起点とした瞬時音響インテンシティベクトルＩｔとして表すことができる。また、他のリスニングポイントＬＰ２について、単位ベクトルｎは、瞬時音響インテンシティベクトルＩｔの起点である他のリスニングポイントＬＰ２から、算出対象リスニングポイントＬＰ１に向いた単位長さの法線ベクトルとして表すことができる。さらに、瞬時音響インテンシティベクトルＩｔと単位ベクトルｎとのなす角をθとして表すことができる。
【００９０】
到来音エネルギー算出手段５０は、瞬時音響インテンシティベクトルＩｔと単位ベクトルｎとの内積によって、算出対象リスニングポイントＬＰ１での到来音エネルギーを算出する。具体的には、到来音エネルギー算出手段５０は、下記の式（１０）に示すように、瞬時音響インテンシティベクトルＩｔの大きさと、単位ベクトルｎの大きさと、ｃｏｓθの値との積によって、算出対象リスニングポイントＬＰ１での到来音エネルギーを算出する。なお、式（１０）では、Ｉ（ｔ）が瞬時音響インテンシティベクトルＩｔである。
【００９１】
【数１０】
0005480748.tif 000011
【００９２】
なお、図７では、他のリスニングポイントＬＰ２が１個として説明したが、到来音エネルギー算出手段５０は、算出対象リスニングポイントＬＰ１以外の全リスニングポイントから、算出対象リスニングポイントＬＰ１への到来音エネルギーを算出することは言うまでもない。
【００９３】
［音響システム用デジタルスコープの動作］
以下、図８を参照して、音響システム用デジタルスコープ１の動作について説明する（適宜図１参照）。
【００９４】
まず、音響システム用デジタルスコープ１は、パラメータ設定手段１０によって、音響システム用デジタルスコープ１に必要なパラメータが入力されると共に、入力されたパラメータをパラメータ記憶手段２０に記憶する（ステップＳ１）。
【００９５】
また、音響システム用デジタルスコープ１は、音源信号入力手段３０によって、Ｎ個の音源信号ファイルを読み込んで、音源信号ファイルのそれぞれに格納された音源信号のデータサイズｄを求める（ステップＳ２）。そして、音響システム用デジタルスコープ１は、音源信号入力手段３０によって、求めたデータサイズｄだけ、音源信号ファイルのそれぞれから、Ｎ個の音源信号を読み込む（ステップＳ３）。
【００９６】
また、音響システム用デジタルスコープ１は、音響インテンシティ算出手段４０によって、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに瞬時音響インテンシティを算出する。このとき、音響システム用デジタルスコープ１は、音響インテンシティ算出手段４０によって、リスニングポイントとスピーカとの位置関係から、スピーカごとに音源信号のＸ軸音圧成分とＹ軸音圧成分とを算出する。そして、音響システム用デジタルスコープ１は、音響インテンシティ算出手段４０によって、スピーカごとのＸ軸音圧成分とＹ軸音圧成分とを累計して、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分とを算出する。さらに、音響システム用デジタルスコープ１は、音響インテンシティ算出手段４０によって、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分とに基づいて粒子速度のＸ軸成分およびＹ軸成分を算出し、スピーカごとの音源信号のＸ軸粒子速度成分とＹ軸粒子速度成分とをそれぞれ累計して、粒子速度のＸ軸成分累計値およびＹ軸成分累計値を算出する。その後、音響システム用デジタルスコープ１は、音響インテンシティ算出手段４０によって、粒子速度のＸ軸成分累計値およびＹ軸成分累計値のそれぞれに音源信号の音圧を乗じて、瞬時音響インテンシティのＸ軸成分およびＹ軸成分を算出する（ステップＳ４）。
【００９７】
また、音響システム用デジタルスコープ１は、到来音エネルギー算出手段５０によって、瞬時音響インテンシティベクトルと、瞬時音響インテンシティベクトルの起点である他のリスニングポイントからリスニングポイントに向いた単位ベクトルとの内積によって、リスニングポイントごとに到来音エネルギーを算出する。つまり、音響システム用デジタルスコープ１は、到来音エネルギー算出手段５０によって、前記した式（１０）を用いて、リスニングポイントごとに到来音エネルギーを算出する（ステップＳ５）。
【００９８】
また、音響システム用デジタルスコープ１は、音響インテンシティ表示手段７１によって、瞬時音響インテンシティをモニタにグラフィック表示する。このとき、音響システム用デジタルスコープ１は、音響インテンシティ表示手段７１によって、瞬時音響インテンシティをリアルタイムでグラフィック表示する（ステップＳ６）。
【００９９】
また、音響システム用デジタルスコープ１は、到来音エネルギー表示手段７３によって、到来音エネルギーをモニタにグラフィック表示する。このとき、音響システム用デジタルスコープ１は、到来音エネルギー表示手段７３によって、到来音エネルギーをリアルタイムでグラフィック表示する（ステップＳ７）。
【０１００】
以上のように、本発明の第１実施形態に係る音響システム用デジタルスコープ１は、２次元のリスニングエリア内の様々な位置（リスニングポイント）において、到来音エネルギーを算出して表示する。これによって、音響システム用デジタルスコープ１は、作業者が、リスニングエリア内の様々な位置において、再生音場を客観的に把握することができる。
【０１０１】
また、音響システム用デジタルスコープ１は、実際の音場における実測データを用いて到来音エネルギーなどの音響情報を算出する従来技術に比べて、実測データを用いることなく、シミュレーションによって到来音エネルギーを算出するため、放送番組の制作現場で利用し易くなる。
【０１０２】
さらに、音響システム用デジタルスコープ１は、音響インテンシティと到来音エネルギーとをリアルタイムで算出してグラフィック表示するため、リスニングエリアにおける再生音場の経時変化を把握することができる。例えば、音源信号として楽曲を入力した場合、音響システム用デジタルスコープ１は、作業者が、楽曲の先頭部分、中間部分および終了部分などで、リスニングエリアにおける再生音場がどのように変化するか把握できる。
【０１０３】
従って、音響システム用デジタルスコープ１は、例えば、放送番組、ラジオ番組および音楽コンテンツの効率的な制作支援を可能とする。この他、音響システム用デジタルスコープ１は、例えば、音響空間や音響システムの客観的な評価にも利用することができる。
【０１０４】
（第２実施形態：３次元）
［音響システム用デジタルスコープの構成］
以下、図１に戻り、本発明の第２実施形態に係る音響システム用デジタルスコープ１Ｂの構成について、第１実施形態と異なる点を説明する。
なお、第２実施形態では、第１軸（Ｘ軸）を実空間の横方向（左右）とし、第２軸（Ｙ軸）を実空間の奥行方向（前後）とし、第３軸（Ｚ軸）を実空間の高さ方向（上下）とした３次元のリスニングエリアの例で説明する。
【０１０５】
音響システム用デジタルスコープ１Ｂは、図１に示すように、パラメータ設定手段（パラメータ入力手段）１０Ｂと、パラメータ記憶手段２０と、音源信号入力手段３０と、音響インテンシティ算出手段４０Ｂと、到来音エネルギー算出手段５０Ｂと、演算結果記憶手段６０と、表示手段７０Ｂとを備える。
【０１０６】
パラメータ設定手段１０Ｂは、音響システム用デジタルスコープ１Ｂに必要なパラメータが入力されると共に、入力されたパラメータをパラメータ記憶手段２０に記憶（設定）する。
【０１０７】
パラメータ設定手段１０Ｂに入力されるパラメータを説明する。
スピーカ位置（Ｘｓｐｉ，Ｙｓｐｉ，Ｚｓｐｉ）は、音響システムが備える各スピーカの位置を、例えば、３次元座標系（ＸＹＺ座標系）で示すパラメータである。
リスニングポイント位置（Ｐｘ，Ｐｙ，Ｐｚ）は、それぞれのリスニングポイントの位置を、例えば、３次元座標系（ＸＹＺ座標系）で示すパラメータである。
【０１０８】
以下、図９を参照し、リスニングエリアＬＡの一例を説明した後、パラメータとしてのリスニングエリア特定情報を説明する（適宜図１参照）。
【０１０９】
図９（ａ）に示すように、リスニングエリアＬＡは、スウィートスポットＳＳを中心とし、立方体に設定することができる。この場合、図示を省略したリスニングポイントは、例えば、上側左前、上側右前、上側左後および上側右後の４点、および、下側左前、下側右前、下側左後および下側右後の４点（合計８点）、に設定することができる。
【０１１０】
また、図９（ｂ）に示すように、リスニングエリアＬＡは、スウィートスポットＳＳを中心とし、Ｘ軸、Ｙ軸およびＺ軸の３方向で同一半径の球体に設定することができる。
なお、リスニングエリアは、これ以外、直方体、楕円体などの任意の形状に設定することができる（不図示）。また、リスニングポイントは、任意の場所に任意の数設定できることは言うまでもない。
【０１１１】
ここで、パラメータとしてのリスニングエリア特定情報について説明する。
リスニングエリア特定情報は、第１実施形態と同様、リスニングエリアＬＡの形状および大きさを特定する情報である。例えば、図９（ａ）のようにリスニングエリアＬＡが立方体の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの上側左前点の座標と、リスニングエリアＬＡの一辺の長さとが含まれる。また、図９（ｂ）のようにリスニングエリアＬＡが球体の場合、リスニングエリア特定情報は、例えば、リスニングエリアＬＡの中心点の座標と、リスニングエリアＬＡの半径とが含まれる。
なお、スピーカ数と、音源信号ファイル名と、メッシュグリッド数と、算出時間間隔と、動画フレームレートとは、第１実施形態と同様のため、説明を省略する。
【０１１２】
以下、図１に戻り、音響システム用デジタルスコープ１Ｂの構成について説明を続ける。
音響インテンシティ算出手段４０Ｂは、音源信号入力手段３０から音源信号が入力されると共に、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに瞬時音響インテンシティを算出する。ここで、音響インテンシティ算出手段４０Ｂは、瞬時音響インテンシティのＸ軸成分、Ｙ軸成分およびＺ軸成分を算出する。そして、音響インテンシティ算出手段４０Ｂは、算出した瞬時音響インテンシティのＸ軸成分、Ｙ軸成分およびＺ軸成分を到来音エネルギー算出手段５０Ｂに出力すると共に、演算結果記憶手段６０に記憶する。
なお、瞬時音響インテンシティの算出については、詳細を後記する。
【０１１３】
到来音エネルギー算出手段５０Ｂは、音響インテンシティ算出手段４０Ｂから瞬時音響インテンシティが入力されると共に、入力された瞬時音響インテンシティを示す瞬時音響インテンシティベクトルと、単位ベクトルとの内積によって、リスニングポイントごとに到来音エネルギーを算出する。ここで、到来音エネルギー算出手段５０Ｂは、到来音エネルギー算出手段５０と同様、前記した式（１０）を用いて、到来音エネルギーを算出する。そして、到来音エネルギー算出手段５０は、算出した到来音エネルギーを演算結果記憶手段６０に記憶する。
【０１１４】
なお、到来音エネルギー算出手段５０Ｂは、瞬時音響インテンシティベクトルが２次元のベクトルであるか、３次元のベクトルあるか以外、到来音エネルギー算出手段５０と同様のため、詳細な説明を省略する。
【０１１５】
表示手段７０Ｂは、演算結果記憶手段６０の演算結果を参照し、図示を省略したモニタに音響情報をグラフィック表示するものであり、音響インテンシティ表示手段７１Ｂと、到来音エネルギー表示手段７３Ｂとを備える。
【０１１６】
音響インテンシティ表示手段７１Ｂは、演算結果記憶手段６０に記憶された瞬時音響インテンシティをモニタにグラフィック表示する。ここで、音響インテンシティ表示手段７１Ｂは、図３の表示画面をＺ軸方向に拡張して、リスニングポイントごとの瞬時音響インテンシティを３次元のベクトルでグラフィック表示する（不図示）。
【０１１７】
到来音エネルギー表示手段７３Ｂは、演算結果記憶手段６０に記憶された到来音エネルギーをモニタにグラフィック表示する。ここで、到来音エネルギー表示手段７３Ｂは、図４の表示画面をＺ軸方向に拡張して、リスニングポイントごとの到来音エネルギーを３次元方向でグラフィック表示する（不図示）。
【０１１８】
＜瞬時音響インテンシティの算出＞
以下、図１０〜図１２を参照し、瞬時音響インテンシティの算出について説明する。
【０１１９】
図１０の例では、聴取者の頭部中心（不図示）を基準として、Ｘ軸、Ｙ軸およびＺ軸を設定した。また、スピーカＳＰは、半径３ｍの球体面上に位置している。図１０では、Ｘ軸方向で半径３ｍの円ｃｒ２、Ｙ軸方向で半径３ｍの円Ｃｒ１、Ｚ軸方向で半径３ｍの円Ｃｒ３を破線で図示した。
【０１２０】
図１０に示すように、リスニングエリアＬＡは、Ｘ軸方向に２ｍ×Ｙ軸方向に２ｍ×Ｚ軸方向に２ｍの大きさの立方体となっている。ここで、説明を簡易にするため、立方体のリスニングエリアＬＡのうち、Ｚ軸の値が聴取者の頭部中心の高さに一致するＸ−Ｙ平面上に該当する部分のみをリスニングエリアＬＡ’とする。
また、スピーカＳＰを１個とし、リスニングポイントＬＰが聴取者の頭部中心（スウィートスポット）であるとする。つまり、リスニングポイントＬＰは、リスニングエリアＬＡ’上に位置する。
【０１２１】
図１０のリスニングエリアＬＡ’において、リスニングポイントＬＰとスピーカＳＰとの位置関係を簡略化すると、図１１に示すようになる。この図１１では、なす角Φは、Ｘ―Ｙ平面と、スピーカＳＰの位置（Ｘｓｐｉ，Ｙｓｐｉ，Ｚｓｐｉ）からリスニングポイントＬＰの位置（Ｐｘ，Ｐｙ，Ｐｚ）までの線分とのなす角である。また、なす角θは、Ｙ−Ｚ平面と、スピーカＳＰの位置（Ｘｓｐｉ，Ｙｓｐｉ，Ｚｓｐｉ）からリスニングポイントＬＰの位置（Ｐｘ，Ｐｙ，Ｐｚ）までの線分とのなす角である。
【０１２２】
まず、音響インテンシティ算出手段４０Ｂは、リスニングポイントＬＰとスピーカＳＰとの位置関係から、スピーカＳＰの音響信号のＺ軸音圧成分（第３軸音圧成分）を算出する。図１１に示すように、音響インテンシティ算出手段４０Ｂは、音源信号の音圧をｐ（ｔ）とすると、スピーカＳＰについて、音響信号のＺ軸音圧成分ｓｉｎΦ＊｜ｐ（ｔ）｜を算出する。ここで、ベクトルＶ１は、スピーカＳＰの音響信号のＺ軸音圧成分を図示したものである。
【０１２３】
また、音響インテンシティ算出手段４０Ｂは、リスニングポイントＬＰとスピーカＳＰとの位置関係から、スピーカＳＰの音響信号のＸ軸音圧成分およびＹ軸音圧成分を算出する。このため、音響インテンシティ算出手段４０Ｂは、スピーカＳＰからＸ―Ｙ平面に到る垂線と、Ｘ―Ｙ平面との交点を求める。そして、音響インテンシティ算出手段４０Ｂは、この交点からリスニングポイントＬＰまでの線分と同じ向きで、このリスニングポイントＬＰを始点としたベクトルＶ２を求める。このベクトルＶ２は、ｃｏｓΦ＊｜ｐ（ｔ）｜と表される。
【０１２４】
そして、図１２に示すように、音響インテンシティ算出手段４０Ｂは、ベクトルＶ２を用いて、スピーカＳＰの音響信号のＸ軸音圧成分およびＹ軸音圧成分を算出する。この場合、音響インテンシティ表示手段７１Ｂは、スピーカＳＰについて、音響信号のＸ軸音圧成分ｃｏｓΦｃｏｓθ＊｜ｐ（ｔ）｜と、音響信号のＹ軸音圧成分ｃｏｓΦｓｉｎθ＊｜Ｐ（ｔ）｜とを算出する。ここで、ベクトルＶ３は、スピーカＳＰの音響信号のＹ軸音圧成分ｃｏｓΦｓｉｎθ＊｜ｐ（ｔ）｜と表される。また、ベクトルＶ４は、スピーカＳＰの音響信号のＸ軸音圧成分ｃｏｓΦｃｏｓθ＊｜ｐ（ｔ）｜と表される。
【０１２５】
次に、音響インテンシティ算出手段４０Ｂは、音響インテンシティ算出手段４０と同様、スピーカＳＰのＸ軸音圧成分、Ｙ軸音圧成分およびＹ軸音圧成分について、その符号（正または負）を求める。ここで、音響システム用デジタルスコープ１Ｂでは、リスニングポイントＬＰへの到来音を扱う。このため、音響インテンシティ算出手段４０Ｂは、Ｚ軸方向では、下方向（図１０，図１１下側）からリスニングポイントＬＰに到来する音を正（＋）の符号とし、上方向（図１０，図１１上側）からリスニングポイントＬＰに到来する音を負（−）の符号として累計する。例えば、音響インテンシティ算出手段４０Ｂは、スピーカＳＰがリスニングポイントＬＰの上側に位置するので、Ｚ軸音圧成分の符号を負（−）とする。その後、音響インテンシティ算出手段４０Ｂは、求めた符号を用いて、スピーカＳＰごとのＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とをそれぞれ累計して、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とを算出する。
【０１２６】
また、音響インテンシティ算出手段４０Ｂは、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分と軸音圧成分とに基づいて粒子速度のＸ軸成分とＸ軸成分とＺ軸成分とを算出し、スピーカＳＰごとの音源信号のＸ軸粒子速度成分とＹ軸粒子速度成分とＺ軸粒子速度成分（第３軸粒子速度成分）とをそれぞれ累計して、粒子速度のＸ軸成分累計値とＹ軸成分累計値とＺ軸成分累計値とを算出する。ここで、スピーカＳＰが１個であることから、音響インテンシティ算出手段４０Ｂは、下記の式（１１）に示すように、粒子速度のＸ軸成分累計値ｕｘ（ｔ）、Ｙ軸成分累計値ｕｙ（ｔ）およびＺ軸成分累計値ｕｚ（ｔ）を算出する。
【０１２７】
【数１１】
0005480748.tif 000012
【０１２８】
その後、音響インテンシティ算出手段４０Ｂは、粒子速度のＸ軸成分累計値ｕｘ（ｔ）、Ｙ軸成分累計値ｕｙ（ｔ）およびＺ軸成分累計値ｕｚ（ｔ）のそれぞれに音源信号の音圧を乗じて、瞬時音響インテンシティのＸ軸成分およびＹ軸成分を算出する。具体的には、音響インテンシティ算出手段４０Ｂは、下記の式（１２）に示すように、粒子速度のＸ軸成分累計値ｕｘ（ｔ）と音源信号の音圧ｐ（ｔ）とを乗じて瞬時音響インテンシティのＸ軸成分Ｉｘ（ｔ）を算出する。また、音響インテンシティ算出手段４０Ｂは、粒子速度のＹ軸成分累計値ｕｙ（ｔ）と音源信号の音圧ｐ（ｔ）とを乗じて瞬時音響インテンシティのＹ軸成分Ｉｙ（ｔ）を算出する。さらに、音響インテンシティ算出手段４０Ｂは、粒子速度のＺ軸成分累計値ｕｚ（ｔ）と音源信号の音圧ｐ（ｔ）とを乗じて瞬時音響インテンシティのＺ軸成分Ｉｚ（ｔ）を算出する。
【０１２９】
【数１２】
0005480748.tif 000013
【０１３０】
なお、スピーカＳＰが複数の場合、第１実施形態と同様に、音響インテンシティ算出手段４０Ｂは、スピーカＳＰごとに音源信号のＸ軸音圧成分、Ｙ軸音圧成分およびＺ軸音圧成分を求める。そして、音響インテンシティ算出手段４０Ｂは、これらＸ軸音圧成分、Ｙ軸音圧成分およびＺ軸音圧成分を累計し、リスニングポイントＬＰごとにＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とを算出することは言うまでもない。
また、Ｚ軸方向におけるリスニングポイントＬＰの位置が、聴取者の頭部の高さであるとして説明したが、これに限定されないことは言うまでもない。
【０１３１】
［音響システム用デジタルスコープの動作］
以下、図１３を参照して、音響システム用デジタルスコープ１Ｂの動作について説明する（適宜図１参照）。
【０１３２】
まず、音響システム用デジタルスコープ１Ｂは、パラメータ設定手段１０Ｂによって、音響システム用デジタルスコープ１Ｂに必要なパラメータが入力されると共に、入力されたパラメータをパラメータ記憶手段２０に記憶する（ステップＳ１１）。
【０１３３】
また、音響システム用デジタルスコープ１Ｂは、音源信号入力手段３０によって、Ｎ個の音源信号ファイルを読み込んで、音源信号ファイルのそれぞれに格納された音源信号のデータサイズｄを求める（ステップＳ１２）。そして、音響システム用デジタルスコープ１Ｂは、音源信号入力手段３０によって、求めたデータサイズｄだけ、音源信号ファイルのそれぞれから、Ｎ個の音源信号を読み込む（ステップＳ１３）。
【０１３４】
また、音響システム用デジタルスコープ１Ｂは、音響インテンシティ算出手段４０Ｂによって、入力された音源信号の音圧に粒子速度を乗じることで、リスニングポイントごとに瞬時音響インテンシティを算出する。このとき、音響システム用デジタルスコープ１Ｂは、音響インテンシティ算出手段４０Ｂによって、リスニングポイントとスピーカとの位置関係から、スピーカごとに音源信号のＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とを算出する。そして、音響システム用デジタルスコープ１Ｂは、音響インテンシティ算出手段４０Ｂによって、スピーカごとのＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とをそれぞれ累計して、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分とＺ軸音圧成分とを算出する。さらに、音響システム用デジタルスコープ１Ｂは、音響インテンシティ算出手段４０Ｂによって、リスニングポイントのＸ軸音圧成分とＹ軸音圧成分と軸音圧成分とに基づいて粒子速度のＸ軸成分とＹ軸成分とＺ軸成分とを算出し、スピーカごとの音源信号のＸ軸粒子速度成分とＹ軸粒子速度成分とＺ軸粒子速度成分とをそれぞれ累計して、粒子速度のＸ軸成分累計値とＹ軸成分累計値とＺ軸成分累計値とを算出する。その後、音響システム用デジタルスコープ１Ｂは、音響インテンシティ算出手段４０Ｂによって、粒子速度のＸ軸成分累計値とＹ軸成分累計値とＺ軸成分累計値とのそれぞれに音源信号の音圧を乗じて、瞬時音響インテンシティのＸ軸成分とＹ軸成分とＺ軸成分とを算出する（ステップＳ１４）。
【０１３５】
また、音響システム用デジタルスコープ１Ｂは、到来音エネルギー算出手段５０Ｂによって、瞬時音響インテンシティベクトルと、瞬時音響インテンシティベクトルの起点である他のリスニングポイントからリスニングポイントに向いた単位ベクトルとの内積によって、リスニングポイントごとに到来音エネルギーを算出する。つまり、音響システム用デジタルスコープ１Ｂは、到来音エネルギー算出手段５０Ｂによって、前記した式（１０）を用いて、リスニングポイントごとに到来音エネルギーを算出する（ステップＳ１５）。
【０１３６】
また、音響システム用デジタルスコープ１Ｂは、音響インテンシティ表示手段７１Ｂによって、瞬時音響インテンシティをモニタにグラフィック表示する。このとき、音響システム用デジタルスコープ１Ｂは、音響インテンシティ表示手段７１Ｂによって、瞬時音響インテンシティをリアルタイムでグラフィック表示する（ステップＳ１６）。
【０１３７】
また、音響システム用デジタルスコープ１Ｂは、到来音エネルギー表示手段７３Ｂによって、到来音エネルギーをモニタにグラフィック表示する。このとき、音響システム用デジタルスコープ１Ｂは、到来音エネルギー表示手段７３Ｂによって、到来音エネルギーをリアルタイムでグラフィック表示する（ステップＳ１７）。
【０１３８】
以上のように、本発明の第２実施形態に係る音響システム用デジタルスコープ１Ｂは、３次元のリスニングエリア内の様々な位置（リスニングポイント）において、到来音エネルギーを算出して表示する。これによって、音響システム用デジタルスコープ１Ｂは、作業者が、リスニングエリア内の様々な位置において、再生音場を客観的に把握することができる。
【０１３９】
また、音響システム用デジタルスコープ１Ｂは、実際の音場における実測データを用いて到来音エネルギーなどの音響情報を算出する従来技術に比べて、実測データを用いることなく、シミュレーションによって到来音エネルギーを算出するため、放送番組の制作現場で利用し易くなる。
【０１４０】
さらに、音響システム用デジタルスコープ１Ｂは、音響インテンシティと到来音エネルギーとをリアルタイムで算出してグラフィック表示するため、リスニングエリアにおける再生音場の経時変化を把握することができる。例えば、音源信号として楽曲を入力した場合、音響システム用デジタルスコープ１Ｂは、作業者が、楽曲の先頭部分、中間部分および終了部分などで、リスニングエリアにおける再生音場がどのように変化するか把握できる。
【０１４１】
なお、第１，２実施形態では、本発明に係る音響システム用デジタルスコープを独立した装置として説明したが、本発明では、一般的なコンピュータのハードウェアを、前記した各手段として協調動作させるプログラムによって実現できる。このプログラムは、通信回線を介して配布しても良く、ＣＤ−ＲＯＭやフラッシュメモリ等の記録媒体に書き込んで配布しても良い。
【符号の説明】
【０１４２】
１，１Ｂ  音響システム用デジタルスコープ
１０，１０Ｂ  パラメータ設定手段（パラメータ入力手段）
２０  パラメータ記憶手段
３０  音源信号入力手段
４０，４０Ｂ  音響インテンシティ算出手段
５０，５０Ｂ  到来音エネルギー算出手段
６０  演算結果記憶手段
７０Ｂ  表示手段
７１，７１Ｂ  音響インテンシティ表示手段
７３，７３Ｂ  到来音エネルギー表示手段
【図１】
0005480748.tif 000014
【図２】
0005480748.tif 000015
【図３】
0005480748.tif 000016
【図４】
0005480748.tif 000017
【図５】
0005480748.tif 000018
【図６】
0005480748.tif 000019
【図７】
0005480748.tif 000020
【図８】
0005480748.tif 000021
【図９】
0005480748.tif 000022
【図１０】
0005480748.tif 000023
【図１１】
0005480748.tif 000024
【図１２】
0005480748.tif 000025
【図１３】
0005480748.tif 000026