本发明涉及一种基于自相关分段迭代法的电离层tec单站预报方法，属于电离层tec预报领域。
背景技术
电离层的参数预报通常分为两种：一种是长期预报，另一种是短期预报。长期预报的预报提前量一般为3个月以上，是指预测将来若干个月的电离层特性的平均情况。该措施是先建立各电离层月中值参量与其相应的表征太阳活动指数的相关关系，然后将预报的太阳活动指数代入相关方程，从而求得电离层参量，实现预报。针对fof2和m(3000)f2，在我国通常用以下两种方法进行预测：1)ccir方法：预测全球上述电离层特征参数的月中值，用itu的ccir340报告建议的方法，以r12为随太阳活动周期变化的参量，并以fof2与r12线性相关为基础。2)亚大方法：亚洲大洋洲部分区域上述特征参数的月中值的预测方法，它采用由重庆站的fof2导出的预报指数ic，在亚大地区，特别是中国区域，fof2与ic的相关性比fof2与r12的相关性更好，在此基础上形成中国参考电离层cri。
电离层短期预报则通常是对电离层参数进行提前1-3天的预报。因短期预报的实用性，且随着研究的深入，短期预报备受关注。在电离层短期预报方面，有人利用自相关分析法首先对中国地区电离层fof2/tec进行了短期预报，还有人提出了适用于全球电离层总电子含量tec的预报方法，构建了中国电离层tec现报模型等。其中，基于自相关法的电离层tec预报方法，结合克里格插值方法，能形成一种适用于中国地区的电离层tec短期预报方法(分为单站预报和区域重构两部分),现已应用于国家气象局空间天气监测预警中心。不过该方案中，基于自相关方法的单站预报方案是制约其预报精度的重要原因之一,提高其预报精度将有利于整体预报精度的提高。
现有的基于自相关分析法的电离层tec短期预报方法中，需假设电离层tec为平稳随机过程。在此基础上获得的结果显示随着预报提前量的增大，预报误差非线性增大。
技术实现要素：
本发明为了解决现有技术中存在的问题，提供一种可以提高电离层tec短期预报精度的基于自相关不等间隔迭代法的电离层tec单站预报方法。
为了达到上述目的，本发明提出的技术方案为：一种基于自相关不等间隔迭代法的电离层tec单站预报方法，其特征在于，包括以下步骤：
步骤一、基于预报时刻前设定时间段的电离层tec数据作为数据集，计算电离层tec相对变化量的自相关系数；
步骤二、将预报提前量t拆分为若干间隔提前量t1～tm，m为间隔提前量个数，t1～tm均为预报量时间间隔的整数倍，且t1～tm的和等于t，至少一个间隔提前量与其他间隔提前量不同；
步骤三、进行提前t1的预报，获取电离层tec的预报值，将提前t1的预报值加入自相关系数计算的数据集，并再开展提前t2的预报；
步骤四、将提前t1和t2的预报值加入自相关系数计算的数据集，然后进行提前t3的预报，以此类推，直至进行提前tm的预报，则提前tm的预报结果即为提前量t的预报值。
对上述技术方案的进一步设计为：电离层tec的相对变化量计算方法为：
其中，tec(t)为预报时刻前设定时间段中t时刻的tec值，
在进行提前预报时，用电离层tec相对变化量作为预报参数。
所述tm至少为预报量时间间隔的2倍。
在进行提前预报时，在设定时间段内选取与当前间隔提前量最近的一定范围，在该范围内挑选前n个最大的自相关系数，并确定n个最大的自相关系数所对应的相对变化量及自相关系数，利用这些相对变化量及自相关系数计算预报值。
将预报值加入自相关系数计算的数据集的规则为，将当前数据集的第一个数据删除，然后将当前预报值补入数据集的末尾。
本发明的有益效果为：
本发明通过电离层tec的相对变化量作为预报参数，采用不等间隔迭代的方法，一方面一定程度上保证了预报参数的“平稳性”要求；另一方面有效改善了现有方法中随着预报提前量增加预报误差非线性增大的问题。
本发明的不等间隔迭代的方法，同时还减弱了等间隔迭代过程中因为迭代次数过多引起的误差累计增大的问题；新方法的预测精度相比于现有方法精度更高。
附图说明
图1为本发明对比例一的曲线图；
图2为本发明对比例二的曲线图；
图3为本发明对比例三的曲线图；
图4为本发明对比例四的曲线图；
图5为本发明对比例五的曲线图；
图6为本发明对比例六的曲线图。
具体实施方式
下面结合附图以及具体实施例对本发明进行详细说明。
实施例
本实施例中预报提前量t＝1h，将t拆分为0.25h、0.25h和0.5h。则本实施例的基于自相关不等间隔迭代法的电离层tec单站预报方法步骤包括：
步骤1、选取15分钟时间间隔、预报时刻前30天的电离层tec观测序列为数据集(记作数据集一)，计算tec的相对变化量
其中，tec(t)为前三十天中t时刻的tec值，
步骤2、根据下式计算电离层z3自相关系数ρ(τ)，下式中的z(t)为z3值
其中：
其中，e、б和ρ分别表示统计参数的数学期望、方差和自相关系数，n为统计参数的个数，s为统计计算相关参数的中间量，p为待统计参数的序号，τ为时间间隔。
步骤3、对步骤2计算获得的电离层tec自相关系数进行函数拟合：
其中a24，a12和τ0为待定系数，通过最小均方拟合的方法来确定。
步骤4、对电离层tec进行提前k＝0.25h预报时，选取最近k开始的4天范围内挑选前48个最大的自相关系数(n＝48)，并确定所对应的z3值及自相关系数带入(7)式，求解系数λj；
其中i＝0,1,2…n-1；μ为拉格朗日乘数因子；ρ是电离层tec的自相关系数，λj是tj时刻的z3的权重系数。
步骤5、将权重系数及相应的z3值带入下式计算z3的预报值，然后基于(1)式换算成电离层tec值，记作tec15。
步骤6、将观测数据集一的第一个数据删除，在最后补入tec15，并将此时的电离层tec数据集记作数据集二。以数据集二为输入，重新按照步骤1-5计算，并将获得的tec值记作tec30；
步骤7、将数据集二的第一个数据删除，在最后一位补入tec30，并将此时的电离层tec序列记作数据集三。以数据集三作为输入，重新按照步骤1-5计算，其中步骤4中的提前量k＝0.5h，计算获得的tec值记作tec1h。该值即为基于电离层tec相对变化量的不等间隔迭代获得电离层tec预报值。
下面根据多组对比例对本发明的预报误差做进一步说明。
对比例一
本对比例以北京站为例，分别将直接预报值和间接预报值和观测值进行比较，从图1中可明显观察到间接预测值与观测值的吻合度更高，误差小。本对比例中直接预报值采用电离层tec作为预报参量，直接进行提前1小时预报，即现有方法；间接预报值采用电离层tec作为预报参量，然后按本发明的不等间隔迭代法先进行两次提前15分钟预报，然后进行提前1小时预报。
对比例二
本对比例以广州站为例，分别将直接预报值和间接预报值和观测值进行比较，从图2中可明显观察到间接预测值与观测值的吻合度更高，误差小。本对比例中直接预报值采用电离层tec作为预报参量，直接进行提前1小时预报，即现有方法；间接预报值采用电离层tec作为预报参量，然后按本发明的不等间隔迭代法先进行两次提前15分钟预报，然后进行提前1小时预报。
结合对比例一和二可说明采用不等间隔迭代的方法，可显著减小预测误差，预测精度高。对比例三
本对比例以北京站为例，分别将直接预报值和间接预报值和观测值进行比较，从图3中可明显观察到间接预测值与观测值的吻合度更高，误差小。本对比例中直接预报值采用电离层tec作为预报参量，直接进行提前1小时预报，即现有方法；间接预报值采用本发明的电离层tec相对变化量作为预报参量，同时按本发明的不等间隔迭代法先进行两次提前15分钟预报，然后进行提前1小时预报。
对比例四
本对比例以广州站为例，分别将直接预报值和间接预报值和观测值进行比较，从图4中可明显观察到间接预测值与观测值的吻合度更高，误差小。本对比例中直接预报值采用电离层tec作为预报参量，直接进行提前1小时预报，即现有方法；间接预报值采用本发明的电离层tec相对变化量作为预报参量，同时按本发明的不等间隔迭代法先进行两次提前15分钟预报，然后进行提前1小时预报。
结合对比例三和四可说明本发明利用电离层tec相对变化量作为预报参量，同时采用不等间隔迭代的方法，得到的预报值精度高。
对比例五
本对比例以广州站为例，将三种不同预报参量下的不同预报方法得到的提前1小时预报误差进行对比，z1的预报参量为电离层tec，z2的预报参量为电离层tec的绝对变化值(tec减去均值)，z3的预报参量是电离层tec的相对变化值；方法1是直接提前1小时预报(现有方法)，方法2是先提前15分钟预报，然后提前45分钟预报，方法3是提前15分钟后再提前15分钟预报，最后提前30分钟预报(不等间隔迭代方法)，方法4为依次重复4次提前15分钟预报(等间隔迭代方法)，结果如图5所示，“z3-方法3”的绝对误差最小。
对比例六
本对比例以北京站为例，将三种不同预报参量下的不同预报方法得到的提前1小时预报误差进行对比，z1的预报参量为电离层tec，z2的预报参量为电离层tec的绝对变化值(tec减去均值)，z3的预报参量是电离层tec的相对变化值；方法1是直接提前1小时预报(现有方法)，方法2是先提前15分钟预报，然后提前45分钟预报，方法3是提前15分钟后再提前15分钟预报，最后提前30分钟预报(不等间隔迭代方法)，方法4为依次重复4次提前15分钟预报(等间隔迭代方法)，结果如图6所示，“z3-方法3”的绝对误差最小。
本发明的技术方案不局限于上述各实施例，凡采用等同替换方式得到的技术方案均落在本发明要求保护的范围内。