1.一种亚音速真空管道列车周围气流参数计算方法，其特征在于，包括如下步骤：S1、确定列车的运行速度以及在该速度下的阻塞比；

S2、计算临界阻塞比，并根据临界阻塞比判断列车是否处于壅塞状态；

S3、若列车处于非壅塞状态，则计算列车周围临界区域的气流参数；若列车处于壅塞状态，则分别计算列车前方、列车周围以及列车尾部的气流参数，具体方式为：若管道内部处于非壅塞状态时，列车周围临界区域的气流参数的计算方式包括：S301、将管道内部区域划分为列车前方区域、临界区域以及列车后方区域，其中，列车前方区域和列车后方区域的气流参数保持一致且均为初始气流参数；

S302、根据临界区域与列车前方区域的截面比值与对应速度的关系求解临界区域的马赫数；

S303、根据临界区域马赫数以及临界区域内气流参数与列车前方区域气流参数之间的关系求解临界区域的气流参数；

若管道内部处于壅塞状态，列车前方、列车周围以及列车尾部的气流参数计算方式为：S311、将管道内部按照列车运行方向划分为区域1、区域2、临界区域、区域3、区域4、区域5和区域6，其中，区域1和区域6中的气流参数为初始值；

S312、根据列车运行速度、前驱激波运行速度、尾部激波运行速度以及管道阻塞比求解管道面积与临界区域喉道面积的比值求解壅塞状态下的区域1的马赫数；

S313、根据步骤S312计算的壅塞状态下的列车马赫数求解区域2中管道与列车运行对向的气体速度，并根据该速度计算前驱激波运行马赫数及区域2中的温度参数；

S314、根据区域1与区域2中的气流参数关系，求解区域2中的气流参数；

S315、根据区域2与临界区域的气流参数关系求解临界区域中的气流参数，并根据临界区域与区域3的气流参数关系求解临界区域3的气流参数；

S316、根据尾部激波运行速度求解区域3相对于尾部激波马赫数，并结合区域3与区域4的气流参数关系求解区域4的气流参数；

S317、根据区域3和区域4的参数关系求解区域4的尾部激波马赫数，并根据区域4与区域5的气流参数关系求解区域5的气流参数。

2.根据权利要求1所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述临界阻塞比的计算方式为：

其中，βt为临界阻塞比，Mat为列车速度对应的马赫数且 ut为列车速度，R为理想气体常数，K为气体绝热指数，T1为管道环境温度。

3.根据权利要求2所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述步骤S2中，若βt≥β，则整个管道处于非壅塞状态；若βt<β，管道内部处于壅塞状态，其中，β为管道阻塞比。

4.根据权利要求1所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述步骤S302中临界区域的马赫数的计算方式为：其中，Ma1为列车运行马赫数且Ma1＝Mat，Macr为临界区域气流相对于列车的马赫数且Macr＜1，A1为列车前方区域的气道截面，Acr为临界区域的气道截面。

5.根据权利要求4所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述步骤S303中临界区域的气流参数计算方式为：其中，Pcr为临界区域的管道压力，ρcr为临界区域的空气密度，Tcr为临界区域的温度，P1为初始环境压力、T1为初始环境温度、ρ1为初始环境密度。

6.根据权利要求1所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述步骤S313中壅塞状态下前驱激波运行马赫数计算方式为：其中，Ma1s为前驱激波运行马赫数，u2s为区域2相对于前驱激波的气流速度且u2s＝u1s‑u2r，u1s为区域1相对于前驱激波的气流速度，u2r为区域2中管道与列车运行对向的气体速度；

区域2中管道与列车运行对向的气体速度计算方式为：其中，Ma1为区域1的马赫数，T2为区域2的区域2中的管道环境温度,ut为列车运行速度。

7.根据权利要求6所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述步骤S316中，区域3处尾部激波马赫数的计算方式为：其中，Ma3s为区域3处尾部激波马赫数，T3为区域3中的管道环境温度，us2为列车尾部激波运行速度；

所述步骤S317中，区域4处的尾部激波马赫数计算方式为：其中，Ma4s为区域4处的尾部激波马赫数。

8.根据权利要求7所述的亚音速真空管道列车周围气流参数计算方法，其特征在于，所述区域2中的气流参数表示为：

其中，T2为区域2中的管道环境温度，T1为区域1中的初始环境温度，P2为区域2的管道压力，P1为区域1的初始环境压力，ρ2为管道2的环境密度，ρ1为区域1的初始环境密度；

临界区域的气流参数计算方式为：

其中，Pcr为临界区域的管道压力，ρcr为临界区域的环境密度，Tcr为临界区域的环境温度；

区域3的气流参数计算方式为：

其中，P3为区域3的管道压力、ρ3为区域3的环境密度、T3为区域3的环境温度；

区域4的气流参数计算方式为：

其中，P4为区域4的管道压力、ρ4为区域4的环境密度、T4为区域4的环境温度；

区域5的气流参数计算方式为：

P5＝P4；

其中，P5为区域5的管道压力，T5为区域5的环境温度，T6为区域6的初始环境温度，u5为相对于管道的速度，c6为声速，ρ5为区域5的环境密度，ρ6为区域6的初始环境密度。