Air-fuel ratio control system for internal combustion engine

Abstract

An allowable range (adaptive allowable range) for limiting a manipulated variable generated in order to converge the difference between an output from an O2 sensor disposed downstream of a catalytic converter and a target value thereof to “0” is sequentially updated depending on how the manipulated variable deviates from the allowable range.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to an air-fuel ratio control system for an internal combustion engine.

2. Description of the Prior Art

The applicant of the present application has proposed an air-fuel ratio control system for controlling the air-fuel ratio of an internal combustion engine, more accurately, the air-fuel ratio of an air-fuel mixture to be combusted by the internal combustion engine, in order to enable a three-way catalytic converter disposed in the exhaust system of the internal combustion engine. See, for example, Japanese laid-open patent publication No. 9-273439 (U.S. Pat. No. 5,845,491).

According to the disclosed technology, an exhaust gas sensor (O 2 sensor) for detecting the concentration of a certain component of an exhaust gas that has passed through a catalytic converter, e.g., the concentration of oxygen, is disposed downstream of the catalytic converter. The air-fuel ratio of the internal combustion engine is manipulated to equalize an output from the exhaust gas sensor (a detected value of the concentration of oxygen) to a predetermined target value (constant value) for thereby enabling the catalytic converter to achieve an optimum purifying capability irrespective of aging of the catalytic converter.

More specifically, a manipulated variable for manipulating the air-fuel ratio of the internal combustion engine is successively generated in given control cycles in order to converge the output of the exhaust gas sensor to its target value based on a feedback control process (specifically, a sliding mode control process) with a manipulated variable generating means implemented by a microcomputer. Based on the manipulated variable, the amount of fuel supplied to the internal combustion engine is regulated to manipulate the air-fuel ratio of the internal combustion engine. In this manner, the internal combustion engine is controlled to emit an exhaust gas for converging the output of the exhaust gas sensor to its target value, and achieving an optimum purifying capability of the catalytic converter.

According to the above publication, the manipulated variable generated by the manipulated variable generating means represents a target value for the difference between a target air-fuel ratio for the internal combustion engine or an air-fuel ratio of the internal combustion engine and a predetermined reference value (which difference is the difference between the target air-fuel ratio and the reference value). However, the manipulated variable may alternatively be a corrective variable for the amount of fuel supplied to the internal combustion engine, or a quantity that can define the air-fuel ratio of the internal combustion engine.

In the above publication, the exhaust gas sensor disposed downstream of the catalytic converter comprises an O 2 sensor. However, the exhaust gas sensor may comprise an NOx sensor, a CO sensor, an HC sensor, or another exhaust gas sensor. It is possible to enable the catalytic converter to have a desired purifying ability by controlling the air-fuel ratio of the internal combustion engine so as to converge the output of such an exhaust gas sensor to a suitable target value.

The manipulated variable generated by the manipulated variable generating means according to the feedback control process to converge the output from the exhaust gas sensor disposed downstream of the catalytic converter to the target value for enabling the catalytic converter to have a desired purifying ability suffers a large variation due to disturbance, e.g., suffers a spike-like variation. If the air-fuel ratio of the internal combustion engine is manipulated directly using the manipulated variable generated by the manipulated variable generating means, then the air-fuel ratio may be excessively be varied, tending to make unstable the operating state of the internal combustion engine.

In view of the above drawback, it has been customary to set an allowable range for the manipulated variable used to actually manipulate the air-fuel ratio of the internal combustion engine to a fixed range. After the generated manipulated variable is limited to a value in the allowable range, the air-fuel ratio of the internal combustion engine is manipulated based on the limited manipulated variable. Specifically, when the generated manipulated variable exceeds the upper limit value or lower limit value of the allowable range, the manipulated variable is forcibly limited to the upper limit value or lower limit value of the allowable range.

The limiting process is effective to avoid excessive changes of the air-fuel ratio of the internal combustion engine depending on the manipulated variable, allowing the internal combustion engine to operate stably.

However, further studies made by the inventors of the present application have revealed that the fixed allowable range for the manipulated variable results in the following shortcomings:

A range (variation range) of values of the manipulated variable generated by the manipulated variable generating means for determining an air-fuel ratio required to converge the output of the exhaust gas sensor disposed downstream of the catalytic converter to a predetermined target value, or a central value in the range is variable by various factors including operating conditions (e.g., the atmospheric pressure, intake air temperature and humidity, rotational speed, etc.) of the internal combustion engine, the nature of the fuel of the internal combustion engine, states (temperature, deterioration, etc.) of the catalytic converter.

Therefore, even when the range of values of the manipulated variable generated by the manipulated variable generating means is stable and less subject to the effect of disturbances, the range may be displaced toward the upper or lower limit of the fixed allowable range, i.e., the central value of the range of values of the generated manipulated variable may be offset from the central value of the allowable range.

In such a case, even though the manipulated variable generated by the manipulated variable generating means poses no problem in stably operating the internal combustion engine, the manipulated variable limited by the above limiting process is highly frequently limited to the upper or lower limit of the allowable range. Occasionally, the manipulated variable may be continuously held to the upper or lower limit of the allowable range.

The manipulated variable forcibly limited by the limiting process differs from the manipulated variable generated by the manipulated variable generating means in order to converge the output from the exhaust gas sensor to the target value, i.e., the manipulated variable defining the air-fuel ratio of the internal combustion engine required to converge the output from the exhaust gas sensor to the target value. Consequently, when the manipulated variable limited by the above limiting process is highly frequently limited to the upper or lower limit of the allowable range, or the manipulated variable is continuously held to the upper or lower limit of the allowable range, the controllability of the output from the exhaust gas sensor at the target value is impaired, i.e., the quick response and stability of the control process for converging the output from the exhaust gas sensor to the target value are lowered, making it difficult to achieve a desired purifying ability of the catalytic converter.

If the range of values of the manipulated variable generated by the manipulated variable generating means is displaced toward the upper or lower limit of the fixed allowable range, then when the manipulated variable contains a spike and is temporarily largely changed to one of the upper and lower limits due to disturbances, the spike manipulated variable often cannot be properly limited by the limiting process. As a consequence, the air-fuel ratio defined by the manipulated variable suffers an excessive change, impairing the stability of operation of the internal combustion engine.

SUMMARY OF THE INVENTION

It is therefore an object of the present invention to provide an air-fuel ratio control system for an internal combustion engine for generating a manipulated variable to manipulate the air-fuel ratio of the internal combustion engine in order to converge an output of an exhaust gas sensor disposed downstream of a catalytic converter to a predetermined target value, limiting the value of the manipulated variable to a predetermined allowable range, and manipulating the air-fuel ratio of the internal combustion engine based on the manipulated variable, the air-fuel ratio control system being capable of well controlling the output of the exhaust gas sensor to the target value while appropriately limiting the manipulated variable to keep operation of the internal combustion engine stable.

To achieve the above object, there is provided in accordance with the present invention an air-fuel ratio control system for an internal combustion engine, comprising an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter, manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from the exhaust gas sensor to a predetermined target value, and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range,for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable, the limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by the manipulated variable generating means.

With the above system, the allowable range is variably set depending on the manipulated variable generated by the manipulated variable generating means. The allowable range can thus be appropriately established which is suitable for the manipulated variable for manipulating the air-fuel ratio of the air-fuel mixture combusted by the internal combustion engine in order to converge the output from the exhaust gas sensor disposed in the downstream of the catalytic converter to the predetermined target value, i.e., the manipulated variable defining the air-fuel ratio of the air-fuel mixture required to converge the output from the exhaust gas sensor to the predetermined target value. As a result, the manipulated variable generated by the manipulated variable generating means (the manipulated variable may also be referred to as a demand manipulated variable in the description of the invention) is prevented from being limited more than necessary or insufficiently by the limiting means. According to the present invention, therefore, it is possible to control the output from the exhaust gas sensor at the target value while properly limiting the demand manipulated variable to achieve the stability of operation of the internal combustion engine.

The demand manipulated variable may comprise a target air-fuel ratio for the air-fuel mixture, a target value for the difference between the air-fuel ratio of the air-fuel mixture and a predetermined reference value (the difference between the target air-fuel ratio and the reference value), or a corrective variable for the amount of fuel supplied to the internal combustion engine.

If the demand manipulated variable represents the target air-fuel ratio or the difference between the target air-fuel ratio and the predetermined reference value, then the air-fuel ratio should preferably be manipulated based on the manipulated variable produced by limiting the demand manipulated variable (which manipulated variable may also be referred to as a command manipulated variable) as follows:

The actual air-fuel ratio of the air-fuel mixture is detected by a suitable sensor, and the amount of fuel supplied to the internal combustion engine is regulated according to a feedback control process in order to converge a detected value (sensor output) to the target air-fuel ratio determined by the command manipulated variable for thereby manipulating the air-fuel ratio into the target air-fuel ratio. However, the air-fuel ratio of the air-fuel mixture can be manipulated into the target air-fuel ratio according to a feed-forward control process by determining the amount of fuel supplied to the internal combustion engine from the target air-fuel ratio determined by the command manipulated variable, using a map or the like.

The limiting means comprises means for, if the value of the manipulated variable generated by the manipulated variable generating means deviates from the allowable range beyond an upper or lower limit thereof, changing the limit beyond which the manipulated variable deviates, in a direction to increase the allowable range for thereby updating the allowable range.

In a situation where a range (variation range) of values of the demand manipulated variable sequentially generated by the manipulated variable generating means for determining an air-fuel ratio required to converge the output of the exhaust gas sensor to the target value is displaced in its entirety toward one of upper and lower limits of the allowable range, the value of the demand manipulated variable tends to deviate from the allowable range beyond one of the limits even if the value of the demand manipulated variable is suitable for converging the output of the exhaust gas sensor to the target value and stabilizing operation of the internal combustion engine. Stated otherwise, the demand manipulated variable is often liable to be forcibly limited by the limiting process, i.e., the frequency at which the command manipulated variable is forcibly set to one of the limits of the allowable ranges is increased. According to the present invention, when the value of the demand manipulated variable deviates from the allowable range, the limit beyond which it deviates is changed in a direction to increase the allowable range for thereby updating the allowable range, so that the value of the demand manipulated variable will easily fall in the updated allowable range. In this manner, the value of the demand manipulated variable is prevented from being frequently forcibly limited by the limiting means, allowing the air-fuel ratio to be manipulated based on the demand manipulated variable itself, i.e., allowing the air-fuel ratio to be easily manipulated based on the command manipulated variable equal to the demand manipulated variable. As a result, the control process for converging the output from the exhaust gas sensor to the target value can well be performed while achieving its quick response and stability.

When the allowable range is updated by changing the limit thereof beyond which the demand manipulated variable deviates from the allowable range, in the direction to increase the allowable range depending on the deviation of the demand manipulated variable from the allowable range, the limit opposite to the limit beyond which the demand manipulated variable deviates may also be changed in the direction to increase the allowable range or may be maintained as it is for the purpose of well performing the control process for converging the output from the exhaust gas sensor to the target value.

However, in a situation where the range of values of the demand manipulated variable is displaced in its entirety to the limit beyond which the demand manipulated variable deviates, the range of values of the demand manipulated variable is spaced in its entirety from the limit opposite to the limit beyond which the demand manipulated variable deviates. Therefore, even when the value of the demand manipulated variable is largely changed as a temporary spike due to disturbance toward the limit opposite to the limit beyond which the demand manipulated variable deviates, the value of the demand manipulated variable is likely to fall in the allowable range. As a consequence, the demand manipulated variable tends to be insufficiently limited at the limit opposite to the limit beyond which the demand manipulated variable deviates.

In the case where the deviation of the value of the demand manipulated variable from the allowable range is temporary, if the limit of the allowable range beyond which the demand manipulated variable deviates is changed in the direction to increase the allowable range and the limit opposite to the limit beyond which the demand manipulated variable deviates is changed in the direction to increase the allowable range or is maintained as it is, then the updated allowable range becomes unnecessarily wider than a normal range of values of the demand manipulated variable. As a result, the limiting process at the upper or lower limit of the allowable range tends to become insufficient.

Therefore, in changing the limit of the allowable range beyond which the demand manipulated variable deviates in the direction to increase the allowable range depending on the deviation of the demand manipulated variable from the allowable range, the limiting means comprises means for, if the value of the manipulated variable generated by the manipulated variable generating means deviates from the allowable range beyond an upper or lower limit thereof, changing the limit opposite to the limit beyond which the manipulated variable deviates, in a direction to decrease the allowable range for thereby updating the allowable range.

Thus, when the value of the demand manipulated variable deviates from the allowable range, the limit beyond which the demand manipulated variable deviates is changed in the direction to increase the allowable range, and the limit opposite to the limit beyond which the demand manipulated variable deviates is changed in the direction to decrease the allowable range. At this time, the allowable range is updated so that the allowable range in its entirety is shifted toward the limit beyond which the demand manipulated variable deviates from the allowable range before it is updated. Therefore, in a situation where the range of values of the demand manipulated variable is displaced toward either one of the upper and lower limits of the present allowable range, if the allowable range is updated depending on the deviation of the value of the demand manipulated variable from the allowable range beyond one of the limits thereof, then the other limit of the updated allowable range is closer to the range of values of the demand manipulated variable. If the deviation of the value of the demand manipulated variable from the allowable range is temporary, then since the limit opposite to the limit beyond which the demand manipulated variable deviates is changed in the direction to decrease the allowable range, so that after the allowable range is updated, the demand manipulated variable is liable to deviate beyond the opposite limit. Depending on the deviation, the limits (both the upper and lower limits) of the allowable range are changed to update the allowable range as described above, for thereby restoring the allowable range to a range suitable for the value of the demand manipulated variable.

Therefore, when the allowable range is updated depending on the deviation of the demand manipulated variable from the allowable range by changing both the upper and lower limits thereof, it is possible to establish an allowable range that matches the range of values of the demand manipulated variable. As a consequence, the manipulated variable limited by the above limiting process is prevented from being frequently limited forcibly, and the limiting process can well be performed. Thus, the control process for converging the output of the exhaust gas sensor to the target value can be carried out well, and the stability of operation of the internal combustion engine is achieved.

The limiting means comprises means for, if the value of the manipulated variable generated by the manipulated variable generating means is present in the allowable range, changing at least one of upper and lower limits of the allowable range in a direction to decrease the allowable range for thereby updating the allowable range.

If the value of the demand manipulated variable does not deviate from the allowable range, then the allowable range is considered to be excessively wider than the range of values of the demand manipulated variable for converging the output of the exhaust gas sensor to the target value. When the demand manipulated variable temporarily deviates from the allowable range due to disturbance or the like, the allowable range may become unnecessarily wider than the range of values of the demand manipulated variable because of the process of updating the allowable range due to the deviation of the demand manipulated variable, particularly, the process of changing the limit beyond which the demand manipulated variable deviates in the direction to increase the allowable range.

According to the present invention, if the manipulated variable is present in the allowable range, then one of the upper and lower limits (preferably both the upper and lower limits) of the allowable range is changed in the direction to decrease the allowable range. Therefore, even if the allowable range becomes unnecessarily wider than the range of values of the demand manipulated variable because of the process of updating the allowable range, the allowable range can quickly be decreased so as to match the range of values of the demand manipulated variable. As a result, the limiting process is carried so as to be suitable for the rang e of values of the demand manipulated variable for reliably preventing the air-fuel ratio from being manipulated based on the demand manipulated variable which is in the form of a temporary spike.

According to the present invention, the process of updating the allowable range depending on the value of the demand manipulated variable should be carried out in a complex manner. Specifically, the limiting means comprises means for changing an upper limit of the allowable range in a direction to increase or decrease the allowable range depending on whether the value of the manipulated variable generated by the manipulated variable generating means is larger than the upper limit or not, and changing a lower limit of the allowable range in a direction to increase or decrease the allowable range depending on whether the value of the manipulated variable generated by the manipulated variable generating means is smaller than the lower limit or not, for thereby updating the allowable range.

When the demand manipulated variable deviates from the allowable range beyond its upper limit, since the value of the demand manipulated variable is greater than the upper limit and the lower limit, the upper limit is changed in the direction to increase the allowable range, and the lower limit is changed in the direction to decrease the allowable range. When the demand manipulated variable deviates from the allowable range beyond its lower limit, since the value of the demand manipulated variable is smaller than the upper limit and the lower limit, the lower limit is changed in the direction to increase the allowable range, and the upper limit is changed in the direction to decrease the allowable range. When the demand manipulated variable is present in the allowable range, since the value of the demand manipulated variable is equal to or smaller than the upper limit and equal to or greater than the lower limit, the upper and lower limits are changed in the direction to decrease the allowable range.

When the allowable range is updated by thus changing the upper and lower limits thereof, an optimum allowable range can adaptively be set up so as to exactly match the range of values of the demand manipulated variable. Consequently, the forced limitation of the demand manipulated variable according to the limiting process is held to a minimum, and the limiting process can be performed more appropriately. Thus, the control process for converging the output of the exhaust gas sensor to the target value can be carried out better, and the stability of operation of the internal combustion engine is achieved more reliably.

In updating the allowable range by changing the upper or lower limit thereof in the direction to increase or decrease the allowable range, the upper or lower limit is preferably changed by the limiting means in the direction to increase the allowable range, by an amount greater than an amount by which the upper or lower limit is changed in the direction to decrease the allowable range.

Since the amount by which the upper or lower limit is changed in the direction to increase the allowable limit is made relatively large, and the amount by which the upper or lower limit is changed in the direction to decrease the allowable range is made relatively small, the value of the demand manipulated variable can be brought quickly into the allowable range. As a result, the control process for converging the output of the exhaust gas sensor to the target value can be carried out better.

When the allowable range is updated as described above, the limiting process can be performed based on the updated allowable range. However, when the demand manipulated variable is temporarily largely changed due to disturbance or the like, in order to reliably prevent the air-fuel ratio from being manipulated based on the demand manipulated variable thus changed, the limiting means preferably comprises means for limiting the manipulated variable to the value in the allowable range before the allowable range is updated.

The limiting means comprises means for limiting the upper or lower limit of the allowable range which is being changed depending on the value of the manipulated variable generated by the manipulated variable generating means, to a value in a predetermined range corresponding to the upper or lower limit of the allowable range.

Therefore, when the allowable range is updated as described above, the upper limit or lower limit of the allowable range is prevented from being excessively large or small, and hence the allowable range for the demand manipulated variable for manipulating the air-fuel ratio of the internal combustion engine into an inappropriate air-fuel ratio for the stable operation of the internal combustion engine is prevented from being set up.

When the manipulated variable is generated by the manipulated variable generating means in given control cycles, the above process of updating the allowable range should preferably be sequentially performed each time the manipulated variable is generated, i.e., in each of the control cycles. The process of updating the allowable range may also be periodically carried out in cycles some-what longer than the control cycles.

The amount by which the upper and lower limits of the allowable range are changed may be of a predetermined fixed value, but may also be determined depending on how the demand manipulated variable deviates from the allowable range.

The manipulated variable generating means comprises means for generating the manipulated variable (demand manipulated variable) according to a sliding mode control process, though the manipulated variable generating means can generate the manipulated variable (demand manipulated variable) according to any of various feedback control processes such as a PI control process.

The sliding mode control process is a variable-structure feedback control process, and is highly stable against disturbances. Therefore, the value of the demand manipulated variable is less subject to excessive changes due to disturbances. As a result, when the allowable range is variable set as described above, it is possible to reliably establish an allowable range suitable for the range of values of the demand manipulated variable without impairing the stability of the operating states of the internal combustion engine.

One sliding mode control process is an adaptive sliding mode control process which is a combination of an ordinary sliding mode control process and a control law referred to as an adaptive control law for eliminating the effect of disturbances as much as possible. It is particularly preferable to employ such an adaptive sliding mode control process.

When the demand manipulated variable is generated according to the sliding mode control process, the limiting means comprises means for sequentially determining whether the output from the exhaust gas sensor is stable or not, and, if the output from the exhaust gas sensor is determined to be unstable, forcibly setting the allowable range to a predetermined range at least while the output from the exhaust gas sensor is continuously being determined to be unstable.

Specifically, in a situation where the output from the exhaust gas sensor is unstable, it is preferable to positively reduce the variation range of the command manipulated variable (a manipulated variable used to manipulate the air-fuel ratio) produced by limiting the demand manipulated variable in order to stabilize the output from the exhaust gas sensor. According to the present invention, at least while the output from the exhaust gas sensor is continuously being determined to be unstable, the allowable range is forcibly set to a predetermined range (which, basically, should preferably be a relatively narrow range). In this manner, the output from the exhaust gas sensor can be stabilized, and the purifying capability of the catalytic converter can be stabilized.

The limiting means comprises means for, if the output from said exhaust gas sensor is determined to be unstable, determining the level of instability in a plurality of stages, and thereby making said predetermined range to which said allowable range is set by said limiting means to be determined smaller as the level of instability determined by said limiting means is higher.

As described above, the level of instability of the output from the exhaust gas sensor is determined, and the range to be set up as the allowable range is made smaller as the level of instability is higher (the output from the exhaust gas sensor is unstabler). Therefore, the output from the exhaust gas sensor is stabilized more reliably. If the level of instability is not so high, the frequency of manipulation of the air-fuel ratio using the demand manipulated variable itself, i.e., manipulation of the air-fuel ratio based on the command manipulated variable that is equal to the demand manipulated variable, is increased as much as possible for converging the output from the exhaust gas sensor to the target value.

The limiting means comprises means for determining whether the output from the exhaust gas sensor is stable or not based on the value of a switching function used in the sliding mode control process.

The limiting means comprises means for sequentially temporarily determining whether the output from the exhaust gas sensor is stable or not based on the value of a switching function used in the sliding mode control process, and determining whether the output from the exhaust gas sensor is stable or not and the level of instability based on the frequency of the temporarily determined instability within a predetermined period of time.

In the sliding mode control process (including the adaptive sliding mode control process), a function referred to as a switching function is used, and it is important to converge the value of the switching function stably to “0” in stably converging the controlled variable (which is the output from the exhaust gas sensor in this invention) to the target value. Thus, the stability of the output from the exhaust gas sensor can be determined based on the value of the switching function.

For example, when the product of the value of the switching function and its rate of change (which product corresponds to the time-differentiated value of a Lyapunov function) is determined, if the product is of a positive value, then the value of the switching function is getting away from “0”, and if the product is of a negative value, then the value of the switching function is getting closely to “0”. Basically, therefore, it is possible to determine whether the output from the exhaust gas sensor is unstable or stable depending on whether the above product is of a positive value or a negative value. Instead, the stability of the output from the exhaust gas sensor can be determined by comparing the magnitude of the value of the switching function or the magnitude of its rate of change with a suitable given value.

In determining the level of instability, the determination of stability (whether the output from the exhaust gas sensor is stable or not) based on the value of the switching function is sequentially temporarily carried out. It is determined whether the output from the exhaust gas sensor is stable or not and the level of instability based on the frequency of the temporarily determined instability within a predetermined period of time. For example, in a situation where the frequency is equal to or smaller than a predetermined value, the output from the exhaust gas sensor is determined as being stable, and in a situation where the frequency is higher than the predetermined value, the output from the exhaust gas sensor is determined as being unstabler as the frequency is higher.

The switching function is represented by a linear function or the like having as a component thereof time-series data of the difference between the controlled variable and the target value, for example.

The manipulated variable (demand manipulated variable) generated by the manipulated variable generating means and the upper and lower limits of the allowable range represent differences with a predetermined reference value for the air-fuel ratio of the air-fuel mixture, further comprising reference value variable setting means for variably setting the reference value depending on the manipulated variable (demand manipulated variable) generated by the manipulated variable generating means.

When the demand manipulated variable and the upper and lower limits of the allowable range represent differences with a predetermined reference value for the air-fuel ratio of the air-fuel mixture, the demand manipulated variable added to or subtracted from the reference value corresponds to an air-fuel ratio (target air-fuel ratio) of the air-fuel mixture required to converge the output from the exhaust gas sensor to the target value. By variably setting the reference value depending on the demand manipulated variable, it is possible to balance the range of values of the target air-fuel ratio of the internal combustion engine which is defined depending on the value of the demand manipulated variable, above and below the reference value as a central value, which is equivalent to balancing the range of values of the demand manipulated variable between positive and negative values about a central value of “0”. As a result, the upper and lower limits of the allowable range can be balanced as positive and negative values, and the demand manipulated variable which is in the form of a temporary spike can appropriately limited by the limiting process.

If the reference value is variably set while the manipulated variable generating means is generating the manipulated variable (demand manipulated variable) according to the sliding mode control process, then the sliding mode control process preferably comprises an adaptive sliding mode control process, and the manipulated variable (demand manipulated variable) generated by the manipulated variable generating means according to the adaptive sliding mode control process includes an adaptive control law component based on an adaptive control law of the adaptive sliding mode control process, the reference value variable setting means comprising means for variably setting the reference value based on the value of the adaptive control law component of the manipulated variable.

The demand manipulated variable generated as a control input to be given to a controlled object according to the adaptive sliding mode control process is given as the sum of a component based on a control law (so-called equivalent control input) for holding the value of a switching function to “0”, a component based on a reaching control law for converging the value of the switching function to “0”, and a component based on an adaptive control law (adaptive algorithm) for eliminating disturbances as much as possible in converging the value of the switching function to “0”. According the finding of the inventors of the present application, by variably setting the reference value depending on the component (adaptive control law component) based on the adaptive control law (adaptive algorithm), the range of values of the target air-fuel ratio for the air-fuel mixture combusted in the internal combustion engine can be balanced above and below the reference value as a central value, which is equivalent to balancing the range of values of the demand manipulated variable between positive and negative values. Specifically, the component based on the adaptive control law is of a value proportional to a time integral of the value of the switching function, for example.

For variably setting the reference value depending on the value of the adaptive control law component of the demand manipulated variable, the reference value variable setting means comprises means for variably setting the reference value by increasing or decreasing the reference value depending on the magnitude of the value of the adaptive control law component of the manipulated variable with respect to a predetermined value or a range close to and containing the predetermined value.

The value of the adaptive control law component is variably set so as to be converged to the predetermined value or a value close thereto. Thus, the reference value can appropriately be variably set depending on the value of the adaptive control law component. Since the reference value serves as a reference for the demand manipulated variable generated according to the adaptive sliding mode control process, if the reference value is varied too frequently, the demand manipulated variable generated according to the adaptive sliding mode control process is adversely affected, tending to impair the stability of the output from the exhaust gas sensor. Therefore, in variably setting the reference value depending on the value of the adaptive control law component, it is preferable to increase or decrease the reference value depending on the magnitude of the value of the adaptive control law component with respect to the predetermined value or the range close to and containing the predetermined value (the reference value is not varied when the value of the adaptive control law component is present in the above close range.

For variably setting the reference value depending on the demand manipulated variable generated according to the sliding mode control process, the reference value variable setting means comprises means for sequentially determining whether the output from the exhaust gas sensor is stable or not, and holding the reference value to a predetermined value irrespective of the manipulated variable if the output from the exhaust gas sensor is unstable.

Specifically, in a situation where the output from the exhaust gas sensor is determined to be unstable, even if the reference value is variably set depending on the demand manipulated variable, the reliability of the reference value is poor. Therefore, the reference value is not variably set depending on the demand manipulated variable, but held to a predetermined value (e.g., the present value or a predetermined fixed value). Thus, the reliability of the reference value as a central value in the range of values of the target air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine can be increased. In the situation where the output from the exhaust gas sensor is determined to be unstable, the process of holding the reference value to the predetermined value is preferably carried out when the reference value is variably set depending on the value of the adaptive control law component.

As with the limiting means for determining whether the output from the exhaust gas sensor is stable, the reference value variable setting means comprises means for determining whether the output from the exhaust gas sensor is stable or not based on the value of a switching function used in the sliding mode control process.

Furthermore, the reference value variable setting means comprises means for determining whether the output from the exhaust gas sensor is substantially converged to the target value or not, and holding the reference value to a predetermined value irrespective of the manipulated variable (demand manipulated variable) if the output from the exhaust gas sensor is not converged to the target value.

Specifically, in a situation where the output from the exhaust gas sensor is not converged to the target value, the generated state of the manipulated variable tends to be unstable. Therefore, the reference value is not variably set, but is held to a predetermined value (e.g., the present value or a predetermined fixed value). In this manner, the reliability of the reference value as a central value in the range of values of the target air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine can be increased.

The process of holding the reference value to the predetermined value in the situation where the output from the exhaust gas sensor is not converged to the target value is preferably carried out when the reference value is variably set depending on the value of the adaptive control law component relative to the sliding mode control process. Whether the output from the exhaust gas sensor is substantially converged to the target value or not can be determined by comparing the magnitude of the difference of the output from the exhaust gas sensor and the target value with a suitable given value. More specifically, if the magnitude of the difference is equal to or smaller than the given value, then the output from the exhaust gas sensor is judged as being substantially converged to the target value, and if the magnitude of the difference is greater than the given value, then the output from the exhaust gas sensor is judged as being not converged to the target value.

The limiting means comprises means for recognizing an operating state of the internal combustion engine, and setting the allowable range depending on the recognized operating state.

By setting the allowable range depending on the operating state of the internal combustion engine, the process of limiting the demand manipulated variable is made suitable for the operating state of the internal combustion engine, making it possible to increase the quick response and stability of the control process for converging the output from the exhaust gas sensor to the target value in various operating states of the internal combustion engine.

More specifically, the exhaust gas sensor comprises an oxygen concentration sensor, and the operating state of the internal combustion engine recognized by the limiting means includes an elapsed time after the internal combustion engine starts to operate, and the limiting means comprises means for setting the allowable range to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease the allowable range beyond a predetermined value until the elapsed time reaches a predetermined time.

Generally, while the internal combustion engine is being started by cranking, air (oxygen) delivered into the catalytic converter is stored in the catalytic converter (so-called O 2 storage). Therefore, the output from the oxygen concentration sensor disposed downstream of the catalytic converter (which output depends on the air-fuel ratio of the air-fuel mixture which is turned into an exhaust gas upon combustion) indicates a leaner air-fuel ratio with respect to the target value. At this time, the demand manipulated variable is of such a value which makes the air-fuel ratio defined thereby richer. In such a situation, for converging the output from the oxygen concentration sensor as the exhaust gas sensor quickly to the target value, it is preferable that the air-fuel ratio of the air-fuel mixture manipulated based on the command manipulated variable produced by limiting the demand manipulated variable be varied so as to be as richer as possible. According to the present invention, therefore, immediately after the internal combustion engine starts until the elapsed time after the start of the internal combustion engine reaches a predetermined value, the allowable range is established to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease the allowable range beyond a predetermined value. Thus, the command manipulated variable produced by limiting the demand manipulated variable can take a relatively large value close to the limit of the allowable range which corresponds to the richer side of the air-fuel ratio of the air-fuel mixture, so that the air-fuel ratio tends to vary toward a richer value. As a result, the output from the oxygen concentration sensor disposed downstream of the catalytic converter can quickly be converged to the target value after the start of the internal combustion engine, and hence the desired purifying capability of the catalytic converter can quickly be achieved after the start of the internal combustion engine.

Furthermore, if the exhaust gas sensor comprises an oxygen concentration sensor, then the operating state of the internal combustion engine recognized by the limiting means includes an elapsed time after the supply of fuel to the internal combustion engine is cut off, and the limiting means comprises means for setting the allowable range to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease the allowable range beyond a predetermined value until the elapsed time reaches a predetermined time.

Generally, immediately after the supply of fuel to the internal combustion engine is cut off, the output from the oxygen concentration sensor tends to indicate a richer air-fuel ratio since the catalytic converter has stored a large amount of oxygen while the supply of fuel is being cut off. Therefore, for converging the output from the oxygen concentration sensor quickly to the target value, it is preferable that the air-fuel ratio of the air-fuel mixture manipulated based on the command manipulated variable produced by limiting the demand manipulated variable be varied so as to be as richer as possible. According to the present invention, therefore, immediately after the supply of fuel to the internal combustion engine is cut off until the elapsed time after the supply of fuel is cut off reaches a predetermined value, the allowable range is established to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease the allowable range beyond a predetermined value. Thus, as immediately after the start of the internal combustion engine, the output from the oxygen concentration sensor disposed down-stream of the catalytic converter can quickly be converged to the target value after the supply of fuel is cut off, and hence the desired purifying capability of the catalytic converter can quickly be achieved after the supply of fuel is cut off.

Moreover, if the exhaust gas sensor comprises an oxygen concentration sensor, then the operating state of the internal combustion engine recognized by the limiting means includes an elapsed time after the internal combustion engine starts driving a load, and the limiting means comprises means for setting the allowable range to inhibit at least one, which corresponds to a leaner side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to increase the allowable range beyond a predetermined value until the elapsed time reaches a predetermined time.

Specifically, immediately after the internal combustion engine starts driving a load (which corresponds to the drive wheels of a vehicle propelled by the internal combustion engine, for example), the actual air-fuel ratio tends to be leaner than the air-fuel ratio determined depending on the demand manipulated variable. In such a case, when the demand manipulated variable has a relatively large value in a direction corresponding to a leaner air-fuel ratio, if the air-fuel ratio of the air-fuel mixture is manipulated according to the demand manipulated variable directly as the command manipulated variable, the output from the oxygen concentration sensor disposed downstream of the catalytic converter is liable to be controlled excessively at a leaner air-fuel ratio than the target value. According to the present invention, therefore, until the elapsed time after the internal combustion engine start driving its load reaches a predetermined value, the allowable range is set to inhibit at least one, which corresponds to a leaner side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to increase the allowable range beyond a predetermined value. Thus, the command manipulated variable produced by limiting the demand manipulated variable is prevented from becoming too large close to the limit value corresponding to a leaner air-fuel ratio. The air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine is prevented from being excessively leaner, avoiding a situation where the output from the oxygen concentration sensor is controlled excessively at a leaner air-fuel ratio than the target value.

According to the present invention, the operating state of the internal combustion engine recognized by the limiting means includes an idling state or an operating state other than the idling state, and the limiting means comprises means for setting a predetermined range as the allowable range irrespective of the manipulated variable (demand manipulated variable) generated by the manipulated variable generating means if the operating state is the idling state, and variably setting the allowable range depending on the manipulated variable (demand manipulated variable) generated by the manipulated variable generating means if the operating state is the operating state other than the idling state.

While the internal combustion engine is in the idling state, it is not preferable for the air-fuel ratio to vary largely for stable operation of the internal combustion engine. If the allowable range is variably set depending on the value of the demand manipulated variable in the idling state, then the allowable range tends to become too large, and the variation of the air-fuel ratio manipulated based on the command manipulated variable produced by limiting the demand manipulated variable tends to be relatively large. According to the present invention, therefore, while the internal combustion engine is in the idling state, the allowable range is set to a predetermined range (relatively small range) irrespective of the demand manipulated variable, and while the internal combustion engine is in the operating state other than the idling state, the allowable range is variably set depending on the demand manipulated variable. When the internal combustion engine is in the idling state, therefore, the variation of the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine is reduced for stably operating the internal combustion engine. When the internal combustion engine is in the operating state other than the idling state, the output from the exhaust gas sensor can well be converged to the target value by variably setting the allowable range depending on the demand manipulated variable.

If the exhaust gas sensor comprises an oxygen concentration sensor, then the internal combustion engine has operation modes including a normal operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated based on the manipulated variable limited by the limiting means, and a lean operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated into a lean air-fuel ratio irrespective of the manipulated variable. The operating state of the internal combustion engine recognized by the limiting means includes an elapsed time from a transition from the lean operation mode to the normal operation mode. The limiting means comprises means for setting the allowable range to inhibit at least one, which corresponds to a leaner side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to increase the allowable range beyond a predetermined value until the elapsed time reaches a predetermined time.

Specifically, in order to reduce fuel consumption, the internal combustion engine may operate on a lean air-fuel mixture where the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine represents less fuel than the stoichiometric air-fuel ratio under suitable operating conditions. In such a situation, since the air-fuel ratio of the air-fuel mixture cannot be manipulated to converge the output from the oxygen concentration sensor as the exhaust gas sensor, the air-fuel ratio is not manipulated based on the command manipulated variable which is produced by limiting the demand manipulated variable generated by the manipulated variable generating means. According to the present invention, therefore, in the case where the internal combustion engine may operate on a lean air-fuel mixture, the internal combustion engine has operation modes including a normal operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated based on the manipulated variable (command manipulated variable) limited by the limiting means, and a lean operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated into a lean air-fuel ratio irrespective of the manipulated variable.

If the operation modes of the internal combustion engine include the normal operation mode and the lean operation mode, then the output from the oxygen concentration sensor represents a leaner air-fuel ratio when the internal combustion engine operates in the lean operation mode. For converging the output from the oxygen concentration sensor quickly to the target value when the lean operation modes changes to the normal operation mode, it is preferable that the air-fuel ratio of the air-fuel mixture manipulated based on the manipulated variable (command manipulated variable) produced by limiting the demand manipulated variable be varied so as to be as richer as possible. According to the present invention, therefore, until the elapsed time from the transition from the lean operation modes to the normal operation mode reaches the predetermined value, i.e., immediately after the lean operation modes changes to the normal operation mode, the allowable range is established to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease the allowable range beyond a predetermined value. Thus, as immediately after the start of the internal combustion engine or immediately after the supply of fuel is cut off, the output from the oxygen concentration sensor disposed downstream of the catalytic converter can quickly be converged to the target value after the end of the lean operation mode of the internal combustion engine, and hence the desired purifying capability of the catalytic converter can quickly be achieved after the end of the lean operation mode of the internal combustion engine.

With the internal combustion engine being operable on a lean air-fuel mixture, the catalytic converter preferably comprises a catalytic converter for absorbing a nitrogen oxide in the exhaust gas when the concentration of oxygen in the exhaust gas passing through the catalytic converter is a concentration corresponding to a lean air-fuel ratio of the air-fuel mixture, and reducing the nitrogen oxide in the exhaust gas when the concentration of oxygen in the exhaust gas passing through the catalytic converter is a concentration corresponding to a rich air-fuel ratio of the air-fuel mixture.

The above catalytic converter absorbs the nitrogen oxide in the exhaust gas when the operation mode of the internal combustion engine is the lean operation mode. When the operation mode of the internal combustion engine changes from the lean operation mode to the normal operation mode, the demand manipulated variable generated by the manipulated variable generating means immediately after the changing or transition of the operation mode is a manipulated variable that changes the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine largely to a richer air-fuel ratio. By setting the allowable range limited by the limiting means as described above immediately after the transition of the operation mode, the demand manipulated variable is less subject to being limited by the limiting process. The command manipulated variable produced by the limiting process is also a manipulated variable that changes the air-fuel ratio of the air-fuel mixture largely to a richer air-fuel ratio. As a result, the air-fuel ratio of the air-fuel mixture is quickly manipulated into a much richer air-fuel ratio immediately after the transition from the lean operation mode to the normal operation mode. Thus, the nitrogen oxide absorbed by the catalytic converter is quickly reduced by the catalytic converter. At the same time, since the air-fuel ratio of the air-fuel mixture is manipulated to converge the output from the oxygen concentration sensor (exhaust gas sensor) to the target value, the desired purifying capability of the catalytic converter can be achieved.

According to the present invention, if the exhaust gas sensor comprises an oxygen concentration sensor, then the target value for the output from the exhaust gas sensor be of a predetermined constant value for the catalytic converter to achieve an optimum purifying capability.

The above and other objects, features, and advantages of the present invention will become apparent from the following description when taken in conjunction with the accompanying drawings which illustrate preferred embodiments of the present invention by way of example.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of an air-fuel ratio control system for an internal combustion engine according to the present invention;

FIG. 2 is a diagram showing output characteristics of an O 2 sensor used in the air-fuel ratio control system shown in FIG. 1 ;

FIG. 3 is a block diagram showing a basic arrangement of a target air-fuel ratio generator in the air-fuel ratio control system shown in FIG. 1 ;

FIG. 4 is a diagram illustrative of a sliding mode control process employed by the air-fuel ratio control system shown in FIG. 1 ;

FIG. 5 is a block diagram of an adaptive controller in the air-fuel ratio control system shown in FIG. 1 ;

FIG. 6 is a flowchart of an engine operation control process carried out by the air-fuel ratio control system shown in FIG. 1 ;

FIG. 7 is a flowchart of a subroutine of the flowchart shown in FIG. 6 ;

FIG. 8 is a flowchart of an overall process carried out by the target air-fuel ratio generator in the air-fuel ratio control system shown in FIG. 1 ;

FIG. 9 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 10 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 11 is a diagram illustrative of partial processing of the flowchart shown in FIG. 10 ;

FIG. 12 is a diagram illustrative of partial processing of the flowchart shown in FIG. 10 ;

FIG. 13 is a flowchart of a subroutine of the flowchart shown in FIG. 10 ;

FIG. 14 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 15 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 16 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 17 is a diagram illustrative of the flowchart shown in FIG. 16 ;

FIG. 18 is a diagram illustrative of an allowable range for a limiting process carried out in a subroutine of the flowchart shown in FIG. 8 ;

FIG. 19 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ;

FIG. 20 is a flowchart of a subroutine of the flowchart shown in FIG. 19 ;

FIG. 21 is a diagram illustrative of the flowchart shown FIG. 19 ;

FIG. 22 is a diagram illustrative of the flowchart shown FIG. 19 ;

FIG. 23 is a diagram illustrative of the flowchart shown FIG. 19 ;

FIG. 24 is a diagram illustrative of the flowchart shown FIG. 19 ;

FIG. 25 is a diagram illustrative of the flowchart shown FIG. 19 ;

FIG. 26 is a flowchart of a subroutine of the flowchart shown in FIG. 8 ; and

FIG. 27 is a diagram illustrative of the flowchart shown FIG. 26 .

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

An air-fuel ratio control system for an internal combustion engine according to the present invention will be described below with reference to FIGS. 1 through 27 .

FIG. 1 shows in block form the plant control system according to the embodiment. As shown in FIG. 1 , a four-cylinder internal combustion engine 1 mounted on a propulsion source on an automobile or a hybrid vehicle, i.e., a drive source for drive wheels thereof. The internal combustion engine 1 generates exhaust gases produced by combustion of an air-fuel mixture in the cylinders. The exhaust gases are collected into a common discharge pipe 2 (exhaust passage) positioned near the internal combustion engine 1 , from which the exhaust gas is discharged into the atmosphere. Two catalytic converters 3 , 4 , each comprising a three-way catalytic converter combined with a NOx absorber (nitrogen oxide absorber), are mounted in the common exhaust pipe 2 at successively downstream locations thereon.

The NOx absorber in each of the catalytic converters 3 , 4 may be of the occlusion type which occludes NOx when the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine 1 is of a lean state (the fuel is lesser than the stoichiometric air-fuel ratio) and the concentration of oxygen in the exhaust gas is relatively high (NOx in the exhaust gas is relatively high), or of the adsorption type which adsorbs NOx in the exhaust gas when the air-fuel ratio of the air-fuel mixture is of a lean state. The occlusion-type NOx absorber may be made of barium oxide (BaO), and the adsorption-type NOx absorber may be made of sodium (Na), titanium (Ti), or strontium (Sr).

When the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine 1 is close to the stoichiometric air-fuel ratio or of a rich state (the fuel is greater than the stoichiometric air-fuel ratio), making the oxygen concentration in the exhaust gas relatively low, the occlusion-type NOx absorber discharges the occluded NOx, which is reduced by HC (hydrocarbon) and CO (carbon monoxide) in the exhaust gas. When the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine 1 is close to the stoichiometric air-fuel ratio or of a rich state (the fuel is greater than the stoichiometric air-fuel ratio), making the oxygen concentration in the exhaust gas relatively low, the adsorptiontype NOx absorber causes the adsorbed NOx to be reduced by HC (hydrocarbon) and CO (carbon monoxide) in the exhaust gas, producing a nitrogen gas that is discharged from the NOx absorber.

Of the catalytic converters 3 , 4 , the catalytic converter related to the present invention is the upstream catalytic converter 3 , and the downstream catalytic converter 4 may be dispensed with.

The plant control system serves to control an air-fuel ratio of the internal combustion engine 1 (more accurately, the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine 1) in order to enable the catalytic converter 3 to achieve optimum exhaust gas purifying performance. The plant control system comprises a wide-range air-fuel ratio sensor 5 mounted on the exhaust pipe 2 upstream of the catalytic converter 3 , or more precisely at a position where exhaust gases from the cylinders of the internal combustion engine 1 are put together, an O 2 sensor (oxygen concentration sensor) 6 mounted on the exhaust pipe 2 downstream of the catalytic converter 3 and upstream of the catalytic converter 4 , and a control unit 7 for carrying out a control process (described later on) based on detected output signals from the sensors 5 , 6 . The control unit 7 is supplied with detected output signals from the sensors 5 , 6 and also detected output signals from various other sensors for detecting operating conditions of the internal combustion engine 1 , including a engine speed sensor, an intake pressure sensor, a coolant temperature sensor, etc. The O 2 sensor 6 corresponds to an exhaust gas sensor.

The wide-range air-fuel ratio sensor 5 is in the form of an O 2 sensor, and outputs a signal having a level depending on the air-fuel ratio of an air-fuel mixture from which an exhaust gas introduced into the catalytic converter 3 is generated by fuel combustion in the internal combustion engine 1 . The air-fuel ratio is represented by the oxygen concentration of the exhaust gas introduced into the catalytic converter. The output signal from the wide-range air-fuel ratio sensor 5 (hereinafter referred to as an LAF sensor 5 ) is processed by a detecting circuit such as a linearizer (not shown) into a signal having an output signal KACT having a level which is proportional to the air-fuel ratio of the air-fuel mixture combusted in the internal combustion engine 1 in a wide range of air-fuel ratios, i.e., an output signal KACT representative of a detected value of the air-fuel ratio. The LAF sensor 5 is disclosed in detail in Japanese laid-open patent publication No. 4-369471, which corresponds to U.S. Pat. No. 5,391,282, and will not be described below.

The O 2 sensor 6 disposed downstream of the catalytic converter 3 generates an output signal VO2/OUT having a level depending on the oxygen concentration in the exhaust gas that has passed through the catalytic converter 3 , i.e., an output signal VO2/OUT representing a detected value of the oxygen concentration in the exhaust gas, as with ordinary O 2 sensors. The output signal VO2/OUT from the O 2 sensor 6 will change with high sensitivity in substantial proportion to the oxygen concentration in the exhaust gas that has passed through the catalytic converter 3 , with the air-fuel ratio recognized by the oxygen concentration in the exhaust gas being in a range Δ close to a stoichiometric air-fuel ratio, as shown in FIG. 2 .

The control unit 7 comprises a microcomputer and is roughly divided into a control unit 7 a for performing in predetermined control cycles a process for successively generating a target air-fuel ratio KCMD (which is also a target value for the output KACT of the LAF sensor 5 ) for the internal combustion engine 1 , and a control unit 7 b for performing in predetermined control cycles a process for controlling the air-fuel ratio of the internal combustion engine 1 at the target air-fuel ratio KCMD. The control unit 7 a will be referred to as an exhaust-side control unit 7 a , and the control unit 7 b will be referred to as an engine-side control unit 7 b.

The target air-fuel ratio KCMD generated by the exhaust-side control unit 7 a is basically an air-fuel ratio of the internal combustion engine 1 required to set the output of the O 2 sensor 6 (the detected value of the oxygen concentration) at a predetermined target value (constant value). The control cycles of the process performed by the exhaust-side control unit 7 a for generating the target air-fuel ratio KCMD are of a constant period (e.g., 30-100 ms) in view of a relatively long dead time of an exhaust system E (described later on) including the catalytic converter 3 and calculating loads.

The process performed by the engine-side control unit 7 b for controlling the air-fuel ratio of the internal combustion engine 1 , more specifically the process for adjusting the amount of fuel supplied to the internal combustion engine 1 , is required to be synchronous with the rotational speed of the internal combustion engine 1 . The control cycles of the process performed by the engine-side control unit 7 b are of a period in synchronism with a crankshaft angle period (so-called TDC) of the internal combustion engine 1 . Output data from various sensors including the LAF sensor 5 and the O 2 sensor 6 are also read in synchronism with the crankshaft angle period (TDC) of the internal combustion engine 1 .

The constant period of the control cycles of the exhaust-side control unit 7 a is longer than the crankshaft angle period (TDC).

The engine-side control unit 7 b has, as its main functional components, a basic fuel injection quantity calculator 8 for determining a basic fuel injection quantity Tim to be injected into the internal combustion engine 1 , a first correction coefficient calculator 9 for determining a first correction coefficient KTOTAL to correct the basic fuel injection quantity Tim, and a second correction coefficient calculator 10 for determining a second correction coefficient KCMDM to correct the basic fuel injection quantity Tim.

The basic fuel injection quantity calculator 8 determines a reference fuel injection quantity for the internal combustion engine 1 from the rotational speed NE and intake pressure PB of the internal combustion engine 1 using a predetermined map, and corrects the determined reference fuel injection quantity depending on the effective opening area of a throttle valve (not shown) of the internal combustion engine 1 , thereby calculating a basic fuel injection quantity Tim.

The first correction coefficient KTOTAL determined by the first correction coefficient calculator 9 serves to correct the basic fuel injection quantity Tim in view of an exhaust gas recirculation ratio of the internal combustion engine 1 , i.e., the proportion of an exhaust gas contained in an air-fuel mixture introduced into the internal combustion engine 1 , an amount of purged fuel supplied to the internal combustion engine 1 when a canister (not shown) is purged, a coolant temperature, an intake temperature, etc. of the internal combustion engine 1 .

The second correction coefficient KCMDM determined by the second correction coefficient calculator 10 serves to correct the basic fuel injection quantity Tim in view of the charging efficiency of an air-fuel mixture due to the cooling effect of fuel flowing into the internal combustion engine 1 depending on a target air-fuel ratio KCMD determined by the exhaust-side control unit 7 a (described later on).

The control unit 7 corrects the basic fuel injection quantity Tim with the first correction coefficient KTOTAL and the second correction coefficient KCMDM by multiplying the basic fuel injection quantity Tim by the first correction coefficient KTOTAL and the second correction coefficient KCMDM, thus producing a demand fuel injection quantity Tcyl for the internal combustion engine 1 .

Specific details of processes for calculating the basic fuel injection quantity Tim, the first correction coefficient KTOTAL, and the second correction coefficient KCMDM are disclosed in detail in Japanese laid-open patent publication No. 5-79374 and U.S. Pat. No. 5,253,630, and will not be described below.

The engine-side control unit 7 b also has, in addition to the above functional components, a feedback controller 14 for adjusting a fuel injection quantity of the internal combustion engine 1 according to a feedback control process so as to converge the output signal KACT of the LAF sensor 5 (the detected air-fuel ratio) to the target air-fuel ratio KCMD which is sequentially determined by the exhaust control unit 7 a.

The feedback controller 14 comprises a general feedback controller 15 for feedback-controlling a total air-fuel ratio for all the cylinders of the internal combustion engine 1 and a local feedback controller 16 for feedback-controlling an air-fuel ratio for each of the cylinders of the internal combustion engine 1 .

The general feedback controller 15 sequentially determines a feedback correction coefficient KFB to correct the demand fuel injection quantity Tcyl (by multiplying the demand fuel injection quantity Tcyl) so as to converge the output signal KACT from the LAF sensor 5 to the target air-fuel ratio KCMD.

The general feedback controller 15 comprises a PID controller 17 for determining a feedback manipulated variable KLAE as the feedback correction coefficient KFB depending on the difference between the output signal KACT from the LAF sensor 5 and the target air-fuel ratio KCMD according to a known PID control process, and an adaptive controller 18 (indicated by “STR” in FIG. 1 ) for adaptively determining a feedback manipulated variable KSTR for determining the feedback correction coefficient KFB in view of changes in operating conditions of the internal combustion engine 1 or characteristic changes thereof from the output signal KACT from the LAF sensor 5 and the target air-fuel ratio KCMD.

In present embodiment, the feedback manipulated variable KLAF generated by the PID controller 17 is of “1” and can be used directly as the feedback correction coefficient KFB when the output signal KACT (the detected air-fuel ratio) from the LAF sensor 5 is equal to the target air-fuel ratio KCMD. The feedback manipulated variable KSTR generated by the adaptive controller 18 becomes the target air-fuel ratio KCMD when the output signal KACT from the LAF sensor 5 is equal to the target air-fuel ratio KCMD. A feedback manipulated variable kstr (=KSTR/KCMD) which is produced by dividing the feedback manipulated variable KSTR by the target air-fuel ratio KCMD with a divider 19 can be used as the feedback correction coefficient KFB.

The feedback manipulated variable KLAF generated by the PID controller 17 and the feedback manipulated variable kstr which is produced by dividing the feedback manipulated variable KSTR from the adaptive controller 18 by the target air-fuel ratio KCMD are selected one at a time by a switcher 20 . A selected one of the feedback manipulated variable KLAF and the feedback manipulated variable KSTR is used as the feedback correction coefficient KFB. The demand fuel injection quantity Tcyl is corrected by being multiplied by the feedback correction coefficient KFB. Details of the general feedback controller 15 (particularly, the adaptive controller 18 ) will be described later on.

The local feedback controller 16 comprises an observer 21 for estimating real air-fuel ratios #nA/F (n=1, 2, 3, 4) of the respective cylinders from the output signal KACT from the LAF sensor 5 , and a plurality of PID controllers 22 (as many as the number of the cylinders) for determining respective feedback correction coefficients #nKLAF for fuel injection quantities for the cylinders from the respective real air-fuel ratios #nA/F estimated by the observer 21 according to a PID control process so as to eliminate variations of the air-fuel ratios of the cylinders.

Briefly stated, the observer 21 estimates a real air-fuel ratio #nA/F of each of the cylinders as follows: A system from the internal combustion engine 1 to the LAF sensor 5 (where the exhaust gases from the cylinders are combined) is considered to be a system for generating an air-fuel ratio detected by the LAF sensor 5 from a real air-fuel ratio #nA/F of each of the cylinders, and is modeled in view of a detection response delay (e.g., a time lag of first order) of the LAF sensor 5 and a chronological contribution of the air-fuel ratio of each of the cylinders to the air-fuel ratio detected by the LAF sensor 5 . Based on the modeled system, a real air-fuel ratio #nA/F of each of the cylinders is estimated from the output signal KACT from the LAF sensor 5 .

Details of the observer 21 are disclosed in Japanese laid-open patent publication No. 7-83094 and U.S. Pat. No. 5,531,208, for example, and will not be described below.

Each of the PID controllers 22 of the local feedback controller 16 divides the output signal KACT from the LAF sensor 5 by an average value of the feedback correction coefficients #nKLAF determined for all the cylinders by the respective PID controllers 22 in a preceding control cycle to produce a quotient value, and uses the quotient value as a target air-fuel ratio for the corresponding cylinder. Each of the PID controllers 22 then determines a feedback correction coefficient #nKLAF for each cylinder in a present control cycle so as to eliminate any difference between the target air-fuel ratio and the corresponding real air-fuel ratio #nA/F determined by the observer 21 .

The local feedback controller 16 multiplies a value, which has been produced by multiplying the demand fuel injection quantity Tcyl by the selected feedback correction coefficient KFB produced by the general feedback controller 15 , by the feedback correction coefficient #nKLAF for each of the cylinders, thereby determining an output fuel injection quantity #nTout (n=1, 2, 3, 4) for each of the cylinders.

The output fuel injection quantity #nTout thus determined for each of the cylinders is corrected for accumulated fuel particles on intake pipe walls of the internal combustion engine 1 by a corresponding fuel accumulation corrector 23 in the engine-side control unit 7 b . The corrected output fuel injection quantity #nTout is applied to each of fuel injectors (not shown) of the internal combustion engine 1 , which injects fuel into each of the cylinders with the corrected output fuel injection quantity #nTout.

The correction of the output fuel injection quantity in view of accumulated fuel particles on intake pipe walls is disclosed in detail in Japanese laid-open patent publication No. 8-21273 and U.S. Pat. No. 5,568,799, for example, and will not be described in detail below.

A sensor output selector 24 shown in FIG. 1 serves to select the output signal KACT from the LAF sensor 5 , which is suitable for the estimation of a real air-fuel ratio #nA/F of each cylinder with the observer 21 , depending on the operating conditions of the internal combustion engine 1 . Details of the sensor output selector 24 are disclosed in detail in Japanese laid-open patent publication No. 7-259588 and U.S. Pat. No. 5,540,209, and will not be described in detail below.

The exhaust-side control unit 7 a has, as its main functional component, a target air-fuel ratio generator 28 for sequentially generating the target air-fuel ratio KCMD in control cycles of the engine-side control unit 7 b using the data of the output signal KACT of the LAF sensor 5 which is given via the engine-side control unit 7 a and the data of the output signal VO2/OUT of the O 2 sensor 6 .

As shown in FIG. 3 , the target air-fuel ratio generator 28 comprises a reference value setting unit 11 (reference value variable setting means) for sequentially variably setting a reference value FLAF/BASE with respect to the air-fuel ratio of the internal combustion engine 1 (hereinafter referred to as an “air-fuel ratio reference value FLAF/BASE”, which is also a reference value with respect to the output signal KACT of the LAF sensor 5 ), a subtractor 12 for determining a difference kact (=KACT-FLAF/BASE) between the output signal KACT from the LAF sensor 5 and the air-fuel ratio reference value FLAF/BASE, and a subtractor 13 for determining a difference VO 2 (=VO2/OUT−VO2/TARGET) between the output signal VO2/OUT from the O 2 sensor 6 and a target value VO2/TARGET therefor. The target value VO2/TARGET for the output signal VO2/OUT from the O 2 sensor 6 is established as a constant value (fixed value) since in the present embodiment, optimum exhaust gas purifying performance of the catalytic converter 3 can be achieved with the output signal VO2/OUT from the O 2 sensor 6 being set at a certain constant value (see FIG. 2 ).

The differences kact, VO2 determined respectively by the subtractors 12 , 13 are referred to as a differential output kact of the LAF sensor 5 and a differential output VO2 of the O 2 sensor 6 , respectively.

The target air-fuel ratio generator 28 also has a manipulated variable generator 29 for sequentially generating a manipulated variable usl for controlling the air-fuel ratio of the internal combustion engine 1 in order to converge the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, i.e., to converge the differential output VO2 to “0”, using the data of the differential output kact of the LAF sensor 5 and the differential output VO2 of the O 2 sensor 6 , a limiter 30 for generating a manipulated variable kcmd by limiting the manipulated variable usl to a value within a predetermined allowable range, and an adder 31 for adding the air-fuel ratio reference value FLAF/BASE established by the reference value setting unit 11 to the manipulated variable kcmd generated by the limiter 30 .

The manipulated variable generator 29 and the limiter 30 correspond respectively to a manipulated variable generating means and a limiting means.

The manipulated variable usl generated by the manipulated variable generator 29 as described later on represents the difference between the actual air-fuel ratio (detected by the LAF sensor 5 ) of the internal combustion engine 1 required to converge the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, and the air-fuel ratio reference value FLAF/BASE. Therefore, basically, in order to converge the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, the sum of the manipulated variable usl (hereinafter referred to as a “demand differential air-fuel ratio usl”) and the air-fuel ratio reference value FLAF/BASE may be generated as the target air-fuel ratio KCMD.

The demand differential air-fuel ratio usl generated by the manipulated variable generator 29 may occasionally suffer relatively large changes due to disturbances. When the actual air-fuel ratio (detected by the LAF sensor 5 ) of the internal combustion engine 1 is controlled at the target air-fuel ratio (=usl+FLAF/BASE) determined depending on the demand differential air-fuel ratio usl, the internal combustion engine 1 may possibly operate unstably.

In order to avoid the above drawback, the limiter 30 in the target air-fuel ratio generator 28 limits the demand differential air-fuel ratio usl to generate the manipulated variable kcmd whose value is limited within the allowable range.

The limiter 30 limits the demand differential air-fuel ratio usl as follows: When the value of the demand differential air-fuel ratio usl is in the allowable range (a normal state), the limiter 30 sets up the demand differential air-fuel ratio usl as the manipulated variable kcmd. When the value of the demand differential air-fuel ratio usl deviates from the allowable range beyond its upper limit value or lower limit value, the limiter 30 forcibly sets the demand differential air-fuel ratio usl to the upper limit value or lower limit value of the allowable range.

The adder 31 in the target air-fuel ratio generator 28 adds the air-fuel ratio reference value FLAF/BASE to the manipulated variable kcmd (hereinafter referred to as a “command differential air-fuel ratio kcmd”) produced by limiting the demand differential air-fuel ratio usl, generating the target air-fuel ratio KCMD, i.e., the target air-fuel ratio KCMD (=kcmd+FLAF/BASE) to be given as a command value for the air-fuel ratio of the internal combustion engine 1 to the engine-side control unit 7 b.

The manipulated variable generator 29 will further be described below. The manipulated variable generator 29 sequentially determines the demand differential air-fuel ratio usl in control cycles of the exhaust-side control unit 7 a as a control input go be given to an exhaust system E (hereinafter referred to as an “object exhaust system E”), which extends from the position of the LAF sensor 5 in the exhaust pipe 2 (upstream of the catalytic converter 3 ) to the position of the O 2 sensor 6 (downstream of the catalytic converter 3 ) and which includes the catalytic converter 3 , in order to converge the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET therefor, i.e., to converge the differential output VO2 of the O 2 sensor 6 to “0”, according to a sliding mode control process, which is a type of feedback control process, in view of a dead time present in the object exhaust system E, dead times present in the internal combustion engine 1 and the engine control unit 7 b , and behavioral changes of the object exhaust system E.

For generating the demand differential air-fuel ratio usl, the object exhaust system E is regarded as a system for generating the differential output VO2 of the O 2 sensor 6 with a dead time and a response delay from the differential output kact of the LAF sensor 5 which corresponds to the difference between the actual air-fuel ratio of the air-fuel mixture, i.e., the air-fuel mixture from which the exhaust gas entering the catalytic converter 3 is generated upon combustion, and the air-fuel ratio reference value FLAF/BASE, and the behavior of the system is modeled in advance. In addition, the system comprising the internal combustion engine 1 and the engine-side control unit 7 b is regarded as a system (hereinafter referred to as an “air-fuel ratio manipulating system”) for generating the differential output kact of the LAF sensor 5 with a dead time from the command differential air-fuel ratio kcmd (which usually agrees with the demand differential air-fuel ratio usl) that represents the difference (=KCMD-FLAF/BASE) between the target air-fuel ratio KCMD and the air-fuel ratio reference value FLAF/BASE, i.e., a system in which the differential output kact of the LAF sensor 5 at each point of time agrees with the command differential air-fuel ratio kcmd prior to the dead time of the air-fuel ratio manipulating system, and the behavior of the system is modeled in advance.

The model representing the behavior of the object exhaust system E (hereinafter referred to as an “exhaust system model”) is expressed as a model of a discretetime system (more specifically, an autoregressive model having a dead time in the differential output kact as an input to the object exhaust system E) according to the following equation (1):

VO2(k+1)=a1·VO2(k)+a2·VO2(k−1)+b1·kact(k−d1)  (1)

where “k” represents the number of a discrete-time control cycle of the exhaust-side control unit 7 a , and “d1” the dead time present in the object exhaust system E as expressed in terms of the number of control cycles. The dead time of the object exhaust system E (the time required until the air-fuel ratio detected by the LAF sensor 5 at each point of time is reflected in the output signal VO2/OUT from the O 2 sensor 6 ) is generally equal to the time of 3-10 control cycles (d1=3-10) if the period (constant) of control cycles of the exhaust-side control unit 7 a ranges from 30 to 100 ms. In present embodiment, a preset constant value (e.g., d1=7 in the present embodiment) equal to or slightly longer than the actual dead time of the object exhaust system E is used as the dead time d1 in the exhaust system model as represented by the equation (1).

The first and second terms of the right side of the equation (1) correspond to a response delay of the object exhaust system E, the first term being a primary autoregressive term and the second term being a secondary autoregressive term. In the first and second terms, “a1”, “a2” represent respective gain coefficients of the primary autoregressive term and the secondary autoregressive term. Stated otherwise, these gain coefficients “a1”, “a2” are coefficients relative to the differential output VO2 of the O 2 sensor 6 in the exhaust system model.

The third term of the right side of the equation (1) represents the differential output kact of the LAF sensor 5 as an input to the object exhaust system E, including the dead time d1 of the object exhaust system E. In the third term, “b1” represents a gain coefficient relative to its input (=the differential output kact of the LAF sensor 5 ). The gain coefficients “a1”, “a2”, “b1” are parameters which are to be set (identified) to values in defining the behavior of the exhaust system model, and are sequentially identified by an identifier which will be described later on.

In the exhaust system model expressed as the discrete time system according to the equation (1), the differential output VO2(k+1) of the O 2 sensor 6 as an output of the object exhaust system E in each control cycle of the exhaust-side control unit 7 a is expressed by a plurality of (two in this embodiment) differential outputs VO2(k), VO(k−1) and a differential output kact(k−d1) of the LAF sensor 5 in past control cycles prior to the control cycle.

The model of the air-fuel ratio manipulating system comprising the internal combustion engine 1 and the engine-side control unit 7 b (hereinafter referred to as an “air-fuel ratio manipulating system model) is expressed as a discrete-time system model according to the following equation (2):

kact(k)=kcmd(k−d2)  (2)

where “d2” represents the dead time (second dead time) of the air-fuel ratio manipulating system in terms of the number of control cycles of the exhaust-side control unit 7 a . The dead time of the air-fuel ratio manipulating system (the time required until the target air-fuel ratio KCMD or the command differential air-fuel ratio kcmd at each point of time is reflected in the output signal KACT or the differential output kact of the LAF sensor 5 ) varies with the rotational speed of the internal combustion engine 1 , and is longer as the rotational speed of the internal combustion engine 1 is lower. In present embodiment, in view of the above characteristics of the dead time of the air-fuel ratio manipulating system expressed by the equation (2), a preset constant value (d2=3 in the present embodiment) which is equal to or slightly longer than the dead time of the actual air-fuel ratio manipulating system at an idling rotational speed of the internal combustion engine 1 (the dead time is a maximum dead time which can be taken by the air-fuel ratio manipulating system at an arbitrary rotational speed of the internal combustion engine 1 ) is used as the value of the dead time d2 in the air-fuel ratio manipulating system model expressed by the equation (2).

In the air-fuel ratio manipulating system model expressed by the equation (2), the differential output kact(k) of the LAF sensor 5 in each control cycle of the exhaust-side control unit 5 a is assumed to agree with the command differential air-fuel ratio kcmd(k−d2) prior to the dead time d2 of the air-fuel ratio manipulating system, and is expressed by the command differential air-fuel ratio kcmd(k−d2).

The air-fuel ratio manipulating system actually includes a response delay of the internal combustion engine 1 in addition to the dead time. Since a response delay of the output KACT or the differential output kact of the LAF sensor 5 with respect to the target air-fuel ratio KCMD or the command differential air-fuel ratio kcmd is basically compensated for by the feedback controller 14 (particularly the adaptive controller 18 ) in the engine-side control unit 7 b , there will arise no problem if the response delay of the internal combustion engine 1 is not taken into account in the air-fuel ratio manipulating system as viewed from the manipulated variable generator 29 in the exhaust-side control unit 7 a.

The manipulated variable generator 29 carries out the process constructed on the basis of the exhaust system model and the air-fuel ratio manipulating system model expressed by the respective equations (1), (2) in the control cycles of the exhaust-side control unit 7 a , to sequentially generate the demand manipulated quantity usl as an input to be given to the object exhaust system E for converging the output VO2/OUT of the O 2 sensor 6 to its target value VO2/TARGET. In order to generate the demand manipulated quantity usl, the manipulated variable generator 29 has its functional components as shown in FIG. 3 .

Specifically, the manipulated variable generator 29 comprises an identifier 25 for sequentially identifying in each control cycle values of the gain coefficients a1, a2, b1 that are parameters to be established for the exhaust system model (the equation (1)), using the data of the differential output kact from the LAF sensor 5 and the differential output VO2 from the O 2 sensor 6 , an estimator 26 for sequentially estimating in each control cycle an estimated value VO2 bar of the differential output VO2 from the O 2 sensor 6 (hereinafter referred to as an “estimated differential output VO2 bar”) after the total dead time d (=d1+d2) which is the sum of the dead time d1 of the object exhaust system E and the dead time d2 of the air-fuel ratio manipulating system, using the data of the differential output kact from the LAF sensor 5 , the data of the differential output VO2 from the O 2 sensor 6 , the data of the command differential air-fuel ratio kcmd (normally, kcmd=usl) produced when the demand differential air-fuel ratio usl determined in the past by a sliding mode controller 27 is limited by the limiter 30 , identified values a1 hat, a2 hat, b1 hat of the gain coefficients a1, a2, b1 that are calculated by the identifier 25 (hereinafter referred to as “identified gain coefficients a1 hat, a2 hat, b1 hat”), and a sliding mode controller 27 for sequentially determining in each control cycle the demand differential air-fuel ratio usl, using the data of the estimated differential output VO2 bar from the O 2 sensor 6 which has been determined by the estimator 26 and the identified gain coefficients a1 hat, a2 hat, b1 hat, according to an adaptive slide mode control process.

The algorithm of a processing operation to be carried out by the identifier 25 , the estimator 26 , and the sliding mode controller 27 is constructed as follows:

The identifier 25 serves to identify the gain coefficients a1, a2, b1 sequentially on a real-time basis for the purpose of minimizing a modeling error of the actual object exhaust system E of the exhaust system model expressed by the equation (1). The identifier 25 carries out its identifying process as follows:

In each control cycle of the exhaust-side control unit 7 a , the identifier 25 determines an identified value VO2(k) hat of the differential output VO2 from the O 2 sensor 6 (hereinafter referred to as an “identified differential output VO2 hat”) on the exhaust system model, using the identified gain coefficients a1 hat, a2 hat, b1 hat of the presently established exhaust system model, i.e., identified gain coefficients a1 hat (k−1), a2 hat (k−1), b1 hat (k−1) determined in a preceding control cycle, and past data of the differential output kact from the LAF sensor 5 and the differential output VO2 from the O 2 sensor 6 , according to the following equation (3):

VÔ2(k)=â1(k−1)·VO2(k−1)+â2(k−1)·VO2(k−2)+{circumflex over (b)}1(k−1)·kact(k−d1−1)  (3)

The equation (3) corresponds to the equation (1) which is shifted into the past by one control cycle with the gain coefficients a1, a2, b1 being replaced with the respective identified gain coefficients a1 hat (k−1), a2 hat (k−1), b1 hat (k−1). The value of the dead time “d1” of the object exhaust system E in the third term of the equation (3) represents a preset constant value (d1=7 in this embodiment) as described above.

If vectors Θ, ξ defined by the following equations (4), (5) are introduced (the letter T in the equations (4), (5) and other equations represents a transposition), then the equation (3) is expressed by the equation (6):

Θ T (k)=[â1(k)â2(k){circumflex over (b)}1(k)]  (4)

ξ T (k)=[VO2(k−1)VO2(k−2)kact(k−d1−1)]  (5)

VÔ2(k)=Θ T (k−1)·ξ(k)  (6)

The identifier 25 also determines a difference id/e between the identified differential output VO2 hat from the O 2 sensor 6 which is determined by the equation (3) or (6) and the present differential output VO2 from the O 2 sensor 6 , as representing a modeling error of the exhaust system model with respect to the actual object exhaust system E (hereinafter the difference id/e will be referred to as an “identified error id/e”), according to the following equation (7):

id/e(k)=VO2(k)−VÔ2(k)  (7)

The identifier 25 further determines new identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat, stated otherwise, a new vector Θ(k) having these identified gain coefficients as elements (hereinafter the new vector Θ(k) will be referred to as an “identified gain coefficient vector Θ”), in order to minimize the identified error id/e, according to the equation (8) given below. That is, the identifier 25 varies the identified gain coefficients al hat (k−1), a2 hat (k−1), b1 hat (k−1) determined in the preceding control cycle by a quantity proportional to the identified error id/e for thereby determining the new identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat.

Θ(k)=Θ(k−1)+Kθ(k)·id/e(k)  (8)

where Kθ represents a cubic vector determined by the following equation (9), i.e., a gain coefficient vector for determining a change depending on the identified error id/e of the identified gain coefficients a1 hat, a2 hat, b1 hat): $$ K ⁢   ⁢ θ ⁡ ( k ) = P ⁡ ( k - 1 ) ⁢ ξ ⁡ ( k ) 1 + ξ T ⁡ ( k ) ⁢ P ⁡ ( k - 1 ) ⁢ ξ ⁡ ( k ) ( 9 )

where P represents a cubic square matrix determined by a recursive formula expressed by the following equation (10): $$ P ⁡ ( k ) = 1 λ 1 ⁡ ( k ) ⁢   [ I - λ 2 ⁡ ( k ) ⁢ P ⁡ ( k - 1 ) ⁢ ξ ⁡ ( k ) ⁢ ξ T ⁡ ( k ) λ 1 ⁡ ( k ) + λ 2 ⁡ ( k ) ⁢ ξ T ⁡ ( k ) ⁢ P ⁡ ( k - 1 ) ⁢ ξ ⁡ ( k ) ] ⁢ P ⁡ ( k - 1 ) ( 10 )

where I represents a unit matrix.

In the equation (10), λ 1 , λ 2 are established to satisfy the conditions 0<λ 1 ≦1 and 0≦λ 2 <2, and an initial value P(0) of P represents a diagonal matrix whose diagonal components are positive numbers.

Depending on how λ 1 , λ 2 in the equation (10) are established, any one of various specific algorithms including a fixed gain method, a degressive gain method, a method of weighted least squares, a method of least squares, a fixed tracing method, etc. may be employed. According to present embodiment, a method of least squares (λ 1 =λ 2 =1), for example, is employed.

Basically, the identifier 25 sequentially determines in each control cycle the identified gain coefficients a1 hat, a2 hat, b1 hat of the exhaust system model in order to minimize the identified error id/e according to the above algorithm (calculating operation). Through this operation, it is possible to sequentially obtain the identified gain coefficients a1 hat, a2 hat, b1 hat which match the actual object exhaust system E.

The calculating operation described above is the basic processing that is carried out by the identifier 25 . In present embodiment, the identifier 25 performs additional processes such as a limiting process, on the identified gain coefficients a1 hat, a2 hat, b1 hat in order to determine them. Such operations of the identifier 25 will be described later on.

The estimator 26 sequentially determines in each control cycle the estimated differential output VO2 bar which is an estimated value of the differential output VO2 from the O 2 sensor 6 after the total dead time d (=d1+d2) in order to compensate for the effect of the dead time d1 of the object exhaust system E and the effect of the dead time d2 of the air-fuel ratio manipulating system for the calculation of the demand differential air-fuel ratio usl with the sliding mode controller 27 as described in detail later on. Since details of the estimator 26 are disclosed in U.S. patent application Ser. No. 09/311353, the estimator 26 will briefly be described below.

If the equation (2) expressing the model of the air-fuel ratio manipulating system is applied to the equation (1) expressing the exhaust system model, then the equation (1) can be rewritten as the following equation (11):

$$ VO2 ⁡ ( k + 2 ) =   ⁢ a1 · VO2 ⁡ ( k ) + a2 · VO2 ⁡ ( k - 1 ) +   ⁢ b1 · kcmd ⁡ ( k - d1 - d2 ) =   ⁢ a1 · VO2 ⁡ ( k ) + a2 · VO2 ⁡ ( k - 1 ) +   ⁢ b1 · kcmd ⁡ ( k - d ) ( 11 )

The equation (11) expresses a system which is a combination of the object exhaust system E and the air-fuel ratio manipulating system as the model of a discrete time system, regarding such a system as a system for generating the differential output VO2 of the O 2 sensor 6 from the command differential air-fuel ratio kcmd with dead times of the object exhaust system E and the air-fuel ratio manipulating system and a response delay of the object exhaust system E.

Using the equation (11), the estimated differential output VO2(k+d) bar which is an estimated value of the differential output VO2 of the O 2 sensor 6 after the total dead time d in each control cycle is expressed using present and past time-series data VO2(k), VO2(k−1) of the differential output VO2 of the O 2 sensor 6 and past time-series data kcmd(k−j) (j=1, 2, . . . , d) of the command differential air-fuel ratio kcmd (=KCMD−FLAF/BASE), according to the following equation (12): $$ VO2 _ ⁡ ( k + d ) =   ⁢ α ⁢ 1 · VO2 ( k ) + α ⁢ 2 · VO2 ( k - 1 ) +   ⁢ ∑ j = 1 d ⁢   ⁢ β j · kcmd ⁡ ( k - j ) ( 12 )

where

α1=the first-row, first-column element of A d ,

α2=the first-row, second-column element of A d ,

βj=the first-row elements of A j−1 ·B $$ A = [ a1 a2 1 0 ] B = [ b1 0 ]

In the equation (12), a1, a2 represent the first-row, first-column element and the first-row, second-column element of the dth power A d (d: total dead time) of the matrix A defined as described above, and βj represents the first-row elements of the product A j−1 ·B of the (j−1)th power A j−1 (j=1, 2, . . . , d) of the matrix A and the vector B defined as described above.

Of the past time-series data kcmd(k−j) (j=1, 2, . . . , d) of the command differential air-fuel ratio kcmd in the equation (12), the past time-series data kcmd(k−d2), kcmd(k−d2−1), . . . , kcmd(k−d) of the command differential air-fuel ratio kcmd before the dead time d2 of the air-fuel ratio manipulating system from the present time can be replaced with data kact(k), kact(k−1), . . . , kact(k−d+2), respectively, prior to the present time, of the differential output kact of the LAF sensor 5 . The replacement of the data results in the following equation (13): $$ VO2 _ ⁡ ( k + d ) =   ⁢ α ⁢ 1 · VO2 ( k ) + α ⁢ 2 · VO2 ( k - 1 ) +   ⁢ ∑ j = 1 d 2 - 1 ⁢   ⁢ β j · kcmd ⁡ ( k - j ) + ∑ i = 0 d - d 2 ⁢   ⁢ β i + d2 · kact ⁡ ( k - i ) =   ⁢ α ⁢ 1 · VO2 ( k ) + α ⁢ 2 · VO2 ( k - 1 ) +   ⁢ ∑ j = 1 d 2 - 1 ⁢   ⁢ β j · kcmd ⁡ ( k - j ) + ∑ i = 0 d 1 ⁢   ⁢ β i + d2 · kact ⁡ ( k - i ) ( 13 )

The above equation (13) is a basic equation for the estimator 26 to calculate the estimated differential output VO2(k+d) bar in this embodiment. Stated otherwise, the estimator 26 determines the estimated differential output VO2 bar of the O 2 sensor 6 according to the equation (13), using the present and past time-series data VO2(k), VO2(k−1) of the differential output VO2 of the O 2 sensor 6 , the past data kcmd(k−j) (j=1, . . . , d2−1) of the command differential air-fuel ratio kcmd which is produced by limiting the demand differential air-fuel ratio usl generated by the sliding mode controller 27 , and the present and past time-series data kact(k−i) (i=0, . . . , d1) of the differential output kact of the LAF sensor 5 .

In this embodiment, the values of the coefficients α1, α2, βj required to calculate the estimated differential output VO2 (k+d) bar according to the equation (13) are basically calculated using the identified gain coefficients a1 hat, a2 hat, b1 hat which are identified values of the gain coefficients a1, a2, b1 (these are elements of the matrix A and the vector B defined above with respect to the equation (12)). For the values of the dead times d1, d2 required for the calculation of the equation (13), the preset values described above are used.

The estimated differential output VO2(k+d) bar may be determined according to the equation (12) without using the data of the differential output kact of the LAF sensor 5 . In such a case, the estimated differential output VO2(k+d) bar is determined using the present and past time-series data VO2(k), VO2(k−1) of the differential output VO2 of the O 2 sensor 6 , the past time-series data kcmd(k−j) (j=1, 2, . . . , d) of the command differential air-fuel ratio kcmd, and the values of the coefficients α1, α2, βj (j=1, 2, . . . , d) determined by the identified gain coefficients a1 hat, a2 hat, b1 hat. For increasing the reliability of the estimated differential output VO2(k+d) bar, it is preferable to calculate the estimated differential output VO2(k+d) bar according to the equation (13) using the data of the differential output kact of the LAF sensor 5 which reflects the actual behavior of the internal combustion engine 1 .

If the dead time d2 of the air-fuel ratio manipulating system may be set to d2=1, then all the past time-series data kcmd(k−j) (j=1, 2, . . . , d) of the command differential air-fuel ratio kcmd in the equation (12) may be replaced with the time-series data kact(k), kact(k−1), . . . , kact(k−d+d2), prior to the present time, of the LAF sensor 5 . In this case, the estimated differential output VO2(k+d) bar can be determined according to the following equation (14) which includes no data of the command differential air-fuel ratio kcmd: $$ VO2 _ ⁡ ( k + d ) =   ⁢ α ⁢ 1 · VO2 ⁡ ( k ) + α ⁢ 2 · VO2 ⁢ ( k - 1 ) +   ⁢ ∑ j = 0 d - 1 ⁢   ⁢ β j + 1 · kact ⁡ ( k - j ) ( 14 )

In the present embodiment, the estimator 26 determines in each control cycle the estimated differential output VO2 bar of the O 2 sensor according to the equation (13) which differs from the equation (12) in that all, prior to the dead time d2 of the air-fuel ratio manipulating system, of the time-series data of the target differential air-fuel ratio kcmd in the equation (12) are replaced with the differential output kact of the LAF sensor 5 . However, the estimated differential output VO2 bar may be determined according to an equation which differs from the equation (12) in that only part of the time-series data of the target differential air-fuel ratio kcmd prior to the dead time d2 in the equation (12) is replaced with the differential output kact of the LAF sensor 5 .

The above calculating process is the basic algorithm for the estimator 26 to determine the estimated differential output VO2(k+d) bar which is an estimated value after the total dead time d of the differential output VO2 of the O 2 sensor in each control cycle.

The sliding mode controller 27 will be described in detail below. Since details of the sliding mode controller 27 are disclosed in U.S. patent application Ser. No. 09/311353, the sliding mode controller 27 will briefly be described below.

The sliding mode controller 27 sequentially determines the demand differential air-fuel ratio usl as a manipulated variable for manipulating the air-fuel ratio of the internal combustion engine 1 in order to converge the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, i.e., to converge the differential output VO2 of the O 2 sensor 6 to “0”, according to an adaptive sliding mode control process which incorporates an adaptive control law for minimizing the effect of a disturbance, in the normal sliding mode control process. An algorithm for carrying out the adaptive sliding mode control process is constructed as follows:

A switching function required for the adaptive sliding mode control process of the sliding mode controller 27 and a hyperplane defined by the switching function (also referred to as a slip plane) will first be described below.

According to a basic concept of the sliding mode control process in the present embodiment, the differential output VO2(k) from the O 2 sensor 6 in each control cycle and the differential output VO2(k−1) in each preceding control cycle are used, and a switching function σ for the sliding mode control process is established according to the following equation (15). The switching function σ is defined by a linear function having as components the present and past time-series data VO2(k), VO2(k−1) of the differential output VO2 of the O 2 sensor 6 . The vector X defined according to the equation (15) as a vector having the differential outputs VO2 (k), VO2(k−1) as its components will hereinafter be referred to as a state quantity X. $$ σ ⁡ ( k ) = s1 · VO2 ⁡ ( k ) + s2 · VO2 ⁡ ( k - 1 ) = S · X ⁢ ⁢ ( S = [ s1 ⁢   ⁢ s2 ] , X = [ VO2 ⁡ ( k ) VO2 ⁡ ( k - 1 ) ] ) ( 15 )

The coefficients s1, s2 relative to the components VO2(k), VO2(k−1) of the switching function σ are established to meet the condition of the following equation (16): $$ - 1 < s2 s1 < 1 ⁢ ⁢ ( when ⁢   ⁢ s1 = 1 , - 1 < s2 < 1 ) ( 16 )  (when s1=1, −1<s2<1)  (16)

In the present embodiment, for the sake of brevity, the coefficient s1 is set to s1=1 (s2/s1=s2), and the coefficient s2 is established to satisfy the condition:

−1<s2<1.

With the switching function σ thus defined, the hyperplane for the sliding mode control process is defined by the equation σ=0. Since the state quantity X is of the second degree, the hyperplane σ=0 is represented by a straight line as shown in FIG. 4 , and, at this time, the hyperplane is called also a switching function.

The time-series data of the estimated differential output VO2 bar determined by the estimator 26 is actually used as the components of the switching function, as described later on.

The adaptive sliding mode control process in this embodiment serves to converge the state quantity X onto the hyperplane σ=0 according to a reaching control law which is a control law for converging the state quantity X=(VO2(k), VO2(k−1)) onto the hyperplane σ=0, i.e., for converging the value of the switching function σ to “0”, and an adaptive control law (adaptive algorithm) which is a control law for compensating for the effect of a disturbance in converging the state quantity X onto the hyperplane σ=0 (mode 1 in FIG. 4 ). While converging the state quantity X onto the hyperplane σ=0 according to an equivalent control input (holding the value of the switching function σ at “0”), the state quantity X is converted to a balanced point on the hyperplane σ=0 where VO2(k)=VO2(k−1)=0, i.e., a point where time-series data VO2/OUT(k), VO2/OUT(k−1) of the output VO2/OUT of the O 2 sensor 6 are equal to the target value VO2/TARGET.

The command differential air-fuel ratio usl to be generated by the sliding mode controller 27 according to the sliding mode control process for converging the state quantity X to the balanced point on the hyperplane σ=0 is expressed as the sum of an equivalent control input ueq to be applied to the object exhaust system E according to the control law for converging the state quantity X onto the hyperplane σ=0, an input urch (hereinafter referred to as a “reaching control law input urch”) to be applied to the object exhaust system E according to the reaching control law, and an input uadp (hereinafter referred to as an “adaptive control law input uadp”) to be applied to the object exhaust system E according to the adaptive control law (see the following equation (17)).

Usl=Ueq+Urch+Uadp  (17)

The equivalent control input ueq, the reaching control law input urch, and the adaptive control law uadp are determined on the basis of the model of the discrete time system expressed by the equation (11) (a model where the differential output kact(k−d1) of the LAF sensor 5 in the equation (1) is replaced with the command differential air-fuel ratio kcmd(k−d) using the total dead time d), as follows:

The equivalent control input ueq which is an input to be applied to the object exhaust system E for converging the state quantity X onto the hyperplane σ=0 is equal to the command differential air-fuel ratio kcmd which satisfies the condition: σ(k+1)=σ(k)=0. Using the equations (11), (15), the equivalent control input ueq which satisfies the above condition is given by the following equation (18): $$ Ueq ⁡ ( k ) =   ⁢ - ( S · B ) - 1 · { S · ( A - 1 ) } · X ⁡ ( k + d ) =   ⁢ - 1 s1b1 · { [ s1 · ( a1 - 1 ) + s2 ] · VO2 ⁡ ( k + d ) +   ⁢ ( s1 · a2 - s2 ) · VO2 ⁡ ( k + d - 1 ) } ( 18 )

The equation (18) is a basic formula for determining the equivalent control input ueq in each control cycle.

According to present embodiment, the reaching control law input urch is basically determined according to the following equation (19): $$ Urch ⁡ ( k ) = - ( S · B ) - 1 · F · σ ⁡ ( k + d ) = - 1 s1b1 · F · σ ⁡ ( k + d ) ( 19 )

Specifically, the reaching control law input urch is determined in proportion to the value σ((k+d) of the switching function σ after the total dead time d, in view of the effect of the total dead time d.

The coefficient F in the equation (19) which determines the gain of the reaching control law is established to satisfy the condition expressed by the following equation (20):

0<F<2  (20)

The value of the switching function σ may possibly vary in an oscillating fashion (so-called chattering) with respect to “0”. In order to suppress such chattering, it is preferable that the coefficient F relative to the reaching control law input urch be established to further satisfy the condition of the following equation (21):

0<F<1  (21)

The adaptive control law input uadp is basically determined according to the following equation (22) (ΔT in the equation (22) represents the period of the control cycles of the exhaust-side control unit 7 a ): $$ Uadp ⁡ ( k ) = - ( S · B ) - 1 · G · ∑ i = 0 k + d ⁢ ( σ ⁡ ( i ) · Δ ⁢   ⁢ T ) = - 1 s1b1 · G · ∑ i = 0 k + d ⁢ ( σ ⁡ ( i ) · Δ ⁢   ⁢ T ) ( 22 )

The adaptive control law input uadp is determined in proportion to an integrated value (which corresponds to an integral of the values of the switching function σ) over control cycles of values of the switching function σ until after the total dead time d, in view of the effect of the total dead time d.

The coefficient G (which determines the gain of the adaptive control law) in the equation (22) is established to satisfy the condition of the following equation (23): $$ G = J · 2 - F Δ ⁢   ⁢ T ⁢ ⁢ ( 0 < J < 2 ) ( 23 )

A specific process of deriving conditions for establishing the equations (16), (20), (21), (23) is described in detail in Japanese patent application No. 11-93741 and U.S. patent application Ser. No. 09/153032, and will not be described in detail below.

The demand differential air-fuel ratio usl as an input to be given to the object exhaust system E for converging the output signal VO2/OUT of the O 2 sensor 6 to its target value VO2/TARGET may basically be determined as the sum (ueq+urch+uadp) of the equivalent control input ueq, the reaching control law input urch, and the adaptive control law uadp determined according to the respective equations (18), (19), (22). However, the differential outputs VO2(K+d), VO2(k+d−1) of the O 2 sensor 6 and the value σ(k+d) of the switching function σ, etc. used in the equations (18), (19), (22) cannot directly be obtained as they are values in the future.

According to present embodiment, therefore, the sliding mode controller 27 uses the estimated differential outputs VO2(k+d) bar, VO2(k+d−1) bar determined by the estimator 26 , instead of the differential outputs VO2(K+d), VO2(k+d−1) from the O 2 sensor 6 for determining the equivalent control input ueq according to the equation (18), and calculates the equivalent control input ueq in each control cycle according to the following equation (24): $$ Ueq ⁡ ( k ) =   ⁢ - 1 s1b1 ⁢ { [ s1 · ( a1 - 1 ) + s2 ] · VO2 _ ⁡ ( k + d ) +   ⁢ ( s1 · a2 - s2 ) · VO2 _ ⁡ ( k + d - 1 ) } ( 24 )

According to present embodiment, furthermore, the sliding mode controller 27 actually uses time-series data of the estimated differential output VO2 bar sequentially determined by the estimator 26 as described as a state quantity to be controlled. That is, the sliding mode controller 27 defines a switching function σ bar according to the following equation (25) (the linear function σ bar corresponds to time-series data of the differential output VO2 in the equation (15) which is replaced with time-series data of the estimated differential output VO2 bar), in place of the switching function σ established according to the equation (15):

{overscore (σ(k))}=s1·{overscore (VO2+L )}(k)+s2·{overscore (VO2+L )}(k−1)  (25)

The sliding mode controller 27 calculates the reaching control law input urch in each control cycle according to the following equation (26), using the value of the switching function σ bar represented by the equation (25), rather than the value of the switching function σ for determining the reaching control law input urch according to the equation (19): $$ Urch ⁡ ( k ) = - 1 s1b1 · F · σ _ ⁡ ( k + d ) ( 26 )

Similarly, the sliding mode controller 27 calculates the adaptive control law input uadp in each control cycle according to the following equation (27), using the value of the switching function a bar represented by the equation (25), rather than the value of the switching function σ for determining the adaptive control law input uadp according to the equation (22): $$ Uadp ⁡ ( k ) = - 1 s1b1 · G · ∑ i = 0 k + d ⁢ ( σ _ ⁡ ( i ) · Δ ⁢   ⁢ T ) ( 27 )

The latest identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat which have been determined by the identifier 25 are basically used as the gain coefficients a1, a2, b1 that are required to calculate the equivalent control input ueq, the reaching control law input urch, and the adaptive control law input uadp according to the equations (24), (26), (27).

The sliding mode controller 27 determines the sum of the equivalent control input ueq, the reaching control law input urch, and the adaptive control law input uadp determined according to the equations (24), (26), (27), as the demand differential air-fuel ratio usl (see the equation (17)). The conditions for establishing the coefficients s1, s2, F, G used in the equations (24), (26), (27) are as described above.

The demand differential air-fuel ratio usl determined by the sliding mode controller 27 as described above signifies an input to be given to the object exhaust system 1 for converting the estimated differential output VO2 bar from the O 2 sensor 6 to “0”, and as a result, for converting the output VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, i.e., a target value for the difference between the air-fuel ratio of the internal combustion engine 1 and the reference value FLAF/BASE.

The above process is a calculating process (algorithm) for generating the demand differential air-fuel ratio usl in each control cycle by the sliding mode controller 27 .

As described above, the limiter 30 (see FIG. 3 ) in the target air-fuel ratio generator 28 of the exhaust-side control unit 7 a limits the demand differential air-fuel ratio usl generated by the sliding mode controller 27 thereby to determine the command differential air-fuel ratio kcmd (normally kcmd=usl). Then, the adder 31 adds the air-fuel ratio reference value FLAF/BASE to the determined command differential air-fuel ratio kcmd for thereby determining the target air-fuel ratio KCMD (=kcmd+FLAF/BASE) to be given to the engine-side control unit 7 b according to the following equation (28):

 KCMD(k)=kcmd(k)+FLAF/BASE  (28)

In this embodiment, the limiter 30 sequentially variably sets up an allowable range for limiting the demand differential air-fuel ratio usl depending on the demand differential air-fuel ratio usl and the operating state of the internal combustion engine 1 , as described later on.

In the present embodiment, furthermore, the air-fuel ratio reference value FLAF/BASE used as a reference for the demand differential air-fuel ratio usl and the command differential air-fuel ratio kcmd is sequentially variably established by the reference value setting unit 11 depending on the adaptive control law uadp which is an adaptive control law component of the demand differential air-fuel ratio usl generated by the sliding mode controller 27 , as described later on.

The general feedback controller 15 of the engine-side control unit 7 b , particularly, the adaptive controller 18 , will be described below.

As shown in FIG. 1 , the general feedback controller 15 effects a feedback control process to converge the output KACT (detected air-fuel ratio) from the LAF sensor 5 to the target air-fuel ratio KCMD as described above. If such a feedback control process were carried out under the known PID control only, it would be difficult keep stable controllability against dynamic behavioral changes including changes in the operating conditions of the internal combustion engine 1 , characteristic changes due to aging of the internal combustion engine 1 , etc.

The adaptive controller 18 is a recursive-type controller which makes it possible to carry out a feedback control process while compensating for dynamic behavioral changes of the internal combustion engine 1 . As shown in FIG. 5 , the adaptive controller 18 comprises a parameter adjuster 32 for establishing a plurality of adaptive parameters using the parameter adjusting law proposed by I. D. Landau, et al., and a manipulated variable calculator 33 for calculating the feedback manipulated variable KSTR using the established adaptive parameters.

The parameter adjuster 32 will be described below. According to the parameter adjusting law proposed by I. D. Landau, et al., when polynomials of the denominator and the numerator of a transfer function B(Z −1 )/A(Z −1 ) of a discrete-system object to be controlled are generally expressed respectively by equations (29), (30), given below, an adaptive parameter θ hat (j) (j indicates the number of a control cycle) established by the parameter adjuster 32 is represented by a vector (transposed vector) according to the equation (31) given below. An input ξ(j) to the parameter adjuster 32 is expressed by the equation (32) given below. In the present embodiment, it is assumed that the internal combustion engine 1 , which is an object to be controlled by the general feedback controller 15 , is considered to be a plant of a first-order system having a dead time d p corresponding to three control cycles (a time corresponding to three combustion cycles of the internal combustion engine 1 ), and m=n=1, d p =3 in the equations (29)-(32), and five adaptive parameters s 0 , r 1 , r 2 , r 3 , b 0 are established (see FIG. 5 ). In the upper and middle expressions of the equation (32), u s , y s generally represent an input (manipulated variable) to the object to be controlled and an output (controlled variable) from the object to be controlled. In the present embodiment, the input is the feedback manipulated variable KSTR and the output from the object (the internal combustion engine 1 ) is the output KACT (detected air-fuel ratio) from the LAF sensor 5 , and the input ξ(j) to the parameter adjuster 32 is expressed by the lower expression of the equation (32) (see FIG. 5 ).

A(Z −1 )=1+a1Z −1 +. . . +anZ −n   (29)

B(Z −1 )=b0+b1Z −1 +. . . +bmZ −m   (30)

$$ θ ^ T ⁡ ( j ) = [ b ^ 0 ⁡ ( j ) , B ^ R ⁡ ( Z - 1 , j ) , S ^ ⁡ ( Z - 1 , j ) ] = [ b 0 ⁡ ( j ) , r 1 ⁡ ( j ) , … ⁢   , r m + d p - 1 ⁡ ( j ) , s o ⁡ ( j ) , … ⁢   , s n - 1 ⁡ ( j ) ] = [ b 0 ⁡ ( j ) , r 1 ⁡ ( j ) , r 2 ⁡ ( j ) , r 3 ⁡ ( j ) , s o ⁡ ( j ) ] ( 31 ) ζ T ⁡ ( j ) =   ⁢ [ u s ⁡ ( j ) , … ⁢   , u s ⁡ ( j - m - dp + 1 ) ,   ⁢ y s ⁡ ( j ) , … ⁢   , y s ⁡ ( j - n + 1 ) ] =   ⁢ [ u s ⁡ ( j ) , u s ⁡ ( j - 1 ) , u s ⁡ ( j - 2 ) , u s ⁡ ( j - 3 ) , y s ⁡ ( j ) ] =   ⁢ [ KSTR ⁡ ( j ) , KSTR ⁡ ( j - 1 ) , KSTR ⁡ ( j - 2 ) ,   ⁢ KSTR ⁡ ( j - 3 ) , KACT ⁡ ( j ) ] ( 32 )

The adaptive parameter θ hat expressed by the equation (31) is made up of a scalar quantity element b 0 hat (Z −1 ,j) for determining the gain of the adaptive controller 18 , a control element B R hat (Z −1 ,j) expressed using a manipulated variable, and a control element S hat (Z −1 ,j) expressed using a controlled variable, which are expressed respectively by the following equations (33)˜(35) (see the block of the manipulated variable calculator 33 shown in FIG. 5 ): $$ b ^ 0 - 1 ⁡ ( j ) = 1 b 0 ( 33 ) B ^ R ⁡ ( Z - 1 , j ) = r 1 ⁢ Z - 1 + r 2 ⁢ Z - 2 + … + r m + d p - 1 ⁢ Z - ( n + dp - 1 ) = r 1 ⁢ Z - 1 + r 2 ⁢ Z - 2 + r 3 ⁢ Z - 3 ( 34 ) S ^ ⁡ ( Z - 1 , j ) = s 0 + s 1 ⁢ Z - 1 + … + s n - 1 ⁢ Z ( n - 1 ) = s 0 ( 35 )

The parameter adjuster 32 establishes coefficients of the scalar quantity element and the control elements, described above, and supplies them as the adaptive parameter θ hat expressed by the equation (31) to the manipulated variable calculator 33 . The parameter adjuster 32 calculates the adaptive parameter θ hat so that the output KACT from the LAF sensor 5 will agree with the target air-fuel ratio KCMD, using time-series data of the feedback manipulated variable KSTR from the present to the past and the output KACT.

Specifically, the parameter adjuster 32 calculates the adaptive parameter θ hat according to the following equation (36):

{circumflex over (θ)}(j)={circumflex over (θ)}(j−1)+Γ(j−1)·ξ(j−d p )·e*(j)  (36)

where Γ (j) represents a gain matrix (whose degree is indicated by m+n+d p ) for determining a rate of establishing the adaptive parameter θ hat, and e*(j) an estimated error of the adaptive parameter θ hat. Γ (j) and e*(j) are expressed respectively by the following recursive formulas (37), (38): $$ Γ ⁡ ( j ) =   ⁢ 1 λ 1 ⁡ ( j ) [ Γ ⁡ ( j - 1 ) -   ⁢ λ 2 ⁡ ( j ) · Γ ⁡ ( j - 1 ) · ζ ⁡ ( j - d p ) · ζ T ⁡ ( j - d p ) · Γ ⁡ ( j - 1 ) λ 2 ⁡ ( j ) + λ 2 ⁡ ( j ) · ζ T ⁡ ( j - d p ) · Γ ⁡ ( j - 1 ) · ζ ⁡ ( j - d p ) ] ⁢ ⁢ where ⁢   ⁢ 0 < λ 1 ⁡ ( j ) ≤ 1 , 0 ≤ λ 2 ⁡ ( j ) < 2 , Γ ⁡ ( 0 ) > 0. ( 37 ) e * ⁡ ( j ) = D ⁡ ( Z - 1 ) · KACT ⁡ ( j ) - θ ^ T ⁡ ( j - 1 ) · ζ ⁡ ( j - d p ) 1 + ζ T ⁡ ( j - d p ) · Γ ⁡ ( j - 1 ) · ζ ⁡ ( j - d p ) ( 38 )

where D(Z −1 ) represents an asymptotically stable polynomial for adjusting the convergence. In the present embodiment, D(Z −1 )=1.

Various specific algorithms including the degressive gain algorithm, the variable gain algorithm, the fixed trace algorithm, and the fixed gain algorithm are obtained depending on how λ 1 (j), λ 2 (j) in the equation (37) are selected. For a time-dependent plant such as a fuel injection process, an air-fuel ratio, or the like of the internal combustion engine 1 , either one of the degressive gain algorithm, the variable gain algorithm, the fixed gain algorithm, and the fixed trace algorithm is suitable.

Using the adaptive parameter θ hat (s 0 , r 1 , r 2 , r 3 , b 0 ) established by the parameter adjuster 32 and the target air-fuel ratio KCMD calculated by the exhaust-side control unit 7 a , the manipulated variable calculator 33 determines the feedback manipulated variable KSTR according to a recursive formula expressed by the following equation (39): $$ KSTR = KCMD ⁡ ( j ) - s 0 · KACT ⁡ ( j ) - r 1 · KSTR ⁡ ( j - 1 ) - r 2 · KSTR ⁡ ( j - 2 ) - r 3 · KSTR ⁡ ( j - 3 ) b 0 ( 39 )

The manipulated variable calculator 33 shown in FIG. 5 represents a block diagram of the calculations according to the equation (39).

The feedback manipulated variable KSTR determined according to the equation (39) becomes the target air-fuel ratio KCMD insofar as the output KACT of the LAF sensor 5 agrees with the target air-fuel ratio KCMD. Therefore, the feedback manipulated variable KSTR is divided by the target air-fuel ratio KCMD by the divider 19 for thereby determining the feedback manipulated variable kstr that can be used as the feedback correction coefficient KFB.

As is apparent from the foregoing description, the adaptive controller 18 thus constructed is a recursive-type controller taking into account dynamic behavioral changes of the internal combustion engine 1 which is an object to be controlled. Stated otherwise, the adaptive controller 18 is a controller described in a recursive form to compensate for dynamic behavioral changes of the internal combustion engine 1 , and more particularly a controller having a recursive-type adaptive parameter adjusting mechanism.

A recursive-type controller of this type may be constructed using an optimum regulator. In such a case, however, it generally has no parameter adjusting mechanism. The adaptive controller 18 constructed as described above is suitable for compensating for dynamic behavioral changes of the internal combustion engine 1 .

The details of the adaptive controller 18 have been described above.

The PID controller 17 , which is provided together with the adaptive controller 18 in the general feedback controller 20 , calculates a proportional term (P term), an integral term (I term), and a derivative term (D term) from the difference between the output KACT of the LAF sensor 5 and the target air-fuel ratio KCMD, and calculates the total of those terms as the feedback manipulated variable KLAF, as is the case with the general PID control process. In the present embodiment, the feedback manipulated variable KLAF is set to “1” when the output KACT of the LAF sensor 5 agrees with the target air-fuel ratio KCMD by setting an initial value of the integral term (I term) to “1”, so that the feedback manipulated variable KLAF can be used as the feedback correction coefficient KFB for directly correcting the fuel injection quantity. The gains of the proportional term, the integral term, and the derivative term are determined from the rotational speed and intake pressure of the internal combustion engine 1 using a predetermined map.

The switcher 20 of the general feedback controller 15 outputs the feedback manipulated variable KLAF determined by the PID controller 17 as the feedback correction coefficient KFB for correcting the fuel injection quantity if the combustion in the internal combustion engine 1 tends to be unstable as when the temperature of the coolant of the internal combustion engine 1 is low, the internal combustion engine 1 rotates at high speeds, or the intake pressure is low, or if the output KACT of the LAF sensor 5 is not reliable due to a response delay of the LAF sensor 5 as when the target air-fuel ratio KCMD changes largely or immediately after the air-fuel ratio feedback control process has started, or if the internal combustion engine 1 operates highly stably as when it is idling and hence no high-gain control process by the adaptive controller 18 is required. Otherwise, the switcher 20 outputs the feedback manipulated variable kstr which is produced by dividing the feedback manipulated variable KSTR determined by the adaptive controller 18 by the target air-fuel ratio KCMD, as the feedback correction coefficient KFB for correcting the fuel injection quantity. This is because the adaptive controller 18 effects a high-gain control process and functions to converge the output KACT of the LAF sensor 5 quickly to the target air-fuel ratio KCMD, and if the feedback manipulated variable KSTR determined by the adaptive controller 18 is used when the combustion in the internal combustion engine 1 is unstable or the output KACT of the LAF sensor 5 is not reliable, then the air-fuel ratio control process tends to be unstable.

Such operation of the switcher 20 is disclosed in detail in Japanese laid-open patent publication Ser. No. 8-105345 and U.S. Pat. No. 5,558,075, and will not be described in detail below.

Operation of the plant control system will be described below.

First, a process, carried out by the engine-side control unit 7 b , of calculating an output fuel injection quantity #nTout (n=1, 2, 3, 4) for each of the cylinders of the internal combustion engine 1 for controlling the air-fuel ratio of the internal combustion engine 1 will be described below with reference to FIG. 6 . The engine-side control unit 7 b calculates an output fuel injection quantity #nTout (n=1, 2, 3, 4) for each of the cylinders in synchronism with a crankshaft angle period of the internal combustion engine 1 as follows:

In FIG. 7 , the engine-side control unit 7 b reads outputs from various sensors including the LAF sensor 5 and the O 2 sensor 6 in STEPa. At this time, the output KACT of the LAF sensor 5 and the output VO2/OUT of the O 2 sensor 6 , including data obtained in the past, are stored in a time-series fashion in a memory (not shown).

Then, the basic fuel injection quantity calculator 8 corrects a fuel injection quantity corresponding to the rotational speed NE and intake pressure PB of the internal combustion engine 1 depending on the effective opening area of the throttle valve, thereby calculating a basic fuel injection quantity Tim in STEPb. The first correction coefficient calculator 9 calculates a first correction coefficient KTOTAL depending on the coolant temperature and the amount by which the canister is purged in STEPc.

The engine control unit 7 b determines an operation mode of the internal combustion engine 1 , and sets a value of a flag f/prism/on which represents whether the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 of the exhaust-side control unit 7 a is to be used to manipulate the air-fuel ratio of the internal combustion engine 1 or not in STEPd.

Operation modes of the internal combustion engine 1 , specifically modes of manipulating the air-fuel ratio of the internal combustion engine 1 , include an operation mode in which the engine-side control unit 7 b uses the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 (which is basically a target air-fuel ratio for converging the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET), and manipulates the air-fuel ratio of the internal combustion engine 1 (i.e., adjust the fuel injection quantity of the internal combustion engine 1 ) in order to the output signal KACT of the LAF sensor (detected air-fuel ratio) to the target air-fuel ratio KCMD, and an operation mode in which the engine-side control unit 7 b manipulates the air-fuel ratio of the internal combustion engine 1 (i.e., adjust the fuel injection quantity of the internal combustion engine 1 ), without using the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 .

The former operation mode is a normal operation mode while the internal combustion engine 1 is in operation (hereinafter referred to as a “normal operation mode”). The latter operation mode represents a plurality of operation modes including an operation mode in which the fuel supplied to the internal combustion engine 1 is cut off (stopped), an operation mode in which the throttle valve (not shown) of the internal combustion engine 1 is fully open, and an operation mode in which the air-fuel ratio of the internal combustion engine 1 is manipulated to generate a lean air-fuel mixture (lean operation mode).

The deciding process of STEPd is a process of determining whether the internal combustion engine 1 is in the normal operation mode or not, i.e., whether the engine-side control unit 7 b is to manipulate the air-fuel ratio of the internal combustion engine 1 using the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 , or not. When the value of the flag f/prism/on is “0”, it means that the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 is not to be used (the internal combustion engine 1 is not in the normal operation mode), and when the value of the flag f/prism/on is “1”, it means that the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 is to be used (the internal combustion engine 1 is in the normal operation mode).

The deciding subroutine of STEPd is shown in detail in FIG. 7 . As shown in FIG. 7 , the engine-side control unit 7 b decides whether the O 2 sensor 6 and the LAF sensor 5 are activated or not respectively in STEPd- 1 , STEPd- 2 . If neither one of the O 2 sensor 6 and the LAF sensor 5 is activated, since detected data from the O 2 sensor 6 and the LAF sensor 5 for use by the target air-fuel ratio generator 28 are not accurate enough, the target air-fuel ratio generator 28 cannot generate an appropriate target air-fuel ratio. Therefore, since the internal combustion engine 1 cannot operate in the normal operation mode, the value of the flag f/prism/on is set to “0” in STEPd- 10 .

Then, the engine-side control unit 7 b decides whether the internal combustion engine 1 is operating with a lean air-fuel mixture or not in STEPd- 3 . The engine control unit 7 b decides whether the ignition timing of the internal combustion engine 1 is retarded for early activation of the catalytic converter 3 immediately after the start of the internal combustion engine 1 or not in STEPd- 4 . The engine control unit 7 b decides whether the throttle valve of the internal combustion engine 1 is fully open or not in STEPd- 5 . The engine control unit 7 b decides whether the supply of fuel to the internal combustion engine 1 is being stopped or not in STEPd- 6 . If either one of the conditions of these steps is satisfied (YES), then since the internal combustion engine 1 is not in the normal operation mode, the value of the flag f/prism/on is set to “0” in STEPd- 10 .

The engine-side control unit 7 b then decides whether the rotational speed NE and the intake pressure PB of the internal combustion engine 1 fall within respective given ranges (normal ranges) or not respectively in STEPd- 7 , STEPd- 8 . If either one of the rotational speed NE and the intake pressure PB does not fall within its given range, then since it is not preferable to control the air-fuel ratio of the internal combustion engine 1 using the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 , the value of the flag f/prism/on is set to “0” in STEPd- 10 .

If the conditions of STEPd- 1 , STEPd- 2 , STEPd- 7 , STEPd- 8 are satisfied, and the conditions of STEPd- 3 , STEPd- 4 , STEPd- 5 , STEPd- 6 are not satisfied (normal), then the engine-side control unit 7 b determines that the internal combustion engine 1 is in the normal operation mode in which the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 is used to manipulate the air-fuel ratio of the internal combustion engine 1 , and sets the value of the flag f/prism/on to “1” in STEPd- 9 .

After the value of the flag f/prism/on is set to “0” in STEPd- 10 , the engine-side control unit 7 b sets, to a given initial value TMSTB, a timer counter tm/stb (count-down timer) for measuring an elapsed time from the time when the value of the flag f/prism/on changes from “0” to “1”, i.e., an elapsed time after the air-fuel ratio of the internal combustion engine 1 starts being manipulated according to the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 , i.e., the control process for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET has started to be carried out, and activates the timer counter tm/stb in STEPd- 11 . The count of the timer counter tm/stb is counted down at a constant period (longer than the period of control cycles of the exhaust-side control unit 7 a ). The initial value TMSTB corresponds to a period in which to inhibit the decision as to the stability of the sliding mode controller 27 .

In FIG. 6 , after the value of the flag f/prism/on has been set, the engine-side control unit 7 b determines the value of the flag f/prism/on in STEPe. If f/prism/on=1, then the engine-side control unit 7 b reads the latest target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 in STEPf. If f/prism/on=0, then the engine-side control unit 7 b sets the target air-fuel ratio KCMD to a predetermined value in STEPg. The predetermined value to be established as the target air-fuel ratio KCMD is determined from the rotational speed NE and intake pressure PB of the internal combustion engine 1 using a predetermined map, for example. Particularly when the engine 1 is in the lean operation mode, the target air-fuel ratio KCMD set in STEPg is an air-fuel ratio in a lean range.

In the local feedback controller 16 , the PID controllers 22 calculates respective feedback correction coefficients #nKLAF in order to eliminate variations between the cylinders, based on actual air-fuel ratios #nA/F of the respective cylinders which have been estimated from the output KACT of the LAF sensor 5 by the observer 21 , in STEPh. Then, the general feedback controller 15 calculates a feedback correction coefficient KFB in STEPi.

Depending on the operating conditions of the internal combustion engine 1 , the switcher 20 selects either the feedback manipulated variable KLAF determined by the PID controller 17 or the feedback manipulated variable kstr which has been produced by dividing the feedback manipulated variable KSTR determined by the adaptive controller 18 by the target air-fuel ratio KCMD (normally, the switcher 20 selects the feedback manipulated variable kstr). The switcher 20 then outputs the selected feedback manipulated variable KLAF or kstr as a feedback correction coefficient KFB.

When switching the feedback correction coefficient KFB from the feedback manipulated variable KLAF from the PID controller 17 to the feedback manipulated variable kstr from the adaptive controller 18 , the adaptive controller 18 determines a feedback manipulated variable KSTR in a manner to hold the correction coefficient KFB to the preceding correction coefficient KFB (=KLAF) as long as in the cycle time for the switching. When switching the feedback correction coefficient KFB from the feedback manipulated variable kstr from the adaptive controller 18 to the feedback manipulated variable KLAF from the PID controller 17 , the PID controller 17 calculates a present correction coefficient KLAF in a manner to regard the feedback manipulated variable KLAF determined by itself in the preceding cycle time as the preceding correction coefficient KFB (=kstr).

After the feedback correction coefficient KFB has been calculated, the second correction coefficient calculator 10 calculates in STEPj a second correction coefficient KCMDM depending on the target air-fuel ratio KCMD determined in STEPf or STEPg.

Then, the engine-side control unit 7 b multiplies the basic fuel injection quantity Tim, determined as described above, by the first correction coefficient KTOTAL, the second correction coefficient KCMDM, the feedback correction coefficient KFB, and the feedback correction coefficients #nKLAF of the respective cylinders, determining output fuel injection quantities #nTout of the respective cylinders in STEPk. The output fuel injection quantities #nTout are then corrected for accumulated fuel particles on intake pipe walls of the internal combustion engine 1 by the fuel accumulation corrector 23 in STEPm. The corrected output fuel injection quantities #nTout are applied to the non-illustrated fuel injectors of the internal combustion engine 1 in STEPn.

In the internal combustion engine 1 , the fuel injectors inject fuel into the respective cylinders according to the respective output fuel injection quantities #nTout.

The above calculation of the output fuel injection quantities #nTout and the fuel injection of the internal combustion engine 1 are carried out in successive cycle times synchronous with the crankshaft angle period of the internal combustion engine 1 for controlling the air-fuel ratio of the internal combustion engine 1 in order to converge the output KACT of the LAF sensor 5 (the detected air-fuel ratio) to the target air-fuel ratio KCMD. While the feedback manipulated variable kstr from the adaptive controller 18 is being used as the feedback correction coefficient KFB, the output KACT of the LAF sensor 5 is quickly converged to the target air-fuel ratio KCMD with high stability against behavioral changes such as changes in the operating conditions of the internal combustion engine 1 or characteristic changes thereof. A response delay of the internal combustion engine 1 is also appropriately compensated for.

Concurrent with the above air-fuel ratio manipulation for the internal combustion engine 1 , the target air-fuel ratio generator 28 of the exhaust-side control unit 7 a executes a main routine shown in FIG. 8 in control cycles of a constant period.

As shown in FIG. 8 , the target air-fuel ratio generator 28 decides whether its own processing (the processing of the identifier 25 , the estimator 26 , the sliding mode controller 27 , and the limiter 30 ) is to be executed or not, and sets a value of a flag f/prism/cal indicative of whether the processing is to be executed or not in STEP 1 . When the value of the flag f/prism/cal is “0”, it means that the processing of the target air-fuel ratio generator 28 is not to be executed, and when the value of the flag f/prism/cal is “1”, it means that the processing of the target air-fuel ratio generator 28 is to be executed.

The deciding subroutine in STEP 1 is shown in detail in FIG. 9 . As shown in FIG. 9 , the target air-fuel ratio generator 28 decides whether the O 2 sensor 6 and the LAF sensor 5 are activated or not respectively in STEP 1 - 1 , STEP 1 - 2 . If neither one of the O 2 sensor 6 and the LAF sensor 5 is activated, since detected data from the O 2 sensor 6 and the LAF sensor 5 for use by the target air-fuel ratio generator 28 are not accurate enough, the value of the flag f/prism/cal is set to “0” in STEP 1 - 6 . Then, in order to initialize the identifier 25 as described later on, the value of a flag f/id/reset indicative of whether the identifier 25 is to be initialized or not is set to “1” in STEP 1 - 7 . When the value of the flag f/id/reset is “1”, it means that the identifier 25 is to be initialized, and when the value of the flag f/id/reset is “0”, it means that the identifier 25 is not to be initialized.

The target air-fuel ratio generator 28 decides whether the internal combustion engine 1 is operating with a lean air-fuel mixture or not in STEP 1 - 3 . The target air-fuel ratio generator 28 decides whether the ignition timing of the internal combustion engine 1 is retarded for early activation of the catalytic converter 3 immediately after the start of the internal combustion engine 1 or not in STEP 1 - 4 . If the conditions of these steps are satisfied, then since the target air-fuel ratio KCMD calculated to converge the output VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET is not used for the fuel control for the internal combustion engine 1 , the value of the flag f/prism/cal is set to “0” in STEP 1 - 6 , and the value of the flag f/id/reset is set to “1” in order to initialize the identifier 25 in STEP 1 - 7 .

If the conditions of STEP 1 - 1 , STEP 1 - 2 are satisfied, and the conditions of STEP 1 - 3 , STEP 1 - 4 are not satisfied, then the value of the flag f/prism/cal is set to “1” in STEP 1 - 5 .

In FIG. 8 , after the above deciding subroutine, the target air-fuel ratio generator 28 decides whether a process of identifying (updating) the gain coefficients a1, a1, b1 with the identifier 25 is to be executed or not, and sets a value of a flag f/id/cal indicative of whether the process of identifying (updating) the gain coefficients a1, a1, b1 is to be executed or not in STEP 2 .

In STEP 2 , the target air-fuel ratio generator 28 decides whether the throttle valve of the internal combustion engine 1 is substantially fully open or not, and whether the supply of fuel to the internal combustion engine 1 is being stopped or not. If either one of the conditions of these steps is satisfied, then since it is difficult to identify the gain coefficients a1, a1, b1 appropriately, the value of the flag f/id/cal is set to “0”. If neither one of the conditions of these steps is satisfied, then the value of the flag f/id/cal is set to “1” to identify (update) the gain coefficients a1, a2, b1 with the identifier 25 .

Referring back to FIG. 8 , the target air-fuel ratio generator 28 calculates the latest differential outputs kact(k) (=KACT−FLAF/BASE), VO2(k) (=VO2/OUT−VO2/TARGET) respectively from the subtractors 12 , 13 in STEP 3 . Specifically, the subtractors 12 , 13 select latest ones of the time-series data read and stored in the non-illustrated memory in STEPa shown in FIG. 6 , and calculate the differential outputs kact(k), VO2(k). The air-fuel ratio reference value FLAF/BASE required to calculate the differential output kact(k) is of a latest nature set by the reference value setting unit 11 as described later on. The data of the differential outputs kact(k), VO2(k), as well as the data thereof calculated in the past, are stored in a time-series manner in a memory (not shown) in the exhaust-side control unit 7 a.

Then, in STEP 4 , the target air-fuel ratio generator 28 determines the value of the flag f/prism/cal set in STEP 1 . If the value of the flag f/prism/cal is “0”, i.e., if the processing of the target air-fuel ratio generator 28 is not to be executed, then the target air-fuel ratio generator 28 forcibly sets the command differential air-fuel ratio kcmd for determining the target air-fuel ratio KCMD to a predetermined value in STEP 13 . The predetermined value may be a fixed value (e.g., “0”) or the value of the command differential air-fuel ratio kcmd determined in a preceding control cycle. After the command differential air-fuel ratio kcmd is set to the predetermined value in STEP 13 , the adder 31 adds the air-fuel ratio reference value FLAF/BASE (the latest one set by the reference value setting unit 11 ) to the command differential air-fuel ratio kcmd for thereby determining a target air-fuel ratio KCMD in the present control cycle in STEP 11 . Then, the reference value setting unit 11 carries out a process of setting the air-fuel ratio reference value FLAF/BASE, as described later on, in STEP 12 , after which the processing in the present control cycle is finished.

If the value of the flag f/prism/cal is “1” in STEP 4 , i.e., if the processing of the target air-fuel ratio generator 28 is to be executed, then the target air-fuel ratio generator 28 effects the processing of the identifier 25 in STEP 5 .

The processing subroutine of STEP 5 is shown in detail in FIG. 10 .

The identifier 25 determines the value of the flag f/id/cal set in STEP 2 in STEP 5 - 1 . If the value of the flag f/id/cal is “0” (the throttle valve of the internal combustion engine 1 is substantially fully open or the fuel supply of the internal combustion engine 1 is being cut off), then since the process of identifying the gain coefficients a1, a2, b1 with the identifier 25 is not carried out, control immediately goes back to the main routine shown in FIG. 8 .

If the value of the flag f/id/cal is “1”, then the identifier 25 determines the value of the flag f/id/reset set in STEP 1 with respect to the initialization of the identifier 25 in STEP 5 - 2 . If the value of the flag f/id/reset is “1”, the identifier 25 is initialized in STEP 5 - 3 . When the identifier 25 is initialized, the identified gain coefficients a1 hat, a2 hat, b1 hat are set to predetermined initial values (the identified gain coefficient vector Θ according to the equation (4) is initialized), and the elements of the matrix P (diagonal matrix) according to the equation (9) are set to predetermined initial values. The value of the flag f/id/reset is reset to “0”.

Then, the identifier 25 calculates the identified differential output VO2(k) hat from exhaust system model (see the equation (3)) which is expressed using the present identified gain coefficients a1(k−1) hat, a2(k−1) hat, b1(k−1) hat, according to the equation (3) or the equation (6) equivalent thereto, using the past data VO2(k−1), VO2(k−2), kact(k−d−1) of the differential outputs VO2, kact calculated in each control cycle in STEP 3 , and the identified gain coefficients a1(k−1) hat, a2(k−1) hat, b1(k−1) hat, in STEP 5 - 4 .

The identifier 25 then calculates the vector Kθ(k)·to be used in determining the new identified gain coefficients a1 hat, a2 hat, b1 hat according to the equation (9) in STEP 5 - 5 . Thereafter, the identifier 25 calculates the identified error id/e, i.e., the difference between the identified differential output VO2 hat from the O 2 sensor 6 in the exhaust system model and the actual differential output VO2 (see the equation (7)), in STEP 5 - 6 .

The identified error id/e obtained in STEP 5 - 6 may basically be calculated according to the equation (7). In the present embodiment, however, a value (=VO2−VO2 hat) calculated according to the equation (7) from the differential output VO2 acquired in each control cycle in STEP 3 (see FIG. 8 ), and the identified differential output VO2 hat calculated in each control cycle in STEP 5 - 4 is filtered with low-pass characteristics to calculate the identified error id/e.

This is because since the object exhaust system E including the catalytic converter 3 generally has low-pass characteristics, it is preferable to attach importance to the low-frequency behavior of the object exhaust system E in appropriately identifying the gain coefficients a1, a2, b1 of the exhaust system model.

Both the differential output VO2 and the identified differential output VO2 hat may be filtered with the same low-pass characteristics. For example, after the differential output VO2 and the identified differential output VO2 hat have separately been filtered, the equation (7) may be calculated to determine the identified error id/e. The above filtering is carried out by a moving average process which is a digital filtering process, for example.

After the identifier 25 has calculated the identified error id/e, the identifier 25 calculates a new identified gain coefficient vector Θ(k), i.e., new identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat, according to the equation (8) using the identified error id/e and Kθ calculated in SETP 5 - 5 in STEP 5 - 7 .

After having calculated the new identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat, the identifier 25 further limits the values of the gain coefficients a1 hat, a2 hat, b1 hat (elements of the identified gain coefficient vector Θ), are limited to meet predetermined conditions in STEP 5 - 8 .

The predetermined conditions for limiting the values of the identified gain coefficients a1 hat, a2 hat, b1 hat include a condition (hereinafter referred to as a “first limiting condition”) for limiting combinations of the values of the identified gain coefficients a1 hat, a2 hat relative to a predetermined combination, and a condition (hereinafter referred to as a “second limiting condition”) for limiting the value of the identified gain coefficient b1 hat.

Prior to describing the first and second limiting conditions and the specific processing details of STEP 5 - 8 , the reasons for limiting the values of the identified gain coefficients a1 hat, a2 hat, b1 hat will be described below.

The inventors of the present invention have found that if the values of the identified gain coefficients a1 hat, a2 hat, b1 hat are not particularly limited, while the output signal VO2/OUT of the O 2 sensor 6 is being stably controlled at the target value VO2/TARGET, there are developed a situation in which the demand differential air-fuel ratio usl determined by the sliding mode controller 27 and the target air-fuel ratio KCMD change smoothly with time, and a situation in which the demand differential air-fuel ratio usl and the target air-fuel ratio KCMD oscillate with time at a high frequency. Neither of these situations poses problems in controlling the output VO2/OUT of the O 2 sensor 6 at the target value VO2/TARGET. However, the situation in which the target air-fuel ratio KCMD oscillates with time at a high frequency is not preferable in smoothly operating the internal combustion engine 1 .

A study of the above phenomenon by the inventors has shown that whether the demand differential air-fuel ratio usl and the target air-fuel ratio KCMD change smoothly or oscillate at a high frequency is affected by the combinations of the values of the identified gain coefficients a1 hat, a2 hat identified by the identifier 25 and the value of the identified gain coefficient b1 hat.

In the present embodiment, the first and second limiting conditions are established appropriately, and the combinations of the values of the identified gain coefficients a1 hat, a2 hat and the value of the identified gain coefficient b1 hat are appropriately limited to eliminate the situation in which the target air-fuel ratio KCMD oscillates at a high frequency.

According to the present embodiment, the first and second limiting conditions are established as follows:

With respect to the first limiting condition for limiting the values of the identified gain coefficients a1 hat, a2 hat, the study by the inventors indicates that obtaining the demand differential air-fuel ratio usl and the target air-fuel ratio KCMD is closely related to combinations of the coefficient values α1, α2 in the equations (12)-(14) which are determined by the values of the gain coefficients a1, a2, i.e., the coefficient values α1, α2 used for the estimator 26 to determine the estimated differential output VO2(k+d) bar (the coefficient values α1, α2 are the first-row, first-column element and the first-row, second-column element of the matrix A d which is a power of the matrix A defined by the equation (12)).

Specifically, as shown in FIG. 11 , when a coordinate plane whose coordinate components are represented by the coefficient values α1, α2 is established, if a point on the coordinate plane which is determined by a combination of the coefficient values α1, α2 lies in a hatched range, which is surrounded by a triangle Q 1 Q 2 Q 3 (including the boundaries) and will hereinafter be referred to as an “estimating coefficient stable range”, then the demand differential air-fuel ratio usl and the target air-fuel ratio KCMD tend to be smooth.

Therefore, the combinations of the values of the gain coefficients a1, a2 identified by the identifier 25 , i.e., the combinations of the values of the identified gain coefficients a1 hat, a2 hat, should be limited such that the point on the coordinate plane shown in FIG. 11 which corresponds to the combination of the coefficient values α1, α2 determined by the values of the gain coefficients a1, a2 or the values of the identified gain coefficients a1 hat, a2 hat will lie within the estimating coefficient stable range.

In FIG. 11 , a triangular range Q 1 Q 4 Q 3 on the coordinate plane which contains the estimating coefficient stable range is a range that determines combinations of the coefficient values α1, α2 which makes theoretically stable a system defined according to the following equation (40), i.e., a system defined by an equation similar to the equation (12) except that VO2(k), VO2(k−1) on the right side of the equation (12) are replaced respectively with VO2(k) bar, VO2(k−1) bar (VO2(k) bar, VO2(k−1) bar mean respectively an estimated differential output determined in each control cycle by the estimator 26 and an estimated differential output determined in a preceding cycle by the estimator 26 ). $$ VO2 _ ⁡ ( k + d ) =   ⁢ α ⁢ 1 · VO2 _ ⁡ ( k ) + α ⁢ 2 · VO2 _ ⁡ ( k - 1 ) +   ⁢ ∑ j = 1 d ⁢ β j · kcmd ⁡ ( k - j ) ( 40 )

The condition for the system defined according to the equation (40) to be stable is that a pole of the system (which is given by the following equation (41)) exists in a unit circle on a complex plane:

Pole of the system according to the equation (40) $$ = α ⁢ 1 ± α ⁢ 1 2 + 4 · α ⁢ 2 2 ( 41 )

The triangular range Q 1 Q 4 Q 3 shown in FIG. 11 is a range for determining the combinations of the coefficient values α1, α2 which satisfy the above condition. Therefore, the estimating coefficient stable range is a range indicative of those combinations where α1≧0 of the combinations of the coefficient values α1, α2 which make stable the system defined by the equation (40).

Since the coefficient values α1, α2 are determined by a combination of the values of the gain coefficients a1, a2, a combination of the values of the gain coefficients a1, a2 is determined by a combination of the coefficient values α1, α2. Therefore, the estimating coefficient stable range shown in FIG. 11 which determines preferable combinations of the coefficient values α1, α2 can be converted into a range on a coordinate plane shown in FIG. 12 whose coordinate components are represented by the gain coefficients a1, a2. Specifically, the estimating coefficient stable range shown in FIG. 11 is converted into a range enclosed by the imaginary lines in FIG. 12 , which is a substantially triangular range having an undulating lower side and will hereinafter be referred to as an “identifying coefficient stable range”, on the coordinate plane shown in FIG. 12 . Stated otherwise, when a point on the coordinate plane shown in FIG. 12 which is determined by a combination of the values of the gain coefficients a1, a2 resides in the identifying coefficient stable range, a point on the coordinate plane shown in FIG. 11 which corresponds to the combination of the coefficient values α1, α2 determined by those values of the gain coefficients a1, a2 resides in the estimating coefficient stable range.

Consequently, the first limiting condition for limiting the values of the identified gain coefficients a1 hat, a2 hat determined by the identifier 25 should preferably be basically established such that a point on the coordinate plane shown in FIG. 12 which is determined by those values of the identified gain coefficients a1 hat, a2 hat reside in the identifying coefficient stable range.

However, since a boundary (lower side) of the identifying coefficient stable range indicated by the imaginary lines in FIG. 12 is of a complex undulating shape, a practical process for limiting the point on the coordinate plane shown in FIG. 12 which is determined by the values of the identified gain coefficients a1 hat, a2 hat is liable to be complex.

For this reason, according to the present embodiment, the identifying coefficient stable range is substantially approximated by a quadrangular range Q 5 Q 6 Q 7 Q 8 enclosed by the solid lines in FIG. 12 , which has straight boundaries and will hereinafter be referred to as an “identifying coefficient limiting range”. As shown in FIG. 12 , the identifying coefficient limiting range is a range enclosed by a polygonal line (including line segments Q 5 Q 6 and Q 5 Q 8 ) expressed by a functional expression |a1|+a2=1, a straight line (including a line segment Q 6 Q 7 ) expressed by a constant-valued functional expression a1=A1L (A1L: constant), and a straight line (including a line segment Q 7 Q 8 ) expressed by a constant-valued functional expression a2=A2L (A2L: constant). The first limiting condition for limiting the values of the identified gain coefficients a1 hat, a2 hat is established such that the point on the coordinate plane shown in FIG. 15 which is determined by those values of the identified gain coefficients a1 hat, a2 hat lies in the identifying coefficient limiting range. Although part of the lower side of the identifying coefficient limiting range deviates from the identifying coefficient stable range, it has experimentally been confirmed that the point determined by the identified gain coefficients a1 hat, a2 hat determined by the identifier 25 does not actually fall in the deviating range. Therefore, the deviating range will not pose any practical problem.

The above identifying coefficient limiting range is given for illustrative purpose only, and may be equal to or may substantially approximate the identifying coefficient stable range, or may be of any shape insofar as most or all of the identifying coefficient limiting range belongs to the identifying coefficient stable range. Thus, the identifying coefficient limiting range may be established in various configurations in view of the ease with which to limit the values of the identified gain coefficients a1 hat, a2 hat and the practical controllability.

For example, while the boundary of an upper portion of the identifying coefficient limiting range is defined by the functional expression |a1|+a2=1 in the illustrated embodiment, combinations of the values of the gain coefficients a1, a2 which satisfy this functional expression are combinations of theoretical stable limits where a pole of the system defined by the equation (41) exists on a unit circle on a complex plane. Therefore, the boundary of the upper portion of the identifying coefficient limiting range may be determined by a functional expression |a1|+a2=r (r is a value slightly smaller than “1” corresponding to the stable limits, e.g., 0.99) for higher control stability.

The above identifying coefficient stable range shown in FIG. 12 as a basis for the identifying coefficient limiting range is given for illustrative purpose only. The identifying coefficient stable range which corresponds to the estimating coefficient stable range shown in FIG. 11 is affected by the dead time d (more precisely, its set value) and has its shape varied depending on the dead time d, as can be seen from the definition of the coefficient values α1, α2 (see the equation (12)). Irrespective of the shape of the identifying coefficient stable range, the identifying coefficient limiting range may be established, as described above, in a manner to match the shape of the identifying coefficient stable range.

In the present embodiment, the second limiting condition for limiting the value of the gain coefficient b1 identified by the identifier 25 , i.e., the value of the identified gain coefficient b1 hat, is established as follows:

The inventors have found that the situation in which the time-depending change of the target air-fuel ratio KCMD is oscillatory at a high frequency tends to happen also when the value of the identified gain coefficient b1 hat is excessively large or small. According to the present embodiment, an upper limit value B 1 H and a lower limit value B1L (B1H>B1L>0) for the identified gain coefficient b1 hat are determined in advance through experimentation or simulation. Then, the second limiting condition is established such that the identified gain coefficient b1 hat is equal to or smaller than the upper limit value B1H and equal to or greater than the lower limit value B1L (B1L≦b1 hat≦B1H).

A process of limiting the values of the identified gain coefficients a1 hat, a2 hat, b1 hat according to the first and second limiting conditions is carried out by in STEP 5 - 8 as follows:

As shown in FIG. 13 , the identifier 25 limits combinations of the identified gain coefficients a1(k) hat, a2(k) hat determined in STEP 5 - 7 shown in FIG. 10 according to the first limiting condition in STEP 5 - 8 - 1 through STEP 5 - 8 - 8 .

Specifically, the identifier 25 decides whether or not the value of the identified gain coefficient a2(k) hat determined in STEP 5 - 7 is equal to or greater than a lower limit value A2L (see FIG. 12 ) for the gain coefficient a2 in the identifying coefficient limiting range in STEP 5 - 8 - 1 .

If the value of the identified gain coefficient a2(k) is smaller than A2L, then since a point on the coordinate plane shown in FIG. 12 , which is expressed by (a1(k) hat, a2(k) hat), determined by the combination of the values of the identified gain coefficients a1(k) hat, a2(k) hat does not reside in the identifying coefficient limiting range, the value of a2(k) hat is forcibly changed to the lower limit value A2L in STEP 5 - 8 - 2 . Thus, the point (a1(k) hat, a2(k) hat) on the coordinate plane shown in FIG. 12 is limited to a point in a region on and above a straight line, i.e., the straight line including the line segment Q 7 Q 8 , expressed by at least a2=A2L.

Then, the identifier 25 decides whether or not the value of the identified gain coefficient a1(k) hat determined in STEP 5 - 7 is equal to or greater than a lower limit value A1L (see FIG. 12 ) for the gain coefficient a1 in the identifying coefficient limiting range in STEP 5 - 8 - 3 , and then decides whether or not the value of the identified gain coefficient a1(k) hat is equal to or smaller than an upper limit value A1H (see FIG. 12 ) for the gain coefficient a1 in the identifying coefficient limiting range in STEP 5 - 8 - 5 . The upper limit value A1H for the gain coefficient a1 in the identifying coefficient limiting range is represented by A1H=1−A2L because it is an al coordinate of the point Q 8 where the polygonal line |a1|+a2=1 (a1>0) and the straight line a2=A2L intersect with each other, as shown in FIG. 12 .

If the value of the identified gain coefficient a1(k) hat is smaller than the lower limit value A1L or greater than the upper limit value A1H, then since the point (a1(k) hat, a2(k) hat) on the coordinate plane shown in FIG. 12 does not reside in the identifying coefficient limiting range, the value of a1(k) hat is forcibly changed to the lower limit value A1L or the upper limit value A1H in STEP 5 - 8 - 4 , STEP 5 - 8 - 6 .

Thus, the point (a1(k) hat, a2(k) hat) on the coordinate plane shown in FIG. 12 is limited to a region on and between a straight line, i.e., the straight line including the line segment Q 6 Q 7 , expressed by a1=A1L, and a straight line, i.e., the straight line passing through the point Q 8 and perpendicular to the al axis, expressed by a1=A1H.

The processing in STEP 5 - 8 - 3 and STEP 5 - 8 - 4 and the processing in STEP 5 - 8 - 5 and STEP 5 - 8 - 6 may be switched around. The processing in STEP 5 - 8 - 1 and STEP 5 - 8 - 2 may be carried out after the processing in STEP 5 - 8 - 3 through STEP 5 - 8 - 6 .

Then, the identifier 25 decides whether the present values of a1(k) hat, a2(k) hat after STEP 5 - 8 - 1 through STEP 5 - 8 - 6 satisfy an inequality |a1+a2≦1 or not, i.e., whether the point (a1(k) hat, a2(k) hat) is positioned on or below or on or above the polygonal line (including line segments Q 5 Q 6 and Q 5 Q 8 ) expressed by the functional expression |a1+a2=1 in STEP 5 - 8 - 7 .

If |a1|+a2≦1, then the point (a1(k) hat, a2(k) hat) determined by the values of a1(k) hat, a2(k) hat after STEP 5 - 8 - 1 through STEP 5 - 8 - 6 exists in the identifying coefficient limiting range (including its boundaries).

If |a1|+a2>1, then since the point (a1(k) hat, a2(k) hat) deviates upwardly from the identifying coefficient limiting range, the value of the a2(k) hat is forcibly changed to a value (1−|a1(k) hat|) depending on the value of a1(k) hat in STEP 5 - 8 - 8 . Stated otherwise, while the value of a1(k) hat is being kept unchanged, the point (a1(k) hat, a2(k) hat) is moved onto a polygonal line expressed by the functional expression |a1|+a2=1, i.e., onto the line segment Q 5 Q 6 or the line segment Q 5 Q 8 which is a boundary of the identifying coefficient limiting range.

Through the above processing in STEP 5 - 8 - 1 through 5 - 8 - 8 , the values of the identified gain coefficients a1(k) hat, a2(k) hat are limited such that the point (a1(k) hat, a2(k) hat) determined thereby resides in the identifying coefficient limiting range. If the point (a1(k) hat, a2(k) hat) corresponding to the values of the identified gain coefficients a1(k) hat, a2(k) hat that have been determined in STEP 5 - 7 exists in the identifying coefficient limiting range, then those values of the identified gain coefficients a1(k) hat, a2(k) hat are maintained.

The value of the identified gain coefficient a1(k) hat relative to the primary autoregressive term of the discrete-system model is not forcibly changed insofar as the value resides between the lower limit value A1L and the upper limit value A1H of the identifying coefficient limiting range. If a1(k) hat<A1L or a1(k) hat >A1H, then since the value of the identified gain coefficient a1(k) hat is forcibly changed to the lower limit value A1L which is a minimum value that the gain coefficient a1 can take in the identifying coefficient limiting range or the upper limit value A1H which is a maximum value that the gain coefficient al can take in the identifying coefficient limiting range, the change in the value of the identified gain coefficient a1(k) hat is minimum. Stated otherwise, if the point (a1(k) hat, a2(k) hat) corresponding to the values of the identified gain coefficients a1(k) hat, a2(k) hat that have been determined in STEP 5 - 7 deviates from the identifying coefficient limiting range, then the forced change in the value of the identified gain coefficient a1(k) hat is held to a minimum.

After having limited the values of the identified gain coefficients a1(k) hat, a2(k) hat, the identifier 25 limits the identified gain coefficient b1(k) hat according to the second limiting condition in STEP 5 - 8 - 9 through STEP 5 - 8 - 12 .

Specifically, the identifier 25 decides whether or not the value of the identified gain coefficient b1(k) hat determined in STEP 5 - 7 is equal to or greater than the lower limit value B1L in STEP 5 - 8 - 9 . If the lower limit value B1L is greater than the value of the identified gain coefficient b1(k) hat, the value of b1(k) hat is forcibly changed to the lower limit value B1L in STEP 5 - 8 - 10 .

The identifier 25 decides whether or not the value of the identified gain coefficient b1(k) hat is equal to or smaller than the upper limit value B1H in STEP 5 - 8 - 11 . If the upper limit value B1H is smaller than the value of the identified gain coefficient b1(k) hat, the value of b1(k) hat is forcibly changed to the upper limit value B1H in STEP 5 - 8 - 12 .

If B1L≦b1(k) hat≦B1H, then the value of the identified gain coefficient b1(k) is maintained as it is.

Through the above processing in STEP 5 - 8 - 9 through 5 - 8 - 12 , the value of the identified gain coefficient b1(k) hat is limited to a range between the lower limit value B1L and the upper limit value B1H.

After the identifier 25 has limited the combination of the values of the identified gain coefficients a1(k) hat, a2(k) hat and the identified gain coefficient b1(k) hat, control returns to the sequence shown in FIG. 10 .

The preceding values a1(k−1) hat, a2(k−1) hat, b1(k−1) hat of the identified gain coefficients used for determining the identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat in STEP 5 - 7 shown in FIG. 10 are the values of the identified gain coefficients limited according to the first and second limiting conditions in STEP 5 - 8 in the preceding control cycle.

After having limited the identified gain coefficients a1(k) hat, a2(k) hat, b1(k) hat as described above, the identifier 25 updates the matrix P(k) according to the equation (10) for the processing of a next control cycle in STEP 5 - 9 , after which control returns to the main routine shown in FIG. 8 .

The processing subroutine of STEP 5 for the identifier 25 has been described above.

In FIG. 8 , after the processing of the identifier 25 has been carried out, the target air-fuel ratio generator 28 determines the gain coefficients a1, a2, b1 in STEP 6 .

In STEP 6 , if the value of the flag f/id/cal set in STEP 2 is “1”, i.e., if the gain coefficients a1, a2, b1 have been identified by the identifier 25 , then the values of the gain coefficients a1, a2, b1 are set to the identified gain coefficients a1 hat, a2 hat, b1 hat (limited in STEP 5 - 8 ) determined by the identifier 25 in STEP 5 . If the value of the flag f/id/cal is “0”, i.e., if the gain coefficients a1, a2, b1 have not been identified by the identifier 25 , then the values of the gain coefficients a1, a2, b1 are set to predetermined values.

Then, the target air-fuel ratio generator 28 effects a processing operation of the estimator 26 , i.e., calculates the estimated differential output VO2 bar, in STEP 7 of the main routine shown in FIG. 8 .

The calculating subroutine of STEP 7 is shown in detail in FIG. 14 . As shown in FIG. 14 , the estimator 26 calculates the coefficients α1, α2, βj (j=1−d) to be used in the equation (13), using the gain coefficients a1, a2, b1 determined in STEP 6 (these values are basically the identified gain coefficients a1 hat, a2 hat, b1 hat which have been limited in STEP 5 - 8 shown in FIG. 10 ) according to the definition accompanying the equation (12) in STEP 7 - 1 .

Then, in STEP 7 - 2 , the estimator 26 calculates the estimated differential output VO2(k+d) bar (estimated value of the differential output VO2 after the total dead time d from the time of the present control cycle) according to the equation (13), using the time-series data VO2(k), VO2(k−1), from before the present control cycle, of the differential output VO2 of the O 2 sensor calculated in each control cycle in STEP 3 shown in FIG. 8 , the time-series data kact(k−j) (j=0−d1), from before the present control cycle, of the differential output kact of the LAF sensor 5 , the time-series data kcmd(k−j) (normally, kcmd(k−j)=usl(k−j) (j=1−d2−1), from before the preceding control cycle, of the command differential air-fuel ratio kcmd (=the demand differential air-fuel ratio usl as limited) determined in each control cycle by the limiter 30 as described later on, and the coefficients α1, α2, βj calculated as described above.

Referring back to FIG. 8 , the target air-fuel ratio generator 28 then calculates the demand differential air-fuel ratio usl with the sliding mode controller 27 in STEP 8 .

The calculating subroutine of STEP 8 is shown in detail in FIG. 15 .

As shown in FIG. 15 , the sliding mode controller 27 calculates a value σ(k+d) bar (corresponding to an estimated value, after the dead time d, of the switching function σ defined according to the equation (15)), after the total dead time d from the present control cycle, of the switching function σ bar defined according to the equation (25), using the time-series data VO2(k+d) bar, VO2(k+d−1) bar (present and preceding values of the estimated differential output VO2 bar) of the estimated differential output VO2 bar determined by the estimator 26 in STEP 7 in STEP 8 - 1 .

If the switching function σ bar is excessively large, then the value of the reaching control law input urch determined depending on the value of the switching function σ bar becomes excessively large, causing the adaptive control law input uadp to change abruptly. Therefore, the demand differential air-fuel ratio usl determined by the sliding mode controller 27 and the target air-fuel ratio KCMD tends to become unstable. According to the present embodiment, the value of the switching function σ bar is set to fall in a predetermined range, and if the value of the σ bar determined according to the equation (25) exceeds an upper or lower limit of the predetermined range, then the value of the σ bar is forcibly set to the upper or lower limit of the predetermined range.

Then, the sliding mode controller 27 accumulates values of the switching function σ bar calculated in respective control cycles in STEP 8 - 1 (more accurately, values produced when the value of the σ bar is multiplied by the period (constant period) of the control cycles of the exhaust-side control unit 7 a ), i.e., adds a value of the σ bar calculated in the present control cycle to the sum determined in the preceding control cycle, thereby calculating an integrated value of the σ bar (which corresponds to the term at the right end of the equation (27)) in STEP 8 - 2 .

In order to prevent the adaptive control law input uadp determined depending on the integrated value of the σ bar, the integrated value of the σ bar is set to fall in a predetermined range, as with STEP 8 - 1 . Specifically, if the integrated value of the σ bar exceeds an upper or lower limit of the predetermined range, then the integrated value of the σ bar is limited to the upper or lower limit.

Then, the sliding mode controller 27 calculates the equivalent control input ueq, the reaching control law input urch, and the adaptive control law uadp according to the respective equations (24), (26), (27) in STEP 8 - 3 , using the time-series data VO2(k+d) bar, VO2(k+d−1) bar of the estimated differential output VO2 bar determined by the estimator 26 in STEP 7 , the value σ(k+d) bar of the switching function and its integrated value which are determined respectively in STEP 8 - 1 and STEP 8 - 2 , the gain coefficients a1, a1, b1 determined in STEP 6 (which are basically the gain coefficients a1 hat, a2 hat, b1 hat limited in STEP 5 - 8 shown in FIG. 10 ).

The sliding mode controller 27 then adds the equivalent control input ueq, the reaching control law input urch, and the adaptive control law uadp determined in STEP 8 - 3 to calculate an input to be applied to the object exhaust system E for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET in STEP 8 - 4 .

The processing operation of the sliding mode controller 27 in STEP 8 has been described above.

In FIG. 8 , the target air-fuel ratio generator 28 carries out a process with the limiter 30 . Prior to limiting the demand differential air-fuel ratio usl calculated by the sliding mode controller 27 , the limiter 30 determines the stability of the status of the output signal VO2/OUT of the O 2 sensor 6 (the output status of the object exhaust system, hereinafter referred to as an “SLD control status”) which is controlled according to the adaptive sliding mode control process carried out by the sliding mode controller 27 in STEP 9 .

A process of determining the stability will briefly be described prior to specifically describing the details of the process of determining the stability.

In the present embodiment, in each control cycle of the exhaust-side control unit 7 a , the sliding mode controller 27 uses a difference Δσ bar (corresponding to a rate of change of the value of the switching function σ bar) between the present value σ(k+d) bar of the switching function σ bar calculated in STEP 8 - 1 and a preceding value σ(k+d−1) bar thereof in a preceding control cycle, and the product σ(k+d) bar·Δσ bar of the switching function σ bar and the present value σ(k+d) bar, as a basic parameter for determining the stability of the SLD control status (the product σ(k+d) bar·Δσ bar will hereinafter be referred to as a “stability determining basic parameter Pstb”).

The stability determining basic parameter Pstb (=σ(k+d) bar·Δσ bar) corresponds to the time-differentiated function of a Lyapunov function σ bar 2 /2 relative to the switching function σ bar. The state in which Pstb≦0 is basically a state in which the value of the switching function σ bar is converged to or converging to “0” (the state quantity comprising the time-series data VO2(k+d) bar, VO2(k+d−1) bar of the estimated differential output VO2 bar is converged to or converging to the hyperplane σ bar=0). The state in which Pstb>0 is basically a state in which the value of the switching function σ bar is getting away from “0” (the state quantity comprising the time-series data VO2(k+d) bar, VO2(k+d−1) bar of the estimated differential output VO2 bar is getting away from the hyperplane σ bar=0).

Therefore, it is possible to determined whether the SLD control status is stable or not based on whether or not the value of the stability determining basic parameter Pstb is equal to or smaller than “0”.

However, if the stability of the SLD control status is determined solely by comparing the value of the stability determining basic parameter Pstb with “0”, then any slight noise contained in the value of the σ bar will affect the determined result of the stability. Furthermore, if the determined result of the stability is established based on the value of the stability determining basic parameter Pstb in each control cycle, then the determined result tends to change frequently.

In this embodiment, in each control cycle of the exhaust-side control unit 7 a , it is temporarily determined whether the SLD control status is stable or not based on whether or not the value of the stability determining basic parameter Pstb is equal or smaller than a predetermined value ε that is of a positive value slightly greater than “0”. In addition, in a predetermined period longer than the control cycles of the exhaust-side control unit 7 a , the frequency cnt/judst at which the SLD control status is temporarily determined as unstable based on the value of the stability determining basic parameter Pstb (more specifically, the number of control cycles in which Pstb>ε in the predetermined period, hereinafter referred to as a “temporarily unstable decision frequency cnt/judst”) is measured. Whether the SLD control status is stable or not is determined by comparing the temporarily unstable decision frequency cnt/judst with predetermined thresholds.

In the present embodiment, the predetermined thresholds with which the temporarily unstable decision frequency cnt/judst is compared include a first threshold SSTB 1 and a second threshold SSTB 2 (SSTB 1 <SSTB 2 ). If the temporarily unstable decision frequency cnt/judst is equal to or smaller than the first threshold SSTB 1 (cnt/judst≦SSTB 1 ), then the SLD control status is determined as stable. If cnt/judst>SSTB 1 , then the SLD control status is determined as unstable. If the SLD control status is determined as unstable (cnt/judst>SSTB 1 ), then when the temporarily unstable decision frequency cnt/judst is equal to or smaller than the second threshold SSTB 2 (cnt/judst≦SSTB 2 ), the level of the instability of the SLD control status is determined as low (such a state will hereinafter be referred to as a “low-level unstable state”), and when cnt/judst>SSTB 2 , the level of the instability of the SLD control status is determined as high (such a state will hereinafter be referred to as a “high-level unstable state”). Thus, if the SLD control status is determined as unstable, then the level of the instability of the SLD control status is also determined depending on the value of the temporarily unstable decision frequency cnt/judst.

Based on the above principles of the determination of the stability of the SLD control status, the process of determining the stability of the SLD control status will be described in greater detail below.

The process of determining the stability of the SLD control status is carried out according to a flowchart shown in FIG. 16 .

The limiter 30 calculates the stability determining basic parameter Pstb (=σ(k+d) bar·Δσ bar) defined above from the present value σ(k+d) and the preceding value σ(k+d−1) of the switching function σ bar determined by the sliding mode controller 27 in STEP 8 - 1 in STEP 9 - 1 in FIG. 16 .

Then, the limiter 30 determines whether the value of the timer counter tm/stb (count-down timer) initialized in STEPd- 11 in FIG. 7 has become “0” or not, i.e., whether the time that has elapsed from the start of the manipulation of the air-fuel ratio of the internal combustion engine 1 based on the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 (the control process for converging the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, which may hereinafter be referred to as an “air-fuel ratio manipulation exhaust system output control process”) has reached a predetermined time expressed by the initial value TMSTB of the timer counter tm/stb, in STEP 9 - 2 .

If tm/stb≠0, and hence the time that has elapsed from the start of the air-fuel ratio manipulation exhaust system output control process has not reached the predetermined time (: TMSTB) (immediately after start of the air-fuel ratio manipulation exhaust system output control process), then the stability of the SLD control status based on the measurement of the temporarily unstable decision frequency cnt/judst, etc. is not determined, and STEP 9 - 3 is carried out, after which control returns to the main routine shown in FIG. 8 .

In STEP 9 - 3 , the value of the timer counter tm/stb (count-down timer) for measuring the predetermined period in which to measure the temporarily unstable decision frequency cnt/judst is set to a predetermined initial value TMJUDST corresponding to the time of the predetermined period. In STEP 9 - 3 , the value of the temporarily unstable decision frequency cnt/judst is initialized to “0”. In STEP 9 - 3 , moreover, the value of a flag f/stb 1 which indicates whether the SLD control status is stable or not with “0” and “1”, respectively, is initialized to “0”, and the value of a flag f/stb 2 which indicates whether the SLD control status is in the low-level unstable state or the high-level unstable state with “0” and “1”, respectively, is initialized to “0”.

If the value of the timer tm/stb is “0”, and hence the time that has elapsed from the start of the air-fuel ratio manipulation exhaust system output control process has reached the predetermined time expressed by the initial value TMSTB of the timer tm/stb, then the limiter 30 compares the stability determining basic parameter Pstb with the predetermined value ε (>0) in STEP 9 - 4 . If Pstb≦ε, then the SLD control status is regarded temporarily stable, and the value of the temporarily unstable decision frequency cnt/judst is maintained at its present value (preceding value) in STEP 9 - 5 . If Pstb>ε, then the SLD control status is regarded as temporarily stable, and the value of the temporarily unstable decision frequency cnt/judst is incremented from the present value (preceding value) by “1” in STEP 9 - 6 .

Then, the limiter 30 compares the value of the temporarily unstable decision frequency cnt/judst determined in the present cycle in STEP 9 - 5 or STEP 9 - 6 with the first threshold SSTB 1 in STEP 9 - 7 . If cnt/judst≦SSTB 1 , then the SLD control status is regarded as stable, and the flag f/stb 1 (hereinafter referred to as a “stability decision flag f/stb 1 ”) is set to “0” in STEP 9 - 8 .

If cnt/judst>SSTB 1 (the SLD control status is unstable), then the limit 30 compares the value of the temporarily unstable decision frequency cnt/judst with the second threshold SSTB 2 in STEP 9 - 9 . If cnt/judst≦SSTB 2 , then the SLD control status is regarded as in the low-level unstable state, and the stability decision flag f/stb 1 is set to “1” in STEP 9 - 10 . If cnt/judst>SSTB 2 , then the SLD control status is regarded as in the high-level unstable state, and the stability decision flag f/stb 1 is set to “1” and the flag f/stb 2 (hereinafter referred to as an “instability level decision flag f/stb 2 ”) is set to “1” in STEP 9 - 11 .

Then, the limiter 30 decrements the value of the timer counter tm/judst from the present value (preceding value) by “1” in STEP 9 - 12 , and thereafter determines whether the updated value of the timer counter tm/judst has reached “0” or not, i.e., whether the predetermined period represented by the initial value TMJUDST of the timer counter tm/judst has elapsed or not, in STEP 9 - 13 .

If tm/judst≠0, and hence the predetermined period (: TMJUDST) has not elapsed (the timer counter tm/judst has not run out), then control returns to the main routine shown in FIG. 8 .

If tm/judst=0, and hence the predetermined period (: TMJUDST) has elapsed (the timer counter tm/judst has run out), then the limiter 30 determines the value of the stability decision flag f/stb 1 in STEP 9 - 14 . If f/stb 1 =1, then the value of the timer counter tm/judst is set to its initial value TMJUDST, and the value of the temporarily unstable decision frequency cnt/judst is reset to “0” in STEP 9 - 16 , after which control returns to the main routine shown in FIG. 8 .

If f/stb 1 =0 in STEP 9 - 14 , i.e., if the SLD control status in the present control cycle is determined as stable, then the value of the instability level decision flag f/stb 2 is reset to “0” in STEP 9 - 15 , and STEP 9 - 16 is executed, after which control returns to the main routine shown in FIG. 8 .

The manner in which the values of the stability determining basic parameter Pstb, the temporarily unstable decision frequency cnt/judst, the stability decision flag f/stb 1 , and the instability level decision flag f/stb 2 vary is illustrated respectively in first, second, third, and fourth stages of FIG. 17 . In the first stage of FIG. 17 , the value of the stability determining basic parameter Pstb is shown as being equal to or greater than “0”. However, the stability determining basic parameter Pstb may actually be of a negative value.

As shown in FIG. 17 , the temporarily unstable decision frequency cnt/judst, i.e., the frequency (the number of times) at which the value of the stability determining basic parameter Pstb is Pstb>ε, is measured in every predetermined period (: TMJUDST) and reset to “0” each time the predetermined period (: TMJUDST) elapses (see STEP 9 - 16 ). As indicated in a period T1 or T2 in FIG. 17 , if the value of the temporarily unstable decision frequency cnt/judst does not reach the first threshold SSTB 1 , i.e., if a situation in which the value of the stability determining basic parameter Pstb exceeds the predetermined value E does not occur or occurs only temporarily, then the SLD control status is regarded as stable, and the value of the stability decision flag f/stb 1 is set to “0” (see STEP 9 - 7 , STEP 9 - 8 ).

As indicated in a period T3 in FIG. 17 , if the situation in which the value of the stability determining basic parameter Pstb exceeds the predetermined value ε in the predetermined period continues for a certain time or occurs frequently, and the value of the temporarily unstable decision frequency cnt/judst exceeds the first threshold SSTB 1 , then the SLD control status is determined as in the low-level unstable state, and the value of the stability decision flag f/stb 1 is set to “1” (see STEP 9 - 9 , STEP 9 - 10 ). In this case, as the temporarily unstable decision frequency cnt/judst is reset in every predetermined period, the value of the stability decision flag f/stb 1 is also reset to “0” in every predetermined period.

As indicated in a period T4 in FIG. 17 , if the situation in which the value of the stability determining basic parameter Pstb exceeds the predetermined value ε in the predetermined period continues for a relatively long time or occurs frequently, and the value of the temporarily unstable decision frequency cnt/judst exceeds the second threshold SSTB 2 , then the SLD control status is determined as in the high-level unstable state, and the value of the stability decision flag f/stb 1 is set to “1” and the value of the instability level decision flag f/stb 2 is set to “1” (see STEP 9 - 11 ). In this case, the value of the instability level decision flag f/stb 2 is kept at “1” regardless of the value of the temporarily unstable decision frequency cnt/judst during a next predetermined period (period T5). Only if the value of the temporarily unstable decision frequency cnt/judst is continuously kept equal to or smaller than the first threshold SSTB 1 and the value of the stability decision flag f/stb 1 is kept at “0” during the next predetermined period (period T5), the instability level decision flag f/stb 2 is reset to “0” at the end of the period T5 (STEP 9 - 14 , STEP 9 - 15 ). Therefore, if the SLD control status is determined as in the high-level unstable state within a certain predetermined period, the value of the instability level decision flag f/stb 2 is kept at “1” indicative of the high-level unstable state within a next predetermined period unless the value of the temporarily unstable decision frequency cnt/judst is continuously kept equal to or smaller than the first threshold SSTB 1 , so that the decision of the high-level unstable state essentially continues.

Referring back to FIG. 8 , after the limiter 30 in the target air-fuel ratio generator 28 determines the stability of the SLD control status as described above, the limiter 30 limits the demand differential air-fuel ratio usl calculated by the sliding mode controller 27 in STEP 8 to determine the command differential air-fuel ratio kcmd in STEP 10 .

Prior to specifically describing the limiting process in detail, allowable ranges used by the limiting process will be described below.

The limiting process carried out by the limiter 30 forcibly sets the command differential air-fuel ratio kcmd to an upper or lower limit of an allowable range if the demand differential air-fuel ratio usl deviates from the allowable range beyond the upper or lower limit thereof, in order to keep the command differential air-fuel ratio kcmd which defines the target air-fuel ratio KCMD to be finally given to the engine-side control unit 7 b . If the demand differential air-fuel ratio usl falls in the allowable range, the limiting process sets up the demand differential air-fuel ratio usl directly as the command differential air-fuel ratio kcmd.

In the present embodiment, the limiting process uses a plurality of allowable ranges as shown in FIG. 18 . Specifically, the allowable ranges include an allowable range for unstable low level used to limit the demand differential air-fuel ratio usl when the determined stability of the SLD control status represents the low-level unstable state (f/stb 1 =1 and f/stb 2 =0) (except when the internal combustion engine 1 is idling) in STEP 9 , and an allowable range for unstable high level used to limit the demand differential air-fuel ratio usl when the determined stability of the SLD control status represents the high-level unstable state (f/stb 2 =1).

The allowable ranges for the limiting process also include an allowable range after FC (fuel cut-off) used to limit the demand differential air-fuel ratio usl immediately after the fuel supply to the internal combustion engine 1 is cut off (specifically, until a predetermined time elapses after the fuel supply to the internal combustion engine 1 is cut off) except when the determined stability of the SLD control status represents the low-level unstable state or the high-level unstable state, an allowable range after start used to limit the demand differential air-fuel ratio usl immediately after the internal combustion engine 1 starts (more specifically, until a predetermined time elapses after the internal combustion engine 1 starts), and an allowable range after operation with lean mixture used to limit the demand differential air-fuel ratio usl immediately after the operation mode of the internal combustion engine 1 changes from the lean operation mode in which the air-fuel ratio of the internal combustion engine 1 is a lean ratio to the normal operation mode in which the air-fuel ratio manipulation exhaust system output control process is performed (specifically, until a predetermined time elapses after the lean operation mode is ended). In this embodiment, the allowable range after FC, the allowable range after start, and the allowable range after operation with lean mixture are identical allowable ranges, and will collectively be referred to as an allowable range after FC/start/operation with lean mixture.

The allowable ranges for limiting the demand differential air-fuel ratio usl further include an allowable range after being driven under load used to limit the demand differential air-fuel ratio usl immediately after a vehicle incorporating the internal combustion engine 1 starts (specifically, until a predetermined time elapses after the internal combustion engine 1 starts to operate drive wheels of the vehicle which serves as a load on the internal combustion engine 1 ), except when the determined stability of the SLD control status represents the low-level unstable state or the high-level unstable state, and except immediately after the fuel supply to the internal combustion engine 1 is cut off, or the internal combustion engine 1 starts, or the internal combustion engine 1 operates with lean mixture, and an allowable range idling used to limit the demand differential air-fuel ratio usl while the internal combustion engine 1 is idling, except when the determined stability of the SLD control status represents the high-level unstable state, and except immediately after the fuel supply to the internal combustion engine 1 is cut off, or the internal combustion engine 1 starts, or the internal combustion engine 1 operates with lean mixture.

The allowable ranges for limiting the demand differential air-fuel ratio usl also include an adaptive allowable range used to limit the demand differential air-fuel ratio usl while the internal combustion engine 1 is operates normally except those states corresponding to the above various allowable ranges.

The allowable range for unstable low level is a fixed allowable range with its upper and lower limits set respectively to predetermined fixed values H, L (the predetermined value H will hereinafter be referred to as an “upper-limit third predetermined value H” and the predetermined value L will hereinafter be referred to as a “lower-limit fourth predetermined value L”). The allowable range for unstable low level has a standard (intermediate) extent at the upper and lower limits thereof, among the various allowable ranges.

The allowable range for unstable high level is a fixed allowable range with its upper and lower limits set respectively to predetermined fixed values STABH, STABL (the predetermined value STABH will hereinafter be referred to as an “upper-limit first predetermined value STABH” and the predetermined value STABL will hereinafter be referred to as a “lower-limit first predetermined value STABL”). The allowable range for unstable high level has a narrowest extent at the upper and lower limits thereof, among the various allowable ranges.

The allowable range for idling is a fixed allowable range with its upper and lower limits set respectively to predetermined fixed values HI, LI (the predetermined value HI will hereinafter be referred to as an “upper-limit second predetermined value HI” and the predetermined value LI will hereinafter be referred to as a “lower-limit third predetermined value LI”). The allowable range for idling has a relatively narrow extent at the upper and lower limits thereof (narrower than the allowable range for unstable low level).

The allowable range after FC/start/operation with lean mixture is an allowable range with its upper limit set to a predetermined fixed value AFCH (hereinafter referred to as an “upper-limit fifth predetermined value AFCH”) and its lower limit sequentially variable (in each control cycle) depending on how the demand differential air-fuel ratio usl deviates from the allowable range, between the lower-limit first predetermined value STABL and a predetermined value LL smaller than the lower-limit first predetermined value STABL (hereinafter referred to as a “lower-limit fifth predetermined value LL”). The lower limit of the allowable range after FC/start/operation with lean mixture is widest among the various allowable ranges. While the internal combustion engine 1 is idling, the lower limit of the allowable range after FC/start/operation with lean mixture is not variable, but set to the lower-limit third predetermined value LI which is the lower limit of the allowable range for idling.

The allowable range after being driven under load is an allowable range with its lower limit set to a predetermined fixed value VSTL (hereinafter referred to as a “lower-limit second predetermined value VSTL”) and its upper limit sequentially variable (in each control cycle) depending on how the demand differential air-fuel ratio usl deviates from the allowable range, between the upper-limit first predetermined value STABH and a predetermined value HH larger than the upper-limit first predetermined value STABH (hereinafter referred to as an “upper-limit fourth predetermined value HH”). The lower limit of the allowable range after being driven under load is relatively narrow (narrower than the lower limit of the allowable range for idling).

The adaptive allowable range is an allowable range with its upper limit sequentially variable (in each control cycle) depending on how the demand differential air-fuel ratio usl deviates from the allowable range, between the upper-limit first predetermined value STABH and the upper-limit fourth predetermined value HH and its lower limit sequentially variable (in each control cycle) depending on how the demand differential air-fuel ratio usl deviates from the allowable range, between the lower-limit first predetermined value STABL and the lower-limit fifth predetermined value LL.

Since the demand differential air-fuel ratio usl to be limited to the various allowable ranges represents the difference between the air-fuel ratio of the internal combustion engine 1 and the air-fuel ratio reference value FLAF/BASE, the upper and lower limits of the allowable ranges represent differences with the air-fuel ratio reference value FLAF/BASE, i.e., positive and negative values with respect to the air-fuel ratio reference value FLAF/BASE represented as “0”. The upper-limit first through fifth predetermined values STABH, HI, H, HH, AFCH are related in magnitude with each other as 0<STABH<HI<H<HH<AFCH as shown in FIG. 18 . Similarly, the lower-limit first through fifth predetermined values STABL, VSTL, LI, L, LL are related in magnitude with each other as 0>STABL>VSTL>LI>L>LL as shown in FIG. 18 .

In view of the above definition of the various allowable ranges, the limiting process in STEP 10 will be described in specific detail below. The limiting process is carried out according to a flowchart shown in FIG. 19 . An allowable range for the limiting process which is established in an idling state of the internal combustion engine 1 , which allowable range is not limited to the above allowable range for idling, has an upper limit AHFI and a lower limit ALFI which will be referred to respectively as an idling upper limit AHFI and an idling lower limit ALFI. An allowable range for the limiting process which is established in an operating state of the internal combustion engine 1 except for the idling state has an upper limit AHF and a lower limit ALF which will be referred to respectively as a non-idling upper limit AHF and a non-idling lower limit ALF. The variable upper and lower limits of the adaptive allowable range will be referred to respectively as an adaptive upper limit ah and an adaptive lower limit al.

As shown in FIG. 19 , the limiter 30 carries out a process of determining an allowable range for limiting the demand differential air-fuel ratio usl in the present cycle in STEP 10 - 1 .

The process of determining an allowable range in STEP 10 - 1 is performed according to a flowchart shown in FIG. 20 .

First, the limiter 30 determines the value of the instability level decision flag f/stb 2 established in the process of determining the stability of the SLD control status in STEP 10 - 1 - 1 . If f/stb 2 =1, i.e., if the SLD control status is in the high-level unstable state, then the limiter 30 sets the idling upper limit AHFI, the non-idling upper limit AHF, and the adaptive upper limit ah (more specifically, a present value ah(k−1) of the adaptive upper limit ah to the upper-limit first predetermined value STABH which is the upper limit of the allowable range for unstable high level (see FIG. 18 ), and also sets the idling lower limit ALFI, the non-idling lower limit ALF, and the adaptive lower limit al (more specifically, a present value al(k−1) of the adaptive lower limit al) to the lower-limit first predetermined value STABL which is the lower limit of the allowable range for unstable high level in STEP 10 - 1 - 2 . After STEP 10 - 1 - 2 , control returns to the routine shown in FIG. 19 .

In STEP 10 - 1 - 2 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range for unstable high level, i.e., the narrowest fixed allowable range, regardless of the operating state of the internal combustion engine 1 . The adaptive allowable range which is made variable as described later on is initialized to the allowable range for unstable high level.

If f/stb 2 =0 in STEP 10 - 1 - 1 , then the limiter 30 determines the stability decision flag f/stb 1 established in the process of determining the stability of the SLD control status in STEP 10 - 1 - 3 . If f/stb 1 =1, i.e., if the SLD control status is in the low-level unstable state, then the limiter 30 sets the idling upper limit AHFI to the upper-limit second predetermined value HI which is the upper limit of the idling allowable range (see FIG. 18 ) and also sets the idling lower limit ALFI to the lower-limit third predetermined value LI which is the lower limit of the idling allowable range in STEP 10 - 1 - 4 . In STEP 10 - 1 - 4 , the limiter 30 also sets the non-idling upper limit AHF and the present value ah(k−1) of the adaptive upper limit ah to the upper-limit third predetermined value H which is the fixed upper limit of the allowable range for unstable low level (see FIG. 18 ), and also sets the non-idling lower limit ALF and the adaptive lower limit al (al(k−1)) to the lower-limit fourth predetermined value L which is the fixed lower limit of the allowable range for unstable low level. After STEP 10 - 1 - 4 , control returns to the routine shown in FIG. 19 .

In STEP 10 - 1 - 4 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range for unstable low level (the fixed standard range) in an operating state of the internal combustion engine 1 except for the idling state, and set to the allowable range for idling (fixed narrow range) which is narrower than the allowable range for unstable low level in the idling state of the internal combustion engine 1 . The adaptive allowable range is initialized to the allowable range for unstable low level.

If f/stb 1 =0 in STEP 10 - 1 - 3 , i.e., if the SLD control status is stable in STEP 9 , then the limiter 30 determines whether the internal combustion engine 1 is in a state immediately after the supply of fuel is cut off or not, i.e., whether the time that has elapsed after the end of the cutting-off of the supply of fuel has not reached a predetermined time or not, in STEP 10 - 1 - 5 .

The decision in STEP 10 - 1 - 5 is carried out as follows: In this embodiment, data indicative of whether the supply of fuel to the internal combustion engine 1 is cut off or not is given from the engine-side control unit 7 b to the exhaust-side control unit 7 a . The limiter 30 activates a timer (not shown) from the time when it recognizes the end of the cutting-off of the supply of fuel based on the data, for thereby measuring the time that has elapsed from the time of the end of the cutting-off of the supply of fuel. The limiter 30 determines a period until the elapsed time reaches a predetermined time (fixed value) as the state immediately after the supply of fuel is cut off.

In this embodiment, while the supply of fuel is being cut off, the manipulation of the air-fuel ratio based on the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 , i.e., the control process of converging the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET, is not performed (see FIG. 7 showing details of STEPd in FIG. 6 ), but the processing in the target air-fuel ratio generator 28 for calculating the demand differential air-fuel ratio usl and the target air-fuel ratio KCMD is carried out. In the present embodiment, the state in which the supply of fuel to the internal combustion engine 1 is being cut off is determined as the state immediately after the supply of fuel is cut off.

If the internal combustion engine 1 is in the state immediately after the supply of fuel is cut off (including the state in which the supply of fuel to the internal combustion engine 1 is being cut off) in STEP 10 - 1 - 5 , then the limiter 30 sets the idling upper limit AHFI and a non-idling upper limit AHF to the upper-limit fifth predetermined value AFCH which is the upper limit of the allowable range after FC/start/operation with lean mixture (see FIG. 18 ) in STEP 10 - 1 - 6 . In STEP 10 - 1 - 6 , the limiter 30 also sets the idling lower limit ALFI to the lower-limit third predetermined value LI which is the lower limit of the allowable range for idling and sets the non-idling lower limit ALF to the present value of the adaptive lower limit al (the adaptive lower limit al(k−1) determined in the preceding control cycle). After STEP 10 - 1 - 6 , control goes back to the routine shown in FIG. 19 .

In STEP 10 - 1 - 6 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after FC/start/operation with lean mixture whose lower limit is fixed to the lower-limit third predetermined value LI which is the lower limit of the allowable range for idling if the internal combustion engine 1 is in the idling state. If the internal combustion engine 1 is in an operating state other than the idling state, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after FC/start/operation with lean mixture whose lower limit is set to the sequentially variable adaptive lower limit al. At any rate, the allowable range for limiting the demand differential air-fuel ratio usl is set to a range which is widest at its upper limit, or more generally, a range whose upper limit is not lower than the upper-limit fifth predetermined value AFCH.

If the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off in STEP 10 - 1 - 5 , then the limiter 30 determines whether the internal combustion engine 1 is in a state immediately after it starts or not, i.e., whether the time that has elapsed after the start of the internal combustion engine 1 (more specifically, after the confirmation of the full fuel combustion in the internal combustion engine 1 ) does not reach a predetermined time or not, in STEP 10 - 1 - 7 .

In this case, data indicative of whether the internal combustion engine 1 starts or not, i.e., whether the full fuel combustion in the internal combustion engine 1 is confirmed or not, is given from the engine-side control unit 7 b to the exhaust-side control unit 7 a . The limiter 30 activates a timer (not shown) from the time when it recognizes the start of the internal combustion engine 1 based on the data, for thereby measuring the time that has elapsed from the start of the internal combustion engine 1 . The limiter 30 determines a period until the elapsed time reaches a predetermined time (fixed value) as the state immediately after the internal combustion engine 1 starts.

If the internal combustion engine 1 is in the state immediately after it starts in STEP 10 - 1 - 7 , then the limiter 30 carries out STEP 10 - 1 - 6 to set the idling upper limit AHFI, the non-idling upper limit AHF, the idling lower limit ALFI, and the non-idling lower limit ALF as described above, after which control returns to the routine shown in FIG. 19 .

If the internal combustion engine 1 is not in the state immediately after it starts in STEP 10 - 1 - 7 , then the limiter 30 determines whether the internal combustion engine 1 is in a state immediately after the end of its lean operation mode or not, i.e., whether the time that has elapsed after the operation mode of the internal combustion engine 1 changes from the lean operation mode to the normal operation mode has reached a predetermined time or not, in STEP 10 - 1 - 8 .

In this case, data indicative of whether the internal combustion engine 1 is in the lean operation mode or not is given from the engine-side control unit 7 b to the exhaust-side control unit 7 a . The limiter 30 activates a timer (not shown) from the time when it recognizes the end of the lean operation mode based on the data, for thereby measuring the time that has elapsed from the start of the timer. The limiter 30 determines a period until the elapsed time reaches a predetermined time (fixed value) as the state immediately after the end of its lean operation mode, more specifically, the state immediately after the operation mode of the internal combustion engine 1 changes from the lean operation mode to the normal operation mode.

If the internal combustion engine 1 is in the state immediately after the end of its lean operation mode in STEP 10 - 1 - 8 , then the limiter 30 carries out STEP 10 - 1 - 6 to set the idling upper limit AHFI, the non-idling upper limit AHF, the idling lower limit ALFI, and the non-idling lower limit ALF as described above, after which control returns to the routine shown in FIG. 19 .

If the internal combustion engine 1 is not in the state immediately after the end of its lean operation mode in STEP 10 - 1 - 8 , then the limiter 30 determines whether the internal combustion engine 1 is in its idling state or not in STEP 10 - 1 - 9 .

In this case, data indicative of whether the internal combustion engine 1 is in its idling state or not is given from the engine-side control unit 7 b to the exhaust-side control unit 7 a . The limiter 30 determines whether the internal combustion engine 1 is in its idling state or not based on the data.

If the internal combustion engine 1 is in the idling state (at this time, the SLD control status is not in the low-level unstable state or the high-level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode), then the limiter 30 sets the idling upper limit AHFI to the upper-limit second predetermined value HI which is the upper limit of the allowable range for idling, and also sets the idling lower limit ALFI to the lower-limit third predetermined value LI which is the lower limit of the allowable range for idling in STEP 10 - 1 - 10 . Thereafter, control returns to the routine shown in FIG. 19 .

In STEP 10 - 1 - 10 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the fixed allowable range for idling.

If the internal combustion engine 1 is not in the idling state (at this time, the SLD control status is not in the low-level unstable state or the high-level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, or the internal combustion engine 1 is not idling), then the limiter 30 performs a process of limiting the present value ah(k−1) (the value determined in the preceding control cycle) of the adaptive upper limit ah to a value in a range between the upper-limit first predetermined value STABH and the upper-limit fourth predetermined value HH (see FIG. 18 ), i.e., a process of limiting the adaptive upper limit ah, in STEP 10 - 1 - 11 through STEP 10 - 1 - 14 described below.

If the present value ah(k−1) of the adaptive upper limit ah is smaller than the upper-limit first predetermined value STABH in STEP 10 - 1 - 11 , then the limiter 30 sets the present value ah(k−1) forcibly to the upper-limit first predetermined value STABH in STEP 10 - 1 - 12 . If the present value ah(k−1) of the adaptive upper limit ah is greater than the upper-limit fourth predetermined value HH in STEP 10 - 1 - 13 , then the limiter 30 sets the present value ah(k−1) forcibly to the upper-limit fourth predetermined value HH in STEP 10 - 1 - 14 . Otherwise, the limiter 30 maintains the present value ah(k−1) of the adaptive upper limit ah.

Similarly, the limiter 30 perform a process of limiting the present value al(k−1) of the adaptive lower limit al to a value in a range between the lower-limit first predetermined value STABL and the lower-limit fifth predetermined value LL (see FIG. 18 ), i.e., a process of limiting the adaptive lower limit al in STEP 10 - 1 - 15 through STEP 10 - 1 - 18 . The values a1, STABL, LL are negative values.

After having limited the ranges for the present upper value ah(k−1) and the present lower value al(k−1), the limiter 30 sets the non-idling upper limit AHF and the non-idling lower limit ALF respectively to the adaptive upper limit ah(k−1) and the adaptive lower limit al(k−1) in STEP 10 - 1 - 19 .

The limiter 30 then determines whether the vehicle is in a state immediately after its start or not, i.e., whether the time that has elapsed after the internal combustion engine 1 starts to operate drive wheels of the vehicle which serves as a load on the internal combustion engine 1 has not reached a predetermined time or not, in STEP 10 - 1 - 20 .

The decision in STEP 10 - 1 - 20 is carried out as follows: In the present embodiment, the data indicative of whether the internal combustion engine 1 is in the idling state and also data indicative of the speed of the vehicle are given from the engine-side control unit 7 b to the exhaust-side control unit 7 a . Based on these data, the limiter 30 recognizes a state in which the internal combustion engine 1 is idling and the vehicle speed is substantially “0”, i.e., a parked or stopped state of the vehicle, and recognizes a start of the vehicle when the vehicle speed exceeds a predetermined value (sufficiently small value) from the parked or stopped state of the vehicle. The limiter 30 measures the time that has elapsed from the recognized start of the vehicle with a timer (not shown), and recognizes a period until the elapsed time reaches a predetermined time as the state immediately after the start of the vehicle.

If the vehicle is in the state immediately after its start in STEP 10 - 1 - 20 , then the limiter 30 sets the non-idling lower limit ALF to the lower-limit second predetermined value VSTL which is the lower limit of the allowable range after being driven under load (see FIG. 18 ) in STEP 10 - 1 - 21 , after which control returns to the routine shown in FIG. 19 .

If the vehicle is not in the state immediately after its start in STEP 10 - 1 - 20 , then the limiter 30 maintains the non-idling upper limit AHF and the non-idling lower limit ALF established in STEP 10 - 1 - 19 , after which control returns to the routine shown in FIG. 19 .

In STEP 10 - 1 - 19 through STEP 10 - 1 - 21 , if the SLD control status is not in the low-level unstable state or the high-level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, or the internal combustion engine 1 is not idling, and if the vehicle is not in the state immediately after its start (in a normal case), the allowable range for limiting the demand differential air-fuel ratio usl is set to the adaptive allowable range whose upper limit (adaptive upper limit ah) and lower limit (adaptive lower limit al) are sequentially varied, as described later on. If the vehicle is in the state immediately after its start, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after being driven under load, i.e., a range whose lower limit is relative narrow, or more generally, a range whose lower limit is not smaller than the lower-limit second predetermined value VSTL.

Referring back to FIG. 19 , after the allowable range for limiting the demand differential air-fuel ratio usl has been determined, the limiter 30 determines whether the internal combustion engine 1 is in its idling state or not in STEP 10 - 2 .

If the internal combustion engine 1 is in its idling state, then the limiter 30 limits the demand differential air-fuel ratio usl to the allowable range (normally the allowable range for idling) determined by the idling upper limit AHFI and the idling lower limit ALFI which have been set in either one of STEP 10 - 1 - 2 , STEP 10 - 1 - 4 , STEP 10 - 1 - 6 , and STEP 10 - 1 - 10 (normally STEP 10 - 1 - 10 ), for thereby determining the command differential air-fuel ratio kcmd in STEP 10 - 3 through STEP 10 - 7 .

Specifically, if the demand differential air-fuel ratio usl generated by the sliding mode controller 27 in STEP 8 (see FIG. 8 ) deviates from the allowable range beyond the idling lower limit ALFI (negative) (usl<ALFI), then the limiter 30 limits the value of the demand differential air-fuel ratio usl to the idling lower limit ALFI in STEP 10 - 3 and STEP 10 - 4 . If the demand differential air-fuel ratio usl deviates from the allowable range beyond the idling upper limit AHFI (positive) (usl>AHFI), then the limiter 30 limits the value of the demand differential air-fuel ratio usl to the idling upper limit AHFI in STEP 10 - 5 and STEP 10 - 6 . If the demand differential air-fuel ratio usl falls in the allowable range (ALFI≦usl≦AHFI), then the limiter 30 sets the demand differential air-fuel ratio usl as the command differential air-fuel ratio kcmd in STEP 10 - 7 .

If the demand differential air-fuel ratio usl deviates from the allowable range between the idling upper limit AHFI and the idling lower limit ALFI in STEP 10 - 3 or STEP 10 - 5 , then in order to prevent the value of the adaptive control law input uadp determined while the sliding mode controller 27 is generating the demand differential air-fuel ratio usl from becoming larger than necessary, the limiter 30 sets the integrated value of the switching function σ bar (see STEP 8 - 2 in FIG. 15 ) forcibly to the integrated value determined in the preceding control cycle in STEP 10 - 8 . If the limiter 30 limits the demand differential air-fuel ratio usl while the internal combustion engine 1 is idling, the present values of the adaptive upper and lower limits ah, al are maintained in STEP 10 - 9 . After STEP 10 - 3 through STEP 10 - 9 , control returns to the main routine shown in FIG. 8 .

If the internal combustion engine 1 is not in its idling state, then the limiter 30 limits the demand differential air-fuel ratio usl to the allowable range (normally the adaptive allowable range) determined by the non-idling upper limit AHF and the non-idling lower limit ALF which have been set in either one of STEP 10 - 1 - 2 , STEP 10 - 1 - 4 , STEP 10 - 1 - 6 , STEP 10 - 1 - 19 , and STEP 10 - 1 - 21 (normally STEP 10 - 1 - 19 ), for thereby determining the command differential air-fuel ratio kcmd in STEP 10 - 10 through STEP 10 - 14 .

Specifically, if the demand differential air-fuel ratio usl generated by the sliding mode controller 27 in STEP 8 (see FIG. 8 ) deviates from the allowable range beyond the non-idling lower limit ALF (negative) (usl<ALF), then the limiter 30 limits the value of the demand differential air-fuel ratio usl to the non-idling lower limit ALF in STEP 10 - 10 and STEP 10 - 11 . If the demand differential air-fuel ratio usl deviates from the allowable range beyond the non-idling upper limit AHF (positive) (usl>AHF), then the limiter 30 limits the value of the demand differential air-fuel ratio usl to the non-idling upper limit AHF in STEP 10 - 12 and STEP 10 - 13 . If the demand differential air-fuel ratio usl falls in the allowable range (ALF≦usl≦AHF), then the limiter 30 sets the demand differential air-fuel ratio usl as the command differential air-fuel ratio kcmd in STEP 10 - 14 .

In case the demand differential air-fuel ratio usl is limited by the allowable range determined by the non-idling upper limit AHF and the non-idling lower limit ALF, the limiter 30 then updates (changes) the values of the adaptive upper and lower limits ah, al relative to the adaptive allowable range in STEP 10 - 15 through STEP 10 - 19 .

Specifically, if the demand differential air-fuel ratio usl falls in the allowable range in STEP 10 - 10 and STEP 10 - 12 , i.e., if STEP 10 - 14 is executed, then the limiter 30 adds a predetermined change ΔDEC (>0, hereinafter referred to as a “decreasing unit change ΔDEC”) to the present value al(k−1) of the adaptive lower limit al for thereby determining a new adaptive lower limit al(k) in STEP 10 - 15 . In STEP 10 - 15 , the limiter 30 also subtracts the decreasing unit change ΔDEC from the present value ah(k−1) of the adaptive upper limit ah for thereby determining a new adaptive upper limit ah(k).

If the demand differential air-fuel ratio usl deviates from the allowable range beyond the non-idling lower limit ALF (negative), i.e., if STEP 10 - 11 is executed, then the limiter 30 determines whether the vehicle is in a state immediately after its start or not in STEP 10 - 16 . If the vehicle is not in the state immediately after its start, then the limiter 30 subtracts a predetermined change ΔINC (>0, hereinafter referred to as an “increasing unit change ΔINC”) from the present value al(k−1) of the adaptive lower limit al for thereby determining a new adaptive lower limit al(k) in STEP 10 - 17 . In STEP 10 - 17 , the limiter 30 also subtracts the decreasing unit change ΔDEC from the present value ah(k−1) of the adaptive upper limit ah for thereby determining a new adaptive upper limit ah(k).

If the vehicle is in the state immediately after its start in STEP 10 - 16 , then the present values of the adaptive upper and lower limits ah, al are maintained in STEP 10 - 18 .

In the present embodiment, the decreasing unit change ΔDEC and the increasing unit change ΔINC are related to each other such that ΔDEC<ΔINC.

If the demand differential air-fuel ratio usl deviates from the allowable range beyond the non-idling upper limit AHF (positive), then the limiter 30 determines whether the internal combustion engine 1 is in the state immediately after the supply of fuel is cut off or not (including while the supply of fuel is being cut off or not), or whether the internal combustion engine 1 is in the state immediately after it starts or not, or whether the internal combustion engine 1 is in the state immediately after the end of its lean operation mode or not in STEP 10 - 19 . If the internal combustion engine 1 is not in any of these states, then the limiter 30 adds the increasing unit change ΔINC to the present value ah(k−1) of the adaptive upper limit ah for thereby determining a new adaptive upper limit ah(k) in STEP 10 - 20 . In STEP 10 - 20 , the limiter 30 also adds the decreasing unit change ΔDEC to the present value al(k−1) of the adaptive lower limit al for thereby determining a new adaptive lower limit al(k).

If the internal combustion engine 1 is in either one of the states in STEP 10 - 19 , then the present values of the adaptive upper and lower limits ah, al are maintained in STEP 10 - 21 .

If the demand differential air-fuel ratio usl deviates from the allowable range determined by the non-idling upper limit AHF and the non-idling lower limit ALF in STEP 10 - 10 or STEP 10 - 12 , then the limiter 30 sets the integrated value of the switching function σ bar (see STEP 8 - 2 in FIG. 15 ) forcibly to the integrated value determined in the preceding control cycle in STEP 10 - 22 . After STEP 10 - 22 , control returns to the main routine shown in FIG. 8 .

In the present embodiment, the initial values of the adaptive upper and lower limits ah, al when the exhaust-side control unit 7 a is activated, i.e., when the vehicle starts operating, are respectively the lower-limit fourth predetermined value L and the upper-limit third predetermined value H. That is, the initial range of the adaptive allowable range is the same as the allowable range for unstable low level.

In the present embodiment, if the SLD control status is in the low-level or high-level unstable state when the internal combustion engine 1 is in an operating state other than the idling state, the adaptive upper and lower limits ah, al are updated in STEP 10 - 15 , STEP 10 - 17 , STEP 10 - 20 . In the above unstable state, the adaptive upper and lower limits ah, al are forcibly set to fixed predetermined values, i.e., the upper-limit first predetermined value STABH, etc. corresponding to the allowable range for unstable low level and the allowable range for unstable high level, by the processing of STEP 10 - 1 - 2 or STEP 10 - 1 - 4 (see FIG. 20 ) in a next control cycle. Therefore, if the SLD control status is in the low-level or high-level unstable state, then the processing of STEP 10 - 15 , STEP 10 - 17 , STEP 10 - 20 may be omitted.

The command differential air-fuel ratio kcmd determined in each control cycle by the limiting process in STEP 10 is stored in a time-series fashion in a memory (not shown), and will be used for the above processing operation of the estimator 26 .

According to the limiting process in STEP 10 , while the internal combustion engine 1 is in the idling state, either one of the fixed allowable ranges including the allowable range for unstable low level, the allowable range for unstable high level, the allowable range after FC/start/operation with lean mixture (specifically, the allowable range after FC/start/operation with lean mixture whose lower limit is fixed to the lower-limit third predetermined value LI (see FIG. 18 )), and the allowable range for idling, is established as the allowable range for limiting the demand differential air-fuel ratio usl (normally the allowable range for idling is established) depending on the SLD control status and the operating state of the internal combustion engine 1 . The demand differential air-fuel ratio usl is limited by the established allowable range to determine the command differential air-fuel ratio kcmd.

When the internal combustion engine 1 is in an operating state other than the idling state, either one of the allowable ranges including the allowable range for unstable low level, the allowable range for unstable high level, the allowable range after FC/start/operation with lean mixture (specifically, the allowable range after FC/start/operation with lean mixture whose lower limit is set to the variable adaptive lower limit al), the allowable range after being driven under load, and the adaptive allowable range is established as the allowable range for limiting the demand differential air-fuel ratio usl (normally the adaptive allowable range is established) depending on the SLD control status and the operating state of the internal combustion engine 1 . The demand differential air-fuel ratio usl is limited by the established allowable range to determine the command differential air-fuel ratio kcmd.

Regardless of whether the internal combustion engine 1 is in an operating state other than the idling state or in the idling state, when the determined stability of the SLD control status is in the high-level unstable state, the allowable range for limiting the demand differential air-fuel ratio usl is set to the narrowest allowable range for unstable high level whose upper limit AHF or AHFI and lower limit ALF or ALFI are set respectively to the fixed upper-limit first predetermined value STABH and the fixed lower-limit first predetermined value STABL in STEP 10 - 1 - 2 .

If the SLD control status is in the low-level unstable state and the internal combustion engine 1 is in an operating state other than the idling state, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range for unstable low level, which has a standard extent, whose upper and lower limits AHF, ALF are set respectively to the upper-limit third predetermined value H and the lower-limit fourth predetermined value L in STEP 10 - 1 - 4 .

If the SLD control status is not in the high-level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, or the vehicle is not in the state immediately after it starts, and the internal combustion engine 1 is in the idling state, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the relatively narrow allowable range for idling whose upper and lower limits AHFI, ALFI are set respectively to the fixed upper-limit second predetermined value HI and the fixed lower-limit third predetermined value LI in STEP 10 - 1 - 4 , STEP 10 - 1 - 10 .

If in the normal state, i.e., if the SLD control status is not in the high-level unstable state or in the low level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, or not in the idl-ing state, or the vehicle is not in the state immediately after it starts, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the adaptive allowable range whose upper and lower limits AHF, ALF are set respectively to the adaptive upper and lower limits ah, al in STEP 10 - 1 - 19 .

In the adaptive allowable range, if the demand differential air-fuel ratio usl is in the adaptive allowable range (at this time, kcmd=usl), then after the limiting process, the adaptive upper and lower limits ah, al of the adaptive allowable range are updated in such a direction that they decrease by the decreasing unit change ΔDEC in each control cycle to respective limits provided by the upper-limit first predetermined value STABH and the lower-limit first predetermined value STABL in STEP 10 - 15 .

If the demand differential air-fuel ratio usl deviates from the adaptive allowable range beyond the upper limit thereof (at this time, kcmd=ah(k−1), then after the limiting process, the adaptive upper limit ah is updated in such a direction that the upper limit of the adaptive allowable range increases, by the increasing unit change ΔINC in each control cycle to a limit provided by the upper-limit fourth predetermined value HH in STEP 10 - 20 . At the same time, the adaptive lower limit al is updated in such a direction that the lower limit of the adaptive allowable range decreases, by the decreasing unit change ΔDEC in each control cycle to a limit provided by the lower-limit first predetermined value STABL.

Similarly, if the demand differential air-fuel ratio usl deviates from the adaptive allowable range beyond the lower limit thereof (at this time, kcmd=al(k−1), then after the limiting process, the adaptive lower limit al is updated in such a direction that the lower limit of the adaptive allowable range increases, by the increasing unit change ΔINC in each control cycle to a limit provided by the lower-limit fifth predetermined value LL in STEP 10 - 17 . At the same time, the adaptive upper limit ah is updated in such a direction that the upper limit of the adaptive allowable range decreases, by the decreasing unit change ΔDEC in each control cycle to a limit provided by the upper-limit first predetermined value STABH.

The manner in which the adaptive allowable range thus changes with time and also the manner in which the demand differential air-fuel ratio usl changes with time are illustrated in FIG. 21 . As shown in FIG. 21 , when the demand differential air-fuel ratio usl lies in the adaptive allowable range (al<usl<ah), both the adaptive upper and lower limits ah, al are updated in such a direction that the adaptive allowable range decreases, for thereby progressively reducing the upper and lower limits of the adaptive allowable range.

If the demand differential air-fuel ratio usl deviates from the adaptive allowable range beyond either one of the adaptive upper and lower limits ah, al (usl>ah or usl<al), one of the adaptive upper and lower limits ah, al beyond which the demand differential air-fuel ratio usl deviates is updated in such a direction that the adaptive allowable range increases, and hence progressively increases. At the same time, the other of the adaptive upper and lower limits ah, al is updated in such a direction that the adaptive allowable range decreases, and hence progressively decreases.

Except when the SLD control status is in the high-level unstable state or in the low-level unstable state, if the internal combustion engine 1 is in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, then the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after FC/start/operation with lean mixture, i.e., the allowable range whose upper limit is widest, whose upper limit AHF or AHFI is set to the fixed upper-limit fifth predetermined value AFCH. Thus, the upper limit AHF or AHFI is inhibited from having a value beyond the upper-limit fifth predetermined value AFCH in the direction to reduce the allowable range in STEP 10 - 1 - 6 .

In the allowable range after FC/start/operation with lean mixture, if the internal combustion engine 1 is in the idling state, then the lower limit ALFI is set to the fixed lower-limit third predetermined value LI which is the lower limit of the allowable range for idling. If the internal combustion engine 1 is in an operating state other than the idling state, then the lower limit ALFI of the allowable range after FC/start/operation with lean mixture is set to the variable adaptive lower limit al. The lower limit ALF (=al) is updated in each control cycle exactly in the same manner as the adaptive lower limit al of the adaptive allowable range depending on how the demand differential air-fuel ratio usl deviates from the allowable range after FC/start/operation with lean mixture, except when the demand differential air-fuel ratio usl deviates from the allowable range after FC/start/operation with lean mixture beyond its upper limit in STEP 10 - 15 , STEP 10 - 17 , STEP 10 - 21 .

The manner in which the allowable range after FC, of the allowable range after FC/start/operation with lean mixture, established as the allowable range changes with time and also the manner in which the demand differential air-fuel ratio usl changes with time are shown in FIG. 22 . Likewise, the manner in which the allowable range after start changes with time and the manner in which the allowable range after operation with lean mixture changes with time are shown in FIGS. 23 and 24 , respectively.

As shown in FIG. 22 , when the supply of fuel to the internal combustion engine 1 starts being cut off, the allowable range after FC is set up as the allowable range for limiting the demand differential air-fuel ratio usl until a predetermined time elapses after the supply of fuel is cut off, and the upper limit AHF thereof is fixedly maintained at the upper-limit fifth predetermined value AFCH. Therefore, the upper limit of the allowable range (the allowable range after FC) is made wide. In the illustrated example, the internal combustion engine 1 is in an operating state other than the idling state, and the lower limit ALF of the allowable range after FC is set to the adaptive lower limit al. Basically, the lower limit ALF of the adaptive lower limit al is sequentially updated in the same manner as with the adaptive allowable range. However, if the demand differential air-fuel ratio usl exceeds the upper limit AHF (=AFCH) of the allowable range after FC when the internal combustion engine 1 is in the state immediately after the supply of fuel is cut off (including the state in which the supply of fuel is being cut off), the lower limit ALF (=al) of the adaptive lower limit al is maintained at a constant value (not updated) according to the processing of STEP 10 - 21 .

In FIG. 22 , the allowable range before the supply of fuel starts being cut off and after elapse of a predetermined time after the supply of fuel is cut off is set to the adaptive allowable range (normal operating state). In this embodiment, since the state in which the supply of fuel is being cut off is included in the state immediately after the supply of fuel is cut off, the allowable range from the start of the cutting-off of the supply of fuel is set to the allowable range after FC and its upper limit AHF is set to the upper-limit fifth predetermined value AFCH. However, while the supply of fuel is being cut off, the air-fuel ratio is not manipulated according to the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 . Therefore, the allowable range while the supply of fuel is being cut off may be any arbitrary allowable range, and may not necessarily be the allowable range after FC. Alternatively, while the supply of fuel is being cut off, the process of limiting the demand differential air-fuel ratio usl may not be carried out.

The allowable range after start is set up as follows: As shown in FIG. 23 , when the processing of the target air-fuel ratio generator 28 starts after the internal combustion engine 1 starts to operate (the processing of the target air-fuel ratio generator 28 starts basically when the activation of both the O 2 sensor 6 and the LAF sensor 5 is confirmed, see FIG. 9 ), the demand differential air-fuel ratio usl starts being calculated. The allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after start. The upper limit AHFI thereof (in the illustrated example, the operating state of the internal combustion engine 1 after it starts is the idling state) is fixedly maintained at the upper-limit fifth predetermined value AFCH (whose upper limit is wide) until a predetermined time elapses after the internal combustion engine 1 starts. Since the internal combustion engine 1 is in the idling state at this time, the lower limit ALFI of the allowable range after start is set to the lower-limit third predetermined value LI.

If the internal combustion engine 1 after it starts is kept in the idling state, then the allowable range after elapse of the predetermined time after the start of the internal combustion engine 1 is basically set to the allowable range for idling (AHFI=HI, ALFI=LI) as shown.

The allowable range after operation with lean mixture is set up as follows: As shown in FIG. 24 , when the operation mode of the internal combustion engine 1 changes from the lean operation mode to the normal operation mode, the processing operation of the target air-fuel ratio generator 28 and the manipulation of the air-fuel ratio according to the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 are resumed. Until a predetermined time elapses after the processing operation of the target air-fuel ratio generator 28 and the manipulation of the air-fuel ratio according to the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 are resumed (immediately after the end of the lean operation mode), the allowable range after operation with lean mixture is set up as the allowable range for limiting the demand differential air-fuel ratio usl, and its upper limit AHF is fixedly maintained at the upper-limit fifth predetermined value AFCH. Therefore, the upper limit of the allowable range (the allowable range after operation with lean mixture) is made wide. In the illustrated example, the operating state of the internal combustion engine 1 after the lean operation mode is a normal operating state other than the idling state, and the lower limit ALF of the allowable range after operation with lean mixture is set to the adaptive lower limit al and sequentially updated in the same manner as with the adaptive allowable range. However, if the demand differential air-fuel ratio usl exceeds the upper limit AHF (=AFCH) of the allowable range after operation with lean mixture, the lower limit ALF (=al) of the allowable range after operation with lean mixture is kept constant according to the processing in STEP 10 - 21 .

In the lean operation mode, inasmuch as the processing operation of the target air-fuel ratio generator 28 is not carried out, the demand differential air-fuel ratio usl and the upper and lower limits AHF, ALF of the allowable range are kept at values immediately prior to the start of the lean operation mode. In the present embodiment, the period in which the allowable range after operation with lean mixture is set up as the allowable range for limiting the demand differential air-fuel ratio usl extends until the predetermined time elapses after the end of the lean operation mode. However, the allowable range after operation with lean mixture may be set up until the estimated differential output VO2 bar sequentially determined by the estimator 26 after the end of the lean operation mode, for example, i.e., the estimated value of the output of the O 2 sensor 6 determined by the estimated differential output VO2 bar becomes substantially equal to the air-fuel ratio reference value FLAF/BASE.

If the SLD control status is not in the low-level unstable state or the high-level unstable state, or the internal combustion engine 1 is not in the state immediately after the supply of fuel is cut off, or immediately after it starts, or immediately after the end of the lean operation mode, or the vehicle with the internal combustion engine 1 is in the state immediately after it starts (at this time, the internal combustion engine 1 is not idling), then the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after being driven under load (whose lower limit is narrow) which has the lower limit ALF set to the lower-limit second predetermined value VSTL, and the lower limit ALF is inhibited from having a value beyond the lower-limit second predetermined value VSTL in the direction to increase the allowable range in STEP 10 - 1 - 21 .

In the allowable range after being driven under load, the upper limit AHF is set to the variable adaptive upper limit ah in STEP 10 - 1 - 19 . The upper limit AHF (=ah) is updated in each control cycle exactly in the same manner as the adaptive upper limit ah of the adaptive allowable range depending on how the demand differential air-fuel ratio usl deviates from the allowable range after being driven under load, except when the demand differential air-fuel ratio usl deviates from the allowable range after being driven under load beyond its lower limit in STEP 10 - 15 , STEP 10 - 18 , STEP 10 - 20 .

The manner in which the allowable range after being driven under load established as the allowable range changes with time and also the manner in which the demand differential air-fuel ratio usl changes with time are shown in FIG. 25 .

As shown in FIG. 25 , when the vehicle starts running from the idling state of the internal combustion engine 1 (the internal combustion engine 1 starts to drive its load), the allowable range is set to the allowable range after being driven under load and has its lower limit ALF fixedly maintained at the lower-limit second predetermined value VSTL until a predetermined time elapses after the vehicle starts running (immediately after the vehicle starts running). Therefore, the lower limit of the allowable range (the allowable range after being driven under load) is made relatively narrow. In this case, the upper limit AHF of the allowable range after being driven under load is set to the adaptive upper limit ah and sequentially updated in the same manner as with the adaptive allowable range.

While the internal combustion engine 1 is idling before the vehicle starts running, the allowable range is normally set to the allowable range for idling (AHFI=HI, ALFI=LI). In FIG. 25 , the allowable range after elapse of the predetermined time after the vehicle starts running is set to the adaptive allowable range (normal state).

The details of the limiting process carried out in STEP 10 have been described above.

Referring back to FIG. 8 , the adder 31 in the target air-fuel ratio generator 28 adds the air-fuel ratio reference value FLAF/BASE (more specifically, the air-fuel ratio reference value FLAF/BASE determined by the reference value setting unit 11 in STEP 12 , described later on, in the preceding control cycle) to the command differential air-fuel ratio kcmd determined by the limiter 30 for thereby determining a target air-fuel ratio KCMD in the present control cycle in STEP 11 .

The target air-fuel ratio KCMD thus determined is stored in a time-series fashion in a memory (not shown) in each control cycle. When the engine-side control unit 7 b uses the target air-fuel ratio KCMD determined by the target air-fuel ratio generator 28 of the exhaust-side control unit 7 a (see STEPf in FIG. 6 ), the latest target air-fuel ratio KCMD stored in the time-series fashion is selected.

Then, the target air-fuel ratio generator 28 carries out a process of setting (updating) the air-fuel ratio reference value FLAF/BASE with the reference value setting unit 11 in STEP 12 , after which the processing in the present control cycle is finished.

In the present embodiment, the air-fuel ratio reference value FLAF/BASE is defined as the sum of a predetermined fixed component flaf/base (hereinafter referred to as a “reference value fixed component flaf/base”) and a variable component flaf/adp (hereinafter referred to as a “reference value variable component flaf/adp”), i.e., FLAF/BASE=flaf/base+flaf/adp. For variably setting up the air-fuel ratio reference value FLAF/BASE, the value of the reference value variable component flaf/adp is adjusted. In the present embodiment, the reference value fixed component flaf/base is regarded as a “stoichiometric air-fuel ratio”.

The processing of STEP 12 is carried out according to a flowchart shown in FIG. 26 .

The reference value setting unit 11 determines whether the differential output VO2 of the O 2 sensor 6 is of value within a predetermined range near “0” (a range whose lower and upper limits are set respectively to fixed values ADL (<0), ADH (>0), hereinafter referred to as a “convergence determining range”) or not in STEP 12 - 1 . Stated otherwise, the reference value setting unit 11 determines whether the output VO2/OUT of the O 2 sensor 6 is converged substantially to its target value VO2/TARGET or not. In the present embodiment, the absolute values of the lower and upper limits ADL, ADH of the convergence determining range are identical to each other (|ADL|=|ADH|).

If the differential output VO2 falls in the convergence determining range (ADL<VO2<ADH) and hence the output VO2/OUT of the O 2 sensor 6 is converged substantially to the target value VO2/TARGET, then the reference value setting unit 11 compares the value of the stability determining basic parameter Pstb (see STEP 9 - 1 shown in FIG. 16 ) determined by the limiter 30 in STEP 9 (the process of determining the stability of the SLD control status) with a predetermined value δ (>0) in STEP 12 - 2 , for thereby determining the stability of the SLD control status as in STEP 9 - 4 shown in FIG. 16 .

The predetermined value δ to be compared with the stability determining basic parameter Pstb in STEP 12 - 2 is smaller than the predetermined value ε used in STEP 9 - 4 shown in FIG. 16 , thus making stricter the condition to determine that the SLD control status is stable.

If Pstb<δ, determining that the SLD control status is stable, then the reference value setting unit 11 adjust the reference value variable component flaf/adp depending on the adaptive control law input uadp (the adaptive control law component of the demand differential air-fuel ratio usl) determined by the sliding mode controller 27 in STEP 8 - 3 shown in FIG. 15 in STEP 12 - 3 through STEP 12 - 7 .

More specifically, the reference value setting unit 11 compares the value of the adaptive control law input uadp with a predetermined range near “0” (a range whose lower and upper limits are set respectively to predetermined (fixed) values NRL (<0), NRH (>0), hereinafter referred to as a “reference value adjusting dead zone”) or not in STEP 12 - 3 , STEP 12 - 5 . In the present embodiment, the absolute values of the lower and upper limits NRL, NRH of the reference value adjusting dead zone are identical to each other (|NRL|=|NRH|).

If the adaptive control law input uadp is smaller than the lower limit NRL of the reference value adjusting dead zone (uadp<NRL), then the reference value setting unit 11 subtracts a predetermined (constant) change Δflaf (>0, hereinafter referred to as a “reference value unit change Δflaf”) from a present value flaf/adp(k−1) (a value determined in the preceding control cycle) of the reference value variable component flaf/adp for thereby determining a new reference value variable component flaf/adp(k) in STEP 12 - 4 . That is, the reference value variable component flaf/adp is reduced by the reference value unit change Δflaf.

If the adaptive control law input uadp is greater than the upper limit NRH of the reference value adjusting dead zone (uadp>NRH), then the reference value setting unit 11 adds the reference value unit change Δflaf to the present value flaf/adp(k−1) of the reference value variable component flaf/adp for thereby determining a new reference value variable component flaf/adp(k) in STEP 12 - 6 . That is, the reference value variable component flaf/adp is increased by the reference value unit change Δflaf.

If the adaptive control law input uadp falls in the reference value adjusting dead zone (NRL≦uadp≦NRH), then the value of the reference value variable component flaf/adp is not changed, but kept at the present value flaf/adp(k−1) in STEP 12 - 7 .

Then, the reference value setting unit 11 adds the value of the reference value variable component flaf/adp(k) determined in either one of STEP 12 - 4 , STEP 12 - 6 , STEP 12 - 7 to the reference value fixed component flaf/base for thereby determining the air-fuel ratio reference value FLAF/BASE to be used to determine the target air-fuel ratio KCMD in STEP 11 in the next control cycle in STEP 12 - 8 . Thereafter, control returns to the main routine shown in FIG. 8 .

If the output VO2/OUT of the O 2 sensor 6 is not converged to its target value VO2/TARGET (VO2≦ADL or VO2≧ADH) in STEP 12 - 1 , or if the SLD control status is unstable (Pstb≧δ) in STEP 12 - 2 , then the reference value variable component flaf/adp is not changed, but STEP 12 - 7 is executed to hold the value of the reference value variable component flaf/adp at the present value flaf/adp(k−1). Then, STEP 12 - 8 is executed to determine the air-fuel ratio reference value FLAF/BASE, after which control returns to the main routine shown in FIG. 8 .

The reference value variable component flaf/adp which is varied depending on the adaptive control law input uadp is stored in a nonvolatile memory (not shown), e.g., an EEPROM, so that it will not be lost when the internal combustion engine 1 is stopped in operation and the exhaust-side control unit 7 a is turned off. The stored value of the reference value variable component flaf/adp will be used as an initial value of the reference value variable component flaf/adp when the internal combustion engine 1 is operated next time. The initial value of the internal combustion engine 1 when the internal combustion engine 1 is operated for the first time is “0”.

The manner in which the adaptive control law input uadp varies in STEP 12 and the manner in which the reference value variable component flaf/adp and the air-fuel ratio reference value FLAF/BASE vary depending on the adaptive control law input uadp are shown respectively in upper and lower stages in FIG. 27 .

As shown in FIG. 27 , when the adaptive control law input uadp falls in the reference value adjusting dead zone (periods T 1 , T 3 , T 5 ), the reference value variable component flaf/adp and the air-fuel ratio reference value FLAF/BASE are not varied, but kept constant. When the adaptive control law input uadp is greater than the upper limit NRH of the reference value adjusting dead zone (a period T2 in FIG. 27 ), the reference value variable component flaf/adp and the air-fuel ratio reference value FLAF/BASE are increased by the reference value unit change Δflaf in each control cycle. When the adaptive control law input uadp is smaller than the lower limit NRL of the reference value adjusting dead zone (a period T4 in FIG. 27 ), the reference value variable component flaf/adp and the air-fuel ratio reference value FLAF/BASE are reduced by the reference value unit change Δflaf in each control cycle. In this manner, the adaptive control law input uadp is finally converted to such a value that the adaptive control law input uadp falls within the reference value adjusting dead zone.

Details of the operation of the plant control system according to the present embodiment have been described above.

The operation of the plant control system is summarized as follows: Basically, the manipulated variable generator 29 (the sliding mode controller 27 , the estimator 26 , and the identifier 25 ) of the target air-fuel ratio generator 28 sequentially generates the demand differential air-fuel ratio usl as an input to be given to the object exhaust system E for converging the output signal VO2/OUT from the O 2 sensor 6 disposed downstream of the catalytic converter 3 to the target value VO2/TARGET, i.e., the difference between the air-fuel ratio of the internal combustion engine 1 and the air-fuel ratio reference value FLAF/BASE, required to converge the output signal VO2/OUT to the target value VO2/TARGET. The air-fuel ratio reference value FLAF/BASE is added to the command differential air-fuel ratio kcmd (basically kcmd=usl) produced by limiting the demand differential air-fuel ratio usl for thereby sequentially determining the target air-fuel ratio KCMD. The engine-side control unit 7 b adjusts the fuel injection quantity of the internal combustion engine 1 in order to converge the output of the LAF sensor 5 (detected air-fuel ratio) to the target air-fuel ratio KCMD for thereby feedback-controlling the air-fuel ratio of the internal combustion engine 1 at the target air-fuel ratio KCMD. In this manner, the output signal VO2/OUT from the O 2 sensor 6 is converged to the target value VO2/TARGET, making the catalytic converter 3 capable of performing an optimum exhaust gas purifying capability regardless of aging thereof.

In the present embodiment, the sliding mode controller 27 performs the adaptive sliding mode control process which is highly stable against the effect of disturbances to calculate the demand differential air-fuel ratio usl. In calculating the demand differential air-fuel ratio usl, the estimator 26 sequentially determines the estimated differential output VO2 bar which is the estimated value of the differential output VO 2 of the O 2 sensor 6 after the total dead time d which is the sum of the dead time d1 of the object exhaust system E including the catalytic converter 3 and the dead time t2 of the air-fuel ratio manipulating system (comprising the internal combustion engine 1 and the engine-side control unit 7 b ). The adaptive sliding mode control process constructed using the estimated differential output VO 2 bar, more specifically, the adaptive sliding mode control process constructed for converging the estimated differential output VO2 bar to “0”, is performed by the sliding mode controller 27 to determine the demand differential air-fuel ratio usl.

In this fashion, it is possible to generate the demand differential air-fuel ratio usl suitable for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET while compensating for the effect of the dead times d1, d2 of the object exhaust system E and the air-fuel ratio manipulating system.

The gain coefficients a1, a2, b1 as parameters of the exhaust system model (the model for expressing the behavior of the object exhaust system E) required for performing the adaptive sliding mode control process with the sliding mode controller 27 and the process of calculating the estimated differential output VO2 bar with the estimator 26 are sequentially identified in a real-time basis by the identifier 25 . The sliding mode controller 27 and estimator 26 determine the demand differential air-fuel ratio usl and the estimated differential output VO2 bar using the identified gain coefficients a1 hat, a2 hat, b1 hat which are the identified values of the gain coefficients a1, a2, b1.

Therefore, the identified gain coefficients a1 hat, a2 hat, b1 hat have their values highly reliable as depending on the behavior of the object exhaust system E, thus increasing the reliability of the demand differential air-fuel ratio usl and the estimated differential output VO2 bar.

In the engine-side control unit 7 b for manipulating the air-fuel ratio of the internal combustion engine 1 , primarily the adaptive controller 18 as a recursive-type controller capable of appropriately compensating for behavioral changes of the internal combustion engine 1 controls the air-fuel ratio of the internal combustion engine 1 at the target air-fuel ratio KCMD. Therefore, the air-fuel ratio of the internal combustion engine 1 can accurately be controlled at the target air-fuel ratio KCMD, and hence the output signal VO2/OUT of the O 2 sensor 6 can stably and accurately be converged to the target value VO2/TARGET.

In the present embodiment, the demand differential air-fuel ratio usl generated by the manipulated variable generator 29 for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET is limited to produce the command differential air-fuel ratio kcmd limited in the allowable range, by which the target air-fuel ratio KCMD is finally determined. Therefore, even when the demand differential air-fuel ratio usl suffers a large variation such as a spike, the target air-fuel ratio KCMD and hence the actual air-fuel ratio of the internal combustion engine 1 are prevented from varying excessively, so that the internal combustion engine 1 can operate stably.

In the present embodiment, the adaptive allowable range is normally set up as the allowable range for limiting the demand differential air-fuel ratio usl. The adaptive upper and lower limits ah, al of the adaptive allowable range are sequentially variably updated depending on how the demand differential air-fuel ratio usl deviates from the allowable range, specifically, depending on the magnitude of the demand differential air-fuel ratio usl with respect to the upper and lower limits (see FIG. 21 ).

Consequently, the adaptive allowable range can be made to match exactly the range of values of the demand differential air-fuel ratio usl as an input to be given to the object exhaust system E for converging the output signal VO2/OUT from the O 2 sensor 6 to the target value VO2/TARGET.

For example, in a situation in which the range of values of the demand differential air-fuel ratio usl that are generated steadily is relatively narrow, i.e., in a situation in which the demand differential air-fuel ratio usl falls steadily in the adaptive allowable range, the adaptive allowable range can be narrowed to match the narrow range. Conversely, in a situation in which the range of values of the demand differential air-fuel ratio usl that are generated steadily is relatively wide, i.e., in a situation in which the demand differential air-fuel ratio usl frequently falls out of the adaptive allowable range, the adaptive allowable range can be widened to match the wide range. In a situation in which the range of values of the demand differential air-fuel ratio usl is shifted toward the upper limit (positive) of the adaptive allowable range, i.e., in a situation in which the demand differential air-fuel ratio usl frequently deviates from the adaptive allowable range beyond its upper limit, the adaptive allowable range can be shifted toward the positive limit to match the shifted range.

Thus, the demand differential air-fuel ratio usl suitable for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET tends to fall in the allowable range (adaptive allowable range), resulting in increased opportunities to determine the target air-fuel ratio KCMD using the demand differential air-fuel ratio usl as the command differential air-fuel ratio kcmd (KCMD=usl+FLAF/BASE). Since the demand differential air-fuel ratio usl is generated according to the adaptive sliding mode control process, the control process for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET can stably be performed.

In a situation in which the demand differential air-fuel ratio usl temporarily varies greatly into a spike due to the effect of a disturbance and deviates from the allowable range (adaptive allowable range), the command differential air-fuel ratio kcmd is limited to a value in the allowable range for thereby avoiding large variations of the target air-fuel ratio KCMD to allow the internal combustion engine 1 to operate stably. In this case, because the adaptive allowable range is updated after limiting the demand differential air-fuel ratio usl, the determination of the target air-fuel ratio KCMD using the above inappropriate target air-fuel ratio KCMD directly and the resulting manipulation of the air-fuel ratio of the internal combustion engine 1 can reliably be eliminated.

For updating the adaptive allowable range, the increasing unit change ΔINC which defines a change in one control cycle of the adaptive upper limit ah or the adaptive lower limit al in a direction to increase the upper limit or the lower limit is greater than the decreasing unit change ΔDEC which defines a change in one control cycle of the adaptive upper limit ah or the adaptive lower limit al in a direction to decrease the upper limit or the lower limit. Therefore, in a situation in which the demand differential air-fuel ratio usl suitable for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET deviates from the adaptive allowable range, the adaptive allowable range can quickly be changed to an allowable range in which the demand differential air-fuel ratio usl falls.

The adaptive upper limit ah of the adaptive allowable range is limited to a range between the upper-limit first predetermined value STABH and the upper-limit fourth predetermined value HH, and the adaptive lower limit al of the adaptive allowable range is limited to a range between the lower-limit first predetermined value STABL and the lower-limit fifth predetermined value LL (see STEP 10 - 1 - 11 through STEP 10 - 1 - 18 shown in FIG. 20 ). Therefore, the target air-fuel ratio KCMD is prevented from becoming an excessively lean or rich value which is not suitable for smoothly operating the internal combustion engine 1 .

In the present embodiment, the allowable range for limiting the demand differential air-fuel ratio usl can be set to the variable adaptive allowable range, and also can be set to any of corresponding allowable ranges if the SLD control status is unstable or the internal combustion engine 1 is in a certain operating state, e.g., in a state immediately after the supply of fuel is cut off or in an idling state, as described. As a result, the following advantages are offered.

If the SLD control status is unstable, the output signal VO2/OUT of the O 2 sensor 6 tends to be unstable and suffer variations. However, according to the present embodiment, the relatively narrow allowable range for unstable low level or the allowable range for unstable high level is established as the allowable range for limiting the demand differential air-fuel ratio usl. Therefore, variations of the command differential air-fuel ratio kcmd and hence the target air-fuel ratio KCMD are suppressed, resulting in stability of the output signal VO2/OUT of the O 2 sensor 6 .

The allowable range established if the SLD control status is unstable is narrower when it is unstable at high level than when it is unstable at low level (see FIG. 18 ). That is, the more unstable the SLD control status, the narrower the allowable range. If the SLD control status is unstable at high level, therefore, variations of the command differential air-fuel ratio kcmd and hence the target air-fuel ratio KCMD can be suppressed as positively as possible, so that the stability of the output signal VO2/OUT of the O 2 sensor 6 can reliably be achieved. If the SLD control status is unstable at low level, while the output signal VO2/OUT of the O 2 sensor 6 is rendered stable, the ability of the output signal VO2/OUT to the target value VO2/TARGET is achieved to some extent.

In this embodiment, for the predetermined time (: TMSTB) after the air-fuel ratio of the internal combustion engine 1 starts being manipulated according to the target air-fuel ratio KCMD generated by the target air-fuel ratio generator 28 , i.e., after the control process of converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET starts, the stability of the SLD control status is not determined, i.e., the SLD control status is regarded as being stable (see STEPd- 11 shown in FIG. 7 and STEP 9 - 2 , STEP 9 - 3 shown in FIG. 16 ), so that the allowable range for unstable low level and the allowable range for unstable high level are prevented from being set up as the allowable range for limiting the demand differential air-fuel ratio usl. Specifically, immediately after the control process of converging the output signal VO2/OUT of the O 2 sensor 6 starts, the output signal VO2/OUT is not converged to the target value VO2/TARGET. In order to accelerate the convergence of the output signal VO2/OUT to the target value VO2/TARGET, limiting the command differential air-fuel ratio kcmd to the upper or lower limit of the allowable range is avoided as much as possible, and the frequency at which the demand differential air-fuel ratio usl becomes the command differential air-fuel ratio kcmd is increased. Thus, the output signal VO2/OUT of the O 2 sensor 6 can quickly be brought closely to the target value VO2/TARGET.

Immediately after the supply of fuel to the internal combustion engine 1 is cut off, as shown in FIG. 22 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after FC with its upper limit being wide. Specifically, since a large amount of oxygen is stored in the catalytic converter 3 while the supply of fuel to the internal combustion engine 1 is being cut off, the output signal VO2/OUT of the O 2 sensor 6 is spaced away from the target value VO2/TARGET toward a leaner value of the air-fuel ratio (the output signal VO2/OUT of the O 2 sensor 6 becomes smaller, see FIG. 2 ). Therefore, immediately after the supply of fuel to the internal combustion engine 1 is cut off, the demand differential air-fuel ratio usl for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET becomes large in a direction to make the air-fuel ratio richer (in the present embodiment, toward the upper limit of the allowable range in a positive direction of the demand differential air-fuel ratio usl), as shown in FIG. 22 . Therefore, the upper limit of the allowable range after FC is set to the upper-limit fifth predetermined value AFCH, making the allowable range wider at the upper limit. It is thus possible for the command differential air-fuel ratio kcmd immediately after the supply of fuel to the internal combustion engine 1 is cut off to have a large value at the upper limit of the allowable range according to the demand differential air-fuel ratio usl, allowing the output signal VO2/OUT of the O 2 sensor 6 to converge quickly to the target value VO2/TARGET.

In the state immediately after the start of the internal combustion engine 1 , as shown in FIG. 23 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after start, with its upper limit made wide in the same manner as immediately after the supply of fuel is cut off. Specifically, since oxygen delivered to the catalytic converter 3 upon cranking of the internal combustion engine 1 to start same is stored in the catalytic converter 3 , immediately after the internal combustion engine 1 starts, the output signal VO2/OUT of the O 2 sensor 6 tends to be spaced away from the target value VO2/TARGET toward a leaner value of the air-fuel ratio. In the present embodiment, therefore, the upper limit (corresponding to the leaner value of the air-fuel ratio) of the allowable range after start is made wide in the same manner as immediately after the supply of fuel is cut off. Therefore, the output signal VO2/OUT of the O 2 sensor 6 can be converged quickly to the target value VO2/TARGET.

In the state immediately after the end of the lean operation mode of the internal combustion engine 1 (immediately after the lean operation mode changes to the normal operation mode), as shown in FIG. 24 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after operation with lean mixture, with its upper limit made wide in the same manner as immediately after the supply of fuel is cut off. Specifically, since the air-fuel ratio of the internal combustion engine 1 is manipulated toward a leaner value during the lean operation mode of the internal combustion engine 1 , the output signal VO2/OUT of the O 2 sensor 6 is spaced away from the target value VO2/TARGET toward a leaner value of the air-fuel ratio. In the present embodiment, therefore, the upper limit of the allowable range after operation with lean mixture is made wide in the same manner as immediately after the supply of fuel is cut off. Therefore, it is possible for the command differential air-fuel ratio kcmd immediately after the end of the lean operation mode to have a large value at the upper limit of the allowable range according to the demand differential air-fuel ratio usl. Consequently, the output signal VO2/OUT of the O 2 sensor 6 can be converged quickly to the target value VO2/TARGET, and the catalytic converter 3 can quickly achieve an appropriate exhaust gas purifying capability.

At this time, inasmuch as the command differential air-fuel ratio kcmd has a large value at the upper limit of the allowable range (corresponding to a richer value of the air-fuel ratio), the actual air-fuel ratio of the internal combustion engine 1 quickly becomes a rich air-fuel ratio after the end of the lean operation mode. Thus, NOx absorbed by the NOx absorber (not shown) contained in the catalytic converters 3 , 4 during the lean operation mode of the internal combustion engine 1 can quickly be reduced. In a next cycle of the lean operation mode, therefore, NOx in the exhaust gases can sufficiently be absorbed by the NOx absorber contained in the catalytic converters 3 , 4 for exhaust gas purification.

In the state immediately after the start of the vehicle, i.e., in the state immediately after the internal combustion engine 1 starts driving its load, as shown in FIG. 25 , the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range after being driven under load, with its upper limit made narrow. Specifically, immediately after the internal combustion engine 1 starts driving its load, the air-fuel ratio of the internal combustion engine 1 tends to change to a leaner value. If such a situation occurs when the air-fuel ratio is changed in a direction toward a leaner value by the demand differential air-fuel ratio usl (in the present embodiment, in a negative direction of the demand differential air-fuel ratio usl toward the lower limit of the allowable range), then the output signal VO2/OUT of the O 2 sensor 6 is likely to change excessively to a leaner air-fuel ratio with respect to the target value VO2/TARGET. According to the present embodiment, the lower limit of the allowable range after being driven under load is set to the lower-limit second predetermined value VSTL, thus making the allowable range narrower at the lower limit. The magnitude (absolute value) of the value of the command differential air-fuel ratio kcmd which can be taken at the lower limit of the allowable range immediately after the start of the vehicle, i.e., immediately after the internal combustion engine 1 starts driving its load, is limited to a relatively small value, thus making the output signal VO2/OUT of the O 2 sensor 6 stable.

In the present embodiment, the demand differential air-fuel ratio usl is defined as a value produced by subtracting the air-fuel ratio reference value FLAF/BASE from the air-fuel ratio to be given to the object exhaust system E for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET. Therefore, the upper limit (positive) of the allowable range corresponds to a richer air-fuel ratio and the lower limit (negative) of the allowable range corresponds to a leaner air-fuel ratio. However, the demand differential air-fuel ratio usl may be defined as a value produced by subtracting the air-fuel ratio to be given to the object exhaust system E for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, from the air-fuel ratio reference value FLAF/BASE. In this case, since the sign of the demand differential air-fuel ratio usl is opposite to the sign in the embodiment, the upper limit (positive) of the allowable range corresponds to a leaner air-fuel ratio and the lower limit (negative) of the allowable range corresponds to a richer air-fuel ratio.

While the internal combustion engine 1 is idling, the allowable range for limiting the demand differential air-fuel ratio usl is set to the allowable range for idling with both upper and lower limits thereof being made relatively narrow. Specifically, in the idling state of the internal combustion engine 1 , if the air-fuel ratio of the internal combustion engine 1 is changed largely, then the stability of the idling state tends to be impaired. According to the present embodiment, therefore, the narrow allowable range for idling is established as the allowable range for limiting the demand differential air-fuel ratio usl. In this manner, the output signal VO2/OUT of the O 2 sensor 6 can be converged to the target value VO2/TARGET while keeping the idling state of the internal combustion engine 1 stable.

According to the present embodiment, furthermore, the air-fuel ratio reference value FLAF/BASE is variably established, as described above, depending on the adaptive control law input uadp which is a component based on the adaptive control law (adaptive algorithm) of the demand differential air-fuel ratio usl generated by the sliding mode controller 27 . In this manner, the air-fuel ratio reference value FLAF/BASE can be of a central value in a range of values of the air-fuel ratio (=usl+FLAF/BASE) that is the sum of the air-fuel ratio reference value FLAF/BASE and the demand differential air-fuel ratio usl, i.e., the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, which air-fuel ratio is basically equal to the target air-fuel ratio KCMD. As a result, values of the demand differential air-fuel ratio usl which are sequentially generated by the sliding mode controller 27 can be balanced between positive and negative values, and the adaptive allowable range that is made variable depending on the demand differential air-fuel ratio usl can be balanced between its upper limit (positive) and its lower limit (negative).

Specifically, with the output signal VO2/OUT of the O 2 sensor 6 steadily converged to the target value VO2/TARGET, as is apparent from the equations (24)-(26), the equivalent control law input ueq and the reaching control law input urch, which are components other than the adaptive control law input uadp of the demand differential air-fuel ratio usl, become “0”, and usl=uadp. The adaptive control law input uadp is thus significant as a central value in the range of values of the demand differential air-fuel ratio usl with the output signal VO2/OUT of the O 2 sensor 6 steadily converged to the target value VO2/TARGET. The sum of the adaptive control law input uadp and the air-fuel ratio reference value FLAF/BASE is significant as a central value of the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, i.e., the target air-fuel ratio KCMD.

By adjusting the target value VO2/TARGET such that the adaptive control law input uadp will be of a value close to “0”, the air-fuel ratio reference value FLAF/BASE can be of a central value of the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, i.e., the target air-fuel ratio KCMD. In this embodiment, therefore, by appropriately changing the air-fuel ratio reference value FLAF/BASE depending on the value of the adaptive control law input uadp, the air-fuel ratio reference value FLAF/BASE is adjusted to keep the value of the adaptive control law input uadp in the reference value adjusting dead zone. In this fashion, the air-fuel ratio reference value FLAF/BASE can be adjusted to a central value of the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, i.e., the target air-fuel ratio KCMD. As a result, values of the demand differential air-fuel ratio usl can be balanced between positive and negative values, and the adaptive allowable range that is made variable depending on the demand differential air-fuel ratio usl can be balanced between its upper limit (positive) and its lower limit (negative). With the adaptive allowable range balanced between its upper limit (positive) and its lower limit (negative), the demand differential air-fuel ratio usl can be limited appropriately in a balance at both the upper and lower limits.

In this embodiment, the air-fuel ratio reference value FLAF/BASE is adjusted (updated) only when the output signal VO2/OUT of the O 2 sensor 6 is substantially converged to the target value VO2/TARGET and it is determined that the SLD control status is stable from the stability determining basic parameter Pstb. Thus, the air-fuel ratio reference value FLAF/BASE is adjusted when the value of the adaptive control law input uadp is stabilized, and the reliability of the air-fuel ratio reference value FLAF/BASE as a central value of the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, i.e., the target air-fuel ratio KCMD, is increased.

In this embodiment, in adjusting the air-fuel ratio reference value FLAF/BASE, the value of the air-fuel ratio reference value FLAF/BASE is not changed while the adaptive control law input uadp is present in the reference value adjusting dead zone. Thus, it is possible to avoid frequent variations of the air-fuel ratio reference value FLAF/BASE for thereby avoiding situations in which the SLD control status is unstable.

The above adjustment of the air-fuel ratio reference value FLAF/BASE offers the following advantages:

By changing the air-fuel ratio reference value FLAF/BASE depending on the adaptive control law input uadp, the quick response of the control process for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET can be increased. Specifically, in the case where the air-fuel ratio reference value FLAF/BASE is constant, e.g., FLAF/BASE=flaf/base, then if there is a steady error between the actual air-fuel ratio of the internal combustion engine 1 and the target air-fuel ratio KCMD, then the adaptive control law input uadp determined by the sliding mode controller 27 finally corresponds to a learned value of the error. If the error is relatively large, then it takes a period of time for the adaptive control law input uadp to finally correspond to a learned value of the error. According to the present invention, since the air-fuel ratio reference value FLAF/BASE is changed depending on the adaptive control law input uadp, the adaptive control law input uadp can be of a sufficiently small value close to “0”. Stated otherwise, the above error can be absorbed by the air-fuel ratio reference value FLAF/BASE, with the result that the quick response of the control process for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET can be increased.

Furthermore, by changing the air-fuel ratio reference value FLAF/BASE depending on the adaptive control law input uadp, the accuracy of the estimated differential output VO2 bar determined by the estimator 26 and the identified gain coefficients a1 hat, a2 hat, b1 hat determined by the identifier 25 can be increased. The reasons for this are as follows: The exhaust system model expressed by the equation (1) with the air-fuel ratio reference value FLAF/BASE used as a reference for the input to the object exhaust system E is a model in which the output KACT (detected air-fuel ratio) of the LAF sensor 5 becomes the air-fuel ratio reference value FLAF/BASE with the output signal VO2/OUT of the O 2 sensor 6 steadily converged to the target value VO2/TARGET. Therefore, the air-fuel ratio reference value FLAF/BASE should be a central value of the air-fuel ratio of the internal combustion engine 1 with the output signal VO2/OUT of the O 2 sensor 6 steadily converged to the target value VO2/TARGET. In this embodiment, by changing the air-fuel ratio reference value FLAF/BASE depending on the adaptive control law input uadp, the air-fuel ratio reference value FLAF/BASE can be adjusted to a central value of the air-fuel ratio required for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET. As a consequence, the behavior of the exhaust system model can be made to better match the actual behavior of the object exhaust system. Therefore, the accuracy of the estimated differential output VO2 bar determined by the estimator 26 based on the exhaust system model can be increased, and the accuracy of the identified gain coefficients a1 hat, a2 hat, b1 hat determined by the identifier 25 as identified values of parameters of the exhaust system model can be increased. With the accuracy of the estimated differential output VO2 bar and the identified gain coefficients a1 hat, a2 hat, b1 being increased, the demand differential air-fuel ratio usl determined by the sliding mode controller 27 using these data can be made optimum for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET. As a result, the accuracy of the control process for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET can be increased.

A plant control system according to another embodiment of the present invention will be described below. The embodiment basically differs from the above embodiment only as to the processing carried out by the estimator 26 . Therefore, the same reference characters as those of the above embodiment will be used in the description of the other embodiment.

In the above embodiment, in order to compensate for the effect of the total dead time d which is the sum of the dead time d1 of the object exhaust system E and the dead time d2 of the air-fuel ratio manipulating system (comprising the internal combustion engine 1 and the engine-side control unit 7 b ), an estimated value of the differential output VO2 of the O 2 sensor 6 (estimated differential output VO2 bar) after the total dead time d is determined. If the dead time d2 of the air-fuel ratio manipulating system is sufficiently small as compared with the dead time d1 of the object exhaust system E, then an estimated value VO2 (k+d1) of the differential output VO2 of the O 2 sensor 6 after the dead time d1 of the object exhaust system E (hereinafter referred to as a “second estimated differential output VO2 bar”) may be determined, and the demand differential air-fuel ratio usl may be determined using the second estimated differential output VO2 bar. In this embodiment, the second estimated differential output VO2 bar is determined to converge the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET.

In this case, the estimator 26 determines the second estimated differential output VO2 bar as follows: Using the equation (1) representing the exhaust system model of the object exhaust system E, the second estimated differential output VO2 bar which represents an estimated value of the differential output VO2 of the O 2 sensor 6 after the dead time d 1 of the object exhaust system E in each control cycle is expressed, using time-series data VO2 (k), VO2 (k−1) of the differential output VO2 of the O 2 sensor 6 and time-series data kact(k−j) (j=1, 2, . . . , d1) of past values of the differential output kact (kact=KACT−FLAF/BASE) of the LAF sensor 5 , according to the following equation (42): $$ VO2 _ ⁡ ( k + d1 ) =   ⁢ α ⁢ 3 · VO2 ( k ) + α ⁢ 4 · VO2 ( k - 1 ) +   ⁢ ∑ j = 1 d1 ⁢   ⁢ γ j · kact ⁡ ( k - j ) ( 42 )

where

α3=the first-row, first-column element of A d1 ,

α4=the first-row, second-column element of A d1 ,

γj=the first-row elements of A j−1 ·B $$ A = [ a1 a2 1 0 ] B = [ b1 0 ]

In the equation (42), α3, α4 represent the first-row, first-column element and the first-row, second-column element of the d 1 th power A d1 (d1: total dead time of the object exhaust system E) of the matrix A defined in the equation (12), and γj (j=1, 2, . . . , d1) represents the first-row elements of the product A j−1 ·B of the (j−1)th power A j−1 (j=1, 2, . . . , d1) of the matrix A and the vector B defined in the equation (12).

The above equation (42) is a basic equation for the estimator 26 to calculate the second estimated differential output VO2 (k+d1) bar in this embodiment. The equation (42) is similar to the equation (12) except that “kcmd” in the equation (12) is replaced with “kact” and “d” in the equation (12) is replaced with “d1”. In this embodiment, the equation (42) is calculated using time-series data VO2(k), VO2(k−1) of the differential output VO2 of the O 2 sensor 6 and time-series data kact(k−j) (j=1, 2, . . . , d1) of past values of the differential output kact (kact=KACT−FLAF/BASE) of the LAF sensor 5 in each control cycle, to determine the second estimated differential output VO2(k+d1) bar of the O 2 sensor 6 .

The values of the coefficients α3, α4 and γj (j=1, 2, . . . , d1) which are necessary to calculate the second estimated differential output VO2(k+d1) bar according to the equation (42) are calculated using the identified gain coefficients a1 hat, a2 hat, b1 hat which are identified values of the gain coefficients a1, a2, b1. The value of the dead time d1 required in the calculation of the equation (42) is the same as in the previous embodiment.

The other processing than the processing described above is basically the same as in the previous embodiment. However, the sliding mode controller 27 determines the equivalent control input ueq, the reaching control law input urch, and the adaptive control law uadp, which are components of the demand differential air-fuel ratio usl according to the equations (24), (26), (27), respectively, where “d” is replaced with “d1”. In the process of limiting combinations of the identified gain coefficients a1 hat, a2 hat, carried out by the identifier 25 , depending on the value of the dead time d1 of the object exhaust system E, the range for limiting the combinations (which corresponds to the identifying coefficient stable range or the identifying coefficient limiting range shown in FIG. 12 ) may differ from the range in the previous embodiment, but may be established in the same manner as in the previous embodiment.

The plant control system according to this embodiment operates on the same manner and offers the same advantages as the with previous embodiment with respect to the process of limiting the demand differential air-fuel ratio usl and the variable setting of the air-fuel ratio reference value FLAF/BASE.

The air-fuel ratio control system for the internal combustion engine is not limited to the above embodiments, but may be modified as follows:

In the above embodiments, the LAF sensor (wide-range air-fuel ratio sensor) 5 is employed for detecting an air-fuel ratio of the internal combustion engine 1 . However, an ordinary O 2 sensor or any of various other types of sensors may be employed insofar as it can detect the air-fuel ratio of the internal combustion engine 1 .

In the above embodiments, the O 2 sensor 6 is employed as the exhaust gas sensor disposed downstream of the catalytic converter. However, the exhaust gas sensor may comprise any of various other types of sensors insofar as it can detect the concentration of a certain component of an exhaust gas downstream of the catalytic converter to be controlled. For example, if carbon monoxide in an exhaust gas downstream of the catalytic converter is to be controlled, the exhaust gas sensor may comprise a CO sensor. If nitrogen oxide (NOx) in an exhaust gas downstream of the catalytic converter is to be controlled, the exhaust gas sensor may comprise an NOx sensor. If hydrocarbon (HC) in an exhaust gas downstream of the catalytic converter is to be controlled, the exhaust gas sensor may comprise an HC sensor. When a three-way catalytic converter is employed, then even if the concentration of any of the above gas components is detected, it may be controlled to maximize the purifying performance of the three-way catalytic converter. If a catalytic converter for oxidation or reduction is employed, then purifying performance of the catalytic converter can be increased by directly detecting a gas component to be purified.

In the above embodiments, after the demand differential air-fuel ratio usl is limited, the air-fuel ratio reference value FLAF/BASE is added to the command differential air-fuel ratio kcmd produced by the limiting process for thereby determining the target air-fuel ratio KCMD. However, the target air-fuel ratio KCMD may be determined by limiting the sum (demand air-fuel ratio) of the demand differential air-fuel ratio usl and the air-fuel ratio reference value FLAF/BASE. In this case, the air-fuel ratio reference value FLAF/BASE is added to the upper and lower limits of the allowable range established as described above with respect to the demand differential air-fuel ratio usl for thereby determining the allowable range with respect to the demand air-fuel ratio. The demand air-fuel ratio is then limited by the allowable range to determine the target air-fuel ratio KCMD. The demand air-fuel ratio may be regarded as the manipulated variable for manipulating the air-fuel ratio of the internal combustion engine 1 in order to converge the output VO2/OUT of the O 2 sensor to the target value VO2/TARGET.

In the above embodiments, the upper limit of the allowable range immediately after the supply of fuel to the internal combustion engine 1 is cut off or not, immediately after the start of the internal combustion engine 1 , or immediately after the end of the lean operation mode, i.e., the allowable range after FC/operation with lean mixture, is fixed to the upper-limit fifth predetermined value AFCH. However, the upper limit of the allowable range may be variably set as with the adaptive upper limit ah of the adaptive allowable range so as not to be smaller than the upper-limit fifth predetermined value AFCH. Furthermore, the upper limits of the allowable range after FC, the allowable range after start, and the allowable range after operation with lean mixture may be established separately from each other.

In the above embodiments, the lower limit of the allowable range for the limiting process immediately after the start of the vehicle, i.e., the lower limit of the allowable range after being driven under load, is fixed to the lower-limit second predetermined value VSTL. However, the lower limit (<0) may be variably set as with the adaptive lower limit al of the adaptive allowable range so as not to be greater than the lower-limit second predetermined value VSTL (<0).

In the above embodiments, the gain coefficients a1, a2, b1 as parameters of the exhaust system model are identified by the identifier 25 . However, the gain coefficients a1, a2, b1 may be fixed to predetermined values, or may be established from a map or the like depending on operating conditions of the internal combustion engine 1 and deteriorated states of the catalytic converter 3 .

In the above embodiments, the estimator 26 and the sliding mode controller 27 use a common exhaust system model of the object exhaust system E. However, the estimator 26 and the sliding mode controller 27 may use respective models. In such a case, an input to the exhaust system model for use in the processing of the sliding mode controller 27 may not necessarily be expressed using the air-fuel ratio reference value FLAF/BASE.

In the above embodiments, the exhaust system model is expressed by a discrete system (discrete-time system). However, it may be expressed a continuous system (continuous-time system), and an algorithm for the processing of the estimator 26 and the sliding mode controller 27 may be constructed based on the model of such a continuous system (continuous-time system).

In the above embodiments, the adaptive sliding mode control process is used as a feedback control process for generating the manipulated variable (the demand differential air-fuel ratio usl) for manipulating the air-fuel ratio of the internal combustion engine 1 for converging the output signal VO2/OUT of the O 2 sensor 6 to the target value VO2/TARGET, using the estimated differential output VO2 bar. However, any of various other feedback control processes (preferably for generating a manipulated variable component corresponding to the adaptive control law input uadp) may be used.

In the above embodiments, the air-fuel ratio of the internal combustion engine 1 is feedback-controlled at the target air-fuel ratio KCMD using the output signal from the LAF sensor 5 . However, it is possible to manipulate the air-fuel ratio of the internal combustion engine 1 into the target air-fuel ratio KCMD by adjusting the amount of fuel supplied to the internal combustion engine 1 under feed-forward control using a map or the like based on the target air-fuel ratio KCMD.

In the above embodiments, in adjusting the air-fuel ratio reference value FLAF/BASE, the value of the air-fuel ratio reference value FLAF/BASE is not changed while the adaptive control law input uadp is present in the reference value adjusting dead zone. However, the air-fuel ratio reference value FLAF/BASE may be changed such that if the adaptive control law input uadp is greater than “0”, then the air-fuel ratio reference value FLAF/BASE is incremented by the reference value unit change Δflaf in each control cycle, and if the adaptive control law input uadp is smaller than “0”, then the air-fuel ratio reference value FLAF/BASE is decremented by the reference value unit change Δflaf in each control cycle. The air-fuel ratio reference value FLAF/BASE can then be adjusted such that the adaptive control law input uadp will be substantially “0”.

In the above embodiments, the air-fuel ratio reference value FLAF/BASE is changed depending on the adaptive control law input uadp according to the adaptive sliding mode control process. However, even if the demand differential air-fuel ratio usl is generated according the normal sliding mode control process which does not include the adaptive control law input uadp, the air-fuel ratio reference value FLAF/BASE can variably be established. Specifically, if the demand differential air-fuel ratio usl is determined as the sum of the equivalent control input ueq and the reaching control law input urch (usl=ueq+urch) according to the normal sliding mode control process which does not employ the adaptive control law, then the reaching control law input urch in a state wherein the value of the stability determining basic parameter Pstb (=σ bar·Δσ bar) or the value of the rate of change of the switching function σ bar is steadily substantially “0” corresponds to the adaptive control law input uadp in the above embodiments. Therefore, the same advantages as in the above embodiments can be achieved by changing the air-fuel ratio reference value FLAF/BASE depending on the reaching control law input urch in the above state in the same manner as in the above embodiments.

In the above embodiments, the air-fuel ratio control system for the internal combustion engine 1 mounted on the vehicle has been described by way of example. However, the present invention is also applicable to an engine for actuating an object other than the vehicle, e.g., an electric generator or the like.

Although certain preferred embodiments of the present invention have been shown and described in detail, it should be understood that various changes and modifications may be made therein without departing from the scope of the appended claims.

Claims

1. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means; wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit beyond which the manipulated variable deviates, in a direction to increase said allowable range for thereby updating the allowable range.

2. An air-fuel ratio control system according to claim 1, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit opposite to the limit beyond which the manipulated variable deviates, in a direction to decrease said allowable range for thereby updating the allowable range.

3. An air-fuel ratio control system according to claim 1, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

4. An air-fuel ratio control system according to claim 2, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

5. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means; wherein said limiting means comprises means for changing an upper limit of said allowable range in a direction to increase or reduce said allowable range depending on whether the value of the manipulated variable generated by said manipulated variable generating means is larger than said upper limit or not, and changing a lower limit of said allowable range in a direction to increase or decrease said allowable range depending on whether the value of the manipulated variable generated by said manipulated variable generating means is smaller than said lower limit or not, for thereby updating the allowable range.

6. An air-fuel ratio control system according to any one of claims 2 through 5, wherein said upper or lower limit is changed by said limiting means in the direction to increase the allowable range, by an amount greater than an amount by which said upper or lower limit is changed in the direction to decrease the allowable range.

7. An air-fuel ratio control system according to any one of claims 1 through 5, wherein said limiting means comprises means for limiting the manipulated variable to the value in the allowable range before the allowable range is updated.

8. An air-fuel ratio control system according to any one of claims 1 through 5, wherein said limiting means comprises means for limiting the upper or lower limit of the allowable range which is being changed depending on the value of the manipulated variable generated by said manipulated variable generating means, to a value in a predetermined range corresponding to the upper or lower limit of the allowable range.

9. An air-fuel ratio control system according to any one of claims 1 through 5, wherein said manipulated variable generating means comprises means for generating said manipulated variable according to a sliding mode control process.

10. An air-fuel ratio control system according to claim 9, wherein said limiting means comprises means for sequentially determining whether the output from said exhaust gas sensor is stable or not, and, if the output from said exhaust gas sensor is determined to be unstable, forcibly setting said allowable range to a predetermined range at least while the output from said exhaust gas sensor is continuously being determined to be unstable.

11. An air-fuel ratio control system according to claim 10, wherein said limiting means comprises means for, if the output from said exhaust gas sensor is determined to be unstable, determining the level of instability in a plurality of stages, and thereby making said predetermined range to which said allowable range is set by said limiting means to be determined smaller as the level of instability determined by said limiting means is higher.

12. An air-fuel ratio control system according to claim 11, wherein said limiting means comprises means for sequentially temporarily determining whether the output from said exhaust gas sensor is stable or not based on the value of a switching function used in said sliding mode control process, and determining whether the output from said exhaust gas sensor is stable or not and the level of instability based on the frequency of the temporarily determined instability within a predetermined period of time.

13. An air-fuel ratio control system according to claim 10, wherein said limiting means comprises means for determining whether the output from said exhaust gas sensor is stable or not based on the value of a switching function used in said sliding mode control process.

14. An air-fuel ratio control system according to claim 9, wherein the manipulated variable generated by said manipulated variable generating means and the upper and lower limits of said allowable range represent differences with a predetermined reference value for the air-fuel ratio of the air-fuel mixture, further comprising reference value variable setting means for variably setting said reference value depending on the manipulated variable generated by said manipulated variable generating means.

15. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means; wherein the manipulated variable generated by said manipulated variable generating means and the upper and lower limits of said allowable range represent differences with a predetermined reference value for the air-fuel ratio of the air-fuel mixture, further comprising reference value variable setting means for variably setting said reference value depending on the manipulated variable generated by said manipulated variable generating means.

16. An air-fuel ratio control system according to claim 15, wherein said reference value variable setting means comprises means for determining whether the output from said exhaust gas sensor is substantially converged to said target value or not, and holding said reference value to a predetermined value irrespective of said manipulated variable if the output from said exhaust gas sensor is not converged to said target value.

17. An air-fuel ratio control system according to claim 15, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit beyond which the manipulated variable deviates, in a direction to increase said allowable range for thereby updating the allowable range.

18. An air-fuel ratio control system according to claim 15, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit opposite to the limit beyond which the manipulated variable deviates, in a direction to decrease said allowable range for thereby updating the allowable range.

19. An air-fuel ratio control system according to claim 15, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

20. An air-fuel ratio control system according to claim 15, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

21. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means, and means for recognizing an operating state of said internal combustion engine, and setting said allowable range depending on the recognized operating state; wherein said exhaust gas sensor comprises an oxygen concentration sensor, and said operating state of the internal combustion engine recognized by said limiting means includes an elapsed time after the internal combustion engine starts to operate, and wherein said limiting means comprises means for setting said allowable range to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease said allowable range beyond a predetermined value until said elapsed time reaches a predetermined time.

22. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means, and means for recognizing an operating state of said internal combustion engine, and setting said allowable range depending on the recognized operating state; wherein said exhaust gas sensor comprises an oxygen concentration sensor, and said operating state of the internal combustion engine recognized by said limiting means includes an elapsed time after the supply of fuel to the internal combustion engine is cut off, and wherein said limiting means comprises means for setting said allowable range to inhibit at least one, which corresponds to a richer side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to decrease said allowable range beyond a predetermined value until said elapsed time reaches a predetermined time.

23. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means, and means for recognizing an operating state of said internal combustion engine, and setting said allowable range depending on the recognized operating state; wherein said exhaust gas sensor comprises an oxygen concentration sensor, and said operating state of the internal combustion engine recognized by said limiting means includes an elapsed time after the internal combustion engine starts driving a load, and wherein said limiting means comprises means for setting said allowable range to inhibit at least one, which corresponds to a leaner side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to increase said allowable range beyond a predetermined value until said elapsed time reaches a predetermined time.

24. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means, and means for recognizing an operating state of said internal combustion engine, and setting said allowable range depending on the recognized operating state; wherein said operating state of the internal combustion engine recognized by said limiting means includes an idling state or an operating state other than the idling state, and wherein said limiting means comprises means for setting a predetermined range as said allowable range irrespective of the manipulated variable generated by said manipulated variable generating means if said operating state is the idling state, and variably setting said allowable range depending on the manipulated variable generated by said manipulated variable generating means if said operating state is the operating state other than the idling state.

25. An air-fuel ratio control system for an internal combustion engine, comprising:an exhaust gas sensor disposed downstream of a catalytic converter in an exhaust passage of the internal combustion engine, for detecting the concentration of a particular component of an exhaust gas having passed through the catalytic converter; manipulated variable generating means for sequentially generating a manipulated variable which manipulates the air-fuel ratio of an air-fuel mixture combusted by the internal combustion engine in order to converge an output from said exhaust gas sensor to a predetermined target value; and limiting means for limiting the manipulated variable generated by said manipulated variable generating means to a value in a predetermined allowable range, for thereby making the air-fuel ratio of the air-fuel mixture to be determined based on the manipulated variable; said limiting means comprising means for variably setting the allowable range depending on the manipulated variable generated by said manipulated variable generating means, and means for recognizing an operating state of said internal combustion engine, and setting said allowable range depending on the recognized operating state; wherein said exhaust gas sensor comprises an oxygen concentration sensor, said internal combustion engine has operation modes including a normal operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated based on the manipulated variable limited by said limiting means, and a lean operation mode in which the air-fuel ratio of the air-fuel mixture is manipulated into a lean air-fuel ratio irrespective of the manipulated variable, and said operating state of the internal combustion engine recognized by said limiting means includes an elapsed time from a transition from said lean operation mode to said normal operation mode, and wherein said limiting means comprises means for setting said allowable range to inhibit at least one, which corresponds to a leaner side of the air-fuel ratio of the air-fuel mixture, of upper and lower limits of the allowable range from becoming a value in a direction to increase said allowable range beyond a predetermined value until said elapsed time reaches a predetermined time.

26. An air-fuel ratio control system according to claim 25, wherein said catalytic converter comprises a catalytic converter for absorbing a nitrogen oxide in the exhaust gas when the concentration of oxygen in the exhaust gas passing through said catalytic converter is a concentration corresponding to a lean air-fuel ratio of the air-fuel mixture, and reducing the nitrogen oxide in the exhaust gas when the concentration of oxygen in the exhaust gas passing through said catalytic converter is a concentration corresponding to a rich air-fuel ratio of the air-fuel mixture.

27. An air-fuel ratio control system according to claims 21-26, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit beyond which the manipulated variable deviates, in a direction to increase said allowable range for thereby updating the allowable range.

28. An air-fuel ratio control system according to claim 27, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

29. An air-fuel ratio control system according to claims 21-26, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means deviates from said allowable range beyond an upper or lower limit thereof, changing the limit opposite to the limit beyond which the manipulated variable deviates, in a direction to decrease said allowable range for thereby updating the allowable range.

30. An air-fuel ratio control system according to claims 21-26, wherein said limiting means comprises means for, if the value of the manipulated variable generated by said manipulated variable generating means is present in said allowable range, changing at least one of upper and lower limits of said allowable range in a direction to decrease said allowable range for thereby updating the allowable range.

31. An air-fuel ratio control system according to claims 17-20, wherein said sliding mode control process comprises an adaptive sliding mode control process, and said manipulated variable generated by said manipulated variable generating means according to said adaptive sliding mode control process includes an adaptive control law component based on an adaptive control law of said adaptive sliding mode control process, said reference value variable setting means comprising means for variably setting said referenced value based on the value of the adaptive control law component of said manipulated variable.

32. An air-fuel ratio control system according to claim 31, wherein said reference value variable setting means comprises means for variably setting said reference value by increasing or decreasing said reference value depending on the magnitude of the value of the adaptive control law component of said manipulated variable with respect to a predetermined value or a range close to and containing said predetermined value.

33. An air-fuel ratio control system according to claims 17-20, wherein said reference value variable setting means comprises means for sequentially determining whether the output from said exhaust gas sensor is stable or not, and holding said reference value to a predetermined value irrespective of said manipulated variable if the output from said exhaust gas sensor is unstable.

34. An air-fuel ratio control system according to claim 33, wherein said reference value variable setting means comprises means for determining whether the output from said exhaust gas sensor is stable or not based on the value of a switching function used in said sliding mode control process.

35. An air-fuel ratio control system according to claims 17-20, wherein said reference value variable setting comprises means for determining whether the output from said exhaust gas sensor is substantially converged to said target value or not, and holding said reference value to a predetermined value irrespective of said manipulated variable if the output from said exhaust gas sensor is not converged to said target value.