JPS61190143A - 内燃機関制御装置 - Google Patents
内燃機関制御装置Info
- Publication number
- JPS61190143A JPS61190143A JP60031265A JP3126585A JPS61190143A JP S61190143 A JPS61190143 A JP S61190143A JP 60031265 A JP60031265 A JP 60031265A JP 3126585 A JP3126585 A JP 3126585A JP S61190143 A JPS61190143 A JP S61190143A
- Authority
- JP
- Japan
- Prior art keywords
- internal combustion
- combustion engine
- amount
- intake air
- target
- 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.)
- Granted
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D31/00—Use of speed-sensing governors to control combustion engines, not otherwise provided for
- F02D31/001—Electric control of rotation speed
- F02D31/002—Electric control of rotation speed controlling air supply
- F02D31/003—Electric control of rotation speed controlling air supply for idle speed control
- F02D31/004—Electric control of rotation speed controlling air supply for idle speed control by controlling a throttle stop
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D43/00—Conjoint electrical control of two or more functions, e.g. ignition, fuel-air mixture, recirculation, supercharging or exhaust-gas treatment
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1413—Controller structures or design
- F02D2041/1415—Controller structures or design using a state feedback or a state space representation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1413—Controller structures or design
- F02D2041/1415—Controller structures or design using a state feedback or a state space representation
- F02D2041/1416—Observer
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1413—Controller structures or design
- F02D2041/1426—Controller structures or design taking into account control stability
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1433—Introducing closed-loop corrections characterised by the control or regulation method using a model or simulation of the system
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
[産業上の利用分野]
本発明は内燃機関を制御する装置であって、特に、内燃
機関の運転に関する系の動的なモデル(ダイナミックモ
デル)に基づいて内燃機関のアイドル回転数を最小の燃
料噴射量で目標回転数に維持するようにした内燃機関制
御装置に関する。
機関の運転に関する系の動的なモデル(ダイナミックモ
デル)に基づいて内燃機関のアイドル回転数を最小の燃
料噴射量で目標回転数に維持するようにした内燃機関制
御装置に関する。
[従来の技術]
従来から、内燃機関のアイドル回転数を制御する装置と
して、例えば特開昭54−76723号公報に示される
如く、機関回転数を回転数検出手段により検出し、この
検出値と冷却水温等に応じて設定された目標アイドル回
転数との偏差をなくすよう、言い換えれば、機関回転数
を目標アイドル回転数に保持するよう、スロットルバル
ブのバイパス通路を通過する空気流量を制御するフィー
ドバック制御によるアイドル回転数制御装置が提案され
ている。
して、例えば特開昭54−76723号公報に示される
如く、機関回転数を回転数検出手段により検出し、この
検出値と冷却水温等に応じて設定された目標アイドル回
転数との偏差をなくすよう、言い換えれば、機関回転数
を目標アイドル回転数に保持するよう、スロットルバル
ブのバイパス通路を通過する空気流量を制御するフィー
ドバック制御によるアイドル回転数制御装置が提案され
ている。
[発明が解決しようとする問題点]
しかしながら、上記のような制御装置によると、アイド
ル運転中にエアコンがオンまたはオフしたことによりそ
れぞれ目標アイドル回転数をそれまでの値から所定値だ
け増加もしくは減少させる場合に、所望の追従性を実現
することが難しく、燃費向上と出力確保といういわば相
反する問題を充分に解決するまでには至らない。
ル運転中にエアコンがオンまたはオフしたことによりそ
れぞれ目標アイドル回転数をそれまでの値から所定値だ
け増加もしくは減少させる場合に、所望の追従性を実現
することが難しく、燃費向上と出力確保といういわば相
反する問題を充分に解決するまでには至らない。
本発明は上記のような問題点を解決することを目的とし
、いわゆる現代制御理論を取り入れ、燃費向上、出力確
保を実現できる充分な追従性を得ると共に燃費が最小と
なる吸入空気量で機関回転数を目標アイドル回転数に維
持することを目的としている。
、いわゆる現代制御理論を取り入れ、燃費向上、出力確
保を実現できる充分な追従性を得ると共に燃費が最小と
なる吸入空気量で機関回転数を目標アイドル回転数に維
持することを目的としている。
[問題点を解決するための手段]
このような目的を達成するための本発明の構成は、第1
図に示すように、 内燃機関M1の吸気管に配設される燃料噴射弁M2、お
よび吸気通路のスロットルバルブM3の開度を調節する
スロットルアクチュエータM4と、 機関回転数検出手段M5、吸入空気流量検出手段M6お
よび機関制御のためのその他の検出手段M7と、 これら検出手段M5、M6、M7からの検出信号N、A
R,・・・に基づき上記燃料噴射弁M2へ燃料噴射弁駆
動信倍FRを、上記スロットルアクチュエータM4ヘス
ロットル間度を規定するアクチュエータ駆動信号θをそ
れぞれ出力する電子制御手段であって、目標アイドル回
転数N*を設定するアイドル回転数設定部M9を備える
ものと、 からなる内燃機関制御装置において、 上記電子制御手段M8を、 上記アイドル回転数設定部M9による目標アイドル回転
数N*の情報を受け、機関回転数Nをこの目標アイドル
回転数N*に維持できる最小の燃料噴射量に対応する吸
入空気11A R*を、目標回転数N本に対応して予め
定めた燃料噴射量・吸入空気量相関パターンに基づいて
定める目標吸入空気量設定部M10と、 この目標吸入空気量設定部M10による目標吸入空気M
AR”の情報、上記アイドル回転数設定部M9による目
標回転数N*の情報および内燃機関M1の制御入出力値
FR,θ、N、ARをそれぞれ受け、内燃機関M1の運
転に関する系の動的なモデルに従って予め定められた機
閣内部状態推定のための演算式および最適フィードバッ
クゲインに基づいて、燃料噴射量およびスロットル開度
のそれぞれについてのフィードバック量を決定し、これ
ら決定されたフィードバック量のそれぞれに基づく燃料
噴射弁駆動信号FRおよびアクチュエータ駆動信号θを
出力する付加積分型レギュレータM11と、から構成し
た ことを要旨とする。
図に示すように、 内燃機関M1の吸気管に配設される燃料噴射弁M2、お
よび吸気通路のスロットルバルブM3の開度を調節する
スロットルアクチュエータM4と、 機関回転数検出手段M5、吸入空気流量検出手段M6お
よび機関制御のためのその他の検出手段M7と、 これら検出手段M5、M6、M7からの検出信号N、A
R,・・・に基づき上記燃料噴射弁M2へ燃料噴射弁駆
動信倍FRを、上記スロットルアクチュエータM4ヘス
ロットル間度を規定するアクチュエータ駆動信号θをそ
れぞれ出力する電子制御手段であって、目標アイドル回
転数N*を設定するアイドル回転数設定部M9を備える
ものと、 からなる内燃機関制御装置において、 上記電子制御手段M8を、 上記アイドル回転数設定部M9による目標アイドル回転
数N*の情報を受け、機関回転数Nをこの目標アイドル
回転数N*に維持できる最小の燃料噴射量に対応する吸
入空気11A R*を、目標回転数N本に対応して予め
定めた燃料噴射量・吸入空気量相関パターンに基づいて
定める目標吸入空気量設定部M10と、 この目標吸入空気量設定部M10による目標吸入空気M
AR”の情報、上記アイドル回転数設定部M9による目
標回転数N*の情報および内燃機関M1の制御入出力値
FR,θ、N、ARをそれぞれ受け、内燃機関M1の運
転に関する系の動的なモデルに従って予め定められた機
閣内部状態推定のための演算式および最適フィードバッ
クゲインに基づいて、燃料噴射量およびスロットル開度
のそれぞれについてのフィードバック量を決定し、これ
ら決定されたフィードバック量のそれぞれに基づく燃料
噴射弁駆動信号FRおよびアクチュエータ駆動信号θを
出力する付加積分型レギュレータM11と、から構成し
た ことを要旨とする。
ここでアイドル回転数設定部M9は、その他検出手段M
7に含まれるアイドルスイッチ、冷却水温センサ、車速
センサ、エアコンスイッチなどからの信号に基づき、例
えば、アイドルスイッチがオンかつ車速がOkm/ h
に近い所定値以下のとき冷却水漏とエアコンスイッチの
オン・オフに応じた目標アイドル回転数を設定する。
7に含まれるアイドルスイッチ、冷却水温センサ、車速
センサ、エアコンスイッチなどからの信号に基づき、例
えば、アイドルスイッチがオンかつ車速がOkm/ h
に近い所定値以下のとき冷却水漏とエアコンスイッチの
オン・オフに応じた目標アイドル回転数を設定する。
目標吸入空気量設定部MIOは次のように所定の目標ア
イドル回転数に対して燃料供給量を最小にする吸入空気
量を設定する。
イドル回転数に対して燃料供給量を最小にする吸入空気
量を設定する。
第2図は内燃vlA関M10目標アイドル回転数N−7
一 本を一定とした時の吸入空気量ARと燃料供給量FRと
の関係を示す定回転数線図である。今、吸入空気量がA
b、燃料供給量がFbの点すで内燃機関が目標アイドル
回転数N*で運転されていたとすると、ここから吸入空
気量をΔAO増加させた点(Aa、、Fa)において燃
料供給1tFaが最小になることがわかる。目標吸入空
気量設定部M10は吸入空気量の目標値AR*に関して
、所定の目標アイドル回転数N*に対して上述のように
燃料供給11FRを最小とする値を与えるよう構成され
るが、一般的には後述する電子iI+制御手段M8の一
部として、マイクロコンピュータ等による処理の一つと
して実現される。
一 本を一定とした時の吸入空気量ARと燃料供給量FRと
の関係を示す定回転数線図である。今、吸入空気量がA
b、燃料供給量がFbの点すで内燃機関が目標アイドル
回転数N*で運転されていたとすると、ここから吸入空
気量をΔAO増加させた点(Aa、、Fa)において燃
料供給1tFaが最小になることがわかる。目標吸入空
気量設定部M10は吸入空気量の目標値AR*に関して
、所定の目標アイドル回転数N*に対して上述のように
燃料供給11FRを最小とする値を与えるよう構成され
るが、一般的には後述する電子iI+制御手段M8の一
部として、マイクロコンピュータ等による処理の一つと
して実現される。
電子制御手段M8は、通常マイクロプロセッサを用いR
OM、RAM等の周辺素子や入出力回路と共に構成され
た電子回路として実現され、予め記憶された処理手順に
従って、目標吸入空気量設定部MIO、アイドル回転数
設定部M9によって設定された目標値と内燃機関M1の
運転状態の諸量とを知って、運転状態が目標に近づくよ
う、予め内燃tII閏M1の運転に関する系の動的なモ
デルに従って定められた最適フィードバックゲインから
定まるフィードバック量により制御するよう構成されて
いる。即ち、制御手段M8は目標吸入空気量設定部M1
0によって設定された目標吸入空気量と内燃機1$1M
1の運転状態の諸量とから最適なフィードバック量を定
める付加積分型最適レギュレータとして構成されている
。
OM、RAM等の周辺素子や入出力回路と共に構成され
た電子回路として実現され、予め記憶された処理手順に
従って、目標吸入空気量設定部MIO、アイドル回転数
設定部M9によって設定された目標値と内燃機関M1の
運転状態の諸量とを知って、運転状態が目標に近づくよ
う、予め内燃tII閏M1の運転に関する系の動的なモ
デルに従って定められた最適フィードバックゲインから
定まるフィードバック量により制御するよう構成されて
いる。即ち、制御手段M8は目標吸入空気量設定部M1
0によって設定された目標吸入空気量と内燃機1$1M
1の運転状態の諸量とから最適なフィードバック量を定
める付加積分型最適レギュレータとして構成されている
。
こうした付加積分型最適レギュレータの構成の手法は、
例えば古田勝久著「線形システム制御理論」 (昭和5
1年NIB晃堂等に詳しいが、ここで実際の構成の手法
について一通の見通しを与えることにする。尚、以下の
説明においてF、X、A。
例えば古田勝久著「線形システム制御理論」 (昭和5
1年NIB晃堂等に詳しいが、ここで実際の構成の手法
について一通の見通しを与えることにする。尚、以下の
説明においてF、X、A。
B、C,v、u、L、Go G、IR,T、Pはベクト
ル!(行列)を示し、ATの如き添字Tは行列の転置を
、A−1の如き添字−1は逆行列を、更にXの如き添字
△はそれが推定値であることを、6の如き記号−は制御
対象の系から変換等により生成された別の系、ここでは
状態観測器(以下、オブザーバと呼ぶ)で扱われている
量であることを、y*の如き記号*は目標値であること
を、各々示している。
ル!(行列)を示し、ATの如き添字Tは行列の転置を
、A−1の如き添字−1は逆行列を、更にXの如き添字
△はそれが推定値であることを、6の如き記号−は制御
対象の系から変換等により生成された別の系、ここでは
状態観測器(以下、オブザーバと呼ぶ)で扱われている
量であることを、y*の如き記号*は目標値であること
を、各々示している。
制御対象、ここでは内燃機関M1の制御において、この
制御対象の動的な振舞は、 X (k ) −=A・X (k−1> 十B−u (
k−1)・・・(1) V (k )−C−X (k ) ・・・
(2)として記述されることが現代制御理論より知られ
ている。ここで式(1)は状態方程式1式(2)は出力
方程式と呼ばれ、X(k)は内燃機11M1の内部状態
を表わす状態変数量であり、(A (k )は内燃i関
M1の制御入力つまり運転条件を示づ金量からなるベク
トル、s (k )は内燃機関M1の1IillIll
出力つまり運転状態を示すliI量からなるベクトルで
ある。又、式(1)、(2>は離散系で記述されており
、添字には現時点であることを、k−1は1回前のサン
プリング時点であることを、各々示している。
制御対象の動的な振舞は、 X (k ) −=A・X (k−1> 十B−u (
k−1)・・・(1) V (k )−C−X (k ) ・・・
(2)として記述されることが現代制御理論より知られ
ている。ここで式(1)は状態方程式1式(2)は出力
方程式と呼ばれ、X(k)は内燃機11M1の内部状態
を表わす状態変数量であり、(A (k )は内燃i関
M1の制御入力つまり運転条件を示づ金量からなるベク
トル、s (k )は内燃機関M1の1IillIll
出力つまり運転状態を示すliI量からなるベクトルで
ある。又、式(1)、(2>は離散系で記述されており
、添字には現時点であることを、k−1は1回前のサン
プリング時点であることを、各々示している。
内燃機II]Mlの内部状態を示す状態変数量X(k)
は、その制御系における未来への影響を予−1〇 − 測するために必要十分な系の履歴に関する情報を示して
いる。従って、内燃機関M1の運転に関する系の動的な
モデルが明らかになり、式(1)。
は、その制御系における未来への影響を予−1〇 − 測するために必要十分な系の履歴に関する情報を示して
いる。従って、内燃機関M1の運転に関する系の動的な
モデルが明らかになり、式(1)。
(2)のベクトルA、IB、Cを定めることができれば
、状態変数量X(k)を用いて内燃機関M1の運転を最
適に制御できることになる。尚、サーボ系においては系
を拡大する必要が生じるが、これについては後述する。
、状態変数量X(k)を用いて内燃機関M1の運転を最
適に制御できることになる。尚、サーボ系においては系
を拡大する必要が生じるが、これについては後述する。
ところが、内燃機関M1のように複雑な対象については
その動的なモデルを理論的に正確に求めることは困難で
あり、何らかの形で実験的に定めることが必要となる。
その動的なモデルを理論的に正確に求めることは困難で
あり、何らかの形で実験的に定めることが必要となる。
これが所謂システム同定と呼ばれるモデル横築の手法で
あって、内燃機関M1が所定の運転状態で運転されてい
る場合、その状態の近傍では線形の近似が成立つとして
、式(1)、(2)の状態方程式に則ってモデルを構築
するのである。従って、内燃機関M1のようにその運転
に関する動的なモデルが非線形のような場合にも、定常
的な複数の運転状態に分離することによって線形な近似
を行なうことができ、個々の動的なモデルを定めること
によって広範な運転域まで拡張することができるのであ
る。
あって、内燃機関M1が所定の運転状態で運転されてい
る場合、その状態の近傍では線形の近似が成立つとして
、式(1)、(2)の状態方程式に則ってモデルを構築
するのである。従って、内燃機関M1のようにその運転
に関する動的なモデルが非線形のような場合にも、定常
的な複数の運転状態に分離することによって線形な近似
を行なうことができ、個々の動的なモデルを定めること
によって広範な運転域まで拡張することができるのであ
る。
ここで、制御対象が比較的容易に物理的なモデルとして
構築できるのものであれば周波数応答法やスペクトル解
析法といった手法によりシステム同定を行なって、動的
な系のモデル(ここではベクトルA、 E3. C)を
定めることができるが、内燃機関M1のような多元系の
制御対象は、ある程度近似のよい物理モデルをつくるこ
とも困難であり、この場合には最小2乗法や補助変数法
あるいはオンライン同定法などにより動的なモデルの構
築を行なう。
構築できるのものであれば周波数応答法やスペクトル解
析法といった手法によりシステム同定を行なって、動的
な系のモデル(ここではベクトルA、 E3. C)を
定めることができるが、内燃機関M1のような多元系の
制御対象は、ある程度近似のよい物理モデルをつくるこ
とも困難であり、この場合には最小2乗法や補助変数法
あるいはオンライン同定法などにより動的なモデルの構
築を行なう。
動的なモデルが定まれば、状態変数量X(k)と運転状
態の諸量s (k )及びその目標値V*(k)からフ
ィードバック量が定まり運転条件の諸量%A(k )が
理論的に最適に定められる。通常内燃機関Ml等では内
燃機関M1の運転に直接関与する諸量として、例えば実
際に吸入されている空気量や燃焼の動的挙動、あるいは
燃焼に関与している混合気中の燃料量や内燃機関の出力
トルクといった量を状態変数量X(k)として扱えばよ
いのであるが、これらの諸量の大部分は直接観測するこ
とが極めて困難である。そこで、こうし゛た場合には、
電子制御手段M8内に状態観測器(オブザーバ)と呼ば
れる手段を構成し、内燃機関M1の運転条件の諸量と運
転状態の諸量とを用いて、内燃機関M1の状態変数量X
(k)を推定することができる。これが所謂、現代制御
理論におけるオブザーバであり、種々のオブザーバとそ
の設計法が知られている。これらは、例えば古田勝久他
著「メカニカルシステム制御」(昭和59年)オーム社
等に詳解されており、適用される制御対象、ここでは内
燃機関M1とその運転状態制御装置との態様に合わせて
最小次元オブザーバや有限整定オブザーバとして設計す
ればよい。
態の諸量s (k )及びその目標値V*(k)からフ
ィードバック量が定まり運転条件の諸量%A(k )が
理論的に最適に定められる。通常内燃機関Ml等では内
燃機関M1の運転に直接関与する諸量として、例えば実
際に吸入されている空気量や燃焼の動的挙動、あるいは
燃焼に関与している混合気中の燃料量や内燃機関の出力
トルクといった量を状態変数量X(k)として扱えばよ
いのであるが、これらの諸量の大部分は直接観測するこ
とが極めて困難である。そこで、こうし゛た場合には、
電子制御手段M8内に状態観測器(オブザーバ)と呼ば
れる手段を構成し、内燃機関M1の運転条件の諸量と運
転状態の諸量とを用いて、内燃機関M1の状態変数量X
(k)を推定することができる。これが所謂、現代制御
理論におけるオブザーバであり、種々のオブザーバとそ
の設計法が知られている。これらは、例えば古田勝久他
著「メカニカルシステム制御」(昭和59年)オーム社
等に詳解されており、適用される制御対象、ここでは内
燃機関M1とその運転状態制御装置との態様に合わせて
最小次元オブザーバや有限整定オブザーバとして設計す
ればよい。
電子制御手段M8は、観測された状態変数量または上記
のオブザーバによって推定された状態変数量X(k)の
他に、目標吸入空気量設定部M10によって設定された
目標吸入空気量、アイドル回転数設定部M9によって設
定された目標アイドル回転数と、実際の吸入空気量、機
関回転数との各々の偏差を累積した累積値を用いて拡大
された系において、両者と、予め定められた最適フィー
ドバックゲインとから最適なフィードバック酸を定め燃
料噴割弁M2およびスロットルアクチュエータM4を制
御する。累積値は運転状態の目標値が内燃機11M1へ
の要求量によって変化することから必要となる社である
。一般にサーボ系の制御においては目標値と実際の制御
値との定常偏差を消去するような制御が必要となり、こ
れは伝達関数において1/S(1次の積分)を含む必要
があるとされる。また、既述したようなシステム同定に
より系の伝達関数を定め、これから状態方程式をたてて
いるような場合には、対ノイズ安定性の上からもこうし
た積分量を含むことが望ましい。
のオブザーバによって推定された状態変数量X(k)の
他に、目標吸入空気量設定部M10によって設定された
目標吸入空気量、アイドル回転数設定部M9によって設
定された目標アイドル回転数と、実際の吸入空気量、機
関回転数との各々の偏差を累積した累積値を用いて拡大
された系において、両者と、予め定められた最適フィー
ドバックゲインとから最適なフィードバック酸を定め燃
料噴割弁M2およびスロットルアクチュエータM4を制
御する。累積値は運転状態の目標値が内燃機11M1へ
の要求量によって変化することから必要となる社である
。一般にサーボ系の制御においては目標値と実際の制御
値との定常偏差を消去するような制御が必要となり、こ
れは伝達関数において1/S(1次の積分)を含む必要
があるとされる。また、既述したようなシステム同定に
より系の伝達関数を定め、これから状態方程式をたてて
いるような場合には、対ノイズ安定性の上からもこうし
た積分量を含むことが望ましい。
本発明においては立−1、即ち一次型の積分を考慮すれ
ばよい。従って、上述の状態変数量X(k)にこの累積
値を加えて系を拡大し、両者と予め定められた最適なフ
ィードバックゲイン「とにより帰還量を定めれば、付加
積分型最適レギュレータとして、制御対象への制御入力
値、即ち内燃ta関M1の運転条件の諸量が定まる。
ばよい。従って、上述の状態変数量X(k)にこの累積
値を加えて系を拡大し、両者と予め定められた最適なフ
ィードバックゲイン「とにより帰還量を定めれば、付加
積分型最適レギュレータとして、制御対象への制御入力
値、即ち内燃ta関M1の運転条件の諸量が定まる。
次に、最適フィードバックゲインについて説明する。上
記の如く積分量を付加した最適レギュレータでは、評価
関数Jを最小とするような制御入力(ここでは内燃機関
M1の運転条件の諸量)の求め方が明らかにされており
、最適フィードバックゲインもリカッチ方程式の解と状
態方程式(1)、出力方程式(2)のA、IB、Cマト
リックス及び評価関数に用いられる重みパラメータ行列
とから求められることがわかっている(前掲書他)。
記の如く積分量を付加した最適レギュレータでは、評価
関数Jを最小とするような制御入力(ここでは内燃機関
M1の運転条件の諸量)の求め方が明らかにされており
、最適フィードバックゲインもリカッチ方程式の解と状
態方程式(1)、出力方程式(2)のA、IB、Cマト
リックス及び評価関数に用いられる重みパラメータ行列
とから求められることがわかっている(前掲書他)。
ここで重みパラメータは当初任意に与えられるものであ
って、評価関数Jが内燃機関M1の運転条件諸量の挙動
を制約する重みを変更するものである。重みパラメータ
を任意に与えて大型コンピュータによるシミュレーショ
ンを行ない、得られた運転状11諸量の挙動から重みパ
ラメータを所定量変更してシミュレーションを繰返し、
最適な値を決定しておくことができる。その結果最適フ
ィードバックゲインFも定められる。
って、評価関数Jが内燃機関M1の運転条件諸量の挙動
を制約する重みを変更するものである。重みパラメータ
を任意に与えて大型コンピュータによるシミュレーショ
ンを行ない、得られた運転状11諸量の挙動から重みパ
ラメータを所定量変更してシミュレーションを繰返し、
最適な値を決定しておくことができる。その結果最適フ
ィードバックゲインFも定められる。
従って、本発明の内燃機関制御装置の電子制御手段M8
は、予めシステム同定等により決定された内燃機関M1
の動的モデルを用いて付加積分型最適レギュレータM1
1として構成され、その内部におけるオブザーバのパラ
メータや最適フィードバックゲイン「などは、全て、予
め内燃機関M1を用いたシミュレーションにより決定さ
れているのである。
は、予めシステム同定等により決定された内燃機関M1
の動的モデルを用いて付加積分型最適レギュレータM1
1として構成され、その内部におけるオブザーバのパラ
メータや最適フィードバックゲイン「などは、全て、予
め内燃機関M1を用いたシミュレーションにより決定さ
れているのである。
尚、以上の説明において状態変数量X(k)は内燃機関
M1の内部状態を表わす量として説明したが、これは実
際の物理量に対応した変数量である必要はなく、内燃機
関M1の状態を表わす適当な次数のベクトル団として設
計することができる。
M1の内部状態を表わす量として説明したが、これは実
際の物理量に対応した変数量である必要はなく、内燃機
関M1の状態を表わす適当な次数のベクトル団として設
計することができる。
[作用]
上記構成を有する本発明の内燃機関制御装置は、それぞ
れ目標吸入空気量設定部M10により内燃11関M1へ
の目標吸入空気量を、アイドル回転数設定部M9により
目標アイドル回転数を算出し、付加積分型最適レギュレ
ータM11により、内燃機関M1の運転状態の諸量が上
記の目標値となるよう最適フィードバック量を求めて燃
料噴射弁M2およびスロットルアクチュエータM4の制
御を行なうよう働く。しかも目標吸入空気量設定部M1
0は内燃機関M1のアイドル回転数を一定にした状態で
燃料消費量が最小となるよう目標吸入空気量を算出する
よう働くので、結果的に本発明の内燃機関制御装置は、
内燃機関M1を目標アイドル回転数でしかも燃料消費量
を最小とするような運転状態に最適に制御する。
れ目標吸入空気量設定部M10により内燃11関M1へ
の目標吸入空気量を、アイドル回転数設定部M9により
目標アイドル回転数を算出し、付加積分型最適レギュレ
ータM11により、内燃機関M1の運転状態の諸量が上
記の目標値となるよう最適フィードバック量を求めて燃
料噴射弁M2およびスロットルアクチュエータM4の制
御を行なうよう働く。しかも目標吸入空気量設定部M1
0は内燃機関M1のアイドル回転数を一定にした状態で
燃料消費量が最小となるよう目標吸入空気量を算出する
よう働くので、結果的に本発明の内燃機関制御装置は、
内燃機関M1を目標アイドル回転数でしかも燃料消費量
を最小とするような運転状態に最適に制御する。
[実施例]
次に、本発明の実施例を図面に基づいて詳細に説明する
。第3図は本発明の実施例における内燃機関とその周辺
装置を表わす概略構成図、第4図は内燃機関の運転状態
を制御する系の制御モデルを示す制−系統図、第5図、
第6図はそれぞれシステム同定の説明に用いるブロック
線図、第7図は電子制御回路において実行される制御の
一例を示すフローチャート、第8図は燃料消費を最小と
する吸入空気量を求める1Ii11御の一例を示すフロ
ーチャート、第9図は本実施例の効果を説明するグラフ
である。以下、この順に説明する。
。第3図は本発明の実施例における内燃機関とその周辺
装置を表わす概略構成図、第4図は内燃機関の運転状態
を制御する系の制御モデルを示す制−系統図、第5図、
第6図はそれぞれシステム同定の説明に用いるブロック
線図、第7図は電子制御回路において実行される制御の
一例を示すフローチャート、第8図は燃料消費を最小と
する吸入空気量を求める1Ii11御の一例を示すフロ
ーチャート、第9図は本実施例の効果を説明するグラフ
である。以下、この順に説明する。
第3図では4気筒4サイクルの内燃i関1のひとつの気
筒を中心に示しているが、吸気系2には上流から図示し
ないエアクリーナ、吸入空気流量AR@測定するエア7
0メータ3、吸気温Thaを検出する吸気温センサ5、
吸入空気量を制御するスロットルバルブ7、サージタン
ク9、電磁式の燃料噴射弁11等が備えられている。ま
た内燃機関1の排気は排気管14により図示しない排気
浄化装置、消音器等を介して外部へ排出される。
筒を中心に示しているが、吸気系2には上流から図示し
ないエアクリーナ、吸入空気流量AR@測定するエア7
0メータ3、吸気温Thaを検出する吸気温センサ5、
吸入空気量を制御するスロットルバルブ7、サージタン
ク9、電磁式の燃料噴射弁11等が備えられている。ま
た内燃機関1の排気は排気管14により図示しない排気
浄化装置、消音器等を介して外部へ排出される。
燃焼室(シリンダ)はピストン15.吸気弁17゜排気
弁199点火プラグ21等から構成されているが、これ
らの作動は周知のものなので説明は省略する。
弁199点火プラグ21等から構成されているが、これ
らの作動は周知のものなので説明は省略する。
内燃機関1にはこの他に冷却水の温度Thwを検出する
冷却水温センサ29やディストリビュータ25内に備え
られ、内燃機関1の回転数Nに応じた周波数のパルス信
号を出力する回転数センサ31や、内燃機関1の1回転
(クランク角の720°)に1発のパルス信号を出力す
る気筒判別セフサ33等が備えられている。又、スロッ
トルバルブ7は直流モータを動ノj源とするスロットル
アクチュエータ35によってその開度θを制御されてい
る。尚、第3図中37はアクセル38の踏込量Aceを
検出するアクセル開度センサである。
冷却水温センサ29やディストリビュータ25内に備え
られ、内燃機関1の回転数Nに応じた周波数のパルス信
号を出力する回転数センサ31や、内燃機関1の1回転
(クランク角の720°)に1発のパルス信号を出力す
る気筒判別セフサ33等が備えられている。又、スロッ
トルバルブ7は直流モータを動ノj源とするスロットル
アクチュエータ35によってその開度θを制御されてい
る。尚、第3図中37はアクセル38の踏込量Aceを
検出するアクセル開度センサである。
以上の構成を有する内燃機関1とその周辺装置において
、その燃料噴射1iFRやスロットルバルブ開度θ等は
電子制御回路40によって制御されている。電子制御回
路40はキースイッチ41を介してバッテリ43より電
力の供給をうけて作動しているが、周知のマイクロプロ
セッサ(MPU)44、ROM45.RAM46.バッ
クアップRAM47.入力ボート49.出力ボート51
等から構成され、上記各素子・ボートは相互にバス53
により接続されている。
、その燃料噴射1iFRやスロットルバルブ開度θ等は
電子制御回路40によって制御されている。電子制御回
路40はキースイッチ41を介してバッテリ43より電
力の供給をうけて作動しているが、周知のマイクロプロ
セッサ(MPU)44、ROM45.RAM46.バッ
クアップRAM47.入力ボート49.出力ボート51
等から構成され、上記各素子・ボートは相互にバス53
により接続されている。
電子制御回路40の入力ポート49は、内燃機関1の運
転状態を示す信号を各センサより入力する。具体的には
、吸入空気流量ARをエアフロメータ3より、吸気IT
haを吸気温センサ5より、機関アイドル状態であるか
ないかを判断するための信号をアイドルスイッチ55、
車速センサ56などにより冷却水温Thwを冷却水温セ
ンサ29より、内燃機関1の回転数Nを回転数センサ3
1より、気筒判別信号を気筒判別センサ33より、エア
コンがオンされているかどうかを判断するための信号を
エアコンスイッチ57より各々入力する入力部で構成さ
れている。
転状態を示す信号を各センサより入力する。具体的には
、吸入空気流量ARをエアフロメータ3より、吸気IT
haを吸気温センサ5より、機関アイドル状態であるか
ないかを判断するための信号をアイドルスイッチ55、
車速センサ56などにより冷却水温Thwを冷却水温セ
ンサ29より、内燃機関1の回転数Nを回転数センサ3
1より、気筒判別信号を気筒判別センサ33より、エア
コンがオンされているかどうかを判断するための信号を
エアコンスイッチ57より各々入力する入力部で構成さ
れている。
一方、出力ポート51は、アクチュエータ35を介して
スロットルバルブ7の開度θを、燃料噴射弁11を開・
閉弁して燃料噴射量FRを、その他、イグナイタ211
を介して点火時期などを、個々に制御する信号を出力す
る。これら電子制御回路40のMPU44による処理に
ついては、点火時期制御などのための処理を除いて、後
に第6図。
スロットルバルブ7の開度θを、燃料噴射弁11を開・
閉弁して燃料噴射量FRを、その他、イグナイタ211
を介して点火時期などを、個々に制御する信号を出力す
る。これら電子制御回路40のMPU44による処理に
ついては、点火時期制御などのための処理を除いて、後
に第6図。
第7図のフローチャー1〜に拠って詳述する。
次に第4図の機能ブロック図に拠って、電子制御回路4
0による処理について説明し、特にシステム同定による
状態方程式(1)、出力方程式(2)等におけるベクト
ルA、B、Cの求め方やこれに基づくオブザーバの求め
方、フィードパツクゲイン「の求め方について、実際に
即して説明する。尚、第4図は機能をブロックとして表
わした図であって、ハード構成を示すものではない。
0による処理について説明し、特にシステム同定による
状態方程式(1)、出力方程式(2)等におけるベクト
ルA、B、Cの求め方やこれに基づくオブザーバの求め
方、フィードパツクゲイン「の求め方について、実際に
即して説明する。尚、第4図は機能をブロックとして表
わした図であって、ハード構成を示すものではない。
又、第4図に示す機能は、実際には第6図のフローチャ
ートに示した一連のプログラムの実行による離散系とし
て実現されるものである。
ートに示した一連のプログラムの実行による離散系とし
て実現されるものである。
第4図に示すように、内燃Il!閏のアイドル回転数N
*は冷却水m T h W 、エアコンスイッチ57の
オン・オフに応じてアイドル回転数設定部P1により設
定される。一方、目標吸入空気量AR*は目標アイドル
回転数N ” +実際の検出された吸入空気量AR,回
転数N、内燃機関1に噴射された燃料噴射量FRから、
後に第8図の説明において詳述する手法により目標吸入
空気量設定部P2により燃料消費量を最小にする値とし
て定められる。積分器P3は目標アイドル回転数N*と
実際の回転数Nとの偏差SNを累積して累積値ZN(k
)を求め、積分器P4は目標吸入空気量AR*と実際の
吸入空気量ARとの偏差SARを累積して累積値ZAR
(k )を求めるものである。
*は冷却水m T h W 、エアコンスイッチ57の
オン・オフに応じてアイドル回転数設定部P1により設
定される。一方、目標吸入空気量AR*は目標アイドル
回転数N ” +実際の検出された吸入空気量AR,回
転数N、内燃機関1に噴射された燃料噴射量FRから、
後に第8図の説明において詳述する手法により目標吸入
空気量設定部P2により燃料消費量を最小にする値とし
て定められる。積分器P3は目標アイドル回転数N*と
実際の回転数Nとの偏差SNを累積して累積値ZN(k
)を求め、積分器P4は目標吸入空気量AR*と実際の
吸入空気量ARとの偏差SARを累積して累積値ZAR
(k )を求めるものである。
P5は、吸入空気量AR,回転数Nについて、定常的な
運転状態での8値(ARa、Na)からの摂動分を抽出
する摂動分抽山部を示している。
運転状態での8値(ARa、Na)からの摂動分を抽出
する摂動分抽山部を示している。
これは、既述したように、非線形なシステムに対して線
形の近似を行なう為に、内燃機関1の運転状態を、個々
の定常的な運転状態の連続とみなし、個々の定常運転近
傍について線形な動的モデルを構築して内燃機関1の運
転に関する広範囲にわたる動的なモデルを構築したこと
によっている。従って、内燃機関1の運転状態の諸1(
AR,N>を、一旦、最も近い定常な運転状態からの摂
動分δAR(=AR−ARa)、δN (=N−Na)
として扱うのである。前記の積分器P3.P4とオブザ
ーバP6とフィードバック量決定部P7とによって求め
られる内燃機関1への制御入力、即ちスロットル開度θ
、燃料噴射IFRに関する操作量も、摂動分δθ、δF
Rとして扱われている。
形の近似を行なう為に、内燃機関1の運転状態を、個々
の定常的な運転状態の連続とみなし、個々の定常運転近
傍について線形な動的モデルを構築して内燃機関1の運
転に関する広範囲にわたる動的なモデルを構築したこと
によっている。従って、内燃機関1の運転状態の諸1(
AR,N>を、一旦、最も近い定常な運転状態からの摂
動分δAR(=AR−ARa)、δN (=N−Na)
として扱うのである。前記の積分器P3.P4とオブザ
ーバP6とフィードバック量決定部P7とによって求め
られる内燃機関1への制御入力、即ちスロットル開度θ
、燃料噴射IFRに関する操作量も、摂動分δθ、δF
Rとして扱われている。
オブザーバP6は、運転状態の摂動分δAR。
δNと制御入力つまり運転条件の摂動分δθ、δFRと
から内燃機関1の□内部状態を表現する状態変数量X(
k)を推定して状態変数推定間貸(k)を求めるもので
あり、この状態変数推定量X(k)と上述の累積値ZT
(k >、 ZAR(k )とに、フィードバック量
決定部P7において、最適フィードバックゲイン「を乗
算し、操作量(δθ、δFR)を求めるのである。この
操作量の組(δθ。
から内燃機関1の□内部状態を表現する状態変数量X(
k)を推定して状態変数推定間貸(k)を求めるもので
あり、この状態変数推定量X(k)と上述の累積値ZT
(k >、 ZAR(k )とに、フィードバック量
決定部P7において、最適フィードバックゲイン「を乗
算し、操作量(δθ、δFR)を求めるのである。この
操作量の組(δθ。
δFR)は摂動分抽山部P5によって選ばれた定常的な
運転状態に対応した運転条件からの摂動分なので、これ
に基準設定値加算部P8によりこの定常的な運転条件に
対応した基準設定値θa、FRaを加えて、内燃機関1
に対する制御入力の諸量、θ、FRを定めるのである。
運転状態に対応した運転条件からの摂動分なので、これ
に基準設定値加算部P8によりこの定常的な運転条件に
対応した基準設定値θa、FRaを加えて、内燃機関1
に対する制御入力の諸量、θ、FRを定めるのである。
以上、簡単にこの制御系の構成について説明したが、こ
うした運転状態(AR,N)及び運転条件(θ、FR)
を実施例として取上げたのは、これらの諸量が内燃機関
1の出力に関与する基本的な量であることによっている
。従って、本実施例では内燃機関1を2人力2出力の多
元系としてとらえた。内燃機関1の出力に関与する量と
しては、この他にも、例えば点火時期や排ガス還流量な
ども考えられ、必要に応じてそれらを加味してシステム
のモデルをたてればよいが、ここでは、内燃機関1の動
的なモデルの横築には上記の2人力2出力のモデルを用
い、系の動的な振舞を変化させるものとして、他に内燃
機関1の冷却水温ThWや吸気温T h aを用いるが
、内燃機関1の冷却水温Thw等は内燃機関1の制御系
の構成を変えるものではなく、その動的な振舞の様子を
変化させるにすぎない。従って、内燃機関1の制御系に
ついてこの動的なモデルを構築する際、状態方程式(1
)、出力方程式(2)のベクトルA、B、Cが内燃機関
1の冷却水iThw等に応じて定められることになる。
うした運転状態(AR,N)及び運転条件(θ、FR)
を実施例として取上げたのは、これらの諸量が内燃機関
1の出力に関与する基本的な量であることによっている
。従って、本実施例では内燃機関1を2人力2出力の多
元系としてとらえた。内燃機関1の出力に関与する量と
しては、この他にも、例えば点火時期や排ガス還流量な
ども考えられ、必要に応じてそれらを加味してシステム
のモデルをたてればよいが、ここでは、内燃機関1の動
的なモデルの横築には上記の2人力2出力のモデルを用
い、系の動的な振舞を変化させるものとして、他に内燃
機関1の冷却水温ThWや吸気温T h aを用いるが
、内燃機関1の冷却水温Thw等は内燃機関1の制御系
の構成を変えるものではなく、その動的な振舞の様子を
変化させるにすぎない。従って、内燃機関1の制御系に
ついてこの動的なモデルを構築する際、状態方程式(1
)、出力方程式(2)のベクトルA、B、Cが内燃機関
1の冷却水iThw等に応じて定められることになる。
以上、内燃機関1のハード的な構成と内燃ll関1の出
力の制御を行なうものとして2人力2出力の系を取り上
げた場合の制御系の構成について説明した。そこで、次
に実際のシステム同定による動的モデルの構築、オブザ
ーバP6の設計、最適フィードバックゲイン「の与え方
について説明する。
力の制御を行なうものとして2人力2出力の系を取り上
げた場合の制御系の構成について説明した。そこで、次
に実際のシステム同定による動的モデルの構築、オブザ
ーバP6の設計、最適フィードバックゲイン「の与え方
について説明する。
まず内燃機関1の動的なモデルを構築する。第5図は2
人力2出力の系として定常運転されている内燃機関1の
系を伝達関数G1(Z)〜G4(2)により書き表わし
た図である。尚、2は入出力信号のサンプル値のZ変換
を示し、G1(z)〜G4(Z)は適当な次数をもつも
のとする。従って、全体の伝達関数行列G(Z)は、で
表わされる。
人力2出力の系として定常運転されている内燃機関1の
系を伝達関数G1(Z)〜G4(2)により書き表わし
た図である。尚、2は入出力信号のサンプル値のZ変換
を示し、G1(z)〜G4(Z)は適当な次数をもつも
のとする。従って、全体の伝達関数行列G(Z)は、で
表わされる。
本実施例の内燃機関1のように、その制御系が2人力2
出力の系であり、入出力の諸量に干渉が存在するような
場合には、物理的なモデルを定めることが極めて困難と
なる。このような場合には、システム同定と呼ばれる一
種のシミュレーションにより伝達関数を求めることがで
きる。
出力の系であり、入出力の諸量に干渉が存在するような
場合には、物理的なモデルを定めることが極めて困難と
なる。このような場合には、システム同定と呼ばれる一
種のシミュレーションにより伝達関数を求めることがで
きる。
システム同定の手法は、例えば相良節夫他著、「システ
ム同定」 (昭和56年)社団法人計測自動−1111
学会等に詳解されているが、ここでは最小2乗法により
同定する。
ム同定」 (昭和56年)社団法人計測自動−1111
学会等に詳解されているが、ここでは最小2乗法により
同定する。
内燃機関1を所定の運転状態で定常運転し、スロットル
開度の変化分δθを零として、燃料供給量の変化分δF
Rとしての適当な試験信号を加え、その時の入力δFR
と出力である回転数の変化分δNのデータをN回に亘っ
てサンプリングする。
開度の変化分δθを零として、燃料供給量の変化分δF
Rとしての適当な試験信号を加え、その時の入力δFR
と出力である回転数の変化分δNのデータをN回に亘っ
てサンプリングする。
これを入力のデータ系列(u(i))−(δFR1)、
出力のデータ系列(y(i))=(δN1)(但し、i
=1.2,3.・・・N)と表わす。この時、系は1
人力1出力とみなすことができ、系の伝達関数Gl(Z
)は、 G1 (z )=B (z−1)/A (z−1)・
(3)即ち、 Gl(Z) −(bo +b1 − z −1+・+bn z −’
) /(1+al ・2−壜+a2・2+・・−+a
n−z)・・・(4) で求められる。尚、ここで z−1は単位推移演算子で
あって、z−1−x (k )=x (k−1)を意味
している。
出力のデータ系列(y(i))=(δN1)(但し、i
=1.2,3.・・・N)と表わす。この時、系は1
人力1出力とみなすことができ、系の伝達関数Gl(Z
)は、 G1 (z )=B (z−1)/A (z−1)・
(3)即ち、 Gl(Z) −(bo +b1 − z −1+・+bn z −’
) /(1+al ・2−壜+a2・2+・・−+a
n−z)・・・(4) で求められる。尚、ここで z−1は単位推移演算子で
あって、z−1−x (k )=x (k−1)を意味
している。
入出力のデータ系列(u (i >)、 N/(i
>)から式(4)のパラメータa1〜an、bO〜bn
を定めれば系の伝達関数Gi(Z)が求められる。最小
2乗法によるシステム同定では、このパラメータa1〜
an、bO〜bnを JO−Σ[(V (k ) +a1 − y (k−1
) +−・・±an−y(k−n))−(bo−u(k
)+b1 ・u(k−1)+・・・ +bn −u (k−n ) )] ・・・(5) が最小となるよう定める。本実施例ではn=2として、
各パラメータを求めた。この場合、系のシグナルフロー
線図は第6図のようになり、状態変数量としr[xl(
k) x2(k)]”をとって、その状態・出力方程
式は、 ・・・ (6) と表わせられる。従って、1人力1出力の系とみなした
場合のシステムパラメータA、 E3. Cを各々AI
”、 IBl +、 C1′とすれば、となる。
>)から式(4)のパラメータa1〜an、bO〜bn
を定めれば系の伝達関数Gi(Z)が求められる。最小
2乗法によるシステム同定では、このパラメータa1〜
an、bO〜bnを JO−Σ[(V (k ) +a1 − y (k−1
) +−・・±an−y(k−n))−(bo−u(k
)+b1 ・u(k−1)+・・・ +bn −u (k−n ) )] ・・・(5) が最小となるよう定める。本実施例ではn=2として、
各パラメータを求めた。この場合、系のシグナルフロー
線図は第6図のようになり、状態変数量としr[xl(
k) x2(k)]”をとって、その状態・出力方程
式は、 ・・・ (6) と表わせられる。従って、1人力1出力の系とみなした
場合のシステムパラメータA、 E3. Cを各々AI
”、 IBl +、 C1′とすれば、となる。
本実施例ではG1(z)についてのパラメータとして、
[al ’a2]=[−1゜91 0.92
3][bObl b2] − [04,86X10 4.73 XIO]
を得た。同様の手法により、伝達関数G2(Z)ないし
G4(z)及び各々についてのシステムパラメータA2
′ないしA4−、[32′ないしB4−、C2′ない
し04 ′が求められる。そこでこれらのシステムパラ
メータから元の2人力2出力の多元系のシステムパラメ
ータ、即ち状態方程式(1)、出力方程式(2)のベク
トルA、 13.0を定めることができる。
3][bObl b2] − [04,86X10 4.73 XIO]
を得た。同様の手法により、伝達関数G2(Z)ないし
G4(z)及び各々についてのシステムパラメータA2
′ないしA4−、[32′ないしB4−、C2′ない
し04 ′が求められる。そこでこれらのシステムパラ
メータから元の2人力2出力の多元系のシステムパラメ
ータ、即ち状態方程式(1)、出力方程式(2)のベク
トルA、 13.0を定めることができる。
こうして本実施例の動的なモデルがシステム同定により
求められたが、この、動的なモデルは、内燃機関1が所
定の状態で運転されている時、この状態の近傍では線形
の近似が成立つという形で定められる。従って、定常的
な複数の運転状態に関して、上記の手法で伝達関数01
(Z)ないしG4(Z)が各々求められ、各々の状態方
程式(1)、出力方程式(2)におけるベクトルA。
求められたが、この、動的なモデルは、内燃機関1が所
定の状態で運転されている時、この状態の近傍では線形
の近似が成立つという形で定められる。従って、定常的
な複数の運転状態に関して、上記の手法で伝達関数01
(Z)ないしG4(Z)が各々求められ、各々の状態方
程式(1)、出力方程式(2)におけるベクトルA。
B、Cが求められ、その入出力の関係は摂動弁δの間に
成立することになる。
成立することになる。
次にオブザーバP6の設計方法について説明する。オブ
ザーバの設計にはゴピナスの設計法などがあって、古田
勝久・佐野昭共著「基礎システム理論」 (昭和53年
)コロナ社等々に詳しいが、本実施例では有限整定オブ
ザーバとして設計する。
ザーバの設計にはゴピナスの設計法などがあって、古田
勝久・佐野昭共著「基礎システム理論」 (昭和53年
)コロナ社等々に詳しいが、本実施例では有限整定オブ
ザーバとして設計する。
オブザーバP6は内燃機関1の運転状態の諸量の摂動弁
(δAR,δN)と運転条件の諸量の摂動弁(δθ、δ
FR)とから内燃機関1の内部の状態変数量X(k)を
推定するものであるが、オブザーバP6によって求めら
れた状態変数推定量X(k)を、内燃機関1の制御にお
いて、実際の状態変数量X(k)として扱うことができ
るという根拠は次の点にある。今、オブザーバP6の出
力X(k)を次式(9)のように構成したとする。
(δAR,δN)と運転条件の諸量の摂動弁(δθ、δ
FR)とから内燃機関1の内部の状態変数量X(k)を
推定するものであるが、オブザーバP6によって求めら
れた状態変数推定量X(k)を、内燃機関1の制御にお
いて、実際の状態変数量X(k)として扱うことができ
るという根拠は次の点にある。今、オブザーバP6の出
力X(k)を次式(9)のように構成したとする。
+旧−u (k−1)+L−v (k−1>・・・(9
) 式(9)においてLは任意に与えられる行列である。式
(1)、(2)、(9)より変形すると、[X (k
) −X (k ) ] −(]A−L−C[X (k−1>−X (k−1)
]・・・(10) を得る。従って(A−L−C)なる行列の固有値が単位
円内にある様に行列りを選択すればに→■でX(k)→
X(k)となり、制御対象の内部の状態変数量X(k)
を入力制御ベクトルu (k )、出力制御ベクトルs
(k )の過去からの系列U(*) 、 y (*)
を用いて正しく推定することができる。
) 式(9)においてLは任意に与えられる行列である。式
(1)、(2)、(9)より変形すると、[X (k
) −X (k ) ] −(]A−L−C[X (k−1>−X (k−1)
]・・・(10) を得る。従って(A−L−C)なる行列の固有値が単位
円内にある様に行列りを選択すればに→■でX(k)→
X(k)となり、制御対象の内部の状態変数量X(k)
を入力制御ベクトルu (k )、出力制御ベクトルs
(k )の過去からの系列U(*) 、 y (*)
を用いて正しく推定することができる。
今、最小2乗法によってシステム同定され定められた状
態方程式(1)、出力方程式(2)のベクトルA、B、
Cはこの系が可観測であることから、正則な行列■を用
いた新たな状態変数量ヌ(k)=T−1・X(k)を考
えて、次の可観測正準形に相似変換することができる。
態方程式(1)、出力方程式(2)のベクトルA、B、
Cはこの系が可観測であることから、正則な行列■を用
いた新たな状態変数量ヌ(k)=T−1・X(k)を考
えて、次の可観測正準形に相似変換することができる。
X(k)−Ao −X(k−1>+B13o−u(k
−1>・・・(11) y(k)=Co−X(k) ・・・(12
)ここでAO−T”・A・工、1BO=T−1・B。
−1>・・・(11) y(k)=Co−X(k) ・・・(12
)ここでAO−T”・A・工、1BO=T−1・B。
C0=C−Tであり、正則なTを適当に選んで、「−[
β1 β2・・・βn]” −・−(14)co=
[oo・・・1] ・・・(15)とするこ
とができる。そこで、式(10)のし行列をL=[−C
1−C2・・・−αn]Tとおいて、式(13)、(1
4)、(15)より、となり、有限整定オブザーバを設
計することができた。ここでI9口、CはA、IB、C
を相似変換をしたものであるが、この操作によっても状
態方程式による制御の正しさは保証されている。
β1 β2・・・βn]” −・−(14)co=
[oo・・・1] ・・・(15)とするこ
とができる。そこで、式(10)のし行列をL=[−C
1−C2・・・−αn]Tとおいて、式(13)、(1
4)、(15)より、となり、有限整定オブザーバを設
計することができた。ここでI9口、CはA、IB、C
を相似変換をしたものであるが、この操作によっても状
態方程式による制御の正しさは保証されている。
以上、システム同定により求めた状態方程式(1)等の
ベクトルA、 El、 CよりオブザーバP6を設計し
たが、以後、このオブザーバの出力を改めてX(k)と
表わすことにする。
ベクトルA、 El、 CよりオブザーバP6を設計し
たが、以後、このオブザーバの出力を改めてX(k)と
表わすことにする。
次に最適フィードバックゲイン「の求め方について説明
するが、最適フィードバックゲイン「を求める手法は、
例えば「線形システム制御理論」(前掲書)等に詳しい
ので、ここでは詳解は略して結果のみを示しておく。
するが、最適フィードバックゲイン「を求める手法は、
例えば「線形システム制御理論」(前掲書)等に詳しい
ので、ここでは詳解は略して結果のみを示しておく。
運転条件の諸量u (k )と運転状態の諸量り(k)
とについて、 δu(k)−u(k)−u(k−1>”・(17)δv
(k ) =v (k ) −v (k−1> ・(
18)とし、次の評価関数Jを最小にする最適制御入力
、即ち運転条件u”(k)を求めることが内燃機関M1
に関する付加積分型最適レギュレータとしての制御問題
を解くことになる。
とについて、 δu(k)−u(k)−u(k−1>”・(17)δv
(k ) =v (k ) −v (k−1> ・(
18)とし、次の評価関数Jを最小にする最適制御入力
、即ち運転条件u”(k)を求めることが内燃機関M1
に関する付加積分型最適レギュレータとしての制御問題
を解くことになる。
J−Σ[δす”(k)・Q・δす(k)IIO
+δu” (k ) l’lδU(k)]・・・ (
19) 尚、ここでG、IRは重みパラメータ行列を、kは制御
I開始時点をOとするサンプル回数を、各々示しており
、式(19)右辺はΩ、Rを対角行列とする所謂2次形
式表現である。
19) 尚、ここでG、IRは重みパラメータ行列を、kは制御
I開始時点をOとするサンプル回数を、各々示しており
、式(19)右辺はΩ、Rを対角行列とする所謂2次形
式表現である。
この時、最適なフィードバックゲイン「はF=−(lR
+E3”−P−B)−1−B” −P−A・・・(20
) として求められる。尚、式(20)におけるλ。
+E3”−P−B)−1−B” −P−A・・・(20
) として求められる。尚、式(20)におけるλ。
酊は各々、
P=A” −P−へ一へT ・P・百・ (百T ・P
・・・ 〈 23) の解である。尚、ここで式(19)の評価関数Jの意味
は内燃機関1に対する制御入力としての運転条件の1t
Wlu (k ) = [δθ δFR]の動きを制約
しつつ、制御出力としての運転状態の諸量y(k)、こ
こでは吸入空気量δAR,回転数δNを含む諸my (
k )の目標値s”(k)からの偏差を最小にしようと
意図したものである。運転条件の諸Mu (k )に対
する制約の重み付けは、重みパラメータ行列G、lRの
値によって変更することができる。従って、すでに求め
ておいた内燃機関1の動的なモデル、即ち行列A、B、
C(ここではA、B、C)を用い、任意の重みパラメー
タ行列Q、Rを選択して式(23)を解いてPを求め、
式(20)により最適フィードバックゲイン「を求めれ
ば、状態変数量X (k )は状態推定m’pi (h
>とじて式(9)より求められるので、u(k)=F
・ [X(k) ZN(k)ZAR(k )
] ・・・ (24)により内燃機関1にとって
の制御入力運転条件の諸り !A (k )を求めるこ
とができる。重みパラメータ行列Ω、Rを変えて最適な
制御特性が1qられるまで以上のシミュレーションを繰
返すことによって、最適フィードバックゲイン「が求め
られた。
・・・ 〈 23) の解である。尚、ここで式(19)の評価関数Jの意味
は内燃機関1に対する制御入力としての運転条件の1t
Wlu (k ) = [δθ δFR]の動きを制約
しつつ、制御出力としての運転状態の諸量y(k)、こ
こでは吸入空気量δAR,回転数δNを含む諸my (
k )の目標値s”(k)からの偏差を最小にしようと
意図したものである。運転条件の諸Mu (k )に対
する制約の重み付けは、重みパラメータ行列G、lRの
値によって変更することができる。従って、すでに求め
ておいた内燃機関1の動的なモデル、即ち行列A、B、
C(ここではA、B、C)を用い、任意の重みパラメー
タ行列Q、Rを選択して式(23)を解いてPを求め、
式(20)により最適フィードバックゲイン「を求めれ
ば、状態変数量X (k )は状態推定m’pi (h
>とじて式(9)より求められるので、u(k)=F
・ [X(k) ZN(k)ZAR(k )
] ・・・ (24)により内燃機関1にとって
の制御入力運転条件の諸り !A (k )を求めるこ
とができる。重みパラメータ行列Ω、Rを変えて最適な
制御特性が1qられるまで以上のシミュレーションを繰
返すことによって、最適フィードバックゲイン「が求め
られた。
以上、最小2乗法によるシステム同定により内燃機関1
の制御系の動的モデルの構築、有限整定オブザーバの設
計、最適フィードバックゲイン「の算出について説明し
たが、これらは予め求めておき、電子制御回路40の内
部ではその結果のみを用いて実際の制御を行なうのであ
る。
の制御系の動的モデルの構築、有限整定オブザーバの設
計、最適フィードバックゲイン「の算出について説明し
たが、これらは予め求めておき、電子制御回路40の内
部ではその結果のみを用いて実際の制御を行なうのであ
る。
そこで、次に、第7図のフローチャートに拠って電子制
御回路40が実際に行なう制御について説明する。尚、
以下の説明では現実の処理において扱われている量を添
字(k)付で、前回に扱われた量を添字(k−1)付で
表わすことにする。MPU44は内燃機関1が起動した
後、繰返しステップ100以下の処理を繰返している。
御回路40が実際に行なう制御について説明する。尚、
以下の説明では現実の処理において扱われている量を添
字(k)付で、前回に扱われた量を添字(k−1)付で
表わすことにする。MPU44は内燃機関1が起動した
後、繰返しステップ100以下の処理を繰返している。
まずステップ100では、内燃機関1の運転状態、即ち
吸入空気流量AR(k−1> 、回転数N(k−4)等
を各センサより読み込む処理を行なう。
吸入空気流量AR(k−1> 、回転数N(k−4)等
を各センサより読み込む処理を行なう。
続くステップ110では、冷却水温Thwなどに基づい
て内燃機関1の目標アイドル回転数N*を算出し、ステ
ップ120では内燃機関1の目標吸入空気量AR*を算
出する。目標吸入空気量AR*は、内燃機関1の燃料消
費量を最小にすべく定められるものであって、その算出
は第8図に拠って後述する制御により行なわれる。これ
らの処理が第4図P1.P2の各設定部に相当する。
て内燃機関1の目標アイドル回転数N*を算出し、ステ
ップ120では内燃機関1の目標吸入空気量AR*を算
出する。目標吸入空気量AR*は、内燃機関1の燃料消
費量を最小にすべく定められるものであって、その算出
は第8図に拠って後述する制御により行なわれる。これ
らの処理が第4図P1.P2の各設定部に相当する。
ステップ130では、この目標アイドル回転数N本と実
際に検出された回転数N(k−1)の偏差SNを、また
、この目標吸入空気量AR”と実際の吸入空気量AR(
k−1>との偏差SARを、各々求める処理が行なわれ
る。続くステップ140では、ステップ130で求めた
各偏差を累積する処理、即ち ZN (k ) −ZN (k−1) +SN (k−
1>により累積値ZN(k)を、一方ZAR(k )=
ZAR(k−1> +SAR(k−1> ニ、にす累積
値ZAR(k)を求める処理が行なわれる。この処理が
第4図の積分器P3.P4に相当する。
際に検出された回転数N(k−1)の偏差SNを、また
、この目標吸入空気量AR”と実際の吸入空気量AR(
k−1>との偏差SARを、各々求める処理が行なわれ
る。続くステップ140では、ステップ130で求めた
各偏差を累積する処理、即ち ZN (k ) −ZN (k−1) +SN (k−
1>により累積値ZN(k)を、一方ZAR(k )=
ZAR(k−1> +SAR(k−1> ニ、にす累積
値ZAR(k)を求める処理が行なわれる。この処理が
第4図の積分器P3.P4に相当する。
続くステップ150では、ステップ100で読み込んだ
内燃機関1の運転状態から、内燃機関1の動的モデルを
構築した際、線形近似が成立つ範囲として取り上げた定
常的な運転状態のうちで最も近い状態(以下、これを定
常点ARa、Naと呼ぶ)を求める処理を行なう。ステ
ップ160では、ステップ100で読み込んだ内燃機関
1の運転状態をこの定常点(ARa、Na>からの摂動
分(δAR,δN)として求める処理を行なう。
内燃機関1の運転状態から、内燃機関1の動的モデルを
構築した際、線形近似が成立つ範囲として取り上げた定
常的な運転状態のうちで最も近い状態(以下、これを定
常点ARa、Naと呼ぶ)を求める処理を行なう。ステ
ップ160では、ステップ100で読み込んだ内燃機関
1の運転状態をこの定常点(ARa、Na>からの摂動
分(δAR,δN)として求める処理を行なう。
この処理が第4図摂動分抽出部P5に相当する。
続くステップ170では、内燃機関1の冷却水mThw
を読み込み、この水IThwに応じて内燃機関1の動的
モデルが変化することから、予め冷却水mThw毎に用
意されたオブザーバ内のパラメータAO,IBO,L及
び最適フィードバックゲイン「を選択する処理を行なう
。
を読み込み、この水IThwに応じて内燃機関1の動的
モデルが変化することから、予め冷却水mThw毎に用
意されたオブザーバ内のパラメータAO,IBO,L及
び最適フィードバックゲイン「を選択する処理を行なう
。
ステップ180では、ステップ170で選択されたAO
,EIO,Lとステップ160で求められた摂動分(δ
AR,δN)と、前回求められた状態変数推定IX (
k−1> = [xl (k−1) X2(k−1
> −x4 (k−1) ]”と、前回求められた
燃料噴射IFR(k−1)とスロットルバルブ開度θ(
k−1)の摂動分δFR(k−1>、δθ(k−1)と
から、次式(25)により新たな状態変数推定11X
(k )を算出する処理が行なわれる。この処理が第4
図オブザーバP6に相当するが、本実施例では既述した
如く、オブザーバP6は有限整定オブザーバとして構成
されている。即ち、X(k) −AO−X (k−1) +Bo ・[δFRδθ]+
L・[δAR(k−1)δN(k−1)]・・・(25
) の計算が行なわれるのである。
,EIO,Lとステップ160で求められた摂動分(δ
AR,δN)と、前回求められた状態変数推定IX (
k−1> = [xl (k−1) X2(k−1
> −x4 (k−1) ]”と、前回求められた
燃料噴射IFR(k−1)とスロットルバルブ開度θ(
k−1)の摂動分δFR(k−1>、δθ(k−1)と
から、次式(25)により新たな状態変数推定11X
(k )を算出する処理が行なわれる。この処理が第4
図オブザーバP6に相当するが、本実施例では既述した
如く、オブザーバP6は有限整定オブザーバとして構成
されている。即ち、X(k) −AO−X (k−1) +Bo ・[δFRδθ]+
L・[δAR(k−1)δN(k−1)]・・・(25
) の計算が行なわれるのである。
続くステップ190では、ステップ180で求めた状態
変数推定fMX (k )と、ステップ140で求めら
れた累積値ZN (k ) 、 ZAR(k )と、予
め用意されステップ170で選択された最適フィードバ
ックゲイン、 をベクトル乗算することにより、即ち[δFR△ (k ) δθ(k ) ] =F・[X (k )
ZN(k ’) ZAR(k ) 1 ”により
操作量の摂動分δFR(k )とδθ(k)とを求める
処理が行なわれる。これが第4図のフィードバック量決
定部P7に相当している。
変数推定fMX (k )と、ステップ140で求めら
れた累積値ZN (k ) 、 ZAR(k )と、予
め用意されステップ170で選択された最適フィードバ
ックゲイン、 をベクトル乗算することにより、即ち[δFR△ (k ) δθ(k ) ] =F・[X (k )
ZN(k ’) ZAR(k ) 1 ”により
操作量の摂動分δFR(k )とδθ(k)とを求める
処理が行なわれる。これが第4図のフィードバック量決
定部P7に相当している。
ステップ200では、ステップ190で求めた操作量の
摂動分δFR(k ) 、δθ(k)と定常点における
各操作量FRa、θaとを加えて、実際に内燃機関1の
燃料噴射弁11及びアクチュエータ35へ出力される操
作量、即ち運転条件FR(k)、θ(k)が求められる
。
摂動分δFR(k ) 、δθ(k)と定常点における
各操作量FRa、θaとを加えて、実際に内燃機関1の
燃料噴射弁11及びアクチュエータ35へ出力される操
作量、即ち運転条件FR(k)、θ(k)が求められる
。
続くステップ210ではサンプリング回数を示ti+1
1kを1だけインクリメントする処理を行ない、上記の
一連の処理、ステップ100ないしステップ210を終
わる。
1kを1だけインクリメントする処理を行ない、上記の
一連の処理、ステップ100ないしステップ210を終
わる。
以上の制御を周期的に継続して行なうことにより、電子
制御回路40は内燃機関1の運転状態を目標アイドル回
転数N*と目標吸入空気量AR”に制御する付加積分型
最適レギュレータとして、最適のフィードバックゲイン
により制御を行なうことになる。
制御回路40は内燃機関1の運転状態を目標アイドル回
転数N*と目標吸入空気量AR”に制御する付加積分型
最適レギュレータとして、最適のフィードバックゲイン
により制御を行なうことになる。
次に、ステップ120の目標吸入空気量AR”を求める
ルーチンについて説明する。このルーチンは第8図のフ
ローチャートに図示する如く、次の手順により、同一の
アイドル回転数N (k )を維持しつつ、燃料消費量
を最小にするように目標吸入空気MAR”を算出する。
ルーチンについて説明する。このルーチンは第8図のフ
ローチャートに図示する如く、次の手順により、同一の
アイドル回転数N (k )を維持しつつ、燃料消費量
を最小にするように目標吸入空気MAR”を算出する。
尚、以下の説明では、本ルーチンにおける前回の目標値
をAR”(k−1)で、今回算出された目標値をAR*
(k )で表わすことがある。
をAR”(k−1)で、今回算出された目標値をAR*
(k )で表わすことがある。
このルーチンはステップ300より開始され、まず第6
図の処理において定められた目標アイドル回転数N本(
k)が前回の値N* (k−1)と、実際の回転数N
(k )が目標アイドル回転数N*(k)と、それぞれ
等しいか否かの判断を行なう。
図の処理において定められた目標アイドル回転数N本(
k)が前回の値N* (k−1)と、実際の回転数N
(k )が目標アイドル回転数N*(k)と、それぞれ
等しいか否かの判断を行なう。
それらのいずれかひとつでも成立しないような場合には
、制御系は平衡に達していないことから、燃料消費量を
最小とするような吸入空気量の探索は行なえないとして
、処理はステップ310へ移行し、内燃機関1の回転数
Nから予め設定されたマツプにより与えられる吸入空気
MAR(N)を目標吸入空気量AR” (k )として
与える処理を行ない、NEXTへ抜けて、本制御ルーチ
ンを終了する。即ち、第7図のフローチャートに還って
言えば、ステップ120において目標吸入空気量AR*
(k )を内燃Il関1が過渡状態にあるとして、マツ
プより定めるのである。
、制御系は平衡に達していないことから、燃料消費量を
最小とするような吸入空気量の探索は行なえないとして
、処理はステップ310へ移行し、内燃機関1の回転数
Nから予め設定されたマツプにより与えられる吸入空気
MAR(N)を目標吸入空気量AR” (k )として
与える処理を行ない、NEXTへ抜けて、本制御ルーチ
ンを終了する。即ち、第7図のフローチャートに還って
言えば、ステップ120において目標吸入空気量AR*
(k )を内燃Il関1が過渡状態にあるとして、マツ
プより定めるのである。
一方ステップ300において、N*(k)=N宰(k−
1)、かつN(k)−N”(k)であれば、内燃i関1
は平衡状態にあるとみなし、燃料消費量を最小とする吸
入空気量の探索を行なうために、処理はステップ320
へ移行する。ステップ320ではフラッグFsが1であ
るか否かの判断を行なうが燃料消費量の探索が開始され
る以前はフラッグFSの値はOに初期セットされている
ので判断は1No」となって処理はステップ330へ進
む。ステップ330では、回転数N (k )を最小の
燃料消費量で目標アイドル回転数N*(k)に維持でき
る吸入空気量を求める探索を開始するとして、フラッグ
Fsを値1に、探索方向を示す係数りつまり吸入空気量
の増加方向または減少方向を指定する係数りを値1に、
本ルーチンの実行回数を示すカウンタC5を値0に、各
々設定する処理を行なう。
1)、かつN(k)−N”(k)であれば、内燃i関1
は平衡状態にあるとみなし、燃料消費量を最小とする吸
入空気量の探索を行なうために、処理はステップ320
へ移行する。ステップ320ではフラッグFsが1であ
るか否かの判断を行なうが燃料消費量の探索が開始され
る以前はフラッグFSの値はOに初期セットされている
ので判断は1No」となって処理はステップ330へ進
む。ステップ330では、回転数N (k )を最小の
燃料消費量で目標アイドル回転数N*(k)に維持でき
る吸入空気量を求める探索を開始するとして、フラッグ
Fsを値1に、探索方向を示す係数りつまり吸入空気量
の増加方向または減少方向を指定する係数りを値1に、
本ルーチンの実行回数を示すカウンタC5を値0に、各
々設定する処理を行なう。
続くステップ340ではカウンタCsの値がOを越えて
いるか否かの判断を行なう。探索開始直後にはカウンタ
C5=Oなので処理はステップ350に移行し、目標吸
入空気WAR*(k ’)を前向の目標値AR” (k
−1)よりD×ΔARだけ増量して定め、続くステップ
360でカウンタCsの値を1だけインクリメントし、
NEXTへ抜けて本ルーチンを終了する。
いるか否かの判断を行なう。探索開始直後にはカウンタ
C5=Oなので処理はステップ350に移行し、目標吸
入空気WAR*(k ’)を前向の目標値AR” (k
−1)よりD×ΔARだけ増量して定め、続くステップ
360でカウンタCsの値を1だけインクリメントし、
NEXTへ抜けて本ルーチンを終了する。
こうして探索が開始された後に、本ルーチンが実行され
るとステップ320.ステップ340での判断は共にr
YESJとなって、処理はステップ370へ進み、燃料
噴射量FR(k )に関し、定常点からの摂動分δFR
(k )が前回の時点でのδFR(k−1>と較べて、
どうなったかを判定する。δFR(k )−δFR(k
−1>の値が所定値−ΔF以下であれば燃料噴射量は更
に減少し得るとして、そのまま探索を継続すべく、処理
はステップ350以下を実行する。このことは丁度第5
図においてb点側からa点側へ近づきつつあることを示
している。
るとステップ320.ステップ340での判断は共にr
YESJとなって、処理はステップ370へ進み、燃料
噴射量FR(k )に関し、定常点からの摂動分δFR
(k )が前回の時点でのδFR(k−1>と較べて、
どうなったかを判定する。δFR(k )−δFR(k
−1>の値が所定値−ΔF以下であれば燃料噴射量は更
に減少し得るとして、そのまま探索を継続すべく、処理
はステップ350以下を実行する。このことは丁度第5
図においてb点側からa点側へ近づきつつあることを示
している。
一方、δFR(k )−δFR(k−1)の値が所定値
61以上であれば、燃料噴射量は増加していることにな
るので探索方向を逆転すべく、ステップ380にて探索
方向係数りの値を−1に設定し、以下、上述のステップ
350,360の処理を行なう。従って、この探索によ
り目標吸入空気IAR*(k)は減少され、スロットル
開度θは閉方向に向う。第2図に即して言えば、図中0
点側からa点側へ探索する場合である。
61以上であれば、燃料噴射量は増加していることにな
るので探索方向を逆転すべく、ステップ380にて探索
方向係数りの値を−1に設定し、以下、上述のステップ
350,360の処理を行なう。従って、この探索によ
り目標吸入空気IAR*(k)は減少され、スロットル
開度θは閉方向に向う。第2図に即して言えば、図中0
点側からa点側へ探索する場合である。
こうして、燃料噴射量を減らす方向での探索が行なわれ
ると、やがて、δFR(k )−δFR(k−1)の値
が所定の偏差±ΔF以内となる点が見い出される。ここ
が同一目標アイドル回転数において燃料消費量を最小と
する吸入空気量となる点である。そこで、探索は一応終
了したとして、ステップ390にてフラッグl”sを値
0に設定し、続くステップ400では、この時得られて
いた目標吸入空気量AR*(k−1>を回転数Nから吸
入空気量を定めるマツプの値として入替える処理を行な
う。即ちAR(N)−AR” (k−1>である。
ると、やがて、δFR(k )−δFR(k−1)の値
が所定の偏差±ΔF以内となる点が見い出される。ここ
が同一目標アイドル回転数において燃料消費量を最小と
する吸入空気量となる点である。そこで、探索は一応終
了したとして、ステップ390にてフラッグl”sを値
0に設定し、続くステップ400では、この時得られて
いた目標吸入空気量AR*(k−1>を回転数Nから吸
入空気量を定めるマツプの値として入替える処理を行な
う。即ちAR(N)−AR” (k−1>である。
続くステン7’410では前回定められた目標吸入空気
量AR” (k−1>を今回も用いるとして、AR”(
k)の値をこれに更新し、NEXTへ抜けて本ルーチン
を終了する。
量AR” (k−1>を今回も用いるとして、AR”(
k)の値をこれに更新し、NEXTへ抜けて本ルーチン
を終了する。
以上で、−回の探索過程を終了し、その後は冒頭の処理
、ステップ320,330.340より再び探索を継続
するのである。
、ステップ320,330.340より再び探索を継続
するのである。
以上、説明した如く、第7図、第8図に示した制御ルー
チンを繰返し実行することにより、本実施例の内燃機関
制御装置は、内燃機関1の回転数を目標アイドル回転数
に制御するのみならず、その燃料消費量を最小にするよ
うに働く。この時、内燃機関1を制御する系はそのフィ
ードバックゲインが最適フィードバックとなる付加積分
型最適レギュレータとなっており、スロットルバルブ開
度θと燃料噴射ff1FRとの制御は、従来実現不可能
であった素早い応答性と安定性とにおいて実現されてい
るり従って、内燃機関1の運転者の運転フィーリグをす
こしも損うことなく、スロットルバルブ7の開度θを変
化させて、燃料噴射量FRを最小にする制御が可能とな
った。
チンを繰返し実行することにより、本実施例の内燃機関
制御装置は、内燃機関1の回転数を目標アイドル回転数
に制御するのみならず、その燃料消費量を最小にするよ
うに働く。この時、内燃機関1を制御する系はそのフィ
ードバックゲインが最適フィードバックとなる付加積分
型最適レギュレータとなっており、スロットルバルブ開
度θと燃料噴射ff1FRとの制御は、従来実現不可能
であった素早い応答性と安定性とにおいて実現されてい
るり従って、内燃機関1の運転者の運転フィーリグをす
こしも損うことなく、スロットルバルブ7の開度θを変
化させて、燃料噴射量FRを最小にする制御が可能とな
った。
しかも本実施例では、内燃機関1の冷却水温Thwに応
じて動的なモデルが変化する為に、冷却水温Thwによ
ってオブザーバのパラメータと最適フィードバックゲイ
ンを切換えて制御を行なっており、内燃機関1の冷却水
温Thwによらず安定した制御を行なうことができる。
じて動的なモデルが変化する為に、冷却水温Thwによ
ってオブザーバのパラメータと最適フィードバックゲイ
ンを切換えて制御を行なっており、内燃機関1の冷却水
温Thwによらず安定した制御を行なうことができる。
こうした優れた応答性と安定性とが実現できて初めて、
内燃機関1の燃料噴射IFRを最小とするような探索が
可能となった。従来のフィードバック制御によってアク
チュエータ35を介してスロットルバルブ7を駆動した
としても、探索を行なうことはできても、応答性・安定
性などの面から現実の使用に耐えるものではないからで
ある。
内燃機関1の燃料噴射IFRを最小とするような探索が
可能となった。従来のフィードバック制御によってアク
チュエータ35を介してスロットルバルブ7を駆動した
としても、探索を行なうことはできても、応答性・安定
性などの面から現実の使用に耐えるものではないからで
ある。
第9図はこれを比較したものであり、一点鎖線rはアイ
ドル回転数の目標値N* (k )を、実線9は本実施
例の制御を行なった場合のm関口転数N (k )の−
例を、破線すは従来のフィードバック制御を行なった場
合の機関回転数N (k )の−例を、各々示している
。図より明らかな如く、付加積分型最適レギュレータと
して構成された本実施例の内燃機関制御a置によれば、
従来のフィードバック制御より早い応答性(立ち上がり
)を実現した上でオーバーシュート、ダウンシュートも
ほとんど生じることなくアイドル回転数の制御がなされ
ている。内燃機関1の機関回転数が平衡に達するまでの
時間で比較すれば1桁以上の改善が実現されており、こ
のことが燃料噴射量を最小とする吸入空気量の探索を現
実的なものにしている。
ドル回転数の目標値N* (k )を、実線9は本実施
例の制御を行なった場合のm関口転数N (k )の−
例を、破線すは従来のフィードバック制御を行なった場
合の機関回転数N (k )の−例を、各々示している
。図より明らかな如く、付加積分型最適レギュレータと
して構成された本実施例の内燃機関制御a置によれば、
従来のフィードバック制御より早い応答性(立ち上がり
)を実現した上でオーバーシュート、ダウンシュートも
ほとんど生じることなくアイドル回転数の制御がなされ
ている。内燃機関1の機関回転数が平衡に達するまでの
時間で比較すれば1桁以上の改善が実現されており、こ
のことが燃料噴射量を最小とする吸入空気量の探索を現
実的なものにしている。
従って、マクロ的に観た場合、常時、内燃機関1の燃料
消費量は最小に制御されている。
消費量は最小に制御されている。
以上の実施例では、内燃機関1を燃料噴@量FRとスロ
ットル開度θを入力とし吸入空気量AR。
ットル開度θを入力とし吸入空気量AR。
回転数Nを出力とする2人力2出力の系としてとらえ、
最小2乗法によるシステム同定を用いて動的なモデルを
構築して付加積分型最適レギュレータを構成しているが
、適用する内燃機関の態様に合わせて、この他の入出力
をも加味し、系の動的なモデルを構築することも、本発
明の要旨を変更することなく行なうことができる。
最小2乗法によるシステム同定を用いて動的なモデルを
構築して付加積分型最適レギュレータを構成しているが
、適用する内燃機関の態様に合わせて、この他の入出力
をも加味し、系の動的なモデルを構築することも、本発
明の要旨を変更することなく行なうことができる。
[発明の効果]
以上詳述したように、本発明の内燃機関の運転状態制御
装置は、目標吸入空気量を、目標アイドル回転数を一定
とした時の吸入空気量と燃料供給量との相関に基づいて
燃料供給量を最小とする吸入空気量に定め、一方、その
制御手段を、内燃機関の運転に関する系の動的なモデル
に従って予め定められた最適フィードバックゲインに基
づいてフィードバック量を定める付加積分型最適レギ1
レータとして構成している。
装置は、目標吸入空気量を、目標アイドル回転数を一定
とした時の吸入空気量と燃料供給量との相関に基づいて
燃料供給量を最小とする吸入空気量に定め、一方、その
制御手段を、内燃機関の運転に関する系の動的なモデル
に従って予め定められた最適フィードバックゲインに基
づいてフィードバック量を定める付加積分型最適レギ1
レータとして構成している。
従って、従来のスロットルアクチュエータ付内燃i閏で
は得られなかった高い応答性と安定性とを実現しながら
、内燃機関の回転数を目標アイドル回転数に制御すると
共に、その燃料消費量を最小にすることができるという
優れた効果を奏する。
は得られなかった高い応答性と安定性とを実現しながら
、内燃機関の回転数を目標アイドル回転数に制御すると
共に、その燃料消費量を最小にすることができるという
優れた効果を奏する。
従って車載の内燃機関に適用すれば、アイドル回転数を
安定に制御し、ドライブフィーリングを十分に好適なも
のとすることができるばかりか、車両の燃費も大幅に改
善されるなど、制御特性を格段に向上させることができ
る。
安定に制御し、ドライブフィーリングを十分に好適なも
のとすることができるばかりか、車両の燃費も大幅に改
善されるなど、制御特性を格段に向上させることができ
る。
第1図は本発明の基本的構成図、第2図は燃料量FRと
吸入空気量ARとの関係を示す定アイドル回転数If図
、第3図は本発明一実施例としての内燃機関とその周辺
装置の構成を示す概略構成図、第4図は同じくその制御
系統図、第5図は実施例の系のモデルを同定するのに用
いたブロック線図、第6図は伝達関数を求める為のシグ
ナルフロー線図、第7図は実施例における付加積分型最
適レギュレータとしての制御を示すフローチャート、第
8図は同じくその燃料消費量を最小にする制御ルーチン
を示すフローチャート、第9図は実施例との制御特性と
従来の制御の一例とを比較するグラフ、である。 1・・・内燃機関 3・・・エア70メータ 7・・・スロットルバルブ 11・・・燃料噴射弁 31・・・回転数センサ 40・・・電子制御回路 44・・・MPU
吸入空気量ARとの関係を示す定アイドル回転数If図
、第3図は本発明一実施例としての内燃機関とその周辺
装置の構成を示す概略構成図、第4図は同じくその制御
系統図、第5図は実施例の系のモデルを同定するのに用
いたブロック線図、第6図は伝達関数を求める為のシグ
ナルフロー線図、第7図は実施例における付加積分型最
適レギュレータとしての制御を示すフローチャート、第
8図は同じくその燃料消費量を最小にする制御ルーチン
を示すフローチャート、第9図は実施例との制御特性と
従来の制御の一例とを比較するグラフ、である。 1・・・内燃機関 3・・・エア70メータ 7・・・スロットルバルブ 11・・・燃料噴射弁 31・・・回転数センサ 40・・・電子制御回路 44・・・MPU
Claims (1)
- 【特許請求の範囲】 1 内燃機関の吸気管に配設される燃料噴射弁、および
吸気通路のスロットルバルブの開度を調節するスロット
ルアクチュエータと、 機関回転数検出手段、吸入空気流量検出手段および機関
制御のためのその他の検出手段と、これら検出手段から
の検出信号に基づき上記燃料噴射弁へ燃料噴射弁駆動信
号を、上記スロットルアクチュエータへスロットル開度
を規定するアクチュエータ駆動信号をそれぞれ出力する
電子制御手段であって、目標アイドル回転数を設定する
アイドル回転数設定部を備えるものと、 からなる内燃機関制御装置において、 上記電子制御手段を、 上記アイドル回転数設定部による目標アイドル回転数の
情報を受け、機関回転数をこの目標アイドル回転数に維
持できる最小の燃料噴射量に対応する吸入空気量を、目
標回転数に対応して予め定めた燃料噴射量・吸入空気量
相関パターンに基づいて定める目標吸入空気量設定部と
、この目標吸入空気量設定部による目標吸入空気量の情
報、上記アイドル回転数設定部による目標回転数の情報
および内燃機関の制御入出力値をそれぞれ受け、内燃機
関の運転に関する系の動的なモデルに従って予め定めら
れた機関内部状態推定のための演算式および最適フィー
ドバックゲインに基づいて、燃料噴射量およびスロット
ル開度のそれぞれについてのフィードバック量を決定し
、これら決定されたフィードバック量のそれぞれに基づ
く燃料噴射弁駆動信号およびアクチュエータ駆動信号を
出力する付加積分型最適レギュレータと、 から構成した ことを特徴とする内燃機関制御装置。 2 前記付加積分型最適レギュレータが、 前記内燃機関の運転に関する系の動的なモデルに基づい
て予め設定されたパラメータを用いて、前記内燃機関の
運転状態と運転条件とから該系の動的な内部状態を表わ
す適当な次数の状態変数量を推定する状態観測部と、 前記目標吸入空気量設定部、アイドル回転数設定部によ
ってそれぞれ定められた目標吸入空気量、目標アイドル
回転数と前記検出された運転状態諸量との各偏差を各々
累積する累積部と、 前記系の動的なモデルに基づいて予め設定されたフィー
ドバックゲインと前記推定された状態変数量と前記累積
値とから、燃料供給量、スロットル開度の各操作量を決
定するフィードバック量決定部と、 から構成された特許請求の範囲第1項記載の内燃機関制
御装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60031265A JPH0612093B2 (ja) | 1985-02-19 | 1985-02-19 | 内燃機関制御装置 |
| DE19863605282 DE3605282A1 (de) | 1985-02-19 | 1986-02-19 | Regelvorrichtung fuer brennkraftmaschinen |
| US06/830,935 US4638778A (en) | 1985-02-19 | 1986-02-19 | Idle speed control apparatus for internal combustion engine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60031265A JPH0612093B2 (ja) | 1985-02-19 | 1985-02-19 | 内燃機関制御装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS61190143A true JPS61190143A (ja) | 1986-08-23 |
| JPH0612093B2 JPH0612093B2 (ja) | 1994-02-16 |
Family
ID=12326511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60031265A Expired - Lifetime JPH0612093B2 (ja) | 1985-02-19 | 1985-02-19 | 内燃機関制御装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US4638778A (ja) |
| JP (1) | JPH0612093B2 (ja) |
| DE (1) | DE3605282A1 (ja) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0645311B2 (ja) * | 1985-02-18 | 1994-06-15 | 日本電装株式会社 | 車速制御装置 |
| JPH0660593B2 (ja) * | 1985-08-05 | 1994-08-10 | 株式会社日立製作所 | 電子式内燃機関制御装置 |
| US4785780A (en) * | 1986-07-08 | 1988-11-22 | Nippondenso Co., Ltd. | Control apparatus |
| JP2674077B2 (ja) * | 1988-04-12 | 1997-11-05 | トヨタ自動車株式会社 | 内燃機関の非線形フィードバック制御方法 |
| JP2748488B2 (ja) * | 1989-01-18 | 1998-05-06 | 株式会社デンソー | スロットル開度制御装置 |
| IT1241215B (it) * | 1990-05-07 | 1993-12-29 | Fiat Auto Spa | Procedimento ed apparato per il controllo della velocita' di rotazione al minimo di un motore a combustione interna. |
| JP3064346B2 (ja) * | 1990-07-24 | 2000-07-12 | 株式会社デンソー | エンジンの回転数制御装置 |
| DE69226709T2 (de) * | 1991-06-10 | 1999-05-06 | Denso Corp., Kariya, Aichi | Vorrichtung zur Regelung der Drehzahl einer Brennkraftmaschine |
| JP3005455B2 (ja) * | 1995-06-14 | 2000-01-31 | トヨタ自動車株式会社 | 内燃機関の回転数制御装置 |
| US5577474A (en) * | 1995-11-29 | 1996-11-26 | General Motors Corporation | Torque estimation for engine speed control |
| DE19741965C1 (de) * | 1997-09-23 | 1999-01-21 | Siemens Ag | Verfahren zur Laufruheregelung |
| US6742497B1 (en) * | 1999-04-06 | 2004-06-01 | Toyota Jidosha Kabushiki Kaisha | Device for controlling rotational speed of internal combustion engine |
| DE102011110607A1 (de) * | 2011-08-18 | 2013-02-21 | Rauch Landmaschinenfabrik Gmbh | Zweischeibendüngerstreuer |
| DE102012003581B3 (de) * | 2012-02-27 | 2013-07-18 | Iav Gmbh Ingenieurgesellschaft Auto Und Verkehr | Leerlaufregler und Verfahren zum Betrieb von Brennkraftmaschinen |
| CN121066725B (zh) * | 2025-09-16 | 2026-07-03 | 江苏奥易克斯汽车电子科技股份有限公司 | 基于动态怠速偏差补偿的发动机怠速转速控制方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5614836A (en) * | 1979-07-13 | 1981-02-13 | Hitachi Ltd | Controlling device for internal combustion engine |
| JPS5641431A (en) * | 1979-09-10 | 1981-04-18 | Hitachi Ltd | Controller for internal combustion engine |
| JPS5912860A (ja) * | 1982-07-13 | 1984-01-23 | Fujitsu Ltd | ワイヤドツトプリンタ用ワイヤガイドの製法 |
| JPS5943943A (ja) * | 1982-09-06 | 1984-03-12 | Nissan Motor Co Ltd | 内燃機関のアイドル回転速度制御方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2903977A1 (de) * | 1979-02-02 | 1980-08-07 | Bosch Gmbh Robert | Einrichtung zum steuern von betriebsparameterabhaengigen und sich wiederholenden vorgaengen bei brennkraftmaschinen |
| DE3268810D1 (en) * | 1981-08-01 | 1986-03-13 | Nippon Denso Co | Method and system for controlling ignition timing in a multicylinder internal combustion engine |
| US4475506A (en) * | 1981-09-25 | 1984-10-09 | Brunswick Corporation | Programmable fuel economy optimizer for an internal combustion engine |
| US4562814A (en) * | 1983-02-04 | 1986-01-07 | Nissan Motor Company, Limited | System and method for controlling fuel supply to an internal combustion engine |
| JPS59188052A (ja) * | 1983-04-08 | 1984-10-25 | Nippon Denso Co Ltd | 内燃機関の空燃比制御方法 |
| US4513713A (en) * | 1983-09-06 | 1985-04-30 | Honda Giken Kogyo Kabushiki Kaisha | Method of controlling operating amounts of operation control means for an internal combustion engine |
-
1985
- 1985-02-19 JP JP60031265A patent/JPH0612093B2/ja not_active Expired - Lifetime
-
1986
- 1986-02-19 US US06/830,935 patent/US4638778A/en not_active Expired - Lifetime
- 1986-02-19 DE DE19863605282 patent/DE3605282A1/de active Granted
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5614836A (en) * | 1979-07-13 | 1981-02-13 | Hitachi Ltd | Controlling device for internal combustion engine |
| JPS5641431A (en) * | 1979-09-10 | 1981-04-18 | Hitachi Ltd | Controller for internal combustion engine |
| JPS5912860A (ja) * | 1982-07-13 | 1984-01-23 | Fujitsu Ltd | ワイヤドツトプリンタ用ワイヤガイドの製法 |
| JPS5943943A (ja) * | 1982-09-06 | 1984-03-12 | Nissan Motor Co Ltd | 内燃機関のアイドル回転速度制御方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE3605282C2 (ja) | 1993-07-01 |
| US4638778A (en) | 1987-01-27 |
| DE3605282A1 (de) | 1986-08-21 |
| JPH0612093B2 (ja) | 1994-02-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS61145339A (ja) | 内燃機関の運転状態制御装置 | |
| JPS61191435A (ja) | 車速制御装置 | |
| JPS62147027A (ja) | 内燃機関駆動車両の加速制御装置 | |
| Hendricks et al. | SI engine controls and mean value engine modelling | |
| EP1643101B1 (en) | Intake air amount control device of internal combustion engine and control device | |
| Powell et al. | Observer-based air fuel ratio control | |
| JP3510021B2 (ja) | 内燃機関の空燃比制御装置 | |
| EP0286104B1 (en) | Method of controlling fuel supply to engine by prediction calculation | |
| JPH0612093B2 (ja) | 内燃機関制御装置 | |
| JPS6349060B2 (ja) | ||
| Kang et al. | Idle speed controller based on active disturbance rejection control in diesel engines | |
| EP0115868A2 (en) | System and method for contolling fuel supply to an internal combustion engine | |
| Takahashi et al. | A simple engine model for idle speed control | |
| EP1645740B1 (en) | Intake airvolume controller of internal combustion engine | |
| US4747385A (en) | Air-fuel ratio control system for an automotive engine | |
| JPS6325173B2 (ja) | ||
| JP2510186B2 (ja) | 制御装置 | |
| JPS6325177B2 (ja) | ||
| Sweet | Automotive applications of modern control theory | |
| JPS6340929B2 (ja) | ||
| JPS6245957A (ja) | 内燃機関の制御装置 | |
| JPS6328226B2 (ja) | ||
| Keller | Two-stage model predictive control for the air path of a turbocharged gasoline engine with exhaust gas recirculation | |
| JPS6345438A (ja) | エンジン制御装置 | |
| Alt et al. | Control of idle engine speed and torque reserve with higher order sliding modes |