一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法技术方案

本发明专利技术涉及一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，与现有技术相比解决了针对于雷达跟踪未知目标无法保证其滤波精度的缺陷。本发明专利技术包括以下步骤：雷达目标跟踪系统的设定；建立虚拟噪声模型；估值计算；自校正阿尔法贝塔滤波器的滤波。本发明专利技术利用拟合虚拟噪声的方法建立一个虚拟噪声模型，用虚拟噪声代替扰动和系统本身的动态噪声，将参数带扰动的系统转化为带稳定的未知参数和已知噪声的滤波问题，生成自校正α‑β滤波器，利用自校正α‑β滤波器对未知目标进行有效滤波，保证了滤波精度。

A perturbation target filtering method for radar target tracking system based on Alphabeta filter

The present invention relates to a perturbation target filtering method for radar target tracking system based on Alphabeta filter. Compared with the existing technology, it solves the defect that the filtering accuracy of radar target tracking system can not be guaranteed. The invention comprises the following steps: setting up a radar target tracking system; establishing a virtual noise model; estimating calculation; and filtering of a self-tuning Alphabeta filter. The invention establishes a virtual noise model by fitting virtual noise, replaces the disturbance and the dynamic noise of the system itself with virtual noise, transforms the system with disturbance parameters into a filtering problem with stable unknown parameters and known noise, generates a self-tuning alpha_beta filter, and effectively filters the unknown target with self-tuning alpha_beta filter to ensure the filtering precision. Degree.

【技术实现步骤摘要】
一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法
本专利技术涉及雷达目标跟踪
，具体来说是一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法。
技术介绍
传统的自校正α-β滤波器可以有效解决雷达目标跟踪系统的稳定未知参数的滤波问题，而面对带有扰动的未知目标，则很难进行有效滤波。现有技术中，有学者提出用保性能滤波器进行滤波，这种滤波器的特点是用扰动的上界代替实际参数。但是这种方法虽然可以保证滤波的稳定性，滤波精度却较差。因此，如何提出一种针对未知目标能够有效滤波的方法已经成为急需解决的技术问题。
技术实现思路
本专利技术的目的是为了解决现有技术中针对于雷达跟踪未知目标无法保证其滤波精度的缺陷，提供一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法来解决上述问题。为了实现上述目的，本专利技术的技术方案如下：一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，包括以下步骤：雷达目标跟踪系统的设定，对雷达目标跟踪系统进行预设定；建立虚拟噪声模型，拟合虚拟噪声，用虚拟噪声代替扰动和雷达目标跟踪系统的动态噪声；估值计算，将虚拟噪声的方差通过建立函数方程的方法求解出对应估算；自校正阿尔法贝塔滤波器的滤波，将虚拟噪声的方差替换自阿尔法贝塔滤波器动态噪声的方差，生成自校正阿尔法贝塔滤波器，通过自校正阿尔法贝塔滤波器进行雷达目标跟踪系统扰动目标的滤波处理。所述雷达目标跟踪系统的设定包括以下步骤：设定雷达目标跟踪系统表达式，其表达式如下：x1(t+1)＝x1(t)+Tx2(t)，x2(t+1)＝ρx2(t)+w(t)，其中，T表示采样周期，x1(t)和x2(t)分别表示t时刻目标的运动目标的位置和速度，ρ表示目标在t时刻和t+1时刻的速度之间的加速度率，w(t)表示系统噪声；设向量X(t)＝[x1(t),x2(t)]T，则有，定义设定雷达目标跟踪系统的观测方程表达式，其表达式如下：y(t)＝[10]X(t)+v(t)，其中，y(t)是雷达对于目标位置的观测信号，观测信号的误差用噪声形式表示为v(t)。所述建立虚拟噪声模型包括以下步骤：假设w(t)和v(t)是零均值、方差各为和的独立的高斯白噪声，ρ为加速度率，设ρ是一个围绕某个值ρ0变化的变量，表示为ρ＝ρ0+m(t)，m(t)＝n0sin(kt)，其中，ρ0为常数，m(t)为随时间变化的扰动变量，将m(t)设为一个正弦函数，n0和k为常数；雷达目标跟踪子系统扰动表示为其中系数阵分别记作：设虚拟噪声如下：ξ(t)＝m(t)x(t)+w(t)，其中，虚拟噪声ξ(t)为一随机序列，即将系统噪声w(t)和m(t)x(t)的和作为一个虚拟噪声；根据雷达目标跟踪子系统扰动方程和观测方程，通过矩阵的求逆和转化过程可得，y(t)＝H(In-q-1Φe)-1q-1ξ(t)+v(t)，其中，q-1为延迟因子；H(In-q-1Φe)-1q-1＝A-1(q-1)B(q-1)；其中，A(q-1)和B(q-1)都是多项式，具有形式X(q-1)＝1+x1q-1+x2q-2+...+xnq-n；y(t)公式转化为：A(q-1)y(t)＝B(q-1)ξ(t)+A(q-1)v(t)；引入r(t)将已知的观测y(t)转化为多项式序列的形式，其形式如下：r(t)＝A(q-1)y(t)。所述估值计算包括以下步骤：r(t)函数递归形式的转化，其转化形式如下：具有一致性，则有：求解R(k)的值，由A(q-1)y(t)的计算公式，在两边求解相关函数，则可得方程组通过该方程组可以求得的估值所述自校正阿尔法贝塔滤波器的滤波包括以下步骤：将该方差替换动态噪声w(t)的方差带入至α-β滤波器，即将虚拟噪声替代真实噪声带入到稳态kalman滤波公式中，其表示如下：Ψf＝(I2-KfH)Φ，则可求得自校正α-β滤波器；自校正α-β滤波器获取雷达目标跟踪系统扰动目标，对其进行滤波处理。有益效果本专利技术的一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，与现有技术相比利用拟合虚拟噪声的方法建立一个虚拟噪声模型，用虚拟噪声代替扰动和系统本身的动态噪声(虚拟噪声的方差通过建立相关函数方程的方法求解出对应估值)，将参数带扰动的系统转化为带稳定的未知参数和已知噪声的滤波问题，生成自校正α-β滤波器，利用自校正α-β滤波器对未知目标进行有效滤波，保证了滤波精度。附图说明图1为本专利技术的方法顺序图；图2为利用本专利技术方法对位置变量x1(t)的滤波结果和真实值的对比图；图3为利用本专利技术方法对速度变量x2(t)的滤波结果和真实值的对比图。具体实施方式为使对本专利技术的结构特征及所达成的功效有更进一步的了解与认识，用以较佳的实施例及附图配合详细的说明，说明如下：如图1所示，本专利技术所述的一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，包括以下步骤：第一步，雷达目标跟踪系统的设定，对雷达目标跟踪系统进行预设定。其具体步骤如下：(1)设定雷达目标跟踪系统表达式，其表达式如下：x1(t+1)＝x1(t)+Tx2(t)，x2(t+1)＝ρx2(t)+w(t)，其中，T表示采样周期，x1(t)和x2(t)分别表示t时刻目标的运动目标的位置和速度，ρ表示目标在t时刻(T周期)和t+1时刻(T+1周期)的速度之间的加速度率，w(t)表示系统噪声。(2)为了方便线性计算，在此，设向量X(t)＝[x1(t),x2(t)]T，则有，(3)为了方便线性运算，便于在后边的公式转化中使用，在此为矩阵和向量设定标号。定义(4)设定雷达目标跟踪系统的观测方程表达式，其表达式如下：y(t)＝[10]X(t)+v(t)，其中，y(t)是雷达对于目标位置的观测信号，由于测量存在误差，该误差可以用噪声形式表示，即观测信号的误差用噪声形式表示为v(t)。第二步，建立虚拟噪声模型。拟合虚拟噪声，用虚拟噪声代替扰动和雷达目标跟踪系统的动态噪声，其具体步骤如下：(1)假设w(t)和v(t)是零均值、方差各为和的独立的高斯白噪声，ρ为加速度率，设ρ是一个围绕某个值ρ0变化的变量，可表示为ρ＝ρ0+m(t)，m(t)＝n0sin(kt)，其中，ρ0为常数，m(t)为随时间变化的扰动变量，将m(t)设为一个正弦函数，n0和k为常数。(2)雷达目标跟踪子系统扰动表示为其中系数阵分别记作：(3)设虚拟噪声如下：ξ(t)＝m(t)x(t)+w(t)，其中，虚拟噪声ξ(t)为一随机序列，即将系统噪声w(t)和m(t)x(t)的和作为一个虚拟噪声，以方便计算。(4)根据雷达目标跟踪子系统扰动方程和观测方程，通过矩阵的求逆和转化过程可得，y(t)＝H(In-q-1Φe)-1q-1ξ(t)+v(t)，其中，q-1为延迟因子；H(In-q-1Φe)-1q-1＝A-1(q-1)B(q-1)，在此引入两个多项式，其中，A(q-1)和B(q-1)都是多项式，具有形式X(q-1)＝1+x1q-1+x2q-2+...+xnq-n；y(t)公式转化为：A(q-1)y(t)＝B(q-1)ξ(t)+A(q-1)v(t)。(5)引入r(t)将已知的观测y(t)转化为多项式序列的形式，其形式如下：r(t)＝A(q-1)y(t)。第三步，估值计算。将虚拟噪声的方差通过建立函数方程的方法求解出对本文档来自技高网...

【技术保护点】
1.一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，其特征在于，包括以下步骤：11)雷达目标跟踪系统的设定，对雷达目标跟踪系统进行预设定；12)建立虚拟噪声模型，拟合虚拟噪声，用虚拟噪声代替扰动和雷达目标跟踪系统的动态噪声；13)估值计算，将虚拟噪声的方差通过建立函数方程的方法求解出对应估算；14)自校正阿尔法贝塔滤波器的滤波，将虚拟噪声的方差替换自阿尔法贝塔滤波器动态噪声的方差，生成自校正阿尔法贝塔滤波器，通过自校正阿尔法贝塔滤波器进行雷达目标跟踪系统扰动目标的滤波处理。

【技术特征摘要】
1.一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，其特征在于，包括以下步骤：11)雷达目标跟踪系统的设定，对雷达目标跟踪系统进行预设定；12)建立虚拟噪声模型，拟合虚拟噪声，用虚拟噪声代替扰动和雷达目标跟踪系统的动态噪声；13)估值计算，将虚拟噪声的方差通过建立函数方程的方法求解出对应估算；14)自校正阿尔法贝塔滤波器的滤波，将虚拟噪声的方差替换自阿尔法贝塔滤波器动态噪声的方差，生成自校正阿尔法贝塔滤波器，通过自校正阿尔法贝塔滤波器进行雷达目标跟踪系统扰动目标的滤波处理。2.根据权利要求1所述的一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，其特征在于，所述雷达目标跟踪系统的设定包括以下步骤：21)设定雷达目标跟踪系统表达式，其表达式如下：x1(t+1)＝x1(t)+Tx2(t)，x2(t+1)＝ρx2(t)+w(t)，其中，T表示采样周期，x1(t)和x2(t)分别表示t时刻目标的运动目标的位置和速度，ρ表示目标在t时刻和t+1时刻的速度之间的加速度率，w(t)表示系统噪声；22)设向量X(t)＝[x1(t),x2(t)]T，则有，23)定义24)设定雷达目标跟踪系统的观测方程表达式，其表达式如下：y(t)＝[10]X(t)+v(t)，其中，y(t)是雷达对于目标位置的观测信号，观测信号的误差用噪声形式表示为v(t)。3.根据权利要求1所述的一种基于阿尔法贝塔滤波器的雷达目标跟踪系统扰动目标滤波方法，其特征在于，所述建立虚拟噪声模型包括以下步骤：31)假设w(t)和v(t)是零均值、方差各为和的独立的高斯白噪声，ρ为加速度率，设ρ是一个围绕某个值ρ0变化的变量，表示为ρ＝ρ0+m(t)，m(t)＝n0sin(kt)，其中，ρ0为常数，m(t)为随时间变...

【专利技术属性】
技术研发人员：李恒，张晨晨，李佩君，王诗兵，刘志远，史晓凤，尹佳来，徐小东，周正港，张炎炎，李啊敏，刘飞彪，陈佳琪，袁冬冬，
申请(专利权)人：阜阳师范学院，
类型：发明
国别省市：安徽,34