Unbiased converted measurement manoeuvering target tracking under maximum correntropy criterion

Cognitive Computation and SystemsVolume 2, Issue 3 p. 125-129 Research ArticleOpen Access Unbiased converted measurement manoeuvering target tracking under maximum correntropy criterion Guoyong Wang, Corresponding Author standardwgy@163.com School of Computer and Information Engineering, Luoyang Institute of Science and Technology, Luoyang, 471000 People's Republic of ChinaSearch for more papers by this authorXiaoliang Feng, College of Electrical Engineering, Henan University of Technology, Zhengzhou, 450001 People's Republic of ChinaSearch for more papers by this author Guoyong Wang, Corresponding Author standardwgy@163.com School of Computer and Information Engineering, Luoyang Institute of Science and Technology, Luoyang, 471000 People's Republic of ChinaSearch for more papers by this authorXiaoliang Feng, College of Electrical Engineering, Henan University of Technology, Zhengzhou, 450001 People's Republic of ChinaSearch for more papers by this author First published: 14 September 2020 https://doi.org/10.1049/ccs.2020.0010AboutSectionsPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onEmailFacebookTwitterLinked InRedditWechat Abstract In this study, the manoeuvering target tracking problem is addressed by using the unbiased converted measurements from a two-dimensional radar system. Due to the fact that radar measurements are usually expressed in polar coordinates while the target motion is described in the Cartesian coordinates, the unbiased converted measurements are utilised to linearise the system model of the manoeuvering target tracking problem in the Cartesian coordinates. The manoeuver acceleration is modelled as the unknown input of the constant velocity kinematic model of the target. First, it is pointed out that the converted measurement noise no longer satisfies Gaussian distribution, even if the raw radar measurement noise is Gaussian noise. In order to solve the manoeuvering target tracking problem with non-Gaussian disturbances, a joint estimation method for the target state and the unknown input is presented under the maximum correntropy criterion. In the simulation, the proposed manoeuvering target tracking method is compared with the one developed on the basis of the traditional Kalman filter. The simulation results verify the effectiveness of the method proposed in this study. 1 Introduction The radar system is a widely used target tracking measurement system in various military and civilian fields [1-4]. As well known radar measurement is usually expressed in polar coordinates, while the target kinematic model is described in the Cartesian coordinates. They are usually converted into the same coordinate system. However, the non-linear conversion possibly makes the Gaussian noise no longer obey Gaussian distribution. It is one of the important issues that need to be solved for target tracking with radar platforms. The target tracking problem based on the radar platform has received extensive attention from researchers, and some important results have been achieved. The existing results mainly include two categories. (i) Re-express the radial distance and the pitch angle in radar measurement with the position of the target, then the target-tracking problems are considered as non-linear filtering problems. On this basis, with the help of extended Kalman filter (EKF), unscented Kalman filter (UKF) and other non-linear filtering methods, the target states were estimated in various target-tracking scenarios [5-8]. (ii) Convert the radial distance and the pitch angle of the radar measurement into the measurement of the position of the target, then the target tracking problems are constructed as linear filtering problems [9-12]. In this context, the famous Kalman filter was utilised to estimate the state of the target. However, the converted measurement noise no longer satisfies Gaussian distribution, after being transformed by a non-linear function. The filtering method based on Gaussian system may reduce the target tracking accuracy. On the other hand, the manoeuvering target tracking problems have also attracted considerable attention in recent years [13-19]. The manoeuvering target tracking methods were studied on the basis of different target manoeuvering models. Based on the Singer model, a manoeuvering target tracking method was presented, in [13], by using the Kalman filter. With the help of three-dimensional constant-turn models, the EKF and the UKF were used to design the manoeuvering target-tracking methods in [14]. By modelling the manoeuvering target with multiple models, the interacting multiple model approaches were studied for the manoeuvering target tracking problems, by considering discrete-time hybrid systems with Gaussian system noises in [15-17]. Taking the target manoeuvering acceleration as the unknown input, the target state and the manoeuvering acceleration were simultaneously estimated in [18, 19]. A common characteristic of these above-mentioned works is that the measurement noises were assumed satisfying Gaussian distribution. However, the converted measurement noise no longer satisfies Gaussian distribution. In addition, the manoeuvering target tracking research with non-Gaussian noise has been seldom reported. In this paper, a manoeuvering target tracking method is studied with the converted measurements from a two-dimensional radar system. The radar measurement in polar coordinates is firstly converted into an unbiased conversion measurement in the Cartesian coordinates. It is shown that the converted measurement noises no longer obey the Gaussian distribution, even if the raw radar measurement noise is Gaussian noise. The target kinematic equation is modelled as a constant velocity motion with an unknown manoeuvering acceleration. Then, with the help of the information-theoretic learning theory [20], the target state and the unknown manoeuvering acceleration are estimated, simultaneously, under the maximum correntropy criterion. The effectiveness of the proposed method is illustrated by numerical simulation, on comparing with the method developed on the basis of the Kalman filter. The main contributions of this paper are two-fold. (i) It is pointed out that the converted measurement noises no longer obey the Gaussian distribution, even if the raw radar measurement noise is Gaussian noise, due to the non-linear measurement conversion. (ii) A non-Gaussian joint state and unknown input estimation method are presented under the maximum correntropy criterion, to estimate the target state and the unknown manoeuvering acceleration. The rest of this paper is organised as follows. Section 2 introduces the manoeuvering target tracking system with unbiased converted measurement and points out that the converted measurement noises no longer satisfy the Gaussian distribution. In Section 3, the target state and the unknown manoeuvering acceleration are estimated under the maximum correntropy criterion. Simulation is provided, in Section 4, to illustrate the effectiveness of the proposed method. Section 5 concludes this paper. 2 Problem formulation Consider the following discrete manoeuvering target tracking model (1)where represents the state of the target at k, including its positions and velocities on the x and y axes. Namely, is the unknown manoeuvering acceleration, the evolution equation of which is given by (2)In (1), and are, respectively, the state transition matrix and the manoeuvering acceleration transition matrix (3)The process noise and satisfy the zero-mean Gaussian distributions with the covariance and . T is the discretisation period. The target is observed by a two-dimensional radar platform. As well known, the measurement of the radar platform is obtained in the polar coordinates. In addition, the measurement function is modelled as follows in this paper (4)where is the radial distance between the radar system and the target and is its measurement obtained by the radar system. is the pitch angle between the radar system and the target, is its measurement, as well as the noises . Converting the measurement in the polar coordinates into the one in the Cartesian coordinates, it yields (5) (6)Using the trigonometric formula, we can obtain (7)in which [21] (8) (9) Remark 1.From (8) and (9), it is noted that the converted measurement noises and no longer obey the Gaussian distribution, even if the raw radar measurement noises are Gaussian noise. Under the condition of the first k measurements (10)and (11)where (12) (13) (14)Denote the augmented state and the augmented process noise as follows: (15)It yields (16)where and (17)The unbiased converted measurement can be further re-written as (18)where (19)As mentioned in the Remark 1, the measurement noises in (8) and (9) do not satisfy the Gaussian distribution, even if are Gaussian noises. Therefore, the estimator of the target state and the manoeuvering acceleration with the above unbiased converted measurements should not be designed on the basis of the famous Kalman filter and its evolution methods, which are proposed for the systems with Gaussian noises. In the next section, without assuming the converted measurement noises satisfy Gaussian distribution, the manoeuvering target tracking approach is studied under the maximum correntropy criterion. 3 Manoeuvering target tracking under maximum correntropy criterion In this section, the manoeuvering target tracking approach is studied with unbiased converted measurement, under the maximum correntropy criterion. Without loss of generality, assume the augmented state estimate and its estimation error covariance at are and . According to (16), the augmented state of the target at k can be predicted by (20)and (21)Combining (20) and (18), it yields (22)with (23)Under the maximum correntropy criterion, the cost function is given by [20, 22] (24)where are the ith row of , L is the dimension of , and is the Gaussian kernel which is given by (25)and (26)Therefore, the optimal estimate of the augmented state is the solution of (27)Thus the optimal estimate is the one satisfying (28)Further, it yields (29)where (30)with Substituting (26) and (30) into (29), we have (31)in which, (32)It should be noted that the augmented state is included in , as shown in (29) and (30). Namely, (31) is a fixed point function with respect to . Taking its prediction as the initial solution , the iterative solution of the fixed point function results into a iterative non-Gaussian filtering algorithm [22]. The tth iterative updating augmented state estimate is (33)The filter gain is obtained by (34)The estimation error covariance is given by (35)Given a small scalar , the final augmented state estimate is obtained by , if the following inequality holds (36)The final augmented state estimate includes the manoeuvering acceleration estimate as its 5th,6th components, and the estimate of the target position and velocity as its first four components. Remark 2.For the method developed from the Kalman filter, the converted measurement noise is assumed to obey the zero-mean Gaussian distribution. Denote the noise variance as . The augmented system state is estimated by the following method (37) where (38) The estimation error covariance is obtained by (39) 4 Simulation Consider the following manoeuvering target tracking system. (40) (41) (42)where The idea position of the radar platform is . The initial estimate of the target manoeuver value is . The initial state of the target is and the initial estimate covariance is . In this simulation, the proposed method is compared with the one developed from Kalman filter as shown in Remark 2. The proposed method is denoted as manoeuvering target maximum correntropy Kalman filter (MTMCKF)-based method, and the one developed from Kalman filter is denoted as KF-based method. The simulation results are shown in Table 1, Figs. 1 and 2. Table 1. Mean of absolution estimation error MTMCKF-based method KF based method MAEE 0.18811 0.22525 As shown in Figs. 1 and 2, the effectiveness of the proposed method is illustrated. In this simulation, the radar measurements on the polar coordinates are converted into the one on the Cartesian coordinates. The converted measurement noise is taken as Gaussian noise in the method developed from Kalman filter, but is not in the proposed tracking method. Therefore, the distance between the real target and the estimate obtained by the proposed method is smaller than the one obtained by the Kalman filter based method. The mean of absolution estimation error of the proposed method is 0.18811, which improves 16.5%, compared with the one obtained by the method developed from Kalman filter, whose mean of absolution estimation error (MAEE) is 0.22525. Fig. 1Open in figure viewerPowerPoint Position curves of the target and its estimates Fig. 2Open in figure viewerPowerPoint Distance curves between the target and its estimates 5 Conclusion In this paper, the manoeuvering target tracking problem is addressed based on the unbiased conversion of the radar measurement. The converted measurement noise no longer obeys the Gaussian distribution, even if the raw measurement noise is Gaussian noise. Therefore, in this paper, a manoeuvering target tracking method is developed for the non-Gaussian measurements. The manoeuvering target tracking system is modelled as a CV motion with the target manoeuvering acceleration which is taken as an unknown input. Under the maximum correntropy criterion, a state and unknown input joint estimation approach is given to solve the manoeuvering target tracking problem. The final simulation illustrates the effectiveness of the proposed approach. How to avoid the singular matrix that appears in the Cholesky decomposition is an important direction in our next work. The square root filtering methods give us meaningful inspiration. 6 Acknowledgments This work was supported in part by the National Natural Science Foundation of China (U1604149). 7 References 1Weinmann, F.: 'EM modelling of radar signatures of targets behind wind farms–A time-gating ray tracking approach'. Proc. of 2016 European Radar Conf., 2016, pp. 9– 12Google Scholar 2Zhang, G., Liang, J., Zhao, H. et al.: 'Sequential Monte Carlo Implementation for infrared/radar maneuvering target tracking'. Proc. of Sixth World Congress on Intelligent Control and Automation, Dalian, China, 2006, pp. 5066– 5069Google Scholar 3Levita, B., Drozdov, M., Naidionova, I. et al.: 'MIMO UWB radar for moving target tracking'. Proc. of 17th Int. Radar Symp., Krakow, Poland, 2016, pp. 1– 4Google Scholar 4Mao, Y., Li, X., Duan, Z. et al.: 'Unbiased measurement model conversion for tracking with multiple radars or sonars'. Proc. of the 31st Chinese Control Conf., Hefei, China, 2012, pp. 897– 903Google Scholar 5Peng, Z., Xia, H., Lu, H. et al.: 'Adaptive robust high-degree cubature filter based on target tracking'. Proc. of Third Advanced Information Technology, Electronic and Automation Control Conf., Chongqing, China, 2018, pp. 552– 557Google Scholar 6Li, B., Qu, C., Zhou, Z. et al.: 'Application of UKF alogrithm in passive localization and tracking'. Proc. of Third Advanced Information Technology, Electronic and Automation Control Conf., Chongqing, China, 2018, pp. 2530– 2534Google Scholar 7Ge, Q., Shao, T., Duan, Z. et al.: 'Performance analysis of the Kalman filter with mismatched noise covariances', IEEE Trans. Autom. Control., 2016, 61, (12), pp. 4014– 4019CrossrefWeb of Science®Google Scholar 8Ge, Q., Shao, T., Yang, Q. et al.: 'Multisensor nonlinear fusion methods based on adaptive ensemble fifth-degree iterated cubature information filter for biome-chatronics', IEEE Trans. Syst. Man Cybern. Syst., 2016, 46, (7), pp. 912– 925CrossrefWeb of Science®Google Scholar 9Mo, L., Song, X., Zhou, Y. et al.: 'Unbiased converted measurements for tracking', IEEE Trans. Aerosp. Electron. Syst., 1998, 34, (3), pp. 1023– 1027CrossrefWeb of Science®Google Scholar 10Duan, Z., Han, C., Li, X.: 'Comments on unbiased converted measurements for tracking', IEEE Trans. Aerosp. Electron. Syst., 2004, 40, (4), pp. 1374– 1377CrossrefWeb of Science®Google Scholar 11Bordonaro, S., Willett, P., Bar-Shalom, Y.: 'Decorrelated unbiased converted measurement Kalman filter', IEEE Trans. Aerosp. Electron. Syst., 2014, 50, (2), pp. 1431– 1444CrossrefWeb of Science®Google Scholar 12Liu, H., Zhou, Z., Yu, L. et al.: 'Two unbiased converted measurement Kalman filtering algorithms with range rate', IET Radar Sonar Navig., 2018, 12, (11), pp. 1217– 1224Wiley Online LibraryWeb of Science®Google Scholar 13Singer, R.A.: 'Estimating optimal tracking filter performance for manned maneuvering targets', IEEE Trans. Aerosp. Electron. Syst., 1970, 6, (4), pp. 473– 483CrossrefWeb of Science®Google Scholar 14Zhan, R., Wan, J.: 'Passive maneuvering target tracking using 3D constant-turn model'. Proc. of 2006 IEEE Conf. on Radar, NY, USA, 2006, DOI:10.1109/RADAR.2006.1631832Google Scholar 15Nadarajah, N., Tharmarasa, R., McDonald, M. et al.: 'IMM forward filtering and backward smoothing for maneuvering target tracking', IEEE Trans. Aerosp. Electron. Syst., 2012, 48, (3), pp. 2673– 2678CrossrefWeb of Science®Google Scholar 16Han, B., Huang, H., Lei, L. et al.: 'An improved IMM algorithm based on STSRCKF for maneuvering target tracking', IEEE Access, 2019, 7, pp. 57795– 57804CrossrefWeb of Science®Google Scholar 17Li, W., Jia, Y.: 'Location of mobile station with maneuvers using an IMM-based cubature Kalman filter', IEEE Trans. Ind. Electron., 2012, 59, (11), pp. 4338– 4348CrossrefWeb of Science®Google Scholar 18Hu, S., Yang, J., Zeng, F. et al. 'Interacting multiple model tracking algorithm with modified input estimation', Acta Electron. Sin., 2009, 37, (12), pp. 2810– 2814Google Scholar 19Lan, H., Liang, Y., Yang, F. et al.: 'Joint estimation and identification for stochastic systems with unknown inputs', IET Control Theory Applic., 2013, 7, (10), pp. 1377– 1386Wiley Online LibraryWeb of Science®Google Scholar 20Principe, J.C.: ' Information theoretic learning: renyis entropy and kernel perspectives' ( Springer: New York, USA, 2010) CrossrefGoogle Scholar 21Han, C., Zhu, H., Duan, Z.: ' Multi-source information fusion' ( Tsinghua University Press, Beijing, 2006) Google Scholar 22Chen, B., Liu, X., Zhao, H. et al.: 'Maximum correntropy Kalman filter', Automatica, 2017, 76, pp. 70– 77CrossrefWeb of Science®Google Scholar Volume2, Issue3September 2020Pages 125-129 FiguresReferencesRelatedInformation