Fuel metering control system and cylinder air flow estimation method in internalcombustion engine

Abstract: Fuel metering control system in an internal combustion engine utilizing adaptive control having an intake manifold wall's fuel adherence plant. In the system, an actual air/fuel ratio in the individual cylinders is accurately estimated using an exhaust manifold model with an observer. Also, an actual cylinder air flow is estimated using a fluid model. Based on them, a desired cylinder fuel flow is determined by dividing the actual cylinder air flow by a desired air/fuel ratio and an actual cylinder fuel flow is determined by dividing the actual cylinder air flow by the estimated actual air/fuel ratio. The adaptive controller operates such that the actual cylinder fuel flow constantly coincides with the desired cylinder fuel flow. In an embodiment, in order to respond the change in wall adherence parameters, a compensator is connected in series with the wall adherence plant, a virtual plant incorporating the compensator is postulated and when the transfer characteristics of the virtual plant is other than 1 or thereabout, the adaptive controller is operated to have a transfer characteristics inverse thereto. At the same time, a method for estimating cylinder air flow inducted in the engine using the aforesaid fluid model is explained.

Description: BACKGROUND OF THE INVENTION Field of the Invention This invention relates to a system for controlling fuel metering in a multicylinder internal combustion engine according to the preamble part of claim 1. Description of the Prior Art During transient engine operation, a cylinder fuel flow is apt to be out of a desired value, and a lean or rich spike occurs in the air/fuel ratio. One cause for this is fuel transport delay caused by the adherence of fuel to the wall of the intake manifold etc. The behavior of the fuel transport delay changes depending on the operating states of the engine, initial manufacturing variance, and time-course changes of the intake manifold or the like owing to the adherence of deposits its wall. With a view to overcoming the problems caused by fuel transport delay, Japanese Laid-open Patent Publication Nos. 2(1990)-173,334 and 3(1991)-26,839 propose that fuel metering in an internal combustion engine be controlled by incorporating adaptive control in which a fuel adherence plant and a parameter adjuster are established such that the plant's output, i.e., an actual cylinder fuel flow coincides with a desired value even in a transient operating state of the engine. For adaptively compensating for fuel transport delay in a multicylinder internal combustion engine, however, it is indispensable to determine the air/fuel ratio at the individual cylinders with high precision so as to be able to estimate accurately the actual cylinder fuel flow inducted into the individual cylinders. Since the prior art controls proposed in the aforementioned references immediately use the air/fuel ratio measured at the exhaust gas confluence point for the whole cylinders, however, it is not able to estimate the actual cylinder fuel flow inducted into the individual cylinders with good accuracy. From U.S. 5,023,795 a system according to the preamble part of claim 1 is known. In this prior art system the fuel amount (wetting fuel amount) adhering to a wall of the air intake passage is estimated on a real-time basis in order to achieve the object of decreasing the fuel transport delay such that the fuel amount injected coincides with the desired (target) fuel amount. However, it is difficult to accurately estimate the wetting fuel amount on a real-time basis because of the effort of evaluation. From JP-A-1276334, a further system for controlling a fuel injection quantity is known, wherein a fuel behavior predicting means predicts the behavior of fuel using a model based on the detected engine speed, the fuel quantity evaporated from the fuel intake pipe wall, the fuel quantity flowing into the engine cylinder and the fuel quantity sprayed out over the fuel injector. The fuel behavior predicting means estimates - based on the predicted result - a parameter of said model. The fuel injection is thus adaptively controlled based on the estimated parameter. Thus, the adapted fuel injection quantity control according to JP-A-1267334 is also designed to estimate the wetting fuel amount on a real-time basis, which also leads to a high evaluation effort. A similar system is knwon from JP-A-1200040. Summary of the Invention It is therefore a primary object of the invention to provide a system for controlling fuel metering in a multi-cylinder internal combustion engine according to the preamble part of claim 1, wherein the actual cylinder fuel flow inducted into each cylinder can be determined with high accuracy and the fuel injection amount can be adaptively controlled on the basis of the so determined actual cylinder fuel flow. The operating states of the engine generating the fuel transport delay includes not only the states defined by engine coolant temperature, intake air temperature or the like that change relatively slow with respect to time, but also the state defined by manifold absolute pressure which varies rapidly. For example, when an accelerator pedal is depressed at a low engine speed, the manifold absolute pressure rises quickly, resulting rapid change in the fuel adherence condition. Since, however, the prior art control observes only the plant's input-output response, it is not able to follow up such a rapid change in the engine operating state. In other words, the actual fuel behavior finishes its change before it appears as the change in the plant's output. The prior art control, nevertheless, estimates the adherence parameter only when the plant's output changes and hence, leaves much to be improved in control response. Another object of the invention is therefore to overcome the aforesaid drawbacks of the prior art and to provide a system for controlling fuel metering in an internal combustion engine wherein the fuel behavior is observed at a real time such that the actual cylinder fuel flow follows a desired value with a better response according to the change in the fuel transport delay. The fuel metering control is usually encountered with a time lag problem. More specifically, it is not possible to immediately detect the air/fuel ratio of a mixture supplied into an engine cylinder. It can only be detected after the mixture burns and resultant exhaust gas reaches an air/fuel ratio sensor provided at the exhaust gas passage and emerges as a chemical-electric output signal. In addition, the lag is enlarged by a time required for fuel metering calculation and other factors such as a timing lag in outputting the calculated value. Even when the fuel metering is conducted through an adaptive control, it is not free from the problem. Thus, without accurately adjusting a timing between the input and output in the control, it is impossible to carry out a correction for the fuel behavior so as to determine a proper manipulated variable (control input) particularly at transient operating state of the engine. In the prior art control, however, although the air/fuel sensor's detection lag is observed, no further attention is made for adjusting individual input and output timings in the adaptive controller. Further object of the invention is therefore to overcome the aforesaid shortcoming of the prior art and to provide a system for controlling fuel metering in an internal combustion engine wherein no timing error occurs between a desired cylinder fuel flow and an output of a fuel adherence plant, i.e., an actual cylinder fuel flow so that an air/fuel ratio accurately converges on a desired value even at transient operating condition of the engine. Aside from the above, various methods have been proposed for measuring or estimating air flow drawn in an engine cylinder including the method for directly measuring the mass air flow or the so-called speed density method which estimates it through manifold absolute pressure. Both methods are, however, not free from the influence from engine transient operating condition, sensor's initial manufacturing variance or degradation in sensor's service life. In view of the above, there are proposed techniques to measure or estimate air flow by Japanese Laid-open Patent Publication No. 2(1990)-5745 and United States Patent No. 4,446,523 utilizing a fluid dynamic model. However, since the former technique proposed by the Japanese reference predicts pressure in the air intake passage and does not detect it directly, it was disadvantageous in accuracy. Further, since the former technique utilizes a recurrence formula, if the pressure be erroneously estimated, the error will then be accumulated and further enlarged. The latter technique proposed by the US patent relates to a mass air flow meter for merely measuring the mass air flow rate passing through a throttle plate and is silent to estimating an actual cylinder air flow. Still further object of the invention is therefore to overcome the aforesaid drawbacks of the prior art and to provide a system for estimating cylinder air flow wherein estimation accuracy is enhanced by directly detecting the pressure in the air flow passage and even if an error occurs, the error will not influence on the next estimation. For solving at least the primary object mentioned above the BRIEF DESCRIPTION OF THE DRAWINGS These and other objects and advantages of the invention will be more apparent from the following description and drawings, in which: Figure 1 is an overall block diagram showing a fuel metering control system according to the invention; Figure 2 is a block diagram focussing on a fuel metering control block redrawn from that illustrated in Figure 1; Figure 3 is a block diagram showing a wall adherence plant referred to in Figure 2; Figure 4 is a block diagram showing that a Model Reference Adaptive Control System is applied for the wall adherence compensation illustrated in Figure 2; Figure 5 is a block diagram showing rearranged configuration illustrated in Figure 4; Figure 6 is a view showing simulation results conducted on the configuration of Figure 5; Figure 7 is a view showing verification conducted on the data of Figure 6; Figure 8 is a block diagram illustrating a configuration provided with dead time in the configuration of Figure 5; Figure 9 is a view showing simulation results conducted on the configuration of Figure 8; Figure 10 is a view showing simulation results showing verification conducted on the data of Figure 9; Figure 11 is a view showing simulation results conducted on the configuration of Figure 5 using a constant gain method; Figure 12 is a view similar to Figure 11, but using a decreasing gain method; Figure 13 is a view similar to Figure 11, but using a variable gain method; Figure 14 is a view similar to Figure 11, but using a constant trace method; Figure 15 is a view showing an air intake system model for estimating cylinder air flow to be used in the fuel metering control illustrated in Figure 1; Figure 16 is a view showing simulation results to the actual cylinder air flow estimated using the model of Figure 15; Figure 17 is a view showing a testing apparatus to be used for the estimation; Figure 18 is a view showing test results using the apparatus of Figure 17; Figure 19 is a view showing test results for identifying the flow rate coefficient with respect to throttle opening; Figure 20 is a view showing estimated values obtained by using the identification results of Figure 19 and illustrated in contrast with measured values; Figure 21 is a view showing values obtained based on the model of Figure 15 and illustrated in contrast with measured values; Figure 22 is a block diagram showing calculation of a throttle effective opening area using the flow rate coefficient etc.; Figure 23 is a view showing the characteristics of a mapped data of a coefficient including the flow rate coefficient set with respect to manifold absolute pressure and throttle opening; Figure 24 is a view showing a control error with respect to throttle opening; Figure 25 is a view showing a control with respect to pressures at upstream and downstream of a throttle valve; Figure 26 is a block diagram showing an air/fuel ratio estimation used in the fuel metering control system of Figure 1; Figure 27 is a block diagram showing a detailed configuration of an EXMN PLANT illustrated in Figure 1; Figure 28 is a block diagram showing the configuration of Figure 27 incorporated with an observer; Figure 29 is a view showing that the fuel metering control system is applied to an actual engine; Figure 30 is a block diagram showing the details of a control unit illustrated in Figure 29; Figure 31 is a flow chart showing the operation of the system of Figure 29; Figure 32 is a block diagram showing a second embodiment of the invention; Figure 33 is a view, similar to Figure 32, but showing a third embodiment of the invention; Figure 34 is a view, similar to Figure 2, but showing a fourth embodiment of the invention; Figure 35 is a view, similar to Figure 8, but showing the fourth embodiment of the invention; and Figure 36 is a view, similar to Figure 22, but showing a fifth embodiment of the invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Figure 1 is an overall block diagram of a fuel metering control system according to the present invention utilizing adaptive control. The control system includes a MAP block comprising predetermined characteristics prepared as a mapped data in a computer memory from which a desired cylinder fuel flow Ti is retrieved using engine speed Ne, manifold absolute pressure Pb and the like as address data, a Gair model block for estimating the dynamic behavior of an actual cylinder air flow Gair from throttle opening TH, manifold absolute pressure Pb etc., and an A/F observer block for estimating an air/fuel ratio of the individual cylinders from the air/fuel ratio measured at the exhaust gas confluence point, and a fuel metering control block for determining an fuel injection amount Tout. In this configuration, the cylinder fuel flow Gfuel at each instant (combustion cycle) is estimated from the estimated (actual) cylinder air flow Gair and air/fuel ratio A/F, and the parameters of the fuel metering control block are adjusted to determine the fuel injection amount Tout such that the actual cylinder fuel flow Gfuel coincides with the desired cylinder fuel flow Ti. Here, the word "mapped" data means a data stored in a computer memory with respect to two parameters. Similarly a word "table" means a look-up table stored in the memory with respect to a single parameter. These will now be explained in detail. The fuel metering control block will be explained first. As regards fuel metering control, Figure 1 can be redrawn as shown in Figure 2. The input parameters are: (1) Desired cylinder fuel flow Ti Value obtained by dividing an actual cylinder air flow Gair estimated using the inputs from the sensors by a desired A/F ratio. (The calculation of the actual cylinder air flow Gair will be explained later.) (2) Actual cylinder fuel flow Gfuel Value obtained by dividing the actual cylinder air flow Gair by an actual air/fuel ratio at the same cylinder calculated from the value measured by an air/fuel ratio sensor. (The calculation of the actual air/fuel ratio at the individual cylinders will be explained later.) (3) Others Various measured and estimated values required by a wall adherence correction compensator (e.g., engine coolant temperature Tw, manifold absolute pressure Pb, engine speed Ne etc.) Specifically, as is clear from the foregoing, the actual cylinder air flow Gair in a combustion cycle at a given time (k-n) is obtained and divided by the desired air/fuel ratio A/F (k-n) to determine the desired cylinder fuel flow Ti (k-n). In addition, the actual cylinder air flow Gair (k-n) in the same combustion cycle is divided by the measured and calculated air/fuel ratio A/F at the same cylinder to determine the actual cylinder fuel flow Gfuel (k-n). Then, a dynamic compensator in an adaptive controller is adjusted so that the actual cylinder fuel flow Gfuel (k-n) constantly coincides with the desired cylinder fuel flow Ti (k-n), whereby the manipulated variable (fuel injection amount) Tout is determined. In order to respond promptly to the aforesaid adherence parameters' change, the aforesaid wall adherence correction compensator is inserted ahead of a wall adherence plant. The transfer function of the wall adherence correction compensator is the inverse of that of the wall adherence plant. The adherence parameters of the wall adherence correction compensator are retrieved from a mapped data prepared beforehand on the basis of their correspondence with the engine operating states. If the adherence parameters of the wall adherence correction compensator are equal to adherence parameters of an actual engine, the transfer function of the two as seen from the outside is 1, namely the product of the transfer functions of the plant and the compensator is 1. Since this means that the actual cylinder fuel flow equals the desired cylinder fuel flow, perfect correction should be obtained. In fact, however, the adherence parameters generally vary complexly depending on the engine operating states, making it difficult to realize perfect coincidence. Moreover, the actual engine experiences initial manufacturing variance and time-course changes owing to the adherence of deposits and the like. If these factors should cause the adherence parameters to vary between the compensator and the actual engine, the value of the transfer function will become something other than 1 or thereabout, i.e, 1.1, 1.2, 0.9, 0.8,.. Since time response therefore occurs, the desired cylinder fuel flow and the actual cylinder fuel flow will not be equal. In view of the above, therefore, a virtual plant incorporating the adherence correction compensator is postulated and when the transfer characteristic of the virtual plant is other than 1 or thereabout, the adaptive controller is operated to have a transfer characteristic inverse thereto. The desired cylinder fuel flow is input to the adaptive controller as a desired value and adaptive parameters are used which vary so that the actual cylinder fuel flow, namely the output of the virtual plant, coincides with the desired value. The parameters of the adaptive controller are calculated by an adaptive parameter adjuster (identifier). The adaptive parameter adjuster (identifier) uses input/output values including past values input to the virtual plant. The adaptive controller also functions to absorb errors in the estimated (actual) cylinder air flow. In other words, since in the end the adaptive parameters are adjusted so that coincidence is constantly maintained between the actual cylinder fuel flow obtained by dividing the cylinder air flow by the measured air/fuel ratio and the desired cylinder fuel flow obtained by dividing the cylinder air flow by the desired air/fuel ratio, any error in the estimated (actual) cylinder air flow can therefore be absorbed. This will be explained in more detail. As the wall adherence plant the first-order model such as expressed in Eq. 1 is used. Here, two parameters are used. Qin(k) = A · Qout(k) + B · Qt(k - 1) Qt(k) = (1 - B) · Qt(k - 1) + (1 - A) · Qout(k) where Qt(k): Wall adherence amount A (0≤A≤1: Direct ratio (cylinder flow ratio) B (0≤B≤1): Carry-off ratio (vaporization ratio) Qin(k): Actual cylinder fuel flow Qout(k): Injector's injection amount Expressed as a discrete transfer function it becomes as shown in Eq. 2. Shown in block diagram it becomes as shown in Figure 3. Qin(z) Qout(z) = Az - (A - B) z - (1 - B) The transfer function of the wall adherence correction compensator is represented by Eq. 3. As mentioned earlier, it is the inverse of the transfer function of the wall adherence plant. F(z) = z - (1 - B ) A z - ( A - B ) The characteristics of the aforesaid direct ratio A and the carry-off ratio B (here both expressed with a circumflex) of the wall adherence correction compensator are stored as the mapped data in advance as functions of the engine operating states, as earlier mentioned, such as engine coolant temperature Tw, manifold absolute pressure Pb, engine speed Ne and the like and are retrieved using the values of these. (In this specification, a value with the circumflex represents an estimated value.) The adaptive controller will then be explained. Among the conditions required of the wall adherence correction are that it constantly work to reduce the transport delay and that it be able to follow variation in the A and B terms in the equations. A well-known system for achieving adaptive control that follows such a time-varying plant is MRACS (Model Reference Adaptive Control System). The configuration when MRACS is applied for wall adherence compensation is shown in Figure 4. In this case, a priori model (model reference) can be taken near the center value of the time-varying plant or it can be taken so as to facilitate control of the wall adherence correction compensator. Since MRACS is effective only for a plant with dead time (delay time), dead time is apparently inserted by delaying input to the adherence plant by one cycle, thus constituting the virtual plant (the word "virtual" is appended to the inserted blocks). It will be noted that the virtual adherence correction compensator and the virtual model reference are connected in series. Therefore, since their transfer functions are the inverse of each other, they can be canceled. This results in a z d = z (d=1) block and D (z -1 ) remaining immediately after the virtual model reference. However, since z is a transfer function which outputs a future value, it cannot exist as it is. Therefore, by defining D (z -1 ) as D (z -1 ) = z -1 , the two can be canceled. Although D (z -1 ) is normally defined as D (z -1 ) = 1 + d 1 z -1 + ... +d n z -n , defining it as D (z -1 ) = z -1 does not cause a problem regarding stability. Thus, rearranging Figure 4 gives the configuration of Figure 5. (As a result, the adaptive controller becomes a controller which handles a regulator problem and is modified to an STR (Self-Tuning Regulator). The adaptive controller receives the coefficient vectors identified by a parameter identifier, thus constituting a feedback compensator. However, since this operation is known (see, for example, the detailed explanation from page 28 to 41 in an article entitled "Digital Adaptive Control" in a magazine "Computrol " No. 27), it will not be explained further here. Figure 6 shows the responses obtained by simulation with the illustrated configuration. From this figure it can be seen that the MRACS parameter identifier operates normally in the aforesaid configuration, but that the behavior of the air/fuel ratio remains jagged. If, in order to verify this microscopically, a desired cylinder fuel flow such as shown in Figure 7(a) is input, the plant output and the air/fuel ratio become as shown in Figures 7(b) and 7(c). It will be noted that the plant output is delayed by one cycle. This delay occurs because the virtual plant was constituted by insertion of dead time. It will also be noted that a lean spike occurs in the air/fuel ratio owing to the fact that during transient engine operation a one-cycle time difference arises between the desired cylinder fuel flow and the plant. Since inserting the dead time z -d before the plant and inserting it after the plant are equivalent when the virtual plant is viewed from the outside, it is here inserted after. The plant output y' (k) is extracted immediately after the plant, i.e. between the plant and the dead time z -d , and the virtual plant output y (k) required by the parameter identifier is extracted from after the dead time z -d . This arrangement ensures that no dead time is present in the path from the input r (k) and the plant output y' (k) and enables the parameter identifier to use the virtual plant output y (k) including the dead time. The configuration is shown in Figure 8. Figure 9 shows the simulation results for the configuration of Figure 8. As shown in Figure 9(c), after convergence, the actual cylinder fuel flow becomes substantially equal to the desired cylinder fuel flow. The air/fuel ratio at this time stays flat in the vicinity of 14.7. Moreover, when the microscopic response after completion of identification is verified in the same scale, it becomes as shown in Figure 10 (in which, for comparison, the response before implementing the dead time is shown by broken line curves). It will be noted that once the identification is completed the response after implementing the dead time is characterized by a very flat air/fuel ratio. This insertion of dead time is not limited to that explained in the foregoing. As shown by the phantom blocks in Figure 8, it can be appropriately inserted in the input and/or output in correspondence with the order of the plant output. This will later be referred in a fourth embodiment. Turning next to the parameter identification laws, when the method proposed by I. D. Landau et al is used in the parameter identifier shown in Figure 4, the gain matrix is represented by Γ(k)= 1 λ1(k) [Γ(k-1)- λ2(k)Γ(k-1)ζ(k-1)ζ T (k-1)Γ(k-1) λ1(k)+λ2(k)ζ T (k-1)Γ(k-1)ζ(k-1) ] Where: 0<λ1(k)≤1, 0<λ2(k)<2, Γ(0)>0 The specific parameter identification laws are determined by how lambda 1(k) and lambda 2(k) are chosen. The typical MRACS identification laws fall in four categories: constant gain method, decreasing gain method (including the method of least squares), variable gain method (including the method of weighted least squares) and the constant trace method. Based on the configuration of Figure 4, simulation was conducted with respect to each under the following conditions. Specifically, the time-varying plant was used since it is apt to be the one involved in application to the actual engine. Figures 11 to 14 show the results of the simulation. As will be understood from these simulation results, in the case of a time-varying plant, when the constant gain method is used (Figure 11) the plant output value exhibits intense hunting centered on the desired value. The hunting is particularly pronounced when the desired value is changing (during transient engine operation). During transient engine operation the difference between the model reference and the plant output value, which is the desired value of the model reference output, becomes large and, therefore, the MRACS parameter identifier attempts to make a sudden great change in the parameter values. As a result, if, for example, the plant variation is too fast, overshooting occurs and causes hunting. In the case of the decreasing gain method (Figure 12), the variable gain method (Figure 13) and the constant trace method (Figure 14), on the other hand, the plant output faithfully follows the model reference constituting the desired value. Although it oscillates in spots, it can be seen to converge on the desired value. Oscillation of this degree can be suppressed by adjusting the parameters, e.g. by varying the gain matrix values or D (z -1 ), without sacrificing the convergence speed. Thus the last-mentioned three identification laws enable faster convergence speed than the constant gain method and can provide faithful following even if the plant is time variable. The estimation of the actual cylinder air flow Gair will now be explained. As was pointed out earlier, for accurately determining the actual cylinder fuel flow Gfuel it is necessary to determine the air mass flow rate with high precision. Conventional methods available for this include the method of measuring mass flow rate of air directly and the speed density method of indirect estimation from the manifold absolute pressure. However, since these known methods operate on the principle of retrieving the air flow rate from a mapped data prepared using parameters having a high degree of correlation with the cylinder air flow, they are powerless with respect to changes of the parameters not taken into account during preparing the mapped data and, therefore, lack toughness in respect of deterioration, variance and aging. Moreover, since preparing the mapped data can intrinsically be conducted only with respect to steady-state engine operating conditions, it cannot express transient engine operating states. This means that for transient engine operation there is no choice other than to have set the cylinder air flow in advance by an engineer's volition. In this invention, therefore, there is applied a fluid dynamic model capable of reflecting variation in the air flow under various air intake system conditions. Notwithstanding that the measurement is more indirect than in the conventional methods, its accuracy is higher owing to the fact that preparation of the mapped data or setting the data by an engineer's volition is eliminated. More specifically, the throttle is viewed as an orifice, the mass of air passing through the throttle is estimated using a fluid dynamic model of the vicinity of the throttle, and the actual air mass flow rate past the throttle is dynamically estimated with consideration to the chamber charging delay. This will now be explained. If the throttle is viewed as an orifice as shown in an air intake system model of Figure 15, it is possible from Eq. 5 (Bernoulli's equation), Eq. 6 (equation of continuity) and Eq. 7 (relational equation of adiabatic process) to derive Eq. 8, which is a standard orifice equation for compressible fluid flow. It is thus possible to determine the air mass flow rate Gth through the throttle valve per unit time. v 1 2 2 + κ κ-1 · P 1 ρ 1 = v 2 2 2 + κ κ-1 · P 2 ρ 1 where the flow is assumed to be the adiabatic process, and P 1 : Absolute pressure on upstream side P 2 : Absolute pressure on downstream side ρ 1 : Air density on upstream side ρ 2 : Air density on downstream side v 1 : Flow velocity on upstream side v 2 : Flow velocity on downstream side κ: Ratio of specific heats ρ 1 · v 1 · A up = ρ 2 · v 2 · S where: A up : Flow passage area on upstream side S: Throttle projection area [= f(TH)] P 1 ρ 1 κ = P 2 ρ 2 κ Gth = ε · ρ 1 · α · S · 2g·(P 1 - P 2 ) γ 1 where: g: Gravitational acceleration γ 1 : Air specific weight on upstream side (=ρ 1 ·g) α: Flow rate coefficient (coefficient of discharge) ε: Correction coefficient (expansion factor of gas) Next, the mass of air in the chamber is calculated from Eq. 9, which is based on the ideal-gas law. The term "chamber" is used here to mean not only the part corresponding to the so-called surge tank but all portions between immediately downstream of the throttle and the intake port. Gb(k) = V RT · P(k) where: V: Chamber volume T: Air temperature R: Gas constant P: Pressure Therefore, the change delta Gb in the mass of air in the chamber in the current cycle can be obtained from the pressure change using Eq. 10. Specifically, under steady-state engine operating conditions it holds that Gth = Gair . On the other hand, under transient engine operation condition the reason that the manifold absolute pressure rises when, for example, the throttle valve is opened suddenly is that the chamber is full of air. This means that if the mass of air charged in the chamber and the mass of air that passed through the throttle valve are known, the mass of air inducted into the cylinder can be found. In other words, if it is assumed that, as is only natural, the mass of air charged in the chamber is not inducted into the cylinder, then the actual cylinder air flow Gair per time unit delta T can be expressed by Eq. 11, whereby it becomes possible to estimate the dynamic behavior of the actual cylinder air flow. Figure 16 shows the results of simulation using this method. Gair = Gth · ΔT - ΔGb The results of a test regarding the foregoing will be set out. The testing apparatus used is shown schematically in Figure 17. The test was conducted by maintaining the throttle opening constant and measuring the change in pressures at upstream and downstream of the throttle when the air flow was varied. Regarding the upstream side of the throttle, the test was conducted for ten different throttle openings. Among these, the results for a throttle opening of 31.6 degrees are shown in Figure 18. From these and the other test results it could be concluded as follows. (1) Moving downstream from the throttle, the pressure drops at a distance of 1D to 2D (D : throttle bore diameter), recovers at 3D to 4D and then gradually decreases from thereon (owing to the contraction, swirling and separation of the flow caused by the throttle valve). (2) It is necessary to calculate the throttle-pass air flow using the recovered pressure value because the pressure difference at upstream and downstream of the throttle appears larger than actual when measured in the pressure drop region. It was further found that the pressure drops just before the throttle on the upstream side. From the foregoing, it was concluded to be preferable to measure the pressure Pthdown (P 2 in Eq. 5) downstream of the throttle at a position in the pressure recovery region (i.e. about 3D (ideally 3D - 4D) from the throttle valve) and to measure the pressure Pthup (P 1 in Eq. 5) upstream of the throttle at a position which is as close to the throttle valve as possible but which is unaffected by the throttle valve (i.e. about 1D or more) from the throttle valve. Since in this sense the pressure downstream of the throttle can be assumed equal to the chamber (surge tank) pressure, as will be explained further later one possible arrangement is to define the detection value of a pressure sensor installed in the surge tank as the pressure Pthdown downstream of the throttle. Taking the flow rate coefficient α and the correction coefficient epsilon to be unknown in Eq. 8, the product of the flow rate coefficient α and the correction coefficient epsilon was identified by the foregoing test (the parameter rho 1 was calculated from the barometric condition at the test). The identification was conducted by using the measured pressures across the throttle to calculate the mass flow rate passing the throttle Gth per unit time (the initial value of which was appropriately set), comparing the calculated value with the measured value, varying the product to bring the calculated and measured values into coincidence, repeating the foregoing to obtain the value involving minimum error and defining this value as the flow rate coefficient. The relationship between the product identified by this method and the throttle opening is shown in Figure 19. The values estimated using the identified product are compared with the measured values in Figure 20 (only for a throttle opening of 31.6 degrees). Figure 21 shows a comparison between measured values and the values calculated by simulation using the product of the flow rate coefficient and correction coefficient obtained in the foregoing manner and the values measured at positions 4D downstream and 1D upstream of the throttle valve. This figure shows the data obtained when the throttle opening was varied between 7 and 20 degrees. The value Pb is the value measured by manifold absolute pressure sensor and the value Gth is the value measured by an air flow meter. In the data illustrated in Figure 21, the values obtained through simulation almost coincided with the measured values. Continuing tests, however, it was found that they were not always equal at all situations. Here, when rewriting Eq. 8, Eq. 12 will be obtained. In Eq. 12, representing the product of the flow rate coefficient α and correction coefficient epsilon by a coefficient C, it was considered, as mentioned earlier, that the coefficient C could solely be determined from the profile of the throttle and dependent on its opening. After conducting tests repeatedly, however, it was found that the coefficient C could not be identified from the throttle opening alone, since a laminar flow or a turbulent flow happened depending on the rate of flow and the state of flow in the vicinity of wall changed by the occurrence of separation or swirling. Namely, it was confirmed that the coefficient C depended, not only on the throttle opening, but also on the flow rate. However, since the coefficient C must be determined in order to determine the flow rate itself, it is not possible to use, as a matter of fact, the flow rate as an input parameter. Instead, therefore, engine load, i.e. manifold absolute parameter Pb was used as a parameter indicative of the state of flow and good results were obtained. And, as illustrated in Eq. 12, the coefficient C thus obtained was multiplied to a throttle projection area S to determine throttle effective opening area A. As a result, it becomes possible to determine the throttle effective opening area A at all engine operating states with accuracy and to estimate the actual cylinder fuel flow precisely. Figure 22 shows the configuration. Here, the throttle projection area is an area generated when throttle valve be projected in the direction parallel to the throttle bore's longitudinal direction. It should be noted that the characteristics of the coefficient C with respect to throttle opening TH and manifold absolute pressure Pb are determined in advance through experiments and is prepared as a mapped data in a computer memory as illustrated in Figure 23. And at the time of preparing the mapped data, an interval between adjacent lattice points should be set to be decreasing with decreasing throttle opening. This is because the change of the coefficient C to the change of the throttle opening becomes large with decreasing throttle opening. Moreover, as illustrated in the same figure, the coefficient C should be set to be at or below 1.0. That is, it is difficult to imagine in the sense of physics that the effective opening area becomes greater than the projection area and the effective opening area is assumed to be increasing monotonously relative to the throttle opening. In addition, since the flow rate coefficient α and correction coefficient epsilon are both found to be related to manifold absolute pressure, they are treated as a whole as explained before. This brings a side effect that an error, if happened, will be lessened when compared with a case in which they are determined separately. Furthermore, as shown in Eq. 12, the pressures P1 and P2 at upstream and downstream of the throttle are represented by atmospheric (barometric) pressure Pa and manifold absolute pressure Pb. And answers in the square root using the pressures are calculated in advance and stored as a mapped data similarly to that shown in Figure 23. Moreover, as illustrated in Figure 22, the throttle projection area S is obtained through a detected throttle opening TH and the coefficient C is multiplied thereto to obtain the throttle opening area A. The relationship between the throttle opening TH and the projection area S is accordingly determined in advance through experiments and stored in a table in a computer memory. The relationship with the sensor's resolving power will next be discussed. Figure 24 is based on measured data, the vertical axis representing the control error for a given measurement error and the horizontal axis representing throttle opening. The figure shows that the control error with respect to a given measurement error increases with decreasing throttle opening. It is therefore preferable to use a sensor whose measurement error decreases with decreasing throttle opening, i.e, one whose resolving power increases with decreasing throttle opening. Figure 25 is based on measured data, the vertical axis again representing control error and the horizontal axis representing the ratio of the pressures on opposite sides of the throttle valve. It will be understood that it is preferable to use a manifold absolute pressure sensor whose resolving power increases with increasing load (toward the atmospheric pressure side indicated by 1 in the figure). In application to an actual engine, therefore, both, or at least one, of the throttle opening sensor and the manifold absolute pressure sensor should exhibit such preferable resolving power characteristics. Some additional comments can be made regarding measurement of the air flow rate. First, the air flow rate is fixed at a prescribed value (e.g. 0.528) when the ratio of the pressures on opposite sides of the throttle valve is lower than a prescribed value since the flow velocity is equal to the sound velocity at such times. Further, for enhancing the calculation accuracy, the intake air temperature sensor is located near the throttle valve on the upstream side. In addition, it is preferable to install a hygrometer and use its output for correcting the air specific weight in Eq. 8. It should further be noted that although the pressure Pb is detected in terms of absolute pressure, it is alternatively possible to detect by gauge pressure. Further, the coefficient C can be determined from the throttle opening TH and a deviation (Pa - Pb) between the manifold absolute pressure Pb and the atmospheric pressure Pa or their ratio (Pb/Pa). Furthermore, the coefficient C may be determined from the throttle opening and any other environmental factor. The detection of the air/fuel ratios at the individual cylinders will now be explained. From the points of cost and durability, multicylinder internal combustion engines are generally equipped with only a single air/fuel ratio sensor mounted at the exhaust gas confluence point. This makes it necessary to determine the air/fuel ratios at the individual cylinders from the air/fuel ratio at the confluence point. In this invention, therefore, the air/fuel ratio behavior at the convergence point is modeled and the air/fuel ratios at the individual cylinders are estimated by numerical calculation from the air/fuel ratio at the convergence point. Here, the air/fuel ratio sensor indicates not the so-called O 2 sensor, but a sensor which can detects an air/fuel ratio varying linearly with the oxygen concentration of the exhaust gas over a broad range extending from the lean direction to the rich direction. As this air-fuel ratio is explained in detail in the assignee's earlier Japanese patent application (Japanese Patent Application No. 3(1991)-169,456 filed June 14, 1991), it will not be discussed further here. First, the response delay of the air/fuel ratio sensor is approximately modeled as a first-order delay, the state equation for this is obtained and the result is discretized for the period delta T, giving Eq. 13. In this equation, LAF stands for the air/fuel sensor output and A/F for the input air/fuel ratio. LAF(k+1) = α LAF(k) + (1 - α )A/F(k) where: α=1+αΔT+(1/2!)α 2 ΔT 2 +(1/3!)α 3 ΔT 3 +(1/4!)α 4 ΔT 4 Z-transforming Eq. 13 to express it as a transfer function gives Eq. 14. In other words, as shown in Figure 26, the air/fuel ratio in the preceding cycle (time k-1) can be obtained by multiplying the sensor output LAF in the current cycle (time k) by the inverse transfer function of Eq. 14. t(z) = (1 - α )/(Z - α ) An explanation will now be given regarding the method used to separate and extract the air/fuel ratios at the individual cylinders from the air/fuel ratio corrected for delay in the foregoing manner. First, the internal combustion engine exhaust system is modeled as shown in Figure 27. This model corresponds to EXMN PLANT in Figure 1. It should be noted that fuel is a controlled variable in this model (plant) so that a fuel/air ratio F/A is used here. The inventors found that the air-fuel ratio at the exhaust confluence point could be expressed as an average weighted to reflect the time-based contribution of the air-fuel ratios of the individual cylinders. That is to say, it can be expressed in the manner of Eq. 15. [Confluence point F/A](k) = C1 · [F/A](k-3) + C2 · [F/A](k-2) + C3 · [F/A](k-1) + C4 · [F/A](k) Eq. 16 expresses the air/fuel ratios at the individual cylinders in the form of a recurrence formula. Since the input U (k) is unknown, if, assuming a four-cylinder engine, a recurrence formula is written for reproducing the air/fuel ratio once every 4 TDC (top dead center), the result becomes Eq. 17. The problem is thus reduced to the ordinary state equation such as expressed by Eq. 18. Therefore, if the time-based degree of contribution C is known, it is possible, by designing a Kalman filter and configuring the observer shown in Figure 28, to estimate X (k) at each instant from Y (k). In other words, an appropriate gain matrix is established for a state equation such as the foregoing, and consideration is given to X circumflex (k) of an equation such as Eq. 19. X (k+1) = (A - KC) · X (k) + K·Y(k) If (A - KC) is a stable matrix, X circumflex (k) becomes X (k), and X (k) (air/fuel ratios at individual cylinders) can be estimated from Y (k) (air/fuel ratio at exhaust confluence point). As this was explained in detail in the assignee's earlier Japanese Patent Application No. 3(1991)-359,340 filed December 27, 1991 (and filed in the Unite States on December 24, 1992 under the number of 997,769 and in the EPO on December 29, 1992 under the number of 92 31 1841.8), it will not be discussed further here. A specific example of the application of the foregoing to an actual engine will now be explained. An overall view of the example is shown in Figure 29. Reference numeral 10 in this figure designates an internal combustion engine. Air drawn in through an air cleaner 14 mounted on the far end of an air intake path 12 is supplied to first to fourth cylinders through a surge tank (chamber) 18 and an intake manifold 20 while the flow thereof is adjusted by a throttle valve 16. An injector 22 for injecting fuel is installed in the vicinity of the intake valve (not shown) of each cylinder. The injected fuel mixes with the intake air to form an air-fuel mixture that is ignited in the associated cylinder by a spark plug (not shown). The resulting combustion of the air-fuel mixture drives down a piston (not shown). The exhaust gas produced by the combustion is discharged through an exhaust valve (not shown) into an exhaust manifold 24, from where it passes through an exhaust pipe 26 to a three-way catalytic converter 28 where it is removed of noxious components before being discharged to the exterior. A crank angle sensor 34 for detecting the piston crank angles is provided in a distributor (not shown) of the internal combustion engine 10, a throttle position sensor 36 is provided for detecting the degree of opening TH of the throttle valve 16, and a manifold absolute pressure sensor 38 is provided for detecting the absolute pressure Pb of the intake air downstream of the throttle valve 16. On the upstream side of the throttle valve 16 are provided an atmospheric pressure sensor 40 for detecting the atmospheric (barometric) pressure Pa, an intake air temperature sensor 42 for detecting the temperature of the intake air and a hygrometer 44 for detecting the humidity of the intake air. The aforesaid air/fuel ratio sensor 46 comprising an oxygen concentration detector is provided in the exhaust system at a point downstream of the exhaust manifold 24 and upstream of a three-way catalytic converter 28, where it detects the air/fuel ratio of the exhaust gas. The outputs of the sensor 34 etc. are sent to a control unit 50. In the foregoing configuration, the atmospheric pressure sensor 40 for detecting the pressure upstream of the throttle is disposed at a position apart from the throttle valve 16 by at least 1D (D : diameter of the intake passage 12) and the manifold absolute pressure sensor 38 for detecting the pressure downstream of the throttle is disposed in the surge tank 18 and the surge tank 18 is disposed at least 3D apart from the throttle valve 16. The intake air temperature sensor 42 and the hygrometer 44 are disposed as close as possible to the throttle valve 16. The resolving power of the throttle position sensor 36 is at least 0.01 degree and that of the manifold absolute pressure sensor 38 at least 0.1 mmHg. Details of the control unit 50 are shown in the block diagram of Figure 30. The output of the air/fuel ratio sensor 46 is received by a detection circuit 52 of the control unit 50, where it is subjected to appropriate linearization processing to obtain an air/fuel ratio (A/F) characterized in that it varies linearly with the oxygen concentration of the exhaust gas over a broad range extending from the lean side to the rich side, as was referred to earlier. The output of the detection circuit 52 is forwarded through an A/D (analog/digital) converter 54 to a microcomputer comprising a CPU (central processing unit) 56, a ROM (read-only memory) 58 and a RAM (random access memory) 60 and is stored in the RAM 58. Similarly, the analogue outputs of the throttle position sensor 36 etc. are input to the microcomputer through a level converter 62, a multiplexer 64 and a second A/D converter 66, while the output of the crank angle sensor 34 is shaped by a waveform shaper 68 and has its output value counted by a counter 70, the result of the count being input to the microcomputer. In accordance with commands stored in the ROM 58, the CPU 56 of the microcomputer computes control values in accordance with the adaptive control method explained earlier and drives the injectors 22 of the individual cylinders via a drive circuit 72. The operation of the control apparatus of Figure 30 will now be explained with reference to the flow chart of Figure 31. The engine speed Ne detected by the crank angle sensor 34 is read in step S10. Control then passes to step S12 in which the atmospheric pressure Pa (same as pressure Pthup or P1 upstream of the throttle), the manifold absolute pressure Pb (same as pressure Pthdown or P2 downstream of the throttle), the throttle opening TH, the air/fuel ratio A/F and the like detected by the atmospheric pressure sensor 40 etc. are read. Program then passes to step S14 in which discrimination is made as to whether or not the engine is cranking, and if it is not, to step S16 in which a discrimination is made as to whether or not the fuel supply has been cut off. If the result of the discrimination is negative, program passes to step S18 in which the desired cylinder fuel flow Ti is calculated by map retrieval as shown in Figure 1 using the engine speed Ne and the manifold absolute pressure Pb as address data, and to step S20 in which the fuel injection amount Tout is calculated in terms of injector's injection period in accordance with the basic mode equation. (The basic mode is a well-known method that does not use the aforesaid adaptive control.) Program then passes to step S22 in which a discrimination is made as to whether or not activation of the air/fuel ratio sensor 46 has been completed, and if it has, to step S24 in which the air/fuel ratios of the individual cylinders are estimated by the method described in the foregoing, to step S26 in which the actual cylinder air flow Gair is estimated, to step S28 in which the actual cylinder fuel flow Gfuel is estimated, to step S30 in which the fuel injection amount Tout is finally determined in accordance with the aforesaid adaptive control, and to step S32 in which the value Tout is output to the injector 22 of the associated cylinder through the drive circuit 72. When it is found in step S14 that the engine is cranking, program passes through steps S34 and S36 for calculating the start mode control value. When step S16 finds that the fuel supply has been cut off, program passes to step S38 in which the value Tout is set to zero. If step S22 finds that the sensor has not been activated, program jumps directly to step S32 and the injector is driven by the basic mode control value. In the foregoing configuration, the actual cylinder fuel flow is estimated with high precision based on the estimated air/fuel ratio at the individual cylinders and the parameters of the controller are adaptively controlled so as to make the actual cylinder fuel flow coincide with the desired value. As a result, it is possible to achieve high-precision adaptive control. In addition, since a compensator with a transfer coefficient that is the inverse of that of the fuel adherence plant is connected in series with the fuel adherence plant, adaptive control to the desired value can be achieved while closely following any variation in the adherence state even in cases where the variation is due to a factor which varies rapidly with time such as the manifold absolute pressure. What is more, since a virtual plant incorporating the adherence compensator is postulated and when the transfer characteristic of the virtual plant is other than 1 or thereabout the adaptive controller is operated to have the inverse transfer characteristic, adaptive control that realizes the desired value can be achieved while closely following any variation that may occur owing to deviation of the preset characteristics from the actual characteristics as a result of aging or the like. Although the invention was explained with reference to the configuration of Figure 1, this is not the only configuration to which it can be applied. Figure 32 shows a second embodiment of the invention. The configuration of the second embodiment does not have the Gair model block for estimating the dynamic behavior of the actual cylinder air flow but instead estimates the actual cylinder air flow Gair by multiplying the mapped value by the stoichiometric air/fuel ratio 14.7 and absorbs the intake system behavior by conducting adaptive control. In other words, as was explained earlier, even error in the estimated actual cylinder air flow can be absorbed. Figure 33 shows a third embodiment of the invention, wherein the desired cylinder fuel flow Ti is not stored as a mapped data but is decided by multiplying the actual cylinder air flow Gair estimated by the Gair model block by 1/14.7. Figure 34 and 35 show a fourth embodiment of the invention. As illustrated in Figure 34, the wall adherence correction compensator is omitted in the configuration of the fourth embodiment in contrast to that in the first embodiment shown in Figure 2. With the arrangement, however, when the transfer characteristic of the virtual plant becomes other than 1 or thereabout, the adaptive controller is also operated such that the transfer characteristic of the virtual plant and the adaptive controller becomes, as a whole, 1 or thereabout, i.e., the adaptive controller operates to have a transfer characteristics inverse thereto. A second characteristic feature of the fourth embodiment is that dead time factors are inserted between the virtual plant and the parameter identifier. Namely, as mentioned earlier, there exist various lags in a fuel metering control such as a lag generated by an air/fuel ratio sensor's detection, a lag generated by sensor outputs' A/D conversion timing, a lag caused by fuel injection amount calculation, a lag due to outputting timing thereof etc. and what is worse, the lags may change depending on the states of engine or fuel metering control system. Therefore, the fourth embodiment aims to conduct a timing adjustment between the plant and the parameter identifier using dead time such that it can cope with the change of the lags. For that purpose, the configuration illustrated in Figure 8 is slightly modified in the fourth embodiment as shown in Figure 35 wherein dead time factors are interposed between the virtual plant and the parameter identifier or the adaptive controller. To be more specific, explaining parameter identification laws in the configuration of Figure 35, the adaptive parameter  ∧ (k) can be expressed as Eq. 20 when using the method proposed by I.D. Landau et al. The identification error signal e* (k) and the gain matrix Γ (k) will be respectively expressed as Eq. 20 and Eq. 21.  (k) =  (k-1) + Γ(k-1)ξ(k-d)e*(k) e*(k)= D(z -1 )y(k) -  T (k-1)ζ(k-d) 1 + ζ T (k-d) Γ(k-1)ζ(k-d) Here, the orders of the  ∧ (k) vector and the gain matrix Γ are solely determined from the order of the virtual plant and the order of the dead time (delay time factor) of the virtual plant. Accordingly, when dead time varies in response to the engine operating states, the orders of the vector and matrix used in the parameter identifier must be varied. Namely, the algorithm itself should be modified. That is not practical when realizing the system actually. As an answer to the problem, the orders of the vector and matrix in the parameter identifier to be used for calculation is set to be possible maximum and dead time factors z -h , z -i and z -j are inserted as illustrated in Figure 35. As a result, if dead time actually becomes shorter than that, they can cope with various time lags existing between the input and output of the virtual plant. More specifically, at a high engine speed various time lags may become greater in total than a calculation cycle of the fuel metering control system so that the order of dead time could be d=4 at the maximum. The parameter identifier and adaptive controller should therefore be configured as d=4. On the other hand, the calculation cycle is relatively long at a low engine speed so that dead time becomes relatively short. If it is presumed that the order be d=2, the values h, i, j in Figure 35 will then be adjusted such that h=2, i=0 and j=2. Consequently, dead time of the virtual plant's output will be apparently d=4 if viewed from the parameter identifier and adaptive controller. Alternatively, the parameter identifier and adaptive controller can be configured in such a manner that dead time is set to be shorter than a possible maximum value. For example, assume that, when the order of dead time be d=4, the parameter identifier be configured to be prepared for a case in which the order of the plant's dead time is d=2. At such instance, if the dead time factors be configured as h=0, i=2 and j=0, the plant's output y(k) includes dead time order d=4 with respect to u(k-2), so that a difference therebetween will be 2. The identifier thus configured with its dead time order as quadratic can operate stably. With the arrangement, since no time error occurs between the desired cylinder fuel flow and the plant output indicative of the actual cylinder fuel flow even during transient engine operating state, an air/fuel ratio can be converted to a desired value. Here, it should be noted that, in order to cope with various time lags existing between the plant's input and output or their variations, dead time can be provided to the plant's input and output in an appropriate manner other than that mentioned above. Figure 36 is a view similar to Figure 22, but shows a fifth embodiment of the invention relating to the determination of the coefficient C used in estimating the actual cylinder fuel flow Gair. When using the testing apparatus illustrated in Figure 17, the throttle effective opening area increases with increasing throttle opening. In the actual engine such as shown in Figure 29, however, there exists a critical value at a certain level at which the effective area becomes maximum. In other words, when viewing the engine air intake system as a whole, resistance at the intake port or the air cleaner becomes greater so that the valve does not function as a throttle. Since an engine is a kind of pump, it has a fully operated area at which no more air will be inducted even if the throttle valve is opened more. At such a fully opened area, if the effective opening area obtained from the testing apparatus of Figure 17 is used in the calculation on an actual engine, a correct air flow will not be obtained. At the fully opened area, the critical value should therefore be used. Since the critical value should separately be determined from individual engine speeds as experienced in the case of the full throttle area, each throttle position corresponding to the fully opening area is obtained as the critical value for respective engine speeds and stored as a table data. A detected throttle opening is then compared with the critical value at the engine speed concerned and if the detected value is found to exceed the critical value, the detected value is replaced with the critical value, and the throttle effective opening area is calculated using the replaced critical value. Figure 34 illustrates this. Further, it should be noted that the cylinder air flow estimation was described in the first and fifth embodiments with reference to the fuel metering control using the adaptive control. This technique is applicable not only to the control disclosed herein but also to ordinary control of fuel metering or to ignition timing control.

Claims: A system for controlling fuel metering in multicylinder internal combustion engine (10), comprising: first means (50, S18) for determining a desired cylinder fuel flow Ti(k-n) in response to operating states (Ne, Pb, Pa, ΘTH) of the engine (10); second means (50, S26) for determining an actual cylinder air flow Gair(k-n) at a combustion cycle (k-n) at or earlier than the last; third means (50, S28) for dividing the value Gair(k-n) by air/fuel ratio A/F (k-n) in the cylinder concerned at a combustion cycle at or earlier than the last to determine an actual cylinder fuel flow Gfuel (k-n) for individual cylinders; and fourth means (50, S30) for determining a fuel injection amount Tout such that the actual cylinder fuel flow Gfuel (k-n) constantly coincides with the desired cylinder fuel flow Ti(k-n) for the individual cylinders of said engine (10), wherein: said fourth means comprises a wall adherence correction compensator for compensating the behavior of fuel adhering to an air intake passage (20) of said engine (10); and wherein said wall adherence correction compensator operates on the basis of predetermined characteristics; characterized in that said fourth means further comprises a dynamic compensator provided independently of and in series with said wall adherence correction compensator and said engine (10), assuming said wall adherence correction compensator and said engine (10) as a single plant such that said dynamic compensator dynamically compensates for the behavior of said engine (10) and of said wall adherence correction compensator, with parameters of the dynamic compensator being adjusted by an adaptive parameter adjuster depending on the actual cylinder fuel flow Gfuel in order to achieve said coincidence between the actual cylinder fuel flow Gfuel (k-n) with the desired cylinder fuel flow Ti(k-n) for the individual cylinders of said engine (10). The system according to claim 1, wherein said third means (50, S28) determines the actual cylinder fuel flow Gfuel(k-n) using the actual cylinder air flow Gair(k-n) at a combustion cycle (k-n) at or earlier than the last and the air/fuel ratio A/F (k-n) at the same combustion cycle. The system according to claim 1 or 2, wherein the air/fuel ratio is determined by the following steps of: A) approximating modeling the response delay in an air/fuel ratio sensor (46) provided at an exhaust system of the engine as a first-order delay to determine a state equation describing the behavior thereof; B) discretizing the state equation for a period delta T to determine a transfer function thereof; C) determining a transfer function inverse thereto and multiplying an output of the sensor (46) to estimate an air/fuel ratio of the mixture inducted in any cylinder of the engine (10); D) assuming that the estimated air/fuel ratio as an average indicative of the individual cylinders' combustion histories weighted by a predetermined value (C 1 , C 2 , C 3 , C 4 ) to establish a model describing the behavior of the exhaust system and determining a state equation using the individual cylinders' air/fuel ratios as variables defining its internal state; E) establishing an observer for observing the internal state to obtain its output, and F) estimating the air/fuel ratio of each cylinder from the observer output. The system according to any of the preceding claims 1 to 3, wherein the actual cylinder air flow Gair is determined from the steps of: A) assuming a throttle valve (16) provided at an air intake passage (12) of the engine (10) as an orifice to establish a fluid dynamic model and based on the model, determining air flow Gth passing therethrough at least using detected pressures at upstream and downstream of the throttle valve (16); B) determining air mass Gb filling a chamber (18) in the passage (12) extending from the throttle valve (16) to an intake port of the cylinder using the ideal-gas law; C) determining the change delta Gb of the air mass Gb in the chamber (18) from pressure change in the chamber (18); and D) estimating a cylinder air flow Gair by subtracting the change delta Gb of the air mass Gb in the chamber (18) from the throttle valve (16) passing air flow Gth. The system according to claim 4, wherein said air flow Gth passing through the throttle valve (16) in response to throttle opening is determined by using an equation based on a fluid dynamic model and defined as: Gth = C · S ·ρ · 2 g ·( P1-P2 ) γ where C: coefficient, S: throttle projection area, ρ: air density at throttle's upstream side, g: gravitational acceleration, γ: air specific weight on throttle's upstream side, P1: pressure on throttle's upstream side, P2: pressure on throttle's downstream side. The system according to claim 4 or 5, wherein the pressure (P1) upstream of the throttle valve (16) is measured at a position away from the throttle valve (16) at least by 1D when the air intake passage diameter is defined as D. The system according to any of preceding claims 4 to 6, wherein the pressure (P2) downstream of the throttle valve (16) is measured at a position away from the throttle valve (16) at least by 3D when the air intake passage diameter is defined as D. The system according to any of preceding claims 4 to 7, wherein the pressure (P2) downstream of the throttle valve (16) is determined from the pressure at the chamber (18). The system according to any of preceding claims 4 to 8, wherein resolving power of a sensor (36) for measuring throttle opening is set to be increased with decreasing throttle opening. The system according to any of preceding claims 4 to 9, wherein resolving power of a sensor (38) for measuring the pressure (P2) at downstream of the throttle valve (16) is increased with increasing pressure. The system according to any of preceding claims 5 to 10, wherein the coefficient (C) is determined from throttle opening and a value indicative of engine load at least one among manifold pressure (Pb), a deviation between manifold pressure (Pb) and atmospheric pressure (Pa) and a ratio of manifold pressure (Pb) to atmospheric pressure (Pa). The system according to claim 11, wherein the coefficient (C) is determined from throttle opening (Θth) and engine load in advance and is stored as mapped data in a computer memory. The system according to claim 12, wherein an interval between adjacent lattice points in the mapped data is set to be smaller with decreasing throttle opening (Θth). The system according to any of preceding claims 11 to 13, wherein a critical throttle opening at which engine load becomes maximum is determined with respect to engine speed (Ne) and when a detected throttle opening exceeds the critical throttle opening, the detected value is replaced with a critical value. The system according to one of claims 11 or 12, wherein the coefficient (C) includes at least the flow rate coefficient. The system according to any of preceding claims 5 to 14, wherein the pressures (P1) and (P2) are respectively represented by atmospheric pressure (Pa) and manifold pressure (Pb) and an answer in the equation's root square is calculated in advance to be stored as a mapped data in a computer memory. The system according to one of the preceding claims, wherein said fourth means further comprises an adaptive controller which uses a ratio between a fuel amount (Qt(k)) adhering to the air intake passage (20) of the engine (10) and that fuel amount (Qin(k)) flowing into a cylinder of the engine (10) as a state variable of a plant and identifies/adjusts at least one of the parameters (A, B) of a transfer function of the plant such that the actual cylinder fuel flow Gfuel coincides with the desired cylinder fuel flow and adjusts at least one of the parameters of the adaptive controller in response to the at least one identified parameter; wherein a dead time factor (Z -d ) corresponding to a delayed time between the timing at which control input (u(k)) is generated and the time at which the plant output (y'(k)) or an estimated value thereof is generated, is provided at least one among a first position between control input (u(k)) to the plant and the parameter identifier/adjuster, a second position between the plant output (y'(k)) or its estimated value and the parameter identifier/adjuster and a third position between the plant output (y'(k)) or its estimated value and the adaptive controller. The system according to claim 17, wherein the order of the dead time factor (Z -d ) is varied in response to at least one of the operating states of the engine and the fuel metering control system. The system according to claim 17 or 18, wherein the parameter identifier/adjuster uses at least one among a decreasing gain method, variable gain method and a constant trace method.

Inventors: AKAZAKI SHUSUKE, HASEGAWA YUSUKE, HIROTA TOSHIAKI, KOMORIYA ISAO, MAKI HIDETAKA, AKAZAKI, SHUSUKE, HASEGAWA, YUSUKE, HIROTA,TOSHIAKI, KOMORIYA, ISAO, MAKI, HIDETAKA
Applicants: HONDA MOTOR CO LTD, HONDA GIKEN KOGYO KABUSHIKI KAISHA
Application year: 1993
Main classification: F02D41
IPC codes: F02D41