JP4783873B2 - 環境状態測定装置及び方法 - Google Patents
環境状態測定装置及び方法 Download PDFInfo
- Publication number
- JP4783873B2 JP4783873B2 JP2000244781A JP2000244781A JP4783873B2 JP 4783873 B2 JP4783873 B2 JP 4783873B2 JP 2000244781 A JP2000244781 A JP 2000244781A JP 2000244781 A JP2000244781 A JP 2000244781A JP 4783873 B2 JP4783873 B2 JP 4783873B2
- Authority
- JP
- Japan
- Prior art keywords
- sound wave
- propagation time
- sound
- measurement
- generator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Air Conditioning Control Device (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は音波発生検出装置、環境状態測定装置及び方法に係り、特に、環境状態の測定に利用可能な音波発生検出装置、前記音波発生検出装置を用いて空間内の温度分布等の環境状態を測定する環境状態測定方法、及び該環境状態測定方法を適用可能な環境状態測定装置に関する。
【０００２】
【従来の技術】
被空調空間内の温度や湿度を所望の値に制御するためには、被空調空間内の温度、湿度を測定する必要があるが、被空調空間内の各箇所における温度、湿度が一定であることは稀であり、特にドーム球場や劇場等の閉鎖された大空間では、空間内の各箇所における温度、湿度は大きくばらついていることが多い。このため、例えばドーム球場や劇場等の大空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、室内の温度や湿度を精密に所望の値に制御する場合には、被空調空間内の温度や湿度の分布を測定する必要がある。
【０００３】
しかし、温度センサや湿度センサによって温度や湿度の分布を直接測定しようとすると、センサを、被空調空間の中央部付近を含む被空調空間内にマトリクス状に多数配設する必要があるが、コストが嵩み景観が損なわれると共に被測定空間の利用も制限されるので、上記のように多数のセンサを配設することは現実的ではない。またドーム球場や劇場等では、被測定空間の周縁部（例えば壁体等）に取付けたセンサにより周縁部の温度等を測定すると共に、風船等でセンサを吊り上げることで被測定空間の中央部付近の温度等を測定することが考えられるが、この方法では、被測定空間内で例えば野球等のイベントが実施されている最中に温度や湿度の分布を測定することは不可能である。
【０００４】
一方、医学の分野では、人体に対して多方向からＸ線を照射し、人体を透過したＸ線量を高感度のセンサで測定し、測定されたＸ線量に基づき人体の横断面に沿った各位置における体内組織のＸ線吸収量をコンピュータにより演算して２次元画像として再構成し、人体の横断面像としてモニタに表示するコンピュータ断層撮影（ＣＴ：Computed Tomography)技術が広く利用されている。このＣＴ技術によれば、測定対象（この場合は人体）内部の状態を非接触で測定することができるので、このＣＴ技術を利用して被測定空間内の温度分布を求めることが提案されている。
【０００５】
具体的には、Ｘ線に代えて、媒質の温度によって伝搬時間（速度）が変化する音波を用い、例えば被測定空間を区画している壁体の多数箇所にスピーカ、マイクロフォン、温度センサを設け、スピーカ及びマイクロフォンによって前記多数箇所の間の音波の伝搬時間を各々測定し、被測定空間内の多数箇所の間の音波の伝搬時間から、被測定空間内の温度分布を演算によって求めている。
【０００６】
【発明が解決しようとする課題】
音波の伝播時間の測定は、詳しくはデジタルの原音声データをＤ／Ａ変換器を介してスピーカに供給することでスピーカから音波を発生させ、被測定空間を伝播した音波がマイクロフォン配設位置に到達することでマイクロフォンから出力された音声信号を、Ａ／Ｄ変換器によってデジタルの検出音声データへ変換し、検出音声データの波形と原音声データの波形を時間的にずらしながら相関性を演算し、相関性が最も高くなるときの時間差を音波の伝播時間とすることによって行われる。
【０００７】
しかしながら、スピーカは音声信号に応じて振動板を振動させることで音を発する構造であり、マイクロフォンは周囲の音に応じて発生する振動板の振動を音声信号に変換する構造であるので、振動板の質量の影響により、スピーカから発せられる音の波形及びマイクロフォンから出力される音声信号の波形は、何れも元の波形に対して鈍っていると共に、周波数が高くなるに従って振動板の振動が追従し切れなくなることで所定周波数以上の高周波成分がカットされる。また、スピーカとマイクロフォンは扱う音の音圧レベルが異なっており、より音圧レベルの高い音を扱うスピーカの振動板は、相応の強度が必要とされるために質量が大きく、上記の現象がより顕著である。
【０００８】
このため、スピーカから発せられた音をマイクロフォンで検出することによって得られる音声信号の波形は、スピーカ及びマイクロフォン（特にスピーカ）の影響により原音声データの波形と相違するので、伝播時間の測定における測定精度低下の原因となっており、温度分布の測定精度の低下も招いていた。また、音声信号に対して振動板の振動のレスポンスの悪いスピーカは、マイクロフォンから出力される音声信号の波形の鈍りを助長し、伝播時間及び温度分布の測定精度の更なる低下に繋がるために使用できないという問題もあった。
【０００９】
また、スピーカの振動板は例えばコーン型やドーム型等の立体的な形状である（平面型は波形再現性が劣るという欠点がある）と共に、振動板上における実際に音が発せられる点（発音点）の位置は、振動板の分割振動により、スピーカから発せられる音の周波数によって相違するため、スピーカの発音点とマイクロフォンとの距離もスピーカから発せられる音の周波数によって若干変化する。このスピーカの発音点とマイクロフォンとの距離の変化は、特に狭い空間内における伝播時間及び温度分布の測定に無視できない影響を及ぼすので、伝播時間及び温度分布の測定のためにスピーカから発する音が、例えば単一周波数の単純音又はそれに近い音に制限されていた。
【００１０】
また、３台のスピーカ及びマイクロフォンを用いて発音点の位置を計測する３スピーカ法等のキャリブレーション方法も提案されているが、計測精度が充分でないという問題があった。
【００１２】
本発明は、被測定空間内の温度分布等の環境状態を高精度に測定できる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１３】
【課題を解決するための手段】
上記目的を達成するために請求項１記載の発明に係る環境状態測定装置は、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられ、被測定空間を横切る複数の伝播時間測定区間の両端部に、前記ブラケットが取付けられた側同士が対向するように各々配置された音波発生検出装置と、伝播時間測定区間の両端部に配置された一対の前記音波発生検出装置のうちの一方の前記音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の前記音波発生検出装置の前記第１音波検出手段、及び他方の前記音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行う制御手段と、前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行う伝播時間演算手段と、前記伝播時間測定手段によって測定された前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算する環境状態演算手段と、を含んで構成されている。
【００１４】
請求項１記載の発明は、音波発生検出装置を備えており、この音波発生検出装置は、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含んで構成されている。
【００１５】
上記の音波発生器としては、例えばスピーカを用いることができ、詳しくは動電形（例えばダイナミックスピーカ）、電磁形（例えばマグネティックスピーカ）、静電形（例えばコンデンサスピーカ）等の各種スピーカを適用することができるが、音波発生器が発生する音波は可聴域外の周波数の音波であってもよく、例えば音波発生器として超音波を発する送波器等を用いることも可能である。
【００１６】
第１音波検出手段及び第２音波検出手段としては、例えばマイクロフォンを用いることができ、詳しくは動電形（例えばダイナミックマイクロフォンやリボンマイクロフォン）、静電形（例えばコンデンサマイクロフォン）、圧電形（例えばクリスタルマイクロフォン）等の各種マイクロフォン等の音波検出器を適用することができるが、音波発生器が発生する音波が超音波であれば、例えば超音波を検出する受波器等の音波検出器を用いることも可能である。
【００１７】
また、第１音波検出手段と第２音波検出手段は別体の音波検出器である必要はない。すなわち、第１音波検出手段と第２音波検出手段は検出対象の音波の到来方向がおおよそ逆方向であるので、例えば第１音波検出手段及び第２音波検出手段を、両指向性の単一の音波検出器によって構成することも可能である。更に、第１音波検出手段と第２音波検出手段を各々独立した単一指向性の音波検出器によって構成する場合には、第１音波検出手段に対応する音波検出器は音波検出方向が音波発生器側を向くように配置し、第２音波検出手段に対応する音波検出器は音波検出方向が音波発生器と反対側を向くように配置すればよい。
【００１８】
また、上記の音波発生検出装置は、第１音波検出手段、第２音波検出手段及び光反射板が、音波発生器と連結されたブラケットに一体に取付けられ、被測定空間を横切る複数の伝播時間測定区間の両端部に、ブラケットが取付けられた側同士が対向するように各々配置されている。これにより、伝播時間測定区間の両端部に各々配置された一対の音波発生検出装置のうち、一方の音波発生検出装置（便宜的に「第１の音波発生検出装置」という）の音波発生器から音波を発生させれば、該発生させた音波を、第１の音波発生検出装置の第１音波検出手段により、前記区間の前記音波発生器側の端部の所定位置（詳しくは前記音波発生器から所定距離隔てた位置：便宜的に「第１の位置」という）で検出すると共に、他方の音波発生検出装置（便宜的に「第２の音波発生検出装置」という）の第２音波検出手段により、前記区間を挟んで前記音波発生器と反対側の端部の所定位置（便宜的に「第２の位置」という）で検出することができる。
【００１９】
ここで、例えば音波発生器が、入力された音声信号に応じた音波を発生する構成であると共に、発生される音波の波形が音声信号の波形から変化する（例えば波形が鈍る）特性を有しているとしても、音波発生器によって発生された音波を検出する第１の音波発生検出装置の第１音波検出手段及び第２の音波発生検出装置の第２音波検出手段は、音波発生器から発せられた音波（音波発生器によって波形が変化した後の音波）を検出するので、第１音波検出手段によって検出される音波の波形と第２音波検出手段によって検出される音波の波形は高い精度で一致する。
【００２０】
また、音波発生器が音波を発生する際の発音点の位置は、音波発生器が発生する音波の周波数に応じて変化するので、この発音点の位置の変化が第２の音波発生検出装置の第２音波検出手段による音波の検出に影響を及ぼし、第１の音波発生検出装置の音波発生器から発生されて第２の音波発生検出装置の第２音波検出手段で検出される音波の伝播距離は、発音点の位置の変化に応じて変化し、この音波の伝播距離の変化に伴って音波の伝播時間も変化する。
【００２１】
しかしながら、第１の音波発生検出装置の第１音波検出手段は第２の音波発生検出装置の第２音波検出手段と同一の音波を検出するので、上述した発音点の位置の移動は、第１の音波発生検出装置の第１音波検出手段による音波の検出に対しても同様の影響を及ぼし、第１の音波発生検出装置の音波発生器から発生されて第１の音波発生検出装置の第１音波検出手段で検出される音波の伝播距離及び伝播時間は、第２の音波発生検出装置の第２音波検出手段による音波の検出と同様に変化する。
【００２２】
従って、第１の音波発生検出装置の第１音波検出手段によって検出された音波の波形と、第２の音波発生検出装置の第２音波検出手段によって検出された音波の波形を照合し（例えば相関性を演算する等）、伝播時間を測定すべき区間の音波の伝播時間として、音波発生器によって発生された音波が第１の音波発生検出装置の第１音波検出手段に到達してから、該音波が第２の音波発生検出装置の第２音波検出手段に到達する迄の時間差、すなわち第１の音波発生検出装置の第１音波検出手段と、第２の音波発生検出装置の第２音波検出手段の間の音波の伝播時間を求めることにより、音波の伝播時間を、音波発生器の特性や発音点の位置の移動の影響を受けることなく高精度に測定することができる。
【００２３】
また、複数の区間の音波の伝播時間を測定することで温度分布等の環境状態を測定するに際しては、伝播時間測定区間の距離を正確に把握することが重要であるが、伝播時間測定区間として音波発生器から音波検出器に至る区間を用いた場合、発音点の位置が変化すると音波の実際の伝播距離が変化するため、伝播時間測定区間の距離を正確に把握することは困難である。これに対し、請求項１記載の発明に係る音波発生検出装置を用いた音波の伝播時間の測定では、第１の音波発生検出装置の第１音波検出手段と第２の音波発生検出装置の第２音波検出手段の間の区間が伝播時間測定区間となるため、発音点の位置が変化しても音波の実際の伝播距離が変化することはない。
【００２４】
請求項１記載の発明において、制御手段は、伝播時間測定区間の両端部に配置された一対の音波発生検出装置のうちの一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、一方の音波発生検出装置の第１音波検出手段、及び他方の音波発生検出装置の第２音波検出手段によって各々検出させることを、複数の伝播時間測定区間について各々行う。
【００２５】
また伝播時間演算手段は、第２音波検出手段によって検出された音波の波形を第１音波検出手段によって検出された音波の波形と照合し、波形の相関性を時間軸に沿って演算した結果に基づいて、第１音波検出手段による音波の検出位置から第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、複数の伝播時間測定区間について各々行う。これにより、音波発生器の特性の影響や発音点の位置が移動することによる影響を受けることなく、複数の伝播時間測定区間について、音波の伝播時間を各々精度良く測定することができる。
【００２６】
また、上記のように音波の伝播時間の測定精度が発音点位置の移動の影響を受けることがないため、音波発生器から発生させる音波として任意の音波（例えば自然音に相当する音波）を用いることが可能となる。これにより、被測定空間内に人間が存在している場合にも、伝播時間の測定にあたり被測定空間内に存在している人間に違和感や不快感を与えることを回避することができる。
【００２７】
そして環境状態演算手段は、伝播時間測定手段によって測定された複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて、被測定空間内の環境状態を演算するので、被測定空間内の温度分布等の環境状態を高精度に測定することができる。
【００２８】
また、音波発生検出装置は、音波発生器と連結されたブラケットに第１音波検出手段、第２音波検出手段及び光反射板が一体に取付けられているので、音波発生検出装置の設置に際し、音波発生器から距離が等しくなるように第１音波検出手段、第２音波検出手段及び光反射板の相対位置を調整する必要がなくなり、音波発生検出装置の設置を容易に行うことができると共に、周囲温度の変化に伴う第１音波検出手段、第２音波検出手段及び光反射板の相対位置の変化を抑制することができる。
【００２９】
更に、音波発生検出装置のブラケットには、音波発生器側と反対側より到来した光を、第１音波検出手段及び第２音波検出手段に対して音波発生器からの距離が等しい位置（音波発生器から所定距離隔てた位置）で反射する光反射板が設けられている。この光反射板は、例えば光反射面が音波発生器と反対側を向くように配置した全反射ミラーによって構成することができる。この光反射板を利用することで、伝播時間測定区間の距離を、例えばレーザ距離計等によって容易かつ正確に測定することができ、音波の伝播時間の測定結果及び伝播時間測定区間の距離の測定結果に基づいて、温度分布等の環境状態を精度良く測定することができる。
【００３０】
なお、環境状態演算手段による被測定空間内の環境状態の演算は、詳しくは、例えば被測定空間を多数の小領域に区分し、各小領域内の音速を変数とし、或る音波発生装置から或る音波検出装置に至る伝搬経路を伝搬する音波の伝搬時間を、前記伝搬経路上に存在する各小領域の音波の通過時間の和で表した方程式を連立方程式として多数生成し、測定した音波の伝搬時間を連立方程式に代入して解くことによっても実現できるが、下記のようにして演算することが好ましい。
【００３１】
すなわち、各々一定の方向に沿いかつ互いに異なる位置で被測定空間を貫く等しい長さの多数の仮想経路を、互いに異なる複数の方向について各々設定しておき、前記各仮想経路について、第１の音波発生検出装置の第１音波検出手段と、第２の音波発生検出装置の第２音波検出手段と、の間に相当する仮想経路上の測定区間の音波の伝搬時間を実際に音波を伝搬させて測定し、仮想経路上の測定区間外の区間の音波の伝搬時間を被測定空間外が一定温度の媒質で満たされていると仮定して推定することにより、仮想経路上を音波が伝搬するときの伝搬時間の推定値を各々求め、所定方向に沿った多数の仮想経路の各々の音波の伝搬時間の推定値に基づき、前記所定方向についての音波の伝搬時間の投影データを求めることを、前記複数の方向について各々行い、前記複数の方向の各々についての音波の伝搬時間の投影データに基づいて被測定空間内の環境状態を演算することが好ましい。
【００３２】
上記では、伝搬時間測定経路の長さを仮想的に一致させることで、音波の伝搬時間の投影データを取得できるので、短時間で演算可能な演算方法（フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等）を適用することができ、演算時間の短縮や伝播時間を測定するための装置構成の簡易化を実現できる。
【００３３】
なお、請求項１記載の発明に係る第１音波検出手段は、検出対象の音波を発生する音波発生器が比較的近い位置に位置しているので、検出対象の音波の音圧レベルが比較的高く、第１音波検出手段によって検出される音波の波形が歪んだり飽和したりする可能性がある。これを考慮すると、請求項２に記載したように、第１音波検出手段による音波検出位置の音波到来側に、到来した音波を減衰させる減衰手段を配置することが好ましい。なお、減衰手段としては、例えば吸音材等を用いることができる。これにより、第１音波検出手段によって検出される音波の波形が歪んだり飽和することを防止することができる。
【００３４】
また、音波検出手段によって検出される音波の波形は、到来した音波の波形に対して若干ではあるが変化する（例えば鈍りが生ずる等）。これを考慮すると、請求項３に記載したように、第１音波検出手段及び第２音波検出手段を、少なくとも類似種（好ましくは同一種）の音波検出器によって音波を検出するよう構成することが好ましい。これにより、第１音波検出手段によって検出される波形及び第２音波検出手段によって検出される波形が、到来した音波の波形に対して同様の変化を示すことになる。
【００３５】
従って、本発明に係る音波発生検出装置を、音波の伝播時間を測定すべき区間の両端部に各々配置する態様において、音波発生器から音波を発生させる第１の音波発生検出装置の第１音波検出手段によって検出される音波の波形と、第２の音波発生検出装置の第２音波検出手段によって検出される音波の波形を高い精度で一致させることができ、音波伝播時間をより精度良く測定することができる。
【００３８】
請求項４記載の発明に係る環境状態測定方法は、被測定空間を横切る複数の伝播時間測定区間の両端部に、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられた音波発生検出装置を、前記ブラケットが取付けられた側同士が対向するように各々配置しておき、伝播時間測定区間の両端部に配置した一対の音波発生検出装置のうちの一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の音波発生検出装置の前記第１音波検出手段、及び他方の音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行い、前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行い、前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算するので、請求項１記載の発明と同様に、被測定空間内の温度分布等の環境状態を高精度に測定することができる。
【００３９】
また、請求項５記載の発明は、請求項４記載の発明において、伝播時間測定区間の両端部に各々配置した一対の音波発生検出装置の一対の光反射板に各々光を照射して前記一対の光反射板の距離を測定することを、前記複数の伝播時間測定区間について各々行い、前記被測定空間内の環境状態の演算を、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して行うことを特徴としている。
【００４０】
請求項５記載の発明では、伝播時間測定区間の両端部に各々配置した一対の音波発生検出装置の一対の光反射板に各々光を照射して一対の光反射板の距離を測定する。なお、この距離測定には例えばレーザ距離計等を用いることができ、光反射板にレーザ光を照射することで一対の光反射板の距離、すなわち伝播時間測定区間の距離を容易かつ正確に測定することができる。
【００４１】
そして請求項５記載の発明では、上記の距離の測定を複数の伝播時間測定区間について各々定期的に行い、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して、被測定空間内の環境状態の演算を行うので、例えば季節変化等に起因する環境温度の変化や経時変化に伴って伝播時間測定区間の距離が変化したとしても、距離変化の影響を受けることなく、被測定空間内の環境状態を精度良く測定することができる。
【００４２】
【発明の実施の形態】
以下、図面を参照して本発明の実施形態の一例を詳細に説明する。図１には、環境状態測定方法が適用された環境状態測定装置１０が示されている。環境状態測定装置１０は、多数のノード１２Ａ、１２Ｂ、１２Ｃ、…と、ホストコンピュータ１４を備えており、これらが伝送媒体１８を介して互いに接続されて構成されている。
【００４３】
ホストコンピュータ１４は、ＣＰＵ１４Ａ、ＲＯＭ１４Ｂ、ＲＡＭ１４Ｃ、入出力ポート１４Ｄを備え、これらがバス１４Ｅを介して互いに接続されて構成されている。入出力ポート１４Ｄはネットワーク伝送部１６を介して伝送媒体１８に接続されている。また、入出力ポート１４Ｄには、各種の情報を表示するためのディスプレイ２０と、オペレータが各種のデータやコマンド等を入力するためのキーボード２２が接続されている。
【００４４】
一方、多数のノード１２は、被測定空間の周縁部に、周縁部の全周に亘って略一定の間隔で配置されている。例として図２には、円筒状の壁体２２によって区画され、上方が屋根（図示省略）によって閉塞されたドーム状の被測定空間２４（例えばドーム球場等）内の環境状態を測定する場合の各ノード１２の配置が示されている。なお、ノード１２の数及び間隔は図２に示した例に限定されるものではないが、ノード１２の数は偶数であることが望ましい。また、各ノード１２を識別するために、各ノードには便宜的にノード番号が付与されている。
【００４５】
ノード１２Ａ、１２Ｂ、１２Ｃ、…は各々同一の構成であるので、以下ではノード１２Ａについて説明する。ノード１２Ａは、被測定空間２４内に向けて配置された音波発生器としてのスピーカ３０と、スピーカ３０へ向けて配置された第１音波検出手段としての第１マイクロフォン３２と、被測定空間２４内に向けて配置された第２音波検出手段としての第２マイクロフォン３４を備えている。スピーカ３０は壁体２２に埋設されており（図示省略）、増幅器３６、Ｄ／Ａ変換器３８を介してデータ処理部４０に接続されている。
【００４６】
図３（Ａ）に示すように、スピーカ３０には、中間部が被測定空間２４側へ突出するようにＵ字型に屈曲されたアーム５４の基部が取付けられている。図３（Ｂ）に示すように、アーム５４の中間部には、スピーカ３０と対向する面に孔５４Ａが穿設されており、第１マイクロフォン３２は、スピーカ３０から発せられた音波をスピーカ３０から所定距離隔てた位置で検出するように、孔５４Ａの内部に埋設されている。なお、以下では第１マイクロフォン３２の配置位置を「第１の位置」と称する。
【００４７】
第１マイクロフォン３２は自ノードのスピーカ３０から発せられた音波を検出するので、音源（スピーカ３０）との距離が小さく、検出対象の音波の音圧レベルが比較的高い。このため、第１マイクロフォン３２は孔５４Ａの底部側に配置されており、孔５４Ａの手前側（第１マイクロフォン３２とスピーカ３０の間）には吸音材５６が設けられている。これにより、第１マイクロフォン３２から出力される音声信号の信号レベルが歪んだり飽和することが防止される。吸音材５６は請求項２に記載の減衰手段に対応している。
【００４８】
また、図３（Ｂ）に示すように、アーム５４の中間部には、孔５４Ａが穿設された面と反対側の面（被測定空間２４側を向いている面）に孔５４Ｂ、５４Ｃが各々穿設されており、第２マイクロフォン３４は、他ノードのスピーカ３０から発せられた音波を、スピーカ３０から所定距離隔てた位置で検出するように、孔５４Ｂの内部に埋設されている。また、以下では第２マイクロフォン３４の配置位置を「第２の位置」と称する。第２マイクロフォン３４は他ノードのスピーカ３０から発せられた音波を検出するので、音源との距離が大きく検出対象の音波の音圧レベルが比較的小さいので、孔５４Ｂには吸音材は設けられていない。
【００４９】
なお、本実施形態では第１マイクロフォン３２及び第２マイクロフォン３４として同一機種のマイクロフォンを用いており、第１マイクロフォン３２及び第２マイクロフォン３４は請求項３に記載の「類似種の音波検出器」に対応している。一方、孔５４Ｃの底部には、スピーカ３０側と反対側より到来した光をスピーカ３０から所定距離隔てた位置で反射するように光反射板５８が取付けられている。光反射板５８は各ノード１２のアーム５４間の距離をレーザ距離計等によって測定する際に用いられる。
【００５０】
なお、上述したスピーカ３０、アーム５４、第１マイクロフォン３２、第２マイクロフォン３４及び光反射板５８は、請求項１に記載の音波発生検出装置に対応している。
【００５１】
第１マイクロフォン３２はフィルタ／増幅器４２、Ａ／Ｄ変換器４４を介してデータ処理部４０に接続されており、第２マイクロフォン３４はフィルタ／増幅器４６、Ａ／Ｄ変換器４８を介してデータ処理部４０に接続されている。データ処理部４０はマイクロコンピュータ等を含んで構成されており、ネットワーク伝送部５０を介して伝送媒体１８に接続されている。
【００５２】
次に本実施形態の作用を説明する。本実施形態では、後述する温度分布演算処理において、音波の伝搬時間の推定を行う経路として、例として図４に示すように、多数の仮想経路を予め設定している。なお、図４では或る一定の方向に沿った仮想経路のみを示しているが、実際には、一例として図６（Ｂ）にも示すように、互いに異なる複数の方向（例えば１０°ずつ異なる方向）について、各々多数の仮想経路が予め設定されている。
【００５３】
仮想経路は、被測定空間２４外に、所定方向に直交するように配置された仮想スピーカ列上に位置している仮想スピーカと、被測定空間２４を挟んで仮想スピーカ列と反対側の被測定空間２４外に一定の間隔Ｌ0 （図５参照）を隔てて仮想スピーカ列と平行に配置された仮想マイクロフォン列上に位置している仮想マイクロフォンと、を結ぶ経路として定義されており、仮想スピーカ及び仮想マイクロフォンは、各仮想経路が、各々所定方向に沿いかつ互いに異なるノード１２が配置されている位置で壁体２２と交差して被測定空間２４を貫くように、仮想スピーカ列及び仮想マイクロフォン列上での位置が定められている。従って、各仮想経路の長さは互いに等しくされている。
【００５４】
各仮想経路のうち、同一の仮想経路上に位置している一対のノード１２の間の区間（図４に実線で示す測定対象区間：詳しくは一対のノード１２の各々のアーム５４の間の区間）については、実際に音波を伝搬させることによって音波の伝搬時間を測定するが、仮想スピーカとノード１２の間、及びノード１２と仮想マイクロフォンの間の区間（伝搬時間推定区間）については、音波の伝搬時間を測定することはできない。このため、ホストコンピュータ１４のＲＯＭ１４Ｂには、各仮想経路に対応して伝搬時間推定区間の長さが各々記憶されている。
【００５５】
なお、本実施形態では、図５に示すように、被測定空間２４の中心を原点Ｏとする直交座標系ｘ−ｙを定義し、仮想経路の各々を、仮想経路に直交し原点Ｏを通る直線Ｓの長さＸ（距離Ｘ）、及び前記直線Ｓとｘ軸との成す角度θ（傾きθ）によって識別する（仮想経路（Ｘ，θ）と称する）。
【００５６】
次に図７のフローチャートを参照し、被測定空間２４内の環境状態（本実施形態では温度分布）を求める場合にホストコンピュータ１４で実行される温度分布演算処理について説明する。
【００５７】
本実施形態では、伝播時間の測定に用いる音波データ（音波をアナログ信号に変換し更にデジタルデータに変換したもの）がホストコンピュータ１４のＲＯＭ１４Ｂ等に予め記憶されており、ステップ１００では、伝送媒体１８を介して前記音波データを全てのノード１２へ転送する。各ノード１２に転送された音波データは、各ノード１２のデータ処理部４０の内蔵メモリ等に記憶される。
【００５８】
なお、この音波データとしては任意の音声を表す音波データを用いることができ、例えば音楽、メッセージ、自然音（鳥のさえずりやドアの開閉音等）を表す音波データを用いることができる。また、被測定空間２４の大きさが大きくなるに従って（すなわち被測定空間１２の中央部を挟んで対向するノード１２間の距離が大きくなるに従って）、音波の周波数帯域が低周波側に偏倚するように音波データを選択してもよい。
【００５９】
次のステップ１０２では音波を発生させるノード１２のノード番号を表すカウント値ｍに１を代入する。ステップ１０４では音波発生対象のノードｍからの音波発生時刻を決定し、次のステップ１０６では、決定した音波発生時刻及び音波発生対象のノード１２のノード番号を全てのノード１２に通知する。次のステップ１０８では、音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。音波発生タイミングが到来すると、ステップ１０８の判定が肯定されてステップ１１０へ移行し、音波発生対象のノードｍに対し、先に転送した音波データを用いて音波を発生するよう指示する。
【００６０】
音波発生対象のノードｍのデータ処理部４０では、図８（Ａ）に示す音波発生処理が行われる。この音波発生処理は、ホストコンピュータ１４から通知された音波発生対象のノード１２のノード番号に基づき、自ノードが音波発生対象のノード１２であることを認識すると実行され、まずステップ１５０で音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。
【００６１】
伝送媒体１８及びネットワーク伝送部４８を介してホストコンピュータ１４から音波の発生が指示されると、ステップ１５０の判定が肯定されてステップ１５２へ移行し、先に転送された音波データをＤ／Ａ変換器３８に出力する。データ処理部４０から出力された音波データは、Ｄ／Ａ変換器３８でアナログの音声信号に変換され、増幅器３６によって所定の増幅率で増幅された後にスピーカ３０に供給される。これにより、ノードｍのスピーカ３０からは、温度分布演算処理（図７）のステップ１０４で判定された音波発生タイミングで、音波データが表す音波が所定の音量で発せられる。
【００６２】
音波発生対象のノードｍのスピーカ３０によって発生された音波はノードｍの第１マイクロフォン３２によって検出され、ノードｍの第１マイクロフォン３２からアナログの音声信号が出力される。ステップ１５２では、スピーカ３０から音波を発生させる処理と並行してＡ／Ｄ変換器４４を所定時間（例えば数秒）作動させ、第１マイクロフォン３２が音波を検出することで第１マイクロフォン３２から出力される音声信号のサンプリングを所定時間実行させる。
【００６３】
これにより、第１マイクロフォン３２から出力された音声信号は、所定時間の間、フィルタ／増幅器４２で雑音が除去されて増幅された後に、Ａ／Ｄ変換器４４によって一定のサンプリング周期でサンプリングされてデジタルの音波データに変換され、所定時間分の比較用音波測定データ（第１マイクロフォン３２で検出された音波の振幅を、所定時間に亘り（ｔWAVE＝０〜ｔLAST：図９（Ｂ）参照）、Ａ／Ｄ変換器４２のサンプリング周期に相当する時間間隔で時系列に表すデータ）としてデータ処理部４０に順次入力される。データ処理部４０は順次入力される比較用音波測定データをメモリ等に記憶する。
【００６４】
ステップ１５４では、上記の処理によって得られた所定時間分の比較用音波測定データを自ノード以外の各ノード１２へ転送し、音波発生処理を終了する。
【００６５】
例として、図９（Ａ）には音波発生対象のノード１２からの音波の発生に用いた音波データが表す音波の波形の一例を、図９（Ｂ）には比較用音波測定データが表す音波の波形の一例を各々示す。比較用音波測定データは、スピーカ３０から発せられた音波を第１マイクロフォン３２によって検出した結果を表しているので、図９（Ａ）の波形と図９（Ｃ）の波形を比較しても明らかなように、比較用音波測定データが表す音波の波形は、スピーカ３０及び第１マイクロフォン３２の特性の影響を受け、元の音波データが表す音波の波形に対して若干変化している（歪んでいる）。
【００６６】
一方、音波発生対象のノード１２以外のノード１２（測定対象のノード１２）のデータ処理部４０では、図８（Ｂ）に示す伝播時間演算処理が行われる。この伝播時間演算処理は、ホストコンピュータ１４から通知された音波発生対象のノード１２のノード番号に基づき、自ノードが音波発生対象のノード１２でないことを認識すると実行され、まずステップ１６０では、ホストコンピュータ１４からノード番号と共に通知された音波発生時刻を記憶する。次のステップ１６２では、記憶した音波発生時刻を現時刻と比較することで音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。音波発生時刻が到来すると、上記判定が肯定されてステップ１６４へ移行する。
【００６７】
先にも説明したように音波発生時刻が到来すると、音波発生対象のノード１２のスピーカから音波が発せられるが、この音波は被測定空間２４全域に放射され、音波測定対象の各ノード１２の第２マイクロフォン３４によって各々検出され、音波測定対象の各ノード１２の第２マイクロフォン３４からアナログの音声信号が出力される。ステップ１６４ではＡ／Ｄ変換器４８を所定時間（例えば数秒）作動させ、第２マイクロフォン３４が音波を検出することで第２マイクロフォン３４から出力される音声信号のサンプリングを所定時間実行させる。
【００６８】
これにより、第２マイクロフォン３２から出力された音声信号は、所定時間の間、フィルタ／増幅器４６で雑音が除去されて増幅された後に、Ａ／Ｄ変換器４８によって一定のサンプリング周期でサンプリングされてデジタルの音波データに変換され、所定時間分の音波測定データ（第２マイクロフォン３２で検出された音波の振幅を、所定時間に亘り、Ａ／Ｄ変換器４２のサンプリング周期に相当する時間間隔で時系列に表すデータ）としてデータ処理部４０に順次入力される。データ処理部４０は順次入力される音波測定データをメモリ等に記憶する。なお、音波測定データのデータ量を削減するために、音波発生時刻が到来してから或る時間（音波発生対象のノード１２からの音波の伝播時間よりも若干短い時間）経過後に音声信号のサンプリングを開始するようにしてもよい。
【００６９】
次のステップ１６６では音波発生対象のノード１２から比較用音波測定データを受信したか否か判定し、判定が肯定される迄待機する。比較用音波測定データを受信すると、ステップ１６６の判定が肯定されてステップ１６８へ移行し、音波測定データが表す波形と音波発生対象のノード１２から受信した比較用音波測定データが表す波形との相関性を演算する。
【００７０】
音波測定データが表す波形と比較用音波測定データが表す波形との相関性は、音波発生対象のノード１２の第１マイクロフォン３２の配設位置（第１の位置）に測定対象音波が到達したタイミングを基準とし、この基準タイミングから、測定対象のノード１２の第２マイクロフォン３４の配設位置（第２の位置）に測定対象音波が到達する迄の時間（音波の伝播時間）として或る値ｔX を仮定し、ｔWAVEが所定値（初期値はｔWAVE＝０：基準タイミング（図９（Ｂ）参照))のときの振幅値を表すデータを比較用音波測定データから取り出すと共に、対応するタイミング（基準タイミングからｔX ＋ｔWAVE経過したタイミング）における振幅値を表すデータを音波測定データから取り出し、両者を乗算することを、ｔWAVEの値を０〜ｔLAST迄順に変化させながら繰り返し、乗算結果を積算することで求めることができる。
【００７１】
上記の演算によって求まる相関性の値は、音波発生対象のノード１２のアーム５４の配設位置（第１マイクロフォン３２の配設位置）に測定対象音波が到達してから時間ｔX が経過したときに、測定対象のノード１２のアーム５４の配設位置（第２マイクロフォン３４の配設位置）に測定対象音波が到達した確からしさを表しており（相関性の値が符号が正でかつ絶対値が大きくなる程、前記確からしさが高いとみなすことができる）、上記の演算処理を、時間ｔX の値を徐々に変化させながら繰り返すことで、例として図９（Ｄ）に示すように、音波測定データが表す波形と比較用音波測定データが表す波形との相関性の値の時間軸に沿った推移を表すデータを得ることができる。
【００７２】
次のステップ１７０では、ステップ１６８における相関性の演算結果に基づいて、相関性の値が最大になるときの時間ｔX の値を求め、音波発生対象のノード１２から測定対象のノード１２（自ノード）への測定対象音波の伝播時間（詳しくは音波発生対象のノード１２の第１の位置から測定対象のノード１２の第２の位置への測定対象音波の伝播時間）として、相関性の値が最大のときの時間ｔX を設定し、設定した伝播時間を、ノードｍからノードｎへの音波の伝播時間データｔmnとしてホストコンピュータ１４へ転送し、伝播時間演算処理を終了する。
【００７３】
例として、図９（Ｃ）には音波測定データが表す音波の波形の一例を示す。音波測定データは、スピーカ３０から発せられた音波を第２マイクロフォン３４によって検出した結果を表しているので、図９（Ａ）の波形と図９（Ｃ）の波形を比較しても明らかなように、音波測定データが表す音波の波形は、スピーカ３０及び第２マイクロフォン３４の特性の影響を受け、元の音波データが表す音波の波形に対して若干変化している（歪んでいる）。
【００７４】
しかしながら、比較用音波測定データは、同一のスピーカ３０から発せられた音波を、第２マイクロフォン３４と同一機種の第１マイクロフォン３２（すなわち、第２マイクロフォン３４と特性が同一のマイクロフォン）によって検出した結果を表しているので、音波測定データが表す音波の波形及び比較用音波測定データが表す音波の波形は、図９（Ｂ）の波形と図９（Ｃ）の波形を比較しても明らかなように、元の音波データが表す音波の波形に対して同様に変化しており、波形の相関度が飛躍的に向上している。
【００７５】
従って、ステップ１６８の演算によって求まる相関性の値は、時間ｔX の値が音波発生対象のノード１２から測定対象のノード１２への測定対象音波の伝播時間に一致したときに正の最大値となり、時間ｔX の値がその前後の値のときにも相関性の値の絶対値は多少大きいものの、時間ｔX の値が前記伝播時間と大きく異なる値のときには相関性の値の絶対値は非常に小さな値となる。これにより、スピーカ３０や第２マイクロフォン３４の特性の影響を受けることなく、伝播時間を高精度に測定・設定することができる。また、各ノード１２のアーム５４の間を伝播速度の測定対象区間としているので、測定対象区間の音波伝播時間の測定結果が、スピーカ３０から発せられる音波の周波数変化に伴う発音点位置の変化の影響を受けて変動することもない。なお、ステップ１６８，１７０は請求項１に記載の伝播時間演算手段に対応している。
【００７６】
上記の伝播時間演算処理が、音波発生対象のノード１２以外の各ノード１２（測定対象のノード１２）で各々行われることにより、音波発生対象のノード１２とそれ以外の各ノード１２との間の測定対象区間（例えば図６（Ａ）において、単一の被測定空間の図に重ねて示した放射状の線に相当する測定対象区間）について、音波発生対象のノード１２から測定対象のノード１２への音波の伝播時間が各々演算され、演算された伝播時間がホストコンピュータ１４へ転送されることになる。
【００７７】
一方、温度分布演算処理（図７）では、音波発生対象のノード１２に対して音波の発生を指示（ステップ１１０）するとステップ１１２へ移行し、全ての測定対象区間の伝播時間の演算が完了したか否か判定し、測定対象のノード１２から音波の伝播時間を受信する毎に、受信した音波の伝播時間をＲＡＭ１４Ｃ等に記憶する処理を繰り返しながら、判定が肯定される迄待機する。
【００７８】
測定対象の全てのノード１２から音波の伝播時間を受信すると、ステップ１１２の判定が肯定されてステップ１１４へ移行し、カウンタｍの値がノード１２の総数ｍMAXに一致したか否か、すなわち全てのノード１２を音波発生対象のノード１２として音波検出処理（音波発生対象のノード１２から音波を発生させ、音波発生対象のノード１２から測定対象のノード１２への音波の伝播時間を演算する処理）を行ったか否かを判定する。判定が否定された場合には、ステップ１１６でカウンタｍの値を１だけインクリメントした後にステップ１０４に戻る。
【００７９】
これにより、ステップ１１４の判定が肯定される迄の間は、音波発生対象のノード１２を順に切り換えてステップ１０４〜１１４が繰り返され（所謂ファンビーム測定：図６（Ａ）参照）、全ての測定対象区間の音波の伝播時間が各々両方向について測定・演算されることになる。
【００８０】
ステップ１１４の判定が肯定されるとステップ１１８へ移行し、ステップ１１８以降で被測定空間２４内の温度分布を演算する。すなわち、ステップ１１８では、所定方向に平行な多数の仮想経路（傾きθが等しい全ての仮想経路）に対応する全ての測定対象区間（図６（Ｂ）において、単一の測定対象区間に重ねて示した互いに平行な線に相当する区間）の伝搬時間のデータを取り込む。
【００８１】
ところで、音波が伝播する空間内に空気流（風）が生じていた場合、音波の伝播時間は空気流の影響（詳しくは、音波の伝播方向に平行な方向の空気流成分の流速の影響）を受けて変化する。本実施形態では、被測定空間２４内に生じている空気流の影響で温度分布の測定精度が低下することを回避するために、全ての測定対象区間の音波の伝播時間を各々両方向について測定しており、次のステップ１２０では、ステップ１１８で取り込んだ伝播時間のデータに基づいて、単一の測定対象区間に対応する両方向の伝播時間のデータ毎にその平均値を演算することを、所定方向に平行な全ての測定対象区間について行う。
【００８２】
例えばノードｍ及びノードｎを各々通る仮想経路に対応する測定対象区間の伝播時間のデータとしては、ノードｍからノードｎへの音波の伝播時間ｔmnのデータと、ノードｎからノードｍへの音波の伝播時間ｔnmのデータが存在しており、これらの平均値ｔm-n（＝（ｔmn＋ｔnm）÷２）を演算する。後述するように、本実施形態ではこの平均値ｔm-nを用いて温度分布を演算するので、空気流の影響による伝播時間の変化分を補正することができる。
【００８３】
なお、本実施形態のようにファンビーム計測を行った場合、一部の測定対象区間については、両方向の音波の伝播時間の測定が時間を空けて行われることになる。被測定空間２４内に生じている空気流が主に空調等に起因するものであれば、空気流の流れる方向や強さは経時的に略一定であることが多く、両方向の伝播時間の測定を時間を空けて行っても影響は小さい。しかし、被測定空間２４内に生じている空気流の流れる方向や強さは経時的に大きく変化する場合には、一方向の伝播時間の測定を行ってから他方向の伝播時間の測定を行う迄の間に空気流の流れる方向や強さが変化することで、空気流の影響による伝播時間の変化に対する補正の精度が不足することも考えられる。
【００８４】
上記のような場合には、ファンビーム計測を行うことに代えて、単一の測定対象区間に対して一方向の伝播時間を測定すると、直ちに他方向の伝播時間を測定することを、全ての測定対象区間について順次行うことが好ましい。これにより、ファンビーム計測と比較して測定に時間はかかるものの、被測定空間２４内に生じている空気流の流れる方向や強さの経時的な変化の影響を受けることなく、被測定空間２４内の温度分布を精度良く求めることができる。
【００８５】
次のステップ１２２では、平均値ｔm-nを演算した全ての測定対象区間について、両方向の音波の伝播時間の差（ｔmn−ｔnm）を各々演算し、演算結果を各測定対象区間における空気流データとしてＲＡＭ１４Ｃ等に記憶する。この空気流データは、温度分布の演算結果を表示する際（後述）に併せて表示するようにしてもよい。
【００８６】
ステップ１２４では、予め設定してＲＯＭ１４Ｂに記憶している伝搬時間推定区間の長さ（仮想スピーカとノード１２（詳しくはノード１２のアーム５４）との間の長さＬSN 、及びノード１２（詳しくはノード１２のアーム５４）と仮想マイクロフォンの間の長さＬNM )を、所定方向に平行な全ての仮想経路について各々取り込む。なお、測定対象区間の距離（各ノード１２のアーム５４間の距離）は光反射板５８を用いてレーザ距離計により正確に計測することができるので、上記の長さＬSN及び長さＬNMについても、測定対象区間の距離の計測結果に基づいて正確に求めることができる。
【００８７】
ステップ１２６では、測定対象区間外が所定温度の媒質（空気）で満たされていると仮定して、所定方向に平行な全ての仮想経路（Ｘ，θ）について、音波の伝搬時間の推定値ｔ（Ｘ，θ）を各々演算する。伝搬時間の推定値ｔ（Ｘ，θ）は次の（１）式により求めることができる。
【００８８】
ｔ（Ｘ，θ）＝（ＬSN÷ｖ0 ）＋ｔm-n＋（ＬNM÷ｖ0 ） …（１）
但し、ｖ0 は空気が所定温度のときの音速であり、ｔm-n、ＬSN及びＬNMは仮想経路（Ｘ，θ）に対応する値である。（１）式における第１項及び第３項は、被測定区間外が一定温度の媒質で満たされていると仮定したときの、仮想経路上の被測定区間外の区間（伝搬時間推定区間）の音波の伝搬時間の推定値に対応している。
【００８９】
ステップ１２８では、所定方向に平行な全ての仮想経路（Ｘ，θ）の音波の伝搬時間の推定値ｔ（Ｘ，θ）に基づいて、投影データｐ（Ｘ，θ）を演算する。
投影データｐ（Ｘ，θ）は次の（２）式によって求めることができる。
【００９０】
ｐ（Ｘ，θ）＝ｔ（Ｘ，θ）−（Ｌ0 ÷ｖ0 ） …（２）
ここで、位置（ｘ，ｙ）における音速をｖ（ｘ，ｙ）とし、音速分布ｆ（ｘ，ｙ）を（温度分布は音速分布ｆ（ｘ，ｙ）から容易に算出できる）、
ｆ（ｘ，ｙ）＝（１／ｖ（ｘ，ｙ))−（１／ｖ0 ） …（３）
上記（３）式のように定義（すなわち、測定対象区間外ではｆ（ｘ，ｙ）＝０となる）すると共に、直線Ｓ（図５参照）に沿う方向をＸ軸、仮想経路に沿う方向をＹ軸とする直交座標系Ｘ−Ｙ（直交座標系ｘ−ｙを原点Ｏ周りに角度θだけ回転させた座標系）を定めると、投影データｐ（Ｘ，θ）と音速分布ｆ（ｘ，ｙ）との関係は、仮想経路（Ｘ，θ）が積分区間となり、
【００９１】
【数１】
【００９２】
となる。（４）式に示したように、投影データｐ（Ｘ，θ）は、仮想経路（Ｘ，θ）に沿った音速分布ｆ（ｘ，ｙ）を投影したものである。前述のステップ１１０では、（２）式の演算を所定方向に平行な全ての仮想経路（Ｘ，θ）（すなわち、距離Ｘが異なり傾きθが等しい全ての仮想経路）について各々行うので、角度θが同一の投影データｐ（Ｘ，θ）の集合（平行な投影データ：parallel projection data）が得られることになる。
【００９３】
次のステップ１３０では全ての方向について投影データｐ（Ｘ，θ）を演算したか否か判定する。判定が否定された場合にはステップ１１８に戻り、ステップ１３０の判定が肯定される迄、ステップ１１８〜ステップ１３０を繰り返す。全ての方向について仮想経路（Ｘ，θ）の投影データｐ（Ｘ，θ）の演算を完了すると（なお全ての方向について仮想経路（Ｘ，θ）の投影データｐ（Ｘ，θ）を演算することは、一般にラドン（Radon)変換と称されている）、これらの投影データ（Ｘ，θ）に基づき、フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等の短時間で演算可能な演算方法を適用することで、被測定空間２４内の温度分布を演算することが可能となる。
【００９４】
このため、ステップ１３０の判定が肯定されるとステップ１３２へ移行し、先のステップ１２８において全ての方向について演算した投影データｐ（Ｘ，θ）に基づいて、被測定空間２４内の温度分布を演算する。このステップ１３２は環境状態演算手段に対応しており、以下では、投影データに基づいて温度分布を演算する演算方法の一例として、フーリエ領域再構成法について説明する。
【００９５】
音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換をＦ（ｕ，ｖ）とすると、
【００９６】
【数２】
【００９７】
フーリエ変換は上記（５）式で定義される。また、Ｆ（ｕ，ｖ）からの２次元フーリエ逆変換は、
【００９８】
【数３】
【００９９】
上記（６）式で定義される。
【０１００】
一方、直交座標系ｘ−ｙと直交座標系Ｘ−Ｙとの間には、以下の（７）式の関係が成り立つ。
【０１０１】
ｘ＝Ｘｃｏｓθ−Ｙｓｉｎθ
ｙ＝Ｘｓｉｎθ＋Ｙｃｏｓθ …（７）
（７）式より、音速分布ｆ（ｘ，ｙ）は次の（８）式で表すことができる。
【０１０２】
なお、ｆ'(Ｘ，Ｙ）はｆ（ｘ，ｙ）を原点Ｏ周りに角度−θだけ回転させることで得られる関数である。（８）式を用いて（５）式のフーリエ変換の定義式を変数変換すると、
【０１０３】
【数４】
【０１０４】
となる。
【０１０５】
また、空間周波数領域においても、実領域と同様に、座標系（ｕ，ｖ）に対して原点周りに角度θだけ回転させることで得られる座標系を座標系（Ｕ，Ｖ）とすると、先の（７）式と同様に、
ｕ＝Ｕｃｏｓθ−Ｖｓｉｎθ
ｖ＝Ｕｓｉｎθ＋Ｖｃｏｓθ …（10）
また
Ｕ＝ｕｃｏｓθ＋Ｖｓｉｎθ
Ｖ＝−ｕｓｉｎθ＋ｖｃｏｓθ …（11）
の関係があり、（10）式よりＦ（ｕ，ｖ）は、
と表すことができる。これらの式を（９）式に代入することにより、次の（13）式を得ることができる。
【０１０６】
【数５】
【０１０７】
（13）式より、実領域において回転した物体（この場合は音速分布ｆ（ｘ，ｙ))のスペクトルは、元の物体のスペクトルを同じ角度だけ回転したものと同じになることがわかる。
【０１０８】
（７）式及び（８）式によれば、投影データｐ（Ｘ，θ）に関する（４）式は次の（14）式のように表すことができる。
【０１０９】
【数６】
【０１１０】
なお、（14）式では音速分布ｆ（ｘ，ｙ）、すなわちｆ'(ｘ，ｙ）が測定対象区間外では０であるため、積分区間を±∞に拡大している。また、投影データｐ（Ｘ，θ）のＸ軸に関するフーリエ変換は、次の（15）式で表すことができる。
【０１１１】
【数７】
【０１１２】
（14）式を（15）式に代入すると、次の（16）式が得られる。
【０１１３】
【数８】
【０１１４】
（16）式は先の（13）式においてＶ＝０としたときの式に一致するので、
Ｐ（Ｕ，θ）＝Ｆ'(Ｕ，０）
が成り立つ。
【０１１５】
このように、傾きθが同一の全ての仮想経路の投影データｐ（Ｘ，θ）に対して、（15）式に従ってＸ軸に関するフーリエ変換（１次元フーリエ変換）を行った結果Ｐ（Ｕ，θ）は、音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換Ｆ（ｕ，ｖ）のθ方向成分（フーリエ変換Ｆ（ｕ，ｖ）のｕ軸に対して角度θだけ傾いた直線（すなわちＵ軸）上の値）に等しい（これを２次元断層定理（projection slice theorem）という）。
【０１１６】
従って、各方向について各々求めた投影データに対して各々フーリエ変換を行って座標系（ｕ，ｖ）上に放射状に並べれば、音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換Ｆ（ｕ，ｖ）を得ることができる。そして、このＦ（ｕ，ｖ）に対し、先の（６）式に従って２次元フーリエ逆変換を行うことにより、音速分布ｆ（ｘ，ｙ）を求めることができる。
【０１１７】
ステップ１３２では、まず上記で説明したフーリエ領域再構成法、或いは公知のフィルタ処理データの逆投影法や２次元フィルタリング法等の短時間で演算可能な演算方法を適用して投影データｐ（Ｘ，θ）から音速分布ｆ（ｘ，ｙ）を求める。次に、求めた音速分布ｆ（ｘ，ｙ）から（３）式に基づいて被測定空間２４内の各位置における音速ｖ（ｘ，ｙ）を演算し、音速ｖと音波が伝搬する媒体（空気）の温度Ｔとの関係を表す次の（17）式に基づき、被測定空間２４内における温度分布Ｔ（ｘ，ｙ）を演算する。
【０１１８】
ｖ＝ 331.45 ＋ 0.607・Ｔ 〔ｍ／秒〕 …（17）
次のステップ１３４では、上記により得られた被測定空間２４内の温度分布を、例として図１０に示すコンターマップ等の形態でディスプレイ２０に表示し、処理を終了する。なお、ディスプレイ２０への表示に代えて、被空調空間２４内の空気調和を行う空調装置に対し、被測定空間２４内の温度分布及び湿度分布を表すデータを出力するようにしてもよい。これにより、空調装置が被測定空間２４に対し、例えば快適性と省エネルギーの双方を満足する空気調和を行うことも可能となる。
【０１１９】
本実施形態では、被測定空間２４の周縁部に多数のノード１２を設け、仮想経路上の一対のノード１２の間の測定対象区間（一対のノード１２の各々のアーム５４の間の区間）については、実際に音波を両方向について各々伝搬させて伝搬時間ｔmn及び伝搬時間ｔnmを測定し、伝播時間の平均値ｔm-nを求め、測定対象区間外の区間（仮想スピーカとノード１２の間、及びノード１２と仮想マイクロフォンの間の伝搬時間推定区間）については、気温が所定温度であると仮定し、空気が所定温度のときの音速ｖ0 を用いて音波の伝搬時間を推定し、伝搬時間の推定値ｔ（Ｘ，θ）を演算しているので、被測定空間２４外にノード１２を多数設けることなく投影データｐ（Ｘ，θ）を得ることができ、この投影データｐ（Ｘ，θ）に基づき、フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等の短時間で演算可能な演算方法を適用して被測定空間内の環境状態を演算することができる。従って、環境状態測定装置１０の構成を簡単にすることができると共に、温度分布の演算に要する時間を短縮することができる。
【０１２０】
なお、上記では第１の位置に到達した音波を第１マイクロフォン３２によって検出する場合を説明したが、音波発生対象のノード１２のスピーカ３０から音波を発生させるための音波データを予め固定的に定めておくと共に、該音波データを用いてスピーカ３０から音波を発生させ、第１の位置に到達した音波をマイクロフォンによって検出し、比較用音波測定データを求めてホストコンピュータ１４のＲＯＭ１４Ｂ、又は各ノード１２のデータ処理部４０の内蔵メモリに記憶しておき、伝播時間の演算時に比較用音波データを読み出して用いることも可能である。
【０１２１】
また、上記では第１マイクロフォン３２及び第２マイクロフォン３４として、同一機種のマイクロフォン（特性が同一のマイクロフォン）を用いていたが、これに限られるものではなく、類似種のマイクロフォン（例えば音波の検出原理が同一等のように特性が類似しているマイクロフォン）であれば第１マイクロフォン３２及び第２マイクロフォン３４としての使用に好適である。
【０１２２】
また、上記では音波の伝搬時間の投影データから温度分布を演算していたが、これに限定されるものではなく、例えば壁体２２の壁面近傍の温度を温度センサによって測定し、壁面近傍の温度の測定値に基づいて前記温度分布を補正するようにしてもよい。
【０１２３】
更に、上記では被測定空間２４を貫くように仮想経路を設定し、仮想経路のうち測定対象区間外の区間については音波の伝播時間を推定することにより、仮想経路上を伝播する音波の伝播時間を求めて温度分布を求めていたが、本発明はこれに限定されるものではなく、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間の音波の伝搬時間と、前記各領域を音波が横切る際の伝搬速度との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の伝搬時間を連立方程式に代入し、伝搬速度が未知の領域について音波が横切る際の伝搬速度を演算し、伝搬速度から各領域の温度を演算することによって温度分布を求めるようにしてもよい。
【０１２４】
また、上記ではドーム球場等のように床面が円形の被測定空間２４内の環境状態を測定する場合を例に説明したが、本発明によれば、音波発生器及び音波検出手段を被測定空間の周縁部に設けることで被測定空間内の温度分布を求めることができるので、任意の形状の空間を被測定空間として採用し、任意の形状の空間内の温度分布を求めることができる。
【０１２５】
また、上記では環境状態として温度分布のみを演算する場合を例に説明したが、これに限定されるものではなく、複数種の周波数の音波について減衰量を各々測定して減衰率分布を求めることにより、被測定空間内の湿度分布も併せて演算したり、他の環境状態の演算も行うようにしてもよい。
【０１２６】
【発明の効果】
請求項１及び請求項４記載の発明は、第１音波検出手段、第２音波検出手段及び光反射板が、音波発生器と連結されたブラケットに一体に取付けられた音波発生検出装置を、被測定空間を横切る複数の伝播時間測定区間の両端部に、ブラケットが取付けられた側同士が対向するように各々配置し、一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、一方の音波発生検出装置の第１音波検出手段、及び他方の音波発生検出装置の第２音波検出手段によって各々検出させることを、複数の伝播時間測定区間について各々行い、検出された音波の波形を照合し、波形の相関性を時間軸に沿って演算した結果に基づいて、第１音波検出手段による音波の検出位置から第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、複数の伝播時間測定区間について各々行い、複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算するので、被測定空間内の温度分布等の環境状態を高精度に測定できる、という優れた効果を有する。
【０１２８】
請求項２記載の発明は、請求項１記載の発明において、第１音波検出手段による音波検出位置の音波到来側に、到来した音波を減衰させる減衰手段を配置するので、上記効果に加え、第１音波検出手段によって検出される音波の波形が歪んだり飽和することを防止することができる、という効果を有する。
【０１２９】
請求項３記載の発明は、請求項１記載の発明において、第１音波検出手段及び第２音波検出手段が少なくとも類似種の音波検出器によって音波を検出するので、音波伝播時間をより精度良く測定することができる、という効果を有する。
【０１３１】
請求項５記載の発明は、請求項４記載の発明において、伝播時間測定区間の両端部の一対の音波発生検出装置の一対の光反射板に各々光を照射して距離を測定することを、複数の伝播時間測定区間について各々定期的に行い、被測定空間内の環境状態の演算を、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して行うので、上記効果に加え、環境温度の変化や経時変化に伴う伝播時間測定区間の距離変化の影響を受けることなく、被測定空間内の環境状態を精度良く測定することができる、という効果を有する。
【図面の簡単な説明】
【図１】 本実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図２】 ノードの配置を示す被測定空間の平面図である。
【図３】 （Ａ）単一のノードにおけるスピーカ及びマイクロフォンの配置を示す側面図、（Ｂ）は（Ａ）におけるマイクロフォン配置位置周辺の拡大図である。
【図４】 測定対象区間及び伝搬時間推定区間から成る仮想経路を、被測定空間の平面図に重ねて示す概念図である。
【図５】 被測定空間に対して設定した直交座標系ｘ−ｙ、各仮想経路を識別するための直線Ｓの距離Ｘ及び傾きθを示す概念図である。
【図６】 （Ａ）及び（Ｂ）は、伝搬時間データの測定順序と、伝搬時間データの処理順序と、の関係を示す概念図である。
【図７】 ホストコンピュータで実行される温度分布演算処理の内容を示すフローチャートである。
【図８】 測定対象の各ノードで実行される音波発生処理及び伝播時間演算処理の内容を示すフローチャートである。
【図９】 （Ａ）は音波データが表す音波の波形の一例、（Ｂ）は比較用音波測定データが表す音波の波形の一例、（Ｃ）は音波測定データが表す音波の波形の一例、（Ｄ）は相関性演算結果の一例を各々示す線図である。
【図１０】 被測定空間内の温度分布の表示例を示すイメージ図である。
【符号の説明】
１０ 環境状態測定装置
１２ ノード
１４ ホストコンピュータ
２４ 被測定空間
３０ スピーカ
３２ 第１マイクロフォン
３４ 第１マイクロフォン
【発明の属する技術分野】
本発明は音波発生検出装置、環境状態測定装置及び方法に係り、特に、環境状態の測定に利用可能な音波発生検出装置、前記音波発生検出装置を用いて空間内の温度分布等の環境状態を測定する環境状態測定方法、及び該環境状態測定方法を適用可能な環境状態測定装置に関する。
【０００２】
【従来の技術】
被空調空間内の温度や湿度を所望の値に制御するためには、被空調空間内の温度、湿度を測定する必要があるが、被空調空間内の各箇所における温度、湿度が一定であることは稀であり、特にドーム球場や劇場等の閉鎖された大空間では、空間内の各箇所における温度、湿度は大きくばらついていることが多い。このため、例えばドーム球場や劇場等の大空間において快適性と省エネルギーの双方を満足する空気調和を行う場合や、室内の温度や湿度を精密に所望の値に制御する場合には、被空調空間内の温度や湿度の分布を測定する必要がある。
【０００３】
しかし、温度センサや湿度センサによって温度や湿度の分布を直接測定しようとすると、センサを、被空調空間の中央部付近を含む被空調空間内にマトリクス状に多数配設する必要があるが、コストが嵩み景観が損なわれると共に被測定空間の利用も制限されるので、上記のように多数のセンサを配設することは現実的ではない。またドーム球場や劇場等では、被測定空間の周縁部（例えば壁体等）に取付けたセンサにより周縁部の温度等を測定すると共に、風船等でセンサを吊り上げることで被測定空間の中央部付近の温度等を測定することが考えられるが、この方法では、被測定空間内で例えば野球等のイベントが実施されている最中に温度や湿度の分布を測定することは不可能である。
【０００４】
一方、医学の分野では、人体に対して多方向からＸ線を照射し、人体を透過したＸ線量を高感度のセンサで測定し、測定されたＸ線量に基づき人体の横断面に沿った各位置における体内組織のＸ線吸収量をコンピュータにより演算して２次元画像として再構成し、人体の横断面像としてモニタに表示するコンピュータ断層撮影（ＣＴ：Computed Tomography)技術が広く利用されている。このＣＴ技術によれば、測定対象（この場合は人体）内部の状態を非接触で測定することができるので、このＣＴ技術を利用して被測定空間内の温度分布を求めることが提案されている。
【０００５】
具体的には、Ｘ線に代えて、媒質の温度によって伝搬時間（速度）が変化する音波を用い、例えば被測定空間を区画している壁体の多数箇所にスピーカ、マイクロフォン、温度センサを設け、スピーカ及びマイクロフォンによって前記多数箇所の間の音波の伝搬時間を各々測定し、被測定空間内の多数箇所の間の音波の伝搬時間から、被測定空間内の温度分布を演算によって求めている。
【０００６】
【発明が解決しようとする課題】
音波の伝播時間の測定は、詳しくはデジタルの原音声データをＤ／Ａ変換器を介してスピーカに供給することでスピーカから音波を発生させ、被測定空間を伝播した音波がマイクロフォン配設位置に到達することでマイクロフォンから出力された音声信号を、Ａ／Ｄ変換器によってデジタルの検出音声データへ変換し、検出音声データの波形と原音声データの波形を時間的にずらしながら相関性を演算し、相関性が最も高くなるときの時間差を音波の伝播時間とすることによって行われる。
【０００７】
しかしながら、スピーカは音声信号に応じて振動板を振動させることで音を発する構造であり、マイクロフォンは周囲の音に応じて発生する振動板の振動を音声信号に変換する構造であるので、振動板の質量の影響により、スピーカから発せられる音の波形及びマイクロフォンから出力される音声信号の波形は、何れも元の波形に対して鈍っていると共に、周波数が高くなるに従って振動板の振動が追従し切れなくなることで所定周波数以上の高周波成分がカットされる。また、スピーカとマイクロフォンは扱う音の音圧レベルが異なっており、より音圧レベルの高い音を扱うスピーカの振動板は、相応の強度が必要とされるために質量が大きく、上記の現象がより顕著である。
【０００８】
このため、スピーカから発せられた音をマイクロフォンで検出することによって得られる音声信号の波形は、スピーカ及びマイクロフォン（特にスピーカ）の影響により原音声データの波形と相違するので、伝播時間の測定における測定精度低下の原因となっており、温度分布の測定精度の低下も招いていた。また、音声信号に対して振動板の振動のレスポンスの悪いスピーカは、マイクロフォンから出力される音声信号の波形の鈍りを助長し、伝播時間及び温度分布の測定精度の更なる低下に繋がるために使用できないという問題もあった。
【０００９】
また、スピーカの振動板は例えばコーン型やドーム型等の立体的な形状である（平面型は波形再現性が劣るという欠点がある）と共に、振動板上における実際に音が発せられる点（発音点）の位置は、振動板の分割振動により、スピーカから発せられる音の周波数によって相違するため、スピーカの発音点とマイクロフォンとの距離もスピーカから発せられる音の周波数によって若干変化する。このスピーカの発音点とマイクロフォンとの距離の変化は、特に狭い空間内における伝播時間及び温度分布の測定に無視できない影響を及ぼすので、伝播時間及び温度分布の測定のためにスピーカから発する音が、例えば単一周波数の単純音又はそれに近い音に制限されていた。
【００１０】
また、３台のスピーカ及びマイクロフォンを用いて発音点の位置を計測する３スピーカ法等のキャリブレーション方法も提案されているが、計測精度が充分でないという問題があった。
【００１２】
本発明は、被測定空間内の温度分布等の環境状態を高精度に測定できる環境状態測定方法及び環境状態測定装置を得ることが目的である。
【００１３】
【課題を解決するための手段】
上記目的を達成するために請求項１記載の発明に係る環境状態測定装置は、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられ、被測定空間を横切る複数の伝播時間測定区間の両端部に、前記ブラケットが取付けられた側同士が対向するように各々配置された音波発生検出装置と、伝播時間測定区間の両端部に配置された一対の前記音波発生検出装置のうちの一方の前記音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の前記音波発生検出装置の前記第１音波検出手段、及び他方の前記音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行う制御手段と、前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行う伝播時間演算手段と、前記伝播時間測定手段によって測定された前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算する環境状態演算手段と、を含んで構成されている。
【００１４】
請求項１記載の発明は、音波発生検出装置を備えており、この音波発生検出装置は、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含んで構成されている。
【００１５】
上記の音波発生器としては、例えばスピーカを用いることができ、詳しくは動電形（例えばダイナミックスピーカ）、電磁形（例えばマグネティックスピーカ）、静電形（例えばコンデンサスピーカ）等の各種スピーカを適用することができるが、音波発生器が発生する音波は可聴域外の周波数の音波であってもよく、例えば音波発生器として超音波を発する送波器等を用いることも可能である。
【００１６】
第１音波検出手段及び第２音波検出手段としては、例えばマイクロフォンを用いることができ、詳しくは動電形（例えばダイナミックマイクロフォンやリボンマイクロフォン）、静電形（例えばコンデンサマイクロフォン）、圧電形（例えばクリスタルマイクロフォン）等の各種マイクロフォン等の音波検出器を適用することができるが、音波発生器が発生する音波が超音波であれば、例えば超音波を検出する受波器等の音波検出器を用いることも可能である。
【００１７】
また、第１音波検出手段と第２音波検出手段は別体の音波検出器である必要はない。すなわち、第１音波検出手段と第２音波検出手段は検出対象の音波の到来方向がおおよそ逆方向であるので、例えば第１音波検出手段及び第２音波検出手段を、両指向性の単一の音波検出器によって構成することも可能である。更に、第１音波検出手段と第２音波検出手段を各々独立した単一指向性の音波検出器によって構成する場合には、第１音波検出手段に対応する音波検出器は音波検出方向が音波発生器側を向くように配置し、第２音波検出手段に対応する音波検出器は音波検出方向が音波発生器と反対側を向くように配置すればよい。
【００１８】
また、上記の音波発生検出装置は、第１音波検出手段、第２音波検出手段及び光反射板が、音波発生器と連結されたブラケットに一体に取付けられ、被測定空間を横切る複数の伝播時間測定区間の両端部に、ブラケットが取付けられた側同士が対向するように各々配置されている。これにより、伝播時間測定区間の両端部に各々配置された一対の音波発生検出装置のうち、一方の音波発生検出装置（便宜的に「第１の音波発生検出装置」という）の音波発生器から音波を発生させれば、該発生させた音波を、第１の音波発生検出装置の第１音波検出手段により、前記区間の前記音波発生器側の端部の所定位置（詳しくは前記音波発生器から所定距離隔てた位置：便宜的に「第１の位置」という）で検出すると共に、他方の音波発生検出装置（便宜的に「第２の音波発生検出装置」という）の第２音波検出手段により、前記区間を挟んで前記音波発生器と反対側の端部の所定位置（便宜的に「第２の位置」という）で検出することができる。
【００１９】
ここで、例えば音波発生器が、入力された音声信号に応じた音波を発生する構成であると共に、発生される音波の波形が音声信号の波形から変化する（例えば波形が鈍る）特性を有しているとしても、音波発生器によって発生された音波を検出する第１の音波発生検出装置の第１音波検出手段及び第２の音波発生検出装置の第２音波検出手段は、音波発生器から発せられた音波（音波発生器によって波形が変化した後の音波）を検出するので、第１音波検出手段によって検出される音波の波形と第２音波検出手段によって検出される音波の波形は高い精度で一致する。
【００２０】
また、音波発生器が音波を発生する際の発音点の位置は、音波発生器が発生する音波の周波数に応じて変化するので、この発音点の位置の変化が第２の音波発生検出装置の第２音波検出手段による音波の検出に影響を及ぼし、第１の音波発生検出装置の音波発生器から発生されて第２の音波発生検出装置の第２音波検出手段で検出される音波の伝播距離は、発音点の位置の変化に応じて変化し、この音波の伝播距離の変化に伴って音波の伝播時間も変化する。
【００２１】
しかしながら、第１の音波発生検出装置の第１音波検出手段は第２の音波発生検出装置の第２音波検出手段と同一の音波を検出するので、上述した発音点の位置の移動は、第１の音波発生検出装置の第１音波検出手段による音波の検出に対しても同様の影響を及ぼし、第１の音波発生検出装置の音波発生器から発生されて第１の音波発生検出装置の第１音波検出手段で検出される音波の伝播距離及び伝播時間は、第２の音波発生検出装置の第２音波検出手段による音波の検出と同様に変化する。
【００２２】
従って、第１の音波発生検出装置の第１音波検出手段によって検出された音波の波形と、第２の音波発生検出装置の第２音波検出手段によって検出された音波の波形を照合し（例えば相関性を演算する等）、伝播時間を測定すべき区間の音波の伝播時間として、音波発生器によって発生された音波が第１の音波発生検出装置の第１音波検出手段に到達してから、該音波が第２の音波発生検出装置の第２音波検出手段に到達する迄の時間差、すなわち第１の音波発生検出装置の第１音波検出手段と、第２の音波発生検出装置の第２音波検出手段の間の音波の伝播時間を求めることにより、音波の伝播時間を、音波発生器の特性や発音点の位置の移動の影響を受けることなく高精度に測定することができる。
【００２３】
また、複数の区間の音波の伝播時間を測定することで温度分布等の環境状態を測定するに際しては、伝播時間測定区間の距離を正確に把握することが重要であるが、伝播時間測定区間として音波発生器から音波検出器に至る区間を用いた場合、発音点の位置が変化すると音波の実際の伝播距離が変化するため、伝播時間測定区間の距離を正確に把握することは困難である。これに対し、請求項１記載の発明に係る音波発生検出装置を用いた音波の伝播時間の測定では、第１の音波発生検出装置の第１音波検出手段と第２の音波発生検出装置の第２音波検出手段の間の区間が伝播時間測定区間となるため、発音点の位置が変化しても音波の実際の伝播距離が変化することはない。
【００２４】
請求項１記載の発明において、制御手段は、伝播時間測定区間の両端部に配置された一対の音波発生検出装置のうちの一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、一方の音波発生検出装置の第１音波検出手段、及び他方の音波発生検出装置の第２音波検出手段によって各々検出させることを、複数の伝播時間測定区間について各々行う。
【００２５】
また伝播時間演算手段は、第２音波検出手段によって検出された音波の波形を第１音波検出手段によって検出された音波の波形と照合し、波形の相関性を時間軸に沿って演算した結果に基づいて、第１音波検出手段による音波の検出位置から第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、複数の伝播時間測定区間について各々行う。これにより、音波発生器の特性の影響や発音点の位置が移動することによる影響を受けることなく、複数の伝播時間測定区間について、音波の伝播時間を各々精度良く測定することができる。
【００２６】
また、上記のように音波の伝播時間の測定精度が発音点位置の移動の影響を受けることがないため、音波発生器から発生させる音波として任意の音波（例えば自然音に相当する音波）を用いることが可能となる。これにより、被測定空間内に人間が存在している場合にも、伝播時間の測定にあたり被測定空間内に存在している人間に違和感や不快感を与えることを回避することができる。
【００２７】
そして環境状態演算手段は、伝播時間測定手段によって測定された複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて、被測定空間内の環境状態を演算するので、被測定空間内の温度分布等の環境状態を高精度に測定することができる。
【００２８】
また、音波発生検出装置は、音波発生器と連結されたブラケットに第１音波検出手段、第２音波検出手段及び光反射板が一体に取付けられているので、音波発生検出装置の設置に際し、音波発生器から距離が等しくなるように第１音波検出手段、第２音波検出手段及び光反射板の相対位置を調整する必要がなくなり、音波発生検出装置の設置を容易に行うことができると共に、周囲温度の変化に伴う第１音波検出手段、第２音波検出手段及び光反射板の相対位置の変化を抑制することができる。
【００２９】
更に、音波発生検出装置のブラケットには、音波発生器側と反対側より到来した光を、第１音波検出手段及び第２音波検出手段に対して音波発生器からの距離が等しい位置（音波発生器から所定距離隔てた位置）で反射する光反射板が設けられている。この光反射板は、例えば光反射面が音波発生器と反対側を向くように配置した全反射ミラーによって構成することができる。この光反射板を利用することで、伝播時間測定区間の距離を、例えばレーザ距離計等によって容易かつ正確に測定することができ、音波の伝播時間の測定結果及び伝播時間測定区間の距離の測定結果に基づいて、温度分布等の環境状態を精度良く測定することができる。
【００３０】
なお、環境状態演算手段による被測定空間内の環境状態の演算は、詳しくは、例えば被測定空間を多数の小領域に区分し、各小領域内の音速を変数とし、或る音波発生装置から或る音波検出装置に至る伝搬経路を伝搬する音波の伝搬時間を、前記伝搬経路上に存在する各小領域の音波の通過時間の和で表した方程式を連立方程式として多数生成し、測定した音波の伝搬時間を連立方程式に代入して解くことによっても実現できるが、下記のようにして演算することが好ましい。
【００３１】
すなわち、各々一定の方向に沿いかつ互いに異なる位置で被測定空間を貫く等しい長さの多数の仮想経路を、互いに異なる複数の方向について各々設定しておき、前記各仮想経路について、第１の音波発生検出装置の第１音波検出手段と、第２の音波発生検出装置の第２音波検出手段と、の間に相当する仮想経路上の測定区間の音波の伝搬時間を実際に音波を伝搬させて測定し、仮想経路上の測定区間外の区間の音波の伝搬時間を被測定空間外が一定温度の媒質で満たされていると仮定して推定することにより、仮想経路上を音波が伝搬するときの伝搬時間の推定値を各々求め、所定方向に沿った多数の仮想経路の各々の音波の伝搬時間の推定値に基づき、前記所定方向についての音波の伝搬時間の投影データを求めることを、前記複数の方向について各々行い、前記複数の方向の各々についての音波の伝搬時間の投影データに基づいて被測定空間内の環境状態を演算することが好ましい。
【００３２】
上記では、伝搬時間測定経路の長さを仮想的に一致させることで、音波の伝搬時間の投影データを取得できるので、短時間で演算可能な演算方法（フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等）を適用することができ、演算時間の短縮や伝播時間を測定するための装置構成の簡易化を実現できる。
【００３３】
なお、請求項１記載の発明に係る第１音波検出手段は、検出対象の音波を発生する音波発生器が比較的近い位置に位置しているので、検出対象の音波の音圧レベルが比較的高く、第１音波検出手段によって検出される音波の波形が歪んだり飽和したりする可能性がある。これを考慮すると、請求項２に記載したように、第１音波検出手段による音波検出位置の音波到来側に、到来した音波を減衰させる減衰手段を配置することが好ましい。なお、減衰手段としては、例えば吸音材等を用いることができる。これにより、第１音波検出手段によって検出される音波の波形が歪んだり飽和することを防止することができる。
【００３４】
また、音波検出手段によって検出される音波の波形は、到来した音波の波形に対して若干ではあるが変化する（例えば鈍りが生ずる等）。これを考慮すると、請求項３に記載したように、第１音波検出手段及び第２音波検出手段を、少なくとも類似種（好ましくは同一種）の音波検出器によって音波を検出するよう構成することが好ましい。これにより、第１音波検出手段によって検出される波形及び第２音波検出手段によって検出される波形が、到来した音波の波形に対して同様の変化を示すことになる。
【００３５】
従って、本発明に係る音波発生検出装置を、音波の伝播時間を測定すべき区間の両端部に各々配置する態様において、音波発生器から音波を発生させる第１の音波発生検出装置の第１音波検出手段によって検出される音波の波形と、第２の音波発生検出装置の第２音波検出手段によって検出される音波の波形を高い精度で一致させることができ、音波伝播時間をより精度良く測定することができる。
【００３８】
請求項４記載の発明に係る環境状態測定方法は、被測定空間を横切る複数の伝播時間測定区間の両端部に、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられた音波発生検出装置を、前記ブラケットが取付けられた側同士が対向するように各々配置しておき、伝播時間測定区間の両端部に配置した一対の音波発生検出装置のうちの一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の音波発生検出装置の前記第１音波検出手段、及び他方の音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行い、前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行い、前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算するので、請求項１記載の発明と同様に、被測定空間内の温度分布等の環境状態を高精度に測定することができる。
【００３９】
また、請求項５記載の発明は、請求項４記載の発明において、伝播時間測定区間の両端部に各々配置した一対の音波発生検出装置の一対の光反射板に各々光を照射して前記一対の光反射板の距離を測定することを、前記複数の伝播時間測定区間について各々行い、前記被測定空間内の環境状態の演算を、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して行うことを特徴としている。
【００４０】
請求項５記載の発明では、伝播時間測定区間の両端部に各々配置した一対の音波発生検出装置の一対の光反射板に各々光を照射して一対の光反射板の距離を測定する。なお、この距離測定には例えばレーザ距離計等を用いることができ、光反射板にレーザ光を照射することで一対の光反射板の距離、すなわち伝播時間測定区間の距離を容易かつ正確に測定することができる。
【００４１】
そして請求項５記載の発明では、上記の距離の測定を複数の伝播時間測定区間について各々定期的に行い、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して、被測定空間内の環境状態の演算を行うので、例えば季節変化等に起因する環境温度の変化や経時変化に伴って伝播時間測定区間の距離が変化したとしても、距離変化の影響を受けることなく、被測定空間内の環境状態を精度良く測定することができる。
【００４２】
【発明の実施の形態】
以下、図面を参照して本発明の実施形態の一例を詳細に説明する。図１には、環境状態測定方法が適用された環境状態測定装置１０が示されている。環境状態測定装置１０は、多数のノード１２Ａ、１２Ｂ、１２Ｃ、…と、ホストコンピュータ１４を備えており、これらが伝送媒体１８を介して互いに接続されて構成されている。
【００４３】
ホストコンピュータ１４は、ＣＰＵ１４Ａ、ＲＯＭ１４Ｂ、ＲＡＭ１４Ｃ、入出力ポート１４Ｄを備え、これらがバス１４Ｅを介して互いに接続されて構成されている。入出力ポート１４Ｄはネットワーク伝送部１６を介して伝送媒体１８に接続されている。また、入出力ポート１４Ｄには、各種の情報を表示するためのディスプレイ２０と、オペレータが各種のデータやコマンド等を入力するためのキーボード２２が接続されている。
【００４４】
一方、多数のノード１２は、被測定空間の周縁部に、周縁部の全周に亘って略一定の間隔で配置されている。例として図２には、円筒状の壁体２２によって区画され、上方が屋根（図示省略）によって閉塞されたドーム状の被測定空間２４（例えばドーム球場等）内の環境状態を測定する場合の各ノード１２の配置が示されている。なお、ノード１２の数及び間隔は図２に示した例に限定されるものではないが、ノード１２の数は偶数であることが望ましい。また、各ノード１２を識別するために、各ノードには便宜的にノード番号が付与されている。
【００４５】
ノード１２Ａ、１２Ｂ、１２Ｃ、…は各々同一の構成であるので、以下ではノード１２Ａについて説明する。ノード１２Ａは、被測定空間２４内に向けて配置された音波発生器としてのスピーカ３０と、スピーカ３０へ向けて配置された第１音波検出手段としての第１マイクロフォン３２と、被測定空間２４内に向けて配置された第２音波検出手段としての第２マイクロフォン３４を備えている。スピーカ３０は壁体２２に埋設されており（図示省略）、増幅器３６、Ｄ／Ａ変換器３８を介してデータ処理部４０に接続されている。
【００４６】
図３（Ａ）に示すように、スピーカ３０には、中間部が被測定空間２４側へ突出するようにＵ字型に屈曲されたアーム５４の基部が取付けられている。図３（Ｂ）に示すように、アーム５４の中間部には、スピーカ３０と対向する面に孔５４Ａが穿設されており、第１マイクロフォン３２は、スピーカ３０から発せられた音波をスピーカ３０から所定距離隔てた位置で検出するように、孔５４Ａの内部に埋設されている。なお、以下では第１マイクロフォン３２の配置位置を「第１の位置」と称する。
【００４７】
第１マイクロフォン３２は自ノードのスピーカ３０から発せられた音波を検出するので、音源（スピーカ３０）との距離が小さく、検出対象の音波の音圧レベルが比較的高い。このため、第１マイクロフォン３２は孔５４Ａの底部側に配置されており、孔５４Ａの手前側（第１マイクロフォン３２とスピーカ３０の間）には吸音材５６が設けられている。これにより、第１マイクロフォン３２から出力される音声信号の信号レベルが歪んだり飽和することが防止される。吸音材５６は請求項２に記載の減衰手段に対応している。
【００４８】
また、図３（Ｂ）に示すように、アーム５４の中間部には、孔５４Ａが穿設された面と反対側の面（被測定空間２４側を向いている面）に孔５４Ｂ、５４Ｃが各々穿設されており、第２マイクロフォン３４は、他ノードのスピーカ３０から発せられた音波を、スピーカ３０から所定距離隔てた位置で検出するように、孔５４Ｂの内部に埋設されている。また、以下では第２マイクロフォン３４の配置位置を「第２の位置」と称する。第２マイクロフォン３４は他ノードのスピーカ３０から発せられた音波を検出するので、音源との距離が大きく検出対象の音波の音圧レベルが比較的小さいので、孔５４Ｂには吸音材は設けられていない。
【００４９】
なお、本実施形態では第１マイクロフォン３２及び第２マイクロフォン３４として同一機種のマイクロフォンを用いており、第１マイクロフォン３２及び第２マイクロフォン３４は請求項３に記載の「類似種の音波検出器」に対応している。一方、孔５４Ｃの底部には、スピーカ３０側と反対側より到来した光をスピーカ３０から所定距離隔てた位置で反射するように光反射板５８が取付けられている。光反射板５８は各ノード１２のアーム５４間の距離をレーザ距離計等によって測定する際に用いられる。
【００５０】
なお、上述したスピーカ３０、アーム５４、第１マイクロフォン３２、第２マイクロフォン３４及び光反射板５８は、請求項１に記載の音波発生検出装置に対応している。
【００５１】
第１マイクロフォン３２はフィルタ／増幅器４２、Ａ／Ｄ変換器４４を介してデータ処理部４０に接続されており、第２マイクロフォン３４はフィルタ／増幅器４６、Ａ／Ｄ変換器４８を介してデータ処理部４０に接続されている。データ処理部４０はマイクロコンピュータ等を含んで構成されており、ネットワーク伝送部５０を介して伝送媒体１８に接続されている。
【００５２】
次に本実施形態の作用を説明する。本実施形態では、後述する温度分布演算処理において、音波の伝搬時間の推定を行う経路として、例として図４に示すように、多数の仮想経路を予め設定している。なお、図４では或る一定の方向に沿った仮想経路のみを示しているが、実際には、一例として図６（Ｂ）にも示すように、互いに異なる複数の方向（例えば１０°ずつ異なる方向）について、各々多数の仮想経路が予め設定されている。
【００５３】
仮想経路は、被測定空間２４外に、所定方向に直交するように配置された仮想スピーカ列上に位置している仮想スピーカと、被測定空間２４を挟んで仮想スピーカ列と反対側の被測定空間２４外に一定の間隔Ｌ0 （図５参照）を隔てて仮想スピーカ列と平行に配置された仮想マイクロフォン列上に位置している仮想マイクロフォンと、を結ぶ経路として定義されており、仮想スピーカ及び仮想マイクロフォンは、各仮想経路が、各々所定方向に沿いかつ互いに異なるノード１２が配置されている位置で壁体２２と交差して被測定空間２４を貫くように、仮想スピーカ列及び仮想マイクロフォン列上での位置が定められている。従って、各仮想経路の長さは互いに等しくされている。
【００５４】
各仮想経路のうち、同一の仮想経路上に位置している一対のノード１２の間の区間（図４に実線で示す測定対象区間：詳しくは一対のノード１２の各々のアーム５４の間の区間）については、実際に音波を伝搬させることによって音波の伝搬時間を測定するが、仮想スピーカとノード１２の間、及びノード１２と仮想マイクロフォンの間の区間（伝搬時間推定区間）については、音波の伝搬時間を測定することはできない。このため、ホストコンピュータ１４のＲＯＭ１４Ｂには、各仮想経路に対応して伝搬時間推定区間の長さが各々記憶されている。
【００５５】
なお、本実施形態では、図５に示すように、被測定空間２４の中心を原点Ｏとする直交座標系ｘ−ｙを定義し、仮想経路の各々を、仮想経路に直交し原点Ｏを通る直線Ｓの長さＸ（距離Ｘ）、及び前記直線Ｓとｘ軸との成す角度θ（傾きθ）によって識別する（仮想経路（Ｘ，θ）と称する）。
【００５６】
次に図７のフローチャートを参照し、被測定空間２４内の環境状態（本実施形態では温度分布）を求める場合にホストコンピュータ１４で実行される温度分布演算処理について説明する。
【００５７】
本実施形態では、伝播時間の測定に用いる音波データ（音波をアナログ信号に変換し更にデジタルデータに変換したもの）がホストコンピュータ１４のＲＯＭ１４Ｂ等に予め記憶されており、ステップ１００では、伝送媒体１８を介して前記音波データを全てのノード１２へ転送する。各ノード１２に転送された音波データは、各ノード１２のデータ処理部４０の内蔵メモリ等に記憶される。
【００５８】
なお、この音波データとしては任意の音声を表す音波データを用いることができ、例えば音楽、メッセージ、自然音（鳥のさえずりやドアの開閉音等）を表す音波データを用いることができる。また、被測定空間２４の大きさが大きくなるに従って（すなわち被測定空間１２の中央部を挟んで対向するノード１２間の距離が大きくなるに従って）、音波の周波数帯域が低周波側に偏倚するように音波データを選択してもよい。
【００５９】
次のステップ１０２では音波を発生させるノード１２のノード番号を表すカウント値ｍに１を代入する。ステップ１０４では音波発生対象のノードｍからの音波発生時刻を決定し、次のステップ１０６では、決定した音波発生時刻及び音波発生対象のノード１２のノード番号を全てのノード１２に通知する。次のステップ１０８では、音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。音波発生タイミングが到来すると、ステップ１０８の判定が肯定されてステップ１１０へ移行し、音波発生対象のノードｍに対し、先に転送した音波データを用いて音波を発生するよう指示する。
【００６０】
音波発生対象のノードｍのデータ処理部４０では、図８（Ａ）に示す音波発生処理が行われる。この音波発生処理は、ホストコンピュータ１４から通知された音波発生対象のノード１２のノード番号に基づき、自ノードが音波発生対象のノード１２であることを認識すると実行され、まずステップ１５０で音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。
【００６１】
伝送媒体１８及びネットワーク伝送部４８を介してホストコンピュータ１４から音波の発生が指示されると、ステップ１５０の判定が肯定されてステップ１５２へ移行し、先に転送された音波データをＤ／Ａ変換器３８に出力する。データ処理部４０から出力された音波データは、Ｄ／Ａ変換器３８でアナログの音声信号に変換され、増幅器３６によって所定の増幅率で増幅された後にスピーカ３０に供給される。これにより、ノードｍのスピーカ３０からは、温度分布演算処理（図７）のステップ１０４で判定された音波発生タイミングで、音波データが表す音波が所定の音量で発せられる。
【００６２】
音波発生対象のノードｍのスピーカ３０によって発生された音波はノードｍの第１マイクロフォン３２によって検出され、ノードｍの第１マイクロフォン３２からアナログの音声信号が出力される。ステップ１５２では、スピーカ３０から音波を発生させる処理と並行してＡ／Ｄ変換器４４を所定時間（例えば数秒）作動させ、第１マイクロフォン３２が音波を検出することで第１マイクロフォン３２から出力される音声信号のサンプリングを所定時間実行させる。
【００６３】
これにより、第１マイクロフォン３２から出力された音声信号は、所定時間の間、フィルタ／増幅器４２で雑音が除去されて増幅された後に、Ａ／Ｄ変換器４４によって一定のサンプリング周期でサンプリングされてデジタルの音波データに変換され、所定時間分の比較用音波測定データ（第１マイクロフォン３２で検出された音波の振幅を、所定時間に亘り（ｔWAVE＝０〜ｔLAST：図９（Ｂ）参照）、Ａ／Ｄ変換器４２のサンプリング周期に相当する時間間隔で時系列に表すデータ）としてデータ処理部４０に順次入力される。データ処理部４０は順次入力される比較用音波測定データをメモリ等に記憶する。
【００６４】
ステップ１５４では、上記の処理によって得られた所定時間分の比較用音波測定データを自ノード以外の各ノード１２へ転送し、音波発生処理を終了する。
【００６５】
例として、図９（Ａ）には音波発生対象のノード１２からの音波の発生に用いた音波データが表す音波の波形の一例を、図９（Ｂ）には比較用音波測定データが表す音波の波形の一例を各々示す。比較用音波測定データは、スピーカ３０から発せられた音波を第１マイクロフォン３２によって検出した結果を表しているので、図９（Ａ）の波形と図９（Ｃ）の波形を比較しても明らかなように、比較用音波測定データが表す音波の波形は、スピーカ３０及び第１マイクロフォン３２の特性の影響を受け、元の音波データが表す音波の波形に対して若干変化している（歪んでいる）。
【００６６】
一方、音波発生対象のノード１２以外のノード１２（測定対象のノード１２）のデータ処理部４０では、図８（Ｂ）に示す伝播時間演算処理が行われる。この伝播時間演算処理は、ホストコンピュータ１４から通知された音波発生対象のノード１２のノード番号に基づき、自ノードが音波発生対象のノード１２でないことを認識すると実行され、まずステップ１６０では、ホストコンピュータ１４からノード番号と共に通知された音波発生時刻を記憶する。次のステップ１６２では、記憶した音波発生時刻を現時刻と比較することで音波発生タイミングが到来したか否か判定し、判定が肯定される迄待機する。音波発生時刻が到来すると、上記判定が肯定されてステップ１６４へ移行する。
【００６７】
先にも説明したように音波発生時刻が到来すると、音波発生対象のノード１２のスピーカから音波が発せられるが、この音波は被測定空間２４全域に放射され、音波測定対象の各ノード１２の第２マイクロフォン３４によって各々検出され、音波測定対象の各ノード１２の第２マイクロフォン３４からアナログの音声信号が出力される。ステップ１６４ではＡ／Ｄ変換器４８を所定時間（例えば数秒）作動させ、第２マイクロフォン３４が音波を検出することで第２マイクロフォン３４から出力される音声信号のサンプリングを所定時間実行させる。
【００６８】
これにより、第２マイクロフォン３２から出力された音声信号は、所定時間の間、フィルタ／増幅器４６で雑音が除去されて増幅された後に、Ａ／Ｄ変換器４８によって一定のサンプリング周期でサンプリングされてデジタルの音波データに変換され、所定時間分の音波測定データ（第２マイクロフォン３２で検出された音波の振幅を、所定時間に亘り、Ａ／Ｄ変換器４２のサンプリング周期に相当する時間間隔で時系列に表すデータ）としてデータ処理部４０に順次入力される。データ処理部４０は順次入力される音波測定データをメモリ等に記憶する。なお、音波測定データのデータ量を削減するために、音波発生時刻が到来してから或る時間（音波発生対象のノード１２からの音波の伝播時間よりも若干短い時間）経過後に音声信号のサンプリングを開始するようにしてもよい。
【００６９】
次のステップ１６６では音波発生対象のノード１２から比較用音波測定データを受信したか否か判定し、判定が肯定される迄待機する。比較用音波測定データを受信すると、ステップ１６６の判定が肯定されてステップ１６８へ移行し、音波測定データが表す波形と音波発生対象のノード１２から受信した比較用音波測定データが表す波形との相関性を演算する。
【００７０】
音波測定データが表す波形と比較用音波測定データが表す波形との相関性は、音波発生対象のノード１２の第１マイクロフォン３２の配設位置（第１の位置）に測定対象音波が到達したタイミングを基準とし、この基準タイミングから、測定対象のノード１２の第２マイクロフォン３４の配設位置（第２の位置）に測定対象音波が到達する迄の時間（音波の伝播時間）として或る値ｔX を仮定し、ｔWAVEが所定値（初期値はｔWAVE＝０：基準タイミング（図９（Ｂ）参照))のときの振幅値を表すデータを比較用音波測定データから取り出すと共に、対応するタイミング（基準タイミングからｔX ＋ｔWAVE経過したタイミング）における振幅値を表すデータを音波測定データから取り出し、両者を乗算することを、ｔWAVEの値を０〜ｔLAST迄順に変化させながら繰り返し、乗算結果を積算することで求めることができる。
【００７１】
上記の演算によって求まる相関性の値は、音波発生対象のノード１２のアーム５４の配設位置（第１マイクロフォン３２の配設位置）に測定対象音波が到達してから時間ｔX が経過したときに、測定対象のノード１２のアーム５４の配設位置（第２マイクロフォン３４の配設位置）に測定対象音波が到達した確からしさを表しており（相関性の値が符号が正でかつ絶対値が大きくなる程、前記確からしさが高いとみなすことができる）、上記の演算処理を、時間ｔX の値を徐々に変化させながら繰り返すことで、例として図９（Ｄ）に示すように、音波測定データが表す波形と比較用音波測定データが表す波形との相関性の値の時間軸に沿った推移を表すデータを得ることができる。
【００７２】
次のステップ１７０では、ステップ１６８における相関性の演算結果に基づいて、相関性の値が最大になるときの時間ｔX の値を求め、音波発生対象のノード１２から測定対象のノード１２（自ノード）への測定対象音波の伝播時間（詳しくは音波発生対象のノード１２の第１の位置から測定対象のノード１２の第２の位置への測定対象音波の伝播時間）として、相関性の値が最大のときの時間ｔX を設定し、設定した伝播時間を、ノードｍからノードｎへの音波の伝播時間データｔmnとしてホストコンピュータ１４へ転送し、伝播時間演算処理を終了する。
【００７３】
例として、図９（Ｃ）には音波測定データが表す音波の波形の一例を示す。音波測定データは、スピーカ３０から発せられた音波を第２マイクロフォン３４によって検出した結果を表しているので、図９（Ａ）の波形と図９（Ｃ）の波形を比較しても明らかなように、音波測定データが表す音波の波形は、スピーカ３０及び第２マイクロフォン３４の特性の影響を受け、元の音波データが表す音波の波形に対して若干変化している（歪んでいる）。
【００７４】
しかしながら、比較用音波測定データは、同一のスピーカ３０から発せられた音波を、第２マイクロフォン３４と同一機種の第１マイクロフォン３２（すなわち、第２マイクロフォン３４と特性が同一のマイクロフォン）によって検出した結果を表しているので、音波測定データが表す音波の波形及び比較用音波測定データが表す音波の波形は、図９（Ｂ）の波形と図９（Ｃ）の波形を比較しても明らかなように、元の音波データが表す音波の波形に対して同様に変化しており、波形の相関度が飛躍的に向上している。
【００７５】
従って、ステップ１６８の演算によって求まる相関性の値は、時間ｔX の値が音波発生対象のノード１２から測定対象のノード１２への測定対象音波の伝播時間に一致したときに正の最大値となり、時間ｔX の値がその前後の値のときにも相関性の値の絶対値は多少大きいものの、時間ｔX の値が前記伝播時間と大きく異なる値のときには相関性の値の絶対値は非常に小さな値となる。これにより、スピーカ３０や第２マイクロフォン３４の特性の影響を受けることなく、伝播時間を高精度に測定・設定することができる。また、各ノード１２のアーム５４の間を伝播速度の測定対象区間としているので、測定対象区間の音波伝播時間の測定結果が、スピーカ３０から発せられる音波の周波数変化に伴う発音点位置の変化の影響を受けて変動することもない。なお、ステップ１６８，１７０は請求項１に記載の伝播時間演算手段に対応している。
【００７６】
上記の伝播時間演算処理が、音波発生対象のノード１２以外の各ノード１２（測定対象のノード１２）で各々行われることにより、音波発生対象のノード１２とそれ以外の各ノード１２との間の測定対象区間（例えば図６（Ａ）において、単一の被測定空間の図に重ねて示した放射状の線に相当する測定対象区間）について、音波発生対象のノード１２から測定対象のノード１２への音波の伝播時間が各々演算され、演算された伝播時間がホストコンピュータ１４へ転送されることになる。
【００７７】
一方、温度分布演算処理（図７）では、音波発生対象のノード１２に対して音波の発生を指示（ステップ１１０）するとステップ１１２へ移行し、全ての測定対象区間の伝播時間の演算が完了したか否か判定し、測定対象のノード１２から音波の伝播時間を受信する毎に、受信した音波の伝播時間をＲＡＭ１４Ｃ等に記憶する処理を繰り返しながら、判定が肯定される迄待機する。
【００７８】
測定対象の全てのノード１２から音波の伝播時間を受信すると、ステップ１１２の判定が肯定されてステップ１１４へ移行し、カウンタｍの値がノード１２の総数ｍMAXに一致したか否か、すなわち全てのノード１２を音波発生対象のノード１２として音波検出処理（音波発生対象のノード１２から音波を発生させ、音波発生対象のノード１２から測定対象のノード１２への音波の伝播時間を演算する処理）を行ったか否かを判定する。判定が否定された場合には、ステップ１１６でカウンタｍの値を１だけインクリメントした後にステップ１０４に戻る。
【００７９】
これにより、ステップ１１４の判定が肯定される迄の間は、音波発生対象のノード１２を順に切り換えてステップ１０４〜１１４が繰り返され（所謂ファンビーム測定：図６（Ａ）参照）、全ての測定対象区間の音波の伝播時間が各々両方向について測定・演算されることになる。
【００８０】
ステップ１１４の判定が肯定されるとステップ１１８へ移行し、ステップ１１８以降で被測定空間２４内の温度分布を演算する。すなわち、ステップ１１８では、所定方向に平行な多数の仮想経路（傾きθが等しい全ての仮想経路）に対応する全ての測定対象区間（図６（Ｂ）において、単一の測定対象区間に重ねて示した互いに平行な線に相当する区間）の伝搬時間のデータを取り込む。
【００８１】
ところで、音波が伝播する空間内に空気流（風）が生じていた場合、音波の伝播時間は空気流の影響（詳しくは、音波の伝播方向に平行な方向の空気流成分の流速の影響）を受けて変化する。本実施形態では、被測定空間２４内に生じている空気流の影響で温度分布の測定精度が低下することを回避するために、全ての測定対象区間の音波の伝播時間を各々両方向について測定しており、次のステップ１２０では、ステップ１１８で取り込んだ伝播時間のデータに基づいて、単一の測定対象区間に対応する両方向の伝播時間のデータ毎にその平均値を演算することを、所定方向に平行な全ての測定対象区間について行う。
【００８２】
例えばノードｍ及びノードｎを各々通る仮想経路に対応する測定対象区間の伝播時間のデータとしては、ノードｍからノードｎへの音波の伝播時間ｔmnのデータと、ノードｎからノードｍへの音波の伝播時間ｔnmのデータが存在しており、これらの平均値ｔm-n（＝（ｔmn＋ｔnm）÷２）を演算する。後述するように、本実施形態ではこの平均値ｔm-nを用いて温度分布を演算するので、空気流の影響による伝播時間の変化分を補正することができる。
【００８３】
なお、本実施形態のようにファンビーム計測を行った場合、一部の測定対象区間については、両方向の音波の伝播時間の測定が時間を空けて行われることになる。被測定空間２４内に生じている空気流が主に空調等に起因するものであれば、空気流の流れる方向や強さは経時的に略一定であることが多く、両方向の伝播時間の測定を時間を空けて行っても影響は小さい。しかし、被測定空間２４内に生じている空気流の流れる方向や強さは経時的に大きく変化する場合には、一方向の伝播時間の測定を行ってから他方向の伝播時間の測定を行う迄の間に空気流の流れる方向や強さが変化することで、空気流の影響による伝播時間の変化に対する補正の精度が不足することも考えられる。
【００８４】
上記のような場合には、ファンビーム計測を行うことに代えて、単一の測定対象区間に対して一方向の伝播時間を測定すると、直ちに他方向の伝播時間を測定することを、全ての測定対象区間について順次行うことが好ましい。これにより、ファンビーム計測と比較して測定に時間はかかるものの、被測定空間２４内に生じている空気流の流れる方向や強さの経時的な変化の影響を受けることなく、被測定空間２４内の温度分布を精度良く求めることができる。
【００８５】
次のステップ１２２では、平均値ｔm-nを演算した全ての測定対象区間について、両方向の音波の伝播時間の差（ｔmn−ｔnm）を各々演算し、演算結果を各測定対象区間における空気流データとしてＲＡＭ１４Ｃ等に記憶する。この空気流データは、温度分布の演算結果を表示する際（後述）に併せて表示するようにしてもよい。
【００８６】
ステップ１２４では、予め設定してＲＯＭ１４Ｂに記憶している伝搬時間推定区間の長さ（仮想スピーカとノード１２（詳しくはノード１２のアーム５４）との間の長さＬSN 、及びノード１２（詳しくはノード１２のアーム５４）と仮想マイクロフォンの間の長さＬNM )を、所定方向に平行な全ての仮想経路について各々取り込む。なお、測定対象区間の距離（各ノード１２のアーム５４間の距離）は光反射板５８を用いてレーザ距離計により正確に計測することができるので、上記の長さＬSN及び長さＬNMについても、測定対象区間の距離の計測結果に基づいて正確に求めることができる。
【００８７】
ステップ１２６では、測定対象区間外が所定温度の媒質（空気）で満たされていると仮定して、所定方向に平行な全ての仮想経路（Ｘ，θ）について、音波の伝搬時間の推定値ｔ（Ｘ，θ）を各々演算する。伝搬時間の推定値ｔ（Ｘ，θ）は次の（１）式により求めることができる。
【００８８】
ｔ（Ｘ，θ）＝（ＬSN÷ｖ0 ）＋ｔm-n＋（ＬNM÷ｖ0 ） …（１）
但し、ｖ0 は空気が所定温度のときの音速であり、ｔm-n、ＬSN及びＬNMは仮想経路（Ｘ，θ）に対応する値である。（１）式における第１項及び第３項は、被測定区間外が一定温度の媒質で満たされていると仮定したときの、仮想経路上の被測定区間外の区間（伝搬時間推定区間）の音波の伝搬時間の推定値に対応している。
【００８９】
ステップ１２８では、所定方向に平行な全ての仮想経路（Ｘ，θ）の音波の伝搬時間の推定値ｔ（Ｘ，θ）に基づいて、投影データｐ（Ｘ，θ）を演算する。
投影データｐ（Ｘ，θ）は次の（２）式によって求めることができる。
【００９０】
ｐ（Ｘ，θ）＝ｔ（Ｘ，θ）−（Ｌ0 ÷ｖ0 ） …（２）
ここで、位置（ｘ，ｙ）における音速をｖ（ｘ，ｙ）とし、音速分布ｆ（ｘ，ｙ）を（温度分布は音速分布ｆ（ｘ，ｙ）から容易に算出できる）、
ｆ（ｘ，ｙ）＝（１／ｖ（ｘ，ｙ))−（１／ｖ0 ） …（３）
上記（３）式のように定義（すなわち、測定対象区間外ではｆ（ｘ，ｙ）＝０となる）すると共に、直線Ｓ（図５参照）に沿う方向をＸ軸、仮想経路に沿う方向をＹ軸とする直交座標系Ｘ−Ｙ（直交座標系ｘ−ｙを原点Ｏ周りに角度θだけ回転させた座標系）を定めると、投影データｐ（Ｘ，θ）と音速分布ｆ（ｘ，ｙ）との関係は、仮想経路（Ｘ，θ）が積分区間となり、
【００９１】
【数１】
となる。（４）式に示したように、投影データｐ（Ｘ，θ）は、仮想経路（Ｘ，θ）に沿った音速分布ｆ（ｘ，ｙ）を投影したものである。前述のステップ１１０では、（２）式の演算を所定方向に平行な全ての仮想経路（Ｘ，θ）（すなわち、距離Ｘが異なり傾きθが等しい全ての仮想経路）について各々行うので、角度θが同一の投影データｐ（Ｘ，θ）の集合（平行な投影データ：parallel projection data）が得られることになる。
【００９３】
次のステップ１３０では全ての方向について投影データｐ（Ｘ，θ）を演算したか否か判定する。判定が否定された場合にはステップ１１８に戻り、ステップ１３０の判定が肯定される迄、ステップ１１８〜ステップ１３０を繰り返す。全ての方向について仮想経路（Ｘ，θ）の投影データｐ（Ｘ，θ）の演算を完了すると（なお全ての方向について仮想経路（Ｘ，θ）の投影データｐ（Ｘ，θ）を演算することは、一般にラドン（Radon)変換と称されている）、これらの投影データ（Ｘ，θ）に基づき、フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等の短時間で演算可能な演算方法を適用することで、被測定空間２４内の温度分布を演算することが可能となる。
【００９４】
このため、ステップ１３０の判定が肯定されるとステップ１３２へ移行し、先のステップ１２８において全ての方向について演算した投影データｐ（Ｘ，θ）に基づいて、被測定空間２４内の温度分布を演算する。このステップ１３２は環境状態演算手段に対応しており、以下では、投影データに基づいて温度分布を演算する演算方法の一例として、フーリエ領域再構成法について説明する。
【００９５】
音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換をＦ（ｕ，ｖ）とすると、
【００９６】
【数２】
フーリエ変換は上記（５）式で定義される。また、Ｆ（ｕ，ｖ）からの２次元フーリエ逆変換は、
【００９８】
【数３】
上記（６）式で定義される。
【０１００】
一方、直交座標系ｘ−ｙと直交座標系Ｘ−Ｙとの間には、以下の（７）式の関係が成り立つ。
【０１０１】
ｘ＝Ｘｃｏｓθ−Ｙｓｉｎθ
ｙ＝Ｘｓｉｎθ＋Ｙｃｏｓθ …（７）
（７）式より、音速分布ｆ（ｘ，ｙ）は次の（８）式で表すことができる。
【０１０２】
【０１０３】
【数４】
となる。
【０１０５】
また、空間周波数領域においても、実領域と同様に、座標系（ｕ，ｖ）に対して原点周りに角度θだけ回転させることで得られる座標系を座標系（Ｕ，Ｖ）とすると、先の（７）式と同様に、
ｕ＝Ｕｃｏｓθ−Ｖｓｉｎθ
ｖ＝Ｕｓｉｎθ＋Ｖｃｏｓθ …（10）
また
Ｕ＝ｕｃｏｓθ＋Ｖｓｉｎθ
Ｖ＝−ｕｓｉｎθ＋ｖｃｏｓθ …（11）
の関係があり、（10）式よりＦ（ｕ，ｖ）は、
【０１０６】
【数５】
（13）式より、実領域において回転した物体（この場合は音速分布ｆ（ｘ，ｙ))のスペクトルは、元の物体のスペクトルを同じ角度だけ回転したものと同じになることがわかる。
【０１０８】
（７）式及び（８）式によれば、投影データｐ（Ｘ，θ）に関する（４）式は次の（14）式のように表すことができる。
【０１０９】
【数６】
なお、（14）式では音速分布ｆ（ｘ，ｙ）、すなわちｆ'(ｘ，ｙ）が測定対象区間外では０であるため、積分区間を±∞に拡大している。また、投影データｐ（Ｘ，θ）のＸ軸に関するフーリエ変換は、次の（15）式で表すことができる。
【０１１１】
【数７】
（14）式を（15）式に代入すると、次の（16）式が得られる。
【０１１３】
【数８】
（16）式は先の（13）式においてＶ＝０としたときの式に一致するので、
Ｐ（Ｕ，θ）＝Ｆ'(Ｕ，０）
が成り立つ。
【０１１５】
このように、傾きθが同一の全ての仮想経路の投影データｐ（Ｘ，θ）に対して、（15）式に従ってＸ軸に関するフーリエ変換（１次元フーリエ変換）を行った結果Ｐ（Ｕ，θ）は、音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換Ｆ（ｕ，ｖ）のθ方向成分（フーリエ変換Ｆ（ｕ，ｖ）のｕ軸に対して角度θだけ傾いた直線（すなわちＵ軸）上の値）に等しい（これを２次元断層定理（projection slice theorem）という）。
【０１１６】
従って、各方向について各々求めた投影データに対して各々フーリエ変換を行って座標系（ｕ，ｖ）上に放射状に並べれば、音速分布ｆ（ｘ，ｙ）の２次元フーリエ変換Ｆ（ｕ，ｖ）を得ることができる。そして、このＦ（ｕ，ｖ）に対し、先の（６）式に従って２次元フーリエ逆変換を行うことにより、音速分布ｆ（ｘ，ｙ）を求めることができる。
【０１１７】
ステップ１３２では、まず上記で説明したフーリエ領域再構成法、或いは公知のフィルタ処理データの逆投影法や２次元フィルタリング法等の短時間で演算可能な演算方法を適用して投影データｐ（Ｘ，θ）から音速分布ｆ（ｘ，ｙ）を求める。次に、求めた音速分布ｆ（ｘ，ｙ）から（３）式に基づいて被測定空間２４内の各位置における音速ｖ（ｘ，ｙ）を演算し、音速ｖと音波が伝搬する媒体（空気）の温度Ｔとの関係を表す次の（17）式に基づき、被測定空間２４内における温度分布Ｔ（ｘ，ｙ）を演算する。
【０１１８】
ｖ＝ 331.45 ＋ 0.607・Ｔ 〔ｍ／秒〕 …（17）
次のステップ１３４では、上記により得られた被測定空間２４内の温度分布を、例として図１０に示すコンターマップ等の形態でディスプレイ２０に表示し、処理を終了する。なお、ディスプレイ２０への表示に代えて、被空調空間２４内の空気調和を行う空調装置に対し、被測定空間２４内の温度分布及び湿度分布を表すデータを出力するようにしてもよい。これにより、空調装置が被測定空間２４に対し、例えば快適性と省エネルギーの双方を満足する空気調和を行うことも可能となる。
【０１１９】
本実施形態では、被測定空間２４の周縁部に多数のノード１２を設け、仮想経路上の一対のノード１２の間の測定対象区間（一対のノード１２の各々のアーム５４の間の区間）については、実際に音波を両方向について各々伝搬させて伝搬時間ｔmn及び伝搬時間ｔnmを測定し、伝播時間の平均値ｔm-nを求め、測定対象区間外の区間（仮想スピーカとノード１２の間、及びノード１２と仮想マイクロフォンの間の伝搬時間推定区間）については、気温が所定温度であると仮定し、空気が所定温度のときの音速ｖ0 を用いて音波の伝搬時間を推定し、伝搬時間の推定値ｔ（Ｘ，θ）を演算しているので、被測定空間２４外にノード１２を多数設けることなく投影データｐ（Ｘ，θ）を得ることができ、この投影データｐ（Ｘ，θ）に基づき、フーリエ領域再構成法やフィルタ処理データの逆投影法、２次元フィルタリング法等の短時間で演算可能な演算方法を適用して被測定空間内の環境状態を演算することができる。従って、環境状態測定装置１０の構成を簡単にすることができると共に、温度分布の演算に要する時間を短縮することができる。
【０１２０】
なお、上記では第１の位置に到達した音波を第１マイクロフォン３２によって検出する場合を説明したが、音波発生対象のノード１２のスピーカ３０から音波を発生させるための音波データを予め固定的に定めておくと共に、該音波データを用いてスピーカ３０から音波を発生させ、第１の位置に到達した音波をマイクロフォンによって検出し、比較用音波測定データを求めてホストコンピュータ１４のＲＯＭ１４Ｂ、又は各ノード１２のデータ処理部４０の内蔵メモリに記憶しておき、伝播時間の演算時に比較用音波データを読み出して用いることも可能である。
【０１２１】
また、上記では第１マイクロフォン３２及び第２マイクロフォン３４として、同一機種のマイクロフォン（特性が同一のマイクロフォン）を用いていたが、これに限られるものではなく、類似種のマイクロフォン（例えば音波の検出原理が同一等のように特性が類似しているマイクロフォン）であれば第１マイクロフォン３２及び第２マイクロフォン３４としての使用に好適である。
【０１２２】
また、上記では音波の伝搬時間の投影データから温度分布を演算していたが、これに限定されるものではなく、例えば壁体２２の壁面近傍の温度を温度センサによって測定し、壁面近傍の温度の測定値に基づいて前記温度分布を補正するようにしてもよい。
【０１２３】
更に、上記では被測定空間２４を貫くように仮想経路を設定し、仮想経路のうち測定対象区間外の区間については音波の伝播時間を推定することにより、仮想経路上を伝播する音波の伝播時間を求めて温度分布を求めていたが、本発明はこれに限定されるものではなく、例えば被測定空間を仮想的に多数の領域（ボクセル）に区切り、或る一対の箇所の間を音波が伝搬した際に前記一対の箇所の間に存在する各領域を音波が横切る距離を求め、前記一対の箇所の間の音波の伝搬時間と、前記各領域を音波が横切る際の伝搬速度との関係を、前記各領域を音波が横切る距離を用いて数式で表すことを前記複数箇所の間について各々行って連立方程式を立て、測定した複数箇所の間の音波の伝搬時間を連立方程式に代入し、伝搬速度が未知の領域について音波が横切る際の伝搬速度を演算し、伝搬速度から各領域の温度を演算することによって温度分布を求めるようにしてもよい。
【０１２４】
また、上記ではドーム球場等のように床面が円形の被測定空間２４内の環境状態を測定する場合を例に説明したが、本発明によれば、音波発生器及び音波検出手段を被測定空間の周縁部に設けることで被測定空間内の温度分布を求めることができるので、任意の形状の空間を被測定空間として採用し、任意の形状の空間内の温度分布を求めることができる。
【０１２５】
また、上記では環境状態として温度分布のみを演算する場合を例に説明したが、これに限定されるものではなく、複数種の周波数の音波について減衰量を各々測定して減衰率分布を求めることにより、被測定空間内の湿度分布も併せて演算したり、他の環境状態の演算も行うようにしてもよい。
【０１２６】
【発明の効果】
請求項１及び請求項４記載の発明は、第１音波検出手段、第２音波検出手段及び光反射板が、音波発生器と連結されたブラケットに一体に取付けられた音波発生検出装置を、被測定空間を横切る複数の伝播時間測定区間の両端部に、ブラケットが取付けられた側同士が対向するように各々配置し、一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、一方の音波発生検出装置の第１音波検出手段、及び他方の音波発生検出装置の第２音波検出手段によって各々検出させることを、複数の伝播時間測定区間について各々行い、検出された音波の波形を照合し、波形の相関性を時間軸に沿って演算した結果に基づいて、第１音波検出手段による音波の検出位置から第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、複数の伝播時間測定区間について各々行い、複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算するので、被測定空間内の温度分布等の環境状態を高精度に測定できる、という優れた効果を有する。
【０１２８】
請求項２記載の発明は、請求項１記載の発明において、第１音波検出手段による音波検出位置の音波到来側に、到来した音波を減衰させる減衰手段を配置するので、上記効果に加え、第１音波検出手段によって検出される音波の波形が歪んだり飽和することを防止することができる、という効果を有する。
【０１２９】
請求項３記載の発明は、請求項１記載の発明において、第１音波検出手段及び第２音波検出手段が少なくとも類似種の音波検出器によって音波を検出するので、音波伝播時間をより精度良く測定することができる、という効果を有する。
【０１３１】
請求項５記載の発明は、請求項４記載の発明において、伝播時間測定区間の両端部の一対の音波発生検出装置の一対の光反射板に各々光を照射して距離を測定することを、複数の伝播時間測定区間について各々定期的に行い、被測定空間内の環境状態の演算を、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して行うので、上記効果に加え、環境温度の変化や経時変化に伴う伝播時間測定区間の距離変化の影響を受けることなく、被測定空間内の環境状態を精度良く測定することができる、という効果を有する。
【図面の簡単な説明】
【図１】 本実施形態に係る環境状態測定装置の概略構成を示すブロック図である。
【図２】 ノードの配置を示す被測定空間の平面図である。
【図３】 （Ａ）単一のノードにおけるスピーカ及びマイクロフォンの配置を示す側面図、（Ｂ）は（Ａ）におけるマイクロフォン配置位置周辺の拡大図である。
【図４】 測定対象区間及び伝搬時間推定区間から成る仮想経路を、被測定空間の平面図に重ねて示す概念図である。
【図５】 被測定空間に対して設定した直交座標系ｘ−ｙ、各仮想経路を識別するための直線Ｓの距離Ｘ及び傾きθを示す概念図である。
【図６】 （Ａ）及び（Ｂ）は、伝搬時間データの測定順序と、伝搬時間データの処理順序と、の関係を示す概念図である。
【図７】 ホストコンピュータで実行される温度分布演算処理の内容を示すフローチャートである。
【図８】 測定対象の各ノードで実行される音波発生処理及び伝播時間演算処理の内容を示すフローチャートである。
【図９】 （Ａ）は音波データが表す音波の波形の一例、（Ｂ）は比較用音波測定データが表す音波の波形の一例、（Ｃ）は音波測定データが表す音波の波形の一例、（Ｄ）は相関性演算結果の一例を各々示す線図である。
【図１０】 被測定空間内の温度分布の表示例を示すイメージ図である。
【符号の説明】
１０ 環境状態測定装置
１２ ノード
１４ ホストコンピュータ
２４ 被測定空間
３０ スピーカ
３２ 第１マイクロフォン
３４ 第１マイクロフォン
Claims (5)
- 音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられ、被測定空間を横切る複数の伝播時間測定区間の両端部に、前記ブラケットが取付けられた側同士が対向するように各々配置された音波発生検出装置と、
伝播時間測定区間の両端部に配置された一対の前記音波発生検出装置のうちの一方の前記音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の前記音波発生検出装置の前記第１音波検出手段、及び他方の前記音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行う制御手段と、
前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行う伝播時間演算手段と、
前記伝播時間測定手段によって測定された前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算する環境状態演算手段と、
を含む環境状態測定装置。 - 前記音波発生検出装置は、前記第１音波検出手段による音波検出位置の音波到来側に配置され、到来した音波を減衰させる減衰手段を更に備えたことを特徴とする請求項１記載の環境状態測定装置。
- 前記音波発生検出装置の前記第１音波検出手段及び第２音波検出手段は、少なくとも類似種の音波検出器によって音波を検出することを特徴とする請求項１記載の環境状態測定装置。
- 被測定空間を横切る複数の伝播時間測定区間の両端部に、音波発生器と、前記音波発生器の音波放射側に配置され、音波発生器から放射された音波を音波発生器から所定距離隔てた位置で検出する第１音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した音波を音波発生器から所定距離隔てた位置で検出する第２音波検出手段と、音波発生器の音波放射側に配置され、音波発生器側と反対側より到来した光を音波発生器から所定距離隔てた位置で反射する光反射板と、を含み、前記第１音波検出手段、第２音波検出手段及び前記光反射板が、前記音波発生器と連結されたブラケットに一体に取付けられた音波発生検出装置を、前記ブラケットが取付けられた側同士が対向するように各々配置しておき、
伝播時間測定区間の両端部に配置した一対の音波発生検出装置のうちの一方の音波発生検出装置の音波発生器から音波を発生させ、発生させた音波を、前記一方の音波発生検出装置の前記第１音波検出手段、及び他方の音波発生検出装置の前記第２音波検出手段によって各々検出させることを、前記複数の伝播時間測定区間について各々行い、
前記第２音波検出手段によって検出された音波の波形を前記第１音波検出手段によって検出された音波の波形と照合し、前記波形の相関性を時間軸に沿って演算した結果に基づいて、前記第１音波検出手段による音波の検出位置から前記第２音波検出手段による音波の検出位置へ伝播する音波の伝播時間を演算することを、前記複数の伝播時間測定区間について各々行い、
前記複数の伝播時間測定区間の音波の伝播時間の測定値に基づいて被測定空間内の環境状態を演算する
環境状態測定方法。 - 伝播時間測定区間の両端部に各々配置した一対の音波発生検出装置の一対の光反射板に各々光を照射して前記一対の光反射板の距離を測定することを、前記複数の伝播時間測定区間について各々定期的に行い、
前記被測定空間内の環境状態の演算を、複数の伝播時間測定区間の各々における一対の光反射板の距離の測定結果も考慮して行う
ことを特徴とする請求項４記載の環境状態測定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000244781A JP4783873B2 (ja) | 2000-08-11 | 2000-08-11 | 環境状態測定装置及び方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000244781A JP4783873B2 (ja) | 2000-08-11 | 2000-08-11 | 環境状態測定装置及び方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002054999A JP2002054999A (ja) | 2002-02-20 |
| JP4783873B2 true JP4783873B2 (ja) | 2011-09-28 |
Family
ID=18735417
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000244781A Expired - Fee Related JP4783873B2 (ja) | 2000-08-11 | 2000-08-11 | 環境状態測定装置及び方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4783873B2 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5945304B2 (ja) * | 2014-08-20 | 2016-07-05 | 株式会社東日製作所 | 締付工具 |
| JP6771681B2 (ja) * | 2017-10-11 | 2020-10-21 | 三菱電機株式会社 | 空調用コントローラ |
| CN107917768B (zh) * | 2017-10-18 | 2020-02-21 | 华北电力大学 | 一种基于低频声波的空气温度测量装置及方法 |
| WO2021145440A1 (ja) * | 2020-01-17 | 2021-07-22 | ピクシーダストテクノロジーズ株式会社 | 計測装置、計測方法、及び、プログラム |
| JP6927449B1 (ja) * | 2020-04-30 | 2021-09-01 | ダイキン工業株式会社 | 環境検出システム |
| CN111693139B (zh) * | 2020-06-19 | 2022-04-22 | 浙江讯飞智能科技有限公司 | 声音强度测量方法、装置、设备及存储介质 |
| CN115950077A (zh) * | 2022-12-23 | 2023-04-11 | 珠海格力电器股份有限公司 | 基于语音空调的温度控制方法、装置、空调及存储介质 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60120246A (ja) * | 1983-12-05 | 1985-06-27 | Fujitsu Ltd | 音響信号検出方法 |
| JPH0715782A (ja) * | 1993-06-28 | 1995-01-17 | Matsushita Electric Ind Co Ltd | スピーカユニット |
| JP3421124B2 (ja) * | 1994-05-17 | 2003-06-30 | 株式会社小野測器 | 音響測定方法および装置 |
| JP3324876B2 (ja) * | 1994-07-15 | 2002-09-17 | ティーオーエー株式会社 | プローブマイクロホン |
| JPH09243367A (ja) * | 1996-03-13 | 1997-09-19 | Mitsubishi Heavy Ind Ltd | 一脚距離計測装置 |
| JP3602951B2 (ja) * | 1997-12-15 | 2004-12-15 | 株式会社竹中工務店 | 環境状態測定方法及び装置 |
| JP3699269B2 (ja) * | 1998-02-04 | 2005-09-28 | 株式会社竹中工務店 | 環境状態測定方法及び装置 |
| JP3613971B2 (ja) * | 1998-03-26 | 2005-01-26 | いすゞ自動車株式会社 | 伝達特性測定装置 |
-
2000
- 2000-08-11 JP JP2000244781A patent/JP4783873B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002054999A (ja) | 2002-02-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5509414A (en) | Apparatus and method for non-contacting detection of respiration | |
| Lehnert et al. | Principles of binaural room simulation | |
| US9194938B2 (en) | Time difference of arrival determination with direct sound | |
| US6673014B2 (en) | Noninvasive methods and apparatuses for measuring the intraocular pressure of a mammal eye | |
| Katz | Acoustic absorption measurement of human hair and skin within the audible frequency range | |
| JP4783873B2 (ja) | 環境状態測定装置及び方法 | |
| JPH051890B2 (ja) | ||
| US7920709B1 (en) | Vector sound-intensity probes operating in a half-space | |
| DK2502035T3 (en) | A method and apparatus for determining coefficients of acoustic and acoustic power | |
| Achdjian et al. | Source localisation in a reverberant plate using average coda properties and early signal strength | |
| US8958269B2 (en) | Transducer for phased array acoustic systems | |
| Stalnov et al. | Auditory detection probability of propeller noise in hover flight in presence of ambient soundscape | |
| US11249187B2 (en) | Method and acoustic device for measuring surface movements | |
| JP3699269B2 (ja) | 環境状態測定方法及び装置 | |
| JP3602951B2 (ja) | 環境状態測定方法及び装置 | |
| JP3977559B2 (ja) | 伝播時間測定方法、環境状態測定方法及び装置 | |
| US6449566B1 (en) | Acoustic scattering measurement and processing for determining variances in multiple features | |
| Jacobsen et al. | The influence of microphone vents on the performance of sound intensity probes | |
| CN111198282A (zh) | 用于校准扬声器的集成体积加速度传感器的方法和系统 | |
| WO2021145440A1 (ja) | 計測装置、計測方法、及び、プログラム | |
| Terashima et al. | Acoustic characteristics of a miniature dynamic speaker driver unit MT006B for measurement of head-related impulse responses by reciprocal method | |
| JP2018523504A (ja) | カテーテルデバイス | |
| WO2022150950A1 (zh) | 评估传声器阵列一致性的方法和装置 | |
| JP2019032446A (ja) | 音響シミュレーション方法、装置、及びプログラム | |
| JP2700530B2 (ja) | 音響伝達特性測定装置および方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070807 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101019 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101220 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110125 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20110224 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110224 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20110225 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20110428 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140722 Year of fee payment: 3 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |