# Stability Precision Error Correction of Photoelectric Detection by Unmanned Aerial Vehicle

Authors: Huajie Hong; Keyan He; Zihao Gan; Guilin Jiang
Journal: Journal of Sensors
 (2021)

---

## Abstract

For getting clear images and overcoming shaking caused by various disturbances, real-time compensation of pointing errors will improve the overall stability performance of photoelectric detection by unmanned aerial vehicle. However, the compensation will be greatly deteriorated by error-causing sources, and the error correction process is of great importance. In this research, the problem of stability precision error correction is comprehensively studied. First, by modeling overall kinematics, error-causing sources, and error compensation, the error correction process is mathematically modeled and simulated. Then, by using simulation data regression, error correction models including the global function model and parametric model are established. The models are validated by carrying out both simulations and validation experiments. At last, the performances of the error correction models are compared and analyzed, which concerns the factors of parameter identification, model simplicity, and final improvement effect. Results show that the final stability precision can be greatly improved over 20%, and the parametric model outperforms the global function model comprehensively. It can be concluded that, either in simulation environment or real application scenarios, the obtained models and related analysis results are effective in improving the system stability performance.


---

## Body

## 1. Introduction

To obtain clear images and overcome disturbances, stability precision (hereinafter, the precision) is a key index for evaluating performances of photoelectric detection devices, especially for lightweight and small-sized carriers like unmanned aerial vehicle (UAV), for those are easily affected by various disturbances and errors [1, 2]. The typical application scenarios of this kind of device are shown in Figure 1.Figure 1
Typical application of UAV photoelectric detection.Researchers have been always focusing on the stability precision control of photoelectric devices. The application scenarios include carriers like unmanned vehicles, aerial carriers, and naval vessels. For instance, Zhuchong et al. [3] developed an inertial-stabilized platform for airborne remote sensing using magnetic bearings. Zhou et al. [4] presented the design of a mechatronic system for a two-axis inertial-stabilized platform in an unmanned helicopter- (UH-) based airborne power line inspection system. Dong et al. [5] proposed a dynamic model and control method for a two-axis initially stabilized platform, and there are plenty of works in other fields concerning the problem of stability precession improvement [6–8]. Therefore, in this research, stability precision is always treated as the key performance index to evaluate the whole system.Except for the casual environmental disturbances caused by carrier movements, winds, noises, etc., the over-all stability precision can also be severely affected by other error-causing sources (ECS) related to the quality of the device itself, like assembly errors, mass eccentricities, and sensor measurement accuracies. For example, Wu et al. [9] studied the device stabilization with total warping variation model. Yu et al. [10] studied video stabilization based on the modeling of motion imaging. Mao et al. [11] proposed the continuous second-order sliding mode control based on disturbance observer for LOS stabilization. Chen et al. [12] studied the disturbance observer-based control and related methods. Dasgupta et al. [13] designed the disturbance observer for nonlinear systems by using a Hirschorn inverse approach. Gao et al. [14] carried out an analysis on the influence of assembly error of off-axis three-mirror optical system on imaging line-of-sight. Liu et al. [15] carried out a machining error analysis of freeform surface off-axis three-mirror system based on optical performance evaluation.However, these studies commonly adopt the methods of off-line modeling and precorrection before the real application. Comparatively, the real-time error correction is considered in this research, which is beneficial for the real-time improvement during the real application process.The real-time compensation of the overall pointing errors will generally dissolve the disturbances coupled to the device base and improve the final precision. Specifically, to stabilize the line of sight (LOS) of the device, the kinematic pair of the azimuth frame and motor, and that of the pitch one, must be both used to compensate simultaneously. In this research, the two pairs are abbreviated as AZ and PI hereafter. However, the compensation performance will be greatly deteriorated by those ECS [16–18]. Therefore, during real applications, the extra error corrections aiming at dissolving the effects of ECS will be needed. The general correction process is shown in Figure 2.Figure 2
General error correction by using algorithms.The correction algorithm will be a parametric model used to predict the effects of ECS. It is used to calculate the pointing angle deviations between when those ECS exist and do not exist. Besides, when there are no ECS existing, the pointing angles are called theoretical pointing angles, which can be obtained by using analysis of general kinematics and compensation, and are often treated as the input parameters. In the real application scenarios, the real-time correction of those ECS effect will be realized by using extra servo control mechanism, including frames, motors, or more advanced structures, like multiaxis linkages, which are assembled on the outer side of AZ and PI. Those extra mechanisms will rotate according to the correction angles generated by using the algorithms above, so the deviations caused by ECS can be dissolved. Therefore, the design of the correction algorithm will be the key factor determining correction performance, which can be analyzed by comparing the final precision with the original precision without any correction. A better algorithm will predict the deviations more accurately and finally achieve a higher precision [19].For instance, Wang et al. [20] proposed the real-time FPGA-based Kalman filter for constant and nonconstant velocity periodic error correction. Xingshun et al. [21] studied the real-time phase correction based on FPGA in the beam position and phase measurement system. Refan et al. [22] proposed the real-time global poisoning stability improvement by utilizing support vector machine. Sen et al. [23] studied the real-time correction of periodic nonlinearity in detection for scanning beam interference. Yu et al. [24] carried out research on dynamic tracking and compensation method for imaging, and there are other similar researches related to the error correction modeling [25–27]. There are plenty of researches concerning real-time error corrections using methods like collective morphing, PID algorithms, and lateral control [28–30]. However, these studies typically focused on only one kind of error like that of beam axis, sensors, motors, or accelerators, without comprehensive modeling and analysis. Besides, they barely focused on the two-axis stabilization mechanism, which is commonly adopted by UAV photoelectric detection.In this research, based on our previous modeling works of the precision and various ECS, different correction algorithms are modeled and introduced into the overall simulation model, and the performances of those correction algorithms are compared. At last, the chosen algorism is validated by carrying out both simulation and experiment using data of a real product. Therefore, the real-time error correction of the UVA stability precision is comprehensively studied. Furthermore, it can be concluded that this research fills the technical deficiency of using a comprehensive method to compensate for the stability error of UAV photoelectric detection and provides a pioneering theoretical guidance for the related engineering applications in the future.The organization of this paper is as follows: the background and literature review are presented in Section1; the modeling and simulations of ECS are presented in Section 2; the error correction modeling is presented in Section 3; the simulation and validation experiments are presented in Section 4; the final conclusions and future work are presented in Section 5.


## 2. Modeling and Simulations of UAV, Stability Precision, and Error-Causing Sources

Based on our previous research of error modeling [31], the necessary predefinitions are shown as follows, and all modeling components are integrated together for the later analysis. Furthermore, because the device base is considered as being fixed to the UAV body, the dynamic modeling and related analysis of UAV can be simplified as that of the device base, which is presented as follows:
(i)
The final precision is denoted as an index ofRMS, and it is defined as the root mean square of the angular deviations Δθsi. The values of Δθsi are the differences between the ideal pointing positions and the real ones, and they are also defined as the pointing errors. The specific method to calculate Δθsi and RMS can be seen in our previous work [31] by using a series of statistic indexes(ii)
The angular velocity of the device base related to the ground is denoted asω→ib=ωibx,ωiby,ωibzT. That of AZ is denoted as ω→ia=ωiax,ωiay,ωiazT; the angle and the velocity of AZ compensation, which are driven by AZ motor, are denoted as θa and θ̇a. Those parameters of PI are denoted as ω→ip=ωipx,ωipy,ωipzT, θp, and θ̇p, respectively(iii)
For the ECS modeling, in this research, the commonly-seen ECS of perpendicularity deviations includingδ→a=δax,δay,0T and δ→p=δpx,0,δpzT, which can be measured by using high precision instruments [32], run-outs [33], including γ→a=γax,γay,0T and γ→p=γpx,0,γpzT, and assembly deviations of gyroscope locations [34] including ε→a=εax,εay,0T and ε→p=εpx,0,εpzT, are considered(iv)
For the kinematics modeling, as mentioned, AZ and PI must be both used to compensate simultaneously. According to kinematics analysis of homogeneous coordinate transformation put forward by [35], the rotation transformation matrix from the device base to AZ and that from AZ to PI is denoted as Rba and Rap. Those of perpendicularity deviations are denoted as Ra′ and Rb′, those of run-outs are denoted as ΔRa and ΔRp, and those of gyroscope location deviations are denoted as Rag and RpgMore details of those ECS and related kinematics analysis can be seen in our previous work [31]. By integrating all the above modules, the overall stability precision is modeled in Simulink, which is shown in Figure 3. The original system base angular velocity ω→ib=ωibx,ωiby,ωibzT is set as the input. All of the modules of the kinematics, components, and ECS are embedded in the middle, and the details of those modules can also be seen in our previous work [31]. The final index of stability precision RMS is set as the output.Figure 3
Overall simulation modeling of the stability precision.


## 3. Error Correction Modeling

### 3.1. Necessary Predefinitions and Preparations before the Modeling

As mentioned, the error correction process is used to dissolve the effects caused by ECS. The correction algorithm model is used to predict the pointing angle deviations between when those ECS exist and do not exist. When there are no ECS existing, those are called the theoretical pointing angles. For the case of the two-axis-two-frame structure of AZ and PI, the theoretical pointing angles are denoted asαi,βi, where αi denotes the AZ pointing angle, βi denotes the PI one, i=1,2,3,¡­,nm, and nm is the total number of the samples used for modeling. αi,βi can be obtained by using methods of general kinematics and compensation and are often treated as the input parameters of the correction model.Respectively, when there are ECS existing, those are called the real pointing angles, which are denoted asα^i,β^i. The deviations between the theoretical and the real are denoted as δαi,δβi. In that case, the general form of the error correction model can be expressed as (1), where fαi,βi and gαi,βi are the model functions, and εi and ηi are the prediction deviations of the model.
(1)δαi=α^i−αi=fαi,βi+εi,δβi=β^i−βi=gαi,βi+ηi.Therefore, the modeling process is to find the most appropriate forms offαi,βi and gαi,βi to obtain the minimum accumulation values of εi and ηi, which is expressed by (2):
(2)min∑i=1nmεi2=min∑i=1nmfαi,βi−δαi2,min∑i=1nmηi2=min∑i=1nmgαi,βi−δβi2.Thenm groups of sampling values of δαi,δβi and αi,βi are the data basis for the modeling process. Based on simulation models presented in Section 2 and by carrying out the simulation with all those ECS set to zeros, the signal values of Δθszi and Δθsyi are taken as the sample values of αi,βi by definition. Similarly, the differences of those signal values, between when those ECS are set to certain values and zeros, are taken as the sample values of δαi,δβi. Additionally, the ECS are all set to be 0.5 degrees considering the actual manufacturing and assembly accuracy. The sampling data process is shown in Figure 4.Figure 4
Data-acquiring of the samples used for the modeling.About more details of the above data-acquiring process, the UAV movements, and other various disturbances coupled to the device base, namely,ω→ib=ωibx,ωiby,ωibzT, are simulated by using back-ward deduction of a power spectrum. In addition, based on real measurement data of UAV movements in outdoor environments, the power spectrum is obtained by using the method of correction function estimation [36]. The backward deduction means that the amplitude of the simulated signal is obtained by inverse Fourier transform of the power spectrum data, the phase is randomly generated, and the three directions of ω→ib are all set equal. Those processes and obtained data are shown in Figure 5, and the relative MATLAB codes are also attached as supplementary files.Figure 5
Simulation of the UAV moves and other disturbances.Based on the above background, 13600 groups ofδαi,δβi and αi,βi are obtained, namely, nm is set as 13600. Those obtained data groups are also attached as supplementary files, which can be used as the basis for the error correction modeling.


### 3.2. Global Function Model

Global function model (GFM) [37] takes the form of a spherical function of Fourier series expansion hα,β as (3), where 0≤α≤2π and 0≤β≤2π denote the rotary angles of a time point, Amα and Bmα denote the Fourier series coefficients, and m denotes the Fourier series number of an expansion term.
(3)hα,β=∑m=0∞Amαcmβ+Bmαsmβ.Then by using the spherical function as the basis,Amα and Bmα are expanded in 0,2π as (4), where Plmcα is the l-order associated Legendre polynomial [38], and Alm and Blm are the coefficients.
(4)Amα=∑l=m∞AlmPlmcα,Bmα=∑l=m∞BlmPlmcα.By introducing (4) into (3), we get (5), which can be further expanded as a form of polynomial shown as (6).
(5)h=∑m=0∞∑l=m∞Almcmβ+BlmsmβPlmcα,(6)hα,β=a0+a1sβ+a2cαcβ+a3sαcβ+a4s2β+a5cαsβcβ+a6sαsβcβ+a7c2αc2β+a8s2αc2β+a9s3β+a10cαs2βcβ+a11sαs2βcβ+a12c2αc2βsβ.Therefore, according to the definition of error correction models shown in (1), the GFM model can be expressed as (7) by introducing (6) into (1) for the AZ and PI, respectively:
(7)δαicβi=fαi,βi=a0+a1sβi+⋯+a12c2αic2βisβi+εi,δβi=gαi,βi=b0+b1sβi+⋯+b12c2αic2βisβi+ηi.Equation (7) is expressed in form of a matrix as in (8):
(8)δαδβ=AαAβx+εη.Those symbols in (8) are explained in (9). By carrying out regression analysis using the 13600 data groups mentioned in Subsection 3.1, we get the result of x=−5.495×103,1.000×10−5,1.000×10−6,1.000×10−6,−1.211×102,1.000×10−5,1.906×104,5.495×103,−4.239×10−1,3.617×10−1,1.098×100,1.000×10−6,1.000×10−5,−1.476×100,2.913×104,3.045×103,−1.052×10−1,1.000×10−6,1.000×10−4,1.000×10−5,1.000×10−6, which means that the GFM model is obtained. Apparently, GFM takes the form originated from general spherical function, and it contains many parameters without specific physical meaning. That is not beneficial for deep understanding of ECS.
(9)δα=δα1,δα2,⋯,δαnn×1T,δβ=δβ1,δβ2,⋯,δβnn×1T,ε=ε1,ε2,⋯,εnn×1T,η=η1,η2,⋯,ηnn×1T,x=a0,a1,⋯,a12,b0,b1,⋯,b1226×1T,Aa=1/cβ1sβ1/cβ1⋯c2α1cβ1sβ101/cβ2sβ2/cβ2⋯c2α2cβ2sβ20⋮⋮⋮⋮⋮1/cβnsβn/cβn⋯c2αncβnsβn0n×26,Aβ=01sβ1⋯c2α1c2β1sβ101sβ2⋯c2α2c2β2sβ2⋮⋮⋮⋮⋮01sβn⋯c2αnc2βnsβnn×26,n=nm=13600.


### 3.3. Parametric Model

Parametric model (PM) is established by using general movement theory and previous analysis results of ECS in Section2. The target position in the base coordinate Cb can be calculated by (10), where Ppgi is one of the sample values of the target position in the Cpg coordinate at the end of the coordinate series, and Pbi=xbi,ybi,zbi is that in the Cb coordinate at the begging of the series. By doing that, the real-time target position in the absolute coordinate can be obtained just by using the gyroscope signals.
(10)Pbi=Ra′ΔRaRbaRp′ΔRpRapRagRpgPpgi.We introduce the contents in Section2 into (10) and then simplify the result by adhering to rules as follows: first, any of the elements Δ belongs to ECS are treated as small angles, so sΔ≈0, and cΔ≈1; second, all of the high-order elements are neglected. Simultaneously, the transfer relation between the rectangular coordinate and spherical coordinate is shown as (11), where d is the distance away from the target.
(11)tanαi+δαi=−ybixbitanβi+δβi=zbid⇒sαi+cαiδαicαi−sαiδαi=−ybixbi,sβi+cβiδβi=zbid.Finally, the PM model can be expressed as (12):
(12)δαi≈εay+εax+δpx−δpz+γaytanβi+γpxsecβi+δaxsαitanβi−δaycαitanβi+εi,δβi≈εpz+εpx+γax+δaxcαi+δaysαi+ηi.Equation (12) is expressed in form of a matrix as (13):
(13)δαδβ=AαAβx+εη.Those symbols in Eq. (13) are explained in Eq. (14), where x1=εay+εax+δpx,x2=−δpz+γay,x3=γpx,x4=δax,x5=−δay,y1=εpz+εpx+γax,y2=δax and y3=δay.
(14)δα=δα1,δα2,⋯,δαnn×1T,δβ=δβ1,δβ2,⋯,δβnn×1T,ε=ε1,ε2,⋯,εnn×1T,η=η1,η2,⋯,ηnn×1T,x=x1,x2,x3,x4,x5,y1,y2,y38×1T,Aα=1tanβ1secβ1sα1tβ1cα1tβ101tanβ2secβ2sα2tβ2cα2tβ20⋮⋮⋮⋮⋮⋮1tanβnsecβnsαntβncαntβn0n×8,Aβ=01cα1sα1⋮⋮⋮⋮01cα2sα201cαnsαnn×8,n=nm=13600.By carrying out regression analysis using the 13600 data groups mentioned in Subsection3.1, we get the result of x=−5.963×103,3.571×10−1,5.963×103,−9.279×102,1.000×10−5,7.977×103,−7.977×103,9.472×10−1, which means that the PM model is obtained. Compared to GFM, PM has few parameters. Those parameters can be explained by referring to the meaning of the predefined ECS, which is beneficial for possible future improvement.


## 3.1. Necessary Predefinitions and Preparations before the Modeling

As mentioned, the error correction process is used to dissolve the effects caused by ECS. The correction algorithm model is used to predict the pointing angle deviations between when those ECS exist and do not exist. When there are no ECS existing, those are called the theoretical pointing angles. For the case of the two-axis-two-frame structure of AZ and PI, the theoretical pointing angles are denoted asαi,βi, where αi denotes the AZ pointing angle, βi denotes the PI one, i=1,2,3,¡­,nm, and nm is the total number of the samples used for modeling. αi,βi can be obtained by using methods of general kinematics and compensation and are often treated as the input parameters of the correction model.Respectively, when there are ECS existing, those are called the real pointing angles, which are denoted asα^i,β^i. The deviations between the theoretical and the real are denoted as δαi,δβi. In that case, the general form of the error correction model can be expressed as (1), where fαi,βi and gαi,βi are the model functions, and εi and ηi are the prediction deviations of the model.
(1)δαi=α^i−αi=fαi,βi+εi,δβi=β^i−βi=gαi,βi+ηi.Therefore, the modeling process is to find the most appropriate forms offαi,βi and gαi,βi to obtain the minimum accumulation values of εi and ηi, which is expressed by (2):
(2)min∑i=1nmεi2=min∑i=1nmfαi,βi−δαi2,min∑i=1nmηi2=min∑i=1nmgαi,βi−δβi2.Thenm groups of sampling values of δαi,δβi and αi,βi are the data basis for the modeling process. Based on simulation models presented in Section 2 and by carrying out the simulation with all those ECS set to zeros, the signal values of Δθszi and Δθsyi are taken as the sample values of αi,βi by definition. Similarly, the differences of those signal values, between when those ECS are set to certain values and zeros, are taken as the sample values of δαi,δβi. Additionally, the ECS are all set to be 0.5 degrees considering the actual manufacturing and assembly accuracy. The sampling data process is shown in Figure 4.Figure 4
Data-acquiring of the samples used for the modeling.About more details of the above data-acquiring process, the UAV movements, and other various disturbances coupled to the device base, namely,ω→ib=ωibx,ωiby,ωibzT, are simulated by using back-ward deduction of a power spectrum. In addition, based on real measurement data of UAV movements in outdoor environments, the power spectrum is obtained by using the method of correction function estimation [36]. The backward deduction means that the amplitude of the simulated signal is obtained by inverse Fourier transform of the power spectrum data, the phase is randomly generated, and the three directions of ω→ib are all set equal. Those processes and obtained data are shown in Figure 5, and the relative MATLAB codes are also attached as supplementary files.Figure 5
Simulation of the UAV moves and other disturbances.Based on the above background, 13600 groups ofδαi,δβi and αi,βi are obtained, namely, nm is set as 13600. Those obtained data groups are also attached as supplementary files, which can be used as the basis for the error correction modeling.


## 3.2. Global Function Model

Global function model (GFM) [37] takes the form of a spherical function of Fourier series expansion hα,β as (3), where 0≤α≤2π and 0≤β≤2π denote the rotary angles of a time point, Amα and Bmα denote the Fourier series coefficients, and m denotes the Fourier series number of an expansion term.
(3)hα,β=∑m=0∞Amαcmβ+Bmαsmβ.Then by using the spherical function as the basis,Amα and Bmα are expanded in 0,2π as (4), where Plmcα is the l-order associated Legendre polynomial [38], and Alm and Blm are the coefficients.
(4)Amα=∑l=m∞AlmPlmcα,Bmα=∑l=m∞BlmPlmcα.By introducing (4) into (3), we get (5), which can be further expanded as a form of polynomial shown as (6).
(5)h=∑m=0∞∑l=m∞Almcmβ+BlmsmβPlmcα,(6)hα,β=a0+a1sβ+a2cαcβ+a3sαcβ+a4s2β+a5cαsβcβ+a6sαsβcβ+a7c2αc2β+a8s2αc2β+a9s3β+a10cαs2βcβ+a11sαs2βcβ+a12c2αc2βsβ.Therefore, according to the definition of error correction models shown in (1), the GFM model can be expressed as (7) by introducing (6) into (1) for the AZ and PI, respectively:
(7)δαicβi=fαi,βi=a0+a1sβi+⋯+a12c2αic2βisβi+εi,δβi=gαi,βi=b0+b1sβi+⋯+b12c2αic2βisβi+ηi.Equation (7) is expressed in form of a matrix as in (8):
(8)δαδβ=AαAβx+εη.Those symbols in (8) are explained in (9). By carrying out regression analysis using the 13600 data groups mentioned in Subsection 3.1, we get the result of x=−5.495×103,1.000×10−5,1.000×10−6,1.000×10−6,−1.211×102,1.000×10−5,1.906×104,5.495×103,−4.239×10−1,3.617×10−1,1.098×100,1.000×10−6,1.000×10−5,−1.476×100,2.913×104,3.045×103,−1.052×10−1,1.000×10−6,1.000×10−4,1.000×10−5,1.000×10−6, which means that the GFM model is obtained. Apparently, GFM takes the form originated from general spherical function, and it contains many parameters without specific physical meaning. That is not beneficial for deep understanding of ECS.
(9)δα=δα1,δα2,⋯,δαnn×1T,δβ=δβ1,δβ2,⋯,δβnn×1T,ε=ε1,ε2,⋯,εnn×1T,η=η1,η2,⋯,ηnn×1T,x=a0,a1,⋯,a12,b0,b1,⋯,b1226×1T,Aa=1/cβ1sβ1/cβ1⋯c2α1cβ1sβ101/cβ2sβ2/cβ2⋯c2α2cβ2sβ20⋮⋮⋮⋮⋮1/cβnsβn/cβn⋯c2αncβnsβn0n×26,Aβ=01sβ1⋯c2α1c2β1sβ101sβ2⋯c2α2c2β2sβ2⋮⋮⋮⋮⋮01sβn⋯c2αnc2βnsβnn×26,n=nm=13600.


## 3.3. Parametric Model

Parametric model (PM) is established by using general movement theory and previous analysis results of ECS in Section2. The target position in the base coordinate Cb can be calculated by (10), where Ppgi is one of the sample values of the target position in the Cpg coordinate at the end of the coordinate series, and Pbi=xbi,ybi,zbi is that in the Cb coordinate at the begging of the series. By doing that, the real-time target position in the absolute coordinate can be obtained just by using the gyroscope signals.
(10)Pbi=Ra′ΔRaRbaRp′ΔRpRapRagRpgPpgi.We introduce the contents in Section2 into (10) and then simplify the result by adhering to rules as follows: first, any of the elements Δ belongs to ECS are treated as small angles, so sΔ≈0, and cΔ≈1; second, all of the high-order elements are neglected. Simultaneously, the transfer relation between the rectangular coordinate and spherical coordinate is shown as (11), where d is the distance away from the target.
(11)tanαi+δαi=−ybixbitanβi+δβi=zbid⇒sαi+cαiδαicαi−sαiδαi=−ybixbi,sβi+cβiδβi=zbid.Finally, the PM model can be expressed as (12):
(12)δαi≈εay+εax+δpx−δpz+γaytanβi+γpxsecβi+δaxsαitanβi−δaycαitanβi+εi,δβi≈εpz+εpx+γax+δaxcαi+δaysαi+ηi.Equation (12) is expressed in form of a matrix as (13):
(13)δαδβ=AαAβx+εη.Those symbols in Eq. (13) are explained in Eq. (14), where x1=εay+εax+δpx,x2=−δpz+γay,x3=γpx,x4=δax,x5=−δay,y1=εpz+εpx+γax,y2=δax and y3=δay.
(14)δα=δα1,δα2,⋯,δαnn×1T,δβ=δβ1,δβ2,⋯,δβnn×1T,ε=ε1,ε2,⋯,εnn×1T,η=η1,η2,⋯,ηnn×1T,x=x1,x2,x3,x4,x5,y1,y2,y38×1T,Aα=1tanβ1secβ1sα1tβ1cα1tβ101tanβ2secβ2sα2tβ2cα2tβ20⋮⋮⋮⋮⋮⋮1tanβnsecβnsαntβncαntβn0n×8,Aβ=01cα1sα1⋮⋮⋮⋮01cα2sα201cαnsαnn×8,n=nm=13600.By carrying out regression analysis using the 13600 data groups mentioned in Subsection3.1, we get the result of x=−5.963×103,3.571×10−1,5.963×103,−9.279×102,1.000×10−5,7.977×103,−7.977×103,9.472×10−1, which means that the PM model is obtained. Compared to GFM, PM has few parameters. Those parameters can be explained by referring to the meaning of the predefined ECS, which is beneficial for possible future improvement.


## 4. Simulation and Validation Experiment of Error Correction

Based on the modeling results in Section3, both simulation and validation experiments are carried out to evaluate the performances of the modeling methods.


### 4.1. Simulations of Error Correction

The simulation process is shown in Figure6, in which the device base movements are simulated by signals generated in Subsection 3.1. The overall integrated precision model in Section 2 is used to generate real-time real pointing angles and theoretical pointing angles. Different correction models, including GFM and PM in Subsections 3.2 and 3.3, and the most commonly used models of ADRC, DAC, and CHADC in [4, 12, 20], are all used to calculate the correction angles. The final stability precisions after correction are also compared to evaluate the performances of the correction models.Figure 6
The simulation process of error correction.The modeling and simulations are all carried out in MATLAB R2019b Simulink, and the MATLAB codes of the simulation works are available online at Supplementary Materials.The necessary parameters of the simulation are listed in Table1, and other parameters for the calculations are set equal to those listed in Table 2 shown later, which is originated from a real pod product. The simulation results are listed in Table 3. The deviation angles Δθsi along the simulation process are presented in Figure 7, which is used to calculate the final RMS.Table 1
Preset parameters of the simulation.
Software versionMATLAB R2019b SimulinkDuration time0-100 sSolver typeVariable-step (ode45)Tolerance of relative and absoluteAutomaticOutput optionRefine outputTable 2
Necessary parameters of the product for the experiment.
PARMDescriptionsValuesUnitsJmaRotary inertia of the azimuth motor0.02kg·m2JLaRotary inertia of the azimuth load0.04kg·m2JmpRotary inertia of the pitch motor0.01kg·m2JLpRotary inertia of the pitch load0.02kg·m2RmTotal resistor of the two armature circuits8.60ΩLmTotal inductance of the armature circuits0.01HDLDamping coefficient of the two motors10/KLElastic coefficient of the two motors100000/CeEMF coefficient of the two motors0.33/CmMotor torque coefficient of the two motors0.33/KpAmplification coefficient of driving circuit20/RsResistor of the driving circuit0.4ΩKgAmplification coefficient of the gyroscopes15/π/ωgNatural frequency of the gyroscopes200πrad/sλgDamping coefficient of the gyroscopes0.7/Table 3
The simulation results of error correction.
Correction model usedStability precisionRMS (μrad)Improvement percentage related to the original (%)The originalRMS without correction169.3880GFM165.7032.175PM135.20720.179ADRC156.2327.766DAC163.4323.516CHADC149.52311.728Figure 7
The deviationsΔθsi along the simulations.From Table3 and Figure 7, it is obvious that PM outperformed GFM and the others in improving the overall stability precision, and the improvement percentage can reach up to over 20%.


### 4.2. Validation Experiment of Error Correction

For further validation of the correction models, validation experiments using a real product assembled on an experiment platform is carried out, as shown in Figure8.Figure 8
The studied product on the experiment platform. (1) High precision servo turntable. (2) Photoelectric system product. (3) PSD position measurement system. (4) Two-axis position platform.The details of the experiment platform are listed as follows:(i)
The high precision servo turntables is used to simulate the UAV movements and other disturbances coupled to the product base, namely,ω→ib=ωibx,ωiby,ωibzT. Those movements are also generated by using signals in Subsection 3.1, and based on which, the theoretical pointing angles αi,βi are calculated in the background software to support the real-time error correction process using GFM or PM. Although all the simulation and experiments are carried out indoors, the movements are simulated according to the real outdoor environments(ii)
The photoelectric system product is a UAV detection pod(iii)
The PSD position measurement system is used to measure the real pointing anglesα^i,β^i and further calculate the final stability precision RMS(iv)
The two-axis position platform is used to adjust the position of the PSD system for better measurementBy referring to real measurement data of the pod product, the necessary parameters, used for the theoretical pointing angles calculation, are preset before the experiments. Those parameters are listed in Table2.GFM, PM, and other commonly used methods are both used to do the correction. The experiment results are listed in Table4. From Table 4, PM also outperformed GFM and the others, and the improvement percentage can reach up to over 15%.Table 4
The validation experiment results of error correction.
Correction model usedStability precisionRMS (μrad)Improvement percentage related to the original (%)The originalRMS without correction357.4530GFM343.2143.983PM302.38215.407ADRC347.4532.798DAC332.6756.932CHADC329.7657.746


## 4.1. Simulations of Error Correction

The simulation process is shown in Figure6, in which the device base movements are simulated by signals generated in Subsection 3.1. The overall integrated precision model in Section 2 is used to generate real-time real pointing angles and theoretical pointing angles. Different correction models, including GFM and PM in Subsections 3.2 and 3.3, and the most commonly used models of ADRC, DAC, and CHADC in [4, 12, 20], are all used to calculate the correction angles. The final stability precisions after correction are also compared to evaluate the performances of the correction models.Figure 6
The simulation process of error correction.The modeling and simulations are all carried out in MATLAB R2019b Simulink, and the MATLAB codes of the simulation works are available online at Supplementary Materials.The necessary parameters of the simulation are listed in Table1, and other parameters for the calculations are set equal to those listed in Table 2 shown later, which is originated from a real pod product. The simulation results are listed in Table 3. The deviation angles Δθsi along the simulation process are presented in Figure 7, which is used to calculate the final RMS.Table 1
Preset parameters of the simulation.
Software versionMATLAB R2019b SimulinkDuration time0-100 sSolver typeVariable-step (ode45)Tolerance of relative and absoluteAutomaticOutput optionRefine outputTable 2
Necessary parameters of the product for the experiment.
PARMDescriptionsValuesUnitsJmaRotary inertia of the azimuth motor0.02kg·m2JLaRotary inertia of the azimuth load0.04kg·m2JmpRotary inertia of the pitch motor0.01kg·m2JLpRotary inertia of the pitch load0.02kg·m2RmTotal resistor of the two armature circuits8.60ΩLmTotal inductance of the armature circuits0.01HDLDamping coefficient of the two motors10/KLElastic coefficient of the two motors100000/CeEMF coefficient of the two motors0.33/CmMotor torque coefficient of the two motors0.33/KpAmplification coefficient of driving circuit20/RsResistor of the driving circuit0.4ΩKgAmplification coefficient of the gyroscopes15/π/ωgNatural frequency of the gyroscopes200πrad/sλgDamping coefficient of the gyroscopes0.7/Table 3
The simulation results of error correction.
Correction model usedStability precisionRMS (μrad)Improvement percentage related to the original (%)The originalRMS without correction169.3880GFM165.7032.175PM135.20720.179ADRC156.2327.766DAC163.4323.516CHADC149.52311.728Figure 7
The deviationsΔθsi along the simulations.From Table3 and Figure 7, it is obvious that PM outperformed GFM and the others in improving the overall stability precision, and the improvement percentage can reach up to over 20%.


## 4.2. Validation Experiment of Error Correction

For further validation of the correction models, validation experiments using a real product assembled on an experiment platform is carried out, as shown in Figure8.Figure 8
The studied product on the experiment platform. (1) High precision servo turntable. (2) Photoelectric system product. (3) PSD position measurement system. (4) Two-axis position platform.The details of the experiment platform are listed as follows:(i)
The high precision servo turntables is used to simulate the UAV movements and other disturbances coupled to the product base, namely,ω→ib=ωibx,ωiby,ωibzT. Those movements are also generated by using signals in Subsection 3.1, and based on which, the theoretical pointing angles αi,βi are calculated in the background software to support the real-time error correction process using GFM or PM. Although all the simulation and experiments are carried out indoors, the movements are simulated according to the real outdoor environments(ii)
The photoelectric system product is a UAV detection pod(iii)
The PSD position measurement system is used to measure the real pointing anglesα^i,β^i and further calculate the final stability precision RMS(iv)
The two-axis position platform is used to adjust the position of the PSD system for better measurementBy referring to real measurement data of the pod product, the necessary parameters, used for the theoretical pointing angles calculation, are preset before the experiments. Those parameters are listed in Table2.GFM, PM, and other commonly used methods are both used to do the correction. The experiment results are listed in Table4. From Table 4, PM also outperformed GFM and the others, and the improvement percentage can reach up to over 15%.Table 4
The validation experiment results of error correction.
Correction model usedStability precisionRMS (μrad)Improvement percentage related to the original (%)The originalRMS without correction357.4530GFM343.2143.983PM302.38215.407ADRC347.4532.798DAC332.6756.932CHADC329.7657.746


## 5. Conclusions and Future Work

The major contributions of this research are as follows. Additionally, the nomenclature is presented in the end.(1)
Stability precision error correction of UAV photoelectric detection is comprehensively defined and analyzed. The analysis results and related data can be valuable guidance for real application(2)
Error correction models including global function model GFM and parametric model PM are established. The models are validated by carrying out both simulations and validation experiments and can be used in real applications(3)
The performances of the studied error correction models are compared and analyzed. PM outperformed GFM and other methods comprehensively, and the improvement percentage can reach up to over 20% for the simulation, and over 15% for the experimentThe limits of this research are listed as follows: (1) So far, only linear models are studied. Other error correction models with nonlinear expressions or fuzzy algorithms should be considered. (2) Other error sources and situations except for assembly errors to further verify the proposed analysis method should be considered, like mass eccentricities, deviations gyro measurements, and various kinds of noises. (3) The limitations, constraints used in this research, and the possibility of uncertainty are not concerned, like the changes of error distributions, target locations, and UAV linear disturbances.As the next step, we will investigate other error correction methods like semiparametric model, BP neural net, and machining learning. The same is true for other error-causing sources and situations, like frictions, structural eccentricity, and inertia coupling. The definition of limitations and constraints should be well-adjusted according to the scenario. For example, the linear disturbances of UAV and target movements may have a great influence on the results. Similar control mechanism of more complex system with more dimensions or new materials will be studied, like two-axis-four-frame and three-axis-four-frame used on larger carriers.



---

*Source: 5564448-2021-04-05.xml*