CN101469640A - 控制装置 - Google Patents

Abstract

本发明提供如下控制装置：在利用表示周期性的变动行为的控制输入对控制量进行控制以使其收敛于目标控制量的情况下，可迅速消除由于控制对象的动态特性的骤变引起的控制输入的变动频率的下降，可提高控制精度。控制装置(1)计算空燃比校正值DKCMD，以使氧浓度传感器(15)的输出VO2收敛于目标输出VO2_TRGT，通过使用应用Δ∑调制算法的算法调制该空燃比校正值DKCMD，来计算调制值DKCMD_DSM。根据排气容积SV_hat计算基准空燃比KCMDBS，使用计算成使空燃比校正值DKCMD为值0的修正系数θ_i来计算模型修正系数KTRQFF，使用KCMDADP＝KCMDBS×KTRQFF的算式计算自适应基准空燃比KCMDADP，使用KCMD＝KCMDADP+DKCMD_DSM的算式计算目标空燃比KCMD。

Description

控制装置

技术领域

本发明涉及控制装置，该控制装置使用规定的调制算法来对使用反馈控制方法以使控制量收敛于目标控制量而计算出的值进行调制，并使用该调制值计算针对控制对象的控制输入。

背景技术

作为控制内燃机的混合气的空燃比的控制装置，本申请人已提出了在专利文献1中记载的控制装置。该控制装置具有LAF传感器、氧浓度传感器、状态预测器、车载同定器、滑模控制器以及目标空燃比计算部等。这些LAF传感器和氧浓度传感器均是用于检测表示内燃机的排气通路内的排气中的氧浓度的值即空燃比的传感器，并设置在排气通路的集合部的下游侧。并且，LAF传感器配置在催化装置的上游侧，氧浓度传感器配置在催化装置的下游侧。

在该控制装置中，作为控制对象模型，使用将LAF传感器的检测空燃比KACT与学习值FLAFBASE之间的偏差(以下称为“空燃比偏差”)DKACT作为输入、并将氧浓度传感器的输出VOUT与规定的目标值VOUT_TARGET之间的偏差(以下称为“输出偏差”)DVO2作为输出的离散时间系模型，按以下所述，计算作为控制输入的目标空燃比KCMD(目标当量比)。

即，在状态预测器中，使用基于上述控制对象模型的规定的预测算法来计算输出偏差DVO2的预测值，在车载同定器中，使用逐次型最小平方法来同定控制对象模型的模型参数。并且，在滑模控制器中，根据输出偏差的预测值和模型参数的同定值，使用滑模控制算法计算操作量Usl，以使输出偏差DVO2收敛于值0。

而且，在目标空燃比计算部中，通过对一定值flafbase加上学习分量flafadp来计算学习值FLAFBASE。当作为操作量Usl的一分量Uadp在规定范围内、而且

成立时，该学习分量flafadp被保持为其上次值，当Uadp不在规定范围内时，通过对学习分量flafadp的上次值加上或减去规定值Δflaf来计算该学习分量flafadp。该学习值FLAFBASE的运算周期即学习速度被设定为较慢的速度，以避免干扰滑模控制器的滑模控制。

然后，通过使按以上计算出的学习值FLAFBASE加上操作量Usl，来计算目标空燃比KCMD。结果，燃料喷射量被计算成使检测空燃比KACT收敛于上述目标空燃比KCMD，从而空燃比被控制成使氧浓度传感器的输出VOUT收敛于规定的目标值VOUT_TARGET。该规定的目标值VOUT_TARGET被设定为这样的值，即：当氧浓度传感器的输出VOUT为该值时，可在催化装置中获得良好的排气净化率，结果，通过以上的控制，可在催化装置中获得良好的排气净化率。

并且，作为控制具有非线性特性的控制对象的控制装置，本申请人已提出了在专利文献2中记载的控制装置。在该专利文献2的图9所示的控制装置中，在控制器51中，计算参考输入，以使排气传感器的输出Vex收敛于其目标值Vex_cmd，在Δ∑调制器52中，通过使用Δ∑调制算法调制该参考输入，来计算作为控制输入的燃料参数Ufuel。由此，可在对发动机和催化装置的响应延迟或变动进行补偿的同时，使排气传感器的输出Vex精度良好地收敛于其目标值Vex_cmd，在催化装置中，可确保良好的排气特性。

【专利文献1】日本特开2000—234550号公报

【专利文献2】日本特开2005—275489号公报

根据上述现有的专利文献1的控制装置，由于根据催化装置的下游侧的氧浓度传感器的输出VOUT来控制空燃比，因而有发生以下所述的问题的可能性。即，一般，在多气缸内燃机中空燃比在气缸间发生变动的情况下，从各气缸排出的排气尽管在排气通路的集合部的下游侧混合，然而其混合程度有限制，当排气流入催化装置时，有时在催化装置内排气的空燃比处于不均匀的状态。因此，例如在氧浓度传感器的输出VOUT收敛于规定的目标值VOUT_TARGET的情况下，即在成立的情况下，有时发生催化装置的一半为浓气氛、而剩余一半为稀气氛的状态，在这种状态长时间持续的情况下，催化装置整体的排气净化率下降，排气特性恶化。

为了解决以上的专利文献1的控制装置的问题，考虑了将专利文献2的控制方法应用于专利文献1的控制装置，并按以下来构成。即，构成为，通过使用Δ∑调制算法来调制所述的操作量Usl，计算调制操作量，通过对该调制操作量加上学习值FLAFBASE，计算目标空燃比KCMD。在这样构成的情况下，目标空燃比KCMD是这样来计算，即：当

成立时，以学习值FLAFBASE为中心以高频率来重复周期性的变动行为。结果，当学习值FLAFBASE达到使

成立的目标空燃比KCMD的最佳值时，可避免由于各气缸的空燃比表示高频率的周期性的变动行为而使催化装置内的气氛长时间维持在不均匀的状态，可确保良好的排气特性。

然而，当学习值FLAFBASE大幅偏离上述目标空燃比KCMD的最佳值时，由于氧浓度传感器的输出VOUT大幅偏离规定的目标值VOUT_TARGET而使目标空燃比KCMD的变动频率下降。本申请人已确认，在这样目标空燃比KCMD的变动频率下降的情况下，尽管当催化装置未恶化时能将催化装置的排气净化率保持在良好的状态，然而当催化装置恶化时，排气净化率显著下降(例如，专利3880861号公报的段落“0202”～“0205”以及图38)。

近年，使用了乙醇混合燃料作为内燃机的燃料，使用了钙钛矿等非金属催化剂或低贵金属保持催化剂等作为催化装置中的催化剂，在这种情况下，伴随内燃机的运转状态的变化(例如排气容积的变化)，使DVO2≒0成立的目标空燃比KCMD的最佳值容易大幅骤变。与此相对，由于学习值FLAFBASE如上所述具有学习速度慢的特性，因而如上所述，在伴随内燃机的运转状态的变化、目标空燃比KCMD的最佳值大幅骤变的情况下，学习值FLAFBASE容易大幅偏离其最佳值，目标空燃比KCMD的变动频率容易显著下降，使得上述问题变得显著。

发明内容

本发明是为了解决上述课题而作成的，本发明的目的是提供一种控制装置，该控制装置在使用表示周期性的变动行为的控制输入对控制量进行控制以使其收敛于目标控制量的情况下，可迅速地消除由于控制对象的动态特性的骤变引起的控制输入的变动频率的下降，由此可提高控制精度。

为了达到上述目的，根据发明1的发明是一种控制装置1、1A，其使用控制输入(目标空燃比KCMD、发动机转矩TRQENG)对控制对象中的控制量(氧浓度传感器15的输出VO2、驱动轮速Ws_act)进行控制，其特征在于，该控制装置具有：控制量检测单元(ECU 2、氧浓度传感器15、左右的后轮速传感器22、23)，其检测控制量；目标控制量设定单元(ECU2、目标车轮速计算部101)，其设定成为控制量的目标的目标控制量(目标输出VO2_TRGT、目标车轮速Ws_cmd)；第1输入值计算单元(ECU 2、SM(滑模)控制器43、车轮速反馈控制器102)，其使用规定的反馈控制算法[式(19)～(23)、(48)～(58)]计算用于对控制量进行反馈控制以使其收敛于目标控制量的第1输入值(空燃比校正值DKCMD、转矩反馈值TRQFB)；调制输入值计算单元(ECU 2、调制器44、103)，其通过使用规定的调制算法对第1输入值进行调制来计算调制输入值(调制值DKCMD_DSM、TRQFB_DSM)，以使得第1输入值的绝对值越小则调制输入值的频率越高；动作状态参数检测单元(ECU 2、空气流量传感器10、进气压传感器11、曲轴转角传感器13、进气温度传感器17)，其检测控制量以外的、表示控制对象的动作状态的动作状态参数(排气容积SV_hat、发动机转速NE)；误差参数计算单元(ECU 2、SM控制器43、模型修正器120)，其计算误差参数(空燃比校正值DKCMD、转矩误差Etrq)，该误差参数表示控制输入相对于被估计为控制量成为目标控制量时的控制输入的估计值(最佳空燃比、估计转矩TRQINST)的误差；修正值计算单元(ECU 2、模型修正器70、120)，其使用规定的控制算法分别计算与将动作状态参数可变化的区域进行划分而成的多个区域对应的多个修正值(修正系数θ_i、基准局部校正值DKTRQFF_BS_ij)，以使误差参数成为规定的目标值(值0)；模型修正单元(ECU 2、自适应基准空燃比计算部50、牵引控制器100)，其使用多个修正值来修正相关关系模型(图6、19)，该相关关系模型表示用于对控制量进行前馈控制的第2输入值与动作状态参数之间的关系；第2输入值计算单元(ECU 2、自适应基准空燃比计算部50、牵引控制器100)，其使用修正后的相关关系模型和动作状态参数来计算第2输入值(自适应基准空燃比KCMDADP、自适应限制转矩TRQFF_LMT)；以及控制输入计算单元(ECU 2、空燃比控制器40、牵引控制器100)，其根据第2输入值和调制输入值的和，将控制输入(目标空燃比KCMD、发动机转矩TRQENG)计算成表现出周期性的变动行为。

根据该控制装置，使用规定的反馈控制算法来计算用于对控制量进行反馈控制以使其收敛于目标值的第1输入值，通过使用规定的调制算法对第1输入值进行调制来计算调制输入值，以使第1输入值的绝对值越小则调制输入值的频率越高，使用修正后的相关关系模型和动作状态参数来计算用于对控制量进行前馈控制的第2输入值，根据第2输入值和调制输入值的和来计算控制输入。由于按以上所述来计算控制输入，因而在由于控制对象的动态特性的骤变而使控制输入相对于控制输入的估计值的误差骤增时，如果使用第1输入值来补偿该误差，由于第1输入值是使用反馈控制算法来计算的，因而需要时间，控制精度下降。

与此相对，第2输入值是通过分别计算与将动作状态参数可变化的区域进行划分而成的多个区域对应的多个修正值以使误差参数成为规定的目标值，并使用利用这些多个修正值进行了修正的相关关系模型和动作状态参数来计算的。即，第2输入值是使用针对多个区域中的每个区域进行了修正的相关关系模型来计算的，以使误差参数成为规定的目标值，因而即使当由于控制对象的动态特性的骤变而使控制输入相对于控制输入的估计值的误差骤增时，也能通过利用修正后的相关关系模型计算出的第2输入值来迅速且精度良好地补偿这种误差。

特别是，作为相关关系模型，使用在前馈控制方法中一般性的、表示动作状态参数与第2输入值之间的相关关系的N(N是自然数)维映射图、或表示两者的相关关系的计算式等，从而与根据利用反馈控制算法计算出的第1输入值来补偿误差参数即相对于控制输入的估计值的误差的情况相比，可迅速地补偿该误差参数。由此，即使在由于控制对象的动态特性变化等而使控制输入相对于控制输入的估计值的误差暂时骤增的条件下，也能将控制量迅速地控制在目标控制量附近，可减小第1输入值的绝对值，从而可迅速地提高调制输入值的频率即控制输入的变动频率。这样，可迅速地消除由于控制对象的动态特性的骤变引起的、控制输入的变动频率的下降，由此可提高控制精度(另外，本说明书中的“相关关系模型”不限于响应曲面模型和数学性模型，包含N(N是自然数)维映射图和规定的计算算法等表示动作状态参数与第2输入值之间的相关关系的全部模型。并且，本说明书中的“动作状态参数的检测”不限于使用传感器直接检测动作状态参数，还包含通过计算或估计得到动作状态参数)。

发明2的发明，其特征在于，在发明1所述的控制装置1、1A中，模型修正单元通过把规定的多个函数(权重函数W_i、WNE_ij)的值乘以误差参数和规定的目标值中的一方和另一方之间的偏差(空燃比校正值DKCMD、转矩误差Etrq)，分别计算多个相乘值(分配误差Ew_i、修正转矩误差Wetrq_ij)，并根据多个相乘值分别计算多个修正值(修正系数θ_i、基准局部校正值DKTRQFF_BS_ij)，多个区域中的相邻区域相互重叠，规定的多个函数分别对应于多个区域，规定的多个函数仅在对应的区域内被设定为值0以外的值且在对应的区域以外被设定为值0，并且在相互重叠的区域内，与重叠区域对应的函数的值的总和的绝对值(值1)被设定成等于函数中的最大值的绝对值(值1)。

根据该控制装置，规定的多个函数分别对应于多个区域，仅在对应的区域内被设定为值0以外的值且在对应的区域以外被设定为值0，并且在相互重叠的区域内，与该重叠的区域对应的函数的值的总和的绝对值被设定成等于函数中的最大值的绝对值。通过把这种规定的多个函数的值乘以误差参数和规定的目标值中的一方和另一方之间的偏差，分别计算多个相乘值，并根据多个相乘值分别计算多个修正值，因而可将误差参数和规定的目标值中的一方和另一方之间的偏差通过规定的多个函数的值分配给多个修正值，由此，可在通过多个修正值中的各个修正值来反映多个区域中的各个区域的误差程度的同时，恰当地修正相关关系模型。而且，由于与相互重叠的区域对应的函数的值的总和的绝对值被设定成等于函数中的最大值的绝对值，因而使用这种多个函数的值计算出的多个修正值表示相互连续的值，从而即使当动作状态参数骤变时，也能将第2输入值计算成平滑且无级差地变化。由此，即使当由于控制对象的动态特性的骤变而使控制输入相对于控制输入的估计值的误差暂时骤增时，也能避免第2输入值的不恰当的骤变或阶梯状的骤变，可提高控制精度和控制稳定性。

发明3的发明，其特征在于，在发明1或2所述的控制装置中，误差参数计算单元具有设定单元(ECU 2、空燃比控制器40)，该设定单元将第1输入值(空燃比校正值DKCMD)设定为误差参数，并将规定的目标值设定为值0。

根据该控制装置，由于第2输入值是使用针对多个区域中的每个区域进行了修正的相关关系模型来计算的，以使第1输入值成为值0，因而即使当由于控制对象的动态特性的骤变而使控制输入大幅偏离估计值、并且控制量大幅偏离目标控制量，从而使得第1输入值骤增时，也能通过利用修正后的相关关系模型计算出的第2输入值将骤增的第1输入值迅速控制为值0。结果，能够将控制量迅速控制为目标控制量，从而可迅速提高调制输入值的频率即控制输入的变动频率，由此可提高控制精度。

发明4的发明，其特征在于，在发明3所述的控制装置1中，控制量是排气浓度传感器(氧浓度传感器15)的输出VO2，该排气浓度传感器检测在内燃机3的排气通路7中的催化装置(第1催化装置8)的下游侧的排气中的规定成分的浓度(氧浓度)，目标控制量是当估计为催化装置的排气净化率成为规定状态时的目标输出VO2_TRGT，动作状态参数是表示内燃机的运转状态的运转状态参数(排气容积SV_hat)，控制输入计算单元计算在内燃机3的混合气的空燃比控制中成为空燃比的目标的目标空燃比KCMD作为控制输入。

根据该控制装置，由于第2输入值是使用针对多个区域中的每个区域进行了修正的相关关系模型和运转状态参数来计算的，以使第1输入值成为值0，并且根据第2输入值和调制输入值的和来计算目标空燃比，因而即使当由于内燃机的运转状态的骤变而使计算出的目标空燃比大幅偏离估计值、并且排气浓度传感器的输出大幅偏离目标输出，从而使得第1输入值骤增时，也能通过利用修正后的相关关系模型计算出的第2输入值将骤增的第1输入值迅速控制为值0。由此，可将排气浓度传感器的输出迅速控制为目标输出，从而可迅速提高调制输入值的频率即目标空燃比的变动频率。结果，可迅速消除由于内燃机的运转状态的骤变引起的、目标空燃比的变动频率的下降，并且通过将该目标输出设定为使催化装置的排气净化率成为良好状态的值，可提高催化装置的排气净化率。

发明5的发明，其特征在于，在发明1或2所述的控制装置1A中，误差参数计算单元具有：估计值计算单元(ECU 2、估计转矩计算部106)，其计算控制输入的估计值(估计转矩TRQINST)；以及误差参数设定单元(ECU 2、转矩误差计算部122)，其将估计值与第2输入值之间的偏差(转矩误差Etrq)设定为误差参数。

根据该控制装置，由于第2输入值是使用针对多个区域中的每个区域进行了修正的相关关系模型来计算的，以使控制输入的估计值与第2输入值之间的偏差成为值0，因而即使当由于控制对象的动态特性的骤变而使控制输入大幅偏离估计值、并且控制量大幅偏离目标控制量，从而使得第1输入值骤增时，也能通过利用修正后的相关关系模型计算出的第2输入值将控制量迅速控制为目标控制量，可迅速减小第1输入值。由此，可迅速提高调制输入值的频率即控制输入的变动频率，可提高控制精度。

发明6的发明，其特征在于，在发明5所述的控制装置1A中，控制量是以内燃机3为动力源的车辆的驱动轮的速度(驱动轮速Ws_act)，目标控制量是用于使车辆的驱动轮相对于非驱动轮的滑移量成为规定量(滑移偏置值OptSlip)的目标速度(目标车轮速Ws_cmd)，动作状态参数是表示内燃机的运转状态的运转状态参数(发动机转速NE)，控制输入计算单元计算在内燃机3的输出控制中成为输出的目标的目标输出(限制转矩TRQTC)作为控制输入。

如该控制装置这样，在将内燃机的目标输出作为控制输入来控制车辆的驱动轮的速度的情况下，在驱动轮发生滑移的前后，轮胎与路面之间的摩擦特性大幅地非线性变化(以下将这种特性称为“抓地行为的非线性特性”)，因而即使当驱动轮的速度收敛于目标速度时，也很有可能驱动轮的滑移量发生波动、或者相对于规定值产生恒定偏差，为了避免这一点，有必要将目标输出计算成表现出高频率的变动行为。与此相对，根据该控制装置，由于第2输入值是使用针对多个区域中的每一个区域进行了修正的相关关系模型和运转状态参数来计算的，以使第1输入值成为值0，并且根据第2输入值和调制输入值的和来计算目标输出，因而即使当由于内燃机的运转状态的骤变而使计算出的目标输出大幅偏离其估计值、并且驱动轮的速度大幅偏离目标速度，从而使得第1输入值骤增时，也能通过利用修正后的相关关系模型计算出的第2输入值将骤增的第1输入值迅速控制为值0。由此，可将驱动轮的速度迅速控制为目标速度，从而可迅速提高调制输入值的频率即目标输出的变动频率。结果，可补偿上述的抓地行为的非线性特性，从而可避免驱动轮的滑移量发生波动、或者相对于规定值产生恒定偏差的情况，并且通过将该规定量设定为极小的值，可确保高水准的驱动轮速度的控制精度。即，可实现高精度的牵引控制。

发明7的发明，其特征在于，在发明6所述的控制装置1A中，控制装置还具有：要求输出计算单元(ECU 2、要求转矩计算部104)，其计算车辆的驾驶者所要求的内燃机3的输出作为要求输出(要求转矩TRQDRV)；以及选择单元(ECU 2、发动机转矩计算部110)，其将要求输出(要求转矩TRQDRV)和目标输出(限制转矩TRQTC)中的小的一方选择为输出控制中的内燃机3的输出目标(发动机转矩TRQENG)。

根据该控制装置，由于将车辆的驾驶者所要求的要求输出和目标输出中的小的一方选择为输出控制中的内燃机的输出的目标，因而在要求输出大于目标输出的情况下，将目标输出选择为目标，因而如上所述，可执行高精度的牵引控制。另一方面，在目标输出大于要求输出的情况下，将要求输出选择为目标，因而可获得驾驶者要求的内燃机的输出，可确保良好的运转性能。

发明8的发明，其特征在于，在发明7所述的控制装置1A中，控制输入值计算单元在第2输入值(自适应限制转矩TRQFF_LMT)小于要求输出(要求转矩TRQDRV)时，将目标输出(限制转矩TRQTC)计算为第2输入值和调制输入值的和(式(89)、(91))，并在第2输入值大于等于要求输出时，将目标输出计算为要求输出和调制输入值的和(式(90)、(91))。

根据该控制装置，当第2输入值小于要求输出时，将目标输出计算为第2输入值和调制输入值的和，并当第2输入值大于等于要求输出时，将目标输出计算为要求输出和调制输入值的和。即，目标输出被计算为第2输入值和要求输出中的小的一方与调制输入值的和，并且目标输出和要求输出中的小的一方被选择为内燃机的输出的目标。因此，在要求输出大于第2输入值而很有可能发生驱动轮的滑移的情况下，把目标输出计算为第2输入值和调制输入值的和，并且把这样计算出的目标输出选择为内燃机的输出的目标，因而可将驱动轮的速度迅速控制为目标车轮速，可执行高精度的牵引控制。

另一方面，在要求输出小于目标输出的情况下，把要求输出选择为内燃机的输出的目标，因而可避免驱动轮的滑移。在该情况下，在使用反馈控制算法来计算第1输入值的关系上，有时第1输入值会成为负值。在该状态下，由于驾驶者的急剧的油门操作等而使要求输出骤增并超过目标输出的情况下，尽管把目标输出选择为内燃机的输出的目标，然而在使用反馈控制算法来计算第1输入值的关系上，其具有响应性低的特性。因此，当第1输入值为负值时，目标输出被抑制为过度小的值，内燃机产生的输出被过度地抑制，从而存在车辆的加速响应性降低、运转性能降低的可能性。

与此相对，由于把目标输出计算为第2输入值和要求输出中的小的一方与调制输入值的和，因而在要求输出骤增并大幅超过第2输入值的情况下，目标输出被计算为可避免驱动轮的滑移的第2输入值和调制输入值的和，从而可将内燃机的输出的目标计算为以第2输入值为中心重复反转的值。由此，可在避免把目标输出抑制为过度小的值的同时，将内燃机产生的输出保持在可避免驱动轮的滑移的输出的上限值附近，从而可在执行高精度的牵引控制的同时，获得良好的加速响应性。

附图说明

图1是示出本发明的第1实施方式的控制装置和应用该控制装置的内燃机的概略结构的示意图。

图2是示出在使用了恶化状态和未恶化状态的第1催化装置的情况下，分别测定出相对于检测空燃比KACT的、两个第1催化装置的HC和NOx的净化率以及O2传感器的输出VO2的结果的一例的图。

图3是示出第1实施方式的控制装置的概略结构的功能框图。

图4是示出空燃比控制器的概略结构的功能框图。

图5是示出自适应基准空燃比计算部的概略结构的功能框图。

图6是示出在基准空燃比KCMDBS的计算中使用的映射图的一例的图。

图7是表示在第1催化装置是新品和恶化品的情况下，排气容积SV_hat与使

成立的混合气的空燃比的最佳值(最佳空燃比)之间的关系的图。

图8(a)是示出在权重函数W_i的计算中使用的映射图的一例的图，图8(b)是示出在第1催化装置是新品和恶化品的情况下的基准空燃比KCMDBS的一例的图。

图9(a)是示出权重函数W_i的计算结果的一例的图，图9(b)是示出相乘值θ_i·W_i的计算结果的一例的图，图9(c)是示出累积值Kff’的计算结果的一例的图，图9(d)是示出模型修正系数Kff的计算结果的一例的图，图9(e)是示出自适应基准空燃比KCMDADP的计算结果的一例的图。

图10是示出空燃比控制处理的流程图。

图11是示出在目标空燃比的起动时用值KCMD_ST的计算中使用的映射图的一例的图。

图12是示出在目标空燃比的催化剂暖机用值KCMD_AST的计算中使用的映射图的一例的图。

图13是示出第1实施方式的控制装置的空燃比控制的模拟结果的一例的时序图。

图14是为了比较而示出当使用现有的控制方法时的空燃比控制的模拟结果的一例的时序图。

图15是示出第2实施方式的控制装置的概略结构的框图。

图16是示出牵引控制器的概略结构的框图。

图17是示出在最大转矩Trq_max和最小转矩Trq_min的计算中使用的映射图的一例的图。

图18是示出在归一化要求驱动力Ktrq_ap的计算中使用的映射图的一例的图。

图19是示出在基准前馈转矩TRQFFM的计算中使用的映射图的一例的图。

图20是示出模型修正器的概略结构的功能框图。

图21是示出在权重函数WNE_ij的计算中使用的映射图的一例的图。

图22(a)是示出权重函数WNE_ij的计算结果的一例的图，图22(b)是示出相乘值Dktrqff_ij的计算结果的一例的图，图22(c)是示出局部校正值DKTRQFF_j的计算结果的一例的图，图22(d)是示出模型修正系数KTRQFF的计算结果的一例的图，图22(e)是示出自适应限制转矩TRQFF_LMT的计算结果的一例的图。

图23是示出第2实施方式的控制装置的牵引控制的模拟结果例的时序图。

图24是为了比较而示出在保持为模型修正系数KTRQFF＝1的情况下的牵引控制的模拟结果例的时序图。

标号说明

1：控制装置；1A：控制装置；2：ECU(控制量检测单元、目标控制量设定单元、第1输入值计算单元、调制输入值计算单元、动作状态参数检测单元、误差参数计算单元、修正值计算单元、模型修正单元、第2输入值计算单元、控制输入计算单元、设定单元、估计值计算单元、误差参数设定单元、要求输出计算单元、选择单元)；3：内燃机；7：排气通路；8：第1催化装置(催化装置)；10：空气流量传感器(动作状态参数检测单元)；11：进气压传感器(动作状态参数检测单元)；13：曲轴转角传感器(动作状态参数检测单元)；15：氧浓度传感器(控制量检测单元、排气浓度传感器)；17：进气温度传感器(动作状态参数检测单元)；22：左后轮速传感器(控制量检测单元)；23：右后轮速传感器(控制量检测单元)；40：空燃比控制器(控制输入计算单元、设定单元)；43：SM控制器(第1输入值计算单元、误差参数计算单元)；44：调制器(调制输入值计算单元)；50：自适应基准空燃比计算部(模型修正单元、第2输入值计算单元)；70：模型修正器(修正值计算单元)；100：牵引控制器(模型修正单元、第2输入值计算单元、控制输入计算单元)；101：目标车轮速计算部(目标控制量设定单元)；102：车轮速反馈控制器(第1输入值计算单元)；103：调制器(调制输入值计算单元)；104：要求转矩计算部(要求输出计算单元)；106：估计转矩计算部(估计值计算单元)；110：发动机转矩计算部(选择单元)；120：模型修正器(误差参数计算单元、修正值计算单元)；122：转矩误差计算部(误差参数设定单元)；KCMD：目标空燃比(控制输入)；VO2：氧浓度传感器的输出(控制量、排气浓度传感器的输出)；VO2_TRGT：目标输出(目标控制量)；DKCMD：空燃比校正值(第1输入值、误差参数、误差参数和规定的目标值中的一方与另一方之间的偏差)；DKCMD_DSM：调制值(调制输入值)；SV_hat：排气容积(动作状态参数、运转状态参数)；θ_i：修正系数(多个修正值)；KCMDADP：自适应基准空燃比(第2输入值)；W_i：权重函数(规定的多个函数)；Ew_i：分配误差(多个相乘值)；TRQENG：发动机转矩(控制输入)；Ws_act：驱动轮速(控制量、驱动轮的速度)；Ws_cmd：目标车轮速(目标控制量、目标速度)；OptSlip：滑移偏置值(规定量)；TRQFB：转矩反馈值(第1输入值)；TRQFB_DSM：调制值(调制输入值)；NE：发动机转速(动作状态参数、运转状态参数)；TRQINST：估计转矩(控制输入的估计值)；Etrq：转矩误差(误差参数、误差参数和规定的目标值中的一方与另一方之间的偏差、估计值与第2输入值之间的偏差)；DKTRQFF_BS_ij：基准局部校正值(多个修正值)；TRQFF_LMT：自适应限制转矩(第2输入值)；WNE_ij：权重函数(规定的多个函数)；Wetrq_ij：修正转矩误差(多个相乘值)；TRQTC：限制转矩(目标输出)；TRQDRV：要求转矩(要求输出)；TRQENG：发动机转矩(内燃机的输出目标)。

具体实施方式

以下，参照附图说明本发明的第1实施方式的控制装置。本实施方式的控制装置是控制内燃机的空燃比的控制装置，图1示出该控制装置1和应用该控制装置1的内燃机(以下称为“发动机”)3的概略结构。如该图所示，该控制装置1具有ECU 2，该ECU 2如后所述，根据发动机3的运转状态控制提供给其气缸内的混合气的空燃比。

该发动机3是安装在具有有档位式自动变速器的车辆(均未作图示)上的直列4缸型汽油发动机，具有第1～第4的4个气缸#1～#4。在该发动机3的进气通路4上，从上游侧起依次设置有空气流量传感器10、节气门机构5以及进气压传感器11。另外，在本实施方式中，空气流量传感器10和进气压传感器11相当于动作状态参数检测单元。

空气流量传感器10由热线式空气流量计构成，检测在进气通路4内流动的空气的流量，将表示该流量的检测信号输出到ECU 2。ECU 2根据空气流量传感器10的检测信号计算空气流量Gth(单位：g/sec)。

并且，节气门机构5具有节气门5a和对节气门5a进行开闭驱动的TH致动器(未作图示)等。TH致动器是将齿轮机构与连接到ECU 2的电动机(均未作图示)组合而成，通过使用来自ECU 2的控制输入信号控制TH致动器，来使节气门5a的开度变化。节气门可自由转动地设置在进气通路4的中途，并根据伴随该转动的开度变化而使空气流量Gth变化。

而且，进气压传感器11由例如半导体压力传感器等构成，检测进气通路4内的绝对压力(以下称为“进气压”)PB，并将表示该绝对压力PB的检测信号输出到ECU 2。

并且，进气通路4经由进气歧管4a的4个分支部4b与4个气缸#1～#4分别连接。在各分支部4b，在各气缸的未作图示的进气端口的上游侧安装有燃料喷射阀6。在发动机3运转时，使用来自ECU 2的控制输入信号控制各燃料喷射阀6的打开时间即燃料喷射量Tout和喷射定时。

另一方面，在发动机3的主体上安装有由例如热敏电阻等构成的水温传感器12。水温传感器12检测在发动机3的气缸体内循环的冷却水的温度即发动机水温TW，并将表示该发动机水温TW的检测信号输出到ECU2。

并且，在发动机3的曲轴(未作图示)上设置有曲轴转角传感器13。该曲轴转角传感器13伴随曲轴的旋转，向ECU 2输出均作为脉冲信号的CRK信号和TDC信号。

CRK信号每隔规定的曲轴转角(例如30°)输出1个脉冲。ECU 2根据该CRK信号计算发动机3的转速(以下称为“发动机转速”)NE。并且，TDC信号是表示各气缸的活塞(未作图示)位于比进气行程的TDC位置稍微靠前的规定的曲轴转角位置处的信号，每隔规定曲轴转角输出1个脉冲。另外，在本实施方式中，曲轴转角传感器13相当于运转状态参数检测单元。

另一方面，在排气通路7的排气歧管7a的下游侧，从上游侧起依次保持间隔地设置有第1催化装置8和第2催化装置9。两个催化装置8、9均是将NOx催化剂和三元催化剂组合而成的，该NOx催化剂未作图示，是将铱催化剂(保持铱的碳化硅晶须粉末和硅的烧制体)覆盖在蜂窝状结构的基材的表面上，并在其上还覆盖钙钛矿型双氧化物(LaCoO₃粉末和硅的烧制体)而成的。

两个催化装置8、9借助NOx催化剂的氧化还原作用来净化稀燃运转时的排气中的NOx，并借助三元催化剂的氧化还原作用来净化稀燃运转以外的运转时的排气中的CO、HC和NOx。另外，两个催化装置8、9不限于将NOx催化剂和三元催化剂组合而成，只要能净化排气中的CO、HC和NOx即可。例如，两个催化装置8、9可以使用钙钛矿型催化剂等非金属催化剂和/或三元催化剂等金属催化剂构成。

在该第1和第2催化装置8、9之间安装有氧浓度传感器(以下称为O2传感器)15。该O2传感器15由氧化锆和铂电极等构成，并将基于第1催化装置8的下游侧的排气中的氧浓度的检测信号输出到ECU 2。该O2传感器15的输出(以下称为“传感器输出”)VO2在比理论空燃比浓的混合气燃烧时成为高电平的电压值(例如0.8V)，当混合气是稀时，成为低电平的电压值(例如0.2V)，并且当混合气在理论空燃比附近时，成为高电平和低电平之间的规定的目标输出VO2_TRGT(例如0.6V)(参照图2)。

另外，在本实施方式中，O2传感器15相当于控制量检测单元和排气浓度传感器，传感器输出VO2相当于控制量和排气浓度传感器的输出，目标输出VO2_TRGT相当于目标控制量。

并且，在第1催化装置8的上游侧的排气歧管7a的集合部附近安装有LAF传感器14。LAF传感器14由氧化锆和铂电极等构成，在从比理论空燃比浓的浓区域到极稀区域的大范围的空燃比的区域中，线性地检测在排气通路7内流动的排气中的氧浓度，并将表示该氧浓度的检测信号输出到ECU 2。ECU 2根据该LAF传感器14的检测信号的值，计算表示排气中的空燃比的检测空燃比KACT。另外，该检测空燃比KACT具体地说被计算为当量比。

下面，参照图2说明第1催化装置8(催化装置)的排气净化率与传感器输出VO2之间的关系。该图表示分别测定出在第1催化装置8处于因长时间使用而使净化能力下降的恶化状态和处于净化能力高的未恶化状态的情况下，当检测空燃比KACT即提供给发动机3的混合气的空燃比在理论空燃比附近变化时的2个第1催化装置8的HC和NOx的净化率以及传感器输出VO2的结果的一例。在该图中，由虚线表示的数据均是在第1催化装置8处于未恶化状态的情况下的测定结果，由实线表示的数据均是在第1催化装置8处于恶化状态的情况下的测定结果。并且表示，检测空燃比KACT越大，则混合气的空燃比在越浓侧的情况。

如该图所示，在第1催化装置8恶化的情况下，与未恶化状态的情况相比，排气净化能力下降，从而当检测空燃比KACT是更稀侧的值KACT1时，传感器输出VO2横切目标输出VO2_TRGT。另一方面，不管第1催化装置8处于恶化状态还是未恶化状态，当传感器输出VO2为目标输出VO2_TRGT时，第1催化装置8具有效率最佳地净化HC和NOx的特性。因此得知，通过控制混合气的空燃比以使传感器输出VO2成为目标输出VO2_TRGT，可利用第1催化装置8效率最佳地净化排气。因此，在后述的空燃比控制中，目标空燃比KCMD被计算成使传感器输出VO2收敛于目标输出VO2_TRGT。

而且，油门开度传感器16和进气温度传感器17等与ECU2连接。该油门开度传感器16检测车辆的未作图示的油门踏板的踩下量(以下称为“油门开度”)AP，并将表示该AP的检测信号输出到ECU 2。并且，进气温度传感器17检测进气温度TA，并将表示该进气温度TA的检测信号输出到ECU 2。另外，在本实施方式中，进气压传感器11相当于动作状态参数检测单元。

ECU 2由微型计算机构成，该微型计算机由CPU、RAM、ROM以及I/O接口(全都未作图示)等所构成，ECU2根据所述的各种传感器10～17的检测信号等，判别发动机3的运转状态，并执行各种控制。具体地说，ECU 2如后所述，根据发动机3的运转状态计算目标空燃比KCMD，并根据该目标空燃比KCMD，针对各气缸计算燃料喷射阀6的燃料喷射量Tout和喷射定时。然后，通过使用基于计算出的燃料喷射量Tout和喷射定时的控制输入信号来驱动燃料喷射阀6，来控制混合气的空燃比。

另外，在本实施方式中，ECU 2相当于：控制量检测单元、目标控制量设定单元、第1输入值计算单元、调制输入值计算单元、动作状态参数检测单元、误差参数计算单元、修正值计算单元、模型修正单元、第2输入值计算单元、控制输入计算单元、设定单元、估计值计算单元、误差参数设定单元以及选择单元。

下面，参照图3说明本实施方式的控制装置1。如该图所示，该控制装置1具有：ST控制器30、基本喷射量计算部31、总校正系数计算部32、乘法器33以及空燃比控制器40，具体地说它们由ECU 2构成。

首先，在空燃比控制器40中，如后所述，根据空气流量Gth、进气压PB、发动机转速NE、传感器输出VO2以及检测空燃比KACT来计算目标空燃比KCMD。该目标空燃比KCMD被计算为当量比。另外，在本实施方式中，空燃比控制器40相当于控制输入计算单元和设定单元，目标空燃比KCMD相当于控制输入。

并且，ST控制器30计算空燃比校正系数KAF，并具有PID控制器和自适应控制器(均未作图示)。在该PID控制器中，使用公知的PID控制方法计算空燃比校正系数KAF，以使检测空燃比KACT收敛于目标空燃比KCMD，在自适应控制器中，使用本申请人在日本特开2000—234550号公报等中已提出的自适应控制方法来计算空燃比校正系数KAF，以使检测空燃比KACT跟随目标空燃比KCMD。然后，根据发动机3的运转状态，将使用2种控制方法分别计算出的值中的一方选择为空燃比校正系数KAF。

而且，在基本喷射量计算部31中，通过根据发动机转速NE和进气压PB来检索未作图示的映射图，来计算基本喷射量Tibs。

另一方面，在总校正系数计算部32中，通过根据发动机水温TW和进气温度TA等表示运转状态的各种参数来检索未作图示的各种映射图，来计算各种校正系数，并通过把这些各种校正系数相乘，来计算总校正系数KTOTAL。

然后，在乘法器33中，根据下式(1)计算燃料喷射量Tout。

【算式1】

Tout(k)＝Tibs(k)·KAF(k)·KT0TAL(k)      ……(1)

在上式(1)中，带有符号(k)的各数据表示是按规定的控制周期ΔT(TDC信号的产生周期)所计算或采样的离散数据，符号k表示各离散数据的控制时刻。例如，符号k表示是按本次的控制定时计算出的本次值，符号k—1表示是按上次的控制定时计算出的上次值。这一点在以下的算式中也是一样。另外，在以下说明中，恰当省略了各离散数据中的符号(k)。

下面，参照图4说明所述的空燃比控制器40。如该图所示，该空燃比控制器40具有：状态预测器41、车载同定器42、SM控制器43、调制器44、排气容积计算部45、加法器46以及自适应基准空燃比计算部50。

在该空燃比控制器40中，在状态预测器41中使用后述的预测算法计算输出偏差DVO2的预测值DVPRE，在车载同定器42中使用后述的逐次型同定算法计算模型参数矢量θb。

并且，在SM控制器43中，如后所述，使用自适应滑模控制算法计算空燃比校正值DKCMD，在调制器44中，如后所述，使用应用Δ∑调制算法的调制算法计算空燃比校正值DKCMD的调制值DKCMD_DSM。

而且，在排气容积计算部45中，使用后述的控制方法计算排气容积SV_hat，在自适应基准空燃比计算部50中，使用后述的控制方法计算自适应基准空燃比KCMDADP。然后，最终在加法器46中，根据下式(2)计算目标空燃比KCMD。

【算式2】

KCMD(k)＝DKCMD_DSM(k)+KCMDADP(k-1)      ……(2)

下面，说明所述的状态预测器41。该状态预测器41使用以下所述的预测算法计算输出偏差DVO2的预测值DVPRE，首先，按下式(3)所示定义输出偏差DVO2。

【算式3】

DVO2(k)＝VO2(k)-VO2_TRGT                ……(3)

然后，按下式(4)所示定义将后述的空燃比校正值DKCMD作为输入、并将输出偏差DVO2作为输出的系统的控制对象模型。

【算式4】

DVO2(k)＝a1·DVO2(k-1)+a2·DVO2(k-2)+b1·DKCMD(k-d)    ……(4)

上式(4)的a1、a2、b1表示模型参数。并且，上式(4)的d表示从由燃料喷射阀6将目标空燃比KCMD的混合气提供给进气系统后到反映到传感器输出VO2上的浪费时间，并按下式(5)所示来定义。

【算式5】

d＝d1+d2……(5)

这里，d1表示从燃料喷射阀6到LAF传感器14的发动机系统的浪费时间，d2表示从LAF传感器14到O2传感器15的排气系统的浪费时间。

并且，预测值DVPRE是对从目标空燃比KCMD的混合气被提供给进气系统后到经过浪费时间d之后的输出偏差DVO2(k+d)进行预测得到的值，当根据上式(4)推导预测值DVPRE的计算式时，得到下式(6)。

【算式6】

$=a1\·\mathrm{DVO}2(k+d-1)+a2\·\mathrm{DVO}2(k+d-2)+b1\·\mathrm{DKCMD}\left(K\right)......\left(6\right)$

在该式(6)的情况下，需要计算与输出偏差DVO2(k)的将来值相当的DVO2(k+d—1)和DVO2(k+d—2)，实际上进行程序化是困难的。因此，使用在车载同定器42中如后所述计算的模型参数的同定值a1(k)、a2(k)、b1(k)，按下式(7)、(8)所示定义矩阵A、B，通过重复使用上式(6)的递推公式对上式(6)进行变形，则得到下式(9)。

【算式7】

$A=\left[\begin{array}{cc}a1\left(k\right)& a2\left(k\right)\\ 1& 0\end{array}\right]\·\·\·\·\·\left(7\right)$

$B=\left[\begin{array}{c}b1\left(k\right)\\ 0\end{array}\right]\·\·\·\·\·\left(8\right)$

$\mathrm{DVPRE}\left(k\right)=\mathrm{\α}1\left(k\right)\·\mathrm{DVO}2\left(k\right)+\mathrm{\α}2\left(k\right)\·\mathrm{DVO}2(k-1)+\underset{i=1}{\overset{d}{\mathrm{\Σ}}}\mathrm{\βi}\left(k\right)\·\mathrm{DKCMD}(k-i)\·\·\·\·\·\left(9\right)$

这里，α1(k)：A^d的第1行第1列要素

α2(k)：A^d的第1行第2列要素

βi(k)：A^i-1B的第1行要素

然后，当按下式(10)所示定义检测空燃比偏差DKACT时，可视为DKACT(k)＝DKCMD(k—d1)，因而当把该关系应用于所述的式(9)时，得到下式(11)。

【算式8】

DKACT(k)＝KACT(k)-KCMDADP(k-1)          ……(10)

$\mathrm{DVPRE}\left(k\right)=\mathrm{\α}1\left(k\right)\·\mathrm{DVO}2\left(k\right)+\mathrm{\α}2\left(k\right)\·\mathrm{DVO}2(k-1)$

$+\underset{i=1}{\overset{d1-1}{\mathrm{\Σ}}}\mathrm{\βi}\left(k\right)\·\mathrm{DKCMD}(k-i)+\underset{j=0}{\overset{d2}{\mathrm{\Σ}}}\mathrm{\βj}+d1\left(k\right)\·\mathrm{DKACT}(k-j)\·\·\·\·\·\left(11\right)$

这里，βj：A^j-1B的第1行要素

这里，上式(10)的KCMDADP是自适应基准空燃比，如后所述，由自适应基准空燃比计算部50计算。在本实施方式的状态预测器41中，使用按以上那样导出的式(11)来计算预测值DVPRE。

下面说明所述的车载同定器42。在该车载同定器42中，使用下式(12)～(18)所示的逐次型同定算法计算模型参数a1、a2、b1的矢量θb。

【算式9】

θb(k)＝θb(k-1)+KQ(k)·Eid(k)           ……(12)

θb(k)^T＝[a1(k)aZ(k)b1(k)]             ……(13)

Eid(k)＝DVO2(k)-DVO2_HAT(k)             ……(14)

DVO2_HAT(k)＝θb(k-1)^Tζ(k)          ……(15)

【算式10】

ζ(k)^T＝[DVO2(k-1) DVO2(k-2) DKACT(k-d2)]   ……(16)

$\mathrm{KQ}\left(k\right)=\frac{P(k-1)\mathrm{\ζ}\left(k\right)}{1+\mathrm{\ζ}{\left(k\right)}^{T}P(k-1)\mathrm{\ζ}\left(k\right)}\·\·\·\·\·\left(17\right)$

$P\left(k\right)=\frac{1}{{\mathrm{\λ}}_1}(1-\frac{{\mathrm{\λ}}_{2}P(k-1)\mathrm{\ζ}\left(k\right)\mathrm{\ζ}{\left(k\right)}^T}{{\mathrm{\λ}}_1+{\mathrm{\λ}}_2\mathrm{\ζ}{\left(k\right)}^{T}P(k-1)\mathrm{\ζ}\left(k\right)})P(k-1)\·\·\·\·\·\left(18\right)$

这里，I是单位矩阵。

在上式(12)中，矢量θb的转置矩阵按式(13)来定义，该式(12)的Eid是根据式(14)计算的跟随误差。该式(14)的DVO2_HAT是输出偏差DVO2的估计值，根据式(15)来计算。该式(15)的ζ是按式(16)来定义其转置矩阵的矢量。

并且，式(12)的KQ是根据式(17)计算的增益系数的矢量，该式(17)的P是按式(18)所示来定义的3次方阵。该式(18)的I表示3次单位矩阵，λ1、λ2分别表示被设定成使0<λ1≤1、0<λ2≤1成立的权重函数。

下面，说明所述的SM控制器43。在该SM控制器43中，使用以下所示的式(19)～(23)的自适应滑模控制算法计算空燃比校正值DKCMD。

【算式11】

σ2(k)＝DVPRE(k)+S2·DVPRE(k-1)               ......(19)

$\mathrm{Urch}2\left(k\right)=\frac{-\mathrm{Krch}2}{b1\left(k\right)}\mathrm{\&sigma;}2\left(k\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(20\right)$

$\mathrm{Uadp}2\left(k\right)=\frac{-\mathrm{Kadp}2}{b1\left(k\right)}\underset{i=0}{\overset{k}{\mathrm{\&Sigma;}}}\mathrm{\&sigma;}2\left(i\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(21\right)$

【算式12】

$\mathrm{Ueq}2\left(k\right)=\frac{-1}{b1\left(k\right)}\{[a1\left(k\right)-1+S2]\&CenterDot;\mathrm{DVPRE}\left(k\right)+[a2\left(k\right)-S2)\&CenterDot;\mathrm{DVPRE}(k-1)\}\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(22\right)$

DKCMD(k)＝Urch2(k)+Uadp2(k)+Ueq2(k)      ……(23)

上式(19)的σ2是切换函数，该式(19)的S2是被设定成使—1<S2<0的关系成立的切换函数设定参数。在该情况下，使用切换函数设定参数S2的设定值来指定预测值DVPRE向值0的收敛速度。并且，如上式(23)所示，空燃比校正值DKCMD被计算为趋近律输入Urch2、自适应律输入Uadp2以及等效控制输入Ueq2的和，该趋近律输入Urch2使用式(20)来计算。在该式(20)中，Krch2表示规定的趋近律增益。

并且，自适应律输入Uadp2使用式(21)来计算，在该式(21)中，Kadp2表示规定的自适应律增益。而且，等效控制输入Ueq2使用式(22)来计算。

如上所述，在SM控制器43中，使用式(19)～(23)的自适应滑模控制算法计算空燃比校正值DKCMD，作为用于使预测值DVPRE收敛于值0的值，即用于使传感器输出VO2收敛于目标输出VO2_TRGT的值。

另外，在本实施方式中，SM控制器43相当于第1输入值计算单元和误差参数计算单元，空燃比校正值DKCMD相当于第1输入值、误差参数、以及误差参数和规定的目标值中的一方与另一方之间的偏差，值0相当于规定的目标值。

下面，说明所述的调制器44。在该调制器44中，如以下所述计算调制值DKCMD_DSM。首先，根据下式(24)～(26)计算小振幅分量值DKCMD_L。

【算式13】

    ·DKCMD(k)<-DKCMD_L_LMT 时

DKCMD_L(k)＝-DKCMD_L_LMT           ……(24)

    ·-DKCMD_L_LMT≦DKCMD(k)≦DKCMD_L_LMT时

DKCMD_L(k)＝DKCMD(k)             ……(25)

    ·DKCMD_L_LMT<DKCMD(k)时

DKCMD_L(k)＝DKCMD_L_LMT            ……(26)

上式(24)～(26)的DKCMD_L_LMT是规定的阈值，如这些式(24)～(26)所示，小振幅分量值DKCMD_L是通过对空燃比校正值DKCMD的振幅实施以DKCMD_L_LMT为上限并以—DKCMD_L_LMT为下限的限制处理来计算的。即，小振幅分量值DKCMD_L相当于当空燃比校正值DKCMD的变动宽度小、且其绝对值处于不超过阈值DKCMD_L_LMT的范围内时的空燃比校正值DKCMD。

并且，大振幅分量值DKCMD_H根据下式(27)～(29)来计算。

【算式14】

    ·DKCMD(k)<-DKCMD_L_LMT时

DKCMD_H(k)＝DKCMD(k)+DKCMD_L_LMT           ……(27)

    ·-DKCMD_L_LMT≦DKCMD(k)≦DKCMD_L_LMT时

DKCMD_H(k)＝0              ……(28)

    ·DKCMD_L_LMT<DKCMD(k)时

DKCMD_H(k)＝DKCMD(k)-DKCMD_L_LMT        ……(29)

如以上的式(27)～(29)所示，当空燃比校正值DKCMD的绝对值不超过阈值DKCMD_L_LMT时，大振幅分量值DKCMD_H被计算为值0，当空燃比校正值DKCMD的绝对值超过阈值DKCMD_L_LMT时，大振幅分量值DKCMD_H被计算为超过部分的值。即，在由于空燃比校正值DKCMD的变动大而要求控制的快速响应性的情况下，大振幅分量值DKCMD_H被计算为用于使这种空燃比校正值DKCMD恰当地反映在目标空燃比KCMD上的值。

而且，通过使用下式(30)～(33)所示的应用Δ∑调制算法的算法来调制小振幅分量值DKCMD_L，来计算调制分量值DKCMD_L_DSM。

【算式15】

δdsm(k)＝DKCMD_L(k)-DKCMD_L_DSM(k-1)        ……(30)

σdsm(k)＝σdsm(k-1)+δdsm(k)              ……(31)

    ·σdsm(k)<0时

DKCMD_L_DSM(k)＝-DKCMD_L_AMP          ……(32)

    ·σdsm(k)≧0时

DKCMD_L_DSM(k)＝DKCMD_L_AMP          ……(33)

如上式(30)所示，偏差δdsm被计算为小振幅分量值DKCMD_L与调制分量值DKCMD_L_DSM的上次值之间的偏差。并且，在式(31)中，σdsm表示偏差δdsm的积分值。并且，式(32)、(33)的DKCMD_L_AMP是规定的振幅值。参照以上的式(30)～(33)可以明白，调制分量值DKCMD_L_DSM被计算为以值0为中心在最小值—DKCMD_L_AMP与最大值DKCMD_L_AMP之间重复反转的值。

然后，最终如下式(34)所示，通过对调制分量值DKCMD_L_DSM加上大振幅分量值DKCMD_H，来计算调制值DKCMD_DSM。

【算式16】

DKCMD_DSM(k)＝DKCMD_L_DSM(k)+DKCMD_H(k)        ……(34)

另外，在本实施方式中，调制器44相当于调制输入值计算单元，调制值DKCMD_DSM相当于调制输入值。

另一方面，在所述的排气容积计算部45中，根据以下的式(35)、(36)计算排气容积SV_hat。

【算式17】

$\mathrm{Gcy}1\left(k\right)=\frac{\mathrm{Gth}\left(k\right)}{\mathrm{NE}\left(k\right)\&CenterDot;\mathrm{SV}\_\mathrm{PRA}}-\frac{\{\mathrm{PB}\left(k\right)-\mathrm{PB}(k-1)\}\&CenterDot;\mathrm{Vb}}{R\&CenterDot;\mathrm{TA}\left(k\right)}\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(35\right)$

SV_hat(k)＝NE(k)·Gcyl(k)·SV_PRA       ……(36)

在上式(35)中，Gcyl表示吸入空气量，Vb表示进气管内体积，R表示规定的气体常数。该式(35)的SV_PRA是根据发动机3的气缸数所设定的换算系数。在该情况下，可以取代上式(35)，而使用根据进气压PB和燃料喷射量Tout来检索映射图的方法来计算吸入空气量Gcyl。而且，也可以将空气流量Gth用作排气容积SV_hat。另外，在本实施方式中，排气容积SV_hat相当于动作状态参数和运转状态参数。

下面说明所述的自适应基准空燃比计算部50。该自适应基准空燃比计算部50使用以下所述的方法计算自适应基准空燃比KCMDADP，并如图5所示，具有基准空燃比计算部60、模型修正器70以及乘法器90。另外，在本实施方式中，自适应基准空燃比计算部50相当于模型修正单元和第2输入值计算单元，自适应基准空燃比KCMDADP相当于第2输入值。

首先，在基准空燃比计算部60中，通过根据排气容积SV_hat来检索图6所示的映射图(一维映射图)，来计算基准空燃比KCMDBS。该基准空燃比KCMDBS针对该时刻的排气容积SV_hat表示使VO2≒VO2_TRGT成立的混合气的空燃比的最佳值(以下称为“最佳空燃比”)。该映射图是根据新品的第1催化装置8中的平均特性来设定的，在本实施方式中，该映射图相当于相关关系模型。

并且，在模型修正器70中，如后所述，计算模型修正系数Kff。然后，在乘法器90中，如下式(37)所示，计算自适应基准空燃比KCMDADP。

【算式18】

KCMDADP(k)＝Kff(k)·KCMDBS(k)         ……(37)

如上式(37)所示，自适应基准空燃比KCMDADP是通过使用模型修正系数Kff来修正(或校正)基准空燃比KCMDBS而计算的。换句话说，自适应基准空燃比KCMDADP相当于使用根据模型修正系数Kff对作为相关关系模型的图6的映射图进行了非线性修正后的映射图计算出的值。

下面，参照图7说明使用上述模型修正系数Kff的原因。图7中的实线所示的曲线表示在第1催化装置8是新品且未恶化的情况下的排气容积SV_hat与所述的最佳空燃比的关系，图7中的虚线所示的曲线表示在第1催化装置8恶化的情况下的两者的关系。

比较2个曲线可以明白，在第1催化装置8恶化的情况下，最佳空燃比不是相对于排气容积SV_hat的全部区域均匀地下降，其下降程度不均匀地变化。由此，在修正基准空燃比KCMDBS时，需要对应于此进行非线性的修正。因此，在模型修正器70中，为了对第1催化装置8未恶化时的基准空燃比KCMDBS进行非线性修正，使用以下所述的方法计算模型修正系数Kff。

如图5所示，模型修正器70具有：3个权重函数计算部71～73，3个SM控制器74～76，3个乘法器77～79，以及2个加法器80、81。

首先，在权重函数计算部71～73中，通过根据排气容积SV_hat来检索图8(a)所示的映射图，来分别计算3个权重函数W_i(i＝1～3)的值。在该图中，SV1～SV4表示被设定成使0<SV1<SV2<SV3<SV4成立的排气容积SV_hat的规定值。

这里，权重函数W_i的下标i表示是与以下所述的排气容积SV_hat的3个区域对应的值，该关系在后述的各种值中也是一样。具体地说，权重函数W₁被设定成对应于被规定为0≤SV_hat<SV2的第1区域，权重函数W₂被设定成对应于被规定为SV1<SV_hat<SV4的第2区域，权重函数W₃被设定成对应于被规定为SV2<SV_hat的第3区域。

并且，如该图所示，3个权重函数W_i在上述的对应的区域内分别被设定为值1以下的正值，且在除此以外的区域内分别被设定为值0，而且相邻的各2个权重函数W_m、W_m+1(m＝1或2)相互交叉，并且交叉部分的两者的和被设定为权重函数W_i的最大值1。

即，相邻的2个权重函数W₁、W₂在规定值SV1和规定值SV2之间交叉，并且在两者的交叉点处的排气容积SV_hat如图8(b)所示成为在其前后基准空燃比KCMDBS的斜率发生变化的斜率变化点。与此相同，相邻的2个权重函数W₂、W₃在规定值SV3处交叉，并且该规定值SV3也成为在其前后基准空燃比KCMDBS的斜率发生变化的斜率变化点。

这样，相邻的2个权重函数W_m、W_m+1(m＝1或2)的交叉点被设定成与基准空燃比KCMDBS的斜率变化点对应是由于以下原因。即，如图8(b)中的虚线所示，当第1催化装置8恶化时，基准空燃比KCMDBS具有上述斜率变化点而且在维持其位置的状态下变化，因而所述的模型修正系数Kff也有必要计算为针对与斜率变化点对应的各区域斜率不同。因此，由于将模型修正系数Kff的斜率设定为针对与基准空燃比KCMDBS的斜率变化点对应的各区域而不同，因而相邻的2个权重函数W_i的交叉点被设定成与基准空燃比KCMDBS的斜率变化点对应。

另外，在基准空燃比KCMDBS的斜率变化点的位置伴随第1催化装置8的恶化而移动的情况下，只要把2个权重函数W_m、W_m+1设定成在包含新品和恶化品的第1催化装置8的基准空燃比KCMDBS的斜率变化点的区域内交叉即可。而且，在设定了多个权重函数W的情况下，可以增加与包含斜率变化点的区域对应的权重函数W的数量，并在该区域内密集配置。

另一方面，在3个滑模控制器74～76中，分别计算3个修正系数θ_i(i＝1～3)。这些修正系数θ_i用于修正3个权重函数W_i，如以下所述，通过应用滑模控制算法的算法[式(38)～(44)]来计算。

首先，根据下式(38)分别计算3个分配误差Ew_i(i＝1～3)。

【算式19】

Ew_i(k)＝W_i(k)·DKCMD(k)           ……(38)

如上式(38)所示，由于3个分配误差Ew_i是通过对3个权重函数W_i乘以空燃比校正值DKCMD来计算的，因而3个分配误差Ew_i被计算为将空燃比校正值DKCMD分别分配给排气容积SV_hat的所述第1～第3区域的值。

然后，最终，如下式(39)～(44)所示，通过利用遗忘系数λ的滑模控制算法来计算3个修正系数θ_i。

【算式20】

σw_i(k)＝Ew_i(k)+S·Ew_i(k-1)        ……(39)

θrch_i(k)＝Krch_w·σw_i(k)         ……(40)

θadp_i(k)＝λ·θadp_i(k-1)+Kadp_w·σw_i(k)       ……(41)

     ·θL≦θ_i(k-1)≦θH时

λ＝1                     ……(42)

     ·θ_i(k-1)<θL或θH<θ_i(k-1)时

λ＝λlmt              ……(43)

θ_i(k)＝θrch_i(k)+θadp_i(k)        ……(44)

上式(39)的σw_i是切换函数，S是被设定成使—1<S<0的关系成立的切换函数设定参数。并且，式(40)的θrch_i是趋近律输入，Krch_i表示规定的趋近律增益。

而且，式(41)的θadp_i是自适应律输入，Kadp_i表示规定的自适应律增益。并且，如式(42)、(43)所示，根据修正系数的上次值θ_i(k—1)与规定的上下限值θH、θL的比较结果把该式(41)的遗忘系数λ设定为值1或规定值λlmt。该上限值θH被设定为正的规定值，下限值θL被设定为负的规定值，并且规定值λlmt被设定为0<λlmt<1成立的值。

而且，如式(44)所示，修正系数θ_i被计算为趋近律输入θrch_i与自适应律输入θadp_i的和。

在以上的修正系数θ_i的计算算法中，使用遗忘系数λ是由于以下原因。即，在由于第1催化装置8中的基准空燃比KCMDBS的变化以外的要因(例如，吹扫气体(purge gas)的导入或者含有假定程度以上的乙醇燃料等)而产生大的分配误差Ew_i的情况下，起因于此，修正系数θ_i因为误自适应而暂时成为不恰当的值，过渡控制性下降。

与此相对，在所述的式(41)中，在修正系数的上次值θ_i(k—1)的绝对值大的情况下，为了避免自适应律输入θadp_i的增大，把被设定为0<λ<1的范围内的值的遗忘系数λ乘以自适应律输入的上次值θadp_i(k—1)。在该情况下，当使用递推公式展开所述的式(41)时，对h(h是2以上的自然数)次前的控制定时时的自适应律输入的上次值θadp_i(k—h)乘以λ^h(≒0)，因而即使当进行了运算处理时，也能避免自适应律输入θadp_i的增大。结果，可避免修正系数θ_i因为误自适应而暂时成为不恰当的值，可提高过渡控制性。

并且，在将遗忘系数λ总是设定为0<λ<1的范围内的值的情况下，当修正系数θ_i的区域改变而使得Ew_i＝0时，或者上述要因被消除而使得

时，由于基于遗忘系数λ的遗忘效果而使修正系数θ_i收敛于值0附近，因而当在这种状态下再次产生大的分配误差Ew_i时，为了消除该分配误差Ew_i需要时间。因此，为了避免这一点，并迅速地消除分配误差Ew_i，即使当修正系数θ_i的绝对值较小时，也有必要将修正系数θi恰当地保持为可迅速补偿分配误差Ew_i的值，因而当θL≤θ_i(k—1)≤θH成立时，为了消除基于遗忘系数λ的遗忘效果，设定为λ＝1。另外，在总是不需要基于遗忘系数λ的遗忘效果的情况下，在式(41)中，不管修正系数的上次值θ_i(k—1)的大小如何，只要设定为λ＝1即可。

如上所述，通过利用遗忘系数λ的滑模控制算法来计算3个修正系数θ_i，以使3个权重函数W_i分别对应的、排气容积SV_hat的3个区域中的空燃比校正值DKCMD收敛于值0。

然后，在3个乘法器77～79中，通过把3个修正系数θ_i和3个权重函数W_i分别相乘，来计算3个相乘值θ_iW_i。

进而，在加法器80中，根据下式(45)计算相乘值之和Kff’。

【算式21】

${\mathrm{Kff}}^{,}\left(k\right)=\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{\mathrm{\&theta;}}_i\left(k\right)\&CenterDot;W_i\left(k\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(45\right)$

然后，在加法器81中，根据下式(46)计算模型修正系数Kff。

【算式22】

$\mathrm{Kff}\left(k\right)=1+{\mathrm{Kff}}^{,}\left(k\right)=1+\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{\mathrm{\&theta;}}_i\left(k\right)\&CenterDot;W_i\left(k\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(46\right)$

如上式(46)所示，模型修正系数Kff是通过对相乘值之和Kff’加上值1来计算的。这是因为，当3个修正系数θ_i均收敛于值0时，在式(37)中KCMDADP＝KCMDBS成立。

另外，在本实施方式中，模型修正器70相当于修正值计算单元，修正系数θ_i相当于多个修正值，权重函数W_i相当于规定的多个函数，分配误差Ew_i相当于多个相乘值。

图9示出自适应基准空燃比计算部50中的各种值的计算例。如该图所示，3个相乘值θ_iW_i被计算为表示互不相同的变化的值，并且模型修正系数Kff被计算为相对于排气容积SV_hat表示表现出非线性变化的值。由此，作为使用模型修正系数Kff对基准空燃比KCMDBS进行了非线性修正后的结果，可计算自适应基准空燃比KCMDADP。

如上所述，在自适应基准空燃比计算部50中，计算修正系数θ_i以使空燃比校正值DKCMD成为值0(即，传感器输出VO2成为目标输出VO2_TRGT)，并且在使用这种修正系数θ_i修正权重函数W_i的同时，计算模型修正系数Kff，而且，通过使用该模型修正系数Kff来修正基准空燃比KCMDBS，来计算自适应基准空燃比KCMDADP。由此，即使当第1催化装置8中的实际的基准空燃比KCMDBS由于所述的各种要因而发生非线性变化时，也能在自适应于这种非线性的基准空燃比KCMDBS的变化的同时，计算自适应基准空燃比KCMDADP。

下面，参照图10说明由ECU 2执行的空燃比控制处理。本处理是计算应从燃料喷射阀6喷射的燃料喷射量Tout的处理，按所述规定的控制周期ΔT来执行。

在该处理中，首先，在步骤1(图中简记为“S1”。下同)中，通过根据发动机转速NE和进气压PB来检索未作图示的映射图，来计算基本喷射量Tibs。

然后，进到步骤2，判别传感器故障标志F_SENNG是否是“1”。在未作图示的判定处理中，当空气流量传感器10和进气压传感器11中的至少一方发生故障时，该传感器故障标志F_SENNG被设定为“1”，在除此以外的情况时被设定为“0”。

当步骤2的判别结果是“否”，即2个传感器10、11均是正常时，进到步骤3，判别发动机起动标志F_ENGSTART是否是“1”。该发动机起动标志F_ENGSTART是通过在未作图示的判定处理中，根据发动机转速NE和点火开关(未作图示)的接通/断开信号来判定是否是发动机起动控制中即转动动力输出(cranking)中而设定的，具体地说，当是发动机起动控制中时被设定为“1”，在除此以外的情况时被设定为“0”。

当步骤3的判别结果是“是”，即是发动机起动控制中时，进到步骤4，通过根据发动机水温TW来检索图11所示的映射图，来计算目标空燃比的起动时用值KCMD_ST。在该映射图中，发动机水温TW越低则起动时用值KCMD_ST被设定为越浓侧的值。这是因为，在发动机水温TW低的情况下，为了提高发动机3的起动性能，有必要将混合气控制为浓侧的值。

在继步骤4之后的步骤5中，将目标空燃比KCMD设定为上述起动时用值KCMD_ST。然后，进到步骤6，将燃料喷射量Tout设定为基本喷射量Tibs和目标空燃比KCMD的积Tibs·KCMD，之后结束本处理。

另一方面，当步骤3的判别结果是“否”，即不是发动机起动控制中时，进到步骤7，判别油门开度AP是否小于规定值APREF。当该判别结果是“是”，即未踩下油门踏板时，进到步骤8，判别起动后定时器的计时值Tast是否小于规定值Tastlmt。

当该判别结果是“是”，即Tast<Tastlmt时，应执行催化剂暖机控制，进到步骤9，通过根据起动后定时器的计时值Tast和发动机水温TW来检索图12所示的映射图，来计算目标空燃比的催化剂暖机用值KCMD_AST。在该图中，TW1～TW3表示TW1<TW2<TW3的关系成立的发动机水温TW的规定值。

在该映射图中，在TW＝TW1成立的低水温区域中，计时值Tast越小，为了加快催化剂的活化，则催化剂暖机用值KCMD_AST被设定为越浓侧的大值。并且，在TW＝TW3成立且催化剂暖机结束的高温区域中，催化剂暖机用值KCMD_AST被设定为相当于理论空燃比的值1。

然后，进到步骤10，将目标空燃比KCMD设定为上述催化剂暖机用值KCMD_AST。之后，在步骤11中，使用所述的ST控制器30中的计算方法计算空燃比校正系数KAF。

在继步骤11之后的步骤12中，使用所述的总校正系数计算部32中的计算方法计算总校正系数KTOTAL。然后，进到步骤13，将燃料喷射量Tout设定为基本喷射量Tibs、空燃比校正系数KAF以及总校正系数KTOTAL的积，之后结束本处理。

另一方面，当步骤7或8的判别结果是“否”，即踩下了油门踏板时，或者Tast≥Tastlmt时，进到步骤14，使用所述的空燃比控制器40中的计算方法计算目标空燃比KCMD。

然后，如上所述，执行步骤11～13之后，结束本处理。

另一方面，当步骤2的判别结果是“是”，即2个传感器10、11中的至少一方发生故障时，进到步骤15，将燃料喷射量Tout设定为基本喷射量Tibs和目标空燃比KCMD的规定的故障时用值KFS的积Tibs·KFS，之后结束本处理。

本实施方式的控制装置1通过以上的空燃比控制处理计算燃料喷射量Tout，尽管未作图示，然而根据该燃料喷射量Tout和发动机转速NE计算燃料喷射定时，并根据基于该燃料喷射量Tout和燃料喷射定时的控制输入信号来驱动燃料喷射阀6，从而控制混合气的空燃比。

下面，参照图13和图14说明本实施方式的控制装置1的空燃比控制的模拟结果(以下称为“控制结果”)。首先，图13示出本实施方式的控制装置1的控制结果例，图14为了比较而示出在使用与专利文献1相同的方法计算出自适应基准空燃比KCMDADP的情况下的控制结果例。

在图13(a)中，虚线所示的曲线表示第1催化装置8的上游侧的NOx量，实线所示的曲线表示第1催化装置8的下游侧的NOx量。并且，该图(b)表示第1催化装置8的NOx净化率，该图(c)的KCMDx表示所述的最佳空燃比，即

成立的混合气的空燃比的最佳值。而且，该图(d)的曲线表示目标空燃比与最佳空燃比之间的偏差KCMD—KCMDx，该图(e)的曲线表示目标空燃比与自适应基准空燃比之间的偏差KCMD—KCMDADP。以上的关系在图14(a)～(e)中也是一样。

首先，参照图14可以明白，在使用与专利文献1相同的方法计算出自适应基准空燃比KCMDADP的情况下，当发动机3的NOx排出量增大时，由于自适应基准空燃比KCMDADP对最佳空燃比KCMDx的跟随性低，因而目标空燃比KCMD大幅偏离最佳空燃比KCMDx，输出偏差DVO2增大(氧浓度传感器的输出VO2大幅偏离规定的目标输出VO2_TRGT)，使得目标空燃比KCMD的变动频率下降(时刻t10～t11，t12～t13)。结果得知，第1催化装置8的NOx净化率下降，第1催化装置8的下游侧的NOx量暂时增大。

与此相对，如图13所示，根据本实施方式的控制装置1，在发动机3的NOx排出量增大的情况下，在控制刚开始之后(时刻t1～t2)，自适应基准空燃比KCMDADP相对于最佳空燃比KCMDx产生少许偏差，使得目标空燃比KCMD的变动频率暂时下降。然而，在此以后，伴随控制执行时间的经过，能使自适应基准空燃比KCMDADP自适应于最佳空燃比KCMDx的特性，能使输出偏差DVO2收敛于值0，从而使得目标空燃比KCMD的变动频率被保持为高频带(例如5Hz以上的频带)。结果得知，可确保第1催化装置8的良好的NOx净化率，可减少第1催化装置8的下游侧的NOx量。

如上所述，根据本实施方式的控制装置1，在SM控制器43中，使用自适应滑模控制算法计算空燃比校正值DKCMD，以使O2传感器的输出VO2收敛于规定的目标输出VO2_TRGT，在调制器44中，通过应用Δ∑调制算法的算法来调制该空燃比校正值DKCMD，由此来计算调制值DKCMD_DSM。由此，当

成立时，成立，使得调制值DKCM]D_DSM被计算为表示在高频率(例如5Hz以上的频率)时的反转行为的值。

并且，在自适应基准空燃比计算部50中，分别计算与划分排气容积SV_hat可变化的区域而得到的3个区域对应的3个修正系数θ_i，以使空燃比校正值DKCMD成为值0(即，使传感器输出VO2成为目标输出VO2TRGT)，对权重函数W_i乘以修正系数θ_i，并对这些相乘值W_iθ_i的总和即相乘值之和Kff’加上值1，从而计算模型修正系数Kff。而且，通过使用该模型修正系数Kff来修正基准空燃比KCMDBS，来计算自适应基准空燃比KCMDADP。由此，即使当由于发动机3的运转状态的骤变而使所述的最佳空燃比针对排气容积SV_hat的所述3个区域发生非线性变化时，也能在迅速地自适应于这种非线性的最佳空燃比的变化的同时，计算自适应基准空燃比KCMDADP。

即，即使在伴随发动机3的运转状态的骤变而使图6中的基准空燃比KCMDBS与排气容积SV_hat的关系相对于实际的最佳空燃比与排气容积SV_hat的关系在所述3个区域中的某个区域内产生局部误差的情况下，或者在两者的实际关系在3个区域间发生变动的情况下，也能使自适应基准空燃比KCMDADP与排气容积SV_hat的关系恰当且迅速地与实际的最佳空燃比与排气容积SV_hat的关系一致。

而且，通过对按以上计算出的自适应基准空燃比KCMDADP加上调制值DKCMD_DSM，来计算目标空燃比KCMD，因而即使当由于发动机3的运转状态的骤变而使目标空燃比KCMD暂时大幅偏离最佳空燃比，并且由于传感器输出VO2大幅偏离目标输出VO2_TRGT而使空燃比校正值DKCMD暂时骤增时，也能使用自适应基准空燃比KCMDADP将骤增的空燃比校正值DKCMD迅速控制为值0。由此，由于能将传感器输出VO2迅速控制为目标输出VO2_TRGT，从而能迅速提高调制值DKCMD_DSM的频率即目标空燃比KCMD的变动频率。结果，可迅速消除由于发动机3的运转状态的骤变引起的、目标空燃比KCMD的变动频率的下降，从而可将第1催化装置中的排气净化率维持在高水准。

并且，在所述的对应的区域内3个权重函数W_i分别被设定为值1以下的正值，且在除此以外的区域内分别被设定为值0，而且，相邻的各2个权重函数W_m、W_m+1(m＝1或2)相互交叉，并且交叉部分的两者的和被设定为权重函数W_i的最大值1，因而当修正基准空燃比KCMDBS时，可针对3个区域连续地修正，使得作为修正结果的自适应基准空燃比KCMDADP不会有不连续点。由此，可避免自适应基准空燃比KCMDADP由于基准空燃比KCMDBS的不连续点而暂时成为不恰当的状态。

另外，第1实施方式是将本发明的控制装置应用于把氧浓度传感器15的输出VO2作为控制量并把目标空燃比KCMD作为控制输入的控制对象的例子，然而本发明的控制装置不限于此，作为控制对象，可应用于把各种产业设备中的输出和输入分别作为控制量和控制输入的系统。

并且，第1实施方式是使用自适应滑模控制算法作为规定的反馈控制算法的例子，然而本发明的规定的反馈控制算法不限于此，只要是能将输入值计算成使控制量收敛于目标控制量的算法即可。例如，可以使用通常的滑模控制算法、PID控制算法、后步进控制算法、将滑模控制算法中的控制对象模型置换为一次系数模型后的响应指定型控制算法、最佳调节器等。

而且，第1实施方式是使用应用Δ∑调制算法的调制算法[式(24)～(34)]作为规定的调制算法的例子，然而本发明的规定的调制算法不限于此，只要是针对规定的调制算法的输入值的绝对值越小、则能将输入值的调制值计算成越高频率的算法即可。例如，作为规定的调制算法，可以使用Δ调制算法或∑Δ调制算法、或者基于这些算法的调制算法。

另一方面，第1实施方式是使用排气容积SV_hat作为动作状态参数的例子，然而本发明的动作状态参数不限于此，只要是表示控制对象的动作状态的参数即可。

并且，第1实施方式是使用排气容积SV_hat作为运转状态参数的例子，然而本发明的运转状态参数不限于此，只要是表示内燃机的运转状态的参数即可。例如，可以使用空气流量Gth作为运转状态参数。

并且，第1实施方式是使用图6的一维映射图作为相关关系模型的例子，然而本发明的相关关系模型不限于此，只要是表示第2输入值和动作状态参数的关系的模型即可。例如，可以使用其他计算式或二维以上的映射图等作为相关关系模型。

而且，第1实施方式是使用氧浓度传感器15作为排气浓度传感器的例子，然而本发明的排气浓度传感器不限于此，只要是检测排气中的规定分量的浓度的传感器即可。例如，可以使用检测排气中的NOx浓度的NOx浓度传感器等作为排气浓度传感器。

并且，第1实施方式是将作为控制输入的目标空燃比KCMD计算为作为第2输入值的自适应基准空燃比KCMDADP与作为调制输入值的调制值DKCMD_DSM之和的例子，然而本发明的控制输入不限于此，只要是根据第2输入值和调制输入值之和计算出的控制输入即可。例如，可以通过对第2输入值和调制输入值之和加上规定值，或者对第2输入值和调制输入值之和乘以规定值来计算控制输入。

下面，说明本发明的第2实施方式的控制装置1A(参照图15)。另外，在以下说明中，关于与第1实施方式相同的结构，附上相同标号，并省略其说明。该控制装置1A是应用于所谓的FR方式的车辆(未作图示)的控制装置，具体地说，是执行该车辆的牵引控制的控制装置，该FR方式的车辆将所述的具有有档位式自动变速器的发动机3安装在前侧，并将均未作图示的后轮和前轮分别用作驱动轮和非驱动轮(从动轮)。

另外，牵引控制是指以下的控制方法，即：在车辆加速时，当发动机3的产生转矩过大、且发生驱动轮相对于非驱动轮空转的状态时，通过降低发动机3的产生转矩来避免空转状态，从而确保车辆的稳定性并提高加速性能。

如图15所示，该控制装置1A具有ECU 2，该ECU 2除了与所述的传感器10～17等连接以外，还连接有左右的前轮速传感器20、21、左右的后轮速传感器22、23、以及档位传感器24。左右的前轮速传感器20、21分别将表示左右的前轮速的检测信号输出到ECU 2，左右的后轮速传感器22、23分别将表示左右的后轮速的检测信号输出到ECU 2。

ECU 2根据左右的前轮速传感器20、21的检测信号计算左右的前轮速，并计算它们的相加平均值作为非驱动轮速Ws_ref。而且，ECU2根据左右的后轮速传感器22、23的检测信号计算左右的后轮速，并计算它们的相加平均值作为驱动轮速Ws_act。

另外，在本实施方式中，ECU 2相当于：控制量检测单元、目标控制量设定单元、第1输入值计算单元、调制输入值计算单元、动作状态参数检测单元、误差参数计算单元、修正值计算单元、模型修正单元、第2输入值计算单元、控制输入计算单元、设定单元、估计值计算单元、误差参数设定单元、要求输出计算单元以及选择单元。而且，左右的后轮速传感器22、23相当于控制量检测单元，驱动轮速Ws_act相当于控制量和驱动轮的速度，曲轴转角传感器13相当于动作状态参数检测单元，发动机转速NE相当于动作状态参数和运转状态参数。

另一方面，有档位式自动变速器具有6个前进档位和1个倒车档位，档位传感器24检测该有档位式自动变速器的档位，将表示该档位的检测信号输出到ECU 2。ECU 2根据该档位传感器24的检测信号设定前进档位值NGEAR。具体地说，该前进档位值NGEAR在1～3前进档位时分别被设定为值1～3，在4～6前进档位时被设定为值4。

并且，如图16所示，控制装置1A具有牵引控制器100。该牵引控制器100如下所述，是计算被估计为能够避免驱动轮的空转状态、并可确保车辆的稳定性和加速性的发动机3的转矩、即发动机转矩TRQENG的控制器，具体地说由ECU 2构成。另外，在本实施方式中，牵引控制器100相当于模型修正单元、第2输入值计算单元以及控制输入计算单元，发动机转矩TRQENG相当于内燃机的输出目标。

如该图所示，牵引控制器100具有：目标车轮速计算部101、车轮速反馈控制器102、调制器103、要求转矩计算部104、基准前馈转矩计算部105、估计转矩计算部106、乘法器107、前馈转矩计算部108、加法器109、发动机转矩计算部110以及模型修正器120。

首先，在目标车轮速计算部101中，根据下式(47)计算目标车轮速Ws_cmd。另外，在本实施方式中，目标车轮速计算部101相当于目标控制量设定单元，目标车轮速Ws_cmd相当于目标控制量和目标速度。

【算式23】

Ws_cmd(k)＝Ws_ref(k)+OptSlip        ……(47)

上式(47)的OptSlip是与在驱动轮和非驱动轮之间可容许的滑移(slip)量相当的规定的滑移偏置值，在本实施方式中，被设定为一定值(例如10km/h)。在该情况下，可以根据规定的参数(例如，非驱动轮速Ws_ref、路面的摩擦阻力系数的估计值、横摆率传感器的检测信号以及车身的滑移角传感器的检测信号等)，通过映射图检索或使用规定的计算式来决定滑移偏置值OptSlip。另外，在本实施方式中，滑移偏置值OptSlip相当于规定量。

下面，说明所述的车轮速反馈控制器102。在该车轮速反馈控制器102中，如下所述，使用将目标值滤波型2自由度滑模控制算法和自适应干扰观测器组合而应用的控制算法[式(48)～(58)]来计算转矩反馈值TRQFB。另外，在本实施方式中，车轮速反馈控制器102相当于第1输入值计算单元，转矩反馈值TRQFB相当于第1输入值。

首先，使用下式(48)所示的一次延迟类型的低通滤波算法计算目标车轮速的滤波值Ws_cmd_f。

【算式24】

Ws_cmd_f(k)＝-Rt·Ws_cmd_f(k-1)+(1+Rt)Ws_cmd(k)      ……(48)

这里，上式(48)的Rt是目标值滤波设定参数，被设定成使—1<Rt<0的关系成立的值。在该情况下，使用目标值滤波设定参数Rt的设定值决定滤波值Ws_cmd_f向目标车轮速Ws_cmd的跟随速度。

然后，使用式(49)～(52)所示的应用滑模控制算法的控制算法来计算趋近律输入Urch_t和非线性输入Unl_t。

【算式25】

Et(k)＝Ws_act(k)-Ws_cmd_f(k)           ……(49)

σt(k)＝Et(k)+St·Et(k-1)                  ……(50)

Urch_t(k)＝-Krch_t·σt(k)                 ……(51)

Unl_t(k)＝-Knl_t·sgn(σt(k))              ……(52)

这里，式(49)的Et是跟随误差，式(50)的σt是切换函数。并且，式(50)的St是切换函数设定参数，被设定成使—1<St<0的关系成立的值。在该情况下，使用切换函数设定参数St的设定值来指定跟随误差Et向值0的收敛速度。而且，式(51)的Krch_t表示规定的趋近律增益，式(52)的Knl_t表示规定的非线性输入用增益。而且，式(52)中的sgn(σt(k))是符号函数，其值被设定成，当σt(k)≥0时为sgn(σt(k))＝1，并当σt(k)<0时为sgn(σt(k))＝—1(另外，也可以设定成，当σt(k)＝0时为sgn(σt(k))＝0)。

然后，使用式(53)～(57)所示的应用自适应干扰观测器的控制算法来计算干扰估计值Uls_t。

【算式26】

σt_hat(k)＝Urch_t(k-1)+Unl_t(k-1)+Uls_t(k-1)       ……(53)

Et_sig(k)＝σt(k)-σt_hat(k)

         ＝σt(k)-Urch_t(k-1)-Un_t(k-1)-Uls_t(k-1)     ……(54)

$\mathrm{Uls}\_t\left(k\right)=\mathrm{\&lambda;t}\&CenterDot;\mathrm{Uls}\_t(k-1)+\frac{\mathrm{Pt}}{1+\mathrm{Pt}}\mathrm{Et}\_\mathrm{sig}\left(k\right).....\left(55\right)$

     ·Uls_t_L<Uls_t(k-1)<Uls_t_H时

λt＝1　　　　　　　　　　　　……(56)

     ·Uls_t(k-1)≦Uls_t_L或Uls_t_H≦Uls_t(k-1)时

λt＝λtlmt                  ……(57)

这里，式(53)的σt_hat是切换函数的估计值，Uls_t是干扰估计值。该干扰估计值Uls_t是使用式(54)、(55)所示的固定增益式的同定算法来计算的。该式(54)的Et_sig表示估计误差，式(55)的Pt是一定值的同定增益。

而且，式(55)的λt是遗忘系数，其值如式(56)、(57)所示，根据干扰估计值的上次值Uls_t(k—1)与规定的上下限值Uls_t_H、Uls_t_L的比较结果被设定为值1或规定值λtlmt。该上限值Uls_t_H被设定为正的规定值，下限值Uls_t_L被设定为负的规定值，并且规定值λtlmt被设定为使0<λtlmt<1成立的值。

然后，最终如下式(58)所示转矩反馈值TRQFB被计算为趋近律输入Urch_t、非线性输入Unl_t以及干扰估计值Uls_t的和。

【算式27】

TRQFB(k)＝Urch_t(k)+Unl_t(k)+Uls_t(k)         ……(58)

如上所述，在车轮速反馈控制器102中，由于转矩反馈值TRQFB是使用式(48)～(58)所示的控制算法来计算的，因而该转矩反馈值TRQFB被计算为用于使驱动轮速Ws_act收敛于目标车轮速的滤波值Ws_cmd_f的值，即用于使驱动轮速Ws_act收敛于目标车轮速Ws_cmd的值。在该情况下，如上所述由于目标车轮速Ws_cmd是通过对非驱动轮速Ws_ref加上滑移偏置值OptSlip来计算的，因而在

的状态下，变成

成立的状态。

并且，由于在使用遗忘系数λt的同时，计算转矩反馈值TRQFB，因而在干扰估计值的上次值Uls_t(k—1)的绝对值大的情况下，可避免由于所述的遗忘效果而使干扰估计值Uls_t即转矩反馈值TRQFB积分式地增大，从而可确保控制系统的过渡响应的稳定性。并且，在干扰估计值的上次值Uls_t(k—1)的绝对值小的情况下，由于遗忘系数λt被设定为值1，因而即使当跟随误差Et收敛于值0时，也能将转矩反馈值TRQFB保持为可迅速补偿跟随误差Et的恰当的值，由此，可提高跟随误差Et开始增大时的响应性。

下面，说明所述的调制器103。在该调制器103中，通过使用以下所述的调制算法来调制转矩反馈值TRQFB，来计算调制值TRQFB_DSM。首先，根据下式(59)～(61)计算小振幅分量值TRQFB_L。

【算式28】

     ·TRQFB(k)<-TRQFB_L_LMT时

TRQFB_L(k)＝-TRQFB_L_LMT          ……(59)

     ·-TRQFB_L_LMT≦TRQFB(k)≦TRQFB_L_LMT时

TRQFB_L(k)＝TRQFB(k)               ……(60)

     ·TRQFB_L_LMT<TRQFB(k)时

TRQFB_L(k)＝TRQFB_L_LMT             ……(61)

上式(59)～(61)的TRQFB_L_LMT是规定的阈值，如这些式(59)～(61)所示，小振幅分量值TRQFB_L是通过对转矩反馈值TRQFB的振幅实施以值TRQFB_L_LMT为上限并以值—TRQFB_L_LMT为下限的限制处理来计算的。即，小振幅分量值TRQFB_L相当于当转矩反馈值TRQFB的变动幅度小、且其绝对值处于不超过阈值TRQFB_L_LMT的范围内时的转矩反馈值TRQFB。

并且，大振幅分量值TRQFB_H根据下式(62)～(64)来计算。

【算式29】

     ·TRQFB(k)<-TRQFB_L_LMT时

TRQFB_H(k)＝TRQFB(k)+TRQFB_L_LMT             ……(62)

     ·-TRQFB_L_LMT≦TRQFB(k)≦TRQFB_L_LMT时

TRQFB_H(k)＝0                             ……(63)

     ·TRQFB_L_LMT<TRQFB(k)时

TRQFB_H(k)＝TRQFB(k)-TRQFB_L_LMT             ……(64)

如以上的式(62)～(64)所示，当转矩反馈值TRQFB的绝对值不超过阈值TRQFB_L_LMT时，大振幅分量值TRQFB_H被计算为值0，当转矩反馈值TRQFB的绝对值超过阈值TRQFB_L_LMT时，大振幅分量值TRQFB_H被计算为超过部分的值。即，在由于转矩反馈值TRQFB的变动大而要求控制的快速响应性的情况下，大振幅分量值TRQFB_H被计算为用于把这种转矩反馈值TRQFB恰当地反映在发动机转矩TRQENG上的值。

而且，通过使用下式(65)～(68)所示的应用Δ∑调制算法的调制算法来调制小振幅分量值TRQFB_L，来计算调制分量值TRQFB_L_DSM。

【算式30】

δdsm2(k)＝TRQFB_L(k)-TRQFB_L_DSM(k-1)            ……(65)

σdsm2(k)＝σdsm2(k-1)+δdsm2(k)                    ……(66)

    ·σdsm2(k)<0时

TRQFB_L_DSM(k)＝-TRQFB_L_AMP                   ……(67)

     ·σdsm2(k)≧0时

TRQFB_L_DSM(k)＝TRQFB_L_AMP                    ……(68)

如上式(65)所示，偏差δdsm2被计算为小振幅分量值TRQFB_L与调制分量值TRQFB_L_DSM的上次值之间的偏差。并且，在式(66)中，σdsm2表示偏差δdsm2的积分值。并且，式(68)的TRQFB_L_AMP是规定的振幅值。参照以上的式(65)～(68)可以明白，调制分量值TRQFB_L_DSM被计算为以值0为中心在最小值—TRQQB_L_AMP与最大值TRQFB_L_AMP之间重复反转的值。

然后，最终如下式(69)所示，通过对调制分量值TRQFB_L_DSM加上大振幅分量值TRQFB_H，来计算调制值TRQFB_DSM。

【算式31】

FRQFB_DSM(k)＝TRQFB_L_DSM(k)+TRQFB_H(k)             ……(69)

另外，在本实施方式中，调制器103相当于调制输入值计算单元，调制值DKCMD_DSM相当于调制输入值。

下面，说明所述的要求转矩计算部104。该要求转矩计算部104使用以下所述的方法计算要求转矩TRQDRV，该要求转矩TRQDRV相当于驾驶者对发动机3所要求的转矩。另外，在本实施方式中，要求转矩计算部104相当于要求输出计算单元，要求转矩TRQDRV相当于要求输出。

首先，通过根据发动机转速Ne来检索图17所示的映射图，来分别计算最大转矩Trq_max和最小转矩Trq_min。该图的NEhigh表示规定的最大容许转速(例如7000rpm)。这些值Trq_max、Trq_min相当于在该发动机转速NE时可达到的发动机转矩的最大值和最小值。并且，在该映射图中，最小转矩Trq_min被设定为负值。这是因为，最小转矩Trq_min相当于在未踩下油门踏板、且减速燃料切断运转中的发动机制动状态下的发动机转矩。

并且，通过根据油门开度AP来检索图18所示的映射图，来计算归一化要求驱动力Ktrq_ap。该图的APmax表示油门开度的最大值(100％)。并且，归一化要求驱动力Ktrq_ap表示以AP＝APmax时的要求驱动力Trq_apmax为基准将根据油门开度AP所决定的要求驱动力Trq_ap归一化后的值，即Ktrq_ap＝Trq_ap÷Trq_apmax成立的值。

然后，最终根据下式(70)计算要求转矩TRQDRV。

【算式32】

TRQDRV(k)＝Ktrq_ap(k){Trq_max(k)-Ttrq_min(k)}+Ttrq_min(k)

                                             ……(70)

另一方面，在所述的基准前馈转矩计算部105中，通过根据发动机转速NE和前进档位值NGEAR来检索图19所示的映射图，来计算基准前馈转矩TRQFFM。该基准前馈转矩TRQFFM相当于在代表性的路面状况和轮胎状况下，驱动轮相对于非驱动轮不会产生滑移、且成立的发动机3的转矩的最大值。

并且，在所述的估计转矩计算部106中，估计转矩TRQINST是使用以吸入空气量Gcyl、进气压PB、检测空燃比KACT、EGR率、燃料喷射量Tout和点火正时等作为输入并以估计转矩TRQINST作为输出的神经网络模型(未作图示)来计算的。该估计转矩TRQINST是估计出在该时刻

成立的发动机3的转矩的上限值而得到的。

在该情况下，估计转矩计算部106中的估计转矩TRQINST的计算方法不限于上述的计算方法，可以使用映射图检索方法来取代神经网络模型。并且，也可以使用检测气缸内的压力的缸内压传感器，根据其检测信号计算图示平均有效压力IMEP，并根据该图示平均有效压力IMEP和预先设定的发动机3的摩擦数据来对估计转矩TRQINST进行估计。而且，也可以把所述的转矩反馈值TRQFB用作估计转矩TRQINST。

另外，在本实施方式中，估计转矩计算部106相当于估计值计算单元，估计转矩TRQINST相当于控制输入的估计值。

下面，参照图20说明所述的模型修正器120。该模型修正器120如以下所述是计算模型修正系数KTRQFF的模型修正器，该模型修正系数KTRQFF如后所述用于修正(校正)基准前馈转矩TRQFFM。另外，在本实施方式中，模型修正器120相当于误差参数计算单元和修正值计算单元。

如该图所示，模型修正器120具有：权重函数计算部121、转矩误差计算部122、档位系数设定部123、修正转矩误差计算部124、基准局部校正值计算部125以及修正系数计算部126。

首先，在权重函数计算部121中，如以下所述，根据发动机转速NE和前进档位值NGEAR来计算权重函数WNE_ij。这里，权重函数WNE_ij的下标i(i＝1～3)表示以下所述的发动机转速NE的3个区域，下标j(j＝1～4)表示前进档位值NGEAR可取的所有的值1～4。即，权重函数WNE_ij被计算为以12个值作为要素的矢量，这些要素分别对应于通过发动机转速NE和前进档位值NGEAR的组合而决定的12个区域。

该权重函数计算部121具有与前进档位值NGEAR对应的4个映射图作为在权重函数WNE_ij的计算中使用的映射图。权重函数WNE_ij是这样来计算的，即：首先，从4个映射图中选择与前进档位值NGEAR对应的映射图，然后根据发动机转速NE检索所选择的映射图。

图21是示出在4个映射图中的NGEAR＝j＝3时使用的映射图，即当自动变速器为前进档位3时所选择的映射图的一例。在该图中，NEx1～6是发动机转速NE的规定值，被设定成使NEx1<NEx2<NEx3<NEx4<NEx5<NEx6(＝NEhigh)的关系成立。

如该图所示，在将发动机转速NE可变化的区域划分为0≤NE<NEx3、NEx1<NE<NEx5、NEx3<NE这3个区域的情况下，3个权重函数WNE_ij分别被设定成与这3个区域对应，并且在对应的区域内被设定为值1以下的正值，在对应区域以外被设定为值0。

具体地说，权重函数WNE_1j在其对应的区域(0≤NE<NEx3)内，以NE≤NEx1时的值1作为最大值，被设定为发动机转速NE越高则越小的正值，并在除此以外的区域内被设定为值0。并且，权重函数WNE_2j在其对应的区域(NEx1<NE<NEx5)内，以NE＝NEx3时的值1作为最大值，被设定为呈三角形的斜边状变化的值，并在除此以外的区域内被设定为值0。

而且，权重函数WNE_3j在其对应的区域(NEx3<NE)内，以NEx5≤NE时的值1作为最大值，被设定为发动机转速NE越高则越大的正值，并在除此以外的区域内被设定为值0。

除了以上之外，3个权重函数WNE_ij各自对应的3个区域被设定为上述那样的相邻区域相互重叠的区域，与这些相互重叠的区域对应的权重函数WNE_ij的值的和被设定成等于各权重函数WNE_ij中的最大值1。例如，在NE＝NEx2的情况下，与该值NEx2对应的2个权重函数WNE_1j、WNE_2j的值分别被设定为0.5，这些权重函数的和WNE_1j+WNE_2j为等于权重函数WNE_ij各自的最大值的值1。而且，在NE＝NEx4的情况下，与规定值NEx4对应的2个权重函数WNE_2j、WNE_3j的和WNE_2j+WNE_3i为等于权重函数WNE_ij各自的最大值的值1。

在该情况下，作为权重函数WNE_ij，可以取代在图21所示的1个映射图中以3个值作为要素的权重函数WNE_ij，而使用在1个映射图中以2个或4个以上的值作为要素的权重函数WNE_ij。在该情况下，可以根据要素的数量将发动机转速NE的区域设定成相互重叠。另外，在本实施方式中，权重函数WNE_ij相当于多个函数。

并且，在转矩误差计算部122中，根据下式(71)～(74)计算转矩误差Etrq。

【算式33】

DT(k)＝TRQFF_LMT(k-1)-TRQNST(k)　　　　　　　　……(71)

     ·DT(k)<ETRQ_L时

Etrq(k)＝ETRQ_L　　　　　　　　　　　　……(72)

     ·ETRQ_L≦DT(k)≦ETRQ_H时

Etrq(k)＝DT(k)　　　　　　　　……(73)

      ·ETRQ_H<DT(k)时

Etrq(k)＝ETRQ_H　　　　　　　　……(74)

上式(71)的DT是偏差，该式(71)的TRQFF_LMT是如后所述计算的自适应限制转矩。并且，式(72)的ETRQ_L是规定的下限值，式(73)的ETRQ_H是规定的上限值。如上所述，转矩误差Etrq是通过对偏差DT实施以2个值ETRQ_L、ETRQ_L作为上下限值的限制处理来计算的。

另外，在本实施方式中，转矩误差计算部122相当于误差参数设定单元，转矩误差Etrq相当于误差参数和规定的目标值中的一方与另一方之间的偏差、以及估计值与第2输入值之间的偏差，值0相当于规定的目标值。

并且，在档位系数设定部123中，根据下式(75)、(76)分别设定4个档位系数δ_j。

【算式34】

    ·NGEAR<k)＝j时

δ_j(NGEAR(k))＝1            ……(75)

    ·NGEAR(k)≠j时

δ_j(NGEAR(k))＝0                 ……(76)

参照上式(75)、(76)可以明白，4个档位系数δ_j中，只有其下标j等于前进档位值NGEAR的档位系数δ_j才被设定为值1，除此以外被设定为值0。例如，当NGEAR＝1时，被设定为δ₁＝1、δ₂～δ₄＝0，当NGEAR＝3时，被设定为δ₃＝1、δ₂＝δ₂＝δ₄＝0。

然后，在修正转矩误差计算部124中，根据下式(77)计算修正转矩误差Wetrq_ij。

【算式35】

Wetrq_ij(k)＝WNE_ij(k)·δ_j(NGEAR(k))·Etrq(k)      ……(77)

参照上式(77)可以明白，修正转矩误差Wetrq_ij被计算为以12个值作为要素的矢量，并且在这些要素中，下标j≠NGEAR的要素全部为值0。另外，在本实施方式中，修正转矩误差Wetrq_ij相当于多个相乘值。

然后，在基准局部校正值计算部125中，使用下式(78)～(84)所示的利用遗忘系数λv的滑模控制算法来计算基准局部校正值DKTRQFF_BS_ij。即，基准局部校正值DKTRQFF_BS_ij被计算为以与通过发动机转速NE和前进档位值NGEAR的组合而决定的12个区域对应的12个值为要素的矢量。

【算式36】

σv_ij(k)＝Wetrq_ij(k)+Sv·Wetrq_ij(k-1)         ……(78)

Urch_v_ij(k)＝-Krch_v·σv_ij(k)                ……(79)

Unl_v_ij(k)＝-Knl_v·sgn(σv_ij(k))             ……(80)

【算式37】

Uadp_v_ij(k)＝λv·Uadp_v_ij(k-1)+Kadp_v·σv_ij(k)       ……(81)

     ·DKTRQ_L<DKTRQFF_BS_ij(k—1)<DKTRQ_H时

λv＝1                              ……(82)

·DKTRQFF_BS_ij(k-1)≦DKTRQ_L

             DKTRQ_H≦DKTRQFF_BS_ij(k-1)时

λv＝λvlmt                           ……(83)

【算式38】

DKTRQFF_BS_ij(k)＝Urch_v_ij(k)+Unl_v_ij(k)+Uadp_v_ij(k)                       ……(84)

在上式(78)中，σv_ij是切换函数，Sv是被设定成使—1<Sv<0的关系成立的切换函数设定参数。在该情况下，使用切换函数设定参数Sv的设定值来指定修正转矩误差Wetrq_ij向值0的收敛速度。并且，在式(79)中，Urch_v_ij是趋近律输入，Krch_v表示规定的趋近律增益。而且，在式(80)中，Unl_v_ij是非线性输入，Knl_v表示规定的非线性输入用增益。并且，式(80)中的sgn(σv_ij(k))是符号函数，其值被设定成，当σv_ij(k)≥0时为sgn(σv_ij(k))＝1，并当σv_ij(k)<0时为sgn(σv_ij(k))＝—1(另外，也可以设定成，当σv_ij(k)＝0时为sgn(σv_ij(k))＝0)。

而且，在式(81)中，Uadp_v_ij是自适应律输入，Kadp_v表示规定的自适应律增益。并且，式(81)的遗忘系数λv如式(82)、(83)所示，根据基准局部校正值的上次值DKTRQFF_BS_ij(k—1)与规定的上下限值DKTRQ_H、DKTRQ_L的比较结果被设定为值1或规定值λvlmt。该上限值DKTRQ_H被设定为正的规定值，下限值DKTRQ_L被设定为负的规定值，并且规定值λvlmt被设定为使0<λvlmt<1成立的值。

使用这种遗忘系数λv是由于以下原因。即，在由于局部的路面变化(例如，路面的局部冻结或局部的路面铺装状态的变化等)而使路面的抓地状态暂时骤变的情况下，转矩误差Etrq的绝对值骤增，伴随于此，自适应律输入Uadp_v_ij的绝对值骤增，后述的自适应限制转矩TRQFF_LMT也骤减或骤增。当在该状态下局部的路面变化消除时，如不使用遗忘系数λv，则自适应律输入Uadp_v_ij被维持在大的值上，使得自适应限制转矩TRQFF_LMT大幅偏离最佳的规定范围，结果，产生加速响应的下降或瞬间的车轮滑移。因此，在本实施方式中，使用遗忘系数λv，以便于当局部的路面变化消除时，使自适应限制转矩TRQFF_LMT迅速恢复到规定范围内。

并且，如式(84)所示，基准局部校正值DKTRQFF_BS_ij被计算为趋近律输入Urch_v_ij、非线性输入Unl_v_ij以及自适应律输入Uadp_v_ij的和。

如上所述，在该基准局部校正值计算部125中，由于基准局部校正值DKTRQFF_BS_ij是使用式(78)～(84)所示的控制算法来计算的，因而该基准局部校正值DKTRQFF_BS_ij分别被计算为用于使修正转矩误差Wetrq_ij收敛于值0的值，换句话说，用于使自适应限制转矩TRQFF_LMT收敛于估计转矩TRQINST的值。

在该情况下，在基准局部校正值DKTRQFF_BS_ij的计算中使用的控制算法不限于所述的式(78)～(84)，只要是能将Wetrq_ij分别收敛于值0的反馈控制算法即可。例如，可以使用PID控制算法、反向步进控制算法、将滑模控制算法中的控制对象模型置换为一次系统的模型后的响应指定型控制算法、最佳调节器等。另外，在本实施方式中，基准局部校正值DKTRQFF_BS_ij相当于多个修正值。

并且，在所述的修正系数计算部126中，根据下式(85)～(87)计算模型修正系数KTRQFF。

【算式39】

Dktrqff_ij(k)＝WNE_ij(k)·δ_j(NGEAR(k))·DKTRQFF_BS_ij(k)               ……(85)

${\mathrm{DKTRQFF}}_j=\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{\mathrm{Dktrqff}}_{\mathrm{ij}}\left(k\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(86\right)$

$\mathrm{KTRQFF}\left(k\right)=1+\underset{j=1}{\overset{4}{\mathrm{\&Sigma;}}}{\mathrm{DKTRQFF}}_j\left(k\right)\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\&CenterDot;\left(87\right)$

上式(85)的Dktrqff_ij表示相乘值，式(86)的DKTRQFF_j表示局部校正值。在该情况下，由于相乘值Dktrqff_ij是通过乘以档位系数δ_j来计算的，因而在4个局部校正值DKTRQFF_j中，下标j≠NGEAR成立的3个局部校正值DKTRQFF_j均被计算为值0。

如上所述，在该模型修正器120中，模型修正系数KTRQFF是通过对4个局部校正值DKTRQFF_j的总和加上值1来计算的。这是因为，如以下所述，由于模型修正系数KTRQFF用作对基准前馈转矩TRQFFM的相乘值，因而在无需校正基准前馈转矩TRQFFM的情况下，使KTRQFF＝1。

即，在所述的乘法器107中，根据下式(88)计算自适应限制转矩TRQFF_LMT。另外，在本实施方式中，自适应限制转矩TRQFF_LMT相当于第2输入值。

【算式40】

TRQFF_LMT(k)＝KTRQFF(k)·TRQFFM(k)            ……(88)

图22示出以上的模型修正器120和乘法器107中的各种值的计算结果例，特别是当前进档位值NGEAR＝j＝3时即自动变速器为前进档位3时的计算结果例。如该图所示，3个相乘值Dktrqff_ij被计算为表示相互不同的变化的值，并且模型修正系数KTRQFF被计算为相对于发动机转速NE表现出非线性变化的值。由此，作为使用模型修正系数KTRQFF对基准前馈转矩TRQFFM进行了非线性修正后的结果，可计算自适应限制转矩TRQFF_LMT。

如上所述，模型修正器120和乘法器107中，计算基准局部校正值DKTRQFF_BS_ij以使转矩误差Etrq成为值0(即，驱动轮速Ws_act成为目标车轮速Ws_cmd)，并且在使用这种基准局部校正值DKTRQFF_BS_ij来修正权重函数WNE_ij的同时，计算模型修正系数KTRQFF，而且，通过使用该模型修正系数KTRQFF修正基准前馈转矩TRQFFM，来计算自适应限制转矩TRQFF_LMT。由此，即使当不会产生驱动轮的滑移的最佳的转矩(以下称为“最佳转矩”)由于所述的各种要因而发生非线性变化时，也能在自适应于这种非线性变化的同时，计算自适应限制转矩TRQFF_LMT。

并且，在所述的前馈转矩计算部108中，根据下式(89)、(90)计算前馈转矩TRQFF。

【算式41】

     ·TRQFF_LMT(k)<TRQDRV(k)时

TRQFF(k)＝TRQFF_LMT(k)              ……(89)

     ·TRQFF_LMT(k)≧TRQDRV(k)时

TRQFF(k)＝TRQDRV(k)                 ……(90)

如上式(89)、(90)所示，前馈转矩TRQFF被设定为自适应限制转矩TRQFF_LMT和要求转矩TRQDRV中的小的一方。其原因在后面描述。

然后，在所述的加法器109中，如下式(91)所示，通过对调制值DKCMD_DSM加上前馈转矩TRQFF，来计算限制转矩TRQTC。

【算式42】

TRQTC(k)＝TRQFB_DSM(k)+TRQFF(k)          ……(91)

然后，在所述的发动机转矩计算部110中，根据下式(92)、(93)最终计算发动机转矩TRQENG。

【算式43】

    ·TRQTC(k)<TRQDRV(k)时

TRQENG(k)＝TRQTC(k)                     ……(92)

    ·TRQTC(k)≧TRQDRV(k)时

TRQRNG(k)＝TRQDRV(k)                    ……(93)

如上式(92)、(93)所示，发动机转矩TRQENG被设定为限制转矩TRQTC和要求转矩TRQDRV中的小的一方。另外，在本实施方式中，发动机转矩计算部110相当于选择单元，限制转矩TRQTC相当于目标输出。

如上所述，在该牵引控制器100中，在发动机转矩计算部110中把2个转矩TRQTC、TRQDRV中的小的一方选择为发动机转矩TRQENG，并且在前馈转矩计算部108中把2个转矩TRQFF_LMT、TRQDRV中的小的一方选择为前馈转矩TRQFF。这是由于以下原因。

即，在驾驶者要求的要求转矩TRQDRV过大且驱动轮很有可能滑移的情况下，在发动机转矩计算部110中，通过对要求转矩TRQDRV实施以限制转矩TRQTC为上限的限制处理，来计算发动机转矩TRQENG，因而可将发动机3产生的转矩恰当地控制为可避免驱动轮的滑移的限制转矩TRQTC，由此可确保运行稳定性。

另一方面，在驾驶者要求的要求转矩TRQDRV小于限制转矩TRQTC的情况下，在发动机转矩计算部110中把发动机转矩TRQENG设定为要求转矩TRQDRV，从而可避免驱动轮的滑移。在该情况下，在车轮速反馈控制器102中，有时转矩反馈值TRQFB为负值。在该状态下，由于驾驶者的急剧的油门踏板操作而使要求转矩TRQDRV骤增并大幅超过限制转矩TRQTC的情况下，在发动机转矩计算部110中，尽管把发动机转矩TRQENG设定为限制转矩TRQTC，然而转矩反馈值TRQFB是使用反馈控制算法来计算的，因而其具有响应性低的特性。因此，当转矩反馈值TRQFB为负值时，由于限制转矩TRQTC被抑制为过度小的值，发动机3产生的转矩被过度地抑制，从而具有车辆的加速响应性能下降、运转性能下降的可能性。

与此相对，在前馈转矩计算部108中，在要求转矩TRQDRV骤增并大幅超过自适应限制转矩TRQFF_LMT的情况下，把前馈转矩TRQFF设定为自适应限制转矩TRQFF_LMT。该自适应限制转矩TRQFF_LMT被计算为可避免驱动轮的滑移的转矩上限值，并且通过对自适应限制转矩TRQFF_LMT加上对转矩反馈值TRQFB进行了调制后的调制值TRQFB_DSM，来计算限制转矩TRQTC，因而可将发动机转矩TRQENG计算为以自适应限制转矩TRQFF_LMT为中心重复反转的值。由此，可避免将限制转矩TRQTC抑制为过度小的值，可将发动机3产生的转矩保持在可避免驱动轮的滑移的转矩上限值附近，从而可获得良好的加速响应性。因此，为了获得以上的效果，在2个计算部108、110中，使用所述的方法计算发动机转矩TRQENG和前馈转矩TRQFF。

在本实施方式的控制装置1A中，如上所述，由牵引控制器100计算发动机转矩TRQENG，并且尽管未作图示，以该发动机转矩TRQENG为目标，执行吸入空气量控制处理、空燃比控制处理以及点火正时控制处理，以使发动机3的实际转矩成为该发动机转矩TRQENG。

下面，参照图23、24说明第2实施方式的控制装置1A的牵引控制的模拟结果(以下称为“控制结果”)。图23示出在摩擦阻力小的路面上重复进行车辆的加速/减速时的本实施方式的控制装置1A的控制结果例，特别是当j＝3且自动变速器为前进档位3时的控制结果例。并且，图24为了比较而示出在把模型修正系数KTRQFF保持为值1的情况下，即直接使用基准前馈转矩TRQFFM作为自适应限制转矩TRQFF_LMT的情况下的控制结果例(以下称为“比较例”)。

首先，参照图24，在该比较例中，由于基准前馈转矩TRQFFM用作自适应限制转矩TRQFF_LMT，因而仅使用转矩反馈值TRQFB来校正限制转矩TRQTC，结果得知，在驱动轮速Ws_act达到目标车轮速Ws_cmd的时刻(时刻t31、t33、t35、t37)以后，驱动轮速Ws_act相对于目标车轮速Ws_cmd发生过冲。并且，由于转矩反馈值TRQFB不收敛于值0，因而调制值TRQFB_DSM的变动频率容易下降，不能对滑移状态的驱动轮赋予微小振动，从而不能补偿所述的抓地行为的非线性特性，结果得知，要使驱动轮的滑移量的恒定偏差(驱动轮速Ws_act相对于目标车轮速Ws_cmd的恒定偏差)收敛于值0需要时间。

与此相对，如图23所示，根据本实施方式的控制装置1A得知，当要求转矩TRQDRV变化为发生滑移的过大值时，在第1次牵引控制的执行期间(时刻t21～t22)中，尽管发生少许过冲，然而该程度比起比较例来已得到抑制，而且，在第2次及第2次以后的牵引控制的执行期间(时刻t23～t24、t25～t26、t27～t28)中，驱动轮速Ws_act收敛于目标车轮速Ws_cmd而不发生过冲。这是因为，通过多次执行牵引控制，来推进模型修正系数KTRQFF的学习，自适应限制转矩TRQFF_LMT自适应于不会发生驱动轮的滑移的最佳转矩，能够使用利用该自适应限制转矩TRQFF_LMT计算出的限制转矩TRQTC来恰当地限制发动机转矩TRQENG。

并且，由于能够将转矩反馈值TRQFB保持在值0附近，从而可将调制值TRQFB_DSM的变动频率保持在高频带。可知能够由此对滑移状态的驱动轮赋予微小振动，可补偿所述的抓地行为的非线性特性，从而可抑制驱动轮的滑移量的恒定偏差。可知，如上所述本实施方式的控制结果例与比较例相比控制精度提高。

如上所述，根据第2实施方式的控制装置1A，计算转矩反馈值TRQFB，以使驱动轮速Ws_act收敛于目标车轮速Ws_cmd，通过使用应用Δ∑调制算法的算法来调制该转矩反馈值TRQFB，来计算调制值TRQFB_DSM。由此，当

成立时，

成立，使得调制值TRQFB_DSM被计算为表示在高频率(例如5Hz以上的频率)下的反转行为的值。

并且，根据发动机转速NE和油门开度AP来计算要求转矩TRQDRV，把自适应限制转矩TRQFF_LMT和要求转矩TRQDRV中的小的一方选择为前馈转矩TRQFF，通过对该前馈转矩TRQFF加上调制值TRQFB_DSM来计算限制转矩TRQTC，并且把限制转矩TRQTC和要求转矩TRQDRV中的小的一方选择为发动机转矩TRQENG。

在该情况下，在模型修正器120中，通过使用模型修正系数KTRQFF来修正(校正)基准前馈转矩TRQFFM来计算自适应限制转矩TRQFF_LMT，该模型修正系数KTRQFF被计算成使修正转矩误差Wetrqij收敛于值0。即，由于模型修正系数KTRQFF被计算成使自适应限制转矩TRQFF_LMT迅速收敛于估计转矩TRQINST，因而使用这种模型修正系数KTRQFF计算出的自适应限制转矩TRQFF_LMT被计算为可使驱动轮速Ws_act收敛于目标车轮速Ws_cmd的最佳转矩。

根据以上结构，在驾驶者要求的要求转矩TRQDRV过大且驱动轮很有可能滑移的情况下，由于选择限制转矩TRQTC作为发动机转矩TRQENG，因而可将发动机3产生的转矩控制成以可避免驱动轮的滑移的最佳转矩为中心、表现出周期性的变动行为。由此，可被保持在

成立的状态，从而可将调制值TRQFB_DSM即限制转矩TRQTC的变动频率保持在高频带。结果，可补偿所述的抓地行为的非线性，可实现高水准的牵引控制。而且，如上所述，可避免限制转矩TRQTC被抑制为过度小的值，可将发动机3产生的转矩保持在可避免驱动轮的滑移的转矩上限值附近，从而可在实现高水准的牵引控制的同时，获得良好的加速响应性。

并且，在模型修正器120中，通过对偏差DT(＝TRQFF_LMT—TRQINST)实施限制处理来计算转矩误差Etrq，修正转矩误差Wetrq_ij是通过把转矩误差Etf、档位系数δ_j以及权重函数WNE_ij相乘来计算的。而且，计算基准局部校正值DKTRQFF_BS_ij以使这种修正转矩误差Wetrq_ij收敛于值0，通过把档位系数δ_j、权重函数WNE_ij以及基准局部校正值DKTRQFF_BS_ij相乘来计算相乘值Dktrqff_ij，并且计算4个局部校正值DKTRQFF_j作为与前进档位值NGEAR可取的4个值1～4分别对应的相乘值Dktrqff_ij的和。然后，通过对这些局部校正值DKTRQFF_j的总和加上值1来计算模型修正系数KTRQFF，并且通过使用该模型修正系数KTRQFF来修正(校正)基准前馈转矩TRQFFM来计算自适应限制转矩TRQFF_LMT。

由此，即使在图19中的基准前馈转矩TRQFFM和发动机转速NE的关系相对于实际的最佳转矩和发动机转速NE的关系在所述3个区域中的某个区域内产生局部误差的情况下，或者在两者的实际关系在3个区域之间发生变动的情况下，也能使基准前馈转矩TRQFFM和发动机转速NE的关系恰当且迅速地与实际的最佳转矩和发动机转速NE的关系一致。结果，能可靠地取得上述的作用效果。

并且，与前进档位值NGEAR的任一值对应的3个权重函数WNE_ij的各个在所述的对应区域内被设定为值1以下的正值，且在除此以外的区域内被设定为值0，而且，相邻的各2个权重函数WNE_mj、WNE_m+1j(m＝1或2)相互交叉，并且交叉部分的两者的和被设定为权重函数WNE_ij的最大值1，因而当修正基准前馈转矩TRQFFM时，可针对3个区域连续地修正，使得作为修正结果的自适应限制转矩TRQFF_LMT不会有不连续点。由此，可避免自适应限制转矩TRQFF_LMT由于基准前馈转矩TRQFFM的不连续点而暂时成为不恰当的状态的情况。

另外，第2实施方式是将本发明的控制装置应用于把驱动轮速Ws_act作为控制量并把发动机转矩TRQENG作为控制输入的控制对象的例子，然而本发明的控制装置不限于此，作为控制对象，可应用于把各种产业设备中的输入输出分别作为控制输入和控制量的系统。

并且，第2实施方式是使用发动机转速NE作为表示内燃机的运转状态的运转状态参数的例子，然而本发明的运转状态参数不限于此，只要是表示内燃机的运转状态的参数即可。例如，可以使用气缸内的压力作为运转状态参数。

Claims (8)

1.一种控制装置，其使用控制输入对控制对象中的控制量进行控制，其特征在于，该控制装置具有：

控制量检测单元，其检测所述控制量；

目标控制量设定单元，其设定成为所述控制量的目标的目标控制量；

第1输入值计算单元，其使用规定的反馈控制算法计算用于对所述控制量进行反馈控制以使其收敛于所述目标控制量的第1输入值；

调制输入值计算单元，其通过使用规定的调制算法对所述第1输入值进行调制来计算调制输入值，以使得所述第1输入值的绝对值越小则所述调制输入值的频率越高；

动作状态参数检测单元，其检测所述控制量以外的、表示所述控制对象的动作状态的动作状态参数；

误差参数计算单元，其计算误差参数，该误差参数表示所述控制输入相对于当估计所述控制量成为所述目标控制量时的所述控制输入的估计值的误差；

修正值计算单元，其使用规定的控制算法分别计算与将所述动作状态参数可变化的区域进行划分而成的多个区域对应的多个修正值，以使所述误差参数成为规定的目标值；

模型修正单元，其使用该多个修正值来修正相关关系模型，该相关关系模型表示用于对所述控制量进行前馈控制的第2输入值与所述动作状态参数之间的关系；

第2输入值计算单元，其使用该修正后的相关关系模型和所述动作状态参数来计算所述第2输入值；以及

控制输入计算单元，其根据所述第2输入值和所述调制输入值的和，将所述控制输入计算成表现出周期性的变动行为。

2.根据权利要求1所述的控制装置，其特征在于，

上述模型修正单元通过把规定的多个函数的值乘以所述误差参数和所述规定的目标值中的一方与另一方之间的偏差，来分别计算多个相乘值，并根据该多个相乘值分别计算所述多个修正值，

所述多个区域中的相邻区域相互重叠，

所述规定的多个函数分别对应于所述多个区域，所述规定的多个函数仅在该对应的区域内被设定为值0以外的值且在该对应的区域以外被设定为值0，并且在所述相互重叠的区域内，与该重叠的区域对应的所述函数的值的总和的绝对值被设定成等于所述函数中的最大值的绝对值。

3.根据权利要求1或2所述的控制装置，其特征在于，所述误差参数计算单元具有设定单元，该设定单元将所述第1输入值设定为所述误差参数，并将所述规定的目标值设定为值0。

4.根据权利要求3所述的控制装置，其特征在于，

所述控制量是排气浓度传感器的输出，该排气浓度传感器检测在内燃机的排气通路中的催化装置的下游侧的排气中的规定成分的浓度，

所述目标控制量是当估计所述催化装置的排气净化率成为规定状态时的目标输出，

所述动作状态参数是表示所述内燃机的运转状态的运转状态参数，

所述控制输入计算单元计算在所述内燃机的混合气的空燃比控制中成为该空燃比的目标的目标空燃比，作为所述控制输入。

5.根据权利要求1或2所述的控制装置，其特征在于，所述误差参数计算单元具有：

估计值计算单元，其计算所述控制输入的所述估计值；以及

误差参数设定单元，其将所述估计值与所述第2输入值之间的偏差设定为所述误差参数。

6.根据权利要求5所述的控制装置，其特征在于，

所述控制量是以内燃机为动力源的车辆的驱动轮的速度，

所述目标控制量是用于使所述车辆的所述驱动轮相对于非驱动轮的滑移量成为规定量的目标速度，

所述动作状态参数是表示所述内燃机的运转状态的运转状态参数，

所述控制输入计算单元计算在所述内燃机的输出控制中成为该输出的目标的目标输出，作为所述控制输入。

7.根据权利要求6所述的控制装置，其特征在于，该控制装置还具有：

要求输出计算单元，其计算所述车辆的驾驶者所要求的所述内燃机的输出作为要求输出；以及

选择单元，其将该要求输出和所述目标输出中的小的一方选择为所述输出控制中的所述内燃机的输出目标。

8.根据权利要求7所述的控制装置，其特征在于，所述控制输入计算单元在所述第2输入值小于所述要求输出时，将所述目标输出计算为所述第2输入值和所述调制输入值的和，并在所述第2输入值大于等于所述要求输出时，将所述目标输出计算为所述要求输出和所述调制输入值的和。

Images