車両内に設置されたセンサによりエンジンの振動を検出する工程と、前記エンジンの諸元値からエンジンの機構モデルを作成する工程と、前記機構モデルに筒内圧のモデル波形を入力してエンジン振動を出力する工程と、前記出力されたエンジン振動と前記センサにより検出されたエンジン振動との振動レベルを比較する工程と、前記比較した振動レベルを基に筒内圧のモデル波形の係数を変化させて、筒内圧を推定する工程と、を含むことを特徴とするエンジンの起振力推定方法。車両内に設置されたセンサによりエンジンの共振周波数を検出する工程と、前記機構モデルから得られるエンジンの共振周波数と前記センサにより検出されたエンジンの共振周波数とを比較して、前記エンジンの諸元値を変更する工程と、を更に含むことを特徴とする請求項1に記載のエンジンの起振力推定方法。しかしながら、エンジンの筒内圧を計測してエンジンの起振力を推定する方法は、実験車もしくはエンジンを用意しなければならず、コストがかかるという問題がある。また、エンジン筒内にセンサを設置するため、エンジンに穴を開ける必要があり、多くの工数がかかるという問題がある。本発明は、このような問題点に鑑みてなされたものであり、本発明の目的は、実験車のエンジンの筒内圧を計測せずに、エンジンの起振力を推定することのできるエンジンの起振力推定方法を提供することにある。上記目的を達成するため、本発明のエンジンの起振力推定方法は、車両内に設置されたセンサによりエンジンの振動を検出する工程と、前記エンジンの諸元値からエンジンの機構モデルを作成する工程と、前記機構モデルに筒内圧のモデル波形を入力してエンジン振動を出力する工程と、前記出力されたエンジン振動と前記センサにより検出されたエンジン振動との振動レベルを比較する工程と、前記比較した振動レベルを基に筒内圧のモデル波形の係数を変化させて、筒内圧を推定する工程と、を含むことを特徴とする。本発明のエンジンの起振力推定方法は、車両内に設置されたセンサによりエンジンの共振周波数を検出する工程と、前記機構モデルから得られるエンジンの共振周波数と前記センサにより検出されたエンジンの共振周波数とを比較して、前記エンジンの諸元値を変更する工程と、を更に含むことが好ましい。本発明は、実測のエンジン振動と計算で求めたエンジン振動の振動レベルを比較し、その振動レベルに基づいて筒内圧の係数を変更してエンジンの起振力を推定するので、実験車のエンジンの筒内圧を計測せずに、エンジンの起振力を推定することができる。また、本発明は、実験車のエンジンの共振周波数と、機構モデルのエンジンの共振周波数を一致させることで、より正確にエンジンの起振力を推定することができる。本発明の実施の形態について、図面を参照して説明する。図1は、本発明のエンジンの起振力推定方法を実施するための起振力推定システムの構成を示す図である。図1に示す起振力推定システムは、実験車のエンジンの共振周波数固有値を検出するための実測固有値取得手段以下、センサという11と、実験車のエンジンの振動を検出するための実測振動取得手段以下、センサという12と、コンピュータ10とからなり、センサ11、12は、コンピュータ10に接続されている。コンピュータ10は、機構モデル作成手段13と、固有値計算手段14と、固有値比較調整手段15と、モデル波形作成手段16と、振動計算手段17と、振動比較調整手段18とを備えており、これらの各手段は、ソフトウェアプログラムを実行することにより実現される。また、コンピュータ10は、データを入力するための入力手段図示せずと、ソフトウェアプログラムを格納するための記憶手段図示せずを備えている。センサ11、12は、起振部に近い箇所であり、エンジンの挙動を代表する箇所であるエンジンブロックの端や、エンジンマウント周辺等に設置される。センサ11、12には、加速度センサ、変位センサ等が用いられ、センサ11、12は、1つのセンサで代用してもよい。センサ11、12は、エンジンの負荷を統一するために、十分にエンジンの暖機運転を行って回転数が安定した状態で、エンジンの加速度または変位とエンジンの回転数を計測する。センサ11実測固有値取得手段は、モーダル解析を用いて共振周波数を求めるのが好ましい。機構モデル作成手段13は、ピストン、コンロッド、クランクシャフト、エンジンブロック、およびエンジンマウントからなるエンジン機構モデルを作成する。固有値計算手段14は、作成した機構モデルに、エンジンの質量、慣性モーメント、重心位置等の諸元値特性を入力して、機構モデルにおけるエンジンの共振周波数を計算する。固有値比較調整手段15は、センサ11より得られた実験車のエンジンの共振周波数と機構モデルのエンジンの共振周波数とを比較して、共振周波数が一致するように機構モデルのエンジンの諸元値特性を調整する。モデル波形作成手段16は、機構モデルを動作させるためのモデル波形を作成する。モデル波形は、機構モデル上でエンジンの筒内圧波形を力要素として設定し、シャフト軸に抵抗となる力要素を設定して作成する。シャフト軸に抵抗となる力要素は、熱力学サイクルをベースとしたPV線図を積分することで、回転トルクを算出し、算出した回転トルクNmと、目標回転数rad/s例えば、実測したアイドル振動をターゲットとする場合、そのときの回転数より負荷係数Nms/radを算出し、負荷係数に現在のエンジン回転数を乗算して算出する。振動計算手段17は、機構モデルにモデル波形を与えて機構モデルにおけるエンジンの振動の時系列データを計算する。振動比較調整手段18は、計算した振動の時系列データと、センサ12より得られた実験車のエンジンの振動の時系列データとを比較して、時系列データが一致するようにエンジンの筒内圧を調整する。ここで、モデル波形筒内圧とエンジン振動との関係を簡単に説明する。エンジンは気筒内の爆発を回転エネルギーに変換する機構であり、その際、様々な不平衡によりエンジンには振動が発生する。剛体近似でのエンジンに働く力は、ピストンが動くことによる慣性力と、エンジンの回転変動により働くトルクである。この反力がエンジンにかかり、エンジンの振動の起振力となる。本発明では、上記トルクに代表される起振力のうち、慣性力に起因する起振力は、機構モデルで自動的に計算されるという前提のもと、残りのエンジンの回転変動によるトルクに起因する起振力は、エンジン内の爆発で発生した内力筒内圧によるトルクの変動が支配的であると想定し、また、筒内圧がエンジンの振動レベルに比例すると近似して、筒内圧を、エンジン振動係数の比較により推定している。次に、本発明のエンジンの起振力推定方法について図2のフローチャートに基づいて説明する。まず、機構モデル作成手段13は、ピストン、コンロッド、クランクシャフト、エンジンブロック、およびエンジンマウントからなるエンジン機構モデルを作成するS101。次に、図示しない入力手段から、機構モデルに、エンジンの質量、慣性モーメント、重心位置等の諸元値特性が入力されるとS102、固有値計算手段14は、機構モデルにおけるエンジンの共振周波数を計算するS103。固有値比較調整手段15は、センサ11実測固有値取得手段から実験車のエンジンの共振周波数を取得するとS104、実験車のエンジンの共振周波数と、計算で求めた機構モデルのエンジンの共振周波数とが一致するか否かを判定するS105。共振周波数が一致しない場合S105でNoの場合は、実験車のエンジンの共振周波数と一致するように、機構モデルの車両諸元値を変更しS106、S103に戻る。共振周波数が一致する場合S105でYesの場合は、モデル波形作成手段16は、図示しない入力手段から、エンジンの筒内圧波形とシャフト軸の抵抗が力要素として入力されると、機構モデルを動作させるためのモデル波形を作成するS107。次に、振動計算手段17は、機構モデルの釣り合いを計算した後、機構モデルにモデル波形を与えて機構モデルにおけるエンジンの振動の時系列データを計算するS108。振動比較調整手段18は、センサ12実測固有値取得手段から実験車のエンジンの振動の時系列データを取得するとS109、実験車の振動の時系列データのピークと、計算で求めた機構モデルの振動の時系列データのピークとが一致するか否かを判定するS110。振動の時系列データのピークが一致しない場合S110でNoの場合は、計算で求めた振動のピークと、実験車の振動のピークとの係数を求め、例えば、計算で求めた振動のピークと、実験車の振動のピークとがA倍ずれていれば、計算で求めた機構モデルのモデル波形に1/Aを掛け、それを筒内圧としS111、S107に戻る。振動の時系列データのピークが一致した場合S110でYesの場合、すなわち、係数が一致した場合は、そのときの筒内圧を起振力と推定するS112。以後、この推定された起振力を用いて車両や車両部品の開発に供する解析モデルを構築する。図3は、エンジンマウント右側に設置したセンサ12により得られた上下方向の加速度の時系列データと、エンジンマウント右側の位置において計算で求めた上下方向の加速度の係数違いの時系列データである。図4は、エンジンマウント右側に設置したセンサ12により得られた左右方向の加速度の時系列データと、エンジンマウント右側の位置おいて計算で求めた左右方向の加速度の係数違いの時系列データである。本発明の起振力推定方法は、図3および図4に示すように、計算で求めた振動データが、実験車の振動データ実測データと係数が異なる場合、係数が一致するように、その係数に基づいて筒内圧を変更してエンジンの起振力を推定するものである。