![雷达数据处理及应用(第四版)](https://wfqqreader-1252317822.image.myqcloud.com/cover/900/47379900/b_47379900.jpg)
3.2.3 卡尔曼滤波的初始化
本节讨论状态估计的初始化问题是运用卡尔曼滤波的一个重要的前提条件,只有进行了初始化,才能利用卡尔曼滤波对目标进行跟踪。
1.二维状态向量估计的初始化
系统的状态方程和量测方程同式(3.31)、式(3.32),此时的状态向量表示为,量测噪声W(k)~N(0,r),且与过程噪声相互独立。这种情况下的状态估计初始化可采用两点差分法,该方法只利用第一和第二时刻的两个量测值Z(0)和Z(1)进行初始化,即初始状态为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_70_3.jpg?sign=1738883662-NWs323N6AszL2tWHjZtXxH92FRBVOeeg-0-fa936a70c2189ecb8d72fbb2e2d43d37)
式中,T为采样间隔。初始协方差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_70_4.jpg?sign=1738883662-nr5E73N4D1Oov18MpvnrAXAfq1QUbCs7-0-75906fae944ca33fa655d41b75040a14)
于是状态估计和滤波从k=2时刻开始。
在对算法进行多次Monte Carlo试验时,则在每次试验中都必须重新产生新的噪声,接着再使用同样的方法初始化。在Monte Carlo试验中,重复使用同样的初始条件将导致有偏估计,所以每次试验初始状态估计应重新随机选择。二维卡尔曼滤波通常在x、y、z轴解耦滤波的情况下使用。
2.四维状态向量估计的初始化
这种情况描述的是两坐标雷达的数据处理问题,此时系统的状态向量若表示为,而直角坐标系下的量测值Z(k)为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_70_6.jpg?sign=1738883662-tjCwvr4kRDAKBVyyrZAsUAVwdBGlioRy-0-cbe152b025a0babae1d5c5d392c24044)
式中,ρ和θ分别为极坐标系下雷达的目标径向距离和方位角量测数据。则系统的初始状态可利用前两个时刻的量值测Z(0)和Z(1)来确定,即
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_70_7.jpg?sign=1738883662-87TbBWoMILZ9PMyNW8phoupviKXT6P45-0-12203f94bcb8fb70cbddeee7c3897001)
k时刻量测噪声在直角坐标系下的协方差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_1.jpg?sign=1738883662-JkJHvIjGDidU0xiUMiMQDswmzkrjtlGO-0-0a1747273484402b9a7fcae6a9ce3211)
式中,和
分别为目标径向距离和方位角测量误差的方差,而
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_4.jpg?sign=1738883662-xzC9sXwWLoMMvmzbSFawTVGy6Z4NqHY5-0-4263f866ef84ee0ea0f508e8b688b26e)
由量测噪声协方差的各元素可得四维状态向量情况下的初始协方差矩阵为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_5.jpg?sign=1738883662-FLpbE8YrfPUbrdLiWDbEtTCAgCaKHqkU-0-be3a16c7f63abafcea4a5c5ee520a412)
并且滤波器从k=2时刻开始工作。
3.六维状态向量估计的初始化
针对三坐标雷达数据处理问题,此时系统状态向量若表示为,而直角坐标系下的量测值Z(k)为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_7.jpg?sign=1738883662-uvVdGKhzQyfqf9YB2kisRFZGz750nra6-0-142c248a4f8d9b42206633e7b5a0ae9c)
式中,ρ和θ的定义同四维状态向量情况,而ε为目标的俯仰角量测数据。此时系统的初始状态仍只需利用前两个时刻的量测值Z(0)和Z(1)来确定,即
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_8.jpg?sign=1738883662-njN7rxZ0iOQigG7PRB64DyrxBxGmNX7z-0-39f1306a529438c88c8203e7fcb808fb)
在这种情况下,k时刻直角坐标系下的量测噪声协方差为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_9.jpg?sign=1738883662-3i31Q4MmhMrJGdilP4aPrI0dn2gdnIpb-0-458dee1deaf6a7f920ed1ce309d32478)
式中,和
的定义同四维状态向量情况,
为俯仰角测量误差的方差,而
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_13.jpg?sign=1738883662-SlWkeWWGW8wkth6KCfyKxBAfyJfOlWoI-0-d91bad3265761ab6e754e45f60a8b1a6)
由量测噪声协方差的各元素可得六维状态向量情况下的初始协方差阵为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_71_14.jpg?sign=1738883662-EBuAVkyKfu3mIeA5KvhrNz9hrYULNEXc-0-700796f3d222783116fe5b8d11c8abb4)
并且滤波器从k=2时刻开始工作。
4.九维状态向量估计的初始化
该情况下系统的状态向量若表示为,它与六维情况相比只是多了加速度项,所以此时直角坐标系下的目标量测值Z(k)、量测噪声协方差R(k)仍和六维情况相同。
由于此时含加速度项,所以系统的初始状态需利用前三个时刻的量测值Z(0)、Z(1)和Z(2)确定,即
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_72_2.jpg?sign=1738883662-UMnw4dwJMkFm3XsLXraIIYsZ3aytoE2A-0-ff334a0f8196ba2ee7e120868993279f)
而初始协方差矩阵为
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_72_3.jpg?sign=1738883662-1H1TO54Ig0DtgbbQKe3DfkEqo6h1E2Wd-0-38ed15867f55ea656b9d654fdaa63b4d)
式中,P11、P12、P13、P22、P23和P33为分块矩阵,且
![](https://epubservercos.yuewen.com/C6293B/26763973309544206/epubprivate/OEBPS/Images/43988_72_4.jpg?sign=1738883662-nvm0xGM7ceScqtEGLeUwdPs3yCyPw8k1-0-e5153bcb799b6e0cdf9a1beca18cddda)
并且滤波器从k=3时刻开始工作。