Abstract
In order to improve the transient stability of frequency in a small stand-alone microgrid (SSM), this paper takes a SSM composed of a direct-drive permanent magnet synchronous generator (D-PMSG) and a micro gas turbine (MGT) as the background and uses wind turbine generator (WTG) virtual inertia (VI) to participate in the primary (short-term) system frequency regulation. First of all, this paper constructs a grid-connected model composed of a WTG and a MGT, analyzes the WTG virtual inertia frequency regulation mechanism, and explains the principle of proportional-differentiation (PD) virtual inertia control (VIC) and its shortcomings. Secondly, the paper introduces the structure principle of -order active disturbance rejection control (ADRC) and deduces the design process of second-order ADRC-VIC. Finally, through the simulation and experimental verification, comparing the frequency perturbation of without-VIC, PD-VIC, and ADRC-VIC, it is concluded that PD-VIC and ADRC-VIC both can use the WTG virtual inertia to participate in the primary frequency regulation. The frequency regulation effect of ADRC-VIC is better than PD-VIC, ADRC-VIC can extend the rotor speed recovery time and avoid overshoot, and its frequency fluctuation amplitude and settling time are obviously improved, and ADRC-VIC can effectively avoid the overshoot phenomenon of the MGT output power.
1. Introduction
The small stand-alone microgrid (SSM) is an effective way to solve the power supply in the mountains, islands, and remote areas. The principle is to use the collected wind energy, solar energy, and other renewable energy to connect with the traditional synchronous generator (SG), thus building a power supply system independent of the large power grid [1–3]. With the progress of wind power technology, its application in the stand-alone microgrid is more and more extensive. Direct-drive permanent magnet synchronous generator (D-PMSG) has the advantages of flexible control, good reliability, and high power generation efficiency, and it has good application prospect [4]. Micro gas turbine (MGT) has the advantages of low emissions, high efficiency, good fuel adaptability, and flexible control, and it is widely used in distributed power generation technology [5].
In this paper, a SSM constructed by a D-PMSG and MGT is used as the research object to analyze the problem of the transient stability of frequency. The master-control power supply of the SSM is MGT, so the stability of frequency mainly relies on MGT. On the one hand, due to economic reasons, the design capacity of microgrid and the demand of load are matched, resulting in the fact that the microgrid capacity is limited and the system inertia is small. On the other hand, due to technical reasons, the system inertia can only damp the system frequency changes, but it can not completely prevent the system frequency changes [6–8]. Therefore, the system requires the WTG to actively participate in the frequency response, so as to improve the transient stability of system frequency. In general, WTG participates in system frequency regulation in two ways [9–11]. One is called deloaded operation of WTG; when the system frequency has a deviation, the WTG releases the standby power to participate in the system frequency regulation. This process takes a long time, but it is at the expense of wind energy utilization and the operating economy of WTG, and it is difficult to promote the application. The other is virtual inertia control (VIC) [12]; since the WTG itself has mechanical inertia, it has available rotor kinetic energy in operation. The rotor kinetic energy can transform into electromagnetic power through the control algorithm, to provide support for the system frequency, and this process takes a short time and “zero power consumption” [13, 14]. In addition, the permeability of wind power is often larger due to the limited capacity of SSM, which ensures the effectiveness of WTG participating in the system transient frequency regulation [15–17]. Many scholars have started the research on how to use the virtual inertia of WTG in the system frequency regulation and improve the transient stability of frequency.
The paper [12] gives smart microgrid hierarchical frequency control ancillary service provision based on virtual inertia concept, both the static and dynamic frequency securities of an islanded microgrid are provided in primary and secondary control levels. In [10, 11], frequency-based VIC is realized based on the rate of change of frequency (ROCOF) or frequency deviation from the nominal value or a combination of both. The proportional (P) term with makes the WTG respond to frequency changes and the differential (D) term with emulates the inertia. The paper [13] introduces a method of PD-VIC to improve the primary frequency contribution of grid-connected variable speed WTG; by adjusting the droop of the WTG in response to wind velocities, the system primary frequency response is significantly improved. In [14], it is concluded that WTG can provide an extra active power production from its stored kinetic energy; the amount of extra injected power is determined according to the frequency deviation and/or the rate of change of frequency .
The above literatures illustrate the feasibility of using virtual inertia to participate in system frequency regulation and PD-VIC has a certain practicality. However, the PD algorithm has its inherent limitations, such as the fact that it is difficult to determine the optimal control parameters and likely to cause an overshoot, and the robustness of the controller is poor.
This paper attempts to design a controller from the view of suppressing the system frequency disturbance, which has strong robustness and can avoid overshoot. In [18], the performance superiority of WTG with the proposed linear ADRC approach over that with the traditional proportional integral and fuzzy-proportional integral-based controllers is validated by the simulation results. The simulation results indicate that the performance of the LADRC-based controller is much better than PI and Fuzzy-PI-based controllers. The result also shows that the proposed LADRC-based controller has an excellent robustness and a good capability of disturbance rejection. The authors of [19, 20] have presented the ADRC for the doubly fed induction generator (DFIG) used in wind energy conversion systems, the proposed control has demonstrated its effectiveness under wind variation, and both generated power and rotational speed are kept in their safety region.
Based on the above research, this paper studies the PD-VIC and second-order ADRC-VIC algorithms, and it analyzes the frequency regulation mechanism of the two algorithms. Under the premise of not exceeding WTG virtual inertia frequency regulation capability, the characteristics of the two algorithms are compared and analyzed by simulation and experiment.
The remainder of this paper is given as follows. Section 2 gives a single WTG grid-connected structure and the mechanism of VIC. In Section 3, the structure of -order ADRC is described first, and then the second-order ADRC is designed according to the frequency dynamic response equation. Simulation analysis and experiment verification are given in Sections 4 and 5, respectively. Finally, conclusions are drawn in Section 6.
2. A WTG Grid-Connected Structure and the Mechanism of VIC
2.1. A WTG Grid-Connected Structure
After a WTG combined to the grid, the system structure is shown in Figure 1. The WTG includes wind wheel, D-PMSG, full-power converter, and its control system.
The wind turbine (WT) is used to capture the wind energy through the wind wheel to drive the D-PMSG to generate electricity. The motion equation of the WTG rotor iswhere denotes the inertia time constant of D-PMSG, the rotor speed of wind wheel, the mechanical power captured by wind turbine, and the electromagnetic power of the WTG.
The expression of is [21]where denotes the air density, the radius of the wind wheel, the wind speed, the power coefficient of WTG, the pitch angle, and the tip speed ratio. The expression of is as follows:
When the pitch angle is constant, is a function of . can reach the maximum value by adjusting the rotor speed and reaches the optimal value at the same time, as shown in Figure 2(a). At this point, it can be seen from (2) that corresponds to only one maximum point () at different wind speeds and rotor speed , and these maximum points are connected to a line, which is the maximum power point tracking (MPPT) curve, as shown in Figure 2(b).
(a) The maximum power coefficient point
(b) The MPPT curve
When the WTG runs according to the MPPT curve, the wind energy obtained by WTG iswhere is the coefficient that makes WTG obtain the maximum wind energy, and its expression is . When the VIC is introduced into the microgrid, as shown by the red dotted line in Figure 1, the expression of iswhere is the output power of the MPPT control and is the auxiliary frequency regulation power of VIC output. As shown in Figure 1, the input of VIC is the system frequency deviation , and , is the system frequency obtained by the grid-side phase-locked loop (PLL), and is the reference frequency 50 Hz. is calculated by the frequency deviation . When the system is in steady-state operation, that is, , , the WTG runs in the state of MPPT control, and virtual inertia does not participate in the system frequency regulation. When the system frequency is disturbed, that is, , , the WTG adjust the generator rotating speed and output power through the full-power converter to use virtual inertia in the system frequency regulation.
2.2. The Response Mechanism of VIC
The operation process of WTG includes MPPT and VIC. Figure 3 shows the characteristic curves of WTG in the processes of MPPT and VIC. The wind speed keeps constant 7 m·s−1, and the system load is increased by 15% after WTG reached the MPPT point (B). and denote the mechanical power and the electromagnetic power, respectively.
The Process of MPPT (). The process by which WTG changes from point () to point is the MPPT process under the constant wind speed. In this process, changes along the curve of MPPT, and of WTG changes along its characteristic curve. As is greater than , is constantly rising, and the system is stable at point ultimately. As can be seen from the curve of in Figure 3, point is the maximum power point of WTG at this wind speed.
The Process of VIC (). When the system frequency falls due to the sudden increase in load, the VIC of WTG makes climb from points to quickly. In this process, due to the inertia of WTG, does not change immediately, so remains at point . At this time, since is less than , starts to decrease. During the decrease process of , changes along the curve , and runs along the curve , when the system finally runs to point , and then reaches the minimum value. At this time, is greater than , starts to increase, of WTG change along the curve , run along the curve , the system is back to point finally, the mechanical power and the electromagnetic power are balanced again, and the rotor speed is stabilized at the optimal speed of point .
In the process of VIC, the principle of auxiliary frequency regulation is mainly to change the power of WTG to reduce the impact of the system load disturbance on the MGT, thereby reducing the changing rate and amplitude of the system frequency.
2.3. The Rotor Speed Recovery and Overshoot
As shown in Figure 4, when the system frequency drops, VIC can quickly reduce the rotor speed and release the rotor kinetic energy. WTG delivers active power to the system for a short period of time to provide support for the system frequency. The inherent drawback of this process is that, at the rapid recovery of the rotor speed, on the one hand, this causes the oscillation and overshoot phenomenon of the WTG output power and the rotor speed (the shadow portion of the time period); on the other hand, during the rotor speed recovery process, the WTG absorbs energy from the grid (the shadow portion of the time period), resulting in a smaller power supply to the microgrid, which is not conducive to frequency recovery stability.
In fact, the rotor speed recovery process is carried out when the frequency is not stable, this process is to reduce the WTG output power and increase the system load, which has played a “superposition effect” in the decline process of the frequency, making the frequency response characteristics deteriorate, and in severe cases can cause the system frequency to fall twice.
In Figure 4, is the time period of VIC. is the time period of the rotor speed recovery and overshoot, which is the time of the rapid decline in frequency regulation capacity, through extending this time period to slow down the speed recovery rate, so as to reduce the rapid decline of frequency regulation capacity and reduce the system frequency deviation.
3. Second-Order ADRC-VIC
ADRC is proposed by Han [22], which is based on the error feedback control theory of traditional PID and the development of modern control theory [23, 24]. ADRC is less dependent on the accuracy of the mathematical model of the controlled object, and the main idea is to estimate the system real-time disturbance by using the extended state observer (ESO) [25] and compensate it. ADRC has the characteristics of fast response, robustness, and adaptability [26]. In this paper, the second-order ADRC is designed according to the frequency dynamic response equation of the microgrid.
MGT is the main control power supply, which can be equivalent to conventional synchronous generator model for analysis, and the system voltage fluctuations caused by the load switching will not be taken into account. Thus, the constant of the synchronous generator model (in p.u.) is as follows:where , , , and are -axis voltage, -axis voltage, -axis current, and -axis current of the generator, respectively. and are the -axis transient reactance and -axis reactance, respectively. denotes the stator resistance. is the rotor speed. and are the rotational inertia and the damping coefficient, respectively. is the mechanical torque.
Laplace transformation of (7) can be expressed as
The speed governor model of the synchronous generator can be expressed aswhere and are the proportional and integral coefficients of PI controller, respectively. is the reference value of the synchronous generator speed. The active power of the microgrid synchronous generator can be expressed as
Substituting (8), (9), and into (8) leads to
It is known from (10) that the VIC system frequency response dynamic equation is
Substituting (10) into (11) leads towhere .
The relationship between the system frequency (p.u.) and the generator rotor speed (p.u.) is expressed as
Substituting (13) into (12) leads to
Both sides of the equation are multiplied by , and then
Namely,
Substituting into (16) leads to
The second-order differential equation of the power system frequency deviation is expressed aswhere the uncertain items of the controlled object can be defined as the sum of internal and external disturbance as follows:
The differential equation of and in (18) indicates that when the ADRC is substituted for the conventional PD controller, the output of the controlled object is the same as that of the conventional PD controller, which is the value of the system frequency deviation . The block diagram of the second-order ADRC is based on (18), as shown in Figure 5.
In Figure 5, is the estimated value of the system frequency deviation value , is the estimated value of , and is the estimated value of the system uncertain item ). denotes the compensation coefficient of the estimated values and . In this paper, and . Since the purpose of the ADRC is to suppress fluctuations in system frequency, the reference values and of variation in status can be directly set to 0, respectively. Thus, the design of TD shown in Figure 5 is left out.
In practical engineering applications, there is a need to discretize the algorithm of each module for the ADRC [22, 23]. The discretized form of the nonlinear ESO can be expressed aswhere is the sampling step size and , and are the observer parameters needed to be determined. According to the characteristics of inherent parameters of the ADRC, , , and , which are used in this paper. In order to avoid the chattering, the saturation function is expressed aswhere is the linear interval length and sgn(·) is the symbolic function and denotes the saturation function factor. When the absolute value of the error is relatively large, its gain will decrease to prevent the overshoot. When the absolute value of the error is relatively small, its gain will increase to speed up the convergence process.
The NLSEF and the discrete control amount produced by disturbance compensation are expressed aswhere and are the output correction coefficients and the parameters are used to mediate the nonlinearity of the function. The ADRC designed in this paper has the following advantages over the conventional PD controller: due to its observation and compensation of power disturbance, its better dynamic characteristic of suppressing the frequency disturbance; increasing the error tracking efficiency by designing the optimized NLSEF and using saturation function to achieve the objective of “small error and large gain, large error and small gain”; its better applicability and robustness because of adopting ESO for the estimation and compensation of the system parameter disturbance.
Although the parameters of ADRC are much more than PD, the advantages of ADRC over PD are as follows. On the one hand, in order to release more WTG virtual inertia, the PD is achieved by increasing the proportional factor and differential factor , but this causes the WTG speed overshoot phenomenon, when the overshoot is serious which will cause the system frequency secondary drop. On the other hand, as more WTG virtual inertia was released，ADRC will not cause the WTG rotor speed overshoot, and its ability to adjust the system frequency will be further improved. And ADRC has a certain ability to resist parameter disturbance. The following is through simulation and experiment to illustrate this.
4. Simulation Analysis
In this paper, a SSM simulation model composed of a WGT and a MGT is built in MATLAB/Simulink. The main parameters of the simulation model are shown in Appendix A. First of all, the estimation effect of ESO and the robustness of antiparameter perturbation of ADRC are verified by simulation.
4.1. The Simulation of ESO Observation Effect
Figure 6 shows the estimated effect of ESO when the load is switched. The red solid line is the p.u. value of the actual system power perturbation, . The blue dotted line is the ESO estimate value. is the estimate amount of . It can be seen from the figure that ESO can realize the accurate tracking and estimation of the system load disturbance, which provides a reliable basis for the compensation and suppression of the system frequency disturbance.
(a) Sudden increase in load
(b) Sudden decrease in load
During the operation of the system, the inertia time constant of the system is in the process of change due to the random fluctuation of the wind resource and the load disturbance. Therefore, it is necessary to examine the impact of parameter on control algorithm performance. When the wind speed fluctuates randomly and takes 1/12, 1/6, or 1/24, respectively, the ADRC control characteristic is shown in Figure 7.
(a) The effect of ESO
(b) The frequency response curve
Figures 7(a) and 7(b) are the estimated effect curves and system frequency response curves of ESO, respectively. As can be seen from Figure 7(b), the frequency response curves of the system are completely coincident with different parameters . The ESO of ADRC has strong robustness, which is inherent in ESO performance. When the parameter is disturbed within a certain range, the ADRC control performance is not affected, which is determined by its robustness. Therefore, the change of the value of will not affect the ESO accurate estimation of the system power disturbance; that is, it will not affect the control effect of ADRC on the system frequency. This shows that ADRC has a certain ability to resist system parameters perturbation, which reduces the ADRC dependence on grid parameters.
4.2. The Simulation of PD-VIC and Second-Order ADRC-VIC
In order to compare the effects of the PD-VIC and second-order ADRC-VIC to suppress the frequency disturbance, the microgrid frequency disturbance is caused by sudden increase load. Figure 8 is the waveforms of simulation during the sudden increase in load, and wind speed is 7 m/s.
(a) Microgrid frequency
(b) WTG rotor speed
(c) WTG output power
(d) MGT output power
The design of simulation scenario: at s, a WTG and MGT are connected to bear the fixed load kW together, the output power of WTG is kW, the rotor speed of WTG is rad/s, and the output power of MGT is kW. At s, the sudden increase in load is kW.
In Figure 8, blue dotted lines are without-VIC simulation waveforms. At s, is put into the system. MGT independently bears the load disturbance, and WTG does not participate in frequency regulation. The frequency nadir is 49.55 Hz and the setting time is at s after two oscillation periods. kW and rad/s remain unchanged. Since the load disturbance is only borne by MGT, increases 12 kW by step.
The red long-dotted lines are PD-VIC simulation waveforms. At s, the WTG uses virtual inertia to participate in frequency regulation. The lowest point drops to is 49.7 Hz and the setting time is at s after one oscillation period. During the sudden increase in load, increases to the maximum value 4 kW and then returns to 3 kW, falls to 13.25 rad/s and then returns to 14 rad/s, and and have the phenomenon of oscillation and overshoot. increases to 12 kW after a certain distortion and appears overshoot phenomenon which is higher than 12 kW.
The black solid lines are the second-order ADRC-VIC simulation waveforms. The lowest point drops to is 49.85 Hz and the setting time is at s. During the sudden increase in load, increases to a maximum point 4.5 kW and then returns to 3 kW, falls to 13 rad/s and then returns to 14rad/s, and and have no phenomenon of oscillation and overshoot. The recovery time of is longer than PD-VIC. Compared with that in PD-VIC, does not appear overshoot phenomenon, can be more smooth added to 12 kW, and second-order ADRC-VIC reduces the impact of load disturbance on MGT.
Based on the above conclusions, when the wind speed is constant, PD-VIC and second-order ADRC-VIC both can use the WTG virtual inertia to participate in the shot-term frequency regulation and effectively reduce the transient fluctuations of the system frequency. Compared with PD-VIC, second-order ADRC-VIC can extend the recovery time of and avoid overshoot, reduce the amplitude of the frequency by 0.15 Hz and the settling time by 2 s, and make more smooth transition during load disturbances.
When the wind speed fluctuates randomly and the system load is constant, the system frequency disturbance is mainly caused by the fluctuation of WTG output power. The effects of different VIC algorithms on suppressing frequency perturbation are compared by simulation in Figure 9. Compared with without-VIC and PD-VIC, second-order ADRC-VIC (blue dash-dotted line) frequency disturbance is significantly smaller. As WTG uses second-order ADRC-VIC to increase the intensity of participating in the frequency regulation, the disturbance range of the rotor speed becomes larger, which indicates that WTG virtual inertia can be more fully released.
(a) Wind speed
(b) Microgrid frequency
(c) WTG output power
(d) WTG rotor speed
5. Experimental Verification
5.1. The Design of Simulative Experiment Platform
In order to further verify the effectiveness and feasibility of the above algorithms, a set of simulative WTG grid-connected experimental systems with capacity of 27.5 kVA is built, which includes a simulative WTG and a virtual synchronous generator (VSG) simulated by a full-power converter. The full-power converter simulates the MGT frequency characteristics according to the algorithm in [27]. Figures 10 and 11 are the schematic diagram of the simulative experiment system and the field of the simulative experimental platform, respectively.
In the experimental system, the rated power of the simulative WTG is kW, the rated rotor speed is rad/s, the radius of the wind turbine is m, and the rated wind speed is 10 m/s. The rated capacity of simulative MGT is kW and inertia time constant is s. The fixed load is kW, and the switchable load is kW. The controller chips of simulative WTG and MGT are DSP F2812 of TI company, the controllers exchange data with LabWindows/CVI of host computers through the RS485 communication protocol.
The design of experimental scenario: at , the WTG and MGT are grid-connected, kW, and kW. At s, the system load was increased by 5 kW. Three groups of experiment are carried out in the order of without-VIC, PD-VIC, and ADRC-VIC. Since the focus of this paper is on the VIC algorithm of the short-time frequency regulation, the wind speed is set to a constant 6.5 m/s during this period.
5.2. The Analysis of Experimental Results
Figure 12 is the experimental waveforms of without-VIC, PD-VIC, and second-order ADRC-VIC during the load increase, respectively.
(a) Microgrid frequency
(b) WTG rotor speed
(c) WTG output power
(d) MGT output power
The analysis of without-VIC experimental waveforms: at s, with sudden increase of the load , the simulative MGT independently bears the load disturbance, and the fluctuation characteristic of frequency is similar to the actual MGT. As shown in Figure 12(a), the maximum frequency excursion is about Hz, the stabilization time of frequency is about s after three cycles of oscillation, and s. In Figures 12(b) and 12(c), the output power of WTG kW and the rotor speed rad/s remain unchanged. In Figure 12(d), the increases 5 kW by step.
The analysis of PD-VIC experimental waveforms: compared with that in without-VIC, the fluctuation amplitude of is decreased, Hz, the settling time is about s after two cycles of oscillation, and . In Figure 12(c), the increased auxiliary power of is the use of virtual inertia to participate in frequency regulation, the extreme value is Δ kW, and and have the phenomenon of oscillation and overshoot. is increased 5 kW with the load increasing, but the WTG participates in the system frequency regulation with virtual inertia at this time, resulting in a distortion and overshoot phenomenon of .
The analysis of second-order ADRC-VIC experimental waveforms: compared with that in without-VIC, the fluctuation amplitude of frequency is obviously reduced, the maximum deviation of the frequency is about Hz, the settling time is at s or so, and s. The extreme of auxiliary power is kW. and have no phenomenon of oscillation and overshoot. Second-order ADRC-VIC makes recovery time longer and avoids overshoot, so that the WTG virtual inertia is released more fully. appears to be no overshoot phenomenon.
Therefore, the PD-VIC and the second-order ADRC-VIC both can effectively improve the transient stability of the system frequency in the same experimental environment. Compared with the PD-VIC, the ADRC-VIC can effectively extend the recovery time and avoid overshoot, and second-order ADRC-VIC reduces the microgrid frequency fluctuation amplitude value by 0.2 Hz, the time of stabilization is reduced by about 1.5 s, and the overshoot phenomenon in can be avoided.
6. Conclusions and Future Works
6.1. The Conclusions of This Paper
Based on the background of a SSM, this paper studies the algorithm of primary frequency regulation by using WTG virtual inertia and focuses on the algorithms of PD-VIC and second-order ADRC-VIC. The paper shows that the PD-VIC is clear in physical meaning and easy to implement, but there is a problem that the setting of the optimal parameter is difficult and it is easy to cause overshoot. ADRC is a controller designed to suppress the system frequency disturbance. This paper presents the design and implementation of the second-order ADRC-VIC.
Through the verification of simulation and experiment, PD-VIC and second-order ADRC-VIC both can effectively use the WTG virtual inertia to participate in the system primary frequency regulation. Compared with PD-VIC, the second-order ADRC-VIC can effectively extend the recovery time and avoid overshoot, and the simulation shows that second-order ADRC-VIC reduces the maximum deviation and the settling time of frequency by 0.15 Hz and 2 s, respectively, the experiment shows that ADRC reduces the maximum deviation and the settling time of frequency by 0.2 Hz and 1.5 s, respectively, and ADRC can avoid the overshoot of the MGT output power. Therefore, ADRC makes WTG virtual inertia release more full, and the frequency regulation effect of second-order ADRC-VIC is better.
6.2. The Future Works
This paper only studies the algorithm of the WTG virtual inertial control but does not fully consider the WTG frequency regulation capability. How to fully release the WTG virtual inertia to participate in frequency regulation is a research that needs to be done in the future. In this paper, a single WTG participating in the frequency regulation is studied only, and multiple WTGs in coordinate frequency regulation are the next research to be carried out.
Appendix
A. Parameters of Simulation Model
MGT rate power : 50 kW MGT inertia time constant : 8 s MGT line-to-line voltage : 380 V a WTG rate power : 20 kW Microgrid steady frequency : 50 Hz Fixed load : 10 kW Switchable load : 5 kW PD coefficients , , Sampling step size : 0.001 s ESO parameters , , and : 1000, 334, and 332 ESO parameters , , , and : 0.5, 0.8, 1, and 1/24 NLSFE parameters , : 70, 10 NLSFE parameters , , , and : 0.5, 2, 1, and 24
B. Parameters of Experimental System
MGT rate power : 20 kW MGT inertia time constant : 10 s MGT line-to-line voltage : 380 V a WTG rate power : 7.5 kW Fixed load : 10 kW Switchable load : 5 kW PD coefficients , , Sampling step size : 0.01 s ESO parameters , , and : 100, 34, and 31 ESO parameters , , , and : 0.7, 0.9, 1, and 1/40 NLSFE parameters , : 90, 30 NLSFE parameters , , , and : 0.7, 5, 3, and 40
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant no. 51377167) and the National Basic Research Program (973 Program) of China (2012GB215103).