[发明专利]一种航天器多遥测参数分析窗口划分方法

权利要求书

1.一种航天器多遥测参数分析窗口划分方法，其特征在于包括下述步骤：

步骤1，对于航天器遥测数据分析中的n个遥测参数p₁、p₂、….、p_n，取得每个遥测参数p_i的采样周期s_i，1≤i≤n，然后针对遥测参数p_i使用时间间隔s_i通过均值或中位数方式进行等时间间隔采样，获得采样后遥测参数序列S_i＝S(p_i)＝[d_1i，d_2i…，d_mi]，m为采样后遥测参数序列长度；

步骤2，对每个S_i进行快速傅里叶变换，得到变化结果F_i＝FFT(S_i)＝[f₀，f₁，…，f_m-1]，将零频直流分量f₀设置为0，然后计算F_i中的最大值max(F_i)，得到遥测参数p_i的近似周期t_i＝(m*s_i)/max(F_i)；循环计算得到所有参数的周期集合[t₁，t₂，…，t_n]；

步骤3，统计[t₁，t₂，…，t_n]中相同长度周期出现的次数，找到出现次数最多的周期t_max，将[t₁，t₂，…，t_n]中所有长度等于t_max的参数列为备选参数集[p₁，p₂，…，p_v]，v≤n；

步骤4，对[p₁，p₂，…，p_v]中的每一个参数p_j，利用时间周期t_max对采样序列S_j进行预划分，得到子序列集合[sub₁，sub₂，…，sub_c]；对子序列集合中每个序列sub_k，1≤k≤c，取得前两个采样周期的采样点对(s_d1，s_d2)_k；利用采样序列S_j的最大值和最小值进行等间隔分箱处理，得到采样点对(s_d1，s_d2)_k对应的分箱值对(e₁，e₂)_k；所有子序列集合计算得到分箱值对集合[(e₁，e₂)₁，(e₁，e₂)₂，…，(e₁，e₂)_c]；

步骤5，在[(e₁，e₂)₁，(e₁，e₂)₂，…，(e₁，e₂)_c]中得到出现次数最多的分箱值对(e₁，e₂)_q，1≤q≤c，将其作为遥测参数p_j的窗口划分特征；利用该窗口划分特征对采样序列S_j重新进行划分，得到新的划分子序列集合[sub_n1，sub_n2，…，sub_nc]；依此在新的子序列集合中计算后一项与前一项的卷积，即corr_ii＝corr(sub_ii*sub_ii-1)，其中2≤ii≤nc，得到[corr₁，corr₂，…，corr_nc-1]，计算该集合的标准差std_j，该值为遥测参数p_j重新划分窗口后得到的标准差，所有备选参数都计算划分窗口后的标准差，得到[std₁，std₂，…，std_v]；

步骤6，在备选参数标准差集[std₁，std₂，…，std_v]中求得具有最小标准差值的参数p_sel作为基准参数，将基准参数p_sel的划分特征(e₁，e₂)作为基准划分特征；

步骤7，设置最小保护时间间隔t_safe，确保每个划分窗口不小于t_safe；使用最小保护时间间隔t_safe和基准划分特征(e₁，e₂)对基准参数p_sel的采样序列S_sel＝[d_1sel，d_2sel…，d_msel]进行时间窗口的划分，将划分出的每个窗口使用窗口开始时间和窗口结束时间进行表示，得到基准划分窗口集为[(ws₁,we₁)，(ws₂,we₂)，…，(ws_h,we_h)]；

步骤8，依次选取待分析的n个遥测参数中的p_i，从1至h顺序从基准划分窗口集[(ws₁,we₁)，(ws₂,we₂)，…，(ws_h,we_h)]中取出划分时间窗口对(ws_jj,we_jj)，1≤jj≤h，将遥测参数p_i中开始时间≥ws_jj，结束时间≤we_jj的遥测值序列作为该窗口对应的划分结果。