## Abstract

It has been observed that changes in the birefringence, which are difficult or impossible to directly measure, can significantly affect mode-locking in a fiber laser. In this work we develop techniques to estimate the effective birefringence by comparing a test measurement of a given objective function against a learned library. In particular, a toroidal search algorithm is applied to the laser cavity for various birefringence values by varying the waveplate and polarizer angles at incommensurate angular frequencies, thus producing a time-series of the objective function. The resulting time series, which is converted to a spectrogram and then dimensionally reduced with a singular value decomposition, is then labelled with the corresponding effective birefringence and concatenated into a library of modes. A sparse search algorithm (*L*_{1}-norm optimization) is then applied to a test measurement in order to classify the birefringence of the fiber laser. Simulations show that the sparse search algorithm performs very well in recognizing cavity birefringence even in the presence of noise and/or noisy measurements. Once classified, the wave plates and polarizers can be adjusted using servo-control motors to the optimal positions obtained from the toroidal search. The result is an efficient, self-tuning laser.

© 2014 Optical Society of America

## 1. Introduction

Mode-locked fiber lasers have continued to make tremendous strides in engineering performance, both in terms of output energy and peak powers, yielding approximately a two orders of magnitude power increase in the last decade [1]. Despite such exceptional experimental achievements, first-principals modeling of fiber lasers has remained, for over two decades now, *qualitative* in nature. The underlying and primary reason which has prevented *quantitative* modeling efforts is the fiber birefringence [2–5] (See also the recent review article by Gordon and Kogelnik [6]). It is well understood from these studies, primarily aimed at fiber optic communications, that the fiber birefringence is stochastic in nature, varying randomly along the length of the fiber laser cavity and highly susceptible (and sensitive) to environmental factors such as bend, twist, anisotropic stress, and ambient conditions such as temperature. With modern data-analysis methods, we propose to optimize mode-locking performance by learning a proxy measure and classification for the fiber birefringence, thus allowing for a self-tuning laser design capable of adapting rapidly to changes in birefringence.

Just as in optical communications, a fiber laser cavity propagates pulses over ultra-long distances in fractions of a second. Signal distortions due to the chromatic dispersion and nonlinearity accumulate after many round trips of the laser cavity, as does the signal distortion due to the fiber birefringence [2–6]. Successful pulsed laser operation is achieved when the linear and nonlinear cavity effects balance each other resulting in stable mode-locked pulses [7, 8]. Although single-mode fibers are typically used for such laser cavities, the so-called single-mode fibers in fact support two modes simultaneously, which are orthogonally polarized. In an idealized circular-core fiber, these two modes will propagate with the same phase velocity. However, practical fibers are not perfectly circularly symmetric. As a result, the two modes propagate with slightly different phase and group velocities due to small differences in the effective index of refraction experienced by each. While this birefringence is small in absolute terms in standard optical fibers, approximately 10^{−7} index of refraction difference in the two modes, the corresponding beat length *L _{B}* is about 10 meters with variations occurring on lengths of 100 meters, which is often on the same order as the dispersive and/or nonlinear length scales. As a result, the birefringence can have a significant impact on mode-locking dynamics.

To illustrate the cavity sensitivity to birefringence, consider one of the most commercially successful mode-locked lasers to date (See Fig. 1): the well-known mode-locked fiber laser that relies on nonlinear polarization rotation (NPR) for achieving saturable absorption using a combination of waveplates and polarizer [8–11]. This NPR based laser concept is more than two decades old and is so successful in part due to its reliance on simple off-the-shelf telecom components, rendering it a highly cost-effective mode-locking source. More recently, tremendous performance advances in this NPR laser have been made in power delivery by using all-normal dispersion fiber cavities [12–14] and/or self-similar pulse evolutions [15–17]. It has also been recently conjectured that multi-NPR sections can be used in the cavity to overcome the multi-pulsing instability and achieve additional performance gains [18–20]. However, such commercial lasers must enforce strict environmental control to maintain performance, i.e., the fiber birefringence is controlled by pinning into place and shielding it from temperature fluctuations. Such system sensitivity has prevented it from major performance advances, limiting power and pulsewidths. Moreover, failure to accurately model the stochastic and sensitive birefringence fluctuations in the cavity have deprived the community for more than two decades of a quantitatively accurate model of this highly successful laser system.

Our objective in this manuscript is to make use of modern data-analysis methods, i.e. machine learning techniques, to help discover a proxy measure for the effective cavity birefringence. Unlike optical communication lines where over long distances a statistical average might be experienced by a pulse, here a single realization of a stochastic variation of the birefringence is what drives the laser cavity dynamics (See Fig. 1(c)). If the cavity is perturbed by bend, twist, anisotropic stress, and/or ambient temperature, then a new realization results. For optimizing performance, it is critical to characterize, or recognize, the fiber birefringence correctly in order to determine the waveplate and polarizer settings, for instance, required to give the best energy performance.

Using pattern learning methods, we demonstrate that spectrogram measurements, which are dimensionally reduced using a singular value decomposition, uniquely characterize the average cavity birefringence. This gives rise to a sparse representation and classification scheme for identifying the dynamic regime of the cavity. Our algorithm allows for efficient self-tuning of the laser cavity when combined with an adaptive controller [21] with servo-driven components [22, 23]. Thus instead of attempting to model the stochastic birefringence fluctuations directly, which can only practically be done in a statistical way, we instead measure and learn the impact of birefringence on mode-locking performance and provide a method by which optimal, self-tuning can be achieved. Although we demonstrate the method on a computational model, the algorithm would ideally apply directly to the laser cavity as there is no need for an underlying model.

The paper is outlined as follows: Sec. 2 outlines the theoretical (computational) model used for the laser cavity along with the adaptive controller and objective function for optimizing the cavity performance. Section 3 develops a toroidal search algorithm which allows for an extensive and rapid search of the best candidate mode-locked states in the system. The construction of spectrograms for various average birefringence values are also shown. Section 4 develops the primary contribution of the paper: an algorithm that is capable of recognizing the average birefringence state of the system. Key to an accurate recognition is the application of a sparse classifier, which is a novel method for identifying dynamical regimes. Section 5 gives the results of the classification methodology and shows the success of the algorithm. An overview of the method and future outlook of its applicability are given in the concluding Sec. 6.

## 2. Fiber laser model and objective function

In previous work [21], an extremum seeking controller (ESC) was designed to achieve and maintain a high-energy, single-pulse state in a mode-locked fiber laser. Given initial parameter values, the ESC enables us to obtain the local maximum of our objective function by varying the fiber laser control parameters, i.e. waveplates and polarizers for a cavity based upon nonlinear polarization rotation. However, the performance of the laser is still limited since the controller itself is only capable of tracking the *local maximum* of the objective function instead of the desired *global maximum*. As a result, we propose a data-driven technique which allows us to build a library of time series of the objective function that globally samples the entire parameter space for different birefringence values. Once a sufficiently large library is obtained, the state of fiber laser system can be then characterized by matching the current system behavior with library entries, hence the global optimal performance can be identified, and then subsequently maintained by ESC. In order to perform such a task, an example model laser system is introduced.

#### 2.1. Laser Cavity Model

In order to demonstrate that the machine learning algorithm provides us with the accurate birefringence characterization, full simulations of a laser cavity are performed. In practice, the algorithm advocated here would sample from the experimental laser dynamics and no theoretical model would be required.

To describe the propagation dynamics in the laser fiber, including the interaction of chromatic dispersion, self-phase modulation, birefringence, cavity attenuation/loss, and bandwidth limited gain and saturation, we use the coupled nonlinear Schrödinger equations (CNLS) [3, 4]:

*u*and

*v*represent two orthogonally polarized electric field envelopes (the fast- and slow-axis respectively) in the fiber cavity of an optical fiber with birefringence

*K*. The variable

*z*denotes the propagation distance which is normalized by the length of the first fiber section and

*t*is the retarded time normalized by the full-width at half-maximum of the pulse. The parameter

*D*is the averaged group velocity dispersion of the fiber section. It is positive for anomalous dispersion and negative for normal dispersion. The nonlinear coupling parameters A (cross-phase modulation) and B (four-wave mixing) are determined by the material of the optical fiber. For axially symmetric fibers

*A*= 2/3 and

*B*= 1/3. The right hand side of the equations, which are dissipative terms, account for the bandwidth limited gain saturation and attenuation, where the operator

*R*of the dissipative terms is defined as follows:

*g*

_{0}and

*e*

_{0}are the nondimensional pumping strength and the saturating energy of the gain. Parameter

*τ*characterizes the bandwidth of the pump, and Γ measures the losses (taken to be distributed) caused by the output coupling and the fiber attenuation.

The waveplates and polarizer are modeled by discrete components in the laser cavity corresponding to Jones matrices. The standard Jones matrices of the quarter-waveplate, half-waveplate and polarizer are given, respectively, by:

*α*(

_{k}*k*= 1, 2, 3,

*p*), the Jones matrices are modified so that where

*W*is one of the given Jones matrices and

*R*is the rotation (alignment) matrix:

The CNLS model together with the Jones matrices provides a full description of pulse propagation and mode-locking dynamics in the laser system. The full simulation of the laser evolves the CNLS equations, using a spectral (Fourier transform) decomposition in the time domain and an adaptive Runge-Kutta method for propagating along the fiber, and a periodic (after every round trip) implementation of the Jones matrices of waveplates and polarizers. The discrete application of Jones matrices after each cavity round trip produces an effective saturable absorption that can be used to control and tune the mode-locking dynamics.

#### 2.2. Objective Function

For any sampling and machine learning method to work, one need to decide how to characterize the state of the laser system effectively. In previous work [21], an objective function *O* was introduced which was obtained by dividing the pulse energy *E* by the spectral kurtosis *M*_{4} (fourth-moment) of the wave form:

*M*

_{4}measures the spread of the waveform.

The objective function enables us to discriminate multi-pulse and chaotic wave forms from the desired single-pulse wave form while providing us with a quantity that favors tight, high-energy, single-pulse wave forms. Figure 2(h) shows the normalized objective function (red), pulse energy (black) and kurtosis of the spectrum (blue) when rotating on the 2-torus of wave-plate *α*_{3} and polarizer *α _{p}*. The objective function selected is ideal for optimizing pulse energy while simultaneously keeping the mode-locking away from instability boundaries (gray regions) [21]. Figure 2(g) shows the wave forms that correspond to the settings marked by the circle, square, diamond, and triangle in Fig. 2(h). The maximal energy occurs away from the regions of single-pulse mode locking (white region), thus tracking energy alone as the objective function would lead to chaotic solutions (grey regions). In contrast, the objective function has local maxima in the single-pulse, mode-locked region. The spectral kurtosis (blue) is small in the single-pulse mode-locking region and is much larger in the grey regions because of multi-pulsing or chaotic wave forms. As a result, dividing energy by the spectral kurtosis penalizes non-mode-locked and multi-pulsing mode-locked states. In the remainder of this paper, we will use this objective function to characterize the state of the fiber laser. It should be noted, however, that any objective function could be picked so that the shortest pulse, broadest spectrum or some other quantity of interest can be selected as a quantity of interest.

## 3. Toroidal search and library building

In this section, a search algorithm will be developed for characterizing the optimal mode-locking performance as a function of birefringence. In practice, this algorithm could take a few to tens of minutes to execute given the fact that mode-locking itself occurs in microseconds. Thus the only limitation to how fast the algorithm can be executed is the speed of the servo-driven controllers [22, 23].

#### 3.1. Toroidal Search

In the experimental setting, the NPR unit consists of the polarizer and three waveplates that can be rotated from 0 to 2*π*, thereby creating a parameter space that is a 4-torus. Multiple NPR sections can be included, let’s say *N* of them, to further enhance performance [18–20]. In order to sample the resulting parameter space (4*N*-torus), a toroidal search algorithm is developed. If we want to sample this 4*N*-dimensional torus, then 4*N* time series are constructed from:

*j*∈ [1,...,4

*N*], where

*θ*

_{j0}are initial parameter values,

*ω*are angular frequencies which are incommensurate, i.e. equation doesn’t have integer solution. In other words,

_{i}*ω*is irrational for any

_{j}/ω_{k}*j*,

*k*∈ [1,...,4

*N*], given

*j*≠

*k*. It is easy to prove that under such conditions, [

*θ*

_{1}(

*t*) ...

*θ*

_{4}

*(*

_{N}*t*)] is dense on the torus [24]. Thus using this method, it is guaranteed that one can sample any torus sufficiently well if sampling for a long enough time or using a high enough sample rate. In the following sections of this paper, we take a 2-torus for a single NPR laser as a simple example case. However, the methods mentioned in this paper can be applied to toroidal parameter spaces of any dimensionality. Figure 2 shows how the toroidal sampling works on a 2-torus comprised of parameters

*α*

_{3}and

*α*. Specifically, the resulting time series of the objective function

_{p}*O*is demonstrated as the 2-torus is sampled. In Fig. 2(a), the torus in under-sampled and aliasing of the objective function occurs. However, as the sampling rate is increased, as shown in Fig. 2(e), the objective function is fully constructed and an evaluation of best performance can be ascertained. Indeed, the red dot in Fig. 2(f) shows the optimal global maxima of the laser cavity. The narrow shaded region around this peak performance is highlighted in Fig. 2(h) where an additional evaluation is made of whether the solution is mode-locked or not.

#### 3.2. The Gábor Transform (spectrogram) and Library Building

For a given birefringence value, toroidal sampling is used to produce a time series of the objective function. Note that we use the entire time series to compute spectrograms in our library building process. Once the library is constructed, optimal parameter settings are kept for future use. In order to develop a robust algorithm that matches the current objective function time series with the library entries, we want to utilize both the temporal and spectral (frequency) signatures of the time series. As a result, we introduce the Gábor transform and construct a spectrogram [25] of the optimal solution.

It is observed that the time series collected from toroidal sampling are comprised of various frequency components that are exhibited at different times. Although the Fourier transform of the signal contains all frequency information, there is no indication of when each frequency occurs in time. Indeed, by definition, the Fourier transform eliminates all time-domain information since we integrate over all time. To circumvent the limitation of the direct application of Fourier transform, Gábor proposed a formal method for keeping information in both time and frequency. His method involved a simple modification of the Fourier transform kernel:

where the filter*g*(

*τ*−

*t*) was introduced with the aim of localizing both time and frequency. The Gábor transform, also known as the short-time Fourier transform is then defined as:

*g*(

*τ*−

*t*) acts as a time filter for localizing the signal and its frequency content over a specific window of time, allowing for the construction of a spectrogram. A spectrogram represents a time series (signal) in both the time and spectral domain, as shown in Fig. 3.

Our key observation is that these spectrograms are unique for varying cavity birefringence. Thus the spectrogram serves as a proxy measure for classifying the underlying cavity birefringence. Unique spectrograms for various birefringence values in the library are shown in Fig 4. By definition the spectrograms are symmetric in frequency, for storage and computation efficiency concerns, we only use the positive frequency part of the spectrograms for classification purposes. These spectrograms serve as the basis of a pattern recognition/classification scheme for determining the value of cavity birefringence.

## 4. Birefringence classification and recognition

In the previous section, we built a spectrogram library using a Gábor transform with Gaussian window. To match a current time series with library entries, the Gábor transform is also applied to the time series sampled from the current system to get the spectrogram. A sparse sampling technique is then applied to recognize the current birefringence from possible library elements.

To start, assume we have computed the spectrogram *S _{k}* for a large number of possible birefringence values where

*k*ranges from 1 to

*M*, and for each

*k*, a singular value decomposition (SVD) is applied to the spectrogram [25]:

*k*value, we keep the first

*m*(

*m*<

*n*) modes (low-rank approximation) of

*U*which has the highest energy and store them in the modes library

_{k}*U*such that where the

_{L}*k*-th sub-library

*Ũ*contains the first

_{k}*m*modes of

*U*: Once we have constructed our dimensionally reduced modes library, we can take a measurement of the laser system (the objective function) and compute the spectrogram. Note that the sampling time does not have to be of the same length as the time series collected when the library was built. We perform an SVD reduction on the measured spectrogram and keep the first

_{k}*m*modes as before, as illustrated in Fig. 5. With the most important (dominant) modes from the measurement in hand, we can do an

*L*

_{1}-norm library search, thus promoting sparsity in our solution [25]. In the

*L*

_{1}-norm search, our objective is to find a vector subject to Here we require the number of library modes to be greater than the dimensionality of the frequency domain. Given this condition, this becomes an underdetermined linear system of equations. The

*L*

_{1}-norm minimization produces a sparse vector

*a*, i.e. only a small portion of the elements are non-zero, as shown in Fig. 6. The non-zero elements of vector

*a*act as a classifier (indicator function) for identifying which sub-library the birefringence falls into. Thus if the largest element falls into the

*i*-th sub-library, the recognized birefringence value is equal to

*K*. This sparsity promoting optimization, when used in conjunction with the unique spectrograms, gives a rapid and accurate classification scheme for the fiber birefringence. Thus birefringence recognition can be easily accomplished. Note that our classification scheme essentially uses the

_{i}*L*

_{1}-norm minimization produce as an indicator function for the correct library elements. More sophisticated sparse classification/recognition strategies can be applied if desired [26], potentially yielding even better recognition results.

## 5. Classification results

The sparse search algorithm is tested using the pre-computed spectrograms with Gaussian random noise added. In our test, the birefringence *K* is varied following a gaussian random walk. For each trial, the spectrogram corresponding to the current birefringence is computed and the *L*_{1}-norm sparse search is executed. Recognition results and errors are showed in Fig 7. In the figure, the recognition algorithm is tested in two scenarios: (i) well-aligned data given the assumption that the servo motors that control the waveplates and polarizers work without error, and (ii) the mis-aligned data that considers the error in the initial angle of the servo motors. In both of these two scenarios, the sparse search works very well. In the well-aligned case, a birefringence recognition (classification) rate of 98% is achieved while in the mis-aligned case, we get a 88% recognition rate. It should also be noted that even when our recognition algorithm fails to find the correct birefringence value, the error between the true birefringence and the recognized value is very small. Thus, even if we use the mis-classified birefringence, is is likely that the predicted optimal parameters will be near the true optimal parameters and the adaptive controller [21] will bring the laser system back to peak performance.

## 6. Conclusions and outlook

The stochastic nature of fiber birefringence has been the major impediment in allowing for *quantitatively* accurate modeling of fiber lasers for optimizing their performance. Indeed, all other physical parameters in the system, such as the Kerr nonlinearity, dispersion characteristics as a function of wavelength, gain and gain bandwidth, can be fairly well characterized in theoretical models. Thus only the birefringence remains unknown and randomly varying. And unlike optical communications, where statistical averaging methods can be used to quantify its effects statistically, a fixed laser cavity represents a single, and unknown, statistical realization of the birefringence which is highly susceptible and sensitive to environmental factors such as bend, twist, anisotropic stress, and ambient conditions such as temperature. Such a system requires new modeling methods which are based upon state-of-the-art data-driven strategies.

We have demonstrated that a toroidal search method forms the basis of a machine learning algorithm for characterizing the performance of a mode-locked laser cavity as a function of birefringence. By constructing a library of objective function spectrograms and their peak performance, a unique signature is given for various cavity birefringence values. Indeed, each unique spectrogram can be reduced using a singular value decomposition and stored in a library for future reference and classification. To evaluate the current cavities birefringence, a time-series sample of the objective function is taken, converted to a spectrogram, reduced using a singular value decomposition, and classified using a sparsity promoting *L*_{1}-norm optimization routine. Accuracies as high as 98% are achieved, thus suggesting the algorithm is highly promising for application purposes. And even when birefringence is misclassified, the results are only off by a small percentage, thus suggesting that cavity tuning can still be effective and efficient. Although the algorithm was demonstrated on an underlying theoretical model, the method can be integrated directly into an experimental laser cavity design, i.e. the advocated method does not rely on an underlying model of the laser dynamics.

Philosophically, the approach taken here does not attempt to construct a better model of birefringence and its stochastic fluctuations. Rather, the cavity birefringence simply is what it is, and the goal is to use state-of-the-art data methods for classifying the birefringence into a previously learned library, i.e. the library and classifier become the expert-in-the-loop. Once classified, the optimal waveplate and polarizer settings are then already known from the toroidal search algorithm. This ultimately leads to an efficient, self-tuning laser since the waveplates and polarizers can be adjusted directly to their optimal positions using servo-control motors [22,23]. Specifics of the physical implementation are given in Fig. 8 with the training and execution algorithms given by the pseudo-codes in Table 1. The flow charts and algorithm architecture, which are color coded with each other, illustrate the practical implementation of the method.

More broadly, such difficulties in quantitative modeling are a hallmark feature of complex systems which display some underlying stochastic dynamics. For instance, the modeling of climate and weather is extremely difficult due to underlying stochastic, micro-scale physics. In such complex systems, data-driven modeling techniques are critical for improving quantitative predictions. In weather forecasting, for instance, data-assimilation methods have been critical in achieving better performance (accuracy) and longer forecast windows [25]. Thus instead of attempting to construct more refined, and typically over-parametrized models for the stochastic effects, the goal is to simply make use of recognized, coherent patterns of activity and direct measurements to inform decisions about the state of the system. The demonstration here shows that such a method is highly effective for modeling fiber lasers.

## Acknowledgments

J. N. Kutz acknowledges support from the U.S. Air Force Office of Scientific Research (AFOSR) ( FA9550-09-0174).

## References and links

**1. **D. J. Richardson, J. Nilsson, and W. A. Clarkson, “High power fiber lasers: current status and future perspectives,” J. Opt. Soc. Am. B **27**, B63–B92 (2010). [CrossRef]

**2. **C. D. Poole and R. E. Wagner, “Phenomenological approach to polarization dispersion in long single-mode fibers,” Electron. Lett. **22**, 1029–1030 (1986). [CrossRef]

**3. **C. R. Menyuk, “Pulse propagation in an elliptically birefringent Kerr media,” IEEE J. Quant. Electron. **25**, 2674–2682 (1989). [CrossRef]

**4. **C. R. Menyuk, “Nonlinear pulse propagation in birefringent optical fibers,” IEEE J. Quant. Electron. **23**, 174–176 (1987). [CrossRef]

**5. **P. K. A. Wai and C. Menyuk, “Polarization mode dispersion, decorrelation, and diffusion in optical fibers with randomly varying birefringence,” J. Light. Tech. **14**, 148–157 (1996). [CrossRef]

**6. **J. P. Gordon and H. Kogelnik, “PMD fundamentals: polarization mode dispersion in optical fibers,” PNAS **97**, 4541–4550 (2000). [CrossRef] [PubMed]

**7. **J. N. Kutz, “Mode-locked soliton lasers,” SIAM Review **48**, 629–678 (2006). [CrossRef]

**8. **H. A. Haus, “Mode-locking of lasers,” IEEE J. Sel. Top. Quant. Elec. **6**, 1173–1185 (2000). [CrossRef]

**9. **K. Tamura, E.P. Ippen, H.A. Haus, and L.E. Nelson, “77-fs Pulse generation from a stretched-pulse mode-locked all-fiber ring laser,” Opt. Lett. **18**, 1080–1082 (1993). [CrossRef] [PubMed]

**10. **K. Tamura and M. Nakazawa, “Optimizing power extraction in stretched pulse fiber ring lasers,” App. Phys. Lett. **67**, 3691–3693 (1995). [CrossRef]

**11. **G. Lenz, K. Tamura, H. A. Haus, and E. P. Ippen, “All-solid-state femtosecond source at 1.55 μm,” Opt. Lett. **20**, 1289–1291 (1995). [CrossRef] [PubMed]

**12. **A. Chong, W. H. Renninger, and F. W. Wise, “Properties of normal-dispersion femtosecond fiber lasers,” J. Opt. Soc. Am. B **25**, 140–148 (2008). [CrossRef]

**13. **A. Chong, J. Buckley, W. Renninger, and F. Wise, “All-normal-dispersion femtosecond fiber laser,” Opt. Express **14**, 10095 (2006). [CrossRef] [PubMed]

**14. **W. Renninger, A. Chong, and F. W. Wise, “Dissipative solitons in normal-dispersion fiber lasers,” Phys. Rev. A **77**, 023814 (2008). [CrossRef]

**15. **F.Ö. Ilday, J. Buckley, and F. W. Wise, “Self-similar evolution of parabolic pulses in a laser cavity,” Phys. Rev. Lett. **92**, 213902 (2004). [CrossRef]

**16. **W. H. Renninger, A. Chong, and F. W. Wise, “Self-similar pulse evolution in an all-normal-dispersion laser,” Phys. Rev. A **82**, 021805 (2010). [CrossRef]

**17. **B. Bale and S. Wabnitz, “Strong spectral filtering for a mode-locked similariton fiber laser,” Opt. Lett. **35**, 2466–2468 (2010). [CrossRef] [PubMed]

**18. **F. Li, P. K. A. Wai, and J. N. Kutz, “Geometrical description of the onset of multi-pulsing in mode-locked laser cavities,” J. Opt. Soc. Am. B **27**, 2068–2077 (2010). [CrossRef]

**19. **F. Li, E. Ding, J. N. Kutz, and P. K. A. Wai, “Dual transmission filters for enhanced energy in mode-locked fiber lasers,” Opt. Express **19**, 23408–23419 (2011). [CrossRef] [PubMed]

**20. **X. Fu and J. N. Kutz, “High-energy mode-locked fiber lasers using multiple transmission filters and a genetic algorithm,” Opt. Express **21**, 6526–6537 (2013). [CrossRef] [PubMed]

**21. **S. L. Brunton, X. Fu, and J. N. Kutz, “Extremum-seeking control of a mode-locked laser,” IEEE J. Quant. Electron. **49**, 852–861 (2013). [CrossRef]

**22. **X. Shen, W. Li, M. Yan, and H. Zeng, “Electronic control of nonlinear-polarization-rotation mode locking in Yb-doped fiber lasers,” Opt. Lett. **37**, 3426–3428 (2012). [CrossRef]

**23. **D. Radnatarov, S. Khripunov, S. Kobtsev, A. Ivanenko, and S. Kukarin, “Automatic electronic-controlled mode locking self-start in fibre lasers with non-linear polarisation evolution,” Opt. Express **21**, 20626–20631 (2013). [CrossRef] [PubMed]

**24. **S. Wiggins, *Introduction to Applied Nonlinear Dynamical Systems and Chaos*, 2, (Springer2003).

**25. **J. N. Kutz, *Data-Driven Modeling and Scientific Computation* (Oxford2013).

**26. **J. Wright, A. Yang, A. Ganesh, S. Sastry, and Y. Ma, “Robust face recognition via sparse representation,” IEEE Trans. Pattern Ana. Mach. Int. **31**, 210–227 (2009). [CrossRef]

**27. **D. Needell and J. A. Tropp, “CoSaMP: iterative signal recovery from incomplete and inaccurate samples,” Comm. of the ACM **53**, 93–100 (2010). [CrossRef]