技术摘要：
本发明公开了一种基于SVM多分类算法的太阳射电暴强度判定方法，所述SVM多分类算法的基本原理为针对一个多分类问题，将其中每两个类型组合成一个二分类器，对于每个二分类器，找出一个最优超平面，将样本分为两类，最终通过投票统计的方法得到分类结果。在此过程中，首  全部
背景技术：
随着GNSS等卫星技术在现代社会越来越广泛应用，太阳射电暴强度对GNSS信号的 影响已成为不可忽视的重要部分。由于射电爆发现象通常发生得十分突然，并且常和太阳 活动区的质子爆发、X射线爆发或耀斑爆发等现象同时发生，其辐射强度大、变化剧烈，当抵 达地球时，会造成一系列的地球物理效应，比如磁暴、激光、通讯干扰等。在以往多起太阳射 电爆发期间，太阳射电暴可能导致卫星载噪比下降，定位误差增大，几何精度因子增加，卫 星失锁、可见星数大幅下降等情况，如此可见，太阳射电暴产生的噪声信号是导航信号的影 响因素之一。 由于不同的太阳射电暴强度，对GNSS信号的干扰是不同的，为了使导航系统在太 阳射电暴发生时能够正常工作，太阳射电暴的分类工作有利于对于太阳射电的研究，对维 持卫星导航系统的正常运行具有重大意义。同时，太阳射电暴强度也影响着天气，太阳射电 暴的分类工作有助于对灾害性天气做出预警，具有很好的实用价值。 传统分类方法一般需要具有专业知识的天文学家对太阳射电暴强度进行手工分 类，这种方法费时费力，并且效率和准确率不高。
技术实现要素：
为了解决上述存在问题。本发明提供一种基于SVM多分类算法的太阳射电暴强度 判定方法，太阳射电暴发生时，对GNSS信号产生的影响进行综合考查，结合机器学习中的 SVM算法，判定太阳射电暴的强度。此方法对太阳射电暴强度分类具有可行性，判定结果能 同时给出多个台站太阳射电暴强度类型。对比于传统方法，结合多种因素，节省时间人力成 本，识别准确度和效率相对提高。为达此目的： 本发明提供一种基于SVM多分类算法的太阳射电暴强度判定方法，具体步骤如下: (1)进行数据预处理，计算观测地不同时刻的特征向量，特征向量包括观测地载噪 比下降值、观测地总定位误差、几何精度GDOP因子和卫星失锁数目； (2)根据射电流量的变化，将太阳射电暴划分为太阳射电暴未发生、强度较弱、强 度较强、强度很强4种类型，并且将4种类型的特征向量对应的标签分别赋值为1、2、3和4； (3)将其中每两个类型进行组合构造一个二分类器，得到6组训练集，对于每组训 练集构造一个未知的非线性SVM二分类模型，并对训练样本进行交叉验证选取最优参数，得 出6组训练好的二分类模型； (4)将待识别的卫星导航信号提取出来的特征向量输入步骤(3)中的分类模型中， 每一个模型输出一个结果，如判定结果为p，则第p类计一票，最终预测样本的判定结果为得 票数最多的一类。 5 CN 111553393 A 说　明　书 2/9 页 作为本发明进一步细化，步骤(1)具体包括： (1.1)输入观测地每颗卫星的载噪比、观测地的定位误差、几何精度GDOP因子和卫 星失锁数目； (1.2)对观测地载噪比进行数据预处理： 取当前每颗卫星载噪比下降值的平均值作为该台站载噪比下降值，记为x1，单位： dBHz； snri＝前8个时刻所测卫星载噪比平均值-该时刻载噪比值 式中，snri为每颗卫星载噪比下降值，N为卫星个数； (1.3)观测地定位误差由东北天坐标系下的X、Y、Z三个方向的定位误差组成，三个 方向取标准差作为台站的定位误差，记为x2，单位：m； 式中，ri(i＝1,2,3)分别表示该时刻观测地X、Y、Z方向上的定位误差； (1.4)GDOP因子记为x3，卫星失锁数目记为x4，单位：个。 作为本发明进一步细化，步骤(2)具体包括： 将所述观测地载噪比下降值x1、观测地总定位误差x2、GDOP因子x3、卫星失锁数目 x4作为特征向量x，样本特征向量即为： 按照下表将太阳射电暴强度进行分类； 表1太阳射电暴强度种类划分 式中，i表示第i个样本，x(i)表示第i个样本的特征向量，y(i)表示第i个样本的标 签，(x(i) ,y(i))即为样本点。 作为本发明进一步细化，步骤(3)具体包括： (3.1)6组训练集分别为标签为1和2、1和3、1和4、2和3、2和4对应的训练样本； (3 .2)假设对标签为k和l对应的训练样本构造的二分类器(k，l＝1、2、3、4且k＜ l)，所述非线性SVM二分类模型具体为： f(x)＝wTΦ(x) b 其中Φ(x)为高维线性映射，将4维特征向量x映射到M维的映射函数，M＞＞4， 6 CN 111553393 A 说　明　书 3/9 页 两者为待求参数； (3.3)解决如下优化问题： 式中，ξi为每个样本i的的松弛变量，ξi越大，代表样本点离群越远，C为超参数，C越 大，表示对离群点带来的目标函数损失越重视，y(i)为样本i的标签，x(i)表示样本i的特征向 量； (3.4)求解拉格朗日乘子： 式中，αi，y(i)，i＝1,2,3...n分别是拉格朗日乘子和数据分类标签，n为样本的个 数； 进一步，得到 其中x(s)为拉格朗日 乘子αi≠0项所对应的样本的特征向量，即支持向量，y(s)为对应标签； (3.5)将计算得到的w0和b0代入非线性SVM分类器模型表达式中得： 式中， 为核函数，取高斯核函数为核函数 σ为核参数； ( 3 .6 ) 将 标 签 为 k 和 l 对 应的 训练 样 本的 样 本 点 输 入 M A T L A B 中 的 ClassificationLearner，选择交叉验证并设定折数k，即随机取矩阵Z的n个样本中的 k＝2,3,4......作为训练样本进行机器学习，剩下的 个样本对学习过的模型 进行测试，选择高斯非线性SVM分类器模型，输入参数C、 得到6组训练集，分别对每组 训练集进行机器学习； (3.7)预设足够大的次数，更改参数C、 的值，返回重复步骤(3.6)，记录平均准 7 CN 111553393 A 说　明　书 4/9 页 确率，直至达到预设次数； (3.8)比较所有参数对应的平均准确率，找出最大准确率对应的超参数C和核参数 σ作为最佳参数，该参数设置下训练所得模型作为最优分类模型。 作为本发明进一步细化，步骤(4)具体包括： (4 .1)从待判定的观测地卫星的载噪比、观测地三个方位的定位误差、几何精度 GDOP因子、卫星失锁数目中提取出特征向量，表示为X＝(x(1) ,x(2) ,...,x(N))，N为待检测样 本总数； (4.2)将X输入得到的6个分类模型中，在每个模型中，对于每个特征向量x(t)，会 输出一个判定结果p(p＝1,2,3,4)，对应不同的太阳射电暴强度类型标签，待6个SVM模型均 分类结束后，标签出现次数最多对应的类型即为该地区该时刻的太阳射电暴强度类型。 本发明与现有技术相比，其显著优点是： 本发明提出了一种基于SVM多分类算法的太阳射电暴强度判定方法。所述方法首 先提取观测地卫星载噪比、三个方向定位误差、几何精度(GDOP)因子、卫星失锁数目进行数 据预处理，得到能有效反应太阳射电暴强度的特征向量。基于SVM本身是二元分类器，对不 同太阳射电暴强度类型采用一对一的方法，将四种强度类型分别标记为1、2、3和4，与特征 向量组成样本。将标签为1、2、3和4的样本两两组合，再对每个组合进行SVM二分类模型训 练，得到相应的最优分类模型，总共6组。对于待判定的观测地数据，在提取其特征向量后输 入到每组最优分类模型中进行分类，每组模型输出一个特征向量所对应的标签，采取投票 形式，标签出现次数最多的类型即为所对应的强度类型。传统太阳射电暴需要具有专业知 识的天文学家进行手工分类，这种方法费时费力。本方法成本低，准确度相对增强，而且能 同时处理大批量数据，提高了检测效率，对维持卫星导航系统的正常运行具有重大意义。 附图说明 图1是本发明的一个实施例的流程示意图。