1.一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于，具体步骤如下：步骤S1：建立包含有不确定因子的随机系统，不确定因子包括隐变量和加性干扰；

步骤S2：设计关于隐变量和加性干扰的Kalman滤波器；

步骤S3：基于步骤S2中的Kalman滤波器和随机系统设计出求解随机系统状态状态变量估计值的自适应滤波器，并根据自适应滤波器对随机系统进行实时滤波。

2.根据权利要求1所述的一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于：在步骤S1中，随机系统的状态模型和观测模型如下：(1)

   (2)

其中， 为离散时间序列， 为系统状态转移矩阵， 为观测矩阵， 和均为系数矩阵， 为状态向量， 为观测向量， 和为系统中的不确定因子，其中 为不确定因子中的隐变量； 表示n维实数集， 表示m维实数集；

和 为不相关的白噪声序列，满足如下公式：    (3)

(4)

   (5)

其中， 表示白噪声 序列的方差； 表示Kronecker积 ；表示离散时间序列中的时刻点； 表示白噪声 序列的方差。

3.根据权利要求2所述的一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于：在步骤S2中，关于隐变量 的Kalman滤波器设计步骤如下：步骤S21a：建立一步滞后的 动态模型，如下：   (6)；

其中， 和 分别表示隐变量的状态转移矩阵和零均值的高斯白噪声；

步骤S22a：为了对隐变量 进行估计，通过随机系统的状态模型简化随机系统的观测模型，简化过程如下：(7)；

统计特性如下：

  (8)

   (9)

其中， 表示为白噪声 的方差； 和 表示为零均值的不相关的白噪声；

步骤S23a：基于隐变量 的动态模型和观测模型，建立 的Kalman滤波器并求解估计值 和估计误差协方差矩阵 。

4.根据权利要求3所述的一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于：步骤S23a具体为：时间更新：

通过公式(6)，得隐变量 预测估计值和预测估计误差协方差矩阵：   (10)

   (11)

其中 为白噪声 的方差；

测量更新：

根据公式(7), 可得隐变量 的测量预测值及测量预测误差，如下：  (12)

  (13)

其中， 为 时刻到 时刻的 的测量预测值， 为时刻到 时刻的 的测量预测误差， 为已知的 时刻最优估计值，已知的 时刻估计误差值， 为 时刻对 时刻的预测值，为 时刻对 时刻的预测误差值； 为 时刻对 时刻的隐变量 估计值， 为 时刻对 时刻的隐变量 估计误差值；

隐变量 的Kalman滤波器设计如下：(14)

相应的估计误差为：

(15)

结合正交原理得到Kalman滤波器中 增益矩阵如下：  (16)

其中， 为  时刻的预测估计值误差协方差矩阵， 为预测误差协方差矩阵；

隐变量 的估计误差协方差矩阵：(17)。

5.根据权利要求4所述的一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于：在步骤S2中，关于加性干扰的Kalman滤波器设计步骤如下：步骤S21b：建立加性干扰的动态模型，如下：  (18)

其中， 表示加性干扰的状态转移矩阵, 表示零均值的白噪声；

步骤S22b：基于加性干扰的状态模型，得到 的预测值及预测误差协方差矩阵：(19)

  (20)

其中， 为白噪声 的方差；

步骤S23b：根据公式(2)得 的测量预测值和测量预测误差如下：    (21)

(22)

的Kalman滤波器设计如下：(23)

相应的估计误差：

   (24)

利用正交原理再结合估计误差以及观测模型可得 增益矩阵如下：   (25)

根据式(24)得到 如下：  (26)。

6.根据权利要求5所述的一种动态系统加性干扰与状态估计的实时滤波方法，其特征在于：步骤S3具体如下：步骤S31：时间更新，根据公式(1)可得 的状态预测值、状态预测误差及预测估计值误差协方差矩阵如下：   (28)

(29)

(30)；

步骤S32：测量更新，根据状态预测值和观测模型可得测量预测值和测量预测误差如下：  (31)

  (32)；

步骤S33：设计状态变量 的自适应滤波器如下：   (33)

相应的状态估计误差如下：  (34)

其中， 为增益矩阵，利用正交原理求解得增益矩阵 如下：    (35)

利用状态估计误差方差阵计算公式，并根据(34)得到 ，  (36)；

其中， 为单位矩阵。