Abstract
The problem of adaptive estimation using spatialtimevarying filtering in dynamic systems with random structure is investigated. The proposed approach of the extension of state estimation in the classical stochastic dynamic systems with deterministic structure to the case of signal processing and parameter identification in stochastic systems with random structure or with switching parameters using the correlationextremum methods and the theory of Markov processes provides the system operation in varying and uncertain external conditions.
1. Introduction
In many practical cases the state estimate of the observable object coordinates cannot be obtained using the optimal filtering algorithms. In particular, this relates to the case of radar tracking of an airborne target when the vehicle dynamics is described by nonlinear differential equations. The known nonlinear filtering algorithms are based on the assumption of linearization of the nonlinear functions in the state and measurements equations relative to estimation errors. These algorithms remain true when the estimation errors are small enough to satisfy a linearization.
Some previous work has already dealt with the linearization problem in the case of great estimation errors. One of the first approaches was taken in [1], where the twofeedback filtering system with changing operating conditions depending on the estimation error value of the observable object coordinates was suggested. But the estimation method discussed in [1] was not a result of a rigorous solution of such system synthesis problem that did not allow the twofeedback system to be considered as an optimal one.
The other filtering problems to be solved are connected with the moving objects parameters estimation, particularly, in the cases of discontinuous jump of the estimated process features, abrupt increasing of the measurements noise, informative signal interruption, and etc.
In most practical applications of filtering theory the uncertain parts of the dynamic system model are sometimes represented by an unknown parameter vector and the variations of these parameters and their identification play an important role.
Many scientific researches have been performed in the class of adaptive filtering scheme of the state vector estimation together with parameter identification (for example, [2, 3]). In [2] the problem of state estimation for discrete systems with parameters which may be switching within a finite set of values was considered for the general case and for the Markov parameter case. The approach to unknown parameters estimation taken in [3] was based on the numerical differentiation method with ${B}^{0}$ splines to assure better noiseresistance.
In the known earlier researches the theory of Markov processes and systems with random structure has been developed for signal processing when the observable signal has been considered only as a function of time and the correlationextremum methods have not been used. This limitation did not allow the system with random structure to use the information of spatialtimevarying signals of different physical nature (obtained by radar, optical, infrared and other image sensors).
The solution of the abovementioned problems applied to the spatialtimevarying signals processing was originally proposed in [4] where the correlationextremum systems theory was extended to systems with random structure and the new suboptimal filtering algorithms were derived. The case when the change of the structure is supposed to be Markov process with two states is considered in this paper and the solution of the linearization problem based on the derived filtering algorithm is proposed, particularly, for the case of great estimation errors.
There are a number of potential military and civil application areas of the above proposed [4] nonlinear information processing algorithms. In this paper the solution for a target tracking system (TTS) is considered as an example.
The problem of target tracking is under consideration since many previous researches (see, for example, [5], where the treatment was based on a closedform quasilinearized solution of the proportional navigation (PN) equations followed by the smallangle approximation only to line of sight (LOS) angle rate). In [5] only the dynamic system state model without noise and measurements was discussed, and the state estimation problem was not solved. The other solutions of the state estimation problem for target tracking systems are based on Kalman filtering schemes (linear or suboptimal extended one) which can process only timevarying signal functions.
A new approach of the extension of Kalman filter adaptive estimation using the theory of Markov processes and systems with random structure or with switching parameters was proposed in [6] to the case of signal processing and parameter identification in stochastic magnetic navigation systems, where the parameters to be identified are the Poisson's coefficients.
The purpose of this paper is to present the solutions of the problem of nonlinear filtering and identification algorithms synthesis and analysis for correlationextremum systems, based on Markov processes and the theory of optimization of stochastic systems with random structure for military and civil applications.
1.1. Problem statement
The problem under consideration is nonlinear estimation for the dynamic state process described by a stochastic differential equation (Eq. (1)):
where $\mathbf{\Lambda}\left(t\right)$ is the state vector, which contains the random, unknown, and timevarying parameters vector ${\mathit{a}}^{T}=({a}_{1},\dots \mathrm{},{a}_{\mathrm{q}})$, $u\left(t\right)$ is the known control vector, which may depends on the state vector components, $l\left(t\right)$ is a stationary Markov process taking values in the set {1, 2,..., $p$} (number of the state). Here $\mathbf{\Lambda}\left(t\right)$ is $n$dimensional vector, in general case, with initial Gaussian value$\mathbf{\Lambda}{(t}_{0})$, ${\mathbf{F}}^{\left(l\right)}(\mathbf{\Lambda},\mathbf{u},\mathrm{}t)$ is the nonlinear deterministic vector function ${\mathbf{F}}^{\left(l\right)}(\mathbf{\Lambda},\mathbf{u},\mathrm{}t)=\Vert {f}_{i}^{\left(l\right)}(\mathbf{\Lambda},\mathbf{u},\mathrm{}t)\Vert \text{,}$$(i=\overline{1,n})$ which satisfies Lipschitz conditions. ${\mathbf{W}}^{\left(l\right)}\left(t\right)$ is a vector process of the state Gaussian white noise with diagonal intensity matrix ${\mathbf{Q}}^{\left(l\right)}\left(t\right)=\Vert {q}_{i}^{\left(l\right)}\left(t\right)\Vert $, $(i=\overline{1,n})$.
The observable signal $\mathbf{r}\left(x,\mathrm{}y,\mathrm{}t\right)$ is the $m$dimensional spatialtimevarying process described by the following measurement equation (Eq. (2)):
where $x$, $y$ are the space variables – space coordinates at any point – $x\in X=[{x}_{0},{x}_{X}]$, $y\in Y=[{y}_{0},{y}_{Y}]$, $t$ is the time variable $t\in T=[{t}_{0},{t}_{T}]$, ${\mathbf{S}}^{\left(l\right)}\left(x,\mathrm{}y,\mathrm{}\mathbf{\Lambda},t\right)$ is the vector of spatialtimevarying signals of different physical nature, ${\mathbf{N}}^{\left(l\right)}\left(x,\mathrm{}y,\mathrm{}t\right)$ is the spatialtimevarying vector process of the measurements Gaussian white noise with diagonal intensity matrix ${\mathbf{C}}_{0}^{\left(l\right)}$ and correlation function $h(\mathrm{\Delta}x,\mathrm{\Delta}y,\mathrm{\Delta}t)={\mathbf{C}}_{0}^{\left(l\right)}\delta (\mathrm{\Delta}x,\mathrm{\Delta}y,\mathrm{\Delta}t)$, $\delta $ is the delta function.
The signal position on the image plane $XOY$ can be determined by parameters vector:
Then the signal may be represented as ${\mathbf{S}}^{\left(l\right)}(x,\mathrm{}y,\mathrm{}\mathbf{\Lambda},t)={\mathbf{S}}^{\left(l\right)}\left(x{\lambda}_{x},y{\lambda}_{y},t\right)\text{.}$ The suboptimal estimator in this case represents a tracker system. The estimation errors are formed by receiving the derivatives of the crosscorrelation function of the observed spatialtimevarying signal $\mathbf{r}(x,\mathrm{}y,\mathrm{}t)$ and the reference image ${\mathbf{S}}^{\left(l\right)}(x,\mathrm{}y,\mathrm{}\widehat{\mathbf{\Lambda}},t)$.
The changes of the structure are supposed to be Markov process with p finite states and the transition intensities ${\nu}_{jl}\left(t\right)$ and ${\nu}_{l}\left(t\right)$, where $j$, $l=\overline{1,p}$. The behavior of the system may be explicated as follows. The system begins in a particular mode of operation, say $l\left(t\right)=\text{1}$, then at a random time the system jumps to one of the other ($p1$) possible modes of operation and may or may not remain in this state. The model has an interpretation as a problem of a maneuvering target tracking interruption or a case of great estimation errors.
The nonlinear filtering or estimation problem defined above for Eqs. (12) is to find the finitedimensional dynamical system whose output is the best minimum variance estimate of the joint Markov process ${(\mathbf{\Lambda}\left(t\right),l\left(t\right))}^{T}$, for $t\ge 0$ given the observed data $\mathbf{r}(x,y,t)$.
1.2. Solution of the algorithms synthesis problem
The synthesis of the filtering algorithms for correlationextremum systems with random structure [4] was based on the generalized FokkerPlankKolmogorovStratonovich equation for the evolution of joint conditional probability density of the state dynamics $\mathbf{\Lambda}\left(t\right)$ and the system structure $l\left(t\right)$$\widehat{\omega}(\mathbf{\Lambda},l,t)={\widehat{\omega}}^{\left(l\right)}(\mathbf{\Lambda},t)$ (Eq. (3)):
$\frac{1}{2}{\widehat{\omega}}^{\left(l\right)}\left(\mathbf{\Lambda},t\right)\left[\mathrm{}{\mathrm{\Phi}}^{\left(l\right)}\left(\mathbf{\Lambda},\mathbf{r},t\right)\underset{\mathrm{\infty}}{\overset{\mathrm{\infty}}{\int}}{\mathrm{\Phi}}^{\left(l\right)}\left(\mathbf{z},\mathbf{r},t\right){\widehat{\mathrm{\omega}}}^{\left(l\right)}\left(\mathbf{z},t\right)d\mathbf{z}\mathrm{}\right]$
$+\sum _{j=1\left(j\ne l\right)}^{p}\mathrm{}{\nu}_{jl}\left(t\right)\frac{{\widehat{P}}_{j}\left(t\right)}{{\widehat{P}}_{l}\left(t\right)}\left[\mathrm{}{\widehat{\omega}}^{\left(j\right)}\left(\mathbf{\Lambda},t\right){\widehat{\omega}}^{\left(l\right)}\left(\mathbf{\Lambda},t\right)\right],\omega \mathrm{}({\mathbf{\Lambda}}_{0},\mathrm{}{t}_{0}),l=\overline{1,p},$
where ${\widehat{\mathbf{\pi}}}^{\left(l\right)}\left(\mathbf{\Lambda},t\right)$ is the probability density flow vector:
where ${\mathbf{K}}_{1}^{\left(l\right)}\left(\mathbf{\Lambda},t\right)$ is the local rate vector, ${\mathbf{K}}_{2}^{\left(l\right)}\left(\mathbf{\Lambda},t\right)$ is the diffusion matrix, ${\widehat{P}}_{l}\left(t\right)$ is the a posteriori probability of the $l$th state, $\omega \left({\mathbf{\Lambda}}_{0},\mathrm{}{t}_{0}\right)$ is the initial value of probability density of the state dynamics $\mathbf{\Lambda}{(t}_{0})$, ${\mathrm{\Phi}}^{\left(l\right)}(\mathbf{\Lambda},\mathbf{r},t)$ is the derivative of the likelihood function logarithm in the $l$th state ($z$ is the variable $z\in \mathbf{\Lambda}$, the sign ^ means the a posteriori function value).
The a posteriori probability density for the whole dynamics process is determined as:
The suboptimal estimate of Markov process $\mathbf{\Lambda}\left(t\right)$ on the assumption of mean square loss function is the conditional mathematical expectation. The optimal estimate of discrete process $l\left(t\right)$ by the a posteriori probability criterion will be such a value of $l$ that makes the value of ${\widehat{P}}_{l}\left(t\right)$ maximum.
The estimate of the state is the probabilistically weighted average $\widehat{\mathbf{\Lambda}}\left(t\right)=\sum _{l=1}^{p}{\widehat{P}}_{l}\left(t\right){\widehat{\mathbf{\Lambda}}}^{\left(l\right)}\left(t\right)$.
The main distinction of signal processing in systems with the random structure consists in the presence of the state probability estimate equations (Eq. (4)) in the filtering algorithms and in the relation between these equations:
$+\frac{1}{2}{\widehat{P}}_{l}\left(t\right)\left\{\underset{\mathrm{\infty}}{\overset{\mathrm{\infty}}{\int}}{\mathrm{\Phi}}^{\left(l\right)}\left(\mathbf{z},\mathbf{r},t\right){\widehat{\omega}}^{\left(l\right)}\left(\mathbf{z},t\right)d\mathbf{z}\mathrm{}\mathrm{}\sum _{k=1}^{p}{\widehat{P}}_{k}\left(t\right)\underset{\mathrm{\infty}}{\overset{\mathrm{\infty}}{\int}}{\mathrm{\Phi}}^{\left(k\right)}\left(\mathbf{z},\mathbf{r},t\right){\widehat{\omega}}^{\left(k\right)}\left(\mathbf{z},t\right)d\mathbf{z}\right\}.$
The a priori state probabilities ${P}_{l}\left(t\right)$ can be found according to Kolmogorov equations.
In [4] the solution of the problem was obtained in Gaussian approximation of the a posteriori probability density. In the algorithms synthesis it was supposed that the signal or image position along one of the axes (e.g., in $y$ direction) was known to simplify the derivation, denoting the signal ${\mathbf{S}}^{\left(l\right)}(x,\mathrm{}\mathbf{\Lambda},t)={\mathbf{S}}^{\left(l\right)}(x{\lambda}_{x},t)$, and the state parameter ${\lambda}_{x}$ without index ${\lambda}_{x}=\lambda $.
For the case when the parameters changes form Markov process with two states ($l=\overline{1,2}$) and the transition intensity $\nu \left(t\right)$ the solution of nonlinear filtering problem for systems with random structure derived in [4] is presented by the following correlationextremum algorithms for computing the a posteriori probabilities of state (Eq. (5)), the state estimates (Eq. (6)), and the covariance (Eq. (7)).
The differential equation for the a posteriori probabilities of state is:
$+\left\{\nu \frac{{\widehat{P}}_{1}\left(t\right)}{{C}_{X}^{\left(2\right)}}\left[{k}_{2}+{B}^{\left(2\right)}\left({\mathrm{\Delta}\lambda}_{\left(2\right)}\right)+\frac{1}{2}{\sigma}_{\left(2\right)}^{2}\left(t\right)\frac{{\partial}^{2}{B}^{\left(2\right)}\left({\mathrm{\Delta}\lambda}_{\left(2\right)}\right)}{\partial {\mathrm{\Delta}\lambda}_{\left(2\right)}^{2}}\right]\right\}\left[1{\widehat{P}}_{1}\left(t\right)\right],$
${\widehat{P}}_{2}\left(t\right)=1{\widehat{P}}_{1}\left(t\right),$
where ${\widehat{P}}_{2}\left(t\right)$ is the a posteriori probability of the second state, ${\mathrm{\Delta}\lambda}_{\left(l\right)}\left(t\right)$ is the state estimate error ${\mathrm{\Delta}\lambda}_{\left(l\right)}\left(t\right)=\lambda \left(t\right){\widehat{\lambda}}^{\left(l\right)}\left(t\right)$, ${\sigma}_{\left(l\right)}^{2}\left(t\right)$ is the variance of the a posteriori probability density function ${\sigma}_{\left(l\right)}^{2}\left(t\right)=\u2329{\left[\right(\lambda \left(t\right){\widehat{\lambda}}^{\left(l\right)}\left(t\right)]}^{2}\u232a$, $l=\overline{1,2}$; ${B}^{\left(l\right)}\left({\mathrm{\Delta}\lambda}_{\left(l\right)},t\right)$ is the spatial correlation function in the $l$th state ${B}^{\left(l\right)}({\mathrm{\Delta}\lambda}_{\left(l\right)},t)=\langle {{\mathbf{S}}^{\left(l\right)}}^{T}(x{\widehat{\lambda}}^{\left(l\right)},t){\mathbf{S}}^{\left(l\right)}(x\lambda ,t)\rangle $, ${C}_{X}^{\left(l\right)}$ is the specific spectral intensity ${C}_{X}^{\left(l\right)}={C}_{0}^{\left(l\right)}/X$ of the spatialtimevarying noise ${\mathbf{N}}^{\left(l\right)}(x,\mathrm{}t)$.
The Eq. (5) has been derived using the assumption of the “unpowered” parameters. This assumption means that the integrals $\underset{X}{\overset{X}{\int}}[{{S}^{\left(l\right)}(x{\widehat{\lambda}}^{\left(l\right)},t)]}^{2}dx$ and $\underset{X}{\overset{X}{\int}}{r}^{2}(x,t)dx\text{,}$ which represent signal energy and are explicitly independent of the estimate parameter, may be included in the ${k}_{1}$ and ${k}_{2}$ coefficients.
The state estimate equation (Eq. (6)) is presented below:
$+\nu \frac{{\widehat{P}}_{j}\left(t\right)}{{\widehat{P}}_{l}\left(t\right)}\left[{\widehat{\lambda}}^{\left(l\right)}\left(t\right){\widehat{\lambda}}^{\left(j\right)}\left(t\right)\right],{\widehat{\lambda}}^{\left(l\right)}\left({t}_{0}\right)={\widehat{\lambda}}_{0}^{\left(l\right)},\left(l,j=\overline{\mathrm{1,2}},j\ne l\right),$
where ${N}_{X}^{\left(l\right)}={\int}_{X}^{X}\frac{\partial {{\mathit{S}}^{\left(l\right)}}^{T}\left(x\u2013{\widehat{\lambda}}^{\left(l\right)},t\right)}{\partial {\widehat{\lambda}}^{\left(l\right)}}{\mathbf{N}}^{\left(l\right)}(x,\mathrm{}t)dx$, (as a remark: in many cases the measurements signals and noises are uncorrelated). The suboptimal state estimate for the whole process can be obtained by using a weighted sum $\widehat{\lambda}\left(t\right)={\widehat{P}}_{1}\left(t\right){\widehat{\lambda}}^{\left(1\right)}\left(t\right)+{\widehat{P}}_{2}\left(t\right){\widehat{\lambda}}^{\left(2\right)}\left(t\right)$.
The variance equation is:
$+\nu \frac{{\widehat{P}}_{j}\left(t\right)}{{\widehat{P}}_{l}\left(t\right)}\left[{\sigma}_{\left(j\right)}^{2}\left(t\right){\sigma}_{\left(l\right)}^{2}\left(t\right)+{\left({\widehat{\lambda}}^{\left(j\right)}\left(t\right){\widehat{\lambda}}^{\left(l\right)}\left(t\right)\right)}^{2}\right],{\sigma}_{\left(l\right)}^{2}\left({t}_{0}\right),\left(l,j=\overline{\mathrm{1,2}},j\ne l\right),$
and the variance of the estimate error for the whole process is ${\sigma}^{2}\left(t\right)={\widehat{P}}_{1}\left(t\right){\sigma}_{\left(1\right)}^{2}\left(t\right)+{\widehat{P}}_{2}\left(t\right){\sigma}_{\left(2\right)}^{2}\left(t\right)$. The correlationextremum filtering algorithm with random structure for estimation of signal position along the $y$ axis has been derived, similarly.
The proposed algorithms (Eqs. (57)) can be applied to the solution of the linearization problem in the case of great estimation errors and provide the estimator adaptive capability due to the extension of the correlationextremum filtering and identification methods for systems with random structure.
The research is continuing in different areas of the proposed algorithms application and, particularly, for the systems with interrupted signal information, and for the case when the measurement noise vector ${\mathbf{N}}^{\left(l\right)}(x,\mathrm{}y,\mathrm{}t)$ is the spatialtimevarying colored GaussianMarkov process new algorithms have been derived.
1.3. Applications to TTS
A guidance kinematic scheme under consideration can be described by a system of nonlinear stochastic differential equations (Eq. (8)) for the state dynamics model:
${\dot{D}}^{\left(l\right)}\left(t\right)={V}_{c}\left(t\right)\mathrm{c}\mathrm{o}\mathrm{s}\left({\epsilon}^{\left(l\right)}\left(t\right){\theta}_{c}^{\left(l\right)}\left(t\right)\right)+{V}_{T}\left(t\right)\mathrm{c}\mathrm{o}\mathrm{s}\left({\epsilon}^{\left(l\right)}\left(t\right){\theta}_{T}^{\left(l\right)}\left(t\right)\right)+{W}_{D}\left(t\right),$
${\dot{\theta}}_{c}^{\left(l\right)}\left(t\right)={N}_{c}^{\left(l\right)}\dot{\widehat{\epsilon}}\left(t\right),{\epsilon}^{\left(l\right)}\left({t}_{0}\right)={\epsilon}_{0}^{\left(l\right)},{D}^{\left(l\right)}\left({t}_{0}\right)={D}_{0}^{\left(l\right)},{\theta}_{c}^{\left(l\right)}\left({t}_{0}\right)={\theta}_{c0}^{\left(l\right)},$
where $\epsilon \left(t\right)$ is a LOS angle relative to reference line, ${\theta}_{c}\left(t\right)$ is an angular position of an aircraft (pursuer) velocity vector ${V}_{c}\left(t\right)$, ${\theta}_{T}\left(t\right)$ is an angular position of a target velocity vector ${V}_{T}\left(t\right)$, $D\left(t\right)$ is a current range from an aircraft to a target, ${W}_{\epsilon}\left(t\right)$ and ${W}_{D}\left(t\right)$ are the system uncorrelated zeromean Gaussian white noises with spectral densities ${q}_{\epsilon}\left(t\right)$ and ${q}_{D}\left(t\right)$, respectively; ${\epsilon}^{\left(l\right)}\left({t}_{0}\right)$, ${D}^{\left(l\right)}\left({t}_{0}\right)$, ${\theta}_{c}^{\left(l\right)}\left({t}_{0}\right)$, ${V}_{c}\left({t}_{0}\right)$, ${V}_{T}\left({t}_{0}\right)$, ${\theta}_{T}^{\left(l\right)}\left({t}_{0}\right)$ are the initial values; ${N}_{c}^{\left(l\right)}$ is the proportional navigation law constant ${N}_{c}$ in the $l$th state. As a remark, the form of the equations description (Eq. (8)) for systems with random structure is first presented in this paper.
The measurements equation (Eq. (9)) represents a spatialtimevarying signal $S(x,\epsilon ,\mathrm{}t)$ observed by means of a coordinator with a narrow field of view which measures small deviations relative to the boresight:
where ${\mathbf{N}}^{\left(l\right)}\left(x,\mathrm{}t\right)$ is the spatialtimevarying noise of white Gaussian type with spectral densities ${C}_{0}^{\left(l\right)}$ in the case of stationary measurements noise, and ${C}_{0}^{\left(l\right)}\left(t\right)$, in the case of nonstationary noise depending on ${D}^{\left(l\right)}\left(t\right)$. The form of the signal which is a nonlinear function of the estimated Markov process $\epsilon \left(t\right)$ is determined by antenna pattern with the beam width of 3°7°. The changes of structure and the corresponding algorithms occur by the Markov process law with the two states $l=\overline{1,2}$ and the transition intensities ${\nu}_{12}\left(t\right)\text{,}$ and ${\nu}_{21}\left(t\right)\text{,}$ taking into account the values of estimation errors of a LOS angle. In a special case the system changes the mode of tracking when the estimation error exceeds the threshold value and the moments of the structure changes are known. There is a special case of the $\epsilon \left(t\right)$ coordinate jumping disturbance, when the system noise ${W}_{\epsilon}\left(t\right)$ is a function of a structure ${W}_{\epsilon}^{\left(l\right)}\left(t\right)$, particularly, when the disturbance is a hit.
The purpose of this part of the present scientific research is to create reliable nonlinear information processing algorithms for a target tracking system as a part of a guidancenavigation complex which provide suboptimal estimation of the range from an aircraft to a target $D\left(t\right)$, the line of sight angle $\epsilon \left(t\right)$, the bearing angle, and the guidance law parameters for different cases and conditions.
The adaptive filtering schemes are considered for 1) the proposed in this paper approach considering the navigation constant ${N}_{c}$ as a parameter to be identified, which changes its structure according to various conditions in contrast to all known previous solutions including [5] where the PN guidance law $\dot{\theta}\left(t\right)={N}_{c}\dot{\epsilon}\left(t\right)$ with the traditionally used navigation constant values ${N}_{c}$, interpreted as the constant values received by experimental approach was applied; 2) the proposed derived guidance law, in which the navigation coefficient ${N}_{c}$ is not a constant, (it is a function of the state estimates as an analytical expression); 3) the direct guidance law.
The estimate of the bearing angle was obtained for its state dynamics as a constant value and in the presence of noise.
The behavior of the system has been investigated with comparison for different first proposed conditions: 1) with adaptive switching from one to the other mode of operation using the information artificial aging techniques on the required parts of the trajectories (when the moment of the structure change depends on the estimation errors value); 2) with changes of the corresponding guidance law parameters, and 3) with switching from one guidance law to the other during tracking.
The adaptation mechanism of the proposed estimation algorithms represents an effective means of switching the field of view of a radar tracker or an optical (infrared image) tracker, and ensures the system operation, particularly, for a wide dynamic range of target maneuvers.
2. Conclusions
The proposed algorithms in comparison with the traditional estimation algorithms first show a relationship between the a posteriori probability density maximum criterion, the covariance matrix minimum criterion, and the crosscorrelation function maximum criterion.
The proposed expansion of mathematical description based on the nonlinear filtering theory of processes and the systems with random structure optimization provides adaptive features for the correlationextremum systems by generating the probabilistically weighted average of the state estimates. The correlationextremum methods combined with the theory of Markov processes and systems with random structure ensure the stochastic dynamic systems reliability and operation in different conditions and under environment influences.
These algorithms allow to solve the filtering and identification problems for complex stochastic dynamic systems equipped with image sensors (e.g., radar, optics, and etc.) using information of spatialtimevarying signals of different nature fields.
References

Baklitski V., Yuriev A. CorrelationExtremum Methods in Navigation. Radio and Communication, Moscow, 1982, (in Russian).

Chang C., Athans M. State estimation for discrete systems with switching parameters. IEEE Transactions on Electronic Systems, Vol. 14, Issue 3, 1978, p. 418425.

Kaminskas V., Sakalauskas E. Identification algorithm for linear systems with distributed parameters. Engineering Cybernetics, Vol. 3, 1987, p. 184189, (in Russian).

Kolosovskaya T. Spatialtimevarying signals processing algorithms in systems with random structure. Mechanical Engineering and Machine Reliability Problems, Vol. 5, 1995, p. 105112, (in Russian).

Shukla U., Mahapatra P. Generalized linear solution of proportional navigation. IEEE Transactions on Aerospace and Electronic Systems, Vol. 24, Issue 3, 1988, p. 231237.

Kolosovskaya T. The application of an extended Kalman filter to synthesis and analysis of estimation algorithms for magnetic navigation systems. 14th International Conference on Aviation and Cosmonautics, Moscow, 2015, p. 306309.
About this article
The author would like to thank the mentioned authors [13], [5] for their researches that stimulate further investigations.