PATH TRACKING CONTROL METHOD AND DEVICE FOR INTELLIGENT ELECTRIC VEHICLES

Abstract

A path tracking control method for an intelligent electric vehicle. The method includes the following steps. A lateral stability state of a vehicle is determined, where the lateral stability state includes a stable state, a critical-destabilized state, and a destabilized state. Path tracking control is performed on the vehicle according to the lateral stability state of the vehicle. This application also provides a path tracking control device for an intelligent electric vehicle, which includes a lateral stability state determination module and a path tracking control module.

Claims

1. A path tracking control method for an intelligent electric vehicle, comprising: determining a lateral stability state of a vehicle, wherein the lateral stability state comprises a stable state, a critical-destabilized state, and a destabilized state; performing path tracking control on the vehicle according to the lateral stability state of the vehicle; determining, when the vehicle is in the stable state, a front wheel steering increment of the vehicle; and determining a front wheel steering of the vehicle based on the front wheel steering increment of the vehicle; determining, when the vehicle is in the critical-destabilized state, a first yaw moment increment and an additional yaw moment weighting coefficient of the vehicle; determining a first additional yaw moment based on the first yaw moment increment and the additional yaw moment weighting coefficient; and determining a first longitudinal force on each wheel of the vehicle based on the first additional yaw moment; and determining, when the vehicle is in the destabilized state, a second yaw moment increment of the vehicle; determining a second additional yaw moment based on the second yaw moment increment; and determining a second longitudinal force on each wheel based on the second additional yaw moment; wherein the stable state is determined by $0<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.??;$ the critical-destabilized state is determined by $?<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.?1;$ and the destabilized state is determined by $.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.>1;$ wherein ? represents a center-of-mass sideslip angle of the vehicle; {dot over (?)} represents a center-of-mass sideslip angular velocity; k represents a slope; c is constant; and A represents a road surface adhesion coefficient.

2. The path tracking control method of claim 1, wherein when the vehicle is in the stable state, the front wheel steering increment of the vehicle is determined through steps of: constructing an objective function J: $\min J={.Math.}_{j=i+1}^{N_p}?\left(t+j.Math.t\right)-{?}_{\mathrm{ref}}\left(t+j.Math.t\right)Q^2+{.Math.}_{i=0}^{N_c-1}?u\left(t+i.Math.t\right)R^2+1{?}^2;$ wherein N.sub.p represents a prediction time domain; N.sub.c represents a control N time domain; t represents a current moment; ?(t+j|t) represents a lateral position coordinate of the vehicle at a (t+j).sup.th moment when the current moment is t; ?.sub.ref(t+j|t) represents a reference value of the lateral position coordinate at the (t+j).sup.th moment when the current moment is t; Q, R, and l represent weighting matrices of an output, a control increment, and a relaxation factor, respectively; ? represents the relaxation factor; and ?u(t+i|t) represents a front wheel steering increment of the vehicle at a (t+i).sup.th moment when the current moment is t; solving the objective function J; and determining a front wheel steering increment ?u(t+i|t) of the vehicle at each moment within the control time domain when the objective function J reaches a minimum value; and acquiring a front wheel steering u(t) of the vehicle at the current moment t based on a front wheel steering increment ?u(t+0|t) at a (t+0).sup.th moment within the control time domain: $u\left(t\right)=u\left(t-1\right)+?u\left(t+0.Math.t\right)$

3. The path tracking control method of claim 1, wherein when the vehicle is in the critical-destabilized state, the first yaw moment increment is determined through steps of: constructing an objective function J: $\min J={.Math.Y\left(k\right)-Y_{\mathrm{ref}}\left(k\right).Math.}_Q^2+{.Math.?U\left(k\right).Math.}_R^2={.Math.}_{j=i+1}^{N_p}{.Math.Y\left(k+j.Math.k\right)-Y_{\mathrm{ref}}\left(k+j.Math.k\right).Math.}_Q^2+{.Math.}_{i=0}^{N_c-1}{.Math.?u\left(k+i.Math.k\right).Math.}_R^2;$ wherein N.sub.p represents a prediction time domain; N.sub.c represents a control time domain; Y(k) represents a predicted value of an output of the vehicle in the prediction time domain N.sub.p; k represents a current moment; Y.sub.ref ( ) represents a reference value of the output; ?U(k) represents a sequence of a yaw moment increment at a current moment k; Y(k+j|k) represents a predicted value of the output at a (k+j).sup.th moment when the current moment is k; Y.sub.ref(k+j|k) represents a reference value of the output at the (k+j).sup.th moment when the current moment is k; ?u(k+i|k) represents a yaw moment increment at a (k+i).sup.th moment when the current moment is k; and Q and R represent weight matrices of the output and a control increment, respectively; solving the objective function J; and determining the yaw moment increment ?u(k+i|k) at each moment within the control time domain when the objective function J takes a minimum value.

4. The path tracking control method of claim 1, wherein when the vehicle is in the critical-destabilized state, the additional yaw moment weighting coefficient is determined through steps of: constructing a coordination controller based on a reinforcement learning deep deterministic policy gradient (DDPG) algorithm, wherein state inputs are yaw velocity w.sub.r, center-of-mass sideslip angle ?, lateral displacement Y, heading angle ?, front wheel steering angle ?.sub.f and lateral acceleration a.sub.y; and an action output is the additional yaw moment weighting coefficient ?; wherein a reward function of the coordination controller is as follows: $r=r_1+r_2+r_3;$ $\mathrm{wherein}{r}_1=-\left(e_y^2+e_{?}^2\right)+C_1;r_2=-.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.+C_2;r_3=-.Math.a_y.Math.+C_3;$ e.sub.y represents a deviation between an actual lateral displacement of the vehicle and a reference path; e.sub.? represents a deviation between the heading angle and a reference heading angle; a.sub.y represents a lateral acceleration of the vehicle; C.sub.1, C.sub.2, and C.sub.3 are constant; c represents a coefficient; k represents a slope; {dot over (?)} represents a center-of-mass sideslip angular velocity of the vehicle; and ? represents a center-of-mass sideslip angle of the vehicle.

5. The path tracking control method of claim 1, wherein the path tracking control is performed through steps of: acquiring a yaw moment u(k) at a current moment k based on a yaw moment increment ?u(k+0|k) at a (k+0).sup.th moment within a control time domain; wherein u(k)=u(k?1)+?u(k+0|k); determining an additional yaw moment ?M according to the yaw moment u(k) and the additional yaw moment weighting coefficient ?; wherein ?M=?u(k); and determining a longitudinal force on each wheel of the vehicle based on the additional yaw moment ?M in a drive-brake-combined manner.

6. The path tracking control method of claim 1, wherein when the vehicle is in the destabilized state, the second yaw moment increment of the vehicle is determined, and the path tracking control is performed based on the second yaw moment increment through steps of: constructing an objective function J: $\min J={.Math.Y\left(k\right)-Y_{\mathrm{ref}}\left(k\right).Math.}_Q^2+{.Math.?U\left(k\right).Math.}_R^2={.Math.}_{j=i+1}^{N_p}{.Math.Y\left(k+j.Math.k\right)-Y_{\mathrm{ref}}\left(k+j.Math.k\right).Math.}_Q^2+{.Math.}_{i=0}^{N_c-1}{.Math.?u\left(k+i.Math.k\right).Math.}_R^2;$ wherein N.sub.p represents a prediction time domain; N.sub.c represents a control time domain; Y(k) represents a predicted value of an output of the vehicle in the prediction time domain N.sub.p; k represents a current moment; Y.sub.ref ( ) represents a reference value of the output; ?U(k) represents a sequence of a yaw moment increment at a current moment k; Y(k+j|k) represents a predicted value of the output at a (k+j).sup.th moment when the current moment is k; Y.sub.ref(k+j|k) represents a reference value of the output at the (k+j).sup.th moment when the current moment is k; ?u(k+i|k) represents a yaw moment increment at a (k+i).sup.th moment when the current moment is k; and Q and R represent weight matrices of the output and a control increment, respectively; solving the objective function J; and determining the yaw moment increment ?u(k+i|k) at each moment within the control time domain when the objective function J takes a minimum value; acquiring a yaw moment u(k) at a current moment k based on a yaw moment increment ?u(k+0|k) at a (k+0).sup.th moment within a control time domain; wherein u(k)=u(k?1)+?u(k+0|k); determining an additional yaw moment ?M according to the yaw moment u(k) and the additional yaw moment weighting coefficient ?; wherein ?M=?u(k); and determining a longitudinal force on each wheel of the vehicle based on the additional yaw moment ?M in a drive-brake-combined manner.

7. A path tracking control device for an intelligent electric vehicle, comprising: a lateral stability state determination module; and a path tracking control module; wherein the lateral stability state determination module is configured to determine a lateral stability state of a vehicle, wherein the lateral stability state comprises a stable state, a critical-destabilized state, and a destabilized state; and the path tracking control module is configured to perform path tracking control on the vehicle according to the lateral stability state of the vehicle, determine a front wheel steering increment of the vehicle when the vehicle is in the stable state, determine a front wheel steering of the vehicle based on the front wheel steering increment of the vehicle, determine a first yaw moment increment and an additional yaw moment weighting coefficient of the vehicle when the vehicle is in the critical-destabilized state, determine a first additional yaw moment based on the first yaw moment increment and the additional yaw moment weighting coefficient, determine a first longitudinal force on each wheel of the vehicle based on the first additional yaw moment, determine a second yaw moment increment of the vehicle when the vehicle is in the destabilized state, determine a second additional yaw moment based on the second yaw moment increment, and determine a second longitudinal force on each wheel based on the second additional yaw moment; wherein the stable state is determined by $0<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.??;$ the critical-destabilized state is determined by $?<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.?1;$ and the destabilized state is determined by $.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.>1;$ wherein ? represents a center-of-mass sideslip angle of the vehicle; {dot over (?)} represents a center-of-mass sideslip angular velocity; k represents a slope; c is constant; and ? represents a road surface adhesion coefficient.

8. A computer-readable storage medium, wherein the computer-readable storage medium is configured to store a computer program; and the computer program is configured to be executed by a processor to implement the path tracking control method of claim 1.

9. A computer program product, comprising: a computer program or a computer instruction; wherein the computer program or the computer instruction is configured to be executed by a processor to implement the path tracking control method of claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0052] The present application may be better understood with reference to the following description in conjunction with the accompanying drawings, which are included in and form a part of this specification.

[0053] FIG. 1 is a flow chart of a path tracking control method for an intelligent electric vehicle according to an embodiment of the present disclosure; and

[0054] FIG. 2 is a block diagram of a structure of a path tracking control device for an intelligent electric vehicle according to an embodiment of the present disclosure.

DETAILED DESCRIPTION OF EMBODIMENTS

[0055] Exemplary embodiments of the present application will be described below with reference to the accompanying drawings. For the sake of clarity and brevity, not all features of the actual embodiments are described in the specification. However, it should be understood that many decisions specific to embodiment can be made during the development of any such practical embodiments to achieve the specific goals of the developer, and these decisions may vary from embodiment to embodiment.

[0056] Furthermore, to avoid obscuring the present application with unnecessary details, only the structure closely related to the embodiments according to the present application is shown in the accompanying drawings, and other details less relevant to the present application have been omitted.

[0057] It should be understood that the present application is not limited to the form of embodiment described with reference to the accompanying drawings below. As far as practicable herein, embodiments may be combined with each other, features can be substituted or borrowed between different embodiments, and one or more features can be omitted in one embodiment.

[0058] FIG. 1 is a flow chart of a path tracking control method for an intelligent electric vehicle. Referring to FIG. 1, the path tracking control method includes the following steps.

[0059] (S110) A lateral stability state of a vehicle is determined, where the lateral stability condition monitoring function of the vehicle is turned on, and the lateral stability state includes a stable state, a critical-destabilized state, and a destabilized state.

[0060] In this step, the lateral stability state of the vehicle is determined through a center-of-mass sideslip angle-sideslip angular velocity phase plane method. Specifically, the stable state is determined by

[00009]$0<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.??;$

the critical-destabilized state is determined by

[00010]$?<.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.?1;$

and the destabilized state is determined by

[00011]$.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.>1;$ [0061] where ? represents a center-of-mass sideslip angle of the vehicle; {dot over (?)} represents a center-of-mass sideslip angular velocity; k represents a slope; c is constant; and A represents a road surface adhesion coefficient.

[0062] (S120) Path tracking control is performed on the vehicle according to the lateral stability state of the vehicle. When the vehicle is in the stable state, a front wheel steering increment of the vehicle is determined, and a front wheel steering of the vehicle is determined based on the front wheel steering increment of the vehicle. When the vehicle is in the critical-destabilized state, a first yaw moment increment and an additional yaw moment weighting coefficient of the vehicle are determined, a first additional yaw moment is determined based on the first yaw moment increment and the additional yaw moment weighting coefficient, and a first longitudinal force on each wheel of the vehicle is determined based on the first additional yaw moment. When the vehicle is in the destabilized state, a second yaw moment increment of the vehicle is determined, a second additional yaw moment is determined based on the second yaw moment increment, and a second longitudinal force on each wheel is determined based on the second additional yaw moment.

[0063] In an embodiment, the stability state of the vehicle is determined, and then according to the different stability states of the vehicle, different control strategies are adopted to maintain the path tracking accuracy while taking into account the stability control of the vehicle, which improves the safety of the intelligent electric vehicle.

[0064] In an embodiment, when the vehicle is in the stable state, the front wheel steering increment of the vehicle is determined, and the front wheel steering is determined based on the front wheel steering increment through the following steps.

[0065] Kinematic equations are established as follows, based on the vehicle kinematics model:

[00012]$\begin{array}{cc}\{\begin{array}{c}\stackrel{.}{X}=v_r\cos ?\\ \stackrel{.}{Y}=v_r\sin ?\\ \stackrel{?}{?}=\frac{v_r\tan {?}_f}{l}\end{array};& \left(1\right)\end{array}$ [0066] where X represents a lateral position coordinate of the vehicle; Y represents a longitudinal position coordinate of the vehicle; ? represents a heading angle of the vehicle; {dot over (X)}, {dot over (Y)} and {dot over (?)} are derivatives of lateral position coordinate, the longitudinal position coordinate and the heading angle, respectively; v.sub.r represents a speed of the vehicle; l represents a wheelbase of the vehicle; and ?.sub.f represents a front wheel steering angle.

[0067] The path tracking model is established by the above kinematic equations and the objective function is constructed to solve.

[0068] The kinematic equations are represented in state space:

[00013]$\begin{array}{cc}\{\begin{array}{c}\stackrel{?}{?}=f\left(?,u\right)\\ ?=C.Math.?\end{array};& \left(2\right)\end{array}$ [0069] where a state variable is represented by ?=[X, Y, ?].sup.T; {dot over (?)} represents a derivative of the state variable ?; ? denotes a lateral position coordinate of the vehicle; and the a control variable is represented by u=[?.sub.f].sup.T; ?.sub.f represents a front wheel steering angel; and

[00014]$C=\left[\begin{array}{ccc}1& 1& 0\\ 0& 1& 1\end{array}\right].$

[0070] The formula (2) is linearized and discretized to obtain a discrete state space equation as follows:

[00015]$\begin{array}{cc}\{\begin{array}{c}{\stackrel{?}{?}}_{k+1}={\stackrel{?}{A}}_{k,t}{\stackrel{?}{?}}_k+{\stackrel{?}{B}}_{k,t}{\tilde{u}}_k\\ ?=C.Math.{\stackrel{?}{?}}_k\end{array};& \left(3\right)\end{array}$

[0071] where ?.sub.k,t represents a matrix of difference coefficients of the state variable at a k.sup.th moment when the current moment is t; B.sub.k,t represents a matrix of difference coefficients of the control variable at the k.sup.th moment when the current moment is t;

[00016]${\stackrel{?}{A}}_{k,t}=\left[\begin{array}{ccc}1& 0& -T{v}_r\sin ?\\ 0& 1& T{v}_r\sin ?\\ 0& 0& 1\end{array}\right];{\stackrel{?}{B}}_{k,t}=\left[\begin{array}{cc}T\cos ?& 0\\ T\sin ?& 0\\ \frac{T\tan {?}_f}{l}& \frac{T{v}_r}{l{\cos }^2{?}_f}\end{array}\right];$

T represents a sampling period; {tilde over (?)}represents a state variable difference, and {tilde over (?)}=???.sub.r; ? represents a control variable difference, and ?=u?u.sub.r; {dot over (?)}.sub.r=f(?.sub.r, u.sub.r), which is a Taylor expansion point at any point; ?.sub.r represents a state variable at any point; u.sub.r represents a control variable at any point; &k represents a state variable difference at the k.sup.th moment; and ?.sub.k represents a control variable difference at the k.sup.th moment.

[0072] The state variable and the control variable are combined to reconstruct a new state space expression ?(k|t) as follows:

[00017]$\begin{array}{cc}?\left(k.Math.t\right)={\left[\stackrel{?}{?}\left(k.Math.t\right)u\left(k-1.Math.t\right)\right]}^T;& \left(4\right)\end{array}$

[0073] The equation (4) is substituted into equation (3) to obtain equation (5) as follows:

[00018]$\begin{array}{cc}\{\begin{array}{c}?\left(k+1.Math.t\right)={?}_{k,t}?\left(k.Math.t\right)+{\stackrel{?}{B}}_{k,t}?u\left(k.Math.t\right)\\ ?\left(k.Math.t\right)={\stackrel{?}{C}}_{k,t}.Math.?\left(k.Math.t\right)\end{array};& \left(5\right)\end{array}$$\mathrm{where}{?}_{k,t}=\left[\begin{array}{cc}{\stackrel{?}{A}}_{k,t}& {\stackrel{?}{B}}_{k,t}\\ O_{m?n}& I_m\end{array}\right];{\stackrel{?}{B}}_{k,t}=\left[\begin{array}{c}{\stackrel{?}{B}}_{k,t}\\ I_m\end{array}\right];$

?u(k|t) represents a control variable increment at a k.sup.th moment when the current moment is t, and ?u(k|t)=u(k|t)?u(k?1|t); ?(k|t) represents a lateral position coordinate of the vehicle at the k.sup.th moment when the current moment is t; ?R=[C Omxn]; m represents a control variable dimension, and m=1; n represents a state variable dimension, and n=3; Im represents a unit matrix of m dimensions; and O.sub.m?n represents a 0 matrix of m?n.

[0074] To simplify the calculation, the prediction time domain is allowed to be N P the control time domain is allowed to be N.sub.c, and assumptions are made as follows:

[00019]$\begin{array}{cc}\{\begin{array}{ccc}{\hat{A}}_{k,t}={\hat{A}}_{t,t},& k=1,2,.Math.,& t+N_P-1\\ {\hat{B}}_{k,t}={\hat{B}}_{t,t},& k=1,2,.Math.,& t+N_P-1\\ {\hat{C}}_{k,t}={\hat{C}}_{t,t},& k=1,2,.Math.,& t+N_P-1\end{array}.& \left(6\right)\end{array}$

[0075] An expression for a predicted output of the system in the prediction time domain N.sub.p is derived as follows:

Y(t)=??(t|t)+??U(t)(7); [0076] in the formula (7),

[00020]$Y\left(k\right)=\left[\begin{array}{c}?\left(t.Math.t\right)\\ ?\left(t+1.Math.t\right)\\ .Math.\\ ?\left(t+N_c.Math.t\right)\\ .Math.\\ ?\left(t+N_p.Math.t\right)\end{array}\right]?\left(t\right)=\left[\begin{array}{c}{\hat{C}}_{t,t}{\hat{A}}_{t,t}\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}^2\\ .Math.\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_c}\\ .Math.\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_p}\end{array}\right]?=\left[\begin{array}{cccc}{\hat{C}}_{t,t}{\hat{B}}_{t,t}& 0& .Math.& 0\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}{\hat{B}}_{t,t}& {\hat{C}}_{t,t}{\hat{B}}_{t,t}& .Math.& 0\\ .Math.& .Math.& ?& .Math.\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_c-1}{\hat{B}}_{t,t}& {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_c-2}{\hat{B}}_{t,t}& .Math.& {\hat{C}}_{t,t}{\hat{B}}_{t,t}\\ {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_c}{\hat{B}}_{t,t}& {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_c-1}{\hat{B}}_{t,t}& .Math.& {\hat{C}}_{t,t}{\hat{A}}_{t,t}{\hat{B}}_{t,t}\\ .Math.& .Math.& ?& .Math.\\ {\hat{C}}_k{\hat{A}}_K^{N_P-1}{\hat{B}}_k& {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_p-2}{\hat{B}}_{t,t}& .Math.& {\hat{C}}_{t,t}{\hat{A}}_{t,t}^{N_p-N_c-1}{\hat{B}}_{t,t}\end{array}\right].$

[0077] ?U(t) represents a sequence of control increment corresponding to the current moment t.

[0078] During the path tracking control of the vehicle, the state error and the control increment of the model prediction controller at each sampling moment are required to be minimized to obtain the best path tracking performance and driving stability of the vehicle. Therefore, the objective function J of the following form is used:

[00021]$\min J={.Math.}_{j=i+1}^{N_p}?\left(t+j.Math.t\right)-{?}_{\mathrm{ref}}\left(t+j.Math.t\right)Q^2+{.Math.}_{i=0}^{N_c-1}?u\left(t+i.Math.t\right)R^2+l{?}^2;$ [0079] where N.sub.p represents a prediction time domain; N.sub.c represents a control time domain; t represents a current moment; ?(t+j|t) represents a lateral position coordinate of the vehicle at a (t+j).sup.th moment when the current moment is t; ?.sub.ref(t+j|t) represents a reference value of the lateral position coordinate at the (t+j).sup.th moment when the current moment is t; Q, R, and l represent weighting matrices of an output, a control increment, and a relaxation factor, respectively; ? represents the relaxation factor; and ?u(t+i|t) represents a front wheel steering increment of the vehicle at a (t+i).sup.th moment when the current moment is t.

[0080] In the actual control process, the control variable and the control increment are required to be limited within a reasonable range with constraint expressions as follows:

[00022]$\begin{array}{cc}\{\begin{array}{c}{u}_{\min }\left(t\right)?u\left(t\right)?u_{\max }\left(t\right)\\ ?u_{\min }\left(t\right)??u\left(t\right)??u_{\max }\left(t\right)\end{array};& \left(9\right)\end{array}$ [0081] where u.sub.min represents the minimum value of the control variable; u.sub.max represents the maximum value of the control variable; ?u.sub.min represents the minimum value of the control increment; and ?u.sub.max represents the maximum value of the control increment.

[0082] After solving the objective function in each control time domain, a series of sequences of increment in the control time domain are obtained, i.e., the front wheel steering increment ?u(t+i|t) of the vehicle when the objective function J takes the minimum value.

[0083] According to the model predictive control (MPC) principle, only the first element in the sequence of increment is taken as the actual control input increment to act on the system each time, i.e., the front wheel steering increment ?u(t+0|t) of the vehicle is selected to act on the system to obtain the front wheel steering u(t) of the vehicle corresponding to the current moment t, so as to realize the path tracking control on the vehicle, where u(t)=u(t?1)+?u(t+0|t).

[0084] In an embodiment, when the vehicle is in the critical-destabilized state, a first yaw moment increment and an additional yaw moment weighting coefficient of the vehicle are determined, a first additional yaw moment is determined based on the first yaw moment increment and the additional yaw moment weighting coefficient, and a first longitudinal force on each wheel of the vehicle is determined based on the first additional yaw moment through the following steps.

[0085] (1) The first yaw moment increment is determined.

[0086] An additional yaw moment controller is established, and the dynamics equations considering the yaw stability control are as follows:

[00023]$\begin{array}{cc}\stackrel{?}{v}=\frac{\cos ?}{m}\left[F_{\mathrm{lf}}\cos {?}_f-F_{\mathrm{sf}}\sin {?}_f+F_{lr}\right]+\frac{\sin ?}{m}\left(F_{\mathrm{sf}}\cos {?}_f+F_{\mathrm{lf}}\sin {?}_f+F_{sr}\right)& \left(10\right)\end{array}$$\stackrel{?}{?}=\frac{\cos ?}{mv}\left(F_{\mathrm{sf}}\cos {?}_f+F_{\mathrm{lf}}\sin {?}_f+F_{sr}\right)-\frac{\sin ?}{mv}\left(F_{\mathrm{lf}}\cos {?}_f-F_{\mathrm{sf}}\sin {?}_f+F_{lr}\right)-?$$\stackrel{?}{?}=\frac{1}{I_z}\left(F_{\mathrm{lf}}a\sin {?}_f+F_{\mathrm{sf}}a\cos {?}_f-F_{\mathrm{sf}}b+?M\right);$ [0087] where ? represents a center-of-mass side slip angle of the vehicle; {dot over (?)} represents a center-of-mass side slip angular velocity of the vehicle; ? represents a yaw angular velocity; {dot over (?)} represents a derivative of the yaw angular velocity; ?M represents a yaw moment; I.sub.z represents a moment of inertia of the vehicle around the z-axis; a and b are distances from the center of mass to the front and rear axes, respectively; F.sub.sf and F.sub.sr are a lateral force of a front axle and a lateral force of the rear axle, respectively; F.sub.lf and F.sub.lr are a longitudinal force of the front axle and a longitudinal force of the rear axle, respectively; ?.sub.f is a steering angle of the front wheel; {dot over (v)} is a longitudinal acceleration of the vehicle; and m represents a weight of the vehicle.

[0088] After linearizing the above dynamics equations, the system state equation is obtained as follows:

[00024]$\begin{array}{cc}\{\begin{array}{c}\stackrel{?}{X}=AX+BU+\mathrm{ED}\\ Y=\mathrm{CX}\end{array};& \left(11\right)\end{array}$ [0089] where a state variable is represented by X=[v, ?, ?].sup.T; an output variable is represented by Y=[?, ?].sup.T; a control variable is represented by U=[?M], an output quantity matrix is expressed as

[00025]$C=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\end{array}\right];D=\left[{?}_f\right];$

? represents a yaw angular velocity; and ? represents a center-of-mass sideslip angle of the vehicle.

[0090] To simplify the calculation and the control system, the linear continuous system is discretized using the forward Euler method, expressed as:

X(k+1)=A.sub.kX(k)+B.sub.kU(k)+E.sub.kD(k)(12); [0091] where X(k) is a state variable at a k.sup.th moment; X(k+1) is a state variable at a (k+1).sup.th moment; U(k) represents a control variable at the k.sup.th moment; D(k) is a front wheel angle at the k.sup.th moment; Ex is a parameter matrix of the vehicle at the k.sup.th moment; A.sub.k=1+TA; B.sub.k=TB; T represents a sampling period; I represents a unit matrix. Since the prediction time domain is small, D(k) is considered as constant in the prediction time domain, and does not affect the construction of the objective function. Therefore, the E.sub.kD(k) term can be ignored.

[0092] The model predictive model equation is constructed let. Letting

[00026]$?\left(k\right)=\left(\begin{array}{c}x\left(k\right)\\ U\left(k-1\right)\end{array}\right)\mathrm{and}?U\left(k\right)=U\left(k\right)-U\left(k-1\right),$

the state space equation of the predictive model is obtained as follows:

[00027]$$\begin{gathered} \begin{array}{cc}?\left(k+1\right)={?}_k?\left(k\right)+{\stackrel{?}{B}}_k?U\left(k\right)+{?}_{k}D\left(k\right) \\ ?\left(k\right)=C_k?\left(k\right);& \left(13\right)\end{array} \end{gathered}$$ [0093] where n represents the number of the state variable; m represents the number of the control variable; ?(k) represents a state variable at a k.sup.th moment; and ?(k) represents a lateral position coordinate of the vehicle at a k.sup.th moment; [0094] where

[00028]${\stackrel{?}{A}}_k={\left(\begin{array}{cc}{A}_k& B_k\\ 0_{m?n}& I_m\end{array}\right)}_{\left(n+m\right)?\left(n+m\right)},{\stackrel{?}{B}}_k={\left(\begin{array}{c}{B}_k\\ I_m\end{array}\right)}_{\left(n+m\right)?m},C_k=\left[\begin{array}{cc}C& 0\end{array}\right],{?}_k={\left(\begin{array}{c}{E}_k\\ I_m\end{array}\right)}_{\left(n+m\right)?m}.$

[0095] In the prediction time domain, the predicted value of the output can be calculated by the following equation:

[00029]$\begin{array}{cc}Y\left(k\right)={?}_k?\left(k\right)+{?}_k?U\left(k\right)+{?}_{k}D\left(k\right);& \left(14\right)\end{array}$ [0096] in the formula (14),

[00030]$Y\left(k\right)={\left[\begin{array}{c}?\left(k+1\right)\\ ?\left(k+2\right)\\ ?\left(k+3\right)\\ .Math.\\ ?\left(k+N_p\right)\end{array}\right]}_{N_p?1}?\left(k\right)={\left[\begin{array}{c}{C}_k{\stackrel{?}{A}}_k\\ C_k{\stackrel{?}{A}}_k^2\\ C_k{\stackrel{?}{A}}_k^3\\ .Math.\\ C_k{\stackrel{?}{A}}_k^{N_p}\end{array}\right]}_{N_p?1}?U_k={\left[\begin{array}{c}?u\left(k\right)\\ ?u\left(k+1\right)\\ ?u\left(k+2\right)\\ .Math.\\ ?u\left(k+N_c-1\right)\end{array}\right]}_{N_c?1}$${?}_k=\left[\begin{array}{ccccc}{C}_k{\stackrel{?}{B}}_k& 0& .Math.& .Math.& 0\\ C_k{\stackrel{?}{A}}_K{\stackrel{?}{B}}_k& C_k{\stackrel{?}{B}}_k& .Math.& .Math.& 0\\ C_k{\stackrel{?}{A}}_K^2{\stackrel{?}{B}}_k& C_k{\stackrel{?}{A}}_K{\stackrel{?}{B}}_k& C_k{\stackrel{?}{B}}_k& .Math.& 0\\ .Math.& .Math.& .Math.& ?& .Math.\\ C_k{\stackrel{?}{A}}_K^{N_P-1}{\stackrel{?}{B}}_k& .Math.& .Math.& C_k{\stackrel{?}{A}}_K^{N_P-N_c}{\stackrel{?}{B}}_k& C_k{\stackrel{?}{B}}_k\end{array}\right]$${?}_k=\left[\begin{array}{ccccc}{C}_k{\stackrel{?}{E}}_k& 0& .Math.& .Math.& 0\\ C_k{\stackrel{?}{A}}_K{\stackrel{?}{E}}_k& C_k{\stackrel{?}{E}}_k& .Math.& .Math.& 0\\ C_k{\stackrel{?}{A}}_K^2{\stackrel{?}{E}}_k& C_k{\stackrel{?}{A}}_K{\stackrel{?}{E}}_k& C_k{\stackrel{?}{E}}_k& .Math.& 0\\ .Math.& .Math.& .Math.& ?& .Math.\\ C_k{\stackrel{?}{A}}_K^{N_P-1}{\stackrel{?}{E}}_k& .Math.& .Math.& C_k{\stackrel{?}{A}}_K^{N_P-N_c}{\stackrel{?}{E}}_k& C_k{\stackrel{?}{E}}_k\end{array}\right]$ [0097] where N.sub.p represents a prediction time domain, N.sub.c represents a control time domain, and N.sub.p?N.sub.c.

[0098] The objective function is set as follows to ensure a good control effect while smoothing the input of the control variable:

[00031]$\begin{array}{cc}\min J={.Math.Y\left(k\right)-Y_{\mathrm{ref}}\left(k\right).Math.}_Q^2+{.Math.?U\left(k\right).Math.}_R^2={.Math.}_{j=i+1}^{N_p}{.Math.Y\left(k+j.Math.k\right)-Y_{\mathrm{ref}}\left(k+j.Math.k\right).Math.}_Q^2+{.Math.}_{i=0}^{N_c-1}{.Math.?u\left(k+i.Math.k\right).Math.}_R^2;& \left(15\right)\end{array}$ [0099] where N.sub.p represents a prediction time domain; N.sub.c represents a control time domain; Y(k) represents a predicted value of an output of the vehicle in the prediction time domain N.sub.p, and

[00032]$Y\left(k\right)={\left[\begin{array}{c}?\left(k+1\right)\\ ?\left(k+2\right)\\ ?\left(k+3\right)\\ .Math.\\ ?\left(k+N_p\right)\end{array}\right]}_{N_p?1};$ k represents a current moment; ? represents a lateral position coordinate of the vehicle; Y.sub.ref (k) represents a reference value of the output; ?U(k) represents a sequence of a yaw moment increment at a current moment k; Y(k+j|k) represents a predicted value of the output at a (k+j).sup.th moment when the current moment is k; Y.sub.ref(k+j|k) represents a reference value of the output at the (k+j).sup.th moment when the current moment is k; ?u(k+i|k) represents a yaw moment increment at a (k+i).sup.th moment when the current moment is k; and Q and R represent weight matrices of the output and a control increment, respectively.

[0100] Considering the constraints on the control variable, the limit value of the control increment and the output, the constraints are organized as follows:

[00033]$\{\begin{array}{c}{U}_{\min }?U\left(k\right)?U_{\max }\\ ?U_{\min }??U\left(k\right)??U_{\max }\\ Y_{\min }?Y\left(k\right)?Y_{\max }\end{array}.$

[0101] In this model, the prediction time domain N.sub.p is 5; the control domain N.sub.c is 3;

[00034]$Q=\left[\begin{array}{cc}10& 0\\ 0& 15\end{array}\right];R=\left[\begin{array}{ccc}500& 0& 0\\ 0& 5000& 0\\ 0& 0& 8000\end{array}\right];$

a simulation step T is 0.01 s; a reference value of the output of the vehicle is represented by Y.sub.ref=[?.sub.ref, ?.sub.ref]; ?.sub.ref is a reference value of the yaw angular velocity; and ?.sub.ref represents a reference value of center-of-mass sideslip angle of the vehicle, which is calculated by the linear two-degree-of-freedom model.

[0102] The objective function J is solved to determine the yaw moment increment ?u(k+i|k) at each moment in the control time domain when the objective function J takes the minimum value.

[0103] (2) The additional yaw moment weighting coefficient is determined.

[0104] A coordination controller is constructed. Since the path tracking performance of the vehicle will be affected after performing the yaw moment control, the coordination controller is constructed based on a reinforcement learning deep deterministic policy gradient (DDPG) algorithm to coordinate the conflict between the trajectory tracking and the stability control, where state inputs are yaw velocity w.sub.r, center-of-mass side slip angle ?, lateral displacement Y, heading angle ?, front wheel steering angle ?.sub.f and lateral acceleration a.sub.y; and an action output is the additional yaw moment weighting coefficient ?(0<?<1) In this case, the path tracking accuracy can be ensured with a stable control on the vehicle. The design principles of a reward function of the coordination controller are described as follows.

[0105] (1) At this time, the vehicle is still not destabilized, and the path tracking requirements should be satisfied first, so let r.sub.1=?(e.sub.y.sup.2+e.sub.?.sup.2)+C.sub.1.

[0106] (2) Since the vehicle is in the critical-destabilized state, to avoid the destabilization of the vehicle, it is necessary to carry out stability control on the vehicle to make it return to the stable state, so let

[00035]$r_2=-.Math.\frac{1}{c}\stackrel{?}{?}+\frac{k}{c}?.Math.+C_2.$

[0107] (3) During stability control, the discomfort to the driver caused by the generation of additional yaw moment should be minimized, so let r.sub.3=?|a.sub.y|+C.sub.3.

[0108] Finally, the reward function of the coordination controller is r=r.sub.1+r.sub.2+r.sub.3; [0109] where e.sub.y represents a deviation between an actual lateral displacement of the vehicle and a reference path; e.sub.? represents a deviation between the heading angle and a reference heading angle; a.sub.y represents a lateral acceleration of the vehicle; C.sub.1, C.sub.2, and C.sub.3 are constant; c represents a coefficient; k represents a slope; {dot over (?)} represents a center-of-mass sideslip angular velocity of the vehicle; and ? represents a center-of-mass sideslip angle of the vehicle.

[0110] (3) The path tracking control is performed on the vehicle based on the yaw moment increment and the coordinated control weight.

[0111] The yaw moment increment ?u(k+0|k) at a (k+0).sup.th moment in the control time domain is selected to act on the system to obtain a yaw moment u(k) at a current moment k, where u(k)=u(k?1)+?u(k+0|k).

[0112] An additional yaw moment ?M is determined according to the yaw moment u(k) and the additional yaw moment weighting coefficient ?, where ?M=?u(k).

[0113] (4) A longitudinal force on each wheel of the vehicle is determined based on the additional yaw moment ?M in a drive-brake-combined manner.

[0114] After obtaining the output additional jaw moment, the drive-brake-combined method is adopted to convert it into longitudinal forces on four wheels of the vehicle. By adopting the combination of driving and braking, the longitudinal force on the four wheels be reasonably utilized, which can avoid the tire wear caused by braking or driving of a single wheel, and ensure that the speed of the vehicle is unchanged during the control process, thereby reducing the burden of the driver. The relationship between the longitudinal forces at the center of each wheel is as follows:

[00036]$\{\begin{array}{c}{F}_{\mathrm{xfl}}+F_{\mathrm{xrl}}=F_{\mathrm{xfr}}+F_{\mathrm{xrr}}\\ F_{\mathrm{xfl}}=F_{\mathrm{xrl}}\\ F_{\mathrm{xfr}}=F_{\mathrm{xrr}}\end{array};$ [0115] where F.sub.xft, F.sub.xrl, F.sub.xfr and F.sub.xrr represent longitudinal forces on a left front wheel, a left rear wheel, a right front wheel and a right rear wheel, respectively, and the longitudinal forces on the four wheels satisfy the following constraints:

[0116] (1)

[00037]$.Math.?M.Math.=.Math.F_{\mathrm{xfl}}+F_{\mathrm{xrl}}.Math.?\frac{d}{2}-.Math.F_{\mathrm{xfr}}+F_{\mathrm{xrr}}.Math.?\frac{d}{2},$

where when ?M>0, F.sub.xft=F.sub.xrl<0 and F.sub.xfr=F.sub.xrr>0; when ?M<0, F.sub.xft=F.sub.xrl>0 and F.sub.xfr=F.sub.xrr<0; and d represents a wheelbase.

[0117] (2) |F.sub.xi|??F.sub.zi, where F.sub.xi represents a longitudinal force of each wheel; i represents a wheel number, corresponding to the left front wheel, the left rear wheel, the right front wheel and the right rear wheel, respectively; F.sub.zi represents a vertical load of each wheel; and ? represents an adhesion coefficient of road surface.

[0118] (3)

[00038]$.Math.F_{xi}.Math.?\frac{T_{\max }}{r},$

where T.sub.max represents the maximum torque of the motor, and r represents a radius of awheel.

[0119] By using the drive-brake combined strategy, the longitudinal force acting on the center of each wheel can be finally obtained to achieve the distribution of the additional yaw moment.

[0120] In an embodiment, when the vehicle is in the destabilized state, the second yaw moment increment of the vehicle is determined, and the path tracking control is performed based on the second yaw moment increment through the following steps.

[0121] An objective function J is constructed:

[00039]$\min J={.Math.Y\left(k\right)-Y_{ref}\left(k\right).Math.}_Q^2+{.Math.?U\left(k\right).Math.}_R^2={.Math.}_{j=i+1}^{N_p}{.Math.Y\left(k+j.Math.k\right)-Y_{ref}\left(k+j.Math.k\right).Math.}_Q^2+{.Math.}_{i=0}^{N_c-1}{.Math.?u\left(k+i.Math.k\right).Math.}_R^2;$ [0122] where N.sub.p represents a prediction time domain; N.sub.c represents a control time domain; Y(k) represents a predicted value of an output of the vehicle in the prediction time domain N.sub.p; k represents a current moment; Y.sub.ref ( ) represents a reference value of the output; ?U(k) represents a sequence of a yaw moment increment at a current moment k; Y(k+j|k) represents a predicted value of the output at a (k+j).sup.th moment when the current moment is k; Y.sub.ref(k+j|k) represents a reference value of the output at the (k+j).sup.th moment when the current moment is k; ?u(k+i|k) represents a yaw moment increment at a (k+i).sup.th moment when the current moment is k; and Q and R represent weight matrices of the output and a control increment, respectively.

[0123] The objective function J; is solved to determine the yaw moment increment ?u(k+i|k) at each moment within the control time domain when the objective function J reaches a minimum value.

[0124] The yaw moment increment ?u(k+0|k) at a (k+0).sup.th moment in the control time domain is selected to act on the system to obtain a yaw moment u(k) at a current moment k, where u(k)=u(k?1)+?u(k+0|k).

[0125] An additional yaw moment ?M is determined according to the yaw moment u(k) and the additional yaw moment weighting coefficient ?, where ?M=?u(k).

[0126] A longitudinal force on each wheel of the vehicle is determined based on the additional yaw moment ?M in a drive-brake-combined manner.

[0127] In this embodiment, the specific calculation of the longitudinal force on each wheel of the vehicle refers to the calculation method mentioned in the aforementioned embodiment where the vehicle is in a critical-destabilized state, and will not be repeated herein.

[0128] Based on the same inventive conception as the aforementioned path tracking control method, the present disclosure also provides a path tracking control device corresponding thereto. FIG. 2 shows a block diagram of a structure of a path tracking control device for an intelligent electric vehicle. The device includes a lateral stability state determination module 210 and a path tracking control module 220.

[0129] The lateral stability state determination module 210 is configured to determine a lateral stability state of a vehicle, where the lateral stability state comprises a stable state, a critical-destabilized state, and a destabilized state.

[0130] The path tracking control module 220 is configured to perform path tracking control on the vehicle according to the lateral stability state of the vehicle, determine a front wheel steering increment of the vehicle when the vehicle is in the stable state, determine a front wheel steering of the vehicle based on the front wheel steering increment of the vehicle, determine a first yaw moment increment and an additional yaw moment weighting coefficient of the vehicle when the vehicle is in the critical-destabilized state, determine a first additional yaw moment based on the first yaw moment increment and the additional yaw moment weighting coefficient, determine a first longitudinal force on each wheel of the vehicle based on the first additional yaw moment, determine a second yaw moment increment of the vehicle when the vehicle is in the destabilized state, determine a second additional yaw moment based on the second yaw moment increment, and determine a second longitudinal force on each wheel based on the second additional yaw moment.

[0131] The present disclosure also provides a computer-readable storage medium. The computer-readable storage medium is configured to store a computer program, and the computer program is configured to be executed by a processor to implement the aforementioned path tracking control method.

[0132] The present disclosure also provides a computer program product, which includes a computer program or a computer instruction, where the computer program or the computer instruction is configured to be executed by a processor to implement the aforementioned path tracking control method.

[0133] In the embodiments provided herein, it should be understood that the disclosed devices and methods may also be realized in other ways. The device embodiments described above are merely schematic. For example, the flow chart and the block diagram in the accompanying drawings show the architecture, functionality, and operation of possible implementations of devices, methods, and computer program products according to a plurality of embodiments of the present application. In this regard, each box in the flow chart or the block diagram may represent a module, a program segment, or a portion of codes, which contains one or more executable instructions for implementing a prescribed logical function. It should also be noted that in some implementations as substitutions, the functions labeled in the boxes may also occur in a different order than those labeled in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, or in reverse order, depending on the involved function. Moreover, each box in the block diagram and/or flowchart, and combinations of boxes in the block diagram and/or flowchart, may be implemented with a specialized hardware-based system that performs the specified function or action, or may be implemented with a combination of specialized hardware and computer instructions.

[0134] In addition, the various function modules in various embodiments of the present application may be integrated together to form a separate part, or the individual modules may stand alone, or two or more modules may be integrated to form a separate part.

[0135] The functions may be stored in a computer-readable storage medium if implemented as software function modules and sold or used as a separate product. Based on this, the technical solutions of the present application may be embodied in the form of a software product, which is essentially or in part a contribution to the prior art, or part of the technical solution may be embodied in the form of a software product. The computer program product is stored in a storage medium, and includes a series of instructions to make a computer device (a personal computer, a server, or a network device, etc.) to carry out all or part of the steps of the method described in the various embodiments of the present application. The aforementioned storage medium can store program codes, and can be a USB flash drive, a removable hard disk, a read-only memory (ROM), a random-access memory (RAM), a diskette or a compact disc.

[0136] Described above are merely various embodiments of the present application, but the scope of protection of the present application is not limited thereto. Any changes or substitutions within the scope of the technology disclosed in the present application can be easily obtained by one of ordinary skilled in the art, and shall be included in the scope of protection of the present application. Therefore, the scope of protection of the present application shall be defined by the claims.