JPS59145338A - Control method of idling speed for internal-combustion engine - Google Patents

Control method of idling speed for internal-combustion engine

Info

Publication number
JPS59145338A
JPS59145338A JP58018100A JP1810083A JPS59145338A JP S59145338 A JPS59145338 A JP S59145338A JP 58018100 A JP58018100 A JP 58018100A JP 1810083 A JP1810083 A JP 1810083A JP S59145338 A JPS59145338 A JP S59145338A
Authority
JP
Japan
Prior art keywords
control
combustion engine
engine
speed
rotation speed
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
Application number
JP58018100A
Other languages
Japanese (ja)
Other versions
JPS6340928B2 (en
Inventor
Takashi Ueno
植野 隆司
Toru Takahashi
徹 高橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP58018100A priority Critical patent/JPS59145338A/en
Publication of JPS59145338A publication Critical patent/JPS59145338A/en
Publication of JPS6340928B2 publication Critical patent/JPS6340928B2/ja
Granted legal-status Critical Current

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D31/00Use of speed-sensing governors to control combustion engines, not otherwise provided for
    • F02D31/001Electric control of rotation speed
    • F02D31/002Electric control of rotation speed controlling air supply
    • F02D31/003Electric control of rotation speed controlling air supply for idle speed control

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)
  • Electrical Control Of Ignition Timing (AREA)
  • Exhaust-Gas Circulating Devices (AREA)

Abstract

PURPOSE:To aim at a distinct improvement in antistallingness, by discriminating an idling speed with a ratio between an idle switch and a car speed and an engine speed. CONSTITUTION:A something modellizing the dynamic characteristics of an internal-combustion engine 12, an actuator and a sensor is stored in a controller, while a combination of a bypass air quantity PA, ignition timing IT, a fuel feed quantity and an exhaust reflux quantity is set down to the control input and also an idling speed is set down to the control input whereby a state variable representing an inner state is estimated from the input and output of control. Then, using an integral value of a deviation SA between a desired value Nr and an actual value N of the idling speed, a control input valve is determined and the idling speed of the internal-combustion engine 12 is feedback-controlled to meet the desired value, through which a judgement for whether the idling speed is right or wrong is attained with a ratio between a throttle valve switch and a car speed plus an engine speed. Thus, an antistallingness at high-speed driving is improved.

Description

【発明の詳細な説明】 (技術分野) この発明は、内燃機関のアイドル時の回転速度の制御方
法に関し、より詳細には、従来一般的なPID(比例積
分微分)制御とは異なり、内燃機関の内部状態を考慮し
て機関をダイナミック(動的)なシステムとして捕え、
内部状態を規定する状態変数によって機関の動的な振舞
いを推定しながら、機関の入力変数を決定する状態変数
制御の手法を用いて、アイドル回転速度を制御する方法
に関する。
Detailed Description of the Invention (Technical Field) The present invention relates to a method for controlling the rotational speed of an internal combustion engine during idling, and more specifically, unlike conventional PID (proportional integral derivative) control, Considering the internal state of the engine, we view the institution as a dynamic system,
The present invention relates to a method of controlling idle rotation speed using a state variable control technique that determines engine input variables while estimating the dynamic behavior of the engine using state variables that define the internal state.

(従来技術) 従来の内燃機関のアイドル回転速度制御方法としては、
例えば第1図に示すようなものがある。
(Prior art) As a conventional method for controlling the idle rotation speed of an internal combustion engine,
For example, there is one shown in FIG.

第1図にお(・て、アイドル回転速度制御用のAACバ
ルブ1は、VCMバルブ2の制御ソレノイド3の駆動パ
ルス幅PAをデユーティ制御することによってリフト量
が変わり、スロソトルノくルブ4の);イパス5を通過
するバイパス空ネ量が変化して、アイドル回転速度が制
御される。
In FIG. 1, the lift amount of the AAC valve 1 for idle rotational speed control is changed by duty-controlling the drive pulse width PA of the control solenoid 3 of the VCM valve 2, and the amount of lift of the throttle valve 4 is changed; By changing the amount of bypass air passing through the air path 5, the idle rotation speed is controlled.

コントロールユニット6は、スロットルバルブスイッチ
(アイドル回転速度とも言う。)7によるアイドル(I
DLE)信号、ニュートラルスイッチ8による二一一ト
ラル(NEUT)信号、車速センサ9による車速(VS
、P)信号などによって機関がアイドル状態にあること
を検知すると、水温センサ10による冷却水温度(Tw
)K応じた1次元テーブルルックアップによって、アイ
ドル回転速度の基本目標値を算出する。そして、エアコ
ンスイッチ11にヨルエアコン(A/C)M号、二一一
トラル(NEUT)信号、バッテリ電圧(VB)信号な
どに応じた補正を行なって最終的に算出されたアイドル
回転速度の目標値Nrに対し、機関の実際のアイドル回
転速度Nとその目標値Nrとの偏差SAが小さくなるよ
うに制御ソレノイド3の駆動パルス幅P^を比例・積分
(PI)のデユーティ制御をして、目標アイドル回転速
度Nrにフィードバック制御する。
The control unit 6 controls the idle (I) speed by controlling the throttle valve switch (also called idle rotation speed) 7.
DLE) signal, 211 tral (NEUT) signal from the neutral switch 8, vehicle speed (VS
, P) When it is detected that the engine is in an idle state by a signal etc., the cooling water temperature (Tw
) Calculate the basic target value of the idle rotation speed by looking up a one-dimensional table according to K. Then, the target idle rotation speed is finally calculated by correcting the air conditioner switch 11 according to the air conditioner (A/C) M number, the 211 tral (NEUT) signal, the battery voltage (VB) signal, etc. With respect to the value Nr, the drive pulse width P^ of the control solenoid 3 is controlled by proportional-integral (PI) duty so that the deviation SA between the actual idle rotation speed N of the engine and its target value Nr becomes small. Feedback control is performed to the target idle rotation speed Nr.

以上の制御方法を流れ図で示したのが第2図である。FIG. 2 shows a flowchart of the above control method.

しかしながら、このような従来の内燃機関のアイドル回
転速度制御方法にあっては、機関、アクチュエータおよ
びセンサの動特性を効果的に用いたPI制御を行なって
いる訳ではな(、さらには、制御手法としてのPI制御
は多入出力システムに対する制御には不向きなものとな
って(・たため、機関が他の運転状態からアイドル状態
に入る時、またはアイドル状態から出る時、さらには種
々の負荷外乱が加わった直後等の、機関がダイナミック
な振舞いを呈する時には、制御追従性すなわち過渡応答
が悪いという問題があった。また、他の制御入力を加え
て制御の自由度を上げ、制御性を高めようとする時には
、PI制御の手法では適用が難しいという問題があった
However, such conventional idle speed control methods for internal combustion engines do not perform PI control that effectively uses the dynamic characteristics of the engine, actuators, and sensors (and furthermore, the control method PI control is not suitable for controlling multi-input/output systems. When the engine exhibits dynamic behavior, such as immediately after a motor is applied, there is a problem of poor control followability, or transient response.In addition, it is necessary to increase the degree of freedom in control by adding other control inputs to improve controllability. When doing so, there was a problem in that it was difficult to apply the PI control method.

そして特に、アイドル回転速度制御を行なうかどうかの
判定をニュートラルスイッチ8、スロットルバルブスイ
ッチ7および車速で行なうと、高速でギヤを二一一トラ
ルにした場合にエンストし易(、また二一一トラルスイ
ノチ8での判定を除いた場合でも、エンジンブレーキで
減速してクラッチを切った場合にエンストし易いという
問題点があった。
In particular, if the determination of whether or not to perform idle speed control is made based on the neutral switch 8, throttle valve switch 7, and vehicle speed, it is easy to stall the engine when the gear is set to 211 tral at high speed (and the 211 tral Even when excluding the judgment in 8, there was a problem in that the engine was likely to stall if the engine brake was used to decelerate and the clutch was released.

(発明の目的) この発明は、このような従来の問題点に着目してなされ
たもので、機関が他の運転状態からアイドル状態へ入る
時、またはアイドル状態から出る時、さらには負荷外乱
が加わった直後等の、機関がダイナミックな振舞いを呈
する時の制御追従性すなわち過渡応答を最適にし、さら
に、多数の制御入力変数を加えて制御自由度を上げ、制
御性を高めることを容易にし、もってより安定なアイド
ル回転速度制御を行なうことを目的とする。そして特に
、アイドル回転速度制御を行なうかどうかの判定を適切
に行なうようにして、高速運転でのギヤ抜差時やエンジ
ンブレーキでの減速時の耐ストール性を増すことを目的
とする。
(Object of the Invention) The present invention has been made by focusing on the conventional problems as described above. It optimizes the control followability, that is, the transient response when the engine exhibits dynamic behavior, such as immediately after engine application, and also increases the degree of control freedom by adding a large number of control input variables, making it easy to improve controllability. The purpose of this is to perform more stable idle rotation speed control. In particular, it is an object of the present invention to appropriately determine whether or not to perform idle rotational speed control, thereby increasing stall resistance when a gear is inserted or removed during high-speed operation or during deceleration by engine braking.

(発明の構成〕 そこでこの発明は、内燃機関、アクチュエータおよびセ
ンサの動特性をモデル化したものをマイクロコンピュー
タ等からなるコントローラに記憶しておき、バイパス空
気量もしくは相当量、点火時期、燃料供給量もしくは相
当量および排気還流(EGR,)量もしくは相当量から
選択されるいずれか1つまたは任意の2つ以上の組合せ
を制御入力とし、かつアイドル回転速度を制御出方とし
、制御入力と制御出力から、ダイナミツクモチルである
内燃機関等の内部状態を代衣する状態変数量を推定し、
その推定値、およびアイドル回転速度の目標値と実際値
の偏差の積分値とを用いて制御入力値を決定し、内燃機
関のアイドル回転速度を目標値にフィードバンク制御す
ることを特徴とする。
(Structure of the Invention) Therefore, the present invention stores models of the dynamic characteristics of the internal combustion engine, actuators, and sensors in a controller consisting of a microcomputer, etc. or an equivalent amount and an exhaust gas recirculation (EGR) amount or an equivalent amount, or any combination of two or more of them as a control input, and the idle rotation speed as a control output method, and the control input and the control output. From this, we estimate the state variable quantity that substitutes for the internal state of the internal combustion engine, etc., which is a dynamic model,
The control input value is determined using the estimated value and the integral value of the deviation between the target value and the actual value of the idle rotation speed, and the idle rotation speed of the internal combustion engine is subjected to feedbank control to the target value.

この制御手法は、従来一般的なPID制御に代わり、多
数の入出力変数を総合的に制御する多変数制御の手法を
用いるものである。
This control method uses a multivariable control method that comprehensively controls a large number of input and output variables in place of the conventional general PID control.

そして特に、アイドル回転速度制御を行なうかどうかの
判定をスロノトルバルブスインチ、および車速と機関回
転速度との比に注目して行なうことを特徴とするもので
ある。
Particularly, the present invention is characterized in that the determination as to whether or not to perform idle rotational speed control is made by paying attention to the throttle valve blink and the ratio between the vehicle speed and the engine rotational speed.

(実施例) 以下、この発明の実施例を図面を参照して説明する。(Example) Embodiments of the present invention will be described below with reference to the drawings.

第3図は、この発明による内燃機関のアイドル回転速度
制御方法の一実施例を実現する装置の構成図である。
FIG. 3 is a block diagram of an apparatus for implementing an embodiment of the method for controlling the idle rotation speed of an internal combustion engine according to the present invention.

同図において、12は制御対象である内燃機関で、アイ
ドル回転速度制御の他、空燃比フィードバック制御を含
む燃料噴射制御その他を行なっている。
In the figure, reference numeral 12 denotes an internal combustion engine to be controlled, which performs not only idle rotational speed control but also fuel injection control including air-fuel ratio feedback control.

制御対象12の制御出力をアイドル回転速度とした場合
、制御入力としては、空気量、点火時期、燃料供給量お
よび排気還流蓋かも選択されたいずれか1つまたは任意
の2つ以上の組合せをとり得る。
When the control output of the controlled object 12 is the idle rotation speed, the control input may be any one of the air amount, ignition timing, fuel supply amount, and exhaust gas recirculation lid, or any combination of two or more of them. obtain.

本実施例では、2制御入力として、アイドル時のバイパ
ス空気量を調整するためのVCMバルブ2の制御ソレノ
イド3(第1図)を駆動するパルス幅PA(すなわちバ
イパス空気量に相当する量)と点火時期ITとをとる。
In this embodiment, the two control inputs are the pulse width PA (that is, the amount corresponding to the amount of bypass air) that drives the control solenoid 3 (Fig. 1) of the VCM valve 2 for adjusting the amount of bypass air at idle; Take the ignition timing IT.

制御出力はアイドル回転速度Nで、■出力である。The control output is the idle rotation speed N, and is an output.

13は、制御対象である内燃機関12のダイナミノクモ
チルを記憶していて、上記3つの制御入出力情報PA、
IT、Nから機関12のダイナミックな内部状態を推定
する状態観測器(オブザーバ)であり、内部状態を代表
する状態変数量x(=x1.1=1.2゜・・・l”o
例えば4つの量XI+ X2+ X3+ ×4のベクト
ル表示)の推定値x(=Xl、■=1.2.・・・、n
)を計算する。
13 stores the dynamics of the internal combustion engine 12 which is the controlled object, and stores the above three control input/output information PA,
It is a state observer that estimates the dynamic internal state of the engine 12 from IT and N, and the state variable quantity x (=x1.1=1.2°...l”o
For example, the estimated value x (=Xl, ■=1.2..., n
).

状態観測器13は制御対象である機関12をシミュレー
ションするもので、ダイナミックな内部状態を状態変数
ixで代表する。制御対象である機関12の内部状態を
表わす状態変数としては、具体的には例えばインテーク
マニホールドの給体圧や吸入負圧、実際にシリンダに吸
入された空気量、燃焼の動的挙動、機関トルク等が挙げ
られる。これらの値をセンサにより検出できれば、その
検出値を用いることによって動的な振舞いを把握し、制
御をより精密に行なうことができる。しかしながら現時
点ではそれらの値を検出できる実用的なセンサはあまり
存在しな見・。そこで機関12の内部状態を状態変数量
Xで代表させるが、但し状態変数量Xは実際の内部状態
を表わす種々の物理量に対応させる必要はな(、全体と
して機関12をシミュレーションさせるものである。状
態変数量Xの次数nは、nが大きい程シミュレーション
カ1精確になるが、反面計算が複雑になる。そこでモデ
ルとしては低次元近似されたものを使用し、近似誤差ま
たは機関個体差による誤差を積分(I)動作で吸収する
。この発明における2人力1出力の場合には、n = 
4程度が適当である。
The state observer 13 simulates the engine 12 that is the object of control, and the dynamic internal state is represented by a state variable ix. Specifically, state variables representing the internal state of the engine 12 to be controlled include, for example, intake manifold supply pressure, suction negative pressure, amount of air actually taken into the cylinder, dynamic behavior of combustion, and engine torque. etc. If these values can be detected by a sensor, dynamic behavior can be grasped by using the detected values, and control can be performed more precisely. However, at present there are not many practical sensors that can detect these values. Therefore, the internal state of the engine 12 is represented by the state variable quantity X, but the state variable quantity X does not need to correspond to various physical quantities representing the actual internal state (the engine 12 as a whole is simulated. As for the order n of the state variable quantity X, the larger n is, the more accurate the simulation becomes, but the computation becomes more complicated. Therefore, a low-dimensional approximation is used as a model, and errors due to approximation errors or individual engine differences are avoided. is absorbed by integral (I) operation.In the case of two people and one output in this invention, n =
A value of about 4 is appropriate.

第3図において、14は積分動作とゲインブロックで、
第4図に詳細な示すよ5に、機関回転速度の指定された
目標(if Nrと実際値Nとの偏差SAを積分した量
および状態観測器13で計算された状態変数量の推定値
Xから、2つの制御入力PAとITの値を計算する。そ
して、下記状態観測器13と積分動作とゲインブロック
14とでコントローラを構成する。
In Fig. 3, 14 is an integral operation and gain block,
As shown in detail in FIG. The values of the two control inputs PA and IT are calculated from the following.A controller is constituted by the state observer 13, integral operation, and gain block 14 described below.

次に作用を説明する。Next, the action will be explained.

制御対象である機関12は2人力1出カシステムで、こ
の入出力間の回転同期サンプル値系のある基準設定値近
辺で求められた線彫近似された伝達関数行列T (z)
から、機関12のダイナミックな内部状態を推定するこ
とが可能である。その1つの手法として状態観測器13
がある。アイドル回転速度近辺の運転条件で、機関12
の伝達関数行列T(z)が実験的に求まり、 となる。但し、2は入出力信号のサンプル値の2−変換
を示し、TI(Z)とT2(zlは例えば202次伝達
関数である。
The engine 12 to be controlled is a two-man power, one-output system, and the transfer function matrix T (z) is a line-shaped approximation obtained near a certain reference setting value of the rotation synchronized sample value system between the input and output.
From this, it is possible to estimate the dynamic internal state of the engine 12. As one method, the state observation device 13
There is. Under operating conditions near idle speed, engine 12
The transfer function matrix T(z) of T(z) is found experimentally, and becomes the following. However, 2 indicates 2-conversion of the sample value of the input/output signal, and TI(Z) and T2(zl are, for example, 202nd-order transfer functions).

入力、出力および伝達関数T+(z)、T2(Z)の関
係を示す機関12のモデル構造を第5図に示す。但し、
入出力はそれぞれ基準設定値からのズレδPA+δIT
FIG. 5 shows a model structure of the engine 12 showing the relationship between input, output, and transfer functions T+(z) and T2(Z). however,
The input and output are each deviation from the reference setting value δPA + δIT
.

δNを用いている。δN is used.

この伝達関数行列T(zlから、次のように状態観測器
13を構成することができる。
From this transfer function matrix T(zl), the state observer 13 can be configured as follows.

先ず、T(z)から機関12の動的な振舞いを記述する
状態変数モデル x(n) −Ax(n−1) 十Bu(n−1)   
  (21y(n−1)=Cx(n−1)      
 (31を導(。ここで、容量のカッコ内のnは現時点
を、n−1は1つ前のサンプル時点を表わす。u(n−
1)は制御入力ベクトルで、ある基準設定値からの線形
近似が成り立つ範囲内での摂動分を表わす、制御ソレノ
イド3(第1図)の駆動パルス幅δPA(nl)と点火
時期δITを要素とする。すなわち、また、y(n−1
)は制御出力で、制御入力ベクトルと同様に、ある基準
回転速度Na(例えば650rpm)からの摂動分を表
わすδN(n−1)を要素とする。
First, from T(z), a state variable model x(n) −Ax(n-1) 10Bu(n-1) describing the dynamic behavior of the engine 12 is created.
(21y(n-1)=Cx(n-1)
(31 is derived (.Here, n in the parentheses of the capacity represents the current time, n-1 represents the previous sample time.
1) is a control input vector, which includes the driving pulse width δPA (nl) of the control solenoid 3 (Fig. 1) and the ignition timing δIT, which represent the perturbation within a range where linear approximation from a certain reference setting value holds. do. That is, also y(n-1
) is a control output, which, like the control input vector, has an element of δN(n-1) representing a perturbation from a certain reference rotational speed Na (for example, 650 rpm).

すなわち、 y(n−1)=δN(n−1)(51 X(・)は状態変数ベクトルであり、行列A、B、Cは
伝達関数行列T(z)の係数から決まる定数行列である
That is, y (n-1) = δN (n-1) (51 .

ここで、次のようなアルゴリズムを持つ状態観測器を構
成する。
Here, we configure a state observer with the following algorithm.

x(n)−(A−GC)x(n−1)+Bu(n−1)
+Gy(n−1)  (6)ここに、Gは任意に与えら
れる行列で、X(・)は機関12の内部状態変数X(・
)の推定値である。(2)(31(6)式より変形する
と、 (x(n)−x(nl〕=(A−GC)〔x(n−1)
−x(n−1)]   f力となり、行列(A−GC)
の固有値が単位円内にあるようにGを選べば、 n→太で   x(0)→X(n)      (8)
となり、内部状態変数量x(n)を入力U(・)と出力
y(・)から推定することができる。また、行列Gを適
当に選び、行列(A−GC)の固有値を全て零にするこ
とも可能で、この時状態観測器13は有限整定状態観測
器となる。
x(n)-(A-GC)x(n-1)+Bu(n-1)
+Gy(n-1) (6) Here, G is an arbitrarily given matrix, and X(・) is the internal state variable X(・
) is the estimated value. (2) (31 Transforming from equation (6), (x(n)-x(nl)=(A-GC)[x(n-1)
-x(n-1)] f force, matrix (A-GC)
If we choose G so that the eigenvalue of is within the unit circle, then n→thick, x(0)→X(n) (8)
Therefore, the internal state variable quantity x(n) can be estimated from the input U(·) and the output y(·). It is also possible to appropriately select the matrix G and make all the eigenvalues of the matrix (A-GC) zero, in which case the state observer 13 becomes a finitely stable state observer.

このようにして推定された状態変数量X(・)と、目標
回転速度Nrと現在の実際の回転速度N(・)との偏差
5A=(N、−N(・))の情報を用いて、制御入力で
ある制御ソレノイド3の駆動パルス幅の基準設定値(P
A)&からの線形近似が成り立つ範囲内での増量分δP
A(・)と、点火時期の基準設定値■Taからの線形近
似が成り立つ範囲内での増量分δIT(・)を決定し、
機関のアイドル回転速度Nの最適レギュレータ制御を行
なう。レギュレータ制御とは、アイドル回転速度Nを一
定値である目標回転速度Nrに合致するように制御する
定値制御を意味する。
Using the state variable amount X(・) estimated in this way and the information of the deviation 5A=(N, -N(・)) between the target rotational speed Nr and the current actual rotational speed N(・), , the reference setting value (P
A) Increase amount δP within the range where linear approximation from &
Determine the amount of increase δIT(・) within the range where linear approximation holds from A(・) and the reference setting value of ignition timing ■Ta,
Optimal regulator control of engine idle speed N is performed. Regulator control means constant value control that controls the idle rotation speed N to match the target rotation speed Nr, which is a constant value.

なおこの発明では、前述したように実験的に求めたモデ
ルが低次元化された近似モデルであるため、その近似誤
差を吸収するための積分(I)動作を付加して(・るが
、ここでは積分動作を含めて最適レギュレータ制御を行
なう。
In this invention, since the experimentally obtained model is a reduced-dimensional approximate model as described above, an integral (I) operation is added to absorb the approximation error. Now, perform optimal regulator control including integral operation.

この発明の制御対象である一関は、前述したように2人
力l出カシステムであり、これを最適にレギュレータ制
御するものであるが、一般的な多変数システムの最適レ
ギュレータ制御アルゴリズムは、例えば、吉田勝久著「
線形システム制御理論」(昭51年)昭晃堂その他に説
明されているので、ここでは詳細な説明は省略する。結
果のみを記述すると、いま、 δu(n)=u(n)−u(n−1)       (
9)δe fn) = Nr  N (n)     
     QO)とし、評価関数Jを、 J= Σ〔δe(k12+δu蜀R!3u(k)〕’ 
 (I I)k=0 とする。ここでRは重みパラメータ行列、【は転置を示
す。kは制御開始時点を0とするサンプル回数で、Rを
対角行列とすると(11)式の右辺第2項は(9)式の
2乗を表わす。また、00式の右辺第2項を(9)式の
ような制御入力の差分の2次形式としているが、これは
第4図のように積分動作を付加したためである。
As mentioned above, the Ichinoseki that is the control target of this invention is a two-man power output system, which is optimally controlled by a regulator, but the optimal regulator control algorithm for a general multivariable system is, for example, Written by Katsuhisa Yoshida “
Since it is explained in "Linear System Control Theory" (1972) by Shokodo et al., detailed explanation will be omitted here. Describing only the results, now δu(n)=u(n)-u(n-1) (
9) δe fn) = Nr N (n)
QO), and the evaluation function J is J= Σ[δe(k12+δu ShuR!3u(k)]'
(I I) Let k=0. Here, R is a weight parameter matrix and [denotes transposition. If k is the number of samples with the control start time being 0, and R is a diagonal matrix, the second term on the right side of equation (11) represents the square of equation (9). Furthermore, the second term on the right side of equation 00 is a quadratic form of the difference in control inputs as shown in equation (9), but this is because an integral action is added as shown in FIG.

制御人力u (klは、 となる。02式で、 K= −(R十B PB)B PA         
   (13)とおくと、Kは最適ゲイン行列である。
Control human power u (kl is as follows. In formula 02, K= - (R 1 B PB) B PA
(13), K is the optimal gain matrix.

また(121式であり、Pは、 のリカノテイ(Riccati )方程式の解である。Also, (formula 121, P is This is the solution of the Riccati equation.

(11)式の評価関数Jの意味は、制御人力U(・)の
動きを制約しつつ、制御出力y(・)であるアイドル回
転速度Nの目標値Nrからの偏差SA (回転変動)を
最小にしようと意図したもので、その制約の重み付けは
重みパラメータ行列Rで変えることができる。従って、
適当なRを選択し、アイドル時の機関のダイナミックモ
デル(状態変数モデル)を用い、(16)式を解いたP
を用いて計算したa′5式の最適ゲイン行列Kをマイク
ロコンビーータに記憶し、アイドル回転速度の目標値N
rと実際値Nの偏差SAの積分値および推定された状態
変数量x(k)から、02式によって最適制御入力値u
 (k)を簡単に決定することができる。また前述した
ように、機関12のダイナミックな状態変数量の推定値
x(k)を求めるには、行列へ、B、C,Gの値をマイ
クロコンビーータに記憶しておき、(6)式により計算
すればよい。
The meaning of the evaluation function J in equation (11) is to limit the deviation SA (rotation fluctuation) of the idle rotation speed N, which is the control output y (・), from the target value Nr while constraining the movement of the control human power U (・). The weighting of the constraints can be varied with the weight parameter matrix R. Therefore,
Select an appropriate R, use a dynamic model (state variable model) of the engine at idle, and solve equation (16).P
The optimum gain matrix K of formula a'5 calculated using
From the integral value of the deviation SA between r and the actual value N and the estimated state variable quantity x(k), the optimal control input value u is determined by formula 02.
(k) can be easily determined. Further, as mentioned above, in order to obtain the estimated value x(k) of the dynamic state variable quantity of the engine 12, the values of B, C, and G are stored in the matrix in the microconbeater, and (6) It can be calculated using the formula.

このような内燃機関のアイドル回転速度制御方法におい
て、前述したように、特に、アイドル回転速度制御を行
なうがどうかの判定(アイドル判定)を、従来のように
二一一トラルスイッチ8、スロットルバルブスインチア
および車速で行なうと、高速回転中にギヤを二一一トラ
ルにした場合にエンストし易いという問題があり、二一
一トラルスイッチによる判定を除いても、エンジンブレ
ーキによる減速時にアイドル回転速度制御を行なってし
まうため、減速途中で実際の機関回転速度Nが目標回転
速度Nrより高い状態で、バイパス空気量を減少させる
方向に制御が行われる。この状態でクラッチを切ると、
機関回転速度が低下し、エンストしてしまうという問題
がある。
In such an idle rotation speed control method for an internal combustion engine, as described above, in particular, the determination of whether or not to perform idle rotation speed control (idle determination) is performed using the conventional 211 tral switch 8 and throttle valve switch. If this is done at the engine brake and vehicle speed, there is a problem that the engine is likely to stall if the gear is set to 211 tral during high speed rotation. Therefore, control is performed in the direction of reducing the amount of bypass air while the actual engine rotational speed N is higher than the target rotational speed Nr during deceleration. If you release the clutch in this condition,
There is a problem that the engine speed decreases and the engine stalls.

従ってこの発明では、アイドル判定をニュートラルスイ
ッチ8では行わず、スロットルバルブスイッチ7、およ
び車速と機関回転速度との比に注目して行なうこととし
た。
Therefore, in the present invention, the idle determination is not performed using the neutral switch 8, but is performed by paying attention to the throttle valve switch 7 and the ratio between the vehicle speed and the engine rotation speed.

以上のアイドル回転速度制御の手順を示したのが第6図
である。手順を説明すると、ステップ(9)では、エア
コンのオン・オフ状態、冷却水温度TrwO値等により
アイドル回転速度の目標値Nrを決める。ステノブ31
テは、スロットルバルブスイッチ7が入っているかどう
かを判定し、入っていなければリターンし、入っていれ
ばステップ32に進む。
FIG. 6 shows the procedure of the above idle rotation speed control. To explain the procedure, in step (9), the target value Nr of the idle rotation speed is determined based on the on/off state of the air conditioner, the cooling water temperature TrwO value, etc. Stenobu 31
In step 3, it is determined whether the throttle valve switch 7 is on, and if it is not on, the process returns, and if it is on, the process proceeds to step 32.

ステップ32では、車速と機関回転速度との比から、ク
ラッチがつながっているかどうかを判定し、つながって
いればリターンし、つながっていなければステップ33
に進む。ステップ33では、初めてアイドル回転速度制
御に入ったかどうかを判定し、初めて入った場合は、ス
テップあでスムーズにアイドル回転速度制御に入れるよ
うに、実際の機関回転速度Nに応じて状態変数量xIn
 x2+ x3+ x4およびDUNに初期値を与え、
ステップ39へ進む。
In step 32, it is determined whether the clutch is engaged based on the ratio of the vehicle speed and the engine rotational speed. If it is engaged, the process returns; if it is not engaged, step 33 is carried out.
Proceed to. In step 33, it is determined whether the idle rotation speed control is entered for the first time, and if it is the first time, the state variable quantity
Give initial values to x2+ x3+ x4 and DUN,
Proceed to step 39.

ステップ33でアイドル回転速度制御が初めて入ったも
のではない場合は、ステハブ35でアイドル回転速度の
目標値Nrと実際値Nの備差SAを計算する。ステップ
36では、制御を始めてから前の周期までの回転速度の
偏差SAを加算していて、結果をDUNというレジスタ
に移す。ステップ37では、回転速度の実際値Nの基準
設定値Na (例えば650rpm )からのズレδN
を計算する。ステップあでは、前の制御周期で推定され
た機関のダイナミックな内部状態を表わす状態変数量x
1〜X3(前回計算値)と、計算された制御入力値δP
AおよびδITと、さらに制御出力値であるδNとを重
み付は加算して、各状態変数量X、〜x4を計算する。
If it is determined in step 33 that this is not the first time that idle rotational speed control has been applied, the step hub 35 calculates the difference SA between the target value Nr and the actual value N of the idle rotational speed. In step 36, the rotational speed deviation SA from the start of control to the previous cycle is added, and the result is transferred to a register called DUN. In step 37, the deviation δN of the actual value N of the rotational speed from the standard set value Na (for example, 650 rpm) is determined.
Calculate. In step A, the state variable quantity x representing the dynamic internal state of the engine estimated in the previous control cycle is
1 to X3 (previously calculated value) and the calculated control input value δP
A, δIT, and δN, which is the control output value, are weighted and added to calculate each state variable amount X, to x4.

但し、(6)式の行列(A−GC月末、 の形で、有限整定オブザーノ;を形成した例である。However, the matrix of equation (6) (A-GC month end, This is an example of forming a finitely settled observer in the form.

なお、(A、B、C)は可観測正準形を用(・て(・る
Note that (A, B, C) uses the observable canonical form (・te(・ru).

ステップ39では、推定された機関のダイナミックな内
部状態変数量X、〜X4とDUNに最適ゲインにの要素
に、jを乗じて加算し、基準設定値(PA)aおよびl
7Ll’、に対し制御入力値をどれだけ増量する力・を
計算する。
In step 39, the estimated dynamic internal state variables of the engine X, ~
7Ll', calculate how much force to increase the control input value.

第6図の係数E)ij + gi+ kij等は、予め
求めておし・テマイクロコンビーータ等に記憶してお(
The coefficients E) ij + gi + kij, etc. in Fig. 6 are calculated in advance and stored in a microcomputer, etc. (
.

以上の手順で、アイドル回転速度が一定の状態における
種々の外乱に対する過渡応答と、アイドル回転速度の目
標値を変更した場合の過渡応答を実験した結果を、従来
のPI制御とこの発明による多変数制御とで比較したの
が第7図ないし第10図である。
Using the above procedure, we compared the results of experiments on the transient response to various disturbances when the idle rotation speed is constant and the transient response when the target value of the idle rotation speed is changed. Figures 7 to 10 show a comparison with the control.

第7図はクラッチ接続時(10点で半クラツチ接続、但
しブレーキを踏んでいる)のアイドル回転速度Nの過渡
応答を示し、(A)は従来のPI制御、(Blはこの発
明の多変数制御の場合である。第8図をまクラッチ遮断
時(10点で遮断)の過渡応答を示し、囚は従来方法、
(B)はこの発明の方法の場合である。
FIG. 7 shows the transient response of the idle rotation speed N when the clutch is engaged (the clutch is half engaged at 10 points, but the brake is depressed), (A) is the conventional PI control, (Bl is the multivariable variable of this invention). In the case of control, Figure 8 shows the transient response when the clutch is disengaged (disengaged at 10 points).
(B) is the case of the method of this invention.

第9図はエアコンをオンし、目標アイドル回転速度Nr
を80Orpmに移行した場合、およびエアコンをオフ
し、目標アイドル回転速度Nrを650rpmに戻した
場合の過渡応答を示し、(A)は従来方法、(Blはこ
の発明の方法の場合である。第10図は無負荷高回転状
態から目標値Nr ” 650rpmにコーステイング
する場合の過渡応答を示し、(A)は従来方法、CB+
はこの発明の方法の場合である。第7図ないし第10図
から明らかなように、いずれの場合もこの発明の方法に
よる方が、過渡応答性が大幅に改善されていることがわ
かる。なお、第7図(A)ではアイドル回転速度が目標
値Nrに整冗しない。
Figure 9 shows the target idle rotation speed Nr when the air conditioner is turned on.
80 rpm, and when the air conditioner is turned off and the target idle rotation speed Nr is returned to 650 rpm, (A) is the conventional method, (Bl is the method of the present invention. Figure 10 shows the transient response when coasting from a no-load high speed state to the target value Nr'' 650rpm, and (A) shows the conventional method, CB+
is the case with the method of this invention. As is clear from FIGS. 7 to 10, it can be seen that in all cases, the method of the present invention significantly improves the transient response. In addition, in FIG. 7(A), the idle rotation speed is not aligned with the target value Nr.

第11図は高速回転中にギヤをニュートラルにした場合
の実験結果で、(Alは二一一トラルスインチをアイド
ル判定に入れた従来方法、(B)はニュートラルスイッ
チをアイドル判定に入れないこの発明の方法をそれぞれ
示す。(Blの方が目標回転速度Nr ”” 650r
pmからの落込みが小さいことがわかる。
Figure 11 shows the experimental results when the gear is set to neutral during high-speed rotation. Each method is shown below. (Bl has a target rotational speed Nr "" 650r
It can be seen that the drop from pm is small.

なお、(B)では機関回転速度が1.1.0Orpmを
通過する点t、で、第6図のステップ34で説明した初
期値をセントしている。
In addition, in (B), the initial value explained in step 34 of FIG. 6 is set at the point t where the engine rotational speed passes 1.1.0 Orpm.

第12図はエンジンブレーキで減速後にクラッチを切っ
た場合の実験結果であり、(Alは従来方法、(Blは
この発明の方法である。期間t。〜t1中にエンジンブ
レーキで減速し、時刻【1でクラッチを切っている。(
B)の方が目標回転速度Nrからの落込みが小さいこと
がわかる。
Figure 12 shows the experimental results when the clutch is disengaged after deceleration by engine braking, (Al is the conventional method, (Bl is the method of the present invention). [The clutch is disengaged at 1. (
It can be seen that B) has a smaller drop from the target rotational speed Nr.

前述したように、この発明における内燃機関の制御を力
をアイドル回転速度とした時に、制御入力としては、空
気量(または相当量)、点火時期、燃料供給量(または
相当量)および排気還流量(または相当量)から選択さ
れたいずれか1つまたは任意の2つ以上の組合せを用い
ることができ、上述の実施例では、バイノくス空気量の
相当量であるV C’Mバルブの制御ソレノイドの駆動
ノ(117幅と点火時期とを制御入力とする場合につい
て説明した。
As mentioned above, when controlling the internal combustion engine in this invention with the power being the idle speed, the control inputs include the air amount (or equivalent amount), ignition timing, fuel supply amount (or equivalent amount), and exhaust gas recirculation amount. (or an equivalent amount) can be used, and in the above embodiment, the control of the V C'M valve, which is an equivalent amount of the binox air amount, can be used. The case where the solenoid drive width (117 width) and ignition timing are used as control inputs has been described.

(発明の効果) 以上説明したように、この発明の内燃機関のアイドル回
転速度制御方法によれば、内燃機関のダイナミックなモ
デルに基づく多変数制御手法を適用してアイドル回転速
度制御を行なし・、しかも内燃機関のダイナミックな状
態を推定する手順を付加し、なおかつ、状態観測器内の
機関モデルを低次元化したものを用いて計算時間を短縮
し、その近似誤差分は積分動作で吸収するようにしたた
め、アイドル状態で問題となる失火外乱や負荷外乱なと
の外乱に対する制御過渡応答を最適にでき、しかも制御
自由度を上げて制御性を高めるために、多変数制御入力
を加えて制御することも容易であり、より安定なアイド
ル回転速度制御が実現できるという効果が得られる。
(Effects of the Invention) As explained above, according to the method for controlling the idle rotation speed of an internal combustion engine of the present invention, the idle rotation speed is controlled by applying a multivariable control method based on a dynamic model of the internal combustion engine. Moreover, it adds a step to estimate the dynamic state of the internal combustion engine, reduces calculation time by using a reduced-dimensional version of the engine model in the state observer, and absorbs the approximation error by integral operation. As a result, it is possible to optimize the control transient response to disturbances such as misfire disturbances and load disturbances, which are problems in the idling state.Moreover, in order to increase the degree of control freedom and improve controllability, multivariable control inputs are added to the control. It is also easy to do so, and the effect that more stable idle rotation speed control can be achieved can be obtained.

そして特に、アイドル判定を、アイドルスイッチおよび
車速と機関回転速度との比に注目して行なう構成とした
ため、高速回転中のギヤ抜き時にエンストしにくく、か
つエンジンブレーキで減速した後の耐ストール性が向上
し、より安定なアイドル運転を実現できるという効果が
得られる。
In particular, since the idle judgment is performed by focusing on the idle switch and the ratio between the vehicle speed and the engine rotation speed, the engine is less likely to stall when out of gear during high-speed rotation, and the stall resistance after decelerating with engine braking is improved. The effect is that it is possible to realize more stable idling operation.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は従来の内燃機関のアイドル回転速度制御装置の
構成図、第2図は従来のアイドル回転速度制御方法を示
すフローチャート、第3図はこの発明による内燃機関の
アイドル回転速度制御方法を実現する制御装置の構成図
、第4図は第3図の積分動作とゲインブロックの詳細な
構成図、第5図は第3図の制御入出力と機関の関係を示
すブロック図、第6図はこの発明による制御方法を説明
するフローチャート、第7図(A)(B)はクラッチ接
続時の過渡応答の実験結果を示す図、第8図(AXBI
はクラッチ遮断時の過渡応答の実験結果を示す図、第9
図(A)(B+はエアコンのオンオフ時の過渡応答の実
験結果を示す図、第10図(A)(B)はコーステイン
グ時の過渡応答の実験結果を示す図、第11図(A)(
B)はアイドル判定に二一一トラルスインチを用いた場
合と用いない場合の高速ギヤ抜き時の実験結果を示す図
、第12図(A)(Blはエンジンブレーキで減速後に
クラッチを切った場合の実験結果を示す図である。 1・・・AACバルブ、  2・・・VCMバルブ、3
・・・制御ソレノイド、  4・・・スロットルバルブ
、5・・・バイパス、7・・・スロットルバルブスイッ
チ、8・・・二一一トラルスイノチ、10・・・水温セ
ンサ、11・・・エアコンスイッチ、12・・・内燃機
関(制御対象)、13・・・状態観測器、14・・・積
分動作とゲインブロック、N・・・アイドル回転速度の
実際値、 Nr・・・アイドル回転速度の目標値、N、・・・アイ
ドル回転速度の基準設定値、SA・・・アイドル回転速
度の目標値と実際値の偏差、PA・・・バイパス空気量
を規定する制御ソレノイドの駆動パルス幅、 IT・・・点火時期、    !(=xi)・・・状態
変数量、X(”Xi)・・・状態変数の推定量。 特許出願人 日産自動車株式会社 特許出願代理人 弁理士  山 本 恵 − 幕3図 竿、4U!J 4 幕7図 (A) g!f藺i(S扛) (8ン 時間f、(Seり 革、8r:!J (A] (8) 時開t(6eり 秦デ図 (A〕 ’B%MIf;(sec) (δン 吟%’T 6 (seす 21E、 10凹 (A) CB)
Fig. 1 is a block diagram of a conventional idle rotation speed control device for an internal combustion engine, Fig. 2 is a flowchart showing a conventional idle rotation speed control method, and Fig. 3 is a realization of the idle rotation speed control method for an internal combustion engine according to the present invention. Fig. 4 is a detailed block diagram of the integral operation and gain block shown in Fig. 3, Fig. 5 is a block diagram showing the relationship between the control input/output and the engine shown in Fig. 3, and Fig. 6 is A flowchart explaining the control method according to the present invention, FIGS.
Figure 9 shows the experimental results of the transient response when the clutch is disengaged.
Figures (A) (B+ are diagrams showing the experimental results of the transient response when turning on and off the air conditioner, Figures 10 (A) and (B) are diagrams showing the experimental results of the transient response during coasting, Figure 11 (A) (
B) is a diagram showing the experimental results during high-speed gear removal with and without using the 211 Torals inch for idle determination, and Figure 12 (A) (Bl is the result when the clutch is disengaged after deceleration with engine braking. It is a diagram showing experimental results. 1... AAC valve, 2... VCM valve, 3
...Control solenoid, 4...Throttle valve, 5...Bypass, 7...Throttle valve switch, 8...211 Toralsinochi, 10...Water temperature sensor, 11...Air conditioner switch, 12... Internal combustion engine (controlled object), 13... Condition observation device, 14... Integral operation and gain block, N... Actual value of idle rotation speed, Nr... Target value of idle rotation speed , N,...Reference setting value of idle rotation speed, SA...Difference between target value and actual value of idle rotation speed, PA...Driving pulse width of control solenoid that defines bypass air amount, IT...・Ignition timing! (=xi)...Amount of state variable, X ("Xi)...Estimated amount of state variable. Patent applicant: Nissan Motor Co., Ltd. Patent application agent Megumi Yamamoto - Act 3 diagram pole, 4U!J 4 Act 7 diagram (A) g!f藺i (S扛) (8n time f, (Sere leather, 8r:!J (A) (8) Time opening t (6e り駦でG) (A) 'B%MIf; (sec) (δngin%'T 6 (sesu21E, 10 concave (A) CB)

Claims (1)

【特許請求の範囲】[Claims] 内燃機関のアイドル時に、アイドル回転速度の目標値N
rと実際値Nの偏差SAに基づいてアイドル回転速度を
フィードバンク制御する方法において、コントローラに
記憶された前記内燃機関のダイナミックモデルに基づき
、前記内燃機関の制御入力値であるスロットルバルブを
バイパスして該内燃機関に供給されるバイパス空気量も
しくは該バイパス空気量に相当する量、該内燃機関の点
火時期、該内燃機関への燃料供給量もしくは該燃料供給
量に相当する量および該内燃機関への排気還流量もしく
は該排気還流量に相当する量から選択されるいずれか1
つまたは任意の2つ以上の組合せと、該内燃機関の制御
出力値であるアイドル回転速度とから、該内燃機関のダ
イナミックな内部状態を代表する適当な次数の状態変数
量xi(i=1.2.・・・、n)を推定し、該推定さ
れた状態変数量x1から前記制御入力値を決定し、さら
に、アイドル回転速度制御を行なうかどうかの判定を、
スロットルバルブスイッチ、車速と機関回転速度との比
および機関回転速度によって行なうことを特徴とする内
燃機関のアイドル回転速度制御方法。
When the internal combustion engine is idling, the target value N of the idle rotation speed
In a method for feedbank controlling an idle rotation speed based on a deviation SA between r and an actual value N, a throttle valve, which is a control input value of the internal combustion engine, is bypassed based on a dynamic model of the internal combustion engine stored in a controller. the amount of bypass air supplied to the internal combustion engine or the amount equivalent to the amount of bypass air; the ignition timing of the internal combustion engine; the amount of fuel supplied to the internal combustion engine or the amount equivalent to the amount of fuel supplied to the internal combustion engine; Any one selected from the exhaust gas recirculation amount or the amount equivalent to the exhaust gas recirculation amount.
or any combination of two or more and the idle rotational speed which is the control output value of the internal combustion engine, a state variable quantity xi (i=1. 2...., n), determine the control input value from the estimated state variable amount x1, and further determine whether to perform idle rotation speed control,
A method for controlling idle rotation speed of an internal combustion engine, characterized in that control is performed using a throttle valve switch, a ratio between vehicle speed and engine rotation speed, and engine rotation speed.
JP58018100A 1983-02-08 1983-02-08 Control method of idling speed for internal-combustion engine Granted JPS59145338A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58018100A JPS59145338A (en) 1983-02-08 1983-02-08 Control method of idling speed for internal-combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58018100A JPS59145338A (en) 1983-02-08 1983-02-08 Control method of idling speed for internal-combustion engine

Publications (2)

Publication Number Publication Date
JPS59145338A true JPS59145338A (en) 1984-08-20
JPS6340928B2 JPS6340928B2 (en) 1988-08-15

Family

ID=11962202

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58018100A Granted JPS59145338A (en) 1983-02-08 1983-02-08 Control method of idling speed for internal-combustion engine

Country Status (1)

Country Link
JP (1) JPS59145338A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63154862A (en) * 1986-12-17 1988-06-28 Mitsubishi Electric Corp Electronic control device for internal combustion engine

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63154862A (en) * 1986-12-17 1988-06-28 Mitsubishi Electric Corp Electronic control device for internal combustion engine

Also Published As

Publication number Publication date
JPS6340928B2 (en) 1988-08-15

Similar Documents

Publication Publication Date Title
JPS6349060B2 (en)
CN1796749B (en) engine control unit
Van Nieuwstadt et al. Coordinated EGR-VGT control for diesel engines: an experimental comparison
JPH01177431A (en) Electronic throttle control device for internal combustion engines
EP3449109B1 (en) Physics-based vehicle turbocharger control techniques
JPS6325174B2 (en)
Gäfvert et al. Control of GDI engines using torque feedback exemplified by simulations
US6947824B1 (en) Engine RPM and torque control transition
JPS59145338A (en) Control method of idling speed for internal-combustion engine
JPS6325177B2 (en)
JPS6330491B2 (en)
JPS6349061B2 (en)
JPS60216045A (en) Intake-air amount controller for internal-combustion engine
JPS597751A (en) Simultaneous control for idle revolution speed and air-fuel ratio in internal-combustion engine
JPS5943943A (en) Control of idle revolution speed for internal-combustion engine
JPS6340929B2 (en)
JPS6343575B2 (en)
JPS6340264B2 (en)
JPS59120753A (en) Method of controlling idling speed of rotation for internal-combustion engine
JPS6328225B2 (en)
JPS6331659B2 (en)
JPH06504603A (en) Method and apparatus for closed-loop control of power provided by an internal combustion engine driving a motor vehicle
Hu et al. Electronic throttle controller design using a triple-step nonlinear method
Cook et al. Automotive control systems
Morelos et al. Fuzzy control strategy for stoichiometric air-fuel mixture in automotive systems