#### Abstract

The vehicle mounted and unbalanced barrel pitching system is a typical kind of high precision electrohydraulic servo system. The key to system control is not only to overcome the influence of external interference such as unbalanced torque and road fluctuation, but also to suppress the influence of internal parameters’ change induced by oil contamination. Firstly, the overall model of barrel pitching system is established and its sensitivity analysis is conducted, and then the correction link and the pressure difference feedback loop are respectively designed to suppress the influence of unbalanced torque and road fluctuation. Secondly, the influence of oil contamination is categorized and the core decay parameters are determined, and then the fractional order fuzzy PID controller is designed. Finally, the simulation model of MATLAB/Simulink is built and the system simulation is carried out. The results show when the given angle is fixed on 25°, the setting time is reduced to 0.2 s, and the maximum steady state error is only 0.000132 rad. Thus, the presented method is effective for barrel pitching system’s control, and the service life of the equipment can be prolonged under the same required accuracy.

#### 1. Introduction

In order to improve the shooting accuracy, the tank is usually equipped with an artillery stabilizer that can effectively isolate the effect of vehicle motion, as shown in Figure 1. The artillery stabilizer, also known as gun control system, can make the barrel unaffected by the pitching and steering of vehicle and then ensure the tank shooting during marching. It can also control the position of the barrel in the vertical direction and the horizontal direction so as to precisely aim [1], in which the vertical pitching motion attracts more attention due to the greater effect of road fluctuation.

The electrohydraulic servo system is more adopted for stability and turnover of the barrel, which has the advantages of large power to weight ratio, fast response, and high control precision. However, the outstanding weakness of this system is its low antipollution capability. According to authoritative statistics, 70%–80% of the hydraulic faults are caused by oil contamination. Also the electrohydraulic servo system has many uncertainties such as external interference, which seriously restricts the reaction speed and control precision [2, 3].

Specific to the barrel pitching system, external interference such as unbalanced torque and road fluctuation can seriously influence the control accuracy. In addition, the “blood of hydraulic system,” that is, the working oil, will be polluted with the prolonged service life under severe working conditions. Oil contamination will induce oil temperature rise and wear aggravate, which causes leakage increase, pressure instability, and action abnormality, and then result in frequent system failure or end the service life because system performance cannot meet requirements. So, oil contamination will greatly increase the control difficulty, deteriorate the control performance of the tank, and ultimately lead to decreasing equipment reliability and shortening durability, even which can cause casualties and loss of initiative during war. Therefore, the control research of vehicle mounted and unbalanced barrel pitching system under oil contamination condition can greatly improve the reliability of tank and enhance the actual combat capability, which is a major strategic issue related to safety and stability.

Around the control of barrel pitching system, scholars have done extensive research, and the result shows that the unbalanced torque is an important reason for the deterioration of control performance. Currently, the method of balancing such imbalances can be attributed to mass balance, balancing machine balance, and adding drive source balance [4]. Among them, the required balance quality for mass balance is large, which is difficult to achieve because of compact and light-weight requirements. Balancing machine balance can achieve light-weight requirements to a certain extent, but it is difficult to achieve compact requirements because of additional mechanical configuration. Therefore, Gao et al. [5] gave a new type of servo system to balance and control the unbalanced barrel, which uses the balance chamber of the three-chamber hydraulic cylinder to compensate the unbalanced force, and the two drive chambers are used for position and interference suppressing, which can achieve precise control of large diameter unbalanced barrel pitching system. Aiming at overcoming road fluctuation’s influence, Zhao [6] established a dynamic stability model of the barrel’s direct stability system, and analyzed the stability accuracy of the system to compensate the road fluctuation, which has important theoretical value.

In addition, the barrel pitching system is sensitive to oil contamination, and the oil contamination degree has strict requirements. Because of the existence of servo valve, the required oil level of barrel pitching system is no lower than NAS 7. In the service process, the hydraulic system is inevitably mixed with solid particles, gases, and liquids, which will reduce the reliability and increase the control difficulty of barrel pitching system. However, the existing methods pay more attention to active control of oil contamination such as purification and contamination detection, and there is no research on adaptive compensation control for system performance degradation [7].

In conclusion, there is no research on the control of barrel pitching system integratedly considering the effects of road fluctuation, unbalanced torque, and oil contamination. In this work, the overall model of barrel pitching system is established and its sensitivity analysis is conducted, and then the correction link is designed to overcome the influence of road fluctuation, and the pressure difference feedback loop is constructed to overcome the influence of changing unbalanced torque during pitching process. After that the oil contamination influence has been categorized based on the difference of contamination source, and then the core decay parameters induced by oil contamination are determined. At last the PID controller, the fuzzy PID controller, and the fractional order fuzzy PID controller are designed to restrain the complex influence of oil contamination, road fluctuation, and unbalanced torque.

#### 2. Working Principle and Modeling of the Barrel Pitching System

##### 2.1. Working Principle

The vehicle mounted and unbalanced barrel pitching system is actually an electro-hydraulic servo control system, the system structure is shown in Figure 2. The speed gyroscope detects pitching speed of barrel in real time and sends it to the controller together with its integral signal. The controller compares it with the given signal and outputs a control voltage signal, and then turns it into a current signal by the servo amplifier to control the proportional valve. The servo valve outputs flow and drives the hydraulic cylinder to move, and it finally drives the barrel to pitch around fulcrum O. In which *T*_{L} is the unbalanced torque generated by the barrel relative to the fulcrum O. The barrel pitching servo system is mounted on the vehicle, which will sway follow the road fluctuation, and the sensor measures the pressure of hydraulic cylinder’s two chambers, which is corrected and sent to servo amplifier for suppressing the influence of unbalanced torque. The gyroscope sensor detects the pitching angle in the vertical direction of the carrier in real time, and then the angle is corrected and sent to servo amplifier for suppressing the influence of road fluctuation.

##### 2.2. Modeling of the Barrel Pitching System

###### 2.2.1. Proportional Servo Valve Model

The servo amplifier converts received control voltage *U*_{1}(*t*) into output current *I*(*t*), and its mathematical model can be approximated to the proportional component shown as follows:

The operating frequency of the selected proportional servo valve is much higher than that of the actuator’s, so its mathematical model can be expressed as follows:where *Q*_{o}(*t*) denotes the no-load flow.

Double closed loop feedback of angular position and angular velocity are adopted. The relationship between given voltage *U*_{r}(*t*) of barrel pitching system and received voltage *U*_{1}(*t*) of servo amplifier can be expressed as follows:where *k*_{l} and *k*_{v} denotes the feedback coefficient of angular position and angular velocity and *θ* denotes the input angle of position sensor and velocity sensor.

###### 2.2.2. Variable Pump and Hydraulic Cylinder Model

Ignoring the smaller return oil pressure, the continuous flow equation for variable pump and hydraulic cylinder is as follows:where *p*_{1} denotes the inlet pressure, *C*_{t} denotes the leakage coefficient, *A* denotes the effective area of piston, *L* denotes the displacement of piston rod, denotes the single chamber volume of cylinder, and *β*_{e} denotes the volume elastic modulus of oil.

In Figure 2, set *θ*_{1} as the pitching angle of barrel relative to the vehicle and *θ*_{2} as the rotation angle of vehicle relative to the inertial space around the rotation axis. The rotation angle of barrel relative to the inertial space around rotation axis is as follows:

Since the pitching angle *θ*_{1} of barrel is the control variable, it is necessary to establish the relationship between the output displacement of piston rod *l* and the pitching angle *θ*_{1}. Set the length of AB in Figure 2 as follows:where *l*_{1} denotes the length of OA, *l*_{2} denotes the length of OB, *l*_{3} denotes the length of AB when the rotation angle of barrel is 0 degree, and *γ* + *θ*_{1} is the angle between OA and OB. Take the derivative of equation (6) about *t*, it can be obtained as

Based on equations (3)–(5) and (7), the continuous flow equation of hydraulic pump and cylinder can be obtained as follows:

The balance equation between cylinder and load is as follows:where *M* denotes the output torque of the hydraulic cylinder, *J* denotes the rotational inertia of barrel relative to the trunnion, *B*_{m} denotes the viscous friction coefficient, and denotes the unbalanced torque.

The output torque *M* of the hydraulic cylinder is the product of output force versus arm of output force relative to barrel rotation center:

###### 2.2.3. State Space Model

Based on equations (8)–(10) and assuming , , and , the state space model of barrel pitching system can be obtained as equation (11), from which the factors affecting control performance include external disturbance and internal parameter disturbance. The former includes unbalanced torque and road fluctuation *θ*_{2}, and the latter includes parameter change induced by working conditions such as oil pollution. Therefore, firstly the sensitivity analysis of each parameter to system performance should be analyzed to find out the more sensitive factors. Secondly, the isolation design of external disturbances should be carried out. Finally, an effective controller should be designed to effectively suppress the sensitive factors and external disturbances.where ,

#### 3. Sensitivity Analysis of the Barrel Pitching System

The sensitivity analysis of the barrel pitching system is carried out using the first order trajectory sensitivity model [8]. Trajectory sensitivity is to study the change of state vector caused by the change of parameters according to the state equation. The state equation of the barrel pitching system can be expressed as follows:where **x** is *N* dimensional state vector, **u** is *R* dimensional input independent of **α**, and *t* is time.

The solution of equation (13) is expressed as follows:

The sensitivity function of vector **x** to parameter **α** and its initial condition are defined as

In this work, **u** and **α** are independent, taking the derivative of equation (13) on both sides of **α** as

Select related vector of barrel pitching system aswhere the state vector parameters of barrel system are , , and ; the input vector parameters are *u*_{1} = *U*_{r}, *u*_{2} = *T*_{g}, and *u*_{3} = *θ*_{2}; and the independent vector parameters are *α*_{1} = *K*_{a}, *α*_{2} = *K*_{1}, *α*_{3} = *A*, *α*_{4} = *L*, *α*_{5} = *C*_{t}, *α*_{6} = *V*_{0}, *α*_{7} = *β*_{e}, *α*_{8} = *B*_{m}, *α*_{9} = *J*, *α*_{10} = *k*_{l}, and *α*_{11} = *k*_{v}.

Then the state space model of barrel pitching system can be expanded into the following equation:

The partial derivative of state equation function to the state vector **x** is conducted, and then the coefficient term of sensitivity function can be obtained as equation (19).

The partial derivative of state equation function to the state vector **α** is conducted, and then the free term of sensitivity function can be obtained as equation (20).

The initial value of state vector **x** is **x**_{0} = 0, and then the initial value of sensitivity function is as follows:

The trajectory sensitivity equation (16) is a first order linear non-homogeneous differential equation with time varying coefficients. There are *n* × (*p* + 1) equations to be solved, including *n* state equations and *n* × *p* sensitivity equations. In this work, only the sensitivity of barrel’s pitching angle *x*_{1} to the parameter vector **α** is considered.

The simulation model of barrel pitching system is established in MATLAB/Simulink. The initial value of parameters is shown in Table 1.

Figure 3 is the step response of barrel pitching system, the setting time is more than 2 s, the steady state curve is sine curve and its magnitude is almost 0.007 rad, and the steady state error is 0.001 rad, which cannot satisfy the requirements.

Based on the simulation model, the differential equation is solved in the command interface of MATLAB using the fourth-order/fifth-order Runge–Kutta algorithm, and then the sensitivity function of step response *x*_{1} to the parameter **α**_{i} (*i* = 1, 2, …, 11) is shown in Figure 4.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(g)**

**(h)**

**(i)**

**(j)**

**(k)**

Comparing Figures 3 and 4, it can be seen that the parameters *α*_{1}, *α*_{2}, *α*_{3}, *α*_{4}, *α*_{6}, *α*_{7,}*α*_{8}, *α*_{9}, and *α*_{11} only have a greater influence during the setting time; the parameters *α*_{5} and *α*_{10} have an effect on the output angle during all sampling time. And the sinusoidal curves with large magnitude are induced by the road fluctuation, which needs further suppression. In order to quantify the influence of each parameter, the percentage of output angle’s change Δ*x*_{1} relative to the angle’s step value *sx*_{1} is shown in Figure 5, and the maximum influence value of each parameter on the output angle is shown in Figure 6. It is clear that the values of *α*_{5} and *α*_{10} are small, although the values of *α*_{1}, *α*_{3}, *α*_{4}, *α*_{6}, *α*_{9,} and *α*_{11} are no less than 5%. Actually, the change of parameter’s value can be neglected after the barrel pitching system is determined, so it is necessary to carry out impact analysis and control of *α*_{2}, *α*_{7,} and *α*_{8}, which have great influence and will change with oil pollution.

#### 4. External Interference and Its Correction

##### 4.1. Road Fluctuation and Its Isolation Model

Based on the model of barrel pitching system, the control diagram can be obtained as shown in Figure 7. *G*_{1}(*s*) is the corrective component, and its function is to compensate the influence of *θ*_{2}*Js*^{2}. The control command *U*_{r} and unbalanced torque *T*_{g} are both set as 0, and the transfer function of output angle *θ* with respect to the road fluctuation signal *θ*_{2} is obtained aswhere *F*_{1}(*s*), *F*_{2}(*s*), and *F*_{3}(*s*) are denoted as follows:

To compensate for the effect of *θ*_{2}*Js*^{2}, we can get the following equation:

##### 4.2. Unbalanced Torque and Its Suppression Model

If *T*_{g0} is the unbalanced torque at 0 degree, the unbalanced torque *T*_{g} varies with the pitching angle:

In Figure 7, *G*_{2}(*s*) is the pressure feedback loop that suppresses the influence of unbalanced torque. The control command *U*_{r} and road fluctuation *θ*_{2} are both set as 0, the transfer function of the output angle *θ* with respect to unbalanced torque *T*_{g} is obtained as equation (26).

Combining equation (26) and the parameters’ value, it is found that the influence of unbalanced torque’s change on steady state accuracy is mainly derived from equation (27).

So the transfer function of pressure feedback loop can be obtained as follows:

#### 5. Influence of Oil Contamination on Barrel Pitching System

The failures of hydraulic system caused by oil contamination can directly affect its reliability and durability, where oil contaminant refers to substances that are mixed into the oil and have a harmful effect on the performance and service life of control system, which mainly includes solid particles, gases, and liquids.

The failures caused by oil contamination can be divided into two classes: sudden failures and gradual failures. The latter is defined as a condition when a performance indicator cannot be maintained within a specified range [9]. Under coupled action of pollutants, the technical state such as viscosity and elastic modulus will change, also the structural parameters such as the throttling edge of spool valve will be degraded. Before final failure, the control method proposed can delay the gradual degradation process.

##### 5.1. Solid Particles

Solid particle contamination is the most common and hazardous pollutant in hydraulic system. The oil flow with solid particle at a certain speed and angle will cause wear on the internal surface of material. The wear is mainly distributed in spool valve of the barrel pitching system, and the throttle model of spool valve is shown in Figure 8.

The radial clearance Δ and the throttle edge radius *r*_{1} and *r*_{2} will increase because of wear, so the equivalent opening amount of the same spool displacement *x*_{v} becomes as large as equation (29), which brings about the flow coefficient increasing sharply and stabilizing at a certain value as the spool displacement continues to increase. So, the flow coefficient *K*_{1} will increase:

##### 5.2. Gases

At present, the influence of gas pollutant in the hydraulic system is often neglected. The gas is mixed into oil, and it will affect the pressure, viscosity, temperature, and density, but mainly affect the elastic modulus of the oil. The experimental data shows the oil entrains 1% of gas; the elastic modulus is reduced to 35.6% of pure oil, which will seriously affect the working performance. The relationship between elastic modulus *K*_{e} and gas content of oil is shown as equation (30), from which we find the elastic modulus of oil drops sharply when the gas content increases.where *K*_{f} denotes the elastic modulus of gas-free oil; *K*_{c} denotes the elastic modulus of hydraulic pipe system; *K*_{b} denotes the elastic modulus of gas contained in oil; *V*_{b} denotes the volume of gas contained in oil; and *V*_{f} denotes the volume of oil and gas mixture.

##### 5.3. Liquids

The performance parameters will not have significant change when the liquid is mixed into oil alone, but can raise the oil temperature together with solid particles and gases, and the viscosity of oil is closely related to oil temperature. After research, Jiang [10] concluded that the relationship between viscosity and oil temperature satisfies the following empirical equation:where *μ*_{x} denotes the viscosity of oil when the pressure is *P* and the temperature is *t*; *μ*_{0} denotes the viscosity of oil when the pressure is standard atmospheric pressure and the temperature is *t*_{0}; *α* denotes the viscosity coefficient; and *λ*_{x} denotes the proportional constants of oil to temperature.

According to the analysis result, when the oil temperature changes from 0° to 150°, the oil viscosity is relatively small, which decreases slowly with the increase of temperature. That is, the oil viscosity will decrease and then lead to the dumping coefficient’s decrease of spool core under the action of liquid together with other pollutants.

#### 6. Fractional Order Fuzzy PID Controller

In view of the parameters’ change induced by oil contamination of barrel pitching system, a fractional order fuzzy PID controller is designed for the control of barrel pitching system. The block diagram of fractional order fuzzy adaptive controller is shown in Figure 9.

##### 6.1. Fuzzification

###### 6.1.1. Definition of Variables and Domains

The input variables of the controller are the pitching angle deviation *e* and the deviation change rate *e*_{c}, and the output variables are the parameters’ correction amount Δ*k*_{p}, Δ*k*_{i}, and Δ*k*_{d}. The domain of deviation *e* and deviation change rate *e*_{c} are [−0.5, 0.5] and [−5, 5], the domain of Δ*k*_{p} is [−500, 500], the basic domain of Δ*k*_{i} is [−5, 5], and the basic domain of Δ*k*_{d} is taken as [−5, 5]. Each variable domain is discretely quantized into 7 levels using NB, NM, NS, ZE, PS, PM, and PB.

###### 6.1.2. Selection of Membership Function

The membership function reflects the relationship between the exact value and the fuzzy value. Both the input and output variables use the triangle membership function.

###### 6.1.3. Fuzzy Rule Base

Based on the knowledge accumulated, the logical relationship between input and output is obtained. Tables 2–4 are the parameter tuning rule bases.

##### 6.2. Fuzzy Inference

The fuzzy reasoning method can be divided into Mamdani reasoning method, Larsen reasoning method, and Zadeh reasoning method. This work selects the most widely used Mamdani reasoning method, and its fuzzy implication relationship is relatively simple, which is obtained by taking the Cartesian product of fuzzy set and .

##### 6.3. Defuzzification

Defuzzification is to convert the fuzzy control quantity obtained by fuzzy inference into the exact quantity actually used for control. In this work, the center of gravity method is used to defuzzify.where *μ*_{Δkpj} denotes the membership degree of Δ*k*_{p} and the exact values of Δ*k*_{i} and Δ*k*_{d} can be obtained.

#### 7. Simulation and Discussion

Based on the given angle signal, the barrel is driven by the electrohydraulic servo system to perform pitching motion, which is detected by the sensor in real time. The PID controller, the fuzzy PID (FPID) controller, and the fractional order fuzzy PID (FOFPID) controller are designed to overcome the influence of unbalanced torque, road fluctuation, and oil contamination and finally achieve high precision pitching control of vehicle mounted and unbalanced barrel. The simulation model is built using MATLAB, and the value of simulation parameters is shown in Table 5.

##### 7.1. Simulation of FOFPID Control under Initial Value of Parameters

Based on the FOFPID controller, the barrel pitching angle’s step response curve under initial value of parameters is shown in Figure 10, in which the given angle is 0.436 (25°) rad, and *λ* and *μ* are set as 0.78 and 0.84. As observed, the dynamic response under initial value of parameters has been improved. The FOFPID controller reduces the setting time to 0.2 s. At the same time, the magnified steady state curve is still sine curve instead of straight line, but the magnitude of sine curve and the steady state error are greatly reduced, which is far less than that without control. To summarize, the FOFPID control method can fully satisfy the high precision control requirements under initial unbalanced torque, road fluctuation, and clean oil condition.

##### 7.2. Contrast Simulation under Bigger Unbalanced Torque

Actually, the unbalanced torque *T*_{g} will increase because of barrel length’s increasing or center of gravity’s changing. The bigger unbalanced torque is set as 5000 Nm, then the PID controller, the FPID controller, and the FOFPID controller are separately applied to the barrel pitching system, and the contrast simulation curves are shown in Figure 11, from which it is clear to find the dynamic characteristics under FPID is slightly better than the PID’s, but the difference is not obvious, and then the FOFPID’s is far better than them. Simultaneously, the magnified steady state curves are all sine curves, but the steady state error is 0.000038 rad, which is much less than the others; that is, the correction link and the FOFPID control can effectively suppress the influence of unbalanced torque’s increase.

##### 7.3. Contrast Simulation under Bigger Magnitude of Road Fluctuation

The magnitude of road fluctuation *θ*_{2} will increase because of road level’s changing. The bigger unbalanced torque is set as 5° (0.087 rad), then the PID controller, the FPID controller, and the FOFPID controller are separately applied to the barrel pitching system, and the contrast simulation curves are shown in Figure 12, from which it is clear to find that the dynamic characteristics under FPID is slightly better than the PID’s but the difference is not obvious, and then the FOFPID’s is far better than them. Simultaneously, the magnified steady state curves are all sine curves, but the steady state error is 0.000132 rad, which is less than the others; that is, the pressure difference negative feedback loop and the FOFPID control can effectively suppress the influence of road fluctuation magnitude’s increase.

##### 7.4. Contrast Simulation under Bigger Flow Coefficient

The flow coefficient *K*_{1} will increase because of oil contamination. The bigger flow coefficient is set as 0.131 m^{3}/(s·A), then the PID controller, the FPID controller, and the FOFPID controller are separately applied to the barrel pitching system, and the contrast simulation curves are shown in Figure 13, from which it is clear to find that the dynamic characteristics under FPID is slightly better than the PID’s but the difference is not obvious, and then the FOFPID’s is far better than them. Simultaneously, the magnified steady state curves are all sine curves, but the steady state error is 0.00003 rad, which is much less than the others; that is, the FOFPID control can effectively suppress the influence of flow coefficient’s increase.

##### 7.5. Contrast Simulation under Smaller Elastic Modulus

The elastic modulus *β*_{e} will decrease because of oil contamination. The smaller elastic modulus is set as 100 MPa, then the PID controller, the FPID controller, and the FOFPID controller are separately applied to the barrel pitching system, and the contrast simulation curves are shown in Figure 14, from which it is not difficult to find that the dynamic characteristics under FPID is far better than the PID’s, and then the FOFPID’s is slightly better than the FPID’s. Similarly, the FOFPID’s steady state performance is slightly better than the FPID’s but far better than the PID’s. Then the steady state error is 0.00009 rad; that is, the FOFPID control can effectively suppress the influence of elastic modulus’s decrease.

##### 7.6. Contrast Simulation under Smaller Viscosity Coefficient

The viscosity coefficient *B*_{m} will decrease because of oil contamination. The smaller viscosity coefficient is set as 20 Nm/(rad/s), then the PID controller, the FPID controller, and the FOFPID controller are separately applied to the barrel pitching system, and the contrast simulation curves are shown in Figure 15, from which it is clear to find that the dynamic characteristics under FPID is slightly better than the PID’s but the difference is not obvious, and then the FOFPID’s is far better than them. Simultaneously, the magnified steady state curves are all sine curves, but the steady state error is 0.00005 rad, which is much less than the others; that is, the FOFPID control can effectively suppress the influence of viscosity coefficient’s decrease.

#### 8. Conclusions

In this work, to accomplish the high precision control, the sensitivity analysis of the barrel pitching system is conducted, and the correction link and pressure difference negative feedback loop are respectively designed for suppressing the influence of unbalanced torque and road fluctuation. Also the fractional order fuzzy PID controller has been presented in order to compensate the core parameters’ decay induced by the oil contamination and the external interference.

As elaborated, in addition to dynamic characteristic improvement, the setting time is less than 0.2 s. The control method leads to the decrease of steady state error so that the maximum value is 0.000132 rad under external interference and oil contamination conditions. And compared with the PID controller and the fuzzy PID controller, the fractional order fuzzy PID can adjust control parameters in order to adapting to different conditions, so the presented control method is suitable for the control of vehicle mounted and unbalanced barrel pitching system.

Finally, based on the results of this work, an approach for the future research should include the considering of the vehicle vibration’s impact on the barrel pitching servo control system.

#### Data Availability

The data used to support the findings of this study have not been made available because the data also forms part of an ongoing study.

#### Conflicts of Interest

The author declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) (grant number 51805403), the project supported by Natural Science Basic Research Plan in Shaanxi Province of China (grant number 2019JQ-836), the Natural Science Foundation of Shaanxi Provincial Department of Education (grant number 19JK0416), and the Dean Fund of Photoelectric Institute (grant number 2017GDYJY01).