1. Introduction
As the core subsystem of the train operation control system, the train positioning system plays an important role in the safe operation of trains. At present, with the rapid development of China’s railway industry, the traditional train positioning methods such as odometers, track circuits, and transponders have gradually exposed various disadvantages, which are difficult to meet the needs of the new train control system. At the same time, with the rapid development of railway modernization and intelligence, it will become a trend to reduce trackside equipment and vehicle-mounted trackside equipment [
1]. Therefore, a new generation of low-cost vehicle-mounted train positioning systems is urgently needed. In June 2020, china successfully launched the last global networking satellite of beidou-3, marking the arrival of the Beidou global networking era. The train positioning system based on satellite navigation will become the development direction of the next generation of train positioning [
2]. Different from the traditional vehicle integrated positioning system based on satellite navigation, the train positioning system has higher requirements for reliability and safety. Therefore, the safety risk analysis is particularly important, and it has important guiding significance for the later system safety design and system safety protection.
At present, a large number of studies on security risk analysis based on fault tree analysis (FTA) and analytic hierarchy process (AHP) have been conducted in China and abroad. The FTA method can trace the source of the causes and influence of the target risk events by establishing a model, in which the failure probability of the basic events is accurately expressed, but in fact, due to the lack of sufficient statistical calculation of events, such accuracy is not accurate. In view of this problem, An et al. [
3] proposed a security risk analysis method based on FTA and fuzzy system combined with fuzzy logic. Mamdikar et al. [
4] proposed a framework for the dynamic reliability assessment using the Fault Tree and the Dynamic Bayesian Network. Aiming at the problem of FTA cannot judge the overall security risk level of the target risk event, Liu et al. [
5] introduced AHP into the security risk assessment, and took the basic cause event information calculated by FTA as the input data for establishing the AHP model, thereby improved the accuracy of risk assessment. In view of the problems of uncertainty and inaccuracy caused by observation and statistics of relevant data in AHP modeling, Bakır et al. [
6] proposed an integrated Fuzzy Analytical Hierarchy Process (F-AHP) and Fuzzy Measurement Alternatives and Ranking according to Compromise Solution (F-MARCOS) approach. In order to effectively solve the problem of consistency of risk judgment matrix in AHP, su et al. [
7] improved AHP and proposed to use neural network technology to build a neural network model corresponding to risk judgment matrix for objective risk assessment. Ma et al. [
8] combined analytic hierarchy process (AHP) and CRITIC to assign weights, and entered the cloud model on this basis, providing a reference for evaluating the train operation control system based communication.
According to the above research contents and the particularity of research objects, this paper introduces the idea of interval number and improves the AHP method, combines FTA and IAHP to propose a safety risk assessment method for a train positioning system based on FTA-IAHP on the basis of the above literature research. It is used to solve the problems when applying one of the above methods for safety risk analysis, so that the analysis results can not only reflect the information of subjective evaluators, but also reasonably reflect the objective facts. This method is applicable to the risk analysis of system insecurity caused by internal fault factors of parts and external uncontrollable environmental factors.
In this paper, combined with the self-developed train positioning system based on satellite navigation and the safety evaluation system of train positioning system can be set up.
Section 2 focuses on presenting the overall structure of train positioning system based on satellite navigation.
Section 3 introduces the improved safety risk analysis method in detail. The safety evaluation system and safety risk analysis are established for the research object in
Section 4, and the effectiveness and availability of the method are verifified through the analyses on actual examples. Finally,
Section 5 concludes this paper with a summary.
2. Train Positioning System Based on Satellite Navigation
Carry out safety risk analysis on the train positioning system of satellite navigation, and the system diagram is shown in
Figure 1.
The whole hardware system is divided into four parts: multi-channel security power subsystem, sensor information acquisition subsystem, data fusion processing subsystem, and 3G wireless communication subsystem. The multi-channel safety power supply subsystem mainly realizes the power supply function for the whole system, and at the same time meets the functions of power short circuit, reverse connection, and overvoltage protection. The sensor information acquisition system mainly realizes the information acquisition and synchronous processing of inertial measurement unit (IMU), satellite navigation system (Beidou, GPS, GALNASS), and the vehicle odometer (ODO). The data fusion processing subsystem mainly completes the combined positioning algorithm processing and two out of two functions of two or three sets of sensor data collected from the bus. Considering the inertial navigation solution and algorithm processing, high-speed DSP is used here. The communication board mainly completing the final analysis of data and to complete the wireless communication function. The whole system adopts the double two out of two security design, which can achieve high-precision and fast positioning and velocity measurement while meeting reliability and security.
Figure 2 and
Figure 3 are the front view of chassis and physical drawing for the independently developed train positioning system, respectively. According to the system risk assessment method proposed in this paper, the security risk of the system will be evaluated.
3. Safety Risk Assessment Method Based on FTA-IAHP
Any reasonable and comprehensive analysis of complex systems is based on information and data of relevant factor, and the judgment matrix is an important form of information. To judge the importance of each factor, we can represent each factor as a numerical result and present the results in matrix form. This paper uses two different methods to construct the judgment matrix: one is FTA basic events based on (), and the other one is expert scoring based on IAHP. The analysis results of the two judgment matrices are synthesized, a new risk analysis method is proposed.
Figure 4 is the flow chart of the risk evaluation method proposed in this paper. The whole method can be divided into two parts: FTA analysis and the improved IAHP method. Firstly, the target risk events of the system are determined, and the FTA model is established to determine the basic cause events. Secondly, the FTA judgment matrix is established by comparing the structural importance of basic cause events; IAHP method is used to further analyze the expert score and establish the interval judgment matrix; at the same time, this paper improves the weight of the traditional IAHP, and proposes new combined weight coefficient solution method by integrating the weight of FTA and IAHP judgment matrix. The comprehensive evaluation of the train positioning system’s target risk events is completed by comparing with the safety risk assessment standards of railway.
The detailed steps are as follows:
3.1. FTA Analysis Method
In order to estimate the security of the militia launch control system(MLCA), Watson first proposed the FTA analysis method in 1961. Later, the FTA method was improved quantitatively by using computers. At present, this method is widely used, such as aviation, nuclear engineering, and safety management. To adopt FTA method, we first need to analyze the system to determine the target risk events, etc. The specific steps are shown in
Table 1.
3.1.1. Building Fault Tree and Calculating Structural Importance
After establishing FTA model according to the specific steps in
Table 1, qualitative or quantitative analysis can be carried out. Structural importance of basic cause events is the most direct indicator of the impact on target risk events [
9]. The structure function of fault tree can be defined as
where,
stands for boolean variable for the state of the base-cause events;
n represents the quantity of all basic cause events obtained from FTA analysis; 0 indicates that the target event did not occur; 1 represents the occurrence of the target risk event. The structural importance of the
ith basic cause event can be expressed as:
Among them,
(
) represents the structural importance of the
ith basic cause event;
indicates that after adding
ith basic cause event to
combination based on cause event, the combination changes from a non-cut set to a cut set, which can be expressed as:
From the above formula, we can get the quantitative influence degree of basic cause events on target risk events.
3.1.2. Establish the Safety Risk Hierarchy Model
The top event of the fault tree is taken as the target layer T (the highest layer) of the hierarchical structure model, and the basic cause event at the bottom of the fault tree is taken as the indicator layer E (the lowest layer) of the hierarchical structure. It is classified and summarized according to certain logic rules to form the criterion layer B (the middle layer).
3.1.3. Judgment Matrix Based on FTA
Formula (2) can provide the
(
) of each basic event, and it is easy to obtain the least common multiple (LCM) of denominator (common characteristic) [
10]. The judgment factor of basic events is calculated as follows:
The judgment factors of basic events reflect the effectiveness of top-level events. Judgment matrix can be constructed by the ratio of two factors. Because each criterion layer contains some index factors, this paper summarizes the index factors of each criterion layer, not the judgment factors. Since each element of the judgment matrix must be an integer, the numbers in each matrix are rounded to the nearest integer. The formula of the criterion layer used to create the judgment matrix E is as follows:
Among them,
m and
n is the number of indicator factors for each criterion,
and
j are the subscripts of the elements of the matrix, which respectively represents
i Line and section
j Column. Similarly, the method for creating the judgment matrix
B of the index layer for each standard layer in the article is given in Formula (6).
According to the (i) of each basic event, it is easy to get the lowest common multiple and the judgment factors of criterion layer and index layer.
3.1.4. Consistency Test and Vector Calculation of Judgment Matrix
The judgment matrix is the basis for solving the weight. A chaotic and untenable matrix can lead to wrong results. Therefore, in order to improve the reliability of the judgment matrix, it is required that the judgment matrix be generally consistent. In consideration of the above, the consistency check of the judgment matrix shall be carried out by the following methods.
Then the consistency test formula is:
If
CR < 0.1 is satisfied, the judgment matrix is considered to have good consistency; otherwise, the matrix needs to be readjusted until the consistency reaches the standard. Where
n is the order of the judgment matrix,
is the maximum eigenvalue of each layer of the matrix, and
RI can be queried from
Table 2. For example, if the judgment matrix is a third-order matrix, the corresponding
RI value is 0.58.
The maximum eigenvalue of the judgment matrix and the corresponding eigenvector can be calculated by MATLAB software (version: R2021b, creator: Jack little and clever moler, USA), the steps are as follows:
3.2. Interval Analytic Hierarchy Process (IAHP)
Analytic hierarchy process (AHP) is a multiple criteria decision-making method. In this method, the measurement theory is introduced into the quantification of subjective judgment of experts, and the combination of qualitative analysis and quantitative analysis is realized. However, due to the incompleteness and uncertainty of information, uncertain subjective judgment often appears in pairwise comparison in practice. Therefore, it is not reasonable to use the point value of traditional AHP method to describe the subjective judgment value. At this time, interval mathematics is most suitable for accurate description of subjective judgment value [
11]. Therefore, IAHP (interval analytic hierarchy process) is proposed.
The IAHP method is based on the traditional AHP method and integrated with interval mathematics. The steps are as follows: the interval number is used to replace the point value to describe the experts’ assessment of the relative importance of various safety factors, and to construct the judgment matrix by pairwise comparison; after checking the consistency of interval, the weight vector of interval number is calculated; the interval comprehensive weight is obtained by calculating the judgment matrix and interval weight vector. Finally, it is sorted, and the specific implementation block diagram is shown in
Figure 5. This method can well describe the uncertainty of judgment, so it can reflect the fuzziness of experts’ subjective judgment on the importance of events, and reduce the subjective influence of evaluation to a certain extent. Therefore, IAHP is used to describe the experts’ assessment of the relative importance of various safety factors. The following is a detailed introduction to the important steps:
3.2.1. Composition of Interval Judgment Matrix
If a real number satisfies
, we call
an interval number. If interval number is used to express the relative importance of each factor in pairwise comparison, the judgment matrix of interval sequence is formed [
12].
Assuming that there are n safety factors in a certain level of the index system, experts use 1~9 scale method to evaluate the relative importance of safety factors.
According to the reciprocal 1–9 scale table [
13], the qualitative description of expert pairwise comparison is scalar, as shown in
Table 3, in which the language description degree of level 2, 4, 6, and 8 is the intermediate value of adjacent judgments in the table.
If the index judgment matrix and = [,] is satisfied, it indicates the importance of evaluation index i relative to evaluation index j. So, the comparison result of is the reciprocal of to represents the result of comparing evaluation index i with itself, we can see that it is the same as that of itself [1,1].
Then, the expert opinions are synthesized, and the interval number judgment matrix is obtained.
3.2.2. Consistency Test of the Judgment Matrix
If there are more than two comparisons in the process of the expert pairwise comparison and judgment, there will be inconsistent judgments. When the amount of the pairwise comparison is large, this kind of inconsistency is more likely to occur, and sometimes totally inconsistent judgment results may be obtained. At this time, the reliability of the judgment information will be reduced. If the result is further calculated, it may lead to wrong results. In view of the above situation, it is necessary to check the consistency of the judgment matrix [
14].
Set up
, in which
, and take:
If and are satisfied, in this case, we consider that the judgment matrix has good consistency; on the contrary, if 1 or 1 is satisfied, it is considered that the consistency of judgment matrix is poor. At this time, we need to feed back to the experts to rejudge until the satisfactory consistency is obtained.
3.2.3. Solving the Weight o Judgment Matrix
Since the 1990s, the research on the weight calculation algorithm of interval number judgment matrix has been carried out successively, including iterative method, random simulation method, interval eigenvalue method [
15], etc. In this paper, the eigenvalue method is used to solve the weight, and the specific steps are as follows:
Step 1: Assuming represents the interval 1–9 scale reciprocal judgment matrix, where = [,], and assuming that = and , A = [,] is available.
Step 2: Using the eigenvalue method to calculate the point value matrix
and
, the weight vectors is obtained as follows:
At this time, the weight vector of interval number is:
where,
= [
,
],
i = 1, 2,
L,
n.
3.3. IAHP-FTA Portfolio Evaluation
3.3.1. Solution of Combination Weight Coefficient
This paper studies the weight of FTA
and IAHP
. The former reflects the influence degree of basic factors on top-level events, while the latter describes experts’ assessment of relative importance of various safety factors. Therefore, it is feasible to combine the two risk assessment methods by using the following formula. Through this solution method of the combined weight coefficient, it can not only integrate the quantitative reflection of the impact of the event importance on the target event in the FTA method, but also weaken the subjectivity of the system evaluation. It is a new reasonable and accurate method of weight calculation.
Among them, α and are weighting coefficients, is the sum of scale factor and consistency factor of FTA judgment matrix, and is the sum of scale factor and consistency factor of interval consistency approximation matrix.
Their calculation method is as follows:
3.3.2. Weighted Coefficient Solution
The consistency factor can be solved according to Formula (14), wherein the solution can be obtained by Formulas (7) and (8):
Supposing there are two sets of number
and
, the square root of two groups of average square error is expressed by RMS (risk matrix proportion) as follows:
Then, the following formula can be used to calculate the scale factor.
where
and
are the maximum
RMS value and the minimum
RMS value corresponding to the scale, and are the values corresponding to the 1–90 scale and the 1–9 scale respectively.
Calculate the
RMS value under different scales, and the scale factor under different scales can be obtained according to Equation (16). The different scale factors are shown in
Table 4. The scaling factor of the intermediate scale can be obtained by interpolation.
3.3.3. Final Overall Assessment
Referring to EN50126 standard and the current research results, when formulating the scoring rules, the safety risk assessment standard of the train positioning system is divided into four levels due to the involvement of the train signal safety system, and the four grade indexes are quantified [
16]. The scoring range is as shown in the
Table 5.
Combined with the fault of the train positioning system and the experience of relevant experts, relevant experts are organized to evaluate and score the safety risk of each factor of the index layer. The evaluation standard is scored according to the score value in
Table 5. Then the safety assessment score of the criterion layer is calculated according to Equation (17).
Among them, refers to the final weight vector of the ith factor in the index layer obtained according to the combined weight method. Through this formula, the safety risk of the factors on the upper layer can be calculated. According to this law, it can be calculated to the target level. Finally, the comprehensive safety risk assessment results of the target layer factors of the train positioning system can be obtained.
4. Safety Risk Assessment of the Train Positioning System Based on Satellite Navigation (Example Verification)
4.1. Establishment and Analysis of the System FTA Model
According to the self-developed train positioning system based on satellite navigation as the research object, the FTA model is established and analyzed by using the above method. Taking the fault of the train positioning system as the target risk event [
17], the FTA model (taking the two-out-of-two system as an example) is established as shown in
Figure 6.
Wherein the structural importance of each basic cause event is obtained according to Equation (2). As shown in
Table 6.
4.2. Establishment of Hierarchical Model for Safety Risk Analysis of Train Positioning System
Taking the train positioning system based on satellite navigation as the target level factor, combined with the analysis results of FTA model, the safety evaluation hierarchy system of train positioning system based on satellite navigation is constructed, which is composed of target layer, criterion layer and index layer. The basic cause events in FTA model are divided into four categories, corresponding to the factors
in the criteria layer of the positioning system model respectively. The correspondence of each factor is shown in
Table 7.
4.3. Establishment of Judgment Matrix and Weight Solution
4.3.1. The Judgment Matrix Based on FTA and Weight Solution
Taking the fault of positioning unit as the target layer factor, combining with the basic cause events obtained by fault tree analysis, it is divided into four categories, corresponding to quasi lateral layer factors respectively. The basic cause events are taken as the index layer factors. According to
Table 6, the
of each basic event can be provided, and the least common multiple (LCM) of denominator (common characteristic) is easily obtained is 32. The judgment factor of each basic cause event is calculated by Formula (4).
Then according to Formulas (5) and (6), the judgment matrix of each layer and the judgment matrix of target layer and the weight vector of each layer factor are constructed respectively. As shown in
Table 8. This paper uses MATLAB software to carry out consistency test according to Formulas (7) and (8), and the judgment matrix has good consistency,
. According to the above method, the weights of the indexes of each layer can be obtained in sequence.
The judgment matrix of each index layer under the criterion layer can be obtained according to Equation (6). The calculation process is described by the index of criterion layer
B1, as shown in
Table 9.
The consistency inspection index
CR =
meets the requirements of
CR < 0.1, that is to say, it meets the consistency test.
B2,
B3 and
B4 weight vectors are obtained by the same method, and the final results are listed in
Table 10.
4.3.2. The Judgment Matrix Based on IAHP Expert Rating
Firstly, 10 experienced experts are invited to compare the safety status of 25 underlying index factors of the evaluation system. Interval judgment matrix is established by layers one by one to check whether the matrix meets the consistency, and the weight vector of interval number is calculated according to Formulas (11) and (12). The calculation process is illustrated by the criterion layer index
B3. Through expert discussion on the index factors E18, E19, E20 and E21 of
B3 criterion layer, the expert score judgment matrix is obtained according to Formula (9):
The judgment matrix of other criterion layers can be obtained in turn.
Finally, the weight of expert on the evaluation of the underlying influencing factors is obtained. It is shown in the table below. The calculation is as follows: the interval number judgment matrix shown in
B3 is divided into two matrices, which are
According to Equation (11), the eigenvector is
According to Formula (10):
= 0.928,
, the interval weight vector calculated by substituting into Equation (12) is:
. The consistency test index
CR =
meets the requirements of
CR < 0.1, that is to say, it meets the consistency test. According to the same method, B1, B2, B4 weight vectors are obtained, and the final results are listed in
Table 10.
4.4. Overall Safety Risk Assessment of Positioning System
According to the weight of the judgment matrix based on FTA and IAHP experts evaluation, combined with the combination weight method proposed in this paper, according to Formulas (14) and (16), matlab is used to get the result
and
. Finally, the final weight factor calculated according to Formula (13) is shown in
Table 10.
According to Formula (17), the comprehensive safety evaluation score of the whole positioning system is 4.95, which meets the second level of the safety risk evaluation standard of the train positioning system. According to reference [
18], most scholars at home and abroad believe that the security goal of GNSS based positioning subsystem is SIL2. However, others say that higher levels of security, such as achieving the target SIL3 or SIL4, may limit the implementation of high precision. Therefore, it can be concluded that the system can meet the requirements of the current train positioning system SIL2 level safety requirements.
5. Conclusions
This paper applies the proposed security risk analysis method based on FTA-IAHP to conduct comprehensive security risk analysis on the independently developed positioning system hardware. It is of great significance to discover the potential risks of the system, optimize the system structure, and ensure the high security of the system. By constructing a hierarchical security risk security model, two judgment matrixes are constructed based on the structural importance of FTA basic events and expert evaluation based on IAHP model, and the difference between the two judgment matrixes is adjusted by weighting factor. The method combines the weights of the two judgment matrices to determine the priorities of the risks that lead to the system failure, and reflect the security risk status of the entire system. The results show that the hardware random safety integrity level of the train positioning system based on GNSS designed in this paper is SIL2, which meets the integrity risk objective required by EN56102. Through the safety risk assessment, we can know the basic causes of the hazard sources that cause the failure of the positioning system, and make it clear that the components that have a great impact on the safety risk status of the positioning system are the power supply board and the communication processing CPU. The weak links of the system structure are found, which provides a basis for the structural optimization of the system. The authors believes that the hazard identification methods in the actual application of the positioning system can also be analyzed in the future research through the methods proposed in this paper, to further provide scientific basis for the follow-up safety risk control. The limitation of this study is that the impact of safety related application requirements on risk assessment parameters in the train control system was not fully considered. In the future, in combination with this study, further consider the impact of the location-based service application requirements on the positioning unit in the train control system, study the safety integrity monitoring algorithm, etc., to further reduce the system safety risk, and better meet the safety requirements of the high-speed railway for the positioning system.
Author Contributions
Conceptualization, Y.Y. and G.C.; methodology, Y.Y. and D.W.; validation, Y.Y. and D.W.; formal analysis, Y.Y. and G.C.; investigation, D.W. and G.C.; data curation, D.W. and G.C.; writing—original draft preparation, Y.Y.; writing—review and editing, G.C. and Y.Y.; supervision, G.C. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported in part by the Youth Program of Science and Technology Foundation of Gansu Province (21JR7RA323), The National Natural Science Foundation of China (62141303), The Science and Technology Guidance Plan of Gansu Province (2020-61-14) and the Gansu Provincial Science and Technology Major Project of China (21ZD4WA018).
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
All datasets generated for this study are included in the article.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Mo, Z. Technical development direction and Prospect of intelligent railway train control system. Railw. Commun. Signal 2022, 58, 1–7. [Google Scholar]
- Wang, J.; Zhou, Z.; Jiang, W.; Cai, B.; Pan, P.F. High precision real-time positioning method of train based on GPS/BDS joint solution. J. Transp. Eng. 2021, 21, 286–296. [Google Scholar]
- An, H.; Huang, A.; An, M.; Fan, L.; Jin, M. Comprehensive risk assessment of falling from height in construction based on fuzzy fault tree. Sci. Technol. Eng. 2022, 22, 8568–8576. [Google Scholar]
- Rao, M.M.; Vinay, K.; Pooja, S. Dynamic reliability analysis framework using fault tree and dynamic Bayesian network: A case study of NPP. Nucl. Eng. Technol. 2022, 54, 1213–1220. [Google Scholar]
- Liu, Y.; Cao, G.; Yang, H.; Feng, J. Research on safety risk assessment of STP wireless shunting locomotive signal and monitoring system. Railw. Transp. Econ. 2018, 40, 91–97. [Google Scholar]
- Bakır, M.; Atalık, Ö. Application of Fuzzy AHP and Fuzzy MARCOS Approach for the Evaluation of E-Service Quality in the Airline Industry. Decis. Mak. Appl. Manag. Eng. 2021, 4, 127–152. [Google Scholar] [CrossRef]
- Su, C.; Li, Y.; Mao, W.; Hu, S. Information Network Risk Assessment Based on AHP and Neural Network. In Proceedings of the 2018 10th International Conference on Communication Software and Networks (ICCSN), Chengdu, China, 6–9 July 2018; pp. 227–231. [Google Scholar]
- Ma, J. Reliability Analysis and Evaluation of Urban Rail Transit CBTC System; Lanzhou Jiaotong University: Lanzhou, China, 2020. [Google Scholar]
- Feng, F.; Yan, M. Risk assessment of Railway Express cargo loss based on BDD. J. Beijing Jiaotong Univ. 2018, 42, 48–54 + 74. [Google Scholar]
- Chen, S.; Xiong, S.; He, C.; Lan, B. Safety risk analysis of subway deep foundation pit collapse accident. J. Saf. Environ. 2020, 20, 52–58. [Google Scholar]
- Cong, L.; Xiao, Z.; Chen, Y.; Guo, Z. State evaluation method of a certain type of anti-ship missile based on improved cloud model. J. Weapons Equip. Eng. 2021, 42, 85–93. [Google Scholar]
- Zhang, Y.; Zhang, F.; Zhu, H.; Guo, P. An Optimization-Evaluation Agricultural Water Planning Approach Based on Interval Linear Fractional Bi-Level Programming and IAHP-TOPSIS. Water 2019, 11, 1094. [Google Scholar] [CrossRef]
- Ding, J.; Shen, S.; Li, C. Evaluation method of intelligent distribution network planning scheme based on interval preference and MADM. Power Syst. Prot. Control. 2018, 46, 18–27. [Google Scholar]
- Huang, R.; Tian, Z.; Lv, Y. Research on consistency of reciprocal judgment matrix of interval rough numbers. Fuzzy Syst. Math. 2019, 33, 124–133. [Google Scholar]
- Duan, L.; Liu, Y.; Liu, S.; Tang, M. Safety evaluation of reciprocating compressor unit based on IEM fuzzy analysis theory. Ind. Saf. Environ. Prot. 2019, 45, 31–35 + 39. [Google Scholar]
- Yao, J. Research on Risk Assessment of CTCS-3 Train Control System; Lanzhou Jiaotong University: Lanzhou, China, 2020. [Google Scholar]
- Wu, B. Research on Evaluation Method of Train Autonomous Positioning Performance Based on GNSS/INS; Beijing Jiaotong University: Beijing, China, 2021. [Google Scholar]
- Beugin, J.; Legrand, C.; Marais, J.; Berbineau, M.; Miloudi, E.K.E. Safety Appraisal of GNSS-Based Localization Systems Used in Train Spacing Control. IEEE Access 2018, 6, 9898–9916. [Google Scholar] [CrossRef]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).