[发明专利]一种基于时间序列的光纤电流互感器随机噪声实时滤波方法

权利要求书

1.一种基于时间序列的光纤电流互感器随机噪声实时滤波方法，其特征在于，包括如下步骤：

(1)采集光纤电流互感器输出的某相电流数据{x0(t)}；

(2)电流序列平均滤波处理；对有效序列进行平均滤波处理，计算每十个数据的均值作为新的序列{x1(t)}；

(3)建立时间序列模型；

(a)平稳性检验；采用单位根检验法则，不满足平稳性要求则进行差分处理，以获取平稳的电流序列{x_n}；

(b)正态性检验；采用偏峰态检验法，电流序列{x_n}具有以下4个表示其总体概率密度函数的参数如下：

均值：

方差：

标准偏度系数：

标准峰度系数：

当g₁和g₂约等于0时，可以认为电流序列{x_n}为正态时序；

(c)周期性检验；周期性检验用来识别光纤电流互感器输出数据中是否包含有随机量以外的周期性分量，这样在解释数据分析的结果时可以避免出现错误；周期性检验的方法是直接考察从输出数据中得到的概率密度函数或自相关函数或功率谱密度函数的图形；

(d)建立时间序列模型；依次选取p和q值计算各模型的AIC值，然后选择最小的AIC值，确定时间序列模型的阶次，即p和q值，采用最小二乘法拟合出模型参数，写出时间序列模型；

(4)写出与确定的时间序列模型相对应的卡尔曼滤波方程的状态空间模型；

状态方程：X_k＝AX_k+BV_k；

输出方程：Y_k＝CX_k+W_k；

其中，V_k和W_k的统计特性为：

系统的状态方程为过程噪声为V_k＝[r_k,0]^T；

对于AR(p)模型，

V_k＝[r_k,0]^T；

对于ARMA(p,q)模型，

V_k＝[r_k,…,r_k-q]^T；

卡尔曼滤波处理，采用以下卡尔曼滤波器对光纤电流互感器的时间序列输出信号进行滤波处理

K_k＝P_k,k-1C^T[CP_k,k-1C^T+R_k]^-1

P_k,k-1＝AP_k,k-1A^T+BQ_k-1,kB^T

P_k,k＝[I-K_kC]P_k,k-1

式中，为滤波状态的进一步估计，为k时刻滤波器的状态，K_k为k时刻滤波器的增益矩阵，R为系统量测噪声误差，Q为系统过程噪声方差，P_k,k为滤波器误差协方差矩阵，为k时刻滤波器的输出。

2.如权利要求1所述的基于时间序列的光纤电流互感器随机噪声实时滤波方法，其特征在于，步骤(2)中采用平均值滤波的方法对原始数据进行处理。

3.如权利要求1所述的基于时间序列的光纤电流互感器随机噪声实时滤波方法，其特征在于，步骤(3)中采用AIC信息准则来确定时间序列模型的阶次，计算各个模型的AIC值，取AIC值最小的模型为使用模型；采用快速算法RLS进行AR模型的参数估计，采用长自归白噪声估计法建立ARMA模型；采用检验模型残差是否为白噪声来确定模型的可用性。

4.如权利要求1所述的基于时间序列的光纤电流互感器随机噪声实时滤波方法，其特征在于，步骤(4)中采用卡尔曼滤波的方法对建立的时间序列模型进行滤波处理。