METHOD AND SYSTEM FOR OBTAINING REFERENCE SIGNALS FOR VEHICLES CONTROL SYSTEMS AND CORRESPONDING CONTROL SYSTEM

Abstract

A method for obtaining reference signals for vehicle control systems, in function of a vehicle geographical position along a travel route, includes providing data relating to the vehicle and data relating to a route to travel, and determining a vehicle driving force reference signal and a vehicle speed reference signal through a first optimisation process configured to optimise the driving force along the travel route. An engaged gear reference signal, in function of the positions of the vehicle along the travel route, is determined through a second optimisation process configured to optimise fuel consumption of the vehicle along the travel route. The second optimisation process is subsequent to the first optimisation process, and the data relating to the travel route, as well as the driving force reference signal and the speed reference signal, is received as input, determined through the first optimisation process.

Claims

1. Method for obtaining reference signals for systems for controlling a vehicle, in function of a geographical position of said vehicle along a route to travel, comprising the operating steps of: providing at least data relating to said vehicle and at least data relating to a route to travel; determining, on the basis of said data, at least one reference signal of a driving force F of said vehicle and at least one reference signal of a speed for said vehicle, in function of the position s of said vehicle along said route to travel, through a first optimisation process configured to optimise the driving force F along said route to travel; and determining at least one reference signal of an engaged gear and/or of gear change u.sub.sh of said vehicle, in function of the position s of said vehicle along said route to travel, through a second optimisation process configured to optimise a fuel consumption of said vehicle along said route to travel; said second optimisation process being subsequent to said first optimisation process and receiving as input said data relating to said vehicle, said data relating to a route to travel, as well as said at least one reference signal of the driving force F and said at least one reference signal of the speed , determined through said first optimisation process.

2. Method according to claim 1, wherein said data of said vehicle comprise information related to a mass, an aerodynamic resistance, a rolling friction, a curve of a torque developed by an engine, one or more fuel flow maps and transmission ratios of said engine, wherein optionally said data are constant or variable during the life of said vehicle and/or of said engine of said vehicle and wherein said variable data can be more optionally estimated in-line or not in-line.

3. Method according to claim 1, wherein said data of said route comprise: an itinerary of said route to travel including a starting geographical position and a destination geographical position, a current geographical position, a slope and/or an elevation along said route to travel, and optionally one or more data selected from the group comprising data relating to speed limits established along said route to travel, an estimated arrival time or a desired maximum travel time, at least one desired average speed, at least one desired maximum speed and at least one desired minimum speed, at least one degree of adherence to said desired average speed, to said at least one desired maximum speed and to said at least one desired minimum speed, data relating to traffic conditions along said route to travel and data relating to weather conditions along said route to travel.

4. Method according to claim 1, wherein said first optimisation process is an analytical process of minimisation of a cost function J(F) of the driving force F required to drive said vehicle having mass m along said route to travel.

5. Method according to claim 4, wherein said minimisation of the cost function J(F) of the driving force F is given by the relationship: $\begin{array}{cc}\underset{F}{\min }.Math.\left\{J\left(F\right)\right\}=\underset{F}{\min }.Math.\left\{{}_0^S.Math.F^2.Math.\mathrm{ds}\right\}& \left(1\right)\end{array}$ with the following state equation: $\begin{array}{cc}m.Math.\frac{\mathrm{dv}}{\mathrm{dt}}=F-F_{\mathrm{grade}}-F_{\mathrm{drag}}-F_{\mathrm{roll}}& \left(2\right)\end{array}$ and the following constraints: $\begin{array}{cc}{v}_{\min }v{v}_{\max }& \left(3\right)\\ F{F}_{\max }& \left(4\right)\\ {}_0^S.Math.\frac{1}{v}.Math.\mathrm{ds}{T}_{\max }=\frac{S}{{\stackrel{\_}{v}}_{\mathrm{des}}}& \left(5\right)\end{array}$ where is the speed of said vehicle, s is the length travelled by said vehicle along said route of overall length S, .sub.min is the desired minimum speed and .sub.max is the desired maximum speed, F.sub.max is the maximum available driving force of said vehicle, T.sub.max is a maximum travel time and .sub.des is the desired average speed, and taking account of exogenous forces as follows:
F.sub.grade=m g sin() (6)
F.sub.drag=.sub.airC.sub.xA.sup.2 (7)
F.sub.roll=m g cos()[C.sub.r0+C.sub.r1.sup.2](8) where m is the mass of said vehicle (V), g is the gravitational acceleration, is an angle corresponding to the slope along said route to travel, .sub.air is an ambient air density along said route to travel, C.sub.x is an aerodynamic resistance coefficient, A is a front area of said vehicle, C.sub.r0 is a first rolling resistance coefficient, and C.sub.r1 is a second rolling resistance coefficient.

6. Method according to claim 5, wherein said analytical process of minimisation of the cost function J(F) of the driving force F performs a minimisation of the integral of the Lagrangian (F), function of the driving force F, given by: $\begin{array}{cc}\left(F\right)={}_0^S.Math.\left\{F^2+a\left(s\right)\left[\frac{\mathrm{dK}}{\mathrm{ds}}-F+m.Math..Math.g.Math..Math.\sin \left(\right)+\frac{{}_{\mathrm{air}}.Math..Math.C_x.Math..Math.A}{m}.Math.K+m.Math..Math.g.Math..Math.\cos \left(\right).Math.\left(C_{r.Math..Math.0}-C_{r.Math..Math.1}.Math.\frac{2.Math..Math.K}{m}\right)\right]+g_1\left(K\right)+g_2\left(K\right)+g_3\left(F,v\left(K\right)\right)+\left(\frac{\mathrm{dt}}{\mathrm{ds}}-\sqrt{\frac{2.Math..Math.K}{m}}\right)\right\}.Math.\mathrm{ds}+b_1\left(K\left(0\right)-K_0\right)+b_2\left(t\left(0\right)-0\right)+{\left(t\left(S\right)-T_{\max }\right)}^2& \left(15\right)\end{array}$ where the functions g.sub.i, g.sub.2 and g.sub.3 are defined by the following equations: $\begin{array}{cc}{g}_1\left(K\right)=\{\begin{array}{cc}{c_1\left(K-K_{\max }\right)}^2& \mathrm{if}.Math..Math.K>K_{\max }\\ 0.Math.& \mathrm{if}.Math..Math.K{K}_{\max }\end{array}& \left(16\right)\\ g_2\left(K\right)=\{\begin{array}{cc}{c_2\left(K_{\min }-K\right)}^2& \mathrm{if}.Math..Math.K<K_{\min }\\ 0.Math.& \mathrm{if}.Math..Math.K{K}_{\min }\end{array}& \left(17\right)\\ g_3\left(F,v\right)=\{\begin{array}{cc}{c_3\left(F-F_{\max }\left(v\right)\right)}^2& \mathrm{if}.Math..Math.F>F_{\max }\\ 0.Math.& \mathrm{if}.Math..Math.F{F}_{\max }\end{array}& \left(18\right)\end{array}$ where: K is a kinetic energy associated to the driving force F, K.sub.min is a kinetic energy associated to the desired minimum speed .sub.min, K.sub.max is a kinetic energy associated to the desired maximum speed .sub.max, c.sub.1, c.sub.2, c.sub.3 are penalisation coefficients, t(s) is the travel time of said vehicle (V), is an amplitude of a quadratic penalty term, (s) is the adjoint field, wherein the following constraints have to be satisfied $\begin{array}{cc}\frac{\mathrm{dK}}{\mathrm{ds}}-F+m.Math..Math.g.Math..Math.\sin \left(\right)+\frac{{}_{\mathrm{air}}.Math.C_x.Math.A}{m}.Math.K-m.Math..Math.g.Math..Math.\cos \left(\right)\left[C_{\mathrm{r0}}-C_{r.Math..Math.1}.Math.\frac{2.Math.K}{m}\right]& \left(9\right)\\ \frac{1}{2}.Math.{\mathrm{mv}}_{\min }^2=K_{\min }K{K}_{\max }=\frac{1}{2}.Math.{\mathrm{mv}}_{\max }^2& \left(10\right)\\ F{F}_{\max }& \left(11\right)\\ \frac{\mathrm{dt}}{\mathrm{ds}}=\sqrt{\frac{m}{2.Math.K}}& \left(13\right)\\ t\left(S\right)=T_{\max }& \left(14\right)\end{array}$ wherein said minimisation of the integral of the Lagrangian is executed by an algorithm of computation of the gradient of the integral of the Lagrangian $\frac{}{.Math..Math.F},$ wherein said maximum available driving force F.sub.max is optionally expressible as an analytical function approximating an envelope of the maximum driving force for each gear.

7. Method according to claim 6, wherein said algorithm of computation of the gradient of the integral of the Lagrangian $\frac{}{.Math..Math.F}$ is an adjoint-based iterative algorithm comprising the following steps: computing a gradient $\frac{}{.Math..Math.F}$ of the Lagrangian (F); iteratively using said gradient $\frac{}{.Math..Math.F}$ of the integral of the Lagrangian (F) to determine a search direction .sub.k along which searching a minimum of said gradient $\frac{}{.Math..Math.F}$ of the integral of the Lagrangian (F); executing a line minimisation; and repeating said iterations until at least one criterion to stop said iterations, indicating that a minimum of the integral of the Lagrangian (F) has been reached, is satisfied wherein the penalisation coefficients c.sub.1, c.sub.2, c.sub.3 are optionally constant or optionally progressively increased in correspondence of successive iterations, and wherein the amplitude of the quadratic penalty term is optionally constant or optionally dynamically modifiable in correspondence of successive iterations.

8. Method according to claim 7, wherein said gradient $\frac{}{.Math..Math.F}$ of the integral of the Lagrangian (F) is given by: $\begin{array}{cc}\frac{}{.Math..Math.F}=-a+2.Math.F+\frac{d.Math..Math.g_3}{\mathrm{dF}}& \left(22\right)\end{array}$ wherein said search direction .sub.k in correspondence of the k-th iteration of said adjoint-based iterative optimisation algorithm is given by: $\begin{array}{cc}\{\begin{array}{cc}{p}_k=-\frac{}{.Math..Math.F}.Math.{|}_k& \mathrm{if}.Math..Math.k=1\\ p_k=-\frac{}{.Math..Math.F}.Math.{|}_k.Math.+{}_k.Math.p_{k-1}& \mathrm{if}.Math..Math.k>1\end{array}& \left(23\right)\end{array}$ where represents the momentum of the search direction in correspondence of the preceding iteration, that is optionally determined through the Polak-Ribire formula: $\begin{array}{cc}{}_k=\max \left(0;\frac{{\left(-\frac{}{.Math..Math.F}.Math.{|}_k\right)}^T.Math.\left(-\frac{}{.Math..Math.F}.Math.{|}_k.Math.+\frac{}{.Math..Math.F}.Math.{|}_{k-1}\right)}{{\left(-\frac{}{.Math..Math.F}.Math.{|}_{k-1}\right)}^T.Math.\left(-\frac{}{.Math..Math.F}.Math.{|}_{k-1}\right)}\right)& \left(24\right)\end{array}$ wherein said line minimisation computes a step length h.sub.k of the driving force F ensuring a maximum decrease of a cost function given by:
F.sub.k+1=F.sub.k+h.sub.k.sub.k (25) iteratively solving a nonlinear equation for the scalar h given by: $\begin{array}{cc}\frac{\left(F+\mathrm{hp}\right)}{h}=\frac{\left(F+\mathrm{hp}\right)}{\left(F+\mathrm{hp}\right)}.Math.\frac{\left(F+\mathrm{hp}\right)}{h}=\left(2.Math.\left(F+\mathrm{hp}\right)-a+\frac{d.Math..Math.g_3}{\mathrm{dF}}\right)=0& \left(26\right)\end{array}$

9. Method according to claim 7, wherein said at least one criterion to stop said iterations is given by: $\begin{array}{cc}\{\begin{array}{c}\frac{J\left(F_k\right)-J\left(F_{k-1}\right)}{J\left(F_{k-1}\right)}{\mathrm{.Math.}}_J\\ \frac{t\left(S\right).Math.{\text{|}}_k-t\left(S\right).Math.{\text{|}}_{k-1}}{t\left(S\right).Math.{\text{|}}_{k-1}}{\mathrm{.Math.}}_t\end{array}& \left(27\right)\end{array}$ where .sub.j and .sub.t are two predetermined threshold values.

10. Method according to claim 1, wherein said second optimisation process is a dynamic process of minimisation of a cost function J.sub.2 of the fuel consumption of said vehicle along said route to travel expressed as sum of M terms, where M is the number of segments having length s.sub.2 forming said route to travel having length S, of the fuel consumption by said vehicle, along a i-th segment of length s.sub.2 of said route of length S, and of penalising factors aimed at penalising gear changes on uphill and downhill sections.

11. Method according to claim 10, wherein said minimisation of the cost function J.sub.2 of the fuel consumption is given by the following relationship: $\begin{array}{cc}\underset{{}_i{}_i}{\min }.Math.\left\{J_2\right\}=\underset{{}_i{}_i}{\min }.Math.\left\{\underset{i=1}{\overset{M}{.Math.}}.Math..Math.\left(\frac{W_{f,i}\left({}_i\right).Math..Math..Math.s_2}{v_i}+{}_{1,i}\left({}_i,m\right).Math..Math.u_{\mathrm{sh},i}.Math.+{}_2.Math..Math.u_{\mathrm{sh},i}.Math.\right)\right\}& \left(28\right)\end{array}$ with the following state equation:
.sub.i+1=.sub.i+u.sub.sh,i (29) and the following constraints:
.sub.i.sub.i (30)
.sub.i:n.sub.e,minn.sub.e(.sub.i, K.sub.i)n.sub.e,max (31)
u.sub.sh,i[2; 1; 0; +1; +2](32) where: $\begin{array}{cc}M=\frac{S}{.Math..Math.s_2}& \left(33\right)\\ {}_{1,i}=\{\begin{array}{cc}m.Math.{}_i& \mathrm{if}.Math..Math.{}_i>0\\ 0& \mathrm{if}.Math..Math.{}_{i}0\end{array}& \left(34\right)\end{array}$ and wherein: u.sub.sh,i is a gear change signal in segment i; M is a number of summing terms of the cost function J.sub.2 of the fuel consumption; W.sub.f,i is a fuel flow rate in segment i; s.sub.2 is a length of each one of the M segments; .sub.1,i is a first penalisation coefficient that is variable in function of the mass m of the vehicle and of the slope .sub.i in segment i; .sub.2 is a second penalisation coefficient that is constant; .sub.i is the gear of the vehicle in segment i, included in a set of gears .sub.i, n.sub.e is the engine speed, ranging between a maximum value and a minimum value n.sub.e,max and n.sub.e,min.

12. Method according to claim 11, wherein fuel flow rate W.sub.f, function of the engine speed n.sub.e and of the torque T.sub.e developed by the engine, is expressed in function of said engaged gear as follows: $\begin{array}{cc}{T}_e\left(F,\right)=\frac{{\mathrm{FR}}_w}{i\left(\right)}& \left(35\right)\\ n_e\left(v,\right)=\frac{\mathrm{vi}\left(\right)}{R_w}& \left(36\right)\end{array}$ where i() is a global transmission ratio, with the engaged gear , and R.sub.w is a wheel radius of said vehicle, and F and are said at least one reference signal of the driving force and said at least one reference signal of a speed determined through said first optimisation process.

13. System configured for obtaining reference signals for systems for controlling a vehicle, in function of a geographical position of said vehicle along a route to travel, comprising: at least one data collection device, configured to output at least data relating to said vehicle and at least data relating to said route to travel, and one or more processing units, configured to receive from said at least one data collection device said data relating to said vehicle and said data relating to said route to travel, said one or more processing units being configured to execute the method for obtaining reference signals for systems for controlling a vehicle according to claim 1.

14. System configured for controlling the speed of a vehicle, comprising at least one system configured for obtaining reference signals according to claim 13 and one or more controlling apparatuses configured to control a speed of said vehicle and an engaged gear of said vehicle along said route to travel on the basis of reference signals received from said at least one system configured for obtaining reference signals.

15. System according to claim 14, wherein said one or more controlling apparatuses comprise one or more apparatuses selected from the group comprising a cruise control system configured to control an engine unit, a unit controlling a gearbox of said vehicle and an apparatus configured to act on at least one braking device of said vehicle.

16. (canceled)

17. Set of one or more computer readable media, having stored thereon one set of one or more computer executable instructions that, when executed by one or more processing units, cause said one or more processing units to perform a method for obtaining reference signals for systems for controlling a vehicle, based on a geographical position of said vehicle along a route to travel, comprising the operating steps of: providing at least data relating to said vehicle and at least data relating to a route to travel: determining, on the basis of said data, at least one reference signal of a driving force F of said vehicle and at least one reference signal of a speed for said vehicle, in function of the position s of said vehicle along said route to travel, through a first optimisation process configured to optimise the driving force F along said route to travel; and determining at least one reference signal of an engaged gear and/or of gear change u.sub.sh of said vehicle, in function of the position s of said vehicle along said route to travel, through a second optimisation process configured to optimise a fuel consumption of said vehicle along said route to travel; said second optimisation process being subsequent to said first optimisation process and receiving as input said data relating to said vehicle, said data relating to a route to travel, as well as said at least one reference signal of the driving force F and said at least one reference signal of the speed , determined through said first optimisation process.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The present invention will be now described, for illustrative but not limiting purposes, according to its preferred embodiments, with particular reference to the drawings in the accompanying Figures, wherein:

[0023] FIG. 1 shows a first flow diagram of the method according to the present invention;

[0024] FIG. 2 is an example block diagram of a system configured to perform the method of FIG. 1;

[0025] FIG. 3 shows an example block diagram of first optimisation process of the method of FIG. 1;

[0026] FIG. 4 is a representation of an approximation of the driving force of a vehicle, according to the method of the present invention;

[0027] FIG. 5 shows an example block diagram of an iterative optimization procedure, carried out during the first optimisation process of the method according to the invention;

[0028] FIG. 6 shows an example block diagram of second optimisation process of the method of FIG. 1;

[0029] FIG. 7 is an exemplary representation of one between various versions of a system configured to implement the method according to the present invention;

[0030] FIG. 8, shows the trend of the speed profiles (b), fuel consumption (c) and gear switching (d), respectively, obtainable along a route shown in (a) with a traditional method (CC) for adjusting the speed of a vehicle and with the method of FIG. 1 (OPT); and

[0031] FIG. 9 shows some experimental results supporting the efficacy of the present invention.

EMBODIMENTS OF THE INVENTION

[0032] With reference to the accompanying Figures, in particular at FIGS. 1 and 2 it will be noted that the method for obtaining reference signals for control systems of a vehicle according to the present invention is generally indicated with reference 1 and comprises at least the following operating steps of: [0033] providing (100) at least data relating to said vehicle V and at least data relating to a route to travel (step 100); [0034] Determine, on the basis of said data, at least one reference signal of the driving force F=F(s) and a reference signal of the speed =(s) for the vehicle V, in function of the position s of the vehicle itself, along the route to travel, through a first optimisation process (step 200) (for which the speed =(s) is generally not a constant cruising speed); and [0035] determine at least one reference signal of the engaged gear (s) (or of the corresponding gear changes u.sub.sh(s)) of the vehicle V, in function of the position s of said vehicle along the route to travel. Through a second optimisation process (step 300).

[0036] The second optimisation process (step 300) is subsequent to the first optimisation process (step 200) and receives as input both the data relating to the vehicle V and the data relating to the route to travel, and the reference signal of the driving force F and the speed reference signal , processed during the first optimisation process (step 200).

[0037] With particular reference to step 100 of the preferred embodiment of the method according to the present invention (see FIG. 3), it envisages to provide information relating to the vehicle, which may include: information relating to its mass, to its aerodynamic resistance and to the rolling friction (for example in the form of corresponding coefficients), to the pattern of the maximum torque developed by the engine of the vehicle and, as will be seen below, information relating to fuel flow maps (for example, as function of the torque developed by the engine, and the speed of the same engine) and the transmission ratios.

[0038] The step 100 of the method according to the present invention also provides for the supply of certain data relating to the route that the vehicle V will travel, including, for example: [0039] a current geographical position of the vehicle, provided for example by means of a suitable Global Positioning System (GPS, GLONASS or other similar system), [0040] a destination geographical position, [0041] information relating to the route between the starting geographical position of the vehicle (which is advantageously the geographical position of departure of the route to travel) and the destination geographical position (including the path of the roadway of the route to travel), such as the slope or the altitude of the roadway that the vehicle must travel, [0042] data corresponding to the speed limits along the route, as well as [0043] a desired arrival time, inserted for example by the driver of the vehicle, and [0044] other optional data such as, for example, an average speed, the maximum and minimum speed desired by the driver as well as the degree of adherence to those speeds.

[0045] The information that can be inputted by the driver are used by the method according to the present invention, as will be seen hereinafter, to take into account also his driving style and the conditions of the road to travel.

[0046] These data are used to perform a first optimisation process (step 200), after suitable treatment (for example comprising a filtering step and re-sampling and/or other treatment of any suitable type), which provides for determining a reference signal of the vehicle speed (indicated in FIG. 7 with reference V) that optimizes the fuel consumption during the route to travel, as well as a reference signal of the driving force F of the vehicle along that route.

[0047] More in particular, said first optimisation process comprises an analytical process of minimization of a first cost function which is a cost function J(F) of the driving force F, required to drive the vehicle V of mass m along a route having global length S.

[0048] The minimization of the cost function J(F) with respect to the drive force F is given by the relation:

[00001]$\begin{array}{cc}\underset{F}{\min }.Math.\left\{J\left(F\right)\right\}=\underset{F}{\min }.Math.\left\{{}_0^S.Math.F^2.Math.\mathrm{ds}\right\}& \left(1\right)\end{array}$

with the following state equation in the time domain:

[00002]$\begin{array}{cc}m.Math.\frac{\mathrm{dv}}{\mathrm{dt}}=F-F_{\mathrm{grade}}-F_{\mathrm{drag}}-F_{\mathrm{roll}}& \left(2\right)\end{array}$

and the constraints:

[00003]$\begin{array}{cc}{v}_{\min }v{v}_{\max }& \left(3\right)\\ F{F}_{\max }& \left(4\right)\\ {}_0^S.Math.\frac{1}{v}.Math.\mathrm{ds}{T}_{\max }=\frac{S}{{\stackrel{\_}{v}}_{\mathrm{des}}}& \left(5\right)\end{array}$

wherein is the vehicle speed, S is, as said above, the global length of the route to travel .sub.min is the minimum speed and .sub.max is the maximum speed of the vehicle (optionally provided by the driver and/or from the data relating to the route the vehicle V will travel, such as those relating the speed limits), F.sub.max is the available maximum driving force of the vehicle (where the maximum driving force F.sub.max is the control variable), T.sub.max is a desired maximum travel time (optionally provided by the driver) and .sub.des is the set desired average speed, even this optionally provided for by the driver of the vehicle V.

[0049] The minimization (1) of the cost function J(F) is aimed at minimising the energy used to drive the vehicle having mass m along a route of finite length S, given the constraints for the state equation (2) relating to the speed , maximum driving force F.sub.max and desired maximum travel time T.sub.max mentioned above in equations (3) to (5). Moreover, the state equation (2) also takes into account the fact that the vehicle V is subject to the following exogenous forces:

F.sub.grade=m g sin() (6)

F.sub.drag=.sub.airC.sub.xA.sup.2 (7)

F.sub.roll=m g cos()[C.sub.r0+C.sub.r1.sup.2](8)

[0050] Due to the slope of the road along which the vehicle is travelling (F.sub.grade), to its aerodynamic resistance F.sub.drag and to the rolling friction F.sub.roll, respectively, which depend, other than the global mass m of the vehicle V, also from the gravitational acceleration g, from the angle corresponding to the slope of the road (which in general can be variable along the route, for which =(s)), from the ambient air density .sub.air (which in general can be variable along the route, for example because the route runs through different altitudes for which .sub.air=.sub.air.sup.(s)), from the aerodynamic drag coefficient C.sub.x, from the frontal area A of the vehicle V, from the first coefficient C.sub.r0 of rolling resistance (which in general can be variable along the route, for example due to the progressive wear of the tires and/or the different type of asphalt or soil and/or the potentially different weather conditions, for which C.sub.r0=C.sub.r0(S)) and from the second coefficient C.sub.r1 of rolling resistance (which in general can be variable along the route, for example as a result of the progressive wear of the tires and/or the different type of asphalt or soil and/or the potentially different weather conditions, so that C.sub.r1=C.sub.r1(s)).

[0051] Since the cost function to J(F) to be minimized and the constraint relating to maximum travel time T.sub.max to be met reflect both objectives in the space domain, the status equation (2) can be rewritten in the space domain (see equation (9) below) and the state variable (i.e. the speed), is replaced by the kinetic energy K so that the equations (2)-(5), using the definitions of exogenous forces acting on the vehicle V, given by equations (6)-(8) can be rewritten as follows:

[00004]$\begin{array}{cc}\frac{\mathrm{dK}}{\mathrm{ds}}=F-m.Math..Math.g.Math..Math.\sin \left(\right)-\frac{{}_{\mathrm{air}}.Math..Math.C_x.Math..Math.A}{m}.Math.K-m.Math..Math.g.Math..Math.\cos \left(\right)\left[C_{r.Math..Math.0}-C_{r.Math..Math.1}.Math.\frac{2.Math..Math.K}{m}\right]& \left(9\right)\\ .Math.\frac{1}{2}.Math.m.Math..Math.v_{\min }^2=K_{\min }K{K}_{\max }=\frac{1}{2}.Math.m.Math..Math.v_{\max }^2& \left(10\right)\\ .Math.F{F}_{\max }& \left(11\right)\\ .Math.{}_0^S.Math.\sqrt{\frac{m}{2.Math..Math.K}}.Math.\mathrm{ds}{T}_{\max }=\frac{S}{{\stackrel{\_}{v}}_{\mathrm{des}}}& \left(12\right)\end{array}$

where K.sub.min is the kinetic energy associated with the minimum speed .sub.min and K.sub.max is the kinetic energy associated with the maximum speed .sub.max.

[0052] The integral constraint (12), to be suitably treated, requires the introduction of the additional state variable t(s), i.e. the time required to travel the space s, thus resulting in the addition of a further constraint equation:

[00005]$\begin{array}{cc}\frac{\mathrm{dt}}{\mathrm{ds}}=\sqrt{\frac{m}{2.Math..Math.K}}& \left(13\right)\end{array}$

according to which, the integral constraint (12) can be explicitly rewritten as:

t(S)=T.sub.max (14)

[0053] This stated, based on the method according to the present invention, the minimization of the cost function J(F), while respecting the constraints described above can be obtained according to the method of the present invention, by means of the minimization of the integral of the Lagrangian, function of the driving force F, given by:

[00006]$\begin{array}{cc}\left(F\right)={}_0^S.Math.\left\{F^2+a\left(s\right)\left[\frac{\mathrm{dK}}{\mathrm{ds}}-F+m.Math..Math.g.Math..Math.\sin \left(\right)+\frac{{}_{\mathrm{air}}.Math..Math.C_x.Math..Math.A}{m}.Math.K+m.Math..Math.g.Math..Math.\cos \left(\right).Math.\left(C_{r.Math..Math.0}-C_{r.Math..Math.1}.Math.\frac{2.Math..Math.K}{m}\right)\right]+g_1\left(K\right)+g_2\left(K\right)+g_3\left(F,v\left(K\right)\right)+\left(\frac{\mathrm{dt}}{\mathrm{ds}}-\sqrt{\frac{2.Math..Math.K}{m}}\right)\right\}.Math.\mathrm{ds}+b_1\left(K\left(0\right)-K_0\right)+b_2\left(t\left(0\right)-0\right)+{\left(t\left(S\right)-T_{\max }\right)}^2.Math.& \left(15\right)\end{array}$

where the penalisation functions g.sub.1, g.sub.2 and g.sub.3 are defined by the following equations:

[00007]$\begin{array}{cc}{g}_1\left(K\right)=\{\begin{array}{cc}{c_1\left(K-K_{\max }\right)}^2& \mathrm{if}.Math..Math.K>K_{\max }\\ 0.Math.& \mathrm{if}.Math..Math.K{K}_{\max }\end{array}& \left(16\right)\\ g_2\left(K\right)=\{\begin{array}{cc}{c_2\left(K_{\min }-K\right)}^2& \mathrm{if}.Math..Math.K<K_{\min }\\ 0.Math.& \mathrm{if}.Math..Math.K{K}_{\min }\end{array}& \left(17\right)\\ g_3\left(F,v\right)=\{\begin{array}{cc}{c_3\left(F-F_{\max }\left(v\right)\right)}^2& \mathrm{if}.Math..Math.F>F_{\max }\\ 0.Math.& \mathrm{if}.Math..Math.F{F}_{\max }\end{array}& \left(18\right)\end{array}$

where c.sub.1, c.sub.2 and c.sub.3 are penalisation coefficients.

[0054] In the equations above, the available maximum driving force F.sub.max, according to the method of the present invention, can be obtained for example by means of an analytical approximation function of the envelope of the maximum driving force for each gear, as shown in FIG. 4, given by the sum of three Gaussian functions according to the relation:

[00008]$\begin{array}{cc}{F}_{\max }\left(v\right)=\underset{j=1}{\overset{3}{.Math.}}.Math..Math.{}_j.Math.e^{-{\left(\frac{v-{}_j}{{}_j}\right)}^2}& \left(19\right)\end{array}$

wherein is the vehicle speed and .sub.1 . . . 3, .sub.1 . . . , .sub.1 . . . 3 are nine constant coefficients.

[0055] It must be considered that the Gaussian functions necessary to the representation of the envelope of the maximum driving force for each gear (i.e. as a function of speed ) could be any other number P greater than or equal to 2 (e.g. 2, 4 or 5), and the number of the constant coefficients .sub.1 . . . P.sub.1 . . . P, .sub.1 . . . P, is equal to 3 times P. Moreover, it must be kept in mind that the available maximum driving force F.sub.max could be represented in analytical form with other basic functions rather than with Gaussian functions, for example by using a piecewise linear function.

[0056] In the integral of the Lagrangian (15), the presence of the constraints of equality (9) and (13), according to the method of the invention requires the introduction of added fields a(s) and ; in addition, in the integral of the Lagrangian (15) is the amplitude of a quadratic penalisation term, to meet the constraint (14), while b.sub.1 and b.sub.2 are Lagrange multipliers, which however do not have any operational significance.

[0057] The method according to the invention performs an algorithm of computation of the Lagrangian gradient (15). In the preferred embodiment of the method according to the invention, that algorithm is adjoint-based iterative and, as shown in FIG. 5, follows the method of the conjugate gradient. That algorithm provides for the computation of the gradient of the integral of the Lagrangian (15)

[00009]$\frac{}{.Math..Math.F}$

(step 250) using the added field, and the minimization of the integral of the Lagrangian is in fact performed by the method of the conjugate gradient by iteratively using the gradient

[00010]$\frac{}{.Math..Math.F}$

thus calculated, to determine the search direction along which to search for the minimum (step 260). In particular, the penalisation coefficients c.sub.1, c.sub.2 and c.sub.3 of equations (16)-(18) can be constant or optionally progressively increased at successive iterations; similarly, the amplitude of the quadratic penalisation term can be a constant value or optionally modifiable, dynamically, at successive iterations. Once the search direction is determined, the algorithm provides to perform the line minimization (step 270) and repeat the iterations until at least one stop criterion of the same iterations is met (step 280), which indicates that a minimum of the Lagrangian integral has been reached.

[0058] Since this is an iterative method, it provides to give likely initial values K(s) and t(s) for the variable F(s), for example by solving equations (9) and (13) (step 251).

[0059] For the gradient computation, the added field is then determined (step 252), by solving the following added equation:

[00011]$\begin{array}{cc}\frac{}{.Math..Math.t}=0.fwdarw.\{\begin{array}{c}\frac{d.Math..Math.}{\mathrm{ds}}=0.Math.\\ \left(S\right)=2.Math.\left(t\left(S\right)-T_{\max }\right)\end{array}& 20)\end{array}$

derived from the first disturbance of the Lagrangian integral (15).

[0060] The solution of the above differential equation is trivial and highlights that must be constant and equal to its final condition.

[0061] At the next step (step 253), the field a(s)is determined, by solving the following the second added equation:

[00012]$\begin{array}{cc}\frac{}{.Math..Math.K}=0.fwdarw.\{\begin{array}{c}\frac{\mathrm{da}}{\mathrm{ds}}=-2.Math.a\left(\frac{1}{2}.Math.\frac{{}_{\mathrm{air}}.Math..Math.C_x.Math..Math.A}{m}+C_{r.Math..Math.1}.Math..Math.g.Math..Math.\cos \left(\right)\right)+\frac{{\mathrm{dg}}_1}{\mathrm{dK}}+\frac{{\mathrm{dg}}_2}{\mathrm{dK}}+\frac{1}{2}.Math..Math.\sqrt{\frac{m}{2}}.Math.\frac{\sqrt{1}}{K^3}.Math.\\ a\left(S\right)=0.Math.\end{array}& \left(21\right)\end{array}$

from which (step 254) the gradient of the Lagrangian integral is obtained based on the relation:

[00013]$\begin{array}{cc}\frac{}{.Math..Math.F}=-a+2.Math.F+\frac{{\mathrm{dg}}_3}{\mathrm{dF}}& \left(22\right)\end{array}$

[0062] The gradient of the Lagrangian integral thereby calculated is used in step (step 260) for determining the search direction .sub.k, in the space of decision variables, wherein, at the k-th iteration, the search direction .sub.k is determined based on the relation:

[00014]$\begin{array}{cc}\{\begin{array}{cc}{p}_k=-\frac{}{.Math..Math.F}.Math.{}_k.Math.& \mathrm{if}.Math..Math.k=1\\ p_k=-\frac{}{.Math..Math.F}.Math.{}_k.Math.+{}_k.Math.p_{k-1}& \mathrm{if}.Math..Math.k>1\end{array}& \left(23\right)\end{array}$

wherein: [0063] in the first iteration (k=1), as search direction the opposite of the gradient, i.e. the steepest descent direction is used, and [0064] In subsequent iterations (k>1) the conjugate gradient is determined, where represents the momentum of the search direction at the previous iteration and is determined by the formula of Polak-Ribire:

[00015]$\begin{array}{cc}{}_k=\max \left(0;\frac{{\left(-\frac{}{.Math..Math.F}.Math.{}_k\right)}^T.Math.\left(-\frac{}{.Math..Math.F}.Math.{}_k.Math.+\frac{}{.Math..Math.F}.Math.{}_{k-1}\right)}{{\left(-\frac{}{.Math..Math.F}.Math.{}_{k-1}\right)}^T.Math.\left(-\frac{}{.Math..Math.F}.Math.{}_{k-1}\right)}\right)& \left(24\right)\end{array}$

[0065] IT is to be considered that can be determined by other formulas, such as for example the Fletcher-Reeves or Hestenes-Stiefel formula.

[0066] Once the search direction has been determined, the line minimization is performed by calculating (step 270) the length of the pitch (h.sub.k) of the driving force F that ensures the maximum reduction of the cost function

F.sub.k+1=F.sub.k+h.sub.k.sub.k (25)

by iteratively solving the following not linear equation, for the scalar h,

[00016]$\begin{array}{cc}\frac{\left(F+h.Math..Math.p\right)}{h}=\frac{\left(F+h.Math..Math.p\right)}{\left(F+\mathrm{hp}\right)}.Math.\frac{\left(F+h.Math..Math.p\right)}{h}=\left(2.Math.\left(F+\mathrm{hp}\right)-a+\frac{{\mathrm{dg}}_3}{\mathrm{dF}}\right)=0& \left(26\right)\end{array}$

for example, by the method of Brent described by R. P. Brent in Chapter 4: An algorithm with Guaranteed Convergence for finding a zero of a Function, Algorithms for minimization without derivatives, Englewood Cliffs, N.J.: Prentice-Hall, 1973, ISBN 0-13-022335-2, or also by other methods such as Dekker algorithm.

[0067] Once the length h.sub.k of the pitch is determined, the control variable F, is updated (step 271) and at step 280 it is checked whether a criterion for stopping the iteration is met or not.

[0068] The stop criterion of the algorithm according to the preferred embodiment of the method of the present invention is given by:

[00017]$\begin{array}{cc}\{\begin{array}{c}\frac{J\left(F_k\right)-J\left(F_{k-1}\right)}{J\left(F_{k-1}\right)}{\mathrm{.Math.}}_J\\ \frac{t\left(S\right).Math.{}_k.Math.-t\left(S\right).Math.{}_{k-1}}{t\left(S\right).Math.{}_{k-1}}{\mathrm{.Math.}}_t\end{array}& \left(27\right)\end{array}$

wherein .sub.j and .sub.t are two predetermined threshold values.

[0069] As shown, the stop criterion is given by two conditions: a first condition, linked to the cost function J(F), which is considered to be minimized when the relative variation of the same between an iteration and the previous one is less than or equal to .sub.j; the second stop criterion, which also must be met in order for the iterations to stop, concerns the duration t(S) of the whole route, whose variation between an iteration and the previous one cannot be greater than .sub.t.

[0070] Based on the calculations mentioned above, the first optimization process (step 200 of FIGS. 1 and 5) provides as output a reference signal of the speed (s) of vehicle V and a reference signal of the driving force F(s), as a function of its position s along the route to travel.

[0071] Based on the results outputted from the first optimisation process, the method according to the present invention provides for performing a second optimisation process (step 300 of FIG. 1) of a second cost function, which is different from the cost function of the first optimisation process, by the method of the dynamic programming (also known as dynamic programming or dynamic optimization), which is a cost function J.sub.2 of fuel consumption. The cost function J.sub.2is expressed as the sum, on M length segments s.sub.2 that make up the route of length S, of the costs due to the fuel consumption along each segment i-th having length s.sub.2 and penalisation factors, aimed at penalising the gear changes, wherein an additional penalisation is applied when a gear change occurs on uphill sections of the route.

[0072] The block diagram of FIG. 6 schematically illustrates the second optimisation process of the preferred embodiment of the method according to the present invention. In particular, the second optimisation process, based on the reference signal of the speed (s) of the vehicle and a reference signal of the driving force F(s) obtained from the first optimization process (step 200 of FIGS. 1 and 5) of the method and the data mentioned above relating to the type of the route to travel, in particular with reference to the slope or altitude of the roadway to travel and the vehicle V (mass, fuel flow maps, transmission ratios or the like).

[0073] More in detail, this second optimisation process performs the minimization of the cost function J.sub.2 as a function of the gear changes , associated with the fuel consumption of the vehicle, and given by the relation:

[00018]$\begin{array}{cc}\underset{{}_i{}_i}{\min }.Math.\left\{J_2\right\}=\underset{{}_i{}_i}{\min }.Math.\left\{\underset{i=1}{\overset{M}{.Math.}}.Math..Math.\left(\frac{W_{f,i}\left({}_i\right).Math..Math..Math.s_2}{v_i}+{}_{1,i}\left({}_i,m\right).Math..Math.u_{\mathrm{sh},i}.Math.+{}_2.Math..Math.u_{\mathrm{sh},i}.Math.\right)\right\}& \left(28\right)\end{array}$

with the following state equation:

.sub.i+1=.sub.i+u.sub.sh,i (29)

and the constraints that follow:

.sub.i.sub.i (30)

.sub.i:n.sub.e,minn.sub.e(.sub.i, K.sub.i)n.sub.e,max (31)

u.sub.sh,i[2; 1, 0; +1; +2](32)

wherein:

[00019]$\begin{array}{cc}M=\frac{S}{.Math..Math.s_2}& \left(33\right)\\ {}_{1,i}=\{\begin{array}{cc}m.Math.{}_i& \mathrm{if}.Math..Math.{}_i>0\\ 0.Math.& \mathrm{if}.Math..Math.{}_{i}0\end{array}& \left(34\right)\end{array}$

and wherein: [0074] u.sub.sh,i is the gear change signal in segment i and represents the control variable; [0075] M is the number of addenda of the cost function J.sub.2 (equal to the number of length segments s.sub.2 that compose the route of length S); [0076] W.sub.f,i Is the fuel flow rate (mass/time) in segment i; [0077] s.sub.2 is the constant length of each one of the M segments (even if in other embodiments of the method according to the invention the length s.sub.2,i of each segment i can be in general different from the lengths of the other segments); [0078] .sub.1,i is a first penalisation coefficient that is variable as a function of the mass m of the vehicle and the slope .sub.i of the segment i; [0079] .sub.2 is a second penalisation coefficient which is constant; [0080] .sub.i Is the gear of the vehicle in the segment i, comprised between the set of gears .sub.i, and represents the state variable (in general, the set of gears .sub.i can change section by section), [0081] n.sub.e Is the speed of the engine, between a maximum value and a minimum value n.sub.e,max and n.sub.e,min.

[0082] It must be considered that, similarly to the prior art systems, when exceptional or emergency conditions occur that require the engagement of gears outside the set .sub.i or the operation of the brakes, both in automatic mode and for direct intervention of the driver, the method according to the invention stops its execution (i.e. is turned off).

[0083] A sufficiently high value must be assigned to penalisation factors .sub.1,i and .sub.2 so that the gear changes that are not needed are properly penalised.

[0084] Moreover, it should be noted that equation (32) limits the set of possible operations between: maintain the gear (0), increase one or two gears (+1, +2) or decrease one or two gears (1, 2).

[0085] The term W.sub.f,i the first contribution of the addenda of the sum of the cost function J.sub.2, can be expressed as a function of the torque T.sub.e developed by the engine and the speed n.sub.e of the same engine, according to the fuel flow map W.sub.f(T.sub.e, n.sub.e) of the engine of the vehicle. As it can be easily understood, other engine parameters can be alternatively or additionally used to determine the fuel flow rate W.sub.f,i, in a totally obvious way for the person skilled in the art.

[0086] According to the present invention, however, the speed n.sub.e and the torque T.sub.e developed by the engine can be expressed as a function of the gear through the relations:

[00020]$\begin{array}{cc}{T}_e\left(F,\right)=\frac{F.Math..Math.R_w}{i\left(\right)}& \left(35\right)\\ n_e\left(v,\right)=\frac{v.Math..Math.i\left(\right)}{R_w}& \left(36\right)\end{array}$

where i() is the global transmission ratio, with the engaged gear , and R.sub.w is the radius of the wheel of the vehicle V.

[0087] From equations (35) and (36) it is therefore clear that, once known F(s) and (s), supplied in output from the first optimisation process of the method according to the present invention, the cost function J.sub.2, different from the cost function of the first optimisation process, associated with the fuel consumption during the route, is a function of only the engaged gear .sub.i in each segment of the route to travel.

[0088] Therefore, in output from the second optimization process of the method according to the present invention, there is obtained a reference signal of the engaged gear (s) (or the gear changes u.sub.sh(s)) of the vehicle V, along the route to travel; advantageously, although the signal of the engaged gear .sub.i refers to the segment i, nevertheless the reference signal of the engaged gear (s) (or the gear changes u.sub.sh(s)) is provided as a control signal at a certain travelled distance s.

[0089] Having said that, in view of the above, it is clear how the method according to the present invention, that implements two successive processes of optimization, each one based on a cost function that is different from that of the other optimization process, is able to obtain reference signals of the speed (s) of a vehicle V and the gear engaged (s) (or the gear changes u.sub.sh(s)) along the entire route travel so that, given the current position s of the vehicle V along the way, it will be possible to determine in a completely unique way the corresponding values of the speed (s) and of the reference gear (s).

[0090] A system for obtaining reference signals for control systems of a vehicle V, as function of the geographic position of the same vehicle along a route travel will be now described purely by way of example and with reference to FIG. 2.

[0091] Such a system 2 comprises at least one device 3 for collecting data, configured to output data relating vehicle V and a route to travel, as well as any data that can be inserted by the driver of the vehicle V and data related to the engine of the vehicle V. In FIG. 2, said at least one device 3 comprises: a memory 3a that stores data related to the vehicle V, such as the total mass m of the vehicle V (the value of which can possibly be updated by weight sensors positioned on the specific components of the vehicle V), an aerodynamic resistance coefficient C.sub.x, frontal area A of the vehicle V, the first and the second coefficient C.sub.r0 and C.sub.r1 rolling resistance of the tires of the vehicle V (which can be updated based on the progressive wear of the tires); a satellite navigator 3b (provided for example with GPS sensor and/or GLONASS) that provides data relating to the position of the vehicle and the route to travel, such as the current position of the vehicle, a geographical destination position, information relating to the route between the position of the vehicle V and the destination (such as the slope or the altitude of the roadway that the vehicle V must travel), data corresponding to the speed limits along the way, and optionally information relating to traffic conditions and the weather condition; an input/output interface 3c through which the driver of the vehicle can enter user data, such as for example a desired arrival time, an average speed, the maximum and minimum speed required by the driver, and optionally a degree of adherence to the speed set by the driver; a 3d system configured to provide in output data relating to the engine, optionally comprising one or more sensors 3d configured to provide in output data related to the torque developed by the engine and/or the speed of the engine itself (such sensors may be those already present in the vehicle and the data can be retrieved from the communication board system, such as the CAN-BUS) and a memory which stores the pattern of the maximum torque developed by the engine of the vehicle, fuel flow maps (for example as function of the torque developed by the engine and the engine speed) and the transmission ratios. Other embodiments of the system according to the invention can comprise additional sensors, for example for detecting a ambient air density and/or weather conditions.

[0092] The system 2 also comprises at least one first module 4 for optimizing the driving force, operatively connected to said at least one device 3, and configured to determine, by means of a first optimisation process such as that described above with reference to FIG. 5, a driving force signal F(s) and a speed signal (s) to the vehicle, based on the position s of the vehicle V, along the route.

[0093] The system 2 further comprises at least one second module 5 for optimizing the fuel consumption, arranged downstream of the first module 4 for optimizing the driving force and operatively connected to said first module 4 and to said at least one device 3 for collecting data. The second optimization module 5 is configured to determine a reference signal of the reference gear (s) (or the gear changes u.sub.sh(s)) of the vehicle V, along the route to be travelled, through a second optimisation process.

[0094] In particular, the first and the second module 4 and 5 can be implemented through two respective processing units, provided for example with at least one microprocessor or micro-controller, or through a single processing unit, provided for example with at least one microprocessor or micro-controller, configured to perform the method according to the invention.

[0095] Said at least one device 3 is configured to provide the first module 4 and the second module 5 with the data illustrated above with reference to the optimization processes of FIGS. 3, 5 and 6.

[0096] The first and the second module 4 and 5 receive from the satellite navigator 3b the information regarding to the route to travel i.e. the slope or the altitude of the roadway along the route to travel, which route is set by the driver by means of the interface input/output 3c, and speed limits provided along that route (and optionally the information regarding the traffic conditions and the weather condition) based on which the first and the second module 4 and 5 perform the first and the second optimization process, respectively, as shown in FIGS. 3 and 6 (wherein such information are cumulatively indicated with the reference number 30).

[0097] The information 30 provided by the satellite navigator 3b may also have been stored beforehand in the same satellite navigator, or may be processed on case by case basis, optionally through a connection by means of a wireless radio communication network (wireless) according to one communication protocol suitable for one or more remote servers or to a remote processing unit.

[0098] FIGS. 3 and 6 also show the information 35 supplied by the memory 3a and by the system 3d to the first and the second module 4 and 5 (even if the information 35 provided by the system 3d is actually used only by the second module 5), i.e. the information comprising data relating to the mass, the aerodynamic resistance, the rolling friction, to pattern of the maximum torque developed by the engine of said vehicle, to the corresponding fuel flow maps, the fuel flow velocity (i.e. the flow rate of fuel) and the gear ratios.

[0099] Now, as it can easily be imagined one system like the one described above can be totally integrated into one vehicle or it can be only in part. Some versions are in fact foreseen, which hereinafter will be presented with reference to FIG. 7.

[0100] According to three main versions, the optimization modules 4 and 5 can be installed on board of the vehicle V or may be external to the vehicle V, for example installed, by means of a suitable software application on a computer device 6 in the cloud, with which they communicate through a suitable communication protocol and/or can be installed on an external device 3, for example a smartphone, tablet, personal computer, operatively connected by means of suitable communication interfaces and in a manner known to a person skilled in the art, to the other components of the system.

[0101] At least one of the above versions provides that the data of the vehicle V concerning in particular the engine are stored on board the vehicle itself, for example on an adjusted memory 7 of the non-volatile type. Such data, as said above, shall be measured or made available by the manufacturers of the vehicle V and/or the engine of the vehicle V as part of the specifications of the vehicle V and/or of the engine itself.

[0102] At least one of the above versions provides that at least the constant data of the vehicle V and/or the engine of the vehicle V (with constant a parameter that does not vary during the life of the vehicle V and/or the engine of the vehicle V is intended), are stored on board, for example by means of an adjusted non-volatile memory 7 suitable for the purpose. The parameters that the vehicle V and/or the engine of the vehicle V which vary during the life of the vehicle V and/or the engine, may instead be estimated in-line through estimation techniques such as Kalman filters, state observers, state estimators and the like. Reference is made in this case, to: [0103] parameters of the vehicle V, such as the aerodynamic resistance coefficient, the rolling friction coefficient, the masses, the inertia, the characteristics of the chassis or other; and [0104] engine parameters, such as the torque and/or power that can be developed, the fuel consumption maps and any other indication of the fuel quantity as a function of the operating conditions.

[0105] Alternatively, the method and the relative system according to the present invention can provide that the parameters of the vehicle V and/or engine, among those previously mentioned, which vary during the life cycle of the vehicle V and/or of the engine can be estimated not in-line and not on board the vehicle itself by an external system 8, for example during maintenance operations of the vehicle V and/or of the engine and then subsequently stored on board in an adjusted non-volatile memory, optionally on a external device 3, such as for example a smartphone, tablet or the like.

[0106] According to another version of the method and the corresponding system of the present invention, both the constant parameters and the ones that are not constant of the vehicle V and/or the engine may be stored on an external device 3 (for example a smartphone or tablet or the like), independently from the fact that they are constants or that they are estimated in-line.

[0107] In view of the above it can be clearly understood that a system like the one described above can be integrated in a control system of a vehicle (which is also an object of the present invention), including, inter alia, one or more control apparatuses configured to control the speed () and the gear engaged () (or the gear change (u.sub.sh) for the vehicle V, operatively connected downstream of the system 2, from which they receive, in use, the reference signals determined according to the method described above.

[0108] Said one or more apparatuses for controlling the speed and the gear engaged (or of the gear changed u.sub.sh) of the vehicle V, include, for example, a cruise control system which regulates a power-train, or a control unit of the gearbox or any other control device that contributes directly or indirectly to the control of the fuel consumption, the reduction of emissions, the safety while driving and similar.

[0109] The method and the system for obtaining reference signals described above fulfil the objects indicated above.

[0110] With particular reference to the method, experimental studies have shown that, under equal conditions with respect to traditional methods, the method according to the present invention provides for reference signals for vehicles control systems that allow to adjust the speed of a vehicle V, improving performance from the consumption point of view and reducing, at the same time, the travel times. In addition, the two optimization processes performed in sequence, and not at the same time, by the method according to the present invention, allows to considerably reduce the computational complexity of the method, thereby making it also more efficient from the computational point of view.

[0111] Below the results are shown of some tests carried out and illustrated in FIG. 8. FIG. 8 shows a comparison between the behaviour of the proposed method (OPT) as compared to the traditional method (CC), in the case of a vehicle with a full load, running the route whose profile in terms of altitude is reported in graph (a).

[0112] The traditional method implemented by a fixed point cruise control system (continuous line), acts on the gas control, on the basis of the residual between a reference speed (constant) and an actual speed. The gear changes are activated by the prediction of the current load, considering the maximum efficiency range of the engine in terms of revolutions per minute (rpm). Such a logic may generate, in some circumstances, an excessive series gear changes, leading to some disadvantages such as the loss of momentum, during the time interval in which the clutch disengages the engine from the wheels, with a consequent increase in fuel consumption on uphill sections of the route, and increased stresses on the transmission.

[0113] The method according to the present invention, on the other hand, (see the dashed lines in (b), (c) and (d) of FIG. 8), dynamically changes the reference speed and the engagement of the gears, based on the topography of the route, adapting to the profiles of the slope and, therefore, saving fuel.

[0114] By way of example, with reference to FIG. 8(b), it should be noted that the speed is slightly increased at approximately 3000 m, at the beginning of an uphill section, and at the same time the gear is predictively downshifted (Errore. L'origine riferimento non stata trovata. (d)). The increased momentum is beneficial to tackle the next section at the higher speed, where the slope of the roadway increases rapidly. Furthermore, the gear change in advance allows to avoid to further change gear during the ascent. At the end of the uphill section, the traditional method provides for the quick return to the reference speed, while the method according to the present invention leaves the vehicle slowing down gently to exploit the next section of descent (4500-6000 m) for regaining speed.

[0115] At the end of the slope, the method according to the present invention maintains a higher speed compared to the traditional method, to compensate for the time lost in the preceding section.

[0116] Since this section is approximately plane, the method according to the present invention provides for a slight increase with reference to consumed fuel, which however is compensated by a gain in terms of journey time.

[0117] In this regard, in FIG. 9 the fuel consumption (FC) is reported together with the variations of travel times (time), in the case of a vehicle of 44t at full load with different reference speeds, both for the traditional method (with reference speed constantwhose graphs are indicated with Traditional Cruise Controller in FIG. 9), both for the method according to the present invention (whose graphs are indicated with Optimized in FIG. 9).

[0118] The variations shown in the graph are related to the fuel consumption and the travel times obtained with a traditional method at reference speed set equal to 80 km/h on a stretch of road actually existing over more than 100 km.

[0119] As it can be noted, the method of the present invention allows, at equal reference speed, to save fuel and time.

[0120] Moreover, it is also possible to increase the reference speed to 82 km/h to save nearly 3% in terms of journey time, while maintaining a saving of fuel consumed of 1.37%. On the contrary, the increase of the reference speed to 82 km/h entails an increase in fuel consumption of 2.43% for the traditional method known as cruise control at a fixed point to, although it has a reduction in journey time of about 2%.

[0121] On the other hand, the reduction of the reference speed to 79 km/h allows for the method according to the present invention to save over the 4.04% of fuel, increasing the duration only of 0.70%.

[0122] In the foregoing the preferred embodiments were described and some modifications of this invention have been suggested, but it should be understood that those skilled in the art can make modifications and changes without departing from the relative scope of protection, as defined by the appended claims.