JP5050789B2 - シミュレーション装置およびプログラム - Google Patents
シミュレーション装置およびプログラム Download PDFInfo
- Publication number
- JP5050789B2 JP5050789B2 JP2007290790A JP2007290790A JP5050789B2 JP 5050789 B2 JP5050789 B2 JP 5050789B2 JP 2007290790 A JP2007290790 A JP 2007290790A JP 2007290790 A JP2007290790 A JP 2007290790A JP 5050789 B2 JP5050789 B2 JP 5050789B2
- Authority
- JP
- Japan
- Prior art keywords
- simulation
- sound
- virtual
- sounding body
- existing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10H—ELECTROPHONIC MUSICAL INSTRUMENTS; INSTRUMENTS IN WHICH THE TONES ARE GENERATED BY ELECTROMECHANICAL MEANS OR ELECTRONIC GENERATORS, OR IN WHICH THE TONES ARE SYNTHESISED FROM A DATA STORE
- G10H1/00—Details of electrophonic musical instruments
- G10H1/0008—Associated control or indicating means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、仮想発音体が発する音を解析・可聴化する技術に関する。
従来、波動音響シミュレーションは、その対象物、対象空間、拘束条件などが単純なモデルの場合（たとえば金属板や木板を打撃した音など）は、実物の現象を精度良く再現可能であった。しかし、対象物が複雑な条件を有する発音体である場合、シミュレーションの解の精度は悪く、定性的な音響特性は類似していても定量的には実物の現象との差異が大きかった。
以上に述べた問題点を解消するために、以下に述べるような技術が提案されている。
特許文献１においては、音響特性解析において、成形品の局所的な物性データを計算し、該物性データを局所的な領域ごとに与えて計算することにより、精度良く音響特性を解析する技術が開示されている。
特許文献２においては、振動・音響解析の音圧スペクトルデータを時系列波形に変換し、その波形をスピーカで再生して評価する技術が開示されている。
特許文献３においては、比較的大規模な音場の音場解析において、防吸音部材の物理特性を理論値や実測値を用いて計算する技術が開示されている。
特開２００３−０９０７５８号公報
特開２００５−３０８７２６号公報
特開２００６−０６５４６６号公報
特許文献１においては、音響特性解析において、成形品の局所的な物性データを計算し、該物性データを局所的な領域ごとに与えて計算することにより、精度良く音響特性を解析する技術が開示されている。
特許文献２においては、振動・音響解析の音圧スペクトルデータを時系列波形に変換し、その波形をスピーカで再生して評価する技術が開示されている。
特許文献３においては、比較的大規模な音場の音場解析において、防吸音部材の物理特性を理論値や実測値を用いて計算する技術が開示されている。
上記特許文献１および２の技術においては、解析で設定する入力条件のうち、構造体の拘束条件、支持条件、結合条件、吸音条件などの境界条件、音源や振動体の加振条件、構造体の振動減衰係数などを正確に与えることは難しく、これらの誤差がシミュレーションの精度を下げるとの問題点があった。
また、このシミュレーションの解をそのまま用いてスピーカ再生した場合、実物が発する音との聴感的な差異が大きく、適切な音の評価ができなかった。
また、シミュレーションにおいて精度の良い周波数帯域のみのデータを用いてスピーカ再生した場合、すなわち帯域を限定して音再生をした場合も、実物が発する音と聴感的な差異が大きく、適切な音の評価ができなかった。
また、特許文献３の技術においては、理論値の利用では、理論モデルの妥当性や解析対象モデルとの整合性、実測値の利用では実音場ではなく材料単体の理想的な条件での実測、等の点から解析精度に限界があった。
また、シミュレーションにおいて精度の良い周波数帯域のみのデータを用いてスピーカ再生した場合、すなわち帯域を限定して音再生をした場合も、実物が発する音と聴感的な差異が大きく、適切な音の評価ができなかった。
また、特許文献３の技術においては、理論値の利用では、理論モデルの妥当性や解析対象モデルとの整合性、実測値の利用では実音場ではなく材料単体の理想的な条件での実測、等の点から解析精度に限界があった。
本発明は上述の課題に鑑みてなされたものであり、対象が複雑な系のモデルであっても、発音特性が精度良くシミュレーションできるシミュレーション装置およびプログラムを提供することを目的とする。
本発明に係るシミュレーション装置は、特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、前記既存発音体シミュレーション手段と前記仮想発音体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第１の差異データ生成手段と、前記差異データを用いて、前記既存発音体測定手段による測定結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第１の特性補正手段とを具備することを特徴とする。
上記の構成において、前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、周波数応答に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定しても良い。また、前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、位相特性に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定しても良い。
上記の構成において、前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、周波数応答に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定しても良い。また、前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、位相特性に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定しても良い。
本発明に係る別の態様のシミュレーション装置は、特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第２の差異データ生成手段と、前記差異データを用いて、前記仮想発音体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第２の特性補正手段とを具備することを特徴とする。
上記の構成において、前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、周波数応答に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定しても良い。また、前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、位相特性に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定しても良い。
上記の構成において、前記特定の発音体の演奏内容を表す音データを受取る受取手段と、前記受取手段が受取った音データを、前記仮想発音体予測データに基づいて変形し出力する音データ変形手段とを具備していても良い。
上記の構成において、前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、周波数応答に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定しても良い。また、前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、位相特性に関する差異を算出し、前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定しても良い。
上記の構成において、前記特定の発音体の演奏内容を表す音データを受取る受取手段と、前記受取手段が受取った音データを、前記仮想発音体予測データに基づいて変形し出力する音データ変形手段とを具備していても良い。
本発明に係るシミュレーション装置は、上記の構成において、特定の構造体に係る音の伝達特性をシミュレートする伝達特性シミュレート手段と、前記伝達特性シミュレート手段によるシミュレーション結果を用いて前記仮想発音体予測データを補正して、前記仮想の発音体が発音した音が前記特定の構造体を伝達した結果の複合体発音特性を予測する第１の予測データ補正手段とを具備していても良い。また、特定の構造体に係る音の伝達特性を測定する伝達特性測定手段と、前記伝達特性測定手段による測定結果を用いて前記仮想発音体予測データを補正して、前記仮想の発音体が発音した音が前記特定の構造体を伝達した結果の複合体発音特性を予測する第２の予測データ補正手段とを具備していても良い。
本発明に係るシミュレーション装置は、上記の構成において、特定の構造体に係る音の伝達特性をシミュレートする既存構造体シミュレーション手段と、仮想の構造体に係る音の伝達特性をシミュレートする仮想構造体シミュレーション手段と、前記特定の構造体が音を伝達する際の伝達特性を測定する既存構造体測定手段と、前記既存構造体シミュレーション手段と前記仮想構造体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第３の差異データ生成手段と、前記差異データを用いて、前記既存構造体測定手段による測定結果を補正し、前記仮想の構造体による音の伝達特性である仮想構造体予測データを生成する第３の特性補正手段と、前記仮想発音体予測データを、前記第３の特性補正手段が生成した仮想構造体予測データに基づいて補正して、前記仮想の発音体が発音した音が前記仮想の構造体を伝達した結果の複合体発音特性を予測する第３の予測データ補正手段とを具備していても良い。また、特定の構造体に係る音の伝達特性をシミュレートする既存構造体シミュレーション手段と、仮想の構造体に係る音の伝達特性をシミュレートする仮想構造体シミュレーション手段と、前記特定の構造体が音を伝達する際の伝達特性を測定する既存構造体測定手段と、前記既存構造体シミュレーション手段によるシミュレーション結果と前記既存構造体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第４の差異データ生成手段と、前記差異データを用いて、前記仮想構造体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の構造体による音の伝達特性である仮想構造体予測データを生成する第４の特性補正手段と、前記仮想発音体予測データを、前記第４の特性補正手段が生成した仮想構造体予測データに基づいて補正して、前記仮想の発音体が発音した音が前記仮想の構造体を伝達した結果の複合体発音特性を予測する第４の予測データ補正手段とを具備していても良い。
本発明に係るプログラムは、コンピュータを、特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、前記既存発音体シミュレーション手段と前記仮想発音体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第１の差異データ生成手段と、前記差異データを用いて、前記既存発音体測定手段による測定結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第１の特性補正手段として機能させることを特徴とする。
本発明に係るプログラムの別の態様は、コンピュータを、特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第２の差異データ生成手段と、前記差異データを用いて、前記仮想発音体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第２の特性補正手段として機能させることを特徴とする。
本発明に係るプログラムの別の態様は、コンピュータを、特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第２の差異データ生成手段と、前記差異データを用いて、前記仮想発音体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第２の特性補正手段として機能させることを特徴とする。
本発明に係るシミュレーション装置およびプログラムによれば、発音特性が精度良くシミュレーションできる。
以下、図面を参照しつつ、本発明を実施する際の最良の形態について説明する。
（Ａ；構成）
図１は、本発明に係る仮想音生成装置１の全体構成を示す図である。仮想音生成装置１は、例えばパーソナルコンピュータで構成される。仮想音生成装置１は、制御部１１、ＲＯＭ（Read Only Memory）１２、ＲＡＭ（Random Access Memory）１３、操作部１４、表示部１５、ＨＤＤ１６（Hard Disk Drive）、音声処理部１７、および音声処理部１７に接続されたスピーカ１８およびマイクロホン１９を有する。上記各部はバスを介して互いに接続されている。
（Ａ；構成）
図１は、本発明に係る仮想音生成装置１の全体構成を示す図である。仮想音生成装置１は、例えばパーソナルコンピュータで構成される。仮想音生成装置１は、制御部１１、ＲＯＭ（Read Only Memory）１２、ＲＡＭ（Random Access Memory）１３、操作部１４、表示部１５、ＨＤＤ１６（Hard Disk Drive）、音声処理部１７、および音声処理部１７に接続されたスピーカ１８およびマイクロホン１９を有する。上記各部はバスを介して互いに接続されている。
制御部１１は、例えばＣＰＵ（Central Processing Unit）であり、ＲＯＭ１２に格納された制御プログラムを実行することにより各部の制御を行う。
ＲＯＭ１２は、制御部１１が実行する制御プログラムを格納している。
ＲＡＭ１３は、制御部１１によってワークエリアとして利用される。
ＲＯＭ１２は、制御部１１が実行する制御プログラムを格納している。
ＲＡＭ１３は、制御部１１によってワークエリアとして利用される。
操作部１４は、例えばキーボードやマウスであり各種の操作子を備える。操作部１４はユーザによる操作内容を表す操作信号を制御部１１へ出力する。
表示部１５は、例えばＬＣＤ（Liquid Crystal Display）などの表示手段である。このディスプレイの画面には、音データの処理に係る各種設定を行うための操作画面が表示される。
ＨＤＤ１６は、大容量の記憶装置である。
表示部１５は、例えばＬＣＤ（Liquid Crystal Display）などの表示手段である。このディスプレイの画面には、音データの処理に係る各種設定を行うための操作画面が表示される。
ＨＤＤ１６は、大容量の記憶装置である。
音声処理部１７は、Ｄ／Ａ（Digital/Analog）変換器とＡ／Ｄ変換器とアンプとからなる。マイクロホン１９から入力された音を表すアナログ信号は、Ａ／Ｄ変換器によりデジタルデータに変換され、制御部１１に出力される。また、制御部１１から受取った音を表すデジタルデータはＤ／Ａ変換によりアナログ信号に変換され、アンプにより振幅が調整された後、スピーカ１８に出力される。
スピーカ１８は、音声処理部１７から受取ったアナログ信号に基づいて音声を放音する。
マイクロホン１９は、音を表すアナログ信号を出力する。
以上が仮想音生成装置１の構成である。
スピーカ１８は、音声処理部１７から受取ったアナログ信号に基づいて音声を放音する。
マイクロホン１９は、音を表すアナログ信号を出力する。
以上が仮想音生成装置１の構成である。
（Ｂ；動作）
次に、上記の構成を有する仮想音生成装置１の動作について説明する。
次に、上記の構成を有する仮想音生成装置１の動作について説明する。
（Ｂ−１；処理の概要）
処理内容の詳細な説明に入る前に、処理の概要について説明する。
本発明に係る仮想音生成装置１は、「実在しない発音体」がどのような音を発音するか、実在する発音体との関連性に基づいてシミュレーションする装置である。後述する動作例においては、現行モデルのキーボード（以下、既存品）に基づいて、開発中のキーボード（以下、開発品）の発音を予測する場合を例にとって説明する。
処理内容の詳細な説明に入る前に、処理の概要について説明する。
本発明に係る仮想音生成装置１は、「実在しない発音体」がどのような音を発音するか、実在する発音体との関連性に基づいてシミュレーションする装置である。後述する動作例においては、現行モデルのキーボード（以下、既存品）に基づいて、開発中のキーボード（以下、開発品）の発音を予測する場合を例にとって説明する。
図２は、本発明に係る仮想音生成装置１が実行する処理の流れを概念的に示した図である。図中、シミュレーション手段１１０、解析手段１１１、差分算出手段１１２、補正手段１１３、合成手段１１４、変形手段１１５は、制御部１１により実現される手段である。以下では、図中の各処理過程に付されたアルファベット（ａ〜ｉ）を参照しながら説明する。
まず、従来のシミュレーションを実行するシミュレーション手段１１０により、既存品Ｅおよび開発品Ｄの発音特性（周波数応答）をシミュレートし（ａ、ｂ）、それぞれのシミュレーション結果を表す既存品シミュレーションスペクトルおよび開発品シミュレーションスペクトルを生成する。そして、差分算出手段１１２により、既存品シミュレーションスペクトルと開発品シミュレーションスペクトルの差分を表す比率データを算出する（ｃ）。一方、既存品Ｅを実際に発音させ（ｄ）、解析手段１１１により該実際の発音を解析することにより既存品の発音特性を解析し（ｅ）、解析結果である既存品実測スペクトルを生成する。次に、補正手段１１３において、上記比率データに基づいて、上記既存品実測スペクトルを補正することにより、新たなスペクトル（特性予測スペクトル）を生成する（ｆ）。該生成された特性予測スペクトルは、既存品Ｅによる実際の発音に基づく発音特性に対して既存品Ｅと開発品Ｄのシミュレーション結果の差分が加味されていることから、開発品Ｄの発音特性を精度良く予測したデータとなる。合成手段１１４において、上記特性予測スペクトルと、既存品実測スペクトルの一部を合成し、合成スペクトルを生成する（ｇ）。そして変形手段１１５において、既存品Ｅにより発音された音（音データ）（ｈ）を、上記生成された合成スペクトルに基づいて変形することで（ｉ）、開発品Ｄの音を模した予測音データを生成する。
まず、従来のシミュレーションを実行するシミュレーション手段１１０により、既存品Ｅおよび開発品Ｄの発音特性（周波数応答）をシミュレートし（ａ、ｂ）、それぞれのシミュレーション結果を表す既存品シミュレーションスペクトルおよび開発品シミュレーションスペクトルを生成する。そして、差分算出手段１１２により、既存品シミュレーションスペクトルと開発品シミュレーションスペクトルの差分を表す比率データを算出する（ｃ）。一方、既存品Ｅを実際に発音させ（ｄ）、解析手段１１１により該実際の発音を解析することにより既存品の発音特性を解析し（ｅ）、解析結果である既存品実測スペクトルを生成する。次に、補正手段１１３において、上記比率データに基づいて、上記既存品実測スペクトルを補正することにより、新たなスペクトル（特性予測スペクトル）を生成する（ｆ）。該生成された特性予測スペクトルは、既存品Ｅによる実際の発音に基づく発音特性に対して既存品Ｅと開発品Ｄのシミュレーション結果の差分が加味されていることから、開発品Ｄの発音特性を精度良く予測したデータとなる。合成手段１１４において、上記特性予測スペクトルと、既存品実測スペクトルの一部を合成し、合成スペクトルを生成する（ｇ）。そして変形手段１１５において、既存品Ｅにより発音された音（音データ）（ｈ）を、上記生成された合成スペクトルに基づいて変形することで（ｉ）、開発品Ｄの音を模した予測音データを生成する。
（Ｂ−２；処理の詳細）
以下では、仮想音生成装置１が実行する処理の詳細を説明する。図３は、仮想音生成装置１の制御部１１が行う処理の流れを示したフローチャートである。
以下では、仮想音生成装置１が実行する処理の詳細を説明する。図３は、仮想音生成装置１の制御部１１が行う処理の流れを示したフローチャートである。
ステップＳＡ１０において、制御部１１は、既存品について、従来のシミュレーションにより発音特性（周波数応答）をシミュレートする。ここで、従来のシミュレーションとは、各種の条件（構造体の拘束条件、支持条件、結合条件、吸音条件などの境界条件、音源や振動体の加振条件、構造体の振動減衰係数など）を入力することにより、既存品による発音の周波数応答を算出するものである。具体的には、例えばホワイトノイズ（すべての周波数成分で同じ強度であるノイズ）の音データを既存品に発音させた場合に、出力される音の周波数帯域ごとの強度を示す周波数スペクトルを生成する。
図４は、既存品について、従来のシミュレーションにより生成された周波数スペクトル（以下、既存品シミュレーションスペクトル）を示す。なお、本実施形態においては、周波数が３００Ｈｚ以下の成分についてシミュレートする。
なお、シミュレーションの対象は、キーボードのスピーカ部分で、寸法が３００ｍｍ（幅）×高さ１００ｍｍ（高さ）×１５０ｍｍ（奥行）である。採用したシミュレーション手法は、要素数１１０５、節点数１０８９とした境界要素法である。
図４は、既存品について、従来のシミュレーションにより生成された周波数スペクトル（以下、既存品シミュレーションスペクトル）を示す。なお、本実施形態においては、周波数が３００Ｈｚ以下の成分についてシミュレートする。
なお、シミュレーションの対象は、キーボードのスピーカ部分で、寸法が３００ｍｍ（幅）×高さ１００ｍｍ（高さ）×１５０ｍｍ（奥行）である。採用したシミュレーション手法は、要素数１１０５、節点数１０８９とした境界要素法である。
ステップＳＡ２０において、制御部１１は、従来のシミュレーションにより開発品の発音特性（周波数応答）をシミュレートする。そのシミュレート方法はステップＳＡ１０と同様である。
図５は、従来のシミュレーションにより開発品について生成された周波数スペクトル（以下、開発品シミュレーションスペクトル）である。なお、開発品シミュレーションスペクトルについても、周波数が３００Ｈｚ以下の成分についてシミュレートする。
図５は、従来のシミュレーションにより開発品について生成された周波数スペクトル（以下、開発品シミュレーションスペクトル）である。なお、開発品シミュレーションスペクトルについても、周波数が３００Ｈｚ以下の成分についてシミュレートする。
図６には、ステップＳＡ１０およびステップＳＡ２０において生成された、既存品および開発品の周波数スペクトルを示す。図６に示したシミュレーション結果から、開発品は、既存品よりも特定の周波数帯域（１００Ｈｚ〜２００Ｈｚと、２１０Ｈｚ〜２８０Ｈｚ）における周波数成分が強調されるとの傾向が予測できる。
ステップＳＡ３０において、制御部１１は、図６に見られたような既存品と開発品との差異を数値化したデータである「比率データ」を生成する。具体的には、周波数帯域ごとに、開発品シミュレーションスペクトルの値を既存品シミュレーションスペクトルの値で除した比の値を算出する。図７は、そのようにして生成された比率データを示す。該比率データにおいて、値が１を超える周波数成分は、開発品において既存品よりも強調されていることを表し、逆に値が１を下回る周波数成分は、既存品において開発品よりも強調されていることを表す。
ステップＳＡ４０において、制御部１１は、既存品に実際にホワイトノイズを入力して発音させ、出力される音の周波数帯域ごとの振幅強度を示す周波数スペクトル（以下、既存品実測スペクトル）を生成する。
具体的には、既存品にホワイトノイズを入力して発音させた音を収音し、該収音された音を表す音データをフーリエ変換し、パワースペクトルを経時的に生成する。そして、該パワースペクトルにおける振幅特性の時間平均を周波数帯域ごとに算出する。なお、既存品実測スペクトルは、周波数３００Ｈｚ以上の周波数帯域についても生成される。
具体的には、既存品にホワイトノイズを入力して発音させた音を収音し、該収音された音を表す音データをフーリエ変換し、パワースペクトルを経時的に生成する。そして、該パワースペクトルにおける振幅特性の時間平均を周波数帯域ごとに算出する。なお、既存品実測スペクトルは、周波数３００Ｈｚ以上の周波数帯域についても生成される。
図８に、ステップＳＡ１０およびステップＳＡ４０において生成された、既存品シミュレーションスペクトルおよび既存品実測スペクトルを示す。なお、既存品実測スペクトルは、周波数３００Ｈｚ以下の周波数帯域について取り出して示す。同図に示されるように、既存品シミュレーションスペクトルは、既存品実測スペクトルの定性的な傾向は捉えているものの、その絶対値には乖離が見られる。ここで見られる乖離は、従来のシミュレーションにおいて各種条件の設定時に生じる誤差や既存品実測スペクトルを求める際に生じる誤差などに由来する。
ステップＳＡ５０において、制御部１１は、ステップＳＡ４０で得られた既存品実測スペクトルに、ステップＳＡ３０において生成された比率データを積算することにより補正し、新たな周波数スペクトル（以下、特性予測スペクトル）を生成する。すなわち、既存品実測スペクトルの各周波数帯域の値に、比率データにおいて対応する周波数帯域の値を積算する。
図９は、そのようにして生成された特性予測スペクトルを示す。該特性予測スペクトルにおいては、図８に示した既存品実測スペクトルにおける、１００Ｈｚ〜２００Ｈｚと、２１０Ｈｚ〜２８０Ｈｚの周波数成分が強調されている。該強調されている周波数成分は、図７に示した比率データにおいて値が１を越える周波数帯域である。
図９は、そのようにして生成された特性予測スペクトルを示す。該特性予測スペクトルにおいては、図８に示した既存品実測スペクトルにおける、１００Ｈｚ〜２００Ｈｚと、２１０Ｈｚ〜２８０Ｈｚの周波数成分が強調されている。該強調されている周波数成分は、図７に示した比率データにおいて値が１を越える周波数帯域である。
ここで、特性予測スペクトルについて簡単にまとめる。図１０には、開発品シミュレーションスペクトルと、開発品実測スペクトルが示されている。開発品実測スペクトルとは、既存品実測スペクトルと同様の周波数スペクトルを開発品について生成したものである。本発明は、そもそも開発品が開発される前に開発品の音響特性を予測することが目的であるため、開発品は未だ実在しないとの状況を想定しているが、ここでは開発品の音響特性を予測した特性予測スペクトルの精度を検証するために、完成した開発品の開発品実測スペクトルを示す。
同図において、開発品シミュレーションスペクトルは定性的な傾向は捉えているものの、開発品実測スペクトルとの乖離が見られる。ここで見られる乖離は、既存品についてシミュレーションと実測との間に見られた乖離（図８参照）と同様に、従来のシミュレーションにおいて各種条件の設定時に生じる誤差や開発品実測スペクトルを求める際に生じる誤差などに由来する。
図１１には、特性予測スペクトルと、開発品実測スペクトルが示されている。同図からも明らかであるように、特性予測スペクトルは、従来のシミュレーションにより生成された開発品シミュレーションスペクトル（図１０参照）よりも、開発品実測スペクトルを精度良く予測できていることがわかる。
さて、図８においては、既存品実測スペクトルの周波数３００Ｈｚ以下の成分について示したが、図１２には、全ての周波数帯域について横軸を対数表示して示す。
ステップＳＡ６０において、制御部１１は、特性予測スペクトルを、図１２に示す既存品実測スペクトルと合成する。すなわち、特性予測スペクトルの周波数３００Ｈｚ以下の部分と、既存品実測スペクトルの周波数３００Ｈｚを上回る部分をつなぎ合わせ、新たな周波数スペクトル（以下、合成スペクトル）を生成する。
図１３は、生成された合成スペクトルを示した図である。図１４には、図１１にも示した開発品実測スペクトルを全ての周波数について示す。両図を比較すると明らかであるように、合成スペクトルは、開発品実測スペクトルの特に低周波数帯域について精度良く予測できていることがわかる。
ステップＳＡ６０において、制御部１１は、特性予測スペクトルを、図１２に示す既存品実測スペクトルと合成する。すなわち、特性予測スペクトルの周波数３００Ｈｚ以下の部分と、既存品実測スペクトルの周波数３００Ｈｚを上回る部分をつなぎ合わせ、新たな周波数スペクトル（以下、合成スペクトル）を生成する。
図１３は、生成された合成スペクトルを示した図である。図１４には、図１１にも示した開発品実測スペクトルを全ての周波数について示す。両図を比較すると明らかであるように、合成スペクトルは、開発品実測スペクトルの特に低周波数帯域について精度良く予測できていることがわかる。
ステップＳＡ７０において、制御部１１は、既存品が内蔵する音源により発音された音を表す音データをマイクロホン１９から受取る。ここでは、既存品であるキーボードにより任意の演奏がなされ、該演奏音が収音されれば良い。
ステップＳＡ８０において、制御部１１は、ステップＳＡ７０において受取った音データに対し、ステップＳＡ６０において生成した合成スペクトルを積算することにより、開発品から発音されると予測される音（予測音）を表す予測音データを生成し、音声処理部１７に出力する。音声処理部１７は、受取った予測音データをアナログ信号に変換し、スピーカ１８から放音させる。該放音された音は、ステップＳＡ７０における既存品による演奏が、あたかも開発品により発音されているように音響特性が変換されている。
（Ｂ−３；まとめ）
以上の処理をまとめる。従来のシミュレーションにおいては、各種条件の設定時に入力されるパラメータに理想値からの誤差が生じていた。該誤差は、既存品および開発品の複雑な構造などに由来するものであり、パラメータを制御することにより該誤差を小さくすることには限界があった。
本発明においては、既存品と開発品それぞれについて、従来のシミュレーションにより音響特性（周波数スペクトル）の差分（比）を算出する処理を行った。該処理は、シミュレーション結果同士で差分を生成するものであるため、上記シミュレーションの設定時に生じる誤差の多くが互いにキャンセルされる。一方、既存品と開発品との間の差異についてはキャンセルされない。その結果、既存品と開発品との差異を精度良く抽出したデータ（比率データ）が生成される。そして、比率データに基づいて既存品実測スペクトルが補正されることにより、開発品実測スペクトルを精度良く予測した特性予測スペクトルが生成される。最終的に、特性予測スペクトルを既存品による音データに対し積算することにより、既存品の音データを開発品の音を模した音データに変換することが可能となる。このようにして生成された音データは、実物が発する音と聴感的な差異が小さく、仮想の発音体の音を適切に評価することが可能である。
なお、本実施形態においては、周波数３００Ｈｚ以下の成分について特性予測スペクトルを生成して開発品の音の予測に用いたが、周波数がより高い成分についても特性予測スペクトルを生成しても良い。
また、シミュレーション手法は、差分法、境界要素法、有限要素法等、種々の手法から、計算精度、計算時間、計算容量等の計算条件に応じて適宜選択すればよい。
以上の処理をまとめる。従来のシミュレーションにおいては、各種条件の設定時に入力されるパラメータに理想値からの誤差が生じていた。該誤差は、既存品および開発品の複雑な構造などに由来するものであり、パラメータを制御することにより該誤差を小さくすることには限界があった。
本発明においては、既存品と開発品それぞれについて、従来のシミュレーションにより音響特性（周波数スペクトル）の差分（比）を算出する処理を行った。該処理は、シミュレーション結果同士で差分を生成するものであるため、上記シミュレーションの設定時に生じる誤差の多くが互いにキャンセルされる。一方、既存品と開発品との間の差異についてはキャンセルされない。その結果、既存品と開発品との差異を精度良く抽出したデータ（比率データ）が生成される。そして、比率データに基づいて既存品実測スペクトルが補正されることにより、開発品実測スペクトルを精度良く予測した特性予測スペクトルが生成される。最終的に、特性予測スペクトルを既存品による音データに対し積算することにより、既存品の音データを開発品の音を模した音データに変換することが可能となる。このようにして生成された音データは、実物が発する音と聴感的な差異が小さく、仮想の発音体の音を適切に評価することが可能である。
なお、本実施形態においては、周波数３００Ｈｚ以下の成分について特性予測スペクトルを生成して開発品の音の予測に用いたが、周波数がより高い成分についても特性予測スペクトルを生成しても良い。
また、シミュレーション手法は、差分法、境界要素法、有限要素法等、種々の手法から、計算精度、計算時間、計算容量等の計算条件に応じて適宜選択すればよい。
（Ｃ；変形例）
以上、本発明の実施形態について説明したが、本発明は以下のように種々の態様で実施することができる。なお、以下に説明する種々の実施形態を、適宜組み合わせて実施することも可能である。
以上、本発明の実施形態について説明したが、本発明は以下のように種々の態様で実施することができる。なお、以下に説明する種々の実施形態を、適宜組み合わせて実施することも可能である。
（１）上記実施形態においては、既存品シミュレーションスペクトルと開発品シミュレーションスペクトルとの差異に基づいて既存品実測スペクトルを補正し、開発品実測スペクトルを予測したデータである特性予測スペクトルを生成する場合について説明した。しかし、既存品シミュレーションスペクトル、開発品シミュレーションスペクトル、および既存品実測スペクトルを以下のようなアルゴリズムで処理することにより特性予測スペクトルを生成するようにしても良い。その場合の処理の概略を以下に説明する。なお、一例として、図１５に示したフローチャートに沿って処理を行う場合について説明する。また、処理の流れを概念的に示した図１６を参照して説明する。なお、図１６において、シミュレーション手段１１００、解析手段１１１０、差分算出手段１１２０、補正手段１１３０、合成手段１１４０、変形手段１１５０は、制御部１１により実現される手段である。
まず、シミュレーション手段１１００は、既存品Ｅの発音特性をシミュレートし、既存品シミュレーションスペクトルを生成する（ステップＳＢ１０）。また、既存品Ｅに実際に音を発音させ、該実際の発音を解析手段１１１０により解析し、既存品実測スペクトルを生成する（ステップＳＢ２０）。ステップＳＢ１０およびステップＳＢ２０は、それぞれ上記実施形態に説明したステップＳＡ１０およびステップＳＡ４０の処理内容と同じであるため、詳細な説明を省略する。
次に、ステップＳＢ３０において、差分算出手段１１２０は、既存品実測スペクトルと既存品シミュレーションスペクトルとを比較し、その差分を表す比率データを算出する。比率データは、周波数帯域ごとに既存品実測スペクトルの値を既存品シミュレーションスペクトルの値で除した比の値に基づいて生成する。
ステップＳＢ４０において、シミュレーション手段１１００は、開発品Ｄの発音特性をシミュレートし、開発品シミュレーションスペクトルを生成する。ステップＳＢ４０は、上記実施形態に説明したステップＳＡ２０の処理内容と同じであるため、詳細な説明を省略する。
ステップＳＢ５０において、補正手段１１３０は、ステップＳＢ４０で生成された開発品シミュレーションスペクトルを、ステップＳＢ３０で生成された比率データを用いて補正し、新たな周波数スペクトル（特性予測スペクトル）を生成する。その補正方法は、上記実施形態に説明したステップＳＡ５０と同様であるため、詳細な説明を省略する。
該生成された特性予測スペクトルには、開発品Ｄのシミュレーションによる発音特性に対して、既存品について見られたシミュレーションと実測との差異（シミュレーションによる誤差）が加味されていることから、開発品実測スペクトルを精度良く予測したスペクトルとなる。
ステップＳＢ６０において、合成手段１１４０は、特性予測スペクトルを既存品実測スペクトルの一部と合成し、新たな周波数スペクトル（合成スペクトル）を生成する。ステップＳＢ７０において既存品を実際に発音させて音データを得る。ステップＳＢ８０において、変形手段１１５０は、上記合成スペクトルに基づいて該音データを変形して予測音データを生成する。
次に、ステップＳＢ３０において、差分算出手段１１２０は、既存品実測スペクトルと既存品シミュレーションスペクトルとを比較し、その差分を表す比率データを算出する。比率データは、周波数帯域ごとに既存品実測スペクトルの値を既存品シミュレーションスペクトルの値で除した比の値に基づいて生成する。
ステップＳＢ４０において、シミュレーション手段１１００は、開発品Ｄの発音特性をシミュレートし、開発品シミュレーションスペクトルを生成する。ステップＳＢ４０は、上記実施形態に説明したステップＳＡ２０の処理内容と同じであるため、詳細な説明を省略する。
ステップＳＢ５０において、補正手段１１３０は、ステップＳＢ４０で生成された開発品シミュレーションスペクトルを、ステップＳＢ３０で生成された比率データを用いて補正し、新たな周波数スペクトル（特性予測スペクトル）を生成する。その補正方法は、上記実施形態に説明したステップＳＡ５０と同様であるため、詳細な説明を省略する。
該生成された特性予測スペクトルには、開発品Ｄのシミュレーションによる発音特性に対して、既存品について見られたシミュレーションと実測との差異（シミュレーションによる誤差）が加味されていることから、開発品実測スペクトルを精度良く予測したスペクトルとなる。
ステップＳＢ６０において、合成手段１１４０は、特性予測スペクトルを既存品実測スペクトルの一部と合成し、新たな周波数スペクトル（合成スペクトル）を生成する。ステップＳＢ７０において既存品を実際に発音させて音データを得る。ステップＳＢ８０において、変形手段１１５０は、上記合成スペクトルに基づいて該音データを変形して予測音データを生成する。
上記実施形態における処理方法と、本変形例（１）に示した処理方法とでは、互いに同じ既存品および開発品について処理を行った場合であっても、生成される特性予測スペクトルや、開発品の予測音データは異なったものとなる。そこで、上記実施形態における処理方法と、本変形例（１）に示した処理方法とをユーザが選択することができるようにしても良い。例えば、既存品や開発品の材質・形状に応じていずれかの処理方法を選択しても良いし、既存品シミュレーションスペクトルや開発品シミュレーションスペクトルの生成に用いた「従来のシミュレーション」における設定内容に応じていずれかの処理方法を選択しても良い。また、両方法により生成された予測音データを共に出力するようにしても良い。
（２）上記実施形態においては、比率データを生成する際、周波数帯域ごとに、既存品シミュレーションスペクトルの値を開発品シミュレーションスペクトルの値で除した比の値を計算する場合について説明した。しかし、上記比率データは、既存品シミュレーションスペクトルおよび開発品シミュレーションスペクトルについて、周波数帯域ごとの値の差を計算しても良い。また、上記各周波数における比の値または差の値に所定の係数を乗じたスペクトルを比率データとしても良い。
（３）上記実施形態においては、マイクロホン１９により収音された既存品による実際の音に基づいて既存品実測スペクトルを生成する場合について説明した。しかし、既存品実測スペクトルが予め得られている場合には、該スペクトルを記憶手段（ＲＯＭ１２、ＲＡＭ１３、ＨＤＤ１６など）に記憶しておき、適宜用いれば良い。
（４）上記実施形態においては、既存品にホワイトノイズを発音させ、該発音された音をマイクロホン１９で収音し解析することにより既存品実測スペクトルを生成する場合について説明した。しかし、ホワイトノイズに限らず、他の音を発音させることにより既存品実測スペクトルを生成しても良い。
（５）上記実施形態においては、本発明を開発中の電子楽器（例えば、キーボード）が発する音の予測に用いる場合について説明した。しかし、本発明は、電子楽器が発音する音の予測に限られない。例えば、電子楽器以外の発音体であるゲーム機、スピーカ装置、携帯電話機などに本発明を適用しても良いし、ゴルフの打撃音、カメラのシャッター音、自動車のドアミラーの風切り音など、発音を目的とされていないものが生じる音について評価するために用いても良い。
（６）上記実施形態においては、本発明に係る仮想音生成装置１に特徴的な機能を実現するための制御プログラムを、ＲＯＭ１２に予め書き込んでおく場合について説明したが、磁気テープ、磁気ディスク、フレキシブルディスク、光記録媒体、光磁気記録媒体、ＲＡＭ、ＲＯＭなどのコンピュータ読み取り可能な記録媒体に上記制御プログラムを記録して配布するとしても良く、インターネット網などの電気通信回線経由のダウンロードにより上記制御プログラムを配布するようにしても良い。
（７）上記実施形態においては、開発品が発する音の周波数応答についてシミュレーションを行う場合について説明した。しかし、以下のようにして、開発品が発する音の他の特性についてシミュレーションしても良い。
例えば、音の振幅の減衰特性についてシミュレーションしても良い。その場合の実施形態を以下に簡潔に説明する。仮想音生成装置１の構成は、ＲＯＭ１２に格納された制御プログラムを除いて、上記実施形態と同様である。
その処理内容は、以下の通りである。まず、既存品および開発品について、従来のシミュレーション方法により振幅の減衰特性をシミュレートする。例えば、既存品のキーボードの筐体の中には大きな空洞があるため音の減衰が比較的遅く、開発品のキーボードの筐体の中には小さい空洞があるため音の減衰が比較的早い、などのようにシミュレートされる。そして、既存品と開発品の間で、シミュレーションによる減衰特性波形の差分（比）を算出する。一方、既存品を実際に発音させ、既存品による発音の減衰特性を測定する。次に、上記減衰特性波形の差分（比）に基づいて、上記既存品による発音の減衰特性を補正することにより、新たな減衰特性波形を生成する。該生成された減衰特性波形には、既存品による発音の減衰特性に対して、シミュレーションにおいて見られた既存品と開発品の間での減衰特性波形の差分が加味されていることから、開発品による発音の減衰特性を精度良く予測したデータとなる。この後、既存品により発音された音の振幅を該生成された減衰特性波形に基づいて変形することで、減衰特性が開発品と類似した音データが生成される。以上に説明した振幅の減衰特性を制御する処理は、上記実施形態における周波数応答に関するシミュレーションと併せて行っても良い。
また、音の位相特性についてシミュレーションしても良い。シミュレーションおよび実測により得られる周波数応答からは、振幅（Amplitude）と位相（Phase）の情報を抽出することができる。すなわち、周波数応答を複素数として表記した場合、振幅は複素数の絶対値として算出され、位相は虚部と実部の比のアークタンジェントとして算出される。既存品および開発品に関するシミュレーション結果、および既存品に関する実測結果から、開発品が発する音の位相特性をシミュレーションし、開発品が発する音と位相特性が類似した音を生成しても良い。
例えば、音の振幅の減衰特性についてシミュレーションしても良い。その場合の実施形態を以下に簡潔に説明する。仮想音生成装置１の構成は、ＲＯＭ１２に格納された制御プログラムを除いて、上記実施形態と同様である。
その処理内容は、以下の通りである。まず、既存品および開発品について、従来のシミュレーション方法により振幅の減衰特性をシミュレートする。例えば、既存品のキーボードの筐体の中には大きな空洞があるため音の減衰が比較的遅く、開発品のキーボードの筐体の中には小さい空洞があるため音の減衰が比較的早い、などのようにシミュレートされる。そして、既存品と開発品の間で、シミュレーションによる減衰特性波形の差分（比）を算出する。一方、既存品を実際に発音させ、既存品による発音の減衰特性を測定する。次に、上記減衰特性波形の差分（比）に基づいて、上記既存品による発音の減衰特性を補正することにより、新たな減衰特性波形を生成する。該生成された減衰特性波形には、既存品による発音の減衰特性に対して、シミュレーションにおいて見られた既存品と開発品の間での減衰特性波形の差分が加味されていることから、開発品による発音の減衰特性を精度良く予測したデータとなる。この後、既存品により発音された音の振幅を該生成された減衰特性波形に基づいて変形することで、減衰特性が開発品と類似した音データが生成される。以上に説明した振幅の減衰特性を制御する処理は、上記実施形態における周波数応答に関するシミュレーションと併せて行っても良い。
また、音の位相特性についてシミュレーションしても良い。シミュレーションおよび実測により得られる周波数応答からは、振幅（Amplitude）と位相（Phase）の情報を抽出することができる。すなわち、周波数応答を複素数として表記した場合、振幅は複素数の絶対値として算出され、位相は虚部と実部の比のアークタンジェントとして算出される。既存品および開発品に関するシミュレーション結果、および既存品に関する実測結果から、開発品が発する音の位相特性をシミュレーションし、開発品が発する音と位相特性が類似した音を生成しても良い。
（８）上記実施形態においては、キーボードという複数のユニットから構成されるシステム全体における発音特性をシミュレーションする場合について説明した。しかし、一旦該発音体を構成する各ユニットについて本発明に係るシミュレーション方法を用いて発音特性をシミュレートし、該シミュレーション結果を用いてシステム全体における発音特性をシミュレートしても良い。以下では、具体例を挙げて説明する。
（例１）キーボードは、入力された音信号に基づいて音を自ら発音するユニットであるスピーカ（発音体）、および自らは発音せず、前記発音体に音響的に接続して発音体の音を伝搬させるユニット（構造体）であるスピーカキャビネットやキーボードキャビネットなど、複数のユニットから構成される。そこで、スピーカについては本発明に係るシミュレーション方法で発音特性をシミュレートし、前記スピーカ以外のユニットについては実測の結果または他のシミュレーション手段により得られる最も良いシミュレーション結果を採用し、それらの結果を組み合わせて評価することにより、キーボードというシステム（上記発音体と構造体とによる複合体）の発音特性をシミュレートしても良い。例えば、実測により、スピーカキャビネットおよびキーボードキャビネットが高周波数成分を低減する傾向があるなどといったことが判明した場合、本発明のシミュレーション方法をスピーカに適用して得られた発音特性（周波数応答）の高周波数部分に変更を加えたものを、仮想のスピーカが内蔵されたキーボード（発音体と構造体から構成される複合体）の発音特性とシミュレートしても良い。
（例２）また、上記組み合わせによるシミュレーション方法を以下のように応用しても良い。例えば、上記実施形態のようにして発音特性が明らかになった仮想のキーボード（発音体）を、既存または仮想の音響空間（構造体）に設置する場合について考える。このような場合、該音響空間が有する音の伝達特性を、実測または他のシミュレーション手段により得ておき、本発明に係るシミュレーション方法により得られたキーボードの発音特性と組み合わせることにより、上記キーボードが設置された音響空間（複合体）から漏れ聞こえるキーボードの漏洩音の特性を評価する、といったことも可能である。例えば、従来のシミュレーション方法により、仮想の音響空間が高周波数成分を低減する傾向があるなどといったことが判明した場合、本発明のシミュレーション方法により得られたキーボードの発音特性（周波数応答）の高周波数部分に変更を加えたものを、仮想のキーボードが設置された仮想の音響空間の発音特性（漏れ聞こえる音の特性）としても良い。
（例３）上記例１および２においては、本発明に係るシミュレーション方法と、実測または他のシミュレーション手段とを組み合わせてシステム全体の発音特性を評価する場合について説明した。しかし、本発明に係るシミュレーション方法を複数組み合わせて新たな発音特性をシミュレートしても良い。例えば、仮想のキーボード（発音体）を仮想の音響空間（構造体）に設置した場合について以下のようにシミュレートすることも可能である。
まず、仮想の音響空間による音の「伝達特性」を、以下に説明する方法Ａまたは方法Ｂに示すように、本発明を応用したシミュレーション方法によりシミュレートしておく。
（方法Ａ）従来のシミュレーション方法で、既存および仮想の音響空間について音の伝達特性（周波数応答）をシミュレートする。一方、既存の音響空間内にホワイトノイズの音源を置き、漏れ聞こえる音を実測するなどして、既存の音響空間の音の伝達特性を測定する。そして、該実測された既存の音響空間の音の伝達特性を、上記仮想と既存の音響空間に関するシミュレーション結果の比で補正することにより、仮想の音響空間の伝達特性をシミュレートする。
（方法Ｂ）従来のシミュレーション方法で、既存および仮想の音響空間について音の伝達特性（周波数応答）をシミュレートする。一方、既存の音響空間内にホワイトノイズの音源を置き、漏れ聞こえる音を実測するなどして、既存の音響空間の音の伝達特性を測定する。そして、仮想の音響空間に関するシミュレーション結果を、実測された既存の音響空間の音の伝達特性と既存の音響空間に関するシミュレーション結果の比で補正することにより、仮想の音響空間の伝達特性をシミュレートする。
また、上記実施形態に記載されたシミュレーション方法によりキーボードの発音特性についても別途シミュレートする。そして、該キーボードの発音特性と上記方法Ａまたは方法Ｂによりシミュレーションされた音響空間の伝達特性とを組み合わせることにより、仮想のキーボードが設置された仮想の音響空間（複合体）から発音される（漏れ聞こえる）音の特性をシミュレーションすることができる。
以上のようなシミュレーション方法によれば、システムを構成する１または複数のユニット（この場合、音響空間とキーボード）について本発明に係るシミュレーションを行えば、その結果に基づいてシステム全体（キーボードが備え付けられた音響空間）の発音特性がシミュレートされる。
なお、上記説明において、「構造体」として、キーボードにおけるスピーカキャビネットやキーボードキャビネット、発音体が設置される音響空間（部屋）などを例示した。これらの構造体は、構造体自体に空隙を有していたり、発音体と複合体を形成した際に該複合体内に空隙が生じたりする。従って、該複合体全体としての発音特性には、空気を介して音が伝搬されることによる効果が含まれていると考えられる。しかし、「構造体」はそのような空気伝搬を行う構造に限られるものではない。発音体から発せられた音を構造体自体の振動により固体伝搬音として伝搬する構造体でも良い。該構造体として、例えば、コンクリートの床や発音体を支持するシャフトなどの構造体が挙げられ、床の上に設置されたキーボード、シャフトで支持されたスピーカ、空間に設置されたゲーム機等のように、発音体の複雑な設置状況や支持状態を含めてシミュレートすることが困難である場合などにも適用することが可能である。
また、システムに含まれる特定のユニットのみを取り替えた場合のシステム全体の音響特性を評価するような場合に、該取り替えたユニットの発音特性をシミュレートしさえすれば、システム全体のシミュレーションを行うことなくシステム全体のシミュレートが可能である。また、特に精密なシミュレーションを要求されるユニット（例えば、スピーカ周辺）などに対して本発明に係るシミュレーション方法を適用し、他のユニットについては簡易なシミュレーション方法を用いる、といったことも可能である。
（例１）キーボードは、入力された音信号に基づいて音を自ら発音するユニットであるスピーカ（発音体）、および自らは発音せず、前記発音体に音響的に接続して発音体の音を伝搬させるユニット（構造体）であるスピーカキャビネットやキーボードキャビネットなど、複数のユニットから構成される。そこで、スピーカについては本発明に係るシミュレーション方法で発音特性をシミュレートし、前記スピーカ以外のユニットについては実測の結果または他のシミュレーション手段により得られる最も良いシミュレーション結果を採用し、それらの結果を組み合わせて評価することにより、キーボードというシステム（上記発音体と構造体とによる複合体）の発音特性をシミュレートしても良い。例えば、実測により、スピーカキャビネットおよびキーボードキャビネットが高周波数成分を低減する傾向があるなどといったことが判明した場合、本発明のシミュレーション方法をスピーカに適用して得られた発音特性（周波数応答）の高周波数部分に変更を加えたものを、仮想のスピーカが内蔵されたキーボード（発音体と構造体から構成される複合体）の発音特性とシミュレートしても良い。
（例２）また、上記組み合わせによるシミュレーション方法を以下のように応用しても良い。例えば、上記実施形態のようにして発音特性が明らかになった仮想のキーボード（発音体）を、既存または仮想の音響空間（構造体）に設置する場合について考える。このような場合、該音響空間が有する音の伝達特性を、実測または他のシミュレーション手段により得ておき、本発明に係るシミュレーション方法により得られたキーボードの発音特性と組み合わせることにより、上記キーボードが設置された音響空間（複合体）から漏れ聞こえるキーボードの漏洩音の特性を評価する、といったことも可能である。例えば、従来のシミュレーション方法により、仮想の音響空間が高周波数成分を低減する傾向があるなどといったことが判明した場合、本発明のシミュレーション方法により得られたキーボードの発音特性（周波数応答）の高周波数部分に変更を加えたものを、仮想のキーボードが設置された仮想の音響空間の発音特性（漏れ聞こえる音の特性）としても良い。
（例３）上記例１および２においては、本発明に係るシミュレーション方法と、実測または他のシミュレーション手段とを組み合わせてシステム全体の発音特性を評価する場合について説明した。しかし、本発明に係るシミュレーション方法を複数組み合わせて新たな発音特性をシミュレートしても良い。例えば、仮想のキーボード（発音体）を仮想の音響空間（構造体）に設置した場合について以下のようにシミュレートすることも可能である。
まず、仮想の音響空間による音の「伝達特性」を、以下に説明する方法Ａまたは方法Ｂに示すように、本発明を応用したシミュレーション方法によりシミュレートしておく。
（方法Ａ）従来のシミュレーション方法で、既存および仮想の音響空間について音の伝達特性（周波数応答）をシミュレートする。一方、既存の音響空間内にホワイトノイズの音源を置き、漏れ聞こえる音を実測するなどして、既存の音響空間の音の伝達特性を測定する。そして、該実測された既存の音響空間の音の伝達特性を、上記仮想と既存の音響空間に関するシミュレーション結果の比で補正することにより、仮想の音響空間の伝達特性をシミュレートする。
（方法Ｂ）従来のシミュレーション方法で、既存および仮想の音響空間について音の伝達特性（周波数応答）をシミュレートする。一方、既存の音響空間内にホワイトノイズの音源を置き、漏れ聞こえる音を実測するなどして、既存の音響空間の音の伝達特性を測定する。そして、仮想の音響空間に関するシミュレーション結果を、実測された既存の音響空間の音の伝達特性と既存の音響空間に関するシミュレーション結果の比で補正することにより、仮想の音響空間の伝達特性をシミュレートする。
また、上記実施形態に記載されたシミュレーション方法によりキーボードの発音特性についても別途シミュレートする。そして、該キーボードの発音特性と上記方法Ａまたは方法Ｂによりシミュレーションされた音響空間の伝達特性とを組み合わせることにより、仮想のキーボードが設置された仮想の音響空間（複合体）から発音される（漏れ聞こえる）音の特性をシミュレーションすることができる。
以上のようなシミュレーション方法によれば、システムを構成する１または複数のユニット（この場合、音響空間とキーボード）について本発明に係るシミュレーションを行えば、その結果に基づいてシステム全体（キーボードが備え付けられた音響空間）の発音特性がシミュレートされる。
なお、上記説明において、「構造体」として、キーボードにおけるスピーカキャビネットやキーボードキャビネット、発音体が設置される音響空間（部屋）などを例示した。これらの構造体は、構造体自体に空隙を有していたり、発音体と複合体を形成した際に該複合体内に空隙が生じたりする。従って、該複合体全体としての発音特性には、空気を介して音が伝搬されることによる効果が含まれていると考えられる。しかし、「構造体」はそのような空気伝搬を行う構造に限られるものではない。発音体から発せられた音を構造体自体の振動により固体伝搬音として伝搬する構造体でも良い。該構造体として、例えば、コンクリートの床や発音体を支持するシャフトなどの構造体が挙げられ、床の上に設置されたキーボード、シャフトで支持されたスピーカ、空間に設置されたゲーム機等のように、発音体の複雑な設置状況や支持状態を含めてシミュレートすることが困難である場合などにも適用することが可能である。
また、システムに含まれる特定のユニットのみを取り替えた場合のシステム全体の音響特性を評価するような場合に、該取り替えたユニットの発音特性をシミュレートしさえすれば、システム全体のシミュレーションを行うことなくシステム全体のシミュレートが可能である。また、特に精密なシミュレーションを要求されるユニット（例えば、スピーカ周辺）などに対して本発明に係るシミュレーション方法を適用し、他のユニットについては簡易なシミュレーション方法を用いる、といったことも可能である。
１…仮想音生成装置、１１…制御部、１２…ＲＯＭ、１３…ＲＡＭ、１４…操作部、１５…表示部、１６…ＨＤＤ、１７…音声処理部、１８…スピーカ、１９…マイクロホン、１１０、１１００…シミュレーション手段、１１１、１１１０…解析手段、１１２、１１２０…差分算出手段、１１３、１１３０…補正手段、１１４、１１４０…合成手段、１１５、１１５０…変形手段
Claims (13)
- 特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、
仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、
前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、
前記既存発音体シミュレーション手段と前記仮想発音体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第１の差異データ生成手段と、
前記差異データを用いて、前記既存発音体測定手段による測定結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第１の特性補正手段と
を具備することを特徴とするシミュレーション装置。 - 特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、
仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、
前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、
前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第２の差異データ生成手段と、
前記差異データを用いて、前記仮想発音体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第２の特性補正手段と
を具備することを特徴とするシミュレーション装置。 - 前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、周波数応答に関する差異を算出し、
前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定する
ことを特徴とする請求項１に記載のシミュレーション装置。 - 前記第１の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記仮想発音体シミュレーション手段によるシミュレーション結果から、位相特性に関する差異を算出し、
前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定する
ことを特徴とする請求項１または３に記載のシミュレーション装置。 - 前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、周波数応答に関する差異を算出し、
前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の周波数応答を測定することを特徴とする請求項２に記載のシミュレーション装置。 - 前記第２の差異データ生成手段は、前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とから、位相特性に関する差異を算出し、
前記既存発音体測定手段は、前記特定の発音体を実際に発音させた際の位相特性を測定することを特徴とする請求項２または５に記載のシミュレーション装置。 - 前記特定の発音体の演奏内容を表す音データを受取る受取手段と、
前記受取手段が受取った音データを、前記仮想発音体予測データに基づいて変形し出力する音データ変形手段と
を具備することを特徴とする請求項１ないし６のいずれかに記載のシミュレーション装置。 - 特定の構造体に係る音の伝達特性をシミュレートする伝達特性シミュレート手段と、
前記伝達特性シミュレート手段によるシミュレーション結果を用いて前記仮想発音体予測データを補正して、前記仮想の発音体が発音した音が前記特定の構造体を伝達した結果の複合体発音特性を予測する第１の予測データ補正手段と
を具備することを特徴とする請求項１ないし７のいずれかに記載のシミュレーション装置。 - 特定の構造体に係る音の伝達特性を測定する伝達特性測定手段と、
前記伝達特性測定手段による測定結果を用いて前記仮想発音体予測データを補正して、前記仮想の発音体が発音した音が前記特定の構造体を伝達した結果の複合体発音特性を予測する第２の予測データ補正手段と
を具備することを特徴とする請求項１ないし７のいずれかに記載のシミュレーション装置。 - 特定の構造体に係る音の伝達特性をシミュレートする既存構造体シミュレーション手段と、
仮想の構造体に係る音の伝達特性をシミュレートする仮想構造体シミュレーション手段と、
前記特定の構造体が音を伝達する際の伝達特性を測定する既存構造体測定手段と、
前記既存構造体シミュレーション手段と前記仮想構造体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第３の差異データ生成手段と、
前記差異データを用いて、前記既存構造体測定手段による測定結果を補正し、前記仮想の構造体による音の伝達特性である仮想構造体予測データを生成する第３の特性補正手段と、
前記仮想発音体予測データを、前記第３の特性補正手段が生成した仮想構造体予測データに基づいて補正して、前記仮想の発音体が発音した音が前記仮想の構造体を伝達した結果の複合体発音特性を予測する第３の予測データ補正手段と
を具備することを特徴とする請求項１ないし７のいずれかに記載のシミュレーション装置。 - 特定の構造体に係る音の伝達特性をシミュレートする既存構造体シミュレーション手段と、
仮想の構造体に係る音の伝達特性をシミュレートする仮想構造体シミュレーション手段と、
前記特定の構造体が音を伝達する際の伝達特性を測定する既存構造体測定手段と、
前記既存構造体シミュレーション手段によるシミュレーション結果と前記既存構造体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第４の差異データ生成手段と、
前記差異データを用いて、前記仮想構造体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の構造体による音の伝達特性である仮想構造体予測データを生成する第４の特性補正手段と、
前記仮想発音体予測データを、前記第４の特性補正手段が生成した仮想構造体予測データに基づいて補正して、前記仮想の発音体が発音した音が前記仮想の構造体を伝達した結果の複合体発音特性を予測する第４の予測データ補正手段と
を具備することを特徴とする請求項１ないし７のいずれかに記載のシミュレーション装置。 - コンピュータを、
特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、
仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、
前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、
前記既存発音体シミュレーション手段と前記仮想発音体シミュレーション手段とによるそれぞれのシミュレーション結果を比較し、その差異を表す差異データを生成する第１の差異データ生成手段と、
前記差異データを用いて、前記既存発音体測定手段による測定結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第１の特性補正手段
として機能させるためのプログラム。 - コンピュータを、
特定の発音体の発音特性をシミュレートする既存発音体シミュレーション手段と、
仮想の発音体の発音特性をシミュレートする仮想発音体シミュレーション手段と、
前記特定の発音体を実際に発音させた際の発音特性を測定する既存発音体測定手段と、
前記既存発音体シミュレーション手段によるシミュレーション結果と前記既存発音体測定手段による測定結果とを比較し、その差異を表す差異データを生成する第２の差異データ生成手段と、
前記差異データを用いて、前記仮想発音体シミュレーション手段によるシミュレーション結果を補正し、前記仮想の発音体の発音特性である仮想発音体予測データを生成する第２の特性補正手段
として機能させるためのプログラム。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007290790A JP5050789B2 (ja) | 2007-11-08 | 2007-11-08 | シミュレーション装置およびプログラム |
| US12/266,496 US8321043B2 (en) | 2007-11-08 | 2008-11-06 | Simulation apparatus and program |
| EP08019523.3A EP2058795B1 (en) | 2007-11-08 | 2008-11-07 | Simulation apparatus and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007290790A JP5050789B2 (ja) | 2007-11-08 | 2007-11-08 | シミュレーション装置およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2009116169A JP2009116169A (ja) | 2009-05-28 |
| JP5050789B2 true JP5050789B2 (ja) | 2012-10-17 |
Family
ID=40380384
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007290790A Expired - Fee Related JP5050789B2 (ja) | 2007-11-08 | 2007-11-08 | シミュレーション装置およびプログラム |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US8321043B2 (ja) |
| EP (1) | EP2058795B1 (ja) |
| JP (1) | JP5050789B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5459709B2 (ja) * | 2010-03-25 | 2014-04-02 | 一般財団法人電力中央研究所 | 電力機器の設計支援方法及び設計支援装置並びに設計支援プログラム |
| JP5987365B2 (ja) * | 2012-03-07 | 2016-09-07 | ヤマハ株式会社 | 伝達関数演算装置、およびプログラム |
| JP5949139B2 (ja) * | 2012-05-18 | 2016-07-06 | ヤマハ株式会社 | シミュレーション装置、およびプログラム |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2681349B2 (ja) | 1986-08-08 | 1997-11-26 | ヤマハ株式会社 | スピーカ再生装置 |
| US5572443A (en) | 1993-05-11 | 1996-11-05 | Yamaha Corporation | Acoustic characteristic correction device |
| US6760451B1 (en) | 1993-08-03 | 2004-07-06 | Peter Graham Craven | Compensating filters |
| US6350943B1 (en) * | 2000-12-28 | 2002-02-26 | Korg, Inc. | Electric instrument amplifier |
| JP2003090758A (ja) | 2001-09-17 | 2003-03-28 | Toray Ind Inc | 成形品の音響特性のシミュレーション方法および装置 |
| JP4674104B2 (ja) | 2004-03-24 | 2011-04-20 | 美津濃株式会社 | サウンドシミュレーションシステム |
| US7630501B2 (en) * | 2004-05-14 | 2009-12-08 | Microsoft Corporation | System and method for calibration of an acoustic system |
| JP2006065466A (ja) | 2004-08-25 | 2006-03-09 | Toyota Motor Corp | Ｆｅｍモデル化方法 |
| US8842847B2 (en) | 2005-01-06 | 2014-09-23 | Harman International Industries, Incorporated | System for simulating sound engineering effects |
| US20100010786A1 (en) * | 2008-07-08 | 2010-01-14 | Cynthia Maxwell | Sound synthesis method and software system for shape-changing geometric models |
-
2007
- 2007-11-08 JP JP2007290790A patent/JP5050789B2/ja not_active Expired - Fee Related
-
2008
- 2008-11-06 US US12/266,496 patent/US8321043B2/en active Active
- 2008-11-07 EP EP08019523.3A patent/EP2058795B1/en not_active Not-in-force
Also Published As
| Publication number | Publication date |
|---|---|
| EP2058795A3 (en) | 2010-12-01 |
| EP2058795A2 (en) | 2009-05-13 |
| US8321043B2 (en) | 2012-11-27 |
| JP2009116169A (ja) | 2009-05-28 |
| EP2058795B1 (en) | 2015-08-12 |
| US20090125135A1 (en) | 2009-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2602787B1 (en) | Signal processing device | |
| JP4716422B2 (ja) | 共鳴音発生装置 | |
| US8447585B2 (en) | System and method for characterizing, synthesizing, and/or canceling out acoustic signals from inanimate sound sources | |
| US8530736B2 (en) | Musical tone signal synthesis method, program and musical tone signal synthesis apparatus | |
| Woodhouse et al. | Interpreting the input admittance of violins and guitars | |
| JP4915773B2 (ja) | 伝達特性測定方法および装置 | |
| US10170094B2 (en) | Musical sound generating device, control method for same, storage medium, and electronic musical instrument | |
| JP2007193129A (ja) | 共鳴音音像生成装置および記憶媒体 | |
| Chatziioannou et al. | Investigating clarinet articulation using a physical model and an artificial blowing machine | |
| US9626949B2 (en) | System of modeling characteristics of a musical instrument | |
| JP4785053B2 (ja) | 共鳴音発生装置 | |
| US20170024495A1 (en) | Method of modeling characteristics of a musical instrument | |
| JP5050789B2 (ja) | シミュレーション装置およびプログラム | |
| WO2006013683A1 (ja) | 電子楽器 | |
| Lee et al. | Mathematical modelling and acoustical analysis of classical guitars and their soundboards | |
| JP5217875B2 (ja) | 音場支援装置、音場支援方法およびプログラム | |
| JP5332296B2 (ja) | 楽音合成装置およびプログラム | |
| Valiente et al. | Modeling the interaction between piano strings and the soundboard | |
| Miranda Valiente et al. | Influence of soundboard modelling approaches on piano string vibration | |
| JP3304443B2 (ja) | 楽音合成装置 | |
| JP5949139B2 (ja) | シミュレーション装置、およびプログラム | |
| JP2023143837A (ja) | 楽音信号合成方法、楽音信号合成装置およびプログラム | |
| JP2012252147A (ja) | 信号補正装置、音響再現システムおよびプログラム | |
| JP5987365B2 (ja) | 伝達関数演算装置、およびプログラム | |
| Bruyns | Sound Synthesis from Shape-Changing Geometric Models |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100921 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120601 |
|
| 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: 20120626 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120709 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5050789 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150803 Year of fee payment: 3 |
|
| LAPS | Cancellation because of no payment of annual fees |