JP2010181824A - Sound simulation apparatus and sound simulation program - Google Patents
Sound simulation apparatus and sound simulation program Download PDFInfo
- Publication number
- JP2010181824A JP2010181824A JP2009027622A JP2009027622A JP2010181824A JP 2010181824 A JP2010181824 A JP 2010181824A JP 2009027622 A JP2009027622 A JP 2009027622A JP 2009027622 A JP2009027622 A JP 2009027622A JP 2010181824 A JP2010181824 A JP 2010181824A
- Authority
- JP
- Japan
- Prior art keywords
- sound
- pressure level
- sound pressure
- room
- calculating
- 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.)
- Pending
Links
Images
Landscapes
- Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
Description
本発明は、室内に設けられる音源から発せられる音が、所定の受音点においてどのような音圧レベルとなるかを演算によって算出する音響シミュレーション装置および音響シミュレーションプログラムに関する。 The present invention relates to an acoustic simulation apparatus and an acoustic simulation program for calculating, by calculation, what kind of sound pressure level a sound emitted from a sound source provided in a room has at a predetermined sound receiving point.
従来、ホールなどの室内にスピーカなどの音源を設置する場合、特定の受音点での音圧レベルや周波数特性などを解析するため、コンピュータを用いた音響シミュレーションが行われている(例えば特許文献1)。このような音響シミュレーションにおける演算手法のひとつとして、音源から受音点に向かって直接伝播する直接音のみを考慮して音圧レベルを算出するものがある。しかし、音源が室内に設置される場合、実際には壁面などで1回以上反射する反射音も受音点に到達するため、直接音のみによるシミュレーションでは、実際の音圧レベルよりも小さなレベルとなって算出される傾向があり、正確な音圧レベルを求めることができない。 Conventionally, when a sound source such as a speaker is installed in a room such as a hall, an acoustic simulation using a computer has been performed in order to analyze a sound pressure level or a frequency characteristic at a specific sound receiving point (for example, Patent Documents). 1). As one of calculation methods in such acoustic simulation, there is a method for calculating a sound pressure level in consideration of only a direct sound propagating directly from a sound source toward a sound receiving point. However, when the sound source is installed indoors, the reflected sound that is reflected once or more on the wall surface or the like actually reaches the sound receiving point. Therefore, in the simulation using only direct sound, the sound level is lower than the actual sound pressure level. Therefore, an accurate sound pressure level cannot be obtained.
一方、音響シミュレーションにおける他の演算手法として、直接音だけでなく、壁面などで1回以上反射する反射音も考慮した演算を行う幾何音響シミュレーションがある。幾何音響シミュレーションでは、音源が設置される室内の三次元形状データと、各壁面などの周波数ごとの吸音率とを予め入力することにより、室内の各壁面や天井、床面などで少なくとも1回反射して受音点に到達する全ての反射音を算出してシミュレーションを行うので、より正確な音圧レベルを求めることができる。 On the other hand, as another calculation method in the acoustic simulation, there is a geometric acoustic simulation that performs calculation in consideration of not only the direct sound but also the reflected sound that is reflected at least once by a wall surface or the like. In geometric acoustic simulation, the three-dimensional shape data of the room where the sound source is installed and the sound absorption coefficient for each frequency such as each wall surface are input in advance, so that reflection is performed at least once on each wall surface, ceiling, floor surface, etc. Since all the reflected sounds reaching the sound receiving point are calculated and the simulation is performed, a more accurate sound pressure level can be obtained.
しかしながら、幾何音響シミュレーションの場合、室内の三次元形状データを入力する際、各壁面の位置や形状などを正確に入力していく必要があるため、この入力作業に多大な工数がかかる。加えて、各壁面や天井、床面などの周波数ごとの吸音率も同時に入力していく必要があるため、シミュレーションを開始するまでの入力作業に長期間を要することもある。例えばこのような入力作業に数日〜数週間程度を要することもあり、幾何音響シミュレーションを行う場合の大きな問題点となっている。 However, in the case of geometric acoustic simulation, when inputting indoor three-dimensional shape data, it is necessary to accurately input the position, shape, etc. of each wall surface. In addition, since it is necessary to simultaneously input the sound absorption coefficient for each frequency such as each wall surface, ceiling, floor surface, etc., it may take a long time to input until the simulation is started. For example, such an input operation may take several days to several weeks, which is a serious problem when performing geometric acoustic simulation.
また幾何音響シミュレーションにおける演算処理では、上記のようにして入力された三次元形状データと各壁面における吸音率とに基づいて一つ一つの反射音を算出していき、最終的には受音点に到達する全ての反射音が算出されることになるが、この場合の演算量が膨大であるため、直接音と全ての反射音とを考慮した受音点での音圧レベルを得るまでに多大な時間を要するという問題がある。 In the calculation process in the geometric acoustic simulation, each reflected sound is calculated based on the three-dimensional shape data input as described above and the sound absorption coefficient on each wall surface, and finally the sound receiving point. However, since the amount of computation in this case is enormous, it is necessary to obtain the sound pressure level at the sound receiving point considering the direct sound and all the reflected sounds. There is a problem that it takes a lot of time.
本発明は、上記従来の問題点を解決することを目的としてなされたものであり、直接音のみによるシミュレーションよりも高精度に音圧レベルを算出できるようにし、更には入力作業に要する手間を軽減しつつ、音圧レベルを算出するまでの演算処理を効率的に行うことができるようにした音響シミュレーション装置および音響シミュレーションプログラムを提供するものである。 The present invention has been made for the purpose of solving the above-described conventional problems, and enables calculation of a sound pressure level with higher accuracy than a simulation using only direct sound, and further reduces labor required for input work. However, it is an object to provide an acoustic simulation apparatus and an acoustic simulation program that can efficiently perform arithmetic processing until a sound pressure level is calculated.
上記目的を達成するため、請求項1にかかる発明は、室内に設けられる音源から音が発せられる場合の前記室内における所定の受音点での音圧レベルを演算によって算出する音響シミュレーション装置であって、前記室内の総表面積と、前記室内における平均吸音率との入力を受け付ける入力受付手段と、前記音源から前記受音点に向かって直接伝播する直接音の音圧レベルを算出する第1の音圧レベル算出手段と、前記入力受付手段によって受け付けられた前記総表面積と前記平均吸音率とに基づいて、前記音源から前記室内に拡散する拡散音の音圧レベルを算出する第2の音圧レベル算出手段と、前記第1の音圧レベル算出手段で算出される直接音の音圧レベルと、前記第2の音圧レベル算出手段で算出される拡散音の音圧レベルとに基づいて前記受音点での音圧レベルを算出する第3の音圧レベル算出手段とを備えることを特徴としている。
In order to achieve the above object, an invention according to
かかる発明によれば、直接音だけでなく、拡散音を考慮して受音点における音圧レベルが算出されるため、直接音のみによるシミュレーションよりも高精度に音圧レベルを算出できるようになる。また拡散音を算出するための情報として、室内の総表面積と前記室内における平均吸音率とを入力する構成であるので、従来の幾何音響シミュレーションの場合のように複雑な三次元形状データや壁面ごとの吸音率などは入力する必要がなく、入力作業に要する手間を軽減することができる。さらに演算処理に要する時間も短縮されるので、効率的にシミュレーションを行うことができるようになる。 According to this invention, since the sound pressure level at the sound receiving point is calculated in consideration of not only the direct sound but also the diffused sound, the sound pressure level can be calculated with higher accuracy than the simulation using only the direct sound. . In addition, as the information for calculating the diffused sound, it is configured to input the total surface area of the room and the average sound absorption coefficient in the room. Therefore, it is not necessary to input the sound absorption coefficient, and the labor required for input work can be reduced. Furthermore, since the time required for the arithmetic processing is shortened, simulation can be performed efficiently.
また請求項2にかかる発明は、請求項1記載の音響シミュレーション装置において、前記平均吸音率に関する複数種類のデータを予め記憶しておく記憶手段をさらに備え、前記入力受付手段は、前記記憶手段に記憶された前記複数種類のデータの中から一のデータの選択を受け付けることを特徴としている。
The invention according to
かかる発明によれば、拡散音を算出するための平均吸音率を入力する際、記憶手段に記憶された複数種類のデータの中から一のデータを選択すれば良いので、入力作業をより一層効率的に行うことができる。 According to this invention, when inputting the average sound absorption rate for calculating the diffused sound, it is only necessary to select one data from a plurality of types of data stored in the storage means, so that the input work is made more efficient. Can be done automatically.
また請求項3にかかる発明は、請求項1又は2記載の音響シミュレーション装置において、前記第3の音圧レベル算出手段で算出される前記受音点での音圧レベルに基づいて、前記音源を制御するための制御パラメータを決定し、前記受音点での音圧レベルが所定レベルとなるように制御する制御手段をさらに備えることを特徴としている。 According to a third aspect of the present invention, in the acoustic simulation apparatus according to the first or second aspect, the sound source is determined based on a sound pressure level at the sound receiving point calculated by the third sound pressure level calculating means. Control means for determining a control parameter for control and controlling the sound pressure level at the sound receiving point to be a predetermined level is further provided.
かかる発明によれば、シミュレーション結果に基づく制御を迅速に、かつ高精度に行うことができるようになる。尚、この場合の制御には、音圧レベルのオーバオール値を制御する場合の他、周波数ごとの音圧レベル(周波数特性)を制御する場合も含まれる。 According to this invention, the control based on the simulation result can be performed quickly and with high accuracy. The control in this case includes not only controlling the overall value of the sound pressure level but also controlling the sound pressure level (frequency characteristic) for each frequency.
また請求項4にかかる発明は、請求項3記載の音響シミュレーション装置において、前記音源は、音の指向性を制御可能なアレイスピーカであり、前記第1および第2の音圧レベル算出手段は、前記アレイスピーカの指向性を制御するための複数の制御パターンに基づいて音圧レベルを算出し、前記第3の音圧レベル算出手段は、前記複数の制御パターンのそれぞれについて前記受音点での音圧レベルを算出し、前記制御手段は、前記複数の制御パターンのそれぞれを適用して算出された前記受音点での音圧レベルが略同一レベルとなるように制御することを特徴としている。 According to a fourth aspect of the present invention, in the acoustic simulation apparatus according to the third aspect, the sound source is an array speaker capable of controlling sound directivity, and the first and second sound pressure level calculating means include: The sound pressure level is calculated based on a plurality of control patterns for controlling the directivity of the array speaker, and the third sound pressure level calculating means is configured to calculate the sound pressure level at each sound receiving point for each of the plurality of control patterns. A sound pressure level is calculated, and the control means controls the sound pressure level at the sound receiving point calculated by applying each of the plurality of control patterns to be substantially the same level. .
かかる発明によれば、アレイスピーカの指向性を制御した場合でも音圧レベルの変化を抑制することができるようになる。 According to this invention, even when the directivity of the array speaker is controlled, a change in the sound pressure level can be suppressed.
また請求項5にかかる発明は、情報処理装置によって実行され、前記情報処理装置に、室内に設けられる音源から音が発せられる場合の前記室内における所定の受音点での音圧レベルを演算する処理を行わせる音響シミュレーションプログラムであって、前記情報処理装置に、前記室内の総表面積と、前記室内における平均吸音率との入力を受け付ける入力受付ステップと、前記音源から前記受音点に向かって直接伝播する直接音の音圧レベルを算出する第1の音圧レベル算出ステップと、前記入力受付ステップにおいて受け付けられた前記総表面積と前記平均吸音率とに基づいて、前記音源から前記室内に拡散する拡散音の音圧レベルを算出する第2の音圧レベル算出ステップと、前記第1の音圧レベル算出ステップで算出される直接音の音圧レベルと、前記第2の音圧レベル算出ステップで算出される拡散音の音圧レベルとを加算することにより前記受音点での音圧レベルを算出する第3の音圧レベル算出ステップとを実行させることを特徴としている。 The invention according to claim 5 is executed by an information processing device, and calculates a sound pressure level at a predetermined sound receiving point in the room when sound is emitted from a sound source provided in the room to the information processing device. An acoustic simulation program for performing processing, wherein the information processing device receives an input of a total surface area of the room and an average sound absorption rate in the room, and from the sound source toward the sound receiving point. Diffusion from the sound source into the room based on the first sound pressure level calculating step for calculating the sound pressure level of the direct sound directly propagating and the total surface area and the average sound absorption rate received in the input receiving step The second sound pressure level calculating step for calculating the sound pressure level of the diffuse sound to be performed, and the direct sound calculated in the first sound pressure level calculating step. A third sound pressure level calculating step for calculating the sound pressure level at the sound receiving point by adding the pressure level and the sound pressure level of the diffused sound calculated in the second sound pressure level calculating step; It is characterized by executing.
本発明によれば、直接音のみをシミュレーションする場合よりも、受音点での音圧レベルを高精度に算出できるようになり、しかもシミュレーションを行うための条件を入力する入力作業は従来の幾何音響シミュレーションの場合よりも軽減されるようになる。また演算時間も短縮されるので、効率的に音響シミュレーションを行うことが可能である。 According to the present invention, the sound pressure level at the sound receiving point can be calculated with higher accuracy than when only the direct sound is simulated, and the input operation for inputting the conditions for performing the simulation is the conventional geometry. This is less than in the case of acoustic simulation. In addition, since the calculation time is shortened, it is possible to perform an acoustic simulation efficiently.
以下、本発明の好ましい一実施形態について図面を参照しつつ詳細に説明する。図1は、本実施形態における音響システム1の一構成例を示すブロック図である。この音響システム1は、音声信号を処理する信号処理部10と、信号処理部10によって処理された音声信号に基づく音を発生させるための音源となるスピーカ装置18と、スピーカ装置18から発せられる音の音響特性を演算によって算出する音響シミュレーション装置2とを備える構成である。
Hereinafter, a preferred embodiment of the present invention will be described in detail with reference to the drawings. FIG. 1 is a block diagram illustrating a configuration example of an
図2は、本実施形態におけるスピーカ装置18の一構成例を示す斜視図である。本実施形態においてスピーカ装置18は、筐体18aに対して複数のスピーカユニット181,182,…,18nが組み込まれたアレイスピーカとなっている。このアレイスピーカは、各スピーカユニットに供給する音声信号を個別に制御することで、スピーカ装置18から放出される音の指向性を制御することが可能である。尚、図例では複数のスピーカユニット181,182,…,18nは、縦方向および横方向に配置された二次元配列となっている場合を示しているが、これに限られるものではなく、一次元配列であっても構わない。また、アレイスピーカとしてスピーカ装置18に配置するスピーカユニットの数は任意であり、図例のものに限定されるものではない。そして本実施形態では、スピーカ装置18は例えば図3に示すようにホールなどの室内8に設置された状態で使用される。
FIG. 2 is a perspective view illustrating a configuration example of the
図1に戻り、信号処理部10は、信号入力部11、レベル調整部12、イコライザ13、アレイ制御用フィルタ141,142,…,14n、および、信号出力部151,152,…,15nを備えている。信号入力部11は、音声信号を入力するインタフェースであり、信号入力部11に入力した音声信号は、レベル調整部12に出力される。レベル調整部12は、音声信号の全体的な音圧レベルを調整するためのものである。このレベル調整部12は、音響シミュレーション装置2によって制御される。つまり、レベル調整部12において音圧レベルを調整するためのパラメータは音響シミュレーション装置2によって設定されるようになっている。そしてレベル調整部12でレベル調整が行われた音声信号は、イコライザ13に出力される。イコライザ13は、スピーカ装置18に対して供給する音声信号の周波数特性を調整するためのものである。このイコライザ13もまた音響シミュレーション装置2によって制御される。つまりイコライザ13において調整される周波数帯域ごとのパラメータは音響シミュレーション装置2によって設定されるようになっている。尚、イコライザ13は、PEQ(パラメトリックイコライザ)などのIIR型フィルタによる構成であっても良いし、またFIR型フィルタによる構成であっても良い。そしてイコライザ13から出力される音声信号は、複数のアレイ制御用フィルタ141,142,…,14nに入力する。
1, the
アレイ制御用フィルタ141,142,…,14nは、スピーカ装置18に設けられた複数のスピーカユニット181,182,…,18nに対応して個別に設けられており、それぞれのスピーカユニット181,182,…,18nに供給する音声信号に対して個別にディレイやゲインを付与することでスピーカ装置18から放出される音の指向性を制御するフィルタである。各アレイ制御用フィルタ141,142,…,14nにおいて音声信号に付与するディレイやゲインは、音響シミュレーション装置2によって設定される。尚、アレイ制御用フィルタ141,142,…,14nは、FIR型フィルタフィルタであっても良い。
The array control filters 141, 142,..., 14 n are individually provided corresponding to the plurality of
信号出力部151,152,…,15nもまた、複数のスピーカユニット181,182,…,18nに対応して個別に設けられており、各アレイ制御用フィルタ141,142,…,14nから出力される音声信号をアンプで増幅させた後、それぞれに対応するスピーカユニット181,182,…,18nに出力する。
.., 15n are also provided individually corresponding to the plurality of
一方、音響シミュレーション装置2は、一般的なコンピュータ(PC)である情報処理装置20によって構成される。この情報処理装置20は、各種演算処理を行うCPU21と、キーボードやマウスなどを含む操作部22と、CRTやLCDなどで構成される表示部23と、各種データやプログラムなどを記憶しておくハードディスク装置24と、CPU21が演算処理を行う際に一時的なデータなどを格納するためのメモリ25とを備えている。ハードディスク装置24には、情報処理装置20を音響シミュレーション装置2として機能させるための音響シミュレーションプログラム25が予めインストールされている。またハードディスク装置24には、音響シミュレーションプログラム25のインストールに伴って、複数種類の平均吸音率データ26が予め記憶されている。そしてCPU21が音響シミュレーションプログラム25を読み出して実行することにより、情報処理装置20は音響シミュレーション装置2として機能する。
On the other hand, the
音響シミュレーション装置2は、図3に示すように、室内8に設けられたスピーカ装置18から発せられる音が、ユーザによって指定される室内8のある受音点9においてどのような音圧レベルとなるかをシミュレーション演算することによって算出する。特に本実施形態におけるシミュレーションでは、スピーカ装置18から受音点9に対して直接伝播する直接音と、スピーカ装置18から室内空間に拡散する拡散音とを算出し、これら直接音と拡散音とを考慮して受音点9での音圧レベルを算出する。そして音響シミュレーション装置2は、その算出した音圧レベルに基づいてレベル調整部12やイコライザ13に設定するための制御パラメータを決定することにより、受音点における音響特性がユーザの所望する特性となるように制御する。
As shown in FIG. 3, the
また音響シミュレーション装置2は、ユーザによって指定されるスピーカ装置18の制御パターンに基づいてアレイ制御用フィルタ141,142,…,14nのそれぞれにディレイやゲイン付与を設定し、アレイスピーカの指向性を制御する。このとき、音響シミュレーション装置2は、スピーカ装置18の制御パターンが変更される前後において音圧レベルが大きく変化しないように、シミュレーション結果に基づいて音圧レベルを調整することもできる。
The
図4は、情報処理装置20が音響シミュレーション装置2として機能する場合のCPU21における機能構成を示すブロック図である。入力受付部31は、操作部22を介してユーザからの入力操作を受け付けるための処理部であり、シミュレーションを行うための各種条件や、受音点における所望の音圧レベルなどの入力を受け付ける。
FIG. 4 is a block diagram illustrating a functional configuration of the
シミュレーションを行うための条件としては、例えば図3に示すような音響環境の場合、スピーカ装置18の構成や特性、スピーカ装置18と受音点9との位置関係などが入力される。また本実施形態では室内8における拡散音を算出するための条件として、室内8の総表面積と、室内8における平均吸音率とが入力される。ここで室内8の平均吸音率は、周波数ごとに異なる吸音率となるのが一般的である。そのため、室内8の平均吸音率は、周波数ごとに入力される。ところが、室内8における周波数ごとの平均吸音率を、ユーザが常に手動操作で入力しなければならないこととすると、シミュレーションの条件を入力するためのユーザの負担が大きくなる。そのため本実施形態では、予め代表的な複数種類の平均吸音率データ26をハードディスク装置24に格納しており、ユーザが室内8の平均吸音率を設定する際、それら複数種類の平均吸音率データ26の中から一の平均吸音率データを選択することにより、周波数ごとの平均吸音率を自動設定することができるようになっている。
As conditions for performing the simulation, for example, in the case of an acoustic environment as shown in FIG. 3, the configuration and characteristics of the
図5は、ハードディスク装置24に記憶される複数種類の平均吸音率データ26の一例を示す図である。図5に示すように、平均吸音率データ26には、スタジオ、教会、バンケットルーム、コンサートホールなどの代表的な音響空間における平均的な周波数ごとの平均吸音率が予め格納されている。また平均吸音率データ26には、ユーザが予め入力した周波数ごとの平均吸音率データがユーザメモリ1,2,3として格納されている。そのため、ユーザは、音響シミュレーション装置2におけるシミュレーションの条件として平均吸音率を入力する際、これら複数種類の平均吸音率データ26の中から、スピーカ装置18が設置された室内8の音響空間に最も近いデータを選択すれば良い。
FIG. 5 is a diagram illustrating an example of a plurality of types of average sound
図6は、音響シミュレーション装置2において室内8の総表面積と平均吸音率との入力受付を行う際、入力受付部31によって表示部23に表示される表示画面の一例を示す図である。図6に示すように、表示部23には、室内8の総表面積を入力するための面積入力欄23aと、平均吸音率を入力するための平均吸音率入力欄23bとが表示される。面積入力欄23aには、例えばユーザが手計算などで算出した室内8の総表面積をテンキーなどで入力する。また平均吸音率入力欄23bの右側にはプルダウンボタン23cが表示されており、ユーザはマウスなどを操作することでプルダウンボタン23cをクリックすると、平均吸音率入力欄23bの下側に複数種類の平均吸音率データ26がプルダウンメニュー23dとして表示される。したがって、ユーザはこのプルダウンメニュー23dの中から室内8の空間に最も適した平均吸音率データを選択することで、平均吸音率を指定することができる。尚、図例では、室内8の総表面積として「2300」が入力され、平均吸音率として「教会(ライブ)」が選択された場合を示している。
FIG. 6 is a diagram illustrating an example of a display screen displayed on the
また入力受付部31は、上記の他にも、シミュレーションを行うための条件として、例えばレベル調整部12やイコライザ13に設定する制御パラメータや、スピーカ装置18の指向性を制御するための制御パターンなどの各種情報の入力も受け付ける。
In addition to the above, the
そしてシミュレーションを行うために必要な全ての情報が入力されると、CPU21はその入力された情報に基づいてシミュレーションのための演算処理を開始する。このときCPU21は、図4に示すように音圧レベル算出部32として機能する。この音圧レベル算出部32は、第1の音圧レベル算出部33、第2の音圧レベル算出部34および第3の音圧レベル算出部35として更に機能する。第1の音圧レベル算出部33は、音源であるスピーカ装置18からユーザによって指定された受音点9に向かって直接伝播する直接音の音圧レベルを算出するための処理部である。また第2の音圧レベル算出部34は、室内8の総表面積と平均吸音率とに基づいて、スピーカ装置18から室内8に拡散する拡散音の音圧レベルを算出するための処理部である。そして第3の音圧レベル算出部35は、第1の音圧レベル算出部33で算出される直接音の音圧レベルと、第2の音圧レベル算出部34で算出される拡散音の音圧レベルとに基づいて受音点9における音圧レベルを算出する処理部である。本実施形態では、これら第1〜第3の音圧レベル算出部33〜35のそれぞれにおいて順次演算処理が行われることにより、直接音と拡散音との双方を含む受音点9での音圧レベルが算出される。以下、さらに詳しく説明する。
When all the information necessary for performing the simulation is input, the
まず、第1の音圧レベル算出部33における演算処理について説明する。第1の音圧レベル算出部33は、スピーカ装置18に含まれる複数のスピーカユニット181,182,…,18nの一つ一つについて受音点9に伝播する直接音の音圧レベルを算出した後、各スピーカユニットからの音圧レベルを加算することで、スピーカ装置18からの直接音による音圧レベルを算出する。
First, calculation processing in the first sound pressure
図7は、i番目(但し、iは1≦i≦nを満たす任意の数)のスピーカユニットと、受音点9の関係を示す図である。図7に示すように、i番目のスピーカユニットの正面方向であるZ軸を基準とし、受音点9に対する方位角をθiとし、仰角をφiとすると共に、受音点9までの距離(m)をriとした場合、i番目のスピーカユニットから受音点9までの伝達関数は次の数1の式で表される。 FIG. 7 is a diagram showing the relationship between the i-th speaker unit (where i is an arbitrary number satisfying 1 ≦ i ≦ n) and the sound receiving point 9. As shown in FIG. 7, the azimuth angle with respect to the sound receiving point 9 is θi, the elevation angle is φi, and the distance to the sound receiving point 9 (m ) Is ri, the transfer function from the i-th speaker unit to the sound receiving point 9 is expressed by the following equation (1).
ここで、fは周波数(Hz)であり、λは波長(m)である。またDi(f,θi,φi)は、i番目のスピーカユニットの指向係数であり、周波数ごとに異なる係数となる。この指向係数は、スピーカユニット固有の係数であり、方位角θi及び仰角φiの方向に向かうエネルギーが基準となるZ軸方向のエネルギーの何倍になるかを示している。このような指向係数は、スピーカ装置18の構成や特性に基づいて決定される。さらに、m(f)は、空気吸収による音の減衰係数(1/m)であり、周波数ごとに異なる係数となる。
Here, f is a frequency (Hz) and λ is a wavelength (m). Di (f, θi, φi) is a directivity coefficient of the i-th speaker unit, and is a different coefficient for each frequency. This directivity coefficient is a coefficient specific to the speaker unit, and indicates how many times the energy in the direction of the azimuth angle θi and the elevation angle φi becomes the reference energy in the Z-axis direction. Such a directivity coefficient is determined based on the configuration and characteristics of the
また図7に示すように、i番目のスピーカユニットの正面側に1m離れたZ軸上の点での音圧P1iを算出する。この音圧P1iは、スピーカ装置18の構成や特性、並びに、信号処理部10における各部(レベル調整部12、イコライザ13、アレイ制御用フィルタ14iおよび信号出力部15i)の設定状態を反映させた公知の演算を行うことにより算出することができる。そして1m点での音圧P1iを算出すると、数1の式の伝達関数により、i番目のスピーカユニットから受音点9に伝播する直接音の音圧Pi(f)は、次の数2の式で表される。
Further, as shown in FIG. 7, the sound pressure P1i at a point on the Z-axis that is 1 m away from the front side of the i-th speaker unit is calculated. This sound pressure P1i reflects the configuration and characteristics of the
上記数2の式に基づく演算を行うことにより、i番目のスピーカユニットから受音点9に直接伝播する直接音の音圧Pi(f)を算出することができる。そしてスピーカ装置18の全体から発せられる直接音による受音点9での音圧Pdは、上記数2の式に基づいて算出される音圧Pi(f)をスピーカユニットごとに算出し、それらを加算することによって求められるため、次の数3の式で表される。尚、数3の式において、Nは、スピーカ装置18におけるスピーカユニットの数(n)である。
By performing the calculation based on the equation (2), the sound pressure Pi (f) of the direct sound that propagates directly from the i-th speaker unit to the sound receiving point 9 can be calculated. The sound pressure Pd at the sound receiving point 9 due to the direct sound emitted from the
第1の音圧レベル算出部32は、上記演算を行うことにより、スピーカ装置18からの直接音による音圧Pdを算出する。尚、音圧Pdは、周波数fの関数であるため、上記数3の式に周波数fを代入して演算を行うことにより、第1の音圧レベル算出部32は、受音点9における直接音の音圧Pdを周波数ごとに算出する。
The first sound pressure
次に、第2の音圧レベル算出部34における演算処理について説明する。第2の音圧レベル算出部34は、スピーカ装置18から発せられる拡散音による音圧レベルを算出するため、スピーカ装置18を中心に設置した仮想的な球体を想定した演算を行う。図8は、この仮想的な球体を示す図である。第2の音圧レベル算出部34は、図8に示すように、仮想的な球体の表面を均等間隔に設定された複数の経線及び緯線によって細分化し、経線と緯線とが交叉する格子点(グリッド)に対して順番に番号を付与する。そして格子点の数がM個である場合、j番目(但し、jは1≦j≦Mを満たす任意の数)のグリッドに着目し、そのグリッドjがスピーカ装置18から十分に遠いと仮定し、次の数4の式に基づく演算を行うことにより、グリッドjにおける音圧を計算する。
Next, calculation processing in the second sound pressure
ここでP1iは、上述したようにi番目のスピーカユニットの正面側に1m離れたZ軸上の点での音圧P1iである。そのため音圧P1iには、信号処理部10における各部(レベル調整部12、イコライザ13、アレイ制御用フィルタ14iおよび信号出力部15i)の設定状態が反映されている。またグリッドjのスピーカ装置18に対する方位角がθjであり、仰角がφjである。またDi(f,θi,φi)は、上述したようにi番目のスピーカユニットの指向係数である。したがって、数4の式に基づく演算を行うことにより、アレイスピーカとして設けられたN個のスピーカユニットから発せられる音のグリッドjでの音圧が求められる。
Here, P1i is the sound pressure P1i at the point on the Z-axis that is 1 m away from the front side of the i-th speaker unit as described above. Therefore, the sound pressure P1i reflects the setting state of each unit (
そして数4の式で求められるグリッドjでの音圧を、仮想的な球体の全体について積分することにより、スピーカ装置18から全方向に向かう音響出力Wが求められる。すなわち、音響出力Wは、次の数5の式で表される。
Then, the sound output W directed in all directions from the
ここでρは空気密度(kg/m3)であり、cは音速(m/s)である。またP0は、基準音圧=2×10-5(Pa)であり、r∞は十分遠い距離(m)を示している。さらにSjは、グリッドjにおける細分化された面積である。尚、数5の式は周波数fの関数となっており、数5の式に基づく演算を周波数ごとに行うことで、周波数ごとの音響出力Wを算出することができる。 Here, ρ is the air density (kg / m 3), and c is the speed of sound (m / s). Further, P0 is the reference sound pressure = 2 × 10 −5 (Pa), and r∞ indicates a sufficiently far distance (m). Furthermore, Sj is a subdivided area in the grid j. Note that the formula 5 is a function of the frequency f, and the sound output W for each frequency can be calculated by performing the calculation based on the formula 5 for each frequency.
そして数5の式によって求められるスピーカ装置18の音響出力Wが、室内8に対して出力されている場合、ユーザによって入力された室内8の平均吸音率をα、室内8の総表面積をSaとすると、室内8に存在する拡散音のエネルギ密度Esは、次の数6の式によって求められる。
When the sound output W of the
さらに数6の式によって求められる拡散音のエネルギ密度Esに基づいて次の数7の式による演算を行うことにより、室内8における拡散音の音圧を算出することができる。 Furthermore, the sound pressure of the diffused sound in the room 8 can be calculated by performing the calculation according to the following expression (7) based on the energy density Es of the diffused sound obtained by the expression (6).
したがって、第2の音圧レベル算出部34は、上記演算を行うことにより、スピーカ装置18から全方向に出力され、室内8に拡散した状態で存在する拡散音の音圧を算出する。尚、数7の式におけるPsもまた、周波数fの関数であるため、上記数7の式に周波数fを代入して演算を行うことにより、第2の音圧レベル算出部34は、室内8における拡散音の音圧を周波数ごとに算出する。
Therefore, the second sound pressure
次に、第3の音圧レベル算出部35における演算処理について説明する。第3の音圧レベル算出部35は、受音点9における直接音の音圧と、拡散音の音圧とに基づいて、直接音と拡散音とを含む受音点9での音圧レベルLを算出する。具体的には、次の数8の式に基づく演算を行うことによって音圧レベルLを算出する。 Next, calculation processing in the third sound pressure level calculation unit 35 will be described. The third sound pressure level calculation unit 35 is based on the sound pressure of the direct sound at the sound receiving point 9 and the sound pressure of the diffused sound, and the sound pressure level at the sound receiving point 9 including the direct sound and the diffused sound. L is calculated. Specifically, the sound pressure level L is calculated by performing an operation based on the following equation (8).
ここで直接音のみによる音圧Pdは第1の音圧レベル算出部33によって算出されており、また拡散音による音圧Psは第2の音圧レベル算出部34によって算出されている。また基準音圧P0は予め定義された既知の値である。したがって、上記数8の式にそれぞれの値を代入して演算を行うことにより、第3の音圧レベル算出部35は、受音点9における直接音と拡散音の双方を考慮した音圧レベルLを算出することができる。尚、数8の式におけるLもまた、周波数fの関数であるため、上記数8の式に周波数fを代入して演算を行うことにより、第3の音圧レベル算出部35は、受音点9における直接音と拡散音との双方を考慮した音圧レベルLを周波数ごとに算出することができる。したがって、受音点9における周波数特性を算出することができると共に、それら周波数ごとの音圧レベルLに基づいてオーバオール値(OA値)を算出することもできる。
Here, the sound pressure Pd due to only the direct sound is calculated by the first sound pressure
図9は、音圧レベル算出部32において上述した演算が行われることによって算出される受音点9での音圧レベル(周波数特性)の一例を示す図である。図例において、音圧レベル特性T1は上述した演算によって直接音と拡散音との双方をシミュレーションすることによって算出される特性を示している。また音圧レベル特性T2は比較例であり、直接音のみをシミュレーションすることによって算出される特定を示している。受音点9での音圧レベルを実際に測定すると、直接音のみをシミュレーションした音圧レベル特性T2よりも測定値の方が大きなレベルとなり、特に周波数の低い低音域になる程そのレベル差は大きくなる。一方、本実施形態のシミュレーションによって算出される音圧レベル特性T1は、直接音のみをシミュレーションした音圧レベル特性T2よりも実際の測定値により近い特性となっており、低周波域で音圧レベル特性T2よりも大きな音圧レベルが算出されている。したがって、本実施形態におけるシミュレーションは、直接音のみをシミュレーションする場合と比較すれば、受音点9での音圧レベルLを高精度に算出することができる。
FIG. 9 is a diagram illustrating an example of a sound pressure level (frequency characteristic) at the sound receiving point 9 calculated by performing the above-described calculation in the sound pressure
またシミュレーションの演算時間に関し、従来の幾何音響シミュレーションと比較すると、本実施形態のシミュレーションに要する演算時間は短くなる。つまり、従来の幾何音響シミュレーションでは、詳細な三次元形状データと各壁面での吸音率とに基づいて一つ一つの反射音を算出していくことが必要であり、受音点に到達する全ての反射音を算出するまでの演算量は膨大なものとなるが、本実施形態では室内8に対してほぼ均一に拡散する拡散音を算出するようにしているので、壁面での反射角度などを計算しながら受音点に到達する一つ一つの反射音を算出する必要がなく、演算量が少なくなるので効率的な演算が可能である。 Further, regarding the calculation time of the simulation, the calculation time required for the simulation of the present embodiment is shorter than that of the conventional geometric acoustic simulation. In other words, in the conventional geometric acoustic simulation, it is necessary to calculate each reflected sound based on detailed three-dimensional shape data and the sound absorption rate at each wall surface, However, in this embodiment, the diffused sound that diffuses almost uniformly in the room 8 is calculated, so the reflection angle on the wall surface and the like can be calculated. It is not necessary to calculate each reflected sound that reaches the sound receiving point while calculating, and the calculation amount is reduced, so that efficient calculation is possible.
さらに本実施形態のシミュレーションでは、反射音ではなく、拡散音を算出するようにしているので、ユーザは室内8の総表面積と、室内8の平均吸音率とを入力すれば良い。つまり、従来の幾何音響シミュレーションのように室内の詳細な三次元形状データを入力する必要がなくなるので、ユーザの入力作業を軽減することが可能である。 Furthermore, in the simulation according to the present embodiment, not the reflected sound but the diffused sound is calculated. Therefore, the user may input the total surface area of the room 8 and the average sound absorption rate of the room 8. That is, since it is not necessary to input detailed indoor three-dimensional shape data as in the conventional geometric acoustic simulation, it is possible to reduce the user's input work.
上記のようにして受音点9での直接音と拡散音との双方を考慮した音圧レベルLを算出すると、CPU21では図4に示す音圧レベル解析部36が機能する。音圧レベル解析部36は、音圧レベル算出部32において算出された受音点9での音圧レベルLを表示部23に表示すると共に、操作部22を介して入力されるユーザが所望する音圧レベルとの比較を行うなど、算出した音圧レベルLの解析を行う。
When the sound pressure level L is calculated in consideration of both the direct sound and the diffused sound at the sound receiving point 9 as described above, the sound pressure
そして算出した音圧レベルLがユーザの所望する音圧レベルと差がある場合には、制御部37が機能し、信号処理部10を制御する。制御部37は、音圧レベルLとユーザが所望する音圧レベルとの差に基づいてレベル調整部12やイコライザ13を調整するための制御パラメータを決定し、その制御パラメータをレベル調整部12やイコライザ13に設定することで受音点9での音圧レベルLがユーザの所望する音圧レベルとなるように制御する。尚、制御部37が決定する制御パラメータを音圧レベル算出部32に再度入力することにより、音圧レベル算出部32がその決定された制御パラメータに基づいて再度シミュレーション演算を行うようにしても良い。
When the calculated sound pressure level L is different from the sound pressure level desired by the user, the
図10は、CPU21において上述した各機能を実現するための処理手順の一例を示すフローチャートである。この処理は、CPU21が音響シミュレーションプログラム25を読み出して実行することによって開始される。処理が開始すると、音響シミュレーション装置2は、まずシミュレーションを行うための条件の入力を受け付ける(ステップS11)。ここでは上述した演算を行うために必要となる各種の情報が入力される。このとき、表示部23には図6に示した画面が表示され、拡散音を算出するために必要となる室内8の総表面積と、室内8の平均吸音率との入力を受け付ける。またこのとき、ユーザが所望する音圧レベルの入力を受け付けるようにしても良い。
FIG. 10 is a flowchart illustrating an example of a processing procedure for realizing each function described above in the
シミュレーションを行うための条件の入力が終了すると、音響シミュレーション装置2は、直接音による音圧レベルの算出演算を行う(ステップS12)。ここではCPU21において上述した第1の音圧レベル算出部33が機能し、受音点9における直接音の音圧レベルPdが算出される。次に音響シミュレーション装置2は、拡散音による音圧レベルの算出演算を行う(ステップS13)。ここではCPU21において上述した第2の音圧レベル算出部34が機能し、室内8における拡散音の音圧レベルPsが算出される。その後、音響シミュレーション装置2は、直接音と拡散音との双方を合算した受音点9における音圧レベルLの算出演算を行う(ステップS14)。ここではCPU21において上述した第3の音圧レベル算出部34が機能し、室内8における受音点9での直接音と拡散音の双方による音圧レベルLが算出される。
When the input of the conditions for performing the simulation is completed, the
そして音響シミュレーション装置2は、他にシミュレーションを行うための条件があるかどうかを判断する(ステップS15)。例えば、上述したステップS12〜S14において特定の周波数での音圧レベルLが算出されるとした場合、他の周波数での音圧レベルLも算出する必要があればYESとなる。またスピーカ装置18の指向性を制御する制御パターンを変化させてシミュレーションを行う場合もYESとなる。そしてシミュレーションを行うための条件を変更して上述したステップS12〜S14の算出処理を繰り返し実行する。一方、他にシミュレーションを行うための条件がない場合にはステップS15においてNOとなる。
Then, the
そして音響シミュレーション装置2は、ステップS14で算出した受音点での音圧レベルLの解析を行う(ステップS16)。このとき、周波数ごとに算出される音圧レベルLに基づいて解析を行っても良いし、また周波数ごとに算出される音圧レベルLのオーバオール値(OA値)に基づいて解析を行っても良い。いずれに基づいて解析を行うかは、ユーザが設定する解析条件による。例えばアレイスピーカであるスピーカ装置18の制御パターンを変化させることによる音圧レベルLの変化を抑制したい場合には、オーバオール値による解析を行うことが好ましい。また受音点9での周波数特性をユーザが所望する特性となるように調整したい場合には、周波数ごとに算出される音圧レベルLに基づいて周波数特性の解析を行うことが好ましい。
The
図11は、スピーカ装置18の制御パターンをパターン1〜5に変化させてシミュレーションを行った場合に算出されるそれぞれの音圧レベルL1〜L5のオーバオール値を示す図である。スピーカ装置18の制御パターンを変化させると、受音点9での音圧レベルL1〜L5は、図11に示すように異なる値となる。そのため、制御パターンの変化に伴う音圧レベルLの変化を抑制したい場合、音響シミュレーション装置2は、例えば図11に示す複数の音圧レベルL1〜L5の最大値と最小値との間に所定レベル以上の差があるか否かを判定する。
FIG. 11 is a diagram illustrating the overall values of the sound pressure levels L1 to L5 calculated when the simulation is performed by changing the control pattern of the
また周波数ごとの音圧レベルLをユーザが所望する音圧レベルに設定したい場合、音響シミュレーション装置2は、図9に示した音圧レベル特性T1がユーザの所望する周波数特性と差があるか否かを判定する。
When it is desired to set the sound pressure level L for each frequency to a sound pressure level desired by the user, the
そして音響シミュレーション装置2は、音圧レベルLの解析結果に基づき、信号処理部10を制御する必要があるか否かを判断し(ステップS17)、制御する必要がない場合には処理を終了する。また制御する必要がある場合には、次のステップに進み、レベル調整部12とイコライザ13との少なくとも一方を制御するための制御パラメータを決定する(ステップS18)。例えば、図11に示したような制御パターンの変化に伴う音圧レベルL1〜L5の変化を抑制したい場合、音響シミュレーション装置2は、複数の制御パターンのそれぞれを適用した場合の複数の音圧レベルL1〜L5が略同一レベルとなるようにレベル調整部12に対する制御パラメータを決定する。このとき、複数の音圧レベルL1〜L5の平均レベルALで略同一となるように各制御パターン適用時の制御パラメータを決定しても良いし、またユーザが入力する所望の音圧レベルをターゲットレベルTLとし、そのターゲットレベルTLで略同一となるように各制御パターン適用時の制御パラメータを決定しても良い。また例えば、周波数ごとの音圧レベルLをユーザが所望する周波数特性に設定したい場合、音響シミュレーション装置2は、ユーザの所望する周波数特性と差に基づいて周波数ごとに音圧レベルを調整するためのイコライザ13に対する制御パラメータを決定する。尚、レベル調整部12とイコライザ13の双方を制御するための制御パラメータを決定しても良い。
Then, the
そして制御パラメータを決定すると、音響シミュレーション装置2は、その決定した制御パラメータをレベル調整部12及び/又はイコライザ13に対して設定することで音圧レベルの制御を行う(ステップS19)。以上で、この処理は終了するが、制御パラメータを設定した状態で再度シミュレーションを行う場合には、上述した処理を繰り返し実行するような処理手順としても良い。
When the control parameter is determined, the
以上のように本実施形態の音響シミュレーション装置2は、室内8においてスピーカ装置18から受音点9に向かって直接伝播する直接音と、スピーカ装置18から室内8に拡散する拡散音との双方に基づいて受音点9での音圧レベルを算出するように構成されており、従来の直接音のみによるシミュレーションよりも高精度に受音点9の音圧レベルを算出することができる。また本実施形態の音響シミュレーション装置2は、拡散音を算出するために、室内8の総表面積と室内8における平均吸音率との2つの情報をユーザによる入力情報としている。そのため、この音響シミュレーション装置2を利用して音圧レベルのシミュレーションを行う場合、ユーザは室内8の複雑な三次元形状データを入力する必要がなく、また壁面ごとに吸音率を入力する必要もない。それ故、シミュレーションを行うための条件の入力作業を速やかに完了することができ、ユーザにかかる負担を低減することができる。
As described above, the
また本実施形態では、代表的な複数種類の平均吸音率データ26を予め記憶しているので、室内8の平均吸音率を入力する際には、その複数種類の平均吸音率データ26の中から室内8の音響空間に最も適したデータを選択することで平均吸音率を設定することが可能である。そのため、スピーカ装置18が設置された室内8における周波数ごとの平均吸音率をユーザが手計算などで算出する必要はなく、この点において入力作業を効率化できる構成となっている。ただし、必ずしも代表的な複数種類の平均吸音率データ26の中から選択しなければならないものではなく、シミュレーションの精度をより一層向上させるためには、スピーカ装置18が設置された室内8の壁面の材質などに基づいてその室内8での周波数ごとの平均吸音率を算出し、それを入力するようにしても良い。
In the present embodiment, representative plural types of average sound
また本実施形態の音響シミュレーション装置2は、従来の幾何音響シミュレーションと比較すると、シミュレーションに要する演算時間が短くなるので、条件を変更して繰り返しシミュレーションを行う場合であっても効率的にその結果を参照することができるようになる。
In addition, the
以上、本発明の一実施形態について説明したが、本発明は上述した内容に限定されるものではない。 As mentioned above, although one Embodiment of this invention was described, this invention is not limited to the content mentioned above.
例えば、上述した実施形態では、スピーカ装置18がアレイスピーカである場合について説明したが、本発明は必ずしもアレイスピーカであることに限定されるものではない。例えばスピーカユニットが1個の場合には、シミュレーションを行う際、上記各数式においてN=1として演算を行うことにより、受音点9での音圧レベルLを算出することができる。またスピーカユニットが複数ある場合であっても、指向性を制御可能なアレイスピーカではなく、2WAYスピーカや3WAYスピーカなどのように帯域が分かれているスピーカ装置の場合には、上記各数式においてN=1とした演算式でシミュレーションを行うことができる。
For example, in the above-described embodiment, the case where the
1 音響システム
2 音響シミュレーション装置
8 室内
9 受音点
10 信号処理部
18 スピーカ(音源)
20 情報処理装置
24 ハードディスク装置(記憶手段)
25 音響シミュレーションプログラム
26 平均吸音率データ
31 入力受付部(入力受付手段)
33 第1の音圧レベル算出部(第1の音圧レベル算出手段)
34 第2の音圧レベル算出部(第2の音圧レベル算出手段)
35 第3の音圧レベル算出部(第3の音圧レベル算出手段)
37 制御部(制御手段)
DESCRIPTION OF
20
25
33 1st sound pressure level calculation part (1st sound pressure level calculation means)
34 2nd sound pressure level calculation part (2nd sound pressure level calculation means)
35 Third sound pressure level calculation unit (third sound pressure level calculation means)
37 Control part (control means)
Claims (5)
前記室内の総表面積と、前記室内における平均吸音率との入力を受け付ける入力受付手段と、
前記音源から前記受音点に向かって直接伝播する直接音の音圧レベルを算出する第1の音圧レベル算出手段と、
前記入力受付手段によって受け付けられた前記総表面積と前記平均吸音率とに基づいて、前記音源から前記室内に拡散する拡散音の音圧レベルを算出する第2の音圧レベル算出手段と、
前記第1の音圧レベル算出手段で算出される直接音の音圧レベルと、前記第2の音圧レベル算出手段で算出される拡散音の音圧レベルとに基づいて前記受音点での音圧レベルを算出する第3の音圧レベル算出手段と、
を備えることを特徴とする音響シミュレーション装置。 An acoustic simulation device for calculating a sound pressure level at a predetermined sound receiving point in the room when sound is emitted from a sound source provided in the room by calculation,
Input receiving means for receiving input of the total surface area in the room and the average sound absorption coefficient in the room;
First sound pressure level calculating means for calculating a sound pressure level of a direct sound that propagates directly from the sound source toward the sound receiving point;
Second sound pressure level calculating means for calculating a sound pressure level of diffused sound diffused from the sound source into the room based on the total surface area and the average sound absorption coefficient received by the input receiving means;
Based on the sound pressure level of the direct sound calculated by the first sound pressure level calculating means and the sound pressure level of the diffuse sound calculated by the second sound pressure level calculating means, Third sound pressure level calculating means for calculating the sound pressure level;
An acoustic simulation apparatus comprising:
前記入力受付手段は、前記記憶手段に記憶された前記複数種類のデータの中から一のデータの選択を受け付けることを特徴とする請求項1記載の音響シミュレーション装置。 A storage means for storing a plurality of types of data relating to the average sound absorption rate in advance;
The acoustic simulation apparatus according to claim 1, wherein the input receiving unit receives selection of one data from the plurality of types of data stored in the storage unit.
前記第1および第2の音圧レベル算出手段は、前記アレイスピーカの指向性を制御するための複数の制御パターンに基づいて音圧レベルを算出し、
前記第3の音圧レベル算出手段は、前記複数の制御パターンのそれぞれについて前記受音点での音圧レベルを算出し、
前記制御手段は、前記複数の制御パターンのそれぞれを適用して算出された前記受音点での音圧レベルが略同一レベルとなるように制御することを特徴とする請求項3記載の音響シミュレーション装置。 The sound source is an array speaker capable of controlling sound directivity,
The first and second sound pressure level calculation means calculate a sound pressure level based on a plurality of control patterns for controlling the directivity of the array speaker,
The third sound pressure level calculating means calculates a sound pressure level at the sound receiving point for each of the plurality of control patterns;
4. The acoustic simulation according to claim 3, wherein the control means controls the sound pressure level at the sound receiving point calculated by applying each of the plurality of control patterns to be substantially the same level. apparatus.
前記情報処理装置に、
前記室内の総表面積と、前記室内における平均吸音率との入力を受け付ける入力受付ステップと、
前記音源から前記受音点に向かって直接伝播する直接音の音圧レベルを算出する第1の音圧レベル算出ステップと、
前記入力受付ステップにおいて受け付けられた前記総表面積と前記平均吸音率とに基づいて、前記音源から前記室内に拡散する拡散音の音圧レベルを算出する第2の音圧レベル算出ステップと、
前記第1の音圧レベル算出ステップで算出される直接音の音圧レベルと、前記第2の音圧レベル算出ステップで算出される拡散音の音圧レベルとを加算することにより前記受音点での音圧レベルを算出する第3の音圧レベル算出ステップと、
を実行させることを特徴とする音響シミュレーションプログラム。 An acoustic simulation program that is executed by an information processing apparatus and causes the information processing apparatus to perform a process of calculating a sound pressure level at a predetermined sound receiving point in the room when sound is emitted from a sound source provided in the room. And
In the information processing apparatus,
An input receiving step for receiving input of the total surface area in the room and the average sound absorption coefficient in the room;
A first sound pressure level calculating step of calculating a sound pressure level of a direct sound that propagates directly from the sound source toward the sound receiving point;
A second sound pressure level calculating step for calculating a sound pressure level of diffused sound diffused from the sound source into the room based on the total surface area and the average sound absorption coefficient received in the input receiving step;
The sound receiving point is obtained by adding the sound pressure level of the direct sound calculated in the first sound pressure level calculating step and the sound pressure level of the diffused sound calculated in the second sound pressure level calculating step. A third sound pressure level calculating step for calculating the sound pressure level at
An acoustic simulation program characterized in that
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009027622A JP2010181824A (en) | 2009-02-09 | 2009-02-09 | Sound simulation apparatus and sound simulation program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009027622A JP2010181824A (en) | 2009-02-09 | 2009-02-09 | Sound simulation apparatus and sound simulation program |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010181824A true JP2010181824A (en) | 2010-08-19 |
Family
ID=42763419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009027622A Pending JP2010181824A (en) | 2009-02-09 | 2009-02-09 | Sound simulation apparatus and sound simulation program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2010181824A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012242597A (en) * | 2011-05-19 | 2012-12-10 | Hitachi Ltd | Acoustic simulator, acoustic consulting device, and processing method for them |
US20160094925A1 (en) * | 2014-09-29 | 2016-03-31 | Bose Corporation | Systems and methods for determining metric for sound system evaluation |
WO2022034700A1 (en) * | 2020-08-11 | 2022-02-17 | ツナガルデザイン株式会社 | Sound-leakage-prevention application program and sound-leakage-prevention device |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001025084A (en) * | 1999-07-07 | 2001-01-26 | Matsushita Electric Ind Co Ltd | Speaker system |
JP2001275198A (en) * | 2000-03-28 | 2001-10-05 | Matsushita Electric Ind Co Ltd | Sound volume controller |
JP2001346299A (en) * | 2000-05-31 | 2001-12-14 | Sony Corp | Sound field correction method and audio unit |
JP2002123262A (en) * | 2000-10-18 | 2002-04-26 | Matsushita Electric Ind Co Ltd | Device and method for simulating interactive sound field, and recording medium with recorded program thereof |
JP2003157090A (en) * | 2001-11-20 | 2003-05-30 | Matsushita Electric Ind Co Ltd | Reverberation sound generating method and reverberation sound simulator |
JP2007026221A (en) * | 2005-07-19 | 2007-02-01 | Yamaha Corp | Acoustic design support device and acoustic design support program |
-
2009
- 2009-02-09 JP JP2009027622A patent/JP2010181824A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001025084A (en) * | 1999-07-07 | 2001-01-26 | Matsushita Electric Ind Co Ltd | Speaker system |
JP2001275198A (en) * | 2000-03-28 | 2001-10-05 | Matsushita Electric Ind Co Ltd | Sound volume controller |
JP2001346299A (en) * | 2000-05-31 | 2001-12-14 | Sony Corp | Sound field correction method and audio unit |
JP2002123262A (en) * | 2000-10-18 | 2002-04-26 | Matsushita Electric Ind Co Ltd | Device and method for simulating interactive sound field, and recording medium with recorded program thereof |
JP2003157090A (en) * | 2001-11-20 | 2003-05-30 | Matsushita Electric Ind Co Ltd | Reverberation sound generating method and reverberation sound simulator |
JP2007026221A (en) * | 2005-07-19 | 2007-02-01 | Yamaha Corp | Acoustic design support device and acoustic design support program |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012242597A (en) * | 2011-05-19 | 2012-12-10 | Hitachi Ltd | Acoustic simulator, acoustic consulting device, and processing method for them |
US20160094925A1 (en) * | 2014-09-29 | 2016-03-31 | Bose Corporation | Systems and methods for determining metric for sound system evaluation |
US9723419B2 (en) * | 2014-09-29 | 2017-08-01 | Bose Corporation | Systems and methods for determining metric for sound system evaluation |
WO2022034700A1 (en) * | 2020-08-11 | 2022-02-17 | ツナガルデザイン株式会社 | Sound-leakage-prevention application program and sound-leakage-prevention device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8332060B2 (en) | Acoustic design support apparatus | |
US8693705B2 (en) | Response waveform synthesis method and apparatus | |
CN111535470B (en) | Method and system for improving indoor language definition | |
WO2008047089A2 (en) | Speaker configuration | |
CN110809215B (en) | Loudspeaker signal feeding method for conference system | |
Jeon et al. | Objective and subjective assessment of sound diffuseness in musical venues via computer simulations and a scale model | |
JP2010181824A (en) | Sound simulation apparatus and sound simulation program | |
US20230306953A1 (en) | Method for generating a reverberation audio signal | |
Aspöck | Validation of room acoustic simulation models | |
TWI640983B (en) | Method for simulating room acoustics effect | |
Jo et al. | Optimizing sound diffusion in a concert hall using scale-model measurements and simulations | |
JP5286739B2 (en) | Sound image localization parameter calculation device, sound image localization control device, sound image localization device, and program | |
US9723419B2 (en) | Systems and methods for determining metric for sound system evaluation | |
Pignatelli et al. | A wooden acoustic shell for open-air chamber music concert | |
JP3482055B2 (en) | High-accuracy sound ray tracking device and high-accuracy sound ray tracking method | |
JP4726174B2 (en) | Indoor acoustic design method | |
Kostek et al. | Modeling and designing acoustical conditions of the interior–case study | |
JP2003042836A (en) | Method for estimating vibration force in construction machine, method and apparatus for estimating vibration and noise, record medium, program, and apparatus for feeling noise in building construction | |
Jo et al. | Objective assessment of a scattered sound field in simulated concert halls and scale model | |
Agus et al. | Energy-Based Binaural Acoustic Modeling | |
Stumpf González et al. | Análisis acústico del Auditório Padre Werner em Unisinos | |
Feistel et al. | Acoustic Modelling–Basics | |
Seeralaselvan | Acoustics Simulation for Stadium Design Using EASE: Analyzing Acoustics and Providing Retrofit Options for the Los Angeles Memorial Coliseum | |
JP2008051935A (en) | Sound insulating structure selecting apparatus | |
Girón et al. | Objective Measures of Spatial Effects in Spanish Concert Halls |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20111216 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130311 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130319 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130515 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20130813 |