一种卫星导航技术领域的用于GPS的基于影响因素特征的脆弱性监测系统及方法，该系统包括：控制计算模块、评估计算模块、空间段脆弱性仿真模块、软件接收机模块和监测站模块。本发明解决了现有卫星导航系统脆弱性影响因素的监测过程复杂且缺乏定量判断标准的技术问题，在不引入可见的计算误差的前提下，大大提高脆弱性影响因素的监测效率，使得全球导航卫星系统的脆弱性监测与判断有系统、快速、准确的实现途径。
1.一种用于GPS的基于影响因素特征的脆弱性监测方法，其特征在于，包括以下步骤：
步骤一、脆弱性监测系统启动，在监测设备上设定需要监测的目标系统信号，得到所有卫星信号的集合M＝{x1(t)，x2(t)，x3(t)，......，xN(t)}，并且得到初始集合的大小N，其中：xi(t)表示在选定频点上不同类型的信号，N表示信号类型总数；
步骤二、根据目标频段的系统信息，从脆弱性监测系统处获得GPS的卫星导航信号的信息列表；
步骤三、对实时导航电文进行电文一致性评估：监测来自脆弱性监测系统实测的实时导航电文，根据其中的轨道信息，进行当前卫星坐标的推算，将得到的卫星坐标推算结果与来自IGS服务器的精确预报星历中的卫星坐标进行对比，得到每颗卫星的当前推算坐标与预测坐标的偏差值 其中：x0、y0、z0分别表示根据实测导航电文推算出的卫星坐标，x、y、z表示从IGS服务器单元精确预报星历中得到的卫星坐标，当偏差值未超过设定阈值则判定导航电文正确并执行步骤四，否则判定导航电文异常并进一步判定为发生了系统电文异常或者产生式欺骗干扰；因为产生式欺骗干扰会先导致接收机信号质量的下降，以至接收机失锁才能完成欺骗，故要结合步骤四的信号质量评估来判断到底是系统电文异常还是产生式欺骗干扰；
步骤四、信号质量评估：监测接收到的最小信号功率是否达到或超过ICD指定数值，当小于等于ICD指定数值则认为出现产生式欺骗，否则认为是系统电文异常；
步骤五、采用IONEX和Klobuchar两个模型结合的方式进行基于电离层异常的大气环境评估，具体为：设定Klobuchar延迟的2倍作为延迟阈值，当IONEX延迟值延迟超出延迟阈值时，即当 时，则判定电离层异常发生，否则认为电离层正常，其中：ΔIONEX表示根据IONEX模型实时计算的电离层延迟，ΔKlobuchar表示根据Klobuchar模型实时计算的电离层延迟；
步骤六、完好性监测评估：软件接收机模块通过卫星间观测数据，利用SAIM算法，从系统内判断卫星故障状况，并得到当前异常工作的卫星信息，同时利用RAIM算法判断卫星通道故障情况，并得到当前测量值错误的通道信息；
步骤七、接收机定位误差监测：软件接收机模块得到当前定位误差并与设定阈值进行比较，当当前定位误差未超过设定阈值，则判定定位误差处于可接受范围判断无脆弱性事件，否则判定定位误差过大，判断肯定存在脆弱性影响因素，最后软件接收机模块根据比较结果生成定位误差结果；
步骤八、将步骤三到步骤七得到的导航电文、信号质量、大气环境、完好性监测及定位误差结果，输入至脆弱性监测系统，得到系统脆弱性信息；
所述的系统包括：控制计算模块、评估计算模块、空间段脆弱性仿真模块、软件接收机模块和监测站模块，其中：控制计算模块与空间段脆弱性仿真模块相连并传输脆弱性因素生成参数，空间段脆弱性仿真模块分别向监测站模块传输电子浓度、对流层模型参数、电离层闪烁参数和/或电离层模型参数，并向软件接收机模块传输星间链路测距值，评估计算模块分别与空间段脆弱性仿真模块、监测站模块和控制计算模块连接，利用空间段脆弱性仿真模块和监测站模块生成的各种脆弱性因素模拟数据计算评估结果并将结果反馈给控制计算模块。
2.根据权利要求1所述的方法，其特征是，所述的空间段脆弱性仿真模块包括：GPS星间链路信号模拟器、电离层和对流层仿真单元，其中：GPS星间链路信号模拟器与评估计算模块相连并传输星间伪距测量值，电离层和对流层仿真单元与监测站模块中的电离层闪烁监测单元相连并传输电离层模型参数或电子浓度、电离层闪烁参数、对流层模型参数。
3.根据权利要求1所述的方法，其特征是，所述的监测站模块包括：分别与软件接收机模块连接的CORS站单元、IGS服务器单元、四模参考站接收机单元和电离层闪烁监测单元，其中：CORS站单元向软件接收机模块输出接收机观测信息，IGS服务器单元向软件接收机模块输出卫星星历和大气环境文件，四模参考站接收机单元向软件接收机模块输出伪距测量值、信号功率、载噪比、导航电文和定位坐标，电离层闪烁监测单元向软件接收机模块输出电离层闪烁监测信息。
4.根据权利要求1所述的方法，其特征是，所述的软件接收机模块包括：分别与评估计算模块相连接的卫星自主完好性监测单元、接收机自主完好性监测单元、定位结果计算单元、电文检测单元和载噪比计算单元、大气延迟计算单元，其中：卫星自主完好性监测单元向评估计算模块输出卫星在空间段的完好性信息，接收机自主完好性监测单元向评估计算模块输出卫星在用户端的完好性信息，定位结果计算单元向评估计算模块输出用户接收机的实际定位信息，电文检测单元向评估计算模块输出电文一致性检测结果，载噪比计算单元向评估计算模块输出实际得到的载噪比数据，大气延迟计算单元向评估计算模块输出大气延迟带来的传播误差。
5.根据权利要求1所述的方法，其特征是，所述的评估计算模块包括：脆弱性综合评估单元以及分别与之相连接的SAIM结果信息单元、RAIM结果信息单元、定位误差评估单元、电文一致评估单元、信号质量评估单元以及大气环境评估单元，其中：脆弱性综合评估单元根据接收到的卫星在空间段的完好性信息、卫星在用户端的完好性信息、用户接收机的实际定位信息、电文一致性检测结果、实际得到的载噪比数据、大气延迟带来的传播误差计算得到系统脆弱性信息。
6.根据权利要求1所述的方法，其特征是，所述的定位误差结果包含：卫星可见性分歧度Δn，导航电文轨道偏差值|Δp|，载噪比偏差值 电离层IONEX模型与Klobuchar模型比值rion，完好性监测结果以及定位误差实测值。
7.根据权利要求1所述的方法，其特征是，所述的步骤四具体是指：监测当前可见星情况，分析与预报星历可见卫星数、卫星号的异同；对于标示为健康的卫星信号，连续不间断地监测接收机的 值，报告 突然降低或突然升高，以及偏离周期性规律值或理论值的情况，采用自适应算法：通过计算一个探测窗范围内的 的平均值，与设定的阈值进行比较，当平均值小于阈值则认为是系统电文异常，否则认为出现产生式欺骗。
8.根据权利要求1所述的方法，其特征是，所述的脆弱性监测系统根据门限进行逐步比对得出当前GPS脆弱与否的结论，并计算出当前脆弱度指数，经各自对应的权重进行加权计算后实现脆弱性因素的确定。
9.根据权利要求7所述的方法，其特征是，所述的与设定的阈值进行比较，通过判决公式实现，具体为：|mean(Z)|＞|mean(Y)|+2·σ(Y)时，认定有干扰发生，其中：mean(Z)表示探测窗范围内的 的平均值，Y表示设定阈值，σ表示阈值波动的标准差。
用于GPS的基于影响因素特征的脆弱性监测方法
技术领域
[0001] 本发明涉及的是一种卫星导航领域的系统及方法，具体是一种用于GPS（Global Positioning System）的基于影响因素特征的脆弱性监测系统及方法。
背景技术
[0002] 全球导航卫星系统脆弱性是指全球导航卫星系统在各类因素影响下服务性能缺失或衰减的程度，是系统的固有属性。其内涵指在各类故障、干扰因素发生时，系统及服务质量受到影响的程度。其外延主要包括：卫星故障失效、星间链路通信干扰、星地链路通信干扰、大气环境影响（电离层/对流层）、多径效应、终端故障等影响因素。随着全球导航卫星系统在世界范围内军事、国民经济、日常生活等几乎所有领域的广泛使用，人类越来越依赖全球卫星导航系统所提供的定位、导航与授时服务。因此对于所有全球卫星导航系统，分析其脆弱性程度，判断其能否提供稳定服务已成为重要研究课题。
[0003] 经对 现有技 术文 献的 检索 发现，Saint Germain en Laye.撰文 GNSS Vulnerability and Mitigation Measures[C]IALA Recommendation R-1292012.“全球卫星导航系统脆弱性和缓解方法”，文中提出了一种全球卫星导航系统的脆弱性主要影响因素及其影响的缓解方法，但该技术针对特定脆弱性因素过程复杂，且没有定量的衡量分类标准，不具有可操作性，因此构建有效的脆弱性监测系统对于全球导航卫星系统脆弱性评估问题具有重要的意义。
发明内容
[0004] 本发明针对现有技术存在的上述不足，提出一种用于GPS的基于影响因素特征的脆弱性监测系统及方法，解决了现有卫星导航系统脆弱性影响因素的监测过程复杂且缺乏定量判断标准的技术问题，在不引入可见的计算误差的前提下，大大提高脆弱性影响因素的监测效率，使得全球导航卫星系统的脆弱性监测与判断有系统、快速、准确的实现途径。
[0005] 本发明是通过以下技术方案实现的：
[0006] 本发明涉及一种用于GPS的基于影响因素特征的脆弱性监测系统，包括：控制计算模块、评估计算模块、空间段脆弱性仿真模块、软件接收机模块和监测站模块，其中：控制计算模块与空间段脆弱性仿真模块相连并传输脆弱性因素生成参数，空间段脆弱性仿真模块分别向监测站模块传输电子浓度、对流层模型参数、电离层闪烁参数和/或电离层模型参数，并向软件接收机模块传输星间链路测距值，评估计算模块分别与空间段脆弱性仿真模块、监测站模块和控制计算模块连接，利用空间段脆弱性仿真模块和监测站模块生成的各种脆弱性因素模拟数据计 算评估结果并将结果反馈给控制计算模块。
[0007] 所述的空间段脆弱性仿真模块包括：GPS星间链路信号模拟器、电离层和对流层仿真单元，其中：GPS星间链路信号模拟器与评估计算模块相连并传输星间伪距测量值，电离层和对流层仿真单元与监测站模块中的电离层闪烁监测单元相连并传输电离层模型参数或电子浓度、电离层闪烁参数、对流层模型参数。
[0008] 所述的监测站模块包括：分别与软件接收机模块连接的CORS(Continuous Operational Reference System连续运 行卫星 定位服 务综合系 统)站单 元、IGS(International GPS)服务器单元、四模参考站接收机单元和电离层闪烁监测单元，其中：CORS站单元向软件接收机模块输出接收机观测信息，IGS服务器单元向软件接收机模块输出卫星星历和大气环境文件，四模参考站接收机单元向软件接收机模块输出伪距测量值、信号功率、载噪比、导航电文和定位坐标，电离层闪烁监测单元向软件接收机模块输出电离层闪烁监测信息。
[0009] 所述的软件接收机模块包括：分别与评估计算模块相连接的卫星自主完好性监测(SAIM，Satellite Autonomous Integrity Monitoring)单元、接收机自主完好性监测(RAIM，Receiver Autonomous Integrity Monitoring)单元、定位结果计算单元、电文检测单元和载噪比计算单元、大气延迟计算单元，其中：卫星自主完好性监测单元向评估计算模块输出卫星在空间段的完好性信息，接收机自主完好性监测单元向评估计算模块输出卫星在用户端的完好性信息，定位结果计算单元向评估计算模块输出用户接收机的实际定位信息，电文检测单元向评估计算模块输出电文一致性检测结果，载噪比计算单元向评估计算模块输出实际得到的载噪比数据，大气延迟计算单元向评估计算模块输出大气延迟带来的传播误差。
[0010] 所述的评估计算模块包括：脆弱性综合评估单元以及分别与之相连接的SAIM结果信息单元、RAIM结果信息单元、定位误差评估单元、电文一致评估单元、信号质量评估单元、大气环境评估单元以及脆弱性综合评估单元，其中：脆弱性综合评估单元根据接收到的卫星在空间段的完好性信息、卫星在用户端的完好性信息、用户接收机的实际定位信息、电文一致性检测结果、实际得到的载噪比数据、大气延迟带来的传播误差计算得到系统脆弱性信息。
[0011] 本发明涉及上述系统的脆弱性监测方法，包括以下步骤：
[0012] 步骤一、脆弱性监测系统启动，在监测设备上设定需要监测的目标系统信号，得到所有卫星信号的集合M={x1(t)，x2(t)，x3(t)，……，xN(t)}，并且得到初始集合的大小N，其中：xi(t)表示在选定频点上不同类型的信号，N表示信号类型总数。
[0013] 步骤二、根据目标频段的系统信息，从脆弱性监测系统处获得GPS的卫星导航信号的信息列表；
[0014] 所述的系统信息是指：当前系统的选定频点，包括GPS的L1、L2、L5频点。
[0015] 所述的信息列表包括：当前可见星信息，当前定位精度信息，当前卫星几何分布信息（GDOP），当前实时导航电文信息，实时导航信号质量（C/No）信息，实时对流层、电离层延迟信息，当前的RAIM/SAIM结果信息，接收机定位误差信息。
[0016] 步骤三、对实时导航电文进行电文一致性评估：监测来自脆弱性监测系统实测的实时导航电文根据其中的轨道信息，进行当前卫星坐标的推算，将得到的卫星坐标推算结果与来自IGS服务器的精确预报星历中的卫星坐标进行对比，得到每颗卫星的当前推算坐标与预测坐标的偏差值 其中：x0、y0、z0分别表示根据实测导航电文推算出的卫星坐标，x、y、z表示从IGS服务器单元精确预报星历中得到的卫星坐标，当偏差值未超过设定阈值则判定导航电文正确并执行步骤四，否则判定导航电文异常并进一步判定为发生了系统电文异常或者产生式欺骗干扰。因为产生式欺骗干扰会先导致接收机信号质量的下降，以至接收机失锁才能完成欺骗，故要结合步骤四的信号质量评估来判断到底是系统电文异常还是产生式欺骗干扰。
[0017] 步骤四、信号质量评估：监测接收到的最小信号功率是否达到或超过ICD指定数值，当小于等于ICD制定数值则认为出现产生式欺骗，否则认为是系统电文异常，具体包括：监测当前可见星情况，分析与预报星历可见卫星数、卫星号的异同；对于标示为健康的卫星信号，连续不间断地监测接收机的 值，报告 突然降低或突然升高，以及偏离周期性规律值或理论值的情况，采用自适应算法：通过计算一个探测窗范围内的 的平均值，与设定的阈值进行比较，当平均值小于阈值则认为是系统电文异常，否则认为出现产生式欺骗。
[0018] 所述的比较通过判决公式实现，具体为：|mean(Z)|>|mean(Y)|+2·σ(Y)时，认定有干扰发生，其中：Z表示探测窗范围内的 的平均值，Y表示设定阈值，σ表示阈值波动的标准差。
[0019] 所述的阈值通过估计的计算 窗两边的 噪声来确定。
[0020] 步骤五、采用IONEX和Klobuchar两个模型结合的方式进行基于电离层异常的大气环境评估，具体为：设定Klobuchar延迟的2倍作为延迟阈值，当IONEX延迟值延迟超出延迟阈值时，即当 时，则判定电离层异常发生，否则认为电离层正常，其中：ΔIONEX表示根据IONEX模型实时计算的电离层延迟，ΔKlobuchar表示根据Klobuchar模型计实时算的电离层延迟。
[0021] 步骤六、完好性监测评估：软件接收机模块通过卫星间观测数据，利用SAIM算法，从系统内判断卫星故障状况，并得到当前异常工作的卫星信息，同时利用RAIM算法判断卫星通道故障情况，并得到当前测量值错误的通道信息。
[0022] 步骤七、接收机定位误差监测：软件接收机模块得到当前定位误差并与设定阈值进行比 较，当当前定位误差未超过设定阈值，则判定定位误差处于可接受范围判断无脆弱性事件，否则判定定位误差过大，判断肯定存在脆弱性影响因素，最后软件接收机模块根据比较结果生成定位误差结果。
[0023] 所述的定位误差结果包含：卫星可见性分歧度Δn，导航电文轨道偏差值|Δp|，载噪比偏差值 电离层IONEX模型与Klobuchar模型比值rion，完好性监测结果以及定位误差实测值。
[0024] 步骤八、将步骤三到步骤七得到的卫星可见性、导航电文、信号质量、大气环境、完好性监测及定位误差结果，输入至脆弱性监测系统，得到系统脆弱性信息。
[0025] 所述的脆弱性监测系统根据门限进行逐步比对得出当前GPS脆弱与否的结论，并计算出当前脆弱度指数，经各自对应的权重进行加权计算后实现脆弱性因素的确定。
[0026] 技术效果
[0027] 与现有技术相比，本发明首先在监测设备上设定需要检测的目标信号，按GPS系统组成的各个分段，对卫星信号进行分类，将待检测卫星信号分别与已分类的已知卫星信号进行脆弱性影响因素的监测，最终得到脆弱性影响因素的判断结果。
[0028] 本发明对卫星导航信号进行脆弱性影响因素的监测，以便快速判定不同卫星导航信号的脆弱性情况。该监测系统是基于卫星导航系统各部分脆弱性影响因素，分别设定脆弱性门限，且实施的每一步骤都有定量的衡量分类标准，在不引入可见的计算误差的前提下，大大提高脆弱性各影响因素的监测效率和系统性，使得全球导航卫星系统的脆弱性监测与判断有快速、准确的实现途径。
附图说明
[0029] 图1为本发明脆弱性监测系统示意图。
[0030] 图2为脆弱性因素影响结果分析图。
[0031] 图3为脆弱性影响因素评估算法流程图。
[0032] 图4为多卫星电离层模型延迟及对比率示意图。
[0033] 图5为载噪比 理论值设定和实际观测值对比示意图。
具体实施方式
[0034] 下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0035] 实施例1
[0036] 如图1所示，本实施例涉及的脆弱性监测系统，包括：控制计算模块、评估计算模块、空间段脆弱性仿真模块、软件接收机模块和监测站模块，其中：控制计算模块与空间段脆弱性仿真模块相连并传输脆弱性因素生成参数，空间段脆弱性仿真模块分别向监测站模块传输电子 浓度、对流层模型参数、电离层闪烁参数和/或电离层模型参数，并向软件接收机模块传输星间链路测距值，评估计算模块分别与空间段脆弱性仿真模块、监测站模块和控制计算模块连接，利用空间段脆弱性仿真模块和监测站模块生成的各种脆弱性因素模拟数据计算评估结果并将结果反馈给控制计算模块。
[0037] 所述的控制计算模块内置综合控制服务器单元。
[0038] 所述的空间段脆弱性仿真模块包括：GPS星间链路信号模拟器、电离层和对流层仿真单元，其中：GPS星间链路信号模拟器与评估计算模块中的卫星自主完好性监测单元相连并传输星间伪距测量值，电离层和对流层仿真单元与监测站模块中的电离层闪烁监测单元相连并传输电离层模型参数或电子浓度、电离层闪烁参数、对流层模型参数。
[0039] 所述的监测站模块包括：分别与软件接收机模块中对应单元连接的CORS站单元、IGS服务器单元、四模参考站接收机单元和电离层闪烁监测单元，其中：CORS站单元与接收机自主完好性监测单元和定位结果计算单元相连并传输接收机观测信息，IGS服务器单元与电文检测单元相连并传输卫星星历和大气环境文件，四模参考站接收机单元与载噪比计算单元相连并传输伪距测量值、信号功率、载噪比、导航电文和定位坐标，电离层闪烁监测单元与大气延迟计算单元相连并传输电离层闪烁监测信息。
[0040] 所述的软件接收机模块包括：分别与评估计算模块中的各个对应单元相连接的卫星自主完好性监测单元、接收机自主完好性监测单元、定位结果计算单元、电文检测单元和载噪比计算单元、大气延迟计算单元，其中：卫星自主完好性监测单元与SAIM结果信息单元相连并传输卫星在空间段的完好性信息，接收机自主完好性监测单元与RAIM结果信息单元相连并传输卫星在用户端的完好性信息，定位结果计算单元与定位误差评估单元相连并传输用户接收机的实际定位信息，电文检测单元与电文一致评估单元相连并传输电文一致性检测结果，载噪比计算单元与信号质量评估单元相连并传输实际得到的载噪比数据，大气延迟计算单元与大气环境评估单元相连并传输大气延迟带来的传播误差。
[0041] 所述的评估计算模块包括：脆弱性综合评估单元以及分别与之相连接的SAIM结果信息单元、RAIM结果信息单元、定位误差评估单元、电文一致评估单元、信号质量评估单元、大气环境评估单元以及脆弱性综合评估单元，其中：脆弱性综合评估单元根据接收到的卫星在空间段的完好性信息、卫星在用户端的完好性信息、用户接收机的实际定位信息、电文一致性检测结果、实际得到的载噪比数据、大气延迟带来的传播误差计算得到GNSS脆弱性信息。
[0042] 本实施例根据上述系统进行脆弱性影响因素的判断，具体包括以下实施条件：设计待检测卫星信号x1(t)工作在L1频段上，x1(t)信号的伪随机码码长Ｌ为1023，信号码片速率fc为1.023MHz、导航数据速率fd为50Hz。可以引起脆弱性的因素有很多，这里不妨假设当前只发 生了电离层异常以及一颗卫星的信号故障。
[0043] 本实施例上述系统，具体通过以下步骤实现脆弱性监测：
[0044] 步骤一：系统启动，在监测设备上设定需要检测的目标频段为L1，得到位于该频段上的所有卫星信号的集合M={x1(t)，GPS CA，GPS L1P(Y)，GPS L1C}，并且得到初始集合的大小N=4；
[0045] 步骤二：根据目标频段的系统信息，获得GPS系统的L1频段的卫星导航信号的信息列表；
[0046] 所述的信息列表包括：当前可见星信息，当前定位精度信息，当前卫星几何分布信息（GDOP），当前实时导航电文信息，实时导航信号质量（C/No）信息，实时对流层、电离层延迟信息，当前的RAIM/SAIM结果信息，接收机定位误差信息。
[0047] 步骤三、对实时导航电文进行电文一致性评估，监测来自脆弱性监测系统实测的实时导航电文，根据其中的轨道信息，进行当前卫星坐标的推算，将得到的卫星坐标推算结果与来自IGS的精确预报星历中的卫星坐标进行对比，得到每颗卫星的当前推算坐标与预测坐标的偏差值 其中：x0、y0、z0表示根据实测导航电文推算出的卫星坐标，x、y、z表示从IGS精确预报星历中得到的卫星坐标，这一偏差值未超过设定阈值则判定导航电文正确。否则判定导航电文异常，需要进一步判定究竟发生了系统电文异常或者产生式欺骗干扰。
[0048] 步骤四、信号质量评估，监测接收到的最小信号功率是否达到或超过ICD上指定的数值；
[0049] 监测当前可见星情况，分析与预报星历可见卫星数、卫星号的异同；对于标示为健康的卫星信号，连续不间断地监测接收机的 值，报告 突然降低或突然升高，以及偏离周期性规律值或理论值的情况，采用自适应算法：通过计算一个探测窗范围内的 的平均值，与设定的阈值进行比较，该阈值通过估计的计算 窗两边的 噪声来确定。判决公式：|mean(Z)|>|mean(Y)|+2·σ(Y)时，认定有干扰发生。监测各频段上信号的码-载波发散（code-carrier divergence）在100-7200秒时间内是否低于ICD规定的数值；监测L1频段上的军用码与其它频段上的民用信号之间的群差分延迟是否低于规定的数值；监测L1和L2频段的差分偏移量（differential bias）的稳定度是否在规定范围内；根据周期性的C/No载噪比评估信号干扰问题，需要历史记录的载噪比数据，而当用户没有记录历史载噪比情况时，可以利用理论C/No载噪比作为阈值来判定实际载噪比异常情况，来实现干扰探测。
[0050] 利用本脆弱性监测系统的观测数据，以2013年1月14日观测的PRN12卫星的载噪比记录为例，进行说明。
[0051] 利用载噪比理论值的计算公式，根据实际参考站接收机观测情况对公式中的参数进行设置——接收机处理损耗：2dB；接收机噪声温度：290K；天线噪声温度：190K；大气损耗参数V：0.125；
[0052] 接收机天线增益按如下公式计算：Crv=-5+15.5*sinEL，其中：，EL为卫星仰角。
[0053] 理论上下限补偿设定：±1dB；已知参考站接收机坐标：[-2840948，4659985，3289597]；计算卫星坐标的快速预报星历文件：igr17231.sp3；载噪比理论的上下限曲线与实际参考站接收机获得的载噪比曲线基本吻合，说明理论载噪比参数设置较符合实际情况。图4中，显示的PRN12卫星载噪比数据时间为GPS时间从32500s至42500s，即当地时间17:01:40至19:48:20。在GPS时间40000s至42500s期间，实际载噪比出现快速下降，略低于设定的理论曲线的下限，即参考站天线对该时刻卫星方向可能出现了遮挡或多径的情况。连续处理10天的参考站GPS L1/L2、BeiDou B1/B2的载噪比数据，载噪比的周期性变化如图所示。
[0054] 对上面记录的载噪比 进行数据分析，可以发现： 的平均差为0.6426dB，标准差为0.7028dB，分析结果表明， 周期性规律非常明显，两天之间的载噪比差异并不大，在设定阈值后，可以用于对信号异常情况的及时发现。此处设定阈值时，需要考虑天气的变化等不确定因素，因此应留一些裕度，2倍的标准差或者平均差作为 变化范围比较合适。
[0055] 步骤五、电离层延迟的评估根据电离层延迟模型和预报星历，结合参考站观测数据来进行监测。采用两个模型结合的方式，其中IONEX延迟值略大于Klobuchar校正延迟，前者比后者获得的误差值比率在1.2。设定Klobuchar延迟的2倍作为阈值，当IONEX延迟值或其他实时模型延迟超出该阈值时，则判定电离层异常发生；
[0056] 步骤六、完好性监测评估：通过卫星间观测数据，利用SAIM算法，从系统内判断卫星故障状况，得到当前哪一颗卫星自身发生了异常；通过接收机观测数据，利用RAIM算法，从接收机角度判断卫星通道故障情况，得到当前哪一个通道的测量值有误。
[0057] 利用SAIM算法，从NASA官网上下载2011年11月22日21:21的RINEX导航文件，文件名为brdc3270.11n，计算星历预测的卫星位置，星历预测的卫星位置与星间双向测距值之和的平均值作差作为真实的轨道残差。轨道残差没有超过相应的最小可检测偏差，说明没有错误；反之，说明有故障。
[0058] 步骤七、接收机定位误差监测，定位误差未超过设定阈值，则判定定位误差处于可接受范围，若超过了设定阈值，则判定定位误差过大，不能够接受。
[0059] 步骤八、将步骤三到步骤七得到的卫星可见性、导航电文、信号质量、大气环境、完好性监测及定位误差结果，输入到脆弱性监测系统，这些结果包含以下参数：导航电文轨道偏差 值|Δp|，载噪比偏差值 电离层IONEX模型与Klobuchar模型比值rion，完好性监测结果以及定位误差实测值，脆弱性监测系统根据由先验信息得到的各个部分的门限进行逐步检验计算得出当前GNSS脆弱与否的结论，并列出具体是什么因素导致当前的脆弱性以及脆弱性的脆弱度指数，并实现影响因素判定。在这一实例中，可以得到导航电文轨道偏差值|Δp|和载噪比偏差值 未超过门限，电离层IONEX模型与Klobuchar模型比值rion超过门限，完好性监测出一个卫星信号的异常以及定位误差超过阈值的结果。因此可以很容易得到当前脆弱性影响因素为电离层异常以及一颗卫星的信号故障，并可以综合电离层IONEX模型与Klobuchar模型比值rion以及定位误差超过阈值的大小来计算当前的脆弱度指数。
[0060] 再进行其它因素的假设判断：假设当前发生了阻塞干扰，则可以发现所有卫星的载噪比偏差值 异常甚至得不到载噪比的实测值，定位误差超过门限值很多，很容易得到发生了阻塞干扰的结论，且脆弱度指数相当高。假设当前是卫星电文发生了异常，那么导航电文轨道偏差值|Δp|会超过门限，可以得到是电文异常的结论。与电文异常的情况比较类似，产生式欺骗干扰也会导致导航电文轨道偏差值|Δp|超过门限，但欺骗干扰同时也会导致短时间或者持续的载噪比异常，这样可以通过噪比偏差值 是否有超过阈值的情况来判断到底是系统电文异常还是产生式欺骗干扰。再假设当前发生的是多径干扰或者转发式欺骗干扰，那么包括载噪比偏差值 电离层IONEX模型与Klobuchar模型比值rion在内的因素都不会超过阈值，完好性监测结果也是正常，只有定位误差会超过门限，此时便可以根据排除法的思想判断出当前的脆弱性影响因素为多径效应或者转发式欺骗干扰。
[0061] 本实施例的脆弱性影响因素监测方法，相比现有分项技术监测方法，是一种系统的脆弱性监测手段，而且得到的脆弱性分离系数误差小很多。
