本发明公开了一种多层次深耦合的导航信号异常高精度监测方法，属于信号质量监测评估技术领域。其通过对零中频信号进行多个层次的指标处理，可对信号接收层和信息接收层的明显异常进行快速告警；此外，根据多个层次的监测结果获取时域、频域、调制域和信息域等的先验信息，依据先验信息实现基准信号的构建，并对基准信号进行多个层次的指标处理，并与实测信号进行特征级深耦合比对，最终利用计算得到关联符合度完成异常的高精度检测。本发明可深入信号内部进行精细化分析，解决了传统监测手段无法兼顾高精度和实时性的技术难题，保证对各类信号异常的精准识别与快速告警。
1.一种多层次深耦合的导航信号异常高精度监测方法，其特征在于，包括以下步骤：(1)利用高增益天线对导航信号进行采集接收，并将信号载波和多普勒剥离，得到零中频信号；具体方式为：利用导航接收机接收导航卫星的全部信号分量，并利用跟踪环路中的多普勒信息和先验的中频载波频率信息完成零中频信号S0(n)的获取：S0(n)＝S(n)×[cos(2πfcn/Fs)+jsin(2πfcn/Fs)]×[cos(2πfdn/Fs)+jsin(2πfdn/Fs)]，其中，Fs为采样率，fc为中频频率，S(n)为原始数据，j为虚数单位，n为时间变量；
(2)从信号原始层、信号接收层和信息接收层对零中频信号进行指标计算处理，并通过信号接收层和信息接收层实时发现的异常进行快速告警；具体方式为：(201)将步骤(1)得到的零中频信号进行信号原始层处理：
利用傅立叶变换的功率谱分析，得到零中频信号的频域观测信息，并根据理想信号的功率谱包络进行计算传输函数H(n)：其中，DFT为离散傅里叶变换，Sideal(n)为根据信号接口控制文件格式生成的离线信号；
对零中频信号的时域波形进行码片整周期脉冲码片剥离，得到时域观测指标以及噪声信号波形信息noise(n)：其中， 为Sideal(n)的共轭，rank为计算信号向量的长度的运算；
将零中频信号的实部和虚部进行二维空间的分布绘制，得到调制域复用映射关系Mult[a1(n),...,ak(n)]，其中a1(n),...,ak(n)代表复用信号中的各个信号分量；
(202)将步骤(1)得到的零中频信号进行信号接收层处理：
将零中频信号与本地伪码进行相关处理，得到相关域观测的指标信息与确定的本地伪码ck(n)，将相关后的相关曲线与码跟踪环路的延迟相位进行测距计算，得到测量域的观测信息以及信号与本地时的时延差τk，当信号与伪码相关后的相关曲线无明显相关峰值时，对信号异常进行实时告警；
(203)将步骤(1)得到的零中频信号进行信息接收层处理：
判决零中频信号中包含的导航电文信息dk(n)，当电文信息存在与标准格式的明显错误时，对信号异常进行实时告警；
(3)根据先验的时域、频域、调制域以及信息域信号特征构建理论信号基准模型，得到基准信号，并从信号原始层、信号接收层和信息接收层对基准信号进行指标计算处理；具体方式为：(301)根据步骤(2)得到的信号原始层、信号接收层和信息接收层的指标构建理论信号基准模型，得到基准信号J(n)：J(n)＝IDFT[H(n)]*{Mult[c1(n‑τ1)d1(n‑τ1),...,ck(n‑τk)dk(n‑τk)]+noise(n)}；
式中，IDFT为傅里叶反变换，c1c2...ck为各个分量的本地伪码，d1d2...dk为各个分量的电文，τ1τ2...τk为各个分量观测信息以及信号与本地时的时延差，Mult为步骤(201)得到调制域复用映射关系；
(302)将基准信号J(n)进行信号原始层处理，得到功率谱谱线、功率谱偏差、时域重构波形、调制映射星座图；此时基准信号的传输函数H(n)恒为1，噪声信号波形信息noise(n)恒为0，将基准信号J(n)的实部和虚部进行二维空间的分布绘制，得到调制域复用映射关系Multideal[a1(n),...,ak(n)]，其中a1(n),...,ak(n)代表复用信号中的各个信号分量；
将基准信号J(n)进行信号接收层处理，得到各分量信号与本地伪码的相关峰曲线、相关峰峰值、相关损失、信号载噪比；
将基准信号与本地伪码进行相关处理，得到相关域观测的指标信息与确定的本地伪码ck(n)，将相关后的相关曲线与码跟踪环路的延迟相位进行测距计算，得到测量域的观测信息与信号与本地时的时延差τk；
将基准信号J(n)进行信息接收层处理，得到电文信息原始位、各电文帧的信息位数值，判决得到基准信号中包含的导航电文信息dk(n)；
(4)将所述基准信号和所述零中频信号的各层级指标进行特征级深耦合，计算各层级指标间的关联符合度；步骤(4)中关联符合度的计算方式为：其中，x为每个信号分量信号原始层、信号接收层和信息接收层的观测指标，μ为基准信号各个维域的观测指标，n代表不同信号分量，k和p代表不同维度的指标数量，i和m代表不同维度的指标排列序号；
(5)将各层级指标间的关联符合度与异常检测判决矩阵进行匹配运算，完成对信号异常的检测；具体方式为：(501)对关联符合度矩阵L和异常检测判决矩阵M进行匹配运算，得到异常检测判别矩阵(502)对异常检测判别矩阵进行门限判决，完成对信号异常的检测。
一种多层次深耦合的导航信号异常高精度监测方法
技术领域
[0001] 本发明涉及GNSS导航卫星信号质量监测评估领域，特别是指一种多层次深耦合的导航信号异常高精度监测方法。
背景技术
[0002] 解决空间信号质量精确监测评估问题，是全球卫星导航系统安全运行的重要环节。因此，空间信号质量监测评估是全球卫星导航系统开放服务和安全运行的前提。
[0003] 传统的信号质量监测评估手段是利用标准仪器和导航接收机完成仪器测试指标和接收观测量的符合性判断，无法兼容对信号异常的高精度监测和快速告警，而且对于深层次的信号异常，无法完成精准发现，极大影响系统服务能力的快速恢复和保持，因此，迫切需要发明新方法，攻克新技术，提供新手段，解决空间信号质量精确评定和快速故障定位的难题，从而保障导航系统的可靠服务。
发明内容
[0004] 本发明的目的在于提供一种多层次深耦合的导航信号异常高精度监测方法，可实现对多个层次指标的综合分析评估，挖掘深层次信号异常，提升异常的检测精度，解决传统手段无法兼顾监测实时性和测量精度的技术难题。
[0005] 为了实现上述目的，本发明采用的技术方案为：
[0006] 一种多层次深耦合的导航信号异常高精度监测方法，包括以下步骤：
[0007] (1)利用高增益天线对导航信号进行采集接收，并将信号载波和多普勒剥离，得到零中频信号；
[0008] (2)从信号原始层、信号接收层和信息接收层对零中频信号进行指标计算处理，并通过信号接收层和信息接收层实时发现的异常进行快速告警；
[0009] (3)根据先验的时域、频域、调制域以及信息域信号特征构建理论信号基准模型，得到基准信号，并从信号原始层、信号接收层和信息接收层对基准信号进行指标计算处理；
[0010] (4)将所述基准信号和所述零中频信号的各层级指标进行特征级深耦合，计算各层级指标间的关联符合度；
[0011] (5)将各层级指标间的关联符合度与异常检测判决模版矩阵进行匹配运算，完成对信号异常的检测。
[0012] 进一步的，步骤(1)的具体方式为：
[0013] 利用导航接收机接收导航卫星的全部信号分量，并利用跟踪环路中的多普勒信息和先验的中频载波频率信息完成零中频信号S0(n)的获取：
[0014] S0(n)＝S(n)×[cos(2πfcn/Fs)+jsin(2πfcn/Fs)]×[cos(2πfdn/Fs)+jsin(2πfdn/Fs)]，
[0015] 其中，Fs为采样率，fc为中频频率，S(n)为原始数据，j为虚数单位，n为时间变量。
[0016] 进一步的，步骤(2)的具体方式为：
[0017] (201)将步骤(1)得到的零中频信号进行信号原始层处理：
[0018] 利用傅立叶变换的功率谱分析，得到零中频信号的频域观测信息，并根据理想信号的功率谱包络进行计算传输函数H(n)：
[0019]
[0020] 其中，DFT为离散傅里叶变换，Sideal(n)为根据信号接口控制文件格式生成的离线信号；
[0021] 对零中频信号的时域波形进行码片整周期脉冲码片剥离，得到时域观测指标以及噪声信号波形信息noise(n)：
[0022]
[0023] 其中， 为Sideal(n)的共轭，rank为计算信号向量的长度的运算；
[0024] 将零中频信号的实部和虚部进行二维空间的分布绘制，得到调制域复用映射关系Mult[a1(n),...,ak(n)]，其中a1(n),...,ak(n)代表复用信号中的各个信号分量；
[0025] (202)将步骤(1)得到的零中频信号进行信号接收层处理：
[0026] 将零中频信号与本地伪码进行相关处理，得到相关域观测的指标信息与确定的本地伪码ck(n)，将相关后的相关曲线与码跟踪环路的延迟相位进行测距计算，得到测量域的观测信息以及信号与本地时的时延差τk，当信号与伪码相关后的相关曲线无明显相关峰值时，对信号异常进行实时告警；
[0027] (203)将步骤(1)得到的零中频信号进行信息接收层处理：
[0028] 判决零中频信号中包含的导航电文信息dk(n)，当电文信息存在与标准格式的明显错误时，对信号异常进行实时告警。
[0029] 进一步的，步骤(3)的具体方式为：
[0030] (301)根据步骤(2)得到的信号原始层、信号接收层和信息接收层的指标构建理论信号基准模型，得到基准信号J(n)：
[0031] J(n)＝IDFT[H(n)]*{Mult[c1(n‑τ1)d1(n‑τ1),...,ck(n‑τk)dk(n‑τk)]+noise(n)}；
[0032] 式中，IDFT为傅里叶反变换，c1c2...ck为各个分量的本地伪码，d1d2...dk为各个分量的电文，τ1τ2...τk为各个分量观测信息以及信号与本地时的时延差，Mult为步骤(201)得到调制域复用映射关系；
[0033] (302)将基准信号J(n)进行信号原始层处理，得到功率谱谱线、功率谱偏差、时域重构波形、调制映射星座图；此时基准信号的传输函数H(n)恒为1，噪声信号波形信息noise(n)恒为0，将基准信号J(n)的实部和虚部进行二维空间的分布绘制，得到调制域复用映射关系Multideal[a1(n),...,ak(n)]，其中a1(n),...,ak(n)代表复用信号中的各个信号分量；
[0034] 将基准信号J(n)进行信号接收层处理，得到各分量信号与本地伪码的相关峰曲线、相关峰峰值、相关损失、信号载噪比；
[0035] 将基准信号与本地伪码进行相关处理，得到相关域观测的指标信息与确定的本地伪码ck(n)，将相关后的相关曲线与码跟踪环路的延迟相位进行测距计算，得到测量域的观测信息与信号与本地时的时延差τk；
[0036] 将基准信号J(n)进行信息接收层处理，得到电文信息原始位、各电文帧的信息位数值，判决得到基准信号中包含的导航电文信息dk(n)。
[0037] 进一步的，步骤(4)中关联符合度的计算方式为：
[0038]
[0039] 其中，x为每个信号分量信号原始层、信号接收层和信息接收层的观测指标，μ为基准信号各个维域的观测指标，n代表不同信号分量，k和p代表不同维度的指标数量，i和m代表不同维度的指标排列序号。
[0040] 进一步的，步骤(5)的具体方式为：
[0041] (501)对关联符合度矩阵L和异常检测判决矩阵M进行匹配运算，得到异常检测判别矩阵
[0042] (502)对异常检测判别矩阵进行门限判决，完成对信号异常的检测。
[0043] 本发明与现有技术相比的有益效果为：
[0044] 1、本发明可兼顾实时和高精度异常监测。
[0045] 2、本发明通过与理想基准信号的多层次监测结果的深耦合，深入信号内部进行分析，可实现对深层次异常的精准识别。
附图说明
[0046] 图1是本发明实施例中多层次深耦合导航信号异常高精度监测方法的流程图。
具体实施方式
[0047] 以下结合附图和具体实施方式对本发明做进一步的说明。
[0048] 如图1所示，一种多层次深耦合的导航信号异常高精度监测方法，该方法首先利用高增益天线完成导航信号接收，对多个层次指标的综合分析评估，其中信号接收层和信息接收层完成对异常的快速告警，信号原始层充分计算信号先验信息，支持理论信号构建，最终利用理论基准信号和实测信号的各层级分析结果进行深耦合匹配比对，完成对信号异常的检测。
[0049] 该方法包括以下步骤：
[0050] (1)利用高增益天线完成导航信号的采集接收，并完成信号载波和多普勒剥离，得到零中频信号；
[0051] 为了保证导航信号的异常监测具备很高的精度，使用高增益定向天线完成对卫星的实时温度跟踪，并利用导航信号的接收设备完成该导航卫星包含的全部信号分量的接收，并利用PLL跟踪环路中的多普勒信息和先验已知的中频载波频率信息完成零中频信号S0(n)的获取，获取方法为：
[0052] S0(n)＝S(n)×[cos(2πfcn/Fs)+jsin(2πfcn/Fs)]×[cos(2πfdn/Fs)+jsin(2πfdn/Fs)]
[0053] 其中Fs为采样率，fc为中频频率，S(n)为原始数据。
[0054] (2)从信号原始层、信号接收层和信息接收层对零中频信号进行指标计算处理，并通过信号接收层和信息接收层实时发现的异常进行快速告警；
[0055] 将上一步得到的零中频信号分别进行多层次处理，得到各个层次的监测指标和信号先验信息，多层次处理采用并行的方式进行，进一步提升异常响应速度。
[0056] 对信号进行信号原始层处理，主要对信号的频域频谱、时域波形以及调制复用方式进行监测处理，包括以下三部分：
[0057] a)利用傅立叶变换的功率谱分析，得到频域观测信息，并根据理想信号的功率谱包络进行传输函数H(n)的计算，计算方法为
[0058] b)对时域波形进行码片整周期脉冲码片剥离，得到时域观测指标以及噪声信号波形信息noise(n)，计算方法为
[0059] c)将零中频信号的实部和虚部进行二维空间的分布绘制，得到调制域复用映射关系Mult[a1(n),...ak(n)]，其中ak(n)代表复用信号中的各个信号分量；
[0060] 对信号进行信号接收层处理，主要是从相关域和测量域的角度对信号进行接收处理，包括与本地伪码相关处理；得到相关域观测的指标信息与确定的本地伪码ck(n)，将相关后的相关曲线与码跟踪环路的延迟相位进行测距计算，得到测量域的观测信息与信号与本地时的时延差τk，当信号与伪码相关后的相关曲线无明显相关峰值时，则完成对信号异常的实时告警；
[0061] 对信号进行信息接收层处理，主要是判决得到零中频信号中包含的导航电文信息dk(n)，当电文信息存在与标准格式的明显错误时，完成对信号异常的实时告警；
[0062] (3)根据先验的时域、频域、调制域以及信息域信号特征构建理论信号基准模型，并完成同样的信号原始层、信号接收层和信息接收层指标计算处理；
[0063] 根据上一步监测得到的各类多维域先验信息，包括传递函数，噪声波形、复用方式、信息电文等，实现对信号基准模型构建，构建信号基准J(n)为：
[0064] J(n)＝IDF[H(n)]*{Mult[c1(n‑τ1)d1(n‑τ1),...,ck(n‑τk)dk(n‑τk)]+noise(n)}；
[0065] 将构建的基准信号J(n)进行与实测信号进行同样的信号原始层、信号接收层和信息接收层指标计算处理，得到各层级各维域的指标监测结果；
[0066] (4)将理论信号和实际信号的各层级指标进行特征级深耦合，计算关联符合度；
[0067] 对各维度监测指标进行特征级深耦合，耦合的思路是对理论信号指标与实测信号指标进行范数计算，获取得到二者之间在矢量空间的差异程度，最终将各维度的指标差异程度组合为矩阵，即为关联符合度矩阵，关联符合度计算方法为：
[0068]
[0069] 其中，x为每个信号分量各个维域的观测指标，μ为基准信号各个维域的观测指标，n代表不同信号分量，k和p代表不同维度的指标数量，i和m代表不同维度的指标排列序号；
[0070] (5)将各层级指标间的关联符合度与异常检测判决模版矩阵进行匹配运算，最终完成对信号异常的检测。
[0071] 在完成关联符合度矩阵计算后，需要对各维域指标的关联符合程度进行判断，判断方法为对将关联符合度矩阵L异常检测判决矩阵M进行匹配运算，得到异常检测判别矩阵最终将判别矩阵与各维域关联的异常门限值进行综合比对判决，最终完成对信号异常的检测。
[0072] 本发明首先剥离载波的零中频信号，然后对零中频信号进行多个层次(信号原始层、信号接收层和信息接收层)的指标处理，可对信号接收层和信息接收层的明显异常进行快速告警。此外，根据多个层次的监测结果获取时域、频域、调制域和信息域等的先验信息，依据先验信息构建基准信号，对基准信号进行与零中频信号相同的多个层次的指标处理，并将处理结果与实测信号进行特征级深耦合比对，最终利用计算得到关联符合度完成异常的高精度检测。本发明可深入信号内部进行精细化分析，解决了传统监测手段无法兼顾高精度和实时性的技术难题，保证对各类信号异常的精准识别与快速告警。
