#### Abstract

A novel method of fault line selection based on IOS is presented. Firstly, the IOS is established by using math model, which adopted TZSC signal to replace built-in signal of duffing chaotic oscillator by selecting appropriate parameters. Then, each line’s TZSC decomposed by db10 wavelet packet to get CFB with the maximum energy principle, and CFB was solved by IOS. Finally, maximum chaotic distance and average chaotic distance on the phase trajectory are used to judge fault line. Simulation results show that the proposed method can accurately judge fault line and healthy line in strong noisy background. Besides, the nondetection zones of proposed method are elaborated.

#### 1. Introduction

When SPG fault occurs in power distribution network, fault line selection method is always study emphasis. For the arc suppression coil grounding system, the fault current is very weak, and judging the fault line accurately becomes more difficult. Therefore, power distribution network in china is also called small current to ground system. For a long time, many scholars did a lot of studies about fault line selection and came up with many novel methods of different theories, which obtained some achievements [1–7].

The signal processing tools, such as wavelet transform (WT), S-transform, Prony method, and Hough Transform, which used to obtain fault signals feature in many methods. Then, many papers establish fault line selection criterion with artificial neural network [8] or support vector machine [9]. In [1], the paper adopted WT to study SPG fault, it used good localization ability in time and frequency domain to extract the characteristic of transient signal on different scales, but WT is easily influenced by noise. Besides, different wavelet-based functions have different results for characteristic of transient signal. In [2], TZSC data of first T/4 is used, and S-transform is carried out to determine the main characteristic frequency. Therefore, fault line selection result is detected by comparing S-transform energies of each feeder line. To avoid TAs density effect of magnetic saturation, paper [3] puts forward a piecewise Prony method; it improved the precision of the traditional Prony method, and the proposed method was used to match TZSC signals to get fault information. In [4], Hough transform was used to achieve fault line selection; actually, Hough transform has an excellent performance in the case of not exceeding two-dimensional parameter space, but if parameter space increased, calculation amount will rise sharply, meanwhile, costly storage space and time-consumption also will be surged. Paper [8] presents a method based on rough membership neural network to classify 10 kinds of fault types about transmission line. In [9], support vector machine (SVM) has advantages of solving small sample size, nonlinear and high dimensional pattern recognition problems; however, its identifying ability is influenced by self-parameters easily.

Chaotic oscillator detection method has some advantages, such as it can effectively extract weak signals under the noise background, it is sensitive to small periodic signals, and it has certain immunity to noise. Therefore, the application areas of the method are very wide., for example, arcing faults protection [10], fault diagnosis of wind turbine gearbox [11], faults line detection [12], and islanding detection [13].

This paper proposed a novel method of fault line selection based on IOS. Firstly, in the aspect of extracting feature signal, TZSC is decomposed by db10 wavelet packet. Next, according to the maximum energy principle, CFB is selected and is inputted to IOS. Finally, fault line selection criterion based on maximum chaotic distance and average chaotic distance of phase trajectory is presented to judge the fault line.

The remaining of this paper is arranged as follows. In Section 2, we analyzed the small current to ground fault. In Section 3, we indicated duffing chaotic oscillator disadvantage that cannot detect the damped oscillation signal. In Section 4, we introduced the IOS work principle. In Section 5, we proposed the fault line selection method and introduced the wavelet packet theory and detection factor setting. In Section 6, we simulated the small current to ground fault and verified the effectiveness of proposed method. In Section 7, the nondetection zone of proposed method is presented. A brief summary concludes the paper in Section 8.

#### 2. Analysis of Small Current to Ground Fault

When SPG fault occurs in small current to ground system, the current flow distribution is shown in Figure 1.

In Figure 1, the power’s neutral point attached an arc suppression coil, so current flow has a change by the coil, but capacitor voltage amplitude and distribution are both the same as nonsolidly grounded system. However, the difference is that inductive current flow from line to ground point. Therefore, the total current is where and were capacitance current to ground and inductance current, respectively. Using to represent inductance, so .

Phase difference between and is 180°, so will decrease due to arc suppression coil compensation. When the coil runs with over compensation stage, the residual current to ground is inductive. Reactive power flow direction of capacitive is from bus line to fault line, it is same as healthy line. Therefore, the fault information is less.

Transient zero-sequence equivalent circuit of SPG fault is used to analyze TZSC. In Figure 2, , are zero-sequence capacitance and zero-sequence inductance, respectively; is transition resistance; , are equivalent resistance and inductance of arc suppression coil respectively; is zero-sequence voltage.

In Figure 2, TZSC is [14] where , are inductive current and capacitive current of TZSC and their initial values are and , respectively (, ), is phase voltage amplitude, is angular frequency of power frequency, , are oscillation angular frequency and attenuation coefficient of TZSC, is decay time constant of inductive current, and is initial phase angel.

When SPG fault occurs, ground point will flow transient capacitive current and transient inductive current; the former has fast attenuation speed and the later has slow speed, and determined the amplitude and frequency of . With power lines length increasing, the self-oscillation frequency decreases, and the free oscillation component amplitude of also decreases; meanwhile, free oscillation duration will be reduced to half power frequency cycle. Moreover, the distribution network is a small current to ground system, so fault signal is weak, and it influenced by noise easily. Therefore, TZSC difference between fault line and healthy line is small, and the fault information is less; it was difficult to achieve fault line selection for distribution networks.

In recent years, duffing chaotic oscillator detection method is used widely, because it is sensitive to small periodic signal and has immunity to noise. Therefore, it can be used to detect weak signals in strong noisy background for power distribution networks.

#### 3. Duffing Chaotic Oscillator Disadvantage

The Holmes-type duffing chaotic oscillator has two characteristics.(1)System motion state is not influenced by large amplitude of irregular noise.(2)System motion state changed by the periodic signal, its frequency is the same as the system, even if the signal is very weak.

Holmes-type duffing equation is given as follows: where , , and are output signal, damping ratio, and built-in signal, respectively, is amplitude factor of built-in signal, and is equal to generally.

If is equal to 0.8 and 0.82625, respectively, we can get chaotic state and periodic state; it is shown in Figure 3.

**(a) Chaotic state**

**(b) Periodic state**

In Figure 3, the system has different states by setting different values. Therefore, duffing chaotic oscillator detecting method usually used system’s different states to detect different signals.

Considering the external signal variability, (3) is changed by scale variation. Then, add measured signal, it is shown as follows: Here, is angular frequency of measured signal, is external signal.

We choose damped oscillation signal to study, it contains a range of frequency components; it was shown in Figure 4; Figure 5 is obtained by calculating (4).

**(a) Original amplitude**

**(b) Amplitude reduced 10 times**

**(c) Amplitude reduced 100 times**

**(d) Amplitude reduced 1000 times**

It showed that duffing oscillator detection system is not influenced by mutative amplitude of external signal in Figure 5. The reason is that external signal has different frequencies, but the system has only one frequency. Therefore, the damped oscillation signal cannot make duffing chaotic oscillator turn to other states, in other words, the duffing oscillator detection system has weak detection ability for damped oscillation signal.

#### 4. IOS Work Principle

TZSC signal of small current to ground system is shown in Figure 6.

In Figure 6, TZSC signal has oscillation and attenuation states, and (2) shows that different TZSC have different oscillation and attenuation states. Therefore, TZSC math model is proposed. Consider the following: Equation (5) consists of four components, where , , , and are component amplitude, natural number, decimal number, attenuation coefficient, and initial phase, respectively, , and . Each component, respectively, represents fundamental component, integer harmonic, noninteger harmonics, and attenuation component.

In order to make up for the disadvantage that cannot detect damped oscillation signal, the built-in signal is replaced by (5); therefore, IOS is proposed and shown as follows: where ), , , , , , is amplitude factor of built-in signal, and is detection factor.

##### 4.1. Amplitude Factor Impact

If , , and is equal to 10^{−4}, 0.025, 0.03, and 0.035, respectively, the IOS phase trajectories are shown in Figure 7.

(a) |

(b) |

(c) |

(d) |

Comparing Figure 3 with Figure 7, in this paper, the state of Figure 7(a) is named chaotic-like state (CLS). Similarly, Figures 7(b) and 7(c) are named critical state (CS), Figure 7(d) is named periodic-like state (PLS).

In Figure 7(a), when IOS is in CLS, points movement are small on phase trajectory, motion tendency is obscure, and it is integrally in chaotic state. When increased to a certain value, the phase trajectory turned to PLS; its movement scale is big and movement tendency is clear; the overall movement state can be thought as circular movement with reduced radius; it is shown in Figure 7(d).

CS is between CLS and PLS, and the overall movement state can be thought as a curve segment keeping away from the origin point; it shows in Figures 7(b) and 7(c).

##### 4.2. Damping Ratio Influence

We used the PLS parameters and set damping ratio as equal to 0.01, 0.05, 0.5, and 1, respectively. By solving IOS, the phase trajectories are shown in Figure 8.

(a) |

(b) |

(c) |

(d) |

According to Figure 8, with damping ratio increasing, movement scale of phase trajectory becomes smaller; overall movement state changed from circular movement with decreasing radius to a clear curve segment. However, when decreased, movement scale became larger and overall movement state is from a clear curve segment turn to a closed curve. In other words, when increased from small to large, PLS phase trajectory followed: .

All in all, rational can make PLS phase trajectory become dense and make the difference between CS and CLS be more obvious; the phenomenon is more beneficial to identifying PLS. But, lower and bigger would make PLS turn to other states. Therefore, choosing the rational damping ratio is very important.

##### 4.3. IOS Example

To demonstrate IOS detection ability, transient zero-sequence signal from paper [15] is chosen and inputted to IOS (, , ). By changing amplitudes and frequencies of transient zero-sequence signal, Figure 9 is obtained.

**(a) Change frequency**

**(b) Change amplitude 1**

**(c) Change amplitude 2**

**(d) No change**

Figure 9(a) results are obtained by the changed frequency signal that comes from paper [15], on the contrary, Figures 9(b) and 9(c) results are obtained by changing the signal amplitude. Similarly, Figure 9(d) results are obtained by importing the signal directly without change.

According to Figure 9, the example signal in the paper can make IOS have different phase trajectories. Therefore, different oscillation damping signals can generate different phase trajectories; it shows that IOS has superior ability to detect oscillation damping signals.

#### 5. Fault Line Selection Method

##### 5.1. Wavelet Packet Transform Theory

Wavelet packet transform (WPT) is an improved method of wavelet transform (WT), and it has better time-frequency resolution, therefore, it can effectively extract fault feature of TZSC signal. WPT is described as follows.

Defining is the closure space of , so is closure space of , and meet the following two-scale equation: where , and have orthogonal relationship; when , (7) can be expressed. Consider the following:

When is equal to , it can generalize , , , orthogonal wavelet packet is determined by basis function: . Since is determined by uniquely, so is called the orthogonal wavelet packet about sequence . About wavelet packet decomposition, it is the band-pass or low-pass filter actually. The filter bandwidth is , where is level number of wavelet decompositions, is the* l*th contacts of wavelet decomposition, and is input signal frequency; WPT decomposition structure is shown in Figure 10.

*A* is low-frequency,* D* is high-frequency. The decomposition has the following equation.

According to power signal feature, WPT can choose appropriate CFB to match signal spectrum adaptively; it improved time-frequency resolution and had a wide prospect. Because daubechies series wavelet has some merits, such as orthogonal characteristic, compact support, and vanishing moments; therefore, db10 wavelet packet is chosen to analyze TZSC signal.

##### 5.2. Detection Factor Setting

Each line’s TZSC is decomposed by db10 wavelet packet; then, the CFB is chosen by maximum energy principle, with the change of grounding resistance, the different CFB are obtained; it is shown in Figure 11.

**(a) Small faults resistance**

**(b) Big faults resistance**

In Figure 11, CFB is also damped oscillation signal, and its maximum amplitude decreased with fault resistance increasing. In other words, when fault situation changed, CFB amplitude will be changed. In order to make IOS have a better applicability, the detection factor should be chosen appropriately. Therefore, CFB was input to IOS, and the equation is proposed as follows: where is CFB signal, , , and . The value ofis chosen by a large number of experiments [12].

##### 5.3. Fault Line Selection Criterion Based on Chaos Distance

Compared with paper [12], it is difficult to accurately distinguish between CS and PLS only by the phase trajectory; the reason is that their phase trajectories are very similar. Therefore, chaotic distance and average chaotic distance are put forward to distinguish the phase trajectory state.

Chaotic distance is square distance from each point to origin point on phase trajectory; it can represent the size of each point movement scale: is the vector which is acquired by solving (4), is the first-order differential vector of , and the maximum of is named .

Average chaotic distance is the average of ; it can represent the overall movement scale: is the number of points on phase trajectory.

By using the chaos distance, a range of CFB signals is chosen and inputted into IOS, by solving IOS, we get Figure 12.

**(a) Chaotic-like state**

**(b) Chaotic-like distance**

**(c) Criticality state**

**(d) Criticality distance**

**(e) Period-like state**

**(f) Period-like distance**

In Figures 12(a) and 12(b), CLS scale is very small, and its points always move near the origin point, so the of CLS is also very small. The of CS increases with time increasing, which complies with curvilinear movement that is away from origin point; it is shown in Figures 12(c) and 12(d). But of PLS shows oscillation damping state, which complies circular movement with reducing radius; it is shown in Figures 12(e) and 12(f). Chaotic distance can clearly indicate movement scale and tendency of each state; it avoids error distinction. and values of Figure 12 are shown in Table 1.

value of CLS is very small, and its magnitude order stabilized near (in this paper, ), but of CS and PLS are both to stabilize near (in the paper, ). For CS and PLS, their has a larger difference. In order to distinguish CS and PLS accurately, is proposed (). While (in the paper, ), CS and PLS states cannot be distinguished by solely. Therefore, it proposed :
where is the maximum value of in all phase trajectories,* j* is the number of phase trajectories.

If (in the paper ), the phase trajectory that has is judged PLS, others are judged CS. If , all phase trajectories are judged CS.

In addition, phase trajectory of magnitude order near the is judged CLS.

From above analysis, CLS, CS, and PLS can be accurately distinguished by and . When SPG fault occurs, CFB with detection factor is inputted to IOS, and it used and to select the fault line. Therefore, the novel criterion of fault line selection method is proposed.

*Step 1. *If only one line has (this line is called* L*), and the phase trajectory magnitude of other lines stabilized near . In other words,* L* makes IOS turn to PLS or CS and other lines make IOS turn to CLS. Therefore,* L* is judged as fault line and other lines are judged as healthy line.

*Step 2. *If many lines have , calculated , If (), the line with maximum makes IOS turn to PLS (this line is called* L*), other lines make IOS turn to CLS or CS. Therefore,* L* is judged as fault line and other lines are judged as healthy line.

*Step 3. *If all lines do not comply with Steps 1 and 2, the fault is judged as bus fault.

#### 6. Simulation and Verification

##### 6.1. Simulation Model

In this paper, SPG fault simulated by ATP-EMTP, and the simulation model is shown in Figure 13. There are 4 lines, such as:* S*_{1},* S*_{2},* S*_{3,} and* S*_{4}; their lengths are 13.5 km, 24 km, 17 km, and 10 km, respectively. Simulation parameters are

*Overhead Line*. Positive-sequence parameters:* R*_{1} = 0.17 Ω/km,* L*_{1} = 1.2 mH/km,* C*_{1} = 9.697 nF/km, zero-sequence parameters:* R*_{0} = 0.23 Ω/km,* L*_{0} = 5.48 mH/km,* C*_{0} = 6 nF/km.

*Cable Line*. Positive-sequence parameters:* R*_{11} = 0.193 Ω/km,* L*_{11} = 0.442 mH/km,* C*_{11} = 143 nF/km, zero-sequence parameters:* R*_{00} = 1.93 Ω/km,* L*_{00} = 5.48 mH/km,* C*_{00} = 143 nF/km.

*Transformer*. 110/10.5 kV, leakage impedances of high voltage side: 0.40 +* j*12.20 Ω, leakage impedances of low voltage side: 0.006 +* j*0.183 Ω, excitation current: 0.672 A, magnetizing flux: 202.2 Wb, magnetic resistance: 400 kΩ.

*Load*. All are delta connected, .

*Arc Suppression Coil*. , .

##### 6.2. SPG Fault Analysis

Set fault resistance and initial angle as equal to 0.0001 Ω, and 0°, respectively; fault line voltage of phases A, B, and C was shown in Figure 14, and each line zero-sequence current was shown in Figure 15.

It shows that phase voltage of healthy phases B and C are increasing and their amplitude is equal, but the phase voltage of fault phase A is equal to 0 approximately in Figure 14. This phenomenon meets the boundary condition of A phase to ground fault. Furthermore, all lines have same voltage value, because each line is linked with the same bus line in parallel connection. Therefore, this paper chose zero-sequence current to study the fault line selection.

In Figure 15, when SPG fault occurs, zero-sequence current value of each line is not equal to 0, its current waveforms have damped oscillation state, and the attenuation speed is fast. The gap of zero-sequence current between fault line and healthy line is smaller after 0.03 s, but the gap of TZSC between fault line and healthy line is larger in 0.02-0.03 s. In addition, healthy line TZSC is similar. Therefore, TZSC is chosen to study fault line selection method in this paper.

##### 6.3. Changing Initial Phase Angle and Resistance

When fault occurred in* S*_{1}, fault location is 5 km away from bus line, and the initial angles are 0°, 30°, 60°, and 90° respectively. Therefore, each TZSC of T/2 decomposed by db10 wavelet packet when fault occurred. According to paper [16], IOS parameters were set as , , and .

To show phase trajectory difference, this paper gives 2 examples, one is that initial angle is equal to 30° and fault resistance is equal to 1100 Ω, to express conveniently, this paper uses (30°, 1100 Ω) to represent the fault condition; the other example is (90°, 2000 Ω). The TZSC and CFB of the two examples were shown in Figures 16 and 17, respectively.

**(a) Transient zero-sequence current**

**(b) Characteristic frequency band**

**(a) Transient zero-sequence current**

**(b) Characteristic frequency band**

It shows that TZSC and CFB are both damped oscillation signal in Figures 16 and 17, different fault conditions have different oscillation damping degrees. Furthermore, CFB can describe time-frequency characteristics of transient signal more clearly in some frequency bands; it provides a basis for fault feature information processing [17]. Therefore, CFB is analyzed by IOS for fault line selection. According to new criterion, fault line selection result is shown in Table 2.

Under the fault condition: (30°, 1100 Ω), according to Table 2, magnitude orders of* S*_{3} and* S*_{4} are both ; it judges* S*_{3},* S*_{4} as healthy line. of* S*_{1} and of* S*_{2} are both less than 0.1; it indicated that* S*_{1},* S*_{2} are PLS or CS, combining , ; therefore,* S*_{1} makes phase trajectory turn to PLS, and* S*_{2} makes phase trajectory turn to CS, which is same as Figure 18. Then,* S*_{1} is judged as fault line; other lines are judged as healthy line.

(a) |

(b) |

(c) |

(d) |

Under the fault condition: (90°, 2000 Ω), the of* S*_{1} meets , because is smaller than 0.1 and the* MD*_{C} magnitude orders of other lines are all 10^{−3}, it judges that* S*_{1} makes IOS turn to PLS and other lines make IOS turn to CLS. Therefore, the* S*_{1} is judged as fault line, and other lines are judged as healthy line. The selection result is the same as the simulation setting. It is as shown in Figure 19.

(a) |

(b) |

(c) |

(d) |

From above analysis, fault line selection results are same as simulation settings, which indicate that the method of this paper can adapt to different resistances and different phase angles.

##### 6.4. Different Lines Fault

Refer to paper [18];* S*_{3} is hybrid line and* S*_{4} is pure cable line, therefore, simulate the different lines occur fault, and the results are shown in Table 3. With cable lines taking part in simulation, TZSC decay process becomes shorter [16], but it does not affect fault line selection results. Table 3 shows that the method has better results when different lines fault occurred.

##### 6.5. Fault with Strong Noise

Considering actual fault situation, fault signal usually has noise; therefore, add 0.5 db random Gaussian white noise to TZSC to validate the method’s noise immunity. The waveform was shown in Figure 20, it shows that the TZSC signal is covered by noise completely. Fault line selection results are shown in Table 4.

Under various fault conditions, the proposed method can also achieve better result, and it indicates that the method has good noise immunity. Furthermore, the fault line selection result is better than paper [19] which adds 20 db noise.

##### 6.6. Different Distances Fault

In fact, fault location has different distances from bus line when fault occurs. In Section 6.6, simulate* S*_{1} line with different distance fault, and the fault resistance is 300 Ω, fault distances are 5 km, 30 km, 35 km, 40 km, and 45 km, respectively. Fault line selection results are shown in Table 5.

According to Table 5, the result is same as simulation setting, which indicates that the proposed method adapts to different distance faults.

#### 7. Nondetection Zone of Proposed Method

##### 7.1. Small Resistance to Ground Fault

Set fault resistance between 0 Ω and 80 Ω. According to proposed method, partial results of fault line selection are shown in Table 6.

The proposed method in this paper does not work for small resistance to ground fault in Table 6. In addition, when fault resistance value is bigger than 10 Ω, cable line and cable hybrids line usually make IOS turn to PLS, and overhead line makes IOS turn to CLS, but when fault resistance value is smaller than 10 Ω, overhead line can also make IOS turn to PLS; for example, the fault condition is (0°, 5 Ω). To express result obviously, phase trajectory of (0°, 10 Ω) is given in Figure 21.

(a) |

(b) |

(c) |

(d) |

It shows that* S*_{1} makes IOS turn to CS,* S*_{2} makes IOS turn to CLS,* S*_{3}, and* S*_{4} make IOS turn to PLS in Figure 21; it contradicts the proposed criterion in this paper (Steps 1, 2, and 3). The reason is that fault signal amplitudes are large and damped oscillation states are obvious in this fault situation, it makes all total driving force of IOS too large. There are two conditions, one is that fault line makes IOS directly cross PLS and turn to CLS, the other is that healthy line makes IOS directly cross CLS and turn to PLS or remain CLS. When fault resistance continues to decrease, fault line and healthy line can both make IOS turn to PLS, such as, fault condition: (0°, 5 Ω). In a word, the most essential thing is that the detection factor is determined by experiment, and it is not suitable for small resistance to ground fault.

##### 7.2. High Resistance to Ground Fault at the End of Line

The length of* S*_{1} is 13.5 km, simulate fault location at the end of* S*_{1}, and fault resistance is set to 1.5 kΩ to 2 kΩ. Its result is shown in Table 7.

According to Table 7, when each line CFB was solved by IOS, the healthy line usually makes IOS turn to PLS or CS, and fault line makes IOS turn to CLS, the phase trajectory of (0°, 1900 Ω) is shown in Figure 22.

(a) |

(b) |

(c) |

(d) |

It is indicated that* S*_{1} makes IOS turn to CLS,* S*_{2} makes IOS turn to CS, and* S*_{3} and* S*_{4} make IOS turn to PLS in Figure 22. This result contradicts the proposed criterion in this paper (Steps 1, 2, and 3). This condition occurs due to the fact that fault location is in the end of line, and the fault resistance is larger. It makes voltage and current decrease and the fault current signal become weaker. Therefore, the proposed method is not suitable for this fault condition.

In conclusion, the proposed method is not suitable for small resistance or high resistance to ground fault in this paper; in other words, detection factor has a range, so how to choose better detection factor to adapt extensive fault conditions is to be further studied.

#### 8. Conclusions

SPG fault line selection method based on IOS is proposed in this paper, and the conclusions are as follows.

IOS has a better recognition for TZSC signal, to a certain extent; it makes up for deficiency of the duffing chaotic oscillator which did not detect oscillation damping signal. In addition, TZSC signal decomposed by db10 wavelet packet, it decreased CFB quantities, and improved IOS solving speed.

PLS phase trajectory can be seen as a circular movement with reduced radius; CS phase trajectory can be seen as a curve segment far away from origin point, but CLS phase trajectory is generally seen as chaotic state of movement. In addition, the of CS is similar to PLS, and their value difference usually is about 0.1. But the between CS and PLS has a significant gap which usually is equal to 1.5 for their ratio. Therefore, phase trajectory state is usually judged by , and of CLS is very small. Besides, magnitude order of is to stabilize at 10^{−3} for CLS.

There are some nondetection zones of the proposed method: detection factor is determined by experiment, it can cause some errors, and the method has problem of insufficient sensitivity when small and high resistance to ground occur fault.

#### Notations

SPG: | Single phase to ground |

TZSC: | Transient zero-sequence current |

CFB: | Characteristic frequency band |

SVM: | Support vector machine |

IOS: | Improved oscillator system |

WT: | Wavelet transform |

WPT: | Wavelet packet transform |

CLS: | Chaotic-like state |

CS: | Critical state |

PLS: | Periodic-like state |

: | Chaotic distance |

: | Maximum of |

: | Average of chaotic distance. |

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this article.

#### Acknowledgments

This work is supported by Science and Technology Research (12B470002, 14A470004) and Control Engineering Lab Project (KG2011-15) of Henan Province, China, and Youth Foundation (Q2012-28, Q2012-43A) of Henan Polytechnic University, China.