Online Energy Management of City Cars with Multi-objective Linear Parameter-varying L 2 -gain Control

This work aims at online regulating transient current out of the batteries of small-sized electric cars that transport people and goods around cities. In a city with heavy traffic, transient current dominates the energy economy and propulsion capability, which are in opposition to each other. In order to manage the trade-off between energy consumption per distance and propulsion capability in transience, the authors improve on previous work on multi-objective linear parameter-varying (LPV) L2-gain control. The observer embedded into this multi-objective controller no longer assumes Kalman-filtering structure, and structural conservatism is thus removed. A full-spectrum set of experiments is performed. The results reveal that the feedback design significantly improves energy-motion management.


Introduction
This work is concerned with online energy management of small-sized electric cars that can carry people and goods around a city with heavy traffic without discharging pollutants or green-house gases.There are currently four types of energy management systems that can be employed in city cars: (a) management of active power sources; (b) assistance with passive energy storage; (c) restructuring of power converters; and (d) implementation of control rules for the energy economy, as discussed in the sequel.
Most power sources available to the car industry suffer energy high ratios, either to cost or to weight.For example, lithium-ion batteries, fuel cells [1], high-capacity ultracapacitors, and ultrahigh-speed flywheels are too costly, and lead-acid batteries are too heavy.Therefore, it has become important for the industry to manage batteries to reduce consumption of electric energy.For example, researchers and practitioners put much effort into online measuring the state-of-charge and/or state-of-health of a diversity of batteries, and then regulating current outputs to decrease energy consumption per distance, while increasing the lifespan of the batteries [2][3][4][5][6][7].Besides, parallel connection of ultracapacitors with batteries can boost energy efficiency.For example, ultracapacitors installed between batteries and power converters can increase the energy harvested in braking phases from mechanical kinetics to electricity [8,9].Since capacitors are passive elements, they can efficiently absorb the impulsive currents out of boost choppers; otherwise, impulsive currents input to batteries merely become heat being dissipated instead of electricity being stored.There is the possibility that supercapacitors could also be employed for fast charging of city cars [10].Furthermore, ultracapacitors deliver initial overshoots of currents to start motors, and thus are indispensable to fuel cell-powered city cars [11,12].However, most city cars are lightweight and possess low horsepower, thus the amount of energy conserved in braking phases is limited, and the initial torque is usually large enough to start the motor.
Restructuring power converters can lead to vital reduction of electricity-to-heat, which has been well investigated by researchers and practitioners in the field of power-electronics control.For example, some designed multiple bi-directional DC/DC power converters to "soften" switching of currents, and RCD snubbers that network diodes with resistors and capacitors to reduce energy consumption from high-frequency switching [13][14][15].Some have embedded microcontrollers into power converters for optimal timing of switching to suppress electricity loss [16][17][18][19].In types (a), (b) and (c), energy management is focused on batteries and power converters, in which additional components are needed for the fulfillment of efficient energy.However, in type (d), a set of dynamic rules is programmed into the controller, an already existing component, to boost energetic economy [20][21][22][23], with the advantage of no additional hardware.This paper discusses type (d).
A city car is usually used for transportation within a heavy-traffic city; therefore, its transient current plays the leading role in economy of propulsion energy and safety of accelerating out of a busy crossing.Controlling current overshoots, which simultaneously dull the motion, can reduce energy dissipation per distance.This paper develops a multi-objective control to manage the trade-off between energy economy and motion dexterousness.
During the modeling process, we found through lab experiments that linear time-invariant (LTI) parameterization of city-car dynamics is unable to capture overshoots and the times the currents out of batteries rise.However, a city car is for transportation inside traffic jams, so its transient current plays the key to energy economy and driving safety.To this, we develop linear parameter-varying (LPV) modeling and identification to remedy such a situation, following the previous works in [24][25][26][27][28][29][30].Here, the tire speed is taken as the slow-time state variable and simultaneously the scheduling parameter.Then, guided by dynamic duality, the parametric identification leads to an LPV plant that is well capable of matching its computed responses with the measured speed and current in both transient and steady states.Therein, the nonlinear dynamics in nature is represented as a second-order LPV plant served for feedback control.
The identified LPV plant is then combined into performance requirements upon energy consumption and motion transience per distance, rather than per time, to form a generalized plant served for multi-objective LPV 2 L -gain control design.It is worth noting that the system is always under tracking of throttle commands and disturbances of gravity, so the energy-motion regulation is not specified by the mixed H objective [31][32][33][34], which, in the last decade, has been more popular than double 2 L -gain objective.Here the double-objective LPV 2 L -gain feedback design is improved from the previous game-theoretic control in [23], such that the embedded observer is free of structure conservatism to yield much more excellent performance on energy-motion regulation.Therein, a feasible set of Luenberger gains is allowed, which are formulated into differential linear matrix inequalities (LMIs) by a slack variable, rather than being assumed with the structure of Kalman-filtering.Moreover, a finite-element method to numerically solve the differential LMIs arising from general rate-bounded LPV control is developed.This finite-element method removes the assumption that the system matrices be affine-dependent on the time-varying parameter for gain scheduling, which is usually made in the polygon or other convex-hull methods [33,34].For the time being, such a finite-element method is the most general and non-conservative solver for differential LMIs.
Even though the sophisticated LMIs algebra was widely applied for feedback design for single 2 L -gain objectives or mixed H objectives [31][32][33][34], we did not apply them to multi-objective 2 L -gain control, for the following reasons.Virtually LMI-based control starts at bounded real-lemmas (BRL), and then arrives at convex formulations of all feasible but unstructured controllers through elimination of decision matrices, e.g., in [35,36], or addition of slack variables, e.g., in [37,38].
The main difficulty of extending these BRL-stemmed strategies into multiple 2 L -gain objectives is that the set of feasible closed-loop Lyapunov matrices is unable to numerically track back to a common output-feedback that achieves all 2 L -gain objectives.Even conservatively assuming an identical Lyapunov matrix for all objectives, this difficulty still exists, although there is no such difficulty in state feedback [39][40][41].For the present, the game-theoretic feedback control in the paper provides a solution to this difficulty.
Here the multiple 2 L -gain control is going on a Nash game, wherein the players consist of several talkers, a worker, a disturber, a negotiator and an observer.The talkers have different 2 L -gain demands, negotiating with one another through an agent.In a real time fashion, the worker delivers the control action to minimize the control storage plus estimation storage, while the disturber exerts plant disturbances, sensing noises and reference commands to maximize the storage.These two players against each other will make a decision at the same Nash equilibrium for all of the talkers' demands, rather than bringing conservatism with decisions.An observer with the maximum disturber in mind comes into play for assisting the worker in reconstructing internal status of dynamics to make optimal decisions.Including this introduction section, this paper is organized into eight sections.Section 2 details the improved version of multi-objective LPV 2 L -gain control and Section 3 provides the finite-element solution of differential LMIs therefrom.Section 4 presents the dynamic modeling and identification of city cars.Section 5 applies the multi-objective control to online manage the energy-motion of city cars.Section 6 demonstrates the procedure of implementing the energy-motion management.Section 7 starts pilot runs in full spectrum and discusses the results therefrom.Section 8 recapitulates the present work.

Multi-Objective Linear Parameter-Varying (LPV) L 2 -Gain Control
For the following pattern of generalized plants: the feedback controllers are to be designed from the Nash-game perspective to achieve multiple 2 L -gain objectives: under zero initial conditions.The matrices ) , ( 1 D B are firstly chosen as orthogonal to each other, i.e., , for canonically quadratic regulation. Therein the measured output y is sent into the to-be-designed feedback controller that real-time regulates the control action u for the multi-objective performance of Equation ( 2) with the metric of 2 L -gains from the exogenous disturbance w to its counterparts i z 's.The exogenous disturbance w can be stacked up by plant disturbances, sensor noises, reference commands, and uncertainty-induced disturbances.The performance variables i z 's are assigned for quadratic regulation of the state vector x and the control vector u, an example of which is ( , , ,ρ) where ( ,ρ) i V x is a positive-definite quadratic storage of the current state, i.e., ( ,ρ) (ρ) Explicitly, the ith Hamiltonian Hi is: where . Similar notations will appear in the sequel.
By "completing the square", the ith Hamiltonian in Equation (5) appears to be: Therefore, the minimum control * i u and the maximum disturbance * i w of the ith Hamiltonian Hi are, respectively: Substitution of Equation (7) into Equation ( 5) yields: * * ( , , ,ρ) for ρ   , then the closed-loop is asymptotically stable and, with the definition in Equation (3): when * i u u  .This is just the ith objective of Equation ( 2).Equation ( 7) reveals that the Nash equilibrium will be achieved if all the Hamiltonians i H 's share the same control storage, i.e., That is, the control and the disturbance: are the minimum control and the minimum disturbance, respectively, of all Hamiltonians simultaneously.Based on Equation ( 9), the multiple 2 L -gain objectives are all achievable if the following differential Riccati inequalities are feasible for 0  X : Now let us continue to the output feedback embedded with a calibrated Luenberger observer: where L is the Luenberger gain and Wcal is the calibration following the maximum disturbance * w in Equation ( 11) and the plant in Equation ( 1), that is: Define the error state by . Subtracting Equation (1) from Equation ( 13) in conjunction with Equation ( 14) yields: wherein all feasible Luenberger gains L are to be determined.Inferred from the state-feedback above, the ith Hamiltonian function Hi assigned to the ith 2 L -gain objective for the output feedback is to be: , , ,ρ) where the ith storage function i V is chosen as control storage plus estimation storage, i.e., ( , ,ρ) (ρ) (ρ) With Equations ( 1) and ( 15), the ith Hamiltonian function Hi becomes: in which: Known from the differential geometry, the minimum control * i u and the maximum disturbance * i w of the ith Hamiltonian Hi can be obtained through: With Equations ( 19)-( 21), it turns out that: If we let Z Z i  for i  , Nash equilibrium is achieved when all Hamiltonians share the same control storage and estimation storage.Thus, in this multi-objective Hamiltonian game, the maximum disturbance is: However, the unavailability of the full information prevents the control u from being chosen as the minimum control, . Instead, we choose: as the best estimation of the minimum control.Therein the temporal trajectory of maximum disturbance * i w is evaluated by the maximum state-feedback in Equation ( 23), which is unable to destabilize the system.On ground of the principle of Equivalent Value in optimal control, the optimal decision to an open loop is identical at any time to the control input in optimal feedback, provided that the open-loop is asymptotically stable.Here the minimum feedback of Equation ( 24) stabilized the closed-loop system.Substituting Equations ( 23) and ( 24) into Equation (18) yields:  27) to be: In summary, the multi-objective 2 L -gain performance is achievable if the following LMIs-like inequalities are feasible for where cal and Y is a slack variable.Any quadruple ( , , , α ) i X Z Y satisfying Equations ( 30) and ( 31) defines a feasible feedback.The Luenberger gain is then calculated by , and the feedback controllers take the form: Extended from the canonical LQ optimal control, this game-theoretic strategy provides a solution to multi-objective 2 L -gain robust control without bringing conservatism into the feasibility.

Finite-Element Solution of the Differential Linear Matrix Inequalities (LMIs)
Partition the parameter space centered at the grid points: . The same notation is also used for the decision function  (33) at which the values of the decision functions are to be solved.If the value of the decision function Z at the point , then the inner-product β Z  at that point can be approximated by: Similar notations become clear for 2  p .
( For brevity of notations, let the differential LMIs of Equations ( 30) and (31) be written, respectively, by: With the above notations, the finite-element version of the differential LMIs of Equations ( 30) and (31) become, respectively: (β ( ), ( ),ρ( )) 0 which spawn qNn and qN LMIs, respectively, where q is the number of vertices of the polygon  that bounds the parameter variation-rates.For a feasible set of 2 L -gain indexes i  's, firstly call 'feasp' in the Matlab-LMI toolbox to solve the feasible solution of Equation (37), which outputs the minimum solution min X over ρ   , that is 0 , ρ   .Secondly, substitute the solution min X into Equation (38), and then call the 'mincx' routine to obtain, if feasible, the solution max Z with minimum integration of trace over ρ   : , which maximizes the Luenberger gain L .This strategy is named here "min-action-max-estimation".

Modeling and Identification of Car Dynamics
The city car with speed under 55 km/h appears as the one in Figure 2.With (A1)-(A8) and electromechanical principles, the city-car dynamics in Figure 2 can be parameterized to second order with electric inductance L, electric resistance R, equivalent inertia of rotation J, mechanical damping B, gear ratio ρ, and torque coefficient  .The input includes the controlled voltage u, being the electric effort across the motor, and the equivalent torque f, being the mechanical effort from the gravity of uphill.The state of dynamics consists of the motor current i, being the fast-time state variable, and the tire speed ω, being the slow-time state.With (A5), Lagrange dynamics formulates the equivalent inertia J and the equivalent torque f to be: where m stands for the total mass of the car plus its driver, t J for the rotation inertia of each tire at the wheel hub (assuming four tires to be identical), t r for the radius of each tire, and φ for the inclination angle of the road surface from the horizon.Parametric identification is further performed in the sequel.Based on the postulation of dynamic duality, we found on experimental account the following rules: (DD-1) Like the mechanical inertia, the electric inertia is a constant.(DD-2) The electric resistance is monotonically increased with the tire speed, while the mechanical damping is monotonically decreased with the speed.(DD-3) The torque coefficient (counter electromagnetic coefficient) is a constant.
With (DD-1), the electric inductance L is taken as a constant.Without the mechanical load, the response of motor current in the neighborhood of initial time is generated by the transfer function , where 0 R is the static resistance.Give the motor a step voltage to measure the overshoot and time constant of rising current, which identify the electric inductance L and the steady-state resistance 0 R , respectively.
Then, vary the step voltages to measure rotation speeds of the motor in steady state, and record the result.It shows a linear relationship between the motor voltage and the motor speed in steady state, when the steady-state current is close to zero.Therefore, the torque coefficient  is a constant, consistent with (DD-3).With (DD-2), Kirchhoff's theorems upon the equivalent circuit in Figure 2 realizes the state space of the car dynamics to be: With this LPV parameterization, we are ready for identification of the speed-dependent electric resistance R and the mechanical damping B .Let the car on the road be maneuvered by the driving of pulse-width modulated voltages u and the braking of road-surface inclinations f , both of which are measured and recorded.A different ) , ( f u makes a different ( ,ω) i , the current and speed in steady state, and there is one-to-one mapping between them.The result is found to be in line with (DD-2), as shown in Figures 3 and 4, respectively.Since the mechanical friction in low-speed region is close to the static friction, the mechanical damping ought to be nearly inversely proportional to the tire speed with an offset representing the dynamic damping in high-speed region.The tire speed represents the mechanical flow but the counter effort in the electric part, so that the electric resistance ought to be nearly linearly proportional to the speed, with an offset representing the static resistance.Experiments joining the computer simulations show that such LPV parameterization has its computed responses in perfect agreement with those measured on the road, whereof a represented result is shown in Figures 5 and 6.For comparison, the responses of the linear LTI model with steady-state resistance 0 R and damping 0 B are also computed.Unlike LPV parameterization, LTI parameterization is unable to capture overshoots and rising times of the motor currents in practice.With road tests on a set of road-surface inclinations and driving-voltage curves, it is found that such a LPV modeling really catches the car dynamics.This indirectly verifies the principle of dynamic duality in LPV parameterization for city cars, as proposed in (DD-1), (DD-2) and (DD-3).The state variables ω and i, the time t, the electric effort u, and the mechanical effort f in Equation ( 40) can be further made dimensionless by: where u stands for the battery voltage (assuming that buck chopper is used), ω and i for the level speed and current, respectively, in steady state, as the motor voltage is kept at u .The m τ stands for the time constant of mechanical parts ) 1 ( / B J .Moreover, dimensionless resistance λ and dimensionless damping b are defined, respectively, by: The order-reduced dynamics is thus LPV-parameterized by: where β and φ indicate mechatronic ratios of AC characteristics and power ratings of electric vehicles, respectively.In the sequel, this identified LPV modeling is combined with energy and motion objectives to form a generalized plant, served for feedback design.

Energy-Motion Management
Referring to the left part of Figure 2, the following symbols are employed to construct a generalized plant.There are four state variables: The motor current i related to the rate of energy dissipated from the electric resistor.(ii) 2 x : The tire speed ω related to the rate of energy dissipated from the mechanical damper.(iii) 3 x : The integration of tracking error of motor current Also defined are the general disturbance and two performance counterparts as well as two performance indexes: The general disturbance consisting of the command reference r , the inclined gravity f , and the measurement contamination θ that is the temporal integration of current-sensor noise after white-noise filtering, as shown in Figure 2. (vi) : The performance variable related to quickness-in-motion (QM).
, The performance variable related to economy-in-energy (EE).(viii) 1 γ 0  is the performance index of QM, when the closed-loop 2 L -gain from the general disturbance q to the QM variable 1 z is set to one: where the 2 L -norms of both sides are of integrations by distance d ω d s t  so that 2 L -gain performance is signified by per-distance rather than by per-time.
(ix) 2 γ 0  is the performance index of EE, when the closed-loop 2 L -gain from the general disturbance q to the EE variable 2 z is also set to one: where the 2 L -norms are weighted by the tire speed ω to emphasize energy consumption per distance.
A smaller QM-index γ1 implies that the distance for reaching the reference speed or torque is shorter, that is, the motion is quicker under servocontrol.A smaller EE-index γ2 means that heat generation from mechanical damper and electric resistor is less in transience per distance.The QM-variable 1 z and EE-variable 2 z include the motor voltage u as an entry to protect both objectives from being achieved by cheap control, whereby to keep the motor voltage within a practical range in real operations.
As the 2 L -gains of Equations ( 44) and ( 45) are per-distance specified, we replace the performance variables, the general disturbance, the motor voltage, and the output (the input to the controller as shown in Figure 2) by: respectively, in order to construct the generalized plant served for the control design in Section 2.
Combination of performance requirements in Equations ( 44) and ( 45) and the nominal plant in Equation ( 43) results in the following LPV generalized plant: where small positive values 1 ε and 2 ε , which indicate relaxed integrations, are chosen to replace zeros to avoid the singularity of 2 L -gain control regarding to DC disturbances.Moreover, a freely chosen parameter 3 ε is also included for additional weighing on torque tracking.
Equations ( 47)-( 50) can be briefly written as: where the tire speed ω [ε,1]  plays as the slow-time gain-scheduling parameter, other than a state variable, in LPV control.In case of ε 0  , the state-space realization of Equation (51) will be uncontrollable and unobservable.In practice, the positive number ε can be set arbitrarily close to zero.The rate-bounds of the tire speed are assumed available at the stage of propulsion design: The LPV generalized plant of Equation ( 51) is served for multi-objective LPV 2 L -gain control design presented in Section 2.
The feedback design in Section 2 in conjunction with its finite-element solver in Section 3 is then employed to solve feasible quadruples 1 2 ( , ,γ ,γ ) X Z that guarantee the performances QM-index 1 γ and EE-index 2 γ .Here, we provide a numerical procedure as follows to figure out the set of feasible region (γ ,γ ) in the 1 γ -2 γ 2D plane.Choose a set of negotiation parameters α and repeat the procedure to obtain the corresponding points.Connect those points into a curve, the left-down region of which must be of infeasibility.
Figure 7 demonstrates in non-dimensional fashion the nature of city cars: the trade-off between energy economics and quick motions.This vehicular nature has been commonly experienced, and is parameterized here.In the Figure 7, the feedbacks with calibrated non-conservative observers embedded are juxtaposed with those with calibrated Kalman-filters embedded to visualize the claimed non-conservatism.

Implementation of Energy-Motion Management
The procedure to implement the controller that fulfills the online energy-motion management is described as follows: Step 1: Identify equivalent parameters of city car.Identify, as in Section 4, the equivalent inductance L, resistance R and torque coefficient κ of electric parts as well as the equivalent damping B and inertia of moment J of mechanical parts.
Step 2: Choose the mode of transient transmission.
Choose a proper mode of transient transmission with negotiation parameter α in the EE-QM regulator of Figure 7, and then calculate the control dynamics of EE-QM regulation with the corresponding pair of performance indexes 2 1 (γ ,γ ) .Plug the control dynamics into the dynamics in Equation ( 43), and then perform computer simulations with the sampling time T .From the computed responses prepare a set of negotiation parameters for desired transmission modes.
Step 3: Implement the transmission mode into a microcontroller.The energy-motion regulator is then implemented into a Microchip-dsPIC chip with the Euler discretization as follows.Continuous-to-digital conversion yields: where k and j are fast-time and slow-time sampling indexes, respectively.The system matrices ) , (   can be calculated online in a slow-time fashion; explicitly: where T is the fast-time sampling period that is short enough to make the above Taylor-series expansion accurate even with few terms, say, 2  N or 3 as T is 0.001 s.In the program of Equations ( 53) and (54), ( A , B , C ) are the system matrices of the integrator in series with the multi-objective LPV 2 L -gain feedback as shown in Figure 2. At the present time k t , the dsPIC merely stores the current state k x , and the update of state from k x to 1  k x at next instant is fulfilled by the DSP-engine performing addition and multiplication of floating numbers according to Equation (53).In this sense, any time can be treated as the initial time, which makes the real-time processing as efficient as possible.Moreover, the intervals of state update are held identical to those in computer simulation, so that the real-timed operation matches the dynamics that has been verified by offline calculation, thus achieving robust implementation.At any instance, a pulse-width modulated (PWM) signal in line with the control signal u is sent to the gate-driving circuit of the switched power converter, as shown in Figure 8. Step 4: Prepare controller circuit board.Figure 9 shows main components of the controller in design.On the circuit board are (1) a current sensor; (2) a dsPIC chip embedded with A/Ds, PWMs and the DSP programs; (3) voltage regulators; (4) switching buck choppers; (5) sensor and transducer circuits; (6) the amplification circuits driving the gates of power MOSFETs in the choppers; and (7) starting and braking relays.Figure 10 plots the auxiliary circuit of the dsPIC chip.Its reliability is known from considerable times of trials-and-errors.

Experimental Results and Discussion
A city car with 96  u V, 52  i A and ω 9200  rpm, installing an ACS754SCB-200 Hall-effect current sensor/transducer (single bias; 0-100A input; 0-5V output) is made for these experiments.A toy DC-motor is pivoted on the axial of the motor for sensing tire speeds, and a USB4711 box is for real-timed signal acquisition of motor currents and tire speeds.The signal conditioning inside the dsPIC chip is as follows: (1) the command signal 0-5 V corresponds to the motor current 0-100 A; (2) the output of current sensor 0 V-5 V corresponds to the motor current 0-100 A; (3) the frequency of the PWM is 20 k Hz with its duty cycle 0%-100% being proportional to motor average voltage 0-96 V; and (4) the sampling time in the MIT program is 001 .0  T s.During the pilot run, we embed two modes of transient transmission into the dsPIC microcontroller to verify the multi-objective control strategy developed in Section 5.One is an economic type with the negotiation parameter α 1.1  , and the other is a quick type with the negotiation parameter being α 2.0  .Figure 11 records the typical response of these two types of controllers.We find that Allegro ACS754SCB-200 Hall-effect current sensor is of wide bandwidth and thus lures high-frequency noises into the loop.Though noises contaminate measurement, we can discern that the multi-objective strategy works.It is more convenient to observe the responses of energy consumption per distance, especially in traffic jams, by changing the independent variable from the time t to the angular displacement θ ( dθ / d ω t  ).Accordingly, we record the responses of energy consumption, motor current and tire speed versus the angular displacement θ of the city car driven in a traffic jam in a short range and long range, respectively.Their responses of energy consumption indicate that the economic type ( α 1.1  ) saves about 50% energy with respect to the quick type ( α 2.0  ).Energy is saved because of reduced electric dissipation resulted from the depression of current overshoots during the acceleration and deceleration.This also protects the battery from overcharging and the DC motor from overheating, and thus saves energy.Letting the car run in the local city and accumulating data therefrom, we calculate the statistical relationship between the amount of energy being saved and the negotiation parameter, which is plotted in Figure 12.It is found that the percentage of energy saved with respect to the quick type ( α 2.0  ), η, is almost inversely proportional to the negotiation parameter.These findings correspond to our experiences: commands of rapid acceleration and deceleration tend to shorten driving ranges; however, it is usually necessary for rapid accelerations and decelerations to safely drive the car in a crowded condition.Therefore, every condition of traffic has its own mode of transient transmission when manufacturers employ our EE-QM regulator in designing controllers of city cars.
Furthermore, we record the responses of motor current and car speed when driving on level ground and uphill, respectively.On the level ground, the economic type significantly suppresses the current overshoots in transience with a bit of retard of speed.This significantly saves energy in traffic jams, when electric energy is unreasonably dissipated due to current overshoots happening during acceleration and deceleration.However, on the uphill the retardation of speed becomes significant.Then, the car should switch to the controller with a negotiation parameter large enough to climb up a steep hill.When the car is back to the level road, it should switch back to economic type to reduce energy consumption.Therefore, the negotiator parameter in transience compares to the gear ratio in steady state for regulating the EE-QM behavior, and the regulation is well captured by Figure 7, obtained by the multi-objective LPV 2 L -gain control.
Next, we design a special road test, arranging a 1000-meter runway equipped with several stop flags.The experimental car stops at each flag, and then starts with maximum throttle up to the level speed at 40 km/h.Then, the car decelerates and stops at the next flag.We recorded the responses of motor current, tire speed, battery voltage, and PWM duty cycles to estimate the power dissipation.It is proven that per-distance control strategy has significantly better EE-QM performance than per-time control strategy, especially in a traffic jam.Experimental results are recorded in Figures 13-15.
It is seen in Figure 13 that the overshoots of motor current are significantly suppressed, especially when the car is accelerated at the times of stop, say at 0 s, 36 s, 72 s, etc.That is, the per-distance control strategy not only saves energy but also protects the battery from being overcharged and power converters from being burned out.As for the small current peaks at other times, they result from small throttles needed to control the car to stop at the flags, at the times the car is running at a low speed.These peaks, large or small, have the values in line with the LPV nature of the city car, as shown in Figure 3: the equivalent resistance is almost proportional to the car speed.We see in Figure 14 that the significant suppression of current overshoots is accompanied by a slight retardation of speed.That is, energy efficiency is achieved in transience nearly without sacrificing the time efficiency.In Figure 14, we calculate the consumption of electric energy versus the driving range, which shows that the LPV per-distance control saves a surprising amount of energy.More energy will be saved in heavier traffic.

Conclusions
This research work presents four prongs of knowledge regarding online energy management of city cars: 1) The current multi-objective LPV 2 L -gain control significantly improves the performance of energy-motion regulation.2) Per-distance control strategy is more suitable than per-time strategy for energy management of city cars.It spawns differential LMIs, which can be numerically solved with the finite-element approximation.
3) The LPV parameterization can perfectly capture the nonlinear nature of DC propulsion, which provides system-level understanding to the common experience that quicker motions always give rise to larger consumption of energy.4) To arrive at sound energy-motion Pareto regulation, the iteration of feedback design and dynamic identification is needed.

Figure 2 .
Figure 2. Modeling and controller design of DC-motorized city cars.

Figure 3 .
Figure 3. Identification of dependence electric resistance on tire speed.

Figure 4 .
Figure 4. Identification of dependence of mechanical damping on tire speed.
command.The former two state variables, x1 and x2, quantify the rate of energy consumption, and the latter two state variables, x3 and x4 indicate the responses of motion.

Figure 7 .
Figure 7. Improved trade-off between energy and motion objectives.

Figure 8 .
Figure 8. Means plus functions of the controller and instrumentation.

Figure 9 .
Figure 9. Sensor and driver circuits in the controller.

Figure 11 .
Figure 11.Measured responses of economic type and quick type.

Figure 13 .
Figure 13.Current responses in LTI and LPV controls.

Figure 14 .
Figure 14.Motion responses in LTI and LPV controls.
Therein the matrices Qi's are positive-definite, and the performance indexes αi's are positive scalars that quantify the negotiation among those individual 2