[发明专利]一种初相稳定性判定方法、存储介质及电子设备

权利要求书

1.一种基于VLBI记录系统的初相稳定性判定方法，其特征在于，包括：

步骤1，生成基准信号S_ref；

步骤2，按预设时间段连续获取VLBI记录系统的记录数据S_n(t)；

步骤3，对所述基准信号S_ref以及所述记录数据S_n(t)进行互相关计算；

步骤4，根据互相关计算结果Scor_n(t)进行残余相位序列计算；

步骤5，根据残余相位序列的计算结果Φ_n(t)计算相位变化率均值φ_ref；

步骤6，根据相位变化率均值φ_ref计算预设时间段端点时刻的差分相位以及基准差分相位Φ_ref；

步骤7，根据所述差分相位以及所述基准差分相位Φ_ref判断VLBI记录系统的初相是否稳定；

其中，步骤1具体为：

通过高稳频率源生成基准信号S_ref，基准信号S_ref的计算公式为：

S_ref＝exp(j·2π·f_ref·t)

其中，f_ref为基准信号频率，t为信号时间，j为虚数单位；

步骤2具体为：

获取VLBI记录系统的记录数据S_n(t)的公式为：

S_n(t)＝exp[j·2π·(f_ref-f_Ln)·t_M]

其中，n表示第n个时间段，f_Ln为n个时间段的本振频率，

t_M为该时间段内的采样点序列，F_n为第n个时间段的采样频率；

步骤3具体为：

对所述基准信号S_ref以及所述记录数据进行互相关计算，得到互相关计算结果Scor_n(t)，公式为：

其中，conj()为复共轭计算，M为进行互相关计算的数据点数，S_n(t_i)为第n个预设时间段的VLBI记录系统的记录数据；

步骤4具体为：

根据互相关计算结果Scor_n(t)进行残余相位序列计算，得到残余相位序列的计算结果Φ_n(t)，计算公式为：

Φ_n(t)＝angle{Scor_n(t)}

其中，angle{}为取幅角计算；

步骤5具体为：

根据残余相位序列的计算结果Φ_n(t)计算相位变化率均值φ_ref，公式为：

其中，N为全部N个时间段，R为每个时间段内互相关计算得到的残余相位数，t_i为第i个时间段的起始时刻，Φ_n(t_i+1)为第n个预设时间内的第i+1个残余相位，Φ_n(t_i)为第n个预设时间内的第i个残余相位；

步骤6具体为：

根据相位变化率均值φ_ref计算预设时间段端点时刻的差分相位以及基准差分相位Φ_ref，基准差分相位Φ_ref计算公式为：

Φ_ref＝φ_ref·(t_n+1,1-t_n,M)

差分相位计算公式为：

其中，t_n+1,1为第n+1个时间段的第一个点的时间，t_n,M为第n个时间段的最后一个点的时间，Φ_n+1(t₁)为第n+1个预设时间内的第一个残余相位，Φ_n(t_M)为第n个预设时间内的第M个残余相位；步骤7具体为：

根据所述差分相位以及所述基准差分相位Φ_ref差值的绝对值与阈值的比较判断VLBI记录系统的初相是否稳定，具体为：

当时，判定VLBI记录系统的初相稳定；

其中，为初相稳定性指标。