(57)【要約】【課題】  基準周波数の数を増やすことなく、より自然な聴感をもたらすインパルス応答データを得る手段を提供する。【解決手段】  補間反射吸音率算定部２０４は、基準周波数ポイントに関する反射吸音率データＤ９０４を用いて、補間周波数ポイントに関する補間反射吸音率データＤ２０８を作成する。ポイントエコーダイアグラム算定部９０２は、反射吸音率データＤ９０４および補間反射吸音率データＤ２０８を用いて、基準周波数ポイントおよび補間周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７を作成する。帯域エコーダイアグラム算定部９０３は、周波数軸上でクロスフェードする時間波形を用いて、ポイントエコーダイアグラムデータＤ９０７から帯域エコーダイアグラムデータＤ９０８を作成する。加算部９０４は、帯域エコーダイアグラムデータＤ９０８を加算してインパルス応答データＤ９０９を作成する。【選択図】    図２
2006030689.tif 000002
【特許請求の範囲】
【請求項１】
音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定する過程と、
前記複数の基準周波数の各々に関するエコーダイアグラムを補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定する過程と、
前記複数の基準周波数および前記１以上の補間周波数の各々に関し、当該基準周波数もしくは当該補間周波数を中心周波数とし、当該中心周波数を含む通過帯域においてのみ正値の周波数成分を有する時間波形を表す係数列を、当該基準周波数もしくは当該補間周波数に対応したエコーダイアグラムに畳み込む過程と、
前記複数の基準周波数および前記１以上の補間周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する過程と
を有するインパルス応答データ生成方法。
【請求項２】
前記特徴データは、前記複数の基準周波数の各々における前記音場の境界面の吸音率を示す情報を含み、
前記インパルス応答データ生成方法は、
前記複数の基準周波数の各々における吸音率を補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々における吸音率を算定する過程を備え、
前記特徴データ、前記音源位置データ、前記受音位置データおよび前記１以上の補間周波数の各々における吸音率を用いて、前記１以上の補間周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記１以上の補間周波数の各々について算定することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定する
ことを特徴とする請求項１に記載のインパルス応答データ生成方法。
【請求項３】
音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定する過程と、
前記複数の基準周波数の各々に関し、当該基準周波数を中心周波数とし、各中心周波数の低域側および高域側に向けてスロープを描いて減衰する周波数成分を有する時間波形を表す係数列を、当該基準周波数に対応したエコーダイアグラムに畳み込む過程と、
前記複数の基準周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する過程と
を有するインパルス応答データ生成方法。
【請求項４】
音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定する過程と、
前記複数の基準周波数の各々に関するエコーダイアグラムを補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定する過程と、
前記複数の基準周波数および前記１以上の補間周波数の各々に関し、当該基準周波数もしくは当該補間周波数を中心周波数とし、各中心周波数の低域側および高域側に向けてスロープを描いて減衰する周波数成分を有する時間波形を表す係数列を、当該基準周波数に対応したエコーダイアグラムに畳み込む過程と、
前記複数の基準周波数および前記１以上の補間周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する過程と
を有するインパルス応答データ生成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音場シミュレーションのためのデータ生成技術に関する。
【背景技術】
【０００２】
音データに対し、シミュレーションにより、特定の音場において音が発せられた場合の残響音等の音場情報を付加する技術がある。そのようなシミュレーションの方法の１つとして、音線法と呼ばれるものがある。
【０００３】
図１２は、この種の音線法により音データに音場情報を付加する装置（以下、「音場シミュレータ９０」と呼ぶ）の機能構成を示した図である。以下、音場シミュレータ９０の動作を説明する。まず、音場シミュレータ９０は、予め以下のデータを記憶部９００に記憶している。
（ａ）形状データＤ９０１：音場の形状を示すデータであり、音場と外界との各境界面の頂点群の座標を示す座標データ列の集まりである。
（ｂ）音源位置データＤ９０２：音場内の音源位置の座標を示すデータである。
（ｃ）受音位置データＤ９０３：音場内の受音位置の座標を示すデータである。
（ｄ）反射吸音率データＤ９０４：基準となる周波数ポイント（以下、「基準周波数ポイント」と呼ぶ）の音が音場と外界との境界面において反射する場合に、反射直前の音の強さと反射直後の音の強さの差を、反射直前の音の強さで除した値を示すデータである。反射吸音率データＤ９０４は、例えば、３１．３Ｈｚ、６２．５Ｈｚ、１２５Ｈｚ、・・・、１６ｋＨｚのようにオクターブ間隔で並んだ複数の基準周波数ポイントの各々に関し準備されている。
（ｅ）音データＤ９０５：音場情報が付加される前の音を示すデータである。
【０００４】
経路算定部９０１は、形状データＤ９０１、音源位置データＤ９０２および受音位置データＤ９０３を用いて、音源位置データＤ９０２により示される音源位置から音が発せされた場合に、その音が受音位置データＤ９０３により示される受音位置に到達するまでの経路を算定し、その経路の始点である音源位置、終点である受音位置および経由点である反射位置の各々を示す座標列からなる経路データを生成する。
【０００５】
経路算定部９０１が算定する音の経路は、音源位置と受音位置を直線で結んで得られる経路に加え、音源位置から発せられた音が、形状データＤ９０１により示される音場と外界との境界面において反射した後、受音位置に到達する経路も含んでいる。そのような経路は無限数あるが、経路算定部９０１は、例えば反射回数が一定数以下、経路の総延長が一定距離以下、といった所定の条件を満たす経路についてのみ、経路データを生成する。経路算定部９０１により生成された複数の経路データは、経路データ群Ｄ９０６として一時記憶部９０６に記憶される。
【０００６】
ポイントエコーダイアグラム算定部９０２は、経路データ群Ｄ９０６および反射吸音率データＤ９０４を用いて、基準周波数ポイントの各々に関し、基準周波数ポイントの音が音源位置において発せられた場合に受音位置において受音される音の時間と強さの関係を示すエコーダイアグラム（以下、周波数ポイントの音に関するエコーダイアグラムを「ポイントエコーダイアグラム」と呼ぶ）を算定し、そのエコーダイアグラムを示すポイントエコーダイアグラムデータＤ９０７を生成する。図１３は、ポイントエコーダイアグラムデータＤ９０７の内容を例示した図である。ポイントエコーダイアグラム算定部９０２により生成されたポイントエコーダイアグラムデータＤ９０７は、一時記憶部９０６に記憶される。以下、例として、ポイントエコーダイアグラム算定部９０２が１２５Ｈｚの基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７を生成する方法を説明する。
【０００７】
ポイントエコーダイアグラム算定部９０２は、経路データ群Ｄ９０６に含まれる経路データの各々に関し、以下の処理を行う。まず、ポイントエコーダイアグラム算定部９０２は、経路データに含まれる隣り合う反射位置の座標間の距離を全て加算し、その経路データにより示される経路の総延長を算出する。また、ポイントエコーダイアグラム算定部９０２は、経路データに含まれる反射位置の数をカウントする。続いて、ポイントエコーダイアグラム算定部９０２は、先に算出した経路の総延長を音の空気中の伝播速度で除することにより、音源位置において発せられた音が経路データにより示される経路をたどって受音位置へ到達するまでの時間（以下、「到達時間」と呼ぶ）を算出する。また、ポイントエコーダイアグラム算定部９０２は、経路の総延長の自乗に４πを乗じた値の逆数を算出する。そのようにして得られる数値は、音源位置において強さ「１．０」で発せられた音が発散により減衰した後、受音位置において受音される際の強さを示す数値（以下、「発散後音強値」と呼ぶ）である。ただし、この発散後音強値には、音の反射に伴う減衰が考慮されていない。以上のように算出される経路の総延長、反射位置の数、到達時間および発散後音強値を示すデータは、周波数により左右されないデータであるため、１２５Ｈｚ以外の基準周波数ポイントに関しても共通して利用される。
【０００８】
続いて、ポイントエコーダイアグラム算定部９０２は、上記の発散後音強値に対し、１２５Ｈｚに対応した反射吸音率データＤ９０４により示される比率を１から減じた数値、すなわち、減衰率を、先にカウントした反射位置の数だけ繰り返し乗じる。その結果得られる数値は、音源位置において「１．０」のレベルで発せられた１２５Ｈｚの音が、経路データにより示される反射点を経由して受音位置に到達する場合の受音位置における音の強さを示す数値（以下、「減衰後音強値」と呼ぶ）である。
【０００９】
ポイントエコーダイアグラム算定部９０２は、上記のようにして１つの経路データにより示される経路に関し、１２５Ｈｚの音の到達時間および減衰後音強値を算出することができる。ポイントエコーダイアグラム算定部９０２は、同様の算出処理を経路データ群Ｄ９０６に含まれる他の全ての経路データに関して繰り返す。そのようにして得られる到達時間と減衰後音強値の組み合わせを示すデータ群が、１２５Ｈｚの基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７である（図１３参照）。ただし、同じ到達時間に対応する複数の減衰後音強値が存在する場合、ポイントエコーダイアグラム算定部９０２はそれらの減衰後音強値の合算値を、その到達時間に対応する減衰後音強値とする。以上が、ポイントエコーダイアグラム算定部９０２がポイントエコーダイアグラムデータＤ９０７を生成する際の処理である。
【００１０】
帯域エコーダイアグラム算定部９０３は、３１．３Ｈｚの基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７により示されるエコーダイアグラムに対し、図１４（１）に示す周波数成分を有する時間波形の係数列を畳み込むことにより、４４．２Ｈｚ以下の周波数帯域に関するエコーダイアグラムを示す帯域エコーダイアグラムデータＤ９０８を生成する。同様に、帯域エコーダイアグラム算定部９０３は、６２．５Ｈｚ、１２５Ｈｚ、・・・の基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７により示されるエコーダイアグラムの各々に対し、図１４（２）、（３）、・・・に示す周波数成分を有する時間波形の係数列を各々畳み込むことにより、それらの基準周波数ポイントを中心周波数ポイントとする周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を生成する。帯域エコーダイアグラム算定部９０３により生成された帯域エコーダイアグラムデータＤ９０８は、一時記憶部９０６に記憶される。
【００１１】
加算部９０４は、帯域エコーダイアグラム算定部９０３により生成された帯域エコーダイアグラムデータＤ９０８の全てを加算する。そのようにして得られるデータは、あらゆる周波数ポイントの音を音源位置において発音した場合に受音位置において受音される音の強さの時間的変化を擬似的に示すデータである。すなわち、加算部９０４により生成されるデータは音源位置においてインパルスを発した場合の受音位置におけるインパルス応答を擬似的に示すデータであるため、以下、これをインパルス応答データＤ９０９と呼ぶ。加算部９０４により生成されたインパルス応答データＤ９０９は、一時記憶部９０６に記憶される。
【００１２】
畳み込み演算部９０５は、音データＤ９０５に対し、インパルス応答データＤ９０９により示される係数列を畳み込むことにより、インパルス応答データＤ９０９により示される音場情報の付加された音を示す音場情報付音データＤ９１０を生成する。すなわち、音場情報付音データＤ９１０を再生した場合に得られる音は、実際に形状データＤ９０１により示される音場内の音源位置データＤ９０２により示される音源位置において、音データＤ９０５に従った音を再生した場合に、受音位置データＤ９０３により示される受音位置において受音される音と類似した音となる。
【００１３】
以上が、音線法により音データに音場情報を付加する手順の説明である。音線法によれば、実際にコンサートホール等において発音および録音を行うことなく、それと類似の結果を演算処理によって得ることができる。その際、音源位置や受音位置を様々に変更した場合の結果を得るのも容易である。さらに、音線法によれば、実在しない音場における音響効果を確認することもできる。
【００１４】
例えば、音線法に関する技術を開示したものとして、特許文献１および２がある。
【特許文献１】特開昭６０−０４６４３０号公報
【特許文献２】特開２０００−２８４７８８号公報
【発明の開示】
【発明が解決しようとする課題】
【００１５】
ところで、以上説明した従来技術により得られる音は、実際の音場において聴感される音と比較して、発音後に一定時間が経過したタイミングで、特定の周波数帯域に含まれる音の成分がカットされてゆくような聴感上の不自然さを伴う。この不自然さが生じる理由は、隣り合う周波数帯域に関するエコーダイアグラムが、周波数軸上で不連続に接合されているためである。例えば、４４．２Ｈｚ以下の周波数帯域に関する帯域エコーダイアグラムデータＤ９０８により示されるエコーダイアグラムの４４．２Ｈｚの周波数成分のレベルと、４４．２Ｈｚ〜８８．４Ｈｚの周波数帯域に関する帯域エコーダイアグラムデータＤ９０８により示されるエコーダイアグラムの４４．２Ｈｚの周波数成分のレベルとの間に、実際には存在しないはずのギャップが生じる。このギャップにより、上述した聴感上の不自然さが生じるのである。
【００１６】
上記のギャップを小さくするために、インパルス応答を得るために用いるエコーダイアグラムの周波数軸上の解像度を高めることが考えられる。すなわち、基準周波数ポイントの数を増やすことにより、各々の帯域エコーダイアグラムデータＤ９０８によりカバーされる周波数帯域の幅を狭くするのである。しかしながら、基準周波数ポイントの数を増やすと、ポイントエコーダイアグラムデータＤ９０７を生成するために必要な反射吸音率データＤ９０４等のデータ量が増大し、それらのデータを準備する手間やコストが増大する。また、反射吸音率データＤ９０４等のデータに関しては、既にオクターブ間隔の所定の基準周波数ポイントに関し、多数のライブラリデータが存在しているにもかかわらず、基準周波数ポイントの数を増やすと、それらのライブラリデータを利用できない、という不都合がある。
【００１７】
上記の状況に鑑み、本発明は、基準周波数ポイントの数を増やすことなく、より自然な聴感をもたらすインパルス応答データを得る手段を提供することを目的とする。
【課題を解決するための手段】
【００１８】
以下、本発明により、インパルス応答を得るために用いるエコーダイアグラムの周波数軸上の解像度を高める方法を説明する。まず、インパルス波形Ｉｍ（ｔ）のフーリエ変換をａ（ｊω）とすると、インパルス応答波形Ｉｍ（ｔ）は次の（式１）のように表すことができる。
【数１】
2006030689.tif 000003
【００１９】
音場における音源位置から受音位置までの周波数伝達特性をＦ（ｊω）とすると、音源位置においてインパルス波形Ｉｍ（ｔ）を発生したときに受音位置において観測されるインパルス応答波形Ｒｅｓ（ｔ）は、次の（式２）のように表すことができる。
【数２】
2006030689.tif 000004
【００２０】
ここで、実際の音場において得られるインパルス応答波形をフーリエ変換すれば、Ｆ（ｊω）を求めることができる。ただし、実際の音場においてインパルス応答波形を実測するのはコストや労力を要し、また存在しない音場においては実測が不可能という問題もある。そこで、コンピュータを用いたシミュレーションにより、インパルス応答波形を近似することが必要となる。そのために、まず、インパルス波形Ｉｍ（ｔ）を以下の（式３）のように近似する。ここで、下記の（式３）においてｎを∞に、Δωを０に近づけた極限が上記（式２）である。
【数３】
2006030689.tif 000005
【００２１】
インパルス波形Ｉｍ（ｔ）を構成する１つのスペクトルａ（ｊｋΔω）ｅjkΔωtに対応した応答波形Ｒｅｓk（ｔ）は次の（式４）のようになる。
【数４】
2006030689.tif 000006
【００２２】
インパルス応答波形Ｒｅｓ（ｔ）は、インパルス波形Ｉｍ（ｔ）を構成する全てのスペクトルａ（ｊｋΔω）ｅjkΔωtに対応した応答波形Ｒｅｓk（ｔ）の総和と考えることができるので、次の（式５）のように表すことができる。
【数５】
2006030689.tif 000007
【００２３】
上記（式３）および（式４）において、Ｆ（ｊｋΔω）は、音源位置から受音位置までの音の伝播経路を各種想定し、各伝播経路における遅延時間と減衰率を算出することにより、例えば以下のように求めることができる。
【数６】
2006030689.tif 000008
【００２４】
従って、音源位置から受音位置までの全伝播経路に対応した伝達関数の総和、すなわち、音源位置から受音位置までの伝達関数は、次の（式６）のように表すことができる。上述したエコーダイアグラムは、この式（６）に相当するものである。
【数７】
2006030689.tif 000009
・・・（式６）
【００２５】
ここで、より自然な音場情報を再現するためには、インパルス波形Ｉｍ（ｔ）の周波数分解能を高くする、すなわち、ａ（ｊｋΔω）ｅjkΔωtの個数ｎを大きくする必要がある。一方、インパルス波形Ｉｍ（ｔ）の周波数分解能に応じて、Ｆ（ｊω）の周波数分解能を高くし、個数ｎの角周波数ωについて上記（式６）を求めることは、減衰率σおよび遅延時間τを算出するためのパラメータを角周波数ωの各々に関し準備する必要があるため容易ではなく、また膨大な演算量を必要とする。従って、インパルス波形Ｉｍ（ｔ）の周波数分解能に対し、Ｆ（ｊω）の周波数分解能を低くする手段が求められる。
【００２６】
従来技術においては、基準となる角周波数ωに関するＦ（ｊｋΔω）をそのまま、近隣の角周波数ωに関するａ（ｊｋΔω）ｅjkΔωtに対応する伝達関数として用いることにより、インパルス波形Ｉｍ（ｔ）の周波数分解能に対し、Ｆ（ｊω）の周波数分解能を低くしていた。例えば、以下の通りである。ただし、「→」は右が左に対応付けられる伝達関数であることを示す。
【数８】
2006030689.tif 000010
この結果、従来技術は、上述した不都合を生んだ。
【００２７】
これに対し、本発明においては、各々の角周波数ωに関するａ（ｊｋΔω）ｅjkΔωtの伝達関数を、その低域側および高域側にそれぞれ位置する基準となる角周波数ωに関するＦ（ｊｋΔω）を補間することにより求めることとした。例えば、以下の通りである。
【数９】
2006030689.tif 000011
【００２８】
より具体的には、本発明では、音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定し、前記複数の基準周波数の各々に関するエコーダイアグラムを補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定し、前記複数の基準周波数および前記１以上の補間周波数の各々に関し、当該基準周波数もしくは当該補間周波数を中心周波数とし、当該中心周波数を含む通過帯域においてのみ正値の周波数成分を有する時間波形を表す係数列を、当該基準周波数もしくは当該補間周波数に対応したエコーダイアグラムに畳み込み、前記複数の基準周波数および前記１以上の補間周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する。
【００２９】
また、好ましい態様においては、１の中心周波数を含む１の通過帯域の上限と、前記１の中心周波数より高域側に位置し前記１の中心周波数と隣り合う他の中心周波数を含む他の通過帯域の下限周波数が一致する。
【００３０】
上記の補間を行う場合の係数、すなわちｂおよびｃの値を決定する好ましい態様として、基準となる角周波数ωに関する減衰率σおよび遅延時間τを算出するためのパラメータを補間し、補間したパラメータにより補間した角周波数ωに関する減衰率σおよび遅延時間τを算出することにより、上記のｂおよびｃを決定する態様が考えられる。
【００３１】
より具体的には、前記複数の基準周波数の各々における前記音場の境界面の吸音率を示す情報を含む前記特徴データを用いて、前記インパルス応答データ生成方法は、前記複数の基準周波数の各々における吸音率を補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々における吸音率を算定し、前記特徴データ、前記音源位置データ、前記受音位置データおよび前記１以上の補間周波数の各々における吸音率を用いて、前記１以上の補間周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記１以上の補間周波数の各々について算定することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定する。
【００３２】
ところで、隣り合う周波数帯域に関するエコーダイアグラムが、周波数軸上で不連続に接合される問題を解決するための他の手段として、周波数軸上でクロスフェードする周波数成分を有する時間波形の係数列を、基準となる角周波数ωに関する伝達関数を時間軸上に展開したエコーダイアグラムに畳み込むことが考えられる。
【００３３】
より具体的には、音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定し、前記複数の基準周波数の各々に関し、当該基準周波数を中心周波数とし、各中心周波数の低域側および高域側に向けてスロープを描いて減衰する周波数成分を有する時間波形を表す係数列を、当該基準周波数に対応したエコーダイアグラムに畳み込み、前記複数の基準周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する。
【００３４】
また、好ましい態様においては、前記畳み込みにおいて用いられる時間波形の周波数成分が正値をとる帯域は、隣り合う帯域間で重複部分を有する。さらに好ましい態様においては、前記重複部分における周波数成分の和は１である。
【００３５】
さらに、上述した補間周波数に関するエコーダイアグラムの算定と、周波数軸上でクロスフェードする周波数成分を有する時間波形の係数列の畳み込みを組み合わせることにより、隣り合う周波数帯域に関するエコーダイアグラムを周波数軸上でより滑らかに接合することが可能となる。
【００３６】
より具体的には、音場の特徴を示す特徴データ、前記音場内の音源位置を示す音源位置データおよび前記音場内の受音位置を示す受音位置データを用いて、複数の基準周波数の各々の音を前記音源位置において発生させた場合に前記受音位置に到達する当該音のレベルの時間的変化を表すエコーダイアグラムを前記複数の基準周波数の各々について算定し、  前記複数の基準周波数の各々に関するエコーダイアグラムを補間することにより、前記複数の基準周波数の間の１以上の補間周波数の各々におけるエコーダイアグラムを算定し、  前記複数の基準周波数および前記１以上の補間周波数の各々に関し、当該基準周波数もしくは当該補間周波数を中心周波数とし、各中心周波数の低域側および高域側に向けてスロープを描いて減衰する周波数成分を有する時間波形を表す係数列を、当該基準周波数に対応したエコーダイアグラムに畳み込み、前記複数の基準周波数および前記１以上の補間周波数の各々に関する前記畳み込みを施されたエコーダイアグラムを加算することにより前記音場における音の伝達特性を表すインパルス応答波形データを生成する。
【００３７】
かかる構成によるインパルス応答データ生成方法によれば、基準周波数の解像度を上げることなく、周波数軸上の不連続性が低減されたインパルス応答データが生成される。
【００３８】
また、本発明は、上記のインパルス応答データの生成方法をデータ処理により実現する装置を提供する。さらに、本発明は、上記のインパルス応答データの生成方法の各々を、コンピュータに所定のデータ処理を行わせることにより実現するための手順を示すプログラムを提供する。
【発明の効果】
【００３９】
本発明にかかる方法によれば、基準周波数ポイントの数を増やすことなく、より自然な聴感をもたらすインパルス応答データを得ることができる。
【発明を実施するための最良の形態】
【００４０】
［１．第１実施形態］
本発明の第１実施形態にかかるデータ処理装置（以下、「音場シミュレータ１０」と呼ぶ）は、音データに対し音場情報を付加するためのインパルス応答データを生成し、生成したインパルス応答データを音データに畳み込む処理を行うことにより、音場情報の付加された音を示す音データを生成する装置である。
【００４１】
［１．１．音場シミュレータの構成］
図１は、音場シミュレータ１０のハードウェア構成を示す図である。音場シミュレータ１０は、一般的なコンピュータと同様に、各種データ処理を行うとともに他の構成部の制御を行うＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１０１、ＣＰＵ１０１に基本的処理を行わせるＢＩＯＳ（Ｂａｓｉｃ  Ｉｎｐｕｔ  Ｏｕｔｐｕｔ  Ｓｙｓｔｅｍ）等を記憶するＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）１０２、ＣＰＵ１０１のデータ処理において一時的にデータを記憶するＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）１０３、音場シミュレータ１０の機能を実現するためのアプリケーションプログラム（以下、「音場シミュレーションプログラム」と呼ぶ）や各種データを記憶するＨＤ（Ｈａｒｄ  Ｄｉｓｋ）１０４、ＣＰＵ１０１の指示に従いユーザに対し音データの選択等を促す画面等を表示するディスプレイ１０５、ユーザの操作に応じた信号を送信するキーボード等の操作部１０６、Ｄ／Ａ（Ｄｉｇｉｔａｌ  ｔｏ  Ａｎａｌｏｇ）コンバータ、アンプおよびスピーカを備え音データに基づき発音を行うサウンドシステム１０７および外部機器との間でデータの送受信を行う入出力Ｉ／Ｆ（Ｉｎｔｅｒｆａｃｅ）１０８を有している。
【００４２】
ＣＰＵ１０１は、音場シミュレーションプログラムを実行することにより、音場シミュレータ１０の各種機能を実現する。図２は、音場シミュレータ１０の機能構成を示したブロック図である。ただし、音場シミュレータ１０がユーザとの間で情報をやりとりするための表示部や操作部、外部機器とのデータ通信を行う通信部等は、一般的な情報処理機器におけるものと同様であるため、図２においては省略されている。また、音場シミュレータ１０の各構成部は、記憶部９００および一時記憶部９０６に含まれる様々なデータを利用するが、その全てを示すと煩雑となるため、図２においては主要なデータの利用関係のみが矢印で示されている。音場シミュレータ１０の機能は、図１２に示した音場シミュレータ９０の機能と共通部分を有している。従って、図２においては、図１２と同様もしくは類似の機能構成部に対しては同じ符号が付されている。
【００４３】
音場シミュレータ１０の記憶部９００は、音場シミュレータ９０の記憶部９００が記憶する各種データと同様もしくは類似のものとして、予め以下のデータを記憶している。ただし、以下の説明において、座標データにより示される数値の単位はメートルであるものとする。また、以下に示すデータの数値は説明のためのものであり、必ずしも実際の数値を意味しない。
（ａ）形状データＤ９０１：音場の形状および音場と外界との各境界面の材質を示すデータであり、境界面の頂点群の座標を示す座標データ列と境界面の材質を示す材質データとを対応付けたデータ（以下、「境界面データ」と呼ぶ）の集まりである。図３に形状データＤ９０１の例を示す。
（ｂ）音源位置データＤ９０２：音場内に含まれる音源位置および音源が音を発する方向を示すデータであり、音源位置の座標データと音を発する方向を示すベクトルデータの組み合わせである。例えば、音源位置データＤ９０２は「（６，５，３），（０，０．８７，０．５）」のデータ形式をとり、この場合、座標データ（６，５，３）で示される音源位置から、ベクトルデータ（０，０．８７，０．５）で示される方向に向けて音が発せされることを示している。
（ｃ）受音位置データＤ９０３：音場内に含まれる受音位置を示す座標データである。
（ｄ）反射吸音率データＤ９０４：各基準周波数ポイントの音が音場と外界との境界面において反射する場合に、反射直前の音の強さと反射直後の音の強さの差を、反射直前の音の強さで除した値を示す吸音率データを、境界面の材質を示す材質データおよび音の入射角を示す入射角データと対応付けたデータ（以下、「材質入射角別吸音率データ」と呼ぶ）の集まりである。図４に反射吸音率データＤ９０４の例を示す。
（ｅ）音データＤ９０５：音場情報が付加される前の音を示すデータである。
【００４４】
音場シミュレータ１０の記憶部９００は、さらに以下のデータを予め記憶している。
（ｆ）音源名データＤ２００：音源の名称を示すデータである。
（ｇ）温度湿度データＤ２０１：音場内の空気の温度および湿度を示すデータである。
（ｈ）補間数テンプレートデータ群Ｄ２０２：隣り合う基準周波数ポイントの間を補間する周波数ポイント（以下、「補間周波数ポイント」と呼ぶ）の数を低域側から順に示す補間数データを、楽曲等の分類名を示すジャンルデータと対応付けたデータ（以下、「補間数テンプレートデータ」と呼ぶ）の集まりである。図５に補間数テンプレートデータ群Ｄ２０２の例を示す。例えば、図５の先頭行に示される補間数テンプレートデータは、「ロック」のジャンルに属する楽曲を示す音データに対し音場情報を付加する場合に適する補間周波数ポイントの数を示しており、補間数データにより示される（１，１，・・・，３）は、順に「３１．３Ｈｚ〜６２．５Ｈｚ」、「６２．５Ｈｚ〜１２５Ｈｚ」、・・・、「８ｋＨｚ〜１６ｋＨｚ」の各周波数帯域に対応している。各補間数テンプレートデータに含まれる補間数データは、楽曲のジャンル等に応じて、それらの特徴を顕著に含む周波数帯域ほど、大きな補間数を示すように作成されている。
（ｉ）指向性データＤ２０３：基準点を中心とし単位距離（例えば１メートル）を半径とする球面上の各点の位置を示すベクトルデータに、基準点において特定の音源により基準周波数ポイントの音を基準となる方向、例えば（０，１，０）に向けて発音した場合に、球面上の位置において受音される音の強さを示す音強データを対応付けたデータ（以下、「指向性要素データ」と呼ぶ）の集まりである。指向性データＤ２０３は、全体として、各基準周波数ポイントの音に関する音源の指向性を示している。指向性データＤ２０３は、複数の音源の各々に関し準備されている。図６に指向性データＤ２０３の例を示す。図６の先頭行の指向性要素データは、ベクトルデータ（０，１，０）により示される方向、すなわち基準となる方向における音源位置から１メートルの位置で、「スピーカＳＰ１」により発音された３１．３Ｈｚの音が、「０．９７」の強さで受音されることを示している。また、図６の第２行の指向性要素データは、基準となる方向からわずかにｘ軸正方向に傾いた方向における音源位置から１メートルの位置で、「スピーカＳＰ１」により発音された３１．３Ｈｚの音が、「０．９５」の強さで受音されることを示している。音源の種類はスピーカに限らず楽器などでもよい。
（ｊ）空気吸音率データＤ２０４：各基準周波数ポイントの音が単位距離（例えば１メートル）だけ伝播する場合に、伝播前の音の強さと伝播後の音の強さの差を、伝播前の音の強さで除した値を示す吸音率データを、空気の温度を示す温度データおよび湿度を示す湿度データと対応付けたデータ（以下、「温度湿度別吸音率データ」と呼ぶ）の集まりである。図７に空気吸音率データＤ２０４の例を示す。
【００４５】
上記の（ａ）〜（ｊ）のデータのうち、（ａ）形状データＤ９０１、（ｂ）音源位置データＤ９０２、（ｃ）受音位置データＤ９０３、（ｅ）音データＤ９０５、（ｆ）音源名データＤ２００および（ｇ）温度湿度データＤ２０１は、ユーザがシミュレーションを行いたいと考える音場の環境等に応じて、自由に変更もしくは選択することができる。また、その他のデータに関しては、ユーザが自ら作成等する他、ライブラリとして既に作成されているものを利用することができる。
【００４６】
［１．２．音場シミュレータの機能および動作］
以下、音場シミュレータ１０が備える機能およびその動作を説明する。まず、経路算定部９０１は、形状データＤ９０１、音源位置データＤ９０２および受音位置データＤ９０３を用いて、経路データ群Ｄ９０６を生成する。経路算定部９０１が経路データ群Ｄ９０６を生成する方法は、音場シミュレータ９０における場合と同様である。
【００４７】
経路データ群Ｄ９０６が一時記憶部９０６に記憶されると、音場シミュレータ１０は続いてデータの補間処理を行う。そのために、音場シミュレータ１０は、ユーザに対し補間数を決定する方法に関する以下の選択肢のいずれかを選択するよう促す画面を表示する。
（ａ）固定数「３」
（ｂ）個別に指定
（ｃ）テンプレートを選択
（ｄ）基準データに基づき自動設定
（ｅ）音データに基づき自動設定
【００４８】
補間数決定部２０１は、上記の選択肢を示す画面に対するユーザの操作に応じて、基準周波数ポイント間を補間する補間周波数ポイントの数を決定する。ユーザにより、上記選択肢（ａ）が選択された場合、補間数決定部２０１は（３，３，・・・，３）という内容の補間数データＤ２０５を作成する。補間数データＤ２０５に含まれる数値は、上述した補間数テンプレートデータ群Ｄ２０２（図５参照）に含まれる補間数データと同様に、それぞれ「３１．３Ｈｚ〜６２．５Ｈｚ」、「６２．５Ｈｚ〜１２５Ｈｚ」、・・・、「８ｋＨｚ〜１６ｋＨｚ」の各周波数帯域を補間する補間周波数ポイントの数を示している。
【００４９】
ユーザにより上記選択肢（ｂ）が選択された場合、音場シミュレータ１０は、さらにユーザに対し、各周波数帯域における補間周波数ポイントの数の入力を促す画面を表示する。その画面に応じて、ユーザが数値を入力すると、補間数決定部２０１はユーザにより入力された数値の列を示すデータを補間数データＤ２０５として作成する。ここで、ユーザは全ての隣り合う基準周波数ポイント間に関し、補間周波数ポイントの数として「０」を入力することもできる。その場合、以下に説明する一連の補間処理はスキップされる。
【００５０】
ユーザにより上記選択肢（ｃ）が選択された場合、音場シミュレータ１０は、さらにユーザに対し、補間数テンプレートデータ群Ｄ２０２に含まれるジャンルデータの内容を示すリストからいずれかを選択するよう促す画面を表示する。その画面に応じて、ユーザがいずれかのジャンルを選択すると、補間数決定部２０１はユーザにより選択されたジャンルを示すジャンルデータを含む補間数テンプレートデータの補間数データのコピーを、補間数データＤ２０５として作成する。
【００５１】
ユーザにより上記選択肢（ｄ）が選択された場合、補間数決定部２０１は指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４のそれぞれについて相違の程度を示す指標を算出し、算出した指標に応じて各基準周波数ポイント間の補間周波数ポイントの数を以下のように決定する。
【００５２】
まず、補間数決定部２０１は、音源名データＤ２００により示される音源名に対応する指向性データＤ２０３（図６参照）のうち、隣り合う基準周波数ポイントに関するものの間の相違の程度を示す指標（以下、「指向性相違指標」と呼ぶ）を算出する。例えば、補間数決定部２０１は、音源名データＤ２００が「スピーカＳＰ１」を示す場合、「スピーカＳＰ１」に対応する３１．３Ｈｚの指向性データＤ２０３と６２．５Ｈｚの指向性データＤ２０３に関し、同じベクトルデータを含む指向性要素データに含まれる音強データが示す値の差の自乗を全ての指向性要素データに関し加算することにより、それらの指向性データＤ２０３の間の指向性相違指標を算出する。補間数決定部２０１は、他の隣り合う基準周波数ポイントの組み合わせに関しても、同様に指向性相違指標を算出する。
【００５３】
続いて、補間数決定部２０１は、各基準周波数ポイントに関する空気吸音率データＤ２０４（図７参照）のうち、温度湿度データＤ２０１により示される温度および湿度を示す温度データおよび湿度データを含む温度湿度別吸音率データの吸音率データを読み出す。補間数決定部２０１は、そのように読み出した吸音率データにより示される吸音率のうち、隣り合う基準周波数ポイントに関するものの間の相違の程度を示す指標（以下、「空気吸音率相違指標」と呼ぶ）を算出する。例えば、補間数決定部２０１は、温度湿度データＤ２０１が「１０℃」および「１０％」を示す場合、３１．３Ｈｚおよび６２．５Ｈｚの空気吸音率データＤ２０４の各々から温度データが「１０℃」であり、湿度データが「１０％」である温度湿度別吸音率データの吸音率データを読み出し、読み出した吸音率データが示す値の差の自乗を空気吸音率相違指標として算出する。補間数決定部２０１は、他の隣り合う基準周波数ポイントの組み合わせに関しても、同様に空気吸音率相違指標を算出する。
【００５４】
続いて、補間数決定部２０１は、隣り合う基準周波数ポイントに関する反射吸音率データＤ９０４（図４参照）の各々から、形状データＤ９０１に含まれる材質データを含む材質入射角別吸音率データを抽出し、抽出した材質入射角別吸音率データに含まれる吸音率データの相違の程度を示す指標（以下、「反射吸音率相違指標」と呼ぶ）を算出する。例えば、補間数決定部２０１は、形状データＤ９０１に「コンクリート」を示す材質データが含まれている場合、３１．３Ｈｚおよび６２．５Ｈｚの反射吸音率データＤ９０４の各々から材質データが「コンクリート」である材質入射角別吸音率データを抽出する。補間数決定部２０１は、３１．３Ｈｚの反射吸音率データＤ９０４から抽出した材質入射角別吸音率データと、６２．５Ｈｚの反射吸音率データＤ９０４から抽出した材質入射角別吸音率データに各々含まれる入射角データの等しいデータ間において、吸音率データにより示される吸音率の差の自乗を算出し、それらを全て加算する。そのようにして算出される値は、「コンクリート」に関する反射吸音率の相違の程度を示す指標である。
【００５５】
補間数決定部２０１は、他の材質データに関しても同様に指標を算出する。さらに、補間数決定部２０１は、形状データＤ９０１に含まれる材質データにより示される材質の各々について、座標データ列により示される境界面の面積を算出する。補間数決定部２０１は、各材質データに関する境界面の面積をウェイトとして、各材質データに関する反射吸音率の相違の程度を示す指標の加重平均値を算出する。そのようにして算出される加重平均値は、３１．３Ｈｚおよび６２．５Ｈｚの反射吸音率データＤ９０４に関する反射吸音率相違指標である。補間数決定部２０１は、他の隣り合う基準周波数ポイントの組み合わせに関しても、同様に反射吸音率相違指標を算出する。
【００５６】
続いて、補間数決定部２０１は隣り合う基準周波数ポイント間のそれぞれに関し、上記のように算出した指向性相違指標、空気吸音率相違指標および反射吸音率相違指標の平均値を算出する。その際、補間数決定部２０１は、指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４のそれぞれが、それらを用いてポイントエコーダイアグラム算定部９０２により算定されるエコーダイアグラムの形状に対し与える影響の大きさに応じたウェイトを用いて、指向性相違指標、空気吸音率相違指標および反射吸音率相違指標の加重平均値を算出するようにしてもよい。このように算出した指標が大きい場合、隣り合う基準周波数ポイントに関する指向性データＤ２０３等に大きな差異があるため、補間数決定部２０１は、その基準周波数ポイント間の補間周波数ポイントの数として大きな数値を決定する。一方、算出した指標が小さい場合、補間数決定部２０１はその基準周波数ポイント間の補間周波数ポイントとして小さな数値を決定する。補間数決定部２０１はそのように決定した数値列を示すデータを補間数データＤ２０５として作成する。
【００５７】
ユーザにより上記選択肢（ｅ）が選択された場合、補間数決定部２０１は音データＤ９０５により示される音の時間波形の全てもしくは一部をフーリエ変換し、その時間波形に含まれる周波数成分を算出する。補間数決定部２０１は、隣り合う基準周波数ポイント間の周波数帯域に含まれる周波数成分の振幅値の大きさに応じて、補間周波数ポイントの数として大きな数値を決定する。補間数決定部２０１は決定した数値列を示すデータを補間数データＤ２０５として作成する。
【００５８】
上記のように補間数決定部２０１により作成された補間数データＤ２０５は一時記憶部９０６に記憶される。補間数データＤ２０５が一時記憶部９０６に記憶されると、補間指向性算定部２０２、補間空気吸音率算定部２０３および補間反射吸音率算定部２０４は、それぞれ指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４により示される基準周波数ポイントにおける数値を、補間数データＤ２０５により示される数だけ補間し、補間指向性データＤ２０６、補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８を作成する。以下、それらの処理を説明する。
【００５９】
補間指向性算定部２０２は、音源名データＤ２００により示される音源名に対応する各々の基準周波数ポイントに関する指向性データＤ２０３（図６参照）から、ベクトルデータが等しい指向性要素データを読み出し、読み出したそれらの指向性要素データに含まれる音強データにより示される数値をｙ座標、対応する基準周波数ポイントをｘ座標とする点を通過するスプライン曲線を算出する。続いて、補間指向性算定部２０２は、算出したスプライン曲線上の点であり、ｘ軸を対数表示した場合に、各基準周波数ポイントを示すｘ座標間を補間数データＤ２０５により示される数で各々内分する点の座標を算出する。
【００６０】
例えば、補間数データＤ２０５により示される最初の数が「３」である場合、補間指向性算定部２０２は、３１．３Ｈｚと６２．５Ｈｚを対数表示において４等分する３点、すなわち３７．２Ｈｚ、４４．２Ｈｚおよび５２．６Ｈｚをｘ座標とするスプライン曲線上の点の座標を求める。これらの点（以下、「補間点」と呼ぶ）のｘ座標が補間周波数ポイントであり、ｙ座標が補間周波数ポイントに関する音強の推定値を示す。補間指向性算定部２０２は、指向性データＤ２０３に含まれる全ての指向性要素データに関し上記のように補間点の座標を算出すると、ｘ座標毎、すなわち補間周波数ポイントごとに、ベクトルデータとｙ座標、すなわち音強の推定値を示す音強データを対応付けたデータの集まりを作成する。そのように作成されるデータが補間指向性データＤ２０６である。補間指向性データＤ２０６の構造は、１つの音源名に対応する指向性データＤ２０３（図６参照）と同様である。補間指向性算定部２０２により作成された補間指向性データＤ２０６は、一時記憶部９０６に記憶される。
【００６１】
ここで、スプライン曲線を用いるのは例示であり、補間指向性算定部２０２は他の種類の曲線や、隣り合う基準周波数ポイントに対応する点を結ぶ直線上の点を補間点として算出するようにしてもよい。この点は、以下の補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８の作成に関しても同様である。
【００６２】
補間空気吸音率算定部２０３は、各基準周波数ポイントに関する空気吸音率データＤ２０４（図７参照）から、温度湿度データＤ２０１により示される温度および湿度を示す温度データおよび湿度データを含む温度湿度別吸音率データを読み出し、読み出したそれらの温度湿度別吸音率データに含まれる吸音率データにより示される数値をｙ座標、対応する基準周波数ポイントをｘ座標とする点を通過するスプライン曲線を算出する。続いて、補間空気吸音率算定部２０３は、上述した補間指向性算定部２０２における場合と同様に、補間数データＤ２０５により示される数の補間点の座標を算出する。そのように算出される座標のｙ座標は補間周波数ポイントに関する空気吸音率の推定値を示す。
【００６３】
補間空気吸音率算定部２０３は、上記のように算出した補間点のｘ座標により示される補間周波数ポイントに関し、補間点のｙ座標を示すデータを補間空気吸音率データＤ２０７として作成する。ここで、各補間周波数ポイントに関する補間空気吸音率データＤ２０７は、温度湿度データＤ２０１により示される温度および湿度における吸音率の推定値を示す吸音率データを１つだけ含んでいる。補間空気吸音率算定部２０３により作成された補間空気吸音率データＤ２０７は、一時記憶部９０６に記憶される。
【００６４】
補間反射吸音率データＤ２０８は、各基準周波数ポイントに関する反射吸音率データＤ９０４（図４参照）から、形状データＤ９０１に含まれる材質データを含む材質入射角別吸音率データを抽出する。補間反射吸音率算定部２０４は、各基準周波数ポイントに関し、抽出した材質入射角別吸音率データから、材質データおよび入射角データが等しいものを読み出し、読み出した材質入射角別吸音率データに含まれる吸音率データにより示される数値をｙ座標、対応する基準周波数ポイントをｘ座標とする点を通過するスプライン曲線を算出する。続いて、補間反射吸音率算定部２０４は、上述した補間指向性算定部２０２における場合と同様に、補間数データＤ２０５により示される数の補間点の座標を算出する。そのように算出される補間点のｙ座標は、補間周波数ポイントの音を材質データにより示される材質の境界面に対し入射角データにより示される入射角で反射させた場合の反射吸音率の推定値を示す。
【００６５】
補間反射吸音率算定部２０４は、先に抽出した全ての材質入射角別吸音率データに関し、上記のように補間点の座標を算出すると、ｘ座標毎、すなわち補間周波数ポイントごとに、材質データと入射角データをｙ座標、すなわち反射吸音率の推定値を示す吸音率データに対応付けたデータの集まりを作成する。そのように作成されるデータが補間反射吸音率データＤ２０８である。補間反射吸音率データＤ２０８の構造は、反射吸音率データＤ９０４（図４参照）と同様である。ただし、補間反射吸音率データＤ２０８には、形状データＤ９０１に含まれる材質データに関する材質入射角別吸音率データしか含まれていない。補間反射吸音率算定部２０４により作成された補間反射吸音率データＤ２０８は、一時記憶部９０６に記憶される。以上が、音場シミュレータ１０における一連の補間処理である。
【００６６】
以上のように、補間指向性データＤ２０６、補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８が一時記憶部９０６に記憶されると、ポイントエコーダイアグラム算定部９０２は、経路データ群Ｄ９０６、音源名データＤ２００および温度湿度データＤ２０１と、各基準周波数ポイントに関する指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４を用いて、各基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７を作成する。以下、ポイントエコーダイアグラム算定部９０２がポイントエコーダイアグラムデータＤ９０７を作成する手順を説明する。
【００６７】
まず、ポイントエコーダイアグラム算定部９０２は、経路データ群Ｄ９０６により示される音の経路の各々に関し、音源位置から発せられる方向に対応する音強を、指向性データＤ２０３（図６参照）に基づき決定する。そのため、ポイントエコーダイアグラム算定部９０２は、指向性データＤ２０３において基準となる方向を示すベクトル、すなわち（０，１，０）と、音源位置データＤ９０２に含まれるベクトルデータにより示されるベクトルとの差ベクトルを求める。続いて、ポイントエコーダイアグラム算定部９０２は、経路データ群Ｄ９０６に含まれる各経路データに関し、経路データが含む最初の座標を始点とし、２番目の座標を終点とする方向ベクトルの単位ベクトルを算出する。さらに、ポイントエコーダイアグラム算定部９０２は、そのように算出した単位ベクトルに、先に求めた差ベクトルを加算する。その結果得られるベクトルは、音源位置データＤ９０２により示される音の発せられる方向を、指向性データＤ２０３における基準となる方向と一致させた場合に、経路データにより示される経路をたどる音が音源位置から発せられる際の方向を示す。従って、ポイントエコーダイアグラム算定部９０２は、そのように算出したベクトルを示すベクトルデータを含む指向性要素データを、音源名データＤ２００により示される音源名に対応する指向性データＤ２０３から検索する。このようにして検索された指向性要素データに含まれる音強データにより示される数値が、その経路における音強の初期値として用いられる。以下、この音強の初期値を「音強値１」と呼ぶ。
【００６８】
続いて、ポイントエコーダイアグラム算定部９０２は、空気吸音率データＤ２０４（図７参照）から、温度湿度データＤ２０１により示される温度および湿度を示す温度データおよび湿度データを含む温度湿度別吸音率データを読み出し、読み出した温度湿度別吸音率データに含まれる吸音率データにより示される数値を、経路データにより示される経路の総延長を示す数値に乗じる。その乗算により得られる数値は、総延長を音が伝播した場合に空気により減衰される音強の比率を示す。従って、ポイントエコーダイアグラム算定部９０２は、乗算により得た数値を１から減じた数値を、先に求めた音強値１に乗じる。その乗算により得られる数値は、指向性および空気吸音率のみを考慮した場合の受音位置における音強を示す。以下、この音強を示す値を「音強値２」と呼ぶ。
【００６９】
続いて、ポイントエコーダイアグラム算定部９０２は、経路データにより示される反射位置の各々について、形状データＤ９０１により示される境界面に対する反射位置における入射角の算出およびその境界面の材質の特定を行う。例えば、経路データが「・・・，（６，９，０），（６，０，９），・・・」であり、形状データＤ９０１が図３に示す境界面データを含む場合、（６，０，９）で示される反射位置は図３の先頭行の境界面データの座標データ列により示される境界面上の反射位置を示し、その入射角は４５度と算出される。また、その反射位置に対応する材質は「コンクリート」であると特定される。
【００７０】
続いて、ポイントエコーダイアグラム算定部９０２は反射吸音率データＤ９０４（図４参照）から、特定した材質を示す材質データおよび算出した入射角を示す入射角データを含む材質入射角別吸音率データを検索し、検索した材質入射角別吸音率データに含まれる吸音率データにより示される反射吸音率を特定する。ポイントエコーダイアグラム算定部９０２は、上記のようにして特定した反射吸音率を１から減じる。ポイントエコーダイアグラム算定部９０２は、各々の反射位置に関し、同様に反射吸音率を特定し、それぞれを１から減じる。ポイントエコーダイアグラム算定部９０２は、そのようにして算出した数値の全てを、先に算出した音強値２に乗じる。その乗算により得られる数値は、指向性、空気吸音率および反射吸音率のみを考慮した場合の受音位置における音強を示す。以下、この音強を示す値を「音強値３」と呼ぶ。
【００７１】
続いて、ポイントエコーダイアグラム算定部９０２は経路データにより示される経路の総延長の自乗に４πを乗じた値の逆数を算出し、その逆数を先に算出した、音強値３に乗じる。その乗算により得られる数値は、指向性、空気吸音率、反射吸音率および発散を全て考慮した場合の受音位置における音強を示す。以下、その数値を「減衰後音強値」と呼ぶ。
【００７２】
続いて、ポイントエコーダイアグラム算定部９０２は、温度湿度データＤ２０１により示される温度を用いて音速を算出する。その際、ポイントエコーダイアグラム算定部９０２は、例えば音速の算出式として、「音速（ｍ／秒）＝３３１＋（温度×０．６）」を用いる。ポイントエコーダイアグラム算定部９０２は、経路データにより示される経路の総延長を、算出した音速で除することにより、音源位置において発せられた音が受音位置に到達するまでの時間（以下、「到達時間」と呼ぶ）を算出する。
【００７３】
ポイントエコーダイアグラム算定部９０２は、各基準周波数ポイントに関し、上記のようにして全ての経路データに関し音の到達時間および減衰後音強値を算出すると、音場シミュレータ９０における場合と同様に、それらの数値の組み合わせを示すデータ群をポイントエコーダイアグラムデータＤ９０７として作成する。
【００７４】
続いて、ポイントエコーダイアグラム算定部９０２は、各補間周波数ポイントに関しても、上記の基準周波数ポイントに関する場合と同様の手順によりポイントエコーダイアグラムデータＤ９０７を作成する。ただし、その場合、ポイントエコーダイアグラム算定部９０２は指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４の代わりに、補間指向性データＤ２０６、補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８を用いる。
【００７５】
以上が、ポイントエコーダイアグラム算定部９０２がポイントエコーダイアグラムデータＤ９０７を作成する手順である。ポイントエコーダイアグラム算定部９０２により作成された基準周波数ポイントおよび補間周波数ポイントの各々に関するポイントエコーダイアグラムデータＤ９０７は、一時記憶部９０６に記憶される。
【００７６】
ポイントエコーダイアグラムデータＤ９０７が一時記憶部９０６に記憶されると、音場シミュレータ１０は、ユーザに対し、帯域エコーダイアグラム算定部９０３が帯域エコーダイアグラムデータＤ９０８を作成する際に用いる時間関数に関する選択肢のうち、いずれかを選択するように促す画面を表示する。音場シミュレータ１０は、例えば以下の選択肢を表示する。
（ａ）矩形（クロスフェードなし）
（ｂ）サイン形
（ｃ）ノコギリ形
（ｄ）台形
【００７７】
ユーザが、上記の画面において選択肢（ａ）を選択した場合、帯域エコーダイアグラム算定部９０３は音場シミュレータ９０における場合と同様に、図１４に示したような、基準周波数ポイントもしくは補間周波数ポイントを中心周波数ポイントとする矩形の周波数成分を含む時間波形の係数列を、各々対応するポイントエコーダイアグラムデータＤ９０７に畳み込むことにより、帯域エコーダイアグラムデータＤ９０８を作成する。
【００７８】
ユーザが、上記の画面において選択肢（ｂ）を選択した場合、帯域エコーダイアグラム算定部９０３は、基準周波数ポイントおよび補間周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７により示されるエコーダイアグラムに対し、図８（１）〜（ｎ）に示す周波数成分を有する時間波形の係数列を各々畳み込むことにより、それらの基準周波数ポイントもしくは補間周波数ポイントを中心周波数ポイントとする周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を生成する。ここで、図８（１）〜（ｎ）に示すｆ1〜ｆnは基準周波数ポイントおよび補間周波数ポイントを低域側から順に並べたものである。この点は、以下の図９および図１０に関しても同様である。
【００７９】
図８（１）、すなわち最も低域側の基準周波数ポイントに畳み込まれる時間波形に含まれる周波数成分の振幅値は、例えば以下の関数により与えられる。
【数１０】
2006030689.tif 000012
【００８０】
図８（２）〜（ｎ−１）に示される周波数成分の振幅値は、例えば以下の関数により与えられる。
【数１１】
2006030689.tif 000013
【００８１】
図８（ｎ）、すなわち最も高域側の基準周波数ポイントに畳み込まれる時間波形に含まれる周波数成分の振幅値は、例えば以下の関数により与えられる。
【数１２】
2006030689.tif 000014
【００８２】
ユーザが、上記の画面において選択肢（ｃ）を選択した場合、帯域エコーダイアグラム算定部９０３は、基準周波数ポイントおよび補間周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７により示されるエコーダイアグラムに対し、図９（１）〜（ｎ）に示す周波数成分を有する時間波形の係数列を各々畳み込むことにより、それらの基準周波数ポイントもしくは補間周波数ポイントを中心周波数ポイントとする周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を生成する。
【００８３】
ユーザが、上記の画面において選択肢（ｄ）を選択した場合、帯域エコーダイアグラム算定部９０３は、基準周波数ポイントおよび補間周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７により示されるエコーダイアグラムに対し、図１０（１）〜（ｎ）に示す周波数成分を有する時間波形の係数列を各々畳み込むことにより、それらの基準周波数ポイントもしくは補間周波数ポイントを中心周波数ポイントとする周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を生成する。
【００８４】
ここで、図８、図９および図１０に示した周波数関数は、以下の特徴を共通して有している。
（ａ）全ての周波数帯域に関する周波数関数を加算すると、全周波数域において周波数成分のレベルが「１」になる。
（ｂ）各中心周波数ポイントにおける周波数成分のレベルは「１」である。
（ｃ）隣り合う周波数帯域間で重なり部分をもち、重なり部分において、各々の周波数関数がクロスフェードしている。
【００８５】
ただし、実際には丸め誤差等により、上記（ａ）および（ｂ）の条件に関し、全周波数域における加算後の周波数成分のレベルが「１」より多少上下したり、また各中心周波数ポイントにおける周波数成分のレベルが「１」より多少上下したりしてもよい。帯域エコーダイアグラム算定部９０３は、上記（ａ）〜（ｃ）の条件を満たす他の形状の周波数関数を、図８〜図１０に示した周波数関数の代わりに用いるようにしてもよい。さらに、特殊な音響効果を付加する目的で、例えば上記（ａ）や（ｂ）の条件を満たさない形状の周波数関数を用いるようにしてもよい。以上のようにして帯域エコーダイアグラム算定部９０３により作成された帯域エコーダイアグラムデータＤ９０８は、一時記憶部９０６に記憶される。
【００８６】
加算部９０４は、一時記憶部９０６に記憶された帯域エコーダイアグラムデータＤ９０８の全てを加算することにより、インパルス応答データＤ９０９を作成する。加算部９０４により生成されたインパルス応答データＤ９０９は、一時記憶部９０６に記憶される。
【００８７】
畳み込み演算部９０５は、音データＤ９０５に対し、インパルス応答データＤ９０９により示される係数列を畳み込むことにより、インパルス応答データＤ９０９により示される音場情報の付加された音を示す音場情報付音データＤ９１０を生成する。音場情報付音データＤ９１０を再生した場合に得られる音は、実際に形状データＤ９０１により示される音場内の音源位置データＤ９０２により示される音源位置において、音データＤ９０５に従った音を再生した場合に、受音位置データＤ９０３により示される受音位置において受音される音と類似した音となる。
【００８８】
ユーザが、上述した補間数の決定方法に関する選択肢のうち、「個別に指定」を選択し、かつ全ての隣り合う基準周波数ポイント間の補間数として「０」を選択した場合を除き、音場シミュレータ１０は、補間処理により、音場シミュレータ９０における場合と比較して、より狭い周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を接合したインパルス応答データＤ９０９を生成する。従って、そのインパルス応答データＤ９０９に含まれる周波数成分の周波数軸上のギャップは音場シミュレータ９０におけるものと比較して小さくなる。その結果、音場シミュレータ１０により作成される音場情報付音データＤ９１０を再生して得られる音は、より自然な聴感をもたらすものとなる。
【００８９】
また、ユーザが、上述した帯域エコーダイアグラムデータＤ９０８を作成する際に用いられる時間関数に関する選択肢のうち、「矩形（クロスフェードなし）」を選択した場合を除き、音場シミュレータ１０は、隣り合う帯域エコーダイアグラムデータＤ９０８に含まれる周波数成分のレベルが互いにクロスフェードするように、帯域エコーダイアグラムデータＤ９０８を作成する。従って、この場合にも、帯域エコーダイアグラムデータＤ９０８を加算して得られるインパルス応答データＤ９０９に含まれる周波数成分の周波数軸上のギャップは音場シミュレータ９０におけるものと比較して小さくなる。その結果、音場シミュレータ１０により作成される音場情報付音データＤ９１０を再生して得られる音は、より自然な聴感をもたらすものとなる。
【００９０】
ところで、本実施形態においては、音場シミュレータ１０は、反射吸音率データＤ９０４として材質および入射角を考慮した反射吸音率を示すデータを用い、さらに指向性データＤ２０３および空気吸音率データＤ２０４を用いることにより、音場シミュレータ９０における場合と比較して、より詳細なシミュレーションを行っている。しかしながら、音場シミュレータ１０は、音場シミュレータ９０と同様に、材質もしくは入射角の両方もしくは一方を考慮しない反射吸音率を示す反射吸音率データＤ９０４を用いたり、指向性データＤ２０３および空気吸音率データＤ２０４の両方もしくは一方を用いないで、上述した補間処理およびクロスフェード処理を行うようにしてもよい。
【００９１】
また、本実施形態においては、音場シミュレータ１０は、ポイントエコーダイアグラムデータＤ９０７に対し時間関数の係数列を畳み込むことにより、帯域エコーダイアグラムデータＤ９０８を作成するものとして説明したが、音場シミュレータ１０は、時間関数であるポイントエコーダイアグラムデータＤ９０７を、フーリエ変換等により周波数関数に変換し、その周波数関数と、図８〜図１０等に示される周波数関数とを乗算して得られる周波数関数を逆フーリエ変換等により時間関数に変換することにより、帯域エコーダイアグラムデータＤ９０８を作成するようにしてもよい。このように、時間関数による演算は、周波数関数による等価な演算を有するため、上述した一連の処理における時間関数間の演算処理の一部もしくは全てを、周波数関数間の等価な演算処理と置き換えてもよい。
【００９２】
また、ポイントエコーダイアグラムデータＤ９０７を用いて帯域エコーダイアグラムデータＤ９０８を作成する処理を、既存のフィルタ装置を用いて行うようにしてもよい。
【００９３】
［２．第２実施形態］
本発明の第２実施形態にかかるデータ処理装置（以下、「音場シミュレータ３０」と呼ぶ）は、音場シミュレータ１０と同様に、音データに対し音場情報を付加するためのインパルス応答データを生成し、生成したインパルス応答データを音データに畳み込む処理を行うことにより、音場情報の付加された音を示す音データを生成する装置である。音場シミュレータ３０は、第１実施形態における音場シミュレータ１０と多くの点で共通しているので、以下、音場シミュレータ３０が音場シミュレータ１０と異なる点を中心に説明する。
【００９４】
図１１は、音場シミュレータ３０の機能構成を示した図である。まず、本実施形態において経路算定部９０１により経路データ群Ｄ９０６が作成される手順や、補間数決定部２０１により補間数データＤ２０５が作成される手順は第１実施形態における場合と同様である。
【００９５】
次に補間処理についてであるが、本実施形態における処理は第１実施形態における処理と異なっている。第１実施形態の音場シミュレータ１０においては、補間数データＤ２０５に基づき、補間指向性算定部２０２、補間空気吸音率算定部２０３および補間反射吸音率算定部２０４が、それぞれ指向性データＤ２０３、空気吸音率データＤ２０４および反射吸音率データＤ９０４を補間するデータとして、補間指向性データＤ２０６、補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８を作成する。これに対し、本実施形態の音場シミュレータ３０は、補間指向性算定部２０２、補間空気吸音率算定部２０３および補間反射吸音率算定部２０４は備えておらず、従って補間指向性データＤ２０６、補間空気吸音率データＤ２０７および補間反射吸音率データＤ２０８は作成されない。すなわち、ポイントエコーダイアグラム算定部９０２は、基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７のみを作成する。
【００９６】
その一方で、本実施形態の音場シミュレータ３０は、補間ポイントエコーダイアグラム算定部３０１を備えている。補間ポイントエコーダイアグラム算定部３０１は、基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７から、補間周波数ポイントに関するエコーダイアグラムを示す補間ポイントエコーダイアグラムデータＤ３０１を生成する構成部である。以下、補間ポイントエコーダイアグラム算定部３０１の処理を説明する。
【００９７】
補間ポイントエコーダイアグラム算定部３０１は、基準周波数ポイントに関するポイントエコーダイアグラムデータＤ９０７が一時記憶部９０６に記憶されると、それらのポイントエコーダイアグラムデータＤ９０７に含まれる、同じ時間に関する音強のレベルを示すデータを読み出す。続いて、補間ポイントエコーダイアグラム算定部３０１は、読み出したデータにより示される音強を示す数値をｙ座標、対応する基準周波数ポイントをｘ座標とする点を全て通過するスプライン曲線を算出する。続いて、補間ポイントエコーダイアグラム算定部３０１は、算出したスプライン曲線上の点であり、ｘ軸を対数表示した場合に、各基準周波数ポイントを示すｘ座標間を補間数データＤ２０５により示される数で各々内分する点の座標を算出する。
【００９８】
例えば、補間数データＤ２０５により示される最初の数が「３」であり、各基準周波数ポイント、すなわち３１．３Ｈｚ、６２．５Ｈｚ、１２５Ｈｚ、・・・のポイントエコーダイアグラムデータＤ９０７が、０．２５ミリ秒における音強のレベルとして、それぞれ「０．８９０」、「０．９３１」、「０．８４２」、・・・を示す場合、補間ポイントエコーダイアグラム算定部３０１はまず、（３１．３，０．８９０）、（６２．５，０．９３１）、（１２５，０．８４２）、・・・の各座標が示す点を全て通過するスプライン曲線を算出する。続いて、補間ポイントエコーダイアグラム算定部３０１は、算出したスプライン曲線上の点であり、対数表示において３１．３と６２．５を４等分する３つの数値、すなわち３７．２、４４．２および５２．６をｘ座標とする点のｙ座標を求める。これらのｙ座標が、補間周波数ポイント、すなわち３７．２Ｈｚ、４４．２Ｈｚおよび５２．６Ｈｚに関する０．２５ミリ秒における音強の推定値を示す。
【００９９】
補間ポイントエコーダイアグラム算定部３０１は、同様に０．２５ミリ秒における他の補間周波数ポイントに関する音強の推定値を算出する。さらに、補間ポイントエコーダイアグラム算定部３０１は、同様に、各補間周波数ポイントに関し他の時間における音強の推定値も算出する。補間ポイントエコーダイアグラム算定部３０１は、そのように算出した各補間周波数ポイントに関する、各時間における音強の推定値を示すデータを、各補間周波数ポイントに関する補間ポイントエコーダイアグラムデータＤ３０１として作成する。以上が、補間ポイントエコーダイアグラム算定部３０１により、補間ポイントエコーダイアグラム算定部３０１が作成される場合の処理である。補間ポイントエコーダイアグラム算定部３０１により作成された補間ポイントエコーダイアグラムデータＤ３０１は、一時記憶部９０６に記憶される。
【０１００】
補間ポイントエコーダイアグラムデータＤ３０１が一時記憶部９０６に記憶されると、音場シミュレータ９０における場合と同様に、帯域エコーダイアグラム算定部９０３は帯域エコーダイアグラムデータＤ９０８を作成する。ただし、本実施形態においては、帯域エコーダイアグラム算定部９０３は各補間周波数ポイントを中心周波数ポイントとする周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を作成する際には、補間ポイントエコーダイアグラム算定部３０１により作成された補間ポイントエコーダイアグラムデータＤ３０１を用いる。その結果、第１実施形態における場合と同様に、各基準周波数ポイントおよび各補間周波数ポイントを中心周波数ポイントとする周波数帯域に関し、帯域エコーダイアグラムデータＤ９０８が作成される。
【０１０１】
帯域エコーダイアグラムデータＤ９０８を用いて、加算部９０４によりインパルス応答データＤ９０９が作成される手順、音データＤ９０５およびインパルス応答データＤ９０９を用いて、畳み込み演算部９０５により音場情報付音データＤ９１０が作成される手順は、第１実施形態における場合と同様である。
【０１０２】
上記のように、音場シミュレータ３０により生成されるインパルス応答データＤ９０９も、音場シミュレータ１０により生成されるインパルス応答データＤ９０９と同様に、音場シミュレータ９０における場合と比較して、より狭い周波数帯域に関する帯域エコーダイアグラムデータＤ９０８を接合したインパルス応答データＤ９０９となる。従って、そのインパルス応答データＤ９０９に含まれる周波数成分の周波数軸上のギャップは音場シミュレータ９０におけるものと比較して小さくなる。その結果、音場シミュレータ１０により作成される音場情報付音データＤ９１０を再生して得られる音は、より自然な聴感をもたらすものとなる。
【図面の簡単な説明】
【０１０３】
【図１】本発明の実施形態にかかる音場シミュレータのハードウェア構成を示した図である。
【図２】本発明の第１実施形態にかかる音場シミュレータの機能構成を示したブロック図である。
【図３】本発明の実施形態にかかる形状データの構造を示した図である。
【図４】本発明の実施形態にかかる反射吸音率データの構造を示した図である。
【図５】本発明の実施形態にかかる補間数テンプレートデータ群の構造を示した図である。
【図６】本発明の実施形態にかかる指向性データの構造を示した図である。
【図７】本発明の実施形態にかかる空気吸音率データの構造を示した図である。
【図８】本発明の実施形態にかかる帯域エコーダイアグラムデータを作成する際に用いられる時間波形の周波数成分を示した図である。
【図９】本発明の実施形態にかかる帯域エコーダイアグラムデータを作成する際に用いられる時間波形の周波数成分を示した図である。
【図１０】本発明の実施形態にかかる帯域エコーダイアグラムデータを作成する際に用いられる時間波形の周波数成分を示した図である。
【図１１】本発明の第２実施形態にかかる音場シミュレータの機能構成を示したブロック図である。
【図１２】従来技術にかかる音場シミュレータの機能構成を示した図である。
【図１３】従来技術にかかるエコーダイアグラムデータを例示した図である。
【図１４】従来技術にかかる帯域エコーダイアグラムデータを作成する際に用いられる時間波形の周波数成分を示した図である。
【符号の説明】
【０１０４】
１０・３０・９０…音場シミュレータ、１０１…ＣＰＵ、１０２…ＲＯＭ、１０３…ＲＡＭ、１０４…ＨＤ、１０５…ディスプレイ、１０６…操作部、１０７…サウンドシステム、１０８…入出力Ｉ／Ｆ、２０１…補間数決定部、２０２…補間指向性算定部、２０３…補間空気吸音率算定部、２０４…補間反射吸音率算定部、３０１…補間ポイントエコーダイアグラム算定部、９００…記憶部、９０１…経路算定部、９０２…ポイントエコーダイアグラム算定部、９０３…帯域エコーダイアグラム算定部、９０４…加算部、９０５…畳み込み演算部、９０６…一時記憶部、Ｄ２００…音源名データ、Ｄ２０１…温度湿度データ、Ｄ２０２…補間数テンプレートデータ群、Ｄ２０３…指向性データ、Ｄ２０４…空気吸音率データ、Ｄ２０５…補間数データ、Ｄ２０６…補間指向性データ、Ｄ２０７…補間空気吸音率データ、Ｄ２０８…補間反射吸音率データ、Ｄ３０１…補間ポイントエコーダイアグラムデータ、Ｄ９０１…形状データ、Ｄ９０２…音源位置データ、Ｄ９０３…受音位置データ、Ｄ９０４…反射吸音率データ、Ｄ９０５…音データ、Ｄ９０６…経路データ群、Ｄ９０７…ポイントエコーダイアグラムデータ、Ｄ９０８…帯域エコーダイアグラムデータ、Ｄ９０９…インパルス応答データ、Ｄ９１０…音場情報付音データ。
【図１】
2006030689.tif 000015
【図２】
2006030689.tif 000016
【図３】
2006030689.tif 000017
【図４】
2006030689.tif 000018
【図５】
2006030689.tif 000019
【図６】
2006030689.tif 000020
【図７】
2006030689.tif 000021
【図８】
2006030689.tif 000022
【図９】
2006030689.tif 000023
【図１０】
2006030689.tif 000024
【図１１】
2006030689.tif 000025
【図１２】
2006030689.tif 000026
【図１３】
2006030689.tif 000027
【図１４】
2006030689.tif 000028