[发明专利]一种电离层闪烁发生概率季节变化相位角翕动预测方法

权利要求书

1.一种电离层闪烁发生概率季节变化相位角翕动预测方法，其特征在于，包括如下步骤：

步骤1，用傅里叶函数表示电离层闪烁发生概率季节变化相位角φ_doy随年积日的变化，其函数形式为：φ_doy＝a₀+a₁*cos(doy*w)+b₁*sin(doy*w)：

步骤1A，首先计算某一年全部的电离层闪烁发生概率季节变化相位角φ_doy，doy＝1,2,3,......,d_max，，d_max为一年总天数，闰年d_max＝366，平年d_max＝365；

步骤1B，构成两个向量：一是作为自变量的输入向量X_doy＝1,2,3,......,d_max，二是作为因变量的目标输出向量

步骤1C，采用非线性函数拟合工具对X_doy和Y_φ进行傅里叶函数拟合，傅里叶函数形式为：φ_doy＝a₀+a₁*cos(doy*w)+b₁*sin(doy*w)，获得傅里叶函数系数a₀,a₁,b₁和w，其同时满足向量函数关系式：Y_φ＝a₀+a₁*cos(X_doy*w)+b₁*sin(X_doy*w)；

步骤2，利用经验模式预测季节变化相位角φ_doy的傅里叶函数角频率的翕动修正量Δw：

步骤2A，预测翕动修正量Δw的经验模式是年平均太阳活动指数Rz的函数，并且年平均太阳活动指数Rz存在适用范围限制，Rz适用范围检查方法如下所示：

其中，R_in为限定值域后的年平均太阳活动指数；

步骤2B，依据限定值域后的年平均太阳活动指数R_in，利用经验模式预测翕动修正量Δw，其计算方法如下：

Δw＝-(a_Δ0+a_Δ1*cos(R_in*w_Δ)+b_Δ1*sin(R_in*w_Δ))，式中，a_Δ0＝0.008056，a_Δ1＝-0.005449，b_Δ1＝0.005821，w_Δ＝-0.01081；

步骤3，利用修正后的傅里叶函数角频率w+Δw计算修正后的电离层闪烁发生概率季节变化相位角具体计算公式如下：