METHOD AND SYSTEM FOR VIRTUALLY COUPLED TRAIN SET CONTROL

Abstract

A method and system for virtually coupled train set (VCTS) control is provided. The method includes following steps: determining whether to execute a backup control strategy based on an actual state for a current cycle of each train unit and a target state sequence for a first preset number of cycles before the current cycle to obtain a first determination result; if the first determination result is yes, executing the backup control strategy to control each train unit; if the first determination result is no, calculating the target state sequence for the current cycle of each train unit according to a position or calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship; and controlling each train unit according to the target state sequence for the current cycle of each train unit, respectively.

Claims

1. A method for virtually coupled train set (VCTS) control, comprising: acquiring an actual state for a current cycle of each train unit in VCTS; determining whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first predetermined number of cycles before the current cycle, to obtain a first determination result; wherein the backup control strategy comprises a control strategy for tracking a recommended driving curve by a first train unit and a control strategy for tracking a i-th train unit by a (i+1)-th train unit, wherein a value of i is greater than or equal to 1; executing the backup control strategy to control each train unit, if the first determination result is yes; executing following operations, if the first determination result is no: determining whether synchronization of each train unit in the VCTS meets a predetermined condition, to obtain a second determination result; calculating a target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit, if the second determination result is yes; calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit, if the second determination result is no; controlling each train unit according to the target state sequence for the current cycle of each train unit, respectively.

2. The method according to claim 1, wherein the determining whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first predetermined number of cycles before the current cycle to obtain a first determination result comprises: determining whether a flag bit of a first cycle before the current cycle is displayed normally, to obtain a third determination result; determining whether a difference between an actual speed for the current cycle of each train unit and a first target speed in a target state sequence for n cycles before the current cycle is less than a speed difference threshold, if the third determination result is yes; wherein if the difference between the speed of for the current cycle of each train unit and the first target speed in the target state sequence in for the n cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and a flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal; determining whether a difference between the actual speed for the current cycle of each train unit and a first target speed in a target state sequence for m cycles before the current cycle is less than the speed difference threshold, if the third determination result is no; wherein if the difference between the speed for of the current cycle of each train unit and the first target speed in the target state sequence in for the m cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and the flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal, wherein n and m are values of the first predetermined number in different situations, and mn.

3. The method according to claim 1, wherein the determining whether synchronization of each train unit in the VCTS meets a predetermined condition to obtain a second determination result comprises: confirming that the second determination result is yes when time index deviations between any two adjacent train units in the VCTS are all less than a time index deviation threshold; confirming that the second determination result is no when the time index deviations between any two adjacent train units in the VCTS are not all less than the time index deviation threshold.

4. The method according to claim 1, wherein the calculating a target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state as: ${\stackrel{}{s}}_{i,0|k}=s_{i,k},{\stackrel{}{v}}_{i,0|k}=V\left({\stackrel{}{s}}_{i,0|k}\right),i=1,2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, v.sub.i,0k is an initial target speed in the target state sequence for the current cycle of the i-th train unit, s.sub.i,k and v.sub.i,k are an actual position and an actual speed for the current cycle of the i-th train unit, respectively, V( ) is a calculation function of a target speed, $V\left({\stackrel{}{s}}_{i,0|k}\right)={\stackrel{}{v}}_{i,q}+\frac{\left({\stackrel{}{s}}_{i,0|k}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,q+1}-{\stackrel{}{s}}_{i,q}},s.t.{\stackrel{}{s}}_{i,q}{\stackrel{}{s}}_{i,0|k}<{\stackrel{}{s}}_{i,q+1},$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on a recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and I is a number of train units in the VCTS; calculating, based on the initial target position and the initial target speed, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1|k}={\stackrel{}{s}}_{i,j|k}+{\stackrel{}{}}_{i,j|k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1|k}=V\left({\stackrel{}{s}}_{i,j+1|k}\right)$ wherein s.sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, N is a number of target states in the target state sequence, is a sampling interval time, V( ) is a calculation function of the target speed, $V\left({\stackrel{}{s}}_{i,j+1|k}\right)={\stackrel{}{v}}_{i,p}+\frac{({\stackrel{}{s}}_{i,j+1|k-{\stackrel{}{s}}_{i,p})\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{i,p}\right)}}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.{\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1|k}<{\stackrel{}{s}}_{i,p+1},$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

5. The method according to claim 1, wherein the calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state of each train unit as: ${\stackrel{}{s}}_{i,0|k}=s_{i,k},{\stackrel{}{v}}_{1,0|k}=V\left({\stackrel{}{s}}_{1,0|k}\right),{\stackrel{}{v}}_{i^{},0|k}=V\left({\stackrel{}{s}}_{i^{},0|k}\right)+v_a\left(0\right),i^{}=2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, s.sub.i,k is an actual position for the current cycle of the i-th train unit, v.sub.1,0|k and v.sub.i,0|k are initial target speeds in target state sequences for the current cycle of the first train unit and a i-th train unit, respectively, v.sub.a(0) is an adjustment amount of the initial target speed in the target state sequence for the current cycle of the i-th train unit, V( ) is a calculation function of a target speed, $\begin{array}{cc}V\left({\stackrel{}{s}}_{i,0|k}\right)={\stackrel{}{}}_{i,q}+\frac{\left({\stackrel{}{s}}_{i,0|k}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,q+1}-{\stackrel{}{s}}_{i,q}},s.t.& {\stackrel{}{s}}_{i,q}{\stackrel{}{s}}_{i,0|k}<{\stackrel{}{s}}_{i,q+1}\end{array},$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and I is a number of train units; calculating, based on the initial target state, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1|k}={\stackrel{}{s}}_{i,j|k}+{\stackrel{}{v}}_{i,j|k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1|k}=V\left({\stackrel{}{s}}_{i,j+1|k}\right)+v_a\left(j\right))$ wherein s.sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, t is a sampling interval time, V( ) is a calculation function of the target speed, $\begin{array}{cc}V\left({\stackrel{}{s}}_{i,j+1|k}\right)={\stackrel{}{v}}_{i,p}+\frac{\left({\stackrel{}{s}}_{i,j+1|k}-{\stackrel{}{s}}_{i,p}\right)\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{i,p}\right)}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.& {\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1|k}<{\stackrel{}{s}}_{i,p+1},\end{array}$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and v.sub.a(j) is an adjustment amount of the j-th target speed in the target state sequence; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

6. The method according to claim 5, wherein a calculation formula of the adjustment amount of the j-th target speed in the target state sequence is: $v_a\left(j\right)$ $=\{\begin{array}{cc}\left(c_a+j\right)a_a,& \left(c_a+j\right){}_a/a_a/\\ {}_a,& {}_a/a_a/\left(c_a+j\right)T_a/-{}_a/a_a/\\ \left(T_a/-\left(c_a+j\right)\right)a_a& \left(c_a+j\right)T_a/-{}_a/a_a/\\ 0& \mathrm{else}\end{array}$ wherein c.sub.a is a number of cycles between the current cycle and a cycle during which adjustment is performed based on the adjustment amount for a first time, a.sub.a is a predetermined adjustment acceleration, is the sampling interval time, v.sub.a is a maximum adjustment speed, and T.sub.a is a predetermined adjustment time.

7. The method according to claim 1, wherein prior to the controlling each train unit according to the target state sequence for the current cycle of each train unit, the method further comprises: executing following operations, when a difference between an actual speed and an emergency braking intervention (EBI) speed for a second predetermined number of cycles between a previous moment and the current cycle of each train unit does not meet a predetermined condition: adjusting a target speed and a target position in the target state sequence for the current cycle of each train unit based on following formulas: ${{\stackrel{}{v}}^{}}_{i,j|k}={\stackrel{}{v}}_{i,j|k}-k_e{v}_e$ ${{\stackrel{}{s}}^{}}_{i,j+1|k}={{\stackrel{}{s}}^{}}_{i,j|k}+{{\stackrel{}{v}}^{}}_{i,j|k}$ ${}_e=\frac{1}{l}\underset{l=k-n+1}{\overset{k}{.Math.}}{}_{i,l}-{\stackrel{}{}}_{i,\mathrm{ebi}}+{}_e$ wherein v.sub.i,j|k is a j-th adjusted target speed in the target state sequence for the current cycle of the i-th train unit, s.sub.i,j|k and s.sub.i,j+1|k are j-th and (j+1)-th adjusted target positions in the target state sequence for the current cycle of the i-th train unit, respectively, v.sub.i,j|k is a j-th target speed in the target state sequence for the current cycle of the i-th train unit, k.sub.e is a compensation proportional coefficient, is a sampling interval time, v.sub.e is a speed compensation amount, v.sub.i,l is an actual speed at a l-th moment before the current cycle, l1 is a second predetermined number, v.sub.i,ebi is an EBI speed of the i-th train unit, and v.sub.e is an EBI speed margin value; performing differential calculation on an adjusted target speed in the target state sequence for the current cycle of each train unit to obtain an adjusted target acceleration in the target state sequence for the current cycle of each train unit.

8. A system for virtually coupled train set (VCTS) control wherein the system is applied to the method according to claim 1, and the system comprises: a state acquiring module, configured to acquire an actual state for a current cycle of each train unit in VCTS; a first determination module, configured to determine whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first predetermined number of cycles before the current cycle to obtain a first determination result; wherein the backup control strategy comprises a control strategy for tracking a recommended driving curve by a first train unit and a control strategy for tracking a i-th train unit by a (i+1)-th train unit, wherein a value of i is greater than or equal to 1; a first control module, configured to execute the backup control strategy to control each train unit, if the first determination result is yes; a second determination module, configured to determine whether synchronization of each train unit in the VCTS meets a predetermined condition to obtain a second determination result, if the first determination result is no; a first target state sequence calculating module, configured to calculate the target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit, if the second determination result is yes; a second state sequence calculating module, configured to calculate the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit, if the second determination result is no; a second control module, configured to control each train unit according to the target state sequence for the current cycle of each train unit, respectively.

9. The system according to claim 8, wherein the determining whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first predetermined number of cycles before the current cycle to obtain a first determination result comprises: determining whether a flag bit of a first cycle before the current cycle is displayed normally, to obtain a third determination result; determining whether a difference between an actual speed for the current cycle of each train unit and a first target speed in a target state sequence for n cycles before the current cycle is less than a speed difference threshold, if the third determination result is yes; wherein if the difference between the speed of for the current cycle of each train unit and the first target speed in the target state sequence in for the n cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and a flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal; determining whether a difference between the actual speed for the current cycle of each train unit and a first target speed in a target state sequence for m cycles before the current cycle is less than the speed difference threshold, if the third determination result is no; wherein if the difference between the speed for of the current cycle of each train unit and the first target speed in the target state sequence in for the m cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and the flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal, wherein n and m are values of the first predetermined number in different situations, and mn.

10. The system according to claim 8, wherein the determining whether synchronization of each train unit in the VCTS meets a predetermined condition to obtain a second determination result comprises: confirming that the second determination result is yes when time index deviations between any two adjacent train units in the VCTS are all less than a time index deviation threshold; confirming that the second determination result is no when the time index deviations between any two adjacent train units in the VCTS are not all less than the time index deviation threshold.

11. The system according to claim 8, wherein the calculating a target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state as: ${\stackrel{}{s}}_{i,0|k}=s_{i,k},{\stackrel{}{v}}_{i,0|k}=V\left({\stackrel{}{s}}_{i,0|k}\right),i=1,2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, v.sub.i,0|k is an initial target speed in the target state sequence for the current cycle of the i-th train unit, s.sub.i,k and v.sub.i,k are an actual position and an actual speed for the current cycle of the i-th train unit, respectively, V( ) is a calculation function of a target speed, $V\left({\stackrel{}{s}}_{i,0|k}\right)={\stackrel{}{v}}_{i,q}+\frac{\left({\stackrel{}{s}}_{i,0|k}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,q+1}-{\stackrel{}{s}}_{i,q}},s.t.{\stackrel{}{s}}_{i,q}{\stackrel{}{s}}_{i,0|k}<{\stackrel{}{s}}_{i,q+1,}$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on a recommended driving curve of the i-th train unit, respectively, v.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and/is a number of train units in the VCTS; calculating, based on the initial target position and the initial target speed, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1|k}={\stackrel{}{s}}_{i,j|k}+{\stackrel{}{v}}_{i,j|k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1|k}=V\left({\stackrel{}{s}}_{i,j+1|k}\right)$ wherein s.sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, Nis a number of target states in the target state sequence, is a sampling interval time, V( ) is a calculation function of the target speeu, $V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)={\stackrel{}{v}}_{i,p}+\frac{\left({\stackrel{}{s}}_{i,j+1.Math.k}-{\stackrel{}{s}}_{i,p}\right)\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{i,p}\right)}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.{\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1.Math.k}<{\stackrel{}{s}}_{i,p+1},$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

12. The system according to claim 8, wherein the calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state of each train unit as: ${\stackrel{}{s}}_{i,0.Math.k}=s_{i,k},{\stackrel{}{v}}_{1,0.Math.k}=V\left({\stackrel{}{s}}_{1,0.Math.k}\right),{\stackrel{}{v}}_{i^{},0.Math.k}=V\left({\stackrel{}{s}}_{i,0.Math.k}\right)+v_a\left(0\right),i^{}=2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, s.sub.i,k is an actual position for the current cycle of the i-th train unit, v.sub.1,0|k and v.sub.i,0|k are initial target speeds in target state sequences for the current cycle of the first train unit and a i-th train unit, respectively, v.sub.a(0) is an adjustment amount of the initial target speed in the target state sequence for the current cycle of the i-th train unit, V( ) is a calculation function of a target speed, $V\left({\stackrel{}{s}}_{i,0.Math.k}\right)={\stackrel{}{v}}_{i,q}+\frac{\left({\stackrel{\_}{s}}_{i,0.Math.k}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,q+1}-{\stackrel{}{s}}_{i,q}},s.t.{\hat{s}}_{i,q}{\stackrel{}{s}}_{i,0.Math.k}<{\stackrel{}{s}}_{i,q+1},$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and I is a number of train units; calculating, based on the initial target state, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1.Math.k}={\stackrel{}{s}}_{i,j.Math.k}+{\stackrel{}{v}}_{i,j.Math.k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1.Math.k}=V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)+v_a\left(j\right)$ wherein .sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, is a sampling interval time, V( ) is a calculation function of the target speed, $V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)={\stackrel{}{v}}_{i,p}+\frac{\left({\stackrel{}{s}}_{i,j+1.Math.k}-{\stackrel{}{s}}_{i,p}\right)\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{i,p}\right)}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.{\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1.Math.k}<{\stackrel{}{s}}_{i,p+1},$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and v.sub.a(j) is an adjustment amount of the j-th target speed in the target state sequence; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

13. The system according to claim 12, wherein a calculation formula of the adjustment amount of the j-th target speed in the target state sequence is: $v_a\left(j\right)=\{\begin{array}{cc}\left(c_a+j\right)a_a,& \left(c_a+j\right)v_a/a_a/\\ v_a,& v_a/a_a/\left(c_a+j\right)T_a/-v_a/a_a/\\ \left(T_a/-\left(c_a+j\right)\right)a_a& \left(c_a+j\right)T_a/-v_a/a_a/\\ 0& \mathrm{else}\end{array}$ wherein c.sub.a is a number of cycles between the current cycle and a cycle during which adjustment is performed based on the adjustment amount for a first time, a.sub.a is a predetermined adjustment acceleration, is the sampling interval time, v.sub.a is a maximum adjustment speed, and T.sub.a is a predetermined adjustment time.

14. The system according to claim 8, wherein prior to the controlling each train unit according to the target state sequence for the current cycle of each train unit, the method further comprises: executing following operations, when a difference between an actual speed and an emergency braking intervention (EBI) speed for a second predetermined number of cycles between a previous moment and the current cycle of each train unit does not meet a predetermined condition: adjusting a target speed and a target position in the target state sequence for the current cycle of each train unit based on following formulas: ${{\stackrel{\_}{v}}^{}}_{i,j.Math.k}={\stackrel{}{v}}_{i,j.Math.k}-k_e{v}_e$ ${{\stackrel{\_}{s}}^{}}_{i,j+1.Math.k}={{\stackrel{\_}{s}}^{}}_{i,j.Math.k}+{{\stackrel{\_}{v}}^{}}_{i,j.Math.k}$ $v_e=\frac{1}{l}\underset{l=k-n+1}{\overset{k}{.Math.}}{v}_{i,l}-{\stackrel{}{v}}_{i,ebi}+v_e$ wherein v.sub.i,j|k is a j-th adjusted target speed in the target state sequence for the current cycle of the i-th train unit, s.sub.i,j|k and s.sub.i,j+1|k are j-th and (j+1)-th adjusted target positions in the target state sequence for the current cycle of the i-th train unit, respectively, v.sub.i,j|k is a j-th target speed in the target state sequence for the current cycle of the i-th train unit, k.sub.e is a compensation proportional coefficient, t is a sampling interval time, v.sub.e is a speed compensation amount, v.sub.i,l is an actual speed at a l-th moment before the current cycle, l1 is a second predetermined number, v.sub.i,ebi is an EBI speed of the i-th train unit, and v.sub.e is an EBI speed margin value; performing differential calculation on an adjusted target speed in the target state sequence for the current cycle of each train unit to obtain an adjusted target acceleration in the target state sequence for the current cycle of each train unit.

15. An electronic device comprising a memory, a processor and a computer program product, comprising a computer readable hardware storage device having computer readable program code stored therein, said program code executable by a processor of a computer system to implement a method stored in the memory and executable on the processor, wherein the processor implements the method according to claim 1 when executing the computer program.

16. The electronic device according to claim 15, wherein the determining whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first predetermined number of cycles before the current cycle to obtain a first determination result comprises: determining whether a flag bit of a first cycle before the current cycle is displayed normally, to obtain a third determination result; determining whether a difference between an actual speed for the current cycle of each train unit and a first target speed in a target state sequence for n cycles before the current cycle is less than a speed difference threshold, if the third determination result is yes; wherein if the difference between the speed of for the current cycle of each train unit and the first target speed in the target state sequence in for the n cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and a flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal; determining whether a difference between the actual speed for the current cycle of each train unit and a first target speed in a target state sequence for m cycles before the current cycle is less than the speed difference threshold, if the third determination result is no; wherein if the difference between the speed for of the current cycle of each train unit and the first target speed in the target state sequence in for the m cycles before a the current cycle is less than the speed difference threshold, the first determination result is yes, and the flag bit of the current cycle is set as normal, otherwise, the first determination result is no, and the flag bit of the current cycle is set as abnormal, wherein n and m are values of the first predetermined number in different situations, and mn.

17. The electronic device according to claim 15, wherein the determining whether synchronization of each train unit in the VCTS meets a predetermined condition to obtain a second determination result comprises: confirming that the second determination result is yes when time index deviations between any two adjacent train units in the VCTS are all less than a time index deviation threshold; confirming that the second determination result is no when the time index deviations between any two adjacent train units in the VCTS are not all less than the time index deviation threshold.

18. The electronic device according to claim 15, wherein the calculating a target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state as: ${\stackrel{}{s}}_{i,0.Math.k}=s_{i,k},{\stackrel{}{v}}_{i,0.Math.k}=V\left({\stackrel{}{s}}_{i,0.Math.k}\right),i=1,2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, v.sub.i,0|k is an initial target speed in the target state sequence for the current cycle of the i-th train unit, s.sub.i,k and v.sub.i,k are an actual position and an actual speed for the current cycle of the i-th train unit, respectively, V( ) is a calculation function of a target speed, $V\left({\stackrel{}{s}}_{i,0.Math.k}\right)={\stackrel{}{v}}_{i,q}+\frac{\left({\stackrel{}{s}}_{\left(i,0.Math.k\right)}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,q+1}-{\stackrel{}{s}}_{i,q}},s.t.{\stackrel{}{\hat{s}}}_{i,q}{\stackrel{}{s}}_{i,0.Math.k}<{\stackrel{}{s}}_{i,q+1},$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on a recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and/is a number of train units in the VCTS; calculating, based on the initial target position and the initial target speed, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1.Math.k}={\stackrel{}{s}}_{i,j.Math.k}+{\stackrel{}{v}}_{i,j.Math.k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1.Math.k}=V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)$ wherein s.sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, Nis a number of target states in the target state sequence, t is a sampling interval time, V( ) is a calculation function of the target speed, $V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)={\stackrel{}{v}}_{i,p}+\frac{\left({\stackrel{}{s}}_{i,j+1.Math.k}-{\stackrel{}{s}}_{i,p}\right)\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{ip}\right)}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.{\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1.Math.k}<{\stackrel{}{s}}_{i,p+1},$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

19. The electronic device according to claim 15, wherein the calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit comprises: determining, according to an actual position and an actual speed for the current cycle of each train unit, an initial target state of each train unit as: ${\stackrel{}{s}}_{i,0.Math.k}=s_{i,k},{\stackrel{}{v}}_{1,0.Math.k}=V\left({\stackrel{}{s}}_{1,0.Math.k}\right),{\stackrel{}{v}}_{i^{},0.Math.k}=V\left({\stackrel{}{s}}_{i,0.Math.k}\right)+v_a\left(0\right),i^{}=2,3,.Math.,I;$ wherein s.sub.i,0|k is an initial target position in a target state sequence for the current cycle of the i-th train unit, s.sub.i,k is an actual position for the current cycle of the i-th train unit, v.sub.1,0|k and v.sub.i,0|k are initial target speeds in target state sequences for the current cycle of the first train unit and a i-th train unit, respectively, v.sub.a(0) is an adjustment amount of the initial target speed in the target state sequence for the current cycle of the i-th train unit, V( ) is a calculation function of a target speed, $V\left({\stackrel{}{s}}_{i,0.Math.k}\right)={\stackrel{}{v}}_{i,q}+\frac{\left({\stackrel{}{s}}_{i,0.Math.k}-{\stackrel{}{s}}_{i,q}\right)\left({\stackrel{}{v}}_{i,q+1}-{\stackrel{}{v}}_{i,q}\right)}{{\stackrel{}{s}}_{i,+1}-{\stackrel{}{s}}_{i,q}},s.t.{\hat{s}}_{i,q}{\stackrel{}{s}}_{i,0.Math.k}<{\stackrel{}{s}}_{i,q+1},$ .sub.i,q and .sub.i,q+1 represent q-th and (q+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,q and {circumflex over (v)}.sub.i,q+1 represent q-th and (q+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and/is a number of train units; calculating, based on the initial target state, a target position and a target speed in the target state sequence for the current cycle of each train unit based on following formulas: ${\stackrel{}{s}}_{i,j+1.Math.k}={\stackrel{}{s}}_{i,j.Math.k}+{\stackrel{}{v}}_{i,j.Math.k},j=0,1,.Math.,N-1$ ${\stackrel{}{v}}_{i,j+1.Math.k}=V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)+v_a\left(j\right)$ wherein s.sub.i,j+1|k and v.sub.i,j+1|k are a (j+1)-th target position and a (j+1)-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, s.sub.i,j|k and v.sub.i,j|k are a j-th target position and a j-th target speed in the target state sequence for the current cycle of the i-th train unit, respectively, t is a sampling interval time, V( ) is a calculation function of the target speed, $V\left({\stackrel{}{s}}_{i,j+1.Math.k}\right)={\stackrel{}{v}}_{i,p}+\frac{\left({\stackrel{}{s}}_{i,j+1.Math.k}-{\stackrel{}{s}}_{i,p}\right)\left({\stackrel{}{v}}_{i,p+1}-{\stackrel{}{v}}_{i,p}\right)}{{\stackrel{}{s}}_{i,p+1}-{\stackrel{}{s}}_{i,p}},s.t.{\stackrel{}{s}}_{i,p}{\stackrel{}{s}}_{i,j+1.Math.k}<{\stackrel{}{s}}_{i,p+1},$ .sub.i,p and .sub.i,p+1 represent p-th and (p+1)-th recommended position values on the recommended driving curve of the i-th train unit, respectively, {circumflex over (v)}.sub.i,p and {circumflex over (v)}.sub.i,p+1 represent p-th and (p+1)-th recommended speed values on the recommended driving curve of the i-th train unit, respectively, and v.sub.a(j) is an adjustment amount of the j-th target speed in the target state sequence; performing differential calculation on the target speed in the target state sequence for the current cycle of each train unit to obtain a target acceleration in the target state sequence for the current cycle of each train unit.

20. A computer-readable storage medium storing a computer program, wherein the computer program, when executed, implements the method according to claim 1.

Description

BRIEF DESCRIPTION

[0057] Some of the embodiments will be described in detail, with reference to the following figures, wherein like designations denote like members, wherein:

[0058] FIG. 1 is a flow chart of the method for VCTS control according to an embodiment of the present disclosure;

[0059] FIG. 2 is a process flow diagram of generating a target state sequence according to an embodiment of the present disclosure;

[0060] FIG. 3 is a graph showing speed and distance of preceding and following train units with time under an ideal control effect according to an embodiment of the present disclosure;

[0061] FIG. 4 is a graph showing speed and distance of preceding and following train units with time under the condition that each train unit still operates according to its own recommended speed curve and the departure of the following train unit is delayed, according to an embodiment of the present disclosure;

[0062] FIG. 5 is a graph showing speed and distance of train units with time under the condition that VCTS is controlled by the present disclosure and the departure of the following train unit is delayed, according to an embodiment of the present disclosure;

[0063] FIG. 6 is a graph showing speed and distance of preceding and following train units with time under the condition that each train unit still operates according to its own recommended speed curve and control errors of the following train unit are accumulated resulting in unsynchronization (the following train unit is relatively slower than the preceding train unit) according to an embodiment of the present disclosure;

[0064] FIG. 7 is a graph showing speed and distance of train units with time under the condition that VCTS is controlled by the present disclosure and the control errors of the following train unit are accumulated resulting in unsynchronization (the following train unit is relatively slower than the preceding train unit) according to an embodiment of the present disclosure;

[0065] FIG. 8 is a graph showing speed and distance of preceding and following train units with time under the condition that each train unit still operates according to its own recommended speed curve and the control errors of the following train unit are accumulated resulting in unsynchronization and overspeed (the following train unit is relatively faster than the preceding train unit) according to an embodiment of the present disclosure; and

[0066] FIG. 9 is a graph showing speed and distance of train units with time under the condition that VCTS is controlled by the present disclosure and the control errors of the following train unit are accumulated resulting in unsynchronization (the following train unit is relatively faster than the preceding train unit) according to an embodiment of the present disclosure.

DETAILED DESCRIPTION

[0067] The technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in the embodiments of the present disclosure. Apparently, the described embodiments are only some embodiments of the present disclosure, rather than all of the embodiments. Based on the embodiment of the present disclosure, all other embodiments obtained by those skilled in the conventional art without creative efforts shall fall within the scope of protection of the present disclosure.

[0068] An objective of embodiments of the present disclosure is to provide a method and a system for VCTS control, so as to ensure that a synchronous operation relationship among all train units can be maintained, while the train units track their respective recommended driving curves.

[0069] To make the above objective, features and advantages of the present disclosure more apparent, the present disclosure will be further explained in detail below with references to the drawings and detailed description.

Embodiment 1

[0070] Embodiment 1 of the present disclosure provides a method for VCTS control. As shown in FIG. 1, the method includes steps of: [0071] acquiring an actual state for a current cycle of each train unit in VCTS; [0072] determining whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first preset number of cycles before the current cycle to obtain a first determination result; where the backup control strategy includes a control strategy for tracking a recommended driving curve by a first train unit and a control strategy for tracking the i-th train unit by the (i+1)-th train unit, where a value of i is greater than or equal to 1.

[0073] If the first determination result is yes, the backup control strategy is executed to control each train unit.

[0074] If the first determination result is no, following operations is executed: [0075] determining whether synchronization of each train unit in VCTS meets a preset condition to obtain a second determination result; [0076] if the second determination result is yes, calculating the target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit; [0077] if the second determination result is no, calculating the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit; [0078] controlling each train unit according to the target state sequence for the current cycle of each train unit, respectively.

[0079] As shown in FIG. 2, the above method specifically includes following steps.

[0080] The position, speed and acceleration of the train unit i at time k are received and represented by s.sub.i,k, v.sub.i,k, and a.sub.i,k, respectively, where i=1 or 2, representing a preceding train unit or a following train unit. First, the operation stage of VCTS is determined. When all the following conditions are true, VCTS is considered to be in the station arrival stage, otherwise, VCTS is considered to be in the inter-station operation stage: (1) the positions of all train units are within the station parking area; (2) the speeds of all train units are zero; (3) the countdown for parking is not zero; (4) no departure signal is received from the station.

[0081] In the case of the station arrival stage, other modules are responsible for handling the station operation function (the module can use the design in the existing research and inventions, and the present disclosure does not involve the design of the module for realizing this part of the function).

[0082] In the case of the inter-station operation stage, the control target sequence of each train unit is calculated.

[0083] The target sequence in the future is represented by s.sub.i,j|k, v.sub.i,j|k and .sub.ij|k, where j=1, 2, . . . , N, N represents a prediction horizon. The recommended speed curve calculated offline is represented by .sub.i,j|k, {circumflex over (v)}.sub.i,j|k and .sub.i,j|k, and the sampling interval time is .

[0084] The process proceeds to function 1 to determine whether the control errors of all train units are within an allowable range. When the control errors of all train units meet the requirement, which means that they can track the recommended driving curve well, the problems of asynchronous operation resulted from the accumulation of the control errors can be effectively alleviated by adjusting the target sequence. However, if the control effect of the train unit is not good enough to keep up with the target sequence, the difference between the actual state of the train and the target sequence is large, and thus it is difficult to determine the influence of adjusting the target sequence on the actual state, and it is difficult to achieve the goal of synchronous operation of train units by adjusting the target sequence. Therefore, first, it is determined whether the train can track the target (speed) well, so as to decide the ways of adjustment.

[0085] The function 1 is executed as follows.

[0086] (1) The difference v.sub.i,kv.sub.i,1|k-1 between the actual speed and the target speed is calculated.

[0087] (2) If the displayed control effect of the flag bit of a previous cycle is normal, it is determined whether the difference in n previous cycles meets the requirements of the threshold V.sub.v.sub.r:v.sub.i,jv.sub.i,1|j-1V.sub.v.sub.r, j=k,k1, . . . , kn+1.

[0088] (3) If the displayed control effect of the flag bit of a previous cycle is abnormal, it is determined whether the difference within m previous cycles meets the requirements of the threshold V.sub.v.sub.r, v.sub.i,jv.sub.i,1|j-1V.sub.v.sub.r, j=k, k1, . . . , km+1.

[0089] (4) If the difference meets the requirement, the flag bit with the normal control effect is output, turning to the main control strategy.

[0090] (5) If the difference does not meet the requirement, the flag bit with the abnormal control effect is output, turning to the backup control strategy.

[0091] The design of the threshold is as follows. V.sub.v.sub.r is used to represent a set:

[00009]$V_{v_r}:=\left\{v:v{v}_r\right\}$ [0092] where v.sub.r can be a speed-related value v.sub.r=f(v.sub.1,k, v.sub.2,k) or a constant value, which is determined according to the varying feature of the difference between the actual speed and the target speed in the experimental results. But it should be noted that v.sub.rv.sub.r and mn. This is to avoid the instability of the algorithm resulting from frequent switching between normal and abnormal determination results.

[0093] Then, proceed to function 2 to check the synchronization relationship of all train units. The synchronization of train units is determined by the time index for the current state of all train units on their respective recommended driving curves.

[0094] The function 2 is executed as follows.

[0095] (1) The time index of two train units is calculated: t.sub.i,k=T(s.sub.i,k), in which T(s.sub.i,k) represents the function of converting the train position into the corresponding time index. A desired calculation method is:

[00010]$\begin{array}{cc}T\left(s_{i,k}\right)=\underset{j}{\min }j+\left(s_{i,k}-{\stackrel{}{s}}_{i,j}\right)/\left({\stackrel{}{s}}_{i,j+1}-{\stackrel{}{s}}_{i,j}\right),& s.t.s_{i,k}{\stackrel{}{s}}_{i,j}\end{array}$

[0096] (2) The time index difference between two train units is calculated: t.sub.k=t.sub.1,kt.sub.2,k.

[0097] (3) If the flag bit of a previous cycle shows that the time index difference between two train units is normal, it is determined whether the time index difference t.sub.k=t.sub.1,kt.sub.2,k between both train units in n cycles meets the requirement t.sub.kT.sub.t.

[0098] (4) If the flag bit of a previous cycle shows that the time index difference between two train units is abnormal, it is determined whether the time index difference t.sub.k=t.sub.1,kt.sub.2,k between both train units in m cycles meets the requirement t.sub.kT.sub.t.

[0099] (5) If the difference meets the requirement, the process turns to function 3.

[0100] (6) If the difference does not meet the requirement and the speeds of all train units are greater than v.sub.a, the process turns to function 4.

[0101] A desired design scheme of the threshold value T.sub.t is as follows:

[00011]$T_t:=\left\{t:tt\right\}$ [0102] where t can be a speed-related value t=f(v.sub.1,k, v.sub.2,k) or a constant value, which is determined according to the variation characteristic of the difference between the target speed and the actual speed in the experimental results. But it should be noted that tt and mn. This is to avoid the instability of the algorithm resulting from frequent switching between normal and abnormal determination results.

[0103] According to the output value of the function 2, the process turns to function 3 to calculate the control target sequence in the future. If the process proceeds to function 3, it means that all train units have good tracking control performances and are still in synchronous operation at present. Therefore, all train units can calculate the target sequence according to their respective recommended driving curves, and still ensure the synchronous operation of the virtual coupling at small intervals.

[0104] The function 3 is executed as follows.

[0105] (1) The target speed for the current position where the train is located is calculated: v.sub.i,0|k=V(.sub.i,0|k), s.sub.i,0|k=s.sub.i,k.

[0106] (2) According to a position and a speed at time j, the target position and speed at the next time j+1 in the future are calculated respectively as follows: s.sub.i,j+1|k=s.sub.i,j|k+V.sub.i,j|k and v.sub.i,j+1|k=V(s.sub.i,j+1|k), j=0, 1, . . . , N1.

[0107] (3) Differential calculation is performed on the target speed to calculate the target acceleration .sub.i,j|k, j=0,1, . . . , N1.

[0108] The calculation function V.sub.r(s.sub.i) of the target speed is defined as

[00012]$V_r\left(s_i\right)={\stackrel{}{v}}_{i,j}+$$\begin{array}{cc}\frac{\left(s_i-{\stackrel{}{s}}_{i,j}\right)\left({\stackrel{}{v}}_{i,j+1}-{\stackrel{}{v}}_{i,j}\right)}{{\stackrel{}{s}}_{i,j+1}-{\stackrel{}{s}}_{i,j}},s.t.& {\stackrel{}{s}}_{i,j}{s}_i<{\stackrel{}{s}}_{i,j+1}.\end{array}$

[0109] According to the output value of the function 2, the process turns to function 4 to calculate the control target sequence in the future. If the process turns to function 4, it means that the VCTS is not good in synchronization at this time, and need to be adjusted in combination with the real-time state. The basic idea of the design function 4 is to find the relative displacement by which the following train unit needs to be adjusted, through the time index of the preceding train unit, that is, to convert the time unsynchronization into the position unsynchronization. Thereafter, it is stipulated that the following train unit needs to adjust back by this displacement within a period of time.

[0110] The function 4 is executed as follows.

[0111] (1) When the process proceeds to the function 4 for the first time, the flag bit F.sub.a=1 and the counter c.sub.a=1 are initialized and adjusted. The distance difference s.sub.a=.sub.2,t.sub.1,k.sub.2,t.sub.2,k between the corresponding positions of the time nodes for two train units searched on the recommended speed curve of the following train unit is calculated, and the maximum adjustment speed v.sub.a=T.sub.aa.sub.a{square root over (a.sub.a.sup.2T.sub.a.sup.2/4a.sub.as.sub.a)} is calculated, where T.sub.a and a.sub.a represent the adjustment time designed in advance and the acceleration in the adjustment process, respectively;

[0112] (2) The target speed of the following train unit at the current position of the recommended speed curve of the following train unit is calculated: v.sub.2,0|k=V(s.sub.2,0|k)+v.sub.a(0), s.sub.2,0|k=s.sub.2,k, where v.sub.a(j) represents the adjustment amount related to the synchronization of train units, which is expressed by the following formula:

[00013]$$\begin{gathered} {}_a\left(j\right)= \\ \{\begin{array}{cc}\left(c_a+j\right)a_a,& \left(c_a+j\right){}_a/a_a/\\ {}_a,& {}_a/a_a/\left(c_a+j\right)T_a/-{}_a/a_a/\\ \left(T_a/-\left(c_a+j\right)\right)a_a& \left(c_a+j\right)T_a/-{}_a/a_a/\\ 0& \mathrm{else}\end{array} \end{gathered}$$ [0113] where c.sub.a is the number of cycles between the current cycle and a cycle during which adjustment is performed based on the adjustment amount for a first time and is used to characterize how many cycles have passed since adjustment for the first time by this equation. For example, when adjustment is performed by this method for the first time at a moment k of 10, c.sub.a=2 at the moment k of 11, and c.sub.a=3 at the moment k of 12, and so on.

[0114] (3) According to a position and a speed at time j, the target position and speed at the next time j+1 in the future are calculated: s.sub.i,j+1|k=s.sub.i,j|k+v.sub.i,j|k and v.sub.i,j+1|k=V(s.sub.i,j+1|k)+v.sub.a(j) j=0,1, . . . , N1.

[0115] (4) Differential calculation is performed on the target speed to calculate the target acceleration .sub.i,j|k, j=0,1, . . . , N1.

[0116] This design can eliminate the difference between the time nodes of two train units in a limited time domain T.sub.a. Moreover, the speed adjustment value is gradually increased first, and then is gradually reduced to zero. The changing trend of the target speed sequence is reflected in the subsequent experimental results.

[0117] The process turns to function 5 according to the output values of the function 1 and the function 2. The implementation of the function 5 can follow the existing control method of tracking operation of VCTS, which will not be described in detail in the example of the present disclosure.

[0118] Finally, the process turns to function 6 to prevent the target speed of the train from being too high, resulting in emergency braking during the tracking operation. Its core is realized by determining the relationship between the current speed and the EBI speed of the train.

[0119] The function 6 is executed as follows.

[0120] (1) The difference v.sub.i,kv.sub.i,ebi between the actual speed and the EBI speed is calculated.

[0121] (2) It is determined whether the difference in n cycles meets the requirement of the threshold V.sub.ebi, v.sub.i,ebiv.sub.i,kV.sub.ebi, j=k, k1, . . . , kn+1, where V.sub.ebi:={v: vv.sub.e}, and v.sub.e is the EBI speed margin value designed in advance.

[0122] (4) If the difference meets the requirement, the calculated sequence of the target position, the speed and the acceleration is normally output.

[0123] (5) If the difference does not meet the requirement, the calculated sequence of the target position, the speed and the acceleration is adjusted as follows: calculating

[00014]$v_e=\frac{1}{n}{}_{j=k-n+1}^k{v}_{i,k}-{\stackrel{}{v}}_{i,ebi}+v_e,$

the following operation is performed on all target speeds: v.sub.i,j|k=v.sub.i,j|kk.sub.ev.sub.e, where k.sub.e is a compensation proportional coefficient, and s.sub.i,j+1|k=s.sub.i,j|k+v.sub.i,j|k, j=0,1, . . . , N1.

[0124] (6) Differential calculation is performed on the target speed to calculate the target acceleration .sub.i,j|k, j=0, 1, . . . , N1.

[0125] Thereafter, according to the calculated target sequence, a control command generation module calculates the control command and acts on the train. Because the target generated at this time is a sequence of target states in the future, a model predictive control method can be selected to calculate the control command.

[0126] In FIG. 3, the speed-time graph and distance-time graph of preceding and following train units under the ideal control effect (when following a recommended driving curve) are given. The arrival time interval between both train units at the target station is 1.6 seconds. Further, the following three scenarios are experimented to show some benefits of the present disclosure.

[0127] In a first scenario, the departure of the following train unit is delayed by 5 seconds. If the control method described in the present disclosure is not used, but each train unit still operates according to its own recommended speed curve, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 4. The arrival time interval between the two train units at the target station is 7.2 seconds, which is 5.6 seconds more than that in a benchmark experiment.

[0128] If the system and the method for VCTS control according to the present disclosure are used, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 5. The arrival time interval between the two train units at the target station is 2.8 seconds. When the speed of the following train unit is higher than 8 m/s, the function 4 is enabled. The calculated control target is to add a speed adjustment amount on the basis of the recommended driving curve. The adjustment time reserved in advance is 40 seconds, which means that the speed of the following train unit after 40 seconds should be consistent with the result of the benchmark experiment. The results in FIG. 3 and FIG. 4 also prove this point.

[0129] In a second scenario, the accumulated value of the control errors of the following train unit causes the two train units to be out of synchronization. If the control method described in the present disclosure is not used, but each train unit still operates according to its own recommended speed curve, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 6. The arrival time interval between the two train units at the target station is 4.2 seconds. It can be seen that in an initial traction stage, the following train unit operates slower than the preceding train unit due to the accumulation of control errors. Because both train units are tracking their recommended speed curves, this error has been accumulated to the final arrival stage, resulting in an increase of the arrival time interval.

[0130] If the system and the method for VCTS control according to the present disclosure are used, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 7. The arrival time interval between both train units at the target station is 1.2 seconds.

[0131] In a third scenario, the control error of the following train unit causes the following train unit to operate faster than the preceding train unit, and the difference between the speed of the following train unit and the EBI speed is smaller than that in the benchmark experiment. If the control method described in the present disclosure is not used, but each train unit still operates according to its own recommended speed curve, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 8. The speed of the following train unit is higher than the EBI speed at the 17th second, so that the emergency braking is carried out.

[0132] If the system and the method for VCTS control according to the present disclosure are used, the speed-time relationship and the distance-time relationship between the preceding and following train units are shown in FIG. 9. It can be seen that the difference between the speed of the following train unit and the EBI speed at this time is smaller than the result of the benchmark experiment in FIG. 3, but larger than the result in FIG. 8, and there is no following train unit carrying out emergency braking due to overspeed.

[0133] By comparing the result in FIG. 4 with the result in FIG. 5, the result in FIG. 6 with the result in FIG. 7, and the result in FIG. 8 with the result in FIG. 9, it can be proved that some embodiments of the present disclosure have the following benefits.

[0134] (1) The present disclosure can effectively reduce the effects from the delay of the departure of the following train unit in VCTS, reduce the following distance of the train units in VCTS, and improve the synchronization of arriving at the station.

[0135] (2) The VCTS controlled by the present disclosure can meet the safety protection constraints, avoid triggering emergency braking, and ensure the stability of operation of the virtual coupling.

Embodiment 2

[0136] Embodiment 2 of the present disclosure provides a system for VCTS control, where the system is applied to the method described above, and the system includes a state acquiring module, a first determination module, a first control module, a second determination module, a first target state sequence calculating module, a second state sequence calculating module, and a second control module.

[0137] The state acquiring module is configured to acquire an actual state for a current cycle of each train unit in VCTS.

[0138] The first determination module is configured to determine whether to execute a backup control strategy based on the actual state for the current cycle of each train unit and a target state sequence for a first preset number of cycles before the current cycle to obtain a first determination result; where the backup control strategy includes a control strategy for tracking a recommended driving curve by a first train unit and a control strategy for tracking the i-th train unit by the (i+1)-th train unit, where the value of i is greater than or equal to 1.

[0139] The first control module is configured to execute the backup control strategy to control each train unit, if the first determination result is yes.

[0140] The second determination module is configured to obtain a second determination result, by determining whether synchronization of each train unit in VCTS meets a preset condition, if the first determination result is no.

[0141] The first target state sequence calculating module is configured to calculate the target state sequence for the current cycle of each train unit based on a position according to the actual state for the current cycle of each train unit, if the second determination result is yes.

[0142] The second state sequence calculating module is configured to calculate the target state sequence for the current cycle of each train unit by using a synchronization relationship according to the actual state for the current cycle of each train unit, if the second determination result is no.

[0143] The second control module is configured to control each train unit according to the target state sequence for the current cycle of each train unit, respectively.

[0144] The present disclosure provides an electronic device including a memory, a processor and a computer program stored in the memory and executable on the processor, where the processor implements the method described above when executing the computer program.

[0145] The present disclosure provides a computer-readable storage medium storing a computer program that is executed to implement the method described above.

[0146] The technical solution of the present disclosure based on the above embodiments has the following beneficial effect.

[0147] According to the present disclosure, the recommended driving curve is taken as the control target under normal conditions, and the real-time control target can be adaptively and dynamically adjusted actively on the basis of the recommended driving curve under disturbance conditions. Thereafter, the control command is calculated and output to the train according to the control target. Therefore, the present disclosure can control the VCTS to arrive in the station synchronously under disturbances and prevent the following train unit from triggering emergency braking due to overspeed.

[0148] Embodiments of the present invention are described in a progressive manner, each embodiment focuses on the difference from other embodiments, and the same and similar parts between the embodiments may refer to each other. Since the system disclosed in an embodiment corresponds to the method disclosed in another embodiment, the description is relatively simple, and reference can be made to the method description.

[0149] Although the invention has been illustrated and described in greater detail with reference to the exemplary embodiments, the invention is not limited to the examples disclosed, and further variations can be inferred by a person skilled in the art, without departing from the scope of protection of the invention.

[0150] For the sake of clarity, it is to be understood that the use of a or an throughout this application does not exclude a plurality, and comprising does not exclude other steps or elements.