JP5366875B2

JP5366875B2 - 音生成装置

Info

Publication number: JP5366875B2
Application number: JP2010093854A
Authority: JP
Inventors: 信治瀬戸; 賢二平工; 洋祐田部; 謙一郎松原; 靖高野
Original assignee: Hitachi Ltd
Current assignee: Hitachi Ltd
Priority date: 2010-04-15
Filing date: 2010-04-15
Publication date: 2013-12-11
Anticipated expiration: 2030-04-15
Also published as: JP2011226790A; US20110255706A1; US8885843B2

Description

本発明は、車両に搭載される部品が動作することで発生する車両音を生成する音生成装置に関する。

車両に搭載される部品が動作することで発生する車両音の生成方法としては、従来、部品がエンジンである場合について、たとえば、単発音の波形データを、吸気音、エンジン爆発音、排気音などといったその単発音を成す要素毎の波形データとして準備しておいて、エンジンの設計条件に応じて選択したそれらの要素毎の波形データに周波数変更や時間、振幅の揺らぎを付加した後に合成して出力することによって、エンジン音のシミュレーションを行う方法が示されている。

特開２００７−２５８９８０号公報

しかしながら、上述した従来の方法では、エンジン以外の部品が起因となって発生する音や複数の部品が起因となって発生する音を生成することはできない。したがって、エンジン以外の部品を新たな部品に交換したときに、その新たな部品が起因となって発生する音等については、実際に製作して車両に取り付けるまで、聞くことはできなかった。

本発明は、上記の点に鑑みてなされたものであり、その目的は、車両に搭載された複数の部品が動作することで発生する複数の音成分のうち、特定の音成分を変更した場合の車両音を生成することができ、これにより、車両音の善し悪しを迅速に判断することができる音生成装置を提供することにある。

上記課題を解決する本発明の音生成装置は、車両から発生する車両音の音圧データを取得し、現行部品が起因となって発生する現行部品起因音を演算し、車両音の音圧データから現行部品起因音を除去した現行部品除去音を演算し、新規部品が起因となって発生する新規部品起因音を演算し、現行部品除去音に新規部品起因音を合成して合成車両音の音圧波形を生成し、その合成車両音の音圧波形から合成車両音を発生させることを特徴としている。

本発明によれば、車両に搭載される特定の現行部品を新規部品に交換したときに車両から発生する車両音を合成した合成車両音を生成することができるので、例えば、車両に搭載された複数の部品が動作することで発生する複数の音成分のうち、特定の音成分を変更した場合の車両音を生成することができ、音の善し悪しを迅速に判断することができる。

本発明の第一の実施形態に係る音生成装置を示す図。本発明の第一の実施形態に係る実測音のフーリエ変換結果の一例を示す図。本発明の第一の実施形態に係る音圧の回転次数成分を示す図。本発明の第一の実施形態に係る音圧の回転次数成分の値の一例を示す図。本発明の第一の実施形態に係る複数条件の音圧の回転次数成分の値の一例を示す図。本発明の第一の実施形態に係るエンジン音算出方法を示す図。本発明の第一の実施形態に係る部品起因音の伝達経路を示す図。本発明の第二の実施形態に係る音生成装置を搭載した車両の概要図。

本発明の音生成装置の第一実施例を図１から図７を用いて説明する。
音生成装置１は、複数の部品が動作することで発生する複数の音成分のうち、特定の音成分を変更した場合の車両音を生成するものである。本実施例では、自動車の車内音のうちエンジンによって駆動される部品のひとつによって発生する音成分を変更した場合の車内音を生成する音生成装置１を例に説明する。

図１において、音生成装置１は、現行音取得手段としての現行音取得部２、現行部品起因音除去手段としての現行部品起因音除去部３、現行部品起因音演算手段としての現行部品起因音演算部４、新規部品起因音演算手段としての新規部品起因音演算部５、合成音生成手段としての合成音生成部６、合成音発生手段としての合成音発生部７とから構成される。

まず、現行音取得部２について説明する。
現行音取得部２では、自動車車内音の音圧データを取得する処理が行われる。この際、いろいろな条件での音圧データがあることが望ましい。例えば、車内音は主にエンジン回転数により変わるため、回転数を徐々に変化させた場合の音圧データを取得しておく。さらに必要に応じてアクセル開度や走行速度などを変化させた場合の音圧データを取得しておくことが望ましい。また、エンジン燃焼の一サイクルは燃料の吸気、圧縮、爆発、排気で行われ、このタイミングと車内音圧には関係があるためこの１サイクルを３６０°としたときの位相（以下、燃焼サイクル）と音圧データは同時に計測しておくことが望ましい。

つぎに、現行部品起因音除去部３について説明する。
自動車には、プーリーやチェーンなどによりエンジン駆動力が伝えられて駆動する回転部品がある。この回転部品の加振力がさまざまな経路を伝播して車内音となる。

現行部品起因音除去部３は、現行音取得部２で取得した音圧データから、回転部品のうちの一つから発生する音（以下、現行部品起因音と称す）を除去した音圧データを作成する部分である。以下、現行部品起因音を除去したデータは主にエンジンが起因となって発生する音であるため、エンジン音と称す。

現行部品起因音除去部３での処理について説明する。まず、現行音取得部２で得られた音圧データを回転数ごとに、フーリエ変換を行う。変換を行った結果の一例を図２に示す。図２は横軸が周波数、縦軸が回転数、色の濃さが音圧を表わす。音圧が高い部分を結ぶと直線となり、回転数に比例して周波数が変わる部分があることがわかる。これは、主にエンジンに起因する音であり、燃焼サイクルの１次成分とその倍数成分（次数成分）である。

この次数成分のうち、燃焼サイクルｎ次成分の車内音Ｐｎに注目すると、図３に示すように、横軸が回転数、縦軸が振幅Ａｎと位相遅れθｎのグラフが得られる。位相遅れは、燃焼サイクルの所定位置（たとえばピストン上死点での燃料吸入開始時）を０゜としたときの遅れを示す（以下、この位相遅れを燃焼サイクル基準位相遅れと称す）。回転数と振幅Ａｎ、燃焼サイクル基準位相遅れθｎの関係をデータとして、例えば、図４に示す形で保存しておく。また、さらにアクセル開度など別の関数となる場合には、例えば、図５に示すように階層的にデータを保存する。

一方、現行部品起因音は、現行部品回転軸の回転により発生するために、燃焼サイクルの次数成分ではなく、現行部品を駆動する部品軸の回転数の次数成分で発生する。

そこで、まずクランクシャフトと現行部品回転軸の回転比から現行部品起因音の燃焼サイクル次数を算出する。たとえば、クランクシャフト回転と現行部品の回転比が４:３で、現行部品回転数の１次、２次、３次、４次・・・で音が発生する部品について考える。クランクシャフト２回転で燃焼サイクル１回が行われる４サイクルエンジンを想定とすると、現行部品起因音が発生する現行部品回転数の1次、２次、３次、４次は、燃焼サイクル次数では、８／３次、１６／３次、８次、３２／３次、・・・となる。一方、エンジン音は燃焼サイクル次数で１次、２次、３次、・・・で発生する。

そのため、８次、１６次、・・・では現行部品起因音とエンジン音の周波数が一致する。そこで、エンジン音は車内音から現行部品起因音を差し引くことで求める。後述する現行部品起因音演算部４により、現行部品起因音Ｐｃの次数成分Ｐｃｎ（振幅Ａｃｎ、燃焼サイクル基準位相遅れθｃｎ）を回転数ごとに計算しておく。

エンジン音Ｐｅｎは車内音Ｐｎと現行部品起因音Ｐｃｎから
［数１］
Ｐｅｎ＝Ｐｎ−Ｐｃｎ・・・（１）
により求められる。

図６に示すようにＰｎ、Ｐｃｎは、それぞれ振幅（Ａｎ、Ａｃｎ）と燃焼サイクル基準位相遅れ（θｎ、θｃｎ）の情報を持っているため、エンジン音Ｐｅｎも振幅Ａｅｎと燃焼サイクル基準位相遅れθｅｎが求まる。これも図４に示す形で保存しておく。

つぎに、現行部品起因音演算部４について説明する。
一般に部品が起因となって発生する部品起因音は、部品からの振動などが複数の伝達経路を伝播することで発生する。図７に示すように、部品１５が３か所で車体１１に取り付けられる場合は、それぞれの取り付け点１７ａ、１７ｂ、１７ｃでの力が伝達経路１０ａ、１０ｂ、１０ｃを伝播し車内音となり、取り付け点１７の力Ｆ１、Ｆ２、Ｆ３と伝達関数Ｇ１、Ｇ２、Ｇ３をかけあわせることで車内音が求まる。一般化すると、すべてのｉに対する力Ｆｉと伝達関数Ｇｉを求めることで車内音が求まる。

そこで、現行部品起因音を求めるには、ｉ番目の取り付け点の力Ｆｉは回転数を徐々に変化させた時の時刻歴データを計測しておく。また同時に燃焼サイクルも測定しておく。このデータを用いて、ｉ番目の取り付け点の力Ｆｉを回転数ごとにフーリエ変換を行い次数成分の抽出を行うとＦｉの燃焼サイクルｎ次成分Ｆｎｉが求まる。Ｆｎｉは振幅と燃焼サイクル基準位相遅れの情報を持つ。したがって、（２）式に示す式により現行部品起因音Ｐｃのｎ次成分Ｐｃｎが求まる。

Ｆｎｉは振幅と燃焼サイクル基準位相遅れ、Ｇｉは入力に対する出力のゲインと入力に対する出力の位相遅れの情報をもつため、Ｐｃｎは振幅Ａｃｎと燃焼サイクル基準位相遅れθｃｎの情報を持つ。

ここではＦｉ、Ｇｉはあらかじめ計測しておいたが、たとえば有限要素法などにより求めてもよい。

つぎに新規部品起因音演算部５について示す。
新規部品起因音とは、車両に搭載されている現行部品を現行部品とは異なる仕様の新規部品に変更した場合や、部品から車内音に至る伝達経路を変更した場合などに、部品が起因で発生する音である。

新規部品起因音は、新規部品を搭載した場合には新規部品の車体へのｉ番目の取り付け点の力Ｆｉｎｅｗを算出し、車体取り付け点から車内音に至る伝達経路を変更した場合には車体取り付け点の力から車内音までの伝達関数Ｇｉｎｅｗを、それぞれ算出することで求まる。

Ｆｉｎｅｗは、新規部品回転軸の回転により発生するため、新規部品回転軸の回転数のｍ次成分Ｆｍｉｎｅｗを振幅と位相遅れの情報として算出しておく。ここでの位相遅れは新規部品回転軸の所定の角度を０°とし、それに対する遅れを示す（以下、この位相遅れを部品回転軸基準位相遅れと称す）。

これにより、新規部品起因音Ｐｃｎｅｗの新規部品回転数のｍ次成分Ｐｃｍｎｅｗが（３）式で表わされる。

Ｐｃｍｎｅｗも振幅Ａｃｍｎｎｅｗと部品回転軸基準位相遅れθｃｍｎｎｅｗとして保存しておく。

ここで、部品取り付け点で発生する力Ｆｍｉｎｅｗの算出方法について述べる。Ｆｍｉｎｅｗは、たとえば、計算機上で求めるものであり、有限要素法などを用いて部品内部で発生する力を入力、取り付け点で発生する力を出力して求める。また、実際に新規部品から発生する力を計測してもよい。

伝達関数Ｇｉｎｅｗは、計測して求めてもよいし、計算機上で有限要素法などを用いて求める方法でもよい。

つぎに、合成音生成部６について説明する。
合成音生成部６は、エンジン音と新規部品起因音を合成した合成車内音を生成するものである。まず再現したいエンジン回転数の時刻歴ωｃｓ（ｔ）を決め、これを積分することで（４）式に示す燃焼サイクルの時刻歴波形を生成する。

これによりエンジン音で発生する燃焼サイクルｎ次成分の音圧の時刻歴波形Ｐｅｎ（ｔ）は（５）式となる。

［数５］
Ｐｅｎ（ｔ）＝Ａｅｎ（ｔ）×Ｓｉｎ (ｎθｅ（ｔ）-θｅｎ（ｔ）) ・・・（５）

ここでＡｅｎ（ｔ）、θｅｎ（ｔ）は（１）式で算出した振幅と燃焼サイクル位相遅れのデータを元に各回転数間でのデータを補間するように作成する。

つぎに、各次数の音を足し合わせてエンジン音の音圧の時刻歴波形Ｐｅ（ｔ）は（６）式となる。

一方、新規部品回転軸の燃焼サイクルに対する回転比をｋとすると新規部品回転軸の角度の時刻歴波形は（７）式で求まる。

［数７］
θｃ（ｔ）=ｋθｅ（ｔ）＋θ０・・・（７）
θ０は燃焼サイクルの基準位置に対する新規部品回転軸の初期角度を表わす。

なお、回転比ｋや初期位相θ０は必要に応じて変更できるようにすることで、回転比や初期角度を変更したときの音の合成が容易にできる。この式から新規部品回転数のｍ次で発生する音圧の時刻歴波形Ｐｃｍｎｅｗ（ｔ）は（８）式で生成される。

［数８］
Ｐｃｍｎｅｗ（ｔ）＝Ａｃｍｎｅｗ（ｔ）×Ｓｉｎ（ｍθｃ（ｔ）−θｃｍｎｅｗ（ｔ））・・・（８）

ここでＡｃｍｎｅｗ（ｔ）が振幅、θｃｍｎｅｗ（ｔ）の値は（３）式で得られた値を元に各回転数間でデータを補間するように作成する。

つぎに、（９）式に示すように各次数の音圧の時刻歴波形を足し合わせることで新規部品起因音の音圧の時刻歴波形Ｐｃ（ｔ）が求められる。

車内音の音圧の時刻歴波形Ｐ（ｔ）はエンジン音の音圧の時刻歴波形Ｐｅ（ｔ）と部品起因音の音圧の時刻歴波形のＰｃ（ｔ）の合成なので（１０）式のように足し合わせて車内音が生成される。

［数１０］
Ｐ（ｔ）＝Ｐｅ（ｔ）＋Ｐｃ（ｔ）・・・（１０）

なお、回転数に依存しない成分（例えば風切り音、ロードノイズ）などは走行条件違いなどでエンジン停止時の車内音を計測しておき、それを（１０）式に足し合わせることで、走行時の音に関して、より現実に近い音を合成することが可能である。

つぎに、合成音発生部７について説明する。
合成音発生部７は、得られた音圧の時刻歴波形Ｐ（ｔ）を電気信号に変換し、スピーカ等で音を発生させるものである。

以上のような構成により、部品変更時や、部品から車内音までの伝達経路変更時の車内音をスピーカで発生させることができる。

本実施例によれば、車両に搭載された複数の部品が動作することで発生する複数の音成分のうち、特定の音成分を変更した場合の車両音を生成することができるため、音の善し悪しを迅速に判断することができる。

したがって、車両に搭載されている現行部品を現行部品とは仕様が異なる新規部品に交換したときの車内音を生成して発生させることができ、例えば、新規部品を製作する部品メーカは、実際に部品を交換することなく、顧客である車両メーカに対して、新規部品を取り付けた状態の車内音を聞かせることができ、事前に音の評価を受けることが可能となる。

なお、本実施例では車内音の生成、発生について記載したが、車外音など別の場所での音についても同様の方法で生成が可能である。また、エンジン以外の部品から発生する音としても同様の方法で生成が可能である。

つぎに、本発明の第二実施例について図８を用いて説明する。
図８は音生成装置１とその一部を搭載した車両２０を示す図である。

音生成装置１は、現行部品起因音演算部４と、新規部品起因音演算部５と、合成音生成部６と、合成音発生部７と、部品回転軸センサ２２とから構成され、そのうち少なくとも部品回転軸センサ２２、合成音生成部６、合成音発生装置７とを車両２０に搭載する。

部品回転軸センサ２２は、車両に搭載されていてエンジン回転駆動力により回転駆動される現行部品２１の回転軸の角度と回転速度を検出するものであり、角度か回転速度の一方から他方を演算できるようなものであればよい。

現行部品起因音演算部４では、あらかじめ部品軸の回転数の次数成分ごとに受音点２３での音圧の振幅Ａｃｍと部品回転軸基準位相遅れθｃｍを計算しておく。これは第一実施例で示した方法と同様に求める。

さらに、合成音発生部７から受音点２３までの伝達関数Ｇｒを考慮し、（１１）式に示すように

を演算し、振幅Ａｃｍと部品回転軸基準位相遅れθｃｍを回転数ごとに求めておくことが望ましい。

新規部品起因音演算部５では、受音点２３での新規部品回転数のｍ次成分の振幅Ａｃｍｎｅｗと部品回転軸基準位相遅れθｃｍｎｅｗを計算しておく。これは第一実施例で示した方法と同様の方法で求める。さらに、合成音発生部７から受音点２３までの伝達関数Ｇｒを考慮し、（１２）式に示すように

を演算しておくことが望ましい。振幅Ａｃｍｎｅｗと部品回転角基準θｃｍｎｅｗを回転数ごとに求めておく。

合成音生成部６では、現行部品音演算部で得られた結果と、新規部品音演算部で得られた結果、さらに部品回転軸センサ２２から得られる結果から部品音の時刻歴波形を（１３）式により生成する。

［数１３］
Ｐｃｍｏ（ｔ）＝Ａｃｍｏ（ｔ）×Ｓｉｎ（ｍθｃ（ｔ）−θｃｍｏ（ｔ））
・・・（１３）

ここで、θｃは部品回転軸センサ２２から得られる部品軸の回転角度であり、振幅Ａｃｍｏ（ｔ）と部品回転軸基準位相遅れθｃｍｏ（ｔ）は同様に部品回転軸センサ２２から得られる回転軸の回転速度の値から補間して求める。また、新規部品起因音の音圧の時刻歴波形を（１４）式により生成する。

［数１４］
Ｐｃｍｏｎｅｗ（ｔ）＝Ａｃｍｏｎｅｗ（ｔ）×Ｓｉｎ（ｍθｃ（ｔ）−θｃｍｏｎｅｗ（ｔ））・・・（１４）

ここで、θｃは部品回転軸センサ２２から得られる部品軸の回転角度であり、振幅Ａｃｍｏｎｅｗ（ｔ）と部品回転軸基準位相遅れθｃｍｏｎｅｗ（ｔ）は同様に部品回転軸センサ２２から得られる回転軸の回転速度の値から補間して求める。

そして（１５）式により合成音発生部７から発生させる音圧の時刻歴波形Ｐｏ（ｔ）を合成する。

この音圧の時刻歴波形Ｐｏ（ｔ）を音発生部７により発生させる。（１５）式の第１項で発生する音はもともとの車両で発生する現行部品起因音と受音点では逆位相の音となるため、現行部品起因音を打ち消すこととなる。第２項が新規部品音として受音点２３で聴こえる。また、部品起因音以外の音はそのまま受音点で聴こえることになる。

本実施例によれば、ある車両に取り付けられた部品を取り換えた場合や伝達経路や伝達関数が変更になった場合に聴こえる車内音を、実際に部品を交換することなく、あるいは伝達経路を変更することなく、聞くことが可能になり、かつ変更した部品以外の音は車両の実際の音をそのままに聴くことができるため、より現実に近い音が聴くことができる。

以上、本発明によれば、車両に搭載された複数の部品が動作することで発生する複数の音成分のうち、特定の音成分を変更した場合の車両音を生成することができるため、音の善し悪しを迅速に判断することができる。

１音生成装置
２現行音取得部
３現行部品起因音除去部
４現行部品起因音演算部
５新規部品起因音演算部
６合成音生成部
７合成音発生部
１０伝達経路
２０車両
２１回転駆動部品
２２部品回転軸センサ
２３受音点

Claims

車両に搭載される特定の現行部品を新規部品に交換したときに車両から発生する車両音を合成した合成車両音を生成する音生成装置であって、
車両から発生する車両音の音圧データを取得する現行音取得手段と、
前記現行部品が起因となって発生する現行部品起因音を演算する現行部品起因音演算手段と、
該車両音の音圧データから前記現行部品起因音を除去した現行部品除去音を演算する現行部品起因音除去手段と、
前記新規部品が起因となって発生する新規部品起因音を演算する新規部品起因音演算手段と、
前記現行部品除去音に前記新規部品起因音を合成して前記合成車両音の音圧波形を生成する合成音生成手段と、
該合成音生成手段により生成された合成車両音の音圧波形から合成車両音を発生させる合成音発生手段と、
を有することを特徴とする音生成装置。
前記現行音取得手段は、回転機械の複数の回転数で、車両で発生する車両音の音圧データを取得し、
前記現行部品起因音演算手段は、前記車両音のうち前記回転機械の駆動力によって駆動される特定の部品が起因で発生する現行部品起因音の音圧データを前記回転機械の回転数に比例した周波数である回転次数成分ごとの音圧次数成分データで演算し、
前記現行部品起因音除去手段は、前記車両音の音圧データから前記回転次数成分ごとの音圧次数成分データに変換し、現行部品起因音の音圧次数成分データを、前記車両音の音圧次数成分データから同じ位相基準で差し引くことで、前記回転次数成分の現行部品除去音の音圧次数成分データを求めることを特徴とする請求項１に記載の音生成装置。
前記音圧次数成分データが音圧の振幅と前記回転機械の所定の位置を基準とした位相遅れから構成されていることを特徴とする請求項２に記載の音生成装置。
新規部品起因音演算手段は、前記回転機械の回転駆動力により部品回転軸が回転されることで部品が起因で発生する新規部品起因音の音圧データを前記部品回転軸の回転数に比例する周波数である部品回転次数成分ごとの音圧次数成分データで演算し、
前記合成音生成手段は、前記回転機械の回転数を入力として、前記部品除去音の前記音圧次数成分データから前記現行部品除去音の音圧波形を生成し、前記部品回転軸の回転数を前記回転機械の回転数から生成し、前記新規部品音の音圧次数成分データから前記新規部品音の音圧波形を生成し、前記現行部品除去音の音圧波形と前記新規部品音の音圧波形を足し合わせて合成音の音圧波形を生成することを特徴とする請求項３に記載の音生成装置。
前記新規部品起因音の音圧次数成分データが振幅と前記部品回転軸の所定の位置を基準とした位相遅れから構成され、
前記合成音生成手段は、前記新規部品音の音圧次数成分データから前記新規部品音の音圧波形を前記部品回転軸の初期位置を任意に変更した音圧波形を生成し、前記現行部品除去音の音圧と前記新規部品音の音圧波形を足し合わせて合成音の音圧波形を生成することを特徴とする請求項４に記載の音生成装置。
前記現行音取得手段は、回転機械駆動時の車両音と回転機械停止時の車両走行時の車内音を取得し、
前記合成音生成手段は、前記現行部品除去音と新規部品起因音と回転機械停止時の車内音を合成して合成車両音の音圧波形を生成することを特徴とする請求項１から請求項５のいずれか一項に記載の音生成装置。
現行部品起因音演算手段と、部品起因音除去手段と、新規部品起因音演算手段と、合成音生成手段と、合成音発生手段と、部品回転角度計測手段と、部品回転速度計測手段とから構成され、
部品回転角度計測手段および前記部品回転速度計測手段は、車両に搭載された特定の部品の回転角度および回転速度を計測し、
前記現行部品起因音演算手段は、回転機械の駆動力によって駆動される特定の部品が起因で発生する現行部品起因音の音圧データを前記回転機械の回転数に比例した周波数である回転次数成分ごとの音圧次数成分データで演算し、
前記新規部品起因音演算手段は、前記回転機械の回転駆動力により部品回転軸が回転されることで部品が起因で発生する新規部品起因音の音圧データを前記部品回転軸の回転数に比例する周波数である部品回転次数成分ごとの音圧次数成分データで演算し、
前記合成音生成手段は、前記部品回転角度計測手段と部品回転速度計測手段とから得られた前記特定の部品の部品回転軸の回転角度および回転速度から現行部品起因音とは逆位相の音圧波形を生成し、前記新規部品音の音圧次数成分データから前記新規部品音の音圧波形を前記部品回転軸の初期位置を任意に変更した音圧波形を生成し、前記現行部品起因音とは逆位相の音圧波形と前記新規部品起因音の音圧波形を足し合わせた合成音の音圧波形を生成し、
前記合成音発生手段は、前記合成音の音圧波形から音を発生することを特徴とする音生成装置。
前記現行部品起因音または前記新規部品起因音が、部品で発生する加振力が複数の伝達経路を伝播して受音点に達するものであり、
前記現行部品起因音演算手段または前記新規部品起因音演算手段は、前記伝達経路への入力と前記伝達経路の伝達関数を乗じることで、前記現行部品起因音または前記新規部品起因音を算出し、前記伝達経路への入力と前記伝達関数を計測して求めることを特徴とする請求項１から請求項７のいずれか一項に記載の音生成装置。
前記現行部品起因音演算手段または前記新規部品起因音演算手段は、前記伝達関数への入力と前記伝達経路の伝達関数のうち少なくとも一方が計算機上で演算されることを特徴とする請求項１から請求項７のいずれか一項に記載の音生成装置。