[发明专利]一种地铁潮流计算中节点电导矩阵的快速计算方法

权利要求书

1.一种地铁潮流计算中节点电导矩阵的快速计算方法，其特征在于，包括：

S100，提供各个元件相对上行列车起点的距离和各个元件的数量，以及电阻参数，其中，所述各个元件包括牵引变电所、上行列车以及下行列车，所述电阻参数包括接触网单位电阻、阻性钢轨单位电阻和轨道对地的泄露单位电阻；

S200，按照相对上行列车起点的距离依次分别对牵引变电所进行排序编号，将列车线路的起点和终点之间的区间划分为多个供电区间，并对供电区间进行编号和长度计算；

S300，构建节点电导矩阵，所述节点电导矩阵包括牵引变电所与牵引变电所之间的节点子电导矩阵Y_SS、牵引变电所与上行列车之间的节点子电导矩阵Y_SU、牵引变电所与下行列车之间的节点子电导矩阵Y_SD、上行列车与牵引变电所之间的节点子电导矩阵Y_US、上行列车与上行列车之间的节点子电导矩阵Y_UU、下行列车与牵引变电所之间的节点子电导矩阵Y_DS、下行列车与下行列车之间的节点子电导矩阵Y_DD，并将Y_SS、Y_SU、Ys_D、Y_US、Y_UU、Y_DS以及Y_DD置零；

S400，按照无列车运行时计算Y_SS；

S500，按照相对上行列车起点的距离依次分别对列车进行排序编号，其中，上行列车的编号与牵引变电所的编号连续，下行列车的编号与上行列车的编号连续；

S600，清零供电区间内上行列车数向量、上行供电区间处理标志向量、上行列车左端元件索引编号以及上行列车右端元件索引编号，设置上行列车循环计数次数的初始值；并根据上行列车循环计数的次数与上行列车辆数的关系，循环计算上行列车左端元件索引编号以及上行列车右端元件索引编号；以及对Y_SS进行修正或修正Y_SS、计算Y_SU以及Y_UU或计算Y_UU或计算Y_US；

S700，清零供电区间内下行列车数向量、下行供电区间处理标志向量、下行列车左端元件索引编号以及下行列车右端元件索引编号，设置下行列车循环计数次数的初始值；并根据下行列车循环计数的次数与下行列车辆数的关系，循环计算下行列车左端元件索引编号以及下行列车右端元件索引编号；以及对Y_SS进行修正、Y_SD以及Y_DD或计算Y_DD或计算Y_DS。

2.根据权利要求1所述的地铁潮流计算中节点电导矩阵的快速计算方法，其特征在于，步骤S400包括，

通过

计算Y_SS，其中，i₁和j₁为节点编号，j₁∈i₁表示直接与i₁相连的节点，r₁为接触网单位电阻，L_t为第t个供电区间的长度，t为供电区间编号，n为牵引变电所总数，r₂为走行钢轨单位电阻，r₃为钢轨对地泄露单位电阻。

3.根据权利要求2所述的地铁潮流计算中节点电导矩阵的快速计算方法，其特征在于，对Y_SS进行修正包括：当上行列车循环计数的次数与上行列车辆数的关系满足i＝m+1，上行列车供电区间的处理标志向量dealup(j)等于0，且上行供电区间的索引编号满足j≠1且j≠n+1时，或当列车循环计数的次数与下行列车辆数的关系满足i＝k+1，下行列车供电区间的处理标志向量dealdown(j)等于0，且供电区间的索引编号满足j≠1且j≠n+1时，

通过对Y_SS进行修正，其中，i为上行列车或下行列车列车循环计数的次数，m为上行列车数量，k为上行列车数量，n为牵引变电所总数，n+1为供电区间的总数，L_j为上行列车或下行列车所属第j个供电区间的长度。