JP2004077352A - Geometric acoustic simulation system, geometric acoustic simulation method, and program - Google Patents
Geometric acoustic simulation system, geometric acoustic simulation method, and program Download PDFInfo
- Publication number
- JP2004077352A JP2004077352A JP2002240085A JP2002240085A JP2004077352A JP 2004077352 A JP2004077352 A JP 2004077352A JP 2002240085 A JP2002240085 A JP 2002240085A JP 2002240085 A JP2002240085 A JP 2002240085A JP 2004077352 A JP2004077352 A JP 2004077352A
- Authority
- JP
- Japan
- Prior art keywords
- reflection
- inner product
- sound
- sound wave
- space
- 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.)
- Granted
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
【０００１】
【発明の属する技術分野】
この発明は、コンサートホールや劇場、スタジオ等の建造物の音場の音響特性を解析するために幾何音響シミュレーションを行うための装置及びその方法に関し、特に、その幾何音響シミュレーションにおける音波の反射方向の演算を高速化する技術に関する。
【０００２】
【従来の技術】
コンサートホール等の建造物の音場である所定の空間内の音響特性を解析することは、残響を付加するためのデータを測定したり、建築前にホールの音響特性をシミュレーションしてその情報を提供するなどのために必要になる。
そのためには、対象となる空間内のある一点の音源から発せられた音が、到達する収音点まで複数の壁や天井などの反響（音の反射）に関与する面をどのような経路で伝わるかを解析することが、必要不可欠となる。
このような解析手法はいくつか知られているが、音の波動性を考慮するか否かで２つに大別できる。このうち、音の波動性を考慮しない幾何学的な手法には虚像法と音線追跡法が知られている。
【０００３】
虚像法は、光の鏡面反射を解析する場合と同様に、音源から反射面に垂線を引き、この垂線の長さと同じ距離だけ反射面から奥に離れた点を虚音像とし、この虚音像と収音点を結んで反射音の経路を得る方法である。
一方、音線追跡法は、図６に示すように音源１を中心に発せられる音を、一定のスピードで直線状に進む線でモデル化した多数の音線２で表現し、各音線２の壁面３などによる反射を１本ずつ追跡して、所定の収音点４を中心とする収音エリア５に到達した音線２のみを反射音として採用する方法である。
虚像法によれば正確な反射経路を求めることができるが、この方法は反射次数（反射回数）が大きくなると演算時間が指数的に増加するため、実用性に乏しい。
【０００４】
そこで、ホール等の音場の音響特性を解析したり、その内部での音響を予測あるいは再現したりするために、長い遅れ時間を持つ反射音まで解析することが必要な場合には、実用的な観点から一般的に音線追跡法が用いられている。
しかし、従来の一般的な音線追跡法では、一定の半径を持つ球形の収音エリアに入った音線を全て反射音として採用していたため、同じ反射経路を通った複数の音線が反射音として採用されてしまい、演算上の誤差が大きくなるという問題があった。
【０００５】
そこで、特開平９−１６６４８２号公報には、この問題を解決する高精度音線追跡装置及び高精度音線追跡方法が提案されている。その音線追跡方法によれば、受音エリア（収音エリア）の大きさを２段階に変化させると共に、反射経路の履歴を照合して同一経路の音線を反射音として複数採用した場合には、１本を残して他の音線を除去するようにしている。
【０００６】
【発明が解決しようとする課題】
しかしながら、上述の公報には、音線１本１本の反射経路の算出に関しては特に記載されておらず、例えば、「壁で反射を繰り返した音線の経路探索」について具体的な方法を開示していないが、これは以下のような一般的に方法で行われているものと思料される。この従来の反射経路探索処理について、図７を用いて説明する。
シミュレーション中の音線の各線分には、その方向を示す音線ベクトルがその起点の位置情報と共に与えられている。そして、音源から発せられた音線又は壁（反射面）で反射された音線が次にどの位置で反射されるかを演算する際に、図７のフローチャートに示す処理を行う。
【０００７】
まず、ステップ▲１▼で、シミュレーション系において定義されている全ての面の中から一つの面を選択する。そして、ステップ▲２▼へ進んでその面の法線ベクトルと演算対象の音線の音線ベクトルとの内積を計算する。ここで、面の法線ベクトルは、音線が到達し得る側（表側：建物内部のシミュレーションの場合には内側）の向きのベクトルを使用するものとする。
次に、ステップ▲３▼でその内積が負かどうか判定し、負であれば音線は表側からその面に向かっていると判断してステップ▲４▼へ進み、音線とその面が交わっているかどうか判断する。交わっていればステップ▲５▼へ進み、その面を反射面の候補として記憶する。ここで、交点の座標や音線の線分の起点から交点までの経路長を共に記憶してステップ▲６▼へ進む。
【０００８】
ステップ▲３▼で負でなかった場合には、音線は裏側（建物内部のシミュレーションの場合には外側）からその面に向かっているものと判断し、実際にはその側から音線が面に向かうことはありえないので、その面を反射面の候補とせずにステップ▲６▼へ進む。ステップ▲４▼で交わっていなかった場合にも、その面は反射面とはなりえないので候補とせずにステップ▲６▼へ進む。
ステップ▲６▼では、全ての面について調べたかどうか判断し、調べていなければステップ▲８▼へ進んで次の面を選択し、ステップ▲２▼に戻って同じ処理を繰り返す。全ての面について調べていた場合にはステップ▲７▼へ進み、記憶した反射面の候補の中で交点までの経路が最も短い面を反射面として採用して反射点の位置や経路長を記憶し、反射後の音線の進行方向を求めて終了する。
この処理によって、音線が次に反射される面とその反射点の位置が求められ、それを基に反射後の音線の起点と音線ベクトルが求められるので、その値を用いて同様の処理を繰り返すことにより、音線の反射経路を追跡することができる。
【０００９】
上述の処理について、図８及び図９を用いて具体的に説明する。まず、図８によって二つのベクトルの内積について説明する。
図８に示すベクトルＡとベクトルＢの内積（スカラー積ともいう）は、その二つのベクトルＡとＢがなす角がθであると、｜Ａ｜・｜Ｂ｜ｃｏｓθである。これを便宜上（ＡＢｃｏｓθ）とする。この値の正負の符号は、ｃｏｓθで決まる。
０≦θ≦π（１８０°）であるから、内積（ＡＢｃｏｓθ）の正負は次のようになる。
０≦θ＜π／２ すなわちθが鋭角であれば正
θ＝π／２ すなわちθが直角であれば０
π／２＜θ≦π すなわちθが鈍角であれば負
このθが鈍角の場合のベクトルＡと角度θを、図８では仮想線のべクトルＡ′と角度θ′で示している。
【００１０】
そこで、図９に示すように、一つの面６を選択し、その内側向きの法線ベクトルＤと演算対象の音線２の音線ベクトルＳとの内積（ＤＳｃｏｓθ）を計算する。このとき音線２が図９に実線で示すように内側から面６に向かっていれば、その音線ベクトルＳと法線ベクトルＤとのなす角θは鈍角になるので、内積（ＤＳｃｏｓθ）は負になる。音線２が図９に仮想線で示すように外側から面６に向かっていれば、その音線ベクトルＳ′と法線ベクトルＤとのなす角θ′は鋭角になるので、内積（ＤＳ′ｃｏｓθ′）は正になる。しかし、実際には音線が外側から面６に向かうことはないので、その面６は反射面の候補としない。θが９０°（内積の計算値が０）ということも反射面としては実際にはありえないので、この場合もその面を反射面の候補としない。したがって、内積の計算値が正であることが、その面を反射面の候補とするために必要な条件である。
また、図９に示すように、音線２と面６とが交わっていることも、その面を反射面の候補とするために必要な条件である。
【００１１】
また、この種の従来技術として、特開平８−２７２３８０号公報に記載されている仮想３次元空間音響の再生方法もある。この方法によると、複数の音響（音線）ベクトルを設定し、仮想空間を多角形の境界から構成し、所定の継続時間内に、音波がその境界により反射して伝播する音響ベクトルの伝播履歴データを算出記憶し、その伝播履歴データをもとに、各音響ベクトルについて、境界に反射する音波とその音波の微少面素によってつくられる観測点への速度ポテンシャルから、観測点に到達する過渡応答を、時系列の数値配列に加算記憶する。
この場合も、音波の反射に寄与する境界であるかどうか（界面が表か裏か）を判定するために、境界面ごとの法線ベクトルと音線ベクトルとの内積の計算が必要であることが、上記公報の段落００３２に記載されている。
【００１２】
このように、従来の音響シミュレーションにおいては、ある境界面が音の反射に寄与する面であるかどうか（境界面が表か裏か）を判定するために、境界面ごとにその法線ベクトルと音線ベクトルとの内積演算が必要であった。そのために、１回の反射を求める毎に膨大な量の演算が必要になり、シミュレーションの高速化を妨げる原因のひとつとなっていた。
この発明は、このような従来の問題を解決するためになされたものであり、所定の空間内の音源から発せられる音波に関して、その空間内における反射方向を演算することによって反射の経路を解析する幾何音響シミュレーションを、その演算量を少なくして高速で行えるようにすることを目的とする。
【００１３】
【課題を解決するための手段】
この発明による幾何音響シミュレーション装置は、所定の空間内の音源から発せられる音波の該空間内における反射方向を反射方向演算手段によって演算することによって上記音波の反射経路を解析する幾何音響シミュレーション装置であり、上記の目的を達成するため、次のように構成したものである。
すなわち、上記空間内又は該空間の境界において上記音波の反射に関与する各面に関する法線ベクトルを記憶する法線ベクトル記憶手段と、予め設定した基準となる複数の基準ベクトルのそれぞれと上記法線ベクトルとの内積を求めて、その内積の符号の組み合わせに応じて上記音波の反射に関与する面を複数のグループに分類（クラスタリング）する面分類手段とを設ける。
【００１４】
そして、上記反射方向演算手段を、上記音波の反射方向を演算する際に、その音波の音線ベクトルと上記複数の基準ベクトルのそれぞれとの内積を求め、上記複数のグループのうち該内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算する手段としたことを特徴とする。これによって、音線追跡時に反射の有無を判定する際、上記グループ（クラスタ）単位で、反射の有無の判定対象となる面を削減できる。
上記所定のグループは、上記音線ベクトルと分類される面の法線ベクトルとの内積が所定の符号になるようなグループであるのがよい。
【００１５】
この発明による幾何音響シミュレーション方法は、所定の空間内の音源から発せられる音波の該空間内における反射方向を演算することによって反射の経路を解析する幾何音響シミュレーション方法であって、上記の目的を達成するため、次のようにシミュレーションを行なう。
上記空間内又は該空間の境界において上記音波の反射に関与する面の法線ベクトルを記憶し、予め設定した基準となる複数の基準ベクトルのそれぞれと上記法線ベクトルとの内積を求めて、その内積の符号の組み合わせに応じて上記音波の反射に関与する面を複数のグループに分類（クラスタリング）する。
【００１６】
そして、上記音波の反射方向を演算する際に、その音波の音線ベクトルと上記複数の基準ベクトルのそれぞれとの内積を求め、上記複数のグループのうち該内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算することを特徴とする。
上記所定のグループを、上記音線ベクトルと分類される面の法線ベクトルとの内積が所定の符号になるグループとするとよい。
これによって、音線追跡時に反射の有無を判定する際、上記グループ（クラスタ）単位で、その判定対象となる面を削減できる。
【００１７】
また、この発明のプログラムは、音波の反射経路を解析する幾何音響シミュレーション装置を制御するコンピュータを、所定の空間内の音源から発せられる音波のその空間内における反射方向を演算する反射方向演算手段と、上記空間内又はその空間の境界において上記音波の反射に関与する各面に関する法線ベクトルを記憶する法線ベクトル記憶手段と、予め設定した基準となる複数の基準ベクトルのそれぞれと上記法線ベクトルとの内積を求めてその内積の符号の組み合わせに応じて上記音波の反射に関与する面を複数のグループに分類する面分類手段として機能させるプログラムであり、上記反射方向演算手段の機能は、上記音波の反射方向を演算する際に、その音波の音線ベクトルと上記複数の基準ベクトルのそれぞれとの内積を求め、上記複数のグループのうちその内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算する機能である。
【００１８】
【発明の実施の形態】
以下、この発明の好ましい実施の形態を図面を参照して説明する。
まず、この発明による幾何音響シミュレーション装置の実施形態について、図１を用いて説明する。図１は、その幾何音響シミュレーション装置の構成を示すブロック図である。
この図１に示す幾何音響シミュレーション装置は、ＣＰＵ１１，プログラムメモリ１２，面情報入力手段１３，ＲＡＭ１４，結果出力手段１５がシステムバス１６によって接続されている。
ＣＰＵ１１は、装置全体の制御及びシミュレーションの演算を行う演算制御部であり、プログラムメモリ１２に格納された制御・演算プログラムを読み出して実行する。
【００１９】
プログラムメモリ１２は、ＲＯＭやハードディスクによって構成され、ＣＰＵ１１で実行するための各種プログラムを記憶する。このプログラムには、ＣＰＵに幾何音響シミュレーションを実行させるためのプログラムも含まれる。
面情報入力手段１３は、シミュレーションを実行する仮想空間内あるいはその空間の境界に配置される音の反射に関与する面の情報を入力する手段である。建物内部の音響についてシミュレーションを行うのであれば、壁面や天井および床などの情報となる。
面の情報としては、位置と法線ベクトルを入力するものとし、曲面については多数の平面で近似して、その各位置と法線ベクトルで表現するものとする。通常、面の位置は座標に基づいた方程式で表わされる。反射率を考慮する場合にはその情報も共に入力するようにしてもよいが、ここでは説明を簡単にするために反射率は考慮しないものとする。入力方法としては、キーボード等の入力手段によってユーザが直接入力するか、所定の形式で作成したファイルを読み込むようにする等の方法が考えられる。
【００２０】
ＲＡＭ１４は、面情報入力手段１３から入力された面情報を記憶したり、ＣＰＵ１１のワークメモリとして使用したり、シミュレーションの結果を記憶したりするために用いる記憶手段である。
結果出力手段１５は、シミュレーションの結果を出力する手段であり、ディスプレイ、プリンタ、ハードディスクや、外部記録媒体への書き込みデバイスなどによって構成することができる。
このような幾何音響シミュレーション装置によって、この発明の幾何音響シミュレーション方法を実行することができるが、この発明による幾何音響シミュレーション装置は、汎用のパーソナルコンピュータに必要なプログラムをインストールして構成することもできる。
【００２１】
次に、以上説明した幾何音響シミュレーション装置の動作、およびその装置を用いて実施するこの発明による幾何音響シミュレーション方法の実施形態について、図２乃至図５を用いて説明する。
図２はその幾何音響シミュレーション方法における基準ベクトル及びその基準ベクトルによるベクトルの分類について説明するための図、図３は上記の幾何音響シミュレーション装置における面の法線ベクトルの分類の処理を示すフローチャート、図４は同じく幾何音響シミュレーションの処理を示すフローチャート、図５はその音線の反射経路追跡処理の詳細を示すフローチャートである。
【００２２】
この発明による幾何音響シミュレーション装置及び幾何音響シミュレーション方法の特徴は、音の反射に関与する各面（ｎ個）の０でない法線ベクトル（Ｖ１，Ｖ２，Ｖ３，……Ｖｎとする）を予め複数のグループに分類（クラスタリング）しておくことにより、音線の経路探索の際に所定のグループの面のみを考慮して探索を行うことができるようにした点であるので、この点についてまず説明する。ただし、実際のシミュレーションは３次元空間を想定して行うが、ここでは説明を簡単にするために２次元のベクトルを用いて説明する。
【００２３】
例えば、図２に示すように、基準となる基準ベクトルｎ１，ｎ２，ｎ３，ｎ４を４５°間隔で定めておく。このとき問題なのは、基準ベクトルｎ１，ｎ２，ｎ３，ｎ４の関係であり、絶対的な向きは都合のよいように決めればよい。
その基準ベクトルｎ１，ｎ２，ｎ３，ｎ４のそれぞれと各面の法線ベクトル（Ｖ１，Ｖ２，Ｖ３，……Ｖｎ）との内積を演算すると、各面はその基準ベクトルと法線ベクトルとのなす角によって決まる内積の符号（＋，−，０）の組み合わせによって、表１に示すように１６のグループに分類することができる。
【００２４】
【表１】
これらの各グループに分類される面の法線ベクトルの向きは、図２に示すＸ−Ｙ座標の原点Ｏを中心に１から１６までの番号を付した向きである。
経路探索の対象となる音線の進行方向を示す０でない音線ベクトル（ｘとする）も、通常は各法線ベクトル（Ｖ１，Ｖ２，Ｖ３，……Ｖｎ）の全てとの内積を計算するが、本実施例では基準ベクトルｎ１，ｎ２，ｎ３，ｎ４の各々との内積を演算し、その符号の組み合わせによって同様な１６のグループに分類することができる。そうすると、音線ベクトルの属するグループと面の属するグループとの関係によっては、音線ベクトルとその面の法線ベクトルとの内積の符号が、実際に内積の演算を行わなくてもわかる場合がある。
例えば、ある面ｉの法線ベクトルｖｉが、表１において全ての基準ベクトルｎ１，ｎ２，ｎ３，ｎ４との内積が正（＋）であるグループ１に分類されていた場合、図２から明らかなように、その法線ベクトルｖｉは、数１のようにベクトルｎ２とベクトルｎ３の正の係数による線形結合で表現することができる。
【００２６】
【数１】
一方、ある音線ベクトルｘが、表１において基準ベクトルｎ１，ｎ２，ｎ３との内積が正で基準ベクトルｎ４との内積が負であるグループ２に分類される場合、図２から明らかなように、ベクトルｘは数２のようにベクトルｎ１とベクトルｎ２の正の係数による線形結合で表現することができる。
すると、法線ベクトルｖｉの内積と音線ベクトルｘとの内積は、数３に示すとおり、具体的に数値を演算するまでもなく正であることがわかる。
【００２８】
【数２】
【数３】
なお、ここで重要なのは基準ベクトルｎ１，ｎ２，ｎ３，ｎ４の向きの関係であり、その絶対的な方向は異なっていても、同様な分類によって同様な符号の推定を行うことができる。
以上、２次元の場合について説明したが、３次元にも同様な手法を拡張することができ、基準ベクトルの各々と法線ベクトルとの内積の符号の組み合わせによって各面を複数のグループに分類しておくと、音線ベクトルと基準ベクトルの各々との内積の符号の組み合わせに応じて、いくつかのグループの面については、実際に内積を演算するまでもなくその法線ベクトルとの内積の符号を確定することができる。
【００３１】
従って、シミュレーションの過程において内積が０以上であることが必要な場合には、内積が負であることが予めわかるグループの面については、実際に内積を演算することなく候補から除外することができる。また、内積が正であることが必要な場合には内積が０又は負であることがわかるグループ、負であることが必要な場合には０又は正であることがわかるグループについては、実際に内積を演算することなく候補から除外することができる。但し、内積の正負がわからない場合は候補から除外できない。
これにより、シミュレーション実行時の内積の演算量を大きく削減することができる。
【００３２】
次に、前述した幾何音響シミュレーション装置の動作及びこの発明の幾何音響シミュレーション方法について説明する。
図１によって前述した幾何音響シミュレーション装置は、シミュレーションを実行する前に、面情報入力手段１３から音の反射に関与する面の面情報（上述のように法線ベクトルの情報を含む）を読み込んでＲＡＭ１４に記憶する。ここでは、ＲＡＭ１４が法線ベクトル記憶手段として機能する。
そして、その後にユーザから事前処理開始の指示があると、ＣＰＵ１１は図３のフローチャートに示す事前処理を開始する。
【００３３】
まず、ステップＳ１で面の一つを最初の面として選択し、ステップＳ２で選択された面の法線ベクトルと所定の複数の基準ベクトルの各々との内積を演算する。そして、ステップＳ３に進んで求めた内積の符号の組み合わせによって面を分類し、その結果をＲＡＭ１４に記憶する。ここでは、ＣＰＵ１１が法線ベクトル分類手段として機能する。また、面の法線ベクトルは、音線の到達し得る側（表側：建物内部のシミュレーションの場合には内側）の向きのベクトルを使用するものとする。
その後、ステップＳ４に進み、全ての面について調べたかどうかを判断し、まだしていなければステップＳ５に進んで次の面を選択してステップＳ２に戻って処理を繰り返し、分類が完了していれば処理を終了する。
【００３４】
この事前処理により、音の反射に関与する全ての面を、その法線ベクトルの向きによって分類（クラスタリング）することができる。
この後、ユーザがシミュレーションの開始を指示すると、ＣＰＵ１１は図４のフローチャートに示す処理を開始する。
まず、ステップＳ１１で音源の位置，収音点の位置，収音エリアの大きさ，音線の刻み角，音線の追跡時間等のシミュレーション条件を設定する。
そして、ステップＳ１２に進んで初めの音線を選択し、ステップＳ１３でその音線の進行方向を示す音線ベクトルを設定し、開始点を音源の位置に設定して音線を発射する処理を行う。
そして、ステップＳ１４に進んでその音線について１段階の反射経路追跡処理を行う。このサブルーチンの処理は、音線が次に反射される位置と反射後の音線の進行方向を求める処理であり、図５にその処理内容をフローチャートで示している。
【００３５】
この音線の反射経路追跡処理時には、ＣＰＵ１１が反射方向演算手段として機能し、まず図５のステップＳ２１で音線ベクトルと基準ベクトルの各々との内積を演算する。そして、ステップＳ２２に進んで求めた内積の符号の組み合わせに従って音線ベクトルを分類する。
そして、ステップＳ２３に進んでその分類結果に従って探索対象の面を限定する。すなわち、この実施形態においては、法線ベクトルは面の表側を向いているものを採用しているため、法線ベクトルと音線ベクトルとの内積が負になる面のみが、追跡対象の音線が表側から到達し得る面である。従って、音線ベクトルが分類されたグループとの関係において内積が正又は０になることがわかっているグループに分類されている面については、予め反射面の候補から除いてしまっても問題ない。なお、上記内積の正負が不明な面は候補から除外できない。
【００３６】
次に、ステップＳ２４で残った面の中から最初に探索する面を選択する。そして、ステップＳ２５で選択された面について、その法線ベクトルと音線ベクトルとの内積を演算する。
そして、ステップＳ２６でその内積が負かどうか判定し、負であれば音線は表側からその面に向かっていると判断してステップＳ２７に進み、音線とその面が交わっているかどうか判断する。交わっていれば、ステップＳ２８に進み、その面を反射面の候補としてＲＡＭ１４に記憶する。ここで、交点の座標や、音線の線分の起点から交点までの経路長を共に記憶してステップＳ２９に進む。
【００３７】
ステップＳ２６で負でなかった場合には、音線は裏側（建物内部のシミュレーションの場合には外側）からその面に向かっているものと判断して、実際にはその側から音線が面に向かうことはありえないため、その面を反射面の候補とせずにステップＳ２９に進む。ステップＳ２７で交わっていなかった場合にも、その面は反射面とはなりえないので候補とせずにステップＳ２９に進む。
ステップＳ２９では、ステップＳ２３で限定した面の全てについて調べたかどうか判断し、調べていなければステップＳ３１に進んで次の面を選択し、ステップＳ２５に戻って処理を繰り返す。調べていた場合には、ステップＳ３０に進み、記憶した反射面の候補の中で交点までの経路が最も短い面を反射面として採用して反射点の位置や経路長を記憶し、反射後の音線の進行方向を求めて図４のステップＳ１５に戻る。
反射後の音線の進行方向を示す音線ベクトルｘ′は、もとの音線ベクトルをｘ，反射面の法線ベクトルを単位ベクトルｖｉとすると、数４によって求めることができる。
【００３８】
【数４】
なお、図５のフローチャートの処理において、法線ベクトルと音線ベクトルとの内積が負になることがわかっているグループに分類された面について探索する場合には、ステップＳ２６を省略し、ステップＳ２５に代えてその面がステップＳ３０で反射面として採用されたあとにその内積を演算するようにしてもよい。このようにすれば、演算量をさらに少なくすることができる。
図４のステップＳ１５に戻ると、音線が求めた収音エリアに到達したかどうか判断し、到達していればステップＳ１６でその音線の情報としてその経路や到達時間をＲＡＭ１４に記憶して、ステップＳ１７に進む。到達していなければそのままステップＳ１７に進む。そして、ステップＳ１７では、音源からの経路長の総計によって反射点までの所要時間を算出し、それが追跡時間内かどうか判断する。追跡時間内であれば、ステップＳ１４に戻って処理を繰り返し、追跡時間内でなければステップＳ１８に進む。
【００４０】
ステップＳ１８では、全ての音線の追跡が終了したかどうか判断し、終了していれば処理を終了する。終了していなければ、ステップＳ１９に進んで次の音線を選択し、ステップＳ１３に戻って処理を繰り返す。
このような処理によってシミュレーションを実行することにより、音線の経路探索時に、音波の反射に関与する面の一部について、その法線ベクトルと音線ベクトルとの内積を実際に求めることなく反射面の候補から除くことができるため、演算を高速化することができる。
【００４１】
【発明の効果】
以上説明してきたように、この発明による幾何音響シミュレーション装置及び幾何音響シミュレーション方法によって、建築物の音場などの所定空間の音響シミュレーションを行なえば、音線の経路探索時に、音波の反射に関与する面の一部については、その法線ベクトルと音線ベクトルとの内積を実際に求めることなく、反射面の候補から除くことができるため、多数の音線の経路情報を得るための演算を高速化することができる。
【図面の簡単な説明】
【図１】この発明の実施形態である幾何音響シミュレーション装置の構成を示すブロック図である。
【図２】この発明の幾何音響シミュレーション方法における複数の基準ベクトルと法線ベクトルの向きについて説明するための図である。
【図３】この発明の幾何音響シミュレーション装置における反射に関与する面の法線ベクトルによる分類の処理を示すフローチャートである。
【図４】同じく幾何音響シミュレーションの処理を示すフローチャートである。
【図５】その音線の反射経路追跡処理の詳細を示すフローチャートである。
【図６】従来の音線追跡法について説明するための図である。
【図７】従来の音線の反射経路追跡処理の一部を示すフローチャートである。
【図８】二つのベクトルの内積についての説明図である。
【図９】選択した面とその法線ベクトルおよびその面に向かう音線の音線ベクトルの関係を示す説明図である。
【符号の説明】
１１…ＣＰＵ、１２…プログラムメモリ、１３…面情報入力手段、１４…ＲＡＭ１５…結果出力手段、１６…システムバス[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an apparatus and a method for performing a geometric acoustic simulation for analyzing acoustic characteristics of a sound field of a building such as a concert hall, a theater, and a studio. The present invention relates to a technique for speeding up computation.
[0002]
[Prior art]
Analyzing the acoustic characteristics in a given space, which is the sound field of a building such as a concert hall, measures data for adding reverberation, or simulates the acoustic characteristics of the hall before construction to obtain that information. It becomes necessary for providing.
To achieve this, the sound emitted from a single sound source in the target space travels through multiple walls, ceilings, and other surfaces involved in reverberation (reflection of sound) through any path to the sound collection point that reaches it. Analyzing the transmission is essential.
Although several such analysis methods are known, they can be broadly classified into two types depending on whether or not the wave nature of sound is considered. Among them, a virtual image method and a ray tracing method are known as geometric methods that do not consider the wave nature of sound.
[0003]
In the virtual image method, as in the case of analyzing the specular reflection of light, a perpendicular line is drawn from the sound source to the reflecting surface, and a point away from the reflecting surface by the same distance as the length of the perpendicular line is regarded as a virtual sound image. This is a method of obtaining a reflected sound path by connecting sound collection points.
On the other hand, in the sound ray tracing method, as shown in FIG. 6, a sound emitted from the
According to the virtual image method, an accurate reflection path can be obtained, but this method is not practical because the calculation time exponentially increases as the reflection order (the number of reflections) increases.
[0004]
Therefore, if it is necessary to analyze the reflected sound with a long delay time to analyze the acoustic characteristics of the sound field such as a hall or to predict or reproduce the sound inside it, From the viewpoint, the ray tracing method is generally used.
However, in the conventional general ray tracing method, all sound rays that enter a spherical sound collection area with a fixed radius are adopted as reflected sounds, so that multiple sound rays passing through the same reflection path are reflected. There is a problem in that it is adopted as a sound, and a calculation error increases.
[0005]
Therefore, Japanese Patent Application Laid-Open No. 9-166482 proposes a high-accuracy sound ray tracing device and a high-accuracy sound ray tracing method that solve this problem. According to the sound ray tracking method, when the size of the sound receiving area (sound collecting area) is changed in two stages, the history of the reflection path is collated, and a plurality of sound rays of the same path are adopted as the reflection sound. Is to remove one sound ray except one.
[0006]
[Problems to be solved by the invention]
However, the above-mentioned publication does not specifically describe the calculation of the reflection path of each sound ray, and discloses, for example, a specific method for “search for a sound ray path that has been repeatedly reflected on a wall”. Although not done, it is believed that this is done in the following general way: This conventional reflection path search processing will be described with reference to FIG.
A sound ray vector indicating the direction is given to each line segment of the sound ray being simulated together with the position information of the starting point. Then, when calculating where the sound ray emitted from the sound source or the sound ray reflected on the wall (reflection surface) will be reflected next, the processing shown in the flowchart of FIG. 7 is performed.
[0007]
First, in step (1), one surface is selected from all surfaces defined in the simulation system. Then, the process proceeds to step (2) to calculate the inner product of the normal vector of the surface and the sound ray vector of the sound ray to be calculated. Here, as the normal vector of the surface, a vector in the direction on the side to which the sound ray can reach (front side: in the case of the simulation inside the building) is used.
Next, in step (3), it is determined whether or not the inner product is negative. If the inner product is negative, it is determined that the sound ray is directed from the front side to the surface, and the process proceeds to step (4), where the sound ray and the surface intersect. Judge whether it is. If they do, the process proceeds to step (5), and the surface is stored as a reflection surface candidate. Here, the coordinates of the intersection and the path length from the starting point of the sound ray line segment to the intersection are both stored, and the process proceeds to step (6).
[0008]
If it is not negative in step (3), it is determined that the sound ray is directed from the back side (outside in the case of the simulation of the interior of the building) to the surface, and in fact, the sound ray is directed from that side to the surface. Therefore, the process proceeds to step (6) without setting the surface as a reflection surface candidate. Even if they do not intersect in step (4), the surface cannot be a reflection surface, so that the process proceeds to step (6) without setting it as a candidate.
In step (6), it is determined whether or not all the surfaces have been checked. If not, the process proceeds to step (8) to select the next surface, and returns to step (2) to repeat the same processing. If all the surfaces have been examined, the process proceeds to step (7), and the surface having the shortest path to the intersection among the stored reflection surface candidates is adopted as the reflection surface, and the position of the reflection point and the path length are stored. Then, the traveling direction of the sound ray after the reflection is obtained, and the processing ends.
By this processing, the surface where the sound ray is reflected next and the position of the reflection point are obtained, and based on this, the starting point of the sound ray after reflection and the sound ray vector are obtained. By repeating the processing, the reflection path of the sound ray can be tracked.
[0009]
The above processing will be specifically described with reference to FIGS. First, the inner product of two vectors will be described with reference to FIG.
The inner product (also called a scalar product) of the vector A and the vector B shown in FIG. 8 is | A | · | B | cos θ when the angle between the two vectors A and B is θ. This is referred to as (ABcos θ) for convenience. The sign of this value is determined by cos θ.
Since 0 ≦ θ ≦ π (180 °), the sign of the inner product (ABcosθ) is as follows.
0 ≦ θ <π / 2, that is, positive if θ is an acute angle
θ = π / 2, that is, 0 if θ is a right angle
π / 2 <θ ≦ π That is, if θ is an obtuse angle, negative
The vector A and the angle θ when the angle θ is an obtuse angle are shown by the vector A ′ and the angle θ ′ of a virtual line in FIG.
[0010]
Therefore, as shown in FIG. 9, one
Also, as shown in FIG. 9, the intersection of the
[0011]
Further, as a conventional technique of this kind, there is a method of reproducing virtual three-dimensional spatial sound described in Japanese Patent Application Laid-Open No. 8-272380. According to this method, a plurality of sound (sound ray) vectors are set, a virtual space is constituted by polygonal boundaries, and within a predetermined duration, a propagation history of acoustic vectors in which sound waves are reflected and propagated by the boundaries. Calculates and stores data, and based on the propagation history data, for each acoustic vector, the transient response that reaches the observation point from the sound potential reflected on the boundary and the velocity potential at the observation point created by the micro-plane element of the sound wave Is added to and stored in the time-series numerical array.
Also in this case, it is necessary to calculate the inner product of the normal vector and sound ray vector for each boundary surface to determine whether the boundary contributes to the reflection of sound waves (whether the interface is front or back). Is described in paragraph 0032 of the above publication.
[0012]
As described above, in the conventional acoustic simulation, in order to determine whether a certain boundary surface is a surface contributing to sound reflection (whether the boundary surface is front or back), the normal vector and the normal vector are determined for each boundary surface. An inner product operation with the sound ray vector was required. For this reason, an enormous amount of calculation is required every time one reflection is obtained, which is one of the causes that hinders the speeding up of the simulation.
The present invention has been made in order to solve such a conventional problem, and analyzes a reflection path of a sound wave emitted from a sound source in a predetermined space by calculating a reflection direction in the space. An object of the present invention is to enable geometric acoustic simulation to be performed at high speed with a reduced amount of calculation.
[0013]
[Means for Solving the Problems]
A geometric acoustic simulation apparatus according to the present invention is a geometric acoustic simulation apparatus that analyzes a reflection path of a sound wave emitted from a sound source in a predetermined space by calculating a reflection direction in the space by a reflection direction calculation unit. In order to achieve the above object, the following configuration is provided.
That is, normal vector storage means for storing a normal vector related to each surface involved in the reflection of the sound wave in the space or at the boundary of the space, and each of a plurality of reference vectors serving as a preset reference and the normal Surface classifying means is provided for obtaining an inner product with a vector and classifying (clustering) surfaces involved in the reflection of the sound wave into a plurality of groups according to a combination of signs of the inner product.
[0014]
When calculating the reflection direction of the sound wave, the reflection direction calculation means obtains an inner product of the sound ray vector of the sound wave and each of the plurality of reference vectors, and calculates a sign of the inner product of the plurality of groups. The reflection direction is calculated without considering the surfaces classified into a predetermined group to be excluded from the candidates determined according to the combination. Thus, when determining the presence or absence of reflection at the time of sound ray tracing, it is possible to reduce the number of surfaces for which the presence or absence of reflection is determined for each group (cluster).
The predetermined group may be a group in which an inner product of the sound ray vector and a normal vector of a surface to be classified has a predetermined code.
[0015]
A geometric acoustic simulation method according to the present invention is a geometric acoustic simulation method for analyzing a reflection path by calculating a reflection direction of a sound wave emitted from a sound source in a predetermined space in the space, and achieves the above object. To do so, a simulation is performed as follows.
The normal vector of the surface involved in the reflection of the sound wave in the space or at the boundary of the space is stored, and the inner product of each of the plurality of reference vectors serving as a preset reference and the normal vector is obtained. The surfaces involved in the reflection of the sound wave are classified (clustered) into a plurality of groups according to the combination of the signs of the inner products.
[0016]
When calculating the reflection direction of the sound wave, an inner product of the sound ray vector of the sound wave and each of the plurality of reference vectors is obtained, and is determined according to a combination of signs of the inner product among the plurality of groups. It is characterized in that the reflection direction is calculated without considering the surface classified into a predetermined group to be excluded from candidates.
The predetermined group may be a group in which an inner product of the sound ray vector and a normal vector of a surface to be classified has a predetermined code.
Thus, when determining the presence or absence of reflection during sound ray tracing, the number of faces to be determined can be reduced for each group (cluster).
[0017]
The program according to the present invention further comprises: a computer that controls a geometric acoustic simulation apparatus that analyzes a reflection path of a sound wave; a reflection direction calculation unit that calculates a reflection direction of a sound wave emitted from a sound source in a predetermined space in the space; Normal vector storage means for storing a normal vector for each surface involved in the reflection of the sound wave in the space or at the boundary of the space, and a plurality of reference vectors each serving as a preset reference and the normal vector Is a program that functions as a surface classification unit that classifies the surface involved in the reflection of the sound wave into a plurality of groups according to a combination of signs of the inner product with the inner product, and the function of the reflection direction calculation unit is When calculating the reflection direction of the sound wave, the inner product of the sound ray vector of the sound wave and each of the plurality of reference vectors is calculated. , The surface that are classified into a given group should be excluded from the candidates determined according to the combination of the sign of the inner product of the plurality of groups is a function of calculating a reflection direction without considering.
[0018]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, preferred embodiments of the present invention will be described with reference to the drawings.
First, an embodiment of a geometric acoustic simulation apparatus according to the present invention will be described with reference to FIG. FIG. 1 is a block diagram showing the configuration of the geometric acoustic simulation device.
In the geometric acoustic simulation apparatus shown in FIG. 1, a CPU 11, a
The CPU 11 is a calculation control unit that performs control of the entire apparatus and calculation of simulation, and reads and executes a control / calculation program stored in the
[0019]
The
The surface
As the surface information, a position and a normal vector are input, and a curved surface is approximated by a large number of planes and is represented by each position and a normal vector. Usually, the position of a surface is represented by an equation based on coordinates. When the reflectance is taken into consideration, the information may be input together, but here, the reflectance is not taken into account for the sake of simplicity. As an input method, a method of directly inputting by a user using an input means such as a keyboard, or a method of reading a file created in a predetermined format can be considered.
[0020]
The
The
The geometric acoustic simulation method according to the present invention can be executed by such a geometric acoustic simulation apparatus. However, the geometric acoustic simulation apparatus according to the present invention can also be configured by installing a necessary program on a general-purpose personal computer. .
[0021]
Next, an operation of the above-described geometric acoustic simulation apparatus and an embodiment of a geometric acoustic simulation method according to the present invention, which is performed using the apparatus, will be described with reference to FIGS.
FIG. 2 is a diagram for explaining a reference vector and a classification of a vector based on the reference vector in the geometric acoustic simulation method. FIG. 3 is a flowchart showing a process of classifying a surface normal vector in the above-described geometric acoustic simulation device. 4 is a flowchart showing the processing of the geometric acoustic simulation, and FIG. 5 is a flowchart showing the details of the reflection path tracking processing of the sound ray.
[0022]
A feature of the geometric acoustic simulation apparatus and the geometric acoustic simulation method according to the present invention is that a plurality of non-zero normal vectors (V1, V2, V3,..., Vn) of each surface (n) involved in sound reflection are previously determined. Is classified (clustered) so that a search can be performed in consideration of only planes of a predetermined group when searching for a sound ray path. First, this point will be described. I do. However, the actual simulation is performed on the assumption of a three-dimensional space. Here, for simplicity, the description will be made using a two-dimensional vector.
[0023]
For example, as shown in FIG. 2, reference vectors n1, n2, n3, and n4 serving as references are defined at 45 ° intervals. The problem at this time is the relationship between the reference vectors n1, n2, n3, and n4, and the absolute direction may be determined as convenient.
When the inner product of each of the reference vectors n1, n2, n3, and n4 and the normal vector (V1, V2, V3,..., Vn) of each surface is calculated, each surface is formed by the reference vector and the normal vector. Depending on the combination of the signs (+,-, 0) of the inner product determined by the angles, they can be classified into 16 groups as shown in Table 1.
[0024]
[Table 1]
The directions of the normal vectors of the surfaces classified into these groups are the directions numbered from 1 to 16 around the origin O of the XY coordinates shown in FIG.
The non-zero sound ray vector (x) indicating the traveling direction of the sound ray to be searched for is normally calculated as the inner product of all the normal vectors (V1, V2, V3,..., Vn). However, in the present embodiment, the inner product of each of the reference vectors n1, n2, n3, and n4 is calculated, and they can be classified into the same 16 groups by the combination of the signs. Then, depending on the relationship between the group to which the sound ray vector belongs and the group to which the surface belongs, the sign of the inner product of the sound ray vector and the normal vector of the surface may be known without actually performing the calculation of the inner product. .
For example, the normal vector v of a certain surface i i Is classified into
[0026]
(Equation 1)
On the other hand, when a certain sound ray vector x is classified into
Then, the normal vector v i It can be understood that the inner product of the sound ray vector x and the inner product of the sound ray vector x is positive, without needing to specifically calculate a numerical value.
[0028]
(Equation 2)
[Equation 3]
What is important here is the relationship between the directions of the reference vectors n1, n2, n3, and n4. Even if the absolute directions are different, similar codes can be estimated by similar classification.
Although the two-dimensional case has been described above, a similar method can be extended to three-dimensional cases. Each surface is classified into a plurality of groups by combining the signs of the inner products of each of the reference vectors and the normal vectors. Depending on the combination of the sign of the inner product of each of the sound ray vector and each of the reference vectors, for some groups of surfaces, the sign of the inner product with its normal vector can be calculated without actually calculating the inner product. Can be determined.
[0031]
Therefore, when the inner product is required to be 0 or more in the course of the simulation, the faces of the group in which the inner product is known in advance to be negative can be excluded from the candidates without actually calculating the inner product. . When the inner product is required to be positive, the group whose inner product is known to be 0 or negative, and when the inner product is required to be negative, the group known to be 0 or positive is actually It can be excluded from the candidates without calculating the inner product. However, if the sign of the inner product is not known, it cannot be excluded from the candidates.
As a result, the amount of calculation of the inner product during the execution of the simulation can be greatly reduced.
[0032]
Next, the operation of the above-described geometric acoustic simulation apparatus and the geometric acoustic simulation method of the present invention will be described.
Before executing the simulation, the geometric acoustic simulation apparatus described with reference to FIG. 1 reads the surface information (including the normal vector information as described above) of the surface involved in sound reflection from the surface
Then, when there is an instruction to start the pre-processing from the user thereafter, the CPU 11 starts the pre-processing shown in the flowchart of FIG.
[0033]
First, one of the surfaces is selected as the first surface in step S1, and an inner product of the normal vector of the selected surface and each of a plurality of predetermined reference vectors is calculated in step S2. Then, the process proceeds to step S3, where the surfaces are classified according to the combination of the signs of the inner products obtained, and the result is stored in the RAM. Here, the CPU 11 functions as a normal vector classification means. As the normal vector of the surface, a vector in the direction of the side to which the sound ray can reach (front side: in the case of simulation inside the building) is used.
Thereafter, the process proceeds to step S4, where it is determined whether or not all the surfaces have been checked. If not, the process proceeds to step S5, the next surface is selected, the process returns to step S2, and the process is repeated. If it is, the process ends.
[0034]
By this preprocessing, all surfaces involved in sound reflection can be classified (clustered) according to the direction of the normal vector.
Thereafter, when the user instructs to start the simulation, the CPU 11 starts the processing shown in the flowchart of FIG.
First, in step S11, simulation conditions such as the position of a sound source, the position of a sound collection point, the size of a sound collection area, the cut angle of a sound ray, and the tracking time of a sound ray are set.
Then, proceeding to step S12, the first sound ray is selected, and in step S13, a sound ray vector indicating the traveling direction of the sound ray is set, and the start point is set to the position of the sound source, and the sound ray is fired. Do.
Then, the process proceeds to step S14 to perform a one-step reflection path tracking process on the sound ray. The processing of this subroutine is processing for obtaining the position where the sound ray is reflected next and the traveling direction of the sound ray after reflection, and FIG. 5 is a flowchart showing the processing contents.
[0035]
At the time of this sound ray reflection path tracking processing, the CPU 11 functions as a reflection direction calculation means, and first calculates the inner product of the sound ray vector and each of the reference vectors in step S21 of FIG. Then, proceeding to step S22, the sound ray vectors are classified according to the combination of the signs of the inner products obtained.
Then, the process proceeds to step S23 to limit the search target surface according to the classification result. That is, in this embodiment, since the normal vector is directed to the front side of the surface, only the surface where the inner product of the normal vector and the sound ray vector is negative is the sound ray to be tracked. Is a surface that can be reached from the front side. Therefore, there is no problem even if a surface classified into a group whose inner product is known to be positive or zero in relation to the group in which the sound ray vector is classified is removed from reflection surface candidates in advance. It should be noted that the plane whose sign of the inner product is unknown cannot be excluded from the candidates.
[0036]
Next, the surface to be searched first is selected from the remaining surfaces in step S24. Then, for the surface selected in step S25, the inner product of the normal vector and the sound ray vector is calculated.
Then, in step S26, it is determined whether or not the inner product is negative. If the inner product is negative, it is determined that the sound ray is directed from the front side to the surface, and the process advances to step S27 to determine whether the sound ray and the surface intersect. . If they do, the process proceeds to step S28, and that surface is stored in the
[0037]
If it is not negative in step S26, it is determined that the sound ray is directed from the back side (outside in the case of the simulation of the inside of the building) to the surface, and the sound ray is actually moved from the side to the surface. Since it is not possible to head, the process proceeds to step S29 without setting the surface as a reflection surface candidate. Even if they do not intersect in step S27, the surface cannot be a reflection surface, so the process proceeds to step S29 without setting it as a candidate.
In step S29, it is determined whether or not all the surfaces limited in step S23 have been checked. If not, the process proceeds to step S31 to select the next surface, and returns to step S25 to repeat the processing. If it has been checked, the process proceeds to step S30, and among the stored reflection surface candidates, the surface having the shortest path to the intersection is adopted as the reflection surface, and the position and path length of the reflection point are stored. The process returns to step S15 in FIG. 4 by determining the traveling direction of the sound ray.
The sound ray vector x 'indicating the traveling direction of the sound ray after reflection is represented by x as the original sound ray vector and unit vector v as the normal vector of the reflection surface. i Then, it can be obtained by
[0038]
(Equation 4)
In the process of the flowchart of FIG. 5, when searching for a surface classified into a group in which the inner product of the normal vector and the sound ray vector is known to be negative, step S26 is omitted, and step S25 is omitted. Alternatively, the inner product may be calculated after the surface is adopted as the reflection surface in step S30. In this way, the amount of calculation can be further reduced.
Returning to step S15 in FIG. 4, it is determined whether or not the sound ray has arrived at the obtained sound collection area. If so, the route and arrival time are stored in the
[0040]
In step S18, it is determined whether or not the tracking of all sound rays has been completed, and if completed, the process is terminated. If not, the process proceeds to step S19 to select the next sound ray, and returns to step S13 to repeat the processing.
By performing a simulation by such processing, the reflection surface can be obtained without actually calculating the inner product of the normal vector and the sound ray vector for a part of the surface involved in sound wave reflection when searching for the sound ray path. Can be eliminated from the candidates, and the calculation can be speeded up.
[0041]
【The invention's effect】
As described above, if the acoustic simulation of a predetermined space such as a sound field of a building is performed by the geometric acoustic simulation device and the geometric acoustic simulation method according to the present invention, the acoustic acoustic path is involved in sound wave reflection when searching for a sound ray path. Since a part of the surface can be excluded from reflection surface candidates without actually calculating the inner product of its normal vector and sound ray vector, the calculation for obtaining path information of many sound rays can be performed at high speed. Can be
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of a geometric acoustic simulation apparatus according to an embodiment of the present invention.
FIG. 2 is a diagram for explaining directions of a plurality of reference vectors and normal vectors in the geometric acoustic simulation method of the present invention.
FIG. 3 is a flowchart showing classification processing by a normal vector of a surface involved in reflection in the geometric acoustic simulation apparatus of the present invention.
FIG. 4 is a flowchart showing processing of a geometric acoustic simulation.
FIG. 5 is a flowchart showing details of the sound ray reflection path tracking processing.
FIG. 6 is a diagram for describing a conventional sound ray tracking method.
FIG. 7 is a flowchart showing a part of a conventional sound ray reflection path tracking process.
FIG. 8 is an explanatory diagram of an inner product of two vectors.
FIG. 9 is an explanatory diagram showing a relationship between a selected plane, a normal vector thereof, and a sound ray vector of a sound ray directed to the plane.
[Explanation of symbols]
11: CPU, 12: Program memory, 13: Surface information input means, 14:
Claims (4)
前記空間内又は該空間の境界において前記音波の反射に関与する各面に関する法線ベクトルを記憶する法線ベクトル記憶手段と、
予め設定した基準となる複数の基準ベクトルのそれぞれと前記法線ベクトルとの内積を求めて該内積の符号の組み合わせに応じて前記音波の反射に関与する面を複数のグループに分類する面分類手段とを設け、
前記反射方向演算手段は、前記音波の反射方向を演算する際に、該音波の音線ベクトルと前記複数の基準ベクトルのそれぞれとの内積を求め、前記複数のグループのうち該内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算する手段であることを特徴とする幾何音響シミュレーション装置。In a geometric acoustic simulation apparatus for analyzing a reflection path of the sound wave by calculating a reflection direction of a sound wave emitted from a sound source in a predetermined space in the space by a reflection direction calculation unit,
Normal vector storage means for storing a normal vector for each surface involved in the reflection of the sound wave in the space or at the boundary of the space,
Surface classifying means for obtaining an inner product of each of a plurality of reference vectors serving as a reference set in advance and the normal vector, and classifying surfaces involved in the reflection of the sound wave into a plurality of groups according to a combination of signs of the inner product. And
When calculating the reflection direction of the sound wave, the reflection direction calculation means obtains an inner product of the sound ray vector of the sound wave and each of the plurality of reference vectors, and a combination of signs of the inner product of the plurality of groups. A geometrical acoustic simulation apparatus for calculating a reflection direction without considering a surface classified into a predetermined group to be excluded from candidates determined according to the following.
前記所定のグループは、前記音線ベクトルと分類される前記面の法線ベクトルとの内積が所定の符号になるようなグループであることを特徴とする幾何音響シミュレーション装置。The geometric acoustic simulation device according to claim 1,
The geometric acoustic simulation apparatus according to claim 1, wherein the predetermined group is a group in which an inner product of the sound ray vector and a normal vector of the surface to be classified has a predetermined code.
前記空間内又は該空間の境界において前記音波の反射に関与する面の法線ベクトルを記憶し、
予め設定した基準となる複数の基準ベクトルのそれぞれと前記法線ベクトルとの内積を求めて、該内積の符号の組み合わせに応じて前記音波の反射に関与する面を複数のグループに分類し、
前記音波の反射方向を演算する際に、該音波の音線ベクトルと前記複数の基準ベクトルのそれぞれとの内積を求め、前記複数のグループのうち該内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算することを特徴とする幾何音響シミュレーション方法。A geometric acoustic simulation method for analyzing a reflection path by calculating a reflection direction of a sound wave emitted from a sound source in a predetermined space in the space,
Store the normal vector of the surface involved in the reflection of the sound wave in the space or at the boundary of the space,
Obtain the inner product of each of the plurality of reference vectors to be a preset reference and the normal vector, classify the surface involved in the reflection of the sound wave into a plurality of groups according to a combination of signs of the inner product,
When calculating the reflection direction of the sound wave, the inner product of each of the sound ray vector of the sound wave and the plurality of reference vectors is obtained, and from among the plurality of groups, a candidate determined according to a combination of signs of the inner product. A geometric acoustic simulation method comprising calculating a reflection direction without considering a surface classified into a predetermined group to be excluded.
所定の空間内の音源から発せられる音波の該空間内における反射方向を演算する反射方向演算手段と、
前記空間内又は該空間の境界において前記音波の反射に関与する各面に関する法線ベクトルを記憶する法線ベクトル記憶手段と、
予め設定した基準となる複数の基準ベクトルのそれぞれと前記法線ベクトルとの内積を求めて該内積の符号の組み合わせに応じて前記音波の反射に関与する面を複数のグループに分類する面分類手段として機能させ、
前記反射方向演算手段の機能は、前記音波の反射方向を演算する際に、該音波の音線ベクトルと前記複数の基準ベクトルのそれぞれとの内積を求め、前記複数のグループのうち該内積の符号の組み合わせに応じて定められる候補から除外すべき所定のグループに分類されている面は考慮せずに反射方向を演算する機能であることを特徴とするプログラム。A computer that controls a geometric acoustic simulation device that analyzes a reflection path of a sound wave,
Reflection direction calculating means for calculating a reflection direction of a sound wave emitted from a sound source in a predetermined space in the space,
Normal vector storage means for storing a normal vector for each surface involved in the reflection of the sound wave in the space or at the boundary of the space,
Surface classifying means for obtaining an inner product of each of a plurality of reference vectors serving as a reference set in advance and the normal vector, and classifying surfaces involved in the reflection of the sound wave into a plurality of groups according to a combination of signs of the inner product. Function as
When calculating the reflection direction of the sound wave, the function of the reflection direction calculation means obtains an inner product of a sound ray vector of the sound wave and each of the plurality of reference vectors, and calculates a sign of the inner product of the plurality of groups. A function of calculating a reflection direction without considering a surface classified into a predetermined group to be excluded from candidates determined according to a combination of the above.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002240085A JP3900041B2 (en) | 2002-08-21 | 2002-08-21 | Geometric acoustic simulation apparatus, geometric acoustic simulation method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002240085A JP3900041B2 (en) | 2002-08-21 | 2002-08-21 | Geometric acoustic simulation apparatus, geometric acoustic simulation method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004077352A true JP2004077352A (en) | 2004-03-11 |
| JP3900041B2 JP3900041B2 (en) | 2007-04-04 |
Family
ID=32022967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002240085A Expired - Fee Related JP3900041B2 (en) | 2002-08-21 | 2002-08-21 | Geometric acoustic simulation apparatus, geometric acoustic simulation method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3900041B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589676A (en) * | 2011-12-21 | 2012-07-18 | 中山大学 | Indoor space partitioning method applied to sound ray tracing |
| WO2023202551A1 (en) * | 2022-04-19 | 2023-10-26 | 北京字跳网络技术有限公司 | Acoustic transmission method and device, and nonvolatile computer readable storage medium |
-
2002
- 2002-08-21 JP JP2002240085A patent/JP3900041B2/en not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102589676A (en) * | 2011-12-21 | 2012-07-18 | 中山大学 | Indoor space partitioning method applied to sound ray tracing |
| CN102589676B (en) * | 2011-12-21 | 2015-08-12 | 中山大学 | A kind of acoustic ray tracing method based on interior space subdivision |
| WO2023202551A1 (en) * | 2022-04-19 | 2023-10-26 | 北京字跳网络技术有限公司 | Acoustic transmission method and device, and nonvolatile computer readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3900041B2 (en) | 2007-04-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7512485B2 (en) | Method for routing multiple paths through polygonal obstacles | |
| US7499053B2 (en) | Real-time precision ray tracing | |
| TWI303039B (en) | Methods, computing apparatus and computing system for ray tracing , and machine-accessible medium for application in ray tracing | |
| Goodrich et al. | Snap rounding line segments efficiently in two and three dimensions | |
| Alton et al. | An ordered upwind method with precomputed stencil and monotone node acceptance for solving static convex Hamilton-Jacobi equations | |
| US7030875B2 (en) | Environmental reasoning using geometric data structure | |
| US20090083007A1 (en) | System, method, and computer program product for determining wall thickness in graphic model | |
| JP2023064082A (en) | Method of constructing three-dimensional map in high-definition map, apparatus, device, and storage medium | |
| US7379710B2 (en) | Radio-wave propagation characteristic forecasting system and its method, and program | |
| US8289332B2 (en) | Apparatus and method for determining intersections | |
| JP3900041B2 (en) | Geometric acoustic simulation apparatus, geometric acoustic simulation method, and program | |
| JP2007514243A (en) | Systems and methods for transparent rendering | |
| CN107567641B (en) | System and method for identifying replicas | |
| JP3482055B2 (en) | High-accuracy sound ray tracking device and high-accuracy sound ray tracking method | |
| Sefi | Pay tracing tools for high frequency electromagnetics simulations | |
| CN114092932A (en) | Method and device for identifying chamfered surface, electronic equipment and storage medium | |
| CN109840945B (en) | Finite element preprocessing method and device | |
| Noichl et al. | Enhancing point cloud semantic segmentation in the data‐scarce domain of industrial plants through synthetic data | |
| JP7435779B2 (en) | Distance measurement error calculation device, distance measurement error calculation method, and program | |
| CA2322071C (en) | Improved trim operation based on enhanced face definition | |
| Jankauskas et al. | Measuring Error between Similar Surfaces | |
| Kovalcikova et al. | Method of Detecting Orientation of Red Blood Cells Based on Video Data | |
| JP4157856B2 (en) | Acoustic reflection path discrimination method, computer program, acoustic reflection path discrimination apparatus, and acoustic simulation apparatus | |
| JP2003233631A (en) | Model transformation program and model transformation method | |
| JPH10207927A (en) | Method and device for calculating area and face retrieving device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050224 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061122 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061205 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061218 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313532 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120112 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130112 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140112 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |