京广铁路运行仿真及导航完好性监测系统

资 源 简 介

本系统专门设计用于模拟从北京站到广州南站的长距离铁路线路运行过程。系统首先基于真实的GIS地理数据建立京广线全线三维轨迹模型，涵盖沿途地形起伏、弯道半径及重要车站分布（如石家庄、郑州、武汉、长沙等）。在仿真过程中，系统采用列车运动学和动力学方程实时解算列车的位移、速度和加速度。该项目的核心功能是针对完好性分析进行的GNSS信号环境模拟，能够根据列车所处位置计算实时卫星可见性、空间几何分布（DOP值）以及功率延迟分布，模拟由于桥梁、建筑或地形遮挡导致的信号中断或多路径干扰。系统集成了完好性监测算法，可计算水

详 情 说 明

京广铁路运行线路仿真及其完好性监测分析系统

项目介绍

本系统是一个高度集成的轨道交通仿真平台，专门针对京广铁路（北京站至广州南站）的高速列车运行过程进行全方位模拟。该系统融合了轨道几何建模、列车动力学仿真、全球卫星导航系统（GNSS）信号环境模拟以及接收机自主完好性监测（RAIM）算法。其核心目标是在模拟的长距离运行环境下，评估卫星定位技术在覆盖、精度及完好性方面的表现，特别是在面对地形遮挡和多路径干扰等复杂电磁环境时的安全预警能力。

功能特性

• 全线三维轨迹建模：通过京广线关键站点的地理坐标（纬度、经度、海拔），利用高精度插值算法生成超过2000公里的连续三维轨道模型。
• 列车运行动力学模拟：实现列车从起步加速、常速巡航到终点减速的全过程运动学仿真，实时计算位置、速度及加速度变化。
• 动态GNSS环境评估：根据列车地理位置实时模拟卫星可见性，并计算衡量定位几何强度的DOP值（HDOP/VDOP）。
• 完好性风险分析（RAIM）：集成基于最小二乘残差法的完好性监测算法，实时计算水平（HPL）和垂直（VPL）保护水平。
• 故障注入与告警验证：支持在特定里程点注入伪距故障偏差，测试系统对定位误差的检测能力及告警限值（AL）的触发逻辑。
• 全维可视化分析报告：提供3D运行轨迹图、卫星质量图谱、完好性风险评估曲线及列车运行状态图。

实现逻辑说明

系统的执行流程遵循严谨的数学仿真步骤，确保了物理特性与导航逻辑的一致性：

1. 系统参数配置：系统首先定义地球物理常数、京广线2298公里的全长参数、300km/h的最高时速限制，以及RAIM算法所需的虚警率、漏警率和安全告警限值（HAL/VAL）。
2. 轨道模型生成：采用分段三次埃尔米特（PCHIP）插值技术，将北京、石家庄、郑州、武汉、长沙、广州等核心车站的GIS坐标连接，模拟真实的京广线走向。同时，通过正弦函数叠加随机噪声的方式，模拟从华北平原到南方丘陵的海拔起伏。
3. 运动状态解算：基于时间步进法，根据列车所处路段（启动区、制动区或巡航区）动态调整加速度，计算列车沿轨道的累积里程，并实时映射回对应的经纬度坐标。
4. GNSS观测模拟：预设24颗卫星组成的星座分布，根据当前列车坐标计算观测矩阵。系统加入了遮挡逻辑：在模拟的高架桥或隧道等特定里程段，增加信号丢失概率，以验证极端环境下的定位可用性。
5. 完好性指标计算：系统利用观测矩阵的投影性质计算各卫星的故障影响斜率（Slope）。通过卡方分布的分位数确定检测门限，从而推导出当前时刻的保护水平（PL值）。
6. 统计与可视化：对比计算出的保护水平与设定的告警限值，统计全行程的系统可用性百分比，并绘制多维图表展示仿真结果。

关键算法与实现细节分析

• 三维插值算法：为了使模拟线路平滑且符合铁路弯道特征，系统使用了PCHIP插值而非简单的线性插值，这保证了轨道在各站点的坐标衔接具有一阶连续性。
• GNSS几何分析：算法通过构建由卫星方位角和仰角组成的几何矩阵G，并计算其伪逆矩阵，进而得出DOP值。这是评估卫星空间分布对定位精度影响的核心手段。
• RAIM保护水平估算：采用了经典的“斜率法”（Slope Method）。算法提取几何矩阵中对水平和垂直分量贡献最大的项，结合预设的统计概率门限，计算出在当前卫星分布下，能够以极高概率包络定位误差的最小半径。
• 自适应分布函数：由于仿真可能在缺乏统计工具箱的环境下运行，系统内部实现了逆卡方分布（chi2inv）和逆误差函数（erfcinv）的数值逼近算法，确保了完好性判决门限计算的准确性。
• 故障检测逻辑：系统在第800及1500仿真采样点人为注入30米的伪距偏差。当受扰后的定位误差超过计算出的HPL时，系统会自动标记完好性故障。

使用方法

1. 环境准备：确保已安装MATLAB（推荐版本R2018b及以上）。
2. 运行仿真：直接在MATLAB命令行窗口输入该主函数名称并回车。
3. 交互查看：

* 观察弹出的Figure窗口，左上图可旋转查看京广线全线3维运行轨迹。 * 右上图分析沿途经过“隧道/桥梁”（卫星数下降区域）时对DOP值的影响。 * 左下图重点关注黑点（实际误差）是否突破红线（保护水平）及红虚线（告警限值）。

1. 读取报告：查看MATLAB命令行输出的“完好性评估报告”，获取总仿真里程、平均卫星数、风险点数量及系统整体可用性指标。

系统要求

• 软件环境：MATLAB 2016a或更高版本。
• 硬件要求：建议内存4GB以上，系统已针对性能优化，支持长距离快速仿真。
• 工具箱依赖：系统自带完整的数学函数实现，不强制依赖Statistics and Machine Learning Toolbox，具备良好的兼容性。