METHOD FOR COMPENSATING FOR VISUAL-MEASUREMENT TIME LAG OF ELECTRO-OPTICAL TRACKING SYSTEM

Abstract

The present invention provides a generalized proportional integral observer-based method for compensating for visual-measurement time lag of an electro-optical tracking system. For visual-measurement time lag present in an electro-optical tracking system, an improved generalized proportional integral observer-based feedback control method is used to mitigate the impact of the measurement time lag on the system and suppress kinematic uncertainty of the system. The core of the method lies in that an observer is used to estimate a state, uncertainty, and a difference of the system at a previous moment, a state and uncertainty of the system at a current moment are then calculated by using these estimated values and a state-space model of the system, and a control input of the system is finally acquired according to the estimated values of the state and uncertainty of the system at the current moment. The method mitigates the adverse impact of visual-measurement time lag on the system and enhances the uncertainty suppression and the tracking precision of the system.

Description

BACKGROUND

Technical Field

The present invention relates to a method for compensating for visual-measurement time lag of an electro-optical tracking system, which belongs to the field of high-precision control technology for electro-optical tracking systems.

Related Art

An electro-optical tracking system is an apparatus that integrates optics, electricity, machinery, and control, and is widely applied to military and civilian fields. When the system is in a tracking mode, the visual axis of a mounted camera needs to remain aimed at a tracked target.

An electro-optical tracking platform is subject to control difficulties such as multi-source disturbance, parameter uncertainty, the maneuver of a tracked target, and visual-measurement time lag. In a two-axis inertially stabilized platform, the optical axis of a camera mounted on an inner frame is susceptible to interference due to the movement of a carrier, leading to considerable difficulty in electro-optical tracking control. In addition to the movement of the carrier, great kinematic uncertainty is also generated in the maneuver of the tracked target in space, especially during the tracking of a fast maneuvering target. Because the camera is used as a sensor in the system to obtain a miss distance of the tracked target, the processing of the image generates visual-measurement time lag, and the time lag reduces the stability and control precision of the system. To implement high-precision tracking control, it is necessary to consider both the visual-measurement time lag and the kinematic uncertainty in the system, and design a corresponding control method.

A Smith predictor is usually used in existing compensation for visual-measurement time lag of an electro-optical tracking system. However, models established in the system need to have very high precision in this method. It is very difficult to meet requirements in practical applications. In addition, only measurement time lag is processed in many methods, but kinematic uncertainty is not processed. As a result, the control precision of the system is reduced.

SUMMARY

To solve the foregoing problems, the present invention discloses a method for compensating for visual-measurement time lag of an electro-optical tracking system. In the method, a compound controller is designed by constructing an improved generalized proportional integral observer and combining the observer with a feedback linearization algorithm, to implement high-precision tracking of the electro-optical tracking system and suppression of kinematic uncertainty of the system. A specific technical solution of the compound controller is as follows:

step 1: installing a camera on an inner pitch frame of a two-axis inertially stabilized platform, connecting the camera to a host computer to form an electro-optical tracking system, shooting, by the camera, a tracked target in real time, extracting, by the host computer, a miss distance of the tracked target according to an image shot by the camera, and generating, by the two-axis inertially stabilized platform, a control input of the system according to the miss distance.

step 2: selecting the miss distance of the tracked target as a state of the system, and establishing a discrete-time state-space model considering visual-measurement time lag and kinematic uncertainty of the system:

$\{\begin{array}{c}X\left(k+1\right)=X\left(k\right)+BU\left(k\right)+\Delta \left(k\right)\\ Y\left(k\right)=X\left(k-d\right)\end{array},$

where X(k) represents the state of the system, U(k) represents the control input of the system, Δ(k) represents the kinematic uncertainty of the system, Y(k) represents a measurement output of the system, d represents the visual-measurement time lag, B represents a control input matrix, and k represents a k^th moment;

step 3: defining H(k)=Δ(k+1)−Δ(k) as a difference of the kinematic uncertainty of the system, and constructing, according to the discrete-time state-space model established in step 2, an improved generalized proportional integral observer:

$\{\begin{array}{c}{Z}_1\left(k+1\right)=Z_1\left(k\right)+BU\left(k-d\right)+Z_2\left(k\right)-L_1\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_2\left(k+1\right)=Z_2\left(k\right)+Z_3\left(k\right)-L_2\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_3\left(k+1\right)=Z_3\left(k\right)-L_3\left(Z_1\left(k\right)-Y\left(k\right)\right)\end{array},$

where Z₁(k), Z₂(k), and Z₃(k) represent states of the observer and are respectively estimated values of X(k−d), Δ(k−d), and Δ(k−d), and L₁, L₂, and L₃ represent gains of the observer.

step 4: calculating a predicted value {circumflex over (X)}(k) of a state of the system and a predicted value {circumflex over (Δ)}(k) of the kinematic uncertainty at a current moment according to the estimated values Z₁(k), Z₂(k), and Z₃(k) in step 3 and the discrete-time state-space model in step 2;

step 5: designing a compound controller U(k)=B⁻¹(−{circumflex over (Δ)}(k)−K{circumflex over (X)}(k)) according to the predicted values {circumflex over (X)}(k) and {circumflex over (Δ)}(k) obtained in step 4 and based on a feedback linearization algorithm, where K represents a parameter of the controller.

X(k), U(k), Δ(k), and Y(k) in the discrete-time state-space model in step 2 are respectively expressed as:

$X\left(k\right)=\left[\begin{array}{c}{x}_1\left(k\right)\\ x_2\left(k\right)\end{array}\right],U\left(k\right)=\left[\begin{array}{c}{u}_1\left(k\right)\\ u_2\left(k\right)\end{array}\right],\Delta \left(k\right)=\left[\begin{array}{c}{\Delta }_1\left(k\right)\\ {\Delta }_2\left(k\right)\end{array}\right],Y\left(k\right)=\left[\begin{array}{c}{y}_1\left(k\right)\\ y_2\left(k\right)\end{array}\right],$

where x₁(k) and x₂(k) respectively represent components of the miss distance of the tracked target in a yaw direction and a pitch direction, u₁(k) and u₂(k) respectively represent a pitch angular velocity and a yaw angular velocity of the camera, Δ₁(k) and Δ₂(k) respectively represent kinematic uncertainty of the system in the yaw direction and kinematic uncertainty of the system in the pitch direction, and y₁(k) and y₂(k) respectively represent measured values of the components of the miss distance of the tracked target in the yaw direction and the pitch direction.

The control input matrix B in the discrete-time state-space model in step 2 is:

$B=\left[\begin{array}{cc}0& -T\lambda m\\ T\lambda m& 0\end{array}\right],$

where T represents a sampling period of discretization of the system, λ represents the focal length of the camera, and m represents a ratio of one meter to the side length of a pixel.

The kinematic uncertainty H(k) and Z₁(k), Z₂(k), and Z₃(k) in the observer in step 3 are respectively expressed as:

$H\left(k\right)=\left[\begin{array}{c}{h}_1\left(k\right)\\ h_2\left(k\right)\end{array}\right],Z_1\left(k\right)=\left[\begin{array}{c}{z}_{11}\left(k\right)\\ z_{12}\left(k\right)\end{array}\right],Z_2\left(k\right)=\left[\begin{array}{c}{z}_{21}\left(k\right)\\ z_{22}\left(k\right)\end{array}\right],Z_3\left(k\right)=\left[\begin{array}{c}{z}_{31}\left(k\right)\\ z_{32}\left(k\right)\end{array}\right],$

where h₁(k) and h₂(k) respectively represent differences of uncertainty in a yaw direction and uncertainty in a pitch direction, z₁₁(k) and z₁₂(k) respectively represent estimated values of state variables x₁(k−d) and x₂(k−d), z₂₁(k) and z₂₂(k) respectively represent estimated values of uncertainty Δ₁(k−d) and uncertainty Δ₂(k−d), and z₃₁(k) and z₃₂(k) respectively represent estimated values of differences h₁(k−d) and h₂(k−d) of uncertainty.

L₁, L₂, and L₃ of the observer in step 3 are respectively:

$L_1=\left[\begin{array}{cc}{l}_{11}& 0\\ 0& l_{12}\end{array}\right],L_2=\left[\begin{array}{cc}{l}_{21}& 0\\ 0& l_{22}\end{array}\right],L_3=\left[\begin{array}{cc}{l}_{31}& 0\\ 0& l_{32}\end{array}\right],$

where I₁₁, I₁₂, I₂₁, I₂₂, I₃₁, and I₃₂ represent parameters of the gains of the observer, a characteristic root of a square matrix

$\left[\begin{array}{ccc}-L_1& I& O\\ -L_2& I& I\\ -L_3& O& I\end{array}\right]$

is within a unit circle, and I and O respectively represent a second-order identity matrix and a second-order zero matrix.

A specific process of calculating the predicted value {circumflex over (Δ)}(k) of the kinematic uncertainty and the predicted value {circumflex over (X)}(k) of the state at the current moment in step 4 is as follows:

$\{\begin{array}{c}\hat{\Delta }\left(k\right)=Z_2\left(k\right)+d·Z_3\left(k\right)\\ \hat{X}\left(k\right)=Z_1\left(k\right)+d·Z_2\left(k\right)+\frac{d\left(d-1\right)}{2}·Z_3\left(k\right)+\sum _{m=0}^{d-1}BU\left(k-1-m\right).\end{array}$

The compound controller in step 5 is specifically expressed as:

$K=\left[\begin{array}{cc}{k}_1& 0\\ 0& k_2\end{array}\right],$

where k₁ and k₂ represent to-be-designed parameters of the controller and satisfy: 0<k₁,k₂<1.

Beneficial effects of the present invention are as follows:

In the present invention, an improved generalized proportional integral observer is used to estimate a state, kinematic uncertainty, and a difference at a previous moment in an electro-optical tracking system, a state and uncertainty of the system at a current moment are predicted according to the obtained estimation information, and a compound controller is designed in combination with a feedback linearization algorithm. The method mitigates the adverse impact of visual-measurement time lag on tracking performance and enhances the uncertainty suppression and the tracking precision of the system.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1(a) is a structural diagram of an electro-optical tracking system according to the present invention;

FIG. 1(b) is a block diagram of a method for compensating for visual-measurement time lag of an electro-optical tracking system of the present invention;

FIG. 2(a) is a comparison diagram of the effect of controlling a miss distance between a DGPIO method that considers measurement time lag and a GPIO method that does not consider measurement time lag under the same visual-measurement time lag and kinematic uncertainty; and

FIG. 2(b) is a comparison diagram of control variables between a DGPIO method that considers measurement time lag and a GPIO method that does not consider measurement time lag under the same visual-measurement time lag and kinematic uncertainty.

DETAILED DESCRIPTION

The present invention is further described below with reference to the accompanying drawings and specific implementation. It should be understood that the following specific implementations are merely intended to describe the present invention rather than to limit the scope of the present invention.

A method for compensating for visual-measurement time lag of an electro-optical tracking system includes specific steps as follows:

Step 1: Install a camera on a pitch inner frame of a two-axis inertially stabilized platform, connect the camera to a host computer to form an electro-optical tracking system, shoot, by the camera, a tracked target in real time, extract, by the host computer, a miss distance of the tracked target according to an image shot by the camera, and generate, by the two-axis inertially stabilized platform, a control input of the system according to the miss distance.

As shown in FIG. 1(a), the two-axis inertially stabilized platform includes a yaw outer frame and a pitch inner frame. The camera is installed on the pitch inner frame. In this system, a desktop computer is used as the host computer. Real-time data acquired by the camera is connected to the host computer through a data cable. The host computer extracts the miss distance of the tracked target according to obtained picture information. The host computer sends the miss distance to a control chip (digital signal processor, DSP) through the data cable. The control chip calculates information about the control input of the system according to encoder information and gyroscope information. A torque motor generates corresponding torque.

Step 2: Select the miss distance of the tracked target as a state of the system, and establish a discrete-time state-space model considering visual-measurement time lag and kinematic uncertainty of the system, for example:

$\{\begin{array}{c}X\left(k+1\right)=X\left(k\right)+BU\left(k\right)+\Delta \left(k\right)\\ Y\left(k\right)=X\left(k-2\right)\end{array},$

where X(k) represents the state of the system, U(k) represents the control input of the system, Δ(k) represents the kinematic uncertainty of the system, Y(k) represents a measurement output of the system, B represents a control input matrix, and k represents a k^th moment.

In the model X(k), U(k), Δ(k), and Y(k) are respectively expressed as:

$X\left(k\right)=\left[\begin{array}{c}{x}_1\left(k\right)\\ x_2\left(k\right)\end{array}\right],U\left(k\right)=\left[\begin{array}{c}{u}_1\left(k\right)\\ u_2\left(k\right)\end{array}\right],\Delta \left(k\right)=\left[\begin{array}{c}{\Delta }_1\left(k\right)\\ {\Delta }_2\left(k\right)\end{array}\right],Y\left(k\right)=\left[\begin{array}{c}{y}_1\left(k\right)\\ y_2\left(k\right)\end{array}\right],$

where x₁(k) and x₂(k) respectively represent components of the miss distance of the tracked target in a yaw direction and a pitch direction, u₁(k) and u₂(k) respectively represent a pitch angular velocity and a yaw angular velocity of the camera, Δ₁(k) and Δ₂(k) respectively represent kinematic uncertainty of the system in the yaw direction and kinematic uncertainty of the system in the pitch direction, y₁(k) and y₂(k) respectively represent measured values of the components of the miss distance of the tracked target in the yaw direction and the pitch direction, and the control input matrix B and the kinematic uncertainty Δ(k) are respectively as follows:

$B=\left[\begin{array}{cc}0& -280\\ 280& 0\end{array}\right],\Delta \left(k\right)=\left[\begin{array}{c}{\Delta }_1\left(k\right)\\ {\Delta }_2\left(k\right)\end{array}\right]=\left[\begin{array}{c}15\sin \left(0.1k\right)+5\\ 15\sin \left(0.1k\right)+5\end{array}\right].$

Step 3: Define H(k)=Δ(k+1)−Δ(k) as a difference of the kinematic uncertainty of the system, and construct, according to the discrete-time state-space model established in step 2, an improved generalized proportional integral observer:

$\{\begin{array}{c}{Z}_1\left(k+1\right)=Z_1\left(k\right)+BU\left(k-2\right)+Z_2\left(k\right)-L_1\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_2\left(k+1\right)=Z_2\left(k\right)+Z_3\left(k\right)-L_2\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_3\left(k+1\right)=Z_3\left(k\right)-L_3\left(Z_1\left(k\right)-Y\left(k\right)\right)\end{array},$

where Z₁(k), Z₂(k), and Z₃(k) represent states of the observer and are respectively estimated values of X(k−d), Δ(k−d), and Δ(k−2), and L₁, L₂, and L₃ represent gains of the observer. Specific selection is as follows:

$L_1=\left[\begin{array}{cc}1/2& 0\\ 0& 1/2\end{array}\right],L_2=\left[\begin{array}{cc}3/4& 0\\ 0& 3/4\end{array}\right],L_3=\left[\begin{array}{cc}1/8& 0\\ 0& 1/8\end{array}\right].$

Step 4: Calculate a predicted value {circumflex over (X)}(k) of a state of the system and a predicted value {circumflex over (Δ)}(k) of the kinematic uncertainty at a current moment according to the estimated values Z₁(k), Z₂(k), and Z₃(k) in step 3 and the discrete-time state-space model in step 2. A specific process of calculating the predicted values is as follows:

$\hat{\Delta }\left(k\right)=Z_2\left(k\right)+2·Z_3\left(k\right)$ $\hat{X}\left(k\right)=Z_1\left(k\right)+2·Z_2\left(k\right)+Z_3\left(k\right)+\sum _{m=0}^1\mathrm{BU}\left(k-1-m\right).$

Step 5: Design a compound controller U(k)=B⁻¹(−{circumflex over (Δ)}(k)−K{circumflex over (X)}(k)) according to the predicted values {circumflex over (X)}(k) and Δ(k) obtained in step 4 and based on a feedback linearization algorithm, where K is specifically as follows:

$K=\left[\begin{array}{cc}\text{1/2}& 0\\ 0& \text{1/2}\end{array}\right].$

In the drawings, DGPIO represents the foregoing method, and REF represents a reference signal of the system.

The effectiveness and advantage of the present invention are described by using a comparative experiment below. Control variables generated in the two methods are basically kept the same for the rationality of comparison. A generalized proportional integral observer-based control method designed for an electro-optical tracking system without considering visual-measurement time lag as follows:

Step 1: Select the miss distance of the tracked target as a state of the system, and establish a discrete-time state-space model of the system that does not consider visual-measurement time lag but considers kinematic uncertainty as follows:

$\hspace{1em}\{\begin{array}{c}X\left(k+1\right)=X\left(k\right)+\mathrm{BU}\left(k\right)+\Delta \left(k\right)\\ Y\left(k\right)=X\left(k\right)\end{array}.$

Step 2: Construct, according to the discrete-time state-space model established in step 1, an improved generalized proportional integral observer:

$\{\begin{array}{c}{Z}_1\left(k+1\right)=Z_1\left(k\right)+\mathrm{BU}\left(k\right)+Z_2\left(k\right)-L_1\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_2\left(k+1\right)=Z_2\left(k\right)+Z_3\left(k\right)-L_2\left(Z_1\left(k\right)-Y\left(k\right)\right)\\ Z_3\left(k+1\right)=Z_3\left(k\right)-L_3\left(Z_1\left(k\right)-Y\left(k\right)\right)\end{array},$

where L₁, L₂, and L₃ represent gains of the observer. Specific selection is as follows:

$L_1=\left[\begin{array}{cc}1/2& 0\\ 0& 1/2\end{array}\right],L_2=\left[\begin{array}{cc}3/4& 0\\ 0& 3/4\end{array}\right],L_3=\left[\begin{array}{cc}1/8& 0\\ 0& 1/8\end{array}\right].$

Step 3: Design a compound controller U(k)=B⁻¹(−Z₂(k)−KZ₁(k)) according to the estimated values Z₁(k) and Z₂(k) in step 2 and based on a feedback linearization algorithm. K is specifically as follows:

$K=\left[\begin{array}{cc}1/2& 0\\ 0& 1/2\end{array}\right].$

In the drawings, GPIO represents the generalized proportional integral observer-based control method designed without considering visual-measurement time lag.

The technical means disclosed in the solutions of the present invention are not limited to the technical means disclosed in the foregoing technical means, but also include technical solutions including any combination of the foregoing technical features.

With the foregoing preferred embodiment according to the present invention as enlightenment, through the foregoing description, a related worker can make various changes and modifications without departing from the scope of the technical idea of the present invention. The technical scope of the present invention is not limited to the content of the description, and the technical scope of the present invention needs to be determined according to the scope of the claims.

Claims

1. A method for compensating for visual-measurement time lag of an electro-optical tracking system, comprising the following steps: step 1: installing a camera on a pitch inner frame of a two-axis inertially stabilized platform, connecting the camera to a host computer to form an electro-optical tracking system, shooting, by the camera, a tracked target in real time, extracting, by the host computer, a miss distance of the tracked target according to an image shot by the camera, and generating, by the two-axis inertially stabilized platform, a control input of the system according to the miss distance;step 2: selecting the miss distance of the tracked target as a state of the system, and establishing a discrete-time state-space model considering visual-measurement time lag and kinematic uncertainty of the system, for example:

2. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein X(k), U(k), Δ(k), and Y(k) in the discrete-time state-space model in step 2 are respectively expressed as:

3. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein the control input matrix B in the discrete-time state-space model in step 2 is:

4. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein the kinematic uncertainty H(k) and Z1(k), Z2(k), and Z3(k) in the observer in step 3 are respectively expressed as:

5. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein the gains L1, L2, and L3 of the observer in step 3 are respectively:

6. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein a specific process of calculating the predicted value {circumflex over (Δ)}(k) of the kinematic uncertainty and the predicted value {circumflex over (X)}(k) of the state at the current moment in step 4 is as follows:

7. The method for compensating for visual-measurement time lag of an electro-optical tracking system according to claim 1, wherein the compound controller in step 5 is specifically expressed as: