Multiple Fusion Based on the CCD and MEMS Accelerometer for the Low-Cost Multi-Loop Optoelectronic System Control

In the charge-coupled device (CCD) and micro-electro-mechanical system (MEMS) accelerometer based low-cost multi-loop optoelectronic control system (OCS), due to accelerometers’ drift and noise in low frequency, the disturbance suppression (DS) is insufficient. Previously, based on the acceleration and position dual-loop control (ADLC), researchers combined a disturbance observer (DOB) with a virtual velocity loop to make some medium-frequency DS exchange for low-frequency performance. However, it is not optimal because the classic DOB based on accelerometers’ inaccurate signals cannot observe accurate disturbance in low frequency and the velocity based on a CCD and accelerometer time-domain fusion carried the CCD’s delay, resulting in the drop of medium-frequency DS. In this paper, considering the CCD’s advantage in low frequency and the accelerometer’s strength in high frequency, we propose to fuse their signals twice with a modified complementary filter method to respectively acquire an acceleration and velocity. The new acceleration with no drift and less noise but lower bandwidth creates a new acceleration model and is only used in fusion DOB (FDOB), while the velocity with little delay is to build an additional velocity loop. Compared with the traditional DOB enhanced by the time-domain fusion velocity loop, experiments verify that the proposed multiple fusion would apparently enhance the system’s DS, especially in low and medium frequency.


Introduction
The CCD-based optoelectronic control system (OCS) is an important component in the high-precision capture and tracking optoelectronic platform, which is commonly used in astronomical observation, laser communication, target tracking and line-of-sight stabilization [1][2][3][4]. With the expansion of generalization and commercialization, the OCS becomes smaller and is more installed on moving platforms, such as vehicles, airplanes, satellites, which are susceptible to external disturbance. Limited by the CCD's delay and low-sampling rate, the disturbance suppression(DS) performance of the outer position loop is extremely restricted [5,6]. Generally, to enhance the DS, we attempted to use inertial sensors to build a high-sampling rate inner loop [7,8]. As a member of the inertial sensors family, the MEMS accelerometer has advantages with low price, small size, high precision and Section 4 discusses the modified complementary filter and the details of the fusion acceleration and velocity. Section 5 is the experiment part which respectively presents the DS improvement from the FDOB and the virtual velocity loop. Concluding remarks are presented in Section 6.

The EDOB Built in the Acceleration Control Loop
DOB as a disturbance compensation feed forward method, has commonly been used in industrial control [20,21]. As it does not need to add an additional sensor and basically does not affect the stability of the system, it is very suitable for the small OCS whose model can be recognized by spectrum fitting and seldom changes in motion. Deng first introduced the DOB disturbance feed forward control into the OCS and later proposed a modified method called EDOB [13,22]. Figure 1 presents the structure of EDOB. Compared with the DOB, EDOB changed the compensation object G a to be the open-loop transfer function C a G a and optimized the controller design. improvement from the FDOB and the virtual velocity loop. Concluding remarks are presented in Section 6.

The EDOB Built in the Acceleration Control Loop
DOB as a disturbance compensation feed forward method, has commonly been used in industrial control [20,21]. As it does not need to add an additional sensor and basically does not affect the stability of the system, it is very suitable for the small OCS whose model can be recognized by spectrum fitting and seldom changes in motion. Deng first introduced the DOB disturbance feed forward control into the OCS and later proposed a modified method called EDOB [13,22]. The closed-loop acceleration is given as follows.
As u is a substitutable factor, after substitution, we get  The closed-loop acceleration is given as follows.
As u is a substitutable factor, after substitution, we get Since we only focus on the line-of-sight stabilization, a re f could be set as 0. The acceleration loop's DS transfer function is as Equation (4).
Theoretically, in Equation (4), if C f is set to be (C a G a ) −1 , E EDOB would become 0. However, due to model mismatch and noise effects, the DS improvement is limited, especially in low frequency. The reason is that the observed disturbance in low frequency from the MEMS accelerometer is not accurate with a big drift and much noise and would lead to poor effects of compensation. In order to get a satisfying low-frequency DS, we first need to acquire a relatively accurate acceleration signal depending on the existing sensors, decreasing the influence of the drift and noise.

The FDOB Control
In OCS, the CCD's low-frequency signal is very valuable. On the contrary, although the accelerometer has a bad low-frequency signal, it has a high bandwidth (around 1000 Hz) nearly with no delay. It is reasonable to combine the CCD with the MEMS accelerometer to get a fusion acceleration. The basic principle of the proposed complementary filter is shown in Figure 2. Theoretically, in Equation (4), if f C is set to be ( ) , EDOB E would become 0. However, due to model mismatch and noise effects, the DS improvement is limited, especially in low frequency. The reason is that the observed disturbance in low frequency from the MEMS accelerometer is not accurate with a big drift and much noise and would lead to poor effects of compensation. In order to get a satisfying low-frequency DS, we first need to acquire a relatively accurate acceleration signal depending on the existing sensors, decreasing the influence of the drift and noise.

The FDOB Control
In OCS, the CCD's low-frequency signal is very valuable. On the contrary, although the accelerometer has a bad low-frequency signal, it has a high bandwidth (around 1000 Hz) nearly with no delay. It is reasonable to combine the CCD with the MEMS accelerometer to get a fusion acceleration. The basic principle of the proposed complementary filter is shown in Figure 2.   Figure 3. From Figure 2, we get The closed-loop transfer function B closed is a low-pass filter which is to extract the low-bandwidth sensor's signal, while the restrain transfer function B restrain is a high-pass one to filter the high-bandwidth sensor's signal. If H low−bandwidth = H high−bandwidth , the fusion output is equal to the real motion status. However, even if H low−bandwidth = H high−bandwidth , it does not matter because we can also calculate the real model with spectrum fitting. This is equivalent to rebuilding the platform's model and producing a virtual sensor.
Suppose the platform's new acceleration model is F a (s) under the complementary filter method, F a (s)/G a (s) represents the transfer characteristic of the virtual accelerometer. Depending on F a (s), we build a new DOB structure called FDOB, which is shown in Figure 3.
The closed-loop acceleration is given as follows After substitution, Similarly, we set a re f = 0, the acceleration loop DS of the FDOB method is as follows.
Compared with Equation (4), the compensation object has changed to be F a (s). Different from directly using the accelerometer's data, the fusion acceleration is better with no drift and less noise in low frequency, which means we can observe a more accurate outer disturbance and the low-frequency DS would be much more improved than by the EDOB method.
Unfortunately, the FDOB as a feedforward way would not solve all the problems. Due to the existing of residual noise, the effect of compensation would be poor and even terrible in lower frequency, which is faced by all the feedforward methods. Considering that a feedback loop could get a converged result and naturally suppress noise's influence, we attempt to build a virtual velocity loop between the acceleration loop and the position loop. The closed-loop acceleration is given as follows After substitution, Similarly, we set 0 ref a = , the acceleration loop DS of the FDOB method is as follows.
Compared with Equation (4), the compensation object has changed to be ( ) a F s  . Different from directly using the accelerometer's data, the fusion acceleration is better with no drift and less noise in

The Fusion Virtual Velocity Loop
After a virtual velocity loop is added, the ADLC system turns to be a triple-loop control(TLC) structure which is exhibited in Figure 4. The fusion method in Reference [14] is based on time domain, which would bring the CCD's time delay to the velocity no matter in low frequency or high frequency and caused a big drop of DS in relatively high frequency. The proposed complementary filter method is based on frequency domain. Since we combine the low-frequency signal of the CCD which is less affected by the delay with the high-frequency signal of the MEMS accelerometer, we minimize the impact of delay on the system. The whole DS without and with a virtual velocity loop is respectively as follows.   Similarly, Generally, we choose to design the open-loop transfer function to be an integral element in order to get an adequate phase and amplitude margin. So if the controllers can be designed ideally, , which means the whole DS would be enhanced in all frequency domain. Similarly, Generally, we choose to design the open-loop transfer function to be an integral element in order to get an adequate phase and amplitude margin. So if the controllers can be designed ideally, obviously, |1/(1 + C v · (1/s))| < 1. And |E TLC−FDOB | would be smaller than |E ADLC−FDOB |, which means the whole DS would be enhanced in all frequency domain.
As we know, the system's robustness is a key indicator of the control performance. And it refers to the system's sensitivity to parameter changes. Since the ADLC and TLC-FDOB both have a same position loop, we just need to examine the velocity sensitivity transfer function.
The sensitivity functions of ADLC and TLC-FDOB are as follows.
where Apparently, according to the above analysis, |S TLC−FDOB | < |S ADLC |, which signifies that the FDOB and additional virtual velocity loop could also improve the robustness of the traditional ADLC structure. If the parameters of the controlled object change a little because of external interference, the system's stability of the TLC-FDOB would suffer less than that of the ADLC.

The Fusion Acceleration Based on the Modified Complementary Filter Method
According to Equation (5), considering the CCD's delay e −τs , we can get the transfer function of the virtual accelerometer as follows. H CCD respects the transfer characteristic of the CCD without delay, H ACC is the transfer characteristic of the MEMS accelerometer.
When in very low frequency, B restrain approaches to 0, a f usion is equal to the acceleration got through the CCD with no drift. Similarly, when in relatively high frequency, B closed−loop approaches to 0, the a f usion is equal to the acceleration got by the MEMS accelerometers with nearly no delay. In medium frequency, a f usion contains both the CCD's and MEMS accelerometers' signals.
However, in fact, although the original complementary filter method can eliminate the low-frequency drift, it still cannot get a very satisfied acceleration because of the MEMS accelerometer's high-frequency noise.
Assume R CCD and R ACC are respectively the CCD's signal without delay and the MEMS accelerometer's signal and they contain both useful signal and interference signal. According to Equation (16), we can deduce the fusion acceleration as follows.
The new fusion acceleration contains two items, the R ACC and the low-frequency data filtered by C low−pass . When given a low-frequency input, R ACC not only contains a zero drift but also contains high-frequency noise. The second item of Equation (17) could help eliminate the drift of R ACC but cannot deal with the high-frequency noise signal, which will decrease the accuracy of the last acceleration. Therefore, in this case the basic complementary filter cannot work well. In order to prevent the pollution of the high-frequency noise, we change C high−pass to be a band-pass filter C band−pass .
Considering that the high-order low-pass filter would bring a large phase lag, we choose a first-order low-pass filter as C low−pass . Equation (18) presents the new fusion acceleration.
In Equation (18), although the double differential of R CCD to get the acceleration would amplify the impact of the noise, since the signal of noise's differential mainly distribute on high-frequency domain, C low−pass will nearly eliminate the effect of the noise. Similarly, the band-pass filter C band−pass could cut down the low-frequency drift and decrease the pollution of the high-frequency noise. Hence, now we could ignore the drift and noise and get R CCD s 2 ≈ R ACC ≈ a real , then In Equation (19), the form of the acceleration transfer function has been fixed. The rest is to determine the handover frequency 1/T 1 . Theoretically, whatever the value of 1/T 1 is, the final acceleration model could be identified by spectrum fitting. Nevertheless, in fact, 1/T 1 will affect the linearity of the object, which will determine the difficulty of the controller design. Figure 5 exhibits the simulation of the acceleration bode response with different 1/T 1 .
In Figure 5, as 1/T 1 increases, the curve fluctuates more and more violently. It is because the smaller handover frequency means less influence of the non-linear link e −τs on the fusion acceleration. Since the curve fluctuations are detrimental to controller design, we should try to choose a small handover frequency. However, since the drift and noise of the MEMS accelerometer have a serious pollution on the useful signal below 2 Hz, to reduce the MEMS accelerometer's proportion in the fusion acceleration, 1/T 1 should be no less than 2 Hz. Therefore, we choose 2 Hz as the handover frequency.
Since the non-linear link e −τs is not convenient for spectrum analysis, we give out its approximate form.
Substituting Equation (20) into Equation (19) a f usion (s) ≈ ( 1  Equation (21) is the approximation of Equation (19). Intuitively, the rebuilt acceleration object a F is to add an inertial element behind the original one. However, in fact, purely adding an inertial element could only filter the high-frequency noise but the complementary filter method would simultaneously eliminate the time-domain drift. Considering that the proposed way would also lead to a decrease of the acceleration's bandwidth, in order to avoid the decrease of the acceleration Equation (21) is the approximation of Equation (19). Intuitively, the rebuilt acceleration object F a is to add an inertial element behind the original one. However, in fact, purely adding an inertial element could only filter the high-frequency noise but the complementary filter method would simultaneously eliminate the time-domain drift. Considering that the proposed way would also lead to a decrease of the acceleration's bandwidth, in order to avoid the decrease of the acceleration closed-loop bandwidth, in the acceleration feedback loop we still use the MEMS accelerometers' original signal. Since an accurate model in high frequency is commonly unreliable, with DOB method we generally choose to compensate low-and medium-frequency disturbance. What's more, as a feed forward method, DOB would not affect the closed-loop performance. Therefore, the relatively low-bandwidth fusion acceleration is very suitable to the DOB method.

The Fusion Velocity
Similarly, according to Equation (5), we can differentiate the CCD's signal and integrate the accelerometer's signal to get the fusion velocity as follows.
In Equation (22), we still choose a one order filter to extract the signals of the CCD and the MEMS accelerometer. Since the fusion velocity is used in closed-loop control, we choose a high-pass filter not a band-pass to extract the MEMS accelerometer's signal, or it will result in the decrease of closed-loop bandwidth. Substitute the filters to Equation (22) as follows.
In Equation (23), we can treat sH CCD ≈ 1 Substitute Equation (20) into Equation (24) From Equation (25), the fusion velocity would have the same form with the platform's real velocity but in fact, due to the influence of the delay, the natural frequency and damping coefficient may change slightly. The handover frequency 1/T 1 could be 2 Hz, which is same with the fusion acceleration's. With the spectrum fitting method, we can identify the detail parameters and give out the accurate transfer function. Figure 6 exhibits the experimental apparatus. As the OCS is a two-axis symmetrical system, we only take one axis into consideration. To simulate the external disturbance, an additional OCS platform is used as the pedestal which can be driven by the dynamic signal analyzer. The disturbance signal is a sine wave with variable frequency. The laser light source emits a light as a reference. The CCD receives the reflected light from the OCS to calculate the offset error to the center of the boresight. Two linear MEMS accelerometers (Model 1221, SILICON DESIGNS, Inc., Kirkland, WA, USA) work in a differential configuration to get the angle acceleration of one direction. The motion of the disturbance platform is measured by the eddy fixed on its bottom. The CCD works in a 100 Hz sampling rate with 0.02 s delay, while the MEMS accelerometers works in a 5000 Hz rate. In the processing of measuring the DS, the disturbance platform works on an open-loop mode, continuously outputting sine signal from 1~100 Hz. The stabilization platform works on a closed-loop mode to stabilize the visual axis. Figure 7 presents the time-domain curves of the MEMS accelerometer and the virtual accelerometer in different frequencies. Below 2 Hz, the MEMS accelerometer's signal contains significant drift and noise, while the virtual accelerometer's signal is much better with no drift and less noise, which means the modified complementary is valid. Above 2 Hz, the signals of the MEMS accelerometer and the virtual one are both good. Although they may have different amplitudes in the same frequency, it does not matter because we can identify the last transfer function with spectrum fitting. Figure 8 exhibits their open-loop bode responses from 1 Hz to 1 KHz. The blue line represents the response measured from the real platform, while the red line represents the fitting one from the math model. Obviously, the matching degree of fitting is very high. And the specific identified transfer functions are as follows.

The FDOB Experiment Based on Acceleration and Position Dual-Loop Control
G a (s) = 0.0022s 2 0.007s 2 + 0.0185s + 1  Figure 7 presents the time-domain curves of the MEMS accelerometer and the virtual accelerometer in different frequencies. Below 2 Hz, the MEMS accelerometer's signal contains significant drift and noise, while the virtual accelerometer's signal is much better with no drift and less noise, which means the modified complementary is valid. Above 2 Hz, the signals of the MEMS accelerometer and the virtual one are both good. Although they may have different amplitudes in the same frequency, it does not matter because we can identify the last transfer function with spectrum fitting. Figure 8 exhibits their open-loop bode responses from 1 Hz to 1 KHz. The blue line represents the response measured from the real platform, while the red line represents the fitting one from the math model. Obviously, the matching degree of fitting is very high. And the specific identified transfer functions are as follows.  G a (s) is the highly approximate model of the original acceleration model based on the MEMS accelerometer, which has a nearly 1 KHz bandwidth and is very suitable for closed-loop feedback control. F a (s) with no more than 100 Hz bandwidth, is the rebuilt acceleration model based on the virtual accelerometer, which is more suitable for DOB control. To measure the closed-loop DS of the OCS, firstly, we should design the acceleration controller a C , which is usually designed as a lag controller.
After the closed-loop acceleration has been designed, we can add the FDOB. According to Equation (12), f C should be designed as follows. To measure the closed-loop DS of the OCS, firstly, we should design the acceleration controller a C , which is usually designed as a lag controller.
After the closed-loop acceleration has been designed, we can add the FDOB. According to Equation (12), f C should be designed as follows. To measure the closed-loop DS of the OCS, firstly, we should design the acceleration controller C a , which is usually designed as a lag controller.
After the closed-loop acceleration has been designed, we can add the FDOB. According to Equation (12), C f should be designed as follows.
From Equation (11), we know the improvement of DS is brought by the numerator, so we focus on its value after compensation.
Equation (31) is a high-pass filter as Figure 9. Compared with the traditional pure ADLC method, the low-frequency improvement is very apparent and as the frequency goes up, the promotion becomes smaller and smaller until there is no improvement, which is in line with forecasts. Figure 10 presents the measured improvement of DS brought by EDOB and FDOB to the pure ADLC method. Obviously, the DS of the pure ADLC is not satisfied. Although the EDOB method could enhance the DS in low and medium frequencies, the low-frequency improvement is insufficient because of the MEMS accelerometer's drift and much noise. Fortunately, the proposed FDOB has solved the problem. Compared with pure ADLC method, FDOB could extremely promote the low-frequency DSA and the maximum improvement could reach −30 dB in 1.8 Hz. Although the increase below 1.8 Hz has slowed down due to the effect of the residual noise, the improvement in low frequency is more apparent than in medium and high frequencies, which is coincident with the above simulation. In order to further promote the system's DS and make full of sensors' potential, we would continuously add the virtual velocity loop between the acceleration loop and the position loop. In Equation (29), since the order of numerator is higher than the one of denominator, it cannot be accomplished in physics. The actually used controller is presented as Equation (30).
From Equation (11), we know the improvement of DS is brought by the numerator, so we focus on its value after compensation.
Equation (31) is a high-pass filter as Figure 9. Compared with the traditional pure ADLC method, the low-frequency improvement is very apparent and as the frequency goes up, the promotion becomes smaller and smaller until there is no improvement, which is in line with forecasts. Figure 10 presents the measured improvement of DS brought by EDOB and FDOB to the pure ADLC method. Obviously, the DS of the pure ADLC is not satisfied. Although the EDOB method could enhance the DS in low and medium frequencies, the low-frequency improvement is insufficient because of the MEMS accelerometer's drift and much noise. Fortunately, the proposed FDOB has solved the problem. Compared with pure ADLC method, FDOB could extremely promote the low-frequency DSA and the maximum improvement could reach −30 dB in 1.8 Hz. Although the increase below 1.8 Hz has slowed down due to the effect of the residual noise, the improvement in low frequency is more apparent than in medium and high frequencies, which is coincident with the above simulation. In order to further promote the system's DS and make full of sensors' potential, we would continuously add the virtual velocity loop between the acceleration loop and the position loop.    Figure 11 presents the time-domain waves of the virtual gyro in different frequencies and we also present a group of signals of the fiber-optic gyro (FOG) (XW-FG70-20, Starneto, Beijing, China) as a comparison. The signal of the virtual gyro has a negligible phase lag and is very close to the real FOG's. It means the frequency-domain fusion is much better than the time-fusion fusion in Reference [14], which would lead to the signal's big phase lag at a relatively high frequency. Figure  12 shows the open-loop bode responses of the velocity from the FOG and virtual gyro. The identified transfer functions are as Equations (32) and (33). Their models are the same but the parameters are somewhat different.   Figure 11 presents the time-domain waves of the virtual gyro in different frequencies and we also present a group of signals of the fiber-optic gyro (FOG) (XW-FG70-20, Starneto, Beijing, China) as a comparison. The signal of the virtual gyro has a negligible phase lag and is very close to the real FOG's. It means the frequency-domain fusion is much better than the time-fusion fusion in Reference [14], which would lead to the signal's big phase lag at a relatively high frequency. Figure 12 shows the open-loop bode responses of the velocity from the FOG and virtual gyro. The identified transfer functions are as Equations (32) and (33). Their models are the same but the parameters are somewhat different.  Figure 11 presents the time-domain waves of the virtual gyro in different frequencies and we also present a group of signals of the fiber-optic gyro (FOG) (XW-FG70-20, Starneto, Beijing, China) as a comparison. The signal of the virtual gyro has a negligible phase lag and is very close to the real FOG's. It means the frequency-domain fusion is much better than the time-fusion fusion in Reference [14], which would lead to the signal's big phase lag at a relatively high frequency. Figure  12 shows the open-loop bode responses of the velocity from the FOG and virtual gyro. The identified transfer functions are as Equations (32) and (33). Their models are the same but the parameters are somewhat different.

The Virtual Velocity Loop with the FDOB
From the open-loop transfer functions, although the resonant frequency and damping coefficient of the virtual gyro has changed, its bandwidth is very high and fits the closed-loop feedback control. Now, the ADLC feedback control has turned to be the TLC structure. a C and f C could be the same as the previous design. Since the inner acceleration loop has improved the platform's characteristics, we can both use a low-pass filter as the velocity controller v C and the position controller p C .  Figure 13 shows the comparison of the DS with different fusion methods. Compared with the ADLC method enhanced by FDOB, the introduction of the virtual velocity loop with the time-domain fusion method in Reference [14] can enhance the low-frequency DS a little but the medium-frequency DS would decrease because of the CCD's time delay. However, the virtual velocity loop based on the frequency-domain fusion could apparently improve the low-and medium-frequency DS. Although the virtual velocity loop cannot enhance the high-frequency DS, it does not matter because it commonly relies on the mechanical design. In total, after fusing the sensors' data twice, the performance of the sensors has been fully utilized and we have obtained a satisfied system with extremely strong DS. Figure 14 presents the DS's comparison of the proposed method and Reference [14]. Figure 15 displays the residual stabilization errors of the two methods in different frequencies. Obviously, whether in frequency domain or time domain, the promotion is very big and could reach 10 dB in low and medium frequency, which means we have succeeded in further enhancing the system's DS. G v (s) = 2.3s 0.0072s 2 + 0.0202s + 1 From the open-loop transfer functions, although the resonant frequency and damping coefficient of the virtual gyro has changed, its bandwidth is very high and fits the closed-loop feedback control. Now, the ADLC feedback control has turned to be the TLC structure. C a and C f could be the same as the previous design. Since the inner acceleration loop has improved the platform's characteristics, we can both use a low-pass filter as the velocity controller C v and the position controller C p .
C v = 0.15 · 1 1 + 0.00077s , C p = 12 · 1 1 + 0.0005s (34) Figure 13 shows the comparison of the DS with different fusion methods. Compared with the ADLC method enhanced by FDOB, the introduction of the virtual velocity loop with the time-domain fusion method in Reference [14] can enhance the low-frequency DS a little but the medium-frequency DS would decrease because of the CCD's time delay. However, the virtual velocity loop based on the frequency-domain fusion could apparently improve the low-and medium-frequency DS. Although the virtual velocity loop cannot enhance the high-frequency DS, it does not matter because it commonly relies on the mechanical design. In total, after fusing the sensors' data twice, the performance of the sensors has been fully utilized and we have obtained a satisfied system with extremely strong DS. Figure 14 presents the DS's comparison of the proposed method and Reference [14]. Figure 15 displays the residual stabilization errors of the two methods in different frequencies. Obviously, whether in frequency domain or time domain, the promotion is very big and could reach 10 dB in low and medium frequency, which means we have succeeded in further enhancing the system's DS.

Conclusions
In the small and cost-efficient OCS platform, to save cost and space, we choose the MEMS accelerometer as the inertial stabilization sensor. At the same time, the proposed complimentary filter, with a small amount of calculation, saves expenses on computing devices. To fully release sensors' potential and enhance the DS as much as possible, this paper attempted to fuse signals of the MEMS accelerometer and CCD twice to get two virtual sensors without extra investment. We improved the fusion method, rebuilt the acceleration model and proposed a new FDOB structure that was not dependent on the original object. What is more, we skillfully combined the disturbance feed forward control and the feedback control to increase the system's robustness. Compared with the previous time-domain once fusion in Reference [14], the proposed frequency-domain fusion minimizes the effect of the accelerometers' defect and the CCD's time delay. These ideas are worth learning on other occasions. Experiments verified that the proposed method is valid.