JP2003241803A - 制御装置 - Google Patents

制御装置

Info

Publication number
JP2003241803A
JP2003241803A JP2002039094A JP2002039094A JP2003241803A JP 2003241803 A JP2003241803 A JP 2003241803A JP 2002039094 A JP2002039094 A JP 2002039094A JP 2002039094 A JP2002039094 A JP 2002039094A JP 2003241803 A JP2003241803 A JP 2003241803A
Authority
JP
Japan
Prior art keywords
value
fuel ratio
air
output
sampling
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
JP2002039094A
Other languages
English (en)
Other versions
JP3973922B2 (ja
Inventor
Yuji Yasui
裕司 安井
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.)
Honda Motor Co Ltd
Original Assignee
Honda 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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2002039094A priority Critical patent/JP3973922B2/ja
Priority to DE10306278A priority patent/DE10306278B4/de
Priority to US10/366,670 priority patent/US7124013B2/en
Publication of JP2003241803A publication Critical patent/JP2003241803A/ja
Priority to US11/503,950 priority patent/US7647157B2/en
Application granted granted Critical
Publication of JP3973922B2 publication Critical patent/JP3973922B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/047Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators the criterion being a time optimal performance criterion
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/02Circuit arrangements for generating control signals
    • F02D41/14Introducing closed-loop corrections
    • F02D41/1438Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1477Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the regulation circuit or part of it,(e.g. comparator, PI regulator, output)
    • F02D41/1479Using a comparator with variable reference
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/24Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means
    • F02D41/2406Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means using essentially read only memories
    • F02D41/2425Particular ways of programming the data
    • F02D41/2429Methods of calibrating or learning
    • F02D41/2451Methods of calibrating or learning characterised by what is learned or calibrated
    • F02D41/2454Learning of the air-fuel ratio control
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/02Circuit arrangements for generating control signals
    • F02D41/14Introducing closed-loop corrections
    • F02D41/1401Introducing closed-loop corrections characterised by the control or regulation method
    • F02D2041/1413Controller structures or design
    • F02D2041/1415Controller structures or design using a state feedback or a state space representation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/02Circuit arrangements for generating control signals
    • F02D41/14Introducing closed-loop corrections
    • F02D41/1401Introducing closed-loop corrections characterised by the control or regulation method
    • F02D2041/1433Introducing closed-loop corrections characterised by the control or regulation method using a model or simulation of the system
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D2200/00Input parameters for engine control
    • F02D2200/02Input parameters for engine control the parameters being related to the engine
    • F02D2200/04Engine intake system parameters
    • F02D2200/0404Throttle position
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D2200/00Input parameters for engine control
    • F02D2200/02Input parameters for engine control the parameters being related to the engine
    • F02D2200/04Engine intake system parameters
    • F02D2200/0406Intake manifold pressure
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D2200/00Input parameters for engine control
    • F02D2200/02Input parameters for engine control the parameters being related to the engine
    • F02D2200/04Engine intake system parameters
    • F02D2200/0414Air temperature
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D2200/00Input parameters for engine control
    • F02D2200/70Input parameters for engine control said parameters being related to the vehicle exterior
    • F02D2200/703Atmospheric pressure
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/02Circuit arrangements for generating control signals
    • F02D41/14Introducing closed-loop corrections
    • F02D41/1438Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1439Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the position of the sensor
    • F02D41/1441Plural sensors
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/02Circuit arrangements for generating control signals
    • F02D41/14Introducing closed-loop corrections
    • F02D41/1438Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor
    • F02D41/1444Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases
    • F02D41/1454Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio
    • F02D41/1456Introducing closed-loop corrections using means for determining characteristics of the combustion gases; Sensors therefor characterised by the characteristics of the combustion gases the characteristics being an oxygen content or concentration or the air-fuel ratio with sensor output signal being linear or quasi-linear with the concentration of oxygen
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/24Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means
    • F02D41/2406Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means using essentially read only memories
    • F02D41/2425Particular ways of programming the data
    • F02D41/2429Methods of calibrating or learning
    • F02D41/2477Methods of calibrating or learning characterised by the method used for learning

Abstract

(57)【要約】 【課題】 位相遅れやむだ時間などが比較的大きい制御
対象を制御する場合において、制御対象の入出力間での
制御タイミングのずれの解消、および制御精度の向上を
いずれも達成することができる制御装置を提供する。 【解決手段】 制御装置1のECU2は、状態予測器2
2およびDSMコントローラ24を備える。状態予測器
22は、予測アルゴリズムに基づき、出力偏差VO2の
予測値PREVO2を所定の算出周期ΔTkで算出し、
DSMコントローラ24は、ΣΔ変調アルゴリズムに基
づき、算出された予測値PREVO2に応じて、目標空
燃比KCMDを、予測値PREVO2の算出周期ΔTk
よりも短い算出周期ΔTmで算出する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、制御対象の出力を
制御するために、制御対象への制御入力をΔΣ変調アル
ゴリズムなどに基づいて算出する制御装置に関する。
【0002】
【従来の技術】従来、この種の制御装置として、例えば
特開2001−154704号公報に記載されたものが
知られている。この制御装置は、制御対象の出力を検出
し、その検出結果をアナログ量の検出信号として出力す
る検出手段と、上位装置から入力されるアナログ量の目
標値と検出信号との偏差を演算する偏差演算手段と、演
算された偏差を1ビットデジタル信号に変換する変換手
段と、変換手段からの1ビットデジタル信号を補償した
後、入力信号として制御対象に入力する補償手段とを備
えている(同公報の図6参照)。
【0003】この制御装置では、偏差演算手段により、
検出信号と目標値との偏差(アナログ量)が演算され、
この演算偏差は、変換手段でのΔΣ変調によって1ビッ
トデジタル信号に変換され、さらに補償手段により補償
された後、入力信号として制御対象に入力される。以上
の構成により、制御対象の出力と目標値との偏差を打ち
消すように、偏差と逆位相の操作量が生成され、制御対
象に入力される。その結果、制御対象の出力が目標値に
収束するようにフィードバック制御される。
【0004】
【発明が解決しようとする課題】上記従来の制御装置に
よれば、制御対象の動特性が比較的大きな位相遅れやむ
だ時間などを有している場合、それに起因して、制御対
象への入力信号の入力後、入力信号を反映する出力信号
が制御対象から出力されるまでに時間がかかり、制御対
象の入出力間での制御タイミングのずれを生じる。その
結果、制御性の低下を招き、制御系が不安定になってし
まうおそれがある。例えば、内燃機関の燃料噴射量を入
力として、内燃機関の排気ガスの空燃比を制御する場
合、燃料噴射が実際に行われてから、排気ガスの空燃比
の状態が実際に変化するまでに時間を要するため、空燃
比制御の安定性および制御性が低下し、触媒で浄化され
た排気ガスの特性が不安定な状態になってしまうことが
ある。
【0005】以上のような、位相遅れやむだ時間などの
大きな制御対象における制御タイミングのずれを補償す
る制御手法として、例えば特開2000−179385
号公報に記載されたものが知られている。この制御手法
では、同定器および状態予測器により、制御対象の出力
として、制御対象の位相遅れやむだ時間などの応答遅れ
を考慮した予測値が算出される。具体的には、目標空燃
比を入力としかつ空燃比センサの出力を出力とする制御
対象を、その応答遅れを考慮した離散時間系モデルとし
てモデル化し、この離散時間系モデルのモデルパラメー
タを同定器により算出するとともに、状態予測器によ
り、空燃比センサの出力と所定の目標値との偏差の予測
値が算出される。
【0006】しかし、このような制御手法を前記従来の
制御装置に適用した場合、以下のような問題がある。す
なわち、従来の制御装置では、ΔΣ変調アルゴリズムを
用いた制御が実行されるため、制御周期をより短く設定
するほど、制御精度が向上する。言い換えれば、制御周
期をより長く設定するほど、制御精度が低下してしま
う。これに対して、同定器および状態予測器を用いる上
記制御手法の場合、同定器および状態予測器の算出周期
を短く設定することは、モデルパラメータの同定値およ
び予測値の算出精度を向上させることに必ずしもつなが
らず、場合によっては、これらの算出精度を低下させる
おそれがある。例えば、制御対象の入出力、または出力
の目標値の変化が小さい状況では、同定器および状態予
測器の算出周期を短く設定すると、これらの算出値に制
御対象の周波数特性が適切に反映されないことで、算出
精度が逆に低下してしまう。
【0007】本発明は、上記課題を解決するためになさ
れたもので、位相遅れやむだ時間などが比較的大きい制
御対象を制御する場合において、制御対象の入出力間で
の制御タイミングのずれの解消、および制御精度の向上
をいずれも達成することができる制御装置を提供するこ
とを目的とする。
【0008】
【課題を解決するための手段】この目的を達成するため
に、請求項1に係る制御装置1は、制御対象の出力(酸
素濃度センサ15の出力Vout)と所定の目標値Vo
pとの偏差(出力偏差VO2)を所定の偏差算出周期
(算出周期ΔTk)で算出する偏差算出手段(ECU
2、状態予測器22)と、Δ変調アルゴリズム、ΔΣ変
調アルゴリズムおよびΣΔ変調アルゴリズムのうちのい
ずれか1つの変調アルゴリズムに基づき、算出された偏
差に応じて、制御対象の出力を目標値に収束させるため
の制御対象への制御入力(目標空燃比KCMD)を、偏
差算出周期よりも短い所定の制御入力算出周期(算出周
期ΔTm)で算出する制御入力算出手段(ECU2、D
SMコントローラ24)と、を備えることを特徴とす
る。
【0009】この制御装置によれば、偏差算出手段によ
り、制御対象の出力と所定の目標値との偏差が所定の偏
差算出周期で算出され、制御入力算出手段により、制御
対象の出力を目標値に収束させるための制御対象への制
御入力が、Δ変調アルゴリズム、ΔΣ変調アルゴリズム
およびΣΔ変調アルゴリズムのうちのいずれか1つの変
調アルゴリズムに基づき、算出された偏差に応じて、偏
差算出周期よりも短い所定の制御入力算出周期で算出さ
れる。このように、制御対象の出力と所定の目標値との
偏差が、制御入力の算出周期よりも長い偏差算出周期で
算出されるので、制御対象の応答遅れ(位相遅れやむだ
時間など)が大きい場合、または制御対象の出力の目標
値の変動周波数が低い場合でも、制御入力算出周期を可
能な限り短く設定することで、制御入力の算出精度を向
上させることが可能になると同時に、同じ理由により、
その場合、制御入力算出周期の短さに起因する高周波ノ
イズが制御対象の出力に混入したときでも、偏差の算出
値に混入するのを回避でき、偏差を精度よく算出でき
る。以上により、制御性を向上させることができる。
【0010】請求項2に係る発明は、請求項1に記載の
制御装置1において、制御入力算出手段は、1つの変調
アルゴリズムに基づき、偏差に応じて、第1の中間値
(DSM信号u''(m)の値)を算出するとともに、算出
した第1の中間値に所定のゲインFdを乗算した値(増
幅DSM信号u(m)の値)に基づき、制御入力(目標空
燃比KCMD)を算出することを特徴とする。
【0011】一般に、ΔΣ変調アルゴリズム、ΣΔ変調
アルゴリズムおよびΔ変調アルゴリズムの各々は、制御
対象のゲインが値1であるとして制御入力を決定するも
のであるため、制御対象の実際のゲインが値1と異なる
場合には、制御入力が適切な値として算出されなくなる
ことで、制御性が低下することがある。例えば、制御対
象の実際のゲインが値1よりも大きい場合には、制御入
力が必要以上に大きい値として算出されるため、オーバ
ーゲインの状態になってしまうおそれがある。これに対
して、この制御装置によれば、制御入力が、1つの変調
アルゴリズムに基づいて算出された第1の中間値に所定
のゲインを乗算した値に基づいて、算出されるので、こ
の所定のゲインを適切に設定することにより、良好な制
御性を確保することができる。
【0012】請求項3に係る発明は、請求項1に記載の
制御装置1において、制御入力算出手段は、1つの変調
アルゴリズムに基づき、偏差に応じて、第2の中間値
(増幅DSM信号u(m)の値)を算出するとともに、算
出した第2の中間値に所定値(基準値FLAFBAS
E)を加算することにより、制御入力(目標空燃比KC
MD)を算出することを特徴とする。
【0013】一般に、Δ変調アルゴリズム、ΔΣ変調ア
ルゴリズムおよびΣΔ変調アルゴリズムはいずれも、値
0を中心とした正負反転型の制御入力しか算出できない
ものである。これに対して、この制御装置によれば、制
御入力算出手段により、制御入力の算出が、1つの変調
アルゴリズムに基づいて算出された第2の中間値に所定
値を加算することにより行われるので、制御入力を、値
0を中心として正負反転する値だけでなく、所定値を中
心として所定幅の増減を繰り返す値として算出すること
ができ、制御の自由度を高めることができる。
【0014】請求項4に係る制御装置1は、予測アルゴ
リズムに基づき、制御対象の出力を表す値(出力偏差V
O2)の予測値PREVO2を所定の予測値算出周期
(算出周期ΔTk)で算出する予測値算出手段(ECU
2、状態予測器22)と、Δ変調アルゴリズム、ΔΣ変
調アルゴリズムおよびΣΔ変調アルゴリズムのうちのい
ずれか1つの変調アルゴリズムに基づき、算出された予
測値PREVO2に応じて、制御対象の出力を制御する
ための制御対象への制御入力(目標空燃比KCMD)
を、予測値算出周期よりも短い制御入力算出周期(算出
周期ΔTm)で算出する制御入力算出手段(ECU2、
DSMコントローラ24)と、を備えることを特徴とす
る。
【0015】この制御装置によれば、予測値算出手段に
より、制御対象の出力を表す値の予測値が、予測アルゴ
リズムに基づき、所定の予測値算出周期で算出されると
ともに、制御入力算出手段により、制御対象の出力を制
御するための制御対象への制御入力が、Δ変調アルゴリ
ズム、ΔΣ変調アルゴリズムおよびΣΔ変調アルゴリズ
ムのうちのいずれか1つの変調アルゴリズムに基づき、
算出された予測値に応じて、予測値算出周期よりも短い
制御入力算出周期で算出される。前述したように、制御
対象の入出力の変化が小さい状況では、予測値算出手段
による予測値の算出周期を短く設定すると、制御対象の
周波数特性が予測値に適切に反映されないことで、予測
値の算出精度が逆に低下してしまうおそれがある。
【0016】これに対して、この制御装置では、予測値
算出周期が制御入力算出周期よりも長く設定されるの
で、位相遅れやむだ時間の大きい制御対象を制御する場
合、予測値算出周期を、制御入力算出周期よりも長くか
つ制御対象の周波数特性に応じた値に設定することによ
り、予測値を、制御対象の位相遅れやむだ時間などの動
特性が反映された適切な値として、精度よく算出するこ
とが可能になる。その結果、そのような予測値に応じて
制御入力を算出することにより、制御対象の入出力間で
の制御タイミングのずれを適切に解消することが可能に
なる。同じ理由により、予測値算出手段および制御入力
算出手段を、例えばマイクロコンピュータなどの演算処
理装置で構成した場合、予測値算出周期を制御入力算出
周期と同じ値に設定したときと比べて、予測値の不要な
算出処理を省略できる分、演算処理装置の演算負荷を低
減でき、それにより、制御入力の算出周期をより短く設
定できることで、その算出精度を向上させることができ
る。以上により、制御性を向上させることができる。
【0017】以下に述べる請求項5および6に係る発明
はそれぞれ、前述した請求項2および3に係る発明とほ
ぼ同様に構成されている。すなわち、請求項5に係る発
明は、請求項4に記載の制御装置1において、制御入力
算出手段は、1つの変調アルゴリズムに基づき、予測値
PREVO2に応じて第1の中間値(DSM信号u''
(m)の値)を算出するとともに、算出した第1の中間値
に所定のゲインFdを乗算した値(増幅DSM信号u
(m)の値)に基づき、制御入力(目標空燃比KCMD)
を算出することを特徴とする。
【0018】また、請求項6に係る発明は、請求項4に
記載の制御装置1において、制御入力算出手段は、1つ
の変調アルゴリズムに基づき、予測値PREVO2に応
じて第2の中間値(増幅DSM信号u(m)の値)を算出
するとともに、算出した第2の中間値に所定値(基準値
FLAFBASE)を加算することにより、制御入力
(目標空燃比KCMD)を算出することを特徴とする。
【0019】以上の請求項5および6に係る発明によれ
ば、請求項2および3に係る発明と同様の作用効果をそ
れぞれ得ることができる。
【0020】請求項7に係る発明は、請求項4ないし6
のいずれかに記載の制御装置1において、制御対象の出
力を表す値(出力偏差VO2)を、制御入力算出周期
(算出周期ΔTm)よりも長い所定の出力サンプリング
周期(サンプリング周期ΔTsk)で、離散データとし
てサンプリングする出力サンプリング手段(ECU2)
をさらに備え、予測アルゴリズムは、サンプリングされ
た制御対象の出力を表す値の離散データ(VO2
(k),VO2(k−1))に応じて、予測値PREV
O2を算出するアルゴリズムであることを特徴とする。
【0021】この制御装置によれば、出力サンプリング
手段により、制御対象の出力を表す値の離散データが、
制御入力算出周期よりも長い所定の出力サンプリング周
期でサンプリングされるとともに、予測値算出手段によ
り、予測アルゴリズムに基づき、サンプリングされた制
御対象の出力を表す値の離散データに応じて、予測値が
算出される。したがって、出力サンプリング周期を、制
御対象の出力を表す値のパワースペクトルが存在する周
波数域に応じた適切な値に設定することにより、制御対
象の入出力の変化が小さい状況でも、制御対象の出力を
表す値の予測値の周波数特性を、制御対象の周波数特性
に一致させることができる。その結果、制御性をさらに
向上させることができる。
【0022】請求項8に係る発明は、請求項7に記載の
制御装置1において、制御入力を表す値(空燃比偏差の
デシメーション値DKCMD’)、および制御対象に入
力された制御入力を反映する値(LAF出力偏差のデシ
メーション値DKACT’)の少なくとも一方を、制御
入力算出周期よりも長い所定の入力サンプリング周期
(サンプリング周期ΔTsk)で、離散データとしてサ
ンプリングする入力サンプリング手段(ECU2、デシ
メーションフィルタ40)をさらに備え、予測アルゴリ
ズムは、サンプリングされた、制御入力を表す値の離散
データ(DKCMD’(k−i))、および制御対象に
入力された制御入力を反映する値の離散データ(DKA
CT’(k−j))の少なくとも一方にさらに応じて、
予測値PREVO2を算出するアルゴリズムであること
を特徴とする。
【0023】この制御装置によれば、入力サンプリング
手段により、制御入力を表す値の離散データおよび/ま
たは制御対象に入力された制御入力を反映する値の離散
データが、制御入力算出周期よりも長い所定の入力サン
プリング周期でサンプリングされるとともに、予測値算
出手段により、予測アルゴリズムに基づき、サンプリン
グされた離散データに応じて、予測値が算出される。こ
のように、制御対象の出力を表す値の予測値が、制御入
力を表す値の離散データおよび/または制御対象に入力
された制御入力を反映する値の離散データにさらに応じ
て算出されるので、位相遅れやむだ時間の大きい制御対
象を制御する場合でも、予測値を、制御入力算出周期の
短さに起因する高周波ノイズが予測値の算出値に混入す
るのを回避しながら、制御対象の位相遅れやむだ時間な
どの動特性が反映された適切な値として、精度よく算出
することができ、制御対象の入出力間の制御タイミング
のずれを精度よく解消することができる。その結果、制
御性をより一層、向上させることができる。
【0024】請求項9に係る発明は、請求項8に記載の
制御装置1において、入力サンプリング手段は、制御入
力を表す値の離散データ(DKCMD’(k−i))、
および制御対象に入力された制御入力を反映する値の離
散データ(DKACT’(k−j))の少なくとも一方
を、入力サンプリング周期(サンプリング周期ΔTs
k)よりも短い所定の第2入力サンプリング周期(サン
プリング周期ΔTsm)でサンプリングした値にデシメ
ーション処理を施すことにより、入力サンプリング周期
でサンプリングすることを特徴とする。
【0025】この制御装置によれば、制御入力を表す値
の離散データおよび/または制御対象に入力された制御
入力を反映する値の離散データが、入力サンプリング周
期よりも短い所定の第2入力サンプリング周期でサンプ
リングされた値に、デシメーション処理を施すことによ
り、入力サンプリング周期でサンプリングされる。これ
により、制御対象の出力を表す値の予測値の算出におい
て、制御入力算出周期の短さに起因する高周波ノイズが
予測値の算出値に混入するのを、確実に回避することが
できる。また、第2入力サンプリング周期を、制御入力
算出周期と同じ値に設定した場合には、演算負荷を低減
することができる。
【0026】請求項10に係る発明は、請求項4ないし
6のいずれかに記載の制御装置1において、予測アルゴ
リズムは、制御入力を表す値(空燃比偏差のデシメーシ
ョン値DKCMD’)および制御対象に入力された制御
入力を反映する値(LAF出力偏差のデシメーション値
DKACT’)の一方と、制御対象の出力を表す値(出
力偏差VO2)とを変数とする制御対象モデル(式(6
0))に基づき、予測値PREVO2を算出するアルゴ
リズム(式(61),(76))であることを特徴とす
る。
【0027】この制御装置によれば、制御入力を表す値
または制御対象に入力された制御入力を反映する値と、
制御対象の出力を表す値とを変数とする制御対象モデル
に基づいて、予測値が算出されるので、この制御対象モ
デルを、制御対象の位相遅れやむだ時間などの動特性が
反映されたものとして定義することにより、予測値を、
制御対象の位相遅れやむだ時間などの動特性が適切に反
映された値として算出することができる。その結果、制
御の安定性を確保することができ、制御性を向上させる
ことができる。
【0028】請求項11に係る発明は、請求項10に記
載の制御装置1において、制御入力を表す値(空燃比偏
差のデシメーション値DKCMD’)と制御対象に入力
された制御入力を反映する値(LAF出力偏差のデシメ
ーション値DKACT’)との少なくとも一方、および
制御対象の出力を表す値(出力偏差VO2)を、制御入
力算出周期(算出周期ΔTm)よりも長い所定の第1サ
ンプリング周期(サンプリング周期ΔTsk)で、離散
データとしてサンプリングするサンプリング手段(EC
U2、デシメーションフィルタ40)をさらに備え、制
御対象モデルは、サンプリングされた、制御入力を表す
値の離散データ(DKCMD’(k−dt))および制
御対象に入力された制御入力を反映する値の離散データ
の一方と、制御対象の出力を表す値の離散データ(VO
2(k−1),VO2(k−2))とを変数とする離散
時間系モデル(式(60))であることを特徴とする。
【0029】この制御装置によれば、制御対象モデル
が、制御入力算出周期よりも長い所定の第1サンプリン
グ周期でサンプリングされた、制御入力を表す値の離散
データまたは制御対象に入力された制御入力を反映する
値の離散データと、制御対象の出力を表す値の離散デー
タとを変数とする離散時間系モデルであるので、第1サ
ンプリング周期を適切に設定することより、この離散時
間系モデルに、制御対象の周波数特性を適切に反映させ
ることができ、その結果、予測値の算出精度を向上させ
ることができる。また、制御対象モデルとして離散時間
系モデルを用いるので、連続時間系モデルを用いる場合
と比べて、例えば最小2乗法などの一般的な同定アルゴ
リズムにより、モデルパラメータの同定を容易に行うこ
とが可能になる。
【0030】請求項12に係る発明は、請求項11に記
載の制御装置1において、サンプリング手段は、制御入
力を表す値の離散データ(DKCMD’(k−d
t))、および制御対象に入力された制御入力を反映す
る値の離散データの少なくとも一方を、第1サンプリン
グ周期(サンプリング周期ΔTsk)よりも短い所定の
第2サンプリング周期(サンプリング周期ΔTsm)で
サンプリングした値にデシメーション処理を施すことに
より、第1サンプリング周期でサンプリングすることを
特徴とする。
【0031】この制御装置によれば、請求項9に係る発
明と同様の作用効果を得ることができる。
【0032】請求項13に係る発明は、請求項11また
は12に記載の制御装置1において、制御対象の出力を
表す値の離散データ(VO2(k−1),VO2(k−
2))に応じて、離散時間系モデルのモデルパラメータ
a1,a2,b1を同定する同定手段(ECU2、オン
ボード同定器23)をさらに備えることを特徴とする。
【0033】この制御装置によれば、同定手段により、
制御対象の出力を表す値の離散データに応じて、離散時
間系モデルのモデルパラメータが同定されるので、制御
対象の動特性が経時変化したり、ばらついたりする場合
でも、それに応じてモデルパラメータを適切に同定する
ことができ、制御対象モデルの動特性を制御対象の実際
の動特性に適合させることができる。その結果、制御性
および制御の安定性を向上させることができる。これに
加えて、離散時間系モデルを用いることにより、同定手
段によるモデルパラメータの同定を容易に行うことがで
きる。
【0034】請求項14に係る発明は、請求項13に記
載の制御装置1において、同定手段は、制御入力を表す
値の離散データ(DKCMD’(k−dt))、および
制御対象に入力された制御入力を反映する値の離散デー
タ(DKACT’(k−d−dd))の一方にさらに応
じて、モデルパラメータを同定することを特徴とする。
【0035】この制御装置によれば、同定手段により、
制御対象の出力を表す値の離散データに加えて、制御入
力を表す値の離散データ、または制御対象に入力された
制御入力を反映する値の離散データにさらに応じて、モ
デルパラメータが同定されるので、制御対象の入出力間
の相関関係を観測しながら、モデルパラメータを同定す
ることができ、その同定精度を高めることができる。そ
の結果、制御性および制御の安定性をさらに向上させる
ことができる。
【0036】請求項15に係る発明は、請求項13また
は14に記載の制御装置1において、同定手段によるモ
デルパラメータの同定周期(算出周期ΔTk)は、制御
入力算出周期(算出周期ΔTm)よりも長い値に設定さ
れていることを特徴とする。
【0037】一般に、演算処理装置においては、同定器
などの同定手段によるモデルパラメータの同定演算が過
大な演算負荷となることで、制御入力の算出周期をより
長く設定しなければならないことがあり、その場合に
は、制御性の低下を招くおそれがある。これに対して、
この制御装置によれば、同定周期が制御入力算出周期よ
りも長い値に設定されているので、そのような演算処理
装置における演算負荷を低減することができ、その分、
制御入力の算出周期をより短く設定できることで、制御
性をより一層、向上させることができる。
【0038】請求項16に係る制御装置1は、内燃機関
3の排気通路(排気管7)の触媒(第1触媒装置8a、
第2触媒装置8b)よりも下流側における排気ガスの空
燃比を表す検出信号を出力する下流側空燃比センサ(酸
素濃度センサ15)と、下流側空燃比センサの出力Vo
utと所定の目標値Vopとの偏差(出力偏差VO2)
を所定の偏差算出周期(算出周期ΔTk)で算出する偏
差算出手段(ECU2、状態予測器22)と、Δ変調ア
ルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アル
ゴリズムのうちのいずれか1つの変調アルゴリズムに基
づき、算出された偏差に応じて、下流側空燃比センサの
出力を目標値に収束させるための、内燃機関に供給すべ
き混合気の目標空燃比KCMDを、偏差算出周期よりも
短い所定の空燃比算出周期(算出周期ΔTm)で算出す
る空燃比算出手段(ECU2、DSMコントローラ2
4)と、算出された目標空燃比KCMDに応じて、内燃
機関に供給される混合気の空燃比を制御する空燃比制御
手段(ECU2)と、を備えることを特徴とする。
【0039】この制御装置によれば、偏差算出手段によ
り、触媒よりも下流側における排気ガスの空燃比を表す
下流側空燃比センサの出力と、所定の目標値との偏差が
所定の偏差算出周期で算出され、制御入力算出手段によ
り、下流側空燃比センサの出力を目標値に収束させるた
めの、内燃機関に供給すべき混合気の目標空燃比が、Δ
変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変
調アルゴリズムのうちのいずれか1つの変調アルゴリズ
ムに基づき、算出された偏差に応じて、偏差算出周期よ
りも短い所定の空燃比算出周期で算出される。このよう
に、下流側空燃比センサの出力と目標値との偏差が、目
標空燃比の算出周期よりも長い偏差算出周期で算出され
るので、内燃機関への目標空燃比の混合気の供給に対し
て、下流側空燃比センサの出力の応答遅れ(位相遅れや
むだ時間など)が大きい空燃比制御の場合や、目標値の
変動周波数が低い場合においても、空燃比算出周期を可
能な限り短く設定することで、目標空燃比の算出精度を
向上させることができる。同じ理由により、空燃比算出
周期を可能な限り短く設定した場合、空燃比算出周期の
短さに起因する高周波ノイズが、下流側空燃比センサの
出力に混入したときでも、偏差の算出値に混入するのを
回避でき、偏差を精度よく算出できる。以上により、目
標空燃比を、下流側空燃比センサの出力を目標値に精度
よく収束させることができる値として算出でき、空燃比
の制御性を向上させることができる。したがって、この
目標値を適切な値に設定することにより、触媒による排
気ガス浄化率を向上させることができる。
【0040】請求項17に係る発明は、請求項16に記
載の制御装置1において、空燃比算出手段は、1つの変
調アルゴリズムに基づき、偏差に応じて、第1の中間値
(DSM信号u''(m)の値)を算出するとともに、算出
した第1の中間値に所定のゲインFdを乗算した値(増
幅DSM信号u(m)の値)に基づき、目標空燃比KCM
Dを算出することを特徴とする。
【0041】一般に、ΔΣ変調アルゴリズム、ΣΔ変調
アルゴリズムおよびΔ変調アルゴリズムの各々は、制御
対象のゲインが値1であるとして制御入力すなわち目標
空燃比を決定するものであるので、制御対象の実際のゲ
インが値1と異なる場合には、目標空燃比が適切な値と
して算出されなくなることで、空燃比の制御性が低下す
ることがある。例えば、制御対象の実際のゲインが値1
よりも大きい場合には、目標空燃比が必要以上に大きい
値として算出されるため、オーバーゲインの状態になっ
てしまうおそれがある。これに対して、この制御装置に
よれば、目標空燃比が、1つの変調アルゴリズムに基づ
いて算出された第1の中間値に所定のゲインを乗算した
値に基づいて、算出されるので、この所定のゲインを適
切に設定することにより、良好な空燃比の制御性を確保
することができる。
【0042】請求項18に係る発明は、請求項16に記
載の制御装置1において、空燃比算出手段は、1つの変
調アルゴリズムに基づき、偏差に応じて、第2の中間値
(増幅DSM信号u(m)の値)を算出するとともに、算
出した第2の中間値に所定値(基準値FLAFBAS
E)を加算することにより、目標空燃比KCMDを算出
することを特徴とする。
【0043】一般に、Δ変調アルゴリズム、ΔΣ変調ア
ルゴリズムおよびΣΔ変調アルゴリズムはいずれも、値
0を中心とした正負反転型の制御入力しか算出できない
ものである。これに対して、この制御装置によれば、空
燃比算出手段により、制御入力としての目標空燃比の算
出が、1つの変調アルゴリズムに基づいて算出された第
2の中間値に所定値を加算することにより行われるの
で、目標空燃比を、値0を中心として正負反転する値だ
けでなく、所定値を中心として所定幅の増減を繰り返す
値として算出することができ、空燃比制御の自由度を高
めることができる。
【0044】請求項19に係る制御装置1は、内燃機関
3の排気通路(排気管7)の触媒(第1触媒装置8a、
第2触媒装置8b)よりも下流側における排気ガスの空
燃比を表す検出信号を出力する下流側空燃比センサ(酸
素濃度センサ15)と、予測アルゴリズム(式(6
1),(76))に基づき、下流側空燃比センサの出力
を表す値(出力偏差VO2)の予測値PREVO2を所
定の予測値算出周期(算出周期ΔTk)で算出する予測
値算出手段(ECU2、状態予測器22)と、Δ変調ア
ルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ変調アル
ゴリズムのうちのいずれか1つの変調アルゴリズムに基
づき、算出された予測値に応じて、下流側空燃比センサ
の出力を制御するための、内燃機関に供給すべき混合気
の目標空燃比を、予測値算出周期よりも短い空燃比算出
周期(算出周期ΔTm)で算出する空燃比算出手段(E
CU2、DSMコントローラ24)と、算出された目標
空燃比KCMDに応じて、内燃機関に供給される混合気
の空燃比を制御する空燃比制御手段(ECU2)と、を
備えることを特徴とする。
【0045】この制御装置によれば、予測値算出手段に
より、下流側空燃比センサの出力を表す値の予測値が、
予測アルゴリズムに基づき、所定の予測値算出周期で算
出されるとともに、空燃比算出手段により、下流側空燃
比センサの出力を制御するための、内燃機関に供給すべ
き混合気の目標空燃比が、Δ変調アルゴリズム、ΔΣ変
調アルゴリズムおよびΣΔ変調アルゴリズムのうちのい
ずれか1つの変調アルゴリズムに基づき、算出された予
測値に応じて、予測値算出周期よりも短い空燃比算出周
期で算出される。前述したように、制御対象の入出力、
すなわち目標空燃比および下流側空燃比センサの出力の
変化が小さい状況では、予測値算出手段による予測値の
算出周期を短く設定すると、制御対象の周波数特性が予
測値に適切に反映されないことで、予測値の算出精度が
逆に低下してしまうおそれがある。
【0046】これに対して、この制御装置では、予測値
算出周期が空燃比算出周期よりも長く設定されるので、
内燃機関に供給する混合気の空燃比によって下流側空燃
比センサの出力を制御する場合、すなわち位相遅れやむ
だ時間の大きい制御対象を制御する場合、予測値算出周
期を、空燃比算出周期よりも長くかつ制御対象の周波数
特性に応じた値に設定することにより、予測値を、空燃
比算出周期が短いことに起因して目標空燃比の算出値に
混入する高周波ノイズの影響を回避しながら、制御対象
の位相遅れやむだ時間などの動特性が反映された適切な
値として、精度よく算出することが可能になる。その結
果、そのような予測値に応じて目標空燃比を算出するこ
とにより、目標空燃比に応じて内燃機関に供給される混
合気の空燃比と、下流側空燃比センサの出力との間での
制御タイミングのずれを適切に解消することが可能にな
る。同じ理由により、予測値算出手段および空燃比算出
手段を、例えばマイクロコンピュータなどの演算処理装
置で構成した場合、予測値算出周期を空燃比算出周期と
同じ値に設定したときと比べて、予測値の不要な算出処
理を省略できる分、演算処理装置の演算負荷を低減で
き、それにより、目標空燃比の算出周期をより短く設定
できることで、その算出精度を向上させることができ
る。以上により、空燃比制御の制御性を向上させること
ができ、排気ガス浄化率を向上させることができる。
【0047】以下に述べる請求項20および21に係る
発明はそれぞれ、前述した請求項17および18に係る
発明とほぼ同様に構成されている。すなわち、請求項2
0に係る発明は、請求項19に記載の制御装置1におい
て、空燃比算出手段は、1つの変調アルゴリズムに基づ
き、予測値に応じて、第1の中間値(DSM信号u''
(m)の値)を算出するとともに、算出した第1の中間値
に所定のゲインFdを乗算した値(増幅DSM信号u
(m)の値)に基づき、目標空燃比KCMDを算出するこ
とを特徴とする。
【0048】また、請求項21に係る発明は、請求項1
9に記載の制御装置1において、空燃比算出手段は、1
つの変調アルゴリズムに基づき、予測値に応じて、第2
の中間値(増幅DSM信号u(m)の値)を算出するとと
もに、算出した第2の中間値に所定値(基準値FLAF
BASE)を加算することにより、目標空燃比KCMD
を算出することを特徴とする。
【0049】以上の請求項20および21に係る発明に
よれば、請求項17および18に係る発明と同様の作用
効果をそれぞれ得ることができる。
【0050】請求項22に係る発明は、請求項19ない
し21のいずれかに記載の制御装置1において、下流側
空燃比センサの出力を表す値(出力偏差VO2)を、空
燃比算出周期(算出周期ΔTm)よりも長い所定の出力
サンプリング周期(サンプリング周期ΔTsk)で、離
散データとしてサンプリングする出力サンプリング手段
(ECU2)をさらに備え、予測アルゴリズムは、サン
プリングされた下流側空燃比センサの出力を表す値の離
散データ(VO2(k),VO2(k−1))に応じ
て、予測値PREVO2を算出するアルゴリズムである
ことを特徴とする。
【0051】この制御装置によれば、出力サンプリング
手段により、下流側空燃比センサの出力を表す値の離散
データが、空燃比算出周期よりも長い所定の出力サンプ
リング周期でサンプリングされるとともに、予測値算出
手段により、予測アルゴリズムに基づき、サンプリング
された下流側空燃比センサの出力を表す値の離散データ
に応じて、予測値が算出される。したがって、出力サン
プリング周期を、下流側空燃比センサの出力を表す値の
パワースペクトルが存在する周波数域に応じて適切に設
定することにより、目標空燃比および下流側空燃比セン
サの出力の変化が小さい状況でも、下流側空燃比センサ
の出力を表す値の予測値の周波数特性を、制御対象の周
波数特性に一致させることができる。その結果、制御性
をさらに向上させることができ、排気ガス浄化率をさら
に向上させることができる。
【0052】請求項23に係る発明は、請求項22に記
載の制御装置1において、目標空燃比を表す値(空燃比
偏差のデシメーション値DKCMD’)を、空燃比算出
周期(算出周期ΔTm)よりも長い所定の空燃比サンプ
リング周期(サンプリング周期ΔTsk)で、離散デー
タとしてサンプリングする空燃比サンプリング手段(E
CU2)をさらに備え、予測アルゴリズムは、サンプリ
ングされた目標空燃比を表す値の離散データ(DKCM
D’(k−i))にさらに応じて、予測値PREVO2
を算出するアルゴリズム(式(76))であることを特
徴とする。
【0053】この制御装置によれば、空燃比サンプリン
グ手段により、目標空燃比を表す値の離散データが、空
燃比算出周期よりも長い所定の空燃比サンプリング周期
でサンプリングされるとともに、予測値算出手段によ
り、予測アルゴリズムに基づき、サンプリングされた目
標空燃比を表す値の離散データにさらに応じて、予測値
が算出される。このように、下流側空燃比センサの出力
を表す値の予測値が、目標空燃比を表す値の離散データ
に応じて算出されるので、内燃機関の空燃比を制御する
場合、すなわち位相遅れやむだ時間の大きい制御対象を
制御する場合でも、空燃比算出周期の短さに起因する高
周波ノイズが、目標空燃比の算出値に混入するのを回避
しながら、制御対象の位相遅れやむだ時間などの動特性
が反映された適切な値として、予測値を精度よく算出す
ることができる、これにより、目標空燃比と下流側空燃
比センサの出力を表す値との間、すなわち制御対象の入
出力間の制御タイミングのずれを精度よく解消すること
ができ、空燃比制御の制御性をより一層、向上させるこ
とができる。
【0054】請求項24に係る発明は、請求項23に記
載の制御装置1において、空燃比サンプリング手段は、
目標空燃比を表す値の離散データ(DKCMD’(k−
i))を、空燃比サンプリング周期(サンプリング周期
ΔTsk)よりも短い所定の第2空燃比サンプリング周
期(サンプリング周期ΔTsm)でサンプリングした値
にデシメーション処理を施すことにより、空燃比サンプ
リング周期でサンプリングすることを特徴とする。
【0055】この制御装置によれば、目標空燃比を表す
値の離散データが、空燃比サンプリング周期よりも短い
所定の第2空燃比サンプリング周期でサンプリングされ
た値にデシメーション処理を施すことにより、空燃比サ
ンプリング周期でサンプリングされる。これにより、下
流側空燃比センサの出力を表す値の予測値の算出におい
て、空燃比算出周期の短さに起因する高周波ノイズが目
標空燃比の算出値に混入するのを、確実に回避すること
ができる。また、第2空燃比サンプリング周期を、空燃
比算出周期と同じ値に設定した場合には、演算負荷を低
減することができる。
【0056】請求項25に係る発明は、請求項19ない
し21のいずれかに記載の制御装置1において、予測ア
ルゴリズムは、目標空燃比を表す値(空燃比偏差のデシ
メーション値DKCMD’)および下流側空燃比センサ
の出力を表す値(出力偏差VO2)を変数とする制御対
象モデル(式(60))に基づき、予測値PREVO2
を算出するアルゴリズムであることを特徴とする。
【0057】この制御装置によれば、目標空燃比を表す
値および下流側空燃比センサの出力を表す値を変数とす
る制御対象モデルに基づき、下流側空燃比センサの出力
を表す値の予測値が算出されるので、この制御対象モデ
ルを、制御対象の位相遅れやむだ時間などの動特性が反
映されたものとして定義することにより、予測値を、制
御対象の位相遅れやむだ時間などの動特性が適切に反映
された値として算出することができる。その結果、空燃
比制御の安定性を確保することができ、その制御性を向
上させることができる。
【0058】請求項26に係る発明は、請求項25に記
載の制御装置1において、目標空燃比を表す値(空燃比
偏差のデシメーション値DKCMD’)および下流側空
燃比センサの出力を表す値(出力偏差VO2)を、空燃
比算出周期(算出周期ΔTm)よりも長い所定の第1サ
ンプリング周期(サンプリング周期ΔTsk)で、離散
データとしてサンプリングするサンプリング手段(EC
U2)をさらに備え、制御対象モデルは、サンプリング
された、目標空燃比を表す値の離散データ(DKCM
D’(k−dt))と、下流側空燃比センサの出力を表
す値の離散データ(VO2(k−1),VO2(k−
2))とを変数とする離散時間系モデル(式(60))
であることを特徴とする。
【0059】この制御装置によれば、制御対象モデル
が、空燃比算出周期よりも長い所定の第1サンプリング
周期でサンプリングされた、目標空燃比を表す値の離散
データと、下流側空燃比センサの出力を表す値の離散デ
ータとを変数とする離散時間系モデルであるので、第1
サンプリング周期を適切に設定することより、この離散
時間系モデルに、目標空燃比を表す値および下流側空燃
比センサの出力を表す値を入出力とする制御対象の周波
数特性を適切に反映させることができ、その結果、予測
値の算出精度を向上させることができる。また、制御対
象モデルとして離散時間系モデルを用いるので、連続時
間系モデルを用いる場合と比べて、例えば最小2乗法な
どの一般的な同定アルゴリズムにより、モデルパラメー
タの同定を容易に行うことが可能になる。
【0060】請求項27に係る発明は、請求項26に記
載の制御装置1において、サンプリング手段は、目標空
燃比を表す値の離散データ(DKCMD’(k−d
t))を、第1サンプリング周期(サンプリング周期Δ
Tsk)よりも短い所定の第2サンプリング周期(サン
プリング周期ΔTsm)でサンプリングした値にデシメ
ーション処理を施すことにより、第1サンプリング周期
でサンプリングすることを特徴とする。
【0061】この制御装置によれば、請求項24に係る
発明と同様の作用効果を得ることができる。
【0062】請求項28に係る発明は、請求項請求項2
6または27に記載の制御装置1において、下流側空燃
比センサの出力を表す値の離散データ(VO2(k−
1),VO2(k−2))に応じて、離散時間系モデル
のモデルパラメータa1,a2,b1を同定する同定手
段(オンボード同定器23)をさらに備えることを特徴
とする。
【0063】この制御装置によれば、同定手段により、
下流側空燃比センサの出力を表す値の離散データに応じ
て、離散時間系モデルのモデルパラメータが同定される
ので、目標空燃比を表す値および下流側空燃比センサの
出力を表す値を入出力とする制御対象において、その動
特性が経時変化したり、ばらついたりする場合でも、そ
れに応じてモデルパラメータを適切に同定することがで
き、制御対象モデルの動特性を制御対象の実際の動特性
に適合させることができる。その結果、空燃比制御の精
度および安定性を向上させることができ、排気ガス浄化
率を向上させることができる。これに加えて、離散時間
系モデルを用いることにより、同定手段によるモデルパ
ラメータの同定を容易に行うことができる。
【0064】請求項29に係る発明は、請求項28に記
載の制御装置1において、同定手段は、目標空燃比を表
す値の離散データ(DKCMD’(k−dt))にさら
に応じて、モデルパラメータを同定することを特徴とす
る。
【0065】この制御装置によれば、同定手段により、
下流側空燃比センサの出力を表す値の離散データに加え
て、目標空燃比を表す値の離散データにさらに応じて、
モデルパラメータが同定されるので、制御対象の入出力
間の相関関係を観測しながら、モデルパラメータを同定
することができ、その同定精度を高めることができる。
その結果、空燃比制御の精度および安定性をさらに向上
させることができ、排気ガス浄化率をさらに向上させる
ことができる。
【0066】請求項30に係る発明は、請求項28また
は29に記載の制御装置1において、同定手段によるモ
デルパラメータの同定周期(算出周期ΔTk)は、空燃
比算出周期(算出周期ΔTm)よりも長い値に設定され
ていることを特徴とする。
【0067】一般に、演算処理装置においては、同定器
などの同定手段によるモデルパラメータの同定演算が過
大な演算負荷となることで、制御入力としての目標空燃
比を表す値の算出周期をより長く設定しなければならな
いことがあり、その場合には、制御性の低下を招くおそ
れがある。これに対して、この制御装置によれば、同定
周期が空燃比算出周期よりも長い値に設定されているの
で、そのような演算処理装置における演算負荷を低減す
ることができ、その分、目標空燃比を表す値の算出周期
をより短く設定できることで、空燃比制御の精度をより
一層、向上させることができ、排気ガス浄化率をより一
層、向上させることができる。
【0068】請求項31に係る発明は、請求項22に記
載の制御装置1において、内燃機関3の排気通路(排気
管7)の触媒(第1触媒装置8a)よりも上流側におけ
る排気ガスの空燃比を表す検出信号を出力する上流側空
燃比センサ(LAFセンサ14)と、目標空燃比を表す
値(空燃比偏差のデシメーション値DKCMD’)と上
流側空燃比センサの出力を表す値(LAF出力偏差のデ
シメーション値DKACT’)との少なくとも一方を、
空燃比算出周期(算出周期ΔTm)よりも長い所定の空
燃比サンプリング周期(サンプリング周期ΔTsk)
で、離散データとしてサンプリングする空燃比サンプリ
ング手段(ECU2)と、をさらに備え、予測アルゴリ
ズムは、サンプリングされた、目標空燃比を表す値の離
散データ(DKCMD’(k−i))、および上流側空
燃比センサの出力を表す値の離散データ(DKACT’
(k−j))の少なくとも一方にさらに応じて、予測値
PREVO2を算出するアルゴリズム(式(61))で
あることを特徴とする。
【0069】この制御装置によれば、下流側空燃比セン
サの出力を表す値の離散データに加えて、目標空燃比を
表す値の離散データにさらに応じて、予測値を算出する
場合には、請求項23に係る発明と同様の作用効果を得
ることができる。また、下流側空燃比センサの出力を表
す値の離散データに加えて、またはこれと目標空燃比を
表す値の離散データとに加えて、上流側空燃比センサの
出力を表す値の離散データにさらに応じて、予測値を算
出する場合には、触媒よりも上流側の空燃比の状態を予
測値に反映させることができ、これにより、予測値の算
出精度をさらに向上させることができる。
【0070】請求項32に係る発明は、請求項31に記
載の制御装置1において、サンプリング手段は、目標空
燃比を表す値の離散データ(DKCMD’(k−
i))、および上流側空燃比センサの出力を表す値の離
散データ(DKACT’(k−j))の少なくとも一方
を、空燃比サンプリング周期(サンプリング周期ΔTs
k)よりも短い所定の第2空燃比サンプリング周期(サ
ンプリング周期ΔTsm)でサンプリングした値にデシ
メーション処理を施すことにより、空燃比サンプリング
周期でサンプリングすることを特徴とする。
【0071】この制御装置によれば、第2空燃比サンプ
リング手段およびデシメーション処理手段により、目標
空燃比を表す値の離散データをサンプリングする場合に
は、請求項24に係る発明と同様の作用効果を得ること
ができる。また、上流側空燃比センサの出力を表す値の
離散データ、またはこれに加えて目標空燃比を表す値の
離散データをサンプリングする場合には、空燃比算出周
期の短さに起因する高周波ノイズが、上流側空燃比セン
サの出力を表す値の離散データに混入するのを、確実に
回避することができ、下流側空燃比センサの出力を表す
値の予測値の算出精度をより一層、向上させることがで
きる。
【0072】請求項33に係る発明は、請求項19ない
し21のいずれかに記載の制御装置1において、内燃機
関3の排気通路(排気管7)の触媒(第1触媒装置8
a)よりも上流側における排気ガスの空燃比を表す検出
信号を出力する上流側空燃比センサ(LAFセンサ1
4)をさらに備え、予測アルゴリズムは、目標空燃比を
表す値(空燃比偏差のデシメーション値DKCMD’)
および上流側空燃比センサの出力を表す値(LAF出力
偏差のデシメーション値DKACT’)の一方と、下流
側空燃比センサの出力を表す値(出力偏差VO2)とを
変数とする制御対象モデルに基づき、予測値PREVO
2を算出するアルゴリズムであることを特徴とする。
【0073】この制御装置によれば、制御対象モデル
を、目標空燃比を表す値と下流側空燃比センサの出力を
表す値とを変数とするモデルとした場合には、請求項2
5に係る発明と同様の作用効果を得ることができる。ま
た、制御対象モデルを、目標空燃比を表す値と上流側空
燃比センサの出力を表す値とを変数とするモデルとした
場合には、触媒よりも上流側の空燃比の状態を制御対象
モデルに反映させることができ、これにより、予測値の
算出精度をさらに向上させることができる。
【0074】請求項34に係る発明は、請求項33に記
載の制御装置1において、目標空燃比を表す値(空燃比
偏差のデシメーション値DKCMD’)と上流側空燃比
センサの出力を表す値(LAF出力偏差のデシメーショ
ン値DKACT’)との少なくとも一方、および下流側
空燃比センサの出力を表す値(出力偏差VO2)を、空
燃比算出周期(算出周期ΔTm)よりも長い所定の第1
サンプリング周期(サンプリング周期ΔTsk)で、離
散データとしてサンプリングするサンプリング手段(E
CU2)をさらに備え、制御対象モデルは、サンプリン
グされた、目標空燃比を表す値の離散データ(DKCM
D’(k−dt))および上流側空燃比センサの出力を
表す値の離散データの一方と、下流側空燃比センサの出
力を表す値の離散データ(VO2(k−1),VO2
(k−2))とを変数とする離散時間系モデル(式(6
0))であることを特徴とする。
【0075】この制御装置によれば、制御対象モデル
を、目標空燃比を表す値の離散データと、下流側空燃比
センサの出力を表す値の離散データとを変数とする離散
時間系モデルとした場合には、請求項26に係る発明と
同様の作用効果を得ることができる。また、制御対象モ
デルを、上流側空燃比センサの出力を表す値の離散デー
タと、下流側空燃比センサの出力を表す値の離散データ
とを変数とする離散時間系モデルとした場合には、触媒
よりも上流側の空燃比の状態を制御対象モデルすなわち
予測値に反映させることができ、これにより、予測値の
算出精度をより一層、向上させることができる。
【0076】請求項35に係る発明は、請求項34に記
載の制御装置1において、サンプリング手段は、目標空
燃比を表す値の離散データ(DKCMD’(k−
i))、および上流側空燃比センサの出力を表す値の離
散データ(DKACT’(k−j))の少なくとも一方
を、第1サンプリング周期(サンプリング周期ΔTs
k)よりも短い所定の第2サンプリング周期(サンプリ
ング周期ΔTsm)でサンプリングした値にデシメーシ
ョン処理を施すことにより、第1サンプリング周期でサ
ンプリングすることを特徴とする。
【0077】この制御装置によれば、請求項32に係る
発明と同様の作用効果を得ることができる。
【0078】請求項36に係る発明は、請求項34また
は35に記載の制御装置1において、下流側空燃比セン
サの出力を表す値の離散データ(VO2(k−1),V
O2(k−2))に応じて、離散時間系モデルのモデル
パラメータa1,a2,b1を同定する同定手段(オン
ボード同定器23)をさらに備えることを特徴とする。
【0079】この制御装置によれば、請求項28に係る
発明と同様の作用効果を得ることができる。
【0080】請求項37に係る発明は、請求項36に記
載の制御装置1において、同定手段は、目標空燃比を表
す値の離散データ、および上流側空燃比センサの出力を
表す値の離散データ(DKACT’(k−d−dd))
の一方にさらに応じて、モデルパラメータを同定するこ
とを特徴とする。
【0081】この制御装置によれば、同定手段により、
モデルパラメータが目標空燃比を表す値の離散データに
さらに応じて同定される場合には、請求項29に係る発
明と同様の作用効果を得ることができる。また、モデル
パラメータが上流側空燃比センサの出力を表す値の離散
データにさらに応じて同定される場合には、目標空燃比
を表す値と比べて、下流側空燃比センサの出力を表す値
との間のむだ時間がより小さい上流側空燃比センサの出
力を表す値の離散データを用い、両者の間の相関関係を
観測しながら、モデルパラメータを同定することがで
き、これにより、モデルパラメータの同定精度をさらに
高めることができる。
【0082】請求項38に係る発明は、請求項36また
は37に記載の制御装置1において、同定手段によるモ
デルパラメータの同定周期(算出周期ΔTk)は、空燃
比算出周期(算出周期ΔTm)よりも長い値に設定され
ていることを特徴とする。
【0083】この制御装置によれば、請求項30に係る
発明と同様の作用効果を得ることができる。
【0084】
【発明の実施の形態】以下、図面を参照しながら、本発
明の第1実施形態に係る制御装置について説明する。こ
の第1実施形態は、制御装置を内燃機関の空燃比を制御
するものとして構成した例であり、図1は、この制御装
置1およびこれを適用した内燃機関3の概略構成を示し
ている。同図に示すように、この制御装置1は、ECU
2を備えており、このECU2は、後述するように、内
燃機関(以下「エンジン」という)3の運転状態に応じ
て、これに供給する混合気の空燃比を制御する。
【0085】このエンジン3は、図示しない車両に搭載
された直列4気筒型ガソリンエンジンであり、第1〜第
4の4つの気筒#1〜#4を備えている。このエンジン
3の吸気管4のスロットル弁5の近傍には、例えばポテ
ンショメータなどで構成されたスロットル弁開度センサ
10が設けられている。このスロットル弁開度センサ1
0は、スロットル弁5の開度(以下「スロットル弁開
度」という)θTHを検出して、その検出信号をECU
2に送る。
【0086】さらに、吸気管4のスロットル弁5よりも
下流側には、吸気管内絶対圧センサ11が設けられてい
る。この吸気管内絶対圧センサ11は、例えば半導体圧
力センサなどで構成され、吸気管4内の吸気管内絶対圧
PBAを検出し、その検出信号をECU2に出力する。
【0087】また、吸気管4は、インテークマニホール
ド4aの4つの分岐部4bを介して4つの気筒#1〜#
4にそれぞれ接続されている。各分岐部4bには、各気
筒の図示しない吸気ポートの上流側に、インジェクタ6
が取り付けられている。各インジェクタ6は、エンジン
3の運転時に、ECU2からの駆動信号によって、その
開弁時間である最終燃料噴射量TOUTおよび噴射タイ
ミングが制御される。
【0088】一方、エンジン3の本体には、例えばサー
ミスタなどで構成された水温センサ12が取り付けられ
ている。水温センサ12は、エンジン3のシリンダブロ
ック内を循環する冷却水の温度であるエンジン水温TW
を検出し、その検出信号をECU2に出力する。
【0089】また、エンジン3のクランクシャフト(図
示せず)には、クランク角センサ13が設けられてい
る。このクランク角センサ13は、クランクシャフトの
回転に伴い、いずれもパルス信号であるCRK信号およ
びTDC信号をECU2に出力する。
【0090】CRK信号は、所定のクランク角(例えば
30゜)ごとに1パルスが出力される。ECU2は、こ
のCRK信号に応じ、エンジン3の回転数(以下「エン
ジン回転数」という)NEを算出する。また、TDC信
号は、各気筒のピストン(図示せず)が吸気行程のTD
C位置よりも若干、手前の所定のクランク角位置にある
ことを表す信号であり、所定クランク角ごとに1パルス
が出力される。
【0091】一方、排気管7(排気通路)のエキゾース
トマニホールド7aよりも下流側には、上流側から順に
第1および第2の触媒装置8a,8b(触媒)が間隔を
存して設けられている。各触媒装置8は、NOx触媒と
3元触媒を組み合わせたものであり、このNOx触媒
は、図示しないが、イリジウム触媒(イリジウムを担持
した炭化ケイ素ウイスカ粉末とシリカの焼成体)をハニ
カム構造の基材の表面に被覆し、その上にペロブスカイ
ト型複酸化物(LaCoO3粉末とシリカの焼成体)を
さらに被覆したものである。触媒装置8は、NOx触媒
による酸化還元作用により、リーンバーン運転時の排気
ガス中のNOxを浄化するとともに、3元触媒の酸化還
元作用により、リーンバーン運転以外の運転時の排気ガ
ス中のCO、HCおよびNOxを浄化する。なお、触媒
装置8は、NOx触媒と3元触媒を組み合わせたものに
限らず、排気ガス中のCO、HCおよびNOxを浄化で
きるものであればよい。例えば、触媒装置8を、ペロブ
スカイト型触媒などの非金属触媒および/または3元触
媒などの金属触媒で構成してもよい。
【0092】これらの第1および第2触媒装置8a,8
bの間には、酸素濃度センサ(以下「O2センサ」とい
う)15が取り付けられている。このO2センサ15
(下流側空燃比センサ)は、ジルコニアおよび白金電極
などで構成され、第1触媒装置8aの下流側の排気ガス
中の酸素濃度に基づく出力VoutをECU2に送る。
このO2センサ15の出力Vout(制御対象の出力)
は、理論空燃比よりもリッチな混合気が燃焼したときに
は、ハイレベルの電圧値(例えば0.8V)となり、混
合気がリーンのときには、ローレベルの電圧値(例えば
0.2V)となるとともに、混合気が理論空燃比付近の
ときには、ハイレベルとローレベルの間の所定の目標値
Vop(例えば0.6V)となる(図2参照)。
【0093】また、第1触媒装置8aよりも上流側のエ
キゾーストマニホールド7aの集合部付近には、LAF
センサ14(上流側空燃比センサ)が取り付けられてい
る。このLAFセンサ14は、O2センサ15と同様の
センサとリニアライザなどの検出回路とを組み合わせる
ことによって構成されており、リッチ領域からリーン領
域までの広範囲な空燃比の領域において排気ガス中の酸
素濃度をリニアに検出し、その酸素濃度に比例する出力
KACTをECU2に送る。この出力KACTは、空燃
比の逆数に比例する当量比として表される。
【0094】次に、図2を参照しながら、第1触媒装置
8aの排気ガスの浄化率とO2センサ15の出力Vou
t(電圧値)との関係について説明する。同図は、第1
触媒装置8aが、長時間の使用により浄化能力が低下し
た劣化状態と、浄化能力の高い未劣化状態の場合におい
て、LAFセンサ14の出力KACTすなわちエンジン
3に供給される混合気の空燃比が理論空燃比の付近で変
化したときの、2つの第1触媒装置8aのHCおよびN
Oxの浄化率と、O2センサ15の出力Voutをそれ
ぞれ測定した結果の一例を示している。同図において、
破線で示すデータはいずれも、第1触媒装置8aが未劣
化状態の場合の測定結果であり、実線で示すデータはい
ずれも、第1触媒装置8aが劣化状態の場合の測定結果
である。また、LAFセンサ14の出力KACTが大き
いほど、混合気の空燃比がよりリッチ側であることを示
している。
【0095】同図に示すように、第1触媒装置8aが劣
化している場合には、未劣化状態の場合と比べて、排気
ガスの浄化能力が低下していることにより、LAFセン
サ14の出力KACTがよりリーン側の値KACT1の
ときに、O2センサ15の出力Voutが目標値Vop
を横切っている。一方、第1触媒装置8aは、その劣化
・未劣化状態にかかわらず、O2センサ15の出力Vo
utが目標値Vopにあるときに、HCおよびNOxを
最も効率よく浄化する特性を有している。したがって、
O2センサ15の出力Voutが目標値Vopになるよ
うに、混合気の空燃比を制御することにより、第1触媒
装置8aによって排気ガスを最も効率よく浄化できるこ
とが判る。このため、後述する空燃比制御では、O2セ
ンサ15の出力Voutが目標値Vopに収束するよう
に、目標空燃比KCMDが制御される。
【0096】さらに、ECU2には、アクセル開度セン
サ16、大気圧センサ17、吸気温センサ18および車
速センサ19などが接続されている。このアクセル開度
センサ16は、車両の図示しないアクセルペダルの踏み
込み量(以下「アクセル開度」という)APを検出し、
その検出信号をECU2に出力する。また、大気圧セン
サ17、吸気温センサ18および車速センサ19はそれ
ぞれ、大気圧PA、吸気温TAおよび車速VPを検出
し、その検出信号をECU2に出力する。
【0097】次に、ECU2(偏差算出手段、制御入力
算出手段、予測値算出手段、出力サンプリング手段、入
力サンプリング手段、サンプリング手段、同定手段、空
燃比算出手段、空燃比制御手段、空燃比サンプリング手
段)について説明する。
【0098】このECU2は、I/Oインターフェー
ス、CPU、RAMおよびROMなどからなるマイクロ
コンピュータから構成されており、前述した各種のセン
サ10〜19の出力に応じて、エンジン3の運転状態を
判別するとともに、ROMに予め記憶された制御プログ
ラムやRAMに記憶されたデータなどに従って、後述す
る適応空燃比制御処理またはマップ検索処理を実行する
ことにより、O2センサ15の出力Voutを目標値V
opに収束させるために、目標空燃比KCMDを算出す
る(なお、この目標空燃比KCMDは、空燃比の逆数に
比例する当量比として算出される)。さらに、後述する
ように、この目標空燃比KCMDに基づいて、インジェ
クタ6の最終燃料噴射量TOUTを気筒ごとに算出し、
この算出した最終燃料噴射量TOUTに基づいた駆動信
号で、インジェクタ6を駆動することにより、LAFセ
ンサ14の出力KACTを目標空燃比KCMDに一致さ
せるように、混合気の空燃比をフィードバック制御す
る。
【0099】図3に示すように、制御装置1は、目標空
燃比KCMD(制御入力)を算出するADSMコントロ
ーラ20およびPRISMコントローラ21を備えてお
り、両コントローラ20,21はいずれも、具体的に
は、ECU2により構成されている。
【0100】以下、ADSMコントローラ20について
説明する。このADSMコントローラ20は、以下に述
べる適応予測型ΔΣ変調制御(Adaptive prediction De
ltaSigma Modulation Control:以下「ADSM」とい
う)処理の制御アルゴリズムにより、O2センサ15の
出力Voutを目標値Vopに収束させるための目標空
燃比KCMDを算出するものであり、状態予測器22、
オンボード同定器23およびDSMコントローラ24に
より構成されている。なお、このADSM処理の具体的
なプログラムについては、後述する。
【0101】まず、状態予測器22について説明する。
この状態予測器22は、以下に述べる予測アルゴリズム
により、出力偏差VO2の予測値PREVO2を予測
(算出)するものである。本実施形態では、制御対象へ
の制御入力を混合気の目標空燃比KCMDとし、制御対
象の出力をO2センサ15の出力Voutとし、インジ
ェクタ6を含むエンジン3の吸気系から、第1触媒装置
8aを含む排気系の第1触媒装置8aの下流側のO2セ
ンサ15までの系を、制御対象と見なすとともに、この
制御対象を、下式(1)に示すように、離散時間系モデ
ルであるARXモデル(auto-regressive model with e
xogeneous input:外部入力を持つ自己回帰モデル)とし
てモデル化する。
【0102】 VO2(k)=a1・VO2(k-1)+a2・VO2(k-2)+b1・DKCMD(k-dt) ……(1) ここで、VO2は、O2センサ15の出力Voutと前
述した目標値Vopとの偏差(Vout−Vop)であ
る出力偏差を表し、DKCMDは、目標空燃比KCMD
(=φop)と基準値FLAFBASEとの偏差(KC
MD−FLAFBASE)である空燃比偏差を表し、記
号kは、サンプリングにより離散化された時間を表して
いる。この基準値FLAFBASEは、所定の一定値に
設定される。また、a1,a2,b1はモデルパラメー
タを表しており、オンボード同定器23により、後述す
るように逐次同定される。
【0103】さらに、上記式(1)のdtは、目標空燃
比KCMDの混合気がインジェクタ6により吸気系に供
給されてから、O2センサ15の出力Voutに反映さ
れるまでの予測時間を表しており、下式(2)のように
定義される。 dt=d+d'+dd ……(2) ここで、dは、LAFセンサ14からO2センサ15ま
での排気系のむだ時間を、d'は、インジェクタ6から
LAFセンサ14までの空燃比操作系のむだ時間を、d
dは、排気系と空燃比操作系との間の位相遅れ時間をそ
れぞれ表している(なお、後述する適応空燃比制御処理
の制御プログラムでは、ADSM処理とPRISM処理
とに切り換えて目標空燃比KCMDを算出する処理を行
うため、位相遅れ時間dd=0に設定されている)。
【0104】以上のように、制御対象モデルを、出力偏
差VO2の離散データ(時系列データ)および空燃比偏
差DKCMDの離散データで構成した理由は以下によ
る。すなわち、一般に、制御対象モデルでは、制御対象
の入出力と所定値との偏差を、入出力を表す変数として
定義した場合の方が、入出力の絶対値を変数として定義
した場合よりも、モデルパラメータをより正確に同定ま
たは定義できることで、制御対象モデルの動特性を制御
対象の実際の動特性に適合させることができるという事
実が知られている。したがって、本実施形態の制御装置
1のように、制御対象モデルを、出力偏差VO2の離散
データおよび空燃比偏差DKCMDの離散データで構成
することにより、O2センサ15の出力Voutおよび
目標空燃比KCMDの絶対値を変数とする場合と比べ
て、制御対象の実際の動特性に対する制御対象モデルの
動特性の適合性を向上させることができ、それにより予
測値PREVO2の算出精度を向上させることができ
る。
【0105】また、予測値PREVO2は、目標空燃比
KCMDの混合気が吸気系に供給されてから予測時間d
tが経過した後の出力偏差VO2(k+dt)を予測し
た値であり、上記式(1)に基づき、予測値PREVO
2の算出式を導出すると、下式(3)が得られる。 PREVO2(k)≒VO2(k+dt) =a1・VO2(k+dt-1)+a2・VO2(k+dt-2)+b1・DKCMD(k) ……(3)
【0106】この式(3)では、出力偏差VO2(k)
の未来値に相当するVO2(k+dt−1),VO2
(k+dt−2)の算出が必要となり、実際にプログラ
ム化するのは困難である。そのため、マトリクスA、B
を、モデルパラメータa1,a2,b1を用いて図4に
示す式(4)、(5)のように定義するとともに、上式
(3)の漸化式を繰り返し用いることにより、上式
(3)を変形すると、図4に示す式(6)が得られる。
予測アルゴリズムすなわち予測値PREVO2の算出式
として、この式(6)を用いた場合、予測値PREVO
2が、出力偏差VO2および空燃比偏差DKCMDによ
り算出される。
【0107】次に、LAF出力偏差DKACTを、LA
Fセンサ14の出力KACT(=φin)と基準値FL
AFBASEとの偏差(KACT−FLAFBASE)
として定義すると、DKACT(k)=DKCMD(k
−d')の関係が成立するので、この関係を図4の式
(6)に適用すると、図4に示す式(7)が得られる。
【0108】以上の式(6)または式(7)により算出
される予測値PREVO2を用い、後述するように目標
空燃比KCMDを算出することによって、制御対象の入
出力間の応答遅れやむだ時間を適切に補償しながら、目
標空燃比KCMDを算出することができる。特に、予測
アルゴリズムとして、上記式(7)を用いた場合、予測
値PREVO2が、出力偏差VO2、LAF出力偏差D
KACTおよび目標空燃比KCMDにより算出されるの
で、第1触媒装置8aに実際に供給される排気ガスの空
燃比の状態が反映された値として、予測値PREVO2
を算出でき、その算出精度すなわち予測精度を上記式
(6)を用いた場合よりも向上させることができる。ま
た、式(7)を用いた場合において、d'≦1と見なせ
るときには、空燃比偏差DKCMDを用いることなく、
出力偏差VO2およびLAF出力偏差DKACTのみに
より、予測値PREVO2を算出できる。本実施形態で
は、LAFセンサ14がエンジン3に設けられているの
で、予測アルゴリズムとして上記式(7)を採用する。
【0109】なお、前述した式(1)の制御対象モデル
は、DKACT(k)=DKCMD(k−d')の関係
を適用することにより、出力偏差VO2およびLAF出
力偏差DKACTを変数とするモデルとして定義するこ
とも可能である。
【0110】次に、オンボード同定器23(同定手段)
について説明する。このオンボード同定器23は、以下
に述べる逐次型同定アルゴリズムにより、前述した式
(1)のモデルパラメータa1,a2,b1を同定(算
出)するものである。具体的には、図5に示す(8),
(9)により、モデルパラメータのベクトルθ(k)を
算出する。同図の式(8)において、KP(k)は、ゲ
イン係数のベクトルであり、ide_f(k)は同定誤
差フィルタ値である。また、式(9)におけるθ(k)
Tは、θ(k)の転置行列を表し、a1'(k)、a2'
(k)およびb1'(k)は、後述するリミット処理を施す
前のモデルパラメータを表している。なお、以下の説明
では、「ベクトル」という表記を適宜、省略する。
【0111】上記式(8)の同定誤差フィルタ値ide
_f(k)は、図5に示す式(11)〜(13)により
算出される同定誤差ide(k)に、図5の式(10)
に示す移動平均フィルタリング処理を施した値である。
図5の式(10)のnは、移動平均フィルタリング処理
のフィルタ次数(1以上の整数)を表しており、式(1
2)のVO2HAT(k)は、出力偏差VO2の同定値
を表している。
【0112】この同定誤差フィルタ値ide_f(k)
を用いる理由は以下による。すなわち、本実施形態の制
御対象は、目標空燃比KCMDを制御入力とし、O2セ
ンサ15の出力Voutを制御対象の出力とするもので
あり、その周波数特性としてはローパス特性を有してい
る。このようなローパス特性を有する制御対象では、オ
ンボード同定器23の同定アルゴリズム、具体的には後
述する重み付き最小2乗法アルゴリズムの周波数重み特
性に起因して、制御対象の高周波特性が強調された状態
で、モデルパラメータが同定されるため、制御対象モデ
ルのゲイン特性が制御対象の実際のゲイン特性よりも低
くなる傾向を示す。その結果、制御装置1によりADS
M処理またはPRISM処理が実行された際、オーバー
ゲイン状態になることで、制御系が発散状態になり、不
安定になる可能性がある。
【0113】したがって、本実施形態では、重み付き最
小2乗法アルゴリズムの周波数重み特性を適切に補正
し、制御対象モデルのゲイン特性を、制御対象の実際の
ゲイン特性に一致させるために、上記同定誤差ide
(k)に移動平均フィルタリング処理を施した同定誤差
フィルタ値ide_f(k)を用いるとともに、後述す
るように、移動平均フィルタリング処理のフィルタ次数
nを、排気ガスボリュームAB_SVに応じて設定して
いる。
【0114】さらに、前述した図5の式(8)のゲイン
係数のベクトルKP(k)は、図5の式(14)により
算出される。この式(14)のP(k)は、図5の式
(15)で定義される3次の正方行列である。
【0115】以上のような同定アルゴリズムでは、式
(15)の重みパラメータλ1、λ2の設定により、以
下の4つの同定アルゴリズムのうちの1つが選択され
る。すなわち、 λ1=1,λ2=0 ;固定ゲインアルゴリズム λ1=1,λ2=1 ;最小2乗法アルゴリズム λ1=1,λ2=λ ;漸減ゲインアルゴリズム λ1=λ,λ2=1 ;重み付き最小2乗法アルゴリズ
ム ただし、λは、0<λ<1に設定される所定値。
【0116】本実施形態では、これらの4つの同定アル
ゴリズムのうちの重み付き最小2乗法アルゴリズムを採
用する。これは、重みパラメータλ1の値をエンジン3
の運転状態、具体的には排気ガスボリュームAB_SV
に応じて設定することにより、同定精度と、モデルパラ
メータの最適値への収束速度とを適切に設定できること
による。例えば、低負荷運転状態のときには、それに応
じて重みパラメータλ1の値を値1に近い値に設定する
ことで、すなわち最小2乗法アルゴリズムに近いアルゴ
リズムに設定することで、良好な同定精度を確保できる
とともに、高負荷運転状態のときには、それに応じて重
みパラメータλ1の値を低負荷運転状態のときよりも小
さい値に設定することにより、モデルパラメータを迅速
に最適値に収束させることができる。以上のように、重
みパラメータλ1の値を排気ガスボリュームAB_SV
に応じて設定することにより、同定精度と、モデルパラ
メータの最適値への収束速度とを適切に設定することが
でき、それにより、触媒装置8a,8bによる排気ガス
浄化率を向上させることができる。
【0117】以上の式(8)〜(15)の同定アルゴリ
ズムにおいて、前述したDKACT(k)=DKCMD
(k−d')の関係を適用すると、図6に示す式(1
6)〜(23)の同定アルゴリズムが得られる。本実施
形態では、LAFセンサ14がエンジン3に設けられて
いるので、これらの式(16)〜(23)を用いる。こ
れらの式(16)〜(23)を用いた場合、前述した理
由により、モデルパラメータを、第1触媒装置8aに実
際に供給される排気ガスの空燃比の状態がより反映され
た値として同定することができ、それにより、上記式
(8)〜(15)の同定アルゴリズムを用いた場合より
も、モデルパラメータの同定精度を向上させることがで
きる。
【0118】また、このオンボード同定器23では、以
上の同定アルゴリズムにより算出されたモデルパラメー
タa1'(k)、a2'(k)およびb1'(k)に、後述する
リミット処理を施すことにより、モデルパラメータa1
(k)、a2(k)およびb1(k)が算出される。さらに、
前述した状態予測器22では、このようにリミット処理
を施した後のモデルパラメータa1(k)、a2(k)およ
びb1(k)に基づき、予測値PREVO2が算出され
る。
【0119】次に、DSMコントローラ24について説
明する。このDSMコントローラ24は、ΔΣ変調アル
ゴリズムを応用した制御アルゴリズムにより、状態予測
器22で算出された予測値PREVO2に基づき、制御
入力φop(k)(=目標空燃比KCMD)を生成(算
出)するとともに、これを制御対象に入力することによ
り、制御対象の出力としてのO2センサ15の出力Vo
utを目標値Vopに収束させるように制御するもので
ある。
【0120】まず、一般的なΔΣ変調アルゴリズムにつ
いて説明する。図7は、ΔΣ変調アルゴリズムを適用し
たコントローラ26により、制御対象27を制御する制
御系の構成を示している。同図に示すように、このコン
トローラ26では、差分器26aにより、参照信号r
(k)と遅延素子26bで遅延されたDSM信号u(k−
1)との偏差として偏差信号δ(k)が生成される。次
に、積分器26cにより、偏差積分値σd(k)が、偏差
信号δ(k)と遅延素子26dで遅延された偏差積分値σ
d(k−1)との和の信号として生成される。次いで、量
子化器26e(符号関数)により、DSM信号u(k)
が、この偏差積分値σd(k)を符号化した信号として生
成される。そして、以上のように生成されたDSM信号
u(k)が制御対象27に入力されることにより、出力信
号 y(k)が制御対象27から出力される。
【0121】以上のΔΣ変調アルゴリズムは、以下の数
式(24)〜(26)で表される。 δ(k)=r(k)−u(k−1) ……(24) σd(k)=σd(k−1)+δ(k) ……(25) u(k)=sgn(σd(k)) ……(26) ただし、符号関数sgn(σd(k))の値は、σd(k)≧
0のときにはsgn(σd(k))=1となり、σd(k)<
0のときにはsgn(σd(k))=−1となる(なお、
σd(k)=0のときに、sgn(σd(k))=0と設定し
てもよい)。
【0122】図8は、以上の制御系の制御シミュレーシ
ョン結果を示している。同図に示すように、正弦波状の
参照信号r(k)を制御系に入力した場合、DSM信号u
(k)が矩形波状の信号として生成され、これを制御対象
27に入力することにより、参照信号r(k)と異なる振
幅で同じ周波数の、ノイズを有するものの全体として同
様の波形の出力信号y(k)が、制御対象27から出力さ
れる。このように、ΔΣ変調アルゴリズムの特性は、参
照信号r(k)から生成されたDSM信号u(k)を制御対
象27に入力したときに、制御対象27の出力y(k)
が、参照信号r(k)に対して、異なる振幅で同じ周波数
の、全体として同様の波形の信号となるような値とし
て、DSM信号u(k)を生成できるという点にある。言
い換えれば、DSM信号u(k)を、参照信号r(k)が制
御対象27の実際の出力y(k)に再現されるような値と
して、生成(算出)できるという点にある。
【0123】DSMコントローラ24は、このようなΔ
Σ変調アルゴリズムの特性を利用し、O2センサ15の
出力Voutを目標値Vopに収束させるための制御入
力φop(k)を算出するものである。その原理について
説明すると、例えば図9に1点鎖線で示すように、出力
偏差VO2が値0に対して揺らいでいる場合(すなわ
ち、O2センサ15の出力Voutが目標値Vopに対
して揺らいでいる場合)、出力偏差VO2を値0に収束
させる(すなわち出力Voutを目標値Vopに収束さ
せる)には、図9に破線で示す、出力偏差VO2を打ち
消すような逆位相波形の出力偏差VO2*が生じるよう
に、制御入力φop(k)を生成すればよい。
【0124】しかし、前述したように、本実施形態の制
御対象では、制御入力φop(k)としての目標空燃比K
CMDが制御対象に入力されてからO2センサ15の出
力Voutに反映されるまでに、予測時間dt分の時間
遅れが発生するため、現在の出力偏差VO2に基づい
て、制御入力φop(k)を算出した場合の出力偏差VO
#は、図9に実線で示すように、出力偏差VO2*に対
して遅れを生じ、それにより、制御タイミングのずれが
生じてしまう。したがって、これを補償するために、本
実施形態のADSMコントローラ20におけるDSMコ
ントローラ24では、出力偏差VO2の予測値PREV
O2を用いることにより、制御入力φop(k)が、制御
タイミングのずれを生じることなく、現在の出力偏差V
O2を打ち消すような出力偏差(逆位相波形の出力偏差
VO2*と同様の出力偏差)を生じさせる信号として生
成される。
【0125】具体的には、このDSMコントローラ24
では、図10に示すように、反転増幅器24aにより、
参照信号r(k)が、値−1、参照信号用のゲインGd
よび予測値PREVO2(k)を互いに乗算した信号とし
て生成される。次に、差分器24bにより、この参照信
号r(k)と遅延素子24cで遅延されたDSM信号u''
(k−1)との偏差として偏差信号δ(k)が生成される。
【0126】次いで、積分器24dにより、偏差積分値
σd(k)が、偏差信号δ(k)と遅延素子24eで遅延さ
れた偏差積分値σd(k−1)との和の信号として生成さ
れ、次に、量子化器24f(符号関数)により、DSM
信号u''(k)が、この偏差積分値σd(k)を符号化した
値として生成される。そして、増幅器24gにより、増
幅DSM信号u(k)がDSM信号u''(k)を所定のゲイ
ンFdで増幅した値として生成され、次に、加算器24
hにより、この増幅DSM信号u(k)を所定の基準値F
LAFBASEに加算した値として、制御入力φop
(k)が生成される。
【0127】以上のDSMコントローラ24の制御アル
ゴリズムは、以下の式(27)〜(32)で表される。 r(k)=−1・Gd・PREVO2(k) ……(27) δ(k)=r(k)−u''(k−1) ……(28) σd(k)=σd(k−1)+δ(k) ……(29) u''(k)=sgn(σd(k)) ……(30) u(k)=Fd・u''(k) ……(31) φop(k)=FLAFBASE+u(k) ……(32) ここで、Gd,Fdはゲインを表す。また、符号関数sg
n(σd(k))の値は、σd(k)≧0のときにはsgn
(σd(k))=1となり、σd(k)<0のときにはsgn
(σd(k))=−1となる(なお、σd(k)=0のとき
に、sgn(σd(k))=0と設定してもよい)。
【0128】このDSMコントローラ24では、以上の
式(27)〜(32)に示す制御アルゴリズムにより、
前述したように、制御入力φop(k)が、制御タイミン
グのずれを生じることなく、出力偏差VO2を打ち消す
ような出力偏差VO2*を生じさせる値として算出され
る。すなわち、制御入力φop(k)が、O2センサ15
の出力Voutを目標値Vopに収束させることができ
る値として算出される。また、制御入力φop(k)が、
増幅DSM信号u(k)を所定の基準値FLAFBASE
に加算した値として算出されるので、制御入力φop
(k)を値0を中心して正負反転する値だけでなく、基準
値FLAFBASEを中心として増減を繰り返す値とし
て算出できる。これにより、通常のΔΣ変調アルゴリズ
ムと比べて、制御の自由度を高めることができる。
【0129】次に、前記PRISMコントローラ21に
ついて説明する。このPRISMコントローラ21は、
以下に述べるオンボード同定型スライディングモード制
御処理(以下「PRISM処理」という)の制御アルゴ
リズムにより、O2センサ15の出力Voutを目標値
Vopに収束させるための目標空燃比KCMDを算出す
るものであり、状態予測器22、オンボード同定器23
およびスライディングモードコントローラ(以下「SL
Dコントローラ」という)25により構成されている。
なお、このPRISM処理の具体的なプログラムについ
ては後述する。
【0130】このPRISMコントローラ21のうちの
状態予測器22およびオンボード同定器23について
は、既に説明したので、ここではSLDコントローラ2
5についてのみ説明する。このSLDコントローラ25
は、スライディングモード制御アルゴリズムに基づいて
スライディングモード制御を行うものであり、以下、一
般的なスライディングモード制御アルゴリズムについて
説明する。このスライディングモード制御アルゴリズム
では、前述した式(1)の離散時間系モデルを制御対象
モデルとして用いるため、切換関数σは、下式(33)
に示すように、出力偏差VO2の時系列データの線形関
数として設定される。 σ(k)=S1・VO2(k)+S2・VO2(k−1) ……(33) ここで、S1,S2は、−1<(S2/S1)<1の関
係が成立するように設定される所定の係数である。
【0131】一般にスライディングモード制御アルゴリ
ズムでは、切換関数σが2つの状態変数(本実施形態で
は出力偏差VO2の時系列データ)で構成されている場
合、2つの状態変数で構成される位相空間は、これらを
それぞれ縦軸および横軸とする2次元の位相平面となる
ため、この位相平面上において、σ=0を満たす2つの
状態変数の値の組み合わせは、切換直線と呼ばれる直線
上に載ることになる。したがって、制御対象への制御入
力を、2つの状態変数の組み合わせが切換直線上に収束
する(載る)ように適切に決定することにより、2つの
状態変数をいずれも、値0になる平衡点に収束(スライ
ディング)させることができる。さらに、スライディン
グモード制御アルゴリズムでは、切換関数σの設定によ
り、状態変数の動特性、より具体的には収束挙動や収束
速度を指定することができる。例えば、本実施形態のよ
うに、切換関数σが2つの状態変数で構成されている場
合には、切換直線の傾きを値1に近づけると、状態変数
の収束速度が遅くなる一方、値0に近づけると、収束速
度が速くなる。
【0132】本実施形態では、前記式(33)に示すよ
うに、切換関数σが出力偏差VO2の2つの時系列デー
タ、すなわち出力偏差VO2の今回値VO2(k)およ
び前回値VO2(k−1)により構成されているので、
これらの今回値VO2(k)および前回値VO2(k−
1)の組み合わせを切換直線上に収束させるように、制
御対象への制御入力すなわち目標空燃比KCMDを設定
すればよい。具体的には、制御量Usl(k)を、基準
値FLAFBASEとの和が目標空燃比KCMDとなる
値として定義すると、今回値VO2(k)および前回値
VO2(k−1)の組み合わせを切換直線上に収束させ
るための制御量Usl(k)は、適応スライディングモ
ード制御アルゴリズムにより、図11に示す式(34)
のように、等価制御入力Ueq(k)、到達則入力Ur
ch(k)および適応則入力Uadp(k)の総和とし
て設定される。
【0133】この等価制御入力Ueq(k)は、出力偏
差VO2の今回値VO2(k)および前回値VO2(k
−1)の組み合わせを切換直線上に拘束しておくための
ものであり、具体的には、図11に示す式(35)のよ
うに定義される。また、到達則入力Urch(k)は、
外乱やモデル化誤差などにより、出力偏差VO2の今回
値VO2(k)および前回値VO2(k−1)の組み合
わせが切換直線上から外れた際に、これらを切換直線上
に収束させるためのものであり、具体的には、図11に
示す式(36)のように定義される。この式(36)に
おいて、Fはゲインを表す。
【0134】さらに、適応則入力Uadp(k)は、制
御対象の定常偏差、モデル化誤差および外乱の影響を抑
制しながら、出力偏差VO2の今回値VO2(k)およ
び前回値VO2(k−1)の組み合わせを、切換超平面
上に確実に収束させるためのものであり、具体的には、
図11に示す式(37)のように定義される。この式
(37)において、Gはゲインを、ΔTは制御周期をそ
れぞれ表す。
【0135】本実施形態のPRISMコントローラ21
のSLDコントローラ25では、前述したように、出力
偏差VO2に代えて予測値PREVO2を用いるので、
PREVO2(k)≒VO2(k+dt)の関係を適用
することにより、以上の式(33)〜(37)のアルゴ
リズムを、図12に示す式(38)〜(42)に書き換
えて用いる。この式(38)におけるσPREは、予測
値PREVO2を用いたときの切換関数(以下「予測切
換関数」という)の値である。すなわち、このSLDコ
ントローラ25では、以上のアルゴリズムで算出される
制御量Usl(k)を基準値FLAFBASEに加算す
ることによって、目標空燃比KCMDが算出される。
【0136】以下、ECU2により実行される燃料噴射
量の算出処理について、図13を参照しながら説明す
る。なお、以下の説明では、今回値であることを示す記
号(k)を適宜、省略する。図13は、この制御処理の
メインルーチンを示しており、本処理は、TDC信号の
入力に同期して割り込み実行される。この処理では、後
述する適応空燃比制御処理、またはマップ検索処理によ
り算出された目標空燃比KCMDを用いることによっ
て、燃料噴射量TOUTが気筒ごとに算出される。
【0137】まず、ステップ1(図では「S1」と略
す。以下同じ)において、前述した各種のセンサ10〜
19の出力を読み込むとともに、読み込んだデータをR
AM内に記憶する。
【0138】次に、ステップ2に進み、基本燃料噴射量
Timを算出する。この処理では、エンジン回転数NE
および吸気管内絶対圧PBAに応じて、図示しないマッ
プを検索することにより、基本燃料噴射量Timを算出
する。
【0139】次いで、ステップ3に進み、総補正係数K
TOTALを算出する。この総補正係数KTOTAL
は、各種の運転パラメータ(例えば吸気温TAや、大気
圧PA、エンジン水温TW、アクセル開度APなど)に
応じて、各種のテーブルやマップを検索することで各種
の補正係数を算出するとともに、これらの各種の補正係
数を互いに乗算することにより、算出される。
【0140】次に、ステップ4に進み、適応制御フラグ
F_PRISMONの設定処理を実行する。この処理の
内容は図示しないが、具体的には、以下の(A)〜
(F)の条件がいずれも成立しているときには、適応空
燃比制御処理で算出された目標空燃比KCMDを使用す
る条件が成立しているとして、それを表すために、適応
制御フラグF_PRISMONが「1」にセットされ
る。一方、(A)〜(F)の条件のうちの少なくとも1
つが成立していないときには、適応制御フラグF_PR
ISMONが「0」にセットされる。 (A)LAFセンサ14およびO2センサ15がいずれ
も活性化していること。 (B)エンジン3がリーンバーン運転中でないこと。 (C)スロットル弁5が全開状態でないこと。 (D)点火時期の遅角制御中でないこと。 (E)フューエルカット運転中でないこと。 (F)エンジン回転数NEおよび吸気管内絶対圧PBA
がいずれも、所定の範囲囲内の値であること。
【0141】次に、ステップ5に進み、ステップ4で設
定された適応制御フラグF_PRISMONが「1」で
あるか否かを判別する。この判別結果がYESのときに
は、ステップ6に進み、目標空燃比KCMDを、後述す
る適応空燃比制御処理で算出された適応目標空燃比KC
MDSLDに設定する。
【0142】一方、ステップ5の判別結果がNOのとき
には、ステップ7に進み、目標空燃比KCMDをマップ
値KCMDMAPに設定する。このマップ値KCMDM
APは、エンジン回転数NEおよび吸気管内絶対圧PB
Aに応じて、図示しないマップを検索することにより、
算出される。
【0143】以上のステップ6または7に続くステップ
8では、オブザーバフィードバック補正係数#nKLA
Fを気筒ごとに算出する。このオブザーバフィードバッ
ク補正係数#nKLAFは、気筒ごとの実際の空燃比の
ばらつきを補正するためのものであり、具体的には、オ
ブザーバによりLAFセンサ14の出力KACTから気
筒ごとの実際の空燃比を推定し、これらの推定した空燃
比に応じて、PID制御により算出される。なお、この
オブザーバフィードバック補正係数#nKLAFの記号
#nは、気筒の番号#1〜#4を表すものであり、これ
は、後述する要求燃料噴射量#nTCYLおよび最終燃
料噴射量#nTOUTにおいても同様である。
【0144】次いで、ステップ9に進み、フィードバッ
ク補正係数KFBを算出する。このフィードバック補正
係数KFBは、LAFセンサ14の出力KACTを目標
空燃比KCMDに一致させるためのものであり、具体的
には、以下のように算出される。すなわち、LAFセン
サ14の出力KACTと目標空燃比KCMDとの偏差に
応じて、PID制御によりフィードバック係数KLAF
を算出する。また、図示しないSelf Tuning Regulator
型の適応制御器によりフィードバック補正係数KSTR
を算出し、これを目標空燃比KCMDで除算することに
より、フィードバック補正係数kstrを算出する。そ
して、エンジン3の運転状態に応じて、これらの2つの
フィードバック係数KLAFおよびフィードバック補正
係数kstrの一方を、フィードバック補正係数KFB
として設定する。
【0145】次いで、ステップ10に進み、補正目標空
燃比KCMDMを算出する。この補正目標空燃比KCM
DMは、空燃比A/Fの変化による充填効率の変化を補
償するためのものであり、前述したステップ6または7
で算出された目標空燃比KCMDに応じて、図示しない
テーブルを検索することにより算出される。
【0146】次に、ステップ11に進み、以上のように
算出した基本燃料噴射量Tim、総補正係数KTOTA
L、オブザーバフィードバック補正係数#nKLAF、
フィードバック補正係数KFB、および補正目標空燃比
KCMDMを用い、下式(43)により、気筒ごとの要
求燃料噴射量#nTCYLを算出する。 #nTCYL=Tim・KTOTAL・KCMDM・KFB・#nKLAF ……(43)
【0147】次に、ステップ12に進み、要求燃料噴射
量#nTCYLを付着補正することにより、最終燃料噴
射量#nTOUTを算出する。この最終燃料噴射量#n
TOUTは、具体的には、今回の燃焼サイクルでインジ
ェクタ6から噴射された燃料が燃焼室の内壁面に付着す
る割合などを、運転状態に応じて算出し、そのように算
出した割合に基づいて、要求燃料噴射量#nTCYLを
補正することにより、算出される。
【0148】次いで、ステップ13に進み、以上のよう
に算出した最終燃料噴射量#nTOUTに基づく駆動信
号を、対応する気筒のインジェクタ6に出力した後、本
処理を終了する。以上により、LAFセンサ14の出力
KACTを目標空燃比KCMDに一致させるように、混
合気の空燃比がフィードバック制御される。
【0149】次に、図14および図15を参照しなが
ら、ADSM処理およびPRISM処理を含む適応空燃
比制御処理について説明する。この処理は、所定の制御
周期ΔT(30msec)で実行される。また、この処
理では、エンジン3の運転状態に応じて、ADSM処
理、PRISM処理、またはスライディングモード制御
量DKCMDSLDを所定値SLDHOLDに設定する
処理により、目標空燃比KCMDが算出される。
【0150】この処理では、まず、ステップ20におい
て、F/C後判定処理を実行する。この処理の内容は図
示しないが、この処理では、フューエルカット運転中
は、それを表すためにF/C後判定フラグF_AFCが
「1」にセットされ、フューエルカット運転の終了後、
所定時間X_TM_AFCが経過したときには、それを
表すためにF/C後判定フラグF_AFCが「0」にセ
ットされる。
【0151】次に、ステップ21に進み、車速VPに基
づいて、エンジン3を搭載した車両が発進したか否かを
判定する発進判定処理を実行する。図16に示すよう
に、この処理では、まず、ステップ49において、アイ
ドル運転フラグF_IDLEが「1」であるか否かを判
別する。このアイドル運転フラグF_IDLEは、アイ
ドル運転中であるときに「1」に、それ以外のときに
「0」にセットされる。
【0152】この判別結果がYESで、アイドル運転中
であるときには、ステップ50に進み、車速VPが所定
車速VSTART(例えば1km/h)より小さいか否
かを判別する。この判別結果がYESで、停車中である
ときには、ステップ51に進み、ダウンカウント式の第
1発進判定タイマのタイマ値TMVOTVSTを第1所
定時間TVOTVST(例えば3msec)に設定す
る。
【0153】次いで、ステップ52に進み、ダウンカウ
ント式の第2発進判定タイマのタイマ値TMVSTを、
上記第1所定時間TVOTVSTよりも長い第2所定時
間TVST(例えば500msec)に設定する。次い
で、ステップ53,54において、第1および第2発進
フラグF_VOTVST,F_VSTをいずれも「0」
にセットした後、本処理を終了する。
【0154】一方、ステップ49または50の判別結果
がNOのとき、すなわちアイドル運転中でないか、また
は車両が発進したときには、ステップ55に進み、第1
発進判定タイマのタイマ値TMVOTVSTが値0より
大きいか否かを判別する。この判別結果がYESで、ア
イドル運転の終了後または車両の発進後、第1所定時間
TVOTVSTが経過していないときには、第1発進モ
ード中であるとして、ステップ56に進み、それを表す
ために第1発進フラグF_VOTVSTを「1」にセッ
トする。
【0155】一方、ステップ55の判別結果がNOで、
アイドル運転の終了後または車両の発進後、第1所定時
間TVOTVSTが経過したときには、第1発進モード
が終了したとして、ステップ57に進み、第1発進フラ
グF_VOTVSTを「0」にセットする。
【0156】ステップ56または57に続くステップ5
8では、第2発進判定タイマのタイマ値TMVSTが値
0より大きいか否かを判別する。この判別結果がYES
で、アイドル運転の終了後または車両の発進後、第2所
定時間TVSTが経過していないときには、第2発進モ
ード中であるとして、ステップ59に進み、それを表す
ために第2発進フラグF_VSTを「1」にセットした
後、本処理を終了する。
【0157】一方、ステップ59の判別結果がNOで、
アイドル運転の終了後または車両の発進後、第2所定時
間TVSTが経過したときには、第2発進モードが終了
したとして、前記ステップ54を実行した後、本処理を
終了する。
【0158】図14に戻り、ステップ21に続くステッ
プ22では、状態変数の設定処理を実行する。図示しな
いが、この処理では、RAM内に記憶されている、目標
空燃比KCMD、LAFセンサ14の出力KACTおよ
び出力偏差VO2の時系列データをいずれも、1サンプ
リングサイクル分ずつ過去側にシフトさせる。その後、
KCMD、KACTおよびVO2の時系列データの最新
の値と、基準値FLAFBASEと、後述する適応補正
項FLAFADPとに基づき、KCMD、KACTおよ
びVO2の今回値を算出する。
【0159】次に、ステップ23に進み、PRISM/
ADSM処理の実行判定処理を行う。この処理は、PR
ISM処理またはADSM処理の実行条件が成立してい
るか否かを判定するものであり、具体的には、図17に
示すフローチャートのように実行される。
【0160】すなわち、図17のステップ60〜63に
おいて、以下の(G)〜(J)の条件がいずれも成立し
ているときには、PRISM処理またはADSM処理を
実行すべき運転状態にあるとして、それを表すために、
ステップ64で、PRISM/ADSM実行フラグF_
PRISMCALを「1」にセットした後、本処理を終
了する。一方、(G)〜(J)の条件の少なくとも1つ
が成立していないときには、PRISM処理またはAD
SM処理を実行すべき運転状態にないとして、それを表
すために、ステップ65で、PRISM/ADSM実行
フラグF_PRISMCALを「0」にセットした後、
本処理を終了する。 (G)O2センサ15が活性化していること。 (H)LAFセンサ14が活性化していること。 (I)エンジン3がリーンバーン運転中でないこと。 (J)点火時期の遅角制御中でないこと。
【0161】図14に戻り、ステップ23に続くステッ
プ24では、同定器演算の実行判定処理を行う。この処
理は、オンボード同定器23によるパラメータ同定の実
行条件が成立しているか否かを判定するものであり、具
体的には、図18に示すフローチャートのように実行さ
れる。
【0162】すなわち、図18のステップ70および7
1の判別結果がいずれもNOのとき、言い換えれば、ス
ロットル弁開度θTHが全開状態でなく、かつフューエ
ルカット運転中でないときには、パラメータ同定を実行
すべき運転状態であるとして、ステップ72に進み、同
定実行フラグF_IDCALを「1」にセットした後、
本処理を終了する。一方、ステップ70または71の判
別結果がYESのときには、パラメータ同定を実行すべ
き運転状態にないとして、ステップ73に進み、同定実
行フラグF_IDCALを「0」にセットした後、本処
理を終了する。
【0163】図14に戻り、ステップ24に続くステッ
プ25では、各種パラメータ(排気ガスボリュームAB
_SVなど)を算出する。この処理の具体的な内容は、
後述する。
【0164】次に、ステップ26に進み、前記ステップ
23で設定されたPRISM/ADSM実行フラグF_
PRISMCALが「1」であるか否かを判別する。こ
の判別結果がYESで、PRISM処理またはADSM
処理の実行条件が成立しているときには、ステップ27
に進み、前記ステップ24で設定された同定実行フラグ
F_IDCALが「1」であるか否かを判別する。
【0165】この判別結果がYESで、オンボード同定
器23によるパラメータ同定を実行すべき運転状態のと
きには、ステップ28に進み、パラメータ初期化フラグ
F_IDRSETが「1」であるか否かを判別する。こ
の判別結果がNOで、RAMに記憶されているモデルパ
ラメータa1,a2,b1の初期化が不要であるときに
は、後述するステップ31に進む。
【0166】一方、この判別結果がYESで、モデルパ
ラメータa1,a2,b1の初期化が必要であるときに
は、ステップ29に進み、モデルパラメータa1,a
2,b1を、それぞれの初期値に設定した後、それを表
すためにステップ30に進み、パラメータ初期化フラグ
F_IDRSETを「0」にセットする。
【0167】このステップ30または28に続くステッ
プ31では、オンボード同定器23の演算を実行し、モ
デルパラメータa1,a2,b1を同定した後、後述す
る図15のステップ32に進む。このオンボード同定器
23の演算の具体的な内容については、後述する。
【0168】一方、ステップ27の判別結果がNOで、
パラメータ同定を実行すべき運転状態でないときには、
以上のステップ28〜31をスキップして、図15のス
テップ32に進む。ステップ27または31に続くステ
ップ32では、モデルパラメータa1,a2,b1とし
て、同定値または所定値を選択する。この処理の内容は
図示しないが、具体的には、前記ステップ24で設定さ
れた同定実行フラグF_IDCALが「1」のときに
は、モデルパラメータa1,a2,b1をステップ31
で同定された同定値に設定する。一方、同定実行フラグ
F_IDCALが「0」のときには、モデルパラメータ
a1,a2,b1を所定値に設定する。
【0169】次に、ステップ33に進み、後述するよう
に、状態予測器22の演算を実行し、予測値PREVO
2を算出する。その後、ステップ34に進み、後述する
ように、制御量Uslを算出する。
【0170】次いで、ステップ35に進み、SLDコン
トローラ25の安定判別を実行する。この処理の内容は
図示しないが、具体的には、予測切換関数σPREの値
に基づき、SLDコントローラ25によるスライディン
グモード制御が安定状態にあるか否かを判別する。
【0171】次に、ステップ36および37において、
後述するように、SLDコントローラ25およびDSM
コントローラ24により、スライディングモード制御量
DKCMDSLDおよびΔΣ変調制御量DKCMDDS
Mをそれぞれ算出する。
【0172】次いで、ステップ38に進み、後述するよ
うに、SLDコントローラ25により算出されたスライ
ディングモード制御量DKCMDSLD、またはDSM
コントローラ24により算出されたΔΣ変調制御量DK
CMDDSMを用いて、適応目標空燃比KCMDSLD
を算出する。この後、ステップ39に進み、後述するよ
うに、適応補正項FLAFADPを算出した後、本処理
を終了する。
【0173】一方、図14に戻り、前記ステップ26の
判別結果がNOで、PRISM処理およびADSM処理
の実行条件がいずれも成立していないときには、ステッ
プ40に進み、パラメータ初期化フラグF_IDRSE
Tを「1」にセットする。次に、図15のステップ41
に進み、スライディングモード制御量DKCMDSLD
を所定値SLDHOLDにセットする。次いで、前述し
たステップ38,39を実行した後、本処理を終了す
る。
【0174】次に、図19を参照しながら、前述したス
テップ25の各種パラメータを算出する処理について説
明する。この処理では、まず、ステップ80において、
下式(44)により、排気ガスボリュームAB_SV
(空間速度の推定値)を算出する。 AB_SV=(NE/1500)・PBA・X_SVPRA ……(44) ここで、X_SVPRAは、エンジン排気量に基づいて
決定される所定の係数である。
【0175】次に、ステップ81に進み、前述した空燃
比操作系のむだ時間KACT_D(=d')、排気系の
むだ時間CAT_DELAY(=d)および予測時間d
tを算出する。具体的には、ステップ80で算出された
排気ガスボリュームAB_SVに応じて、図20に示す
テーブルを検索することにより、むだ時間KACT_
D,CAT_DELAYをそれぞれ算出するとともに、
これらの和(KACT_D+CAT_DELAY)を予
測時間dtとして設定する。すなわち、この制御プログ
ラムでは、位相遅れ時間ddが値0に設定される。
【0176】このテーブルでは、排気ガスボリュームA
B_SVが大きいほど、むだ時間KACT_D,CAT
_DELAYがより小さい値に設定されている。これ
は、排気ガスボリュームAB_SVが大きいほど、排気
ガスの流速が大きくなることで、むだ時間KACT_
D,CAT_DELAYが短くなることによる。以上の
ように、むだ時間KACT_D,CAT_DELAYお
よび予測時間dtが、排気ガスボリュームに応じて算出
されるので、これらを用いて算出した出力偏差VO2の
予測値PREVO2に基づき、後述する適応目標空燃比
KCMDSLDを算出することにより、制御対象の入出
力間の制御タイミングのずれを解消することができる。
また、モデルパラメータa1,a2,b1が、上記むだ
時間CAT_DELAYを用いて同定されるので、制御
対象モデルの動特性を、制御対象の実際の動特性に適合
させることができ、それにより、制御対象の入出力間の
制御タイミングのずれをさらに解消することができる。
【0177】次に、ステップ82に進み、同定アルゴリ
ズムの重みパラメータλ1,λ2の値を算出する。具体
的には、重みパラメータλ2を値1に設定すると同時
に、重みパラメータλ1を、排気ガスボリュームAB_
SVに応じて、図21に示すテーブルを検索することに
より算出する。
【0178】このテーブルでは、排気ガスボリュームA
B_SVが大きいほど、重みパラメータλ1がより小さ
い値に設定されており、言い換えれば、排気ガスボリュ
ームAB_SVが小さいほど、重みパラメータλ1がよ
り大きくかつ値1により近い値に設定されている。これ
は、排気ガスボリュームAB_SVが大きいほど、言い
換えれば高負荷運転状態であるほど、モデルパラメータ
の同定をより迅速に行う必要があるので、重みパラメー
タλ1をより小さく設定することによって、モデルパラ
メータの最適値への収束速度を高めるためである。これ
に加えて、排気ガスボリュームAB_SVが小さいほ
ど、すなわち低負荷運転状態であるほど、空燃比が変動
しやすくなり、排気ガス浄化率が不安定になりやすいこ
とで、モデルパラメータの良好な同定精度を確保する必
要があるので、重みパラメータλ1を値1に近づける
(最小2乗法アルゴリズムに近づける)ことによって、
モデルパラメータの同定精度をより高めるためである。
【0179】次に、ステップ83に進み、モデルパラメ
ータa1,a2の値を制限するための下限値X_IDA
2Lと、モデルパラメータb1の値を制限するための下
限値X_IDB1Lおよび上限値X_IDB1Hとを、
排気ガスボリュームAB_SVに応じて、図22に示す
テーブルを検索することにより算出する。
【0180】このテーブルでは、下限値X_IDA2L
は、排気ガスボリュームAB_SVが大きいほど、より
大きい値に設定されている。これは、排気ガスボリュー
ムAB_SVの変化に応じたむだ時間の増減に伴い、制
御系が安定状態となるモデルパラメータa1,a2の組
み合わせが変化することによる。また、下限値X_ID
B1Lおよび上限値X_IDB1Hも、排気ガスボリュ
ームAB_SVが大きいほど、より大きい値に設定され
ている。これは、排気ガスボリュームAB_SVが大き
いほど、触媒前空燃比(第1触媒装置8aよりも上流側
の排気ガスの空燃比)がO2センサ15の出力Vout
に及ぼす影響の度合、すなわち制御対象のゲインがより
大きくなることによる。
【0181】次いで、ステップ84に進み、移動平均フ
ィルタリング処理のフィルタ次数nを算出した後、本処
理を終了する。この処理では、フィルタ次数nを、排気
ガスボリュームAB_SVに応じて、図23に示すテー
ブルを検索することにより、算出する。
【0182】このテーブルでは、排気ガスボリュームA
B_SVが大きいほど、フィルタ次数nがより小さい値
に設定されている。これは、以下の理由による。すなわ
ち、前述したように、排気ガスボリュームAB_SVが
変化すると、制御対象の周波数特性、特にゲイン特性が
変化するので、制御対象モデルのゲイン特性を、制御対
象の実際のゲイン特性に一致させるためには、重み付き
最小2乗法アルゴリズムの周波数重み特性を、排気ガス
ボリュームAB_SVに応じて適切に補正する必要があ
る。したがって、移動平均フィルタリング処理のフィル
タ次数nを、上記テーブルのように排気ガスボリューム
AB_SVに応じて設定することにより、排気ガスボリ
ュームAB_SVの変化にかかわらず、一定の同定重み
を同定アルゴリズムにおいて確保できるとともに、制御
対象モデルと制御対象との間で互いのゲイン特性を一致
させることができ、これにより、同定精度を向上させる
ことができる。
【0183】次に、図24を参照しながら、前記ステッ
プ31のオンボード同定器23の演算処理について説明
する。同図に示すように、この処理では、まず、ステッ
プ90において、前述した式(22)より、ゲイン係数
KP(k)を算出する。次に、ステップ91に進み、前
述した式(20)より、出力偏差VO2の同定値VO2
HAT(k)を算出する。
【0184】次いで、ステップ92に進み、前述した式
(18)(19)より、同定誤差フィルタ値ide_f
(k)を算出する。次に、ステップ93に進み、前述し
た式(16)より、モデルパラメータのベクトルθ
(k)を算出した後、ステップ94に進み、モデルパラ
メータのベクトルθ(k)の安定化処理を実行する。こ
の処理については後述する。
【0185】次いで、ステップ95に進み、前述した式
(23)より、正方行列P(k)の次回値P(k+1)
を算出する。この次回値P(k+1)は、次回のループ
での算出において、正方行列P(k)の値として用いら
れる。
【0186】以下、図25を参照しながら、上記ステッ
プ94におけるモデルパラメータのベクトルθ(k)の
安定化処理について説明する。同図に示すように、ま
ず、ステップ100で、3つのフラグF_A1STA
B,F_A2STAB,F_B1STABをいずれも
「0」にセットする。
【0187】次に、ステップ101に進み、後述するよ
うに、a1'&a2'のリミット処理を実行する。次い
で、ステップ102で、後述するように、b1'のリミ
ット処理を実行した後、本処理を終了する。
【0188】以下、図26を参照しながら、上記ステッ
プ101のa1'&a2'のリミット処理について説明す
る。同図に示すように、まず、ステップ110におい
て、前記ステップ93で算出したモデルパラメータの同
定値a2'が、前記図19のステップ83で算出された
下限値X_IDA2L以上であるか否かを判別する。こ
の判別結果がNOのときには、ステップ111に進み、
制御系を安定化させるために、モデルパラメータa2を
下限値X_IDA2Lに設定すると同時に、モデルパラ
メータa2の安定化を実行したことを表すために、フラ
グF_A2STABを「1」にセットする。一方、この
判別結果がYESで、a2'≧X_IDA2Lのときに
は、ステップ112に進み、モデルパラメータa2を同
定値a2'に設定する。
【0189】これらのステップ111または112に続
くステップ113では、前記ステップ93で算出したモ
デルパラメータの同定値a1'が、所定の下限値X_I
DA1L(例えば値−2以上で値0より小さい一定値)
以上であるか否かを判別する。この判別結果がNOのと
きには、ステップ114に進み、制御系を安定化させる
ために、モデルパラメータa1を下限値X_IDA1L
に設定すると同時に、モデルパラメータa1の安定化を
実行したことを表すために、フラグF_A1STABを
「1」にセットする。
【0190】一方、ステップ113の判別結果がYES
のときには、ステップ115に進み、同定値a1'が、
所定の上限値X_IDA1H(例えば値2)以下である
か否かを判別する。この判別結果がYESで、X_ID
A1L≦a1'≦X_IDA1Hのときには、ステップ
116に進み、モデルパラメータa1を同定値a1'に
設定する。一方、この判別結果がNOで、X_IDA1
H<a1'のときには、ステップ117に進み、モデル
パラメータa1を上限値X_IDA1Hに設定すると同
時に、モデルパラメータa1の安定化を実行したことを
表すために、フラグF_A1STABを「1」にセット
する。
【0191】これらのステップ114、116または1
17に続くステップ118では、以上のように算出した
モデルパラメータa1の絶対値と、モデルパラメータa
2との和(|a1|+a2)が、所定の判定値X_A2
STAB(例えば値0.9)以下であるか否かを判別す
る。この判別結果がYESのときには、モデルパラメー
タa1,a2の組み合わせが、制御系の安定性を確保で
きる範囲(図27にハッチングで示す規制範囲)内にあ
るとして、そのまま本処理を終了する。
【0192】一方、ステップ118の判別結果がNOの
ときには、ステップ119に進み、モデルパラメータa
1が、判定値X_A2STABから下限値X_IDA2
Lを減算した値(X_A2STAB−X_IDA2L)
以下であるか否かを判別する。この判別結果がYESの
ときには、ステップ120に進み、モデルパラメータa
2を、判定値X_A2STABからモデルパラメータa
1の絶対値を減算した値(X_A2STAB−|a1
|)に設定すると同時に、モデルパラメータa2の安定
化を実行したことを表すために、フラグF_A2STA
Bを「1」にセットした後、本処理を終了する。
【0193】一方、ステップ119の判別結果がNO
で、a1>(X_A2STAB−X_IDA2L)のと
きには、ステップ121に進み、制御系を安定化させる
ために、モデルパラメータa1を、判定値X_A2ST
ABから下限値X_IDA2Lを減算した値(X_A2
STAB−X_IDA2L)に設定し、モデルパラメー
タa2を下限値X_IDA2Lに設定する。これと同時
に、モデルパラメータa1,a2の安定化を実行したこ
とを表すために、フラグF_A1STAB,F_A2S
TABをいずれも「1」にセットする。その後、本処理
を終了する。
【0194】前述したように、逐次型の同定アルゴリズ
ムでは、制御対象の入出力が定常状態になると、自己励
起条件の不足化に起因して、同定されたモデルパラメー
タの絶対値が増大する、いわゆるドリフト現象が発生し
やすくなることで、制御系が不安定になったり、振動状
態になったりすることがある。また、その安定限界も、
エンジン3の運転状態に応じて変化する。例えば、低負
荷運転状態のときには、排気ガスボリュームAB_SV
が小さくなることで、供給された混合気に対する排気ガ
スの応答遅れやむだ時間などが大きくなり、それによ
り、O2センサ15の出力Voutが振動状態になりや
すい。
【0195】これに対して、以上のa1'&a2'のリミ
ット処理では、モデルパラメータa1,a2の組み合わ
せが、図27にハッチングで示す規制範囲内の値に収ま
るように設定されるとともに、この規制範囲を決定する
下限値X_IDA2Lが、排気ガスボリュームAB_S
Vに応じて設定されるので、この規制範囲をエンジン3
の運転状態の変化、すなわち制御対象の動特性の変化に
伴う安定限界の変化が反映された適切な安定限界の範囲
として設定することができ、そのような規制範囲内に収
まるように規制されたモデルパラメータa1,a2を用
いることにより、上記ドリフト現象の発生を回避でき、
制御系の安定性を確保することができる。これに加え
て、モデルパラメータa1,a2の組み合わせを、制御
系の安定性を確保できる上記規制範囲内の値として設定
することにより、モデルパラメータa1およびモデルパ
ラメータa2を単独で規制した場合における、制御系の
不安定な状態の発生を回避できる。以上により、制御系
の安定性を向上させることができ、排気ガス浄化率を向
上させることができる。
【0196】次に、図28を参照しながら、前記ステッ
プ102のb1'のリミット処理について説明する。同
図に示すように、この処理では、ステップ130におい
て、前記ステップ93で算出されたモデルパラメータの
同定値b1'が、前記図19のステップ83で算出され
た下限値X_IDB1L以上であるか否かを判別する。
【0197】この判別結果がYESで、b1'≧X_I
DB1Lのときには、ステップ1311に進み、モデル
パラメータの同定値b1'が、前記図19のステップ8
3で算出された上限値X_IDB1H以下であるか否か
を判別する。この判別結果がYESで、X_IDB1L
≦b1'≦X_IDB1Hのときには、ステップ132
に進み、モデルパラメータb1を同定値b1'に設定し
た後、本処理を終了する。
【0198】一方、ステップ131の判別結果がNO
で、b1'>X_IDB1Hのときには、ステップ13
3に進み、モデルパラメータb1を上限値X_IDB1
Hに設定すると同時に、それを表すためにフラグF_B
1LMTを「1」にセットした後、本処理を終了する。
【0199】一方、ステップ130の判別結果がNO
で、b1'<X_IDB1Lのときには、ステップ13
4に進み、モデルパラメータb1を下限値X_IDB1
Lに設定すると同時に、それを表すためにフラグF_B
1LMTを「1」にセットした後、本処理を終了する。
【0200】以上のb1'のリミット処理を実行するこ
とにより、モデルパラメータb1を、X_IDB1L以
上かつX_IDB1H以下の規制範囲内の値に制限する
ことができ、それにより、逐次型の同定アルゴリズムに
よるドリフト現象の発生を回避できる。さらに、前述し
たように、これらの上下限値X_IDB1H,X_ID
B1Lが、排気ガスボリュームAB_SVに応じて設定
されるので、規制範囲をエンジン3の運転状態の変化、
すなわち制御対象の動特性の変化に伴う安定限界の変化
が反映された適切な安定限界の範囲として、設定するこ
とができ、そのような規制範囲内に規制されたモデルパ
ラメータb1を用いることにより、制御系の安定性を確
保することができる。以上により、制御系の安定性を向
上させることができ、排気ガス浄化率を向上させること
ができる。
【0201】次に、図29を参照しながら、前述したス
テップ33の状態予測器22の演算処理について説明す
る。この処理では、まず、ステップ140において、前
述した式(7)の行列要素α1,α2,βi,βjを算
出する。次いで、ステップ141に進み、ステップ14
0で算出した行列要素α1,α2,βi,βjを式
(7)に適用することにより、出力偏差VO2の予測値
PREVO2を算出した後、本処理を終了する。
【0202】次に、図30を参照しながら、前述したス
テップ34の制御量Uslを算出する処理について説明
する。この処理では、まず、ステップ150において、
前述した図12の式(38)により、予測切換関数σP
REを算出する。
【0203】次に、ステップ151に進み、予測切換関
数σPREの積算値SUMSIGMAを算出する。この
処理では、図31に示すように、まず、ステップ160
において、下記の3つの条件(L)〜(N)のうちの少
なくとも1つが成立しているか否かを判別する。 (L)適応制御フラグF_PRISMONが「1」であ
ること。 (M)後述する積算値保持フラグF_SS_HOLDが
「0」であること。 (N)後述するADSM実行フラグF_KOPRが
「0」であること。
【0204】このステップ160の判別結果がYESの
とき、すなわち積算値SUMSIGMAの算出条件が成
立しているときには、ステップ161に進み、積算値S
UMSIGMAの今回値SUMSIGMA(k)を、前
回値SUMSIGMA(k−1)に、制御周期ΔTと予
測切換関数σPREとの積を加算した値[SUMSIG
MA(k−1)+ΔT・σPRE]に設定する。
【0205】次いで、ステップ162に進み、ステップ
161で算出した今回値SUMSIGMA(k)が所定
の下限値SUMSLより大きいか否かを判別する。この
判別結果がYESのときには、ステップ162に進み、
今回値SUMSIGMA(k)が所定の上限値SUMS
Hより小さいか否かを判別する。この判別結果がYES
で、SUMSL<SUMSIGMA(k)<SUMSH
のときには、そのまま本処理を終了する。
【0206】一方、ステップ163の判別結果がNO
で、SUMSIGMA(k)≧SUMSHのときには、
ステップ164に進み、今回値SUMSIGMA(k)
を上限値SUMSHに設定した後、本処理を終了する。
一方、ステップ162の判別結果がNOで、SUMSI
GMA(k)≦SUMSLのときには、ステップ165
に進み、今回値SUMSIGMA(k)を下限値SUM
SLに設定した後、本処理を終了する。
【0207】一方、ステップ160の判別結果がNOの
とき、すなわち3つの条件(L)〜(N)がいずれも不
成立で、積算値SUMSIGMAの算出条件が不成立で
あるときには、ステップ166に進み、今回値SUMS
IGMA(k)を前回値SUMSIGMA(k−1)に
設定する。すなわち、積算値SUMSIGMAをホール
ドする。この後、本処理を終了する。
【0208】図30に戻り、ステップ151に続くステ
ップ152〜154において、前述した図12の式(4
0)〜(42)により、等価制御入力Ueq、到達則入
力Urchおよび適応則入力Uadpをそれぞれ算出す
る。
【0209】次に、ステップ155に進み、これらの等
価制御入力Ueq、到達則入力Urchおよび適応則入
力Uadpの和を、制御量Uslとして設定した後、本
処理を終了する。
【0210】次に、図32,33を参照しながら、前述
した図15のステップ36のスライディングモード制御
量DKCMDSLDの算出処理について説明する。この
処理では、まず、ステップ170において、制御量Us
lのリミット値算出処理を実行する。この処理では、そ
の詳細は説明は省略するが、前述したステップ35のコ
ントローラの安定判別処理の判別結果と、後述する制御
量Uslの適応上下限値Usl_ah,Usl_alと
に基づいて、非アイドル運転用の上下限値Usl_ah
f,Usl_alfと、アイドル運転用の上下限値Us
l_ahfi,Usl_alfiとをそれぞれ算出す
る。
【0211】次いで、ステップ171に進み、アイドル
運転フラグF_IDLEが「0」であるか否かを判別す
る。この判別結果がYESで、アイドル運転中でないと
きには、ステップ171に進み、前述した図30の処理
で算出された制御量Uslが、非アイドル運転用の下限
値Usl_alf以下であるか否かを判別する。
【0212】この判別結果がNOで、Usl>Usl_
alfのときには、ステップ173に進み、制御量Us
lが非アイドル運転用の上限値Usl_ahf以上であ
るか否かを判別する。この判別結果がNOで、Usl_
alf<Usl<Usl_ahfのときには、ステップ
174に進み、スライディングモード制御量DKCMD
SLDを制御量Uslに設定すると同時に、積算値保持
フラグF_SS_HOLDを「0」にセットする。
【0213】次いで、ステップ175に進み、適応下限
値の今回値Usl_al(k)を、前回値Usl_al
(k−1)に所定の減少側値X_AL_DECを加算し
た値[Usl_al(k−1)+X_AL_DEC]に
設定すると同時に、適応上限値の今回値Usl_ah
(k)を、前回値Usl_ah(k−1)から所定の減
少側値X_AL_DECを減算した値[Usl_al
(k−1)−X_AL_DEC]に設定した後、本処理
を終了する。
【0214】一方、ステップ173の判別結果がYES
で、Usl≧Usl_ahfのときには、ステップ17
6に進み、スライディングモード制御量DKCMDSL
Dを非アイドル運転用の適応上限値Usl_ahfに設
定すると同時に、積算値保持フラグF_SS_HOLD
を「1」にセットする。
【0215】次いで、ステップ177に進み、始動後タ
イマのタイマ値TMACRが所定時間X_TMAWAS
Tより小さいこと、またはF/C後判定フラグF_AF
Cが「1」であることが成立しているか否かを判別す
る。この始動後タイマは、エンジン3の始動後の経過時
間を計時するアップカウント式のタイマである。
【0216】この判別結果がYESのとき、すなわち、
エンジン始動後、所定時間X_TMAWASTが経過し
ていないか、またはフューエルカット運転の終了後、所
定時間X_TM_AFCが経過していないときには、そ
のまま本処理を終了する。
【0217】一方、ステップ177の判別結果がNOの
とき、すなわち、エンジン始動後、所定時間X_TMA
WASTが経過し、かつフューエルカット運転の終了
後、所定時間X_TM_AFCが経過したときには、ス
テップ178に進み、適応下限値の今回値Usl_al
(k)を、前回値Usl_al(k−1)に減少側値X
_AL_DECを加算した値[Usl_al(k−1)
+X_AL_DEC]に設定すると同時に、適応上限値
の今回値Usl_ah(k)を、前回値Usl_ah
(k−1)に所定の増大側値X_AL_INCを加算し
た値[Usl_ah(k−1)+X_AL_INC]に
設定した後、本処理を終了する。
【0218】一方、ステップ172の判別結果がYES
で、Usl≦Usl_alfのときには、ステップ17
9に進み、スライディングモード制御量DKCMDSL
Dを非アイドル運転用の適応下限値Usl_alfに設
定すると同時に、積算値保持フラグF_SS_HOLD
を「1」にセットする。
【0219】次いで、ステップ180に進み、第2発進
フラグF_VSTが「1」であるか否かを判別する。こ
の判別結果がYESで、車両の発進後、第2所定時間T
VSTが経過しておらず、第2発進モード中であるとき
には、そのまま本処理を終了する。
【0220】一方、ステップ180の判別結果がNO
で、車両の発進後、第2所定時間TVSTが経過し、第
2発進モードが終了したときには、ステップ181に進
み、適応下限値の今回値Usl_al(k)を、前回値
Usl_al(k−1)から増大側値X_AL_INC
を減算した値[Usl_al(k−1)−X_AL_I
NC]に設定すると同時に、適応上限値の今回値Usl
_ah(k)を、前回値Usl_ah(k−1)から減
少側値X_AL_DECを減算した値[Usl_ah
(k−1)−X_AL_DEC]に設定する。その後、
本処理を終了する。
【0221】一方、ステップ171の判別結果がNO
で、アイドル運転中であるときには、図33のステップ
182に進み、制御量Uslが、アイドル運転用の下限
値Usl_alfi以下であるか否かを判別する。この
判別結果がNOで、Usl>Usl_alfiのときに
は、ステップ183に進み、制御量Uslがアイドル運
転用の上限値Usl_ahfi以上であるか否かを判別
する。
【0222】この判別結果がNOで、Usl_alfi
<Usl<Usl_ahfiのときには、ステップ18
4に進み、スライディングモード制御量DKCMDSL
Dを制御量Uslに設定すると同時に、積算値保持フラ
グF_SS_HOLDを「0」にセットした後、本処理
を終了する。
【0223】一方、ステップ183の判別結果がYES
で、Usl≧Usl_ahfiのときには、ステップ1
85に進み、スライディングモード制御量DKCMDS
LDをアイドル運転用の上限値Usl_ahfiに設定
すると同時に、積算値保持フラグF_SS_HOLDを
「1」にセットした後、本処理を終了する。
【0224】一方、ステップ182の判別結果がYES
で、Usl≦Usl_alfiのときには、ステップ1
86に進み、スライディングモード制御量DKCMDS
LDをアイドル運転用の下限値Usl_alfiに設定
すると同時に、積算値保持フラグF_SS_HOLDを
「1」にセットした後、本処理を終了する。
【0225】次に、図34を参照しながら、前述した図
15のステップ37のΔΣ変調制御量DKCMDDSM
を算出する処理について説明する。同図に示すように、
この処理では、まず、ステップ190において、RAM
に記憶されている、前回のループで算出されたDSM信
号値の今回値DSMSGNS(k)[=u''(k)]
を、前回値DSMSGNS(k−1)[=u''(k−
1)]として設定する。
【0226】次に、ステップ191に進み、RAMに記
憶されている、前回のループで算出された偏差積分値の
今回値DSMSIGMA(k)[=σd(k)]を、前
回値DSMSIGMA(k−1)[=σd(k−1)]
として設定する。
【0227】次いで、ステップ192に進み、出力偏差
の予測値PREVO2(k)が値0以上であるか否かを
判別する。この判別結果がYESのときには、エンジン
3が混合気の空燃比をリーン側に変更すべき運転状態に
あるとして、ステップ193に進み、参照信号値用のゲ
インKRDSM(=Gd)を、リーン化用の値KRDS
MLに設定した後、後述するステップ195に進む。
【0228】一方、ステップ192の判別結果がNOの
ときには、エンジン3が混合気の空燃比をリッチ側に変
更すべき運転状態にあるとして、ステップ194に進
み、参照信号値用のゲインKRDSMを、リーン化用の
値KRDSMLよりも大きいリッチ化用の値KRDSM
Rに設定した後、ステップ195に進む。
【0229】このように、リーン化用の値KRDSML
およびリッチ化用の値KRDSMRが互いに異なる値に
設定されている理由は、以下による。すなわち、混合気
の空燃比をリーン側に変更する際には、第1触媒装置8
aのNOx浄化率を確保すべく、リーンバイアスによる
NOx排出量の抑制効果を得るために、リーン化用の値
KRDSMLをリッチ化用の値KRDSMRよりも小さ
い値に設定することで、O2センサ15の出力Vout
の目標値Vopへの収束速度がリッチ側への変更時より
も遅くなるように、空燃比を制御する。一方、混合気の
空燃比をリッチ側に変更する際には、第1および第2触
媒装置8a,8bのNOx浄化率を十分に回復させるた
め、リッチ化用の値KRDSMRをリーン化用の値KR
DSMLよりも大きい値に設定することで、O2センサ
15の出力Voutの目標値Vopへの収束速度がリー
ン側への変更時よりも速くなるように、空燃比を制御す
る。以上により、混合気の空燃比をリッチ側およびリー
ン側に変更する際、良好な排気ガス浄化率を確保するこ
とができる。
【0230】ステップ193または194に続くステッ
プ195では、値−1、参照信号値用のゲインKRDS
Mおよび予測値の今回値PREVO2(k)を互いに乗
算した値から、上記ステップ190で算出したDSM信
号値の前回値DSMSGNS(k−1)を減算した値
[−1・KRDSM・PREVO2(k)−DSMSG
NS(k−1)]を、偏差信号値DSMDELTA[=
δ(k)]として設定する。この処理は、前述した式
(27),(28)に相当する。
【0231】次いで、ステップ196に進み、偏差積分
値の今回値DSMSIGMA(k)を、ステップ191
で算出した前回値DSMSIGMA(k−1)と、ステ
ップ195で算出した偏差信号値DSMDELTAとの
和[DSMSIGMA(k−1)+DSMDELTA]
に設定する。この処理は、前述した式(29)に相当す
る。
【0232】次に、ステップ197〜199において、
ステップ196で算出した偏差積分値の今回値DSMS
IGMA(k)が値0以上のときには、DSM信号値の
今回値DSMSGNS(k)を値1に設定し、偏差積分
値の今回値DSMSIGMA(k)が値0よりも小さい
ときには、DSM信号値の今回値DSMSGNS(k)
を値−1に設定する。以上のステップ197〜199の
処理は、前述した式(30)に相当する。
【0233】次いで、ステップ200において、排気ガ
スボリュームAB_SVに応じて、図35に示すテーブ
ルを検索することにより、DSM信号値用のゲインKD
SM(=Fd)を算出する。同図に示すように、このゲ
インKDSMは、排気ガスボリュームAB_SVが小さ
いほど、より大きな値に設定されている。これは、排気
ガスボリュームAB_SVが小さいほど、すなわちエン
ジン3の運転負荷が小さい状態であるほど、O2センサ
15の出力Voutの応答性が低下するので、それを補
償するためである。このようにゲインKSDMを設定す
ることにより、ΔΣ変調制御量DKCMDDSMを、例
えばオーバーゲイン状態などを回避しながら、エンジン
3の運転状態に応じて適切に算出することができ、それ
により、排気ガス浄化率を向上させることができる。
【0234】なお、このゲインKDSMの算出に用いる
テーブルは、ゲインKDSMが排気ガスボリュームAB
_SVに応じて設定されている上記テーブルに限らず、
エンジン3の運転負荷状態を表すパラメータ(例えば基
本燃料噴射時間Tim)に応じてゲインKDSMが予め
設定されているものであればよい。また、触媒装置8
a,8bの劣化判別器が設けられている場合には、この
劣化判別器で判別された触媒装置8a,8bの劣化度合
が大きいほど、ゲインDSMをより小さい値に補正する
ようにしてもよい。
【0235】次に、ステップ201に進み、ΔΣ変調制
御量DKCMDDSMを、DSM信号値用のゲインKD
SMと、DSM信号値の今回値DSMSGNS(k)と
を互いに乗算した値[KDSM・DSMSGNS
(k)]に設定した後、本処理を終了する。この処理
が、前述した式(31)に相当する。
【0236】次に、図36を参照しながら、前述した図
15のステップ38の適応目標空燃比KCMDSLDを
算出する処理について説明する。同図に示すように、こ
の処理では、まず、ステップ210において、アイドル
運転フラグF_IDLEが「1」であること、およびア
イドル時ADSM実行フラグF_SWOPRIが「1」
であることがいずれも成立しているか否かを判別する。
このアイドル時ADSM実行フラグF_SWOPRI
は、エンジン3がアイドル運転中で、かつADSM処理
を実行すべき運転状態のときに「1」に、それ以外のと
きに「0」にセットされる。
【0237】この判別結果がYESのとき、すなわちエ
ンジン3がアイドル運転中でADSM処理により適応目
標空燃比KCMDSLDを算出すべき運転状態のときに
は、ステップ211に進み、適応目標空燃比KCMDS
LDを、基準値FLAFBASEにΔΣ変調制御量DK
CMDDSMを加算した値[FLAFBASE+DKC
MDDSM]に設定する。この処理が、前述した式(3
2)に相当する。
【0238】次いで、ステップ212に進み、ADSM
処理を実行したことを表すために、ADSM実行済みフ
ラグF_KOPRを「1」に設定した後、本処理を終了
する。
【0239】一方、ステップ210の判別結果がNOの
ときには、ステップ213に進み、触媒/O2センサフ
ラグF_FCATDSMが「1」であるか否かを判別す
る。この触媒/O2センサフラグF_FCATDSM
は、以下の4つの条件(O)〜(R)のうちの少なくと
も1つが成立しているときに「1」に、それ以外は
「0」にセットされる。 (O)第1触媒装置8aの触媒容量が所定値以上である
こと。 (P)第1触媒装置8aの貴金属含有量が所定値以上で
あること。 (Q)LAFセンサ14がエンジン3の排気管7に設け
られていないこと。 (R)O2センサ15が第2触媒装置8bよりも下流に
設けられていること。
【0240】この判別結果がYESのときには、ステッ
プ214に進み、第1発進フラグF_VOTVST、お
よび発進後ADSM実行フラグF_SWOPRVSTが
いずれも「1」であるか否かを判別する。この発進後A
DSM実行フラグF_SWOPRVSTは、車両の発進
後で、かつエンジン3がADSM処理を実行すべき運転
状態のときに「1」に、それ以外のときに「0」にセッ
トされる。
【0241】この判別結果がYESのとき、すなわち車
両の発進後、第1所定時間TVOTVSTが経過し、か
つADSM処理を実行すべき運転状態のときには、前述
したように、ステップ211,212を実行した後、本
処理を終了する。
【0242】一方、ステップ214の判別結果がNOの
ときには、ステップ215に進み、排気ガスボリューム
AB_SVが所定値OPRSVH以下であること、およ
び小排気時ADSM実行フラグF_SWOPRSVが
「1」であることがいずれも成立しているか否かを判別
する。小排気時ADSM実行フラグF_SWOPRSV
は、エンジン3の排気ガスボリュームAB_SVが小さ
い状態で、かつエンジン3がADSM処理を実行すべき
運転状態のときに「1」に、それ以外のときに「0」に
セットされる。
【0243】この判別結果がYESのとき、すなわち排
気ガスボリュームAB_SVが小さく、かつエンジン3
がADSM処理を実行すべき運転状態のときには、前述
したように、ステップ211,212を実行した後、本
処理を終了する。
【0244】一方、ステップ215の判別結果がNOの
ときには、エンジン3がPRISM処理を実行すべき運
転状態であるとして、ステップ216に進み、適応目標
空燃比KCMDSLDを、基準値FLAFBASEに適
応補正項FLAFADPおよびスライディングモード制
御量DKCMDSLDを加算した値[FLAFBASE
+FLAFADP+DKCMDSLD]に設定する。次
いで、ステップ217に進み、PRISM処理を実行し
たことを表すために、ADSM実行済みフラグF_KO
PRを「0」にセットした後、本処理を終了する。
【0245】一方、前記ステップ213の判別結果がN
Oのとき、すなわち前述した4つの条件(O)〜(R)
がいずれも成立していないときには、ステップ214,
215をスキップし、前述したステップ216,217
を実行した後、本処理を終了する。以上のように、この
適応目標空燃比KCMDSLDの算出処理では、適応目
標空燃比KCMDSLDが、エンジン3の運転状態に応
じて、ADSM処理またはPRISM処理に切り換えて
算出される。
【0246】次に、図37を参照しながら、図15のス
テップ39の適応補正項FLAFADPの算出処理につ
いて説明する。同図に示すように、この処理では、ま
ず、ステップ220において、出力偏差VO2が所定の
範囲(ADL<VO2<ADH)内の値であるか否かを
判別する。この判別結果がYESのとき、すなわち出力
偏差VO2が小さく、O2センサ15の出力Voutが
目標値Vopの近傍にあるときには、ステップ221に
進み、適応則入力Uadpが所定の下限値NRLより小
さいか否かを判別する。
【0247】この判別結果がNOで、Uadp≧NRL
のときには、ステップ222に進み、適応則入力Uad
pが所定の上限値NRHより大きいか否かを判別する。
この判別結果がNOで、NRL≦Uadp≦NRHのと
きには、ステップ223に進み、適応補正項の今回値F
LAFADP(k)を前回値FLAFADP(k−1)
に設定する。すなわち、適応補正項FLAFADPの値
をホールドする。この後、本処理を終了する。
【0248】一方、ステップ222の判別結果がYES
で、Uadp>NRHのときには、ステップ224に進
み、適応補正項の今回値FLAFADP(k)を、前回
値FLAFADP(k−1)に所定の更新値X_FLA
FDLTを加算した値[FLAFADP(k−1)+X
_FLAFDLT]に設定した後、本処理を終了する。
【0249】一方、ステップ221の判別結果がYES
で、Uadp<NRLのときには、ステップ225に進
み、適応補正項の今回値FLAFADP(k)を、前回
値FLAFADP(k−1)から所定の更新値X_FL
AFDLTを減算した値[FLAFADP(k−1)−
X_FLAFDLT]に設定した後、本処理を終了す
る。
【0250】以上のように、第1実施形態の制御装置1
によれば、目標空燃比KCMDを制御入力とし、O2セ
ンサ15の出力Voutを出力とする、位相遅れやむだ
時間などが比較的大きい動特性を有する制御対象におい
て、制御対象の入出力間での制御タイミングのずれを適
切に解消することができ、それにより、制御の安定性お
よび制御性を向上させることができ、排気ガス浄化率を
向上させることができる。
【0251】以下、本発明の第2〜第10実施形態に係
る制御装置について説明する。なお、以下の各実施形態
の説明では、上述した第1実施形態と同じまたは同等の
構成要素については、同一の参照番号を付し、その説明
は適宜、省略するものとする。
【0252】まず、図38を参照しながら、第2実施形
態の制御装置について説明する。同図に示すように、こ
の第2実施形態の制御装置1は、第1実施形態の制御装
置1と比べて、オンボード同定器23のみが異なってい
る。具体的には、第1実施形態のオンボード同定器23
では、KACT、Voutおよびφop(KCMD)に
基づいて、モデルパラメータa1,a2,b1が算出さ
れるのに対して、本実施形態のオンボード同定器23で
は、Voutおよびφopに基づいて、モデルパラメー
タa1,a2,b1が算出される。
【0253】すなわち、このオンボード同定器23で
は、第1実施形態の図6の式(16)〜(23)に示す
同定アルゴリズムに代えて、前述した図5の式(8)〜
(15)に示す同定アルゴリズムにより、モデルパラメ
ータの同定値a1',a2',b1'が算出されるととも
に、これらに前述した図26,28のリミット処理を施
すことにより、モデルパラメータa1,a2,b1が算
出される。このオンボード同定器23の演算処理の具体
的なプログラムは、図示しないが、第1実施形態ものと
ほぼ同様に構成される。以上のような本実施形態の制御
装置1によれば、第1実施形態の制御装置1と同様の効
果を得ることができる。
【0254】次に、図39を参照しながら、第3実施形
態の制御装置について説明する。同図に示すように、こ
の第3実施形態の制御装置1は、第1実施形態の制御装
置1と比べて、状態予測器22のみが異なっている。具
体的には、第1実施形態の状態予測器22では、a1、
a2、b1、KACT、Voutおよびφop(KCM
D)に基づいて、予測値PREVO2が算出されるのに
対して、本実施形態のオンボード同定器23では、a
1、a2、b1、Voutおよびφopに基づいて、予
測値PREVO2が算出される。
【0255】すなわち、この状態予測器22では、第1
実施形態の図4の式(7)に示す予測アルゴリズムに代
えて、同図の式(6)に示す予測アルゴリズムにより、
出力偏差VO2の予測値PREVO2が算出される。こ
の状態予測器22の演算処理の具体的なプログラムは、
図示しないが、第1実施形態のものとほぼ同様に構成さ
れる。この制御装置1によれば、第1実施形態の制御装
置1と同様の効果を得ることができる。
【0256】次に、図40を参照しながら、第4実施形
態の制御装置について説明する。同図に示すように、こ
の第4実施形態の制御装置1は、第1実施形態の制御装
置1と比べると、ADSMコントローラ20、PRIS
Mコントローラ21およびオンボード同定器23に代え
て、スケジュール型DSMコントローラ20A、スケジ
ュール型状態予測スライディングモードコントローラ2
1Aおよびパラメータスケジューラ28を用いること
で、モデルパラメータa1,a2,b1を算出する点の
みが異なっている。
【0257】このパラメータスケジューラ28では、ま
ず、前述した式(44)により、エンジン回転数NEお
よび吸気管内絶対圧PBAに基づいて、排気ガスボリュ
ームAB_SVが算出される。次いで、図41に示すテ
ーブルにより、排気ガスボリュームAB_SVに応じ
て、モデルパラメータa1,a2,b1が算出される。
【0258】このテーブルでは、モデルパラメータa1
は、排気ガスボリュームAB_SVが大きいほど、より
小さい値に設定されており、これとは逆に、モデルパラ
メータa2,b1は、排気ガスボリュームAB_SVが
大きいほど、より大きい値に設定されている。これは、
排気ガスボリュームAB_SVの増大に伴い、制御対象
の出力すなわちO2センサ15の出力Voutが安定化
する一方、排気ガスボリュームAB_SVの減少に伴
い、O2センサ15の出力Voutが振動的になること
による。
【0259】スケジュール型DSMコントローラ20A
は、以上のように算出されたモデルパラメータa1,a
2,b1を用い、前述した第1実施形態と同様のDSM
コントローラ24により目標空燃比KCMDを算出す
る。また、スケジュール型状態予測スライディングモー
ドコントローラ21Aも、以上のように算出されたモデ
ルパラメータa1,a2,b1を用い、前述した第1実
施形態と同様のSLDコントローラ25により目標空燃
比KCMDを算出する。
【0260】この制御装置1によれば、第1実施形態の
制御装置1と同様の効果を得ることができる。これに加
えて、パラメータスケジューラ28を用いることによ
り、オンボード同定器23を用いる場合と比べて、モデ
ルパラメータa1,a2,b1を、より迅速に算出する
ことができる。これにより、制御の応答性を向上させる
ことができ、良好な排気ガス浄化率をより迅速に確保す
ることができる。
【0261】次に、図42を参照しながら、第5実施形
態の制御装置について説明する。この第5実施形態の制
御装置1は、第1実施形態の制御装置1のDSMコント
ローラ24に代えて、SDMコントローラ29を用いる
点のみが異なっている。このSDMコントローラ29
は、ΣΔ変調アルゴリズムを適用した制御アルゴリズム
により、予測値PREVO2(k)に基づいて、制御入
力φop(k)を算出するものである。
【0262】すなわち、同図に示すように、このSDM
コントローラ29では、反転増幅器29aにより、参照
信号r(k)が、値−1、参照信号用のゲインGdおよび
予測値PREVO2(k)を互いに乗算した信号として生
成される。次に、積分器29bにより、参照信号積分値
σdr(k)が、遅延素子29cで遅延された参照信号積
分値σdr(k−1)と参照信号r(k)との和の信号とし
て生成される。一方、積分器29dにより、SDM信号
積分値σdu(k)が、遅延素子29eで遅延されたSD
M信号積分値σdu(k−1)と、遅延素子29jで遅延
されたSDM信号u''(k−1)との和の信号として生成
される。そして、差分器29fにより、参照信号積分値
σdr(k)とSDM信号積分値σdu(k)との偏差信号
δ''(k)が生成される。
【0263】次いで、量子化器29g(符号関数)によ
り、SDM信号u''(k)が、この偏差信号δ''(k)を符
号化した値として生成される。そして、増幅器29hに
より、増幅SDM信号u(k)がSDM信号u''(k)を所
定のゲインFdで増幅した値として生成され、次に、加
算器29iにより、この増幅SDM信号u(k)を所定の
基準値FLAFBASEに加算した値として、制御入力
φop(k)が生成される。
【0264】以上のSDMコントローラ29の制御アル
ゴリズムは、以下の式(45)〜(51)で表される。 r(k)=−1・Gd・PREVO2(k) ……(45) σdr(k)=σdr(k−1)+r(k) ……(46) σdu(k)=σdu(k−1)+u''(k−1) ……(47) δ''(k)=σdr(k)−σdu(k) ……(48) u''(k)=sgn(δ''(k)) ……(49) u(k)=Fd・u''(k) ……(50) φop(k)=FLAFBASE+u(k) ……(51) ここで、Gd,Fdはゲインを表す。また、符号関数sg
n(δ''(k))の値は、δ''(k)≧0のときにはsgn
(δ''(k))=1となり、δ''(k)<0のときにはsg
n(δ''(k))=−1となる(なお、δ''(k)=0のと
きに、sgn(δ''(k))=0と設定してもよい)。
【0265】以上のSDMコントローラ29の制御アル
ゴリズムにおけるΣΔ変調アルゴリズムの特性は、ΔΣ
変調アルゴリズムと同様に、SDM信号u(k)を、これ
を制御対象に入力した際、参照信号r(k)が制御対象の
出力に再現されるような値として、生成(算出)できる
という点にある。すなわち、SDMコントローラ29
は、前述したDSMコントローラ24と同様の制御入力
φop(k)を生成できるという特性を備えている。し
たがって、このSDMコントローラ29を用いる本実施
形態の制御装置1によれば、第1実施形態の制御装置1
と同様の効果を得ることができる。なお、SDMコント
ローラ29の具体的なプログラムは図示しないが、DS
Mコントローラ24とほぼ同様に構成される。
【0266】次に、図43を参照しながら、第6実施形
態の制御装置について説明する。この第6実施形態の制
御装置1は、第1実施形態の制御装置1のDSMコント
ローラ24に代えて、DMコントローラ30を用いる点
のみが異なっている。このDMコントローラ30は、Δ
変調アルゴリズムを適用した制御アルゴリズムにより、
予測値PREVO2(k)に基づいて、制御入力φop
(k)を算出するものである。
【0267】すなわち、同図に示すように、このDMコ
ントローラ30では、反転増幅器30aにより、参照信
号r(k)が、値−1、参照信号用のゲインGdおよび予
測値PREVO2(k)を互いに乗算した信号として生成
される。一方、積分器30bにより、DM信号積分値σ
du(k)が、遅延素子30cで遅延されたDM信号積分
値σdu(k−1)と、遅延素子30hで遅延されたDM
信号u''(k−1)との和の信号として生成される。そし
て、差分器30dにより、参照信号r(k)とDM信号積
分値σdu(k)との偏差信号δ''(k)が生成される。
【0268】次いで、量子化器30e(符号関数)によ
り、DM信号u''(k)が、この偏差信号δ''(k)を符号
化した値として生成される。そして、増幅器30fによ
り、増幅DM信号u(k)がDM信号u''(k)を所定のゲ
インFdで増幅した値として生成され、次に、加算器3
0gにより、この増幅DM信号u(k)を所定の基準値F
LAFBASEに加算した値として、制御入力φop
(k)が生成される。
【0269】以上のDMコントローラ30の制御アルゴ
リズムは、以下の式(52)〜(57)で表される。 r(k)=−1・Gd・PREVO2(k) ……(52) σdu(k)=σdu(k−1)+u''(k−1) ……(53) δ''(k)=r(k)−σdu(k) ……(54) u''(k)=sgn(δ''(k)) ……(55) u(k)=Fd・u''(k) ……(56) φop(k)=FLAFBASE+u(k) ……(57) ここで、Gd,Fdはゲインを表す。また、符号関数sg
n(δ''(k))の値は、δ''(k)≧0のときにはsgn
(δ''(k))=1となり、δ''(k)<0のときにはsg
n(δ''(k))=−1となる(なお、δ''(k)=0のと
きに、sgn(δ''(k))=0と設定してもよい)。
【0270】以上のDMコントローラ30の制御アルゴ
リズムすなわちΔ変調アルゴリズムの特性は、ΔΣ変調
アルゴリズムおよびΣΔ変調アルゴリズムと同様に、D
M信号u(k)を制御対象に入力した際、参照信号r(k)
が制御対象の出力に再現されるような値として、DM信
号u(k)を生成(算出)できるという点にある。すなわ
ち、DMコントローラ30は、前述したDSMコントロ
ーラ24およびSDMコントローラ29と同様の制御入
力φop(k)を生成できるという特性を備えている。
したがって、このDMコントローラ30を用いる本実施
形態の制御装置1によれば、第1実施形態の制御装置1
と同様の効果を得ることができる。なお、DMコントロ
ーラ30の具体的なプログラムは図示しないが、DSM
コントローラ24とほぼ同様に構成される。
【0271】次に、図44および図45を参照しなが
ら、第7実施形態の制御装置について説明する。図44
に示すように、この第7実施形態の制御装置1は、第1
実施形態の制御装置1と比べて、LAFセンサ14がエ
ンジン3に設けられていないとともに、O2センサ15
が第2触媒装置8bよりも下流側に設けられている点の
みが異なっている。
【0272】また、LAFセンサ14を備えていないた
め、この制御装置1では、図45に示すように、オンボ
ード同定器23により、O2センサ15の出力Vout
および制御入力φop(目標空燃比KCMD)に基づい
て、モデルパラメータa1,a2,b1が算出される。
すなわち、このオンボード同定器23では、前述した図
5の式(8)〜(15)に示す同定アルゴリズムによ
り、モデルパラメータの同定値a1',a2',b1'が
算出されるとともに、これらに前述したリミット処理を
施すことにより、モデルパラメータa1,a2,b1が
算出される。
【0273】さらに、状態予測器22により、モデルパ
ラメータa1,a2,b1、O2センサ15の出力Vo
utおよび制御入力φopに基づいて、出力偏差VO2
の予測値PREVO2が算出される。すなわち、この状
態予測器22では、図4の式(6)に示す予測アルゴリ
ズムにより、出力偏差VO2の予測値PREVO2が算
出される。なお、これらの状態予測器22およびオンボ
ード同定器23の演算処理の具体的なプログラムは、図
示しないが、第1実施形態のものとほぼ同様に構成さ
れ、それら以外のプログラムも、第1実施形態のものと
同様に構成される。
【0274】以上のような本実施形態の制御装置1によ
れば、第1実施形態の制御装置1と同様の効果を得るこ
とができる。特に、前述したように、図34のステップ
192〜194において、参照信号値用のゲインKRD
SMを、排気ガスをリーン側に制御する場合と、リッチ
側に制御する場合とで互いに異なる値に設定し、目標空
燃比KCMDの目標値Vopへの収束速度を変更するこ
とにより、本実施形態のようなO2センサ15のみで空
燃比を制御する場合においても、混合気の空燃比をリッ
チ側およびリーン側に変更する際、良好な排気ガス浄化
率を確実に得ることができる。これに加えて、LAFセ
ンサ14を用いることなく、良好な排気ガス浄化率を確
保できるので、その分、製造コストを削減することがで
きる。
【0275】次に、図46を参照しながら、第8実施形
態の制御装置について説明する。同図に示すように、こ
の第8実施形態の制御装置1は、上記第7実施形態の制
御装置1において、ADSMコントローラ20、PRI
SMコントローラ21およびオンボード同定器23を、
前記第4実施形態のスケジュール型DSMコントローラ
20A、スケジュール型状態予測スライディングモード
コントローラ21Aおよびパラメータスケジューラ28
に置き換えたものであり、これらのコントローラ20
A,21Aおよびパラメータスケジューラ28は、第4
実施形態のものと同様に構成されている。この制御装置
1によれば、上記第7実施形態の制御装置1と同様の効
果を得ることができる。これに加えて、パラメータスケ
ジューラ28を用いることにより、オンボード同定器2
3を用いる場合と比べて、モデルパラメータa1,a
2,b1を、より迅速に算出することができる。これに
より、制御の応答性を向上させることができ、良好な排
気ガス浄化率をより迅速に確保することができる。
【0276】次に、図47〜図49を参照しながら、第
9実施形態の制御装置について説明する。図47に示す
ように、この第2実施形態の制御装置1は、第1実施形
態の制御装置1に、2つのデシメーションフィルタ4
0,40(入力サンプリング手段、サンプリング手段)
と、オーバーサンプラ41とを付加したものであり、そ
れに伴い、以下に述べるように、DSMコントローラ2
4による目標空燃比KCMDの算出周期ΔTm(制御入
力算出周期、空燃比算出周期)と、状態予測器22およ
びオンボード同定器23による、予測値PREVO2お
よびモデルパラメータa1,a2,b1の算出周期ΔT
k(偏差算出周期、予測値算出周期)とが、互いに異な
る値に設定されている(図49参照)。
【0277】これらのデシメーションフィルタ40,4
0の一方は、DSMコントローラ24により算出された
目標空燃比KCMDにデシメーション処理を施すととも
に、そのデシメーション値FD(KCMD)と基準値F
LAFBASEとにより、空燃比偏差のデシメーション
値DKCMD’(k)(図49(h)に黒丸で示す値)を
サンプリングする。具体的には、まず、目標空燃比KC
MDを、所定のサンプリング周期ΔTsm(10mse
c)でサンプリングする。次に、これらのサンプリング
値(図49(g)に黒丸および白抜きの丸で示す値)に
ローパスフィルタリング処理を施した後、ΔTsmの4
倍のサンプリング周期ΔTsk(40msec)でダウ
ンサンプリングすることにより、目標空燃比のデシメー
ション値FD(KCMD)(図49(g)に黒丸で示す
値)をサンプリングする。次いで、このデシメーション
値FD(KCMD)と基準値FLAFBASEとの偏差
を、空燃比偏差のデシメーション値DKCMD’(k)と
してサンプリングする。
【0278】また、もう一方のデシメーションフィルタ
40は、LAFセンサ14の出力KACTに上記と同じ
デシメーション処理を施すとともに、そのデシメーショ
ン値FD(KACT)と基準値FLAFBASEとによ
り、LAF出力偏差のデシメーション値DKACT’
(k)(図49(c)に黒丸で示す値)をサンプリングす
る。具体的には、LAFセンサ14の出力KACTを、
所定のサンプリング周期ΔTsm(10msec)でサ
ンプリングし、次に、これらのサンプリング値(図49
(b)に黒丸および白抜きの丸で示す値)にローパスフ
ィルタリング処理を施した後、ΔTsmの4倍のサンプ
リング周期ΔTsk(40msec)でダウンサンプリ
ングすることにより、デシメーション値FD(KAC
T)(図49(b)に黒丸で示す値)をサンプリングす
る。次いで、このデシメーション値FD(KACT)と
基準値FLAFBASEとの偏差を、LAF出力偏差の
デシメーション値DKACT’(k)としてサンプリング
する。なお、以上のデシメーション処理において、サン
プリング周期ΔTskでのダウンサンプリング処理と、
ローパスフィルタリング処理とを、上記とは逆の順序で
実行してもよい。
【0279】さらに、図示しないサンプラにより、O2
センサ15の出力のサンプリング値Vout(k)がサン
プリング周期ΔTskでサンプリングされる(すなわち
図49(a)に黒丸で示す値がサンプリングされる)。
また、これらの3つの値DKCMD’(k),DKAC
T’(k),Vout(k)は、互いに同じタイミングでサ
ンプリングされる(図49参照)。
【0280】オンボード同定器23(同定手段)は、上
記3つの値DKCMD’(k),DKACT’(k),Vo
ut(k)のサンプリングタイミングに同期して、サンプ
リングされたサンプリング値Vout(k)およびデシメ
ーション値DKACT’(k)に基づき、出力偏差VO2
およびモデルパラメータa1,a2,b1を算出(同
定)する(図49(d)参照)。すなわち、オンボード
同定器23によるモデルパラメータa1,a2,b1の
算出周期ΔTkは、サンプリング周期ΔTskと同じ値
(40msec)に設定されている。
【0281】また、状態予測器22(予測値算出手段、
偏差算出手段)は、モデルパラメータa1,a2,b1
の同定値および上記3つの値DKCMD’(k),DKA
CT’(k),Vout(k)に基づき、上記サンプリング
タイミングに同期して、出力偏差VO2および予測値P
REVO2を算出する(図49(e)参照)。すなわ
ち、予測値PREVO2も、上記算出周期ΔTkで算出
される。
【0282】次に、オーバーサンプラ41は、予測値P
REVO2のオーバーサンプル値Fo(PREVO2)
をサンプリングする。具体的には、予測値PREVO2
を、前述したサンプリング周期ΔTsmでオーバーサン
プリングすることにより、オーバーサンプル値Fo(P
REVO2)をサンプリングする(すなわち図49
(f)に白抜きの丸で示す値をサンプリングする)。
【0283】次いで、DSMコントローラ24(制御入
力算出手段、空燃比算出手段)は、このオーバーサンプ
ル値Fo(PREVO2)に基づき、そのサンプリング
タイミングに同期して、目標空燃比KCMDを、前記サ
ンプリング周期ΔTsmと同じ値(10msec)の算
出周期ΔTmで算出する。
【0284】次に、この制御装置1における状態予測器
22の予測アルゴリズムについて説明する。前述したよ
うに、空燃比偏差のデシメーション値DKCMD’(k)
およびLAF出力偏差のデシメーション値DKACT’
(k)はそれぞれ、下式(58),(59)に示すように
定義される。 DKCMD’(k)=Fd(φop(m))−FLAFBASE =Fd(KCMD(m))−FLAFBASE ……(58) DKACT’(k)=Fd(φin(m))−FLAFBASE =Fd(KACT(m))−FLAFBASE ……(59) ここで、記号mは、サンプリング周期ΔTsmでのサン
プリングにより離散化された時間を表し、記号kは、サ
ンプリング周期ΔTskでのサンプリングにより離散化
された時間を表す。また、KCMD(m)およびKAC
T(m)はそれぞれ、目標空燃比KCMDおよびLAF
センサ14の出力KACTの離散データを表す。
【0285】また、制御対象モデルは、下式(60)の
ように定義される。 VO2(k)=a1・VO2(k-1)+a2・VO2(k-2)+b1・DKCMD'(k-dt) ……(60) なお、この式(60)に、DKACT’(k)=DKC
MD’(k−d')の関係を適用することにより、制御
対象モデルを、VO2とDKACT’の関係で表すこと
も可能である。
【0286】さらに、上記式(60)の制御対象モデル
に基づき、第1実施形態と同様の手法を用いると、予測
値PREVO2の算出式は、図48に示す式(61)と
して定義される。
【0287】一方、オンボード同定器23の同定アルゴ
リズムは、図48に示す式(62)〜(69)のように
定義される。
【0288】また、DSMコントローラ24の制御アル
ゴリズムは、以下の式(70)〜(75)で表される。 r(m)=−1・Gd・PREVO2(m) ……(70) δ(m)=r(m)−u''(m−1) ……(71) σd(m)=σd(m−1)+δ(m) ……(72) u''(m)=sgn(σd(m)) ……(73) u(m)=Fd・u''(m) ……(74) φop(m)=KCMD(m)=FLAFBASE+u(m) ……(75) ここで、符号関数sgn(σd(m))の値は、σd(m)≧
0のときにはsgn(σd(m))=1となり、σd(m)<
0のときにはsgn(σd(m))=−1となる(なお、
σd(m)=0のときに、sgn(σd(m))=0と設定し
てもよい)。
【0289】以上のように、本実施形態の制御装置1で
は、状態予測器22およびオンボード同定器23によ
り、予測値PREVO2およびモデルパラメータa1,
a2,b1が、互いに同じ所定のサンプリング周期ΔT
sk(40msec)でサンプリングされた値DKCM
D’,DKACT’,Vout(k)を用い、互いに同じ
所定の算出周期ΔTk(40msec)で算出される。
これに加えて、DSMコントローラ24により、目標空
燃比KCMDが、算出周期ΔTkの1/4の算出周期Δ
Tm(10msec)で算出される。これは以下の理由
による。
【0290】すなわち、DSMコントローラ24では、
算出周期ΔTmが短いほど、目標空燃比KCMDの算出
精度は向上する。これに対して、一般に、O2センサ1
5の出力Voutのパワースペクトルは、0〜5Hzの
周波数域に存在するため、サンプリング周期ΔTskお
よび算出周期ΔTkが20msecよりも短い場合に
は、制御対象の周波数特性が、モデルパラメータa1,
a2,b1の同定値および制御対象モデルに適切に反映
されなくなるとともに、予測値PREVO2の算出値に
も適切に反映されなくなることがある。それにより、モ
デルパラメータa1,a2,b1および予測値PREV
O2の算出精度が低下してしまうことがある。したがっ
て、本実施形態の制御装置1によれば、サンプリング周
期ΔTskおよび算出周期ΔTkが上記の値に設定され
ることにより、制御対象の周波数特性を制御対象モデル
に適切に反映させることができ、モデルパラメータa
1,a2,b1および予測値PREVO2の算出値に、
制御対象の周波数特性を適切に反映させることができ
る。その結果、モデルパラメータa1,a2,b1およ
び予測値PREVO2の算出精度を、第1実施形態の例
よりもさらに向上させることができると同時に、算出周
期ΔTmが算出周期ΔTkよりも短い値に設定されるこ
とにより、目標空燃比KCMDの良好な算出精度も確保
することができる。また、本実施形態の制御装置1によ
れば、前述した第1実施形態の制御装置1と同様の効果
も得ることができる。
【0291】なお、目標空燃比KCMDおよびLAFセ
ンサ14の出力KACTをサンプリング周期ΔTskで
サンプリングするための構成として、第9実施形態で用
いたデシメーションフィルタ40に代えて、これらをサ
ンプリング周期ΔTskでサンプリングするサンプラな
どを用いてもよい。また、DSMコントローラ24に代
えて、前述したSDMコントローラ29またはDMコン
トローラ30を用いてもよい。さらに、各算出周期ΔT
m,ΔTkおよびサンプリング周期ΔTsm,ΔTsk
は、前述した値に限らず、適切な値に設定することが可
能である。例えば、算出周期ΔTmおよびサンプリング
周期ΔTsmを、前述した値(10msec)よりも短
い値に、算出周期ΔTkおよびサンプリング周期ΔTs
kを、前述した値(40msec)よりも長い値にそれ
ぞれ設定してもよい。
【0292】次に、第10実施形態の制御装置について
説明する。この制御装置1は、前述した図44の例と同
様に、LAFセンサ14がエンジン3に設けられていな
いとともに、O2センサ15が第2触媒装置8bよりも
下流側に設けられているエンジン3に適用した例であ
る。そのため、この制御装置1では、図50に示すよう
に、O2センサ15の出力Voutのサンプリング値
と、目標空燃比のデシメーション値DKCMD’とを用
いて、状態予測器22よる予測値PREVO2の算出、
およびオンボード同定器23によるモデルパラメータa
1,a2,b1の算出が行われる。
【0293】また、状態予測器22の予測アルゴリズム
では、前述した式(60)の制御対象モデルに基づき、
予測値PREVO2の算出式は、図51に示す式(7
6)として定義される。さらに、オンボード同定器23
の同定アルゴリズムは、図51に示す式(77)〜(8
4)で表される。また、DSMコントローラ24の制御
アルゴリズムは、前述した式(70)〜(75)で表さ
れる。
【0294】以上の制御装置1によれば、上記第9実施
形態の制御装置1と同様の効果を得ることができる。こ
れに加えて、そのような制御装置1をLAFセンサ14
を用いることなく、構成することができ、その分、製造
コストを削減することができる。
【0295】また、この第10実施形態の制御装置1で
も、第9実施形態の制御装置1と同様に、デシメーショ
ンフィルタ40に代えて、目標空燃比KCMDをサンプ
リング周期ΔTskでサンプリングするサンプラなどを
用いてもよく、DSMコントローラ24に代えて、SD
Mコントローラ29またはDMコントローラ30を用い
てもよい。
【0296】なお、以上の各実施形態は、本発明の制御
装置を内燃機関3の空燃比を制御するものとして構成し
た例であるが、本発明はこれに限らず、他の任意の制御
対象を制御する制御装置に広く適用可能であることは言
うまでもない。また、ADSMコントローラ20および
PRISMコントローラ21を、実施形態のプログラム
に代えて、電気回路により構成してもよい。
【0297】
【発明の効果】以上のように、本発明の制御装置によれ
ば、位相遅れやむだ時間などが比較的大きい制御対象を
制御する場合において、制御対象の入出力間での制御タ
イミングのずれの解消、および制御精度の向上をいずれ
も達成することができる。
【図面の簡単な説明】
【図1】本発明の第1実施形態に係る制御装置およびこ
れを適用した内燃機関の概略構成を示す図である。
【図2】劣化状態および未劣化状態の第1触媒装置を用
いた場合において、LAFセンサの出力KACTに対す
る、両第1触媒装置のHCおよびNOxの浄化率と、O
2センサ15の出力Voutとをそれぞれ測定した結果
の一例を示す図である。
【図3】第1実施形態の制御装置のADSMコントロー
ラおよびPRISMコントローラの構成を示すブロック
図である。
【図4】状態予測器の予測アルゴリズムの数式の一例を
示す図である。
【図5】オンボード同定器の同定アルゴリズムの数式の
一例を示す図である。
【図6】オンボード同定器の同定アルゴリズムの数式の
他の一例を示す図である。
【図7】ΔΣ変調を実行するコントローラおよびこれを
備えた制御系の構成を示すブロック図である。
【図8】図7の制御系の制御結果の一例を示すタイミン
グチャートである。
【図9】第1実施形態のADSMコントローラによる適
応予測型ΔΣ変調制御の原理を説明するためのタイミン
グチャートである。
【図10】ADSMコントローラのうちのDSMコント
ローラの構成を示すブロック図である。
【図11】スライディングモード制御アルゴリズムの数
式を示す図である。
【図12】PRISMコントローラのスライディングモ
ード制御アルゴリズムの数式を示す図である。
【図13】内燃機関の燃料噴射制御処理を示すフローチ
ャートである。
【図14】適応空燃比制御処理を示すフローチャートで
ある。
【図15】図14の続きを示すフローチャートである。
【図16】図14のステップ21における発進判定処理
を示すフローチャートである。
【図17】図14のステップ23におけるPRISM/
ADSM処理の実行判定処理を示すフローチャートであ
る。
【図18】図14のステップ24における同定器演算の
実行判定処理を示すフローチャートである。
【図19】図14のステップ25における各種パラメー
タの算出処理を示すフローチャートである。
【図20】むだ時間CAT_DELAY,KACT_D
の算出に用いるテーブルの一例を示す図である。
【図21】重みパラメータλ1の算出に用いるテーブル
の一例を示す図である。
【図22】モデルパラメータa1,a2,b1の値を制
限するリミット値X_IDA2L,X_IDB1L,X
_IDB1Hの算出に用いるテーブルの一例を示す図で
ある。
【図23】フィルタ次数nの算出に用いるテーブルの一
例を示す図である。
【図24】図14のステップ31における同定器の演算
処理を示すフローチャートである。
【図25】図24のステップ94におけるθ(k)の安
定化処理を示すフローチャートである。
【図26】図25のステップ101におけるa1'&a
2'のリミット処理を示すフローチャートである。
【図27】図26の処理によりa1'&a2'の組み合わ
せが規制される規制範囲を示す図である。
【図28】図25のステップ102におけるb1'のリ
ミット処理を示すフローチャートである。
【図29】図15のステップ33の状態予測器の演算処
理を示すフローチャートである。
【図30】図15のステップ34の制御量Uslの算出
処理を示すフローチャートである。
【図31】図30のステップ151の予測切換関数σP
REの積算値算出処理を示すフローチャートである。
【図32】図15のステップ36のスライディングモー
ド制御量DKCMDSLDの算出処理を示すフローチャ
ートである。
【図33】図32の続きを示すフローチャートである。
【図34】図15のステップ37のΔΣ変調制御量DK
CMDDSMの算出処理を示すフローチャートである。
【図35】KDSMの算出に用いるテーブルの一例を示
す図である。
【図36】図15のステップ38の適応目標空燃比KC
MDSLDの算出処理を示すフローチャートである。
【図37】図15のステップ39の適応補正項FLAF
ADPの算出処理を示すフローチャートである。
【図38】第2実施形態の制御装置の概略構成を示すブ
ロック図である。
【図39】第3実施形態の制御装置の概略構成を示すブ
ロック図である。
【図40】第4実施形態の制御装置の概略構成を示すブ
ロック図である。
【図41】第4実施形態の制御装置のパラメータスケジ
ューラにおいて、モデルパラメータの算出に用いるテー
ブルの一例を示す図である。
【図42】第5実施形態の制御装置のSDMコントロー
ラの概略構成を示すブロック図である。
【図43】第6実施形態の制御装置のDMコントローラ
の概略構成を示すブロック図である。
【図44】第7実施形態に係る制御装置およびこれを適
用した内燃機関の概略構成を示す図である。
【図45】第7実施形態の制御装置の構成を示すブロッ
ク図である。
【図46】第8実施形態の制御装置の構成を示すブロッ
ク図である。
【図47】第9実施形態の制御装置の構成を示すブロッ
ク図である。
【図48】第9実施形態の制御装置における状態予測器
の予測アルゴリズムおよびオンボード同定器の同定アル
ゴリズムの数式の一例を示す図である。
【図49】第9実施形態の制御装置における各算出値の
算出タイミングおよび各サンプリング値のサンプリング
タイミングを示すタイミングチャートである。
【図50】第10実施形態の制御装置の構成を示すブロ
ック図である。
【図51】第10実施形態の制御装置における状態予測
器の予測アルゴリズムおよびオンボード同定器の同定ア
ルゴリズムの数式の一例を示す図である。
【符号の説明】
1 制御装置 2 ECU(偏差算出手段、制御入力算出手段、予測
値算出手段、出力サンプリング手段、入力サンプリング
手段、サンプリング手段、同定手段、空燃比算出手段、
空燃比制御手段、空燃比サンプリング手段) 3 内燃機関 7 排気管(排気通路) 8a 第1触媒装置(触媒) 8b 第2触媒装置(触媒) 14 LAFセンサ(上流側空燃比センサ) 15 酸素濃度センサ(下流側空燃比センサ) 22 状態予測器(予測値算出手段、偏差算出手段) 23 オンボード同定器(同定手段) 24 DSMコントローラ(制御入力算出手段、空燃
比算出手段) 40 デシメーションフィルタ(入力サンプリング手
段、サンプリング手段) ΔTk 算出周期(偏差算出周期、予測値算出周期) ΔTm 算出周期(制御入力算出周期、空燃比算出周
期) ΔTsk サンプリング周期(出力サンプリング周
期、入力サンプリング周期、第1サンプリング周期、同
定周期、空燃比サンプリング周期) ΔTsm サンプリング周期(サンプリング周期、第
2入力サンプリング周期、第2サンプリング周期、第2
空燃比サンプリング周期) KCMD 目標空燃比(制御入力) DKCMD’ 空燃比偏差のデシメーション値(制御入
力を表す値、制御入力と基準値との偏差、目標空燃比を
表す値) DKCMD'(k-dt) 空燃比偏差の離散データ DKCMD'(k-i) 空燃比偏差の離散データ KACT LAFセンサの出力(制御対象に入力され
た制御入力を反映する値) DKACT’ LAF出力偏差のデシメーション値(制
御対象に入力された制御入力を反映する値と基準値との
偏差、制御対象に入力された制御入力を反映する値、上
流側空燃比センサの出力を表す値) DKACT'(k-j) LAF出力偏差の離散データ DKACT'(k-d-dd) LAF出力偏差の離散データ Vout 酸素濃度センサの出力(制御対象の出力) Vop 目標値 VO2 出力偏差(制御対象の出力と目標値との偏
差、制御対象の出力を表す値、下流側空燃比センサの出
力を表す値) VO2(k)〜VO2(k-2) 出力偏差の離散データ PREVO2 出力偏差の予測値 a1 モデルパラメータ a2 モデルパラメータ b1 モデルパラメータ u''(m) DSM信号(第1の中間値を表す信号) u(m) 増幅DSM信号(第1の中間値にゲインを乗
算した値を表す信号、第2の中間値) Fd ゲイン FLAFBASE 基準値(所定値)
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) G05B 11/36 G05B 11/36 M 13/04 13/04 21/02 21/02 Z H03M 3/02 H03M 3/02 Fターム(参考) 3G084 AA03 BA09 BA13 CA03 CA04 DA04 DA10 EA05 EA11 EB06 EB12 EB25 EC01 EC03 FA01 FA02 FA05 FA10 FA11 FA20 FA30 FA33 FA38 3G301 HA01 HA06 JA11 KA06 LB02 MA01 MA11 NA01 NA03 NA04 NA05 NA06 NA08 NA09 NB02 NB13 NC02 ND05 ND18 ND45 NE01 NE06 NE17 NE19 NE21 PA07Z PA09Z PA10Z PA11Z PB08Z PD09A PD09Z PE01Z PE03Z PE08Z PF01Z PF03Z 5H004 GA10 GA18 GA40 GB12 HA13 HB01 HB03 HB04 HB07 HB08 JB07 JB24 KA05 KA33 KA74 KB02 KB04 KB06 KB21 KB28 KC26 KC28 KC35 KC38 KC45 KC46 LA02 LA03 MA11 5J064 BA02 BA03 BB01 BB03 BC08 BC10 BC15 BD01

Claims (38)

    【特許請求の範囲】
  1. 【請求項1】 制御対象の出力と所定の目標値との偏差
    を所定の偏差算出周期で算出する偏差算出手段と、 Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ
    変調アルゴリズムのうちのいずれか1つの変調アルゴリ
    ズムに基づき、前記算出された偏差に応じて、前記制御
    対象の出力を前記目標値に収束させるための前記制御対
    象への制御入力を、前記偏差算出周期よりも短い所定の
    制御入力算出周期で算出する制御入力算出手段と、 を備えることを特徴とする制御装置。
  2. 【請求項2】 前記制御入力算出手段は、前記1つの変
    調アルゴリズムに基づき、前記偏差に応じて、第1の中
    間値を算出するとともに、当該算出した第1の中間値に
    所定のゲインを乗算した値に基づき、前記制御入力を算
    出することを特徴とする請求項1に記載の制御装置。
  3. 【請求項3】 前記制御入力算出手段は、前記1つの変
    調アルゴリズムに基づき、前記偏差に応じて、第2の中
    間値を算出するとともに、当該算出した第2の中間値に
    所定値を加算することにより、前記制御入力を算出する
    ことを特徴とする請求項1に記載の制御装置。
  4. 【請求項4】 予測アルゴリズムに基づき、制御対象の
    出力を表す値の予測値を所定の予測値算出周期で算出す
    る予測値算出手段と、 Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ
    変調アルゴリズムのうちのいずれか1つの変調アルゴリ
    ズムに基づき、前記算出された予測値に応じて、前記制
    御対象の出力を制御するための前記制御対象への制御入
    力を、前記予測値算出周期よりも短い制御入力算出周期
    で算出する制御入力算出手段と、 を備えることを特徴とする制御装置。
  5. 【請求項5】 前記制御入力算出手段は、前記1つの変
    調アルゴリズムに基づき、前記予測値に応じて第1の中
    間値を算出するとともに、当該算出した第1の中間値に
    所定のゲインを乗算した値に基づき、前記制御入力を算
    出することを特徴とする請求項4に記載の制御装置。
  6. 【請求項6】 前記制御入力算出手段は、前記1つの変
    調アルゴリズムに基づき、前記予測値に応じて第2の中
    間値を算出するとともに、当該算出した第2の中間値に
    所定値を加算することにより、前記制御入力を算出する
    ことを特徴とする請求項4に記載の制御装置。
  7. 【請求項7】 前記制御対象の出力を表す値を、前記制
    御入力算出周期よりも長い所定の出力サンプリング周期
    で、離散データとしてサンプリングする出力サンプリン
    グ手段をさらに備え、 前記予測アルゴリズムは、当該サンプリングされた前記
    制御対象の出力を表す値の離散データに応じて、前記予
    測値を算出するアルゴリズムであることを特徴とする請
    求項4ないし6のいずれかに記載の制御装置。
  8. 【請求項8】 前記制御入力を表す値、および前記制御
    対象に入力された制御入力を反映する値の少なくとも一
    方を、前記制御入力算出周期よりも長い所定の入力サン
    プリング周期で、離散データとしてサンプリングする入
    力サンプリング手段をさらに備え、 前記予測アルゴリズムは、当該サンプリングされた、前
    記制御入力を表す値の離散データ、および前記制御対象
    に入力された制御入力を反映する値の離散データの少な
    くとも一方にさらに応じて、前記予測値を算出するアル
    ゴリズムであることを特徴とする請求項7に記載の制御
    装置。
  9. 【請求項9】 前記入力サンプリング手段は、前記制御
    入力を表す値の離散データ、および前記制御対象に入力
    された制御入力を反映する値の離散データの少なくとも
    一方を、前記入力サンプリング周期よりも短い所定の第
    2入力サンプリング周期でサンプリングした値にデシメ
    ーション処理を施すことにより、前記入力サンプリング
    周期でサンプリングすることを特徴とする請求項8に記
    載の制御装置。
  10. 【請求項10】 前記予測アルゴリズムは、前記制御入
    力を表す値および前記制御対象に入力された制御入力を
    反映する値の一方と、前記制御対象の出力を表す値とを
    変数とする制御対象モデルに基づき、前記予測値を算出
    するアルゴリズムであることを特徴とする請求項4ない
    し6のいずれかに記載の制御装置。
  11. 【請求項11】 前記制御入力を表す値と前記制御対象
    に入力された制御入力を反映する値との少なくとも一
    方、および前記制御対象の出力を表す値を、前記制御入
    力算出周期よりも長い所定の第1サンプリング周期で、
    離散データとしてサンプリングするサンプリング手段を
    さらに備え、 前記制御対象モデルは、当該サンプリングされた、前記
    制御入力を表す値の離散データおよび前記制御対象に入
    力された制御入力を反映する値の離散データの一方と、
    前記制御対象の出力を表す値の離散データとを変数とす
    る離散時間系モデルであることを特徴とする請求項10
    に記載の制御装置。
  12. 【請求項12】 前記サンプリング手段は、前記制御入
    力を表す値の離散データ、および前記制御対象に入力さ
    れた制御入力を反映する値の離散データの少なくとも一
    方を、前記第1サンプリング周期よりも短い所定の第2
    サンプリング周期でサンプリングした値にデシメーショ
    ン処理を施すことにより、前記第1サンプリング周期で
    サンプリングすることを特徴とする請求項11に記載の
    制御装置。
  13. 【請求項13】 前記制御対象の出力を表す値の離散デ
    ータに応じて、前記離散時間系モデルのモデルパラメー
    タを同定する同定手段をさらに備えることを特徴とする
    請求項11または12に記載の制御装置。
  14. 【請求項14】 前記同定手段は、前記制御入力を表す
    値の離散データ、および前記制御対象に入力された制御
    入力を反映する値の離散データの一方にさらに応じて、
    前記モデルパラメータを同定することを特徴とする請求
    項13に記載の制御装置。
  15. 【請求項15】 前記同定手段による前記モデルパラメ
    ータの同定周期は、前記制御入力算出周期よりも長い値
    に設定されていることを特徴とする請求項13または1
    4に記載の制御装置。
  16. 【請求項16】 内燃機関の排気通路の触媒よりも下流
    側における排気ガスの空燃比を表す検出信号を出力する
    下流側空燃比センサと、 当該下流側空燃比センサの出力と所定の目標値との偏差
    を所定の偏差算出周期で算出する偏差算出手段と、 Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ
    変調アルゴリズムのうちのいずれか1つの変調アルゴリ
    ズムに基づき、前記算出された偏差に応じて、前記下流
    側空燃比センサの出力を前記目標値に収束させるため
    の、前記内燃機関に供給すべき混合気の目標空燃比を、
    前記偏差算出周期よりも短い所定の空燃比算出周期で算
    出する空燃比算出手段と、 当該算出された目標空燃比に応じて、前記内燃機関に供
    給される混合気の空燃比を制御する空燃比制御手段と、 を備えることを特徴とする制御装置。
  17. 【請求項17】 前記空燃比算出手段は、前記1つの変
    調アルゴリズムに基づき、前記偏差に応じて、第1の中
    間値を算出するとともに、当該算出した第1の中間値に
    所定のゲインを乗算した値に基づき、前記目標空燃比を
    算出することを特徴とする請求項16に記載の制御装
    置。
  18. 【請求項18】 前記空燃比算出手段は、前記1つの変
    調アルゴリズムに基づき、前記偏差に応じて、第2の中
    間値を算出するとともに、当該算出した第2の中間値に
    所定値を加算することにより、前記目標空燃比を算出す
    ることを特徴とする請求項16に記載の制御装置。
  19. 【請求項19】 内燃機関の排気通路の触媒よりも下流
    側における排気ガスの空燃比を表す検出信号を出力する
    下流側空燃比センサと、 予測アルゴリズムに基づき、前記下流側空燃比センサの
    出力を表す値の予測値を所定の予測値算出周期で算出す
    る予測値算出手段と、 Δ変調アルゴリズム、ΔΣ変調アルゴリズムおよびΣΔ
    変調アルゴリズムのうちのいずれか1つの変調アルゴリ
    ズムに基づき、前記算出された予測値に応じて、前記下
    流側空燃比センサの出力を制御するための、前記内燃機
    関に供給すべき混合気の目標空燃比を、前記予測値算出
    周期よりも短い空燃比算出周期で算出する空燃比算出手
    段と、 当該算出された目標空燃比に応じて、前記内燃機関に供
    給される混合気の空燃比を制御する空燃比制御手段と、 を備えることを特徴とする制御装置。
  20. 【請求項20】 前記空燃比算出手段は、前記1つの変
    調アルゴリズムに基づき、前記予測値に応じて、第1の
    中間値を算出するとともに、当該算出した第1の中間値
    に所定のゲインを乗算した値に基づき、前記目標空燃比
    を算出することを特徴とする請求項19に記載の制御装
    置。
  21. 【請求項21】 前記空燃比算出手段は、前記1つの変
    調アルゴリズムに基づき、前記予測値に応じて、第2の
    中間値を算出するとともに、当該算出した第2の中間値
    に所定値を加算することにより、前記内燃機関に供給す
    べき混合気の目標空燃比を算出することを特徴とする請
    求項19に記載の制御装置。
  22. 【請求項22】 前記下流側空燃比センサの出力を表す
    値を、前記空燃比算出周期よりも長い所定の出力サンプ
    リング周期で、離散データとしてサンプリングする出力
    サンプリング手段をさらに備え、 前記予測アルゴリズムは、当該サンプリングされた前記
    下流側空燃比センサの出力を表す値の離散データに応じ
    て、前記予測値を算出するアルゴリズムであることを特
    徴とする請求項19ないし21のいずれかに記載の制御
    装置。
  23. 【請求項23】 前記目標空燃比を表す値を、前記空燃
    比算出周期よりも長い所定の空燃比サンプリング周期
    で、離散データとしてサンプリングする空燃比サンプリ
    ング手段をさらに備え、 前記予測アルゴリズムは、当該サンプリングされた前記
    目標空燃比を表す値の離散データにさらに応じて、前記
    予測値を算出するアルゴリズムであることを特徴とする
    請求項22に記載の制御装置。
  24. 【請求項24】 前記空燃比サンプリング手段は、 前記目標空燃比を表す値の離散データを、前記空燃比サ
    ンプリング周期よりも短い所定の第2空燃比サンプリン
    グ周期でサンプリングした値にデシメーション処理を施
    すことにより、前記空燃比サンプリング周期でサンプリ
    ングすることを特徴とする請求項23に記載の制御装
    置。
  25. 【請求項25】 前記予測アルゴリズムは、前記目標空
    燃比を表す値および前記下流側空燃比センサの出力を表
    す値を変数とする制御対象モデルに基づき、前記予測値
    を算出するアルゴリズムであることを特徴とする請求項
    19ないし21のいずれかに記載の制御装置。
  26. 【請求項26】 前記目標空燃比を表す値および前記下
    流側空燃比センサの出力を表す値を、前記空燃比算出周
    期よりも長い所定の第1サンプリング周期で、離散デー
    タとしてサンプリングするサンプリング手段をさらに備
    え、 前記制御対象モデルは、当該サンプリングされた、前記
    目標空燃比を表す値の離散データと、前記下流側空燃比
    センサの出力を表す値の離散データとを変数とする離散
    時間系モデルであることを特徴とする請求項25に記載
    の制御装置。
  27. 【請求項27】 前記サンプリング手段は、 前記目標空燃比を表す値の離散データを、前記第1サン
    プリング周期よりも短い所定の第2サンプリング周期で
    サンプリングした値にデシメーション処理を施すことに
    より、前記第1サンプリング周期でサンプリングするこ
    とを特徴とする請求項26に記載の制御装置。
  28. 【請求項28】 前記下流側空燃比センサの出力を表す
    値の離散データに応じて、前記離散時間系モデルのモデ
    ルパラメータを同定する同定手段をさらに備えることを
    特徴とする請求項26または27に記載の制御装置。
  29. 【請求項29】 前記同定手段は、前記目標空燃比を表
    す値の離散データにさらに応じて、前記モデルパラメー
    タを同定することを特徴とする請求項28に記載の制御
    装置。
  30. 【請求項30】 前記同定手段による前記モデルパラメ
    ータの同定周期は、前記空燃比算出周期よりも長い値に
    設定されていることを特徴とする請求項28または29
    に記載の制御装置。
  31. 【請求項31】 前記内燃機関の前記排気通路の前記触
    媒よりも上流側における排気ガスの空燃比を表す検出信
    号を出力する上流側空燃比センサと、 前記目標空燃比を表す値と前記上流側空燃比センサの出
    力を表す値との少なくとも一方を、前記空燃比算出周期
    よりも長い所定の空燃比サンプリング周期で、離散デー
    タとしてサンプリングする空燃比サンプリング手段と、
    をさらに備え、 前記予測アルゴリズムは、当該サンプリングされた、前
    記目標空燃比を表す値の離散データ、および前記上流側
    空燃比センサの出力を表す値の離散データの少なくとも
    一方にさらに応じて、前記予測値を算出するアルゴリズ
    ムであることを特徴とする請求項22に記載の制御装
    置。
  32. 【請求項32】 前記サンプリング手段は、 前記目標空燃比を表す値の離散データ、および前記上流
    側空燃比センサの出力を表す値の離散データの少なくと
    も一方を、前記空燃比サンプリング周期よりも短い所定
    の第2空燃比サンプリング周期でサンプリングした値に
    デシメーション処理を施すことにより、前記空燃比サン
    プリング周期でサンプリングすることを特徴とする請求
    項31に記載の制御装置。
  33. 【請求項33】 前記内燃機関の前記排気通路の前記触
    媒よりも上流側における排気ガスの空燃比を表す検出信
    号を出力する上流側空燃比センサをさらに備え、 前記予測アルゴリズムは、前記目標空燃比を表す値およ
    び前記上流側空燃比センサの出力を表す値の一方と、前
    記下流側空燃比センサの出力を表す値とを変数とする制
    御対象モデルに基づき、前記予測値を算出するアルゴリ
    ズムであることを特徴とする請求項19ないし21のい
    ずれかに記載の制御装置。
  34. 【請求項34】 前記目標空燃比を表す値と前記上流側
    空燃比センサの出力を表す値との少なくとも一方、およ
    び前記下流側空燃比センサの出力を表す値を、前記空燃
    比算出周期よりも長い所定の第1サンプリング周期で、
    離散データとしてサンプリングするサンプリング手段を
    さらに備え、 前記制御対象モデルは、当該サンプリングされた、前記
    目標空燃比を表す値の離散データおよび前記上流側空燃
    比センサの出力を表す値の離散データの一方と、前記下
    流側空燃比センサの出力を表す値の離散データとを変数
    とする離散時間系モデルであることを特徴とする請求項
    33に記載の制御装置。
  35. 【請求項35】 前記サンプリング手段は、 前記目標空燃比を表す値の離散データ、および前記上流
    側空燃比センサの出力を表す値の離散データの少なくと
    も一方を、前記第1サンプリング周期よりも短い所定の
    第2サンプリング周期でサンプリングした値にデシメー
    ション処理を施すことにより、前記第1サンプリング周
    期でサンプリングすることを特徴とする請求項34に記
    載の制御装置。
  36. 【請求項36】 前記下流側空燃比センサの出力を表す
    値の離散データに応じて、前記離散時間系モデルのモデ
    ルパラメータを同定する同定手段をさらに備えることを
    特徴とする請求項34または35に記載の制御装置。
  37. 【請求項37】 前記同定手段は、前記目標空燃比を表
    す値の離散データ、および前記上流側空燃比センサの出
    力を表す値の離散データの一方にさらに応じて、前記モ
    デルパラメータを同定することを特徴とする請求項36
    に記載の制御装置。
  38. 【請求項38】 前記同定手段による前記モデルパラメ
    ータの同定周期は、前記空燃比算出周期よりも長い値に
    設定されていることを特徴とする請求項36または37
    に記載の制御装置。
JP2002039094A 2002-02-15 2002-02-15 制御装置 Expired - Fee Related JP3973922B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2002039094A JP3973922B2 (ja) 2002-02-15 2002-02-15 制御装置
DE10306278A DE10306278B4 (de) 2002-02-15 2003-02-14 Regelvorrichtung und Regelverfahren
US10/366,670 US7124013B2 (en) 2002-02-15 2003-02-14 Control device, control method, control unit, and engine control unit
US11/503,950 US7647157B2 (en) 2002-02-15 2006-08-15 Control device, control method, control unit, and engine control unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002039094A JP3973922B2 (ja) 2002-02-15 2002-02-15 制御装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2006306602A Division JP4277959B2 (ja) 2006-11-13 2006-11-13 制御装置

Publications (2)

Publication Number Publication Date
JP2003241803A true JP2003241803A (ja) 2003-08-29
JP3973922B2 JP3973922B2 (ja) 2007-09-12

Family

ID=27678220

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002039094A Expired - Fee Related JP3973922B2 (ja) 2002-02-15 2002-02-15 制御装置

Country Status (3)

Country Link
US (2) US7124013B2 (ja)
JP (1) JP3973922B2 (ja)
DE (1) DE10306278B4 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005038221A1 (ja) * 2003-10-20 2005-04-28 Honda Motor Co., Ltd. 内燃機関の吸入空気量制御装置
WO2005091088A1 (ja) * 2004-03-23 2005-09-29 Honda Motor Co., Ltd. 変調アルゴリズムを用いたプラントの制御装置および制御方法
WO2006059789A1 (en) * 2004-12-02 2006-06-08 Honda Motor Co., Ltd. Air/fuel ratio control apparatus of an internal combustion engine
US7190503B2 (en) 2003-10-28 2007-03-13 Seiko Epson Corporation Electro-optical device, method of manufacturing the same, and electronic apparatus
JP2008257741A (ja) * 2008-05-28 2008-10-23 Honda Motor Co Ltd Δς変調アルゴリズムを用いてプラントを制御する制御装置
US7654252B2 (en) 2006-10-16 2010-02-02 Toyota Jidosha Kabushiki Kaisha Air-fuel ratio control system and method for internal combustion engine
JP4763713B2 (ja) * 2004-11-18 2011-08-31 ユニバーシティ オブ レスター 予測診断のためのスライディングモード方法

Families Citing this family (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3973922B2 (ja) * 2002-02-15 2007-09-12 本田技研工業株式会社 制御装置
DE102004046874A1 (de) * 2004-09-28 2006-04-13 Robert Bosch Gmbh Verfahren zum Betreiben eines Verwaltungssystems von Funktionsmodulen
US7389773B2 (en) 2005-08-18 2008-06-24 Honeywell International Inc. Emissions sensors for fuel control in engines
JP4639166B2 (ja) * 2006-05-18 2011-02-23 本田技研工業株式会社 制御装置
JP4329799B2 (ja) * 2006-09-20 2009-09-09 トヨタ自動車株式会社 内燃機関の空燃比制御装置
US8060290B2 (en) 2008-07-17 2011-11-15 Honeywell International Inc. Configurable automotive controller
US7930147B2 (en) * 2009-02-03 2011-04-19 Deere & Company Sensor response time accelerator
US8140248B2 (en) * 2009-04-07 2012-03-20 General Electric Company System and method for obtaining an optimal estimate of NOx emissions
US8620461B2 (en) 2009-09-24 2013-12-31 Honeywell International, Inc. Method and system for updating tuning parameters of a controller
US8224557B2 (en) * 2009-09-30 2012-07-17 GM Global Technology Operations LLC Control systems and methods using geometry based exhaust mixing model
US8583349B2 (en) * 2009-11-05 2013-11-12 GM Global Technology Operations LLC Systems and methods for diagnosing oxygen sensors and catalytic converters of exhaust systems
US8504175B2 (en) 2010-06-02 2013-08-06 Honeywell International Inc. Using model predictive control to optimize variable trajectories and system control
US20120283849A1 (en) * 2011-05-06 2012-11-08 Kureemun Ridwan Sensor system having time lag compensation
CN102354112B (zh) * 2011-07-22 2013-04-03 杭州和迈科技有限公司 一种船平航状态下调距桨螺距和主机转速的控制算法
CN102385314B (zh) * 2011-07-22 2013-02-20 杭州和迈科技有限公司 一种船启动时调距桨螺距和主机转速匹配的控制算法
US9677493B2 (en) 2011-09-19 2017-06-13 Honeywell Spol, S.R.O. Coordinated engine and emissions control system
US9650934B2 (en) 2011-11-04 2017-05-16 Honeywell spol.s.r.o. Engine and aftertreatment optimization system
US20130111905A1 (en) 2011-11-04 2013-05-09 Honeywell Spol. S.R.O. Integrated optimization and control of an engine and aftertreatment system
JP5883140B2 (ja) * 2012-07-17 2016-03-09 本田技研工業株式会社 内燃機関の制御装置
GB2516035B (en) * 2013-07-08 2017-03-29 Jaguar Land Rover Ltd Adaptive powertrain control for optimized performance
JP6269371B2 (ja) * 2014-07-28 2018-01-31 トヨタ自動車株式会社 内燃機関
EP3051367B1 (en) 2015-01-28 2020-11-25 Honeywell spol s.r.o. An approach and system for handling constraints for measured disturbances with uncertain preview
EP3056706A1 (en) 2015-02-16 2016-08-17 Honeywell International Inc. An approach for aftertreatment system modeling and model identification
EP3091212A1 (en) 2015-05-06 2016-11-09 Honeywell International Inc. An identification approach for internal combustion engine mean value models
EP3125052B1 (en) 2015-07-31 2020-09-02 Garrett Transportation I Inc. Quadratic program solver for mpc using variable ordering
US10272779B2 (en) 2015-08-05 2019-04-30 Garrett Transportation I Inc. System and approach for dynamic vehicle speed optimization
US10415492B2 (en) 2016-01-29 2019-09-17 Garrett Transportation I Inc. Engine system with inferential sensor
FR3049084B1 (fr) * 2016-03-15 2022-11-11 Fraunhofer Ges Forschung Dispositif de codage pour le traitement d'un signal d'entree et dispositif de decodage pour le traitement d'un signal code
US10124750B2 (en) 2016-04-26 2018-11-13 Honeywell International Inc. Vehicle security module system
US10036338B2 (en) 2016-04-26 2018-07-31 Honeywell International Inc. Condition-based powertrain control system
CN105892295B (zh) * 2016-04-28 2018-09-11 上海交通大学 一种提高拟滑模域特性的离散滑模控制方法
US20180058350A1 (en) * 2016-08-31 2018-03-01 GM Global Technology Operations LLC Method and apparatus for controlling operation of an internal combustion engine
US11199120B2 (en) 2016-11-29 2021-12-14 Garrett Transportation I, Inc. Inferential flow sensor
US11057213B2 (en) 2017-10-13 2021-07-06 Garrett Transportation I, Inc. Authentication system for electronic control unit on a bus
US11829935B2 (en) * 2018-04-03 2023-11-28 Florida A&M University Application of a multi-objective optimization model for automatic vessel scheduling in liner shipping
DE102018207703A1 (de) * 2018-05-17 2019-11-21 Robert Bosch Gmbh Verfahren und Vorrichtung zum Betreiben einer Abgasnachbehandlungseinrichtung eines Motorsystems mit einem Verbrennungsmotor
CN115598983B (zh) * 2022-10-28 2023-06-13 湖南大学无锡智能控制研究院 一种考虑时变时延的无人车横纵向协同控制方法及装置

Family Cites Families (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01216047A (ja) * 1988-02-24 1989-08-30 Hitachi Ltd エンジンの空燃比制御方法および装置
US5070846A (en) * 1990-11-26 1991-12-10 General Motors Corporation Method for estimating and correcting bias errors in a software air meter
US5266907A (en) * 1991-06-25 1993-11-30 Timeback Fll Continuously tuneable frequency steerable frequency synthesizer having frequency lock for precision synthesis
US5363647A (en) * 1992-10-13 1994-11-15 Mitsubishi Denki Kabushiki Kaisha Dual-sensor type air fuel ratio control system for internal combustion engine and catalytic converter diagnosis apparatus for the same
JP3206357B2 (ja) * 1994-04-19 2001-09-10 トヨタ自動車株式会社 内燃機関の燃料噴射量制御装置
EP0682304B1 (fr) * 1994-05-11 1999-12-29 CSEM Centre Suisse d'Electronique et de Microtechnique S.A. - Recherche et Développement Microsystème à faible consommation d'énergie
JP3315814B2 (ja) * 1994-06-09 2002-08-19 トヨタ自動車株式会社 車両用制御装置,クラッチのスリップ制御装置およびアイドル回転数制御装置
US5467185A (en) * 1994-07-15 1995-11-14 General Electric Company Emissions control for internal combustion engine
JP3156604B2 (ja) * 1996-02-28 2001-04-16 トヨタ自動車株式会社 内燃機関の空燃比制御装置
US5992383A (en) * 1996-05-28 1999-11-30 U.S. Philips Corporation Control unit having a disturbance predictor, a system controlled by such a control unit, an electrical actuator controlled by such a control unit, and throttle device provided with such an actuator
JP3581762B2 (ja) * 1996-06-20 2004-10-27 トヨタ自動車株式会社 内燃機関の空燃比制御装置
JP3521790B2 (ja) * 1998-03-25 2004-04-19 株式会社デンソー 内燃機関の制御装置
JP3484074B2 (ja) * 1998-05-13 2004-01-06 本田技研工業株式会社 プラントの制御装置
JP2000043609A (ja) 1998-07-30 2000-02-15 Mazda Motor Corp 車両用定速走行装置
JP3572961B2 (ja) * 1998-10-16 2004-10-06 日産自動車株式会社 エンジンの排気浄化装置
US6594613B1 (en) * 1998-12-10 2003-07-15 Rosemount Inc. Adjustable bandwidth filter for process variable transmitter
JP3484088B2 (ja) * 1998-12-17 2004-01-06 本田技研工業株式会社 プラントの制御装置
US6535153B1 (en) * 1999-02-04 2003-03-18 Med-El Electromedizinische Gerate Ges.M.B.H. Adaptive sigma-delta modulation with one-bit quantization
JP3868693B2 (ja) * 1999-03-03 2007-01-17 日産自動車株式会社 内燃機関の空燃比制御装置
JP2000332611A (ja) 1999-05-17 2000-11-30 Matsushita Electric Ind Co Ltd 音声信号処理装置
JP3655145B2 (ja) 1999-10-08 2005-06-02 本田技研工業株式会社 多気筒内燃機関の空燃比制御装置
JP2001154704A (ja) * 1999-11-26 2001-06-08 Prime Motion:Kk フィードバック制御システム
JP4354068B2 (ja) * 2000-02-02 2009-10-28 本田技研工業株式会社 内燃機関の排ガスの空燃比制御装置
JP3680217B2 (ja) * 2000-06-26 2005-08-10 トヨタ自動車株式会社 内燃機関の空燃比制御装置
JP3467455B2 (ja) * 2000-07-17 2003-11-17 本田技研工業株式会社 内燃機関の気筒別空燃比推定装置
JP3922980B2 (ja) * 2001-07-25 2007-05-30 本田技研工業株式会社 制御装置
JP3904923B2 (ja) * 2001-12-28 2007-04-11 本田技研工業株式会社 制御装置
US6633793B2 (en) * 2001-08-13 2003-10-14 Promos Technologies Method to reduce lot-to-lot variation of array threshold voltage in a DRAM device
JP3973922B2 (ja) * 2002-02-15 2007-09-12 本田技研工業株式会社 制御装置
JP4028334B2 (ja) * 2002-09-12 2007-12-26 本田技研工業株式会社 制御装置

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005038221A1 (ja) * 2003-10-20 2005-04-28 Honda Motor Co., Ltd. 内燃機関の吸入空気量制御装置
CN100434670C (zh) * 2003-10-20 2008-11-19 本田技研工业株式会社 内燃机的吸入空气量控制装置
US7568454B2 (en) 2003-10-20 2009-08-04 Honda Motor Co., Ltd. Intake air amount control system for internal combustion engine
US7190503B2 (en) 2003-10-28 2007-03-13 Seiko Epson Corporation Electro-optical device, method of manufacturing the same, and electronic apparatus
WO2005091088A1 (ja) * 2004-03-23 2005-09-29 Honda Motor Co., Ltd. 変調アルゴリズムを用いたプラントの制御装置および制御方法
US7502681B2 (en) 2004-03-23 2009-03-10 Honda Motor Co., Ltd. Plant control device and control method using modulation algorithm
JP4763713B2 (ja) * 2004-11-18 2011-08-31 ユニバーシティ オブ レスター 予測診断のためのスライディングモード方法
WO2006059789A1 (en) * 2004-12-02 2006-06-08 Honda Motor Co., Ltd. Air/fuel ratio control apparatus of an internal combustion engine
US7680580B2 (en) 2004-12-02 2010-03-16 Honda Motor Co., Ltd. Air/fuel ratio control apparatus of an internal combustion engine
US7654252B2 (en) 2006-10-16 2010-02-02 Toyota Jidosha Kabushiki Kaisha Air-fuel ratio control system and method for internal combustion engine
JP2008257741A (ja) * 2008-05-28 2008-10-23 Honda Motor Co Ltd Δς変調アルゴリズムを用いてプラントを制御する制御装置
JP4630354B2 (ja) * 2008-05-28 2011-02-09 本田技研工業株式会社 Δς変調アルゴリズムを用いてプラントを制御する制御装置

Also Published As

Publication number Publication date
US7647157B2 (en) 2010-01-12
DE10306278B4 (de) 2013-09-05
JP3973922B2 (ja) 2007-09-12
US7124013B2 (en) 2006-10-17
DE10306278A1 (de) 2003-10-30
US20030154953A1 (en) 2003-08-21
US20060282211A1 (en) 2006-12-14

Similar Documents

Publication Publication Date Title
JP2003241803A (ja) 制御装置
JP3922980B2 (ja) 制御装置
JP3904923B2 (ja) 制御装置
JP3484088B2 (ja) プラントの制御装置
JP3926703B2 (ja) 制御装置
JP3621839B2 (ja) プラントの制御装置
JP3880861B2 (ja) 内燃機関の空燃比制御装置
JP2003003892A (ja) 内燃機関の空燃比制御装置
JP4028334B2 (ja) 制御装置
JP2000234550A (ja) 内燃機関の空燃比制御装置
JP4439508B2 (ja) 制御装置
JP4277958B2 (ja) 内燃機関の空燃比制御装置
JP4277959B2 (ja) 制御装置
JP4240325B2 (ja) 制御装置
EP1279820B1 (en) Control apparatus, control method, and engine control unit
JP2003074399A (ja) 内燃機関の空燃比制御装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041130

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060825

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060912

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20061113

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20070130

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070329

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20070412

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20070515

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070613

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100622

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110622

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110622

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130622

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130622

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140622

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees