(57)【特許請求の範囲】
【請求項１】
機関負荷と機関回転数によって定まる内燃機関の各運転状態において該内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であって、
上記制御パラメータの値を様々な値に設定して該制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて該特性パラメータと上記制御パラメータとの関係を全運転状態に渡って求め、少なくとも該関係に基づいて各運転状態における上記制御パラメータの適合値の存在範囲を推定し、次いで該存在範囲内において上記制御パラメータの値を様々な値に設定して該制御パラメータの各値毎に上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて該特性パラメータと上記制御パラメータとの各運転状態における関係を求め、該関係に基づいて各運転状態における上記制御パラメータの適合値を求める、内燃機関の制御パラメータの適合値を求める方法。
【請求項２】
機関負荷と機関回転数によって定まる内燃機関の各運転状態において該内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であって、
上記制御パラメータの値を様々な値に設定して該制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータ、並びに機関負荷と機関回転数をパラメータとして含むモデル式である全域モデルのモデル式を求め、該全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定することと、
内燃機関の各運転状態毎に、内燃機関の運転状態がその運転状態になる場合において上記制御パラメータの値を上記適合値の存在範囲内で様々な値に設定した場合における上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータをパラメータとして含む一方、機関負荷と機関回転数をパラメータとして含まないモデル式である局所モデルのモデル式を求め、該局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求めることと、を含む内燃機関の制御パラメータの適合値を求める方法。
【請求項３】
上記制御パラメータの適合値の存在範囲が、上記全域モデルのモデル式から得られる特性パラメータの値の他、上記全域モデルのモデル式を求める際に上記特性パラメータの計測を実施した上記制御パラメータの範囲と、上記全域モデルのモデル式を求める際のモデル化誤差と、上記制御パラメータを制御する精度と、上記制御パラメータを計測する精度と、上記特性パラメータを計測する精度とのうちの少なくとも一つに基づいて推定される、請求項２に記載の内燃機関の制御パラメータの適合値を求める方法。
【請求項４】
機関負荷と機関回転数によって定まる内燃機関の各運転状態において該内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であり、
上記制御パラメータの値を様々な値に設定して該制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて、特性パラメータの値を表すモデル式であって上記制御パラメータ、並びに機関負荷と機関回転数をパラメータとして含むモデル式である全域モデルのモデル式を求め、該全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求めると共に、
内燃機関の運転状態のうちの予め定められた特定運転状態については、内燃機関の運転状態が上記特定運転状態になる場合において上記制御パラメータの値を様々な値に設定した場合における上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて、特性パラメータの値を表すモデル式であって上記制御パラメータをパラメータとして含む一方、機関負荷と機関回転数をパラメータとして含まないモデル式である局所モデルのモデル式を求め、該局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を再度求めるようになっていて、
上記特定運転状態以外の運転状態に関しては上記全域モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値を最終的な適合値とする一方、上記特定運転状態に関しては上記局所モデルのモデル式から得られる特性パラメータの値に基づいて再度求められた上記制御パラメータの適合値を最終的な適合値とする、制御パラメータの適合値を求める方法であって、
上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲は、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて決定される、内燃機関の制御パラメータの適合値を求める方法。
【請求項５】
上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲が、上記全域モデルのモデル式から得られる特性パラメータの値の他、上記全域モデルのモデル式を求める際に上記特性パラメータの計測を実施した上記制御パラメータの範囲と、上記全域モデルのモデル式を求める際のモデル化誤差と、上記制御パラメータを制御する精度と、上記制御パラメータを計測する精度と、上記特性パラメータを計測する精度とのうちの少なくとも一つに基づいて決定される、請求項４に記載の内燃機関の制御パラメータの適合値を求める方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は内燃機関の適合方法、すなわち内燃機関の制御において用いられる制御パラメータの適合値を求める方法に関する。
【背景技術】
【０００２】
一般に、内燃機関の制御は、トルク、排気エミッション及び燃費等の内燃機関の特性についての要求条件を満たすようにスロットル開度、点火時期、吸気弁又は排気弁の開閉弁特性、燃料噴射量等の制御パラメータの値を変化させることによって行われる。そしてこれらの制御パラメータには、そのときの内燃機関に対する要求に応じて、例えば機関負荷（トルク）及び機関回転数により定まる運転状態毎に最適な値が存在するため、このような運転状態毎の制御パラメータの最適な値を目標値として予め設定しておく必要がある。
【０００３】
このような運転状態毎の制御パラメータの最適な値は、一般に、上記制御パラメータを様々な値に設定し、そのときの発生トルクやＮＯX排出量等の内燃機関の特性を表す特性パラメータを計測してその結果から運転状態毎の各制御パラメータの最適な値（すなわち、適合値）を求める作業、いわゆる適合作業によって求められる。
【０００４】
斯かる適合作業においては、定常運転時における特性パラメータの値を計測する必要があることから、制御パラメータの値を変化させて設定される各計測点毎に運転状態が安定するまで待ってから、すなわち例えばトルク、機関回転数がほぼ一定の値に収束するまで待ってから計測が行われる。このため、各計測点において特性パラメータの計測値を得るまでに長い時間を要する。また、適合精度を高いものとするためには、多くの計測点において計測が必要であり、場合によっては計測点数が数千〜数十万点にも及ぶ。このため、適合作業全体の計測工数は膨大なものとなる。
【０００５】
そこで、各計測点間の間隔を広げ、すなわち各計測点間の制御パラメータの値の差を大きくして計測点数を低減すると共に、特性パラメータの計測値に基づいてモデル式を求め、斯かるモデル式に基づいて計測点間の各制御パラメータの値に対する特性パラメータの値を推定することが提案されている（特許文献１）。このようにすると、計測工数を低減して各制御パラメータの最適な値（すなわち、適合値）を求めることが可能となる。
【０００６】
【特許文献１】特開２００２−２０６４５６号公報
【特許文献２】特開２００４−６８７２９号公報
【特許文献３】特開２０００−２４８９９１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
ところで、上記特性パラメータの値を表すモデル式として各制御パラメータに加えて機関負荷及び機関回転数もパラメータとして含むモデル式を用い、実際に運転が行われると想定される全ての運転状態（以下、単に「全運転状態」という）を一つのモデル（以下、「全域モデル」という）で表現するようにする場合がある。そしてこのようにする場合には比較的少ない計測工数で全運転状態における各制御パラメータの適合値を求めることが可能となる。
【０００８】
ところが、この場合に上記各制御パラメータの適合値として求められる値は、実際にはその精度が充分でない場合がある。この原因の一つは上記全域モデルを構築するためには全運転状態に渡って計測を実施する必要があるため、計測期間が相当長期間となり、その結果気象条件や内燃機関の状態等の計測条件に変化が生じてモデル精度が低下することにある。
【０００９】
一方、上記特性パラメータの値を表すモデル式として各運転状態毎に機関負荷と機関回転数をパラメータとして含まないモデル式（すなわち、機関負荷と機関回転数を固定した場合のモデル式）を求めるようにし、このモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求める方法もある。すなわち、この場合には機関負荷と機関回転数をパラメータとして含まない（すなわち、機関負荷と機関回転数を固定した場合の）モデル（以下、「局所モデル」という）が構築されて、各制御パラメータの適合値が求められる。
【００１０】
この局所モデルを用いた場合には、一般に上記全域モデルを用いた場合に比べて上記制御パラメータの適合値を精度よく求めることが出来るのであるが、各運転状態毎にモデル式を求めるのに充分な数の特性パラメータの計測値が必要となるため、計測工数は上記全域モデルを用いた場合に比べて相当に大きくなってしまう。
【００１１】
本発明は以上のような点に鑑みてなされたものであり、その目的は計測工数の増大を抑えつつ、内燃機関の制御パラメータの適合値を精度よく求める方法を提供することである。
【課題を解決するための手段】
【００１２】
本発明は、上記課題を解決するための手段として、特許請求の範囲の各請求項に記載された内燃機関の制御パラメータの適合値を求める方法を提供する。
【００１３】
請求項１に記載の発明は、機関負荷と機関回転数によって定まる内燃機関の各運転状態において同内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であって、上記制御パラメータの値を様々な値に設定して同制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて同特性パラメータと上記制御パラメータとの関係を全運転状態に渡って求め、少なくとも同関係に基づいて各運転状態における上記制御パラメータの適合値の存在範囲を推定し、次いで同存在範囲内において上記制御パラメータの値を様々な値に設定して同制御パラメータの各値毎に上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて同特性パラメータと上記制御パラメータとの各運転状態における関係を求め、同関係に基づいて各運転状態における上記制御パラメータの適合値を求める、内燃機関の制御パラメータの適合値を求める方法を提供する。
【００１４】
請求項１に記載の発明のようにすれば、まず上記制御パラメータの適合値の存在範囲を推定し、その後、その存在範囲内において集中的に特性パラメータの計測を行うことになるので、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めることが可能となる。すなわち、制御パラメータの適合値の精度を効率的に向上することが出来る。
【００１５】
請求項２に記載の発明は、機関負荷と機関回転数によって定まる内燃機関の各運転状態において同内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であって、上記制御パラメータの値を様々な値に設定して同制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータ、並びに機関負荷と機関回転数をパラメータとして含むモデル式である全域モデルのモデル式を求め、同全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定することと、内燃機関の各運転状態毎に、内燃機関の運転状態がその運転状態になる場合において上記制御パラメータの値を上記適合値の存在範囲内で様々な値に設定した場合における上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータをパラメータとして含む一方、機関負荷と機関回転数をパラメータとして含まないモデル式である局所モデルのモデル式を求め、同局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求めることと、を含む内燃機関の制御パラメータの適合値を求める方法を提供する。
【００１６】
先にも述べたように上記全域モデルを用いて制御パラメータの適合値を求める場合には、計測工数が少ない反面、求められる適合値の精度が充分でない場合がある。他方、上記局所モデルを用いれば求められる適合値の精度を向上することができるが、その場合には計測工数が相当増大してしまう。
【００１７】
これに対し、請求項２に記載の発明では、まず上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定し、次いでその存在範囲内でのみ制御パラメータを様々な値に設定して特性パラメータの値の計測を行い局所モデルのモデル式を求めるようになっている。このようにすると局所モデルのモデル式を求めるための計測工数が大幅に低減できるので、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めることが可能となる。
【００１８】
請求項３に記載の発明では請求項２に記載の発明において、上記制御パラメータの適合値の存在範囲が、上記全域モデルのモデル式から得られる特性パラメータの値の他、上記全域モデルのモデル式を求める際に上記特性パラメータの計測を実施した上記制御パラメータの範囲と、上記全域モデルのモデル式を求める際のモデル化誤差と、上記制御パラメータを制御する精度と、上記制御パラメータを計測する精度と、上記特性パラメータを計測する精度とのうちの少なくとも一つに基づいて推定されるようになっている。
【００１９】
請求項３に記載の発明のようにすれば、上記制御パラメータの適合値の存在範囲（すなわち、局所モデルを作成する範囲）を適切に推定することができ、計測工数の不必要な増大を招くことなくより高精度な適合値を求めることができる。
【００２０】
請求項４に記載の発明は、機関負荷と機関回転数によって定まる内燃機関の各運転状態において同内燃機関の特性が要求条件を満たすようにする制御パラメータの値である制御パラメータの適合値を求める方法であり、上記制御パラメータの値を様々な値に設定して同制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて、特性パラメータの値を表すモデル式であって上記制御パラメータ、並びに機関負荷と機関回転数をパラメータとして含むモデル式である全域モデルのモデル式を求め、同全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求めると共に、内燃機関の運転状態のうちの予め定められた特定運転状態については、内燃機関の運転状態が上記特定運転状態になる場合において上記制御パラメータの値を様々な値に設定した場合における上記特性パラメータの値を計測し、それら特性パラメータの計測値に基づいて、特性パラメータの値を表すモデル式であって上記制御パラメータをパラメータとして含む一方、機関負荷と機関回転数をパラメータとして含まないモデル式である局所モデルのモデル式を求め、同局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を再度求めるようになっていて、上記特定運転状態以外の運転状態に関しては上記全域モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値を最終的な適合値とする一方、上記特定運転状態に関しては上記局所モデルのモデル式から得られる特性パラメータの値に基づいて再度求められた上記制御パラメータの適合値を最終的な適合値とする、制御パラメータの適合値を求める方法であって、上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲は、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて決定される、内燃機関の制御パラメータの適合値を求める方法を提供する。
【００２１】
上述したように上記全域モデルを用いて制御パラメータの適合値を求める場合には、計測工数が少ない反面、求められる適合値の精度が充分でない場合がある。他方、上記局所モデルを用いれば求められる適合値の精度を向上することができるが、その場合には計測工数が相当増大してしまう。
【００２２】
これに対し請求項４に記載の発明では、局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲は、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて決定されるようになっている。このようにすると局所モデルのモデル式を求めるための計測工数が大幅に低減できるので、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めることが可能となる。
【００２３】
また、上記制御パラメータの適合値の精度の問題は、実際の運転で使用される頻度の高い運転状態に対して設定された制御パラメータの適合値の精度が充分でない場合に特に影響が大きく、良好な制御を実現できなくなる恐れがある。
【００２４】
これに対し、請求項４に記載の発明では、予め定められた特定運転状態以外の運転状態に関しては上記全域モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値を最終的な適合値とする一方、上記特定運転状態に関しては局所モデルのモデル式を求め、同局所モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値を最終的な適合値とするようになっている。つまり、請求項４に記載の発明では、例えば実際の運転で使用される頻度の高い運転状態等を上記特定運転状態とすることで、予め定めた内燃機関の制御上重要な運転状態に対してのみ局所モデルを構築し、同局所モデルを用いてより高精度の制御パラメータの適合値が求められるようになっている。そしてこのようにすることで、計測工数の増大をより一層抑制しつつ、制御パラメータの適合値の精度の向上をより効果的に行うことができる。
【００２５】
請求項５に記載の発明では請求項４に記載の発明において、上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲が、上記全域モデルのモデル式から得られる特性パラメータの値の他、上記全域モデルのモデル式を求める際に上記特性パラメータの計測を実施した上記制御パラメータの範囲と、上記全域モデルのモデル式を求める際のモデル化誤差と、上記制御パラメータを制御する精度と、上記制御パラメータを計測する精度と、上記特性パラメータを計測する精度とのうちの少なくとも一つに基づいて決定されるようになっている。
【００２６】
請求項５に記載の発明のようにすれば、上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲（すなわち、局所モデルを作成する範囲）を適切に決定することができ、計測工数の不必要な増大を招くことなくより高精度な適合値を求めることができる。
【発明の効果】
【００２７】
各請求項に記載の発明は、計測工数の増大を抑えつつ内燃機関の制御パラメータの適合値を精度よく求めることができるという共通の効果を奏する。
【発明を実施するための最良の形態】
【００２８】
以下、図面を参照して本発明について詳細に説明する。図１は後述する適合作業の対象となる内燃機関、並びに当該適合作業に用いられる各計測装置及び制御装置を示している。
【００２９】
図１を参照すると１は機関本体、２はシリンダブロック、３はシリンダブロック２内で往復動するピストン、４はシリンダブロック２上に固定されたシリンダヘッド、５はピストン３とシリンダヘッド４との間に形成された燃焼室、６は吸気弁、７は吸気ポート、８は排気弁、９は排気ポートをそれぞれ示す。図１に示したようにシリンダヘッド４の内壁面の中央部には点火プラグ１０が配置され、シリンダヘッド４内壁面周辺部には燃料噴射弁１１が配置される。またピストン３の頂面上には燃料噴射弁１１の下方から点火プラグ１０の下方まで延びるキャビティ１２が形成されている。
【００３０】
各気筒の吸気ポート７はそれぞれ対応する吸気枝管１３を介してサージタンク１４に連結され、サージタンク１４は吸気管１５に連結される。吸気管１５内にはステップモータ１７によって駆動されるスロットル弁１８が配置される。一方、各気筒の排気ポート９は排気マニホルド１９に連結される。また、吸気弁６には吸気弁６の開閉弁特性、すなわち位相角及び作用角を変更するための可変動弁機構２０が取付けられている。
【００３１】
一般に、図１に示したような内燃機関の制御は、内燃機関の運転中に変化するトルク、排気エミッション及び燃費等についての要求条件を満たすように、すなわち実際のトルク、排気エミッション及び燃費等が目標トルク、目標排気エミッション及び目標燃費等となるように、内燃機関の運転状態に影響を与える制御可能なパラメータ（すなわち、制御パラメータ）の値を変化させることによって行われる。
【００３２】
このような制御パラメータには、そのときの内燃機関に対する要求に応じて、例えば、機関負荷及び機関回転数により定まる運転状態毎に最適な値が存在する。例えば、点火プラグ１０による点火時期については、内燃機関のトルク、燃費や失火等を考慮すると、一般に、トルクが最も大きくなるような最小進角時期、いわゆるＭＢＴ（Minimum Advance for Best Torque）付近で点火を行うのが好ましい。このＭＢＴは、全ての運転状態に対して同じではなく、例えば機関回転数が異なると、ＭＢＴも異なる時期となる。また、一方で、内燃機関の排気浄化のために内燃機関の排気通路に設けられた排気浄化触媒（図示せず）を高温にする必要があるような場合には、機関本体１から排出される排気ガスの温度（以下、「排気温度」と称す）を高めるために上記ＭＢＴよりも或る程度進角側の時期に点火を行うのが好ましい。
【００３３】
このような内燃機関に対する要求に応じた運転状態毎の各制御パラメータの最適な値（すなわち、適合値）は、数値計算のみから算出することは困難であるため、通常、内燃機関の形式毎に適合作業によって求められる。ここで、適合作業とは、特定の制御パラメータを様々な値に設定し、各制御パラメータの値毎に特性パラメータ（制御パラメータの値を変更することによりその値が変わり得るパラメータであって内燃機関の特性を表すパラメータ）を計測し、これら特性パラメータの計測値から各運転状態に対する制御パラメータの最適な値（すなわち、適合値）を求める作業を意味する。
【００３４】
図１には、適合作業の対象となる内燃機関に加えて、この内燃機関の特性パラメータを計測する各計測装置（各種センサ）、並びに制御装置４０が示されている。図示したように、適合作業の対象となる内燃機関に対しては、スロットル弁１８の開度を計測するためのスロットル開度センサ３１がスロットル弁１８に取付けられ、また、吸気管１５内を流れる空気の流量を計測するエアフロメータ３２がスロットル弁１８上流側の吸気管１５内に取付けられる。さらに、機関本体１から排出された排気ガスの温度を計測する排気温度センサ３３及び機関本体１から排出された排気ガスの空燃比を計測する空燃比センサ３４が排気ポート又は排気マニホルド１９に取付けられる。さらに、機関本体１のクランクシャフト（図示せず）には内燃機関による駆動力であるトルクを検出するためのトルクセンサ（図示せず）が取付けられる。これらセンサ３１〜３４は、制御装置４０に接続されている。制御装置４０ではこれらセンサ３１〜３４によって計測された各特性パラメータの値が記憶されると共に演算処理され、各制御パラメータの適合値が求められるようになっている。なお、本実施形態ではこの制御装置４０において後に詳述するような方法で各制御パラメータの適合値が求められるようになっている。
【００３５】
なお、上述したスロットル弁駆動用のステップモータ１７、燃料噴射弁１１及び点火プラグ１０も制御装置４０に接続されており、これらステップモータ１７等は制御装置４０によって駆動、制御される。すなわち、制御装置４０によって制御パラメータの値が変更される。
【００３６】
ところで、一般に適合作業においては、定常運転時における特性パラメータの値を計測する必要があることから、制御パラメータの値を変化させて設定される各計測点毎に運転状態が安定するまで待ってから、すなわち例えばトルク、機関回転数がほぼ一定の値に収束するまで待ってから計測が行われる。このため、各計測点において特性パラメータの計測値を得るまでに長い時間を要する。また、適合精度を高いものとするためには、多くの計測点において計測が必要であり、場合によっては計測点数が数千〜数十万点にも及ぶ。このため、適合作業全体の計測工数は膨大なものとなってしまう。
【００３７】
そこで、近年、各計測点間の間隔を広げ、すなわち各計測点間の制御パラメータの値の差を大きくして計測点数を低減すると共に、特性パラメータの計測値に基づいてモデル式を求め、斯かるモデル式に基づいて計測点間の各制御パラメータの値に対する特性パラメータの値を推定することが提案されている。このようにすると、計測工数を低減して各制御パラメータの適合値を求めることが可能となる。
【００３８】
そして上記特性パラメータの値を表すモデル式として各制御パラメータに加えて機関負荷及び機関回転数もパラメータとして含むモデル式を用い、全運転状態（すなわち、実際に運転が行われると想定される全ての運転状態）を一つの全域モデルで表現するようにする場合がある。このようにする場合には比較的少ない計測工数で全運転状態における各制御パラメータの適合値を求めることが可能となる。
【００３９】
ところが、この場合に上記各制御パラメータの適合値として求められる値は、実際にはその精度が充分でない場合がある。この原因の一つは上記全域モデルを構築するためには全運転状態に渡って計測を実施する必要があるため、依然として計測期間が相当長期間となり、その結果気象条件や内燃機関の状態等の計測条件に変化が生じてモデル精度が低下することにある。
【００４０】
一方、上記特性パラメータの値を表すモデル式として各運転状態毎に機関負荷と機関回転数をパラメータとして含まないモデル式（すなわち、機関負荷と機関回転数を固定した場合のモデル式）を求めるようにし、このモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値を求める方法もある。すなわち、この場合には機関負荷と機関回転数をパラメータとして含まない局所モデル（すなわち、機関負荷と機関回転数を固定した場合のモデル）が構築されて、各制御パラメータの適合値が求められる。
【００４１】
この局所モデルを用いた場合には、一般に上記全域モデルを用いた場合に比べて上記制御パラメータの適合値を精度よく求めることが出来るのであるが、各運転状態毎にモデル式を求めるのに充分な数の特性パラメータの計測値が必要となるため、計測工数は上記全域モデルを用いた場合に比べて相当に大きくなってしまう。
【００４２】
そこで本実施形態では、以上のような点を踏まえ、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めるために、以下で説明するような方法で各制御パラメータの適合値を求めるようにしている。
【００４３】
すなわち、本実施形態では、まず上記制御パラメータの値が全運転状態に渡って様々な値に設定されて同制御パラメータの各値毎に内燃機関の特性を表す特性パラメータの値が計測され、それら特性パラメータの計測値に基づいて同特性パラメータと上記制御パラメータとの関係が全運転状態に渡って求められ、少なくともこの関係に基づいて各運転状態における上記制御パラメータの適合値の存在範囲が推定されるようになっている。
【００４４】
そして次に、同存在範囲内において集中的に上記制御パラメータの値が様々な値に設定されて同制御パラメータの各値毎に上記特性パラメータの値が計測され、それら特性パラメータの計測値に基づいて同特性パラメータと上記制御パラメータとの各運転状態におけるより詳細な関係が求められ、同関係に基づいて各運転状態における上記制御パラメータの適合値が求められる。
【００４５】
このようにすると、まず上記制御パラメータの適合値の存在範囲を推定し、その後、その存在範囲内において集中的に特性パラメータの計測を行うことになるので、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めることが可能となる。すなわち、制御パラメータの適合値の精度を効率的に向上することが出来る。
【００４６】
より詳細には本実施形態では、まず上記制御パラメータの値を全運転状態に渡って様々な値に設定して計測した特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータ、並びに機関負荷と機関回転数をパラメータとして含むモデル式である全域モデルのモデル式が求められる。そして、同全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲が推定される。
【００４７】
次いで、内燃機関の各運転状態毎に、内燃機関の運転状態が同一の運転状態になる場合において上記制御パラメータの値を上記適合値の存在範囲内で様々な値に設定した場合における上記特性パラメータの値が計測され、それら特性パラメータの計測値に基づいて特性パラメータの値を表すモデル式であって上記制御パラメータをパラメータとして含む一方、機関負荷と機関回転数をパラメータとして含まないモデル式である局所モデルのモデル式が求められる。そして同局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値が求められる。
【００４８】
このように、本実施形態ではまず上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定し、次いでその存在範囲内でのみ制御パラメータを様々な値に設定して特性パラメータの値の計測を行い局所モデルのモデル式を求めるようになっている。このようにすると局所モデルのモデル式を求めるための計測工数が大幅に低減できるので、計測工数の増大を抑えつつ制御パラメータの適合値を精度よく求めることが可能となる。なおここで、全域モデルのモデル式で用いられる機関負荷を表すパラメータとしては発生トルクや筒内空気充填率（筒内充填空気質量を、１気圧、２５℃の状態において排気量に相当する容積を占める空気の質量で除算した値）等を用いることができる。
【００４９】
また、他の実施形態においては上記全運転状態について上記局所モデルのモデル式が求められ、同局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値が求められるようになっていてもよいが、本実施形態では、上記全運転状態のうちの予め定められた特定運転状態についてのみ局所モデルのモデル式が求められ、同局所モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値が求められるようになっている。
【００５０】
すなわち、本実施形態では、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定する際に、上記全域モデルのモデル式から得られる特性パラメータの値のみに基づいて上記制御パラメータの適合値Ｇｐが求められるようになっている。そして、上記特定運転状態に関しては更に上述したように上記局所モデルのモデル式を求めて、それから得られる特性パラメータの値に基づいて上記制御パラメータの適合値Ｌｐを再度求めるようになっている。
【００５１】
そして本実施形態では、上記特定運転状態以外の運転状態に関しては上記全域モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値Ｇを最終的な適合値とする一方、上記特定運転状態に関しては上記局所モデルのモデル式から得られる特性パラメータの値に基づいて再度求められた上記制御パラメータの適合値Ｌｐを最終的な適合値とするようにしている。
【００５２】
このようにすると、計測工数の増大をより一層抑制しつつ、制御パラメータの適合値の精度の向上をより効果的に行うことができる。すなわち、上述したような制御パラメータの適合値の精度の問題は、実際の運転で使用される頻度の高い運転状態に対して設定された制御パラメータの適合値の精度が充分でない場合に特に影響が大きく、良好な制御を実現できなくなる恐れがある。
【００５３】
これに対し、本実施形態では上述したように、予め定められた特定運転状態以外の運転状態に関しては上記全域モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値Ｇｐを最終的な適合値とする一方、上記特定運転状態に関しては局所モデルのモデル式を求め、同局所モデルのモデル式から得られる特性パラメータの値に基づいて求められた上記制御パラメータの適合値Ｌｐを最終的な適合値とするようになっている。つまり、本実施形態では、例えば実際の運転で使用される頻度の高い運転状態等を上記特定運転状態とすることで、予め定めた内燃機関の制御上重要な運転状態に対してのみ局所モデルを構築し、同局所モデルを用いてより高精度の制御パラメータの適合値が求められるようになっている。したがって、本実施形態のようにすることで、計測工数の増大をより一層抑制しつつ、制御パラメータの適合値の精度の向上をより効果的に行うことができるのである。
【００５４】
次に、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記制御パラメータの適合値の存在範囲を推定する方法、すなわち上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記局所モデルのモデル式を求める際に特性パラメータの計測を行う時の制御パラメータの値を設定する範囲（以下、「局所モデル作成範囲」という）を決定する方法ついて説明する。
【００５５】
図２は本実施形態において上記局所モデル作成範囲を決定するために実施される方法について説明するためのフローチャートである。また、図３は図２のフローチャートで示される制御フローを実施した場合の一例について説明するための図である。以下、これらの図を参照しつつ説明する。なお、ここでは説明及び理解を容易にするために一つの特性パラメータＹに対し一つの制御パラメータＸを用いて説明する。すなわち図３の横軸は制御パラメータＸを示し、縦軸は特性パラメータＹを示している。また図３の曲線Ｙａは、上記全域モデルのモデル式から得られる特性パラメータＹの値を示したものである。
【００５６】
図２の制御フローがスタートすると、まずステップ１０１において上記全域モデルのモデル式から得られる特性パラメータＹの値（すなわち上記曲線Ｙａ）に基づいて特性パラメータＹの値の最良値Ｙｐが求められる。またこれと同時に、上記特性パラメータＹの値が最良値Ｙｐとなる時の上記制御パラメータＸの値、すなわち上記全域モデルに基づいた上記制御パラメータＸの適合値Ｇｐも求められる。なお、事前に上記特性パラメータＹの最良値Ｙｐ及び制御パラメータＸの適合値Ｇｐが求められている場合には、ステップ１０１においてこれらの値の読込みが行われる。
【００５７】
ステップ１０１において上記特性パラメータＹの最良値Ｙｐ及び制御パラメータＸの適合値Ｇｐが求められると（もしくは読込まれると）ステップ１０３に進む。ステップ１０３では後のステップにおいて上記局所モデル作成範囲を決定するために使用される上記特性パラメータＹの最良値Ｙｐからの落ち幅ΔＹが設定される。この落ち幅ΔＹは、予め定めた方法によって設定されるようになっていても良いし、予め定めた一定値とされても良い。すなわち、例えば上記落ち幅ΔＹを上記全域モデルのモデル化誤差を表す値（例えば、残差二乗和平均（ＲＭＳＥ））としても良い。この場合、上記全域モデルのモデル化の精度が低い程、上記落ち幅ΔＹは大きくなる。
【００５８】
ステップ１０３において上記落ち幅ΔＹが設定されると、ステップ１０５に進み、下記（１）式が設定される。
Ｙ（ｋ＋１）＝Ｙ（ｋ）−ΔＹ          …  （１）
そして続くステップ１０７において、ｋ＝０、Ｙ（０）＝Ｙｐという初期値が与えられ、ステップ１０９に進んで上記（１）式によりＹ（ｋ＋１）が算出される。
【００５９】
ステップ１０９において上記Ｙ（ｋ＋１）が算出されるとステップ１１１に進む。ステップ１１１では、上記Ｙ（ｋ＋１）に対応する制御パラメータＸの範囲Ｒ（ｋ＋１）が求められる。この制御パラメータＸの範囲Ｒ（ｋ＋１）は、図３からも明らかなように、上記全域モデルのモデル式に基づいた場合に、上記特性パラメータＹの値がステップ１０９で求めたＹ（ｋ＋１）となる上記制御パラメータＸの値のうち、上記適合値Ｇｐよりも大きい側及び小さい側のそれぞれで最も上記適合値Ｇｐに近い値の間で表される範囲である。
【００６０】
ステップ１１１において上記制御パラメータＸの範囲Ｒ（ｋ＋１）が求められると、ステップ１１３に進む。ステップ１１３ではステップ１１１で求められた上記制御パラメータＸの範囲Ｒ（ｋ＋１）の大きさを制御パラメータＸの設定可能な全範囲Ｒｓ（例えば、全域モデルのモデル式を求める際に特性パラメータの計測を行う時に上記制御パラメータＸの値を設定した範囲）の大きさで除算した値（＝｜Ｒ（ｋ＋１）｜／｜Ｒｓ｜）が予め定めた基準比率Ｃｒ以上であるか否かが判定される。ここで、この基準比率Ｃｒは例えば上記制御パラメータＸの制御精度に基づいて設定される。すなわち、上記制御パラメータＸの制御精度が低い程、上記基準比率Ｃｒとして大きな値が設定される。
【００６１】
また、ここでは説明上、制御パラメータが一つであるとしているためステップ１１３で判定される比率は長さの比率となっているが、当然のことながら、制御パラメータが二つの場合には面積の比率となり、三つの場合には体積の比率となる。
【００６２】
ステップ１１３において、上記制御パラメータＸの範囲Ｒ（ｋ＋１）の大きさを制御パラメータＸの設定可能な全範囲Ｒｓの大きさで除算した値が上記基準比率Ｃｒ以上であると判定された場合にはステップ１１５に進み、その時に求められている上記制御パラメータＸの範囲Ｒ（ｋ＋１）が上記局所モデル作成範囲Ｒｍとされて制御フローが終了する。一方、ステップ１１３において上記制御パラメータＸの範囲Ｒ（ｋ＋１）の大きさを制御パラメータＸの設定可能な全範囲Ｒｓの大きさで除算した値が上記基準比率Ｃｒ未満であると判定された場合にはステップ１１７に進み、ｋにｋ＋１が代入されて再度ステップ１０９からの制御が繰り返される。図３に示した例では、ステップ１０９からの制御が３度繰り返されてＹ（３）及びＲ（３）まで求められ、Ｒ（３）が上記局所モデル作成範囲Ｒｍとされている。
【００６３】
以上の説明からも明らかなように、この方法が実施された場合には、上記制御パラメータＸの適合値Ｇｐの付近において上記制御パラメータＸの値の変化に対して上記特性パラメータＹの値の変化が緩やかである場合には、その変化が急峻である場合に比べて上記局所モデル作成範囲が広くなるようになっている。つまり、この方法が実施された場合には、制御パラメータの値の変化に対する特性パラメータの値の変化の傾向に応じた上記局所モデル作成範囲の設定が可能となる。そしてその結果、計測工数の不必要な増大を招くことなく、より高精度な適合値が求められるように局所モデルを作成する範囲を適切に設定することができる。
【００６４】
次に、上記全域モデルのモデル式から得られる特性パラメータの値に基づいて上記局所モデル作成範囲を決定する別の方法について説明する。図４はこの方法で上記局所モデル作成範囲を決定する場合について説明するためのフローチャートであり、図５は図４のフローチャートで示される制御フローを実施した場合の一例について説明するための図である。以下、これら図４及び図５を参照しつつ説明する。なお、ここでも説明及び理解を容易にするために一つの特性パラメータＹに対し一つの制御パラメータＸを用いて説明する。すなわち図５の横軸及び縦軸は、図３の場合と同様、それぞれ制御パラメータＸ及び特性パラメータＹを示している。図５中、曲線Ｙａは上記全域モデルのモデル式から得られる特性パラメータＹの値を示したものである。
【００６５】
また、図５中の一点鎖線Ｙｒｊは上記全域モデルのモデル式から得られる特性パラメータＹの値の信頼性区間の上限値を示したものであり、一点鎖線Ｙｒｋは上記全域モデルのモデル式から得られる特性パラメータＹの値の信頼性区間の下限値を示したものである。より詳細には、図５に示された例では上記二つの一点鎖線ＹｒｊとＹｒｋとの間は上記全域モデルのモデル式から得られる特性パラメータＹの値の９５％信頼性区間を示している。すなわち、特性パラメータＹの真値は９５％の確率で上記二つの一点鎖線Ｙｒｊ及びＹｒｋの間に存在する。
【００６６】
図４の制御フローがスタートすると、まずステップ２０１において上記全域モデルのモデル式から得られる特性パラメータＹの値（すなわち上記曲線Ｙａ）に基づいて上記制御パラメータＸの適合値Ｇｐが求められる。この適合値Ｇｐは、すなわち上記全域モデルのモデル式から得られる特性パラメータＹの値（すなわち上記曲線Ｙａ）に基づいて特性パラメータＹの値が最良値Ｙｐとなる時の上記制御パラメータＸの値である。ここでは、上記制御パラメータＸの適合値Ｇｐが求められるときに上記特性パラメータＹの最良値Ｙｐも同時に求められるようになっている。なお、事前に上記制御パラメータＸの適合値Ｇｐ及び特性パラメータＹの最良値Ｙｐが求められている場合には、ステップ２０１においてこれらの値の読込みが行われる。
【００６７】
ステップ２０１において上記制御パラメータＸの適合値Ｇｐ及び特性パラメータＹの最良値Ｙｐが求められると（もしくは読込まれると）ステップ２０３に進む。ステップ２０３では最小確保範囲Ｒａが求められる。この最小確保範囲Ｒａの大きさは上記制御パラメータＸの制御上もしくは計測上の最小分解能（すなわち、それより小さい大きさでは充分な精度をもって制御したり計測することができない大きさ）に基づいて求められる。また、この最小確保範囲Ｒａは図５に示されているように、上記制御パラメータＸの適合値Ｇｐを中心として設定される。
【００６８】
ステップ２０３において最小確保範囲Ｒａが求められるとステップ２０５に進む。ステップ２０５では信頼性区間考慮範囲Ｒｂが求められる。この信頼性区間考慮範囲Ｒｂは、上記全域モデルのモデル式から得られる特性パラメータＹの値（すなわち上記曲線Ｙａ）のモデル化誤差を考慮して定められる。図５に示された例では、上記全域モデルのモデル式から得られる特性パラメータＹの値の９５％信頼性区間の上限値を示した一点鎖線Ｙｒｊにおいて上記特性パラメータＹの値が上記最良値Ｙｐとなる時の二つの制御パラメータＸの値Ｒｂ１、Ｒｂ２の間の範囲を上記信頼性区間考慮範囲Ｒｂとしている。
【００６９】
ステップ２０５において信頼性区間考慮範囲Ｒｂが求められるとステップ２０７に進む。ステップ２０７では計測誤差考慮範囲Ｒｃが求められる。この計測誤差考慮範囲Ｒｃは、上記全域モデルのモデル式を求める際に行った特性パラメータＹの計測における計測誤差（すなわち、計測機器の誤差や計測条件の変化に伴う計測値のバラツキ等）を考慮して定められる。図５に示された例では、上記全域モデルのモデル式から得られる特性パラメータＹの値（すなわち上記曲線Ｙａ）に基づいて求められた上記特性パラメータＹの最良値Ｙｐの計測誤差としてΔｅを想定し、上記曲線Ｙａにおいて上記特性パラメータＹの値が上記最良値Ｙｐから上記計測誤差Δｅ分だけ減算した値Ｙｅとなる時の二つの制御パラメータＸの値Ｒｃ１、Ｒｃ２の間の範囲を上記計測誤差考慮範囲Ｒｃとしている。
【００７０】
ステップ２０７において計測誤差考慮範囲Ｒｃが求められるとステップ２０９に進む。ステップ２０９ではステップ２０５で求められた信頼性区間考慮範囲Ｒｂとステップ２０７で求められた計測誤差考慮範囲Ｒｃとの論理和となる範囲（すなわち、信頼性区間考慮範囲Ｒｂと計測誤差考慮範囲Ｒｃの少なくとも一方に含まれる制御パラメータＸの範囲）Ｒｂｃが求められる。図５に示された例では、上記信頼性区間考慮範囲Ｒｂが上記計測誤差考慮範囲Ｒｃ内に含まれているので、上記計測誤差考慮範囲Ｒｃがそのまま上記範囲Ｒｂｃとなる。
【００７１】
ステップ２０９において上記範囲Ｒｂｃが求められるとステップ２１１に進む。ステップ２１１では全域モデル計画範囲Ｒｄが読込まれる。この全域モデル計画範囲Ｒｄは、上記全域モデルのモデル式を求めるために行われた特性パラメータＹの計測実験の計画範囲である。図５に示された例では、上記全域モデル計画範囲Ｒｄは、後述する制御パラメータ制御可能範囲Ｒｇから範囲Ｄｘを除いた部分となっている。この例において上記範囲Ｄｘは、上記制御パラメータＸの値がその範囲Ｄｘ内にある場合には、例えば機関の運転上の理由等で実際に実験を実施して計測値を得るのが困難な範囲である。そのため、このような実験を実施できない範囲は、たとえ制御パラメータＸが制御可能な範囲であっても上記全域モデル計画範囲Ｒｄからは除外されている。なお、この範囲Ｄｘにおいては実際に計測した値が得られないため、上記全域モデルのモデル式に基づいた外挿によって上記特性パラメータＹの値が求められることになる。そのため、この範囲Ｄｘ内においては上述した信頼性区間が広がっている。
【００７２】
ステップ２１１において上記全域モデル計画範囲Ｒｄが読込まれるとステップ２１３に進む。ステップ２１３では全域モデル制約考慮範囲Ｒｅが求められる。この全域モデル制約考慮範囲Ｒｅは、上記全域モデルのモデル式を求めるために行われた特性パラメータＹの計測実験において制約となった様々な条件について考慮した上記制御パラメータＸの範囲である。図５に示された例では、上記全域モデル制約考慮範囲Ｒｅは、後述する制御パラメータ制御可能範囲Ｒｇから範囲Ｅｘ１及びＥｘ２を除いた部分となっている。この例において上記範囲Ｅｘ１及びＥｘ２はそれぞれが異なる制約を表す範囲である（すなわち、例えば上記制御パラメータＸの値が上記範囲Ｅｘ１内にある場合及び上記範囲Ｅｘ２内にある場合には何れも精度の良い計測値を得ることが困難となるが、その原因は異なっている）。つまり、上記全域モデル制約考慮範囲Ｒｅは、例えば上記制御パラメータＸの制御可能な範囲Ｒｇからこのような制約を表す範囲Ｅｘ１、Ｅｘ２を総て除外することにより求めることができる。また、この全域モデル制約考慮範囲Ｒｅと上記全域モデル計画範囲Ｒｄとの論理積となる範囲（すなわち、上記全域モデル制約考慮範囲Ｒｅと上記全域モデル計画範囲Ｒｄの両方に含まれる制御パラメータＸの範囲）が、上記全域モデルのモデル式を求める際に上記特性パラメータＹの計測を実際に実施した範囲である。
【００７３】
ステップ２１３において上記全域モデル制約考慮範囲Ｒｅが求められるとステップ２１５に進む。ステップ２１５では制御パラメータＸの制御可能な範囲である制御パラメータ制御可能範囲Ｒｇが読込まれる。ステップ２１５に続いてステップ２１７に進むと、そこで上記範囲Ｒｂｃ、Ｒｄ、Ｒｅ及びＲｇの論理積となる範囲（すなわち、上記範囲Ｒｂｃ、Ｒｄ、Ｒｅ及びＲｇの総てに含まれる制御パラメータＸの範囲）Ｒｈが求められる。図５に示された例では、上記範囲Ｒｂｃ（＝上記範囲Ｒｃ）がその他の範囲Ｒｄ、Ｒｅ及びＲｇ内に含まれているので、上記範囲Ｒｂｃ（＝上記範囲Ｒｃ）がそのまま上記範囲Ｒｈとなる。
【００７４】
ステップ２１７において上記範囲Ｒｈが求められるとステップ２１９に進む。ステップ２１９ではステップ２０３で求められた最小確保範囲Ｒａとステップ２１７で求められた上記範囲Ｒｈとの論理和となる範囲（すなわち、最小確保範囲Ｒａと上記範囲Ｒｈの少なくとも一方に含まれる制御パラメータＸの範囲）Ｒｍが求められる。この範囲Ｒｍが上記局所モデル作成範囲となる。ステップ２１９において上記範囲Ｒｍが求められると本制御フローが終了する。図５に示された例では、上記最小確保範囲Ｒａが上記範囲Ｒｈ（＝上記範囲Ｒｃ）内に含まれているので、上記範囲Ｒｈ（＝上記範囲Ｒｃ）がそのまま局所モデル作成範囲Ｒｍとなる。
【００７５】
以上の説明からも明らかなように、この方法では上記全域モデルのモデル式を求める際の情報（すなわち、例えば上記全域モデルのモデル化誤差や上記全域モデルのモデル式を求める際に行った特性パラメータＹの計測における計測誤差等）を利用して、局所モデルを用いて求められる上記制御パラメータの適合値Ｌｐの存在範囲である上記局所モデル作成範囲が決定されるようになっている。そして、このようにすることによっても計測工数の不必要な増大を招くことなく、より高精度な適合値が求められるように局所モデルを作成する範囲を適切に設定することができる。
【００７６】
なお、以上の説明からも明らかであると思われるが、上述した各種範囲のうち上記最小確保範囲Ｒａ、信頼性区間考慮範囲Ｒｂ及び計測誤差考慮範囲Ｒｃは、局所モデル作成範囲Ｒｍが含むべき範囲である。一方、上述した各種範囲のうち全域モデル計画範囲Ｒｄ、全域モデル制約考慮範囲Ｒｅ及び制御パラメータ制御可能範囲Ｒｇは、局所モデル作成範囲Ｒｍが含まれるべき範囲である。そして、上述した方法では上記範囲Ｒａ、Ｒｂ及びＲｃの総てを含み、且つ上記範囲Ｒｄ、Ｒｅ及びＲｇの総てに含まれる最小範囲を求め、上記局所モデル作成範囲Ｒｍとしていると言える。したがって、このような考え方にしたがえば図４に示された制御フローの各ステップの順序を入れ替えても同様に上記局所モデル作成範囲Ｒｍが求められることは容易に理解できる。
【００７７】
また、上記全域モデル制約考慮範囲Ｒｅ等については、その範囲設定において信頼性区間や誤差を考慮するようにしても良い。上記全域モデル制約考慮範囲Ｒｅの場合、制約のある範囲（例えば、上記範囲Ｅｘ１やＥｘ２）の信頼性区間や誤差についての判断によって上記全域モデル制約考慮範囲Ｒｅが広く設定されたり、狭く設定されたりすることになり、その結果的にとして上記局所モデル作成範囲Ｒｍが広く設定されたり、狭く設定されたりすることになる。局所モデルの精度の向上を重視する場合には、上記局所モデル作成範囲Ｒｍが広く設定されるようにすると有効であり、絞込みの度合いを高めて計測工数を低減することを重視する場合には上記局所モデル作成範囲Ｒｍが狭く設定されるようにすると有効である。
【図面の簡単な説明】
【００７８】
【図１】図１は、適合作業の対象となる内燃機関、並びに適合作業に用いられる各計測装置及び制御装置を示す図である。
【図２】図２は、局所モデル作成範囲を決定するために実施される方法について説明するためのフローチャートである。
【図３】図３は、図２のフローチャートで示される制御フローを実施した場合の一例について説明するための図である。
【図４】図４は、局所モデル作成範囲を決定するために実施される別の方法について説明するためのフローチャートである。
【図５】図５は、図４のフローチャートで示される制御フローを実施した場合の一例について説明するための図である。
【符号の説明】
【００７９】
１    機関本体
６    吸気弁
８    排気弁
１０    点火プラグ
１１    燃料噴射弁
１８    スロットル弁
３１    スロットル開度センサ
３２    エアフロメータ
３３    排気温度センサ
３４    空燃比センサ
４０    制御装置
【図１】
0004419918.tif 000002
【図２】
0004419918.tif 000003
【図３】
0004419918.tif 000004
【図４】
0004419918.tif 000005
【図５】
0004419918.tif 000006