Ignition timing control method for an internal combustion engine and apparatus therefor

Abstract: An engine control apparatus for shifting an ignition timing signal to an ignition coil (48) so that the engine produces optimum torque generates a pseudo random signal (26) which is preferably an M series signal and applies the pseudo random signal to a parameter of the engine having a noise frequency component such as the engine rotational speed y(t). The M series signal is a signal having a impulse type self correlation function. A mutual correlation function of the number of revolutions of the internal combustion engine and of the M series signal is computed (30-44) and an optimum ignition advance angle is supplied to an adder (46). The adder (46) combines the optimised signal ignition advance angle with a basic ignition timing signal (ϑ advB ) to provide an ignition timing signal (ϑig) to the ignition coil (48).

Description: Background of the Invention This invention relates to an ignition timing control method for an internal combustion engine and to an apparatus for performing said method. The invention has particular suitablility for controlling and maintaining an internal combustion engine ignition timing at an optimum position to obtain improved fuel consumption. Description of Related Art Under the same operating conditions such as amount of fuel, number of engine revolutions, load, fuel properties, etc., an internal combustion engine changes its generating torque when the ignition timing is adjusted, and generates the maximum torque at an optimum ignition timing. Accordingly, it is clear that the fuel consumption rate of the internal combustion engine will be improved if the ignition timing is continously varied so as to yield the maximum torque under different operating conditions. A present date technique is to hold ignition timing information in a read only memory (ROM) and it has been proposed to control the ignition timing of an internal combustion engine in accordance with the map data, the map data being prepared to indicate the ignition timing providing maximum output in response to the number of engine revolutions and load on the internal combustion engine. However, the optimum ignition timing fluctuates with behaviour of individual engines and due to ageing caused by carbon deposits, sensor drift, actuator drift, and in the use of fuels having different octane numbers. It has therefore been extremely difficult to control the engine in proper response to such fluctuating conditions. An attempt at overcoming such difficulties is disclosed in an article published in the SAE PAPER (SAE) 870083 (February 1982) pp. 43 - 50 which discloses a method for predicting ignition timing giving the maximum torque output from the detected rate of change of rotation of an internal combustion engine when the rotational speed is changed slightly by increasing or decreasing the ignition timing while the internal combustion engine is running. This method moved the ignition timing advance angle in proportion to the gradient of the output torque of the internal combustion engine, i.e. in accordance with the characteristic of output torque against ignition advance angle. Thus, denoting the output torque of an internal combustion engine by T, denoting the number of engine revolutions by N, and denoting the ignition advance angle by , then the following formula applies:- Δ T Δ = Δ T Δ N · Δ N Δ ≒ K Δ N Δ An optimum control is therefore achieved by applying the method termed the "hill climbing" method in the said article; that is to say instead of determining the change gradient of output torque to ignition advance angle (ΔT/Δ), a change gradient of the number of revolutions of the internal combustion engine to ignition advance angle (ΔN/Δ) is determined, and the amount of the ignition advance angle is moved in proportion to the gradient of the characteristic ΔN/Δ. To detect the optimum ignition timing with this method, it is necessary to determine the gradient of the number of engine revolutions to the ignition timing as described above and the optimum ignition timing is normally taken to be when the gradient change of the characteristic is zero. However, in the known prior art device, due to the low signal/noise ratio in the engine revolutions sensor, a sufficient determination cannot be obtained unless a considerable change is given to the number of revolutions of the internal combustion engine, which results in the vehicle jerking and reducing driving comfort. EP-A-0275169 discloses determining the slope of engine speed with respect to spark advance by perturbing the spark advance angle or, alternatively, determining the slope of engine speed with respect to a parameter relating to the fuel mixture by perturbing the fuel mixture. It is an object of this invention to provide a method and apparatus which detects an optimum ignition internal combustion engine timing without upsetting the normal operation of the internal combustion engine. Summary of the Invention According to one aspect of this invention there is provided an engine control apparatus as claimed in claim 1 herein. According to another aspect of this invention there is provided a method as claimed in claim 6 herein. Various optional features of the invention are defined in the dependent claims. It will, therefore, be understood that the present invention changes the ignition timing of an internal combustion engine in accordance with a search signal which may be a signal known as an "M series signal" having an impulse-type self correlation function and to control the ignition timing on the basis of the rate of change of the number of revolutions of the engine at that time. Brief Description of the Drawings The invention will now be described by way of example with reference to the accomanying drawings in which:- Figure 1(a) shows a block schematic diagram of a simplified form of an engine control apparatus in accordance with this invention, Figure 1(b) shows in simplified block schematic form apparatus for effecting corrective action, Figure 2 shows a block schematic diagram of one embodiment of the invention related to controlling ignition timing, Figure 3(a) shows an M series signal used in this invention, Figure 3(b) a self correlation function of the M series signal shown in Figure 3(a), Figure 3(c) shows a torque against ignition timing characteristic used in the embodiment of the invention, Figure 4 and Figure 5 show flow diagrams used in performing the invention, Figure 6(a) shows a register count for the M series signal, Figure 6(b) shows in further detail the M series signal, Figures 7 and 8 show further flow diagrams used in the operation of the invention, Figures 9(a)-9(h) show control timing diagrams for the programme used in the invention, Figures 10(a) and 10(b) are explanatory wave forms of the M series signal distribution used in the invention, Figure 11 shows a further embodiment in accordance with this invention used to control ignition timing of an internal combustion engine, Figure 12 shows a flow diagram of the operation of the embodiment shown in Figure 11, Figures 13(a)-13(e) show graphical wave forms of the results of a simulation using the invention, Figure 14 shows further aspects of the M series signal overlapping onto an ignition timing signal, Figure 15 shows the results of overlapping the M series signal on the acceleration of the vehicle, Figures 16(a)-16(c) shows the result of overlapping an M series signal onto an ignition timing signal with relation to the engine rotational speed and torque gradient, and Figures 17(a)-(c) show test results actually overlapping the M series signal on an ignition control showing the optimised ignition timing signal and change in engine rotational speed. In the figures like reference numerals denote like parts. Description of the Preferred Embodiments An apparatus of the invention in simplified form is shown in Figure 1a in which an internal combustion engine 1 has a parameter thereof which is to be sensed or corrected. The parameter to be sensed or corrected has a noise signal component thereon which makes detection of a small incremental change of the parameter impossible to detect in the prior art. In the present invention a psuedo random signal is generated which is herein termed a search signal and which may be an M series signal known per se , briefly described later herein but described more fully in "CORRELATION FUNCTION AND SPECTRUM - Measurement and Application thereof", edited by Takashi Isobe, published by Tokyo-Daigaku-Shuppansha. The search signal is generated by a generator 2 which outputs the search signal to be added by adder 7 into the parameter to be sensed or corrected, e.g. fuel injection time and/or ignition timing. The parameter with the search signal added thereto and, of course including the noise signal component, is applied as one input to a correlator 3, the other input of the correlator 3 being provided from the search signal generator 2. The correlator 3 thus is able to detect the combined engine parameter and the search signal and to correlate the search signal so that, effectively, the noise and the search signal are removed in the correlation process. The output of the correlator is applied to an integrator 4 which acts as a digital filter, the output of the intergrator 4 providing a correction factor γ 1 . The correction factor γ 1 as shown in Figure 1(b) is applied to a unit 5 for correcting or diagnosing sensor and/or engine parameters such as fuel injection time and/or ignition timing. The output of the unit 5 may also be applied to an alarm 6 shown in broken lines. An embodiment of the invention will now be described in more detail with reference to Figure 2. The speed, i.e. number of revolutions N, of an internal combustion engine is detected by a crank angle sensor 12, and air quantity Q a sucked into a cylinder of the engine is detected by an air flow sensor 14. The crank angle sensor 12 also supplies reference signals REF shown in Figure 9(a) and position indicating signals at, for example, 1° intervals as shown in Figure 9(b). The signals REF are generated for example at 110° before top dead center (TDC) and the REF signals and POS signals are supplied to a control unit (not shown in Figure 2). A divider 16 calculates the ratio of the air quantity Q a versus number of revolutions of the internal combustion engine N, and generates a signal corresponding to load L (L=Q a /N). An air fuel ratio compensator 18 receives signals representative of load L and engine revolutions N (and optionally the output of an air fuel ratio sensor O 2 ) and supplies a signal representative thereof to one input of a control unit 20, shown diagrammatically as a multiplier. The other input of the unit 20 is a standard injection timing signal T p derived by unit 19 that is proportional to load L. The unit 20 provides a fuel injection time signal T i for a cylinder of the internal combustion engine, the signal being provided to a fuel injector 22. The control unit 20 provides as an output the actual fuel injection time T i either by adding the injection time calculated by the air/fuel ratio compensator 18 to the standard fuel injection time T p (as shown in Figure 2) or by multiplying the standard time T p by a correction coefficient (not shown). The fuel injector 12 injects the fuel into the cylinder of the engine for the injection period T i . An ignition timing control unit 24 generates a basic ignition advance angle  advB which is determined in response to the number of revolutions of the internal combustion engine and load and read from a map held in ROM in the unit 14. A pseudo random signal having varying mark (Δ)/space ratio which repeats after time period nΔ where n is an integer is used as a correction signal and the pseudo random signal, search signal is an M series signal generated by a generator 26 located within a single-chip micro computer 50. The M series signal x and(t) is based on the data shown in Figfure 6(a) and 6(b) described later herein. The M series signal has a signal component of ignition advance angle Δ advM . A parameter y(t) representative of the number of engine revolutions is combined with the M series signal, the resultant is detected and the M series signal correlated by unit 28. The output of the correlator γ is applied to an integral control 30 and to a proportional control 32. The output of the proportional control 32, kγ is multiplied in unit 34 to provide a signal Δ advPI which is applied to an adder 36. The signal kγ is applied to a unit 38 evaluating 1-β and thence to unit 40 producing e -LS which is a lag operator where L is lag time and S is Laplace operator d/dt =∫ dt to provide an output Δ adv , P, L which is applied as another input to adder 36. The output of adder 36 is signal Δ advP which is applied to one input of another adder 42 having as another input the output Δ advI from unit 30. The resultant output from adder 42 is Δ adv O which is the optimised signal component of the ignition advance angle and given by x and 1 (t). The output from adder 42 is applied to another adder 44 which receives as another input the output Δ advM from the M series signal generator 26. The combined signals are applied to another adder 46 having input x and O (t) which is the basic ignition advance angle  advB to provide an output x(t) which is the ignition timing signal ig applied to the ignition coil 48. Thus a mutual correlation function is performed of the M series signal x and(t) and engine RPM N represented by y (t), the correlation function is integrated and the output γ is a phase integral which is itself integrated and after combining with the M series signal in adder 44 so the signal is overlapped onto the basic ignition advance angle  advB denoted by x and O (t) and the ignition timing supplied to the coil 48. The M series signal, as shown in Figure 3(a) has an amplitude a, a minumum pulse width Δ, period nΔ where n is the maximum sequence and in the present embodiment equals 31. However n equals 5 and 15 are integers found by trial and experiment for different engine types which may also be used. The M series signal has an impulse type self correlation function as defined by  x x (α) = n+1 n a 2 As will be described later herein the M series signal x and(t) is generated with the amplitude a in a range causing a change in the number of engine revolutions that will not be perceived by a driver of the vehicle and the M series signal is overlapped on the ignition timing signal ig. The mutual correlation function and the phase shift integral of this M series signal x and(t) versus the number of revolutions y (t) of the engine are calculated to determine output torque gradient γ. To set the ignition timing forward or backward from its basic position together with changing the magnitude thereof in response to the positive or negative sign and the magnitude of the output torque gradient γ, the output torque gradient is integrated and combined on the original engine timing signal. In such a manner the ignition timing is controlled so as to be maintained at the optimum position without the necessity as in the prior art of excessively fluctuating the engine speed of revolution so that jerking of the vehicle is not noted by a driver with use of this invention. Although discrete units 28-44 have been shown, in the preferred embodiment, the functions performed by these units are effected under software control in the micro computer 50. The theory behind the present invention will now be described. The impulse response g(α), when an M series signal x and(t) is used as the input signal in the method of the present invention is determined by calculating the mutual correlation function  x andy(α) of input x and(t) and the output y(t) based on the input signal x and(t). Accordingly if the following relation holds, x (t) = x o (t) + x 1 (t) the equations (1) and (2) below hold. Because x (t) changes more slowly than x and(t) it can be regarded as a DC component so that the signal x and(t) can be regarded as a signal varying about a mean level x (t). Similarly for the engine RPM y (t) may be considered as the DC component about which the engine RPM y and(t) fluctuates. x(t) = x (t) + x (t) y(t) = y (t) + y (t) The characteristics shown in Figure 3(c) of torque (which is equivalent to engine speed of rotation) against ignition timing (ig) has a characteristic which is initially positive and then negative. The optimum ignition angle  optimum is taken when the slope of the torque/ignition angle curve is zero. The slope of the curve is thus constantly monitored for the optimum angle at which torque is maximum. In this invention the slope of the characteristic is measured over the search signal amplitude and it will be seen from Figure 3(c) that if the search signal x and(t) has a sufficiently small amplitude then the torque/ignition angle characteristic is substantially linear. Therefore the relationship between the search signal x and(t) and the output component y and(t) which corresponds to the search signal (in other words the relationship between the ignition timing and the number of revolutions of the engine) is represented by the following equations (3) to (5) nΔ : One cycle of of M series signal Δ : Minimum pulse width of M series signal n : Sequence number of M series signal α : Phase shift between successive nΔ periods of x and(t). Further, the mutual correlation function x andy and(α) for the search signal x and(t) and the output signal y and(t) is represented by the following equation (6). From control theory output Here, x andx and(α) is a self correlation function for M-series signals, and is given by the following formula; Because the search signal x and(t), is an M series signal which includes all frequency components, its power spectrum density function x andx and(α) is constant, accordingly. Φ x x (α)=Φ x x (o) As a result, the self correlation function,  x andx and(t-α), which appears in equation (6), is represented by equation (8) using delta function δ;  x x (t-α)=Φ x x (o)·δ(t-α) Hence, the mutual correlation function  x andy and(α) shown in equation (6) is tranformed as follows; As is evident from the above, the impulse response g(α) is given by equation (10) below using the mutual correlation function x andy and(α) between x and(t) and y and(t). g(α)= x y (α) / Φ x x (o) where, x andx and(O) corresponds to the integrated value of self correlation function x andx and, and is given by the following equation; Φ x x (o) = (n+1)Δ a 2 /n = Z (constant) where a:amplitude of the M series signal The mutual correlation function  x andy and(α) is transformed as shown below using equation (2); Thus, g(α)={ x y(α) -  x y (α)}/Z where the second term of equation (13) x and y (α) is the mutual correlation function between the M-series signal x and(t) and the DC component of output y (t). The first term x andy(α) is a mutual correlation function between the M-series signal input x and(t) and output y(t). y(t) is composed of fluctuating components due to the influence of M-series signal x and(t), and the DC component from x(t) ; however, it is difficult to separate and detect these components so that a directly obtainable function is a mutual correlation function x andy shown by the following equation. However, the value of x and y (α) approximates to the value of x andy(α) if the value of α is large because the value of x and y (α) is no longer influenced by x and(t). Therefore, x and y (α) can be approximately substituted by the average value of x andy(α) in the interval between α1 and α2. where, α1 and α2 are bias compensation terms and they are selected to have values close to n . Δ. Further, output torque slope with respect to ignition angle γ, i.e. change of torque with respect to change of ignition angle ig, in the interval between α S - α L is given by equation (15), α S is the starting time of the integration in consideration of the leading adge of the impulse response due to the pseudo-white noise of the M-series signal. α L is the ending time of the integration inverval for impulse response integration. This is set in advance, in accordance with the impulse response characteristics. This output torque slope γ corresponds to the change in number of revolutions of the internal combustion engine, when the ignition timing is changed by a unit quantity by the search signal, and this is called the output torque gradient. In the embodiment of the present invention shown in Figure 2, the optimum ignition timing is more smoothly achieved by overlapping the further integration of the above-mentioned output torque gradient γ on the ignition timing signal ig. The invention will now be described by way of example in conjuction with various flow diagrams where the operation of micro computer 50 is described. Referring to Figure 4 the timing is started with a reference signal REF and a basic ignition advance angle  advB which is a function of engine revolutions N and load (Q a /N) is determined in step 401. In an optimised control routine 402 under the flag ON condition an M series ignition advance angle setting routine 403 is set to start, the M series routine being described more fully with reference to Figure 5. In an ignition advance angle routine 404 the ignition advance ig is determined using equation (16)  ig = advB + Δ advM + Δ adv O where, ig : ignition advance angle,  advB : basic ignition advance angle, Δ advM : M series signal component of the ignition advance angle, Δ adv O : optimised signal component of the ignition advance angle. It will be understood that the processing step 404 is shown schematically in Figure 2 as the summation of inputs of adders 44 and 46. In an ignition energising start timing routine 405, the power is supplied to the ignition coil 48. The M series signal is shown in Figure 6(b) having a abscissa of time against ignition angle ig and the mean amplitude of the M series signal is at  advB with the amplitude of the M series signal being represented by Δ advM . Referring now to Figure 5, the M series signal component ignition advance angle routine is shown in detail. In this routine the M series signals are generated by successive readout of bit data from previously set M series signal x and(t) data. The data is read out of a register, schematically shown in Figure 6(a) having register positions 2 O to 2 n-1 and having a total count herein shown as 2 MCNT . The first run of data is determined in step 510 and if affirmative a counter MCNT is cleared at step 511. Retrievals for the M series signal bit data are then performed at step 512, then, using equation: - where a is the amplitude of the M series signal. The M series signal component ignition advance angle Δ advM is thus generated at steps 513, 514 and 515. Next, at steps 516, 517 and 518, the counter/register MCNT is updated in accordance with following: where, n-1 : period of M-series signal. Figure 7 is a flow diagram of the slope of the torque/ignition timing slope detection and correlation performed by unit 28. The M series signal x and(t) is sampled together with the number of revolutions of the internal combustion engine N at step 701 and both sets of data are input and memorised in the micro computer. When data input corresponding to one cycle of M series signal is completed, that is to say when m. n samples are taken, where m is the number of samples in period Δ and n is the number of stages in the register shown in Figure 6(a) (typically m = 5 and n = 31) then the self correlation function of equation 13' is performed at step 703, it being remembered that N is equivalent to y(t). The output torque gradient γ is computed at step 704 in accordance with equations 14 and 15. In the above calculations m is an interger as will be described later herein. At the step 702 it is determined whether the data m and n is stored in memory and if it is then the calculation for γ can be performed. In the flow diagram of Figure 8 an optimised signal component of the ignition advance angle Δ adv O is determined. At step 801 Δ advI present is calculated from Δ advI present = Δ advI previous + G.γ where Δ advI present is the integral of the present value Δ adv , Δ advI previous is the integral of the previous value of Δ adv , G is the integration coefficient. The equation (19a) evaluated at step 801 is the output of unit 30 shown in Figure 2. At step 802 the following equation is calculated Δ advPI = k.β.γ where k is a proportional coefficient and β is an immediate factor of proportional output which in the present embodiment = 0.3. The terms used in this calculation are shown in Figure 14(a) and the calculation at step 802 is the output of unit 34 shown in Figure 2. The integration control gain coefficient G is a coefficient representing the relationship between an output torque gradient and an optmised ignition timing and is set in accordance with individual internal combustion engines. As shown in Figure 14 the optimised signal component ignition advance angle Δ adv O is preferably formed in two steps in conformity with the immediate factor of proportional output β and in accordance with a phase delay time L. Thus at step 803 it is determined whether time delay L has occurred and if it has then the progranme steps to step 804 and Δ advPL = k (1-β)γ is calculated. The output of step 804 is the output of the units 38 and 40 schematically shown in Figure 2 and the equation Δ adv O = Δ advI +Δ advP where Δ advP = Δ advPI + Δ advPL is calculated at step 805 and this is the optimised signal component of the ignition signal which is the output of adder 44 shown in Figure 2. By arranging that the optimised signal component of the ignition advance angle Δ adv O is given in two steps as shown in Figure 14, a sudden change in ignition timing is suppressed. Figure 9 shows the relative timings on which the respective calculation routines operate. As previously described, Figure 9(a) shows the ignition timing setting routine to start with the timing of reference signals REF which are generated for each cylinder. The reference signals are generated at 110° before top dead centre (TDC) of each cylinder. For a six cylinder engine, reference signals REF are generated every 120°, that is, three pulses are generated per revolution, i.e. two revolutions are performed in one cycle so that six reference signals REF are generated during one cycle. Reference signals R 1 to R 3 correspond to the first cylinder to the third cylinder only and the period T ref of the reference signal REF becomes smaller as the number of engine revolutions increases. The position signals are shown in Figure 9(b) and are generated at 1° intervals at rotation of the crank shaft. The ignition coil current is controlled in response to the calculated results of this routine, and an ignition pulse is generated at the previously determined ignition timing. The flow time for the ignition coil current is determined upon conditions such as output voltage of the battery and number of revolutions of the internal combustion engine, and the starting time for the current T S is adjusted to a value computed by the ignition advance angle setting routine, as shown in Figure 9(d) and 9(e). For example, if an M series signal such as the one shown in Figure 9(c) is produced, and the ignition advance angle is changed with ±a then the starting time for the current T S is modified with ±a resulting in an adjustment T f of the ignition coil current and the ignition pulse shown in Figure 9(f) is generated at time T f . Independently from the ignition timing setting routine which is set to start synchronously with reference signal REF, an optimised control routine starts at an optimised control timing which is determined by dividing the reference signal REF into 1/m, where m is predetermined integer. Figures 9(g) and 9(h) show an embodiment where m = 5. As the timing period T ref /m at which the optimised control routine is set to start is proportional to the reference signal REF, the number of revolutions of the internal combustion engine is detected by measuring the interval of the optmised control timing operation. Since the detected number of revolutions has the same value within the period from one optimised control timing pulse generation to the next timing pulse generation (such as an interval T shown in Figure 9(g)), the optimised control routine is set to start at anywhere within the interval T. Any number from 1 to 5 can be selected as the value for integer m, however, even if a larger number of m is selected, the detected number of revolutions is virtually the same at low speed running and such a larger number will only result in increasing a burden on the micro-computer. In practice, values such as 1 or 2 are adequate. If the ignition advance angle setting routine and the optimised control routine are independently controlled as described above, both routines are not always sychronised and, moreover, priority may be given with regard to either of the processings. As a result, the optimised control routine may be run on a time basis; further if there is insufficient processing time, the processing of the ignition advance angle setting routine may be given priority so that the control can be made certain. Additionally, as shown in Figure 15, the processing may be separately executed during the measuring period for obtaining an output torque gradient in every period of the M series signal T ref .N and during the control output period so as to control the ignition timing at an optimised value. The minimum pulse widths Δ of the M series signal is an integer as large as the number of combustion strokes of the internal combustion engine. In the case of a six cylinder engine, for example, a reference signal REF is generated at every 120°, that is to say, six signals for every 2 revolutions, and the minimum pulse width Δ is set at an integer as large as the period T ref of the reference signal REF. For example, with an M series signal, if the minimum pulse width Δ is set at the same magnitude as the number of combustion strokes, then the result is as shown in Figure 10(a), and if the minimum pulse width is set to be six times as large as the number of combustion strokes then the result is as shown in Figure 10(b). If the minimum pulse width is set at the number of combustion strokes, all the cylinders are given the same ignition timing signal. If the minimum pulse width his set as a magnitude less than the number of combustion strokes it may happen that two or more ignition timing commands are given simultaneously to one cylinder or the M series signal falls into disorder. This minimum pulse width is set at a small magnitude with an increasing number of engine revolutions. Another embodiment of the present invention will now be described, which embodiment facilitates a decreased memory capacity with respect to the embodiment shown in Figure 2. In the calculations for the indicial (step) response y, the equation is transformed into the form of equation (20) below by exchanging the time integral in the mutual correlation function with the integral of the above phase α : - where Z : power of auto correlation of M series signal, i.e. area of triangle having base 2Δ of Figure 3(b), and X(t) is a function corresponding to the intergration by parts of the signal x and(t) represented by equation (21) below, and depends on x and(t) only, with no relation to the response signal y(t) of an internal combustion engine control system. From equation (12): Reforming the above, the indicial response γ is represented by: k 1 = 2 N+1 ·Δ·a k 2 = α L -α S α 2 -α 1 X(t), which is given by equation (24), is the function which corresponds to the partially integrated value of the search signal x and(t), and which is called the correlation signal. Not all the data of this correlation signal X(t) needs to be stored in a memory, provided the initial value X(0) is first determined and the difference is calculated at each timing. Now, when a sampling period is denoted by T S , the following equations are used for the determination. X(t) - χ(t-Ts) = Ts [ x (Ts+ Δ) - x (t+Δ - (p+1)Ts) -K 2 {x(t-α 1 )-x(t-α 1 -(q+1)Ts}] provided = αL-α S Ts , q = α 2 -α 1 Ts If the time integral in equation (28) is approximated by a moving average, i.e. a time shifted repetitive calculation of integration of a portion of X(t), the necessary memory capacity for the integral calculation will be greatly reduced. Referring to Figure 11, the micro processor 150 in this embodiment has an M series signal generator 26 supplying signal x and(t) to an adder 44 added and to a correlation signal generator 128. The output from the correlation signal generator is signal X(t) which is supplied as one input to a multiplier 129, the other input of the multiplying being signal y(t) representative of the engine rotational speed. Output from multiplier 129 is fed to a time integrater 130 the output of which is signal γ supplied to an integral control unit 131. Output from unit 131 is signal Δ adv O which is supplied to another input to the adder 44. The remaining part of the circuit shown in Figure 11 is similar to that shown in Figure 2. The correlation function X(t) is calculated in accordance with equation (28) sychronously with the M series signal and is stored in a memory. The optimised control programme for execution by a micro computer is shown in the flow diagram of Figure 12. In a data input 1201 the number of revolutions of the internal combustion engine y is sampled and the correlation signal is generated at step 1202 to provide X(t) as shown in equation (21). At step 1203 the following equation is calculated γ present = γ previous + Xy When the above processing is executed for one cycle period of M series signals determined by step 1204, the optimimum signal component ignition advance angle Δ adv O is determined using equation (19) (c) at step 1205. The output torque gradient γ is reset at step 1206 to prepare for the next cycle calculation. In this embodiment the correlation function is calculated in a step-by-step method. Therefore there is no need to store the M series data signals x(t) and the number of revolutions of the internal combustion engine y over one cycle of the M series signal, thus the memory capacity required is greatly reduced. Because, the integral by phase α is executed in advance, real time execution is limited to the time integral only and the time required for calculations is therefore greatly reduced. Figure 13 shows the results of some simulations of the present invention when applied to a six cylinder internal combustion engine. Figure 13(a) shows M series signal changes of plus or minus 1° which are overlapped onto an ignition timing signal and the number of engine revolutions changed with the M series signal by approximately 50 RPM as shown by Figure 13(b). A mutual correlation function between the detected number of revolutions of the engine was calculated for each period of the M series signal to provide an output torque gradient as shown in Figure 13(c) and as a result of the overlap of the integrated value of the output torque gradient on the ignition timing signal, so the ignition timing moved from its initial position of 20° before TDC to a new position of 28° before TDC (the optimum position) in about 4 seconds as shown by Figure 13(d). At this movement, the acceleration of the vehicle in the direction of travel was within ±0.03 G, which is in a range that would not be perceived by a driver as shown in Figure 13(e). Referring again to Figure 15, the solid line shows the crank shaft rotation where two steps are taken to reach Δ adv O whereas if a single step is taken then the larger change in engine rotation causes the crank oscillate as shown by the broken line in Figure 15. In Figure 14 there is shown an example when the M series signal is successively overlapped on the ignition timing signal. If the M series signal is given a change of ±2° as shown in Figure 16(a) then the number of revolutions of the crank shaft changes by approximately ±30 RPM as shown in Figure 16(b). When the M series signal was overlapped for approximately 600 msec, the torque gradient γ changed by about 6.5 RPM/° as shown in Figure 16(c). Considering such a change in relation to the embodiment of Figure 2, the torque gradient is determined in such a way that the mutual correlation function between the M series signal x and(t) and the output y(t) is calculated using equation (13'), then using this correlation function, the torque gradient was determined with equations (14) and (15). Figure 17 shows the test results using an actual car, where the M series signal overlapped for 620 msec, as shown in Figure 17(a). As a result, the ignition timing was corrected by about 10° from 23° to 30° as shown in Figure 17(a) and 17(b) and the delay time L was 140 msec. As a result, the revolution speed exhibited a hill climbing characteristic as shown in Figure 17(c) and the ignition timing moved to the optimum position. Having described the invention, it will be understood that the invention renders ignition timing control of an internal combustion engine smoother, providing very small velocity changes for change in ignition timing. Although the subject invention has been described in relation to modifying the ignition timing, the invention is not intended to be so restricted. In this respect the M series signal could be used to incrementally modify fuel injection pulse width or could be used to detect an output from an O 2 sensor to control fuel injection to thereby provide a cleaner exhaust gas.

Claims: An engine system control apparatus for controlling an engine system in response to a control signal, including means (2) for generating a pseudo random signal which fluctuates an engine system control parameter in order to examine a behaviour and responsiveness characteristic of said engine system, said pseudo random signal being applied to said control signal which is produced according to said engine system control parameter, means (3) for detecting an output signal corresponding to said pseudo random signal which reflects an operational condition depending on said responsiveness characteristic of said engine system, means (3) for deriving a mutual correlation function from a relationship existing between said pseudo random signal and said output signal, wherein said mutual correlation function is obtained through processing including a phase shift of at least one of said pseudo random signal and said output signal, forming plural multiplications therebetween and summing said plural multiplications, whereby a correct control input to the engine system is obtained, and means (5) for correcting said control signal in accordance with said correct control input. An apparatus as claimed in claim 1 wherein said pseudo random signal has an impulse type self correlation function. An apparatus as claimed in claim 1 wherein said correction means (5) for said control signal uses a plurality of the results of said multiplication, said control signal including an ignition timing signal, and a fuel injection signal. An apparatus as claimed in claim 1 wherein said output signal from the engine reflecting the operational condition thereof is indicative of the number of revolutions of the engine. An apparatus as claimed in any preceding claim wherein there is provided an air quantity sensor (14) for measuring air quantity Q 3 supplied to the internal combustion engine, an injector (22) for supplying fuel to said engine, micro-computer means for controlling said injector and said ignition coil (48), said micro-computer being arranged to generate a fuel injection period signal T i in dependence upon a load on the engine (L=Q a /N) determined by the ratio of outputs of the air quantity sensor and the engine revolution sensor. A method of deriving a mutual correlation function from an input signal to an engine system and an output signal therefrom in order to examine a behaviour and responsiveness characteristic of said engine system, characterised by comprising the steps of: applying a pseudo random signal to said input signal, said pseudo random signal fluctuating said input signal within a permissible level which is marginal; detecting a component associated with said pseudo random signal in said output signal which is subjected to actual behaviour of said engine system; and calculating a mutual correlation function by: shifting a phase of at least one of said pseudo random signal and said output signal; forming plural multiplications of said pseudo random signal and said output signal upon shifting the phase of at least one of them; and adding up said plural multiplications to formulate the mutual correlation function. A method as claimed in claim 6 wherein the pseudo random signal has an impulse type self correlation function. A method as claimed in claim 6 wherein the control signal is an ignition timing control for an internal combustion engine and the result of said correction is to provide an optimum ignition timing signal. A method as claimed in claim 6 wherein said pseudo random signal has an impulse type self correlation function and is overlapped onto an ignition timing signal, determining an output torque gradient relative to ignition timing caused by the pseudo random signal from said mutual correlation and the number of revolutions of the internal combustion engine detected by a revolution sensor, generating a correction signal obtained by integration of said output torque gradient and correcting said ignition timing by the correction signal to produce an up/down manipulation of said ignition timing signal to thereby vary the number of revolutions of the engine. A method as claimed in any of claims 6 to 9 wherein said pseudo random signal is an M series signal having binary values. A method as claimed in any of claims 6 to 10 wherein said pseudo random signal has a minimum pulse width defined by the number of combustion strokes of the internal combustion engine. A method as claimed in any of claims 6 to 11 wherein said pseudo random signal has a predetermined period which is arranged to become shorter with increasing number of revolutions of the engine. A method as claimed in claim 6 wherein the number of engine revolutions N are detected, an air quantity Q a to be supplied to the engine is detected, a micro-computer generates a fuel injection period signal T i in dependence upon the load on the engine (L=Q a /N) determined by the ratio of air quantity to engine revolutions and said micro-computer also generates a basic ignition timing signal in dependence upon load on the engine and number of engine revolutions ( advB ) the basic ignition timing signal being combined with the result of integration for application to an ignition coil.

Inventors: ISHII JUNICHI, KANEYASU MASAYOSHI, KURIHARA NOBUO, ISHII, JUNICHI, KANEYASU, MASAYOSHI, KURIHARA, NOBUO
Applicants: HITACHI LTD, HITACHI, LTD.
Application year: 1990
Main classification: F02P5
IPC codes: F02P5, F02D37, F02D41, F02D45