(57)【特許請求の範囲】
【請求項１】
エンジンを制御する複数の制御パラメータにおいて、当該制御パラメータの限界値によって構成される制御可能領域の範囲内における、当該制御パラメータの組み合わせである計測点を算出する計測点算出装置であって、
乱数を発生する乱数発生手段と、
前記制御可能領域の範囲内に、第１計測点を決定する第１計測点決定手段と、
前記乱数発生手段によって発生させた乱数に基づいて前記制御可能領域の範囲内に仮計測点を算出する仮計測点算出手段と、
前記第１計測点決定手段によって決定された前記第１計測点から、前記仮計測点算出手段によって算出された前記仮計測点へ結んだ直線の延長線と、前記制御可能領域との交点を算出する交点算出手段と、
前記乱数発生手段によって発生させた乱数に基づいて、前記交点算出手段によって結ばれた前記直線上の点であって前記仮計測点と前記交点との間の第２計測点を算出する第２計測点算出手段と、
前記第２計測点算出手段によって算出された前記第２計測点を前記第１計測点として、前記仮計測点算出手段と、前記交点算出手段と、前記第２計測点算出手段とが行う処理を繰り返す繰り返し手段と、
を備える計測点算出装置。
【請求項２】
前記第１計測点決定手段は、前記制御可能領域の重心を算出し、第１計測点とする請求項１に記載の計測点算出装置。
【請求項３】
前記第２計測点算出手段によって算出された計測点におけるエンジンの性能を測定する測定手段をさらに備える請求項１又は２に記載の計測点算出装置。
【請求項４】
エンジンを制御する複数の制御パラメータにおいて、当該制御パラメータの限界値によって構成される制御可能領域の範囲内における、当該制御パラメータの組み合わせである計測点を、乱数発生手段を有する計測点算出装置が算出する方法であって、
前記制御可能領域の範囲内に、第１計測点を決定する第１計測点決定ステップと、
前記乱数発生手段によって発生させた乱数に基づいて前記制御可能領域の範囲内に仮計測点を算出する仮計測点算出ステップと、
前記第１計測点決定ステップによって決定された前記第１計測点から、前記仮計測点算出ステップによって算出された前記仮計測点へ結んだ直線の延長線と、前記制御可能領域との交点を算出する交点算出ステップと、
前記乱数発生手段によって発生させた乱数に基づいて、前記交点算出ステップによって結ばれた前記直線上の点であって前記仮計測点と前記交点との間の第２計測点を算出する第２計測点算出ステップと、
前記第２計測点算出ステップによって算出された前記第２計測点を前記第１計測点として、前記仮計測点算出ステップと、前記交点算出ステップと、前記第２計測点算出ステップとが行う処理を繰り返す繰り返しステップと、
を備える方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計測点算出装置及び方法に関する。
【背景技術】
【０００２】
近年、エンジンは、高性能化により、多くの制御パラメータによって制御されて動作する。この制御パラメータによって制御されるエンジンの特性は、実験計画法（ＤＯＥ）によって配置された所定数の計測点における計測データに基づいて、モデル化されている。このような計測データに基づいたモデルを使用して、エンジンの性能を最適にする制御パラメータの適合値が算出される（非特許文献１及び非特許文献２参照）。
【０００３】
例えば、空間充填法（Ｓｐａｃｅ  Ｆｉｌｌｉｎｇ）を応用して、制御パラメータの適合値について、実験計画を生成する発明を開示する特許文献１が知られている。
【０００４】
特許文献１で開示する適合装置は、所定数の計測点を実験計画法（ＤＯＥ）のスペースフィリング（Ｓｐａｃｅ  Ｆｉｌｌｉｎｇ）を使用して配置した後、複数の制御パラメータのうちの少なくとも１つの制御パラメータに関して計測点の位置を当該制御パラメータの中心値から境界方向へ拡張し、この拡張により所定数の計測点のうち境界を超えたものを境界上に再配置するようにしている。これにより、全体の計測点数を増加させなくても、境界付近のモデル誤差を目標レベルに収めることが可能となり、モデル精度向上の要求と適合時間（計測時間）短縮化の要求とを共に満たすことができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００６−２４４０４２号公報
【非特許文献】
【０００６】
【非特許文献１】Ｒ．Ｌ．Ｓｍｉｔｈ，Ｅｆｆｉｃｉｅｎｔ  Ｍｏｎｔｅ−Ｃａｒｌｏ  ｐｒｏｃｅｄｕｒｅｓ  ｆｏｒ  ｇｅｎｅｒａｔｉｎｇ  ｐｏｉｎｔｓ  ｕｎｉｆｏｒｍｌｙ  ｄｉｓｔｒｉｂｕｔｅｄ  ｏｖｅｒ  ｂｏｕｎｄｅｄ  ｒｅｇｉｏｎｓ，Ｏｐｅｒａｔｉｏｎｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ３２，ｐｐ．１２９６−１３０８，１９８４．
【非特許文献２】Ｌ．Ｌｏｖａｓｚ，Ｈｉｔ−ａｎｄ−ｒｕｎ  ｍｉｘｅｓ  ｆａｓｔ，  Ｍａｔｈｅｍａｔｉｃａｌ  Ｐｒｏｇｒａｍｍｉｎｇ，ｖｏｌ．８６，ｐｐ．４４３−４６１，１９９８．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１に開示された発明における空間充填法は、多くの制御パラメータ（例えば、１０以上）があり、かつ複雑な凸境界において所定数の計測点を配置する場合、多大な時間がかかり、実用的でない。
【０００８】
そこで、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる装置が求められている。
【０００９】
本発明は、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる計測点算出装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明では、以下のような解決手段を提供する。
【００１１】
（１）  エンジンを制御する複数の制御パラメータにおいて、当該制御パラメータの限界値によって構成される制御可能領域の範囲内における当該制御パラメータの組み合わせである計測点を算出する計測点算出装置であって、乱数を発生する乱数発生手段と、前記制御可能領域の範囲内に、第１計測点を決定する第１計測点決定手段と、前記乱数発生手段によって発生させた乱数に基づいて前記制御可能領域の範囲内に仮計測点を算出する仮計測点算出手段と、前記第１計測点決定手段によって決定された前記第１計測点から、前記仮計測点算出手段によって算出された前記仮計測点へ結んだ直線の延長線と、前記制御可能領域との交点を算出する交点算出手段と、前記乱数発生手段によって発生させた乱数に基づいて、前記交点算出手段によって結ばれた前記直線上の点であって前記仮計測点と前記交点との間の第２計測点を算出する第２計測点算出手段と、前記第２計測点算出手段によって算出された前記第２計測点を前記第１計測点として、前記仮計測点算出手段と、前記交点算出手段と、前記第２計測点算出手段とが行う処理を繰り返す繰り返し手段と、を備える計測点算出装置。
【００１２】
（１）の構成によれば、本発明に係る計測点算出装置は、乱数を発生する乱数発生手段を有する。そして、本発明に係る計測点算出装置は、制御可能領域の範囲内に、第１計測点を決定し、乱数発生手段によって発生させた乱数に基づいて制御可能領域の範囲内に仮計測点を算出し、決定された第１計測点から、算出された仮計測点の方向へ結んだ直線の延長線と、制御可能領域との交点を算出する。次に、本発明に係る計測点算出装置は、乱数発生手段によって発生させた乱数に基づいて、交点と結んだ直線上の点であって仮計測点と交点との間の第２計測点を算出する。そして、本発明に係る計測点算出装置は、算出した第２計測点を第１計測点として、仮計測点の算出と、交点の算出と、第２計測点の算出とを繰り返し、計測点を算出する。
【００１３】
すなわち、本発明に係る計測点算出装置は、第１計測点から、乱数に基づく仮計測点に基づいて第２計測点を算出し、第２計測点の算出を繰り返して計測点を算出する。よって、本発明に係る計測点算出装置は、第１計測点からランダムに直線の方向を変えながら、制御可能領域内の計測点をランダムに算出する。したがって、本発明に係る計測点算出装置は、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる。
【００１４】
（２）  前記第１計測点決定手段は、前記制御可能領域の重心を算出し、第１計測点とする（１）に記載の計測点算出装置。
【００１５】
したがって、（２）に係る計測点算出装置は、制御可能領域内において、偏ることなく計測点を算出することができる。
【００１６】
（３）  前記第２計測点算出手段によって算出された計測点におけるエンジンの性能を測定する測定手段をさらに備える（１）又は（２）に記載の計測点算出装置。
【００１７】
したがって、（３）に係る計測点算出装置は、妥当な時間内で算出された、十分な数の計測点により、エンジンの性能を測定することができる。
【００１８】
（４）  エンジンを制御する複数の制御パラメータにおいて、当該制御パラメータの限界値によって構成される制御可能領域の範囲内における、当該制御パラメータの組み合わせである計測点を、乱数発生手段を有する計測点算出装置が算出する方法であって、前記制御可能領域の範囲内に、第１計測点を決定する第１計測点決定ステップと、前記乱数発生手段によって発生させた乱数に基づいて前記制御可能領域の範囲内に仮計測点を算出する仮計測点算出ステップと、前記第１計測点決定ステップによって決定された前記第１計測点から、前記仮計測点算出ステップによって算出された前記仮計測点へ結んだ直線の延長線と、前記制御可能領域との交点を算出する交点算出ステップと、前記乱数発生手段によって発生させた乱数に基づいて、前記交点算出ステップによって結ばれた前記直線上の点であって前記仮計測点と前記交点との間の第２計測点を算出する第２計測点算出ステップと、前記第２計測点算出ステップによって算出された前記第２計測点を前記第１計測点として、前記仮計測点算出ステップと、前記交点算出ステップと、前記第２計測点算出ステップとが行う処理を繰り返す繰り返しステップと、を備える方法。
【００１９】
（４）の構成によれば、本発明に係る方法は、第１計測点から、乱数に基づく仮計測点に基づいて第２計測点を算出し、第２計測点の算出を繰り返して計測点を算出する。よって、本発明に係る方法は、第１計測点からランダムに直線の方向を変えながら、制御可能領域内の計測点をランダムに算出する。したがって、本発明に係る方法は、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる。
【発明の効果】
【００２０】
本発明によれば、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる計測点算出装置及び方法を提供することができる。
【図面の簡単な説明】
【００２１】
【図１】本発明の特徴を示す図である。
【図２】本発明の一実施形態に係る計測点算出装置の機能を示す機能ブロック図である。
【図３】本発明の一実施形態に係る計測点算出装置のハードウェア構成の一例を示す図である。
【図４】本発明の一実施形態に係る計測点算出装置の処理内容を示すフローチャートである。
【図５】本発明の一実施形態に係る計測点算出装置の測定処理の内容を示すフローチャートである。
【図６】本発明の一実施形態に係る計測点算出装置による計測点算出の例を示す図である。
【発明を実施するための形態】
【００２２】
以下、本発明の実施形態について図を参照しながら説明する。
【００２３】
図１は、本発明の特徴を示す図である。計測点算出装置１０は、エンジンを制御する複数の制御パラメータにおいて、制御パラメータの限界値に対応する限界点１０１，１０２，１０３，１０４，１０５，１０６によって構成される制御可能領域１９０の範囲内における、制御パラメータの組み合わせである計測点を算出する。例えば、計測点算出装置１０は、制御可能領域１９０の範囲内に、第１計測点２０１を決定し、乱数発生部２１（後述する図２参照）によって発生させた乱数に基づいて制御可能領域１９０の範囲内に仮計測点３０１を算出し、決定された第１計測点２０１から、算出された仮計測点３０１へ結んだ直線１５１と、制御可能領域１９０との交点１１１を算出する。次に、計測点算出装置１０は、乱数発生部２１によって発生させた乱数に基づいて、第１計測点２０１と交点１１１とを結んだ直線上の点であって仮計測点３０１と交点１１１との間の第２計測点２０２を算出する。そして、計測点算出装置１０は、算出された第２計測点２０２を第１計測点として、仮計測点３０２の算出と、直線１５２と制御可能領域１９０との交点１１２の算出と、第２計測点２０３の算出とを繰り返し、計測点を算出する。このように、計測点算出装置１０は、第１計測点からランダムに直線の方向を変えながら、制御可能領域１９０内の計測点をランダムに算出する。したがって、計測点算出装置１０は、多くの制御パラメータによる複雑な凸境界で構成された制御可能領域１９０であっても、妥当な時間内に、十分な数の計測点を算出することができる。
【００２４】
図２は、本発明の一実施形態に係る計測点算出装置１０の機能を示す機能ブロック図である。計測点算出装置１０は、乱数発生手段としての乱数発生部２１と、第１計測点決定手段としての第１計測点決定部１１と、仮計測点算出手段としての仮計測点算出部１２と、交点算出手段としての交点算出部１３と、第２計測点算出手段としての第２計測点算出部１４と、繰り返し手段としての繰り返し部１５と、測定手段としての測定部１６とを備える。そして、計測点算出装置１０は、エンジンを制御する複数の制御パラメータにおいて、制御パラメータの限界値によって構成される制御可能領域の範囲内における制御パラメータの組み合わせである計測点を算出する。このような計測点算出装置１０について、各部ごとに詳述する。
【００２５】
乱数発生部２１は、乱数を発生する。乱数発生部２１は、信号を元に乱数を生成するハードウェア乱数生成器や、ＣＰＵ１０１０（後述する図３参照）によって計算によって作成される擬似乱数を生成する擬似乱数生成器であってよい。
【００２６】
第１計測点決定部１１は、制御可能領域の範囲内に、第１計測点を決定する。制御可能領域は、制御パラメータの限界値による条件式等で表わされる凸境界によって構成される。最初の第１計測点は、例えば、制御可能領域内の重心や、制御可能領域内の中心点等である。
【００２７】
仮計測点算出部１２は、乱数発生部２１によって発生させた乱数に基づいて制御可能領域の範囲内に仮計測点を算出する。例えば、仮計測点算出部１２は、発生させた乱数を中心角と半径とに対応付けて、第１計測点を起点とする仮計測点を算出する。算出した点が制御可能領域外である場合、仮計測点算出部１２は、制御可能領域内の点に変換して制御可能領域内の仮計測点を算出する。
【００２８】
交点算出部１３は、第１計測点決定部１１によって決定された第１計測点から、仮計測点算出部１２によって算出された仮計測点へ結んだ直線の延長線と、制御可能領域との交点を算出する。例えば、交点算出部１３は、第１計測点から仮計測点へ結んだ直線の方程式と、制御可能領域を表わす条件式とに基づいて交点の解を算出する。
【００２９】
第２計測点算出部１４は、乱数発生部２１によって発生させた乱数に基づいて、交点算出部１３によって結ばれた直線上の点であって仮計測点と交点との間の第２計測点を算出する。例えば、第２計測点算出部１４は、発生させた乱数を仮計測点と交点とを内分する比率に対応付け、仮計測点と交点との間を内分する第２計測点を算出する。
【００３０】
繰り返し部１５は、第２計測点算出部１４によって算出された第２計測点を第１計測点として、仮計測点算出部１２と、交点算出部１３と、第２計測点算出部１４とが行う処理を繰り返し、計測点を算出する。すなわち、計測点算出装置１０は、算出された計測点を起点として、ランダムに方向を変えながら、制御可能領域内の計測点をランダムに次々と算出する。計測点算出装置１０は、算出した計測点のうち、数点おきに抽出して測定を行う計測点としてもよい。
【００３１】
測定部１６は、第２計測点算出部１４によって算出された計測点におけるエンジンの性能を測定する。例えば、測定部１６は、エンジンを制御するＥＣＵ（Ｅｌｅｃｔｒｏｎｉｃ  Ｃｏｎｔｒｏｌ  Ｕｎｉｔ）制御システムに、算出した計測点に対応する制御パラメータの値を送信し、ＥＣＵ制御システムから測定結果を受信する。
【００３２】
図３は、本発明の一実施形態に係る計測点算出装置１０のハードウェア構成の一例を示す図である。計測点算出装置１０は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１０１０、バスライン１００５、通信Ｉ／Ｆ１０４０、メインメモリ１０５０、ＢＩＯＳ（Ｂａｓｉｃ  Ｉｎｐｕｔ  Ｏｕｔｐｕｔ  Ｓｙｓｔｅｍ）１０６０、Ｉ／Ｏコントローラ１０７０、キーボード／マウス１１００及び表示装置１０２２を備える。
【００３３】
Ｉ／Ｏコントローラ１０７０には、ハードディスク１０７４、半導体メモリ１０７８、等の記憶手段を接続することができる。
【００３４】
ＢＩＯＳ１０６０は、計測点算出装置１０の起動時にＣＰＵ１０１０が実行するブートプログラムや、計測点算出装置１０のハードウェアに依存するプログラム等を格納する。
【００３５】
ハードディスク１０７４は、計測点算出装置１０が本発明の機能を実行するためのプログラムを記憶しており、さらに、各種データベースを構成可能である。
【００３６】
計測点算出装置１０に提供されるプログラムは、ハードディスク１０７４、又はメモリカード等の記録媒体に格納されて提供される。このプログラムは、Ｉ／Ｏコントローラ１０７０を介して、記録媒体から読み出され、又は通信Ｉ／Ｆ１０４０を介してダウンロードされることによって、計測点算出装置１０にインストールされ実行されてもよい。
【００３７】
前述のプログラムは、専用通信回線に接続されたサーバシステムに設けたハードディスク又は光ディスクライブラリ等の記憶装置を記録媒体として使用し、通信回線を介して計測点算出装置１０に提供されるとしてもよい。
【００３８】
ここで、表示装置１０２２は、計測点算出装置１０による演算処理結果の画面を表示したりするものであり、ブラウン管表示装置（ＣＲＴ）、液晶表示装置（ＬＣＤ）等のディスプレイ装置を含む。
【００３９】
通信Ｉ／Ｆ１０４０は、計測点算出装置１０を専用ネットワークを介して端末（例えば、測定対象の装置やＥＣＵ制御システム等）と接続できるようにするためのネットワーク・アダプタである。
【００４０】
図４は、本発明の一実施形態に係る計測点算出装置１０の処理内容を示すフローチャートである。なお、本処理は、例えば、プログラム開始指令を受け付けて開始し、プログラム終了指令又は終了条件により終了する。
【００４１】
ステップＳ１０１において、ＣＰＵ１０１０は、第１計測点を決定する。より具体的には、ＣＰＵ１０１０は、制御可能領域を構成する限界点の座標に基づいて重心を算出し、第１計測点とする。その後、ＣＰＵ１０１０は、処理をステップＳ１０２に移す。
【００４２】
ステップＳ１０２において、ＣＰＵ１０１０は、乱数により仮計測点を算出する。より具体的には、ＣＰＵ１０１０は、乱数発生部２１から、例えば、４桁の自然数を取得する。そして、例えば制御パラメータが２パラメータの場合、ＣＰＵ１０１０は、０〜９９９９を０ラジアン〜２πラジアンに対応付け、取得した乱数から中心角を算出し、０〜９９９９を第１計測点と、第１計測点から最大距離の限界点とに対応付け、取得した乱数から半径を算出し、仮計測点の座標を算出する。そして、ＣＰＵ１０１０は、算出した仮計測点の座標が制御可能領域外である場合、制御可能領域内の座標に変換するために、例えば、制御可能領域の境界で半径の方向を反対方向に折り返した点を仮計測点とする。その後、ＣＰＵ１０１０は、処理をステップＳ１０３に移す。
【００４３】
ステップＳ１０３において、ＣＰＵ１０１０は、交点を算出する。より具体的には、ＣＰＵ１０１０は、第１計測点から仮計測点へ結んだ直線の方程式と、制御可能領域を表わす条件式とに基づいて交点の解を算出する。その後、ＣＰＵ１０１０は、処理をステップＳ１０４に移す。
【００４４】
ステップＳ１０４において、ＣＰＵ１０１０は、乱数により第２計測点を算出する。より具体的には、ＣＰＵ１０１０は、乱数発生部２１から、例えば、４桁の自然数を取得する。そして、ＣＰＵ１０１０は、０〜９９９９を仮計測点と交点とを内分する比率に対応付け、取得した乱数ｕにより仮計測点と交点との間をｕ：（９９９９−ｕ）に内分する第２計測点を算出する。その後、ＣＰＵ１０１０は、処理をステップＳ１０５に移す。
【００４５】
ステップＳ１０５において、ＣＰＵ１０１０は、算出した計測点の個数が十分か否かを判断する。すなわち、ＣＰＵ１０１０は、算出した計測点の個数が所定の個数に達したか否かを判断する。この判断がＹＥＳの場合、ＣＰＵ１０１０は、処理をステップＳ１０６に移し、ＮＯの場合、ＣＰＵ１０１０は、第２計測点を第１計測点として処理をステップＳ１０２に移す。
【００４６】
ステップＳ１０６において、ＣＰＵ１０１０は、算出した計測点に基づいて測定処理（後述する図５参照）を行う。その後、ＣＰＵ１０１０は、処理を終了する。
【００４７】
図５は、本発明の一実施形態に係る計測点算出装置１０の測定処理の内容を示すフローチャートである。
【００４８】
ステップＳ４０１において、ＣＰＵ１０１０は、制御パラメータの値を設定する。より具体的には、ＣＰＵ１０１０は、エンジンを動作させるシステム（例えば、ＥＣＵ制御システム）に、算出した計測点に基づく制御パラメータを送信する。その後、ＣＰＵ１０１０は、処理をステップＳ４０２に移す。
【００４９】
ステップＳ４０２において、ＣＰＵ１０１０は、測定を行う。より具体的には、ＣＰＵ１０１０は、エンジンを動作させるシステム（例えば、ＥＣＵ制御システム）から、動作結果としての測定結果を受信する。その後、ＣＰＵ１０１０は、処理をステップＳ４０３に移す。
【００５０】
ステップＳ４０３において、ＣＰＵ１０１０は、測定値が所定の動作条件を超えたか否かを判断する。この判断がＹＥＳの場合、ＣＰＵ１０１０は、処理をステップＳ４０５に移し、ＮＯの場合、ＣＰＵ１０１０は、処理をステップＳ４０４に移す。
【００５１】
ステップＳ４０４において、ＣＰＵ１０１０は、測定値が安定したか否かを判断する。すなわち、ＣＰＵ１０１０は、受信した測定結果が変化しないで一定の値を示しているか否かを判断する。この判断がＹＥＳの場合、ＣＰＵ１０１０は、処理をステップＳ４０５に移し、ＮＯの場合、ＣＰＵ１０１０は、処理をステップＳ４０２に移す。
【００５２】
ステップＳ４０５において、ＣＰＵ１０１０は、全ての計測点を測定したか否かを判断する。この判断がＹＥＳの場合、ＣＰＵ１０１０は、処理をステップＳ４０６に移し、ＮＯの場合、ＣＰＵ１０１０は、処理をステップＳ４０１に移す。
【００５３】
ステップＳ４０６において、ＣＰＵ１０１０は、測定値を返す。その後、ＣＰＵ１０１０は、処理を終了し、本処理に移るステップの次のステップに処理を戻す。
【００５４】
図６は、本発明の一実施形態に係る計測点算出装置１０による計測点算出の例を示す図である。図６において、計測点算出装置１０は、第１計測点２０１から、乱数に基づいて直線１５１上に第２計測点２０２を算出し、第２計測点２０２を次の起点として乱数に基づいて直線１５２上に次の計測点２０３を算出し、算出を繰り返して直線１５３，１５４，１５５，１５６上に計測点２０４，２０５，２０６，２０７を算出する。図６の例は、このようにして計測点算出装置１０が、計測点の算出を繰り返し、多くの制御パラメータによる複雑な凸境界で構成された制御可能領域１９０であっても、妥当な時間内に、領域内に一様に分布する、十分な数の計測点を算出することを示す例である。
【００５５】
本実施形態によれば、計測点算出装置１０は、乱数を発生する乱数発生部２１を有する。そして、計測点算出装置１０は、制御可能領域の範囲内に、第１計測点を決定し、乱数発生部２１によって発生させた乱数に基づいて制御可能領域の範囲内に仮計測点を算出し、決定された第１計測点から、算出された仮計測点の方向へ結んだ直線の延長線と、制御可能領域との交点を算出する。次に、計測点算出装置１０は、乱数発生手段によって発生させた乱数に基づいて、交点と結んだ直線上の点であって仮計測点と交点との間の第２計測点を算出する。そして、計測点算出装置１０は、算出した第２計測点を第１計測点として、仮計測点の算出と、交点の算出と、第２計測点の算出とを繰り返し、計測点を算出する。したがって、計測点算出装置１０は、多くの制御パラメータによる複雑な凸境界であっても、妥当な時間内に、十分な数の計測点を算出することができる。そして、計測点算出装置１０は、算出した十分な数の計測点におけるエンジンの性能を測定する。
【００５６】
以上、本発明の実施形態について説明したが、本発明は上述した実施形態に限るものではない。また、本発明の実施形態に記載された効果は、本発明から生じる最も好適な効果を列挙したに過ぎず、本発明による効果は、本発明の実施形態に記載されたものに限定されるものではない。
【符号の説明】
【００５７】
１０  計測点算出装置
１１  第１計測点決定部
１２  仮計測点算出部
１３  交点算出部
１４  第２計測点算出部
１５  繰り返し部
１６  測定部
２１  乱数発生部
【図１】
0005443283.tif 000002
【図２】
0005443283.tif 000003
【図３】
0005443283.tif 000004
【図４】
0005443283.tif 000005
【図５】
0005443283.tif 000006
【図６】
0005443283.tif 000007